Publications

Microorganisms are critical in mediating carbon (C) and nitrogen (N) cycling processes in soils. Yet, it has long been debated whether the processes underlying biogeochemical cycles are affected by the composition and diversity of the soil microbial community or not. The composition and diversity of soil microbial communities can be influenced by various environmental factors, which in turn are known to impact biogeochemical processes. The objectives of this study were to test effects of multiple edaphic drivers individually and represented as the multivariate soil environment interacting with microbial community composition and diversity, and concomitantly on multiple soil functions (i.e. soil enzyme activities, soil C and N processes). We employed high-throughput sequencing (Illumina MiSeq) to analyze bacterial/archaeal and fungal community composition by targeting the 16S rRNA gene and the ITS1 region of soils collected from three land uses (cropland, grassland and forest) deriving from two bedrock forms (silicate and limestone). Based on this data set we explored single and combined effects of edaphic variables on soil microbial community structure and diversity, as well as on soil enzyme activities and several soil C and N processes. We found that both bacterial/archaeal and fungal communities were shaped by the same edaphic factors, with most single edaphic variables and the combined soil environment representation exerting stronger effects on bacterial/archaeal communities than on fungal communities, as demonstrated by (partial) Mantel tests. We also found similar edaphic controls on the bacterial/archaeal/fungal richness and diversity. Soil C processes were only directly affected by the soil environment but not affected by microbial community composition. In contrast, soil N processes were significantly related to bacterial/archaeal community composition and bacterial/archaeal/fungal richness/diversity but not directly affected by the soil environment. This indicates direct control of the soil environment on soil C processes and indirect control of the soil environment on soil N processes by structuring the microbial communities. The study further highlights the importance of edaphic drivers and microbial communities (i.e. composition and diversity) on important soil C and N processes.

Diazotrophic microorganisms introduce biologically available nitrogen (N) to the global N cycle through the activity of the nitrogenase enzyme. The genetically conserved dinitrogenase reductase (nifH) gene is phylogenetically distributed across four clusters (I-IV) and is widely used as a marker gene for N₂ fixation, permitting investigators to study the genetic diversity of diazotrophs in nature and target potential participants in N₂ fixation. To date there have been limited, standardized pipelines for the nifH functional gene, which is in stark contrast to the rRNA gene. Here we present a bioinformatics pipeline for processing nifH amplicon datasets – NifMAP (“NifH MiSeq Illumina amplicon Analysis Pipeline”), which as a novel aspect uses Hidden-Markov models to filter out homologous genes to nifH. By using this pipeline, we evaluated the broadly inclusive primer pairs (Ueda19F-R6, IGK3-DVV, F2-R6) that target the nifH gene. To evaluate any systematic biases, the nifH gene was amplified with the aforementioned primer pairs in a diverse collection of environmental samples (soils, rhizosphere and roots samples, biological soil crusts and estuarine samples), in addition to a nifH mock community consisting of six phylogenetically diverse members. We noted that all primer pairs co-amplified nifH homologs to varying degrees; up to 90% of the amplicons were nifH homologs with IGK3-DVV in some samples (rhizosphere and roots from tall oat-grass). In regards to specificity, we observed some degree of bias across the primer pairs. For example, primer pair F2-R6 discriminated against cyanobacteria (amongst others), yet captured many sequences from subclusters IIIE and IIIL-N. These aforementioned subclusters were largely missing by the primer pair IGK3-DVV, which also tended to discriminate against Alphaproteobacteria, but amplified sequences within clusters IIIC (affiliated with Clostridia) and clusters IVB and IVC. Primer pair Ueda19F-R6 exhibited the least bias and successfully captured diazotrophs in cluster I and subclusters IIIE, IIIL, IIIM and IIIN, but discriminated against Firmicutes and subcluster IIIC. Taken together, our newly established bioinformatics pipeline, NifMAP, along with our systematic evaluations of nifH primer pairs permit more robust, high-throughput investigations of diazotrophs in diverse environments. 

Sulfur-cycling microorganisms impact organic matter decomposition in wetlands and consequently greenhouse gas emissions from these globally relevant environments. However, their identities and physiological properties are largely unknown. By applying a functional metagenomics approach to an acidic peatland, we recovered draft genomes of seven novel Acidobacteria species with the potential for dissimilatory sulfite (dsrAB, dsrC, dsrD, dsrN, dsrT, dsrMKJOP) or sulfate respiration (sat, aprBA, qmoABC plus dsr genes). Surprisingly, the genomes also encoded DsrL, which so far was only found in sulfur-oxidizing microorganisms. Metatranscriptome analysis demonstrated expression of acidobacterial sulfur-metabolism genes in native peat soil and their upregulation in diverse anoxic microcosms. This indicated an active sulfate respiration pathway, which, however, might also operate in reverse for dissimilatory sulfur oxidation or disproportionation as proposed for the sulfur-oxidizing Desulfurivibrio alkaliphilus. Acidobacteria that only harbored genes for sulfite reduction additionally encoded enzymes that liberate sulfite from organosulfonates, which suggested organic sulfur compounds as complementary energy sources. Further metabolic potentials included polysaccharide hydrolysis and sugar utilization, aerobic respiration, several fermentative capabilities, and hydrogen oxidation. Our findings extend both, the known physiological and genetic properties of Acidobacteria and the known taxonomic diversity of microorganisms with a DsrAB-based sulfur metabolism, and highlight new fundamental niches for facultative anaerobic Acidobacteria in wetlands based on exploitation of inorganic and organic sulfur molecules for energy conservation.

Members of the phylum Acidobacteria are abundant and ubiquitous across soils. We performed the largest (to date) comparative genome analysis spanning subdivisions 1, 3, 4, 6, 8, and 23 (n=24) with the goal to identify features to help explain their prevalence in soils and understand their ecophysiology. In contrast to earlier studies, our analysis revealed that bacteriophage integration events along with transposable and mobile elements influenced the structure and plasticity of these genomes. Low- and high-affinity respiratory oxygen reductases were detected in multiple genomes, suggesting the capacity for growing across different oxygen gradients. Amongst many genomes, the capacity to use a diverse collection of carbohydrates, as well as inorganic and organic N sources (such as extracellular peptidases), were detected – both advantageous traits in environments with fluctuating nutrient environments. We also identified multiple soil acidobacteria with the potential to scavenge atmospheric concentrations of H₂, now encompassing mesophilic soil strains within the subdivision 1 and 3, in addition to a previously identified thermophilic strain in subdivision 4. This large-scale acidobacteria genome analysis reveals traits that provide genomic, physiological and metabolic versatility, presumably allowing flexibility and versatility in the challenging and fluctuating soil environment.

Investigating active participants in the fixation of dinitrogen gas is vital as N is often a limiting factor for primary production. Biological nitrogen fixation (BNF) is performed by a diverse guild of bacteria and archaea (diazotrophs), which can be free-living or symbionts. Free-living diazotrophs are widely distributed in the environment, yet our knowledge about their identity and ecophysiology is still limited. A major challenge in investigating this guild is inferring activity from genetic data as this process is highly regulated. To address this challenge, we evaluated and improved several 15N-based methods for detecting N2 fixation activity (with a focus on soil samples) and studying active diazotrophs. We compared the acetylene reduction assay and the 15N2 tracer method and demonstrated that the latter is more sensitive in samples with low activity. Additionally, tracing 15N into microbial RNA provides much higher sensitivity compared to bulk soil analysis. Active soil diazotrophs were identified with a 15N-RNA-SIP approach optimized for environmental samples and benchmarked to 15N-DNA-SIP. Lastly, we investigated the feasibility of using SIP-Raman microspectroscopy for detecting 15N-labelled cells. Taken together, these tools allow identifying and investigating active free-living diazotrophs in a highly sensitive manner in diverse environments, from bulk to the single-cell level.

The combined approach of incubating environmental samples with stable isotope-labeled substrates followed by single-cell analyses through high-resolution secondary ion mass spectrometry (NanoSIMS) or Raman microspectroscopy provides insights into the in situ function of microorganisms. This approach has found limited application in soils presumably due to the dispersal of microbial cells in a large background of particles. We developed a pipeline for the efficient preparation of cell extracts from soils for subsequent single-cell methods by combining cell detachment with separation of cells and soil particles followed by cell concentration. The procedure was evaluated by examining its influence on cell recoveries and microbial community composition across two soils. This approach generated a cell fraction with considerably reduced soil particle load and of sufficient small size to allow single-cell analysis by NanoSIMS, as shown when detecting active N₂-fixing and cellulose-responsive microorganisms via ¹⁵N₂ and ¹³C-UL-cellulose incubations, respectively. The same procedure was also applicable for Raman microspectroscopic analyses of soil microorganisms, assessed via microcosm incubations with a ¹³C-labeled carbon source and deuterium oxide (D₂O, a general activity marker). The described sample preparation procedure enables single-cell analysis of soil microorganisms using NanoSIMS and Raman microspectroscopy, but should also facilitate single-cell sorting and sequencing.

Biological nitrogen fixation is the primary supply of N to most ecosystems, yet there is considerable uncertainty about how N-fixing bacteria will respond to global change factors such as increasing atmospheric CO₂ and N deposition. Using the nifH gene as a molecular marker, we studied how the community structure of N-fixing soil bacteria from temperate pine, aspen, and sweet gum stands and a brackish tidal marsh responded to multiyear elevated CO₂conditions. We also examined how N availability, specifically, N fertilization, interacted with elevated CO₂ to affect these communities in the temperate pine forest. Based on data from Sanger sequencing and quantitative PCR, the soil nifHcomposition in the three forest systems was dominated by species in the Geobacteraceae and, to a lesser extent, Alphaproteobacteria. The N-fixing-bacterial-community structure was subtly altered after 10 or more years of elevated atmospheric CO₂, and the observed shifts differed in each biome. In the pine forest, N fertilization had a stronger effect on nifH community structure than elevated CO₂ and suppressed the diversity and abundance of N-fixing bacteria under elevated atmospheric CO₂ conditions. These results indicate that N-fixing bacteria have complex, interacting responses that will be important for understanding ecosystem productivity in a changing climate.

Microbial communities that deconstruct plant biomass have broad relevance in biofuel production and global carbon cycling. Biomass pretreatments reduce plant biomass recalcitrance for increased efficiency of enzymatic hydrolysis. We exploited these chemical pretreatments to study howthermophilic bacterial consortia adapt to deconstruct switchgrass (SG) biomass of varying compositions. Microbial communities were adapted to untreated, ammonium fiber expansion (AFEX)-pretreated, and ionic liquid (IL)-pretreated SG under aerobic, thermophilic conditions using green waste compost as the inoculum to study biomass deconstruction by microbial consortia. After microbial cultivation, gravimetric analysis of the residual biomass demonstrated that both AFEX- and IL-pretreatment enhanced the deconstruction of the SG biomass by approximately 2-fold. 2D-NMR experiments and acetyl bromide-reactive lignin analysis indicated that polysaccharide hydrolysis was the dominant process occurring during microbial biomass deconstruction and lignin remaining in the residual biomass was largely unmodified. SSU rRNA gene amplicon libraries revealed that although the dominant taxa across these chemical pretreatments were consistently represented by members of the Firmicutes, Bacteroidetes, and Deinococcus-Thermus phyla, the abundance of select OTUs varied suggesting adaptations to the different substrates. Combining the observations of differences in the community structure and the chemical and physical structure of the biomass, we hypothesize specific roles for individual community members in biomass deconstruction.

Enzymatic hydrolysis of cellulose is a key process in the global carbon cycle and the industrial conversion of biomass to biofuels. In natural environments, cellulose hydrolysis is predominately performed by microbial communities. However, detailed understanding of bacterial cellulose hydrolysis is primarily confined to a few model isolates. Developing models for cellulose hydrolysis by mixed microbial consortia will complement these isolate studies and may reveal new mechanisms for cellulose deconstruction. Microbial communities were adapted to microcrystalline cellulose under aerobic, thermophilic conditions using green waste compost as the inoculum to study cellulose hydrolysis in a microbial consortium. This adaptation selected for three dominant taxa – the Firmicutes, Bacteroidetes and Thermus. A high-resolution profile of community development during the enrichment demonstrated a community transition from Firmicutes to a novel Bacteroidetes population that clusters in the Chitinophagaceae family. A representative strain of this population, strain NYFB, was successfully isolated, and sequencing of a nearly full-length 16S rRNA gene demonstrated that it was only 86% identical compared with other validated strains in the phylum Bacteroidetes. Strain NYFB grew well on soluble polysaccharide substrates, but grew poorly on insoluble polysaccharide substrates. Similar communities were observed in companion thermophilic enrichments on insoluble wheat arabinoxylan, a hemicellulosic substrate, suggesting a common model for deconstruction of plant polysaccharides. Combining observations of community dynamics and the physiology of strain NYFB, a cooperative successional model for polysaccharide hydrolysis by the Firmicutes and Bacteroidetes in the thermophilic cellulolytic consortia is proposed.

Taxonomic and phylogenetic fingerprinting based on sequence analysis of gene fragments from the large-subunit rRNA (LSU) gene or the internal transcribed spacer (ITS) region is becoming an integral part of fungal classification. The lack of an accurate and robust classification tool trained by a validated sequence database for taxonomic placement of fungal LSU genes is a severe limitation in taxonomic analysis of fungal isolates or large data sets obtained from environmental surveys. Using a hand-curated set of 8,506 fungal LSU gene fragments, we determined the performance characteristics of a naïve Bayesian classifier across multiple taxonomic levels and compared the classifier performance to that of a sequence similarity-based (BLASTN) approach. The naïve Bayesian classifier was computationally more rapid (>460-fold with our system) than the BLASTN approach, and it provided equal or superior classification accuracy. Classifier accuracies were compared using sequence fragments of 100 bp and 400 bp and two different PCR primer anchor points to mimic sequence read lengths commonly obtained using current high-throughput sequencing technologies. Accuracy was higher with 400-bp sequence reads than with 100-bp reads. It was also significantly affected by sequence location across the 1,400-bp test region. The highest accuracy was obtained across either the D1 or D2 variable region. The naïve Bayesian classifier provides an effective and rapid means to classify fungal LSU sequences from large environmental surveys. The training set and tool are publicly available through the Ribosomal Database Project (http://rdp.cme.msu.edu/classifier/classifier.jsp).

Many bacteria and fungi are known to degrade cellulose in culture, but their combined response to cellulose in different soils is unknown. Replicate soil microcosms amended with [¹³C]cellulose were used to identify bacterial and fungal communities responsive to cellulose in five geographically and edaphically different soils. The diversity and composition of the cellulose-responsive communities were assessed by DNA-stable isotope probing combined with Sanger sequencing of small-subunit and large-subunit rRNA genes for the bacterial and fungal communities, respectively. In each soil, the ¹³C-enriched, cellulose-responsive communities were of distinct composition compared to the original soil community or ¹²C-nonenriched communities. The composition of cellulose-responsive taxa, as identified by sequence operational taxonomic unit (OTU) similarity, differed in each soil. When OTUs were grouped at the bacterial order level, we found that members of the Burkholderiales, Caulobacteriales, Rhizobiales, Sphingobacteriales, Xanthomonadales, and the subdivision 1 Acidobacteria were prevalent in the ¹³C-enriched DNA in at least three of the soils. The cellulose-responsive fungi were identified as members of the Trichocladium, Chaetomium, Dactylaria, and Arthrobotrys genera, along with two novel Ascomycota clusters, unique to one soil. Although similarities were identified in higher-level taxa among some soils, the composition of cellulose-responsive bacteria and fungi was generally unique to a certain soil type, suggesting a strong potential influence of multiple edaphic factors in shaping the community.

Six terrestrial ecosystems in the USA were exposed to elevated atmospheric CO2 in single or multifactorial experiments for more than a decade to assess potential impacts. We retrospectively assessed soil bacterial community responses in all six-field experiments and found ecosystem-specific and common patterns of soil bacterial community response to elevated CO2. Soil bacterial composition differed greatly across the six ecosystems. No common effect of elevated atmospheric CO2 on bacterial biomass, richness and community composition across all of the ecosystems was identified, although significant responses were detected in individual ecosystems. The most striking common trend across the sites was a decrease of up to 3.5-fold in the relative abundance of Acidobacteria Group 1 bacteria in soils exposed to elevated CO2 or other climate factors. The Acidobacteria Group 1 response observed in exploratory 16S rRNA gene clone library surveys was validated in one ecosystem by 100-fold deeper sequencing and semi-quantitative PCR assays. Collectively, the 16S rRNA gene sequencing approach revealed influences of elevated CO2 on multiple ecosystems. Although few common trends across the ecosystems were detected in the small surveys, the trends may be harbingers of more substantive changes in less abundant, more sensitive taxa that can only be detected by deeper surveys.

Bacteria modulate glycoside hydrolase expression in response to the changes in the composition of lignocellulosic biomass. The response of switchgrass-adapted thermophilic bacterial consortia to perturbation with a variety of biomass substrates was characterized to determine if bacterial consortia also responded to changes in biomass composition. Incubation of the switchgrass-adapted consortia with these alternative substrates produced shifts in glycoside hydrolase activities and bacterial community composition. Substantially increased endoglucanase activity was observed upon incubation with microcrystalline cellulose and trifluororacetic acid-pretreated switchgrass. In contrast, culturing the microbial consortia with ionic liquid-pretreated switchgrass increased xylanase activity dramatically. Microbial community analyses of these cultures indicated that the increased endoglucanase activity correlated with an increase in bacteria related to Rhodothermus marinus. Inclusion of simple organic substrates in the culture medium abrogated glycoside hydrolase activity and enriched for bacteria related to Thermus thermophilus. These results demonstrate that the composition of biomass substrates influences the glycoside hydrolase activities and community composition of biomass-deconstructing bacterial consortia.

Members of the phylum Acidobacteria are among the most abundant bacteria in soil. Although they have been characterized as versatile heterotrophs, it is unclear if the types and availability of organic resources influence their distribution in soil. The potential for organic resources to select for different acidobacteria was assessed using molecular and cultivation-based approaches with agricultural and managed grassland soils in Michigan. The distribution of acidobacteria varied with the carbon content of soil: the proportion of subdivision 4 sequences was highest in agricultural soils (ca. 41%) that contained less carbon than grassland soils, where the proportions of subdivision 1, 3, 4, and 6 sequences were similar. Either readily oxidizable carbon or plant polymers were used as the sole carbon and energy source to isolate heterotrophic bacteria from these soils. Plant polymers increased the diversity of acidobacteria cultivated but decreased the total number of heterotrophs recovered compared to readily oxidizable carbon. Two phylogenetically novel Acidobacteria strains isolated on the plant polymer medium were characterized. Strains KBS 83 (subdivision 1) and KBS 96 (subdivision 3) are moderate acidophiles with pH optima of 5.0 and 6.0, respectively. Both strains grew slowly (μ = 0.01 h⁻¹) and harbored either 1 (strain KBS 83) or 2 (strain KBS 96) copies of the 16S rRNA encoding gene—a genomic characteristic typical of oligotrophs. Strain KBS 83 is a microaerophile, growing optimally at 8% oxygen. These metabolic characteristics help delineate the niches that acidobacteria occupy in soil and are consistent with their widespread distribution and abundance.

Members of the bacterial phylum Acidobacteria are widespread in soils and sediments worldwide, and are abundant in many soils. Acidobacteria are challenging to culture in vitro,and many basic features of their biology and functional roles in the soil have not been determined. Candidatus Solibacter usitatus strain Ellin6076 has a 9.9 Mb genome that is approximately 2–5 times as large as the other sequenced Acidobacteria genomes. Bacterial genome sizes typically range from 0.5 to 10 Mb and are influenced by gene duplication, horizontal gene transfer, gene loss and other evolutionary processes. Our comparative genome analyses indicate that the Ellin6076 large genome has arisen by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, and widespread small-scale gene duplications, resulting in an increased number of paralogs. Low amino acid sequence identities among functional group members, and lack of conserved gene order and orientation in regions containing similar groups of paralogs, suggest that most of the paralogs are not the result of recent duplication events. The genome sizes of additional cultured Acidobacteria strains were estimated using pulsed-field gel electrophoresis to determine the prevalence of the large genome trait within the phylum. Members of subdivision 3 had larger genomes than those of subdivision 1, but none were as large as the Ellin6076 genome. The large genome of Ellin6076 may not be typical of the phylum, and encodes traits that could provide a selective metabolic, defensive and regulatory advantage in the soil environment.

Bacteria in the phylum Acidobacteria are widely distributed and abundant in soils, but their ecological roles are poorly understood, owing in part to a paucity of cultured representatives. In a molecular survey of acidobacterial diversity at the Michigan State University Kellogg Biological Station Long-Term Ecological Research site, 27% of acidobacterial 16S rRNA gene clones in a never-tilled, successional plant community belonged to subdivision 1, whose relative abundance varied inversely with soil pH. Strains of subdivision 1 were isolated from these never-tilled soils using low-nutrient medium incubated for 3 to 4 weeks under elevated levels of carbon dioxide, which resulted in a slightly acidified medium that matched the pH optima of the strains (between 5 and 6). Colonies were approximately 1 mm in diameter and either white or pink, the latter due to a carotenoid(s) that was synthesized preferentially under 20% instead of 2% oxygen. Strains were gram-negative, aerobic, chemo-organotrophic, nonmotile rods that produced an extracellular matrix. All strains contained either one or two copies of the 16S rRNA encoding gene, which along with a relatively slow doubling time (10 to 15 h at ca. 23°C) is suggestive of an oligotrophic lifestyle. Six of the strains are sufficiently similar to one another, but distinct from previously named Acidobacteria, to warrant creation of a new genus, Terriglobus, with Terriglobus roseus defined as the type species. The physiological and nutritional characteristics of Terriglobus are consistent with its potential widespread distribution in soil.

An integrative approach was used to obtain pure cultures of previously uncultivated members of the divisions Acidobacteria and Verrucomicrobia from agricultural soil and from the guts of wood-feeding termites. Some elements of the cultivation procedure included the following: the use of agar media with little or no added nutrients; relatively long periods of incubation (more than 30 days); protection of cells from exogenous peroxides; and inclusion of humic acids or a humic acid analogue (anthraquinone disulfonate) and quorum-signaling compounds (acyl homoserine lactones) in growth media. The bacteria were incubated in the presence of air and in hypoxic (1 to 2% O₂ [vol/vol]) and anoxic atmospheres. Some bacteria were incubated with elevated concentrations of CO₂(5% [vol/vol]). Significantly more Acidobacteria were found on isolation plates that had been incubated with 5% CO₂. A simple, high-throughput, PCR-based surveillance method (plate wash PCR) was developed. This method greatly facilitated detection and ultimate isolation of target bacteria from as many as 1,000 colonies of nontarget microbes growing on the same agar plates. Results illustrate the power of integrating culture methods with molecular techniques to isolate bacteria from phylogenetic groups underrepresented in culture.