Using thermodynamics to obtain geochemical information from genomes

Thermodynamic characterization of the relative stabilities of chemical compounds is a pillar of conceptual models in various fields of geosciences. Analogous models applied to genomes can yield new information about the relationship between genomes and their geochemical environments. In this perspective article, we present a chemical and thermodynamic analysis of prokaryotic lineages that have been the target of previous phylogenomic studies of evolutionary adaptation to varying redox conditions. The thermodynamic model development begins by quantifying the effects of hydrogen activity (aH₂) and temperature on the relative stabilities of organic compounds with different carbon oxidation state. When applied to proteins instead of metabolites, the same techniques can be used to identify combinations of aH₂ and temperature at which reference proteomes for Class I or Class II methanogens are relatively stable. The calculated aH₂ values are compatible with reported measurements for habitats of methanogens ranging from highly reducing submarine hydrothermal systems to less reducing environments including methanogenic sediments. In contrast to the transition between the two classes of methanogenic archaea, that between basal and terrestrial groups of Thaumarchaeota (denoting the origin of ammonia-oxidizing archaea) occurs at a less-reducing redox boundary. These examples reveal the consequences of energy minimization driving evolution and show how geochemical calculations involving biomolecules can be used to quantify and better understand the coevolution of the geosphere and biosphere.     

Quantifying the energetics of metabolic reactions in diverse biogeochemical systems: electron flow and ATP synthesis

Classical solution chemistry and thermodynamics have been used to quantify the energetics of two fundamental bioenergetic processes over a range of environmental conditions in which life is known to thrive: (1) the flow of electrons originating from an electron donor through the nicotinamide adenine dinucleotide (NAD) redox couple to an electron acceptor, and (2) the synthesis of ATP from ADP and aqueous phosphate. The approach taken explicitly accounts for the chemical formulas, charge states, complexation of ADP and ATP with magnesium, and the thermodynamic properties of individual biochemical species in stoichiometric and charge‐balanced reactions among biomolecules and other compounds. Because these species are represented in the reactions by explicit formula units, the chemical and thermodynamic consequences of the reactions can be evaluated as a function of temperature, pH, and bulk composition. To illustrate the utility of this approach, the energetics of the oxidation of glucose and hydrogen and the reduction of oxygen and sulfate coupled to NAD were characterized as functions of pH and temperature. The thermodynamic drive (chemical affinity, A ) for glucose to reduce NAD increases as the temperature increases, whereas the opposite is true for hydrogen. Also, at lower pHs, the chemical affinity is lower for these two electron donors to reduce NAD than at higher pHs. Similarly, the chemical affinity for oxidation of NAD by oxygen decreases by more than 2 kcal mol⁻¹ as temperature increases from 0 to 125 °C but the chemical affinity for oxidation of NAD by sulfate decreases by less than 1 kcal mol⁻¹. Calculations were also carried out to quantify the energetics of the synthesis of ATP for different bulk compositions, pHs, and temperatures. The chemical affinity for the synthesis of MgATP from MgADP and aqueous monophosphate at pH 5 minimizes with increasing temperature from 0 to 125 °C, which is not the case at pH 7 or 9. The procedures employed in these various calculations can be used to better understand how different environmental variables influence biogeochemical interactions. In addition, they help constrain the minimum energy required to sustain a particular microbial population and provide the means to determine why certain types of metabolism occur in the environments in which they do.     

Nitrospira-dominated biofilm within a thermal artesian spring: a case for nitrification-driven primary production in a geothermal setting

Water chemistry, energetic modeling, and molecular analyses were combined to investigate the microbial ecology of a biofilm growing in a thermal artesian spring within Hot Springs National Park, AR. This unique fresh water spring has a low dissolved chemical load and is isolated from both light and direct terrestrial carbon input - resulting in an oligotrophic ecosystem limited for fixed carbon and electron donors. Evaluation of energy yields of lithotrophic reactions putatively linked to autotrophy identified the aerobic oxidation of methane, hydrogen, sulfide, ammonia, and nitrite as the most exergonic. Small subunit (SSU) rRNA gene libraries from biofilm revealed a low-diversity microbial assemblage populated by bacteria and archaea at a gene copy ratio of 45:1. Members of the bacterial family 'Nitrospiraceae', known for their autotrophic nitrite oxidation, dominated the bacterial SSU rRNA gene library (approximately 45%). Members of the Thaumarchaeota ThAOA/HWCGIII (>96%) and Thaumarchaeota Group I.1b (2.5%), which both contain confirmed autotrophic ammonia oxidizers, dominated the archaeal SSU rRNA library. Archaea appear to dominate among the ammonia oxidizers, as only ammonia monooxygenase subunit A (amoA) genes belonging to members of the Thaumarchaeota were detected. The geochemical, phylogenetic, and genetic data support a model that describes a novel thermophilic biofilm built largely by an autotrophic nitrifying microbial assemblage. This is also the first observation of 'Nitrospiraceae' as the dominant organisms within a geothermal environment.     

Comparison of Metagenomic DNA Extraction Methods for Soil Sediments of High Elevation Puga Hot Spring in Ladakh,India to Explore Bacterial Diversity

Extraction of good-quality metagenomic DNA from extreme environments is quite challenging, particularly from high elevation hot spring sediments. Low microbial load, high humic acid content and other contaminants complicate the process of extraction of metagenomic DNA from hot spring sediments. In the present study, efficacy of five manual DNA extraction protocols with modifications has been evaluated for metagenomic DNA extraction from boron–sulfur rich high elevation Puga hot spring sediments. Best suited protocol was identified based on the cell lysis efficiency, DNA yield, humic acid content, PCR reproducibility and representation of bacterial diversity. Quantity as well as quality of crude metagenomic DNA differed remarkably between various protocols used and were not pure enough to give PCR amplification using 16S rRNA bacterial and archaeal primers. Crude metagenomic DNA extracted using five different DNA extraction protocols was purified using spin column based purification method. Even after purification, only three protocols C, D and E yielded metagenomic DNA that could be amplified using both archaeal and bacterial primers. To evaluate the degree of microbial diversity represented by protocols C, D and E, phylogenetic genes amplified were subjected to amplified ribosomal DNA restriction analysis (ARDRA) and denaturing gradient gel electrophoresis analysis (DGGE) analysis. ARDRA banding pattern of amplicons generated for all the three extraction protocols, i.e., C, D and E were found to be similar. DGGE of protocol E derived amplicons resulted in the similar number of dominant bands but a greater number of non-dominant bands, i.e., the highest microbial diversity in comparison to protocols C and D, respectively. In the present study, protocol E developed from Yeates et al. protocol has been found to be best in terms of DNA yield, DNA purity and bacterial diversity depiction associated with boron–sulfur rich sediment of high elevation hot springs.     

Archaeal and bacterial communities in geochemically diverse hot springs of Yellowstone National Park, USA

Microbiological and geochemical surveys were conducted at three hot springs (Obsidian Pool, Sylvan Spring, and ‘Bison Pool’) in Yellowstone National Park (Wyoming, USA). Microbial community structure was investigated by polymerase chain reaction (PCR) amplification of 16S rRNA gene sequences from DNA extracted from sediments of each hot spring, followed by molecular cloning. Both bacterial and archaeal DNA was retrieved from all samples. No Euryarchaea were found, but diverse Crenarchaea exist in all three pools, particularly affiliating with deep‐branching, but uncultivated organisms. In addition, cloned DNA affiliating with the Desulphurococcales and Thermoproteales was identified, but the distribution of taxa differs in each hot spring. The bacterial community at all three locations is dominated by members of the Aquificales and Thermodesulfobacteriales, indicating that the ‘knallgas’ reaction (aerobic hydrogen oxidation) may be a central metabolism in these ecosystems. To provide geochemical context for the microbial community structures, energy‐yields for a number of chemolithoautotrophic reactions are provided for >80 sampling sites in Yellowstone, including Obsidian Pool, Sylvan Spring, and ‘Bison Pool’. This energy profile shows that the knallgas reaction is just one of many exergonic reactions in the Yellowstone hot springs, that energy‐yields for certain reactions can vary substantially from one site to the next, and that few of the demonstrated exergonic reactions are known to support microbial metabolism.     

On the energetics of chemolithotrophy in nonequilibrium systems: case studies of geothermal springs in Yellowstone National Park

Chemolithotrophic micro‐organisms are important primary producers in high‐temperature geothermal environments and may catalyse a number of different energetically favourable redox reactions as a primary energy source. Analysis of geochemical constituents followed by chemical speciation and subsequent calculation of reaction free energies (ΔG_rxn) is a useful tool for evaluating the thermodynamic favourability and potential energy available for microbial metabolism. The primary goal of this study was to examine relationships among geochemical gradients and microbial population distribution, and to evaluate the utility of energetic approaches for predicting microbial metabolism from free‐energy calculations, utilizing as examples, several geothermal habitats in Yellowstone National Park where thorough geochemical and phylogenetic analyses have been performed. Acidic (pH ～ 3) and near‐neutral (pH ～ 6–7) geothermal springs were chosen for their range in geochemical properties. Aqueous and solid phase samples obtained from the source pools and the outflow channels of each spring were characterized for all major chemical constituents using laboratory and field methods to accurately measure the concentrations of predominant oxidized and reduced species. Reaction free energies (ΔG_rxn) for 33 oxidation–reduction reactions potentially important to chemolithotrophic micro‐organisms were calculated at relevant spring temperatures after calculating ion activities using an aqueous equilibrium model. Free‐energy values exhibit significant variation among sites for reactions with pH dependence. For example, free‐energy values for reactions involving Fe³⁺ are especially variable across sites due in large part to the pH dependence of Fe³⁺ activity, and exhibit changes of up to 40 kJ mol^?1 electron from acidic to near neutral geothermal springs. Many of the detected 16S rRNA gene sequences represent organisms whose metabolisms are consistent with exergonic processes. However, sensitivity analyses demonstrated that reaction free energies do not generally represent the steep gradients in local geochemical conditions resulting from air–water gas exchange and solid phase deposition that are important in defining microbial habitats and 16S rRNA gene sequence distribution within geothermal outflow channels.     

Calculation of total free energy yield as an alternative approach for predicting the importance of potential chemolithotrophic reactions in geothermal springs

To inform hypotheses regarding the relative importance of chemolithotrophic metabolisms in geothermal environments, we calculated free energy yields of 26 chemical reactions potentially supporting chemolithotrophy in two US Great Basin hot springs, taking into account the effects of changing reactant and product activities on the Gibbs free energy as each reaction progressed. Results ranged from 1.2 × 10(-5) to 3.6 J kg(-1) spring water, or 3.7 × 10(-5) to 11.5 J s(-1) based on measured flow rates, with aerobic oxidation of CH(4) or NH4 + giving the highest average yields. Energy yields calculated without constraining pH were similar to those at constant pH except for reactions where H(+) was consumed, which often had significantly lower yields when pH was unconstrained. In contrast to the commonly used normalization of reaction chemical affinities per mole of electrons transferred, reaction energy yields for a given oxidant varied by several orders of magnitude and were more sensitive to differences in the activities of products and reactants. The high energy yield of aerobic ammonia oxidation is consistent with previous observations of significant ammonia oxidation rates and abundant ammonia-oxidizing archaea in sediments of these springs. This approach offers an additional lens through which to view the thermodynamic landscape of geothermal springs.     

Protein folding stabilities are a major determinant of oxidation rates for buried methionine residues

The oxidation of protein-bound methionines to form methionine sulfoxides has a broad range of biological ramifications, making it important to delineate factors that influence methionine oxidation rates within a given protein. This is especially important for biopharmaceuticals, where oxidation can lead to deactivation and degradation. Previously, neighboring residue effects and solvent accessibility have been shown to impact the susceptibility of methionine residues to oxidation. In this study, we provide proteome-wide evidence that oxidation rates of buried methionine residues are also strongly influenced by the thermodynamic folding stability of proteins. We surveyed the Escherichia coli proteome using several proteomic methodologies and globally measured oxidation rates of methionine residues in the presence and absence of tertiary structure, as well as the folding stabilities of methionine-containing domains. These data indicated that buried methionines have a wide range of protection factors against oxidation that correlate strongly with folding stabilities. Consistent with this, we show that in comparison to E. coli, the proteome of the thermophile Thermus thermophilus is significantly more stable and thus more resistant to methionine oxidation. To demonstrate the utility of this correlation, we used native methionine oxidation rates to survey the folding stabilities of E. coli and T. thermophilus proteomes at various temperatures and propose a model that relates the temperature dependence of the folding stabilities of these two species to their optimal growth temperatures. Overall, these results indicate that oxidation rates of buried methionines from the native state of proteins can be used as a metric of folding stability.     

Occurrence and Carbon Metabolism of Green Nonsulfur-like Bacteria in Californian and Nevada Hot Spring Microbial mats as Revealed by Wax Ester Lipid Analysis

Green nonsulfur-like bacteria (GNSLB) in Yellowstone hot spring microbial mats have been extensively studied and are thought to operate both as photoheterotrophs and photoautotrophs. Here we studied the occurrence and carbon metabolisms of GNSLB by analyzing the distribution and isotopic composition of their characteristic wax ester lipids in four Californian and Nevada hot spring microbial mats at a range of temperatures (37–96°C). The distribution of wax esters varied strongly with temperature. At temperatures between 50–60°C the wax ester composition in each of the four hot spring microbial mats was dominated by C₃₀ to C₃₆ wax esters, consisting of mixtures of C₁₅-C₁₈ n-alkyl and branched fatty acids and alcohols, typical for GNSLB. Stable carbon isotopic analysis showed that these wax esters were only depleted by 5 to 10‰ compared to dissolved inorganic carbon in the overlying water, suggesting that these GNSLB were mainly autotrophic. However, analysis of different depth layers of one microbial mat showed that these GNSLB wax esters were increasingly depleted in ¹³C with depth, suggesting that photoautotrophy mainly occurred in the top layer of the mat. ¹³C-depleted C₃₆-C₄₄ wax esters were found in one hot spring at high temperatures (77–96°C) and are likely derived from allochtonous plant waxes. At several lower temperature sites (35–40°C) the wax esters were predominantly composed of C₂₈, C₃₀ and C₃₂ wax esters consisting of mixtures of C₁₄-C₁₆ fatty acids and n-alkanols and were depleted in ¹³C by 15–20‰ relative to dissolved inorganic carbon, suggesting they may be derived from heterotrophic organisms. Our results indicate that autotrophic GNSLB occur widely in hot springs and that diverse groups of organisms contribute to the pool of wax ester lipids in hot spring environments.     

Energy landscapes shape microbial communities in hydrothermal systems on the Arctic Mid-Ocean Ridge

Methods developed in geochemical modelling combined with recent advances in molecular microbial ecology provide new opportunities to explore how microbial communities are shaped by their chemical surroundings. Here, we present a framework for analyses of how chemical energy availability shape chemotrophic microbial communities in hydrothermal systems through an investigation of two geochemically different basalt-hosted hydrothermal systems on the Arctic Mid-Ocean Ridge: the Soria Moria Vent field (SMVF) and the Loki''s Castle Vent Field (LCVF). Chemical energy landscapes were evaluated through modelling of the Gibbs energy from selected redox reactions under different mixing ratios between seawater and hydrothermal fluids. Our models indicate that the sediment-influenced LCVF has a much higher potential for both anaerobic and aerobic methane oxidation, as well as aerobic ammonium and hydrogen oxidation, than the SMVF. The modelled energy landscapes were used to develop microbial community composition models, which were compared with community compositions in environmental samples inside or on the exterior of hydrothermal chimneys, as assessed by pyrosequencing of partial 16S rRNA genes. We show that modelled microbial communities based solely on thermodynamic considerations can have a high predictive power and provide a framework for analyses of the link between energy availability and microbial community composition.     

Metagenomic resolution of microbial functions in deep-sea hydrothermal plumes across the Eastern Lau Spreading Center

Microbial processes within deep-sea hydrothermal plumes affect ocean biogeochemistry on global scales. In rising hydrothermal plumes, a combination of microbial metabolism and particle formation processes initiate the transformation of reduced chemicals like hydrogen sulfide, hydrogen, methane, iron, manganese and ammonia that are abundant in hydrothermal vent fluids. Despite the biogeochemical importance of this rising portion of plumes, it is understudied in comparison to neutrally buoyant plumes. Here we use metagenomics and bioenergetic modeling to describe the abundance and genetic potential of microorganisms in relation to available electron donors in five different hydrothermal plumes and three associated background deep-sea waters from the Eastern Lau Spreading Center located in the Western Pacific Ocean. Three hundred and thirty one distinct genomic ‘bins'' were identified, comprising an estimated 951 genomes of archaea, bacteria, eukarya and viruses. A significant proportion of these genomes is from novel microorganisms and thus reveals insights into the energy metabolism of heretofore unknown microbial groups. Community-wide analyses of genes encoding enzymes that oxidize inorganic energy sources showed that sulfur oxidation was the most abundant and diverse chemolithotrophic microbial metabolism in the community. Genes for sulfur oxidation were commonly present in genomic bins that also contained genes for oxidation of hydrogen and methane, suggesting metabolic versatility in these microbial groups. The relative diversity and abundance of genes encoding hydrogen oxidation was moderate, whereas that of genes for methane and ammonia oxidation was low in comparison to sulfur oxidation. Bioenergetic-thermodynamic modeling supports the metagenomic analyses, showing that oxidation of elemental sulfur with oxygen is the most dominant catabolic reaction in the hydrothermal plumes. We conclude that the energy metabolism of microbial communities inhabiting rising hydrothermal plumes is dictated by the underlying plume chemistry, with a dominant role for sulfur-based chemolithoautotrophy.     

NMR investigations on residue level unfolding thermodynamics in DLC8 dimer by temperature dependent native state hydrogen exchange

Understanding protein stability at residue level detail in the native state ensemble of a protein is crucial to understanding its biological function. At the same time, deriving thermodynamic parameters using conventional spectroscopic and calorimetric techniques remains a major challenge for some proteins due to protein aggregation and irreversibility of denaturation at higher temperature values. In this regard, we describe here the NMR investigations on the conformational stabilities and related thermodynamic parameters such as local unfolding enthalpies, heat capacities and transition midpoints in DLC8 dimer, by using temperature dependent native state hydrogen exchange; this protein aggregates at high (>65°C) temperatures. The stability (free energy) of the native state was found to vary substantially with temperature at every residue. Significant differences were found in the thermodynamic parameters at individual residue sites indicating that the local environments in the protein structure would respond differently to external perturbations; this reflects on plasticity differences in different regions of the protein. Further, comparison of this data with similar data obtained from GdnHCl dependent native state hydrogen exchange indicated many similarities at residue level, suggesting that local unfolding transitions may be similar in both the cases. This has implications for the folding/unfolding mechanisms of the protein. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Phase transition in superdense hydrogen and deuterium

Results are presented from numerical simulations of the thermodynamic properties of superdense hydrogen and deuterium plasmas by the Monte Carlo method and from calculations by a multicomponent chemical model. The results obtained reveal the anomalous behavior of the thermodynamic functions and composition of molecular gas plasmas in the submegabar and megabar pressure ranges. Such behavior is interpreted as a dissociative phase transition. The results of calculations by the chemical model are compared with the experimental data on the equation of state and conductivity of hydrogen and deuterium plasmas.     

Coordinating environmental genomics and geochemistry reveals metabolic transitions in a hot spring ecosystem

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.     

Role of a reaction coordinate in free energy calculations

Abstract

The free energy calculation method emerges as a viable technique for ‘in-silico’ calorimetry. Efficient sampling techniques and the good choice of a reaction path connecting the reactant and the product state enable accurate computations of the free energy differences. We argue that in many cases the thermodynamic integration technique has the lowest variance when the transformation between the reactant and the product state proceeds along the natural path of the studied chemical reaction. We provide examples of free energy calculations for the fragmentation of the charged clusters and the swapping reaction of oligomer formation in proteins that follow a tentative reaction mechanism.     

Field- and Lab-Based Potentiometric Titrations of Microbial Mats from the Fairmont Hot Spring,Canada

Potentiometric titrations are an effective tool to constrain the protonation constants and site concentrations for microbial surface ligands. Protonation models developed from these experiments are often coupled with data from metal adsorption experiments to calculate microbial ligand-metal binding constants. Ultimately, the resulting surface complexation models can be used to predict metal immobilization behavior across diverse chemical conditions. However, most protonation and metal-ligand thermodynamic constants have been generated in laboratory experiments that use cultured microbes which may differ in their chemical reactivity from environmental samples. In this study, we investigate the use of in situ field potentiometric titrations of microbial mats at a carbonate hot spring located at Fairmont Hot Springs, British Columbia, with the aim to study microbial reactivities in a natural field system. We found that authigenic carbonate minerals complicated the potentiometric titration process due to a “carbonate spike” introduced by the contribution of inorganic carbonate mineral dissolution and subsequent carbonate speciation changes during the transition from low to high pH. This inhibits the determination of microbial surface ligand variety and concentrations. Our preliminary study also highlights the need for developing novel probes to quantify in situ microbial mat reactivity in future field investigations.     

Comparative oxidation studies of methionine residues reflect a structural effect on chemical kinetics in rhG-CSF

The effect of protein conformation on the rate of chemical degradation is poorly understood. To address the role of structure on chemical degradation kinetics, comparative oxidation studies of methionine residues in recombinant human granulocyte colony-stimulating factor (rhG-CSF) were performed. The kinetics of oxidation of methionine residues by hydrogen peroxide (H2O2) in rhG-CSF and corresponding chemically synthesized peptides thereof was measured at different temperatures. To assess structural effects, equilibrium denaturation experiments also were conducted on rhG-CSF, yielding the free energy of unfolding as a function of temperature. A comparison of the relative rates of oxidation of methionine residues in short peptides with those of corresponding methionine residues in rhG-CSF yields an understanding of how protein tertiary structure affects oxidation reactions. For the temperature range that was studied, 4-45 degrees C, the oxidation rate constants followed an Arrhenius equation quite well, suggesting the lack of temperature-induced local structural perturbations that affect chemical degradation rates. One of the four methionine residues, Met 122, exhibited an activation energy significantly different from that of the corresponding peptide. Extrapolation of kinetic data predicts non-Arrhenius behavior around the melting temperature. Three phenomenological models based on different mechanisms are discussed, and an application to shelf life prediction of pharmaceuticals is presented.     

Microbial Community Structure and Arsenic Biogeochemistry in an Acid Vapor-Formed Spring in Tengchong Geothermal Area,China

Arsenic biogeochemistry has been studied extensively in acid sulfate-chloride hot springs, but not in acid sulfate hot springs with low chloride. In this study, Zhenzhuquan in Tengchong geothermal area, a representative acid sulfate hot spring with low chloride, was chosen to study arsenic geochemistry and microbial community structure using Illumina MiSeq sequencing. Over 0.3 million 16S rRNA sequence reads were obtained from 6-paired parallel water and sediment samples along its outflow channel. Arsenic oxidation occurred in the Zhenxhuquan pool, with distinctly high ratios of arsenate to total dissolved arsenic (0.73–0.86). Coupled with iron and sulfur oxidation along the outflow channel, arsenic accumulated in downstream sediments with concentrations up to 16.44 g/kg and appeared to significantly constrain their microbial community diversity. These oxidations might be correlated with the appearance of some putative functional microbial populations, such as Aquificae and Pseudomonas (arsenic oxidation), Sulfolobus (sulfur and iron oxidation), Metallosphaera and Acidicaldus (iron oxidation). Temperature, total organic carbon and dissolved oxygen significantly shaped the microbial community structure of upstream and downstream samples. In the upstream outflow channel region, most microbial populations were microaerophilic/anaerobic thermophiles and hyperthermophiles, such as Sulfolobus, Nocardia, Fervidicoccus, Delftia, and Ralstonia. In the downstream region, aerobic heterotrophic mesophiles and thermophiles were identified, including Ktedonobacteria, Acidicaldus, Chthonomonas and Sphingobacteria. A total of 72.41–95.91% unassigned-genus sequences were derived from the downstream high arsenic sediments 16S rRNA clone libraries. This study could enable us to achieve an integrated understanding on arsenic biogeochemistry in acid hot springs.     

On tautomerism of the cytosine molecule

Tautomerism of the cytosine molecule is discussed in connection with recent experimental matrix-isolation infrared spectroscopic measurements and recent ab initio calculations of relative stabilities of tautomers and of IR spectra for different tautomeric forms of the compound. Experimental IR spectra in the N-H and O-H stretching regions and in the C = O stretching region are presented for cytosine and for its several derivatives considered as model compounds. This experimental evidence, as well as the quantum-mechanical calculations (including both electron correlation and zero-point vibrational contributions), clearly indicate that two tautomers of cytosine, i.e. the amino-hydroxy and amino-oxo forms with the hydrogen atom at the N(1) position, exist in equilibrium when the cytosine molecule is isolated in an inert environment. The effect of the environment on the relative stabilities of several tautomers is also discussed briefly.     

Genomics, environmental genomics and the issue of microbial species

A microbial species concept is crucial for interpreting the variation detected by genomics and environmental genomics among cultivated microorganisms and within natural microbial populations. Comparative genomic analyses of prokaryotic species as they are presently described and named have led to the provocative idea that prokaryotes may not form species as we think about them for plants and animals. There are good reasons to doubt whether presently recognized prokaryotic species are truly species. To achieve a better understanding of microbial species, we believe it is necessary to (i) re-evaluate traditional approaches in light of evolutionary and ecological theory, (ii) consider that different microbial species may have evolved in different ways and (iii) integrate genomic, metagenomic and genome-wide expression approaches with ecological and evolutionary theory. Here, we outline how we are using genomic methods to (i) identify ecologically distinct populations (ecotypes) predicted by theory to be species-like fundamental units of microbial communities, and (ii) test their species-like character through in situ distribution and gene expression studies. By comparing metagenomic sequences obtained from well-studied hot spring cyanobacterial mats with genomic sequences of two cultivated cyanobacterial ecotypes, closely related to predominant native populations, we can conduct in situ population genetics studies that identify putative ecotypes and functional genes that determine the ecotypes' ecological distinctness. If individuals within microbial communities are found to be grouped into ecologically distinct, species-like populations, knowing about such populations should guide us to a better understanding of how genomic variation is linked to community function.