Contemporary and historic factors influence differently genetic differentiation and diversity in a tropical palm

Population genetics theory predicts loss in genetic variability because of drift and inbreeding in isolated plant populations; however, it has been argued that long-distance pollination and seed dispersal may be able to maintain gene flow, even in highly fragmented landscapes. We tested how historical effective population size, historical migration and contemporary landscape structure, such as forest cover, patch isolation and matrix resistance, affect genetic variability and differentiation of seedlings in a tropical palm (Euterpe edulis) in a human-modified rainforest. We sampled 16 sites within five landscapes in the Brazilian Atlantic forest and assessed genetic variability and differentiation using eight microsatellite loci. Using a model selection approach, none of the covariates explained the variation observed in inbreeding coefficients among populations. The variation in genetic diversity among sites was best explained by historical effective population size. Allelic richness was best explained by historical effective population size and matrix resistance, whereas genetic differentiation was explained by matrix resistance. Coalescence analysis revealed high historical migration between sites within landscapes and constant historical population sizes, showing that the genetic differentiation is most likely due to recent changes caused by habitat loss and fragmentation. Overall, recent landscape changes have a greater influence on among-population genetic variation than historical gene flow process. As immediate restoration actions in landscapes with low forest amount, the development of more permeable matrices to allow the movement of pollinators and seed dispersers may be an effective strategy to maintain microevolutionary processes.

Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru

The major Neotropical malaria vector, Anopheles darlingi, was reintroduced into the Iquitos, Loreto, Peru area during the early 1990s, where it displaced other anophelines and caused a major malaria epidemic. Since then, case numbers in Loreto have fluctuated, but annual increases have been reported since 2012. The population genetic structure of An. darlingi sampled before and after the introduction of long-lasting insecticidal nets (LLINs) was investigated to test the hypothesis of temporal population change (2006 vs. 2012). Current samples of An. darlingi were used to test the hypothesis of ecological adaptation to human modified (highway) compared with wild (riverine) habitat, linked to forest cover. In total, 693 An. darlingi from nine localities in Loreto, Peru area were genotyped using 13 microsatellite loci. To test the hypothesis of habitat differentiation in An. darlingi biting time patterns, HBR and EIR, four collections of An. darlingi from five localities (two riverine and three highway) were analysed. Analyses of microsatellite loci from seven (2006) and nine settlements (2012-2014) in the Iquitos area detected two distinctive populations with little overlap, although it is unclear whether this population replacement event is associated with LLIN distribution or climate. Within the 2012-2014 population two admixed subpopulations, A and B, were differentiated by habitat, with B significantly overrepresented in highway, and both in near-equal proportions in riverine. Both subpopulations had a signature of expansion and there was moderate genetic differentiation between them. Habitat and forest cover level had significant effects on HBR, such that Plasmodium transmission risk, as measured by EIR, in peridomestic riverine settlements was threefold higher than in peridomestic highway settlements. HBR was directly associated with available host biomass rather than forest cover. A population replacement event occurred between 2006 and 2012-2014, concurrently with LLIN distribution and a moderate El Niño event, and prior to an increase in malaria incidence. The likely drivers of this replacement cannot be determined with current data. The present-day An. darlingi population is composed of two highly admixed subpopulations, which appear to be in an early stage of differentiation, triggered by anthropogenic alterations to local habitat.

Social cohesion among kin, gene flow without dispersal and the evolution of population genetic structure in the killer whale (Orcinus orca)

In social species, breeding system and gregarious behavior are key factors influencing the evolution of large-scale population genetic structure. The killer whale is a highly social apex predator showing genetic differentiation in sympatry between populations of foraging specialists (ecotypes), and low levels of genetic diversity overall. Our comparative assessments of kinship, parentage and dispersal reveal high levels of kinship within local populations and ongoing male-mediated gene flow among them, including among ecotypes that are maximally divergent within the mtDNA phylogeny. Dispersal from natal populations was rare, implying that gene flow occurs without dispersal, as a result of reproduction during temporary interactions. Discordance between nuclear and mitochondrial phylogenies was consistent with earlier studies suggesting a stochastic basis for the magnitude of mtDNA differentiation between matrilines. Taken together our results show how the killer whale breeding system, coupled with social, dispersal and foraging behaviour, contributes to the evolution of population genetic structure.

Evidence for genetic differentiation of Octopus vulgaris (Mollusca, Cephalopoda) fishery populations from the southern coast of Brazil as revealed by microsatellites

Octopus vulgaris is a cephalopod species in several oceans and commonly caught by artisanal and industrial fisheries. In Brazil, O. vulgaris populations are mainly distributed along the southern coast and have been subjected to intensive fishing during recent years. Despite the importance of this marine resource, no genetic study has been carried out to examine genetic differences among populations along the coast of Brazil. In this study, 343 individuals collected by commercial vessels were genotyped at six microsatellite loci to investigate the genetic differences in O. vulgaris populations along the southern coast of Brazil. Genetic structure and levels of differentiation among sampling sites were estimated via a genotype assignment test and F-statistics. Our results indicate that the O. vulgaris stock consists of four genetic populations with an overall significant analogous F(ST). (phi(CT) = 0.10710, P<0.05) value. The genetic diversity was high with an observed heterozygosity of Ho = 0.987. The negative values of F(IS) found for most of the loci examined suggested a possible bottleneck process. These findings are important for further steps toward more sustainable octopus fisheries, so that this marine resource can be preserved for long-term utilization. (C) 2011 Elsevier B.V. All rights reserved.

Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti

Aedes aegypti is the most important vector of dengue viruses in tropical and subtropical regions. Because vaccines are still under development, dengue prevention depends primarily on vector control. Population genetics is a common approach in research involving Ae. aegypti. In the context of medical entomology, wing morphometric analysis has been proposed as a strong and low-cost complementary tool for investigating population structure. Therefore, we comparatively evaluated the genetic and phenotypic variability of population samples of Ae. aegypti from four sampling sites in the metropolitan area of Sao Paulo city, Brazil. The distances between the sites ranged from 7.1 to 50 km. This area, where knowledge on the population genetics of this mosquito is incipient, was chosen due to the thousands of dengue cases registered yearly. The analysed loci were polymorphic, and they revealed population structure (global F-ST = 0.062; p < 0.05) and low levels of gene flow (Nm = 0.47) between the four locations. Principal component and discriminant analyses of wing shape variables (18 landmarks) demonstrated that wing polymorphisms were only slightly more common between populations than within populations. Whereas microsatellites allowed for geographic differentiation, wing geometry failed to distinguish the samples. These data suggest that microevolution in this species may affect genetic and morphological characters to different degrees. In this case, wing shape was not validated as a marker for assessing population structure. According to the interpretation of a previous report, the wing shape of Ae. aegypti does not vary significantly because it is stabilised by selective pressure. (C) 2011 Elsevier B.V. All rights reserved.

Geographic differentiation of polymorphism in the Plasmodium falciparum malaria vaccine candidate gene SERA5

SERA5 is regarded as a promising malaria vaccine candidate of the most virulent human malaria parasite Plasmodium falciparum. SERA5 is a 120 kDa abundantly expressed blood-stage protein containing a papain-like protease. Since substantial polymorphism in blood-stage vaccine candidates may potentially limit their efficacy, it is imperative to fully investigate polymorphism of the SERA5 gene (sera5). In this study, we performed evolutionary and population genetic analysis of sera5. The level of inter-species divergence (kS = 0.076) between P. falciparum and Plasmodium reichenowi, a closely related chimpanzee malaria parasite is comparable to that of housekeeping protein genes. A signature of purifying selection was detected in the proenzyme and enzyme domains. Analysis of 445 near full-length P. falciparum sera5 sequences from nine countries in Africa, Southeast Asia, Oceania and South America revealed extensive variations in the number of octamer repeat (OR) and serine repeat (SR) regions as well as substantial level of single nucleotide polymorphism (SNP) in non-repeat regions (2562 bp). Remarkably, a 14 amino acid sequence of SERA5 (amino acids 59-72) that is known to be the in vitro target of parasite growth inhibitory antibodies was found to be perfectly conserved in all 445 worldwide isolates of P. falciparum evaluated. Unlike other major vaccine target antigen genes such as merozoite surface protein-1, apical membrane antigen-1 or circumsporozoite protein, no strong evidence for positive selection was detected for SNPs in the non-repeat regions of sera5. A biased geographical distribution was observed in SNPs as well as in the haplotypes of the sera5 OR and SR regions. In Africa, OR- and SR-haplotypes with low frequency (<5%) and SNPs with minor allele frequency (<5%) were abundant and were mostly continent-specific. Consistently, significant genetic differentiation, assessed by the Wright's fixation index (FST) of inter-population variance in allele frequencies, was detected for SNPs and both OR- and SR-haplotypes among almost all parasite populations. The exception was parasite populations between Tanzania and Ghana, suggesting frequent gene flow in Africa. The present study points to the importance of investigating whether biased geographical distribution for SNPs and repeat variants in the OR and SR regions affect the reactivity of human serum antibodies to variants. (C) 2011 Elsevier Ltd. All rights reserved.

Microsatellite dna reveals low ganetic differentiation among chub mackerel ("Scomber colias) sampled in Atlantic and Mediterranea

trabajo realizado por Medina Alcaraz, C., Castro, J.J., Sosa, P. A.

Peripheral arterial disease and diabetes: effects on endothelial progenitor cell differentiation and nitric oxide metabolism

In the recent years it is emerged that peripheral arterial disease (PAD) has become a growing health problem in Western countries. This is a progressive manifestation of atherothrombotic vascular disease, which results into the narrowing of the blood vessels of the lower limbs and, as final consequence, in critical leg ischemia. PAD often occurs along with other cardiovascular risk factors, including diabetes mellitus (DM), low-grade inflammation, hypertension, and lipid disorders. Patients with DM have an increased risk of developing PAD, and that risk increases with the duration of DM. Moreover, there is a growing population of patients identified with insulin resistance (IR), impaired glucose tolerance, and obesity, a pathological condition known as “metabolic syndrome”, which presents increased cardiovascular risk. Atherosclerosis is the earliest symptom of PAD and is a dynamic and progressive disease arising from the combination of endothelial dysfunction and inflammation. Endothelial dysfunction is a broad term that implies diminished production or availability of nitric oxide (NO) and/or an imbalance in the relative contribution of endothelium-derived relaxing factors. The secretion of these agents is considerably reduced in association with the major risks of atherosclerosis, especially hyperglycaemia and diabetes, and a reduced vascular repair has been observed in response to wound healing and to ischemia. Neovascularization does not only rely on the proliferation of local endothelial cells, but also involves bone marrow-derived stem cells, referred to as endothelial progenitor cells (EPCs), since they exhibit endothelial surface markers and properties. They can promote postnatal vasculogenesis by homing to, differentiating into an endothelial phenotype, proliferating and incorporating into new vessels. Consequently, EPCs are critical to endothelium maintenance and repair and their dysfunction contributes to vascular disease. The aim of this study has been the characterization of EPCs from healthy peripheral blood, in terms of proliferation, differentiation and function. Given the importance of NO in neovascularization and homing process, it has been investigated the expression of NO synthase (NOS) isoforms, eNOS, nNOS and iNOS, and the effects of their inhibition on EPC function. Moreover, it has been examined the expression of NADPH oxidase (Nox) isoforms which are the principal source of ROS in the cell. In fact, a number of evidences showed the correlation between ROS and NO metabolism, since oxidative stress causes NOS inactivation via enzyme uncoupling. In particular, it has been studied the expression of Nox2 and Nox4, constitutively expressed in endothelium, and Nox1. The second part of this research was focused on the study of EPCs under pathological conditions. Firstly, EPCs isolated from healthy subject were cultured in a hyperglycaemic medium, in order to evaluate the effects of high glucose concentration on EPCs. Secondly, EPCs were isolated from the peripheral blood of patients affected with PAD, both diabetic or not, and it was assessed their capacity to proliferate, differentiate, and to participate to neovasculogenesis. Furthermore, it was investigated the expression of NOS and Nox in these cells. Mononuclear cells isolated from peripheral blood of healthy patients, if cultured under differentiating conditions, differentiate into EPCs. These cells are not able to form capillary-like structures ex novo, but participate to vasculogenesis by incorporation into the new vessels formed by mature endothelial cells, such as HUVECs. With respect to NOS expression, these cells have high levels of iNOS, the inducible isoform of NOS, 3-4 fold higher than in HUVECs. While the endothelial isoform, eNOS, is poorly expressed in EPCs. The higher iNOS expression could be a form of compensation of lower eNOS levels. Under hyperglycaemic conditions, both iNOS and eNOS expression are enhanced compared to control EPCs, as resulted from experimental studies in animal models. In patients affected with PAD, the EPCs may act in different ways. Non-diabetic patients and diabetic patients with a higher vascular damage, evidenced by a higher number of circulating endothelial cells (CECs), show a reduced proliferation and ability to participate to vasculogenesis. On the other hand, diabetic patients with lower CEC number have proliferative and vasculogenic capacity more similar to healthy EPCs. eNOS levels in both patient types are equivalent to those of control, while iNOS expression is enhanced. Interestingly, nNOS is not detected in diabetic patients, analogously to other cell types in diabetics, which show a reduced or no nNOS expression. Concerning Nox expression, EPCs present higher levels of both Nox1 and Nox2, in comparison with HUVECs, while Nox4 is poorly expressed, probably because of uncompleted differentiation into an endothelial phenotype. Nox1 is more expressed in PAD patients, diabetic or not, than in controls, suggesting an increased ROS production. Nox2, instead, is lower in patients than in controls. Being Nox2 involved in cellular response to VEGF, its reduced expression can be referable to impaired vasculogenic potential of PAD patients.

Influence of pollination and seed dispersal on the genetic population structure of two Commiphora species in Madagascar and South Africa

Pollination and seed dispersal are important ecological processes for the regeneration of plant populations and both vectors for gene exchange between plant populations. For my thesis, I studied the pollination ecology of the South African tree Commiphora harveyi (Burseraceae) and compared it with C. guillauminii from Madagascar. Both species have low visitation rates and a low number of pollinating insect species, resulting in a low fruit set. While their pollination ecology is very similar, they differ in their seed dispersal with a low seed dispersal rate in the Malagasy and a high seed dispersal rate in the South African species. This should be reflected in a stronger genetic differentiation among populations in the Malagasy than in the South African species. My results, based on AFLP markers, contradict these expectations, the overall differentiation was lower in the Malagasy (FST = 0.05) than in the South African species (FST = 0.16). However, at a smaller spatial scale (below 3 km), the Malagasy species was genetically more strongly differentiated than the South African species, which was reflected by the high inter-population variance within the sample site (C. guillauminii: 72.2 - 85.5 %; C. harveyi: 8.4 - 14.5 %). This strong differentiation could arise from limited gene flow, which was confirmed by spatial autocorrelation analyses. The shape of the autocorrelogram suggested that gene exchange between individuals occurred only up to 3 km in the Malagasy species, whereas up to 30 km in the South African species. These results on the genetic structure correspond to the expectations based on seed dispersal data. Thus, seed dispersal seems to be a key factor for the genetic structure in plant populations on a local scale.

Genetic characterisation and global comparison of ancient Mediterranean Great White Sharks population (Carcharodon carcharias, L.1758)

This study poses as its objective the genetic characterization of the ancient population of the Great White shark, Carcharodon carcharias, L.1758, present in the Mediterranean Sea. Using historical evidence, for the most part buccal arches but also whole, stuffed examples from various national museums, research institutes and private collections, a dataset of 18 examples coming from the Mediterranean Sea has been created, in order to increase the informations regarding this species in the Mediterranean. The importance of the Mediterranean provenance derives from the fact that a genetic characterization of this species' population does not exist, and this creates gaps in the knowledge of this species in the Mediterranean. The genetic characterization of the individuals will initially take place by the extraction of the ancient DNA and the analysis of the variations in the sequence markers of the mitochondrial DNA. This approach has allowed the genetic comparison between ancient populations of the Mediterranean and contemporary populations of the same geographical area. In addition, the genetic characterization of the population of white sharks of the Mediterranean, has allowed a genetic comparison with populations from global "hot spots", using published sequences in online databases (NCBI, GenBank). Analyzing the variability of the dataset, both in terms space and time, I assessed the evolutionary relationships of the Mediterranean population of Great Whites with the global populations (Australia/New Zealand, South Africa, Pacific USA, West Atlantic), and the temporal trend of the Mediterranean population variability. This method based on the sequencing of two portions of mitochondrial DNA genes, markers showed us how the population of Great White Sharks in the Mediterranean, is genetically more similar to the populations of the Australia Pacific ocean, American Pacific Ocean, rather than the population of South Africa, and showing also how the population of South Africa is abnormally distant from all other clusters. Interestingly, these results are inconsistent with the results from tagging of this species. In addition, there is evidence of differences between the ancient population of the Mediterranean with the modern one. This differentiation between the ancient and modern population of white shark can be the result of events impacting on this species occurred over the last two centuries.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

Wide-scale population genomics of Atlantic bluefin tuna (Thunnus thynnus) inferred by novel high-throughput technology

My PhD project was focused on Atlantic bluefin tuna, Thunnus thynnus, a fishery resource overexploited in the last decades. For a better management of stocks, it was necessary to improve scientific knowledge of this species and to develop novel tools to avoid collapse of this important commercial resource. To do this, we used new high throughput sequencing technologies, as Next Generation Sequencing (NGS), and markers linked to expressed genes, as SNPs (Single Nucleotide Polymorphisms). In this work we applied a combined approach: transcriptomic resources were used to build cDNA libreries from mRNA isolated by muscle, and genomic resources allowed to create a reference backbone for this species lacking of reference genome. All cDNA reads, obtained from mRNA, were mapped against this genome and, employing several bioinformatics tools and different restricted parameters, we achieved a set of contigs to detect SNPs. Once a final panel of 384 SNPs was developed, following the selection criteria, it was genotyped in 960 individuals of Atlantic bluefin tuna, including all size/age classes, from larvae to adults, collected from the entire range of the species. The analysis of obtained data was aimed to evaluate the genetic diversity and the population structure of Thunnus thynnus. We detect a low but significant signal of genetic differentiation among spawning samples, that can suggest the presence of three genetically separate reproduction areas. The adult samples resulted instead genetically undifferentiated between them and from the spawning populations, indicating a presence of panmictic population of adult bluefin tuna in the Mediterranean Sea, without different meta populations.

Impact of anthropogenic stressors on the population genetics of the common bioindicator Dreissena polymorpha (Pallas, 1771)

Toxicant inputs from agriculture, industry and human settlements have been shown to severely affect freshwater ecosystems. Pollution can lead to changes in population genetic patterns through various genetic and stochastic processes. In my thesis, I investigated the impact of anthropogenic stressors on the population genetics of the zebra mussel Dreissena polymorpha. In order to analyze the genetics of zebra mussel populations, I isolated five new highly polymorphic microsatellite loci. Out of those and other already existing microsatellite markers for this species, I established a robust marker set of six microsatellite loci for D. polymorpha. rnMonitoring the biogeographical background is an important requirement when integrating population genetic measures into ecotoxicological studies. I analyzed the biogeographical background of eleven populations in a section of the River Danube (in Hungary and Croatia) and some of its tributaries, and another population in the River Rhine as genetic outgroup. Moreover, I measured abiotic water parameters at the sampling sites and analyzed if they were correlated with the genetic parameters of the populations. The genetic differentiation was basically consistent with the overall biogeographical history of the populations in the study region. However, the genetic diversity of the populations was not influenced by the geographical distance between the populations, but by the environmental factors oxygen and temperature and also by other unidentified factors. I found strong evidence that genetic adaptation of zebra mussel populations to local habitat conditions had influenced the genetic constitution of the populations. Moreover, by establishing the biogeographical baseline of molecular variance in the study area, I laid the foundation for interpreting population genetic results in ecotoxicological experiments in this region.rnIn a cooperation project with the Department of Zoology of the University of Zagreb, I elaborated an integrated approach in biomonitoring with D. polymorpha by combining the analysis techniques of microsatellite analysis, Comet assay and micronucleus test (MNT). This approach was applied in a case study on freshwater contamination by an effluent of a wastewater treatment plant (WWTP) in the River Drava (Croatia) and a complementary laboratory experiment. I assessed and compared the genetic status of two zebra mussel populations from a contaminated and a reference site. Microsatellite analysis suggested that the contaminated population had undergone a genetic bottleneck, caused by random genetic drift and selection, whereas a bottleneck was not detected in the reference population. The Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that the contaminated population had an increased percentage of micronuclei in hemocytes in comparison to the reference population. The laboratory experiment with mussels exposed to municipal wastewater revealed that mussels from the contaminated site had a lower percentage of tail DNA and a higher percentage of micronuclei than the reference population. These differences between populations were probably caused by an overall decreased fitness of mussels from the contaminated site due to genetic drift and by an enhanced DNA repair mechanism due to adaptation to pollution in the source habitat. Overall, the combination of the three biomarkers provided sufficient information on the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization.rnIn my thesis, I could show that the newly established marker set of six microsatellite loci provided reliable and informative data for population genetic analyses of D. polymorpha. The adaptation of the analyzed zebra mussel populations to the local conditions of their habitat had a strong influence on their genetic constitution. We found evidence that the different genetic constitutions of two populations had influenced the outcome of our ecotoxicological experiment. Overall, the integrated approach in biomonitoring gave comprehensive information about the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization and was well practicable in a first case study.

Temporal variation of genetic diversity and differentiation of albacore, Thunnus alalunga, in the Mediterranean Sea

This study is on albacore (Thunnus alalunga, Bonnaterre 1788), an epi- and mesopelagic oceanic tuna species cosmopolitan in the tropical and temperate waters of all oceans including the Mediterranean Sea, extending in a broad band between 40°N and 40°S. What it’s known about albacore population structure is based on different studies that used fisheries data, RFLP, mtDNA control region and nuDNA markers, blood lectins analysis, individual tags and microsatellite. At the moment, for T. alalunga six management units are recognized: the North Pacific, South Pacific, Indian, North Atlantic, South Atlantic and Mediterranean stocks. In this study I have done a temporal and spatial comparison of genetic variability between different Mediterranean populations of Thunnus alalunga matching an historical dataset ca. from 1920s composed of 43 individuals divided in 3 populations (NADR, SPAIN and CMED) with a modern dataset composed of 254 individuals and 7 populations (BAL, CYP, LIG, TYR, TUR, ADR, ALB). The investigation was possible using a panel of 94 nuclear SNPs, built specifically for the target species at the University of Basque Country UPV/EHU. First analysis done was the Hardy-Weinberg, then the number of clusters (K) was determined using STRUCTURE and to assess the genetic variability, allele frequencies, the average number of alleles per locus, expected (He) and observed (Ho) heterozygosis, and the index of polymorphism (P) was used the software Genetix. Historical and modern samples gives different results, showing a clear loss of genetic diversity over time leading to a single cluster in modern albacore instead of the two found in historical samples. What this study reveals is very important for conservation concerns, and additional research endeavours are needed.

Inference of Raja miraletus population structure using mitochondrial and nuclear DNA: comparing the resolution power of molecular markers in exploring species and population boundaries

This study focused on the role of oceanographic discontinuities and the presence of transitional areas in shaping the population structure and the phylogeography of the Raja miraletus species complex, coupled with the test of the effective occurrence of past speciation events. The comparisons between the Atlantic African and the North-Eastern Atlantic-Mediterranean geographic populations were unravelled using both Cytochrome Oxidase I and eight microsatellite loci. This approach guaranteed a robust dataset for the identification of a speciation event between the Atlantic African clade, corresponding to the ex Raja ocellifera nominal species, and the NE Atlantic-Mediterranean R. miraletus clade. As a matter of fact, the origin of the Atlantic Africa and the NE Atlantic-Mediterranean deep split dated about 11.74MYA and was likely due to the synergic influence currents and two upwelling areas crossing the Western African Waters. Within the Mediterranean Sea, particular attention was also paid to the transitional area represented by Adventura and Maltese Bank, that might have contributed in sustaining the connectivity of the Western and the Eastern Mediterranean geographical populations. Furthermore, the geology of the easternmost part of Sicily and the geo-morphological depression of the Calabrian Arc could have driven the differentiation of the Eastern Mediterranean Sea. Although bathymetric and oceanographic discontinuity could represent barriers to dispersal and migration between Eastern and Western Mediterranean samples, a clear and complete genetic separation among them was not detected. Results produced by this work identified a speciation event defining Raja ocellifera and R. miraletus as two different species, and describing the R. miraletus species complex as the most ancient cryptic speciation event in the family Rajidae, representing another example of how strictly connected the environment, the behavioural habits and the evolutionary and ecologic drivers are.