Clay swelling is a colloidal phenomenon that has a large influence on flow and solute migration in soils and sediments. While models for clay swelling have been proposed over many years, debate remains as to the interaction forces that combine to produce the observed swelling behavior. Using cryogenic transmission electron microscopy (cryo-TEM) and small-angle X-ray scattering, we study the influence of salinity, in combination with layer charge, interlayer cation, and particle size, on montmorillonite swelling. We observe a decrease in swelling with increased layer charge, increased cation charge, and decreased cation hydration, each indicative of the critical influence of Coulombic attraction between the negatively charged layers and interlayer cations. Cryo-TEM images of individual montmorillonite particles also reveal that swelling is dependent upon the number of layers in a particle. Calculations of the van der Waals (vdW) interaction based on new measurements of Hamaker coefficients confirm that long-range vdW interactions extend beyond near-neighbor layer interactions and result in a decrease in layer spacing with a larger number of layers. This work clarifies the short- and long-range attractive interactions that govern clay structure and ultimately the stability and permeability of hydrated clays in the environment.

Remediation of industrial wastewater is important for preventing environmental contamination and enabling water reuse. Biological treatment for one industrial contaminant, thiocyanate (SCN–), relies upon microbial hydrolysis, but this process is sensitive to high loadings. To examine the activity and stability of a microbial community over increasing SCN– loadings, we established and operated a continuous-flow bioreactor fed increasing loadings of SCN–. A second reactor was fed ammonium sulfate to mimic breakdown products of SCN–. Biomass was sampled from both reactors for metagenomics and metaproteomics, yielding a set of genomes for 144 bacteria and one rotifer that constituted the abundant community in both reactors. We analyzed the metabolic potential and temporal dynamics of these organisms across the increasing loadings. In the SCN– reactor, Thiobacillus strains capable of SCN– degradation were highly abundant, whereas the ammonium sulfate reactor contained nitrifiers and heterotrophs capable of nitrate reduction. Key organisms in the SCN– reactor expressed proteins involved in SCN– degradation, sulfur oxidation, carbon fixation, and nitrogen removal. Lower performance at higher loadings was linked to changes in microbial community composition. This work provides an example of how meta-omics can increase our understanding of industrial wastewater treatment and inform iterative process design and development.

Metagenomic studies recently uncovered form II/III RubisCO genes, originally thought to only occur in archaea, from uncultivated bacteria of the candidate phyla radiation (CPR). There are no isolated CPR bacteria and these organisms are predicted to have limited metabolic capacities. Here we expand the known diversity of RubisCO from CPR lineages. We report a form of RubisCO, distantly similar to the archaeal form III RubisCO, in some CPR bacteria from the Parcubacteria (OD1), WS6 and Microgenomates (OP11) phyla. In addition, we significantly expand the Peregrinibacteria (PER) II/III RubisCO diversity and report the first II/III RubisCO sequences from the Microgenomates and WS6 phyla. To provide a metabolic context for these RubisCOs, we reconstructed near-complete (>93%) PER genomes and the first closed genome for a WS6 bacterium, for which we propose the phylum name Dojkabacteria. Genomic and bioinformatic analyses suggest that the CPR RubisCOs function in a nucleoside pathway similar to that proposed in Archaea. Detection of form II/III RubisCO and nucleoside metabolism gene transcripts from a PER supports the operation of this pathway in situ. We demonstrate that the PER form II/III RubisCO is catalytically active, fixing CO2 to physiologically complement phototrophic growth in a bacterial photoautotrophic RubisCO deletion strain. We propose that the identification of these RubisCOs across a radiation of obligately fermentative, small-celled organisms hints at a widespread, simple metabolic platform in which ribose may be a prominent currency.

Marine microbial communities experience daily fluctuations in light and temperature that can have important ramifications for carbon and nutrient cycling. Elucidation of such short time scale community-wide dynamics is hindered by system complexity. Hypersaline aquatic environments have lower species richness than marine environments and can be well-defined spatially, hence they provide a model system for diel cycle analysis. We conducted a 3-day time series experiment in a well-defined pool in hypersaline Lake Tyrrell, Australia. Microbial communities were tracked by combining cultivation-independent lipidomic, metagenomic and microscopy methods. The ratio of total bacterial to archaeal core lipids in the planktonic community increased by up to 58% during daylight hours and decreased by up to 32% overnight. However, total organism abundances remained relatively consistent over 3 days. Metagenomic analysis of the planktonic community composition, resolved at the genome level, showed dominance by Haloquadratum species and six uncultured members of the Halobacteriaceae. The post 0.8 μm filtrate contained six different nanohaloarchaeal types, three of which have not been identified previously, and cryo-transmission electron microscopy imaging confirmed the presence of small cells. Notably, these nano-sized archaea showed a strong diel cycle, with a pronounced increase in relative abundance over the night periods. We detected no eukaryotic algae or other photosynthetic primary producers, suggesting that carbon resources may derive from patchily distributed microbial mats at the sediment-water interface or from surrounding land. Results show the operation of a strong community-level diel cycle, probably driven by interconnected temperature, light abundance, dissolved oxygen concentration and nutrient flux effects.

A fundamental question in microbial ecology relates to community structure, and how this varies across environment types. It is widely believed that some environments, such as those at very low pH, host simple communities based on the low number of taxa, possibly due to the extreme environmental conditions. However, most analyses of species richness have relied on methods that provide relatively low ribosomal RNA (rRNA) sampling depth. Here we used community transcriptomics to analyze the microbial diversity of natural acid mine drainage biofilms from the Richmond Mine at Iron Mountain, California. Our analyses target deep pools of rRNA gene transcripts recovered from both natural and laboratory-grown biofilms across varying developmental stages. In all, 91.8% of the ~254 million Illumina reads mapped to rRNA genes represented in the SILVA database. Up to 159 different taxa, including Bacteria, Archaea and Eukaryotes, were identified. Diversity measures, ordination and hierarchical clustering separate environmental from laboratory-grown biofilms. In part, this is due to the much larger number of rare members in the environmental biofilms. Although Leptospirillum bacteria generally dominate biofilms, we detect a wide variety of other Nitrospira organisms present at very low abundance. Bacteria from the Chloroflexi phylum were also detected. The results indicate that the primary characteristic that has enabled prior extensive cultivation-independent ‘omic’ analyses is not simplicity but rather the high dominance by a few taxa. We conclude that a much larger variety of organisms than previously thought have adapted to this extreme environment, although only few are selected for at any one time.

Microbial communities that underpin global biogeochemical cycles will likely be influenced by elevated temperature associated with environmental change. Here, we test an approach to measure how elevated temperature impacts the physiology of individual microbial groups in a community context, using a model microbial-based ecosystem. The study is the first application of tandem mass tag (TMT)-based proteomics to a microbial community. We accurately, precisely and reproducibly quantified thousands of proteins in biofilms growing at 40, 43 and 46 °C. Elevated temperature led to upregulation of proteins involved in amino-acid metabolism at the level of individual organisms and the entire community. Proteins from related organisms differed in their relative abundance and functional responses to temperature. Elevated temperature repressed carbon fixation proteins from two Leptospirillum genotypes, whereas carbon fixation proteins were significantly upregulated at higher temperature by a third member of this genus. Leptospirillum group III bacteria may have been subject to viral stress at elevated temperature, which could lead to greater carbon turnover in the microbial food web through the release of viral lysate. Overall, these findings highlight the utility of proteomics-enabled community-based physiology studies, and provide a methodological framework for possible extension to additional mixed culture and environmental sample analyses.

Intimate interconnection of crystal growth, (oriented) aggregation and phase transformation seem common in the formation of nano- and microcrystalline materials from solutions. Yet, the mechanistic linkages between the different processes have not been fully understood. In this work, we studied the hydrothermal growth of akaganeite (β-FeOOH) nanorods and their transformation to micron-sized hematite (α-Fe2O3) spindles using high-resolution cryogenic transmission electron microscopy (cryo-TEM). Only akaganeite particles and hematite spindles were detected in the samples. Further, cryo-electron 3D tomograms show that akaganeite nanorods were aggregated into loose three-dimensional networks with some embedded hematite spindles. Based on our cryo-TEM and additional X-ray diffraction, electron microscopy, and chemical data, we propose the following mechanism: first, formation of the early-stage hematite spindles is driven by phase stability change due to increase in size caused by oriented aggregation of akaganeite. Then, akaganeite particles continue to transform to hematite upon contact with and recrystallization onto hematite surfaces, making hematite grow with a constant aspect ratio and forming micron-sized nano-porous single-crystal spindles. Our growth model interprets experimental observations well and it resolves previous long-time debate over whether the hematite spindles are formed via classical Ostwald ripening or by oriented aggregation of hematite nanoparticles. Possibly, this aggregation-based concurrent growth and transformation model may also be applicable to crystal growth and phase transformation in other systems.

Crystal growth is one of the most fundamental processes in nature. Understanding of crystal growth mechanisms has changed dramatically over the past two decades. One significant advance has been the recognition that growth does not only occur atom by atom, but often proceeds via attachment and fusion of either amorphous or crystalline particles. Results from recent experiments and calculations can be integrated to develop a simple, unified conceptual description of attachment-based crystal growth. This enables us to address three important questions: What are the driving forces for attachment-based growth? For crystalline particles, what enables the particles to achieve crystallographic coalignment? What determines the surface on which attachment occurs? We conclude that the extent of internal nanoparticle order controls the degree of periodicity and anisotropy in the surrounding electrostatic field. For crystalline particles, the orienting force stemming from the electrostatic field can promote oriented attachment events, although solvent-surface interactions modulate this control. In cases where perfect crystallographic alignment is not achieved, misorientation gives rise to structural defects that can fundamentally modify nanomaterial properties.Keywords: aggregation driving force; amorphous particle; electrostatic field; oriented aggregation; oriented attachment; particle orienting force; particle-mediated growth

Fermentation-based metabolism is an important ecosystem function often associated with environments rich in organic carbon, such as wetlands, sewage sludge and the mammalian gut. The diversity of microorganisms and pathways involved in carbon and hydrogen cycling in sediments and aquifers and the impacts of these processes on other biogeochemical cycles remain poorly understood. Here we used metagenomics and proteomics to characterize microbial communities sampled from an aquifer adjacent to the Colorado River at Rifle, CO, USA, and document interlinked microbial roles in geochemical cycling. The organic carbon content in the aquifer was elevated via acetate amendment of the groundwater occurring over 2 successive years. Samples were collected at three time points, with the objective of extensive genome recovery to enable metabolic reconstruction of the community. Fermentative community members include organisms from a new phylum, Melainabacteria, most closely related to Cyanobacteria, phylogenetically novel members of the Chloroflexi and Bacteroidales, as well as candidate phyla genomes (OD1, BD1-5, SR1, WWE3, ACD58, TM6, PER and OP11). These organisms have the capacity to produce hydrogen, acetate, formate, ethanol, butyrate and lactate, activities supported by proteomic data. The diversity and expression of hydrogenases suggests the importance of hydrogen metabolism in the subsurface. Our proteogenomic data further indicate the consumption of fermentation intermediates by Proteobacteria can be coupled to nitrate, sulfate and iron reduction. Thus, fermentation carried out by previously unknown members of sediment microbial communities may be an important driver of nitrogen, hydrogen, sulfur, carbon and iron cycling.

Spontaneous formation of curved nanorods is generally unexpected, since curvature introduces strain energy. However, electron microscopy shows that under hydrothermal conditions, some nanorods grown by oriented attachment of small anatase particles on {101} surfaces are curved and dislocation free. Molecular dynamics simulations show that the lattice energy of a curved anatase rod is actually lower than that of a linear rod due to more attractive long-range interatomic Coulombic interactions among atoms in the curved rod. The thermodynamic driving force stemming from lattice energy could be harnessed to produce asymmetric morphologies unexpected from classical Ostwald ripening with unusual shapes and properties.

The fractionation of hydrogen stable isotopes during lipid biosynthesis is larger in autotrophic than in heterotrophic microorganisms, possibly due to selective incorporation of hydrogen from water into NAD(P)H, resulting in D-depleted lipids. An analogous fractionation should occur during amino acid biosynthesis. Whereas these effects are traditionally measured using gas-phase isotope ratio on 1H-1H and 1H-2H, using an electrospray mass spectrometry-based technique on the original biomolecular structure and fitting of isotopic patterns we measured the hydrogen isotope compositions of proteins from an acidophilic microbial community with organism specificity and compared values with those for lipids. We showed that lipids were isotopically light by −260 ‰ relative to water in the growth solution; alternatively protein isotopic composition averaged −370 ‰. This difference suggests that steps in addition to NAD(P)H formation contribute to D/H fractionation. Further, autotrophic bacteria sharing 94% 16S rRNA gene sequence identity displayed statistically significant differences in protein hydrogen isotope fractionation, suggesting different metabolic traits consistent with distinct ecological niches or incorrectly annotated gene function. In addition, it was found that heterotrophic, archaeal members of the community had isotopically light protein (−323 ‰) relative to growth water and were significantly different from coexisting bacteria. This could be attributed to metabolite transfer from autotrophs and unknown aspects of fractionation associated with iron reduction. Differential fractionation of hydrogen stable isotopes into metabolites and proteins may reveal trophic levels of members of microbial communities. The approach developed here provided insights into the metabolic characteristics of organisms in natural communities and may be applied to analyze other systems.

Vanadium is a commercially important metal that is released into the environment by fossil fuel combustion and mining. Despite its prevalence as a contaminant, the potential for vanadium bioremediation has not been widely studied. Injection of acetate (as a carbon source) directly into an aquifer to biostimulate contaminated sediments in Colorado, United States, resulted in prolonged removal of aqueous vanadium for a period of at least two years. To further investigate this process, we simultaneously added acetate and vanadate (V5+) to columns that were packed with aquifer sediment and inserted into groundwater wells installed on the Colorado River floodplain. This allowed evaluation of the microbial response to amendments in columns that received an influx of natural groundwater. Our results demonstrate the removal of up to 99% of the added V5+(aq) and suggest microbial mediation. Most probable number measurements demonstrate up to a 50-fold increase in numbers of V5+-reducing cells in vanadium-amended columns compared to controls. 16S rRNA gene sequencing indicates decreased diversity and selection for specific taxa in columns that received vanadate compared to those that did not. Overall, our results demonstrate that acetate amendment can be an effective strategy for V removal, and that V bioremediation may be a viable technology.

Sediments are massive reservoirs of carbon compounds and host a large fraction of microbial life. Microorganisms within terrestrial aquifer sediments control buried organic carbon turnover, degrade organic contaminants, and impact drinking water quality. Recent 16S rRNA gene profiling indicates that members of the bacterial phylum Chloroflexi are common in sediment. Only the role of the class Dehalococcoidia, which degrade halogenated solvents, is well understood. Genomic sampling is available for only six of the approximate 30 Chloroflexi classes, so little is known about the phylogenetic distribution of reductive dehalogenation or about the broader metabolic characteristics of Chloroflexi in sediment.We used metagenomics to directly evaluate the metabolic potential and diversity of Chloroflexi in aquifer sediments. We sampled genomic sequence from 86 Chloroflexi representing 15 distinct lineages, including members of eight classes previously characterized only by 16S rRNA sequences. Unlike in the Dehalococcoidia, genes for organohalide respiration are rare within the Chloroflexi genomes sampled here. Near-complete genomes were reconstructed for three Chloroflexi. One, a member of an unsequenced lineage in the Anaerolinea, is an aerobe with the potential for respiring diverse carbon compounds. The others represent two genomically unsampled classes sibling to the Dehalococcoidia, and are anaerobes likely involved in sugar and plant-derived-compound degradation to acetate. Both fix CO2 via the Wood-Ljungdahl pathway, a pathway not previously documented in Chloroflexi. The genomes each encode unique traits apparently acquired from Archaea, including mechanisms of motility and ATP synthesis.Chloroflexi in the aquifer sediments are abundant and highly diverse. Genomic analyses provide new evolutionary boundaries for obligate organohalide respiration. We expand the potential roles of Chloroflexi in sediment carbon cycling beyond organohalide respiration to include respiration of sugars, fermentation, CO2 fixation, and acetogenesis with ATP formation by substrate-level phosphorylation.

The premature infant gut has low individual but high inter-individual microbial diversity compared with adults. Based on prior 16S rRNA gene surveys, many species from this environment are expected to be similar to those previously detected in the human microbiota. However, the level of genomic novelty and metabolic variation of strains found in the infant gut remains relatively unexplored.To study the stability and function of early microbial colonizers of the premature infant gut, nine stool samples were taken during the third week of life of a premature male infant delivered via Caesarean section. Metagenomic sequences were assembled and binned into near-complete and partial genomes, enabling strain-level genomic analysis of the microbial community.We reconstructed eleven near-complete and six partial bacterial genomes representative of the key members of the microbial community. Twelve of these genomes share >90% putative ortholog amino acid identity with reference genomes. Manual curation of the assembly of one particularly novel genome resulted in the first essentially complete genome sequence (in three pieces, the order of which could not be determined due to a repeat) for Varibaculum cambriense (strain Dora), a medically relevant species that has been implicated in abscess formation.During the period studied, the microbial community undergoes a compositional shift, in which obligate anaerobes (fermenters) overtake Escherichia coli as the most abundant species. Other species remain stable, probably due to their ability to either respire anaerobically or grow by fermentation, and their capacity to tolerate fluctuating levels of oxygen. Metabolic predictions for V. cambriense suggest that, like other members of the microbial community, this organism is able to process various sugar substrates and make use of multiple different electron acceptors during anaerobic respiration. Genome comparisons within the family Actinomycetaceae reveal important differences related to respiratory metabolism and motility.Genome-based analysis provided direct insight into strain-specific potential for anaerobic respiration and yielded the first genome for the genus Varibaculum. Importantly, comparison of these de novo assembled genomes with closely related isolate genomes supported the accuracy of the metagenomic methodology. Over a one-week period, the early gut microbial community transitioned to a community with a higher representation of obligate anaerobes, emphasizing both taxonomic and metabolic instability during colonization.

Stimulation of subsurface microorganisms to induce reductive immobilization of metals is a promising approach for bioremediation, yet the overall microbial community response is typically poorly understood. Here we used proteogenomics to test the hypothesis that excess input of acetate activates complex community functioning and syntrophic interactions among autotrophs and heterotrophs. A flow-through sediment column was incubated in a groundwater well of an acetate-amended aquifer and recovered during microbial sulfate reduction. De novo reconstruction of community sequences yielded near-complete genomes of Desulfobacter (Deltaproteobacteria), Sulfurovum- and Sulfurimonas-like Epsilonproteobacteria and Bacteroidetes. Partial genomes were obtained for Clostridiales (Firmicutes) and Desulfuromonadales-like Deltaproteobacteria. The majority of proteins identified by mass spectrometry corresponded to Desulfobacter-like species, and demonstrate the role of this organism in sulfate reduction (Dsr and APS), nitrogen fixation and acetate oxidation to CO2 during amendment. Results indicate less abundant Desulfuromonadales, and possibly Bacteroidetes, also actively contributed to CO2 production via the tricarboxylic acid (TCA) cycle. Proteomic data indicate that sulfide was partially re-oxidized by Epsilonproteobacteria through nitrate-dependent sulfide oxidation (using Nap, Nir, Nos, SQR and Sox), with CO2 fixed using the reverse TCA cycle. We infer that high acetate concentrations, aimed at stimulating anaerobic heterotrophy, led to the co-enrichment of, and carbon fixation in Epsilonproteobacteria. Results give an insight into ecosystem behavior following addition of simple organic carbon to the subsurface, and demonstrate a range of biological processes and community interactions were stimulated.

Iron-reducing bacteria (FeRB) play key roles in anaerobic metal and carbon cycling and carry out biogeochemical transformations that can be harnessed for environmental bioremediation. A subset of FeRB require direct contact with Fe(III)-bearing minerals for dissimilatory growth, yet these bacteria must move between mineral particles. Furthermore, they proliferate in planktonic consortia during biostimulation experiments. Thus, a key question is how such organisms can sustain growth under these conditions. Here we characterized planktonic microbial communities sampled from an aquifer in Rifle, Colorado, USA, close to the peak of iron reduction following in situ acetate amendment. Samples were cryo-plunged on site and subsequently examined using correlated two- and three-dimensional cryogenic transmission electron microscopy (cryo-TEM) and scanning transmission X-ray microscopy (STXM). The outer membranes of most cells were decorated with aggregates up to 150 nm in diameter composed of ~3 nm wide amorphous, Fe-rich nanoparticles. Fluorescent in situ hybridization of lineage-specific probes applied to rRNA of cells subsequently imaged via cryo-TEM identified Geobacter spp., a well-studied group of FeRB. STXM results at the Fe L2,3 absorption edges indicate that nanoparticle aggregates contain a variable mixture of Fe(II)–Fe(III), and are generally enriched in Fe(III). Geobacter bemidjiensis cultivated anaerobically in the laboratory on acetate and hydrous ferric oxyhydroxides also accumulated mixed-valence nanoparticle aggregates. In field-collected samples, FeRB with a wide variety of morphologies were associated with nano-aggregates, indicating that cell surface Fe(III) accumulation may be a general mechanism by which FeRB can grow while in planktonic suspension.

Natural microbial communities are extremely diverse and contain uncharacterized but functionally important small molecules. By coupling a deuterium (D) labeling technique to high mass accuracy untargeted liquid chromatography-electrospray ionization-mass spectrometry (LC–ESI–MS) metabolomic analysis, we found that natural acidophilic microbial biofilms dominated by bacteria of the genus Leptospirillum contained unusual lyso phosphatidylethanolamine (PE) lipids in high abundance (more than 10 nmol/mg of dry biomass). The unusual polar head group structure of these lipids is similar to lipids found in phylogenetically unrelated acidophilic chemoautolithotrophs and may be related to the affinity of these lipids for iron and calcium ions. Correlations of lyso phospholipid and proteome abundance patterns suggest a link between the lyso phospholipids and the UBA-type substrain of Leptospirillum group II. By combining untargeted metabolomics with D exchange we demonstrate the ability to identify cryptic but biologically functional small molecules in mixed microbial communities.

Novel nanomaterials properties can stem from unusual morphologies that arise via oriented attachment (OA)-based crystal growth. However, the currently understood thermodynamic driving force of surface energy reduction for OA cannot predict experimental observation of attachment on relatively low energy surfaces in some crystals or sequential OA events that result in nanocrystals with morphology not predicted by symmetry. In this work, using molecular energetic calculations, we show that orientation-specific long-range interatomic interactions, in addition to surface energy reduction, predict morphology development and explain how OA produces crystals with lower symmetry than the initial material. Results also show that Coulombic interactions, rather than van der Waals interactions, control OA of ionic nanocrystals. Our computational approach and results can guide the selection or design of nanomaterials to achieve desired morphology as well as new physical and chemical properties.Keywords: crystal growth; molecular simulation; morphology development; nanomaterial; oriented aggregation; oriented attachment; thermodynamic driving force;

The intestinal microbiome is a critical determinant of human health. Alterations in its composition have been correlated with chronic disorders, such as obesity and inflammatory bowel disease in adults, and may be associated with neonatal necrotizing enterocolitis in premature infants. Increasing evidence suggests that strain-level genomic variation may underpin distinct ecological trajectories within mixed populations, yet there have been few strain-resolved analyses of genotype–phenotype connections in the context of the human ecosystem. Here, we document strain-level genomic divergence during the first 3 wk of life within the fecal microbiota of an infant born at 28-wk gestation. We observed three compositional phases during colonization, and reconstructed and intensively curated population genomic datasets from the third phase. The relative abundance of two Citrobacter strains sharing ~99% nucleotide identity changed significantly over time within a community dominated by a nearly clonal Serratia population and harboring a lower abundance Enterococcus population and multiple plasmids and bacteriophage. Modeling of Citrobacter strain abundance suggests differences in growth rates and host colonization patterns. We identified genotypic variation potentially responsible for divergent strain ecologies, including hotspots of sequence variation in regulatory genes and intergenic regions, and in genes involved in transport, flagellar biosynthesis, substrate metabolism, and host colonization, as well as differences in the complements of these genes. Our results demonstrate that a community genomic approach can elucidate gut microbial colonization at the resolution required to discern medically relevant strain and species population dynamics, and hence improve our ability to diagnose and treat microbial community-mediated disorders.

Recovery of ribosomal small subunit genes by assembly of short read community DNA sequence data generally fails, making taxonomic characterization difficult. Here, we solve this problem with a novel iterative method, based on the expectation maximization algorithm, that reconstructs full-length small subunit gene sequences and provides estimates of relative taxon abundances. We apply the method to natural and simulated microbial communities, and correctly recover community structure from known and previously unreported rRNA gene sequences. An implementation of the method is freely available at https://github.com/csmiller/EMIRGE.

Extensive genomic characterization of multi-species acid mine drainage microbial consortia combined with laboratory cultivation has enabled the application of quantitative proteomic analyses at the community level. In this study, quantitative proteomic comparisons were used to functionally characterize laboratory-cultivated acidophilic communities sustained in pH 1.45 or 0.85 conditions. The distributions of all proteins identified for individual organisms indicated biases for either high or low pH, and suggests pH-specific niche partitioning for low abundance bacteria and archaea. Although the proteome of the dominant bacterium, Leptospirillum group II, was largely unaffected by pH treatments, analysis of functional categories indicated proteins involved in amino acid and nucleotide metabolism, as well as cell membrane/envelope biogenesis were overrepresented at high pH. Comparison of specific protein abundances indicates higher pH conditions favor Leptospirillum group III, whereas low pH conditions promote the growth of certain archaea. Thus, quantitative proteomic comparisons revealed distinct differences in community composition and metabolic function of individual organisms during different pH treatments. Proteomic analysis revealed other aspects of community function. Different numbers of phage proteins were identified across biological replicates, indicating stochastic spatial heterogeneity of phage outbreaks. Additionally, proteomic data were used to identify a previously unknown genotypic variant of Leptospirillum group II, an indication of selection for a specific Leptospirillum group II population in laboratory communities. Our results confirm the importance of pH and related geochemical factors in fine-tuning acidophilic microbial community structure and function at the species and strain level, and demonstrate the broad utility of proteomics in laboratory community studies.

Acid mine drainage (AMD), an extreme environment characterized by low pH and high metal concentrations, can support dense acidophilic microbial biofilm communities that rely on chemoautotrophic production based on iron oxidation. Field determined production rates indicate that, despite the extreme conditions, these communities are sufficiently well adapted to their habitats to achieve primary production rates comparable to those of microbial communities occurring in some non-extreme environments. To enable laboratory studies of growth, production and ecology of AMD microbial communities, a culturing system was designed to reproduce natural biofilms, including organisms recalcitrant to cultivation. A comprehensive metabolic labeling-based quantitative proteomic analysis was used to verify that natural and laboratory communities were comparable at the functional level. Results confirmed that the composition and core metabolic activities of laboratory-grown communities were similar to a natural community, including the presence of active, low abundance bacteria and archaea that have not yet been isolated. However, laboratory growth rates were slow compared with natural communities, and this correlated with increased abundance of stress response proteins for the dominant bacteria in laboratory communities. Modification of cultivation conditions reduced the abundance of stress response proteins and increased laboratory community growth rates. The research presented here represents the first description of the application of a metabolic labeling-based quantitative proteomic analysis at the community level and resulted in a model microbial community system ideal for testing physiological and ecological hypotheses.

Similar to virtually all components of natural environments, microbial systems are inherently complex and dynamic. Advances in cultivation-independent molecular methods have provided a route to study microbial consortia in their natural surroundings and to begin resolving the community structure, dominant metabolic processes and inter-organism interactions. However, the utility of these methods generally scales inversely with community complexity. By applying genomics-enabled methods to the study of natural microbial communities with reduced levels of species richness, a relatively comprehensive understanding of the metabolic networks and evolutionary processes within these communities can be attained. In such well-defined model systems, it is also possible to link emergent ecological patterns to their molecular and evolutionary underpinnings, facilitating construction of predictive ecosystem models. In this study, we review over a decade of research on one such system—acid mine drainage biofilm communities. We discuss the value and limitations of tractable model microbial communities in developing molecular methods for microbial ecology and in uncovering principles that may explain behavior in more complex systems.

Metagenomics has provided access to genomes of as yet uncultivated microorganisms in natural environments, yet there are gaps in our knowledge—particularly for Archaea—that occur at relatively low abundance and in extreme environments. Ultrasmall cells (<500 nm in diameter) from lineages without cultivated representatives that branch near the crenarchaeal/euryarchaeal divide have been detected in a variety of acidic ecosystems. We reconstructed composite, near-complete ~1-Mb genomes for three lineages, referred to as ARMAN (archaeal Richmond Mine acidophilic nanoorganisms), from environmental samples and a biofilm filtrate. Genes of two lineages are among the smallest yet described, enabling a 10% higher coding density than found genomes of the same size, and there are noncontiguous genes. No biological function could be inferred for up to 45% of genes and no more than 63% of the predicted proteins could be assigned to a revised set of archaeal clusters of orthologous groups. Some core metabolic genes are more common in Crenarchaeota than Euryarchaeota, up to 21% of genes have the highest sequence identity to bacterial genes, and 12 belong to clusters of orthologous groups that were previously exclusive to bacteria. A small subset of 3D cryo-electron tomographic reconstructions clearly show penetration of the ARMAN cell wall and cytoplasmic membranes by protuberances extended from cells of the archaeal order Thermoplasmatales. Interspecies interactions, the presence of a unique internal tubular organelle [Comolli, et al. (2009) ISME J 3:159–167], and many genes previously only affiliated with Crenarchaea or Bacteria indicate extensive unique physiology in organisms that branched close to the time that Cren- and Euryarchaeotal lineages diverged.

Bacterial species concepts are controversial. More widely accepted is the need to understand how differences in gene content and sequence lead to ecological divergence. To address this relationship in ecosystem context, we investigated links between genotype and ecology of two genotypic groups of Leptospirillum group II bacteria in comprehensively characterized, natural acidophilic biofilm communities. These groups share 99.7% 16S rRNA gene sequence identity and 95% average amino acid identity between their orthologs. One genotypic group predominates during early colonization, and the other group typically proliferates in later successional stages, forming distinct patches tens to hundreds of micrometers in diameter. Among early colonizing populations, we observed dominance of five genotypes that differed from each other by the extent of recombination with the late colonizing type. Our analyses suggest that the specific recombinant variant within the early colonizing group is selected for by environmental parameters such as temperature, consistent with recombination as a mechanism for ecological fine tuning. Evolutionary signatures, and strain-resolved expression patterns measured via mass spectrometry–based proteomics, indicate increased cobalamin biosynthesis, (de)methylation, and glycine cleavage in the late colonizer. This may suggest environmental changes within the biofilm during development, accompanied by redirection of compatible solutes from osmoprotectants toward metabolism. Across 27 communities, comparative proteogenomic analyses show that differential regulation of shared genes and expression of a small subset of the ∼15% of genes unique to each genotype are involved in niche partitioning. In summary, the results show how subtle genetic variations can lead to distinct ecological strategies.

Climate change impacts on soil microbial communities could alter the structure of terrestrial ecosystems and biogeochemical cycles of the Earth. We used 16S rRNA gene microarrays to evaluate changes in the composition of grassland soil microbial communities under rainfall amendments simulating alternative climate change scenarios, and to compare these to responses of overlying plants and invertebrates. Following 5 years of rainfall manipulation, soil bacteria and archaea in plots where natural rain was supplemented differed little from ambient controls, despite profound treatment-related changes in the overlying grassland. During the sixth and seventh year, seasonal differences in bacterial and archaeal assemblages emerged among treatments, but only when watering exacerbated or alleviated periods of particularly aberrant conditions in the ambient climate. In contrast to effects on plants and invertebrates, effects on bacteria and archaea did not compound across seasons or years, indicating that soil microbial communities may be more robust than associated aboveground macroorganisms to certain alterations in climate.

Fully understanding the biology of acid mine drainage (AMD) is central to our ability to control and manipulate its environmental impact. Although genomics and biogeochemical methods are relatively well established in the field, their combination with high-resolution imaging of intact members of microbial biofilm communities has not yet reached its full potential. Here, we used three-dimensional (3D) cryogenic electron tomography to determine the size and ultrastructure of intact ARMAN cells, a novel ultra-small archaeon, and sought evidence for their interactions with other members of its community. Within acid mine drainage biofilms, apparently free-living ARMAN cells from a deeply branched archaeal lineage have volumes of 0.009–0.04 μm3 (mean ~0.03±0.01 μm3), only ~92 ribosomes, yet are frequent hosts for replicating viruses. Organization within the periplasm and partitioning of ribosomes to the inner surface of the cytoplasmic membrane may be factors in size minimization. Most cells contain enigmatic tubular structures of unknown function. The low ribosome copy number per unit volume, indicative of slow growth rates and targeting of cells by diverse viruses may account for the low abundance of ARMAN cells compared with other biofilm community members. Our results provide the first 3D analysis of structural features of these novel and enigmatic cells and their interactions with at least two types of viruses. Our findings also emphasize that new biological phenomena remain to be discovered among lower abundance organisms from novel uncultivated lineages.

Enhanced biological phosphorus removal (EBPR) selects for polyphosphate accumulating microorganisms to achieve phosphate removal from wastewater. We used high-resolution community proteomics to identify key metabolic pathways in ‘Candidatus Accumulibacter phosphatis’ (A. phosphatis)-mediated EBPR and to evaluate the contributions of co-existing strains within the dominant population. Overall, 702 proteins from the A. phosphatis population were identified. Results highlight the importance of denitrification, fatty acid cycling and the glyoxylate bypass in EBPR. Strong similarity in protein profiles under anaerobic and aerobic conditions was uncovered (only 3% of A. phosphatis-associated proteins exhibited statistically significant abundance differences). By comprehensive genome-wide alignment of 13 930 orthologous proteins, we uncovered substantial differences in protein abundance for enzyme variants involved in both core-metabolism and EBPR-specific pathways among the A. phosphatis population. These findings suggest an essential role for genetic diversity in maintaining the stable performance of EBPR systems and, hence, demonstrate the power of integrated cultivation-independent genomics and proteomics for the analysis of complex biotechnological systems.

Microbes comprise the majority of extant organisms, yet much remains to be learned about the nature and driving forces of microbial diversification. Our understanding of how microorganisms adapt and evolve can be advanced by genome-wide documentation of the patterns of genetic exchange, particularly if analyses target coexisting members of natural communities. Here we use community genomic data sets to identify, with strain specificity, expressed proteins from the dominant member of a genomically uncharacterized, natural, acidophilic biofilm. Proteomics results reveal a genome shaped by recombination involving chromosomal regions of tens to hundreds of kilobases long that are derived from two closely related bacterial populations. Inter-population genetic exchange was confirmed by multilocus sequence typing of isolates and of uncultivated natural consortia. The findings suggest that exchange of large blocks of gene variants is crucial for the adaptation to specific ecological niches within the very acidic, metal-rich environment. Mass-spectrometry-based discrimination of expressed protein products that differ by as little as a single amino acid enables us to distinguish the behaviour of closely related coexisting organisms. This is important, given that microorganisms grouped together as a single species may have quite distinct roles in natural systems1, 2, 3 and their interactions might be key to ecosystem optimization. Because proteomic data simultaneously convey information about genome type and activity, strain-resolved community proteomics is an important complement to cultivation-independent genomic (metagenomic) analysis4, 5, 6 of microorganisms in the natural environment.