Piphillin: Improved Prediction of Metagenomic Content by Direct Inference from Human Microbiomes

The role of the gut microbiota in the dietary niche expansion of fishing bats

Background

Due to its central role in animal nutrition, the gut microbiota is likely a relevant factor shaping dietary niche shifts. We analysed both the impact and contribution of the gut microbiota to the dietary niche expansion of the only four bat species that have incorporated fish into their primarily arthropodophage diet.

Results

We first compared the taxonomic and functional features of the gut microbiota of the four piscivorous bats to that of 11 strictly arthropodophagous species using 16S rRNA targeted amplicon sequencing. Second, we increased the resolution of our analyses for one of the piscivorous bat species, namely Myotis capaccinii, and analysed multiple populations combining targeted approaches with shotgun sequencing. To better understand the origin of gut microorganisms, we also analysed the gut microbiota of their fish prey (Gambusia holbrooki). Our analyses showed that piscivorous bats carry a characteristic gut microbiota that differs from that of their strict arthropodophagous counterparts, in which the most relevant bacteria have been directly acquired from their fish prey. This characteristic microbiota exhibits enrichment of genes involved in vitamin biosynthesis, as well as complex carbohydrate and lipid metabolism, likely providing their hosts with an enhanced capacity to metabolise the glycosphingolipids and long-chain fatty acids that are particularly abundant in fish.

Conclusions

Our results depict the gut microbiota as a relevant element in facilitating the dietary transition from arthropodophagy to piscivory.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00137-w.

Cross-Talk Between Intestinal Microbiota and Host Gene Expression in Gilthead Sea Bream (Sparus aurata) Juveniles: Insights in Fish Feeds for Increased Circularity and Resource Utilization

New types of fish feed based on processed animal proteins (PAPs), insect meal, yeast, and microbial biomasses have been used with success in gilthead sea bream. However, some drawback effects on feed conversion and inflammatory systemic markers were reported in different degrees with PAP- and non-PAP-based feed formulations. Here, we focused on the effects of control and two experimental diets on gut mucosal-adherent microbiota, and how it correlated with host transcriptomics at the local (intestine) and systemic (liver and head kidney) levels. The use of tissue-specific PCR-arrays of 93 genes in total rendered 13, 12, and 9 differentially expressed (DE) genes in the intestine, liver, and head kidney, respectively. Illumina sequencing of gut microbiota yielded a mean of 125,350 reads per sample, assigned to 1,281 operational taxonomic unit (OTUs). Bacterial richness and alpha diversity were lower in fish fed with the PAP diet, and discriminant analysis displayed 135 OTUs driving the separation between groups with 43 taxa correlating with 27 DE genes. The highest expression of intestinal pcna and alpi was achieved in PAP fish with intermediate values in non-PAP, being the pro-inflammatory action of alpi associated with the presence of Psychrobacter piscatorii. The intestinal muc13 gene was down-regulated in non-PAP fish, with this gene being negatively correlated with anaerobic (Chloroflexi and Anoxybacillus) and metal-reducing (Pelosinus and Psychrosinus) bacteria. Other inflammatory markers (igm, il8, tnfα) were up-regulated in PAP fish, positively correlating the intestinal igm gene with the inflammasome activator Escherichia/Shigella, whereas the systemic expression of il8 and tnfα was negatively correlated with the Bacilli class in PAP fish and positively correlated with Paracoccus yeei in non-PAP fish. Overall changes in the expression pattern of il10, galectins (lgals1, lgals8), and toll-like receptors (tlr2, tlr5, tlr9) reinforced the anti-inflammatory profile of fish fed with the non-PAP diet, with these gene markers being associated with a wide range of OTUs. A gut microbiota-liver axis was also established, linking the microbial generation of short chain fatty acids with the fueling of scd1- and elovl6-mediated lipogenesis. In summary, by correlating the microbiome with host gene expression, we offer new insights in the evaluation of fish diets promoting gut and metabolism homeostasis, and ultimately, the health of farmed fish.

Effects of Dietary Fishmeal Replacement by Poultry By-Product Meal and Hydrolyzed Feather Meal on Liver and Intestinal Histomorphology and on Intestinal Microbiota of Gilthead Seabream (Sparus aurata)

The effects on liver and intestinal histomorphology and on intestinal microbiota in gilthead seabream (Sparus aurata) fed diets that contained poultry by-product meal (PBM) and hydrolyzed feather meal (HFM) as fishmeal replacements were studied. Fish fed on a series of isonitrogenous and isoenergetic diets, where fishmeal protein of the control diet (FM diet) was replaced by either PBM or by HFM at 25%, 50% and 100% without amino acid supplementation (PBM25, PBM50, PBM100, HFM25, HFM50 and HFM100 diets) or supplemented with lysine and methionine (PBM25+, PBM50+, HFM25+ and HFM50+ diets). The use of PBM and HFM at 25% fishmeal replacement generated a similar hepatic histomorphology to FM-fed fish, indicating that both land animal proteins are highly digestible at low FM replacement levels. However, 50% and 100% FM replacement levels by either PBM or HFM resulted in pronounced hepatic alterations in fish with the latter causing more severe degradation of the liver. Dietary amino acid supplementation delivered an improved tissue histology signifying their importance at high FM replacement levels. Intestinal microbiota was dominated by Proteobacteria (58.8%) and Actinobacteria (32.4%) in all dietary groups, but no specific pattern was observed among them at any taxonomic level. This finding was probably driven by the high inter-individual variability observed.

The effect of benzo[a]pyrene on the gut microbiota of Nile tilapia (Oreochromis niloticus)

Benzo[a]pyrene (BaP) is a highly toxic and carcinogenic polycyclic aromatic hydrocarbon (PAH) whose toxicological effects in the gut microbiota of aquatic organisms have not yet been fully revealed. Therefore, in this study, we used high-throughput 16S rRNA gene sequencing to evaluate the effects of BaP in the gut microbiome of Oreochromis niloticus, including its possible participation in the process of detoxification and its ability to recover. The fish were injected with a single intraperitoneal dose of 20 mg kg^-1 of BaP, and the effects in the microbiome were evaluated at 24, 72, and 120 h post-injection. The results indicate a clear dysbiosis (in composition, relative abundance, diversity, and interaction networks) of the gut microbiota during 24 h post-injection, dominated by Fusobacteria and Bacteroidetes and a decrease in Proteobacteria and Spirochaetae. Interestingly, a slight recovery of the microbiome begins at 72 h and stabilises at 120 h post-injection. Pathway analysis revealed the participation of the gut microbiome in PAH degradation mainly at 24 h post-injection. This study provides new insights in the toxicology of BaP in O. niloticus and the first evidence of the ability of the gut microbiome to recovery after a chemical disturbance. KEY POINTS: • Benzo[a]pyrene caused a dysbiosis in the gut microbiota of Oreochromis niloticus. • We observed an enrichment of bacteria involved in the metabolism of xenobiotics. • The gut microbiota was recovered after exposure to benzo[a]pyrene.

Social environment drives sex and age-specific variation in Drosophila melanogaster microbiome composition and predicted function

Social environments influence multiple traits of individuals including immunity, stress and ageing, often in sex-specific ways. The composition of the microbiome (the assemblage of symbiotic microorganisms within a host) is determined by environmental factors and the host's immune, endocrine and neural systems. The social environment could alter host microbiomes extrinsically by affecting transmission between individuals, probably promoting homogeneity in the microbiome of social partners. Alternatively, intrinsic effects arising from interactions between the microbiome and host physiology (the microbiota-gut-brain axis) could translate social stress into dysbiotic microbiomes, with consequences for host health. We investigated how manipulating social environments during larval and adult life-stages altered the microbiome composition of Drosophila melanogaster fruit flies. We used social contexts that particularly alter the development and lifespan of males, predicting that any intrinsic social effects on the microbiome would therefore be sex-specific. The presence of adult males during the larval stage significantly altered the microbiome of pupae of both sexes. In adults, same-sex grouping increased bacterial diversity in both sexes. Importantly, the microbiome community structure of males was more sensitive to social contact at older ages, an effect partially mitigated by housing focal males with young rather than coaged groups. Functional analyses suggest that these microbiome changes impact ageing and immune responses. This is consistent with the hypothesis that the substantial effects of the social environment on individual health are mediated through intrinsic effects on the microbiome, and provides a model for understanding the mechanistic basis of the microbiota-gut-brain axis.

Gut microbiota alterations in response to sleep length among African-origin adults

Sleep disorders are increasingly being characterized in modern society as contributing to a host of serious medical problems, including obesity and metabolic syndrome. Changes to the microbial community in the human gut have been reportedly associated with many of these cardiometabolic outcomes. In this study, we investigated the impact of sleep length on the gut microbiota in a large cohort of 655 participants of African descent, aged 25-45, from Ghana, South Africa (SA), Jamaica, and the United States (US). The sleep duration was self-reported via a questionnaire. Participants were classified into 3 sleep groups: short (<7hrs), normal (7-<9hrs), and long (≥9hrs). Forty-seven percent of US participants were classified as short sleepers and 88% of SA participants as long sleepers. Gut microbial composition analysis (16S rRNA gene sequencing) revealed that bacterial alpha diversity negatively correlated with sleep length (p<0.05). Furthermore, sleep length significantly contributed to the inter-individual beta diversity dissimilarity in gut microbial composition (p<0.01). Participants with both short and long-sleep durations exhibited significantly higher abundances of several taxonomic features, compared to normal sleep duration participants. The predicted relative proportion of two genes involved in the butyrate synthesis via lysine pathway were enriched in short sleep duration participants. Finally, co-occurrence relationships revealed by network analysis showed unique interactions among the short, normal and long duration sleepers. These results suggest that sleep length in humans may alter gut microbiota by driving population shifts of the whole microbiota and also specific changes in Exact Sequence Variants abundance, which may have implications for chronic inflammation associated diseases. The current findings suggest a possible relationship between disrupted sleep patterns and the composition of the gut microbiota. Prospective investigations in larger and more prolonged sleep researches and causally experimental studies are needed to confirm these findings, investigate the underlying mechanism and determine whether improving microbial homeostasis may buffer against sleep-related health decline in humans.

Microbiome and oral squamous cell carcinoma: a possible interplay on iron metabolism and its impact on tumor microenvironment

There is increasing evidence showing positive association between changes in oral microbiome and the occurrence of oral squamous cell carcinoma (OSCC). Alcohol- and nicotine-related products can induce microbial changes but are still unknown if these changes are related to cancerous lesion sites. In an attempt to understand how these changes can influence the OSCC development and maintenance, the aim of this study was to investigate the oral microbiome linked with OSCC as well as to identify functional signatures and associate them with healthy or precancerous and cancerous sites. Our group used data of oral microbiomes available in public repositories. The analysis included data of oral microbiomes from electronic cigarette users, alcohol consumers, and precancerous and OSCC samples. An R-based pipeline was used for taxonomic and functional prediction analysis. The Streptococcus spp. genus was the main class identified in the healthy group. Haemophilus spp. predominated in precancerous lesions. OSCC samples revealed a higher relative abundance compared with the other groups, represented by an increased proportion of Fusobacterium spp., Prevotella spp., Haemophilus spp., and Campylobacter spp. Venn diagram analysis showed 52 genera exclusive of OSCC samples. Both precancerous and OSCC samples seemed to present a specific associated functional pattern. They were menaquinone-dependent protoporphyrinogen oxidase pattern enhanced in the former and both 3',5'-cyclic-nucleotide phosphodiesterase (purine metabolism) and iron(III) transport system ATP-binding protein enhanced in the latter. We conclude that although precancerous and OSCC samples present some differences on microbial profile, both microbiomes act as "iron chelators-like" potentially contributing to tumor growth.

Identification of Biomolecules Involved in the Adaptation to the Environment of Cold-Loving Microorganisms and Metabolic Pathways for Their Production

Cold-loving microorganisms of all three domains of life have unique and special abilities that allow them to live in harsh environments. They have acquired structural and molecular mechanisms of adaptation to the cold that include the production of anti-freeze proteins, carbohydrate-based extracellular polymeric substances and lipids which serve as cryo- and osmoprotectants by maintaining the fluidity of their membranes. They also produce a wide diversity of pigmented molecules to obtain energy, carry out photosynthesis, increase their resistance to stress and provide them with ultraviolet light protection. Recently developed analytical techniques have been applied as high-throughoutput technologies for function discovery and for reconstructing functional networks in psychrophiles. Among them, omics deserve special mention, such as genomics, transcriptomics, proteomics, glycomics, lipidomics and metabolomics. These techniques have allowed the identification of microorganisms and the study of their biogeochemical activities. They have also made it possible to infer their metabolic capacities and identify the biomolecules that are parts of their structures or that they secrete into the environment, which can be useful in various fields of biotechnology. This Review summarizes current knowledge on psychrophiles as sources of biomolecules and the metabolic pathways for their production. New strategies and next-generation approaches are needed to increase the chances of discovering new biomolecules.

Essential Hypertension Is Associated With Changes in Gut Microbial Metabolic Pathways: A Multisite Analysis of Ambulatory Blood Pressure

[Figure: see text].

Microbiota from young mice counteracts selective age-associated behavioral deficits

The gut microbiota is increasingly recognized as an important regulator of host immunity and brain health. The aging process yields dramatic alterations in the microbiota, which is linked to poorer health and frailty in elderly populations. However, there is limited evidence for a mechanistic role of the gut microbiota in brain health and neuroimmunity during aging processes. Therefore, we conducted fecal microbiota transplantation from either young (3-4 months) or old (19-20 months) donor mice into aged recipient mice (19-20 months). Transplant of a microbiota from young donors reversed aging-associated differences in peripheral and brain immunity, as well as the hippocampal metabolome and transcriptome of aging recipient mice. Finally, the young donor-derived microbiota attenuated selective age-associated impairments in cognitive behavior when transplanted into an aged host. Our results reveal that the microbiome may be a suitable therapeutic target to promote healthy aging.

A large-scale metagenomic study for enzyme profiles using the focused identification of the NGS-based definitive enzyme research (FINDER) strategy

Excavating the molecular details of many diverse enzymes from metagenomes remains challenging in agriculture, food, health, and environmental fields. We present a versatile method that accelerates metabolic enzyme discovery for highly selective gene capture in metagenomes using next-generation sequencing. Culture-independent enzyme mining of environmental DNA is based on a set of short identifying degenerate sequences specific for a wide range of enzyme superfamilies, followed by multiplexed DNA barcode sequencing. A strategy of 'focused identification of next-generation sequencing-based definitive enzyme research' enabled us to generate targeted enzyme datasets from metagenomes, resulting in minimal hands-on obtention of high-throughput biological diversity and potential function profiles, without being time-consuming. This method also provided a targeted inventory of predicted proteins and molecular features of metabolic activities from several metagenomic samples. We suggest that the efficiency and sensitivity of this method will accelerate the decryption of microbial diversity and the signature of proteins and their metabolism from environmental samples.

Deficiency of interleukin-17 receptor A1 induces microbiota disruption in the intestine of Japanese medaka, Oryzias latipes

The mutual relationship between the intestinal immune system and the gut microbiota has received a great deal of attention. In mammals, interleukin-17A and F (IL-17A/F) are inflammatory cytokines and key regulators of the gut microbiota. However, in teleosts, the function of IL-17A/F in controlling the gut microbiota is poorly understood. We attempted to elucidate the importance of teleost IL-17 signaling in controlling gut microbiota. We previously established a knockout (KO) of IL-17 receptor A (RA) 1, a receptor for IL-17A/F, in the Japanese medaka (Oryzias latipes) using the CRISPR-Cas9 system and performed 16S rRNA-based metagenomic analyses using the anterior and posterior sections of the intestinal tract. The number of observed OTUs in the anterior intestine was significantly decreased in IL-17RA1 KO medaka compared to that in the wild-type (WT). Furthermore, β-diversity analysis (weighted UniFrac) revealed considerably different bacterial composition in the anterior intestine of IL-17RA1 KO compared to WT, with similar findings in α-diversity. Notably, the pathogen Plesiomonas shigelloides was significantly increased in the posterior intestine of IL-17RA1 KO medaka. These findings indicate that signaling via IL-17RA1 is required to maintain a healthy gut microbiota in teleosts and mammals. The involvement of IL-17RA1 in controlling the gut microbiota has been demonstrated, resulting in microbiome dysbiosis in IL-17RA1 KO medaka.

Transcriptome Profiling and Metagenomic Analysis Help to Elucidate Interactions in an Inflammation-Associated Cancer Mouse Model

Simple Summary

Colitis-associated colorectal cancer is the third most significant condition that increases the overall risk of developing colorectal cancer. In this study, we examined normal colonic mucosa of tumor-bearing mice in the DSS/AOM mouse model by gene expression profiling and fecal samples by 16s rDNA amplicon sequencing. Gene set enrichment analysis revealed that genes associated with fatty acid metabolism, oxidative phosphorylation, and the PI3K-Akt-mTOR pathways were enriched colonic mucosa of DSS/AOM mice. Additionally, enrichment of the sphingolipid signal and lipoarabinomannan biosynthetic pathways were inferred from fecal microbial composition. Our findings provide insights into altered transcriptome and microbiome in a mouse model of colitis-induced carcinogenesis.

Abstract

Colitis is a risk factor for colorectal cancer (CRC) and can change the dynamics of gut microbiota, leading to dysbiosis and contributing to carcinogenesis. The functional interactions between colitis-associated CRC and microbiota remain unknown. In this study, colitis and CRC were induced in BALB/c mice by the administration of dextran sodium sulfate (DSS) and/or azoxymethane (AOM). Whole transcriptome profiling of normal colon was then performed, and gene set enrichment analysis (GSEA) revealed enriched fatty acid metabolism, oxidative phosphorylation, and PI3K-Akt-mTOR signaling in the tissues from DSS/AOM mice. Additionally, immunohistochemical staining showed increased expression levels of phosphorylated S6 ribosomal protein, a downstream target of the PI3K-Akt-mTOR pathway in the inflamed mucosa of DSS/AOM mice. Fecal microbes were characterized using 16S rDNA gene sequencing. Redundancy analysis demonstrated a significant dissimilarity between the DSS/AOM group and the others. Functional analysis inferred from microbial composition showed enrichments of the sphingolipid signal and lipoarabinomannan biosynthetic pathways. This study provides additional insights into alterations associated with DSS/AOM-induced colitis and associates PI3K-Akt-mTOR, sphingolipid-signaling and lipoarabinomannan biosynthetic pathways in mouse DSS/AOM-induced colitis.

Metabarcoding under Brine: Microbial Ecology of Five Hypersaline Lakes at Rottnest Island (WA, Australia)

Hypersaline ecosystems—aquatic environments where concentration of salt exceeds 35 g L−1—host microbial communities that are highly specialised to cope with these extreme conditions. However, our knowledge on the taxonomic diversity and functional metabolisms characterising microbial communities in the water columns of hypersaline ecosystems is still limited, and this may compromise the future preservation of these unique environments. DNA metabarcoding provides a reliable and affordable tool to investigate environmental dynamics of aquatic ecosystems, and its use in brine can be highly informative. Here, we make use of bacterial 16S metabarcoding techniques combined with hydrochemical analyses to investigate the microbial patterns (diversity and functions) from five hypersaline lakes located at Rottnest Island (WA). Our results indicate lake-driven microbial aquatic assemblages that are characterised by taxonomically and functionally moderately to extremely halophilic groups, with TDS (total dissolved solids) and alkalinity amongst the most influential parameters driving the community patterns. Overall, our findings suggest that DNA metabarcoding allows rapid but reliable ecological assessment of the hypersaline aquatic microbial communities at Rottnest Island. Further studies involving different hypersaline lakes across multiple seasons will help elucidate the full extent of the potential of this tool in brine.

Investigation into the bacterial diversity of sediment samples obtained from Berg River, Western Cape, South Africa

This study used conventional culturing and 16S rRNA metagenomics analyses to assess the diversity of bacterial communities in sediment samples obtained from the Berg River, Western Cape, South Africa. Samples were collected from six points: a residential and recreational area, an industrial area, an informal residential settlement, a point next to a wastewater treatment plant (WWTP), a pumping station, and a residential and agricultural farming area along the river. High bacterial counts recorded on general selective and differential culture media signify substantial microbial contamination along the sampling sites. The most prevalent bacterial phyla detected (through metagenomics analyses) along the sampling sites were Proteobacteria (61%), Planctomycetes (9.5%), Firmicutes (7.8%), Bacteroidetes (5%), Acidobacteria (4.6%), and Actinobacteria (4.6%). Some members of the identified predominant bacterial phyla, genera, and classes are important public health bacteria that have been implicated in human diseases and outbreaks, while some others are metal or hydrocarbon tolerant, indicating possible significant environmental pollution. Notable human pathogenic genera such as Bacillus, Clostridium, Shigella, Legionella, Mycobacterium, and Pseudomonas were identified in varying percentages at five of the six sampling areas. Fecal contamination was particularly rife at all residential areas, with the informal housing area being the most notably polluted. Diverse functional pathways were predicted for identified bacteria, such as those associated with different chronic and infectious human diseases as well as those related to hydrocarbon and metal remediation. The point next to a WWTP contained vastly diverse groups of bacterial contaminants as well as the most abundant pathway identities and titles.

Symbiotic nitrogen fixation in the reproductive structures of a basidiomycete fungus

Nitrogen (N) fixation is a driving force for the formation of symbiotic associations between N₂-fixing bacteria and eukaryotes.¹ Limited examples of these associations are known in fungi, and none with sexual structures of non-lichenized species.^2-6 The basidiomycete Guyanagaster necrorhizus is a sequestrate fungus endemic to the Guiana Shield.⁷ Like the root rot-causing species in its sister genera Armillaria and Desarmillaria, G. necrorhizus sporocarps fruit from roots of decaying trees (Figures 1A-1C),⁸ and genome sequencing is consistent with observations that G. necrorhizus is a white-rotting decomposer. This species also represents the first documentation of an arthropod-dispersed sequestrate fungus. Numerous species of distantly related wood-feeding termites, which scavenge for N-rich food, feed on the mature spore-bearing tissue, or gleba, of G. necrorhizus. During feeding, mature spores adhere to termites for subsequent dispersal.⁹ Using chemical assays, isotope analysis, and high-throughput sequencing, we show that the sporocarps harbor actively N₂-fixing Enterobacteriaceae species and that the N content within fungal tissue increases with maturation. Untargeted proteomic profiling suggests that ATP generation in the gleba is accomplished via fermentation. The use of fermentation-an anaerobic process-indicates that the sporocarp environment is anoxic, likely an adaptation to protect the oxygen-sensitive nitrogenase enzyme. Sporocarps also have a thick outer covering, possibly to limit oxygen diffusion. The enriched N content within mature sporocarps may offer a dietary inducement for termites in exchange for spore dispersal. These results show that the flexible metabolic capacity of fungi may facilitate N₂-fixing associations, as well as higher-level organismal associations.

Sasso Pisano Geothermal Field Environment Harbours Diverse Ktedonobacteria Representatives and Illustrates Habitat-Specific Adaptations

The hydrothermal steam environment of Sasso Pisano (Italy) was selected to investigate the associated microbial community and its metabolic potential. In this context, 16S and 18S rRNA gene partial sequences of thermophilic prokaryotes and eukaryotes inhabiting hot springs and fumaroles as well as mesophilic microbes colonising soil and water were analysed by high-throughput amplicon sequencing. The eukaryotic and prokaryotic communities from hot environments clearly differ from reference microbial communities of colder soil sites, though Ktedonobacteria showed high abundances in various hot spring samples and a few soil samples. This indicates that the hydrothermal steam environments of Sasso Pisano represent not only a vast reservoir of thermophilic but also mesophilic members of this Chloroflexi class. Metabolic functional profiling revealed that the hot spring microbiome exhibits a higher capability to utilise methane and aromatic compounds and is more diverse in its sulphur and nitrogen metabolism than the mesophilic soil microbial consortium. In addition, heavy metal resistance-conferring genes were significantly more abundant in the hot spring microbiome. The eukaryotic diversity at a fumarole indicated high abundances of primary producers (unicellular red algae: Cyanidiales), consumers (Arthropoda: Collembola sp.), and endoparasite Apicomplexa (Gregarina sp.), which helps to hypothesise a simplified food web at this hot and extremely nutrient-deprived acidic environment.

Probiotic Yeasts and Vibrio anguillarum Infection Modify the Microbiome of Zebrafish Larvae

The host microbiome plays an essential role in health and disease. Microbiome modification by pathogens or probiotics has been poorly explored especially in the case of probiotic yeasts. Next-generation sequencing currently provides the best tools for their characterization. Debaryomyces hansenii 97 (D. hansenii 97) and Yarrowia lipolytica 242 (Y. lipolytica 242) are yeasts that protect wildtype zebrafish (Danio rerio) larvae against a Vibrio anguillarum (V. anguillarum) infection, increasing their survival rate. We investigate the effect of these microorganisms on the microbiome and neutrophil response (inflammation) in zebrafish larvae line Tg(Bacmpx:GFP)ⁱ¹¹⁴. We postulated that preinoculation of larvae with yeasts would attenuate the intestinal neutrophil response and prevent modification of the larval microbiome induced by the pathogen. Microbiome study was performed by sequencing the V3-V4 region of the 16S rRNA gene and prediction of metabolic pathways by Piphillin in conventionally raised larvae. Survival and the neutrophil response were both evaluated in conventional and germ-free conditions. V. anguillarum infection resulted in higher neutrophil number in the intestinal area compared to non-infected larvae in both conditions. In germ-free conditions, infected larvae pre-inoculated with yeasts showed fewer neutrophil numbers than infected larvae. In both conditions, only D. hansenii 97 increased the survival of infected larvae. Beta diversity of the microbiota was modified by V. anguillarum and both yeasts, compared to non-inoculated larvae. At 3 days post-infection, V. anguillarum modified the relative abundance of 10 genera, and pre-inoculation with D. hansenii 97 and Y. lipolytica 242 prevented the modification of 5 and 6 of these genera, respectively. Both yeasts prevent the increase of Ensifer and Vogesella identified as negative predictors for larval survival (accounting for 40 and 27 of the variance, respectively). In addition, yeast pre-inoculation prevents changes in some metabolic pathways altered by V. anguillarum’s infection. These results suggest that both yeasts and V. anguillarum can shape the larval microbiota configuration in the early developmental stage of D. rerio. Moreover, modulation of key taxa or metabolic pathways of the larval microbiome by yeasts can be associated with the survival of infected larvae. This study contributes to the understanding of yeast–pathogen–microbiome interactions, although further studies are needed to elucidate the mechanisms involved.

Endogenous bacteria inhabiting the Ophiocordyceps highlandensis during fruiting body development

Background

The genus Ophiocordyceps, which includes Ophiocordyceps sinensis, has been demonstrated to be one of the most valuable medicinal taxa. The low rate of larval infection and slow development that characterize the cultivation of this genus should be urgently addressed. To identify potential bioinoculants that stimulate the growth of Ophiocordyceps, O. highlandensis was selected as a model system, and a total of 72 samples were collected to systematically compare the microbial communities present during fruiting body development. By applying high-throughput 16S and ITS2 amplicon sequencing technology, the bacterial and fungal communities were identified in O. highlandensis and its surrounding soil, and the functional dynamics of the bacteria were explored.

Results

The results indicate that the most abundant bacteria across all the samples from O. highlandensis were Proteobacteria, Firmicutes and Bacteroidetes, while members of Ascomycota were detected among the fungi. The pathways enriched in the developmental stages were associated with carbohydrate degradation, nucleotides and pyridoxal biosynthesis, and the TCA cycle. Compared with that in the fungal community, an unexpectedly high taxonomic and functional fluctuation was discovered in the bacterial community during the maturation of O. highlandensis. Furthermore, bipartite network analysis identified four potential supercore OTUs associated with O. highlandensis growth.

Conclusions

All the findings of this study suggest unexpectedly high taxonomic and functional fluctuations in the bacterial community of O. highlandensis during its maturation. O. highlandensis may recruit different endogenous bacteria across its life cycle to enhance growth and support rapid infection. These results may facilitate Ophiocordyceps cultivation and improve the development of strategies for the identification of potential bioinoculant resources.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02227-w.

Changes in the predicted function of the rumen bacterial community of Japanese Black beef cattle during the fattening stages according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses

We investigated changes in the predicted functions of the rumen bacterial community in Japanese Black beef cattle during fattening. Nine cattle were fed a high-concentrate diet during the early, middle, and late fattening stages consecutively (10–14, 15–22, and 23–30 months of age, respectively). The rumen fluid and solid samples collected at each stage were subjected to sequencing analyses. The sequencing results were clustered and classified into operational taxonomic units (OTUs). Representative sequences and a raw counting table for each OTU were submitted to the Piphillin website. The predicted functions were revealed by the Kyoto Encyclopedia of Genes and Genomes database as the ratio of the total sequence. In the early stage, “Biosynthesis of secondary metabolites” was significantly higher in the fluid fraction than in the solid fraction. “Two-component system” in the middle stage was significantly lower and “Purine metabolism” in the late stage was significantly higher in the fluid fraction than those in the solid fraction. The fluid fraction was significantly correlated with acetic acid, propionic acid, and bacterial metabolism, such as “Biosynthesis of secondary metabolites” and “Sugar metabolism.” Moreover, the solid fraction was correlated with “Purine metabolism” and “Biosynthesis of secondary metabolism”. These results suggest that the rumen bacterial community in Japanese Black beef cattle adapts to changes in rumen conditions by altering their functions in response to a long-term high-grain diet.

Influence of seasonal changes and salinity on spinach phyllosphere bacterial functional assemblage

The phyllosphere is the aerial part of plants that is exposed to different environmental conditions and is also known to harbor a wide variety of bacteria including both plant and human pathogens. However, studies on phyllosphere bacterial communities have focused on bacterial composition at different stages of plant growth without correlating their functional capabilities to bacterial communities. In this study, we examined the seasonal effects and temporal variabilities driving bacterial community composition and function in spinach phyllosphere due to increasing salinity and season and estimated the functional capacity of bacterial community16S V4 rRNA gene profiles by indirectly inferring the abundance of functional genes based on metagenomics inference tool Piphillin. The experimental design involved three sets of spinach (Spinacia oleracea L., cv. Racoon) grown with saline water during different seasons. Total bacteria DNA from leaf surfaces were sequenced using MiSeq® Illumina platform. About 66.35% of bacteria detected in the phyllosphere were dominated by four phyla- Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria. Permutational analysis of variance (PERMANOVA) showed that phyllosphere microbiomes were significantly (P < 0.003) affected by season, but not salinity (P = 0.501). The most abundant inferred functional pathways in leaf samples were the amino acids biosynthesis, ABC transporters, ribosome, aminoacyl-tRNA biosynthesis, two-component system, carbon metabolism, purine metabolism, and pyrimidine metabolism. The photosynthesis antenna proteins pathway was significantly enriched in June leaf samples, when compared to March and May. Several genes related to toxin co-regulated pilus biosynthesis proteins were also significantly enriched in June leaf samples, when compared to March and May leaf samples. Therefore, planting and harvesting times must be considered during leafy green production due to the influence of seasons in growth and proliferation of phyllosphere microbial communities.

General Unified Microbiome Profiling Pipeline (GUMPP) for Large Scale, Streamlined and Reproducible Analysis of Bacterial 16S rRNA Data to Predicted Microbial Metagenomes, Enzymatic Reactions and Metabolic Pathways

General Unified Microbiome Profiling Pipeline (GUMPP) was developed for large scale, streamlined and reproducible analysis of bacterial 16S rRNA data and prediction of microbial metagenomes, enzymatic reactions and metabolic pathways from amplicon data. GUMPP workflow introduces reproducible data analyses at each of the three levels of resolution (genus; operational taxonomic units (OTUs); amplicon sequence variants (ASVs)). The ability to support reproducible analyses enables production of datasets that ultimately identify the biochemical pathways characteristic of disease pathology. These datasets coupled to biostatistics and mathematical approaches of machine learning can play a significant role in extraction of truly significant and meaningful information from a wide set of 16S rRNA datasets. The adoption of GUMPP in the gut-microbiota related research enables focusing on the generation of novel biomarkers that can lead to the development of mechanistic hypotheses applicable to the development of novel therapies in personalized medicine.

Beyond taxonomy: Validating functional inference approaches in the context of fish-farm impact assessments

Characterization of microbial assemblages via environmental DNA metabarcoding is increasingly being used in routine monitoring programs due to its sensitivity and cost-effectiveness. Several programs have recently been developed which infer functional profiles from 16S rRNA gene data using hidden-state prediction (HSP) algorithms. These might offer an economic and scalable alternative to shotgun metagenomics. To date, HSP-based methods have seen limited use for benthic marine surveys and their performance in these environments remains unevaluated. In this study, 16S rRNA metabarcoding was applied to sediment samples collected at 0 and ≥1,200 m from Norwegian salmon farms, and three metabolic inference approaches (Paprica, Picrust2 and Tax4Fun2) evaluated against metagenomics and environmental data. While metabarcoding and metagenomics recovered a comparable functional diversity, the taxonomic composition differed between approaches, with genera richness up to 20× higher for metabarcoding. Comparisons between the sensitivity (highest true positive rates) and specificity (lowest true negative rates) of HSP-based programs in detecting functions found in metagenomic data ranged from 0.52 and 0.60 to 0.76 and 0.79, respectively. However, little correlation was observed between the relative abundance of their specific functions. Functional beta-diversity of HSP-based data was strongly associated with that of metagenomics (r ≥ 0.86 for Paprica and Tax4Fun2) and responded similarly to the impact of fish farm activities. Our results demonstrate that although HSP-based metabarcoding approaches provide a slightly different functional profile than metagenomics, partly due to recovering a distinct community, they represent a cost-effective and valuable tool for characterizing and assessing the effects of fish farming on benthic ecosystems.

Taxonomic Characterization and Microbial Activity Determination of Cold-Adapted Microbial Communities in Lava Tube Ice Caves from Lava Beds National Monument, a High-Fidelity Mars Analogue Environment

Martian lava tube caves resulting from a time when the planet was still volcanically active are proposed to contain deposits of water ice, a feature that may increase microbial habitability. In this study, we taxonomically characterized and directly measured metabolic activity of the microbial communities that inhabit lava tube ice from Lava Beds National Monument, an analogue environment to martian lava tubes. We investigated whether this environment was habitable to microorganisms by determining their taxonomic diversity, metabolic activity, and viability using both culture-dependent and culture-independent techniques. With 16S rRNA gene sequencing, we recovered 27 distinct phyla from both ice and ice-rock interface samples, primarily consisting of Actinobacteria, Proteobacteria, Bacteroidetes, Firmicutes, and Chloroflexi. Radiorespiration and Biolog EcoPlate assays found these microbial communities to be metabolically active at both 5°C and -5°C and able to metabolize diverse sets of heterotrophic carbon substrates at each temperature. Viable cells were predominantly cold adapted and capable of growth at 5°C (1.3 × 10⁴ to 2.9 × 10⁷ cells/mL), and 24 of 38 cultured isolates were capable of growth at -5°C. Furthermore, 14 of these cultured isolates, and 16 of the 20 most numerous amplicon sequences we recovered were most closely related to isolates and sequences obtained from other cryophilic environments. Given these results, lava tube ice appears to be a habitable environment, and considering the protections martian lava tubes offer to microbial communities from harsh surface conditions, similar martian caves containing ice may be capable of supporting extant, active microbial communities.

Intestinal Microbiome in Preterm Infants Influenced by Enteral Iron Dosing

OBJECTIVES

The aim of the study was to compare the intestinal microbiome in very low birth weight (VLBW) infants who received different enteral iron supplementation (EIS) doses.

STUDY DESIGN

Longitudinal stool collection in 80 VLBW infants were conducted up to 2 months postnatally in a prospective study. The 16S rRNA regions V4 was used to calculate microbiome compositions and the Piphillin software was used for bacterial functional prediction. Linear mixed effect models and Wilcoxon rank-sum tests were performed to examine the relationships between initial EIS dosage and stool microbiome and bacterial functional potential.

RESULTS

There were 105 samples collected before and 237 collected after EIS started from infants with birth gestational age and weight of 28.1 ± 2.4 weeks and 1103 ± 210 g, respectively. The average postnatal age at start of EIS was 17.9 ± 6.9 days and the average initial EIS dose was 4.8 ± 1.1 mg · kg-1 · day-1. Infants who were started on ≥6 mg · kg-1 · day-1 had higher abundances of Proteus and Bifidobacterium and a lower alpha diversity than those started on lower doses (P < 0.05). Infants given higher EIS doses had higher bacterial predicted functional potentials for ferroptosis and epithelial invasion after 2 weeks post EIS.

CONCLUSIONS

Higher EIS dosage is linked to higher abundances of Proteus and Bifidobacterium, and a less diverse microbiome and higher predicted potential of bacterial epithelial invasion. These observational findings should be further studied in a randomized study to elucidate the optimal dosage of EIS in VLBW infants.

From the Andes to the desert: 16S rRNA metabarcoding characterization of aquatic bacterial communities in the Rimac river, the main source of water for Lima, Peru

The Rimac river is the main source of water for Lima, Peru's capital megacity. The river is constantly affected by different types of contamination including mine tailings in the Andes and urban sewage in the metropolitan area. In this work, we aim to produce the first characterization of aquatic bacterial communities in the Rimac river using a 16S rRNA metabarcoding approach which would be useful to identify bacterial diversity and potential understudied pathogens. We report a lower diversity in bacterial communities from the Lower Rimac (Metropolitan zone) in comparison to other sub-basins. Samples were generally grouped according to their geographical location. Bacterial classes Alphaproteobacteria, Bacteroidia, Campylobacteria, Fusobacteriia, and Gammaproteobacteria were the most frequent along the river. Arcobacter cryaerophilus (Campylobacteria) was the most frequent species in the Lower Rimac while Flavobacterium succinicans (Bacteroidia) and Hypnocyclicus (Fusobacteriia) were the most predominant in the Upper Rimac. Predicted metabolic functions in the microbiota include bacterial motility and quorum sensing. Additional metabolomic analyses showed the presence of some insecticides and herbicides in the Parac-Upper Rimac and Santa Eulalia-Parac sub-basins. The dominance in the Metropolitan area of Arcobacter cryaerophilus, an emergent pathogen associated with fecal contamination and antibiotic multiresistance, that is not usually reported in traditional microbiological quality assessments, highlights the necessity to apply next-generation sequencing tools to improve pathogen surveillance. We believe that our study will encourage the integration of omics sciences in Peru and its application on current environmental and public health issues.

Children with Autism and Their Typically Developing Siblings Differ in Amplicon Sequence Variants and Predicted Functions of Stool-Associated Microbes

The existence of a link between the gut microbiome and autism spectrum disorder (ASD) is well established in mice, but in human populations, efforts to identify microbial biomarkers have been limited due to a lack of appropriately matched controls, stratification of participants within the autism spectrum, and sample size. To overcome these limitations, we crowdsourced the recruitment of families with age-matched sibling pairs between 2 and 7 years old (within 2 years of each other), where one child had a diagnosis of ASD and the other did not. Parents collected stool samples, provided a home video of their ASD child's natural social behavior, and responded online to diet and behavioral questionnaires. 16S rRNA V4 amplicon sequencing of 117 samples (60 ASD and 57 controls) identified 21 amplicon sequence variants (ASVs) that differed significantly between the two cohorts: 11 were found to be enriched in neurotypical children (six ASVs belonging to the Lachnospiraceae family), while 10 were enriched in children with ASD (including Ruminococcaceae and Bacteroidaceae families). Summarizing the expected KEGG orthologs of each predicted genome, the taxonomic biomarkers associated with children with ASD can use amino acids as precursors for butyragenic pathways, potentially altering the availability of neurotransmitters like glutamate and gamma aminobutyric acid (GABA).IMPORTANCE Autism spectrum disorder (ASD), which now affects 1 in 54 children in the United States, is known to have comorbidity with gut disorders of a variety of types; however, the link to the microbiome remains poorly characterized. Recent work has provided compelling evidence to link the gut microbiome to the autism phenotype in mouse models, but identification of specific taxa associated with autism has suffered replicability issues in humans. This has been due in part to sample size that sufficiently covers the spectrum of phenotypes known to autism (which range from subtle to severe) and a lack of appropriately matched controls. Our original study proposes to overcome these limitations by collecting stool-associated microbiome on 60 sibling pairs of children, one with autism and one neurotypically developing, both 2 to 7 years old and no more than 2 years apart in age. We use exact sequence variant analysis and both permutation and differential abundance procedures to identify 21 taxa with significant enrichment or depletion in the autism cohort compared to their matched sibling controls. Several of these 21 biomarkers have been identified in previous smaller studies; however, some are new to autism and known to be important in gut-brain interactions and/or are associated with specific fatty acid biosynthesis pathways.

An Introduction to Next Generation Sequencing Bioinformatic Analysis in Gut Microbiome Studies

The gut microbiome is a microbial ecosystem which expresses 100 times more genes than the human host and plays an essential role in human health and disease pathogenesis. Since most intestinal microbial species are difficult to culture, next generation sequencing technologies have been widely applied to study the gut microbiome, including 16S rRNA, 18S rRNA, internal transcribed spacer (ITS) sequencing, shotgun metagenomic sequencing, metatranscriptomic sequencing and viromic sequencing. Various software tools were developed to analyze different sequencing data. In this review, we summarize commonly used computational tools for gut microbiome data analysis, which extended our understanding of the gut microbiome in health and diseases.

Breast milk flora plays an important role in infantile eczema: cohort study in Northeast China

AIMS

Infantile eczema, usually coupled with a range of hypersensitive phenotypes, has come into notice with its rising prevalence and unclear pathogenesis. Recent studies show close ties between eczema and an infant's intestinal flora. To gain a further understanding of the interactions between microbiota and eczema, we studied the breast milk flora as a new factor and present the links among breast milk flora, infant intestinal flora and infantile eczema through a cohort study in Northeast China.

METHODS AND RESULTS

Fifty-two families were recruited with either an eczema or healthy infant younger than 6 months. Analysis and predictions using amplicon sequencing of microbiota found that Bifidobacterium and Bacteroidetes were enriched in healthy and eczema infant stools, respectively, consistent with previous reports. For breast milk flora, more 'positive' bacteria such as Akkermansia were enriched in breast milk from healthy infants' mothers. Further, higher bacterial delivery efficiencies were found in pairs of breast milk flora and infants' stool flora of families with eczema infants compared with families with healthy infants. Bacteroidetes, a widely known indicator of eczema, was found delivered more in eczema pairs. Further metagenomic predictions revealed that the breast milk microbiota participated significantly less in metabolism and immune system pathways, particularly in antigen processing and presentation and in Th17 cell-related pathways.

CONCLUSIONS

In conclusion, as with other components of breast milk, the breast milk microbiota closely associates with infants' health via mother-infant bacterial delivery and metabolic functions.

SIGNIFICANCE AND IMPACT OF THE STUDY

Our research aimed to fill the gap between the eczema and breast milk flora and describe the connections among breast milk and intestinal flora and eczema.

Mining microbes for mental health: Determining the role of microbial metabolic pathways in human brain health and disease

There is increasing knowledge regarding the role of the microbiome in modulating the brain and behaviour. Indeed, the actions of microbial metabolites are key for appropriate gut-brain communication in humans. Among these metabolites, short-chain fatty acids, tryptophan, and bile acid metabolites/pathways show strong preclinical evidence for involvement in various aspects of brain function and behaviour. With the identification of neuroactive gut-brain modules, new predictive tools can be applied to existing datasets. We identified 278 studies relating to the human microbiota-gut-brain axis which included sequencing data. This spanned across psychiatric and neurological disorders with a small number also focused on normal behavioural development. With a consistent bioinformatics pipeline, thirty-five of these datasets were reanalysed from publicly available raw sequencing files and the remainder summarised and collated. Among the reanalysed studies, we uncovered evidence of disease-related alterations in microbial metabolic pathways in Alzheimer's Disease, schizophrenia, anxiety and depression. Amongst studies that could not be reanalysed, many sequencing and technical limitations hindered the discovery of specific biomarkers of microbes or metabolites conserved across studies. Future studies are warranted to confirm our findings. We also propose guidelines for future human microbiome analysis to increase reproducibility and consistency within the field.

Zebrafish and water microbiome recovery after oxytetracycline exposure

Oxytetracycline (OTC) is a broad-spectrum antibiotic widely used in aquaculture, resulting in contamination of aquatic environments. In a previous study, we observed significant effects of OTC sublethal concentrations in zebrafish, its microbiome and the water bacterial community. Here we assessed the extent to which these effects are reversible after a recovery period. Zebrafish adults were exposed to OTC (10,000 μg/L) via water exposure. Effects were analyzed at 5 days (5 dE) and 2 months (2 mE) of exposure and recovery was assessed at 5 days (5dPE) and 1 month (1mPE) after exposure Impacts were observed in fish energetic reserves and in fish and water microbiomes structure, being significant even at 5 dE. At energetic reserves level, the effect in cellular energy allocation (CEA) was dependent on the exposure time: initially CEA increased while after 2 mE CEA decreased. At microbiome level, diversity was not affected but the richness of the water microbiome significantly decreased at 2 mE. Regarding the post-exposure period, at CEA level, organisms seem to recover. In water and gut microbiomes OTC effects were also attenuated after exposure ceases, indicating a recovery. Even so, the structure of water exposed community remained significantly different towards the control, while richness of this community significantly increased at 1mPE. During exposure the relative abundance of 11 and 16 genera was significantly affected in the gut and water microbiomes, respectively, though these numbers decreased to 4 and 8 genera in the post-exposure period. At functional level during exposure 12 and 13 pathways were predicted to be affected in zebrafish gut and water microbiomes respectively, while post-exposure few pathways remained significantly affected. Hence, our results suggest a recovery of the fish fitness as well as of the water and intestine microbiomes after exposure ceases. Even so, some of the effects caused by OTC remain significant after this recovery period.

Volatility as a Concept to Understand the Impact of Stress on the Microbiome

The microbiome-gut-brain-axis is a complex phenomenon spanning several dynamic systems in the body which can be parsed at a molecular, cellular, physiological and ecological level. A growing body of evidence indicates that this axis is particularly sensitive to the effects of stress and that it may be relevant to stress resilience and susceptibility. Although stress-induced changes in the composition of the microbiome have been reported, the degree of compositional change over time, which we define as volatility, has not been the subject of in-depth scrutiny. Using a chronic psychosocial stress paradigm in male mice, we report that the volatility of the microbiome significantly correlated with several readouts of the stress response, including behaviour and corticosterone response. We then validated these findings in a second independent group of stressed mice. Additionally, we assessed the relationship between volatility and stress parameters in a cohort of health volunteers who were undergoing academic exams and report similar observations. Finally, we found inter-species similarities in the microbiome stress response on a functional level. Our research highlights the effects of stress on the dynamic microbiome and underscores the informative value of volatility as a parameter that should be considered in all future analyses of the microbiome.

Maternal Diet Shapes the Breast Milk Microbiota Composition and Diversity: Impact of Mode of Delivery and Antibiotic Exposure

BACKGROUND

Breast milk is a complex biofluid that provides nutrients and bioactive agents, including bacteria, for the development of the infant gut microbiota. However, the impact of maternal diet and other factors, such as mode of delivery and antibiotic exposure, on the breast milk microbiota has yet to be understood.

OBJECTIVES

This study aimed to examine the association between maternal diet and breast milk microbiota and to ascertain the potential role of mode of delivery and antibiotic exposure.

METHODS

In a cross-sectional study of the MAMI cohort, breast milk microbiota profiling was assessed in 120 samples from healthy mothers by 16S rRNA gene sequencing. Maternal dietary information was recorded through an FFQ, and clinical characteristics, including mode of delivery, antibiotic exposure, and exclusive breastfeeding, were collected.

RESULTS

Maternal diet was grouped into 2 clusters: Cluster I (high intake of plant protein, fiber, and carbohydrates), and Cluster II (high intake of animal protein and lipids). Breast milk microbiota was shaped by maternal dietary clusters. Staphylococcus and Bifidobacterium were associated with carbohydrate intake whereas the Streptococcus genus was associated with intakes of the n-3 PUFAs [EPA and docosapentaenoic acid (22:5ω-3)]. Mode of delivery and antibiotic exposure influenced breast milk microbiota in a diet cluster-dependent manner. Differences between/among the maternal dietary clusters were found in the milk microbiota of the cesarean-section (C-section)/antibiotic group, whereas no differences were observed in vaginal births. Lower abundances of Lactobacillus, Bacteroides, and Sediminibacterium genera were observed in Cluster II/C-section/antibiotic exposure compared with the other groups.

CONCLUSIONS

Maternal diet shapes the composition and diversity of breast milk microbiota, with the most important contributions coming from dietary fiber and both plant and animal protein intakes. The relation between the maternal diet and the milk microbiota needs further research because it has a key impact on infant microbiota development and contributes to infant health outcomes in the short and long term.This trial was registered at clinicaltrials.gov as NCT03552939.

Fermented Soybean Meal Replacement in the Diet of Lactating Holstein Dairy Cows: Modulated Rumen Fermentation and Ruminal Microflora

This study was conducted to examine the influences of replacing soybean meal (SBM) with fermented soybean meal (FSBM) in the diet of lactating Holstein cattle on rumen fermentation and ruminal bacterial microbiome. Twenty-four lactating Chinese Holstein dairy cattle were assigned to each of the two treatments in a completely randomized design: the SBM group [the basal total mixed ration (TMR) diet containing 5.77% SBM] and the FSBM group (the experimental TMR diet containing 5.55% FSBM). This trial lasted for 54 days (14 days for adjustment and 40 days for data and sample collection), and samples of rumen liquid were collected on 34 d and 54 d, respectively. The results showed that replacing SBM with FSBM significantly increased the molar percentages of propionate (P < 0.01) and valerate (P < 0.05), but reduced the total volatile fatty acid (TVFA) concentration (P < 0.05), butyrate molar proportion (P < 0.05), and the acetate to propionate ratio (P < 0.01). The copy numbers of total bacteria (P < 0.05), Fibrobacter succinogenes (P < 0.01), Selenomonas ruminantium (P < 0.01), and Prevotella spp. (P < 0.05) in the FSBM group were greater, while the density of Prevotella ruminicola (P < 0.05) was lower than those in the SBM treatment. Additionally, Succiniclasticum ruminis and Saccharofermentans acetigenes were significantly enriched (P < 0.05) in the rumen fluid of FSBM-fed cows, despite the fact that there was no remarkable difference in the Alpha diversity indexes, structure and KEGG pathway abundances of the bacterial community across the two treatments. It could hence be concluded that the substitution of FSBM for SBM modulated rumen fermentation and rumen bacterial microbiota in lactating Holstein dairy cows. Further research is required to elucidate the relevant mechanisms of FSBM, and provide more insights into the application of FSBM in dairy cattle.

Lignin intermediates lead to phenyl acid formation and microbial community shifts in meso- and thermophilic batch reactors

BACKGROUND

Lignin intermediates resulting from lignocellulose degradation have been suspected to hinder anaerobic mineralisation of organic materials to biogas. Phenyl acids like phenylacetate (PAA) are early detectable intermediates during anaerobic digestion (AD) of aromatic compounds. Studying the phenyl acid formation dynamics and concomitant microbial community shifts can help to understand the microbial interdependencies during AD of aromatic compounds and may be beneficial to counteract disturbances.

RESULTS

The length of the aliphatic side chain and chemical structure of the benzene side group(s) had an influence on the methanogenic system. PAA, phenylpropionate (PPA), and phenylbutyrate (PBA) accumulations showed that the respective lignin intermediate was degraded but that there were metabolic restrictions as the phenyl acids were not effectively processed. Metagenomic analyses confirmed that mesophilic genera like Fastidiosipila or Syntrophomonas and thermophilic genera like Lactobacillus, Bacillus, Geobacillus, and Tissierella are associated with phenyl acid formation. Acetoclastic methanogenesis was prevalent in mesophilic samples at low and medium overload conditions, whereas Methanoculleus spp. dominated at high overload conditions when methane production was restricted. In medium carbon load reactors under thermophilic conditions, syntrophic acetate oxidation (SAO)-induced hydrogenotrophic methanogenesis was the most important process despite the fact that acetoclastic methanogenesis would thermodynamically be more favourable. As acetoclastic methanogens were restricted at medium and high overload conditions, syntrophic acetate oxidising bacteria and their hydrogenotrophic partners could step in for acetate consumption.

CONCLUSIONS

PAA, PPA, and PBA were early indicators for upcoming process failures. Acetoclastic methanogens were one of the first microorganisms to be impaired by aromatic compounds, and shifts to syntrophic acetate oxidation coupled to hydrogenotrophic methanogenesis occurred in thermophilic reactors. Previously assumed associations of specific meso- and thermophilic genera with anaerobic phenyl acid formation could be confirmed.

Alterations of stool metabolome, phenome, and microbiome of the marine fish, red sea bream, Pagrus major, following exposure to phenanthrene: A non-invasive approach for exposure assessment

The present study aimed to assess the impact of phenanthrene (Phe) on fish health by addressing the alteration of fecal characteristics, in lieu of collecting biomarkers that often involves injurious or even fatal sampling of organisms. The marine fish red sea bream, Pagrus major, was exposed to Phe at a concentration of 18 μg/L for 16 days followed by depuration for 13 days. We collected feces from Phe-exposed or control (Phe-free) fish and then analyzed the fecal metabolite profile (metabolome), carbon utilization of microbiota (phenome), and bacterial 16s rRNA gene sequence (microbiome). Along with the increase in physiological stress markers (SOD and EROD) in serum and liver, we noted the possible role of intestine as a Phe reservoir. Furthermore, abnormal fecal appearance (green coloration) and remarkable changes in fecal characteristics were observed. These changes include alterations of cholesterol and putrescine metabolism and the enhanced utilization of putrescine as a carbon source. Phe also altered the microbial community, with an increase in Phe-degrading bacteria such as Pseudomonas. Interestingly, these enteric impairments were ameliorated by depuration. Taken together, our findings suggest that these alterations in feces were associated with adaptive responses to environmentally relevant Phe exposure scenarios, and that stool samples are potential candidates for exposure assessment in fish.

MetGEMs Toolbox: Metagenome-scale models as integrative toolbox for uncovering metabolic functions and routes of human gut microbiome

Investigating metabolic functional capability of a human gut microbiome enables the quantification of microbiome changes, which can cause a phenotypic change of host physiology and disease. One possible way to estimate the functional capability of a microbial community is through inferring metagenomic content from 16S rRNA gene sequences. Genome-scale models (GEMs) can be used as scaffold for functional estimation analysis at a systematic level, however up to date, there is no integrative toolbox based on GEMs for uncovering metabolic functions. Here, we developed the MetGEMs (metagenome-scale models) toolbox, an open-source application for inferring metabolic functions from 16S rRNA gene sequences to facilitate the study of the human gut microbiome by the wider scientific community. The developed toolbox was validated using shotgun metagenomic data and shown to be superior in predicting functional composition in human clinical samples compared to existing state-of-the-art tools. Therefore, the MetGEMs toolbox was subsequently applied for annotating putative enzyme functions and metabolic routes related in human disease using atopic dermatitis as a case study.

The rectal mucosal but not fecal microbiota detects subclinical ulcerative colitis

Ulcerative colitis (UC), a subtype of inflammatory bowel disease, is characterized by repetitive remission and relapse. Gut microbiome is critically involved in pathogenesis of UC. The shifts in microbiome profile during disease remission remain under-investigated. Recent studies revealed that UC pathogenesis is likely to originate in the mucosal barrier. Therefore, we investigated the effectiveness of mucosal tissue microbiomes to differentiate patients with subclinical UC from healthy individuals. The microbiomes of cecal and rectal biopsies and feces were characterized from 13 healthy individuals and 45 patients with subclinical UC. Total genomic DNA was extracted from the samples, and their microbial communities determined using next-generation sequencing. We found that changes in relative abundance of subclinical UC were marked by a decrease in Proteobacteria and an increase in Bacteroidetes phyla in microbiome derived from rectal tissues but not cecal tissue nor feces. Only in the microbiome of rectal tissue had significantly higher community richness and evenness in subclinical UC patients than controls. Twenty-seven operational taxonomic units were enriched in subclinical UC cohort with majority of the taxa from the Firmicutes phylum. Inference of putative microbial functional pathways from rectal biopsy microbiome suggested a differential increase in interleukin-17 signaling and T-helper cell differentiation pathways. Rectal biopsy tissue was suggested to be more suitable than fecal samples for microbiome assays to distinguish patients with subclinical UC from healthy adults. Assessment of the rectal biopsy microbiome may offer clinical insight into UC disease progression and predict relapse of the diseases.

Composition and geographic variation of the bacterial microbiota associated with the coelomic fluid of the sea urchin Paracentrotus lividus

In the present work, culture-based and culture-independent investigations were performed to determine the microbiota structure of the coelomic fluid of Mediterranean sea urchin Paracentrotus lividus individuals collected from two distinct geographical sites neighboring a high-density population bay and a nature reserve, respectively. Next Generation Sequencing analysis of 16S rRNA gene (rDNA) showed that members of the Proteobacteria, Bacteroidetes and Fusobacteria phyla, which have been previously reported to be commonly retrieved from marine invertebrates, dominate the overall population of microorganisms colonizing this liquid tissue, with minority bacterial genera exhibiting remarkable differences among individuals. Our results showed that there is a correlation between microbiota structure and geographical location of the echinoderm collection site, highlighting over-representation of metagenomic functions related to amino acid and bioactive peptides metabolism in specimens inhabiting the nature reserve. Finally, we also described the developmental delay and aberrations exhibited by sea urchin embryos exposed to distinct bacterial isolates, and showed that these defects rely upon hydrophilic compound(s) synthesized by the bacterial strains assayed. Altogether, our findings lay the groundwork to decipher the relationships of bacteria with sea urchins in their aquatic environment, also providing an additional layer of information to understand the biological roles of the coelomic fluid.

A Comparative Evaluation of Tools to Predict Metabolite Profiles From Microbiome Sequencing Data

Metabolomic analyses of human gut microbiome samples can unveil the metabolic potential of host tissues and the numerous microorganisms they support, concurrently. As such, metabolomic information bears immense potential to improve disease diagnosis and therapeutic drug discovery. Unfortunately, as cohort sizes increase, comprehensive metabolomic profiling becomes costly and logistically difficult to perform at a large scale. To address these difficulties, we tested the feasibility of predicting the metabolites of a microbial community based solely on microbiome sequencing data. Paired microbiome sequencing (16S rRNA gene amplicons, shotgun metagenomics, and metatranscriptomics) and metabolome (mass spectrometry and nuclear magnetic resonance spectroscopy) datasets were collected from six independent studies spanning multiple diseases. We used these datasets to evaluate two reference-based gene-to-metabolite prediction pipelines and a machine-learning (ML) based metabolic profile prediction approach. With the pre-trained model on over 900 microbiome-metabolome paired samples, the ML approach yielded the most accurate predictions (i.e., highest F1 scores) of metabolite occurrences in the human gut and outperformed reference-based pipelines in predicting differential metabolites between case and control subjects. Our findings demonstrate the possibility of predicting metabolites from microbiome sequencing data, while highlighting certain limitations in detecting differential metabolites, and provide a framework to evaluate metabolite prediction pipelines, which will ultimately facilitate future investigations on microbial metabolites and human health.

CoMA - an intuitive and user-friendly pipeline for amplicon-sequencing data analysis

In recent years, there has been a veritable boost in next-generation sequencing (NGS) of gene amplicons in biological and medical studies. Huge amounts of data are produced and need to be analyzed adequately. Various online and offline analysis tools are available; however, most of them require extensive expertise in computer science or bioinformatics, and often a Linux-based operating system. Here, we introduce "CoMA-Comparative Microbiome Analysis" as a free and intuitive analysis pipeline for amplicon-sequencing data, compatible with any common operating system. Moreover, the tool offers various useful services including data pre-processing, quality checking, clustering to operational taxonomic units (OTUs), taxonomic assignment, data post-processing, data visualization, and statistical appraisal. The workflow results in highly esthetic and publication-ready graphics, as well as output files in standardized formats (e.g. tab-delimited OTU-table, BIOM, NEWICK tree) that can be used for more sophisticated analyses. The CoMA output was validated by a benchmark test, using three mock communities with different sample characteristics (primer set, amplicon length, diversity). The performance was compared with that of Mothur, QIIME and QIIME2-DADA2, popular packages for NGS data analysis. Furthermore, the functionality of CoMA is demonstrated on a practical example, investigating microbial communities from three different soils (grassland, forest, swamp). All tools performed well in the benchmark test and were able to reveal the majority of all genera in the mock communities. Also for the soil samples, the results of CoMA were congruent to those of the other pipelines, in particular when looking at the key microbial players.

Genetic selection for growth drives differences in intestinal microbiota composition and parasite disease resistance in gilthead sea bream

BACKGROUND

The key effects of intestinal microbiota in animal health have led to an increasing interest in manipulating these bacterial populations to improve animal welfare. The aquaculture sector is no exception and in the last years, many studies have described these populations in different fish species. However, this is not an easy task, as intestinal microbiota is composed of very dynamic populations that are influenced by different factors, such as diet, environment, host age, and genetics. In the current study, we aimed to determine whether the genetic background of gilthead sea bream (Sparus aurata) influences the intestinal microbial composition, how these bacterial populations are modulated by dietary changes, and the effect of selection by growth on intestinal disease resistance. To that aim, three different groups of five families of gilthead sea bream that were selected during two generations for fast, intermediate, or slow growth (F3 generation) were kept together in the same open-flow tanks and fed a control or a well-balanced plant-based diet during 9 months. Six animals per family and dietary treatment were sacrificed and the adherent bacteria from the anterior intestinal portion were sequenced. In parallel, fish of the fast- and slow-growth groups were infected with the intestinal parasite Enteromyxum leei and the disease signs, prevalence, intensity, and parasite abundance were evaluated.

RESULTS

No differences were detected in alpha diversity indexes among families, and the core bacterial architecture was the prototypical composition of gilthead sea bream intestinal microbiota, indicating no dysbiosis in any of the groups. The plant-based diet significantly changed the microbiota in the intermediate- and slow-growth families, with a much lower effect on the fast-growth group. Interestingly, the smaller changes detected in the fast-growth families potentially accounted for more changes at the metabolic level when compared with the other families. Upon parasitic infection, the fast-growth group showed significantly lower disease signs and parasite intensity and abundance than the slow-growth animals.

CONCLUSIONS

These results show a clear genome-metagenome interaction indicating that the fast-growth families harbor a microbiota that is more flexible upon dietary changes. These animals also showed a better ability to cope with intestinal infections. Video Abstract.

PROVIT: Supplementary Probiotic Treatment and Vitamin B7 in Depression-A Randomized Controlled Trial

Gut microbiota are suspected to affect brain functions and behavior as well as lowering inflammation status. Therefore, an effect on depression has already been suggested by recent research. The aim of this randomized double-blind controlled trial was to evaluate the effect of probiotic treatment in depressed individuals. Within inpatient care, 82 currently depressed individuals were randomly assigned to either receive a multistrain probiotic plus biotin treatment or biotin plus placebo for 28 days. Clinical symptoms as well as gut microbiome were analyzed at the begin of the study, after one and after four weeks. After 16S rRNA analysis, microbiome samples were bioinformatically explored using QIIME, SPSS, R and Piphillin. Both groups improved significantly regarding psychiatric symptoms. Ruminococcus gauvreauii and Coprococcus 3 were more abundant and β-diversity was higher in the probiotics group after 28 days. KEGG-analysis showed elevated inflammation-regulatory and metabolic pathways in the intervention group. The elevated abundance of potentially beneficial bacteria after probiotic treatment allows speculations on the functionality of probiotic treatment in depressed individuals. Furthermore, the finding of upregulated vitamin B6 and B7 synthesis underlines the connection between the quality of diet, gut microbiota and mental health through the regulation of metabolic functions, anti-inflammatory and anti-apoptotic properties. Concluding, four-week probiotic plus biotin supplementation, in inpatient individuals with a major depressive disorder diagnosis, showed an overall beneficial effect of clinical treatment. However, probiotic intervention compared to placebo only differed in microbial diversity profile, not in clinical outcome measures.

Structure and predictive metabolic contribution of intestinal microbiota of Longfin yellowtail (Seriola rivoliana) juveniles in aquaculture systems

Seriola rivoliana intestinal microbiota (IM) was characterised under aquaculture conditions through 16S rRNA amplicon sequencing. Specimens of 30 days after hatching (DAH) were maintained in three tanks and fed under the same environmental conditions for characterisation 15 days prior to sampling. Three fish were randomly taken from each tank; total DNA extraction of the gut microbiota was performed to characterise microbial composition and its metabolic prediction. The V3 hypervariable region of the 16S rRNA was amplified and sequenced with Illumina pair-end technology. The prokaryotic components in the S. rivoliana intestine were dominated mainly by the phyla Proteobacteria, Firmicutes, Bacteroidetes, Cyanobacteria and Actinobacteria. No significant differences in beta diversity were detected in the three samples (tanks). However in alpha diversity, they were detected in juveniles of the same cohort within the same group, as exemplified by enrichment of certain bacterial groups, mainly of the Clostridia class, which were specific in each fish within the same tank. The metabolic prediction analyses suggested that S. rivoliana IM contribute to the metabolism of amino acids, carbohydrates, lipids, and immune system. This study provides the first IM characterisation under rearing conditions of S. rivoliana-a species with broad economic potential-and contributes to novel information for potential use of probiotics in future trials.

Multiomic features associated with mucosal healing and inflammation in paediatric Crohn's disease

BACKGROUND

The gastrointestinal microbiota has an important role in mucosal immune homoeostasis and may contribute to maintaining mucosal healing in Crohn's disease (CD).

AIM

To identify changes in the microbiota, metabolome and protease activity associated with mucosal healing in established paediatric CD METHODS: Twenty-five participants aged 3-18 years with CD, disease duration of over 6 months, and maintenance treatment with biological therapy were recruited. They were divided into a low calprotectin group (faecal calprotectin <100 μg/g, "mucosal healing," n = 11), and a high calprotectin group (faecal calprotectin >100 μg/g, "mucosal inflammation," n = 11). 16S gene-based metataxonomics, ¹ H-NMR spectroscopy-based metabolic profiling and protease activity assays were performed on stool samples.

RESULTS

Relative abundance of Dialister species was six-times greater in the low calprotectin group (q = 0.00999). Alpha and beta diversity, total protease activity and inferred metagenomic profiles did not differ between groups. Pentanoate (valerate) and lysine were principal discriminators in a machine-learning model which differentiated high and low calprotectin samples using NMR spectra (R² 0.87, Q² 0.41). Mean relative concentration of pentanoate was 1.35-times greater in the low calprotectin group (95% CI 1.03-1.68, P = 0.036) and was positively correlated with Dialister. Mean relative concentration of lysine was 1.54-times greater in the high calprotectin group (95% CI 1.05-2.03, P = 0.028).

CONCLUSIONS

This multiomic study identified an increase in Dialister species and pentanoate, and a decrease in lysine, in patients with "mucosal healing." It supports further investigation of these as potential novel therapeutic targets in CD.

PUMAA: A Platform for Accessible Microbiome Analysis in the Undergraduate Classroom

Improvements in high-throughput sequencing makes targeted amplicon analysis an ideal method for the study of human and environmental microbiomes by undergraduates. Multiple bioinformatics programs are available to process and interpret raw microbial diversity datasets, and the choice of programs to use in curricula is largely determined by student learning goals. Many of the most commonly used microbiome bioinformatics platforms offer end-to-end data processing and data analysis using a command line interface (CLI), but the downside for novice microbiome researchers is the steep learning curve often required. Alternatively, some sequencing providers include processing of raw data and taxonomy assignments as part of their pipelines. This, when coupled with available web-based or graphical user interface (GUI) analysis and visualization tools, eliminates the need for students or instructors to have extensive CLI experience. However, lack of universal data formats can make integration of these tools challenging. For example, tools for upstream and downstream analyses frequently use multiple different data formats which then require writing custom scripts or hours of manual work to make the files compatible. Here, we describe a microbial ecology bioinformatics curriculum that focuses on data analysis, visualization, and statistical reasoning by taking advantage of existing web-based and GUI tools. We created the Program for Unifying Microbiome Analysis Applications (PUMAA), which solves the problem of inconsistent files by formatting the output files from several raw data processing programs to seamlessly transition to a suite of GUI programs for analysis and visualization of microbiome taxonomic and inferred functional profiles. Additionally, we created a series of tutorials to accompany each of the microbiome analysis curricular modules. From pre- and post-course surveys, students in this curriculum self-reported conceptual and confidence gains in bioinformatics and data analysis skills. Students also demonstrated gains in biologically relevant statistical reasoning based on rubric-guided evaluations of open-ended survey questions and the Statistical Reasoning in Biology Concept Inventory. The PUMAA program and associated analysis tutorials enable students and researchers with no computational experience to effectively analyze real microbiome datasets to investigate real-world research questions.

Chytrid fungi shape bacterial communities on model particulate organic matter

Microbial colonization and degradation of particulate organic matter (POM) are important processes that influence the structure and function of aquatic ecosystems. Although POM is readily used by aquatic fungi and bacteria, there is a limited understanding of POM-associated interactions between these taxa, particularly for early-diverging fungal lineages. Using a model ecological system with the chitin-degrading freshwater chytrid fungus Rhizoclosmatium globosum and chitin microbeads, we assessed the impacts of chytrid fungi on POM-associated bacteria. We show that the presence of chytrids on POM alters concomitant bacterial community diversity and structure, including differing responses between chytrid life stages. We propose that chytrids can act as ecosystem facilitators through saprotrophic feeding by producing ‘public goods’ from POM degradation that modify bacterial POM communities. This study suggests that chytrid fungi have complex ecological roles in aquatic POM degradation not previously considered, including the regulation of bacterial colonization, community succession and subsequent biogeochemical potential.

Microbial ecology and predicted metabolic pathways in various oral environments from patients with acute endodontic infections

AIM

To assess in a cross-sectional clinical study the effect of antibiotics on the diversity, structure and metabolic pathways of bacterial communities in various oral environments in patients with acute primary infections.

METHODOLOGY

Samples of saliva (SA), supragingival biofilm (SB) and from the pulp cavity (PC) were collected from teeth with acute primary infections and then grouped according to previous use of antibiotics (NoAtb = no antibiotics [n = 6]; Atb = antibiotics [n = 6]). DNA sequencing was conducted using MiSeq (Illumina, San Diego, CA, USA). The V1-V3 hyper-variable region of the 16S rRNA gene was amplified. A custom Mothur pipeline was used for 16S rRNA processing. Subsequent analyses of the sequence dataset were performed in R (using vegan, phyloseq and ggplot2 packages) or QIIME.

RESULTS

Twelve patients aged from 22 to 56 years were recruited. Participants in the Atb group had taken the beta-lactamics amoxicillin (5/6) or cephalexin (1/6) for 2-3 days. A total of 332 bacterial taxa (OTUs) were identified, belonging to 120 genera, 60 families and nine phyla. Firmicutes (41%) and Bacteroidetes (38%) were the most abundant phyla in all samples. Taxa clustered significantly by oral site (PCoA analysis; P < 0.05, ANOSIM). Use of antibiotics had little effect on this clustering. However, SA, SB and PC had different degrees of richness, diversity and evenness. The greatest diversity was observed in SB samples and the least diversity was observed in PC samples. Metabolic prediction identified 163 pathways and previous use of antibiotics had a major effect on the estimated functional clustering in SA and PC samples.

CONCLUSION

The ecological niche had a strong influence on the bacterial content of samples from various oral sites. Previous exposure to antibiotics may exert an effect on the phylogenetic composition of SA. Metabolic pathways appear to be modulated by antimicrobial agents in SA and PC samples. The dynamics of host/microbial interactions in the apical region and the functional ecology of the infected pulp cavity should be revisited.

Sediment-associated bacterial community and predictive functionalities are influenced by choice of 16S ribosomal RNA hypervariable region(s): An amplicon-based diversity study

Meta-omics approaches such as high-throughput sequencing of 16S hypervariable region(s) [HVR(s)] is extensively applied for profiling microbial community. Several studies have deciphered the influence of HVR(s) on bacterial diversity; most of these were devoted to human body habitats. Extent to which targeted HVR(s) influences the diversity estimates of environmental samples is rather unclear. Here, we evaluated the performance of five widely used universal primer pairs spanning V1-V3, V3-V4, V4, V5-V6 and V7-V9 HVRs to characterize bacterial diversity and predictive functionality of complex marine sediments. Obtained results revealed that the HVR(s) V4 and V5-V6 represented the higher species richness than others while, V1-V3 and V7-V9 were unsuccessful to detect Bacteroidetes and Planctomycetes. Further, PICRUSt analysis showed that the selected HVR(s) also had significant impact on the predictive functional profile. Conclusively, this study proved that HVR selection has a profound effect on overall results and thus should be selected with utmost caution.

Maternal Diet and Infant Feeding Practices Are Associated with Variation in the Human Milk Microbiota at 3 Months Postpartum in a Cohort of Women with High Rates of Gestational Glucose Intolerance

BACKGROUND

Human milk contains a diverse community of bacteria believed to play a role in breast health and inoculation of the infant's gastrointestinal tract. The role of maternal nutrition and infant feeding practices on the human milk microbiota remains poorly understood.

OBJECTIVE

Our aim was to explore the associations between maternal diet (delivery to 3 mo postpartum), infant feeding practices, and the microbial composition and predicted function in milk from women with varied metabolic status.

METHODS

This was an exploratory analysis of a previously completed prospective cohort study of women with varying degrees of gestational glucose intolerance (NCT01405547). Milk samples (n = 93 mothers) were collected at 3 mo postpartum. Maternal dietary information (validated food-frequency questionnaire) and infant feeding practices (human milk exclusivity, frequency of direct breastfeeding per day) were collected. V4-16S ribosomal RNA gene sequencing (Illumina MiSeq) was conducted to determine microbiota composition.

RESULTS

Intake of polyunsaturated fat [β estimate (SE): 0.036 (0.018), P = 0.047] and fiber from grains [0.027 (0.013), P = 0.048] were positively associated with ɑ-diversity (Shannon index) of human milk. Overall microbial composition of human milk clustered based on human milk exclusivity (weighted UniFrac R2 = 0.034, P = 0.015; Bray-Curtis R2 = 0.041, P = 0.007), frequency of direct breastfeeding per day (Bray-Curtis R2 = 0.057, P = 0.026), and maternal fiber intake from grains (Bray-Curtis R2 = 0.055, P = 0.040). Total fiber, fiber from grains, dietary fat, and infant feeding practices were also associated with a number of differentially abundant taxa. The overall composition of predicted microbial functions was associated with total fiber consumption (Bray-Curtis R2 = 0.067, P = 0.036) and human milk exclusivity (Bray-Curtis R2 = 0.041, P = 0.013).

CONCLUSIONS

Maternal consumption of fiber and fat, as well as mother's infant feeding practices, are important determinants of the human milk microbiota. Understanding whether these microbial changes impact an infant's overall health and development requires future study.