Brute force prey metabarcoding to explore the diets of small invertebrates.

Prey metabarcoding has become a popular tool in molecular ecology for resolving trophic interactions at high resolution, from various sample types and animals. To date, most predator-prey studies of small-sized animals (<1 mm) have met the problem of overabundant predator DNA in dietary samples by adding blocking primers/peptide nucleic acids. These primers aim to limit the PCR amplification and detection of the predator DNA but may introduce bias to the prey composition identified by interacting with sequences that are similar to those of the predator. Here we demonstrate the use of an alternative method to explore the prey of small marine copepods using whole-body DNA extracts and deep, brute force metabarcoding of an 18S rDNA fragment. After processing and curating raw data from two sequencing runs of varying depths (0.4 and 5.4 billion raw reads), we isolated 1.3 and 52.2 million prey reads, with average depths of ~15,900 and ~120,000 prey reads per copepod individual, respectively. While data from both sequencing runs were sufficient to distinguish dietary compositions from disparate seasons, locations, and copepod species, greater sequencing depth led to better separation of clusters. As computation and sequencing are becoming ever more powerful and affordable, we expect the brute force approach to become a general standard for prey metabarcoding, as it offers a simple and affordable solution to consumers that is impractical to dissect or unknown to science.

Seasonality in land-ocean connectivity and local processes control sediment bacterial community structure and function in a High Arctic tidal flat.

Climate change is altering patterns of precipitation, cryosphere thaw, and land-ocean influxes, affecting understudied Arctic estuarine tidal flats. These transitional zones between terrestrial and marine systems are hotspots for biogeochemical cycling, often driven by microbial processes. We investigated surface sediment bacterial community composition and function from May to September along a river-intertidal-subtidal-fjord gradient. We paired metabarcoding of in situ communities with in vitro carbon-source utilization assays. Bacterial communities differed in space and time, alongside varying environmental conditions driven by local seasonal processes and riverine inputs, with salinity emerging as the dominant structuring factor. Terrestrial and riverine taxa were found throughout the system, likely transported with runoff. In vitro assays revealed sediment bacteria utilized a broader range of organic matter substrates when incubated in fresh and brackish water compared to marine water. These results highlight the importance of salinity for ecosystem processes in these dynamic tidal flats, with the highest potential for utilization of terrestrially derived organic matter likely limited to tidal flat areas (and times) where sediments are permeated by freshwater. Our results demonstrate that intertidal flats must be included in future studies on impacts of increased riverine discharge and transport of terrestrial organic matter on coastal carbon cycling in a warming Arctic.

Linking extreme seasonality and gene expression in Arctic marine protists.

At high latitudes, strong seasonal differences in light availability affect marine organisms and regulate the timing of ecosystem processes. Marine protists are key players in Arctic aquatic ecosystems, yet little is known about their ecological roles over yearly cycles. This is especially true for the dark polar night period, which up until recently was assumed to be devoid of biological activity. A 12 million transcripts catalogue was built from 0.45 to 10 µm protist assemblages sampled over 13 months in a time series station in an Arctic fjord in Svalbard. Community gene expression was correlated with seasonality, with light as the main driving factor. Transcript diversity and evenness were higher during polar night compared to polar day. Light-dependent functions had higher relative expression during polar day, except phototransduction. 64% of the most expressed genes could not be functionally annotated, yet up to 78% were identified in Arctic samples from Tara Oceans, suggesting that Arctic marine assemblages are distinct from those from other oceans. Our study increases understanding of the links between extreme seasonality and biological processes in pico- and nanoplanktonic protists. Our results set the ground for future monitoring studies investigating the seasonal impact of climate change on the communities of microbial eukaryotes in the High Arctic.

Seasonality of the bacterial and archaeal community composition of the Northern Barents Sea.

The Barents Sea is a transition zone between the Atlantic and the Arctic Ocean. The ecosystem in this region is highly variable, and a seasonal baseline of biological factors is needed to monitor the effects of global warming. In this study, we report the results from the investigations of the bacterial and archaeal community in late winter, spring, summer, and early winter along a transect through the northern Barents Sea into the Arctic Ocean east of Svalbard using 16S rRNA metabarcoding. Winter samples were dominated by members of the SAR11 clade and a community of nitrifiers, namely Cand. Nitrosopumilus and LS-NOB (Nitrospinia), suggest a prevalence of chemoautotrophic metabolisms. During spring and summer, members of the Gammaproteobacteria (mainly members of the SAR92 and OM60(NOR5) clades, Nitrincolaceae) and Bacteroidia (mainly Polaribacter, Formosa, and members of the NS9 marine group), which followed a succession based on their utilization of different phytoplankton-derived carbon sources, prevailed. Our results indicate that Arctic marine bacterial and archaeal communities switch from carbon cycling in spring and summer to nitrogen cycling in winter and provide a seasonal baseline to study the changes in these processes in response to the effects of climate change.

The mystery of the ice cold rose-Microbiome of an Arctic winter frost flower.

Under very cold conditions, delicate ice-crystal structures called frost flowers emerge on the surface of newly formed sea ice. These understudied, ephemeral structures include saline brine, organic material, inorganic nutrients, and bacterial and archaeal communities in their brine channels. Hitherto, only a few frost flowers have been studied during spring and these have been reported to be dominated by Rhizobia or members of the SAR11 clade. Here we report on the microbiome of frost flowers sampled during the winter and polar night in the Barents Sea. There was a distinct difference in community profile between the extracted DNA and RNA, but both were dominated by members of the SAR11 clade (78% relative abundance and 41.5% relative activity). The data further suggested the abundance and activity of Cand. Nitrosopumilus, Nitrospinia, and Nitrosomonas. Combined with the inference of marker genes based on the 16S rRNA gene data, this indicates that sulfur and nitrogen cycling are likely the major metabolism in these ephemeral structures.

Multiomics in the central Arctic Ocean for benchmarking biodiversity change.

Multiomics approaches need to be applied in the central Arctic Ocean to benchmark biodiversity change and to identify novel species and their genes. As part of MOSAiC, EcoOmics will therefore be essential for conservation and sustainable bioprospecting in one of the least explored ecosystems on Earth.

High Concentrations of Perfluoroalkyl Acids in Arctic Seawater Driven by Early Thawing Sea Ice.

Poly- and perfluoroalkyl substances are synthetic chemicals that are widely present in the global environment including the Arctic. However, little is known about how these chemicals (particularly perfluoroalkyl acids, PFAA) enter the Arctic marine system and cycle between seawater and sea ice compartments. To evaluate this, we analyzed sea ice, snow, melt ponds, and near-surface seawater at two ice-covered stations located north of the Barents Sea (81 °N) with the aim of investigating PFAA dynamics in the late-season ice pack. Sea ice showed high concentrations of PFAA particularly at the surface with snow-ice (the uppermost sea ice layer strongly influenced by snow) comprising 26-62% of the total PFAA burden. Low salinities (<2.5 ppt) and low δ18OH20 values (<1‰ in snow and upper ice layers) in sea ice revealed the strong influence of meteoric water on sea ice, thus indicating a significant atmospheric source of PFAA with subsequent transfer down the sea ice column in meltwater. Importantly, the under-ice seawater (0.5 m depth) displayed some of the highest concentrations notably for the long-chain PFAA (e.g., PFOA 928 ± 617 pg L-1), which were ≈3-fold higher than those of deeper water (5 m depth) and ≈2-fold higher than those recently measured in surface waters of the North Sea infuenced by industrial inputs of PFAAs. The evidence provided here suggests that meltwater arising early in the melt season from snow and other surface ice floe components drives the higher PFAA concentrations observed in under-ice seawater, which could in turn influence the timing and extent of PFAA exposure for organisms at the base of the marine food web.

Proton-pumping rhodopsins are abundantly expressed by microbial eukaryotes in a high-Arctic fjord.

Proton-pumping rhodopsins provide an alternative pathway to photosynthesis by which solar energy can enter the marine food web. Rhodopsin genes are widely found in marine bacteria, also in the Arctic, and were recently reported from several eukaryotic lineages. So far, little is known about rhodopsin expression in Arctic eukaryotes. In this study, we used metatranscriptomics and 18S rDNA tag sequencing to examine the mid-summer function and composition of marine protists (size 0.45-10 µm) in the high-Arctic Billefjorden (Spitsbergen), especially focussing on the expression of microbial proton-pumping rhodopsins. Rhodopsin transcripts were highly abundant, at a level similar to that of genes involved in photosynthesis. Phylogenetic analyses placed the environmental rhodopsins within disparate eukaryotic lineages, including dinoflagellates, stramenopiles, haptophytes and cryptophytes. Sequence comparison indicated the presence of several functional types, including xanthorhodopsins and a eukaryotic clade of proteorhodopsin. Transcripts belonging to the proteorhodopsin clade were also abundant in published metatranscriptomes from other oceanic regions, suggesting a global distribution. The diversity and abundance of rhodopsins show that these light-driven proton pumps play an important role in Arctic microbial eukaryotes. Understanding this role is imperative to predicting the future of the Arctic marine ecosystem faced by a changing light climate due to diminishing sea-ice.

Dead or Alive; or Does It Really Matter? Level of Congruency Between Trophic Modes in Total and Active Fungal Communities in High Arctic Soil.

Describing dynamics of belowground organisms, such as fungi, can be challenging. Results of studies based on environmental DNA (eDNA) may be biased as the template does not discriminate between metabolically active cells and dead biomass. We analyzed ribosomal DNA (rDNA) and ribosomal RNA (rRNA) coextracted from 48 soil samples collected from a manipulated snow depth experiment in two distinct vegetation types in Svalbard, in the High Arctic. Our main goal was to compare if the rDNA and rRNA metabarcoding templates produced congruent results that would lead to consistent ecological interpretation. Data derived from both rDNA and rRNA clustered according to vegetation types. Different sets of environmental variables explained the community composition based on the metabarcoding template. rDNA and rRNA-derived community composition of symbiotrophs and saprotrophs, unlike pathotrophs, clustered together in a similar way as when the community composition was analyzed using all OTUs in the study. Mean OTU richness was higher for rRNA, especially in symbiotrophs. The metabarcoding template was more important than vegetation type in explaining differences in richness. The proportion of symbiotrophic, saprotrophic and functionally unassigned reads differed between rDNA and rRNA, but showed similar trends. There was no evidence for increased snow depth influence on fungal community composition or richness. Our findings suggest that template choice may be especially important for estimating biodiversity, such as richness and relative abundances, especially in Helotiales and Agaricales, but not for inferring community composition. Differences in study results originating from rDNA or rRNA may directly impact the ecological conclusions of one's study, which could potentially lead to false conclusions on the dynamics of microbial communities in a rapidly changing Arctic.

Microbial Eukaryotes in an Arctic Under-Ice Spring Bloom North of Svalbard.

Microbial eukaryotes can play prominent roles in the Arctic marine ecosystem, but their diversity and variability is not well known in the ice-covered ecosystems. We determined the community composition of microbial eukaryotes in an Arctic under-ice spring bloom north of Svalbard using metabarcoding of DNA and RNA from the hypervariable V4 region of 18S nrDNA. At the two stations studied, the photosynthetic biomass was dominated by protists >3 µm and was concentrated in the upper 70-80 m, above the thermocline and halocline. Hierarchical cluster analyses as well as ordination analyses showed a distinct clustering of the microbial eukaryote communities according to a combination of water mass and local environmental characteristics. While samples collected in the surface mixed layer differed distinctly between the two sites, the deeper communities collected in Atlantic Water were fairly similar despite being geographically distant. The differentiation of the microbial eukaryote communities of the upper mixed water was probably driven by local development and advection, while the lack of such differentiation in the communities of Atlantic Water reflects the homogenizing effect of water currents on microbial communities.

Strong Seasonality of Marine Microbial Eukaryotes in a High-Arctic Fjord (Isfjorden, in West Spitsbergen, Norway).

The Adventfjorden time series station (IsA) in Isfjorden, West Spitsbergen, Norway, was sampled frequently from December 2011 to December 2012. The community composition of microbial eukaryotes (size, 0.45 to 10 µm) from a depth of 25 m was determined using 454 sequencing of the 18S V4 region amplified from both DNA and RNA. The compositional changes throughout the year were assessed in relation to in situ fjord environmental conditions. Size fractionation analyses of chlorophyll a showed that the photosynthetic biomass was dominated by small cells (<10 µm) most of the year but that larger cells dominated during the spring and summer. The winter and early-spring communities were more diverse than the spring and summer/autumn communities. Dinophyceae were predominant throughout the year. The Arctic Micromonas ecotype was abundant mostly in the early-bloom and fall periods, whereas heterotrophs, such as marine stramenopiles (MASTs), Picozoa, and the parasitoid marine alveolates (MALVs), displayed higher relative abundance in the winter than in other seasons. Our results emphasize the extreme seasonality of Arctic microbial eukaryotic communities driven by the light regime and nutrient availability but point to the necessity of a thorough knowledge of hydrography for full understanding of their succession and variability.

In vivo expression of a group I intron HEG from the antisense strand of Didymium ribosomal DNA.

Two different isolates of the myxomycete Didymium iridis harbour homing endonuclease genes that are expressed from group I introns inserted into identical sites within the small subunit ribosomal DNA. The homing endonuclease proteins are related in sequence, and their gene structures share similar features such as the presence of small spliceosomal introns and functional polyadenylation sites. However, they are transcribed from opposite strands of the ribosomal DNA and presumable by different RNA polymerases. We have previously described the in vivo expression of the I-DirI homing endonuclease from within the ribosomal RNA precursor. In this paper, we demonstrate the in vivo expression of the I-DirII homing endonuclease from the opposite strand of the Didymium rRNA gene. A comparison of the expression strategies of the two genes demonstrates the feasibility of antisense expression and provides insight into nucleolar gene expression.

The recent transfer of a homing endonuclease gene.

The myxomycete Didymium iridis (isolate Panama 2) contains a mobile group I intron named Dir.S956-1 after position 956 in the nuclear small subunit (SSU) rRNA gene. The intron is efficiently spread through homing by the intron-encoded homing endonuclease I-DirI. Homing endonuclease genes (HEGs) usually spread with their associated introns as a unit, but infrequently also spread independent of introns (or inteins). Clear examples of HEG mobility are however sparse. Here, we provide evidence for the transfer of a HEG into a group I intron named Dir.S956-2 that is inserted into the SSU rDNA of the Costa Rica 8 isolate of D.iridis. Similarities between intron sequences that flank the HEG and rDNA sequences that flank the intron (the homing endonuclease recognition sequence) suggest that the HEG invaded the intron during the recent evolution in a homing-like event. Dir.S956-2 is inserted into the same SSU site as Dir.S956-1. Remarkably, the two group I introns encode distantly related splicing ribozymes with phylogenetically related HEGs inserted on the opposite strands of different peripheral loop regions. The HEGs are both interrupted by small spliceosomal introns that must be removed during RNA maturation.

The group I-like ribozyme DiGIR1 mediates alternative processing of pre-rRNA transcripts in Didymium iridis.

During starvation induced encystment, cells of the myxomycete Didymium iridis accumulate a 7.5-kb RNA that is the result of alternative processing of pre-rRNA. The 5' end corresponds to an internal processing site cleaved by the group I-like ribozyme DiGIR1, located within the twin-ribozyme intron Dir.S956-1. The RNA retains the majority of Dir.S956-1 including the homing endonuclease gene and a small spliceosomal intron, the internal transcribed spacers ITS1 and ITS2, and the large subunit rRNA lacking its two group I introns. The formation of this RNA implies cleavage by DiGIR1 in a new RNA context, and presents a new example of the cost to the host of intron load. This is because the formation of the 7.5-kb RNA is incompatible with the formation of functional ribosomal RNA from the same transcript. In the formation of the 7.5-kb RNA, DiGIR1 catalysed cleavage takes place without prior splicing performed by DiGIR2. This contrasts with the processing order leading to mature rRNA and I-DirI mRNA in growing cells, suggesting an interplay between the two ribozymes of a twin-ribozyme intron.