Microbial-mat communities in the effluent channels of Octopus and Mushroom Springs within the Lower Geyser Basin at Yellowstone National Park have been studied for nearly 50 years. The emphasis has mostly focused on the chlorophototrophic bacterial organisms of the phyla Cyanobacteria and Chloroflexi. In contrast, the diversity and metabolic functions of the heterotrophic community in the microoxic/anoxic region of the mat are not well understood. In this study we analyzed the orange-colored undermat of the microbial community of Mushroom Spring using metagenomic and rRNA-amplicon (iTag) analyses. Our analyses disclosed a highly diverse community exhibiting a high degree of unevenness, strongly dominated by a single taxon, the filamentous anoxygenic phototroph, Roseiflexus spp. The second most abundant organisms belonged to the Thermotogae, which have been hypothesized to be a major source of H 2 from fermentation that could enable photomixotrophic metabolism by Chloroflexus and Roseiflexus spp. Other abundant organisms include two members of the Armatimonadetes (OP10); Thermocrinis sp.; and phototrophic and heterotrophic members of the Chloroflexi. Further, an Atribacteria (OP9/JS1) member; a sulfate-reducing Thermodesulfovibrio sp.; a Planctomycetes member; a member of the EM3 group tentatively affiliated with the Thermotogae, as well as a putative member of the Arminicenantes (OP8) represented ≥1% of the reads. Archaea were not abundant in the iTag analysis, and no metagenomic bin representing an archaeon was identified. A high microdiversity of 16S rRNA gene sequences was identified for the dominant taxon, Roseiflexus spp. Previous studies demonstrated that highly similar Synechococcus variants in the upper layer of the mats represent ecological species populations with specific ecological adaptations. This study suggests that similar putative ecotypes specifically adapted to different niches occur within the undermat community, particularly for Roseiflexus spp.

Dynamic environmental factors such as light, nutrients, salt, and temperature continuously affect chlorophototrophic microbial mats, requiring adaptive and acclimative responses to stabilize composition and function. Quantitative metabolomics analysis can provide insights into metabolite dynamics for understanding community response to such changing environmental conditions. In this study, we quantified volatile organic acids, polar metabolites (amino acids, glycolytic and citric acid cycle intermediates, nucleobases, nucleosides, and sugars), wax esters, and polyhydroxyalkanoates, resulting in the identification of 104 metabolites and related molecules in thermal chlorophototrophic microbial mat cores collected over a diel cycle in Mushroom Spring, Yellowstone National Park. A limited number of predominant taxa inhabit this community and their functional potentials have been previously identified through metagenomic and metatranscriptomic analyses and in situ metabolisms, and metabolic interactions among these taxa have been hypothesized. Our metabolomics results confirmed the diel cycling of photorespiration (e.g., glycolate) and fermentation (e.g., acetate, propionate, and lactate) products, the carbon storage polymers polyhydroxyalkanoates, and dissolved gasses (e.g., H2 and CO2) in the waters overlying the mat, which were hypothesized to occur in major mat chlorophototrophic community members. In addition, we have formulated the following new hypotheses: (1) the morning hours are a time of biosynthesis of amino acids, DNA, and RNA; (2) photo-inhibited cells may also produce lactate via fermentation as an alternate metabolism; (3) glycolate and lactate are exchanged among Synechococcus and Roseiflexus spp.; and (4) fluctuations in many metabolite pools (e.g., wax esters) at different times of day result from species found at different depths within the mat responding to temporal differences in their niches.

Genomes were obtained for three closely related strains of Synechococcus that are representative of putative ecotypes (PEs) that predominate at different depths in the 1 mm-thick, upper-green layer in the 60°C mat of Mushroom Spring, Yellowstone National Park, and exhibit different light adaptation and acclimation responses. The genomes were compared to the published genome of a previously obtained, closely related strain from a neighboring spring, and differences in both gene content and orthologous gene alleles between high-light-adapted and low-light-adapted strains were identified. Evidence of genetic differences that relate to adaptation to light intensity and/or quality, CO2uptake, nitrogen metabolism, organic carbon metabolism, and uptake of other nutrients were found between strains of the different putative ecotypes. In situ diel transcription patterns of genes, including genes unique to either low-light-adapted or high-light-adapted strains and different alleles of an orthologous photosystem gene, revealed that expression is fine-tuned to the different light environments experienced by ecotypes prevalent at various depths in the mat. This study suggests that strains of closely related PEs have different genomic adaptations that enable them to inhabit distinct ecological niches while living in close proximity within a microbial community.

Closely related strains of thermophilic Synechococcus were cultivated from the microbial mats found in the effluent channels of Mushroom Spring, Yellowstone National Park (YNP). These strains have identical or nearly identical 16S rRNA sequences but are representative of separate, predicted putative ecotype (PE) populations, which were identified by using the more highly resolving psaA locus and which predominate at different vertical positions within the 1-mm-thick upper-green layer of the mat. Pyrosequencing confirmed that each strain contained a single, predominant psaA genotype. Strains differed in growth rate as a function of irradiance. A strain with a psaA genotype corresponding to a predicted PE that predominates near the mat surface grew fastest at high irradiances, whereas strains with psaA genotypes representative of predominant subsurface populations grew faster at low irradiance and exhibited greater sensitivity to abrupt shifts to high light. The high-light-adapted and low-light-adapted strains also exhibited differences in pigment content and the composition of the photosynthetic apparatus (photosystem ratio) when grown under different light intensities. Cells representative of the different strains had similar morphologies under low-light conditions, but under high-light conditions, cells of low-light-adapted strains became elongated and formed short chains of cells. Collectively, the results presented here are consistent with the hypothesis that closely related, but distinct, ecological species of Synechococcus occupy different light niches in the Mushroom Spring microbial mat and acclimate differently to changing light environments.

Recent studies of bacterial speciation have claimed to support the biological species concept—that reduced recombination is required for bacterial populations to diverge into species. This conclusion has been reached from the discovery that ecologically distinct clades show lower rates of recombination than that which occurs among closest relatives. However, these previous studies did not attempt to determine whether the more-rapidly recombining close relatives within the clades studied may also have diversified ecologically, without benefit of sexual isolation. Here we have measured the impact of recombination on ecological diversification within and between two ecologically distinct clades (A and B') of Synechococcus in a hot spring microbial mat in Yellowstone National Park, using a cultivation-free, multi-locus approach. Bacterial artificial chromosome (BAC) libraries were constructed from mat samples collected at 60°C and 65°C. Analysis of multiple linked loci near Synechococcus 16S rRNA genes showed little evidence of recombination between the A and B' lineages, but a record of recombination was apparent within each lineage. Recombination and mutation rates within each lineage were of similar magnitude, but recombination had a somewhat greater impact on sequence diversity than mutation, as also seen in many other bacteria and archaea. Despite recombination within the A and B' lineages, there was evidence of ecological diversification within each lineage. The algorithm Ecotype Simulation identified sequence clusters consistent with ecologically distinct populations (ecotypes), and several hypothesized ecotypes were distinct in their habitat associations and in their adaptations to different microenvironments. We conclude that sexual isolation is more likely to follow ecological divergence than to precede it. Thus, an ecology-based model of speciation appears more appropriate than the biological species concept for bacterial and archaeal diversification.