Evaluating connectivity in the brooding brittle star Astrotoma agassizii across the drake passage in the Southern Ocean

Surface Bacterioplankton Community Structure Crossing the Antarctic Circumpolar Current Fronts

The Antarctic Circumpolar Current (ACC) is the major current in the Southern Ocean, isolating the warm stratified subtropical waters from the more homogeneous cold polar waters. The ACC flows from west to east around Antarctica and generates an overturning circulation by fostering deep-cold water upwelling and the formation of new water masses, thus affecting the Earth's heat balance and the global distribution of carbon. The ACC is characterized by several water mass boundaries or fronts, known as the Subtropical Front (STF), Subantarctic Front (SAF), Polar Front (PF), and South Antarctic Circumpolar Current Front (SACCF), identified by typical physical and chemical properties. While the physical characteristics of these fronts have been characterized, there is still poor information regarding the microbial diversity of this area. Here we present the surface water bacterioplankton community structure based on 16S rRNA sequencing from 13 stations sampled in 2017 between New Zealand to the Ross Sea crossing the ACC Fronts. Our results show a distinct succession in the dominant bacterial phylotypes present in the different water masses and suggest a strong role of sea surface temperatures and the availability of Carbon and Nitrogen in controlling community composition. This work represents an important baseline for future studies on the response of Southern Ocean epipelagic microbial communities to climate change.

Development of two new sets of PCR primers for eDNA metabarcoding of brittle stars (Echinodermata, Ophiuroidea)

Brittle stars (class Ophiuroidea) are marine invertebrates comprising approximately 2,100 extant species, and are considered to constitute the most diverse taxon of the phylum Echinodermata. As a non-invasive method for monitoring biodiversity, we developed two new sets of PCR primers for metabarcoding environmental DNA (eDNA) from brittle stars. The new primer sets were designed to amplify 2 short regions of the mitochondrial 16S rRNA gene, comprising a conserved region (111–115 bp, 112 bp on average; named “16SOph1”) and a hyper-variable region (180–195 bp, 185 bp on average; named “16SOph2”) displaying interspecific variation. The performance of the primers was tested using eDNA obtained from two sources: a) rearing water of an 2.5 or 170 L aquarium tanks containing 15 brittle star species and b) from natural seawater collected around Misaki, the Pacific coast of central Japan, at depths ranging from shallow (2 m) to deep (> 200 m) sea. To build a reference library, we obtained 16S rRNA sequences of brittle star specimens collected from around Misaki and from similar depths in Japan, and sequences registered in International Nucleotide Sequence Database Collaboration. As a result of comparison of the obtained eDNA sequences with the reference library 37 (including cryptic species) and 26 brittle star species were detected with certain identities by 16SOph1 and 16SOph2 analyses, respectively. In shallow water, the number of species and reads other than the brittle stars detected with 16SOph1 was less than 10% of the total number. On the other hand, the number of brittle star species and reads detected with 16SOph2 was less than half of the total number, and the number of detected non-brittle star metazoan species ranged from 20 to 46 species across 6 to 8 phyla (only the reads at the “Tank” were less than 0.001%). The number of non-brittle star species and reads at 80 m was less than 10% with both of the primer sets. These findings suggest that 16SOph1 is specific to the brittle star and 16SOph2 is suitable for a variety of marine metazoans. It appears, however, that further optimization of primer sequences would still be necessary to avoid possible PCR dropouts from eDNA extracts. Moreover, a detailed elucidation of the brittle star fauna in the examined area, and the accurate identification of brittle star species in the current DNA databank is required.

Multiple Doris “ kerguelenensis ” (Nudibranchia) species span the Antarctic Polar Front

Despite strong historical biogeographical links between benthic faunal assemblages of the Magellan region of South America and the Antarctic Peninsula, very few studies have documented contemporary movement and gene flow in or out of the Southern Ocean, especially across the Antarctic Polar Front (APF). In fact, oceanographic barriers such as the APF and Antarctica's long geologic isolation have substantially separated the continents and facilitated the evolution of endemic marine taxa found within the Antarctic region. The Southern Ocean benthic sea slug complex, Doris “kerguelenensis,” are a group of direct‐developing, simultaneous hermaphrodites that lack a dispersive larval stage. To date, there are 59 highly divergent species known within this complex. Here, we provide evidence to show intraspecific genetic connectivity occurs across the APF for multiple species within the D. “kerguelenensis” nudibranch species complex. We addressed questions of genetic connectivity by examining the phylogeographic structure of the three best‐sampled D. “kerguelenensis” species and another three trans‐APF species using the protein coding mtDNA gene, cytochrome oxidase I. We also highlight alternative refugia uses among species with the same life history traits (i.e., benthic and direct developers) and for some species, extremely large distributions are established (e.g., circumpolarity). By improving our sampling of these nudibranchs, we gain better insight into the population structure and connectivity of the Antarctic region. This work also demonstrates how difficult it is to make generalizations across Antarctic marine species, even among ecologically‐similar, closely related species.

Origin, diversity, and biogeography of Antarctic scale worms (Polychaeta: Polynoidae): a wide-scale barcoding approach

The Antarctic marine environment hosts diversified and highly endemic benthos owing to its unique geologic and climatic history. Current warming trends have increased the urgency of understanding Antarctic species history to predict how environmental changes will impact ecosystem functioning. Antarctic benthic lineages have traditionally been examined under three hypotheses: (1) high endemism and local radiation, (2) emergence of deep‐sea taxa through thermohaline circulation, and (3) species migrations across the Polar Front. In this study, we investigated which hypotheses best describe benthic invertebrate origins by examining Antarctic scale worms (Polynoidae). We amassed 691 polynoid sequences from the Southern Ocean and neighboring areas: the Kerguelen and Tierra del Fuego (South America) archipelagos, the Indian Ocean, and waters around New Zealand. We performed phylogenetic reconstructions to identify lineages across geographic regions, aided by mitochondrial markers cytochrome c oxidase subunit I (Cox1) and 16S ribosomal RNA (16S). Additionally, we produced haplotype networks at the species scale to examine genetic diversity, biogeographic separations, and past demography. The Cox1 dataset provided the most illuminating insights into the evolution of polynoids, with a total of 36 lineages identified. Eunoe sp. was present at Tierra del Fuego and Kerguelen, in favor of the latter acting as a migration crossroads. Harmothoe fuligineum, widespread around the Antarctic continent, was also present but isolated at Kerguelen, possibly resulting from historical freeze–thaw cycles. The genus Polyeunoa appears to have diversified prior to colonizing the continent, leading to the co‐occurrence of at least three cryptic species around the Southern and Indian Oceans. Analyses identified that nearly all populations are presently expanding following a bottleneck event, possibly caused by habitat reduction from the last glacial episodes. Findings support multiple origins for contemporary Antarctic polynoids, and some species investigated here provide information on ancestral scenarios of (re)colonization. First, it is apparent that species collected from the Antarctic continent are endemic, as the absence of closely related species in the Kerguelen and Tierra del Fuego datasets for most lineages argues in favor of Hypothesis 1 of local origin. Next, Eunoe sp. and H. fuligineum, however, support the possibility of Kerguelen and other sub‐Antarctic islands acting as a crossroads for larvae of some species, in support of Hypothesis 3. Finally, the genus Polyeunoa, conversely, is found at depths greater than 150 m and may have a deep origin, in line with Hypothesis 2. These “non endemic” groups, nevertheless, have a distribution that is either north or south of the Antarctic Polar Front, indicating that there is still a barrier to dispersal, even in the deep sea.

Mitochondrial phylogeny of the brittle star genus Ophioderma

We reconstructed the mitochondrial phylogeny of the species of the brittle star genus Ophioderma, using sequences of the Cytochrome Oxidase I gene (COI) to address four questions: (i) Are the species of Ophioderma described on morphological evidence reflected in mitochondrial genealogy? (ii) Which species separated from which? (iii) When did speciation events occur? (iv) What is the rate of COI evolution in ophiuroids? We found that most of the 22 described species we sampled coincide with monophyletic clusters of COI sequences, but there are exceptions, most notably in the eastern Pacific, in which three undescribed species were indicated. The COI phylogeny lacks resolution in the deeper nodes, but it does show that there are four species pairs, the members of which are found on either side of the central American Isthmus. Two pairs with a genetic distance of ~ 4% between Atlantic and Pacific members were probably split during the final stages of Isthmus completion roughly 3 million years ago. The rate of divergence provided by these pairs allowed the calibration of a relaxed molecular clock. Estimated dates of divergence indicate that the lineages leading to extant species coalesce at times much older than congeneric species in other classes of echinoderms, suggesting that low extinction rates may be one of the reasons that ophiuroids are species-rich. The mean rate of COI substitution in Ophioderma is three times slower than that of echinoids. Conclusions of previous mitochondrial DNA studies of ophiuroids that relied on echinoid calibrations to determine divergence times need to be revised.

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.

Divergence time and species delimitation of microbivalves in the Southern Ocean: the case of Kidderia species

The systematics of Subantarctic and Antarctic near-shore marine benthic invertebrates requires major revision and highlights the necessity to incorporate additional sources of information in the specimen identification chart in the Southern Ocean (SO). In this study, we aim to improve our understanding of the biodiversity of Kidderia (Dall 1876) through molecular and morphological comparisons of Antarctic and Subantarctic taxa. The microbivalves of the genus Kidderia are small brooding organisms that inhabit intertidal and shallow subtidal rocky ecosystems. This genus represents an interesting model to test the vicariance and dispersal hypothesis in the biogeography of the SO. However, the description of Kidderia species relies on a few morphological characters and biogeographic records that raise questions about the true diversity in the group. Here we will define the specimens collected with genetic tools, delimiting their respective boundaries across provinces of the SO, validating the presence of two species of Kidderia. Through the revision of taxonomic issues and species delimitation, it was possible to report that the Antarctic species is Kidderia subquadrata and the species recorded in the Subantarctic islands Diego Ramirez, South Georgia and the Kerguelen Archipelago is Kidderia minuta. The divergence time estimation suggests the origin and diversification of Kidderia lineages are related to historical vicariant processes probably associated with the separation of the continental landmasses close to the late Eocene.

An approach using ddRADseq and machine learning for understanding speciation in Antarctic Antarctophilinidae gastropods

Sampling impediments and paucity of suitable material for molecular analyses have precluded the study of speciation and radiation of deep-sea species in Antarctica. We analyzed barcodes together with genome-wide single nucleotide polymorphisms obtained from double digestion restriction site-associated DNA sequencing (ddRADseq) for species in the family Antarctophilinidae. We also reevaluated the fossil record associated with this taxon to provide further insights into the origin of the group. Novel approaches to identify distinctive genetic lineages, including unsupervised machine learning variational autoencoder plots, were used to establish species hypothesis frameworks. In this sense, three undescribed species and a complex of cryptic species were identified, suggesting allopatric speciation connected to geographic or bathymetric isolation. We further observed that the shallow waters around the Scotia Arc and on the continental shelf in the Weddell Sea present high endemism and diversity. In contrast, likely due to the glacial pressure during the Cenozoic, a deep-sea group with fewer species emerged expanding over great areas in the South-Atlantic Antarctic Ridge. Our study agrees on how diachronic paleoclimatic and current environmental factors shaped Antarctic communities both at the shallow and deep-sea levels, promoting Antarctica as the center of origin for numerous taxa such as gastropod mollusks.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Phylogeography and genetic diversity of the microbivalve Kidderia subquadrata, reveals new data from West Antarctic Peninsula

It is well established that Antarctic biodiversity has been strongly influenced by rapid climatic fluctuations during the Quaternary. Marine invertebrates from Antarctica constitute an interesting lens through which to study the impacts of the last glacial periods as glaciation impacted the distribution and intraspecific genetic variation of these animals. However, the impact on the spatial genetic distribution and historical demography of local processes in areas adjacent to the West Antarctic Peninsula (WAP) is less clear. Here we present new genetic information on the bivalve Kidderia subquadrata, a small mollusk that inhabits intertidal rocky island ecosystems throughout the WAP. Using a phylogeographical approach, we examined the spatial patterns of genetic diversity in this brooder species to test the hypothesis of strong genetic structure in incubating organisms and the hypothesis of glacial refugia in organisms with limited dispersion. We found evidence of strong genetic structure among populations of the WAP and a recent expansion in the South Shetland Islands. Our findings are concordant with the predictions that incubating organisms, abundant in Antarctica, present a strong genetic structure among their populations and also support the hypothesis of glacial refugia in organisms with limited dispersion. The effect of the coastal current pattern in the WAP is suggested as a driver to the local spatial dynamics of the genetic diversity distribution. Although genetic information about this microbivalve is still scarce, the knowledge reported here has increased our understanding of the evolutionary patterns of this organism that is endemic to the Southern Ocean.

Cold Adaptation in Antarctic Notothenioids: Comparative Transcriptomics Reveals Novel Insights in the Peculiar Role of Gills and Highlights Signatures of Cobalamin Deficiency

Far from being devoid of life, Antarctic waters are home to Cryonotothenioidea, which represent one of the fascinating cases of evolutionary adaptation to extreme environmental conditions in vertebrates. Thanks to a series of unique morphological and physiological peculiarities, which include the paradigmatic case of loss of hemoglobin in the family Channichthyidae, these fish survive and thrive at sub-zero temperatures. While some of the distinctive features of such adaptations have been known for decades, our knowledge of their genetic and molecular bases is still limited. We generated a reference de novo assembly of the icefish Chionodraco hamatus transcriptome and used this resource for a large-scale comparative analysis among five red-blooded Cryonotothenioidea, the sub-Antarctic notothenioid Eleginops maclovinus and seven temperate teleost species. Our investigation targeted the gills, a tissue of primary importance for gaseous exchange, osmoregulation, ammonia excretion, and its role in fish immunity. One hundred and twenty genes were identified as significantly up-regulated in Antarctic species and surprisingly shared by red- and white-blooded notothenioids, unveiling several previously unreported molecular players that might have contributed to the evolutionary success of Cryonotothenioidea in Antarctica. In particular, we detected cobalamin deficiency signatures and discussed the possible biological implications of this condition concerning hematological alterations and the heavy parasitic loads typically observed in all Cryonotothenioidea.

Larval Dispersal Modeling Suggests Limited Ecological Connectivity Between Fjords on the West Antarctic Peninsula

Larval dispersal is a key process for community assembly and population maintenance in the marine environment, yet it is extremely difficult to measure at ecologically relevant spatio-temporal scales. We used a high-resolution hydrodynamic model and particle-tracking model to explore the dispersal of simulated larvae in a hydrographically complex region of fjords on the West Antarctic Peninsula. Modeled larvae represented two end members of dispersal potential observed in Antarctic benthos resulting from differing developmental periods and swimming behavior. For simulations of low dispersing larvae (pre-competency period = 8 days, settlement period = 15 days, swimming downward) self-recruitment within fjords was important, with no larval settlement occurring in adjacent fjords <50 km apart. For simulations of highly dispersing organisms (pre-competency period = 35-120 days, settlement period = 30-115 days, no swimming behavior), dispersal between fjords occurred when larvae were in the water column for at least 35 days, but settlement was rarely successful even for larvae spending up to 150 days in the plankton. The lack of ecological connectivity between fjords within a single spawning event suggests that these fjords harbor ecologically distinct populations in which self-recruitment may maintain populations, and genetic connectivity between fjords is likely achieved through stepping-stone dispersal. Export of larvae from natal fjord populations to the broader shelf region (>100 km distance) occurred within surface layers (<100 m depth) and was enhanced by episodic katabatic wind events that may be common in glaciomarine fjords worldwide.

Contrasting biogeographical patterns in Margarella (Gastropoda: Calliostomatidae: Margarellinae) across the Antarctic Polar Front

Members of the trochoidean genus Margarella (Calliostomatidae) are broadly distributed across Antarctic and sub-Antarctic ecosystems. Here we used novel mitochondrial and nuclear gene sequences to clarify species boundaries and phylogenetic relationships among seven nominal species distributed on either side of the Antarctic Polar Front (APF). Molecular reconstructions and species-delimitation analyses recognized only four species: M. antarctica (the Antarctic Peninsula), M. achilles (endemic to South Georgia), M. steineni (South Georgia and Crozet Island) and the morphologically variable M. violacea (=M. expansa, M. porcellana and M. pruinosa), with populations in southern South America, the Falkland/Malvinas, Crozet and Kerguelen Islands. Margarella violacea and M. achilles are sister species, closely related to M. steineni, with M. antarctica sister to all these. This taxonomy reflects contrasting biogeographic patterns on either side of the APF in the Southern Ocean. Populations of Margarella north of the APF (M. violacea) showed significant genetic variation but with many shared haplotypes between geographically distant populations. By contrast, populations south of the APF (M. antarctica, M. steineni and M. achilles) exhibited fewer haplotypes and comprised three distinct species, each occurring across a separate geographical range. We hypothesize that the biogeographical differences may be the consequence of the presence north of the APF of buoyant kelps - potential long-distance dispersal vectors for these vetigastropods with benthic-protected development - and their near-absence to the south. Finally, we suggest that the low levels of genetic diversity within higher-latitude Margarella reflect the impact of Quaternary glacial cycles that exterminated local populations during their maxima.

Gene flow in the Antarctic bivalve Aequiyoldia eightsii (Jay, 1839) suggests a role for the Antarctic Peninsula Coastal Current in larval dispersal

The Antarctic Circumpolar Current (ACC) dominates the open-ocean circulation of the Southern Ocean, and both isolates and connects the Southern Ocean biodiversity. However, the impact on biological processes of other Southern Ocean currents is less clear. Adjacent to the West Antarctic Peninsula (WAP), the ACC flows offshore in a northeastward direction, whereas the Antarctic Peninsula Coastal Current (APCC) follows a complex circulation pattern along the coast, with topographically influenced deflections depending on the area. Using genomic data, we estimated genetic structure and migration rates between populations of the benthic bivalve Aequiyoldia eightsii from the shallows of southern South America and the WAP to test the role of the ACC and the APCC in its dispersal. We found strong genetic structure across the ACC (between southern South America and Antarctica) and moderate structure between populations of the WAP. Migration rates along the WAP were consistent with the APCC being important for species dispersal. Along with supporting current knowledge about ocean circulation models at the WAP, migration from the tip of the Antarctic Peninsula to the Bellingshausen Sea highlights the complexities of Southern Ocean circulation. This study provides novel biological evidence of a role of the APCC as a driver of species dispersal and highlights the power of genomic data for aiding in the understanding of the influence of complex oceanographic processes in shaping the population structure of marine species.

Benthic ecoregionalization based on echinoid fauna of the Southern Ocean supports current proposals of Antarctic Marine Protected Areas under IPCC scenarios of climate change

The Southern Ocean (SO) is among the regions on Earth that are undergoing regionally the fastest environmental changes. The unique ecological features of its marine life make it particularly vulnerable to the multiple effects of climate change. A network of Marine Protected Areas (MPAs) has started to be implemented in the SO to protect marine ecosystems. However, considering future predictions of the Intergovernmental Panel on Climate Change (IPCC), the relevance of current, static, MPAs may be questioned under future scenarios. In this context, the ecoregionalization approach can prove promising in identifying well-delimited regions of common species composition and environmental settings. These so-called ecoregions are expected to show similar biotic responses to environmental changes and can be used to define priority areas for the designation of new MPAs and the update of their current delimitation. In the present work, a benthic ecoregionalization of the entire SO is proposed for the first time based on abiotic environmental parameters and the distribution of echinoid fauna, a diversified and common member of Antarctic benthic ecosystems. A novel two-step approach was developed combining species distribution modeling with Random Forest and Gaussian Mixture modeling from species probabilities to define current ecoregions and predict future ecoregions under IPCC scenarios RCP 4.5 and 8.5. The ecological representativity of current and proposed MPAs of the SO is discussed with regard to the modeled benthic ecoregions. In all, 12 benthic ecoregions were determined under present conditions, they are representative of major biogeographic patterns already described. Our results show that the most dramatic changes can be expected along the Antarctic Peninsula, in East Antarctica and the sub-Antarctic islands under both IPCC scenarios. Our results advocate for a dynamic definition of MPAs, they also argue for improving the representativity of Antarctic ecoregions in proposed MPAs and support current proposals of Conservation of Antarctic Marine Living Resources for the creation of Antarctic MPAs.

Dwarf brooder versus giant broadcaster: combining genetic and reproductive data to unravel cryptic diversity in an Antarctic brittle star

Poecilogony, or multiple developmental modes in a single species, is exceedingly rare. Several species described as poecilogenous were later demonstrated to be multiple (cryptic) species with a different developmental mode. The Southern Ocean is known to harbor a high proportion of brooders (Thorson's Rule) but with an increasing number of counter examples over recent years. Here we evaluated poecilogony vs. crypticism in the brittle star Astrotoma agassizii across the Southern Ocean. This species was initially described from South America as a brooder before some pelagic stages were identified in Antarctica. Reproductive and mitochondrial data were combined to unravel geographic and genetic variation of developmental modes. Our results indicate that A. agassizii is composed of seven well-supported and deeply divergent clades (I: Antarctica and South Georgia; II: South Georgia and Sub-Antarctic locations including Kerguelen, Patagonian shelf, and New Zealand; III-VI-VII: Patagonian shelf, IV-V: South Georgia). Two of these clades demonstrated strong size dimorphism when in sympatry and can be linked to differing developmental modes (Clade V: dwarf brooder vs. Clade I: giant broadcaster). Based on their restricted geographic distributions and on previous studies, it is likely that Clades III-VI-VII are brooders. Clade II is composed of different morphological species, A. agassizii and A. drachi, the latter originally used as the outgroup. By integrating morphology, reproductive, and molecular data we conclude that the variation identified in A. agassizii is best described as crypticism rather than poecilogony.

Challenging Dogma Concerning Biogeographic Patterns of Antarctica and the Southern Ocean

Antarctica is enormous, cold, remote, and particularly sensitive to climate change. Most biological research below 60°S has focused on the isolated nature of the biota and how organisms have adapted to the cold and ice. However, biogeographic patterns in Antarctica and the Southern Ocean, and the processes explaining how those patterns came about, still await adequate explanation. Both terrestrial and marine organisms have been influenced by climatic change (e.g., glaciation), physical phenomena (e.g., oceanic currents), and/or potential barriers to gene flow (e.g., steep thermal gradients). Whereas the Antarctic region contains diverse and complex marine communities, terrestrial systems tend to be comparatively simple with limited diversity. Here, we challenge the current dogma used to explain the diversity and biogeographic patterns present in the Antarctic. We assert that relatively modern processes within the last few million years, rather than geo-logical events that occurred in the Eocene and Miocene, account for present patterns of biodiversity in the region. Additionally, reproductive life history stages appear to have little influence in structuring genetic patterns in the Antarctic, as currents and glacial patterns are noted to be more important drivers of organismal patterns of distribution. Finally, we highlight the need for additional sampling, high-throughput genomic approaches, and broad, multinational cooperation for addressing outstanding questions of Antarctic biogeography and biodiversity.

Distinct Oceanic Microbiomes From Viruses to Protists Located Near the Antarctic Circumpolar Current

Microbes occupy diverse ecological niches and only through recent advances in next generation sequencing technologies have the true microbial diversity been revealed. Furthermore, lack of perceivable marine barriers to genetic dispersal (i.e., mountains or islands) has allowed the speculation that organisms that can be easily transported by currents and therefore proliferate everywhere. That said, ocean currents are now commonly being recognized as barriers for microbial dispersal. Here we analyzed samples collected from a total of six stations, four located in the Indian Ocean, and two in the Southern Ocean. Amplicon sequencing was used to characterize both prokaryotic and eukaryotic plankton communities, while shotgun sequencing was used for the combined environmental DNA (eDNA), microbial eDNA (meDNA), and viral fractions. We found that Cyanobacteria dominated the prokaryotic component in the South-West Indian Ocean, while γ-Proteobacteria dominated the South-East Indian Ocean. A combination of γ- and α-Proteobacteria dominated the Southern Ocean. Alveolates dominated almost exclusively the eukaryotic component, with variation in the ratio of Protoalveolata and Dinoflagellata depending on station. However, an increase in haptophyte relative abundance was observed in the Southern Ocean. Similarly, the viral fraction was dominated by members of the order Caudovirales across all stations; however, a higher presence of nucleocytoplasmic large DNA viruses (mainly chloroviruses and mimiviruses) was observed in the Southern Ocean. To our knowledge, this is the first that a statistical difference in the microbiome (from viruses to protists) between the subtropical Indian and Southern Oceans. We also show that not all phylotypes can be found everywhere, and that meDNA is not a suitable resource for monitoring aquatic microbial diversity.

Genetic structure and demographic inference of the regular sea urchin Sterechinus neumayeri (Meissner, 1900) in the Southern Ocean: The role of the last glaciation

One of the most relevant characteristics of the extant Southern Ocean fauna is its resiliency to survive glacial processes of the Quaternary. These climatic events produced catastrophic habitat reductions and forced some marine benthic species to move, adapt or go extinct. The marine benthic species inhabiting the Antarctic upper continental shelf faced the Quaternary glaciations with different strategies that drastically modified population sizes and thus affected the amount and distribution of intraspecific genetic variation. Here we present new genetic information for the most conspicuous regular sea urchin of the Antarctic continental shelf, Sterechinus neumayeri. We studied the patterns of genetic diversity and structure in this broadcast-spawner across three Antarctic regions: Antarctic Peninsula, the Weddell Sea and Adélie Land in East Antarctica. Genetic analyses based on mitochondrial and nuclear markers suggested that S. neumayeri is a single genetic unit around the Antarctic continent. The species is characterized by low levels of genetic diversity and exhibits a typical star-like haplotype genealogy that supports the hypothesis of a single in situ refugium. Based on two mutation rates standardized for this genus, the Bayesian Skyline plot analyses detected a rapid demographic expansion after the Last Glacial Maximum. We propose a scenario of rapid postglacial expansion and recolonization of Antarctic shallow areas from a less ice-impacted refugium where the species survived the LGM. Considering the patterns of genetic diversity and structure recorded in the species, this refugium was probably located in East Antarctica.

Closely related octopus species show different spatial genetic structures in response to the Antarctic seascape

Determining whether comparable processes drive genetic divergence among marine species is relevant to molecular ecologists and managers alike. Sympatric species with similar life histories might be expected to show comparable patterns of genetic differentiation and a consistent influence of environmental factors in shaping divergence. We used microsatellite loci to quantify genetic differentiation across the Scotia Arc in three species of closely related benthic octopods, Pareledone turqueti, P. charcoti, and Adelieledone polymorpha. The relative importance of environmental factors (latitude, longitude, depth, and temperature) in shaping genetic structure was investigated when significant spatial genetic structure was uncovered. Isolated populations of P. turqueti and A. polymorpha at these species' range margins were genetically different to samples close to mainland Antarctica; however, these species showed different genetic structures at a regional scale. Samples of P. turqueti from the Antarctic Peninsula, Elephant Island, and Signy Island were genetically different, and this divergence was associated primarily with sample collection depth. By contrast, weak or nonsignificant spatial genetic structure was evident across the Antarctic Peninsula, Elephant Island, and Signy Island region for A. polymorpha, and slight associations between population divergence and temperature or depth (and/or longitude) were detected. Pareledone charcoti has a limited geographic range, but exhibited no genetic differentiation between samples from a small region of the Scotia Arc (Elephant Island and the Antarctic Peninsula). Thus, closely related species with similar life history strategies can display contrasting patterns of genetic differentiation depending on spatial scale; moreover, depth may drive genetic divergence in Southern Ocean benthos.

Invertebrate population genetics across Earth's largest habitat: The deep-sea floor

Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non-chemosynthetic ecosystems on the deep-sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow-water species. Generally, populations at similar depths were well connected over 100s-1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s-1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean-wide (under 4%), and 48% were Atlantic-focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, "ecosystem engineers" and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single-locus mitochondrial genes revealing a common pattern of non-neutrality, consistent with demographic instability or selective sweeps; similar to deep-sea hydrothermal vent fauna. The absence of a clear difference between vent and non-vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single-locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta-analyses where broad inferences about deep-sea ecology could be made.

Crossing the Divide: Admixture Across the Antarctic Polar Front Revealed by the Brittle Star Astrotoma agassizii

The Antarctic Polar Front (APF) is one of the most well-defined and persistent oceanographic features on the planet and serves as a barrier to dispersal between the Southern Ocean and lower latitudes. High levels of endemism in the Southern Ocean have been attributed to this barrier, whereas the accompanying Antarctic Circumpolar Current (ACC) likely promotes west-to-east dispersal. Previous phylogeographic work on the brittle star Astrotoma agassizii Lyman, 1875 based on mitochondrial genes suggested isolation across the APF, even though populations in both South American waters and the Southern Ocean are morphologically indistinguishable. Here, we revisit this finding using a high-resolution 2b-RAD (restriction-site-associated DNA) single-nucleotide polymorphism (SNP)-based approach, in addition to enlarged mitochondrial DNA data sets (16S rDNA, COI, and COII), for comparison to previous work. In total, 955 biallelic SNP loci confirmed the existence of strongly divergent populations on either side of the Drake Passage. Interestingly, genetic admixture was detected between South America and the Southern Ocean in five individuals on both sides of the APF, revealing evidence of recent or ongoing genetic contact. We also identified two differentiated populations on the Patagonian Shelf with six admixed individuals from these two populations. These findings suggest that the APF is a strong but imperfect barrier. Fluctuations in location and strength of the APF and ACC due to climate shifts may have profound consequences for levels of admixture or endemism in this region of the world.

Locked in the icehouse: Evolution of an endemic Epimeria (Amphipoda, Crustacea) species flock on the Antarctic shelf

The Antarctic shelf's marine biodiversity has been greatly influenced by the climatic and glacial history of the region. Extreme temperature changes led to the extinction of some lineages, while others adapted and flourished. The amphipod genus Epimeria is an example of the latter, being particularly diverse in the Antarctic region. By reconstructing a time-calibrated phylogeny based on mitochondrial (COI) and nuclear (28S and H3) markers and including Epimeria species from all oceans, this study provides a temporal and geographical framework for the evolution of Antarctic Epimeria. The monophyly of this genus is not supported by Bayesian Inference, as Antarctic and non-Antarctic Epimeria form two distinct well-supported clades, with Antarctic Epimeria being a sister clade to two stilipedid species. The monophyly of Antarctic Epimeria suggests that this clade evolved in isolation since its origin. While the precise timing of this origin remains unclear, it is inferred that the Antarctic lineage arose from a late Gondwanan ancestor and hence did not colonize the Antarctic region after the continent broke apart from the other fragments of Gondwanaland. The initial diversification of the clade occurred 38.04Ma (95% HPD [48.46Ma; 28.36Ma]) in a cooling environment. Adaptation to cold waters, along with the extinction of cold-intolerant taxa and resulting ecological opportunities, likely led to the successful diversification of Epimeria on the Antarctic shelf. However, there was neither evidence of a rapid lineage diversification early in the clade's history, nor of any shifts in diversification rates induced by glacial cycles. This suggests that a high turnover rate on the repeatedly scoured Antarctic shelf could have masked potential signals of diversification bursts.

A comparison of genetic connectivity in two deep sea corals to examine whether seamounts are isolated islands or stepping stones for dispersal

Ecological processes in the deep sea are poorly understood due to the logistical constraints of sampling thousands of metres below the ocean’s surface and remote from most land masses. Under such circumstances, genetic data provides unparalleled insight into biological and ecological relationships. We use microsatellite DNA to compare the population structure, reproductive mode and dispersal capacity in two deep sea corals from seamounts in the Southern Ocean. The solitary coral Desmophyllum dianthus has widespread dispersal consistent with its global distribution and resilience to disturbance. In contrast, for the matrix-forming colonial coral Solenosmilia variabilis asexual reproduction is important and the dispersal of sexually produced larvae is negligible, resulting in isolated populations. Interestingly, despite the recognised impacts of fishing on seamount communities, genetic diversity on fished and unfished seamounts was similar for both species, suggesting that evolutionary resilience remains despite reductions in biomass. Our results provide empirical evidence that a group of seamounts can function either as isolated islands or stepping stones for dispersal for different taxa. Furthermore different strategies will be required to protect the two sympatric corals and consequently the recently declared marine reserves in this region may function as a network for D. dianthus, but not for S. variabilis.

Positive selection on sperm ion channels in a brooding brittle star: consequence of life-history traits evolution

Closely related species are key models to investigate mechanisms leading to reproductive isolation and early stages of diversification, also at the genomic level. The brittle star cryptic species complex Ophioderma longicauda encompasses the sympatric broadcast-spawning species C3 and the internal brooding species C5. Here, we used de novo transcriptome sequencing and assembly in two closely related species displaying contrasting reproductive modes to compare their genetic diversity and to investigate the role of natural selection in reproductive isolation. We reconstructed 20 146 and 22 123 genes for C3 and C5, respectively, and characterized a set of 12 229 orthologs. Genetic diversity was 1.5-2 times higher in C3 compared to C5, confirming that species with low parental investment display higher levels of genetic diversity. Forty-eight genes were the targets of positive diversifying selection during the evolution of the two species. Notably, two genes (NHE and TetraKCNG) are sperm-specific ion channels involved in sperm motility. Ancestral sequence reconstructions show that natural selection targeted the two genes in the brooding species. This may result from an adaptation to the novel environmental conditions surrounding sperm in the brooding species, either directly affecting sperm or via an increase in male/female conflict. This phenomenon could have promoted prezygotic reproductive isolation between C3 and C5. Finally, the sperm receptors to egg chemoattractants differed between C3 and C5 in the ligand-binding region. We propose that mechanisms of species-specific gamete recognition in brittle stars occur during sperm chemotaxis (sperm attraction towards the eggs), contrary to other marine invertebrates where prezygotic barriers to interspecific hybridization typically occur before sperm-egg fusion.

Geographic structure in the Southern Ocean circumpolar brittle star Ophionotus victoriae (Ophiuridae) revealed from mtDNA and single-nucleotide polymorphism data

Marine systems have traditionally been thought of as “open” with few barriers to gene flow. In particular, many marine organisms in the Southern Ocean purportedly possess circumpolar distributions that have rarely been well verified. Here, we use the highly abundant and endemic Southern Ocean brittle star Ophionotus victoriae to examine genetic structure and determine whether barriers to gene flow have existed around the Antarctic continent. Ophionotus victoriae possesses feeding planktotrophic larvae with presumed high dispersal capability, but a previous study revealed genetic structure along the Antarctic Peninsula. To test the extent of genetic differentiation within O. victoriae, we sampled from the Ross Sea through the eastern Weddell Sea. Whereas two mitochondrial DNA markers (16S rDNA and COI) were employed to allow comparison to earlier work, a 2b‐RAD single‐nucleotide polymorphism (SNP) approach allowed sampling of loci across the genome. Mitochondrial data from 414 individuals suggested three major lineages, but 2b‐RAD data generated 1,999 biallelic loci that identified four geographically distinct groups from 89 samples. Given the greater resolution by SNP data, O. victoriae can be divided into geographically distinct populations likely representing multiple species. Specific historical scenarios that explain current population structure were examined with approximate Bayesian computation (ABC) analyses. Although the Bransfield Strait region shows high diversity possibly due to mixing, our results suggest that within the recent past, dispersal processes due to strong currents such as the Antarctic Circumpolar Current have not overcome genetic subdivision presumably due to historical isolation, questioning the idea of large open circumpolar populations in the Southern Ocean.

Marked phylogeographic structure of Gentoo penguin reveals an ongoing diversification process along the Southern Ocean

Two main hypotheses have been debated about the biogeography of the Southern Ocean: (1) the Antarctic Polar Front (APF), acting as a barrier between Antarctic and sub-Antarctic provinces, and (2) the Antarctic Circumpolar Current (ACC), promoting gene flow among sub-Antarctic areas. The Gentoo penguin is distributed throughout these two provinces, separated by the APF. We analyzed mtDNA (HVR1) and 12 microsatellite loci of 264 Gentoo penguins, Pygoscelis papua, from 12 colonies spanning from the Western Antarctic Peninsula and the South Shetland Islands (WAP) to the sub-Antarctic Islands (SAI). While low genetic structure was detected among WAP colonies (mtDNA Ф_ST=0.037-0.133; microsatellite F_ST=0.009-0.063), high differentiation was found between all SAI and WAP populations (mtDNA Ф_ST=0.678-0.930; microsatellite F_ST=0.110-0.290). These results suggest that contemporary dispersal around the Southern Ocean is very limited or absent. As predicted, the APF appears to be a significant biogeographical boundary for Gentoo penguin populations; however, the ACC does not promote connectivity in this species. Our data suggest demographic expansion in the WAP during the last glacial maximum (LGM, about 20kya), but stability in SAI. Phylogenetic analyses showed a deep divergence between populations from the WAP and those from the SAI. Therefore, taxonomy should be further revised. The Crozet Islands resulted as a basal clade (3.57Mya), followed by the Kerguelen Islands (2.32Mya) as well as a more recent divergence between the Falkland/Malvinas Islands and the WAP (1.27Mya). Historical isolation, local adaptation, and past climate scenarios of those Evolutionarily Significant Units may have led to different potentials to respond to climate changes.

EmpPrior: using outside empirical data to inform branch-length priors for Bayesian phylogenetics

Background

Branch-length parameters are a central component of phylogenetic models and of intrinsic biological interest. Default branch-length priors in some Bayesian phylogenetic software can be unintentionally informative and lead to branch- and tree-length estimates that are unreasonable. Alternatively, priors may be uninformative, but lead to diffuse posterior estimates. Despite the widespread availability of relevant datasets from other groups, biologists rarely leverage outside information to specify branch-length priors that are specific to the analysis they are conducting.

Results

We developed the software package EmpPrior to facilitate the collection and incorporation of relevant, outside information when setting branch-length priors for phylogenetics. EmpPrior efficiently queries TreeBASE to find data that are similar to focal data, in terms of taxonomic and genetic sampling, and uses them to inform branch-length priors for the focal analysis. EmpPrior consists of two components: EmpPrior-search, written in Java to query TreeBASE, and EmpPrior-fit, written in R to parameterize branch-length distributions. In an example analysis, we show how the use of relevant, outside data is made possible by EmpPrior and improves tree-length estimates from a focal dataset.

Conclusion

EmpPrior is easy to use, fast, and improves both the accuracy and precision of branch-length estimates in many circumstances. While EmpPrior’s focus is on branch lengths, the strategy it employs could easily be extended to address other prior parameterization problems in phylogenetics.

The Antarctic Circumpolar Current as a diversification trigger for deep-sea octocorals

BACKGROUND

Antarctica is surrounded by the Antarctic Circumpolar Current (ACC), the largest and strongest current in the world. Despite its potential importance for shaping biogeographical patterns, the distribution and connectivity of deep-sea populations across the ACC remain poorly understood. In this study we conducted the first assessment of phylogeographical patterns in deep-sea octocorals in the South Pacific and Southern Ocean, specifically a group of closely related bottlebrush octocorals (Primnoidae: Tokoprymno and Thourella), as a test case to study the effect of the ACC on the population structure of brooding species. We assessed the degree to which the ACC constitutes a barrier to gene flow between northern and southern populations and whether the onset of diversification of these corals coincides with the origin of the ACC (Oligocene-Miocene boundary).

RESULTS

Based on DNA sequences of two nuclear genes from 80 individuals and a combination of phylogeographic model-testing approaches we found a phylogenetic break corresponding to the spatial occurrence of the ACC. We also found significant genetic structure among our four regional populations. However, we uncovered shared haplotypes among certain population pairs, suggesting long-distance, asymmetrical migration. Our divergence time analyses indicated that the separation of amphi-ACC populations took place during the Middle Miocene around 12.6 million years ago, i.e., after the formation of the ACC.

CONCLUSION

We suggest that the ACC constitutes a semi-permeable barrier to these deep-sea octocorals capable of separating and structuring populations, while allowing short periods of gene flow. The fluctuations in latitudinal positioning of the ACC during the Miocene likely contributed to the diversification of these octocorals. Additionally, we provide evidence that the populations from each of our four sampling regions could actually constitute different species.

Widespread intra-specific genetic homogeneity of coastal Antarctic benthic foraminifera

Benthic foraminifera are a major component of the Antarctic biota. Coastal foraminiferal morphospecies are widely distributed in Antarctic waters. The question is whether these morphotypes are genetically identical or, rather, they represent a cohort of cryptic species. Here, we compared genetically nine benthic foraminiferal morphospecies from Admiralty Bay (South Shetlands) and the western Ross Sea (McMurdo Sound, Terra Nova Bay), separated by a distance of ~4500 km. Additionally, for three of these morphospecies, we included specimens from Rothera (Marguerite Bay), which is located between the two main areas of interest. Our study, based on SSU and ITS rDNA sequence data, shows that all examined morphospecies share the same genotypes despite the presence of considerable intra-individual genetic variability.

Genetic data, reproduction season and reproductive strategy data support the existence of biological species in Ophioderma longicauda

Cryptic species are numerous in the marine environment. The brittle star Ophioderma longicauda is composed of six mitochondrial lineages, encompassing brooders, which form a monophyletic group, and broadcasters, from which the brooders are derived. To clarify the species limits within O. longicauda, we compared the reproductive status of the sympatric lineages L1 and L3 (defined after sequencing a portion of the mitochondrial gene COI) during the month of May in Greece. In addition, we genotyped a nuclear marker, intron i51. Each L3 female was brooding, whereas all L1 specimens displayed full gonads, suggesting temporal pre-zygotic isolation between brooders and broadcasters. Statistical differences were found among lineages in morphology and bathymetric distribution. Finally, the intron i51 was polymorphic in L1 (60 individuals), but monomorphic in L3 (109 individuals), confirming the absence of gene flow between brooders and broadcasters. In conclusion, the broadcasting lineage L1 and the brooding lineage L3 are different biological species.

A hierarchical classification of benthic biodiversity and assessment of protected areas in the Southern Ocean

An international effort is underway to establish a representative system of marine protected areas (MPAs) in the Southern Ocean to help provide for the long-term conservation of marine biodiversity in the region. Important to this undertaking is knowledge of the distribution of benthic assemblages. Here, our aim is to identify the areas where benthic marine assemblages are likely to differ from each other in the Southern Ocean including near-shore Antarctica. We achieve this by using a hierarchical spatial classification of ecoregions, bathomes and environmental types. Ecoregions are defined according to available data on biogeographic patterns and environmental drivers on dispersal. Bathomes are identified according to depth strata defined by species distributions. Environmental types are uniquely classified according to the geomorphic features found within the bathomes in each ecoregion. We identified 23 ecoregions and nine bathomes. From a set of 28 types of geomorphic features of the seabed, 562 unique environmental types were classified for the Southern Ocean. We applied the environmental types as surrogates of different assemblages of biodiversity to assess the representativeness of existing MPAs. We found that 12 ecoregions are not represented in MPAs and that no ecoregion has their full range of environmental types represented in MPAs. Current MPA planning processes, if implemented, will substantially increase the representation of environmental types particularly within 8 ecoregions. To meet internationally agreed conservation goals, additional MPAs will be needed. To assist with this process, we identified 107 spatially restricted environmental types, which should be considered for inclusion in future MPAs. Detailed supplementary data including a spatial dataset are provided.