Mitochondrial DNA Introgression in the European Abalone Haliotis tuberculata tuberculata: Evidence for Experimental mtDNA Paternal Inheritance and a Natural Hybrid Sequence

Two subspecies of the European abalone have been morphologically recognized: Haliotis tuberculata tuberculata, present in the North Atlantic, and Haliotis tuberculata coccinea, present in the Canary Islands. Among the different nuclear markers used to differentiate these two subspecies, the sperm lysin gene was the most reliable, leading to a 2.2% divergence. Concerning the subunit I of the mitochondrial cytochrome oxydase gene (COI), we observed a difference of 3.3% between the two subspecies. In the North Atlantic, an introgression of mitochondrial DNA from H. tuberculata coccinea to H. tuberculata tuberculata was evident in around 30% of individuals. Due to this difference, we were able to experimentally detect the transfer of paternal mitochondrial DNA (mtDNA) by specific quantitative polymerase chain reaction measurements. The presence of the two mtDNA signatures was also detected in 20% of individuals tested in the field. Moreover, one mt DNA hybrid sequence was identified. The sequencing of this mitochondrial DNA hybrid revealed a mosaic structure with many specific mutations. The origin of this hybrid sequence is discussed.     

Multilocus approach reveals an incipient differentiation process in the Stone-curlew, <Emphasis Type="Italic">Burhinus oedicnemus</Emphasis> around the Mediterranean basin

The Stone-curlew Burhinus oedicnemus, a steppe bird species, is mainly distributed in the Mediterranean and Macaronesian regions, which are considered hotspots of biodiversity with priority for animal and plant species richness conservation. In this study, we investigated the genetic diversity of the Stone-curlew in the Mediterranean basin and in the Canary Islands by applying a multilocus approach. We analysed mitochondrial and nuclear markers in order to evaluate the genetic structure and the congruence between morphological subspecies and geographic samples. We found a significant level of genetic differentiation between Mediterranean and Canary Island populations with all markers. Both in the Mediterranean basin and in the Canary Islands, we found a significant level of genetic diversity with nuclear markers only. We identified seven population groups, including insular populations. The four subspecies described for the Western Palaearctic were confirmed with some changes in distribution range. In spite of habitat fragmentation and negative population trend, the Stone-curlew showed a significant level of genetic diversity and gene flow among continental populations. However, islands constitute important reservoirs of genetic diversity and a potential for the evolution of the species.     

Phylogeography of the common ragworm Hediste diversicolor (Polychaeta: Nereididae) reveals cryptic diversity and multiple colonization events across its distribution

Previous studies on the common ragworm Hediste diversicolor (Polychaeta: Nereididae) revealed a marked genetic fragmentation across its distribution and the occurrence of sibling taxa in the Baltic Sea. These results suggested that the phylogeographic patterns of H. diversicolor could reflect interactions between cryptic differentiation and multiple colonization events. This study aims to describe the large-scale genetic structuring of H. diversicolor and to trace the phylogeographic origins of the genetic types described in the Baltic Sea. Samples of H. diversicolor (2 <  n  < 28) were collected at 16 locations across the NE Atlantic coasts of Europe and Morocco and in the Mediterranean, Black and Caspian Seas and sequenced at two mitochondrial gene fragments (COI and cyt b , 345 and 290 bp, respectively). Bayesian analyses revealed deep phylogeographic splits yielding three main clades corresponding to populations (i) from the NE Atlantic coasts (from Germany to Morocco) and from part of the Western Mediterranean, (ii) from the Mediterranean Sea, and (iii) from the Black and Caspian Seas. These clades are further divided in well-supported subclades including populations from different regions of NE Atlantic and Mediterranean (i.e. Portugal/Morocco, Western Mediterranean, Adriatic Sea). The Baltic Sea comprises three sympatric lineages sharing a common evolutionary history with populations from NE Atlantic, Western Mediterranean and Black/Caspian Seas, respectively. Hence, the current patterns of genetic structuring of H. diversicolor appear as the result of allopatric isolation, multiple colonization events and possible adaptation to local environmental conditions.     

Phylogeography of the Atlanto-Mediterranean sea cucumber Holothuria (Holothuria) mammata: the combined effects of historical processes and current oceanographical pattern

We assessed the genetic structure of populations of the widely distributed sea cucumber Holothuria (Holothuria) mammata Grube, 1840, and investigated the effects of marine barriers to gene flow and historical processes. Several potential genetic breaks were considered, which would separate the Atlantic and Mediterranean basins, the isolated Macaronesian Islands from the other locations analysed, and the Western Mediterranean and Aegean Sea (Eastern Mediterranean). We analysed mitochondrial 16S and COI gene sequences from 177 individuals from four Atlantic locations and four Mediterranean locations. Haplotype diversity was high (H=0.9307 for 16S and 0.9203 for COI), and the haplotypes were closely related (π=0.0058 for 16S and 0.0071 for COI). The lowest genetic diversities were found in the Aegean Sea population. Our results showed that the COI gene was more variable and more useful for the detection of population structure than the 16S gene. The distribution of mtDNA haplotypes, the pairwise F(ST) values and the results of exact tests and amova revealed: (i) a significant genetic break between the population in the Aegean Sea and those in the other locations, as supported by both mitochondrial genes, and (ii) weak differentiation of the Canary and Azores Islands from the other populations; however, the populations from the Macaronesian Islands, Algarve and West Mediterranean could be considered to be a panmictic metapopulation. Isolation by distance was not identified in H. (H.) mammata. Historical events behind the observed findings, together with the current oceanographic patterns, were proposed and discussed as the main factors that determine the population structure and genetic signature of H. (H.) mammata.     

The distribution and molecular diversity of the Eastern Atlantic and Mediterranean chthamalids (Crustacea, Cirripedia)

The three chthamalids Chthamalus stellatus , C. montagui and Euraphia depressa are common inhabitants of the intertidal zone in the Eastern Atlantic, Mediterranean Sea and Black Sea. In this study, we investigated the occurrence of these barnacles in a wide range of their distribution. Population divergences of these two species have been inferred using three molecular markers — internal transcribed spacer (ITS), elongation factor 1α (EF-1α) and cytochrome oxidase subunit I (COI). ITS sequences of C. stellatus were identical throughout the species range, whereas ITS sequences of C. montagui indicated that the Black Sea and Mediterranean populations are isolated from the Atlantic population. The COI and EF-1α sequences were the most variable and informative. They indicated a high genetic divergence between Atlantic, Mediterranean and Black Sea populations for C. montagui . In addition significant genetic structure was found among the populations of C. stellatus based on EF-1α but not COI. Interestingly, our molecular dating analysis correlated the pattern of diversification in C. montagui to major geological changes that occurred in the Mediterranean during the end of the Messinian and Pleiocene periods. We suggest that palaeohistory shaped the divergences between Chthamalus populations that have probably been maintained by current hydrographic conditions. Finally, COI phylogenetic analysis placed the genus Euraphia within the Chthamalus clade, suggesting the need for a taxonomic revision of Euraphia . This study represents the most detailed phylogeographical analysis of intertidal Mediterranean species to date, and shows that geological events have strongly shaped the current diversity pattern of this fauna.     

EVOLUTIONARY AND ECOLOGICAL DIFFERENTIATION IN THE PANTROPICAL TO WARM-TEMPERATE SEAWEED DIGENEA SIMPLEX (RHODOPHYTA)1

Genetic differentiation among geographic isolates of the pantropical to warm-temperate red alga Digenea simplex (Wulfen) C. Agardh was investigated using random amplified polymorphic DNA (RAPD) markers, crossing studies, and temperature tolerances experiments. Eleven isolates representing populations from the Caribbean, eastern Atlantic, and Indo-West Pacific were compared. RAPD analysis clearly revealed an Indo-West Pacific group, a Caribbean/Cape Verde Islands group, and a Canary Islands group. Crossing studies showed different levels of inter fertility. In most crosses between Western Australian and Atlantic isolates, no hybrid tetrasporophytes were formed. In crosses between Caribbean and Cape Verde Islands isolates, tetrasporophytes developed, but the viability of tetraspores was reduced. Full sexual compatibility was observed among Cape Verde Islands isolates and among isolates from Bonaire. Temperature tolerance studies indicate that Pacific isolates have a broader temperature survival range than Atlantic isolates, which may be correlated to local temperature extremes. Despite the reduced level of sexual compatibility between Caribbean and Cape Verde Islands isolates, their shared position in the RAPD analysis and similar temperature responses suggest trans-Atlantic dispersal in the near geological past. In addition to their discrete position in the RAPD distance analysis, the Canary Islands isolates were significantly more cold-tolerant than the other Atlantic isolates. This finding is consistent with the hypothesis that the Canary Islands were recolonized from cold-adapted eastern Mediterranean populations after the last Pleistocene glaciation.     

Eutrophication and restoration of shallow lakes – the concept of stable equilibria revisited

We studied the genetic structure of populations of the Atlanto-Mediterranean ascidian Clavelina lepadiformis (Müller, 1776). A 369 bp segment of the COI mitochondrial gene was sequenced in Mediterranean and Atlantic populations from inside harbours, marinas and fjords (interior populations), and from the open rocky littoral (exterior populations). Previous work identified genetic differences between C. lepadiformis inhabiting Mediterranean harbours and the Mediterranean rocky littoral, however, the origin of these two clades remained speculative. Here we compared the Mediterranean populations with four Atlantic populations (two interior and two exterior). Gene differentiation and maximum likelihood analyses showed that the Atlantic forms were not divided into interior and exterior clades, and were closely related to the interior clade in the Mediterranean. The results support the hypothesis that both clades evolved allopatrically in the two seas, and that a recent colonisation of Mediterranean marinas from the Atlantic was caused by ship-hull transport. Colonisation of habitats by new genetic variants, morphologically indistinguishable from local populations, may be common among benthic invertebrates, and only genetic tools can uncover these cryptic invasions.     

Duplication and selection on abalone sperm lysin in an allopatric population

While gene duplication is a major source of evolutionary novelty, the importance of this process in reproductive protein evolution has not been widely investigated. Here, we report the first known case of gene duplication of abalone sperm lysin in an allopatric subspecies found in the Eastern Atlantic, Haliotis tuberculata coccinea. Mass spectrometry identified both copies of the lysin protein in testis tissue, and 3-dimensional structural modeling suggests that both proteins remain functional. We also detected positive selection acting on both paralogs after duplication and found evidence of a recent selective sweep. Because H. t. coccinea occurs in geographic isolation from other abalone species, these findings suggest that the evolution of lysin is not driven to create reproductive barriers to unfit hybrid formation with an overlapping species. Instead, sexual selection or sexual conflict acting during abalone fertilization could be responsible for the recent positive selection on this protein. The presence of multiple, rapidly evolving lysin genes in H. tuberculata presents an opportunity to study the early stages of diversification of a protein whose function is well understood.     

Forensic genetic identification of abalone (<Emphasis Type="Italic">Haliotis</Emphasis> spp.) of the northeastern Pacific Ocean

International trade in abalone meats, exclusive of their shells, is a lucrative business based upon both legally and illegally harvested abalone from many jurisdictions. The inability to visually identify abalone meat to species in the absence of the shell impedes enforcement efforts to reduce the illegal exploitation of the world’s abalone resources. We describe species-specific DNA sequences for the gamete recognition proteins, sperm lysin and vitelline egg receptor for lysin, and their use in forensic species identification among abalone of the northeastern Pacific Ocean. Some commercially relevant abalone species from the Northern and Southern hemispheres can be differentiated on the basis of the length of the second intron of the sperm lysin gene. The seven North American species of abalone that occupy the waters of Mexico, the USA and Canada can be distinguished based on sequence differentiation in the first three repeats of the vitelline receptor gene.     

Historical population dynamics and demography of the Eastern Atlantic pomacentrid Chromis limbata (Valenciennes, 1833)

Recent studies have focused on the relationship between the marine fauna of the Eastern Atlantic and the Mediterranean Sea, but within the Atlantic, little is known about genetic relationships between populations of the Macaronesian islands. In this study, we tested whether the paleo-climatology and paleo-oceanography of the region could predict the genetic relationships among three Eastern Atlantic populations (Azores, Madeira, and Canaries) of a damselfish, Chromis limbata, and compared our results with its Mediterranean and adjacent Atlantic sister species, Chromis chromis. We combined phylogeographic and coalescent approaches using the fast evolving mitochondrial control region gene. No population structure was found for the three archipelagos. The coalescence time estimated for C. limbata (0.857-1.17 Mya) was much greater than that estimated for C. chromis. We propose that this difference reflects differences in glaciating extents in the Northeastern Atlantic and the Mediterranean. Diversity indexes (Hd and genetic distances) together with historical demographic parameters of C. limbata (Theta and g) revealed a more stable population history when compared to C. chromis. Our results suggest that the Macaronesian populations of C. limbata have probably been less affected by the last glaciation than the Mediterranean populations of C. chromis. Migration across the three archipelagos was estimated and a prevailing northwest trend was detected. This result supports the idea of a colonization of the Azores by warm water fish from Madeira or the westernmost Canary islands which acted as major glacial refugia for the tropical and subtropical marine fauna during the glaciations.     

Phylogeography of the Calonectris shearwaters using molecular and morphometric data

We investigated phylogenetic relationships and the biogeographic history of the Calonectris species complex, using both molecular and biometric data from one population of the Cape Verde shearwater Calonectris edwardsii (Cape Verde Islands), one from the streaked shearwater C. leucomelas (western Pacific Ocean) and 26 from Cory's shearwater populations distributed across the Atlantic (C. d. borealis) and the Mediterranean (C. d. diomedea). The streaked shearwater appeared as the most basal and distant clades, whereas the genetic divergences among the three main clades within the Palearctic were similar. Clock calibrations match the first speciation event within Calonectris to the Panama Isthmus formation, suggesting a vicariant scenario for the divergence of the Pacific and the Palearctic clades. The separation between the Atlantic and Mediterranean clades would have occurred in allopatry by range contraction followed by local adaptation during the major biogeographic events of the Pleistocene. The endemic form from Cape Verde probably evolved as a result of ecological divergence from the Mediterranean subspecies. Finally, one Mediterranean population (Almeria) was unexpectedly grouped into the Atlantic subspecies clade, both by genetic and by morphometric analyses, pointing out the Almeria-Oran oceanographic front (AOOF) as the actual divide between the two Cory's shearwater subspecies. Our results highlight the importance of oceanographic boundaries as potentially effective barriers shaping population and species phylogeographical structure in pelagic seabirds.     

Patterns of diversity of the north-eastern Atlantic blenniid fish fauna (Pisces: Blenniidae)

This paper presents an analysis of the distributional patterns of blenniids (Pisces: Blenniidae) in the north‐eastern Atlantic. Two peaks of species diversity were found, both in terms of number of species and number of endemics: one in the tropical African coast and another in the Mediterranean Sea. A cluster analysis of similarity values (Jaccard coefficient) among the eastern Atlantic zoogeographical areas, revealed the following groups: a north temperate group, a tropical group formed by the tropical African coast and Mauritania, another group formed by the islands of Cape Verde, a south temperate group (South Africa), and a southern Atlantic group formed by the islands of Ascension and St Helena. Within the north temperate group, the subgroups with higher similarities were: Azores and Madeira, Canary Islands and Morocco, and the Mediterranean and the Atlantic coast of the Iberian Peninsula. Based on affinity indices, the probable directions of faunal flows were inferred. The tropical coast of Africa and the Mediterranean emerged from this analysis as probable speciation centres of the north‐eastern Atlantic blenniid fauna. The Mediterranean may have also acted as a refuge during glacial periods.     

Distributional and demographic consequences of Pleistocene climate fluctuations for a marine demersal fish in the north-eastern Atlantic

Aim The Pleistocene glaciations were the most significant historical event during the evolutionary life span of most extant species. However, little is known about the consequences of these climate changes for the distribution and demography of marine animals of the north‐eastern Atlantic. The present study focuses on the phylogeographic and demographic patterns of the sand goby, Pomatoschistus minutus (Teleostei: Gobiidae), a small marine demersal fish. Location North‐eastern Atlantic, Mediterranean, Irish, North and Baltic seas. Methods Analysis was carried out by sequencing the mtDNA cytochrome b gene of sand gobies from 12 localities throughout the species’ range, and using this information in combination with published data of allozyme markers and mtDNA control region sequences. Several phylogenetic methods and a network analysis were used to explore the phylogeographic pattern. The historical demography of P. minutus was studied through a mismatch analysis and a Bayesian skyline plot. Results Reciprocal monophyly was found between a Mediterranean Sea (MS) clade and an Atlantic Ocean (AO) clade, both with a Middle Pleistocene origin. The AO Clade contains two evolutionary significant units (ESUs): the Iberian Peninsula (IB) Group and the North Atlantic (NA) Group. These two groups diverged during Middle Pleistocene glacial cycles. For the NA Group there is evidence for geographic sorting of the ancestral haplotypes with recent radiations in the Baltic Sea, Irish Sea, North Sea and Bay of Biscay. The demographic histories of the Mediterranean Clade and the two Atlantic ESUs were influenced mainly by expansions dated as occurring during the Middle Pleistocene glaciations and post‐Eem, respectively. Main conclusions The pre‐LGM (Last Glacial Maximum) subdivision signals were not erased for P. minutus during the LGM. Middle Pleistocene glaciations yielded isolated and differently evolving sets of populations. In contrast to the case for most other taxa, only the northern Atlantic group contributed to the post‐glacial recolonization. The historical demography of Mediterranean sand gobies was influenced mainly by Middle Pleistocene glaciations, in contrast to that of the Atlantic populations, which was shaped by Late Pleistocene expansions.     

Phylogeography of the planktonic chaetognath Sagitta setosa reveals isolation in European seas

Numerous planktonic species have disjunct distribution patterns in the world's oceans. However, it is unclear whether these are truly unconnected by gene flow, or whether they are composed of morphologically cryptic species. The marine planktonic chaetognath Sagitta setosa Müller has a discontinuous geographic distribution over the continental shelf in the northeastern Atlantic, Mediterranean Sea, and Black Sea. Morphological variation between these populations has been described, but overlaps and is therefore unsuitable to determine the degree of isolation between populations. To test whether disjunct populations are also genetically disjunct, we sequenced a 504-bp fragment of mitochondrial DNA comprising the cytochrome oxidase II region of 86 individuals. Sequences were highly variable; each represented a different haplotype. Within S. setosa, sequence divergence ranged from 0.2 to 8.1% and strong phylogeographic structure was found, with four main groups corresponding to the northeastern Atlantic, Mediterranean Sea (including Ligurian Sea, Tyrrhenian Sea and Gulf of Gabes), Adriatic Sea, and Black Sea. Two of these (Atlantic and Black Sea) were resolved as monophyletic clades, thus gene flow between disjunct populations of S. setosa has been extremely limited and lineage sorting has taken place. The deepest divergence was between Atlantic and Mediterranean/Black Sea populations followed by a split between Mediterranean and Black Sea populations. The Mediterranean/Black Sea clade comprised three groups, with the Adriatic Sea as the most likely sister clade of the Black Sea. These data are consistent with a colonization of the Black Sea from the Mediterranean. Furthermore, a possible cryptic species was found in the Black Sea with 23.1% sequence divergence from S. setosa. Two possibilities for the evolutionary origin of this species are proposed, namely, that it represents a relict species from the ancient Paratethys, or that it represents another chaetognath species that colonized the Black Sea more recently. Even though the exact timing of disjunction of S. setosa populations remains unclear, on the basis of the geological and paleoclimatic history of the European basins and our estimates of net nucleotide divergence, we suggest that disjunct populations arose through vicariance resulting from the cyclical changes in temperature and sea levels during the Pleistocene. We conclude that these populations have remained disjunct, not because of limited dispersal ability, but because of the inability to maintain viable populations in suboptimal, geographically intermediate areas.     

Complex population genetic structure in the endemic Canary Island pine revealed using chloroplast microsatellite markers

The Canary archipelago, located on the northwestern Atlantic coast of Africa, is comprised of seven islands aligned from east to west, plus seven minor islets. All the islands were formed by volcanic eruptions and their geological history is well documented providing a historical framework to study colonization events. The Canary Island pine (Pinus canariensis C. Sm.), nowadays restricted to the westernmost Canary Islands (Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro), is considered an old (Lower Cretaceous) relic from an ancient Mediterranean evolutionary centre. Twenty seven chloroplast haplotypes were found in Canary Island pine but only one of them was common to all populations. The distribution of haplotypic variation in P. canariensis suggested the colonization of western Canary Islands from a single continental source located close to the Mediterranean Basin. Present-day populations of Canary Island pine retain levels of genetic diversity equivalent to those found in Mediterranean continental pine species, Pinus pinaster and Pinus halepensis. A hierarchical analysis of variance (AMOVA) showed high differentiation among populations within islands (approximately 19%) but no differentiation among islands. Simple differentiation models such as isolation by distance or stepping-stone colonization from older to younger islands were rejected based on product-moment correlations between pairwise genetic distances and both geographic distances and population-age divergences. However, the distribution of cpSSR diversity within the islands of Tenerife and Gran Canaria pointed towards the importance of the role played by regional Pliocene and Quaternary volcanic activity and long-distance gene flow in shaping the population genetic structure of the Canary Island pine. Therefore, conservation strategies at the population level are strongly recommended for this species.Communicated by D.B. NealeA. Gómez and S.C. González-Martínez as joint authors     

Genetic diversity and species identification in the endangered white abalone (<Emphasis Type="Italic">Haliotis sorenseni</Emphasis>)

In 2001, the white abalone Haliotis sorenseni became the first marine invertebrate in United States waters to receive federal protection as an endangered species. Prior to the endangered species listing, 20 abalone were collected as potential broodstock for a captive rearing program. Using DNA from these animals, we have developed genetic markers, including five nuclear microsatellite loci and partial sequences of one nuclear (VERL) and two mitochondrial (COI and CytB) genes, to assess genetic variability in the species, aid in species identification, and potentially track the success of future outplanting of captive-reared animals. All five microsatellite loci were polymorphic and followed expectations of simple Mendelian inheritance in laboratory crosses. Each of the wild-caught adult abalone exhibited a unique composite microsatellite genotype, suggesting that significant genetic variation remains in natural populations. A combination of nuclear and mitochondrial gene sequencing demonstrated that one of the original wild-caught animals was, in fact, not a white abalone, but H. kamtschatkana (possibly subspecies assimilis). Similarly, another animal of uncertain identity accidentally collected by dredging was also shown to be H. kamtschatkana. Inclusion of these two animals as broodstock could have resulted in unintentional hybridizations detrimental to the white abalone recovery program. Molecular genetic identifications will be useful both in preventing broodstock contamination and as markers for future restocking operations.     

Mitochondrial phylogeography of the Dalmatian wall lizard, Podarcis melisellensis (Lacertidae)

A 903 bp section of the mitochondrial cytochrome b gene was sequenced from 73 specimens of Podarcis melisellensis collected at 52 localities distributed over the major part of the species’ range. In addition, parts of the 12S (about 470 bp) and 16S rRNA (about 500 bp) genes were analysed for 11 representative samples leading to a congruent phylogeny. Our study includes representatives of all 20 subspecies recognized today. The phylogenetic analysis of the sequence data revealed three main clades: mainland with nearby islands, Vis archipelago, and Lastovo archipelago. The degree of mitochondrial DNA divergence among these clades suggests a separation of the respective population groups during the earliest Pleistocene. The phylogenetic pattern observed within the species is in sharp contrast to the actual taxonomic division into subspecies. A correlation between genetic diversity of P. melisellensis populations and paleogeography of the regions they inhabit is discussed.     

Cryptic divergence and strong population structure in the colonial invertebrate Pycnoclavella communis (Ascidiacea) inferred from molecular data

We studied sequence variation in the mitochondrial gene cytochrome c oxidase subunit I (COI) for 135 individuals from eight Mediterranean populations of the colonial ascidian Pycnoclavella communis across most of its presently known range of distribution in the Mediterranean. Three haplotypes from Atlantic locations were also included in the study. Phylogenetic, phylogeographic and population genetic analyses were used to unravel the genetic variability within and between populations. The study revealed 32 haplotypes for COI, 29 of them grouped within two Mediterranean lineages of P. communis (mean nucleotide divergence between lineages was 8.55%). Phylogenetic and network analyses suggest the possible existence of cryptic species corresponding to these two lineages. Population genetic analyses were restricted to the five populations belonging to the main genetic lineage, and for these localities we compared the information gleaned from COI sequence data and from eight microsatellite loci. A high genetic divergence between populations was substantiated using both kinds of markers (COI, global Fst=0.343; microsatellite loci, global Fst=0.362). There were high numbers of private haplotypes (COI) and alleles (microsatellites) in the populations studied. Restricted gene flow and inbreeding occur in the present range of distribution of the species. Microsatellite loci showed a strong incidence of failed amplifications, which we attribute to the marked intraspecies variability that hampered the application of these highly specific markers. Our results show important genetic variability at all levels studied, from within populations to between basins, possibly coupled to speciation processes. This variability is attributable to restricted gene flow among populations due to short-distance dispersal of the larvae.     

Genetic structure in the European endemic seabird,Phalacrocorax aristotelis,shaped by a complex interaction of historical and contemporary,physical and nonphysical drivers

Geographically separated populations tend to be less connected by gene flow, as a result of physical or nonphysical barriers preventing dispersal, and this can lead to genetic structure. In this context, highly mobile organisms such as seabirds are interesting because the small effect of physical barriers means nonphysical ones may be relatively more important. Here, we use microsatellite and mitochondrial data to explore the genetic structure and phylogeography of Atlantic and Mediterranean populations of a European endemic seabird, the European shag, Phalacrocorax aristotelis, and identify the primary drivers of their diversification. Analyses of mitochondrial markers revealed three phylogenetic lineages grouping the North Atlantic, Spanish/Corsican and eastern Mediterranean populations, apparently arising from fragmentation during the Pleistocene followed by range expansion. These traces of historical fragmentation were also evident in the genetic structure estimated by microsatellite markers, despite significant contemporary gene flow among adjacent populations. Stronger genetic structure, probably promoted by landscape, philopatry and local adaptation, was found among distant populations and those separated by physical and ecological barriers. This study highlights the enduring effect of Pleistocene climatic changes on shag populations, especially within the Mediterranean Basin, and suggests a role for cryptic northern refugia, as well as known southern refugia, on the genetic structure of European seabirds. Finally, it outlines how contemporary ecological barriers and behavioural traits may maintain population divergence, despite long‐distance dispersal triggered by extreme environmental conditions (e.g. population crashes).