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

BACKGROUND: With the emergence of Oxford Nanopore technology, now the on-site sequencing of 16S rRNA from environments is available. Due to the error level and structure, the analysis of such data demands some database of reference sequences. However, many taxa from complex and diverse environments, have poor representation in publicly available databases. In this paper, we propose the METASEED pipeline for the reconstruction of full-length 16S sequences from such environments, in order to improve the reference for the subsequent use of on-site sequencing. RESULTS: We show that combining high-precision short-read sequencing of both 16S and full metagenome from the same samples allow us to reconstruct high-quality 16S sequences from the more abundant taxa. A significant novelty is the carefully designed collection of metagenome reads that matches the 16S amplicons, based on a combination of uniqueness and abundance. Compared to alternative approaches this produces superior results. CONCLUSION: Our pipeline will facilitate numerous studies associated with various unknown microorganisms, thus allowing the comprehension of the diverse environments. The pipeline is a potential tool in generating a full length 16S rRNA gene database for any environment.

Spatial fractionation of phosphorus accumulating biofilm: stratification of polyphosphate accumulation and dissimilatory nitrogen metabolism.

The spatial distribution of microorganisms represents a critical issue in understanding biofilm function. The aim of the current work was to develop a method for biofilm fractionation, facilitating the analysis of individual spatial biofilm layers using metagenomic approaches. Phosphorus accumulating biofilm applied in an enhanced biological phosphorus removal wastewater treatment plant, were fractionated, and analyzed. The fractionated biofilm revealed a surprising spatial distribution of bacteria and genes, where potential polyphosphate accumulating organisms account for ∼ 47% of the inner layer microbiome. A spatial distribution of genes involved in dissimilatory nitrogen reduction was observed, indicating that different layers of the biofilm were metabolically active during the anoxic and aerobic zones of the phosphorus removal process. The physical biofilm separation into individual fractions unveiled functional layers of the biofilm, which will be important for future understanding of the phosphorus removal wastewater process.

Fecal Microbiota Nutrient Utilization Potential Suggests Mucins as Drivers for Initial Gut Colonization of Mother-Child-Shared Bacteria.

The nutritional drivers for mother-child sharing of bacteria and the corresponding longitudinal trajectory of the infant gut microbiota development are not yet completely settled. We therefore aimed to characterize the mother-child sharing and the inferred nutritional utilization potential for the gut microbiota from a large unselected cohort. We analyzed in depth gut microbiota in 100 mother-child pairs enrolled antenatally from the general population-based Preventing Atopic Dermatitis and Allergies in Children (PreventADALL) cohort. Fecal samples collected at gestational week 18 for mothers and at birth (meconium), 3, 6, and 12 months for infants were analyzed by reduced metagenome sequencing to determine metagenome size and taxonomic composition. The nutrient utilization potential was determined based on the Virtual Metabolic Human (VMH, www.vmh.life) database. The estimated median metagenome size was ∼150 million base pairs (bp) for mothers and ∼20 million bp at birth for the children. Longitudinal analyses revealed mother-child sharing (P < 0.05, chi-square test) from birth up to 6 months for 3 prevalent Bacteroides species (prevalence, >25% for all age groups). In a multivariate analysis of variance (ANOVA), the mother-child-shared Bacteroides were associated with vaginal delivery (1.7% explained variance, P = 0.0001). Both vaginal delivery and mother-child sharing were associated with host-derived mucins as nutrient sources. The age-related increase in metagenome size corresponded to an increased diversity in nutrient utilization, with dietary polysaccharides as the main age-related factor. Our results support host-derived mucins as potential selection means for mother-child sharing of initial colonizers, while the age-related increase in diversity was associated with dietary polysaccharides.IMPORTANCE The initial bacterial colonization of human infants is crucial for lifelong health. Understanding the factors driving this colonization will therefore be of great importance. Here, we used a novel high-taxonomic-resolution approach to deduce the nutrient utilization potential of the infant gut microbiota in a large longitudinal mother-child cohort. We found mucins as potential selection means for the initial colonization of mother-child-shared bacteria, while the transition to a more adult-like microbiota was associated with dietary polysaccharide utilization potential. This knowledge will be important for a future understanding of the importance of diet in shaping the gut microbiota composition and development during infancy.

A game theory model for gut bacterial nutrient utilization strategies during human infancy.

Despite the fact that infant gut colonization patterns have been extensively studied, we have limited knowledge about the underlying ecological processes. This particularly relates to the ecological choice of nutrient utilization strategies. The aim of the current study was therefore to compare empirically determined nutrient utilization strategies with that expected from a combinatorial game theory model. Observational analyses for 100 mother-child pairs suggested mother-child transmission of specialists with the potential to use few nutrients. Generalists, on the other hand, with the potential to use many nutrients, peaked at three months of age for the children. The level of generalists was gradually replaced with specialists up to 12 months of age. Game theory simulation revealed a competitive advantage of generalists in an expanding population, while more specialized bacteria were favoured with the maturation of the population. This suggests that the observed increase in generalists in the three-month-old children could be due to an immature, expanding gut microbiota population while the increase of specialists at 12 months could be due to population maturation. The simulated and empirical data also correspond with respect to an increased α diversity and a decreased ß diversity with the number of simulations and age, respectively. Taken together, game theory simulation of nutrient utilization strategies can therefore provide novel insight into the maturation of the human gut microbiota during infancy.

Experimental support for multidrug resistance transfer potential in the preterm infant gut microbiota.

BACKGROUND: There is currently a lack of experimental evidence for horizontal gene transfer (HGT) mechanisms in the human gut microbiota. The aim of this study was therefore to experimentally determine the HGT potential in the microbiota of a healthy preterm infant twin pair and to evaluate the global occurrence of the mobilized elements. METHODS: Stool samples were collected. Both shotgun metagenome sequencing and bacterial culturing were done for the same samples. A range of experimental conditions were used to test DNA transfer for the cultured isolates. Searches for global distribution of transferable elements were done for the ~120,000 metagenomic samples in the Sequence Read Archive (SRA) database. RESULTS: DNA transfer experiments demonstrated frequent transmission of an ESBL encoding IncI1 plasmid, a high copy number ColEI plasmid, and bacteriophage P1. Both IncI1 and ColE1 were abundant in the stool samples. In vitro competition experiments showed that transconjugants containing IncI1 plasmids outcompeted the recipient strain in the absence of antibiotic selection. The SRA searches indicated a global distribution of the mobilizable elements, with chicken identified as a possible reservoir for the IncI1 ESBL encoding plasmid. CONCLUSION: Our results experimentally support a major horizontal transmission and persistence potential of the preterm infant gut microbiota mobilome involving genes encoding ESBL.

Microbial ecological processes in MBBR biofilms for biological phosphorus removal from wastewater.

Phosphorus is both a major environmental pollutant and a limiting resource. Although enhanced biological phosphorus removal (EBPR) is used worldwide for phosphorus removal, the standard activated sludge-based EBPR process shows limitations with stability and efficiency. Recently, a new EBPR moving bed biofilm reactor (MBBR) process has been developed at HIAS (Hamar, Norway), enabling a phosphorus removal stability above 90% during a whole year cycle. To increase the knowledge of the HIAS (MBBR) process the aim of the current work was to characterize the MBBR microbiota and operational performance weekly for the operational year. Surprisingly, we found a major succession of the microbiota, with a five-fold increase in phosphorus accumulating organisms (PAOs), and major shifts in eukaryote composition, despite a stable phosphorus removal. Temperature was the only factor that significantly affected both phosphorus removal and the microbiota. There was a lower phosphor removal during the winter, coinciding with a higher microbiota alpha diversity, and a lower beta diversity. This differs from what is observed for activated sludge based EBPR. Taken together, the knowledge gained from the current microbiota study supports the efficiency and stability of MBBR-based systems, and that knowledge from activated sludge-based EBPR approaches cannot be translated to MBBR systems.

Stable Core Gut Microbiota across the Freshwater-to-Saltwater Transition for Farmed Atlantic Salmon.

Gut microbiota associations through habitat transitions are fundamentally important yet poorly understood. One such habitat transition is the migration from freshwater to saltwater for anadromous fish, such as salmon. The aim of the current work was therefore to determine the freshwater-to-saltwater transition impact on the gut microbiota in farmed Atlantic salmon, with dietary interventions resembling freshwater and saltwater diets with respect to fatty acid composition. Using deep 16S rRNA gene sequencing and quantitative PCR, we found that the freshwater-to-saltwater transition had a major association with the microbiota composition and quantity, while diet did not show significant associations with the microbiota. In saltwater there was a 100-fold increase in bacterial quantity, with a relative increase of Firmicutes and a relative decrease of both Actinobacteria and Proteobacteria Irrespective of an overall shift in microbiota composition from freshwater to saltwater, we identified three core clostridia and one Lactobacillus-affiliated phylotype with wide geographic distribution that were highly prevalent and co-occurring. Taken together, our results support the importance of the dominating bacteria in the salmon gut, with the freshwater microbiota being immature. Due to the low number of potentially host-associated bacterial species in the salmon gut, we believe that farmed salmon can represent an important model for future understanding of host-bacterium interactions in aquatic environments.IMPORTANCE Little is known about factors affecting the interindividual distribution of gut bacteria in aquatic environments. We have shown that there is a core of four highly prevalent and co-occurring bacteria irrespective of feed and freshwater-to-saltwater transition. The potential host interactions of the core bacteria, however, need to be elucidated further.

Allergen Immunization Induces Major Changes in Microbiota Composition and Short-Chain Fatty Acid Production in Different Gut Segments in a Mouse Model of Lupine Food Allergy.

BACKGROUND: The incidence of food allergies in western countries has increased in recent decades. OBJECTIVES: To study the association between gut bacterial microbiota composition, short-chain fatty acids (SCFAs) and food allergy in a mouse model. METHODS: After oral immunizations with the human food allergen lupine with the adjuvant cholera toxin (CT) (or buffer in controls), sensitization and anaphylactic responses were determined. Gastrointestinal content was collected from the distal ileum, cecum, colon, and fecal pellets, and the bacterial diversity and composition was determined by deep sequencing of the 16S rRNA gene. SCFAs in gastrointestinal content supernatants were determined by gas chromatography. RESULTS: The microbiota signatures were profoundly affected by allergen immunization. Ten operational taxonomic units (OTUs) were significantly different between immunized and control animals for at least one of the intestinal segments; eight of these OTUs belonged to the Clostridia class. Although consistent across all four gut segments, the colon showed the highest number of OTUs significantly associated with allergic immunization. SCFA levels in the cecum were also altered by immunization. CONCLUSIONS: Allergen immunization with CT in the present food allergy model induced profound changes in the microbiome composition and SCFA production. The result suggests that the colon may be the most sensitive gut segment for investigating changes in the gut microbiome.

Is amniotic fluid of women with uncomplicated term pregnancies free of bacteria?

BACKGROUND: The "sterile womb" paradigm is debated. Recent evidence suggests that the offspring's first microbial encounter is before birth in term uncomplicated pregnancies. The establishment of a healthy microbiota early in life might be crucial for reducing the burden of diseases later in life. OBJECTIVE: We aimed to investigate the presence of a microbiota in sterilely collected amniotic fluid in uncomplicated pregnancies at term in the Preventing Atopic Dermatitis and Allergies in children (PreventADALL) study cohort. STUDY DESIGN: Amniotic fluid was randomly sampled at cesarean deliveries in pregnant women in 1 out of 3 study sites included in the PreventADALL study. From 65 pregnancies at term, where amniotic fluid was successfully sampled, we selected 10 from elective (planned, without ongoing labor) cesarean deliveries with intact amniotic membranes and all 14 with prior rupture of membranes were included as positive controls. Amniotic fluid was analyzed by culture-independent and culture-dependent techniques. RESULTS: The median (min-max) concentration of prokaryotic DNA (16S rRNA gene copies/mL; digital droplet polymerase chain reaction) was low for the group with intact membranes [664 (544-748)]-corresponding to the negative controls [596 (461-679)], while the rupture of amniotic membranes group had >10-fold higher levels [7700 (1066-251,430)] (P = .0001, by Mann-Whitney U test). Furthermore, bacteria were detected in 50% of the rupture of amniotic membranes samples by anaerobic culturing, while none of the intact membranes samples showed bacterial growth. Sanger sequencing of the rupture of amniotic membrane samples identified bacterial strains that are commonly part of the vaginal flora and/or associated with intrauterine infections. CONCLUSION: We conclude that fetal development in uncomplicated pregnancies occurs in the absence of an amniotic fluid microbiota and that the offspring microbial colonization starts after uterine contractions and rupture of amniotic membrane.

Swarm and UNOISE outperform DADA2 and Deblur for denoising high-diversity marine seafloor samples.

The performance of sequence variant resolution analytic tools for metabarcoding has not yet been adequately benchmarked for high-diversity environmental samples. We therefore evaluated the sequence variant tools DADA2, Deblur, Swarm, and UNOISE, using high-diversity seafloor samples, resulting in comparisons of 1800 sequence variant tables. The evaluation was based on 30 sediment grab samples, for which 3 replica samples were collected. Each replica sample was extracted using 5 common DNA extraction kits, resulting in 450 DNA extracts which were 16S rRNA gene sequenced (V3-V4), using Illumina. Assessments included variation across replica samples, extraction kits, and denoising methods, in addition to applying prior knowledge about alpha diversity correlations toward the cosmopolitan marine archaeon Nitrosopumilus with high diversity and the sulfide oxidizing Sulfurovum with low diversity. DADA2 displayed the highest variance between replicates (Manhattan distance 1.14), while Swarm showed the lowest variance (Manhattan distance 0.93). For the analysis based on prior biological knowledge, UNOISE displayed the highest alpha diversity (Simpson's D) correlation toward Nitrosopumilus (Spearman rho = 0.85), while DADA2 showed the lowest (Spearman rho = 0.10). Deblur completely eliminated Nitrosopumilus from the dataset. For Sulfurovum, on the other hand, all the methods showed comparable results. In conclusion, our evaluations show that Swarm and UNOISE performed better than DADA2 and Deblur for high-diversity seafloor samples.

Reduced metagenome sequencing for strain-resolution taxonomic profiles.

BACKGROUND: Studies of shifts in microbial community composition has many applications. For studies at species or subspecies levels, the 16S amplicon sequencing lacks resolution and is often replaced by full shotgun sequencing. Due to higher costs, this restricts the number of samples sequenced. As an alternative to a full shotgun sequencing we have investigated the use of Reduced Metagenome Sequencing (RMS) to estimate the composition of a microbial community. This involves the use of double-digested restriction-associated DNA sequencing, which means only a smaller fraction of the genomes are sequenced. The read sets obtained by this approach have properties different from both amplicon and shotgun data, and analysis pipelines for both can either not be used at all or not explore the full potential of RMS data. RESULTS: We suggest a procedure for analyzing such data, based on fragment clustering and the use of a constrained ordinary least square de-convolution for estimating the relative abundance of all community members. Mock community datasets show the potential to clearly separate strains even when the 16S is 100% identical, and genome-wide differences is < 0.02, indicating RMS has a very high resolution. From a simulation study, we compare RMS to shotgun sequencing and show that we get improved abundance estimates when the community has many very closely related genomes. From a real dataset of infant guts, we show that RMS is capable of detecting a strain diversity gradient for Escherichia coli across time. CONCLUSION: We find that RMS is a good alternative to either metabarcoding or shotgun sequencing when it comes to resolving microbial communities at the strain level. Like shotgun metagenomics, it requires a good database of reference genomes and is well suited for studies of the human gut or other communities where many reference genomes exist. A data analysis pipeline is offered, as an R package at https://github.com/larssnip/microRMS . Video abstract.

Naturalizing laboratory mice by housing in a farmyard-type habitat confers protection against colorectal carcinogenesis.

Living in a farm environment in proximity to animals is associated with reduced risk of developing allergies and asthma, and has been suggested to protect against other diseases, such as inflammatory bowel disease and cancer. Despite epidemiological evidence, experimental disease models that recapitulate such environments are needed to understand the underlying mechanisms. In this study, we show that feralizing conventional inbred mice by continuous exposure to a livestock farmyard-type environment conferred protection toward colorectal carcinogenesis. Two independent experimental approaches for colorectal cancer induction were used; spontaneous (Apc Min/+ mice on an A/J background) or chemical (AOM/DSS). In contrast to conventionally reared laboratory mice, the feralized mouse gut microbiota structure remained stable and resistant to mutagen- and colitis-induced neoplasia. Moreover, the feralized mice exhibited signs of a more mature immunophenotype, indicated by increased expression of NK and T-cell maturation markers, and a more potent IFN-γ response to stimuli. In our study, hygienically born and raised mice subsequently feralized post-weaning were protected to a similar level as life-long exposed mice, although the greatest effect was seen upon neonatal exposure. Collectively, we show protective implications of a farmyard-type environment on colorectal cancer development and demonstrate the utility of a novel animal modeling approach that recapitulates realistic disease responses in a naturalized mammal.

Ecological Processes Affecting Long-Term Eukaryote and Prokaryote Biofilm Persistence in Nitrogen Removal from Sewage.

The factors affecting long-term biofilm stability in sewage treatment remain largely unexplored. We therefore analyzed moving bed bioreactors (MBBRs) biofilm composition and function two years apart from four reactors in a nitrogen-removal sewage treatment plant. Multivariate ANOVA revealed a similar prokaryote microbiota composition on biofilm carriers from the same reactors, where reactor explained 84.6% of the variance, and year only explained 1.5%. Eukaryotes showed a less similar composition with reactor explaining 56.8% of the variance and year 9.4%. Downstream effects were also more pronounced for eukaryotes than prokaryotes. For prokaryotes, carbon source emerged as a potential factor for deterministic assembly. In the two reactors with methanol as a carbon source, the bacterial genus Methylotenera dominated, with M. versatilis as the most abundant species. M. versatilis showed large lineage diversity. The lineages mainly differed with respect to potential terminal electron acceptor usage (nitrogen oxides and oxygen). Searches in the Sequence Read Archive (SRA) database indicate a global distribution of the M. versatilis strains, with methane-containing sediments as the main habitat. Taken together, our results support long-term prokaryote biofilm persistence, while eukaryotes were less persistent.

Pre- and Post-Race Intestinal Microbiota in Long-Distance Sled Dogs and Associations with Performance.

Although our understanding of the role of the gut microbiota in different diseases is improving, our knowledge regarding how the gut microbiota affects functioning in healthy individuals is still limited. Here, we hypothesize that the gut microbiota could be associated with sled dog endurance-race performance. We investigated the gut microbiota in 166 fecal samples from 96 Alaskan Huskies, representing 16 teams participating in the 2016 Femund Race (400 km) in Norway, relating the microbiota composition to performance and metadata derived from questionnaires. For 16S rRNA gene sequencing-derived compositional data, we found a strong negative association between Enterobacteriaceae (dysbiosis-associated) and Clostridium hiranonis (normobiosis-associated). The teams with the best performances showed both the lowest levels of dysbiosis-associated bacteria prior to the race and the lowest change (decrease) in these bacteria after the race. Taken together, our results support the hypothesis that normobiosis-associated bacteria are involved in resilience mechanisms, potentially preventing growth of Enterobacteriaceae during the race.

De novo species identification using 16S rRNA gene nanopore sequencing.

Nanopore sequencing is rapidly becoming more popular for use in various microbiota-based applications. Major limitations of current approaches are that they do not enable de novo species identification and that they cannot be used to verify species assignments. This severely limits applicability of the nanopore sequencing technology in taxonomic applications. Here, we demonstrate the possibility of de novo species identification and verification using hexamer frequencies in combination with k-means clustering for nanopore sequencing data. The approach was tested on the human infant gut microbiota of 3-month-old infants. Using the hexamer k-means approach we identified two new low abundant species associated with vaginal delivery. In addition, we confirmed both the vaginal delivery association for two previously identified species and the overall high levels of bifidobacteria. Taxonomic assignments were further verified by mock community analyses. Therefore, we believe our de novo species identification approach will have widespread application in analyzing microbial communities in the future.

Association between diet and rumen microbiota in wild roe deer.

The association between diet and the rumen microbiota for wild animals remains largely unexplored. Here, we explored this association using a combination of 16S rRNA gene sequencing to determine the prokaryote microbiota and 18S rRNA gene sequencing to determine the dietary components for wild roe deer. These analyses revealed a wide diversity of dietary components, with over-representation of Bacteroidetes for the diet-correlating bacteria. Ruminococcus, on the other hand, dominated the stable diet-independent part of the microbiota. Taken together, the combination of 16S and 18S rRNA gene analyses provide novel insight into rumen microbiota ecology.

Low Maternal Microbiota Sharing across Gut, Breast Milk and Vagina, as Revealed by 16S rRNA Gene and Reduced Metagenomic Sequencing.

The maternal microbiota plays an important role in infant gut colonization. In this work we have investigated which bacterial species are shared across the breast milk, vaginal and stool microbiotas of 109 women shortly before and after giving birth using 16S rRNA gene sequencing and a novel reduced metagenomic sequencing (RMS) approach in a subgroup of 16 women. All the species predicted by the 16S rRNA gene sequencing were also detected by RMS analysis and there was good correspondence between their relative abundances estimated by both approaches. Both approaches also demonstrate a low level of maternal microbiota sharing across the population and RMS analysis identified only two species common to most women and in all sample types (Bifidobacterium longum and Enterococcus faecalis). Breast milk was the only sample type that had significantly higher intra- than inter- individual similarity towards both vaginal and stool samples. We also searched our RMS dataset against an in silico generated reference database derived from bacterial isolates in the Human Microbiome Project. The use of this reference-based search enabled further separation of Bifidobacterium longum into Bifidobacterium longum ssp. longum and Bifidobacterium longum ssp. infantis. We also detected the Lactobacillus rhamnosus GG strain, which was used as a probiotic supplement by some women, demonstrating the potential of RMS approach for deeper taxonomic delineation and estimation.

Comparison of reduced metagenome and 16S rRNA gene sequencing for determination of genetic diversity and mother-child overlap of the gut associated microbiota.

Use of the 16S rRNA gene in microbiota studies is limited by the lack of taxonomic and functional resolution. High resolution analyses are particularly important for understanding transmission and persistence of bacteria. The aim of our work was therefore to compare a novel reduced metagenome sequencing (RMS) approach with 16S rRNA gene sequencing to determine both the metagenome genetic diversity and the mother-to-child sharing of the microbiota in a cohort of 17 mother-child pairs. We found that although both approaches gave comparable results with respect to sample separation and taxonomy, RMS gave higher resolution and the potential for genomic-/functional assignment. Using RMS we estimated that the metagenome size increased from about 60 Mbp for 4-day-old children to about 225 Mbp for mothers. The 4-day-old children shared 7% of the metagenome sequences with the mothers, while the metagenome sequence sharing was >30% among the mothers. We found 15 genomes shared across >50% of the mothers, of which 10 belonged to Clostridia. Only Bacteroides showed a direct mother-child association, with B. vulgatus being abundant in both 4-day-old children and mothers. For the functional assignments, we identified a significant association between antibiotic usage during labor, and quantity of Fosfomycin resistance genes. In conclusion, our results show a higher functional and taxonomic resolution for RMS compared to 16S rRNA gene sequencing, where RMS enabled a detailed description of mother to child gut microbiota transmission - supporting a late recruitment of most gut bacteria and an effect of antibiotic treatment during labor on infant antibiotic resistance gene patterns.