Comparison in diversity of eukaryotic algae in surface sediments from diff erent functional sea areas of Qingdao coast, the Yellow Sea, China: a metabarcoding approach*

Zhaohui WANG , Mingdan LEI , Shuanghui JI , Changliang XIE , Jiazhuo CHEN ,Weiguo LI , Tao JIANG

1 College of Life Science and Technology, Jinan University, Guangzhou 510632, China

2 Engineering Research Center of Tropical and Subtropical Aquatic Ecological Engineering, Ministry of Education, Guangzhou 510632, China

3 School of Ocean, Yantai University, Yantai 264005, China

Abstract Surface sediment samples were collected in three diff erent functional sea areas in Qingdao coast, East China, including the inner Jiaozhou Bay, the Laoshan Coast, and the Amphioxus Reserve area.Diversity and community structure of eukaryotic algae especially those of phytoplankton resting stages were assessed by metabarcoding V4 region of the 18S rDNA. Biogenic elements including total organic carbon(TOC), organic matter (OM), total nitrogen (TN), total phosphorus (TP), and biogenic silicon (BSi) were analyzed. A total of 1 496 eukaryotic operational taxonomic units (OTUs) were measured, including 207 algal OTUs, which contributed to 13.84% of the total OTUs. Ninety-eight species in 8 phyla, 24 classes of eukaryotic algae were detected. Among them, 47 species have been reported to form resting stages, and 12 species are f irstly recorded in Chinese coastal waters. Dinof lagellates dominated in both DNA reads and OTU richness, which contributed to 73.02% and 61.35% of the eukaryotic algal sequences and OTU richness, respectively. DNA reads, OTU richness and alpha diversity indexes of eukaryotic algae were higher in the Laoshan Coast, and lower in Jiaozhou Bay. Eukaryotic algal community diff ered in the three sea areas, which was dominated by chrysophytes in Jiaozhou Bay, by dinof lagellates in the Laoshan Coast,and co-dominated by dinof lagellates and chrysophytes in the Amphioxus Reserve area. Clustering analysis showed that the Laoshan Coast and the Amphioxus Reserve area are clustered together, while Jiaozhou Bay is clustered separately. Thirty-six harmful algal bloom (HAB) species were detected, and 10 species have been reported to form blooms in Jiaozhou Bay and the Qingdao coast before. Some of these species occurred widely and dominantly in this study, suggesting high potential risk of HABs in the Qingdao coastal area.

Keyword: benthic microalgae; 18S rDNA; metabarcoding; resting stages; Jiaozhou Bay; Yellow Sea

1 INTRODUCTION

Marine sediments provide luxuriant nutrients for benthic organisms, and are considered as “seed banks” of biodiversity in the seawater, out of which species can reemerge into the water column in appropriate opportunity (Massana et al., 2015). As the main producers of benthic eukaryotic communities,eukaryotic algae are major components in microbial food webs, playing important roles in energy f low and biogeochemical circle in the benthic ecosystems (Bik et al., 2012). Benthic microalgae have been regarded as an important ecological indicator for assessing anthropogenic eff ects on the marine environments(Blanco et al., 2008; Fierro et al., 2019).

Many phytoplankton species have the ability to experience the dormant phase as a part of their life circle under the stresses of adverse environmental conditions in the forms of resting spores and resting cysts (hereafter referred to as resting stages) (Head,1996). The resting stages thicken the cell wall and sink to the sea f loor, and then become a part of benthic community (Piredda et al., 2017). The resting stages have thick resistant wall, and can withstand low temperature, high pressure, hypoxic and dark benthic environment, and can survive in sediments for hundreds of years (reviewed by Ellegaard and Ribeiro,2018). When the environmental conditions are suitable, the resting stages germinate and release vegetative cells into the water column (McQuoid,2002). Phytoplankton resting stages are regarded as the “seed banks” of algal blooms, and play important role in initiating and dissipation of harmful algal blooms (HABs) (Smayda, 2002; Genovesi et al.,2009).

Identif ication of benthic microalgae has traditionally been based on microscopic observations and cell counting after separating cells from sediments. However, this approach has signif icant limitations, such as the recovery rate during separation and concentration, biases in the recovered groups,and diffi culties associated with morphological identif ication (Piredda et al., 2017). These limitations have been overcome with the implementation of modern molecular based metabarcoding techniques,particular for the utilization of high-throughput sequencing (HTS) techniques (Penna et al., 2017).Metabarcoding allows identif ication of microorganisms without taxonomic expertise, and only need to match the HTS-derived DNA fragments to the reference sequence libraries (Keeley et al.,2018). Standardized protocols can be developed and the results are auditable (Ji et al., 2013). HTS based metabarcoding has become a powerful and rapid technique to identify and quantify benthic organisms for biomonitoring purposes (Keeley et al., 2018).

The coastal zone is the most active area of human activity, and provides abundant resources and have high ecological value and environmental functions.However, many coastal environments have been suff ering from problems of environmental degradation,decline in biodiversity, and species invasion as the rapid economic development, population expansion,and over-exploitation of marine resources (Yu et al.,2019). Qingdao is located in the eastern China,southeast of Shandong Peninsula and east of the Yellow Sea. It is the economic center of Shandong Province as the important international port, the modern marine industry zone, the international shipping hub in Northeast Asia, and the maritime sports base. Qingdao City covered 758.16 km2area with the population of 9.50 million in 2019. Due to the rapid increase of economic development and population size as well as the increase of aquaculture,the pollution level in Qingdao coastal waters has been increasing recently (Yuan et al., 2017). The degradation of water quality results in the frequent occurrence of algal blooms. A total of 36 algal blooms occurred in the coastal areas of Qingdao between 1990 and 2017,in which more than 80% occurred in Jiaozhou Bay(Zhou et al., 2020). Moreover, sediments in Jiaozhou Bay have been polluted to a certain extent (Gu et al.,2021; Wang et al., 2021).

Intensive surveys have been carried out to investigate phytoplankton community in Jiaozhou Bay during the past f ifty years, and results have shown rich phytoplankton diversity and great variations of the dominant phytoplankton taxa with time (reviewed by Liu and Chen, 2021). Meanwhile, metabarcoding studies revealed much more phytoplankton species in Jiaozhou Bay than the routine microscopic observations (Liu et al., 2020b; Chen et al., 2021).However, there are few studies on phytoplankton resting stages in sediments, with only one report on dinof lagellate cysts in Jiaozhou Bay based on microscopic analysis (Li et al., 2017), which cannot fully ref lect the diversity and distribution of phytoplankton resting stages in sediments. In order to understand the distribution of phytoplankton resting stages and the potential of HABs in the Qingdao coast, surface sediments were collected from three diff erent function sea areas of the Qingdao coast in this study, including a eutrophic enclosed bay (the inner Jiaozhou Bay), a conservation area (the Qingdao Amphioxus Reserve area), and a scenic and aquacultural area (the Laoshan Coast). HTS-based metabarcoding was used to investigate diversity and community structure of eukaryotic algae. Biogenic elements including total organic carbon (TOC),organic matter (OM), total nitrogen (TN), total phosphorus (TP), and biogenic silicon (BSi) were analyzed, and the relationships between algal community and biogenic elements were discussed.The purpose of this study was to compare benthic microalgae composition among the three sea areas with diff erent degrees of human activities, and to discuss the distribution of resting stages of HAB species in sediments from the Qingdao coast. This study represents a f irst attempt to study eukaryotic algae in sediments from the Qingdao coast using the metabarcoding approach.

Fig.1 Sampling stations in Qingdao coastal sea area

2 MATERIAL AND METHOD

2.1 Study areas and sediment collection

Qingdao is located in the southeast of the Shandong Peninsula, which borders the Yellow Sea in the east and the south. Surface sediments were collected in three sea areas along the Qingdao coast, i.e., the inner Jiaozhou Bay, the Laoshan Coast, and the Qingdao Amphioxus Reserve area (Fig.1). Jiaozhou Bay, in the northwest of Qingdao city center, is a large semiclosed bay connected to the Yellow Sea by a narrow bay mouth. Jiaozhou Bay has multiple functions, such as port, transportation, and aquaculture. The Laoshan Coast locates in the eastern part of Qingdao City and surround the famous tourist attraction the Laoshan Scenic Area. There are many aquaculture farms in the Laoshan Coast. The Amphioxus Reserve area is located in the southern part of Qingdao City, just outside the mouth of the Jiaozhou Bay and close to the Qingdao port. The Amphioxus Reserve area was established in August 2004 to protect wild amphioxus population and its habitat environment.cold (4 °C) dark condition for six months before DNA extraction to inactivate most of the microalgal vegetative cells, and to make sure that most of the survives were resting stages (Piredda et al., 2017).

2.2 DNA extraction, PCR amplification, and sequencing

Genomic DNA of ca. 1.0 g of sediment sample was extracted using DNeasy PowerSoil Kit (Qiagen,Germany) according to the manufacturer’s instructions. Polymerase chain reaction (PCR)amplif ication of ca. 450-bp fragment of the hypervariable V4 region of the 18S rRNA gene(rDNA) was carried out with the universal eukaryotic primers 3NDf 5ʹ-GGCAAGTCTGGTGCCAG-3ʹ(Cavalier-Smith et al., 2009) and V4_euk_R2 5ʹ-ACGGTATCTATCTCTTCG-3ʹ (Bråte et al.,2010). PCR products were then sequenced on an Illumina MisSeq PE300 platform (Illumina, San Diego, CA). Sequencing was performed by Majorbio Bio-pharm Technology Co., Ltd. (Shanghai, China).

Surface sediments were collected in triplicate from seventeen stations in the three sea areas using a Peterson grab between November and December 2017 (Fig.1). The top 2 cm of sediments were sampled with a polyethylene spatula, and the three replicate samples were mixed well together as one sample and sealed in a polyethylene bag, and then stored in the

2.3 Sequence data cleaning, f iltering, and OTU analysis

Raw data from the Illumina libraries were assigned to individual samples by their barcodes. Paired-end reads were truncated to remove unique barcode and primer sequences, then assembled using FLASH(Version 1.2.7). The obtained f iltered sequences with error were removed via using QIIME (QIIME, version 1.8, http://qiime.org) software. The Operational Taxonomic Unit (OTU) was the sequence similarity threshold, and sequences with similarities above the threshold were merged into the same OTU at 97%sequence similarity by using UPARSE (version 7.1,http://drive5.com/uparse/). Singletons and doubletons were removed from the dataset, and the chimeras were removed using USEARCH (version 8.1.1861,Http://www.drive5.com/Usearch/). The taxonomic assignment of each OTU was classif ied against the PR2 database (version 4.12.0, https://pr2-database.org/, searched on 27 March 2021) using the Ribosomal Database Project (RDP) classif ier (Caporaso et al.,2010) with a bootstrap conf idence of 80%. Those OTUs, which were not assigned to the species level in PR2 database, were further annotated against the NCBI database (https://www.ncbi.nlm.nih.gov/,searched on 17 April 2021) with the highest similarity( ≥ 97%). The raw data were deposited into the National Center for Biotechnology Information (NCBI)Sequence Read Archive (SRA) (https://www.ncbi.nlm.nih.gov/sra/) with the accession number PRJNA722324. The classif ication of algal taxa is based on taxonomic information in AlgaeBase, a global algal database (Guiry and Guiry, 2021, https://www.algaebase.org/). A taxon was annotated as an HAB species if it is categorized as harmful microalgae in the online version of the Intergovernmental Oceanographic Commission (IOC) Taxonomic Reference List of Harmful Micro Algae (Lundholm et al., 2009 onwards) or in “Manual on Harmful Marine Microalgae” (Hallegraeff et al., 2003) or had been reported as an HAB species in previous studies.

2.4 Analyses of biogenic elements

Sediments for biogenic elements analysis were dried in an oven at 40 °C until a constant weight was reached. The dried sediments were ground gently with an agate mortar and pestle, sieved through a 100-μm mesh for homogenization and stored in sealed plastic bags. The biogenic elements were measured in triplicate. TOC and TN was measured by a Perkin-Elmer 2400 Series II CHNS/O Analyzer (Perkin Elmer Inc., USA). TP was measured by potassium persulfate digestion method (Thien and Myers, 1992). OM was determined by ignition loss in a muffl e furnace (SGXL1200, Honglang, Shanghai, China). BSi was measured by the molybdate blue spectrophotometric method after removing the carbonates and organics by 1-mol/L HCl and 10% H2O2and digested using 0.5-mol/L Na2CO3solution (Mortlock and Froelich,1989). The quality assurance/quality control (QA/QC)was assessed by the analyses of blank reagents and f ive replicates of the certif ied reference material(Off shore Marine Sediment, GBW 07314). The analytical precision was controlled to within 5% for biogenic elements.

2.5 Statistical analyses

Rarefaction curves were generated to assess the degree of sample saturation using the R 4.0.0 package vegan. The alpha diversity of eukaryotic organisms and eukaryotic algae, including OTU richness,Margalef, Simpson diversity, Shannon diversity,Pielou index, Chao1 and ACE, were analyzed after rarefying each sample to the number of DNA reads in the sample with the fewest reads (35 984 reads in this study). The Bray-Curtis (BC) index was used as a measure of similarity between the samples (betadiversity). A distance matrix was computed with the BC index, and a hierarchical cluster tree was constructed using the un-weighted pair group method with arithmetic mean (UPGMA) with bootstrap support. Venn diagrams were calculated and plotted based on diff erent samples and sea areas. All these analyses were conducted using the R 4.0.0 package vegan. A phylogenetic tree was constructed to illustrate the relationships between the algal communities in surface sediments from the Qingdao coast using the Neighbor-Joining algorithm method with the maximum likelihood model (MEGA v. 10.0 software) with bootstrap values of 1 000 replicate runs. A bubble chart of the top 25 abundance OTUs was drawn with the R package ggplot2, reshape2, and formats according to their abundances, and the DNA reads were logarithmically transformed. The redundancy analysis (RDA) was performed to revel the impacts of biogenic elements in explaining the distribution patterns of algal communities in seventeen surface sediments using Canoco4.5. All data were logarithmically transformed to obtain the equal weight of the elements in RDA analysis.

3 RESULT

3.1 Sequence data

Fig.2 Rarefaction curves of DNA reads based on OTUs at 97% sequence similarity

A total of 947 886 eukaryotic DNA sequences and 1 497 OTUs were obtained in 17 surface sediment samples from the Qingdao coast. Rarefaction curves suggested adequate coverage for all datasets (Fig.2).DNA reads in each sample ranged from 35 984 at station L1 to 72 247 at station L3, and OTU richness were between 175 and 684 OTUs, with the lowest at station J3 and the highest at station L2, respectively(Figs.2–3). There were 171 722 algal sequences and 207 algal OTUs. The number of algal sequences at each station varied greatly (Fig.3a), ranging from 424 to 28 082 reads with the highest at station L6 and the lowest at station J3. The proportions of algal reads to the total eukaryotic reads ranged from 0.59% at station J3 to 41.25% at station L6. The algal OTU richness of each sample ranged from 20 OTUs at station J3 to 103 OTUs at station L6 with the proportions of 8.93%–20.80% (Fig.3b).

3.2 Eukaryotic algal diversity and community structure

The eukaryotic algae detected in this study included 75 genera in 24 classes of 8 phyla, i.e., Bacillariophyceae and Mediophyceae in Bacillariophyta, Chlorophyceae,Chlorodendrophyceae, Chloropicophyceae,Mamiellophyceae, Palmophyllophyceae,Pedinophyceae, Pyramimonadophyceae,Trebouxiophyceae, and Ulvophyceae in Chlorophyta,Cryptophyceae in Cryptophyta, Dinophyceae, and Syndiniophyceae in Dinophyta, Prymnesiophyceae in Haptophyta, Chrysophyceae, Eustigmatophyceae,Pelagophyceae, Phaeophyceae, Raphidophyceae, and Xanthophyceae in Ochrophyta, Prasinodermatophyceae in Prasinodermatophyta, and Florideophyceae and Rhodellophyceae in Rhodophyta. There were 124 OTUs (59.9%) identif ied at the species level after BLAST against the PR2 and the NCBI databases, and 161 OTUs (77.8%) at the genus level, and other 46 OTUs (22.2%) only into class level. Sometimes more than one or two OTUs were annotated as the same taxa, and altogether 98 species were detected in this study, including 12 Bacillariophyta (13 OTUs), 23 Chlorophyta (27 OTUs), 2 Cryptophyta (2 OTUs), 45 Dinophyta (66 OTUs), 1 Haptophyta (1 OTU), 13 Ochrophyta (13 OTUs), and 2 Rhodophyta (3 OTUs)species (Supplementary Table S1). There were 47, 89,and 61 species recorded in samples from the inner Jiaozhou Bay, the Laoshan Coast, and the Amphioxus Reserve area, respectively. Most of the species or genera identif ied in surface sediments in this study were marine species, some were brackish or marine/freshwater species, and terrestrial and benthic species were also detected. About 47 taxa have been reported to form resting stages, and 12 species are f irstly recorded in the Chinese coastal waters (Supplementary Table S1).

Fig.3 DNA reads (a) and OTU richness (b) of eukaryotic organisms and eukaryotic algae in surface sediments from 17 sites of the Qingdao coastal

Based on DNA reads at the class level (Fig.4a),Dinophyceae dominated in the eukaryotic algal assemblages in surface sediments from the Qingdao coast, with the total of 124 516 reads accounting for 72.51% of total algal sequences. The second dominant class was Chrysophyceae, with 41 137 reads and proportion of 23.96%. Diatoms, including two classes, Mediophyceae and Bacillariophyceae, in this study, ranked the third with 2 197 reads and proportion of 1.27%. There were 1 178 reads of raphidophytes with the proportion of 0.69%. Syndiniophytes contributed 0.51% to the total algal reads. Sequences of other classes were less than 500 reads and contributed low than 0.5% of the algal reads.

Fig.4 Overall relative abundances of DNA reads (a) and OTU richness (b) of eukaryotic algae at the class level

Dinof lagellates (Dinophyta, including Dinophyceae and Syndiniophyceae) were the most diverse algal group in this study, with 106 OTUs and 21 OTUs,respectively, together which represented 61.35% of the total algae OTUs (Fig.4b). However, none of the syndiniophytes were identif ied at the species level.Thirty-three OTUs (15.94%) were recorded in nine classes of Chlorophyta. Among them, 10 OTUs(4.83%) were recorded in the marine green algal class Mamiellophyceae (Fig.4b). Diatoms were the third diverse group, with 12 OTUs (5.80%) and 6 OTUs(2.90%) for Mediophyceae and Bacillariophyceae,respectively.

3.3 Comparison of eukaryotic algal community structure in surface sediments from the three sea areas of the Qingdao coast

To further understand the relationships among eukaryotic algae, a phylogenetic tree was constructed based on all algal OTUs detected in this study(Fig.5). All OTUs in Dinophyta cluster into a big group, while OTUs in the other seven phyla form up another big group. The second group is separated into two large clades: one including Chlorophyta and the other three phyla (Rhodophyta, Cryptophyta,and Haptophyta). However, Ochrophyta and Bacillariophyta cluster to form another clade. OTUs in Chlorophyta join together as a large sub-clade,and those in Rhodophyta, Cryptophyta, and Haptophyta cluster. Generally, OTUs of the same phylum clustered. However, an unclassif ied diatom OTU (OTU241) and three unclassif ied Ochrophyta OTUs cluster into a small branch in the diatom clade. In addition, an OTU of Prasinodermatophyta,a new algal phylum which was recently separated from Chlorophyta (Li et al., 2020), is within the large branch of Chlorophyta.

The composition of algal community was quite diff erent in the three sea areas of the Qingdao coast(Fig.6). The eukaryotic algal community in the inner Jiaozhou Bay (JZ) was dominated by chrysophytes(Fig.6a), which contributed to 77.40% of the total algal sequences; followed by dinof lagellates(including Dinophyceae and Syndiniophyceae)accounting for 17.79% of the algal sequences.However, dinof lagellates predominated in the Laoshan Coast (LS) with a relative percentage of 95.25%. Dinof lagellates and chrysophytes codominated in the Amphioxus Reserve area (AR) with the percentages of 62.98% and 32.28%, respectively.The compositions of OTU richness were comparable in the three sea areas (Fig.6b), in which higher diversities were observed in dinof lagellates(Dinophyceae and Syndiniophyceae), diatoms(Mediophyceae and Bacillariophyceae), and Mamiellophyceae of Chlorophyta. Dinof lagellates accounted for 47.76%, 63.64%, and 66.12% of the algal OTUs in JZ, LS, and AR, respectively.

The Venn diagrams highlighted the diff erences of algal communities among the seventeen samples in the three sea areas (Fig.7). Only 44 OTUs (21.3%)shared among the three sea areas, and 2 OTUs shared among all samples. Particularly, the Laoshan Coast(LS) had the highest number of OTUs, and 66 unique OTUs (31.9%) were detected in this sea area (Fig.7a).The numbers of unique OTUs were comparable in the inner Jiaozhou Bay (JZ) and the Amphioxus Reserve area (AR), which accounted for 4.3% of the total OTUs (Fig.7a). The reason for the high OTU richness in the Laoshan Coast was partly due to more samples collected (9 samples), while only 4 samples in the other two sea areas, respectively. In addition, the diversity of eukaryotic algae per sample was higher in the Laoshan Coast as well, and the top two highest diverse samples (L6 and L1) were all located in the Laoshan Coast (Fig.7b).

Fig.5 A phylogenetic tree on OTU lever based on eukaryotic algal sequences in surface sediments from the Qingdao coast using the Neighbor-Joining algorithm method with the maximum likelihood model (MEGA v.10.0)

Fig.6 Relative abundance of DNA reads (a) and OTU richness (b) of eukaryotic algae in surface sediments from the three sea areas of the Qingdao coast

Fig.7 Venn diagrams highlighting the degree of overlap of algal OTUs among the three sea areas (a) and among the seventeen samples (b)

The OTU richness, Chao1, ACE, Pielou, Shannon and Margalef indices were signif icantly lower in the inner Jiaozhou Bay than those in samples from the other two sea areas (P<0.05 for all), however the Simpson index was higher in the inner Jiaozhou Bay(Supplementary Table S2). The results suggested lower diversity and evenness, and higher dominance in the benthic eukaryotic algal community in the inner Jiaozhou Bay. Stations J3 and J4, located in the eastern part of the inner Jiaozhou Bay and close to Qingdao City, had the lowest diversity, and DNA reads and proportions of algal sequences were low in the two stations as well. The small phagotrophic chrysophyteParaphysomonaslucasiwas predominated at station J1, which contributed to 96.59% of the algal sequences, leading to very low Shannon and Pielou indices and high Simpson index at the station (Supplementary Table S2). There was no signif icant diff erence in diversity indices between samples from the Laoshan Coast and the Amphioxus Reserve area (P>0.05). The overall diversity indices in the Laoshan Coast were higher than those in the Amphioxus Reserve area due to more samples analyzed in the Laoshan Coast.

Fig.8 Cluster analysis of seventeen sediment samples from the three sea areas of Qingdao coast based on Bray-Curtis distance

Cluster analysis showed that thirteen samples from the Laoshan Coast (L1–L9) and the Amphioxus Reserve area (A1–A4) were clustered into one large group (Fig.8). Two samples from the inner Jiaozhou Bay (J1 and J2) were clustered into one group, while J3 and J4 were ungrouped. The results indicated the benthic algal communities were quite diff erent between the inner Jiaozhou Bay and the other two sea areas, and also among diff erent stations in the inner Jiaozhou Bay. Jiaozhou Bay is semi-enclosed sea area with only a narrow mouth open to the Yellow Sea,which makes it a specif ic ecosystem and microalgal community. While the Laoshan Coast and the Amphioxus Reserve area are connected to the Yellow Sea, which share similar microalgal compositions.

3.4 Distribution of dominant species

Figure 9 illustrates the DNA reads of the top 25 abundant OTUs. Most of the top abundant OTUs belong to dinof lagellates, except for two OTUs in Chrysophyceae, one OTU in diatoms, and two OTUs in Raphidophyceae. Most of the identif ied species with high abundance were those conf irmed to form resting stages (Supplementary Table S1), while those unidentif ied at the genus level were cyst/spore forming genera such asGonyaulaxandChaetoceros.Only two OTUs were present at all stations(Biecheleriahalophila, OTU443, OTU1061), four OTUs occurred at 16 stations, i.e.,Scrippsiellaacuminata(OTU1441),Gyrodinium spirale(OTU994),Gymnodiniummicroreticulatum(OTU1001), andHeterosigmaakashiwo(OTU1455).In addition, most of the other dominant OTUs appeared at more than half of the stations. Two dinof lagellate OTUs (OTU737, OTU874) occurred abundantly in a few samples however only classif ied into the class lever.

The small thin-walled and cyst-forming dinof lagellate,Biecheleriahalophila, was the f irst dominant species in the surface sediments from the Qingdao coastal areas. Two OTUs ofB.halophila(OTU443 and OTU1061) were recorded in this study based on 97% similarity, which were the f irst and fourth dominant OTUs, respectively, and both OTUs were presented abundantly at all stations (Fig.9).MeanwhileBiecheleriawas the most diverse genus in this study, and altogether eighteen OTUs were detected intoBiecheleria, ten of which were detected into four species (B.brevisulcata,B.cincta,B.halophila, andB.tirezensis).B.tirezensiswas f irstly recorded in the Chinese coastal waters(Supplementary Table S1).

The phagotrophic chrysophyteParaphysomonaslucasi, a newly recorded species in Chinese coasts,was the second dominant species in this study. It was detected in 12 stations in this study and was the predominant species in the inner Jiaozhou Bay and two stations in the Amphioxus Reserve area, especially at station J1, in which accounted for 96.93% of the total algal reads. The massive occurrence ofParaphysomonaslucasiresulted in chrysophytedominant algal community in the inner Jiaozhou Bay,and chrysophyte and dinophyte co-dominant community in the Amphioxus Reserve area (Fig.6a).

Scrippsiellaacuminata(formerlyS.trochoidea)was the third richest and most widely distributed species in this study.Scrippsiellais a dinof lagellate genus, which forms calcareous cysts (Matsuoka and Fukuyo, 2000). Four OTUs (OTU1441, OTU1122,OTU232, and OTU1108) were classif ied into genusScrippsiella, and three of them were identif ied at the species level (S. acuminata,S. erinaceus, andS.masanensis, Supplementary Table S1).

Gymnodiniumwas one of the most common dominant genera. Five OTUs were classif ied into genusGymnodiniumin this study, and four were identif ied at the species level (Gy.catenatum,Gy.corollarium,Gy.impudicum, andGy.microreticulatum). All these species form microreticulate cysts, which are diffi cult to be identif ied under microscopic observations (Gu et al., 2013a). Three species (Gy.corollarium,Gy.microreticulatum, andGy.impudicum) were recorded as the top 25 abundant OTUs (Fig.9).

Gonyaulaxwas one of the most diverse and abundant genera in this study. Most species ofGonyaulaxform ovoid cysts with numerous membranaceous processes, which belong to palynological genusSpiniferites. Thirteen OTUs were attributed toGonyaulax, and six OTUs were identif ied into four species (Go.cochlea,Go.fragilis,Go.spinifera, andGo.whaseongensis). Four OTUs includingGo.spinifera(OTU420, OTU414),Gonyaulaxsp. (OTU1104), andGo.whaseongensis(OTU191) were recorded as the top 25 abundance.

Azadiniumoccurred dominantly in this study as well, and three of the four OTUs ofAzadiniumwere detected at the species level, i.e.,Az.dalianense,Az.poporum, andAz.trinitatum.Az.trinitatum(OTU1437) andAz.dalianense(OTU157) were the top 8thand top 17thabundant OTUs, respectively. In addition,Az.trinitatumdistributed in 15 stations, andAz.dalianensewas recorded in all stations except for four stations in the inner Jiaozhou Bay. While the azaspiracid shellf ish poisoning (AZP) producerAz.poporum(Krock et al., 2012) occurred in only three samples with low reads.

Chaetoceroswas the most abundant and diverse diatom genus in this study, and eight OTUs were detected including six OTUs identif ied into f ive species (Chaetocerosdebilis,Chae.decipiens,Chae.muelleri,Chae.radicans, andChae.wighamii). Only one diatom OTU (Chaetocerossp., OTU1154) was within the top 25 abundances, and other diatoms occurred in low reads.

Fig.9 The bubble chart showing DNA reads of the top 25 abundant OTUs in each station

Fig.10 Redundancy analysis (RDA) of biogenic elements and the eukaryotic algal community in surface sediments from the Qingdao coasts, showing algal distribution on class level and the directions of biogenic elements to the f rist two RDA axes, axis 1 (61.9%) and axis 2 (20.7%)

3.5 Occurrence of harmful algal bloom species

Notably, 36 potentially toxic/harmful and/or bloom species were detected in surface sediments from the Qingdao coast (Supplementary Table S1), including the paralytic shellf ish poisoning (PSP) producersAlexandriumminutum,Al.andersonii, andGymnodiniumcatenatum, the yessotoxin (YTX)producersGonyaulaxspiniferaandProtoceratiumreticulatum, the AZA producerAzadiniumpoporum,the karlotoxin producerKarlodiniumvenef icum, the ichthyotoxic speciesAureococcusanophageff erens,Chattonellamarina,Fibrocapsajaponica,Heterosigmaakashiwo,Polykrikoshartmannii, andPo.geminatus. Sixteen bloom species occurred in all of the three sea areas, including three bloom causative diatoms, ten HAB dinof lagellates, the brown tide speciesAu.anophageff erens, and two ichthyotoxic raphidophytesChat.marinaandH.akashiwo. All of the 36 bloom species were recorded in the Laoshan Coast, and 20 and 27 bloom species were detected in samples from the inner Jiaozhou Bay and the Amphioxus Reserve area, respectively. Some of these species occurred widely and dominantly in this study,such asB.halophila,S.acuminata,Gy.corollarium,Go.spinifera,Levanderinaf issa,Gy.microreticulatum,Po.geminatus,Gy.impudicum,H.akashiwo,K.venef icum, andChat.marina(Fig.9).

3.6 RDA analyses of benthic algal community and biogenic elements

RDA analysis was conducted based on biogenic elements and benthic eukaryotic algal community(Fig.10). Axis 1 and Axis 2 explained 61.9% and 20.7% of the environmental and biological variables,respectively. Most of the biogenic elements and eukaryotic algal classes were located in the negative direction of Axis 1, and BSi had a higher negative contribution to Axis 1. Dinof lagellates (Dinophyceae and Syndiniophyceae) and diatoms (Bacillariophyceae and Mediophyceae) showed narrow intersection angles with BSi and TP, indicating signif icant inf luences of the two biogenic elements on the distribution of dinof lagellates and diatoms (Fig.10a).TOC, OM, and TN had a negative contribution to Axis 2. Two classes of green algae, Pyramimonadophyceae and Mamiellophyceae, had a small angle with these elements, which indicated that organic matters had a certain inf luence on the distribution of green algae.Other algal classes scattered in the coordinates far away from the biogenic elements, indicating few eff ects of biogenic elements on their distribution.

The ordination sampling station showed that samples in the three sea areas were separately grouped(Fig.10b). Samples in the inner Jiaozhou Bay were clustered along the positive axis of Axis 1, which indicated the low BSi content in this area. Samples in the Amphioxus Reserve area distributed along the positive axis of Axis 2, indicating the low organic matters (TOC, OM, and TP) in this area. Meanwhile,the arrows of the biogenic elements distributed within the group of samples of the Laoshan Coast, indicating the high contents of biogenic elements in this area.Station J3 was far away from other stations, which might be accounted by the inf luences of intensive human activities as it is very close to land. Furthermore,only 424 algal reads were recorded in J3, which was too low to form a strong correlation with biogenic elements.

4 DISCUSSION

4.1 Composition of eukaryotic algal community in surface sediments from the Qingdao coast

Dinof lagellates dominated in the benthic algal community structure in this study, followed by chrysophytes, diatoms, and chlorophytes.Dinof lagellates also had the highest OTU richness.The benthic algal community structure is comparable to those reported in sediments from the other coastal sea areas, which were dominated by dinof lagellates in both DNA reads and OTU richness (Dzhembekova et al., 2018; Liu et al., 2020a). Dinof lagellates have much higher DNA copies (up to thousands) than the other protists (dozens to hundreds) (Lin, 2011), which may lead to their over-representation in metabarcoding studies. In addition, about 20% dinof lagellates have the ability to form resting cysts, and dinof lagellate cysts constitute a major component of sedimentary assemblages of microalgal resting stages (reviewed by Ellegaard and Ribeiro, 2018). Actually, most dominant dinof lagellates detected at the species level in this study were cyst-forming taxa (Fig.9).

Chrysophyceae has been detected abundantly in the eukaryotic sequences in both sediment and water samples (de Vargas et al., 2015; Liu et al., 2020a).Cyst formation is a characteristic feature of chrysophytes, which form stomatocysts with siliceous wall (Findenig et al., 2010). Therefore, the chrysophyte community has been recorded relatively completely in metagenomic monitoring of sediments.Paraphysomonaslucasipredominated in the inner Jiaozhou Bay and also widely and abundantly distributed in the sediments from the Amphioxus Reserve area in this study. It is a new recorded species in the Chinese coastal sea area as far as our knowledge.Wide and massive occurrences of this species in this study suggested its ubiquitous distribution and cystformation characteristics. Meanwhile, chrysophytes are also an important group in phytoplankton community (Liu et al., 2020b; Pan et al., 2020).However, they are generally diffi cult to be distinguished and identif ied under the microscopic observations due to their fragility or their small size(LeRoi and Hallegraeff , 2006; Findenig et al., 2010).More and more chrysophyte species have been recorded in phytoplankton communities as the application of HTS-based metabarcoding (Scoble and Cavalier-Smith, 2014), for exampleParaphysomonasforaminifera,Pa.imperforate, andPedospumellaencystansin Jiaozhou Bay, the Yellow Sea (Liu et al.,2020b), andParaphysomonasvestitain Liaodong Bay, the Bohai Sea (Song et al., 2019). The results of these studies and our study suggested that the abundance and diversity of the nano-sized chrysophytes might be much greater than expected in the Yellow Sea and the Bohai Sea of China.

The OTU richness and abundance of benthic diatoms were greatly reduced in this study compared with the phytoplankton metabarcoding studies in Jiaozhou Bay and the other sea areas (Xu et al., 2017;Song et al., 2019; Liu et al., 2020b). We focused on phytoplankton resting stages in this study. However,only a few central diatom species have been reported to form resting spores (McQuoid and Hobson 1996).Therefore, most diatom cells existing in water column could not be detected in this study. Diatoms recorded in this study were those reported to form resting spores, such as species in generaChaetoceros,Skeletonema, andThalassiosiraetc. (Montresor et al.,2013; Piredda et al., 2017). Meanwhile, the abundance and species diversity of diatoms in the phytoplankton community in Jiaozhou Bay and its adjacent waters showed a decreasing trend due to the decrease of silicon f lux (Liu et al., 2010). Some freshwater and terrestrial species were detected in this study, which might be carried into the sea by surface runoff .

4.2 Comparison of benthic algal communities in the three sea areas

Although the three sea areas are not far apart, the benthic algal community structures were quite diff erent. Chrysophytes and dinof lagellates dominated in the inner Jiaozhou Bay and the Laoshan sea area,respectively. While, dinof lagellates and chrysophytes co-dominated in the Amphioxus Reserve area. Only 44 OTUs (21.3%) shared among the three sea areas,and 2 OTUs shared among all samples in this study.The cluster analysis showed that samples from Jiaozhou Bay were separately clustered. Jiaozhou Bay is an enclosed inner bay with little exchange with outside water bodies. The four stations of Jiaozhou Bay in this study are located in the inner Jiaozhou Bay. The enclosed environment results in a unique phytoplankton and benthic algal community structure.The predominant chrysophyteParaphysomonaslucasirecorded in this study is a ubiquitous phagotroph, which primarily feeds on bacteria (Scoble and Cavalier-Smith, 2014). Therefore, nutrient enrichment and high organic matters in Jiaozhou Bay(Li et al., 2020) provide suitable conditions for the phagotrophic chrysophytes. Actually, chrysophytes occurred abundantly in the phytoplankton community in Jiaozhou Bay as well, and the proportion of chrysophytes was greater in the inner Jiaozhou Bay(Liu et al., 2020b).

Two hundred and seven eukaryotic algae OTUs were obtained from 17 sediment samples from the Qingdao coast in this study, which was comparable to other metabarcoding studies on dinof lagellate cysts and phytoplankton resting stages (Dzhembekova et al., 2018; Liu et al., 2020a). However, compared with benthic microeukaryote studies in the Yellow Sea and European coastal areas (Gong et al., 2015; Massana et al., 2015), lower diversity of eukaryotic algae community was recorded in Qingdao coastal waters in this study. This study focused on the phytoplankton resting stages, and samples were stored at 4 °C for 6 months to inactivate most of the vegetative cells(Piredda et al., 2017). Only some parts of dinof lagellates, diatoms, chrysophytes, chlorophytes,haptophytes, and raphidophytes form resting stages,which can stay alive in sediments for a long time(reviewed by Ellegaard and Ribeiro, 2018). In addition, the inner Jiaozhou Bay had the lowest diversity and algal sequence number among the three sea areas, especially at station J3, in which only ca.400 algal reads were recorded. While eukaryotic sequence numbers measured in all the four samples from Jiaozhou Bay were more than 50 000 reads(Fig.2) with suffi cient sequencing depth.Eutrophication in Jiaozhou Bay has greatly increased recently, and the nutrient status has changed, which resulted in frequent algal blooms and decreased phytoplankton biodiversity (Shi et al., 2020). In addition, Jiaozhou Bay is an important sea area for shellf ish culturing, and the feeding of shellf ish reduces the diversity of eukaryotic microalgae to a certain extent (Yu et al., 2019).

4.3 Neglected taxa under routine microscopic observation

One remarkable diff erence between our study and the previous traditional phytoplankton and resting stage studies is that diverse phyla and classes were recorded in this study, including 24 classes of 8 phyla.However, diatoms and dinof lagellates were mostly recorded in the previous phytoplankton surveys based on traditional microscopic observations, and other phyla such as Chlorophyta, Cryptophyta, Haptophyta,Ochrophyta, Prasinodermatophyta, and Rhodophyta were rarely reported (Yang et al., 2014; Guo et al.,2019). Metabarcoding has a strong ability to detect the small-sized, morphologically similar, and the fragile and rare species (Salonen et al., 2019). A large number of small-sized species were detected in this study, for example, species in generaAzadinium,Biecheleria, andParagymnodinium(<20 μm),calcareous nannofossilsAchradinapulchra(<20 μm)andBraarudosphaerabigelowii(<5 μm), plenty of nano-sized chlorophytes and chrysophytes (<10 μm)(Supplementary Table S1). Liu et al. (2020b) analyzed phytoplankton community in Jiaozhou Bay based on metabarcoding, and detected 70 species which had never been reported in Jiaozhou Bay in the previous investigations, and majority of these species missed by previous studies are relatively small in size(<10 μm), mostly belonging to Cryptophyta,Chlorophyta, Ochrophyta, and Haptophyta, which are diffi cult to detect and identify using light microscopy.Meanwhile, Pigments from eight possible phytoplankton classes (Diatoms, Dinof lagellates,Chlorophyceae, Prasinophyceae, Chrysophyceae,Haptophyceae, Cryptophyceae, and Cyanobacteria)were detected in surface water by high performance liquid chromatography (HPLC) (Yao et al., 2010).These studies also suggested high diversities of phytoplankton community in Jiaozhou Bay.

Twelve taxa are f irstly recorded in the Chinese coastal waters in this study after reviewing the literatures based on both microscopic and metabarcoding studies (Supplementary Table S1).Many new recorded species occurred abundantly and widely in this study, and might be overlooked in previous microscopic studies due to the small size,fragile cell wall, and indistinguishable morphological characteristics. Furthermore, remarkable diverse chlorophytes were recorded in our study, including 34 OTUs in 10 classes. Most of these chlorophytes were recorded in previous phytoplankton metabarcoding studies (Xu et al., 2017; Song et al., 2019; Liu et al.,2020b). Three species in chlorophytes together with the other nine species have not been detected in previous phytoplankton metabarcoding studies in Chinese coasts (Xu et al., 2017; Song et al., 2019; Liu et al., 2020b), however were recorded in our study on sediments. It is well known that the resting stages in sediments are the assemblages of the deposition of phytoplankton resting stages over a period, while results of phytoplankton surveys ref lect the phytoplankton community in the water column at the sampling moment. Actually, many these new recorded species are resting stage forming or have the thick calcareous wall, which can be the fossil assemblages in sediments (Supplementary Table S1).

4.4 Occurrence of HAB species and potential risks of HABs in the Qingdao coast

Thirty-six potentially toxic/harmful and/or bloom species were detected in surface sediments from the Qingdao coast, among which ten species, includingGonyaulaxspinifera,Noctilucascintillans,Prorocentrummicans,Scrippsiellaacuminata,Aureococcusanophageff erens,Heterosigmaakashiwo,Chattonellamarina(Ministry of Ecology and Environment of the People's Republic of China,2010–2017),Chaetocerosdebilis(Yan et al., 2002),Margalef idiniumfulvescens(Lin et al., 2020), andNannochloropsisgranulate(Zhang et al., 2015), have been reported to form blooms in Jiaozhou Bay and the Qingdao coast before. Some of them are toxic species,which may produce toxins penitently harmful to human health. Most HAB species found in this study can form resting stages (29 to 36 species,Supplementary Table S1), and some of them occurred widely and abundantly in the sediments (Fig.9).Notably all of the 36 HAB species occurred in sediment samples from the Laoshan Coast. The result suggested the high potential risk of HABs in the Qingdao coastal area, especially in the Laoshan Coast.

Dinof lagellates were the most important HAB group identif ied in this study. Twenty-two species of dinof lagellates were identif ied to produce toxins and/or form blooms (Supplementary Table S1). The potential PSP producers,Alexandriumandersonii,Al.minutum, andGymnodiniumcatenatum(Lundholm et al., 2009 onwards), were present in sediments but in low DNA reads, indicating that the low historical abundance of these species in the water column. Liu et al. (2020b) did not identify any of these PSP producers in the phytoplankton survey in January 2019. The potential YTX producer,GonyaulaxspiniferaandProtoceratiumreticulatum(Lundholm et al., 2009 onwards), were detected in this study.Three OTUs were attributed toG.spinifera, two of which belonged to the top abundant OTUs (Fig.9).Blooms ofG.spiniferaoccurred at outside of Jiaozhou Bay in 2003 and 2007, respectively (Zhou et al.,2020). As a cyst-forming species,G.spiniferacan form variety types of cysts, mostly belonging to the cyst genusSpiniferiteswith wide morphological variations (Rochon et al., 2009). The great occurrence ofG.spiniferain the Laoshan Coast (Fig.9) indicated high numbers of its vegetative cells in the water column and the high risk of its blooms. AnotherGanyaulaxspecies,G.fragiliswhich has been associated with the massive mucilage formations in the Marmara Sea, Turkey (Balkis et al., 2011), rarely distributed at one station of the Laoshan Coast in this study.

Three species ofAzadiniumwere detected in this study, onlyAz.poporumhas the potential to produce AZP (Krock et al., 2012).Az.poporumwas reported to distribute widely in surface sediments from the Chinese coasts (Gu et al., 2013b, Liu et al., 2020a),and 13 out of 16 strains ofAz.poporumfrom diff erent geographic locations along the Chinese coastline contained azaspiracids (Krock et al., 2014). DNA sequence number ofAz.poporumwas relatively low,on the other hand, the nontoxicAz.trinitatumandAz.dalianensewidely distributed in high DNA reads in this study (Fig.9).

The karlotoxin producer,Karlodiniumvenef icum,occurred in most samples (Fig.9).K.venef icumis a small sized naked donif lagellate, which produces a groups of polyketide toxins (karlotoxins) exerting hemolytic, ichthyotoxic, and cytotoxic eff ects on f ishes and aquatic animals (Place et al., 2012). Since the f irst bloom ofK.venef icumreported in the South China Sea in 2003 (Liu et al., 2020c), its blooms occurred in the Chinese coastal sea areas frequently(Wang et al., 2011; Dai et al., 2014). Cysts of this species were detected widely in sediments from the China Sea as well (Liu et al., 2020c). Considering the association of its blooms and massive f ish killings(Place et al., 2012), the wide distribution ofK.venef icumis worthy of attention.

Biecheleriawas rarely observed in the previous routine phytoplankton surveys due to the small size and frangible thin wall.Biecheleriahas become more frequently detected in the phytoplankton communities as the development of molecular biological techniques(Sundström et al., 2010). Moreover, taxa inBiecheleriagenerally form resting cysts (Moestrup et al., 2009),which made them as the common dominant eukaryotes in sediments (Dzhembekova et al., 2018; Liu et al.,2020a; Rhodes et al., 2020). Ten OTUs ofBiecheleriawere detected in sediments from the Qingdao coast,andB.halophilawas the f irst dominant species detected in this study andB.tirezensiswas newly recorded in the Chinese coastline.B.halophilahas been reported to form a co-occurring bloom withScrippsiellahangoeiin the Baltic Sea (Kremp et al.,2005). The results suggested that species inBiecheleriawere probably of much wider distribution,however might be ignored during the microscopic observation because of their fragility and small size.

Other ichthyotoxic species, whose blooms have been associated with massive f ish kills, included the dinof lagellatesPolykrikoshartmannii, andPo.geminatus, the pelagophyteAureococcusanophageff erens, and the raphidophytesChattonellamarina,Fibrocapsajaponica, andHeterosigmaakashiwo(Supplementary Table S1). Among them,blooms ofAu.anophageff erens,Chat.marina, andH.akashiwooccurred in Jiaozhou Bay and the nearby sea areas (Ministry of Ecology and Environment of the People’s Republic of China, 2010–2017), and brown tides caused byAu.anophageff erenshave recurred frequently in the northwestern Bohai Sea since 2009 (Huang et al., 2020).Po.geminatusis a common bloom species in the Pearl River Estuary of the South China Sea, and its blooms have frequently occurred since 2006 (Shen et al., 2012; Dong et al.,2020).Po.geminatuswas detected in sediment samples from eight stations in this study, and ranked the top 12thabundant OTU (Fig.9). To the best of our knowledge, it is the f irst report of this species present in the northern Chinese coastline, i.e., the Yellow sea and the Bohai Sea. Coincidentally, another bloom species in the Pearl River Estuary,Levanderinaf issa,occurred widely and abundantly in this study (Fig.9).

4.5 Potential of metabarcoding for microalgal resting stage monitoring

Microscopic observations require professional taxonomic expertise for distinguishing the tiny morphological diff erences among species from an assemblage of organisms (Shang et al., 2019).Furthermore, only a small proportion of phytoplankton resting stages are known for their corresponding vegetative cells, and vice versa. Misidentif ication is therefore common and almost unavoidable due to these diffi culties (Gómez et al., 2017; Shang et al.,2019). Metabarcoding can avoid these problems, and standardized protocols can be developed to make results more comparable among laboratories (Pochon et al., 2015). Metabarcoding has the ability to detect small-sized taxa that were not observed in previous microscopic surveys, such as species inAzadinium,Biecheleria, andParagymnodinium, calcareous nannofossilsAchradinapulchraandBraarudosphaerabigelowii, and the nano-sized chlorophytes in this study. The results suggest that metabarcoding provides comprehensive insights into benthic communities, especially for organisms that are diffi cult to identify because of their small size, cryptic habits (such as the parasitic and symbiotic species),occurrence in the form of resting stages or lack of traditional identif ication keys (Salonen et al., 2019).

The f irst problem in sediment metabarcoding analysis is the less effi ciency of genomic DNA extraction for resting stages than for vegetative cells,due to the presence of the thick resistant walls. In addition, DNA extraction is particularly diffi cult from sediment samples containing high concentrations of clay and humic material. Therefore, extraction and purif ication of DNA are crucial steps in the analysis of marine sediments, where the presence of high organic matter, organic pollutants and heavy metals can interfere with the recovery of DNA and reduce the eff ectiveness of PCR (Piredda et al., 2017). This might be part of the reasons the lower diversity in the sediment sample than in phytoplankton. We applied the standardized protocol suggested by DNeasy PowerSoil Kit for DNA isolation, which included complete vortex of the soil samples in the PowerBead tube with Solution C1 to break the tough cells such as resting stages. The kit has been widely used and shown to be well optimized for DNA extraction from a variety of soils and plant tissues (reviewed by Lear et al., 2018). Actually, high diversity of eukaryotic algae (207 OTUs) was obtained in this study, and the cyst/spore-forming taxa such as the dinof lagellatesBiecheleria,Scrippsiella,Gymnodinium, the ChrysophyteParaphysomonas, and the diatomChaetoceroswere the dominant ribotypes in the sediments. The results indicated that DNA was eff ectively extracted from these thick-wall spores/cysts, and thus suggested the effi ciency of DNA extraction.

Our samples were stored in a cold dark condition for six months in order to ensure that most survives in the sediments were resting stages (Piredda et al., 2017).About half taxa that were identif ied to the species level in this study have been reported to form resting stages(Supplementary Table S1), and most of the top abundant OTUs were reported to form resting stages (Fig.9).Other taxa such asParaphysomonaslucasi,Gonyaulaxsp., andChaetocerossp. occurred abundantly in this study (Fig.9). However, whether they could form resting stages was still unknown. Indeed, cyst/spore formation is an important characteristic in genusParaphysomonas(Coradeghini and Vigna, 2008),Gonyaulax(Matsuoka and Fukuyo, 2000), andChaetoceros(Montresor et al., 2013). Microscopic observations were undertaken to observe dinof lagellate cysts in the same sediments synchronously, and 23 taxa of dinof lagellate cysts were observed, and only 13 taxa were identif ied to the species level (Lei, 2021). While 127 dinof lagellate OTUs were detected by metabarcroding in this study, and 66 OTUs were attributed to 45 species of dinof lagellates, including 32 cyst-forming species (Supplementary Table S1). The results suggested much higher biodiversity of dinof lagellate cysts was revealed by metabarcroding than that by microscopy. However, we did not analyze resting stages of the other phytoplankton groups by microscopy due to the lack of morphological references.Those taxa that do not produce resting stages were generally detected in low DNA reads, some of which might be the benthic vegetative cells, and some of which might represent DNA fragments of dead cells or even DNA relics. Actually, it is hard to distinguish the dead and partially degraded phytoplankton cells from authentic sediment inhabiting cells or resting cells during HTS based metabarcoding, which may increase the OTU richness to some extent. However, DNA fragments would be presented as extremely rare OTUs after long time storage under cold dark condition(Shang et al., 2019). The rare OTUs with 1–2 reads(singletons and doubletons) were removed in this study,which reduced the possibility of increased OTU richness due to DNA fragments. On the other hand,storing sediment samples for six months could turn benthic microalgae into resting stages, which might increase OTU richness and diversity of microalgal resting stages in our study. However, majority of cyst/spore forming taxa are phytoplankton species (reviewed by Ellegaard and Ribeiro, 2018), and only a few benthic algae have been reported to form resting stages up until now (Limoges et al., 2015; Mertens et al., 2017).Furthermore, most of the species identif ied in this study were marine phytoplankton (Supplementary Table S1).Therefore, the benthic algae had little contribution to the diversity in this study.

The choice of the targeted region has a potentially signif icant inf luence in metabarcoding surveys. The V4 and V9 regions of 18S rDNA are two of the most commonly used regions for metabarcoding eukaryotic organisms. The V4 region is the longest 18S region in eukaryotes and features high variability, which increases its resolution relative to other 18S rRNA regions (reviewed by Lear et al., 2018). The V4 region is one of the best choices for metabarcoding the entire eukaryotic community from environmental DNA(reviewed by Lear et al., 2018) and has been used in many studies focusing on marine microalgae (Kataoka et al., 2017; Penna et al., 2017; Xu et al., 2017; Song et al., 2019; Liu et al., 2020b; Cui et al., 2021). Gaonkar et al. (2018) compared the 18S rDNA V4 and V9 hypervariable region sequences of taxa in the planktonic diatom family Chaetocerotaceae, and found that the 18S V4 region f it well with target taxa,and the V4-tree topology was similar to that of the whole 18S rDNA. The primer sites used in our study produce amplicons of lengths (450 bp) just right ideal length for sequencing on the Illumina MiSeq platform.Results from the phylogenetic tree showed that OTUs in the same phylum were generally clustered together(Fig.5), suggesting the credibility of species annotation.

In addition, the reference database used potentially has a large impact on the eukaryote community obtained. Of the 207 generated algal OTUs, only 124 OTUs (59.9%) were identif ied at the species level, and some unclassif ied OTUs were the common predominant sequences in our study. The lack of coverage at the species level for many taxa may hinder correct interpretation of metabarcoding data(Forster et al., 2016). However, as new sequencing data are continuously contributing to the existing databases, the situation is likely to improve in the future. On the other hand, the identif ication and enumeration of microalgae based on microscopic methods are restricted to organisms with welldocumented morphological characteristics (Piredda et al., 2017). While the morphologies of resting stages of most microalgae are still unknown, and it is diffi cult to analyze them under microscopy. The capacity of HTS to detect morphologically cryptic, small-sized,and rare taxa makes this method a suffi ciently sensitive detection approach for assessing the composition of resting stages in sediments.

5 CONCLUSION

This study enriched our understanding of eukaryotic algae in surface sediments from the Qingdao Coast based on 18S rDNA metabarcoding.Diverse eukaryotic algae were recorded, including 207 OTUs belonged to 75 genera in 24 classes of 8 phyla. Within 98 taxa identif ied to the species level,12 species were f irstly recorded in the Chinese coastal waters. Dinophyta generally dominated in both DNA reads and OTU richness. The benthic algal community structures were quite diff erent in three sea areas of the Qingdao coast, which was dominated by chrysophytes in Jiaozhou Bay, by dinof lagellates in the Laoshan Coast, and co-dominated by dinof lagellates and chrysophytes in the Amphioxus Reserve area.Furthermore, eukaryotic algae OTU richness and diversity were higher in the Laoshan Coast, and lower in Jiaozhou Bay. A large number of small-sized species were detected in this study, such asAzadiniumspp.,Biecheleriaspp., andParagymnodiniumspp.,and nano-sized chlorophytes and chrysophytes. The results suggested the strong ability of metabarcoding to detect the small-sized, morphologically similar,and the fragile and rare species. Notably, 36 HAB species were detected in this study, and some of these species occurred widely and dominantly in this study,suggesting high potential risk of HABs in the Qingdao coastal area.

6 DATA AVAILABILITY STATEMENT

The data in this study are available from the National Center for Biotechnology Information(NCBI) Sequence Read Archive (SRA) (http://www.ncbi.nlm.nih.gov/Traces/sra) with the accession number PRJNA722324.