61urn:lsid:arphahub.com:pub:F9A654F3-34B8-5338-9F9C-D38A394E7BA2Metabarcoding and MetagenomicsMBMG2534-9708Pensoft Publishers10.3897/mbmg.6.7818178181Research ArticleAlgaeBacteriaCyanobacteriaProtozoaDNA barcodingPopulations & CommunitiesJapanPacific OceanAnnual dynamics of eukaryotic and bacterial communities revealed by 18S and 16S rRNA metabarcoding in the coastal ecosystem of Sagami Bay, JapanSogawaSayaka1ssogawa@affrc.go.jpConceptualizationData curationFormal analysisFunding acquisitionInvestigationVisualizationWriting - original draftTsuchiyaKenji2Formal analysisVisualizationWriting - original draftNagaiSatoshi1Data curationFunding acquisitionMethodologySoftwareWriting - review and editingShimodeShinji3InvestigationResourcesWriting - review and editingKuwaharaVictor S.4Funding acquisitionInvestigationProject administrationSupervisionWriting - review and editingJapan Fisheries Research and Education Agency, 2-12-4 Fukuura, Kanazawa, Yokohama, Kanagawa 236-8648, JapanJapan Fisheries Research and Education AgencyYokohamaJapanNational Institute for Environmental Studies, 16-2 Onogawa, Tsukuba, Ibaraki 305-8506, JapanNational Institute for Environmental StudiesTsukubaJapanManazuru Marine Center for Environmental Research and Education, Graduate School of Environment and Information Sciences, Yokohama National University, 61 Iwa, Manazuru, Kanagawa 259-0202, JapanYokohama National UniversityManazuruJapanGraduate School of Science & Engineering, Soka University, 1-236 Tangi, Hachioji, Tokyo 192-8577, JapanSoka UniversityHachiojiJapan
2022280220226e78181BA40649D-F3E0-5CBE-8A7C-14C5F83079481911202114022022Sayaka Sogawa, Kenji Tsuchiya, Satoshi Nagai, Shinji Shimode, Victor S. KuwaharaThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Sagami Bay, Japan is influenced by both the warm Kuroshio Current and the cold Oyashio Current and rich nutrients are supplied from multiple river sources and the deep-sea, forming a dynamic ecosystem. The aim of the present study was to investigate eukaryotic and bacterial communities in the coastal waters of Sagami Bay, using 16S rRNA and 18S rRNA sequencing and to assess the seasonal and vertical dynamics in relation to physicochemical and biological conditions. Eukaryotic and bacterial communities showed synchronous seasonal and vertical changes along with environmental variability. Diversity of plankton community suspended in the surface was lower than those at the subsurface layers in both the eukaryotes and bacteria communities; however, community diversity showed different characteristics in the subsurface where the eukaryotic community at the deeper layer (100 m) was as low as the surface and highest in intermediate depth layers (10–50 m), while that of bacterial community was highest in the deeper layer (100 m). The annual variability of the coastal microbial communities was driven, not only by the seasonal changes of abiotic and biotic factors and short-term rapid changes by river water inflow and phytoplankton blooms, but also largely influenced by deep-seawater upwellings due to the unique seafloor topography.
Marine pelagic ecosystems are interconnected by two types of food chains: grazing food chains and microbial food chains (microbial loop). Classical grazing food chains begin with phytoplankton photosynthesis (producers), which become food for zooplankton, such as herbivorous copepods and further preyed upon by larger organisms, such as fishes (Fenchel 1988; Azam 1998). The microbial loop is a relatively new food chain concept beginning with bacteria introducing dissolved organic material (DOM) to the classical food chain via heterotrophic protist and microzooplankton (Pomeroy 1974; Azam et al. 1983; Pomeroy et al. 2007). In the marine ecosystem, physical factors, such as currents, temperature, density and light, chemical factors, such as nutrients and oxygen and biological factors related directly or indirectly with marine organisms, regulate the system. Several studies have shown that diversity, abundance and the community structure of marine plankton correspond with changes in physical and chemical factors and interaction between marine plankton (Sin et al. 1999; Hays et al. 2005; Sogawa et. al. 2013; Signori et al. 2014; Vajravelu et al. 2018; Ibarbalz et al. 2019). Further, interactions between primary producers and bacteria influence one another, having significant impact on biogeochemical cycles and form the basis of marine ecosystem community and diversity (Azam 2007; Little et al. 2008; Amin et al. 2015). The dynamics of coastal marine microbial communities are also affected by seasonal change in abiotic and biotic factors and short-term effects of river water and phytoplankton blooms (Lucas et al. 2015). Thus, a comprehensive understanding of plankton communities in marine ecosystems and revealing changing mechanisms are essential for sustainable use and effective management of marine resources (Ministry of Environment of Japan 2011; Sogawa et al. 2018). However, it is difficult to grasp the diverse plankton community by morphological identification using microscopy as many of the marine plankton are small and amorphous (Herndl and Peduzzi 1988; Bochdansky and Herndl 1992; Hirai et al. 2017) and require advanced techniques and experience for identification (Abad et al. 2016, 2017). The advances of molecular biological techniques, such as DNA metabarcoding using high-throughput sequencing, have become an efficient and powerful tool for monitoring and assessing marine communities (Valentini et al. 2016; Ruppert et al. 2019; Suter et al. 2020).
Sagami Bay, located on the Pacific Ocean side of middle Honshu Island, Japan, is known for its variety of aquatic landform, containing rivers (terrestrial runoff), coastal zones and deep seas (deeper than 1,000 m). Abundant nutrients are supplied from multiple river sources and deep ocean upwelling, forming rich harvest grounds containing a wide variety of migratory and local fishery catch. The Bay is influenced by both warm Kuroshio and cold Oyashio currents: Kuroshio, the western boundary current of the North Pacific subtropical gyre, transports warm and saline subtropical waters as well as a diverse assortment of organisms and Oyashio, the western boundary current of the western subarctic gyre, branches and flows from the opposite direction of Kuroshio in the mesopelagic layer and flows into Sagami Bay. Consequently, the Bay would be an appropriate location for monitoring diverse plankton communities and survey of the area could provide indispensable information towards understanding the coastal aquatic ecosystem. Previous studies reported temporal/seasonal variation of bacterio-, phyto- and zooplankton in Sagami Bay and their ecology in relation to environmental factors. (Takasuka et al. 2003; Ara et al. 2006; Miyaguchi et al. 2006; Shimode et al. 2006; Hashihama et al. 2008; Baek et al. 2009; Baki et al. 2009; Tsuchiya et al. 2013, 2016; Sugai et al. 2016, 2018, 2020). A few studies analysed diversity of plankton using 18S and 16S rRNA gene (Kita-Tsukamoto et al. 2006; Kok et al. 2012a, b). However, no studies have conducted comprehensive seasonal analysis on plankton communities in Sagami Bay that address environmental change and interactions between organisms. In this study, a continuous survey over a year was conducted at a long-term monitoring site of Sagami Bay and 18S-RNA and 16S-rRNA genes were analysed to reveal the annual vertical dynamics of eukaryotic and bacterial communities in the temperate coastal pelagic ecosystem.
MethodsField sampling and DNA extraction
Seasonal observations were carried out at the shelf station (120 m depths, latitude 35°09.75'N, 139°10.00'E longitude), located in the western part of Sagami Bay (Fig. 1). Continuous surveys were conducted 23 times over a one-year period from October 2018 to February 2020 (Table 1) on board R/V “Tachibana” of the Manazuru Marine Center for Environmental Research and Education (MMCER), Graduate School of Environment and Information Sciences, Yokohama National University, Yokohama National University.
Survey dates in Sagami Bay. Asterisks represent sampling only on the surface.
Year
Month and day
2018
Oct (18th), Nov (15th), Dec (13th)
2019
Feb (14th), Mar (17th), Apr (19th), May (17th, 24th* and 31st*), Jun (14th), Jul (9th*, 12th and 31st*), Aug (6th*), Sep (30th*), Nov (15th and 25th*), Dec (10th* and 20th)
2020
Jan (6th* and 31st*), Feb (4th* and 20th)
085708FC-66C5-5B89-BAD1-1A33458C8783
Schematic illustration of surface current around Japan (left) and location of survey site (right, red circle) in Sagami Bay, Japan.
https://binary.pensoft.net/fig/651865
Vertical profiles of temperature, salinity, density, chlorophyll a fluorescence and dissolved oxygen were measured continuously using a RINKO-Profiler (ASTD102, JFE Advantech Co., Ltd). Mixed layer depths (MLDs) were calculated from density with the following equation: MLD: Δσt(z) = σt(z) – σt(5m) > 0.125 kg/m3 (Monterey and Levitus 1997). The euphotic zone (1% light level; EZ1%) was measured by Profiling Reflectance Radiometer (PRR) / Compact Optical Profiling System (C-OPS). Seawater samples were collected at several depth layers (maximum depth 100 m) in Niskin bottles and directly from the surface layer in sampling bottles. Seawater samples of 1litre for 16S and 18S rRNA analysis were filtered through Sterivex filter (0.22 μm) and DNA was extracted using 5% Chelex buffer (Nagai et al. 2012; Tanabe et al. 2016). Sampling bottles and filtering equipment were washed with a 10% commercial bleach solution and rinsed with distilled water before usage. Seawater samples for discrete vertical profiles of nutrient (nitrate + nitrite, silicate, phosphate) concentrations were collected at the same depths and measured using a continuous flow injection analyser (QuAAtro, BL-Tech).
For detection of eukaryotic and bacterial communities, respectively primer pairs of V7–9 region in 18S-rRNA gene (18S-V7F: TGGAGYGATHTGTCTGGTTDATTCCG and 18S-V9R: TCACCTACGGAWACCTTGTTACG; Sildever et al. 2019) and V3–4 region in 16S-rRNA gene (16S-341F: CCTACGGGNGGCWGCAG and 16S-805R: GACTACHVGGGTATCTAATCC) (V3-V4 region, Sinclair et al. 2015) were used as universal primers for metabarcoding analysis. Two-step PCR for the construction of paired-end libraries and HTS on Illumina Miseq 300 PE platform (Illumina, San Diego, CA, USA) followed the protocol described by Dzhembekova et al. (2017). The first PCR was performed using a thermal cycler (PC-808; ASTEC, Fukuoka, Japan) in a 25 μl reaction mixture containing 1.0 μl template DNA (< 1 ng), 0.2 mM of each dNTP, 1× PCR buffer, 1.5 mM Mg2+, 1.0 U KOD-Plus-ver.2 (Toyobo, Osaka, Japan) and 1.0 μM of each primer. PCR conditions were first denaturation at 94 °C for 3 min, followed by 30–32 cycles (18S-rRNA gene) or 28 cycles (16S-rRNA gene) at 94 °C for 15 s, 56 °C for 30 s and 68 °C for 40 s. PCR amplification was verified by 1.5% agarose gel electrophoresis. PCR products were purified using Agencourt AMPure XP (Beckman Coulter, Brea, CA, USA) and eluted in 25 μl TE buffer according to the manufacturer’s protocol. The second PCR was performed using the same conditions for the first PCR, except the reaction mixture volume was 50 μl with the addition of 2.0 μl diluted PCR product from the first round as a template. The PCR cycling conditions were as follows: initial denaturation at 94 °C for 3 min, followed by 10–12 cycles at 94 °C for 15 s, 56 °C for 30 s and 68 °C for 40 s. PCR amplification was again verified by agarose gel electrophoresis and the PCR products were purified using an Agencourt AMPure XP (Beckman Coulter, Brea, CA, USA). The amplified PCR products were quantified, mixed equally and stored at –30 °C until used for sequencing.
Treatment processes and operational taxonomic unit picking
Nucleotide sequences were demultiplexed depending on the 5′-multiplex identifier tag and primer sequences according to the default format in MiSeq. Sequences longer than 300 bp were truncated to 300 bp by trimming the 30 tails. The trimmed sequences shorter than 250 bp were filtered out. Demultiplexing and trimming were performed using Trimmomatic version 0.35. The remaining sequences were merged into paired reads using Usearch version 8.0.1517. In addition, singletons were removed. Sequences were then aligned using Clustal Omega version 1.2.0. Multiple sequences were aligned with each other and only sequences that were contained in more than 75% of the read positions were extracted. Filtering and part of the multiple alignment process were performed as described in the Miseq standard operating procedure using Mothur (Schloss et al. 2011). Erroneous and chimeric sequences were detected and removed using Mothur (Edgar et al. 2011). Demultiplexed and filtered, but untrimmed sequence data are deposited into the DDBJ Sequence Read Archive (access no. DRA013426). The sequences were clustered to select representative sequences at a 0.99 level of sequence identity by CD-HIT-EST version 4.6.8 (Li and Godzik 2006).
Taxonomic identification of the OTUs
The sequence database used to identify OTUs was downloaded from GenBank on 20 August 2020. The taxonomic identification of each OTU was performed by the BLAST search (Cheung et al. 2010) with NCBI BLAST+ 2.10.1+ (Camacho et al. 2009). The taxonomic information was obtained from the BLAST hit with the top bitscores for each query sequence. The OTUs of the same top-hit were merged and the OTUs were also merged when an OTU matched several data with the same bitscore and top-hit similarity. Several data, therefore, were sometimes merged at an OTU. The removal of sequences containing errors was imperfect after the successive processes of MPS data treatment; sequences containing different types of errors derived from the original ones remained in the following analytical steps. Therefore, these sequences were detected as artificial OTUs with the same top BLAST hit name, but with slight differences. To avoid overestimation of the OTUs, these artificial OTUs were merged into a single OTU with the greatest similarity score. The phylum and class/order (prokaryotes) level taxonomic grouping were also confirmed by the feature-classifier command in Qiime2 (https://docs.qiime2.org/2020.8/) (Bokulich et al. 2018) using the Silva database (https://docs.qiime2.org/2020.8/data-resources/common/silva-138-99-nb-classifier.qza) (Yilmaz et al. 2014) as the taxonomic database.
Statistical analysis
All multivariate analyses of eukaryotic community structure and diversity were performed using PRIMER version 7 with PERMANOVA+ add-on software (Anderson et al. 2008; Clarke and Gorley 2015). For multivariate analyses, Bray-Curtis similarities (Bray and Curtis 1957) were calculated amongst samples, based on the log-transformed abundance data of sequence reads. The OTU hit as Chloroplast was excluded from the 16S analysis. A non-metric multidimensional scaling (NMDS) ordination was used to visualise the differences in the community amongst samples and analysis of similarity (ANOSIM) tested differences amongst water depths and months (Clarke 1993). A distance-based linear model (DISTLM) and distance-based redundancy analysis (DBRDA) were conducted to explain changes in eukaryotic and bacterial community with respect to environmental variables. Environmental data were normalised and square root transformed prior to the analysis. Multicollinearity was tested prior to DISTLM and DBRDA and variables with correlations stronger than 0.95 were removed. Plankton diversity was measured according to the total number of OTUs together with the Shannon-Wiener Diversity Index (Shannon and Weaver 1949) and the Evenness (Pielou 1966).
ResultsPhysical and chemical oceanographic environment
Vertical profiles of temperature, salinity, density, chlorophyll fluorescence, dissolved oxygen and nutrients in Sagami Bay during the observation period are shown in Fig. 2. Water temperature in the sampling layer ranged from 26.9 °C at the surface in August 2019 to 13.9 °C at 100 m depth layer in April 2019. Salinity varied from 34.6 to 32.8 and the highest value was observed at 100 m depth layer in February 2019 and the lowest value was observed at the surface in July 2019. Density was highest at 100 m depth layer in April 2019 (25.8 σt) and lowest at the surface in July (21.2 σt). The mixed layer depths (MLDs) were the deepest in February when vertical mixing was strong and < 20 m in May to November when the seasonal thermocline was well developed. Chlorophyll a varied from 2.24 μgl-1 to 0.11 μgl-1 with highest at 10 m depth layer in April 2019 and the lowest at 100 m depth layer in October 2018. Nutrient concentration ranged from 14.53 μM to 0.04 μM for nitrate + nitrite and 1.012 μM to 0.054 μM for phosphate with the maximum value observed at 100 m depth layer in April 2019 and the minimum value at 10 m depth in May 2019 for nitrate + nitrite and in July for phosphate. The maximum value of silicate was observed at the surface in July 2019 (22.60 μM) and the minimum value at the surface in September 2019 (0.98 μM).
BC028AF4-A9F4-5F9A-B9AF-75C10DD0CF90
Vertical profiles of temperature, salinity, density, chlorophyll fluorescence, dissolved oxygen (DO) and nutrients in Sagami Bay during observation period. Crosses show water sampling layer. The Mixed Layer Depths (MLDs), calculated
https://binary.pensoft.net/fig/651866Eukaryotic and bacterial community diversity
In the study area, 1,660 OTUs (985 OTUs from the surface) of eukaryotic plankton were detected. Supergroup Alveolata was the largest of 658 OTUs (377 OTUs from the surface), followed by Stramenopiles of 368 OTUs (259 OTUs), Opisthokonta of 262 OTUs (131 OTUs), Rhizaria of 160 OTUs (77 OTUs), Hacrobia of 90 OTUs (66 OTUs) and Viridiplantae of 77 OTUs (50 OTUs). The number of OTUs detected in the study area were ca. three times larger in the bacterial community (4477 OTUs and 2838 OTUs from the surface). Phylum Proteobacteria was the largest of 2483 OTUs (1573 OTUs from the surface), followed by Bacteroidota of 1067 (786 OTUs), Verrucomicrobiota of 112 (66 OTUs), Actinobacteriota of 87 (70 OTUs), Bdellovibrionota of 78 OTUs (31 OTUs) and Cyanobacteria of 71 (56 OTUs).
Diversity analysis was conducted for eukaryotic community amongst depths (Fig. 3). The mean numbers of OTUs were lower in the surface (175) and the bottom depth layer (100 m; 169) and highest in the intermediate depth layer (30 m; 211). Pielou’s Evenness was the lowest in the surface (J′ = 0.64) and higher in the intermediate depth layers (10–50 m; J′ = 0.69–0.71). The Shannon-Wiener Diversity Index was also lower in the surface (H′ = 3.3) and the bottom layer (H′ = 3.4) and higher in the intermediate depth layers (10–50 m; H′ = ca. 3.7). The number of OTUs were higher in intermediate depths layers (10–50 m) and lower in the surface and bottom depth layer, although rarefaction curves of all depths did not reach asymptote (Fig. 4).
C768686A-EE05-50B3-AB1E-FA7EC9B7A8E8
Species richness, Shannon-Wiener Diversity (H′) and Pielou’s Evenness (J′) of eukaryotic and bacterial communities in Sagami Bay comparing amongst sampling depth layers.
Sample-based rarefaction curves of each sampling depth layer of eukaryotic and bacterial community.
https://binary.pensoft.net/fig/651868
Diversity analysis was also conducted for bacterial community amongst depths (Fig. 3). Characteristics of the mean numbers of OTUs were contrary to the eukaryotic community that were higher in the surface (321) and the bottom depth layer (100 m; 318) and lower in the intermediate depth layer (30 m and 50 m; 282 and 289, respectively). Pielou’s Evenness was the lowest in the surface (J′ = 0.61), the same as the eukaryotic community; however, they gradually increased with increasing depth and reached the highest in the bottom depth layer (100 m; J′ = 74). The Shannon-Wiener Diversity Index was also lowest in the surface (H′ = 3.5) and highest in the bottom depth layers (100 m; H′ = 4.2). Rarefaction curves of all depths have not reached asymptote, but were higher in the surface and lower in the intermediate depth layers (30 m and 50 m), contrary to the eukaryotic community (Fig. 4). Diversity of plankton community suspended in the surface was lower than those at the subsurface layers in both the eukaryotes and bacteria communities; however, the diversity of the eukaryotic community at the bottom layer (100 m) was as low as the surface, while that of the bacterial community was highest, contrary to the surface.
Annual dynamics of eukaryotic and bacterial community
The ANOSIM test revealed significant differences amongst months for both eukaryotic and bacterial community (global test: R = 0.395, P = 0.001 for 18S and R = 0.701, P = 0.001 for 16S). On the contrary, there were no significant differences amongst depth layers for both the eukaryotic and bacterial communities (ANOSIM global test: R = 0.009, P = 0.389 for 18S and R = 0.005, P = 0.470 for 16S). NMDS ordination represented samples collected in similar seasons (months) and depth layers are plotted closer, particularly during winter (Fig. 5). In winter, plankton communities were more vertically homogeneous as samples are plotted closer compared to spring to autumn.
FD5B8E53-EC60-5F5B-9F25-C1420987B8A9
Non-metric multidimensional scaling (NMDS) ordination of eukaryotic (upper figures) and bacterial communities (lower figures) in Sagami Bay. Coloured area showed plot ranges of each depth layer (left figures) and each month (right figures). Numbers in parenthesis after month represent year of sampling; “18” as 2018, “19” as 2019 and “20” as 2020.
https://binary.pensoft.net/fig/651869
Distance-based multivariate linear model (DISTLM) analysis revealed significant relationships (p < 0.01) between eight out of nine variables and eukaryotic and bacterial community structures, respectively (Table 2). One variable, phosphate, was excluded prior to analysis as it was strongly correlated (|r| > 0.95) with nitrate + nitrite. MLDs accounted for the largest proportion of fitted variance for eukaryotic community, followed by nitrate + nitrite, density, chlorophyll a, EZ1% and temperature (11.5–4.2%) (Table 2). Axis 1 of the distance-based redundancy analysis (DBRDA) accounted for 46.3% of the fitted model (11.9% of the total variation) and strongly correlated with MLDs (Fig. 6). Axis 2 of DBRDA accounted for 28.6% of the fitted model (7.3% of the total variation) and strongly correlated with nitrate + nitrite.
Percentage of variation explained in a distance-based multivariate linear model (DISTLIM) of eukaryotes and bacterial community.
Marginal Tests (18S)
Marginal Tests (16S)
Variable
p
Variation explained (%)
Variable
p
Variation explained (%)
MLD
0.0001
11.48
MLD
0.0001
11.60
Nitrate + nitrite
0.0001
8.25
Density
0.0001
9.95
Density
0.0001
7.44
EZ1%
0.0001
9.75
Chlorophyll a
0.0001
6.46
Nitrate + nitrite
0.0001
9.31
EZ1%
0.0001
6.45
Chlorophyll a
0.0001
8.71
Temp
0.0001
4.25
Month
0.0001
8.04
Month
0.0003
3.82
Temp
0.0001
5.15
Silica
0.0009
3.45
Silica
0.0002
4.75
Sal
0.5616
1.07
Sal
0.4044
1.37
54B20DDB-D5DB-53F8-A29F-DEC43125902F
Distance-based redundancy analysis (DBRDA) ordination of eukaryotic and bacterial community. The distance-based redundancy analysis was constrained by the best-fit explanatory variables from a distance-based multivariate linear model (DISTLM).
https://binary.pensoft.net/fig/651870
In the top 10 OTU of eukaryotic communities, Amobophyra (Alveolata, Dinophyceae) and Astrosphaera (Rhizaria, Polycystinea), which were abundant at 100 m depths, were positively related to nitrate + nitrite and negatively related to chlorophyll a (Fig. 6). On the contrary, Paracalanus (Opisthokonta, Metazoa), which were abundant on the surface (0 m), was positively correlated to chlorophyll a. Karlodiniumveneficum, Heterocapsarotundata (Alveolata, Dinophyceae) and Leptocylindrusconvexus (Stramenopiles, Bacillariophyta), which were abundant in spring-autumn, were negatively correlated to MLDs. On the other hand, Ostreococcus, Bathycoccusprasinos and Micromonaspusilla (Viridiplantae, Chlorophyta), which were abundant in winter, were positively correlated to MLDs.
MLDs also accounted for the largest proportion of fitted variance for bacterial community, followed by density, EZ1%, nitrate + nitrite and chlorophyll a (11.6–8.7%) (Table 2). Axis 1 of the DBRDA accounted for 33.6% of the fitted model (14.3% of the total variation) and strongly correlated with MLDs, EZ1%, density and chlorophyll a (Fig. 6). Axis 2 of DBRDA accounted for 24.7% of the fitted model (10.5% of the total variation) and strongly correlated with density and MLDs. These two axes were strongly related to water mass structure and seasonality.
In the top 10 OTU of bacterial communities, Flavobacteriales, Cellvibrionales and Rhodobacterales, which were abundant in spring and summer, were positively related to chlorophyll a and negatively related to MLDs (Fig. 6). Synechococcus (CC9902) was located in the second quadrant in summer and autumn and positively correlated to temperature and negatively correlated to density. SAR11 subclade IV was abundant in autumn, where chlorophyll a was relatively low and negatively related to Flavobacteriales and Rhodobacterales. SAR11 subclade Ia was negatively correlated to Flavobacteriales and Rhodobacterales, as well as SAR11 subclade IV.
Cluster analysis of bacterial community
Bacterial communities were classified into four clusters at 30% of the Bray-Curtis Similarity in cluster analysis (Fig. 7, Suppl. material 1: Fig. S1). Moreover, clusters of 3 and 4 were further classified into 3 (3-1, 3-2 and 3-3) and 4 groups (4-1, 4-2, 4-3 and 4-4) at 39% and 41% of the Bray-Curtis Similarity, respectively (Fig. 7, Suppl. material 1: Fig. S1). In cluster 1, limited in deep samples, SAR11_clade was most dominant (14 ± 8%), followed by Flavobacteriales (8.7 ± 2.4%), Sphingomonadales (6.9 ± 6.2%) and Oceanospirillales (6.2 ± 3.2%) (Fig. 8, Suppl. material 1: Fig. S2). SAR324_clade and SAR406_clade were also relatively abundant (5.0 ± 1.4% and 4.9 ± 1.6%, respectively). Cluster 2 consisted of surface samples from July (190709) to November (191115). In this cluster, Favobacteriales (32.7%), Rhodobacterales (11.6%), SAR11_clade (10.8%), Synechococcales (10.6%) and Cellvibrionales (11.7%) were dominant. In cluster 3-1, Synechococcales dominated the communities (32.0%), followed by SAR11_clade (21.3%) and Flabovacteriales (9.4%). In cluster 3-2, SAR11_clade (17.9%), Flavobacteriales (17.3%) and Alteromonadales (10.8%) were dominant. Cluster 3-3 consisted of samples from November (191125) to February (200220) and was characterised by a dominance of SAR11_clade (26.1%) and Flavobacteriales (22.3%). In cluster 4-1, Flavobacteriales (14.6%), SAR11_clade (14.6%) and Parvibaculales (13.9%) were dominant, followed by Alteromonadales (8.2%) and Thiomicrospirales (7.0%). In cluster 4-2, Flavobacteriales (39.6%) and Rhodobacterales (27.4%) were highly dominant. Cluster 4-3 consisted of the sample at 100 m depth in May and Thiomicrospirales (19.3%) was dominant, following Flavobacteriales (25.9%). In cluster 4-4, Flavobacteriales (39.9%) was the most dominant taxa, followed by Rhodobacterales (14.3%), SAR11_clade (10.0%) and Parvibaculales (6.4%).
40303B1A-695D-5CBD-961C-0AE4E8B118B9
Spatiotemporal distributions of bacterial community clusters.
Taxonomic relative contributions of top 30 abundant Order (on average) to the bacterioplankton faction for each cluster.
https://binary.pensoft.net/fig/651872
SAR11 subclade Ia was dominant in SAR11 clade (53.1 ± 16.1%, Suppl. material 1: Fig. S3) and the contribution of this subclade to the bacterial community accounted for 7.6 ± 5.2% (up to 21.7%, Fig. 9). This subclade was ubiquitously found with relatively high contributions, except for the surface samples from spring to summer, when Flavobacteriales and Rhodobacterales were abundant (Suppl. material 1: Fig. S4). Subclade IV was the second dominant taxa in this clade (23.6 ± 18.6%, Suppl. material 1: Fig. S3) and showed relatively high contributions to the bacterial community in cluster 1 and 3-1, 3-2 and 3-3 (autumn to winter, 8.2 ± 4.6%, up to 18.5%, Fig. 9, Suppl. material 1: Fig. S4), while the relative contribution of subclade IV was small in the other clusters (1.3 ± 1.6%). Subclade II was the third dominant in this clade (15.5 ± 10.0%, Fig. 9). Except for clusters 4-2, 4-3 and 4-4 (0.89 ± 0.26%), subclade II showed relatively stable contributions to the bacterial community in the other clusters (1.8 ± 1.1%). The occurrence of subclade II was positively correlated to that of subclade Ia (n = 74, r = 0.79, p < 0.001, figure not shown).
855AF268-0139-52AB-898A-754DC7130A3C
Relative abundance of SAR11 subclades in bacterial community in each cluster.
https://binary.pensoft.net/fig/651873Cluster analysis of eukaryotic (phytoplankton) community
The eukaryotic community was classified into seven clusters, which represent clear seasonal and vertical patterns (Fig. 10, Suppl. material 1: Fig. S1). Cluster 7 was further classified into four clusters (7-1, 7-2, 7-3 and 7-4). Cluster 1, 3 and 4 composed of samples from 100 m depth of October, Spring (April–May) and February, respectively. Cluster 2 composed mostly from surface samples in late spring to autumn (May–September). Cluster 5 composed of samples from all vertical depths in November and the subsurface layer of different months. Cluster 6 composed of samples in all vertical depths of late spring to summer (May–July). Cluster 7 composed of samples in all vertical depths of winter to early spring (December–April).
BDD33D86-E4A6-52B6-8F71-D94EB9456158
Spatiotemporal distributions of eukaryotic community clusters.
https://binary.pensoft.net/fig/651874
The taxonomic groups of phytoplankton (diatom, dinoflagellates, haptophytes and green algae) were represented by heatmap at class level (Fig. 11). Noctilucales (dinoflagellates) showed a sudden increase in the surface in May and July. The abundance of reads in Peridinales (dinoflagellates) increased on the surface in May and deeper than 30 m in November, whereas Gymnodiniales (dinoflagellates) increased in the shallower depths (less than 10 m) in July and November to February. Syndiniales (dinoflagellates) showed large abundance in sequence reads in the deeper depths (deeper than 30 m with the largest at 100 m depth) of spring to autumn. Amongst diatoms, Coscinodiscophyceae showed large abundance in sequence reads from spring to autumn in the shallower depths near the surface (≥ 10 m) and especially large in May from the surface to 100 m depth. Mediophyceae (diatoms) showed a sudden increase on the surface in April. Amongst green algae, Mamiellophyceae showed large abundance in reads during winter in all depths and in April with peak in 30 m depth. Amongst haptophytes, Prymnesiales showed the largest abundance in reads in April with a peak in the surface (≥ 30 m) and during winter at all depths.
24F7183F-1589-5D13-B7FD-5BE56C8C5BFB
Heatmap of taxonomic groups that include phytoplankton community (diatom, dinoflagellates, haptophytes and green algae) at class levels identified by 18S rRNA in Sagami Bay. The colour scale of the heatmap indicates relative abundance in reads between the samples, colour red indicates larger abundance and blue indicates smaller abundance.
https://binary.pensoft.net/fig/651875
Additional cluster analysis was conducted for each taxonomic group that include phytoplankton community (diatom 223 OTUs, dinoflagellates 489 OTUs, green algae 77 OTUs and haptophytes 60 OTUs) and relative abundance top 15 OTUs classified to species level (Suppl. material 1: Figs S5–S8). Amongst diatoms, Leptocylindrusconvexus (16.3%) most dominated in relative reads abundance, followed by Minutocelluspolymorphusin (14.2%), Skeletonemamenzelii (5.6%), Thalassiosiraoceanica (4.2%), Brockmanniellabrockmannii (3.0%), Stephanopyxisturris (2.8%), Pseudo-nitzschia delicatissima (2.3%), Thalassiosiratenera (1.4%), Cylindrothecaclosterium (1.4%) and Thalassiosiraallenii (1.3%). Leptocylindrusconvexus was abundant in May from the surface to 100 m depths and June–July in the subsurface layer. Thalassiosiraallenii dominated from the surface to 30 m depths in June–August. Cylindrothecaclosterium and C.rostratus was abundant in July and especially C.rostratus dominated in the surface layer. Stephanopyxisturris dominated in the subsurface layer during winter (October–March) and most dominated in October at 100 m depth. Skeletonemamenzelii dominated in winter near the surface. Minutocelluspolymorphus, Thalassiosiraoceanica and Brockmanniellabrockmannii were abundant mainly in winter. Pseudo-nitzschia australis dominated locally on the surface in May and the subsurface layer in June.
Amongst dinoflagellates, Karlodiniumveneficum (12.6%) most dominated in relative reads abundance, followed by Heterocapsarotundata (8.0%), Noctilucascintillans (5.2%), Pentapharsodiniumtyrrhenicum (4.4%), Takayamacf.pulchellum (2.3%), Prorocentrummicans (1.8%) and Lepidodiniumviride (1.3%). Karlodiniumveneficum was relatively abundant throughout the survey period, especially dominant at 50 m depths in March–June. Heterocapsarotundata dominated at the surface in May–September and at the subsurface layer (30–100 m) in November. Noctilucascintillans was locally abundant at the surface in May and July and near the surface in November. Pentapharsodiniumtyrrhenicum was relatively abundant in the subsurface layer throughout the survey period, especially at 100 m depth and Parvodiniuminconspicuum at 100 m depths in February.
Sequence reads abundance of haptophytes had an overall peak in April at depths shallower than 30 m. Amongst haptophytes, Chrysochromulinascutellum most dominated in relative reads abundance (13.5%), followed by Prymnesiumpigrum (13.3%), Chrysochromulinasimplex (10.9%), Phaeocystisglobose (10.6%), Prymnesiumkappa (3.6%), Phaeocystiscordata (3.0%), Chrysochromulinacampanulifera (1.7%), Prymnesiumpalpebrale (1.5%), Haptolinafragaria (1.3%) and Chrysochromulinastrobilus (1.3%). Chrysochromulinascutellum occurred throughout the survey period. Prymnesiumpigrum was relatively abundant in winter at depths shallower than 50 m. Chrysochromulinasimplex also occurred throughout the survey period, but dominated in May–July and P.globosa at the 100 m depths.
Sequence reads abundance of green algae overall peaked in winter. Amongst green algae, Ostreococcustauri (%) most dominated in relative reads abundance (27.0%), followed by Micromonaspusilla (20.5%), Bathycoccusprasinos (14.6%), Mantoniellasquamata (5.1%), Pycnococcusprovasolii (3.3%), Pyramimonasdisomata (0.4%) and Pyramimonasaustralis (0.3%). Ostreococcustauri dominated in December–June. Micromonaspusilla and Bathycoccusprasinos dominated in October–July and October–March, respectively.
DiscussionEukaryotic community
Astrosphaerahexagonalis/ Arachnosphaeramyriacantha (radiolarian, Rhizaria) dominated in cluster 1 (October 2018 at 100 m depth) of eukaryotic community (75.3%). These species form a monophyletic group by molecular phylogenetic analysis using 18S rRNA gene (Yuasa et al. 2009). They were recorded in the warm surface water in the Pacific Ocean including the water area around the Kuroshio Current and surface sediment (Suzuki and Sugiyama 2001; Yuasa et al. 2009; Boltovskoy et al. 2010). The surface layer of Sagami Bay is strongly influenced by the oligotrophic Kuroshio (Kawabe and Yoneno 1987; Iwata and Matsuyama 1989; Hashihama et al. 2008; Baek et al. 2009). Furthermore, the same species was also most abundant in sequence reads amongst Rhizaria in the waters around the Kuroshio current south of Japan, with its distribution peak at 500–2000 m depths (Sogawa et al., submitted manuscript, 2022). These species were located in the direction of high nitrate + nitrite concentration, based on the DBRDA results, suggesting their coupled distribution to high nutrient deep-sea water.
The genus Amoebophrya is known to parasitise on a wide range of other dinoflagellates (Gunderson et al. 2002; Park et al. 2004). It also parasitises red tide-forming species (Taylor 1968) and is detrimental to harmful bloom algae (Miller et al. 2012). Parasitism was reported to occur during the bloom events in the surface or sub-surface layers (Nishitani et al. 2021). However, in our study, the parasitic syndinian (order Sindiniales) Amoebophrya sp. (dinoflagellates) dominated in cluster 4 (February at 100 m depth) of the eukaryotic community. No significant difference was observed between 100 m depths in October 2018 and February 2020 in the cluster analysis of only dinoflagellates. Amoebophrya sp. showed positive correlation with three dinoflagellates species; another parasitic, but non-toxic dinoflagellate Pentapharsodiniumtyrrhenicum that forms resting cysts (r = 0.52), Euduboscquella sp. (r = 0.55) and unarmoured dinoflagellates Takayamacf.pulchellum (synonym Gymnodiniumpulchellum) that causes a red tide bloom which may lead to fish kills (r = 0.48). Many species of dinoflagellates form resting cysts and have different life cycles that consist of a cyst and motile stage. The genus Euduboscquella parasitises Tintinnids and produces cysts (Dolan et al. 2013; Shang et al. 2019; Liu et al. 2020). The temporal dominance of Amoebophrya sp. could be considered caused by turbulence from bottom sediment from sudden upwelling events.
In the present study area, early summer is the rainy season and direct typhoon passes occasionally occur during the summer. High precipitation and abundant nutrients supplied by terrestrial run-off cause an increase in chlorophyll a concentration a few days later (Baek et al. 2009; Tsuchiya et al. 2015). Subsequently, abundance of harmful algal bloom (HAB) species of dinoflagellates (Triposfurca, previously described as Ceratiumfurca, Triposfusus, previously described as Ceratiumfusus and Noctilucascintillans) and diatoms (Rhizosoleniadelicatula, Hemiaulussinensis and Navicula spp.) increase and exhaust nitrogen and phosphorus supply, whereas silicate concentrations are not limited (Miyaguchi et al. 2006; Baek et al. 2008, 2009). It can be inferred that the same succession occurred in our study period, based on nutrient and salinity results. Dinoflagellate species are common bloom-forming species in coastal waters; however, our study showed increased abundances for N.scintillans, but not for T.furca and T.fusus. Furthermore, no recent blooms of the two Tripos species occurred at the monitoring site (Shimode’s periodical observation, data not shown). On the other hand, small mixotrophic dinoflagellates Heterocapsarotundata (3–10 μm) and Karlodiniumveneficum dominated. The former species, widely distributed in tropical and temperate coastal areas, is known to form large blooms in temperate estuaries (Millette et al. 2017). The latter species is ichthyotoxic and is distributed worldwide causing massive fish kills in coastal areas (Müller et al. 2019). The Kuroshio Current formed a large meander, which began in August 2017, continuing for a long period and was also observed during the survey periods (Japan Meteorological Agency, 20 May 2021). The surface layer of Sagami Bay is strongly influenced by the oligotrophic Kuroshio as it changes inflow to the Bay, based on the Kuroshio path (Kawabe and Yoneno 1987; Iwata and Matsuyama 1989; Hashihama et al. 2008; Baek et al. 2009). Therefore, the long-term large meandering might have caused changes in the plankton community and occurrence of bloom-forming species in Sagami Bay.
Bacterial community
The bacterial communities transitioned seasonally in the surface waters: from autumn–winter in 2018 (cluster 3-1) to winter–spring (cluster 4-4), rainy season before summer (cluster 4-2) and summer–autumn (cluster 2) in 2019 and to winter in 2019–2020 (cluster 3-3). The succession could be derived from the change of dominant taxa. Synechococcales and SAR11_clade were dominant in cluster 3-1, Flavobacteriales and Rhodobacterales became dominant in clusters 4-4 and 4-2. In cluster 2, Flabovacteriales, Rhodobacteriales, SAR11_clade, Synechococcales and cellvibrionales were evenly dominant. After that, SAR11_clade and Flavobacteriales became dominant. Flavobacteriales and Rhodobacterales are well known as copiotrophic taxa and respond quickly to the flux of fresh labile DOM derived from phytoplankton, subsequently outcompeting other community members as observed during phytoplankton blooms (Buchan et al. 2014; Luria et al. 2017). In this study, samples clustered in 2, 4-2 and 4-4 were located in the direction of high chlorophyll a concentration, based on the DBRDA results, suggesting Flavobacteriales and Rhodobacterales may have responded quickly to the DOM supply released from phytoplankton.
SAR11 was negatively correlated to Flavobacteriales and Rhodobacterales in the present study. The inverse correlation between those taxa was also observed in the coastal Southern Ocean (Williams et al. 2013). Flavobacteriia is essential for the breakdown of organic matter and they make labile organic compounds available to Proteobacteria, including SAR11 (Williams et al. 2013). SAR11 cells oxidise many labile, low-molecular-weight organic compounds and their small cells with periplasms and many substrate-biding proteins (SBPs) are likely to be superior competitors for low-molecular-weight DOM in oligotrophic ocean environments (Giovannoni 2017). Indeed, dissolved organic carbon (DOC) concentration increases during spring to summer and decreases during winter in the mixed layer of Sagami Bay (Sugai et al. 2016), suggesting that the DOC dynamics and the succession from Flavobacteriales and Rhodobacterales to SAR11 are consistent. Their different ecological roles could separate their occurrences spatially and seasonally.
SAR11 can be divided into several subclades, each of which has a different response to the environment (Giovannoni 2017; Bolanos et al. 2021). In this study, three SAR11 subclades (Ia, II and IV) were dominant. Subclade Ia was found with relatively regular contributions in all clusters, whereas subclade IV sporadically increased in clusters 1, 3-1, 3-2 and 3-3. Subclade Ia (Ia.1 in cold environment and Ia.3 in temperate and tropical environments) was characterised by surface ecotypes (Giovannoni 2017) and dominated in subpolar regions in the Atlantic Ocean (Bolanos et al. 2021) and poly- and euhaline sites in San Francisco Bay (Rasmussen et al. 2021), which are consistent with the present study. Since there was no distinct spatiotemporal variation in subclade Ia, DBRDA analysis did not show a relationship between environmental factors and subclade Ia occurrence. At BATS site, subclade IV was found at lower abundance in the surface layer, predominantly around 80 m depth (Vergin et al. 2013). Besides, subclade IV was varying only slightly and not dominant throughout the seasons in the North Atlantic (Bolanos et al. 2021). Subclade II was relatively low (0.89 ± 0.26%) in clusters 4-2, 4-3 and 4-4 (mainly in spring) compared to that in the other clusters (2.2 ± 0.4%). There are seasonal blooms in subclade IIa on the surface and IIb in the mesopelagic during spring (Carlson et al. 2009; Carlson et al. 2010; Vergin et al. 2013), which disagreed with the results of the present study. The function of SAR11 subclades can differ depending on the ocean area and environments, although there is a presumed relationship between ecotype classification and taxonomy in SAR11 subclades (Giovannoni 2017).
Cyanobacteria were relatively abundant in clusters 2 and 3-1, accounting for up to 61% and DBRDA ordination revealed that cyanobacteria (Synechococcus_CC9902) were highly correlated with water temperature. The sporadic increase of cyanobacteria from summer to autumn is consistent with previous studies in the same Sagami Bay area (Mitbavkar et al. 2009; Ara et al. 2019). Cyanobacterial abundance and growth rate are dependent principally on the temperature (Chen and Liu 2010; Chen et al. 2014; Wang et al. 2020). In Sagami Bay, Synechococcus biomass (216 ± 43 mg C m-2 on average annually) was much higher than that of Prochlorococcus (11.5 ± 1.7 mg C m-2 on average annually) (Mitbavkar et al. 2009). This is consistent with the present results that the ratio of Synechococcus read numbers to those of Cyanobacteria was 97 ± 7%. These two cyanobacterial genera show different ecological distribution: Synechococcus increases in nutrient-rich waters (generally coastal waters), whereas Prochlorococcus is most abundant in offshore nutrient-depleted areas (Partensky et al. 1999a). The global distribution prediction for Prochlorococcus and Synechococcus indicate that the abundance of Prochlorococcus and Synechococcus will increase overall by the end of the 21st century under the RCP4.5 scenario; however, the community structure will change significantly in the mid-latitudes (around 35°N or 35°S) and the abundance of Synechococcus may decrease (Flombaum et al. 2013). Flombaum et al. (2013) found that temperature is the primary control of the biogeography of both taxa and nutrient concentration had a limited influence on it. However, Prochlorococcus has been reported to be in a competitive disadvantage in elevated nutrient environments (Partensky et al. 1999b) and its future prediction in coastal waters is still debatable. Therefore, it is essential to understand the dynamics of both Prochlorococcus and Synechococcus in the coastal areas to predict the dynamics of primary production in the future since marine cyanobacteria could be responsible for ~ 25% of ocean net primary productivity (Flombaum et al. 2013).
SAR324_clade increased at deep waters (100 m depth) in summer and autumn, especially in October (5.4%), November (3.3%), December (6.1%) 2018 and July 2019 (5.8%). The results suggested that SAR324 gradually increased in their relative abundance several months after the MLD exceeded 100 m and vertical mixing occurs to the deeper layers. The seasonal variation agreed with a previous study reporting that SAR324 increased in the upper mesopelagic in autumn in the North Atlantic (Bolanos et al. 2021). SAR324 has been associated with a chemolithotrophic lifestyle and their distribution was well documented in oxygen-deficient waters (Zaikova et al. 2010), hydrothermal plumes (Sheik et al. 2015) and aphotic layers (160–3000 m, Wright et al. 1997; López-García et al. 2001; Bolanos et al. 2021). Draft genome analysis has revealed that SAR324 is also involved in new metabolic features, such as aromatic compound degradation and C1 metabolism, such as methanotrophy (Cao et al. 2016). This group might have contributed to carbon cycling in the stratified, deep layers in temperate coastal waters.
Overall plankton community and diversity dynamics
The present study revealed clear seasonal and vertical variations in the eukaryotic and bacterial communities consistent with changes in the environmental parameters. Eukaryotic and bacterial communities showed synchronous seasonal and vertical changes corresponding to the physical change; summer when the thermocline develops, winter when the mixed layer depth (MLD) is deep and transition periods that MLD gradually shallows (spring) and gradually deepens (autumn). In particular, DBRDA results suggested that the succession of eukaryotic and bacterial communities is closely related to the variations of the MLD, followed by nutrients for the eukaryotic community and temperature and chlorophyll a for the bacterial community.
The planktonic community was also vertically uniform in winter when vertical mixing formed a uniform water column and abiotic environment, although there was a slight time lag between the eukaryotic community (November–February) and the bacterial community (December–April). Productivity and biomass of bacterial communities are known to depend on eukaryote-derived resources (Williams et al. 2013; Buchan et al. 2014; Tsuchiya et al. 2015; Sugai et al. 2016; Luria et al. 2017; Tsuchiya et al. 2019); therefore, the bacterial community would also be formed dependent on the eukaryotic community. Significantly different clusters formed during different seasons at 100 m depths (clusters 1, 3 and 4 for eukaryotic community and clusters 1 and 4-3 for bacterial community). Upwelling of deep seawater supplies rich nutrients to the coastal areas from Sagami Trough (more than 1000 m deep), located in the centre of the Bay (Takahashi et al. 1980; Taira and Teramoto 1985; Kawabe and Yoneo 1987). In our study, deep-waters, rich in nutrients, were supplied to deeper layers (100 m) when vertical mixing weakened. It is interesting to note that the different eukaryotic communities formed during each upwelling event (cluster 1 in October 2018, cluster 3 in April–May 2019 and cluster 4 in February 2020). The annual dynamics of the plankton and bacterial communities are driven not only by the seasonal changes of abiotic and biotic factors and short-term rapid changes by river water inflow and phytoplankton blooms, but also largely influenced by deep-seawater upwelling due to the unique seafloor topography.
Vertical characteristics of community diversity were different between the eukaryotic and bacterial communities; the lowest diversity was observed in the surface water in both communities, whereas the highest diversity was observed in intermediate depth layers (10–50 m) for the eukaryotic community and in the deeper layer (100 m) for the bacterial community. The number of OTUs and the Evenness were both low at the surface and 100 m layer in the eukaryotic community. The higher OTU number in the intermediate depth layers could be explained by seasonally changing MLD; the subsurface chlorophyll maxima (SCM), formed at the optimal depth where light availability and upward nutrient flux are sufficient to sustain phytoplankton growth in the water column, is characterised by a high species diversity of phytoplankton (Furuya and Marumo 1983). It is also well known that ultraviolet radiation is harmful to many organisms and the surface is illuminated most with ultraviolet light (Kuwahara et al. 2015). The Evenness could be very low in the surface and 100 m layers due to phytoplankton and dinoflagellate blooms caused by short-term rapid environmental change by river water inflow and upwelling events as noted above. In the bacterial community, the number of OTUs does not vary much amongst depths, but the Evenness gradually increased with increasing depth and reached the highest in the 100 m layer along with the Diversity Index. This could be explained by low oxygen concentration in the deeper layers (except temporal vertical mixing in winter) and higher concentration near the surface as reported in the other water areas (Stevens et al. 2008; Signori et al. 2014).
Acknowledgements
This research was supported by the Japan Aerospace Exploration Agency (JAXA) Global Change Observation Mission – Climate (GCOM-C) project grant JX-PSPC-524336 and a Grant-in-Aid (Establishing a network of environment and fisheries information) by the Ministry of Agriculture, Forestry and Fisheries of Japan. We would like to thank all members of Team Manazuru in Yokohama National University and Soka University for observation support and members of Fisheries Research Institute, Japan Fisheries Research and Education Agency for measurement support.
ReferencesAbadDAlbainaAAguirreMLaza‐MartínezAUriarteIIriarteAVillateFEstonbaA (2016) Is metabarcoding suitable for estuarine plankton monitoring? A comparative study with microscopy.163: 1–13. https://doi.org/10.1007/s00227-016-2920-0AbadDAlbainaAAguirreMEstonbaA (2017) 18S V9 metabarcoding correctly depicts plankton estuarine community drivers.584: 31–43. https://doi.org/10.3354/meps12373AminSAHmeloLRVan TolHMDurhamBPCarlsonLTHealKRMoralesRLBerthiaumeCTParkerMSDjunaediBIngallsAEParsekMRMoranMAArmbrustEV (2015) Interaction and signalling between a cosmopolitan phytoplankton and associated bacteria.522: 98–101. https://doi.org/10.1038/nature14488AndersonMJGorleyRNClarkeKR (2008) PRIMER-E, Plymouth, 214 pp.AraKHiromiJNotesA (2006) Temporal variability in primary and copepod production in Sagami Bay, Japan. Journal of Plankton Research 29: i85–i96. https://doi.org/10.1093/plankt/fbl069AraKFukuyamaSOkutsuTNagasakaSShiomotoA (2019) Seasonal variability in phytoplankton carbon biomass and primary production, and their contribution to particulate carbon in the neritic area of Sagami Bay, Japan.14: 224–250. https://doi.org/10.3800/pbr.14.224AzamFFenchelTFieldJGGrayJJMeyer-ReilLAThingstadF (1983) The ecological role of water-column microbes in the sea.10: 257–263. https://doi.org/10.3354/meps010257AzamF (1998) Microbial control of oceanic carbon flux: The plot thickens.280: 694–696. https://doi.org/10.1126/science.280.5364.694AzamFMalfattiF (2007) Microbial structuring of marine ecosystems.5: 782–791. https://doi.org/10.1038/nrmicro1747BaekSHShimodeSHanMSKikuchiT (2008) Population development of the dinoflagellates Ceratiumfurca and Ceratiumfusus during spring and early summer in Iwa Harbor, Sagami Bay, Japan.43: 49–59. https://doi.org/10.1007/BF03022431BaekSHShimodeSKimHHanMSKikuchiT (2009) Strong bottom-up effects on phytoplankton community caused by a rainfall during spring and summer in Sagami Bay, Japan.75: 253–264. https://doi.org/10.1016/j.jmarsys.2008.10.005BakiMAMotegiCShibataAFukudaHShimodeSKikuchiT (2009) Temporal changes in chlorophyll a concentrations and bacterial, viral, and heterotrophic nanoflagellate abundances in the coastal zone of Sagami Bay, Japan: implications of top-down and bottom-up effects.32: 1–10.BochdanskyABHerndlGJ (1992) Ecology of amorphous aggregations (marine snow) in the Northern Adriatic Sea. 111. Zooplankton interactions with marine snow.87: 135–146. https://doi.org/10.3354/meps087135BokulichNAKaehlerBDRideoutJRDillonMBolyenEKnightRHuttleyGACaporasoJG (2018) Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 6: e90. https://doi.org/10.1186/s40168-018-0470-zBoltovskoyDKlingSATakahashiKBjørklundK (2010) World atlas of distribution of recent Polycystina (Radiolaria). Palaeontological Electronica 13.3.18A: 1–230.BolanosLMChoiCJWordenAZBaetgeNCarlsonCAGiovannoniS (2021) Seasonality of the microbial community composition in the North Atlantic.8: 1–16. https://doi.org/10.3389/fmars.2021.624164BrayJRCurtisJT (1957) An Ordination of the Upland Forest Communities of Southern Wisconsin.27: 325–349. https://doi.org/10.2307/1942268BuchanALeCleirGRGulvikCAGonzálezJM (2014) Master recyclers: features and functions of bacteria associated with phytoplankton blooms.12: 686–698. https://doi.org/10.1038/nrmicro3326.CamachoCCoulourisGAvagyanVMaNPapadopoulosJBealerKMaddenTL (2009) BLAST+: architecture and applications. BMC Bioinformatics 10: e421. https://doi.org/10.1186/1471-2105-10-421CaoHDongCBougouffaSLiJZhangWShaoZBajicVBQianPY (2016) Delta-proteobacterial SAR324 group in hydrothermal plumes on the South Mid-Atlantic Ridge. Scientific Report 6: e22842. https://doi.org/10.1038/srep22842CarlsonCAMorrisRParsonsRTreuschAHGiovannoniSJVerginK (2009) Seasonal dynamics of SAR11 populations in the euphotic and mesopelagic zones of the northwestern Sargasso Sea.3: 283–295. https://doi.org/10.1038/ismej.2008.117CarlsonCHansellDANelsonNSiegelDSmethieWKhatiwalaSMeyersMMHalewoodER (2010) Dissolved organic carbon export and subsequent remineralization in the mesopelagic and bathypelagic realms of the North Atlantic Basin.57: 1433–1445. https://doi.org/10.1016/j.dsr2.2010.02.013ChenBLiuH (2010) Relationships between phytoplankton growth and cell size in surface oceans: Interactive effects of temperature, nutrients, and grazing.55: 965–972. https://doi.org/10.4319/lo.2010.55.3.0965ChenBLiuHHuangBWangJ (2014) Temperature effects on the growth rate of marine picoplankton.505: 37–47. https://doi.org/10.3354/meps10773CheungKLYHuenJHouryWAOrtegaJ (2010) Comparison of the multiple oligomeric structures observed for the Rvb1 and Rvb2 proteins.8: 77–88. https://doi.org/10.1139/o09-159ClarkeKR (1993) Non‐parametric multivariate analyses of changes in community structure.18: 117–143. https://doi.org/10.1111/j.1442-9993.1993.tb00438.xClarkeKRGorleyRN (2015) PRIMER-E, Plymouth, 296 pp.DangHTiegangLMingnaCGuiqiaoH (2008) Cross-Ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters.74: 52–60. https://doi.org/10.1128/AEM.01400-07DolanJRMontagnesDJSAgathaSDCoatsDWStoeckerDK (2013) John Wiley & Sons, Ltd., 296 pp. https://doi.org/10.1002/9781118358092DzhembekovaNUrusizakiSMonchevaSIvanovaPNagaiS (2017) Applicability of massively parallel sequencing on monitoring harmful algae at Varna Bay in the Black Sea.68: 40–51. https://doi.org/10.1016/j.hal.2017.07.004EdgarRCHaasBJClementeJCQuinceCKnightR (2011) UCHIME improves sensitivity and speed of chimera detection.27: 2194–2200. https://doi.org/10.1093/bioinformatics/btr381FenchelT (1988) Marine Plankton Food Chains.19: 19–38. https://doi.org/10.1146/annurev.es.19.110188.000315FlombaumPFlombaumPGallegosJLGordilloRARincónJZabalaLLJiaoNKarlDMLiWKWLomasMWVenezianoDVeraCSVrugtJAMartinyAC (2013) Present and future global distributions of the marine CyanobacteriaProchlorococcus and Synechococcus.110: 9824–9829. https://doi.org/10.1073/pnas.1307701110FuruyaKMarumoR (1983) The structure of the phytoplankton community in the subsurface chlorophyll maxima in the western North Pacific Ocean.5: 393–406. https://doi.org/10.1093/plankt/5.3.393GiovannoniSJ (2017) SAR11 Bacteria: The Most Abundant Plankton in the Oceans.9: 231–255. https://doi.org/10.1146/annurev-marine-010814-015934GundersonJHJohnSABomanWCCoatsDW (2002) Multiple strains of the parasitic dinoflagellate Amoebophrya exist in Chesapeake Bay.49: 469–474. https://doi.org/10.1111/j.1550-7408.2002.tb00230.xHashihamaFHorimotoNKandaJFuruyaKIshimaruTSainoT (2008) Temporal variation in phytoplankton composition related to water mass properties in the central part of Sagami Bay.64: 23–37. https://doi.org/10.1007/s10872-008-0002-8HaysGCRichardsonAJRobinsoneC (2005) Climate change and marine plankton.20: 337–344. https://doi.org/10.1016/j.tree.2005.03.004HerndlJGPeduzziP (1988) The ecology of amorphous aggregations (marine snow) in the northern Adriatic Sea.9: 79–90. https://doi.org/10.1111/j.1439-0485.1988.tb00199.xHiraiJKatakuraSKasaiHNagaiS (2017) Cryptic zooplankton diversity revealed by a metagenetic approach to monitoring metazoan communities in the coastal waters of the Okhotsk Sea, northeastern Hokkaido. Frontiers in Marine Science 4: е379. https://doi.org/10.3389/fmars.2017.00379IbarbalzFMHenryNBrandãoMCMartiniSBusseniGByrneHCoelhoLPEndoHGasolJMGregoryACMahéFRigonatoJRoyo-LlonchMSalazarGSanz-SáezIScalcoESoviadanDZayedAAZingoneALabadieKFerlandJMarecCKandelsSPicheralMDimierCPoulainJPisarevSCarmichaelMPesantSBabinMBossEIudiconeDJaillonOAcinasSGOgataHPelletierEStemmannLSullivanMBSunagawaSBoppLDe VargasCKarp-BossLWinckerPLombardFBowlerCZingerLAcinasSGBabinMBorkPBossEBowlerCCochraneGDe VargasCFollowsMGorskyGGrimsleyNGuidiLHingampPIudiconeDJaillonOKandelsSKarp-BossLKarsentiENotFOgataHPesantSPoultonNRaesJSardetCSpeichSStemmannLSullivanMBSunagawaSWinckerP (2019) Global trends in marine plankton diversity across Kingdoms of life.179: 1084–1097. https://doi.org/10.1016/j.cell.2019.10.008IwataSMatsuyamaM (1989) Surface circulation in Sagami Bay: the response to variations of the Kuroshio axis.45: 310–320. https://doi.org/10.1007/BF02123485KawabeMYonenoM (1987) Water and flow variations in Sagami Bay under the influence of the Kuroshio path.43: 283–294. https://doi.org/10.1007/BF02108696Kita-TsukamotoKYaoKKamiyaAYoshizawaSUchiyamaNKogureKWadaM (2006) Rapid identification of marine bioluminescent bacteria by amplified 16S ribosomal RNA gene restriction analysis.256: 298–303. https://doi.org/10.1111/j.1574-6968.2006.00129.xKokSPKikuchiTTodaTKurosawaN (2012a) Diversity analysis of protistan microplankton in Sagami Bay by 18S rRNA gene clone analysis using newly designed PCR primers.68: 599–613. https://doi.org/10.1007/s10872-012-0121-0KokSPKikuchiTTodaTKurosawaN (2012b) Diversity and community dynamics of protistan microplankton in Sagami Bay revealed by 18S rRNA gene clone analysis.7: 75–86. https://doi.org/10.3800/pbr.7.75KuwaharaVSNozakiSNakanoJTodaTKikuchiTTaguchiS (2015) 18-year variability of ultraviolet radiation penetration in the mid-latitude coastal waters of the western boundary Pacific.160: 1–9. https://doi.org/10.1016/j.ecss.2015.03.029LiWGodzikA (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences.22: 1658–1659. https://doi.org/10.1093/bioinformatics/btl158LittleAEFRobinsonCJPetersonSBRaffaKFHandelsmanJ (2008) Rules of engagement: Interspecies interactions that regulate microbial communities. Annual Review of Microbiology 62: 375–401. https://doi.org/10.1146/annurev.micro.030608.101423LiuLWangZLuS (2020) Diversity and geographical distribution of resting stages of eukaryotic algae in the surface sediments from the southern Chinese coastline based on metabarcoding partial 18S rDNA sequences.41: 1–17. https://doi.org/10.1111/maec.12585López-GarcíaPLópez-LópezAMoreiraDRodríguez-ValeraF (2001) Diversity of free-living prokaryotes from a deep-sea site at the Antarctic Polar Front.36: 193–202. https://doi.org/10.1111/j.1574-6941.2001.tb00840.xLucasJWichelsATeelingHChafeeMScharfeMGerdtsG (2015) Annual dynamics of North Sea bacterioplankton: seasonal variability superimposes short-term variation. FEMS Microbiology Ecology 91: fiv099. https://doi.org/10.1093/femsec/fiv099LuriaCMAmaral-ZettlerLADucklowHWRepetaDJRhyneALRichJJ (2017) Seasonal shifts in bacterial community responses to phytoplankton-derived dissolved organic matter in the western Antarctic Peninsula. Frontiers in Microbiology 8: е2117. https://doi.org/10.3389/fmicb.2017.02117MillerJJDelwicheCFCoatsEW (2012) Ultrastructure of Amoebophrya sp. and its changes during the course of infection.163: 720–745. https://doi.org/10.1016/j.protis.2011.11.007MilletteNCPiersonJJAcevesAStoeckerDK (2017) Mixotrophy in Heterocapsarotundata: A mechanism for dominating the winter phytoplankton.62: 836–845. https://doi.org/10.1002/lno.10470Ministry of Environment of Japan (2011) Marine Biodiversity Conservation Strategy. 58 pp. https://www.env.go.jp/nature/biodic/kaiyo-hozen/pdf/pdf_pf_eng_all.pdfMiyaguchiHFujikiTKikuchiTKuwaharaVSTodaT (2006) Relationship between the bloom of Noctilucascintillans and environmental factors in the coastal waters of Sagami Bay, Japan.28: 313–324. https://doi.org/10.1093/plankt/fbi127MitbavkarSSainoTHorimotoNKandaJIshimaruT (2009) Role of environment and hydrography in determining the picoplankton community structure of Sagami Bay, Japan.65: 195–208. https://doi.org/10.1007/s10872-009-0019-7MontereyGLevitusS (1997) Seasonal variability of mixed layer depth for the world ocean. NOAA Atlas NESDIS 14, U.S. Government Printing Office, Washington D.C. 32(12): е87 https://doi.org/10.1029/2005GL022463MüllerMNDorantes-ArandaJJSegerABotanaMTBrandiniFPHallegraeffGM (2019) Ichthyotoxicity of the Dinoflagellate Karlodiniumveneficum in response to changes in seawater pH.26: 1–6. https://doi.org/10.3389/fmars.2019.00082NagaiSYamamotoKHataNItakuraS (2012) Study of DNA extraction methods for use in loop-mediated isothermal amplification detection of single resting cysts in the toxic dinoflagellates Alexandriumtamarense and A.catenella.7: 51–56. https://doi.org/10.1016/j.margen.2012.03.002NishitaniGYamamotoKNakajimaMShibataYSato-OkoshiWYamaguchiM (2021) A novel parasite strain of Amoebophrya sp. infecting the toxic dinoflagellate Alexandriumcatenella (Group I) and its effect on the host bloom in Osaka Bay, Japan. Harmful Algae 110: е102123. https://doi.org/10.1016/j.hal.2021.102123ParkMGYhiWCoatsDW (2004) Parasities and phytoplankton, with special emphasis on dinoflagellate infection.51: 145–155. https://doi.org/10.1111/j.1550-7408.2004.tb00539.xPartenskyFBlanchotJVaulotD (1999a) Differential distribution and ecology of Prochlorococcus and Synechococcus in oceanic waters: a review. Bulletin de l’Institut océanographique de Monaco Numero Special 19: 457–476.PartenskyFHessWRVaulotD (1999b) Prochlorococcus, a marine photosynthetic prokaryote of global significance.63: 106–127. https://doi.org/10.1128/MMBR.63.1.106-127.1999PielouEC (1966) The measurement of diversity in different types of biological collections.13: 131–144. https://doi.org/10.1016/0022-5193(66)90013-0PomeroyLR (1974) The Ocean’s Food Web, A Changing Paradigm.24: 499–504. https://doi.org/10.2307/1296885PomeroyLRWilliamsPJ leBAzamFHobbieJE (2007) The microbial loop.20: 28–33. https://doi.org/10.5670/oceanog.2007.45RasmussenANDamashekJEloe-FadroshEAFrancisCA (2021) In-depth spatiotemporal characterization of planktonic archaeal and bacterial communities in north and south San Francisco Bay.81: 601–616. https://doi.org/10.1007/s00248-020-01621-7RuppertKMKlineRJRahmanMS (2019) Past, present, and future perspectives of environmental DNA (eDNA) metabarcoding: A systematic review in methods, monitoring, and applications of global eDNA.17: 1–29. https://doi.org/10.1016/j.gecco.2019.e00547SchlossPDGeversDWestcottSL (2011) Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLоS ONE 6: e27310. https://doi.org/10.1371/journal.pone.0027310ShangLHuZDengYLiuYZhaiXChaiZLiuXZhanZDobbsFCTangYZ (2019) Metagenomic Sequencing identifies highly diverse assemblages of dinoflagellate cysts in sediments from ships’ ballast tanks. Microorganisms 7: e250. https://doi.org/10.3390/microorganisms7080250ShannonCEWeaverW (1949) Urbana, IL: University of Illinois Press, 131 pp.SheikCSAnantharamanKBreierJASylvanJBEdwardsKJDickGJ (2015) Spatially resolved sampling reveals dynamic microbial communities in rising hydrothermal plumes across a back-arc basin.9: 1434–1445. https://doi.org/10.1038/ismej.2014.228ShimodeSTodaTKikuchiT (2006) Spatio-temporal changes in diversity and community structure of planktonic copepods in Sagami Bay, Japan.148: 581–597. https://doi.org/10.1007/s00227-005-0093-3SignoriCNThomasFEnrich-PrastAPolleryRCGSievertSM (2014) Microbial diversity and community structure across environmental gradients in Bransfield Strait, Western Antarctic Peninsula.5: 1–12. https://doi.org/10.3389/fmicb.2014.00647SildeverSKawakamiYKannoNKasaiHShiomotoAKatakuraANagaiS (2019) Toxic HAB species from the Sea of Okhotsk detected by a metagenetic approach, seasonality and environmental drivers. Harmful Algae 87: e101631. https://doi.org/10.1016/j.hal.2019.101631SinYWetzelLRAndersonCI (1999) Spatial and temporal characteristic of nutrient and phytoplankton dynamics in the York River Estuary, Virginia. Analysis of long-term data.22: 260–275. https://doi.org/10.4319/lo.2006.51.3.1410SinclairLOsmanOABertilssonSEileA (2015) Microbial community composition and diversity via 16S rRNA gene amplicons: Evaluating the Illumina Platform. PLoS ONE 10(2): e0116955. https://doi.org/10.1371/journal.pone.0116955SogawaSSugisakiHSaitoHOkazakiYShimodeSKikuchiT (2013) Congruence between euphausiid community and water region in the northwestern Pacific: particularly in the Oyashio-Kuroshio Mixed Water Region.69: 71–85. https://doi.org/10.1007/s10872-012-0158-0SogawaSHidakaKKamimuraYTakahashiMSaitoHOkazakiYShimizuYSetouT (2018) Environmental characteristics of spawning and nursery grounds of Japanese sardine and mackerels in the Kuroshio and Kuroshio Extension area.28: 454–467. https://doi.org/10.1111/fog.12423Sogawa S, Nakamura Y, Nagai S, Nishi N, Hidaka K, Shimizu Y, Setou T (undated) Vertical variation and hidden diversity of Alveolata and Rhizaria in the western Pacific revealed by environmental DNA metabarcoding. submitted to Deep-Sea Research Part I.StevensHUlloaO (2008) Bacterial diversity in the oxygen minimum zone of the eastern tropical South Pacific.10: 1244–1259. https://doi.org/10.1111/j.1462-2920.2007.01539.xSugaiYTsuchiyaKKuwaharaVSShimodeSKomatsuKImaiATodaT (2016) Bacterial growth rate and the relative abundance of bacteria to heterotrophic nanoflagellates in the euphotic and disphotic layers in temperate coastal waters of Sagami Bay, Japan.72: 577–587. https://doi.org/10.1007/s10872-016-0352-6SugaiYTsuchiyaKShimodeSTodaT (2018) Seasonal variations in microbial abundance and transparent exopolymer particle concentration in the sea surface microlayer of temperate coastal waters.81: 201–211. https://doi.org/10.3354/ame01869SugaiYTsuchiyaKShimodeSTodaT (2020) Photochemical production and biological consumption of CO in the SML of temperate coastal waters and their implications for air-sea CO exchange. Journal of Geophysical Research: Oceans 125(4): e2019JC015505. https://doi.org/10.1029/2019JC015505SuterLPolanowskiAClarkeLKitchenerJDeagleB (2020) Capturing open ocean biodiversity: comparing environmental DNA metabarcoding to the continuous plankton recorder.30: 3140–3157. https://doi.org/10.1111/mec.15587SuzukiNSugiyamaK (2001) Regular axopodial activity of Diplosphaerahexagonalis Haeckel (spheroidal spumellarian, Radiolaria).5: 131–140. https://doi.org/10.2517/prpsj.5.131TairaKTeramotoT (1985) Bottom currents in Nankai Trough and Sagami Trough.41: 388–398. https://doi.org/10.1007/BF02109033TakahashiMKoikeIIshimaruTSainoTFuruyaKFujitaYHattoriAIchimuraS (1980) Upwelling in Sagami Bay and adjacent water around the Izu Islands, Japan.36: 209–216. https://doi.org/10.1007/BF02070334TakasukaAAokiIMitaniI (2003) Evidence of growth-selective predation on larval Japanese anchovy Enbraulisjaponicus in Sagami Bay.252: 223–238. https://doi.org/10.3354/MEPS252223TanabeASNagaiSHidaKYasuikeMFujiwaraANakamuraYTakanoYKatakuraS (2016) Comparative study of the validity of three regions of the 18S-rRNA gene for massively parallel sequencing-based monitoring of the planktonic eukaryote community.16: 402–414. https://doi.org/10.1111/1755-0998.12459TaylorFJR (1968) Parasitism of the toxin-producing dinoflagellate Gonyaulaxcatenella by the endoparasitic dmoflagellate Amoebophryaceratii.25: 2241–2245. https://doi.org/10.1139/f68-197TsuchiyaKYoshikiTNakajimaRMiyaguchiHKuwaharaVSTaguchiSKikuchiTTodaT (2013) Typhoon-driven variations in primary production and phytoplankton assemblages in Sagami Bay, Japan: A case study of typhoon Mawar (T0511).8: 74–87. https://doi.org/10.3800/pbr.8.74TsuchiyaKKuwaharaVSHamasakiKTadaYIchikawaTYoshikiTNakajimaRImaiAShimodeSTodaT (2015) Typhoon-induced response of phytoplankton and bacteria in temperate coastal waters. Estuarine, Coastal and Shelf Science 167, Part B, 20: 458–465. https://doi.org/10.1016/j.ecss.2015.10.026TsuchiyaKSanoTTomiokaNKohzuAKomatsuKShinoharaRTakamuraNNakagawaMSugaiYKuwaharaVSTodaTFukudaHImaA (2019) Seasonal variability and regulation of bacterial production in a shallow eutrophic lake.64: 2441–2454. https://doi.org/10.1002/lno.11196VajraveluMMartinYAyyappanSMayakrishnanM (2018) Seasonal influence of physico-chemical parameters onphytoplankton diversity, community structure andabundance at Parangipettai coastal waters, Bay of Bengal, South East Coast of India.60: 114–127. https://doi.org/10.1016/j.oceano.2017.08.003ValentiniATaberletPMiaudCCivadeRHerderJThomsenPFBellemainEBesnardACoissacEBoyerFGaboriaudCJeanPPouletNRosetNCoppGHGeniezPPontDArgilleierCBaudoinJPerouxTCrivelliAJOlivierAAcquebergeMBrunMLMøllerPRWillerslevEDejeanT (2016) Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding.25: 929–942. https://doi.org/10.1111/mec.13428VerginKLVerginKLBeszteriBMonierAThrashJCTempertonBTreuschAHKilpertFWordenAZGiovannoniSJ (2013) High-resolution SAR11 ecotype dynamics at the Bermuda Atlantic Time-series Study site by phylogenetic placement of pyrosequences.7: 1322–1332. https://doi.org/10.1038/ismej.2013.32WangHZhangCChenFKanJ (2020) Spatial and temporal variations of bacterioplankton in the Chesapeake Bay: A re‐examination with high‐throughput sequencing analysis.65: 3032–3045. https://doi.org/10.1002/lno.11572WilliamsTJWilkinsDLongEEvansFDeMaereMZRafteryMJCavicchioliR (2013) The role of planktonic Flavobacteria in processing algal organic matter in coastal East Antarctica revealed using metagenomics and metaproteomics.15: 1302–1317. https://doi.org/10.1111/1462-2920.12017WrightTDVerginKLBoydPWGiovannoniSJ (1997) A novel delta-subdivision proteobacterial lineage from the lower ocean surface layer.63: 1441–1448. https://doi.org/10.1128/aem.63.4.1441-1448.1997YilmazPParfreyLWYarzaPGerkenJPruesseEQuastCSchweerTPepliesJLudwigWGlöcknerFO (2014) The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Research 42: D643–D648. https://doi.org/10.1093/nar/gkt1209YuasaTTakahashiOMayamaS (2009) A simple technique for extracting small subunit ribosomal DNAs from both the host and symbiont of a single radiolarian specimen.14: 1–4.ZaikovaEWalshDAStilwellCPMohnWWTortellPDHallamSJ (2010) Microbial community dynamics in a seasonally anoxic fjord: Saanich Inlet, British Columbia.12: 172–191. https://doi.org/10.1111/j.1462-2920.2009.02058.xSupplementary materials10.3897/mbmg.6.78181.suppl1A16C3C83-D384-5B64-8C70-462FE6A7A26C
Figures S1–S8
images
Explanation note: Figure S1. Cluster analysis results (group average) of eukaryotic (upper) and bacterial community (lower) in Sagami Bay. (A) Red lines represented results of SIMPROF test. (B) Each point represents each sample and different color represent each cluster. Figure S2. Taxonomic relative contributions to the bacterioplankton faction for each sample in Order level. Figure S3. Relative reads abundance of SAR11 subclades in SAR11 clade in each cluster. Figure S4. Spatiotemporal distributions of SAR11 subclades Ia, II and IV shown by their relative reads abundances in total bacterial community. Figure S5. Cluster analysis of diatom and relative reads abundance of top 15 OTUs classified to species level. Figure S6. Cluster analysis of dinoflagellates and relative reads abundance of top 15 OTUs classified to species level. Figure S7. Cluster analysis of green algae and relative reads abundance of top 15 OTUs classified to species level. Figure S8. Cluster analysis of haptophytes and relative reads abundance of top 15 OTUs classified to species level.
https://binary.pensoft.net/file/651876This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Author: Sayaka Sogawa, Kenji Tsuchiya, Satoshi Nagai, Shinji Shimode, Victor S. Kuwahara