Skip to main content

Insights into the prokaryotic communities of the abyssal-hadal benthic-boundary layer of the Kuril Kamchatka Trench



The Kuril–Kamchatka Trench (maximum depth 9604 m), located in the NW Pacific Ocean, is among the top seven deepest hadal trenches. The work aimed to investigate the unexplored abyssal-hadal prokaryotic communities of this fascinating, but underrated environment.


As for the bacterial communities, we found that Proteobacteria (56.1–74.5%), Bacteroidetes (6.5–19.1%), and Actinobacteria (0.9–16.1%) were the most represented bacterial phyla over all samples. Thaumarchaeota (52.9–91.1%) was the most abundant phylum in the archaeal communities. The archaeal diversity was highly represented by the ammonia-oxidizing Nitrosopumilus, and the potential hydrocarbon-degrading bacteria Acinetobacter, Zhongshania, and Colwellia were the main bacterial genera. The α-diversity analysis evidenced that both prokaryotic communities were characterized by low evenness, as indicated by the high Gini index values (> 0.9). The β-diversity analysis (Redundancy Analysis) indicated that, as expected, the depth significantly affected the structure of the prokaryotic communities. The co-occurrence network revealed seven prokaryotic groups that covaried across the abyssal-hadal zone of the Kuril–Kamchatka Trench. Among them, the main group included the most abundant archaeal and bacterial OTUs (Nitrosopumilus OTU A2 and OTU A1; Acinetobacter OTU B1), which were ubiquitous across the trench.


This manuscript represents the first attempt to characterize the prokaryotic communities of the KKT abyssal-hadal zone. Our results reveal that the most abundant prokaryotes harbored by the abyssal-hadal zone of Kuril–Kamchatka Trench were chemolithotrophic archaea and heterotrophic bacteria, which did not show a distinctive pattern distribution according to depth. In particular, Acinetobacter, Zhongshania, and Colwellia (potential hydrocarbon degraders) were the main bacterial genera, and Nitrosopumilus (ammonia oxidizer) was the dominant representative of the archaeal diversity.


The oceans represent the greatest ecosystem of our planet providing a plethora of different habitats, hosting a huge diversity of organisms adapted to different environmental conditions [1]. Based on depth, the vertical profile of the oceans is divided into five zones: epipelagic (0–200 m), mesopelagic (200–1000 m), bathypelagic (1000–4000 m), abyssopelagic (4000–6000 m) and hadopelagic (> 6000 m) [1, 2]. These zones show quite diverse physicochemical features (such as salinity, temperature, light, pressure and nutrient concentration) that create distinct conditions affecting microbial communities [3]. Therefore, the characterization of microbial communities across the oceanic zones is important to understand their vertical distribution and ecological role in marine processes, including biogeochemical cycles.

The hadal zones, the deepest part of the oceans composed almost exclusively of trenches, are among the least studied habitats on the planet for their microbial assemblages. Various studies have regarded the isolation of piezophilic and hyperpiezophilic microorganisms from trenches worldwide [4,5,6,7,8]. These include species of archaea, bacteria and fungi including some recently described novel taxa: Muricauda abyssi, Paraoceanicella profunda, Abyssibius alkaniclasticus [9,10,11] As for the profiling of prokaryotic communities, the investigation has been principally focused on the Mariana Trench [3, 12,13,14,15,16,17,18,19,20,21,22,23,24,25,26] and to a lesser extent on other trenches (e.g., Puerto Rico, Japan, Kermadec, Izu-Ogasawara, Atacama, New Britain and Yap trenches) [16, 17, 27,28,29,30,31,32,33,34]. If compared to the literature available for the trench fauna, works that have reported the prokaryote intra-trench variation along the water column are scarce [14, 15, 17, 20,21,22]. Still fewer are the works addressing the inter-trench variability of the prokaryotic diversity; however, most of the studies analyzed the vertical/spatial intra-trench variation of prokaryotes in the Mariana Trench [e.g., 16,17,19,23,24,33]. Meanwhile, the prokaryotic variability within the other trenches or inter-trenches remains poorly investigated. Differences in prokaryotic diversity may exist at different depths within the same trench, but also among different trenches, which are geographically separated habitats with distinctive environmental conditions [4, 29, 35]. Thus, more studies investigating the intra- and inter-trench variation of the prokaryotic communities are needed to extend the knowledge of the community variability of these unique extreme habitats.

The Kuril–Kamchatka Trench (KKT) or Kuril Trench is among the top seven deepest (> 9000 m) hadal trenches [36, 37]. It is a linear (slightly convex) depression of ca. 2089 km in length (surface area of 130,985 km2 and maximum depth of 9604 m) located in the northwest Pacific Ocean, which joins the Japan Trench (creating a hadal continuum) in the south-west and meets the Aleutian Trench in the north-east [36, 38]. The KKT region is affected by various oceanic circulations having different and sometimes opposite directions, resulting also in interactions between cold and warm waters; the trench upper waters are mainly affected by the Oyashio current, whereas two bottom currents affect the deep waters [39, 40]. Actually, water circulation involves only the layers above the hadal region of the trench, whereas in the deeper V-shaped depression no current has been recorded; the Lower Circumpolar Deep Water is the deepest current (~ 4000 m depth) reaching the KKT region [41, 42]. The trench, together with the interconnected Kuril Basin is placed in one of the most productive oceanic areas in the world [43].

Although the KKT region has been investigated since 1949 (by the soviet Vitjaz expeditions; [44]), the studies were principally focused on characterizing its oceanographic, chemical and physical features and on investigating the productivity, plankton and benthic fauna distribution [45,46,47]. Studies carried out on the KKT microbial diversity are definitely scant. One work has been focused on bacteria from hadal bivalve gills [48], and, to the best of our knowledge, only two works (involving sampling in the KKT region or neighboring areas) have investigated the prokaryotic diversity of the site [49, 50]. However, these studies limited the investigation to the upper layers (from the surface to the bathyal zone) considering a maximum depth of 3000 m. Other studies evidenced that many macrofaunal species of the trench and nearby areas floors host chemosynthetic endosymbiotic, indicating that the bottom communities are sustained by the bacterial chemosynthetic organic carbon production [43, 51]: this seems to be a common feature in bottom oceanic fauna [52].

The Aim of this work was an in-depth investigation of the unexplored prokaryotic communities of the abyssal-hadal zone (depth range 5146–9540 m) of the Kuril-Kamchatka Trench. The study was carried out on the benthic-boundary layer waters sampled in various sites of the trench and of the adjacent abyssal plan.

Materials and methods

Sample collection and characterization

The eighteen water samples used for the survey were collected in various abyssal and hadal areas within the Kuril-Kamchatka Trench region (Northwest Pacific Ocean) between August 16 and September 26, 2016, aboard the deep-ocean Research Vessel (RV) Sonne, during the German-Russian expedition KuramBio II (Kurile-Kamchatka Biodiversity Studies II) (Fig. 1).

Fig. 1
figure 1

Kuril-Kamchatka Trench map with indication of the sites (red points) where samples were collected. The map was generated using Google Earth Pro version 7.3.1. and graphically edited

The sampling involved different sites within the trench and in the nearest zones of the adjacent abyssal plain (depth range 5146–9540 m) (Table 1). The samples were representative of the benthic boundary layer (water layers that are directly adjacent to the ocean bottom and can contain particles derived from the resuspension of surface sediments; [15]). For each depth, one sample was taken from the water overlaying (~ 30 cm) one of the cores collected by the OKTOPUS® multicorer (GmbH, Kiel, Germany) equipped with 12 transparent acryl-glass-cylinders (length 60 cm, inner diameter 9.5 cm).

Table 1 Sampling site coordinates and related depths

For each sample, 1 L of seawater was filtered on sterile membranes, which were maintained frozen until DNA extraction (that was carried out on land).

DNA extraction, amplicon libraries and sequencing

Total DNA was extracted from filters using the ZR Fungal/Bacterial DNA MiniPrep kit (Zymo Research Corp., Irvine, CA, USA) according to the manufacturer’s instructions.

The multiplexed amplicon libraries for Bacteria and Archaea (which were prepared and sequenced separately) were obtained using a dual PCR amplification protocol.

The first PCR was performed using primers modified by adding a 6-bp barcodes to allow the parallel processing of multiple samples, as previously reported [53, 54]. For both amplicon libraries, the first PCR was carried out in 3 × 80 µL volume reactions using the Phusion® High-Fidelity DNA Polymerase (NEB Inc., Ipswich, MA, USA) and 1 µM of each primer. For the bacterial community analysis, the V5–V6 hypervariable regions of the 16 S rRNA gene were targeted using the primer set 783 F/1046R [55, 56]. The PCR cycling conditions were: 94 °C for 4 min; 28 cycles at 94 °C for 50 s, 47 °C for 30 s, and 72 °C for 30 s; and a final extension at 72 °C for 5 min. For the archaeal community analysis, the V3 hypervariable region of the 16 S rRNA gene was targeted using the primer set Arch349F/571R [57, 58]. The PCR cycling conditions were: 98 °C for 30 s; 10 cycles at 98 °C for 10 s, 68 °C for 20 s, and 72 °C for 15 s; 35 cycles at 98 °C for 10 s, 58 °C for 15 s, and 72 °C for 15 s; and a final extension at 72 °C for 2 min. The amplicons were purified by Wizard SV Gel and PCR Clean-up System (Promega Corporation, Madison, WI, USA) and quantified using Qubit 2.0. After the purification, nine amplicon samples (identifiable by a unique barcode pair combination) were pooled to obtain a single library.

The second step of the library preparation, consisting of a second PCR for the addition of the standard Nextera indexes (Illumina, Inc., San Diego, CA, USA), the DNA quality evaluation and quantification, and the subsequent amplicon library sequencing were carried out at Nuova Genetica Italiana SRL (Monza-Brianza, Italy). The amplicon libraries were sequenced by Illumina MiSeq (Illumina Inc., San Diego, CA, USA) using a 2 × 250 bp paired-end protocol.

Sequence processing and data analysis

Sequence processing and data analysis were carried out as previously reported [59, 60]. The raw reads obtained from the sequencing were demultiplexed according to the indices and internal barcodes, and processed using the UPARSE pipeline [61]. Forward and reverse reads were merged with perfect overlapping and quality filtered using the default parameters. Suspected chimeras and singletons (unique sequences on the whole data set) were removed. Operational Taxonomic Units (OTUs) were defined on the whole data set clustering the sequences at 97% of similarity, and a representative sequence was defined for each cluster. The OTU abundances per sample were estimated by mapping the sequences to the OTU representative sequences (at 97% of similarity). The taxonomic analysis was performed using RDP classifier [62], applying the 50% confidence cut-off as suggested for short sequences [63].

The sequence data have been submitted to the GenBank database under the accession numbers PRJNA799316 and PRJNA799327.

Statistical methods

Sequence processing and data analysis were carried out as previously reported [64]. The coverage of each sample was estimated by the Good’s method [65]. For both Bacteria and Archaea libraries, rarefaction curves were calculated and plotted to visualize OTU richness as a function of increasing sequencing depth. The Shannon [66] and Gini inequality [67] indices (calculated on rarefied data) were used to describe the α-diversity. The analyses of α-diversity indices were performed using linear models after checking by visual inspection of residual plots (details not shown) that no relevant deviation from model assumption had occurred.

β-diversity analyses were performed on the non-rarefied samples [68] using for each OTU the number of the relating sequences retrieved in each sample as estimation of its abundance in the sample. OTU abundances per sample were transformed using the Hellinger method [69]. The Redundancy Analysis (RDA) was performed to investigate the variation of community structure in relation to the depth (fundamental parameter for the trench environments); the analysis significance was assessed through 9999 permutations. For both datasets (Bacteria and Archaea), the clustering analysis (complete linkage) was performed based on Hellinger distance data. The abundance variation of the most abundant genera (having at least 20,000 (bacterial taxa) and 9,000 (archaeal taxa) sequences in the whole dataset) in relation to the depth was investigated by Generalized Linear Models (GLMs).

The GLMs were performed assuming a Poisson distribution and corrected for overdispersion. We accounted for multiple statistical tests by correcting GLM significance (P-value) according to the False Discovery Rate (FDR) procedure [70]. All the statistical analyses were performed in R environment (R 3.4.2, R Core Team, 2018) using VEGAN, BIODIVERSITYR, MULTTEST, and MULTCOMP packages.

The Functional Annotation of Prokaryotic Taxa software program (FAPROTAX, Version 1.2.6) [71] was used to predict the potential functions of the KKT abyssal-hadal communities, based on the functional annotation information of the cultured prokaryotes.

The network analysis was used to investigate the co-occurrence patterns among bacterial and archaeal OTUs which mainly contribute to the community structure. To reduce the networking complexity, the analysis was carried out only on OTUs with relative abundance (Ra) ≥ 1% at least in one sample. The Spearman correlation matrix (Spearman’s correlation coefficient (ρ) ≥ |0.8|; q values < 0.01 after Bonferroni correction) was calculated on the compositional dataset using the statistical software Systat 8.0 (Systat Software Inc., Point Richmond, CA, USA). The Visone 2.21 free software ( was used for the network exploration and visualization [72].

Results and discussions

This work was carried out within the German-Russian expedition KuramBio II which was specifically planned to increase the knowledge about this oceanic ecosystem, started with a previous oceanic campaign (Kurambio). The main aim of both KURAMBIO expeditions was to study the animal species composition and distribution at the trench bottom, and to evaluate the faunal relationships in the studied deep-sea ecosystems [43]. Samples from KKT, which is practically unknown under the microbiological point of view, were selected for this study due to extreme environmental conditions of the trench and, in particular, its high hydrostatic pressure. During the oceanographic expedition, the environmental parameters related to the benthic-boundary layer samples were not registered due to both technical and logistical reasons.

Bacterial community profiling

After suspected chimera (n = 12,762) and singleton (n = 635) removal, a total of 1,205,906 bacterial mapping sequences were obtained; the Good’s index values (97.9–99.7) indicated a good coverage for the libraries. Overall, 4628 OTUs were defined through the clustering sequence process.

The bacterial OTUs were assigned to 27 phyla, 58 classes, 108 orders, 234 families and 453 genera. Unassigned OTUs ranges were 2.4–9.1%, 5.6–14.1%, 24.9–41.7%, 33.9–50.7%, and 44.9–68.5% at phylum, class, order, family, and genus level, respectively.

At the phylum level, the major bacterial taxa (Ra ≥ 1% in at least one sample) were Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Planctomycetes, Acidobacteria and Marinimicrobia (Fig. 2A). Proteobacteria, Bacteroidetes and Actinobacteria were the most represented taxa over all the KKT samples, and except for Actinobacteria in KKT13 (6065 m depth), they were always recorded with Ra ≥ 1%. Proteobacteria (56.1–74.5%; mean abundance of 66.42%) dominated the abyssal-hadal KKT bacterial communities, followed by Bacteroidetes (6.5–19.1%; mean abundance of 11.69%) and Actinobacteria (0.9–16.1%; mean abundance of 8.11%). Firmicutes (0.3–9.6%) had Ra ≥ 1% in most of the samples (15 out of 18 samples); its mean abundance was 2.9% and it showed the highest abundances (≥ 8.6%) in KKT9 (~ 5146 m depth) and KKT26 (~ 9540 m depth). Planctomycetes (0.8–2.3%) and Acidobacteria (0.3–2.1%) showed mean abundance of 1.5 and 1% and they were recorded among major taxa in 78% and 50% of samples. Marinimicrobia was detected usually as minor taxon, except in KKT14 and KKT26, sampled at ~ 6516 and ~ 9540 m, respectively.

Fig. 2
figure 2

Bacterial composition of the KKT samples at the phylum (A) and genus (B) level. The bar plots show the distribution of the major taxa (Ra ≥ 1% in at least one sample); taxa with Ra < 1% were gathered in “Others”

As expected, the KKT bacterial communities of the abyssal-hadal zone studied in our work differed in both phyla presence and abundance from those of the shallower layers reported by Li et al. [50]. These authors carried out a sampling campaign with stations located around the Sea of Okhotsk and North Western Subarctic Pacific Ocean, including the KKT area affected by the East Kamchatka Current. Samples were mainly taken from the surface and two of them were collected from the mesopelagic (1000 m) and bathypelagic (3000 m) waters. Meso- and bathy-pelagic waters differ from those of the abyssal and hadal ones; they are more unstable and influenced by the currents.

The sample taken at 1000 m, that is generally considered the beginning of piezosphere, should represent a depth where piezophilic microorganisms are considered to acquire advantages over the piezosensitive. Here the waters, should be characterized mainly by the presence of piezotolerant microorganisms than piezophiles, which have their optimal growth conditions at a depth of around 4000–6000 m. The sample taken at 3000 m should not contain piezosensitive micoorganisms but only piezotolerant and piezophilic ones. In any case, in both samples, no extreme-piezophiles should be present [73, 74]. Moreover, a certain taxon could be found at different depths but it could include different species or even strains showing different adaptations to the hydrostatic pressure. For the mentioned reasons the comparison between the present work and that of Li et al. [50] could be considered useful for an intra-trench prokaryotic biodiversity evaluation. Also these authors observed a preponderance of Proteobacteria, which seems to be the main bacterial phylum at almost all depths being its Ra always higher than 40%. In their surface waters, there was also a high presence of Bacteroidetes and Cyanobacteria and sometimes of Verrucomicrobia. In the deeper waters, the biodiversity at the phylum level was higher than that of the surface waters and characterized by the presence also of Planctomycetes, Marinimicrobia, Actinobacteria, and Acidobacteria.

Since our water samples are representative of the benthic boundary layer (presenting also microorganisms from the sediment surface resuspension; [15]) a useful comparison could be done with the prokaryotic communities of both waters and sediments of other trenches (e.g., Mariana Trench, Japan Trench, Puerto Rico Trench, Yap Trench).

Various studies agreed with the general predominance of Proteobacteria in the abyssal-hadal zones [e.g.,14,16,18,22,27,29,31], whereas others reported dominance of different taxa. For example, in the work of Nunoura et al. [12], carried out in the Challenger Deep in the Mariana Trench, Proteobacteria dominated in the bathyal zone (1000–2000 m) and in the hadal zone at 10,000 m, while Bacteroidetes dominated in the abyssal-hadal zone from 6000 to 9000 m. In the same area, Nunoura et al. [13] observed a very strong abundance of Protobacteria in hadal waters at 10,257 m, whereas Chloroflexi and into a second extent Bacteroidetes, Planctomycetes, Marinimicrobia and Gemmatimonadetes were the most abundant taxa in the sediments. In a later work [21], characterizing the Challenger Deep prokaryotic communities, in most of the samples taken from the hadal waters, Marinimicrobia was the most abundant bacterial phylum. Dominance of Marinimicrobia was found also in hadal waters and sediments from other sites [12, 13, 20, 30, 75]. In our case, except for KKT14 (~ 7137 m) and KKT26 (~ 9540 m), this phylum was always a rare taxon. Overall, an inter-trench diversification of the bacterial community assemblages is evident already at the phylum level.

However, a possible comparison of data related to different trenches obtained in different periods (years and/or months), would be more meaningful if referred to samples taken in stable abyssal-hadal zones rather than with those obtained from the upper less stable layers subject to various perturbations such as those caused by oceanic currents.

At the genus level, a very high abundance of unassigned OTUs (from 44.9 to 68.5%) occurred (Fig. 2B). As underlined by others [76, 77], a high proportion of unassigned microbial sequences in deep seas suggests that the majority of the microbial diversity remains largely unknown. However, the sequence length does not allow to confidently state whether the unassigned sequences could be attributed to known genera or to putative new taxa. To overcome this problem the complete 16S rRNA gene could be sequenced in a subsequent study.

Even if the subsequent observations are based on relative abundance data calculated from the sole OTUs assigned at genus level, it is still possible to make some general considerations.

The bacterial OTUs were assigned to 453 genera; among them, 26 taxa showed Ra ≥ 1% in at least one sample. The most represented genus was Acinetobacter, followed by Zhongshania and Colwellia.

Acinetobacter was the sole genus representing a major taxon in all samples, with no significant abundance variation in relation to the depth, as suggested by the GLM analysis (PFDR > 0.05). It was recorded with abundances spanning from 5.4% (in KKT19; ~9499 m depth) to 27.3% (in KKT9; ~5147 m depth). Moreover, its highest abundances (≥ 20%) were recorded both in the upper layers of the abyssal zone (~ 5147 and ~ 5217 m depth) and in the deep layers of the hadal zone (~ 8221 and ~ 9540 m depth). Acinetobacter has been found in a wide array of environments and includes various species having different ecological roles, from pathogeny to biodegradation of xenobiotics [78]. It has been revealed in marine sediments and waters, from the surface to the deepest zones, including some oceanic trenches [31, 79,80,81,82,83,84,85,86,87,88]. The metabolic features and the potential ecological functions of some deep-sea members of Acinetobacter indicated that they can be involved in sulphur and metal oxidation, and hydrocarbon degradation [80, 84, 88,89,90].

Zhongshania and Colwellia were detected in all samples, but they showed Ra < 1% in six (KKT9, KKT10, KKT19 KKT21, KKT23, and KKT26) and two (KKT9 and KKT19) samples, respectively. Also for these genera, a significant abundance variation according to the depth was not evidenced (PFDR > 0.05). Zongshania ranged from 0.1 to 11.7%, and it showed the highest abundances (> 9%) at ~ 6065 m and ~ 8255 m. Colwellia showed abundances in the range 0.4–5.5% with the highest abundance at ~ 6065 m. Zhongshania, usually detected as major taxon in the KKT abyssal and hadal samples, has various psychrophilic and psychrotolerant members isolated from diverse marine environments (from Antarctic coastal waters to tidal sediments of the Yellow Sea), but only up to 500 m depth [91,92,93,94]. No information concerning its presence in the deepest oceanic zones is available and there is scarce knowledge about its behavior in relation to increasing hydrostatic pressure. Only one work tested Zhongshania strains (isolated from coastal seawater samples) at increasing hydrostatic pressures, evidencing their piezotolerant phenotype [95]. Hence, as far as we know, our work is the first reporting the presence of these bacteria in the deepest oceanic zones (5146–9540 m depth) possibly confirming the existence of piezophilic/piezotolerant members. It is interesting to note that various members of the genus Zhongshania are hydrocarbon-degrading bacteria [93, 95,96,97]. In addition, a sequential cooperation among oil-degrading strains of Zhongshania and Colwellia was reported: Zhongshania plays a fundamental role in the initial n-alkanes biotransformation that is completed by Colwellia, through different metabolic pathways including also the beta‑oxidation [96].

Therefore, it is possible to speculate that the most represented bacteria found in the KKT abyssal-hadal zone are potentially hydrocarbon-degraders. In addition, the presence of well-known obligate hydrocarbonoclastic bacteria (Alcanivorax, Oleibacter and Oleispira; [98]) in some KKT samples, although as minor taxa (Fig. 3), might reinforce this hypothesis, indicating also that hydrocarbons could represent a significant carbon source in these deep environments.

Fig. 3
figure 3

Occurrence of Acinetobacter, Zhongshania, Colwellia and genera of obligate hydrocarbonoclastic bacteria (Oleispira, Oleibacter and Alcanivorax) across the KKT abyssal and hadal samples

Presence of hydrocarbons at the trench bottom could be due both to xenobiotic (e.g., oil spill falls) and biotic (from heterotrophic and autotrophic organisms) sources revealed in sinking particles [18]. Thus, the occurrence of the hydrocarbonoclastic bacteria in the KKT samples is consistent with what found for other trench environments: the deepest zones support the growth of these heterotrophic microorganisms that locally contribute to hydrocarbon recycle at the bottom of the ocean [3, 99].

A list of the bacterial genera found in the abyssal-hadal zones of KKT and other main trenches was supplied to provide a more detailed comparison (Additional file 2: Table S1).

Archaeal community profiling

After suspected chimera (n = 3613) and singleton (n = 977) removal, a total of 1,226,044 archaeal mapping sequences were obtained; the Good’s index values (99.4–99.8) indicated a good coverage for the libraries. Overall, 3362 OTUs were defined through the clustering sequence process.

Since for some archaeal genera the lineage is currently not fully defined, we report only the summary of OTU assignation obtained at phylum and genus levels. The archaeal OTUs were assigned to 5 phyla and 21 genera. Unassigned OTUs ranged from 0.3 to 3.6% and from 6 to 33.3% at phylum and genus level, respectively. The archaeal phyla detected were Thaumarchaeota, Woesearchaeota, Euryarchaeota, Crenarchaeota, and Diapherotrites (Fig. 4A). Among them, Diapherotrites was the sole phylum representing always a minor taxon (it was present in 13 out of 18 samples with Ra in the range 0.001-0.2%). Crenarchaeota was detected as a major genus only in two samples (1.7 and 2.6% in KKT19 and KKT13, respectively). Members of this phylum have been detected in various oceanic environments, often representing an abundant archaeal fraction in oxygenated deep waters and surface sediments [100].

Fig. 4
figure 4

Archaeal composition of the KKT samples at the phylum (A) and genus (B) level. The bar plots show the distribution of the major taxa (Ra ≥ 1% in at least one sample); taxa with Ra < 1% were gathered in “Others”

Thaumarchaeota (52.9–91.1%), Woesearchaeota (3.3–31.1%), and Euryarchaeota (0.6–39.1%) were the most represented taxa over all samples; except for Euryarchaeota in KKT23 (~ 6518 m depth), they were always detected with Ra ≥ 1%.

The archaeal community assemblages revealed for the KKT samples are consistent with those recognized by Jing et al. [15] in the benthic boundary layer of the Mariana Trench. As already observed for the bacteria, also the archaeal communities of the KKT abyssal-hadal zone differ from those of the shallower layers in phyla presence/abundance [50]. However, Li et al. [50] outcomes confirmed Thaumarchaeota and Euryarchaeota as the main archaeal phyla in the upper layers (1000 and 3000 m depths) as well.

The current study evidenced that Thaumarchaeota dominated the archaeal communities of the KKT abyssal-hadal zone, showing a notable preponderance (Ra ≥ 80%) in one third of the samples. Other studies revealed a dominance of this taxon in the archaeal communities of both waters and sediments across various deep oceanic habitats, where they are followed by Woesearchaeota (mainly in sediments) and Euryarchaeota (mainly in waters) [13, 14, 16, 21,22,23,24, 28, 29, 31, 32].

Marine Thaumarchaeota are well-known chemolithoautotrophs (spanning from mesopelagic to hadal environments) involved in dark primary production, with key roles in nitrogen and carbon biogeochemical cycles [101].

Marine Euryarchaeota has been divided in three major groups (MG-II, MG-III, and MG-IV). MG-II members represent the dominant fraction of the archaeal communities within the euphotic zone, but they have also been found in deep-sea waters. MG-III and MG-IV are considered to be rare components of deep-sea communities [102]. Little is known about the ecological partitioning and physiology of these groups; in addition, the scarce knowledge of their metabolic competencies means their impact on biogeochemical cycles still speculative. However, MGII and MGIII have been suggested to be involved in the degradation of high-molecular-weight compounds, therefore affecting primary the carbon biogeochemical cycle [103]. Woesearchaeota has been found in different environments, including marine waters/sediments and hydrothermal vents [104]; however, its biodiversity, distribution and metabolisms remain largely unknown [105]. Nevertheless, it has been suggested that most of Woesearchaeota might have a mainly anaerobic or fermentation-based lifestyle, also under potential syntrophic/mutualistic partnership with other microbial groups [104, 105]. It has been demonstrated that Woesearchaeota tend to co-occur in a consortium with anaerobic methanogens of Methanomicrobia and Methanobacteria, and show a possible role in the hydrogen cycle [105]. Also our survey revealed that Woesearchaeota and Methanomicrobia were concurrently detected at all depths (Fig. 5); thus, a potential syntrophic partnership between these two taxa could be reasonably suggested.

The archaeal genera detected as major taxon in at least one KKT sample were Nitrosopumilus, Halococcus, Woesearchaeota Incertae Sedis AR16, Woesearchaeota Incertae Sedis AR20, Halorubrum, Woesearchaeota Incertae Sedis AR15, Salinigranum and Halobellus (Fig. 4B).

Fig. 5
figure 5

Occurrence of Woesearchaeota and classes of methanogens (Methanomicrobia and Methanobacteria) across the KKT abyssal and hadal samples

Among them, the most represented taxa were Nitrosopumilus (52.9-91-1%), Halococcus (0.1–36%), and Woesearchaeota Incertae Sedis AR16 (0.7-9%). All the most abundant genera did not show significant differences in abundances in relation to the depth (PFDR > 0.05), indicating a general ability to cope with different hydrostatic pressure conditions. Woesearchaeota Incertae Sedis AR16 was among the major taxa in 16 samples out of 18, it had mean abundances of 3.3% with the highest abundance (9%) in KKT25 (~ 9539 m depth). Halococcus was recorded with Ra ≥ 1% in eight samples; it showed high abundance (30.4 and 36%) in KKT9 and KKT11 collected at ~ 5147 m and ~ 8224 m depth, respectively.

Hence, the archaeal communities of the KKT abyssal-hadal zone were dominated by Nitrosopumilus, which accounted at least for 99.98% of Thaumarchaeota that was the most abundant archaeal phylum. Therefore, Nitrosopumilus represented the main fraction of the whole KKT archaeal diversity.

Nitrosopumilus is a globally abundant marine genus with high biogeochemical significance playing a critical function in global element cycling particularly for carbon (due to their chemolithoautotrophy) and nitrogen (for their role in ammonia oxidation) [3, 14, 16, 17, 19,20,21, 23, 34, 106]. These chemolithoautotrophs and ammonia-oxidizing archaea are considered key primary producers principally when photoautotrophy is not possible [107, 108]. Moreover, the ammonia-oxidizing archaea are generally found in co-presence with nitrite-oxidizing bacteria in most ecosystems [109, 110]. In the oceans, ammonia is mainly oxidized by ammonia oxidizing archaea (among them Nitrosopumilus), nitrite is further oxidized to nitrate by nitrite oxidizing bacteria, primarily belonging to the phylum Nitrospinae and to a lesser extent to the genera Nitrococcus and Nitrospira [109]. The presence of nitrite oxidizing bacteria of Nitrospina (phylum Nitrospinae) and Nitrospira was revealed (although with Ra ≤ 1%) in the KKT samples (Fig. 6), suggesting a potential partnership with the ammonia oxidizing archaeon Nitrosopumilus.

Fig. 6
figure 6

Occurrence of Nitrosopumilus and genera of nitrite oxidizing bacteria (Nitrospira and Nitrospina) across the KKT abyssal and hadal samples

A list of the archaeal genera found in the abyssal-hadal zones of KKT and other main trenches was supplied to provide a more detailed comparison (Additional file 2: Table S2).

Alpha- and Beta-diversity of the bacterial and archaeal communities

Rarefaction curves were separately evaluated for Bacteria and Archaea datasets (Additional file 1: Fig. S1) to perform dataset rarefaction to the minimum library depth.

The α-diversity indices of the bacterial and archaeal communities of the KKT abyssal and hadal samples are reported in Table 2. Summarizing, the Number of OTUs ranged from 716 to 2163 and from 243 to 1031 for the bacterial and archaeal libraries, respectively. The Shannon diversity index ranged from 2.686 to 5.818 and from 3.155 to 4.764 for Bacteria and Archaea, respectively. The Gini inequality index ranged from 0.927 to 0.992 and from 0.959 to 0.993 for Bacteria and Archaea, respectively. The high Gini values obtained for both prokaryotic communities evidenced that they were characterized by low evenness, and in turn by the presence of few dominant OTUs.

Table 2 Alpha-diversity indices of the KKT bacterial and archaeal amplicon libraries

The depth and the hydrostatic pressure are key interrelated parameters, which are intrinsically fundamental in deep oceanic environments. To understand if differences in hydrostatic pressure exerted significant effects on KKT α-diversity, the GLMs were run to analyze whether the variation of indices according to the depth was significant. The analysis evidenced that for both bacterial and archaeal communities the α-diversity indices did not show significant variation (p > 0.05) in relation to this key-parameter.

As expected, the RDA analysis revealed that the hydrostatic pressure clearly affected the whole structure of the prokaryotic communities. The KKT abyssal-hadal prokaryotic communities significantly differed (Bacteria: z = 3.373, P = 0.001, adjusted-R2 = 0.122; Archaea: z = 2.693, P = 0.001, adjusted-R2 = 0.091) according to the depth, which explained the 5.42% and 4.29% of data variability for the bacterial and archaeal community, respectively (Fig. 7).

Fig. 7
figure 7

Redundancy analysis (RDA) ordination biplot showing the relationship between the KKT samples (based on bacterial (A) and archaeal (B) OTU profiles) and the depth

Although no data are available for the considered ecosystem, it is possible to hypothesize that other variables in addition to the depth, such as nutrients, can explain the structuring of the KKT abyssal-hadal prokaryotic communities. Instead, as generally known, since the trench abyssal-hadal environments are quite stables, with temperature, light, and salinity that are always constant, it is reasonable to think that these environmental parameters do not represent key drivers for the community structure.

For Bacteria, the Cluster Analysis (based on Hellinger distance) identified five groups (Additional file1: Fig. S2). Based on their bacterial communities, unexpectedly, samples KKT25 (9538.6 m) and KKT26 (9540.2 m) grouped together with the samples KKT9 and KKT10, which were collected at the less deep layers (5146.6-5146.7 m). The samples collected in-between 8224 and 8254.7 m depth (KKT11, KKT12, and KKT22) clustered with those taken at 7134.6-7136.7 m (KKT14 and KKT15). Instead, KKT18 (8221.1 m) grouped with less deep samples (KKT16, KKT17, KKT20, KKT21, KKT23, and KKT24) that were collected at depths 5211.3-6517.5 m. Compared to these clusters, KKT13 (6065 m) and KKT19 (9498.9 m) showed the most divergent bacterial communities.

For Archaea, three main groups were identified through the Cluster Analysis (Additional file1: Fig. S3). Overall, KKT19 showed the most divergent archaeal community. A cluster included samples collected at the deepest layers (8224-9540.2 m), which grouped unexpectedly together with KKT21 (taken at 5217.4 m). Finally, a larger cluster included the samples collected at depths in-between 5146.6 and 8221.1 m.

Bacterial and archaeal ecological function prediction

The functional diversity analysis (FAPROTAX) showed that for both Bacteria and Archaea the predicted ecological functions were mainly associated with the carbon, nitrogen, and sulfur cycles (Figs. 8 and 9).

Fig. 8
figure 8

Ecological functions predicted for the KKT abyssal-hadal bacterial communities

Fig. 9
figure 9

Ecological functions predicted for the KKT abyssal-hadal archaeal communities

Chemoheterotrophy was the dominant metabolic strategy revealed in the bacterial communities and associated with 219 out of the 453 identified genera, including the aerobic chemoheterotrophic Acinetobacter, Colwellia, Shewanella, Pseudomonas, Candidatus Pelagibacter, and Pseudoalteromonas, which were among the major taxa found in KKT samples (Fig. 8; Additional file 2: Table S3; Fig. 2).

A rather high fraction, accounting for 5.4–27.3% of the whole bacterial communities, was characterized by aromatic compound degradative ability, related principally to Acinetobacter and to a lesser extent to Gordonia, Nocardioides, Rhodococcus, and Sporobacter (Fig. 8; Additional file 2: Table S3).

Among nitrogen cycle-related functions, nitrate reduction, nitrate respiration, and nitrogen respiration were the most represented predicted functions (in particular in KKT24, collected at 6515.8 m), largely associated with Shewanella; Achromobacter, Aquabacterium, Azospira, Desulfuromusa, Paracoccus, and Stenotrophomonas were involved to a lesser degree (Fig. 8; Additional file 2: Table S3).

As for the sulfur cycle, respiration of sulfur compounds and in particular sulfate respiration were the most predicted metabolisms, involving just minor genera such as Desulfatiglans, Desulfatitalea, Desulfobacula, and Desulfobulbus (Fig. 8; Additional file 2: Table S3).

Regarding the archaeal communities, the aerobic ammonia oxidation (the first step of the nitrification process) was the most represented ecological function, principally related to Nitrosopumilus and to a lesser extent to Nitrososphaera members (Fig. 9; Additional file 2: Table S4).

Also nitrate reduction (step of the denitrification process) was among the most predicted functions related to the nitrogen cycle, involving Halorubrum, Natronomonas, and Thermodiscus (Fig. 9; Additional file 2: Table S4).

A rather high fraction of the archaeal communities (in particular in KKT9 and KKT11 samples) was characterized by a chemoheterotrophyc metabolism, relating to Halococcus, Methanomassiliicoccus, Natronomonas, and Thermodiscus members (Fig. 9; Additional file 2: Table S4).

Metabolic functions related to the sulfur cycle were less represented and uniquely associated with the genus Thermodiscus (Fig. 9; Additional file 2: Table S4).

Bacterial and archaeal OTU co-occurrence network

The co-occurrence analysis was carried out to understand the interactions among the major bacteria and archaeal OTUs (Ra ≥ 1% in at least one sample) found in the benthic boundary layer of the KKT abyssal-hadal zone. The co-occurrence network (Fig. 10) showed seven groups (A-G) consisting of 2–34 OTUs (54 OTUs in total) and 1–60 edges.

Fig. 10
figure 10

Co-occurrence network between bacterial and archaeal OTUs. Nodes represent the archaeal (circles) and bacterial (squares) OTUs, and edges (grey lines) represent the statistically significant positive correlations between each OTU pair; negative correlations are not shown. The node size is proportional to the mean relative abundance of OTUs in the dataset. The nodes are colored according to the taxonomic assignation (confidence ≥ 0.50)

There is a main group (group A) subdivided in 3 sub-groups linked by a connection node constituted by Nitrosopumilus OTU A24. This group included the principal (most abundant and ubiquitous in the KKT abyssal-hadal communities) archaeal and bacterial OTUs: Nitrosopumilus OTU A2 and OTU A1, (in subgroup A-I and A-II, respectively), Acinetobacter OTU B1 (in A-II). Except for two bacterial OTUs, belonging to Colwelliaceae (B128) and Oceanospirillales (B33), the main subgroup (A-I) showed a strong predominance of Nitrosopumilus OTUs. This subgroup was the highest connected module, with Nitrosopumilus OTU A2 and OTU A16 as the most interconnected nodes. Subgroup A-II consisted of various Nitrosopumilus OTUs connected with two bacterial OTUs. Among them, OTU B2 (class Actinobacteria) represented a node linking a peripheral core of bacterial OTUs. The subgroup A-III was the sole consisting of all bacterial OTUs. The groups B, C and E included only archaeal OTUs; among them, the most interconnected was Group B constituted by six Halococcus OTUs. Group D included three Nitrosopumilus and one Planktomarina OTUs. The groups F and G presented only two bacterial OTUs and connected respectively members of Clostridium with Shewanella, and Gammaprotoeobacteria with Candidatus Pelagibacter.

As expected (based on the main ecological functions predicted through FAPROTAX analysis), overall, the groups comprised ammonia-oxidizing archaea (Nitrosopumilus) and chemoheterotrophic bacteria and archaea (Acinetobacter, having also aromatic compound degradation abilities; Shewanella, involved also in nitrate reduction, nitrate and nitrogen respiration; Candidatus Pelagibacter; Oceanospirillales and Colwelliaceae members; Halococcus) (Additional file 2: Table S3 and Table S4).

It is not possible to find a strict relation between the various groups of the network and depth, being almost all OTUs detected at all depths. However, it is worth noting that, within the various groups/subgroups, most of the OTUs were recorded as major or even highly abundant at certain depths (Additional file 2: Table S5). For example, OTU A2 was always present with Ra > 4% and recorded with the highest abundances particularly in the range 5741-7135 m. All the A-II OTUs showed Ra ≥ 1% at the maximum depth (~ 9540 m), and most of them showed highest abundances in the hadal zone. Group C OTUs showed their highest abundances mainly in the KKT abyssal zone; whereas those of group D were major only in KKT19 (~ 9499 m). Group E was represented only by the archaeal OTUs A19 e A3, which showed their highest abundances in KKT22 (8224 m). The bacterial OTUs B8 and B70 in Group F showed highest abundances in KKT24 and KKT26 (~ 6516 m and ~ 9540 m, respectively). Group G gathers the bacterial OTUs B17 and B48, which were major just in one sample (KKT19, ~ 9499 m) with quite similar abundances.


Summarizing, the present work allowed portraying the bacterial and archaeal communities of the KKT abyssal-hadal zone, which were never investigated before. Although the availability of environmental parameters to relate with the studied prokaryotic diversity would have provided a more informative ecological overview, the work represents an important contribution to the scientific community, especially if taking into account that no further explorations of the area are likely to be organized probably for many years to come. The survey evidenced that this extreme environment harbored communities that differed from those of other marine trenches. The highly uneven KKT bacterial and archaeal communities were mainly represented by heterotrophic and chemolithotrophic genera, respectively, which did not show a distinctive zonation but were detected at all depths. In particular, Acinetobacter, Zhongshania, and Colwellia (potential hydrocarbon degraders) were the main bacterial genera, and Nitrosopumilus (ammonia oxidizer) was the dominant archaeal representative of the KKT deepest areas.

Data Availability

The datasets generated during and/or analyzed during the current study are available in the Sequence Read Archive (SRA; repository (SRA accession numbers: PRJNA799316 and PRJNA799327).



Kuril–Kamchatka Trench


Operational Taxonomic Unit


Redundancy Analysis


Generalized Linear Model


False Discovery Rate


Relative abundance


  1. Kennish MJ. Practical handbook of marine science. Boca Raton: CRC press; 2019.

    Book  Google Scholar 

  2. Speight MR, Henderson PA. Marine ecology: concepts and applications. New York: John Wiley & Sons; 2013.

    Google Scholar 

  3. Xue C, Liu J, Lea-Smith DJ, Rowley G, Lin H, Zheng Y, et al. Insights into the vertical stratification of microbial ecological roles across the deepest seawater column on Earth. Microorganisms. 2020;8(9):1309.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Liu R, Wang L, Wei Y, Fang J. The hadal biosphere: recent insights and new directions. Deep-Sea Res II: Top Stud Oceanogr. 2018;155:11–8.

    Article  Google Scholar 

  5. Nogi Y, Kato C, Horikoshi K. Psychromonas kaikoae sp. nov., a novel from the deepest piezophilic bacterium cold-seep sediments in the Japan Trench. Int J Syst Evol Microbiol. 2002;52:1527–32.

    Article  CAS  PubMed  Google Scholar 

  6. Pathom-Aree W, Nogi Y, Sutcliffe IC, Ward AC, Horikoshi K, Bull AT, et al. Dermacoccus abyssi sp. nov., a piezotolerant actinomycete isolated from the Mariana Trench. Int J Syst Evol Microbiol. 2006;56:1233–7.

    Article  CAS  PubMed  Google Scholar 

  7. Cao Y, Chastain RA, Eloe EA, Nogi Y, Kato C, Bartletta DH. Novel psychropiezophilic oceanospirillales species Profundimonas piezophila gen. nov., sp. nov., isolated from the deep-sea environment of the puerto rico trench. Appl Environ Microbiol. 2014;80:54–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Wang H, Zhang Y, Bartlett DH, Xiao X. Transcriptomic analysis reveals common adaptation mechanisms under different stresses for moderately piezophilic Bacteria. Microb Ecol. 2021;81:617–29.

    Article  CAS  PubMed  Google Scholar 

  9. Wang D, Wu Y, Liu Y, Liu B, Gao Y, Yang Y, et al. Muricauda abyssi sp. nov., a marine bacterium isolated from deep seawater of the Mariana Trench. Int J Syst Evol Microbiol. 2022;72:005615.

    Article  CAS  Google Scholar 

  10. Liu P, Ding W, Lai Q, Liu R, Wei Y, Wang L, et al. Physiological and genomic features of Paraoceanicella profunda gen. nov., sp. nov., a novel piezophile isolated from deep seawater of the Mariana Trench. Microbiologyopen. 2020;9:e966.

    Article  CAS  PubMed  Google Scholar 

  11. Wei TT, Fan XB, Quan ZX. Abyssibius alkaniclasticus gen. nov., sp. nov., a novel member of the family Rhodobacteraceae, isolated from the Mariana Trench. Int J Syst Evol Microbiol. 2023;73:005715.

    Article  CAS  Google Scholar 

  12. 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 Nat Acad Sci. 2015;112(11):E1230–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. 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(2):186–94.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 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  PubMed  PubMed Central  Google Scholar 

  15. Jing H, Zhu W, Liu H, Zheng L, Zhang Y. Particle-attached and free-living archaeal communities in the benthic boundary layer of the Mariana Trench. Front Microbiol. 2018;9:2821.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 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(4):e0195102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Peoples LM, Grammatopoulou E, Pombrol M, Xu X, Osuntokun O, Blanton J, et al. Microbial community diversity within sediments from two geographically separated hadal trenches. Front Microbiol. 2019;10:347.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Liu J, Zheng Y, Lin H, Wang X, Li M, Liu Y, et al. Proliferation of hydrocarbon-degrading microbes at the bottom of the Mariana Trench. Microbiome. 2019;7(1):1–13.

    Article  Google Scholar 

  19. Liu R, Wang Z, Wang L, Li Z, Fang J, Wei X, et al. Bulk and active sediment prokaryotic communities in the Mariana and Mussau Trenches. Front Microbiol. 2020;11:521.

    Article  Google Scholar 

  20. Cui G, Li J, Gao Z, Wang Y. Spatial variations of microbial communities in abyssal and hadal sediments across the Challenger Deep. PeerJ. 2019;7:e6961.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Gao ZM, Huang JM, Cui GJ, Li WL, Li J, Wei ZF, et al. In situ meta-omic insights into the community compositions and ecological roles of hadal microbes in the Mariana Trench. Environ Microbiol. 2019;21(11):4092–108.

    Article  CAS  PubMed  Google Scholar 

  22. Wang Y, Gao ZM, Li J, He LS, Cui GJ, Li WL, et al. Hadal water sampling by in situ microbial filtration and fixation (ISMIFF) apparatus. Deep Sea Res Part I Oceanogr Res. 2019;144:132–7.

    Article  Google Scholar 

  23. Wang Z, Wang L, Liu R, Li Z, Wu J, Wei X, et al. Community structure and activity potentials of archaeal communities in hadal sediments of the Mariana and Mussau trenches. Mar Life Sci Technol. 2022;4(1):150–61.

    Article  CAS  PubMed  Google Scholar 

  24. Hiraoka S, Hirai M, Matsui Y, Makabe A, Minegishi H, Tsuda M, et al. Microbial community and geochemical analyses of trans-trench sediments for understanding the roles of hadal environments. ISME J. 2020;14(3):740–56.

    Article  CAS  PubMed  Google Scholar 

  25. Guo R, Liang Y, Xin Y, Wang L, Mou S, Cao C, et al. Insight into the Pico- and Nano-Phytoplankton Communities in the Deepest Biosphere, the Mariana Trench. Front Microbiol. 2018;9:403251.

    Article  Google Scholar 

  26. Wang W, Sun J, Hao J. Spatial variability of bacterial community compositions in the Mariana Trench. 2022;68:633–42.

  27. Eloe EA, Fadrosh DW, Novotny M, Allen LZ, Kim M, Lombardo M-J, et al. Going deeper: metagenome of a hadopelagic microbial community. PLoS ONE. 2011;6(5):e20388.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Eloe EA, Shulse CN, Fadrosh DW, Williamson SJ, Allen EE, Bartlett DH. Compositional differences in particle-associated and free-living microbial assemblages from an extreme deep-ocean environment. Environ Microbiol Rep. 2011;3(4):449–58.

    Article  PubMed  Google Scholar 

  29. Nunoura T, Hirai M, Yoshida-Takashima Y, Nishizawa M, Kawagucci S, Yokokawa T, et al. Distribution and niche separation of planktonic microbial communities in the water columns from the surface to the hadal waters of the Japan Trench under the eutrophic ocean. Front Microbiol. 2016;7:1261.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Liu R, Wang L, Liu Q, Wang Z, Li Z, Fang J, et al. Depth-resolved distribution of particle-attached and free-living bacterial communities in the water column of the New Britain Trench. Front Microbiol. 2018;9:625.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Zhang X, Xu W, Liu Y, Cai M, Luo Z, Li M. Metagenomics reveals microbial diversity and metabolic potentials of seawater and surface sediment from a hadal biosphere at the Yap Trench. Front Microbiol. 2018;9:2402.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Fu L, Li D, Mi T, Zhao J, Liu C, Sun C, et al. Characteristics of the archaeal and bacterial communities in core sediments from Southern Yap Trench via in situ sampling by the manned submersible Jiaolong. Sci Total Environ. 2020;703:134884.

    Article  CAS  PubMed  Google Scholar 

  33. Schauberger C, Glud RN, Hausmann B, Trouche B, Maignien L, Poulain J, et al. Microbial community structure in hadal sediments: high similarity along trench axes and strong changes along redox gradients. ISME J. 2021;15(12):3455–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Zhao X, Luo H, He S, Yang B, Wei T, Hu Y, et al. Vertical distribution of size-fractionated bacterial communities in the water column of the Atacama Trench. Reg Stud Mar Sci. 2022;55:102470.

    Article  Google Scholar 

  35. Wenzhöfer F, Oguri K, Middelboe M, Turnewitsch R, Toyofuku T, Kitazato H, et al. Benthic carbon mineralization in hadal trenches: Assessment by in situ O2 microprofile measurements. Deep Sea Res Part I. 2016;116:276–86.

    Article  CAS  Google Scholar 

  36. Dreutter S, Steffen M, Arbizu PM, Brandt A. Will the “top five” deepest trenches lose one of their members? Prog Oceanogr. 2020;181:102258.

    Article  Google Scholar 

  37. Stern RJ. Ocean Trenches. Encyclopedia of Geology. Elsevier; 2021. 845–54.

  38. Jamieson AJ, Stewart HA. Hadal zones of the Northwest Pacific Ocean. Prog Oceanogr. 2021;190:102477.

    Article  Google Scholar 

  39. Mitsuzawa K, Holloway G. Characteristics of deep currents along trenches in the northwest Pacific. J Geophys Res Oceans. 1998;103(C6):13085–92.

    Article  Google Scholar 

  40. Qiu B. Kuroshio and Oyashio currents. Ocean currents: a derivative of the Encyclopedia of Ocean Sciences. Elsevier; 2001. 61–72.

  41. Rella SF, Uchida M. A Southern Ocean trigger for Northwest Pacific ventilation during the Holocene? Sci Rep. 2014;4(1):1–11.

    Article  CAS  Google Scholar 

  42. Fuhr M, Laukert G, Yu Y, Nürnberg D, Frank M. Tracing water mass mixing from the equatorial to the North Pacific Ocean with dissolved neodymium isotopes and concentrations. Front Ma Sci. 2021;7:1261.

    Article  Google Scholar 

  43. Kamenev GM, Mordukhovich VV, Alalykina IL, Chernyshev AV, Maiorova AS. Macrofauna and nematode abundance in the abyssal and hadal zones of interconnected deep-sea ecosystems in the Kuril Basin (Sea of Okhotsk) and the Kuril-Kamchatka Trench (Pacific Ocean). Front Mar Sci. 2022;9:812464.

    Article  Google Scholar 

  44. Wolff T. The hadal community, an introduction. Deep Sea Res (1953). 1959;6:95–124.

    Article  Google Scholar 

  45. Bogorov VR. Fauna of the Kurile–Kamchatka Trench and its environment, based on data of the 38th cruise of the R/V “Vityaz. Jerusalem: Israel Program for Scientific Translations; 1972.

    Google Scholar 

  46. Zenkevich LA. Biology of the Seas of the USSR. London: George Allen & Unwin LTD; 1963.

    Book  Google Scholar 

  47. Fischer V, Brandt A. Composition of the abyssal infauna of the Kuril–Kamchatka area (NW Pacific) collected with a box corer. Deep-Sea Res II: Top Stud Oceanogr. 2015;111:26–33.

    Article  Google Scholar 

  48. Krylova EM, Drozdov AL, Mironov AN. Presence of bacteria in gills of hadal bivalve “Vesicomyasergeevi Filatova, 1971. Ruthenica. 2000;10:76–9.

    Google Scholar 

  49. Jing H, Cheung S, Xia X, Suzuki K, Nishioka J, Liu H. Geographic distribution of ammonia-oxidizing archaea along the Kuril Islands in the western Subarctic Pacific. Front Microbiol. 2017;8:1247.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Li Y, Jing H, Xia X, Cheung S, Suzuki K, Liu H. Metagenomic insights into the microbial community and nutrient cycling in the western subarctic Pacific Ocean. Front Microbiol. 2018;9:623.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Karaseva N, Gantsevich M, Obzhirov AI, Shakirov RB, Starovoytov A, Smirnov R, et al. Siboglinids (Annelida. Siboglinidae) as possible indicators of carbohydrates on the case of the Sea of Okhotsk. Dokl Akad Nauk. 2019;486(1):127–30.

    Article  Google Scholar 

  52. Levin LA. Oxygen minimum zone benthos: adaptation and community response to hypoxia. Oceanography and Marine Biology, an Annual Review, Vol. 41. CRC Press; 2003. 1–45.

  53. Gagliano AL, Tagliavia M, D’alessandro W, Franzetti A, Parello F, Quatrini P. So close, so different: geothermal flux shapes divergent soil microbial communities at neighbouring sites. Geobiology. 2016;14(2):150–62.

    Article  CAS  PubMed  Google Scholar 

  54. Gorrasi S, Franzetti A, Ambrosini R, Pittino F, Pasqualetti M, Fenice M. Spatio-temporal variation of the bacterial communities along a salinity gradient within a thalassohaline environment (saline di Tarquinia Salterns, Italy). Molecules. 2021;26(5):1338.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Huber JA, Welch DBM, Morrison HG, Huse SM, Neal PR, Butterfield DA, et al. Microbial population structures in the deep marine biosphere. Science. 2007;318(5847):97–100.

    Article  CAS  PubMed  Google Scholar 

  56. Wang Y, Qian PY. Conservative fragments in bacterial 16S rRNA genes and primer design for 16S ribosomal DNA amplicons in metagenomic studies. PLoS ONE. 2009;4(10):e7401.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Takai KEN, Horikoshi K. Rapid detection and quantification of members of the archaeal community by quantitative PCR using fluorogenic probes. Appl Environ Microbiol. 2000;66(11):5066–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Baker GC, Smith JJ, Cowan DA. Review and re-analysis of domain-specific 16S primers. J Microbiol Methods. 2003;55(3):541–55.

    Article  CAS  PubMed  Google Scholar 

  59. Gorrasi S, Pasqualetti M, Franzetti A, Gonzalez-Martinez A, Gonzalez-Lopez J, Muñoz-Palazon B, et al. Persistence of Enterobacteriaceae drawn into a marine saltern (saline di Tarquinia, Italy) from the adjacent coastal zone. Water. 2021;13(11):1443.

    Article  Google Scholar 

  60. Gorrasi S, Pasqualetti M, Braconcini M, Muñoz-Palazon B, Fenice M. Could Pontimonas Harbour Halophilic Members able to withstand very broad salinity variations? Microorganisms. 2022;10(4):790.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.

    Article  CAS  PubMed  Google Scholar 

  62. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Claesson MJ, O’Sullivan O, Wang Q, Nikkilä J, Marchesi JR, Smidt H, et al. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS ONE. 2009;4(8):e6669.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Gorrasi S, Pasqualetti M, Franzetti A, Pittino F, Fenice M. Vibrio communities along a salinity gradient within a marine saltern hypersaline environment (saline di Tarquinia, Italy). Environ Microbiol. 2020;22(10):4356–66.

    Article  CAS  PubMed  Google Scholar 

  65. Good IJ. The population frequencies of species and the estimation of population parameters. Biometrika. 1953;40(3–4):237–64.

    Article  Google Scholar 

  66. Shannon CE. A mathematical theory of communication. Bell Sys Tech J. 1948;27:379–423.

    Article  Google Scholar 

  67. Gini C. Variabilità e mutabilità. Contributo allo studio delle distribuzioni e delle relazioni statistiche. Bologna:Tipografia di Paolo Cuppini; 1912.

  68. McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10(4):e1003531.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Legendre P, Legendre L. Numerical Ecology. Amsterdam: Elsevier; 2012.

    Google Scholar 

  70. Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001;29(4):1165–88.

    Article  Google Scholar 

  71. Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science (1979). 2016;353:1272–7.

  72. Brandes U, Wagner D. Visone - analysis and visualization of social networks. Graph drawing Software. Springer; 2012. 321–40.

  73. Scoma A. Functional groups in microbial ecology: updated definitions of piezophiles as suggested by hydrostatic pressure dependence on temperature. ISME J 2021. 2021;15:7.

    Article  Google Scholar 

  74. Herndl GJ, Bayer B, Baltar F, Reinthaler T. Prokaryotic Life in the Deep Ocean’s Water Column. Ann Rev Mar Sci. 2023;15:461–83.

    Article  PubMed  Google Scholar 

  75. Huang JM, Wang Y. Genomic differences within the phylum Marinimicrobia: from waters to sediments in the Mariana Trench. Mar Genomics. 2019;50:100699.

    Article  PubMed  Google Scholar 

  76. Corinaldesi C. New perspectives in benthic deep-sea microbial ecology. Front Mar Sci. 2015;2:17.

    Article  Google Scholar 

  77. Zamkovaya T, Foster JS, de Crécy-Lagard V, Conesa A. A network approach to elucidate and prioritize microbial dark matter in microbial communities. ISME J. 2021;15(1):228–44.

    Article  PubMed  Google Scholar 

  78. Towner KJ. Biology of Acinetobacter spp. Acinetobacter: Microbiology, Epidemiology, Infections, Management. CRC Press; 2020. 13–36.

  79. Do HK, Kogure K, Simidu U. Identification of deep-sea-sediment bacteria which produce tetrodotoxin. Appl Environ Microbiol. 1990;56(4):1162–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Durand P, Benyagoub A, Prieur D. Numerical taxonomy of heterotrophic sulfur-oxidizing bacteria isolated from southwestern Pacific hydrothermal vents. Can J Microbiol. 1994;40(8):690–7.

    Article  CAS  Google Scholar 

  81. Takami H, Inoue A, Fuji F, Horikoshi K. Microbial flora in the deepest sea mud of the Mariana Trench. FEMS Microbiol Lett. 1997;152(2):279–85.

    Article  CAS  PubMed  Google Scholar 

  82. Maruyama A, Honda D, Yamamoto H, Kitamura K, Higashihara T. Phylogenetic analysis of psychrophilic bacteria isolated from the Japan Trench, including a description of the deep-sea species Psychrobacter pacificensis sp. nov. Int J Syst Evol Microbiol. 2000;50(2):835–46.

    Article  CAS  PubMed  Google Scholar 

  83. Bharathi PL, Nair S. Rise of the dormant: simulated disturbance improves culturable abundance, diversity, and functions of deep-sea bacteria of Central Indian Ocean Basin. Mar Georesources Geotechnol. 2005;23(4):419–28.

    Article  Google Scholar 

  84. Kumar AG, Vijayakumar L, Joshi G, Peter DM, Dharani G, Kirubagaran R. Biodegradation of complex hydrocarbons in spent engine oil by novel bacterial consortium isolated from deep sea sediment. Bioresour Technol. 2014;170:556–64.

    Article  CAS  Google Scholar 

  85. Zhang L, Kang M, Xu J, Xu J, Shuai Y, Zhou X, et al. Bacterial and archaeal communities in the deep-sea sediments of inactive hydrothermal vents in the Southwest India Ridge. Sci Rep. 2016;6(1):1–11.

    Article  CAS  Google Scholar 

  86. Nitahara S, Kato S, Usui A, Urabe T, Suzuki K, Yamagishi A. Archaeal and bacterial communities in deep-sea hydrogenetic ferromanganese crusts on old seamounts of the northwestern Pacific. PLoS ONE. 2017;12(2):e0173071.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Zhao X, Liu J, Zhou S, Zheng Y, Wu Y, Kogure K, et al. Diversity of culturable heterotrophic bacteria from the Mariana Trench and their ability to degrade macromolecules. Mar Life Sci Technol. 2020;2(2):181–93.

    Article  Google Scholar 

  88. Ma M, Gao W, Li Q, Han B, Zhu A, Yang H, et al. Biodiversity and oil degradation capacity of oil-degrading bacteria isolated from deep-sea hydrothermal sediments of the South Mid-Atlantic Ridge. Mar Pollut Bull. 2021;171:112770.

    Article  CAS  PubMed  Google Scholar 

  89. Zhang J, Sun QL, Zeng ZG, Chen S, Sun L. Microbial diversity in the deep-sea sediments of Iheya North and Iheya Ridge, Okinawa Trough. Microbiol Res. 2015;177:43–52.

    Article  PubMed  Google Scholar 

  90. Kong W, Zhao C, Gao X, Wang L, Tian Q, Liu Y, et al. Characterization and transcriptome analysis of a long-chain n-alkane-degrading strain Acinetobacter pittii SW-1. Int J Environ Res Public Health. 2021;18(12):6365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Jang GI, Hwang CY, Choi HG, Kang SH, Cho BC et al. Description of Spongiibacter borealis sp. nov., isolated from Arctic seawater, and reclassification of Melitea salexigens Urios 2008 as a later heterotypic synonym of Spongiibacter marinus Graeber 2008 with emended descriptions of the genus Spongiibacter and Spongiibacter marinus. Int J Syst Evol Microbiol. 2011;61(12):2895–2900.

  92. Li HJ, Zhang XY, Chen CX, Zhang YJ, Gao ZM, Yu Y, et al. Zhongshania antarctica gen. nov., sp. nov. and Zhongshania guokunii sp. nov., gammaproteobacteria respectively isolated from coastal attached (fast) ice and surface seawater of the Antarctic. Int J Syst Evol Microbiol. 2011;61(9):2052–7.

    Article  CAS  PubMed  Google Scholar 

  93. Naysim LO, Kang HJ, Jeon CO. Zhongshania aliphaticivorans sp. nov., an aliphatic hydrocarbon-degrading bacterium isolated from marine sediment, and transfer of Spongiibacter borealis Jang et al. 2011 to the genus Zhongshania as Zhongshania borealis comb. nov. Int J Syst Evol Microbiol. 2014;64(11):3768–74.

    Article  CAS  Google Scholar 

  94. On YG, Oh JS, Roh DH. Zhongshania marina sp. nov., isolated from deep-sea water. Int J Syst Evol Microbiol. 2019;69(2):542–546.

  95. Barbato M, Scoma A. Mild hydrostatic-pressure (15 MPa) affects the assembly, but not the growth, of oil-degrading coastal microbial communities tested under limiting conditions (5°C, no added nutrients). FEMS Microbiol Ecol. 2020;96(11):fiaa160.

    Article  CAS  PubMed  Google Scholar 

  96. Ribicic D, Netzer R, Hazen TC, Techtmann SM, Drabløs F, Brakstad OG. Microbial community and metagenome dynamics during biodegradation of dispersed oil reveals potential key-players in cold norwegian seawater. Mar Pollut Bull. 2018;129(1):370–8.

    Article  CAS  PubMed  Google Scholar 

  97. Murphy SM, Bautista MA, Cramm MA, Hubert CR. Diesel and crude oil biodegradation by cold-adapted microbial communities in the Labrador Sea. Appl Environ Microbiol. 2021;87(20):e00800–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Gutierrez T. Occurrence and roles of the obligate hydrocarbonoclastic bacteria in the ocean when there is no obvious hydrocarbon contamination. Taxonomy, Genomics and Ecophysiology of Hydrocarbon-Degrading Microbes. Springer; 2019. 337–52.

  99. Fan S, Wang M, Ding W, Li YX, Zhang YZ, Zhang W. Scientific and technological progress in the microbial exploration of the hadal zone. Mar Life Sci Technol. 2021;4(1):127–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  100. Gugliandolo C, Maugeri TL. Phylogenetic diversity of archaea in shallow hydrothermal vents of Eolian Islands, Italy. Diversity. 2019;11(9):156.

    Article  CAS  Google Scholar 

  101. Zhong H, Lehtovirta-Morley L, Liu J, Zheng Y, Lin H, Song D, et al. Novel insights into the Thaumarchaeota in the deepest oceans: their metabolism and potential adaptation mechanisms. Microbiome. 2020;8(1):1–16.

    Article  CAS  Google Scholar 

  102. Haro-Moreno JM, Rodriguez-Valera F, López-García P, Moreira D, Martin-Cuadrado AB. New insights into marine group III Euryarchaeota, from dark to light. ISME J. 2017;11(5):1102–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  103. Santoro AE, Richter RA, Dupont CL. Planktonic marine archaea. Ann Rev Mar Sci. 2019;11(1):131–58.

    Article  PubMed  Google Scholar 

  104. Huang WC, Liu Y, Zhang X, Zhang CJ, Zou D, Zheng S, et al. Comparative genomic analysis reveals metabolic flexibility of Woesearchaeota. Nat Commun. 2021;12(1):1–14.

    Article  CAS  Google Scholar 

  105. Liu X, Li M, Castelle CJ, Probst AJ, Zhou Z, Pan J, et al. Insights into the ecology, evolution, and metabolism of the widespread woesearchaeotal lineages. Microbiome. 2018;6(1):1–16.

    Article  Google Scholar 

  106. León-Zayas R, Novotny M, Podell S, Shepard CM, Berkenpas E, Nikolenko S, et al. Single cells within the Puerto Rico Trench suggest hadal adaptation of microbial lineages. Appl Environ Microbiol. 2015;81:8265–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  107. Zhao X, Luo H, He S, Yang B, Wei T, Hu Y, et al. Vertical distribution of size-fractionated bacterial communities in the water column of the Atacama Trench. Reg Stud Mar Sci. 2022;55:102470.

  108. Bayer B, Hansman RL, Bittner MJ, Noriega-Ortega BE, Niggemann J, Dittmar T, et al. Ammonia‐oxidizing archaea release a suite of organic compounds potentially fueling prokaryotic heterotrophy in the ocean. Environ Microbiol. 2019;21(11):4062–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  109. Kim JG, Gazi KS, Awala SI, Jung MY, Rhee SK. Ammonia-oxidizing archaea in biological interactions. J Microbiol. 2021;59(3):298–310.

    Article  CAS  PubMed  Google Scholar 

  110. Kitzinger K, Marchant HK, Bristow LA, Herbold CW, Padilla CC, Kidane AT, et al. Single cell analyses reveal contrasting life strategies of the two main nitrifiers in the ocean. Nat Commun. 2020;11(1):1–12.

    Article  CAS  Google Scholar 

  111. Bertagnolli AD, McCalmont D, Meinhardt KA, Fransen SC, Strand S, Brown S, et al. Agricultural land usage transforms nitrifier population ecology. Environ Microbiol. 2016;18:1918–29.

    Article  CAS  PubMed  Google Scholar 

Download references


The material was collected and sorted within the framework of the KuramBio II project supported by the PTJ (German Ministry for Science and Education), grant 03G0250A to Prof. Dr. Angelika Brandt, University of Hamburg, now Senckenberg Museum, Frankfurt, Germany.

We thank the crew of the RV Sonne for their help on board and all student helpers and technicians for support and help with sorting of the extensive expedition material.


Not applicable.

Author information

Authors and Affiliations



SG, AF and MF conceptualized the study. AB and MF coordinated the research activity planning and execution. AF, AB and MF were involved in funding acquisition, provided the resources and supervised the research activities. UM carried out the sampling. SG, AF, MP and MF carried out the investigation and the formal analysis. SG and MF wrote the original draft. All authors edited and corrected the final version of the manuscript and they approved the submitted version.

Corresponding authors

Correspondence to Susanna Gorrasi or Massimiliano Fenice.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have read the manuscript and consented to the publication.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.


Additional file 1: Fig. S1. Rarefaction curves of KKT abyssal-hadal samples for Bacteria (A) and Archaea (B) libraries; Fig. S2. Cluster Analysis (Bacteria dataset) of the KKT abyssal-hadal samples based on Hellinger distance; Fig. S3. Cluster Analysis (Archaea dataset) of the KKT abyssal-hadal samples based on Hellinger distance.


Additional file 2: Table S1. Bacterial genera found in the abyssal-hadal zones of the KKT and other trenches; Table S2. Archaeal genera found in the abyssal-hadal zones of the KKT and other trenches; Table S3. Bacterial ecological functions predicted through FAPROTAX analysis; Table S4. Archaeal ecological functions predicted through FAPROTAX analysis; Table S5. Abundance of co-occurrent bacterial and archaeal OTUs (network analysis) across the KKT abyssal and hadal samples.

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

Gorrasi, S., Franzetti, A., Brandt, A. et al. Insights into the prokaryotic communities of the abyssal-hadal benthic-boundary layer of the Kuril Kamchatka Trench. Environmental Microbiome 18, 67 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Bacterial communities
  • Archaeal communities
  • Deep Sea
  • Abyssal-hadal zone
  • Kuril Kamchatka Trench
  • 16S rRNA gene amplicon sequencing