Nutritional niches of endemic, facultatively anaerobic heterotrophs from an isolated Antarctic terrestrial hydrothermal refugium elucidated through metagenomics

doi:10.21203/rs.3.rs-4805162/v1

Download PDF

Research Article

Nutritional niches of endemic, facultatively anaerobic heterotrophs from an isolated Antarctic terrestrial hydrothermal refugium elucidated through metagenomics

https://doi.org/10.21203/rs.3.rs-4805162/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background

Tramway Ridge, a geothermal Antarctic Specially Protected Area (1) (elevation 3340 m) located near the summit of Mount Erebus, is home to a unique community composed of cosmopolitan surface-associated micro-organisms and abundant, poorly understood subsurface-associated microorganisms (2–5). Here, we use shotgun metagenomics to compare the functional capabilities of this community to those found elsewhere on Earth and to infer endemism and metabolic capabilities of abundant subsurface taxa.

Results

We found that the functional potential in this community is most similar to that found in terrestrial hydrothermal environments (hot springs, sediments) and that the two dominant organisms in the subsurface are primarily endemic. They were found to be facultative anaerobic heterotrophs that likely share a pool of nitrogenous organic compounds while specializing in different carbon compounds.

Conclusions

Metagenomic insights have provided a detailed understanding of the microbe-based ecosystem found in geothermally heated fumaroles at Tramway Ridge. This approach enabled us to compare Tramway Ridge with other microbial systems, identify endemic taxa and elucidate the key metabolic pathways that may enable specific organisms to dominate the ecosystem.

endemic

hydrothermal

geothermal

volcano

fumarole

metabolism

sediments

Mt. Erebus, Victoria Land, Antarctica, is the highest, most southern, isolated geothermal feature on the planet (6). Since its original submarine eruption ~ 1.3Ma and through its multiple building phases (7, 8), an ever-present Antarctic Circumpolar current has kept Mt. Erebus relatively isolated from volcanoes found elsewhere on Earth (9). Its biogeographical isolation is supported by previous research on the soil microbial communities at Tramway Ridge, a small geothermal feature on the summit plateau, where unique, deep-branching, and potentially endemic lineages of Bacteria and Archaea were found within highly thermally stratified fumaroles (Fig. 1) (2–4, 10, 11).

The fumaroles of Tramway Ridge differ in many ways from other terrestrial hydrothermal features. They are characterized by hot (65°C) CO₂-rich steam venting through slightly alkaline (pH 8) hydrothermally altered mineral soils (12, 13). Oxygen levels within fumarolic sediments are approximately 1 mg L^− 1 (roughly 25% saturated for the temperature) indicative of subsurface hypoxia (14, 15). All Erebus hydrothermal features are driven by a phonolite magmatic source (Sims et al. 2021; Noell et al. 2022). Off-gasses tend to be composed primarily of steam, with low levels of sulfur and elevated concentrations of methane, hydrogen and CO (16, 17). The high pH, moderate temperature, lack of standing water, and hypoxia distinguish them from other well-studied terrestrial hot springs and mud volcanoes and provide a unique geothermally-driven range of micro-environments for resident microbiota. Steam is vented through concentrated hotspots, generating steep temperature (-20 to 65°C) and pH (3.5–8) gradients over less than a meter that are major determinants of the composition of the thick cyanobacterial mats and associated microbial communities observed on the surface (3, 4). Within fumaroles, the temperature is a relatively constant 65°C, even at > 5 cm depth, but can decrease suddenly to less than 20°C and stay low for 24 hours or more at a time (3, 15). Beyond its unique physical characteristics, the hot fumarolic soils of Tramway Ridge also have an extremely low total C:N ratio (ranging from 1:3 to 3:1) (2, 3), which suggests that the microbial community experiences continual carbon-limitation relative to nitrogen (18).

The isolation of Tramway Ridge, its unique geochemical environment, and microbiota make it an exciting site for evaluating potential endemism and for identifying novel metabolic pathways. Driven by the novelty of the taxa encountered there and the lack of information regarding their metabolic potential, we launched a detailed metagenomic study of high-altitude Antarctic fumaroles. First, we used functional profiles to contextualize the functional repertoire of this particular community with respect to other types of microbial communities. Second, we improved on our previous effort to identify endemic taxa (2), where we relied on matching partial 16S rRNA gene amplicon sequences with database entries by developing a novel metric that is based on synonymous polymorphisms within reconstructed environmental populations to define/circumscribe endemic taxa. Finally, we sought to elucidate possible novel metabolic processes encoded by abundant and endemic taxa that are specifically localized to the fumarolic subsurface (depths > 2 cm).

Sample collection

Soil samples were collected within the Tramway Ridge Antarctic Specially Protected Area (ASPA 130) in February 2009 from two sites (site A – 77° 31.103' S, 167° 6.682' S and site B – 77° 31.306' S, 167° 6.668' E). Sites were chosen based on measuring a surface temperature of 65°C with a stainless steel Checktemp1 temperature probe (Hanna Instruments, Rhode Island, USA) sterilized with 70% ethanol immediately prior to use. Temperature measurements were repeated for each layer sampled. Surface soil crusts were carefully set aside prior to collecting samples. Samples were collected by carefully removing 2 cm of soil in an approximately 5 cm x 5 cm square area using an autoclaved stainless steel spatula wiped with 70% ethanol just prior to sampling. Soil was placed into a fresh 50 mL Falcon tube and immediately frozen at -20°C. Sampling continued with the collection of a second layer (2–4 cm depth).

DNA extraction, library preparation and sequencing

DNA was extracted from samples using a modified CTAB (cetyltrimethylammonium bromide) bead-beating protocol (19) and quantified using the Quant-IT dsDNA BR Assay Kit (Invitrogen, Carlsbad, CA, USA). A portion of extracted metagenomic DNA was frozen and sent to sequencing facilities at the University of California-Los Angeles (USA) and the University of Waikato (Hamilton, New Zealand). At each location, samples were processed and sequenced using standard protocols for the 454-Ti platform (Roche 454 Life Sciences, Branford, CT, USA). Additional DNA fractions were sent to the sequencing facility at University of California, Davis, where paired-end libraries were prepared and sequenced using the Illumina Hi-Seq platform (Illumina, San Diego, CA, USA).

Metagenome assembly and binning

Assembly of metagenomic data was carried out using appropriate assemblers for the two different sequencing platforms. 454 reads were assembled from original sff files using Newbler (Roche). Reads from each dataset were assembled independently and were pooled for an additional assembly. Newbler assembly used default parameters except that minimum overlap was set to 100 nucleotides and overlap identity was set to 98% identity. Paired-end Illumina reads were pre-processed by removing any paired-end set for which identifying tags had at least one mismatch or for which the paired-end tags were not identical. Identifying tags were removed, and reads were end-trimmed at the first incidence of a quality score below 14. Any reads that were shorter than 20 nucleotides after tag removal and end-trim were removed from the dataset. 454 datasets were assembled in Newbler (-mi 98 -ml 100 -minlen 45 -a 500 -l 2000). Illumina assemblies were carried out using two sets of options (--careful –only-assembler; --cov-cutoff auto --careful -k 25,55,65,75 --only-assembler) in Spades (20), a single setting (--meta --only-assembler) in MetaSpades (21) and a single setting (--k-min 27 --k-max 127 --k-step 10 --min-contig-len 500 --prune-level 3 --no-mercy --min-count 3 --no-local) in Megahit (22). BBMap (https://sourceforge.net/projects/bbmap/) was used to map all individual read sets (454 and Illumina) against all assemblies. Large (> 2 kb) contigs and scaffolds were clustered into Metagenome-assembled genomes by oligonucleotide frequency and read coverage using Maxbin 2 (23) and Metabat 2 (24). Redundant bins were subsequently dereplicated and evaluated using dRep (25) with a completeness cutoff of 40%, contamination cutoff of 10% and a minimum genome size of 200kb. Final MAGs (Table 1) were classified using GTDB-Tk v2.1 (26) with genome database release R214 (27) and annotated with the NCBI Prokaryotic Genome Annotation Pipeline (28). MAGs discussed in the text were additionally annotated with eggNOG mapper (29), Cytochrome c oxidases were

Table 1

Species-level representatives of metagenome-assembled genomes (MAGs) from Tramway Ridge fumarolic soils. Lineage was assigned using GTDB-Tk v.2.1 (26) with Release R214 of the Genome Taxonomy Database (27). Completeness and contamination estimates were calculated in CheckM (50). Assembly quality follows the recommendations outlined previously (55). For a comprehensive list of genome quality attributes for all binned strains, see supplementary Table 3.
Proposed name	species group	# of strains in species group	Division	assembly quality	GC	Genome size (bp)
Blastocatellia bacterium ERB_27	27	2	Acidobacteria	high	50.9	2698697
Pyrinomonas sp. ERB_32	32	1	Acidobacteria	high	60	3512547
Acidimicrobiia bacterium ERB_23	23	3	Actinobacteria	medium	55.1	1963147
Acidimicrobiia bacterium ERB_8	8	1	Actinobacteria	medium	69.4	1352013
Thermoleophilia bacterium ERB_19	19	3	Actinobacteria	high	69.5	2326210
Armatimonadetes bacterium ERB_24	24	1	Armatimonadetes	medium	60.8	1563163
Armatimonadetes bacterium ERB_33	33	1	Armatimonadetes	high	60.5	3260803
Armatimonadetes bacterium ERB_34	34	2	Armatimonadetes	high	58	2190352
Armatimonadetes bacterium ERB_6	6	3	Armatimonadetes	high	61.1	2058555
Chitinophagaceae sp. ERB_2	2	1	Bacteroidetes	high	40.5	2908196
Chitinophagaceae sp. ERB_3	3	2	Bacteroidetes	medium	39.3	2810122
Candidatus GAL15 bacterium ERB_18	18	2	candidate division GAL15	high	70.1	1986527
Candidatus Dadabacteria bacterium ERB_12	12	1	Candidatus Dadabacteria	medium	44.2	3057813
Candidatus Fervidibacteria erebusii ERB_15	15	1	Candidatus Fervidibacteria	high	56.3	2795547
Candidatus Lakebacteria bacterium ERB_C1	C1	3	Candidatus Lakebacteria	medium	27.3	481934
Chloroflexi bacterium ERB_10	10	1	Chloroflexi	medium	64.7	2882324
Chloroflexi bacterium ERB_11	11	1	Chloroflexi	medium	68.5	2692019
Chloroflexi bacterium ERB_20	20	3	Chloroflexi	high	63.6	3392850
Candidatus Nitrocaldera therma ERB_22	22	2	Chloroflexi	high	65.1	2707038
Chloroflexi bacterium ERB_25	25	1	Chloroflexi	Low-quality draft	64.8	1709950
Chloroflexi bacterium ERB_7	7	1	Chloroflexi	medium	70.2	1489492
Chloroflexi bacterium ERB_9	9	2	Chloroflexi	medium	62.8	3793921
Thermoflexus sp. ERB_21	21	3	Chloroflexi	high	69.6	2384148
Leptolyngbya sp. ERB_1	1	2	Cyanobacteria	medium	47.2	4634638
Mastigocladus sp. ERB_26	26	2	Cyanobacteria	high	41.2	5965913
Meiothermus sp. ERB_29	29	2	Deinococcus-Thermus	medium	66.3	2588060
Meiothermus sp. ERB_30	30	2	Deinococcus-Thermus	medium	61.5	3765706
Meiothermus sp. ERB_31	31	1	Deinococcus-Thermus	medium	68.5	2746960
Thermus sp. ERB_17	17	5	Deinococcus-Thermus	high	65.2	1806162
Ignavibacteria bacterium ERB_28	28	1	Ignavibacteriae	high	55.2	3309624
Nitrospiraceae bacterium ERB_14	14	2	Nitrospirae	medium	56.2	3371732
Gemmataceae bacterium ERB_16	16	2	Planctomycetes	high	60.1	3796795
Rhodanobacteraceae bacterium ERB_4	4	1	Proteobacteria	medium	62.8	1332200
Candidatus Australlarchaeum erebusii ERB_5	5	1	Thaumarchaeota	medium	63.6	1213694
Nitrososphaera sp. ERB_13	13	1	Thaumarchaeota	medium	56.6	1122400

checked with the HCO classifier (30), hydrogenases were checked with HydDB (31), and CAZymes were checked with dbCAN3 (32).

Endemicity Index Calculation

Quality-trimmed Illumina reads were mapped to each MAG using BBmap and further filtered using a hard 97% identity cutoff where identity = number of matches/length of alignment, with the additional requirement that at least 50 nucleotides mapped. Raw diversity was compiled using Samtools mpileup using option: -d 1000000 (33) for each MAG and each dataset independently. SNPs were determined using Varscan2 pileup2snp (options: --min-var-freq 0.01 --p-value 0.05) (34) and further filtered using the Benjamini-Hochberg multiple testing correction FDR = 0.01. SnpEff (35) was used to classify SNPs as synonymous with gff files produced with Prodigal (36). The density of synonymous SNPs (D_SynSNP) for a MAG was calculated as the number of Synonymous SNPs (N_SynSNP) divided by MAG length in Mb (L_MAG): D_SynSNP = N_SynSNP / L_MAG. Because sensitivity increases as read depth increases, D_SynSNP was corrected for read coverage of the MAG (C_MAG) which was calculated as the number of reads mapped (N_M) divided by MAG length (L_MAG): C_MAG = N_M / L_MAG. Endemicity Index (EI) was then the density of synonymous SNPs (D_SynSNP) divided by read coverage (C_MAG): EI = D_SynSNP / C_MAG. EI was calculated for each MAG/dataset combination only if N_SynSNP ≥ 5 and reported EI values are the average of calculated values over four Illumina read sets (SRA accessions SRR6519253, SRR6519254,SRR6519255, SRR6519256).

Functional profile comparisons

Pfam (37) profiles were used to compare functional similarity between the metagenomic assembly from Tramway Ridge and publicly available metagenomic assemblies in the Integrated Microbial Genomes (IMG) database (38). A list of all assembled, published, and “unrestricted” environmental metagenomic datasets available through IMG was downloaded on 31 March, 2023. IMG-generated pfam profiles (counts of pfams present in metagenomic assembly) were downloaded for each available metagenome, resulting in 7652 total pfam profiles. The number of metagenomic datasets was reduced by removing datasets with total assembly length less 5x10⁷ or greater than 1x10⁹ bases. This subset was further reduced by removing metagenomic datasets with fewer than 3500 or greater than 7500 unique pfams. These filtering criteria resulted in a dataset of 4513 publicly sourced metagenomic datasets for comparative analysis. Profiles were analyzed in R 4.2.3. Jaccard dissimilarity was calculated using the vegdist() function from vegan 2.6-4 (39) based on presence/absence of pfams. Principle coordinate analysis (PCoA) was carried out using the pcoa() function from Version 5.7-1 of ape (40) and plotted in three dimensions using the plot3d() function from rgl 1.2.1 (41). t-distributed stochastic neighbor embedding (tSNE) was calculated using the Rtsne() function from Rtsne 0.16 (42) and plotted with ggplot2 3.4.3 (43). Rtsne settings were as follows: Rtsne(X, dims = 2,initial_dims = 5, perplexity = 300, theta = 0.5, check_duplicates = FALSE, pca = TRUE, partial_pca = FALSE, max_iter = 5000000, verbose = getOption("verbose",FALSE), is_distance = TRUE, Y_init = NULL, pca_center = TRUE, pca_scale = FALSE, normalize = TRUE, stop_lying_iter = 500000, mom_switch_iter = 500000, momentum = 0.5, final_momentum = 0.8, eta = 100, exaggeration_factor = 12).

Phylogenomic analysis

A collection of reference genomes for comparative phylogenomic analysis was assembled from representative species defined in release 214 (April, 2023) of the Genome Taxonomy Database (27). GTDB classifications of MAGs from the current study were used to select from GTDB species representatives for an informative tree. For instance, in the case where a MAG was classified into an order but not a family, one representative taxon from each family was included, the phylogeny was calculated, and if the MAG was associated with a particular family, the process was repeated using genus representatives from that family. These genomes were supplemented with additional genomes from thermophilic environments (44, 45), thermophilic nitrifying enrichment cultures (46, 47), additional Nitrospirota (48) and additional Nitrososphaeria (49). All genomes were downloaded and processed using CheckM (50) to generate concatenated alignments of 34 universal marker genes (43 marker HMMs). To be included in phylogenomic reconstruction, reference genomes were required to be at least 60% complete with less than 5% contamination and to have at least 4500 ungapped characters in the concatenated alignment (6988 total positions). All genomes used for phylogenetic analysis are listed in Supplementary Table 8. Phylogenetic reconstruction with IQ-TREE 2 (51) included model selection with ModelFinder (52) and calculation of bootstraps with UFboots (53).

We found that the functional profiles from Tramway Ridge metagenomes resembled those from other thermally-influenced environments, in particular those from terrestrial hydrothermal systems. To learn this, we compared functional profiles of assembled Tramway Ridge metagenomes to 4513 publicly available and assembled environmental metagenomes broadly categorized as terrestrial hydrothermal / non-hydrothermal, freshwater, and marine hydrothermal / non-hydrothermal (Supplementary Table 1). Functional profiles for metagenomes were constructed based on the presence-absence of Pfam protein family domains (54), and dissimilarities were calculated using the Jaccard index. A post-hoc Tukey’s HSD test (Fig. 2a) comparing Jaccard dissimilarity grouped according to the broad categories listed above showed that metagenomes from Tramway Ridge were most similar to terrestrial hydrothermal environments (Tukey category a) and least similar to metagenomes from marine non-hydrothermal and freshwater environments (Tukey category d). We continued with principal coordinate analysis (PCoA) and t-distributed Stochastic Neighbor Embedding (tSNE) visualizations to explore relationships that may not have been clear from grouped pairwise comparisons (Fig. 2bc). In both, Tramway Ridge microbial communities clustered loosely with microbial communities sourced from both terrestrial and marine hydrothermal environments to the exclusion of non-hydrothermal environments. T-SNE visualization (Fig. 2c) recovered several distinct but associated clusters from hydrothermal systems, one of which was composed exclusively of profiles from Tramway Ridge.

Metagenome-assembled genomes (MAGs) were constructed from DNA extracted from soils at two 65°C fumaroles at Tramway Ridge (supplementary Table 2). We recovered 63 MAGs at the strain level (99% average nucleotide identity, ANI), which clustered into 35 species-level representatives (> 96.5% average nucleotide identity) (Fig. 3, Table 1, supplementary Table 3). A total of 16 MAGs met the criteria for high-quality drafts, and 18 MAGs met the criteria for medium-level drafts (55); these include nearly complete (75–96%) genome bins for a novel order of Archaea within the Nitrososphaeria (Fig. 3b) and novel lineages of Bacteria including Armatimonadota, Chloroflexota, Actinobacteriota, and Candidate division CSP1-3. Additional MAGs of interest include those belonging to the Candidatus Patescibacteria phylum (Fig. 3c), Mastigocladus genus (Fig. 3d) and Candidatus Fervidibacter genus (Fig. 3e). MAGs were classified using GTDB-Tk v.2.1 (26) with Release R214 of the Genome Taxonomy Database (27), and novel taxa of note were given proposed names (Table 1, supplementary Table 3).

Tramway Ridge occupies a unique location on Earth that makes it ideal for the study of endemism due to its highly isolated nature and the geothermal enrichment that prevents intrusion of nearby psychrophilic and mesophilic microorganisms. Here, we quantified endemism for each of the MAGs we generated. We defined the degree of endemism of a given taxon as being proportional to the in situ diversity of that taxon, reasoning that endemic species have had an extended opportunity to diversify on-site, as opposed to recent arrivals that would be subject to founder effects. We developed a simple metric, which we called the endemicity index (EI) (Fig. 3F) that measured the frequency of synonymous mutations in a MAG accounting for read depth. High values (e.g. 10^− 2) of EI indicate high diversity whereas low values (e.g. 10^− 6) indicate low diversity.

The lowest median EI value (1 x 10^− 5) for a single species was calculated for ERB_26, (Fig. 3f), which represented two individual strains of a cyanobacterium belonging to the Fischerella / Mastigocladus genus (supplementary Table 3) and dominated the near-surface (5–10% of near-surface reads). The highest median EI value (1 x 10^− 2) for a single species was calculated for ERB_C1 (Fig. 3f), which represented three individual strains of Candidatus Patescibacteria (supplementary Table 3). Ca. Lakebacteria bacterium ERB_C1 classified into the HRBIN35 genus in the HR35 family (Fig. 3c). Based on the percentage of reads recruited, it was slightly more abundant in the near-surface (0–2 cm depth) than in the subsurface (> 2 cm, Fig. 3g). Due to this distribution pattern and the fact that an accelerated evolutionary process is a hallmark feature of members of this lineage (27, 56, 57), we did not pursue further analysis of this genome, which would be best accomplished in a comprehensive comparison with other Ca. Patescibacteria.

Two MAGs, ERB_15_1 and ERB_5_1, with relatively high EI values (3 x 10^− 3 and 1 x 10^− 3, respectively) dominated the subsurface community at depths greater than 2 cm (Fig. 3f, Fig. 3g). ERB_5_1 shared less than 40% average amino acid identity (AAI) with any other Nitrososphaeria (formerly Thaumarchaeota) and formed a singular deep branch of the Nitrososphaeria in phylogenetic trees (FIgure 3b). ERB_15_1 shared 69% AAI and 77% ANI with Candidatus Fervidibacteria sacchari, establishing it as a novel species of this genus (Fig. 3c). Here we propose the names Candidatus Fervidibacteria erebusii (ERB_15_1) and Candidatus Australlarchaeum erebusii (ERB_5_1) to reflect their first observation and endemism at Mt. Erebus. We picked these two MAGs for further analyses based on their high EI values, subsurface abundance, and taxonomic novelty.

Candidatus A. erebusii (ERB_5_1) branches near the base of the Nitrososphaeria (Fig. 3b) and, like other deeply diverging lineages, it lacked marker genes for cobalamin biosynthesis (cob/cbi/bluB) and ammonia oxidation (archaeal amoABC, HAO). Consistent with previous observations of deeply diverging Nitrososphaeria, the genome of Ca. A. erebusii encoded genes for the beta oxidation of fatty acids and peptide degradation and several genes that may be used for amino acid-based oxidation (Supplementary Table 5). Ca. A. erebusii encoded aminotransferases and glutamate dehydrogenase that could be used to produce NAD(P)H by converting glutamate to 2-oxoglutaric acid. A 2-oxoglutarate:ferredoxin oxidoreductase could then convert the 2-oxoglutaric acid to succinyl-CoA, which could then be converted to succinate with the production of ATP via ADP-forming succinate-CoA ligase (Supplementary Table 5). Unlike other deep-branching Nitrososphaeria, Ca. A. erebusii also additionally encoded two aa3-type (low-affinity) cytochrome C-oxidases (Supplementary Table 5). We identified a putative aerobic carbon monoxide dehydrogenase (CODH) cluster of genes (coxMLS) that may confer the ability to use CO as an additional energy source. However, the large subunit of the CODH in Ca. A. erebusii contained a hallmark motif, AYXGAGR, of type II CODH, which oxidizes CO only very slowly and possibly incidentally (58). No predicted carbon fixation pathways were identified.

Candidatus Fervidibacteria erebusii (ERB_15_1) shared 69% AAI and 77% ANI with Candidatus Fervidibacteria sacchari, establishing it as a novel species of this genus (Fig. 3c). Based on conservative annotations, the Ca. F. erebusii genome encoded 144 CAZymes, including 69 glycosyl hydrolase (GH) domains spanning 32 GH families (Supplementary Table 6, Supplementary Table 7), and an impressive 15 variants of the unusual and poorly characterized GH109 family. Furthermore, when we group the Ca. F. erebusii CAZymes by substrate utilization capacity, the GH109s belong to the largest cohort (35.25%) which appear to be involved in depolymerisation cascades of N-acetylhexosamide glycans such as chitins and chitosan (Supplementary Table 7). Other major polysaccharide utilization cohorts include cellulose and hemicellulose (19.8%), pectin (14.3%) and ß-glucans (19.8%), with minor cohorts for ɑ-glucans, ɑ-mannans and ɑ-fucans. Ca. Fervidibacter genomes, including Ca. F. erebusii additionally encoded numerous aminotransferases, glutamate dehydrogenase, indole-pyruvate:ferrodoxin oxidoreductase and ADP-forming succinate-CoA ligase (Supplementary Table 6) which may explain the ability of Ca. Fervidibacter sacchari to grow solely on casamino acids (59). The Ca. F. erebusii genome also encoded Group 2a NiFe hydrogenases which are high-affinity hydrogenases that enable the survival of soil heterotrophs under carbon limitation (60) and enable hydrogenotrophic growth of autotrophic nitrite-oxidizing bacteria under nitrite limitations (61). The Ca. F. erebusii genome also encoded a cluster of putative CODH genes (coxMLS), which contained the hallmark motif, AYXGAGR, of type II CODH, indicating a low, incidental activity (58). No carbon fixation pathways were found.

Metagenome Functional Profiles

Mt. Erebus has evolved over time, starting with seafloor rifting and growing as a subaerial volcano into a modern-day stratovolcano (7, 8). The current conditions at Tramway Ridge differ remarkably from other extant non-Antarctic terrestrial and marine hydrothermal systems, being primarily driven by the unique phonolite magmatic source resulting in alkaline fumaroles with low sulphidic content. It is unknown how long geothermal features such as the fumaroles at Tramway Ridge have been present on Mt. Erebus; however, volcanism was once widespread across the Ross Island massif (6). Extant geothermal features may represent a once widespread ecosystem of similar features. Given this complex history of Mt. Erebus, we questioned whether the Tramway Ridge community has retained a legacy signature of its origin on the seafloor or if the modern-day Tramway Ridge community better resembles other terrestrial hydrothermal sites.

To answer this question, we compared functional profiles of assembled Tramway Ridge metagenomes to a large set of publicly available metagenomes. Metagenomes from Tramway Ridge were distinct from all others (Fig. 2c) but showed the most similarity to terrestrial hydrothermal environments such as hot springs and associated sediments. This indicates that the legacy of a sea floor origin is less important than the extant conditions at Tramway ridge in defining the microbial community. This is also reflected in the prevalence of taxa at Tramway Ridge that are found in other terrestrial hydrothermal environments but not in seafloor hydrothermal systems, such as the genera Mastigocladus, Caldithermus, Meiothermus and phyla such as the Chloroflexota, Armatimonadota and Actinobacteriota.

The unique nature of Tramway Ridge metagenomic profiles likely reflects the site’s relatively low diversity (2) and the type of functional profile analysis used here. First, pfams were used, which are very coarse functional units. Second, gene abundance was not available for most metagenomes, so we were forced to use presence-absence as data. Third, we found that the use of pfam presence-absence was sensitive to metagenome assembly size as well as the total number of unique pfams detected. In combination, these features had the potential to obfuscate any meaningful relationships. Our analysis was optimized to be as permissive as possible, allowing as many metagenomes into the analysis while still comparing like against like. Using this strategy, we were able to compare Tramway Ridge metagenomic datasets to 4513 out of 7652 publicly available and unrestricted metagenomes (see 10.5281/zenodo.10928974).

Within continental Antarctica, only three known surface-expressed active geothermal areas exist (Mt. Rittmann, Mt. Melbourne and Mt Erebus), each of which is separated by vast ice fields (15). It is thought that intercontinental transport of microbe-bearing particulates into Antarctica occurs much less frequently than intracontinental transport (62), suggesting that the introduction of exogenous microbes is relatively infrequent. Recent studies have used database searches to identify endemic species (63), and a lack of sequence identity to database entries has been used in the past to suggest that novel sequences indicate endemism (2). However, defining endemism based on whether sequence matches exist within a database can be problematic as this definition is sensitive to database composition. Conclusions drawn may not withstand the inevitable growth in database size and the diversity it holds. For this study, we took an alternative approach and defined the degree of endemism of a given taxon as being proportional to the in situ diversity of that taxon. For these inferences, we assumed that rare colonization events have been limited to single clones due to the extreme isolation of Mt. Erebus. Therefore, populations arising from recently introduced taxa would be expected to exhibit relatively low levels of genetic polymorphism and endemic microbial populations would be expected to show high levels of genetic polymorphism. In particular, we focused on accumulations of synonymous mutations, which were assumed to be under reduced selection pressure.

We developed the endemicity index (EI) (Fig. 3F) to assess the diversity of a microbial population represented by a MAG. Similar calculations have been successfully applied to approximate effective population size and genomic fluidity (64). High values (e.g. 10^− 2) of EI indicate high diversity which we interpreted as reflecting a neutral evolutionary process occurring under minimal contemporary selection pressures. Low values (e.g. 10^− 6) indicate low diversity, which we interpreted as possible evidence of a relatively recent arrival, a local population bottleneck, or a recent selective sweep.

A cyanobacterium belonging to either the Fischerella or Mastigocladus genus recovered the lowest median EI value (1 x 10^− 5) for a single species. This species dominated the near-surface. Although difficult to distinguish these two genera based on 16S rRNA gene sequence and GTDB classifies all members of the genera Fischerella and Mastigocladus as Fisherella, the distinction of Mastigocladus is recognized as a distinct genus by the List of Prokaryotic names with Standing in Nomenclature (LPSN). Therefore we classified this MAG as a member of the Mastigocladus genus to be consistent with the classification of specimens collected from Tramway Ridge in the past (65). The low endemicity index observed for this taxon was consistent with previous studies that showed that the global phylogeography of Mastigocladus reflects a geologically recent radiation from Yellowstone National Park, USA (66) and that the surface-associated microbial community at Tramway Ridge is likely dominated by aeolian-distributed cosmopolitan members of non-Antarctic temperate and terrestrial hydrothermal soil communities (2).

We identified Candidatus Australlarchaeum erebusii (ERB_5_1) and Candidatus Fervidibacteria erebusii (ERB_15_1) as two MAGs with relatively high EI values and abundance in the subsurface at depths greater than 2 cm (Fig. 3f, Fig. 3g). In our earlier amplicon-based study, these two taxa were similarly identified as dominant, endemic, and associated with the subsurface, but were referred to as “Thaumarchaeota-like archaeon” and OCtSpA1-106 respectively (2).

The most abundant endemic, subsurface-associated organism recovered was Candidatus A. erebusii (ERB_5_1), a member of the Nitrososphaeria (formerly Thaumarchaeota). This class of Archaea is a globally distributed group that is best known for the chemolithoautotrophic oxidation of ammonia (67) and for the apparent universal synthesis of cobalamin (68). However, several deeply divergent lineages of Nitrososphaeria identified through the construction of MAGs (45, 69, 70) and a single cultivated species, Candidatus Conexivisphaera calidus NAS-02 (71) have been shown to lack these hallmark attributes. These deeply diverging lineages have been predicted to be predominantly anaerobic heterotrophs (69, 70) capable of beta oxidation of fatty acids and protein/peptide degradation (71, 72). Candidatus A. erebusii (ERB_5_1) is one of the deepest-branching members of the Nitrososphaeria (Fig. 3b) and like other deeply diverging lineages, it encodes genes for the beta oxidation of fatty acids and peptide degradation while lacking marker genes for cobalamin biosynthesis and ammonia oxidation (Supplementary Table 5). Candidatus A. erebusii may also employ amino acid-based oxidation, similar to a pathway used by Thermococcus kodakarensis (73) and a proposed alternative metabolism for Ca. Nitrosocaldus islandicus (74), a representative of thermophilic ammonia oxidizing archaea (AOA). In this proposed metabolism, glutamate could be utilized to generate both reducing power and ATP to power the cell.

However, Ca. A. erebusii is also predicted to respire oxygen, a unique prediction among its closest, presumably anaerobic relatives. It encodes two aa3-type (low-affinity) cytochrome C-oxidases which could presumably drive beta oxidation of fatty acids. However it is unclear whether aerobic respiration could be coupled with an amino acid degradation pathway, which is typically thought to be an anaerobic metabolism (73, 74). It is unclear if these energy-generating pathways are mutually exclusive and operate under specific oxidative conditions. At least during summer, oxygen levels in the subsurface are around 30% saturation (2) and therefore, no organisms inhabiting the fumaroles are likely to be obligate anaerobes. However, it is also reasonable to assume that the wet, steamy subsurface experiences anoxia at least at small spatial scales. Therefore, we hypothesize that Ca. A. erebusii switches between metabolic pathways depending on the oxic environment, using beta-oxidation of fatty acids when oxygen levels are sufficiently high and peptide fermentation when oxygen levels are low. CO metabolism is likely to be used for maintenance during times of nutritional stress, as has previously been shown for Antarctic soil microbes (49).

Another abundant and endemic subsurface MAG examined in detail belongs to the thermophilic genus Candidatus Fervidibacteria, which was first discovered in Octopus Spring, Yellowstone National Park (clone OctSpA1-106, (75). The genus is named after Ca. Fervidibacter sacchari, which encodes a large repertoire of carbohydrate-active enzymes (CAZymes) (76) and which has recently been isolated and described (59). Like other Ca. Fervidibacter, Ca. F. erebusii (ERB_15_1) encodes a significant number and diversity of CAZymes including a large, diverse cohort of the unusual and poorly characterized GH109 family. Previous findings suggest that these novel enzymes are involved in extracellular polysaccharide metabolism unique to thermophilic systems (59, 77, 78) with the diversity of Ca. F. erebusii GH109s suggesting a diverse and unique polysaccharide utilization profile typical of this genus. Interestingly, Ca. Fervidibacter genomes, including Ca. F. erebusii appear to also encode mechanisms that enable growth on casamino acids (59). The encoded hydrogenase and CODH are likely maintenance mechanisms to enable survival under carbon limitation with growth primarily supported through aerobic heterotrophic growth on saccharides and anaerobic growth on amino acids.

In the current study, we have used metagenomics to assess several aspects of the microbial community inhabiting the fumarolic soils of Tramway Ridge, Mt. Erebus, Antarctica. We observed a shared functional repertoire between Tramway Ridge and other geothermal systems, specifically terrestrial hydrothermal systems such as hot springs. We then assessed metagenome-assembled genomes (MAGs) using a novel endemicity measure to identify two highly endemic taxa that are more abundant in the Tramway Ridge subsurface (> 2 cm depth) than the near-surface (< 2 cm depth). We named these two taxa Candidatus Australlarchaeum erebusii and Candidatus Fervidibacter erebusii to reflect that they were first observed as dominant, unique and endemic to Mt. Erebus. A close examination of the metabolic repertoire of these taxa revealed that they are likely both facultative anaerobic heterotrophs that specialize in using different carbon sources under aerobic conditions, but that use similar organic compounds during anaerobic growth. Like other deep-branching non-AOA Nitrososphaeria, Ca. A. erebusii possesses a putative pathway for the beta-oxidation of fatty acids. Like other Candidatus Fervidibacter, Ca. F. erebusii is predicted to utilize sugars and scavenge hydrogen gas under aerobic conditions (76, 79). Under oxygen-limited conditions, both may utilize similar peptides and amino acids for energy and carbon acquisition. We hypothesize that this pattern of metabolic utilization may reflect the extreme and carbon-limiting C:N ratios (1:3 to 3:1) encountered at Tramway Ridge. Dominant taxa may share nitrogen-rich compounds such as peptides and amino acids since the demand for such compounds may not be as high as for carbon-rich compounds. Instead, different taxa specialize in utilizing specific carbon compounds (sugars vs. fatty acids) through selective exclusion, allowing both to co-exist by carving out specific nutritional niches. Each is also equipped with form II aerobic carbon monoxide dehydrogenase genes that may provide maintenance energy during times of starvation. Together, these insights provide an unprecedented view into the dominant metabolic processes that may sustain life in this harsh, isolated environment.

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Availability of data and material

All raw metagenomic sequence data and metagenome assembled genomes used in this manuscript are available under bioproject accession PRJNA431961 (https://www.ncbi.nlm.nih.gov/bioproject/431961) . Data and scripts necessary to reproduce comparative metagenomic work is available at 10.5281/zenodo.10928974

Competing interests

Not applicable

Funding

Financial support was provided by grant UOW0802 from the New Zealand Marsden Fund to SCC and IRM and a CRE award from the National Geographic Society to S.C.C.

Author contributions

Planning and field sampling by CWH, CJV, SCC and IRM. Data analysis by CJV and CWH. All authors contributed to the interpretation of data and writing of the manuscript.

Acknowledgements

Antarctic logistic support for Event K-023 was provided by Antarctica New Zealand.

Parties of the Antarctic Treaty. Antarctic Treaty database [Internet]. ASPA 175 - High Altitude Geothermal sites of the Ross Sea region. 2014 [cited 2018 Oct 25]. https://www.ats.aq/devAS/info_measures_listitem.aspx?lang=e&id=573
Herbold CW, Lee CK, McDonald IR, Cary SC. Evidence of global-scale aeolian dispersal and endemism in isolated geothermal microbial communities of Antarctica. Nat Commun. 2014;5:3875.
Soo RM, Wood SA, Grzymski JJ, McDonald IR, Cary SC. Microbial biodiversity of thermophilic communities in hot mineral soils of Tramway Ridge, Mount Erebus, Antarctica. Environ Microbiol. 2009;11(3):715–28.
Broady PA. Taxonomic and ecological investigations of algae on steam-warmed soil on Mt Erebus, Ross Island, Antarctica. Phycologia. 1984;23(3):257–71.
Hudson JA, Daniel RM. Enumeration of thermophilic heterotrophs in geothermally heated soils from mount erebus, ross island, antarctica. Appl Environ Microbiol. 1988;54(2):622–4.
Kyle PR, McIntosh WC, Schmidt-Thomé M, Mueller P, Tessensohn F, Noll MR, et al. McMurdo Volcanic Group Western Ross embayment. In: LeMasurier WE, Thomson JW, Baker PE, Kyle PR, Rowley PD, Smellie JL, et al. editors. Volcanoes of the antarctic plate and southern oceans. Washington, D. C.: American Geophysical Union; 1990. pp. 18–145.
Harpel CJ, Kyle PR, Esser RP, McIntosh WC, Caldwell DA. 40Ar/39Ar dating of the eruptive history of Mount Erebus, Antarctica: summit flows, tephra, and caldera collapse. Bull Volcanol. 2004;66(8):687–702.
Esser RP, Kyle PR, McIntosh WC. 40Ar/39Ar dating of the eruptive history of Mount Erebus, Antarctica: volcano evolution. Bull Volcanol. 2004;66(8):671–86.
Barker PF, Thomas E. Origin, signature and palaeoclimatic influence of the Antarctic Circumpolar Current. Earth Sci Rev. 2004;66(1–2):143–62.
Ugolini FC, Starkey RL. Soils and Micro-organisms from Mount Erebus. Antarctica Nat. 1966;211(5047):440–1.
Vickers CJ, Herbold CW, Cary SC, Mcdonald IR. Insights into the metabolism of the high temperature microbial community of Tramway Ridge, Mount Erebus, Antarctica. Antarct Sci. 2016;28(4):241–9.
Wardell LJ, Kyle PR, Campbell AR. Carbon dioxide emissions from fumarolic ice towers, Mount Erebus volcano, Antarctica. Geological Society, London, Special Publications. 2003;213(1):231–46.
Ugolini FC. Soils of mount erebus, antarctica. NZ J Geol Geophys. 1967;10(2):431–42.
Vaquer-Sunyer R, Duarte CM. Thresholds of hypoxia for marine biodiversity. Proc Natl Acad Sci USA. 2008;105(40):15452–7.
Herbold CW, McDonald IR, Cary SC. Microbial ecology of geothermal habitats in antarctica. In: Cowan DA, editor. Antarctic Terrestrial Microbiology. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. pp. 181–215.
Sweeney D, Kyle PR, Oppenheimer C. Sulfur dioxide emissions and degassing behavior of Erebus volcano, Antarctica. J Volcanol Geoth Res. 2008;177(3):725–33.
Ilanko T, Fischer TP, Kyle P, Curtis A, Lee H, Sano Y. Modification of fumarolic gases by the ice-covered edifice of Erebus volcano, Antarctica. J Volcanol Geoth Res. 2019;381:119–39.
Mooshammer M, Wanek W, Hämmerle I, Fuchslueger L, Hofhansl F, Knoltsch A, et al. Adjustment of microbial nitrogen use efficiency to carbon:nitrogen imbalances regulates soil nitrogen cycling. Nat Commun. 2014;5:3694.
Barrett JE, Virginia RA, Wall DH, Cary SC, Adams BJ, Hacker AL, et al. Co-variation in soil biodiversity and biogeochemistry in northern and southern Victoria Land, Antarctica. Antarct Sci. 2006;18(4):535–48.
Nurk S, Bankevich A, Antipov D, Gurevich A, Korobeynikov A, Lapidus A, et al. Assembling Genomes and Mini-metagenomes from Highly Chimeric Reads. In: Deng M, Jiang R, Sun F, Zhang X, editors. Research in computational molecular biology. Berlin, Heidelberg: Springer Berlin Heidelberg; 2013. pp. 158–70.
Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34.
Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31(10):1674–6.
Wu Y-W, Simmons BA, Singer SW. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics. 2016;32(4):605–7.
Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;3:e1165.
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 2017;11(12):2864–8.
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2019;36(6):1925–7.
Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36(10):996–1004.
Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, et al. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016;44(14):6614–24.
Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, von Mering C, et al. Fast Genome-Wide Functional Annotation through Orthology Assignment by eggNOG-Mapper. Mol Biol Evol. 2017;34(8):2115–22.
Sousa FL, Alves RJ, Pereira-Leal JB, Teixeira M, Pereira MM. A bioinformatics classifier and database for heme-copper oxygen reductases. PLoS ONE. 2011;6(4):e19117.
Søndergaard D, Pedersen CNS, Greening C. HydDB: A web tool for hydrogenase classification and analysis. Sci Rep. 2016;6:34212.
Zheng J, Ge Q, Yan Y, Zhang X, Huang L, Yin Y. dbCAN3: automated carbohydrate-active enzyme and substrate annotation. Nucleic Acids Res. 2023;51(W1):W115–21.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.
Hyatt D, Chen G-L, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.
Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, et al. The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 2016;44(D1):D279–85.
Chen I-MA, Chu K, Palaniappan K, Ratner A, Huang J, Huntemann M, et al. The IMG/M data management and analysis system v.7: content updates and new features. Nucleic Acids Res. 2023;51(D1):D723–32.
Oksanen J, Simpson GL, Blanchet FG, Kindt R, Legendre P, Minchin PR et al. vegan: Community Ecology Package. 2022.
Paradis E, Schliep K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics. 2019;35(3):526–8.
Murdoch D, Adler D. rgl: 3D Visualization Using OpenGL. 2023.
Krijthe JH, Van der Maaten L, Rtsne. T-distributed stochastic neighbor embedding using Barnes-Hut implementation. R package version 013, URL https://githubcom/jkrijthe/Rtsne. 2015.
Wickham H. ggplot2: Elegant Graphics for Data Analysis (Use R!). 2nd ed. Cham: Springer; 2016.
Lai D, Hedlund BP, Mau RL, Jiao J-Y, Li J, Hayer M, et al. Resource partitioning and amino acid assimilation in a terrestrial geothermal spring. ISME J. 2023;17(11):2112–22.
Beam JP, Jay ZJ, Kozubal MA, Inskeep WP. Niche specialization of novel Thaumarchaeota to oxic and hypoxic acidic geothermal springs of Yellowstone National Park. ISME J. 2014;8(4):938–51.
Kato S, Sakai S, Hirai M, Tasumi E, Nishizawa M, Suzuki K, et al. Long-Term Cultivation and Metagenomics Reveal Ecophysiology of Previously Uncultivated Thermophiles Involved in Biogeochemical Nitrogen Cycle. Microbes Environ. 2018;33(1):107–10.
Spieck E, Spohn M, Wendt K, Bock E, Shively J, Frank J, et al. Extremophilic nitrite-oxidizing Chloroflexi from Yellowstone hot springs. ISME J. 2020;14(2):364–79.
Mueller AJ, Daebeler A, Herbold CW, Kirkegaard RH, Daims H. Cultivation and genomic characterization of novel and ubiquitous marine nitrite-oxidizing bacteria from the Nitrospirales. ISME J. 2023;17(11):2123–33.
Ji M, Greening C, Vanwonterghem I, Carere CR, Bay SK, Steen JA, et al. Atmospheric trace gases support primary production in Antarctic desert surface soil. Nature. 2017;552(7685):400–3.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55.
Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4.
Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9.
Minh BQ, Nguyen MAT, von Haeseler A. Ultrafast approximation for phylogenetic bootstrap. Mol Biol Evol. 2013;30(5):1188–95.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42(Database issue):D222–30.
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35(8):725–31.
Brown CT, Hug LA, Thomas BC, Sharon I, Castelle CJ, Singh A, et al. Unusual biology across a group comprising more than 15% of domain Bacteria. Nature. 2015;523(7559):208–11.
Nelson WC, Stegen JC. The reduced genomes of Parcubacteria (OD1) contain signatures of a symbiotic lifestyle. Front Microbiol. 2015;6:713.
King GM, Weber CF. Distribution, diversity and ecology of aerobic CO-oxidizing bacteria. Nat Rev Microbiol. 2007;5(2):107–18.
Hedlund B, Nou N, Covington J, Lai D, Mayali X, Seymour C et al. Genome-guided isolation of Fervidibacter sacchari, an aerobic, hyperthermophilic polysaccharide-degrading specialist. Res Sq. 2024.
Greening C, Berney M, Hards K, Cook GM, Conrad R. A soil actinobacterium scavenges atmospheric H2 using two membrane-associated, oxygen-dependent [NiFe] hydrogenases. Proc Natl Acad Sci USA. 2014;111(11):4257–61.
Koch H, Galushko A, Albertsen M, Schintlmeister A, Gruber-Dorninger C, Lücker S, et al. Growth of nitrite-oxidizing bacteria by aerobic hydrogen oxidation. Science. 2014;345(6200):1052–4.
Archer SDJ, Lee KC, Caruso T, Maki T, Lee CK, Cary SC, et al. Airborne microbial transport limitation to isolated Antarctic soil habitats. Nat Microbiol. 2019;4(6):925–32.
Power JF, Carere CR, Welford HE, Hudson DT, Lee KC, Moreau JW, et al. A genus in the bacterial phylum Aquificota appears to be endemic to Aotearoa-New Zealand. Nat Commun. 2024;15(1):179.
Andreani NA, Hesse E, Vos M. Prokaryote genome fluidity is dependent on effective population size. ISME J. 2017;11(7):1719–21.
Melick DR, Broady PA, Rowan KingsleyS. Morphological and physiological characteristics of a non-heterocystous strain of the cyanobacterium Mastigocladus laminosus Cohn from fumarolic soil on Mt Erebus, Antarctica. Polar Biol. 1991;11(2).
Miller SR, Castenholz RW, Pedersen D. Phylogeography of the thermophilic cyanobacterium Mastigocladus laminosus. Appl Environ Microbiol. 2007;73(15):4751–9.
Pester M, Schleper C, Wagner M. The Thaumarchaeota: an emerging view of their phylogeny and ecophysiology. Curr Opin Microbiol. 2011;14(3):300–6.
Doxey AC, Kurtz DA, Lynch MDJ, Sauder LA, Neufeld JD. Aquatic metagenomes implicate Thaumarchaeota in global cobalamin production. ISME J. 2015;9(2):461–71.
Ren M, Feng X, Huang Y, Wang H, Hu Z, Clingenpeel S et al. Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution. ISME J. 2019.
Hua Z-S, Qu Y-N, Zhu Q, Zhou E-M, Qi Y-L, Yin Y-R, et al. Genomic inference of the metabolism and evolution of the archaeal phylum Aigarchaeota. Nat Commun. 2018;9(1):2832.
Kato S, Itoh T, Yuki M, Nagamori M, Ohnishi M, Uematsu K, et al. Isolation and characterization of a thermophilic sulfur- and iron-reducing thaumarchaeote from a terrestrial acidic hot spring. ISME J. 2019;13(10):2465–74.
Lin X, Handley KM, Gilbert JA, Kostka JE. Metabolic potential of fatty acid oxidation and anaerobic respiration by abundant members of Thaumarchaeota and Thermoplasmata in deep anoxic peat. ISME J. 2015;9(12):2740–4.
Yokooji Y, Sato T, Fujiwara S, Imanaka T, Atomi H. Genetic examination of initial amino acid oxidation and glutamate catabolism in the hyperthermophilic archaeon Thermococcus kodakarensis. J Bacteriol. 2013;195(9):1940–8.
Daebeler A, Herbold CW, Vierheilig J, Sedlacek CJ, Pjevac P, Albertsen M, et al. Cultivation and Genomic Analysis of Candidatus Nitrosocaldus islandicus, an Obligately Thermophilic, Ammonia-Oxidizing Thaumarchaeon from a Hot Spring Biofilm in Graendalur Valley, Iceland. Front Microbiol. 2018;9:193.
Blank CE, Cady SL, Pace NR. Microbial composition of near-boiling silica-depositing thermal springs throughout Yellowstone National Park. Appl Environ Microbiol. 2002;68(10):5123–35.
Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng J-F, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013;499(7459):431–7.
Liu QP, Sulzenbacher G, Yuan H, Bennett EP, Pietz G, Saunders K, et al. Bacterial glycosidases for the production of universal red blood cells. Nat Biotechnol. 2007;25(4):454–64.
Teze D, Shuoker B, Chaberski EK, Kunstmann S, Fredslund F, Nielsen TS et al. The Catalytic Acid–Base in GH109 Resides in a Conserved GGHGG Loop and Allows for Comparable α-Retaining and β-Inverting Activity in an N -Acetylgalactosaminidase from Akkermansia muciniphila. ACS Catal. 2020;3809–19.
Hedlund BP, Dodsworth JA, Murugapiran SK, Rinke C, Woyke T. Impact of single-cell genomics and metagenomics on the emerging view of extremophile microbial dark matter. Extremophiles. 2014;18(5):865–75.

No competing interests reported.

Download PDF

Editorial decision: Revision requested
09 Sep, 2024
Reviews received at journal
02 Sep, 2024
Reviews received at journal
29 Aug, 2024
Reviewers agreed at journal
29 Aug, 2024
Reviewers agreed at journal
13 Aug, 2024
Reviewers agreed at journal
07 Aug, 2024
Reviewers invited by journal
07 Aug, 2024
Editor assigned by journal
01 Aug, 2024
Submission checks completed at journal
27 Jul, 2024
First submitted to journal
25 Jul, 2024

You are reading this latest preprint version

Nutritional niches of endemic, facultatively anaerobic heterotrophs from an isolated Antarctic terrestrial hydrothermal refugium elucidated through metagenomics

Status:

Version 1

Abstract

Background

Results

Conclusions

Figures

Background

Methods

Sample collection

DNA extraction, library preparation and sequencing

Metagenome assembly and binning

Endemicity Index Calculation

Functional profile comparisons

Phylogenomic analysis

Results

Discussion

Metagenome Functional Profiles

Conclusion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1