Genetic variability and structuring of deep Corallium rubrum (Linnaeus, 1758) populations

ABSTRACT Given the decline of shallow-water red coral populations resulting from over-exploitation and mass mortality events, deeper populations below 50 metres depth (mesophotic populations) are currently the most harvested; unfortunately, very little is known about their biology and ecology. The persistence of these populations is tightly linked to their adult density, reproductive success, larval dispersal and recruitment. Moreover, for their conservation, it is paramount understand processes such as connectivity within and among populations. Here, for the first time, genetic variability and structuring of Corallium rubrum populations collected in the Tyrrhenian Sea ranging from 58 to 118 metres were analyzed using ten microsatellite loci and two mitochondrial markers (mtMSH and MtC). The aims of the work were 1) to examine patterns of genetic diversity within each geographic area (Elba, Ischia and Praiano) and 2) to define population structuring at different spatial scales (from tens of metres to hundreds of kilometres). Based on microsatellite data set, significant deviations from Hardy-Weinberg equilibrium due to elevated heterozygote deficiencies were detected in all samples, probably related to the presence of null alleles and/or inbreeding, as was previously observed in shallow-water populations. Moreover, significant levels of genetic differentiation were observed at all spatial scale, suggesting a recent isolation of populations. Biological factors which act at small spatial scale and/or abiotic factors at larger scale (e.g. summer gyres or absence of suitable substrata for settlement) could determine this genetic isolation. Using mitochondrial markers, significant differences were found only at wider scale (between Tuscany and Campania regions). These results could be related to the different mutation rate of the molecular makers or to the occurrence of some historical links within regions. A significant isolation by distance pattern was then observed using both data sets, confirming the restricted larval dispersal capability of the species. Therefore, the hypothesis that deeper populations may act as a source of larvae helping recovery of threatened shallow-water populations is not proved. Conservation strategies have to take into account these results, and management plans of deep and currently harvested populations have to be defined at a regional or sub regional level, similarly to shallow-water populations. Nevertheless, further investigations should be needed to understand better the genetic structuring of this species in the mesophotic zone, e.g. extending studies to other Mediterranean deep-water populations.

Genetic structure and connectivity between populations of two common Mediterranean sessile invertebrates

Population genetic and phylogeography of two common mediterranean species were studied in 10 localities located on the coasts of Toscana, Puglia and Calabria. The aim of the study was to verify the extent of genetic breaks, in areas recognized as boundaries between Mediterranean biogeographic sectors. From about 100 sequences obtained from the mitochondrial Cytochrome Oxidase subunit I (COI) gene of Halocynthia papillosa and Hexaplex trunculus genetic diversity, genetic structure at small and large distances and demographic history of both specieswere analyzed. No evidences of genetic breaks were found for the two species in Toscana and Puglia. The genetic structure of H. trunculus evidences the extent of a barrier to gene flow localized in Calabria, which could be represented by the Siculo-Tunisian Strait and the Strait of Messina. The observed patterns showed similar level of gene flow at small distances in both species, although the two species have different larval ecology. These results suggest that other factors, such as currents, local dynamics and seasonal temperatures, influence the connectivity along the Italian peninsula. The geographic distribution of the haplotypes shows that H. papillosacould represent a single genetic pool in expansion, whereas H. trunculus has two distinct genetic pools in expansion. The demographic pattern of the two species suggests that Pleistocene sea level oscillations, in particular of the LGM, may have played a key role in shaping genetic structure of the two species. This knowledge provides basic information, useful for the definition of management plans, or for the design of a network of marine protected areas along the Italian peninsula.

Genetic structuring of Eunicella singularis populations along a bathymetric gradient

Phenotypic plasticity refers to the ability of an organism to express different morphologies depending on the abiotic and biotic environment. Depth integrating many variables (e.g. temperature, light and hydrodynamics), may affect population structure and dynamics of the populations, as well as connectivity patterns and genetic diversity. Eunicella singularis is a Mediterranean arborescent gorgonian who plays an important rule as engineer species providing biomass and complexity to coralligenous habitats. It has a wide bathymetric distribution ranging from shallow rocky bottoms to deep sublittoral reefs. The species shows two depth-related morphotypes which taxonomic status is not yet clarified. The aim of the study is to analyses genetic variability and/or structuring along a vertical gradient to test the presence of the two morphotypes. Furthermore, a preliminary analyses of the phylogenetic relationship among species of the genus Eunicella has been done. Six populations of Eunicella singularis were sampled from 10 to 60 m depth in Cap de Creus and individuals belonging to Eunicella cavolinii, E. verrucosa, E. racemosa and E. stricta aphyta were collected. The genetic analyses were carried out using five microsatellite loci and ITS-1 sequence polymorphism. The results showed a reduction of genetic variability along the vertical gradient. A threshold in connectivity was observed across 30 - 40 m depth, confirming the presence of two different Eunicella singularis morphotypes. The two morphological forms could be due to phenotypic plasticity, which allowed populations to suit different environmental conditions, or to a break in gene flow that determined the isolation of the two populations and an accumulation of genetic differences. The molecular markers used were not able to clarify the phylogenetic relationship among Eunicella species and the systematic position of the two morphotypes, conversely they risen the question on the existence of single species of Mediterranean Eunicella.

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

Are there genetic breaks between Atlantic and Pacific Yellowfin tuna (Thunnus albacares) populations? A preliminary study based on microsatellites gene variation

Yellowfin tuna (Thunnus albacares, YFT, Bonnaterre 1788) is one of the most important market tuna species in the world. The high mortality of juveniles is in part caused by their bycatch. Indeed, if unregulated, it could permanently destabilize stocks health. For this reason investigating and better knowing the stock boundaries represent a crucial concern. Aim of this thesis was to preliminary investigate the YFT population structure within and between Atlantic and Pacific Oceans through the analysis of genetic variation at eight microsatellite loci and assess the occurrence of barriers to the gene flow between Oceans. For this propouse we collected 4 geographical samples coming from Atlantic and Pacific Ocean and selected a panel of 8 microsatellites loci developped by Antoni et al., (2014). Samples 71-2-Y and 77-2-Y, came from rispectively west central pacific ocean (WCPO) and east central pacific ocean (ECPO), instead samples 41-1-Y and 34-2-Y derive from west central atlantic ocean (WCAO) and east central atlantic ocean (ECAO). Total 160 specimens were analyzed (40 per sample) and were carried out several genetic information as allele frequencies, allele number, allelic richness, HWE (using He and Ho) and pairwise Fst genetic distance. Results obtained, may support the panmictic theory of this species, only one of pairwise Fst obtained is statistically significant (Fst= 0.00927; pV= 0.00218) between 41-1-Y and 71-2-Y samples. Results suggest low genetic differentiation and consequent high level of gene flow between Atlantic and Pacific populations. Furthermore, we performed an analysis of molecular taxonomy through the use of ATCO (the flaking region between ATPse6 and cytochrome oxidase subunit III genes mt DNA, to discriminate within the gener Thunnus two of the related species (Yellofin and bigeye tuna) according with their difficult recognition at certain size (<40 cm). ATCO analysis in this thesis, has provided strong discriminate evidence between the target species proving to be one of the most reliable genetic tools capable to indagate within the genus Thunnus. Thus, our study has provided useful information for possible use of this protocol for conservation plans and management of this fish stocks.

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.

Morphological-molecular taxonomy and population genetic structure of common and Egyptian soles of the Adriatic Sea

A total of 352 specimens were analyzed to achieve the different aims of this thesis. 255 central-northern Adriatic specimens of S. solea and S. aegyptiaca were molecularly analysed using microsatellite locus Sos(AC)40 and 205 also morphologically due to evaluate the abundance and the distribution of the cryptic species S. aegyptiaca and to confirm morphologic analyses. Morphological and molecular analyses comparated show a correspondence of 96%. A combined morphologic approach could be proposed to apply multiple criteria on the analyzed external morphological keys. The Adriatic Egyptian soles may lives in shallow waters (up 30 m) and in brackish lagoon. 127 samples of Adriatic common sole added to 326 samples of previous studies showed, using mitochondrial marker (CytB), that the Adriatic Sea as contact zone between Tyrrhenian and Aegean Sea, the divergence within the Adriatic Sea is low but significant between central-north and south, with a longitudinal strong gene flow in central-northern side. It’s also showed as in the Adriatic Sea two near-panmictic populations of common sole exist.

Genetic variability and differentiation in the Mediterranean endemic starry ray Raja asterias (Delaroche, 1809)

The present study deal with the population structure and connectivity of the Mediterranean endemic starry ray Raja asterias (Delaroche, 1809) in the Western and Eastern Mediterranean basin. A panel of eight microsatellite loci which cross-amplify in Rajidae (El Nagar, 2010) was used to assess population connectivity and structure. Those aims were investigated by analyzing the genetic variation of 9 population sample for a total of 185 individuals collected during past scientific surveys (MEDITS, GRUND), commercial trawling and also directly at fish markets. The purpose of this thesis is to estimate the genetic divergence occurring between the Mediterranean populations and, in particular, to assess the presence of any barrier (geographic, hydrogeological and biological) to gene flow for this species. Different statistical approaches were performed to reach this aim evaluating both the genetic diversity (nucleotide diversity, allelic richness, observed and expected heterozygosity and Hardy-Weinberg equilibrium test) and the population differentiation patterns (pairwise Fst estimated and population structure analysis). The results obtained from the analysis of the microsatellite dataset suggest a geographic and genetic separation between the starry ray populations of the Mediterranean basin into three or four distinct groups: Western and Eastern Mediterranean basins and Sicilian coast always clustering as an independent group and Algeria which could be or not considered another separate group. The data were discussed from both an evolutionary and a conservation point of view and in relation to previous results obtained by the analysis of mitochondrial marker. A comparison with other Mediterranean demersal skate species was performed in order to better contextualise our results. Finally, our results could offer useful information to protect vulnerable species as R. asterias and developing effective conservation plans in the Mediterranean.

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.

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.

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.