Species identification, distributioin and reproductive interactions of common smoot-hound and blackspotted smooth-hound shark (Genus mustelus) in the Adriatic Sea

Longstanding taxonomic ambiguity and uncertainty exist in the identification of the common (M. mustelus) and blackspotted (M. punctulatus) smooth-hound in the Adriatic Sea. The lack of a clear and accurate method of morphological identification, leading to frequent misidentification, prevents the collation of species-specific landings and survey data for these fishes and hampers the delineation of the distribution ranges and stock boundaries of the species. In this context, adequate species-specific conservation and management strategies can not be applied without risks of population declining and local extinction. In this thesis work I investigated the molecular ecology of the two smooth-hound sharks which are abundant in the demersal trawl surveys carried out in the NC Adriatic Sea to monitor and assess the fishery resources. Ecological and evolutionary relationships were assessed by two molecular tests: a DNA barcoding analysis to improve species identification (and consequently the knowledge of their spatial ecology and taxonomy) and a hybridization assay based on the nuclear codominant marker ITS2 to evaluate reproductive interactions (hybridization or gene introgression). The smooth-hound sharks (N=208) were collected during the MEDITS 2008 and 2010 campaigns along the Italian and Croatian coasts of the Adriatic Sea, in the Sicilian Channel and in the Algerian fisheries. Since the identification based on morphological characters is not strongly reliable, I performed a molecular identification of the specimens producing for each one the cytochrome oxidase subunit 1 (COI) gene sequence (ca. 640 bp long) and compared them with reference sequences from different databases (GenBank and BOLD). From these molecular ID data I inferred the distribution of the two target species in the NC Adriatic Sea. In almost the totality of the MEDITS hauls I found no evidence of species sympatry. The data collected during the MEDITS survey showed an almost different distribution of M. mustelus (confined along the Italian coasts) and M. punctulatus (confined along the Croatian coasts); just one sample (Gulf of Venice, where probably the ranges of the species overlap) was found to have catches of both the species. Despite these data results suggested no interaction occurred between my two target species at least during the summertime (the period in which MEDITS survey is carried out), I still wanted to know if there were inter-species reproductive interactions so I developed a simple molecular genetic method to detect hybridization. This method is based on DNA sequence polymorphism among species in the nuclear ribosomal Internal Transcribed Spacer 2 locus (ITS2). Its application to the 208 specimens collected raised important questions regarding the ecology of this two species in the Adriatic Sea. In fact results showed signs of hybridization and/or gene introgression in two sharks collected during the trawl survey of 2008 and one collected during the 2010 one along the Italian and Croatian coasts. In the case that it will be confirmed the hybrid nature of these individuals, a spatiotemporal overlapping of the mating behaviour and ecology must occur. At the spatial level, the northern part of the Adriatic Sea (an area where the two species occur with high frequency of immature individuals) could likely play the role of a common nursery area for both species.

Spatial and temporal variability and ecological processes in the epibenthic assemblages of the northern Adriatic Sea

Several coralligenous reefs occur in the soft bottoms of the northern Adriatic continental shelf. Mediterranean coralligenous habitats are characterised by high species diversity and are intrinsically valuable for their biological diversity and for the ecological processes they support. The conservation and management of these habitats require quantifying spatial and temporal variability of their benthic assemblages. This PhD thesis aims to give a relevant contribution to the knowledge of the structure and dynamics of the epibenthic assemblages on the coralligenous subtidal reefs occurring in the northern Adriatic Sea. The epibenthic assemblages showed a spatial variation larger compared to temporal changes, with a temporal persistence of reef-forming organisms. Assemblages spatial heterogeneity has been related to morphological features and geographical location of the reefs, together with variation in the hydrological conditions. Manipulative experiments help to understand the ecological processes structuring the benthic assemblages and maintaining their diversity. In this regards a short and long term experiment on colonization patterns of artificial substrata over a 3-year period has been performed in three reefs, corresponding to the three main types of assemblages detected in the previous study. The first colonisers, largely depending by the different larval supply, played a key role in determining the heterogeneity of the assemblages in the early stage of colonisation. Lateral invasion, from the surrounding assemblages, was the driver in structuring the mature assemblages. These complex colonisation dynamics explained the high heterogeneity of the assemblages dwelling on the northern Adriatic biogenic reefs. The buildup of these coralligenous reefs mainly depends by the bioconstruction-erosion processes that has been analysed through a field experiment. Bioconstruction, largely due to serpulid polychaetes, prevailed on erosion processes and occurred at similar rates in all sites. Similarly, the total energy contents in the benthic communities do not differ among sites, despite being provided by different species. Therefore, we can hypothesise that both bioconstruction processes and energetic storage may be limited by the availability of resources. Finally the major contribution of the zoobenthos compared to the phytobenthos to the total energetic content of assemblages suggests that the energy flow in these benthic habitats is primarily supported by planktonic food web trough the filter feeding invertebrates.

The stratigraphic record of the quaternary sea level fluctuations and the impact of the post-glacial sea level rise (Termination I) in the Adriatic basin (Mediterranean sea)

The modern stratigraphy of clastic continental margins is the result of the interaction between several geological processes acting on different time scales, among which sea level oscillations, sediment supply fluctuations and local tectonics are the main mechanisms. During the past three years my PhD was focused on understanding the impact of each of these process in the deposition of the central and northern Adriatic sedimentary successions, with the aim of reconstructing and quantifying the Late Quaternary eustatic fluctuations. In the last few decades, several Authors tried to quantify past eustatic fluctuations through the analysis of direct sea level indicators, among which drowned barrier-island deposits or coral reefs, or indirect methods, such as Oxygen isotope ratios (δ18O) or modeling simulations. Sea level curves, obtained from direct sea level indicators, record a composite signal, formed by the contribution of the global eustatic change and regional factors, as tectonic processes or glacial-isostatic rebound effects: the eustatic signal has to be obtained by removing the contribution of these other mechanisms. To obtain the most realistic sea level reconstructions it is important to quantify the tectonic regime of the central Adriatic margin. This result has been achieved integrating a numerical approach with the analysis of high-resolution seismic profiles. In detail, the subsidence trend obtained from the geohistory analysis and the backstripping of the borehole PRAD1.2 (the borehole PRAD1.2 is a 71 m continuous borehole drilled in -185 m of water depth, south of the Mid Adriatic Deep - MAD - during the European Project PROMESS 1, Profile Across Mediterranean Sedimentary Systems, Part 1), has been confirmed by the analysis of lowstand paleoshorelines and by benthic foraminifera associations investigated through the borehole. This work showed an evolution from inner-shelf environment, during Marine Isotopic Stage (MIS) 10, to upper-slope conditions, during MIS 2. Once the tectonic regime of the central Adriatic margin has been constrained, it is possible to investigate the impact of sea level and sediment supply fluctuations on the deposition of the Late Pleistocene-Holocene transgressive deposits. The Adriatic transgressive record (TST - Transgressive Systems Tract) is formed by three correlative sedimentary bodies, deposited in less then 14 kyr since the Last Glacial Maximum (LGM); in particular: along the central Adriatic shelf and in the adjacent slope basin the TST is formed by marine units, while along the northern Adriatic shelf the TST is represented by costal deposits in a backstepping configuration. The central Adriatic margin, characterized by a thick transgressive sedimentary succession, is the ideal site to investigate the impact of late Pleistocene climatic and eustatic fluctuations, among which Meltwater Pulses 1A and 1B and the Younger Dryas cold event. The central Adriatic TST is formed by a tripartite deposit bounded by two regional unconformities. In particular, the middle TST unit includes two prograding wedges, deposited in the interval between the two Meltwater Pulse events, as highlighted by several 14C age estimates, and likely recorded the Younger Dryas cold interval. Modeling simulations, obtained with the two coupled models HydroTrend 3.0 and 2D-Sedflux 1.0C (developed by the Community Surface Dynamics Modeling System - CSDMS), integrated by the analysis of high resolution seismic profiles and core samples, indicate that: 1 - the prograding middle TST unit, deposited during the Younger Dryas, was formed as a consequence of an increase in sediment flux, likely connected to a decline in vegetation cover in the catchment area due to the establishment of sub glacial arid conditions; 2 - the two-stage prograding geometry was the consequence of a sea level still-stand (or possibly a fall) during the Younger Dryas event. The northern Adriatic margin, characterized by a broad and gentle shelf (350 km wide with a low angle plunge of 0.02° to the SE), is the ideal site to quantify the timing of each steps of the post LGM sea level rise. The modern shelf is characterized by sandy deposits of barrier-island systems in a backstepping configuration, showing younger ages at progressively shallower depths, which recorded the step-wise nature of the last sea level rise. The age-depth model, obtained by dated samples of basal peat layers, is in good agreement with previous published sea level curves, and highlights the post-glacial eustatic trend. The interval corresponding to the Younger Dyas cold reversal, instead, is more complex: two coeval coastal deposits characterize the northern Adriatic shelf at very different water depths. Several explanations and different models can be attempted to explain this conundrum, but the problem remains still unsolved.

Effects of temperature and acidity on Mediterranean benthic calcifiers, in the face of globale climate change

The Mediterranean Sea is expected to react faster to global change compared to the ocean and is already showing more pronounced warming and acidification rates. A study performed along the Italian western coast showed that porosity of the skeleton increases with temperature in the zooxanthellate (i.e. symbiotic with unicellular algae named zooxanthellae) solitary scleractinian Balanophyllia europaea while it does not vary with temperature in the solitary non-zooxanthellate Leptopsammia pruvoti. These results were confirmed by another study that indicated that the increase in porosity was accompanied by an increase of the fraction of the largest pores in the pore-space, perhaps due to an inhibition of the photosynthetic process at elevated temperatures, causing an attenuation of calcification. B. europaea, L. pruvoti and the colonial non-zooxanthellate Astroides calycularis, transplanted along a natural pH gradient, showed that high temperature exacerbated the negative effect of lowered pH on their mortality rates. The growth of the zooxanthellate species did not react to reduced pH, while the growth of the two non-zooxanthellate species was negatively affected. Reduced abundance of naturally occurring B. europaea, a mollusk, a calcifying and a non-calcifying macroalgae were observed along the gradient while no variation was seen in the abundance of a calcifying green alga. With decreasing pH, the mineralogy of the coral and mollusk did not change, while the two calcifying algae decreased the content of aragonite in favor of the less soluble calcium sulphates and whewellite (calcium oxalate), possibly as a mechanism of phenotypic plasticity. Increased values of porosity and macroporosity with CO2 were observed in B. europaea specimens, indicating reduces the resistance of its skeletons to mechanical stresses with increasing acidity. These findings, added to the negative effect of temperature on various biological parameters, generate concern on the sensitivity of this zooxanthellate species to the envisaged global climate change scenarios.

Demersal communities in the Mediterranean Sea: a case study of Triglidae (Osteichthyes, Scorpaeniformes) on the conservation and sustainable use of marine resources

An appropriate management of fisheries resources can only be achieved with the continuous supply of information on the structure and biology of populations, in order to predict the temporal fluctuations. This study supports the importance of investigating the bio-ecology of increasingly exploited and poorly known species, such as gurnards (Osteichthyes, Triglidae) from Adriatic Sea (Mediterranean), to quantify their ecological role into marine community. It also focuses on investigate inter and intra-specific structuring factor of Adriatic population. These objectives were achieved by: 1) investigating aspects of the population dynamics; 2) studying the feeding biology through the examination of stomach contents; 3) using sagittal otoliths as potential marker of species life cycle; 4) getting preliminary data on mDNA phylogeny. Gurnards showed a specie-specific “critical size” coinciding with the start of sexual maturity, the tendency to migrate to greater depths, a change of diet from crustaceans to fish and an increase of variety of food items eaten. Distribution of prey items, predator size range and depth distribution were the main dimensions that influence the breadth of trophic niche and the relative difference amongst Adriatic gurnards. Several feeding preferences were individuated and a possible impact among bigger-size gurnards and other commercial fishes (anchovy, gadoids) and Crustacea (such as mantis prawn and shrimps) were to be necessary considered. Otolith studies showed that gurnard species have a very fast growth despite other results in other areas; intra-specific differences and the increase in the variability of otolith shape, sulcus acusticus shape, S:O ratios, sulcus acusticus external crystals arrangement were shown between juveniles and adults and were linked to growth (individual genetic factors) and to environmental conditions (e.g. depth and trophic niche distribution). In order to facilitate correct biological interpretation of data, molecular data were obtained for comparing morphological distance to genetic ones.

Species identification of archived fish bones collected from the Mediterranean Sea 100 years ago using molecular techniques

With the discovery that DNA can be successfully recovered from museum collections, a new source of genetic information has been provided to extend our comprehension of the evolutionary history of species. However, historical specimens are often mislabeled or report incorrect information of origin, thus accurate identification of specimens is essential. Due to the highly damaged nature of ancient DNA many pitfalls exist and particular precautions need to be considered in order to perform genetic analysis. In this study we analyze 208 historical remains of pelagic fishes collected in the beginning of the 20th century. Through the adaptation of existing protocols, usually applied to human remains, we manage to successfully retrieve valuable genetic material from almost all of the examined samples using a guanidine and silica column-based approach. The combined use of two mitochondrial markers cytochrome-oxidase-1(mtDNA COI) and Control Region (mtDNA CR), and the nuclear marker first internal transcriber space (ITS1) allowed us to identify the majority of the examined specimens using traditional PCR and Sanger sequencing techniques. The creation of primers capable of amplifying heavily degraded DNA have great potential for future uses, both in ancient and in modern investigation. The methodologies developed in this study can in fact be applied for other ancient fish specimens as well as cooked or canned samples.

Assessment of marine litter in the Western Adriatic Sea and its impact on marine biota

Due to its environmental, safety, health and socio-economic impacts, marine litter has been recognized as a 21st century global challenge, so that it has been included in Descriptor 10 of the EU MSFD. For its morphological features and anthropogenic pressures, the Adriatic Sea is very sensitive to the accumulation of debris, but data are inconsistent and fragmented. This thesis, in the framework of DeFishGear project, intents to assess marine litter on beaches and on seafloor in the Western Adriatic sea, and test if debris ingestion by fish occurs. Three beaches were sampled during two surveys in 2015. Benthic litter monitoring was carried out in the FAO GSA17 during fall 2014, using a rapido trawl. Litter ingestion was investigated through gut contents analysis of 260 fish belonging to 8 commercial species collected in Western Gulf of Venice. Average litter density on beaches was 1.5 items/m2 during spring, and decreased to 0.8 items/m2 in summer. Litter composition was heterogeneous and varied among sites, even if no significant differences were found. Most of debris consisted of plastic sheets, fragments, polystyrene pieces, mussels nets and cottons bud sticks, showing that sources are many and include aquaculture, land-based activities and local users of beaches. Average density of benthic litter was 913 items/Km2 (82 Kg/Km2). Plastic dominated in terms of numbers and weight, and consisted mainly of bags, sheets and mussel nets. The highest density was found close to the coast, and sources driving the major differences in litter distribution were mussel farms and shipping lanes. Litter ingestion occurred in 47% of examined fish, mainly consisting of fibers. Among species, S. pilchardus swallowed almost all debris categories. Findinds may provide a baseline to set the necessary measures to manage and minimize marine litter in the Western Adriatic region and to protect aquatic life from plastic pollution, even accounting the possible implications on human health.

Temporal dynamics of benthic assemblages along a gradient of ocean acidification at a CO2 vent’s system of the Ischia Island

Global warming and ocean acidification, due to rising atmospheric levels of CO2, represent an actual threat to terrestrial and marine environments. Since Industrial Revolution, in less of 250 years, pH of surface seawater decreased on average of 0.1 unit, and is expected to further decreases of approximately 0.3-0.4 units by the end of this century. Naturally acidified marine areas, such as CO2 vent systems at the Ischia Island, allow to study acclimatation and adaptation of individual species as well as the structure of communities, and ecosystems to OA. The main aim of this thesis was to study how hard bottom sublittoral benthic assemblages changed trough time along a pH gradient. For this purpose, the temporal dynamics of mature assemblages established on artificial substrates (volcanic tiles) over a 3 year- period were analysed. Our results revealed how composition and dynamics of the community were altered and highly simplified at different level of seawater acidification. In fact, extreme low values of pH (approximately 6.9), affected strongly the assemblages, reducing diversity both in terms of taxa and functional groups, respect to lower acidification levels (mean pH 7.8) and ambient conditions (8.1 unit). Temporal variation was observed in terms of species composition but not in functional groups. Variability was related to species belonging to the same functional group, suggesting the occurrence of functional redundancy. Therefore, the analysis of functional groups kept information on the structure, but lost information on species diversity and dynamics. Decreasing in ocean pH is only one of many future global changes that will occur at the end of this century (increase of ocean temperature, sea level rise, eutrophication etc.). The interaction between these factors and OA could exacerbate the community and ecosystem effects showed by this thesis.

Ecology of larval fishes and large zooplankton in the Keweenaw Current region of Lake Superior, with special focus on lake herring, Coregonus artedi

I assessed the influence of the Keweenaw Current and spring thermal bar on the distribution of larval fishes and large zooplankton in Lake Superior. In 1998 and 1999, samples were collected from inshore (0.2 – 3.0 km from shore) and offshore (5.0 – 9.0 km from shore) locations on three transects off the western coast of the Keweenaw Peninsula, Michigan. For larval fishes, density and size distribution patterns of lake herring (Coregonus artedi), rainbow smelt (Osmerus mordax), burbot (Lota lota), deepwater sculpin (Myoxocephalus thompsoni), and spoonhead sculpin (Cottus ricei) suggest a seasonal inshore to offshore movement. For zooplankton, seasonal warming appeared to be the major factor that limited planktonic catches of the primarily benthic Mysisrelicta and Diporeia spp., while simultaneously stimulated growth and reproduction of the cladocerans Daphnia spp., Holopedium gibberum, and Bythotrephes cederstroemi. In contrast, calanoid copepods as a group were abundant throughout the entire sampling season. The greatest abundances of zooplankton were generally encountered offshore, even for the cladocerans, which apparently expanded from inshore to offshore locations with seasonal warming. In 2000, sampling efforts focused on lake herring. Samples were collected from surface waters at 0.1 – 17.0 km from shore on two transects. Lake herring larvae were also reared in the laboratory from eggs in order to validate the use of otolith microstructure for aging. Increment deposition was not statistically different from a daily rate starting from 28 days after hatching, near the time of yolk-sac absorption, but larvae with lower growth rates could not be aged as accurately. In Lake Superior, lake herring tended to be slightly more abundant, larger, and older at inshore locations, but a dense patch of younger larvae was also encountered 7 – 13 km from shore. The distribution iiipatterns suggest that larvae were transported by prevailing currents into the study region, possibly from the more productive spawning regions in western Lake Superior. Growth rates were suppressed at offshore locations where temperatures were less than 8°C. These results indicate that lake herring larvae may be transported far distances from spawning concentrations by longshore currents, and water temperatures may largely control their growth.

Investigation of Quaternary diatoms in sediment core CRP-1 from the Ross Sea, Antarctica

In the first season of drilling, the Cape Roberts Project (CRP) recovered one drillcore (CRP-l) from Roberts Ridge in western McMurdo Sound, Ross Sea, Antarctica Diatom biostratigraphy places the upper six lithostratigraphic units (Units 1.1, 2.1, 2.2, 2.3, 3.1, and 4.1) of CRP-l (0.0 to 43.15 mbsf) within the Quaternary. Both non-marine and marine Quaternary diatoms occur in variable abundance in the Quaternary interval of CRP- 1 Biostratigraphic data resolve two Quaternary time slices or events within CRP-1. Marine diatom assemblages in Units 4.1 and 3.1 represent sedimentation within the diatom Actinocyclus ingens Zone (1.35 to 0.66 Ma). Further refinement of the age of Unit 3.l places deposition in the interval 1.15 to 0.75 Ma based on the common occurrence of Thalassiosira elliptipora and correlation to the Southern Ocean acme of this taxon The absence of ActiActinocyclus ingens and the presence ot Thalassiosira antarctica in Unit 2.2 require a younger zonal assignment for this interval, within the diatom Thalassiosira lentiginosa Zone (0.66 to 0.0 Ma). A new diatom species. Rouxia leventerae, is described from marine assemblages of Units 2.2, 2.3, 3.1, and 4.l. Lithostratigraphic Unit 3.1 (33.82 to 31.89 mbsf) is a bryozoan-dominated skeletal-carbonate facies. Low abundance of Fragilariopsis curta and Fragilariopsis cylindrus within this unit combined with the relatively high abundance of species associated with open water indicates deposition in waters that remained ice free for much or all of the year Diatom assemblages suggest carbonate deposition in Unit 3.1 is linked to a significant early Pleistocene event in McMurdo Sound, when elevated surface-water temperatures inhibited the formation of sea ice.

Diversity and faunal composition of benthic foraminifera in Hole 44-390A on the Blake Nose Plateau

Benthic foraminiferal assemblages are a widespread tool to understand changes in organic matter flux and bottom-water oxygenation and their relation to paleoceanographic changes in the Upper Cretaceous oceans. In this study, assemblage data (diversity, total number, and number per species and gram) from Deep Sea Drilling Project (DSDP) Site 390 (Blake Nose, western North Atlantic) were processed for the lower Maastrichtian (Globotruncana falsostuarti - Gansserina gansseri Planktic Foraminiferal Zone). These data document significant changes in nutrient flux to the sea floor as well as bottom-water oxygenation during this time interval. Parallel to the observed changes in the benthic foraminiferal assemblages the number of inoceramid shells decreases, reflecting also a significant increase in bottom-water oxygenation. We speculate, that these data could reflect the onset of a shift from warmer low-latitude to cooler high-latitude deep-water sources. This speculation will predate the major reorganization of the oceanic circulation resulting in a circulation mode similar to today at the Early/Late Maastrichtian boundary by ~1 Ma and therefore improves our understanding of Late Cretaceous paleoceanography.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.