Journal of Oceanology and Limnology   2021, Vol. 39 issue(5): 1718-1729     PDF
Institute of Oceanology, Chinese Academy of Sciences

Article Information

SHI Junfeng, LEI Yanli, LI Haotian, LI Tiegang
NGS-metabarcoding revealing novel foraminiferal diversity in the Western Pacific Magellan Seamount sediments
Journal of Oceanology and Limnology, 39(5): 1718-1729

Article History

Received Nov. 11, 2020
accepted in principle Jan. 14, 2021
accepted for publication Aug. 16, 2021
NGS-metabarcoding revealing novel foraminiferal diversity in the Western Pacific Magellan Seamount sediments
Junfeng SHI1,3,7, Yanli LEI1,4,5,6,7, Haotian LI1,6,7, Tiegang LI2     
1 Laboratory of Marine Organism Taxonomy and Phylogeny, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China;
2 Key Laboratory of Marine Sedimentology and Environmental Geology, First Institute of Oceanography, Ministry of Natural Resources (MNR), Qingdao 266061, China;
3 Clinical Research Center, Affiliated Hospital of Weifang Medical University, Weifang 261031, China;
4 Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China;
5 Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China;
6 University of Chinese Academy of Sciences, Beijing 100049, China;
7 Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
Abstract: Seamount is a unique deep-sea ecosystem widely distributed in the world. Its biodiversity is vibrant due to its specific geographical and hydrological conditions. However, the diversity and features of foraminifera in such an environment have rarely been studied. We extracted environmental DNA (eDNA) in sediments and amplified the partial small subunit ribosomal DNA (SSU rDNA) of foraminifera to understand the foraminiferal diversity from four sites in Magellan Seamount (Western Pacific Ocean). Partial SSU rDNA sequencing was conducted and 912 979 foraminiferal reads were obtained and gathered into 266 operational taxonomic units (OTUs). In the available dataset, a high proportion of rare OTUs and low identity OTUs in each studied sample showed that the Magellan Seamount foraminiferal community might have a high genetic novelty. The relative abundance of foraminifera varied between replicates probably due to the existence of bias in amplification process and patchiness of the deep-sea floor. It showed that the Magellan Seamount has a relatively high benthic foraminiferal diversity characterized mainly by monothalamiids (76.37% of total reads) in association with rotaliids (19.03% of total reads), including planktic foraminiferal sequences (38.58% of rotaliids; 7.36% of total reads). The remaining reads were assigned to miliolids (0.83% of total reads) and textulariids (0.66% of total reads), and 3.11% of total reads are unassigned to a specific family. The comparative analysis with foraminiferal assemblages from coastal and deep-sea environments indicated that seamounts could aggregate species from a nearby deep-sea.
Keywords: metabarcoding    Magellan Seamount    foraminifera    Western Pacific Ocean    

Foraminifera are commonly defined as marine testate protists among the most abundant, diverse, and widely distributed deep-sea microfauna. As significant components of deep-sea biotic communities, foraminifera play pivotal roles in ecological processes and the flow of energy and materials through ecosystems (Cornelius and Gooday, 2004; Gooday and Jorissen, 2012; Lei et al., 2014). Benthic foraminifera are sensitive to biotic and abiotic factors, and their abundance and test morphology allowed the recovery of paleoecological changes (Mojtahid et al., 2010; Bejaoui et al., 2019). Their assemblages can reveal information about salinity, temperature, food supply, oxygen, and bathymetry (Lei et al., 2017). Planktic foraminifera are considered as good evaluating indicators for open marine environments (Li and Keller, 1998; Saïdi and Zaghbib-Turki, 2016). Their community structures tend to be affected by abiotic factors and indicate the water hydrography (Vincent and Berger, 1981; Zhang, 1985; Naidu and Malmgren, 1996; Li and Keller, 1998; Saïdi and Zaghbib-Turki, 2016; Jian et al., 2020). These reconstructions are established from the chemical composition and morphology of foraminiferal fossil test (Li and Keller, 1998; Katz et al., 2010).

Furthermore, the most part of living foraminifera ecological studies focused on shallow water environments (Lei et al., 2014, 2017). Little is known about the biology and the diversity of living benthic foraminifera at the deep-sea environments as few works have been carried out on them (Lee and Anderson, 1991; Lecroq et al., 2011; André et al., 2014, Angeles et al., 2020), especially in the seamounts where research data on foraminifera are very few (Ohkushi and Natori, 2001; Shi et al., 2020). In addition, deep-sea sampling is challenging due to the vastness and remoteness of the deep seabed. Insufficient sampling limits to obtain a comprehensive assessment of foraminiferal diversity in the deep-sea (Lecroq et al., 2011; Goineau and Gooday, 2019).

A seamount corresponds to an uplifted terrain with a height of more than 1 000 m (without reaching the sea surface) in the deep-sea (Epp and Smoot, 1989). It is a unique habitat in the deep-sea, which is globally distributed in the oceans (Clark et al., 2010; Yesson et al., 2011). The seamount has a high biodiversity due to its specific geographical and hydrological characteristics (Taylor Column) (Parker and Tunnicliffe, 1994; Stocks, 2004a, b). All studies on foraminiferal diversity in seamount are based on their morphology (Nienstedt and Arnold, 1988; Heinz et al., 2004; Gooday and Jorissen, 2012). However, most deep-sea species are tiny, fragile, and challenging to identify (Lecroq et al., 2011). A large spectrum of deep-sea foraminifera belongs to early lineages characterized by simple, single-chambered (monothalamiids), organic-walled, or agglutinated tests and are poorly preserved in the fossil record (Gooday et al., 2001, 2004; Lecroq et al., 2011). Thus, it is challenging and even impossible to decipher living foraminiferal diversity in seamounts based only on test morphology. In the past decade, molecular studies have revealed the diversity and geographic distribution of benthic foraminifera in deep-sea sediments (Pawlowski and Holzmann, 2007; Lecroq et al., 2011; Lecroq, 2014). However, most of these studies relied on clone libraries, limiting the number of sequences available for analysis (Lecroq et al., 2011). Next-generation sequencing (NGS) used for environmental rDNA sequencing reduced clone libraries' limitations and made the decoding application possible (Lecroq et al., 2011; Frontalini et al., 2020). This approach made it possible to analyze the living foraminiferal diversity that is tiny, fragile, or rare to identify using classic methods in the deep-sea comprising seamount habitats (Logares et al., 2014; Egge et al., 2015; Piredda et al., 2018; Shi et al., 2020).

Meanwhile, the total environmental DNA (eDNA) contains living cell and extracellular DNA (Pietramellara et al., 2009), it is characterized by mixing the genomic DNA of different biota present in the studied environment and has enormous potential to address fundamental research questions related to molecular biology, ecology, paleontology, and environmental sciences (Thomsen and Willerslev, 2015). This provided an opportunity to study the sequence of foraminifera ex situ.

This study aims to (1) identify the taxonomic composition and diversity of foraminifera in the Magellan Seamount (Western Pacific Ocean) sediments, (2) reveal variations between different replications or sites, and (3) determine if there are foraminiferal sequences from other sources in sediments except in-situ benthic foraminiferal data.

2 MATERIAL AND METHOD 2.1 Location and sample collection

Twelve samples have been collected from four sites (with three sample replicates from each site) from the Magellan Seamount, which lies adjacent to the southeast of the Mariana Trench in the Western Pacific Ocean and comprises two concomitant top flat seamounts (Fig. 1). To compare foraminiferal diversity from the four sites, sediments were collected at water depths ranging from 1 308 to 1 968 m using a push-core operated by a remotely operated vehicle (ROV), Faxian (Discovery in Chinese), equipped by the R/V Kexue (Science in Chinese) in April 2018 (Fig. 1). The ROV was only launched four times in four chosen sites as it was challenging to obtain the sediment samples from the Magellan Seamount. After collection, sediments were immediately frozen at -80 ℃ (sample information is presented in Table 1).

Fig.1 Location map of the studied sites of Magellan Seamount in the Western Pacific Ocean a. location map of the Magellan Seamount in the Western Pacific Ocean; b. detailed overview of sampling sites, each of the red circles represent geographic location of the samples.
Table 1 Summary of sampling sites, numbers of foraminiferal read and OTUs
2.2 eDNA extraction, amplification, and sequencing

The eDNA was extracted using DNeasy PowerSoil kit (QIAGEN, Germany) for sediment samples. Approximately 0.25 g of the sediment was required for each sample under the extraction process. Afterward, the eDNA extracts were stored at -20 ℃. The foraminiferal 37f hypervariable region of SSU rDNA was amplified using a High-Fidelity Master Mix (TRANSGEN, China) by nested polymerase chain reaction (PCR). We used foraminiferal-specific primer pair, s14F3 (5'-ACGCAMGTGTGAAACTTG-3') and s17 (5'-CGGTCACGTTCGTTGC-3'), and a long fragment (300–500 bp) was amplified in the first reaction. A shorter fragment (210–510 bp) was re-amplified using the s14F1 (5'-AAGGGCACCACA-AGAACGC-3') and s17 primer pairs following Lejzerowicz et al. (2010). In the first amplification, 34 cycles were programmed at 50 ℃ for annealing temperature; in the second amplification, the program was set at 24 cycles and 52 ℃. The PCR system contained 12.5 μL of two times High-Fidelity Master Mix, 0.5 μL of each primer (10 μmol/L), 2 μL of eDNA template, and distillation-distillation H2O (ddH2O) was added to reach a final volume of 25 μL.

The PCR products were purified using the Gel Extraction Kit (OMEGA, USA). Sequencing was performed in Novogene Bioinformatics Technology Co., Ltd. The sequence libraries were generated using the TruSeq® DNA PCR-Free Sample Preparation Kit for Illumina (NEB, USA) following the manufacturer's recommendations. The libraries' quality was conducted using the Agilent Bioanalyzer 2100 system. The library was sequenced using an Illumina HiSeq2500 platform, and 250-bp paired-end reads were generated.

2.3 Data analysis methods

Paired-end reads (raw reads) were split into samples based on unique barcodes and merged using FLASH (V1.2.7,; Magoč and Salzberg, 2011). FastQC (, an adopted application in quality control of NGS following Brown et al. (2017), was used to assess this dataset's read quality. Low-quality reads (quality value of < 20) were removed according to the QIIME (V 1.9.1, (Caporaso et al., 2010; Bokulich et al., 2013) after filtering, and singletons were excluded. OTUs that occurred in less than three samples were excluded to avoid bias due to the presence of extracellular DNA. As PCR errors often generated rare sequences, reads occurring in more than three samples and showing a minimal abundance of ten in all samples are obtained.

The filtered reads (clean reads) were generated to operational taxonomic units (OTUs) using the Unoise3 (developed by Edgar (2016)) pipeline with an identified threshold of 100%. After the Unoise3 pipeline generation, each OTU—an artificially defined taxonomic unit based on sequence similarity used to distinguish taxonomic groups—represented a unique sequence. OTUs with less than 10 reads occurring in less than three samples were excluded from the analysis. Representative sequences were extracted from the OTU table and taxonomically assigned based on the PR2 database (Protist Ribosomal Reference Database; Guillou et al., 2012) using BLAST (Version 2.7.1); non-foraminifera (Rank 3 in PR2 database) OTUs were removed. OTUs were defined as abundant when the number of sequences exceeds 1% of all sequences in a sample. Rare OTUs were defined when the number of sequences was less than 0.1% of all sequences in a sample, as was the case for Pedrós-Alió (2006) and Fuhrman (2009). α-diversity (within-habitat diversity) can reflect the richness and evenness of species within habitat. Thus, the Simpson, Shannon, and Margalef α-diversity indices were calculated in this work using QIIME v.1.91.

3 RESULT 3.1 Overview of NGS dataset

The foraminiferal SSU rDNA dataset contained 1 059 636 raw reads. Based on the Illumina sequencing platform, an average of 98.20% reads in the analyzed samples had a high quality score (≥30). The preserved reads were compared with the PR2 database, and seven unique sequences with 1 516 reads were assigned to Excavata (amitochondriate, Protozoa without mitochondria, Rank 2 in PR2 database). These non-foraminiferal sequences were discarded from the dataset. After quality filtering, 912 979 foraminiferal reads were retained and the percentage of filtered reads varied from 4.25% (FXD174.1) to 22.20% (FXD175.3).

The rarefaction curves (Fig. 2a) show that all samples were fully sequenced as they tended to saturation. In this study, the numbers of effective reads in all samples were more than 60 000 (see Reads after filtered; Table 1). The reads were clustered into 266 distinct OTUs based on 100% sequence similarity. To avoid the threshold affecting clustering accuracy and achieve the best phenotype resolution, we defined each sequence as a separate OTU. The OTUs for each sample ranged from 107 to 176, with an average of 139. The details of the samples and foraminiferal reads used in our analysis are shown in Table 1.

Fig.2 Comparison of α-diversity from the four studied sites a. the rarefaction curves. The curves of 4 sites were integrated from 12 samples; b. the histogram showed α-diversity indices of 4 sites.
3.2 Taxonomic composition

The taxonomic assignment of the obtained 266 OTUs revealed 38 genera and 45 species among planktic and benthic foraminifers. The OTUs were assigned to a reference database with low similarity in most cases, and only 91 OTUs (34.21%, corresponding to 20.94% of the reads) reached a similarity of more than 97% (identity-value > 0.97, see Supplementary Table S1) with targeted sequences in the PR2 database. The average identity-value of all detected OTUs was 93.59. The Sankey diagram (Fig. 3a) showed the flow of NGS reads (left) and OTUs (right) from the four studied locations, which are assigned to different foraminiferal families (Rank 5 in PR2 database). Approximately, half of the OTUs (119 out of 266) were with less than 1 000 reads (Fig. 3b & c). An OTUs with maximum numbers of reads were assigned to Toxisarcon sp. (74 346 reads), Elphidium frigidum (33 875 reads), Siphonaperta sp. (7 524 reads), and Trochammia hadai (2 318 reads), becoming the most abundant species of monothalamiids, rotaliids, miliolids, and textulariids, respectively. The 15 most abundant OTUs approximately reached half of the reads. Only 13 OTUs occurred in all the studied samples.

Fig.3 The relationship between OTUs and reads a. next-generation sequencing (NGS) reads (left), classification (middle), and OTUs abundance (right). The violin plots show the number of OTUs in each replication (b) and site (c).

Figure 4a shows the distribution of foraminifera in the analyzed samples, including all reads in sediments. Most OTUs were assigned to (Rank 5 in PR2 database) monothalamiids (single-chambered and soft-walled) (191 out of 266), and 76.37% of the OTU reads were distributed between the samples under study (697 226 out of 912 979). These monothalamiids correspond to more than 70% of the reads from eight samples (Fig. 4b). Six OTUs with 28 392 reads were unassigned to a definite family (Rank 5 named foraminifera_XX; Supplementary Table S1), which represented 3.11% of all reads. The remaining 69 OTUs were assigned to rotaliids (58 OTUs), miliolids (2 OTUs), and textulariids (9 OTUs), representing 19.03%, 0.83%, and 0.66% of the foraminiferal reads (including planktic and benthic foraminifera), respectively.

Fig.4 Relative abundances of the number of OTUs (a) and sequences (b) related to the major families of foraminifera detected by DNA sequencing in different sites

With calcareous tests, planktic foraminifera can be placed within the radiation of the rotaliids (Ujiié et al., 2008). In this dataset, 23 OTUs (out of 266) were assigned to planktic foraminifera (globigerinacea, Rank 6 in PR2 database), representing 67 239 reads (7.36% of the total dataset; Supplementary Table S1). Planktic foraminiferal sequences were annotated to six species: Globigerinita uvula, Globigerinita glutinata, Globigerina bulloides, Candeina nitida, Pulleniatina obliquiloculata, and Globigerinella siphonifera. Among these species, Globigerinita uvula had 14 OTUs and the most number of reads (53 346), representing approximately 79.34% of the planktic foraminiferal sequences. The number of OTUs ranged 14 (in FXD178) – 23 (in FXD175), representing 7.25% (in FXD178) – 10.31% (in FXD175) of each site's total OTUs (Fig. 5a). In the four sites, 4 296 (in FXD172) – 27 132 (in FXD175) planktic foraminiferal reads were retained, representing 1.99% (in FXD172) – 13.32% (in FXD175) of the total sample sequence (Fig. 5b). The absolute and relative abundance of planktic reads varied from site to site, but the number and proportion of OTUs were stable (Fig. 5).

Fig.5 Occurrences of planktic foraminiferal eDNA in the abyssal sediments a. the number and proportion of planktic foraminiferal OTUs in each studied site; b. the number and proportion of planktic foraminiferal reads and their proportions in each studied site.
3.3 Comparative analyses

The benthic foraminiferal diversity variation in sediments was observed among different sites and replicates. The number of OTUs among sites was the highest in FXD172 (on average of 164.7) and the lowest in FXD178 (on average of 128.7; Table 1). The rarefaction curves (Fig. 2a) show that FXD174 has the highest value (Margalef value) and number of observed OTUs but the lowest Shannon-Wiener value. Furthermore, the values of the Margalef and Simpson indices have no remarkable changes among the four stations. The number of shared OTUs in all sites reached 131. Only 26.6% of all OTUs occurred in three replicates. Besides, the proportion of these OTUs in each site was 39.1% (FXD172), 16.5% (FXD174), 18.8% (FXD175), and 40% (FXD178) (Fig. 6ad). The average percentage of the shared OTUs in a site was about 32.8% (46.2% in FXD172, 18.6% in FXD174, 22.4% in FXD175, and 44% in FXD178). The number of exclusive OTUs in FXD174.3 and FXD175.3 exceeded the shared ones in all samples of one site (Fig. 6bc). To mine the dataset more deeply, we divided OTUs into abundant and rare. Moreover, statistical work was conducted separately in two families (Table 1). The number of abundant and rare OTUs varied between sites, whereas the difference between replicates among each site seems to be more significant.

Fig.6 The number of unique and shared OTUs in different replications and sites a. FXD172; b. FXD174; c. FXD175; d. FXD178; e. the number of unique and shared OTUs in the different studied sites.
4 DISCUSSION 4.1 Deficiency of metabarcoding in foraminiferal research

In theory, the more sequences in a dataset, the better it reflects the sample's actual composition. However, superfluous sequences meant a greater investment in the sequencing process, so we had to compromise between the sequencing data amount and sequencing costs. Approximately 20 000 effective reads per sample seemed sufficient for a diversity study (Kozich et al., 2013), covering more than 95% of the OTUs (Lecroq et al., 2011). There is no rationale for defining the optimal size of environmental sequence datasets (Ni et al., 2013; Lejzerowicz et al., 2014). In this dataset, all the studied samples have more than 60 000 effective reads, which means the data are sufficient for ecological research.

The PCR amplification process is biased but it influences the estimation of species relative abundances (Acinas et al., 2005; Clarke et al., 2014; Pawluczyk et al., 2015). Predictions of the species relative abundances often depend on the number of sequences. However, equal amounts of template DNA from different species did not yield equal proportions of amplicons (Elbrecht and Leese, 2015; present study). The proportion of different families (including planktic and benthic foraminifera) among samples was similar regarding the number of OTUs (Fig. 4a). In contrast, the proportion of different families among samples had extremely high variations regarding the relative abundance of sequences (Fig. 4). This indicates that the relative abundance of the sequence was inaccurate for biomass evaluation. The efficiency's inequality can explain the result in amplifying multi-species mixed templates using the universal primer (Polz and Cavanaugh, 2008; Piñol et al., 2015).

4.2 High genetic novelty of foraminifera in Magellan Seamount

The low identity-values of OTUs in our dataset may show a high genetic novelty as it was reported in other seamounts (Shi et al., 2020). The application of NGS contributed substantially to the detection of broader protist diversity. However, related research is also facing some challenges (Pawlowski et al., 2014). Many researchers have used a 97% sequence identity threshold to distinguish different species in the previous works (Ciardo et al., 2006; Westcott and Schloss, 2017). This threshold was first proposed by Stackebrandt and Goebel (1994) when 16S rDNA sequences were considered suitable to estimate the species composition of a community. However, Edgar (2018) proved that the optimal identity threshold of a high-quality set of 16S rRNA sequences from known species was higher than 97%, and the optimal threshold was data-dependent. Previous studies have indicated that foraminifera have a high rDNA variation, leading to overestimated molecular diversity in environmental samples. Meanwhile, the database's lack of reference sequences may also cause the pseudomorph of high genetic novelty. The extremely high genetic novelty recorded in our dataset may be a pseudomorph generated due to the foregoing circumstances or the embodiment of high biodiversity in the seamount area.

4.3 Variations of benthic foraminifera among samples and replicates

Variations of benthic foraminiferal relative abundances occurred in the four studied sites but were the most evident in the site FXD178 (Fig. 4). The water depth in the four stations is approximately about 1 300–2 000 m and cannot be regarded as a controlling factor of these abundances' changes. However, our results infer that heterogeneity of surface sediments can be responsible for the recorded variations in benthic foraminiferal relative abundances (Frontalini et al., 2020; Shi et al., 2020). FXD178 was in a relatively different habitat compared with other stations, and these habitat variations led to even more-significant differences in diversity than that of water depths.

Not all significant families appeared in the three replications (related to the four sampled sites). For example, miliolids were absent in FXD172.2 and FXD172.3 but have been recorded in FXD172.1 (Fig. 4a). With a considerable number of rare OTUs (Table 1) in each sample, heterogeneity occurred in almost every replication. For example, FXD174.3 and FXD175.3 have more exclusive OTUs than the other two replications. This heterogeneity among replications can be explained by the surface sediments "patchiness", which was characterized by vertical and horizontal distribution variations of microfauna in microhabitats even occurred on the centimeter scale and was widespread in the deep-sea (Jorissen et al., 2007; Lejzerowicz et al., 2014; Bejaoui et al., 2019). The heterogeneity of deep-sea microhabitats affected the diversity of replications, and this phenomenon appeared in the study of nematodes (Danovaro et al., 2013), ciliate (Zhao et al., 2017), and foraminifera (Lejzerowicz et al., 2014).

4.4 Planktic foraminiferal assemblages in sediments

In accordance with Morard et al. (2017) and Shi et al. (2020), we infer that the DNA of planktic foraminifera was preserved in the sediments. For example, the planktic foraminiferal reads accounted for 1.94% of the total reads in the seamount (Shi et al., 2020). An organism's DNA can be adsorbed to solid particles (Herndl and Reinthaler, 2013), such as marine snow, fecal pellets, and massive falls of phytoplankton material after surface blooms (Torti et al., 2015). The solid particles may agglomerate particles and enable the absorbed DNA to sink rapidly and become a record of extracellular DNA (Corinaldesi et al., 2011; Torti et al., 2015; Morard et al., 2017). The extracellular DNA could reconstruct planktic eukaryotic groups (Capo et al., 2015) and become a significant archive of present eukaryotic gene sequences (Corinaldesi et al., 2011). Planktic foraminiferal sequences in sediments provide an opportunity to understand the upper water column's sedimentation process. Although the number of planktic sequences fluctuated considerably between different sites or samples, the number of OTUs kept the same order of magnitude. This enables to reconstruct species composition in the water column and determine the traceability of planktic sequences in sediments (Angeles et al., 2020). According to this track, we can determine the dynamic sedimentation pattern of plankton undercurrent. This will provide a basis for using sediment data to reconstruct the water environment column.

4.5 Characteristics of foraminiferal communities

The dominance of monothalamiids is consistent with previously published results based on morphology (Gooday, 2002; Gooday et al., 2004) and molecular biology (Lecroq et al., 2011; Lejzerowicz et al., 2014; Shi et al., 2020). Some classifications inferred from eDNA data (i.e., monothalamiids-clade-E in Supplementary Table S1) are not the most abundant family in our dataset, different from results of Lejzerowicz et al. (2014). The composition of foraminiferal assemblages in the seamounts had significant proximity to the ordinary deep-sea and coastal area characterized by abundant monothalamiids and rotaliids taxa as recorded by Lecroq et al. (2011), Pawlowski et al. (2011), and Li et al. (2020).

As for the diversity of foraminifera, Magellan Seamount had a relatively high Shannon-Wiener index (> 3; Fig. 2b) compared with adjacent abyssal plains (< 2.6; Shi et al., 2020). In the typical deep-sea, the Shannon-Wiener index of foraminifera (including planktic and benthic) exhibited a lower value (approximately 2; Corliss et al., 2009). In the coastal Jiaozhou Bay, the Shannon-Wiener index of foraminifera was higher (superior to 5; Li et al., 2020) than the Shannon-Wiener index recorded in the studied Magellan Seamount. From the coastal shallow water area to the ordinary deep-sea and then to the seamounts, the diversity of foraminifera first declined and then rose. Nevertheless, the composition of foraminiferal assemblages did not change and still be dominated by monothalamiids and rotaliids. The foraminiferal communities in the seamounts seem more similar to the distant coastal area than the nearby deep-sea and seamount sections. Furthermore, the specific hydrologic conditions (Taylor Column) (Parker and Tunnicliffe, 1994) of the seamounts favored species gathering. Thus, we speculate that the foraminiferal community in the deep-sea was from the foraminiferal species diffused from the coastal area. The foraminiferal planktic larvae from the deep-sea aggregated around the seamounts, which caused the high foraminiferal diversity in seamounts and similar species composition to the deep-sea and coastal areas.


Magellan Seamount may have a high genetic novelty due to the high proportion of rare and low identity OTUs. The relative abundance of foraminiferal families varied between replicates. The planktic foraminiferal sequences were detected in the sedimentary environment DNA dataset. The monothalamiids and rotaliids were the prominent foraminiferal families in the Magellan Seamount.

After a comparative analysis of the foraminiferal communities in the coastal area and deep-sea, we speculate that seamounts could aggregate species.


Sequence data that support the findings of this study were submitted to the Sequence Read Archive under accession number PRJNA675069.


The data and samples were collected onboard R/V Kexue. We are also very grateful to Dr. Ezzedine SAÏDI for his useful comments and suggestions.

Electronic supplementary material

Supplementary material (Supplementary Table S1) is available in the online version of this article at

Acinas S G, Sarma-Rupavtarm R, Klepac-Ceraj V, Polz M F. 2005. PCR-induced sequence artifacts and bias: insights from comparison of two 16S rRNA clone libraries constructed from the same sample. Applied and Environmental Microbiology, 71(12): 8 966-8 969. DOI:10.1128/AEM.71.12.8966-8969.2005
André A, Quillévéré F, Morard R, Ujiié Y, Escarguel G, De Vargas C, De Garidel-Thoron T D, Douady C J. 2014. SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications. PLoS One, 9(8): e104641. DOI:10.1371/journal.pone.0104641
Angeles I B, Lejzerowicz F, Cordier T, Scheplitz J, Kucera M, Ariztegui D, Pawlowski J, Morard R. 2020. Planktonic foraminifera eDNA signature deposited on the seafloor remains preserved after burial in marine sediments. Scientific Reports, 10: 20351. DOI:10.1038/s41598-020-77179-8
Bejaoui A, Saïdi E, Zaghbib-Turki D. 2019. Small benthic foraminiferal biostratigraphy and palaeoecology during the Campanian Maastrichtian transition in north-western Tunisia. Turkish Journal of Earth Sciences, 28(4): 500-530. DOI:10.3906/yer-1807-49
Bokulich N A, Subramanian S, Faith J J, Gevers D, Gordon J I, Knight R, Mills D A, Caporaso J G. 2013. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nature Methods, 10(1): 57-59. DOI:10.1038/nmeth.2276
Brown J, Pirrung M, Mccue L A. 2017. FQC dashboard: integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics, 33(19): 3 137-3 139. DOI:10.1093/bioinformatics/btx373
Capo E, Debroas D, Arnaud F, Domaizon I. 2015. Is planktonic diversity well recorded in sedimentary DNA? Toward the reconstruction of past protistan diversity. Microbial Ecology, 70(4): 865-875. DOI:10.1007/s00248-015-0627-2
Caporaso J G, Kuczynski J, Stombaugh J, Bittinger K, Bushman F D, Costello E K, Fierer N, Peña A G, Goodrich J K, Gordon J I, Huttley G A, Kelley S T, Knights D, Koenig J E, Ley R E, Lozupone C A, McDonald D, Muegge B D, Pirrung M, Reeder J, Sevinsky J R, Turnbaugh P J, Walters W A, Widmann J, Yatsunenko T, Zaneveld J, Knight R. 2010. QⅡME allows analysis of high-throughput community sequencing data. Nature Methods, 7(5): 335-336. DOI:10.1038/nmeth.f.303
Ciardo D E, Schär G, Böttger E C, Altwegg M, Bosshard P P. 2006. Internal transcribed spacer sequencing versus biochemical profiling for identification of medically important yeasts. Journal of Clinical Microbiology, 44(1): 77-84. DOI:10.1128/JCM.44.1.77-84.2006
Clark M R, Rowden A A, Schlacher T, Williams A, Consalvey M, Stocks K I, Rogers A D, O'Hara T D, White M, Shank T M, Hall-Spencer J M. 2010. The ecology of seamounts: structure, function, and human impacts. Annual Review of Marine Science, 2(1): 253-278. DOI:10.1146/annurev-marine-120308-081109
Clarke L J, Soubrier J, Weyrich L S, Cooper A. 2014. Environmental metabarcodes for insects: in silico PCR reveals potential for taxonomic bias. Molecular Ecology Resources, 14(6): 1 160-1 170. DOI:10.1111/1755-0998.12265
Corinaldesi C, Barucca M, Luna G M, Dell'Anno A. 2011. Preservation, origin and genetic imprint of extracellular DNA in permanently anoxic deep-sea sediments. Molecular Ecology, 20(3): 642-654. DOI:10.1111/j.1365-294X.2010.04958.x
Corliss B H, Brown C W, Sun X, Showers W J. 2009. Deep-sea benthic diversity linked to seasonality of pelagic productivity. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 56(5): 835-841. DOI:10.1016/j.dsr.2008.12.009
Cornelius N, Gooday A J. 2004. 'Live' (stained) deep-sea benthic foraminiferans in the western Weddell Sea: trends in abundance, diversity and taxonomic composition along a depth transect. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 51(14-16): 1 571-1 602. DOI:10.1016/j.dsr2.2004.06.024
Danovaro R, Carugati L, Corinaldesi C, Gambi C, Guilini K, Pusceddu A, Vanreusel A. 2013. Multiple spatial scale analyses provide new clues on patterns and drivers of deep-sea nematode diversity. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 92: 97-106. DOI:10.1016/j.dsr2.2013.03.035
Edgar R C. 2016. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. BioRxiv: 081257. DOI:10.1101/081257
Edgar R C. 2018. Updating the 97% identity threshold for 16S ribosomal RNA OTUs. Bioinformatics, 34(14): 2 371-2 375. DOI:10.1093/bioinformatics/bty113
Egge E S, Johannessen T V, Andersen T, Eikrem W, Bittner L, Larsen A, Sandaa R A, Edvardsen B. 2015. Seasonal diversity and dynamics of haptophytes in the Skagerrak, Norway, explored by high-throughput sequencing. Molecular ecology, 24(12): 3 026-3 042. DOI:10.1111/mec.13160
Elbrecht V, Leese F. 2015. Can DNA-based ecosystem assessments quantify species abundance? Testing primer bias and biomass-sequence relationships with an innovative metabarcoding protocol. PLoS One, 10(7): e0130324. DOI:10.1371/journal.pone.0130324
Epp D, Smoot N C. 1989. Distribution of seamounts in the North Atlantic. Nature, 337(6204): 254-257. DOI:10.1038/337254a0
Frontalini F, Cordier T, Balassi E, Du Chatelet E A, Cermakova K, Apothéloz-Perret-Gentil L, Martins M V A, Bucci C, Scantamburlo E, Treglia M, Bonamin V, Pawlowski J. 2020. Benthic foraminiferal metabarcoding and morphology-based assessment around three offshore gas platforms: congruence and complementarity. Environment International, 144: 106049. DOI:10.1016/j.envint.2020.106049
Fuhrman J A. 2009. Microbial community structure and its functional implications. Nature, 459(7244): 193-199. DOI:10.1038/nature08058
Goineau A, Gooday A J. 2019. Diversity and spatial patterns of foraminiferal assemblages in the eastern Clarion-Clipperton zone (abyssal eastern equatorial Pacific). Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 149: 103036. DOI:10.1016/j.dsr.2019.04.014
Gooday A J, Hori S, Todo Y, Okamoto T, Kitazato H, Sabbatini A. 2004. Soft-walled, monothalamous benthic foraminiferans in the Pacific, Indian and Atlantic Oceans: aspects of biodiversity and biogeography. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 51(1): 33-53. DOI:10.1016/j.dsr.2003.07.002
Gooday A J, Jorissen F J. 2012. Benthic foraminiferal biogeography: controls on global distribution patterns in deep-water settings. Annual Review of Marine Science, 4(1): 237-262. DOI:10.1146/annurev-marine-120709-142737
Gooday A J, Kitazato H, Hori S, Toyofuku T. 2001. Monothalamous soft-shelled foraminifera at an abyssal site in the north Pacific: a preliminary report. Journal of Oceanography, 57(3): 377-384. DOI:10.1023/A:1012447015771
Gooday A J. 2002. Organic-walled allogromiids: aspects of their occurrence, diversity and ecology in marine habitats. Journal of Foraminiferal Research, 32(4): 384-399. DOI:10.2113/0320384
Guillou L, Bachar D, Audic S, Bass D, Berney C, Bittner L, Boutte C, Burgaud G, De Vargas C, Decelle J, Del Campo J, Dolan J R, Dunthorn M, Edvardsen B, Holzmann M, Kooistra W H C F, Lara E, Le Bescot N, Logares R, Mahé F, Massana R, Montresor M, Morard R, Not F, Pawlowski J, Probert I, Sauvadet A L, Siano R, Stoeck T, Vaulot D, Zimmermann P, Christen R. 2012. The Protist Ribosomal Reference database (PR2): a catalog of unicellular eukaryote small sub-unit rRNA sequences with curated taxonomy. Nucleic Acids Research, 41(D1): D597-D604. DOI:10.1093/nar/gks1160
Heinz P, Ruepp D, Hemleben C. 2004. Benthic foraminifera assemblages at great meteor seamount. Marine Biology, 144(5): 985-998. DOI:10.1007/s00227-003-1257-7
Herndl G J, Reinthaler T. 2013. Microbial control of the dark end of the biological pump. Nature Geoscience, 6(9): 718-724. DOI:10.1038/ngeo1921
Jian Z M, Wang Y, Dang H W, Lea D W, Liu Z Y, Jin H Y, Yin Y Q. 2020. Half-precessional cycle of thermocline temperature in the western equatorial Pacific and its bihemispheric dynamics. Proceedings of the National Academy of Sciences of the United States of America, 117(13): 7 044-7 051. DOI:10.1073/pnas.1915510117
Jorissen F J, Fontanier C, Thomas E. 2007. Paleoceanographical proxies based on deep-sea benthic foraminiferal assemblage characteristics. Developments in Marine Geology, 1: 263-325. DOI:10.1016/S1572-5480(07)01012-3
Katz M E, Cramer B S, Franzese A, Hönisch B, Miller K G, Rosenthal Y, Wright J D. 2010. Traditional and emerging geochemical proxies in foraminifera. Journal of Foraminiferal Research, 40(2): 165-192. DOI:10.2113/gsjfr.40.2.165
Kozich J J, Westcott S L, Baxter N T, Highlander S K, Schloss P D. 2013. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Applied and Environmental Microbiology, 79(17): 5 112-5 120. DOI:10.1128/AEM.01043-13
Lecroq B, Lejzerowicz F, Bachar D, Christen R, Esling P, Baerlocher L, Østerås M, Farinelli L, Pawlowski J. 2011. Ultra-deep sequencing of foraminiferal microbarcodes unveils hidden richness of early monothalamous lineages in deep-sea sediments. Proceedings of the National Academy of Sciences of the United States of America, 108(32): 13 177-13 182. DOI:10.1073/pnas.1018426108
Lecroq B. 2014. Molecular assessment of benthic foraminiferal diversity. In: Kitazato H, Bernhard J M eds. Approaches to Study Living Foraminifera. Tokyo: Springer, 91-102.
Lee J J, Anderson O R. 1991. Biology of Foraminifera. Academic Press, London. p. 1-368.
Lei Y L, Li T G, Jian Z M, Nigam R. 2017. Taxonomy and distribution of benthic foraminifera in an intertidal zone of the Yellow Sea, PR China: correlations with sediment temperature and salinity. Marine Micropaleontology, 133: 1-20. DOI:10.1016/j.marmicro.2017.04.005
Lei Y L, Stumm K, Wickham S A, Berninger U G. 2014. Distributions and biomass of benthic ciliates, foraminifera and amoeboid protists in marine, brackish, and freshwater sediments. Journal of Eukaryotic Microbiology, 61(5): 493-508. DOI:10.1111/jeu.12129
Lejzerowicz F, Esling P, Pawlowski J. 2014. Patchiness of deep-sea benthic Foraminifera across the Southern Ocean: insights from high-throughput DNA sequencing. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 108: 17-26. DOI:10.1016/j.dsr2.2014.07.018
Lejzerowicz F, Pawlowski J, Fraissinet-Tachet L, Marmeisse R. 2010. Molecular evidence for widespread occurrence of foraminifera in soils. Environmental Microbiology, 12(9): 2518-2526. DOI:10.1111/j.1462-2920.2010.02225.x
Li L Q, Keller G. 1998. Maastrichtian climate, productivity and faunal turnovers in planktic foraminifera in south Atlantic DSDP sites 525A and 21. Marine Micropaleontology, 33(1-2): 55-86. DOI:10.1016/S0377-8398(97)00027-3
Li Q X, Lei Y L, Morard R, Li T G, Wang B D. 2020. Diversity hotspot and unique community structure of foraminifera in the world's deepest marine blue hole-Sansha Yongle Blue Hole. Scientific Reports, 10(1): 10257. DOI:10.1038/s41598-020-67221-0
Logares R, Audic S, Bass D, Bittner L, Boutte C, Christen R, Claverie J M, Decelle J, Dolan J R, Dunthorn M, Edvardsen B, Gobet A, Kooistra W H C F, Mahé F, Not F, Ogata H, Pawlowski J, Pernice M C, Romac S, Shalchian-Tabrizi K, Simon N, Stoeck T, Santini S, Siano R, Wincker P, Zingone A, Richards T A, De Vargas C, Massana R. 2014. Patterns of rare and abundant marine microbial eukaryotes. Current Biology, 24(8): 813-821. DOI:10.1016/j.cub.2014.02.050
Magoč T, Salzberg S L. 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics, 27(21): 2 957-2 963. DOI:10.1093/bioinformatics/btr507
Mojtahid M, Griveaud C, Fontanier C, Anschutz P, Jorissen F J. 2010. Live benthic foraminiferal faunas along a bathymetrical transect (140-4800 m) in the Bay of Biscay (NE Atlantic). Revue de Micropaléontologie, 53(3): 139-162. DOI:10.1016/j.revmic.2010.01.002
Morard R, Lejzerowicz F, Darling K F, Lecroq-Bennet B, Pedersen M W, Orlando L, Pawlowski J, Mulitza S, De Vargas C, Kucera M. 2017. Planktonic foraminifera-derived environmental DNA extracted from abyssal sediments preserves patterns of plankton macroecology. Biogeosciences, 14(11): 2. DOI:10.5194/bg-14-2741-2017
Naidu P D, Malmgren B A. 1996. A high-resolution record of late Quaternary upwelling along the Oman Margin, Arabian Sea based on planktonic foraminifera. Paleoceanography, 11(1): 129-140. DOI:10.1029/95PA03198
Ni J J, Yan Q Y, Yu Y H. 2013. How much metagenomic sequencing is enough to achieve a given goal?. Scientific Reports, 3(1): 1 968. DOI:10.1038/srep01968
Nienstedt J C, Arnold A J. 1988. The distribution of benthic foraminifera on seamounts near the East Pacific Rise. Journal of Foraminiferal Research, 18(3): 237-249. DOI:10.2113/gsjfr.18.3.237
Ohkushi K, Natori H. 2001. Living benthic foraminifera of the Hess Rise and Suiko Seamount, central North Pacific. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 48(5): 1 309-1 324. DOI:10.1016/S0967-0637(00)00081-9
Parker T, Tunnicliffe V. 1994. Dispersal strategies of the biota on an oceanic seamount: implications for ecology and biogeography. The Biological Bulletin, 187(3): 336-345. DOI:10.2307/1542290
Pawlowski J, Esling P, Lejzerowicz F, Cedhagen T, Wilding T A. 2014. Environmental monitoring through protist next-generation sequencing metabarcoding: assessing the impact of fish farming on benthic foraminifera communities. Molecular Ecology Resources, 14(6): 1 129-1 140. DOI:10.1111/1755-0998.12261
Pawlowski J, Fontaine D, da Silva A A, Guiard J. 2011. Novel lineages of Southern Ocean deep-sea foraminifera revealed by environmental DNA sequencing. Deep Sea Research Part Ⅱ: Topical Studies in Oceanography, 58(19-20): 1 996-2 003. DOI:10.1016/j.dsr2.2011.01.009
Pawlowski J, Holzmann M. 2007. Diversity and geographic distribution of benthic foraminifera: a molecular perspective. In: Foissner W, Hawksworth D L eds. Protist Diversity and Geographical Distribution. Dordrecht: Springer, p. 83-94.
Pawluczyk M, Weiss J, Links M G, Aranguren M E, Wilkinson M D, Egea-Cortines M. 2015. Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples. Analytical and Bioanalytical Chemistry, 407(7): 1 841-1 848. DOI:10.1007/s00216-014-8435-y
Pedrós-Alió C. 2006. Marine microbial diversity: can it be determined?. Trends in Microbiology, 14(6): 257-263. DOI:10.1016/j.tim.2006.04.007
Pietramellara G, Ascher J, Borgogni F, Ceccherini M T, Guerri G, Nannipieri P. 2009. Extracellular DNA in soil and sediment: fate and ecological relevance. Biology and Fertility of Soils, 45(3): 219-235. DOI:10.1007/s00374-008-0345-8
Piñol J, Mir G, Gomez-Polo P, Agustí N. 2015. Universal and blocking primer mismatches limit the use of high-throughput DNA sequencing for the quantitative metabarcoding of arthropods. Molecular Ecology Resources, 15(4): 819-830. DOI:10.1111/1755-0998.12355
Piredda R, Claverie J M, Decelle J, De Vargas C, Dunthorn M, Edvardsen B, Eikrem W, Forster D, Kooistra W H C F, Logares R, Massana R, Montresor M, Not F, Ogata H, Pawlowski J, Romac S, Sarno D, Stoeck T, Zingone A. 2018. Diatom diversity through HTS-metabarcoding in coastal European seas. Scientific Reports, 8(1): 18059. DOI:10.1038/s41598-018-36345-9
Polz M F, Cavanaugh C M. 1998. Bias in template-to-product ratios in multitemplate PCR. Applied and Environmental Microbiology, 64(10): 3 724-3 730. DOI:10.1128/AEM.64.10.3724-3730.1998
Saïdi E, Zaghbib-Turki D. 2016. Planktonic foraminiferal biostratigraphy and quantitative analysis during the Campanian-Maastrichtian transition at the Oued Necham section (Kalâat Senan, central Tunisia). Turkish Journal of Earth Sciences, 25: 538-572. DOI:10.3906/yer-1602-13
Shi J F, Lei Y L, Li Q X, Lyu M, Li T G. 2020. Molecular diversity and spatial distribution of benthic foraminifera of the seamounts and adjacent abyssal plains in the tropical Western Pacific Ocean. Marine Micropaleontology, 156: 101850. DOI:10.1016/j.marmicro.2020.101850
Stackebrandt E, Goebel B M. 1994. Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. International Journal of Systematic and Evolutionary Microbiology, 44(4): 846-849. DOI:10.1099/00207713-44-4-846
Stocks K I. 2004a. Seamount invertebrates: composition and vulnerability to fishing. In: Morato T, Pauly D eds. Fiseries Centre Research Reports: Seamounts: Biodiversity and Fisheries. Vancouver, Canada: Fisheries Centre, University of British Columbia, p. 17-24.
Stocks K I. 2004b. SeamountsOnline: an online information system for seamount biology. In: Vanden Berghe E ed. Proceedings of 'The Colour of Ocean Data'. World Wide Web Electronic Publication, 77-89.
Thomsen P F, Willerslev E. 2015. Environmental DNA - an emerging tool in conservation for monitoring past and present biodiversity. Biological Conservation, 183: 4-18. DOI:10.1016/j.biocon.2014.11.019
Torti A, Lever M A, Jørgensen B B. 2015. Origin, dynamics, and implications of extracellular DNA pools in marine sediments. Marine Genomics, 24: 185-196. DOI:10.1016/j.margen.2015.08.007
Ujiié Y, Kimoto K, Pawlowski J. 2008. Molecular evidence for an independent origin of modern triserial planktonic foraminifera from benthic ancestors. Marine Micropaleontology, 69(3-4): 334-340. DOI:10.1016/j.marmicro.2008.09.003
Vincent E, Berger W H. 1981. Planktonic foraminifera and their use in paleoceanography. In: Emiliani C ed. The Oceanic Lithosphere, vol 7, The Sea. New York: John Wiley, p. 1 025-1 119.
Westcott S L, Schloss P D. 2017. OptiClust, an improved method for assigning amplicon-based sequence data to operational taxonomic units. mSphere, 2(2): e00073-17. DOI:10.1128/mSphereDirect.00073-17
Yesson C, Clark M R, Taylor M L, Rogers A D. 2011. The global distribution of seamounts based on 30 arc seconds bathymetry data. Deep Sea Research Part Ⅰ: Oceanographic Research Papers, 58(4): 442-453. DOI:10.1016/j.dsr.2011.02.004
Zhang J J N. 1985. Living planktonic foraminifera from the eastern Arabian Sea. Deep Sea Research Part A. Oceanographic Research Papers, 32(7): 789-798. DOI:10.1016/0198-0149(85)90115-3
Zhao F, Filker S, Stoeck T, Xu K D. 2017. Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean. BMC Microbiology, 17(1): 192. DOI:10.1186/s12866-017-1103-6