Phylogeography and population structure of thornback rays (Raja clavata L., Rajidae)

Rapid diversification of grey mangroves (Avicennia marina) driven by geographic isolation and extreme environmental conditions in the Arabian Peninsula

Biological systems occurring in ecologically heterogeneous and spatially discontinuous habitats provide an ideal opportunity to investigate the relative roles of neutral and selective factors in driving lineage diversification. The grey mangroves (Avicennia marina) of Arabia occur at the northern edge of the species' range and are subject to variable, often extreme, environmental conditions, as well as historic large fluctuations in habitat availability and connectivity resulting from Quaternary glacial cycles. Here, we analyse fully sequenced genomes sampled from 19 locations across the Red Sea, the Arabian Sea and the Persian/Arabian Gulf (PAG) to reconstruct the evolutionary history of the species in the region and to identify adaptive mechanisms of lineage diversification. Population structure and phylogenetic analyses revealed marked genetic structure correlating with geographic distance and highly supported clades among and within the seas surrounding the Arabian Peninsula. Demographic modelling showed times of divergence consistent with recent periods of geographic isolation and low marine connectivity during glaciations, suggesting the presence of (cryptic) glacial refugia in the Red Sea and the PAG. Significant migration was detected within the Red Sea and the PAG, and across the Strait of Hormuz to the Arabian Sea, suggesting gene flow upon secondary contact among populations. Genetic-environment association analyses revealed high levels of adaptive divergence and detected signs of multi-loci local adaptation driven by temperature extremes and hypersalinity. These results support a process of rapid diversification resulting from the combined effects of historical factors and ecological selection and reveal mangrove peripheral environments as relevant drivers of lineage diversity.

Measurement of the Growth of the Main Commercial Rays (Raja clavata, Raja brachyura, Torpedo marmorata, Dipturus oxyrinchus) in European Waters Using Intercalibration Methods

The intercalibration of age readings represents a crucial step in the ageing procedure; the use of different sampling methods, structures, preparation techniques, and ageing criteria can significantly affect age and growth data. This study evaluated the precision and accuracy of ageing for the most important North Atlantic (NA) and Mediterranean (M) ray species, Raja clavata, Raja brachyura, Torpedo marmorata, and Dipturus oxyrinchus, through exchange exercises carried out by readers from different laboratories. In addition, growth parameters were estimated from the obtained data. A total of 663 individual batoids were analysed. R. clavata and R. brachyura samples were obtained from both the NA and the M, while vertebral centra of T. marmorata and D. oxyrinchus were only available for the M. High reading variability was observed for all four evaluated species in terms of CV, APE, and PA. D. oxyrinchus and T. marmorata showed relatively slow growth and the von Bertalanffy model with fixed t0 and Gompertz's model were, respectively, the most precise models for each of these species. In R. brachyura, females had a faster growth rate compared to combined sexes. The vbt0p proved the most precise model for describing growth in this species, and no statistical differences were found between the NO and the M. For R. clavata, the best-fitting model was the vbt0p for females and males in the NO and for females from the M, while the best-fitting model for males from the M and sexes combined for both areas was log.p. Distinct growth patterns were observed between the two study areas.

Matching an Old Marine Paradigm: Limitless Connectivity in a Deep-Water Fish over a Large Distance

Investigations of population structuring in wild species are fundamental to complete the bigger picture defining their ecological and biological roles in the marine realm, to estimate their recovery capacity triggered by human disturbance and implement more efficient management strategies for fishery resources. The Blackspot Seabream (Pagellus bogaraveo, Brünnich 1768) is a commercially valuable deep-water fish highly exploited over past decades. Considering its exploitation status, deepening the knowledge of intraspecific variability, genetic diversity, and differentiation using high-performing molecular markers is considered an important step for a more effective stock assessment and fishery management. With one of the largest efforts conceived of and completed by countries overlooking the Atlantic and Mediterranean coasts in recent years, a total of 320 individuals were collected from different fishing grounds in the Mediterranean Sea and Atlantic Ocean and analysed using 29 microsatellite loci. We applied multiple statistical approaches to investigate the species' connectivity and population structure across most of its described distribution area. Considering the incomplete knowledge regarding the migratory behaviour of adults, here we suggest the importance of egg and larval dispersal in sustaining the observed genetic connectivity on such a large geographical scale.

To Be, or Not to Be: That Is the Hamletic Question of Cryptic Evolution in the Eastern Atlantic and Mediterranean Raja miraletus Species Complex

Despite a high species diversity, skates (Rajiformes) exhibit remarkably conservative morphology and ecology. Limited trait variations occur within and between species, and cryptic species have been reported among sister and non-sister taxa, suggesting that species complexes may be subject to stabilising selection. Three sibling species are currently recognised in the Raja miraletus complex: (i) R. miraletus occurring along the Portuguese and Mediterranean coasts, (ii) R. parva in the Central-Eastern Atlantic off West Africa and (iii) R. ocellifera in the Western Indian Ocean off South Africa. In the present study, the genetic variation at mitochondrial and nuclear markers was estimated in the species complex by analysing 323 individuals sampled across most of its geographical distribution area to test the hypothesis that restricted gene flow and genetic divergence within species reflect known climate and bio-oceanographic discontinuities. Our results support previous morphological studies and confirm the known taxonomic boundaries of the three recognised species. In addition, we identified multiple weakly differentiated clades in the Northeastern Atlantic Ocean and Mediterranean, at least two additional cryptic taxa off Senegal and Angola, a pronounced differentiation of ancient South African clades. The hidden genetic structure presented here may represent a valuable support to species' conservation action plans.

Trans-Arctic vicariance in Strongylocentrotus sea urchins

The sea urchins Strongylocentotus pallidus and S. droebachiensis first invaded the Atlantic Ocean from the Pacific following the opening of the Bering seaway in the late Miocene. While trans-Arctic dispersal during the Pleistocene is thought to have maintained species' integrity, a recent genomic analysis identified a reproductively isolated cryptic species within S. droebachiensis. Based on previous studies, the distribution of one of these lineages (S. droebachiensis W) includes the shallow water habitats of the northwest Atlantic and Pacific, while the other (S. droebachiensis E) is found throughout the shallow habitat in the northeast but is mostly restricted to deep habitats (>65 m) in the northwest Atlantic. However, since genetic variation within S. droebachiensis has been largely unstudied in the north Pacific and Arctic oceans, the biogeography of the cryptic species is not well known, and it is difficult to identify the mechanisms driving population subdivision and speciation. Here we use population genetic analyses to characterize the distribution of each species, and to test hypotheses about the role of vicariance in the evolution of systematic and genomic divergence within the genus. We collected individuals of all three Strongylocentrotus species (n = 365) from 10 previously unsampled locations in the northeast Pacific and north Atlantic (Labrador Sea and Norway), and generated mtDNA sequence data for a 418 bp fragment of cytochrome c oxidase subunit I (COI). To assess the biogeography of all three species, we combined our alignment with five previously published data sets (total n = 789) and used statistical parsimony and maximum likelihood to identify species and characterize their distribution within and among oceans. Patterns of haplotype sharing, pairwise F _ST , and hierarchical analyses of molecular variance (AMOVA) identified trans-Arctic dispersal in S. pallidus and S. droebachiensis W, but other than 5 previously reported singletons we failed to detect additional mtDNA haplotypes of S. droebachiensis E in the north Pacific. Within the Atlantic, patterns of habitat segregation suggests that temperature may play a role in limiting the distribution of S. droebachiensis E, particularly throughout the warmer coastal waters along the coast of Nova Scotia. Our results are consistent with the cycles of trans-Arctic dispersal and vicariance in S. pallidus and S. droebachiensis W, but we suggest that the evolution of Atlantic populations of S. droebachiensis E has been driven by persistent trans-Arctic vicariance that may date to the initial invasion in the late Pliocene.

Biological Knowledge of Thornback Ray (Raja clavata) from the Azores: Improving Scientific Information for the Effectiveness of Species-Specific Management Measures

Simple Summary

Sharks, rays, and skates are increasingly being recognized as endangered due to their life-history characteristics, fishing pressure, and habitat degradation. The thornback ray Raja clavata is one of the most commercially important skates in the seas of Northwest Europe. However, due to a lack of biological knowledge about this species in Azorean waters, the types of stock evaluations that can be performed are restricted. This study expands current knowledge on vertical distribution, size-frequency distributions, growth patterns, sex ratios, mortality rates, and reproduction of this species, and provides a baseline for further fishing monitoring.

Abstract

Elasmobranchs are globally recognized as vulnerable due to their life-history characteristics, fishing pressure, and habitat degradation. Among the skates and rays caught by commercial fisheries, the thornback ray Raja clavata is one of the most economically important in Northwest European seas. However, the scarcity of biological knowledge about this species in Azorean waters has limited the stock assessment types that can be conducted. To improve information on its habitat preferences, spatial distribution and movement pattern, growth, sex ratio, mortality, and reproduction, as well as to investigate long-term changes in abundance and size, this study analyzed approximately 25 years of fishery-dependent and independent data from the Azores. Raja clavata was mainly caught at depths up to 250 m. Most of the tagged fish were recaptured near the release point. A larger–deeper trend was found, and females were larger and more abundant than males. Life-history parameters showed that R. clavata has a long lifespan, large size, slow growth, and low natural mortality. The sustainability of its population is of concern to fisheries management and, while our findings suggested a relatively healthy stock in the Azores, a thorough increase in data quality is required to better understand the stock condition and prevent overexploitation.

No leading-edge effect in North Atlantic harbor porpoises: Evolutionary and conservation implications

Understanding species responses to past environmental changes can help forecast how they will cope with ongoing climate changes. Harbor porpoises are widely distributed in the North Atlantic and were deeply impacted by the Pleistocene changes with the split of three subspecies. Despite major impacts of fisheries on natural populations, little is known about population connectivity and dispersal, how they reacted to the Pleistocene changes, and how they will evolve in the future. Here, we used phylogenetics, population genetics, and predictive habitat modeling to investigate population structure and phylogeographic history of the North Atlantic porpoises. A total of 925 porpoises were characterized at 10 microsatellite loci and one quarter of the mitogenome (mtDNA). A highly divergent mtDNA lineage was uncovered in one porpoise off Western Greenland, suggesting that a cryptic group may occur and could belong to a recently discovered mesopelagic ecotype off Greenland. Aside from it and the southern subspecies, spatial genetic variation showed that porpoises from both sides of the North Atlantic form a continuous system belonging to the same subspecies (Phocoena phocoena phocoena). Yet, we identified important departures from random mating and restricted dispersal forming a highly significant isolation by distance (IBD) at both mtDNA and nuclear markers. A ten times stronger IBD at mtDNA compared with nuclear loci supported previous evidence of female philopatry. Together with the lack of spatial trends in genetic diversity, this IBD suggests that migration-drift equilibrium has been reached, erasing any genetic signal of a leading-edge effect that accompanied the predicted recolonization of the northern habitats freed from Pleistocene ice. These results illuminate the processes shaping porpoise population structure and provide a framework for designing conservation strategies and forecasting future population evolution.

Phylogeography in an "oyster" shell provides first insights into the genetic structure of an extinct Ostrea edulis population

The historical phylogeography of Ostrea edulis was successfully depicted in its native range for the first time using ancient DNA methods on dry shells from museum collections. This research reconstructed the historical population structure of the European flat oyster across Europe in the 1870s-including the now extinct population in the Wadden Sea. In total, four haplogroups were identified with one haplogroup having a patchy distribution from the North Sea to the Atlantic coast of France. This irregular distribution could be the result of translocations. The other three haplogroups are restricted to narrow geographic ranges, which may indicate adaptation to local environmental conditions or geographical barriers to gene flow. The phylogenetic reconstruction of the four haplogroups suggests the signatures of glacial refugia and postglacial expansion. The comparison with present-day O. edulis populations revealed a temporally stable population genetic pattern over the past 150 years despite large-scale translocations. This historical phylogeographic reconstruction was able to discover an autochthonous population in the German and Danish Wadden Sea in the late nineteenth century, where O. edulis is extinct today. The genetic distinctiveness of a now-extinct population hints at a connection between the genetic background of O. edulis in the Wadden Sea and for its absence until today.

A globally threatened shark, Carcharias taurus, shows no population decline in South Africa

Knowledge about the demographic histories of natural populations helps to evaluate their conservation status, and potential impacts of natural and anthropogenic pressures. In particular, estimates of effective population size obtained through molecular data can provide useful information to guide management decisions for vulnerable populations. The spotted ragged-tooth shark, Carcharias taurus (also known as the sandtiger or grey nurse shark), is widely distributed in warm-temperate and subtropical waters, but has suffered severe population declines across much of its range as a result of overexploitation. Here, we used multilocus genotype data to investigate the demographic history of the South African C. taurus population. Using approximate Bayesian computation and likelihood-based importance sampling, we found that the population underwent a historical range expansion that may have been linked to climatic changes during the late Pleistocene. There was no evidence for a recent anthropogenic decline. Together with census data suggesting a stable population, these results support the idea that fishing pressure and other threats have so far not been detrimental to the local C. taurus population. The results reported here indicate that South Africa could possibly harbour the last remaining, relatively pristine population of this widespread but vulnerable top predator.

Phylogeographic history of flat periwinkles, Littorina fabalis and L. obtusata

Background

The flat periwinkles, Littorina fabalis and L. obtusata, are two sister species widely distributed throughout the Northern Atlantic shores with high potential to inform us about the process of ecological speciation in the intertidal. However, whether gene flow has occurred during their divergence is still a matter of debate. A comprehensive assessment of the genetic diversity of these species is also lacking and their main glacial refugia and dispersal barriers remain largely unknown. In order to fill these gaps, we sequenced two mitochondrial genes and two nuclear fragments to perform a phylogeographic analysis of flat periwinkles across their distribution range.

Results

We identified two main clades largely composed by species-specific haplotypes corresponding to L. obtusata and L. fabalis, with moderate to strong support, respectively. Importantly, a model of divergence with gene flow between the two species (from L. obtusata to L. fabalis) was better supported, both in Iberia and in northern-central Europe. Three mitochondrial clades were detected within L. fabalis and two within L. obtusata, with strong divergence between Iberia and the remaining populations. The largest component of the genetic variance within each species was explained by differences between geographic regions associated with these clades. Our data suggests that overall intraspecific genetic diversity is similar between the two flat periwinkle species and that populations from Iberia tend to be less diverse than populations from northern-central Europe.

Conclusions

The phylogeographic analysis of this sister-species pair supports divergence with gene flow. This system thus provides us with the opportunity to study the contribution of gene flow and natural selection during diversification. The distribution of the different clades suggests the existence of glacial refugia in Iberia and northern-central Europe for both species, with a main phylogeographic break between these regions. Although the genetic diversity results are not fully conclusive, the lower diversity observed in Iberia could reflect marginal conditions at the southern limit of their distribution range during the current interglacial period.

Phylogeographic and evolutionary history analyses of the warty crab Eriphia verrucosa (Decapoda, Brachyura, Eriphiidae) unveil genetic imprints of a late Pleistocene vicariant event across the Gibraltar Strait, erased by postglacial expansion and admixture among refugial lineages

Background

The Pleistocene cyclic sea-level fluctuations are thought to have markedly affected the distribution and genetic architecture of Atlanto-Mediterranean biota. Despite the acknowledged key role played by these historical events in shaping population genetic structure of marine species, little is still known about the processes involved in shaping the spatial distribution of genetic variation within intertidal species. We intended in this study to reconstruct the phylogeography of a common and widely distributed coastal species across the East Atlantic and Mediterranean Sea (the warty crab Eriphia verrucosa), aiming to unravel potential microevolutionary processes likely involved in shaping its genetic polymorphism. For this purpose, a total of 155 specimens of E. verrucosa from 35 locations across the entire distribution range were analyzed by comparing a 453 basepairs region of the mitochondrial gene cytochrome oxidase subunit 1 (Cox1).

Results

Our results unveiled the prevalence of high genetic connectivity among East Atlantic and Mediterranean populations, with noticeable genetic distinctiveness of the peripheral population from the Azores. Spatio-temporal patterns of genetic diversification and demographic history allowed retrieving genetic imprints of late Pleistocene vicariant event across the Gibraltar Strait followed by subsequent postglacial expansion events for both the East Atlantic and Mediterranean regions. Integrative evidences from the outcomes of comparison of regional genetic diversification, as well as evolutionary and biogeographic histories reconstructions, support the existence of potential glacial refugia for E. verrucosa in the East Atlantic and western Mediterranean. Our results also revealed low levels of genetic variability along with recent demographic and spatial expansion events for eastern Mediterranean warty crabs, suggesting that the eastern areas within the distribution range of the species might have been recently colonized from putative glacial refugia.

Conclusions

These findings provide new insights into the phylogeography and evolutionary history of a common but poorly studied Atlanto-Mediterranean decapod species. Specifically, they contribute to the understanding of the impact of historical processes on shaping contemporary population genetic structure and diversity in intertidal marine species.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1423-2) contains supplementary material, which is available to authorized users.

The hidden side of a major marine biogeographic boundary: a wide mosaic hybrid zone at the Atlantic-Mediterranean divide reveals the complex interaction between natural and genetic barriers in mussels

The Almeria-Oran Front (AOF) is a recognised hotspot of genetic differentiation in the sea, with genetic discontinuities reported in more than 50 species. The AOF is a barrier to dispersal and an ecological boundary; both can determine the position of these genetic breaks. However, the maintenance of genetic differentiation is likely reinforced by genetic barriers. A general drawback of previous studies is an insufficient density of sampling sites at the transition zone, with a conspicuous lack of samples from the southern coastline. We analysed the fine-scale genetic structure in the mussel Mytilus galloprovincialis using a few ancestry-informative loci previously identified from genome scans. We discovered a 600-km-wide mosaic hybrid zone eastward of the AOF along the Algerian coasts. This mosaic zone provides a new twist to our understanding of the Atlantic-Mediterranean transition because it demonstrates that the two lineages can live in sympatry with ample opportunities to interbreed in a large area, but they hardly do so. This implies that some form of reproductive isolation must exist to maintain the two genetic backgrounds locally cohesive. The mosaic zone ends with an abrupt genetic shift at a barrier to dispersal in the Gulf of Bejaia, Eastern Algeria. Simulations of endogenous or exogenous selection in models that account for the geography and hydrodynamic features of the region support the hypothesis that sister hybrid zones could have been differentially trapped at two alternative barriers to dispersal and/or environmental boundaries, at Almeria in the north and Bejaia in the south. A preponderantly unidirectional north-south gene flow next to the AOF can also maintain a patch of intrinsically maintained genetic background in the south and the mosaic structure, even in the absence of local adaptation. Our results concur with the coupling hypothesis that suggests that natural barriers can explain the position of genetic breaks, while their maintenance depends on genetic barriers.

Range-wide genomic data synthesis reveals transatlantic vicariance and secondary contact in Atlantic cod

Recent advances in genetic and genomic analysis have greatly improved our understanding of spatial population structure in marine species. However, studies addressing phylogeographic patterns at oceanic spatial scales remain rare. In Atlantic cod (Gadus morhua), existing range-wide examinations suggest significant transatlantic divergence, although the fine-scale contemporary distribution of populations and potential for secondary contact are largely unresolved. Here, we explore transatlantic phylogeography in Atlantic cod using a data-synthesis approach, integrating multiple genome-wide single-nucleotide polymorphism (SNP) datasets representative of different regions to create a single range-wide dataset containing 1,494 individuals from 54 locations and genotyped at 796 common loci. Our analysis highlights significant transatlantic divergence and supports the hypothesis of westward post-glacial colonization of Greenland from the East Atlantic. Accordingly, our analysis suggests the presence of transatlantic secondary contact off eastern North America and supports existing perspectives on the phylogeographic history of Atlantic cod with an unprecedented combination of genetic and geographic resolution. Moreover, we demonstrate the utility of integrating distinct SNP databases of high comparability.

Meta-analysis of northeast Atlantic marine taxa shows contrasting phylogeographic patterns following post-LGM expansions

Background

Comparative phylogeography enables the study of historical and evolutionary processes that have contributed to shaping patterns of contemporary genetic diversity across co-distributed species. In this study, we explored genetic structure and historical demography in a range of coastal marine species across the northeast Atlantic to assess whether there are commonalities in phylogeographic patterns across taxa and to evaluate whether the timings of population expansions were linked to the Last Glacial Maximum (LGM).

Methods

A literature search was conducted using Web of Science. Search terms were chosen to maximise the inclusion of articles reporting on population structure and phylogeography from the northeast Atlantic; titles and abstracts were screened to identify suitable articles within the scope of this study. Given the proven utility of mtDNA in comparative phylogeography and the availability of these data in the public domain, a meta-analysis was conducted using published mtDNA gene sequences. A standardised methodology was implemented to ensure that the genealogy and demographic history of all mtDNA datasets were reanalysed in a consistent and directly comparable manner.

Results

Mitochondrial DNA datasets were built for 21 species. The meta-analysis revealed significant population differentiation in 16 species and four main types of haplotype network were found, with haplotypes in some species unique to specific geographical locations. A signal of rapid expansion was detected in 16 species, whereas five species showed evidence of a stable population size. Corrected mutation rates indicated that the majority of expansions were estimated to have occurred after the earliest estimate for the LGM (∼26.5 Kyr), while few expansions were estimated to have pre-dated the LGM.

Conclusion

This study suggests that post-LGM expansion appeared to be common in a range of marine taxa, supporting the concept of rapid expansions after the LGM as the ice sheets started to retreat. However, despite the commonality of expansion patterns in many of these taxa, phylogeographic patterns appear to differ in the species included in this study. This suggests that species-specific evolutionary processes, as well as historical events, have likely influenced the distribution of genetic diversity of marine taxa in the northeast Atlantic.

Natural history and molecular evolution of demersal Mediterranean sharks and skates inferred by comparative phylogeographic and demographic analyses

Background

The unique and complex paleoclimatic and paleogeographic events which affected the Mediterranean Sea since late Miocene deeply influenced the distribution and evolution of marine organisms and shaped their genetic structure. Following the Messinian salinity crisis and the sea-level fluctuations during the Pleistocene, several Mediterranean marine species developed deep genetic differentiation, and some underwent rapid radiation. Here, we consider two of the most prioritized groups for conservation in the light of their evolutionary history: sharks and rays (elasmobranchs). This paper deals with a comparative multispecies analysis of phylogeographic structure and historical demography in two pairs of sympatric, phylogenetically- and ecologically-related elasmobranchs, two scyliorhinid catsharks (Galeus melastomus, Scyliorhinus canicula) and two rajid skates (Raja clavata, Raja miraletus). Sampling and experimental analyses were designed to primarily test if the Sicilian Channel can be considered as effective eco-physiological barrier for Mediterranean demersal sympatric elasmobranchs.

Methods

The phylogeography and the historical demography of target species were inferred by analysing the nucleotide variation of three mitochondrial DNA markers (i.e., partial sequence of COI, NADH2 and CR) obtained from a total of 248 individuals sampled in the Western and Eastern Mediterranean Sea as well as in the adjacent northeastern Atlantic Ocean. Phylogeographic analysis was performed by haplotype networking and testing spatial genetic differentiation of samples (i.e., analysis of molecular variance and of principal components). Demographic history of Mediterranean populations was reconstructed using mismatch distribution and Bayesian Skyline Plot analyses.

Results

No spatial genetic differentiation was identified in either catshark species, while phylogeographic structure of lineages was identified in both skates, with R. miraletus more structured than R. clavata. However, such structuring of skate lineages was not consistent with the separation between Western and Eastern Mediterranean. Sudden demographic expansions occurred synchronously during the upper Pleistocene (40,000–60,000 years ago) in both skates and G. melastomus, likely related to optimal environmental conditions. In contrast, S. canicula experienced a slow and constant increase in population size over the last 350,000 years.

Discussion

The comparative analysis of phylogeographic and historical demographic patterns for the Mediterranean populations of these elasmobranchs reveals that historical phylogeographic breaks have not had a large impact on their microevolution. We hypothesize that interactions between environmental and ecological/physiological traits may have been the driving force in the microevolution of these demersal elasmobranch species in the Mediterranean rather than oceanographic barriers.

Parapatric genetic divergence among deep evolutionary lineages in the Mediterranean green crab, Carcinus aestuarii (Brachyura, Portunoidea, Carcinidae), accounts for a sharp phylogeographic break in the Eastern Mediterranean

BACKGROUND

Recently, population genetic studies of Mediterranean marine species highlighted patterns of genetic divergence and phylogeographic breaks, due to the interplay between impacts of Pleistocene climate shifts and contemporary hydrographical barriers. These factors markedly shaped the distribution of marine organisms and their genetic makeup. The present study is part of an ongoing effort to understand the phylogeography and evolutionary history of the highly dispersive Mediterranean green crab, Carcinus aestuarii (Nardo, 1847), across the Mediterranean Sea. Recently, marked divergence between two highly separated haplogroups (genetic types I and II) of C. aestuarii was discerned across the Siculo-Tunisian Strait, suggesting an Early Pleistocene vicariant event. In order to better identify phylogeographic patterns in this species, a total of 263 individuals from 22 Mediterranean locations were analysed by comparing a 587 basepair region of the mitochondrial gene Cox1 (cytochrome oxidase subunit 1). The examined dataset is composed of both newly generated sequences (76) and previously investigated ones (187).

RESULTS

Our results unveiled the occurrence of a highly divergent haplogroup (genetic type III) in the most north-eastern part of the Mediterranean Sea. Divergence between the most distinct type III and the common ancestor of both types I and II corresponds to the Early Pleistocene and coincides with the historical episode of separation between types I and II. Our results also revealed strong genetic divergence among adjacent regions (separating the Aegean and Marmara seas from the remaining distribution zone) and confirmed a sharp phylogeographic break across the Eastern Mediterranean. The recorded parapatric genetic divergence, with the potential existence of a contact zone between both groups in the Ionian Sea and notable differences in the demographic history, suggest the likely impact of paleoclimatic events, as well as past and contemporary oceanographic processes, in shaping genetic variability of this species.

CONCLUSIONS

Our findings not only provide further evidence for the complex evolutionary history of the green crab in the Mediterranean Sea, but also stress the importance of investigating peripheral areas in the species' distribution zone in order to fully understand the distribution of genetic diversity and unravel hidden genetic units and local patterns of endemism.

Insights into the genetic structure of the rabbitfish Chimaera monstrosa (Holocephali) across the Atlantic-Mediterranean transition zone

A comparison of the genetic structure of Chimaera monstrosa populations from the Atlantic Ocean and the Mediterranean Sea was carried out using mitochondrial DNA analysis. Results indicate high and significant pairwise Φ_ST values with no shared haplotypes between the two areas. Furthermore, migration rate estimates suggested absence of gene flow between the two basins. These findings, coupled with the species vertical distribution, suggest that the shallow depth of the Strait of Gibraltar may act as a barrier limiting the dispersal capabilities of these populations, which likely became separated at the end of the middle Pleistocene.

Population genetics of the deep-sea bluntnose sixgill shark, Hexanchus griseus, revealing spatial genetic heterogeneity

Hexanchus griseus is a globally distributed deep-water shark species. It inhabits tropical and temperate waters throughout the world, including the Mediterranean Sea where it is by-caught by small-scale fisheries in the region. In this study, we analysed the genetic variation of H. griseus specimens collected from different areas within and outside the Mediterranean region, to assess its genetic connectivity. The mitochondrial DNA (mtDNA) sequence analysed in this study ranged from cytochrome b to 16S rRNA genes including the control region, the 12S rRNA gene and the interspersed tRNA genes in the region, covering a total of 3731 to 3914 nucleotides. Results have shown that this species exhibits geographically distinct maternal lineages, indicating population structure along geographical ranges. These findings reveal population subdivisions not only between the Pacific Ocean and the Atlantic Ocean, but also within the oceans and on a smaller scale within the Mediterranean Sea. This highlights the need to consider each population subdivision separately when designing management plans for the conservation of this species.

Historical gene flow constraints in a northeastern Atlantic fish: phylogeography of the ballan wrasse Labrus bergylta across its distribution range

The distribution and demographic patterns of marine organisms in the north Atlantic were largely shaped by climatic changes during the Pleistocene, when recurrent glacial maxima forced them to move south or to survive in northern peri-glacial refugia. These patterns were also influenced by biological and ecological factors intrinsic to each species, namely their dispersion ability. The ballan wrasse (Labrus bergylta), the largest labrid fish along Europe's continental margins, is a target for fisheries and aquaculture industry. The phylogeographic pattern, population structure, potential glacial refugia and recolonization routes for this species were assessed across its full distribution range, using mitochondrial and nuclear markers. The existence of a marked population structure can reflect both recolonization from three distinct glacial refugia and current and past oceanographic circulation patterns. Although isolated in present times, shared haplotypes between continental and Azores populations and historical exchange of migrants in both directions point to a common origin of L. bergylta. This situation is likely to be maintained and/or accentuated by current circulation patterns in the north Atlantic, and may lead to incipient speciation in the already distinct Azorean population. Future monitoring of this species is crucial to evaluate how this species is coping with current environmental changes.

Improving the Conservation of Mediterranean Chondrichthyans: The ELASMOMED DNA Barcode Reference Library

Cartilaginous fish are particularly vulnerable to anthropogenic stressors and environmental change because of their K-selected reproductive strategy. Accurate data from scientific surveys and landings are essential to assess conservation status and to develop robust protection and management plans. Currently available data are often incomplete or incorrect as a result of inaccurate species identifications, due to a high level of morphological stasis, especially among closely related taxa. Moreover, several diagnostic characters clearly visible in adult specimens are less evident in juveniles. Here we present results generated by the ELASMOMED Consortium, a regional network aiming to sample and DNA-barcode the Mediterranean Chondrichthyans with the ultimate goal to provide a comprehensive DNA barcode reference library. This library will support and improve the molecular taxonomy of this group and the effectiveness of management and conservation measures. We successfully barcoded 882 individuals belonging to 42 species (17 sharks, 24 batoids and one chimaera), including four endemic and several threatened ones. Morphological misidentifications were found across most orders, further confirming the need for a comprehensive DNA barcoding library as a valuable tool for the reliable identification of specimens in support of taxonomist who are reviewing current identification keys. Despite low intraspecific variation among their barcode sequences and reduced samples size, five species showed preliminary evidence of phylogeographic structure. Overall, the ELASMOMED initiative further emphasizes the key role accurate DNA barcoding libraries play in establishing reliable diagnostic species specific features in otherwise taxonomically problematic groups for biodiversity management and conservation actions.

Resolving taxonomic uncertainty in vulnerable elasmobranchs: are the Madeira skate (Raja maderensis) and the thornback ray (Raja clavata) distinct species?

Skates and rays constitute the most speciose group of chondrichthyan fishes, yet are characterised by remarkable levels of morphological and ecological conservatism. They can be challenging to identify, which makes monitoring species compositions for fisheries management purposes problematic. Owing to their slow growth and low fecundity, skates are vulnerable to exploitation and species exhibiting endemism or limited ranges are considered to be the most at risk. The Madeira skate Raja maderensis is endemic and classified as ‘Data Deficient’ by the IUCN, yet its taxonomic distinctiveness from the morphologically similar and more wide-ranging thornback ray Raja clavata is unresolved. This study evaluated the sequence divergence of both the variable control region and cytochrome oxidase I ‘DNA barcode’ gene of the mitochondrial genome to elucidate the genetic differentiation of specimens identified as R. maderensis and R. clavata collected across much of their geographic ranges. Genetic evidence was insufficient to support the different species designations. However regardless of putative species identification, individuals occupying waters around the Azores and North African Seamounts represent an evolutionarily significant unit worthy of special consideration for conservation management.

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep-sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep-sea ecosystems is scarce, especially for deep-sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep-sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497-bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ(ST) = -0.0091, F(ST) = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ(ST) = 0.0744, P < 0.01; F(ST) = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ(ST) = 0.3808, F(ST) = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.

Strong population genetic structure and contrasting demographic histories for the small-spotted catshark (Scyliorhinus canicula) in the Mediterranean Sea

Coastal and demersal chondrichthyans, such as the small-spotted catshark, are expected to exhibit genetic differentiation in areas of complex geomorphology like the Mediterranean Basin because of their limited dispersal ability. To test this hypothesis, we used a fragment of the mitochondrial cytochrome c oxidase subunit I gene and 12 nuclear microsatellite loci in order to investigate the genetic structure and historical demography of this species, and to identify potential barriers to gene flow. Samples were collected from the Balearic Islands, the Algerian Basin, the Ionian Sea, the Corinthian Gulf and various locations across the Aegean Sea. Additional sequences from the Atlantic and the Levantine Basin retrieved from GenBank were included in the mitochondrial DNA analysis. Both mitochondrial and nuclear microsatellite DNA data revealed a strong genetic subdivision, mainly between the western and eastern Mediterranean, whereas the Levantine Basin shared haplotypes with both areas. The geographic isolation of the Mediterranean basins seems to enforce the population genetic differentiation of the species, with the deep sea acting as a strong barrier to its dispersal. Contrasting historical demographic patterns were also observed in different parts of the species' distribution, most notably a population growth trend in the western Mediterranean/Atlantic area and a slight decreasing one in the Aegean Sea. The different effects of the Pleistocene glacial periods on the habitat availability may explain the contrasting demographic patterns observed. The current findings suggest that the small-spotted catshark exhibits several genetic stocks in the Mediterranean, although further study is needed.

Contemporary population structure and post-glacial genetic demography in a migratory marine species, the blacknose shark, Carcharhinus acronotus

Patterns of population structure and historical genetic demography of blacknose sharks in the western North Atlantic Ocean were assessed using variation in nuclear-encoded microsatellites and sequences of mitochondrial (mt)DNA. Significant heterogeneity and/or inferred barriers to gene flow, based on microsatellites and/or mtDNA, revealed the occurrence of five genetic populations localized to five geographic regions: the southeastern U.S Atlantic coast, the eastern Gulf of Mexico, the western Gulf of Mexico, Bay of Campeche in the southern Gulf of Mexico and the Bahamas. Pairwise estimates of genetic divergence between sharks in the Bahamas and those in all other localities were more than an order of magnitude higher than between pairwise comparisons involving the other localities. Demographic modelling indicated that sharks in all five regions diverged after the last glacial maximum and, except for the Bahamas, experienced post-glacial, population expansion. The patterns of genetic variation also suggest that the southern Gulf of Mexico may have served as a glacial refuge and source for the expansion. Results of the study demonstrate that barriers to gene flow and historical genetic demography contributed to contemporary patterns of population structure in a coastal migratory species living in an otherwise continuous marine habitat. The results also indicate that for many marine species, failure to properly characterize barriers in terms of levels of contemporary gene flow could in part be due to inferences based solely on equilibrium assumptions. This could lead to erroneous conclusions regarding levels of connectivity in species of conservation concern.

A tale of two seas: contrasting patterns of population structure in the small-spotted catshark across Europe

Elasmobranchs represent important components of marine ecosystems, but they can be vulnerable to overexploitation. This has driven investigations into the population genetic structure of large-bodied pelagic sharks, but relatively little is known of population structure in smaller demersal taxa, which are perhaps more representative of the biodiversity of the group. This study explores spatial population genetic structure of the small-spotted catshark (Scyliorhinus canicula), across European seas. The results show significant genetic differences among most of the Mediterranean sample collections, but no significant structure among Atlantic shelf areas. The data suggest the Mediterranean populations are likely to have persisted in a stable and structured environment during Pleistocene sea-level changes. Conversely, the Northeast Atlantic populations would have experienced major changes in habitat availability during glacial cycles, driving patterns of population reduction and expansion. The data also provide evidence of male-biased dispersal and female philopatry over large spatial scales, implying complex sex-determined differences in the behaviour of elasmobranchs. On the basis of this evidence, we suggest that patterns of connectivity are determined by trends of past habitat stability that provides opportunity for local adaptation in species exhibiting philopatric behaviour, implying that resilience of populations to fisheries and other stressors may differ across the range of species.

Reproductive isolation and patterns of genetic differentiation in a cryptic butterfly species complex

Molecular studies of natural populations are often designed to detect and categorize hidden layers of cryptic diversity, and an emerging pattern suggests that cryptic species are more common and more widely distributed than previously thought. However, these studies are often decoupled from ecological and behavioural studies of species divergence. Thus, the mechanisms by which the cryptic diversity is distributed and maintained across large spatial scales are often unknown. In 1988, it was discovered that the common Eurasian Wood White butterfly consisted of two species (Leptidea sinapis and Leptidea reali), and the pair became an emerging model for the study of speciation and chromosomal evolution. In 2011, the existence of a third cryptic species (Leptidea juvernica) was proposed. This unexpected discovery raises questions about the mechanisms preventing gene flow and about the potential existence of additional species hidden in the complex. Here, we compare patterns of genetic divergence across western Eurasia in an extensive data set of mitochondrial and nuclear DNA sequences with behavioural data on inter- and intraspecific reproductive isolation in courtship experiments. We show that three species exist in accordance with both the phylogenetic and biological species concepts and that additional hidden diversity is unlikely to occur in Europe. The Leptidea species are now the best studied cryptic complex of butterflies in Europe and a promising model system for understanding the formation of cryptic species and the roles of local processes, colonization patterns and heterospecific interactions for ecological and evolutionary divergence.

Phylogeography of bivalve Cyclina sinensis: testing the historical glaciations and Changjiang River outflow hypotheses in northwestern Pacific

Background

The marginal seas of northwestern Pacific are characterized by unique topography and intricate hydrology. Two hypotheses have been proposed to explain genetic patterns of marine species inhabiting the region: the historical glaciations hypothesis suggests population genetic divergence between sea basins, whereas the Changjiang River outflow hypothesis suggests genetic break in line with the Changjiang Estuary. Here the phylogeography of bivalve Cyclina sinensis was investigated to test the validity of these two hypotheses for intertidal species in three marginal seas—the East China Sea (ECS), the South China Sea (SCS), and the Japan Sea (JPS).

Methodology/Principal Findings

Fragments of two markers (mitochondrial COI and nuclear ITS-1) were sequenced for 335 individuals collected from 21 populations. Significant pairwise Φ_ST were observed between different marginal sea populations. Network analyses illustrated restricted distribution of haplogroups/sub-haplogroups to sea basins, with a narrow secondary contact zone between the ECS and SCS. Demographic expansion was inferred for ECS and SCS lineages using mismatch distributions, neutral tests, and extended Bayesian Skyline Plots. Based on a molecular clock method, the divergence times among COI lineages were estimated dating from the Pleistocene.

Conclusions

The phylogeographical break revealed for C. sinensis populations is congruent with the historical isolation of sea basins rather than the putative Changjiang River outflow barrier. The large land bridges extending between seas during glaciation allowed accumulation of mutations and subsequently gave rise to deep divergent lineages. The low-dispersal capacity of the clam and coastal oceanography may facilitate the maintenance of the historical patterns as barriers shift. Our study supports the historical glaciations hypothesis for interpreting present-day phylogeographical patterns of C. sinensis, and highlights the importance of historical glaciations for generating genetic structure of marine coastal species especially those with low-dispersal abilities in northwestern Pacific.

Unexpected high genetic diversity at the extreme northern geographic limit of Taurulus bubalis (Euphrasen, 1786)

The longspined bullhead (Taurulus bubalis, Euphrasen 1786) belongs to the family Cottidae and is a rocky shore species that inhabits the intertidal zones of the Eastern Atlantic since Iceland, southward to Portugal and also the North Sea and Baltic, northward to the Gulf of Finland, with some occurrences in the northern Mediterranean coasts eastward to the Gulf of Genoa. We analysed the phylogeographic patterns of this species using mitochondrial and nuclear markers in populations throughout most of its distributional range in west Europe. We found that T. bubalis has a relatively shallow genealogy with some differentiation between Atlantic and North Sea. Genetic diversity was homogeneous across all populations studied. The possibility of a glacial refugium near the North Sea is discussed. In many, but not all, marine temperate organisms, patterns of diversity are similar across the species range. If this phenomenon proves to be most common in cold adapted species, it may reflect the availability of glacial refugia not far from their present-day northern limits.

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.

Molecular markers: progress and prospects for understanding reproductive ecology in elasmobranchs

Application of modern molecular tools is expanding the understanding of elasmobranch reproductive ecology. High-resolution molecular markers provide information at scales ranging from the identification of reproductively isolated populations in sympatry (i.e. cryptic species) to the relationships among parents, offspring and siblings. This avenue of study has not only augmented the current understanding of the reproductive biology of elasmobranchs but has also provided novel insights that could not be obtained through experimental or observational techniques. Sharing of genetic polymorphisms across ocean basins indicates that for some species there may be gene flow on global scales. The presence, however, of morphologically similar but genetically distinct entities in sympatry suggests that reproductive isolation can occur with minimal morphological differentiation. This review discusses the recent findings in elasmobranch reproductive biology like philopatry, hybridization and polyandry while highlighting important molecular and analytical techniques. Furthermore, the review examines gaps in current knowledge and discusses how new technologies may be applied to further the understanding of elasmobranch reproductive ecology.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.

Population structure and phylogeography of the short-tailed stingray, Dasyatis brevicaudata (Hutton 1875), in the Southern Hemisphere

There is accumulating evidence that the degree of vagility explains little of the extent of population subdivision found within elasmobranch species. Instead, patterns of gene flow in elasmobranchs appear more closely linked to the presence of dispersal barriers, either physical or biological. Here, we investigate the potential role of some of these isolating mechanisms in shaping the population structure of a member of the stingray family Dasyatidae (Dasyatis brevicaudata) at various scales (southern hemisphere vs. coastal New Zealand). Analyses of the mitochondrial DNA control region from 176 individuals revealed significant genetic structure between South Africa, Australia, and New Zealand populations (analysis of molecular variance [AMOVA], overall Ф(ST) = 0.67, P < 0.001), although New Zealand and Australia shared some haplotypes. Surprisingly, significant population differentiation was found among several coastal New Zealand locations (AMOVA, overall Ф(ST) = 0.05, P < 0.05). However, data did not support the genetic differentiation between individuals from an offshore breeding area and mainland individuals. Comparisons suggest that these stingrays exhibit similar levels of population differentiation as other coastal elasmobranchs, with high divergence across oceanic basins and lower differentiation along continuous coastal habitats. Differences in coastal population structuring in elasmobranch species studied to date may be attributed to species-specific preferences for coastal habitats, which may be linked to life history functions (e.g., feeding and pupping).