The rate of caesarean section delivery (CSD) is increasing worldwide. It remains unclear whether disruption of mother-to-neonate transmission of microbiota through CSD occurs and whether it affects human physiology. Here we perform metagenomic analysis of earliest gut microbial community structures and functions. We identify differences in encoded functions between microbiomes of vaginally delivered (VD) and CSD neonates. Several functional pathways are over-represented in VD neonates, including lipopolysaccharide (LPS) biosynthesis. We link these enriched functions to individual-specific strains, which are transmitted from mothers to neonates in case of VD. The stimulation of primary human immune cells with LPS isolated from early stool samples of VD neonates results in higher levels of tumour necrosis factor (TNF-α) and interleukin 18 (IL-18). Accordingly, the observed levels of TNF-α and IL-18 in neonatal blood plasma are higher after VD. Taken together, our results support that CSD disrupts mother-to-neonate transmission of specific microbial strains, linked functional repertoires and immune-stimulatory potential during a critical window for neonatal immune system priming.

Background: Next generation sequencing (NGS) technologies can be applied in complex microbial ecosystems for metatranscriptome analysis by employing direct cDNA sequencing, which is known as RNA sequencing (RNA-seq). RNA-seq generates large datasets of great complexity, the comprehensive interpretation of which requires a reliable bioinformatic pipeline. In this study, we focus on the development of such a metatranscriptome pipeline, which we validate using Illumina RNA-seq datasets derived from the small intestine microbiota of two individuals with an ileostomy. Results: The metatranscriptome pipeline developed here enabled effective removal of rRNA derived sequences, followed by confident assignment of the predicted function and taxonomic origin of the mRNA reads. Phylogenetic analysis of the small intestine metatranscriptome datasets revealed a strong similarity with the community composition profiles obtained from 16S rDNA and rRNA pyrosequencing, indicating considerable congruency between community composition (rDNA), and the taxonomic distribution of overall (rRNA) and specific (mRNA) activity among its microbial members. Reproducibility of the metatranscriptome sequencing approach was established by independent duplicate experiments. In addition, comparison of metatranscriptome analysis employing single- or paired-end sequencing methods indicated that the latter approach does not provide improved functional or phylogenetic insights. Metatranscriptome functional-mapping allowed the analysis of global, and genus specific activity of the microbiota, and illustrated the potential of these approaches to unravel syntrophic interactions in microbial ecosystems. Conclusions: A reliable pipeline for metatransciptome data analysis was developed and evaluated using RNA-seq datasets obtained for the human small intestine microbiota. The set-up of the pipeline is very generic and can be applied for (bacterial) metatranscriptome analysis in any chosen niche.

Background: Massive high-throughput sequencing of short, hypervariable segments of the 16S ribosomal RNA (rRNA) gene has transformed the methodological landscape describing microbial diversity within and across complex biomes. However, several studies have shown that the methodology rather than the biological variation is responsible for the observed sample composition and distribution. This compromises meta-analyses, although this fact is often disregarded. Results: To facilitate true meta-analysis of microbiome studies, we developed NG-Tax, a pipeline for 16S rRNA gene amplicon sequence analysis that was validated with different mock communities and benchmarked against QIIME as a frequently used pipeline. The microbial composition of 49 independently amplified mock samples was characterized by sequencing two variable 16S rRNA gene regions, V4 and V5-V6, in three separate sequencing runs on Illumina’s HiSeq2000 platform. This allowed for the evaluation of important causes of technical bias in taxonomic classification: 1) run-to-run sequencing variation, 2) PCR–error, and 3) region/primer specific amplification bias. Despite the short read length (~140 nt) and all technical biases, the average specificity of the taxonomic assignment for the phylotypes included in the mock communities was 97.78%. On average 99.95% and 88.43% of the reads could be assigned to at least family or genus level, respectively, while assignment to ‘spurious genera’ represented on average only 0.21% of the reads per sample. Analysis of α- and β-diversity confirmed conclusions guided by biology rather than the aforementioned methodological aspects, which was not achieved with QIIME. Conclusions: Different biological outcomes are commonly observed due to 16S rRNA region-specific performance. NG-Tax demonstrated high robustness against choice of region and other technical biases associated with 16S rRNA gene amplicon sequencing studies, diminishing their impact and providing accurate qualitative and quantitative representation of the true sample composition. This will improve comparability between studies and facilitate efforts towards standardization.

By modulating the human gut microbiome, prebiotics and probiotics (combinations of which are called synbiotics) may be used to treat diseases such as colorectal cancer (CRC). Methodological limitations have prevented determining the potential combinatorial mechanisms of action of such regimens. We expanded our HuMiX gut-on-a-chip model to co-culture CRC-derived epithelial cells with a model probiotic under a simulated prebiotic regimen, and we integrated the multi-omic results with in silico metabolic modeling. In contrast to individual prebiotic or probiotic treatments, the synbiotic regimen caused downregulation of genes involved in procarcinogenic pathways and drug resistance, and reduced levels of the oncometabolite lactate. Distinct ratios of organic and short-chain fatty acids were produced during the simulated regimens. Treatment of primary CRC-derived cells with a molecular cocktail reflecting the synbiotic regimen attenuated self-renewal capacity. Our integrated approach demonstrates the potential of modeling for rationally formulating synbiotics-based treatments in the future.

The exposure of fish to environmental free-living microbes and its effect on early colonization in the gut have been studied in recent years. However, little is known regarding how the host and environment interact to shape gut communities during early life. Here, we tested whether the early microbial exposure of tilapia larvae affects the gut microbiota at later life stages. The experimental period was divided into three stages: axenic, probiotic and active suspension. Axenic tilapia larvae were reared either under conventional conditions (active suspension systems) or exposed to a single strain probiotic (Bacillus subtilis) added to the water. Microbial characterization by Illumina HiSeq sequencing of 16S rRNA gene amplicons showed the presence of B. subtilis in the gut during the seven days of probiotic application. Although B. subtilis was no longer detected in the guts of fish exposed to the probiotic after day 7, gut microbiota of the exposed tilapia larvae remained significantly different from that of the control treatment. Compared with the control, fish gut microbiota under probiotic treatment was less affected by spatial differences resulting from tank replication, suggesting that the early probiotic contact contributed to the subsequent observation of low inter-individual variation.