# Phylogeny and lipid profiles of snow-algae isolated from Norwegian red-snow microbiomes

Hirono Suzuki, Alexandre Détain, Youngjin Park, Kiron Viswanath, René H Wijffels, Nathalie Leborgne-Castel, Lenka Procházková, Chris J Hulatt

► To cite this version:

Hirono Suzuki, Alexandre Détain, Youngjin Park, Kiron Viswanath, René H Wijffels, et al.. Phylogeny and lipid profiles of snow-algae isolated from Norwegian red-snow microbiomes. FEMS Microbiology Ecology, 2023, 99 (6), pp.art. /fiad057. 10.1093/femsec/fiad057. hal-04163713

HAL Id: hal-04163713

https://hal.inrae.fr/hal-04163713

Submitted on 17 Jul 2023

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
# Phylogeny and lipid profiles of snow-algae isolated from Norwegian red-snow microbiomes

Hirono Suzuki<sup>1,1</sup>, Alexandre Détain<sup>1,1</sup>, Youngjin Park<sup>1</sup>, Kiron Viswanath<sup>1</sup>, René H. Wijffels<sup>1,2</sup>, Nathalie Leborgne-Castel<sup>1,3</sup>, Lenka Procházková<sup>1,4,5</sup>, Chris J. Hulatt<sup>1,6</sup>,

## Abstract

Snow algae blooms often form green or red coloured patches in melting alpine and polar snowfields worldwide, yet little is known about their biology, biogeography, and species diversity. We investigated eight isolates collected from red snow in northern Norway, using a combination of morphology, 18S rRNA gene and internal transcribed spacer 2 (ITS2) genetic markers. Phylogenetic and ITS2 rRNA secondary structure analyses assigned six isolates to the species *Raphidonema nivale*, *Deuterostichococcus epilithicus*, *Chloromonas reticulata*, and *Xanthonema bristolianum*. Two novel isolates belonging to the family Stichococcaceae (ARK-S05-19) and the genus *Chloromonas* (ARK-S08-19) were identified as potentially new species. In laboratory cultivation, differences in the growth rate and fatty acid profiles were observed between the strains. Chlorophyta were characterized by abundant C18:3n-3 fatty-acids with increases in C18:1n-9 in the stationary phase, whilst *Xanthonema* (Ochrophyta) was characterized by a large proportion of C20:5n-3, with increases in C16:1n-7 in the stationary phase. In a further experiment, lipid droplet formation was studied in *C. reticulata* at the single-cell level using imaging flow cytometry. Our study establishes new cultures of snow algae, reveals novel data on their biodiversity and biogeography, and provides an initial characterization of physiological traits that shape natural communities and their ecophysiological properties.

**Keywords:** 18S rRNA, fatty acids, imaging flow cytometry, ITS2 rRNA, microalgae, phylogeny

## Introduction

Microalgae inhabiting terrestrial snow environments are key primary producers in extreme ecosystems and are important for their biodiversity as well as interactions with climate warming (Gray et al. 2020, Krug et al. 2020). Snow algae blooms during spring and summer may form green, orange, or red coloured patches covering small (cm<fnref n="2" />) to large (~km<fnref n="2" />) areas that can be observed in melting polar and alpine snowfields worldwide (Hoham and Remias 2020). Recently, snow algae have attracted a lot of attention due to their contribution to reducing the albedo of snow and ice surfaces, consequently accelerating snowmelt and causing positive feedbacks with a warming climate (Lutz et al. 2016, Cook et al. 2020). Despite their general macroscopic appearance, studies have revealed a substantial amount of global species diversity, much of which has not been characterized.

Red-snow communities are typically dominated by the Chlorophytean algae *Chloromonas*, *Chlainomonas*, and the newly described genus *Sanguina* (Procházková et al. 2019a). Amongst these communities, *Raphidonema* (Trebouxiophyceae) are also commonly observed (Lutz et al. 2016, Segawa et al. 2018, Engstrom et al. 2020). To study the biodiversity and biogeography of snow algae, metabarcoding approaches using common molecular markers have been applied. Lutz et al. (2016) used partial 18S ribosomal RNA (rRNA) gene-based analysis of samples from differ-

ent regions of the Arctic and their data suggested that many groups of snow algae are cosmopolitan. Studies based on faster-evolving internal transcribed spacer 2 (ITS2) rDNA barcodes have also shown that polar red snow communities are typically dominated by a few cosmopolitan species, but also identified additional diverse species that are endemic to particular regions (Segawa et al. 2018). Although amplicon sequencing gives us insight into community structures, accurately distinguishing between closely related species or genera and evaluating biogeographical patterns of snow algae dispersal remain challenging, primarily due to lack of sample site coverage and limited phenotype/genotype reference data for high-throughput sequence analysis (Lutz et al. 2019). Alternatively, isolation and cultivation techniques allow us to characterize cell morphological and physiological properties together with sequence identity at high resolution. In the past, identification of algal species was often limited to morphological observation, but the combined use of molecular markers (18S, ITS2, and *rbcL*) as well as light and electron microscopy have greatly improved our understanding of cryospheric microflora by re-examining isolates and describing novel lineages and species (Hoham and Remias 2020). However, only a fraction of snow algae species, which are often in cyst form in the field, have been successfully cultured (Procházková et al. 2018a).

<footnote n="1"><sup>1</sup>Faculty of Biosciences and Aquaculture, Nord University, 8049 Bodø, Norway</footnote>

<footnote n="2"><sup>2</sup>Bioprocess Engineering, AlgaePARC, Wageningen University, PO Box 16 Wageningen, 6700, AA, The Netherlands</footnote>

<footnote n="3"><sup>3</sup>Agroécologie, Institut Agro Dijon, CNRS, INRAE, Univ. Bourgogne Franche-Comté, F-21000 Dijon, France</footnote>

<footnote n="4"><sup>4</sup>Department of Ecology, Faculty of Science, Charles University, Viničná 7, 12844 Prague, Czech Republic</footnote>

<footnote n="5"><sup>5</sup>Centre for Phycology, Institute of Botany of the Czech Academy of Sciences, Dukelská 135, 37982 Třeboň, Czech Republic</footnote>

<footnote><sup>*</sup>Corresponding author. Faculty of Biosciences and Aquaculture, Nord University, 8049 Bodø, Norway. Tel: +47 90098564; E-mail: christopher.j.hulatt@nord.no</footnote>

<footnote>**Editor:** [Liane G. Benning]</footnote>
A few studies have described natural snow algae communities in mainland Norway (Kol 1963, Procházková et al. 2019a, Tucker and Brown 2022) and neighbouring Sweden (Kol 1974, Lutz et al. 2016). Kol (1963) observed snow algae communities at Finse (southwest Norway), where she found cells representing *Sanguina* (previously assigned to *Chlamydomonas nivalis* Wille), and *Raphidonema nivale* Lagerh. Procházková et al., (2019a) and Tucker and Brown (2022) recently characterized the dispersal of the *Sanguina* population worldwide, including environmental samples from southern Norway. In total, though, very little is known about snow algae biodiversity in this region, and no cultures are established from continental Norway.

Snow algae have diverse adaptations to the harsh, compound stress conditions found in melting snowfields (Leya 2013). In nature, they are often exposed to very high (inhibitory) levels of photosynthetically active radiation (PAR) coupled with elevated UV light levels, low nutrient concentrations, freeze-thaw cycles, and very low temperatures (Řezanka et al. 2008). Lipids are key molecules that perform structural, metabolic, and signalling functions in cell environmental adaptation, provide taxonomic markers, and distinguish the physiological status of communities in the field (Goold et al. 2015, Lutz et al. 2015, Davey et al. 2019, Lupette and Marechal 2020). The stress conditions prevalent in snow habitats are known to trigger neutral lipid and secondary carotenoid accumulation in cytoplasmic lipid droplets (LDs) in many algal species, including the characteristic red-coloured encysted snow-algae cells observed in nature (Leya et al. 2009, Procházková et al. 2019a). LDs in general are compartments that serve as a transient reserve of acyl-chains, and carotenoid-pigmented LDs may additionally have photoprotective roles by absorbing excessive light and scavenging/quenching reactive oxygen species (Britton et al. 2008, Li-Beisson et al. 2015). Assessment of snow algae lipids and their constituent fatty acids has received only limited research, yet lipids connect taxonomic variation with adaptive ecophysiological processes and more widely to the metabolic characteristics of whole snow algae communities.

In 2019, we obtained samples of red snow from melting summer snowfields near to Svartisen glacier, adjacent to the Arctic Circle in northern Norway, a region in which there are no earlier reports of snow-algae blooms. Our aim in this work was to isolate culturable eukaryotic microalgae from these field samples, to characterize the species present, and to establish them as laboratory cultures. Our primary questions were: (i) are culturable species found in this region globally dispersed cosmopolitan strains, or are they endemic to this region? (ii) how do lipid profiles and intrinsic growth rates, features that determine the construction and metabolic properties of snow algae communities, vary between isolates in diverse conditions?

# Materials and methods

## Sampling and algal strain isolation

Red-coloured snow was collected ~1 km from Svartisen glacier in Northern Norway on June 23rd 2019 at an elevation of ~540 m (66°45′03.9 N 14°05′55.0E). The sampling site is shown in Fig. 1. We have observed snow algae blooms on annual snowpack in this area each summer from 2019–2022, which appear in smaller patches (meters scale) dispersed around a wide area (several km²). Samples were collected in sterile plastic bags and transported to the laboratory within 4 h. Snow samples were thawed gently in a refrigerator at 4°C. Using a microscope, we observed a large pro-

portion of red cells, interspersed with smaller numbers of diverse cell types (Fig. S1). To ensure we could obtain sufficient viable cells two isolation methods were used; (i) enrichment cultures followed by plating on agar and (ii) direct plating on agar. Enrichment cultures were prepared by filtering 500 ml of the thawed liquid samples onto ~1.0 µm pore size 47 mm glass fiber filters (VWR, Norway). Filters were then added to a 500 ml Erlenmeyer flask containing autoclaved Bold's Basal Medium (BBM, Bischoff 1963, pH 6.5) and incubated at 6°C and 50 µmol·m⁻²·s⁻¹ photons PAR in a plant growth chamber. After visible (green) cultures were established, culture broth from the flask was streaked onto 1% agar plates. When colonies appeared, each was transferred onto new plates. This procedure was repeated several times to obtain unialgal cultures. In parallel, thawed snow samples were spread directly on agar plates and colonies were selected as above. All isolates were then compared using microscopy (Olympus BX-43 microscope, Olympus Europa GmbH, Hamburg, Germany) to reduce redundancy and retain apparently unique strains based on cell size, morphology, the presence of cilia, chloroplast structure, and whether cells formed chains. Brightfield images were captured using × 1000 magnification and chlorophyll autofluorescence was imaged by excitation with blue LED light and an Olympus DP28 camera.

We selected eight isolates (named ARK-Snn-19, where n is a number) to study based on their morphological uniqueness. Cultivation for DNA extraction was performed in 250 ml Erlenmeyer flasks with BBM at 6°C at 20–40 µmol·m⁻²·s⁻¹ photons PAR. Strains were also deposited at the Culture Collection of Cryophilic Algae, at the Fraunhofer Institute IZI-BB in Potsdam, Germany (http://cccryo.fraunhofer.de/web/strains/) (Table 1).

## DNA extraction, polymerase chain reaction (PCR), and sequencing

For molecular identification of each strain, we attempted to sequence both the 18S and ITS2 rDNA. The 18S rRNA gene analysis, including extraction and amplification by PCR was performed by Roscoff Culture Collection (Roscoff, France). Genomic DNA was extracted with the Nucleospin Tissue Kit (Macherey-Nagel, Düren, Germany). The quality and quantity of DNA were measured by NanoDrop ND-1000 (Thermo Scientific, MI, USA). PCRs were conducted using 1 µL of DNA extract and 29 µL of a mixture of HotStarTaq Plus Master Mix Kit (Qiagen, Hilden, Germany), Coral Load (from the kit), primers, and H₂O (15/3/0.5/0.5/10, v/v/v/v/v/v). The small subunit (18S) rRNA gene was amplified using universal eukaryotic primers (Euk-A, forward): 5'-AACCTGGTTGATCCTGCCAGT-3', (EukBr, reverse): 5'-TGATCCTTCTGCAGGTTCACCTAC-3' (Medlin et al. 1988), (NSF573, forward): 5'-CGCGGTAATTCCAGCTCCA-3' (Hendriks et al. 1989), and (S69R, reverse): 5'-CCGTCADTTCCCTTTTRAGDTT-3' (Probert et al. 2014). Initial denaturation at 95 °C for 5 min was followed by 35 cycles of denaturation at 95°C for 30 s, annealing at 52°C for 30 s, extension at 72°C for 1.5 min, and final extension at 72 °C for 15 min. PCR products were purified using ExoSAP-IT™ PCR Product Cleanup Reagent (Thermo Scientific, Waltham, United States) and sequenced by Macrogen (Amsterdam, Netherlands).

For ITS2 rDNA sequencing, DNA was extracted from wet cell pellets using an E.Z.N.A.® HP Plant DNA Kit (Omega Bio-tek, Georgia, USA). Bead milling (4000 r/m × 5 min, Precellys 24, Bertin Technologies, Montigny le Bretonneux, France) and 0.1 mm glass beads were used to assist cell lysis. The DNA amount and quality were determined by Nanodrop One (Thermo Scientific, MI, USA). Amplicon library preparation and Illumina sequenc-
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Sampling sites of macroscopically visible red-pigmented snow algae communities, Northern Norway. A 30 cm scale bar is shown (right).</figcaption></figure>

<table><caption>**Table 1.** Eight newly established algal isolates from snow assigned to species level, including NCBI GenBank accession numbers for the 18S rRNA gene and 5.8S rDNA partial + ITS2 rDNA + 28S rDNA partial sequences. Asterisk indicates the potentially new species. The isolates are deposited as corresponding ID in the Culture collection of Cryophilic Algae (CCCryo) in Germany.</caption><thead><tr><th rowspan="2">Isolate number</th><th rowspan="2">Species designation</th><th rowspan="2">CCCryo strain ID</th><th colspan="2">GenBank accession number (sequence length, bp)</th></tr><tr><th>18S rRNA gene</th><th>5.8S rDNA partial + ITS2 rDNA + 26S rDNA partial</th></tr></thead><tbody><tr><td>ARK-S01-19</td><td>Raphidonema nivale</td><td>558-22</td><td>OM729978 (2173)</td><td>OM729984 (380)</td></tr><tr><td>ARK-S02-19</td><td>Raphidonema nivale</td><td>559-22</td><td>OM729979 (1649)</td><td>OM729985 (380)</td></tr><tr><td>ARK-S05-19</td><td>Stichococcaceae sp. ARK-S05-19 *</td><td>560-22</td><td>–</td><td>OM729986 (428)</td></tr><tr><td>ARK-S10-19</td><td>Deuterostichococcus epilithicus</td><td>561-22</td><td>–</td><td>OM729987 (395)</td></tr><tr><td>ARK-S08-19</td><td>Chloromonas sp. ARK-S08-19 *</td><td>562-22</td><td>OM729980 (1687)</td><td>OM729987 (437)</td></tr><tr><td>ARK-S11-19</td><td>Chloromonas reticulata</td><td>–</td><td>OM729981 (1603)</td><td>OM729989 (404)</td></tr><tr><td>ARK-S12-19</td><td>Chloromonas reticulata</td><td>563-22</td><td>OM729982 (1631)</td><td>OM729990 (404)</td></tr><tr><td>ARK-S13-19</td><td>Xanthonema bristolianum</td><td>564-22</td><td>OM729983 (1708)</td><td>–</td></tr></tbody></table>

ing were performed by BGI (Guangdong, China) using 300 bp paired end reads on a HiSeq2500 instrument. The forward primer ITS3: 5'-GCATCGATGAAGAACGCAGC-3' and reverse primer ITS4: 5'-TCCTCCGCTTATTGATATGC-3' (White et al. 1990) were used for fragment amplification. The overlapping read pairs from each fragment were fused with FLASH (Magoč and Salzberg 2011), yielding the 5.8 rDNA partial + ITS2 + 28S rDNA partial sequence from seven out of the eight isolates (Table 1). The sequences were clustered with 'cd-hit-est' at setting '-c 1.00 -sc 1' to quantify the unique sequences and eliminate the possibility of other eukaryotic contaminants (Li and Godzik 2006). The 5.8 rDNA partial + ITS2 + 28S rDNA partial, and 18S rRNA gene sequences were deposited in GenBank (Table 1).

## Phylogenetic analysis

Phylogenetic trees of *Raphidonema*-like strains were constructed based on the sequence alignment of the 18S rRNA gene and alignment of ITS2 rRNA secondary structure, respectively, whilst that of *Stichococcus*-like strains was constructed based only on alignment of ITS2 rRNA sequence structure. Two phylogenetic trees of *Chloromonas*-like strains were constructed based on the sequence alignment of 18S rRNA gene and ITS2 rDNA, whilst that of *Xanthonema*-like strains was constructed based only on 18S rRNA gene sequences. The 18S and ITS2 rDNA sequences of related strains were obtained from NCBI GenBank using BLASTn queries (http://blast.ncbi.nlm.nih.gov, Altschul et al. 1990). The sequences were aligned with MUSCLE (multiple sequence alignment by log-expectation) in MEGAX 10.1.8 (Kumar et al. 2018), then visualized

and trimmed to obtain the same lengths for all sequences. The 18S rRNA gene alignment of *Raphidonema*, *Chloromonas*, and *Xanthonema* contained 45 (1661 bp), 52 (1553 bp) and 23 sequences (1730 bp), respectively. The ITS2 rRNA matrix (secondary structure alignment) of *Raphidonema* and *Stichococcus*-like species consisted of 28 (354 bp) and 41 sequences (396 bp), respectively. The ITS2 rDNA alignment of the *Chloromonas* tree contained 16 sequences (321 bp). Phylogenetic trees were constructed using maximum likelihood (ML) with the best model predicted based on BIC in MEGAX v10.1.8. Outgroups and the selected optimal model are specified in each figure or the legend of the corresponding phylogenetic tree. The bootstrapping values (1000 replicates) and Bayesian posterior probabilities were calculated by MEGAX 10.1.8 and MrBayes 3.2.6 (Ronquist and Huelsenbeck 2003), respectively.

## ITS2 rRNA secondary structure analysis

The secondary structure of the ITS2 rRNA was analyzed to examine species boundaries between *Raphidonema* strains, between *Stichococcus*-like strains, and between ARK-S08-19 and *Chloromonas augustae* SAG 5.73. The ITS2, 5.8S, and large subunit (LSU) stem region of rDNA was annotated using the ITS2 database (Keller et al. 2009). The ITS2 secondary structures were computed using a minimum free energy prediction model with the Mfold server (Zuker 2003). Predicted structures were visualized and only sequence-structures containing the specific features of ITS2 rRNA such as four helices with UGG and GGU motifs and U-U mismatch (Coleman 2007, Caisová et al. 2013) were selected for further analysis. The selected sequence-structure was aligned using ClustalW, and
the compensatory base pair changes (CBCs) were computed and visualized in 4SALE 1.7 (Wolf et al. 2014). The aligned sequence-structure data were exported to MEGAX 10.1.8. The secondary structure was drawn by Visualization Applet for RNA (Varna) version 3.8 (Darty et al. 2009) or 4SALE 1.7, followed by editing with Inkscape 0.92 (http://www.inkscape.org/).

## Cultivation experiments

Three different cultivation experiments were performed to target growth rates, fatty acids and LD accumulation. Each of the experiments was conducted in a temperature-controlled plant growth incubator fitted with cool white fluorescent lights (Termaks AS, Norway). Cultivation vessels and nutrient medium were sterilized by autoclaving (121°C, 20 min).

## Cell density

The cell density was measured by the optical density (OD), which was calibrated against the dry weight (DW). The OD was measured at 540 nm, which avoided chlorophyll absorbance peaks. A 1 cm cuvette was used for OD measurement, with samples diluted below absorbance 1.0 for linear response. In Experiments 1 and 3, the DW was measured at the beginning and end of cultivation by filtering 3–10 ml of culture broth through pre-weighted ~1.0 µm pore size 47 mm glass fiber filters (VWR, Oslo, Norway). Filters were dried at 100°C for 24 h, then re-weighed and the DW (g·L⁻¹) was calculated. For each species, linear calibration was used to calculate the DW from the OD (Fig. S2).

## Experiment 1. Growth rates

The eight isolates were cultivated at 2°C and at 10°C in glass tubes (n = 3 at each temperature), containing 90 ml of 3N-BBM (BBM with 3-fold nitrate, Bischoff 1963, pH 6.8) and illuminated with 109 ±5 µmol·m⁻²·s⁻¹ photons PAR. Each tube was sparged with air enriched with 1.0% CO₂ at 62.5 ml·min⁻¹, supplied by a gas mixing system (Photon Systems Instruments GMS 150, Drásov, Czech Republic). OD measurements were taken every 1–2 days and at the end of experiments cells were pelleted and stored at –40°C for further analysis.

The growth rate of each strain was established by fitting the cell density measurements to a logistic growth equation using the R package 'Growthcurver' (Sprouffske and Wagner 2016) in R software version 3.1.2 (Team, 2018).

$$N_t = \frac{K}{1 + \left( \frac{K - N_0}{N_0} \right) e^{-rt}}, \quad (1)$$

where, N_t is the population size (cell density, g·L⁻¹) at time t (days), N₀ is the initial cell density (g·L⁻¹), K is the carrying capacity (maximum cell density, g·L⁻¹) and r is the population growth rate (µ·d⁻¹). Growth curves were fitted independently to each replicate cultivation, and the mean ± SD of n = 3 is presented for each of the eight isolates at each temperature.

## Experiment 2. Fatty-acid profiling

To determine the fatty-acid profile of each isolate under standardized non-stressed conditions, cells were cultured at 6°C under 50 µmol·m⁻²·s⁻¹ photons PAR in 500 ml Erlenmeyer flasks (n = 4 each) containing 200 ml of BBM (pH 6.8). Cells were harvested by centrifugation (3,000 g, 3 min) when the OD reached absorbance 0.25, ensuring the cell density was well below the carrying capacity of the nutrient medium and each culture received equiva-

lent light supply at the time of collection. Samples were stored at –40°C.

## Experiment 3. LD accumulation in ARK-S12-19

Strain ARK-S12-19 was cultivated at 10°C in 27 mm diameter glass bubble tubes containing 380 ml of 3N-BBM (pH 6.8). The tubes were supplied with air containing 1.0% CO₂ at 85 ml·min⁻¹ and 145 ±8 µmol·m⁻²·s⁻¹ photons PAR. The cells were first grown for 9 days, then the cultures were gently centrifuged (1920 g, 15 min) and resuspended in nitrate (N)-free media to initiate the start of the experiment (day 0), providing an accurate onset of N-starvation. The subsequent experimental cultivation period in N-free medium was 11 days, with biomass, imaging flow cytometry and lipid samples collected at day 0, 1, 2, 4, 7, and 11.

## Imaging flow cytometry

An ImageStream®X Mk II Imaging Flow Cytometer (Luminex Corporation, Austin, TX, USA) was used to measure cell dimensions, LD accumulation and chlorophyll content in Experiment 3. Cell LDs were labelled with Bodipy® 493/503 (Molecular Probes, Eugene, OR, USA) dissolved in 100% DMSO (v/v) at a final concentration of 0.075 µg mL⁻¹, mixed and incubated for 15 mins. A 488 nm laser operated at 1.0 mW was used for excitation of the probe and chlorophylls. Detector wavelengths for the different channels were 505–560 nm (Bodipy® label, Channel 2); 595–642 nm (brightfield image, Channel 4); 642–745 nm (chlorophyll autofluorescence, Channel 5). A total of up to 10000 particles were analyzed at a rate of ~50 objects/second at magnification × 40. IDEAS v.6.1.823.0 software was used for downstream analysis, including gating to select the main population of single cells and exclude debris or doublets. Image masking was conducted to measure the length, width, area and fluorescence intensity of each individual cell. Cell dimensions and area were measured from the brightfield images (Channel 4) by applying an image mask, and the fluorescence intensity of Channel 2 and Channel 5 were determined from the same mask area (Fig. S13). The biovolume of each cell was calculated from the length/width dimensions of the image using the formula for a prolate spheroid shape (Hillebrand et al. 1999).

$$V = \frac{\pi}{6} \cdot d^2 \cdot h, \quad (2)$$

where, V is the volume (µm³), d is diameter (µm), and h is height (µm) of the cell. Both the cell area and cell volume were tested as normalization values, by dividing Bodipy® fluorescence by either the area or volume.

## Extraction, isolation of lipids, and methylation

For Experiment 1, we directly extracted and derivatized the fatty acids to fatty-acid methyl esters (FAMES) from cells cultivated at 10°C as described by Radakovits et al. (2012). For Experiments 2 and 3, total lipids were extracted from freeze-dried cell pellets using chloroform and methanol as described by Breuer et al. (2012). For Experiment 3, neutral lipids and polar lipids were additionally separated by solid-phase extraction (Sep-Pak 1 g silica cartridges, 6 ml, Waters, MO, USA) followed by the separation of TAG (Triacylglycerol) from neutral lipids by high-performance thin-layer chromatography (HPTLC). Fatty acids in each lipid fraction were derivatized by adding 3 ml of 5% H₂SO₄ in methanol and incubation at 70°C for 3 h. After washing with NaCl-saturated water, FAMES were recovered in hexane. For a full description of the experiments, see "Supplementary Material" section.
## FAME analysis

FAMES were quantified by a Gas Chromatograph (GC) equipped with a Flame Ionization Detector (FID, SCION 436, Bruker, UK) and an Agilent CP-Wax 52CB column (Agilent Technologies, USA). Supelco 37 component standards (Sigma-Aldrich, MA, USA) were used for identification and quantification of the FAMES, with additional standards for unsaturated C16 and C18 fatty acids. To identify unusual fatty acids, the FAMES structures of ARK-S01-19 and ARK-S13-19 were verified by Helsinki University Lipidomics Unit (HiLIPID). The analyses were performed based on their electron-impact mass spectra recorded by GC-MS-QP2010 Ultra (Shimadzu Scientific Instruments, Kyoto, Japan), compared to the spectra of analytical standard mixtures of FAMES (Supelco® Analytical Products, Merck) and published reference mass spectra (Christie 2021). The GC-mass spectrometer (GC-MS) was equipped with a Zebron ZB-wax capillary column (30 m, 0.25 mm ID and film thickness 0.25 µm; Phenomenex, Torrence CA, USA), which gave the FAMES similar retention time patterns as the column used for the quantitative GC-FID analyses.

## Results

### Snow-algae identification: cell morphology and phylogeny

Morphological features (Fig. 2, Fig. S3) analyzed using light microscopy together with sequence data of molecular markers were used to identify eight snow algae isolates. We constructed different phylogenetic trees to place the genus *Raphidonema* (2 isolates), *Chloromonas* (3 isolates), *Xanthonema* (1 isolate), and the family *Stichococcaceae* (2 isolates) using the sequences presented in Table 1.

### Raphidonema-like species

The cells of ARK-S01-19 were filiform or elongated fusiform, straight 8.1–45.4 µm long and 2.0–3.9 µm wide with bent aristate or rotund apices (Fig. 2 a1–a5). Cells were usually unicellular, but sometimes remained united in short chains up to four cells, with ribbon-like chloroplasts (Fig. S3a). Filaments with branched acute ends were observed. ARK-S02-19 cells were filiform, straight 15.1–47.6 µm in length and 2.3–3.8 µm in width with acute apices or rotund apices (Fig. 2 b1 and b2). The cells commonly formed unbranched filaments consisting of more than four and up to sixteen cells with plate-like chloroplasts (Fig. S3b).

The ITS2 rDNA sequence of ARK-S01-19 was identical to R. nivale (Lagerheim) CCCryo 381-11 and CCCryo375-11 isolated from a red snowfield (77°0′N 16°37′E) and on top of a steep snowfield on Svalbard (77°0′N 16°22′E), respectively. The 18S rRNA gene sequence of ARK-S02-19 was identical to R. nivale strain CCAP 470/4 isolated from a mountain at an altitude of 1295 m in Washington State (USA), with only 1 bp difference in ITS2 rDNA. We constructed a ML tree based on 18S rRNA gene sequences, which revealed that each isolate formed a different subclade together with other *Raphidonema* within a clade with full statistical support (ML/BL:100/1.00), being a sister to the strains identified as *Pseudochlorella* (Fig. S4). ITS2 rRNA sequence-structure analysis was further used for species-level resolution, where ML analysis showed that our two isolates fell into a clade together with other R. nivale (ML:82) (Fig. 3). The worldwide distribution of *Raphidonema* strains in the phylogenetic tree (Fig. 3) is mapped in Figure S5, indicating the geographic location of the closest relatives to the new isolates and highlighting the known distribution of R. nivale

in maritime Antarctica (James Ross), Western USA, Svalbard, continental Norway, and Switzerland.

### Stichococcus-like species

The two *Stichococcus*-like isolates ARK-S05-19 and ARK-S10-19 are both cylindrical cells with rounded ends (Fig. 2C and D). The cells of ARK-S05-19 are 4.1–8.9 µm in length and 2.3–3.6 µm in width, and the length to width ratio was around 2.2 (Fig. 2 c1–c3). The cells were usually unicellular, but sometimes occurred as unbranched two-cell filaments with plate-like chloroplasts (Fig. S1c). In contrast, ARK-S10-19 cells were 5.2–14.8 µm long and 2.2–3.2 µm wide, resulting in the length to width ratio of 3.0, and usually unicellular but occasionally appearing as short unbranched filaments of 2 to 8 cells with plate-like chloroplasts (Fig. 2 d1–d3 and Fig. S3d).

The 18S rRNA gene sequences of *Stichococcus*-like species were unfortunately unobtainable due to unsuccessful amplification by PCR (Table 1). However, the ITS2 rRNA sequence-structure based phylogenetic tree showed that ARK-S10-19 belonged to a *Deuterostichococcus* clade (Pröschold and Darienko 2020). Furthermore, its 5.8S rDNA partial + ITS2 + 28S rDNA partial sequences were identical to that of *D. epilithicus* FG2/4.2 (KM020048), for which the isolation/sample location or habitat are unknown. *Deuterostichococcus epilithicus* have been found worldwide, many of which are isolated from soil or subaerial habitats (Fig. S6). In contrast, the closest ITS2 rDNA of ARK-S05-19 in the NCBI database was the unclassified isolate Prasiolales sp. S2RM26, originating from the surface of snowfall samples in Sweden (Tesson and Šantl-Temkiv 2018), revealing a difference of 49 nucleotides out of 318 (85% identity). ARK-S05-19 and the unclassified isolate Prasiolales sp. S2RM26 together formed a well-supported (93%) independent clade distinct from the sister lineages *Desmococcus olivaceus* and *S. bacillaris* (Fig. 4). Nevertheless, one CBC difference was detected in helix I between ARK-S05-19 and Prasiolales sp. S2RM26 (Fig. 5). Several strains of *S. bacillaris* and five strains of *D. olivaceus* also provided the highest BLAST matches for the ITS2 rDNA sequence of ARK-S05-19, but their rRNA secondary structures differed by three and four CBCs, respectively (Table S1, Table S2).

### Chloromonas-like species

The cells of isolates ARK-S11-19 and ARK-S12-19 were both ellipsoid-shaped, 12–13 µm in length and 6–7 µm in width (Fig. 3F and G), having a prominent hemispherical shaped papilla (Fig. 2 f6 and g6) and a cup-shaped chloroplast without a pyrenoid (Fig. 2 f1, f2, g1, g2 and Fig. S3 f, g). The eyespot was in the anterior 1/2–1/3 of the cell (Fig. 2 f1 and g1) and two contractile vacuoles were located on anterior side (Fig. 2 f5 and g5). The 18S rRNA gene sequence of ARK-S11-19 and ARK-S12-19 (1603 and 1631 bp long, respectively), and the 5.8S rDNA partial + ITS2 + 28S rDNA partial sequences of both isolates (404 bp long each) were identical to the type species *C. reticulata* CCAP 11/110 (= CCCryo 213-05; SAG 29.83; UTEX 1970), a snow alga isolated from Todd lake, Oregon, USA. The phylogeny based on the 18S rRNA gene sequences using the ML method showed that both isolates belong to core *Chloromonas*, located within a subclade containing the type species and corresponding only to the *C. reticulata* clade (Pröschold et al. 2001) within clade 1 (Hoham et al. 2002) (Fig. 6A). Further phylogenetic analysis based on the ITS2 rDNA provided higher resolution and both isolates were clustered together with *C. reticulata*, within a well-supported (ML/BI: 99/0.91) subclade of the *Reticulata* group (Matsuzaki et al. 2012, Fig. 6B). The biogeographical location and habitat of selected *Chloromonas* strains (shown in Fig. 6A and B) is
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Microphotographs of eight algal isolates: R. nivale ARK-S01-19 (a1, a2, a3, a4, a5); R. nivale ARK-S02-19 (b1, b2); Stichococcaceae sp. ARK-S05-19 (c1, c2, c3); D. epilithicus ARK-S10-19 (d1, d2, d3); *Chloromonas* sp. ARK-S08-19 (e1, e2, e3, e4, e5, e6, e7, e8); C. *reticulata* ARK-S11-19 (f1, f2, f3, f4, f5, f6); C. *reticulata* ARK-S12-19 (g1, g2, g3, g4, g5, g6, g7, g8); X. *bristolianum* ARK-S13-19 (h1, h2, h3, h4, h5). Arrowheads indicate the papillae and arrows indicate contractile vacuoles. Scale bars are either 2 µm (with indication of the scale size near the bar) or 5 µm (no indication of the scale size). p = pyrenoid.</figcaption></figure>

presented in Fig. S7, showing a multi-continental distribution of
*C. reticulata* species.

Microscopic analysis revealed that the vegetative cells of ARK-
S08-19 are round or ellipsoidal in shape, 10–20 µm in length and 7–
8 µm in width, have a small hemispherical shaped anterior papilla
(Fig. 2 e1), a nucleus and asteroid chloroplast with centrally lo-
cated pyrenoid (Fig. 2 e1, e2, e8 and Fig. S1e) and two contractile
vacuoles on the anterior side (Fig. 2e7). The vegetative cells have

12–14 µm long pairs of cilia and the eyespots are located on the
anterior 1/3 of the cell (Fig. 2 e1 and e6). The non-motile cells of
ARK-S08-19 are oval or spherical, up to 28 µm in diameter with
pyrenoid (Fig. S8). The formation of two to sixteen zoospore cells
within the parental cell wall was observed. Using the 18S rRNA
gene sequence phylogeny, isolate ARK-S08-19 was located in a
clade together with *C. augustae* within clade 1 (Hoham et al. 2002,
Fig. 6A). Based on available GenBank data, the 18S rRNA gene se-
<figure><img src="image_4.png" /><figcaption>**Figure 3.** ITS2 rRNA sequence-structure ML tree of *Raphidonema* with the new isolates (sequence) shown in bold. An asterisk indicates the type strain of the species *R. nivale* and the species names are as proposed by Yakimovich et al. (2021). The best model was JTT+G+I+F calculated by MEGAX 10.1.8. Numbers next to branches indicate statistical support values [ML bootstraps (1000 replicates)/Bayesian posterior probabilities]. The bootstrap support values above 70% and Bayesian posterior probabilities above 0.90 are shown. Thick lines indicate branches with full statistical support (ML/BI:100/1.00). *Stichococcus bacillaris* SAG 249.80 (MT078156) and *D. epilithicus* SAG 10.97 (MT078169) serve as the outgroup.</figcaption></figure>

quence of ARK-S08-19 was most closely related to *C. augustae* SAG 13.89 isolated from a salt marsh pool in the United Kingdom, with 13 bp difference out of 1687 bp (99.23% similarity), followed by *C. augustae* SAG 5.73 (99.17% similarity), a freshwater alga isolated from the former Czechoslovakia, differing in 14 bp out of 1687 bp (99.17% similarity). The closest and only BLASTn hit of the ITS2 rDNA sequence of ARK-S08-19 against the NCBI database was *C. augustae* SAG 5.73, yet it differed by 32 bp out of 136 bp of aligned sequence (40% query cover, 76.5% similarity). Further ITS2 rRNA sequence-structure analysis showed that the alignable sequences were mainly found in helix II and near the 5' side of helix III (Fig. S9).

## Xanthonema-like species

The ARK-S13-19 cells were oblong with rounded corners, ~8.4 µm in length 5.5 µm in width (Fig. 2 h1 to h5). The cells usually form unbranched filaments consisting of 2 to 30 cells, but it can sometimes be unicellular (Fig. 2 h4), having four to 12 plate-like chloroplasts (Fig. S3h). Eyespots were occasionally visible in some cells (Fig. 2h3 and h5) indicating the development of motile cells. Phylogenetic analysis of the 18S rRNA gene sequence using ML methods showed that ARK-S13-19 belongs to the fully statistically supported *Xanthonema* 1 clade (ML/BI:100/1.00) within the Tribonematales (Fig. 7) (Maistro et al. 2009). The 18S rRNA gene sequence of ARK-S13-19 (1708 bp long) was identical to *X. bristolianum* CCALA 516 isolated from snow in Belianske Tatras, Slovakia. The closely related strain of *X. bristolianum* ACSSI 290 was isolated from soil in Russia.

## Growth rate

The growth patterns of the eight isolates at 2°C and 10°C were modelled (Equation 1) and the results are presented in Fig. 8 and

Table S3. At 10°C, the growth rate $r$ varied from 0.45 ±0.01 to 0.86 ±0.10 d⁻¹ amongst all isolates, with the highest growth rates recorded in *C. reticulata* ARK-S11-19 and ARK-S12-19. *X. bristolianum* ARK-S13-19 exhibited the highest carrying capacity (maximum cell mass production, K) reaching 7.5 ±0.25 g·L⁻¹ whilst the lowest maximum cell production was achieved by *Chloromonas* sp. ARK-S08-19 with a K estimate of 1.3 ±0.08 g·L⁻¹. At lower 2°C, close to the temperature of melting snow, all strains exhibited an extended lag phase and the growth rate of seven strains was substantially reduced compared to growth at 10°C. At 2°C maximum growth rates ranged from 0.20 ±0.04 to 0.56 ±0.10 d⁻¹, with *R. nivale* ARK-S01-19 exhibiting the highest specific growth rate amongst the isolates, comparable to that obtained at 10°C. In contrast, the related isolate *R. nivale* ARK-S02-19 recorded protracted and substantially reduced growth rates of 0.23 ±0.02 d⁻¹ at 2°C. *C. reticulata* ARK-S11-19 and ARK-S12-19 also recorded an extended lag phase, but eventually established intermediate growth rates of 0.34 ±0.03 and 0.36 ±0.01 d⁻¹, respectively. At 2°C *X. bristolianum* ARK-S13-19 showed minimal growth over the whole cultivation period. Interestingly, ARK-S08-19 reached a higher cell density (K) of 2.0 ±0.21 g·L⁻¹ at 2°C than that of 1.3 ±0.08 g·L⁻¹ obtained at 10°C, and its growth rate at low temperature was not suppressed to the same extent as most of the other strains.

## Fatty-acid profiles

The total fatty acid (TFA) contents of each isolate cultivated under standardized conditions with excess nutrients (Experiment 2) ranged from 54.3 to 93.0 mg·g⁻¹ DW, with total polyunsaturated fatty acids accounting for between 57.3% and 70.4% of TFA (Table 2). The fatty acid profiles of both *R. nivale* isolates, ARK-S01-19 and ARK-S02-19, were distinguished by substantial amounts of C18:3n-3, accounting for 25.2 and 19.1% of TFAs
<table><caption>**Table 2.** TFA and fatty acid composition of eight isolates in nutrient-replete conditions (Experiment 2) (mg.g⁻¹ DW). Mean values (±SD) of triplicates per isolate are given in the table.</caption><thead><tr><th></th><th>R. nivale<br>ARK-S01-19</th><th>R. nivale<br>ARK-S02-19</th><th>Stichococcaceae<br>sp. ARK-S05-19</th><th>Chloromonas sp.<br>ARK-S08-19</th><th>D. epilithicus<br>ARK-S10-19</th><th>C. reticulata<br>ARK-S11-19</th><th>C. reticulata<br>ARK-S12-19</th><th>X. bristolianum<br>ARK-S13-19</th></tr></thead><tbody><tr><td>C14:0</td><td>2.9 ± 0.05</td><td>2.6 ± 0.5</td><td>5.1 ± 0.1</td><td>2.2 ± 0.1</td><td>2.8 ± 0.1</td><td>2.2 ± 0.3</td><td>2.4 ± 0.1</td><td>9.9 ± 1.2</td></tr><tr><td>C16:0</td><td>6.7 ± 0.7</td><td>6.9 ± 1.1</td><td>4.6 ± 0.4</td><td>7.7 ± 0.2</td><td>3.9 ± 0.6</td><td>10.1 ± 0.7</td><td>9.8 ± 0.5</td><td>4.5 ± 0.7</td></tr><tr><td>C16:1n-7</td><td>1.2 ± 0.1</td><td>2.5 ± 0.5</td><td>1.8 ± 0.1</td><td>1.7 ± 0.02</td><td>1.0 ± 0.03</td><td>1.0 ± 0.1</td><td>0.7 ± 0.5</td><td>12.1 ± 1.6</td></tr><tr><td>3t-16:1</td><td>1.6 ± 0.1</td><td>1.5 ± 0.3</td><td>2.6 ± 0.1</td><td>1.4 ± 0.03</td><td>1.5 ± 0.1</td><td>1.3 ± 0.1</td><td>1.4 ± 0.1</td><td>3.3 ± 0.4</td></tr><tr><td>C16:2n-6</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.5 ± 0.1</td></tr><tr><td>C16:2n-4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>5.7 ± 0.7</td></tr><tr><td>C16:3n-4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>9.3 ± 1.2</td></tr><tr><td>C16:3n-3</td><td>1.8 ± 0.1</td><td>1.1 ± 0.1</td><td>12.2 ± 0.3</td><td>3.3 ± 0.1</td><td>6.5 ± 0.5</td><td>4.8 ± 0.3</td><td>4.6 ± 0.5</td><td>—</td></tr><tr><td>C16:4n-3</td><td>4.3 ± 0.1</td><td>4.8 ± 1.1</td><td>—</td><td>13.6 ± 0.3</td><td>—</td><td>16.6 ± 0.9</td><td>16.3 ± 1.0</td><td>—</td></tr><tr><td>C18:0</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C18:1n-9</td><td>3.2 ± 0.2</td><td>2.9 ± 0.2</td><td>—</td><td>13.2 ± 0.9</td><td>2.3 ± 0.4</td><td>7.8 ± 0.7</td><td>8.4 ± 01.4</td><td>—</td></tr><tr><td>C18:1n-7</td><td>4.3 ± 0.2</td><td>5.9 ± 1.0</td><td>7.3 ± 0.2</td><td>10.5 ± 0.2</td><td>3.6 ± 0.1</td><td>9.3 ± 0.6</td><td>9.0 ± 0.6</td><td>—</td></tr><tr><td>C18:2n-6</td><td>0.6 ± 0.1</td><td>1.9 ± 0.2</td><td>2.3 ± 0.1</td><td>1.3 ± 0.05</td><td>2.4 ± 0.4</td><td>0.7 ± 0.05</td><td>0.5 ± 0.2</td><td>2.0 ± 0.3</td></tr><tr><td>NMI-18:3</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>0.5 ± 0.1</td></tr><tr><td>C18:3n-3</td><td>15.1 ± 1.4</td><td>10.6 ± 1.9</td><td>26.6 ± 0.7</td><td>38.2 ± 0.5</td><td>17.7 ± 1.0</td><td>37.3 ± 2.3</td><td>36.1 ± 0.31</td><td>—</td></tr><tr><td>C18:3n-4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>1.6 ± 0.2</td></tr><tr><td>C18:4n-4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>4.0 ± 0.6</td></tr><tr><td>NMI-18:4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>2.8 ± 0.3</td></tr><tr><td>C18:4n-3</td><td>2.0 ± 0.1</td><td>4.1 ± 0.8</td><td>—</td><td>—</td><td>1.3 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C20:0</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C20:1n-9</td><td>0.7 ± 0.03</td><td>1.2 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C20:3n-6</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C20:4n-6</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C20:4n-3</td><td>1.7 ± 0.1</td><td>—</td><td>2.5 ± 0.1</td><td>—</td><td>5.0 ± 0.9</td><td>—</td><td>—</td><td>3.1 ± 1.1</td></tr><tr><td>C20:5n-3</td><td>11.8 ± 0.3</td><td>9.2 ± 1.7</td><td>5.0 ± 0.2</td><td>—</td><td>1.6 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C22:0</td><td>0.1 ± 0.3</td><td>—</td><td>—</td><td>—</td><td>3.7 ± 0.2</td><td>—</td><td>—</td><td>27.3 ± 3.3</td></tr><tr><td>C24:0</td><td>0.4 ± 0.2</td><td>0.1 ± 0.1</td><td>0.6 ± 0.1</td><td>—</td><td>0.8 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C24:0</td><td>1.6 ± 2.0</td><td>—</td><td>—</td><td>—</td><td>0.3 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>C22:6n-3</td><td>38.2 ± 0.9</td><td>29.7 ± 5.4</td><td>43.8 ± 1.1</td><td>55.2 ± 0.9</td><td>30.8 ± 1.8</td><td>58.7 ± 3.4</td><td>57.1 ± 4.5</td><td>2.5 ± 0.3</td></tr><tr><td>Omega-3</td><td>0.6 ± 0.1</td><td>1.9 ± 0.2</td><td>4.7 ± 0.1</td><td>1.3 ± 0.05</td><td>7.4 ± 1.0</td><td>0.7 ± 0.04</td><td>0.5 ± 0.2</td><td>30.3 ± 3.7</td></tr><tr><td>Omega-6</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>5.6 ± 1.3</td></tr><tr><td>Omega-4</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>23.5 ± 3.1</td></tr><tr><td>SFA</td><td>10.1 ± 0.3</td><td>9.5 ± 1.5</td><td>10.1 ± 0.4</td><td>9.6 ± 0.4</td><td>7.8 ± 0.7</td><td>11.8 ± 0.9</td><td>11.9 ± 0.7</td><td>14.2 ± 1.9</td></tr><tr><td>MUFA</td><td>11.0 ± 0.3</td><td>14.0 ± 1.9</td><td>11.7 ± 0.3</td><td>26.9 ± 1.0</td><td>8.3 ± 0.4</td><td>19.4 ± 1.5</td><td>19.5 ± 1.5</td><td>15.4 ± 2.0</td></tr><tr><td>PUFA</td><td>38.8 ± 0.9</td><td>31.7 ± 5.5</td><td>48.6 ± 1.1</td><td>56.4 ± 0.8</td><td>38.2 ± 1.9</td><td>59.8 ± 3.5</td><td>57.6 ± 4.6</td><td>59.4 ± 7.4</td></tr><tr><td>TFAs</td><td>59.9 ± 4.1</td><td>55.1 ± 3.2</td><td>70.4 ± 7.4</td><td>93.0 ± 11.3</td><td>54.3 ± 4.2</td><td>91.1 ± 11.1</td><td>88.9 ± 10.8</td><td>89.0 ± 6.9</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 4.** ITS2 rRNA sequence and secondary structure-based ML phylogeny of *Stichococcus*-like species with the new isolates (sequences) shown in bold. The best model was WAG+G+F, calculated by MEGAX 10.1.8. The asterisk shows the type strain *Deuterostichococcus epilithicus* according to Pröschold and Darienko (2020). Numbers next to branches indicate statistical support values [ML bootstraps (1000 replicates)/Bayesian posterior probabilities]. The bootstrap support values above 70% and Bayesian posterior probabilities above 0.9 are shown. Thick lines indicate the branches with full statistical support (ML/BI:100/1.00). *Pseudostichococcus monaillantoides* SAG 379–4 (MT078184) and *P. monaillantoides* SAG 380–1 (KM020066) serve as the outgroup.</figcaption></figure>

(15.1 and 10.6 mg·g⁻¹ DW) respectively, and C20:5n-3 (eicosapentaenoic acid, EPA) amounting to 19.7 and 16.6% of TFAs (11.8 and 9.2 mg·g⁻¹ DW), respectively. C16:0 was also found at intermediate levels (11.2 and 12.5% TFAs, 6.7 and 6.9 mg·g⁻¹ DW, respectively). The two R. nivale strains differed slightly, with ARK-S01-19 containing additional trace amounts of C22:0, C20:4n-3, and C22:6n-3 (docosahexaenoic acid; DHA), which were absent in ARK-S02-19 (Table 2).

Stichococcaceae sp. ARK-S05-19 and D. epilithicus ARK-S10-19 were characterized by high amounts of C18:3n-3, accounting for 37.7% and 32.5% of TFAs (26.6 and 17.7 mg·g⁻¹ DW), respectively. C16:3n-3 was also abundant, accounting for 17.4% and 11.9% TFAs (12.2 and 6.5 mg·g⁻¹ DW), respectively.

The FAME profiles of *Chloromonas* isolates ARK-S08-19, ARK-S11-19 and ARK-S12-19 revealed the dominance of C18:3n-3,

which accounted for 40.6%–41.1% of TFAs, with additional fatty acids including C16:4n-3, C16:0, C18:1n-7, and C18:1n-9. Although profiles were broadly comparable, some differences in fatty acid proportions were observed, for example *Chloromonas* sp. ARK-S08-19 had higher proportions of C18:1n-9 (14.2% TFAs) than *C. reticulata* ARK-S11-19 and ARK-S12-19 (8.6 and 9.4% TFAs, respectively).

In contrast to the other isolates investigated in this study, *X. bristolianum* was characterized by the presence of the omega-4 fatty acids C16:2n-4, C16:3n-4, C18:3n-4, and C18:4n-4. In addition, *X. bristolianum* contained unusual (non-methylene interrupted) NMI-18:3 and NMI-18:4 fatty-acids (Fig. S10), of which the mass spectrum was not included in one of the best-known existing libraries of C18 fatty acids (Christie 2021). The most abundant FA in *X. bristolianum* was EPA (30.7% TFAs) followed by C16:1n-7 and C14:0.
<figure><img src="image_6.png" /><figcaption>**Figure 5.** The consensus structure of ITS2 rRNA sequence secondary structure between ARK-S05-19 (OM729986) and Prasiolales sp. S2RM26 (MK005091). Partial 5.8S and partial LSU stem regions of rRNA are included. One compensatory base change (CBC) between these two species within helix I is indicated in pink. Nucleotides highlighted in circles outside the yellow mark indicates a CBC between ARK-S05-19 (white circles) vs Prasiolales sp. S2RM26 (black circles). The conserved regions are indicated in blue with nucleotides (A, U, G, or C) and non-conserved nucleotides are shown in yellow. The ITS2 secondary structure and sequences were synchronously aligned and visualized using 4SALE.</figcaption></figure>

FAME profiles were also examined in the algal samples collected from the stationary phase of the growth rate experiment at 10°C (Experiment 1, Fig. S11). The amount of TFAs was higher in all isolates (155.3–411.6 mg·g⁻¹ DW) than that in Experiment 2 and the FA profile of each isolate was characterized by accumulation of a few fatty-acids. Substantial proportions of C18:1n-9 were observed in *Raphidonema*, *Stichococcus*, and *Chloromonas* strains accounting for the majority (48.7–57.2%) of TFAs. *Xanthonema* sp. ARK-S13-19 accumulated primarily C16:1n-7 reaching 50.2% of TFAs, followed by C16:0 (17.3% TFAs), C14:0 (16.8% TFAs), and only a small amount of C18:1n-9.

## LD dynamics in C. reticulata

We selected cosmopolitan *C. reticulata* ARK-S12-19 as a model to study LD formation under N-starvation, using both chemical lipid analysis by GC-FID and image-based labelling with imaging flow cytometry (Experiment 3). Over the 11-day experimental period, neutral and polar lipids were separated by solid-phase extraction (Fig. S12), and TAG was further separated from other neutral lipids by thin layer chromatography (TLC, Fig. 9A). Polar lipids peaked at 127.6 mg·g⁻¹ DW at day 7, and decreased to 75.0 mg·g⁻¹ DW after

11 days of N removal, indicating limited membrane lipid degradation towards the end of the experiment (Fig. S12). During the same period, total neutral lipids increased from 14.2 to 208.7 mg·g⁻¹ DW. TLC analysis revealed that C16:4n-3 and C18:3n-3 fatty acids dominated TAG at the start of the experiment, but by day 11 fatty acids C18:1n-9 and C18:1n-7 became the most abundant (Fig. 9A). Using imaging flow cytometry (IFC) we were able to distinguish different cell types including oval shaped bicilliates (27.5% ±5.2%; R4) and cells lacking cilia (35.6% ±6.5%; R3) (Fig. S13). Bodipy® probe fluorescence showed an increase in neutral lipids over time, which was confirmed by visualizing the size, position and number of LDs within the cells (Fig. 9B). As variation in cell size could influence the quantitation of neutral lipids by IFC, we tested whether the normalization of Bodipy® fluorescence to the image bright field area (µm²) and biovolume of individual cells (µm³) improved the correlation between IFC (Bodipy®) and GC quantification (SPE and TLC) of fatty acids (Fig. S14 and S15). Normalizing the Bodipy® 493/503 fluorescence to either 2D cell image area or cell biovolume resulted in improved correlation with TFAs, neutral lipids, and TAG compared to non-normalized Bodipy® fluorescence (Fig. S15). The correlation coefficients were slightly higher for neutral lipids vs Bodipy® (r² = 0.97) than TAG vs Bodipy® (r² = 0.95) (Fig. S15).

## Discussion

Snow-algae communities harbour substantial biodiversity, yet their phylogeny and global distribution are not well understood. This study is the first to isolate and cultivate snow algae from mainland Norway. Our use of relatively conserved (18S rRNA gene) and highly variable (ITS2 rDNA) genetic markers combined with microscopy and lipid profiling will assist detection of these species widely in environmental samples.

## Phylogeny using molecular markers

*Raphidonema* is a cosmopolitan, globally dispersed genus (Segawa et al. 2018) commonly found in alpine or glacier ecosystems (Leya et al. 2009, Yakimovich et al. 2021, Novis et al. 2023, Segawa et al. 2023). Using morphology alone, *Raphidonema* species are difficult to distinguish from other closely related Trebouxiophycean genera, especially *Koliella*, due to their similar cell appearance combined with phenotypic plasticity under different conditions (Hoham 1973, Stibal and Elster 2005). The phylogeny of *Koliella* and *Raphidonema* was recently resolved with molecular markers, where strains designated as *Raphidonema* or *Koliella* were clearly separated based on *rcbL* sequences (Yakimovich et al. 2021). The authors further revised *Raphidonema* using the ITS2 rRNA sequence-structure and proposed five species within the genus. Our analyses showed that the 18S rRNA gene was not adequate in delineating the isolates at species level, but ITS2 rRNA sequence-structure analysis provided high resolution and places ARK-S01-19 and ARK-S02-19 as different strains of *R. nivale*. Most of the existing *Raphidonema* cultures in collections were isolated from snow and glacier habitats (CCCryo culture collection, the Culture Collection of Algae at the University of Texas at Austin, UTEX), but Stibal and Elster (2005) reasonably suggest that *Raphidonema* is a soil-dweller transported onto snow by wind or meltwater. Although *R. nivale* has not yet been cultured from soil, closely related species, including *R. pyrenoidifera* and *R. monicae* have been isolated from this habitat (Yakimovich et al. 2021). Further investigation of summer Arctic and Alpine soil communities (e.g. by amplicon se-
<figure><img src="image_7.png" /><figcaption>**Figure 6. (A)** 18S rRNA gene-based ML phylogenetic tree of *Chloromonas*. The asterisk shows the authentic strain of *C. reticulata*, the type species of the genus *Chloromonas*. Our new isolates (sequence) are shown in bold. *Chloromonas* clade 1,2 and 3 are delimited according to Hoham et al. (2002). (B) ITS2 rDNA ML phylogenetic tree of the *Reticulata* group showing the phylogenetic relationship of ARK-S11-19 and ARK-S12-19 with other *C. reticulata*. For both phylogenetic trees, the best model was K2+I calculated by MEGAX 10.1.8. Numbers next to branches indicate statistical support value [ML bootstraps (1000 replicates)/Bayesian posterior probabilities]. The bootstrap support values above 70% and Bayesian posterior probabilities above 0.90 are shown. Thick lines indicate the branches with full statistical support (ML/BI:100/1.00).</figcaption></figure>

quencing) could improve our understanding of the lifestyle and global distribution of *Raphidonema*.

Stichococcus-like strains belonging to the Prasiola clade are widely distributed in various habitats, from tropical to polar regions, freshwater to marine, and from hot acidic springs to snow and cryoconite hole environments (Butcher 1952, Stibal et al. 2006, Hodač et al. 2016). They often change their morphotype, making it difficult to designate strains at genus and species levels without sequence information (Pröschold and Darienko 2020). Recently, the genus *Stichococcus* was divided into seven lineages using a multigene approach by Pröschold and Darienko (2020).

According to our ITS2 rRNA sequence-structure based phylogeny, isolate ARK-S10-19 clearly belongs to one of the seven lineages under *Deuterostichococcus* designated as *D. epilithicus* (Pröschold and Darienko 2020). However, ARK-S05-19 formed a distinct and unique clade with the isolate Prasiolales sp. S2RM26 as an undescribed sister group to *Stichococcus* and *Desmococcus*. According to Coleman (2003), the closely related species could be delineated by the presence of CBCs in the conserved region of the ITS2 rRNA secondary structure, as it correlates with the separation of biological species. Although a CBC in helix I of ITS2 rRNA between ARK-S05-19 and S2RM26 was found in a position assumed to be outside
<figure><img src="image_8.png" /><figcaption>**Figure 7.** A phylogenetic tree based on the alignment of 18S rRNA gene sequences using the ML method for *Xanthonema*, the new isolate (sequence) is shown in bold. The best model was TN93 + G + I calculated by MEGAX 10.1.8. Numbers next to branches indicate statistical support values [ML bootstraps (1000 replicates)/Bayesian posterior probabilities]. The bootstrap support values above 70% and Bayesian posterior probabilities above 0.9 are shown. Thick lines indicate the branches with full statistical support (ML/BI:100/1.00). *Pylaiella littoralis* CCAP 1330/3 (AY032606) and *Nannochloropsis limnetica* SAG 18.99 (AF251496) serve as the outgroup.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** Growth patterns of eight microalgal isolates at 2°C and 10°C (Experiment 1). The carrying capacity (maximum cell mass production, K, g·L⁻¹) and the intrinsic growth rate (r, d⁻¹) are shown in the gray box. Mean values (±standard error) of triplicates per isolate are given in the table. The lines show the predicted growth models at 2°C (blue dashed line) and at 10°C (grey line). The coloured points indicate the measured biomass concentrations over 20–45 days (n = 3), where lighter coloured points for each isolate indicates the measurements at 2°C. The asterisks indicate the mean values (±standard error) of duplicates per isolate due to an outlier that did not reach the stationary phase.</figcaption></figure>

the conserved region as defined by Pröschold and Darienko (2020), the two algae were clearly phylogenetically distinct (Fig. 4), suggesting that this clade may contain hidden species diversity. ARK-S05-19 displayed *Stichococcus*-like morphology rather than that of *Desmococcus*, which can be easily distinguished by the unique

characteristics of *Desmococcus* including formation of package-like structures with branched filaments and production of akinetes with ornamented cell walls (Broady and Ingerfeld 1993). Unfortunately, morphological comparison with Prasiolales sp. S2RM26 could not be performed, and ARK-S05-19 is tentatively classified
<figure><img src="image_10.png" /><figcaption>**Figure 9.** Neutral lipid accumulation in *C. reticulata* ARK-S12-19 under nitrogen starvation (Experiment 3). (A) Analysis of TAG accumulation (bottom left) and the average TAG fatty-acid composition (top left) from *n* = 3 replicates in *C. reticulata*. 3t-16 : 1 indicates 3-trans-hexadecenoic acid. TAG content is determined by high performance thin-layer chromatography (HPTLC) together with standards (TAG; triacylglycerol, FFA; free fatty acids, DAG; diacylglycerol). The plate is labelled with Primuline and visualized under 254 nm UV light. (B) Flow-cytometric analysis of cell biovolume-normalized Bodipy® fluorescence over 11 days (*n* = 3). The white vertical lines indicate the median of the frequency and the grey shaded points indicate individual cell data. Multispectral images show the accumulation of LDs in the cells. The Bodipy® fluorescence (green) and chlorophyll autofluorescence (red) are shown. R, region; R2; round shaped cells, R4; oval shaped cells.</figcaption></figure>

as Stichococcaceae sp. Further morphological comparison with closely related strains and sequencing of 18S rRNA gene and *rbcL* genes within the clade could in theory facilitate establishment of new species and improve our understanding of the biogeography and ecophysiology of Stichococcaceae and their appearance in melting snowfields worldwide.

The most commonly found snow algae are typically members of the *Chloromonadinia* clade within the Volvocales (Barcytė et al. 2018a, Hoham and Remias 2020). The *Chloromonadinia* originally comprised solely the genus *Chloromonas* based on *C. reticulata* (Gobi) as the type species, but further investigation led to the discovery of new clades, and the genera *Chlainomonas* (Christen 1959, Novis et al. 2008), *Gloeomonas* (Klebs 1886, Nozaki et al. 2010), *Ostravamonas* (Hoham and Remias 2020, Susanti et al. 2020), and *Ixipapillifera* (Nakada et al. 2016) are currently placed within the same phylogroup. Our sequence analyses show that ARK-S11-19 and ARK-S12-19 are identical to the type species *C. reticulata* Gobi CCAP 11/110 (= CCCryo 213-05; SAG 29.83; UTEX 1970) placed in the *Reticulata* group (Matsuzaki et al. 2012), within the core *Chloromonas* (part of clade 1) that were suggested to be a true *Chloromonas* genus by Barcytė et al. (2018a). Both isolates also exhibited the key morphological features of *C. reticulata* which distinguishes them from three closely related species in the *Reticulata* clade (Matsuzaki et al. 2012). Whilst *C. reticulata* was previously evaluated as a robust cryotolerant mesophile (CCALA 970 in Lukeš et al. 2014 was isolated from red snow and belongs to *C. reticulata*, Barcytė et al. 2018b), it likely does not form snow-blooms, though it is commonly found in snow samples using culture-dependent techniques (Novakovskaya et al. 2018). In contrast, clade 2 encompasses psychrophilic microorganisms forming snow blooms, supported by Sanger sequences from field-collected cells in blooms that are conspecific with those obtained from laboratory strains (Matsuzaki et al. 2019). Still, *C. reticulata* has been isolated from

multiple continents; CCAP 11/110 (= CCCryo 213-05; SAG 29.83; UTEX 1970) from snow in the USA, *C. reticulata* SAG 32.86 from moss in Antarctica, CCCryo 338-08 from snow in France, SAG 26.90 a Norwegian lake and SYKOA Ch-054-11 from snow in Russia (Fig. S7). Our analysis further supports the ecological versatility and cosmopolitan distribution of *C. reticulata* in cold habitats, including snow in northern Norway.

On the other hand, the 18S rRNA gene sequence designated ARK-S08-19 as *Chloromonas* sp. and placed it within a clade containing *C. augustae*, yet its partial 5.8S rDNA + ITS2 rDNA + partial 28S rDNA had only one BLASTn hit to *C. augustae* SAG 5.73 with very low query cover (40%). Such low coverage of ITS2 rDNA sequences together with relatively low similarity of the 18S rRNA gene suggests that ARK-S08-19 may be an independent novel species. Yet, morphologically, the vegetative cells resemble those of *C. augustae* (Skuja 1943, Pröschold et al. 2001). *Chloromonas* augustae has been isolated from freshwater habitats including paddy soil, swamp, and salt marsh pools (Pröschold et al. 2001), though there are no previous reports from melting snow. Sequencing the ITS2 rDNA of related strains and the use of additional markers such as *rbcL* genes could provide more robust species-level identification, whilst ultrastructural analysis may also aid characterization of ARK-S08-19.

Lastly, the genus *Xanthonema* comprise freshwater algae that have been found in soil, snow and freshwaters from alpine, polar, and temperate environments (Broady et al. 1997, Andersen 2004, Broady 2005). The most recent taxonomic assessment of Tribonematalean algae was based on *rbcL* and 18S rRNA gene sequences, where *Xanthonema* clade 1 was defined as the genus *Xanthonema*, and clade 2 and 3 were defined as new genera awaiting taxonomic revision (Maistro et al. 2009, Rybalka et al. 2009). Our analyses showed that ARK-S13-19 is *X. bristolianum* and identical to CCALA 516, placed within *Xanthonema* clade 1 (Rybalka et al. 2009). Molec-
ular evidence was supported by the cell and chloroplast mor-
phologies of ARK-S13-19, which resembled that of X. bristolianum
CCALA 516 as described by Pascher (1939), although we occasion-
ally observed eyespots in some cells that were not mentioned in
the original description. Considering the reports of Xanthonema
in colder habitats including soil and snow (Elster et al. 1999, Ry-
balka et al. 2009, Schulz et al. 2016, Borchhardt et al. 2017a,b), X.
bristolianum is probably widely dispersed across the globe, and its
absence in snow algae communities in previous metabarcoding
studies (Segawa et al. 2018, Lutz et al. 2019) suggests that Xan-
thonema is an opportunistic soil alga rather than a typical snow-
inhabiting species.

A limitation of our methodology is that the most abundant red-
pigmented cells in the snow samples may not be readily cultur-
able. Researchers have reported only limited success in recover-
ing carotenoid-rich encysted red stages into green, viable cultures
(e.g. Procházková et al. 2018a,b). However, Raymond et al. (2022)
showed that orange-coloured cells recognized as Sanguina auran-
tia in field samples were conspecific with their isolates of green
biciated cells that are characteristic of Volvocine algae. Their cul-
tures could only be kept in liquid medium, implying that related
Sanguina sp. and presumably similar algae causing red snow may
have specific cultivation requirements. Similarly, Procházková et
al. (2019b) obtained cultivable green cells, which were conspe-
cific with field collected orange cysts of Chloromonas hindakii. The
authors allowed the cysts to recover and produce a motile, cili-
ated cell culture in deionized water under laboratory conditions
by simulating the diurnal temperature regime of melting snow.
Chloromonas hindakii has been reported to propagate faster in low
nutrient medium (0.6 N BBM) (CCCryo culture collection), and the
use of standard media in our study may have selected strains
that are better adapted to high-nutrient conditions. Expanding the
range of isolation and cultivation conditions could potentially in-
crease the number of unique strains obtained, whilst comparing
isolates to whole community level metabarcode (amplicon) se-
quences from the same red snow microbiomes may help to estab-
lish the relative contribution of our isolates to in-situ snow algae
communities.

Snow-algae dispersal, community construction,
and physiological traits

Snow algae communities in northern Europe have not been exten-
sively studied and until recently little was known about their dis-
tribution and how their constituent species compared with visibly
similar red-snow microbiomes around the world. Our isolates in-
clude cosmopolitan species together with new, unique strains that
are potentially endemic to the region. Assuming that microbial
dispersal is not a major barrier to initial colonization, what shapes
snow-algae microbiomes in different locations? Recent data im-
plicate a multitude of factors driving variation in emergent com-
munities that, for example, vary with altitude in the same region
(Yakimovich et al. 2020, Stewart et al. 2021). The interplay be-
tween temperature, underlying soil structure, and snowmelt dy-
namics thus likely mediate the annual propagation of alternative
species and community structures at different sites. Especially,
the growth rate of microbes is a fundamental life history trait that
determines their ecological success, the outcome of competitive
interactions and ultimately the composition of natural microbial
communities and their biogeochemical functions (Pianka 1970,
Dorodnikov et al. 2009, Brown et al. 2016). At 10°C the growth rates
of our isolates varied two-fold from 0.45 to 0.86 d⁻¹. Chloromonas
reticulata ARK-S11-19 and ARK-S12-19 presented the highest spe-

cific growth rates, whilst novel *Chloromonas* sp. ARK-S08-19 was
characterized by substantially lower growth rates and carrying
capacity, reflecting acute competitive differences between related
strains. At lower 2°C temperatures close to that of melting snow,
the growth rates of all strains except R. *nivale* ARK-S01-19 were
suppressed, and differences between strains were more conspic-
uous. The higher growth rate of R. *nivale* strain ARK-S01-19 com-
pared with ARK-S02-19 at 2°C highlights how ecotypic variation
and selection on thermal tolerance may reshape snow algae com-
munities. Although they are the same species, the closest rela-
tive of ARK-S01-19 was isolated from Svalbard (77°N), suggest-
ing that ARK-S01-19 may have acquired physiological adaptation
to lower temperatures compared with ARK-S02-19. *Chloromonas*
isolates ARK-S08-19, ARK-S11-19, ARK-S12-19 exhibited strong
growth rates (0.34–0.53 d⁻¹) and reached consistent, high cell den-
sities at 2°C, supporting the widespread global distribution of this
genus in melting snow. The ephemeral and dynamic nature of
melting snowfields implies that opportunistic *r*-strategists with
higher growth rates, and possibly high motility, should prevail in
these conditions. Nevertheless, snow and ice are multi-stress envi-
ronments characterized by high light, freeze-thaw cycles, and lim-
ited nutrients that may impose additional selective constraints on
the microbes capable of surviving there that may explain the co-
existence and succession of diverse species (Blagodatskaya et al.
2007, Krug et al. 2020, Tucker and Brown 2022). Further, in-situ
growth rates and the recruitment of snow microbial communities
from the soil are probably further curbed by access to light and
nutrients.

Our data show that fatty acid profiles are species and culture age specific. Under standard high-nutrient conditions (Experiment 2) the Chlorophyta *Raphidonema*, Stichococcaceae, and *Chloromonas* were characterized by abundant n-3 PUFAs, predominantly C18:3n-3, with increases in C18:1n-9 in the stationary phase. C18:3n-3 is commonly found in plastidial membrane lipids monogalactosyldiacylglycerol, digalactosyldiacylglycerol, sulfoquinovosyldiacylglycerol, and phosphatidyl-glycerol (Li-Beisson et al. 2015), whilst C18:1n-9 is known to accumulate in TAG under growth limiting conditions in Chlorophyta (Solovchenko 2012, Andeden et al. 2021, Wan Afifudeen et al. 2021). C18:1n-9 has been reported as a dominant fatty acid in many isolated snow algae strains in N-deprived or stationary phase conditions, including *Raphidonema* and *Chloromonas* (Spijkerman et al. 2012, Hulatt et al. 2017, Suzuki et al. 2019), as well as in field samples of red snow composed of *Sanguina nivaloides* (Procházková et al. 2019a) and other red, orange, or green snow (Bidigare et al. 1993, Spijkerman et al. 2012, Davey et al. 2019). On the other hand, monospecific snow algae blooms caused by C. *kaweckae* (green, Procházková et al. 2022), C. *nivalis* subsp.tatrae, subsp. (brownish-red, Procházková et al. 2018a), *Chloromonas krienitzii* (orange, Procházková et al. 2020), C. *hindakii* (orange, Procházková et al. 2019b), *Chlainomonas* sp. (red, Procházková et al. 2018b) contained largely C18:3n-3 fatty acids rather than C18:1n-9. The fatty acid C16:4n-3 was also a major feature of genus *Chloromonas* and secondly *Raphidonema*, which is found exclusively in plastidal membrane monogalactosyldiacylglycerol in *Chlamydomonas reinhardtii*, and occasionally present in the Trebouxiophyceae (Lang et al. 2011). Both *Raphidonema* strains produced moderate proportions of the long-chain PUFA C20:5n-3, which has previously been described amongst related strains from Svalbard (Spijkerman et al. 2012) and the Antarctic (Suzuki et al. 2019).

The fatty acid composition of X. *bristolianum* ARK-S13-19 (Ochrophyta) differed substantially from that of the Chlorophyta. X. *bristolianum* was characterized by a large proportion of C20:5n-3
under high-nutrient conditions, and by large amounts of C16:1n-7 in the stationary phase. The latter FA is not generally abundant in the Chlorophyta and has been detected only at low concentrations in red, orange, and green snow algal blooms in the field (Spijkerman et al. 2012, Davey et al. 2019). X. *bristolianum* was also characterized by multiple omega-4 fatty acids, including C16:2n-4 and C16:3n-4 which have previously been reported within the kingdom Chromista including the related species *Tribonema ultriculsum* (Chen et al. 2021), diatoms including Phaeodactylum tricor nutum and Chaetoceros (Bacillariophyta) (Miller et al. 2014, Qiao et al. 2016), and Isochrysis aff. *galbana* (Haptophyta) (Huerlimann et al. 2014). Both C16:2n-4 and C16:3n-4 are reportedly located pri- marily in the thylakoid membrane lipid monogalactodiacylglyc- erol (Remize et al. 2020). However, to our knowledge the C18:3n-4, C18:4n-4, NMI-18:3 and NMI-18:4 fatty-acids found in ARK-S13- 19 do not seem to have been reported in environmental samples nor other microalgae, and further lipidomic and genomic investi- gation may elucidate their biosynthetic pathways and functional roles.

Lipids can provide information on the physiological status of snow algae communities (Lutz et al. 2015). Davey et al. (2019) showed that red snow in Antarctica could be clearly distinguished from green snow by the elevated neutral lipids (TAGs) in red snow communities. TAGs are particularly relevant in snow-algae be- cause LDs are produced under stress, especially the low-nutrient conditions often found in water films infiltrating melting snow crystals (Spijkerman et al. 2012). The fatty acid composition of TAG varies between species, reflecting the availability of alterna- tive fatty acids and the activity of a range of acyltransferase en- zymes (Xin et al. 2019). TAG in C. *reticulata* ARK-S12-19 under nitro- gen starvation contained comparatively less saturated C16:0 fatty acids than related model alga C. *reinhardtii*, although the PUFA contents of C. *reticulata* (26%) was more comparable. PUFA-rich TAG has been reported in other algae and hypothesized to provide a reservoir for rapid reconstruction of membrane lipids, enabling cells to recover efficiently from environmental changes (Bigogno et al. 2002), although metabolic validation is needed. IFC analy- sis provided sensitive detection of the onset of LD development in C. *reticulata* and quantitatively supported neutral lipid/TAG inside single cells. Because the most abundant Chlamydomon- adalean snow algae comprise multiple cell types, IFC-based anal- ysis could be further extended for rapidly characterizing cell- cell heterogeneity and the ecophysiology of communities in nature.

Conclusion

In this study, eight algae strains were isolated and characterized
from red snow in northern Norway, adjacent to the Arctic Circle.
Using molecular and morphological analyses the isolates were
assigned to the species R. *nivale*, D. *epilithicus*, C. *reticulata*, and
X. *bristolianum* as well as potentially novel species belonging to
the Stichococcaceae family and *Chloromonas*. Cultivation revealed
variation in growth rates and lipids, ecologically important traits
that shape snow algae communities and their physiological prop-
erties in the natural environment. The accumulation of TAG in
LDs was also observed in cosmopolitan C. *reticulata*, illustrating
the adaptive metabolic changes that snow algae undergo to sur-
vive in harsh environments. The new cultures will provide further
opportunities to explore the stress ecophysiology and molecular
biology of snow algae in the laboratory.

Author Contributions

Hirono Suzuki (Conceptualization, Data curation, Formal analy-
sis, Investigation, Methodology, Validation, Visualization, Writing –
original draft, Writing – review & editing), Alexandre Détain (Data
curation, Formal analysis, Investigation, Writing – review & edit-
ing), Youngjin Park (Data curation, Formal analysis, Methodology,
Visualization, Writing – review & editing), Viswanath Kiron (Re-
sources, Supervision, Writing – review & editing), René H. Wijf-
fels (Resources, Supervision, Writing – review & editing), Nathalie
Leborgne-Castel (Funding acquisition, Supervision, Writing – re-
view & editing), Lenka Procházková (Formal analysis, Methodol-
ogy, Supervision, Validation, Visualization, Writing – review & edit-
ing), and Chris J. Hulatt (Conceptualization, Data curation, Formal
analysis, Funding acquisition, Investigation, Methodology, Project
administration, Resources, Supervision, Validation, Visualization,
Writing – original draft, Writing – review & editing).

Acknowledgements

We thank Anjana Palihawadana from Nord University for his sup-
port in fatty acids analysis. We appreciate support from Priscil-
lia Gourvil from Roscoff culture collection for 18S rRNA gene se-
quencing, and Dr. Reijo Käkelä as well as Dr. Hanna Ruhanen
from Helsinki University Lipidomics Unit (HiLIPID), Helsinki Insti-
tute of Life Science (HiLIFE), and Biocenter Finland for perform-
ing EI-mass spectrometry to verify FAME structures. We thank Dr.
Thomas Leya for providing insight into the morphological obser-
vation of our isolates and curation at CCCryo.

**Supplementary Data**

Supplementary data is available at FEMSEConline.

**Conflict of Interest statement.** The authors do not have any competing interests in the outcomes of the study.

Funding

This work was supported by the Norwegian Research Council (Au-
rora project [No. 319352]) to H.S., A.D., N.L.C., and C.J.H., a Euro-
pean Commission Marie Skłodowska-Curie individual fellowship
[No. 749910] to C.J.H., the Charles University Research Centre pro-
gram [No. 204069], the INTER-EXCELLENCE project [grant number
LTAIN19139], and the institutional long-term research plan of the
Institute of Botany of Czech Academy of Sciences, [RVO 67985939]
to LP.

**References**

Altschul SF, Gish W, Miller W et al. Basic local alignment search tool.
*J Mol Biol* 1990;**215**:403-10.

Andeden EE, Ozturk S, Aslim B. Effect of alkaline pH and nitrogen starvation on the triacylglycerol (TAG) content, growth, biochemical composition, and fatty acid profile of *auxenochlorella protothecoides* KP7. *J Appl Phycol* 2021;**33**:211-25.

Andersen RA. Biology and systematics of heterokont and haptophyte
algae. *Am J Bot* 2004;**91**:1508-22.

Barcyté D, Hodač L, Nedbalová L et al. *Chloromonas arctica* sp. nov., a psychrotolerant alga from snow in the High Arctic (Chlamydomonadales, Chlorophyta). *Int J Syst Evol Microbiol* 2018b;**68**: 851-9.
Barcytè D, Hodač L, Nedbalová L et al. *Chloromonas svalbardensis* n. sp. with insights into the phylogroup *Chloromonadinia* (Chlorophyceae). *J Eukaryot Microbiol* 2018a;**65**:882-92.

Bidigare RR, Ondrusek ME, Kennicutt MC et al. Evidence a ph toprotective for secondary carotenoids of snow algae. *J Phycol* 1993;**29**:427-34.

Bigogno C, Khozin-Goldberg I, Cohen Z. Accumulation of arachidonic acid-rich triacylglycerols in the microalga *Parietochloris incisa* (Trebuxiophyceae, Chlorophyta). *Phytochemistry* 2002;**60**:135-43.

Bischoff H. Phycological Studies IV. Some soil algae from enchanted rock and related algal species. *U. Texas Pub.* 1963;**6318**:1-95.

Blagodatskaya EV, Blagodatsky SA, Anderson TH et al. Priming effects in Chernozem induced by glucose and N in relation to microbial growth strategies. *Appl. Soil Ecol.* 2007;**37**:95-105.

Borchhardt N, Baum C, Mikhailuk T et al. Biological soil crusts of Arctic Svalbard-water availability as potential controlling factor for microalgal biodiversity. *Front Microbiol* 2017a;**8**:1485-.

Borchhardt N, Schiefelbein U, Abarca N et al. Diversity of algae and lichens in biological soil crusts of Ardley and King George islands, Antarctica. *Antarct Sci* 2017b;**29**:229-37.

Breuer G, Lamers PP, Martens DE et al. the impact of nitrogen star vation on the dynamics of triacylglycerol accumulation in nine microalgae strains. *Bioresour Technol* 2012;**124**:217-26.

Britton G, Liaaen-Jensen S, Pfander H. *Carotenoids, Vol. 4: Natural Functions*, Birkhäuser Basel. Basel, Switzerland. 2008.

Broady PA, Ingerfeld M. Three new species and a new record of chaetophoracean (Chlorophyta) algae from terrestrial habitats in Antarctica. *Eur J Phycol* 1993;**28**:25-31.

Broady PA, Ohtani S, Ingerfeld M. a comparison of strains of *xanthomonema* (= *Heterothrix*, Tribonematales, Xanthophyceae) from Antarctica, Europe, and New Zealand. *Phycologia* 1997;**36**:164-71.

Broady PA. The distribution of terrestrial and hydro-terrestrial algal associations at three contrasting locations in southern Victoria Land, Antarctica. *Algol Stud* 2005;**118**:95-112.

Brown SP, Ungerer MC, Jumpponen A. A community of clones: snow algae are diverse communities of spatially structured clones. *Int J Plant Sci* 2016;**177**:432-9.

Butcher R. Contributions to our knowledge of the smaller marine algaes. *J Mar Biol Ass* 1952;**31**:175-91.

Caisová L, Marin B, Melkonian M. A consensus secondary structure of ITS2 in the Chlorophyta identified by phylogenetic reconstruction. *Protist* 2013;**164**:482-96.

Chen W, Luo L, Han D et al. Effect of dietary supplementation with filamentous microalga *Tribonema ultriculosum* on growth performance, fillet quality and immunity of rainbow trout *oncorhynchus mykiss*. *Aquacult Nutr* 2021;**27**:1232-43.

Christen H. Flagellaten aus dem Schützenweiher bei Veltheim. *Mitteilungen Der Naturwissenschaftlichen Gesellschaft in Winterthur* 1959;**29**:167-89.

Christie WW. LIPID MAPS Lipidomics Gateway, the LipidWeb, methyl esters of fatty acids, archive of mass spectra. 2021. https://www.lipidmaps.org/resources/lipidweb/index.php?pag e=ms/methesters/me-arch/index.htm.

Coleman AW. ITS2 is a double-edged tool for eukaryote evolutionary comparisons. *Trends Genet* 2003;**19**:370-5.

Coleman AW. Pan-eukaryote ITS2 homologies revealed by RNA secondary structure. *Nucleic Acids Res* 2007;**35**:3322-9.

Cook JM, Tedstone AJ, Williamson C et al. Glacier algae accelerate melt rates on the south-western Greenland ice sheet. *The Cryosphere* 2020;**14**:309-30.

Darty K, Denise A, Ponty Y. VARNA: interactive drawing and editing of the RNA secondary structure. *Bioinformatics* 2009;**25**:1974.

Davey MP, Norman L, Sterk P et al. Snow algae communities in Antarctica: metabolic and taxonomic composition. *New Phytol* 2019;**222**:1242-55.

Dorodnikov M, Blagodatskaya E, Blagodatsky S et al. Stimulation of r - vs. K -selected microorganisms by elevated atmospheric CO2 depends on soil aggregate size. *FEMS Microbiol Ecol* 2009;**69**:43-52.

Elster J, Lukesová A, Svoboda J et al. Diversity and abundance of soil algae in the polar desert, Sverdrup Pass, central Ellesmere Island. *Polar Rec* 1999;**35**:231-54.

Engstrom CB, Yakimovich KM, Quarmby LM. Variation in snow alga e blooms in the coast range of British Columbia. *Front Microbiol* 2020;**11**:569.

Goold H, Beisson F, Peltier G et al. Microalgal lipid droplets: composition, diversity, biogenesis and functions. *Plant Cell Rep* 2015;**34**:545-55.

Gray A, Krolikowski M, Fretwell P et al. Remote sensing reveals Antarctic green snow algae as important terrestrial carbon sink. *Nat Commun* 2020;**11**:1-9.

Hendriks L, Goris A, Neefs J-M et al. the nucleotide sequence of the small ribosomal subunit RNA of the yeast *Candida albicans* and the evolutionary position of the fungi among the eukaryotes. *Syst Appl Microbiol* 1989;**12**:223-9.

Hillebrand H, Dürselen CD, Kirschtel D et al. Biovolume calculation for pelagic and benthic microalgae. *J Phycol* 1999;**35**:403-24.

Hodač L, Hallmann C, Spitzer K et al. Widespread green algae *Chlorella* and *Stichococcus* exhibit polar-temperate and tropical-temperate biogeography. *FEMS Microbiol Ecol* 2016;**92(8)**:fiw122.

Hoham RW, Bonome TA, Martin CW et al. A combined 18S rDNA and *rbcL* phylogenetic analysis of *Chloromonas* and *Chlamydomonas* (Chlorophyceae, Volvocales) emphasizing snow and other cold- temperature habitats. *J Phycol* 2002;**38**:1051-64.

Hoham RW, Remias D. Snow and glacial algae: a review. *J Phycol* 2020;**56**:264-82.

Hoham RW. Pleiomorphism in the snow alga, *Raphidonema nivale* Lagerh. (Chlorophyta), and a revision of the genus *Raphidonema* Lagerh. *Syesis* 1973;**6**:255-63.

Huerlimann R, Steinig EJ, Loxton H et al. Effects of growth phase and nitrogen starvation on expression of fatty acid desaturases and fatty acid composition of *Isochrysis* aff. *Galbana* (TISO). *Gene* 2014;**545**:36-44.

Hulatt CJ, Berecz O, Egeland ES et al. Polar snow algae as a valuable source of lipids? *Bioresour Technol* 2017;**235**:338-47.

Keller A, Schleicher T, Schultz J et al. 8S-28S rRNA interaction and HMM-based ITS2 annotation. *Gene* 2009;**430**:50-7.

Klebs G. Über die Organisation der Gallerte bei Einigen Al gen und Flagellaten. *Untersuchungen Botanische Institut Tübingen* 1886;2:333-417, pls III, IV.

Kol E. on the red snow of Finse (Norway). *Ann Hist Musei Natl Hungarici* 1963;**55**:155-60.

Kol E. Trochiscia (Chlorophyta) red snow from Swedish Lapland. *Ann Hist Musei Natl Hungarici* 1974;**66**:59-63.

Krug L, Erlacher A, Markut K et al. the microbiome of alpine snow algae shows a specific inter-kingdom connectivity and algae-bacteria interactions with supportive capacities. *ISME J* 2020;**14**:2197-210.

Kumar S, Stecher G, Li M et al. MEGA X: molecular evolution-ary genetics analysis across computing platforms. *Mol Biol Evol* 2018;**35**:1547.

Lang I, Hodac L, Friedl T et al. Fatty acid profiles and their distribu- tion patterns in microalgae: a comprehensive analysis of more than 2000 strains from the SAG culture collection. *BMC Plant Biol* 2011;**11**:124.
Leya T, Rahn A, Lütz C et al. Response of Arctic snow and permafrost algae to high light and nitrogen stress by changes in pigment composition and applied aspects for biotechnology. *FEMS Microbiol Ecol* 2009;**67**:432–43.

Leya T. Snow algae: adaptation strategies to survive on snow and ice. In: ,Seckbach J Oren A, Stan-Lotter H (eds.). *Polyextremophiles: Life under Multiple Forms of Stress*, vol 27. 401–23. Springer Netherlands, Dordrecht. 2013.

Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. *Bioinformatics* 2006;**22**:1658–9.

Li-Beisson Y, Beisson F, Riekhof W. Metabolism of acyl-lipids in *Chlamydomonas reinhardtii*. *Plant J* 2015;**82**:504–22.

Lukeš M, LProcházková, Shmidt V, Nedbalová L, Kaftan D. Temperature dependence of photosynthesis and thylakoid lipid composition in the red snow alga *Chlamydomonas* cf. *nivalis* (Chlorophyceae). *FEMS Microbiology Ecology* 2014;**89**:2 303–315.

Lupette J, Marechal E. The puzzling conservation and diversification of lipid droplets from bacteria to eukaryotes. In: M. Kloc. *Symbiosis: Cellular, Molecular, Medical and Evolutionary Aspects*, vol 69. 281–334. Springer, Cham: Switzerland. 2020.

Lutz S, Anesio AM, Field K et al. Integrated 'Omics', targeted metabolite and single-cell analyses of Arctic snow algae functionality and adaptability. *Front Microbiol* 2015;**6**:1323

Lutz S, Anesio AM, Raiswell R et al. the biogeography of red snow microbiomes and their role in melting Arctic glaciers. *Nat Commun* 2016;**7**:11968.

Lutz S, Procházková L, Benning LG et al. Evaluating high-throughput sequencing data of microalgae living in melting snow: improvements and limitations. *Fottea* 2019;**19**:115.

Magoč T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. *Bioinformatics* 2011;**27**:2957–63.

Maistro S, Broady PA, Andreoli C et al. Phylogeny and taxonomy of xanthophyceae (Stramenopiles, Chromalveolata). *Protist* 2009;**160**:412–26.

Matsuzaki R, Hara Y, Nozaki H. A taxonomic revision of *Chloromonas reticulata* (Volvocales, Chlorophyceae), the type species of the genus *Chloromonas*, based on multigene phylogeny and comparative light and electron microscopy. *Phycologia* 2012;**51**:74–85.

Matsuzaki R, Nozaki H, Takeuchi N et al. Taxonomic re-examination of "*Chloromonas* nivalis (Volvocales, Chlorophyceae) zygotes" from Japan and description of *C. muramotoi* sp. nov. *PLoS One* 2019;**14**:e0210986.

Medlin L, Elwood HJ, Stickel S et al. the characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. *Gene* 1988;**71**:491–9.

Miller MR, Quek SY, Staehler K et al. Changes in oil content, lipid class and fatty acid composition of the microalga *Chaetoceros calcitrans* over different phases of batch culture. *Aquac Res* 2014;**45**:1634–47.

Nakada T, Tomita M, Wu JT et al. Taxonomic revision of *chlamydomonas* subg. *Amphichloris* (Volvocales, Chlorophyceae), with restoration of the genus *Dangeardinia* and descriptions of *ixipapillifera* gen. nov. and *rhysamphichloris* gen. nov. *J Phycol* 2016;**52**:283–304.

Novakovskaya IV, Patova EN, Boldina ON et al. Molecular phylogenetic analyses, ecology and morphological characteristics of *Chloromonas reticulata* (Goroschankin) Gobi which causes red blooming of snow in the subpolar Urals. *Cryptogamie, Algologie* 2018;**39**:199–213.

Novis PM, Hoham RW, Beer T et al. Two snow species of the quadriflagellate green alga *chlainomonas* (Chlorophyta, Volvocales): ultrastructure and phylogenetic position within the *Chloromonas* clade. *J Phycol* 2008;**44**:1001–12.

Novis PM, Podolyan A, Kodner R. New Zealand isolates from snow of the widespread algal genus *Raphidonema* assigned to a single species. *N Z J Bot* 2023;**61**:1–12.

Nozaki H, Nakada T, Watanabe S. Evolutionary origin of *Gloeomonas* (Volvocales, Chlorophyceae), based on ultrastructure of chloroplasts and molecular phylogeny *J Phycol* 2010;**46**:195–201.

Pascher W. *Heterokonten*. Akademische Verlagsgesellschaft MBH, Leipzig. 1939.

Pianka ER. On r-and K-selection. *Am Nat* 1970;**104**:592–7.

Probert I, Siano R, Poirier C et al. *Brandtodinium* gen. nov. and B. nutricula comb. nov. (Dinophyceae), a dinoflagellate commonly found in symbiosis with polycystine radiolarians. *J Phycol* 2014;**50**: 388–99.

Procházková L, Leya T, Křížková H et al. *Sanguina nivaloides* and *Sanguina aurantia* gen. Et spp. nov. (Chlorophyta): the taxonomy, phylogeny, biogeography and ecology of two newly recognised algae causing red and orange snow. *FEMS Microbiol Ecol* 2019a;**95**(6).fiz064.

Procházková L, Matsuzaki R, Řezanka T et al. the snow alga *Chloromonas kaweekae* sp. nov. (Volvocales, Chlorophyta) Causes green surface blooms in the high tatras (Slovakia) and tolerates high irradiance. *J Phycol* 2022;**59**(1):236–48.

Procházková L, Remias D, Bilger W et al. Cysts of the snow alga *Chloromonas krienitzii* (Chlorophyceae) show increased tolerance to ultraviolet radiation and elevated visible light. *Front Plant Sci* 2020;**11**:14–27.

Procházková L, Remias D, Holzinger A et al. Ecophysiological and morphological comparison of two populations of *Chlainomonas* sp. (Chlorophyta) causing red snow on ice-covered lakes in the High Tatras and Austrian Alps. *Eur J Phycol* 2018b;**53**:230–43.

Procházková L, Remias D, Řezanka T et al. *Chloromonas nivalis* subsp. Tatrae, subsp. nov. (Chlamydomonadales, Chlorophyta): re-examination of a snow alga from the High Tatra Mountains (Slovakia). *Fottea* 2018a;**18**:1–18.

Procházková L, Remias D, Řezanka T et al. Ecophysiology of *Chloromonas* hindakii sp. nov. (Chlorophyceae), causing orange snow blooms at different light conditions. *Microorganisms* 2019b;**7**:434.

Pröschold T, Darienko T. the green puzzle *Stichococcus* (Trebouxiophyceae, Chlorophyta): new generic and species concept among this widely distributed genus. *Phytotaxa* 2020;**441**:113–42.

Pröschold T, Marin B, Schlösser UG et al. Molecular phylogeny and taxonomic revision of *Chlamydomonas* (Chlorophyta). I. Emendation of *Chlamydomonas* Ehrenberg and *Chloromonas* Gobi, and description of *Oogamochlamys* gen. nov. and *Lobochlamys* gen. nov. *Protist* 2001;**152**:265–300.

Qiao H, Cong C, Sun C et al. Effect of culture conditions on growth, fatty acid composition and DHA/EPA ratio of *phaedodactylum tricornutum*. *Aquaculture* 2016;**452**:311–7.

Radakovits R, Jinkerson RE, Fuerstenberg SI et al. Draft genome sequence and genetic transformation of the oleaginous alga *Nanochloropsis gaditana*. *Nat Commun* 2012;**3**:686.

Raymond BB, Engstrom CB, Quarmby LM. the underlying green biciilate morphology of the orange snow alga *sanguina aurantia*. *Curr Biol* 2022;**32**(2):R68–9.

Remize M, Planchon F, Loh AN et al. Study of synthesis pathways of the essential polyunsaturated fatty acid 20:5n-3 in the diatom *Chaetoceros muelleri* using 13C-isotope labeling. *Biomolecules* 2020;**10**:797.

Řezanka T, Nedbalová L, Sigler K et al. Identification of astaxanthin diglucoside diesters from snow alga *Chlamydomonas nivalis* by liquid chromatography-atmospheric pressure chemical ionization mass spectrometry. *Phytochemistry* 2008;**69**:479–90.
Ronquist F, Huelsenbeck JP. MrBayes 3: bayesian phylogenetic inference under mixed models. *Bioinformatics* 2003;**19**:1572-4.

Rybalka N, Andersen RA, Kostikov I et al. Testing for endemism, genotypic diversity and species concepts in Antarctic terrestrial microalgae of the Tribonemataceae (Stramenopiles, Xanthophyceae). *Environ Microbiol* 2009;**11**:554-65.

Schulz K, Mikhailuk T, Dresler M et al. Biological soil crusts from coastal dunes at the Baltic Sea: cyanobacterial and algal biodiversity and related soil properties. *Microb Ecol* 2016;**71**:178-93.

Segawa T, Matsuzaki R, Takeuchi N et al. Bipolar dispersal of red-snow algae. *Nat Commun* 2018;**9**:3094.

Segawa T, Yonezawa T, Matsuzaki R et al. Evolution of snow algae, from cosmopolitans to endemics, revealed by DNA analysis of ancient ice. *ISME J* 2023;1-11.

Skuja H. Ein fall von fakultativer symbiose zwischen operculatem discomycet und einer *chlamydomonade*. Archiv für protistenkunde. *Nova Hedwigia* 1943;**96**:365-76.

Solovchenko A. Physiological role of neutral lipid accumulation in eukaryotic microalgae under stresses. *Russ J Plant Physiol* 2012;**59**:167-76.

Spijkerman E, Wacker A, Weithoff G et al. Elemental and fatty acid composition of snow algae in Arctic habitats. *Front Microbio* 2012;**3**:1-15.

Sprouffske K, Wagner A. Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves. *BMC Bioinf* 2016;**17**:172.

Stewart A, Rioux D, Boyer F et al. Altitudinal zonation of green algae biodiversity in the French Alps. *Front Plant Sci* 2021;**12**:38-54.

Stibal M, Elster J. Growth and morphology variation as a response to changing environmental factors in two Arctic species of *Raphidonema* (Trebouxiophyceae) from snow and soil. *Polar Biol* 2005;**28**:558-67.

Stibal M, Šabacká M, Kaštovská K. Microbial communities on glacier surfaces in Svalbard: impact of physical and chemical properties on abundance and structure of cyanobacteria and algae. *Microb Ecol* 2006;**52**:644.

Susanti H, Yoshida M, Nakayama T et al. a taxonomic reassessment of *Chlamydomonas meslinii* (Volvocales, Chlorophyceae) with a description of *Paludistella* gen. nov. *Phytotaxa* 2020;**432**:65-80.

Suzuki H, Hulatt CJ, Wijffels RH et al. Growth and LC-PUFA production of the cold-adapted microalga *Koliella antarctica* in photobioreactors. *J Appl Phycol* 2019;**31**:981-97.

Tesson SVM, Šantl-Temkiv T. Ice nucleation activity and aeolian dispersal success in airborne and aquatic microalgae. *Front Microbiol* 2018;**9**:2681.

Tucker AE, Brown SP. Sampling a gradient of red snow algae bloom density reveals novel connections between microbial communities and environmental features. *Sci Rep* 2022;**12**:10536.

Wan Affudeen C-L, Loh SH, Aziz A et al. Double-high in palmitic and oleic acids accumulation in a non-model green microalga, *Mesastrum gracile* SE-MC4 under nitrate-repletion and -starvation cultivations. *Sci Rep* 2021;**11**:381.

White TJ, Bruns T, Lee S et al. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: *PCR Protocols: a Guide to Methods and Applications*. London: Academic Press, 1990; 315-22.

Wolf M, Koetschan C, Müller T. ITS2, 18S, 16S or any other RNA—Simply aligning sequences and their individual secondary structures simultaneously by an automatic approach. *Gene* 2014;**546**:145-9.

Xin Y, Shen C, She Y et al. Biosynthesis of triacylglycerol molecules with a tailored PUFA profile in industrial microalgae. *Molecular Plant* 2019;**12**:474-88.

Yakimovich KM, Engstrom CB, Quarmby LM. Alpine snow algae microbiome diversity in the coast range of British Columbia. *Front Microbiol* 2020;**11**:1721.

Yakimovich KM, Gauthier NP, Engstrom CB et al. A molecular analysis of microalgae from around the globe to revise *Raphidonema* (Trebouxiophyceae, Chlorophyta). *J Phycol* 2021;**57**: 1419-32.

Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. *Nucleic Acids Res* 2003;**31**:3406-15.