A practical guide to amplicon and metagenomic analysis of microbiome data

Analysis and Interpretation of metagenomics data: an approach

Advances in next-generation sequencing technologies have accelerated the momentum of metagenomic studies, which is increasing yearly. The metagenomics field is one of the versatile applications in microbiology, where any interaction in the environment involving microorganisms can be the topic of study. Due to this versatility, the number of applications of this omics technology reached its horizons. Agriculture is a crucial sector involving crop plants and microorganisms interacting together. Hence, studying these interactions through the lenses of metagenomics would completely disclose a new meaning to crop health and development. The rhizosphere is an essential reservoir of the microbial community for agricultural soil. Hence, we focus on the R&D of metagenomic studies on the rhizosphere of crops such as rice, wheat, legumes, chickpea, and sorghum. These recent developments are impossible without the continuous advancement seen in the next-generation sequencing platforms; thus, a brief introduction and analysis of the available sequencing platforms are presented here to have a clear picture of the workflow. Concluding the topic is the discussion about different pipelines applied to analyze data produced by sequencing techniques and have a significant role in interpreting the outcome of a particular experiment. A plethora of different software and tools are incorporated in the automated pipelines or individually available to perform manual metagenomic analysis. Here we describe 8-10 advanced, efficient pipelines used for analysis that explain their respective workflows to simplify the whole analysis process.

Geographical patterns of Fejervarya limnocharis gut microbiota by latitude along mainland China’s coastline

The gut microbiota affects many aspects of host biology and plays key roles in the coevolutionary association with its host. Geographical gradients may play a certain role on gut microbiota variation in the natural environment. However, the distribution pattern of amphibian gut microbiota in the latitudinal gradient remains largely unexplored. Here, we sampled six natural populations of Fejervarya limnocharis along the eastern coastline of mainland China (spanning 20°–30° N = 1,300 km) using 16S rRNA amplicon sequencing to characterize the gut microbiota. First of all, a significant correlation between gut microbial diversity and latitude was observed in our research system. Second, we discovered that latitude influenced the composition of the gut microbiota of F. limnocharis. Finally, we detected that geographical distance could not determine gut microbiota composition in F. limnocharis. These results indicate that latitude can play an important role in shaping the gut microbial diversity of amphibian. Our study offers the first evidence that gut microbial diversity of amphibian presents a latitudinal pattern and highlights the need for increased numbers of individuals to be sampled during microbiome studies in wild populations along environmental gradients.

The microbial community associated with Parascaris spp. infecting juvenile horses

Background

Parasitic nematodes, including large roundworms colloquially known as ascarids, affect the health and well-being of livestock animals worldwide. The equine ascarids, Parascaris spp., are important parasites of juvenile horses and the first ascarids to develop widespread anthelmintic resistance. The microbiota has been shown to be an important factor in the fitness of many organisms, including parasitic nematodes, where endosymbiotic Wolbachia have been exploited for treatment of filariasis in humans.

Methods

This study used short-read 16S rRNA sequences and Illumina sequencing to characterize and compare microbiota of whole worm small intestinal stages and microbiota of male and female intestines and gonads. Diversity metrics including alpha and beta diversity, and the differential abundance analyses DESeq2, ANCOM-BC, corncob, and metagenomeSeq were used for comparisons.

Results

Alpha and beta diversity of whole worm microbiota did not differ significantly between groups, but Simpson alpha diversity was significantly different between female intestine (FI) and male gonad (MG) (P= 0.0018), and Shannon alpha diversity was significantly different between female and male gonads (P = 0.0130), FI and horse jejunum (HJ) (P = 0.0383), and FI and MG (P= 0.0001). Beta diversity (Fig. 2B) was significantly different between female and male gonads (P = 0.0006), male intestine (MI) and FG (P = 0.0093), and MG and FI (P = 0.0041). When comparing organs, Veillonella was differentially abundant for DESeq2 and ANCOM-BC (p < 0.0001), corncob (P = 0.0008), and metagenomeSeq (P = 0.0118), and Sarcina was differentially abundant across four methods (P < 0.0001). Finally, the microbiota of all individual Parascaris spp. specimens were compared to establish shared microbiota between groups.

Conclusions

Overall, this study provided important information regarding the Parascaris spp. microbiota and provides a first step towards determining whether the microbiota may be a viable target for future parasite control options.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05533-y.

Effects of different colored polyethylene mulching films on bacterial communities from soil during enrichment incubation

Studies have shown that mulching agricultural fields with plastic residues can influence microbial communities in the environment, but few studies have investigated the differences in the soil microbial communities in distinct areas under mulching with different colored plastic products. Thus, in this study, we explored how different colored polyethylene mulching films (PMFs) might affect soil bacterial communities during enrichment incubation. We found significant differences in the bacterial communities under different colored PMFs after incubation. Treatment with the same colored PMF obtained more similar bacterial community compositions. For instance, at the class level, Gammaproteobacteria and Bacteroidia were most abundant with black PMF, whereas Actinobacteria and Bacteroidia were most abundant with white PMF. The most abundant genera were Acinetobacter and Chryseobacterium with black PMF but Rhodanobacter and Paenarthrobacter with white PMF. Polyethylene- and hydrocarbon-degrading bacteria were the core members detected under both treatments, and the bacterial communities were predicted to have the potential for the biodegradation and metabolism of xenobiotics after enrichment culture according to the Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) tool. In addition, the bacterial communities in soil from Xinjiang treated with white PMF and in soil from Yangling treated with black PMF were strongly correlated and stable. Our results suggest that the color of the PMF applied affected the soil bacterial communities, where plastics with the same color may have recruited similar species of microorganisms, although the origins of these microorganisms were not the same.

Visualization and quantification of carbon "rusty sink" by rice root iron plaque: Mechanisms, functions, and global implications

Paddies contain 78% higher organic carbon (C) stocks than adjacent upland soils, and iron (Fe) plaque formation on rice roots is one of the mechanisms that traps C. The process sequence, extent and global relevance of this C stabilization mechanism under oxic/anoxic conditions remains unclear. We quantified and localized the contribution of Fe plaque to organic matter stabilization in a microoxic area (rice rhizosphere) and evaluated roles of this C trap for global C sequestration in paddy soils. Visualization and localization of pH by imaging with planar optodes, enzyme activities by zymography, and root exudation by ¹⁴ C imaging, as well as upscale modeling enabled linkage of three groups of rhizosphere processes that are responsible for C stabilization from the micro- (root) to the macro- (ecosystem) levels. The ¹⁴ C activity in soil (reflecting stabilization of rhizodeposits) with Fe²⁺ addition was 1.4-1.5 times higher than that in the control and phosphate addition soils. Perfect co-localization of the hotspots of β-glucosidase activity (by zymography) with root exudation (¹⁴ C) showed that labile C and high enzyme activities were localized within Fe plaques. Fe²⁺ addition to soil and its microbial oxidation to Fe³⁺ by radial oxygen release from rice roots increased Fe plaque (Fe³⁺ ) formation by 1.7-2.5 times. The C amounts trapped by Fe plaque increased by 1.1 times after Fe²⁺ addition. Therefore, Fe plaque formed from amorphous and complex Fe (oxyhydr)oxides on the root surface act as a "rusty sink" for organic matter. Considering the area of coverage of paddy soils globally, upscaling by model revealed the radial oxygen loss from roots and bacterial Fe oxidation may trap up to 130 Mg C in Fe plaques per rice season. This represents an important annual surplus of new and stable C to the existing C pool under long-term rice cropping.

The response of bacterial community to UVB was significantly different between immature periphyton and mature periphyton, but not for physiological indicators

Periphyton has important ecological functions. It can even exist in environments with strong ultraviolet radiation. However, knowledge of periphyton under ultraviolet is limited, which limits the understanding and application of periphyton in environments with high ultraviolet radiation. In this study, immature periphyton (IMP) and mature periphyton (MP) under ultraviolet B (UVB) irradiation were investigated and compared in terms of physiological characteristics and bacterial community. Analysis of the physiological characteristics showed that the response patterns of IMP and MP to UVB were similar. IMP and MP could adapt to UVB of 1 W/m² well. However, high-intensity UVB (2 and 3 W/m²) reduced the periphyton biomass, inhibited photosynthesis and antioxidant enzyme activity and caused severe lipid peroxidation in both IMP and MP. Integrated Biological Response (IBR) analysis and toxicological model fitting showed that the ED₅₀ values of UVB for IMP and MP were 1.25 and 1.50 W/m², respectively. 16 S rRNA gene analysis showed that in both IMP and MP, bacterial community composition, assembly and function were affected by UVB. In addition, the response of the bacterial community in IMP to UVB was stronger than that in MP. The diversity of the IMP community was inhibited by UVB, but that of the MP community was not. Proteobacteria and Deinococcus-Thermus are key microorganisms responsible for tolerance to UVB stress. Neutral community model fitting showed that both UVB and the development process caused the determinism of bacterial succession. However, UVB may weaken the deterministic process caused by development. Functional prediction showed that many metabolic functions of periphyton were inhibited by UVB in IMP and MP. However, UVB caused different changes (enhancement or inhibition) of some ecological functions in them. This study provides valuable information for understanding periphyton in environments with UVB radiation, which may be used to improve the application of periphyton in these environments.

Understanding of mercury and methylmercury transformation in sludge composting by metagenomic analysis

Municipal sewage especially the produced sewage sludge is a significant source releasing mercury (Hg) to the environment. However, the Hg speciation especially methylmercury (MeHg) transformation in sewage sludge treatment process remains poorly understood. This study investigated the transformation of Hg speciation especially MeHg in sludge composting. The distribution of Hg transformation related gene pairs hgcAB and merAB, and their putative microbial hosts were comprehensively analyzed. Both Hg (from 3.16±0.22 mg/kg to 3.20±0.19 mg/kg) and MeHg content (from 4.77±0.64 ng/g to 4.36±0.37 ng/g) were not obviously changed before and after composting, but about 19.69% of Hg and 27.36% of MeHg were lost according to mass balance calculation. The metagenomic analysis further revealed that anaerobes (Desulfobacterota and Euryarchaeota) were the mainly putative Hg methylators especially carrying high abundance of hgcA gene in the initial periods of composting. Among the 151 reconstructed metagenome-assembled genomes (MAGs), only 4 hgcA gene carriers (Myxococcota, Firmicutes, Cyclobacteriaceae, and Methanothermobacter) and 16 merB gene carriers were identified. But almost all of the MAGs carried hgcB gene and merA gene. The merA gene was widely distributed in genomes, which indicated the widespread functionality of microbes for reducing Hg(II) to Hg(0). The hgcA carrying microbes tends to present the similar metabolic pathways including methanogenesis and sulfur metabolism. Besides, both the irregular distribution of hgcA in various species (including Actinobacteria, Archaea, Bacteroidetes, Desulfobacterota, Euryarchaeota, and Nitrospirae, etc.) and opposite evolution trends between hgcA gene abundance and its host genome abundance can be an indication of horizontal gene transfer or gene deletions of hgcA during composting. Our findings thus revealed that sludge composting is not only a hotspot for Hg speciation transformation, but also a potential hotspot for MeHg transformation.

CRAMdb: a comprehensive database for composition and roles of microbiome in animals

CRAMdb (a database for composition and roles of animal microbiome) is a comprehensive resource of curated and consistently annotated metagenomes for non-human animals. It focuses on the composition and roles of the microbiome in various animal species. The main goal of the CRAMdb is to facilitate the reuse of animal metagenomic data, and enable cross-host and cross-phenotype comparisons. To this end, we consistently annotated microbiomes (including 16S, 18S, ITS and metagenomics sequencing data) of 516 animals from 475 projects spanning 43 phenotype pairs to construct the database that is equipped with 9430 bacteria, 278 archaea, 2216 fungi and 458 viruses. CRAMdb provides two main contents: microbiome composition data, illustrating the landscape of the microbiota (bacteria, archaea, fungi, and viruses) in various animal species, and microbiome association data, revealing the relationships between the microbiota and various phenotypes across different animal species. More importantly, users can quickly compare the composition of the microbiota of interest cross-host or body site and the associated taxa that differ between phenotype pairs cross-host or cross-phenotype. CRAMdb is freely available at (http://www.ehbio.com/CRAMdb).

The Role of the Cutaneous Mycobiome in Atopic Dermatitis

Atopic dermatitis is a chronic inflammatory skin disorder characterized by eczematous lesions, itch, and a significant deterioration in the quality of life. Recently, microbiome dysbiosis has been implicated in the pathogenesis of atopic dermatitis. Changes in the fungal microbiome (also termed mycobiome) appear to be an important factor influencing the clinical picture of this entity. This review summarizes the available insights into the role of the cutaneous mycobiome in atopic dermatitis and the new research possibilities in this field. The prevalence and characteristics of key fungal species, the most important pathogenesis pathways, as well as classic and emerging therapies of fungal dysbiosis and infections complicating atopic dermatitis, are presented.

Rumen Microbiome Reveals the Differential Response of CO2 and CH4 Emissions of Yaks to Feeding Regimes on the Qinghai-Tibet Plateau

Shifts in feeding regimes are important factors affecting greenhouse gas (GHG) emissions from livestock farming. However, the quantitative values and associated drivers of GHG emissions from yaks (Bos grunniens) following shifts in feeding regimes have yet to be fully described. In this study, we aimed to investigate CH4 and CO2 emissions differences of yaks under different feeding regimes and their potential microbial mechanisms. Using static breathing chamber and Picarro G2508 gas concentration analyzer, we measured the CO2 and CH4 emissions from yaks under traditional grazing (TG) and warm-grazing and cold-indoor feeding (WGCF) regimes. Microbial inventories from the ruminal fluid of the yaks were determined via Illumina 16S rRNA and ITS sequencing. Results showed that implementing the TG regime in yaks decreased their CO2 and CH4 emissions compared to the WGCF regime. The alpha diversity of ruminal archaeal community was higher in the TG regime than in the WGCF regime. The beta diversity showed that significant differences in the rumen microbial composition of the TG regime and the WGCF regime. Changes in the rumen microbiota of the yaks were driven by differences in dietary nutritional parameters. The relative abundances of the phyla Neocallimastigomycota and Euryarchaeota and the functional genera Prevotella, Ruminococcus, Orpinomyces, and Methanobrevibacter were significantly higher in the WGCF regime than in the TG regime. CO2 and CH4 emissions from yaks differed mainly because of the enrichment relationship of functional H2- and CO2-producing microorganisms, hydrogen-consuming microbiota, and hydrogenotrophic methanogenic microbiota. Our results provided a view that it is ecologically important to develop GHG emissions reduction strategies for yaks on the Qinghai-Tibet Plateau based on traditional grazing regime.

Pesticide soil microbial toxicity: setting the scene for a new pesticide risk assessment for soil microorganisms (IUPAC Technical Report)

Pesticides constitute an integral part of modern agriculture. However, there are still concerns about their effects on non-target organisms. To address this the European Commission has imposed a stringent regulatory scheme for new pesticide compounds. Assessment of the aquatic toxicity of pesticides is based on a range of advanced tests. This does not apply to terrestrial ecosystems, where the toxicity of pesticides on soil microorganisms, is based on an outdated and crude test (N mineralization). This regulatory gap is reinforced by the recent methodological and standardization advances in soil microbial ecology. The inclusion of such standardized tools in a revised risk assessment scheme will enable the accurate estimation of the toxicity of pesticides on soil microorganisms and on associated ecosystem services. In this review we (i) summarize recent work in the assessment of the soil microbial toxicity of pesticides and point to ammonia-oxidizing microorganisms (AOM) and arbuscular mycorrhizal fungi (AMF) as most relevant bioindicator groups (ii) identify limitations in the experimental approaches used and propose mitigation solutions, (iii) identify scientific gaps and (iv) propose a new risk assessment procedure to assess the effects of pesticides on soil microorganisms.

Microbial communities associated with thermogenic gas hydrate-bearing marine sediments in Qiongdongnan Basin, South China Sea

Biogenic and thermogenic gas are two major contributors to gas hydrate formation. Methane hydrates from both origins may have critical impacts on the ecological properties of marine sediments. However, research on microbial diversity in thermogenic hydrate-containing sediments is limited. This study examined the prokaryotic diversity and distributions along a sediment core with a vertical distribution of thermogenic gas hydrates with different occurrences obtained from the Qiongdongnan Basin by Illumina sequencing of 16S rRNA genes as well as molecular and geochemical techniques. Here, we show that gas hydrate occurrence has substantial impacts on both microbial diversity and community composition. Compared to the hydrate-free zone, distinct microbiomes with significantly higher abundance and lower diversity were observed within the gas hydrate-containing layers. Gammaproteobacteria and Actinobacterota dominated the bacterial taxa in all collected samples, while archaeal communities shifted sharply along the vertical profile of sediment layers. A notable stratified distribution of anaerobic methanotrophs shaped by both geophysical and geochemical parameters was also determined. In addition, the hydrate-free zone hosted a large number of rare taxa that might perform a fermentative breakdown of proteins in the deep biosphere and probably respond to the hydrate formation.

The Application of Metagenomics to Study Microbial Communities and Develop Desirable Traits in Fermented Foods

The microbial communities present within fermented foods are diverse and dynamic, producing a variety of metabolites responsible for the fermentation processes, imparting characteristic organoleptic qualities and health-promoting traits, and maintaining microbiological safety of fermented foods. In this context, it is crucial to study these microbial communities to characterise fermented foods and the production processes involved. High Throughput Sequencing (HTS)-based methods such as metagenomics enable microbial community studies through amplicon and shotgun sequencing approaches. As the field constantly develops, sequencing technologies are becoming more accessible, affordable and accurate with a further shift from short read to long read sequencing being observed. Metagenomics is enjoying wide-spread application in fermented food studies and in recent years is also being employed in concert with synthetic biology techniques to help tackle problems with the large amounts of waste generated in the food sector. This review presents an introduction to current sequencing technologies and the benefits of their application in fermented foods.

Detection of Acquired Antibiotic Resistance Genes in Domestic Pig (Sus scrofa) and Common Carp (Cyprinus carpio) Intestinal Samples by Metagenomics Analyses in Hungary

The aim of this study was metagenomics analyses of acquired antibiotic-resistance genes (ARGs) in the intestinal microbiome of two important food-animal species in Hungary from a One Health perspective. Intestinal content samples were collected from 12 domestic pigs (Sus scrofa) and from a common carp (Cyprinus carpio). Shotgun metagenomic sequencing of DNA purified from the intestinal samples was performed on the Illumina platform. The ResFinder database was applied for detecting acquired ARGs in the assembled metagenomic contigs. Altogether, 59 acquired ARG types were identified, 51 genes from domestic pig and 12 genes from the carp intestinal microbiome. ARG types belonged to the antibiotic classes aminoglycosides (27.1%), tetracyclines (25.4%), β-lactams (16.9%), and others. Of the identified ARGs, tet(E), a bla_OXA-_48-like β-lactamase gene, as well as cphA4, ampS, aadA2, qnrS2, and sul1, were identified only in carp but not in swine samples. Several of the detected acquired ARGs have not yet been described from food animals in Hungary. The tet(Q), tet(W), tet(O), and mef(A) genes detected in the intestinal microbiome of domestic pigs had also been identified from free-living wild boars in Hungary, suggesting a possible relationship between the occurrence of acquired ARGs in domestic and wild animal populations.

Improving analysis of the vaginal microbiota of women undergoing assisted reproduction using nanopore sequencing

PURPOSE

Subclinical alterations of the vaginal microbiome have been described to be associated with female infertility and may serve as predictors for failure of in vitro fertilization treatment. While large prospective studies to delineate the role of microbial composition are warranted, integrating microbiome information into clinical management depends on economical and practical feasibility, specifically on a short duration from sampling to final results. The currently most used method for microbiota analysis is either metagenomics sequencing or amplicon-based microbiota analysis using second-generation methods such as sequencing-by-synthesis approaches (Illumina), which is both expensive and time-consuming. Thus, additional approaches are warranted to accelerate the usability of the microbiome as a marker in clinical praxis.

METHODS

Herein, we used a set of ten selected vaginal swabs from women undergoing assisted reproduction, comparing and performing critical optimization of nanopore-based microbiota analysis with the results from MiSeq-based data as a quality reference.

RESULTS

The analyzed samples carried varying community compositions, as shown by amplicon-based analysis of the V3V4 region of the bacterial 16S rRNA gene by MiSeq sequencing. Using a stepwise procedure to optimize adaptation, we show that a close approximation of the microbial composition can be achieved within a reduced time frame and at a minimum of costs using nanopore sequencing.

CONCLUSIONS

Our work highlights the potential of a nanopore-based methodical setup to support the feasibility of interventional studies and contribute to the development of microbiome-based clinical decision-making in assisted reproduction.

Weizmannia coagulans BC2000 Plus Ellagic Acid Inhibits High-Fat-Induced Insulin Resistance by Remodeling the Gut Microbiota and Activating the Hepatic Autophagy Pathway in Mice

(1) Background: Ellagic acid (EA) acts as a product of gut microbiota transformation to prevent insulin resistance, which is limited by high-fat diet (HFD)-induced dysbiosis. The aim of this study was to investigate the synergistic effects and mechanisms of supplementation with the probiotic Weizmannia coagulans (W. coagulans) on the prevention of insulin resistance by EA; (2) Methods: C57BL/6J mice were divided into five groups (n = 10/group): low-fat-diet group, high-fat-diet group, EA intervention group, EA + W. coagulans BC77 group, and EA + W. coagulans BC2000 group; (3) Result: W. coagulans BC2000 showed a synergistic effect on EA's lowering insulin resistance index and inhibiting high-fat diet-induced endotoxemia. The combined effect of BC2000 and EA activated the autophagy pathway in the mouse liver, a urolithin-like effect. This was associated with altered β-diversity of gut microbiota and increased Eggerthellaceae, a potential EA-converting family. Ellagic acid treatment alone and the combined use of ellagic acid and W. coagulans BC77 failed to activate the hepatic autophagy pathway; (4) Conclusions: W. coagulans BC2000 can assist EA in its role of preventing insulin resistance. This study provides a basis for the development of EA-rich functional food supplemented with W. coagulans BC2000.

Quantification and characterization of mouse and human tissue-resident microbiota by qPCR and 16S sequencing

The tissue-resident microbiota is an integral component of multiple tumor types, but it remains challenging to characterize its abundance and composition due to its low biomass. Here, we describe an optimized protocol for quantification and profiling of tissue-resident microbiota. The major optimized steps include DNA extraction, qPCR, 16S library construction, and bioinformatics analysis. This protocol enables robust and accurate characterization of the dynamics of normal and tumor tissue-resident microbiota at its physiological abundance from both mouse and human origins. For complete details on the use and execution of this protocol, please refer to Fu et al. (2022).

Machine learning and deep learning applications in microbiome research

The many microbial communities around us form interactive and dynamic ecosystems called microbiomes. Though concealed from the naked eye, microbiomes govern and influence macroscopic systems including human health, plant resilience, and biogeochemical cycling. Such feats have attracted interest from the scientific community, which has recently turned to machine learning and deep learning methods to interrogate the microbiome and elucidate the relationships between its composition and function. Here, we provide an overview of how the latest microbiome studies harness the inductive prowess of artificial intelligence methods. We start by highlighting that microbiome data - being compositional, sparse, and high-dimensional - necessitates special treatment. We then introduce traditional and novel methods and discuss their strengths and applications. Finally, we discuss the outlook of machine and deep learning pipelines, focusing on bottlenecks and considerations to address them.

Comparative study on the microbiota of colostrum and nipple skin from lactating mothers separated from their newborn at birth in China

Increasing studies have found breast milk (BM) contains its own microbiota. However, the route through which microbes enter the BM is still unclear. In order to verify the entero-mammary pathway of BM, we designed a rigorous study that prevented oral bacteria from contaminating the breast and nipple skin (NS) during baby nursing. Thirty-one healthy, postpartum mothers living in southern China who were immediately separated from their newborn after delivery were enrolled in this study. Using an aseptic protocol for sampling, sterile water was used to wash the NS and was then collected. Then the first drop of BM was discarded and colostrum was collected manually. Amplicon sequencing was performed targeting the V3–V4 region of the bacterial 16S rRNA gene, and the differences between the microbiota of the colostrum and NS were analyzed. Additionally, the effects of environmental factors, such as the delivery mode and intrapartum antibiotic exposure, on the diversity of the colostrum microbiota were also analyzed. We found significant differences in the α diversity and richness between the BM and NS as evidenced by richness, Chao1, and Simpson indices. There were 170 operational taxonomic units (OTUs) shared by colostrum and NS, while 111 and 87 OTUs were unique, respectively, as well as a clear distinction in OTUs was observed by unifrac binary analysis between them. Linear discriminant analysis effect size analysis found that anaerobes, such as Bifidobacterium and Pantoea at the genus level and enterobacteria including Enterobacteriaceae at the family level, were predominant in the colostrum, while the predominant bacteria on the NS were Bacteroides, Staphylococcus, and Parabacteroides at the genus level. BM is colonized by bacteria prior to baby suckling, and the diversity of the colostrum microbiota differs from that of the NS. The predominant microbiota taxa in BM indicated that they were likely to be transferred to the breast through the intestinal tract. Our study provides direct evidence for the revolutionary active migration hypothesis. Additionally, factors like intrapartum antibiotic exposure did not significantly affect the diversity of the microbiota in the BM. Therefore, it is suggested that mothers continue to provide BM for their newborns during separation.

Fluoride induced leaky gut and bloom of Erysipelatoclostridium ramosum mediate the exacerbation of obesity in high-fat-diet fed mice

INTRODUCTION

Fluoride is widely presented in drinking water and foods. A strong relation between fluoride exposure and obesity has been reported. However, the potential mechanisms on fluoride-induced obesity remain unexplored. Objectives and methods The effects of fluoride on the obesity were investigated using mice model. Furthermore, the role of gut homeostasis in exacerbation of the obesity induced by fluoride was evaluated. Results The results showed that fluoride alone did not induce obesity in normal diet (ND) fed mice, whereas, it could trigger exacerbation of obesity in high-fat diet (HFD) fed mice. Fluoride impaired intestinal barrier and activated Toll-like receptor 4 (TLR4) signaling to induce obesity, which was further verified in TLR4^-/- mice. Furthermore, fluoride could deteriorate the gut microbiota in HFD mice. The fecal microbiota transplantation from fluoride-induced mice was sufficient to induce obesity, while the exacerbation of obesity by fluoride was blocked upon gut microbiota depletion. The fluoride-induced bloom of Erysipelatoclostridium ramosum was responsible for exacerbation of obesity. In addition, a potential strategy for prevention of fluoride-induced obesity was proposed by intervention with polysaccharides from Fuzhuan brick tea. Conclusion Overall, these results provide the first evidence of a comprehensive cross-talk mechanism between fluoride and obesity in HFD fed mice, which is mediated by gut microbiota and intestinal barrier. E. ramosum was identified as a crucial mediator of fluoride induced obesity, which could be explored as potential target for prevention and treatment of obesity with exciting translational value.

Metagenomic evidence of suppressed methanogenic pathways along soil profile after wetland conversion to cropland

Wetland conversion to cropland substantially suppresses methane (CH4) emissions due to the strong suppression of methanogenesis, which consists of various pathways. In this study, we evaluated the cultivation impacts on four predominant CH4 production pathways, including acetate, carbon dioxide (CO2), methylamines, and methanol, in a wetland and cultivated cropland in northeastern China. The results showed significant suppression of CH4 production potential and the abundance of genes for all four methanogenic pathways in cropland. The consistency between CH4 production and methanogenesis genes indicates the robustness of genomic genes in analyzing methanogenesis. The suppression effects varied across seasons and along soil profiles, most evident in spring and 0 to 30 cm layers. The acetate pathway accounted for 55% in wetland vs. 70% in the cropland of all functional genes for CH4 production; while the other three pathways were stronger in response to cultivation, which presented as stronger suppressions in both abundance of functional genes (declines are 52% of CO2 pathway, 68% of methanol pathway, and 62% of methylamines pathway, vs. 19% of acetate pathway) and their percentages in four pathways (from 20 to 15% for CO2, 15 to 9% for methylamines, and 10 to 6% for methanol pathway vs. 55 to 70% for acetate pathway). The structural equation models showed that substrate availability was most correlated with CH4 production potential in the wetland, while the positive correlations of acetate, CO2, and methylamine pathways with CH4 production potential were significant in the cropland. The quantitative responses of four CH4 production pathways to land conversion reported in this study provide benchmark information for validating the CH4 model in simulating CH4 cycling under land use and land cover change.

ggClusterNet: An R package for microbiome network analysis and modularity-based multiple network layouts

The network analysis has attracted increasing attention and interest from ecological academics, thus it is of great necessity to develop more convenient and powerful tools. For that reason, we have developed an R package, named "ggClusterNet," to complete and display the network analysis in an easier manner. In that package, ten network layout algorithms are designed to better display the modules of microbiome network (randomClusterG, PolygonClusterG, PolygonRrClusterG, ArtifCluster, randSNEClusterG, PolygonModsquareG, PolyRdmNotdCirG, model_Gephi.2, model_igraph, and model_maptree). For the convenience of the users, many functions related to microbial network analysis, such as corMicor(), net_properties(), node_properties(), ZiPiPlot(), random_Net_compate(), are integrated to complete the network mining. Furthermore, the pipeline function named network.2() and corBionetwork() are also added for the quick achievement of the network or bipartite network analysis as well as their in-depth mining. The ggClusterNet is publicly available via GitHub (https://github.com/taowenmicro/ggClusterNet/) or Gitee (https://gitee.com/wentaomicro/ggClusterNet) for users' access. A complete description of the usages can be found on the manuscript's GitHub page (https://github.com/taowenmicro/ggClusterNet/wiki).

Microbial community roles and chemical mechanisms in the parasitic development of Orobanche cumana

Orobanche cumana Wallr. is a holoparasite weed that extracts water and nutrients from its host the sunflower, thereby causing yield reductions and quality losses. However, the number of O. cumana parasites in the same farmland is distinctly different. The roots of some hosts have been heavily parasitized, while others have not been parasitized. What are the factors contributing to this phenomenon? Is it possible that sunflower interroot microorganisms are playing a regulatory role in this phenomenon? The role of the microbial community in this remains unclear. In this study, we investigated the rhizosphere soil microbiome for sunflowers with different degrees of O. cumana parasitism, that is, healthy, light infection, moderate infection, and severe infection on the sunflower roots. The microbial structures differed significantly according to the degree of parasitism, where Xanthomonadaceae was enriched in severe infections. Metagenomic analyses revealed that amino acid, carbohydrate, energy, and lipid metabolism were increased in the rhizosphere soils of severely infected sunflowers, which were attributed to the proliferation of Lysobacter. Lysobacter antibioticus (HX79) was isolated and its capacity to promote O. cumana seed germination and increase the germ tube length was confirmed by germination and pot experiments. Cyclo(Pro-Val), an active metabolite of strain HX79, was identified and metabolomic and molecular docking approaches confirmed it was responsible for promoting O. cumana seed germination and growth. And we found that Pseudomonas mandelii HX1 inhibited the growth of O. cumana in the host rhizosphere soil. Our findings clarify the role of rhizosphere microbiota in regulating the parasite O. cumana to possibly facilitate the development of a new weed suppression strategy.

Towards standardized and reproducible research in skin microbiomes

Skin is a complex organ serving a critical role as a barrier and mediator of interactions between the human body and its environment. Recent studies have uncovered how resident microbial communities play a significant role in maintaining the normal healthy function of the skin and the immune system. In turn, numerous host-associated and environmental factors influence these communities' composition and diversity across the cutaneous surface. In addition, specific compositional changes in skin microbiota have also been connected to the development of several chronic diseases. The current era of microbiome research is characterized by its reliance on large data sets of nucleotide sequences produced with high-throughput sequencing of sample-extracted DNA. These approaches have yielded new insights into many previously uncharacterized microbial communities. Application of standardized practices in the study of skin microbial communities could help us understand their complex structures, functional capacities, and health associations and increase the reproducibility of the research. Here, we overview the current research in human skin microbiomes and outline challenges specific to their study. Furthermore, we provide perspectives on recent advances in methods, analytical tools and applications of skin microbiomes in medicine and forensics.

Sangerbox: A comprehensive, interaction-friendly clinical bioinformatics analysis platform

In recent decades, with the continuous development of high-throughput sequencing technology, data volume in medical research has increased, at the same time, almost all clinical researchers have their own independent omics data, which provided a better condition for data mining and a deeper understanding of gene functions. However, for these large amounts of data, many common and cutting-edge effective bioinformatics research methods still cannot be widely used. This has encouraged the establishment of many analytical platforms, a portion of databases or platforms were designed to solve the special analysis needs of users, for instance, MG RAST, IMG/M, Qiita, BIGSdb, and TRAPR were developed for specific omics research, and some databases or servers provide solutions for special problems solutions. Metascape was designed to only provide functional annotations of genes as well as function enrichment analysis; BioNumerics and RidomSeqSphere+ perform multilocus sequence typing; CARD provides only antimicrobial resistance annotations. Additionally, some web services are outdated, and inefficient interaction often fails to meet the needs of researchers, such as our previous versions of the platform. Therefore, the demand to complete massive data processing tasks urgently requires a comprehensive bioinformatics analysis platform. Hence, we have developed a website platform, Sangerbox 3.0 (http://vip.sangerbox.com/), a web-based tool platform. On a user-friendly interface that also supports differential analysis, the platform provides interactive customizable analysis tools, including various kinds of correlation analyses, pathway enrichment analysis, weighted correlation network analysis, and other common tools and functions, users only need to upload their own corresponding data into Sangerbox 3.0, select required parameters, submit, and wait for the results after the task has been completed. We have also established a new interactive plotting system that allows users to adjust the parameters in the image; moreover, optimized plotting performance enables users to adjust large-capacity vector maps on the web site. At the same time, we have integrated GEO, TCGA, ICGC, and other databases and processed data in batches, greatly reducing the difficulty to obtain data and improving the efficiency of bioimformatics study for users. Finally, we also provide users with rich sources of bioinformatics analysis courses, offering a platform for researchers to share and exchange knowledge.

Integrating Multiple Omics Identifies Phaeoacremonium rubrigenum Acting as Aquilaria sinensis Marker Fungus to Promote Agarwood Sesquiterpene Accumulation by Inducing Plant Host Phosphorylation

The present study aimed to explore the factors that promote persistent agarwood accumulation. To this end, we first investigated the morphological changes and volatile compound distribution in five layers of "Guan Xiang" agarwood. The agarwood-normal transition layer (TL), an essential layer of persistent agarwood accumulation, showed clear metabolic differences by microscopy and GC-MS analysis. Microbiome analysis revealed that Phaeocremonium rubrigenum was the predominant biomarker fungus in the TL of "Guan Xiang" agarwood samples. Among the seven isolated fungi, P. rubrigenum exhibited a significantly heightened ability to induce the production in Aquilaria sinensis seedlings, especially for sesquiterpene. Tracing the proteome profile changes in P. rubrigenum-induced A. sinensis calli for 18 ds showed that the fungus-induced sesquiterpene biosynthesis increased mainly through the mevalonate (MVA) pathway. Specifically, the phosphorylation modification level, instead of the protein abundance of transcription factors (TFs), showed corresponding changes during sesquiterpene biosynthesis, thus indicating that induced phosphorylation is the key reason for enhanced sesquiterpene production. IMPORTANCE Agarwood is an expensive resinous portion derived from Aquilaria plants and has been widely used as medicine, incense, and perfume. The factors involved in steady agarwood accumulation remain elusive. Our current study suggests that as a TL marker fungus, P. rubrigenum could persistently promote agarwood sesquiterpene accumulation by inducing phosphorylation of the TFs-MVA network in A. sinensis. Moreover, our work provides strategies to improve agarwood industry management and sheds light on the potential molecular mechanisms of plant adaptation to native microbial conditions.

Convergent and Divergent Age Patterning of Gut Microbiota Diversity in Humans and Nonhuman Primates

The gut microbiome has significant effects on healthy aging and aging-related diseases, whether in humans or nonhuman primates. However, little is known about the divergence and convergence of gut microbial diversity between humans and nonhuman primates during aging, which limits their applicability for studying the gut microbiome's role in human health and aging. Here, we performed 16S rRNA gene sequencing analysis for captive rhesus macaques (Macaca mulatta) and compared this data set with other freely available gut microbial data sets containing four human populations (Chinese, Japanese, Italian, and British) and two nonhuman primates (wild lemurs [Lemur catta] and wild chimpanzees [Pan troglodytes]). Based on the consistent V4 region of the 16S rRNA gene, beta diversity analysis suggested significantly separated gut microbial communities associated with host backgrounds of seven host groups, but within each group, significant gut microbial divergences were observed, and indicator bacterial genera were identified as associated with aging. We further discovered six common anti-inflammatory gut bacteria (Prevotellamassilia, Prevotella, Gemmiger, Coprococcus, Faecalibacterium, and Roseburia) that had butyrate-producing potentials suggested by pangenomic analysis and that showed similar dynamic changes in at least two selected host groups during aging, independent of distinct host backgrounds. Finally, we found striking age-related changes in 66 plasma metabolites in macaques. Two highly changed metabolites, hydroxyproline and leucine, enriched in adult macaques were significantly and positively correlated with Prevotella and Prevotellamassilia. Furthermore, genus-level pangenome analysis suggested that those six common indicator bacteria can synthesize leucine and arginine as hydroxyproline and proline precursors in both humans and macaques. IMPORTANCE This study provides the first comprehensive investigation of age patterning of gut microbiota of four human populations and three nonhuman primates and found that Prevotellamassilia, Prevotella, Gemmiger, Coprococcus, Faecalibacterium, and Roseburia may be common antiaging microbial markers in both humans and nonhuman primates due to their potential metabolic capabilities for host health benefits. Our results also provide key support for using macaques as animal models in studies of the gut microbiome's role during human aging.

Gut microbiota on admission as predictive biomarker for acute necrotizing pancreatitis

Background

Acute necrotizing pancreatitis (NP), a severe form of acute pancreatitis (AP), has higher mortality and worse outcome than non-necrotizing pancreatitis (non-NP). Infected NP is a devastating subgroup of NP. To date neither NP nor infected NP has robust prediction strategies, which may delay early recognition and timely intervention. Recent studies revealed correlations between disturbed gut microbiota and AP severity. Some features of intestinal microbiota have the potential to become biomarkers for NP prediction.

Methods

We performed 16S rRNA sequencing to analyze gut microbiota features in 20 healthy controls (HC), and 58 AP patients on hospital admission. The AP patients were later classified into NP and non-NP groups based on subsequent diagnostic imaging features. Random forest regression model and ROC curve were applied for NP and infected NP prediction. PIRCUSt2 was used for bacterial functional pathway prediction analysis.

Results

We found that the three groups (HC, NP, and non-NP) had distinct microorganism composition. NP patients had reduced microbial diversity, higher abundance of Enterobacteriales, but lower abundance of Clostridiales and Bacteroidales compared with the non-NP group. Correlation analyses displayed that intestine bacterial taxonomic alterations were related to severity, ICU admission, and prognosis. By pathway prediction, species more abundant in NP patients had positive correlation with synthesis and degradation of ketone bodies, and benzoate degradation. Enterococcus faecium (ASV2) performed best in discriminating NP and non-NP patients. Finegoldia magna (ASV3) showed the maximal prediction capacity among all ASVs and had comparable accuracy with Balthazar CT to detect patients with infected NP.

Conclusions

Our study suggests that NP patients have distinct intestinal microbiota on admission compared to non-NP patients. Dysbiosis of intestinal microbiota might influence NP progression through ketone body or benzoate metabolism. Enterococcus faecium and Finegoldia magna are potential predictors for NP and infected NP. Our findings explore biomarkers which may inform clinical decision-making in AP and shed light on further studies on NP pathophysiology and management.

Application of Microbiome in Forensics

Recent advances in next-generation sequencing technology and improvements in bioinformatics have expanded the scope of microbiome analysis as a forensic tool. Microbiome research is concerned with the study of the compositional profile and diversity of microbial flora as well as the interactions between microbes, hosts, and the environment. It has opened up many new possibilities for forensic analysis. In this review, we discuss various applications of microbiomes in forensics, including identification of individuals, geolocation inference, post-mortem interval (PMI) estimation, and others.

Interaction between Plasma Metabolomics and Intestinal Microbiome in db/db Mouse, an Animal Model for Study of Type 2 Diabetes and Diabetic Kidney Disease

Evidence has demonstrated that either metabolites or intestinal microbiota are involved in the pathogenesis of type 2 diabetes (T2D) and diabetic kidney disease (DKD). To explore the interaction between plasma metabolomics and intestinal microbiome in the progress of T2D-DKD, in the current study, we analyzed metabolomics in the plasma of db/db mice with liquid chromatography–mass spectrometry and also examined intestinal prokaryotes and entire gut microbiome dysbiosis at the genus level with both 16S rDNA and metagenomic sequencing techniques. We found that Negativibacillus and Rikenella were upregulated, while Akkermansia, Candidatus, Erysipelatoclostridium and Ileibacterium were downregulated in the colon of db/db mice compared with non-diabetic controls. In parallel, a total of 91 metabolites were upregulated, while 23 were downregulated in the plasma of db/db mice. The top five upregulated metabolites included D-arabinose 5-phosphate, estrone 3-sulfate, L-theanine, 3′-aenylic acid and adenosine 5′-monophosphate, and the five most significantly downregulated metabolites were aurohyocholic acid sodium salt, calcium phosphorylcholine chloride, tauro-alpha-muricholic acid sodium salt, galactinol and phosphocholine. These plasma metabolites were interacted with intestinal microbiomes, which are mainly involved in the pathways related to the biosynthesis of unsaturated fatty acids, fatty acid elongation, steroid biosynthesis, and D-arginine and D-ornithine metabolism. In the differential metabolites, N-acetyl-L-ornithine, ornithine and L-kyn could be metabolized by the correspondingly differential ontology genes in the intestinal metagenome. The current study thereby provides evidence for a gut–metabolism–kidney axis in the metabolism of db/db mice, in which the gut microbiome and circulating metabolomics interact, and suggests that information from this axis may contribute to our understanding of T2D and DKD pathogenesis.

A droplet-based microfluidic approach to isolating functional bacteria from gut microbiota

Metabolic interactions within gut microbiota play a vital role in human health and disease. Targeting metabolically interacting bacteria could provide effective treatments; however, obtaining functional bacteria remains a significant challenge due to the complexity of gut microbiota. Here, we developed a facile droplet-based approach to isolate and enrich functional gut bacteria that could utilize metabolites from an engineered butyrate-producing bacteria (EBPB) of anti-obesity potential. This involves the high throughput formation of single-bacteria droplets, followed by culturing “droplets” on agar plates to form discrete single-cell colonies. This approach eliminates the need for sophisticated s instruments to sort droplets and thus allows the operation hosted in a traditional anaerobic chamber. In comparison to the traditional culture, the droplet-based approach obtained a community of substantially higher diversity and evenness. Using the conditioned plates containing metabolites from the EBPB supernatant, we obtained gut bacteria closely associated or interacting with the EBPB. These include anaerobic Lactobacillus and Bifidobacterium, which are often used as probiotics. The study illustrates the potential of our approach in the search for the associated bacteria within the gut microbiota and retrieving those yet-to-be cultured.

The Role of Thermokarst Lake Expansion in Altering the Microbial Community and Methane Cycling in Beiluhe Basin on Tibetan Plateau

One of the most significant environmental changes across the Tibetan Plateau (TP) is the rapid lake expansion. The expansion of thermokarst lakes affects the global biogeochemical cycles and local climate regulation by rising levels, expanding area, and increasing water volumes. Meanwhile, microbial activity contributes greatly to the biogeochemical cycle of carbon in the thermokarst lakes, including organic matter decomposition, soil formation, and mineralization. However, the impact of lake expansion on distribution patterns of microbial communities and methane cycling, especially those of water and sediment under ice, remain unknown. This hinders our ability to assess the true impact of lake expansion on ecosystem services and our ability to accurately investigate greenhouse gas emissions and consumption in thermokarst lakes. Here, we explored the patterns of microorganisms and methane cycling by investigating sediment and water samples at an oriented direction of expansion occurred from four points under ice of a mature-developed thermokarst lake on TP. In addition, the methane concentration of each water layer was examined. Microbial diversity and network complexity were different in our shallow points (MS, SH) and deep points (CE, SH). There are differences of microbial community composition among four points, resulting in the decreased relative abundances of dominant phyla, such as Firmicutes in sediment, Proteobacteria in water, Thermoplasmatota in sediment and water, and increased relative abundance of Actinobacteriota with MS and SH points. Microbial community composition involved in methane cycling also shifted, such as increases in USCγ, Methylomonas, and Methylobacter, with higher relative abundance consistent with low dissolved methane concentration in MS and SH points. There was a strong correlation between changes in microbiota characteristics and changes in water and sediment environmental factors. Together, these results show that lake expansion has an important impact on microbial diversity and methane cycling.

Contribution of omics to biopreservation: Toward food microbiome engineering

Biopreservation is a sustainable approach to improve food safety and maintain or extend food shelf life by using beneficial microorganisms or their metabolites. Over the past 20 years, omics techniques have revolutionised food microbiology including biopreservation. A range of methods including genomics, transcriptomics, proteomics, metabolomics and meta-omics derivatives have highlighted the potential of biopreservation to improve the microbial safety of various foods. This review shows how these approaches have contributed to the selection of biopreservation agents, to a better understanding of the mechanisms of action and of their efficiency and impact within the food ecosystem. It also presents the potential of combining omics with complementary approaches to take into account better the complexity of food microbiomes at multiple scales, from the cell to the community levels, and their spatial, physicochemical and microbiological heterogeneity. The latest advances in biopreservation through omics have emphasised the importance of considering food as a complex and dynamic microbiome that requires integrated engineering strategies to increase the rate of innovation production in order to meet the safety, environmental and economic challenges of the agri-food sector.

MRG Chip: A High-Throughput qPCR-Based Tool for Assessment of the Heavy Metal(loid) Resistome

Bacterial metal detoxification mechanisms have been well studied for centuries in pure culture systems. However, profiling metal resistance determinants at the community level is still a challenge due to the lack of comprehensive and reliable quantification tools. Here, a novel high-throughput quantitative polymerase chain reaction (HT-qPCR) chip, termed the metal resistance gene (MRG) chip, has been developed for the quantification of genes involved in the homeostasis of 9 metals. The MRG chip contains 77 newly designed degenerate primer sets and 9 published primer sets covering 56 metal resistance genes. Computational evaluation of the taxonomic coverage indicated that the MRG chip had a broad coverage matching 2 kingdoms, 29 phyla, 64 classes, 130 orders, 226 families, and 382 genera. Temperature gradient PCR and HT-qPCR verified that 57 °C was the optimal annealing temperature, with amplification efficiencies of over 94% primer sets achieving 80-110%, with R² > 0.993. Both computational evaluation and the melting curve analysis of HT-qPCR validated a high specificity. The MRG chip has been successfully applied to characterize the distribution of diverse metal resistance determinants in natural and human-related environments, confirming its wide scope of application. Collectively, the MRG chip is a powerful and efficient high-throughput quantification tool for exploring the microbial metal resistome.

Contamination with multiple heavy metals decreases microbial diversity and favors generalists as the keystones in microbial occurrence networks

Soil contamination with multiple heavy metals poses threats to human health and ecosystem functioning. Using the Nemerow pollution index, which considers the effects of multiple heavy metals, we compared the diversity and composition of bacteria, fungi and protists and their potential interactions in response to a multi-metal contamination gradient. Multi-metal contamination significantly altered the community composition of bacteria, fungi and protists, and the degree of alteration increased with increasing severity of contamination. The alpha-diversity of bacteria, fungi and protists significantly decreased with increasing contamination level. The dominant generalists, found in all soil samples, were Gammaproteobacteria, Chloroflexi and Bacillus sp, whereas the dominant specialists were Anaerolineaceae, Entoloma sp. and Sandonidae_X sp. The relative abundances of generalists were positively correlated, whereas those of specialists were negatively correlated, with the Nemerow pollution index. In addition, the complexity of the microbial co-occurrence network increased with increasing contamination level. Generalists, rather than specialists, were the keystones in the microbial co-occurrence network and played a crucial role in adaptation to multi-metal contamination through enhanced potential interactions within the entire microbiome. Our results provide insights into the ecological effects of multi-metal contamination on the soil microbiome and will help to develop bio-remediation technologies for contaminated soils.

Association analyses of host genetics, root-colonizing microbes, and plant phenotypes under different nitrogen conditions in maize

The root-associated microbiome (rhizobiome) affects plant health, stress tolerance, and nutrient use efficiency. However, it remains unclear to what extent the composition of the rhizobiome is governed by intraspecific variation in host plant genetics in the field and the degree to which host plant selection can reshape the composition of the rhizobiome. Here we quantify the rhizosphere microbial communities associated with a replicated diversity panel of 230 maize (Zea mays L.) genotypes grown in agronomically relevant conditions under high N (+N) and low N (-N) treatments. We analyze the maize rhizobiome in terms of 150 abundant and consistently reproducible microbial groups and we show that the abundance of many root-associated microbes is explainable by natural genetic variation in the host plant, with a greater proportion of microbial variance attributable to plant genetic variation in -N conditions. Population genetic approaches identify signatures of purifying selection in the maize genome associated with the abundance of several groups of microbes in the maize rhizobiome. Genome-wide association study was conducted using the abundance of microbial groups as rhizobiome traits, and identified n = 622 plant loci that are linked to the abundance of n = 104 microbial groups in the maize rhizosphere. In 62/104 cases, which is more than expected by chance, the abundance of these same microbial groups was correlated with variation in plant vigor indicators derived from high throughput phenotyping of the same field experiment. We provide comprehensive datasets about the three-way interaction of host genetics, microbe abundance, and plant performance under two N treatments to facilitate targeted experiments towards harnessing the full potential of root-associated microbial symbionts in maize production.

The Response of Ruminal Microbiota and Metabolites to Different Dietary Protein Levels in Tibetan Sheep on the Qinghai-Tibetan Plateau

Ruminal microbiota and metabolites play crucial roles in animal health and productivity. Exploring the dynamic changes and interactions between microbial community composition and metabolites is important for understanding ruminal nutrition and metabolism. Tibetan sheep (Ovis aries) are an important livestock resource on the Qinghai-Tibetan Plateau (QTP), and the effects of various dietary protein levels on ruminal microbiota and metabolites are still unknown. The aim of this study was to investigate the response of ruminal microbiota and metabolites to different levels of dietary protein in Tibetan sheep. Three diets with different protein levels (low protein 10.1%, medium protein 12.1%, and high protein 14.1%) were fed to Tibetan sheep. 16S rRNA gene sequencing and gas chromatography coupled with time-of-flight mass spectrometry (GC-TOF-MS) were used to study the profile changes in each group of ruminal microbes and metabolites, as well as the potential interaction between them. The rumen microbiota in all groups was dominated by the phyla Bacteroidetes and Firmicutes regardless of the dietary protein level. At the genus level, Prevotella_1, Rikenellaceae_RC9_gut_group and Prevotellaceae_UCG-001 were dominant. Under the same forage-to-concentrate ratio condition, the difference in the dietary protein levels had no significant impact on the bacterial alpha diversity index and relative abundance of the major phyla and genera in Tibetan sheep. Rumen metabolomics analysis revealed that dietary protein levels altered the concentrations of ruminal amino acids, carbohydrates and organic acids, and significantly affected tryptophan metabolism (p < 0.05). Correlation analysis of the microbiota and metabolites revealed positive and negative regulatory mechanisms. Overall, this study provides detailed information on rumen microorganisms and ruminal metabolites under different levels of dietary protein, which could be helpful in subsequent research for regulating animal nutrition and metabolism through nutritional interventions.

Advances and challenges in cataloging the human gut virome

The human gut virome, which is often referred to as the "dark matter" of the gut microbiome, remains understudied. A better understanding of the composition and variations of the gut virome across populations is critical for exploring its impact on diseases and health. A series of advances in the characterization of human gut virome have unveiled high genetic diversity and various functional potentials of gut viruses. Here, we summarize the recently available human gut virome databases and discuss their features, procedures, and challenges with the intention to provide a reference to researchers to use while choosing a profiling database. We also propose a "best practice" for cataloging the viral population.

Contrasting Diversity and Composition of Human Colostrum Microbiota in a Maternal Cohort With Different Ethnic Origins but Shared Physical Geography (Island Scale)

Colostrum represents an important source for the transfer of important commensal bacteria from mother to newborn and has a strong impact on the newborn’s health after birth. However, the composition of the colostrum microbiome is highly heterogeneous due to geographic factors and ethnicity (maternal, cultural, and subsistence factors). By analyzing the colostrum 16S rRNA gene full-length sequencing dataset in 97 healthy mothers (60 from Han, 37 from Li) from the Hainan island of China, we showed that the ethnic differences of the colostrum microbiome in a maternal cohort with different ethnic origins shared physical geography. Results indicated that the richness of microbial community in colostrum of Han women was higher than that of Li women, but there was no significant difference in Shannon index and invsimpson index between the two groups. Visualization analysis based on the distance showed an obvious ethnicity-associated structural segregation of colostrum microbiota. The relative abundance of Firmicutes was higher in the microbiota of the Han group than in Li’s, while Proteobacteria was on the contrary. At the genus level, the most dominant members of the Han and Li ethnic groups were Acinetobacter and Cupriavidus, two common environmental bacteria, respectively, although skin-derived Staphylococcus and Streptococcus were still subdominant taxa. Cupriavidus lacunae was the most dominant species in the Li group, accounting for 26.10% of the total bacterial community, but only 3.43% for the Han group with the most dominant Staphylococcus petrasii (25.54%), indicating that human colostrum microbiome was more susceptible to local living environmental factors. Hence, the ethnic origin of individuals may be an important factor to consider in human milk microbiome research and its potential clinical significance during the perinatal period in ethnic-diverse societies, even within a small geographic scale.

Changes in the Mucosa-Associated Microbiome and Transcriptome across Gut Segments Are Associated with Obesity in a Metabolic Syndrome Porcine Model

Obesity is a major risk factor for metabolic syndrome, which is the most common cause of death worldwide, especially in developed countries. The link between obesity and gut mucosa-associated microbiota is unclear due to challenges associated with the collection of intestinal samples from humans.

Effect of carbendazim on honey bee health: Assessment of survival, pollen consumption, and gut microbiome composition

Gut microbiota and nutrition play major roles in honey bee health. Recent reports have shown that pesticides can disrupt the gut microbiota and cause malnutrition in honey bees. Carbendazim is the most commonly used fungicide in China, but it is not clear whether carbendazim negatively affects the gut microbes and nutrient intake levels in honey bees. To address this research gap, we assessed the effects of carbendazim on the survival, pollen consumption, and sequenced 16 S rRNA gene to determine the bacterial composition in the midgut and hindgut. Our results suggest that carbendazim exposure does not cause acute death in honey bees even at high concentrations (5000 mg/L), which are extremely unlikely to exist under field conditions. Carbendazim does not disturb the microbiome composition in the gut of young worker bees during gut microbial colonization and adult worker bees with established gut communities in the mid and hindgut. However, carbendazim exposure significantly decreases pollen consumption in honey bees. Thus, exposure of bees to carbendazim can perturb their beneficial nutrition homeostasis, potentially reducing honey bee immunity and increasing their susceptibility to infection by pathogens, which influence effectiveness as pollinators, even colony health.

Fungal diversity in shade-coffee plantations in Soconusco, Mexico

Background

As forested natural habitats disappear in the world, traditional, shade-coffee plantations offer an opportunity to conserve biodiversity and ecosystem services. Traditional coffee plantations maintain a diversity of tree species that provide shade for coffee bushes and, at the same time, are important repositories for plants and animals that inhabited the original cloud forest. However, there is still little information about shade-coffee plantation's fungal diversity despite their relevance for ecosystem functioning as decomposers, symbionts and pathogens. Specifically, it is unknown if and what mycorrhizae-forming fungi can be found on the branches and trunks of coffee bushes and trees, which hold a diversity of epiphytes. Here, we evaluate fungal communities on specific plant microsites on both coffee bushes and shade trees. We investigate the ecological roles played by this diversity, with a special focus on mycorrhizae-forming fungi that may enable the establishment and development of epiphytic plants.

Methods

We collected 48 bark samples from coffee bushes and shade trees (coffee; tree), from four plant microsites (upper and lower trunks, branches and twigs), in two shade-coffee plantations in the Soconusco region in southern Mexico, at different altitudes. We obtained ITS amplicon sequences that served to estimate alpha and beta diversity, to assign taxonomy and to infer the potential ecological role played by the detected taxa.

Results

The bark of shade trees and coffee bushes supported high fungal diversity (3,783 amplicon sequence variants). There were no strong associations between community species richness and collection site, plant type or microsite. However, we detected differences in beta diversity between collection sites. All trophic modes defined by FUNGuild database were represented in both plant types. However, when looking into guilds that involve mycorrhizae formation, the CLAM test suggests that coffee bushes are more likely to host taxa that may function as mycorrhizae.

Discussion

We detected high fungal diversity in shade-coffee plantations in Soconusco, Chiapas, possibly remnants of the original cloud forest ecosystem. Several mycorrhiza forming fungi guilds occur on the bark of coffee bushes and shade trees in this agroecosystem, with the potential of supporting epiphyte establishment and development. Thus, traditional coffee cultivation could be part of an integrated strategy for restoration and conservation of epiphytic populations. This is particularly relevant for conservation of threatened species of Orchidaceae that are highly dependent on mycorrhizae formation.

Bacterial Signatures of Cerebral Thrombi in Large Vessel Occlusion Stroke

It is important to understand the microbial features of the cerebral thrombus and its clinical relevance in stroke patients, of which data were scarce. We aimed to investigate the microbial features of cerebral thrombi retrieved via thrombectomy in stroke patients with large vessel occlusion (LVO) and their correlations with 3-month mortality. In a prospective cohort study, thrombus samples were collected during mechanical thrombectomy in LVO stroke patients with successful revascularization at a tertiary hospital. Oral, fecal, and isolated plasma samples were collected within 12 h of admission. The microbial compositions of all samples were compared using 16S rRNA gene amplicon next-generation sequencing. Fluorescent in situ hybridization (FISH) was used to detect bacteria in thrombus samples. The primary outcome was 3-month mortality. Perioperative adverse events (AEs) within 48 h were also recorded. Bacterial DNA was detected in 96.2% of thrombus samples from 104 patients, and clusters of bacterial signals were seen in the thrombi with FISH. Compared with fecal and oral samples, the thrombus microbiota was mainly characterized by excessive enrichment of Proteobacteria, mainly originating from plasma. The bacterial concentrations, dominant bacteria, and distribution patterns differed in thrombi obtained from cardioembolic and large-artery atherosclerotic strokes. Higher abundances of Acinetobacter and Enterobacteriaceae were associated with a higher risk of perioperative AEs, and a higher abundance of Acinetobacter was independently associated with a higher risk of 90-day mortality. This study demonstrated the presence of bacteria in cerebral thrombi retrieved with thrombectomy in LVO strokes, with some bacteria associated with patients’ prognoses.

All Set before Flowering: A 16S Gene Amplicon-Based Analysis of the Root Microbiome Recruited by Common Bean (Phaseolus vulgaris) in Its Centre of Domestication

Plant roots recruit most prokaryotic members of their root microbiota from the locally available inoculum, but knowledge on the contribution of native microorganisms to the root microbiota of crops in native versus non-native areas remains scarce. We grew common bean (Phaseolus vulgaris) at a field site in its centre of domestication to characterise rhizosphere and endosphere bacterial communities at the vegetative, flowering, and pod filling stage. 16S r RNA gene amplicon sequencing of ten samples yielded 9,401,757 reads, of which 8,344,070 were assigned to 17,352 operational taxonomic units (OTUs). Rhizosphere communities were four times more diverse than in the endosphere and dominated by Actinobacteria, Bacteroidetes, Crenarchaeota, and Proteobacteria (endosphere: 99% Proteobacteria). We also detected high abundances of Gemmatimonadetes (6%), Chloroflexi (4%), and the archaeal phylum Thaumarchaeota (Candidatus Nitrososphaera: 11.5%): taxa less frequently reported from common bean rhizosphere. Among 154 OTUs with different abundances between vegetative and flowering stage, we detected increased read numbers of Chryseobacterium in the endosphere and a 40-fold increase in the abundances of OTUs classified as Rhizobium and Aeromonas (equivalent to 1.5% and over 6% of all reads in the rhizosphere). Our results indicate that bean recruits specific taxa into its microbiome when growing 'at home'.

Distinct Patterns of Rhizosphere Microbiota Associated With Rice Genotypes Differing in Aluminum Tolerance in an Acid Sulfate Soil

Rhizosphere microbes are important for plant tolerance to various soil stresses. Rice is the most aluminum (Al)-tolerant small grain cereal crop species, but the link between rice Al tolerance and rhizosphere microbiota remains unclear. This study aimed to investigate the microbial community structure of aluminum-sensitive and Al-tolerant rice varieties in acid sulfate soil under liming and non-liming conditions. We analyzed the rice biomass and mineral element contents of rice plants as well as the chemical properties and microbial (archaea, bacteria, and fungi) communities of rhizosphere and bulk soil samples. The results showed that the Al-tolerant rice genotype grew better and was able to take up more phosphorus from the acid sulfate soil than the Al-sensitive genotype. Liming was the main factor altering the microbial diversity and community structure, followed by rhizosphere effects. In the absence of liming effects, the rice genotypes shifted the community structure of bacteria and fungi, which accounted for the observed variation in the rice biomass. The Al-tolerant rice genotype recruited specific bacterial and fungal taxa (Bacillus, Pseudomonas, Aspergillus, and Rhizopus) associated with phosphorus solubilization and plant growth promotion. The soil microbial co-occurrence network of the Al-tolerant rice genotype was more complex than that of the Al-sensitive rice genotype. In conclusion, the bacterial and fungal community in the rhizosphere has genotype-dependent effects on rice Al tolerance. Aluminum-tolerant rice genotypes recruit specific microbial taxa, especially phosphorus-solubilizing microorganisms, and are associated with complex microbial co-occurrence networks, which may enhance rice growth in acid sulfate soil.

Does bacterial community succession within the polyethylene mulching film plastisphere drive biodegradation?

Agricultural fields are severely contaminated with polyethylene mulching film (PMF) and this plastic in the natural environment can be colonized by biofilm-forming microorganisms that differ from those in the surrounding environment. In this study, we investigated the succession of the soil microbial communities in the PMF plastisphere using an artificial micro-ecosystem as well as exploring the degradation of PMF by plastisphere communities. The results indicated a significant and gradual decrease in the alpha diversity of the bacterial communities in the plastisphere and surrounding liquid. The community compositions in the plastisphere and surrounding liquid differed significantly from that in agricultural soil. Phyla and genera with the capacity to degrade polyethylene and hydrocarbon were enriched in the plastisphere, and some of these microorganisms were core members of the plastisphere community. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) analysis detected increases in metabolism pathways for PMF plastisphere Xenobiotics Biodegradation and Metabolism, thereby suggesting the possibility of polyethylene degradation in the plastisphere. Observations by scanning electron microscopy (SEM) and confocal laser scanning microscopy demonstrated the formation of biofilms on the incubated PMF. SEM, atomic force microscopy, Fourier transform infrared spectroscopy and water contact angle detected significant changes in the surface microstructure, chemical composition and hydrophobicity change of the films, thereby suggesting that the plastisphere community degraded PMF during incubation. In conclusion, this study provides insights into the changes in agricultural soil microorganisms in the PMF plastisphere and the degradation of PMF.

Obese Individuals With and Without Phlegm-Dampness Constitution Show Different Gut Microbial Composition Associated With Risk of Metabolic Disorders

Background

Obesity is conventionally considered a risk factor for multiple metabolic diseases, such as dyslipidemia, type 2 diabetes, hypertension, and cardiovascular disease (CVD). However, not every obese patient will progress to metabolic disease. Phlegm-dampness constitution (PDC), one of the nine TCM constitutions, is considered a high-risk factor for obesity and its complications. Alterations in the gut microbiota have been shown to drive the development and progression of obesity and metabolic disease, however, key microbial changes in obese patients with PDC have a higher risk for metabolic disorders remain elusive.

Methods

We carried out fecal 16S rRNA gene sequencing in the present study, including 30 obese subjects with PDC (PDC), 30 individuals without PDC (non-PDC), and 30 healthy controls with balanced constitution (BC). Metagenomic functional prediction of bacterial taxa was achieved using PICRUSt.

Results

Obese individuals with PDC had higher BMI, waist circumference, hip circumference, and altered composition of their gut microbiota compared to non-PDC obese individuals. At the phylum level, the gut microbiota was characterized by increased abundance of Bacteroidetes and decreased levels of Firmicutes and Firmicutes/Bacteroidetes ratio. At the genus level, Faecalibacterium, producing short-chain fatty acid, achieving anti-inflammatory effects and strengthening intestinal barrier functions, was depleted in the PDC group, instead, Prevotella was enriched. Most PDC-associated bacteria had a stronger correlation with clinical indicators of metabolic disorders rather than more severe obesity. The PICRUSt analysis demonstrated 70 significantly different microbiome community functions between the two groups, which were mainly involved in carbohydrate and amino acid metabolism, such as promoting Arachidonic acid metabolism, mineral absorption, and Lipopolysaccharide biosynthesis, reducing Arginine and proline metabolism, flavone and flavonol biosynthesis, Glycolysis/Gluconeogenesis, and primary bile acid biosynthesis. Furthermore, a disease classifier based on microbiota was constructed to accurately discriminate PDC individuals from all obese people.

Conclusion

Our study shows that obese individuals with PDC can be distinguished from non-PDC obese individuals based on gut microbial characteristics. The composition of the gut microbiome altered in obese with PDC may be responsible for their high risk of metabolic diseases.