Climate-driven genetic divergence of limpets with different life histories across a southeast African marine biogeographic disjunction: different processes, same outcome

Mind the Gap: A Review of Disjunctions in Coastal Marine Species

Many coastal marine species have discontinuous distributions or genetic breakpoints throughout their geographical ranges. These spatial and genetic disjunctions occur in species that span limited to broad dispersal potential. Thus, the mechanisms that underlie these disjunctions remain speculative or incompletely known, particularly on small spatial scales where long-term historical processes are unlikely to be the only mechanism contributing to disjunction. Rather, ecological or oceanographic factors may be important. To identify key drivers of coastal disjunctions, we reviewed publications investigating spatial and genetic disjunctions in coastal marine species and visually summarized where and why they are thought to occur. The most frequently cited mechanisms implicated in causing disjunctions include historical processes, oceanographic features, heterogeneous habitat, species introductions, and limited larval dispersal capacities. However, the relative importance of each of these processes varies depending on the spatial scales investigated. Furthermore, locations associated with disjunctions for a suite of species are typically associated with multiple processes that maintain these disjunctions. This study provides a non-exhaustive synthesis of disjunctions in coastal marine species by visualizing where they occur, exploring underlying mechanisms, and investigating biases in how the scientific community studies this phenomenon.

Identifying Whitemouth Croaker (Micropogonias furnieri) Populations along the Rio de Janeiro Coast, Brazil, through Microsatellite and Otolith Analyses

The inshore area of the Southwestern Atlantic between 22 °S and 29 °S (South Brazilian Bight) is a transitional climatic zone, where the tropical and warm temperate provinces mix. In its northern part, i.e., in the coastal waters of Rio de Janeiro, Brazil, local oceanographic conditions, such as upwelling in the north, and great bays with different degrees of anthropogenic influences in the center and south can determine the population structure of several fish stocks. The Whitemouth croaker (Micropogonias furnieri) is one the most heavily exploited fishing resources in this area, but there are still some doubts about its population structure. In this study, through combined analyses using nuclear genetic markers and morphological and geochemical signatures of otoliths, a divergence of individuals between two populations was identified using microsatellites, while a finer spatial structure with three populations (north, center and south, respectively) was found based on otolith shapes and elemental signatures. This regional population structure may have direct implications for rational fisheries management and conservation of the species.

Quaternary ice sheets and sea level regression drove divergence in a marine gastropod along Eastern and Western coasts of South America

The southern coastline of South America is a remarkable area to evaluate how Quaternary glacial processes impacted the demography of the near-shore marine biota. Here we present new phylogeographic analyses in the pulmonate Siphonaria lessonii across its distribution, from northern Chile in the Pacific to Uruguay in the Atlantic. Contrary to our expectations, populations from the southwestern Atlantic, an area that was less impacted by ice during glacial maxima, showed low genetic diversity and evidence of recent expansion, similar to the patterns recorded in this study across heavily ice-impacted areas in the Pacific Magellan margin. We propose that Atlantic and Pacific shallow marine hard-substrate benthic species were both affected during the Quaternary in South America, but by different processes. At higher latitudes of the southeast Pacific, ice-scouring drastically affected S. lessonii populations compared to non-glaciated areas along the Chile-Peru province where the species was resilient. In the southwest Atlantic, S. lessonii populations would have been dramatically impacted by the reduction of near-shore rocky habitat availability as a consequence of glacio-eustatic movements. The increase of gravelly and rocky shore substrates in the southwest Atlantic supports a hypothesis of glacial refugia from where the species recolonized lower latitudes across the Atlantic and Pacific margins. Our results suggest that current patterns of genetic diversity and structure in near-shore marine benthic species do not solely depend on the impact of Quaternary glacial ice expansions but also on the availability of suitable habitats and life-history traits, including developmental mode, bathymetry and the likelihood of dispersal by rafting.

Genetic and morphological divergence at a biogeographic break in the beach-dwelling brooder Excirolana hirsuticauda Menzies (Crustacea, Peracarida)

BACKGROUND

There is a biogeographic break located at 30°S in the southeast Pacific, in a coastal area of strong environmental discontinuities. Several marine benthic taxa with restricted dispersal have a coincident phylogeographic break at 30°S, indicating that genetic structure is moulded by life history traits that limit gene flow and thereby promote divergence and speciation. In order to evaluate intraspecific divergence at this biogeographic break, we investigated the genetic and morphological variation of the directly developing beach isopod Excirolana hirsuticauda along 1900 km of the southeast Pacific coast, across 30°S.

RESULTS

The COI sequences and microsatellite data both identified a strong discontinuity between populations of E. hirsuticauda to the north and south of 30°S, and a second weaker phylogeographic break at approximately 35°S. The three genetic groups were evidenced by different past demographic and genetic diversity signatures, and were also clearly distinguished with microsatellite data clustering. The COI sequences established that the genetic divergence of E. hirsuticauda at 30°S started earlier than divergence at 35°. Additionally, the three groups have different past demographic signatures, with probable demographic expansion occurring earlier in the southern group (south of 35°S), associated with Pleistocene interglacial periods. Interestingly, body length, multivariate morphometric analyses, and the morphology of a fertilization-related morphological character in males, the appendix masculina, reinforced the three genetic groups detected with genetic data.

CONCLUSIONS

The degree of divergence of COI sequences, microsatellite data, and morphology was concordant and showed two geographic areas in which divergence was promoted at differing historical periods. Variation in the appendix masculina of males has probably promoted reproductive isolation. This variation together with gene flow restrictions promoted by life history traits, small body size, oceanographic discontinuities and sandy-beach habitat continuity, likely influenced species divergence at 30°S in the southeast Pacific coast. The degree of genetic and morphological differentiation of populations to the north and south of 30°S suggests that E. hirsuticauda harbours intraspecific divergence consistent with reproductive isolation and an advanced stage of speciation. The speciation process within E. hirsuticauda has been shaped by both restrictions to gene flow and a prezygotic reproductive barrier.

Thermal selection as a driver of marine ecological speciation

Intraspecific genetic structure in widely distributed marine species often mirrors the boundaries between temperature-defined bioregions. This suggests that the same thermal gradients that maintain distinct species assemblages also drive the evolution of new biodiversity. Ecological speciation scenarios are often invoked to explain such patterns, but the fact that adaptation is usually only identified when phylogenetic splits are already evident makes it impossible to rule out the alternative scenario of allopatric speciation with subsequent adaptation. We integrated large-scale genomic and environmental datasets along one of the world's best-defined marine thermal gradients (the South African coastline) to test the hypothesis that incipient ecological speciation is a result of divergence linked to the thermal environment. We identified temperature-associated gene regions in a coastal fish species that is spatially homogeneous throughout several temperature-defined biogeographic regions based on selectively neutral markers. Based on these gene regions, the species is divided into geographically distinct regional populations. Importantly, the ranges of these populations are delimited by the same ecological boundaries that define distinct infraspecific genetic lineages in co-distributed marine species, and biogeographic disjunctions in species assemblages. Our results indicate that temperature-mediated selection represents an early stage of marine ecological speciation in coastal regions that lack physical dispersal barriers.

Seascape genetics of the spiny lobster Panulirus homarus in the Western Indian Ocean: Understanding how oceanographic features shape the genetic structure of species with high larval dispersal potential

This study examines the fine-scale population genetic structure and phylogeography of the spiny lobster Panulirus homarus in the Western Indian Ocean. A seascape genetics approach was used to relate the observed genetic structure based on 21 microsatellite loci to ocean circulation patterns, and to determine the influence of latitude, sea surface temperature (SST), and ocean turbidity (KD490) on population-level processes. At a geospatial level, the genetic clusters recovered corresponded to three putative subspecies, P. h. rubellus from the SW Indian Ocean, P. h. megasculptus from the NW Indian Ocean, and P. h. homarus from the tropical region in-between. Virtual passive Lagrangian particles advected using satellite-derived ocean surface currents were used to simulate larval dispersal. In the SW Indian Ocean, the dispersion of particles tracked over a 4-month period provided insight into a steep genetic gradient observed at the Delagoa Bight, which separates P. h. rubellus and P. h. homarus. South of the contact zone, particles were advected southwestwards by prevailing boundary currents or were retained in nearshore eddies close to release locations. Some particles released in southeast Madagascar dispersed across the Mozambique Channel and reached the African shelf. Dispersal was characterized by high seasonal and inter-annual variability, and a large proportion of particles were dispersed far offshore and presumably lost. In the NW Indian Ocean, particles were retained within the Arabian Sea. Larval retention and self-recruitment in the Arabian Sea could explain the recent genetic divergence between P. h. megasculptus and P. h. homarus. Geographic distance and minimum SST were significantly associated with genetic differentiation in multivariate analysis, suggesting that larval tolerance to SST plays a role in shaping the population structure of P. homarus.

Mitochondrial DNA is unsuitable to test for isolation by distance

Tests for isolation by distance (IBD) are the most commonly used method of assessing spatial genetic structure. Many studies have exclusively used mitochondrial DNA (mtDNA) sequences to test for IBD, but this marker is often in conflict with multilocus markers. Here, we report a review of the literature on IBD, with the aims of determining (a) whether significant IBD is primarily a result of lumping spatially discrete populations, and (b) whether microsatellite datasets are more likely to detect IBD when mtDNA does not. We also provide empirical data from four species in which mtDNA failed to detect IBD by comparing these with microsatellite and SNP data. Our results confirm that IBD is mostly found when distinct regional populations are pooled, and this trend disappears when each is analysed separately. Discrepancies between markers were found in almost half of the studies reviewed, and microsatellites were more likely to detect IBD when mtDNA did not. Our empirical data rejected the lack of IBD in the four species studied, and support for IBD was particularly strong for the SNP data. We conclude that mtDNA sequence data are often not suitable to test for IBD, and can be misleading about species' true dispersal potential. The observed failure of mtDNA to reliably detect IBD, in addition to being a single-locus marker, is likely a result of a selection-driven reduction in genetic diversity obscuring spatial genetic differentiation.

Evidence for panmixia despite barriers to gene flow in the hooked mussel, Ischadium recurvum (Mytilidae; Brachidontinae) along the North American coastline

The discovery of phylogeographic patterns within broadly distributed marine species can be particularly challenging because absolute physical barriers to dispersal can be inconspicuous. Genetic boundaries often lie where ocean currents meet, forming sharp physical and ecological gradients, which may act as barriers to successful migrants. In eastern North America, coastal species often show phylogeographic differentiation associated with two recognized genetic barriers: the Gulf/Atlantic and the Virginia/Carolina discontinuities. We examined 185 specimens of the intertidal hooked mussel Ischadium recurvum collected from 15 locations along the eastern coastline of North America to examine phylogeographic, migration and historical demographic patterns associated climate change associated with Pleistocene glacial patterns. Hypothesis testing using Bayes factors in Migrate-n rejected the presence of phylogeographic breaks consistent with either maritime discontinuity and favoured a panmictic population model. The migration rate from the Gulf to the Atlantic was approximately three times higher than the migration from the Atlantic to the Gulf whereas the Carolina-Virginia migration rates were nearly equal. The summary statistics (Tajima's D, Fu's Fs) were significant and the demographic analyses (mismatch distributions, Bayesian skyline plot) were consistent with patterns of population expansion following glacial retreat during the Pleistocene epoch.

Causations of phylogeographic barrier of some rocky shore species along the Chinese coastline

BACKGROUND

Substrate, ocean current and freshwater discharge are recognized as important factors that control the larval dispersal and recruitment of intertidal species. Life history traits of individual species will determine the differential responses to these physical factors, and hence resulting in contrasting phylogeography across the same biogeographic barrier. To determine how these factors affect genetic structure of rocky shore species along the China coast, a comparative phylogeographic study of four intertidal and subtidal species was conducted using mitochondrial and nuclear DNA by combining new sequences from Siphonaria japonica with previously published sequences from three species (Cellana toreuma, Sargassum horneri and Atrina pectinata).

RESULTS

Analysis of molecular variance and pairwise ΦST revealed significant genetic differences between the Yellow Sea (YS) and the other two marginal seas (East China Sea, ECS and South China Sea, SCS) for rocky-shore species (S. japonica, C. toreuma, S. horneri), but not for muddy-shore species Atrina pectinata. Demographic history analysis proved that the population size of all these four species were persistent though the Last Glacial Maximum (LGM, ~20 ka BP). Migration analysis revealed that gene flow differentiated northward and southward migration for these four species. However, the inferred direction of gene flow using alternatively mitochondrial or nuclear markers was contradictory in S. japonica.

CONCLUSIONS

It is concluded that there is a phylogeographical break at the Yangtze River estuary for the rocky shore species and the causation of the barrier is mainly due to the unsuitable substratum and freshwater discharge. All four intertidal and subtidal species appear to have persisted through the LGM in China, indicating the lower impact of LGM on intertidal and subtidal species than generally anticipated. The imbalanced gene flow between YS and ESCS groups for these four species could be explained by historical refugia. The discordance between mitochondrial and nuclear markers in the MIGRATE analysis of S. japonica prove the importance of employing multi-locus data in biogeographic study. Climate change, land reclamation and dam construction, which are changing substrate and hydrological conditions around Yangtze River estuary, will consequently affect the biogeographic pattern of intertidal species.

Genetic structuring across marine biogeographic boundaries in rocky shore invertebrates

Biogeography investigates spatial patterns of species distribution. Discontinuities in species distribution are identified as boundaries between biogeographic areas. Do these boundaries affect genetic connectivity? To address this question, a multifactorial hierarchical sampling design, across three of the major marine biogeographic boundaries in the central Mediterranean Sea (Ligurian-Tyrrhenian, Tyrrhenian-Ionian and Ionian-Adriatic) was carried out. Mitochondrial COI sequence polymorphism of seven species of Mediterranean benthic invertebrates was analysed. Two species showed significant genetic structure across the Tyrrhenian-Ionian boundary, as well as two other species across the Ionian Sea, a previously unknown phylogeographic barrier. The hypothesized barrier in the Ligurian-Tyrrhenian cannot be detected in the genetic structure of the investigated species. Connectivity patterns across species at distances up to 800 km apart confirmed that estimates of pelagic larval dispersal were poor predictors of the genetic structure. The detected genetic discontinuities seem more related to the effect of past historical events, though maintained by present day oceanographic processes. Multivariate statistical tools were used to test the consistency of the patterns across species, providing a conceptual framework for across-species barrier locations and strengths. Additional sequences retrieved from public databases supported our findings. Heterogeneity of phylogeographic patterns shown by the 7 investigated species is relevant to the understanding of the genetic diversity, and carry implications for conservation biology.

Can novel genetic analyses help to identify low-dispersal marine invasive species?

Genetic methods can be a powerful tool to resolve the native versus introduced status of populations whose taxonomy and biogeography are poorly understood. The genetic study of introduced species is presently dominated by analyses that identify signatures of recent colonization by means of summary statistics. Unfortunately, such approaches cannot be used in low-dispersal species, in which recently established populations originating from elsewhere in the species' native range also experience periods of low population size because they are founded by few individuals. We tested whether coalescent-based molecular analyses that provide detailed information about demographic history supported the hypothesis that a sea squirt whose distribution is centered on Tasmania was recently introduced to mainland Australia and New Zealand through human activities. Methods comparing trends in population size (Bayesian Skyline Plots and Approximate Bayesian Computation) were no more informative than summary statistics, likely because of recent intra-Tasmanian dispersal. However, IMa2 estimates of divergence between putatively native and introduced populations provided information at a temporal scale suitable to differentiate between recent (potentially anthropogenic) introductions and ancient divergence, and indicated that all three non-Tasmanian populations were founded during the period of European settlement. While this approach can be affected by inaccurate molecular dating, it has considerable (albeit largely unexplored) potential to corroborate nongenetic information in species with limited dispersal capabilities.

Phylogeographic structure in benthic marine invertebrates of the southeast Pacific coast of Chile with differing dispersal potential

The role of dispersal potential on phylogeographic structure, evidenced by the degree of genetic structure and the presence of coincident genetic and biogeographic breaks, was evaluated in a macrogeographic comparative approach along the north-central coast of Chile, across the biogeographic transition zone at 30°S. Using 2,217 partial sequences of the mitochondrial Cytochrome Oxidase I gene of eight benthic invertebrate species along ca. 2,600 km of coast, we contrasted dispersal potential with genetic structure and determined the concordance between genetic divergence between biogeographic regions and the biogeographic transition zone at 30°S. Genetic diversity and differentiation highly differed between species with high and low dispersal potential. Dispersal potential, sometimes together with biogeographic region, was the factor that best explained the genetic structure of the eight species. The three low dispersal species, and one species assigned to the high dispersal category, had a phylogeographic discontinuity coincident with the biogeographic transition zone at 30°S. Furthermore, coalescent analyses based on the isolation-with-migration model validate that the split between biogeographic regions north and south of 30°S has a historic origin. The signatures of the historic break in high dispersers is parsimoniously explained by the homogenizing effects of gene flow that have erased the genetic signatures, if ever existed, in high dispersers. Of the four species with structure across the break, only two had significant albeit very low levels of asymmetric migration across the transition zone. Historic processes have led to the current biogeographic and phylogeographic structure of marine species with limited dispersal along the north-central coast of Chile, with a strong lasting impact in their genetic structure.

Triton's trident: cryptic Neogene divergences in a marine clam (Lasaea australis) correspond to Australia's three temperate biogeographic provinces

The southern coast of Australia is composed of three distinct biogeographic provinces distinguished primarily by intertidal community composition. Several ecological mechanisms have been proposed to explain their formation and persistence, but no consensus has been reached. The marine clam Lasaea australis is arguably the most common bivalve on southern Australian rocky shores and occurs in all three provinces. Here, we tested if this species exhibits cryptic genetic structuring corresponding to the provinces and if so, what mechanisms potentially drove its divergence. Variation in two mitochondrial genes (16S and COIII) and one nuclear gene (ITS2) was assayed to test for genetic structuring and to reconstruct the clam's phylogenetic history. Our results showed that L. australis is comprised of three cryptic mitochondrial clades, each corresponding almost perfectly to one of the three biogeographic provinces. Divergence time estimates place their cladogenesis in the Neogene. The trident-like topology and Neogene time frame of L. australis cladogenesis are incongruent with Quaternary vicariance predictions: a two-clade topology produced by Pleistocene Bass Strait land bridge formation. We hypothesize that the interaction of the Middle Miocene Climate Transition with the specific geography of the southern coastline of Australia was the primary cladogenic driver in this clam lineage. Additional in-depth studies of the endemic southern Australian marine biota across all three provinces are needed to establish the generality of this proposed older framework for regional cladogenesis.