Search

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

Metabarcoding is an important tool for understanding fungal communities. The internal transcribed spacer (ITS) rDNA is the accepted fungal barcode but has known problems. The large subunit (LSU) rDNA has also been used to investigate fungal communities but available LSU metabarcoding primers were mostly designed to target Dikarya (Ascomycota + Basidiomycota) with little attention to early diverging fungi (EDF). However, evidence from multiple studies suggests that EDF comprise a large portion of unknown diversity in community sampling. Here, we investigate how DNA marker choice and methodological biases impact recovery of EDF from environmental samples. We focused on one EDF lineage, Zoopagomycota, as an example. We evaluated three primer sets (ITS1F/ITS2, LROR/LR3, and LR3 paired with new primer LR22F) to amplify and sequence a Zoopagomycota mock community and a set of 146 environmental samples with Illumina MiSeq. We compared two taxonomy assignment methods and created an LSU reference database compatible with AMPtk software. The two taxonomy assignment methods recovered strikingly different communities of fungi and EDF. Target fragment length variation exacerbated PCR amplification biases and influenced downstream taxonomic assignments, but this effect was greater for EDF than Dikarya. To improve identification of LSU amplicons we performed phylogenetic reconstruction and illustrate the advantages of this critical tool for investigating identified and unidentified sequences. Our results suggest much of the EDF community may be missed or misidentified with "standard" metabarcoding approaches and modified techniques are needed to understand the role of these taxa in a broader ecological context.

The gut mycobiome plays an important role in the health and disease of the human gut, but its exact function is still under investigation. While there is a wealth of information available on the bacterial community of the human gut microbiome, research on the fungal community is still relatively limited. In particular, technical methodologies for mycobiome analysis, especially the DNA extraction method for human faecal samples, varied in different studies. In the current study, two commercial kits commonly used in DNA extraction, the QIAamp® Fast DNA Stool Mini Kit and DNeasy PowerSoil Pro Kit, and one manual method, the International Human Microbiome Standards Protocol Q, were compared. Furthermore, the effectiveness of two different bead-beating machines, the Mini-Beadbeater-16 and FastPrep-24TM 5G, was compared in parallel. A mock fungal community with a known composition of fungal strains was also generated and included to compare different DNA extraction methods. Our results suggested that the method using the DNeasy PowerSoil Pro Kit and Mini-Beadbeater-16 provides the best results to extract DNA from human faecal samples. Based on our data, we propose a standard operating procedure for DNA extraction from human faecal samples for mycobiome analysis.

BackgroundLong-read shotgun metagenomic sequencing is gaining in popularity and offers many advantages over short-read sequencing. The higher information content in long reads is useful for a variety of metagenomics analyses, including taxonomic classification and profiling. The development of long-read specific tools for taxonomic classification is accelerating, yet there is a lack of information regarding their relative performance. Here, we perform a critical benchmarking study using 11 methods, including five methods designed specifically for long reads. We applied these tools to several mock community datasets generated using Pacific Biosciences (PacBio) HiFi or Oxford Nanopore Technology sequencing, and evaluated their performance based on read utilization, detection metrics, and relative abundance estimates.ResultsOur results show that long-read classifiers generally performed best. Several short-read classification and profiling methods produced many false positives (particularly at lower abundances), required heavy filtering to achieve acceptable precision (at the cost of reduced recall), and produced inaccurate abundance estimates. By contrast, two long-read methods (BugSeq, MEGAN-LR & DIAMOND) and one generalized method (sourmash) displayed high precision and recall without any filtering required. Furthermore, in the PacBio HiFi datasets these methods detected all species down to the 0.1% abundance level with high precision. Some long-read methods, such as MetaMaps and MMseqs2, required moderate filtering to reduce false positives to resemble the precision and recall of the top-performing methods. We found read quality affected performance for methods relying on protein prediction or exact k-mer matching, and these methods performed better with PacBio HiFi datasets. We also found that long-read datasets with a large proportion of shorter reads (

Abstract The small subunit (SSU) ribosomal RNA (rRNA) is the most commonly used marker for the identification of microbial taxa, but its full‐length reconstruction from high‐throughput sequencing (HTS) data remains challenging. Metagenomics and total RNA sequencing (total RNA‐Seq) are target‐PCR‐free HTS methods that are used to characterize microbial communities and simultaneously reconstruct SSU rRNA sequences. However, more testing is required to determine and improve their effectiveness. We processed metagenomics and total RNA‐Seq data retrieved from a commercially available mock microbial community and an aquarium sample using 112 combinations of data processing tools. We determined the SSU rRNA reconstruction completeness of both sequencing methods for both samples and analysed the impact of data processing tools on SSU rRNA completeness. In contrast to metagenomics, total RNA‐Seq allowed for the complete or near‐complete reconstruction of all mock community SSU rRNA sequences and generated up to 438 SSU rRNA sequences with ≥80% completeness from the aquarium sample using only 1/5 of an Illumina MiSeq run. SSU rRNA completeness of metagenomics significantly correlated with the genome size of mock community species. Data processing tools impacted SSU rRNA completeness, in particular the utilized assemblers. These results are promising for the high‐throughput reconstruction of novel full‐length SSU rRNA sequences and could advance the simultaneous application of multiple ‐omics approaches in routine environmental assessments to allow for more holistic assessments of ecosystems.

The small subunit (SSU) ribosomal RNA (rRNA) is the most commonly used marker for the identification of microbial taxa, but its full‐length reconstruction from high‐throughput sequencing (HTS) data remains challenging. Metagenomics and total RNA sequencing (total RNA‐Seq) are target‐PCR‐free HTS methods that are used to characterize microbial communities and simultaneously reconstruct SSU rRNA sequences. However, more testing is required to determine and improve their effectiveness.We processed metagenomics and total RNA‐Seq data retrieved from a commercially available mock microbial community and an aquarium sample using 112 combinations of data processing tools. We determined the SSU rRNA reconstruction completeness of both sequencing methods for both samples and analysed the impact of data processing tools on SSU rRNA completeness.In contrast to metagenomics, total RNA‐Seq allowed for the complete or near‐complete reconstruction of all mock community SSU rRNA sequences and generated up to 438 SSU rRNA sequences with ≥80% completeness from the aquarium sample using only 1/5 of an Illumina MiSeq run. SSU rRNA completeness of metagenomics significantly correlated with the genome size of mock community species. Data processing tools impacted SSU rRNA completeness, in particular the utilized assemblers.These results are promising for the high‐throughput reconstruction of novel full‐length SSU rRNA sequences and could advance the simultaneous application of multiple ‐omics approaches in routine environmental assessments to allow for more holistic assessments of ecosystems. [ABSTRACT FROM AUTHOR]

Abstract Receiving the same results from repeated analysis of the same sample is a basic principle in science. The inability to reproduce previously published results has led to discussions of a reproducibility crisis within science. For studies of microbial communities, the problem of reproducibility is more pronounced and has, in some fields, led to a discussion on the very existence of a constantly present microbiota. In this study, DNA from 44 bovine milk samples were extracted twice and the V3–V4 region of the 16S rRNA gene was sequenced in two separate runs. The FASTQ files from the two data sets were run through the same bioinformatics pipeline using the same settings and results from the two data sets were compared. Milk samples collected maximally 2 h apart were used as replicates and permitted comparisons to be made within the same run. Results show a significant difference in species richness between the two sequencing runs although Shannon and Simpson's diversity was the same. Multivariate analyses of all samples demonstrate that the sequencing run was a driver for variation. Direct comparison of similarity between samples and sequencing run showed an average similarity of 42%–45% depending on whether binary or abundance‐based similarity indices were used. Within‐run comparisons of milk samples collected maximally 2 h apart showed an average similarity of 39%–47% depending on the similarity index used and that similarity differed significantly between runs. We conclude that repeated DNA extraction and sequencing significantly can affect the results of a low microbial biomass microbiota study.

Recent advances in new molecular biology methods and next-generation sequencing (NGS) technologies have revolutionized metabarcoding studies investigating complex microbial communities from various environments. The inevitable first step in sample preparation is DNA extraction which introduces its own set of biases and considerations. In this study, we assessed the influence of five DNA extraction methods [B1: phenol/chloroform/isoamyl extraction, B2 and B3: isopropanol and ethanol precipitations, respectively--both modifications of B1, K1: DNeasy PowerWater Kit (QIAGEN), K2: modified DNeasy PowerWater Kit (QIAGEN) and direct PCR approach (P) that completely circumvents this step on community composition and DNA yield of mock and marine sample communities from the Adriatic Sea]. B1--B3 methods generally produced higher DNA yields and more similar microbial communities, but with higher interindividual variability. Each method demonstrated significant differences in a specific community structure, where rare taxa seem to play a crucial role. There was not one superior method closest to the theoretically expected mock community composition, they all demonstrated skewed ratios, but in a similar way which might be attributed to other factors, such as primer bias or 16S rRNA gene count for specific taxa. Direct PCR represents an interesting approach when high throughput in sample processing is required. We emphasize the importance of making a cautious decision about the choice of the extraction method or direct PCR approach, but even more importantly its consistent application throughout the study. [ABSTRACT FROM AUTHOR]

Introduction: The gut microbiome is an integral partner in host health and plays a role in immune development, altered nutrition, and pathogen prevention. The mycobiome (fungal microbiome) is considered part of the rare biosphere but is still a critical component in health. Next generation sequencing has improved our understanding of fungi in the gut, but methodological challenges remain. Biases are introduced during DNA isolation, primer design and choice, polymerase selection, sequencing platform selection, and data analyses, as fungal reference databases are often incomplete or contain erroneous sequences. Methods: Here, we compared the accuracy of taxonomic identifications and abundances from mycobiome analyses which vary among three commonly selected target gene regions (18S, ITS1, or ITS2) and the reference database (UNITE - ITS1, ITS2 and SILVA - 18S). We analyze multiple communities including individual fungal isolates, a mixed mock community created from five common fungal isolates found in weanling piglet feces, a purchased commercial fungal mock community, and piglet fecal samples. In addition, we calculated gene copy numbers for the 18S, ITS1, and ITS2 regions of each of the five isolates from the piglet fecal mock community to determine whether copy number affects abundance estimates. Finally, we determined the abundance of taxa from several iterations of our in-house fecal community to assess the effects of community composition on taxon abundance. Results: Overall, no marker-database combination consistently outperformed the others. Internal transcribed space markers were slightly superior to 18S in the identification of species in tested communities, but Lichtheimia corymbifera, a common member of piglet gut communities, was not amplified by ITS1 and ITS2 primers. Thus, ITS based abundance estimates of taxa in piglet mock communities were skewed while 18S marker profiles were more accurate. Kazachstania slooffiae displayed the most stable copy numbers (83-85) while L. corymbifera displayed significant variability (90-144) across gene regions. Discussion: This study underscores the importance of preliminary studies to assess primer combinations and database choice for the mycobiome sample of interest and raises questions regarding the validity of fungal abundance estimates. [ABSTRACT FROM AUTHOR]

Recent advances in new molecular biology methods and next-generation sequencing (NGS) technologies have revolutionized metabarcoding studies investigating complex microbial communities from various environments. The inevitable first step in sample preparation is DNA extraction which introduces its own set of biases and considerations. In this study, we assessed the influence of five DNA extraction methods [B1: phenol/chloroform/isoamyl extraction, B2 and B3: isopropanol and ethanol precipitations, respectively—both modifications of B1, K1: DNeasy PowerWater Kit (QIAGEN), K2: modified DNeasy PowerWater Kit (QIAGEN) and direct PCR approach (P) that completely circumvents this step on community composition and DNA yield of mock and marine sample communities from the Adriatic Sea]. B1–B3 methods generally produced higher DNA yields and more similar microbial communities, but with higher interindividual variability. Each method demonstrated significant differences in a specific community structure, where rare taxa seem to play a crucial role. There was not one superior method closest to the theoretically expected mock community composition, they all demonstrated skewed ratios, but in a similar way which might be attributed to other factors, such as primer bias or 16S rRNA gene count for specific taxa. Direct PCR represents an interesting approach when high throughput in sample processing is required. We emphasize the importance of making a cautious decision about the choice of the extraction method or direct PCR approach, but even more importantly its consistent application throughout the study.

Microbial communities contain a broad phylogenetic diversity of organisms; however, the majority of methods center on describing bacteria and archaea. Fungi are important symbionts in many ecosystems and are potentially important members of the human microbiome, beyond those that can cause disease. To expand our analysis of microbial communities to include data from the fungal internal transcribed spacer (ITS) region, five candidate DNA extraction kits were compared against our standardized protocol for describing bacteria and archaea using 16S rRNA gene amplicon- and shotgun metagenomics sequencing. The results are presented considering a diverse panel of host-associated and environmental sample types and comparing the cost, processing time, well-to-well contamination, DNA yield, limit of detection and microbial community composition among protocols. Across all criteria, the MagMAX Microbiome kit was found to perform best. The PowerSoil Pro kit performed comparably but with increased cost per sample and overall processing time. The Zymo MagBead, NucleoMag Food and Norgen Stool kits were included.

IntroductionThe gut microbiome is an integral partner in host health and plays a role in immune development, altered nutrition, and pathogen prevention. The mycobiome (fungal microbiome) is considered part of the rare biosphere but is still a critical component in health. Next generation sequencing has improved our understanding of fungi in the gut, but methodological challenges remain. Biases are introduced during DNA isolation, primer design and choice, polymerase selection, sequencing platform selection, and data analyses, as fungal reference databases are often incomplete or contain erroneous sequences.MethodsHere, we compared the accuracy of taxonomic identifications and abundances from mycobiome analyses which vary among three commonly selected target gene regions (18S, ITS1, or ITS2) and the reference database (UNITE - ITS1, ITS2 and SILVA - 18S). We analyze multiple communities including individual fungal isolates, a mixed mock community created from five common fungal isolates found in weanling piglet feces, a purchased commercial fungal mock community, and piglet fecal samples. In addition, we calculated gene copy numbers for the 18S, ITS1, and ITS2 regions of each of the five isolates from the piglet fecal mock community to determine whether copy number affects abundance estimates. Finally, we determined the abundance of taxa from several iterations of our in-house fecal community to assess the effects of community composition on taxon abundance.ResultsOverall, no marker-database combination consistently outperformed the others. Internal transcribed space markers were slightly superior to 18S in the identification of species in tested communities, but Lichtheimia corymbifera, a common member of piglet gut communities, was not amplified by ITS1 and ITS2 primers. Thus, ITS based abundance estimates of taxa in piglet mock communities were skewed while 18S marker profiles were more accurate. Kazachstania slooffiae displayed the most stable copy numbers (83-85) while L. corymbifera displayed significant variability (90-144) across gene regions.DiscussionThis study underscores the importance of preliminary studies to assess primer combinations and database choice for the mycobiome sample of interest and raises questions regarding the validity of fungal abundance estimates.

It is known that data from both 16S and shotgun metagenomics studies are subject to biases that cause the observed relative abundances of taxa to differ from their true values. Model community analyses, in which the relative abundances of all taxa in the sample are known by construction, seem to offer the hope that these biases can be measured. However, it is unclear whether the bias we measure in a mock community analysis is the same as we measure in a sample in which taxa are spiked in at known relative abundance, or if the biases we measure in spike-in samples is the same as the bias we would measure in a real (e.g., biological) sample. Here, we consider these questions in the context of 16S rRNA measurements on three sets of samples: the commercially available Zymo cells model community; the Zymo model community mixed with Swedish Snus, a smokeless tobacco product that is virtually bacteria-free; and a set of commercially available smokeless tobacco products. Each set of samples was subject to four different extraction protocols. The goal of our analysis is to determine whether the patterns of bias observed in each set of samples are the same, i.e., can we learn about the bias in the commercially available smokeless tobacco products by studying the Zymo cells model community? [ABSTRACT FROM AUTHOR]

Microbial community analysis based on the 16S rRNA-gene is used to investigate both beneficial and harmful microorganisms in various fields and environments. Recently, the next-generation sequencing (NGS) technology has enabled rapid and accurate microbial community analysis. Despite these advantages of NGS based metagenomics study, sample transport, storage conditions, amplification, library preparation kits, sequencing, and bioinformatics procedures can bias microbial community analysis results. In this study, eight mock communities were pooled from genomic DNA of Lactobacillus acidophilus KCTC 3164T, Limosilactobacillus fermentum KCTC 3112T, Lactobacillus gasseri KCTC 3163T, Lacticaseibacillus paracasei subsp. paracasei KCTC 3510T, Limosilactobacillus reuteri KCTC 3594T, Lactococcus lactis subsp. lactis KCTC 3769T, Bifidobacterium animalis subsp. lactis KCTC 5854T, and Bifidobacterium breve KCTC 3220T. The genomic DNAs were quantified by droplet digital PCR (ddPCR) and were mixed as mock communities. The mock communities were amplified with various 16S rRNA gene universal primer pairs and sequenced by MiSeq, IonTorrent, MGIseq-2000, Sequel II, and MinION NGS platforms. In a comparison of primer-dependent bias, the microbial profiles of V1-V2 and V3 regions were similar to the original ratio of the mock communities, while the microbial profiles of the V1-V3 region were relatively biased. In a comparison of platform-dependent bias, the sequence read from short-read platforms (MiSeq, IonTorrent, and MGIseq-2000) showed lower bias than that of long-read platforms (Sequel II and MinION). Meanwhile, the sequences read from Sequel II and MinION platforms were relatively biased in some mock communities. In the data of all NGS platforms and regions, L. acidophilus was greatly underrepresented while Lactococcus lactis subsp. lactis was generally overrepresented. In all samples of this study, the bias index (BI) was calculated and PCA was performed for comparison. The samples with biased relative abundance showed high BI values and were separated in the PCA results. In particular, analysis of regions rich in AT and GC poses problems for genome assembly, which can lead to sequencing bias. According to this comparative analysis, the development of reference material (RM) material has been proposed to calibrate the bias in microbiome analysis. [ABSTRACT FROM AUTHOR]

Microbial community analysis based on the 16S rRNA-gene is used to investigate both beneficial and harmful microorganisms in various fields and environments. Recently, the next-generation sequencing (NGS) technology has enabled rapid and accurate microbial community analysis. Despite these advantages of NGS based metagenomics study, sample transport, storage conditions, amplification, library preparation kits, sequencing, and bioinformatics procedures can bias microbial community analysis results. In this study, eight mock communities were pooled from genomic DNA of Lactobacillus acidophilus KCTC 3164T, Limosilactobacillus fermentum KCTC 3112T, Lactobacillus gasseri KCTC 3163T, Lacticaseibacillus paracasei subsp. paracasei KCTC 3510T, Limosilactobacillus reuteri KCTC 3594T, Lactococcus lactis subsp. lactis KCTC 3769T, Bifidobacterium animalis subsp. lactis KCTC 5854T, and Bifidobacterium breve KCTC 3220T. The genomic DNAs were quantified by droplet digital PCR (ddPCR) and were mixed as mock communities. The mock communities were amplified with various 16S rRNA gene universal primer pairs and sequenced by MiSeq, IonTorrent, MGIseq-2000, Sequel II, and MinION NGS platforms. In a comparison of primer-dependent bias, the microbial profiles of V1-V2 and V3 regions were similar to the original ratio of the mock communities, while the microbial profiles of the V1-V3 region were relatively biased. In a comparison of platform-dependent bias, the sequence read from short-read platforms (MiSeq, IonTorrent, and MGIseq-2000) showed lower bias than that of long-read platforms (Sequel II and MinION). Meanwhile, the sequences read from Sequel II and MinION platforms were relatively biased in some mock communities. In the data of all NGS platforms and regions, L. acidophilus was greatly underrepresented while Lactococcus lactis subsp. lactis was generally overrepresented. In all samples of this study, the bias index (BI) was calculated and PCA was performed for comparison. The samples with biased relative abundance showed high BI values and were separated in the PCA results. In particular, analysis of regions rich in AT and GC poses problems for genome assembly, which can lead to sequencing bias. According to this comparative analysis, the development of reference material (RM) material has been proposed to calibrate the bias in microbiome analysis.

ABSTRACT The bovine udder is colonized by a huge quantity of microorganisms that constitute the intramammary ecosystem, with a specific role in modulating not only udder homeostasis and mastitis susceptibility, but also the quality of the dairy products. However, generating high-quality bacterial DNA can be critical, especially starting from a complex biological matrix like milk, characterized by high fat, protein, and calcium contents. Here, bacterial DNA was recovered from a commercial ultra-high-temperature (UHT) milk sample artificially spiked with a predetermined mock community composition and from three bulk tank milk (raw milk) samples. The DNA was isolated using three different protocols to evaluate the effect of the extraction procedures on the milk microbiota composition. In the mock community experiment, the bacterial profiles generated by the three DNA extraction protocols were profoundly different, with the genera Staphylococcus, Lactobacillus, Listeria, and Salmonella underestimated by all the protocols. Only one protocol revealed values close to the expected abundances for Escherichia/Shigella spp., Bacillus spp., Enterococcus spp., and Pseudomonas spp. On the other hand, the nonspiked UHT milk sample exhibited a similar microbiota composition, revealing the prevalence of Acinetobacter spp., for all the DNA extraction protocols. For the raw milk samples, the three DNA extraction kits performed differently, revealing significant separations in both the microbial richness (alpha diversity) and composition (beta diversity). Our study highlights the presence of significant differences among these procedures, probably due to the different DNA extracting capacities and to the different properties of the milk samples, revealing that the selection of DNA extraction protocol is a critical point. IMPORTANCE The advance of high-throughput technologies has increased our knowledge of the world of microorganisms, especially of microbial populations inhabiting living animals. This study provides evidence that milk, as other complex sources, could be critical for generating high-quality DNA for microbiota analysis. In addition, it demonstrates that the microbial population highlighted by metagenomic studies changes in relation to different DNA extraction procedures, revealing that attention should be paid especially when comparing different studies.

It is known that data from both 16S and shotgun metagenomics studies are subject to biases that cause the observed relative abundances of taxa to differ from their true values. Model community analyses, in which the relative abundances of all taxa in the sample are known by construction, seem to offer the hope that these biases can be measured. However, it is unclear whether the bias we measure in a mock community analysis is the same as we measure in a sample in which taxa are spiked in at known relative abundance, or if the biases we measure in spike-in samples is the same as the bias we would measure in a real (e.g., biological) sample. Here, we consider these questions in the context of 16S rRNA measurements on three sets of samples: the commercially available Zymo cells model community; the Zymo model community mixed with Swedish Snus, a smokeless tobacco product that is virtually bacteria-free; and a set of commercially available smokeless tobacco products. Each set of samples was subject to four different extraction protocols. The goal of our analysis is to determine whether the patterns of bias observed in each set of samples are the same, i.e., can we learn about the bias in the commercially available smokeless tobacco products by studying the Zymo cells model community?

Abstract High‐throughput sequencing has the potential to describe biological communities with high efficiency yet comprehensive assessment of diversity with species‐level resolution remains one of the most challenging aspects of metabarcoding studies. We investigated the utility of curated ribosomal and mitochondrial nematode reference sequence databases for determining phylum‐specific species‐level clustering thresholds. We compiled 438 ribosomal and 290 mitochondrial sequences which identified 99% and 94% as the species delineation clustering threshold, respectively. These thresholds were evaluated in HTS data from mock communities containing 39 nematode species as well as environmental samples from Vietnam. We compared the taxonomic description of the mocks generated by two read‐merging and two clustering algorithms and the cluster‐free Dada2 pipeline. Taxonomic assignment with the RDP classifier was assessed under different training sets. Our results showed that 36/39 mock nematode species were identified across the molecular markers (18S: 32, JB2: 19, JB3: 21) in UClust_ref OTUs at their respective clustering thresholds, outperforming UParse_denovo and the commonly used 97% similarity. Dada2 generated the most realistic number of ASVs (18S: 83, JB2: 75, JB3: 82), collectively identifying 30/39 mock species. The ribosomal marker outperformed the mitochondrial markers in terms of species and genus‐level detections for both OTUs and ASVs. The number of taxonomic assignments of OTUs/ASVs was highest when the smallest reference database containing only nematode sequences was used and when sequences were truncated to the respective amplicon length. Overall, OTUs generated more species‐level detections, which were, however, associated with higher error rates compared to ASVs. Genus‐level assignments using ASVs exhibited higher accuracy and lower error rates compared to species‐level assignments, suggesting that this is the most reliable pipeline for rapid assessment of alpha diversity from environmental samples.

Abstract Background Since the analysis of a large number of metagenomic sequences costs heavy computing resources and takes long time, we examined a selected small part of metagenomic sequences as “sample”s of the entire full sequences, both for a mock community and for 10 different existing metagenomics case studies. A mock community with 10 bacterial strains was prepared, and their mixed genome were sequenced by Hiseq. The hits of BLAST search for reference genome of each strain were counted. Each of 176 different small parts selected from these sequences were also searched by BLAST and their hits were also counted, in order to compare them to the original search results from the full sequences. We also prepared small parts of sequences which were selected from 10 publicly downloadable research data of MG-RAST service, and analyzed these samples with MG-RAST. Results Both the BLAST search tests of the mock community and the results from the publicly downloadable researches of MG-RAST show that sampling an extremely small part from sequence data is useful to estimate brief taxonomic information of the original metagenomic sequences. For 9 cases out of 10, the most annotated classes from the MG-RAST analyses of the selected partial sample sequences are the same as the ones from the originals. Conclusions When a researcher wants to estimate brief information of a metagenome’s taxonomic distribution with less computing resources and within shorter time, the researcher can analyze a selected small part of metagenomic sequences. With this approach, we can also build a strategy to monitor metagenome samples of wider geographic area, more frequently.

Next-generation sequencing has revolutionized our ability to investigate the microbiota composition of diverse and complex environments. However, a number of factors can affect the accuracy of microbial community assessment, such as the DNA extraction method, the hypervariable region of 16S rRNA gene targeted, or the PCR primers used for amplification. The aim of this study was to assess the influence of commercially available DNA extraction kits and different primer pairs to provide a non-biased vision of the composition of bacterial communities present in olive xylem sap. For that purpose, branches from “Picual” and “Arbequina” olive cultivars were used for xylem sap extraction using a Scholander chamber device. The DNA extraction protocol significantly affected xylem sap bacterial community assessment. That resulted in significant differences in alpha (Richness) and beta diversity (UniFrac distances) metrics among DNA extraction protocols, with the 12 DNA extraction kits evaluated being clustered in four groups behaving differently. Although the core number of taxa detected by all DNA extraction kits included four phyla, seven classes, 12 orders, 16 or 21 families, and 12 or 14 genera when using the Greengenes or Silva database for taxonomic assignment, respectively, some taxa, particularly those identified at low frequency, were detected by some DNA extraction kits only. The most accurate depiction of a bacterial mock community artificially inoculated on sap samples was generated when using the PowerPlant DNA extraction kit, the combination of 799F/1193R primers amplifying the hypervariable V5–V7 region, and the Silva 132 database for taxonomic assignment. The DESeq2 analysis displayed significant differences among genera abundance between the different PCR primer pairs tested. Thus, Enterobacter, Granulicatella, Prevotella, and Brevibacterium presented a significant higher abundance in all PCR protocols when compared with primer pair 799F/1193R, while the opposite was true for Pseudomonas and Pectobacterium. The methodological approach followed in this study can be useful to optimize plant-associated microbiome analysis, especially when exploring new plant niches. Some of the DNA extraction kits and PCR primers selected in this study will contribute to better characterize bacterial communities inhabiting the xylem sap of olives or other woody crop species.

Next-generation sequencing has revolutionized our ability to investigate the microbiota composition of diverse and complex environments. However, a number of factors can affect the accuracy of microbial community assessment, such as the DNA extraction method, the hypervariable region of 16S rRNA gene targeted, or the PCR primers used for amplification. The aim of this study was to assess the influence of commercially available DNA extraction kits and different primer pairs to provide a non-biased vision of the composition of bacterial communities present in olive xylem sap. For that purpose, branches from "Picual" and "Arbequina" olive cultivars were used for xylem sap extraction using a Scholander chamber device. The DNA extraction protocol significantly affected xylem sap bacterial community assessment. That resulted in significant differences in alpha (Richness) and beta diversity (UniFrac distances) metrics among DNA extraction protocols, with the 12 DNA extraction kits evaluated being clustered in four groups behaving differently. Although the core number of taxa detected by all DNA extraction kits included four phyla, seven classes, 12 orders, 16 or 21 families, and 12 or 14 genera when using the Greengenes or Silva database for taxonomic assignment, respectively, some taxa, particularly those identified at low frequency, were detected by some DNA extraction kits only. The most accurate depiction of a bacterial mock community artificially inoculated on sap samples was generated when using the PowerPlant DNA extraction kit, the combination of 799F/1193R primers amplifying the hypervariable V5–V7 region, and the Silva 132 database for taxonomic assignment. The DESeq2 analysis displayed significant differences among genera abundance between the different PCR primer pairs tested. Thus, Enterobacter , Granulicatella , Prevotella , and Brevibacterium presented a significant higher abundance in all PCR protocols when compared with primer pair 799F/1193R, while the opposite was true for Pseudomonas and Pectobacterium. The methodological approach followed in this study can be useful to optimize plant-associated microbiome analysis, especially when exploring new plant niches. Some of the DNA extraction kits and PCR primers selected in this study will contribute to better characterize bacterial communities inhabiting the xylem sap of olives or other woody crop species. [ABSTRACT FROM AUTHOR]

The gut mycobiome plays an important role in the health and disease of the human gut, but its exact function is still under investigation. While there is a wealth of information available on the bacterial community of the human gut microbiome, research on the fungal community is still relatively limited. In particular, technical methodologies for mycobiome analysis, especially the DNA extraction method for human faecal samples, varied in different studies. In the current study, two commercial kits commonly used in DNA extraction, the QIAamp® Fast DNA Stool Mini Kit and DNeasy PowerSoil Pro Kit, and one manual method, the International Human Microbiome Standards Protocol Q, were compared. Furthermore, the effectiveness of two different bead-beating machines, the Mini-Beadbeater-16 and FastPrep-24 TM 5G, was compared in parallel. A mock fungal community with a known composition of fungal strains was also generated and included to compare different DNA extraction methods. Our results suggested that the method using the DNeasy PowerSoil Pro Kit and Mini-Beadbeater-16 provides the best results to extract DNA from human faecal samples. Based on our data, we propose a standard operating procedure for DNA extraction from human faecal samples for mycobiome analysis., Competing Interests: We declare we have no competing interests., (© 2024 The Authors.)

Summary: Recent studies have questioned the use of high‐throughput sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region to derive a semi‐quantitative representation of fungal community composition. However, comprehensive studies that quantify biases occurring during PCR and sequencing of ITS amplicons are still lacking.We used artificially assembled communities consisting of 10 ITS‐like fragments of varying lengths and guanine‐cytosine (GC) contents to evaluate and quantify biases during PCR and sequencing with Illumina MiSeq, PacBio RS II and PacBio Sequel I technologies.Fragment length variation was the main source of bias in observed community composition relative to the template, with longer fragments generally being under‐represented for all sequencing platforms. This bias was three times higher for Illumina MiSeq than for PacBio RS II and Sequel I. All 10 fragments in the artificial community were recovered when sequenced with PacBio technologies, whereas the three longest fragments (> 447 bases) were lost when sequenced with Illumina MiSeq. Fragment length bias also increased linearly with increasing number of PCR cycles but could be mitigated by optimization of the PCR setup. No significant biases related to GC content were observed.Despite lower sequencing output, PacBio sequencing was better able to reflect the community composition of the template than Illumina MiSeq sequencing. [ABSTRACT FROM AUTHOR]

Abstract Background Advances in sequencing technologies and bioinformatics have made the analysis of microbial communities almost routine. Nonetheless, the need remains to improve on the techniques used for gathering such data, including increasing throughput while lowering cost and benchmarking the techniques so that potential sources of bias can be better characterized. Methods We present a triple-index amplicon sequencing strategy to sequence large numbers of samples at significantly lower c ost and in a shorter timeframe compared to existing methods. The design employs a two-stage PCR protocol, incorpo rating three barcodes to each sample, with the possibility to add a fourth-index. It also includes heterogeneity spacers to overcome low complexity issues faced when sequencing amplicons on Illumina platforms. Results The library preparation method was extensively benchmarked through analysis of a mock community in order to assess biases introduced by sample indexing, number of PCR cycles, and template concentration. We further evaluated the method through re-sequencing of a standardized environmental sample. Finally, we evaluated our protocol on a set of fecal samples from a small cohort of healthy adults, demonstrating good performance in a realistic experimental setting. Between-sample variation was mainly related to batch effects, such as DNA extraction, while sample indexing was also a significant source of bias. PCR cycle number strongly influenced chimera formation and affected relative abundance estimates of species with high GC content. Libraries were sequenced using the Illumina HiSeq and MiSeq platforms to demonstrate that this protocol is highly scalable to sequence thousands of samples at a very low cost. Conclusions Here, we provide the most comprehensive study of performance and bias inherent to a 16S rRNA gene amplicon sequencing method to date. Triple-indexing greatly reduces the number of long custom DNA oligos required for library preparation, while the inclusion of variable length heterogeneity spacers minimizes the need for PhiX spike-in. This design results in a significant cost reduction of highly multiplexed amplicon sequencing. The biases we characterize highlight the need for highly standardized protocols. Reassuringly, we find that the biological signal is a far stronger structuring factor than the various sources of bias.

Abstract Background Metagenomics is the study of the microbial genomes isolated from communities found on our bodies or in our environment. By correctly determining the relation between human health and the human associated microbial communities, novel mechanisms of health and disease can be found, thus enabling the development of novel diagnostics and therapeutics. Due to the diversity of the microbial communities, strategies developed for aligning human genomes cannot be utilized, and genomes of the microbial species in the community must be assembled de novo. However, in order to obtain the best metagenomic assemblies, it is important to choose the proper assembler. Due to the rapidly evolving nature of metagenomics, new assemblers are constantly created, and the field has not yet agreed on a standardized process. Furthermore, the truth sets used to compare these methods are either too simple (computationally derived diverse communities) or complex (microbial communities of unknown composition), yielding results that are hard to interpret. In this analysis, we interrogate the strengths and weaknesses of five popular assemblers through the use of defined biological samples of known genomic composition and abundance. We assessed the performance of each assembler on their ability to reassemble genomes, call taxonomic abundances, and recreate open reading frames (ORFs). Results We tested five metagenomic assemblers: Omega, metaSPAdes, IDBA-UD, metaVelvet and MEGAHIT on known and synthetic metagenomic data sets. MetaSPAdes excelled in diverse sets, IDBA-UD performed well all around, metaVelvet had high accuracy in high abundance organisms, and MEGAHIT was able to accurately differentiate similar organisms within a community. At the ORF level, metaSPAdes and MEGAHIT had the least number of missing ORFs within diverse and similar communities respectively. Conclusions Depending on the metagenomics question asked, the correct assembler for the task at hand will differ. It is important to choose the appropriate assembler, and thus clearly define the biological problem of an experiment, as different assemblers will give different answers to the same question.

In food safety the detection of food contaminations with pathogenic microorganisms is a race against time and often outpaced by error-prone epidemiological approaches. For evidence-based outbreak investigations fast and reliable techniques and procedures are required to identify the source of infection. Metagenomics has the potential to become a powerful tool in the field of modern food safety, since it allows the detection, identification and characterization of a broad range of pathogens in a single experiment without pre-cultivation within a couple of days. Nevertheless, sample handling, sequencing and data analysis are challenging and can introduce errors and biases into the analysis. In order to evaluate the potential of metagenomics in food safety, we generated a mock community containing DNA of foodborne bacteria. Herewith, we compare the aptitude of the two prevalent approaches – 16S rDNA amplicon sequencing and whole genome shotgun sequencing – for the detection of foodborne bacteria using different parameters during sample preparation, sequencing and data analysis. 16S rDNA sequencing did not only result in high deviations from the expected sample composition on genus and species level, but more importantly lacked the detection of several pathogenic species. While shotgun sequencing is more suitable for species detection, abundance estimation, genome assembly and species characterization, the performance can vary depending on the library preparation kit, which was confirmed for a naturally Francisella tularensis contaminated game meat sample. The application of the Nextera XT DNA Library Preparation Kit for shotgun sequencing did not only result in lower reference genome recovery and coverage, but also in distortions of the mock community composition. For data analysis, we propose a publicly available workflow for pathogen detection and characterization and demonstrate its benefits on the usability of metagenomic sequencing in food safety by analyzing an authentic metagenomic sample.

The amount of host DNA poses a major challenge to metagenome analysis. However, there is no guidance on the levels of host DNA, nor on the depth of sequencing needed to acquire meaningful information from whole metagenome sequencing (WMS). Here, we evaluated the impact of a wide range of amounts of host DNA and sequencing depths on microbiome taxonomic profiling using WMS. Synthetic samples with increasing levels of host DNA were created by spiking DNA of a mock bacterial community, with DNA from a mouse-derived cell line. Taxonomic analysis revealed that increasing proportions of host DNA led to decreased sensitivity in detecting very low and low abundant species. Reduction of sequencing depth had major impact on the sensitivity of WMS for profiling samples with 90% host DNA, increasing the number of undetected species. Finally, analysis of simulated datasets with fixed depth of 10 million reads confirmed that microbiome profiling becomes more inaccurate as the level of host DNA increases in a sample. In conclusion, samples with high amounts of host DNA coupled with reduced sequencing depths, decrease WMS coverage for characterization of the microbiome. This study highlights the importance of carefully considering these aspects in the design of WMS experiments to maximize microbiome analyses.

Invasive fungal infections are an increasingly important cause of human morbidity and mortality. We generated a next-generation sequencing (NGS)-based method designed to detect a wide range of fungi and applied it to analysis of the fungal microbiome (mycobiome) of the lung during fungal infection. Internal transcribed spacer 1 (ITS1) amplicon sequencing and a custom analysis pipeline detected 96% of species from three mock communities comprised of potential fungal lung pathogens with good recapitulation of the expected species distributions (Pearson correlation coefficients r = 0.63, p = 0.004; r = 0.71, p < 0.001; r = 0.62, p = 0.002). We used this pipeline to analyze mycobiomes of bronchoalveolar lavage (BAL) specimens classified as culture-negative (n = 50) or culture-positive (n = 39) for Blastomyces dermatitidis/gilchristii, the causative agent of North America blastomycosis. Detected in 91.4% of the culture-positive samples, Blastomyces dominated (>50% relative abundance) the mycobiome in 68.6% of these culture-positive samples but was absent in culture-negative samples. To overcome any bias in relative abundance due to between-sample variation in fungal biomass, an abundance-weighting calculation was used to normalize the data by accounting for sample-specific PCR cycle number and PCR product concentration data utilized during sample preparation. After normalization, there was a statistically significant greater overall abundance of ITS1 amplicon in the Blastomyces-culture-positive samples versus culture-negative samples. Moreover, the normalization revealed a greater biomass of yeast and environmental fungi in several Blastomyces-culture-positive samples than in the culture-negative samples. Successful detection of Coccidioides, Scedosporium, Phaeoacremonium, and Aspergillus in 6 additional culture-positive BALs by ITS1 amplicon sequencing demonstrates the ability of this method to detect a broad range of fungi from clinical specimens, suggesting that it may be a potentially useful adjunct to traditional fungal microbiological testing for the diagnosis of respiratory mycoses.

Metaproteomics, the study of protein expression in microbial communities, is a versatile tool for environmental microbiology. Achieving sufficiently high metaproteome coverage to obtain a comprehensive picture of the activities and interactions in microbial communities is one of the current challenges in metaproteomics. An essential step to maximize the number of identified proteins is peptide separation via liquid chromatography (LC) prior to mass spectrometry (MS). Thorough optimization and comparison of LC methods for metaproteomics are, however, currently lacking. Here, we present an extensive development and test of different 1D and 2D-LC approaches for metaproteomic peptide separations. We used fully characterized mock community samples to evaluate metaproteomic approaches with very long analytical columns (50 and 75 cm) and long gradients (up to 12 h). We assessed a total of over 20 different 1D and 2D-LC approaches in terms of number of protein groups and unique peptides identified, peptide spectrum matches (PSMs) generated, the ability to detect proteins of low-abundance species, the effect of technical replicate runs on protein identifications and method reproducibility. We show here that, while 1D-LC approaches are faster and easier to set up and lead to more identifications per minute of runtime, 2D-LC approaches allow for a higher overall number of identifications with up to >10,000 protein groups identified. We also compared the 1D and 2D-LC approaches to a standard GeLC workflow, in which proteins are pre-fractionated via gel electrophoresis. This method yielded results comparable to the 2D-LC approaches, however with the drawback of a much increased sample preparation time. Based on our results, we provide recommendations on how to choose the best LC approach for metaproteomics experiments, depending on the study aims.

Receiving the same results from repeated analysis of the same sample is a basic principle in science. The inability to reproduce previously published results has led to discussions of a reproducibility crisis within science. For studies of microbial communities, the problem of reproducibility is more pronounced and has, in some fields, led to a discussion on the very existence of a constantly present microbiota. In this study, DNA from 44 bovine milk samples were extracted twice and the V3-V4 region of the 16S rRNA gene was sequenced in two separate runs. The FASTQ files from the two data sets were run through the same bioinformatics pipeline using the same settings and results from the two data sets were compared. Milk samples collected maximally 2 h apart were used as replicates and permitted comparisons to be made within the same run. Results show a significant difference in species richness between the two sequencing runs although Shannon and Simpson's diversity was the same. Multivariate analyses of all samples demonstrate that the sequencing run was a driver for variation. Direct comparison of similarity between samples and sequencing run showed an average similarity of 42%-45% depending on whether binary or abundance-based similarity indices were used. Within-run comparisons of milk samples collected maximally 2 h apart showed an average similarity of 39%-47% depending on the similarity index used and that similarity differed significantly between runs. We conclude that repeated DNA extraction and sequencing significantly can affect the results of a low microbial biomass microbiota study., (© 2023 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd.)

In food safety the detection of food contaminations with pathogenic microorganisms is a race against time and often outpaced by error-prone epidemiological approaches. For evidence-based outbreak investigations fast and reliable techniques and procedures are required to identify the source of infection. Metagenomics has the potential to become a powerful tool in the field of modern food safety, since it allows the detection, identification and characterization of a broad range of pathogens in a single experiment without pre-cultivation within a couple of days. Nevertheless, sample handling, sequencing and data analysis are challenging and can introduce errors and biases into the analysis. In order to evaluate the potential of metagenomics in food safety, we generated a mock community containing DNA of foodborne bacteria. Herewith, we compare the aptitude of the two prevalent approaches – 16S rDNA amplicon sequencing and whole genome shotgun sequencing – for the detection of foodborne bacteria using different parameters during sample preparation, sequencing and data analysis. 16S rDNA sequencing did not only result in high deviations from the expected sample composition on genus and species level, but more importantly lacked the detection of several pathogenic species. While shotgun sequencing is more suitable for species detection, abundance estimation, genome assembly and species characterization, the performance can vary depending on the library preparation kit, which was confirmed for a naturally Francisella tularensis contaminated game meat sample. The application of the Nextera XT DNA Library Preparation Kit for shotgun sequencing did not only result in lower reference genome recovery and coverage, but also in distortions of the mock community composition. For data analysis, we propose a publicly available workflow for pathogen detection and characterization and demonstrate its benefits on the usability of metagenomic sequencing in food safety by analyzing an authentic metagenomic sample. [ABSTRACT FROM AUTHOR]

The amount of host DNA poses a major challenge to metagenome analysis. However, there is no guidance on the levels of host DNA, nor on the depth of sequencing needed to acquire meaningful information from whole metagenome sequencing (WMS). Here, we evaluated the impact of a wide range of amounts of host DNA and sequencing depths on microbiome taxonomic profiling using WMS. Synthetic samples with increasing levels of host DNA were created by spiking DNA of a mock bacterial community, with DNA from a mouse-derived cell line. Taxonomic analysis revealed that increasing proportions of host DNA led to decreased sensitivity in detecting very low and low abundant species. Reduction of sequencing depth had major impact on the sensitivity of WMS for profiling samples with 90% host DNA, increasing the number of undetected species. Finally, analysis of simulated datasets with fixed depth of 10 million reads confirmed that microbiome profiling becomes more inaccurate as the level of host DNA increases in a sample. In conclusion, samples with high amounts of host DNA coupled with reduced sequencing depths, decrease WMS coverage for characterization of the microbiome. This study highlights the importance of carefully considering these aspects in the design of WMS experiments to maximize microbiome analyses. [ABSTRACT FROM AUTHOR]

Invasive fungal infections are an increasingly important cause of human morbidity and mortality. We generated a next-generation sequencing (NGS)-based method designed to detect a wide range of fungi and applied it to analysis of the fungal microbiome (mycobiome) of the lung during fungal infection. Internal transcribed spacer 1 (ITS1) amplicon sequencing and a custom analysis pipeline detected 96% of species from three mock communities comprised of potential fungal lung pathogens with good recapitulation of the expected species distributions (Pearson correlation coefficients r = 0.63, p = 0.004; r = 0.71, p < 0.001; r = 0.62, p = 0.002). We used this pipeline to analyze mycobiomes of bronchoalveolar lavage (BAL) specimens classified as culture-negative (n = 50) or culture-positive (n = 39) for Blastomyces dermatitidis/gilchristii , the causative agent of North America blastomycosis. Detected in 91.4% of the culture-positive samples, Blastomyces dominated (>50% relative abundance) the mycobiome in 68.6% of these culture-positive samples but was absent in culture-negative samples. To overcome any bias in relative abundance due to between-sample variation in fungal biomass, an abundance-weighting calculation was used to normalize the data by accounting for sample-specific PCR cycle number and PCR product concentration data utilized during sample preparation. After normalization, there was a statistically significant greater overall abundance of ITS1 amplicon in the Blastomyces -culture-positive samples versus culture-negative samples. Moreover, the normalization revealed a greater biomass of yeast and environmental fungi in several Blastomyces -culture-positive samples than in the culture-negative samples. Successful detection of Coccidioides , Scedosporium , Phaeoacremonium , and Aspergillus in 6 additional culture-positive BALs by ITS1 amplicon sequencing demonstrates the ability of this method to detect a broad range of fungi from clinical specimens, suggesting that it may be a potentially useful adjunct to traditional fungal microbiological testing for the diagnosis of respiratory mycoses. [ABSTRACT FROM AUTHOR]

Metaproteomics, the study of protein expression in microbial communities, is a versatile tool for environmental microbiology. Achieving sufficiently high metaproteome coverage to obtain a comprehensive picture of the activities and interactions in microbial communities is one of the current challenges in metaproteomics. An essential step to maximize the number of identified proteins is peptide separation via liquid chromatography (LC) prior to mass spectrometry (MS). Thorough optimization and comparison of LC methods for metaproteomics are, however, currently lacking. Here, we present an extensive development and test of different 1D and 2D-LC approaches for metaproteomic peptide separations. We used fully characterized mock community samples to evaluate metaproteomic approaches with very long analytical columns (50 and 75 cm) and long gradients (up to 12 h). We assessed a total of over 20 different 1D and 2D-LC approaches in terms of number of protein groups and unique peptides identified, peptide spectrum matches (PSMs) generated, the ability to detect proteins of low-abundance species, the effect of technical replicate runs on protein identifications and method reproducibility. We show here that, while 1D-LC approaches are faster and easier to set up and lead to more identifications per minute of runtime, 2D-LC approaches allow for a higher overall number of identifications with up to >10,000 protein groups identified. We also compared the 1D and 2D-LC approaches to a standard GeLC workflow, in which proteins are pre-fractionated via gel electrophoresis. This method yielded results comparable to the 2D-LC approaches, however with the drawback of a much increased sample preparation time. Based on our results, we provide recommendations on how to choose the best LC approach for metaproteomics experiments, depending on the study aims. [ABSTRACT FROM AUTHOR]

High‐throughput sequencing has the potential to describe biological communities with high efficiency yet comprehensive assessment of diversity with species‐level resolution remains one of the most challenging aspects of metabarcoding studies. We investigated the utility of curated ribosomal and mitochondrial nematode reference sequence databases for determining phylum‐specific species‐level clustering thresholds. We compiled 438 ribosomal and 290 mitochondrial sequences which identified 99% and 94% as the species delineation clustering threshold, respectively. These thresholds were evaluated in HTS data from mock communities containing 39 nematode species as well as environmental samples from Vietnam. We compared the taxonomic description of the mocks generated by two read‐merging and two clustering algorithms and the cluster‐free Dada2 pipeline. Taxonomic assignment with the RDP classifier was assessed under different training sets. Our results showed that 36/39 mock nematode species were identified across the molecular markers (18S: 32, JB2: 19, JB3: 21) in UClust_ref OTUs at their respective clustering thresholds, outperforming UParse_denovo and the commonly used 97% similarity. Dada2 generated the most realistic number of ASVs (18S: 83, JB2: 75, JB3: 82), collectively identifying 30/39 mock species. The ribosomal marker outperformed the mitochondrial markers in terms of species and genus‐level detections for both OTUs and ASVs. The number of taxonomic assignments of OTUs/ASVs was highest when the smallest reference database containing only nematode sequences was used and when sequences were truncated to the respective amplicon length. Overall, OTUs generated more species‐level detections, which were, however, associated with higher error rates compared to ASVs. Genus‐level assignments using ASVs exhibited higher accuracy and lower error rates compared to species‐level assignments, suggesting that this is the most reliable pipeline for rapid assessment of alpha diversity from environmental samples. Curated DNA barcode data were used to delineate nematode‐specific inter‐specific operational taxonomic unit (OTUs) clustering thresholds, which were subsequently tested on a artificial nematode community and environmental samples. Open‐reference clustering (UClust_ref) at these empirically derived thresholds outperformed denovo (UParse) as well as cluster‐free approaches (DADA2), resulting in the most accurate description of the nematode community. [ABSTRACT FROM AUTHOR]

Background: Since the analysis of a large number of metagenomic sequences costs heavy computing resources and takes long time, we examined a selected small part of metagenomic sequences as "sample"s of the entire full sequences, both for a mock community and for 10 different existing metagenomics case studies. A mock community with 10 bacterial strains was prepared, and their mixed genome were sequenced by Hiseq. The hits of BLAST search for reference genome of each strain were counted. Each of 176 different small parts selected from these sequences were also searched by BLAST and their hits were also counted, in order to compare them to the original search results from the full sequences. We also prepared small parts of sequences which were selected from 10 publicly downloadable research data of MG-RAST service, and analyzed these samples with MG-RAST. Results: Both the BLAST search tests of the mock community and the results from the publicly downloadable researches of MG-RAST show that sampling an extremely small part from sequence data is useful to estimate brief taxonomic information of the original metagenomic sequences. For 9 cases out of 10, the most annotated classes from the MG-RAST analyses of the selected partial sample sequences are the same as the ones from the originals. Conclusions: When a researcher wants to estimate brief information of a metagenome's taxonomic distribution with less computing resources and within shorter time, the researcher can analyze a selected small part of metagenomic sequences. With this approach, we can also build a strategy to monitor metagenome samples of wider geographic area, more frequently. [ABSTRACT FROM AUTHOR]

The European ash (Fraxinus excelsior) is threatened by the introduced ascomycete Hymenoscyphus fraxineus , the causal agent of ash dieback. Endophytic fungi are known to modulate their host's resistance against pathogens. To understand possible consequences of ash dieback on the endophytic mycobiome, F. excelsior leaves were collected in naturally regenerated forests and the fungal communities analyzed by classic culture and Illumina amplicon sequencing using a newly developed and validated fungal-specific primer. Collections were done in the area infested by ash dieback north of the Alps, and in the disease free area on the south side. Sycamore maple (Acer pseudoplatanus) was additionally collected, as well as the flowering ash (F. ornus), which occurs naturally in the south and shows tolerance to ash dieback. Both cultivation and amplicon sequencing revealed characteristic endophytic fungal communities dominated by several strictly host specific Venturia species. On A. pseudoplatanus , a hitherto undescribed Venturia species was identified. Due to its dominance on F. excelsior , V. fraxini is unlikely to go extinct in case of reduced host densities. A majority of species was not strictly host specific and is therefore likely less affected by ash dieback in the future. Still, shifts in community structure and loss of genetic diversity cannot be excluded. The potentially endangered endophyte Hymenoscyphus albidus was rarely found. In addition to host specificity, species with preferences for leaf laminae or petioles were found. We also detected considerable geographical variation between sampling sites and clear differences between the two sides of the Alps for endophytes of F. excelsior , but not A. pseudoplatanus. Since sycamore maple is not affected by an epidemic, this could point toward an influence of ash dieback on ash communities, although firm conclusions are not possible because of host preferences and climatic differences. Furthermore, the mycobiota of F. excelsior trees with or without dieback symptoms were compared, but no clear differences were detected. Besides methodical refinement, our study provides comprehensive data on the ash mycobiome that we expect to be subject to changes caused by an emerging disease of the host tree. [ABSTRACT FROM AUTHOR]

ABSTRACT Mock communities have been used in microbiome method development to help estimate biases introduced in PCR amplification and sequencing and to optimize pipeline outputs. Nevertheless, the strong value of routine mock community analysis beyond initial method development is rarely, if ever, considered. Here we report that our routine use of mock communities as internal standards allowed us to discover highly aberrant and strong biases in the relative proportions of multiple taxa in a single Illumina HiSeqPE250 run. In this run, an important archaeal taxon virtually disappeared from all samples, and other mock community taxa showed >2-fold high or low abundance, whereas a rerun of those identical amplicons (from the same reaction tubes) on a different date yielded “normal” results. Although obvious from the strange mock community results, we could have easily missed the problem had we not used the mock communities because of natural variation of microbiomes at our site. The “normal” results were validated over four MiSeqPE300 runs and three HiSeqPE250 runs, and run-to-run variation was usually low. While validating these “normal” results, we also discovered that some mock microbial taxa had relatively modest, but consistent, differences between sequencing platforms. We strongly advise the use of mock communities in every sequencing run to distinguish potentially serious aberrations from natural variations. The mock communities should have more than just a few members and ideally at least partly represent the samples being analyzed to detect problems that show up only in some taxa and also to help validate clustering. IMPORTANCE Despite the routine use of standards and blanks in virtually all chemical or physical assays and most biological studies (a kind of “control”), microbiome analysis has traditionally lacked such standards. Here we show that unexpected problems of unknown origin can occur in such sequencing runs and yield completely incorrect results that would not necessarily be detected without the use of standards. Assuming that the microbiome sequencing analysis works properly every time risks serious errors that can be detected by the use of mock communities.

The rumen microbiome scientific community has utilized amplicon sequencing as an aid in identifying potential community compositional trends that could be used as an estimation of various production and performance traits including methane emission, animal protein production efficiency, and ruminant health status. In order to translate rumen microbiome studies into executable application, there is a need for experimental and analytical concordance within the community. The objective of this study was to assess these factors in relation to selected currently established methods for 16S phylogenetic community analysis on a microbial community standard (MC) and a DNA standard (DS; ZymoBIOMICSTM). DNA was extracted from MC using the RBBC method commonly used for microbial DNA extraction from rumen digesta samples. 16S rRNA amplicon libraries were generated for the MC and DS using primers routinely used for rumen bacterial and archaeal community analysis. The primers targeted the V4 and V3–V4 region of the 16S rRNA gene and samples were subjected to both 20 and 28 polymerase chain reaction (PCR) cycles under identical cycle conditions. Sequencing was conducted using the Illumina MiSeq platform. As the bacteria contained in the microbial mock community were well-classified species, and for ease of explanation, we used the results of the Basic Local Alignment Search Tool classification to assess the DNA, PCR cycle number, and primer type. Sequence classification methodology was assessed independently. Spearman’s correlation analysis indicated that utilizing the repeated bead beating and column method for DNA extraction in combination with primers targeting the 16S rRNA gene using 20 first-round PCR cycles was sufficient for amplicon sequencing to generate a relatively accurate depiction of the bacterial communities present in rumen samples. These results also emphasize the requirement to develop and utilize positive mock community controls for all rumen microbiomic studies in order to discern errors which may arise at any step during a next-generation sequencing protocol.

Profiling of microbial community composition is frequently performed by partial 16S rRNA gene sequencing on benchtop platforms following PCR amplification of specific hypervariable regions within this gene. Accuracy and reproducibility of this strategy are two key parameters to consider, which may be influenced during all processes from sample collection and storage, through DNA extraction and PCR based library preparation to the final sequencing. In order to evaluate both the reproducibility and accuracy of 16S rRNA gene based microbial profiling using the Ion Torrent PGM platform, we prepared libraries and performed sequencing of a well-defined and validated 20-member bacterial DNA mock community on five separate occasions and compared results with the expected even distribution. In general the applied method had a median coefficient of variance of 11.8% (range 5.5–73.7%) for all 20 included strains in the mock community across five separate sequencing runs, with underrepresented strains generally showing the largest degree of variation. In terms of accuracy, mock community species belonging to Proteobacteria were underestimated, whereas those belonging to Firmicutes were mostly overestimated. This could be explained partly by premature read truncation, but to larger degree their genomic GC-content, which correlated negatively with the observed relative abundances, suggesting a PCR bias against GC-rich species during library preparation. Increasing the initial denaturation time during the PCR amplification from 30 to 120 s resulted in an increased average relative abundance of the three mock community members with the highest genomic GC%, but did not significantly change the overall evenness of the community distribution. Therefore, efforts should be made to optimize the PCR conditions prior to sequencing in order to maximize accuracy.

High-depth sequencing of universal marker genes such as the 16S rRNA gene is a common strategy to profile microbial communities. Traditionally, sequence reads are clustered into operational taxonomic units (OTUs) at a defined identity threshold to avoid sequencing errors generating spurious taxonomic units. However, there have been numerous bioinformatic packages recently released that attempt to correct sequencing errors to determine real biological sequences at single nucleotide resolution by generating amplicon sequence variants (ASVs). As more researchers begin to use high resolution ASVs, there is a need for an in-depth and unbiased comparison of these novel ''denoising'' pipelines. In this study, we conduct a thorough comparison of three of the most widely-used denoising packages (DADA2, UNOISE3, and Deblur) as well as an open-reference 97% OTU clustering pipeline on mock, soil, and host-associated communities. We found from the mock community analyses that although they produced similar microbial compositions based on relative abundance, the approaches identified vastly different numbers of ASVs that significantly impact alpha diversity metrics. Our analysis on real datasets using recommended settings for each denoising pipeline also showed that the three packages were consistent in their per-sample compositions, resulting in only minor differences based on weighted UniFrac and Bray Curtis dissimilarity. DADA2 tended to find more ASVs than the other two denoising pipelines when analyzing both the real soil data and two other host-associated datasets, suggesting that it could be better at finding rare organisms, but at the expense of possible false positives. The open-reference OTU clustering approach identified considerably more OTUs in comparison to the number of ASVs from the denoising pipelines in all datasets tested. The three denoising approaches were significantly different in their run times, with UNOISE3 running greater than 1,200 and 15 times faster than DADA2 and Deblur, respectively. Our findings indicate that, although all pipelines result in similar general community structure, the number of ASVs/OTUs and resulting alpha-diversity metrics varies considerably and should be considered when attempting to identify rare organisms from possible background noise. [ABSTRACT FROM AUTHOR]

The rumen microbiome scientific community has utilized amplicon sequencing as an aid in identifying potential community compositional trends that could be used as an estimation of various production and performance traits including methane emission, animal protein production efficiency, and ruminant health status. In order to translate rumen microbiome studies into executable application, there is a need for experimental and analytical concordance within the community. The objective of this study was to assess these factors in relation to selected currently established methods for 16S phylogenetic community analysis on a microbial community standard (MC) and a DNA standard (DS; ZymoBIOMICSTM). DNA was extracted from MC using the RBBC method commonly used for microbial DNA extraction from rumen digesta samples. 16S rRNA amplicon libraries were generated for the MC and DS using primers routinely used for rumen bacterial and archaeal community analysis. The primers targeted the V4 and V3–V4 region of the 16S rRNA gene and samples were subjected to both 20 and 28 polymerase chain reaction (PCR) cycles under identical cycle conditions. Sequencing was conducted using the Illumina MiSeq platform. As the bacteria contained in the microbial mock community were well-classified species, and for ease of explanation, we used the results of the Basic Local Alignment Search Tool classification to assess the DNA, PCR cycle number, and primer type. Sequence classification methodology was assessed independently. Spearman’s correlation analysis indicated that utilizing the repeated bead beating and column method for DNA extraction in combination with primers targeting the 16S rRNA gene using 20 first-round PCR cycles was sufficient for amplicon sequencing to generate a relatively accurate depiction of the bacterial communities present in rumen samples. These results also emphasize the requirement to develop and utilize positive mock community controls for all rumen microbiomic studies in order to discern errors which may arise at any step during a next-generation sequencing protocol. [ABSTRACT FROM AUTHOR]

Profiling of microbial community composition is frequently performed by partial 16S rRNA gene sequencing on benchtop platforms following PCR amplification of specific hypervariable regions within this gene. Accuracy and reproducibility of this strategy are two key parameters to consider, which may be influenced during all processes from sample collection and storage, through DNA extraction and PCR based library preparation to the final sequencing. In order to evaluate both the reproducibility and accuracy of 16S rRNA gene based microbial profiling using the Ion Torrent PGM platform, we prepared libraries and performed sequencing of a well-defined and validated 20-member bacterial DNA mock community on five separate occasions and compared results with the expected even distribution. In general the applied method had a median coefficient of variance of 11.8% (range 5.5-73.7%) for all 20 included strains in the mock community across five separate sequencing runs, with underrepresented strains generally showing the largest degree of variation. In terms of accuracy, mock community species belonging to Proteobacteria were underestimated, whereas those belonging to Firmicutes were mostly overestimated. This could be explained partly by premature read truncation, but to larger degree their genomic GC-content, which correlated negatively with the observed relative abundances, suggesting a PCR bias against GCrich species during library preparation. Increasing the initial denaturation time during the PCR amplification from 30 to 120 s resulted in an increased average relative abundance of the three mock community members with the highest genomic GC%, but did not significantly change the overall evenness of the community distribution. Therefore, efforts should be made to optimize the PCR conditions prior to sequencing in order to maximize accuracy. [ABSTRACT FROM AUTHOR]

Inclusion and investigation of technical controls in microbiome sequencing studies is important for understanding technical biases and errors. Here, we present chkMocks , a general R-based tool that allows researchers to compare the composition of mock communities that are processed along with samples to their theoretical composition. A visual comparison between experimental and theoretical community composition and their correlation is provided for researchers to assess the quality of their sample processing workflows., Competing Interests: All authors declared that there are no conflicts of interest. While we use ZymoBiomics data, we, the developers of chkMocks, are not associated with the manufacturers and this work should not be considered as an endorsement for the said product., (© The Author(s) 2023.)