Distinct and rich assemblages of giant viruses in Arctic and Antarctic lakes

Giant viruses (GVs) are key players in ecosystem functioning, biogeochemistry, and eukaryotic genome evolution. GV diversity and abundance in aquatic systems can exceed that of prokaryotes, but their diversity and ecology in lakes, especially polar ones, remain poorly understood. We conducted a comprehensive survey and meta-analysis of GV diversity across 20 lakes, spanning polar to temperate regions, combining our extensive lake metagenome database from the Canadian Arctic and subarctic with publicly available datasets. Leveraging a novel GV genome identification tool, we identified 3304 GV metagenome-assembled genomes, revealing lakes as untapped GV reservoirs. Phylogenomic analysis highlighted their dispersion across all Nucleocytoviricota orders. Strong GV population endemism emerged between lakes from similar regions and biomes (Antarctic and Arctic), but a polar/temperate barrier in lacustrine GV populations and differences in their gene content could be observed. Our study establishes a robust genomic reference for future investigations into lacustrine GV ecology in fast changing polar environments.

A Tale of Two Seasons: Distinct Seasonal Viral Communities in a Thermokarst Lake

Thermokarst lakes are important features of subarctic landscapes and are a substantial source of greenhouse gases, although the extent of gas produced varies seasonally. Microbial communities are responsible for the production of methane and CO₂ but the "top down" forces that influence microbial dynamics (i.e., grazers and viruses) and how they vary temporally within these lakes are still poorly understood. The aim of this study was to examine viral diversity over time to elucidate the seasonal structure of the viral communities in thermokarst lakes. We produced virus-enriched metagenomes from a subarctic peatland thermokarst lake in the summer and winter over three years. The vast majority of vOTUs assigned to viral families belonged to Caudovirales (Caudoviricetes), notably the morphological groups myovirus, siphovirus and podovirus. We identified two distinct communities: a dynamic, seasonal community in the oxygenated surface layer during the summer and a stable community found in the anoxic water layer at the bottom of the lake in summer and throughout much of the water column in winter. Comparison with other permafrost and northern lake metagenomes highlighted the distinct composition of viral communities in this permafrost thaw lake ecosystem.

Structural shifts in sea ice prokaryotic communities across a salinity gradient in the subarctic

Current knowledge of the processes that shape prokaryotic community assembly in sea ice across polar ecosystems is scarce. Here, we coupled culture-dependent (bacterial isolation on R2A medium) and culture-independent (high-throughput 16S rRNA gene sequencing) approaches to provide the first comprehensive assessment of prokaryotic communities in the late winter ice and its underlying water along a natural salinity gradient in coastal Hudson Bay, an iconic cryo-environment that marks the ecological transition between Canadian Subarctic and Arctic biomes. We found that prokaryotic community assembly processes in the ice were less selective at low salinity since typical freshwater taxa such as Frankiales, Burkholderiales, and Chitinophagales dominated both the ice and its underlying water. In contrast, there were sharp shifts in community structure between the ice and underlying water samples at sites with higher salinity, with the orders Alteromonadales and Flavobacteriales dominating the ice, while the abovementioned freshwater taxa dominated the underlying water communities. Moreover, primary producers including Cyanobium (Cyanobacteria, Synechococcales) may play a role in shaping the ice communities and were accompanied by known Planctomycetes and Verrucomicrobiae taxa. Culture-dependent analyses showed that the ice contained pigment-producing psychrotolerant or psychrophilic bacteria from the phyla Proteobacteria, Actinobacteriota, and Bacteroidota, likely favored by the combination of low temperatures and the seasonal increase in sunlight. Our findings suggest that salinity, photosynthesis and dissolved organic matter are the main drivers of prokaryotic community structure in the late winter ice of coastal Hudson Bay, the ecosystem with the fastest sea ice loss rate in the Canadian North.

Size-Fractionated Microbiome Structure in Subarctic Rivers and a Coastal Plume Across DOC and Salinity Gradients

Little is known about the microbial diversity of rivers that flow across the changing subarctic landscape. Using amplicon sequencing (rRNA and rRNA genes) combined with HPLC pigment analysis and physicochemical measurements, we investigated the diversity of two size fractions of planktonic Bacteria, Archaea and microbial eukaryotes along environmental gradients in the Great Whale River (GWR), Canada. This large subarctic river drains an extensive watershed that includes areas of thawing permafrost, and discharges into southeastern Hudson Bay as an extensive plume that gradually mixes with the coastal marine waters. The microbial communities differed by size-fraction (separated with a 3-μm filter), and clustered into three distinct environmental groups: (1) the GWR sites throughout a 150-km sampling transect; (2) the GWR plume in Hudson Bay; and (3) small rivers that flow through degraded permafrost landscapes. There was a downstream increase in taxonomic richness along the GWR, suggesting that sub-catchment inputs influence microbial community structure in the absence of sharp environmental gradients. Microbial community structure shifted across the salinity gradient within the plume, with changes in taxonomic composition and diversity. Rivers flowing through degraded permafrost had distinct physicochemical and microbiome characteristics, with allochthonous dissolved organic carbon explaining part of the variation in community structure. Finally, our analyses of the core microbiome indicated that while a substantial part of all communities consisted of generalists, most taxa had a more limited environmental range and may therefore be sensitive to ongoing change.

Chemical Composition of the Unexplored Volatile Fraction of Betula glandulosa, a Prevalent Shrub in Nunavik, Quebec

The volatile fraction of the leaves of Betula glandulosa Michx. has been investigated for its secondary metabolite composition by GC-MS and GC/FID. The rapid expansion of this shrub species in subarctic landscapes, like the ones found in Nunavik (Northern Québec, Canada), highly impacts ecosystem dynamics. Yet, despite its abundance, few phytochemical investigations have yet been conducted on this species. In this study, we present the first phytochemical investigation of the volatile metabolites of B. glandulosa leaves. Although no essential oil was isolated, volatile compounds were extracted from the hydrosol by steam distillation. The main metabolites observed were linalool (14.6 - 19.0%), C 6 oxylipins (known as green leaf volatiles, GLV; total of 18.2 - 40.2%), eugenol (1.6 - 8.6%) and α -terpineol (3.3 - 4.8%). Dwarf birch is an important food source for insects and herbivores, so knowledge of its metabolite composition could help understand parts of its functional role in subarctic ecosystems. The composition of the volatile fraction could serve as marker for differentiating B. glandulosa from other dwarf birch species like Betula nana L.

Functionally explicit partitioning of plant β-diversity reveal soil fungal assembly in the subarctic tundra

Metabarcoding technologies for soil fungal DNA pools have enabled to capture the diversity of fungal community and the agreement of their β-diversity with plant β-diversity. However, processes underlying the synchrony of the aboveground-belowground biodiversity is still unclear. By using partitioning methods for plant β-diversity, this study explored the process driving synchrony in tundra ecosystems, in which drastic vegetation shifts are observed with climate warming. Our methods based on Baselga's partitioning enabled the division of plant β-diversity into two phenomena and three functional components. Correlation of fungal β-diversity with the components of plant β-diversity showed that the spatial replacement of fungi was promoted by plant species turnover, in particular, plant species turnover with functional exchange. In addition, spatial variety of graminoid or forbs species, rather than shrubs, enhanced fungal β-diversity. These results suggest the importance of small-scale factors such as plant-fungal interactions or local environments modified by plants for the fungal community assemblage. The process-based understanding of community dynamics of plants and fungi allows us to predict the ongoing shrub encroachment in the Arctic region, which could weaken the aboveground-belowground synchrony.

Community composition of aquatic fungi across the thawing Arctic

Thermokarst activity at permafrost sites releases considerable amounts of ancient carbon to the atmosphere. A large part of this carbon is released via thermokarst ponds, and fungi could be an important organismal group enabling its recycling. However, our knowledge about aquatic fungi in thermokarstic systems is extremely limited. In this study, we collected samples from five permafrost sites distributed across circumpolar Arctic and representing different stages of permafrost integrity. Surface water samples were taken from the ponds and, additionally, for most of the ponds also the detritus and sediment samples were taken. All the samples were extracted for total DNA, which was then amplified for the fungal ITS2 region of the ribosomal genes. These amplicons were sequenced using PacBio technology. Water samples were also collected to analyze the chemical conditions in the ponds, including nutrient status and the quality and quantity of dissolved organic carbon. This dataset gives a unique overview of the impact of the thawing permafrost on fungal communities and their potential role on carbon recycling.

Re-activation of landslide in sub-Arctic areas due to extreme rainfall and discharge events (the mouth of the Great Whale River, Nunavik, Canada)

Climate change is impacting surficial geomorphic processes, especially in sensitive areas such as the sub-Arctic. One of the most common examples involves landslides, which often develop in glacio-isostatically raised marine clays in northeastern Canada. One of these sites, an expansive area of complex landslide terrain located at the mouth of the Great Whale River in Nunavik, has already been studied due to its age and morphology. We present new data, based on the multidisciplinary research including geomorphic, dendrochronological, and hydroclimatological analyses, allowing us to determine how contemporary climate change has affected landslide reactivation during the last 80 years. Our research included collecting 60 cores from Picea glauca trees, growing on the marginal zone of a landslide deposit, as well as from a reference site. The tilted trees formed eccentric growth-ring patterns, which provided us with reliable dates on the landslide events. In addition to these dendrochronological data, we studied these landslides using repeated aerial photography, which showed changes in river channel constrictions in the period 1969-2019. Based on the eccentricity index of the tree ring data, we recognized disturbance events due to landslides. We compared these data with the hydroclimatological conditions and found clearly visible correlations between heavy rainfall and discharge (>95th percentile) of the Great Whale River. The increased landslide activity over the past several years can be linked to an increase in extreme summertime rainfall events. Increased landslide activity poses a real threat, through its input of large amounts of fine-grained sediment to the river, causing it to narrow.

Structural Characterization of Dissolved Organic Matter in Permafrost Peatland Lakes

Thermokarst lakes result from the thawing of ice-rich permafrost and are widespread across northern landscapes. These waters are strong emitters of methane, especially in permafrost peatland regions, where they are stained black by high concentrations of dissolved organic matter (DOM). In the present study, we aimed to structurally characterize the DOM from a set of peatland thermokarst lakes that are known to be intense sites of microbial decomposition and methane emission. Samples were collected at different depths from three thermokarst lakes in the Sasapimakwananisikw (SAS) River valley near the eastern Hudson Bay community of Kuujjuarapik–Whapmagoostui (Nunavik, Canada). Samples were analyzed by spectrofluorometry, Fourier-transform infrared spectroscopy (FTIR), nuclear magnetic resonance spectroscopy (NMR), and elemental analysis. Fluorescence analyses indicated considerable amounts of autochthonous DOM in the surface waters of one of SAS 1A, indicating a strong bioavailability of labile DOM, and consequently a greater methanogenic potential. The three lakes differed in their chemical composition and diversity, suggesting various DOM transformations phenomena. The usefulness of complementary analytical approaches to characterize the complex mixture of DOM in permafrost peatland waters cannot be overlooked, representing a first step towards greater comprehension of the organic geochemical properties of these permafrost-derived systems.

Hotspots of biotic compositional change in lakes along vast latitudinal transects in northern Canada

Ecotones mark zones of rapid change in ecological structure at various spatial scales. They are believed to be particularly susceptible to shifts caused by environmental transformation, making them key regions for studying the effects of global change. Here, we explored the variation in assemblage structure of aquatic primary producer and consumer communities across latitudinal transects in northeastern North America (Québec-Labrador) to identify spatial patterns in biodiversity that indicated the location of transition zones across the landscape. We analyzed species richness and the cumulative rate of compositional change (expressed as beta-diversity) of diatoms and chironomids to detect any abrupt shifts in the rate of spatial taxonomic turnover. We used principal coordinates analysis to estimate community turnover with latitude, then applied piecewise linear regression to assess the position of ecotones. Statistically significant changes in assemblage composition occurred at 52 and 55°N, corresponding to the transition between closed- and open-crown forest, and to the southern onset of the forest tundra (i.e., the forest limit), respectively. The spatial distribution of ecotones was most strongly related to air temperature for chironomids and to vegetation- and soil-related chemical attributes of lake water for diatoms, including dissolved organic carbon content and water color. Lakes at mid- to high-latitudes currently face pressures from rapidly rising temperatures, accompanied by large increases in organic carbon inputs from their catchments, often leading to browning and its associated effects. The biota at the base of food webs in lakes located in transition zones are disproportionately affected by the cascading effects of these multi-factorial changes, concurrent with pronounced terrestrial greening observed in these regions. Similar patterns of biotic shifts have been observed along alpine aquatic transects, indicating the potential for widespread restructuring of cold, high-altitude and high-latitude freshwater communities due to global change.

Identification of a Threshold Minimum Area for Reflectance Retrieval from Thermokarst Lakes and Ponds Using Full-Pixel Data from Sentinel-2

Thermokarst waterbodies caused by permafrost thawing and degradation are ubiquitous in many subarctic and Arctic regions. They are globally important components of the biogeochemical carbon cycle and have potential feedback effects on climate. These northern waters are mostly small lakes and ponds, and although they may be mapped using very high-resolution satellites or aerial photography, these approaches are generally not suitable for monitoring purposes, due to the cost and limited availability of such images. In this study we evaluated the potential use of widely available high-resolution imagery from Sentinel-2 (S2) for the characterization of the spectral reflectance of thermokarst lakes and ponds. Specifically, we aimed to define the minimum lake area that could be reliably imaged, and to identify challenges and solutions for remote sensing of such waters in the future. The study was conducted in subarctic Canada, in the vicinity of Whapmagoostui-Kuujjuarapik (Nunavik, Québec), an area in the sporadic permafrost zone with numerous thermokarst waterbodies that vary greatly in size. Ground truthing lake reflectance data were collected using an Unmanned Aerial System (UAS) fitted with a multispectral camera that collected images at 13 cm resolution. The results were compared with reflectance from Sentinel-2 images, and the effect of lake area on the reflectance response was assessed. Our results show that Sentinel-2 imagery was suitable for waterbodies larger than 350 m2 once their boundaries were defined, which in the two test sites would allow monitoring from 11% to 30% of the waterbodies and 73% to 85% of the total lake area. Challenges for remote sensing of small lakes include the confounding effects of water reflection (both direct radiation and diffuse), wind and shadow. Given the small threshold area and frequent revisit time, Sentinel-2 provides a valuable approach towards the continuous monitoring of waterbodies, including ponds and small lakes such as those found in thermokarst landscapes. UASs provide a complementary approach for ground truthing and boundary definition.

Diversity and potential activity of methanotrophs in high methane-emitting permafrost thaw ponds

Lakes and ponds derived from thawing permafrost are strong emitters of carbon dioxide and methane to the atmosphere, but little is known about the methane oxidation processes in these waters. Here we investigated the distribution and potential activity of aerobic methanotrophic bacteria in thaw ponds in two types of eroding permafrost landscapes in subarctic Québec: peatlands and mineral soils. We hypothesized that methanotrophic community composition and potential activity differ regionally as a function of the landscape type and permafrost degradation stage, and locally as a function of depth-dependent oxygen conditions. Our analysis of pmoA transcripts by Illumina amplicon sequencing and quantitative PCR showed that the communities were composed of diverse and potentially active lineages. Type I methanotrophs, particularly Methylobacter, dominated all communities, however there was a clear taxonomic separation between the two landscape types, consistent with environmental control of community structure. In contrast, methanotrophic potential activity, measured by pmoA transcript concentrations, did not vary with landscape type, but correlated with conductivity, phosphorus and total suspended solids. Methanotrophic potential activity was also detected in low-oxygen bottom waters, where it was inversely correlated with methane concentrations, suggesting methane depletion by methanotrophs. Methanotrophs were present and potentially active throughout the water column regardless of oxygen concentration, and may therefore be resilient to future mixing and oxygenation regimes in the warming subarctic.

Rare earth elements in freshwater, marine, and terrestrial ecosystems in the eastern Canadian Arctic

Few ecotoxicological studies exist for rare earth elements (REEs), particularly field-based studies on their bioaccumulation and food web dynamics. REE mining has led to significant environmental impacts in several countries (China, Brazil, U.S.), yet little is known about the fate and transport of these contaminants of emerging concern. Northern ecosystems are potentially vulnerable to REE enrichment from prospective mining projects at high latitudes. To understand how REEs behave in remote northern food webs, we measured REE concentrations and carbon and nitrogen stable isotope ratios (∂¹⁵N, ∂¹³C) in biota from marine, freshwater, and terrestrial ecosystems of the eastern Canadian Arctic (N = 339). Wildlife harvesting and tissue sampling was partly conducted by local hunters through a community-based monitoring project. Results show that REEs generally follow a coherent bioaccumulation pattern for sample tissues, with some anomalies for redox-sensitive elements (Ce, Eu). Highest REE concentrations were found at low trophic levels, especially in vegetation and aquatic invertebrates. Terrestrial herbivores, ringed seal, and fish had low total REE levels in muscle tissue (∑REE for 15 elements <0.1 nmol g^-1), yet accumulation was an order of magnitude higher in liver tissues. Age- and length-dependent REE accumulation also suggest that REE uptake is faster than elimination for some species. Overall, REE bioaccumulation patterns appear to be species- and tissue-specific, with limited potential for biomagnification. This study provides novel data on the behaviour of REEs in ecosystems and will be useful for environmental impact assessment of REE enrichment in northern regions.

Different parts, different stories: climate sensitivity of growth is stronger in root collars vs. stems in tundra shrubs

Shrub densification has been widely reported across the circumpolar arctic and subarctic biomes in recent years. Long-term analyses based on dendrochronological techniques applied to shrubs have linked this phenomenon to climate change. However, the multi-stemmed structure of shrubs makes them difficult to sample and therefore leads to non-uniform sampling protocols among shrub ecologists, who will favor either root collars or stems to conduct dendrochronological analyses. Through a comparative study of the use of root collars and stems of Betula glandulosa, a common North American shrub species, we evaluated the relative sensitivity of each plant part to climate variables and assessed whether this sensitivity is consistent across three different types of environments in northwestern Québec, Canada (terrace, hilltop and snowbed). We found that root collars had greater sensitivity to climate than stems and that these differences were maintained across the three types of environments. Growth at the root collar was best explained by spring precipitation and summer temperature, whereas stem growth showed weak and inconsistent responses to climate variables. Moreover, sensitivity to climate was not consistent among plant parts, as individuals having climate-sensitive root collars did not tend to have climate-sensitive stems. These differences in sensitivity of shrub parts to climate highlight the complexity of resource allocation in multi-stemmed plants. Whereas stem initiation and growth are driven by microenvironmental variables such as light availability and competition, root collars integrate the growth of all plant parts instead, rendering them less affected by mechanisms such as competition and more responsive to signals of global change. Although further investigations are required to determine the degree to which these findings are generalizable across the tundra biome, our results indicate that consistency and caution in the choice of plant parts are a key consideration for the success of future dendroclimatological studies on shrubs.

Environmental selection of planktonic methanogens in permafrost thaw ponds

The warming and thermal erosion of ice-containing permafrost results in thaw ponds that are strong emitters of methane to the atmosphere. Here we examined methanogens and other Archaea, in two types of thaw ponds that are formed by the collapse of either permafrost peat mounds (palsas) or mineral soil mounds (lithalsas) in subarctic Quebec, Canada. Using high-throughput sequencing of a hypervariable region of 16S rRNA, we determined the taxonomic structure and diversity of archaeal communities in near-bottom water samples, and analyzed the mcrA gene transcripts from two sites. The ponds at all sites were well stratified, with hypoxic or anoxic bottom waters. Their archaeal communities were dominated by Euryarchaeota, specifically taxa in the methanogenic orders Methanomicrobiales and Methanosarcinales, indicating a potentially active community of planktonic methanogens. The order Methanomicrobiales accounted for most of the mcrA transcripts in the two ponds. The Archaeal communities differed significantly between the lithalsa and palsa ponds, with higher alpha diversity in the organic-rich palsa ponds, and pronounced differences in community structure. These results indicate the widespread occurrence of planktonic, methane-producing Archaea in thaw ponds, with environmental selection of taxa according to permafrost landscape type.

High Methylmercury in Arctic and Subarctic Ponds is Related to Nutrient Levels in the Warming Eastern Canadian Arctic

Permafrost thaw ponds are ubiquitous in the eastern Canadian Arctic, yet little information exists on their potential as sources of methylmercury (MeHg) to freshwaters. They are microbially active and conducive to methylation of inorganic mercury, and are also affected by Arctic warming. This multiyear study investigated thaw ponds in a discontinuous permafrost region in the Subarctic taiga (Kuujjuarapik-Whapmagoostui, QC) and a continuous permafrost region in the Arctic tundra (Bylot Island, NU). MeHg concentrations in thaw ponds were well above levels measured in most freshwater ecosystems in the Canadian Arctic (>0.1 ng L(-1)). On Bylot, ice-wedge trough ponds showed significantly higher MeHg (0.3-2.2 ng L(-1)) than polygonal ponds (0.1-0.3 ng L(-1)) or lakes (<0.1 ng L(-1)). High MeHg was measured in the bottom waters of Subarctic thaw ponds near Kuujjuarapik (0.1-3.1 ng L(-1)). High water MeHg concentrations in thaw ponds were strongly correlated with variables associated with high inputs of organic matter (DOC, a320, Fe), nutrients (TP, TN), and microbial activity (dissolved CO2 and CH4). Thawing permafrost due to Arctic warming will continue to release nutrients and organic carbon into these systems and increase ponding in some regions, likely stimulating higher water concentrations of MeHg. Greater hydrological connectivity from permafrost thawing may potentially increase transport of MeHg from thaw ponds to neighboring aquatic ecosystems.

Bacterial community structure across environmental gradients in permafrost thaw ponds: methanotroph-rich ecosystems

Permafrost thawing leads to the formation of thermokarst ponds that potentially emit CO₂ and CH₄ to the atmosphere. In the Nunavik subarctic region (northern Québec, Canada), these numerous, shallow ponds become well-stratified during summer. This creates a physico-chemical gradient of temperature and oxygen, with an upper oxic layer and a bottom low oxygen or anoxic layer. Our objective was to determine the influence of stratification and related limnological and landscape properties on the community structure of potentially active bacteria in these waters. Samples for RNA analysis were taken from ponds in three contrasting valleys across a gradient of permafrost degradation. A total of 1296 operational taxonomic units were identified by high throughput amplicon sequencing, targeting bacterial 16S rRNA that was reverse transcribed to cDNA. β-proteobacteria were the dominant group in all ponds, with highest representation by the genera Variovorax and Polynucleobacter. Methanotrophs were also among the most abundant sequences at most sites. They accounted for up to 27% of the total sequences (median of 4.9% for all samples), indicating the importance of methane as a bacterial energy source in these waters. Both oxygenic (cyanobacteria) and anoxygenic (Chlorobi) phototrophs were also well-represented, the latter in the low oxygen bottom waters. Ordination analyses showed that the communities clustered according to valley and depth, with significant effects attributed to dissolved oxygen, pH, dissolved organic carbon, and total suspended solids. These results indicate that the bacterial assemblages of permafrost thaw ponds are filtered by environmental gradients, and are complex consortia of functionally diverse taxa that likely affect the composition as well as magnitude of greenhouse gas emissions from these abundant waters.