Enhancing diversity analysis by repeatedly rarefying next generation sequencing data describing microbial communities

Ecological processes of bacterial microbiome assembly in healthy and dysbiotic strawberry farms

The bacterial microbiome plays crucial role in plants' resistance to diseases, nutrient uptake and productivity. We examined the microbiome characteristics of healthy and unhealthy strawberry farms, focusing on soil (bulk soil, rhizosphere soil) and plant (roots and shoots). The relative abundance of most abundant taxa were correlated with the chemical soil properties and shoot niche revealed the least amount of significant correlations between the two. While alpha and beta diversities did not show differences between health groups, we identified a number of core taxa (16-59) and marker bacterial taxa for each healthy (Unclassified Tepidisphaerales, Ohtaekwangia, Hydrocarboniphaga) and dysbiotic (Udaeobacter, Solibacter, Unclassified Chitinophagales, Unclassified Nitrosomonadaceae, Nitrospira, Nocardioides, Tardiphaga, Skermanella, Pseudomonas, Allorhizobium-Neorhizobium-Pararhizobium-Rhizobium, Curtobacterium) niche. We also revealed selective pressure of strawberry rhizosphere soil and roots plants in unhealthy plantations increased stochastic ecological processes of bacterial microbiome assembly in shoots. Our findings contribute to understanding sustainable agriculture and plant-microbiome interactions.

A multi-bin rarefying method for evaluating alpha diversities in TCR sequencing data

MOTIVATION

T cell receptors (TCRs) constitute a major component of our adaptive immune system, governing the recognition and response to internal and external antigens. Studying the TCR diversity via sequencing technology is critical for a deeper understanding of immune dynamics. However, library sizes differ substantially across samples, hindering the accurate estimation/comparisons of alpha diversities. To address this, researchers frequently use an overall rarefying approach in which all samples are sub-sampled to an even depth. Despite its pervasive application, its efficacy has never been rigorously assessed.

RESULTS

In this paper, we develop an innovative "multi-bin" rarefying approach that partitions samples into multiple bins according to their library sizes, conducts rarefying within each bin for alpha diversity calculations, and performs meta-analysis across bins. Extensive simulations using real-world data highlight the inadequacy of the overall rarefying approach in controlling the confounding effect of library size. Our method proves robust in addressing library size confounding, outperforming competing normalization strategies by achieving better-controlled type-I error rates and enhanced statistical power in association tests.

AVAILABILITY AND IMPLEMENTATION

The code is available at https://github.com/mli171/MultibinAlpha. The datasets are freely available at https://doi.org/10.21417/B7001Z and https://doi.org/10.21417/AR2019NC.

Comparative analysis of gut microbiota in children with obstructive sleep apnea: assessing the efficacy of 16S rRNA gene sequencing in metabolic function prediction based on weight status

Background

Analyzing bacterial microbiomes consistently using next-generation sequencing (NGS) is challenging due to the diversity of synthetic platforms for 16S rRNA genes and their analytical pipelines. This study compares the efficacy of full-length (V1-V9 hypervariable regions) and partial-length (V3-V4 hypervariable regions) sequencing of synthetic 16S rRNA genes from human gut microbiomes, with a focus on childhood obesity.

Methods

In this observational and comparative study, we explored the differences between these two sequencing methods in taxonomic categorization and weight status prediction among twelve children with obstructive sleep apnea.

Results

The full-length NGS method by Pacbio® identified 118 genera and 248 species in the V1-V9 regions, all with a 0% unclassified rate. In contrast, the partial-length NGS method by Illumina® detected 142 genera (with a 39% unclassified rate) and 6 species (with a 99% unclassified rate) in the V3-V4 regions. These approaches showed marked differences in gut microbiome composition and functional predictions. The full-length method distinguished between obese and non-obese children using the Firmicutes/Bacteroidetes ratio, a known obesity marker (p = 0.046), whereas the partial-length method was less conclusive (p = 0.075). Additionally, out of 73 metabolic pathways identified through full-length sequencing, 35 (48%) were associated with level 1 metabolism, compared to 28 of 61 pathways (46%) identified through the partial-length method. The full-length NGS also highlighted complex associations between body mass index z-score, three bacterial species (Bacteroides ovatus, Bifidobacterium pseudocatenulatum, and Streptococcus parasanguinis ATCC 15912), and 17 metabolic pathways. Both sequencing techniques revealed relationships between gut microbiota composition and OSA-related parameters, with full-length sequencing offering more comprehensive insights into associated metabolic pathways than the V3-V4 technique.

Conclusion

These findings highlight disparities in NGS-based assessments, emphasizing the value of full-length NGS with amplicon sequence variant analysis for clinical gut microbiome research. They underscore the importance of considering methodological differences in future meta-analyses.

Longitudinal multicompartment characterization of host-microbiota interactions in patients with acute respiratory failure

Critical illness can significantly alter the composition and function of the human microbiome, but few studies have examined these changes over time. Here, we conduct a comprehensive analysis of the oral, lung, and gut microbiota in 479 mechanically ventilated patients (223 females, 256 males) with acute respiratory failure. We use advanced DNA sequencing technologies, including Illumina amplicon sequencing (utilizing 16S and ITS rRNA genes for bacteria and fungi, respectively, in all sample types) and Nanopore metagenomics for lung microbiota. Our results reveal a progressive dysbiosis in all three body compartments, characterized by a reduction in microbial diversity, a decrease in beneficial anaerobes, and an increase in pathogens. We find that clinical factors, such as chronic obstructive pulmonary disease, immunosuppression, and antibiotic exposure, are associated with specific patterns of dysbiosis. Interestingly, unsupervised clustering of lung microbiota diversity and composition by 16S independently predicted survival and performed better than traditional clinical and host-response predictors. These observations are validated in two separate cohorts of COVID-19 patients, highlighting the potential of lung microbiota as valuable prognostic biomarkers in critical care. Understanding these microbiome changes during critical illness points to new opportunities for microbiota-targeted precision medicine interventions.

Zebrafish (Danio rerio) behavioral phenotypes not underscored by different gut microbiota

Different animal behavioral phenotypes maintained and selectively bred over multiple generations may be underscored by dissimilar gut microbial community compositions or not have any significant dissimilarity in community composition. Operating within the microbiota-gut-brain axis framework, we anticipated differences in gut microbiome profiles between zebrafish (Danio rerio) selectively bred to display the bold and shy personality types. This would highlight gut microbe-mediated effects on host behavior. To this end, we amplified and sequenced a fragment of the 16S rRNA gene from the guts of bold and shy zebrafish individuals (n=10) via Miseq. We uncovered no significant difference in within-group microbial diversity nor between-group microbial community composition of the two behavioral phenotypes. Interestingly, though not statistically different, we determined that the gut microbial community of the bold phenotype was dominated by Burkholderiaceae, Micropepsaceae, and Propionibacteriaceae. In contrast, the shy phenotype was dominated by Beijerinckaceae, Pirelullacaeae, Rhizobiales_Incertis_Sedis, and Rubinishaeraceae. The absence of any significant difference in gut microbiota profiles between the two phenotypes would suggest that in this species, there might exist a stable "core" gut microbiome, regardless of behavioral phenotypes, and or possibly, a limited role for the gut microbiota in modulating this selected-for host behavior. This is the first study to characterize the gut microbial community of distinct innate behavioral phenotypes of the zebrafish (that are not considered dysbiotic states) and not rely on antibiotic or probiotic treatments to induce changes in behavior. Such studies are crucial to our understanding of the modulating impacts of the gut microbiome on normative animal behavior.

Tracking arbuscular mycorrhizal fungi to their source: active inoculation and passive dispersal differentially affect community assembly in urban soils

Communities of arbuscular mycorrhizal (AM) fungi assemble passively over time via biotic and abiotic mechanisms. In degraded soils, AM fungal communities can assemble actively when humans manage mycorrhizas for ecosystem restoration. We investigated mechanisms of urban AM fungal community assembly in a 2-yr green roof experiment. We compared AM fungal communities in inoculated and uninoculated trays to samples from two potential sources: the inoculum and air. Active inoculation stimulated more distinct and diverse AM fungal communities, an effect that intensified over time. In the treatment trays, 45% of AM fungal taxa were detected in the inoculum, 2% were detected in aerial samples, 23% were detected in both inoculum and air, and 30% were not detected in either source. Passive dispersal of AM fungi likely resulted in the successful establishment of a small number of species, but active inoculation with native AM fungal species resulted in an immediate shift to a diverse and unique fungal community. When urban soils are constructed or modified by human activity, this is an opportunity for intervention with AM fungi that will persist and add diversity to that system.

Virus Quasispecies Rarefaction: Subsampling with or without Replacement?

In quasispecies diversity studies, the comparison of two samples of varying sizes is a common necessity. However, the sensitivity of certain diversity indices to sample size variations poses a challenge. To address this issue, rarefaction emerges as a crucial tool, serving to normalize and create fairly comparable samples. This study emphasizes the imperative nature of sample size normalization in quasispecies diversity studies using next-generation sequencing (NGS) data. We present a thorough examination of resampling schemes using various simple hypothetical cases of quasispecies showing different quasispecies structures in the sense of haplotype genomic composition, offering a comprehensive understanding of their implications in general cases. Despite the big numbers implied in this sort of study, often involving coverages exceeding 100,000 reads per sample and amplicon, the rarefaction process for normalization should be performed with repeated resampling without replacement, especially when rare haplotypes constitute a significant fraction of interest. However, it is noteworthy that different diversity indices exhibit distinct sensitivities to sample size. Consequently, some diversity indicators may be compared directly without normalization, or instead may be resampled safely with replacement.

Enhanced gut microbiome supplementation of essential amino acids in Diploptera punctata fed low-protein plant-based diet

Introduction

Building on our previous work, we investigate how dietary shifts affect gut microbial essential amino acid (EAA) provisioning in the lactating cockroach Diploptera punctata.

Method

To that end, we fed cockroaches three distinct diets: a plant-only Gari diet composed of starchy and granulated root tuber Yucca (Manihot esculenta), a dog food diet (DF), and a cellulose-amended dog food (CADF) diet. We anticipated that the high carbohydrate, low protein Gari would highlight increased microbial EAA supplementation to the host.

Results

By day 28, we observed distinct profiles of 14 bacterial families in the insect gut microbiomes of the three dietary groups. CADF-fed insects predominantly harbored cellulolytic and nitrogen-fixing bacteria families Streptococcaceae and Xanthomonadaceae. In contrast, Gari-fed insects were enriched in anaerobic lignocellulolytic bacteria families Paludibacteraceae and Dysgonomonadaceae, while DF-fed insects had a prevalence of proteolytic anaerobes Williamwhitmaniaceae and sulfate-reducing bacteria Desulfovibrionaceae. Furthermore, we confirmed significantly higher EAA supplementation in Gari-fed insects than in non-Gari-fed insects based on δ13C-EAA offsets between insect and their diets. The δ13C-EAA offsets between DF and CADF were nearly indistinguishable, highlighting the relevance of using the plant-based Gari in this experiment to unequivocally demonstrate this function in this insect. These results were underscored by lower standard metabolic rate (SMR) relative to the DF insect in Gari-fed (intermediate SMR and dietary quality) and CADF (least SMR and dietary quality) insects.

Discussion

The influence of diet on EAA provisioning and SMR responses in insects underscores the need for further exploration into the role of gut microbial functions in modulating metabolic responses.

Metal-oxide precipitation influences microbiome structure in hyporheic zones receiving acid rock drainage

Streams impacted by historic mining activity are characterized by acidic pH, unique microbial communities, and abundant metal-oxide precipitation, all of which can influence groundwater-surface water exchange. We investigate how metal-oxide precipitates and hyporheic mixing mediate the composition of microbial communities in two streams receiving acid-rock and mine drainage near Silverton, Colorado, USA. A large, neutral pH hyporheic zone facilitated the precipitation of metal particles/colloids in hyporheic porewaters. A small, low pH hyporheic zone, limited by the presence of a low-permeability, iron-oxyhydroxide layer known as ferricrete, led to the formation of steep geochemical gradients and high dissolved-metal concentrations. To determine how these two hyporheic systems influence microbiome composition, we installed well clusters and deployed in situ microcosms in each stream to sample porewaters and sediments for 16S rRNA gene sequencing. Results indicated that distinct hydrogeochemical conditions were present above and below the ferricrete in the low pH system. A positive feedback loop may be present in the low pH stream where microbially mediated precipitation of iron-oxides contributes to additional clogging of hyporheic pore spaces, separating abundant, iron-oxidizing bacteria (Gallionella spp.) above the ferricrete from rare, low-abundance bacteria below the ferricrete. Metal precipitates and colloids that formed in the neutral pH hyporheic zone were associated with a more diverse phylogenetic community of nonmotile, nutrient-cycling bacteria that may be transported through hyporheic pore spaces. In summary, biogeochemical conditions influence, and are influenced by, hyporheic mixing, which mediates the distribution of micro-organisms and, thus, the cycling of metals in streams receiving acid-rock and mine drainage.

IMPORTANCE

In streams receiving acid-rock and mine drainage, the abundant precipitation of iron minerals can alter how groundwater and surface water mix along streams (in what is known as the "hyporheic zone") and may shape the distribution of microbial communities. The findings presented here suggest that neutral pH streams with large, well-mixed hyporheic zones may harbor and transport diverse microorganisms attached to particles/colloids through hyporheic pore spaces. In acidic streams where metal oxides clog pore spaces and limit hyporheic exchange, iron-oxidizing bacteria may dominate and phylogenetic diversity becomes low. The abundance of iron-oxidizing bacteria in acid mine drainage streams has the potential to contribute to additional clogging of hyporheic pore spaces and the accumulation of toxic metals in the hyporheic zone. This research highlights the dynamic interplay between hydrology, geochemistry, and microbiology at the groundwater-surface water interface of acid mine drainage streams.

Sporadic diurnal fluctuations of cyanobacterial populations in oligotrophic temperate systems can prevent accurate characterization of change and risk in aquatic systems

Cyanobacteria increasingly threaten recreational water use and drinking water resources globally. They require dynamic monitoring to account for variability in their distribution arising from diel cycles associated with oscillatory vertical migration. While this has been discussed in marine and eutrophic freshwater contexts, reports of diurnal vertical migration of cyanobacteria in oligotrophic freshwater lakes are scant. Typical monitoring protocols do not reflect these dynamics and frequently focus only on surface water sampling approaches, and either ignore sampling time or recommend large midday timeframes (e.g., 10AM-3PM), thereby preventing accurate characterization of cyanobacterial community dynamics. To evaluate the impact of diurnal migrations and water column stratification on cyanobacterial abundance and composition, communities were characterized in a shallow well-mixed lake interconnected to a thermally stratified lake in the Turkey Lakes Watershed (Ontario, Canada) using amplicon sequencing of the 16S rRNA gene across a multi-time point sampling series in 2018 and 2022. This work showed that cyanobacteria are present in oligotrophic lakes and their community structure varies (i) diurnally, (ii) across the depth of the water column, (iii) interannually within the same lake and (iv) between different lakes that are closely interconnected within the same watershed. It underscored the need for integrating multi-timepoint, multi-depth discrete sampling guidance into lake and reservoir monitoring programs to describe cyanobacteria community dynamics and signal change to inform risk management associated with the potential for cyanotoxin production. Ignoring variability in cyanobacterial community dynamics (such as that reported herein) and reducing sample numbers can lead to a false sense of security and missed opportunities to identify and mitigate changes in trophic status and associated risks such as toxin or taste and odor production, especially in sensitive, oligotrophic systems.

Associations between dietary habits, socio-demographics and gut microbial composition in adolescents

The relationship between dietary habits and microbiota composition during adolescence has not been well examined. This is a crucial knowledge gap to fill considering that diet-microbiota interactions influence neurodevelopment, immune system maturation and metabolic regulation. This study examined the associations between diet and the gut microbiota in a school-based sample of 136 adolescents (Mage = 12·1 years; age range 11-13 years; 48 % female; 47 % Black, 38 % non-Hispanic White, 15 % Hispanic or other minorities) from urban, suburban and rural areas in the Southeast USA. Adolescents completed the Rapid Eating Assessment for Participants and provided stool samples for 16S ribosomal RNA gene sequencing. Parents reported their child and family socio-demographic characteristics. The associations between diet and socio-demographics with gut microbiota diversity and abundance were analysed using multivariable regression models. Child race and ethnicity, sex, socio-economic status and geographic locale contributed to variation within microbiota composition (β-diversity). Greater consumption of processed meat was associated with a lower microbial α-diversity after adjusting for socio-demographic variables. Multi-adjusted models showed that frequent consumption of nutrient-poor, energy-dense foods (e.g. sugar-sweetened beverages, fried foods, sweets) was negatively associated with abundances of genera in the family Lachnospiraceae (Anaerostipes, Fusicatenibacter and Roseburia), which are thought to play a beneficial role in host health through their production of short-chain fatty acids (SCFAs). These results provide new insights into the complex relationships among socio-demographic factors, diet and gut microbiota during adolescence. Adolescence may represent a critical window of opportunity to promote healthy eating practices that shape a homoeostatic gut microbiota with life-long benefits.

Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals

Cycads are known to host symbiotic cyanobacteria, including Nostocales species, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includes Aulosira spp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Considerations for the inclusion of metabarcoding data in the plant protection product risk assessment of the soil microbiome

There is increasing interest in further developing the plant protection product (PPP) environmental risk assessment, particularly within the European Union, to include the assessment of soil microbial community composition, as measured by metabarcoding approaches. However, to date, there has been little discussion as to how this could be implemented in a standardized, reliable, and robust manner suitable for regulatory decision-making. Introduction of metabarcoding-based assessments of the soil microbiome into the PPP risk assessment would represent a significant increase in the degree of complexity of the data that needs to be processed and analyzed in comparison to the existing risk assessment on in-soil organisms. The bioinformatics procedures to process DNA sequences into community compositional data sets currently lack standardization, while little information exists on how these data should be used to generate regulatory endpoints and the ways in which these endpoints should be interpreted. Through a thorough and critical review, we explore these challenges. We conclude that currently, we do not have a sufficient degree of standardization or understanding of the required bioinformatics and data analysis procedures to consider their use in an environmental risk assessment context. However, we highlight critical knowledge gaps and the further research required to understand whether metabarcoding-based assessments of the soil microbiome can be utilized in a statistically and ecologically relevant manner within a PPP risk assessment. Only once these challenges are addressed can we consider if and how we should use metabarcoding as a tool for regulatory decision-making to assess and monitor ecotoxicological effects on soil microorganisms within an environmental risk assessment of PPPs. Integr Environ Assess Manag 2024;20:337-358. © 2023 SETAC.

Aberrant colon metabolome and the sudden infant death syndrome

BACKGROUND

The Sudden Infant Death Syndrome (SIDS) has been associated with increased peripheral serotonin and an abnormal colonic microbiome, suggesting the colonic metabolome may also be abnormal. This study addresses this potential correlation by comparing colonic autopsy tissue from SIDS to age-matched non-SIDS controls.

METHODS

Untargeted metabolomic analysis by mass spectrometry is used to assess human colonic metabolomic differences including serotonin. Expression of genes associated with colonic serotonin synthesis and transport (TPH1, TPH2, DDC, SCL6A4) is measured by qRT-PCR. Microbiome analysis is performed to compare the SIDS and non-SIDS colonic microbiome.

RESULTS

Unsupervised hierarchical cluster and principal component analyses of metabolomic data shows increased variability in the SIDS cohort and separation of SIDS cases from the non-SIDS controls. There is a trend toward increased serotonin in the SIDS cohort but there is no significant difference in expression of the serotonin synthesis and transport genes between SIDS and non-SIDS control cohorts. Microbiome analysis shows no significant difference between the SIDS and non-SIDS control cohorts.

CONCLUSIONS

This study demonstrates increased variability in the colonic metabolome and a trend towards increased colonic serotonin in SIDS. The underlying cause of colon metabolomic variability, and its potential role in SIDS pathogenesis, warrants further investigation.

IMPACT STATEMENT

The key message of this article is that SIDS is associated with an aberrant colonic metabolome. This is a novel observation suggesting another component in the pathophysiology underlying SIDS. Investigation of why the colonic metabolome is aberrant may offer new insights to SIDS pathogenesis and new strategies to reduce risk.

Waste not, want not: revisiting the analysis that called into question the practice of rarefaction

In 2014, McMurdie and Holmes published the provocatively titled "Waste not, want not: why rarefying microbiome data is inadmissible." The claims of their study have significantly altered how microbiome researchers control for the unavoidable uneven sequencing depths that are inherent in modern 16S rRNA gene sequencing. Confusion over the distinction between the definitions of rarefying and rarefaction continues to cloud the interpretation of their results. More importantly, the authors made a variety of problematic choices when designing and analyzing their simulations. I identified 11 factors that could have compromised the results of the original study. I reproduced the original simulation results and assessed the impact of those factors on the underlying conclusion that rarefying data is inadmissible. Throughout, the design of the original study made choices that caused rarefying and rarefaction to appear to perform worse than they truly did. Most important were the approaches used to assess ecological distances, the removal of samples with low sequencing depth, and not accounting for conditions where sequencing effort is confounded with treatment group. Although the original study criticized rarefying for the arbitrary removal of valid data, repeatedly rarefying data many times (i.e., rarefaction) incorporates all the data. In contrast, it is the removal of rare taxa that would appear to remove valid data. Overall, I show that rarefaction is the most robust approach to control for uneven sequencing effort when considered across a variety of alpha and beta diversity metrics.IMPORTANCEOver the past 10 years, the best method for normalizing the sequencing depth of samples characterized by 16S rRNA gene sequencing has been contentious. An often cited article by McMurdie and Holmes forcefully argued that rarefying the number of sequence counts was "inadmissible" and should not be employed. However, I identified a number of problems with the design of their simulations and analysis that compromised their results. In fact, when I reproduced and expanded upon their analysis, it was clear that rarefaction was actually the most robust approach for controlling for uneven sequencing effort across samples. Rarefaction limits the rate of falsely detecting and rejecting differences between treatment groups. Far from being "inadmissible", rarefaction is a valuable tool for analyzing microbiome sequence data.

Anthropic disturbances impact the soil microbial network structure and stability to a greater extent than natural disturbances in an arid ecosystem

Growing pressure from climate change and agricultural land use is destabilizing soil microbial community interactions. Yet little is known about microbial community resistance and adaptation to disturbances over time. This hampers our ability to determine the recovery latency of microbial interactions after disturbances, with fundamental implications for ecosystem functioning and conservation measures. Here we examined the response of bacterial and fungal community networks in the rhizosphere of Haloxylon salicornicum (Moq.) Bunge ex Boiss. over the course of soil disturbances resulting from a history of different hydric constraints involving flooding-drought successions. An anthropic disturbance related to past agricultural use, with frequent successions of flooding and drought, was compared to a natural disturbance, i.e., an evaporation basin, with yearly flooding-drought successions. The anthropic disturbance resulted in a specific microbial network topology characterized by lower modularity and stability, reflecting the legacy of past agricultural use on soil microbiome. In contrast, the natural disturbance resulted in a network topology and stability close to those of natural environments despite the lower alpha diversity, and a different community composition compared to that of the other sites. These results highlighted the temporality in the response of the microbial community structure to disturbance, where long-term adaptation to flooding-drought successions lead to a higher stability than disturbances occurring over a shorter timescale.

Evaluating the Potential Safety Risk of Plant-Based Meat Analogues by Analyzing Microbial Community Composition

Plant-based meat analogues offer an environmentally and scientifically sustainable option as a substitute for animal-derived meat. They contribute to reducing greenhouse gas emissions, freshwater consumption, and the potential risks associated with zoonotic diseases linked to livestock production. However, specific processing methods such as extrusion or cooking, using various raw materials, can influence the survival and growth of spoilage and pathogenic microorganisms, resulting in differences between plant-based meat analogues and animal meat. In this study, the microbial communities in five different types of plant-based meat analogues were investigated using high-throughput sequencing. The findings revealed a diverse range of bacteria, including Cyanobacteria, Firmicutes, Proteobacteria, Bacteroidota, Actinobacteriota, and Chloroflexi, as well as fungi such as Ascomycota, Basidiomycota, Phragmoplastophyta, Vertebrata, and Mucoromycota. Additionally, this study analyzed microbial diversity at the genus level and employed phenotype prediction to evaluate the relative abundance of various bacterium types, including Gram-positive and Gram-negative bacteria, aerobic, anaerobic, and facultative anaerobic bacteria, as well as potential pathogenic bacteria. The insights gained from this study provide valuable information regarding the microbial communities and phenotypes of different plant-based meat analogues, which could help identify effective storage strategies to extend the shelf-life of these products.

Early life exposure to broccoli sprouts confers stronger protection against enterocolitis development in an immunological mouse model of inflammatory bowel disease

IMPORTANCE

To our knowledge, IL-10-KO mice have not previously been used to investigate the interactions of host, microbiota, and broccoli, broccoli sprouts, or broccoli bioactives in resolving symptoms of CD. We showed that a diet containing 10% raw broccoli sprouts increased the plasma concentration of the anti-inflammatory compound sulforaphane and protected mice to varying degrees against disease symptoms, including weight loss or stagnation, fecal blood, and diarrhea. Younger mice responded more strongly to the diet, further reducing symptoms, as well as increased gut bacterial richness, increased bacterial community similarity to each other, and more location-specific communities than older mice on the diet intervention. Crohn's disease disrupts the lives of patients and requires people to alter dietary and lifestyle habits to manage symptoms. The current medical treatment is expensive with significant side effects, and a dietary intervention represents an affordable, accessible, and simple strategy to reduce the burden of symptoms.

Microbial community structure and functional prediction in five full-scale industrial park wastewater treatment plants

The development of industrial parks has become an important global trend contributing significantly to economic and industrial growth. However, this growth comes at a cost, as the treatment of multisource industrial wastewater generated in these parks can be difficult owing to its complex composition. Microorganisms play a critical role in pollutant removal during industrial park wastewater treatment. Therefore, our study focused on the microbial communities in five full-scale industrial park wastewater treatment plants (WWTPs) with similar treatment processes and capacities. The results showed that denitrifying bacteria were dominant in almost every process section of all the plants, with heterotrophic denitrification being the main pathway. Moreover, autotrophic sulfur denitrification and methane oxidation denitrification may contribute to total nitrogen (TN) removal. In plants where the influent had low levels of COD and TN, dominant bacteria included oligotrophic microorganisms like Prosthecobacter (2.88 % ~ 10.02 %) and hgcI_clade (2.05 % ~ 9.49 %). Heavy metal metabolizing microorganisms, such as Norank_f__PHOS-HE36 (3.96 % ~ 5.36 %) and Sediminibacterium (1.86 % ~ 5.34 %), were prevalent in oxidation ditch and secondary settling tanks in certain plants. Functional Annotation of Prokaryotic Taxa (FAPROTAX) revealed that microbial communities in the regulation and hydrolysis tanks exhibited higher potential activity in the nitrogen (N) and sulfur (S) cycles than those in the oxidation ditch. Sulfate/sulfite reduction was common in most plants, whereas the potential occurrence of sulfide compounds and thiosulfate oxidation tended to be higher in plants with a relatively high sulfate concentration and low COD content in their influent. Our study provides a new understanding of the microbial community in full-scale industrial park WWTPs and highlights the critical role of microorganisms in the treatment of industrial wastewater.

Vaginal microbial dynamics and pathogen colonization in a humanized microbiota mouse model

Vaginal microbial composition is associated with differential risk of urogenital infection. Although Lactobacillus spp. are thought to confer protection against infection, the lack of in vivo models resembling the human vaginal microbiota remains a prominent barrier to mechanistic discovery. Using 16S rRNA amplicon sequencing of C57BL/6J female mice, we found that vaginal microbial composition varies within and between colonies across three vivaria. Noting vaginal microbial plasticity in conventional mice, we assessed the vaginal microbiome of humanized microbiota mice (HMbmice). Like the community structure in conventional mice, HMbmice vaginal microbiota clustered into community state types but, uniquely, HMbmice communities were frequently dominated by Lactobacillus or Enterobacteriaceae. Compared to conventional mice, HMbmice were less susceptible to uterine ascension by urogenital pathobionts group B Streptococcus (GBS) and Prevotella bivia. Although Escherichia and Lactobacillus both correlated with the absence of uterine GBS, vaginal pre-inoculation with exogenous HMbmouse-derived E. coli, but not Ligilactobacillus murinus, reduced vaginal GBS burden. Overall, HMbmice serve as a useful model to elucidate the role of endogenous microbes in conferring protection against urogenital pathogens.

Microbe-associated molecular pattern recognition receptors have little effect on endophytic Arabidopsis thaliana microbiome assembly in the field

Plant microbiome structure affects plant health and productivity. A limited subset of environmental microbes successfully establishes within plant tissues, but the forces underlying this selectivity remain poorly characterized. Transmembrane pattern recognition receptors (PRRs), used by plants to detect microbe-associated molecular patterns (MAMPs), are strong candidates for achieving this selectivity because PRRs can potentially interact with many members of the microbiome. Indeed, MAMPs found in many microbial taxa, including beneficials and commensals, can instigate a robust immune response that affects microbial growth. Surprisingly, we found that MAMP-detecting PRRs have little effect on endophytic bacterial and fungal microbiome structure in the field. We compared the microbiomes of four PRR knockout lines of Arabidopsis thaliana to wild-type plants in multiple tissue types over several developmental stages and detected only subtle shifts in fungal, but not bacterial, β-diversity in one of the four PRR mutants. In one developmental stage, lore mutants had slightly altered fungal β-diversity, indicating that LORE may be involved in plant-fungal interactions in addition to its known role in detecting certain bacterial lipids. No other effects of PRRs on α-diversity, microbiome variability, within-individual homogeneity, or microbial load were found. The general lack of effect suggests that individual MAMP-detecting PRRs are not critical in shaping the endophytic plant microbiome. Rather, we suggest that MAMP-detecting PRRs must either act in concert and/or are individually maintained through pleiotropic effects or interactions with coevolved mutualists or pathogens. Although unexpected, these results offer insights into the role of MAMP-detecting PRRs in plant-microbe interactions and help direct future efforts to uncover host genetic elements that control plant microbiome assembly.

Ecological diversity profiles of non-vaccine-targeted HPVs after gender-based community vaccination efforts

The long-term effect of population-level human papillomavirus (HPV) vaccination on the viral ecology of the untargeted HPVs is poorly understood. We performed an 8-year follow-up of 33 communities randomized to gender-neutral HPV16/18 vaccination, girls-only HPV16/18 vaccination, and control communities without HPV vaccination. The 1992/93 and 1994 birth cohorts were invited in school years 2007/8 and 2008/9. Follow-up cervico-vaginal sampling at 18 and 22 years of age, 4 and 8 years post-vaccination, respectively, were attended by 11,396 and 5,602 participants. HPV types 6/11/16/18/31/33/35/39/45/51/52/56/58/59/66/68 were genotyped and used for the community-level ecological diversity estimations. Gender-neutral vaccination communities with a stronger herd immunity than girls-only vaccination communities show a significantly increased HPV α-diversity (p = 1.1 × 10-8) from 4 to 8 years post-vaccination, despite the clearance of the vaccine-targeted HPVs in these communities. This likely sign of niche occupation by the non-vaccine-targeted HPVs will potentially affect the future cervical cancer screening programs but should not interfere with the WHO mission to eliminate cervical cancer.

Effect of Ginger Root Powder on Gastrointestinal Bacteria Composition, Gastrointestinal Symptoms, Mental Health, Fatigue, and Quality of Life: A Double-Blind Placebo-Controlled Trial

BACKGROUND

Despite compositional alterations in gastrointestinal microbiota being purported to underpin some of the therapeutic effects of ginger, the effect of a standardized ginger supplement on gut microbiota has not been tested in humans.

OBJECTIVES

To determine the effect of a standardized ginger (Zingiber officinale) root powder, compared to placebo, on gastrointestinal bacteria and associated outcomes in healthy adults.

METHODS

A randomized double-blind placebo-controlled trial allocated participants aged 18 to 30 y to ginger or microcrystalline cellulose (MCC) placebo. The intervention comprised 1.2 g/d of ginger (4 capsules per day totaling 84 mg/d of active gingerols/shogaols) for 14 d following a 1-wk run-in period. Primary outcomes were gastrointestinal community composition, alpha and beta diversity, and differential abundance, measured using 16S rRNA gene sequencing of fecal samples. Secondary outcomes were gastrointestinal symptoms, bowel function, depression, anxiety, stress, fatigue, quality of life, and adverse events.

RESULTS

Fifty-one participants were enrolled and analyzed (71% female; mean age 25 ± 3 y; ginger: n = 29, placebo: n = 22). There was a greater increase in relative abundance of phylum, Actinobacteria, observed following ginger supplementation compared to placebo (U: 145.0; z: -2.1; P = 0.033). Ginger was associated with a greater abundance of the genera Parabacteroides, Bacillus, Ruminococcaceae incertae sedis, unclassified Bacilli, families Defluviitaleaceae, Morganellaceae, and Bacillaceae as well as lower abundance of the genus Blautia and family Sphingomonadaceae (P < 0.05). An improvement in indigestion symptoms was observed with ginger supplementation (U: 196.0; z: -2.4; P = 0.015). No differences between ginger and placebo groups were found for alpha and beta diversity or other secondary outcomes. No moderate or severe adverse events were reported.

CONCLUSIONS

Supplementation with ginger root powder was safe and altered aspects of gastrointestinal bacteria composition; however, it did not change alpha- or beta diversity, bowel function, gastrointestinal symptoms, mood, or quality of life in healthy adults. These results provide further understanding regarding the mechanisms of action of ginger supplementation. This trial was registered in the Australia New Zealand Clinical Trials Registry as ACTRN12620000302954p and the Therapeutic Goods Administration as CT-2020-CTN-00380-1.

Monospecies probiotic preparation and administration with downstream analysis of sex-specific effects on gut microbiome composition in mice

Dysbiosis of the gut microbiome is implicated in the growing burden of non-communicable chronic diseases, including neurodevelopmental disorders, and both preclinical and clinical studies highlight the potential for precision probiotic therapies in their prevention and treatment. Here, we present an optimized protocol for the preparation and administration of Limosilactobacillus reuteri MM4-1A (ATCC-PTA-6475) to adolescent mice. We also describe steps for performing downstream analysis of metataxonomic sequencing data with careful assessment of sex-specific effects on microbiome composition and structure. For complete details on the use and execution of this protocol, please refer to Di Gesù et al.1.

Antibiotic resistance indicator genes in biofilm and planktonic microbial communities after wastewater discharge

The spread of bacteria with antibiotic resistance genes (ARGs) in aquatic ecosystems is of growing concern as this can pose a risk of transmission to humans and animals. While the impact of wastewater treatment plant (WWTP) effluent on ARG abundance in surface waters has been studied extensively, less is known about the fate of ARGs in biofilms. The proximity and dense growth of microorganisms in combination with the accumulation of higher antibiotic concentrations in biofilms might render biofilms a reservoir for ARGs. Seasonal parameters such as water temperature, precipitation, and antibiotic concentrations should be considered as well, as they may further influence the fate of ARGs in aquatic ecosystems. Here we investigated the effect of WWTP effluent on the abundance of the sulfonamide resistance genes sul1 and sul2, and the integrase gene intI1 in biofilm and surface water compartments of a river in Germany with a gradient of anthropogenic impact using quantitative PCR. Furthermore, we analyzed the bacterial community structure in both compartments via 16S rRNA gene amplicon sequencing, following the river downstream. Additionally, conventional water parameters and sulfonamide concentrations were measured, and seasonal aspects were considered by comparing the fate of ARGs and bacterial community diversity in the surface water compartment between the summer and winter season. Our results show that biofilm compartments near the WWTP had a higher relative abundance of ARGs (up to 4.7%) than surface waters (<2.8%). Sulfonamide resistance genes were more persistent further downstream (>10 km) of the WWTP in the hot and dry summer season than in winter. This finding is likely a consequence of the higher proportion of wastewater and thus wastewater-derived microorganisms in the river during summer periods. We observed distinct bacterial communities and ARG abundance between the biofilm and surface water compartment, but even greater variations when considering seasonal and spatiotemporal parameters. This underscores the need to consider seasonal aspects when studying the fate of ARGs in aquatic ecosystems.

Beneficial microbial consortium improves winter rye performance by modulating bacterial communities in the rhizosphere and enhancing plant nutrient acquisition

The beneficial effect of microbial consortium application on plants is strongly affected by soil conditions, which are influenced by farming practices. The establishment of microbial inoculants in the rhizosphere is a prerequisite for successful plant-microorganism interactions. This study investigated whether a consortium of beneficial microorganisms establishes in the rhizosphere of a winter crop during the vegetation period, including the winter growing season. In addition, we aimed for a better understanding of its effect on plant performance under different farming practices. Winter rye plants grown in a long-time field trial under conventional or organic farming practices were inoculated after plant emergence in autumn with a microbial consortium containing Pseudomonas sp. (RU47), Bacillus atrophaeus (ABi03) and Trichoderma harzianum (OMG16). The density of the microbial inoculants in the rhizosphere and root-associated soil was quantified in autumn and the following spring. Furthermore, the influence of the consortium on plant performance and on the rhizosphere bacterial community assembly was investigated using a multidisciplinary approach. Selective plating showed a high colonization density of individual microorganisms of the consortium in the rhizosphere and root-associated soil of winter rye throughout its early growth cycle. 16S rRNA gene amplicon sequencing showed that the farming practice affected mainly the rhizosphere bacterial communities in autumn and spring. However, the microbial consortium inoculated altered also the bacterial community composition at each sampling time point, especially at the beginning of the new growing season in spring. Inoculation of winter rye with the microbial consortium significantly improved the plant nutrient status and performance especially under organic farming. In summary, the microbial consortium showed sufficient efficacy throughout vegetation dormancy when inoculated in autumn and contributed to better plant performance, indicating the potential of microbe-based solutions in organic farming where nutrient availability is limited.

Profiles of bacterial communities and environmental factors associated with proliferation of malaria vector mosquitoes within the Kenyan Coast

Background

Since Anopheles mosquitoes which transmit and maintain the malaria parasite breed in the outdoor environment, there is an urgent need to manage these mosquito breeding sites. In order to elaborate more on the ecological landscape of mosquito breeding sites, the bacterial community structure and their interactions with physicochemical factors in mosquito larval habitats was characterised in Kwale County (Kenya), where malaria is endemic.

Methods

The physical characteristics and water physicochemical parameters of the habitats were determined and recorded. Water samples were also collected from the identified sites for total metagenomic DNA extraction in order to characterise the bacterial communities within the breeding sites.

Results and Discussion

Sites where mosquito larvae were found were described as positive and those without mosquito larvae as negative. Electrical conductivity, total dissolved solids, salinity and ammonia were lower in the rainy season than in the dry season, which also coincided with a high proportion of positive sites. Pseudomonadota was the most common phyla recovered in all samples followed by Bacteroidota and then Actinomycetota. The presence or absence of mosquito larvae in a potential proliferation site was not related to the bacterial community structure in the sampled sites, but was positively correlated with bacterial richness and evenness.

Conclusion

Generally, the presence of Anopheles mosquito larvae was found to be positively correlated with rainy season, bacterial richness and evenness, and negatively correlated with electrical conductivity, total dissolved solids, salinity and ammonia. The findings of this study have implications for predicting the potential of environmental water samples to become mosquito proliferation sites.

Systematic review of cnidarian microbiomes reveals insights into the structure, specificity, and fidelity of marine associations

Microorganisms play essential roles in the health and resilience of cnidarians. Understanding the factors influencing cnidarian microbiomes requires cross study comparisons, yet the plethora of protocols used hampers dataset integration. We unify 16S rRNA gene sequences from cnidarian microbiome studies under a single analysis pipeline. We reprocess 12,010 cnidarian microbiome samples from 186 studies, alongside 3,388 poriferan, 370 seawater samples, and 245 cultured Symbiodiniaceae, unifying ~6.5 billion sequence reads. Samples are partitioned by hypervariable region and sequencing platform to reduce sequencing variability. This systematic review uncovers an incredible diversity of 86 archaeal and bacterial phyla associated with Cnidaria, and highlights key bacteria hosted across host sub-phylum, depth, and microhabitat. Shallow (< 30 m) water Alcyonacea and Actinaria are characterized by highly shared and relatively abundant microbial communities, unlike Scleractinia and most deeper cnidarians. Utilizing the V4 region, we find that cnidarian microbial composition, richness, diversity, and structure are primarily influenced by host phylogeny, sampling depth, and ocean body, followed by microhabitat and sampling date. We identify host and geographical generalist and specific Endozoicomonas clades within Cnidaria and Porifera. This systematic review forms a framework for understanding factors governing cnidarian microbiomes and creates a baseline for assessing stress associated dysbiosis.

Oral microbiome and preterm birth

Background

The etiology of preterm birth (PTB) is heterogeneous and not yet well known. Maternal periodontal disease has been investigated for decades and is a known risk factor for adverse pregnancy outcomes. However, no particular bacterial species or higher taxonomic order has been found as causative of PTB, leading to studies of the whole oral microbiome. In order to determine if and how the composition of the oral microbiome is associated with PTB, we performed a large case-control study including women with term (TB) and PTB.

Methods

We compared oral microbiomes in PTB to TB, to examine differences in the microbial richness, diversity, and differential abundance of specific taxa. We obtained oral swab samples from 152 Caucasian pregnant women who were classified as either PTB (≤36 6/7 weeks, n = 61) or TB (≥38 0/7 weeks, n = 91) in exclusion of any other major medical or obstetric conditions. The oral microbiomes of these women were characterized by 16S ribosomal RNA (rRNA) gene sequencing of the V3-V4 region on the MiSeq platform.

Results

The dominant microorganisms at the phylum level in all pregnant women regardless of birth week outcomes as belonging to Firmicutes, Proteobacteria, Bacteroidetes, Fusobacteria, and Actinobacteria. The phyla Firmicutes and Bacteroidetes were relatively more abundant in women with a PTB than in women with a TB, while Proteobacteria was less prevalent in women with a PTB. At the genus level, Veillonella, Prevotella, and Capnocytophaga were enriched in the PTB, and while many of the members of these genera could not be resolved to the species level, Veillonella massillensis was shown to be increased in the PTB group.

Conclusion

We identified the genera Veillonella, Prevotella, and Capnocytophaga in the maternal oral microbiome as being associated with PTB independently of clinically apparent infection, uterine anomalies, and other pregnancy complications, including placenta previa, and placental abruption. The clarification of the role of those taxa in the etiology of PTB merits further research.

Machine learning and metagenomics reveal shared antimicrobial resistance profiles across multiple chicken farms and abattoirs in China

China is the largest global consumer of antimicrobials and improving surveillance methods could help to reduce antimicrobial resistance (AMR) spread. Here we report the surveillance of ten large-scale chicken farms and four connected abattoirs in three Chinese provinces over 2.5 years. Using a data mining approach based on machine learning, we analysed 461 microbiomes from birds, carcasses and environments, identifying 145 potentially mobile antibiotic resistance genes (ARGs) shared between chickens and environments across all farms. A core set of 233 ARGs and 186 microbial species extracted from the chicken gut microbiome correlated with the AMR profiles of Escherichia coli colonizing the same gut, including Arcobacter, Acinetobacter and Sphingobacterium, clinically relevant for humans, and 38 clinically relevant ARGs. Temperature and humidity in the barns were also correlated with ARG presence. We reveal an intricate network of correlations between environments, microbial communities and AMR, suggesting multiple routes to improving AMR surveillance in livestock production.

Early life exposure to broccoli sprouts confers stronger protection against enterocolitis development in an immunological mouse model of inflammatory bowel disease

Crohn's Disease (CD) is a presentation of Inflammatory Bowel Disease (IBD) that manifests in childhood and adolescence, and involves chronic and severe enterocolitis, immune and gut microbiome dysregulation, and other complications. Diet and gut-microbiota-produced metabolites are sources of anti-inflammatories which could ameliorate symptoms. However, questions remain on how IBD influences biogeographic patterns of microbial location and function in the gut, how early life transitional gut communities are affected by IBD and diet interventions, and how disruption to biogeography alters disease mediation by diet components or microbial metabolites. Many studies on diet and IBD use a chemically induced ulcerative colitis model, despite the availability of an immune-modulated CD model. Interleukin-10-knockout (IL-10-KO) mice on a C57BL/6 background, beginning at age 4 or 7 weeks, were fed a control diet or one containing 10% (w/w) raw broccoli sprouts, which was high in the sprout-sourced anti-inflammatory sulforaphane. Diets began 7 days prior to, and for 2 weeks after inoculation with Helicobacter hepaticus, which triggers Crohn's-like symptoms in these immune-impaired mice. The broccoli sprout diet increased sulforaphane in plasma; decreased weight stagnation, fecal blood, and diarrhea associated; and increased microbiota richness in the gut, especially in younger mice. Sprout diets resulted in some anatomically specific bacteria in younger mice, and reduced the prevalence and abundance of pathobiont bacteria which trigger inflammation in the IL-10-KO mouse, for example; Escherichia coli and Helicobacter. Overall, the IL-10-KO mouse model is responsive to a raw broccoli sprout diet and represents an opportunity for more diet-host-microbiome research.

Evaluating the role of bacterial diversity in supporting soil ecosystem functions under anthropogenic stress

Ecosystem functions and services are under threat from anthropogenic global change at a planetary scale. Microorganisms are the dominant drivers of nearly all ecosystem functions and therefore ecosystem-scale responses are dependent on responses of resident microbial communities. However, the specific characteristics of microbial communities that contribute to ecosystem stability under anthropogenic stress are unknown. We evaluated bacterial drivers of ecosystem stability by generating wide experimental gradients of bacterial diversity in soils, applying stress to the soils, and measuring responses of several microbial-mediated ecosystem processes, including C and N cycling rates and soil enzyme activities. Some processes (e.g., C mineralization) exhibited positive correlations with bacterial diversity and losses of diversity resulted in reduced stability of nearly all processes. However, comprehensive evaluation of all potential bacterial drivers of the processes revealed that bacterial α diversity per se was never among the most important predictors of ecosystem functions. Instead, key predictors included total microbial biomass, 16S gene abundance, bacterial ASV membership, and abundances of specific prokaryotic taxa and functional groups (e.g., nitrifying taxa). These results suggest that bacterial α diversity may be a useful indicator of soil ecosystem function and stability, but that other characteristics of bacterial communities are stronger statistical predictors of ecosystem function and better reflect the biological mechanisms by which microbial communities influence ecosystems. Overall, our results provide insight into the role of microorganisms in supporting ecosystem function and stability by identifying specific characteristics of bacterial communities that are critical for understanding and predicting ecosystem responses to global change.

Metagenomic comparisons reveal a highly diverse and unique viral community in a seasonally fluctuating hypersaline microbial mat

In spring 2016, a shallow hypersaline pond (50×25 m) was found in the Cuatro Cienegas Basin (CCB). This pond, known as Archaean Domes (AD) because of its elastic microbial mats that form dome-shaped structures due to the production of reducing gases reminiscent of the Archaean eon, such as methane and hydrogen sulfide, harbour a highly diverse microbial community, rich in halophilic and methanogenic archaea. AD is a seasonally fluctuating hypersaline site, with salinity ranging from low hypersaline (5.3%) during the wet season to high hypersaline (saturation) during the dry season. To characterize the viral community and to test whether it resembles those of other hypersaline sites (whose diversity is conditioned by salinity), or if it is similar to other CCB sites (with which it shares a common geological history), we generated 12 metagenomes from different seasons and depths over a 4 year period and compared them to 35 metagenomes from varied environments. Haloarchaeaviruses were detected, but were never dominant (average of 15.37 % of the total viral species), and the viral community structure and diversity were not affected by environmental fluctuations. In fact, unlike other viral communities at hypersaline sites, AD remained more diverse than other environments regardless of season. β-Diversity analyses show that AD is closely related to other CCB sites, although it has a unique viral community that forms a cluster of its own. The similarity of two surface samples to the 30 and 50 cm depth samples, as well as the observed increase in diversity at greater depths, supports the hypothesis that the diversity of CCB has evolved as a result of a long time environmental stability of a deep aquifer that functions as a 'seed bank' of great microbial diversity that is transported to the surface by sporadic groundwater upwelling events.

Evaluating the Reproducibility of Amplicon Sequencing Data Derived from Deep-Sea Cold Seep Sediment-Associated Microbiota

Benefiting from the rapid developments and wide applications of high-throughput sequencing, great advancements have been made in investigating microbiota, which are highly diverse and play key roles in both element cycling and the energy flow of ecosystems. There have been inherent limitations of amplicon sequencing that could introduce uncertainty and raise concerns about the accuracy and reproducibility of this technology. However, studies focusing on the reproducibility of amplicon sequencing are limited, especially in characterizing microbial communities in deep-sea sediments. To evaluate reproducibility, 118 deep-sea sediment samples were used for 16S rRNA gene sequencing in technical replicates (repeated measurements of the same sample) that demonstrate the variability of amplicon sequencing. The average occurrence-based overlaps were 35.98% and 27.02% between two and among three technical replicates, respectively, whereas their abundance-based overlaps reached 84.88% and 83.16%, respectively. Although variations of alpha and beta diversity indices were found between/among technical replicates, alpha diversity indices were similar across samples, and the average beta diversity indices were much smaller for technical replicates than among samples. Moreover, clustering methods (i.e., operational taxonomic units [OTUs] and amplicon sequence variants [ASVs]) were shown to have little impact on the alpha and beta diversity patterns of microbial communities. Taken together, although there are variations between/among technical replicates, amplicon sequencing is still a powerful tool with which to reveal diversity patterns of microbiota in deep-sea sediments. IMPORTANCE The reproducibility of amplicon sequencing is vital for whether the diversities of microbial communities could be accurately estimated. Thus, reproducibility influences the drawing of sound ecological conclusions. Nevertheless, few studies have focused on the reproducibility of microbial communities that are characterized by amplicon sequencing, and studies focusing on microbiota in deep-sea sediments have been especially lacking. In this study, we evaluated the reproducibility of amplicon sequencing targeting microbiota in deep-sea sediments of cold seep. Our results revealed that there were variations between/among technical replicates and that amplicon sequencing was still a powerful tool with which to characterize the diversities of microbial communities in deep-sea sediments. This study provides valuable guidelines for the reproducibility evaluation of future work in experimental design and interpretation.

Association between lipid-A-producing oral bacteria of different potency and fractional exhaled nitric oxide in a Norwegian population-based adult cohort

BACKGROUND

Lipid A is the primary immunostimulatory part of the lipopolysaccharide (LPS) molecule. The inflammatory response of LPS varies and depends upon the number of acyl chains and phosphate groups in lipid A which is specific for a bacterial species or strain. Traditional LPS quantification assays cannot distinguish between the acylation degree of lipid A molecules, and therefore little is known about how bacteria with different inflammation-inducing potencies affect fractional exhaled nitric oxide (F_eNO). We aimed to explore the association between pro-inflammatory hexa- and less inflammatory penta-acylated LPS-producing oral bacteria and F_eNO as a marker of airway inflammation.

METHODS

We used data from a population-based adult cohort from Norway (n = 477), a study center of the RHINESSA multi-center generation study. We applied statistical methods on the bacterial community- (prediction with MiRKAT) and genus-level (differential abundance analysis with ANCOM-BC) to investigate the association between the oral microbiota composition and F_eNO.

RESULTS

We found the overall composition to be significantly associated with increasing F_eNO levels independent of covariate adjustment, and abundances of 27 bacterial genera to differ in individuals with high F_eNO vs. low F_eNO levels. Hexa- and penta-acylated LPS producers made up 2.4% and 40.8% of the oral bacterial genera, respectively. The Bray-Curtis dissimilarity within hexa- and penta-acylated LPS-producing oral bacteria was associated with increasing F_eNO levels independent of covariate adjustment. A few single penta-acylated LPS producers were more abundant in individuals with low F_eNO vs. high F_eNO, while hexa-acylated LPS producers were found not to be enriched.

CONCLUSIONS

In a population-based adult cohort, F_eNO was observed to be associated with the overall oral bacterial community composition. The effect of hexa- and penta-acylated LPS-producing oral bacteria was overall significant when focusing on Bray-Curtis dissimilarity within each of the two communities and F_eNO levels, but only penta-acylated LPS producers appeared to be reduced or absent in individuals with high F_eNO. It is likely that the pro-inflammatory effect of hexa-acylated LPS producers is counteracted by the dominance of the more abundant penta-acylated LPS producers in this population-based adult cohort involving mainly healthy individuals.

Pathobionts in the tumour microbiota predict survival following resection for colorectal cancer

BACKGROUND AND AIMS

The gut microbiota is implicated in the pathogenesis of colorectal cancer (CRC). We aimed to map the CRC mucosal microbiota and metabolome and define the influence of the tumoral microbiota on oncological outcomes.

METHODS

A multicentre, prospective observational study was conducted of CRC patients undergoing primary surgical resection in the UK (n = 74) and Czech Republic (n = 61). Analysis was performed using metataxonomics, ultra-performance liquid chromatography-mass spectrometry (UPLC-MS), targeted bacterial qPCR and tumour exome sequencing. Hierarchical clustering accounting for clinical and oncological covariates was performed to identify clusters of bacteria and metabolites linked to CRC. Cox proportional hazards regression was used to ascertain clusters associated with disease-free survival over median follow-up of 50 months.

RESULTS

Thirteen mucosal microbiota clusters were identified, of which five were significantly different between tumour and paired normal mucosa. Cluster 7, containing the pathobionts Fusobacterium nucleatum and Granulicatella adiacens, was strongly associated with CRC (P_FDR = 0.0002). Additionally, tumoral dominance of cluster 7 independently predicted favourable disease-free survival (adjusted p = 0.031). Cluster 1, containing Faecalibacterium prausnitzii and Ruminococcus gnavus, was negatively associated with cancer (P_FDR = 0.0009), and abundance was independently predictive of worse disease-free survival (adjusted p = 0.0009). UPLC-MS analysis revealed two major metabolic (Met) clusters. Met 1, composed of medium chain (MCFA), long-chain (LCFA) and very long-chain (VLCFA) fatty acid species, ceramides and lysophospholipids, was negatively associated with CRC (P_FDR = 2.61 × 10^-11); Met 2, composed of phosphatidylcholine species, nucleosides and amino acids, was strongly associated with CRC (P_FDR = 1.30 × 10^-12), but metabolite clusters were not associated with disease-free survival (p = 0.358). An association was identified between Met 1 and DNA mismatch-repair deficiency (p = 0.005). FBXW7 mutations were only found in cancers predominant in microbiota cluster 7.

CONCLUSIONS

Networks of pathobionts in the tumour mucosal niche are associated with tumour mutation and metabolic subtypes and predict favourable outcome following CRC resection. Video Abstract.

Unraveling the Dysbiosis of Vaginal Microbiome to Understand Cervical Cancer Disease Etiology-An Explainable AI Approach

Microbial Dysbiosis is associated with the etiology and pathogenesis of diseases. The studies on the vaginal microbiome in cervical cancer are essential to discern the cause and effect of the condition. The present study characterizes the microbial pathogenesis involved in developing cervical cancer. Relative species abundance assessment identified Firmicutes, Actinobacteria, and Proteobacteria dominating the phylum level. A significant increase in Lactobacillus iners and Prevotella timonensis at the species level revealed its pathogenic influence on cervical cancer progression. The diversity, richness, and dominance analysis divulges a substantial decline in cervical cancer compared to control samples. The β diversity index proves the homogeneity in the subgroups' microbial composition. The association between enriched Lactobacillus iners at the species level, Lactobacillus, Pseudomonas, and Enterococcus genera with cervical cancer is identified by Linear discriminant analysis Effect Size (LEfSe) prediction. The functional enrichment corroborates the microbial disease association with pathogenic infections such as aerobic vaginitis, bacterial vaginosis, and chlamydia. The dataset is trained and validated with repeated k-fold cross-validation technique using a random forest algorithm to determine the discriminative pattern from the samples. SHapley Additive exPlanations (SHAP), a game theoretic approach, is employed to analyze the results predicted by the model. Interestingly, SHAP identified that the increase in Ralstonia has a higher probability of predicting the sample as cervical cancer. New evidential microbiomes identified in the experiment confirm the presence of pathogenic microbiomes in cervical cancer vaginal samples and their mutuality with microbial imbalance.

Evaluation of the rbcL marker for metabarcoding of marine diatoms and inference of population structure of selected genera

Diatoms are one of the most important phytoplankton groups in the world's oceans. There are responsible for up to 40% of the photosynthetic activity in the Ocean, and they play an important role in the silicon and carbon cycles by decoupling carbon from atmospheric interactions through sinking and export. These processes are strongly influenced by the taxonomic composition of diatom assemblages. Traditionally, these have been assessed using microscopy, which in some cases is not reliable or reproducible. Next-generation sequencing enabled us to study diversity in a high-throughput manner and uncover new distribution patterns and diversity. However, phylogenetic markers used for this purpose, such as various 18S rDNA regions, are often insufficient because they cannot distinguish between some taxa. In this work, we demonstrate the performance of the chloroplast-encoded rbcL marker for metabarcoding marine diatoms compared to microscopy and 18S-V9 metabarcoding using a series of monthly samples from the Gulf of Trieste (GoT), northern Adriatic Sea. We demonstrate that rbcL is able to detect more taxa compared to 18S-V9 metabarcoding or microscopy, while the overall structure of the diatom assemblage was comparable to the other two methods with some variations, that were taxon dependent. In total, 6 new genera and 22 new diatom species for the study region were identified. We were able to spot misidentification of genera obtained with microscopy such as Pseudo-nitzschia galaxiae, which was mistaken for Cylindrotheca closterium, as well as genera that were completely overlooked, such as Minidiscus and several genera from the Cymatosiraceae family. Furthermore, on the example of two well-studied genera in the region, namely Chaetoceros and particularly Pseudo-nitzschia, we show how the rbcL method can be used to infer even deeper phylogenetic and ecologically significant differences at the species population level. Despite a very thorough community analysis obtained by rbcL the incompleteness of reference databases was still evident, and we shed light on possible improvements. Our work has further implications for studies dealing with taxa distribution and population structure, as well as carbon and silica flux models and networks.

Unique Features of the Gut Microbiome Characterized in Animal Models of Angelman Syndrome

A large subset of patients with Angelman syndrome (AS) suffer from concurrent gastrointestinal (GI) issues, including constipation, poor feeding, and reflux. AS is caused by the loss of ubiquitin ligase E3A (UBE3A) gene expression in the brain. Clinical features of AS, which include developmental delays, intellectual disability, microcephaly, and seizures, are primarily due to the deficient expression or function of the maternally inherited UBE3A allele. The association between neurodevelopmental delay and GI disorders is part of the increasing evidence suggesting a link between the brain and the gut microbiome via the microbiota-gut-brain axis. To investigate the associations between colonization of the gut microbiota in AS, we characterized the fecal microbiome in three animal models of AS involving maternal deletions of Ube3A, including mouse, rat, and pig, using 16S rRNA amplicon sequencing. Overall, we identified changes in bacterial abundance across all three animal models of AS. Specific bacterial groups were significantly increased across all animal models, including Lachnospiraceae Incertae sedis, Desulfovibrios sp., and Odoribacter, which have been correlated with neuropsychiatric disorders. Taken together, these findings suggest that specific changes to the local environment in the gut are driven by a Ube3a maternal deletion, unaffected by varying housing conditions, and are prominent and detectable across multiple small and large animal model species. These findings begin to uncover the underlying mechanistic causes of GI disorders in AS patients and provide future therapeutic options for AS patients. IMPORTANCE Angelman syndrome (AS)-associated gastrointestinal (GI) symptoms significantly impact quality of life in patients. In AS models in mouse, rat, and pig, AS animals showed impaired colonization of the gut microbiota compared to wild-type (healthy) control animals. Common changes in AS microbiomes across all three animal models may play a causal effect for GI symptoms and may help to identify ways to treat these comorbidities in patients in the future.

Insights on the particle-attached riverine archaeal community shifts linked to seasons and to multipollution during a Mediterranean extreme storm event

Even if Archaea deliver important ecosystem services and are major players in global biogeochemical cycles, they remain poorly understood in freshwater ecosystems. To our knowledge, no studies specifically address the direct impact of xenobiotics on the riverine archaeome. Using environmental DNA metabarcoding of the 16S ribosomal gene, we previously demonstrated bacterial communities significant shifts linked to pollutant mixtures during an extreme flood in a typical Mediterranean coastal watercourse. Here, using the same methodology, we sought to determine whether archaeal community shifts coincided with the delivery of environmental stressors during the same flood. Further, we wanted to determine how archaea taxa compared at different seasons. In contrast to the bacteriome, the archaeome showed a specific community in summer compared to winter and autumn. We also identified a significant relationship between in situ archaeome shifts and changes in physicochemical parameters along the flood, but a less marked link to those parameters correlated to river hydrodynamics than bacteria. New urban-specific archaeal taxa significantly related to multiple stressors were identified. Through statistical modeling of both domains, our results demonstrate that Archaea, seldom considered as bioindicators of water quality, have the potential to improve monitoring methods of watersheds.

Maximising signal-to-noise ratios in environmental DNA-based monitoring

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

Impact of Data and Study Characteristics on Microbiome Volatility Estimates

The human microbiome is a dynamic community of bacteria, viruses, fungi, and other microorganisms. Both the composition of the microbiome (the microbes that are present and their relative abundances) and the temporal variability of the microbiome (the magnitude of changes in their composition across time, called volatility) has been associated with human health. However, the effect of unbalanced sampling intervals and differential read depth on the estimates of microbiome volatility has not been thoroughly assessed. Using four publicly available gut and vaginal microbiome time series, we subsampled the datasets to several sampling intervals and read depths and then compared additive, multiplicative, centered log ratio (CLR)-based, qualitative, and distance-based measures of microbiome volatility between the conditions. We find that longer sampling intervals are associated with larger quantitative measures of change (particularly for common taxa), but not with qualitative measures of change or distance-based volatility quantification. A lower sequencing read depth is associated with smaller multiplicative, CLR-based, and qualitative measures of change (particularly for less common taxa). Strategic subsampling may serve as a useful sensitivity analysis in unbalanced longitudinal studies investigating clinical associations with microbiome volatility.

Diagnostic Significance of Targeted Next-Generation Sequencing in Central Nervous System Infections in Neurosurgery of Pediatrics

Background

Cerebrospinal fluid (CSF) pathogen culture suffers from the drawbacks of prolonged cycle time and a low positivity rate in diagnosing intracranial infections in children. This study aims to investigate the diagnostic potential of targeted next-generation sequencing (tNGS) in pediatric neurosurgery for central nervous system (CNS) infections.

Methods

A retrospective study was conducted on children under 14 with suspected intracranial infections following craniocerebral trauma or surgery between November 2018 and August 2020. Routine, biochemical, smear, and pathogen culture tests were performed on CSF during treatment. The main parameters of CSF analysis encompassed white blood cells (WBC, ×10⁶/L) count, percentage of multinucleated cells (%), protein levels (g/L), glucose concentration (GLU, mmol/L), chloride levels (mmol/L), and pressure (mmH2O). The outcomes of tNGS were assessed through the Receiver Operating Characteristic (ROC) curve and pertinent diagnostic parameters.

Results

Among the 35 included pediatric patients, 22 were clinically diagnosed with CNS infection in neurosurgery, tNGS was confirmed in 18 cases. The sensitivity and specificity of tNGS were 81.8% and 76.9%, respectively, while the traditional method of CSF cultures and smears exhibited a sensitivity of 13.6% and a specificity of 100%. ROC curve analysis indicated an area under the curve (AUC) of 0.794 for tNGS and 0.568 for the CSF cultures and smears. CSF analysis indicated that the two groups exhibited statistically significant differences in terms of WBC count [330.0 (110.00-2639.75) vs 14.00 (4.50-26.50), P<0.001] and percentage of multinuclear cells (%) [87.50 (39.75-90.00) vs 0 (0-10.00), P<0.001]. However, the remaining parameters did not statistically significant differences between the groups (all P>0.05).

Conclusion

tNGS demonstrates a high degree of diagnostic accuracy when detecting infections within the CNS of pediatric neurosurgery patients. tNGS can effectively establish for diagnosing CNS infections by detecting pathogenic microorganisms and their corresponding virulence and/or resistance genes within the test samples.

Bioinformatic Approaches for Describing the Oral Microbiota

Cultivation-independent (molecular) analysis of the oral microbiota can provide a comprehensive picture of microbial community composition, yet there is an at-times bewildering array of approaches that can be employed. This chapter introduces some of the key considerations when undertaking microbiota research and describes two alternative bioinformatic pipelines for conducting such studies. The descriptions are based on analysis of bacterial 16S ribosomal RNA gene sequences, but can be easily adapted for analysis of other microbial taxa such as fungi.

Strong restructuration of skin microbiota during captivity challenges ex-situ conservation of amphibians

In response to the current worldwide amphibian extinction crisis, conservation instances have encouraged the establishment of ex-situ collections for endangered species. The resulting assurance populations are managed under strict biosecure protocols, often involving artificial cycles of temperature and humidity to induce active and overwintering phases, which likely affect the bacterial symbionts living on the amphibian skin. However, the skin microbiota is an important first line of defense against pathogens that can cause amphibian declines, such as the chytrid Batrachochytrium dendrobatidis (Bd). Determining whether current husbandry practices for assurance populations might deplete amphibians from their symbionts is therefore essential to conservation success. Here, we characterize the effect of the transitions from the wild to captivity, and between aquatic and overwintering phases, on the skin microbiota of two newt species. While our results confirm differential selectivity of skin microbiota between species, they underscore that captivity and phase-shifts similarly affect their community structure. More specifically, the translocation ex-situ is associated with rapid impoverishment, decrease in alpha diversity and strong species turnover of bacterial communities. Shifts between active and overwintering phases also cause changes in the diversity and composition of the microbiota, and on the prevalence of Bd-inhibitory phylotypes. Altogether, our results suggest that current husbandry practices strongly restructure the amphibian skin microbiota. Although it remains to be determined whether these changes are reversible or have deleterious effects on their hosts, we discuss methods to limit microbial diversity loss ex-situ and emphasize the importance of integrating bacterial communities to applied amphibian conservation.

Analysis of Bacterial Microbiota of Aerated Compost Teas and Effect on Tomato Growth

Mature composts and their water-based extracts, known as aerated compost teas (ACTs), are biofertilizers that share bioactive effects like soil restoration and plant health promotion, widely used for sustainable agriculture. Bioactive effects of compost and ACTs could be associated with their physicochemical and biological characteristics, like carbon/nitrogen (C/N) ratio and microbiota structure respectively. In our study, we elaborated ACTs using mature homemade compost, wheat bran, and grass clippings, following the C/N ratio criteria. Irrigation of tomato plantlets with ACT whose C/N ratio was close to the expected C/N ratio for mature compost evidenced plant growth promotion. Exploring the bacterial microbiota of elaborated ACTs and origin compost revealed significant structural differences, including phyla involved in N mineralization and free-living N-fixing bacteria. Therefore, ACTs harbor diverse bacterial microbiota involved in the N cycle, which would enrich plant and soil bacterial communities at the taxonomic and functional levels. Furthermore, ACTs are considered a part of agroecological and circular economy approaches.

Bacterial diversity across four drinking water distribution systems in Croatia: impacts of water management practices and disinfection by-products

Several factors may impact bacterial diversity in drinking water distribution systems (DWDSs) including the origin of the raw water, the water treatment technologies, and the disinfection practices applied. 16S rRNA metabarcoding was used for the in-depth characterization of bacterial communities in the four studied Croatian DWDSs (A, B, C, D) two of which had residual disinfectant (A, B) and two were without (C, D), while only B utilized the conventional water treatment technology. Significantly higher diversity and species richness were evidenced in non-disinfected DWDSs (p<0.05) compared to disinfected DWDSs. The phylum Proteobacteria was the most abundant in all the DWDSs, being proportionately higher in non-disinfected systems (p<0.05). The most abundant genera in DWDS-A Mycobacterium and Sphingomonas both positively correlated, whereas Lactobacillus negatively correlated with the concentration of disinfection by-products (DBPs) as a sum of haloacetic acids (HAAs). Conversely, the genus Ralstonia positively correlated with the individual DBP dichloroacetic acid. These results indicate that genera Sphingomonas, Mycobacterium, Lactobacillus and Ralstonia could have an effect on promoting the formation of DBPs, in a similar manner to how negatively correlated taxa may influence their degradation.

A longitudinal census of the bacterial community in raw milk correlated with Staphylococcus aureus clinical mastitis infections in dairy cattle

BACKGROUND

Staphylococcus aureus is a common cause of clinical mastitis (CM) in dairy cattle. Optimizing the bovine mammary gland microbiota to resist S. aureus colonization is a growing area of research. However, the details of the interbacterial interactions between S. aureus and commensal bacteria, which would be required to manipulate the microbiome to resist infection, are still unknown. This study aims to characterize changes in the bovine milk bacterial community before, during, and after S. aureus CM and to compare bacterial communities present in milk between infected and healthy quarters.

METHODS

We collected quarter-level milk samples from 698 Holstein dairy cows over an entire lactation. A total of 11 quarters from 10 cows were affected by S. aureus CM and milk samples from these 10 cows (n = 583) regardless of health status were analyzed by performing 16S rRNA gene amplicon sequencing.

RESULTS

The milk microbiota of healthy quarters was distinguishable from that of S. aureus CM quarters two weeks before CM diagnosis via visual inspection. Microbial network analysis showed that 11 OTUs had negative associations with OTU0001 (Staphylococcus). A low diversity or dysbiotic milk microbiome did not necessarily correlate with increased inflammation. Specifically, Staphylococcus xylosus, Staphylococcus epidermidis, and Aerococcus urinaeequi were each abundant in milk from the quarters with low levels of inflammation.

CONCLUSION

Our results show that the udder microbiome is highly dynamic, yet a change in the abundance in certain bacteria can be a potential indicator of future S. aureus CM. This study has identified potential prophylactic bacterial species that could act as a barrier against S. aureus colonization and prevent future instances of S. aureus CM.

Hidden Tenants: Microbiota of the Rhizosphere and Phyllosphere of Cordia dodecandra Trees in Mayan Forests and Homegardens

During domestication, the selection of cultivated plants often reduces microbiota diversity compared with their wild ancestors. Microbiota in compartments such as the phyllosphere or rhizosphere can promote fruit tree health, growth, and development. Cordia dodecandra is a deciduous tree used by Maya people for its fruit and wood, growing, to date, in remnant forest fragments and homegardens (traditional agroforestry systems) in Yucatán. In this work, we evaluated the microbiota's alpha and beta diversity per compartment (phyllosphere and rhizosphere) and per population (forest and homegarden) in the Northeast and Southwest Yucatán regions. Eight composite DNA samples (per compartment/population/region combination) were amplified for 16S-RNA (bacteria) and ITS1-2 (fungi) and sequenced by Illumina MiSeq. Bioinformatic analyses were performed with QIIME and phyloseq. For bacteria and fungi, from 107,947 and 128,786 assembled sequences, 618 and 1092 operating taxonomic units (OTUs) were assigned, respectively. The alpha diversity of bacteria and fungi was highly variable among samples and was similar among compartments and populations. A significant species turnover among populations and regions was observed in the rhizosphere. The core microbiota from the phyllosphere was similar among populations and regions. Forests and homegarden populations are reservoirs of the C. dodecandra phyllosphere core microbiome and significant rhizosphere biodiversity.