Mitochondrial DNA reveals a mosaic pattern of phylogeographical structure in Atlantic and Mediterranean populations of anchovy (Engraulis encrasicolus)

Santamaría MT, Faria R, Machado AM, Fonseca MM, Froufe E, C. Castro LF. Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Mediterranean fish communities are struggling to adapt to global warming. Evidence from the western coast of Italy

Climate change has significant impacts on marine ecosystems, resulting in disruptions in biological interactions, shifts in community composition, and changes in the physiology of fish and other marine organisms. In this study conducted in the central Mediterranean Sea, the mean temperature of the catch (MTC) was employed as an indicator to investigate the climatological factors influencing the fish community. The MTC, which utilizes species-preferred temperatures, was calculated using bottom temperature (BT) data weighted against scientific catches. The estimated MTC increasing rates were 0.01 °C year-1 for the entire community, 0.017 °C year-1 for the shelf break, and 0.004 °C year-1 for the continental slope assemblage. We found that MTC is increasing at a lower rate compared to BT, suggesting a progressive under-adaptation of the fish community that seems not fully able to keep up with the ongoing pace of warming. The study identified sea surface temperature and bottom temperature as key drivers of changes in fish community composition. Notably, the fish community composition exhibited drastic changes over the studied period, and we suggest that the MTC can be a useful index to monitor such changes within the context of the EU's climate change adaptation strategy.

Mitonuclear genetic patterns of divergence in the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) along the Turkish seas

Biogeographical transition zones present good opportunities for studying the effect of the past ice ages on genetic structure of species because secondary contact zones of post-glacial lineages can be formed. In this study, we investigated the population genetic structure of the marbled rock crab, Pachygrapsus marmoratus along the coasts of Turkey. We genotyped 334 individuals from the Black Sea, the Turkish Straits System (TSS), the Aegean, and the Eastern Mediterranean basins. In order to reveal its evolutionary history and its population connectivity, we used mitochondrial CO1 region and five microsatellite loci. CO1 analyzes also included 610 additional samples from Genbank, which covered most of its distribution range. Both microsatellites and mtDNA showed decreased diversity in sampling sites of the TSS and the Black Sea as compared to those along the Aegean and the Levantine coasts. There is an especially strong geographical pattern in distributions of haplotypes in mtDNA, most probably as a result of genetic drift in the Black Sea and the Sea of Marmara (SoM). Microsatellite data analyses revealed two genetically distinct clusters of P. marmoratus (clusters C and M). While individuals belonging to cluster C are present in all the sampling locations, those belonging to cluster M are only detected along the Mediterranean coasts including the Aegean and the Levantine basins. These clusters shared similar haplotypes in the Mediterranean. Haplotypes of two sympatric clusters could be similar due to incomplete lineage sorting of ancestral polymorphisms. In order to retrieve the complex demographic history and to investigate evolutionary processes resulting in sympatric clusters in the Aegean Sea and the Levantine basin, mitochondrial markers with faster mutation rates than CO1 and/or SNP data will be useful.

Genetic hypervariability of a Northeastern Atlantic venomous rockfish

BACKGROUND

Understanding the interplay between climate and current and historical factors shaping genetic diversity is pivotal to infer changes in marine species range and communities' composition. A phylogeographical break between the Atlantic and the Mediterranean has been documented for several marine organisms, translating into limited dispersal between the two basins.

METHODS

In this study, we screened the intraspecific diversity of 150 individuals of the Madeira rockfish (Scorpaena maderensis) across its distributional range (seven sampling locations in the Atlantic and Mediterranean basins) using the mitochondrial control region and the nuclear S7 first intron.

RESULTS

The present work is the most comprehensive study done for this species, yielding no genetic structure across sampled locations and no detectable Atlantic-Mediterranean break in connectivity. Our results reveal deep and hyper-diverse bush-like genealogies with large numbers of singletons and very few shared haplotypes. The genetic hyper-diversity found for the Madeira rockfish is relatively uncommon in rocky coastal species, whose dispersal capability is limited by local oceanographic patterns. The effect of climate warming on the distribution of the species is discussed.

Signals from the deep-sea: Genetic structure, morphometric analysis, and ecological implications of Cyclothone braueri (Pisces, Gonostomatidae) early life stages in the Central Mediterranean Sea

Cyclothone braueri (Stomiiformes, Gonostomatidae) is a widely distributed fish inhabiting the mesopelagic zone of marine tropical and temperate waters. Constituting one of the largest biomasses of the ocean, C. braueri is a key element in most of the ecological processes occurring in the twilight layer. We focused on the ecological processes linked to early life stages in relation to marine pelagic environmental drivers (temperature, salinity, food availability and geostrophic currents) considering different regions of the Central Mediterranean Sea. A multivariate morphometric analysis was carried out using six parameters with the aim of discerning different larval morphotypes, while a fragment of 367 bp representing the 12S ribosomal RNA gene was used to perform molecular analyses aimed at determining the intraspecific genetic variability. Analysis highlighted two geographically distinct morphotypes not genetically discernible and related to the different nutritional conditions due to spatial heterogeneities in terms of temperature and food availability. The body depth (BD) emerged as an appropriate morphometric parameter to detect the larval condition in this species. Molecular analysis highlighted a moderate genetic divergence in the fish population, showing the recurrence of two phylogroups not geographically separated.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Population Genetic Diversity of Two Marine Gobies (Gobiiformes: Gobiidae) from the North-Eastern Atlantic and the Mediterranean Sea

Gobies (Gobiiformes: Gobiidae) are the most species-rich family of fishes in general, and the most abundant fish group in the European seas. Nonetheless, our knowledge on many aspects of their biology, including the population genetic diversity, is poor. Although barriers to gene flow are less apparent in the marine environment, the ocean is not a continuous habitat, as has been shown by studies on population genetics of various marine biota. For the first time, European marine goby species which cannot be collected by common fishery techniques were studied. The population genetic structure of two epibenthic species, Gobius geniporus and Gobius cruentatus, from seven localities across their distribution ranges was assessed, using one mitochondrial (cytochrome b) and one nuclear gene (first intron of ribosomal protein gene S7). Our results showed that there is a great diversity of haplotypes of mitochondrial gene cytochrome b in both species at all localities. Global fixation indices (FST) indicated a great differentiation of populations in both studied gobies. Our results did not show a geographic subdivision to individual populations. Instead, the data correspond with the model of migration which allow divergence and recurrent migration from the ancestral population. The estimated migration routes coincide with the main currents in the studied area. This matches well the biology of the studied species, with adults exhibiting only short-distance movements and planktonic larval stages.

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.

Genetic divergence in natural and farm populations of Pengba fish, Osteobrama belangeri (Valenciennes, 1844), an endemic fish of North-East India derived from mtDNA ATPase 6/8 gene

Pengba fish, Osteobrama belangeri is a freshwater inhabitant, highly endemic, threatened and economically important minor carp for its food value. In the present investigation, population structure of O. belangeri was examined using mitochondrial ATPase 6/8 gene from geographically distinct locations along the Indo-Burma biodiversity hotspot. A total of 80 individuals were collected belonging to natural and farm populations. The hierarchical analysis of molecular variance and conventional Fst values (0.825 in ATPase 6/8, p < .05) indicated significant genetic structure among populations. The result showed that ATPase 6/8 genes are potential marker in determining the genetic divergence between natural and farm populations of O. belangeri from North-East India.

Spatial genetic structure and body size divergence in endangered Gymnogobius isaza in ancient Lake Biwa

Gymnogobius isaza is a freshwater goby endemic to ancient Lake Biwa, the largest lake in Japan. The species is now listed as 'Critically Endangered' in the Red Data Book of Japan. Nevertheless, it remains subject to fishing without any specific management strategies. Previous studies using mitochondrial DNA markers showed that this fish species has two cryptic lineages. However, little is known about spatial genetic structure and ecological differences across the broad lakescape. In this study, we collected fish samples at nine locations along the lakeshore during the breeding season and tested for the presence of spatial heterogeneity in the lineage's composition while measuring body size as the most fundamental biological trait. The results showed that the major lineage dominated all the sampling locations whereas the minor lineage consisted of only 11% (16/143) of samples. Furthermore, although their spatial distributions overlapped (i.e. the two lineages may be well mixed), we found it possible that the minor lineage may have a potentially narrower distribution than the major lineage. In addition, we found that the two lineages differ in body size; specifically, the minor lineage is smaller in size. From the viewpoint of genetic diversity conservation and sustainable resource use, this fish should be managed as two genetic stocks and spatial and/or body size-based fishery management is desirable, with particular attention to the minor (smaller sized) lineage.

Insights on the drivers of genetic divergence in the European anchovy

Anchovies represent the largest world's marine fish catches and the current threats on their populations impose a sustainable exploitment based on sound scientific information. In the European anchovy (Engraulis encrasicolus), the existence of several populations has been proposed but a global view is missing. Using a multidisciplinary approach, here we assessed the divergence among different ecotypes and its possible causes. SNPs have revealed two functionally distinct ecotypes overlapping in the Central Mediterranean, with one ecotype confined near the river estuaries. The same SNPs outliers also segregated two distinct populations in the near Atlantic, despite their large spatial distance. In addition, while most studies suggested that adaptation to low salinity is key to divergence, here we show that the offshore ecotype has higher environmental tolerance and an opportunistic feeding behaviour, as assessed by the study of environmental conditions, anchovy diet and trophic levels, and passive egg dispersal. These results provide insights into the anchovy evolutionary history, stressing the importance of behaviour in shaping ecotypes.

Wandering behaviour prevents inter and intra oceanic speciation in a coastal pelagic fish

Small pelagic fishes have the ability to disperse over long distances and may present complex evolutionary histories. Here, Old World Anchovies (OWA) were used as a model system to understand genetic patterns and connectivity of fish between the Atlantic and Pacific basins. We surveyed 16 locations worldwide using mtDNA and 8 microsatellite loci for genetic parameters, and mtDNA (cyt b; 16S) and nuclear (RAG1; RAG2) regions for dating major lineage-splitting events within Engraulidae family. The OWA genetic divergences (0-0.4%) are compatible with intra-specific divergence, showing evidence of both ancient and contemporary admixture between the Pacific and Atlantic populations, enhanced by high asymmetrical migration from the Pacific to the Atlantic. The estimated divergence between Atlantic and Pacific anchovies (0.67 [0.53-0.80] Ma) matches a severe drop of sea temperature during the Günz glacial stage of the Pleistocene. Our results support an alternative evolutionary scenario for the OWA, suggesting a coastal migration along south Asia, Middle East and eastern Africa continental platforms, followed by the colonization of the Atlantic via the Cape of the Good Hope.

Contrasted levels of genetic diversity in a benthic Mediterranean octocoral: Consequences of different demographic histories?

Understanding the factors explaining the observed patterns of genetic diversity is an important question in evolutionary biology. We provide the first data on the genetic structure of a Mediterranean octocoral, the yellow gorgonian Eunicella cavolini, along with insights into the demographic history of this species. We sampled populations in four areas of the Mediterranean Sea: continental France, Algeria, Turkey, and the Balearic and Corsica islands. Along French coasts, three sites were sampled at two depths (20 and 40 m). We demonstrated a high genetic structure in this species (overall F_ST = 0.13), and most pairwise differentiation tests were significant. We did not detect any difference between depths at the same site. Clustering analyses revealed four differentiated groups corresponding to the main geographical areas. The levels of allelic richness and heterozygosity were significantly different between regions, with highest diversity in Algeria and lowest levels in Turkey. The highest levels of private allelic richness were observed in Algeria followed by Turkey. Such contrasted patterns of genetic diversity were not observed in other Mediterranean octocorals and could be the result of different evolutionary histories. We also provide new empirical evidence of contrasting results between tests and model‐based studies of demographic history. Our results have important consequences for the management of this species.

Parallel genetic divergence among coastal-marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome-wide data set of RAD-derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal-marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine-scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.

Contemporary genetic structure and postglacial demographic history of the black scorpionfish, Scorpaena porcus, in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.

The European anchovy, a genetically highly diverse species displays null within-sample haplotype diversity on a single study?

The European anchovy has been the focus of numerous population genetic studies, most of which exposing high levels of haplotype diversity. However, Keskin and Atar (2012) revealed rather singular results of null haplotype diversities. We therefore call for caution when considering these findings.

Is there a distinct harbor porpoise subpopulation in the Marmara Sea?

Genetic population structure of geographically isolated endangered Black Sea harbor porpoise (Phocoena phocoena relicta) is little known in Turkish waters, especially in the Turkish Straits System (TSS- Marmara Sea, Bosphorus and Dardanelles), which connects the Black Sea and the Aegean Sea. Mitochondrial DNA sequences of 70 new individuals sampled in the Turkish Black Sea, TSS and Aegean Sea, revealed five new haplotypes from the Black Sea. The findings support the idea that harbor porpoises from the Black Sea dispersed into the Aegean through the TSS. Considering signatures of population expansion, all subpopulations showed a signature of population expansion. The network data and the Фst calculations indicated that the Marmara Sea subpopulation was significantly differentiated from all of the other subpopulations, and supports the notion of its isolated. The finding of a potential management unit (MU) within an already heavily impacted subpopulation as a whole suggests that the individuals of P. p. relicta inhabiting the Marmara Sea require a very rigorous conservation strategy to ensure the survival of this subpopulation, represented by its unique haplotype.

A COI Nonsynonymous Mutation as Diagnostic Tool for Intraspecific Discrimination in the European Anchovy Engraulis encrasicolus (Linnaeus)

The European anchovy, Engraulis encrasicolus, is currently one of the principal target species for commercial fisheries in Europe. In this study, the mitochondrial Control Region (CR) and the Cytochrome Oxidase I (COI) mitochondrial gene were analyzed in 74 specimens of E. encrasicolus from four localities in the central Mediterranean. In both populations, the two markers revealed the presence of two main haplogroups, A and B, already detected in previous investigations of different classes of molecular markers. Both CR and COI markers consistently identified two haplogroups. The COI sequence analysis identified a non-synonymous transversion (T to G) at position 116 of the translated sequence, resulting in an amino acid change. All COI sequences of haplogroup A had an amino acid sequence with alanine in this position, while serine was present in the same position in haplogroup B. The two haplogroups A and B were also discriminated by the variable number of TACA elements at the 5’-end of the mitochondrial CR. The selection tests applied to the COI dataset revealed that codon 116 was not under positive selection, that seven amino acid changes were under purifying selection, and that two amino acids were under episodic positive selection.

Thermal adaptation and clinal mitochondrial DNA variation of European anchovy

Natural populations of widely distributed organisms often exhibit genetic clinal variation over their geographical ranges. The European anchovy, Engraulis encrasicolus, illustrates this by displaying a two-clade mitochondrial structure clinally arranged along the eastern Atlantic. One clade has low frequencies at higher latitudes, whereas the other has an anti-tropical distribution, with frequencies decreasing towards the tropics. The distribution pattern of these clades has been explained as a consequence of secondary contact after an ancient geographical isolation. However, it is not unlikely that selection acts on mitochondria whose genes are involved in relevant oxidative phosphorylation processes. In this study, we performed selection tests on a fragment of 1044 bp of the mitochondrial cytochrome b gene using 455 individuals from 18 locations. We also tested correlations of six environmental features: temperature, salinity, apparent oxygen utilization and nutrient concentrations of phosphate, nitrate and silicate, on a compilation of mitochondrial clade frequencies from 66 sampling sites comprising 2776 specimens from previously published studies. Positive selection in a single codon was detected predominantly (99%) in the anti-tropical clade and temperature was the most relevant environmental predictor, contributing with 59% of the variance in the geographical distribution of clade frequencies. These findings strongly suggest that temperature is shaping the contemporary distribution of mitochondrial DNA clade frequencies in the European anchovy.

Population genetic structure of chub mackerel Scomber japonicus in the Northwestern Pacific inferred from microsatellite analysis

Marine pelagic fishes are usually characterized by subtle but complex patterns of genetic differentiation, which are influenced by both historical process and contemporary gene flow. Genetic population differentiation of chub mackerel, Scomber japonicus, was examined across most of its range in the Northwestern Pacific by screening variation of eight microsatellite loci. Our genetic analysis detected a weak but significant genetic structure of chub mackerel, which was characterized by areas of gene flow and isolation by distance. Consistent with previous estimates of stock structure, we found genetic discontinuity between Japan and China samples. Local-scale pattern of genetic differentiation was observed between samples from the Bohai Sea and North Yellow Sea and those from the East China Sea, which we ascribed to differences in spawning time and migratory behavior. Furthermore, the observed homogeneity among collections of chub mackerel from the East and South China Seas could be the result of an interaction between biological characteristics and marine currents. The present study underlies the importance of understanding the biological significance of genetic differentiation to establish management strategies for exploited fish populations.

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.

SNP discovery in European anchovy (Engraulis encrasicolus, L) by high-throughput transcriptome and genome sequencing

Increased throughput in sequencing technologies has facilitated the acquisition of detailed genomic information in non-model species. The focus of this research was to discover and validate SNPs derived from the European anchovy (Engraulis encrasicolus) transcriptome, a species with no available reference genome, using next generation sequencing technologies. A cDNA library was constructed from four tissues of ten fish individuals corresponding to three populations of E. encrasicolus, and Roche 454 GS FLX Titanium sequencing yielded 19,367 contigs. Additionally, the European anchovy genome was sequenced for the same ten individuals using an Illumina HiSeq2000. Using a computational pipeline for combining transcriptome and genome information, a total of 18,994 SNPs met the necessary minor allele frequency and depth filters. A series of further stringent filters were applied to identify those SNPs likely to succeed in genotyping assays, and for filtering of those in potential duplicated genome regions. A novel method for detecting potential intron-exon boundaries in areas of putative SNPs has also been applied in silico to improve genotyping success. In all, 2,317 filtered putative transcriptome SNPs suitable for genotyping primer design were identified. From those, a subset of 530 were selected, with the genotyping results showing the highest reported conversion and validation rates (91.3% and 83.2%, respectively) reported to date for a non-model species. This study represents a promising strategy to discover genotypable SNPs in the exome of non-model organisms. The genomic resource generated for E. encrasicolus, both in terms of sequences and novel markers, will be informative for research into this species with applications including traceability studies, population genetic analyses and aquaculture.

Multiple SNP markers reveal fine-scale population and deep phylogeographic structure in European anchovy (Engraulis encrasicolus L.)

Geographic surveys of allozymes, microsatellites, nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) have detected several genetic subdivisions among European anchovy populations. However, these studies have been limited in their power to detect some aspects of population structure by the use of a single or a few molecular markers, or by limited geographic sampling. We use a multi-marker approach, 47 nDNA and 15 mtDNA single nucleotide polymorphisms (SNPs), to analyze 626 European anchovies from the whole range of the species to resolve shallow and deep levels of population structure. Nuclear SNPs define 10 genetic entities within two larger genetically distinctive groups associated with oceanic variables and different life-history traits. MtDNA SNPs define two deep phylogroups that reflect ancient dispersals and colonizations. These markers define two ecological groups. One major group of Iberian-Atlantic populations is associated with upwelling areas on narrow continental shelves and includes populations spawning and overwintering in coastal areas. A second major group includes northern populations in the North East (NE) Atlantic (including the Bay of Biscay) and the Mediterranean and is associated with wide continental shelves with local larval retention currents. This group tends to spawn and overwinter in oceanic areas. These two groups encompass ten populations that differ from previously defined management stocks in the Alboran Sea, Iberian-Atlantic and Bay of Biscay regions. In addition, a new North Sea-English Channel stock is defined. SNPs indicate that some populations in the Bay of Biscay are genetically closer to North Western (NW) Mediterranean populations than to other populations in the NE Atlantic, likely due to colonizations of the Bay of Biscay and NW Mediterranean by migrants from a common ancestral population. Northern NE Atlantic populations were subsequently established by migrants from the Bay of Biscay. Populations along the Iberian-Atlantic coast appear to have been founded by secondary waves of migrants from a southern refuge.

Phylogeography of isolated freshwater three-spined stickleback Gasterosteus aculeatus populations in the Adriatic Sea basin

Analyses of mitochondrial (mt) DNA and microsatellite variation were carried out to examine the relationships between 10 freshwater populations of three-spined sticklebacks Gasterosteus aculeatus along the eastern coast of the Adriatic Sea. Partial sequences of the mtDNA control region and cytochrome b gene, in addition to 15 microsatellite loci, were used to analyse populations from four isolated river catchments. Results uncovered an Adriatic lineage that was clearly divergent from the European lineage, and confirmed that the most divergent and ancient populations are located within the Adriatic lineage as compared with other European populations. Two northern Adriatic populations formed independent clades within the European mitochondrial lineage, suggesting different colonization histories of the different Adriatic populations. Nuclear marker analyses also indicated deep divergence between Adriatic and European populations, albeit with some discordance between the mtDNA phylogeny of the northern Adriatic populations, further highlighting the strong differentiation among the Adriatic populations. The southern populations within the Adriatic lineage were further organized into distinct clades corresponding to respective river catchments and sub-clades corresponding to river tributaries, reflecting a high degree of population structuring within a small geographic region, concurrent with suggestions of existence of several microrefugia within the Balkan Peninsula. The highly divergent clades and haplotypes unique to the southern Adriatic populations further suggest, in accordance with an earlier, more limited survey, that southern Adriatic populations represent an important reservoir for ancient genetic diversity of G. aculeatus.

Genetic architecture of the marbled goby Pomatoschistus marmoratus (Perciformes, Gobiidae) in the Mediterranean Sea

The marbled goby Pomatoschistus marmoratus, a species inhabiting coastal Mediterranean lagoons, has been studied by measuring its mitochondrial DNA variation. This analysis revealed a Mediterranean west vs east split and, subsequently, an eastern differentiation among the Libyan-Tunisian Gulf, the Adriatic Sea and the Aegean Sea. The high cohesion between the samples collected in the vast area of western Mediterranean contrasts with the genetic mosaic of the more sub-structured eastern Mediterranean. This western homogeneity can not yet be fully explained even if a human-mediated migratory flow, due to a maritime traffic, has been posited. The pattern in the eastern basin revealed a genetic architecture possibly due to the non-migratory habit of the gobid. Within this perspective, the role of the Mediterranean lagoon habitat should be related to how much it amplifies the effects of historical (e.g. past sea-level changes) and environmental (e.g. present-day hydrographic regime) processes as regards the genetic structure of the inhabiting species.

Mitochondrial DNA variation in the caramote prawn Penaeus (Melicertus) kerathurus across a transition zone in the Mediterranean Sea

In this study we analysed mitochondrial DNA variation in Penaeus kerathurus prawns collected from seven locations along a transect across the Siculo-Tunisian region in order to verify if any population structuring exists over a limited geographical scale and to delineate the putative transition zone with sufficient accuracy. Partial DNA sequences of COI and 16S genes were analysed. In contrast to the highly conservative 16S gene, the COI sequences exhibited sufficient diversity for population analysis. The COI gene revealed low levels of haplotype and nucleotide diversities. The size of the annual landings of this commercial species suggests large population sizes. Hence, the low genetic diversity detected in this study could indicate a possible reduction in effective population sizes in the past. We detected significant genetic differentiation between eastern and western populations likely due to restricted gene flow across the Siculo-Tunisian boundary. We discuss the different evolutionary forces that may have shaped the genetic variation and suggest that the genetic divide is probably maintained by present-day dispersal limitation.

Relative role of life-history traits and historical factors in shaping genetic population structure of sardines (Sardina pilchardus)

Background

Marine pelagic fishes exhibit rather complex patterns of genetic differentiation, which are the result of both historical processes and present day gene flow. Comparative multi-locus analyses based on both nuclear and mitochondrial genetic markers are probably the most efficient and informative approach to discerning the relative role of historical events and life-history traits in shaping genetic heterogeneity. The European sardine (Sardina pilchardus) is a small pelagic fish with a relatively high migratory capability that is expected to show low levels of genetic differentiation among populations. Previous genetic studies based on meristic and mitochondrial control region haplotype frequency data supported the existence of two sardine subspecies (S. p. pilchardus and S. p. sardina).

Results

We investigated genetic structure of sardine among nine locations in the Atlantic Ocean and Mediterranean Sea using allelic size variation of eight specific microsatellite loci. Bayesian clustering and assignment tests, maximum likelihood estimates of migration rates, as well as classical genetic-variance-based methods (hierarchical AMOVA test and R_ST pairwise comparisons) supported a single evolutionary unit for sardines. These analyses only detected weak but significant genetic differentiation, which followed an isolation-by-distance pattern according to Mantel test.

Conclusion

We suggest that the discordant genetic structuring patterns inferred based on mitochondrial and microsatellite data might indicate that the two different classes of molecular markers may be reflecting different and complementary aspects of the evolutionary history of sardine. Mitochondrial data might be reflecting past isolation of sardine populations into two distinct groupings during Pleistocene whereas microsatellite data reveal the existence of present day gene flow among populations, and a pattern of isolation by distance.

Pillars of Hercules: is the Atlantic-Mediterranean transition a phylogeographical break?

The geological history of the Mediterranean Sea, its hydrography and connection with the Atlantic Ocean have been well documented. Despite a wealth of historical and oceanographic data, the Atlantic-Mediterranean transition remains controversial at the biological level as there are discordant results regarding the biogeographical separation between the Atlantic and Mediterranean biota. The opening of the Strait of Gibraltar at the end of the Messinian Salinity Crisis (some 5.33 million years ago), removed the land barrier that impeded the marine biota allowing it to disperse freely into the Mediterranean Sea. However, present day genetic patterns suggest a limitation to gene flow for some marine species, preventing population admixture. In the last few years, a large number of studies have challenged the hypothesis of the Strait of Gibraltar representing a phylogeographical break. A review of more than 70 papers reveals no obvious relationship between either dispersal ability or life history, and observed patterns of partial or complete genetic isolation between Atlantic and Mediterranean populations. We re-analysed a selection of this large body of data (20 studies in total) in order to provide a homogeneous and coherent view on the generality of the phylogeographical patterns and the presence of a phylogeographical barrier. This offered the opportunity to summarize the state of the art on this matter and reach some general conclusions on the evolutionary history across the Atlantic-Mediterranean range. Geographically, some species in the transition zone showed step changes of allele frequencies associated with the Almeria-Oran Front rather than with the Strait of Gibraltar itself. A major part of the data describe evolutionary events well within the time frame of the Quaternary age as very few taxa pre-date closure of the Tethys Sea. Results point to a combined signature of vicariance, palaeoclimate fluctuation and life-history traits on the Atlantic-Mediterranean phylogeographical patterns. Principal component analysis failed to show any particular association between biological traits and genetic variables. It would argue that organismal determinism may play a far less significant role than marine biogeographers have generally believed.

Demographic history and population structure of the Antarctic silverfish Pleuragramma antarcticum

The Antarctic silverfish Pleuragramma antarcticum (Nototheniidae) is the most abundant pelagic fish inhabiting Antarctic waters. In this study we investigated, through partial sequencing of the D-loop mitochondrial region, samples collected at four different locations in the Southern Ocean, three in the Atlantic and one in the Pacific sector. Sampling was replicated in two different years at two locations. Sequence analysis showed a remarkably high polymorphism, with 110 haplotypes over the 256 investigated specimens, and about 80% of haplotypes occurring only once. Neutrality tests indicated that all samples were not at mutation-drift equilibrium, and suggested a past population expansion. This result was supported by the presence of a star-like topology in the D-loop gene tree, and by results of mismatch distribution. The start of the expansion was dated, using a specifically calibrated clock, between 111 and 126 thousand years ago. This value corresponds to the start of the cooling period that led to the last glaciation peak, and is in close agreement with a recently suggested range expansion for pelagic Antarctic ecosystems. Analysis of molecular variation indicated a small, though highly significant, value of differentiation between samples. This result, together with the lack of association between clades and geographical locations, indicates a weak population structure for the species.