# Need for focus on microbial species following ice melt and changing freshwater regimes in a Janus Arctic Gateway

Nathalie Joli, Michel Gosselin, Mathieu Ardyna, Marcel Babin, Deo Florence Onda, Jean-Éric Tremblay, Connie Lovejoy

► To cite this version:

Nathalie Joli, Michel Gosselin, Mathieu Ardyna, Marcel Babin, Deo Florence Onda, et al.. Need for focus on microbial species following ice melt and changing freshwater regimes in a Janus Arctic Gateway. Scientific Reports, 2018, 8, pp.9405. 10.1038/s41598-018-27705-6. hal-01829749

HAL Id: hal-01829749

https://hal.sorbonne-universite.fr/hal-01829749

Submitted on 4 Jul 2018

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
SCIENTIFIC REPORTS

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

OPEN

# Need for focus on microbial species following ice melt and changing freshwater regimes in a Janus Arctic Gateway

Nathalie Joli<fnref n="1" />,², Michel Gosselin³, Mathieu Ardyna⁴,⁵, Marcel Babin<fnref n="1" />, Deo Florence Onda<fnref n="1" />,², Jean-Éric Tremblay<fnref n="1" /> & Connie Lovejoy<fnref n="1" />,²

Received: 18 January 2018

Accepted: 8 June 2018

Published online: 20 June 2018

Oceanic gateways are sensitive to climate driven processes. By connecting oceans, they have a global influence on marine biological production and biogeochemical cycles. The furthest north of these gateways is Nares Strait at the top of the North Water between Greenland and Ellesmere Island (Canada). This gateway is globally beneficial, first by supporting high local mammal and bird populations and second with the outflow of phosphate-rich Arctic waters fueling the North Atlantic spring bloom. Both sides of the North Water are hydrologically distinct with counter currents that make this Arctic portal a Janus gateway, after Janus, the Roman god of duality. We examined oceanographic properties and differences in phytoplankton and other protist communities from the eastern and western sides of the North Water (latitude 76.5°N) and found that species differed markedly due to salinity stratification regimes and local hydrography. Typical Arctic communities were associated with south flowing currents along the Canadian side, while potentially noxious *Pseudo-nitzschia* spp. were dominant on the Greenland side and associated with greater surface freshening from ice melt. This susceptibility of the Greenland side to *Pseudo-nitzschia* spp. blooms suggest that monitoring species responses to climate mediated changes is needed.

Oceanic gateways have been critical in mediating past climatic changes and modulate heat and salinity exchange between ocean basins<fnref n="1" />. Nutrients and living organisms are also transported across these gateways and can influence local and global biogeochemical cycles²,³. The Arctic Ocean and surrounding seas are under threat due to rapid climate change⁴ with effects on net primary production and carbon fluxes⁵,⁶. The gateways to Baffin Bay (Nares Strait, Jones Sound and Lancaster Sound) between Canada and Greenland are of special interest as net Arctic outflows. These outflows influence nutrient balances⁷,⁸ and species composition in the North Atlantic⁹. Nares Strait is the deepest and the most northerly of all Arctic gateways and is contiguous with the historically productive North Water (locally referred to as Pikialasorsuaq), where human settlement dates to up to 1000 years before present (BP)¹⁰. The North Water Polynya (open water surrounded by sea ice) usually occurs in winter due to the formation of an ice bridge across Nares Strait¹¹,¹². The area supports large populations of marine mammals and seabirds¹³ partly due to a longer open water season enabling early season access to resources by migrating species, but also aided by estuarine type circulation that results in high productivity¹⁴. A northward flow along the west coast of Greenland (West Greenland Current) transports both surface and deeper warmer saltier nitrate-rich Atlantic Water as far north as Smith Sound at the southern end of Nares Strait (Fig. 1A), before being forced south, driven by cyclonic circulation¹⁵. As well on the Greenland side, non-sea ice sourced freshwater has been detected in the Polar Mixed Layer (PML) that lies between ca. 5 m and the first major halocline in the Arctic water column¹⁶. On the Canadian side of the North Water, Arctic Ocean water flows directly southward¹², transporting

<footnote n="1">¹Département de biologie, Québec Océan and Takuvik Joint International Laboratory (UMI 3376), Université Laval (Canada) - CNRS (France), Université Laval, Québec, QC, G1V 0A6, Canada. ²Institut de biologie intégrative et des systèmes (IBIS), Université Laval, Québec, QC, G1V 0A6, Canada. ³Institut des sciences de la mer de Rimouski, Université du Québec à Rimouski, 310 Allée des Ursulines, Rimouski, QC, G5L 3A1, Canada. ⁴Sorbonne Universités, UPMC Paris 06, INSU-CNRS, Laboratoire d'Océanographie de Villefranche, 181 Chemin du Lazaret, 06230, Villefranche-sur-mer, France. ⁵Department of Earth System Science, Stanford University, Stanford, CA, 94305, USA. Correspondence and requests for materials should be addressed to N.J. (email: joli@biologie.ens.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** North Water sampling sites and temperature-salinity plots at the Canadian (Stn 101) and Greenland (Stn 115) repeated sampling sites. (A) Dots on the map corresponds to the drift sites where samples were collected (see Supplementary Table S1). (B) Potential temperature – salinity (T–S) diagrams from CTD profiles to 100 m. The T–S points are color coded for sampling depth. The 20 m samples are indicated by sequential sampling time (T₁–T₇) from the Canadian (black labels) and Greenland (red labels) sides.</figcaption></figure>

cold phosphate rich water from the Central Arctic and Lincoln Sea, along with first year and multi-year sea-ice into the North Atlantic¹⁷.

Data from satellites indicate that the timing and extent of the spring-summer surface chlorophyll signal in the North Water region is highly variable but tending to occur further North over time¹⁸. While valuable, synoptic views from space lack information on chlorophyll concentrations below the surface¹⁹,²⁰, and are blind to most species changes. Recent increasing northern penetration of Atlantic Water along the West Coast of Greenland and melting of the Greenland Ice Sheet (GrIS) could alter the nature of the Nares Strait gateway at multiple levels. For example, local assemblages of the eukaryotic phytoplankton and other protists likely have differing carbon and nitrogen cycle regimes²¹. In addition, since microbial eukaryotes are the base of the food chain, species changes could influence higher trophic levels¹³. Given sensitivity of the Arctic to climate change and the key role of Northern Baffin Bay in the Arctic-Atlantic system, an understanding of the species variability within this gateway is needed. An earlier microscopy study based on samples collected from a four-month series of transects from April to July 1998 in the North Water, found that microbial eukaryotic communities in both the surface and the upper waters that make up the Polar Mixed Layer (PML), differed regionally²². Since the Canadian and Greenland sides of the North Water region are influenced by different hydrological processes, our working hypothesis was that the microbial eukaryotic communities would differ in August, despite similar light levels at 20 m and the post bloom low nutrient conditions that follow spring and early summer high productivity and nutrient depletion¹⁸,²³. The repeated transects of the 1998 study²⁴ did not resolve short-term variability, in space or time, and the lack of closely spaced replicates, could be considered a shortcoming. Multiple profiles and sample collection at the same location is extremely rare in this region of the Arctic and one profile per station is routine¹⁸. In addition, previous expeditions to the North Water have sampled stations at more or less random times over the 24-hour cycle due to logistic constraints on multidisciplinary missions²⁵. For these reasons, we repeatedly sampled the bottom of the PML in two one-day Lagrangian drifts, with one each of the Canadian and Greenland sides of the North Water (Fig. 1A, Supplementary Table S1). We spaced the sample collection approximately every 4 hours over 24 hours to account for any differences between the two sides that may have been associated with solar or tidal cycles. To investigate microbial eukaryotic community variability, we focused on the V4 region of 18S rRNA using high throughput sequencing (HTS). Given that microbial eukaryotes are sensitive indicators of water masses²⁶, the aim was to capture short term temporal and spatial variability in biological, chemical and physical parameters and link these to species composition and community changes.

## Results

**Oceanographic characteristics.** At the time of sampling in mid-august, the temperature-salinity (TS) profiles showed fresher surface waters and more saline deep waters on the Greenland side compared to the Canadian side (Fig. 1B). On the Canadian side, there was a slight increase in salinity between the surface and the weak sub-surface chlorophyll maximum (SCM) (31.45 to 32.35). This was in contrast to fresher surface waters (salinity of 30.51) and a sharp halocline on the Greenland side with salinity of 33.32 at the SCM (Table 1, Fig. 1B). The surface chlorophyll (Chl *a*) concentrations from satellite data were similar on the two sides and both sides had been ice free since the end of June (Supplementary Fig. S1). Total Chl *a* from samples collected at the surface were also similar on the two sides (Supplementary Fig. S2; ANOVA, p-value = 0.304). At the bottom of the PML starting at 20 m depth, the two sides differed in nearly all respects. On the Canadian side, *in situ* chlorophyll fluorescence (Fig. 2, upper panel) and total Chl *a* concentrations (Supplementary Fig. S2A) changed little down the water
<table><caption>**Table 1.** Environmental data for T₁ to T₇ from 20 m and the surface (Surf) and subsurface chlorophyll maxima (SCM) at T₁ for Stn 101 and at T₂ for Stn 115; Temperature (T) in °C, salinity (Sₚ), dissolved oxygen (DO) in µM, relative fluorescence of colored dissolved organic matter (fCDOM), nitrite (NO₂), nitrate (NO₃), phosphate (PO₄) and silicate (Si) in µM. Ratio of silicate to nitrite + nitrate (Si:N).</caption><thead><tr><th>Stn</th><th>Depth</th><th>Sample</th><th>T</th><th>S<sub>p</sub></th><th>DO</th><th>fCDOM</th><th>NO<sub>2</sub></th><th>NO<sub>3</sub></th><th>PO<sub>4</sub></th><th>Si</th><th>Si:N</th></tr></thead><tbody><tr><td rowspan="10">101</td><td>20</td><td>T<sub>1</sub></td><td>2.65</td><td>32.18</td><td>353.3</td><td>3.13</td><td>0.08</td><td>0.15</td><td>0.56</td><td>0.58</td><td>2.52</td></tr><tr><td>20</td><td>T<sub>2</sub></td><td>3.97</td><td>31.70</td><td>341.6</td><td>2.84</td><td>0.07</td><td>0.12</td><td>0.47</td><td>0.52</td><td>2.71</td></tr><tr><td>20</td><td>T<sub>3</sub></td><td>3.44</td><td>32.11</td><td>350.1</td><td>3.19</td><td>0.09</td><td>0.15</td><td>0.49</td><td>0.93</td><td>3.94</td></tr><tr><td>20</td><td>T<sub>4</sub></td><td>3.65</td><td>32.01</td><td>341.2</td><td>2.99</td><td>0.08</td><td>0.14</td><td>0.48</td><td>0.89</td><td>4.13</td></tr><tr><td>20</td><td>T<sub>5</sub></td><td>3.90</td><td>31.88</td><td>343.9</td><td>2.99</td><td>0.09</td><td>0.14</td><td>0.47</td><td>0.95</td><td>4.32</td></tr><tr><td>20</td><td>T<sub>6</sub></td><td>3.24</td><td>32.20</td><td>332.7</td><td>3.19</td><td>0.09</td><td>0.20</td><td>0.48</td><td>0.97</td><td>3.41</td></tr><tr><td>20</td><td>T<sub>7</sub></td><td>3.55</td><td>32.14</td><td>334.1</td><td>3.20</td><td>0.09</td><td>0.28</td><td>0.53</td><td>0.93</td><td>2.50</td></tr><tr><td>2</td><td>Surf</td><td>4.27</td><td>31.45</td><td>336.7</td><td>2.98</td><td>0.08</td><td>0.21</td><td>0.46</td><td>0.61</td><td>2.09</td></tr><tr><td>38</td><td>SCM</td><td>1.33</td><td>32.35</td><td>352.8</td><td>3.79</td><td>0.10</td><td>0.72</td><td>0.71</td><td>0.33</td><td>0.41</td></tr><tr><td rowspan="10">115</td><td>20</td><td>T<sub>1</sub></td><td>-0.13</td><td>33.03</td><td>414.4</td><td>2.67</td><td>0.11</td><td>0.21</td><td>0.33</td><td>0.66</td><td>2.09</td></tr><tr><td>20</td><td>T<sub>2</sub></td><td>-0.41</td><td>33.11</td><td>410.9</td><td>3.02</td><td>0.10</td><td>0.18</td><td>0.32</td><td>0.52</td><td>1.87</td></tr><tr><td>20</td><td>T<sub>3</sub></td><td>2.60</td><td>32.46</td><td>374.3</td><td>2.41</td><td>0.11</td><td>0.19</td><td>0.32</td><td>1.00</td><td>3.35</td></tr><tr><td>20</td><td>T<sub>4</sub></td><td>1.37</td><td>32.52</td><td>373.8</td><td>2.18</td><td>0.09</td><td>0.26</td><td>0.23</td><td>1.10</td><td>3.08</td></tr><tr><td>20</td><td>T<sub>5</sub></td><td>2.85</td><td>32.63</td><td>367.1</td><td>2.26</td><td>0.10</td><td>0.18</td><td>0.30</td><td>0.31</td><td>1.14</td></tr><tr><td>20</td><td>T<sub>6</sub></td><td>3.90</td><td>32.20</td><td>317.1</td><td>2.19</td><td>0.11</td><td>0.12</td><td>0.26</td><td>0.72</td><td>3.20</td></tr><tr><td>20</td><td>T<sub>7</sub></td><td>2.42</td><td>32.61</td><td>361.3</td><td>2.20</td><td>0.11</td><td>0.15</td><td>0.24</td><td>0.80</td><td>3.02</td></tr><tr><td>2</td><td>Surf</td><td>3.47</td><td>30.51</td><td>332.3</td><td>2.19</td><td>0.11</td><td>0.20</td><td>0.20</td><td>1.30</td><td>4.25</td></tr><tr><td>30</td><td>SCM</td><td>-0.91</td><td>33.32</td><td>333.6</td><td>3.63</td><td>0.27</td><td>3.97</td><td>0.61</td><td>4.26</td><td>1.01</td></tr></tbody></table>

column, with only a very weak SCM (0.69 mg Chl *a* m⁻³) detected at ca. 40 m. On the Greenland side, the SCM at 30 m was pronounced with concentrations reaching 2.56 mg Chl *a* m⁻³. Dissolved oxygen (DO) concentrations were also greater on the Greenland side (Table 1 and Supplementary Table S2) compared to the Canadian side. For nutrients, phosphate concentrations were higher down the water column on the Canadian side. Nitrate concentrations were depleted to at least 20 m on both sides and increased with depth below the PML with higher concentrations on the Greenland side. There was a nitrite maximum detected on the Greenland side at 40 m (Fig. 2).

Our repeated sampling indicated that nitrate and silicate concentrations at 20 m were similar on both sides of the North Water (Supplementary Table S2). However, phosphate concentrations and fluorescent colored dissolved organic matter (fCDOM) showed significant differences between the two sides at 20 m. On the Greenland side, phosphate (avg. 0.29 µM) and fCDOM (avg. 2.42 fluorescence units) were significantly lower compared to the Canadian side (phosphate avg. 0.49 µM and fCDOM avg. 3.08 units) (Table 1, Supplementary Table S2).

**Microbial communities.** For T₁ and T₄, both DNA and RNA templates were sequenced. The relative proportions of taxa inferred from the DNA (rDNA, Supplementary Fig. S3) was similar to proportions of the same taxa from the samples inferred from the RNA (rRNA, see bars for T₁ and T₄ in Fig. 3). A general exception was a higher proportions of dinoflagellates in the rDNA (Supplementary Fig. S3, Supplementary Table S3). Both templates showed the marked difference in the microbial eukaryote communities on the two sides of the North Water (Fig. 3, Supplementary Fig. S3). Unweighted clustering of the rRNA dominants Operational Taxonomic Units (OTUs with >100 read occurrences) also clearly separated communities on the two sides (Supplementary Fig. S4). Environmental variables (Supplementary Table S2) that were significantly correlated with Principle Coordinates Analysis (PCoA) community values were used for the Redundancy Analysis (RDA). The Greenland and Canadian sides clearly separated along the RDA1 axis explaining 20.9% of the total variance (Fig. 4). Higher phosphate, fCDOM and temperature were associated with the Canadian communities, with nitrite, DO and salinity associated with the Greenland communities. Analysis of variance (ANOVA) repeated measures and Kruskal-Wallis tests supported the RDA results (Supplementary Table S2).

Over the ca. 24 hours of 20 m repeated sampling, the Canadian side communities varied substantially, with the T₁ community predominantly the small diatom, *Chaetoceros* aff. *gelidus* (up to 72% of reads at T₁). *Chaetoceros* had largely disappeared 4 hours later; replaced by *Phaeocystis*, which accounted for 43% of the small fraction reads at T₂. Subsequently we detected a mix of small photosynthetic flagellates, which co-occurred with another Hylochaete *Chaetoceros* species; *C. neogracile* (Fig. 3, Supplementary Table S4A). The communities on the Greenland side at 20 m were less variable (Supplementary Table S5), with all samples containing the potentially toxic pennate diatom *Pseudo-nitzschia* (up to 52% of rRNA reads at T₁). The two dominant diatom genera *Chaetoceros* and *Pseudo-nitzschia* were genetically diverse with multiple OTUs found on both the Canadian and Greenland sides (Supplementary Fig. S5). The presence of both genera was confirmed microscopically (Supplementary Fig. S6 plates a,b for *Chaetoceros*, and Fig. S6 plates g,h,i,j for *Pseudo-nitzschia*).

Alveolates represented the second most abundant major group from rRNA on both sides of the North Water. Ciliates made up 28% of the total reads on the Canadian side and 30% of the total reads on the Greenland side (Fig. 3, Supplementary Table S4). Most of the ciliate OTUs were in the order Oligotrichida, with contributions from Choreotrichida, Cyclotrichida (Askenasia), Prorodontida (Urotricha), Haptorida (Monodinium) and other
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Vertical profiles of the upper 60 m of the water column on the Canadian side (Stn 101) and Greenland side (Stn 115) of the North Water. Upper panels) attenuation of photosynthetically available radiation (PAR), nitrate concentration (in µM), and water density (σ<sub>t</sub> in kg m<sup>-3</sup>). Lower panels) phosphate and nitrite concentrations (in µM), and chlorophyll fluorescence (Chl in arbitrary units). The single dotted line for PAR attenuation was from T<sub>2</sub> for Stn 101 and from T<sub>3</sub> for Stn 115. Central solid lines are the average of 7 profiles (time T<sub>1</sub> to T<sub>7</sub>), with exterior dotted lines showing the standard error.</figcaption></figure>

unidentified Litostomatea detected on both sides. *Urotricha* was common in the large fraction on the Canadian side but less on the Greenland side. In contrast, the genus *Monodinium* was more abundant on the Greenland side compared to the Canadian side (Supplementary Table S4). Core dinoflagellates accounted for ca. 30% of total reads on the Canadian side and 21% of total reads on the Greenland side. For both sides, unclassified dinoflagellates accounted for 12% to 18% of reads (Supplementary Table S4). The most common dinoflagellate OTUs on the two sides were *Gymnodiniales*. At the species level *Gyrodinium spirale* was common on the Canadian side and *Torodinium robustum*, which was also identified microscopically (Supplementary Fig. S6 plates c,d) was common on both sides. The majority of dinoflagellate taxa were morphologically diverse gymnodinoids and gyrodinoids between 6–50 µm and not identified further using microscopy. From HTS, *Syndiniales* were predominantly in the *Duboscquellidae* in Group I marine alveolates (MALV I)²⁷ and the still poorly known environmental MALV III, associated with the environmental clone OLI11005²⁸ (GenBank AJ402349).

The relative abundances of smaller flagellated taxa differed between the two sides. For example, several hap-
tophyte OTUs were more common on the Canadian side with *Phaeocystis* accounting for 8.8% of reads in the
small size fraction, and *Chrysochromulina* with 1.2% in the large fraction (Fig. 3, Supplementary Table S4). The
Arctic *Micromonas* CCMP 2099 was relatively more abundant on the Canadian side with 6.1% of the small frac-
tion reads, compared to the Greenland side with 3.8% of the reads (Fig. 3, Supplementary Table S4). A second
chlorophyte, *Pyramimonas* sp. (Supplementary Fig. S6 plates i,j), accounted for nearly 3% of small fraction reads
on the Greenland but only 0.1% on the Canadian side. Two heterotrophic flagellate taxa; Picozoa and group 2
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Barplot of the most abundant species per station per fraction. (A) Stn 101 small fraction. (B) Stn 115 small fraction (C) Stn 101 large fraction. (D) Stn 115 large fraction.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Redundancy analysis from (Canadian) Stn 101 (black dots) and (Greenland) Stn 115 (red dots) based on unweighted UNIFRAC values from the entire microbial eukaryote community within each sample.</figcaption></figure>

uncultivated marine stramenopiles (MAST 2), were common on the Greenland side with each accounting for 1.9% of reads in the small fraction. There were marked differences in the dominant chrysophytes on the two sides (Fig. 3 and Supplementary Table S4) with *Dinobryon balticum* identified from the Canadian side (Supplementary Fig. S6, plates e,f) and *D. belgica* on the Greenland side (Supplementary Fig. S6 plates k,l).

## Discussion

Within the PML the microbial eukaryotic communities from the two sides of the North Water differed taxonom- ically, with *Pseudo-nitzschia* dominated communities persistently found on the Greenland side, and more mixed communities on the Canadian side. The Canadian communities changed several times over the nearly 24 hours of sampling (Fig. 3). The high similarity among samples collected on the Greenland side (Supplementary Table S5) suggested that a single sample could have been used to represent the community in mid-August. In contrast, no single typical community was detected on the Canadian side. The active communities detected using rRNA,
were consistent with the physical oceanography of the respective sides. Salinity determines Arctic Ocean physical oceanographic structure²⁹ and on the Greenland side, Atlantic waters contribute to the higher salinity and sensible heat at depth¹¹, which leads to shoaling of the upper mixed layer³⁰. During our study the fresh surface waters on the Greenland side originated from first year and multi-year sea ice melt and also likely from iceberg and glaciier melt from around Greenland¹⁶,³¹. A high external input of surface fresh water originating from the Greenland Ice Sheet (GrIS) was previously suggested as an explanation for anomalous values for alkalinity and delta O¹⁸ data from the same station on the same cruise¹⁶. The Canadian side reflected the more general freshening of the Arctic Ocean that is due to multi-year sea ice melt³². This overall freshening on the Canadian side could be seen in the weak halocline and lower average water column salinity (Fig. 1B, Fig. 2). The higher phosphate and fCDOM concentrations (Table 1, Table S2), on the Canadian side were also characteristic of Arctic Water³³, in keeping with the Nares Strait Gateway directly contributing excess Arctic Ocean phosphate to the Western Labrador Sea and North Atlantic⁷,⁸.

At the surface, the two sides of Northern Baffin Bay were similar in terms of ice cover and surface bloom phenology (Supplementary Fig. S1), with surface chlorophyll concentrations on both sides within the long-term norm for mid-August¹⁸. However, below the surface, along with the marked differences in physical structure, differences in biological indicators, such as the prominent SCM and overall higher Chl *a* concentrations at depth on the Greenland side, were striking. The formation of the Arctic SCM below the PML depends on the availability of both light and nutrients³⁴ and the high nitrate concentrations below 30 m on the Canadian side, suggest that phytoplankton may have been light limited at that depth. If so, light limitation could be explained by water column instability associated with mixing below the halocline, moving phytoplankton out of the euphotic zone to depths below irradiance levels sufficient for net growth. Alternatively, or in addition, the lack of formation of a prominent SCM on the Canadian side may have been due to the hydrologic complexity of the region, where multiple water-masses with similar densities and different recent histories could converge and interleave³⁵. Such interleaving would be consistent with temperature-salinity (TS) profiles (Fig. 1B) and rapid community changes (Fig. 3) at the single depth sampled. In contrast, there was a prominent SCM on the Greenland side, which was associated with a well stratified water column and presumably continuous nutrient supply to the bottom of the euphotic zone beneath the sharp halocline. Saline Atlantic water below the fresher PML is a characteristic of the Greenland side³⁰, and such ideal conditions for SCM formation³⁶ (Fig. 2, Supplementary Fig. S2) would be expected, suggesting annual formation of an SCM on the Greenland side.

The RDA (Fig. 4) supported the notion that the distinct communities on the two sides were associated with the physical oceanography of the region¹⁷. Higher phosphate and fCDOM concentrations, characteristic of Arctic surface waters²³ were associated with samples collected from the Canadian side. In contrast, Greenland side communities were associated with greater DO concentrations consistent with local upward diffusion of DO from the photosynthetic activity at the much more robust SCM.

The active community from the RNA template (Fig. 3, Supplementary Table S4), documented by microscopy (Supplementary Fig. S6), revealed the presence of Hyalochaete *Chaetoceros* predominantly at T₁ on the Canadian side, when the colder temperatures and the high fCDOM indicated the highest contribution of unmodified Arctic waters³³. These small diatoms were earlier reported (identified as *C. socialis*) as ubiquitous in the summer phytoplankton in the North Water²⁴. The abundant *Chaetoceros* populations suggested that the 20 m sample at T₁ may have originated nearer the surface under high irradiances³⁷. At T₂ on the Canadian side, the community changed, with a high proportion of reads from *Phaeocystis*, which is another common polar species. For example, *Phaeocystis* is reported from the Ross Sea, Antarctica³⁸. Several studies suggest that *Phaeocystis* rapidly responds to increased light and nutrient inputs³⁹,⁴⁰ consistent with changing conditions on the Canadian side of the Nares Strait gateway. Another small flagellate *Micromonas polaris*, which was recently described as an endemic arctic species⁴¹, and almost universally reported from Arctic 18 S rRNA surveys as *Micromonas* CCMP 2099⁴²,⁴³, was retrieved on both sides but with highest relative abundance after T₅ on the Canadian side (ca 14% in the small fraction; Fig. 3). Other typical Arctic flora on the Canadian side included chrysophytes such as *Dinobryon balticum*, which contributes to carbon flux in cold waters⁴⁴. In contrast, the morphospecies *Dinobryon belgica*, which has smaller colonies (Supplementary Fig. S6 plate e versus l) was noted on the Greenland side. The 30-fold difference in the abundance of *Pyramimonas* on the Greenland versus Canadian side may have been facilitated by the low salinity surface waters. For example, Daugbjerg⁴⁵ indicated that some *Pyramimonas* species are able to tolerate a wide range of salinity and *Pyramimonas* was also reported as a dominant nanoflagellate in surface waters with high freshwater input from the Mackenzie River and surface ice melt in the Beaufort Sea⁴³. In addition, *Pyramimonas* is associated with meltwater pools on sea ice¹⁶ and found during the winter-spring transition in the fjord and ice influenced Disko Bay in Greenland south of our study region⁴⁷.

The molecular data (Fig. 3, Supplementary Table S3, Supplementary Table S3, Supplementary Table S4) verified using microscopy (Supplementary Fig. S6) showed that a major difference between the two sides was the overwhelming dominance of *Pseudo-nitzschia* spp. on the Greenland side, accounting up to 52% of reads, but <0.1% of reads on the Canadian side. The presence of *Pseudo-nitzschia* is of concern, since the genus can produce the neurological toxin, domoic acid (DA)⁴⁸. In contrast, more typical summer Arctic communities⁴⁹ were found on Canadian side, with variability consistent with changes in water masses with different recent histories and species³⁵. *Pseudo-nitzschia* has long been present in this region, but high concentrations have not been previously reported. During the 1998 North Water Polynya (NOW) study, Lovejoy *et al.*²² observed a diverse diatom community on the Greenland side in July, including *Pseudo-nitzschia* cf. *seriata*. A *Pseudo-nitzschia* strain (CCMP 2093) now classified as *P. arctica*⁵⁰ was isolated from the North Water (78°35.87′N, 74°29.53′W) in 1998 and other *Pseudo-nitzschia* species have been isolated from Greenland fjords south of our study region²¹. *Pseudo-nitzschia* have wide salinity tolerances, are able to use a variety of nitrogen sources, and are favored under conditions of high dissolved Fe concentrations⁵². The nitrite peak at 40 m on the Greenland side suggests high biological activity and an additional nitrogen source for *Pseudo-nitzschia*⁵³. Overall, we suggest that GrIS inflows could favor
blooms of *Pseudo-nitzschia* spp. and potentially DA production⁴⁸. DA producing *Pseudo-nitzschia* blooms are detrimental to marine ecosystems since DA causes memory impairment and is harmful to marine mammals⁵⁴. To date, DA has not been reported from the North Water, but there have been few or no DA surveys in the region, to our knowledge. Locally isolated *Pseudo-nitzschia* strains from Disko Bay (Greenland) produce DA when fed to Arctic copepodites⁵¹ and *Pseudo-nitzschia* OTUs here clustered with species such as *P. brasiliana* (Supplementary Fig. S5A), which is known to produce DA⁵², suggesting the potential of DA production. Recently, large-scale phytoplankton blooms in the Labrador Sea were reported to coincide with freshwater discharge from glaciers and meltwater from the GrIS, however information on species was lacking²⁵, and given that the bloom was attributed to increased Fe and organic nutrients, it would not be surprising if *Pseudo-nitzschia* was present, and such phenomena warrant investigation.

In sum, although classified as a single eco-region for modeling studies⁵⁶, The Greenland side of the gateway was more saline below PML, but fresher on the surface (Fig. 1B), consistent with recent increases in melt water from the Greenland Ice Sheet (GrIS)⁵⁷. Freshwater from the GrIS, glacial melt and multi-year sea-ice melt enter the North Water via Nares Strait as surface water along the Greenland side¹⁷. Additional GrIS melt from Northeast Greenland⁵⁸ could also contribute to surface freshwater on the Greenland side of the North Water since it flows southward along the east coast of Greenland (via the East Greenland Current) and at the southern tip of Greenland, becomes entrained into the northward coastal flow (West Greenland Current)⁵⁷. Any increase in GrIS meltwater and associated high iron (Fe) and other micronutrient concentrations⁵⁹ has implications for the biogeochemistry of Baffin Bay⁵⁵ and would create conditions that could favor blooms of *Pseudo-nitzschia*. This in turn, suggests a threat to the historic role of the North Water in supporting marine bird and mammal populations including endangered Greenland populations of Bowhead whales, which migrate to this latitude in summer and graze on copepods⁶⁰. Finally, given the increasing cyclonic circulation in Northern Baffin Bay, low-phosphate waters and potential toxic species might flow south to the North Atlantic along the Baffin Bay meridional ridge¹⁷, making this a true Janus gateway. This work highlights the need to consider drivers that select for species and the biological diversity of gateways, if we are to fully understand impacts of climate change in the Arctic and elsewhere.

## Materials and Methods

**Satellite-derived products.** Satellite-derived level-3 data sets of Chl *a* concentration (mg m⁻³) and particulate back-scattering coefficients at 443 nm (bₚₚ, m⁻¹) (Supplementary Fig. S1) were obtained from the European Space Agency's GlobColour project (http://www.globcolour.info). Eight-day composite Chl *a* concentration and bₚₚ coefficients were determined, respectively using standard Case 1 water algorithms⁶¹ and the semi-analytical Garver Siegel Maritorena (GSM) merging algorithm⁶². The sea-ice concentration was derived from the Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E) sensor, and made available by National Snow and Ice Data Center (NSIDC; https://nsidc.org). Actual sea ice conditions northward and just prior to sampling were further verified by inspecting a range of additional satellite products (https://worldview.earthdata.nasa.gov/).

**Field sampling.** Oceanographic data was collected aboard the CCGS *Amundsen* using a rosette system equipped with a conductivity, temperature, depth (CTD) profiler (Sea-Bird SBE-911 CTD), relative nitrate (*In-Situ* Ultraviolet Spectrometer, ISUS, Satlantic), oxygen (Seabird SBE-43), chlorophyll fluorescence (Seapoint), fluorescent colored dissolved organic matter (*f*CDOM; Wetlabs ECO) and photosynthetically available radiation (PAR, 400–700 nm; Biospherical Instruments QDP2300) sensors. The oxygen sensor was calibrated onboard against Winkler titrations³⁶.

The ship followed a buoy with a drogue suspended to 20 m, with the goal of following the bottom of the PML and sampled approximately every 4 hours beginning at solar time 05:15 for Stn 101 (Canadian side) and 06:00 for Stn 115 (Greenland side), ending after 23 hours at and 22.5 hours, respectively (Supplementary Table S1). At Stn 101, the drift started ($T_1$) at lat 76°23.242'N, long 077°23.412' W, and ended ($T_7$) at lat 76°17.372'N, long 077°46.264'W. The Stn 115 drift started ($T_1$) at lat 76°20.297'N, long 071°11.514'W and ended at ($T_7$) at 76°29.858'N, long 071°26.615'W (Fig. 1A and Supplementary Table S1). Water samples were collected on the upcast from 12–L Niskin-type bottles mounted on the rosette. Nutrient samples were collected every 10 m from the surface to 60 m. Nutrients were analyzed on board within 4 h after collection as using a Bran-Luebbe 3 auto-analyzer⁶³. Optical depths were derived from vertical PAR profiles taken from the shadow-free side of the ship using a PNF-300 radiometer (Biospherical Instruments). At time $T_1$ for Stn 101 and time $T_3$ for Stn 115, water samples for duplicate Chl *a* were collected at seven optical depths (100, 50, 30, 15, 5, 1, and 0.2% of surface irradiance), the SCM determined from the downward cast of the CTD, and at 80 m and 100 m. Chl *a* samples were also collected and analyzed from 20 m along with nucleic acids (see below) during the drift. Samples for total Chl *a* were filtered onto 25 mm Whatman GF/F filters (TChl *a*; >0.7 µm) and onto 5 µm pore size polycarbonate (PC) membrane filters (Nuclepore™) to estimate the large fraction (L-Chl *a*; >5 µm) and analyzed on board as in Parsons *et al*.⁶⁴. Concentrations in the smaller fraction (S-Chl *a*; 0.7–5 µm) were from subtraction. For nucleic acids, five independent samples from 20 m were collected from each the Canadian and Greenland sides (Supplementary Table S1). Six L of sample water prefiltered through a 50 µm nylon mesh, and then sequentially filterthrough a 47–mm diameter 3 µm pore size PC membrane filter, and a 0.2 µm Sterivex™ Unit (Millipore). Material on the filters was preserved in RNALater™ (ThermoFisher) then frozen at –80 °C. The different size fractionation protocols for Chl *a* and nucleic acids was due to logistic constraints aboard the ship, and small and large designations should be considered as indicative only.

Whole water samples for fluorescence microscopy were collected at 20 m at $T_1$ for Stn 101 and $T_2$ for Stn 115. Single aliquots of 48 mL were preserved and filtered onto black 0.8 µm PC filters, stained with 4,6-diamidino-2-phenylindole (DAPI) following⁶⁵. The slides were inspected at 1000X magnification using an
Olympus IX71 microscope equipped with UV and blue excitation filter blocks. Images were captured with the integrated QImaging Retiga 2000R CCD Camera and QCapture software version 2.9.11.

**Laboratory procedures.** DNA and RNA were extracted from the same filters using the All-Prep DNA/RNA Minikit (Qiagen). DNA was used to verify the overall community including dead and dormant cells at two time points 12 hours apart on either side (T₁ and T₂). RNA was converted to cDNA using the High Capacity Reverse Transcription Kit (Applied Biosystems). The V4 region of 18 S rRNA gene and rRNA (from cDNA) was amplified using the eukaryote specific primers E572 (forward) and E1009 (reverse), see Comeau *et al.*⁶⁶, coupled with a MiSeq® specific linking primer. To decrease potential PCR bias, 1, 5 and 10-fold diluted template was used for PCRs for each sample. The PCR products of the three dilutions were pooled together and purified using the Axygen® PCR cleanup kit (Axygen) and then quantified spectrophotometrically with the Nanodrop 1000™ (ThermoFisher Scientific). Unique pairs of barcodes (tags) were added to the sample amplicons using the TruSeq® and Nextera® (both Illumina) barcode sets in a nested PCR as described in Comeau *et al.*⁶⁷. Equimolar concentrations of the sample barcoded amplicons were sequenced on the Illumina MiSeq at the Plate-forme d'Analyses Génomiques (IBIS, Université Laval, Québec, QC, Canada). All reads are deposited in NCBI GenBank Sequence Read Archive (SRA) under the BioProject number PRJNA383398 (GenBank: SRX2745611 to SRX2745642).

**Data and statistical analysis.** Paired end reads were processed with UPARSE⁶⁸, with microbial eukaryote operational taxonomic units OTUs retained following Wu *et al.*⁶⁹. Based on a 99% similarity the reads clustered into 2069 OTUs. We then used an arctic specific reference database for taxonomic assignment of OTUs⁷⁰. The resulting OTU classifications are available at https://zenodo.org/record/1205261#.WrOAMJPwb-Y and their associated reads at https://zenodo.org/record/1205255#.WrN_RJPwb-Y.

Data sets were further processed in QIIME⁷¹. Rarefied datasets (11,500 sequences per sample) were used as input for the UniFrac unweighted distances⁷², with beta diversity and principal coordinates (PCoA) as output. Environmental variables that were significantly correlated with PCoA community values were used for the redundancy analysis (RDA)⁷³ computed using the *rda* function in the Rstudio package ‘vegan’ (v2.4-1)⁷⁴. Parameters that best explained variability in the PCoA were selected using the ordistep function in the ‘vegan’ package to build an optimal model (highest adjusted coefficient determination). One-way ANOVA (parametric) and Kruskal-Wallis (non-parametric) tests were carried out in PAST v3.0⁷⁵ to test for differences between the two sides of the North Water. ANOVA repeated measures tests were carried out using PAST v3.0⁷⁵. Statistical tests of physical and nutrient data and profile plots were carried out in Sigma Plot (v11).

With the aim to focus on more abundant organisms, we selected OTUs that had a minimum of 100 occurrences in at least one rRNA sample from both sides of the North Water Clusters and heatmaps were constructed using the heatmap.2 function in gplots package in R. Rarefied rRNA read counts were first log transformed and then clustered by site (Y axis) and by OTUs (X axis) using Euclidean distances. Dendograms were generated using hierarchical clustering UPGMA.

To explore the diversity of the two dominant diatom genera in our samples, reference sequences related to *Pseudo-nitzschia* were retrieved from Percopo *et al.*⁵⁰ and those related to *Chaetoceros* were retrieved from Chamnansin *et al.*⁷⁶. The V4 18 S rRNA sequences were aligned separately with MAFFT (v7)⁷⁷.

**References**

1. Karas, C. *et al.* Pliocene oceanic seaways and global climate. *Sci. Reports* 7, 39842, https://doi.org/10.1038/srep39842 (2017).

2. Takahashi, K. The Bering and Okhotsk Seas: modern and past paleoceanographic changes and gateway impact. *J. Asian Earth Sci.* 16, 49–58, https://doi.org/10.1016/s0743-9547(97)00048-2 (1998).

3. Wassmann, P. *et al.* The contiguous domains of Arctic Ocean advection: Trails of life and death. *Prog.Oceanogr.* 139, 42–65, https://doi.org/10.1016/j.pocean.2015.06.011 (2015).

4. Notz, D. & Stroeve, J. Observed Arctic sea-ice loss directly follows anthropogenic CO₂ emission. *Science* 354, 747–750, https://doi.org/10.1126/science.aag2345 (2016).

5. Tremblay, J. É. *et al.* Global and regional drivers of nutrient supply, primary production and CO₂ drawdown in the changing Arctic Ocean. *Prog.Oceanogr.* 139, 171–196, https://doi.org/10.1016/j.pocean.2015.08.009 (2015).

6. Arrigo, K. R., van Dijken, G. & Pabi, S. Impact of a shrinking Arctic ice cover on marine primary production. *Geophys. Res. Lett.* 35, https://doi.org/10.1029/2008gl035028 (2008).

7. Yamamoto-Kawai, M., Carmack, E. & McLaughlin, F. Nitrogen balance and Arctic throughflow. *Nature* 443, 43–43, https://doi.org/10.1038/443043a (2006).

8. Torres-Valdes, S. *et al.* Export of nutrients from the Arctic Ocean. *J. Geophys. Res.-Oceans* 118, 1625–1644, https://doi.org/10.1002/jgrc.20063 (2013).

9. Luddington, I. A., Lovejoy, C. & Kaczmarska, I. Species-rich meta-communities of the diatom order Thalassiosirales in the Arctic and northern Atlantic Ocean. *J. Plank. Res.* 38, 781–797, https://doi.org/10.1093/plankt/fbw030 (2016).

10. Morrison, D. Radiocarbon Dating Thule Culture. *Arctic Anthropol.* 26, 48–77 (1989).

11. Båcle, J., Carmack, E. C. & Ingram, R. G. Water column structure and circulation under the North Water during spring transition: April-July 1998. *Deep-Sea Res. Pt. II* 49, 4907–4925, https://doi.org/10.1016/s0967-0645(02)00170-4 (2002).

12. Preusser, A., Heinemann, G., Willmes, S. & Paul, S. Multi-decadal variability of polynya characteristics and ice production in the North Water Polynya by means of passive microwave and thermal infrared satellite imagery. *Remote Sens.* 7, 15844–15867, https://doi.org/10.3390/rs71215807 (2015).

13. Heide-Jorgensen, M. P. *et al.* The significance of the North Water Polynya to arctic top predators. *Ambio* 42, 596–610, https://doi.org/10.1007/s13280-012-0357-3 (2013).

14. Beszczynska-Moller, A., Woodgate, R. A., Lee, C., Melling, H. & Karcher, M. A Synthesis of exchanges through the main oceanic gateways to the Arctic Ocean. *Oceanogr.* 24, 82–99, https://doi.org/10.5670/oceanog.2011.59 (2011).

15. Dumont, D., Gratton, Y. & Arbetter, T. E. Modeling wind-driven circulation and landfast ice-edge processes during polynya events in Northern Baffin Bay. *J. Phys. Oceanogr.* 40, 1356–1372, https://doi.org/10.1175/2010JPO4292.1 (2010).

16. Burgers, T. M. *et al.* Surface water pCO₂ variations and sea-air CO₂ fluxes during summer in the Eastern Canadian Arctic. *J. Geophys. Res.-Oceans* 122, 9663–9678, https://doi.org/10.1002/2017jc013250 (2017).
17. Munchow, A., Falkner, K. K. & Melling, H. Baffin Island and West Greenland Current systems in Northern Baffin Bay. *Prog. Oceanogr.* **132**, 305–317, https://doi.org/10.1016/j.pocean.2014.04.001 (2015).

18. Marchese, C. et al. Changes in phytoplankton bloom phenology over the North Water (NOW) polynya: a response to changing environmental conditions. *Polar Biol.* **40**, 1721–1737, https://doi.org/10.1007/s00300-017-2095-2 (2017).

19. Ardyna, M. et al. Parameterization of vertical chlorophyll *a* in the Arctic Ocean: impact of the subsurface chlorophyll maximum on regional, seasonal, and annual primary production estimates. *Biogeosciences*. **10**, 4383–4404, https://doi.org/10.5194/bg-10-4383-2013 (2013).

20. Babin, M. et al. Estimation of primary production in the Arctic Ocean using ocean colour remote sensing and coupled physical-biological models: strengths, limitations and how they compare. *Prog. Oceanogr.* **139**, 197–220, https://doi.org/10.1016/j.pocean.2015.08.008 (2015).

21. Ardyna, M. et al. Shelf-basin gradients shape ecological phytoplankton niches and community composition in the coastal Arctic Ocean (Beaufort Sea). *Limnol. Oceanogr.* **62**, 2113–2132, https://doi.org/10.1002/lno.10554 (2017).

22. Lovejoy, C., Legendre, L., Martineau, M. J., Bâcle, J. & von Quillfeldt, C. H. Distribution of phytoplankton and other protists in the North Water. *Deep-Sea Res. Pt. II* **49**, 5027–5047, https://doi.org/10.1016/s0967-0645(02)00176-5 (2002).

23. Tremblay, J. É., Gratton, Y., Carmack, E. C., Payne, C. D. & Price, N. M. Impact of the large-scale Arctic circulation and the North Water Polynya on nutrient inventories in Baffin Bay. *J. Geophys. Res.-Oceans* **107**, 3112, https://doi.org/10.1029/2000jc000595 (2002).

24. Booth, B. C. et al. Dynamics of *Chaetoceros socialis* blooms in the North Water. *Deep-Sea Res. Pt. II* **49**, 5003–5025, https://doi.org/10.1016/s0967-0645(02)00175-3 (2002).

25. Pedneault, E., Galand, P. E., Polvin, M., Tremblay, J. E. & Lovejoy, C. Archaeal *amoA* and *ureC* genes and their transcriptional activity in the Arctic Ocean. *Sci. Rep* **4**, 4661, https://doi.org/10.1038/srep04661 (2014).

26. Monier, A. et al. Upper Arctic Ocean water masses harbor distinct communities of heterotrophic flagellates. *Biogeosciences* **10**, 4273–4286, https://doi.org/10.5194/bg-10-4273-2013 (2013).

27. Harada, A., Ohtsuka, S. & Horiguchi, T. Species of the parasitic genus *Duboscquella* are members of the enigmatic Marine Alveolate Group I. *Protist* **158**, 337–347, https://doi.org/10.1016/j.protis.2007.03.005 (2007).

28. Guillou, L. et al. Widespread occurrence and genetic diversity of marine parasitoids belonging to Syndiniales (Alveolata). *Environ. Microbiol.* **10**, 3349–3365, https://doi.org/10.1111/j.1462-2920.2008.01731.x (2008).

29. Carmack, E. C. et al. Freshwater and its role in the arctic marine system: sources, disposition, storage, export, and physical and biogeochemical consequences in the Arctic and Global Oceans. *J. Geophys. Res.-Biogeosci.* **121**, 675–717, https://doi.org/10.1002/2015jg003140 (2016).

30. Mei, Z. P. et al. Physical control of spring-summer phytoplankton dynamics in the North Water, April-July 1998. *Deep-Sea Res. Pt. II* **49**, 4959–4982, https://doi.org/10.1016/s0967-0645(02)00173-x (2002).

31. Gillard, L. C., Hu, X. M., Myers, P. G. & Bamber, J. L. Meltwater pathways from marine terminating glaciers of the Greenland ice sheet. *Geophys. Res. Lett* **43**, 10873–10882, https://doi.org/10.1002/2016gl070969 (2016).

32. Munchow, A. Volume and freshwater flux observations from nares strait to the west of Greenland at daily time scales from 2003 to 2009. *J. Phys. Oceanogr.* **46**, 141–157, https://doi.org/10.1175/jpo-d-15-0093.1 (2016).

33. Goncalves-Araujo, R. et al. Using fluorescent dissolved organic matter to trace and distinguish the origin of Arctic surface waters. *Sci. Rep* **6**, 12, https://doi.org/10.1038/srep33978 (2016).

34. Tremblay, J. É. et al. Vertical stability and the annual dynamics of nutrients and chlorophyll fluorescence in the coastal, southeast Beaufort Sea. *J. Geophys. Res.-Oceans* **113**, https://doi.org/10.1029/2007jc004547 (2008).

35. Lovejoy, C., Carmack, E. C., Legendre, L. & Price, N. M. Water column interleaving: A new physical mechanism determining protist communities and bacterial states. *Limnol. Oceanogr.* **47**, 1819–1831 (2012).

36. Martin, J. et al. Prevalence, structure and properties of subsurface chlorophyll maxima in Canadian Arctic waters. *Mar. Ecol. Prog. Ser.* **412**, 69–84, https://doi.org/10.3354/meps08666 (2010).

37. Blais, M. et al. Contrasting interannual changes in phytoplankton productivity and community structure in the coastal Canadian Arctic Ocean. *Limnol. Oceanogr.* **62**, 2480–2497, https://doi.org/10.1002/lno.10581 (2017).

38. Arrigo, K. R., Weiss, A. M. & Smith, W. O. Physical forcing of phytoplankton dynamics in the southwestern Ross Sea. *J. Geophys. Res.-Oceans* **103**, 1007–1021, https://doi.org/10.1029/97jc02326 (1998).

39. Engel, A. et al. Inter-annual variability of transparent exopolymer particles in the Arctic Ocean reveals high sensitivity to ecosystem changes. *Sci. Rep* **7**, 4129, https://doi.org/10.1038/s41598-017-04106-9 (2017).

40. Pavlov, A. K. et al. Altered inherent optical properties and estimates of the underwater light field during an Arctic under-ice bloom of *Phaeocystis poucheti*. *J. Geophys. Res.-Oceans* **122**, 4939–4961, https://doi.org/10.1002/2016jc012471 (2017).

41. Simon, N. et al. Revision of the genus *Micromonas* Manton et Parke (Chlorophyta, Mamiellophyceae), of the type species *M. pusilla* (Butcher) Manton & Parke and of the species *M. commoda* van Baren, Bachy and Worden and description of two new species based on the genetic and phenotypic characterization of cultured isolates. *Protist* **168**, 612–635, https://doi.org/10.1016/j.protis.2017.09.002 (2017).

42. Terrado, R. et al. Protist community composition during spring in an Arctic flaw lead polynya. *Polar Biol.* **34**, 1901–1914, https://doi.org/10.1007/s00300-011-1039-5 (2011).

43. Balzano, S., Marie, D., Gourvil, P. & Vaulot, D. Composition of the summer photosynthetic pico and nanoplankton communities in the Beaufort Sea assessed by T-RFLP and sequences of the 18S rRNA gene from flow cytometry sorted samples. *ISME J.* **6**, 1480–1498, https://doi.org/10.1038/ismej.2011.213 (2012).

44. Olli, K. et al. Seasonal variation in vertical flux of biogenic matter in the marginal ice zone and the central Barents Sea. *J. Mar. Syst.* **38**, 189–204, https://doi.org/10.1016/s0924-7963(02)00177-x (2002).

45. Daugbjerg, N., Hansen, G., Larsen, J. & Moestrup, O. Phylogeny of some of the major genera of dinoflagellates based on ultrastructure and partial LSU rDNA sequence data, including the erection of three new genera of unarmoured dinoflagellates. *Phycologia* **39**, 302–317, https://doi.org/10.2216/i0031-8884-39-4-302.1 (2000).

46. Gradinger, R. Occurrence of an algal bloom under Arctic pack. *Mar. Ecol. Prog. Ser.* **131**, 301–305, https://doi.org/10.3354/meps131301 (1996).

47. Haroardottir, S., Lundholm, N., Moestrup, O. & Nielsen, T. G. Description of *Pyramimonas* *diskoicola* sp. nov. and the importance of the flagellate *Pyramimonas* (Prasinophyceae) in Greenland sea ice during the winter-spring transition. *Polar Biol.* **37**, 1479–1494, https://doi.org/10.1007/s00300-014-1538-2 (2014).

48. McKibben, S. M. et al. Climatic regulation of the neurotoxin domoic acid. *Proc., Natl. Acad. Sci. USA* **114**, 239–244, https://doi.org/10.1073/pnas.1606798114 (2017).

49. Lovejoy, C. et al. Plankton- chapter 3.2. In *State of Arctic Marine Biodiversity Report* (by Conservation of Arctic Flora and Fauna) CAFF, Akureyri, Iceland, (2017).

50. Percopo, I. et al. *Pseudo-nitzschia arctica* sp nov., a new cold-water cryptic *Pseudo-nitzschia* species within the *P. pseudodelicatissima* complex. *J. Phycol.* **52**, 184–199, https://doi.org/10.1111/jpy.12395 (2016).

51. Haroardottir, S. et al. Dangerous relations in the Arctic marine food web: Interactions between toxin producing *Pseudo-nitzschia* diatoms and *Calanus copepodites*. *Mar. Drugs* **13**, 3809–3835, https://doi.org/10.3390/md13063809 (2015).

52. Trainer, V. L. et al. *Pseudo-nitzschia* physiological ecology, phylogeny, toxicity, monitoring and impacts on ecosystem health. *Harmful Algae* **14**, 271–300, https://doi.org/10.1016/j.hal.2011.10.025 (2012).
53. Auro, M. E. & Cochlan, W. P. Nitrogen utilization and toxin production by two diatoms of the *Pseudo-nitzschia pseudodelicatissima* complex: *P. cuspidata* and *P. fryxelliana*. *J. Phycol.* **49**, 156–169, https://doi.org/10.1111/jpy.12033 (2013).

54. Lefebvre, K. A. et al. Prevalence of algal toxins in Alaskan marine mammals foraging in a changing arctic and subarctic environment. *Harmful Algae* **55**, 13–24, https://doi.org/10.1016/j.hal.2016.01.007 (2016).

55. Arrigo, K. R. et al. Melting glaciers stimulate large summer phytoplankton blooms in southwest Greenland waters. *Geophys. Res. Lett.* **44**, 6278–6285, https://doi.org/10.1002/2017GL073583 (2017).

56. Lawrence, J., Popova, E., Yool, A. & Srokosz, M. On the vertical phytoplankton response to an ice-free Arctic Ocean. *J. Geophys. Res.-Oceans* **120**, 8571–8582, https://doi.org/10.1002/2015jc011180 (2015).

57. de la Guardia, L. C., Hu, X. M. & Myers, P. G. Potential positive feedback between Greenland Ice Sheet melt and Baffin Bay heat content on the West Greenland Shelf. *Geophys. Res. Lett.* **42**, 4922–4930, https://doi.org/10.1002/2015gl064626 (2015).

58. Sejr, M. K. et al. Evidence of local and regional freshening of Northeast Greenland coastal waters. *Sci. Reports* **7**, 13183, https://doi.org/10.1038/s41598-017-10610-9 (2017).

59. Cameron, K. A. et al. Meltwater export of prokaryotic cells from the Greenland Ice Sheet. *Environ. Microbiol.* **19**, 524–534, https://doi.org/10.1111/1462-2920.13483 (2017).

60. Nielsen, N. H., Laidre, K., Larsen, R. S. & Heide-Jorgensen, M. P. Identification of potential foraging areas for Bowhead whales in Baffin Bay and adjacent waters. *Arctic* **68**, 169–179 (2015).

61. O'Reilly, J. E. et al. Ocean color chlorophyll algorithms for SeaWiFS. *J. Geophys. Res.-Oceans* **103**, 24937–24953, https://doi.org/10.1029/98jc02160 (1998).

62. Maritorena, S. & Siegel, D. A. Consistent merging of satellite ocean color data sets using a bio-optical model. *Remote Sens. Environ.* **94**, 429–440, https://doi.org/10.1016/j.rse.2004.08.014 (2005).

63. Grasshoff, K., Kremling, K. & Ehrhardt, M. *Methods of seawater analysis*, 3rd edition. (Wiley-VCH, 1999).

64. Parsons, T. R., Maita, Y. & Lalli, C. M. *A manual of chemical and biological methods for seawater analysis*. Pergamon Press (1984).

65. Vaque, D. et al. Seasonal changes in planktonic bacterivory rates under the ice-covered coastal Arctic Ocean. *Limnol. Oceanogr.* **53**, 2427–2438, https://doi.org/10.4319/lo.2008.53.6.2427 (2008).

66. Comeau, A. M., Li, W. K. W., Tremblay, J. É., Carmack, E. C. & Lovejoy, C. Arctic Ocean microbial community structure before and after the 2007 record sea ice minimum. *PLoS One* **6**, e27492, https://doi.org/10.1371/journal.pone.0027492 (2011).

67. Comeau, A. M., Vincent, W. F., Bernier, L. & Lovejoy, C. Novel chytrid lineages dominate fungal sequences in diverse marine and freshwater habitats. *Sci. Reports* **6**, e30120, https://doi.org/10.1038/srep30120 (2016).

68. Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. *Nat. Methods* **10**, 996–998, https://doi.org/10.1038/nmeth.2604 (2013).

69. Wu, W., Logares, R., Huang, B. & Hsieh, C. H. Abundant and rare picoeukaryotic sub-communities present contrasting patterns in the epipelagic waters of marginal seas in the northwestern Pacific Ocean. *Environ. Microbiol.* **19**, 287–300 (2017).

70. Lovejoy, C., Comeau, A. & Thaler, M. Curated reference database of SSU rRNA for northern marine and freshwater communities of Archaea, Bacteria and microbial eukaryotes, v. 1.1 NordicanaD23, https://doi.org/10.5885/45409XD-79A199B76BCC4110 (2002–2008).

71. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. *Nat. Methods* **7**, 335–336, https://doi.org/10.1038/nmeth.f303 (2010).

72. Lozupone, C. & Knight, R. UniFrac: A new phylogenetic method for comparing microbial communities. *Appl. Environ. Microbiol.* **71**, 8228–8235, https://doi.org/10.1128/aem.71.12.8228-8235.2005 (2005).

73. Legendre, P. & Gallagher, E. D. Ecologically meaningful transformations for ordination of species data. *Oecologia* **129**, 271–280, https://doi.org/10.1007/s004420100716 (2001).

74. Oksanen, J. et al. *Vegan: Community Ecology Package*. *R. package version 2.0-10 CRAN* (2013).

75. Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST: PAleontological STatistics software packages for education and data analysis. *Palaeontologia Electronica* **4**, 1–9 (2001).

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

77. Katoh, K. & Frith, M. C. Adding unaligned sequences into an existing alignment using MAFFT and LAST. *Bioinformatics* **28**, 3144–3146, https://doi.org/10.1093/bioinformatics/bts578 (2012).

## Acknowledgements

This study is a contribution to ArcticNet, a Network of Centers of Excellence (Canada). NJ received scholarships from Université Laval and funds from the Canadian Excellence Research Chair (CERC) "Remote Sensing of Canada's New Arctic Frontier" to MB. Research support came from a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery and Northern Supplement grants to CL and MG and the Fonds de recherche du Québec Nature et Technologies (FRQNT) supporting Québec-Océan. We thank D. Colatriano for help during sample collection; P. Guillot for processing CTD data, J. Gagnon and G. Deslongchamps for nutrient analysis and M. Blais for chlorophyll data. L. Benestan and R. Logares provided advice on statistics and bioinformatics. We also acknowledge sequencing support from the Plate-forme d'Analyses Génomiques (IBIS, Université Laval). Data collected by the CCGS Amundsen was made available by the Amundsen Science program, supported by the Canada Foundation for Innovation (CFI) and NSERC.

## Author Contributions

N.J. and C.L. designed the study, N.J. carried out the field and lab work, N.J., D.F.O. and C.L. analyzed the data. M.B., J.E.T., M.G. and M.A. contributed data. N.J. and C.L. wrote the manuscript, all authors commented on and contributed to the various versions of the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-27705-6.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018