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.

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.

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.

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.

Molecular phylogenesis of Mediterranean Octocorals

In order to support the conservation of the Mediterranean octocorals improvements on information regarding their taxonomic units and phylogenetic relationships are strongly needed. In the present thesis work, phylogenetic analyses based on the mitochondrial mtMSH and 16S genes were performed including 15 Mediterranean octocorals species on the 56 recognized to date. Moreover, an extended datasets with Atlanto/Pacific congeners Octocorallia species was implemented to clarify their phylogenetic relationships and estimate the divergence times of the Mediterranean species. Results indicated that: 1) there are similarity and differences among molecular and morphological traits depending on the taxonomical level considered; 2) the molecular phylogeny of the Mediterranean octocorals retrace the previous relationships based on wide octocorals analyses; and 3) the divergence time among Mediterranean and Atlanto/Pacific species varies depending on analysed taxa. At higher taxonomic level, the Mediterranean trees supported the division of the Mediterranean Octocorallia into one major clade (Alcyoniina-Holaxonia) plus two unresolved branch including the single species available of Scleraxonia and Stolonifera respectively. This topology was better supported including the Atlanto/Pacific congeners species. The molecular evidence suggested that Alcyonium palmatum and Corallium rubrum species are the youngest with a divergence time estimated around 4 MYA. Particularly, C. rubrum results were in agreement with the hypothesis that recent orogenesis process of the Mediterranean Sea promoted the allopatric speciation of this specie. Increasing the sample design and implementing the emerging next-generation genomic-sequencing technologies, further studies would be able to improve the understanding of the Mediterranean octocorals phylogenetic relationships and evolution.

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.

Thermal structure and dynamical modeling of a Mediterranean tropical-like cyclone

In the present work, a detailed analysis of a Mediterranean TLC occurred in January 2014 has been conducted. The author is not aware of other studies regarding this particular event at the publication of this thesis. In order to outline the cyclone evolution, observational data, including weather-stations data, satellite data, radar data and photographic evidence, were collected at first. After having identified the cyclone path and its general features, the GLOBO, BOLAM and MOLOCH NWP models, developed at ISAC-CNR (Bologna), were used to simulate the phenomenon. Particular attention was paid on the Mediterranean phase as well as on the Atlantic phase, since the cyclone showed a well defined precursor up to 3 days before the minimum formation in the Alboran Sea. The Mediterranean phase has been studied using different combinations of GLOBO, BOLAM and MOLOCH models, so as to evaluate the best model chain to simulate this kind of phenomena. The BOLAM and MOLOCH models showed the best performance, by adjusting the path erroneously deviated in the National Centre for Environmental Prediction (NCEP) and ECMWF operational models. The analysis of the cyclone thermal phase shown the presence of a deep-warm core structure in many cases, thus confirming the tropical-like nature of the system. Furthermore, the results showed high sensitivity to initial conditions in the whole lifetime of the cyclone, while the Sea Surface Temperature (SST) modification leads only to small changes in the Adriatic phase. The Atlantic phase has been studied using GLOBO and BOLAM model and with the aid of the same methodology already developed. After tracing the precursor, in the form of a low-pressure system, from the American East Coast to Spain, the thermal phase analysis was conducted. The parameters obtained showed evidence of a deep-cold core asymmetric structure during the whole Atlantic phase, while the first contact with the Mediterranean Sea caused a sudden transition to a shallow-warm core structure. The examination of Potential Vorticity (PV) 3-dimensional structure revealed the presence of a PV streamer that individually formed over Greenland and eventually interacted with the low-pressure system over the Spanish coast, favouring the first phase of the cyclone baroclinic intensification. Finally, the development of an automated system that tracks and studies the thermal phase of Mediterranean cyclones has been encouraged. This could lead to the forecast of potential tropical transition, against with a minimum computational investment.

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.

The nematode parasite fauna in three odontocete species: Phocoena phocoena, Lagenorhynchus albirostris, Lagenorhynchus acutus from Northern Scotland and its relationships with diet composition

Food items and nematode parasites were identified from the stomachs of 42 individuals of Phocoena phocoena, 6 of Lagenorhynchus acutus and 8 of L. albirostris stranded off the coastal waters of Northern Scotland between 2004 and 2014. Post-mortem examinations have revealed heavy parasitic worm burdens. Four nematode species complex as Anisakis spp., Contracaeucum spp., Pseudoterronova spp., and Hysterothylacium spp. were recorded. Data on presence of the anisakid species in cetaceans, reported a significative relationship between the presence of Hysterothylacium and the month of host stranding; suggesting a decrease of larval H. aduncum abundance in the period between April and August due to a seasonal effect related to prey availability. Similarly, the parasite burden of the all anisakid genera was related to the year fraction of stranding, and a relationship statistically significant was found just for L. albirostris with an increase between April and October. This finding is explained by a seasonality in occurrence of white-beaked dolphins, with a peak during August, that might be related to movements of shared prey species and competition with other species (Tursiops truncatus). Geographical differences were observed in parasites number of all anisakid species, which was the highest in cetaceans from the East area and lowest in the North coast. The parasites number also increased significantly with the length of the animal and during the year, but with a significant seasonal pattern only for P. phocoena. Regarding diet composition, through a data set consisting of 34 harbour porpoises and 1 Atlantic white-sided dolphins, we found a positive association between parasite number and the cephalopods genus Alloteuthis. This higher level of parasite infection in squid from this area, is probably due to a quantitative distribution of infective forms in squid prey, an abundance of the final host and age or size maturity of squid.

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.

Spatio-temporal diversity of Megistozoobenthos in the Antalya Gulf and relationships with environmental features

Antalya Gulf is situated in the Levantine Sea, the second biggest and most eastern basin in the Mediterranean Sea. This area is an ultra-oligotrophic basin, strongly affected by anthropogenic inputs, in particular in the fishing areas. For this characteristic, in the Levantine Sea, there is a strong pressure on the natural resources and benthic assemblages. Furthermore, many alien species enter from Suez Canal and are well established in the area. All these pressures are leading to a degradation of the Levantine Sea. For this reason it is important to have tools to study and monitoring the functioning of the marine ecosystem. Benthic organisms are superior to many other biological groups for their response to environmental stresses. The variability of benthic assemblages on a site can reflect, in an integrative mode, the entire functioning of the marine ecosystem. In this study, that wants to analyze the spatial and temporal distribution of the benthic macrofaunal assemblages of Antalya Gulf, 90 benthic species divided in 8 taxa (Annelida, Cnidaria, Echinodermata, Echiura, Mollusca, Porifera, Sipunculida and Tunicata) were found. All the analyses conducted on the entire benthic class and later on Mollusca and Echinodermata separately highlighted the importance of depth on structuring benthic community.

On the Mediterranean conveyor belt system

The Mediterranean Sea is a semi-enclosed sea, connected to the Atlantic Ocean through the Gibraltar Strait and subdivided in two different regions by the Sicily Strait. The geographical extension of the basin, the surface heat flux, and the water inflow from the Gibraltar Strait are some of the basic factors determining its horizontal and vertical circulation. In the Mediterranean strong salinity and temperature zonal gradients contribute to maintain the zonal-vertical circulation, while meridional-vertical cells are equally forced by winds and deep water mass formation in three regions, the Gulf of Lyon, the southern Adriatic and the Cretan Sea areas. The objective of this thesis is to study how these cells combine together to form the Mediterranean conveyor belt system. This has never been attempted before so the conclusions are necessarily preliminary. In the first part we discuss the vertical zonal and meridional circulation by reconstructing the Wust Mediterranean vertical salinity and temperature structures in an attempt to evaluate the water mass structure consistent with modern data. Our results confirm that Wust depicted vertical circulation from scarce data is reproduced by the past 27 years observations. The structure of both meridional and zonal circulations was discussed using velocity vertical streamfunctions with two different methods. The first one, eulerian, allowed us to observe vertical structures that were already reported in the literature. Recent studies in the Atlantic Ocean have shown that gyres and eddies have an important influence in the isopycnal vertical circulation. This is called the residual circulation which was computed in this study for the first time. A possible interpretation of horizontal connection between the meridional and zonal cells was discussed using horizontal streamfunction. In the last part of the thesis we have been developing an idealized numerical model to study the vertical circulation in the Mediterranean.

The routes of the Mediterranean Sea circulation

The Mediterranean Sea is a semi-enclosed basin connected to the Atlantic Ocean through the narrow and shallow Strait of Gibraltar and further subdivided in two different sub-basins, the Eastern Mediterranean and the Western Mediterranean, connected through the Stait of Sicily. On annual basis, a net heat budget of −7 W/m2, combined with exceeding evaporation over precipation and runoff together with wind stress, is responsible for the antiestuarine character of the zonal thermoaline circulation. The outflow at Gibraltar Strait is mainly composed of Levantine Intermediate Water (LIW) and deep water masses formed in the Western Mediterranean Sea. The aim of this thesis is to validate and quantitatively assess the main routes of water masses composing the ouflow at Gibraltar Strait, using for the first time in the Mediterranean Sea a lagrangian interpretation of the eulerian velocity field produced from an eddy-resolving reanalysis dataset, spanning from 2000 to 2012. A lagrangian model named Ariane is used to map out three-dimensional trajectories in order to describe the pathways of water mass transport from the Strait of Sicily, the Gulf of Lyon and the Northern Tyrrhenian Sea to the Gibraltar Strait. Numerical experiments were carried out by seeding millions of particles in the Strait of Gibraltar and following them backwards in time to track the origins of water masses and transport exchanged between the different sections of the Mediterranean. Finally, the main routes of the intermediate and deep water masses are reconstructed from virtual particles trajectories, which highlight the role of the Western Mediterranean Deep Water (WMDW) as the main contributor to the Gibraltar Strait outflow. For the first time, the quantitative description of the flow of water masses coming from the Eastern Mediterranean towards the Gibraltar Strait is provided and a new route that directly links the Northern Tyrrhenian Sea to Gibralatr Strait has been detected.