Molecular phylogeny and historical zoogeography of four genera of Eastern Atlantic skates (Rajiformes, Rajidae)

The aim of this study was to reconstruct a solid phylogeny of four genera of the Rajidae family (Chondrichthyans: Batoidea) using a concatenated alignment of mtDNA genes. Then use the resultant tree to estimate divergence time between taxa based on molecular clock and fossil calibration and conduct biogeographic analysis. The intent was to prove that the actual distribution of species of Eastern Atlantic and Mediterranean skates is due to a series of vicariant events. The species considered belongs to two different tribe: Rajini (Raja and Dipturus) and Amblyrajini (Leucoraja and Rajella). The choice of this genera is due to their high presence in the area of interest and to the richness of endemic species. The results show that despite the ancient origin of Rajidae (97 MYA), the Eastern Atlantic and Mediterranean faunas originated more recently, during Middle Miocene-Late Pliocene, after the closure of connection between these areas and the Indo-Pacific ocean (15 MYA). The endemic species of the Mediterranean (Raja asterias, R. radula, R. polystigma and Leucoraja melitensis) originated after the Messinian salinity crisis (7-5 MYA), when the recolonization of the basin occurred, and are still maintained in allopatric distribution by the presence of biogeographic barriers. Moreover from 4 to 2.6 MYA we can observe the formation of sister species for Raja, Leucoraja and Rajella, one of which has a Northern distribution, and the other has a Southern distribution (R. clavata vs R. straeleni, L. wallacei vs L. naevus, R. fyllae vs R. caudaspinosa and R. kukujevi vs R. leopardus + R. barnardi). The Quaternary and present oceanographic discontinuities that occur along the western African continental shelf (e.g., Cape Blanc and the Angola–Benguela Front) might contribute to the maintenance of low or null levels of gene flow between these closely related siblings species. Also sympatric speciation must be invoked to explain the evolution of skates, for example for the division between R. leopardus and R. barnardi. The speciation processes followed a south-to-north pathways for Dipturus and a north-to-south pathways for Raja, Leucoraja and Rajella underling that the evolution of the genera occurred independently. In the end, it is conceivable that the evolutionary pathways of the tribes followed the costal line during the gondwana fragmentation. The results demonstrate that the evolution of this family is characterized by a series of parallel and independent speciation events, strictly correlated to the tectonic movement of continental masses and paleogeographic and paleoclimatic events and so can be explained by a panbiogeographical (vicariance) model.

Evolution and historical biogeography of the Eastern Atlantic skates

In my thesis, I tested the hypothesis that the diversification of the Eastern Atlantic skate faunas arose through vicariance rather than dispersal, using combined approach of molecular phylogeny reconstruction and zoogeography (namely historical biogeography). This analyses have been carried out independently on four Rajidae genera belonging to two different tribes: Rajini (Raja and Dipturus) and Amblyrajini (Rajella and Leucoraja). These taxa were selected because they displayed high species diversity and richness of endemic species in the Eastern Atlantic and Mediterranean. The verification of this hypothesis was carried out by reconstructing the best phylogenetic relationships among four genera and 26 species (including several endemism) based on mtDNA and nuDNA gene variation and several statistical approaches. Divergence times of taxa have been estimated based on molecular clock and fossil calibration to explain evolutionary patterns in the context of geological framework. Main issues are (i) the evidence that Eastern Atlantic skate evolution and displacement of species diversity occurred from pulsed geographical speciation (i.e. repeated series of parallel and independent speciation events) started in the Late Eocene-Early Miocene and they have occurred prevalently during Miocene; (ii) such relatively ancient origin of diversification has been allowed the sympatric displacement and evolution of several congeneric taxa likely because they have accumulated huge differences in the genomic and physiological/behavioural phenotypic traits; (iii) recently diverged sister species and taxa showed allopatric or parapatric evolution by the presence of oceanographic or hydrogeographical barriers which likely prevent large mixing between parapatric sister species.

First assessment on genetic structure and phylogeography of Mediterranean blue shark (prionace glauca, L. 1758) population using mitochondrial gene variation: a comparison with the Atlantic.

The blue shark, Prionace glauca, is one of the most vagile shark species worldwide distributed. The particular body shape allows blue sharks make transoceanic movements, leading to a circumglobal distribution. Due to its reproductive cycle, an extraordinarily high number of specimens is globally registered but, even if it is still a major bycatch of longline fishery rather than a commercial target, it is characterized by a high vulnerability. In this perspective it is important to increase the amount of informations regarding its population extent in the different worldwide areas, evaluating the possible phylogeographic patterns between different locations. This study, included in the "MedBlueSGen" European project, aims exactly at filling a gap in knowledges regarding the genetic population structure of the Mediterranean blue sharks, which has never been investigated before, with a comparison with the North-Eastern Atlantic blue shark population. To reach this objective, we used a dataset of samples from different Mediterranean areas implementing it with some samples from North-Eastern Atlantic. Analyzing the variability of the two mitochondrial markers control region and cytochrome b, with the design of new species-specific primer pairs, we assessed the mitochondrial genetic structure of Mediterranean and North-Eastern Atlantic samples, focusing on the analysis of their possible connectivity, and we tried to reconstruct their demographic history and population size. Data analyses highlighted the absence of a genetic structuring within the Mediterranean and among it and North-Eastern Atlantic, suggesting that the Strait of Gibraltar doesn't represent a phylogeographic barrier. These results are coherent to what has been found in similar investigations on other worldwide blue shark populations. Analysis of the historical demographic trend revealed a general stable pattern for the cytochrome-b and a slightly population expansion for the control region marker.

Analysis of the Variability of Carbonate System Parameters in the North-East Atlantic.

The purpose of this thesis is to analyse the spatial and temporal variability of the aragonite saturation state (ΩAR), commonly used as an indicator of ocean acidification, in the North-East Atlantic. When the aragonite saturation state decreases below a certain threshold, ΩAR <1, calcifying organisms (i.e. molluscs, pteropods, foraminifera, crabs, etc.) are subject to dissolution of shells and aragonite structures. This objective agrees with the challenge 'Ocean, climate change and acidification' of the EU COST Ocean Governance for Sustainability project, which aims to combine the information collected on the state of health of the oceans. Two open-sources data products, EMODnet and GLODAPv2, have been integrated and analysed for the first time in the North-East Atlantic region. The integrated dataset contains 1038 ΩAR vertical profiles whose time distribution spans from 1970 to 2014. The ΩAR has been computed from CO2SYS software considering different combinations of input parameters, pH, Total Alkalinity (TAlk) and Dissolved Inorganic Carbon (DIC), associated with Temperature, Salinity and Pressure at in situ conditions. A sensitivity analysis has been performed to better understand the data consistency of ΩAR computed from the different combinations of pH, Talk and DIC and to verify the difference among observed TAlk and DIC parameters and their output values from the CO2SYS tool. Maps of ΩAR have been computed with the best data coverage obtained from the two datasets, at different levels of depth in the area of investigation and they have been compared to the work of Jiang et al. (2015). The results are consistent and show similar horizontal and vertical patterns. The study highlights some aragonite undersaturated values (ΩAR <1) below 500 meters depth, suggesting a potential effect of acidification in the considered time period. This thesis aims to be a preliminary work for future studies that will be able to design the ΩAR variability on a decadal distribution based on the extended time-series acquired in this work.

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.

Genetic and morphological variation and differentiation of Raja clavata populations in the European seas: marker development and preliminary data

In this study the population structure and connectivity of the Mediterranean and Atlantic Raja clavata (L., 1758) were investigated by analyzing the genetic variation of six population samples (N = 144) at seven nuclear microsatellite loci. The genetic dataset was generated by selecting population samples available in the tissue databases of the GenoDREAM laboratory (University of Bologna) and of the Department of Life Sciences and Environment (University of Cagliari), all collected during past scientific surveys (MEDITS, GRUND) from different geographical locations in the Mediterranean basin and North-east Atlantic sea, as North Sea, Sardinian coasts, Tuscany coasts and Cyprus Island. This thesis deals with to estimate the genetic diversity and differentiation among 6 geographical samples, in particular, to assess the presence of any barrier (geographic, hydrogeological or biological) to gene flow evaluating both the genetic diversity (nucleotide diversity, observed and expected heterozygosity, Hardy- Weinberg equilibrium analysis) and population differentiation (Fst estimates, population structure analysis). In addition to molecular analysis, quantitative representation and statistical analysis of morphological individuals shape are performed using geometric morphometrics methods and statistical tests. Geometric coordinates call landmarks are fixed in 158 individuals belonging to two population samples of Raja clavata and in population samples of closely related species, Raja straeleni (cryptic sibling) and Raja asterias, to assess significant morphological differences at multiple taxonomic levels. The results obtained from the analysis of the microsatellite dataset suggested a geographic and genetic separation between populations from Central-Western and Eastern Mediterranean basins. Furthermore, the analysis also showed that there was no separation between geographic samples from North Atlantic Ocean and central-Western Mediterranean, grouping them to a panmictic population. The Landmark-based geometric morphometry method results showed significant differences of body shape able to discriminate taxa at tested levels (from species to populations).

Genetic patterns in recovered specimens of Caretta caretta from the Adriatic Sea

The Adriatic Sea is considered a feeding and developmental area for Mediterranean loggerhead turtles, but this area is severely threatened by human impacts. In the Adriatic Sea loggerhead turtles are often found stranded or floating, but they are also recovered as by-catch from fishing activities. Nevertheless, information about population structuring and origin of individuals found in the Adriatic Sea are still limited. Cooperation with fishermen and a good network of voluntary collaborators are essential for understanding their distribution, ecology and for developing conservation strategies in the Adriatic Sea. In this study, a comparative analysis of biometric data and DNA sequence polymorphism of the long fragment of the mitochondrial control region was carried out on ninety-three loggerheads recovered from three feeding areas in the Adriatic Sea: North-western, North-eastern and South Adriatic. Differences in turtles body sizes (e.g. Straight Carapace Length) among the three recovery areas and relationship between SCL and the type of recovery were investigated. The origin of turtles from Mediterranean rookeries and the use of the Adriatic feeding habitats by loggerheads in different life-stages were assessed to understand the migratory pathway of the species. The analysis of biometric data revealed a significant difference in turtle sizes between the Southern and the Northern Adriatic. Moreover, size of captured turtles resulted significantly different from the size of stranded and floating individuals. Actually, neritic sub-adults and adults are more affected by incidental captures than juveniles because of their feeding behavior. The Bayesian mixed-stock analysis showed a strong genetic relationship between the Adriatic aggregates and Mediterranean rookeries, while a low pro¬portion of individuals of Atlantic origin were detected in the Adriatic feeding grounds. The presence of migratory pathways towards the Adriatic Sea due to the surface current system was reinforced by the finding of individuals bearing haplotypes endemic to the nesting populations of Libya, Greece and Israel. A relatively high contribution from Turkey and Cyprus to the Northwest and South Adriatic populations was identified when the three sampled areas were analyzed independently. These results have to be taken in account in a conservative perspective, since coastal hazards, affecting the population of turtles feeding in the Adriatic Sea may also affect the nesting populations of the Eastern Mediterranean with a unique genetic pattern.

Phylogeny and phylogeography of the family Hyalidae (Crustacea: Amphipoda) along the northeast Atlantic coasts

The family Hyalidae comprises more than one hundred species, distributed worldwide. They are common and abundant in the littoral and shallow sublittoral habitats and they play an important role in the coastal food chain. Most studies about this family were dealing with taxonomy and ecology, while very little is known about phylogenetic relationship among genera and species. In the present study we aim to achieve the first approach of the phylogenetic patterns of this family in NE Atlantic Ocean and Mediterranean Sea, and to perform the first insight into the phylogeography Apohyale prevostii along both the North Atlantic coasts. In order to do that, eight species belonging to the genera Apohyale, Hyale, Serejohyale and Protohyale were investigated using the mitochondrial COI-5P barcode region. Specimens were collected along European and Moroccan Atlantic rocky shores, including Iceland, the British Isles, Macaronesia and in the Mediterranean Sea. Sequences of A. prevostii, from the NW Atlantic Ocean, available in BOLD and GenBank, were retrieved. As expected, phylogenetic analyses showed highly-divergent clades, clearly discriminating among different species clusters, confirming their morphology-based identifications. Although, within A. perieri, A. media, A. stebbingi, P. (Protohyale) schmidtii and S. spinidactylus, high genetic diversity was found, revealing putative cryptic species. The clade of A. prevostii and A. stebbingi appears well supported and divided from the other two congeneric species, and P. (Protohyale) schmidtii shows a basal divergence. The north-western Atlantic coasts were recently colonized by A. prevostii after the last glacial maximum from the European populations showing also a common haplotype in every population analysed. The use of the COI-5P as DNA barcode provided a good tool to underline the necessity of a revision of this emblematic family, as well as to discern taxonomically the possible new species flagged with this molecular device.

Dynamics of inorganic carbon and carbonate saturation in the Gulf of Cádiz (SW Spain)

The study of inorganic carbon chemistry of the coastal ocean is conducted in the Gulf of Cádiz (GoC). Here we describe observations obtained during 4 sampling cruises in March, June, September and November 2015. The primary data set consists of state-of-the-art measurements of the keystone parameters of the marine CO2 system: Total Alkalinity (TA), pH, dissolved inorganic carbon (DIC). We have then calculated aragonite and calcite saturation state. The distribution of inorganic carbon system parameters in the north eastern shelf of the Gulf of Cádiz showed temporal and spatial variability. River input, mixing, primary production, respiration and remineralization were factors that controlled such distributions. Data related to carbonate saturation of calcite and aragonite reveal the occurrence of a supersaturated water; in any case, both species increased with distance and decreased with depth. The carbon system parameters present a different behaviour close to the coast to offshore ad at deeper water. In this area six water masses are clearly identified by their different chemical properties: Surface Atlantic Water, North Atlantic Central Water (NACW) and Mediterranean Water (MOW). Moreover, with this work the measurement of calcium in seawater is optimize, allowing a better quantification for future work of the saturation state of CaCO3.