Analysis of large 16S rRNA Illumina data sets: Impact of singleton read filtering on microbial community description

Environmental metabolomics characterization of modern stromatolites and annotation of ibhayipeptolides

Lithified layers of complex microbial mats known as microbialites are ubiquitous in the fossil record, and modern forms are increasingly identified globally. A key challenge to developing an understanding of microbialite formation and environmental role is how to investigate complex and diverse communities in situ. We selected living, layered microbialites (stromatolites) in a peritidal environment near Schoenmakerskop, Eastern Cape, South Africa to conduct a spatial survey mapping the composition and small molecule production of the microbial communities from environmental samples. Substrate core samples were collected from nine sampling stations ranging from the upper point of the freshwater inflow to the lower marine interface where tidal overtopping takes place. Substrate cores provided material for parallel analyses of microbial community diversity by 16S rRNA gene amplicon sequencing and metabolomics using LC-MS2. Species and metabolite diversities were correlated, and prominent specialized metabolites were targeted for preliminary characterization. A new series of cyclic hexadepsipeptides, named ibhayipeptolides, was most abundant in substrate cores of submerged microbialites. These results demonstrate the detection and identification of metabolites from mass-limited environmental samples and contribute knowledge about microbialite chemistry and biology, which facilitates future targeted studies of specialized metabolite function and biosynthesis.

Benchmarking MicrobIEM - a user-friendly tool for decontamination of microbiome sequencing data

BACKGROUND

Microbiome analysis is becoming a standard component in many scientific studies, but also requires extensive quality control of the 16S rRNA gene sequencing data prior to analysis. In particular, when investigating low-biomass microbial environments such as human skin, contaminants distort the true microbiome sample composition and need to be removed bioinformatically. We introduce MicrobIEM, a novel tool to bioinformatically remove contaminants using negative controls.

RESULTS

We benchmarked MicrobIEM against five established decontamination approaches in four 16S rRNA amplicon sequencing datasets: three serially diluted mock communities (108-103 cells, 0.4-80% contamination) with even or staggered taxon compositions and a skin microbiome dataset. Results depended strongly on user-selected algorithm parameters. Overall, sample-based algorithms separated mock and contaminant sequences best in the even mock, whereas control-based algorithms performed better in the two staggered mocks, particularly in low-biomass samples (≤ 106 cells). We show that a correct decontamination benchmarking requires realistic staggered mock communities and unbiased evaluation measures such as Youden's index. In the skin dataset, the Decontam prevalence filter and MicrobIEM's ratio filter effectively reduced common contaminants while keeping skin-associated genera.

CONCLUSIONS

MicrobIEM's ratio filter for decontamination performs better or as good as established bioinformatic decontamination tools. In contrast to established tools, MicrobIEM additionally provides interactive plots and supports selecting appropriate filtering parameters via a user-friendly graphical user interface. Therefore, MicrobIEM is the first quality control tool for microbiome experts without coding experience.

Impact of phosphorus limitation on medium-chain-length polyhydroxyalkanoate production by activated sludge

For a sustainable economy, biodegradable biopolymers polyhydroxyalkanoates (PHA) are desirable substitutes to petroleum-based plastics that contaminate our environment. Medium-chain-length (MCL) PHA bioplastics are particularly interesting due to their thermoplastic properties. To hamper the high cost associated to PHA production, the use of bacterial mixed cultures cultivated in open systems and using cheap resources is a promising strategy. Here, we studied the operating conditions favouring direct MCL accumulation by activated sludge, using oleic acid as a model substrate and phosphorus limitation in fed-batch bioreactors. Our results confirm the presence of PHA-accumulating organisms (PHAAO) in activated sludge able to accumulate MCL from oleic acid. A positive correlation between phosphorus (P) limitation and PHA accumulation was demonstrated, allowing up to 26% PHA/total biomass accumulation, and highlighted its negative impact on the MCL/PHA fraction in the polymer. Diversity analysis through 16S rRNA amplicon sequencing revealed a differential selection of PHAAO according to the P-limitation level. A differential behaviour for the orders Pseudomonadales and Burkholderiales at increasing P-limitation levels was revealed, with a higher abundance of the latter at high levels of P-limitation. The PHA accumulation observed in activated sludge open new perspectives for MCL-PHA production system based on P-limitation strategy applied to mixed microbial communities. KEY POINTS: • Direct accumulation of MCL-PHA in activated sludge was demonstrated. • MCL-PHA content is negatively correlated with P-limitation. • Burkholderiales members discriminate the highest P-limitation levels.

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.

Characteristics of gastric cancer gut microbiome according to tumor stage and age segmentation

With the development of 16S rRNA technology, gut microbiome evaluation has been performed in many diseases, including gastrointestinal tumors. Among these cancers, gastric cancer (GC) exhibits high morbidity and mortality and has been extensively studied in its pathogenesis and diagnosis techniques. The current researches have proved that the gut microbiome may have the potential to distinguish GC patients from healthy patients. However, the change of the gut microbiome according to tumor node metastasis classification (TNM) has not been clarified. Besides, the characteristics of gut microbiome in GC patients and their ages of onset are also ambiguous. To address the above shortcomings, we investigated 226 fecal samples and divided them according to their tumor stage and onset age. The findings revealed that surgery and tumor stage can change the characteristic of GC patients' gut microbiota. In specific, the effect of surgery on early gastric cancer (EGC) was greater than that on advanced gastric cancer (AGC), and the comparison of postoperative microflora with healthy people indicated that EGC has more differential bacteria than AGC. Besides, we found that Collinsella, Blautia, Anaerostipes, Dorea, and Lachnospiraceae_ND3007_group expressed differently between EGC and AGC. More importantly, it is the first time revealed that the composition of gut microbiota in GC is different between different onset ages. KEY POINTS: •Gut microbiota of gastric cancer (GC) patients are either highly associated with TNM stage and surgery or not. It shows surgery has more significant changes in early gastric cancer (EGC) than advanced gastric cancer (AGC). •There existed specific gut microbiota between EGC and AGC which may have potential to distinguish the early or advanced GC. •Onset age of GC may influence the gut microbiota: the composition of gut microbiota of early-onset gastric cancer (EOGC) and late-onset gastric cancer (LOGC) is significantly different.

The Role of the Gut Microbiota in the Effects of Early-Life Stress and Dietary Fatty Acids on Later-Life Central and Metabolic Outcomes in Mice

Early-life stress (ELS) leads to increased vulnerability for mental and metabolic disorders. We have previously shown that a low dietary ω-6/ω-3 polyunsaturated fatty acid (PUFA) ratio protects against ELS-induced cognitive impairments. Due to the importance of the gut microbiota as a determinant of long-term health, we here study the impact of ELS and dietary PUFAs on the gut microbiota and how this relates to the previously described cognitive, metabolic, and fatty acid profiles. Male mice were exposed to ELS via the limited bedding and nesting paradigm (postnatal day (P)2 to P9 and to an early diet (P2 to P42) with an either high (15) or low (1) ω-6 linoleic acid to ω-3 alpha-linolenic acid ratio. 16S rRNA was sequenced and analyzed from fecal samples at P21, P42, and P180. Age impacted α- and β-diversity. ELS and diet together predicted variance in microbiota composition and affected the relative abundance of bacterial groups at several taxonomic levels in the short and long term. For example, age increased the abundance of the phyla Bacteroidetes, while it decreased Actinobacteria and Verrucomicrobia; ELS reduced the genera RC9 gut group and Rikenella, and the low ω-6/ω-3 diet reduced the abundance of the Firmicutes Erysipelotrichia. At P42, species abundance correlated with body fat mass and circulating leptin (e.g., Bacteroidetes and Proteobacteria taxa) and fatty acid profiles (e.g., Firmicutes taxa). This study gives novel insights into the impact of age, ELS, and dietary PUFAs on microbiota composition, providing potential targets for noninvasive (nutritional) modulation of ELS-induced deficits. IMPORTANCE Early-life stress (ELS) leads to increased vulnerability to develop mental and metabolic disorders; however, the biological mechanisms leading to such programming are not fully clear. Increased attention has been given to the importance of the gut microbiota as a determinant of long-term health and as a potential target for noninvasive nutritional strategies to protect against the negative impact of ELS. Here, we give novel insights into the complex interaction between ELS, early dietary ω-3 availability, and the gut microbiota across ages and provide new potential targets for (nutritional) modulation of the long-term effects of the early-life environment via the microbiota.

Comparison of traditional and DNA metabarcoding samples for monitoring tropical soil arthropods (Formicidae, Collembola and Isoptera)

The soil fauna of the tropics remains one of the least known components of the biosphere. Long-term monitoring of this fauna is hampered by the lack of taxonomic expertise and funding. These obstacles may potentially be lifted with DNA metabarcoding. To validate this approach, we studied the ants, springtails and termites of 100 paired soil samples from Barro Colorado Island, Panama. The fauna was extracted with Berlese-Tullgren funnels and then either sorted with traditional taxonomy and known, individual DNA barcodes ("traditional samples") or processed with metabarcoding ("metabarcoding samples"). We detected 49 ant, 37 springtail and 34 termite species with 3.46 million reads of the COI gene, at a mean sequence length of 233 bp. Traditional identification yielded 80, 111 and 15 species of ants, springtails and termites, respectively; 98%, 37% and 100% of these species had a Barcode Index Number (BIN) allowing for direct comparison with metabarcoding. Ants were best surveyed through traditional methods, termites were better detected by metabarcoding, and springtails were equally well detected by both techniques. Species richness was underestimated, and faunal composition was different in metabarcoding samples, mostly because 37% of ant species were not detected. The prevalence of species in metabarcoding samples increased with their abundance in traditional samples, and seasonal shifts in species prevalence and faunal composition were similar between traditional and metabarcoding samples. Probable false positive and negative species records were reasonably low (13-18% of common species). We conclude that metabarcoding of samples extracted with Berlese-Tullgren funnels appear suitable for the long-term monitoring of termites and springtails in tropical rainforests. For ants, metabarcoding schemes should be complemented by additional samples of alates from Malaise or light traps.

LANDMark: an ensemble approach to the supervised selection of biomarkers in high-throughput sequencing data

Background

Identification of biomarkers, which are measurable characteristics of biological datasets, can be challenging. Although amplicon sequence variants (ASVs) can be considered potential biomarkers, identifying important ASVs in high-throughput sequencing datasets is challenging. Noise, algorithmic failures to account for specific distributional properties, and feature interactions can complicate the discovery of ASV biomarkers. In addition, these issues can impact the replicability of various models and elevate false-discovery rates. Contemporary machine learning approaches can be leveraged to address these issues. Ensembles of decision trees are particularly effective at classifying the types of data commonly generated in high-throughput sequencing (HTS) studies due to their robustness when the number of features in the training data is orders of magnitude larger than the number of samples. In addition, when combined with appropriate model introspection algorithms, machine learning algorithms can also be used to discover and select potential biomarkers. However, the construction of these models could introduce various biases which potentially obfuscate feature discovery.

Results

We developed a decision tree ensemble, LANDMark, which uses oblique and non-linear cuts at each node. In synthetic and toy tests LANDMark consistently ranked as the best classifier and often outperformed the Random Forest classifier. When trained on the full metabarcoding dataset obtained from Canada’s Wood Buffalo National Park, LANDMark was able to create highly predictive models and achieved an overall balanced accuracy score of 0.96 ± 0.06. The use of recursive feature elimination did not impact LANDMark’s generalization performance and, when trained on data from the BE amplicon, it was able to outperform the Linear Support Vector Machine, Logistic Regression models, and Stochastic Gradient Descent models (p ≤ 0.05). Finally, LANDMark distinguishes itself due to its ability to learn smoother non-linear decision boundaries.

Conclusions

Our work introduces LANDMark, a meta-classifier which blends the characteristics of several machine learning models into a decision tree and ensemble learning framework. To our knowledge, this is the first study to apply this type of ensemble approach to amplicon sequencing data and we have shown that analyzing these datasets using LANDMark can produce highly predictive and consistent models.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04631-z.

Detecting Flavobacterial Fish Pathogens in the Environment via High-Throughput Community Analysis

Diseases caused by the fish pathogens Flavobacterium columnare and Flavobacterium psychrophilum are major contributors of preventable losses in the aquaculture industry. The persistent and difficult-to-control infections caused by these bacteria make timely intervention and prophylactic elimination of pathogen reservoirs important measures to combat these disease-causing agents. In this study, we present two independent assays for detecting these pathogens in a range of environmental samples. Natural water samples were inoculated with F. columnare and F. psychrophilum over 5 orders of magnitude, and pathogen levels were detected using Illumina MiSeq sequencing and droplet digital PCR. Both detection methods accurately identified pathogen-positive samples and showed good agreement in quantifying each pathogen. Additionally, the real-world application of these approaches was demonstrated using environmental samples collected at a rainbow trout (Oncorhynchus mykiss) aquaculture facility. These results show that both methods can serve as useful tools for surveillance efforts in aquaculture facilities, where the early detection of these flavobacterial pathogens may direct preventative measures to reduce disease occurrence. IMPORTANCE Early detection of a deadly disease outbreak in a population can be the difference between mass mortality or mitigated effects. In the present study, we evaluated and compared two molecular techniques for detecting economically impactful aquaculture pathogens. We demonstrate that one of these techniques, 16S rRNA gene sequencing using Illumina MiSeq technology, provides the ability to accurately detect two freshwater fish pathogens, F. columnare and F. psychrophilum, while simultaneously profiling the native microbial community. The second technique, droplet digital PCR, is commonly used for pathogen detection, and the results obtained using the assays we designed with this method served to validate those obtained using the MiSeq method. These two methods offer distinct advantages. The MiSeq method pairs pathogen detection and microbial community profiling to answer immediate and long-term fish health concerns, while the droplet digital PCR method provides fast and highly sensitive detection that is useful for surveillance and rapid clinical responses.

Saccharomyces Cerevisiae Var Boulardii CNCM I-1079 Reduces Expression of Genes Involved in Inflammatory Response in Porcine Cells Challenged by Enterotoxigenic E. Coli and Influences Bacterial Communities in an In Vitro Model of the Weaning Piglet Colon

Enterotoxigenic Escherichia coli (ETEC) is the main infectious agent responsible for piglet post-weaning diarrhea with high mortality rates. Antimicrobials represent the current principal strategy for treating ETEC infections in pig farms, but the occurrence of multi-resistant bacterial strains has considerably increased in the last decades. Thus, finding non-antibiotic alternatives becomes a real emergency. In this context, we investigated the effect of a live yeast strain, Saccharomyces cerevisiae var boulardii CNCM I-1079 (SB) in an in vitro model of the weaning piglet colon implemented with a mucus phase (MPigut-IVM) inoculated with ETEC and coupled with an intestinal porcine cell line IPI-2I. We showed that SB was able to modulate the in vitro microbiota through an increase in Bacteroidiaceae and a decrease in Prevotellaceae families. Effluents collected from the SB treated bioreactors were able to mitigate the expression level of genes encoding non-gel forming mucins, tight junction proteins, innate immune pathway, and pro-inflammatory response in IPI-2I cells. Furthermore, SB exerted a significant protective effect against ETEC adhesion on porcine IPEC-J2 intestinal cells in a dose-dependent manner and showed a positive effect on ETEC-challenged IPEC-J2 by lowering expression of genes involved in pro-inflammatory immune responses. Our results showed that the strain SB CNCM I-1079 could prevent microbiota dysbiosis associated with weaning and protect porcine enterocytes from ETEC infections by reducing bacterial adhesion and modulating the inflammatory response.

Microbial Diversity Profiling of Gut Microbiota of Macropus giganteus Using Three Hypervariable Regions of the Bacterial 16S rRNA

Animal faecal contamination of surface waters poses a human health risk, as they may contain pathogenic bacteria or viruses. Of the numerous animal species residing along surface waterways in Australia, macropod species are a top contributor to wild animals’ faecal pollution load. We characterised the gut microbiota of 30 native Australian Eastern Grey Kangaroos from six geographical regions (five kangaroos from each region) within South East Queensland in order to establish their bacterial diversity and identify potential novel species-specific bacteria for the rapid detection of faecal contamination of surface waters by these animals. Using three hypervariable regions (HVRs) of the 16S rRNA gene (i.e., V1–V3, V3–V4, and V5–V6), for their effectiveness in delineating the gut microbial diversity, faecal samples from each region were pooled and microbial genomic DNA was extracted, sequenced, and analysed. Results indicated that V1-V3 yielded a higher taxa richness due to its larger target region (~480 bp); however, higher levels of unassigned taxa were observed using the V1-V3 region. In contrast, the V3–V4 HVR (~569 bp) attained a higher likelihood of a taxonomic hit identity to the bacterial species level, with a 5-fold decrease in unassigned taxa. There were distinct dissimilarities in beta diversity between the regions, with the V1-V3 region displaying the highest number of unique taxa (n = 42), followed by V3–V4 (n = 11) and V5–V6 (n = 8). Variations in the gut microbial diversity profiles of kangaroos from different regions were also observed, which indicates that environmental factors may impact the microbial development and, thus, the composition of the gut microbiome of these animals.

Pathogen Challenge and Dietary Shift Alter Microbiota Composition and Activity in a Mucin-Associated in vitro Model of the Piglet Colon (MPigut-IVM) Simulating Weaning Transition

Enterotoxigenic Escherichia coli (ETEC) is the principal pathogen responsible for post-weaning diarrhea in newly weaned piglets. Expansion of ETEC at weaning is thought to be the consequence of various stress factors such as transient anorexia, dietary change or increase in intestinal inflammation and permeability, but the exact mechanisms remain to be elucidated. As the use of animal experiments raise more and more ethical concerns, we used a recently developed in vitro model of piglet colonic microbiome and mucobiome, the MPigut-IVM, to evaluate the effects of a simulated weaning transition and pathogen challenge at weaning. Our data suggested that the tested factors impacted the composition and functionality of the MPigut-IVM microbiota. The simulation of weaning transition led to an increase in relative abundance of the Prevotellaceae family which was further promoted by the presence of the ETEC strain. In contrast, several beneficial families such as Bacteroidiaceae or Ruminococcaceae and gut health related short chain fatty acids like butyrate or acetate were reduced upon simulated weaning. Moreover, the incubation of MPigut-IVM filtrated effluents with porcine intestinal cell cultures showed that ETEC challenge in the in vitro model led to an increased expression of pro-inflammatory genes by the porcine cells. This study provides insights about the etiology of a dysbiotic microbiota in post-weaning piglets.

Handling of spurious sequences affects the outcome of high-throughput 16S rRNA gene amplicon profiling

16S rRNA gene amplicon sequencing is a popular approach for studying microbiomes. However, some basic concepts have still not been investigated comprehensively. We studied the occurrence of spurious sequences using defined microbial communities based on data either from the literature or generated in three sequencing facilities and analyzed via both operational taxonomic units (OTUs) and amplicon sequence variants (ASVs) approaches. OTU clustering and singleton removal, a commonly used approach, delivered approximately 50% (mock communities) to 80% (gnotobiotic mice) spurious taxa. The fraction of spurious taxa was generally lower based on ASV analysis, but varied depending on the gene region targeted and the barcoding system used. A relative abundance of 0.25% was found as an effective threshold below which the analysis of spurious taxa can be prevented to a large extent in both OTU- and ASV-based analysis approaches. Using this cutoff improved the reproducibility of analysis, i.e., variation in richness estimates was reduced by 38% compared with singleton filtering using six human fecal samples across seven sequencing runs. Beta-diversity analysis of human fecal communities was markedly affected by both the filtering strategy and the type of phylogenetic distances used for comparison, highlighting the importance of carefully analyzing data before drawing conclusions on microbiome changes. In summary, handling of artifact sequences during bioinformatic processing of 16S rRNA gene amplicon data requires careful attention to avoid the generation of misleading findings. We propose the concept of effective richness to facilitate the comparison of alpha-diversity across studies.

Weaning-associated feed deprivation stress causes microbiota disruptions in a novel mucin-containing in vitro model of the piglet colon (MPigut-IVM)

BACKGROUND

Risk factors for the etiology of post-weaning diarrhea, a major problem in swine industry associated with enormous economic losses, remain to be fully elucidated. In concordance with the ethical concerns raised by animal experiments, we developed a new in vitro model of the weaning piglet colon (MPigut-IVM) including a mucin bead compartment to reproduce the mucus surface from the gut to which gut microbes can adhere.

RESULTS

Our results indicated that the MPigut-IVM is able to establish a representative piglet archaeal and bacterial colon microbiota in terms of taxonomic composition and function. The MPigut-IVM was consequently used to investigate the potential effects of feed deprivation, a common consequence of weaning in piglets, on the microbiota. The lack of nutrients in the MPigut-IVM led to an increased abundance of Prevotellaceae and Escherichia-Shigella and a decrease in Bacteroidiaceae and confirms previous in vivo findings. On top of a strong increase in redox potential, the feed deprivation stress induced modifications of microbial metabolite production such as a decrease in acetate and an increase in proportional valerate, isovalerate and isobutyrate production.

CONCLUSIONS

The MPigut-IVM is able to simulate luminal and mucosal piglet microbiota and represent an innovative tool for comparative studies to investigate the impact of weaning stressors on piglet microbiota. Besides, weaning-associated feed deprivation in piglets provokes disruptions of MPigut-IVM microbiota composition and functionality and could be implicated in the onset of post-weaning dysbiosis in piglets.

Evolution of fungal community associated with ready-to-eat pineapple during storage under different temperature conditions

The international market of fresh-cut products has witnessed dramatic growth in recent years, stimulated by consumer's demand for healthy, nutritious and convenient foods. One of the main challenging issues for the quality and safety of these products is the potential microbial spoilage that can significantly reduce their shelf-life. The complete identification of fresh-cut product microbiota together with the evaluation of environmental factors impact on microbial composition is of primary importance. We therefore assessed the fungal communities associated with the spoilage of ready-to-eat (RTE) pineapple using a metagenetic amplicon sequencing approach, based on the ITS2 region. Our results revealed a significant variability on fungal species composition between the different batches of RTE pineapple. The initial microbiota composition was the main influencing factor and determined the progress of spoilage. Temperature and storage time were the secondary factors influencing spoilage and their impact was depending on the initial prevalent fungal species, which showed different responses to the various modifications. Our results strongly suggest that further large-scale sampling of RTE pineapple production should be conducted in order to assess the full biodiversity range of fungal community involved in the spoilage process and for unravelling the impact of important environmental factors shaping the initial microbiota.

Microbiome Composition and Borrelia Detection in Ixodes scapularis Ticks at the Northwestern Edge of Their Range

Lyme disease-causing Borrelia burgdorferi has been reported in 10–19% of Ixodes ticks from Alberta, Canada, where the tick vector Ixodes scapularis is at the northwestern edge of its range. However, the presence of Borrelia has not been verified independently, and the bacterial microbiome of these ticks has not been described. We performed 16S rRNA bacterial surveys on female I. scapularis from Alberta that were previously qPCR-tested in a Lyme disease surveillance program. Both 16S and qPCR methods were concordant for the presence of Borrelia. The 16S studies also provided a profile of associated bacteria that showed the microbiome of I. scapularis in Alberta was similar to other areas of North America. Ticks that were qPCR-positive for Borrelia had significantly greater bacterial diversity than Borrelia-negative ticks, on the basis of generalized linear model testing. This study adds value to ongoing tick surveillance and is a foundation for deeper understanding of tick microbial ecology and disease transmission in a region where I. scapularis range expansion, induced by climate and land use changes, is likely to have increasing public health implications.

BIOCOM-PIPE: a new user-friendly metabarcoding pipeline for the characterization of microbial diversity from 16S, 18S and 23S rRNA gene amplicons

Background

The ability to compare samples or studies easily using metabarcoding so as to better interpret microbial ecology results is an upcoming challenge. A growing number of metabarcoding pipelines are available, each with its own benefits and limitations. However, very few have been developed to offer the opportunity to characterize various microbial communities (e.g., archaea, bacteria, fungi, photosynthetic microeukaryotes) with the same tool.

Results

BIOCOM-PIPE is a flexible and independent suite of tools for processing data from high-throughput sequencing technologies, Roche 454 and Illumina platforms, and focused on the diversity of archaeal, bacterial, fungal, and photosynthetic microeukaryote amplicons. Various original methods were implemented in BIOCOM-PIPE to (1) remove chimeras based on read abundance, (2) align sequences with structure-based alignments of RNA homologs using covariance models, and (3) a post-clustering tool (ReClustOR) to improve OTUs consistency based on a reference OTU database. The comparison with two other pipelines (FROGS and mothur) and Amplicon Sequence Variant definition highlighted that BIOCOM-PIPE was better at discriminating land use groups.

Conclusions

The BIOCOM-PIPE pipeline makes it possible to analyze 16S, 18S and 23S rRNA genes in the same packaged tool. The new post-clustering approach defines a biological database from previously analyzed samples and performs post-clustering of reads with this reference database by using open-reference clustering. This makes it easier to compare projects from various sequencing runs, and increased the congruence among results. For all users, the pipeline was developed to allow for adding or modifying the components, the databases and the bioinformatics tools easily, giving high modularity for each analysis.

Bacterial community profile after the lethal infection of Steinernema-Xenorhabdus pairs into soil-reared Tenebrio molitor larvae

The host microbiota may have an impact on pathogens. This is often studied in laboratory-reared hosts but rarely in individuals whose microbiota looks like that of wild animals. In this study, we modified the gut microbiota of the insect Tenebrio molitor by rearing larvae in soil sampled from the field. We showed by high throughput sequencing methods that this treatment modifies the gut microbiota so that it is more diversified than that of laboratory-reared insects, and closely resembled the one of soil-dwelling insects. To describe what the entomopathogenic bacterial symbiont Xenorhabdus (Enterobacteriaceae), vectored by the soil-dwelling nematode Steinernema, might experience in natural conditions, we studied the infestation of the soil-reared T. molitor larvae with three Steinernema-Xenorhabdus pairs. We performed the infestation at 18°C, which delays the emergence of new infective juveniles (IJs), the soil-dwelling nematode forms, but which is a temperature compatible with natural infestation. We analyzed by high throughput sequencing methods the composition of the bacterial community within the insect cadavers before the first emergences of IJs. These bacterial communities were generally characterized by one or two non-symbiont taxa. Even for highly lethal Steinernema-Xenorhabdus pairs, the symbiont does not dominate the bacterial community within the insect cadaver.

Microbial community and geochemical analyses of trans-trench sediments for understanding the roles of hadal environments

Hadal trench bottom (>6000 m below sea level) sediments harbor higher microbial cell abundance compared with adjacent abyssal plain sediments. This is supported by the accumulation of sedimentary organic matter (OM), facilitated by trench topography. However, the distribution of benthic microbes in different trench systems has not been well explored yet. Here, we carried out small subunit ribosomal RNA gene tag sequencing for 92 sediment subsamples of seven abyssal and seven hadal sediment cores collected from three trench regions in the northwest Pacific Ocean: the Japan, Izu-Ogasawara, and Mariana Trenches. Tag-sequencing analyses showed specific distribution patterns of several phyla associated with oxygen and nitrate. The community structure was distinct between abyssal and hadal sediments, following geographic locations and factors represented by sediment depth. Co-occurrence network revealed six potential prokaryotic consortia that covaried across regions. Our results further support that the OM cycle is driven by hadal currents and/or rapid burial shapes microbial community structures at trench bottom sites, in addition to vertical deposition from the surface ocean. Our trans-trench analysis highlights intra- and inter-trench distributions of microbial assemblages and geochemistry in surface seafloor sediments, providing novel insights into ultradeep-sea microbial ecology, one of the last frontiers on our planet.

Bacteria and fungi in acute cholecystitis. A prospective study comparing next generation sequencing to culture

OBJECTIVES

Guidelines for antibiotic treatment of acute cholecystitis are based on studies using culture techniques for microbial identification. Microbial culture has well described limitations and more comprehensive data on the microbial spectrum may support adjustments of these recommendations. We used next generation sequencing to conduct a thorough microbiological characterization of bile-samples from patients with moderate and severe acute cholecystitis.

METHODS

We prospectively included patients with moderate and severe acute cholecystitis, undergoing percutaneous or perioperative drainage of the gall bladder. Bile samples were analyzed using both culture and deep sequencing of bacterial 16S rRNA and rpoB genes and the fungal ITS2-segment. Clinical details were evaluated by medical record review.

RESULTS

Thirty-six patients with moderate and severe acute cholecystitis were included. Bile from 31 (86%) of these contained bacteria (29) and/or fungi (5) as determined by sequencing. Culture identified only 40 (38%) of the 106 microbes identified by sequencing. In none of the 15 polymicrobial samples did culture detect all present microbes. Frequently identified bacteria often missed by culture included oral streptococci, anaerobic bacteria, enterococci and Enterobacteriaceae other than Klebsiella spp. and Escherichia coli.

CONCLUSIONS

Culture techniques display decreased sensitivity for the microbial diagnostics of acute cholecystitis leaving possible pathogens undetected.

Microbiota Composition and Functional Profiling Throughout the Gastrointestinal Tract of Commercial Weaning Piglets

Dietary, environmental, and social stresses induced by weaning transition in pig production are associated with alterations of gut microbiota, diarrhea, and enteric infections. With the boom of -omic technologies, numerous studies have investigated the dynamics of fecal bacterial communities of piglets throughout weaning but much less research has been focused on the composition and functional properties of microbial communities inhabiting other gastrointestinal segments. The objective of the present study was to bring additional information about the piglet bacterial and archaeal microbiota throughout the entire digestive tract, both at the structural level by using quantitative PCR and high-throughput sequencing, and on functionality by measurement of short-chain fatty acids and predictions using Tax4Fun tool. Our results highlighted strong structural and functional differences between microbial communities inhabiting the fore and the lower gut as well as a quantitatively important archaeal community in the hindgut. The presence of opportunistic pathogens was also noticed throughout the entire digestive tract and could trigger infection emergence. Understanding the role of the intestinal piglet microbiota at weaning could provide further information about the etiology of post-weaning infections and lead to the development of effective preventive solutions.

Improving eDNA-based protist diversity assessments using networks of amplicon sequence variants

Effective and precise grouping of highly similar sequences remains a major bottleneck in the evaluation of high-throughput sequencing datasets. Amplicon sequence variants (ASVs) offer a promising alternative that may supersede the widely used operational taxonomic units (OTUs) in environmental sequencing studies. We compared the performance of a recently developed pipeline based on the algorithm DADA2 for obtaining ASVs against a pipeline based on the algorithm SWARM for obtaining OTUs. Illumina-sequencing of 29 individual ciliate species resulted in up to 11 ASVs per species, while SWARM produced up to 19 OTUs per species. To improve the congruency between species diversity and molecular diversity, we applied sequence similarity networks (SSNs) for second-level sequence grouping into network sequence clusters (NSCs). At 100% sequence similarity in SWARM-SSNs, NSC numbers decreased from 7.9-fold overestimation without abundance filter, to 4.5-fold overestimation when an abundance filter was applied. For the DADA2-SSN approach, NSC numbers decreased from 3.5-fold to 3-fold overestimation. Rand index cluster analyses predicted best binning results between 97% and 94% sequence similarity for both DADA2-SSNs and SWARM-SSNs. Depending on the ecological questions addressed in an environmental sequencing study with protists we recommend ASVs as replacement for OTUs, best in combination with SSNs.

Rumen bacterial community responses to DPA, EPA and DHA in cattle and sheep: A comparative in vitro study

The role of marine lipids as modulators of ruminal biohydrogenation of dietary unsaturated fatty acids may be explained by the effects of their n-3 polyunsaturated fatty acids (PUFA) on the bacterial community. However, the impact of individual PUFA has barely been examined, and it is uncertain which bacteria are truly involved in biohydrogenation. In addition, despite interspecies differences in rumen bacterial composition, we are not aware of any direct comparison of bovine and ovine responses to dietary PUFA. Therefore, rumen fluid from cannulated cattle and sheep were used as inocula to examine in vitro the effect of 20:5n-3 (EPA), 22:5n-3 (DPA), and 22:6n-3 (DHA) on the bacterial community. Amplicon 16 S rRNA sequencing suggested that EPA and DHA had a greater contribution to the action of marine lipids than DPA both in cattle and sheep. Certain effects were exclusive to each ruminant species, which underlines the complexity of rumen microbial responses to dietary fatty acids. Based on changes in bacterial abundance, Barnesiella, Prevotella, Paraprevotella, Hallela, Anaerovorax, Succiniclasticum, Ruminococcus and Ruminobacter may be involved in the ruminal response in biohydrogenation to the addition of marine lipids, but further research is necessary to confirm their actual role in ruminal lipid metabolism.

Temporal Stability of the Ruminal Bacterial Communities in Beef Steers

Nutritional studies involving ruminants have traditionally relied on relatively short transition or wash-out periods between dietary treatments, typically two to four weeks. However, little is known about adequate adaptation periods required to reach stabilization of the rumen microbiome that could provide more accurate results from nutritional studies in ruminants. This study determined the rumen bacterial communities and rumen environment parameters over ten weeks following transition from a forage-based to concentrate-based diet. Several α-diversity metrics, including observed OTUs and Simpson’s Evenness fluctuated throughout the trial, but were typically either greatest (observed OTUs) or lowest (Simpson’s) at week 5 of the trial contrasted from weeks 1 and 10 (P < 0.05). At week 4, several orders associated with the shift to the final bacterial community composition, including Pasteurellales, Aeromonadales, and Bacteroidales. At week 5, rumen pH was correlated with α-diversity (P = 0.005) and predictive of the rumen microbiome signature at week 10 (R² = 0.48; P = 0.04). Rumen microbiome stability did not occur until approximately 9 weeks following adaptation to the diet and was associated with changes in specific bacterial populations and rumen environment. The results of this study suggest that adaptation and wash-out periods must be re-evaluated in order to accommodate necessary rumen microbiome acclimation.

Composition and dynamics of the bacterial communities present in the post-slaughter environment of farmed Atlantic salmon (Salmo salar L.) and correlations to gelatin degrading activity

Background

Microbial analyses performed in connection with the post-slaughter environment of farmed Atlantic salmon (Salmo salar L.) have mostly focused on specific bacteria that may have negative effects on the health of consumers. However, bacteria may also affect other quality variables. The objective of this study was to provide general knowledge about composition and dynamics of the bacterial communities present at slaughter and cold storage of farmed Atlantic salmon, as well as reveal any possible correlations to gelatinase activity, which may affect fillet quality. Thus, these data may provide a basis for optimization opportunities in the aquaculture industry.

Methods

Samples were taken from the digestive system harvested from 15 salmon immediately after slaughter. Another 17 salmon were taken from the processing line just before the final cleaning stage; of these eight were distributed in three iced storage boxes while the other nine were rinsed an extra time with industrial water before being distributed into another three storage boxes. In the following 6 days, samples were taken of skin mucus, liquids in the abdominal cavity and the storage ice. The compositions of the bacterial communities were analyzed by next-generation sequencing and gelatinase activity was measured in all samples except the storage ice.

Results

The bacterial communities in the digestive tract samples were dominated by the family Mycoplasmataceae. The genus Aliivibrio was also relatively abundant. Bacterial communities in the abdominal cavity were generally more diverse than the intestinal samples. However, all of the abdominal samples from storage box no. 3 had a high relative abundance of Mycoplasmataceae, and could not be distinguished from the intestinal samples (Q = 1.27, p = 0.633) while being significantly different from the other abdominal samples (Q = 9.02, p = 0.01). In addition, the abdominal samples from storage box no. 3 had a significantly higher gelatin degrading activity (Q = 9.43, p = 0.001) than those from the other storage boxes and similar to the high gelatinase activity in the intestinal samples. This indicated that in storage box no. 3 there was a transfer of intestinal fluids to the abdominal cavities, which was not removed by the cleaning procedure. There was a significant difference of the major phyla detected in the skin mucus of salmon rinsed an additional time, as these salmon had a higher relative amount of Firmicutes (F = 4.76, p = 0.04) and lower amount of Proteobacteria (F = 4.41, p = 0.047).

Conclusions

The study showed a correlation between intestinal fluids and bacteria left in the abdominal cavity and gelatinase activity. This suggested that intestinal fluids and/or bacteria could enhance the degradation of connective tissue in the abdominal cavity and hence negatively affect the fillet quality. In addition, the study provided general knowledge of the composition and dynamics of bacterial communities present.

A longitudinal characterization of the Non-Cystic Fibrosis Bronchiectasis airway microbiome

A diverse microbiota exists within the airways of individuals with non-cystic fibrosis bronchiectasis (nCFB). How the lung microbiome evolves over time, and whether changes within the microbiome correlate with future disease progression is not yet known. We assessed the microbial community structure of 133 serial sputa and subsequent disease course of 29 nCFB patients collected over a span of 4–16 years using 16S rRNA paired-end sequencing. Interestingly, no significant shifts in the microbial community of individuals were observed during extended follow-up suggesting the microbiome remains relatively stable over prolonged periods. Samples that were Pseudomonas aeruginosa culture positive displayed markedly different microbial community structures compared to those that were positive for Haemophilus influenzae. Importantly, patients with sputum of lower microbial community diversity were more likely to experience subsequent lung function decline as defined by annual change in ≥−1 FEV₁% predicted. Shannon diversity values <1 were more prevalent in patients with FEV₁ decline (P = 0.002). However, the relative abundance of particular core microbiota constituents did not associate with risk of decline. Here we present data confirming that the microbiome of nCFB individuals is generally stable, and that microbiome-based measurements may have a prognostic role as biomarkers for nCFB.

Reduction in Methane Emissions From Acidified Dairy Slurry Is Related to Inhibition of Methanosarcina Species

Liquid dairy manure treated with sulfuric acid was stored in duplicate pilot-scale storage tanks for 120 days with continuous monitoring of CH₄ emissions and concurrent examination of changes in the structure of bacterial and methanogenic communities. Methane emissions were monitored at the site using laser-based Trace Gas Analyzer whereas quantitative real-time polymerase chain reaction and massively parallel sequencing were employed to study bacterial and methanogenic communities using 16S rRNA and methyl-coenzyme M Reductase A (mcrA) genes/transcripts, respectively. When compared with untreated slurries, acidification resulted in 69–84% reductions of cumulative CH₄ emissions. The abundance, activity, and proportion of bacterial communities did not vary with manure acidification. However, the abundance and activity of methanogens (as estimated from mcrA gene and transcript copies, respectively) in acidified slurries were reduced by 6 and 20%, respectively. Up to 21% reduction in mcrA transcript/gene ratios were also detected in acidified slurries. Regardless of treatment, Methanocorpusculum predominated archaeal 16S rRNA and mcrA gene and transcript libraries. The proportion of Methanosarcina, which is the most metabolically-diverse methanogen, was the significant discriminant feature between acidified and untreated slurries. In acidified slurries, the relative proportions of Methanosarcina were ≤ 10%, whereas in untreated slurries, it represented up to 24 and 53% of the mcrA gene and transcript libraries, respectively. The low proportions of Methanosarcina in acidified slurries coincided with the reductions in CH₄ emissions. The results suggest that reduction of CH₄ missions achieved by acidification was due to an inhibition of the growth and activity of Methanosarcina species.

Detection of an amplification bias associated to Leuconostocaceae family with a universal primer routinely used for monitoring microbial community structures within food products

Objectives

Sequencing of 16S rDNA V3–V4 region is widely applied for food community profiling. However, two different universal forward primers (named here MUYZER-primer1 and KLINDWORTH-primer2) targeting an identical conservative sequence upstream of the V3 region of 16S rRNA gene, and only distinguished by a single mismatch are both used. This study was carried out to compare whether the accuracy of food microbiota analysis would depend on the choice of one of these two primers.

Results

Alignment of both primers with common food-borne bacteria 16S sequences revealed that the mismatch between both primers might specifically affect the amplification of Leuconostoc, Oenococcus and Fructobacillus species but not Weissella species. Food products containing either Leuconostoc and/or Weissella were selected for a detection test. As expected from our in silico analysis, our study showed that this mismatch induced a strong biased amplification specifically associated to the OTUs belonging to the genus Leuconostoc but not to the genus Weissella. In presence of Muyzer-primer1, none of the sequences expected for Leuconostoc genus was detected whereas those sequences were correctly amplified with Klindworth-primer2. Since Leuconostoc is an important genus in food, agro-environments and in digestive tract of animals, we recommend that Muyzer-primer1 should thus be abandoned for the bacterial characterization of their associated microbiota.

Anaerobic lignocellulolytic microbial consortium derived from termite gut: enrichment, lignocellulose degradation and community dynamics

BACKGROUND

Lignocellulose is the most abundant renewable carbon resource that can be used for biofuels and commodity chemicals production. The ability of complex microbial communities present in natural environments that are specialized in biomass deconstruction can be exploited to develop lignocellulose bioconversion processes. Termites are among the most abundant insects on earth and play an important role in lignocellulose decomposition. Although their digestive microbiome is recognized as a potential reservoir of microorganisms producing lignocellulolytic enzymes, the potential to enrich and maintain the lignocellulolytic activity of microbial consortia derived from termite gut useful for lignocellulose biorefinery has not been assessed. Here, we assessed the possibility of enriching a microbial consortium from termite gut and maintaining its lignocellulose degradation ability in controlled anaerobic bioreactors.

RESULTS

We enriched a termite gut-derived consortium able to transform lignocellulose into carboxylates under anaerobic conditions. To assess the impact of substrate natural microbiome on the enrichment and the maintenance of termite gut microbiome, the enrichment process was performed using both sterilized and non-sterilized straw. The enrichment process was carried out in bioreactors operating under industrially relevant aseptic conditions. Two termite gut-derived microbial consortia were obtained from Nasutitermes ephratae by sequential batch culture on raw wheat straw as the sole carbon source. Analysis of substrate loss, carboxylate production and microbial diversity showed that regardless of the substrate sterility, the diversity of communities selected by the enrichment process strongly changed compared to that observed in the termite gut. Nevertheless, the community obtained on sterile straw displayed higher lignocellulose degradation capacity; it showed a high xylanase activity and an initial preference for hemicellulose.

CONCLUSIONS

This study demonstrates that it is possible to enrich and maintain a microbial consortium derived from termite gut microbiome in controlled anaerobic bioreactors, producing useful carboxylates from raw biomass. Our results suggest that the microbial community is shaped both by the substrate and the conditions that prevail during enrichment. However, when aseptic conditions are applied, it is also affected by the biotic pressure exerted by microorganisms naturally present in the substrate and in the surrounding environment. Besides the efficient lignocellulolytic consortium enriched in this study, our results revealed high levels of xylanase activity that can now be further explored for enzyme identification and overexpression for biorefinery purposes.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Uncovering the Potential of Termite Gut Microbiome for Lignocellulose Bioconversion in Anaerobic Batch Bioreactors

Termites are xylophages, being able to digest a wide variety of lignocellulosic biomass including wood with high lignin content. This ability to feed on recalcitrant plant material is the result of complex symbiotic relationships, which involve termite-specific gut microbiomes. Therefore, these represent a potential source of microorganisms for the bioconversion of lignocellulose in bioprocesses targeting the production of carboxylates. In this study, gut microbiomes of four termite species were studied for their capacity to degrade wheat straw and produce carboxylates in controlled bioreactors. All of the gut microbiomes successfully degraded lignocellulose and up to 45% w/w of wheat straw degradation was observed, with the Nasutitermes ephratae gut-microbiome displaying the highest levels of wheat straw degradation, carboxylate production and enzymatic activity. Comparing the 16S rRNA gene diversity of the initial gut inocula to the bacterial communities in lignocellulose degradation bioreactors revealed important changes in community diversity. In particular, taxa such as Spirochaetes and Fibrobacteres that were highly abundant in the initial gut inocula were replaced by Firmicutes and Proteobacteria at the end of incubation in wheat straw bioreactors. Overall, this study demonstrates that termite-gut microbiomes constitute a reservoir of lignocellulose-degrading bacteria that can be harnessed in artificial conditions for biomass conversion processes that lead to the production of useful molecules.