16S rRNA Gene Copies Research Articles

Fungal Abundance and Diversity in the Mariana Trench, the Deepest Ecosystem on Earth.

Hadal trenches host abundant and diversified benthic prokaryotic assemblages, but information on benthic fungi is still extremely limited. We investigated the fungal abundance and diversity in the Challenger Deep (at ca. 11,000 m depth) and the slope of the Mariana Trench in comparison with three sites of the adjacent abyssal plain. Our results indicate that trench sediments are a hotspot of fungal abundance in terms of the 18S rRNA gene copy number. The fungal diversity (as the number of amplicon sequence variants, ASVs) was relatively low at all sites (10-31 ASVs) but showed a high turnover diversity among stations due to the presence of exclusive fungal taxa belonging to Aspergillaceae, Trichosphaeriaceae, and Nectriaceae. Fungal abundance and diversity were closely linked to sediment organic matter content and composition (i.e., phytopigments and carbohydrates), suggesting a specialization of different fungal taxa for the exploitation of available resources. Overall, these findings provide new insights into the diversity of deep-sea fungi and the potential ecological role in trench sediments and pave the way for a better understanding of their relevance in one of the most extreme ecosystems on Earth.

Microeukaryotic plankton evolutionary constraints in a subtropical river explained by environment and bacteria along differing taxonomic resolutions.

Microeukaryotic plankton communities are keystone components for keeping aquatic primary productivity. Currently, variations in microeukaryotic plankton diversity have often been explained by local ecological factors but not by evolutionary constraints. We used amplicon sequencing of 100 water samples across five years to investigate the ecological preferences of the microeukaryotic plankton community in a subtropical riverine ecosystem. We found that microeukaryotic plankton diversity was less associated with bacterial abundance (16S rRNA gene copy number) than bacterial diversity. Further, environmental effects exhibited a larger influence on microeukaryotic plankton community composition than bacterial community composition, especially at fine taxonomic levels. The evolutionary constraints of microeukaryotic plankton community increased with decreasing taxonomic resolution (from 97% to 91% similarity levels), but not significant change from 85% to 70% similarity levels. However, compared with the bacterial community, the evolutionary constraints were shown to be more affected by environmental variables. This study illustrated possible controlling environmental and bacterial drivers of microeukaryotic diversity and community assembly in a subtropical river, thereby indirectly reflecting on the quality status of the water environment by providing new clues on the microeukaryotic community assembly.

Degradation, altered microbial community composition, and protein expression in bacterial consortium/fungus inoculated crude oil contaminated loamy soil

Accidental discharge of crude oil due to rupturing of old oil pipelines and blasts in oil refineries leads to contamination of agricultural land/adjacent areas and changes in microbial community and protein expression in the contaminated soil. Our previous study identified a bacterial consortium (Bacillus amyloliquefaciens MW532755+ Pseudomonas sp. MW444887) and a fungus (Aspergillus sydowii KY614299) with crude oil degrading ability. These microorganisms were used in the present investigation to degrade total petroleum hydrocarbon in naturally contaminated loamy soil (25% w/w) collected from the Assam region of India. The half-life of crude oil in uninoculated (T1), bacterial consortium (T2), and fungus-inoculated (T3) soils were 135.9, 47.5, and 58.2 days, respectively. The bacterial consortium and the fungus treatments altered the soil bacterial community and metabolic pathways, evidently from the analyses using metagenomics and predicted functional proteomics. The sequencing analysis using Ion Torrent showed the presence of 287 genera in the oil-polluted and treated soil samples. At the Genus level, abundances of Pseudomonas and Bacilli were observed, and the reads representing those of the Domain Bacteria for uninoculated (T1), bacterial consortium (T2), and fungus (T3) treated soils were 214525, 584054 and 489423, respectively. The predicted functional proteomics analysis for the protein-regulating genes, functional pathways, and protein families showed the presence of diversified pathways responsible for higher metabolic activities (q < 0.05). The quantitative polymerase chain reaction (qPCR) analysis showed that, compared to fungus treatment, bacterial consortium increased the abundance of 16S rRNA gene copies of Archaea, Alphaproteobacteria, and Bacteroidetes and also the gene copies of nitrogen cycling such as nifH (106 to 108) were highest, followed by nirK (104 to 108), amoA (104 to 105), and nirS (104 to 106)The present study suggested that the bacterial consortium as the hydrocarbon degraders was better than fungus and can be used in the future for oil remediation purposes.

Biosolids-Derived Biochar Improves Biomethane Production in the Anaerobic Digestion of Chicken Manure

Anaerobic digestion has attracted great interest for use in the management of organic wastes and the production of biomethane. However, this process is facing challenges, such as a high concentration of ammonia nitrogen, which affects the methanogenesis process and, thus, the production of methane. This study investigates the use of biosolid-derived biochar for mitigating ammonia stress and improving methane production during the anaerobic digestion of chicken manure, using both pristine biochar and biochar modified with a potassium hydroxide (KOH) solution. Batch mesophilic anaerobic digestion (37 °C) was carried out over 18 days. When compared to chicken-manure-only controls, a significant increase in methane formation was observed in the digesters amended with biochar and KOH-modified biochar, producing 220 L kg−1 volatile solids (VSs) and 262 L kg−1 VSs of methane, respectively, compared to 139 L kg−1 VSs from the control digesters. The use of biochar and KOH-modified biochar resulted in a significant reduction of 8 days in the lag phase. The total ammonia nitrogen (TAN) concentration was reduced in the digesters with biochar and KOH-modified biochar by 25% and 35.5%, respectively. The quantitative polymerase chain reaction (QPCR) data revealed that the number of 16S rRNA gene copies was around 50,000 and 41,000 times higher in the biochar and KOH-modified biochar digesters, respectively, compared to the control digesters on day 18. The taxonomic profiles indicated that the BC and KOH-BC digesters contained a mixture of methanogenic pathways, including acetoclastic (Methanosaetaceae), hydrogenotrophic (Methanosarcinaceae), and methylation (Methanofastidiosaceae). This mix of pathways suggests a more robust archaeal community and, hence, more efficient methanogenesis. The results show that the addition of biosolids biochar enhances anaerobic digestion, mitigates ammonia stress to methanogens, and significantly increases biogas production.

Cultivating Resilience in Dryland Soils: An Assisted Migration Approach to Biological Soil Crust Restoration.

Land use practices and climate change have driven substantial soil degradation across global drylands, impacting ecosystem functions and human livelihoods. Biological soil crusts, a common feature of dryland ecosystems, are under extensive exploration for their potential to restore the stability and fertility of degraded soils through the development of inoculants. However, stressful abiotic conditions often result in the failure of inoculation-based restoration in the field and may hinder the long-term success of biocrust restoration efforts. Taking an assisted migration approach, we cultivated biocrust inocula sourced from multiple hot-adapted sites (Mojave and Sonoran Deserts) in an outdoor facility at a cool desert site (Colorado Plateau). In addition to cultivating inoculum from each site, we created an inoculum mixture of biocrust from the Mojave Desert, Sonoran Desert, and Colorado Plateau. We then applied two habitat amelioration treatments to the cultivation site (growth substrate and shading) to enhance soil stability and water availability and reduce UV stress. Using marker gene sequencing, we found that the cultivated mixed inoculum comprised both local- and hot-adapted cyanobacteria at the end of cultivation but had similar cyanobacterial richness as each unmixed inoculum. All cultivated inocula had more cyanobacterial 16S rRNA gene copies and higher cyanobacterial richness when cultivated with a growth substrate and shade. Our work shows that it is possible to field cultivate biocrust inocula sourced from different deserts, but that community composition shifts toward that of the cultivation site unless habitat amelioration is employed. Future assessments of the function of a mixed inoculum in restoration and its resilience in the face of abiotic stressors are needed to determine the relative benefit of assisted migration compared to the challenges and risks of this approach.

Achieving Discharge Limits in Single-Stage Domestic Wastewater Treatment by Combining Urban Waste Sources and Phototrophic Mixed Cultures.

This work shows the potential of a new way of co-treatment of domestic wastewater (DWW) and a liquid stream coming from the thermal hydrolysis of the organic fraction of municipal solid waste (OFMSW) mediated by a mixed culture of purple phototrophic bacteria (PPB) capable of assimilating carbon and nutrients from the medium. The biological system is an open single-step process operated under microaerophilic conditions at an oxidative reduction potential (ORP) < 0 mV with a photoperiod of 12/24 h and fed during the light stage only so the results can be extrapolated to outdoor open pond operations by monitoring the ORP. The effluent mostly complies with the discharge values of the Spanish legislation in COD and p-values (<125 mg/L; <2 mg/L), respectively, and punctually on values in N (<15 mg/L). Applying an HRT of 3 d and a ratio of 100:7 (COD:N), the presence of PPB in the mixed culture surpassed 50% of 16S rRNA gene copies, removing 78% of COD, 53% of N, and 66% of P. Furthermore, by increasing the HRT to 5 d, removal efficiencies of 83% of COD, 65% of N, and 91% of P were achieved. In addition, the reactors were further operated in a membrane bioreactor, thus separating the HRT from the SRT to increase the specific loading rate. Very satisfactory removal efficiencies were achieved by applying an HRT and SRT of 2.3 and 3 d, respectively: 84% of COD, 49% of N, and 93% of P despite the low presence of PPB due to more oxidative conditions, which step-by-step re-colonized the mixed culture until reaching >20% of 16S rRNA gene copies after 49 d of operation. These results open the door to scaling up the process in open photobioreactors capable of treating urban wastewater and municipal solid waste in a single stage and under microaerophilic conditions by controlling the ORP of the system.

Hyperarid soil microbial community response to simulated rainfall.

The exceptionally long and protracted aridity in the Atacama Desert (AD), Chile, provides an extreme, terrestrial ecosystem that is ideal for studying microbial community dynamics under hyperarid conditions. Our aim was to characterize the temporal response of hyperarid soil AD microbial communities to ex situ simulated rainfall (5% g water/g dry soil for 4 weeks) without nutrient amendment. We conducted replicated microcosm experiments with surface soils from two previously well-characterized AD hyperarid locations near Yungay at 1242 and 1609 masl (YUN1242 and YUN1609) with distinct microbial community compositions and average soil relative humidity levels of 21 and 17%, respectively. The bacterial and archaeal response to soil wetting was evaluated by 16S rRNA gene qPCR, and amplicon sequencing. Initial YUN1242 bacterial and archaeal 16S rRNA gene copy numbers were significantly higher than for YUN1609. Over the next 4 weeks, qPCR results showed significant increases in viable bacterial abundance, whereas archaeal abundance decreased. Both communities were dominated by 10 prokaryotic phyla (Actinobacteriota, Proteobacteria, Chloroflexota, Gemmatimonadota, Firmicutes, Bacteroidota, Planctomycetota, Nitrospirota, Cyanobacteriota, and Crenarchaeota) but there were significant site differences in the relative abundances of Gemmatimonadota and Chloroflexota, and specific actinobacterial orders. The response to simulated rainfall was distinct for the two communities. The actinobacterial taxa in the YUN1242 community showed rapid changes while the same taxa in the YUN1609 community remained relatively stable until day 30. Analysis of inferred function of the YUN1242 microbiome response implied an increase in the relative abundance of known spore-forming taxa with the capacity for mixotrophy at the expense of more oligotrophic taxa, whereas the YUN1609 community retained a stable profile of oligotrophic, facultative chemolithoautotrophic and mixotrophic taxa. These results indicate that bacterial communities in extreme hyperarid soils have the capacity for growth in response to simulated rainfall; however, historic variations in long-term hyperaridity exposure produce communities with distinct putative metabolic capacities.

The resuscitation-promoting factor (Rpf) from Micrococcus luteus and its putative reaction product 1,6-anhydro-MurNAc increase culturability of environmental bacteria.

Dormant bacterial cells do not divide and are not immediately culturable, but they persist in a state of low metabolic activity, a physiological state having clinical relevance, for instance in latent tuberculosis. Resuscitation-promoting factors (Rpfs) are proteins that act as signalling molecules mediating growth and replication. In this study we aimed to test the effect of Rpfs from Micrococcus luteus on the number and diversity of cultured bacteria using insect and soil samples, and to examine if the increase in culturability could be reproduced with the putative reaction product of Rpf, 1,6-anhydro-N-acetylmuramic acid (1,6-anhydro-MurNAc). The rpf gene from Micrococcus luteus was amplified and cloned into a pET21b expression vector and the protein was expressed in Escherichia coli BL21(DE3) cells and purified by affinity chromatography using a hexa-histidine tag. 1,6-Anhydro-MurNAc was prepared using reported chemical synthesis methods. Recombinant Rpf protein or 1,6-anhydro-MurNAc were added to R2A cultivation media, and their effect on the culturability of bacteria from eight environmental samples including four cockroach guts and four soils was examined. Colony-forming units, 16S rRNA gene copies and Illumina amplicon sequencing of the 16S rRNA gene were measured for all eight samples subjected to three different treatments: Rpf, 1,6-anhydro-MurNAc or blank control. Both Rpf and 1,6-anhydro-MurNAc increased the number of colony-forming units and of 16S rRNA gene copies across the samples although the protein was more effective. The Rpf and 1,6-anhydro-MurNAc promoted the cultivation of a diverse set of bacteria and in particular certain clades of the phyla Actinomycetota and Bacillota . This study opens the path for improved cultivation strategies aiming to isolate and study yet undescribed living bacterial organisms.

Validation of 16S rRNA gene sequencing and metagenomics for evaluating microbial immigration in a methanogenic bioreactor.

To quantitatively evaluate the impact of microbial immigration from an upstream community on the microbial assembly of a downstream community, an ecological genomics (ecogenomics)-based mass balance (EGMB) model coupled with 16S rRNA gene sequencing was previously developed. In this study, a mock community was used to further validate the EGMB models and demonstrate the feasibility of using metagenome-based EGMB model to reveal both microbial activity and function. The mock community consisting of Aeromonas, Escherichia, and Pseudomonas was fed into a lab-scale methanogenic bioreactor together with dissolved organic substrate. Using qPCR, 16S rRNA gene, 16S rRNA gene copy number normalization (GCN), and metagenome, results showed highly comparable community profiles in the feed. In the bioreactor, Aeromonas and Pseudomonas exhibited negative growth rates throughout the experiment by all approaches. Escherichia's growth rate was negative by most biomarkers but was slightly positive by 16S rRNA gene. Still, all approaches showed a decreasing trend toward negative in the growth rate of Escherichia as reactor operation time increased. Uncultivated populations of phyla Desulfobacterota, Chloroflexi, Actinobacteriota, and Spirochaetota were observed to increase in abundance, suggesting their contribution in degrading the feed biomass. Based on metabolic reconstruction of metagenomes, these populations possessed functions of hydrolysis, fermentation, fatty acid degradation, or acetate oxidation. Overall results supported the application of both 16S rRNA gene- and metagenome-based EGMB models to measure the growth rate of microbes in the bioreactor, and the latter had advantage in providing insights into the microbial functions of uncultivated populations.

Concatenated 16S rRNA sequence analysis improves bacterial taxonomy.

Background: Microscopic, biochemical, molecular, and computer-based approaches are extensively used to identify and classify bacterial populations. Advances in DNA sequencing and bioinformatics workflows have facilitated sophisticated genome-based methods for microbial taxonomy although sequencing of the 16S rRNA gene is widely employed to identify and classify bacterial communities as a cost-effective and single-gene approach. However, the 16S rRNA sequence-based species identification accuracy is limited because of the occurrence of multiple copies of the 16S rRNA gene and higher sequence identity between closely related species. The availability of the genomes of several bacterial species provided an opportunity to develop comprehensive species-specific 16S rRNA reference libraries. Methods: Sequences of the 16S rRNA genes were retrieved from the whole genomes available in the Genome databases. With defined criteria, four 16S rRNA gene copy variants were concatenated to develop a species-specific reference library. The sequence similarity search was performed with a web-based BLAST program, and MEGA software was used to construct the phylogenetic tree. Results: Using this approach, species-specific 16S rRNA gene libraries were developed for four closely related Streptococcusspecies ( S. gordonii, S. mitis, S. oralis, and S. pneumoniae). Sequence similarity and phylogenetic analysis using concatenated 16S rRNA copies yielded better resolution than single gene copy approaches. Conclusions: The approach is very effective in classifying genetically closely related bacterial species and may reduce misclassification of bacterial species and genome assemblies.

The Effects of Long-Term Application of Stabilized and Coated Urea on Soil Chemical Properties, Microbial Community Structure, and Functional Genes in Paddy Fields

The addition of fertilizers is indispensable in agricultural production, and currently, there is a wide variety of new types of fertilizers available. For example, commonly used are stabilized fertilizers with inhibitors and coated slow-release fertilizers, among others. However, the long-term effects of these fertilizers, when applied continuously are still uncertain. This study will provide scientific and theoretical support for the development and promotion of these fertilizers. A 16-year paddy field with brown soil treated with different urease and nitrification inhibitors, sulfur-coated urea (SCU), and resin-coated urea (PCU) was studied. The study showed that long-term use of conventional urea nitrogen fertilizer resulted in a significant reduction in soil total phosphorus (TP). Long-term application of NBPT and conventional urea significantly increased soil organic matter (SOM). Moreover, except for HQ and NBPT+DMPP, the prolonged application of new urea fertilizers also significantly enhanced soil total potassium (TK). Application of SCU fertilizer in brown soil type paddy fields resulted in a significant decrease in soil pH over time. However, changes in pH had no effect on the abundance of ammonia-oxidizing bacteria (AOB), as AOB was mainly affected by soil-available N. DMPP, HQ+DCD, NBPT+DMPP, SCU, and PCU significantly reduced the 16S rRNA gene copy number of soil bacteria, with the greatest effect of coated urea fertilizer (SCU and PCU). Long-term application of stable urea fertilizer with HQ significantly reduced the bacterial community in paddy soil. Conversely, HQ+DCD-stabilizede urea fertilizer significantly increased the population structure and abundance of Basidiomycota fungi while decreasing the population structure and abundance of Rozellomycota fungi. DMPP-stabilized urea fertilizer notably increased the population structure and abundance of Ascomycota fungi while decreasing the population structure and abundance of Rozellomycota and Chytridiomycota fungi. Furthermore, HQ-stabilized urea fertilizer significantly reduced the population structure and abundance of Chytridiomycota fungi.

Responses of syntrophic microbial communities and their interactions with polystyrene nanoplastics in a microbial electrolysis cell

Microbial electrochemical technologies are promising for simultaneous energy recovery and wastewater treatment. Although the inhibitory effects of emerging pollutants, particularly micro/nanoplastics (MPs/NPs), on conventional wastewater systems have been extensively studied, the current understanding of their impact on microbial electrochemical systems is still quite limited. Microplastics are plastic particles ranging from 1 μm to 5 mm. However, nanoplastics are smaller plastic particles ranging from 1 to 100 nm. Due to their smaller size and greater surface area, they can penetrate deeper into biofilm structures and cell membranes, potentially disrupting their integrity and leading to changes in biofilm composition and function. This study first reports the impact of polystyrene nanoplastics (PsNPs) on syntrophic anode microbial communities in a microbial electrolysis cell. Low concentrations of PsNPs (50 and 250 μg/L) had a minimal impact on current density and hydrogen production. However, 500 μg/L of PsNPs decreased the maximum current density and specific hydrogen production rate by ∼43 % and ∼48 %, respectively. Exposure to PsNPs increased extracellular polymeric substance (EPS) levels, with a higher ratio of carbohydrates to proteins, suggesting a potential defense mechanism through EPS secretion. The downregulation of genes associated with extracellular electron transfer was observed at 500 μg/L of PsNPs. Furthermore, the detrimental impact of 500 μg/L PsNPs on the microbiome was evident from the decrease in 16S rRNA gene copies, microbial diversity, richness, and relative abundances of key electroactive and fermentative bacteria. For the first time, this study presents the inhibitory threshold of any NPs on syntrophic electroactive biofilms within a microbial electrochemical system.

Assembly and co-occurrence patterns of rhizosphere bacterial communities are closely linked to soil fertility during continuous cropping of cut chrysanthemum (Chrysanthemum morifolium Ramat).

Continuous cropping is known to have profound effects on the soil microbial community in different planting systems. However, we lack an understanding of how different years of continuous cropping affects rhizosphere soil bacterial community co-occurrence pattern and assembly processes in the cut chrysanthemum (Chrysanthemum morifolium Ramat.) field. We collected the soils from cut chrysanthemum rhizospheres with planting for 1 year (PY1) and continuous cropping for 6 years (CY6) and 12 years (CY12). Real-time quantitative PCR and flow cytometry (FCM) techniques were used to test the 16S rRNA gene copy number and bacterial cell count, respectively. The bacterial community structure was analysed by using high-throughput sequencing technology. The CY12 had a significantly decreased soil fertility index and rhizosphere bacterial living cell counts and gene copy numbers compared to CY6 and PY1 (P <0.05). The rhizosphere bacterial community dissimilarity increased as the continuous cropping years increased. Three main ecological clusters (modules #1, #2, and #3) were observed in the bacterial co-occurrence network across all samples, and only the relative abundance of module #1 (enriched in the CY12) was significantly correlated with soil fertility (P <0.05). Moreover, the rhizosphere bacterial community assembly was primarily governed by the deterministic process under 12 years of continuous cropping. Soil fertility decline correlates with ecological network modularization and the deterministic assembly process of the rhizosphere bacterial community of cut chrysanthemum during continuous cropping.

Different in root exudates and rhizosphere microorganisms effect on nitrogen removal between three emergent aquatic plants in surface flow constructed wetlands

Swine wastewater contains high concentration of nitrogen (N), causing pollution of surrounding water bodies. Constructed wetlands (CWs) are considered as an effective ecological treatment measure to remove nitrogen. Some emergent aquatic plants could tolerate high ammonia, and play a crucial part in CWs to treat high concentration N wastewater. However, the mechanism of root exudates and rhizosphere microorganisms of emergent plants on nitrogen removal is still unclear. Effects of organic and amino acids on rhizosphere N cycle microorganisms and environmental factors across three emergent plants were investigated in this study. The highest TN removal efficiency were 81.20% in surface flow constructed wetlands (SFCWs) plant with Pontederia cordata. The root exudation rates results showed that organic and amino acids were higher in 56 d than that in 0 d in SFCWs plants with Iris pseudacorus and P. cordata. The highest ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) gene copy numbers were found in I. pseudacorus rhizosphere soil, while the highest nirS, nirK, hzsB and 16S rRNA gene copy numbers were detected in P. cordata rhizosphere soil. Regression analysis results demonstrated that organic and amino acids exudation rates were positive related to rhizosphere microorganisms. These results indicated that organic and amino acids secretion could stimulate growth of emergent plants rhizosphere microorganisms in SFCWs for swine wastewater treatment. In addition, the EC, TN, NH4+-N and NO3−-N were negatively correlated with organic and amino acids exudation rates, and abundances of rhizosphere microorganisms via Pearson correlation analysis. These results imply that organic and amino acids, and rhizosphere microorganisms synergically affected on the nitrogen removal in SFCWs.

A novel Gardnerella, Prevotella, and Lactobacillus standard that improves accuracy in quantifying bacterial burden in vaginal microbial communities.

Bacterial vaginosis (BV) is the most common vaginal dysbiosis. In this condition, a polymicrobial biofilm develops on vaginal epithelial cells. Accurately quantifying the bacterial burden of the BV biofilm is necessary to further our understanding of BV pathogenesis. Historically, the standard for calculating total bacterial burden of the BV biofilm has been based on quantifying Escherichia coli 16S rRNA gene copy number. However, E. coli is improper for measuring the bacterial burden of this unique micro-environment. Here, we propose a novel qPCR standard to quantify bacterial burden in vaginal microbial communities, from an optimal state to a mature BV biofilm. These standards consist of different combinations of vaginal bacteria including three common BV-associated bacteria (BVAB) Gardnerella spp. (G), Prevotella spp. (P), and Fannyhessea spp. (F) and commensal Lactobacillus spp. (L) using the 16S rRNA gene (G:P:F:L, G:P:F, G:P:L and 1G:9L). We compared these standards to the traditional E. coli (E) reference standard using known quantities of mock vaginal communities and 16 vaginal samples from women. The E standard significantly underestimated the copy numbers of the mock communities, and this underestimation was significantly greater at lower copy numbers of these communities. The G:P:L standard was the most accurate across all mock communities and when compared to other mixed vaginal standards. Mixed vaginal standards were further validated with vaginal samples. This new G:P:L standard can be used in BV pathogenesis research to enhance reproducibility and reliability in quantitative measurements of BVAB, spanning from the optimal to non-optimal (including BV) vaginal microbiota.

Accounting for 16S rRNA copy number prediction uncertainty and its implications in bacterial diversity analyses

16S rRNA gene copy number (16S GCN) varies among bacterial species and this variation introduces potential biases to microbial diversity analyses using 16S rRNA read counts. To correct the biases, methods have been developed to predict 16S GCN. A recent study suggests that the prediction uncertainty can be so great that copy number correction is not justified in practice. Here we develop RasperGade16S, a novel method and software to better model and capture the inherent uncertainty in 16S GCN prediction. RasperGade16S implements a maximum likelihood framework of pulsed evolution model and explicitly accounts for intraspecific GCN variation and heterogeneous GCN evolution rates among species. Using cross-validation, we show that our method provides robust confidence estimates for the GCN predictions and outperforms other methods in both precision and recall. We have predicted GCN for 592605 OTUs in the SILVA database and tested 113842 bacterial communities that represent an exhaustive and diverse list of engineered and natural environments. We found that the prediction uncertainty is small enough for 99% of the communities that 16S GCN correction should improve their compositional and functional profiles estimated using 16S rRNA reads. On the other hand, we found that GCN variation has limited impacts on beta-diversity analyses such as PCoA, NMDS, PERMANOVA and random-forest test.

The microbiome of Riccia liverworts is an important reservoir for microbial diversity in temporary agricultural crusts

BackgroundThe microbiota of liverworts provides an interesting model for plant symbioses; however, their microbiome assembly is not yet understood. Here, we assessed specific factors that shape microbial communities associated with Riccia temporary agricultural crusts in harvested fields by investigating bacterial, fungal and archaeal communities in thalli and adhering soil from different field sites in Styria and Burgenland, Austria combining qPCR analyses, amplicon sequencing and advanced microscopy.ResultsRiccia spec. div. was colonized by a very high abundance of bacteria (1010 16S rRNA gene copies per g of thallus) as well as archaea and fungi (108 ITS copies per g of thallus). Each Riccia thallus contain approx. 1000 prokaryotic and fungal ASVs. The field type was the main driver for the enrichment of fungal taxa, likely due to an imprint on soil microbiomes by the cultivated crop plants. This was shown by a higher fungal richness and different fungal community compositions comparing liverwort samples collected from pumpkin fields, with those from corn fields. In contrast, bacterial communities linked to liverworts are highly specialized and the soil attached to them is not a significant source of these bacteria. Specifically, enriched Cyanobacteria, Bacteroidetes and Methylobacteria suggest a symbiotic interaction. Intriguingly, compared to the surrounding soil, the thallus samples were shown to enrich several well-known bacterial and fungal phytopathogens indicating an undescribed role of liverworts as potential reservoirs of crop pathogens.ConclusionsOur results provide evidence that a stable bacterial community but varying fungal communities are colonizing liverwort thalli. Post-harvest, temporary agricultural biocrusts are important reservoirs for microbial biodiversity but they have to be considered as potential reservoirs for pathogens as well.

Host-diet-gut microbiome interactions influence human energy balance: a randomized clinical trial

The gut microbiome is emerging as a key modulator of human energy balance. Prior studies in humans lacked the environmental and dietary controls and precision required to quantitatively evaluate the contributions of the gut microbiome. Using a Microbiome Enhancer Diet (MBD) designed to deliver more dietary substrates to the colon and therefore modulate the gut microbiome, we quantified microbial and host contributions to human energy balance in a controlled feeding study with a randomized crossover design in young, healthy, weight stable males and females (NCT02939703). In a metabolic ward where the environment was strictly controlled, we measured energy intake, energy expenditure, and energy output (fecal and urinary). The primary endpoint was the within-participant difference in host metabolizable energy between experimental conditions [Control, Western Diet (WD) vs. MBD]. The secondary endpoints were enteroendocrine hormones, hunger/satiety, and food intake. Here we show that, compared to the WD, the MBD leads to an additional 116 ± 56 kcals (P < 0.0001) lost in feces daily and thus, lower metabolizable energy for the host (89.5 ± 0.73%; range 84.2-96.1% on the MBD vs. 95.4 ± 0.21%; range 94.1-97.0% on the WD; P < 0.0001) without changes in energy expenditure, hunger/satiety or food intake (P > 0.05). Microbial 16S rRNA gene copy number (a surrogate of biomass) increases (P < 0.0001), beta-diversity changes (whole genome shotgun sequencing; P = 0.02), and fermentation products increase (P < 0.01) on an MBD as compared to a WD along with significant changes in the host enteroendocrine system (P < 0.0001). The substantial interindividual variability in metabolizable energy on the MBD is explained in part by fecal SCFAs and biomass. Our results reveal the complex host-diet-microbiome interplay that modulates energy balance.

Comparison of two methods for bioaerosol sampling and characterization in a low-biomass chamber environment

Bioaerosols are emitted from various sources into the indoor environment and can positively and negatively impact human health. Humans are the major source of bioaerosol emissions indoors, specifically for bacteria. However, efficient sampling to guarantee successful downstream analyses can be challenging due to the relatively low bioaerosol concentrations in many indoor spaces and variable susceptibility of bioaerosols to sampling stress. Establishing standard procedures for collecting bacteria in low biomass indoor environments can help advance the field. We compared the performance of two approaches (sampling with a personal environmental monitor (PEM) and a two-stage dry cyclone sampler) to capture the bacterial emission from human participants in a controlled chamber environment. The comparison was based on quantifying the DNA yield and characterizing the bacterial community diversity by metabarcoding of the 16S rRNA gene. We found notable differences in the performance of the samplers. The cyclone sampler collected significantly more DNA and 16S rRNA gene copies including Gram-negative bacteria than the PEM sampler (p < 0.001). The bacterial barcode sequencing revealed a significant difference in the diversity of bacteria captured by the two samplers (p < 0.05), with a higher diversity of Gram-negative bacteria captured by PEM than the cyclone sampler. Overall, this study showed that both sampling approaches efficiently sampled in low biomass indoor environments, however with substantial differences in captured DNA concentrations, bacterial concentrations, and bacterial diversity. Our results indicate that bioaerosol sampler choice is highly research question dependent, and this study provides data support to make informed choices.

Variable Region Sequences Influence 16S rRNA Performance.

16S rRNA gene sequences are commonly analyzed for taxonomic and phylogenetic studies because they contain variable regions that can help distinguish different genera. However, intra-genus distinction using variable region homology is often impossible due to the high overall sequence identities among closely related species, even though some residues may be conserved within respective species. Using a computational method that included the allelic diversity within individual genomes, we discovered that certain Escherichia and Shigella species can be distinguished by a multi-allelic 16S rRNA variable region single nucleotide polymorphism (SNP). To evaluate the performance of 16S rRNAs with altered variable regions, we developed an in vivo system that measures the acceptance and distribution of variant 16S rRNAs into a large pool of natural versions supporting normal translation and growth. We found that 16S rRNAs containing evolutionarily disparate variable regions were underpopulated both in ribosomes and in active translation pools, even for an SNP. Overall, this study revealed that variable region sequences can substantially influence the performance of 16S rRNAs and that this biological constraint can be leveraged to justify refining taxonomic assignments of variable region sequence data. IMPORTANCE This study reevaluates the notion that 16S rRNA gene variable region sequences are uninformative for intra-genus classification and that single nucleotide variations within them have no consequence to strains that bear them. We demonstrated that the performance of 16S rRNAs in Escherichia coli can be negatively impacted by sequence changes in variable regions, even for single nucleotide changes that are native to closely related Escherichia and Shigella species; thus, biological performance is likely constraining the evolution of variable regions in bacteria. Further, the native nucleotide variations we tested occur in all strains of their respective species and across their multiple 16S rRNA gene copies, suggesting that these species evolved beyond what would be discerned from a consensus sequence comparison. Therefore, this work also reveals that the multiple 16S rRNA gene alleles found in most bacteria can provide more informative phylogenetic and taxonomic detail than a single reference allele.