Explosive extrusion of cold material from the interior of icy bodies, or cryovolcanism, has been observed on Enceladus and, perhaps, Europa, Triton, and Ceres. It may explain the observed evidence for a young surface on Charon (Pluto’s surface is masked by frosts). Here, we evaluate prerequisites for cryovolcanism on dwarf planet-class Kuiper belt objects (KBOs). We first review the likely spatial and temporal extent of subsurface liquid, proposed mechanisms to overcome the negative buoyancy of liquid water in ice, and the volatile inventory of KBOs. We then present a new geochemical equilibrium model for volatile exsolution and its ability to drive upward crack propagation. This novel approach bridges geophysics and geochemistry, and extends geochemical modeling to the seldom-explored realm of liquid water at subzero temperatures. We show that carbon monoxide (CO) is a key volatile for gas-driven fluid ascent; whereas CO₂ and sulfur gases only play a minor role. N₂, CH₄, and H₂ exsolution may also drive explosive cryovolcanism if hydrothermal activity produces these species in large amounts (a few percent with respect to water). Another important control on crack propagation is the internal structure: a hydrated core makes explosive cryovolcanism easier, but an undifferentiated crust does not. We briefly discuss other controls on ascent such as fluid freezing on crack walls, and outline theoretical advances necessary to better understand cryovolcanic processes. Finally, we make testable predictions for the 2015 New Horizons flyby of the Pluto-Charon system.

Oligonucleotide signatures, especially tetranucleotide signatures, have been used as method for homology binning by exploiting an organism’s inherent biases towards the use of specific oligonucleotide words. Tetranucleotide signatures have been especially useful in environmental metagenomics samples as many of these samples contain organisms from poorly classified phyla which cannot be easily identified using traditional homology methods, including NCBI BLAST. This study examines oligonucleotide signatures across 1,424 completed genomes from across the tree of life, substantially expanding upon previous work. A comprehensive analysis of mononucleotide through nonanucleotide word lengths suggests that longer word lengths substantially improve the classification of DNA fragments across a range of sizes of relevance to high throughput sequencing. We find that, at present, heptanucleotide signatures represent an optimal balance between prediction accuracy and computational time for resolving taxonomy using both genomic and metagenomic fragments. We directly compare the ability of tetranucleotide and heptanucleotide world lengths (tetranucleotide signatures are the current standard for oligonucleotide word usage analyses) for taxonomic binning of metagenome reads. We present evidence that heptanucleotide word lengths consistently provide more taxonomic resolving power, particularly in distinguishing between closely related organisms that are often present in metagenomic samples. This implies that longer oligonucleotide word lengths should replace tetranucleotide signatures for most analyses. Finally, we show that the application of longer word lengths to metagenomic datasets leads to more accurate taxonomic binning of DNA scaffolds and have the potential to substantially improve taxonomic assignment and assembly of metagenomic data.

Uncovering the chemical and physical links between natural environments and microbial communities is becoming increasingly amenable owing to geochemical observations and metagenomic sequencing. At the hot spring known as Bison Pool in Yellowstone National Park, the cooling of the water in the outflow channel is associated with an increase in oxidation potential estimated from multiple field-based measurements. Representative groups of proteins whose sequences were derived from metagenomic data also exhibit an increase in average oxidation state of carbon in the protein molecules with distance from the hot-spring source. The energetic requirements of reactions to form selected proteins used in the model were computed using amino-acid group additivity for the standard molal thermodynamic properties of the proteins, and the relative chemical stabilities of the proteins were investigated by varying temperature, pH and oxidation state, expressed as activity of dissolved hydrogen. The relative stabilities of the proteins were found to track the locations of the sampling sites when the calculations included a function for hydrogen activity that increases with temperature and is higher, or more reducing, than values consistent with measurements of dissolved oxygen, sulfide and oxidation-reduction potential in the field. These findings imply that spatial patterns in the amino acid compositions of proteins can be linked, through energetics of overall chemical reactions representing the formation of the proteins, to the environmental conditions at this hot spring, even if microbial cells maintain considerably different internal conditions. Further applications of the thermodynamic calculations are possible for other natural microbial ecosystems.

We have constructed a conceptual model of biogeochemical cycles and metabolic and microbial community shifts within a hot spring ecosystem via coordinated analysis of the “Bison Pool” (BP) Environmental Genome and a complementary contextual geochemical dataset of ∼75 geochemical parameters. 2,321 16S rRNA clones and 470 megabases of environmental sequence data were produced from biofilms at five sites along the outflow of BP, an alkaline hot spring in Sentinel Meadow (Lower Geyser Basin) of Yellowstone National Park. This channel acts as a >22 m gradient of decreasing temperature, increasing dissolved oxygen, and changing availability of biologically important chemical species, such as those containing nitrogen and sulfur. Microbial life at BP transitions from a 92°C chemotrophic streamer biofilm community in the BP source pool to a 56°C phototrophic mat community. We improved automated annotation of the BP environmental genomes using BLAST-based Markov clustering. We have also assigned environmental genome sequences to individual microbial community members by complementing traditional homology-based assignment with nucleotide word-usage algorithms, allowing more than 70% of all reads to be assigned to source organisms. This assignment yields high genome coverage in dominant community members, facilitating reconstruction of nearly complete metabolic profiles and in-depth analysis of the relation between geochemical and metabolic changes along the outflow. We show that changes in environmental conditions and energy availability are associated with dramatic shifts in microbial communities and metabolic function. We have also identified an organism constituting a novel phylum in a metabolic “transition” community, located physically between the chemotroph- and phototroph-dominated sites. The complementary analysis of biogeochemical and environmental genomic data from BP has allowed us to build ecosystem-based conceptual models for this hot spring, reconstructing whole metabolic networks in order to illuminate community roles in shaping and responding to geochemical variability.

Over 100 hot spring sediment samples were collected from 28 sites in 12 areas/regions, while recording as many coincident geochemical properties as feasible (>60 analytes). PCR was used to screen samples for Korarchaeota 16S rRNA genes. Over 500 Korarchaeota 16S rRNA genes were screened by RFLP analysis and 90 were sequenced, resulting in identification of novel Korarchaeota phylotypes and exclusive geographical variants. Korarchaeota diversity was low, as in other terrestrial geothermal systems, suggesting a marine origin for Korarchaeota with subsequent niche-invasion into terrestrial systems. Korarchaeota endemism is consistent with endemism of other terrestrial thermophiles and supports the existence of dispersal barriers. Korarchaeota were found predominantly in >55°C springs at pH 4.7–8.5 at concentrations up to 6.6×10⁶ 16S rRNA gene copies g^-1 wet sediment. In Yellowstone National Park (YNP), Korarchaeota were most abundant in springs with a pH range of 5.7 to 7.0. High sulfate concentrations suggest these fluids are influenced by contributions from hydrothermal vapors that may be neutralized to some extent by mixing with water from deep geothermal sources or meteoric water. In the Great Basin (GB), Korarchaeota were most abundant at spring sources of pH<7.2 with high particulate C content and high alkalinity, which are likely to be buffered by the carbonic acid system. It is therefore likely that at least two different geological mechanisms in YNP and GB springs create the neutral to mildly acidic pH that is optimal for Korarchaeota. A classification support vector machine (C-SVM) trained on single analytes, two analyte combinations, or vectors from non-metric multidimensional scaling models was able to predict springs as Korarchaeota-optimal or sub-optimal habitats with accuracies up to 95%. To our knowledge, this is the most extensive analysis of the geochemical habitat of any high-level microbial taxon and the first application of a C-SVM to microbial ecology.

Streamer biofilm communities (SBC) are often observed within chemosynthetic zones of Yellowstone hot spring outflow channels, where temperatures exceed those conducive to photosynthesis. Nearest the hydrothermal source (75–88°C) SBC comprise thermophilic Archaea and Bacteria, often mixed communities including Desulfurococcales and uncultured Crenarchaeota, as well as Aquificae and Thermus, each carrying diagnostic membrane lipid biomarkers. We tested the hypothesis that SBC can alternate their metabolism between autotrophy and heterotrophy depending on substrate availability. Feeding experiments were performed at two alkaline hot springs in Yellowstone National Park: Octopus Spring and “Bison Pool,” using various ¹³C-labeled substrates (bicarbonate, formate, acetate, and glucose) to determine the relative uptake of these different carbon sources. Highest ¹³C uptake, at both sites, was from acetate into almost all bacterial fatty acids, particularly into methyl-branched C₁₅, C₁₇ and C₁₉ fatty acids that are diagnostic for Thermus/Meiothermus, and some Firmicutes as well as into universally common C_16:0 and C_18:0 fatty acids. ¹³C-glucose showed a similar, but a 10–30 times lower uptake across most fatty acids. ¹³C-bicarbonate uptake, signifying the presence of autotrophic communities was only significant at “Bison Pool” and was observed predominantly in non-specific saturated C₁₆, C₁₈, C₂₀, and C₂₂ fatty acids. Incorporation of ¹³C-formate occurred only at very low rates at “Bison Pool” and was almost undetectable at Octopus Spring, suggesting that formate is not an important carbon source for SBC. ¹³C-uptake into archaeal lipids occurred predominantly with ¹³C-acetate, suggesting also that archaeal communities at both springs have primarily heterotrophic carbon assimilation pathways. We hypothesize that these communities are energy-limited and predominantly nurtured by input of exogenous organic material, with only a small fraction being sustained by autotrophic growth.

Background: The metabolic strategies employed by microbes inhabiting natural systems are, in large part, dictated by the physical and geochemical properties of the environment. This study sheds light onto the complex relationship between biology and environmental geochemistry using forty-three metagenomes collected from geochemically diverse and globally distributed natural systems. It is widely hypothesized that many uncommonly measured geochemical parameters affect community dynamics and this study leverages the development and application of multidimensional biogeochemical metrics to study correlations between geochemistry and microbial ecology. Analysis techniques such as a Markov cluster-based measure of the evolutionary distance between whole communities and a principal component analysis (PCA) of the geochemical gradients between environments allows for the determination of correlations between microbial community dynamics and environmental geochemistry and provides insight into which geochemical parameters most strongly influence microbial biodiversity.

Results: By progressively building from samples taken along well defined geochemical gradients to samples widely dispersed in geochemical space this study reveals strong links between the extent of taxonomic and functional diversification of resident communities and environmental geochemistry and reveals temperature and pH as the primary factors that have shaped the evolution of these communities. Moreover, the inclusion of extensive geochemical data into analyses reveals new links between geochemical parameters (e.g. oxygen and trace element availability) and the distribution and taxonomic diversification of communities at the functional level. Further, an overall geochemical gradient (from multivariate analyses) between natural systems provides one of the most complete predictions of microbial taxonomic and functional composition.

Conclusions: Clustering based on the frequency in which orthologous proteins occur among metagenomes facilitated accurate prediction of the ordering of community functional composition along geochemical gradients, despite a lack of geochemical input. The consistency in the results obtained from the application of Markov clustering and multivariate methods to distinct natural systems underscore their utility in predicting the functional potential of microbial communities within a natural system based on system geochemistry alone, allowing geochemical measurements to be used to predict purely biological metrics such as microbial community composition and metabolism.

Carbon can be a major constituent of crustal and mantle fluids, occurring both as dissolved ionic species (e.g., carbonate ions or organic acids) and molecular species (e.g., CO₂, CO, CH₄, and more complex organic compounds). The chemistry of dissolved carbon changes dramatically with pressure (P) and temperature (T). In aqueous fluids at low P and T, molecular carbon gas species such as CO₂ and CH₄ saturate at low concentration to form a separate phase. With modest increases in P and T, these molecular species become fully miscible with H₂O, enabling deep crustal and mantle fluids to become highly concentrated in carbon. At such high concentrations, carbon species play an integral role as solvent components and, with H₂O, control the mobility of rock-forming elements in a wide range of geologic settings. The migration of carbon-bearing crustal and mantle fluids contributes to Earth’s carbon cycle; however, the mechanisms, magnitudes, and time variations of carbon transfer from depth to the surface remain least understood parts of the global carbon budget (Berner 1991, 1994; Berner and Kothavala 2001).

Microvirga sp. BSC39 was isolated from a biological soil crust near Moab, Utah. The strain appears to be capable of chemotaxis and exopolysaccharide synthesis for biofilm adhesion. The BSC39 genome contains iron siderophore uptake and hydrolysis enzymes; however, it lacks siderophore synthesis pathways, suggesting the uptake of siderophores produced by neighboring microbes.