# Diversity and biogeography of planktonic diatoms in Svalbard fjords: The role of dispersal and Arctic endemism in phytoplankton community structuring

Luka Šupraha, Kerstin Klemm, Sandra Gran-Stadniczeñko, Cora Hörstmann, Daniel Vaulot, Bente Edvardsen, Uwe John

► To cite this version:

Luka Šupraha, Kerstin Klemm, Sandra Gran-Stadniczeñko, Cora Hörstmann, Daniel Vaulot, et al. Diversity and biogeography of planktonic diatoms in Svalbard fjords: The role of dispersal and Arctic endemism in phytoplankton community structuring. Elementa: Science of the Anthropocene, 2022, 10 (1), 10.1525/elementa.2021.00117. insu-03784351

HAL Id: insu-03784351

https://insu.hal.science/insu-03784351

Submitted on 23 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Diversity and biogeography of planktonic diatoms in Svalbard fjords: The role of dispersal and Arctic endemism in phytoplankton community structuring

Luka Šupraha<fnref n="1" />*, Kerstin Klemm<fnref n="2" />, Sandra Gran-Stadniczeňko<fnref n="1" />, Cora Hörstmann<fnref n="2" />, Daniel Vaulot<fnref n="3" />, Bente Edvardsen<fnref n="1" />, and Uwe John<fnref n="2" />,<fnref n="4" />,*

Understanding the processes that shape the community structure of Arctic phytoplankton is crucial for predicting responses of Arctic ecosystems to the ongoing ocean warming. In particular, little is known about the importance of phytoplankton dispersal by the North Atlantic Current and the prevalence and maintenance of Arctic endemism. We investigated the diversity and biogeography of diatoms from five Svalbard fjords and the Hausgarten observatory (Fram Strait) by combining diatom cultivation and 18S rRNA gene metabarcoding. In total, 50 diatom strains were isolated from the area during the HE492 cruise in August 2017. The strains were identified taxonomically using molecular and morphological approaches, and their biogeographic distribution was mapped using the local metabarcoding dataset and a global compilation of published metabarcoding datasets. Biogeographic analysis was also conducted for the locally most abundant diatom metabarcoding amplicon sequence variants. The biogeographic analyses demonstrated that Arctic diatoms exhibit three general biogeographic distribution types: Arctic, Arctic-temperate, and cosmopolitan. At Hausgarten and in outer Isfjorden on the west coast of Svalbard, the communities were dominated by genotypes with Arctic-temperate and cosmopolitan distribution. Diatom communities in nearby Van Mijenfjorden, inner Isfjorden and Kongsfjorden were dominated by genotypes with Arctic-temperate distribution, and cosmopolitan species were less abundant. The genotypes endemic to the Arctic had lower abundance on the west coast of Svalbard. The two northernmost fjords (Woodfjorden and Wijdefjorden) had a higher abundance of genotypes endemic to the Arctic. Our results demonstrate that the diatom communities in the Svalbard area consist of genotypes endemic to the Arctic, and genotypes with broader biogeographic distribution, all of which are further structured by local environmental gradients. Finer biogeographic patterns observed within Arctic-temperate and cosmopolitan genotypes suggest that certain genotypes can be used as indicators of increasing influence of Atlantic waters on the phytoplankton community structure in the Svalbard area.

**Keywords:** Arctic diatoms, Protist diversity, Metabarcoding, Biogeography, Atlantification, Endemism

## Introduction

Diatoms are important primary producers in Arctic marine ecosystems. They dominate both pelagic and sympagic (ice-associated) primary production (Eamer et al., 2013) and have a crucial role in Arctic food web dynamics (Lovejoy et al., 2002) as well as the biogeochemical cycling of

carbon (Le Moigne et al., 2015) and silica (Krause et al., 2018). With more than 1,000 described morphospecies (Poulin et al., 2011), Arctic diatoms exhibit a remarkable taxonomic diversity, along with a high degree of morphological, physiological and functional diversity (Fragoso et al., 2018). Accordingly, the abundance, community composition and distribution of diatoms in the Arctic realm ultimately shape ecosystem functioning (Wiedmann et al., 2020).

The Arctic today is a hotspot of climate change (Overland et al., 2019), with ocean warming, decreasing sea-ice cover (Kwok, 2018), and ocean freshening (Haine et al., 2015) likely to continue under projected carbon emission scenarios (Adakudlu et al., 2019). In this context, understanding the processes that are shaping Arctic diatom communities is essential for predicting the response of Arctic ecosystems to global warming. The structuring of Arctic diatom communities by seasonal and latitudinal

<footnote n="1">¹Section for Aquatic Biology and Toxicology, Department of Biosciences, University of Oslo, Oslo, Norway</footnote>

<footnote n="2">²Ecological Chemistry, Department of Biosciences, Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany</footnote>

<footnote n="3">³Sorbonne Université, CNRS, UMR 7144, Station Biologique de Roscoff, Roscoff, France</footnote>

<footnote n="4">⁴Helmholtz Institute for Functional Marine Biodiversity, Oldenburg, Germany</footnote>

<footnote>* Corresponding author:
Emails: luka.supraha@ibv.uio.no; uwe.john@awi.de</footnote>
patterns in light, nutrient availability, sea-ice and oceanographic conditions is generally well-understood (Eamer et al., 2013; Weydmann et al., 2013; Leu et al., 2015; Leeuwe et al., 2018). However, the fine-scale, genotype-level biogeographic patterns underlying the distribution of Arctic diatoms, particularly the relative importance of endemic Arctic species and of species that are dispersed with Atlantic circulation, are still unknown (Polyakov et al., 2017; Oziel et al., 2020).

Endemism is seemingly rare among marine protists, especially among planktonic species with large population sizes and the potential for global dispersal (Finlay, 2002). This paradigm has been challenged by evidence for endemism among marine planktonic protists, especially in polar areas (Darling et al., 2007; Lovejoy et al., 2007; Vyverman et al., 2010; Ribeiro et al., 2020; Yau et al., 2020). A certain degree of biogeographically restricted distribution has been reported for diatoms in the Arctic, notably for species associated with sea ice (Quillfeldt, 2000; Poulin et al., 2011). A recent, morphology-based analysis of the biogeography of selected cultured Arctic diatom species (Balzano et al., 2017) suggested that Arctic diatoms exhibit diverse biogeographic distribution patterns, which often extend beyond the Arctic realm into temperate and even warm latitudes.

In recent years, several studies have reported an apparent northward expansion of temperate phytoplankton, most notably bloom-forming haptophytes such as *Emiliania huxleyi* and *Phaeocystis* spp. (Orkney et al., 2020; Oziel et al., 2020), into the Atlantic sector of the Arctic. The phenomenon, termed "Atlantification", has been correlated with the increasing inflow of warm Atlantic waters and warming of Arctic waters and is seen as an indicator of global warming (Polyakov et al., 2017; Neukermans et al., 2018). The relevance of dispersed temperate or low-latitude diatom species for community structuring in the Arctic is unknown. Based on available morphology-based data, cosmopolitan species constitute a significant portion of the Arctic diatom flora (Quillfeldt, 2000; Poulin et al., 2011). Unclear is whether they represent warm-water species dispersed with Atlantic circulation, truly cosmopolitan species also adapted to the Arctic environment, or cold-adapted cryptic species and genotypes restricted to the Arctic. Identifying these fine-scale biogeographic patterns is crucial for understanding processes that are shaping Arctic phytoplankton communities, and for developing tools for tracking changes in Arctic microalgal biodiversity related to the warming climate.

Exploring the biogeographic distribution of Arctic diatoms is challenging due to several methodological constraints. The morphology-based datasets are often limited by insufficient and inconsistent taxonomic identification. These issues are especially relevant for species complexes with cosmopolitan biogeographic distribution, which often include cryptic, cold-adapted species (Chamnansinp et al., 2013; Percopo et al., 2016) and cold-adapted genotypes (Stock et al., 2019) that are usually indistinguishable and easily overlooked in morphology-based surveys. On the other hand, combining morphological and molecular approaches (Balzano et al., 2017;

Ribeiro et al., 2020) provides a more robust taxonomic framework for identifying Arctic diatoms and accounting for their cryptic diversity. When combined with ever-increasing environmental metabarcoding datasets, this approach allows for an in-depth analysis of the global distribution of morphologically-verified diatom genotypes using species-level molecular markers such as 18S rRNA gene sequences (De Luca et al., 2019).

This study aimed to reveal the diversity and biogeographic distribution of Arctic planktonic diatoms and understand the processes that shape the diatom community structure in the Svalbard area, particularly the interplay of dispersal by Atlantification, Arctic endemism and local environmental forcing at eastern Fram Strait and along five Svalbard fjords. The work followed three main lines: 1) morphological and molecular characterisation of cultured diatoms strains isolated from the Svalbard area; 2) exploration of the global and local distribution patterns of the cultured diatom genotypes and the most abundant metabarcoding ASVs using a global compilation of 18S V4 rRNA gene metabarcoding datasets; and 3) assessment of the contribution of genotypes with distinct biogeographic distribution to late-summer diatom communities on the west coast of Svalbard.

## Material and methods

### *Studied area*

The samples were collected onboard R/V Heincke between July 29 and August 18, 2017, in the frame of the cruise number HE492 led by Alfred Wegener Institute (Germany). In total, 28 stations were sampled, forming transsects of 4 to 6 stations within five Svalbard fjords (south to north: Van Mijenfjorden, Isfjorden, Kongsfjorden, Woodfjorden and Wijdefjorden) and a 3-station transect at the long-term ecological research observatory Hausgarten, located in eastern Fram Strait (Figure 1A).

The studied area is a unique natural laboratory for exploring how biogeographic processes and environmental gradients are shaping Arctic diatom communities (Figure 1B and C). The Hausgarten stations are located within the West Spitzbergen Current, an extension of the North Atlantic current that brings warmer and nutrient-rich Atlantic waters into the Arctic (Soltwedel et al., 2016). Isfjorden and Kongsfjorden, on the west coast of Svalbard, are increasingly influenced by inflowing Atlantic waters (Nilsen et al., 2016; Menze et al., 2020), as evidenced by the long-term trends in hydrography (Pavlov et al., 2013; Tverberg et al., 2019) and primary production (Hegseth and Tverberg, 2013). In contrast, the characteristic topography of the southernmost Van Mijenfjorden, with an island blocking the entrance to the fjord, reduces the inflow of coastal waters, and the fjord is consequently less influenced by the Atlantic water inflow (Skarðhamar and Svendsen, 2010; Nilsen et al., 2016). Finally, the two northernmost fjords, Wijdefjorden and Woodfjorden, are located at the northern edge of the West Spitzbergen Current and are the least affected by the North Atlantic Current of all the fjords studied here (Menze et al., 2020). Along with different levels of exposure to advection of warmer water masses, the studied fjords exhibit distinct environmental
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Sampling stations and hydrographic conditions at the time of sampling. (A) Sampling stations and transects along the five Svalbard fjords and at Hausgarten area; (B) sea surface temperature and (C) salinity based on ferrybox and CTD data along with the satellite data derived from Copernicus ARCTIC_MULTIYEAR PHY_002_003 (last access: 2022-06-17), averaged over 20170803-20170816 with the Europe coastline derived from the European Environment Agency (https://www.eea.europa.eu/ds_resolveuid/06227e40310045408ac8be0d469e1189) and the map shown in the NSIDC Sea Ice Polar Stereographic North projection using qGIS v3.24.2.</figcaption></figure>

gradients defined by their geography and freshwater and nutrient input by glaciers and runoff (Cottier et al., 2010).

### CTD profiling and water sampling

A conductivity-temperature-depth (CTD) probe (Seabird sbe911+, Seabird, USA) was deployed at each station to determine the hydrographical profile. The CTD was equipped with additional sensors measuring turbidity, fluorescence and oxygen levels. All CTD data collected during the cruise have been deposited in the Pangea data repository (https://doi.pangaea.de/10.1594/PANGAEA.881296). The CTD probe was equipped with a rosette with

10-L Niskin bottles, which were deployed to collect water samples from discrete depths after the initial onboard examination of the profile.

### Culture isolation

Clonal cultures of planktonic diatoms were isolated from stations 6 and 7 located in Wijdefjorden and stations 20, 21 and 22 located at Hausgarten. A concentrated phytoplankton sample was collected from 0–30 m by a vertical net haul using a plankton net with mesh size of 20 µm. Upon arrival to deck the net haul sample was diluted to 1 L with filtered seawater, and 40 mL of this volume was
transferred to a Nunclon culturing flask with added 10 mL of algal growth medium (IMR ½ medium with salinity 34; Eppley et al., 1967). In addition, water samples for metabarcoding collected by the Niskin bottles on the CTD from 3, 15, and 30 m were pooled. Some of this material was pre-filtered through a 200-µm mesh sieve and cells were collected on a 20-µm mesh sieve, rinsed off the sieve with seawater and transferred to a culturing flask to which was added 1/10 parts of IMR ½ medium with salinity 34. Both types of samples were stored onboard at 4°C in an illuminated incubator until the isolation of cultures onboard and later in the lab at UiO.

The larger diatom species were obtained using capillary isolation. Concentrated samples were placed in a Petri dish and observed under a Zeiss Axiovert A1 inverted microscope (Zeiss, Oberkochen, Germany). Individual cells were isolated using a custom-made glass capillary connected to a rubber tube with a mouthpiece. Each cell was washed by repeated transfers to drops of filter-sterilized seawater and was finally transferred to an individual well on a 96-well culture plate (Thermo Fisher Scientific, USA) filled with approximately 300 µL of IMR ½ medium with salinity 34. Cells were kept in a culture room at 4°C, under medium light intensity (approximately 10–20 µmol m⁻² s⁻¹) and 16 h/8 h light/dark cycle, and were inspected regularly for growth and contamination. When the presence of a single, exponentially growing culture was confirmed, the culture was transferred to a larger, 12-well culture plate (Thermo Fisher Scientific, USA) in 4 mL of the growth medium. For the smaller species (< 20 µm) dilution cultures were prepared from total integrated water samples collected by Niskin bottles (from depths of 3, 15 and 30 m), with a 6-step dilution series along the rows of the 96-well culture plate (Thermo Fisher Scientific, USA) with a final dilution of x10⁶. Well plates were incubated in the same culture room as the capillary isolates, and the monocultures detected in the dilution series followed the same treatment as the capillary isolates.

When a culture reached high concentration in 4-mL wells, a drop (approximately 200 µL) of culture was inoculated into a glass cultivation tube with 10 mL medium and assigned a strain code. Subsamples of cultures were collected for live observation by light microscopy and imaging using a Zeiss Axioscope light microscope (Zeiss, Oberkochen, Germany) equipped with a Leica MC170 HD Camera (Leica, Germany). In addition, 1 mL of dense culture material was fixed with EM-grade glutaraldehyde (final concentration 1%; Electron Microscopy Sciences, Hatfield, PA, USA) and stored at 4°C for future morphological analysis. Cultures were maintained in 10-mL glass tubes and transferred to a fresh medium every 5–7 weeks. Subsequently, cultures were deposited to the Norwegian Culture Collection of Algae (NORCCA), where they were assigned a new strain code and made publicly available (Table 1; www.norcca.scrol.net).

## Molecular characterisation of cultures

The main aim of the sequencing effort was to obtain 18S and 28S rRNA gene sequences for all cultured genotypes,

with the particular focus on sequencing the V4 region of the 18S rRNA gene, which was used for mapping the biogeographic distribution of cultured genotypes. A 1.5 mL volume of a dense culture was collected in an Eppendorf tube and centrifuged for 15 min at 6000 g (Eppendorf centrifuge 5424 R) to form a pellet. The supernatant was discarded, and the pellet was stored at –20°C until DNA extraction, which was performed using Qiagen DNeasy blood and tissue kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol.

For amplification of the 18S rRNA gene, the initial PCR reaction was run using 1F as a forward primer and 1528R (Medlin et al., 1988) as reverse primers (Table S1). The PCR reaction steps were as follows: initial denaturation (94°C, 3.5 min), followed by 35 cycles of denaturation (94°C, 50 s), annealing (55°C, 50 s) and elongation (72°C, 2 min), and the final elongation step (72°C, 10 min). A variety of sequencing primers for the 18S rRNA gene was used (1F, 300R, 528F, 850F, 1055R, 1147R, and 1528R), depending on the strain (Table S1). For amplifying the partial 28S rRNA gene sequence, the D1R-F and D2C-R primer pair (Scholin et al., 1994; Table S1) was used in a PCR reaction with the following steps: initial denaturation at (94°C, 2 min), then 35 cycles of denaturation (94°C, 45 s), annealing (55°C, 45 s) and elongation (72°C, 1 min), followed by the final elongation step (72°C, 10 min). The 28S rRNA gene fragment was sequenced using the same primers as in the PCR reaction. Sequencing data (sequences and chromatograms) were quality-checked, and contigs were assembled into consensus sequences using the De Novo Assemble tool, trimmed, edited and aligned in the Geneious software (v. 11.1.5). Sequences were submitted to NCBI, with accession numbers OK147649-147680 for 18S rRNA gene sequences, and OK147681-147729 for 28S rRNA gene sequences (Table 1).

## Metabarcoding of field samples

A total of 60 L or 20 L of seawater were collected at three different depths with Niskin bottles (depths of 3 m, 15 m, and 30 m) and subsequently pooled. The 60 L were gravity-filtered over 200-µm and finally 20-µm mesh size sieves. For the smaller size fraction, after pre-filtration over 20 µm, an in-line filtration followed through 3-µm and 0.2-µm polycarbonate filters (Millipore, USA) with a Millipore Tripod unit (147-mm diameter; Millipore, USA) and a peristaltic pump during a maximum of 30 min. These filtrations yielded three size fractions: microplankton (20–200 µm), nanoplankton (3–20 µm), and picoplankton (0.2–3 µm). Filters were then added to warm lysis buffer, fast-frozen in liquid nitrogen, and kept at –80°C until further processing.

DNA extraction was performed with the NucleoSpin® Soil kit (Macherey-Nagel, Düren, Germany) in SL1 buffer according to the manufacturer's protocol. The 18S rRNA gene V4 region was targeted for amplification using forward and reverse primers (Bradley et al., 2016; Table S2) with overhang adapters attached. The Illumina overhang nucleotide sequences were added to the 18S amplicon PCR forward and reverse primer. The workflow for preparation of 18S rRNA gene amplicons for the Illumina MiSeq
<table><caption>Table 1. List of diatom strains isolated during this study</caption><thead><tr><th rowspan="2">Taxon</th><th rowspan="2">Strain Code (NORCCA)<sup>a</sup></th><th rowspan="2">Station<sup>b</sup></th><th colspan="2">GenBank Accession Number<sup>c</sup></th></tr><tr><th>18S</th><th>28S</th></tr></thead><tbody><tr><td colspan="5"><strong>Bacillariaceae</strong></td></tr><tr><td>Cylindrotheca closterium</td><td>HE492-63 (UiO 469)</td><td>7 (10)</td><td>OK147662</td><td>OK147703</td></tr><tr><td>Fragilariopsis sp.<sup>d</sup></td><td>HE492-50</td><td>21 (28)</td><td>OK147664</td><td>OK147701</td></tr><tr><td>Pseudo-nitzschia granii</td><td>HE492-40 (UiO 453)</td><td>20 (27)</td><td>OK147667</td><td>OK147699</td></tr><tr><td>Pseudo-nitzschia turgidula</td><td>HE492-37 (UiO 451)</td><td>20 (27)</td><td>OK147665</td><td>OK147698</td></tr><tr><td>Pseudo-nitzschia turgidula</td><td>HE492-47 (UiO 460)</td><td>21 (28)</td><td>OK147666</td><td>OK147697</td></tr><tr><td colspan="5"><strong>Attheyaceae</strong></td></tr><tr><td>Attheya longicornis</td><td>HE492-59 (UiO 467)</td><td>20 (27)</td><td>OK147651</td><td>OK147726</td></tr><tr><td colspan="5"><strong>Cymatosiraceae</strong></td></tr><tr><td>Arcocellulus cornucervis</td><td>HE492-38</td><td>20 (27)</td><td>OK147650</td><td>OK147729</td></tr><tr><td>Arcocellulus cornucervis</td><td>HE492-43 (UiO 456)</td><td>21 (28)</td><td>OK147649</td><td>OK147728</td></tr><tr><td>Arcocellulus cornucervis</td><td>HE492-71 (UiO 471)</td><td>6 (9)</td><td>na<sup>e</sup></td><td>OK147727</td></tr><tr><td>Undetermined Cymatosiraceae<sup>d</sup></td><td>HE492-46 (UiO 459)</td><td>21 (28)</td><td>OK147663</td><td>OK147702</td></tr><tr><td colspan="5"><strong>Chaetocerotaceae</strong></td></tr><tr><td>Chaetoceros contortus</td><td>HE492-06 (UiO 434)</td><td>20 (27)</td><td>OK147652</td><td>OK147725</td></tr><tr><td>Chaetoceros convolutus</td><td>HE492-64</td><td>7 (10)</td><td>OK147653</td><td>na</td></tr><tr><td>Chaetoceros debilis (Clade 1)</td><td>HE492-61</td><td>20 (27)</td><td>na</td><td>OK147722</td></tr><tr><td>Chaetoceros debilis (Clade 2)</td><td>HE492-19</td><td>20 (27)</td><td>OK147654</td><td>OK147724</td></tr><tr><td>Chaetoceros debilis (Clade 2)</td><td>HE492-28</td><td>22 (29)</td><td>na</td><td>OK147723</td></tr><tr><td>Chaetoceros diadema (Clade 1)</td><td>HE492-13 (UiO 440)</td><td>20 (27)</td><td>na</td><td>OK147720</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-12 (UiO 439)</td><td>20 (27)</td><td>OK147655</td><td>OK147721</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-15 (UiO 441)</td><td>20 (27)</td><td>na</td><td>OK147719</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-16 (UiO 442)</td><td>20 (27)</td><td>na</td><td>OK147718</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-17 (UiO 443)</td><td>20 (27)</td><td>na</td><td>OK147717</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-18 (UiO 444)</td><td>20 (27)</td><td>OK147656</td><td>OK147716</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-21</td><td>7 (10)</td><td>OK147657</td><td>OK147715</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-58 (UiO 446)</td><td>20 (27)</td><td>OK147658</td><td>OK147714</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-65</td><td>22 (29)</td><td>na</td><td>OK147713</td></tr><tr><td>Chaetoceros diadema (Clade 2)<sup>d</sup></td><td>HE492-75 (UiO 475)</td><td>20 (27)</td><td>na</td><td>OK147712</td></tr><tr><td>Chaetoceros neogracilis (Clade 1)</td><td>HE492-73 (UiO 473)</td><td>6 (9)</td><td>OK147659</td><td>OK147711</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-36 (UiO 450)</td><td>6 (9)</td><td>OK147661</td><td>OK147705</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-02</td><td>20 (27)</td><td>na</td><td>OK147710</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-03</td><td>20 (27)</td><td>OK147660</td><td>OK147709</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-05 (UiO 433)</td><td>20 (27)</td><td>na</td><td>OK147708</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-11</td><td>20 (27)</td><td>na</td><td>OK147707</td></tr><tr><td>Chaetoceros sp.<sup>d</sup></td><td>HE492-27</td><td>22 (29)</td><td>na</td><td>OK147706</td></tr><tr><td colspan="5"><strong>Leptocylindraceae</strong></td></tr><tr><td>Leptocylindrus sp.</td><td>HE492-14</td><td>20 (27)</td><td>na</td><td>OK147700</td></tr><tr><td colspan="5"><strong>Skeletonemataceae</strong></td></tr><tr><td>Skeletonema marinoi</td><td>HE492-39 (UiO 452)</td><td>20 (27)</td><td>na</td><td>OK147691</td></tr><tr><td>Skeletonema marinoi</td><td>HE492-42 (UiO 455)</td><td>20 (27)</td><td>na</td><td>OK147690</td></tr></tbody></table>

(continued)
<table><caption>**Table 1. (continued)**</caption><thead><tr><th rowspan="2">Taxon</th><th rowspan="2">Strain Code (NORCCA)<sup>a</sup></th><th rowspan="2">Station<sup>b</sup></th><th colspan="2">GenBank Accession Number<sup>c</sup></th></tr><tr><th>18S</th><th>28S</th></tr></thead><tbody><tr><td colspan="5"><strong>Thalassiosiraceae</strong></td></tr><tr><td>Shionodiscus bioculatus</td><td>HE492-10 (UiO 438)</td><td>20 (27)</td><td>OK147668</td><td>OK147696</td></tr><tr><td>Shionodiscus bioculatus</td><td>HE492-23</td><td>7 (10)</td><td>OK147669</td><td>OK147695</td></tr><tr><td>Shionodiscus bioculatus</td><td>HE492-26</td><td>7 (10)</td><td>OK147670</td><td>OK147694</td></tr><tr><td>Shionodiscus bioculatus</td><td>HE492-31</td><td>6 (9)</td><td>OK147671</td><td>OK147693</td></tr><tr><td>Shionodiscus bioculatus</td><td>HE492-33</td><td>6 (9)</td><td>na</td><td>OK147692</td></tr><tr><td>Thalassiosira gravida</td><td>HE492-07 (UiO 435)</td><td>20 (27)</td><td>OK147672</td><td>OK147689</td></tr><tr><td>Thalassiosira gravida</td><td>HE492-08 (UiO 436)</td><td>20 (27)</td><td>OK147673</td><td>OK147688</td></tr><tr><td>Thalassiosira gravida</td><td>HE492-30 (UiO 447)</td><td>6 (9)</td><td>OK147674</td><td>OK147687</td></tr><tr><td>Thalassiosira gravida</td><td>HE492-32 (UiO 448)</td><td>6 (9)</td><td>OK147675</td><td>OK147686</td></tr><tr><td>Thalassiosira gravida</td><td>HE492-34</td><td>6 (9)</td><td>OK147676</td><td>OK147685</td></tr><tr><td>Thalassiosira oceanica</td><td>HE492-54</td><td>22 (29)</td><td>OK147677</td><td>OK147684</td></tr><tr><td>Thalassiosira sp.<sup>d</sup></td><td>HE492-45 (UiO 458)</td><td>21 (28)</td><td>OK147678</td><td>OK147683</td></tr><tr><td>Thalassiosira sp.<sup>d</sup></td><td>HE492-48 (UiO 461)</td><td>21 (28)</td><td>OK147679</td><td>OK147682</td></tr><tr><td>Thalassiosira sp.<sup>d</sup></td><td>HE492-52 (UiO 463)</td><td>22 (29)</td><td>OK147680</td><td>OK147681</td></tr><tr><td colspan="5"><strong>Corethraceae</strong></td></tr><tr><td>Corethron sp.</td><td>HE492-01</td><td>20 (27)</td><td>na</td><td>OK147704</td></tr></tbody></table>

<sup>a</sup>Strains deposited in the Norwegian Culture Collection of Algae (NORCCA) have been assigned an additional code (starting with "UiO").

<sup>b</sup>Station names in parentheses were used during the cruise and correspond to station names in the PANGEA publication (John and Wisotzki, 2017).

<sup>c</sup>All obtained 18S and 28S rRNA gene sequences have been deposited in the NCBI database and were assigned an accession number.

<sup>d</sup>Potentially novel undescribed species.

<sup>e</sup>Not available.

system was derived from the "16S Metagenomic Sequencing Library Preparation B" document (part no.15044223; Rev.B) distributed by Illumina with modifications for preparation of 18S ribosomal gene amplicons. Paired end Illumina sequencing (MSC 2.5.0.5/RTA 1.18.54, 2 x 300 bp) of 152 plankton samples was performed on a MiSeq platform (Illumina, United States). The sequence data were processed with the DADA2 pipeline for detecting amplicon sequence variants (ASVs) in a sample from the library of noisy reads generated by amplicon sequencing (Rosen et al., 2012; Callahan et al., 2016) as described in Elferink et al. (2020).

For further analysis, all diatom ASVs (Class = "Bacillariophyta"; N = 230) were extracted from the original ASV table based on the taxonomy assigned during the bioinformatic processing by DADA2 pipeline using PR² version 4.10.0 (https://github.com/pr2database/pr2database/releases/tag/4.10.0; Guillou et al., 2013) as the reference sequence database. The reduced, diatom-only ASV table and the corresponding taxonomical and environmental metadata were then incorporated into a phyloseq object (McMurdie and Holmes, 2013) in R used

in all of the downstream analysis and visualisation of the metabarcoding data.

### Phylogenetic analysis of diatom cultures

Sequences obtained from each of the 50 strains were first compared with the GenBank using the NCBI BLAST algorithm tool in Geneious (v. 11.1.5). After assigning the initial taxonomy based on the morphological observations and molecular (BLAST) results, the sequences of all the strains were aligned using the L-INS-i algorithm in mafft (v7.427; Stamatakis, 2014) to determine the number of unique 18S and 28S rRNA genotypes. For the 18S marker gene, the phylogenetic tree was prepared using selected morphologically verified reference sequences from PR² (v. 4.10.0) and recently published Arctic-diatom reference sequences (Balzano et al., 2017; Ribeiro et al., 2020) and reference sequences of *Chae-toceros* from Gaonkar et al. (2018). In addition, 18S rRNA sequences of cultured genotypes were compared with the local metabarcoding dataset to identify ASVs that were 100% identical to these strains. These ASVs were added to the 18S rRNA alignment using -add fragments algorithm in mafft and included in the 18S rRNA gene phylogenetic tree.
The alignment for the 18S rRNA gene tree included 32 sequences obtained from cultures in this study, 17 ASVs from the local metabarcoding dataset and 81 reference sequences (a total of 130 sequences and 2084 nucleotide positions). For the 28S rRNA tree, we used morphologically verified reference sequences from GenBank (a total of 117 sequences with 700 positions). For both marker genes, sequences were aligned using the L-INS-i algorithm in mafft, and Maximum Likelihood phylogenetic trees were made using RAxML (v8.2.12; GTRGAMMA model with 1000 bootstrap replications).

## Phylogenetic placement of diatom ASVs

To determine the phylogenetic placement of cultured genotypes and diatom ASVs obtained through metabarcoding, we compiled a diatom reference library with an improved representation of Arctic taxa. The diatom sequences (class=="Bacillariophyceae") were extracted from the PR² database (v.4.10.0; Guillou et al., 2013) and all the environmental ("ENV") diatom sequences were removed, leaving the total of 1965 sequences. To these sequences, we added 218 recently published 18S rRNA reference sequences of the genus *Chaetoceros* (Gaonkar et al., 2018), novel morphologically verified Arctic diatom reference sequences published by Balzano et al. (2017) (19 sequences) and Ribeiro et al. (2020) (33 sequences) (all recently included in the PR² v. 4.14.0), as well as 32 sequences newly obtained in this study. The initial reference alignment of the 2267 reference diatom sequences and three outgroup sequences (belonging to class Bolidophyceae, genus *Triparma*) was created using the FFT-NS-i algorithm with a maximum of 1000 iterations in mafft. After inspection of the alignment and trimming of the ends extending beyond the 1F-1528R primer positions (Medlin et al., 1988), the 230 diatom ASVs obtained in this study were placed in the alignment using the -add fragments algorithm in mafft. The final alignment containing reference sequences, our cultured genotypes and ASVs was used to build a Maximum Likelihood phylogenetic tree using the RAxML program (GTRGAMMA model with 1000 bootstrap replications). The final phylogenetic tree (Figure S1; interactive tree provided in ITOL: https://itol.embl.de/tree/193157255 105100381657613640) was then used to assign a phylogeny-based taxonomic identity at the lowest possible taxonomic level to each of the diatom ASVs, following best practice and guidelines for the phylogeny-based taxonomic assignment (Zimmermann et al., 2015; Gaonkar et al., 2020).

## Biogeographic analysis

For biogeographic analysis, we used a compilation of published metabarcoding datasets from different surveys that employed environmental metabarcoding targeting the V4 region of the 18S rRNA gene. The compiled dataset (metaPR²; Vaulot et al., 2022) has global coverage of coastal areas as it includes the samples from the Ocean Sampling Day survey (Kopf et al., 2015), notably in the northern hemisphere, as well as several oceanic transects based on the metabarcoding data collected within, e.g. the Tara

Ocean (de Vargas et al., 2015) and Malaspina expeditions (Duarte, 2015). The dataset has a good coverage of the Arctic realm, especially the Atlantic sector and Greenland, with both planktonic and sea-ice habitats. List of all datasets included in the compiled dataset and related metadata are provided in the GitHub repository (https://github.com/vaulot/metapr2_HE492). The map of the sampling stations included in the combined dataset is shown in **Figures 2–4**.

Metadata and sequence processing are detailed in Vaulot et al. (2022). In brief, raw sequence files and metadata were downloaded from NCBI SRA web site (https://www.ncbi.nlm.nih.gov/Traces/studies) when available or obtained directly from the investigators. For each study, raw sequences files were processed independently de novo. Primer sequences were removed using *cutadapt* (Martin, 2011) and amplicon processing was performed under the R software (R 2013) using the DADA2 package (Callahan et al., 2016). Taxonomic assignation of ASVs was performed against the PR² database (Guillou et al., 2013) version 4.12 (https://pr2-database.org; https://github.com/pr2database/pr2database/releases/tag/v4.12.0). ASV information was retrieved from final database available at https://shiny.metapr2.org.

For the biogeographic analysis of cultured genotypes, the 18S rRNA gene sequences of each genotype were compared with the ASVs found in the compiled dataset using the local BLAST tool in Geneious. After finding the ASVs in the dataset that were identical (100% pairwise identity) to the cultured genotypes, their distribution was mapped using an R script (available at https://vaulot.github.io/metapr2_HE492/metaPR2_HE492_diatoms.html), along with the relative contribution of each ASV to diatom assemblage at the sampling site. The same procedure was followed for mapping the distribution of the most abundant ASVs from the local metabarcoding dataset. After producing and examining the distribution maps, each genotype or ASV from the local metabarcoding dataset was assigned a descriptive biogeographic distribution type based on its global distribution pattern (**Tables 2 and S3**), namely Arctic (**Figure 2**), Arctic-temperate (**Figure 3**) or cosmopolitan (**Figure 4**).

## Results

### Overview of cultured diatom strains

In total, 50 diatom strains were isolated from Hausgarten (37 strains) and Wijdefjorden (13 strains) stations, most of them belonging to the families Chaetocerotaceae (22 strains) and Thalassiosiraceae (14 strains; **Table 1**). The remaining strains were isolated from families Bacillariaceae (5 strains), Cymatosiraceae (4 strains) and Skeleton-mataceae (2 strains), while Attheyaceae, Leptocylindraceae and Corethraceae were each represented with a single strain. The molecular analysis yielded a total of 49 28S rRNA gene sequences and 32 18S rRNA gene sequences. The phylogenetic analysis confirmed that the strains belonged to 23 distinct diatom genotypes. At least one marker gene was sequenced for each genotype, and the V4 region of the 18S rRNA gene was sequenced for 19 distinct genotypes (**Table 1**).
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Examples of genotypes with Arctic biogeographic distribution. Maps showing examples of Arctic biogeographic distribution based on the metaPR² database (right) along with distribution of example genotypes in the local metabarcoding dataset (left): (A) *Chaetoceros diadema* (ASV_786); (B) *Chaetoceros convolutus* (ASV_109); and (C) *Fragilariopsis* sp. (ASV_186). Symbol size is related to the relative abundance of genotype in respect to total diatom reads within each sample of the compiled metabarcoding dataset. Cross symbol (+) indicates samples in which the genotypes were not found.</figcaption></figure>

## Taxonomy and biogeographic distribution patterns of cultured diatom genotypes

Family Bacillariaceae

**Cylindrotheca closterium** strain HE492-63 (Wijdefjorden). Cells solitary, elongated, with two chloroplasts. Typical widening of the central part of the valve (**Figure 5A**). 18S rRNA gene sequence without a 100% match in GenBank (**Figure 6**). A 1 bp difference in the 18S rRNA gene sequence compared with the closest cultured strain KMMCC B-181 (South Korea). Sister clade to recently sequenced Arctic strains RCC5303, RCC5206 (Ribeiro et al., 2020) and RCC1985 (Balzano et al., 2017) (**Figures 6 and S1**). The 28S rRNA gene sequence identical to Strain KD10 (Kongsfjorden, Svalbard), a cold-adapted ecotype (Stock et al., 2019; **Figure 7**). Genotype found both in the local (ASV_766) and in the compiled dataset. Patchy local distribution. Higher read abundances in the innermost parts of all fjords except Isfjorden. Absent from the outermost Hausgarten stations (Figure S2A). Cosmopolitan distribution. Common in the Arctic realm, including

sea-ice communities in the Arctic Basin (Figure S2A). The 28S rRNA signature suggests an Arctic, cold-adapted population or a cryptic Arctic genotype that cannot be distinguished using the 18S rRNA V4 marker and light microscopy (LM)-based morphology.

*Fragilariopsis* sp. strain HE492-50 (Hausgarten). Cells solitary, not chain-forming, with two chloroplasts (**Figure 5B**). Valve lanceolate in trans-valval view and rectangular in girdle view. Characteristic row of costae clearly visible at the valve margin. The LM morphology insufficient for species-level identification. No 100% match in GenBank for either of the two markers. Most similar to various *Fragilariopsis* spp. and *Neodenticula seminae* (based on 28S rRNA; **Figure 7**) sequences. The ASV_186 matches the V4 sequence of this genotype, forming a clade with two more fairly abundant ASVs likely representing other unidentified *Fragilariopsis* species (Figure S1). ASV ubiquitous within the local dataset. Abundant in the northern parts of the studied area (Figure S2B). Globally, restricted to Arctic. No records in the Arctic Basin or in sea-ice samples.
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Examples of genotypes with Arctic-temperate biogeographic distribution. Maps showing examples of Arctic-temperate biogeographic distribution based on the metaPR² database (right) along with distribution of example genotypes in the local metabarcoding dataset (left): (A) *Skeletonema marinoi* (ASV_40); (B) *Thalassiosira gravida* (ASV_243); and (C) *Thalassiosira* sp. (ASV_117). Symbol size is related to the relative abundance of genotype in respect to total diatom reads within each sample of the compiled metabarcoding dataset. Cross symbol (+) indicates samples in which the genotypes were not found.</figcaption></figure>

**Pseudo-nitzschia granii** strain HE492-40 (Hausgarten). Relatively small, lanceolate cells, long chains with overlapping valve-tips, often solitary in culture (**Figure 5C**). Both 18S rRNA (**Figure 6**) and 28S rRNA (**Figure 7**) identical to sequences described by Balzano et al. (2017) from Beaufort Sea and UNC 1102 from North West Pacific. A 100% match with ASV_54. Relatively abundant at all stations in the local dataset (Figure S2C). Globally cosmopolitan distribution in coastal and oceanic areas. Common in the Atlantic sector of the Arctic and in western Greenland, but not in the Arctic Basin.

**Pseudo-nitzschia turgidula** strains HE492-37 and HE492-47 (Hausgarten). Cells relatively large, lanceolate. Chains typical for the genus *Pseudo-nitzschia* (**Figure 5D**). First record of this species in the Arctic. The 28S rRNA gene sequence identical to strain NWFSC-255 (Northeastern Pacific Ocean; Lundholm et al., 2012; **Figure 7**) identified as *P-n.* cf. *turgidula*. Two 100% matching 18S rRNA gene sequences (unknown origin) in the PR²

reference database identified as *P-n. turgidula* (**Figures 6** and S1). Matching ASV_179 most abundant at Hausgarten and present in Isfjorden and Kongsfjorden. Mostly absent from other fjords (Figure S2D). Globally, cosmopolitan but rare. Tropical coastal and oceanic stations and southern temperate and cold areas. Absent from the Arctic Basin. Could represent a warm-water genotype dispersed with Atlantic circulation.

**Attheya longicornis** strain HE492-59 (Hausgarten). Cells solitary, rectangular. Long undulating spines typical for this species projecting outwards diagonally from each valve end (**Figure 5E**). Phylogenetic placement within the *Attheya* clade. Unresolved distinction between *A. septentrionalis* and *A. longicornis*, possibly due to morphological misidentification of reference sequences. The 18S rRNA gene sequence identical to RCC2042 (Beaufort Sea; Balzano et al., 2017) and UNBF-P25C1 (Bay of Fundy, North Atlantic) identified as *A. septentrionalis* (**Figures 6** and S1). The 28S rRNA gene sequence matches several
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Examples of genotypes with cosmopolitan biogeographic distribution. Maps showing examples of cosmopolitan biogeographic distribution based on the metaPR² database (right) along with distribution of example genotypes in the local metabarcoding dataset (left): (A) *Shionodiscus bioculatus* (ASV_110); (B) *Pseudo-nitzschia granii* (ASV_54); and (C) *Thalassiosira oceanica* (ASV_187). Symbol size is related to the relative abundance of genotype in respect to total diatom reads within each sample of the compiled metabarcoding dataset. Cross symbol (+) indicates samples in which the genotypes were not found.</figcaption></figure>

<table><caption>**Table 2.** Descriptions of biogeographic distribution types as defined in this study</caption><thead><tr><td>Biogeographic Distribution Type</td><td>Description</td></tr></thead><tbody><tr><td>Arctic</td><td>Genotype has a biogeographic distribution mostly restricted to the Arctic circle (above 66°N; <strong>Figure 2</strong>).</td></tr><tr><td>Arctic-temperate</td><td>Genotype is present in the Arctic and at temperate latitudes (30°N–66°N and 30°S–66°S; <strong>Figure 3</strong>).</td></tr><tr><td>Cosmopolitan</td><td>Genotype is present in the Arctic and has a cosmopolitan distribution (<strong>Figure 4</strong>).</td></tr></tbody></table>

GenBank entries, including the above-mentioned RCC2042 from the Arctic (Figure 7). Matching ASV_570 is the most abundant ASV from the *Attheya* clade in the

local metabarcoding dataset. Abundant in Woodfjoden and Wijdefjorden, rare in the southern fjords and Hausgarten (Figure S2E). Arctic-temperate distribution. Clear Arctic presence including sea-ice samples.

**Family Cymatosiraceae**

*Arcocellulus cornucervis* strains HE492-38, -43 (Hausgarten) and HE492-71 (Wijdefjord). Saddle-like cells in girdle view. Characteristic pili extending inwards from the valve ends (Figure 5F). 18S (Figure 6) and 28S rRNA (Figure 7) gene sequences identical to strain RCC2270 (Beaufort Sea; Balzano et al., 2017), strain RCC5530 (Baffin Bay; Ribeiro et al., 2020), and strain NIES-3970 identified as *Minuto-cellus polymorphus* (unknown origin). In the reference phylogenetic tree, the 18S rRNA gene sequences and four related ASVs form a well-supported clade, indicating additional genotypes within this species present in the studied area (Figure S1). The most abundant ASV (ASV_24) from
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Light microscopy images of cultured genotypes. Scale bar = 10 µm in all images. (A) *Cylindrotheca closterium* (HE492-63); (B) *Fragilariopsis* sp. (HE492-50); (C) *Pseudo-nitzschia granii* (HE492-40); (D) *Pseudo-nitzschia turgidula* (HE492-37); (E) *Attheya longicornis* (HE492-59); (F) *Arcocellulus cornucervis* (HE492-43); (G) undetermined Cymatosiraceae (HE492-46); (H) *Chaetoceros contortus* (HE492-06); (I) *Chaetoceros convolutus* (HE492-64); (J) *Chaetoceros debilis* Clade 1 (HE492-61); (K) *Chaetoceros debilis* Clade 2 (HE492-28); (L) *Chaetoceros diadema* Clade 1 (HE492-13); (M) *Chaetoceros diadema* Clade 2 (HE492-17); (N) *Chaetoceros* sp. (HE492-05); (O) *Chaetoceros* sp. (HE492-36); (P) *Chaetoceros neogracilis* (HE492-73); (Q) *Skeletonema marinoi* (HE492-42); (R) *Leptocylindrus* sp. (HE492-14); (S) *Corethron* sp. (HE492-01); (T) *Thalassiosira oceanica* (HE492-54); (U) *Thalassiosira* sp. (HE492-45); (V) *Thalassiosira gravida* (HE492-07); and (W) *Shionodiscus bioculatus* (HE492-23).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 6.** The 18S rRNA phylogeny of cultured strains. Maximum-likelihood tree was generated in RAxML using the GTRGAMMAI model with 1000 bootstrap replications. Only bootstrap values > 70 are shown. The cultured strains are marked in red, ASVs obtained from the local metabarcoding dataset are marked in blue, and reference diatom sequences are marked in black. Numbers next to ASVs indicate their identifier and total read abundance in the metabarcoding dataset.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 7.** The 28S rRNA phylogeny of cultured strains. Maximum-likelihood tree was generated in RAxML using the GTRGAMMAI model with 1000 bootstrap replications. Only bootstrap values > 70 are shown. The cultured strains are marked in red, and reference diatom sequences are marked in black.</figcaption></figure>

the *Arcocellulus* clade identical to the cultured strains. Widely distributed. Higher abundance in the southern fjords (Van Mijendfjorden, Isfjorden). Lower abundance in the northern fjords and at Hausgarten (Figure S3A).

Arctic-temperate distribution. Present in Atlantic coastal areas, the Atlantic sector of the Arctic, and the Arctic Basin including sea-ice communities. Relatively rare beyond the polar front.
Unidentified strain belonging to the Cymatosiraceae family, HE492-46 (Hausgarten). Cells solitary. Elongated frustules with pointed ends, asymmetrical widening of valve along the apical axis (Figure 5G). Both 18S and 28S rRNA gene sequences have no match in GenBank (Figures 6 and 7). Novel clade within the Cymatosiraceae family. Three distinct ASVs from local metabarcoding dataset placed in the clade (Figure S1). ASV_46 identical to the cultured strain. Wide distribution in all fjords. Highest read abundance in Isfjorden (Figure S3B). Arctic-temperate distribution. Both coastal and open ocean presence. Not found in the Arctic Basin or the sea-ice samples.

## Family Chaetocerotaceae

Chaetoceros contortus strain HE492-06 (Hausgarten). Cells rectangular. Straight chains with setae emerging closer to the centre of the valve. Fused interlocking setae of neighbouring cells forming wide, hexagonal apertures (Figure 5H). Some apertures in a chain rotated by 90 degrees, dividing chains into segments with differently rotated apertures. Both 18S and 28S rRNA gene phylogeny placing the strain into the *C. contortus* clade. Both 28S and 18S rRNA gene sequences identical to strains Ch9A2, Ch12A4, Ch8A1 and Ch2B2 (Chile waters; Figures 6 and 7). Six ASVs found in the studied area placed in the *C. contortus* clade. Most abundant ASV_1089 identical to the cultured genotype (Figure S1). Higher read abundance in the inner parts of the northern fjords. Lower abundance in the outer Hausgarten stations (Figure S3C). Arctic-temperate distribution, prevalent in coastal samples.

Chaetoceros convolutus strain HE492-64 (Wijdefjorden). Large, heavily silicified species. Cells rectangular. Convex anterior valve, longer in the perivalvar axis than the apical axis. Chloroplasts in setae, tightly interlocking neighbouring cells leaving almost no apertures (Figure 5I). Chains slightly curved and twisted, the posterior terminal valve often has one large seta. No 100% matching 18S rRNA gene sequences in GenBank (Figure 6). Phylogenetical placement with other heavily silicified Chaetoceros species with chloroplasts in setae (subgenus Phaeoceros) such as *C. peruvianus* and *C. rostratus* (Figure S1). Matching ASV_1097 abundant in the northern fjords and at Hausgarten. Absent in the inner Isfjord and Van Mijenfjorden (Figure S3D). Predominantly found in the high Arctic, but also present in the Antarctic.

Chaetoceros debilis strain HE492-61(Hausgarten) belongs to Clade 1 and strains HE492-19 and -28 (Hausgarten) to Clade 2, based on 18S and 28S rRNA gene phylogenies (Figures 6 and 7), following Gaonkar et al. (2018). The Clade 1 strain has rectangular cells, almost identical length of pervalvar and apical axis. Thin setae emerging from the corners of the cells, slightly closer to the valve centre, fusing with neighbouring setae to form a wide hexagonal aperture (Figure 5J). Chains gently curved, spiralling and twisting along the chain axis. The Clade 2 strains have rectangular cells, shorter pervalvar axis compared with the apical axis, relatively narrow apertures formed by fusion of long basal parts of neighbouring setae (Figure 5K). In contrast to Clade 1 strain, the chains have pronounced spiral shape and are twisted along the

chain axis. Despite isolating two genotypes based on the 28S rRNA gene marker, a single ASV from this clade was found in the metabarcoding dataset, with identical sequence to the Clade 1 genotype (Figure S1). Locally high abundance at all stations except Isfjorden (Figure S3E). Arctic-temperate, including Arctic Basin.

Two genotypes of *Chaetoceros diadema* were isolated. Strain HE492-13 (Hausgarten) is placed in Clade 1 of this species based on the 28S rRNA phylogeny. The other nine strains (HE492-12, -15, -16, -17, -18, -21, -58, -65, -and -75), all but one (HE492-21 from Wijdefjorden) isolated at Hausgarten, form a new clade within the *C. diadema* species complex, based on both 28S and 18S rRNA gene phylogeny (Gaonkar et al., 2018; Figures 6 and 7). The Clade 1 strain has rectangular cells, pervalvar axis approximately twice as long as the apical axis (Figure 5L). Relatively long basal parts of the setae emerging from the corners of the cells, interlocking neighbouring setae forming rectangular, wide apertures. Chains relatively short, gently curved. Terminal setae extending diagonally outwards from the short basal parts. Short convex protuberance on the terminal valves. The novel clade has rectangular cells, apical axis longer than the pervalvar axis, as opposed to Clade 1 (Figure 5M). Valves slightly convex. Basal parts of neighbouring setae interlocking to form narrow windows. Chains relatively short and straight. Terminal setae projecting outwards diagonally then running parallel to each other. A single ASV_786 matching the novel clade was retrieved in the local metabarcoding dataset. Predominantly in the northern parts of the studied area (Figure S4A). Globally, mostly restricted to the Arctic, few minor observations below the Arctic circle.

Chaetoceros neogracilis strain HE492-73 (Wijdefjorden). Small single cells, rectangular (Figure 5P), short straight setae emerging diagonally from the corners of valves. 18S (Figure 6) and 28S rRNA gene sequences (Figure 7) placed into Clade 1 of the *C. neogracilis* complex (Balzano et al., 2017). Three ASVs from the local metabarcoding dataset placed in the *C. neogracilis* clade. ASV_121 identical to cultured genotype, but also to the Clade 2 reference sequences (Figure S1), such that the V4 region is identical for the two clades. Present at all stations. Highest abundance in Van Mijenfjorden in the south and the two northern fjords. Lower abundance at Hausgarten, Kongsfjorden and Isfjorden (Figure S4B). Arctic-temperate distribution. Abundant in the Arctic, including Arctic Basin and sea-ice.

Chaetoceros sp. strain HE492-36 (Wijdefjorden). Morphologically similar to *C. lauderi* and *C. teres*. Cells rectangular. Pervalval and apical axes of equal lengths (Figure 5N). Cells densely connected in the chain, overlapping setae, almost no apertures. Terminal setae extending outwards with a short basal part. Chains straight, not twisted along the chain axis. The 28S rRNA gene sequence has no 100% match in GenBank (Figure 7). The 18S rRNA gene sequence 99.9% identical to the strain CCAP 1010/16 (unknown origin; Figure 6). The phylogenetic placement indicates an undescribed, novel species of *Chaetoceros* affiliated to the *C. brevis/C. teres/C. lauderi* clade, with the highest similarity to *C. lauderi*
(Figure S1). A single ASV_1173 matched the cultured strain. Very rare, most abundant in the northern fjords and at two Hausgarten stations. Arctic-temperate distribution (Figure S4C).

**Chaetoceros** sp. strains HE492-02, -03, -05, -11 and -27 (Hausgarten). A novel species and a new clade within the *Chaetoceros* genus. Cells rectangular and elongated (**Figure 50**). Setae from neighbouring cells merging at the relatively short basal parts, leaving small, elliptical apertures. Chains straight, long, terminal setae extending almost laterally in respect to the chain axis. The 28S rRNA marker gene has no 100% match in the GenBank. The 18S rRNA sequence identical to a *Chaetoceros* sp. strain (coast of Japan; **Figures 6 and 7**). Phylogenetic placement as a sister group to the *C. anastomosans/C. vixvisibilis* clade, though with low bootstrap support (Figure S1). ASV_1516 identical to cultured strains, Only detected at Hausgarten (Figure S4D). Arctic-temperate distribution. Rare in the Arctic Basin.

### Family Skeletonemataceae

**Skeletonema marinoi**, strains HE492-39 and HE492-42 (Hausgarten). Cells relatively small, rectangular in girdle view, connected with interlocking processes growing from the edges of valves, forming long characteristic chains (**Figure 5Q**). Only 28S rRNA gene sequences obtained, identical to *S. marinoi* (**Figure 7**). A single ASV (ASV_40) belonging to *Skeletonema marinoi* was found in the metabarcoding dataset, but not certain is whether it matches the 18S rRNA gene sequences of the cultured genotype, as it may represent other *Skeletonema* species, too (Figure S1). Read abundance highest in Kongsfjorden, but high in other fjords as well. Slightly lower abundance at Hausgarten (Figure S4E). Globally very common, Arctic-temperate distribution, including sea-ice samples.

### Family Thalassiosiraceae

**Shionodiscus bioculatus** strains HE492-10, -23, -26, -31 and -33 (Hausgarten and Wijdefjorden). Cells cylindrical, often solitary, sometimes forming shorter chains of cells closely linked with numerous organic threads (**Figure 5W**). The 18S (**Figure 6**) and 28S (**Figure 7**) rRNA marker genes all match *S. bioculatus* strains RCC5532 and RCC1991 (Beaufort Sea and Baffin Bay, respectively; Balzano et al., 2017). ASV_110 had an identical sequence to the cultured strain. Present at all stations. Higher abundance at Hausgarten and outer parts of fjords (Figure S5A). Cosmopolitan, including high Arctic and Antarctic.

**Thalassiosira gravida** strains HE492-07, -08, -30, -32 and -34 (Hausgarten and Wijdefjorden). Cells chain-forming, cylindrical, rectangular in girdle view with slightly curved edges of valves (**Figure 5V**). Single thick thread connecting cells in a chain. Numerous strutted processes visible in girdle view. Both 18S (**Figures 6 and 1** S1) and 28S (**Figure 7**) rRNA marker genes match *T. gravida* and *T. rotula* genotypes. These were considered the same species with temperature-induced morphological variability (Sar et al., 2011), though some studies suggest that they are genetically distinct based on their ITS sequences (Whittaker et al., 2012). A single ASV_243 was

identical to the cultured strains. Present at all stations. Most abundant in the innermost parts of Woodfjorden, Wijdefjorden and Van Mijenfjorden (Figure S5B). Arctic-temperate distribution, few records in Arctic Basin.

**Thalassiosira oceanica** strain HE492-54 (Hausgarten). Small and cylindrical, mostly solitary cells, slightly rounded valve edges in girdle view (**Figure 5T**). The 18S (**Figure 6**) and 28S (**Figure 7**) rRNA marker genes match previously sequenced representatives of this species. Phylogenetic placement in 18S rRNA tree is within a well-supported *T. oceanica* clade (Figure S1). The ASV_187 identical to the cultured genotype. Found at almost all stations. Highest read abundance along the Hausgarten-Kongsfjorden transect (Figure S5C). Cosmopolitan distribution, rare in the Arctic Basin.

The strains HE492-45, -48, and -52 (Hausgarten) belong to a new clade within the **Thalassiosira** genus. Cells exclusively solitary, relatively small and elongated, cylindrical with roundish valve edges (**Figure 5U**). Single, long organic thread extending diagonally from each valve. No matching 28S or 18S rRNA gene sequences in the GenBank. Phylogenetic placement based on 18S rRNA gene suggests that the novel clade is a sister group to *T. concaviuscula*, though with low bootstrap support (Figure S1). The matching ASV_117 from the metabarcoding dataset found in almost all samples. Most abundant in Hausgarten, outer Kongsfjorden and outer stations of northern fjords (Figure S5D). Arctic-temperate distribution, excluding Arctic Basin.

### Family Corethraceae

**Corethron** sp. strain HE492-01 (Hausgarten). Cells with rounded valves, elliptical appearance in girdle view. Numerous long spines projecting outwards from bases of both valves (**Figure 5S**). Numerous small round chloroplasts. Only 28S rRNA gene sequence obtained (**Figure 7**). A single 100% matching sequence in the NCBI, belonging to *Corethron* sp. (Bay of Fundy). Only one ASV was retrieved from the *Corethron* clade in the local metabarcoding dataset, which could not be linked with certainty to the cultured strain due to lack of 18S rRNA sequence. The ASV was present in all fjords. Arctic-temperate distribution excluding Arctic Basin (Figure S5E).

### Family Leptocylindraceae

**Leptocylindrus** sp. strain HE492-14 (Hausgarten). Cells cylindrical, elongated, typical morphology for the genus, either solitary or chain-forming (**Figure 5R**). Numerous small round chloroplasts. Only 28S rRNA sequence obtained. No matching sequences in NCBI (**Figure 7**). A total of 21 ASVs placed within the *Leptocylindrus* clade, 20 of which were forming a clade together with *L. minimus* reference sequences, suggesting either a high undescribed diversity within this genus or a high degree of haplotype diversity within one population. As the 28S sequence of the cultured strain did not place within the *L. danicus/L. aporus* or *L. hargreavesi* clades (**Figure 7**), the strain may belong to the *L. minimus* clade for which there is no 28S rRNA reference sequence available, but which had a relatively high read abundance in the area. Further studies
and 18S rRNA sequencing will be needed to identify this
strain and link it with the environmental sequences.

Diatom diversity and distribution in Svalbard fjords
based on the metabarcoding data

A total of 230 diatom ASVs were obtained after the
processing of the raw metabarcoding data. Following
the phylogenetic placement of the ASVs into the newly
compiled reference library and the construction of the
RAxML tree, 110 ASVs representing 48% of total diatom
ASVs could be assigned taxonomically down to the spe-
cies level with high certainty (Figure S1). The remaining
52% (120 ASVs) could not be placed at the species level
due to a lack of reference sequences. On the other hand,
genus-level taxonomic identity could be assigned to
87 additional ASVs, thus in total 197 (86%) of diatom
ASVs, with only 33 ASVs not being assigned to the
genus level. These mainly belonged to different raphid
and araphid pennate diatoms (21 ASVs) and polar-
centric Mediophyceae (11). Overall, the taxonomic
assignment to species level was above average within
ecologically important and taxonomically well-studied
genera of large centric diatoms such as *Chaetoceros*,
where 16 out of 26 ASVs (62%) were assigned to species
level, or *Thalassiosira*, with 13 out of 19 (68%) ASVs
assigned to the species level.

Overall, diatoms had relatively low contribution to the
protist community in the studied area based on the
metabarcoding data, reaching the maximum of approxi-
mately 25% contribution of sequence reads in Isfjorden
(Figure S11). Their diversity and distribution patterns,
inferred from the metabarcoding data, varied signifi-
cantly among different fjords (Figures 8A and S12). At
Hausgarten stations, the diatom community was domi-
nated by the two *Pseudo-nitzschia* species (*P-n. turgidula*
and *P-n. granii*), along with *S. bioculatus*, *Thalassiosira*
spp. and *Eucampia groenlandica*. In Van Mijenfjorden,
the community was dominated by *A. cornucervis*, as well
as *Skeletonema* spp. (most notably *S. marinoi*), *Thalassiosira*
spp. and *Chaetoceros* spp. The community in Isfjor-
den was unique within the investigated area, with a high
abundance of different ASVs belonging to the *L. minimus*
clade, except at the outermost station, which was dom-
inated by *Pseudo-nitzschia* spp. *A. cornucervis* and *Skele-
tonema* spp. were also significant components of the
Isfjorden community. In Kongsfjorden, the community
was primarily dominated by *Skeletonema* spp., along with
*A. cornucervis* and *Thalassiosira* spp. Finally, the two
northernmost fjords, Woodfjorden and Wijdefjorden,
exhibited similar communities dominated by heavily
silicified *Chaetoceros* species (e.g., *C. convolutus*), *Thalas-
siosira* spp. (especially *T. hispida*), *E. groenlandica* and
*Actinocyclus* sp.

Biogeographic distribution patterns of the most
abundant non-cultured ASVs from the local
metabarcoding dataset

To explore the biogeographic patterns of diatoms at the
community level using the local metabarcoding dataset,
we performed the additional biogeographic analysis of the

25 most abundant ASVs in the dataset that did not corre-
spond to our cultured genotypes (Figures S6–S10; Table
S3). When combined with ASVs that matched the cultured
strains 100%, these accounted for 97% of total diatom
read abundance in the dataset. Of the 25 most abundant
non-cultured ASVs, 11 could not be assigned taxonomi-
cally at the species level, and one raphid pennate ASV
could not be identified at genus level. The biogeographic
patterns observed among non-cultured ASVs were similar
to those of the cultured strains. The most common bio-
geographic distribution type was Arctic-temperate (18),
followed by Arctic (7). No genotypes from this group
exhibited a cosmopolitan distribution.

The ASVs exhibiting Arctic-temperate biogeographic
distribution belonged to *L. minimus* (6 ASVs), *E. groenlan-
dica*, *Thalassiosira concaviuscula*, *Chaetoceros brevis*, *Cha-
toceros gelidus*, *Proboscia alata*, unidentified species of
*Actinocyclus*, *Chaetoceros*, *Fragilariopsis*, *Proboscia*, *Rhizo-
solenia*, *Skeletonema* and an unknown raphid pennate. The
ASVs with Arctic distribution were placed both within
well-known Arctic diatom species such as *T. hispida*, *Fra-
gilariopsis cylindrus* and *C. convolutus*, and within uniden-
tified species belonging to genera *Skeletonema*,
*Chaetoceros* and *Bacterosira*.

Biogeographic patterns at the community level

The biogeographic types assigned to 20 cultured dia-
tom genotypes (ASVs) and 25 most abundant diatom
ASVs from the local metabarcoding dataset were
explored at the community level to explore the possible
relationship between degree of Atlantic influence and
prevalence of biogeographic distribution types in the
Svalbard area (Figure 8B). The Hausgarten stations had
a unique composition of ASVs, mainly dominated by
genotypes with cosmopolitan (approximately 25–40%)
and Arctic-temperate (approximately 25–50%) distribu-
tion. The prevalence of Arctic genotypes at this station
was relatively low (up to 20%). The outermost station in
Isfjorden showed a pattern that was similar to Hausgar-
ten stations, with cosmopolitan genotypes accounting
for approximately 50% of total read abundance. The
relative read abundance of Arctic (up to about 10%)
and cosmopolitan (up to 5%) genotypes was generally
low in Isfjorden, where Arctic-temperate genotypes
made up > 90% of reads. The cosmopolitan genotypes
were generally absent from the innermost parts of
Isfjorden.

A very similar contribution of different biogeographic
types was observed in Van Mijenfjorden and Kongsfjor-
den. Here, the most abundant group were ASVs with
Arctic-temperate distribution (up to approximately
80%), followed by the species that are biogeographically
restricted to the Arctic. A slight gradient was observed
along these fjords in terms of biogeographic distribution
types, with Arctic genotypes increasing, and cosmopoli-
tan genotypes decreasing in contribution towards the
innermost stations. Finally, the two northernmost fjords,
namely Woodfjorden and Wijdefjorden, had a large pro-
portion of genotypes with Arctic distribution compared
with the other studied fjords and Hausgarten stations.
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Diversity and biogeography of diatoms in the Svalbard fjords. (A) Diatom diversity in the studied area shown at the genus level. Only ASVs that constituted more than 10% of total diatom ASVs at each station are included in the plot; (B) Distribution of diatom genotypes in the Svalbard fjords based on their biogeographic type. Each bar segment represents relative abundance of an individual ASV. The ASVs labeled as NA have not been included in the biogeographic analysis.</figcaption></figure>

The contribution of ASVs restricted to the Arctic was increasing towards the inner parts of the fjords (reaching up to 50% of total read abundance), followed by a decrease in the prevalence of the species with Arctic-temperate and cosmopolitan biogeographic distribution. The species with cosmopolitan distribution had a relatively minor contribution (up to about 10%), and were nearly absent from the innermost stations in Woodfjorden and Wijdefjorden.

## Discussion

*Bridging the gap between morphological and molecular approaches in studying Arctic diatom diversity*

A well-established taxonomic framework, integrating both morphological and molecular approaches, is a prerequisite for studying the composition and structure of Arctic diatom communities and understanding how they shape the Arctic ecosystem function. In this work, we used
morphological and molecular characterisation of diatom cultures to link robust taxonomic identification to local and global metabarcoding datasets. The precise identification at the species and genotype level allowed for detecting general biogeographic patterns that are underlying diatom community structure in the Svalbard area.

Large-scale cultivation efforts are crucial for bridging the gap between traditional morphology-based identification of diatoms and state-of-the-art molecular metabarcoding approaches. The phylogenetic analysis of cultured and environmental genotypes shows that Arctic diatom diversity is not yet fully represented in morphological and molecular databases, as demonstrated in recent large-scale diatom cultivation studies in the Arctic (Balzano et al., 2017; Ribeiro et al., 2020). Here, we have found seven diatom strains that could not be identified to species level using light microscopy or available molecular reference libraries. Probable novel species were found even within the taxonomically well studied planktonic genera such as *Chaetoceros* and *Thalassiosira*. We also found a cryptic species and a novel sub-clade within a globally distributed *C. diadema* species complex that is biogeographically restricted to the Arctic realm. These findings underline the importance of further taxonomic studies based on molecular methods and suggest that Arctic diatom communities harbour a significant undescribed diversity of species or genotypes that are specifically adapted to the extreme Arctic environment and endemic to the Arctic.

The taxonomy assignment of ASVs from the local metabarcoding datasets confirmed the notion that the interpretative power of metabarcoding is directly linked to the representation of taxa in reference libraries (Zimmermann et al., 2015; Gaonkar et al., 2020). In our analysis, about half of the diatom ASVs could be identified to the species level, and this level of identification was only possible for groups with high-quality reference sequences such as the genus *Chaetoceros* (Gaonkar et al., 2018) or other well-studied planktonic genera like *Pseudo-nitzschia* and *Thalassiosira*. Overall, Arctic diatoms are underrepresented in reference databases, and large-scale cultivation is needed to link morphological and molecular information. Notable gaps in reference databases could be observed for pennate diatoms, where a significant portion of the ASVs could not be identified to species or even genus level. In particular, our study did not include sea-ice or benthic communities, where pennate diatoms are especially ecologically relevant and diverse, and where gaps in reference databases are likely to be even higher. The 23 cultured genotypes accounted for about a third (35%) of the total diatom sequence reads in the metabarcoding dataset, showing that combining cultivation with metabarcoding surveys can significantly improve the subsequent interpretation of metabarcoding data.

We used two common genetic markers for species-level identification and phylogenetic placement of diatom cultures, namely the 18S and 28S rRNA gene sequences. The two showed overall similar taxonomic resolution and proved to be sufficient for discerning different diatom morphotypes, as was demonstrated earlier (Zimmermann et al., 2015; Gaonkar et al., 2020). The clear advantage of

18S rRNA, apart from being a well-established marker for metabarcoding of diatoms (Zimmermann et al., 2011), is in the availability of reference databases such as PR² (Guillou et al., 2013) and of large metabarcoding datasets such as Tara Ocean (de Vargas et al., 2015), the Ocean Sampling Day (Kopf et al., 2015) and Malaspina (Duarte, 2015) that allow for global biogeographic mapping of ASVs, as shown here and in similar studies on diatoms (De Luca et al., 2019) and other protists (Ichinomiya et al., 2016; Lopes dos Santos et al., 2017). In some cases, the 18S rRNA marker was not sufficient for discerning morphospecies, as was shown for the two strains of *C. debilis* that could only be separated using the 28S rRNA sequences. Similarly, the *T. gravida* cultures had identical 18S and partial 28S rRNA gene sequences to *T. rotula* sequences, as the two can only be discerned using more variable genes such as ITS rRNA (Whittaker et al., 2012).

In this study, the V4 region of the 18S rRNA gene marker appeared to have sufficient resolution for species-level identification of diatoms in most cases. This sufficiency means that the local and global compiled 18S rRNA metabarcoding datasets can be used to determine biogeographic distribution patterns at the genotype level, likely reflecting the ecological niche of a genotype. Clearly, the use of more variable markers such as ITS rRNA, *rbcL* and *psbC* that can differentiate additional diatom genotypes may reveal finer biogeographic patterns. However, metabarcoding datasets based on these markers are lacking, and 18S rRNA gene datasets presently provide the best global coverage.

In some cases, we found a high number of ASVs within certain terminal diatom clades. This finding was most notable within the *L. minimus* clade, but also within other clades such as *C. contortus*, *A. cornucervis* or the new strain from the Cymatosiraceae family. This high diversity of ASVs may indicate high intragenomic haplotype diversity of the 18S rRNA gene, as was shown in environmental metabarcoding data from the genus *Chaetoceros* (Gaonkar et al., 2020) and within clonal strains of the genus *Chaetoceros* (De Luca et al., 2021), where dominant haplotype is accompanied with peripheral haplotypes with order of magnitude lower read numbers. High intragenomic haplotype diversity was also reported for the genus *Skeletonema* (Alverson and Kolnick, 2005). Another argument for the common genomic origin of dominant and peripheral ASVs within *L. minimus* and other clades is the similarity of distribution patterns among the most abundant ASVs in *L. minimus* clade, both within the local dataset and globally. In light of these observations, applying a 100% pairwise identity when linking cultures with metabarcoding datasets likely reveals only the patterns of the dominant haplotype which is also the one obtained by Sanger sequencing of cultures. The dominant haplotype is nevertheless highly representative of a studied genotype, and the analysis is not likely losing significant biogeographic information from minor haplotypes. On the other hand, differences in the global distribution among some of the peripheral ASVs would suggest that they indeed represent distinct species that are currently lacking in the reference databases.
## Biogeographic patterns of Arctic diatoms

Marine planktonic protists, including diatoms, have generally been considered cosmopolitan (globally distributed), with unlimited potential for dispersal through oceanic circulation, and their community structure primarily regulated by environmental forcing (Finlay, 2002). In particular, the population structure of marine planktonic diatoms has been argued to be defined almost exclusively by environmental gradients with no long-term limitations to dispersal that could result in biogeographically restricted or endemic planktonic diatom taxa (Cermeño and Falkowski, 2009). Despite recent findings suggesting that endemic diatoms are common in freshwater systems (Vanormelingen et al., 2009), notably in polar regions (Vyverman et al., 2010), there has been no strong evidence for endemism in marine planktonic diatoms. Still, early morphology-based studies (Hasle, 1976) have shown that planktonic diatoms can exhibit certain biogeographic patterns such as cold-water, warm-water, and cosmopolitan distribution.

In recent years, findings of cryptic cold-adapted, planktonic diatom species within presumably cosmopolitan species complexes (Chamnansin et al., 2013; Percopo et al., 2016) have suggested that the morphology-based taxonomy and morphology-based datasets may, in many cases, prove insufficient for detecting biogeographic patterns and potential endemism or restricted distribution in diatoms. Integrating morphology-based records with global molecular metabarcoding datasets, such as Tara Ocean (Karlusich et al., 2020) and Ocean Sampling Day (Kopf et al., 2015), revealed that genotypes of *Chaetoceros* show both endemic and biogeographically restricted distribution patterns, notably in polar areas (De Luca et al., 2019). Building on this research, we demonstrate here that genotypes of Arctic planktonic diatoms exhibit a variety of biogeographic patterns, ranging from endemic (or Arctic-restricted distribution) to cosmopolitan distribution.

After analysing biogeographic patterns of 45 diatom genotypes from the Svalbard area, we detected three distinct general biogeographic distribution types, namely Arctic, Arctic-temperate and cosmopolitan. Overall, there was no clear link between phylogeny and biogeography, and species within the most relevant diatom genera such as *Chaetoceros* and *Thalassiosira* exhibited a wide range of biogeographic distribution types. These distribution types generally correspond well with the global biogeographic distribution of Arctic diatoms defined by Balzano et al. (2017) and show some similarities with the earlier study of Hasle (1976), both of which used morphology-based biogeographic analyses. For example, the biogeographic patterns delineated by Balzano et al. (2017) and this study were the same for *A. cornucervis* (Arctic-temperate), *C. closterium* (cosmopolitan), *E. groenlandica* (Arctic-temperate), *T. gravida* (Arctic-temperate), *C. gelidus* (Arctic-temperate) and *C. neogracilis* (Arctic-temperate). On the other hand, slightly different biogeographic distribution was observed in this study compared with that of Balzano et al. (2017) for *P.-n. granii* (cosmopolitan in this study versus Arctic-temperate), *T. hispida* (Arctic versus Arctic-

temperate), *S. bioculatus* (cosmopolitan versus Arctic-temperate). Observed differences in biogeography of individual species between this study and solely morphology-based studies can have different causes. They may stem from the finer, genotype-level approach used here that could distinguish cold-adapted cryptic species, but also in some cases from insufficient resolution of 18S rRNA marker when distinguishing morphologically distinct but closely related species, as is often the case in the *Pseudonitzschia* genus (Lundholm et al., 2012).

The most common biogeographic distribution type among the studied genotypes was Arctic-temperate, which included 30 genotypes that are present both at temperate latitudes (30–66°N) and above the polar circle (66°N). Within this general biogeographic group that was dominating diatom communities in the fjords on the west coast of Svalbard (at Hausgarten, Van Mijenfjorden, Isfjorden and Kongsfjorden), we observed a range of possible finer biogeographic patterns that were difficult, however, to delineate objectively based on 18S V4 rRNA metabarcoding datasets used in this study. For example, some Arctic-temperate genotypes (e.g., *C. neogracilis* and *A. longicornis*) have clearly established communities in the Arctic Basin (beyond the polar front) and are also found in sea-ice communities. On the other hand, Arctic-temperate genotypes such as *C. contortus*, *Corethron* sp., or some genotypes of *L. minimus*, are common in the Svalbard area, but are not present beyond the polar front. These Arctic-temperate genotypes, absent from the Arctic Basin, are of particular interest as they are seemingly mainly associated with relatively warmer Atlantic waters (> 3 °C), and their northwards expansion can indicate Atlantification or ocean warming in the area.

Similar patterns could be observed among the cosmopolitan genotypes. Out of the five cosmopolitan genotypes, only two were present continuously at all latitudes between the Arctic Basin and Antarctica (*S. bioculatus*, and *C. closterium*). The remaining three, *P.-n. granii*, *P.-n. turgidula* and *T. oceanica*, were generally absent from the Arctic Basin. In particular, *P.-n. turgidula* was mostly found at Hausgarten, suggesting that it represents a warm-water genotype associated with Atlantic waters. These observations also indicate that finer biogeographic patterns and additional cryptic diversity would likely be revealed for these genotypes using finer genetic markers, as demonstrated by detection of cold-adapted *C. closterium* ecotypes based on 28S rRNA gene (Stock et al., 2019).

Ten genotypes were endemic to the Arctic. These belonged to species which have already been considered as restricted to the Arctic (and sometimes to southern polar areas as well) based on morphology-based studies, such as *C. convolutus* and *F. cylindrus* (Lundholm and Hasle, 2010; Balzano et al., 2017), but also to some cryptic genotypes within species complexes with wider biogeographic distribution such as *C. diadema*, *T. hispida* and *S. costatum*. This work demonstrates that biogeographically restricted distribution or Arctic endemism not only exists but is relatively common among planktonic diatoms in the Svalbard area. Similar studies of sea
ice-associated diatom genotypes would undoubtedly
reveal additional endemism among Arctic diatoms,
which has currently been suspected or assumed for
a number of sea ice-associated pennate diatoms
(Quillfeldt, 2000; Poulin et al., 2011).

Examination of biogeographic distribution types on the
community level within five Svalbard fjords and Hausgar-
ten area revealed distinct patterns in the abundance of
genotypes with specific distribution types. Overall, diatom
communities were shaped by the three main groups of
genotypes: 1) genotypes restricted to the Arctic, that dom-
inated in the northernmost fjords (Woodfjorden and Wij-
defjorden); 2) genotypes with Arctic-temperate
distribution, that dominated fjords on the west coast of
Svalbard (Van Mijenfjorden, Isfjorden and Kongsfjorden),
and 3) genotypes with cosmopolitan distribution, that
dominated at Hausgarten stations. Within groups 2 and
3, there seemed to exist additional distribution patterns,
with differing degrees of expansion within the Arctic
realm. This finding suggests that the diatom communities
in the area are structured by an interplay of several pro-
cesses: dispersal and advection of warm water genotypes
and their ability to survive the extreme Arctic environ-
ment (especially beyond the polar front), ability of geno-
types restricted to the Arctic to expand southwards into
warmer waters, and an underlying stock of genotypes
which seem to be well-adapted both to warmer waters
and the Arctic.

Advection is a well-documented source of highly pro-
ductive, warm, and nutrient-rich Atlantic waters into Fram
Strait and the western Svalbard area (Randelhoff et al.,
2015; Wassmann et al., 2015; Wassmann et al., 2019).
An important component of advection is the transport
of planktonic organisms into the Arctic realm, a phenom-
enon that has been well-documented for temperate meso-
zooplankton (Wassmann et al., 2019). To date, the main
evidence for the transport of phytoplankton with Atlantic
waters was the northward expansion of the haptophyte,
*Emiliania huxleyi*, reported in several studies from Fram
Strait and the Barents Sea (Hegseth and Sundfjord, 2008;
Neukermans et al., 2018; Oziel et al., 2020). Another
example is the northward expansion in the Atlantic Arctic
of the cyanobacterium *Synechococcus* (Paulsen et al.,
2016). The results from this study demonstrate for the first
time that diatom genotypes with broad Arctic-temperate
and cosmopolitan distribution are abundant in the Sval-
bard area, especially in the areas with stronger Atlantic
influence.

While the local metabarcoding data shown here pro-
vide only a snapshot of the late Arctic summer diatom
community, they clearly demonstrate that species with
wider biogeographic distribution dominated the
Atlantic-influenced waters, and genotypes endemic to the
Arctic dominated colder waters. Despite the clear general
trend, the finer scale genotype-level patterns provided
a much more complex picture. For example, the cosmo-
politan *P.-n. turgidula* genotype was not detected in Van
Mijenfjorden and was nearly absent from the northern-
most fjords, while another cosmopolitan genotype *P.-n.
granii* was common in all fjords. Moreover, while the

majority of the Arctic-temperate genotypes were present
in all fjords, most of the Arctic-temperate *L. minimus* gen-
otypes were nearly absent from the innermost parts of the
northern fjords. On the same note, the *F. cylindrus* geno-
type, endemic for the Arctic, was uniformly distributed
throughout the area. These patterns point out that
genotype-level ecophysiology is also important for shap-
ing diatom community structure, as it likely determines
the local distribution and abundance of individual geno-
types with respect to specific hydrography of each fjord
and seasonal variation in environmental drivers (Basedow
et al., 2018; Menze et al., 2020).

Conclusion

This work underscores the importance of combining mor-
phological and molecular approaches in studying biodiver-
versity and biogeography of Arctic diatoms. By providing
novel morphology-based data and taxonomic descriptions
as well as novel reference sequences of 18S and 28S rRNA
marker genes, it represents an important contribution to
future studies on Arctic diatom biodiversity. Moreover, we
demonstrated that Arctic diatom flora harbours a signifi-
cant undescribed biodiversity, and that new taxa are also
found within morphologically well-described genera of
planktonic diatoms. Future studies should focus on for-
mally describing these new taxa and expanding on the
intensive cultivation efforts to further reveal the biodiver-
sity of Arctic diatoms. Our biogeographic analysis shows
that diatom genotypes from the Svalbard area exhibit
a range of biogeographic patterns, from endemic geno-
types restricted to the Arctic, to genotypes with wider,
Arctic-temperate and cosmopolitan biogeographic distri-
bution. Using finer taxonomic markers (e.g., ITS rDNA),
capable of delineating cryptic Arctic species and cold-
adapted Arctic genotypes, would likely reveal additional,
finer scale biogeographic patterns, providing better
insight into biogeographic processes that are shaping the
Arctic community structure.

Data accessibility statement

CTD data from the cruise is available on Pangea link: John
and Wisotzki (2017): Physical oceanography during
HEINCKE cruise HE492; https://doi.org/10.1594/
PANGAEA.881306. Raw Ilumina metabarcoding data from
the HE492 cruise have been deposited to GenBank SRA
under accession number: PRJEB49358. Sequences of the
strains are available at NCBI with accession numbers
OK147649-147680 for 18S rRNA sequences, and
OK147681-147729 for 28S rRNA sequences (Table 1).
Scripts and processed data are available from https://
github.com/vaulot/metapr2_HE492.

**Supplemental files**

The supplemental files for this article can be found as
follows:

Figure S1 is provided both as a supplemental file and as
an interactive phylogenetic tree in ITOL: https://itol.embl.
de/tree/193157255105100381657613640. Supplemen-
tary PDF document including 11 supplementary figures
(Figure S2–S12) and three supplementary tables (Table S1–S3) is submitted with the document.

## Acknowledgments

We are thankful to the crew on board R/V Heincke for their assistance during the HE492 cruise in Svalbard fjords. Jan Heuschele (UiO) provided valuable help with coding in R and data analysis. Bioinformatic analyses were performed on the UiO UNINETT Saga2 computer cluster, project nn9307k, with help from Anders Kristian Krabberod (UiO). We thank the Roscoff ABIMS platform of the FR2424 (CNRS, Sorbonne Université) for bioinformatics resources.

## Funding

This study was financially supported by The Research Council of Norway (RCN) to the projects "TaxMArc: Taxonomy of marine protists in a changing Arctic" no. 268286 (LS, BE, UJ, DV) and by ERA-MBT initiated project "PROMISE" (SGS and BE no. 282936). UJ was supported by the "PROMISE" project no. 161B0561A from DLR/BMBF. UJ and KK were supported by ERA4CS, an ERA-NET initiated by JPI Climate, and funded DLR/BMBF (DE): 01LS1714A with co-funding by the European Union (Grant 690462). UJ was DFG supported to SPP1991, "Taxonomics" project GZ: JO 702/13-1 AOBJ: 641030. UJ received funding and support for the HE492 cruise within the Helmholtz research program "Changing Earth, Sustaining our Future" (Sub-topic 6.2 Adaptation of marine life) of the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Germany. DV was supported by ANR contract "PhytoPol" (ANR-15-CE02-0007).

## Competing interests

The authors have declared that no competing interests exist.

## Author contributions

Contributed to conception and design: LS, BE, UJ.

Contributed to acquisition of data: LS, KK, SGS, CH, BE, UJ.

Contributed to analysis and interpretation of data: LS, SGS, DV, BE, UJ.

Drafted and/or revised the article: LS, KK, SGS, CH, DV, BE, UJ.

Approved the submitted version for publication: LS, KK, SGS, CH, DV, BE, UJ.

## References

Adakudlu, M, Andresen, J, Bakke, J, Beldring, S, Benestad, R, Bilt, W, van der Bogen, J, Borstad, CP, Breili, K, Breivik, Ø, Børsheim, KY, Christiansen, HH, Dobler, A, Engeset, R, Frauenfelder, R, Gerland, S, Gjelten, HM, Gundersen, J, Isaksen, K, Jaedicke, C, Kierulf, H, Kohler, J, Li, H, Lutz, J, Melvold, K, Mezghani, A, Nilsen, F, Nilsen, IB, Nilsen, JEØ, Pavlova, O, Ravndal, O, Risebrobakken, B, Saloranta, T, Sandven, S, Schuler, T, Simpson, MJR, Skogen, MD, Smedsrud, LH, Sund, M, Vikhamar-Schuler, D, Westermann, S, Wong,

WK. 2019. Climate in Svalbard 2100. The Norwegian Centre for Climate Services (NCCS). Available at https://bora.uib.no/bora-xmlui/handle/1956/19136. Accessed 19 March 2021.

Alverson, AJ, Kolnick, L. 2005. Intragenomic nucleotide polymorphism among small subunit (18S) rDNA paralogs in the diatom genus Skeletonema (Bacillariophyta). *Journal of Phycology* **41**(6): 1248–1257. DOI: http://dx.doi.org/10.1111/j.1529-8817.2005.00136.x.

Balzano, S, Percopo, I, Siano, R, Gourvil, P, Chanoine, M, Marie, D, Vaulot, D, Sarno, D. 2017. Morphological and genetic diversity of Beaufort Sea diatoms with high contributions from the *Chaetoceros neo-gracilis* species complex. *Journal of Phycology* **53**(1): 161–187. DOI: http://dx.doi.org/10.1111/jpy.12489.

Basedow, SL, Sundfjord, A, Von Appen, W-J, Halvorsen, E, Kwasniewski, S, Reigstad, M. 2018. Seasonal variation in transport of zooplankton into the Arctic Basin through the Atlantic gateway, Fram Strait. *Frontiers in Marine Science* **5**. DOI: http://dx.doi.org/10.3389/fmars.2018.00194.

Bradley, IM, Pinto, AJ, Guest, JS. 2016. Design and evaluation of Illumina MiSeq-compatible, 18S rRNA gene-specific primers for improved characterization of mixed phototrophic communities. *Applied Environmental Microbiology* **82**(19): 5878–5891. DOI: http://dx.doi.org/10.1128/AEM.01630-16.

Callahan, BJ, McMurdie, PJ, Rosen, MJ, Han, AW, Johnson, AJA, Holmes, SP. 2016. DADA2: High-resolution sample inference from Illumina amplicon data. *Nature Methods* **13**(7): 581–583. DOI: http://dx.doi.org/10.1038/nmeth.3869.

Cermeño, P, Falkowski, PG. 2009. Controls on diatom biogeography in the ocean. *Science* **325**(5947): 1539–1541. DOI: http://dx.doi.org/10.1126/science.1174159.

Chamnansinp, A, Li, Y, Lundholm, N, Moestrup, O. 2013. Global diversity of two widespread, colony-forming diatoms of the marine plankton, *Chaetoceros socialis* (syn. *C. radians*) and *Chaetoceros gelidus* sp. nov. *Journal of Phycology* **49**(6): 1128–1141. DOI: http://dx.doi.org/10.1111/jpy.12121.

Cottier, FR, Nilsen, F, Skogseth, R, Tverberg, V, Skarðhamar, J, Svendsen, H. 2010. Arctic fjords: A review of the oceanographic environment and dominant physical processes. *Geological Society, London, Special Publications* **344**(1): 35–50. DOI: http://dx.doi.org/10.1144/SP344.4.

Darling, KF, Kucera, M, Wade, CM. 2007. Global molecular phylogeography reveals persistent Arctic circumpolar isolation in a marine planktonic protist. *Proceedings of the National Academy of Sciences* **104**(12): 5002–5007. DOI: http://dx.doi.org/10.1073/pnas.0700520104.

De Luca, D, Kooistra, WHCF, Sarno, D, Biffali, E, Pir edda, R. 2021. Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus *Chaetoceros*. *Scientific Reports* **11**(1):
807. DOI: http://dx.doi.org/10.1038/s41598-020-80829-6.

De Luca, D, Kooistra, WHCF, Sarno, D, Gaonkar, CC, Piredda, R. 2019. Global distribution and diversity of *Chaetoceros* (Bacillariophyta, Mediophyceae): Integration of classical and novel strategies. *PeerJ* 7: e7410. DOI: http://dx.doi.org/10.7717/peerj.7410.

De Vargas, C, Audic, S, Henry, N, Decelle, J, Mahé, F, Logares, R, Lara, E, Berney, C, Le Bescot, N, Probert, I, Carmichael, M, Poulain, J, Romac, S, Colin, S, Aury, J-M, Bittner, L, Chaffron, S, Dunthorn, M, Engelen, S, Flegontova, O, Guidi, L, Horák, A, Jaillon, O, Lima-Mendez, G, Lukeš, J, Malviya, S, Morard, R, Mulot, M, Scalco, E, Siano, R, Vincent, F, Zingone, A, Dimier, C, Picheral, M, Searson, S, Kandels-Lewis, S, Acinas, S G, Bork, P, Bowler, C, Gorsky, G, Grimsley, N, Hingamp, P, Iudicone, D, Not, F, Ogata, H, Pesant, S, Raes, J, Sieracki, M E, Speich, S, Stemmann, L, Sunagawa, S, Weissenbach, J, Wincker, P, Karsenti, E. 2015. Eukaryotic plankton diversity in the sunlit ocean. *Science* **348**(6237): 1261605. DOI: http://dx.doi.org/10.1126/science.1261605.

Duarte, CM. 2015. Seafaring in the 21St Century: The Malaspina 2010 Circumnavigation expedition. *Limnology and Oceanography Bulletin* **24**(1): 11–14. DOI: http://dx.doi.org/https://doi.org/10.1002/lob.10008.

Eamer, J, Donaldson, GM, Gaston, J, Kosobokova, KN, Lárusson, KF, Melnikov, I, Reist, JD, Richardson, E, Staples, L, Von Quillfeldt, CH. 2013. Life linked to ice: A guide to sea-ice-associated biodiversity in this time of rapid change. CAFF Assessment Series No. 10 Report No.: 9789935431257. Available at https://oaarchive.arctic-council.org/bitstream/handle/11374/239/Life_Linked_to_Ice_Oct_2013.pdf?sequence=1&isAllowed=yhttp://www.caff.is/assessment-series/254-life-linked-to-ice-a-guide-to-sea-ice-associated-biodiversity-in-this-time-of-ra/download. Accessed 19 March 2021.

Edvardsen, B, Shalchian-Tabrizi, K, Jakobsen, KS, Medlin, LK, Dahl, E, Brubak, S, Paasche, E. 2003. Genetic variability and molecular phylogeny of *Dinophysis* species (Dinophyceae) from Norwegian waters inferred from single cell analyzes of ribosomal DNA. *Journal of Phycology* **39**(2): 395–408. DOI: http://dx.doi.org/10.1046/j.1529-8817.2003.01252.x.

Elferink, S, John, U, Neuhaus, S, Wohlrab, S. 2020. Functional genomics differentiate inherent and environmentally influenced traits in dinoflagellate and diatom communities. *Microorganisms* **8**(4): 567. DOI: http://dx.doi.org/10.3390/microorganisms8040567.

Eppley, RW, Holmes, RW, Strickland, JDH. 1967. Sink- ing rates of marine phytoplankton measured with a fluorometer. *Journal of Experimental Marine Biol- ogy and Ecology* **1**(2): 191–208. DOI: http://dx.doi. org/10.1016/0022-0981(67)90014-7.

Finlay, BJ. 2002. Global dispersal of free-living microbial eukaryote species. *Science* **296**(5570): 1061–1063. DOI: http://dx.doi.org/10.1126/science.1070710.

Fragoso, GM, Poulton, AJ, Yashayaev, IM, Head, EJH, Johnsen, G, Purdie, DA. 2018. Diatom biogeogra- phy from the Labrador Sea revealed through a trait- based approach. *Frontiers in Marine Science* **5**. DOI: http://dx.doi.org/10.3389/fmars.2018.00297.

Gaonkar, CC, Piredda, R, Minucci, C, Mann, DG, Montresor, M, Sarno, D, Kooistra, WHCF. 2018. Anno- tated 18S and 28S rDNA reference sequences of taxa in the planktonic diatom family Chaetocerotaceae. *PLoS One* **13**(12): e0208929. DOI: http://dx.doi. org/10.1371/journal.pone.0208929.

Gaonkar, CC, Piredda, R, Sarno, D, Zingone, A, Montresor, M, Kooistra, WHCF. 2020. Species detec- tion and delineation in the marine planktonic diatoms *Chaetoceros* and *Bacteriastrum* through metabarcoding: Making biological sense of haplo- type diversity. *Environmental Microbiology* **22**(5): 1917–1929. DOI: http://dx.doi.org/10.1111/1462- 2920.14984.

Guillou, L, Bachar, D, Audic, S, Bass, D, Berney, C, Bittner, L, Boutte, C, Burgaud, G, de Vargas, C, Decelle, J, del Campo, J, Dolan, JR, Dunthorn, M, Edvardsen, B, Holzmann, M, Kooistra, WHC F, Lara, E, Le Bescot, N, Logares, R, Mahé, F, Mas- sana, R, Montresor, M, Morard, R, Not, F, Paw- lowski, J, Probert, I, Sauvadet, AL, Siano, R, Stoeck, T, Vaulot, D, Zimmermann, P, Christen, R. 2013. The protist ribosomal reference database (PR²): A catalog of unicellular eukaryote small sub- unit rRNA sequences with curated taxonomy. *Nucleic Acids Research* **41**(Database issue): 597–604. DOI: http://dx.doi.org/10.1093/nar/ gks1160.

Haine, TWN, Curry, B, Gerdes, R, Hansen, E, Karcher, M, Lee, C, Rudels, B, Sreen, G, de Steur, L, Stew- art, KD, Woodgate, R. 2015. Arctic freshwater export: Status, mechanisms, and prospects. *Global and Planetary Change* **125**: 13–35. DOI: http://dx. doi.org/10.1016/j.gloplacha.2014.11.013.

Hasle, GR. 1976. The biogeography of some marine plank- tonic diatoms. *Deep Sea Research and Oceanographic Abstracts* **23**: 319–338.

Hegseth, EN, Sundfjord, A. 2008. Intrusion and bloom- ing of Atlantic phytoplankton species in the high Arctic. *Journal of Marine Systems* **74**(1–2): 108–119. DOI: http://dx.doi.org/10.1016/j.jmarsys. 2007.11.011.

Hegseth, EN, Tverberg, V. 2013. Effect of Atlantic water inflow on timing of the phytoplankton spring bloom in a high Arctic fjord (Kongsfjorden, Sval- bard). *Journal of Marine Systems* **113–114**: 94–105. DOI: http://dx.doi.org/10.1016/j.jmarsys. 2013.01.003.

Hop, H, Vihtakari, M, Bluhm, BA, Assmy, P, Poulin, M, Gradinger, R, Peeken, I, von Quillfeldt, C, Olsen, LM, Zhitina, L, Melnikov, IA. 2020. Changes in sea- ice protist diversity with declining sea ice in the
Arctic Ocean from the 1980s to 2010s. *Frontiers in Marine Science* **7**. DOI: http://dx.doi.org/10.3389/fmars.2020.00243.

Ichinomiya, M, Dos Santos, AL, Gourvil, P, Yoshikawa, S, Kamiya, M, Ohki, K, Audic, S, De Vargas, C, Noël, M-H, Vaulot, D, Kuwata, A. 2016. Diversity and oceanic distribution of the Parmales (Bolidophyceae), a picoplanktonic group closely related to diatoms. *The ISME Journal* **10**(10): 2419–2434. DOI: http://dx.doi.org/10.1038/ismej.2016.38.

John, U, Wisotzki, A. 2017. *Physical oceanography during HEINCKE cruise HE492*. Bremerhaven, Pangaea: Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research. DOI: http://dx.doi.org/10.1594/PANGAEA.881306.

Karlusich, JJP, Ibarbalz, FM, Bowler, C. 2020. Phytoplankton in the Tara Ocean. *Annual Review of Marine Science* **12**(1): 233–265. DOI: http://dx.doi.org/10.1146/annurev-marine-010419-010706.

Kopf, A, Bicak, M, Kottmann, R, Schnetzer, J, Kostadinov, I, Lehmann, K, Fernandez-Guerra, A, Jeanthon, C, Rahav, E, Ullrich, M, Wichels, A, Gerdts, G, Polymenakou, P, Kotoulas, G, Siam, R, Abdallah, RZ, Sonnenschein, EC, Cariou, T, O’Gara, F, Jackson, S, Orlic, S, Steinke, M, Busch, J, Duarte, B, Caçador, I, Canning-Clode, J, Bobrova, O, Marteinsson, V, Reynisson, E, Loureiro, CM, Luna, GM, Quero, GM, Löscher, CR, Kremp, A, DeLorenzo, ME, Øvreås, L, Tolman, J, LaRoche, J, Penna, A, Frischer, M, Davis, T, Katherine, B, Meyer, CP, Ramos, S, Magalhães, C, Jude-Lemeilleur, F, Aguirre-Macedo, ML, Wang, S, Poulton, N, Jones, S, Collin, R, Fuhrman, JA, Conan, P, Alonso, C, Stambler, N, Goodwin, K, Yakimov, MM, Baltar, F, Bodrossy, L, Van De Kamp, J, Frampton, DM, Ostrowski, M, Van Ruth, P, Malthouse, P, Claus, S, Deneudt, K, Mortelmans, J, Pitois, S, Wallom, D, Salter, I, Costa, R, Schroeder, DC, Kandil, MM, Amaral, V, Biancalana, F, Santana, R, Pedrotti, ML, Yoshida, T, Ogata, H, Ingleton, T, Munnik, K, Rodriguez-Ezpeleta, N, Berteaux-Lecellier, V, Wecker, P, Cancio, I, Vaulot, D, Bienhold, C, Ghazal, H, Chaouni, B, Essayeh, S, Ettamimi, S, Zaid, EH, Boukhatem, N, Bouali, A, Chahboune, R, Barrijal, S, Timinouni, M, El Otmani, F, Bennani, M, Mea, M, Todorova, N, Karamfilov, V, ten Hoopen, P, Cochrane, G, L’Haridon, S, Bizsel, KC, Vezzi, A, Lauro, FM, Martin, P, Jensen, RM, Hinks, J, Gebbels, S, Rosselli, R, De Pascale, F, Schiavon, R, dos Santos, A, Villar, E, Pesant, S, Cataletto, B, Malfatti, F, Edirisinghe, R, Silveira, JAH, Barbier, M, Turk, V, Tinta, T, Fuller, W J, Salihoglu, I, Serakinci, N, Ergoren, MC, Bresnan, E, Iriberri, J, Nyhus, PAF, Bente, E, Karlsen, HE, Golyshin, PN, Gasol, JM, Moncheva, S, Dzhembekova, N, Johnson, Z, Sinigalliano, CD, Gidley, ML, Zingone, A, Danovaro, R, Tsiamis, G, Clark, MS, Costa, AC, El Bour, M, Martins, A M, Collins, RE, Ducluzeau, AL, Martinez, J, Costello, MJ, Amaral-Zettler, LA, Gilbert, JA, Davies, N, Field, D, Glöckner, FO.

2015. The ocean sampling day consortium. *GigaScience* **4**(1): 27. DOI: http://dx.doi.org/10.1186/s13742-015-0066-5.

Krause, JW, Duarte, CM, Marquez, IA, Assmy, P, Fernández-Méndez, M, Wiedmann, I, Wassmann, P, Kristiansen, S, Agustí, S. 2018. Biogenic silica production and diatom dynamics in the Svalbard region during spring. *Biogeosciences* **15**(21): 6503–6517. DOI: http://dx.doi.org/10.5194/bg-15-6503-2018.

Kwok, R. 2018. Arctic sea ice thickness, volume, and multiyear ice coverage: Losses and coupled variability (1958–2018). *Environmental Research Letters* **13**(10): 105005. DOI: http://dx.doi.org/10.1088/1748-9326/aae3ec.

Le Moigne, FAC, Poulton, AJ, Henson, SA, Daniels, CJ, Fragoso, GM, Mitchell, E, Richier, S, Russell, BC, Smith, HEK, Tarling, GA, Young, JR, Zubkov, M. 2015. Carbon export efficiency and phytoplankton community composition in the Atlantic sector of the Arctic Ocean. *Journal of Geophysical Research: Oceans* **120**(6): 3896–3912. DOI: http://dx.doi.org/10.1002/2015JC010700.

Leeuwe, M, Tedesco, L, Arrigo, KR, Assmy, P, Campbell, K, Meiners, KM, Rintala, JM, Selz, V, Thomas, DN, Stefels, J. 2018. Microalgal community structure and primary production in Arctic and Antarctic sea ice: A synthesis. *Elementa: Science of the Anthropocene* **6**(1). DOI: http://dx.doi.org/10.1525/elementa.267.

Leu, E, Mundy, CJ, Assmy, P, Campbell, K, Gabrielsen, TM, Gosselin, M, Juul-Pedersen, T, Gradinger, R. 2015. Arctic spring awakening—Steering principles behind the phenology of vernal ice algal blooms. *Progress in Oceanography* **139**: 151–170. DOI: http://dx.doi.org/10.1016/j.pocean.2015.07.012.

Lopes dos Santos, A, Gourvil, P, Tragin, M, Noël, M-H, Decelle, J, Romac, S, Vaulot, D. 2017. Diversity and oceanic distribution of prasinophytes clade VII, the dominant group of green algae in oceanic waters. *The ISME Journal* **11**(2): 512–528. DOI: http://dx.doi.org/10.1038/ismej.2016.120.

Lovejoy, C, Legendre, L, Price, N. 2002. Prolonged diatom blooms and microbial food web dynamics: Experimental results from an Arctic polynya. *Aquatic Microbial Ecology* **29**: 267–278. DOI: http://dx.doi.org/10.3354/ame029267.

Lovejoy, C, Vincent, WF, Bonilla, S, Roy, S, Martineau, MJ, Terrado, R, Potvin, M, Massana, R, Pedrós-Alió, C. 2007. Distribution, phylogeny, and growth of cold-adapted picoprasinophytes in Arctic Sea. *Journal of Phycology* **43**(1): 78–89. DOI: http://dx.doi.org/10.1111/j.1529-8817.2006.00310.x.

Lundholm, N, Bates, SS, Baugh, KA, Bill, BD, Connell, LB, Léger, C, Trainer, VL. 2012. Cryptic and pseudo-cryptic diversity in diatoms—with descriptions of *Pseudo-nitzschia hasleana* sp. nov. and *P. fryxelliana* sp. nov. 1. *Journal of Phycology* **48**(2): 436–454. DOI: http://dx.doi.org/https://doi.org/10.1111/j.1529-8817.2012.01132.x.
Lundholm, N, Hasle, GR. 2010. Fragilariopsis (Bacillariophyceae) of the Northern Hemisphere—Morphology, taxonomy, phylogeny and distribution, with a description of *F. pacifica* sp. nov. *Phycologia* **49**(5): 438–460. DOI: http://dx.doi.org/10.2216/09-97.1.

Martin, M. 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. *EMBnet. Journal* **17**(1): 10. DOI: http://dx.doi.org/10.14806/ej.17.1.200.

McMurdie, PJ, Holmes, S. 2013. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. *PLoS ONE* **8**(4): e61217. DOI: http://dx.doi.org/10.1371/journal.pone.0061217.

Medlin, L, Elwood, HJ, Stickel, S, Sogin, ML. 1988. The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. *Gene* **71**(2): 491–499. DOI: http://dx.doi.org/10.1016/0378-1119(88)90066-2.

Menze, S, Ingvaldsen, RB, Nikolopoulos, A, Hattermann, T, Albretsen, J, Gjøsæter, H. 2020. Producitive detours—Atlantic water inflow and acoustic backscatter in the major troughs along the Svalbard shelf. *Progress in Oceanography* **188**: 102447. DOI: http://dx.doi.org/10.1016/j.pocean.2020.102447.

Neukermans, G, Oziel, L, Babin, M. 2018. Increased intrusion of warming Atlantic water leads to rapid expansion of temperate phytoplankton in the Arctic. *Global Change Biology* **24**(6): 2545–2553. DOI: http://dx.doi.org/10.1111/gcb.14075.

Nilsen, F, Skogseth, R, Vaardal-Lunde, J, Inall, M. 2016. A simple shelf circulation model: Intrusion of Atlantic Water on the West Spitsbergen Shelf. *Journal of Physical Oceanography* **46**(4): 1209–1230. DOI: http://dx.doi.org/10.1175/JPO-D-15-0058.1.

Orkney, A, Platt, T, Narayanaswamy, BE, Kostakis, I, Bouman, HA. 2020. Bio-optical evidence for increasing *Phaeocystis* dominance in the Barents Sea. *Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences* **378**(2181): 20190357. DOI: http://dx.doi.org/10.1098/rsta.2019.0357.

Overland, J, Dunlea, E, Box, JE, Corell, R, Forsius, M, Kattsov, V, Olsen, MS, Pawlak, J, Reiersen, L-O, Wang, M. 2019. The urgency of Arctic change. *Polar Science* **21**: 6–13. DOI: http://dx.doi.org/10.1016/j.polar.2018.11.008.

Oziel, L, Baudena, A, Ardyna, M, Massicotte, P, Randelhoff, A, Sallée, JB, Ingvaldsen, RB, Devred, E, Babin, M. 2020. Faster Atlantic currents drive poleward expansion of temperate phytoplankton in the Arctic Ocean. *Nature Communications* **11**(1): 1–8. DOI: http://dx.doi.org/10.1038/s41467-020-15485-5.

Paulsen, ML, Doré, H, Garczarek, L, Seuthe, L, Müller, O, Sandaa, RA, Bratbak, G, Larsen, A. 2016. Synechococcus in the Atlantic Gateway to the Arctic Ocean. *Frontiers in Marine Science* **3**: 191. DOI: http://dx.doi.org/10.3389/fmars.2016.00191.

Pavlov, AK, Tverberg, V, Ivanov, BV, Nilsen, F, Falk-Petersen, S, Granskog, MA. 2013. Warming of Atlantic water in two west Spitsbergen fjords over the last century (1912–2009). *Polar Research* **32**(1): 11206. DOI: http://dx.doi.org/10.3402/polar.v32i0.11206.

Percopo, I, Ruggiero, MV, Balzano, S, Gourvil, P, Lundholm, N, Siano, R, Tammilehto, A, Vaulot, D, Sarno, D. 2016. *Pseudo-nitzschia arctica* sp. nov., a new cold-water cryptic *Pseudo-nitzschia* species within the *P. pseudodelicatissima* complex. *Journal of Phycology* **52**(2): 184–199. DOI: http://dx.doi.org/https://doi.org/10.1111/jpy.12395.

Polyakov, IV, Pnyushkov, AV, Alkire, MB, Ashik, IM, Baumann, TM, Carmack, EC, Goszczko, I, Guthrie, J, Ivanov, VV, Kanzow, T, Krishfield, R, Kwok, R, Sundfjord, A, Morison, J, Rember, R, Yulin, A. 2017. Greater role for Atlantic inflows on sea-ice loss in the Eurasian Basin of the Arctic Ocean. *Science* **356**(6335): 285–291. DOI: http://dx.doi.org/10.1126/science.aai8204.

Poulin, M, Daugbjerg, N, Gradinger, R, Ilyash, L, Ratkova, T, von Quillfeldt, C. 2011. The pan-Arctic biodiversity of marine pelagic and sea-ice unicellular eukaryotes: A first-attempt assessment. *Marine Biology* **41**(1): 13–28. DOI: http://dx.doi.org/10.1007/s12526-010-0058-8.

Quillfeldt, CHV. 2000. Common diatom species in Arctic spring blooms: Their distribution and abundance. *Botanica Marina* **43**(6): 499–516. DOI: http://dx.doi.org/10.1515/BOT.2000.050.

Randelhoff, A, Sundfjord, A, Reigstad, M. 2015. Seasonal variability and fluxes of nitrate in the surface waters over the Arctic shelf slope. *Geophysical Research Letters* **42**(9): 3442–3449. DOI: http://dx.doi.org/https://doi.org/10.1002/2015GL063655.

Ribeiro, CG, Dos Santos, AL, Gourvil, P, Le Gall, F, Marie, D, Tragin, M, Probert, I, Vaulot, D. 2020. Culturable diversity of Arctic phytoplankton during pack ice melting. *Elementa: Science of the Anthropocene* **8**(6). DOI: http://dx.doi.org/10.1525/elementa.401.

Ribeiro, CG, Dos Santos AL, Probert, I, Vaulot, D, Edvardsen, B. 2020. Taxonomic reassignment of *Pseudohaptolina birgeri comb. nov.* (Haptophyta). *Phycologia* **59**(6): 606–615. DOI: http://dx.doi.org/10.1080/00318884.2020.1830255.

Rosen, MJ, Callahan, BJ, Fisher, DS, Holmes, SP. 2012. Denoising PCR-amplified metagenome data. *BMC Bioinformatics* **13**(1): 283. DOI: http://dx.doi.org/10.1186/1471-2105-13-283.

Sar, EA, Sunesen, I, Lavigne, AS, Lofeudo, S. 2011. *Thalassiosira rotula*, a heterotypic synonym of *Thalassiosira gravida*: Morphological evidence. *Diatom Research* **26**(1): 109–119. DOI: http://dx.doi.org/10.1080/0269249X.2011.573691.

Scholin, CA, Herzog, M, Sogin, M, Anderson, DM. 1994. Identification of group- and strain-specific genetic markers for globally distributed *Alexandrium* (Dinophyceae). II. Sequence analysis of a fragment of the
LSU rRNA gene. *Journal of Phycology* **30**(6): 999-1011. DOI: http://dx.doi.org/10.1111/j.0022-3646.1994.00999.x.

Skarołhamar, J, Svendsen, H. 2010. Short-term hydro-
graphic variability in a stratified Arctic fjord. *Geolog-
ical Society, London, Special Publications* **344**(1):
51–60. DOI: http://dx.doi.org/10.1144/SP344.5.

Soltwedel, T, Bauerfeind, E, Bergmann, M, Bracher, A,
Budaeva, N, Busch, K, Cherkasheva, A, Fahl, K,
Grzelak, K, Hasemann, C, Jacob, M, Kraft, A,
Lalande, C, Metfies, K, Nöthig, EM, Meyer, K,
Quéric, NV, Schewe, I, Włodarska-Kowalczuk,
M, Klages, M. 2016. Natural variability or
anthropogenically-induced variation? Insights from
15 years of multidisciplinary observations at the
Arctic marine LTER site Hausgarten. *Ecological Indi-
cators* **65**: 89–102. DOI: http://dx.doi.org/10.1016/
j.ecolind.2015.10.001.

Stamatakis, A. 2014. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. *Bioinformatics* **30**(9): 1312–3. DOI: http://dx.doi.org/10.1093/bioinformatics/btu033.

Stock, W, Vanselslander, B, Rüdiger, F, Sabbe, K, Vyver-
man, W, Karsten, U. 2019. Thermal niche differen-
tiation in the benthic diatom *Cylindrotheca
closterium* (Bacillariophyceae) Complex. *Frontiers in
Microbiology* **10**: 1395. DOI: http://dx.doi.org/10.
3389/fmicb.2019.01395.

Thomas, MK, Kremer, CT, Klausmeier, CA, Litchman, E. 2012. A global pattern of thermal adaptation in marine phytoplankton. *Science* **338**(6110): 1085-1088. DOI: http://dx.doi.org/10.1126/science.1224836.

Tverberg, V, Skogseth, R, Cottier, F, Sundfjord, A, Walczowski, W, Inall, ME, Falck, E, Pavlova, O, Nilsen, F. 2019. The Kongsfjorden transect: Seasonal and inter-annual variability in hydrography. In Hop H, Wiencke C eds. *The ecosystem of Kongsfjorden, Svalbard*. Cham, Switzerland: Springer International Publishing: 49–104. DOI: http://dx.doi.org/10.1007/978-3-319-46425-1_3.

van Leeuwe, MA, Tedesco, LMV, Arrigo, KR, Assmy, P
Campbell, K, Meiners, KM, Rintala, JM, Selz, V,
Thomas, DN, Stefels, J. 2018. Microalgal commu-
nity structure and primary production in Arctic and
Antarctic sea ice: A synthesis. *Elementa: Science of
the Anthropocene* 6(1): 4. DOI: http://dx.doi.org/10.
1525/elementa.267.

Vanormelingen, P, Verleyen, E, Vyverman, W. 2009. The diversity and distribution of diatoms: From cosmopolitanism to narrow endemism. In Foissner, W, Hawksworth, DL eds., *Protist diversity and geographical distribution*. Dordrecht, the Netherlands: Springer: 159–171. DOI: http://dx.doi.org/10.1007/978-90-481-2801-3_12.

Vaulot, D, Sim, CWH, Ong, D, Teo, B, Biwer, C, Jamy, M, Lopes dos Santos, A. 2022. MetaPR²: A database of eukaryotic 18S rRNA metabarcodes with an emphasis on protists. *Molecular Ecology Resources* **00**:

1-14. DOI: http://dx.doi.org/10.1111/1755-0998.
13674.

Vyverman, W, Verleyen, E, Wilmotte, A, Hodgson, DA,
Willems, A, Peeters, K, Van de Vijver, B, De
Wever, A, Leliaert, F, Sabbe, K. 2010. Evidence for
widespread endemism among Antarctic micro-
organisms. *Polar Science* **4**(2): 103–113. DOI:
http://dx.doi.org/10.1016/j.polar.2010.03.006.

Wassmann, P, Kosobokova, KN, Slagstad, D, Drinkwater, KF, Hopcroft, RR, Moore, SE, Ellingsen, I, Nelson, RJ, Carmack, E, Popova, E, Berge, J. 2015. The contiguous domains of Arctic Ocean advection: Trails of life and death. *Progress in Oceanography* **139**: 42–65. DOI: http://dx.doi.org/10.1016/j.pcean.2015.06.011.

Wassmann, P, Slagstad, D, Ellingsen, I. 2019. Advection of mesozooplankton into the Northern Svalbard shelf region. *Frontiers in Marine Science* **6**. DOI: http://dx.doi.org/10.3389/fmars.2019.00458.

Weydmann, A, Søreide, JE, Kwaśniewski, S, Leu, E, Falk-Petersen, S, Berge, J. 2013. Ice-related seasonality in zooplankton community composition in a high Arctic fjord. *Journal of Plankton Research* **35**(4): 831–842. DOI: http://dx.doi.org/10.1093/plankt/fbt031.

Whittaker, KA, Rignanese, DR, Olson, RJ, Rynearson, TA. 2012. Molecular subdivision of the marine diatom *Thalassiosira rotula* in relation to geographic distribution, genome size, and physiology. *BMC Evolutionary Biology* **12**(1): 209. DOI: http://dx.doi.org/10.1186/1471-2148-12-209.

Wiedmann, I, Ceballos-Romero, E, Villa-Alfageme, M,
Renner, AHH, Dybwad, C, Jagt, H, Svensen, C,
Assmy, P, Wiktor, JM, Tatarek, A, Różańska-
Pluta, M, Iversen, MH. 2020. Arctic observations
identify phytoplankton community composition as
driver of carbon flux attenuation. *Geophysical
Research Letters* **47**(14). DOI: http://dx.doi.org/10.
1029/2020GL087465.

Yau, S, Lopes dos Santos, A, Eikrem, W, Gérikas
Ribeiro, C, Gourvil, P, Balzano, S, Escande, M,
Moreau, H, Vaulot, D. 2020. *Mantoniella beaufortii* and *Mantoniella baffinensis* sp. nov. (Mamiellales, Mamiellophyceae), two new green algal species from the high Arctic. *Journal of Phycology* **56**(1): 37–51. DOI: http://dx.doi.org/10.1111/jpy.12932.

Zimmermann, J, Glückner, G, Jahn, R, Enke, N, Gemeinholzer, B. 2015. Metabarcoding vs. morphological identification to assess diatom diversity in environmental studies. *Molecular Ecology Resources* **15**(3): 526–542. DOI: http://dx.doi.org/10.1111/1755-0998.12336.

Zimmermann, J, Jahn, R, Gemeinholzer, B. 2011. Bar-coding diatoms: Evaluation of the V4 subregion on the 18S rRNA gene, including new primers and protocols. *Organisms Diversity and Evolution* **11**(3): 173–192. DOI: http://dx.doi.org/10.1007/s13127-011-0050-6.
**How to cite this article:** Šupraha, L, Klemm, K, Gran-Stadniczeňko, S, Hörstmann, C, Vaulot, D, Edvardsen, B, John, U. 2022. Diversity and biogeography of planktonic diatoms in Svalbard fjords: The role of dispersal and Arctic endemism in phyto-plankton community structuring. *Elementa: Science of the Anthropocene* 10(1). DOI: https://doi.org/10.1525/elementa.2021.00117

**Domain Editor-in-Chief:** Jody W. Deming, University of Washington, Seattle, WA, USA

**Associate Editor:** Kevin R. Arrigo, Department of Earth System Science, Stanford University, Stanford, CA, USA

**Knowledge Domain:** Ocean Science

**Published:** August 9, 2022

**Accepted:** July 12, 2022

**Submitted:** December 9, 2021

**Copyright:** © 2022 The Author(s). This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License (CC-BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. See http://creativecommons.org/licenses/by/4.0/.

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>