Genomic characterization and computational phenotyping of nitrogen-fixing bacteria isolated from Colombian sugarcane fields.

Previous studies have shown the sugarcane microbiome harbors diverse plant growth promoting microorganisms, including nitrogen-fixing bacteria (diazotrophs), which can serve as biofertilizers. The genomes of 22 diazotrophs from Colombian sugarcane fields were sequenced to investigate potential biofertilizers. A genome-enabled computational phenotyping approach was developed to prioritize sugarcane associated diazotrophs according to their potential as biofertilizers. This method selects isolates that have potential for nitrogen fixation and other plant growth promoting (PGP) phenotypes while showing low risk for virulence and antibiotic resistance. Intact nitrogenase (nif) genes and operons were found in 18 of the isolates. Isolates also encode phosphate solubilization and siderophore production operons, and other PGP genes. The majority of sugarcane isolates showed uniformly low predicted virulence and antibiotic resistance compared to clinical isolates. Six strains with the highest overall genotype scores were experimentally evaluated for nitrogen fixation, phosphate solubilization, and the production of siderophores, gibberellic acid, and indole acetic acid. Results from the biochemical assays were consistent and validated computational phenotype predictions. A genotypic and phenotypic threshold was observed that separated strains by their potential for PGP versus predicted pathogenicity. Our results indicate that computational phenotyping is a promising tool for the assessment of bacteria detected in agricultural ecosystems.

Microbiota-directed fibre activates both targeted and secondary metabolic shifts in the distal gut.

Beneficial modulation of the gut microbiome has high-impact implications not only in humans, but also in livestock that sustain our current societal needs. In this context, we have tailored an acetylated galactoglucomannan (AcGGM) fibre to match unique enzymatic capabilities of Roseburia and Faecalibacterium species, both renowned butyrate-producing gut commensals. Here, we test the accuracy of AcGGM within the complex endogenous gut microbiome of pigs, wherein we resolve 355 metagenome-assembled genomes together with quantitative metaproteomes. In AcGGM-fed pigs, both target populations differentially express AcGGM-specific polysaccharide utilization loci, including novel, mannan-specific esterases that are critical to its deconstruction. However, AcGGM-inclusion also manifests a "butterfly effect", whereby numerous metabolic changes and interdependent cross-feeding pathways occur in neighboring non-mannanolytic populations that produce short-chain fatty acids. Our findings show how intricate structural features and acetylation patterns of dietary fibre can be customized to specific bacterial populations, with potential to create greater modulatory effects at large.

Isolation and Characterization of Extracellular Vesicles Secreted In Vitro by Porcine Microbiota.

The secretion of extracellular vesicles, EVs, is a common process in both prokaryotic and eukaryotic cells for intercellular communication, survival, and pathogenesis. Previous studies have illustrated the presence of EVs in supernatants from pure cultures of bacteria, including Gram-positive and Gram-negative glycan-degrading gut commensals. However, the isolation and characterization of EVs secreted by a complex microbial community have not been clearly reported. In a recent paper, we showed that wood-derived, complex ß-mannan, which shares a structural similarity with conventional dietary fibers, can be used to modulate the porcine gut microbiota composition and activity. In this paper, we investigated the production, size, composition, and proteome of EVs secreted by pig fecal microbiota after 24 h enrichment on complex ß-mannan. Using transmission electron microscopy and nanoparticle tracking analysis, we identified EVs with an average size of 165 nm. We utilized mass spectrometry-based metaproteomic profiling of EV proteins against a database of 355 metagenome-assembled genomes (MAGs) from the porcine colon and thereby identified 303 proteins. For EVs isolated from the culture grown on ß-mannan, most proteins mapped to two MAGs, MAG53 and MAG272, belonging to the orders Clostridiales and Bacilli, respectively. Furthermore, the MAG with the third-most-detected protein was MAG 343, belonging to the order Enterobacteriales. The most abundant proteins detected in the ß-mannan EVs proteome were involved in translation, energy production, amino acid, and carbohydrate transport, as well as metabolism. Overall, this proof-of-concept study demonstrates the successful isolation of EVs released from a complex microbial community; furthermore, the protein content of the EVs reflects the response of specific microbes to the available carbohydrate source.

Effects of H2:CO2 ratio and H2 supply fluctuation on methane content and microbial community composition during in-situ biological biogas upgrading.

BACKGROUND: Commercial biogas upgrading facilities are expensive and consume energy. Biological biogas upgrading may serve as a low-cost approach because it can be easily integrated with existing facilities at biogas plants. The microbial communities found in anaerobic digesters typically contain hydrogenotrophic methanogens, which can use hydrogen (H2) as a reducing agent for conversion of carbon dioxide (CO2) into methane (CH4). Thus, biological biogas upgrading through the exogenous addition of H2 into biogas digesters for the conversion of CO2 into CH4 can increase CH4 yield and lower CO2 emission. RESULTS: The addition of 4 mol of H2 per mol of CO2 was optimal for batch biogas reactors and increased the CH4 content of the biogas from 67 to 94%. The CO2 content of the biogas was reduced from 33 to 3% and the average residual H2 content was 3%. At molar H2:CO2 ratios > 4:1, all CO2 was converted into CH4, but the pH increased above 8 due to depletion of CO2, which negatively influenced the process stability. Additionally, high residual H2 content in these reactors was unfavourable, causing volatile fatty acid accumulation and reduced CH4 yields. The reactor microbial communities shifted in composition over time, which corresponded to changes in the reactor variables. Numerous taxa responded to the H2 inputs, and in particular the hydrogenotrophic methanogen Methanobacterium increased in abundance with addition of H2. In addition, the apparent rapid response of hydrogenotrophic methanogens to intermittent H2 feeding indicates the suitability of biological methanation for variable H2 inputs, aligning well with fluctuations in renewable electricity production that may be used to produce H2. CONCLUSIONS: Our research demonstrates that the H2:CO2 ratio has a significant effect on reactor performance during in situ biological methanation. Consequently, the H2:CO2 molar ratio should be kept at 4:1 to avoid process instability. A shift toward hydrogenotrophic methanogenesis was indicated by an increase in the abundance of the obligate hydrogenotrophic methanogen Methanobacterium.

Feeding frequency influences process performance and microbial community composition in anaerobic digesters treating steam exploded food waste.

In anaerobic digestion, studies of feeding frequency have produced conflicting results. Hence, the effect of feeding frequency on process variables and microbial community structure was investigated by comparing a laboratory-scale digester fed steam exploded food waste 10 times daily vs. one fed an equivalent amount once daily. The Frequently Fed Digester (FFD) produced on average 20% more methane and had lower effluent concentrations of long-chain fatty acids. Greater daily fluctuations in acetate, pH and biogas production rate could explain the lower specific methane yield and ß-oxidation. Feeding frequency also influenced the microbial community whereby Tenericutes (42%) dominated in FFD but Firmicutes (31%) was most abundant in the Daily Fed Digester (DFD). Feeding frequency effects are therefore postulated to occur more often in digesters fed labile feedstocks at high organic loading rates.

Characterization of the fecal and mucosa-associated microbiota in dogs with colorectal epithelial tumors.

Colorectal epithelial tumors occur spontaneously in dogs, and the pathogenesis seems to parallel that of humans. The development of human colorectal tumorigenesis has been linked to alterations in the composition of the intestinal microbiota. This study characterized the fecal- and mucosa-associated microbiota in dogs with colorectal epithelial tumors (n = 10). The fecal microbiota was characterized by 16S rDNA analysis and compared with that of control dogs (n = 13). We also determined the mucosa-associated microbiota composition in colonic tumor tissue (n = 8) and in adjacent non-tumor tissue (n = 5) by 16S rDNA- and rRNA profiling. The fecal microbial community structure in dogs with tumors was different from that of control samples and was distinguished by oligotypes affiliated with Enterobacteriaceae, Bacteroides, Helicobacter, Porphyromonas, Peptostreptococcus and Streptococcus, and lower abundance of Ruminococcaceae, Slackia, Clostridium XI and Faecalibacterium. The overall community structure and populations of mucosal bacteria were not different based on either the 16S rDNA or the 16S rRNA profile in tumor tissue vs. adjacent non-tumor tissue. However, the proportion of live, potentially active bacteria appeared to be higher in non-tumor tissue compared with tumor tissue and included Slackia, Roseburia, unclass. Ruminococcaeceae, unclass. Lachnospiraceae and Oscillibacter. Colorectal tumors are rarely diagnosed in dogs, but despite this limitation, we were able to show that dogs with colorectal tumors have distinct fecal microbiota profiles. These initial results support the need for future case-control studies that are adequately powered, as well as age-matched and breed-matched, in order to evaluate the influence of bacteria on colorectal cancer etiopathogenesis and to determine whether the bacteria may have potential as biomarkers in clinical settings.

Diazotroph Community Characterization via a High-Throughput nifH Amplicon Sequencing and Analysis Pipeline.

The dinitrogenase reductase gene (nifH) is the most widely established molecular marker for the study of nitrogen-fixing prokaryotes in nature. A large number of PCR primer sets have been developed for nifH amplification, and the effective deployment of these approaches should be guided by a rapid, easy-to-use analysis protocol. Bioinformatic analysis of marker gene sequences also requires considerable expertise. In this study, we advance the state of the art for nifH analysis by evaluating nifH primer set performance, developing an improved amplicon sequencing workflow, and implementing a user-friendly bioinformatics pipeline. The developed amplicon sequencing workflow is a three-stage PCR-based approach that uses established technologies for incorporating sample-specific barcode sequences and sequencing adapters. Based on our primer evaluation, we recommend the Ando primer set be used with a modified annealing temperature of 58°C, as this approach captured the largest diversity of nifH templates, including paralog cluster IV/V sequences. To improve nifH sequence analysis, we developed a computational pipeline which infers taxonomy and optionally filters out paralog sequences. In addition, we employed an empirical model to derive optimal operational taxonomic unit (OTU) cutoffs for the nifH gene at the species, genus, and family levels. A comprehensive workflow script named TaxADivA (TAXonomy Assignment and DIVersity Assessment) is provided to ease processing and analysis of nifH amplicons. Our approach is then validated through characterization of diazotroph communities across environmental gradients in beach sands impacted by the Deepwater Horizon oil spill in the Gulf of Mexico, in a peat moss-dominated wetland, and in various plant compartments of a sugarcane field.IMPORTANCE Nitrogen availability often limits ecosystem productivity, and nitrogen fixation, exclusive to prokaryotes, comprises a major source of nitrogen input that sustains food webs. The nifH gene, which codes for the iron protein of the nitrogenase enzyme, is the most widely established molecular marker for the study of nitrogen-fixing microorganisms (diazotrophs) in nature. In this study, a flexible sequencing/analysis pipeline, named TaxADivA, was developed for nifH amplicons produced by Illumina paired-end sequencing, and it enables an inference of taxonomy, performs clustering, and produces output in formats that may be used by programs that facilitate data exploration and analysis. Diazotroph diversity and community composition are linked to ecosystem functioning, and our results advance the phylogenetic characterization of diazotroph communities by providing empirically derived nifH similarity cutoffs for species, genus, and family levels. The utility of our pipeline is validated for diazotroph communities in a variety of ecosystems, including contaminated beach sands, peatland ecosystems, living plant tissues, and rhizosphere soil.

The effect of temperature and retention time on methane production and microbial community composition in staged anaerobic digesters fed with food waste.

BACKGROUND: Food waste is a large bio-resource that may be converted to biogas that can be used for heat and power production, or as transport fuel. We studied the anaerobic digestion of food waste in a staged digestion system consisting of separate acidogenic and methanogenic reactor vessels. Two anaerobic digestion parameters were investigated. First, we tested the effect of 55 vs. 65 °C acidogenic reactor temperature, and second, we examined the effect of reducing the hydraulic retention time (HRT) from 17 to 10 days in the methanogenic reactor. Process parameters including biogas production were monitored, and the microbial community composition was characterized by 16S amplicon sequencing. RESULTS: Neither organic matter removal nor methane production were significantly different for the 55 and 65 °C systems, despite the higher acetate and butyrate concentrations observed in the 65 °C acidogenic reactor. Ammonium levels in the methanogenic reactors were about 950 mg/L NH4+ when HRT was 17 days but were reduced to 550 mg/L NH4+ at 10 days HRT. Methane production increased from ~ 3600 mL/day to ~ 7800 when the HRT was decreased. Each reactor had unique environmental parameters and a correspondingly unique microbial community. In fact, the distinct values in each reactor for just two parameters, pH and ammonium concentration, recapitulate the separation seen in microbial community composition. The thermophilic and mesophilic digesters were particularly distinct from one another. The 55 °C acidogenic reactor was mainly dominated by Thermoanaerobacterium and Ruminococcus, whereas the 65 °C acidogenic reactor was initially dominated by Thermoanaerobacterium but later was overtaken by Coprothermobacter. The acidogenic reactors were lower in diversity (34-101 observed OTU0.97, 1.3-2.5 Shannon) compared to the methanogenic reactors (472-513 observed OTU0.97, 5.1-5.6 Shannon). The microbial communities in the acidogenic reactors were > 90% Firmicutes, and the Euryarchaeota were higher in relative abundance in the methanogenic reactors. CONCLUSIONS: The digestion systems had similar biogas production and COD removal rates, and hence differences in temperature, NH4+ concentration, and pH in the reactors resulted in distinct but similarly functioning microbial communities over this range of operating parameters. Consequently, one could reduce operational costs by lowering both the hydrolysis temperature from 65 to 55 °C and the HRT from 17 to 10 days.

The Use of Degenerate Primers in qPCR Analysis of Functional Genes Can Cause Dramatic Quantification Bias as Revealed by Investigation of nifH Primer Performance.

The measurement of functional gene abundance in diverse microbial communities often employs quantitative PCR (qPCR) with highly degenerate oligonucleotide primers. While degenerate PCR primers have been demonstrated to cause template-specific bias in PCR applications, the effect of such bias on qPCR has been less well explored. We used a set of diverse, full-length nifH gene standards to test the performance of several universal nifH primer sets in qPCR. We found significant template-specific bias in all but the PolF/PolR primer set. Template-specific bias caused more than 1000-fold mis-estimation of nifH gene copy number for three of the primer sets and one primer set resulted in more than 10,000-fold mis-estimation. Furthermore, such template-specific bias will cause qPCR estimates to vary in response to beta-diversity, thereby causing mis-estimation of changes in gene copy number. A reduction in bias was achieved by increasing the primer concentration. We conclude that degenerate primers should be evaluated across a range of templates, annealing temperatures, and primer concentrations to evaluate the potential for template-specific bias prior to their use in qPCR.

A comprehensive aligned nifH gene database: a multipurpose tool for studies of nitrogen-fixing bacteria.

We describe a nitrogenase gene sequence database that facilitates analysis of the evolution and ecology of nitrogen-fixing organisms. The database contains 32 954 aligned nitrogenase nifH sequences linked to phylogenetic trees and associated sequence metadata. The database includes 185 linked multigene entries including full-length nifH, nifD, nifK and 16S ribosomal RNA (rRNA) gene sequences. Evolutionary analyses enabled by the multigene entries support an ancient horizontal transfer of nitrogenase genes between Archaea and Bacteria and provide evidence that nifH has a different history of horizontal gene transfer from the nifDK enzyme core. Further analyses show that lineages in nitrogenase cluster I and cluster III have different rates of substitution within nifD, suggesting that nifD is under different selection pressure in these two lineages. Finally, we find that that the genetic divergence of nifH and 16S rRNA genes does not correlate well at sequence dissimilarity values used commonly to define microbial species, as stains having <3% sequence dissimilarity in their 16S rRNA genes can have up to 23% dissimilarity in nifH. The nifH database has a number of uses including phylogenetic and evolutionary analyses, the design and assessment of primers/probes and the evaluation of nitrogenase sequence diversity. Database URL: http://www.css.cornell.edu/faculty/buckley/nifh.htm.

A comprehensive evaluation of PCR primers to amplify the nifH gene of nitrogenase.

The nifH gene is the most widely sequenced marker gene used to identify nitrogen-fixing Bacteria and Archaea. Numerous PCR primers have been designed to amplify nifH, but a comprehensive evaluation of nifH PCR primers has not been performed. We performed an in silico analysis of the specificity and coverage of 51 universal and 35 group-specific nifH primers by using an aligned database of 23,847 nifH sequences. We found that there are 15 universal nifH primers that target 90% or more of nitrogen fixers, but that there are also 23 nifH primers that target less than 50% of nifH sequences. The nifH primers we evaluated vary in their phylogenetic bias and their ability to recover sequences from commonly sampled environments. In addition, many of these primers will amplify genes that do not mediate nitrogen fixation, and thus it would be advisable for researchers to screen their sequencing results for the presence of non-target genes before analysis. Universal primers that performed well in silico were tested empirically with soil samples and with genomic DNA from a phylogenetically diverse set of nitrogen-fixing strains. This analysis will be of great utility to those engaged in molecular analysis of nifH genes from isolates and environmental samples.

A global census of nitrogenase diversity.

The global diversity of nitrogen-fixing microorganisms was assessed through construction and analysis of an aligned database of 16,989 nifH sequences. We conclude that the diversity of diazotrophs is still poorly described and that many organisms remain to be discovered. Our analyses indicate that diversity is not distributed evenly across phylogenetic groups or across environments and that some of the most diverse assemblages and environments remain the most poorly characterized. The majority of OTUs were rare, falling in the long tail of the frequency distribution. The most dominant OTUs fell into either the Cyanobacteria or the α, ß, and γ Proteobacteria, and five of these dominant OTUs do not have any representatives cultivated in isolation. Soils contained the greatest diversity of nifH sequences of all of the environments surveyed. Cluster III, which is dominated by nifH sequences from obligate anaerobes, was found to contain the greatest diversity of all nifH lineages and is also the group for which diversity is the least sampled. Our findings provide context for ongoing efforts to explore diazotroph diversity, indicating specific groups and environments that remain poorly characterized.