Skip to main content

Comparison of prokaryotes between Mount Everest and the Mariana Trench



Mount Everest and the Mariana Trench represent the highest and deepest places on Earth, respectively. They are geographically separated, with distinct extreme environmental parameters that provide unique habitats for prokaryotes. Comparison of prokaryotes between Mount Everest and the Mariana Trench will provide a unique perspective to understanding the composition and distribution of environmental microbiomes on Earth.


Here, we compared prokaryotic communities between Mount Everest and the Mariana Trench based on shotgun metagenomic analysis. Analyzing 25 metagenomes and 1176 metagenome-assembled genomes showed distinct taxonomic compositions between Mount Everest and the Mariana Trench, with little taxa overlap, and significant differences in genome size, GC content, and predicted optimal growth temperature. However, community metabolic capabilities exhibited striking commonality, with > 90% of metabolic modules overlapping among samples of Mount Everest and the Mariana Trench, with the only exception for CO2 fixations (photoautotrophy in Mount Everest but chemoautotrophy in the Mariana Trench). Most metabolic pathways were common but performed by distinct taxa in the two extreme habitats, even including some specialized metabolic pathways, such as the versatile degradation of various refractory organic matters, heavy metal metabolism (e.g., As and Se), stress resistance, and antioxidation. The metabolic commonality indicated the overall consistent roles of prokaryotes in elemental cycling and common adaptation strategies to overcome the distinct stress conditions despite the intuitively huge differences in Mount Everest and the Mariana Trench.


Our results, the first comparison between prokaryotes in the highest and the deepest habitats on Earth, may highlight the principles of prokaryotic diversity: although taxa are habitat-specific, primary metabolic functions could be always conserved.

Video abstract.


Extreme environments are widespread on Earth, such as hot springs, hydrothermal vents, deserts, glaciers, deep sea, and hadal trenches, which have a wide variety of harsh conditions and provide unique and diverse ecological niches for life on Earth [1]. Prokaryotes are the primary life forms in extreme environments, teaching us about life’s limits [2,3,4]. Prokaryotes are the most widespread form of life, encompassing incredibly diverse taxa and encoding the vast majority of biological functions, playing crucial roles in changing, shaping, and sustaining habitability on Earth [5,6,7]. The extraordinary diversity of prokaryotes is believed to be a product of ~ 3.8 billion years of evolution, 2 billion years longer than that of eukaryotic organisms [8]. It is generally accepted that prokaryotic species are closely related to the environmental conditions of their habitats [9] and most species-level genes (> 95% nucleotide identity) are specific to a single habitat [6]. However, the functional redundancy that multiple prokaryotic taxa encode the same biological function [10] and cross-stress adaptation behaviors for coping with distinct environmental conditions [3] are also reported in many studies. These two seemingly contradictory results bring us a new hypothesis: the prokaryotic taxa tend to adapt to the environment; in contrast, most metabolic pathways (including some adaptation mechanisms) are still conserved, even in two completely different extreme environments. To test this hypothesis, here we chose the highest habitats and the deepest habitats for prokaryotes on Earth, Mount Everest (ME) and the Mariana Trench (MT), respectively, for comparing the taxonomy and metabolic capability of the prokaryotic community and gaining new insights into the prokaryotic biodiversity in two isolated extreme environments.

Mount Everest (ME) is the highest area on Earth, with an 8.8-km elevation summit, and is covered by snow and ice [11]. Observations in situ reveal that > 6 km above sea level (a.s.l) on ME are environments with intense ultraviolet (UV) radiation, as well as the low but dramatically changeable temperatures that vary with seasons and day/night cycles [12]. The automatic meteorological station at 6.5 km a.s.l. recorded the highest net radiation, 392.52 W/m2, and temperature ranged from −19.3 to 9.4 °C, with an average of −5.5 °C from May to July 2021 (unpublished data). In addition, the high elevation of ME also causes lower pressure than the standard atmospheric pressure (0.0033 MPa at the summit of Mount Everest) [13] and lower nutrient concentration than normal environments (e.g., in soils) [11, 14, 15]. These extreme environmental parameters challenge all life forms living in ME. Bacterial abundances in surface snow, moraine lakes, glacier streams, and meltwaters of ME above 6 km were reported with a low concentration of cells, approximately 104–105 cells/mL [11, 16,17,18]. Seasonal variation has been observed for bacterial community structure in snow samples [19]. Both culture and 16S rRNA gene amplicon-based studies reveal that the bacterial community in ME is relatively simple, commonly predominated by Proteobacteria, Actinobacteria, Firmicutes, Cyanobacteria, and Deinococcus-Thermus [11, 14, 16, 17, 20]. Thus far, the recent published study on genome and gene catalog of glacier microbiomes is the only one publication that included the prokaryotic metagenome at 4.6–6.7 km on the Tibetan Plateau, while the metabolic potentials and adaptative mechanisms of prokaryotes in ME remain unclear which needs further investigation [4].

In contrast, the Mariana Trench (MT) is the deepest region on Earth. The 11-km water depth below sea level (b.s.l) of the Challenger Deep, the deepest area in MT, is far beyond areas reached by sunlight. For this reason, the MT is characterized by a stable and near-freezing temperature (~ 2°C), total darkness, and poor nutrient availability, but under an ultrahigh hydrostatic pressure (HHP) up to ~ 115 MPa [21, 22]. Cellular life living in MT usually needs unique mechanisms to cope with the HHP, such as intercellular accumulation of compatible solutes (e.g., TMAO) in both macro- and microorganisms [23,24,25]. The biomass density at the bottom of MT is significantly lower than that at the marine surface, 103–104 cells/mL in near-bottom water and 105–107 cells/mL in sediments [26, 27]. Microbiomes in MT samples have been reported to be with high novelties, especially in bottom-axis sediment samples [28]. In addition, the V-shaped topography of MT enables the collection of particulate organic matters (POMs) and heavy metals, e.g., arsenic (As) and selenium (Se), from the overlying water column, abyssal seafloor, and Earth’s upper crust [28, 29]. These accumulation effects result in a higher content of organic matters (OMs) and more intense microbial degradation activities at the bottom-axis than in adjacent slope sites in MT [30], as well as a higher concentration (> 2-fold) of As and Se in MT sediment samples than those at non-hadal sites [28]. The only two previous metagenomic studies of the sediments in the Challenger Deep of MT suggested the essential roles of predominant heterotrophs in recycling macromolecules and utilizing various carbon sources (e.g., peptides, carbohydrates, hydrocarbons, and aromatic compounds) [21, 28].

In this study, we compare the taxonomy and metabolic capability of prokaryotes between ME and MT by using nine samples of ice/snow at 6–7 km a.s.l. from the East Rongbuk Glacier of ME and 16 samples of water/sediment with depths of 7–11 km b.s.l from the Challenger Deep in MT (Fig. 1). Metagenomic analysis between ME and MT revealed a surprisingly opposite scenario in terms of taxonomy and metabolic capability. The taxonomic diversity between ME and MT was highly variable, with significant differences in community composition, genome size, GC content, and predicted optimal growth temperature between these two habitats. In contrast, the diversity of metabolic capabilities was much more consistent between ME and MT. Most metabolic modules (> 90%) were shared across ME and MT. Near-complete metagenome-assembled genomes (MAG)-based analysis suggested the functional redundancy that distinct prokaryotic taxa can carry the same metabolic functions in ME and MT. These results provide extreme examples for prokaryotic comparison from the habitats with the greatest height differences on Earth, providing a new perspective to understanding prokaryotic biodiversity in terms of taxonomy and metabolism and their specificity to different extreme habitats.

Fig. 1
figure 1

Sampling sites in Mount Everest (ME) and the Mariana Trench (MT). A Location of the ME and MT on the map of Google Earth. Detailed sampling sites in the ME (B) and the MT (C). D Schematic diagram of the samples used in this study


Sampling and preprocessing

Since the environmental conditions of ME and MT are so different, it is almost impossible to collect exactly the same type of sample to compare. Here we selected the most representative samples of ME or MT, including 9 samples (snow and ice) from ME and 16 samples (water and sediment) from MT (Fig. 1, Table 1). Among them, snow and ice represent the typical habitats that harbored diverse prokaryotes on glacier surface of Mount Everest [31]; similarly, water and sediment are the main habitats of prokaryotes in trenches [32]. For ME samples, 1 snow sample and 8 ice samples were collected at 8 sites on northern slope, the East Rongbuk Glacier, of Mount Everest with different elevations ranging from 6280 to 6675 m. The snow sample was collected at ~5-cm depth by a sterile steel spoon, and eight ice samples were collected with pre-clean ice axe, then surface 2 mm of the ice samples was scraped with a sterile scalpel for decontamination. Detailed information of sample sites is shown in Table 1. About 2L of each sample was placed in WhirlPak bags (WhirlPak®, Nasco, USA), and all samples were placed in icebox during transport from sampling site to the laboratory. Snow and ice samples were melted at 10 °C in the lab and filtered onto a 0.22-μm polycarbonate filter membrane (47-mm diameter; Millipore) for further DNA extraction.

Table 1 Sampling sites in this study

For MT samples, 10 water samples and six sediment samples were collected at 13 sites near the Challenger Deep in the Mariana Trench aboard the R/V “Tansuo-01,” with varying water depth ranging from 7344 m to ~ 11 km b.s.l, under hydrostatic pressure from 75 to 115 MPa. Among these 16 samples, 13 were collected from the bottom of the Challenger Deep with nearly the deepest water depth (> 10 km) on Earth. The temperature of these 16 MT sampling sites ranged from 1.8 to 2.6 °C and varies with the water depth (Table 1). All seawater samples were collected using Niskin bottles from the Mariana Trench during TS-09 cruise in September and October 2018. Samples from station TY044 and TY041 were filtered through 0.22-μm polycarbonate membranes (142 mm, Millipore) directly, and other 8 samples were filtered serially through 3-μm and 0.22-μm polycarbonate membranes aboard. After adding RNAlaterTM Stabilization Solution (Thermo Scientific, Wilmington, DE, USA), all filters were frozen at −80°C immediately until further analysis. In situ hydrographical parameters (e.g., salinity, temperature, and depth) were measured by a conductivity-temperature-depth (CTD) rosette system (Sea-Bird Electronics). Sediment samples were collected by box-core at stations TY038 and TY042 during TS-09 cruise, and by gravity-core at station GT04 during the cruise TS01-14. The box-core samples were mixed and frozen at −80 °C immediately until further analysis, while the gravity-core samples were cut at each 2–4 cm and then frozen at −80°C until further analysis. Detailed information of sampling is shown in Table 1.

DNA extraction and sequencing

Metagenomic DNA of snow, ice, and seawater samples was extracted from 0.22-μm filters with a PowerSoil DNA Isolation Kit (Qiagen, Germantown, USA) following its protocol. The concentration of DNA was quantified by Qubit DNA Assay Kit with a Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA), and the quality was checked by gel electrophoresis. The sequencing library was prepared using the KAPA hyper Prep Kit (Roche, Shanghai, China). The shotgun sequencing was performed on Illumina NovaSeq 6000 Platform PE150 (Illumina, San Diego, CA, USA) with 150 bp paired-end reads. Genomic DNA of sediment samples from TY038 and TY042 stations was extracted from ~ 5 g by modified SDS method as descript in a previous study [33]. DNA of sediment samples from GT04 station was extracted from ~ 0.3 g of sediment per sample using PowerSoil® DNA Isolation Kit (QIAGEN, Germany). The sequence library was built in Beijing Genome Institution (BGI, Shenzhen, China) and sequenced on the BGI MGISEQ-2000 platform, which has been reported can obtain comparable results as Illumina platform [34].

Quality control and reads analysis

Metagenomic raw reads are adapter-trimmed using TrimGalore (version 0.6.6) and filtered using sickle (version 1.33) with parameters (--length-threshold 90 -t sanger -g) [35]. Filtered read quality is estimated and reported by fastQC (v0.11.4) and multiQC (version 1.9) [36]. Small-subunit rRNA (SSU rRNA) genes are assembled from clean reads using phyloFlash (v3.4) with parameters (-almosteverything), depending on SILVA database version 138.1 [37, 38]. To assess the alpha diversity, Shannon and richness indexes were evaluated based on counts of SSU sequences extracted from all clean reads of the metagenomes of ME and MT samples, from domain to species levels. Metagenomic contigs are assembled from clean reads using megahit (v1.2.9) with default parameters and --min-contig-len 500 [39]. Clean reads were then mapped to contigs using (last modified on February 13, 2020) with parameters (nodisk k=13 minid=0.95 keepnames=t minaveragequality=5 trimreaddescriptions=t pairlen=350 rescuedist=650) [40]. Genes are predicted by prodigal (V2.6.3) with parameter (-p meta) [41, 42]. Genes not shorter than 99 bp (33 aa) are then annotated by KofamKOALA (version 1.3.0, database build 2021-01-04) and eggNOG (version 2.1.3, database 5.0.2) to assign KEGG KO to genes [43]. We also download a description of a pathway to module (ko00002), pathways, and module for further analysis. After gene annotation, gene coverage is calculated using featureCounts (v2.0.1) [44].

Binning and MAG-level analysis

Metagenomic binning and refinement were performed on contigs using a combination of metabat2 (version 2:2.15), maxbin2 (version 2.2.7), concoct (version 1.1.0), and DAS Tool (version 1.1.2) (Additional file 1: Fig. S1) [45,46,47]. Megahit was run with nine combinations of parameters (--maxP with 60, 75, 90, and –minS with 60, 75, 90). Metabat2 was used with two combinations of parameters (-markerset with 107 and 40). Binnig results of the total 12 parallel methods are then refined using DASTool with parameters (--search_engine diamond –score_threshold 0). Next, binning results from all samples are collected and representative MAGs are choired using dRep (v3.0.0) with parameters (-comp 50 -con 10 -pa 0.9 -sa 0.95) [48]. Finally, 1176 MAGs with completeness ≥ 50 and contaminant ≤ 10 including 411 MAGs with completeness ≥ 90 and contaminant ≤ 5 are clustered into 648 clusters [49]. We used GTDB-Tk (v1.6.0) to taxonomically classify representative MAGs of all clusters with the GTDB release 202 [50, 51]. Next, we use prodigal without parameter (-p meta) and annotate gene function as descript above. Reads are mapped to genes of all MAG genomes by minimap (version 2.1) equivalent to the method used by coverM (version 0.6.1) [52].

Statistical analysis

mOTUs are first filtered and only keep prokaryotes, and analyzed at phylum to species level, any taxon that exists in fewer than two samples or maps lower than ten reads are filtered. Mapped reads of genes are grouped by KO or KEGG module for each sample [6]. Distances between samples were calculated by vegdist in R package vegan, when calculating the binary Jaccard distance, we use binary format of reads table, and the Bray-Curtis distance between samples is calculated using total reads [53]. A method equivalent to TPM was used to normalize gene abundance data for metagenomic comparation. Briefly, reads mapped to genes were normalized by gene length and total mapped number following the formula below:

$$\textrm{TP}{\textrm{M}}_{\textrm{k}}=\frac{{\textrm{r}}_{\textrm{k}}/{\textrm{l}}_{\textrm{k}}}{\sum \left({\textrm{r}}_{\textrm{i}}/{\textrm{l}}_{\textrm{i}}\right)}\times {10}^6$$

in which ri and li refer to the reads mapped to gene i and the length of gene i. The significant difference between different locations was calculated using adonis method in R package Vegan. Significance of the MAG-based comparisons of predicted G+C content, genome size, optimal growth temperature, and minimum generation times between different locations are using Wilcox test in ggsignif package in R [54]. Co-occurrence networks of genes were constructed based on their presence/absence across ME or MT-specific MAGs based on binary Jaccard distance. Edges in the network represent a co-occurrence level of > 0.3. A full record of all statistical analysis is included as Additional file 4 and was created using the knitr package in R [55].

Phylogenomic analysis of MAGs

A phylogenetic tree of reconstructed representative genomes above middle-quality was constructed using the protein sequences of 40 universal gene markers. Approximate maximum-likelihood trees were generated using FastTree (version 2.1.10) and bootstrap values of both phylogenomic trees were evaluated based on 1000 replicates. Other 100 reference genomes from GTDB database were chosen according to fastANI similarity. Genes of each genome are predicted by prodigal in “single” mode, and fetchMGs (v1.2) was used to extract 40 conserved single-copy genes of prokaryotes [56]. Genomes with less than 20 marker genes extracted were filtered. Marker genes are aligned by MAFFT (v7.487) with parameters (--maxiterate 1000 --localpair) and concatenated to trim by trimAl (v1.4.rev15) with the automated1 mode. FastTree (version 2.1.10) was used to construct the tree with parameter (-gamma). All trees were further polished and visualized using the interactive tree of life (iTOL) v5 [57].

MAG-based physiological predictions

Estimation of optimal growth temperature was performed for medium-high-quality MAGs using Growthpred by searching highly expressed gene (ribosomal protein gene) in genomes [58]. We discarded MAGs that belong to the Archaea domain and the number of highly expressed gene blow 10, then estimate minimum generation times using with parameters (-r -t -c 0) [59].


Description of prokaryotes from ME and MT

Total of the 9 ME and 16 MT samples yielded 4,846,075,106 clean raw read pairs (814,407,384 for ME, 4,031,667,722 for MT) after being trimmed and filtered, composing 98% of the raw reads (see the “Methods” section). Since eukaryotic reads account for only 3.45% in average of 25 metagenomes from ME and MT, we only focus on the prokaryotes in this study. Assembly of the above 25 metagenomes generated 16,958,325 scaffolds (average length 1238.34 bp). Binning of these contigs resulted in 1176 draft genomes with qualities above the middle-quality threshold (≥ 50% completeness and ≤ 10% contamination), among which 411 were above the high-quality threshold (≥ 90% completeness and ≤ 5% contamination) according to a widely accepted standard of MAGs [49]. Quality-controlled reads from all metagenomes had relatively high mapping rates (51.5%) to all MAGs above the middle-quality threshold, suggesting an acceptable representation and quality for the use of these MAGs to reflect the prokaryotic community in further analyses.

Distinct taxonomic composition of prokaryotic communities between ME and MT

A total of 1,841,892 reads were mapped to 4084 molecular operational taxonomic unit (mOTU) sequences at the species level or higher based on the MiTAG method using phyloFlash (v3.4) (Methods, Additional file 1: Fig. S1A, Additional file 2: Table S1). A high proportion of microbial dark matter was observed in samples from these two extreme environments. For example, ~50% of family-level and ~80% of genus-level mOTUs were unclassified according to the SILVA 138.1 database [38]. A similar conclusion could be obtained from MAGs, in which 34.5% of high-quality MAGs were unclassified against the Genome Taxonomy database (GTDB) database at the species level.

We independently calculated the alpha and beta diversity of prokaryotic metagenomes from ME and MT samples to reflect and compare the prokaryotic diversity and community dissimilarity (see the “Methods” section). The richness index of the prokaryotic community of ME was significantly lower than that in MT samples at the phylum level, while the opposite pattern was observed at the species level (Welch’s t test; p value < 0.05) (Additional file 1: Fig. S2). Beta diversity was assessed to compare the dissimilarity of prokaryotic communities between ME and MT. This dissimilarity between the two extreme environments increased from the domain to species level (Fig. 2C) and significantly differed at levels below class (Adonis test; p value < 0.05), regardless of the presence (binary Jaccard distance) or abundance (Bary-Cruit distance) of taxa (Fig. 2A, Additional file 1: Fig. S3). These comparisons of alpha and beta diversity reveal apparent differences in the taxonomic composition of prokaryotic communities between ME and MT. The relative abundance of taxa to class level among ME and MT samples calculated by mOTUs also shows the distinct community composition between ME and MT (Additional file 1: Fig. S1B).

Fig. 2
figure 2

Taxonomic and functional diversity of microbiome in Mount Everest (ME) and the Mariana Trench (MT). A Taxonomic composition between ME and MT. Calculation was performed on species abundancy against the clean reads of metagenome. B Functional diversity between ME and MT. Calculation was performed on KO abundancy against the clean reads of metagenome. C Comparison of dissimilarities of function and taxonomy within ME, within MT and across two habitats. The dissimilarities of functions and taxonomy were reflected by the Bray-Curtis distance. Statistic significances were calculated by Wilcox test (*p-value < 0.05; **p-value < 0.001; ***p-value < 1e−5), while no significance was represented by “NS”. D Sankey network of taxonomy and function between ME and MT. Calculation was performed on 648 representative MAGs (above middle-quality, after removing redundancy). We used class level to reflect taxonomy and KEGG metabolic modules to reflect function. ME-specific taxa and functions were represented in red, and MT-specific were represented in blue, while common taxa and functions were represented in purple. More detailed information is shown in Additional file 1: Fig. S1 and Fig. S6

High environmental specificity in prokaryotic taxa

In addition to the analysis of clean reads, the MAG-based Sankey network also revealed the distinct prokaryotic taxa between ME and MT (see the “Methods” section; Fig. 2D). We further constructed a phylogenetic tree by using 648 representative MAGs (dereplicated at 95% average nucleotide identity (ANI)) and 100 reference genomes from the GTDB [50] (see the “Methods” section, Fig. 3A). Except for two unclassified MAGs, the remaining MAGs were distributed into 47 phyla. Among them, seven phyla were commonly identified at both ME and MT. Regardless of ME or MT, most MAGs belonged to the phylum Proteobacteria, followed by the phyla Chloroflexota, Planctomycetota, Bacteroidota, Actinobacteriota, Marinisomatota, and Gemmatimonadota. Specific phyla from ME samples were Cyanobacteria and Deinococcota. No archaea group was observed in metagenomes from ME, which matched the previous observation based on 16S rRNA gene [11]. Conversely, 38 specific phyla were observed only in MT samples, including four phyla belonging to the archaea domain, i.e., Nanoarchaeota, Thermoplasmatota, Iainarchaeota, and Thermoproteota (formerly known as Thaumarchaeota), and 34 bacterial phyla, such as Chloroflexota, Marinisomatota, Nitrospirota, Planctomycetota, Elusimicrobiota, Hydrogenedentota, and Omnitrophota.

Fig. 3
figure 3

Phylogeny, taxonomy, and comparison genomic and predicted physiological features of MAGs between ME and MT. A Phylogenetic tree of MAGs from ME and MT in the species level. This phylogeny was generated using 648 genomes after dereplication clustered into species level (ANI ≥ 95). Bar plots are shown completeness (red) and contamination (brown) of MAGs. Colorstrip in outer ring represents the occurrence of MAGs between ME and MT in the genus level, such as red for ME, blue for MT, and purpose for present in both (shared genus labeled in outer line). Bootstraps are shown in the purple triangle and are based on 1000 replicated trees. Comparison of B genome size, C GC content, and D predicted optimal growth temperature of MAGs between ME and MT. Red dots and boxes represent the data from ME, while blue dots and boxes represent the data from MT. Significant difference between ME and MT is represented by “*” and “***” where the p-value is less than 0.05 and 1e−5 (Wilcox test), representatively

Additionally, ME and MT shared a very small number of taxa, especially at a finer level (at or below the family) (Additional file 1: Fig. S6). For ~92% MAGs above middle-quality, the ANIs between genomes from ME and MT were all less than 80%, lower than the commonly accepted within-family cutoff of genome similarity [60, 61] (Additional file 1: Fig. S4). Furthermore, only eight genera were shared across ME and MT, including the genera Acinetobacter (15 MAGs), Sphingomonas (9 MAGs), Brevundimonas (5 MAGs), Comamonas (6 MAGs), Phenylobacterium (3 MAGs), Pedobacter (3 MAGs), and Nocardioides (9 MAGs) (Fig. 3A, Additional file 1: Fig. S6). At the species level, 29 and 99 species were unique to ME and MT, respectively (Additional file 1: Fig. S6). The only shared species was annotated as Comamonas tsuruhatensis, and the ANI was 96% among C. tsuruhatensis MAGs in ME and MT (Additional file 1: Fig. S4). The above results fully demonstrated that the taxonomic diversity of prokaryotes varied considerably between ME and MT and that most taxa at or below the family level are specific to either ME or MT.

Distinct prokaryotic genomic and physiological characteristics between ME and MT

We further compared the genomic and physiological characteristics of prokaryotes between ME and MT based on all MAGs above middle quality (Fig. 3B–D, Additional file 1: Fig. S8A-C). The medium size of prokaryotic genomes from ME was ~ 4 Mbp, significantly larger than those from MT at ~ 3 Mbp (Wilcox test; p value < 1e−5) (Fig. 3B). The medium GC content of genomes in ME was ~ 57%, which is slightly higher than that in MT of ~ 56% (Wilcox test; p value < 0.05) (Fig. 3C). We also predicted the optimal growth temperature (OGT) according to MAGs above middle-quality as an indicator to represent the physiological features of prokaryotes at ME and MT (see the “Methods” section). These predictions showed that prokaryotes from ME had a significantly lower OGT than those from MT (Wilcox test; p value < 0.05), of which medium OGTs at 25 °C and 40 °C, respectively (Fig. 3D). The above results suggested significant differences in the genomic and physiological features of prokaryotes between these two distinct extreme environments.

Similar metabolic capabilities between the communities in ME and MT

In contrast to the significant differences in taxonomic characteristics between ME and MT, the metabolic capabilities of prokaryotic communities were more conserved. The two groups for ME and MT were clearly separated in the nonmetric multidimensional scaling (NMDS) ordination of taxonomy below the family level (Fig. 2A, Additional file 1: Fig. S3). However, the NMDS ordination of KEGG orthologs (KO) overlapped between ME and MT (Fig. 2B). The dissimilarity of KO composition was significantly lower than the dissimilarity of taxonomic composition (Wilcox test, p value < 1e−5) (Fig. 2C). For metabolic modules, the dissimilarity between the two habitats was even as low as the cross-sample variation in ME (Fig. 2C). Similar trends were also revealed by MAG-based analysis (Additional file 1: Fig. S6, Fig. S7). The overlapped KOs between ME and MT occupied 73% of all annotated KOs from MAGs above middle quality, whereas the ME-specific and MT-specific KOs were only 10% and 17%, respectively. After mapping KOs to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway modules, 91% of mapped function modules were conserved, while only 1 and 8% were specific to ME and MT, respectively. These results imply that almost all metabolic capabilities are shared between ME and MT.

Metabolic capabilities of prokaryotes in ME and MT

The MAG-based Sankey network of taxonomy and metabolism showed that distantly related microorganisms exhibited similar metabolic capabilities in ME and MT (Fig. 2D). To explore the relationship between metabolic capabilities and the taxa that perform them, we first mapped metabolic genes of all MAGs above middle-quality to metabolic pathways. Furthermore, we searched representative genes of target pathways against MAGs to obtain their taxonomic distribution in ME and MT (see the “Methods” section). We found surprising commonalities of metabolic capabilities between ME and MT, and distinct but various taxa could perform most metabolic pathways. We describe the essential functions that could be related to elemental cycling directly and adaptation to extreme stresses in the following paragraphs.

Carbon cycling

Degradation of refractory organic matters (ROMs)

Prokaryotes in both ME and MT exhibited metabolic versatility to utilize complex organic matters, especially for those forms that are generally considered to be refractory, including complex polymers (e.g., cellulose, lignin, and chitin), aromatic compounds (e.g., benzene, aromatic hydrocarbons, aromatic amino acids, and other derivatives), alkanes (short and long chain), and abiotic enantiomers (D-amino acids and L-sugars). Utilization of these ROMs relied on enzymes to break the complex or rigid structures of ROMs at the first step and then through a series of enzymes to produce intermediates that are metabolic nodes linked to central carbon metabolism (Fig. 4A). Most representative genes involved in ROM degradation were not restricted to certain species or clades but had a wide distribution among taxa in both ME and MT (Fig. 5A, Additional file 1: Results, Additional file 3: Table S2). Notably, for D-amino acid (D-AA) metabolism, the catabolic capacities for nine D-AAs, including two typical D-AAs in bacterial cell walls (i.e., D-Ala, D-Glu) and seven non-canonical D-AAs (i.e., D-Gln, D-Ser, D-Thr, D-Arg, D-Pro, D-Phe and D-Cys), were found in MAGs from both ME and MT. Among them, dadA, the typical gene for a broad-spectrum D-AA dehydrogenase in D-AA catabolism [62], was abundant and found in 445 (of 1176) MAGs (59% in ME and 35% in MT). Almost all identified class-level clades (81/84) contain at least one D-AA metabolic gene(s) (Fig. 5A, Additional file 1: Fig. S9A, Additional file 3: Table S2). Like the observations that the same metabolic genes occurred in various taxa, the degradation of complex polymers, aromatic compounds, alkanes, and L-sugars were also found among varied taxa (Figs. 4A and 5A, Additional file 3: Table S2). This result suggested a wide distribution of ROM utilization capacity in prokaryotes of both ME and MT.

Fig. 4
figure 4

Schematic diagram of prokaryotic metabolic versatile in ME and MT. A Refractory organic matter utilization. B Metabolism of nitrogen, sulfur, heavy metal, and ROS. C Distinct pathways in ME and MT. Arrows represent biochemical reactions between two compounds. Double arrows indicate that multiple reactions are involved in the conversion between two shown compounds. Key genes in each pathway are shown in italic, which are corresponding to the gene IDs in Additional file 3: Table S2

Fig. 5
figure 5

Distribution and co-occurrence of selected metabolic genes represented various biofunctions in MAGs. A Percentages of selected metabolic genes across class level in MAGs found in ME only, MT only, and both. Percentages of genomes of each class were shown in four levels: > 50% in dark red, 20–50% in orange, < 20% in light orange, while none in white. All MAGs used are above middle-quality (> 50% completeness and < 10% contamination). Key genes are corresponding to the gene IDs in Fig. 4 and Additional file 3: Table S2. B Co-occurrence network of selected metabolic genes in MAGs of ME or MT. All nodes with edges that binary-Jaccard distance > 0.3 are presented with circles. Colors of each node are corresponding to the color bar of biofunctions for selected metabolic genes in A

CO2 fixation

Although evidence for autotrophic CO2 fixation to produce organic matters (OMs) was detected in both ME and MT, CO2 fixation-associated pathways were specific in ME and MT, (Figs. 4C, 5A, and 6, Additional file 3: Table S2). Genes involved in the typical chemoautotrophic pathways were only detected in metagenomes and MAGs of MT, including the Wood–Ljungdahl (WL) pathway (cdhA and cdhC), reductive Krebs (rTCA) cycle (aclA, aclB, ccsA, ccsB, and ccl), and 3-hydroxypropionate/4-hydroxybutyrate (3HP-4HB) cycle (3hpcs and 4hbcl) (Additional file 1: Fig. S9A, Additional file 3: Table S2). Although gene for the Calvin–Benson–Bassham (CBB) cycle (prkB) was observed in both ME and MT, a series of genes in the photosystems I and II (psaA-X and psbA-Z) and allophycocyanin (AP, apcA-F) system, which is essential for photosynthesis, were only relatively complete in the ME-specific clade Cyanobacteria (Additional file 3: Table S2).

Fig. 6
figure 6

Schematic representation of major prokaryotic biofunctions compared between the ME and MT. ME-specific pathways are in red, while MT-specific pathways are in blue. OM, organic matter. ROM, refractory organic matter

Nitrogen and sulfur cycling

Although our understanding of prokaryote-mediated nitrogen and sulfur cycling is still incomplete for both ME and MT, our results show that most of the known valence transition processes between the essential compounds of nitrogen and sulfur can be observed in the prokaryotic MAGs of both ME and MT (Fig. 4B). Detailed processes in nitrogen and sulfur metabolism and the taxa that perform them at ME and MT are described as follows.

Nitrogen metabolism

The pathways of dissimilatory nitrate reduction to ammonia (napA/nasA/narB and nirA/nirB/nrfA) and denitrification (nirK, norB and nosZ) are both complete in prokaryotic MAGs from ME and MT, which suggested the stable capacities of prokaryotic communities to use nitrate and nitrite as electron acceptors in ME and MT (Fig. 4B, Supplemental Text). The dissimilatory nitrate reduction was complete in 57% of the total MAGs in ME and 19% in MT, and the denitrification was only complete in 2 MAGs in ME and 7 MAGs in MT. These nitrate/nitrite respiratory genes are widely distributed among ME-specific, MT-specific, and cross-habitat class-level clades (Fig. 5A, Additional file 1: Fig. S9B, Additional file 3: Table S2). In addition, the representative gene of the cleavage system from organic nitrogen to ammonium, the gcvT gene for aminomethyltransferase, and the ammonium transporter gene amt had relatively high abundance in both ME and MT. Both gcvT and amt can be observed in more than 70% of MAGs covered in 65 class-level clades, which indicates the potential importance of recycling and exchanging nitrogen sources between cells and inorganic ammonium in both ME and MT.

Compared to prokaryotic functions in MT, the nitrogen cycling in ME is relatively incomplete, i.e., the aerobic ammonia oxidation pathway is missing according to both read-based and MAG-based analyses (Figs. 4B, 5A, and Fig. 6, Additional file 1: Fig. S9B). This is caused by the absence of the typical ammonia oxidation archaea (AOA), Thaumarchaeota (class Nitrososphaeria), from the prokaryotic community of ME. The representative gene pmoA-amoA for ammonia oxidation and genes 3hpcs and 4hbcl for the 3HP-4HB cycle, the specific CO2 fixation pathway in Thaumarchaeota, were only found in the MAGs of Thaumarchaeota from MT (Fig. 5A). Interestingly, the hao gene encoding hydroxylamine dehydrogenase was not only associated with the AOA but also found in 25 other MT-specific classes, which supported the previous assumption that the hao gene can be involved in other nitrogen processes in additional to aerobic ammonia oxidation [63].

Sulfur metabolism

Most genes in known sulfur metabolic processes can be found in MAGs of both ME and MT, including pathways of sulfate reduction (sat, cysND, cysC, cysH for assimilation, and aprA for dissimilation), assimilatory sulfite reduction (cysJ, cysI, and sir), sulfur oxidation (SOX) (soxA, soxC, soxY, and glpE), thiosulfate reduction (phsA), and sulfide oxidation to sulfur (sqr and sdh) (Fig. 4B, Additional file 1: Results, Fig. S9B).

However, the dsrA gene involved in converting sulfite to sulfide in dissimilatory sulfate reduction was only present in 12 MT-specific MAGs, restricted to two classes, Desulfobacteria and Thermodesulfovibrionia (Figs. 4B and 5A, Additional file 3: Table S2). Desulfobacteria is the typical sulfate-reducing bacteria (SRB) in subseafloor sediments [64], while Thermodesulfovibrionia is a recently reported sulfur-nitrogen coupling bacteria [65]. These results indicated a potentially essential role in the dissimilatory sulfate reduction in MT, but further investigation is still needed.

Heavy metal metabolism

Heavy metals, such as arsenate-As (V), arsenite-As (III), and selenium (Se), were reported to accumulate together with the sinking of POM by V-topology in MT [28, 66] and were also determined to be present in snow collected from the East Rongbuk Glacier in ME [67]. Therefore, it is reasonable to suppose that prokaryotes in both ME and MT possibly contained metabolic capabilities for heavy metal utilization and/or detoxication. Our MAG-based results supported this assumption. Indeed, many MAGs in both ME and MT with various class-level taxa harbored arsenic-related genes encoding arsenic detoxication (e.g., arsC1 and arsC2), transport (e.g., ArsB and Acr3), and metabolism (e.g., arsM and aoxA) (Figs. 4B and 5A, Additional file 3: Table S2). Among them, representative genes for arsenate reduction from As(V) to As(III), arsC1 (glutaredoxin) and arsC2 (thioredoxin), were the most abundant among As-related genes in both ME and MT (Additional file 1: Fig. S9C). These two arsC genes were found in ~50% of total MAGs in either ME or MT. More than 80% of class-level clades harbored at least one arsC gene(s), indicating that the detoxication of As(V) is essential for prokaryotes in both ME and MT (Figs. 4B and 5A, Additional file 3: Table S2). Additionally, arsenic metabolism genes, e.g., arsM and aoxA (named as aioA in some studies), which are predicted to be less common than the arsC genes but highly endemic according to global soil microbiome analysis [68], were also found in both ME and MT. The aoxA gene for arsenite oxidation from As(III) to As(V) was found in a similar proportion (~ 10%) of total MAGs in both ME and MT. However, the assimilatory gene arsM for producing methylarsenite from arsenite [69] was found to account for higher coverage of total MAGs in MT (51%) than ME (8%) (Additional file 1: Fig. S9C, Additional file 3: Table S2).

Selenium (Se) is a trace element involved in the biosynthesis of selenium-containing amino acids, e.g., selenocysteine (Sec), most of which serve oxidoreductase functions [70]. Biosynthesis of Sec is not a common trait of prokaryotes [71]. Considering that the accumulation of Se (up to ~ 200 ng/g sample) has been reported in sediments of the bottom-axis of MT in a previous study [28], we found a relatively high abundance of representative Se-related genes among reads in MT as expected. Unexpectedly, Sec-related genes were also observed in ME, with a similar relative abundance to that in MT (Additional file 1: Fig. S9C), although the Se concentration in ME was measured at trace levels (45 ng/g sample) [72]. The selenophosphate synthetase gene selD, which is used to define the overall Se utilization from environments [71], was both overserved in 40% of total MAGs in ME and 47% in MT, being widely distributed in > 50% of class-level clades. A similar abundance of selA, selB, and selU genes encoding the synthesis of selenium-containing amino acids was observed in both ME and MT (Figs. 4B and 5A, Additional file 3: Table S2).

Stress resistance

Compatible solutes

The uptake and metabolism of compatible solutes is a typical strategy for stress adaptation across three domains of life. We found that these functional genes were shared in prokaryotes in both ME and MT (Figs. 4C and 5A, Additional file 1: Fig. S9C, Additional file 3: Table S2). Among all types of compatible solutes, glycine betaine (GB) is probably the most frequently used across three domains of life. The As-containing homolog of GB, arsenobetaine (As-GB), was reported to share the uptake and biosynthesis genes with GB and play a cytoprotective role against extremes in osmotic pressure and low temperature [73]. In our results, ~40% of total MAGs in ME and ~25% in MT harbored the representative gene opuC encoding transporters for both GB and As-GB, and ~28% of MAGs in ME and 33% in MT contained the representative gene betB encoding (arseno)betaine-aldehyde dehydrogenase to synthesize GB or/and As-GB. Both opuC and betB were distributed in various class-level clades (36%), including ME-specific, MT-specific and cross-habitat clades, indicating that genes for both the uptake and biosynthesis of GB and As-GB are widely distributed in prokaryotes in ME and MT.

In addition to the common compatible solutes mentioned above, prokaryotes in marine environments and deep-sea environments can synthesize or utilize some special compatible solutes for preserving macromolecules from high salinity and HHP environments, usually called as osmolytes or piezolytes. Trimethylamine N-oxide (TMAO) is a general piezolyte that commonly functions in both macro- and microorganisms from the deep sea [24, 74]. Interestingly, both biosynthesis and utilization genes for TMAO were found in MAGs from ME and MT (Figs. 4C, 5A, and 6, Additional file 1: Fig. S9C, Additional file 3: Table S2). The TMAO-producing gene tmm was encoded by 36 MAGs of MT (Alphaproteobacteria, Gammaproteobacteria and UBA9042) but was only found in one MAG of ME (Alphaproteobacteria). The relative abundance of the tmm gene among all clean reads in MT was significantly higher than that in ME (Wilcox’s test, p value < 0.001). The TMAO utilization gene torZ required to use TMAO as an anaerobic electron acceptor to produce trimethylamine (TMA) was also found in both ME and MT, with similar abundance (Additional file 1: Fig. S9C). However, the torZ gene had a wider distribution than the tmm gene, which can be found in 11 class-level clades. Among these classes, > 50% of MAGs belong to Fimbriimonadia, Thermoleophilia, and SORT01.

In contrast, dimethylsulfoniopropionate (DMSP)-related genes were restricted to prokaryotes in MT but not found in either the metagenome or MAGs from ME. DMPS is one of the most abundant organosulfur compounds in marine environments and has been reported to be produced in considerable amounts by bacterial cells in MT and act as a compatible solute in bacterial cells [75,76,77]. Representative genes involved in the biosynthesis (dsyB), demethylation (dmdA), and cleavage (dddL) of DMSP were only observed in MAGs of MT (Figs. 4C, 5A, and 6, Additional file 1: Fig. S9C, Additional file 3: Table S2). Among them, nine MAGs belonging to Alphaproteobacteria contained the DMSP-producing gene dsyB; 43 MAGs belonging to Alphaproteobacteria, Gammaproteobacteria, and Bacterioidia contained the DMSP demethylation gene dmdA; 14 MAGs belonging to Alphaproteobacteria and Gammaproteobacteria contained the DMSP cleavage gene dddL [78]. This result suggested a specific essential role of DMSP in MT ecosystems.

Antioxidation of reactive oxygen species (ROS)

Antioxidation was another shared pathway for stress resistance in prokaryotes across the ME and MT. Genes for typical ROS antioxidant enzymes, including superoxide dismutases, i.e., sod1 (Cu-Zn family), sod2 (Fe-Mn family), and sodN (Ni family), and catalase, i.e., katE and katG, were present in both ME and MT (Figs. 4B and 5A, Additional file 3: Table S2). Among them, the relative abundance of the katE and sod2 genes in ME was significantly higher than that in MT; however, the Ni-sodN in MT was slightly higher than that in ME (Wilcox’s test, p value < 0.05) (Additional file 1: Fig. S9C). All of them were widely distributed in MAGs from ME and MT. The sod2 gene had the highest coverage among all these tested ROS antioxidant genes, at ~91% of MAGs in ME and ~56% of MAGs in MT. Tested genes for ROS antioxidation were present at least one in ~80% of identified class-level clades. Our observations illustrated that ROS antioxidation is common and essential for prokaryotes regardless of whether they live at ME or MT.

Distinct cooccurrence of metabolic functions across MAGs compared between ME and MT

Interestingly, the cooccurrence of different metabolic genes was quite different across prokaryotic genomes in ME and MT (see the “Methods” section, Fig. 5B, Additional file 1: Fig. S10). Generally, a closer correlation of the central clustered genes was observed in ME, while the cooccurrence of different genes in MT was more loosely correlated. Notably, in the cooccurrence network of ME, the degradation genes involved in aromatic compounds (benA, catA, pcaH), alkane (ladA), and D-AAs (dsdA) were correlated with each other and linked with the genes related to the universal compatible solute GB (gbsA and opuC) as nodes, and the degradation genes of complex polymers were correlated with functional genes involved in nitrogen, sulfur, and heavy metal metabolism. In contrast, there seems to be no such a cooccurrence of genes related to ROM utilization (except D-AA-related genes) among themselves or genes involved in other tested metabolic functions. In the cooccurrence network of MT, genes related to chemotrophic CO2 fixations were separated from the central cluster. For example, the gene pmoA-amoA for ammonia oxidation and genes 3hpcs and 4hbcl for the 3HP-4HB pathway, specific to Thaumarchaeota from MT, formed a separated cluster. A similar situation occurred for the cdhC and cdhA genes for the WL pathway in MT. The above results indicated that the combination of different metabolic modules, especially those related to different elemental cycling, may be specific to ME or MT. Despite the overall distinct cooccurrence network, regardless of ME or MT, the metabolic genes related to nitrogen, sulfur, As, and Se were highly correlated. This result indicated potential crosstalk among cycling of different elements in both ME and MT. In addition, the ROS antioxidation genes and GB biosynthesis gene (gbsA) were also observed in the central cooccurrence cluster, correlated with genes related to N, S, As, and Se, which suggested that prokaryotes in both ME and MT play an essential role in ROS antioxidation and biosynthesis of the universal compatible solute GB.


In this study, we investigated and compared prokaryotic communities by performing metagenome analysis between two geographically isolated extreme environments, ME and MT, the highest and deepest places on Earth, respectively. We observed significant differences in taxonomy between these two habitats, with many specific taxa that were only found in ME or MT. For taxonomy, our result that taxa tend to adapt and be specific to their habitats is consistent with previous global studies [6]. As the community composition is shown (Additional file 1: Fig. S1B), Gammaproteobacteria were the most abundant group in most samples of both ME and MT, which were reported to account for ~ 30% of glaciers [79, 80], ~60% in hadal trenches [81,82,83,84], and widely distributed in other cold environments, such as the Antarctic lakes and salt-cones [85, 86]. The shared genera between ME and MT, such as Acinetobacter and Sphingomonas, were also reported as the dominant group in cold environments, including worldwide glacier environments [87,88,89,90,91], clouds [92], and seawater [93]. In addition, the taxa specific to ME or MT we observed in this study are generally consistent with the typical prokaryotic clades in glaciers or trench environments, respectively. For example, the ME-specific phylum, Cyanobacteria, was reported as the pioneer colonizers of the glaciers, inhabiting ice surfaces [94], snow [95], and other glacier habitats [96, 97], playing a keystone role as the leading primary producer in the microbial food web in glacier environments (e.g., photosynthetic carbon and nitrogen fixation) [98,99,100,101]. Another ME-specific phylum we found in this study is Deinococcota, which is usually reported as radiation-resistant bacteria [21], and has been observed in ME-related glaciers under high UV conditions [12]. Similarly, the prokaryotic community composition observed in this study is quite comparable to other studies on the microbiome in hadal trenches. The MT-specific archaeal phylum, Nitrososphaeria (AOA), was the most abundance archaea in this study, which were also reported as the dominant archaea group with up to 30% in both water and sediment samples in hadal trenches [27, 102, 103].

This taxonomic variability revealed in this study provides clues to key environmental factors that affect prokaryotic life mostly in ME and MT. Our comparison revealed a significant difference in the genomic and predicted growth features of prokaryotes between ME and MT (Fig. 3, Additional file 1: Fig. S7). Although both ME and MT are low-temperature environments, the in situ temperature in MT is relatively stable (approximately 2 °C), while the temperature of ME, as mentioned above, is much more changeable (span 29 °C, ranging from −19.3 to 9.4 °C in summer) and varies with seasons and day/night cycles, which is also observed in other glacial ecosystem [104, 105]. The tolerance to temperature variations of prokaryotes in ME possibly requires special mechanisms with higher energy costs than adaptation to a stable low-temperature environment in MT, which may cause different genome size, GC contents, and growth features of prokaryotes. Moreover, growth temperature was reported to be closely related to mutation rate and population size in previous studies, especially under temperature fluctuations [106, 107]. Therefore, the difference between stable versus changeable in situ temperatures may cause distinct specialization or environmental selection mechanisms and processes in prokaryotes of ME and MT, which may be further reflected in the taxonomic diversity and composition of prokaryotic communities. Furthermore, bacterial genomes are commonly considered to tend to increase in size by aggregating adaptive gene modules to provide greater metabolic flexibility [108]. As alpha diversity showed above, prokaryotic communities in ME had lower diversity at a high taxonomic level (phylum) but had higher diversity at a low taxonomic level (species) than those in MT. The significantly higher genome size of prokaryotes in ME may be an explanation for why the prokaryotic community in ME has a lower taxonomic diversity than MT, but further investigation is needed. Finally, the cooccurrence network of versatile metabolic genes between ME and MT indicated that although the metabolic capabilities of the prokaryotic community are similar, distinct combinations of metabolic pathways are performed by different taxa, which could be a specific strategy for adaptation to distinct habitats. Although most metabolic pathways are shared between ME and MT and distributed among various taxa, some specialized functions are still restricted within certain taxa (i.e., AOA and SRB) and are only found in MT. These specific functions and the taxa that perform them indicated a potential essential role for ammonia oxidation in the nitrogen cycle and dissimilatory sulfate reduction in the sulfur cycle in MT.

In contrast to the taxonomy, to our surprise, prokaryotic metabolic capabilities exhibited more commonalities across these two extreme habitats, with > 90 metabolic modules overlapping. The Sankey network and taxonomic distribution of representative functional genes confirmed the inference of functional redundancy in which different taxa possibly harbor the same metabolic functions in ME and MT (Figs. 2D and 6). Although this functional redundancy phenomenon has been observed in some host-associated environments and bioreactors [11, 109,110,111,112], here is the first observation across such two distinguished and isolated extreme environments. The high species variability and metabolic commonalities between ME and MT make us rethink the relationship between taxonomy, metabolism, diversity, and adaptation of prokaryotes. These results may shed light on the principle of prokaryotic diversity: although taxa are specific to their habitats, primary metabolic functions could be still conserved.

These surprising metabolic commonalities between prokaryotes in ME and MT indicated underlying common strategies for prokaryotic adaptation to different extreme environmental conditions, although the geological topography and environmental conditions (e.g., sunlight, pressure, source and composition of nutrients) are intuitively different. For example, the stress of intense UV radiation in ME and the HHP in MT are seemingly incomparable at first glance. However, both UV and HHP can cause peroxidation damage for cells [22, 113], which makes antioxidants essential for prokaryotic adaptation in both ME and MT. This can explain why ROS antioxidation is common and widespread among diverse taxa of prokaryotes from both ME and MT (Fig. 5). Another example is the utilization of arsenic and selenic compounds, which is common in both ME and MT (Figs. 4, 5, and 6). A recent study reported a relatively high concentration (μg level) of As and Se in sediment samples from the bottom-axis of MT and pointed out the potentially important role of prokaryotes in redox balancing [28]. However, the concentration of As and Se in ME was reported to be two orders of magnitude lower than that in MT [72], and a similar relative abundance of genes enabling utilization of As and Se compounds was observed in ME in this study, which indicates a possible general requirement of cells for these two heavy metals and functions in their cycling [114,115,116]. Furthermore, both ME and MT can be considered nutrient-poor environments. The versatile metabolic genes encoding the ROM utilization enzymes in both ME and MT suggested a common strategy to survive under nutrient-poor environments where the ROM usually composes a more significant proportion of the organic matter than under normal conditions (e.g., in soil and surface seawater) [117]. Generally, a similar metabolic backbone and the basic biological rules are the basis of the common strategies for prokaryotes regardless of where they live.


Although it is too early to draw a clean conclusion on how prokaryotic diversity formed and which factors influenced it, our results already suggest a possible relationship between prokaryotic taxonomy and metabolism: habitats tend to select the taxonomic composition of the prokaryotic community, while their metabolic functions could be very similar. This conclusion is valid at least in the comparison of prokaryotes between the highest and deepest habitats on Earth, and whether this conclusion could be extended to a broader scale still requires more investigations with a combination of a massive number of environmental parameter measurements, geochemical data, and multiple omics data for prokaryotes.

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available in the NCBI repository in the accession number of projects PRJNA813429 and PRJNA859662. A full record of all statistical analysis is included as Additional file 4 and was created using the knitr package in R. All original scripts and result tables are available in GitHub ( Corresponding taxonomic classifications and information of metagenomic assembled genomes with annotations have all been included as Additional files 2 and 3, respectively.



D-Glycerate 2-phosphate


3-Hydroxypropionate/4-hydroxybutyrate pathway




Acetyl phosphate




Average nucleotide identity


Ammonia oxidation archaea




Beta-D-fructose 6-phosphate


Calvin-Benson-Bassham cycle


D-Amino acid


Dihydroxyacetone phosphate


Dimethylarsinic acid




Beta-D-fructose 1:6-bisphosphate


D-Glycerate 1:3-diphosphate


D-Glucose 6-phosphate


Glyceraldehyde 3-phosphate


Glycine betaine


Genome Taxonomy database


High hydrostatic pressure


Kyoto Encyclopedia of Genes and Genomes


KEGG orthology


Mount Everest


Monomethylarsonic acid


Molecular Operational Taxonomic Unit


The Mariana trench


Nonmetric multidimensional scaling


Optimal growth temperature


Organic matter




Particle organic matter


Practical Salinity Unit


Refractory organic matter


Reactive oxidative species


Reductive tricarboxylic acid cycle


Sulfate-reducing bacteria


Sulfur oxidation




Trimethylamine N-oxide






  1. Shu WS, Huang LN. Microbial diversity in extreme environments. Nat Rev Microbiol. 2021;20:219–35.

    Article  Google Scholar 

  2. Rothschild LJ, Mancinelli RL. Life in extreme environments. Nature. 2001;409:1092–101.

    Article  CAS  Google Scholar 

  3. Schmid AK, Allers T. DiRuggiero J. SnapShot: microbial extremophiles. Cell. 2020;180:818.

    Article  CAS  Google Scholar 

  4. Liu YQ, Ji MK, Yu T, Zaugg J, Anesio AM, Zhang ZH, et al. A genome and gene catalog of glacier microbiomes. Nat Biotechnol. 2022.

  5. Louca S, Polz MF, Mazel F, Albright MBN, Huber JA, O’Connor MI, et al. Function and functional redundancy in microbial systems. Nat Ecol Evol. 2018;2:936–43.

    Article  Google Scholar 

  6. Coelho LP, Alves R, del Río ÁR, Myers PN, Cantalapiedra CP, Giner-Lamia J, et al. Towards the biogeography of prokaryotic genes. Nature. 2022;601:252–6.

    Article  CAS  Google Scholar 

  7. Hug LA, Baker BJ, Anantharaman K, Brown CT, Probst AJ, Castelle CJ, et al. A new view of the tree of life. Nat Microbiol. 2016;1:16048.

    Article  CAS  Google Scholar 

  8. Torsvik V, Øvreås L, Thingstad TF. Prokaryotic diversity--magnitude, dynamics, and controlling factors. Science. 2002;296:1064–6.

    Article  CAS  Google Scholar 

  9. von Mering C, Hugenholtz P, Raes J, Tringe SG, Doerks T, Jensen LJ, et al. Quantitative phylogenetic assessment of microbial communities in diverse environments. Science. 2007;315:1126–30.

    Article  Google Scholar 

  10. Maistrenko OM, Mende DR, Luetge M, Hildebrand F, Schmidt TSB, Li SS, et al. Disentangling the impact of environmental and phylogenetic constraints on prokaryotic within-species diversity. ISME J. 2020;14:1247–59.

    Article  Google Scholar 

  11. Liu YQ, Yao TD, Kang SC, Jiao NZ, Zeng YH, Huang SJ, et al. Microbial community structure in major habitats above 6000 m on Mount Everest. Chin Sci Bull. 2007;52:2350–7.

    Article  Google Scholar 

  12. Moore K, Semple J, Cristofanelli P, Bonasoni P, Stocchi P. Environmental conditions at the South Col of Mount Everest and their impact on hypoxia and hypothermia experienced by mountaineers. Extrem Physiol Med. 2012;1:2–2.

    Article  Google Scholar 

  13. Merino N, Aronson HS, Bojanova DP, Feyhl-Buska J, Wong ML, Zhang S, et al. Living at the extremes: extremophiles and the limits of life in a planetary context. Front Microbiol. 2019;10:780.

    Article  Google Scholar 

  14. Zhang SH, Hou SG, Yang GL, Wang JH. Bacterial community in the East Rongbuk Glacier, Mt. Qomolangma (Everest) by culture and culture-independent methods. Microbiol Res. 2010;165:336–45.

    Article  CAS  Google Scholar 

  15. Liu YQ, Fang PC, Guo BX, Ji MK, Liu PF, Mao GN, et al. A comprehensive dataset of microbial abundance, dissolved organic carbon, and nitrogen in Tibetan Plateau glaciers. Earth Syst Sci Data. 2022;14:2303–14.

    Article  Google Scholar 

  16. Liu YQ, Yao TD, Jiao NZ, Tian LD, Hu AY, Yu WS, et al. Microbial diversity in the snow, a moraine lake and a stream in Himalayan glacier. Extremophiles. 2011;15:411–21.

    Article  Google Scholar 

  17. Liu YQ, Yao TD, Jiao NZ, Kang SC, Zeng YH, Huang SJ. Microbial community structure in moraine lakes and glacial meltwaters, Mount Everest: Microbial diversity in lakes of Mount Everest. FEMS Microbiol Lett. 2006;265:98–105.

    Article  CAS  Google Scholar 

  18. Morvan T, Lemoine C, Gaillard F, Hamelin G, Trinkler B, Carteaux L, et al. A comprehensive dataset on nitrate, Nitrite and dissolved organic carbon leaching losses from a 4-year Lysimeter study. Data Brief. 2020;32:106029.

    Article  Google Scholar 

  19. Liu YQ, Yao TD, Kang SC, Jiao NZ, Zeng YH, Shi Y, et al. Seasonal variation of snow microbial community structure in the East Rongbuk glacier. Mt Everest Sci Bull. 2006;51:1476–86.

    Article  CAS  Google Scholar 

  20. Battista JR. AGAINST ALL ODDS: The survival strategies of Deinococcus radiodurans. Annu Rev Microbiol. 1997;51:203–24.

    Article  CAS  Google Scholar 

  21. Chen P, Zhou H, Huang YY, Xie Z, Zhang MJ, Wei YL, et al. Revealing the full biosphere structure and versatile metabolic functions in the deepest ocean sediment of the Challenger Deep. Genome Biol. 2021;22:207.

    Article  CAS  Google Scholar 

  22. Xiao X, Zhang Y, Wang FP. Hydrostatic pressure is the universal key driver of microbial evolution in the deep ocean and beyond. Environ Microbiol Rep. 2021;13:68–72.

    Article  Google Scholar 

  23. Yancey PH. Cellular responses in marine animals to hydrostatic pressure. J Exp Zool A Ecol Integr Physiol. 2020;333:398–420.

    Article  CAS  Google Scholar 

  24. Qin QL, Wang ZB, Su HN, Chen XL, Miao J, Wang XJ, et al. Oxidation of trimethylamine to trimethylamine N -oxide facilitates high hydrostatic pressure tolerance in a generalist bacterial lineage. Sci Adv. 2021;7:eabf9941.

    Article  CAS  Google Scholar 

  25. Yin QJ, Zhang WJ, Qi XQ, Zhang SD, Jiang T, Li XG, et al. High hydrostatic pressure inducible trimethylamine N-oxide reductase improves the pressure tolerance of piezosensitive bacteria Vibrio fluvialis. Front Microbiol. 2017;8:2646.

    Article  Google Scholar 

  26. Nunoura T, Nishizawa M, Hirai M, Shimamura S, Harnvoravongchai P, Koide O, et al. Microbial diversity in sediments from the bottom of the Challenger Deep, the Mariana Trench. Microbes Environ. 2018;33:186–94.

    Article  Google Scholar 

  27. Nunoura T, Takaki Y, Hirai M, Shimamura S, Makabe A, Koide O, et al. Hadal biosphere: insight into the microbial ecosystem in the deepest ocean on earth. Proc Natl Acad Sci U S A. 2015;112:E1230–6.

    Article  CAS  Google Scholar 

  28. Zhou YL, Mara P, Cui GJ, Edgcomb VP, Wang Y. Microbiomes in the Challenger Deep slope and bottom-axis sediments. Nat Commun. 2022;13:1515.

    Article  CAS  Google Scholar 

  29. Luo M, Gieskes J, Chen LY, Shi XF, Chen DF. Provenances, distribution, and accumulation of organic matter in the southern Mariana Trench rim and slope: implication for carbon cycle and burial in hadal trenches. Mar Geol. 2017;386:98–106.

    Article  CAS  Google Scholar 

  30. Glud RN, Wenzhöfer F, Middelboe M, Oguri K, Turnewitsch R, Canfield DE, et al. High rates of microbial carbon turnover in sediments in the deepest oceanic trench on Earth. Nat Geosci. 2013;6:284–8.

    Article  CAS  Google Scholar 

  31. Hodson A, Anesio AM, Tranter M, Fountain A, Osborn M, Priscu J, et al. Glacial ecosystems. Ecol Monogr. 2008;78:41–67.

    Article  Google Scholar 

  32. Takai K. Recent topics on deep-sea microbial communities in microbes and environments. Microbes Environ. 2019;34:345–6.

    Article  Google Scholar 

  33. Natarajan VP, Zhang XX, Morono Y, Inagaki F, Wang FP. A modified sds-based DNA extraction method for high quality environmental DNA from seafloor environments. Front Microbiol. 2016;7:986.

    Article  Google Scholar 

  34. Lang JD, Zhu RR, Sun X, Zhu SY, Li TB, Shi XL, et al. Evaluation of the MGISEQ-2000 sequencing platform for Illumina target capture sequencing libraries. Front Genet. 2021;12:730519.

    Article  CAS  Google Scholar 

  35. Joshi NA, Fass JN. Sickle: A sliding-window, adaptive, quality-based trimming tool for FastQ files (Version 1.33). 2011.

  36. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32:3047–8.

    Article  CAS  Google Scholar 

  37. Gruber-Vodicka HR, Seah BKB, Pruesse E. phyloFlash: rapid small-subunit rRNA profiling and targeted assembly from metagenomes. MSystems. 2020;5:e00920.

    Article  CAS  Google Scholar 

  38. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.

    Article  CAS  Google Scholar 

  39. Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.

    Article  CAS  Google Scholar 

  40. Bushnell B. BBMap: a fast, accurate, splice-aware aligner. 2014.

  41. Hyatt D, Chen GL, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.

    Article  Google Scholar 

  42. Hyatt D, LoCascio PF, Hauser LJ, Uberbacher EC. Gene and translation initiation site prediction in metagenomic sequences. Bioinformatics. 2012;28:2223–30.

    Article  CAS  Google Scholar 

  43. Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36:2251–2.

    Article  CAS  Google Scholar 

  44. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.

    Article  CAS  Google Scholar 

  45. Kang DD, Li F, Kirton ES, Thomas A, Egan RS, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies; 2019.

    Book  Google Scholar 

  46. Wu YW, Simmons BA, Singer SW. MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics. 2016;32:605–7.

    Article  CAS  Google Scholar 

  47. Sieber CMK, Probst AJ, Sharrar A, Thomas BC, Hess M, Tringe SG, et al. Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat Microbiol. 2018;3:836–43.

    Article  CAS  Google Scholar 

  48. 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:2864–8.

    Article  CAS  Google Scholar 

  49. 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:725–31.

    Article  CAS  Google Scholar 

  50. Parks DH, Chuvochina M, Chaumeil P-A, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol. 2020;38:1079–86.

    Article  CAS  Google Scholar 

  51. Chaumeil PA, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2019;36(6):1925–7.

  52. Woodcroft BJ. CoverM: read coverage calculator for metagenomics. 2021.

  53. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community ecology package. 2015.

    Google Scholar 

  54. Ahlmann-Eltze C, Patil I. ggsignif: R package for displaying significance brackets for ‘ggplot2’. 2021. PsyArXiv. 2021.

  55. Xie YH. knitr: a comprehensive tool for reproducible research in R. In: Stodden V, Leisch F, Peng RD, editors. Implementing reproducible research. Chapman: Hall/CRC; 2014. p. 3–31.

  56. Sorek R, Zhu YW, Creevey CJ, Francino MP, Bork P, Rubin EM. Genome-wide experimental determination of barriers to horizontal gene transfer. Science. 2007;318:1449–52.

    Article  CAS  Google Scholar 

  57. Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293–6.

    Article  CAS  Google Scholar 

  58. Vieira-Silva S, Rocha EPC. The systemic imprint of growth and its uses in ecological (meta)genomics. PLoS Genet. 2010;6:e1000808.

    Article  Google Scholar 

  59. Weissman JL, Hou SW, Fuhrman JA. Estimating maximal microbial growth rates from cultures, metagenomes, and single cells via codon usage patterns. Proc Natl Acad Sci U S A. 2021;118:e2016810118.

    Article  CAS  Google Scholar 

  60. Louca S. The rates of global bacterial and archaeal dispersal. ISME J. 2022;16:159–67.

    Article  Google Scholar 

  61. Rodriguez-R LM, Konstantinidis KT. Bypassing cultivation to identify bacterial species: culture-independent genomic approaches identify credibly distinct clusters, avoid cultivation bias, and provide true insights into microbial species. Microbe Magazine. 2014;9:111–8.

    Article  Google Scholar 

  62. Wang XY, Yang Y, Lv YX, Xiao X, Zhao WS. The capability of utilizing abiotic enantiomers of amino acids by Halomonas sp. LMO_D1 derived from the Mariana Trench. Front Astronomy Space Sci. 2021;8:741053.

  63. Maalcke WJ, Reimann J, de Vries S, Butt JN, Dietl A, Kip N, et al. Characterization of anammox hydrazine dehydrogenase, a key N2-producing enzyme in the global nitrogen cycle*. J Biol Chem. 2016;291:17077–92.

    Article  CAS  Google Scholar 

  64. Zhao R, Biddle JF. Helarchaeota and co-occurring sulfate-reducing bacteria in subseafloor sediments from the Costa Rica Margin. ISME Commun. 2021;1:25.

    Article  Google Scholar 

  65. Umezawa K, Kojima H, Kato Y, Fukui M. Dissulfurispira thermophila gen. nov., sp. nov., a thermophilic chemolithoautotroph growing by sulfur disproportionation, and proposal of novel taxa in the phylum Nitrospirota to reclassify the genus Thermodesulfovibrio. Syst Appl Microbiol. 2021;44:126184.

    Article  CAS  Google Scholar 

  66. Sun RY, Yuan JJ, Sonke JE, Zhang YX, Zhang T, Zheng W, et al. Methylmercury produced in upper oceans accumulates in deep Mariana Trench fauna. Nat Commun. 2020;11:3389.

    Article  CAS  Google Scholar 

  67. Lee K, Hur SD, Hou S, Hong S, Qin X, Ren J, et al. Atmospheric pollution for trace elements in the remote high-altitude atmosphere in central Asia as recorded in snow from Mt. Qomolangma (Everest) of the Himalayas. Sci Total Environ. 2008;404:171–81.

    Article  CAS  Google Scholar 

  68. Dunivin TK, Yeh SY, Shade A. A global survey of arsenic-related genes in soil microbiomes. BMC Biol. 2019;17:45.

    Article  CAS  Google Scholar 

  69. Neff JM. Ecotoxicology of arsenic in the marine environment. Environ Toxicol Chem. 1997;16:917–27.

    CAS  Google Scholar 

  70. Labunskyy VM, Hatfield DL, Gladyshev VN. Selenoproteins: molecular pathways and physiological roles. Physiol Rev. 2014;94:739–77.

    Article  CAS  Google Scholar 

  71. Peng T, Lin J, Xu YZ, Zhang Y. Comparative genomics reveals new evolutionary and ecological patterns of selenium utilization in bacteria. ISME J. 2016;10:2048–59.

    Article  CAS  Google Scholar 

  72. Kang SC, Zhang QG, Kaspari S, Qin DH, Cong ZY, Ren JW, et al. Spatial and seasonal variations of elemental composition in Mt. Everest (Qomolangma) snow/firn. Atmos Environ. 2007;41:7208–18.

    Article  CAS  Google Scholar 

  73. Hoffmann T, Warmbold B, Smits SHJ, Tschapek B, Ronzheimer S, Bashir A, et al. Arsenobetaine: an ecophysiologically important organoarsenical confers cytoprotection against osmotic stress and growth temperature extremes: Stress protection by uptake and synthesis of glycine betaine. Environ Microbiol. 2018;20:305–23.

    Article  CAS  Google Scholar 

  74. Wang K, Shen YJ, Yang YZ, Gan XN, Liu GC, Hu K, et al. Morphology and genome of a snailfish from the Mariana Trench provide insights into deep-sea adaptation. Nat Ecol Evol. 2019;3:823833.

    Article  Google Scholar 

  75. Liu JL, Liu J, Zhang SH, Liang JC, Lin HY, Song DL, et al. Novel insights into bacterial dimethylsulfoniopropionate catabolism in the East China Sea. Front Microbiol. 2018;9:3206.

    Article  Google Scholar 

  76. Zhang XH, Liu J, Liu J, Yang G, Xue CX, Curson ARJ, et al. Biogenic production of DMSP and its degradation to DMS—their roles in the global sulfur cycle. Sci China Life Sci. 2019;62:1296–319.

    Article  CAS  Google Scholar 

  77. Zheng YF, Wang JY, Zhou S, Zhang YH, Liu J, Xue CX, et al. Bacteria are important dimethylsulfoniopropionate producers in marine aphotic and high-pressure environments. Nat Commun. 2020;11:4658.

    Article  CAS  Google Scholar 

  78. Teng ZJ, Wang P, Chen XL, Guillonneau R, Li CY, Zou SB, et al. Acrylate protects a marine bacterium from grazing by a ciliate predator. Nat Microbiol. 2021;6:1351–6.

    Article  CAS  Google Scholar 

  79. Liu YQ, Yao TD, Jiao NZ, Kang SC, Xu BQ, Zeng YH, et al. Bacterial diversity in the snow over Tibetan Plateau Glaciers. Extremophiles. 2009;13:411–23.

    Article  CAS  Google Scholar 

  80. Cameron KA, Hagedorn B, Dieser M, Christner BC, Choquette K, Sletten R, et al. Diversity and potential sources of microbiota associated with snow on western portions of the Greenland Ice Sheet. Environ Microbiol. 2015;17:594–609.

    Article  CAS  Google Scholar 

  81. Tarn J, Peoples LM, Hardy K, Cameron J, Bartlett DH. Identification of free-living and particle-associated microbial communities present in hadal regions of the Mariana Trench. Front Microbiol. 2016;7:665.

    Article  Google Scholar 

  82. Peoples LM, Donaldson S, Osuntokun O, Xia Q, Nelson A, Blanton J, et al. Vertically distinct microbial communities in the Mariana and Kermadec trenches. PLoS One. 2018;13:e0195102.

    Article  Google Scholar 

  83. Li WL, Huang JM, Zhang PW, Cui GJ, Wei ZF, Wu YZ, et al. Periodic and spatial spreading of alkanes and Alcanivorax bacteria in deep waters of the Mariana Trench. Appl Environ Microbiol. 2019;85:e02089–18.

    Article  CAS  Google Scholar 

  84. Liu J, Zheng YF, Lin HY, Wang XC, Li M, Liu Y, et al. Proliferation of hydrocarbon-degrading microbes at the bottom of the Mariana Trench. Microbiome. 2019;7:47.

    Article  CAS  Google Scholar 

  85. Lo Giudice A, Conte A, Papale M, Rizzo C, Azzaro M, Guglielmin M, et al. Prokaryotic diversity and metabolically active communities in brines from two perennially ice-covered Antarctic lakes. Astrobiology. 2021;21:551–65.

    Article  CAS  Google Scholar 

  86. Azzaro M, Papale M, Rizzo C, Forte E, Lenaz D, Guglielmin M, et al. Antarctic salt-cones: an oasis of microbial life? The example of Boulder Clay Glacier (Northern Victoria Land). Microorganisms. 2022;10:1753.

    Article  CAS  Google Scholar 

  87. Mikucki JA, Priscu JC. Bacterial diversity associated with blood falls, a subglacial outflow from the Taylor Glacier, Antarctica. Appl Environ Microbiol. 2007;73:4029–39.

    Article  CAS  Google Scholar 

  88. Miteva VI, Brenchley JE. Detection and isolation of ultrasmall microorganisms from a 120,000-year-old Greenland glacier ice core. Appl Environ Microbiol. 2005;71:7806–18.

    Article  CAS  Google Scholar 

  89. Yang GL, Hou SG, Baoge RL, Li ZG, Xu H, Liu YP, et al. Differences in bacterial diversity and communities between glacial snow and glacial soil on the Chongce Ice Cap, West Kunlun Mountains. Sci Rep. 2016;6:36548.

    Article  Google Scholar 

  90. Liu Q, Liu HC, Zhang JL, Zhou YG, Xin YH. Sphingomonas psychrolutea sp nov., a psychrotolerant bacterium isolated from glacier ice. Int J Syst Evol Microbiol. 2015;65:2955–9.

    Article  CAS  Google Scholar 

  91. Zhong ZP, Tian FN, Roux S, Gazitua MC, Solonenko NE, Li YF, et al. Glacier ice archives nearly 15,000-year-old microbes and phages. Microbiome. 2021;9:160.

    Article  CAS  Google Scholar 

  92. DeLeon-Rodriguez N, Lathem TL, Rodriguez-R LM, Barazesh JM, Anderson BE, Beyersdorf AJ, et al. Microbiome of the upper troposphere: Species composition and prevalence, effects of tropical storms, and atmospheric implications. Proc Natl Acad Sci U S A. 2013;110:2575–80.

    Article  CAS  Google Scholar 

  93. Laskin AI, White DC. Preface to special issue on Sphingomonas. J Ind Microbiol. 1999;23:231.

    CAS  Google Scholar 

  94. Yallop ML, Anesio AM, Perkins RG, Cook J, Telling J, Fagan D, et al. Photophysiology and albedo-changing potential of the ice algal community on the surface of the Greenland ice sheet. ISME J. 2012;6:2302–13.

    Article  CAS  Google Scholar 

  95. Lutz S, Anesio AM, Edwards A, Benning LG. Microbial diversity on Icelandic glaciers and ice caps. Front Microbiol. 2015;6:307.

    Article  Google Scholar 

  96. Edwards A, Anesio AM, Rassner SM, Sattler B, Hubbard B, Perkins WT, et al. Possible interactions between bacterial diversity, microbial activity and supraglacial hydrology of cryoconite holes in Svalbard. ISME J. 2011;5:150–60.

    Article  Google Scholar 

  97. Conte A, Papale M, Amalfitano S, Mikkonen A, Rizzo C, De Domenico E, et al. Bacterial community structure along the subtidal sandy sediment belt of a high Arctic fjord (Kongsfjorden, Svalbard Islands). Sci Total Environ. 2018;619-620:203–11.

    Article  CAS  Google Scholar 

  98. Stibal M, Sabacka K, Kastovska K. Microbial communities on glacier surfaces in Svalbard: Impact of physical and chemical properties on abundance and structure of cyanobacteria and algae. Microb Ecol. 2006;52:644–54.

    Article  Google Scholar 

  99. Murakami T, Takeuchi N, Mori H, Hirose Y, Edwards A, Irvine-Fynn T, et al. Metagenomics reveals global-scale contrasts in nitrogen cycling and cyanobacterial light-harvesting mechanisms in glacier cryoconite. Microbiome. 2022;10:50.

    Article  CAS  Google Scholar 

  100. Hodson AJ, Mumford PN, Kohler J, Wynn PM. The High Arctic glacial ecosystem: new insights from nutrient budgets. Biogeochemistry. 2005;72:233–56.

    Article  CAS  Google Scholar 

  101. Franzetti A, Tagliaferri I, Gandolfi I, Bestetti G, Minora U, Mayer C, et al. Light-dependent microbial metabolisms drive carbon fluxes on glacier surfaces. ISME J. 2016;10:2984–8.

    Article  CAS  Google Scholar 

  102. Jing HM, Xiao X, Zhang Y, Li ZY, Jian HH, Luo YF, et al. Composition and ecological roles of the core microbiome along the abyssal-hadal transition zone sediments of the Mariana Trench. Microbiol Spectr. 2022;10(3):e0198821.

  103. Zhao JL, Jing HM, Wang ZM, Wang L, Jian HH, Zhang R, et al. Novel viral communities potentially assisting in carbon, nitrogen, and sulfur metabolism in the upper slope sediments of Mariana Trench. MSystems. 2022;7:e01358–21.

    Article  CAS  Google Scholar 

  104. Cruaud P, Vigneron A, Fradette MS, Dorea CC, Culley AI, Rodriguez MJ, et al. Annual bacterial community cycle in a seasonally ice-covered river reflects environmental and climatic conditions. Limnol Oceanogr. 2020;65:S21–37.

    Article  Google Scholar 

  105. Bourgeois S, Kerherve P, Calleja ML, Many G, Morata N. Glacier inputs influence organic matter composition and prokaryotic distribution in a high Arctic fjord (Kongsfjorden, Svalbard). J Mar Syst. 2016;164:112–27.

    Article  Google Scholar 

  106. Gu JH, Wang XJ, Ma XP, Sun Y, Xiao X, Luo HW. Unexpectedly high mutation rate of a deep-sea hyperthermophilic anaerobic archaeon. ISME J. 2021;15:1862–9.

    Article  CAS  Google Scholar 

  107. Chen ZY, Wang XJ, Song Y, Zeng QL, Zhang Y, Luo HW. Prochlorococcus have low global mutation rate and small effective population size. Nat Ecol Evol. 2022;6:183–94.

    Article  Google Scholar 

  108. Bobay L-M, Ochman H. The evolution of bacterial genome architecture. Front Genet. 2017;8:72.

    Article  Google Scholar 

  109. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, et al. A core gut microbiome in obese and lean twins. Nature. 2009;457:480–4.

    Article  CAS  Google Scholar 

  110. Burke C, Steinberg P, Rusch D, Kjelleberg S, Thomas T. Bacterial community assembly based on functional genes rather than species. Proc Natl Acad Sci U S A. 2011;108:14288–93.

    Article  CAS  Google Scholar 

  111. Wittebolle L, Vervaeren H, Verstraete W, Boon N. Quantifying community dynamics of nitrifiers in functionally stable reactors. Appl Environ Microbiol. 2008;74:286–93.

    Article  CAS  Google Scholar 

  112. Vanwonterghem I, Jensen PD, Rabaey K, Tyson GW. Genome-centric resolution of microbial diversity, metabolism and interactions in anaerobic digestion: Genome-centric resolution through deep metagenomics. Environ Microbiol. 2016;18:3144–58.

    Article  CAS  Google Scholar 

  113. Kumar R, Acharya V, Mukhia S, Singh D, Kumar S. Complete genome sequence of Pseudomonas frederiksbergensis ERDD5:01 revealed genetic bases for survivability at high altitude ecosystem and bioprospection potential. Genomics. 2019;111:492–9.

    Article  CAS  Google Scholar 

  114. Hattori K, Takahashi Y, Guillot S, Johanson B. Occurrence of arsenic (V) in forearc mantle serpentinites based on X-ray absorption spectroscopy study. Geochim Cosmochim Acta. 2005;69:5585–96.

    Article  CAS  Google Scholar 

  115. Mueller B. Preliminary trace element analysis of arsenic in Nepalese groundwater may pinpoint its origin. Environ. Earth Sci. 2018;77:35 s12665-017-7154-z.

    Article  Google Scholar 

  116. Imran Baloch M, Akhtar M, Khan K, Khalid A, Mehmood A, Rukh S, et al. Total and extractable soil selenium contents variation within and across the parent materials. J Biodiversity Environ Sci. 2016;9:175–86.

    Google Scholar 

  117. Jiao NZ, Cai RH, Zheng Q, Tang K, Liu JH, Jiao FLE, et al. Unveiling the enigma of refractory carbon in the ocean. Natl Sci Rev. 2018;5:459–63.

    Article  CAS  Google Scholar 

Download references


We thank the Second Tibetan Plateau Scientific Expedition and Research Program (STEP) 2021QZKK0100 for sampling from the Mount Everest, and thank the captain, crew, and scientists of the R/V “Tansuo-01” during the TS-09 cruise for sampling from the Mariana Trench. We acknowledge the Center for High Performance Computing at Shanghai Jiao Tong University for computing resources.


This study was financially supported by the following funding: the Natural Science Foundation of China (grant numbers 41921006, 42106087 and U21A20176), the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University (project number SL2021PT103), and New Faculty Startup Program of Shanghai Jiao Tong University (project number 21X010500755).

Author information

Authors and Affiliations



Conceptualization: YL, XX, WZ; research design: WZ, XX, YL, ZZ, AH, JW; sample preparation: YL, XX, HJ, KL; metagenomic analysis: WZ, AH, ZZ, JW, XX, YL. Figure preparation: WZ, AH, ZZ, JW. WZ, YL, MJ, XX, ZZ, and AH wrote the manuscript. All authors reviewed the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xiang Xiao or Weishu Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Results and figures.

Additional file 2.

Dataset Table S1.

Additional file 3.

Dataset Table S2.

Additional file 4.


Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Zhang, Z., Ji, M. et al. Comparison of prokaryotes between Mount Everest and the Mariana Trench. Microbiome 10, 215 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: