Dilution-to-Extinction Platform for the Isolation of Marine Bacteria-Producing Antitumor Compounds.

Gram-negative marine bacteria are an underexplored source of new chemical entities for a wide range of applications. Even though, some have shown a high antitumor activity. This chapter describes an isolation and screening protocol based on the Dilution-to-Extinction approach coupled with an antiproliferative test oriented to the discovery of new cytotoxic compounds synthesized by marine bacteria. In addition to the discovery of new bioactive secondary metabolites, this protocol provides a high-throughput isolation and screening platform for discarding no bioactive strains during the first steps of the drug discovery process.

Identification of trans-AT polyketide clusters in two marine bacteria reveals cryptic similarities between distinct symbiosis factors.

Glutarimide-containing polyketides are known as potent antitumoral and antimetastatic agents. The associated gene clusters have only been identified in a few Streptomyces producers and Burkholderia gladioli symbiont. The new glutarimide-family polyketides, denominated sesbanimides D, E and F along with the previously known sesbanimide A and C, were isolated from two marine alphaproteobacteria Stappia indica PHM037 and Labrenzia aggregata PHM038. Structures of the isolated compounds were elucidated based on 1D and 2D homo and heteronuclear NMR analyses and ESI-MS spectrometry. All compounds exhibited strong antitumor activity in lung, breast and colorectal cancer cell lines. Subsequent whole genome sequencing and genome mining revealed the presence of the trans-AT PKS gene cluster responsible for the sesbanimide biosynthesis, described as sbn cluster. Strikingly, the modular architecture of downstream mixed type PKS/NRPS, SbnQ, revealed high similarity to PedH in pederin and Lab13 in labrenzin gene clusters, although those clusters are responsible for the production of structurally completely different molecules. The unexpected presence of SbnQ homologues in unrelated polyketide gene clusters across phylogenetically distant bacteria, raises intriguing questions about the evolutionary relationship between glutarimide-like and pederin-like pathways, as well as the functionality of their synthetic products.

Detection of a pederin-like compound using a dilution-to-extinction-based platform for the isolation of marine bacteria in drug discovery strategies.

The continued development of culturing technologies for the discovery of new molecules from marine microbes is of paramount importance for drug discovery. Coupled with this, the use of the high-throughput approach shows promise for increasing the number of Gram-negative and non-filamentous bacteria cultures that can be surveyed, since they show a lower potential of bioactivity. In this work, we propose a new strategy of high-throughput cultivation of bacteria inspired by a dilution-to-extinction (DTE) methodology for the isolation of, and screening for, new cytotoxic compound producing marine bacteria. A marine sponge tissue was directly used as inoculum and the results were compared with the data obtained through the direct plating isolation method. Enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR) genomic fingerprinting indicated the isolation of four bioactive strains, three of them producers of a pederin-like compound, and the fourth one able to synthesize a different compound, still unidentified, rendered by the DTE approach, in comparison with one bioactive strain identified through the plating method. Analyses based on the 16S rRNA gene data showed the existence of two different species belonging to the genus Labrenzia. The efficiency and diversity ratio in the number of isolates and compounds are discussed. In view of the results, the proposed DTE approach proved to be efficient for the isolation of new cytotoxic compounds of marine origin and pave the way for future potential applications.

Phylogeography and Population Genetic Analyses in the Iberian Toothcarp (Aphanius iberus Valenciennes, 1846) at Different Time Scales.

Secondary freshwater fish species inhabiting fluctuating and extreme environments are susceptible to changes in dispersion, effective population size, and genetic structure. The Iberian toothcarp Aphanius iberus is an endemic cyprinodontid of the Iberian Peninsula restricted to brackish water of salt marshes and coastal lagoons on the eastern Spanish Mediterranean coast. In this study, we analyzed mitochondrial cytochrome b (cyt b) DNA and microsatellite variation to evaluate ways in which the processes of extinction, dispersal, and colonization of A. iberus across its geographic distribution have affected its population genetic structure over time and space. The A. iberus network reconstruction indicated subtle levels of phylogeographic structuring. This, combined with substantial mitochondrial DNA (mtDNA) genetic diversity, suggests that Pleistocene glaciations had a lesser effect on the demographic structure of its populations than was the case for Iberian freshwater species with a similar distribution. Haplotype network, hierarchical analysis of molecular variance, and pairwise ΦST comparisons involving some Levantine samples showed a relatively high degree of mtDNA differentiation, which could be explained by historical isolation of the Villena Lagoon population. Conversely, significant genetic differentiation that follows an isolation-by-distance pattern, and a reduction in Ne though time was detected with microsatellites, suggesting extensive habitat fragmentation on the Mediterranean coast of the Iberian Peninsula over the past hundreds of years. At a smaller geographical scale (Mar Menor Lagoon), habitat fragmentation, probably due to human activity, appears to have resulted in substantially reduced migration and increased genetic drift, as shown by expanded genetic differentiation of populations.

Phylogenetic relationships of freshwater fishes of the genus Capoeta (Actinopterygii, Cyprinidae) in Iran.

The Middle East contains a great diversity of Capoeta species, but their taxonomy remains poorly described. We used mitochondrial history to examine diversity of the algae-scraping cyprinid Capoeta in Iran, applying the species-delimiting approaches General Mixed Yule-Coalescent (GMYC) and Poisson Tree Process (PTP) as well as haplotype network analyses. Using the BEAST program, we also examined temporal divergence patterns of Capoeta. The monophyly of the genus and the existence of three previously described main clades (Mesopotamian, Anatolian-Iranian, and Aralo-Caspian) were confirmed. However, the phylogeny proposed novel taxonomic findings within Capoeta. Results of GMYC, bPTP, and phylogenetic analyses were similar and suggested that species diversity in Iran is currently underestimated. At least four candidate species, Capoeta sp4, Capoeta sp5, Capoeta sp6, and Capoeta sp7, are awaiting description. Capoeta capoeta comprises a species complex with distinct genetic lineages. The divergence times of the three main Capoeta clades are estimated to have occurred around 15.6-12.4 Mya, consistent with a Mio-Pleistocene origin of the diversity of Capoeta in Iran. The changes in Caspian Sea levels associated with climate fluctuations and geomorphological events such as the uplift of the Zagros and Alborz Mountains may account for the complex speciation patterns in Capoeta in Iran.

Genetic and demographic recovery of an isolated population of brown bear Ursus arctos L., 1758.

The brown bear Ursus arctos L., 1758 population of the Cantabrian Mountains (northwestern Spain) became isolated from other bear populations in Europe about 500 years ago and has declined due to hunting and habitat degradation. At the beginning of the 20th century, the Cantabrian population split into eastern and western subpopulations, and genetic exchange between them ceased. In the early 1990s, total population size was estimated to be < 100 bears. Subsequently, reduction in human-caused mortality has brought about an increase in numbers, mainly in the western subpopulation, likely promoting male-mediated migration and gene flow from the western nucleus to the eastern. To evaluate the possible genetic recovery of the small and genetically depauperate eastern subpopulation, in 2013 and 2014 we genotyped hair and faeces samples (116 from the eastern subpopulation and 36 from the western) for 18 microsatellite markers. Data from the annual count of females with cubs of the year (COY) during the past twenty-six years was used to analyze demographic changes. The number of females with COY fell to a minimum of seven in the western and three in eastern subpopulations in the biennium 1993-1994 and reached a respective maximum of 54 and 10 individuals in 2013-2014. We also observed increased bear dispersal and gene flow, mainly from the western to the eastern subpopulation. Of the 26 unique genotypes detected in the eastern subpopulation, 14 (54%) presented an admixture composition, and seven (27%) were determined to be migrants from the western subpopulation. Hence, the two separated and clearly structured subpopulations identified in the past currently show some degree of genetic admixture. This research shows the partial demographic recovery and a change in genetic composition due to migration process in a population of bears that has been isolated for several centuries.

Contrasting demographic history and gene flow patterns of two mangrove species on either side of the Central American Isthmus.

Comparative phylogeography offers a unique opportunity to understand the interplay between past environmental events and life-history traits on diversification of unrelated but co-distributed species. Here, we examined the effects of the quaternary climate fluctuations and palaeomarine currents and present-day marine currents on the extant patterns of genetic diversity in the two most conspicuous mangrove species of the Neotropics. The black (Avicennia germinans, Avicenniaceae) and the red (Rhizophora mangle, Rhizophoraceae) mangroves have similar geographic ranges but are very distantly related and show striking differences on their life-history traits. We sampled 18 Atlantic and 26 Pacific locations for A. germinans (N = 292) and R. mangle (N = 422). We performed coalescence simulations using microsatellite diversity to test for evidence of population change associated with quaternary climate fluctuations. In addition, we examined whether patterns of genetic variation were consistent with the directions of major marine (historical and present day) currents in the region. Our demographic analysis was grounded within a phylogeographic framework provided by the sequence analysis of two chloroplasts and one flanking microsatellite region in a subsample of individuals. The two mangrove species shared similar biogeographic histories including: (1) strong genetic breaks between Atlantic and Pacific ocean basins associated with the final closure of the Central American Isthmus (CAI), (2) evidence for simultaneous population declines between the mid-Pleistocene and early Holocene, (3) asymmetric historical migration with higher gene flow from the Atlantic to the Pacific oceans following the direction of the palaeomarine current, and (4) contemporary gene flow between West Africa and South America following the major Atlantic Ocean currents. Despite the remarkable differences in life-history traits of mangrove species, which should have had a strong influence on seed dispersal capability and, thus, population connectivity, we found that vicariant events, climate fluctuations and marine currents have shaped the distribution of genetic diversity in strikingly similar ways.

Patterns of genetic variation in the endangered European mink (Mustela lutreola L., 1761).

BACKGROUND: The European mink (Mustela lutreola, L. 1761) is a critically endangered mustelid, which inhabits several main river drainages in Europe. Here, we assess the genetic variation of existing populations of this species, including new sampling sites and additional molecular markers (newly developed microsatellite loci specific to European mink) as compared to previous studies. Probabilistic analyses were used to examine genetic structure within and between existing populations, and to infer phylogeographic processes and past demography. RESULTS: According to both mitochondrial and nuclear microsatellite markers, Northeastern (Russia, Estonia and Belarus) and Southeastern (Romania) European populations showed the highest intraspecific diversity. In contrast, Western European (France and Spain) populations were the least polymorphic, featuring a unique mitochondrial DNA haplotype. The high differentiation values detected between Eastern and Western European populations could be the result of genetic drift in the latter due to population isolation and reduction. Genetic differences among populations were further supported by Bayesian clustering and two main groups were confirmed (Eastern vs. Western Europe) along with two contained subgroups at a more local scale (Northeastern vs. Southeastern Europe; France vs. Spain). CONCLUSIONS: Genetic data and performed analyses support a historical scenario of stable European mink populations, not affected by Quaternary climate oscillations in the Late Pleistocene, and posterior expansion events following river connections in both North- and Southeastern European populations. This suggests an eastern refuge during glacial maxima (as already proposed for boreal and continental species). In contrast, Western Europe was colonised more recently following either natural expansions or putative human introductions. Low levels of genetic diversity observed within each studied population suggest recent bottleneck events and stress the urgent need for conservation measures to counteract the demographic decline experienced by the European mink.

Island survivors: population genetic structure and demography of the critically endangered giant lizard of La Gomera, Gallotia bravoana.

BACKGROUND: The giant lizard of La Gomera (Gallotia bravoana), is an endemic lacertid of this Canary Island that lives confined to a very restricted area of occupancy in a steep cliff, and is catalogued as Critically Endangered by IUCN. We present the first population genetic analysis of the wild population as well as of captive-born individuals (for which paternity data are available) from a recovery center. Current genetic variability, and inferred past demographic changes were determined in order to discern the relative contribution of natural versus human-mediated effects on the observed decline in population size. RESULTS: Genetic analyses indicate that the only known natural population of the species shows low genetic diversity and acts as a single evolutionary unit. Demographic analyses inferred a prolonged decline of the species for at least 230 generations. Depending on the assumed generation time, the onset of the decline was dated between 1200-13000 years ago. Pedigree analyses of captive individuals suggest that reproductive behavior of the giant lizard of La Gomera may include polyandry, multiple paternity and female long-term sperm retention. CONCLUSIONS: The current low genetic diversity of G. bravoana is the result of a long-term gradual decline. Because generation time is unknown in this lizard and estimates had large credibility intervals, it is not possible to determine the relative contribution of humans in the collapse of the population. Shorter generation times would favor a stronger influence of human pressure whereas longer generation times would favor a climate-induced origin of the decline. In any case, our analyses show that the wild population has survived for a long period of time with low levels of genetic diversity and a small effective population size. Reproductive behavior may have acted as an important inbreeding avoidance mechanism allowing the species to elude extinction. Overall, our results suggest that the species retains its adaptive potential and could restore its ancient genetic diversity under favorable conditions. Therefore, management of the giant lizard of La Gomera should concentrate efforts on enhancing population growth rates through captive breeding of the species as well as on restoring the carrying capacity of its natural habitat.

Genetic diversity and population history of the endangered killifish Aphanius baeticus.

The secondary freshwater fish fauna of the western-Iberian Peninsula basin is primarily restricted to local coastal streams, and man-made salt evaporation ponds, etc., which are susceptible to periodical flood and drought events. Despite its uniqueness in ecological adaptation to high saltwater tolerance, very little is known about this fauna's population dynamics and evolutionary history. The killifish, Aphanius baeticus (Cyprinodontidae) is an endemic species restricted to river basins on Spain's southern Atlantic coastline, considered as "Endangered." In this study, the genetic structure, diversity and historical demography of A. baeticus were analyzed using mitochondrial (cytochrome b, N=131) and nuclear (4 out of 19 microsatellites tested, N=288) markers across its distribution range. The phylogenetic and networking reconstruction revealed subtle phylogeographic structuring. A scattered expansion at the beginning of the interglacial periods, coupled with posterior events of extinction and colonization caused by periodical cycles of flooding, could explain the absence of well-defined phylogenetic relationships among populations. Moreover, very low genetic diversity values and a weak population differentiation were detected. We proposed that dispersals allowed by periodic floods connecting river drainages may have promoted a wide genetic exchange among populations and could have contributed to the current genetic relatedness of these populations.

Microsatellite DNA capture from enriched libraries.

Microsatellites are DNA sequences of tandem repeats of one to six nucleotides, which are highly polymorphic, and thus the molecular markers of choice in many kinship, population genetic, and conservation studies. There have been significant technical improvements since the early methods for microsatellite isolation were developed, and today the most common procedures take advantage of the hybrid capture methods of enriched-targeted microsatellite DNA. Furthermore, recent advents in sequencing technologies (i.e., next-generation sequencing, NGS) have fostered the mining of microsatellite markers in non-model organisms, affording a cost-effective way of obtaining a large amount of sequence data potentially useful for loci characterization. The rapid improvements of NGS platforms together with the increase in available microsatellite information open new avenues to the understanding of the evolutionary forces that shape genetic structuring in wild populations. Here, we provide detailed methodological procedures for microsatellite isolation based on the screening of GT microsatellite-enriched libraries, either by cloning and Sanger sequencing of positive clones or by direct NGS. Guides for designing new species-specific primers and basic genotyping are also given.

Evolutionary history of the genus Trisopterus.

The group of small poor cods and pouts from the genus Trisopterus, belonging to the Gadidae family, comprises four described benthopelagic species that occur across the North-eastern Atlantic, from the Baltic Sea to the coast of Morocco, and the Mediterranean. Here, we combined molecular data from mitochondrial (cytochrome b) and nuclear (rhodopsin) genes to confirm the taxonomic status of the described species and to disentangle the evolutionary history of the genus. Our analyses supported the monophyly of the genus Trisopterus and confirmed the recently described species Trisopterus capelanus. A relaxed molecular clock analysis estimated an Oligocene origin for the group (~30 million years ago; mya) indicating this genus as one of the most ancestral within the Gadidae family. The closure and re-opening of the Strait of Gibraltar after the Messinian Salinity Crisis (MSC) probably triggered the speciation process that resulted in the recently described T. capelanus.

Population proteomics of the European Hake (Merluccius merluccius).

We report the novel use of proteomics to investigate protein variation among populations of the European hake (Merluccius merluccius). The liver and brain extracts of 18 hake (N = 36) captured in the Mediterranean Sea, Cantabrian Sea, and Atlantic Ocean were examined by 2D/DIGE and mass spectrometry. Significant differences in protein expression among populations were revealed by 84 spots obtained in the gels for the liver and 145 spots for the brain. Population groups of samples were defined by multivariate analysis (PCA and hierarchical clustering). According to protein expression levels and the functions of the 55 candidate protein spots identified, which showed significant expression differences, highest population discrimination was rendered by brain proteins involved in cell signaling and metabolism/energy and by liver proteins involved in protein fate. Finally, we present a statistically robust framework to accurately classify individuals according to their population of origin. Thus, purposely identified protein isoforms were found to be competent at discriminating populations. These results suggest the possibility of identifying protein biomarkers related to environmental changes in a nonmodel species such as the hake and pave the way for more extensive research on protein variation among populations of marine fishes.

Genetic structuring and migration patterns of Atlantic bigeye tuna, Thunnus obesus (Lowe, 1839).

BACKGROUND: Large pelagic fishes are generally thought to have little population genetic structuring based on their cosmopolitan distribution, large population sizes and high dispersal capacities. However, gene flow can be influenced by ecological (e.g. homing behaviour) and physical (e.g. present-day ocean currents, past changes in sea temperature and levels) factors. In this regard, Atlantic bigeye tuna shows an interesting genetic structuring pattern with two highly divergent mitochondrial clades (Clades I and II), which are assumed to have been originated during the last Pleistocene glacial maxima. We assess genetic structure patterns of Atlantic bigeye tuna at the nuclear level, and compare them with mitochondrial evidence. RESULTS: We examined allele size variation of nine microsatellite loci in 380 individuals from the Gulf of Guinea, Canary, Azores, Canada, Indian Ocean, and Pacific Ocean. To investigate temporal stability of genetic structure, three Atlantic Ocean sites were re-sampled a second year. Hierarchical AMOVA tests, RST pairwise comparisons, isolation by distance (Mantel) tests, Bayesian clustering analyses, and coalescence-based migration rate inferences supported unrestricted gene flow within the Atlantic Ocean at the nuclear level, and therefore interbreeding between individuals belonging to both mitochondrial clades. Moreover, departures from HWE in several loci were inferred for the samples of Guinea, and attributed to a Wahlund effect supporting the role of this region as a spawning and nursery area. Our microsatellite data supported a single worldwide panmictic unit for bigeye tunas. Despite the strong Agulhas Current, immigration rates seem to be higher from the Atlantic Ocean into the Indo-Pacific Ocean, but the actual number of individuals moving per generation is relatively low compared to the large population sizes inhabiting each ocean basin. CONCLUSION: Lack of congruence between mt and nuclear evidences, which is also found in other species, most likely reflects past events of isolation and secondary contact. Given the inferred relatively low number of immigrants per generation around the Cape of Good Hope, the proportions of the mitochondrial clades in the different oceans may keep stable, and it seems plausible that the presence of individuals belonging to the mt Clade I in the Atlantic Ocean may be due to extensive migrations that predated the last glaciation.

As the raven flies: using genetic data to infer the history of invasive common raven (Corvus corax) populations in the Mojave Desert.

Common raven (Corvus corax) populations in Mojave Desert regions of southern California and Nevada have increased dramatically over the past five decades. This growth has been attributed to increased human development in the region, as ravens have a commensal relationship with humans and feed extensively at landfills and on road-killed wildlife. Ravens, as a partially subsidized predator, also represent a problem for native desert wildlife, in particular threatened desert tortoises (Gopherus agassizii). However, it is unclear whether the more than 15-fold population increase is due to in situ population growth or to immigration from adjacent regions where ravens have been historically common. Ravens were sampled for genetic analysis at several local sites within five major areas: the West Mojave Desert (California), East Mojave Desert (southern Nevada), southern coastal California, northern coastal California (Bay Area), and northern Nevada (Great Basin). Analyses of mtDNA control region sequences reveal an increased frequency of raven 'Holarctic clade' haplotypes from south to north inland, with 'California clade' haplotypes nearly fixed in the California populations. There was significant structuring among regions for mtDNA, with high F(ST) values among sampling regions, especially between the Nevada and California samples. Analyses of eight microsatellite loci reveal a mostly similar pattern of regional population structure, with considerably smaller, but mostly significant, values. The greater mtDNA divergences may be due to lower female dispersal relative to males, lower N(e), or effects of high mutation rates on maximal values of F(ST). Analyses indicate recent population growth in the West Mojave Desert and a bottleneck in the northern California populations. While we cannot rule out in situ population growth as a factor, patterns of movement inferred from our data suggest that the increase in raven populations in the West Mojave Desert resulted from movements from southern California and the Central Valley. Ravens in the East Mojave Desert are more similar to ones from northern Nevada, indicating movement between those regions. If this interpretation of high gene flow into the Mojave Desert is correct, then efforts to manage raven numbers by local control may not be optimally effective.

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 RST 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.

Genetic diversity and historical demography of Atlantic bigeye tuna (Thunnus obesus).

Bigeye (Thunnus obesus) is a large, pelagic, and migratory species of tuna that inhabits tropical and temperate marine waters worldwide. Previous studies based on mitochondrial RFLP data have shown that bigeye tunas from the Atlantic Ocean are the most interesting from a genetic point of view. Two highly divergent mitochondrial haplotype clades (I and II) coexist in the Atlantic Ocean. One is almost exclusive of the Atlantic Ocean whereas the other is also found in the Indo-Pacific Ocean. Bigeye tuna from the Atlantic Ocean is currently managed as a single stock, although this assumption remains untested at the genetic level. Therefore, genetic diversity was determined at the mitochondrial control region to test the null hypothesis of no population structure in bigeye tuna from the Atlantic Ocean. A total of 331 specimens were sampled from four locations in the Atlantic Ocean (Canada, Azores, Canary Islands, and Gulf of Guinea), and one in the Indian and Pacific Oceans, respectively. The reconstructed neighbor-joining phylogeny confirmed the presence of Clades I and II throughout the Atlantic Ocean. No apparent latitudinal gradient of the proportions of both clades in the different collection sites was observed. Hierarchical AMOVA tests and pairwise phi(ST) comparisons involving Atlantic Ocean Clades I and II were consistent with a single stock of bigeye tuna in the Atlantic Ocean. Population genetic analyses considering phylogroups independently supported gene flow within Clade II throughout the Atlantic Ocean, and within Clade I between Atlantic and Indo-Pacific Oceans. The latter result suggests present uni-directional gene flow from the Indo-Pacific into the Atlantic Ocean. Moreover, mismatch analyses dated divergence of Clades I and II during the Pleistocene, as previously proposed. In addition, migration rates were estimated using coalescent methods, and showed a net migration from Atlantic Ocean feeding grounds towards the Gulf of Guinea, the best-known spawning ground of Atlantic bigeye tuna.

Signature of an early genetic bottleneck in a population of Moroccan sardines (Sardina pilchardus).

Fishery assessment models meant to determine sustainability of commercial marine fish failed to predict recent stock collapses due to overexploitation. One flaw of assessment models is that they strongly rely on catch and age-composition statistics, but largely ignore the genetic background of the studied populations. We examined population genetic structure of sardine (Sardina pilchardus) in the centraleastern and northeastern Atlantic Ocean and Mediterranean Sea to aid fishery management of this heavily fished small pelagic species. We found that sardine has a striking mitochondrial control region, and sequenced a fragment of 387 bp of its 5'-end in 261 individuals collected off the coasts of Morocco (Dakhla, Tantan, Safi, Larache, and Nador), Portugal (Quarteira), Spain (Pasajes, Barcelona), and Greece (Kavala). High levels of haplotypic diversity rendered a rather unresolved NJ phylogeny. The recovered tree had no phylogeographic structuring except for the clustering of 13 individuals of Safi. In contrast, individuals grouped together according to the presence or absence of a 13-bp insertion in the sequence. Phi(ST) pairwise comparisons and molecular variance analyses supported genetic differentiation between the population of Pasajes (Bay of Biscay), and those of the Mediterranean Sea and Moroccan coast, with a contact zone around the Strait of Gibraltar. This result confirms the existence of two subspecies, S. pilchardus pilchardus and S. pilchardus sardina that were previously identified based on meristics and morphometry. Mismatch distribution analysis showed that sardine populations are expanding since the Pleistocene. Surprisingly, the population of Safi showed strong and statistically significant levels of genetic differentiation that could be related with isolation and genetic drift. Comparative analysis of the Safi population versus the rest including mismatch distributions, and a Bayesian skyline plot suggest that the Safi population likely underwent an early genetic bottleneck. The genetic singularity of the Safi population could have been responsible for the historical collapse of this sardine stock in the 1970s.