Sulfolobus islandicus meta-populations in Yellowstone National Park hot springs

A genus in the bacterial phylum Aquificota appears to be endemic to Aotearoa-New Zealand

Allopatric speciation has been difficult to examine among microorganisms, with prior reports of endemism restricted to sub-genus level taxa. Previous microbial community analysis via 16S rRNA gene sequencing of 925 geothermal springs from the Taupō Volcanic Zone (TVZ), Aotearoa-New Zealand, revealed widespread distribution and abundance of a single bacterial genus across 686 of these ecosystems (pH 1.2-9.6 and 17.4-99.8 °C). Here, we present evidence to suggest that this genus, Venenivibrio (phylum Aquificota), is endemic to Aotearoa-New Zealand. A specific environmental niche that increases habitat isolation was identified, with maximal read abundance of Venenivibrio occurring at pH 4-6, 50-70 °C, and low oxidation-reduction potentials. This was further highlighted by genomic and culture-based analyses of the only characterised species for the genus, Venenivibrio stagnispumantis CP.B2T, which confirmed a chemolithoautotrophic metabolism dependent on hydrogen oxidation. While similarity between Venenivibrio populations illustrated that dispersal is not limited across the TVZ, extensive amplicon, metagenomic, and phylogenomic analyses of global microbial communities from DNA sequence databases indicates Venenivibrio is geographically restricted to the Aotearoa-New Zealand archipelago. We conclude that geographic isolation, complemented by physicochemical constraints, has resulted in the establishment of an endemic bacterial genus.

Genomic insight and physiological characterization of thermoacidophilic Alicyclobacillus isolated from Yellowstone National Park

Introduction

Alicyclobacillus has been isolated from extreme environments such as hot springs, volcanoes, as well as pasteurized acidic beverages, because it can tolerate extreme temperatures and acidity. In our previous study, Alicyclobacillus was isolated during the enrichment of methane oxidizing bacteria from Yellowstone Hot Spring samples.

Methods

Physiological characterization and genomic exploration of two new Alicyclobacillus isolates, AL01A and AL05G, are the main focus of this study to identify their potential relationships with a thermoacidophilic methanotroph (Methylacidiphilum) isolated from the same hot spring sediments.

Results and discussion

In the present study, both Alicyclobacillus isolates showed optimal growth at pH 3.5 and 55°C, and contain ω-alicyclic fatty acids as a major lipid (ca. 60%) in the bacterial membrane. Genomic analysis of these strains revealed specific genes and pathways that the methanotroph genome does not have in the intermediary carbon metabolism pathway such as serC (phosphoserine aminotransferase), comA (phosphosulfolactate synthase), and DAK (glycerone kinase). Both Alicyclobacillus strains were also found to contain transporter systems for extracellular sulfate (ABC transporter), suggesting that they could play an important role in sulfur metabolism in this extreme environment. Genomic analysis of vitamin metabolism revealed Alicyclobacillus and Methylacidiphilum are able to complement each other's nutritional deficiencies, resulting in a mutually beneficial relationship, especially in vitamin B1(thiamin), B3 (niacin), and B7 (biotin) metabolism. These findings provide insights into the role of Alicyclobacillus isolates in geothermal environments and their unique metabolic adaptations to these environments.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Temporal dynamics of geothermal microbial communities in Aotearoa-New Zealand

Microbial biogeography studies, in particular for geothermal-associated habitats, have focused on spatial patterns and/or individual sites, which have limited ability to describe the dynamics of ecosystem behaviour. Here, we report the first comprehensive temporal study of bacterial and archaeal communities from an extensive range of geothermal features in Aotearoa-New Zealand. One hundred and fifteen water column samples from 31 geothermal ecosystems were taken over a 34-month period to ascertain microbial community stability (control sites), community response to both natural and anthropogenic disturbances in the local environment (disturbed sites) and temporal variation in spring diversity across different pH values (pH 3, 5, 7, 9) all at a similar temperature of 60–70°C (pH sites). Identical methodologies were employed to measure microbial diversity via 16S rRNA gene amplicon sequencing, along with 44 physicochemical parameters from each feature, to ensure confidence in comparing samples across timeframes. Our results indicated temperature and associated groundwater physicochemistry were the most likely parameters to vary stochastically in these geothermal features, with community abundances rather than composition more readily affected by a changing environment. However, variation in pH (pH ±1) had a more significant effect on community structure than temperature (±20°C), with alpha diversity failing to adequately measure temporal microbial disparity in geothermal features outside of circumneutral conditions. While a substantial physicochemical disturbance was required to shift community structures at the phylum level, geothermal ecosystems were resilient at this broad taxonomic rank and returned to a pre-disturbed state if environmental conditions re-established. These findings highlight the diverse controls between different microbial communities within the same habitat-type, expanding our understanding of temporal dynamics in extreme ecosystems.

Longitudinal analysis of the Five Sisters hot springs in Yellowstone National Park reveals a dynamic thermoalkaline environment

Research focused on microbial populations of thermoalkaline springs has been driven in a large part by the lure of discovering functional enzymes with industrial applications in high-pH and high temperature environments. While several studies have focused on understanding the fundamental ecology of these springs, the small molecule profiles of thermoalkaline springs have largely been overlooked. To better understand how geochemistry, small molecule composition, and microbial communities are connected, we conducted a three-year study of the Five Sisters (FS) springs that included high-resolution geochemical measurements, 16S rRNA sequencing of the bacterial and archaeal community, and mass spectrometry-based metabolite and extracellular small molecule characterization. Integration of the four datasets facilitated a comprehensive analysis of the interwoven thermoalkaline spring system. Over the course of the study, the microbial population responded to changing environmental conditions, with archaeal populations decreasing in both relative abundance and diversity compared to bacterial populations. Decreases in the relative abundance of Archaea were associated with environmental changes that included decreased availability of specific nitrogen- and sulfur-containing extracellular small molecules and fluctuations in metabolic pathways associated with nitrogen cycling. This multi-factorial analysis demonstrates that the microbial community composition is more closely correlated with pools of extracellular small molecules than with the geochemistry of the thermal springs. This is a novel finding and suggests that a previously overlooked component of thermal springs may have a significant impact on microbial community composition.

Viruses in Subsurface Environments

Over the past 20 years, our knowledge of virus diversity and abundance in subsurface environments has expanded dramatically through application of quantitative metagenomic approaches. In most subsurface environments, viral diversity and abundance rival viral diversity and abundance observed in surface environments. Most of these viruses are uncharacterized in terms of their hosts and replication cycles. Analysis of accessory metabolic genes encoded by subsurface viruses indicates that they evolved to replicate within the unique features of their environments. The key question remains: What role do these viruses play in the ecology and evolution of the environments in which they replicate? Undoubtedly, as more virologists examine the role of viruses in subsurface environments, new insights will emerge.

Intraspecific antagonism through viral toxin encoded by chronic Sulfolobus spindle-shaped virus

Virus-host interactions evolve along a symbiosis continuum from antagonism to mutualism. Long-term associations between virus and host, such as those in chronic infection, will select for traits that drive the interaction towards mutualism, especially when susceptible hosts are rare in the population. Virus-host mutualism has been demonstrated in thermophilic archaeal populations where Sulfolobus spindle-shaped viruses (SSVs) provide a competitive advantage to their host Sulfolobus islandicus by producing a toxin that kills uninfected strains. Here, we determine the genetic basis of this killing phenotype by identifying highly transcribed genes in cells that are chronically infected with a diversity of SSVs. We demonstrate that these genes alone confer growth inhibition by being expressed in uninfected cells via a Sulfolobus expression plasmid. Challenge of chronically infected strains with vector-expressed toxins revealed a nested network of cross-toxicity among divergent SSVs, with both broad and specific toxin efficacies. This suggests that competition between viruses and/or their hosts could maintain toxin diversity. We propose that competitive interactions among chronic viruses to promote their host fitness form the basis of virus-host mutualism. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Genome Sequence of a Thermoacidophilic Methanotroph Belonging to the Verrucomicrobiota Phylum from Geothermal Hot Springs in Yellowstone National Park: A Metagenomic Assembly and Reconstruction

Verrucomicrobiotal methanotrophs are thermoacidophilic methane oxidizers that have been isolated from volcanic and geothermal regions of the world. We used a metagenomic approach that entailed obtaining the whole genome sequence of a verrucomicrobiotal methanotroph from a microbial consortium enriched from samples obtained from Nymph Lake (89.9 °C, pH 2.73) in Yellowstone National Park in the USA. To identify and reconstruct the verrucomicrobiotal genome from Illumina NovaSeq 6000 sequencing data, we constructed a bioinformatic pipeline with various combinations of de novo assembly, alignment, and binning algorithms. Based on the marker gene (pmoA), we identified and assembled the Candidatus Methylacidiphilum sp. YNP IV genome (2.47 Mbp, 2392 ORF, and 41.26% GC content). In a comparison of average nucleotide identity between Ca. Methylacidiphilum sp. YNP IV and Ca. Methylacidiphilum fumariolicum SolV, its closest 16S rRNA gene sequence relative, is lower than 95%, suggesting that Ca. Methylacidiphilum sp. YNP IV can be regarded as a different species. The Ca. Methylacidiphilum sp. YNP IV genome assembly showed most of the key genes for methane metabolism, the CBB pathway for CO₂ fixation, nitrogen fixation and assimilation, hydrogenases, and rare earth elements transporter, as well as defense mechanisms. The assembly and reconstruction of a thermoacidophilic methanotroph belonging to the Verrucomicrobiota phylum from a geothermal environment adds further evidence and knowledge concerning the diversity of biological methane oxidation and on the adaptation of this geochemically relevant reaction in extreme environments.

The biology of thermoacidophilic archaea from the order Sulfolobales

Thermoacidophilic archaea belonging to the order Sulfolobales thrive in extreme biotopes, such as sulfuric hot springs and ore deposits. These microorganisms have been model systems for understanding life in extreme environments, as well as for probing the evolution of both molecular genetic processes and central metabolic pathways. Thermoacidophiles, such as the Sulfolobales, use typical microbial responses to persist in hot acid (e.g. motility, stress response, biofilm formation), albeit with some unusual twists. They also exhibit unique physiological features, including iron and sulfur chemolithoautotrophy, that differentiate them from much of the microbial world. Although first discovered >50 years ago, it was not until recently that genome sequence data and facile genetic tools have been developed for species in the Sulfolobales. These advances have not only opened up ways to further probe novel features of these microbes but also paved the way for their potential biotechnological applications. Discussed here are the nuances of the thermoacidophilic lifestyle of the Sulfolobales, including their evolutionary placement, cell biology, survival strategies, genetic tools, metabolic processes and physiological attributes together with how these characteristics make thermoacidophiles ideal platforms for specialized industrial processes.

High Potential for Biomass-Degrading Enzymes Revealed by Hot Spring Metagenomics

Enzyme stability and activity at elevated temperatures are important aspects in biotechnological industries, such as the conversion of plant biomass into biofuels. In order to reduce the costs and increase the efficiency of biomass conversion, better enzymatic processing must be developed. Hot springs represent a treasure trove of underexplored microbiological and protein chemistry diversity. Herein, we conduct an exploratory study into the diversity of hot spring biomass-degrading potential. We describe the taxonomic diversity and carbohydrate active enzyme (CAZyme) coding potential in 71 publicly available metagenomic datasets from 58 globally distributed terrestrial geothermal features. Through taxonomic profiling, we detected a wide diversity of microbes unique to varying temperature and pH ranges. Biomass-degrading enzyme potential included all five classes of CAZymes and we described the presence or absence of genes encoding 19 glycosyl hydrolases hypothesized to be involved with cellulose, hemicellulose, and oligosaccharide degradation. Our results highlight hot springs as a promising system for the further discovery and development of thermo-stable biomass-degrading enzymes that can be applied toward generation of renewable biofuels. This study lays a foundation for future research to further investigate the functional diversity of hot spring biomass-degrading enzymes and their potential utility in biotechnological processing.

The network structure and eco-evolutionary dynamics of CRISPR-induced immune diversification

As a heritable sequence-specific adaptive immune system, CRISPR-Cas is a powerful force shaping strain diversity in host-virus systems. While the diversity of CRISPR alleles has been explored, the associated structure and dynamics of host-virus interactions have not. We explore the role of CRISPR in mediating the interplay between host-virus interaction structure and eco-evolutionary dynamics in a computational model and compare the results with three empirical datasets from natural systems. We show that the structure of the networks describing who infects whom and the degree to which strains are immune, are respectively modular (containing groups of hosts and viruses that interact strongly) and weighted-nested (specialist hosts are more susceptible to subsets of viruses that in turn also infect the more generalist hosts with many spacers matching many viruses). The dynamic interplay between these networks influences transitions between dynamical regimes of virus diversification and host control. The three empirical systems exhibit weighted-nested immunity networks, a pattern our theory shows is indicative of hosts able to suppress virus diversification. Previously missing from studies of microbial host-pathogen systems, the immunity network plays a key role in the coevolutionary dynamics.

Duplication of leucyl-tRNA synthetase in an archaeal extremophile may play a role in adaptation to variable environmental conditions

Aminoacyl-tRNA synthetases (aaRSs) are ancient enzymes that play a fundamental role in protein synthesis. They catalyze the esterification of specific amino acids to the 3'-end of their cognate tRNAs and therefore play a pivotal role in protein synthesis. Although previous studies suggest that aaRS-dependent errors in protein synthesis can be beneficial to some microbial species, evidence that reduced aaRS fidelity can be adaptive is limited. Using bioinformatics analyses, we identified two distinct leucyl-tRNA synthetase (LeuRS) genes within all genomes of the archaeal family Sulfolobaceae. Remarkably, one copy, designated LeuRS-I, had key amino acid substitutions within its editing domain that would be expected to disrupt hydrolytic editing of mischarged tRNA^Leu and to result in variation within the proteome of these extremophiles. We found that another copy, LeuRS-F, contains canonical active sites for aminoacylation and editing. Biochemical and genetic analyses of the paralogs within Sulfolobus islandicus supported the hypothesis that LeuRS-F, but not LeuRS-I, functions as an essential tRNA synthetase that accurately charges leucine to tRNA^Leu for protein translation. Although LeuRS-I was not essential, its expression clearly supported optimal S. islandicus growth. We conclude that LeuRS-I may have evolved to confer a selective advantage under the extreme and fluctuating environmental conditions characteristic of the volcanic hot springs in which these archaeal extremophiles reside.

Diversified local CRISPR-Cas immunity to viruses of Sulfolobus islandicus

The population diversity and structure of CRISPR-Cas immunity provides key insights into virus-host interactions. Here, we examined two geographically and genetically distinct natural populations of the thermophilic crenarchaeon Sulfolobus islandicus and their interactions with Sulfolobus spindle-shaped viruses (SSVs) and S. islandicus rod-shaped viruses (SIRVs). We found that both virus families can be targeted with high population distributed immunity, whereby most immune strains target a virus using unique unshared CRISPR spacers. In Kamchatka, Russia, we observed high immunity to chronic SSVs that increases over time. In this context, we found that some SSVs had shortened genomes lacking genes that are highly targeted by the S. islandicus population, indicating a potential mechanism of immune evasion. By contrast, in Yellowstone National Park, we found high inter- and intra-strain immune diversity targeting lytic SIRVs and low immunity to chronic SSVs. In this population, we observed evidence of SIRVs evolving immunity through mutations concentrated in the first five bases of protospacers. These results indicate that diversity and structure of antiviral CRISPR-Cas immunity for a single microbial species can differ by both the population and virus type, and suggest that different virus families use different mechanisms to evade CRISPR-Cas immunity. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Microbial biogeography of 925 geothermal springs in New Zealand

Geothermal springs are model ecosystems to investigate microbial biogeography as they represent discrete, relatively homogenous habitats, are distributed across multiple geographical scales, span broad geochemical gradients, and have reduced metazoan interactions. Here, we report the largest known consolidated study of geothermal ecosystems to determine factors that influence biogeographical patterns. We measured bacterial and archaeal community composition, 46 physicochemical parameters, and metadata from 925 geothermal springs across New Zealand (13.9–100.6 °C and pH < 1–9.7). We determined that diversity is primarily influenced by pH at temperatures <70 °C; with temperature only having a significant effect for values >70 °C. Further, community dissimilarity increases with geographic distance, with niche selection driving assembly at a localised scale. Surprisingly, two genera (Venenivibrio and Acidithiobacillus) dominated in both average relative abundance (11.2% and 11.1%, respectively) and prevalence (74.2% and 62.9%, respectively). These findings provide an unprecedented insight into ecological behaviour in geothermal springs, and a foundation to improve the characterisation of microbial biogeographical processes.

Power et al. catalogue the microbial biodiversity and physicochemistry of around 1000 hotsprings across New Zealand, providing insights into the ecological conditions that drive community assembly in these ecosystems.

Structured Populations of Sulfolobus acidocaldarius with Susceptibility to Mobile Genetic Elements

The impact of a structured environment on genome evolution can be determined through comparative population genomics of species that live in the same habitat. Recent work comparing three genome sequences of Sulfolobus acidocaldarius suggested that highly structured, extreme, hot spring environments do not limit dispersal of this thermoacidophile, in contrast to other co-occurring Sulfolobus species. Instead, a high level of conservation among these three S. acidocaldarius genomes was hypothesized to result from rapid, global-scale dispersal promoted by low susceptibility to viruses that sets S. acidocaldarius apart from its sister Sulfolobus species. To test this hypothesis, we conducted a comparative analysis of 47 genomes of S. acidocaldarius from spatial and temporal sampling of two hot springs in Yellowstone National Park. While we confirm the low diversity in the core genome, we observe differentiation among S. acidocaldarius populations, likely resulting from low migration among hot spring “islands” in Yellowstone National Park. Patterns of genomic variation indicate that differing geological contexts result in the elimination or preservation of diversity among differentiated populations. We observe multiple deletions associated with a large genomic island rich in glycosyltransferases, differential integrations of the Sulfolobus turreted icosahedral virus, as well as two different plasmid elements. These data demonstrate that neither rapid dispersal nor lack of mobile genetic elements result in low diversity in the S. acidocaldarius genomes. We suggest instead that significant differences in the recent evolutionary history, or the intrinsic evolutionary rates, of sister Sulfolobus species result in the relatively low diversity of the S. acidocaldarius genome.