Emu: species-level microbial community profiling of full-length 16S rRNA Oxford Nanopore sequencing data

Characterising the allergic fungal rhinosinusitis microenvironment using full-length 16S rRNA gene amplicon sequencing and fungal ITS sequencing

INTRODUCTION

Allergic fungal rhinosinusitis (AFRS) is a severe phenotype of chronic rhinosinusitis with nasal polyposis (CRSwNP), characterised by localised and exaggerated type 2 inflammation. While fungal antigenic stimulation of unregulated Th2-mediated inflammation is the core pathophysiological mechanism, the direct and synergistic role of bacteria in disease modification is a pervasive hypothesis. We set out to define the microenvironment of AFRS to elucidate virulent organisms that may be implicated in the pathophysiology of AFRS.

METHODOLOGY

We undertook a cross-sectional study of AFRS patients and non-fungal CRSwNP patients. Demographics, disease severity, culture and microbiome sequences were analysed. Multimodality microbiome sequencing included short-read next-generation sequencing (NGS) on the Illumina Miseq (16S rRNA and ITS) and full-length 16S rRNA sequencing on the Oxford Nanopore Technologies GridION (ONT).

RESULTS

Thirty-two AFRS and 29 non-fungal CRSwNP patients (NF) were included in this study. Staphylococcus aureus was the dominant organism cultured and sequenced in both AFRS and NF groups (AFRS 27.54%; NF 18.04%; p = .07). Streptococcus pneumoniae (AFRS 12.31%; NF 0.98%; p = .03) and Haemophilus influenzae (AFRS 15.03%; NF 0.24%; p = .005) were significantly more abundant in AFRS. Bacterial diversity (Shannon's index) was considerably lower in AFRS relative to NF (AFRS 0.6; NF 1.0, p = .008). Aspergillus was the most cultured fungus in AFRS (10/32, 31.3%). The AFRS sequenced mycobiome was predominantly represented by Malassezia (43.6%), Curvularia (18.5%) and Aspergillus (16.8%), while the NF mycobiome was nearly exclusively Malassezia (84.2%) with an absence of Aspergillus or dematiaceous fungi.

CONCLUSION

A low diversity, dysbiotic microenvironment dominated by Staphylococcus aureus, Streptococcus pneumoniae and Haemophilus influenzae characterised the bacterial microbiome of AFRS, with a mycobiome abundant in Malassezia, Aspergillus and Curvularia. While Staphylococcus aureus has been previously implicated in AFRS through enterotoxin superantigen potential, Streptococcus pneumoniae and Haemophilus influenzae are novel findings that may represent alternate cross-kingdom pathophysiological mechanisms.

METASEED: a novel approach to full-length 16S rRNA gene reconstruction from short read data

BACKGROUND

With the emergence of Oxford Nanopore technology, now the on-site sequencing of 16S rRNA from environments is available. Due to the error level and structure, the analysis of such data demands some database of reference sequences. However, many taxa from complex and diverse environments, have poor representation in publicly available databases. In this paper, we propose the METASEED pipeline for the reconstruction of full-length 16S sequences from such environments, in order to improve the reference for the subsequent use of on-site sequencing.

RESULTS

We show that combining high-precision short-read sequencing of both 16S and full metagenome from the same samples allow us to reconstruct high-quality 16S sequences from the more abundant taxa. A significant novelty is the carefully designed collection of metagenome reads that matches the 16S amplicons, based on a combination of uniqueness and abundance. Compared to alternative approaches this produces superior results.

CONCLUSION

Our pipeline will facilitate numerous studies associated with various unknown microorganisms, thus allowing the comprehension of the diverse environments. The pipeline is a potential tool in generating a full length 16S rRNA gene database for any environment.

Lightweight taxonomic profiling of long-read sequenced metagenomes with Lemur and Magnet

Taxonomic profiling is a ubiquitous task in the analysis of clinical and environmental microbiomes. The advent of long-read sequencing of microbiomes necessitates the development of new taxonomic profilers tailored to long-read shotgun metagenomic datasets. Here, we introduce Lemur and Magnet, a pair of tools optimized for lightweight and accurate taxonomic profiling from long-read shotgun metagenomic datasets. Lemur is a marker-gene based method that leverages an EM algorithm to reduce false positive calls while preserving true positives; Magnet makes detailed presence/absence calls for bacterial genomes based on whole-genome read mapping. The tools work in sequence: Lemur estimates abundances conservatively, and Magnet operates on the genomes of identified organisms to filter out likely false positive taxa. The result is an increase in precision of as much as 70%, which far exceeds competing methods. By operating only on marker genes, Lemur is a comparatively lightweight software. We demonstrate that it can run in minutes to hours on a laptop with 32 GB of RAM, even for large inputs - a crucial feature given the portability of long-read sequencing machines. Furthermore, the marker gene database used by Lemur is only 4 GB and contains information from over 300,000 RefSeq genomes. The reference is available at https://zenodo.org/records/10802546, and the software is open-source and available at https://github.com/treangenlab/lemur.

Comparative evaluation of analytical pipelines for illumina short- and nanopore long-read 16S rRNA gene amplicon sequencing with mock microbial communities

Utility of a recently developed long-read pipeline, Emu, was assessed using an expectation-maximization algorithm for accurate read classification. We compared it to conventional short- and long-read pipelines, using well-characterized mock bacterial samples. Our findings highlight the necessity of appropriate data-processing for taxonomic descriptions, expanding our understanding of the precise microbiome.

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.

Unveiling microbial diversity: harnessing long-read sequencing technology

Long-read sequencing has recently transformed metagenomics, enhancing strain-level pathogen characterization, enabling accurate and complete metagenome-assembled genomes, and improving microbiome taxonomic classification and profiling. These advancements are not only due to improvements in sequencing accuracy, but also happening across rapidly changing analysis methods. In this Review, we explore long-read sequencing's profound impact on metagenomics, focusing on computational pipelines for genome assembly, taxonomic characterization and variant detection, to summarize recent advancements in the field and provide an overview of available analytical methods to fully leverage long reads. We provide insights into the advantages and disadvantages of long reads over short reads and their evolution from the early days of long-read sequencing to their recent impact on metagenomics and clinical diagnostics. We further point out remaining challenges for the field such as the integration of methylation signals in sub-strain analysis and the lack of benchmarks.

Impact of nocturnal hypoxia on glycaemic control, appetite, gut microbiota and inflammation in adults with type 2 diabetes mellitus: A single-blind cross-over trial

High altitude residents have a lower incidence of type 2 diabetes mellitus (T2DM). Therefore, we examined the effect of repeated overnight normobaric hypoxic exposure on glycaemic control, appetite, gut microbiota and inflammation in adults with T2DM. Thirteen adults with T2DM [glycated haemoglobin (HbA1c): 61.1 ± 14.1 mmol mol-1; aged 64.2 ± 9.4 years; four female] completed a single-blind, randomised, sham-controlled, cross-over study for 10 nights, sleeping when exposed to hypoxia (fractional inspired O2 [F I O 2 ${{F}_{{\mathrm{I}}{{{\mathrm{O}}}_{\mathrm{2}}}}}$ ] = 0.155; ∼2500 m simulated altitude) or normoxic conditions (F I O 2 ${{F}_{{\mathrm{I}}{{{\mathrm{O}}}_{\mathrm{2}}}}}$  = 0.209) in a randomised order. Outcome measures included: fasted plasma [glucose]; [hypoxia inducible factor-1α]; [interleukin-6]; [tumour necrosis factor-α]; [interleukin-10]; [heat shock protein 70]; [butyric acid]; peak plasma [glucose] and insulin sensitivity following a 2 h oral glucose tolerance test; body composition; appetite indices ([leptin], [acyl ghrelin], [peptide YY], [glucagon-like peptide-1]); and gut microbiota diversity and abundance [16S rRNA amplicon sequencing]. During intervention periods, accelerometers measured physical activity, sleep duration and efficiency, whereas continuous glucose monitors were used to assess estimated HbA1c and glucose management indicator and time in target range. Overnight hypoxia was not associated with changes in any outcome measure (P > 0.05 with small effect sizes) except fasting insulin sensitivity and gut microbiota alpha diversity, which exhibited trends (P = 0.10; P = 0.08 respectively) for a medium beneficial effect (d = 0.49; d = 0.59 respectively). Ten nights of overnight moderate hypoxic exposure did not significantly affect glycaemic control, gut microbiome, appetite, or inflammation in adults with T2DM. However, the intervention was well tolerated and a medium effect-size for improved insulin sensitivity and reduced alpha diversity warrants further investigation. KEY POINTS: Living at altitude lowers the incidence of type 2 diabetes mellitus (T2DM). Animal studies suggest that exposure to hypoxia may lead to weight loss and suppressed appetite. In a single-blind, randomised sham-controlled, cross-over trial, we assessed the effects of 10 nights of hypoxia (fractional inspired O2 ∼0.155) on glucose homeostasis, appetite, gut microbiota, inflammatory stress ([interleukin-6]; [tumour necrosis factor-α]; [interleukin-10]) and hypoxic stress ([hypoxia inducible factor 1α]; heat shock protein 70]) in 13 adults with T2DM. Appetite and inflammatory markers were unchanged following hypoxic exposure, but an increased insulin sensitivity and reduced gut microbiota alpha diversity were associated with a medium effect-size and statistical trends, which warrant further investigation using a definitive large randomised controlled trial. Hypoxic exposure may represent a viable therapeutic intervention in people with T2DM and particularly those unable or unwilling to exercise because barriers to uptake and adherence may be lower than for other lifestyle interventions (e.g. diet and exercise).

Sterile sentinels and MinION sequencing capture active soil microbial communities that differentiate crop rotations

BACKGROUND

Soil microbial communities are difficult to measure and critical to soil processes. The bulk soil microbiome is highly diverse and spatially heterogeneous, which can make it difficult to detect and monitor the responses of microbial communities to differences or changes in management, such as different crop rotations in agricultural research. Sampling a subset of actively growing microbes should promote monitoring how soil microbial communities respond to management by reducing the variation contributed by high microbial spatial and temporal heterogeneity and less active microbes. We tested an in-growth bag method using sterilized soil in root-excluding mesh, "sterile sentinels," for the capacity to differentiate between crop rotations. We assessed the utility of different incubation times and compared colonized sentinels to concurrently sampled bulk soils for the statistical power to differentiate microbial community composition in low and high diversity crop rotations. We paired this method with Oxford Nanopore MinION sequencing to assess sterile sentinels as a standardized, fast turn-around monitoring method.

RESULTS

Compared to bulk soil, sentinels provided greater statistical power to distinguish between crop rotations for bacterial communities and equivalent power for fungal communities. The incubation time did not affect the statistical power to detect treatment differences in community composition, although longer incubation time increased total biomass. Bulk and sentinel soil samples contained shared and unique microbial taxa that were differentially abundant between crop rotations.

CONCLUSIONS

Overall, compared to bulk soils, the sentinels captured taxa with copiotrophic or ruderal traits, and plant-associated taxa. The sentinels show promise as a sensitive, scalable method to monitor soil microbial communities and provide information complementary to traditional soil sampling.

A comparison between full-length 16S rRNA Oxford nanopore sequencing and Illumina V3-V4 16S rRNA sequencing in head and neck cancer tissues

Describing the microbial community within the tumour has been a key aspect in understanding the pathophysiology of the tumour microenvironment. In head and neck cancer (HNC), most studies on tissue samples have only performed 16S rRNA short-read sequencing (SRS) on V3-V5 region. SRS is mostly limited to genus level identification. In this study, we compared full-length 16S rRNA long-read sequencing (FL-ONT) from Oxford Nanopore Technology (ONT) to V3-V4 Illumina SRS (V3V4-Illumina) in 26 HNC tumour tissues. Further validation was also performed using culture-based methods in 16 bacterial isolates obtained from 4 patients using MALDI-TOF MS. We observed similar alpha diversity indexes between FL-ONT and V3V4-Illumina. However, beta-diversity was significantly different between techniques (PERMANOVA - R2 = 0.131, p < 0.0001). At higher taxonomic levels (Phylum to Family), all metrics were more similar among sequencing techniques, while lower taxonomy displayed more discrepancies. At higher taxonomic levels, correlation in relative abundance from FL-ONT and V3V4-Illumina were higher, while this correlation decreased at lower levels. Finally, FL-ONT was able to identify more isolates at the species level that were identified using MALDI-TOF MS (75% vs. 18.8%). FL-ONT was able to identify lower taxonomic levels at a better resolution as compared to V3V4-Illumina 16S rRNA sequencing.

Sensitivity and consistency of long- and short-read metagenomics and epicPCR for the detection of antibiotic resistance genes and their bacterial hosts in wastewater

Wastewater surveillance is a powerful tool to assess the risks associated with antibiotic resistance in communities. One challenge is selecting which analytical tool to deploy to measure risk indicators, such as antibiotic resistance genes (ARGs) and their respective bacterial hosts. Although metagenomics is frequently used for analyzing ARGs, few studies have compared the performance of long-read and short-read metagenomics in identifying which bacteria harbor ARGs in wastewater. Furthermore, for ARG host detection, untargeted metagenomics has not been compared to targeted methods such as epicPCR. Here, we 1) evaluated long-read and short-read metagenomics as well as epicPCR for detecting ARG hosts in wastewater, and 2) investigated the host range of ARGs across the wastewater treatment plant (WWTP) to evaluate host proliferation. Results highlighted long-read revealed a wider range of ARG hosts compared to short-read metagenomics. Nonetheless, the ARG host range detected by long-read metagenomics only represented a subset of the hosts detected by epicPCR. The ARG-host linkages across the influent and effluent of the WWTP were characterized. Results showed the ARG-host phylum linkages were relatively consistent across the WWTP, whereas new ARG-host species linkages appeared in the WWTP effluent. The ARG-host linkages of several clinically relevant species found in the effluent were identified.

Tomato domestication rather than subsequent breeding events reduces microbial associations related to phosphorus recovery

Legacy phosphorus (P) is a reservoir of sparingly available P, and its recovery could enhance sustainable use of nonrenewable mineral fertilizers. Domestication has affected P acquisition, but it is unknown if subsequent breeding efforts, like the Green Revolution (GR), had a similar effect. We examined how domestication and breeding events altered P acquisition by growing wild, traditional (pre-GR), and modern (post-GR) tomato in soil with legacy P but low bioavailable P. Wild tomatoes, particularly accession LA0716 (Solanum pennellii), heavily cultured rhizosphere P solubilizers, suggesting reliance on microbial associations to acquire P. Wild tomato also had a greater abundance of other putatively beneficial bacteria, including those that produce chelating agents and antibiotic compounds. Although wild tomatoes had a high abundance of these P solubilizers, they had lower relative biomass and greater P stress factor than traditional or modern tomato. Compared to wild tomato, domesticated tomato was more tolerant to P deficiency, and both cultivated groups had a similar rhizosphere bacterial community composition. Ultimately, this study suggests that while domestication changed tomato P recovery by reducing microbial associations, subsequent breeding processes have not further impacted microbial P acquisition mechanisms. Selecting microbial P-related traits that diminished with domestication may therefore increase legacy P solubilization.

Human nasal microbiota shifts in healthy and chronic respiratory disease conditions

BACKGROUND

An increasing number of studies investigate various human microbiotas and their roles in the development of diseases, maintenance of health states, and balanced signaling towards the brain. Current data demonstrate that the nasal microbiota contains a unique and highly variable array of commensal bacteria and opportunistic pathogens. However, we need to understand how to harness current knowledge, enrich nasal microbiota with beneficial microorganisms, and prevent pathogenic developments.

RESULTS

In this study, we have obtained nasal, nasopharyngeal, and bronchoalveolar lavage fluid samples from healthy volunteers and patients suffering from chronic respiratory tract diseases for full-length 16 S rRNA sequencing analysis using Oxford Nanopore Technologies. Demographic and clinical data were collected simultaneously. The microbiome analysis of 97 people from Lithuania suffering from chronic inflammatory respiratory tract disease and healthy volunteers revealed that the human nasal microbiome represents the microbiome of the upper airways well.

CONCLUSIONS

The nasal microbiota of patients was enriched with opportunistic pathogens, which could be used as indicators of respiratory tract conditions. In addition, we observed that a healthy human nasal microbiome contained several plant- and bee-associated species, suggesting the possibility of enriching human nasal microbiota via such exposures when needed. These candidate probiotics should be investigated for their modulating effects on airway and lung epithelia, immunogenic properties, neurotransmitter content, and roles in maintaining respiratory health and nose-brain interrelationships.

Effect of the consumption of brazzein and monellin, two recombinant sweet-tasting proteins, on rat gut microbiota

Sweet-tasting proteins (SPs) are proteins of plant origin initially isolated from tropical fruits. They are thousands of times sweeter than sucrose and most artificial sweeteners. SPs are a class of proteins capable of causing a sweet taste sensation in humans when interacting with the T1R2/T1R3 receptor. SP thaumatin has already been introduced in the food industry in some countries. Other SPs, such as monellin and brazzein, are promising products. An important stage in researching SPs, in addition to confirming the absence of toxicity, mutagenicity, oncogenicity, and allergenic effects, is studying their influence on gut microbiota. In this paper we describe changes in the composition of rat gut microbiota after six months of consuming one of two recombinant SPs-brazzein or monellin. A full length 16S gene sequencing method was used for DNA library barcoding. The MaAsLin2 analysis results showed noticeable fluctuations in the relative abundances of Anaerocella delicata in brazzein-fed rat microbiota, and of Anaerutruncus rubiinfantis in monellin-fed rat microbiota, which, however, did not exceed the standard deviation. The sucrose-fed group was associated with an increase in the relative abundance of Faecalibaculum rodentium, which may contribute to obesity. Overall, prolonged consumption of the sweet proteins brazzein and monellin did not significantly change rat microbiota and did not result in the appearance of opportunistic microbiota. This provides additional evidence for the safety of these potential sweeteners.

Bioaccessibility of lead and cadmium in soils around typical lead-acid power plants and their effect on gut microorganisms

Potentially toxic elements (Pb and Cd) contamination of soil can adversely affect human health. Moreover, these metal ions interact with the gut microbiota after entering the human digestive system. Based on the physiologically based extraction test and the simulator of human intestinal microbial ecosystem, the bioaccessibility of Pb and Cd in soils contaminated with lead-acid power plants was assessed. The gastric stage exhibited the greatest average bioaccessibility of lead and cadmium (63.39% and 57.22%), followed by the small intestinal stage (6.86% and 36.29%); due to gut microorganisms, the bioaccessibility of lead and cadmium was further reduced in the colon stage (1.86% and 4.22%). Furthermore, to investigate soil contamination's effects on gut microbes, 16S rRNA high-throughput sequencing was used to identify the gut microbial species after the colon period. Due to Pb and Cd exposure, the relative abundance of Firmicutes and unidentified_Bacteria decreased, while the relative abundance of Proteobacteria, Synergistota, and Bacteroidota increased. The relationship between environmental factors and the number of microbial species in the gut was also examined using Spearman correlation analysis. Pb and Cd exposure has been found to affect the composition and structure of the gut microbiota.

Grazing-driven shifts in soil bacterial community structure and function in a typical steppe are mediated by additional N inputs

Herbivore grazing and nitrogen (N) fertilization affect soil microbial diversity and community composition both in direct and indirect pathways (e.g., via alterations in soil microenvironment and plant communities); however, their combination effects are still largely unexplored. We carried out a field study to investigate how soil abiotic properties, plant community composition and functional traits altered soil bacterial community structure and function in response to a long-term herbivore grazing (17-year sheep grazing with four stocking rates) and anthropogenic N inputs (6-year N addition with four levels) experiment. We show that a high stocking rate of 8.7 sheep ha-1 (SR8.7) decreased soil bacterial α- and β-diversity, while α- and β-diversity showed hump-shaped and saddle-shaped responses, respectively, with increasing N addition rate, reaching tipping points at the N application rate of 10 g N m-2 year-1 (N10). The synergistic effects of grazing and N addition induced the highest soil bacterial α-diversity at SR2.7 with N10. The contrasting effects of grazing and N addition induced higher soil bacterial β-diversity at SR8.7 with N20. Plant factors (e.g., aboveground biomass of Stipa bungeana and community-weighted mean carbon [CWM_C]), edaphic factors (e.g., soil moisture, pH, NO3--N, and C:nutrients ratios) and their interactions were the most significant factors affecting the diversity and community composition of bacteria. Our structure equation model (SEM) shows that grazing-induced negative effects on soil pH and plant community composition indirectly increased the β-diversity of soil bacteria, while grazing-induced decreased CWM_C had positive effects on bacterial α-diversity and community structure. However, N addition indirectly increased β-diversity of soil bacteria via changes in soil NO3--N and plant community composition, while N addition had negative impacts on bacterial α-diversity and community structure via variations in CWM_C. The interaction of grazing and N addition increased the complexity and stability of the bacterial network. Based on the KEGG database, grazing and N addition could accelerate the soil functional potential of C and N cycling. Our findings suggest that N application at a rate of <10 g N m-2 year-1 with a stocking rate of <5.3 sheep ha-1 could maintain the development of soil bacteria in supporting the most important ecosystem functions and services. Complex responses of soil microbes to grazing and N addition indicate the need for deeper investigations of the impacts of global change on microbial involvement in biogeochemical cycles.

An integrated method for targeted Oxford Nanopore sequencing and automated bioinformatics for the simultaneous detection of bacteria, fungi, and ARG

AIMS

The use of metagenomics for pathogen identification in clinical practice has been limited. Here we describe a workflow to encourage the clinical utility and potential of NGS for the screening of bacteria, fungi, and antimicrobial resistance genes (ARGs).

METHODS AND RESULTS

The method includes target enrichment, long-read sequencing, and automated bioinformatics. Evaluation of several tools and databases was undertaken across standard organisms (n = 12), clinical isolates (n = 114), and blood samples from patients with suspected bloodstream infections (n = 33). The strategy used could offset the presence of host background DNA, error rates of long-read sequencing, and provide accurate and reproducible detection of pathogens. Eleven targets could be successfully tested in a single assay. Organisms could be confidently identified considering ≥60% of best hits of a BLAST-based threshold of e-value 0.001 and a percent identity of >80%. For ARGs, reads with percent identity of >90% and >60% overlap of the complete gene could be confidently annotated. A kappa of 0.83 was observed compared to standard diagnostic methods. Thus, a workflow for the direct-from-sample, on-site sequencing combined with automated genomics was demonstrated to be reproducible.

CONCLUSION

NGS-based technologies overcome several limitations of current day diagnostics. Highly sensitive and comprehensive methods of pathogen screening are the need of the hour. We developed a framework for reliable, on-site, screening of pathogens.

Identification of consensus head and neck cancer-associated microbiota signatures: a systematic review and meta-analysis of 16S rRNA and The Cancer Microbiome Atlas datasets

Introduction. Multiple reports have attempted to describe the tumour microbiota in head and neck cancer (HNSC).Gap statement. However, these have failed to produce a consistent microbiota signature, which may undermine understanding the importance of bacterial-mediated effects in HNSC.Aim. The aim of this study is to consolidate these datasets and identify a consensus microbiota signature in HNSC.Methodology. We analysed 12 published HNSC 16S rRNA microbial datasets collected from cancer, cancer-adjacent and non-cancer tissues to generate a consensus microbiota signature. These signatures were then validated using The Cancer Microbiome Atlas (TCMA) database and correlated with the tumour microenvironment phenotypes and patient's clinical outcome.Results. We identified a consensus microbial signature at the genus level to differentiate between HNSC sample types, with cancer and cancer-adjacent tissues sharing more similarity than non-cancer tissues. Univariate analysis on 16S rRNA datasets identified significant differences in the abundance of 34 bacterial genera among the tissue types. Paired cancer and cancer-adjacent tissue analyses in 16S rRNA and TCMA datasets identified increased abundance in Fusobacterium in cancer tissues and decreased abundance of Atopobium, Rothia and Actinomyces in cancer-adjacent tissues. Furthermore, these bacteria were associated with different tumour microenvironment phenotypes. Notably, high Fusobacterium signature was associated with high neutrophil (r=0.37, P<0.0001), angiogenesis (r=0.38, P<0.0001) and granulocyte signatures (r=0.38, P<0.0001) and better overall patient survival [continuous: HR 0.8482, 95 % confidence interval (CI) 0.7758-0.9273, P=0.0003].Conclusion. Our meta-analysis demonstrates a consensus microbiota signature for HNSC, highlighting its potential importance in this disease.

Changes of putative pathogenic species within the water bacterial community in large-scale drinking water treatment and distribution systems

Although the management of microbes in drinking water is of paramount importance for public health, there remain challenges in comprehensively examining pathogenic bacteria in the water supply system at the species level. In this study, high-throughput sequencing of nearly full-length 16S rRNA genes was performed to investigate the changes of the water bacterial community in three large-scale drinking water treatment plants (DWTPs) and their corresponding distribution systems during winter and summer. Our findings revealed significant differences in the bacterial community structure between winter and summer water samples for each DWTP and its distribution management area (DMA). In the groundwater-fed DWTP, selective enrichment of mycobacterial species was observed in both seasons, and the subsequent DMA also exhibited strong selection for specific mycobacterial species. In one of the surface water-fed DWTPs, certain Legionella species present in the source water in winter were selectively enriched in the bacterial community after pre-oxidation, although they were susceptible to the subsequent purification steps. A variety of putative pathogenic species (n = 83) were identified based on our pathogen identification pipeline, with the dominant species representing opportunistic pathogens commonly found in water supply systems. While pathogen removal primarily occurred during the purification processes of DWTPs, especially for surface water-fed plants, the relative abundance of pathogenic bacteria in the DMA water flora was lower than that in the DWTP effluent flora, indicating a diminished competitiveness of pathogens within the DMA ecosystem.

The pan-microbiome profiling system Taxa4Meta identifies clinical dysbiotic features and classifies diarrheal disease

Targeted metagenomic sequencing is an emerging strategy to survey disease-specific microbiome biomarkers for clinical diagnosis and prognosis. However, this approach often yields inconsistent or conflicting results owing to inadequate study power and sequencing bias. We introduce Taxa4Meta, a bioinformatics pipeline explicitly designed to compensate for technical and demographic bias. We designed and validated Taxa4Meta for accurate taxonomic profiling of 16S rRNA amplicon data acquired from different sequencing strategies. Taxa4Meta offers significant potential in identifying clinical dysbiotic features that can reliably predict human disease, validated comprehensively via reanalysis of individual patient 16S data sets. We leveraged the power of Taxa4Meta's pan-microbiome profiling to generate 16S-based classifiers that exhibited excellent utility for stratification of diarrheal patients with Clostridioides difficile infection, irritable bowel syndrome, or inflammatory bowel diseases, which represent common misdiagnoses and pose significant challenges for clinical management. We believe that Taxa4Meta represents a new "best practices" approach to individual microbiome surveys that can be used to define gut dysbiosis at a population-scale level.

The clinical utility of Nanopore 16S rRNA gene sequencing for direct bacterial identification in normally sterile body fluids

The prolonged incubation period of traditional culture methods leads to a delay in diagnosing invasive infections. Nanopore 16S rRNA gene sequencing (Nanopore 16S) offers a potential rapid diagnostic approach for directly identifying bacteria in infected body fluids. To evaluate the clinical utility of Nanopore 16S, we conducted a study involving the collection and sequencing of 128 monomicrobial samples, 65 polymicrobial samples, and 20 culture-negative body fluids. To minimize classification bias, taxonomic classification was performed using 3 analysis pipelines: Epi2me, Emu, and NanoCLUST. The result was compared to the culture references. The limit of detection of Nanopore 16S was also determined using simulated bacteremic blood samples. Among the three classifiers, Emu demonstrated the highest concordance with the culture results. It correctly identified the taxon of 125 (97.7%) of the 128 monomicrobial samples, compared to 109 (85.2%) for Epi2me and 102 (79.7%) for NanoCLUST. For the 230 cultured species in the 65 polymicrobial samples, Emu correctly identified 188 (81.7%) cultured species, compared to 174 (75.7%) for Epi2me and 125 (54.3%) for NanoCLUST. Through ROC analysis on the monomicrobial samples, we determined a threshold of relative abundance at 0.058 for distinguishing potential pathogens from background in Nanopore 16S. Applying this threshold resulted in the identification of 107 (83.6%), 117 (91.4%), and 114 (91.2%) correctly detected samples for Epi2me, Emu, and NanoCLUST, respectively, in the monomicrobial samples. Nanopore 16S coupled with Epi2me could provide preliminary results within 6 h. However, the ROC analysis of polymicrobial samples exhibited a random-like performance, making it difficult to establish a threshold. The overall limit of detection for Nanopore 16S was found to be about 90 CFU/ml.

KombOver: Efficient k-core and K-truss based characterization of perturbations within the human gut microbiome

The microbes present in the human gastrointestinal tract are regularly linked to human health and disease outcomes. Thanks to technological and methodological advances in recent years, metagenomic sequencing data, and computational methods designed to analyze metagenomic data, have contributed to improved understanding of the link between the human gut microbiome and disease. However, while numerous methods have been recently developed to extract quantitative and qualitative results from host-associated microbiome data, improved computational tools are still needed to track microbiome dynamics with short-read sequencing data. Previously we have proposed KOMB as a de novo tool for identifying copy number variations in metagenomes for characterizing microbial genome dynamics in response to perturbations. In this work, we present KombOver (KO), which includes four key contributions with respect to our previous work: (i) it scales to large microbiome study cohorts, (ii) it includes both k-core and K-truss based analysis, (iii) we provide the foundation of a theoretical understanding of the relation between various graph-based metagenome representations, and (iv) we provide an improved user experience with easier-to-run code and more descriptive outputs/results. To highlight the aforementioned benefits, we applied KO to nearly 1000 human microbiome samples, requiring less than 10 minutes and 10 GB RAM per sample to process these data. Furthermore, we highlight how graph-based approaches such as k-core and K-truss can be informative for pinpointing microbial community dynamics within a myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) cohort. KO is open source and available for download/use at: https://github.com/treangenlab/komb.

Direct Determination of the Structure of Single Biopolymer Molecules Using Nanopore Sequencing

This review highlights operational principles, features, and modern aspects of the development of third-generation sequencing technology of biopolymers focusing on the nucleic acids analysis, namely the nanopore sequencing system. Basics of the method and technical solutions used for its realization are considered, from the first works showing the possibility of creation of these systems to the easy-to-handle procedure developed by Oxford Nanopore Technologies company. Moreover, this review focuses on applications, which were developed and realized using equipment developed by the Oxford Nanopore Technologies, including assembly of whole genomes, methagenomics, direct analysis of the presence of modified bases.

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.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

Application of Nanopore Sequencing in the Diagnosis and Treatment of Pulmonary Infections

This review provides an in-depth discussion of the development, principles and utility of nanopore sequencing technology and its diverse applications in the identification of various pulmonary pathogens. We examined the emergence and advancements of nanopore sequencing as a significant player in this field. We illustrate the challenges faced in diagnosing mixed infections and further scrutinize the use of nanopore sequencing in the identification of single pathogens, including viruses (with a focus on its use in epidemiology, outbreak investigation, and viral resistance), bacteria (emphasizing 16S targeted sequencing, rare bacterial lung infections, and antimicrobial resistance studies), fungi (employing internal transcribed spacer sequencing), tuberculosis, and atypical pathogens. Furthermore, we discuss the role of nanopore sequencing in metagenomics and its potential for unbiased detection of all pathogens in a clinical setting, emphasizing its advantages in sequencing genome repeat areas and structural variant regions. We discuss the limitations in dealing with host DNA removal, the inherent high error rate of nanopore sequencing technology, along with the complexity of operation and processing, while acknowledging the possibilities provided by recent technological improvements. We compared nanopore sequencing with the BioFire system, a rapid molecular diagnostic system based on polymerase chain reaction. Although the BioFire system serves well for the rapid screening of known and common pathogens, it falls short in the identification of unknown or rare pathogens and in providing comprehensive genome analysis. As technological advancements continue, it is anticipated that the role of nanopore sequencing technology in diagnosing and treating lung infections will become increasingly significant.

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.

Exploring the Feasibility of 16S rRNA Short Amplicon Sequencing-Based Microbiota Analysis for Microbiological Safety Assessment of Raw Oyster

16S rRNA short amplicon sequencing-based microbiota profiling has been thought of and suggested as a feasible method to assess food safety. However, even if a comprehensive microbial information can be obtained by microbiota profiling, it would not be necessarily sufficient for all circumstances. To prove this, the feasibility of the most widely used V3-V4 amplicon sequencing method for food safety assessment was examined here. We designed a pathogen (Vibrio parahaemolyticus) contamination and/or V. parahaemolyticus-specific phage treatment model of raw oysters under improper storage temperature and monitored their microbial structure changes. The samples stored at refrigerator temperature (negative control, NC) and those that were stored at room temperature without any treatment (no treatment, NT) were included as control groups. The profiling results revealed that no statistical difference exists between the NT group and the pathogen spiked- and/or phage treated-groups even when the bacterial composition was compared at the possible lowest-rank taxa, family/genus level. In the beta-diversity analysis, all the samples except the NC group formed one distinct cluster. Notably, the samples with pathogen and/or phage addition did not form each cluster even though the enumerated number of V. parahaemolyticus in those samples were extremely different. These discrepant results indicate that the feasibility of 16S rRNA short amplicon sequencing should not be overgeneralized in microbiological safety assessment of food samples, such as raw oyster.

CuReSim-LoRM: A Tool to Simulate Metabarcoding Long Reads

Metabarcoding DNA sequencing has revolutionized the study of microbial communities. Third-generation sequencing producing long reads had opened up new perspectives. Obtaining the full-length ribosomal RNA gene would permit one to reach a better taxonomic resolution at the species or the strain level. However, Oxford Nanopore Technologies (ONT) sequencing produces reads with high error rates, which introduces biases in analysis. Understanding the biases introduced during the analysis allows one to better interpret the biological results and take care of conclusions drawn from metabarcoding experiments. To benchmark an analysis process, the ground truth, i.e., the real composition of the microbial community, has to be known. In addition to artificial mock communities, simulated data are often used to evaluate the biases and performances of the bioinformatics analysis step. Currently, no specific tool has been developed to simulate metabarcoding long reads, mimic the error rate and the length distribution, and allow one to benchmark the analysis process. Here, we introduce CuReSim-LoRM, for the customized read simulator to generate long reads for metabarcoding. We showed that CuReSim-LoRM is able to produce reads with varying error rates and length distributions by mimicking the real data very well.

Altered Respiratory Microbiomes, Plasma Metabolites, and Immune Responses in Influenza A Virus and Methicillin-Resistant Staphylococcus aureus Coinfection

Influenza A virus (IAV)-methicillin-resistant Staphylococcus aureus (MRSA) coinfection causes severe respiratory infections. The host microbiome plays an important role in respiratory tract infections. However, the relationships among the immune responses, metabolic characteristics, and respiratory microbial characteristics of IAV-MRSA coinfection have not been fully studied. We used specific-pathogen-free (SPF) C57BL/6N mice infected with IAV and MRSA to build a nonlethal model of IAV-MRSA coinfection and characterized the upper respiratory tract (URT) and lower respiratory tract (LRT) microbiomes at 4 and 13 days postinfection by full-length 16S rRNA gene sequencing. Immune response and plasma metabolism profile analyses were performed at 4 days postinfection by flow cytometry and liquid chromatography-tandem mass spectrometry (LC-MS/MS). The relationships among the LRT microbiota, the immune response, and the plasma metabolism profile were analyzed by Spearman's correlation analysis. IAV-MRSA coinfection showed significant weight loss and lung injury and significantly increased loads of IAV and MRSA in bronchoalveolar lavage fluid (BALF). Microbiome data showed that coinfection significantly increased the relative abundances of Enterococcus faecalis, Enterobacter hormaechei, Citrobacter freundii, and Klebsiella pneumoniae and decreased the relative abundances of Lactobacillus reuteri and Lactobacillus murinus. The percentages of CD4+/CD8+ T cells and B cells in the spleen; the levels of interleukin-9 (IL-9), interferon gamma (IFN-γ), tumor necrosis factor alpha (TNF-α), IL-6, and IL-8 in the lung; and the level of mevalonolactone in plasma were increased in IAV-MRSA-coinfected mice. L. murinus was positively correlated with lung macrophages and natural killer (NK) cells, negatively correlated with spleen B cells and CD4+/CD8+ T cells, and correlated with multiple plasma metabolites. Future research is needed to clarify whether L. murinus mediates or alters the severity of IAV-MRSA coinfection. IMPORTANCE The respiratory microbiome plays an important role in respiratory tract infections. In this study, we characterized the URT and LRT microbiota, the host immune response, and plasma metabolic profiles during IAV-MRSA coinfection and evaluated their correlations. We observed that IAV-MRSA coinfection induced severe lung injury and dysregulated host immunity and plasma metabolic profiles, as evidenced by the aggravation of lung pathological damage, the reduction of innate immune cells, the strong adaptation of the immune response, and the upregulation of mevalonolactone in plasma. L. murinus was strongly correlated with immune cells and plasma metabolites. Our findings contribute to a better understanding of the role of the host microbiome in respiratory tract infections and identified a key bacterial species, L. murinus, that may provide important references for the development of probiotic therapies.

Dysbiotic Co-Factors in Cervical Cancer. How the Microbiome Influences the Development of Cervical Intraepithelial Neoplasia (CIN)

Human papillomavirus (HPV) infection is a necessary but not sufficient condition for the development of cervical cancer. The dysbiotic shift in the cervicovaginal microbiome appears to be a major co-factor in carcinogenesis. New analytical methods, such as next-generation sequencing (NGS), can be used to detect all of the vaginal microorganisms present and therefore identify individual therapeutic options. The relationship of bacterial vaginosis and carcinogenesis, as well as possible indications for the use of microbiome analysis, will be discussed.

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.

Commonalities between the Atacama Desert and Antarctica rhizosphere microbial communities

Plant-microbiota interactions have significant effects on plant growth, health, and productivity. Rhizosphere microorganisms are involved in processes that promote physiological responses to biotic and abiotic stresses in plants. In recent years, the interest in microorganisms to improve plant productivity has increased, mainly aiming to find promising strains to overcome the impact of climate change on crops. In this work, we hypothesize that given the desertic environment of the Antarctic and the Atacama Desert, different plant species inhabiting these areas might share microbial taxa with functions associated with desiccation and drought stress tolerance. Therefore, in this study, we described and compared the composition of the rhizobacterial community associated with Deschampsia antarctica (Da), Colobanthus quitensis (Cq) from Antarctic territories, and Croton chilensis (Cc), Eulychnia iquiquensis (Ei) and Nicotiana solanifolia (Ns) from coastal Atacama Desert environments by using 16S rRNA amplicon sequencing. In addition, we evaluated the putative functions of that rhizobacterial community that are likely involved in nutrient acquisition and stress tolerance of these plants. Even though each plant microbial rhizosphere presents a unique taxonomic pattern of 3,019 different sequences, the distribution at the genus level showed a core microbiome with a higher abundance of Haliangium, Bryobacter, Bacillus, MND1 from the Nitrosomonadaceae family, and unclassified taxa from Gemmatiamonadaceae and Chitinophagaceae families in the rhizosphere of all samples analyzed (781 unique sequences). In addition, species Gemmatirosa kalamazoonesis and Solibacter usitatus were shared by the core microbiome of both Antarctic and Desert plants. All the taxa mentioned above had been previously associated with beneficial effects in plants. Also, this microbial core composition converged with the functional prediction related to survival under harsh conditions, including chemoheterotrophy, ureolysis, phototrophy, nitrogen fixation, and chitinolysis. Therefore, this study provides relevant information for the exploration of rhizospheric microorganisms from plants in extreme conditions of the Atacama Desert and Antarctic as promising plant growth-promoting rhizobacteria.

Comprehensive Assessment of 16S rRNA Gene Amplicon Sequencing for Microbiome Profiling across Multiple Habitats

The 16S rRNA gene works as a rapid and effective marker for the identification of microorganisms in complex communities; hence, a huge number of microbiomes have been surveyed by 16S amplicon-based sequencing. The resolution of the 16S rRNA gene is always considered only at the genus level; however, it has not been verified on a wide range of microbes yet. To fully explore the ability and potential of the 16S rRNA gene in microbial profiling, here, we propose Qscore, a comprehensive method to evaluate the performance of amplicons by integrating the amplification rate, multitier taxonomic annotation, sequence type, and length. Our in silico assessment by a "global view" of 35,889 microbe species across multiple reference databases summarizes the optimal sequencing strategy for 16S short reads. On the other hand, since microbes are unevenly distributed according to their habitats, we also provide the recommended configuration for 16 typical ecosystems based on the Qscores of 157,390 microbiomes in the Microbiome Search Engine (MSE). Detailed data simulation further proves that the 16S amplicons produced with Qscore-suggested parameters exhibit high precision in microbiome profiling, which is close to that of shotgun metagenomes under CAMI metrics. Therefore, by reconsidering the precision of 16S-based microbiome profiling, our work not only enables the high-quality reusability of massive sequence legacy that has already been produced but is also significant for guiding microbiome studies in the future. We have implemented the Qscore as an online service at http://qscore.single-cell.cn to parse the recommended sequencing strategy for specific habitats or expected microbial structures. IMPORTANCE 16S rRNA has long been used as a biomarker to identify distinct microbes from complex communities. However, due to the influence of the amplification region, sequencing type, sequence processing, and reference database, the accuracy of 16S rRNA has not been fully verified on a global range. More importantly, the microbial composition of different habitats varies greatly, and it is necessary to adopt different strategies according to the corresponding target microbes to achieve optimal analytical performance. Here, we developed Qscore, which evaluates the comprehensive performance of 16S amplicons from multiple perspectives, thus providing the best sequencing strategies for common ecological environments by using big data.

Comparison of Oxford Nanopore Technologies and Illumina MiSeq sequencing with mock communities and agricultural soil

Illumina MiSeq is the current standard for characterizing microbial communities in soil. The newer alternative, Oxford Nanopore Technologies MinION sequencer, is quickly gaining popularity because of the low initial cost and longer sequence reads. However, the accuracy of MinION, per base, is much lower than MiSeq (95% versus 99.9%). The effects of this difference in base-calling accuracy on taxonomic and diversity estimates remains unclear. We compared the effects of platform, primers, and bioinformatics on mock community and agricultural soil samples using short MiSeq, and short and full-length MinION 16S rRNA amplicon sequencing. For all three methods, we found that taxonomic assignments of the mock community at both the genus and species level matched expectations with minimal deviation (genus: 80.9-90.5%; species: 70.9-85.2% Bray-Curtis similarity); however, the short MiSeq with error correction (DADA2) resulted in the correct estimate of mock community species richness and much lower alpha diversity for soils. Several filtering strategies were tested to improve these estimates with varying results. The sequencing platform also had a significant influence on the relative abundances of taxa with MiSeq resulting in significantly higher abundances Actinobacteria, Chloroflexi, and Gemmatimonadetes and lower abundances of Acidobacteria, Bacteroides, Firmicutes, Proteobacteria, and Verrucomicrobia compared to the MinION platform. When comparing agricultural soils from two different sites (Fort Collins, CO and Pendleton, OR), methods varied in the taxa identified as significantly different between sites. At all taxonomic levels, the full-length MinION method had the highest similarity to the short MiSeq method with DADA2 correction with 73.2%, 69.3%, 74.1%, 79.3%, 79.4%, and 82.28% of the taxa at the phyla, class, order, family, genus, and species levels, respectively, showing similar patterns in differences between the sites. In summary, although both platforms appear suitable for 16S rRNA microbial community composition, biases for different taxa may make the comparison between studies problematic; and even with a single study (i.e., comparing sites or treatments), the sequencing platform can influence the differentially abundant taxa identified.

Cefquinome shows a higher impact on the pig gut microbiome and resistome compared to ceftiofur

Cephalosporins are licensed for treatment of severe bacterial infections in different species. However, the effect of these antimicrobials on the fecal microbiome and potential spread of resistance-associated genes causes great concern. This highlights the need to understand the impact of cephalosporins on the porcine fecal microbiome and resistome. A combination of long-read 16S rRNA gene and shotgun metagenomic sequencing was applied to investigate the effect of conventional treatment with either ceftiofur (3 mg.kg^-1 intramuscular, 3 consecutive days) or cefquinome (2 mg.kg^-1 intramuscular, 5 consecutive days) on the porcine microbiome and resistome. Fecal samples were collected from 17 pigs (6 ceftiofur treated, 6 cefquinome treated, 5 control pigs) at four different timepoints. Treatment with ceftiofur resulted in an increase in Proteobacteria members on microbiome level, while on resistome level selection in TetQ containing Bacteroides, CfxA6 containing Prevotella and bla_TEM-1 containing Escherichia coli was observed. Cefquinome treatment resulted in a decline in overall species richness (α-diversity) and increase in Proteobacteria members. On genus level, administration of cefquinome significantly affected more genera than ceftiofur (18 vs 8). On resistome level, cefquinome resulted in a significant increase of six antimicrobial resistance genes, with no clear correlation with certain genera. For both antimicrobials, the resistome levels returned back to the control levels 21 days post-treatment. Overall, our study provides novel insights on the effect of specific cephalosporins on the porcine gut microbiome and resistome after conventional intramuscular treatment. These results might contribute to better tailoring of the most ideal treatment strategy for some bacterial infections.

The role of microbiome-host interactions in the development of Alzheimer´s disease

Alzheimer`s disease (AD) is the most prevalent cause of dementia. It is often assumed that AD is caused by an aggregation of extracellular beta-amyloid and intracellular tau-protein, supported by a recent study showing reduced brain amyloid levels and reduced cognitive decline under treatment with a beta-amyloid-binding antibody. Confirmation of the importance of amyloid as a therapeutic target notwithstanding, the underlying causes of beta-amyloid aggregation in the human brain, however, remain to be elucidated. Multiple lines of evidence point towards an important role of infectious agents and/or inflammatory conditions in the etiology of AD. Various microorganisms have been detected in the cerebrospinal fluid and brains of AD-patients and have thus been hypothesized to be linked to the development of AD, including Porphyromonas gingivalis (PG) and Spirochaetes. Intriguingly, these microorganisms are also found in the oral cavity under normal physiological conditions, which is often affected by multiple pathologies like caries or tooth loss in AD patients. Oral cavity pathologies are mostly accompanied by a compositional shift in the community of oral microbiota, mainly affecting commensal microorganisms and referred to as 'dysbiosis'. Oral dysbiosis seems to be at least partly mediated by key pathogens such as PG, and it is associated with a pro-inflammatory state that promotes the destruction of connective tissue in the mouth, possibly enabling the translocation of pathogenic microbiota from the oral cavity to the nervous system. It has therefore been hypothesized that dysbiosis of the oral microbiome may contribute to the development of AD. In this review, we discuss the infectious hypothesis of AD in the light of the oral microbiome and microbiome-host interactions, which may contribute to or even cause the development of AD. We discuss technical challenges relating to the detection of microorganisms in relevant body fluids and approaches for avoiding false-positives, and introduce the antibacterial protein lactoferrin as a potential link between the dysbiotic microbiome and the host inflammatory reaction.

Competition between skin antimicrobial peptides and commensal bacteria in type 2 inflammation enables survival of S. aureus

During inflammation, the skin deploys antimicrobial peptides (AMPs) yet during allergic inflammation it becomes more susceptible to Staphylococcus aureus. To understand this contradiction, single-cell sequencing of Il4ra^-/- mice combined with skin microbiome analysis reveals that lower production of AMPs from interleukin-4 receptor α (IL-4Rα) activation selectively inhibits survival of antibiotic-producing strains of coagulase-negative Staphylococcus (CoNS). Diminished AMPs under conditions of T helper type 2 (Th2) inflammation enable expansion of CoNS strains without antibiotic activity and increase Staphylococcus aureus (S. aureus), recapitulating the microbiome on humans with atopic dermatitis. This response is rescued in Camp^-/- mice or after topical steroids, since further inhibition of AMPs enables survival of antibiotic-producing CoNS strains. In conditions of Th17 inflammation, a higher expression of host AMPs is sufficient to directly inhibit S. aureus survival. These results show that antimicrobials produced by the host and commensal bacteria each act to control S. aureus on the skin.

Real-Time Nanopore Q20+ Sequencing Enables Extremely Fast and Accurate Core Genome MLST Typing and Democratizes Access to High-Resolution Bacterial Pathogen Surveillance

Next-generation whole-genome sequencing is essential for high-resolution surveillance of bacterial pathogens, for example, during outbreak investigations or for source tracking and escape variant analysis. However, current global sequencing and bioinformatic bottlenecks and a long time to result with standard technologies demand new approaches. In this study, we investigated whether novel nanopore Q20+ long-read chemistry enables standardized and easily accessible high-resolution typing combined with core genome multilocus sequence typing (cgMLST). We set high requirements for discriminatory power by using the slowly evolving bacterium Bordetella pertussis as a model pathogen. Our results show that the increased raw read accuracy enables the description of epidemiological scenarios and phylogenetic linkages at the level of gold-standard short reads. The same was true for our variant analysis of vaccine antigens, resistance genes, and virulence factors, demonstrating that nanopore sequencing is a legitimate competitor in the area of next-generation sequencing (NGS)-based high-resolution bacterial typing. Furthermore, we evaluated the parameters for the fastest possible analysis of the data. By combining the optimized processing pipeline with real-time basecalling, we established a workflow that allows for highly accurate and extremely fast high-resolution typing of bacterial pathogens while sequencing is still in progress. Along with advantages such as low costs and portability, the approach suggested here might democratize modern bacterial typing, enabling more efficient infection control globally.

SituSeq: an offline protocol for rapid and remote Nanopore 16S rRNA amplicon sequence analysis

Microbiome analysis through 16S rRNA gene sequencing is a crucial tool for understanding the microbial ecology of any habitat or ecosystem. However, workflows require large equipment, stable internet, and extensive computing power such that most of the work is performed far away from sample collection in both space and time. Performing amplicon sequencing and analysis at sample collection would have positive implications in many instances including remote fieldwork and point-of-care medical diagnoses. Here we present SituSeq, an offline and portable workflow for the sequencing and analysis of 16S rRNA gene amplicons using Nanopore sequencing and a standard laptop computer. SituSeq was validated by comparing Nanopore 16S rRNA gene amplicons, Illumina 16S rRNA gene amplicons, and Illumina metagenomes, sequenced using the same environmental DNA. Comparisons revealed consistent community composition, ecological trends, and sequence identity across platforms. Correlation between the abundance of taxa in each taxonomic level in Illumina and Nanopore data sets was high (Pearson's r > 0.9), and over 70% of Illumina 16S rRNA gene sequences matched a Nanopore sequence with greater than 97% sequence identity. On board a research vessel on the open ocean, SituSeq was used to analyze amplicon sequences from deep sea sediments less than 2 h after sequencing, and 8 h after sample collection. The rapidly available results informed decisions about subsequent sampling in near real-time while the offshore expedition was still underway. SituSeq is a portable and user-friendly workflow that helps to bring the power of microbial genomics and diagnostics to many more researchers and situations.

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.

Inoculation with black soldier fly larvae alters the microbiome and volatile organic compound profile of decomposing food waste

The black soldier fly (BSF; Hermetia illucens) is used in sustainable processing of many types of organic waste. However, organic waste being decomposed by BSF produces strong odors, hindering more widespread application. The odor components and how they are produced have yet to be characterized. We found that digestion of food waste by BSF significantly alters the microbial flora, based on metagenomic analyses, and the odor components generated, as shown by thermal desorption gas chromatography mass spectrometry analysis. Inoculation with BSF significantly decreased production of volatile organic sulfur compounds (dimethyl disulfide and dimethyl trisulfide), which are known to be released during methionine and cysteine metabolism by Lactobacillus and Enterococcus bacteria. BSF inoculation significantly changed the abundance of Lactobacillus and Enterococcus and decreased microbial diversity overall. These findings may help in optimizing use of BSF for deodorization of composting food waste.

Characterisation of the Gut Bacteria of Cultured and Wild Spiny Lobster Panulirus ornatus

The commercial onshore aquaculture of the spiny lobster Panulirus ornatus, while in its infancy, has progressed rapidly from the enabling research that continues at the University of Tasmania. The development of lobster feeds, both fresh and manufactured, has been critical to the success of this emerging aquaculture sector. Fresh feeds derived from mussel represent the gold standard in terms of the growth performance of juvenile lobsters. Nonetheless, concerns regarding availability, sustainability, and potential biosecurity issues of fresh feeds highlight the importance of developing manufactured feeds for lobster aquaculture. Wild lobsters are assumed to have a balanced natural diet that allows for standard growth and development, and as such natural diets are often used as a reference for feed development. Similarly, the gut microbiota associated with a natural diet is assumed to reflect a healthy microbial assemblage. The aim of this study was to compare the microbiota of the hindgut and hepatopancreas of cultured P. ornatus fed with a commercial prawn pellet or mussel to that of wild spiny lobster juveniles. Gut samples were analysed using Oxford Nanopore 16S rRNA gene sequencing. Based on principal coordinate analysis, the gut bacteria of cultured lobsters were different from the wild juveniles. The core microbiota of the hindgut and hepatopancreas libraries were phyla Proteobacteria (Gamma, Alpha) and Bacteroidetes. Vibrio was the most dominant genus in both organs. The differences in bacterial relative abundance were mainly between cultured (pellet-, mussel-fed) and wild lobsters. In conclusion, bacteria in the cultured lobsters had significantly different profiles to that of the wild juveniles, indicating that current onshore aquaculture practices alter the gut microbiota. A number of different feeding and culture practices may be required if the aim of closed culture practices is to attain a gut microbiota in cultured animals that is representative of that found in wild spiny lobsters.

Unexpected finding of Fusobacterium varium as the dominant Fusobacterium species in cattle rumen: potential implications for liver abscess etiology and interventions

Fusobacterium varium has been generally overlooked in cattle rumen microbiome studies relative to the presumably more abundant liver abscess-causing Fusobacterium necrophorum. However, F. varium was found to be more abundant in the rumen fluid of cattle and under culture conditions tailored to enrich F. necrophorum. Using near-full length 16S ribosomal ribonucleic acid sequencing, we demonstrate that F. varium grows under restrictive conditions commonly used to enumerate F. necrophorum, suggesting that previous F. necrophorum abundance assessment may have been inaccurate and that F. varium may be an underestimated member of the ruminal bacterial community. Fusobacterium varium were not as susceptible as F. necrophorum to in-feed antibiotics conventionally used in feedlots. Exposure to tylosin, the current gold standard for liver abscess reduction strategies in cattle, consistently hindered growth of the F. necrophorum strains tested by over 67% (P < 0.05) relative to the unexposed control. In contrast, F. varium strains were totally or highly resistant (0%-13% reduction in maximum yield, P < 0.05). Monensin, an ionophore antibiotic, had greater inhibitory activity against F. necrophorum than F. varium. Finally, preliminary genomic analysis of two F. varium isolates from the rumen revealed the presence of virulence genes related to those of pathogenic human F. varium isolates associated with active invasion of mammalian cells. The data presented here encourage further investigation into the ecological role of F. varium within the bovine rumen and potential role in liver abscess development, and proactive interventions.

The Neolithic site "La Marmotta": DNA metabarcoding to identify the microbial deterioration of waterlogged archeological wood

Introduction

The evaluation of biological degradation of waterlogged archeological wood is crucial to choose the conservative and protective treatments to be applied to the wooden material. The waterlogged environmental conditions are characterized by oxygen scarcity, only allowing the growth of adapted microbes capable to degrade the organic wooden material, mainly erosion bacteria and soft-rot fungi. In this work, we characterized and evaluated the biodegradation state and the microbial communities of wooden fragments preserved in storage tanks. These were preserved by waterlogging within the Neolithic village "La Marmotta," currently found under the Bracciano Lake (Lazio, Italy).

Methods

The waterlogged wood samples were first identified taxonomically with an optical microscope, also allowing an evaluation of their preservation state. The microbial community was then evaluated through the sequencing of Internal Transcribed Spacer sequences for fungi and 16S for bacteria with the Oxford Nanopore Technologies (ONT) MinION platform.

Results

The identified microbial community appears to be consistent with the waterlogged samples, as many bacteria attributable to the erosion of wood and ligninolytic fungi have been sequenced.

Discussion

The reported results highlight the first use of targeted metabarcoding by ONT applied to study the biodeterioration of waterlogged archeological wood.

Identification of microbial community in the urban environment: The concordance between conventional culture and nanopore 16S rRNA sequencing

Introduction

Microbes in the built environment have been implicated as a source of infectious diseases. Bacterial culture is the standard method for assessing the risk of exposure to pathogens in urban environments, but this method only accounts for <1% of the diversity of bacteria. Recently, full-length 16S rRNA gene analysis using nanopore sequencing has been applied for microbial evaluations, resulting in a rise in the development of long-read taxonomic tools for species-level classification. Regarding their comparative performance, there is, however, a lack of information.

Methods

Here, we aim to analyze the concordance of the microbial community in the urban environment inferred by multiple taxonomic classifiers, including ARGpore2, Emu, Kraken2/Bracken and NanoCLUST, using our 16S-nanopore dataset generated by MegaBLAST, as well as assess their abilities to identify culturable species based on the conventional culture results.

Results

According to our results, NanoCLUST was preferred for 16S microbial profiling because it had a high concordance of dominant species and a similar microbial profile to MegaBLAST, whereas Kraken2/Bracken, which had similar clustering results as NanoCLUST, was also desirable. Second, for culturable species identification, Emu with the highest accuracy (81.2%) and F1 score (29%) for the detection of culturable species was suggested.

Discussion

In addition to generating datasets in complex communities for future benchmarking studies, our comprehensive evaluation of the taxonomic classifiers offers recommendations for ongoing microbial community research, particularly for complex communities using nanopore 16S rRNA sequencing.

Topography of respiratory tract and gut microbiota in mice with influenza A virus infection

Introduction

Influenza A virus (IAV)-induced dysbiosis may predispose to severe bacterial superinfections. Most studies have focused on the microbiota of single mucosal surfaces; consequently, the relationships between microbiota at different anatomic sites in IAV-infected mice have not been fully studied.

Methods

We characterized respiratory and gut microbiota using full-length 16S rRNA gene sequencing by Nanopore sequencers and compared the nasopharyngeal, oropharyngeal, lung and gut microbiomes in healthy and IAV-infected mice.

Results

The oropharyngeal, lung and gut microbiota of healthy mice were dominated by Lactobacillus spp., while nasopharyngeal microbiota were comprised primarily of Streptococcus spp. However, the oropharyngeal, nasopharyngeal, lung, and gut microbiota of IAV-infected mice were dominated by Pseudomonas, Escherichia, Streptococcus, and Muribaculum spp., respectively. Lactobacillus murinus was identified as a biomarker and was reduced at all sites in IAV-infected mice. The microbiota composition of lung was more similar to that of the nasopharynx than the oropharynx in healthy mice.

Discussion

These findings suggest that the main source of lung microbiota in mice differs from that of adults. Moreover, the similarity between the nasopharyngeal and lung microbiota was increased in IAV-infected mice. We found that IAV infection reduced the similarity between the gut and oropharyngeal microbiota. L. murinus was identified as a biomarker of IAV infection and may be an important target for intervention in post-influenza bacterial superinfections.

Improving analysis of the vaginal microbiota of women undergoing assisted reproduction using nanopore sequencing

PURPOSE

Subclinical alterations of the vaginal microbiome have been described to be associated with female infertility and may serve as predictors for failure of in vitro fertilization treatment. While large prospective studies to delineate the role of microbial composition are warranted, integrating microbiome information into clinical management depends on economical and practical feasibility, specifically on a short duration from sampling to final results. The currently most used method for microbiota analysis is either metagenomics sequencing or amplicon-based microbiota analysis using second-generation methods such as sequencing-by-synthesis approaches (Illumina), which is both expensive and time-consuming. Thus, additional approaches are warranted to accelerate the usability of the microbiome as a marker in clinical praxis.

METHODS

Herein, we used a set of ten selected vaginal swabs from women undergoing assisted reproduction, comparing and performing critical optimization of nanopore-based microbiota analysis with the results from MiSeq-based data as a quality reference.

RESULTS

The analyzed samples carried varying community compositions, as shown by amplicon-based analysis of the V3V4 region of the bacterial 16S rRNA gene by MiSeq sequencing. Using a stepwise procedure to optimize adaptation, we show that a close approximation of the microbial composition can be achieved within a reduced time frame and at a minimum of costs using nanopore sequencing.

CONCLUSIONS

Our work highlights the potential of a nanopore-based methodical setup to support the feasibility of interventional studies and contribute to the development of microbiome-based clinical decision-making in assisted reproduction.