Hongxi Zhang
Taoshu Wei1
Lisheng He
Yong Wang- 1Institute of Deep Sea Science and Engineering, Chinese Academy of Sciences, Sanya, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Sansha Track Ocean Coral Reef Conservation Research Institute Co., Ltd., Sansha, China
- 4Institute for Ocean Engineering, Shenzhen International Graduate School, Shenzhen, China
- 5Shenzhen Key Laboratory of Advanced Technology for Marine Ecology, Shenzhen International Graduate School, Tsinghua University, Shenzhen, China
The Yongle blue hole (YBH), situated in the South China Sea, represents a compelling subject of study in marine microbiology due to its unique redox-layered microbial ecosystems. However, the diversity and ecology of microbial eukaryotes within the YBH remains underexplored. This study endeavors to bridge this gap through the application of the in situ microbial filtration and fixation (ISMIFF) device to collect 0.22–30 μm microbial samples from 21 water layers of YBH. Subsequent extraction of 18S rRNA metagenomic reads of 21 metagenomes and 10 metatranscriptomes facilitated a comprehensive analysis of community structures. Findings revealed a pronounced superiority in the diversity and richness of eukaryotic microorganisms in the oxic zone compared to its suboxic and anoxic counterparts. Notably, Dinophyceae and Maxillopoda emerged as the predominant taxa based on the analysis of the 18S rRNA reads for the V4 and V9 regions, which showed stratification In their relative abundance and suggested their potential role in the thermo-halocline boundaries and oxic-anoxic interface. Specifically, In these eukaryotic microbial communities, Dinophyceae exhibited significant abundance at 20 m (20.01%) and 105 m (26.13%) water depths, while Maxillopoda was prevalent at 40 m (22.84%), 80 m (23.19%), and 100 m (15.42%) depths. A part of these organisms, identified as larvae and protists, were likely attracted by swarming chemosynthetic bacterial prey prevailing at the thermo-halocline boundaries and oxic-anoxic interface. Furthermore, the phylogenetic relationships of the major 18S operational taxonomic units (OTUs) showed a close adjacency to known species, except for three Dinophyceae OTUs. In conclusion, this study provides critical insights into the vertical distribution and transcriptional activity of <30-μm eukaryotic microbes, shedding light on the taxonomic novelty of eukaryotic microorganisms within the semi-enclosed blue holes.
1 Introduction
The Yongle blue hole (YBH), situated in the Xisha Islands of the South China Sea, presents a stratified aqueous environment comprising a superficial oxic layer, a relatively stable oxic-anoxic interface, and a deep anoxic layer (Yao et al., 2018; Xie et al., 2019), with a water depth of 300.89 m (Yao et al., 2020). Water below approximately 100 m, exhibiting seasonal fluctuations in depth, and contains a significant amount of hydrogen sulfide and methane, likely resulting from biodegradation (Xie et al., 2019). The unique topographical structure of the cave acts as a barrier to vertical mixing, preventing the transport of external seawater and particulate organic matter to the cave’s bottom (Li et al., 2018, 2020; Gao et al., 2022). Temperature and salinity within the water column exhibit distinct stratification, and remain relatively constant in the deeper layer (Xie et al., 2019). This results in a unique physical–chemical characteristics, such as a strong thermo-halocline, a highly stratified water column, and thick anoxic and hydrogen sulfided layers. Such conditions are believed to support the development of a distinct ecosystem within the YBH. The stable anoxic environment provided by the YBH serves as an important model for studying ancient anoxic oceans and the early evolutionary history of eukaryotes (Li et al., 2010; Mentel and Martin, 2010; Zhang et al., 2024).
Eukaryotic organisms residing beneath the redox layer encounter environmental conditions markedly distinct from those in the oxygen-rich water layer (Behnke et al., 2006). The formation of symbiotic associations among the domains Bacteria, Archaea, and Eukarya in anoxic water columns is reported to be prevalent, positing a potential adaptation mechanism for eukaryotic microbes to exploit ecologically challenging habitats (Orsi et al., 2012). Surveys of marine blue holes have revealed numerous new species, including various types of crustaceans and foraminifera (Lejeusne and Chevaldonné, 2010; Boxshall and Jaume, 2012). Crustaceans have been found to dominate in some blue holes in the Bahamas and the coast of Mexico with a significant biodiversity, with many species unique to blue holes and not found in the surrounding waters (Kornicker et al., 2007; Iliffe and Kornicker, 2009). The interior of the holes contains a large number of planktons that are significantly different from those found on the surrounding outer reef slope, exhibiting notably higher biodiversity and abundance (Chen et al., 2018). Due to the challenges of studying the morphological and ultrastructural characteristics of microbial eukaryotes, analytical methods based on 18S rRNA gene clone libraries have been used for taxonomic analysis in different marine ecosystems (Zuendorf et al., 2006; Behnke et al., 2010; Wang et al., 2014). Using amplicons of 18S rRNA genes, the microbial community structures in the water column at different depths of YBH were preliminarily analyzed (Liu et al., 2019). It was observed that the abundance of eukaryotic microbes in the cave is significantly higher than that in the outer reef slope, alongside a large number of endemic species have been found, indicating the need for further study (Liu et al., 2019). It is challenging to recover intact eukaryotic microbial cells from the deep anoxic water layer (Wylezich and Jürgens, 2011) and, furthermore, there have been observations of biases introduced in the analysis of the community structure when using18S rRNA gene amplicons (Yeh et al., 2021). To address these challenges, metagenomics, based on the extraction of 18S ribosomal DNA reads from metagenomes, has been used to study microbial diversity and overcome these issues (Logares et al., 2014). This method allows the detection of low-abundance taxa and has been applied in the analysis of bathypelagic microbial eukaryotes (Pernice et al., 2016). It has been noted that the microbial eukaryote communities of the anoxic and suboxic water column are primarily influenced by oxygen concentration, as well as the sampling method (Duret et al., 2015). The sampling method employed could have affected the microorganisms in the deep layers of YBH with low oxygen and light density (Torres-Beltran et al., 2019). An in situ microbial filtration and fixation (ISMIFF) apparatus, which is mounted on manned submersibles, remotely operated vehicles (ROVs), and landers, has the capability to operate automatically in the full ocean depth (Gao et al., 2019; Wang et al., 2019; Wei et al., 2020). This technology was utilized to enrich and fix microorganisms in YBH in situ, allowing for a better understanding of the community structure and in situ metabolism of prokaryotic microbes without being impacting by environmental changes during sampling (Zhang et al., 2023, 2024).
The pico- and nano-sized (0.22–20 μm) eukaryotes (Ceccarelli et al., 2013) play a crucial role in marine plankto ecosystems as both primary producers and major bacterial grazers (Li et al., 2023). Also, the 0.2–30 μm diameter fraction includes pico- (0.2–2 μm), nano- (2–20 μm), and small microplankton (20–30 μm), representing a wide range of protistan species that are studied to analyze the structure of planktonic microbial eukaryotic communities (David et al., 2021). Nanoeukaryotes form part of the diet for some plankton larvae, impacting the community structure of eukaryotic microorganisms (Lindeque et al., 2015). In this study, we employed 18S rRNA metagenomic Illumina reads (miTags) to investigate the eukaryotic communities in the water column of YBH, which provides evidence for the uniqueness and distribution of eukaryotic species in suboxic and anoxic zones of blue holes.
2 Materials and methods
2.1 Sampling
All samples were obtained from YBH (16°31′30″N; 111°46′05″E) in the South China Sea. Water samples were collected by the ISMIFF apparatus (Wang et al., 2019) at 21 depths from the surface to the bottom layer inside the YBH. In the April cruise, 16 samples were collected in the depth range of 1–290 m (Supplementary Table S1). In the September cruise, we collected samples from five additional suboxic and anoxic water layers that were not sampled in April, including depths of 96.5 m, 98.5 m, 99.5 m, 100.5 m, and 250 m. ISMIFF was used to filter ~30–60 L seawater by a 0.22-μm pore size polycarbonate membrane (Millipore, Massachusetts, United States) with a prefiltration through a 30-μm mesh net. Immediately after ISMIFF was recovered, the membrane was taken out from the filtration chamber and put into a 15 mL centrifuge tube containing 3 mL RNAlater (Ambion, Carlsbad, CA, United States). The tube was then frozen in liquid nitrogen immediately and stored at −80°C until further processing. All operations were completed within 30 min. An SBE 37SMP-ODO conductivity-temperature-depth (CTD) unit and a dissolved oxygen (DO) sensor (Sea-Bird, Bellevue, WA, United States) were used to detect the temperature, salinity, and DO. pH was measured directly with Hach model HQ 40d portable combination meter (HACH, Loveland, CO, United States).
2.2 Nucleic acids extraction, metagenome and metatranscriptome sequencing
Total genomic DNA and RNA were extracted from 21 ISMIFF samples using the DNA/RNA co-extraction kit (Tiangen, Beijing, China) following the manufacturer’s instructions. DNA and RNA concentrations were measured using a Qubit™ 2.0 Fluorometer (Invitrogen, Carlsbad, CA, United States). 100 ng DNA was used as a template for the construction of a metagenomic high-throughput sequencing library. In the process of library construction, the long DNA fragments were ultrasonically broken by the Covaris M220 instrument (Covaris, Woburn, Massachusetts, USA), and the length of the inserted fragments in the library was about 350 bp. Then, the metagenomic high-throughput library was constructed by using the VAHTS Universal DNA Library Prep Kit for Illumina V3 (Vazyme, Nanjing, China) according to the instructions. 5 μL constructed library was used for 1% agarose gel electrophoresis to verify the library quality. The DNA library was sequenced using an Illumina Novaseq 6000 platform (2×150 bp).
Using the same method, we extracted total RNA from the 11 ISMIFF samples in-situ filtered with 0.22 μm pore size membranes. Qubit™ 2.0 Fluorescence meter was employed to measure RNA concentration. RNA extracts were treated with DNase I (Tiangen, Beijing, China) as instructed by the manufacturer to ensure DNA removal. PCR was performed to detect DNA residues using the universal primers 341F (5′-CCTAYGGGRBGCASCAG-3′) and 802R (5′-TACNVGGGTA TCTAATCC-3′), and the RNA extracts as a template. Each 50 μL PCR reaction contains 25 μL of 2 × Taq PCR premixed reagent (Tiangen, Beijing, China) and 1 μL of 10 μM primers each. All PCR reaction systems were first heated at 94°C for 3 min, followed by 30 cycles of amplification: 94°C (30 s), 55°C (30 s), and 72°C (1 min), and maintenance at 72°C for 5 min at the end of the cycle. 1% agarose gel electrophoresis was used to confirm that DNA was removed completely. Then, according to instructions for the VAHTS® Universal V8 RNA-seq Library Prep Kit for Illumina (Vazyme, Nanjing, China) kit, the total RNA after DNA removal was used as a template to synthesize double-stranded cDNA and construct a metatranscriptome library. Among them, rRNA was not removed from the total RNA. We then fragmented RNA using conditions of 85°C for 5 min according to the instructions of the reagent kit. 5 μL constructed library was used for 1% agarose gel electrophoresis to verify the library quality. The cDNA library was sequenced using the Illumina Novaseq 6000 platform (2 × 150 bp).
2.3 Community composition analysis
Raw Illumina sequencing data of the 21 metagenomes and 11 metatranscriptomes were evaluated using FastQC (v0.11.8).1 They were then trimmed to remove adapters and then filtered using fastp (v0.23.1) (Chen et al., 2018) with parameters (-w 24 -c -q 20 -u 20 -g -W 5–3 -l 50). Reads with low quality (assigned by a quality score < 20 for >20% of the read length), those shorter than 50 bp, and unpaired reads were removed. Data with a high level of duplications (>5%) were pre-processed with FastUniq (v1.1) (Xu et al., 2012) to eliminate duplicated reads.
Ribosomal RNA miTags (5.8S, 18S, and 28S) were extracted from the clean metagenomic and metatranscriptomic reads using rna_hmm3.py (Huang et al., 2009), which employed HMMER (v.3.1b2) (Mistry et al., 2013) to predict ribosomal RNA gene fragments from both forward and reverse metagenomic reads. 18S miTags (≥100 bp) were identified using an in-house Python script (Zhou et al., 2022) and matched to the V4 and V9 variable regions of the 18S rRNA gene using HMMsearch software (Eddy, 2011) against V4 and V9 Hidden Markov Models. The 18S miTags were imported into QIIME 2 (v.2022.2) (Bolyen et al., 2019) with the setting of --type “SampleData [Sequences]” and clustered into operational taxonomic units (OTUs) with ≥ 97% similarity using the VSEARCH pipeline (Rognes et al., 2016). The representative miTags were classified using the classify-sklearn command in QIIME 2 in reference to the SILVA SSU 138.1 database (Quast et al., 2013; Glockner et al., 2017). Furthermore, mTAGs (Salazar et al., 2021) were used for additional classification of these miTags, referencing the PR2 (Protist Ribosomal Reference) database (Guillou et al., 2013). After annotation, to account for the differences in sequencing depth between samples, the q2-feature-table command in QIIME2 was used to normalize the feature table, obtaining relative abundance at different taxonomic levels. Three metrics, including Good’s coverage, Chao1, and Shannon index at a 3% dissimilarity were calculated using the QIIME 2 package. The relative abundance of the 18S miTags in 21 metagenomes was imported to the “vegan” package (v.2.5–4) in R to estimate Bray–Curtis dissimilarity index for Principal Coordinate Analysis (PCoA). The dissimilarity of eukaryotic communities in the studied samples of different depths was also analyzed by the unweighted pair group method with arithmetic average (UPGMA) clustering, using the upgma function in the R Vegan package.
2.4 Phylogenetic analysis
The reference 18S rRNA gene sequences of Dinophyceae and Maxillopoda were downloaded from the National Center for Biotechnology Information (NCBI) database and the most abundant Dinophyceae and Maxillopoda OTUs in this study were used to construct a phylogenetic tree. The comparison was made using MAFFT (v7.505, setting: –maxiterate 1,000 –localpair) (Katoh and Standley, 2013) and then trimAl (v1.4, setting: -automated1-phylip) (Capella-Gutierrez et al., 2009) was used to construct the alignment results. The maximum likelihood phylogenetic tree of the 18S rRNA genes was constructed using IQ-TREE 2 (v2.2.0, setting: -m MFP -B 1,000 -alrt 1,000 -T AUTO) (Minh et al., 2020). Phylogenetic trees were visualized with iTOL.2
3 Results
3.1 Environmental parameters in sampling stations
When sampling from YBH in 2021, we used sensors to measure temperature, pH, salinity, chlorophyll, and dissolved oxygen concentration (DO) in the blue hole (Figure 1). Except for the increased dissolved oxygen concentration in the September cruise, other environmental parameters were detected in the April cruise. The temperature, pH, and dissolved oxygen obtained in this study showed that the readings of these sensors decreased with increasing depth across the water column, and the salinity showed an opposite trend. The temperature decreases with depth from 28.43°C on the surface to 16.30°C at about 90 m and tends to be stable (Figure 1A). The salinity gradually increased from 33.87 at the surface to 34.45 at 95 m and then basically remained stable in the deeper water as previously reported (Xie et al., 2019; He et al., 2020). There are several possible thermo-haloclines due to lower temperatures and fluctuating changes in salinity, and they are located at 20–30 m, 40–50 m, and 80–100 m (Figure 1A). Chlorophyll value peaked at nearly 30 m, reaching 2.88 μg/L. It decreased rapidly with an increment of water depth, and a secondary peak (1.37 μg/L) occurred at ~45 m and then decreased to nearly 0 at ~55 m. Dissolved oxygen concentrations varied seasonally in YBH (Figure 1B). Suboxic layer (0–0.5 mg/L dissolved oxygen) occurred at 105–110 m and 92–101 m in April and September, respectively. Dissolved oxygen in April began to decrease rapidly at about 80 m (~3.37 mg/L), while that in September showed a generally declining trend. Altogether, the YBH water column contained more dissolved oxygen in April than in September. Moreover, in September, YBH had a stronger temperature and redox gradients than in April. Previously, it has been reported that there are multiple thermoclines in YBH, while their depth varies seasonally and the conditions of the anoxic layer are quite stable (Yao et al., 2018; Xie et al., 2019).
Figure 1. Hydrological parameter profiles in YBH. Profiles of seawater temperature, salinity (A), pH, Chlorophyll, and DO were measured in April and September of 2021 (B) inside YBH. The layers with a blue shade represent the thermo- and halocline layers (A).
During the two cruises, ISMIFF was used to obtain in-situ filtration membranes from 16 water layers in a depth range of 1–290 m in April of 2021 (Spring sample), and from five water layers of 96.5–100.5 m and 250 m in September of 2021 (Autumn sample). Among them, the four water layers from 96.5–100.5 m were combined with suboxic samples. A total of 162.48 Gbp and 23.61 Gbp data for metagenomes and metatranscriptomes, respectively, were obtained for all the samples (Supplementary Table S1).
3.2 Community composition of eukaryotes in YBH
A total of 16,906 18S miTags assigned for the V4 and V9 regions were extracted from clean Illumina reads of the 21 metagenomes to decipher the eukaryotic microbial communities of YBH. Since the number of 18S V4 and V9 miTags in some samples of the suboxic and anoxic layers was less than 100 (Supplementary Table S1), these 18S miTags of 96.5–100.5 m, 110–150 m, and 180–290 m metagenomes were combined into three individual groups for analysis. To further evaluate the microbial activity, 18S V9 miTags were successfully extracted from clean Illumina reads from 10 metatranscriptomes, which resulted in a total of 381,446 miTags (Supplementary Table S1). The alpha-diversity indices showed that the diversity and richness of the eukaryotic microorganisms in the oxic layer of YBH were significantly higher than those in the suboxic and anoxic layers (U test, p < 0.001) (Supplementary Table S1). PCoA plotting showed the divergence of microbial communities at the genus and order levels between different layers (Figure 2; Supplementary Figure S2). The first two principal components accounted for 54.30 and 17.86% of the total variation in the microbial community composition, respectively. In the oxic zone (1–105 m), depth was the main factor driving the division of eukaryotic microbial community structures, explaining 13.35% of the variance in the communities. In contrast, the suboxic (96.5–100.5 m) and anoxic layers (110–290 m) are relatively aggregated, indicating that the community structures resembled each other in the water layers with low oxygen content.
Figure 2. PCoA analysis of eukaryotic community structures. The relative abundances of genera based on the classification of 18S miTags extracted from 21 metagenomes were used for the PCoA analysis. The genera represented by less than 2% in the OTUs were ignored. The sample IDs are referred to Supplementary Table S1.
In total, 20 major eukaryotic phyla (relative abundance >5%) including Arthropoda, Dinoflagellata, Chlorophyta, Porifera, and Protalveolata were consistently detected in all the samples, with the removal of unassigned OTUs (~6.34% of all samples) and Eukaryota OTUs that were not classified into phyla (~47.98% of all samples) in the plots (Supplementary Figure S3). The samples were overwhelmed by Arthropoda (~10.49%), followed by Dinoflagellata (~8.41%). The eukaryotes were further divided into 25 major orders (Supplementary Figure S4) and 28 major genera (relative abundance >5%) (Figure 3A). The dominant phyla in the metatranscriptomes were roughly the same as those in the metagenomes, with the exclusion of unassigned OTUs (~2.69%) and Eukaryota OTUs that were not classified into phyla (~41.50%) in the plots (Figure 3B; Supplementary Figure S3B). In addition, we searched the OTUs against the PR2 database, which showed similar results with the classification based on the Silva database (Supplementary Figure S5). Additionally, we showed the structure of the protozoan community after the removal of the metazoans (Arthropoda, Porifera, Annelida, Cnidaria, Xenacoelomorpha, and Mollusca), in which Dinoflagellata was overwhelmingly dominant (~39.23% of all samples) (Supplementary Figure S6). However, copepod groups have been found to be in higher abundance in the hole (Chen et al., 2018), and 18S rRNA gene amplicon data have been used to assess the metazoan diversity in the pico-nanoplanktonic samples (López-Escardó et al., 2018). Therefore, in order to comprehensively understand the diversity of eukaryotic microorganisms, subsequent analyses included data that encompassed metazoan.
Figure 3. Eukaryotic community structure at genus level based on 18S miTags in metagenomes and metatranscriptomes. (A) Eukaryotic community structure in 21 metagenomes. (B) Eukaryotic community structure in 10 metatranscriptomes. 18S miTags for the V9 region were extracted from the metagenomes and metatranscriptomes for the generation of OTUs with a 97% similarity. The genera that accounted for less than 5 and 2% of the OTUs for metagenomes and metatranscriptomes, respectively, were grouped into “minor group.” Sampling depth (m) is indicated by the number or number range in the IDs that are described in Supplementary Table S1.
The presence of short 18S miTags in this study may have contributed to the potential inaccuracy of eukaryotic taxonomic classification. To address this, we compared the 18S V4 miTags with V9 miTags in the consistency of eukaryotic microbial community structures. Our results obtained using the V4 miTags of the 18S rRNA genes (Supplementary Figure S7) were similar to those based on the V9 ones. Both regions demonstrated the same dominant taxa (Dinophyceae and Maxillopoda). The hypervariable V9 region has been used to independently evaluate eukaryotic microbial diversity at the population level with the availability of a good reference database (Obiol et al., 2020).
Maxillopoda (Arthropoda) were dominant in 40 m (22.84%), 80 m (23.19%), and 100 m (15.42%) water depths, respectively. In the anoxic zone of this study, a certain abundance of Maxillopoda (4.62%) was identified in the metagenomes and was transcriptionally active (about 10.66% of the metatranscriptome-based communities) (Figure 3B). Dinophyceae was the main group of eukaryotic algae with the highest abundance in 20-m (20.01%) and 105-m (26.13%) depth metagenomes, and dominated the transcriptomes from 105 m depth (15.35%) (Figure 3B). The class Sordariomycetes (Ascomycota) prevailed at 98.5 m depth, as they accounted for 16.74% of the metagenome and 18.12% of the transcriptome of the suboxic samples (Figure 3B). Similar to the surface layer, Chlorophyta (11.59%) and Porifera (6.54%) were also rich in the deep anoxic layers (180–290 m depths) (Figure 3A). Interestingly, Porifera showed transcriptional activity at 290 m depth (genus Haplosclerida, occupying about 0.06% of the transcribed 18S miTag V9 reads), and their distribution is currently known to be restricted in shallow water (Redmond et al., 2011). Moreover, the taxonomic status of the Porifera remains to be determined.
UPGMA clustering was applied to compare the community similarity and dissimilarity of different depths based on metagenomic 18S V9 miTags. Our result showed that the samples from different oxygenated layers could be divided into two clades, the first consisting of those sampled from 40 m, 80 m, and 100 m depths in April of 2021 and the second comprising the rest (Figure 4). Compared with suboxic and anoxic layers, the community structures of the oxic layers were more similar. It suggests that oxygen availability plays a significant role in shaping eukaryotic community structures, and the eukaryotic communities in the oxic layers exhibit similar composition patterns. At the sampling depths (40 m, 80 m, and 100 m) with a boundary of thermo-halocline (Figure 1A), the eukaryotic communities resembled the prevalence of Arthropoda in these layers, probably accounted for the thermo-halocline and the oxic-anoxic interface that hinder the vertical movement of eukaryotes (Lovejoy et al., 2006; Weber et al., 2014; Edgcomb, 2016).
Figure 4. UPGMA clustering of the eukaryote communities in YBH. The clustering was based on the relative abundance of eukaryote communities at genus level in 21 metagenomes. The genera represented by less than 2% in the OTUs were ignored. The p-values (%) were calculated using the pvclust function, where the red number represents the AU (approximately unbiased) p value of the cluster, and the green number represents the BP (bootstrap probability) of the cluster. The sample IDs are referred to Supplementary Table S1.
3.3 Phylogenomics of Dinophyceae and Maxillopoda
To learn a more detailed taxonomic affiliation of the most abundant Dinophyceae OTUs of 18S V9 miTags, we reconstructed phylogenetic trees using their representative joint reads from the 20 m and 105 m metagenomes, respectively (Figure 5). The reference sequences of the cultured strains were selected from the NCBI database. These Dinophyceae OTUs were almost all adjacent to known strains with high bootstrap values except for OTU1, OTU2, and OTU4 which occupied, on average, 2.93, 2.93, and 1.17% of the metagenome-based communities at 20 m and 105 m depths. These OTUs (0.92% ~ 2.93%) were scattered into different families of Dinophyceae. Syndiniales group I has been found to survive in anoxic and suboxic ecosystems (Guillou et al., 2010), and the OTUs from this study were not assigned to Syndiniales. Similarly, the most abundant Maxillopoda OTUs showed a phylogenetical affinity with the references from subclass Copepoda that were prevalent in YBH (Chen et al., 2018), indicating that this taxon is a dominant, highly active member in the oxic zone. Although most of the OTUs could be found in an adjacent position of Calanoida, Harpacticoida, and Cyclopoida in the phylogenetic tree (Figure 6), OTU2 (~3.16%) was aggregated with uncultured species. In our phylogenetic tree, Arthropoda OTU1 (~7.38%), OTU3 (~2.06%), and OTU5 (~2.30%) were proximate to the orders from Harpacticoida and Poecilostomatoida.
Figure 5. Maximum-likelihood phylogenetic tree of dominant OTUs of Dinophyceae class in YBH samples. The representative sequences of the dominant OTUs based on 18S V9 miTags were aligned with the closest references from NCBI for the construction of the ML tree. Bootstrap values above 50 based on 1,000 replicates are shown as dots with different sizes. The highest relative abundance of these OTUs in themetagenomes for 20-m and 105-m depths is indicated as a percentage in the parentheses.
Figure 6. Maximum-likelihood phylogenetic tree of dominant OTUs of Maxillopoda class in YBH samples. The representative sequences of the dominant OTUs based on 18S V9 miTags were aligned with the closest references from NCBI for the construction of ML tree. Bootstrap values above 50 based on 1,000 replicates are shown as dots with different sizes. The highest relative abundance of these OTUs in the metagenomes for 40-m, 80-m, and 100-m depths is indicated as a percentage in the parentheses.
4 Discussion
In the present study, we investigated the stratified distribution and diversity of eukaryotic microbes in the water column of YBH using metagenomes and metatranscriptomes. It was observed that plankton in the cave exhibited diurnal vertical movement, particularly in the water layer above 90 m. This area has a strong thermo-halocline that may have hindered the zooplankton’s vertical migration (Chen et al., 2018), possibly explaining the higher abundance of arthropods at the depths of 80 m and 100 m compared to 90 m. The Arthropoda phylum showed high transcriptional activity at depths of 110–290 m (15.79%), indicating the presence of arthropods in the anoxic layer. Copepods with high abundance have been observed through microscopy (Chen et al., 2018), supporting the existence of Maxillopoda in the anoxic zone. However, it has been reported that the oligotrophic condition in the deep anoxic environment of the Cariaco Basin limits the survival rate of arthropods (Stoeck et al., 2003). Copepoda is a group of small Arthropoda with prosome adult length less than 2 mm (Roura et al., 2018). The smallest Copepoda is about 20 μm in size (Steinberg and Landry, 2017). Previous studies on Copepoda microbes from the ocean used 64-μm and even larger mesh nets (Shoemaker et al., 2019; Barth-Jensen et al., 2022; Sugai et al., 2023). In this study, we used two sets of filtration membranes to capture Arthropoda in sizes of 0.22–30 μm. However, from our knowledge, the adults of Harpacticoida and Poecilostomatoida were mostly larger than 30 μm (El-Rashidy and Boxshall, 1999; Dahms and Qian, 2004). Therefore, we believe that some of the Arthropoda captured were larvae harvesting bacterial preys in YBH. Furthermore, we cannot exclude the possibility of encountering novel Arthropoda species smaller than 30 μm in their adult stage. They were likely attracted by swarming chemosynthetic bacterial preys prevailing at the thermo-halocline boundaries and oxic-anoxic interface (Chen et al., 2018).
The growth of eukaryotic algae is affected by light intensity. However, excessive light can inhibit algae growth, as previously reported (Takahashi et al., 1973; Liu et al., 2019). The algae bloom at 20 m depth might be related to optimal light intensity in this study. The higher transcriptional activity of Dinophyceae may indicate an unknown branch capable of growing in lightless and low-oxygen environments. Dinophyceae that have been reported were largely found in the photic zone (Luo et al., 2022). Dinoflagellates are believed to be capable of phagotrophic nutrition (Sherr et al., 2007). Most of the Dinophyceae sequences were closely related to known cultured representatives or environmental clones from other hypoxic marine systems (Mariager Fjord) and an oligotrophic mountain lake (Lake Tovel) (Di Giuseppe and Dini, 2004; Zuendorf et al., 2006). Further genomic and transcriptomic studies of these OTUs will provide evidence for their metabolic activities in YBH.
The marine fungus Sordariomycetes is commonly found in marine sediments and is often associated with corals and sponges or parasitic for algae (Jones et al., 2019). Its prevalence at a depth of 98.5 m might be due to the abundance of organic particles for degradation or its parasitism with eukaryotic algae. Chlorophyta, a class of photosynthetic green algae, was classified as class Clade_VII and Ulvophyceae, which are typically found in euphotic environments (Hodgson and McDermid, 2000). However, their low abundance in the metatranscriptomes of deep anoxic layers (180–290 m) suggests that they may be less active in these conditions. Further investigation is needed to determine whether the Chlorophyta were permanent inhabitants of the YBH’s deep layer. The SCM38C39 marine group (phylum Euglenozoa, order Diplonemea) was found in the 100.5 m suboxic layer. While it accounted for a small percentage of the community (<2%) based on the metagenomics reads, its high activity (18.45%) in based on the metatranscriptomics reads. Diplonemea, known for preying on bacteria, is one of the most abundant groups of marine planktonic protists (Flegontova et al., 2020). It possesses enzymes for anaerobic adaptation (NADP+ oxidoreductase, fumarate reductase, enoyl- coenzyme A reductase, wax ester synthesis enzymes), which enables it to thrive in low-oxygen environments (Škodová-Sveráková et al., 2021). Diplonemea have been found in the anoxic Cariaco Basin (Schoenle et al., 2021)and might be more active in suboxic zones compared to other protists, as indicated by high transcriptomic levels.
Eukaryotic microalgae, specifically the Dinophyceae, dominated the water column at depths of 20 m and 105 m, and they might contribute significantly to carbon input in the YBH ecosystem. These microorganisms, along with Cyanobacteria on the surface, produce organic carbon, which provides important energy for the entire food web (He et al., 2020; Yao et al., 2020). Eukaryotic microorganisms, including Dinophyceae and Maxillopoda, play a role in nutrient cycling within the YBH. The presence of swarming chemosynthetic bacterial prey at the thermo-halocline boundaries and oxic-anoxic interface attracts larvae and protists, including Maxillopoda. These eukaryotic microorganisms likely feed on the bacterial preys (Pernthaler, 2005), contributing to the transfer of energy and nutrients within the ecosystem. Additionally, Maxillopoda may serve as a food source for higher trophic levels, such as fish and larger invertebrates (Suárez-Morales, 2015). The presence of different species at specific depths suggests that the blue hole provides distinct ecological niches for these organisms. This diversity of eukaryotes contributes to the overall biodiversity of the ecosystem and can influence species interactions, such as competition and predation, which shape the structure and dynamics of the community.
Other blue holes and typical anoxic water columns have similar sulfided environments such as that in the Black Sea (Edgcomb et al., 2014; Cabello-Yeves et al., 2021). There are abundant ciliates and flagellate taxa (mainly stramenopiles and dinoflagellates) in anoxic, especially sulfidic water layers, and the species diversity in sulfide layers is higher than that in hypoxic samples (Wylezich and Jürgens, 2011). By contrast, the main dominant phyla of the suboxic and anoxic sulfide water column in YBH are Arthropoda and Dinoflagellata as shown by this study. In YBH, the anoxic layer had lower species diversity compared to the suboxic layer. In the anoxic layer of YBH, ciliates mainly existed at 290 m depth (6.06%), belonging to the class Oligohymenophorea. However, the genus Cyclidium (Oligohymenophorea) which is the dominant bacterial group in the sulfidic layer of the Black Sea, may be a potential bacterial consumer and preys on sulfur-reducing bacteria. Compared with the Black Sea, ciliates are not eukaryotic microorganisms that dominate YBH. In summary, the eukaryotic microorganisms found in the anoxic and sulfided environment of YBH provide new evidence for the adaptation and evolution of early eukaryotes to anoxic habitats (Porter et al., 2018).
There are a large number of eukaryotes that have not been classified into phyla using current references in public databases, which can indicate that there are more novel species in YBH. We suggest more efforts for the exploration of the anaerobic eukaryotic microbes in the anoxic layer of YBH. In addition, although it is a fact that the size of Maxillopoda in general is very small, they are multicellular organisms. Their dominance in the present study might be caused by having multiple copies of the rRNA gene in their genome (Hu et al., 2016; Ibarbalz et al., 2019). However, their high abundance has been identified through morphological observation, a method that is beyond doubt for confirming their presence (Chen et al., 2018). Moreover, we also filtered the samples from over 30 L to avoid sampling drift. Nevertheless, this is an unavoidable limitation of the metagenomics method for marker gene analysis. Since we only analyzed samples <30-μm and included only protists and some larval communities of large eukaryotes (e.g., Maxillopoda, etc.), there was a lack of larger eukaryotes to show the dependence between prokaryotes and eukaryotes in the water column of the hole. Additional analysis will be conducted on eukaryotic communities within samples larger than 30-μm, in order to achieve a more systematic and comprehensive examination of the distinct ecosystem within the YBH.
5 Conclusion
This study analyzed the eukaryotic communities in various chemical gradient layers of YBH using metagenomics and metatranscriptomics and showed the dominant position of Dinophyceae and Maxillopoda in different layers. Dinophyceae were the most abundant at the 20 m and 105 m water depths, while Maxillopoda were dominant at the 40 m, 80 m, and 100 m. Some of these organisms were likely larvae and protists attracted by the chemosynthetic bacterial prey at the thermo-halocline boundaries and oxic-anoxic interface. The environmental barrier of protozoa in the hole hinders the vertical movement of protozoa and forms a stable ecological niche, which provides a microbial oasis for the stability of the material cycle in the hole. The phylogenetic relationships of the major 18S OTUs were mostly adjacent to known species, except for three Dinophyceae OTUs, pointing to the potential taxonomic novelty of these organisms. Consequently, our study has provided a detailed look at the vertical distribution of <30-μm eukaryotic microbes in the YBH, underscoring the taxonomic novelty and ecological importance of eukaryotic microorganisms in these unique semi-enclosed blue holes. The YBH hosts numerous endemic species that require further studied in terms of taxonomy and metabolic activities. In summary, this study presents important evidence for existence of eukaryotic microorganisms in suboxic and anoxic layers of marine blue hole ecosystems. Future efforts may further explore the metabolic pathways, energy acquisition strategies, and material cycling mechanisms of these eukaryotic microorganisms. This will ultimately provide in-depth understanding of the adaptability, diversity, and origin of eukaryotic organisms.
Data availability statement
18S miTags of YBH are available at the NCBI Sequence Read Archive (SRA) under BioProject PRJNA900714 (BioSample accession numbers SAMN31699490–SAMN31699510).
Author contributions
HZ: Conceptualization, Formal analysis, Visualization, Writing – original draft. TW: Formal analysis, Resources, Writing – review & editing. QL: Formal analysis, Resources, Writing – review & editing. LF: Resources, Writing – review & editing. ML: Writing – review & editing. LH: Resources, Writing – review & editing. YW: Conceptualization, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study is supported by National Natural Science Foundation of China (42376149) and Shenzhen Key Laboratory of Advanced Technology for Marine Ecology (ZDSYS20230626091459009).
Acknowledgments
We are grateful to Sansha Track Ocean Coral Reef Conservation Research Institute Co., Ltd. for their help in sample collection.
Conflict of interest
LF was employed by Sansha Track Ocean Coral Reef Conservation Research Institute Co., Ltd.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1420899/full#supplementary-material
Footnotes
References
Barth-Jensen, C., Daase, M., Ormańczyk, M. R., Varpe, Ø., Kwaśniewski, S., and Svensen, C. (2022). High abundances of small copepods early developmental stages and nauplii strengthen the perception of a non-dormant Arctic winter. Polar Biol. 45, 675–690. doi: 10.1007/s00300-022-03025-4
Behnke, A., Barger, K. J., Bunge, J., and Stoeck, T. (2010). Spatio-temporal variations in protistan communities along an O2/H2S gradient in the anoxic Framvaren Fjord (Norway). FEMS Microbiol. Ecol. 72, 89–102. doi: 10.1111/j.1574-6941.2010.00836.x
Behnke, A., Bunge, J., Barger, K., Breiner, H.-W., Alla, V., Stoeck, T. J. A., et al. (2006). Microeukaryote community patterns along an O2/H2S gradient in a supersulfidic anoxic fjord (Framvaren, Norway). Appl. Environ. Microbiol. 72, 3626–3636. doi: 10.1128/AEM.72.5.3626-3636.2006
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 37:1091. doi: 10.1038/s41587-019-0252-6
Boxshall, G. A., and Jaume, D. (2012). Centropages orsinii Giesbrecht, 1889 (Copepoda, Calanoida, Centropagidae) from an anchialine cave in Vanuatu. Zoosystema 34, 377–387. doi: 10.5252/z2012n2a11
Cabello-Yeves, P. J., Callieri, C., Picazo, A., Mehrshad, M., Haro-Moreno, J. M., Roda-Garcia, J. J., et al. (2021). The microbiome of the Black Sea water column analyzed by shotgun and genome centric metagenomics. Environ. Microbiome 16:5. doi: 10.1186/s40793-021-00374-1
Capella-Gutierrez, S., Silla-Martinez, J. M., and Gabaldon, T. (2009). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973. doi: 10.1093/bioinformatics/btp348
Ceccarelli, D. M., McKinnon, A. D., Andréfouët, S., Allain, V., Young, J., Gledhill, D. C., et al. (2013). “Chapter four - the Coral Sea: physical environment, ecosystem status and biodiversity assets” in Adv Mar Biol. ed. M. Lesser (66: Academic Press), 213–290.
Chen, C., Fu, L., Bi, N. S., Ge, R., Liu, G., Zhuang, Y., et al. (2018). Zooplankton community composition and diel vertical distribution in the Yongle blue hole, Xisha Islands, South China Sea. Oceanol. Limnol. Sin. 49, 594–603. doi: 10.11693/hyhz20171200328
Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560
Dahms, H.-U., and Qian, P.-Y. (2004). Life histories of the Harpacticoida (Copepoda, Crustacea): a comparison with meiofauna and macrofauna. J. Nat. Hist. 38, 1725–1734. doi: 10.1080/0022293031000156321
David, G. M., Moreira, D., Reboul, G., Annenkova, N. V., Galindo, L. J., Bertolino, P., et al. (2021). Environmental drivers of plankton protist communities along latitudinal and vertical gradients in the oldest and deepest freshwater lake. Environ. Microbiol. 23, 1436–1451. doi: 10.1111/1462-2920.15346
Di Giuseppe, G., and Dini, F. (2004). Identità dei dinoflagellati proliferanti in strutture sperimentali in situ nel Lago di Tovel (Nord Italia). Studi Trent. Sci. Nat. Acta Biol. 81, 427–438.
Duret, M. T., Pachiadaki, M. G., Stewart, F. J., Sarode, N., Christaki, U., Monchy, S., et al. (2015). Size-fractionated diversity of eukaryotic microbial communities in the eastern tropical North Pacific oxygen minimum zone. FEMS Microbiol. Ecol. 91:fiv037. doi: 10.1093/femsec/fiv037
Eddy, S. R. (2011). Accelerated profile HMM searches. PLoS Comput. Biol. 7:e1002195. doi: 10.1371/journal.pcbi.1002195
Edgcomb, V. P. (2016). Marine protist associations and environmental impacts across trophic levels in the twilight zone and below. Curr. Opin. Microbiol. 31, 169–175. doi: 10.1016/j.mib.2016.04.001
Edgcomb, V. P., Bernhard, J. M., Summons, R. E., Orsi, W., Beaudoin, D., and Visscher, P. T. (2014). Active eukaryotes in microbialites from Highborne cay, Bahamas, and Hamelin Pool (Shark Bay), Australia. ISME J. 8, 418–429. doi: 10.1038/ismej.2013.130
El-Rashidy, H., and Boxshall, G. A. (1999). Ergasilid copepods (Poecilostomatoida) from the gills of primitive Mugilidae (grey mullets). Syst. Parasitol. 42, 161–168. doi: 10.1023/A:1006075223683
Flegontova, O., Flegontov, P., Londoño, P. A. C., Walczowski, W., Šantić, D., Edgcomb, V. P., et al. (2020). Environmental determinants of the distribution of planktonic diplonemids and kinetoplastids in the oceans. Environ. Microbiol. 22, 4014–4031. doi: 10.1111/1462-2920.15190
Gao, J., Fu, T., Zhao, M., Zhang, H., and Tian, L. (2022). Preliminary study on formation process of Sansha Yongle blue hole. J. Trop. Oceanogr. 41, 171–183. doi: 10.11978/2021077
Gao, Z. M., Huang, J. M., Cui, G. J., Li, W. L., Li, J., Wei, Z. F., et al. (2019). In situ meta-omic insights into the community compositions and ecological roles of hadal microbes in the Mariana trench. Environ. Microbiol. 21, 4092–4108. doi: 10.1111/1462-2920.14759
Glockner, F. O., Yilmaz, P., Quast, C., Gerken, J., Beccati, A., Ciuprina, A., et al. (2017). 25 years of serving the community with ribosomal RNA gene reference databases and tools. J. Biotechnol. 261, 169–176. doi: 10.1016/j.jbiotec.2017.06.1198
Guillou, L., Bachar, D., Audic, S., Bass, D., Berney, C., Bittner, L., et al. (2013). The Protist ribosomal reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Res. 41, D597–D604. doi: 10.1093/nar/gks1160
Guillou, L., Viprey, M., Chambouvet, A., Welsh, R. M., and Worden, A. Z. (2010). Widespread occurrence and genetic diversity of marine parasitoids belonging to Syndiniales (Alveolata). Environ. Microbiol. 10, 3349–3365. doi: 10.1111/j.1462-2920.2008.01731.x
He, P., Xie, L., Zhang, X., Li, J., and Li, J. (2020). Microbial diversity and metabolic potential in the stratified Sansha Yongle blue hole in the South China Sea. Sci. Rep. 10:5949. doi: 10.1038/s41598-020-62411-2
Hodgson, L., and McDermid, K. (2000). Marine plants of Pohnpei and ant atoll: Chlorophyta. Phaeophyta Magnoliophyta 32, 289–308.
Hu, Y. O. O., Karlson, B., Charvet, S., and Andersson, A. F. (2016). Diversity of Pico- to Mesoplankton along the 2000 km salinity gradient of the Baltic Sea. Front. Microbiol. 7:679. doi: 10.3389/fmicb.2016.00679
Huang, Y., Gilna, P., and Li, W. (2009). Identification of ribosomal RNA genes in metagenomic fragments. Bioinformatics 25, 1338–1340. doi: 10.1093/bioinformatics/btp161
Ibarbalz, F. M., Henry, N., Brandão, M. C., Martini, S., Busseni, G., Byrne, H., et al. (2019). Global trends in marine plankton diversity across kingdoms of life. Cell 179, 1084–1097.e21. doi: 10.1016/j.cell.2019.10.008
Iliffe, T. M., and Kornicker, L. S. Worldwide diving discoveries of living fossil animals from the depths of anchialine and marine caves. Smithsonian Contributions to Marine Sciences (2009); 38, 269–280. Available at: http://www.tamug.edu/cavebiology/reprints/Reprint-195.pdf.
Jones, E., Devadatha, B., Abdel-Wahab, M. A., Dayarathne, M. C., and Cai, L. (2019). Phylogeny of new marine Dothideomycetes and Sordariomycetes from mangroves and deep-sea sediments. Bot. Mar. 63, 151–181. doi: 10.1515/bot-2019-0014
Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Kornicker, L. S., Iliffe, T. M., and Harrison-Nelson, E. (2007). Ostracoda (Myodocopa) from Anchialine caves and ocean blue holes. Zootaxa 1565, 1–151. doi: 10.11646/zootaxa.1565.1.1
Lejeusne, C., and Chevaldonné, P. (2010). Brooding crustaceans in a highly fragmented habitat: the genetic structure of Mediterranean marine cave-dwelling mysid populations. Mol. Ecol. 15, 4123–4140. doi: 10.1111/j.1365-294X.2006.03101.x
Li, T., Feng, A., Liu, Y., Li, Z., Guo, K., Jiang, W., et al. (2018). Three-dimensional (3D) morphology of Sansha Yongle blue hole in the South China Sea revealed by underwater remotely operated vehicle. Sci. Rep. 8:17122. doi: 10.1038/s41598-018-35220-x
Li, B., Liu, J., Zhou, S., Fu, L., and Zhang, X. H. (2020). Vertical variation in Vibrio community composition in Sansha Yongle blue hole and its ability to degrade macromolecules. J. Mar. Sci. Technol. 2, 60–72. doi: 10.1007/s42995-019-00003-4
Li, C., Love, G. D., Lyons, T. W., Fike, D. A., Sessions, A. L., and Chu, X. (2010). A stratified redox model for the Ediacaran Ocean. Science 328, 80–83. doi: 10.1126/science.1182369
Li, G., Wang, Y., Li, H., Zhang, X., and Gong, J. (2023). Quantifying relative contributions of biotic interactions to bacterial diversity and community assembly by using community characteristics of microbial eukaryotes. Ecol. Indic. 146:109841. doi: 10.1016/j.ecolind.2022.109841
Lindeque, P. K., Dimond, A., Harmer, R. A., Parry, H. E., Pemberton, K. L., and Fileman, E. S. (2015). Feeding selectivity of bivalve larvae on natural plankton assemblages in the Western English Channel. Mar. Biol. 162, 291–308. doi: 10.1007/s00227-014-2580-x
Liu, Y., He, H., Fu, L., Liu, Q., Yang, Z., and Zhen, Y. (2019). Environmental DNA sequencing reveals a highly complex eukaryote community in Sansha yongle blue hole, xisha, South China Sea. Microorganisms 7:624. doi: 10.3390/microorganisms7120624
Logares, R., Sunagawa, S., Salazar, G., Cornejo-Castillo, F. M., Ferrera, I., Sarmento, H., et al. (2014). Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. Environ. Microbiol. 16, 2659–2671. doi: 10.1111/1462-2920.12250
López-Escardó, D., Paps, J., de Vargas, C., Massana, R., Ruiz-Trillo, I., and del Campo, J. (2018). Metabarcoding analysis on European coastal samples reveals new molecular metazoan diversity. Sci. Rep. 8:9106. doi: 10.1038/s41598-018-27509-8
Lovejoy, C., Massana, R., and Pedrós-Alió, C. (2006). Diversity and distribution of marine microbial eukaryotes in the Arctic Ocean and adjacent seas. Appl. Environ. Microbiol. 72, 3085–3095. doi: 10.1128/AEM.72.5.3085-3095.2006
Luo, Z., Zhang, H., Li, Q., Wang, L., Mohamed, H. F., Lü, S., et al. (2022). Characterization of Amphidinium (Amphidiniales, Dinophyceae) species from the China Sea based on morphological, molecular, and pigment data. J. Oceanol. Limnol. 40, 1191–1219. doi: 10.1007/s00343-021-1049-2
Mentel, M., and Martin, W. (2010). Anaerobic animals from an ancient, anoxic ecological niche. BMC Biol. 8:32. doi: 10.1186/1741-7007-8-32
Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015
Mistry, J., Finn, R. D., Eddy, S. R., Bateman, A., and Punta, M. (2013). Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 41:e121. doi: 10.1093/nar/gkt263
Obiol, A., Giner, C. R., Sánchez, P., Duarte, C. M., Acinas, S. G., and Massana, R. (2020). A metagenomic assessment of microbial eukaryotic diversity in the global ocean. Mol. Ecol. Resour. 20, 718–731. doi: 10.1111/1755-0998.13147
Orsi, W., Charvet, S., Vd'ačný, P., Bernhard, J. M., and Edgcomb, V. P. (2012). Prevalence of partnerships between bacteria and ciliates in oxygen-depleted marine water columns. Front. Microbiol. 3:341. doi: 10.3389/fmicb.2012.00341
Pernice, M. C., Giner, C. R., Logares, R., Perera-Bel, J., Acinas, S. G., Duarte, C. M., et al. (2016). Large variability of bathypelagic microbial eukaryotic communities across the world’s oceans. ISME J. 10, 945–958. doi: 10.1038/ismej.2015.170
Pernthaler, J. J. N. R. M. (2005). Predation on prokaryotes in the water column and its ecological implications. Nat. Rev. Microbiol. 3, 537–546. doi: 10.1038/nrmicro1180
Porter, S. M., Agić, H., and Riedman, L. A. (2018). Anoxic ecosystems and early eukaryotes. Emerg. Top. Life Sci. 2, 299–309. doi: 10.1042/ETLS20170162
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Redmond, N. E., Raleigh, J., van Soest, R. W. M., Kelly, M., Travers, S. A. A., Bradshaw, B., et al. (2011). Phylogenetic relationships of the marine Haplosclerida (phylum Porifera) employing ribosomal (28S rRNA) and mitochondrial (cox1, nad1) gene sequence data. PLoS One 6:e24344. doi: 10.1371/journal.pone.0024344
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahe, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Roura, Á., Strugnell, J. M., Guerra, Á., González, Á. F., and Richardson, A. J. (2018). Small copepods could channel missing carbon through metazoan predation. Ecol. Evol. 8, 10868–10878. doi: 10.1002/ece3.4546
Salazar, G., Ruscheweyh, H.-J., Hildebrand, F., Acinas, S. G., and Sunagawa, S. (2021). mTAGs: taxonomic profiling using degenerate consensus reference sequences of ribosomal RNA genes. Bioinformatics (Oxford, England) 38, 270–272. doi: 10.1093/bioinformatics/btab465
Schoenle, A., Hohlfeld, M., Hermanns, K., Mahé, F., de Vargas, C., Nitsche, F., et al. (2021). High and specific diversity of protists in the deep-sea basins dominated by diplonemids, kinetoplastids, ciliates and foraminiferans. Commun. Biol. 4:501. doi: 10.1038/s42003-021-02012-5
Sherr, B. F., Sherr, E. B., Caron, D. A., Vaulot, D., and Worden, A. Z. (2007). Oceanic Protists. Oceanography 20, 130–134. doi: 10.5670/oceanog.2007.57
Shoemaker, K. M., Duhamel, S., and Moisander, P. H. (2019). Copepods promote bacterial community changes in surrounding seawater through farming and nutrient enrichment. Environ. Microbiol. 21, 3737–3750. doi: 10.1111/1462-2920.14723
Škodová-Sveráková, I., Záhonová, K., Juricová, V., Danchenko, M., Moos, M., Baráth, P., et al. (2021). Highly flexible metabolism of the marine euglenozoan protist Diplonema papillatum. BMC Biol. 19:251. doi: 10.1186/s12915-021-01186-y
Steinberg, D. K., and Landry, M. R. (2017). Zooplankton and the ocean carbon cycle. Annu. Rev. Mar. Sci. 9, 413–444. doi: 10.1146/annurev-marine-010814-015924
Stoeck, T., Taylor, G. T., and Epstein, S. S. (2003). Novel eukaryotes from the permanently anoxic Cariaco Basin (Caribbean Sea). Appl. Environ. Microbiol. 69, 5656–5663. doi: 10.1128/AEM.69.9.5656-5663.2003
Suárez-Morales, E. (2015). “Class maxillopoda” in Thorp and Covich's freshwater invertebrates. ed. J. H. Thorp and D. C. Roger (Elsevier), p. 709–755.
Sugai, Y., Natori, N., Tsuchiya, K., Nakagawa, M., Honda, M. C., Shimode, S., et al. (2023). Ingestion rate estimated from food concentration and predatory role of copepod nauplii in the microbial food web of temperate embayment waters. J. Plankton Res. 45, 325–337. doi: 10.1093/plankt/fbad002
Takahashi, M., Fujii, K., and Parsons, T. R. (1973). Simulation study of phytoplankton photosynthesis and growth in the Fraser River estuary. Mar. Biol. 19, 102–116. doi: 10.1007/BF00353581
Torres-Beltran, M., Mueller, A., Scofield, M., Pachiadaki, M. G., Taylor, C., Tyshchenko, K., et al. (2019). Sampling and processing methods impact microbial community structure and potential activity in a seasonally Anoxic Fjord: Saanich inlet, British Columbia. Front. Mar. Sci. 6:6. doi: 10.3389/fmars.2019.00132
Wang, Y., Gao, Z. M., Li, J., He, L. S., Cui, G. J., Li, W. L., et al. (2019). Hadal water sampling by microbial filtration and fixation (ISMIFF) apparatus. Deep Sea Res. Part I Oceanogr. Res. Pap. 144, 132–137. doi: 10.1016/j.dsr.2019.01.009
Wang, Y., Zhang, W. P., Cao, H. L., Shek, C. S., Tian, R. M., Wong, Y. H., et al. (2014). Diversity and distribution of eukaryotic microbes in and around a brine pool adjacent to the Thuwal cold seeps in the Red Sea. Front. Microbiol. 5:37. doi: 10.3389/fmicb.2014.00037
Weber, F., Anderson, R., Mylnikov, A., and Jürgens, K. (2014). Morphological and molecular approaches reveal highly stratified protist communities in Baltic Sea pelagic redox gradients. Aquat. Microb. Ecol. 73, 1–16. doi: 10.3354/ame01702
Wei, Z.-F. L. W.-L., Li, J., Chen, J., Xin, Y.-Z., He, L.-S., and Wang, Y. (2020). Multiple in situ nucleic acid collections (MISNAC) from Deep-Sea waters. Front. Mar. Sci. 7:81. doi: 10.3389/fmars.2020.00081
Wylezich, C., and Jürgens, K. (2011). Protist diversity in suboxic and sulfidic waters of the Black Sea. Environ. Microbiol. 13, 2939–2956. doi: 10.1111/j.1462-2920.2011.02569.x
Xie, L. P., Wang, B. D., Pu, X. M., Xin, M., He, P. Q., Li, C. X., et al. (2019). Hydrochemical properties and chemocline of the Sansha Yongle blue hole in the South China Sea. Sci. Total Environ. 649, 1281–1292. doi: 10.1016/j.scitotenv.2018.08.333
Xu, H., Luo, X., Qian, J., Pang, X., Song, J., Qian, G., et al. (2012). FastUniq: a fast de novo duplicates removal tool for paired short reads. PLoS One 7:e52249. doi: 10.1371/journal.pone.0052249
Yao, P., Chen, L., Fu, L., Yang, Z., Bi, N., Wang, L., et al. (2018). Controls on vertical nutrient distributions in the Sansha Yongle blue hole, South China Sea. Chin. Sci. Bull. 63, 2393–2402. doi: 10.1360/N972018-00155
Yao, P., Wang, X. C., Bianchi, T. S., Yang, Z. S., Fu, L., Zhang, X. H., et al. (2020). Carbon cycling in the World's deepest blue hole. J. Geophys. Res. Biogeosci. 125:e2019JG005307. doi: 10.1029/2019JG005307
Yeh, Y. C., McNichol, J., Needham, D. M., Fichot, E. B., Berdjeb, L., and Fuhrman, J. A. (2021). Comprehensive single-PCR 16S and 18S rRNA community analysis validated with mock communities, and estimation of sequencing bias against 18S. Environ. Microbiol. 23, 3240–3250. doi: 10.1111/1462-2920.15553
Zhang, H., Wei, T., Li, Q., Fu, L., He, L., and Wang, Y. (2023). Metagenomic 16S rDNA reads of in situ preserved samples revealed microbial communities in the Yongle blue hole. PeerJ 11, e16257–e16218. doi: 10.7717/peerj.16257
Zhang, H., Wei, T., Zhang, J., Li, Q., Fu, L., He, L., et al. (2024). Microbial mediation of cryptic methane cycle in semiclosed marine water column. TIG 2024:100082. doi: 10.59717/j.xinn-geo.2024.100082
Zhou, Y. L., Mara, P., Cui, G. J., Edgcomb, V. P., and Wang, Y. (2022). Microbiomes in the challenger deep slope and bottom-axis sediments. Nat. Commun. 13:1515. doi: 10.1038/s41467-022-29144-4
Keywords: Yongle blue hole, 18S rRNA gene, eukaryotic communities, anoxic water, metagenomics
Citation: Zhang H, Wei T, Li Q, Fu L, Li M, He L and Wang Y (2024) Metagenomic 18S rDNA reads revealed zonation of eukaryotic communities in the Yongle blue hole. Front. Microbiol. 15:1420899. doi: 10.3389/fmicb.2024.1420899
Edited by:
Memory Tekere, University of South Africa, South AfricaReviewed by:
Ilunga Kamika, University of South Africa, South AfricaKedibone Masenya, Agricultural Research Council of South Africa (ARC-SA), South Africa
Grace Nkechinyere Ijoma, University of South Africa, South Africa
Copyright © 2024 Zhang, Wei, Li, Fu, Li, He and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Wang, d2FuZ3lvbmdAc3oudHNpbmdodWEuZWR1LmNu