Exploring the Microdiversity Within Marine Bacterial Taxa: Toward an Integrated Biogeography in the Southern Ocean

Most of the microbial biogeographic patterns in the oceans have been depicted at the whole community level, leaving out finer taxonomic resolution (i.e., microdiversity) that is crucial to conduct intra-population phylogeographic study, as commonly done for macroorganisms. Here, we present a new approach to unravel the bacterial phylogeographic patterns combining community-wide survey by 16S rRNA gene metabarcoding and intra-species resolution through the oligotyping method, allowing robust estimations of genetic and phylogeographic indices, and migration parameters. As a proof-of-concept, we focused on the bacterial genus Spirochaeta across three distant biogeographic provinces of the Southern Ocean; maritime Antarctica, sub-Antarctic Islands, and Patagonia. Each targeted Spirochaeta operational taxonomic units were characterized by a substantial intrapopulation microdiversity, and significant genetic differentiation and phylogeographic structure among the three provinces. Gene flow estimations among Spirochaeta populations support the role of the Antarctic Polar Front as a biogeographic barrier to bacterial dispersal between Antarctic and sub-Antarctic provinces. Conversely, the Antarctic Circumpolar Current appears as the main driver of gene flow, connecting sub-Antarctic Islands with Patagonia and maritime Antarctica. Additionally, historical processes (drift and dispersal limitation) govern up to 86% of the spatial turnover among Spirochaeta populations. Overall, our approach bridges the gap between microbial and macrobial ecology by revealing strong congruency with macroorganisms distribution patterns at the populational level, shaped by the same oceanographic structures and ecological processes.

Introduction

Biogeography has traditionally investigated the geographic distribution of macroorganisms in the Eukaryota domain. Yet, the last two decades have witnessed a growing number of studies focused on the biogeography of microorganisms, and repeatedly reporting non-random community assemblages of various prokaryotic microorganisms (Martiny et al., 2006; Wilkins et al., 2013; Karimi et al., 2018). Unlike contemporary driving factors (i.e., environmental selection) that have been extensively studied (Gilbert et al., 2012; Stegen et al., 2012; van der Gast, 2015), the role of historical processes – past ecological and evolutionary events – onto the present-day distribution patterns of microorganisms remains poorly investigated. Initially, the consensus was that the rapid and widespread dispersal of microbes should erase any signal of past historical events (Martiny et al., 2006). Nevertheless, it is now clear that historical processes, such as the dispersal barriers and geographic distance, might substantially contribute to microbes’ biogeography instead of environmental filtering (Hanson et al., 2012, 2019). Similarly to larger organisms, biogeographic regionalization, isolation, and endemism have been reported for microbes, reflecting the predominant effect of geographic distance over environmental variations (Papke et al., 2003; Whitaker et al., 2003).

To date, most of the microbial biogeographic patterns have been depicted at the whole community level, taking advantage of the extensive survey capacity provided by the next-generation sequencing (NGS) technologies (Almasia et al., 2016; Karimi et al., 2018; Wu et al., 2019). Nevertheless, as observed in various empirical studies, a finer taxonomic scale generally allows better detection of geographic patterns (Hanson et al., 2012; Chase and Martiny, 2018; Bay et al., 2020). Moreover, the ecological processes driving the biogeographic patterns at the community-level would mainly result from the accumulation of micro-evolutive processes, i.e., mechanisms contributing to the genetic composition and diversity within populations, and how they vary in space and time (Hanson et al., 2012; Larkin and Martiny, 2017). Hence the comprehensive description of these micro-evolutive processes requires considering the intra-population diversity, as commonly applied in phylogeographic studies of macroorganisms. In other words, microbial assembly processes need to be investigated at a finer taxonomic resolution than usually done by microbial biogeographic surveys and consider the “microdiversity” within groups (Larkin and Martiny, 2017; Chase and Martiny, 2018).

The oceans have been considered among the most challenging environments to test hypotheses about microbial biogeography, mainly due to the speculated transport of organisms over large distances by marine currents and the absence of perceivable marine barriers impeding potential dispersal events (Fenchel and Finlay, 2004). However, oceanic fronts separating different water masses have been recently identified as major microbial dispersal barriers (Flaviani et al., 2018). The Southern Ocean (SO) is a vast region representing approximately 20% of the world ocean surface. It surrounds Antarctica, and its northern limit is the Subtropical Front (Gruber et al., 2019). Two main oceanographic structures shape the SO biogeography; the Antarctic Polar Front (APF) and the Antarctic Circumpolar Current (ACC). The APF is classically considered a harsh north–south obstacle for dispersing marine organisms due to the brutal change in water temperature and salinity (Griffiths, 2010; Halanych and Mahon, 2018). Phylogenetic reconstruction achieved on various vertebrate and invertebrate taxa clearly supports the role of the APF on their respective diversification processes (Díaz et al., 2011; Poulin et al., 2014; Hüne et al., 2015; González−Wevar et al., 2017; González-Wevar et al., 2019). Accordingly, and based on the described distribution of species, biogeographers have traditionally recognized Antarctica and sub-Antarctica as the two main biogeographic provinces of the SO, even if several provinces have been proposed within each of them (Koubbi et al., 2014). Contrarily, outside the APF, the ACC is generally described as the driver of genetic connection across the sub-Antarctic zone due to its eastward circulation (Nikula et al., 2010; Cumming et al., 2014; Moon et al., 2017; González−Wevar et al., 2018). Intraspecific genetic and phylogeographic studies of macroorganisms have demonstrated the ACC’s role in connecting geographically distant sub-Antarctic provinces (Gérard et al., 2008; Fraser et al., 2009; González−Wevar et al., 2018; Frugone et al., 2019). The marine biota distribution in this region has been synthesized in the Biogeographic Atlas of the SO, providing updated biogeographic information of a wide range of benthic and pelagic taxa from Metazoan, macroalgae, and phytoplankton (De Broyer et al., 2014). Despite being the most abundant and diverse domains on Earth, Bacteria, and Archaea are not included in the SO Atlas (Shade et al., 2018). Studies conducted at the whole community-level support (1) the role of ACC as a likely efficient mechanism of circumpolar microbial transport and dispersal (Murray and Grzymski, 2007; Wilkins et al., 2013) and (2) the role of APF as the main dispersal barrier separating Antarctic and sub-Antarctic microbial assemblages (Wilkins et al., 2013; Flaviani et al., 2018; Raes et al., 2018). However, even when geographic distributions of marine microbial communities have been characterized in the region, the underlying evolutionary processes remain unclear, and their comprehensive understanding may rely on a higher taxonomic perspective, exploring bacterial populations’ microdiversity.

Targeting the intraspecific microdiversity using NGS data requires specific computational methods to discriminate the stochastic noise caused by random sequencing errors from those associated with biologically significant diversity (Hanson et al., 2012; Edgar and Flyvbjerg, 2015; Chase and Martiny, 2018). For this purpose, an algorithm called “Minimum Entropy Decomposition” (MED) relying on the oligotyping method has been proposed by Eren et al. (2015). This algorithm allows to identify true sequence variants (i.e., oligotypes) within the “operational taxonomic units” (OTUs), classically defined at 97% identity of the bacterial 16S rRNA gene. This approach has already been successfully used to unravel fine-grained biogeographic patterns of bacterial microdiversity in Arctic sediments, such as variations in oligotype distribution according to spatial and environmental parameters (Buttigieg and Ramette, 2014). Moreover, focusing on the sulfate-reducing genus Desulfotomaculum in Arctic marine sediments, Hanson et al. (2019) showed clear biogeographic patterns – attributed to historical factors associated with past environments – were only evident at the microdiversity level achieved with the oligotyping method. However, the microevolutionary processes causing the microdiversity were not assayed, and the study did not encompass large-scale distribution among different biogeographic provinces, as it was restricted to the west coast of Spitsbergen, Svalbard in the Arctic Ocean.

In the present proof-of-concept study, we aim to elucidate the evolutionary processes driving microbial biogeography across different provinces of the SO by combining (1) community-wide surveying provided by the high-throughput sequencing of the 16S rRNA gene, (2) intra-species microdiversity resolution obtained through the oligotyping method implemented in the MED pipeline, and (3) phylogeographic analysis as traditionally developed for macroorganisms as models. Considering the SO as an outstanding idoneous frame, we investigated the geographic distribution of genetic diversity of marine bacterial taxa across three biogeographic provinces: maritime Antarctica (King George Island, South Shetland Islands, West Antarctic Peninsula), sub-Antarctic Islands of the Indian Ocean (Kerguelen archipelago), and southern South America (Patagonia). We selected sites separated by the APF, i.e., maritime Antarctica and Patagonia, and others connected through the ACC, i.e., Patagonia and Kerguelen archipelago.

As the contribution of geography to biological diversity patterns (i.e., dispersal limitation) is stronger on habitat-specialists (i.e., taxa found in habitat with high selective strength) (Lindström and Langenheder, 2012; Szekely et al., 2013), and emphasized within homogeneous habitats distributed across large spatial scales (Astorga et al., 2012; Hanson et al., 2012; Langenheder and Lindström, 2019), we focused our study on the bacterial community associated to a specific habitat: the gut of Abatus irregular sea urchins. The Abatus genus is distributed across the SO and gathers various sibling species homologous in ecology and habitat, such as Abatus cavernosus in southern South America, Abatus cordatus in Kerguelen Islands, and Abatus agassizii in maritime Antarctica (Poulin and Feral, 1995; David et al., 2005; Díaz et al., 2012; Guillaumot et al., 2018). Since these species lack specialized respiratory structure, they are restricted to the well-oxygenated coarse sediments found at shallow depth (typically 1–3 m depth) in sheltered bays, protected from the swell (Poulin and Feral, 1995). Within the Abatus hosts, we focused on a specific micro-environment – the gut tissue – previously described to act as a selective filter of the external sediment microbiota, as illustrated by the reduction of bacterial diversity at both taxonomic and functional levels (Schwob et al., 2020). Working on the gut community with supposedly more limited dispersal capacity, through a high sequencing depth, is expected to (1) provide robust coverage of the bacterial diversity, (2) minimize the relevance of environmental filtering between provinces, (3) emphasize the contribution of geographic and oceanographic factors, and therefore (4) enhance the detection of phylogeographic signals across the SO (Hanson et al., 2012). As a model taxon to explore bacterial phylogeography in the SO, we chose the Spirochaeta genus (phylum Spirochaetes). Spirochaeta bacteria are recognized as the most prevalent and abundant genus in the Abatus gut tissue (Schwob et al., 2020). Moreover, spirochaetes are classically reported in marine benthic sediments (Bowman and McCuaig, 2003; Beiruti et al., 2017) and, to a lesser extent, in the water column (Ocean Barcode Atlas¹). Thus, due to its ease of detection and ubiquity across biogeographic provinces, Spirochaeta represents an illustrative model to validate our methodology and explore marine bacteria’s spatial genetic patterns, from genus to population level. We hypothesized that the strong biogeographic barrier between South America and maritime Antarctica classically observed in the literature for macroorganisms (i.e., vicariance process) also affects the fine-scale genetic structure and the phylogeographic patterns within Spirochaeta OTUs. Besides, the ACC-mediated connectivity among sub-Antarctic provinces should be reflected by a greater genetic homogeneity of Spirochaeta populations between South American sites and the Kerguelen Islands, rather than with maritime Antarctica. Alternatively, the potential high dispersal capacity of Spirochaeta taxa may result in the absence of genetic and phylogeographic structure across the SO.

Materials and Methods

Sampling Collection, DNA Extraction, and 16S rRNA Gene Library Preparation

Adult Abatus individuals were sampled from four localities across the SO, including two sites in Patagonia, southern South America (Possession Bay, PAT1 and Puerto Deseado, PAT2), one site in Kerguelen Islands (Port-aux-Français, KER), and one site in the West Antarctic Peninsula (King George Island, KGI) (Figure 1 and Table 1). Marine surface sediments (0–5 cm, referred here as “external sediment”) were also sampled in each Abatus population’s immediate vicinity as the ingested food source of the sea urchins. Due to logistic constraints, it was not possible to collect external sediment in the PAT2 site. All individuals were dissected under sterile conditions to collect the whole digestive tract minus the caecum (identified as “gut tissue”) following Schwob et al. (2020). Gut tissue samples were gently rinsed with nuclease-free sterile water to remove the content (i.e., in sediment) and were then individually homogenized using mortar and pestle under a laminar-flow cabinet. Genomic DNA was extracted from both external sediments and the totality of the homogenized gut tissue samples using the DNeasy PowerSoil^® Kit (Qiagen, Hilden, Germany) following the manufacturer’s recommendations.

FIGURE 1

Figure 1. Sampling localities across the Southern Ocean, encompassing Possession Bay and Puerto Deseado in Atlantic Patagonia (PAT1 and PAT2, respectively), King George Island in Maritime Antarctica (KGI), and Port-aux-Français in Kerguelen Islands (KER). The Antarctic Polar Front (APF) and the Southern Antarctic Circumpolar Front (SACCF) are represented by the limit between the white and the blue areas, and by the dashed line, respectively.

TABLE 1

Table 1. Experimental design and sequencing data.

A metabarcoding approach was used to assess the bacterial community composition in the external sediment and Abatus gut tissue samples. Briefly, extracted genomic DNA was used as the template for PCR amplification using the primers 515F 5′-GTGYCAGCMGCCGCGGTA-3′ and 926R 5′-CCCCGYCAATTCMTTTRAGT-3′ (Parada et al., 2016). The PCR conditions and 16S rRNA gene library preparation were the same as described in Schwob et al. (2020).

Metabarcoding Data Processing

External sediment and gut tissue amplicons were sequenced using the paired-end sequencing technology (2 × 250 bp) on the Illumina MiSeq platform at the UWBC DNA Sequencing Facility (University of Wisconsin–Madison, United States). Reads of 16S rRNA were processed using the open-source software Mothur v1.44.0. Briefly, 3′ and 5′ reads were paired and trimmed according to their length and quality as described in Schwob et al. (2020). Chimeric sequences were removed using Uchime implemented in Mothur (Edgar et al., 2011). Reads were clustered into OTUs at 97% identity threshold, and a filter of relative abundance at >0.0001% was applied, as recommended by Bokulich et al. (2013). Following this, a taxonomic classification was performed with the classify.otu function and the SILVA database v138 implemented in Mothur. An OTU table of Spirochaeta was edited (i.e., all OTUs assigned to the genus Spirochaeta), and converted into a presence/absence matrix to detect the shared or exclusive OTUs across the four selected sites, discarding abundance variations resulting from short-term environmental conditions. Bray–Curtis and Unweighted Unifrac distances were calculated from the OTU presence/absence matrix and used to perform a non-metric multidimensional scaling (NMDS) with the metaMDS function of the ade4 package (Chessel et al., 2009). The NMDS was plotted through a scatter diagram using the s.class function implemented in ade4. The locality’s contribution to the Spirochaeta OTUs composition in gut tissue samples was tested with permutational multivariate analysis of variance (permanova), using adonis and pairwise.adonis functions implemented in vegan and pairwiseAdonis R packages, respectively (Oksanen et al., 2011; Martinez Arbizu, 2017).

A subset of the three most abundant Spirochaeta OTUs present in the four localities was retained for further analysis (Supplementary File 2). All the sequences assigned to the selected Spirochaeta OTUs were retrieved using the bin.seqs command in Mothur. Finally, the resulting fasta files were processed independently through the Minimum Entropy Decomposition pipeline following Eren et al. (2015).

Minimum Entropy Decomposition

Minimum Entropy Decomposition pipeline was used to identify nucleotidic polymorphism at fine-scale resolution (<3% of differences) within the 16S rRNA gene sequences from Spirochaeta OTUs. Briefly, MED employs the Shannon entropy algorithm to discriminate biologically meaningful variations of closely related sequences from the stochastic noise caused by random sequencing errors, focusing on informative-rich variable nucleotide positions (Eren et al., 2013, 2015). The resulting taxonomic units will be referred to as Spirochaeta oligotypes. Unsupervised oligotyping was carried out individually on Spirochaeta OTUs using the default dynamically computed threshold from which entropy is considered as zero (−m). Additionally, each identified oligotype had to have a default minimum relative abundance of 2% in the OTU sequences dataset (Eren et al., 2015). Accumulation curves of oligotypes’ richness were computed for each Spirochaeta OTU at a 95% confidence interval using the package iNext (Hsieh et al., 2016) in r v3.6.0 (R Core Team, 2012). Pie charts of the relative site contributions in the total abundance of the Spirochaeta OTUs oligotypes were performed with the pie function in the package graphics in r v3.6.0.

Genetic Diversity and Structure of Spirochaeta Populations

The number of oligotypes (k), the oligotype diversity (H), the number of discriminant sites (S) and the pairwise difference between sequences (Φ) were estimated individually for each Spirochaeta OTU using the packages pegas (Paradis, 2010) and ape v5.3.0 (Paradis et al., 2015) in R v3.6.0. For comparative purposes among sites with unequal sample sizes, a composite bootstrapping script was written to rarefy the sequence datasets to the minimum number of sequences per site and repeat the rarefaction with 1,000 re-samplings. Confidence intervals at 95% of genetic diversity indices were then calculated using these iteration values. The genetic differentiation (F_st and Φ_st) among Spirochaeta populations was analyzed using the software arlequin v3.5.2 (Schneider et al., 2000) with 1,000 permutations and a significance threshold at 0.05. Phylogeographic differentiation was also estimated with the nearest-neighbor statistic Snn (Hudson, 2000), and the significance of Snn estimates was tested with a permutation test through DnaSP v5.10.01 (Librado and Rozas, 2009). The reconstruction of oligotype networks was performed using the Median Joining method with the software Populational Analysis with Reticulate Trees v1.7.0 in PopART (Leigh and Bryant, 2015). Because of the differences in Spirochaeta OTUs’ abundances across the four sites, oligotypes frequencies were calculated for each site and then summed per OTU to reconstruct the oligotype networks.

Quantification of Selection, Dispersal, and Drift

The relative contribution of stochastic (i.e., dispersal, drift) and deterministic (i.e., selection) processes, on Spirochaeta oligotype assembly was measured for the selected OTU, following the analytical framework described in Stegen et al. (2015) and illustrated by Feng et al. (2018). Briefly, the approach relies on the comparison of the phylogenetic turnover between communities across samples (β mean nearest-taxon distance, βMNTD) to a null distribution of βMNTD, and denoted as the β-nearest taxon index (βNTI). βNTI values indicate that taxa between two communities are more (i.e., βNTI < −2) or less (i.e., βNTI > +2) phylogenetically related than expected by chance, thus suggesting that communities experience homogenizing or variable selection, respectively (Stegen et al., 2015). βNTI values ranging from −2 to +2 indicate a limited selection effect and point to dispersal limitation and ecological drift out as possible community composition drivers.

Then, the respective effect of dispersal limitation and ecological drift were disentangled using the pairwise Bray–Curtis-based Raup–Crick dissimilarity index (RC_Bray) among sites (Chase et al., 2011), weighted by oligotype abundance (Stegen et al., 2013). The RC_Bray values < −0.95 and > +0.95 correspond to communities that have –respectively– more or fewer taxa in common than expected by chance, and therefore indicate that community turnover is driven by homogenizing dispersal (RC_Bray < −0.95) or dispersal limitation plus drift (RC_Bray > +0.95). On the contrary, RC_Bray values > −0.95 and < +0.95 are indicative of ecological drift (Ramoneda et al., 2020).

The βMNTD/βNTI and RC_Bray matrices, and the respective contributions of the four ecological processes, were calculated using an optimized version of the initial script of Stegen et al. (2013), developed by Richter-Heitmann et al. (2020). The phylogenetic trees required for the βMNTD/βNTI matrix were generated using PhyML v3.0 (Guindon et al., 2010), and the oligotype sequences of Spirochaeta previously aligned with MUSCLE (Edgar, 2004). Comparisons with βNTI and RC_Bray null models included 999 randomizations.

Testing for Isolation by Distance and Environment

To disentangle the relative effect of geographic distance and abiotic environmental differences on the Spirochaeta oligotype composition between samples, we used the distance-based multiple matrix regression with randomization (MMRR) approach (Wang, 2013). We extracted a set of 9 environmental variables for each of our sampling site from the Bio-ORACLE database (Assis et al., 2018), including pH, the means of nitrate, silicate, and phosphate concentrations, and the means at the mean depth of seawater salinity, dissolved oxygen concentration, seawater temperature, seawater temperature range and chlorophyll concentration. All environmental variables were standardized (($x_i-\overline{x})/sd(x))$, and were then analyzed through principal components analysis (PCA). As a high percentage of the variation among localities was explained by the first component (PC1, >91%, Supplementary Files 6, 7), we transformed the scores of PC1 into Euclidean distance using the vegdist function in the vegan package in R to use it as the environmental distance matrix further. The longitude and latitude coordinates were converted into kilometers using the earth.dist function implemented in the fossil package in R (Vavrek, 2011). The geographic distances were transformed using the Hellinger method through the decostand function of the vegan package in R. The dissimilarity matrix of Spirochaeta oligotype composition among samples was also created from Bray–Curtis distances using the vegdist function of the R package vegan. Finally, to evaluate the relative weight of environmental and geographic distance matrices, an MMRR was performed using the R package PopGenReport (Adamack and Gruber, 2014), and the correlation coefficients and their significance were estimated based on 10,000 random permutations.

Connectivity Among Spirochaeta Populations

The amount and direction of gene flow among Spirochaeta populations were estimated using the coalescent-based program Lamarc v2.1.10 (Kuhner, 2006). A total of 10 runs was performed for each Spirochaeta OTU, consisting of likelihood searches of 20 initial and 2 final chains, with a minimum of 500 and 10,000 recorded trees, respectively, and sampling every 20 generations after a burn-in of 1,000 genealogies. The effective number of migrants per generation (Nm) among Spirochaeta populations was calculated by multiplying the maximum likelihood estimates (MLE) of the mutation parameter (θ) by the migration parameter (M), both estimated through the Lamarc program. We present the mean and standard deviation of the estimated Nm values obtained from the 10 runs for each Spirochaeta OTU.

Results

Sequencing Performance and OTU-Based Analysis

A total of 4,184,226 raw reads was generated from the 91 samples of external sediments and gut tissues. Once processed, 2,542,038 cleaned sequences distributed into 727 OTUs were obtained from the external sediment and gut tissue samples (details provided in Table 1). Out of this condensed dataset, 425,613 sequences associated with the Spirochaeta genus were retrieved, representing a total of 10 OTUs.

Both Bray–Curtis and Unweighted Unifrac distance methods did not reveal any difference in Spirochaeta OTU composition between the Patagonian sites (PAT1 and PAT2) (Figures 2A,B and Supplementary File 1). Conversely, Kerguelen Islands (KER) and martime Antarctic (KGI) sites were significantly different from each other in terms of Spirochaeta OTU composition and with the Patagonian ones (Supplementary File 1).

FIGURE 2

Figure 2. Non-metric multidimensional scaling scatter diagram of the Spirochaeta OTUs composition in gut tissue samples across the localities. Presence/absence matrix converted in Bray–Curtis distances (A) and Unweighted Unifrac distance (B). Colors are assigned to the locality.

The relative abundance analyses among the 10 Spirochaeta OTUs (Figure 3) showed four of them were shared among all the SO’s sampled provinces. Three OTUs (OTU6, OTU7, and OTU532) were more abundant in maritime Antarctica (KGI), four (OTU40, OTU42, OTU278 and OTU561) were more abundant in the Patagonian locality PAT1, two (OTU23 and OTU221) were predominantly found in the Patagonian locality PAT2, and a single one (OTU349) was predominant in Kerguelen Islands (Figure 3).

FIGURE 3

Figure 3. Clustering of Spirochaeta OTUs based on their relative abundances in each site. Clustering of Bray–Curtis distances matrix using the complete linkage method. The size of the circles indicates the repartition of a given OTU sequences among the four localities. The percentages indicate the OTUs’ proportions in the Spirochaeta dataset. Framed labels and black arrows indicate the selected OTUs that were selected to be process through the Minimum Decomposition Entropy pipeline (MED).

To test the genetic and phylogeographic structure of Spirochaeta at the broadest geographic scale available through our dataset, we selected the three most abundant Spirochaeta OTUs which were detected within the four localities of the dataset (i.e., OTU6, OTU7, and OTU40). These three selected co-distributed OTUs are good targets to constitute a metapopulation, which is a meaningful ecological unit of distinct local populations separated by gaps in habitats and interconnected to some extent via dispersal events of individuals (Halsey et al., 2015). The relative abundance in sample types, the closest sequence retrieved from Blast analysis, and the distribution of these OTUs among the localities are provided in the Supplementary Files 2, 3.

Microdiversity Within Spirochaeta OTUs

A total of 48, 96, and 48 oligotypes were defined for OTU6, OTU7, and OTU40, respectively (Supplementary File 4). Accumulation curves of OTU6, OTU7, and OTU40 oligotypes reached saturation in almost all localities indicating that the overall majority of Spirochaeta microdiversity has been found within the analyzed samples (Supplementary File 5). Diversity indices measured as N, k, S, h, and Π for each OTU in each locality are provided in Table 2. The genetic diversity (H) ranged from 0.0681 (OTU40 in KER) to 0.8036 (OTU7 in PAT1) across localities (Table 2). Patagonian sites exhibited higher oligotype and nucleotide diversity for OTU6 and OTU7 than maritime Antarctica and Kerguelen Islands localities. In contrast, the genetic diversity of the OTU40 oligotypes was higher in the maritime Antarctic site and lower for the Kerguelen population.

TABLE 2

Table 2. Summary of oligotypes number and genetic indices per OTU and per site for the three most abundant Spirochaeta OTUs found in all sampling localities.

Populations Differentiation and Phylogeographic Structure of Spirochaeta Oligotypes

Independently of the OTU considered, the genetic (F_st) and phylogeographic (Φ_st) structures between the two closest localities (i.e., Patagonian sites PAT1 and PAT2) were extremely to moderately low. In the case of the OTU40, the genetic diversity and frequencies of Spirochaeta oligotypes were fully homogenous between PAT1 and PAT2, as indicated by the non-significant values of F_st and Φ_st comparisons (Table 3). Contrarily, higher values of F_st and Φ_stT comparisons were recorded among the three provinces considered in this study (Patagonia, PAT1/PAT2; maritime Antarctica, KGI; Kerguelen Islands, KER) (Table 3). Consistently, the distribution of the Spirochaeta oligotypes was geographically discontinuous across the localities, with various province-specific oligotypes (Supplementary File 4). Two exceptions were observed, in the case of maritime Antarctica and Kerguelen Islands (KGI and KER) for OTU6, and in the case of Kerguelen Islands (KER) and Patagonia (PAT1) for OTU7, with relatively lower values of F_st and Φ_st (Table 3). The Snn test for phylogeographic structure among all sites was significant with statistic values ≥ 0.5 (OTU6; Snn = 0.50, p-value < 0.0001, OTU7; Snn = 0.57, p-value < 0.0001, OTU40; Snn = 0.50, p-value < 0.0001). All in all, these results indicate the existence of both genetically and geographically differentiated Spirochaeta populations across the three biogeographic provinces sampled.

TABLE 3

Table 3. Genetic (F_ST) and phylogeographic structure (Φ_ST) of the Spirochaeta populations among localities.

Within the 48 oligotypes identified in the OTU6, 11 (∼23%) were exclusive to one of the three provinces, and more than half were exclusive to Kerguelen Islands (Supplementary File 4). Maritime Antarctic and Kerguelen Islands (KGI and KER) shared 27 (∼66%) of their oligotypes. The Patagonian localities (PAT1/PAT2) shared 18 of the 23 total oligotypes (∼78%) observed in this province (Figure 4 and Supplementary File 2). A total of five (∼10%) oligotypes were broadly distributed across all localities. Oligotype network of OTU6 showed short genealogies and the presence of at least five dominant oligotypes. The dominant oligotype in Patagonia (PAT1/PAT2), and the dominant oligotype in maritime Antarctica and Kerguelen Islands (KGI/KER), were separated by a single substitution (Figure 4).

FIGURE 4

Figure 4. Median-joining oligotype networks of the three selected OTUs of Spirochaeta. Each circle represents a distinct oligotype. Colors indicate the locality of provenance. Circle size is scaled on the sum of the oligotype frequencies in each locality, within the corresponding OTU dataset. Raw abundances are provided in Supplementary File 2.

In the case of OTU7, the percentage of exclusive oligotypes was almost the same as the OTU6, with 21 (∼22%) oligotypes exclusive to one of the provinces (Supplementary File 4). The dominant oligotype was different between maritime Antarctic and Kerguelen Islands localities (KGI and KER) (Figure 4). While most oligotypes from Kerguelen Islands (KER) were detected in at least one of the Patagonian sites (PAT2 or PAT1) (∼94%), fewer oligotypes from maritime Antarctica (KGI) were observed in Kerguelen Islands (∼68%) and even fewer in the Patagonian localities (∼62%) (Figure 4).

For OTU40, we recorded a predominant group of oligotypes specific to the Patagonian sites representing 65% of the oligotypes identified within the OTU40 (Supplementary File 4 and Figure 4). A clear separation was observed in the oligotype network between the KGI/KER and PAT1/PAT2 localities (Figure 4), with only four shared oligotypes (∼8%) (Supplementary File 2). In maritime Antarctica (KGI), three of the four oligotypes were exclusive, whereas the dominant one from Kerguelen Islands (KER) was shared with maritime Antarctica (KGI) (∼8%) (Figure 4 and Supplementary File 4).

Gene Flow Under a Migration–Drift Equilibrium Model

All the analyzed OTUs showed high genetic similarities between the analyzed Patagonian populations (PAT1 and PAT2). Gene flow analyses identified a bidirectional pattern from PAT2 to PAT1 (effective number of migrants per generation, Nm > 4) and from PAT1 to PAT2 (Nm > 14) (Table 4). The connectivity between the Patagonian and maritime Antarctic localities showed relatively low values of Nm, ranging from 0.001 (OTU40, from PAT2 to KGI) to 0.64 (OTU7, from KGI to PAT1) (Table 4). The OTU6 and OTU7 were both characterized by a substantial gene flow between maritime Antarctica and Kerguelen Islands that was stronger in the direction KGI to KER (OTU6, Nm = 9.81 and OTU7, Nm = 1.39) than in the direction KER to KGI (OTU6, Nm = 2.63 and OTU7, Nm = 0.41) (Table 4). Contrarily, an unidirectional and low gene flow from KER to KGI (Nm = 0.87) was recorded for the OTU40 (Table 4). Finally, the connectivity between Patagonian (PAT1 and PAT2) and Kerguelen Islands (KER) localities was illustrated by three distinct patterns; a low-intensity flow (Nm < 0.5) predominant in the direction PAT1/PAT2 to KER for the OTU6, a substantial flow (Nm values from 2.00 up to 24.41) predominant in the direction KER to PAT1/PAT2 for the OTU7, and an absence of connectivity (Nm < 0.03) in the case of the OTU40 (Table 4). The gene flows are summarized in Figure 5.

TABLE 4

Table 4. Effective numbers of migrants per generation (Nm) among Spirochaeta populations of the OTU6, OTU7, and OTU40.

FIGURE 5

Figure 5. Gene flows summary and proposed dispersal routes across the Southern Ocean for each Spirochaeta taxa. Only the gene flows with Nm values > 0.25 are graphically represented. Discontinuous arrows represent Nm values > 0.25 and <1, continuous arrow represent Nm values > 1. Continuous arrows’ width is proportional to Nm values.

Contribution of Contemporary Selection Versus Historical Processes in Shaping the Spirochaeta Microdiversity

For each of the three selected Spirochaeta OTUs, neutral ecological processes were essential in shaping the population composition turnover in Abatus gut membrane. According to the quantitative parsing of ecological processes, the composition of Spirochaeta population was mostly driven by ecological drift, ranging from 50% (OTU40) to 74% (OTU6) of turnover, followed by dispersal limitation, ranging from 12% (OTU6) to 20% (OTU40) of turnover, and homogenizing selection, ranging from 9% (OTU6) to 19% (OTU40) of turnover (Table 5). Overall, deterministic processes (i.e., homogeneous and variable selection) did not account for more than 10% of the populations’ turnover.

TABLE 5

Table 5. Quantitative parsing of ecological processes driving populations turnover within Spirochaeta OTUs.

The MMRR approach was used to disentangle the relative effect of geographic distance environmental abiotic differences on the Spirochaeta oligotype compositions between samples. The geographic distance matrix was linearly correlated to the abundance-based similarity matrix of Spirochaeta population composition for OTU7 and OTU40, explaining about 31 and 67% of the observed variation, respectively (Table 6). In contrast, the geographic distance did not significantly impact Spirochaeta oligotype composition for OTU6 (Table 6). Whatever the OTU considered, the environment distance had a significant but slight contribution (<4% of the observed variation) to the Spirochaeta population composition (Table 6), and the global R² of the model (environmental and geographic distance) did not exceed 24%.

TABLE 6

Table 6. Multiple matrix regression with randomization (MMRR) to quantify the relative effects of isolation by distance (IBD) and isolation by environment (IBE) on oligotypes assemblage within Spirochaeta OTUs.

Discussion

In this study, we coupled 16S rRNA metabarcoding and oligotyping algorithm to reveal the microdiversity within three bacterial OTUs affiliated to the Spirochaeta genus, and co-distributed across Patagonia, Kerguelen Islands, and maritime Antarctic provinces of the SO. Through this innovative approach, we identified numerous oligotypes within each of the Spirochaeta OTUs. These oligotypes, corresponding to Spirochaeta sub-taxa, were characterized by contrasting geographic distribution and high levels of 16S rRNA gene similarity (>97%). Taking advantage of the populational taxonomic resolution provided by the oligotype definition, we depicted the Spirochaeta biogeographic patterns across the analyzed provinces in the SO, using various tools adapted from population genetics classically applied in phylogeographic study of macroorganisms’ models. Despite its low substitution rate [approximately 1% in 50 million years (Espejo and Plaza, 2018)], our study demonstrates that the 16S rRNA gene shows value in the evaluation of the microdiversification, as it offers the best compromise between an informative genetic signal, and robust screening of global microbial diversity at intra-OTU level, in a wide range of barely unknown habitats (Chase and Martiny, 2018).

Unlike the studies with macroorganisms, which are usually more demanding in terms of individual sampling effort, we benefit here from the high sequencing depth provided by the metabarcoding of a low-diversity habitat (i.e., the Abatus gut tissue). This allows a robust coverage of the Spirochaeta diversity (up to 180,000 sequences per OTU), and thus high precision of the oligotypes frequencies. Our methodology echoes the “metaphylogeographic” approach recently proposed by Turon et al. (2020) to investigate eukaryotic intraspecies diversity through COI (cytochrome c oxidase subunit I) gene amplicon-sequencing and an oligotyping-like cleaning protocol of the reads based on entropy variation. We propose to expand the concept of “metaphylogeography” to the prokaryotes since it permits phylogeographic inferences of uncultured microbes from a wide range of habitats.

The β-diversity analysis performed at the Spirochaeta genus level revealed that each of the three geographic provinces might host specific Spirochaeta OTUs representing distinct phylogenetic lineages. We also reported a non-random distribution trend with contrasting patterns of Spirochaeta OTU compositions across the localities. Nevertheless, about half of the Spirochaeta OTUs exhibited a broad distribution encompassing Patagonia, maritime Antarctica, and the Kerguelen Islands located more than 7,000 km to the east. This result suggests that despite being mostly detected in Abatus gut, and to a lesser extent in marine benthic sediments, some Spirochaeta representatives would disperse through the SO currents. Concordantly, previous campaigns of high-throughput sequencing of the ocean water column have consistently reported the presence of free-living Spirochaeta OTUs in surface to mesopelagic water, away from the coastlines (Pesant et al., 2015).

For each of the three assessed OTUs, the Spirochaeta populations were expected to be remarkably homogeneous between the two Patagonian sites due to the geographic vicinity and the absence of an evident oceanographic barrier. Consistently with this assumption, and for each of the three OTUs, most of the Spirochaeta oligotypes were shared, the lowest genetic and phylogeographic structures were reported, and high levels of gene flow were recorded between these two sites. Similarly, low or absent differentiation patterns along the Atlantic coast of Patagonia were previously reported for marine Patagonian macroorganisms, including notothenioid fishes (Ceballos et al., 2016), scorched mussels (Trovant et al., 2015), and pulmonate gastropods (Iriarte et al., 2020) presumably due to their high dispersal potential and the ecological continuum of the sampled localities that may conform a same biogeographic province connected through the equator-ward Falkland current (Phillpot, 1985; Arkhipkin et al., 2004). Further phylogeographic studies focusing on microbial taxa of additional sampling sites from Atlantic Patagonia should confirm the microbial biogeographic consistency of this province.

Between Patagonian and maritime Antarctic provinces, Spirochaeta populations exhibited strong genetic and phylogeographic structures, and low levels of gene flow were estimated between these two provinces. These results corroborate our hypothesis that the APF hinders individual dispersion and genetic homogeneity among bacterial populations and suggest that the geographically structured Spirochaeta populations from these two provinces are genetically diverging over time (Nagylaki, 1980; Charlesworth, 2003). Previous studies focusing on diverse macroorganisms taxonomic groups of the SO have evidenced the critical role of the APF on biogeographic patterns, as an open-ocean barrier inducing a genetic break between South America and Antarctica (e.g., ribbon worms (Thornhill et al., 2008); brittle stars (Hunter and Halanych, 2008); notothenioid fishes (Hüne et al., 2015); limpets (González-Wevar et al., 2012, 2019; González−Wevar et al., 2017); sea urchins (Díaz et al., 2011). Regarding the microbial distribution patterns, significant β-diversity differences between prokaryotes assembly from both sides of the APF have been reported in the past, but most of the studies focused on global community in the water column, at high taxonomic resolution [summarized in Flaviani et al. (2018)]. Here, we extend this discontinuity in bacterial diversity to a fine taxonomic resolution (i.e., intra-OTU), revealing province-restricted oligotypes and strong genetic and phylogeographic structure between Patagonian and maritime Antarctic Spirochaeta populations.

Contrarily, and despite the substantial geographic distance separating the sub-Antarctic Kerguelen Islands and the Patagonian and Antarctic sites (>6,500 km), population genetic analyses suggest the existence of some level of connectivity between Kerguelen and the other sites. These findings support a potential dispersion of Spirochaeta taxa from Patagonia and maritime Antarctica to the Kerguelen Islands, and contrariwise, from the Kerguelen Islands to Patagonia and maritime Antarctica. As evidenced by the numerous shared oligotypes, such connectivity would maintain a sufficient gene flow among these provinces to partially counteract the genetic divergence driven by selection, mutation, and genetic drift, inducing oligotypes mixing, and limiting the spatial differentiation of oligotypes assembly (Martiny et al., 2006; Orsini et al., 2013). Moreover, we suggest that this gene flow is not bidirectional, but governed by exclusively eastward oriented dispersion routes (Figure 5), following the major and constant flow of the ACC (Güller et al., 2020). Under this scenario, Spirochaeta individuals from Kerguelen Islands may seed toward Patagonia following the ACC eastward flow around Antarctica. Such ACC-mediated connectivity among sub-Antarctic provinces (Patagonia and Kerguelen Islands) is well known in a wide range of benthic macroorganisms populations, such as buoyant kelps Durvillaea antarctica and Macrocystis pyrifera (Macaya and Zuccarello, 2010), and several kelp-associated macroinvertebrates (Leese et al., 2010; Nikula et al., 2010; Cumming et al., 2014; González−Wevar et al., 2018; Güller et al., 2020). Occasionally, Spirochaeta individuals from Kerguelen Islands may also be able to reach the maritime Antarctic province. Such pattern has been recently reported for the southern bull kelp D. antarctica, a typical sub-Antarctic macroalgae, which is transported by rafting to as far as the West Antarctic Peninsula coasts, pushed by the circumpolar flow of the ACC or by storms leading to the occasional crossing of the APF (Fraser et al., 2020).

Several studies have provided evidence of a high dispersal capacity of marine bacteria by comparing community composition mostly at high taxonomic resolution (e.g., class, genus, or OTU) among various water masses and oceanic regions (Sul et al., 2013; Sunagawa et al., 2015; Milici et al., 2017; Logares et al., 2020). Particularly, the most abundant marine bacteria are supposed to migrate between adjacent regions through passive transport (Sul et al., 2013). An innovative conceptual framework called “Microbial Conveyor Belt” (MCB) has been proposed by Mestre and Höfer (2020), to emphasize that the marine microorganisms’ dispersion would not merely rely on passive and stochastic dispersal, but instead on the adaptation of life-history traits (e.g., dormancy stage). These traits would allow microorganisms to successfully and recurrently disperse in unfavorable habitats through specific dispersion avenues (Locey, 2010). Here, we provided empiric results from Spirochaeta population based on genetic data supporting a partial MCB in the SO driven by the ACC. Unfortunately, details about the benthic Spirochaeta taxa’s ecology are scarce, with a single isolated strain from subseafloor sediment (Miyazaki et al., 2014). Thus, the life-history traits of Spirochaeta, as the sporulation capacity, remain to be investigated to further understand its distribution pattern in the SO. Nevertheless, in order to disperse, we propose that Spirochaeta individuals (enriched in the digestive tract) could be released from the host gut toward the surrounding benthic sediments through fecal pellets. Such enrichment of the digesta with taxa from the host microbiota, as well as the presence of Spirochaeta within the fecal material, have been demonstrated in the sea urchin species Lytechinus variegatus (Hakim, 2015). The released Spirochaeta individuals may be resuspended in the water column through the action of one or several processes such as upwelling, bioturbation by the benthic deposit-feeders, or water column mixing during winter (Meysman et al., 2006; Petro et al., 2017; Mestre and Höfer, 2020). Once in the water column, these Spirochaeta individuals may disperse over large geographic scales, transported through oceanographic features (e.g., currents, punctual meteorological events) (Mestre and Höfer, 2020). The attachment to suspended particulate matter, either biotic [e.g., hitchhiking on zooplankton (Grossart et al., 2010) and seaweed (Serebryakova et al., 2018)] or abiotic [e.g., microplastics (Bowley et al., 2020), known to have a long-distance dispersion potential], may also contribute to the bacterial spreading in the oceans (Milici et al., 2017; Mestre and Höfer, 2020).

The marine prokaryote communities are usually considered widely dispersed and mainly shaped by contemporary ecological processes such as environmental filtering (Sunagawa et al., 2015; Louca et al., 2016). By applying the ecological framework developed by Stegen et al. (2013) to oligotypes data, we found contrarily that ecological drift was the predominant stochastic mechanism shaping intra-populations turnover within Spirochaeta taxa across the SO. Our previous study of the Abatus gut microbiota showed that non-neutral processes drove the bacterial community at the OTU level in the host gut tissue (Schwob et al., 2020). While deterministic processes are usually prevalent in structuring microbial communities’ assembly at a higher taxonomic resolution (Martiny et al., 2011; Wang, 2013; Sintes et al., 2015; Larkin and Martiny, 2017), the stochastic mechanisms tend to have a more significant contribution at finer taxonomic scales (Logares et al., 2020), since niche overlapping and functional redundancy enhance the susceptibility of populations to drift (Zhou and Ning, 2017). Thus, the biogeographic structure observed among Spirochaeta populations might result from stochastic birth, death, disturbance, emigration, and immigration events rather than oligotype-sorting through the biotic and abiotic environmental variations (Martiny et al., 2006; Nemergut et al., 2013; Vellend, 2016). Consistently, the MMRR analysis revealed that the isolation-by-environment (IBE) model might account for a low percentage of the Spirochaeta oligotypes turnover. Altogether, these results tend to validate the strategy applied in our study, that is, to focus on specialist bacterial taxa hosted in sibling sea urchin species with the same habitat preferences, in order to homogenize the environment, to reduce the diversity, to soften the deterministic selection driven by environmental variations, thus leading to the maximization of the detection of neutral micro-evolutive processes associated with biogeography (Langenheder and Lindstrom, 2019).

By analogy with the genetic drift, whereby changes in gene frequencies occur solely by chance in a population (Chase et al., 2011), our result suggests that the microdiversity observed within the Spirochaeta taxa would be mostly generated by genetic drift without any adaptive implications. An earlier study reported that microdiversity observed in the 16S rRNA gene of marine coastal Vibrio splendidus isolates was ecologically neutral (Thompson et al., 2005). Nevertheless, we cannot discard that, while the microdiversity within the V4–V5 of 16S rRNA gene-targeted here is likely to be acquired through neutral processes (Chust et al., 2016), it may also be associated with substantial modifications in niche-defining traits and functional attributes specific of the Spirochaeta strains, driven by deterministic processes, in order to cope with local conditions (Zure et al., 2017; Chase and Martiny, 2018). Further studies will need to focus on other loci (e.g., functional genes), potentially under selection, as they are expected to display a higher degree of differentiation among populations and to provide an insight into the ecology of the Spirochaeta sub-taxa (White et al., 2010).

Given the highly specific targeted habitat (i.e., gut tissue) and the intricate association between Spirochaeta and Abatus, it remains unclear whether or not the host contributes to the Spirochaeta genetic and phylogeographic pattern reported in the gut populations. Indeed, the Abatus guts may have individual colonization histories also contributing to the microbial community assembly process at the population level. However, due to the limited sequencing coverage of Spirochaeta diversity outside the host, we cannot evaluate to what extent the most abundant oligotypes from the gut tissue, also detected in the external sediment samples, would be representative of the composition of free-living populations, thus preventing a comprehensive understanding of the host contribution to Spirochaeta oligotype assembly. Ultimately, conclusion about the predominant effect of stochastic ecological processes onto Spirochaeta oligotype composition need to be nuanced, as potential historical contingencies related to the host can eventually shape the final state of microbiota composition through both stochastic and/or deterministic processes (Obadia et al., 2017). In strong instance, host early-life microbial colonization by early colonizers benefiting from a “priority effect” may influence the order and timing of the historical sequence of species arrival within the local community (Fukami, 2015). Under the hypothesis of a predominantly stochastic gut colonization by random microbes, diverging chronologies of the ensuing microbial species would be expected, leading toward different stable states and a noisy community assembly across Abatus populations (Obadia et al., 2017). Contrastingly, we reported patterns of high dominance and prevalence of Spirochaeta across Abatus individuals, and high homogeneity of oligotypes composition within populations (data not shown), rather suggesting that Spirochaeta may benefit from a priority effect during Abatus early-life, and that the subsequent interactions may deterministically constrain gut colonization success (Fukami, 2015; Higgins et al., 2020). Moreover, while there is no evidence of vertical transmission in Abatus, the transovarial transmission of a bacterial endosymbiont was reported in the regular sea urchin species Heliocidaris erythrogramma (Carrier et al., 2021) and maternal-inherited bacterial communities were detected in unfertilized eggs of Strongylocentrotus purpuratus (Carrier and Reitzel, 2019). The presence of preexisting Spirochaeta strains within the gut may also limit the ensuing horizontal colonization and constrain the microbiota composition toward a homogenous and stabilized state (Renelies-Hamilton et al., 2021). The presence of vertically transmitted Spirochaeta, the compositional transition states of Abatus gut microbiota, and the biogeographic structure of free-living Spirochaeta populations in the external sediment remain to be further investigated.

Notwithstanding the consistency of the global phylogeographic and connectivity patterns depicted across the three tested OTUs, we also reported some differences according to the taxa considered, which might be related to different ecotypes with distinct ecological niches or different dispersal capacity. For instance, various marine bacterial taxa, such as the cyanobacteria Synechococcus or the Vibrio populations, demonstrate fine-tuning of their physiology by accumulating microdiversity in functional genes through duplication events, SNPs, and allelic variants (Everroad and Wood, 2012; Shapiro et al., 2012). Alternatively, these differences may also be related to the Spirochaeta OTU abundance, since the more relatively abundant the Spirochaeta populations were (i.e., higher number of sequences retrieved from the gut tissue through the metabarcoding approach), the more they tend to exhibit cosmopolitan oligotypes (i.e., detected across each of the four localities). It is not unreasonable to infer that a larger population may have more chance to migrate and successfully reach a suitable habitat, while small-size populations may be more likely diluted along the dispersal route with no/too few dispersive particles to establish in the new habitat (Sul et al., 2013).

The diversity units defined by 16S rRNA gene sequences are generally considered as insensitive to diversification resulting from dispersal limitation (Hanson et al., 2012). Contrastingly, we reported that the dispersal limitation was the second most crucial ecological factor driving the turnover of Spirochaeta oligotypes, and by extension, their genetic divergence. Dispersal limitation is classically considered as a historical factor since current oligotypes assemblage results from past dispersal limitations (Martiny et al., 2006). Our result indicates that the potentially suitable habitats are too distant (Heino et al., 2015), or inaccessible due to the existence of oceanic currents (Sul et al., 2013; Logares et al., 2018), hence limiting the homogenization of Spirochaeta oligotypes’ frequencies across populations and allowing the neutral genetic divergence of genomic regions overtime via genetic drift (Nosil et al., 2008; Hanson et al., 2012). Note that our results obtained from distinct methodologies [i.e., the genetic differentiation and phylogeographic structure, the contribution of dispersal limitation from Stegen et al. framework (Stegen et al., 2013), and the contribution of the geographic distance isolation by distance (IBD) from the MMRR analysis] were highly consistent with each other, and across the three selected Spirochaeta OTUs. For instance, the OTU40 that harbored the overall highest value of genetic divergence was also characterized by the highest estimated contribution of geographic distance and dispersal limitation, thus supporting the interrelation between genetic divergence and oligotypes population turnover, and the overall consistency of the approach implemented.

Conclusion

Our study highlights the application of V4–V5 16S rRNA gene metabarcoding and oligotyping approach as rapid, robust, and resolutive enough to unravel marine bacterial phylogeographic patterns and detect genetic connectivity among the SO provinces. Taken together, the three Spirochaeta OTUs analyzed evidence three consistent phylogeographic patterns, classically observed in the studies involving benthic macroinvertebrates across the SO: (1) a high populational and genetic homogeneity within the Patagonia province, (2) a strong barrier to dispersal between Patagonia and maritime Antarctica due to the APF, resulting in a high differentiation of Spirochaeta populations, and (3) the existence of connectivity between sub-Antarctic provinces of the Kerguelen Islands and Patagonia, and from Kerguelen Islands to the maritime Antarctic, due to the ACC-mediated connectivity. Nevertheless, as connected as these provinces are, the gene flow does not seem to be strong enough to prevent the ongoing intraspecific microdiversification of the Spirochaeta taxa. The microdiversity of Spirochaeta, underlying these biogeographic patterns, is essentially driven by historical processes, such as ecological and genetic drift, and dispersal limitation related to the SO’s oceanographic features. In the future, extending this framework to other localities and taxonomic groups will contribute to the comprehensive understanding of the SO microbiota.

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 below: https://www.ncbi.nlm.nih.gov/ PRJNA658980, https://www.ncbi.nlm.nih.gov/ PRJNA590493, and https://www.ncbi.nlm.nih.gov/ PRJNA659050.

Author Contributions

EP, JO, and LC designed the study. EP and JO organized the sampling missions. EP, JO, LC, and GS collected samples. GS extracted the DNA, organized sequencing and managed data mining and analyses. NS contributed to the gene-flow and the MMRR analyses and designed the illustrative maps. GS, JO, NS, CG-W, and EP interpreted the results. GS wrote the manuscript. All authors contributed substantially to manuscript revisions and read and approved the final manuscript.

Funding

This work was financially financed by the project ANID/CONICYT PIA ACT 172065. Additionally, this research was supported by the post-doctoral projects ANID/CONICYT FONDECYT 3200036 (GS) and 3190482 (NS), and the regular project ANID/CONICYT FONDECYT 1211672 (LC). Field access facilities to Kerguelen Islands were supported by the French Polar Institute-IPEV (program No. 1044 PROTEKER).

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 Peter Beerli and Franz-Sebastian Krah for their technical supports. We also thank Jonathan Flores for providing Abatus samples from Argentinian Patagonia, Karin Gérard for managing the sampling logistic in Chilean Patagonia, and Ignacio Garrido, Paulina Brüning, and Miguel Pardo of the IDEAL Centre for their help in Abatus sampling in King George Island. We recognize the Chilean Antarctic Institute (INACH) for the logistic support during the Chilean Antarctic Expeditions (ECA 55 and 56).

Supplementary Material

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

Supplementary File 1 | Pairwise PERMANOVA on Spirochaeta OTUs composition dissimilarities among localities. p-values are adjusted using the default Bonferroni method implemented in the pairwiseAdonis R package and are considered as significant < 0.05.

Supplementary File 2 | Abundance and closest sequence retrieved from Blast analysis for each of the three OTUs analyzed through the MED pipeline. (1) Bowman JP, McCuaig RD: Biodiversity, community structural shifts, and biogeography of prokaryotes within Antarctic continental shelf sediment. Appl Environ Microbiol 2003, 69:2463–2483. (2) Acosta-González A, Rosselló-Móra R, Marqués S: Characterization of the anaerobic microbial community in oil-polluted subtidal sediments: aromatic biodegradation potential after the Prestige oil spill. Environmental microbiology 2013, 15:77–92.

Supplementary File 3 | Relative contribution of each locality in the total abundance of OTU6, OTU7, and OTU40 sequences. Colors are assigned to the different localities.

Supplementary File 4 | Summary of number of oligotypes in each locality and per OTU of Spirochaeta.

Supplementary File 5 | Accumulation curves of OTU6 (A), OTU7 (B), and OTU40 (C) oligotypes richness. Colors are assigned to each locality. Extrapolation is calculated from Hill numbers of richness (q = 0).

Footnotes

1. ^ http://oba.mio.osupytheas.fr/ocean-atlas/

References

Adamack, A. T., and Gruber, B. (2014). PopGenReport: simplifying basic population genetic analyses in R. Methods Ecol. Evol. 5, 384–387.

Google Scholar

Almasia, R., Carú, M., Handford, M., and Orlando, J. (2016). Environmental conditions shape soil bacterial community structure in a fragmented landscape. Soil Biol. Biochem. 103, 39–45. doi: 10.1016/j.soilbio.2016.08.004

CrossRef Full Text | Google Scholar

Arkhipkin, A. I., Middleton, D., Sirota, A., and Grzebielec, R. (2004). The effect of Falkland current inflows on offshore ontogenetic migrations of the squid Loligo gahi on the southern shelf of the Falkland Islands. Estuar. Coast. Shelf Sci. 60, 11–22. doi: 10.1016/j.ecss.2003.11.016

CrossRef Full Text | Google Scholar

Assis, J., Tyberghein, L., Bosch, S., Verbruggen, H., Serrão, E. A., and De Clerck, O. (2018). Bio−ORACLE v2. 0: extending marine data layers for bioclimatic modelling. Glob. Ecol. Biogeogr. 27, 277–284. doi: 10.1111/geb.12693

CrossRef Full Text | Google Scholar

Astorga, A., Oksanen, J., Luoto, M., Soininen, J., Virtanen, R., and Muotka, T. (2012). Distance decay of similarity in freshwater communities: do macro−and microorganisms follow the same rules? Glob. Ecol. Biogeogr. 21, 365–375. doi: 10.1111/j.1466-8238.2011.00681.x

CrossRef Full Text | Google Scholar

Bay, S. K., Mcgeoch, M. A., Gillor, O., Wieler, N., Palmer, D. J., Baker, D. J., et al. (2020). Soil bacterial communities exhibit strong biogeographic patterns at fine taxonomic resolution. mSystems 5:e00540-20.

Google Scholar

Beiruti, Z. A., Bhatnagar, S., Tegetmeyer, H. E., Geelhoed, J. S., Strous, M., and Ruff, S. E. (2017). Transient exposure to oxygen or nitrate reveals ecophysiology of fermentative and sulfate-reducing benthic microbial populations. bioRxiv [Preprint]. doi: 10.1101/146886

CrossRef Full Text | Google Scholar

Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10, 57–59. doi: 10.1038/nmeth.2276

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowley, J., Baker-Austin, C., Porter, A., Hartnell, R., and Lewis, C. (2020). Oceanic hitchhikers–assessing pathogen risks from marine microplastic. Trends Microbiol. 29, 107–116. doi: 10.1016/j.tim.2020.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowman, J. P., and McCuaig, R. D. (2003). Biodiversity, community structural shifts, and biogeography of prokaryotes within Antarctic continental shelf sediment. Appl. Environ. Microbiol. 69, 2463–2483. doi: 10.1128/aem.69.5.2463-2483.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Buttigieg, P. L., and Ramette, A. (2014). A guide to statistical analysis in microbial ecology: a community-focused, living review of multivariate data analyses. FEMS Microbiol Ecol 90, 543–550. doi: 10.1111/1574-6941.12437

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrier, T. J., Leigh, B. A., Deaker, D. J., Devens, H. R., Wray, G. A., Bordenstein, S. R., et al. (2021). Microbiome reduction and endosymbiont gain from a switch in sea urchin life history. Proc. Natl. Acad. Sci. U.S.A. 118:e2022023118. doi: 10.1073/pnas.2022023118

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrier, T. J., and Reitzel, A. M. (2019). Bacterial community dynamics during embryonic and larval development of three confamilial echinoids. Mar. Ecol. Prog. Ser. 611, 179–188. doi: 10.3354/meps12872

CrossRef Full Text | Google Scholar

Ceballos, S. G., Lessa, E. P., Licandeo, R., and Fernandez, D. A. (2016). Genetic relationships between Atlantic and Pacific populations of the notothenioid fish Eleginops maclovinus: the footprints of quaternary glaciations in Patagonia. Heredity (Edinb.) 116, 372–377. doi: 10.1038/hdy.2015.106

PubMed Abstract | CrossRef Full Text | Google Scholar

Charlesworth, D. (2003). Effects of inbreeding on the genetic diversity of populations. Philos. Trans. R. Soc. Lond. SerB Biol. Sci. 358, 1051–1070.

Google Scholar

Chase, A. B., and Martiny, J. B. (2018). The importance of resolving biogeographic patterns of microbial microdiversity. Microbiol. Aust. 39, 5–8. doi: 10.1071/ma18003

CrossRef Full Text | Google Scholar

Chase, J. M., Kraft, N. J., Smith, K. G., Vellend, M., and Inouye, B. D. (2011). Using null models to disentangle variation in community dissimilarity from variation in α−diversity. Ecosphere 2, 1–11.

Google Scholar

Chessel, D., Dufour, A.-B., and Dray, S. (2009). Package ‘ade4’: Analysis of Ecological Data: Exploratory and Euclidean Methods in Environmental Sciences. R-package Version: 1.5-1.

Google Scholar

Chust, G., Villarino, E., Chenuil, A., Irigoien, X., Bizsel, N., Bode, A., et al. (2016). Dispersal similarly shapes both population genetics and community patterns in the marine realm. Sci. Rep. 6:28730.

Google Scholar

Cumming, R., Nikula, R., Spencer, H., and Waters, J. (2014). Transoceanic genetic similarities of kelp−associated sea slug populations: long−distance dispersal via rafting? J. Biogeogr. 41, 2357–2370. doi: 10.1111/jbi.12376

CrossRef Full Text | Google Scholar

David, B., Choné, T., Mooi, R., and De Ridder, C. (2005). “Antarctic echinoidea,” in Synopses of the Antarctic benthos, eds J. W. Wägele and J. Sieg (Ruggell: ARG Gantner).

Google Scholar

De Broyer, C., Koubbi, P., Griffiths, H., Raymond, B., d’Acoz, U., Van De Putte, A., et al. (2014). Biogeographic Atlas of the Southern Ocean. Cambridge: Scientific Committee on Antarctic Research.

Google Scholar

Díaz, A., Féral, J.-P., David, B., Saucède, T., and Poulin, E. (2011). Evolutionary pathways among shallow and deep-sea echinoids of the genus Sterechinus in the Southern Ocean. Deep Sea Res. 2 Top. Stud. Oceanogr. 58, 205–211. doi: 10.1016/j.dsr2.2010.10.012

CrossRef Full Text | Google Scholar

Díaz, A., González-Wevar, C. A., Maturana, C. S., Palma, A. T., Poulin, E., and Gerard, K. (2012). Restricted geographic distribution and low genetic diversity of the brooding sea urchin Abatus agassizii (Spatangoidea: Schizasteridae) in the South Shetland Islands: a bridgehead population before the spread to the northern Antarctic Peninsula? Rev. Chil. Hist. Nat. 85, 457–468. doi: 10.4067/s0716-078x2012000400008

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., and Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics 31, 3476–3482. doi: 10.1093/bioinformatics/btv401

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381

PubMed Abstract | CrossRef Full Text | Google Scholar

Eren, A. M., Maignien, L., Sul, W. J., Murphy, L. G., Grim, S. L., Morrison, H. G., et al. (2013). Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol. Evol. 4, 1111–1119. doi: 10.1111/2041-210x.12114

PubMed Abstract | CrossRef Full Text | Google Scholar

Eren, A. M., Morrison, H. G., Lescault, P. J., Reveillaud, J., Vineis, J. H., and Sogin, M. L. (2015). Minimum entropy decomposition: unsupervised oligotyping for sensitive partitioning of high-throughput marker gene sequences. ISME J. 9, 968–979. doi: 10.1038/ismej.2014.195

PubMed Abstract | CrossRef Full Text | Google Scholar

Espejo, R. T., and Plaza, N. (2018). Multiple ribosomal RNA operons in bacteria; their concerted evolution and potential consequences on the rate of evolution of their 16S rRNA. Front. Microbiol. 9:1232. doi: 10.3389/fmicb.2018.01232

PubMed Abstract | CrossRef Full Text | Google Scholar

Everroad, R. C., and Wood, A. M. (2012). Phycoerythrin evolution and diversification of spectral phenotype in marine Synechococcus and related picocyanobacteria. Mol. Phylogenet. Evol. 64, 381–392. doi: 10.1016/j.ympev.2012.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenchel, T., and Finlay, B. J. (2004). The ubiquity of small species: patterns of local and global diversity. Bioscience 54, 777–784. doi: 10.1641/0006-3568(2004)054[0777:tuossp]2.0.co;2

CrossRef Full Text | Google Scholar

Feng, Y., Chen, R., Stegen, J. C., Guo, Z., Zhang, J., Li, Z., et al. (2018). Two key features influencing community assembly processes at regional scale: initial state and degree of change in environmental conditions. Mol. Ecol. 27, 5238–5251. doi: 10.1111/mec.14914

PubMed Abstract | CrossRef Full Text | Google Scholar

Flaviani, F., Schroeder, D. C., Lebret, K., Balestreri, C., Highfield, A. C., Schroeder, J. L., et al. (2018). Distinct oceanic microbiomes from viruses to protists located near the Antarctic Circumpolar Current. Front. Microbiol. 9:1474. doi: 10.3389/fmicb.2018.01474

PubMed Abstract | CrossRef Full Text | Google Scholar

Fraser, C. I., Morrison, A., and Rojas, P. O. (2020). “Biogeographic Processes Influencing Antarctic and sub-Antarctic Seaweeds,” in Antarctic Seaweeds, eds I. Gómez and P. Huovinen (Cham: Springer), 43–57. doi: 10.1007/978-3-030-39448-6_3

CrossRef Full Text | Google Scholar

Fraser, C. I., Nikula, R., Spencer, H. G., and Waters, J. M. (2009). Kelp genes reveal effects of subantarctic sea ice during the Last Glacial Maximum. Proc. Natl. Acad. Sci. U.S.A. 106, 3249–3253. doi: 10.1073/pnas.0810635106

PubMed Abstract | CrossRef Full Text | Google Scholar

Frugone, M. J., López, M. E., Segovia, N. I., Cole, T. L., Lowther, A., Pistorius, P., et al. (2019). More than the eye can see: genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean. Mol. Phylogenet. Evol. 139:106563. doi: 10.1016/j.ympev.2019.106563

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukami, T. (2015). Historical contingency in community assembly: integrating niches, species pools, and priority effects. Annu. Rev. Ecol. Evol. Syst. 46, 1–23. doi: 10.1146/annurev-ecolsys-110411-160340

CrossRef Full Text | Google Scholar

Gérard, K., Bierne, N., Borsa, P., Chenuil, A., and Féral, J.-P. (2008). Pleistocene separation of mitochondrial lineages of Mytilus spp. mussels from Northern and Southern Hemispheres and strong genetic differentiation among southern populations. Mol. Phylogenet. Evol. 49, 84–91. doi: 10.1016/j.ympev.2008.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilbert, J. A., Steele, J. A., Caporaso, J. G., Steinbruck, L., Reeder, J., Temperton, B., et al. (2012). Defining seasonal marine microbial community dynamics. ISME J. 6, 298–308. doi: 10.1038/ismej.2011.107

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Wevar, C. A., Hüne, M., Cañete, J. I., Mansilla, A., Nakano, T., and Poulin, E. (2012). Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet Nacella magellanica (Gmelin, 1791). BMC Evol. Biol. 12:139. doi: 10.1186/1471-2148-12-139

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Wevar, C. A., Hüne, M., Rosenfeld, S., Nakano, T., Saucède, T., Spencer, H., et al. (2019). Systematic revision of Nacella (Patellogastropoda: Nacellidae) based on a complete phylogeny of the genus, with the description of a new species from the southern tip of South America. Zool. J. Linnean Soc. 186, 303–336. doi: 10.1093/zoolinnean/zly067

CrossRef Full Text | Google Scholar

González−Wevar, C. A., Hüne, M., Segovia, I., Nakano, T., Spencer, H. G., Chown, S. L., et al. (2017). Following the Antarctic Circumpolar Current: patterns and processes in the biogeography of the limpet Nacella (Mollusca: Patellogastropoda) across the Southern Ocean. J. Biogeogr. 44, 861–874. doi: 10.1111/jbi.12908

CrossRef Full Text | Google Scholar

González−Wevar, C. A., Segovia, I., Rosenfeld, S., Ojeda, J., Hüne, M., Naretto, J., et al. (2018). Unexpected absence of island endemics: long−distance dispersal in higher latitude sub−Antarctic Siphonaria (Gastropoda: Euthyneura) species. J. Biogeogr. 45, 874–884. doi: 10.1111/jbi.13174

CrossRef Full Text | Google Scholar

Griffiths, H. J. (2010). Antarctic marine biodiversity–what do we know about the distribution of life in the Southern Ocean? PLoS One 5:e11683. doi: 10.1371/journal.pone.0011683

PubMed Abstract | CrossRef Full Text | Google Scholar

Grossart, H.-P., Dziallas, C., Leunert, F., and Tang, K. W. (2010). Bacteria dispersal by hitchhiking on zooplankton. Proc. Natl. Acad. Sci. 107, 11959–11964. doi: 10.1073/pnas.1000668107

PubMed Abstract | CrossRef Full Text | Google Scholar

Gruber, N., Landschützer, P., and Lovenduski, N. S. (2019). The variable Southern Ocean carbon sink. Annu. Rev. Mar. Sci. 11, 159–186. doi: 10.1146/annurev-marine-121916-063407

PubMed Abstract | CrossRef Full Text | Google Scholar

Guillaumot, C., Fabri−Ruiz, S., Martin, A., Eléaume, M., Danis, B., Féral, J. P., et al. (2018). Benthic species of the Kerguelen Plateau show contrasting distribution shifts in response to environmental changes. Ecol. Evol. 8, 6210–6225. doi: 10.1002/ece3.4091

PubMed Abstract | CrossRef Full Text | Google Scholar

Guindon, S., Dufayard, J.-F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010

PubMed Abstract | CrossRef Full Text | Google Scholar

Güller, M., Puccinelli, E., and Zelaya, D. G. (2020). The Antarctic Circumpolar Current as a dispersive agent in the Southern Ocean: evidence from bivalves. Mar. Biol. 167:143.

Google Scholar

Hakim, J. A. (2015). Comparison of Gut Microbiomes in Laboratory Cultured Sea Urchins Revealing Selective Attributes of Microbial Composition Based Upon Their Feed and Surroundings. Ph. D. thesis. Birmingham, AL: The University of Alabama at Birmingham.

Google Scholar

Halanych, K. M., and Mahon, A. R. (2018). Challenging dogma concerning biogeographic patterns of Antarctica and the Southern Ocean. Annu. Rev. Ecol. Evol. Syst. 49, 355–378. doi: 10.1146/annurev-ecolsys-121415-032139

CrossRef Full Text | Google Scholar

Halsey, S. J., Bell, T. J., Mceachern, K., and Pavlovic, N. B. (2015). Comparison of reintroduction and enhancement effects on metapopulation viability. Restor. Ecol. 23, 375–384. doi: 10.1111/rec.12191

CrossRef Full Text | Google Scholar

Hanson, C. A., Fuhrman, J. A., Horner-Devine, M. C., and Martiny, J. B. (2012). Beyond biogeographic patterns: processes shaping the microbial landscape. Nat. Rev. Microbiol. 10, 497–506. doi: 10.1038/nrmicro2795

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanson, C. A., Müller, A. L., Loy, A., Dona, C., Appel, R., Jørgensen, B. B., et al. (2019). Historical factors associated with past environments influence the biogeography of thermophilic endospores in Arctic marine sediments. Front. Microbiol. 10:245. doi: 10.3389/fmicb.2019.00245

PubMed Abstract | CrossRef Full Text | Google Scholar

Heino, J., Melo, A. S., Siqueira, T., Soininen, J., Valanko, S., and Bini, L. M. (2015). Metacommunity organisation, spatial extent and dispersal in aquatic systems: patterns, processes and prospects. Freshw. Biol. 60, 845–869. doi: 10.1111/fwb.12533

CrossRef Full Text | Google Scholar

Higgins, S. A., Panke−Buisse, K., and Buckley, D. H. (2020). The biogeography of Streptomyces in New Zealand enabled by high−throughput sequencing of genus−specific rpoB amplicons. Environ. Microbiol. 23, 1452–1468. doi: 10.1111/1462-2920.15350

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsieh, T., Ma, K., and Chao, A. (2016). iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers). Methods Ecol. Evol. 7, 1451–1456. doi: 10.1111/2041-210x.12613

CrossRef Full Text | Google Scholar

Hudson, R. R. (2000). A new statistic for detecting genetic differentiation. Genetics 155, 2011–2014. doi: 10.1093/genetics/155.4.2011

CrossRef Full Text | Google Scholar

Hüne, M., González-Wevar, C., Poulin, E., Mansilla, A., Fernández, D. A., and Barrera-Oro, E. (2015). Low level of genetic divergence between Harpagifer fish species (Perciformes: Notothenioidei) suggests a quaternary colonization of Patagonia from the Antarctic Peninsula. Polar Biol. 38, 607–617. doi: 10.1007/s00300-014-1623-6

CrossRef Full Text | Google Scholar

Hunter, R. L., and Halanych, K. M. (2008). Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the Drake Passage in the Southern Ocean. J. Hered. 99, 137–148. doi: 10.1093/jhered/esm119

PubMed Abstract | CrossRef Full Text | Google Scholar

Iriarte, P. F., González-Wevar, C., Segovia, N., Rosenfeld, S., Hüne, M., Fainburg, L., et al. (2020). Quaternary ice sheets and sea level regression drove divergence in a marine gastropod along Eastern and Western coasts of South America. Sci. Rep. 10, 1–15.

Google Scholar

Karimi, B., Terrat, S., Dequiedt, S., Saby, N. P., Horrigue, W., Lelièvre, M., et al. (2018). Biogeography of soil bacteria and archaea across France. Sci. Adv. 4:eaat1808. doi: 10.1126/sciadv.aat1808

PubMed Abstract | CrossRef Full Text | Google Scholar

Koubbi, P., De Broyer, C., Griffiths, H., Raymond, B., D’udekem D’acoz, C., Van De Putte, A., et al. (2014). “Conclusions: Present and future of Southern Ocean biogeography,” in Biogeographic Atlas of the Southern Ocean, eds C. De Broyer and P. Koubbi (Cambridge: Scientific Committee on Antarctic Research).

Google Scholar

Kuhner, M. K. (2006). LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics 22, 768–770. doi: 10.1093/bioinformatics/btk051

PubMed Abstract | CrossRef Full Text | Google Scholar

Langenheder, S., and Lindstrom, E. S. (2019). Factors influencing aquatic and terrestrial bacterial community assembly. Environ. Microbiol. Rep. 11, 306–315.

Google Scholar

Langenheder, S., and Lindström, E. S. (2019). Factors influencing aquatic and terrestrial bacterial community assembly. Environ. Microbiol. Rep. 11, 306–315. doi: 10.1111/1758-2229.12731

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkin, A. A., and Martiny, A. C. (2017). Microdiversity shapes the traits, niche space, and biogeography of microbial taxa. Environ. Microbiol. Rep. 9, 55–70. doi: 10.1111/1758-2229.12523

PubMed Abstract | CrossRef Full Text | Google Scholar

Leese, F., Agrawal, S., and Held, C. (2010). Long-distance island hopping without dispersal stages: transportation across major zoogeographic barriers in a Southern Ocean isopod. Naturwissenschaften 97, 583–594. doi: 10.1007/s00114-010-0674-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Leigh, J. W., and Bryant, D. (2015). POPART: full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110–1116. doi: 10.1111/2041-210x.12410

CrossRef Full Text | Google Scholar

Librado, P., and Rozas, J. (2009). DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451–1452. doi: 10.1093/bioinformatics/btp187

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindström, E. S., and Langenheder, S. (2012). Local and regional factors influencing bacterial community assembly. Environ. Microbiol. Rep. 4, 1–9. doi: 10.1111/j.1758-2229.2011.00257.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Locey, K. J. (2010). Synthesizing traditional biogeography with microbial ecology: the importance of dormancy. J. Biogeogr. 37, 1835–1841.

Google Scholar

Logares, R., Deutschmann, I. M., Junger, P. C., Giner, C. R., Krabberod, A. K., Schmidt, T. S. B., et al. (2020). Disentangling the mechanisms shaping the surface ocean microbiota. Microbiome 8:55.

Google Scholar

Logares, R., Tesson, S. V. M., Canback, B., Pontarp, M., Hedlund, K., and Rengefors, K. (2018). Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes. Environ. Microbiol. 20, 2231–2240. doi: 10.1111/1462-2920.14265

PubMed Abstract | CrossRef Full Text | Google Scholar

Louca, S., Parfrey, L. W., and Doebeli, M. (2016). Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272–1277. doi: 10.1126/science.aaf4507

PubMed Abstract | CrossRef Full Text | Google Scholar

Macaya, E. C., and Zuccarello, G. C. (2010). DNA Barcoding and genetic divergence in the Giant Kelp Macrocystis (Laminariales) 1. J. Phycol. 46, 736–742. doi: 10.1111/j.1529-8817.2010.00845.x

CrossRef Full Text | Google Scholar

Martinez Arbizu, P. (2017). pairwiseAdonis: Pairwise Multilevel Comparison Using Adonis. R Package Version 1.

Google Scholar

Martiny, J. B., Bohannan, B. J., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341

PubMed Abstract | CrossRef Full Text | Google Scholar

Martiny, J. B., Eisen, J. A., Penn, K., Allison, S. D., and Horner-Devine, M. C. (2011). Drivers of bacterial β-diversity depend on spatial scale. Proc. Natl. Acad. Sci. 108, 7850–7854. doi: 10.1073/pnas.1016308108

PubMed Abstract | CrossRef Full Text | Google Scholar

Mestre, M., and Höfer, J. (2020). The microbial conveyor belt: connecting the globe through dispersion and dormancy. Trends Microbiol. 29, 482–492. doi: 10.1016/j.tim.2020.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Meysman, F. J., Middelburg, J. J., and Heip, C. H. (2006). Bioturbation: a fresh look at Darwin’s last idea. Trends Ecol. Evol. 21, 688–695. doi: 10.1016/j.tree.2006.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Milici, M., Vital, M., Tomasch, J., Badewien, T. H., Giebel, H. A., Plumeier, I., et al. (2017). Diversity and community composition of particle−associated and free−living bacteria in mesopelagic and bathypelagic Southern Ocean water masses: evidence of dispersal limitation in the Bransfield Strait. Limnol. Oceanogr. 62, 1080–1095. doi: 10.1002/lno.10487

CrossRef Full Text | Google Scholar

Miyazaki, M., Sakai, S., Yamanaka, Y., Saito, Y., Takai, K., and Imachi, H. (2014). Spirochaeta psychrophila sp. nov., a psychrophilic spirochaete isolated from subseafloor sediment, and emended description of the genus Spirochaeta. Int. J. Syst. Evol. Microbiol. 64, 2798–2804. doi: 10.1099/ijs.0.062463-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Moon, K. L., Chown, S. L., and Fraser, C. I. (2017). Reconsidering connectivity in the sub−Antarctic. Biol. Rev. 92, 2164–2181. doi: 10.1111/brv.12327

PubMed Abstract | CrossRef Full Text | Google Scholar

Murray, A. E., and Grzymski, J. J. (2007). Diversity and genomics of Antarctic marine micro-organisms. Philos. Trans. R. Soc. Lond. B Biol. Sci. 362, 2259–2271. doi: 10.1098/rstb.2006.1944

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagylaki, T. (1980). The strong-migration limit in geographically structured populations. J. Math. Biol. 9, 101–114. doi: 10.1007/bf00275916

PubMed Abstract | CrossRef Full Text | Google Scholar

Nemergut, D. R., Schmidt, S. K., Fukami, T., O’neill, S. P., Bilinski, T. M., Stanish, L. F., et al. (2013). Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356.

Google Scholar

Nikula, R., Fraser, C., Spencer, H., and Waters, J. (2010). Circumpolar dispersal by rafting in two subantarctic kelp-dwelling crustaceans. Mar. Ecol. Prog. Ser. 405, 221–230. doi: 10.3354/meps08523

CrossRef Full Text | Google Scholar

Nosil, P., Egan, S. P., and Funk, D. J. (2008). Heterogeneous genomic differentiation between walking-stick ecotypes: “isolation by adaptation” and multiple roles for divergent selection. Evolution 62, 316–336. doi: 10.1111/j.1558-5646.2007.00299.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Obadia, B., Guvener, Z. T., Zhang, V., Ceja-Navarro, J. A., Brodie, E. L., Ja, W. W., et al. (2017). Probabilistic invasion underlies natural gut microbiome stability. Curr. Biol. 27, 1999–2006.e8.

Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., O’hara, R., Simpson, G. L., et al. (2011). Vegan: Community Ecology Package, 2011. R package version, 1.17-18.

Google Scholar

Orsini, L., Vanoverbeke, J., Swillen, I., Mergeay, J., and De Meester, L. (2013). Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol. Ecol. 22, 5983–5999. doi: 10.1111/mec.12561

PubMed Abstract | CrossRef Full Text | Google Scholar

Papke, R. T., Ramsing, N. B., Bateson, M. M., and Ward, D. M. (2003). Geographical isolation in hot spring cyanobacteria. Environ. Microbiol. 5, 650–659. doi: 10.1046/j.1462-2920.2003.00460.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023

PubMed Abstract | CrossRef Full Text | Google Scholar

Paradis, E. (2010). pegas: an R package for population genetics with an integrated–modular approach. Bioinformatics 26, 419–420. doi: 10.1093/bioinformatics/btp696

PubMed Abstract | CrossRef Full Text | Google Scholar

Paradis, E., Blomberg, S., Bolker, B., Brown, J., Claude, J., Cuong, H. S., et al. (2015). Package ‘ape’. Analyses of Phylogenetics and Evolution, version 2, 4-1.

Google Scholar

Pesant, S., Not, F., Picheral, M., Kandels-Lewis, S., Le Bescot, N., Gorsky, G., et al. (2015). Open science resources for the discovery and analysis of Tara Oceans data. Sci. Data 2, 1–16.

Google Scholar

Petro, C., Starnawski, P., Schramm, A., and Kjeldsen, K. U. (2017). Microbial community assembly in marine sediments. Aquat. Microb. Ecol. 79, 177–195. doi: 10.3354/ame01826

CrossRef Full Text | Google Scholar

Phillpot, H. (1985). “Physical geography–climate,” in Key Environments, Antarctica, eds W. N. Bonner and D. W. H. Walton (Oxford: Pergamon Press).

Google Scholar

Poulin, E., and Feral, J.-P. (1995). Pattern of spatial distribution of a brood-protecting schizasterid echinoid, Abatus cordatus, endemic to the Kerguelen Islands. Mar. Ecol. Prog. Ser. 118, 179–186. doi: 10.3354/meps118179

CrossRef Full Text | Google Scholar

Poulin, E., González-Wevar, C., Díaz, A., Gérard, K., and Hüne, M. (2014). Divergence between Antarctic and South American marine invertebrates: what molecular biology tells us about Scotia Arc geodynamics and the intensification of the Antarctic Circumpolar Current. Glob. Planet. Change 123, 392–399. doi: 10.1016/j.gloplacha.2014.07.017

CrossRef Full Text | Google Scholar

R Core Team (2012). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Raes, E. J., Bodrossy, L., Van De Kamp, J., Bissett, A., Ostrowski, M., Brown, M. V., et al. (2018). Oceanographic boundaries constrain microbial diversity gradients in the South Pacific Ocean. Proc. Natl. Acad. Sci. U.S.A. 115, E8266–E8275.

Google Scholar

Ramoneda, J., Le Roux, J. J., Frossard, E., Frey, B., and Gamper, H. A. (2020). Experimental assembly reveals ecological drift as a major driver of root nodule bacterial diversity in a woody legume crop. FEMS Microbiol. Ecol. 96:fiaa083.

Google Scholar

Renelies-Hamilton, J., Germer, K., Sillam-Dussès, D., Bodawatta, K. H., and Poulsen, M. (2021). Disentangling the relative roles of vertical transmission, subsequent colonizations, and diet on cockroach microbiome assembly. mSphere 6, e01023-20.

Google Scholar

Richter-Heitmann, T., Hofner, B., Krah, F. S., Sikorski, J., Wust, P. K., Bunk, B., et al. (2020). Stochastic dispersal rather than deterministic selection explains the spatio-temporal distribution of soil bacteria in a temperate grassland. Front. Microbiol. 11:1391. doi: 10.3389/fmicb.2020.01391

PubMed Abstract | CrossRef Full Text | Google Scholar

Schneider, S., Roessli, D., and Excoffier, L. (2000). Arlequin: a software for population genetics data analysis. User Manual Ver 2, 2496–2497.

Google Scholar

Schwob, G., Cabrol, L., Poulin, E., and Orlando, J. (2020). Characterization of the gut microbiota of the Antarctic heart urchin (Spatangoida) Abatus agassizii. Front. Microbiol. 11:308. doi: 10.3389/fmicb.2020.00308

PubMed Abstract | CrossRef Full Text | Google Scholar

Serebryakova, A., Aires, T., Viard, F., Serrao, E. A., and Engelen, A. H. (2018). Summer shifts of bacterial communities associated with the invasive brown seaweed Sargassum muticum are location and tissue dependent. PLoS One 13:e0206734. doi: 10.1371/journal.pone.0206734

PubMed Abstract | CrossRef Full Text | Google Scholar

Shade, A., Dunn, R. R., Blowes, S. A., Keil, P., Bohannan, B. J., Herrmann, M., et al. (2018). Macroecology to unite all life, large and small. Trends Ecol. Evol. 33, 731–744. doi: 10.1016/j.tree.2018.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Shapiro, B. J., Friedman, J., Cordero, O. X., Preheim, S. P., Timberlake, S. C., Szabó, G., et al. (2012). Population genomics of early events in the ecological differentiation of bacteria. Science 336, 48–51. doi: 10.1126/science.1218198

PubMed Abstract | CrossRef Full Text | Google Scholar

Sintes, E., De Corte, D., Ouillon, N., and Herndl, G. J. (2015). Macroecological patterns of archaeal ammonia oxidizers in the Atlantic Ocean. Mol. Ecol. 24, 4931–4942. doi: 10.1111/mec.13365

PubMed Abstract | CrossRef Full Text | Google Scholar

Stegen, J. C., Lin, X., Fredrickson, J. K., Chen, X., Kennedy, D. W., Murray, C. J., et al. (2013). Quantifying community assembly processes and identifying features that impose them. ISME J. 7, 2069–2079. doi: 10.1038/ismej.2013.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Stegen, J. C., Lin, X., Fredrickson, J. K., and Konopka, A. E. (2015). Estimating and mapping ecological processes influencing microbial community assembly. Front. Microbiol. 6:370.

Google Scholar

Stegen, J. C., Lin, X., Konopka, A. E., and Fredrickson, J. K. (2012). Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J. 6, 1653–1664. doi: 10.1038/ismej.2012.22

PubMed Abstract | CrossRef Full Text | Google Scholar

Sul, W. J., Oliver, T. A., Ducklow, H. W., Amaral-Zettler, L. A., and Sogin, M. L. (2013). Marine bacteria exhibit a bipolar distribution. Proc. Natl. Acad. Sci. U.S.A. 110, 2342–2347. doi: 10.1073/pnas.1212424110

PubMed Abstract | CrossRef Full Text | Google Scholar

Sunagawa, S., Coelho, L. P., Chaffron, S., Kultima, J. R., Labadie, K., Salazar, G., et al. (2015). Structure and function of the global ocean microbiome. Science 348:1261359.

Google Scholar

Szekely, A. J., Berga, M., and Langenheder, S. (2013). Mechanisms determining the fate of dispersed bacterial communities in new environments. ISME J. 7, 61–71. doi: 10.1038/ismej.2012.80

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. R., Pacocha, S., Pharino, C., Klepac-Ceraj, V., Hunt, D. E., Benoit, J., et al. (2005). Genotypic diversity within a natural coastal bacterioplankton population. Science 307, 1311–1313. doi: 10.1126/science.1106028

PubMed Abstract | CrossRef Full Text | Google Scholar

Thornhill, D. J., Mahon, A. R., Norenburg, J. L., and Halanych, K. M. (2008). Open-ocean barriers to dispersal: a test case with the Antarctic Polar Front and the ribbon worm Parborlasia corrugatus (Nemertea: Lineidae). Mol. Ecol. 17, 5104–5117. doi: 10.1111/j.1365-294x.2008.03970.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Trovant, B., Orensanz, J. M., Ruzzante, D. E., Stotz, W., and Basso, N. G. (2015). Scorched mussels (BIVALVIA: MYTILIDAE: BRACHIDONTINAE) from the temperate coasts of South America: phylogenetic relationships, trans-Pacific connections and the footprints of quaternary glaciations. Mol. Phylogenet. Evol. 82(Pt A), 60–74. doi: 10.1016/j.ympev.2014.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Turon, X., Antich, A., Palacín, C., Præbel, K., and Wangensteen, O. S. (2020). From metabarcoding to metaphylogeography: separating the wheat from the chaff. Ecol. Appl. 30:e02036.

Google Scholar

van der Gast, C. J. (2015). Microbial biogeography: the end of the ubiquitous dispersal hypothesis? Environ. Microbiol. 17, 544–546. doi: 10.1111/1462-2920.12635

PubMed Abstract | CrossRef Full Text | Google Scholar

Vavrek, M. J. (2011). Fossil: palaeoecological and palaeogeographical analysis tools. Palaeontol. Electron. 14:16.

Google Scholar

Vellend, M. (2016). The Theory of Ecological Communities. Princeton, NJ: Princeton University Press.

Google Scholar

Wang, I. J. (2013). Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution 67, 3403–3411. doi: 10.1111/evo.12134

PubMed Abstract | CrossRef Full Text | Google Scholar

Whitaker, R. J., Grogan, D. W., and Taylor, J. W. (2003). Geographic barriers isolate endemic populations of hyperthermophilic archaea. Science 301, 976–978. doi: 10.1126/science.1086909

PubMed Abstract | CrossRef Full Text | Google Scholar

White, T. A., Stamford, J., and Rus Hoelzel, A. (2010). Local selection and population structure in a deep−sea fish, the roundnose grenadier (Coryphaenoides rupestris). Mol. Ecol. 19, 216–226. doi: 10.1111/j.1365-294x.2009.04446.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, D., Lauro, F. M., Williams, T. J., Demaere, M. Z., Brown, M. V., Hoffman, J. M., et al. (2013). Biogeographic partitioning of Southern Ocean microorganisms revealed by metagenomics. Environ. Microbiol. 15, 1318–1333.

Google Scholar

Wu, L., Ning, D., Zhang, B., Li, Y., Zhang, P., Shan, X., et al. (2019). Global diversity and biogeography of bacterial communities in wastewater treatment plants. Nat. Microbiol. 4, 1183–1195.

Google Scholar

Zhou, J., and Ning, D. (2017). Stochastic community assembly: does it matter in microbial ecology? Microbiol. Mol. Biol. Rev. 81:e00002-17.

Google Scholar

Zure, M., Fernandez-Guerra, A., Munn, C. B., and Harder, J. (2017). Geographic distribution at subspecies resolution level: closely related Rhodopirellula species in European coastal sediments. ISME J. 11, 478–489.

Google Scholar