Comparison of Two 16S rRNA Primers (V3–V4 and V4–V5) for Studies of Arctic Microbial Communities

Fadeev, Eduard; Cardozo-Mino, Magda G.; Rapp, Josephine Z.; Bienhold, Christina; Salter, Ian; Salman-Carvalho, Verena; Molari, Massimiliano; Tegetmeyer, Halina E.; Buttigieg, Pier Luigi; Boetius, Antje

doi:10.3389/fmicb.2021.637526

ORIGINAL RESEARCH article

Front. Microbiol., 16 February 2021

Sec. Aquatic Microbiology

Volume 12 - 2021 | https://doi.org/10.3389/fmicb.2021.637526

This article is part of the Research Topic Marine Microbiomes: Towards Standard Methods and Best Practices View all 10 articles

Comparison of Two 16S rRNA Primers (V3–V4 and V4–V5) for Studies of Arctic Microbial Communities

$\r\nEduard Fadeev,*&#x;$ Eduard Fadeev^1,2*†

Magda G. Cardozo-Mino^1,2

Josephine Z. Rapp³

Christina Bienhold^1,2

Ian Salter^1,4

Verena Salman-Carvalho⁵

Massimiliano Molari²

Halina E. Tegetmeyer^2,6

Pier Luigi Buttigieg^1,2

Antje Boetius^1,2,7

¹Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany
²Max Planck Institute for Marine Microbiology, Bremen, Germany
³School of Oceanography, University of Washington, Seattle, WA, United States
⁴Faroe Marine Research Institute, Tórshavn, Faroe Islands
⁵Department of Microbiology, Morrill Science Center IVN, University of Massachusetts, Amherst, MA, United States
⁶Center for Biotechnology, Bielefeld University, Bielefeld, Germany
⁷MARUM, University of Bremen, Bremen, Germany

Microbial communities of the Arctic Ocean are poorly characterized in comparison to other aquatic environments as to their horizontal, vertical, and temporal turnover. Yet, recent studies showed that the Arctic marine ecosystem harbors unique microbial community members that are adapted to harsh environmental conditions, such as near-freezing temperatures and extreme seasonality. The gene for the small ribosomal subunit (16S rRNA) is commonly used to study the taxonomic composition of microbial communities in their natural environment. Several primer sets for this marker gene have been extensively tested across various sample sets, but these typically originated from low-latitude environments. An explicit evaluation of primer-set performances in representing the microbial communities of the Arctic Ocean is currently lacking. To select a suitable primer set for studying microbiomes of various Arctic marine habitats (sea ice, surface water, marine snow, deep ocean basin, and deep-sea sediment), we have conducted a performance comparison between two widely used primer sets, targeting different hypervariable regions of the 16S rRNA gene (V3–V4 and V4–V5). We observed that both primer sets were highly similar in representing the total microbial community composition down to genus rank, which was also confirmed independently by subgroup-specific catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) counts. Each primer set revealed higher internal diversity within certain bacterial taxonomic groups (e.g., the class Bacteroidia by V3–V4, and the phylum Planctomycetes by V4–V5). However, the V4–V5 primer set provides concurrent coverage of the archaeal domain, a relevant component comprising 10–20% of the community in Arctic deep waters and the sediment. Although both primer sets perform similarly, we suggest the use of the V4–V5 primer set for the integration of both bacterial and archaeal community dynamics in the Arctic marine environment.

Introduction

The Arctic Ocean is the most rapidly changing marine region on the planet due to its fast warming causing substantial sea-ice loss (Peng and Meier, 2018; Dai et al., 2019), as well as increasing pollution (Peeken et al., 2018). To assess the impact of global climate change on marine food web dynamics and elemental cycles, it is important to monitor variations in microbial community structure with time (Karl and Church, 2014; Fuhrman et al., 2015; Buttigieg et al., 2018). However, the Arctic Ocean is generally under-sampled in ice-covered regions and in winter (Wassmann et al., 2011), particularly with regard to assessments of its microbial communities and their biogeochemical functions (Boetius et al., 2015). Until recently, microbial monitoring efforts in the deep Arctic Ocean consisted of 1 year-round long-term time series at the HAUSGARTEN observatory in the Fram Strait (Soltwedel et al., 2005, 2015), as well as a few other process studies (e.g., Kirchman et al., 2007; Alonso-Sáez et al., 2008; Nikrad et al., 2012; Wilson et al., 2017; Müller et al., 2018).

The Arctic Ocean features substantial vertical structure that may select for specific microbial types in the sea ice (Boetius et al., 2015; Rapp et al., 2018), in the ice-free and the ice-covered highly stratified surface waters (Wilson et al., 2017; Fadeev et al., 2018), in the sinking particles (further addressed as “marine snow”; Fadeev et al., 2020), as well as in the water and sediments of the deep-sea where temperatures are year-round close to freezing point temperatures (Bienhold et al., 2012; Hoffmann et al., 2017; Wilson et al., 2017; Rapp et al., 2018; Fadeev et al., 2020). Throughout the annual cycle, Arctic surface waters bacterial and archaeal communities exhibit pronounced fluctuations of the dominant taxonomic groups (Alonso-Sáez et al., 2008; Wilson et al., 2017; Müller et al., 2018), which are strongly associated with presence of sea ice and the seasonal phytoplankton blooms (Kirchman et al., 2007; Nikrad et al., 2012; Fadeev et al., 2018; Cardozo-Mino et al., 2020). In winter, as well as under ice-covered conditions, the communities are dominated by the bacterial classes Alphaproteobacteria (mainly the SAR11 clade), Dehalococcoidia (mainly SAR202 clade), and the archaeal class Nitrososphaeria (Alonso-Sáez et al., 2008; Wilson et al., 2017; Müller et al., 2018). In the summer, and under ice-free conditions, the communities are dominated by the bacterial classes Bacteroidia (mainly the order Flavobacteriales) and Gammaproteobacteria (mainly the orders Alteromonadales and Oceanospirillales; Wilson et al., 2017; Fadeev et al., 2018). During the summer, differences between ice-covered and ice-free communities also affect the microbial diversity of the deep ocean and the seafloor via alterations of microbial communities on marine snow (Fadeev et al., 2020).

In the framework of the FRAM Microbial Observatory (FRontiers in Arctic marine Monitoring), we are aiming to develop a standardized methodology for long-term observations of microbial communities in these highly diverse Arctic Ocean environments, which will be also comparable to other long-term microbial time series locations (e.g., HOT and BATS). Unlike other time series sites of the world, the ice-cover and the harsh conditions of the Arctic Ocean are limiting the accessibility of the sampling sites to the summer months. Sampling campaigns during the winter (when microbial biomass is low; Kirchman et al., 2007; Alonso-Sáez et al., 2008) are rare and have only recently been achieved using autonomous samplers with limited sampling capacities (Liu et al., 2020). Therefore, the unique conditions and the currently available technologies constrain year-round microbial observations to PCR-based approaches (i.e., 16S rRNA gene amplicon sequencing), which can be realized with low concentrations of DNA (Thomas et al., 2012). Metagenomics approaches suggest that the functional capacity of marine microbial communities is strongly linked to their taxonomic composition (Galand et al., 2018; McNichol et al., 2020). Thus, when supported by curated taxonomic databases (e.g., SILVA 16S rRNA gene reference; Quast et al., 2013), 16S rRNA gene amplicon sequencing provides an affordable high-throughput tool for addressing traditional community ecology questions, especially under the constrained sampling conditions of the Arctic marine environment.

A critical step in 16S rRNA gene sequencing studies is the selection of PCR primers for DNA amplification (Armougom, 2009; Wang and Qian, 2009). Throughout the years, many primer sets were designed for diversity studies of specific taxonomic groups (e.g., SAR11 clade; Apprill et al., 2015), and attempts have been made to develop a more universal 16S rRNA gene primer sets that could cover close to the entire diversity of a natural microbial community (e.g., Earth Microbiome Project; Caporaso et al., 2012; Gilbert et al., 2014). The development of primer sets for the amplification of 16S rRNA genes is conducted in silico using reference databases (e.g., Klindworth et al., 2013). The Arctic Ocean is the smallest and shallowest of all five oceans, representing 4% of the area and 1% of the volume of the global ocean. Nevertheless, it plays an important role in global processes that are strongly affected by the ongoing climatic changes and is considered relevant for several Earth System tipping points (Wassmann and Reigstad, 2011; Lenton et al., 2019). Furthermore, being the coldest among the oceans, with strong stratification and only limited deep-water exchange, the Arctic Ocean is likely to contain unique endemic microbial diversity that drives its biogeochemical cycles (Kirchman et al., 2009; Ghiglione et al., 2012; Pedrós-Alió et al., 2015). An example for such locally adapted Arctic diversity was recently found with Arctic specific members of the ubiquitous SAR11 clade (Kraemer et al., 2020). Furthermore, despite its global importance, sampling effort in the Arctic Ocean is low, especially in and under the sea ice and in the deep basin, as well as generally during the wintertime (Wassmann et al., 2011; Royo-Llonch et al., 2020). Thus, the reference databases are likely lacking proper coverage of the complexity and dynamics of the Arctic Ocean microbiomes that may result in biased representations of them by currently available 16S rRNA gene primers.

One of the most extensively used primer set for the investigation of bacterial diversity in various environments is the 341F/785R (targeting the V3–V4 hypervariable regions of the 16S rRNA gene) that was developed by Klindworth et al. (2013). For the investigation of marine microbiomes, an alternative primer set 515F-Y/926R (targeting the V4–V5 hypervariable regions of the 16S rRNA gene), which is also able to capture the diversity of the archaeal communities, has been developed by Parada et al. (2016). Currently, both V3–V4 and V4–V5 primer sets are widely used in studies of marine microbial communities and were extensively tested using mock and natural communities of temperate waters (e.g., Wear et al., 2018; Willis et al., 2019; McNichol et al., 2020). However, no study has systematically tested the performance of these primer sets on microbial communities of the Arctic Ocean.

In an attempt to select the most suitable primer set for the long-term monitoring of Arctic microbial communities as part of the FRAM Molecular Observatory, we present here a performance comparison of the 16S rRNA gene primer sets V3–V4 (341F/785R) and V4–V5 (515F-Y/926R). Our hypothesis was that due to relatively low representation of Arctic microbial communities in public databases (due to low number of existing studies), the 16S rRNA gene primer sets may capture different parts of microbial diversity in these unique environments. To test this hypothesis, we have conducted a direct comparison of the taxonomic coverage and potential biases of the two primer sets in 37 field samples collected from various environments of the Arctic Ocean, including sea-ice, surface and deep water column, marine snow, and deep-sea sediment. As an independent line of validation, we performed cell counting of five key taxonomic subgroups in a subset of the field samples via CARD-FISH (catalyzed reporter deposition-fluorescence in situ hybridization).

Materials and Methods

Sample Collection

The samples included in this study were collected at the long-term ecological research (LTER) site HAUSGARTEN in Fram Strait and the central Arctic Ocean (Supplementary Figure 1 and Supplementary Table 1). The samples were collected as follows:

• The sea-ice cores were collected using an ice corer (9 cm diameter; Kovacs Enterprise, Roseburg, OR, United States) and broken into subsections to facilitate quicker melting. The lower 30–50 cm of the sea ice (depending on total core length) was melted in plastic containers (rinsed with ethanol and ultrapure water) at 4°C in the dark. The melting of the sea ice took ∼24 h and the samples were immediately filtered on 0.22 μm Sterivex^TM membranes as soon as the last piece of sea ice melted. Additional samples for microscopy counts were filtered onto 0.22 μm polycarbonate membranes (Whatman Nucleopore, Buckinghamshire, United Kingdom), with sterile filtered formalin at a final concentration of 2% and stored at −20°C.

• The water sampling was carried out using 12 L Niskin bottles mounted on a CTD rosette (Sea-Bird Electronics Inc., SBE 911 plus probe, Bellevue, WA, United States) and filtered on 0.22 μm Sterivex^TM membranes. The Sterivex^TM membranes were then stored at −20°C until further processing. Additional samples for microscopy counts were filtered onto 0.22 μm polycarbonate membranes (Whatman Nucleopore, Buckinghamshire, United Kingdom), with sterile filtered formalin at a final concentration of 2% and stored at −20°C.

• The deep-sea sediment cores were retrieved by a TV-guided multicorer, and subsamples of the uppermost centimeter of the cores were collected with syringes and immediately stored at −20°C until further processing.

• The marine snow samples were collected using sediment traps of the long-term moorings at the LTER site HAUSGARTEN (Bauerfeind et al., 2009; Lalande et al., 2013). Collection cups (400 ml) were filled with filtered seawater, adjusted to a salinity of 40 and poisoned with HgCl₂ (0.14% final solution) to preserve samples during deployment and after recovery (Metfies et al., 2017). After recovery, samples were stored at +4°C, swimmers were removed and samples were split by a wet splitting procedure (Bodungen et al., 2013). In this study, we used 1/32 splits of the original trap sample. Sinking particles from the sediment trap samples were collected on 0.22 μm Sterivex filters and stored at −20°C.

All metadata of the samples are accessible via the Data Publisher for Earth and Environmental Science PANGAEA¹, the PANGAEA event IDs are listed in Supplementary Table 1. Sampling map was produced using Ocean Data View v5.2.1 (Schlitzer, 2018).

DNA Isolation and 16S rRNA Gene Amplicon Sequencing

Genomic DNA was isolated in a combined chemical and mechanical procedure using the PowerWater DNA Isolation Kit for sea ice, water, and sediment traps and using the PowerSoil DNA Isolation Kit for sediment samples (MO BIO Laboratories, Inc., Carlsbad, CA, United States). Prior to DNA isolation, the 0.22 μm Sterivex^TM membrane cartridges of the seawater and sea ice samples were cracked open in order to place the filters into the kit-supplied bead beating tubes. The isolation was continued according to the manufacturer’s instructions, and DNA was stored at −20°C. Library preparation was performed according to the standard instructions of the 16S Metagenomic Sequencing Library Preparation protocol (Illumina^TM, Inc., San Diego, CA, United States). Two different hypervariable regions of the bacterial 16S rRNA gene were amplified using aliquots of the isolated DNA from each sample. The V3–V4 region was amplified using the S-D-Bact-0341-b-S-17 (5′-CCTACGGGNGGCWGCAG-3′) and the S-D-Bact-0785-a-A-21 (5′-GACTACHVGGGTATCTAATCC-3′) primers (Klindworth et al., 2013). The V4–V5 regions was amplified using the 515F-Y (5′-GTGYCAGCMGCCGCGGTAA-3′) and the 926R (5′-CCGYCAATTYMTTTRAGTTT-3′) primers (Parada et al., 2016). Sequences were obtained on the Illumina MiSeq^TM platform in a 2 × 300 bp paired-end run and for surface water samples on the Illumina HiSeq^TM platform in a 2 × 250 bp paired-end run (CeBiTec, Bielefeld, Germany), following the standard instructions of the 16S Metagenomic Sequencing Library Preparation protocol.

Raw paired-end, primer-trimmed reads were deposited in the European Nucleotide Archive (ENA; Harrison et al., 2019) under accession number PRJEB31938. The data were archived using the brokerage service of the German Federation for Biological Data (GFBio; Diepenbroek et al., 2014).

Bioinformatics and Statistical Analyses

The raw paired-end reads were primer-trimmed using cutadapt (Martin, 2011). Further analyses were conducted using R v4.0.0² in RStudio v1.2.5042³. The libraries were processed using DADA2 v1.16 (Callahan et al., 2016a), following the suggested workflow (Callahan et al., 2016b). The reads in MiSeq libraries were truncated at 255 bp length for forward reads and at 200 bp length for reverse reads, to facilitate the technical quality drop at the end of the reads. Reads in both MiSeq and HiSeq were then trimmed for low-quality bases and merged based on a minimum overlap of 10 bp. Chimeras and amplicon sequence variants (ASVs) that were observed in only one sample were filtered out. The representative sequences were taxonomically classified against SILVA 16S rRNA gene reference database release 138 (Quast et al., 2013; Yilmaz et al., 2014). The ASVs that were taxonomically unclassified at phylum rank or were not assigned to bacterial or archaeal lineages were excluded from further analysis. Furthermore, all ASVs that were taxonomically assigned to mitochondria and chloroplast were removed from the dataset.

Sample data matrices were managed using the R package “phyloseq” v1.32 (McMurdie and Holmes, 2013), and plots were generated using the R package “ggplot2” v3.3.0 (Gómez-Rubio, 2017). The sample rarefaction analyses were conducted using the R package “iNEXT” v2.0.20 (Hsieh et al., 2016). To test the effect of the different primer sets on the taxonomic composition of the microbial communities, as well as to test for differences between microbial communities of different types of samples, a two-way permutation multivariate analysis of variance (“Two-way PERMANOVA”) of Jensen–Shannon Divergence distance matrix was conducted (using the function “adonis2” in the R package “vegan” v.2.5.6; Oksanen et al., 2007).

Scripts for processing data can be accessed at https://github.com/edfadeev/Arctic-16S-Primers-comparison/.

Catalyzed Reporter Deposition Fluorescence in situ Hybridization

Both sea ice and seawater samples were directly fixed in 4% formalin for 4 h at 4°C, filtered onto 0.22 μm polycarbonate Track-Etched^TM membranes (Whatman Nucleopore, Buckinghamshire, United Kingdom), and stored at −20°C. The CARD-FISH was applied based on the protocol established by Pernthaler et al. (2002), using horseradish-peroxidase (HRP)–labeled oligonucleotide probes (⁴ Ulm, Germany; Supplementary Table 4). All probes were checked for specificity and coverage of their target groups against the SILVA 16S rRNA gene reference. All filters were embedded in 0.2% low-gelling-point agarose and treated with 10 mg mL^–1 lysozyme solution (Sigma-Aldrich Chemie GmbH, Hamburg, Germany) for 1 h at 37°C. Subsequently, endogenous peroxidases were inactivated by submerging the filter pieces in 0.15% H₂O₂ in methanol for 30 min, before rinsing in Milli-Q water and dehydration in 96% ethanol. Then, the filters were covered in a hybridization buffer and a probe concentration of 0.2 ng μL^–1. Hybridization was performed at 46°C for 2.5 h, followed by washing in a pre-warmed washing buffer at 48°C for 10 min, and 15 min in 1x PBS. Signal amplification was carried out for 45 min at 46°C with an amplification buffer containing either tyramide-bound Alexa 488 (1 μg/mL^–1) or Alexa 594 (0.33 μg mL^–1). Afterward, the cells were counterstained in 1 μg/mL^–1 DAPI (4′,6-diamidino-2-phenylindole; Thermo Fisher Scientific GmbH, Bremen, Germany) for 10 min at 46°C. After rinsing with Milli-Q water and 96% ethanol, the filter pieces were embedded in a 4:1 mix of Citifluor (Citifluor Ltd., London, United Kingdom) and Vectashield (Vector Laboratories, Inc., Burlingame, United States) and stored overnight at −20°C for later microscopy evaluation.

Automated Image Acquisition and Cell Counting

The filters were evaluated microscopically under a Zeiss Axio Imager.Z2 stand (Carl Zeiss MicroImaging GmbH, Jena, Germany), equipped with a multipurpose fully automated microscope imaging system (MPISYS), a Colibri LED light source illumination system, and a multi-filter set 62HE (Carl Zeiss MicroImaging GmbH, Jena, Germany). Pictures were taken via a cooled charged-coupled-device (CCD) camera (AxioCam MRm; Carl Zeiss AG, Oberkochen, Germany) with a 63x oil objective, a numerical aperture of 1.4, and a pixel size of 0.1016 μm/pixel, coupled to the AxioVision SE64 Rel.4.9.1 software (Carl Zeiss AG, Oberkochen, Germany) as described by Bennke et al. (2016). Exposure times were adjusted after manual inspection with the AxioVision Rel.4.8 software coupled to the SamLoc 1.7 software (Zeder et al., 2011), which was also used to define the coordinates of the filters on the slides. For image acquisition, channels were defined with the MPISYS software, and a minimum of 55 fields of view with a minimum distance of 0.25 mm were acquired of each filter piece by recording a z-stack of seven images in autofocus.

Cell enumeration was performed with the software Automated Cell Measuring and Enumeration Tool (ACMETool3, 2018-11-09; M. Zeder, Technobiology GmbH, Buchrain, Switzerland). Cells were counted as objects according to manually defined parameters separately for the DAPI and FISH channels.

Results and Discussion

In this study, aliquots of 37 DNA samples from different environments in the Arctic Ocean (sea ice, surface and deep ocean water, marine snow, and seafloor sediment; Supplementary Table 1) were sequenced using two common primers sets that amplify either the V3–V4 or the V4–V5 hypervariable regions in the 16S rRNA gene and were subjected to the same bioinformatic workflow. Both primer sets showed a similar decrease in the number of sequences throughout the workflow, with 62 ± 13% and 68 ± 9% of sequences retained per sample, respectively. The final datasets consisted of 3,318,649 sequences in the V3–V4 dataset that were assigned to 12,045 ASVs and 3,340,628 sequences in the V4–V5 dataset that were assigned to 14,505 ASVs (Supplementary Table 2). In addition, the ASVs which were taxonomically assigned to eukaryotic, mitochondrial or chloroplast sequences, as well as ASVs unclassified at phylum rank, were also removed from further analysis (ca. 9% and ca. 17% of sequences in V3–V4 and V4–V5 datasets, respectively). In both datasets, an asymptotic extrapolation of the rarefaction curves did not further increase the number of observed ASVs (Supplementary Figure 2). Although, most likely further microbial diversity remains to be uncovered in all sampled environments, the rarefaction curves suggest that our samples contained most of the potential community richness covered by both primer sets. In sea ice, surface water (<30 m depth) and marine snow, both primer sets showed similar community richness (Figure 1). However, in the deep-water communities (>600 m depth), richness was significantly different between the primer sets (Wilcoxon Signed-Rank Test; p < 0.01), with ca. 40% more bacterial ASVs in the V3–V4. In contrast, the sediment community richness was significantly higher in the V4–V5 dataset (Wilcoxon signed-rank test; p < 0.01), with up to double the amount of bacterial ASVs compared to the V3–V4 dataset. The main taxonomic groups, typically observed in the Arctic marine environment, such as the classes Alphaproteobacteria, Bacteroidia, and Gammaproteobacteria, dominated both datasets (each comprising 10–30% of sequences in V3–V4 and V4–V4 datasets, respectively). However, within these groups significant differences between datasets in the number of observed ASVs were detected.

FIGURE 1

Figure 1. Chao1 richness estimates in the different sample types. Different primer sets represented by colors and shapes. Please note the differences of y-axis between the panels.

In the V3–V4 dataset, the Bacteroidia and Gammaproteobacteria showed the highest differences in number of observed ASVs within each class (i.e., type richness) compared to the V4–V5 dataset (Supplementary Table 2). The family Flavobacteriaceae (class Bacteroidia) comprised 18% of all sequences in both datasets; however, in the V3–V4 dataset, it consisted of one third more ASVs compared to the V4–V5 dataset (total of 278 and 196 ASVs, respectively; Figure 2). This difference in the number of observed ASVs was mainly associated with ASVs of the genus Polaribacter (total of 28 and 14 ASVs, respectively), a key heterotrophic bacterium that responds to phytoplankton blooms in mid- and high-latitudes (Gómez-Pereira et al., 2010; Fadeev et al., 2018; Avcı et al., 2020). The orders Alteromonadales, Cellvibrionales, and Oceanospirillales (all within the class Gammaproteobacteria), which comprised 4–6% of all sequences in the V3–V4 dataset and 3% of all sequences in the V4–V5 dataset, also showed differences between datasets in the number of observed ASVs (Supplementary Table 3). Each of these Gammaproteobacteria orders contained two times more ASVs in the V3–V4 dataset, compared to the V4–V5 dataset (the largest difference was in the order Alteromonadales, with total of 113 and 49 ASVs, respectively). These taxonomic groups are typically associated with organic matter degradation (Buchan et al., 2014), and were previously shown to dominate sea ice microbial communities associated with algal aggregates (Rapp et al., 2018), as well as surface waters during phytoplankton blooms (Fadeev et al., 2018). Furthermore, the family Woeseiaceae (class Gammaproteobacteria) also consisted of ca. 30% more ASVs in the V3–V4 dataset, compared to the V4–V5 dataset (total of 127 and 98 ASVs, respectively; Figure 2). This bacterial family is abundant in deep-sea sediments around the globe, including the Arctic Ocean (Bienhold et al., 2016; Hoffmann et al., 2020).

FIGURE 2

Figure 2. Major taxonomic families in V3–V4 and V4–V5 datasets. The x-axis represents the total sequence proportion of each family in V3–V4 (left panel) and V4–V5 (right panel) datasets. The numbers at each column represent the number of observed ASVs affiliated with each taxonomic family. Different taxonomic classes are represented by color code. Only families that comprised at least 1% of sequences in at least one of the datasets were included in the visualization.

Compared to the V3–V4 dataset, the V4–V5 dataset consisted of at least one third more ASVs in the classes Phycisphaerae (total of 206 and 117 ASVs, respectively) and Planctomycetes (total of 299 and 244 ASVs, respectively). This difference in the number of observed ASVs was mainly associated with the families Pirellulaceae that comprised ca. 2% of all sequences in both datasets (Figure 2), as well as Phycisphaeraceae that comprised less than 1% of all sequences (Supplementary Table 3) in both datasets. These taxonomic groups have been previously shown to be associated with sinking particles in the deep ocean and are also abundant in Arctic deep-sea sediments (Fadeev et al., 2020). Furthermore, the archaeal class Nitrososphaeria was almost absent from the V3–V4 dataset, with only a few sequences associated with four ASVs, compared to 168 ASVs in the V4–V5 dataset that comprised 7% of the total sequences (Figure 2). Marine members of the Archaea in general, and the class Nitrososphaeria in particular, are abundant in the Arctic marine environment and can reach up to one fifth of the cells in Arctic microbial communities (Müller et al., 2018; Cardozo-Mino et al., 2020). Taken together, these observations suggest that on ASV level the diversity of different taxonomic groups are captured differently by the two primer sets. This is potentially a result of differences in the regional hypervariability of the 16S rRNA gene within different taxonomic groups (Yang et al., 2016; Kerrigan et al., 2019). In addition, as was previously shown for various taxonomic groups, such as the SAR11 clade, differences in captured diversity may occur also due to specificity differences of the primer sets to the targeted 16S rRNA gene region (Parada et al., 2016).

Despite the observed differences on an ASV level, the overall taxonomic composition was consistent between the datasets (Figure 3). Sampled sea ice, surface water, and marine snow communities were dominated by heterotrophic bacteria of the classes Bacteroidia (mainly the genus Polaribacter) and Gammaproteobacteria (mainly the genera in the order Alteromonadales), with equivalent relative sequence abundances to those described in previous reports (Bowman et al., 2012; Eronen-Rasimus et al., 2016; Hatam et al., 2016; Wilson et al., 2017; Fadeev et al., 2018, 2020; Rapp et al., 2018). At depth, pelagic communities were dominated by sequences of the class Alphaproteobacteria, SAR324 clade, and the archaeal class Nitrososphaeria, all of which were previously observed to dominate deep Arctic waters, as well as surface communities during the Arctic winter (Wilson et al., 2017; Fadeev et al., 2020). The sediment communities, which have previously been shown to harbor the highest taxonomic diversity among the described Arctic environments by far (Bienhold et al., 2012; Hoffmann et al., 2017; Rapp et al., 2018), were dominated in sequence abundance of Gammaproteobacteria.

FIGURE 3

Figure 3. Taxonomic compositions of the microbial communities. Different taxonomic classes are represented by color code. Classes with sequence proportion below 2% were classified as “Other taxa”.

In order to compare the differences in representation of taxonomic groups between the primer sets, we combined sequence abundances of all ASVs according to their taxonomic affiliation at genus rank (i.e., the highest possible shared between the datasets taxonomic resolution). In the V3–V4 dataset, the ASVs were merged into 306 different genera and 279 lineages that were affiliated to higher taxonomic ranks (i.e., were unclassified on a genus rank). In the V4–V5 dataset, the ASVs were merged into 280 different genera and 299 lineages that were affiliated to higher taxonomic ranks. Overall, 489 (72% of the total) lineages were observed in both datasets at this level of taxonomic resolution. In the V3–V4 dataset there were 96 (14% of the total) lineages that were absent from the V4–V5 dataset, but together they comprised less than 1% of the sequences in the V3–V4 dataset. On the other hand, in the V4–V5 dataset there were 90 (13% of the total) lineages that were absent from the V3–V4 dataset, and together they comprised 5% of the sequences in the V4–V5 dataset. In addition, the dissimilarity of community compositions in merged V3–V4 and V4–V5 datasets revealed consistent and significant difference between the microbiomes captured by both primer sets (Two-way PERMANOVA test; F_4,64 = 86.29, R² = 0.83, p value < 0.001; Figure 4). Only a small fraction of the total variance was associated with the difference between the primer sets (Two-way PERMANOVA test; F_1,64 = 7.59, R² = 0.02, p value < 0.001). No significant combined effect of different primer sets on different sample types was observed (Two-way PERMANOVA test; p value > 0.05). Taken together, these results confirm that, even though the primer sets showed different sensitivity to diversity at the ASV level, both of them reflect similar taxonomic composition down to the genus rank.

FIGURE 4

Figure 4. Non-metric Multi-dimensional Scaling (NMDS) of the microbial communities in merged V3–V4 and V4–V5 datasets, based on Jensen-Shannon Divergence. The different types of samples are represented by colors, and the primer set are represented by shapes. Ellipses encompass clustering of each microbiome type with normal confidence of 0.95.

The research at the FRAM Microbial Observatory is focused on the seasonal and interannual dynamics of the Arctic Ocean associated with changes in sea ice extent and primary production in the surface ocean (e.g., Metfies et al., 2017; Fadeev et al., 2018, 2020). To further evaluate the performance of the two primer sets in these long-term monitored environments, we compared the sequence representation of selected taxonomic groups, which are associated with distinct stages of seasonal dynamics (Fadeev et al., 2018), to microscopically counted cells using CARD-FISH combined with an automated image acquisition (Cardozo-Mino et al., 2020). Fluorescence in situ hybridization techniques have the advantage of providing absolute abundances of (viable) cells that can be directly compared between samples. In microscopy counts of both sea ice and surface water communities the highest observed cell abundance was of the class Bacteroidia (up to 35 and 26% of the total microbial community, respectively), which was consistent with the representation of this taxonomic group by both primer sets. In surface water communities, high levels of consistency between the microscopy counts and both primer sets were observed also in the representation of Alteromonadales and Polaribacter (Table 1). On the other hand, the representation, by both primer sets, of the class Gammaproteobacteria in sea ice and surface water communities was 2–4 times higher in comparison to the proportion observed in microscopy counts (up to 9 and 18%, respectively). In contrast, the proportional abundance of the SAR11 clade was 5–10 times higher in the microscopy counts, compared to its representation by both primer sets (Table 1). Our results suggest that at least for some taxonomic groups (i.e., Polaribacter), both primer sets may provide a consistent semi-quantitative representation. However, microscopy results must be interpreted under several methodological caveats, knowing that low cellular ribosome content or low efficiency of the probe may alter the representation of individual taxa in our cell counts (Amann and Fuchs, 2008). Therefore, the observed inconsistency in the representation of some taxonomic groups (i.e., SAR11 clade) may also result from these limitations. In order to further investigate the quantitative performance of the primer sets, further investigation, using techniques such as mock communities (Yeh et al., 2019) or metagenomics (McNichol et al., 2020), is required.

TABLE 1

Table 1. Overview of cell abundances and sequence proportions range in selected taxonomic groups.

Concluding Remarks

To understand the links between the rapid environmental changes in the Arctic region and the dynamics of microbial communities in the Arctic Ocean, there is a need for robust methods addressing changes in diversity and relative abundance. In order to conduct such observations using a 16S rRNA gene tag-sequencing approach, optimally a similar extraction method and a single PCR primer set should be selected, which can be applied to all environments of the Arctic Ocean (sea ice, water column, and deep-sea sediment). The most suitable primer set for 16S rRNA amplification and sequencing from environmental samples should produce high-quality amplicon libraries and cover with minimum biases the variety of present organisms, as well as their relative abundances. We have found that at all taxonomic ranks down to genus, both primer sets represent the overall richness of the major bacterial taxonomic groups at comparable levels across the different Arctic Ocean biomes. The relative sequence abundance of some dominant taxonomic groups, such as the Polaribacter, corresponds with their proportional representation via microscopic cell counts. Other taxonomic groups such as the SAR11 clade strongly differ between the molecular and the microscopical representations. However, this discrepancy may be due to limitations of the microscopical quantification. On an ASV level, both primer sets capture the diversity within the most abundant taxonomic groups differently, and thus, the use of each primer set may depend on the target groups. However, the main advantage of the V4–V5 primer set is its additional coverage of the archaeal domain, without compromising the detection of other taxonomic groups. Members of the Archaea comprise a substantial fraction of Arctic marine microbial communities, particularly during the dark season and in deep waters. Thus, given the demonstrated similarities and differences, we endorse the use of the V4–V5 primer set for capturing comprehensive insights into microbial community dynamics of the Arctic marine environment.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

EF and AB designed the study. CB, IS, MM, and JR provided the environmental samples for the study. HT conducted the sequencing of the samples. MC-M and VS-C conducted the CARD-FISH counts. EF analyzed the data and wrote the manuscript. All authors contributed to the final version of the manuscript.

Funding

This work was conducted in the framework of the HGF Infrastructure Program FRAM of the Alfred-Wegener-Institute Helmholtz Center for Polar and Marine Research. The project has received funding from the European Research Council (ERC) under the European Union’s Seventh Framework Program (FP7/2007-2013) research project ABYSS (grant agreement no. 294757) awarded to AB. Additional funding came from the Austrian Science Fund (FWF) M-2797 to EF. This publication is Eprint ID 53010 of the Alfred Wegener Institute for Polar and Marine Research, Bremerhaven, Germany.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Daniel Sher for a critical review of the manuscript. We would also like to thank Jana Bäger, Jakob Barz, and Theresa Hargesheimer for DNA extractions and library preparations.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.637526/full#supplementary-material

Footnotes

References

Alonso-Sáez, L., Sánchez, O., Gasol, J. M., Balagué, V., and Pedrós-Alio, C. (2008). Winter-to-summer changes in the composition and single-cell activity of near-surface Arctic prokaryotes. Environ. Microbiol. 10, 2444–2454. doi: 10.1111/j.1462-2920.2008.01674.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Amann, R., and Fuchs, B. M. (2008). Single-cell identification in microbial communities by improved fluorescence in situ hybridization techniques. Nat. Rev. Microbiol. 6, 339–348. doi: 10.1038/nrmicro1888

PubMed Abstract | CrossRef Full Text | Google Scholar

Apprill, A., McNally, S., Parsons, R., and Weber, L. (2015). Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 75, 129–137. doi: 10.3354/ame01753