Skip to main content

Giant viral signatures on the Greenland ice sheet

Abstract

Background

Dark pigmented snow and glacier ice algae on glaciers and ice sheets contribute to accelerating melt. The biological controls on these algae, particularly the role of viruses, remain poorly understood. Giant viruses, classified under the nucleocytoplasmic large DNA viruses (NCLDV) supergroup (phylum Nucleocytoviricota), are diverse and globally distributed. NCLDVs are known to infect eukaryotic cells in marine and freshwater environments, providing a biological control on the algal population in these ecosystems. However, there is very limited information on the diversity and ecosystem function of NCLDVs in terrestrial icy habitats.

Results

In this study, we investigate for the first time giant viruses and their host connections on ice and snow habitats, such as cryoconite, dark ice, ice core, red and green snow, and genomic assemblies of five cultivated Chlorophyta snow algae. Giant virus marker genes were present in almost all samples; the highest abundances were recovered from red snow and the snow algae genomic assemblies, followed by green snow and dark ice. The variety of active algae and protists in these GrIS habitats containing NCLDV marker genes suggests that infection can occur on a range of eukaryotic hosts. Metagenomic data from red and green snow contained evidence of giant virus metagenome-assembled genomes from the orders Imitervirales, Asfuvirales, and Algavirales.

Conclusion

Our study highlights NCLDV family signatures in snow and ice samples from the Greenland ice sheet. Giant virus metagenome-assembled genomes (GVMAGs) were found in red snow samples, and related NCLDV marker genes were identified for the first time in snow algal culture genomic assemblies; implying a relationship between the NCLDVs and snow algae. Metatranscriptomic viral genes also aligned with metagenomic sequences, suggesting that NCLDVs are an active component of the microbial community and are potential “top-down” controls of the eukaryotic algal and protistan members. This study reveals the unprecedented presence of a diverse community of NCLDVs in a variety of glacial habitats dominated by algae.

Introduction

Snow and glacier ice algae thrive on the surface of glaciers and ice sheets worldwide during the summer melt season [1,2,3,4,5,6,7,8,9,10,11,12], producing landscape-wide blooms visible on satellite imagery [13]. Red snow patches on the Greenland ice sheet (GrIS) are dominated by Chloromonas spp. and Chlamydomonas spp. (Chlorophyta), while glacier ice algal blooms are dominated by Ancylonema alaskanum and Ancylonema nordenskioeldii (Streptophyta) species [2, 14]. These algae belong to different taxonomic groups, but they both decrease the surface albedo of the snow and ice, which in turn accelerates melting [15,16,17,18,19,20]. Recent extensive efforts to expand knowledge on the ecology, physiology, and phylogeny of these primary producers’ have so far produced relatively little about their life cycle, including the top-down controls that influence their expansion.

Viruses are abundant and ubiquitous across the whole biosphere [21, 22], including cold [23,24,25,26,27] and polar regions [28,29,30,31,32,33,34,35]. Viruses play an essential role in influencing microbial communities through lysis, metabolic reprogramming, and horizontal gene transfer [36, 37]. The viral shunt within aquatic ecosystems significantly influences the structure of algal blooms and eukaryotic communities [38], thereby impacting local, regional, and global biogeochemical cycles [39] and playing a central role in the termination of marine algal blooms [40, 41]. Most of the nucleocytoplasmic large DNA viruses (NCLDV) investigations typically address marine and freshwater environments [42] and only a few from other environments, including polar regions [43].

NCLDVs (Nucleocytoviricota phylum), also called giant viruses, are a supergroup of double-stranded DNA viruses that infect eukaryotes, possessing large virions (up to 1.2 μm in Pithoviridae [44]) and genome sizes (up to 2.5 Mb in Pandoraviridae [45]). They present a set of signature genes used for phylogenetic analyses but encode genes typical of cellular life, such as tRNA and genes involved in protein biosynthesis [46]. The infection strategies of NCLDVs vary considerably, although similarities in how these viruses enter and exit the host cell can be found [47]. NCLDVs are found as free-living particles in environmental samples, and partial or complete viral genomes have been found to be endogenized in several green algae and other hosts genomes [48, 49]. Nucleocytoviricota was revised after the discovery of unclassifiable families, with the addition of new taxonomic ranks, partitioning them into 6 orders (Chitovirales, Asfuvirales, Pimascovirales, Pandoravirales, Algavirales, and Imitervirales), 32 families, and 344 genera [50]. Recently, taxonomic updates were adopted within the order Imitervirales [50]. This level of viral diversity presents challenges when characterizing environmental samples because of the inherent difficulty of culturing virus-host systems. However, diverse environmental metagenomic studies have emphasized their distribution and diversity, demonstrating their presence in oceans, freshwater, and soil [51, 52], as well as extreme habitats, such as the bathypelagic deep sea ocean [27], marine waters, lakes in Antarctica [53,54,55] and marine waters, cryoconite holes, and an epishelf lake in the Arctic [25, 29, 35].

In this study, we demonstrate through analysis of both metagenomic and metatranscriptomic data that NCLDVs are a key constituent of environmental snow and ice microbial communities from the Greenland ice sheet (GrIS). Habitats analyzed include the following: cryoconite (n = 1), ice core (n = 3), green snow (n = 2), red snow (n = 5), and dark ice (n = 8), including the analysis of one metavirome (< 0.2-µm fraction) from dark ice samples. Furthermore, we assess NCLDVs endogenization within cultivated snow algae genomic assemblies (Chlorophyta). Environmental samples were evaluated for the presence of 10 NCLDV marker genes, encoding for factors for maturation of the viral capsid (MCPs), packaging ATPase (A32), DNA polymerase elongation subunit family B (PolB), D5-like helicase-primase (D5), mRNA-capping enzyme (mRNAc), RNA polymerase large and small subunit (RNApl, RNAps), DNA or RNA helicases of superfamily II (RNR, SFII) and poxvirus late transcription factor VLTF3 like (VLTF3), and their clustering with known viral families. Retrieval of 10 giant virus metagenome-assembled genomes (GVMAGs) that were assigned to the Imitervirales, Asfuvirales, and Algavirales was undertaken for comparison with metagenomic and metatranscriptomic viral genes present in these environmental samples to assess the potential viral influence on snow and glacier ice algal blooms.

Results and discussion

We highlight the unprecedented presence of NCLDV marker genes in microbial communities within Greenland ice sheet surface environments, including cryoconite, dark ice, ice core, and red and green snow, and within the genomic assemblies of five cultivated Chlorophyta snow algae (Fig. 1, Tables S1 and S2).

Fig. 1
figure 1

Greenland 2019 and 2020 sampling campaigns (GrIS19/Mit19 and GrIS20, respectively) for environmental samples. One location on the south side of the Greenland ice sheet (inset 1, bottom left). Three locations on the east side of the Greenland ice sheet: Bruckner Glacier (inset 2, top right), Heim Glacier (inset 3, top right), and Mittivakkat Glacier (inset 4, bottom right). Sample types include the following: cryoconite sediment, ice core, dark surface ice, and green and red snow. Circle sizes indicate the metagenome library’s average coverage depth. (Sample information can be found in Supplementary Table 1)

To reduce false positives of NCLDV marker gene identification, fragmented hits were removed and verified against the NCBI nonredundant (nr) database (April 2023), bolstering the quality of the remaining genes used for phylogenetic comparisons. Inconclusive matches occurred with all tested marker genes, mainly due to the presence of hypothetical proteins generated from poorly annotated bacterial MAGs and unknown endogenous viruses within eukaryotic genomes in the database. This resulted in a total of 879 marker genes, 387 from red snow, 298 from snow algae genomic assemblies, 87 from green snow, and 82 from dark ice (Fig. 2, Figs. S1–S7, Tables S2–S4).

Fig. 2
figure 2

Quality-controlled counts of unfragmented NCLDV marker genes after homology searches against the NCBI nr reference database for each sample of this study. Analysis was carried out in 19 environmental metagenomes (MG) and 18 environmental metatranscriptomes (pooled) obtained from samples of cryoconite (n = 1), ice core (n = 3), green snow (n = 2), red snow (n = 5), dark ice (n = 8)), 1 metavirome (dark ice), and 5 snow algae genomic assemblies from the CCCryo collection. The points represent the total number of each marker gene in the samples with “total” indicating the overall count of marker genes in that sample. The “md” (more depth) notation following selected samples represents those that were re-sequenced with higher metagenomic coverage. Colored symbols on the left of the sample names represent the sample types

To better speculate on the total number of NCLDV’s present, MCP genes were summed and used as a proxy for the number of NCLDVs as they are considered bona fide viral genes [50]. Before stringent quality control measures were applied, there were 211 MCP genes; after quality control, only 67 remained. While this stringent quality control method may limit the detection of novel NCLDVs (the original amount of MCP genes was nearly three times higher), it still emphasizes the potential diversity and abundance of the NCLDVs across the Arctic. The 19 environmental metagenome libraries varied in library coverage and size; however, these values did not relate to the number of NCLDV marker genes identified (Fig. 1, Table 1, Tables S1, S2–S4). NCLDV genes were absent in the three ice cores or in four out of nine dark ice metagenomes (Fig. 2, Table S4).

Table 1 Environmental sample name, sample type, location, filed campaign and year, nucleic acid extraction method (coextraction with PowerLyzer PowerSoil DNA and RNeasy PowerSoil Total RNA kit (denoted with DNA/RNA) or cetyltrimethyl ammonium bromide (denoted with CTAB), or DNA purification resin (denoted with resin)), sequencing platform, and assembly size in base pairs for each environmental sample. The “md” (more depth) notation following selected samples represents those that were re-sequenced with higher metagenomic coverage

The D5, RNApl, and RNAps marker genes were the most abundant in all metagenome samples making up > 50% of the marker genes (Fig. 2, Figs. S1–S3, Tables S2–S3). Similarly, RNApl and RNAps genes in the pooled metatranscriptomes had the most individual counts (35% and 32%, respectively) of the transcribed gene (Tables S2–S3; e-value = 1 × 10−10). However, none of these transcribed sequences matched known NCLDV members on NCBI nr and was excluded from further analyses. In most cases, the signatures of these single NCLDVs marker genes found both in metagenomic and metatranscriptomic data were found in short contigs, impeding a deeper investigation of genomic context. Transcribed MCP genes were the third most abundant (18%) in the metatranscriptomes (Figs. S4 and S8, Tables S2–S4) after RNApl and RNAps counts, confirming high expression in the environment [56]. All the sequences of the transcribed MCP genes were similar to known NCLDV families and clustered closely to the related metagenomic sequence within the phylogenetic tree (Fig. S4). Generally, the marker gene sequences were closely related and clustered on shared phylogenetic tree nodes, despite originating from different metagenomes or genomes. For example, MCP genes from red snow samples, MG12 and MG3, were > 97% identical (Fig. S5) and from an environmental red snow sample (MG28) and Chloromonas remiasii 005–99 or 047–99 were up to 84% identical. Furthermore, the dark ice samples, MG32, MG31, MG19, MG8, and the metavirome, each contained a PolB sequence with a high percentage of identity, > 99.5% (Fig. 3, Table S5).

Fig. 3
figure 3

Maximum-likelihood phylogenetic tree of the NCLDV core gene DNA polymerase (PolB). Sequences recovered from the environmental samples are presented in bold at the tree node. Environmental sample types are specified in correspondence of each sequence. Branches are color-coded by order-level taxonomy. Viral families are specified in the colored ranges. The dark dots at the nodes represent the bootstrap support value of > 70

Red snow samples MG12 and MG3 from the GrIS contained PolB genes that were > 99.4% identical (Fig. 3). The mRNAc genes from dark ice samples, MG32, MG31, MG19, and the metavirome, were > 99.7% identical (Fig. S6). This gene similarity suggests a degree of relatedness with the NCLDVs identified at each environment, despite unique sample types. The eukaryotic diversity and composition in each of the locations are generally composed of the same members (Figs. 4 and S9), which corroborates identifying similar NCLDVs marker genes.

Fig. 4
figure 4

18S rRNA diversity of the 18 environmental samples from TotalRNA. Some eukaryotic phyla are made up of more than one individual. Bacterial phyla are not displayed but make up the empty space above each bar. Calculated relative abundance percentages can be found in Supplementary Table S9

Additionally, D5 marker genes from red snow samples and two of the algae genomic assemblies clustered closely with giant endogenous viral elements (GEVEs) previously found in diverse green algae genomes [48] (Fig. S1). The co-clustering of marker genes found in snow algae genomic assemblies with those from environmental red snow samples (Fig. 3, S1–S7) and GEVEs observed in other green algae (Fig. S1) strongly suggests that chlorophytes may serve as hosts in this environment and have endogenized viral genes.

In total, 10 GVMAGs (Table S6) and 29 individual PolB sequences (Table S4) were retrieved from the 31 different samples investigated, 7 GVMAGs and 22 PolB originating from the environment, and 3 GVMAGs and 7 PolB from the snow algal genome assemblies (Fig. 2, Table S4). Since PolB is the only marker gene typically found as single copy, it is used for phylogenetic placement within known NCLDV families [57]. In red snow and dark ice samples, a few PolB sequences had similar or identical residues (Table S5). Between the 10 GVMAGs, there were 8 unique genome pairs with > 99.1% ANI (Table S7). The origin of these identical GVMAGs were red snow samples from Mittivakkat Glacier (MG12 and MG3) and two Chloromonas remiasii cultures (005–99 and 047–99), further showing giant virus links to red snow algae. However, the GVMAGs retrieved here are not an exhaustive representation of the NCLDVs present in these Greenland environments. There were five MAGs that had less than five NCLDV marker genes, smaller than 100-kbp genome size, and therefore were not considered further as GVMAG. One was from the cryoconite sample (MG30), one from green snow (MG27), two from red snow (MG3_md and MG12), and one from dark ice (MG32). Although these are poor representative, they still indicate potential GVMAG diversity in other habitats in the GrIS. The functional annotation of the 10 GVMAGs highlighted the presence of genes associated with eukaryotic photosynthesis, such as heliorhodopsin, Rubisco LSMT substrate binding, bestrophin chloride channel, and copper amine oxidase [42]. These annotations were found within one Algavirales and three Asfuvirales GVMAGs (MG12_md_6, MG3_12, MG12_md_5, and MG12_2, respectively) (Table S8). These genes are often found endogenized in host genomes, and finding these within GVMAGs from environmental snow samples further indicates a potential host-viral relationship.

Individual phylogenetic trees of marker genes were built to examine the phylogenetic relationships (PolB, Fig. 3) and phylogenetic diversity [50] (D5, RNAps, RNApl, MCP, mRNAc, A32, SFII, VLTF3, and RNR, Figs. S1–S7) between the proteins found in the metagenomic, metatranscriptomic, metaviromic, and genomic contigs in comparison with known viruses. The maximum-likelihood phylogenetic tree of the NCLDV marker gene DNA polymerase (PolB) showed clustering with four known viral families (Allomimiviridae, Pithoviridae, Asfarviridae, Algavirales AG-04), with a clear separation in terms of NCLDV groups based on the sample type (Fig. 3). PolB sequences from red snow samples (three identified in MG12 and one in MG3) and green snow samples (1 sequence from MG27) grouped together with Asfuvirales reference sequences, which is a globally distributed group in the ocean known to infect photosynthetic dinoflagellates, as well as protozoans [58]. The rest of the red snow (MG28, four sequences) and sequences originating from the snow algae genomic assemblies (C. remiasii, 3 sequences, and Microglena cf. sp., 1 sequence) clustered with the Imitervirales, which is the widest order infecting a variety of hosts, including green algae [52, 59]. Signatures found in green snow samples (one sequence retrieved from MG27 and one from MG26) formed a sister group with the Heterosigma akashiwo virus 01 (Algavirales), which has been used as microbiological agent for red tide control in the ocean [60]. Sequences from dark ice were assigned to Pithoviridae (two sequences from MG8, two from MG19, one from MG31, one from MT31, and three from MG32), which mostly infect species of the amoebozoan genus Acanthamoeba [47]. Overall, PolB phylogeny shows a wide diversity of NCLDVs and reveals the potential top-down interactions affecting a diverse eukaryotic host community (algae and protists) on the GrIS.

Different samples of red snow (MG28, MG12, and MG3) contained NCLDV signatures belonging to different families. The concatenated maximum-likelihood tree assigned all the GVMAGs generated from the GrIS2020 red snow sample (MG28), together with the snow algae C. remiasii GVMAGs (Fig. 5), to the family Allomimiviridae, confirming results obtained through PolB phylogeny.

Fig. 5
figure 5

Maximum-likelihood parsimony phylogenetic tree with 1171 external genomes from previously published GVMAGs and 10 GVMAGs from this study. According to the tree, the retrieved GVMAGs cluster in correspondence of the Asfuvirales (3), Imitervirales (6), and Algavirales (1) orders. VGMAGs originated from this study are highlighted by the corresponding environmental sample type symbol. Branches are color-coded by order-level taxonomy. Cultured isolate virus references of interest are labeled in their approximate location along the branches with the following abbreviations: African swine fever virus (ASFV), Heterosigma akashiwo Virus 01 isolate HaV53 (HaV53), Tetraselmis Virus (TetV), Pyramimonas orientalis Virus 01b (PoV-01b), Phaeocystis globosa Virus (PgV)

This family contains the recently cultivated Oceanusvirus kaneohense [61], formerly known as Tetraselmis Virus (TetV), which infects the marine green algae Tetraselmis (Chlorodendrophyceae) [62]. Members of this genus are ubiquitous and commonly found in nutrient-rich marine and fresh waters, although the first TetV-specific host was initially isolated from an oligotrophic habitat (open ocean) [62]. GVMAG_MG28_md_2 and Chloromonas_remiasii_005-99_3 fell within the cluster formed by TetV. Another member of the Allomimiviridae family is the species Heliosvirus raunefjordenense, formerly known as Pyramimonas orientalis Virus 01b (PoV-01b), also infecting chlorophytes [63]. One GVMAG retrieved from green snow (MG27) was assigned to the family IM_18 of the Imitervirales order. This family is represented only by genomes derived from cultivation-independent approaches retrieved from freshwater and marine sources and does not include isolated members at present. One GVMAG originated from red snow (Mittivakkat 2019) was assigned to the order Algavirales (family incertae sedis), an NCLDV order encompassing several well-studied algal viruses [50]. Three GVMAGs originating from red snow sampled from the Mittivakkat Glacier in SE Greenland were assigned to the Asfuvirales family AF_2, and one was assigned to a cluster with uncertain taxonomy. Generally, members of the Asfuvirales infect a mixture of metazoan and protist hosts and are broadly distributed in marine systems [58, 64]. The presence of Pithoviridae signatures in dark ice and their likely associations with protists suggest that the GVMAGs from red snow assigned to the Asfuvirales are also probably associated with protists hosts. The unassigned GVMAG emphasizes the complexity of NCLDV taxonomy, which is constantly growing from metagenomic data but unfortunately lacks additional cultured isolate reference genomes.

The active 18S rRNA eukaryotic community contained algal and protistan members. Dark ice was dominated by the phylum Streptophyta, mainly from the class Zygnematophyceae (7–37% throughout the seven dark ice samples), but also with the presence of chlorophytes, specifically from two classes, Chlorophyceae (5–21%) and Trebouxiophyceae (9–27%) (Fig. 4 and Tables S9–S10). One dark ice sample was used in an attempt to sequence and assemble a draft genome of the Streptophyta glacier ice algae. The final assembly had more than Streptophyta contigs (Table S11), so it was considered as the 19th environmental sample (MG32, Fig. 4) despite different extraction, sequencing, and assembly methods used on the other metagenomic samples. Green and red snow were dominated by algae belonging to the phylum Chlorophyta (20–22% and 28–69%, respectively). The active protistan community included the cercozoa Glissomonadida and ciliate Stokesia, which are commonly found in glaciers, snow, and sea ice [12, 65,66,67]. The order Glissomonadida includes biflagellate gliding bacterivores and algivorous amoeboflagellates [68, 69], which were mainly present in green (4–10%) and red (1–2%) snow samples (Table S6) but also in dark ice (0–1%) and ice core (1%). Stokesia is a large (more than 100 μm) ciliate containing endosymbiotic green algae commonly found within spring phytoplankton blooms in oligo-mesotrophic lakes [70, 71], which was present and active in a dark ice sample from Heim Glacier (MG19, 33%, Table S6). The variety of active algae and protists in these GrIS habitats containing NCLDV marker genes suggests that infection can occur on a range of eukaryotic hosts.

The presence of active members of the community in all sample types was confirmed by the read recruitment analysis showing that reads of the metagenomic samples recruit to the corresponding metatranscriptomic sample. Most of the reads from each sample mainly mapped to their respective assemblies, however also mapped to other environmental types (10 GVMAGS, 23 metagenomes, 1 metavirome, 18 metatranscriptomes (Table S12, Fig. S10)). For example, the Streptophyta-dominated environmental sample (MG32, Fig. S9) mapped 30% of the reads to its own assembly, 4% mapping to red snow or dark ice assemblies, and 8% mapping to the cryoconite assembly (Fig. S10, Table S12). Furthermore, the red snow sample MG28 mainly mapped (27% reads) to a metatranscriptome assembly from green snow (MG27, Fig. S10). This pattern, where one sample type (e.g., red snow) maps at least 5% of the reads to another sample type assembly (e.g., green snow), demonstrates the community overlap between Greenland ice sheet habitat types. This is further seen within the hierarchical clustering groups through the shared read recruitment analysis, where different sample types share similar read recruitment pattern (Fig. S10). The similarities in shared mapping are better underscored by the compared diversity of the 18S rRNA from the metagenomes and metatranscriptomes (Fig. S11). These glacial samples observed diversity was above 325 in most metatranscriptomic samples, except 2 red snow samples (MT12 and MT28, Fig. S11A, Table S1). These two red snow samples, and the deeply sequence MG28, also have a lower observed diversity in 18S rRNA genes from the metagenomes (Fig. S11C). The Shannon index highlighted a high diversity in the samples (Fig. S11A), excluding few exceptions such as the red snow samples MT28, MT12, and MT22 that seemed to harbor a lower diversity (2, 2.9, and 3.3, respectively) and two dark ice samples MT19 and MT31 (2.9, and 2.6, respectively). Overall, the inverse Simpson index showed a low evenness of the samples (ranging between 5.2 and 47) that therefore appeared to be dominated by few taxa (Fig. S11A). The nonmetric multidimensional scaling (NMDS) analysis revealed clustering of the dark ice, ice core, and red snow samples based on the location over a strong association with sample type for both metagenome and metatranscriptomic samples (Figs. S11B and D). Generally, samples also clustered based on sample type, with the exception of the red snow sample MG28, which appeared significantly dissimilar from the others (Fig. S11B).

Overall, metagenomic evidence reveals diverse NCLDV signature genes in snow and ice habitats, highlighting the presence of potential viral controls on these algal communities. Furthermore, the presence of viral genes in Chloromonas spp., Microglena sp., and Sphaerocystis sp. genomic assemblies can be most likely considered a result of past viral DNA integration, as seen already in other non-snow hosted green algae (Chlorophyta) [48] and highlighted by genomic evidence in a comprehensive survey of giant virus DNA integration into genomes of algae and protists [49]. Integration of endogenous viruses in algal genomes is not present in all algal groups and appears to be highly host specific [49]. Nevertheless, the presence of endogenized viruses does impacts the algal genome evolution and potentially the ecological success of these algae [47]. Viruses would not be endogenized in the first place without active viral-host interactions. The co-clustering of metagenomic and endogenized signatures on the phylogenetic trees indicates that these NCDLVs are likely close relatives and allows the host to be inferred. It suggests that the Allomimiviridae group of NCLDV signature genes in these red snow samples is from algal-infecting viruses.

These diverse environmental sample types offer valuable insights into the prevalence of NCLDVs within microbial communities on the Greenland ice sheet. They are primarily associated with snow algae (Chlorophyceae) in red snow, while other signatures, such as Asfarviridae and Pithoviridae, are linked to protists in dark glacier ice algae-dominated habitats. Collectively, these findings suggest that pigmented supraglacial snow and ice habitats contain a diverse array of NCLDVs linked to various eukaryotic hosts. Furthermore, the detection of transcribed viral marker genes that taxonomically identify with NCLDV metagenomic sequences implies an active NCLDV influence on the snow and ice algal community, potentially serving as regulators of colored snow blooms.

Material and methods

Samples collection and preparation

Samples were collected during two fieldwork campaigns in July 2019 and July 2020. Samples in 2019 were collected from three locations in the SE of the GrIS. Mittivakkat Glacier is an independent glacier separated from the GrIS, located on Ammassalik Island, in South-East Greenland, below the Arctic Circle (65.69°N; 37.83°W) (Fig. 1). The samples from Mittivakkat Glacier were collected along a west-sloping transect from the accumulation zone (two red snow samples) to the ablation area (four dark ice samples). The second location was on the GrIS, across the fjord from Mittivakkat Glacier, and one sample of red snow and two dark ice samples were collected from Bruckner and Heim glaciers (65.99°N; 38.44°W, and 65.95°N; 38.53°W, respectively).

In 2020, a set of environmental samples were collected from the GrIS, close to the QAS_U and QAS_M PROMICE stations (61.08°N; 46.83°W and 61.18°N; 46.82°W) in S-Greenland. The samples included ice core (1-m core that included snow and ice transition section, n = 3), cryoconite hole sediment (n = 1), dark ice (n = 1), dark ice sample only for viral fractionation concentration (n = 1) and for the purpose of creating a draft genome of the Ancylonema ice algae (n = 1), green snow algae biofilm (n = 1), green snow (n = 1), and red snow ((n = 2; Fig. 1). Coordinates and details for each sampling site and sample are reported in Table 1 and supplementary material (Table S1). Dark ice represents glacial surface ice that is visually dark as compared to white ice and contains a high abundance (104 cell/ml) of dark pigmented glacier algae, typically dominated by the class Zygnematophyceae [5]. Green and red snow also is visually colored and contains a high abundance of green and red snow algae, both within the class Chlorophyceae [5].

All samples were collected with sterile nitrile gloves and tools and stored in sterile Whirl–Pak® bags. Samples were melted at ambient temperature in the field and filtered through 0.22-μm mixed cellulose ester membrane filters (Sartorius, Germany), which were immediately frozen and transported to the home laboratory in a cryo-shipper at liquid nitrogen temperatures, where they were stored at − 80 °C until further processing. Total DNA and RNA were extracted from the filters with the PowerLyzer PowerSoil DNA isolation kit and the RNeasy PowerSoil Total RNA kit (Qiagen, Germany), respectively, following the manufacturer’s instructions. Nineteen DNA libraries were generated with the NEBNext® Ultra™ II FS DNA Library Prep Kit (Illumina), with 8 rounds of PCR amplification. RNA samples were treated with the DNase Max® kit to remove remaining DNA (Qiagen, Germany) following the manufacturer’s instructions. Eighteen RNA libraries were prepared with the TotalRNA NEBNext® Ultra™ II RNA Library Kit, with 8 rounds of PCR amplification. Sequencing was performed in-house using the NextSeq 500 platform and the 300 cycle v2.5 chemistry (151-bp pair-end reads). The reconstructed, full-length rRNA small subunit (SSU) genes in the 18 environmental transcripts were taxonomically identified with Silva 138.1, BLAST, and CREST4 [72] as part of our in-house TotalRNA workflow (DOI: 10.5281/zenodo.7656004). Chloroplast and mitochondria sequences made up a total of 0.4% to 19.9% of initial sequences and were removed from further analysis. Statistical comparisons of the assembly diversity were analyzed with phyloseq (v 1.44) [73] in R Studio (v 3.17) [74].

Snow algal genomic amplification and assembly

The algal strains Microglena cf. sp. 002b-99, Chloromonas remiasii 047–99 and 005–99, cf. Sphaerocystis sp. 101–99, and Raphidonema sempervirens 011a-99, commonly present on pigmented snow and surface ice of the GrIS and other glaciers, were obtained from the Culture Collection of Cryophilic Algae (CCCryo) at the Fraunhofer IZI-BB Institute (Table 2).

Table 2 Snow algae cultures sequenced and assembled from the CCCryo Culture Collection, all extracted with PowerSoil DNA Isolation Kit for DNA sequencing

They were grown at 10 °C in liquid triple-concentrated Bold’s Basal Medium [75] (pH 5.5) under axenic conditions and continuous illumination as per the CCCryo guidelines. DNA was extracted using the PowerSoil DNA Isolation kit (QIAGEN, Germany) following the manufacturer’s instructions. Each strain was sequenced on a MiSeq flowcell using the 500 cycles v2 chemistry (250-bp pair-end reads) at the Genome Analysis Centre (Earlham Institute, UK). In addition, high-molecular-weight DNA of Sphaerocystis sp. 101–99 was extracted using the QIAGEN genomic-tips extraction kit and sequenced using one PacBio Sequel SMRT Cell (2.0 chemistry) at NERC Biomolecular Analysis Facility — Liverpool. The genomes of all five strains were de novo assembled by the Earlham Institute using the Illumina 250-bp paired-end reads. Quality control of the raw data was done using FastQC (fastqc-0.11.2, http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Preprocessing of the raw reads was done by the Earlham Institute (https://github.com/TGAC/kontaminant) using the pipeline Kontaminant. ABySS (v.1.9.0) [76] was used to perform the de novo assembly of each strain. The PacBio Sequel reads of the strain 101–99 were de novo assembled using Flye (v.2.3.3) [77] using a minimum subread length of 5000 bp and an estimated genome size of 120 Mb. Sample MG32 was extracted with the CTAB (cetyltrimethyl ammonium bromide) method [78] and sequenced on two platforms. Illumina libraries were prepared using the NEBNext® Ultra™ II FS DNA Library Prep Kit (New England Biolabs) and sequenced on a NextSeq 500 instrument with the 300 cycles v2.5 chemistry. Nanopore libraries were prepared using the Ligation Sequencing Kit (LSK-109) and sequenced on a MinION (Oxford Nanopore Technologies, Oxford, UK) with a FLO-MIN106 flow cell, controlled using MinKNOW (19.10.1). Raw nanopore fast5 reads were basecalled with GPU Guppy (3.2.6 + afc8e14). The Illumina reads were quality filtered using trim-galore under default settings (https://github.com/FelixKrueger/TrimGalore). The raw nanopore reads were corrected with the trimmed Illumina reads using LoRDEC [79] with default settings. The corrected long reads were used for de novo whole genome assemblies with Flye [77] under default settings utilizing the − nano-corr flag [77]. under default settings utilizing the − nano-corr flag. The dark ice environmental sample (MG32) containing a high abundance (104 cell/ml) of Ancylonema sp. was taken in an attempt of producing a draft genome of the Streptophyta glacier ice algae. The overall appearance of the sample under the microscope gave the misleading impression that this mixed culture would contain primarily Ancylonema sp. and prokaryotes. The idea was then to remove prokaryotic contigs and have a representative Ancylonema genome. Further analysis on the resulting assembly with BARRNAP (BAsic Rapid Ribosomal RNA Predictor, https://github.com/tseemann/barrnap) revealed the presence of a diverse eukaryotic community. Nevertheless, the sample was kept in the study as it provided another dark ice environmental sample and gave nanopore long reads.

Metagenome, metatranscriptome, and GVMAG assembly

Illumina reads were quality filtered to remove low-quality reads and trimmed with fastp [80] (version 0.20.0) using default options. Trimmed Illumina reads were assembled with metaSPAdes [81] (v3.15.1) specifying the –only-assembler pipeline. Metatranscriptome reads were quality cleaned and trimmed with trim-galore (https://www.bioinformatics.babraham.ac.uk, v0.6.6) using default options. Raw reads were assembled both singularly (each sample) and pooled together (co-assembly) with Trinity assembler [82] (v2.6.6) including the following options: –normalize_by_read_set. Results of the co-assembly are presented as metatranscriptome — pooled.

Giant virus metagenome-assembled genomes (GVMAGs) were created by binning contigs with MetaBAT2 (v2.12.1) [83] using >  = 5000 base-pair contigs. Resulting bins were analyzed for NCLDV marker genes using ViralRecall (v2) [57] and were considered a GVMAG if they had five or more of the marker genes, a genome larger than 100 kbp, and taxonomic placement within other NCLDV genomes [51]. CoverM (v 0.6.1) (https://github.com/wwood/CoverM) was used to assess the read recruitment between all 57 generated assemblies and the environmental sample reads (18 metagenomes, 19 metatranscriptomes, and 1 metavirome). GVMAGs functional annotations were assessed with InterPro [84] and GVOGs [50].

Metavirome construction

Five liters of dark ice from the GrIS 2020 location was prefiltered with 3-μm nitrocellulose membrane filters (Sartorius) to remove large particles and subsequently filtered through 0.2-μm VacuCap™ devices (Pall Corporation), retaining the viral fraction (< 0.2 μm). Viruses were further concentrated from the filtrate using iron chloride flocculation [85] followed by storage at 4 °C. After resuspension in ascorbic-EDTA buffer (0.1-M EDTA, 0.2-M MgCl2, 0.2-M ascorbic acid, pH 6.0), viral particles were concentrated using Amicon Ultra 100-kD centrifugal devices (Millipore) and extracted as previously described [86]. Briefly, viral particle suspensions were treated with Wizard Polymerase Chain Reaction Preps DNA Purification Resin (Promega, Fitchburg, WI, USA) at a ratio of 1-mL sample to 1-mL resin and eluted with TE buffer (10-mM Tris, pH 7.5, 1-mM EDTA) using Wizard Minicolumns. The DNA library was prepared following the NEBNext® Ultra™ II FS DNA Library Prep Kit (Illumina), with 8 rounds of PCR amplification. Sequencing was performed in-house using the NextSeq 500 platform and the 300 cycle v2.5 chemistry (151 -bp pair-end reads). This sample was processed with a small filter size (< 0.2 µm) and treated as the environmental viral fraction. It is important to note that the small filter size will decrease the amount of NCLDV signatures retrieved.

Identification of NCLDVs signatures in metagenomic data

ViralRecall was used to identify NCLDV-like sequences and viral-like regions in all the metagenome, metatranscriptomes, metavirome, and pure algal culture. Options used were as follows: -db marker -c. The “marker” option was used to only search against 10 NCLDV marker genes, encoding for factors for maturation of the viral capsid (MCPs), packaging ATPase (A32), DNA polymerase elongation subunit family B (PolB), D5-like helicase-primase (D5), mRNA-capping enzyme (mRNAc), RNA polymerase large and small subunit (RNApl, RNAps), DNA or RNA helicases of superfamily II (RNR, SFII), and poxvirus late transcription factor VLTF3 like (VLTF3). All resulting hits with an e-value less than e^-10 were used further. These genes are universal NCLDV marker genes and hence are routinely assessed for identification of signatures of NCLDVs in different ecosystems [51]. PolB is the only marker gene typically found as single copy and is therefore used for phylogenetic placement within known NCLDV families [57].

To confirm that virus-like regions belonged to NCLDV families, blastp function against NCBI nr was used, and 50 top hits were verified for each sequence classified as possible NCLDV gene by ViralRecall. A gene was considered from NCLDV when it had NCLDV results within the top 10 hits. The total abundance of the 10 NCLDV core genes in each sample was calculated before and after verification with NCBI nr by summating the marker genes with an e-value cutoff of 1 × 10−10 and normalizing to the total library size. Four of the 19 environmental samples with the highest relative presence of viral marker genes (MG3, MG8, MG12, and MG28; Fig. 2) were chosen to be re-sequenced at a greater depth to provide higher sequencing coverage and increase the chances of assembling GVMAGs.

Phylogeny of unbinned GV marker genes and transcriptomes

MAFFT [87] (v7.475) was used to align the viral regions from sequenced data against the reference sequences using the − auto option to select the appropriate option (L-INS-I, FFT-NS-2, or FFT-NS-i) for each alignment according to the size of input data (options: –maxiterate 1000). Only sequences of marker genes that had an e-value <  = 1 × 10−10 and had a length comparable to the reference sequences (> = 300 aa) were subsequently kept in the tree. Fragmented signatures (< 300 aa) were not included in the phylogenetic placement. For each gene, maximum likelihood phylogenetic trees were built using IQ-TREE [88] v2.0.3. According to BIC scores, LG + F + I + G4 (PolB) was the best model by the “-m TEST” ModelFinder option [89]. IQ-TREE was run with 1000 ultrafast bootstraps (-alrt 1000 -B 1000) to assess confidence [90].

Phylogeny of the GVMAGs against 1171 external Nucleocytoviricota genomes

External Nucleocytoviricota genomes were downloaded from previously published studies [50, 91]. All 1171 external genomes and our 10 GVMAGs were aligned using ncldv_markersearch.py (last update 21Apr2022, github.com/faylward/ncldv_markersearch). A maximum-likelihood phylogenetic tree was constructed using IQ-TREE with the LG + F + I + G4 model with -B ultrafast 1000 bootstraps [90]. Phylogeny assignment was assigned based on previous literature [51].

Availability of data and materials

All metagenomes, metatranscriptomes, and culture genomic assemblies can be found under NCBI BioProject: PRJNA1011216 and BioSamples within. Culture sequenced reads can be found under NCBI BioProject PRJNA1036577.

References

  1. Yallop ML, Anesio AM, Perkins RG, Cook J, Telling J, Fagan D, et al. Photophysiology and albedo-changing potential of the ice algal community on the surface of the Greenland ice sheet. ISME J. 2012;6:2302–13. Available from: https://doi.org/10.1038/ismej.2012.107.

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

    Article  PubMed  PubMed Central  Google Scholar 

  3. Tanaka S, Takeuchi N, Miyairi M, Fujisawa Y, Kadota T, Shirakawa T, et al. Snow algal communities on glaciers in the Suntar-Khayata Mountain Range in eastern Siberia, Russia. Polar Sci. 2016;10:227–38. Available from: https://doi.org/10.1016/j.polar.2016.03.004.

  4. Uetake J, Tanaka S, Hara K, Tanabe Y, Samyn D, Motoyama H, et al. Novel biogenic aggregation of moss gemmae on a disappearing african glacier. PLoS One. 2014;9(11):e112510.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Williamson CJ, Anesio AM, Cook J, Tedstone A, Poniecka E, Holland A, et al. Ice algal bloom development on the surface of the Greenland ice sheet. FEMS Microbiol Ecol. 2018;94(3):fiy025.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lutz S, Anesio AM, Jorge Villar SE, Benning LG. Variations of algal communities cause darkening of a Greenland glacier. FEMS Microbiol Ecol. 2014;89:402–14.

    Article  CAS  PubMed  Google Scholar 

  7. Takeuchi N, Dial R, Kohshima S, Segawa T, Uetake J. Spatial distribution and abundance of red snow algae on the Harding Icefield, Alaska derived from a satellite image. Geophys Res Lett. 2006;33.10. https://doi.org/10.1029/2006GL027819.

  8. Remias D, Jost S, Boenigk J, Wastian J, Lütz C. Hydrurus-related golden algae (Chrysophyceae) cause yellow snow in polar summer snowfields. Phycol Res. 2013;61:277–85.

    Article  CAS  Google Scholar 

  9. Takeuchi N, Tanaka S, Konno Y, Irvine-Fynn TDL, Rassner SME, Edwards A. Variations in phototroph communities on the ablating bare-ice surface of glaciers on Brøggerhalvøya, Svalbard. Front Earth Sci. 2019;7:1–10.

    Google Scholar 

  10. Procházková L, Leya T, Krížková H, Nedbalová L. 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. 2019;95:1–21.

    Article  Google Scholar 

  11. Lutz S, Anesio AM, Edwards A, Benning LG. Linking microbial diversity and functionality of arctic glacial surface habitats. Environ Microbiol. 2017;19:551–65.

    Article  CAS  PubMed  Google Scholar 

  12. Lutz S, Anesio AM, Edwards A, Benning LG. Microbial diversity on icelandic glaciers and ice caps. Front Microbiol. 2015;6:307.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Gray A, Krolikowski M, Fretwell P, Convey P, Peck LS, Mendelova M, et al. Remote sensing phenology of Antarctic green and red snow algae using WorldView satellites. Front Plant Sci. 2021;12:1–16.

    Article  Google Scholar 

  14. Lutz S, McCutcheon J, McQuaid JB, Benning LG. The diversity of ice algal communities on the Greenland ice sheet as revealed by oligotyping. Microb genomics. 2018;4:1–10.

    Article  CAS  Google Scholar 

  15. Stibal M, Box JE, Cameron KA, Langen PL, Yallop ML, Mottram RH, et al. Algae drive enhanced darkening of bare ice on the Greenland ice sheet. Geophys Res Lett. 2017;44:11,463-71.

    Article  Google Scholar 

  16. Cook JM, Tedstone AJ, Williamson C, McCutcheon J, Hodson AJ, Dayal A, et al. Glacier algae accelerate melt rates on the south-western Greenland ice sheet. Cryosphere. 2020;14:309–30.

    Article  Google Scholar 

  17. Cook JM, Tedstone AJ, Williamson C, McCutcheon J, Hodson AJ, Dayal A, et al. Glacier algae accelerate melt rates on the western Greenland ice sheet. Cryosph Discuss. 2019;1–31. https://doi.org/10.5194/tc-2019-58.

  18. Williamson CJ, Cook J, Tedstone A, Yallop M, McCutcheon J, Poniecka E, et al. Algal photophysiology drives darkening and melt of the Greenland ice sheet. Proc Natl Acad Sci. 2020;201918412. Available from: http://www.pnas.org/lookup/doi/10.1073/pnas.1918412117

  19. Lutz S, Anesio AM, Raiswell R, Edwards A, Newton RJ, Gill F, et al. The biogeography of red snow microbiomes and their role in melting arctic glaciers. Nat Commun. 2016;7:1–9.

    Article  Google Scholar 

  20. Williamson CJ, Cameron KA, Cook JM, Zarsky JD, Stibal M, Edwards A. Glacier algae: a dark past and a darker future. Front Microbiol. 2019;10:524.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhang YZ, Shi M, Holmes EC. Using metagenomics to characterize an expanding virosphere. Cell. 2018;172:1168–72. Available from: https://doi.org/10.1016/j.cell.2018.02.043

  22. Graham EB, Paez-Espino D, Brislawn C, Hofmockel KS, Wu R, Kyrpides NC, et al. Untapped viral diversity in global soil metagenomes. bioRxiv. 2019;583997. https://doi.org/10.1101/583997.

  23. Lara E, Vaqué D, Sà EL, Boras JA, Gomes A, Borrull E, et al. Unveiling the role and life strategies of viruses from the surface to the dark ocean. Sci Adv. 2017;3.

  24. Lara E, Roux S, Sullivan MB, Luna GM, Acinas SG, Vaqué D, et al. An inside look at bathypelagic viruses. 2015;54:6237

  25. Gao C, Xia J, Zhou X, Liang Y, Jiang Y, Wang M, et al. Viral characteristics of the warm Atlantic and cold Arctic water masses in the Nordic seas. Appl Environ Microbiol. 2021;87(22):e0116021.

    Article  PubMed  Google Scholar 

  26. Hingamp P, Grimsley N, Acinas SG, Clerissi C, Subirana L, Poulain J, et al. Exploring nucleo-cytoplasmic large DNA viruses in Tara Oceans microbial metagenomes. ISME J. 2013;7:1678–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bäckström D, Yutin N, Jørgensen SL, Dharamshi J, Homa F, Zaremba-Niedwiedzka K, et al. Virus genomes from deep sea sediments expand the ocean megavirome and support independent origins of viral gigantism. MBio. 2019;10(2):e02497-18.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Anesio AM, Mindl B, Laybourn-Parry J, Hodson AJ, Sattler B. Viral dynamics in cryoconite holes on a high Arctic glacier (Svalbard). J Geophys Res Biogeosciences. 2007;112:G4. https://doi.org/10.1029/2006JG000350.

  29. Bellas CM, Anesio AM, Barker G. Analysis of virus genomes from glacial environments reveals novel virus groups with unusual host interactions. Front Microbiol. 2015;6:1–14.

    Article  Google Scholar 

  30. Bellas CM, Anesio AM, Telling J, Stibal M, Tranter M, Davis S. Viral impacts on bacterial communities in Arctic cryoconite. Environ Res Lett. 2013;8(4):045021. https://doi.org/10.1088/1748-9326/8/4/045021.

    Article  Google Scholar 

  31. Bellas CM, Anesio AM. High diversity and potential origins of T4-type bacteriophages on the surface of Arctic glaciers. Extremophiles. 2013;17:861–70.

    Article  CAS  PubMed  Google Scholar 

  32. Liu Y, Jiao N, Xu Zhong K, Zang L, Zhang R, Xiao X, et al. Diversity and function of mountain and polar supraglacial DNA viruses. Sci Bull. 2023;68(20):2418–33.

    Article  CAS  Google Scholar 

  33. Barno AR, Green K, Rohwer F, Silveira CB. Viral and bacterial ecogenomics in globally expanding red snow blooms. bioRxiv. 2023. Available from: https://www.biorxiv.org/content/https://doi.org/10.1101/2023.01.07.523113v1.full#page

  34. Clouthier SC, Vanwalleghem E, Copeland S, Klassen C, Hobbs G, Nielsen O, et al. A new species of nucleo-cytoplasmic large DNA virus (NCLDV) associated with mortalities in Manitoba lake sturgeon Acipenser fulvescens. Dis Aquat Organ. 2013;102:195–209.

    Article  CAS  PubMed  Google Scholar 

  35. Labbé M, Thaler M, Pitot TM, Rapp JZ, Vincent WF, Culley AI. Climate-endangered Arctic epishelf lake harbors viral assemblages with distinct genetic repertoires. Appl Environ Microbiol. 2022;88:1–15.

    Article  Google Scholar 

  36. Mann NH, Cook A, Millard A, Bailey S, Clokie M. Bacterial photosynthesis genes in a virus. Nature. 2003;424:741.

    Article  CAS  PubMed  Google Scholar 

  37. Williamson SJ, Rusch DB, Yooseph S, Halpern AL, Heidelberg KB, Glass JI, et al. The sorcerer II global ocean sampling expedition: metagenomic characterization of viruses within aquatic microbial samples. PLoS One. 2008;3(1):e1456.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Kuhlisch C, Schleyer G, Shahaf N, Vincent F, Schatz D, Vardi A. Viral infection of algal blooms leaves a unique metabolic footprint on the dissolved organic matter in the ocean. Sci Adv. 2021;7:1–14.

    Article  Google Scholar 

  39. Suttle CA. Marine viruses - major players in the global ecosystem. Nat Rev Microbiol. 2007;5:801–12.

    Article  CAS  PubMed  Google Scholar 

  40. Van Etten JL, Dunigan DD, Nagasaki K, Schroeder DC, Grimsley N, Brussaard CPD, et al. Phycodnaviruses (Phycodnaviridae). Encycl. Virol. Elsevier Ltd.; 2021. Available from: https://doi.org/10.1016/B978-0-12-809633-8.21291-0

  41. Brussaard CPD, Kuipers B, Veldhuis MJW. A mesocosm study of Phaeocystis globosa population dynamics: I. Regulatory role of viruses in bloom control. Harmful Algae. 2005;4:859–74.

    Article  Google Scholar 

  42. Schulz F, Roux S, Paez-Espino D, Jungbluth S, Walsh DA, Denef VJ, et al. Giant virus diversity and host interactions through global metagenomics. Nature. 2020;578:432–6. Available from: https://doi.org/10.1038/s41586-020-1957-x

  43. Zhong ZP, Vik D, Rapp JZ, Zablocki O, Maughan H, Temperton B, et al. Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines. Microbiome. 2023;11:1–18.

    Article  Google Scholar 

  44. Legendre M, Bartoli J, Shmakova L, Jeudy S, Labadie K, Adrait A, et al. Thirty-thousand-year-old distant relative of giant icosahedral DNA viruses with a pandoravirus morphology. Proc Natl Acad Sci U S A. 2014;111:4274–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Philippe N, Legendre M, Doutre G, Couté Y, Poirot O, Lescot M, et al. Pandoraviruses: Amoeba viruses with genomes up to 25 Mb reaching that of parasitic eukaryotes. Science(80- ). 2013;341:281–6.

    CAS  Google Scholar 

  46. Van Etten JL, Lane LC, Dunigan DD. DNA viruses: the really big ones (giruses). Annu Rev Microbiol. 2010;64:83–99.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Schulz F, Abergel C, Woyke T. Giant virus biology and diversity in the era of genome-resolved metagenomics. Nat Rev Microbiol. 2022;20:721–36.

    Article  CAS  PubMed  Google Scholar 

  48. Moniruzzaman M, Weinheimer AR, Martinez-Gutierrez CA, Aylward FO. Widespread endogenization of giant viruses shapes genomes of green algae. Nature. 2020;588:141–5.

    Article  CAS  PubMed  Google Scholar 

  49. Gallot-Lavallée L, Blanc G. A glimpse of nucleo-cytoplasmic large DNA virus biodiversity through the eukaryotic genomicswindow. Viruses. 2017;9(1):17.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Aylward FO, Moniruzzaman M, Ha AD, Koonin E V. A phylogenomic framework for charting the diversity and evolution of giant viruses. PLoS Biol. 2021;19:1–18. Available from: https://doi.org/10.1371/journal.pbio.3001430

  51. Moniruzzaman M, Martinez-Gutierrez CA, Weinheimer AR, Aylward FO. Dynamic genome evolution and complex virocell metabolism of globally-distributed giant viruses. Nat Commun. 2020;11:1–11. Available from: https://doi.org/10.1038/s41467-020-15507-2

  52. Ha AD, Moniruzzaman M, Aylward FO. Assessing the biogeography of marine giant viruses in four oceanic transects. ISME Commun. 2023;3:1–13.

    Article  Google Scholar 

  53. Andrade ACDSP, Arantes TS, Rodrigues RAL, Machado TB, Dornas FP, Landell MF, et al. Ubiquitous giants: a plethora of giant viruses found in Brazil and Antarctica. Virol J. 2018;15:1–10.

    Article  Google Scholar 

  54. López-Bueno A, Tamames J, Velázquez D, Moya A, Quesada A, Alcamí A. High diversity of the viral community from an Antarctic lake. Science(80-). 2009;326:858–61.

    Google Scholar 

  55. Yau S, Lauro FM, DeMaere MZ, Brown MV, Thomas T, Raftery MJ, et al. Virophage control of antarctic algal host-virus dynamics. Proc Natl Acad Sci U S A. 2011;108:6163–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Legendre M, Audic S, Poirot O, Hingamp P, Seltzer V, Byrne D, et al. mRNA deep sequencing reveals 75 new genes and a complex transcriptional landscape in Mimivirus. Genome Res. 2010;20:664–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Aylward FO, Moniruzzaman M. Viralrecall—a flexible command-line tool for the detection of giant virus signatures in ‘omic data. Viruses. 2021;13:15–7.

    Article  Google Scholar 

  58. Karki S, Moniruzzaman M, Aylward FO. Comparative genomics and environmental distribution of large dsDNA viruses in the family Asfarviridae. Front Microbiol. 2021;12:1–13.

    Article  Google Scholar 

  59. Moniruzzaman M, Erazo-Garcia MP, Aylward FO. Endogenous giant viruses contribute to intraspecies genomic variability in the model green alga Chlamydomonas reinhardtii. Virus Evol. 2022;8(2):veac102.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Nagasaki K, Tarutani K, Yamaguchi M. Growth characteristics of Heterosigma akashiwo virus and its possible use as a microbiological agent for red tide control. Appl Environ Microbiol. 1999;65:898.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Aylward FO, Abrahão JS, Brussaard CPD, Fischer MG, Moniruzzaman M, Ogata H, et al. Taxonomic update for giant viruses in the order Imitervirales (phylum Nucleocytoviricota). Arch Virol. 2023;168:283. Available from: http://www.ncbi.nlm.nih.gov/pubmed/37904060

  62. Schvarcz CR, Steward GF. A giant virus infecting green algae encodes key fermentation genes. Virology. 2018;518:423–33. Available from: https://doi.org/10.1016/j.virol.2018.03.010

  63. Endo H, Blanc-Mathieu R, Li Y, Salazar G, Henry N, Labadie K, et al. Biogeography of marine giant viruses reveals their interplay with eukaryotes and ecological functions. Nat Ecol Evol. 2020;4:1639–49. Available from: https://doi.org/10.1038/s41559-020-01288-w

  64. Ha AD, Moniruzzaman M, Aylward FO. High transcriptional activity and diverse functional repertoires of hundreds of giant viruses in a coastal marine system. mSystems. 2021;6(4):e0029321.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  66. Luo W, Ding H, Li H, Ji Z, Huang K, Zhao W, et al. Molecular diversity of the microbial community in coloured snow from the Fildes Peninsula (King George Island, Maritime Antarctica). Polar Biol. 2020;43:1391–405. Available from: https://doi.org/10.1007/s00300-020-02716-0

  67. Hardge K, Peeken I, Neuhaus S, Lange BA, Stock A, Stoeck T, et al. The importance of sea ice for exchange of habitat-specific protist communities in the Central Arctic Ocean. J Mar Syst. 2017;165:124–38.

    Article  Google Scholar 

  68. Howe AT, Bass D, Chao EE, Cavalier-Smith T. New genera, species, and improved phylogeny of Glissomonadida (Cercozoa). Protist. 2011;162:710–22. Available from: https://doi.org/10.1016/j.protis.2011.06.002

  69. Hess S, Melkonian M. The mystery of clade X: Orciraptor gen. nov. and viridiraptor gen. nov. are highly specialised, algivorous amoeboflagellates (Glissomonadida, Cercozoa). Protist. 2013;164:706–47. Available from: https://doi.org/10.1016/j.protis.2013.07.003

  70. Przytulska A, Comte J, Crevecoeur S, Lovejoy C, Laurion I, Vincent WF. Phototrophic pigment diversity and picophytoplankton in permafrost thaw lakes. Biogeosciences. 2016;13:13–26.

    Article  CAS  Google Scholar 

  71. Posch T, Eugster B, Pomati F, Pernthaler J, Pitsch G, Eckert EM. Network of interactions between ciliates and phytoplankton during spring. Front Microbiol. 2015;6:1–14.

    Article  Google Scholar 

  72. Lanzén A, Jørgensen SL, Huson DH, Gorfer M, Grindhaug SH, Jonassen I, et al. CREST - Classification Resources for Environmental Sequence Tags. PLoS One. 2012;7(11):e49334.

    Article  PubMed  PubMed Central  Google Scholar 

  73. McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8(4):e61217.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Team RC. R: a language and environment for statistical computing. 2023.

  75. Bischoff HW, Bold HC. Phycological studies IV. Some soil algae from Enchanted Rock and related algal species. Univ Texas Publ. 1963;6318:1–95.

    Google Scholar 

  76. Jackman SD, Vandervalk BP, Mohamadi H, Chu J, Yeo S, Hammond SA, et al. ABySS 2 . 0: resource-efficient assembly of large genomes using a bloom filter effect of bloom filter false positive rate. Genome Res. 2017;27:768–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Kolmogorov M, Yuan J, Lin Y, Pevzner PA. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol. 2019;37:540–6. Available from: https://doi.org/10.1038/s41587-019-0072-8

  78. Cheng S, Xian W, Fu Y, Marin B, Keller J, Wu T, et al. Genomes of subaerial Zygnematophyceae provide insights into land plant evolution. Cell. 2019;179:1057-1067.e14.

    Article  CAS  PubMed  Google Scholar 

  79. Salmela L, Rivals E. LoRDEC: accurate and efficient long read error correction. Bioinformatics. 2014;30:3506–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Chen S, Zhou Y, Chen Y, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34:i884–90.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. MetaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Grabherr MG, Haas Brian J., Yassour Moran, Levin Joshua Z., Thompson Dawn A., Amit Ido, Adiconis Xian, Fan Lin, Raychowdhury Raktima, Zeng Qiandong, Chen Zehua, Mauceli Evan, Hacohen Nir, Gnirke Andreas, Rhind Nicholas, Palma Federica di, W. N Bruce, Friedman, AR. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2013;29:644–52.

    Article  Google Scholar 

  83. Kang DD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. PeerJ. 2015;2015:1–15.

    Google Scholar 

  84. Paysan-Lafosse T, Blum M, Chuguransky S, Grego T, Pinto BL, Salazar GA, et al. InterPro in 2022. Nucleic Acids Res. 2023;51:D418–27.

    Article  CAS  PubMed  Google Scholar 

  85. John SG, Mendez CB, Deng L, Poulos B, Kauffman AKM, Kern S, et al. A simple and efficient method for concentration of ocean viruses by chemical flocculation. Environ Microbiol Rep. 2011;3:195–202.

    Article  CAS  PubMed  Google Scholar 

  86. Hurwitz BL, Sullivan MB. The Pacific Ocean virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. Thompson F, editor. PLoS One. 2013;8:e57355.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Katoh K, Kuma KI, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37:1530–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Kalyaanamoorthy S, Minh BQ, Wong TKF, Von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Hoang DT, Chernomor O, Von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35:518–22.

    Article  CAS  PubMed  Google Scholar 

  91. Gilbert NE, LeCleir GR, Pound HL, Strzepek RF, Ellwood MJ, Twining BS, et al. Giant virus infection signatures are modulated by euphotic zone depth strata and iron regimes of the subantarctic southern ocean. mSystems. 2023;8(2):e0126022.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Special thank you for help with sample collection and processing go to Daniel Rissi and Christopher Trivedi (GFZ) as well as Marie Bolander Jensen (Århus University). We thank the scientific field teams for their support in 2019 and 2020 (James A. Bradley, Matthias Winkel, Eva L. Doting, Laura Halbach, and Rey Mourot). We particularly thank Thomas Leya (Fraunhofer Institute for Cell Therapy and Immunology, Germany) who manager the CCCRyo Culture Collection (http://cccryo.fraunhofer.de/web/infos/welcome/) for culturing and providing the five snow algal cultures from which the genomic assemblies were sequenced.

Funding

This study was financially supported by the European Research Council (ERC) Synergy Grant DEEP PURPLE under the European Union’s Horizon 2020 Research and Innovation Program (Grant Number 856416) awarded to M. T, L. G. B., and A. A. L. G. B. acknowledges financial support from an H2020 EU-funded INTERACT project (AirMiMic, grant agreement No. 730938) through which the fieldwork in SE Greenland was possible.

Author information

Authors and Affiliations

Authors

Contributions

LP and AMA conceived and designed the study. LGB and RM performed sampling and processing of the samples in the field. LP processed the environmental samples in the lab, and AZ generated sequencing data. SL generated the snow algal assemblies. LP, KS, and AMA performed the data analysis and interpretation of the results, with support from CB, MM, and AZ. LP and KS drafted the manuscript. AZ, CB, SL, MM, LGB, MT, and AMA reviewed and edited the manuscript. AMA, LGB, and MT acquired funding for the study.

Corresponding author

Correspondence to Laura Perini.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

40168_2024_1796_MOESM1_ESM.zip

Additional file 1. Figure S1-S7: Maximum-likelihood phylogenetic tree of the NCLDV core gene D5, RNAps, RNApl, MCP, mRNAc, A32, SFII, VLTF3, RNR. Sequences recovered from the environmental samples are presented in bold. Environmental sample types are specified in correspondence of each sequence. Branches are color-coded by order-level taxonomy. Dark dots at the nodes represent the bootstrap support value of >70. Figure S8: Normalized counts of the NCLDV marker genes by each assembly size for all the samples in this study. Analysis was carried out in 19 environmental metagenomes (MG) and 18 environmental metatranscriptomes (pooled) obtained from samples of cryoconite (n=1), ice core (n=3), green snow (n=2), red snow (n=5), dark ice (n=8), one metavirome (dark ice) and five snow algae culture genomic assemblies from the CCCryo collection. The ‘md’ (more depth) notation following select samples are those that were re-sequenced with higher metagenomic coverage. Symbols represent the sample types. Figure S9: Metagenomic SSU relative abundance. Sample descriptions can be found in Table 1 and Supplementary Table 1. The blank space above each bar is comprised of bacterial phyla. The ‘md’ (more depth) following a sample name marks those that were sequenced with a high average library coverage. The full abundance table can be found in Supplementary Table S10. Figure S10: Read mapping percentage transformed into log scale + 0.01 for comparable scales. Purple values are 0 reads mapped. Assemblies (left) and sample read files (right), MG are metagenomes, MG with ‘_2’ note the four samples that were sequenced with higher library coverage, and ‘MT’ are metatranscriptomes. The only culture to have reads recruit was Raphidonema_sempervirens_LIV13260. Sample types are labeled in the same way as Fig. 2.

40168_2024_1796_MOESM2_ESM.xlsx

Additional file 2. Table S1: Samples info. Table S2: ViralRecall GV Count. Table S3: Normalized VRGV Counts. Table S4: NCBI QC GV Counts. Table S5: % identity PolB. Table S6: GVMAG. Table S7: GVMAG ANI %. Table S8: GV MAGs funct ann. Table S9: Percentages 18S rRNA MT. Table S10: Abun Table MetagenomeSSU. Table S11: BARRNAP results. Table S12: CoverM-Mapping

Additional file 3.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Perini, L., Sipes, K., Zervas, A. et al. Giant viral signatures on the Greenland ice sheet. Microbiome 12, 91 (2024). https://doi.org/10.1186/s40168-024-01796-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40168-024-01796-y