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

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

Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

Ocean-Atmosphere Interactions: Linking Oceanic Biological Activity to Marine Aerosol Physico-Chemical and Cloud Properties

In this work, in-situ measurements of aerosol chemical composition, particle number size distribution, cloud-relevant properties and ground-based cloud observations were combined with high-resolution satellite sea surface chlorophyll-a concentration and air mass back-trajectory data to investigate the impact of the marine biota on aerosol physico-chemical and cloud properties. Studies were performed over the North-Eastern Atlantic Ocean, the central Mediterranean Sea, and the Arctic Ocean, by deploying both multi-year datasets and short-time scale observations. All the data were chosen to be representative of the marine atmosphere, reducing to a minimum any anthropogenic input. A relationship between the patterns of marine biological activity and the time evolution of marine aerosol properties was observed, under a variety of aspects, from chemical composition to number concentration and size distribution, up to the most cloud‐relevant properties. At short-time scales (1-2 months), the aerosol properties tend to respond to biological activity variations with a delay of about one to three weeks. This delay should be considered in model applications that make use of Chlorophyll-a to predict marine aerosol properties at high temporal resolution. The impact of oceanic biological activity on the microphysical properties of marine stratiform clouds is also evidenced by our analysis, over the Eastern North Atlantic Ocean. Such clouds tend to have a higher number of smaller cloud droplets in periods of high biological activity with respect to quiescent periods. This confirms the possibility of feedback interactions within the biota-aerosol-cloud climate system. Achieving a better characterization of the time and space relationships linking oceanic biological activity to marine aerosol composition and properties may significantly impact our future capability of predicting the chemical composition of the marine atmosphere, potentially contributing to reducing the uncertainty of future climate predictions, through a better understanding of the natural climate system.

Dynamics of terricolous alpine lichen communities of the Alps and Mediterranean Mountains in a climate change perspective

This thesis focuses on the impact of climate change in alpine ecosystems stressing the response of high elevation terricolous lichen communities. In fact, despite the strong sensitivity of cryptogams to changes in climatic factors, information is still scanty.We collected records in 154 plots placed in the summit area of the Majella Massif. In Following a multitaxon approach, Chapter 1 includes cryptogams and vascular plants. We analysed patterns in species richness, beta diversity and functional composition. In Chapter 2, we analysed the relationships between climatic variables and phylogenetic diversity and structure indices. Chapter 3 provides a long-term response relative to the consequences of climate change on a representative terricolous lichen genus across the Alps. Chapter 4 explores the relationships between the species richness and the functional composition of lichen growing on two types of substrates (carbonatic and siliceous soils) along different elevation gradients in the Eastern Alps. Climate change could affect cryptogams and lichens much more than vascular plants in Mediterranean mountains. Contrasting species-climate and traits-climate relationships were found between lichens and bryophytes, suggesting that each group may be sensitive to different components of climate change. Ongoing climate change may also lead to a loss of genetic diversity at high elevation ranges in the Mediterranean mountains, pauperising the life history richness of lichens. Alpine results forecasted that moderate range loss dynamics will occur at low elevation and in peripheral areas of the alpine chain. Results also support the view that range dynamics could be associated with functional traits mainly related to water-use strategies, dispersal, and establishment ability. We also highlighted the importance of substrates as a main driver of both species’ richness and functional traits composition. A “trade-off” also occurs between stress tolerance and the competitive response of communities of terricolous lichens that grow above siliceous and carbonatic soils.

Moving toward stock units identification based on spatial population structure of marine species in the Mediterranean Sea and adjacent waters through multidisciplinary and holistic approaches for the sustainability of fisheries resources

Investigating stock identity of marine species in a multidisciplinary holistic approach can reveal patterns of complex spatial population structure and signatures of potential local adaptation. The population structure of common sole (Solea solea) in the Mediterranean Sea was delineated using genomic and otolith data, including single nucleotide polymorphisms (SNPs) markers and otolith data. SNPs were correlated with environmental and spatial variables to evaluate the impact of these features on the actual genetic population structure. Integrated holistic approach was applied to combine the tracers with different spatio-temporal scales. SNPs data was also used to illustrate the population structure of European hake (Merluccius merluccius) within the Alboran Sea, extending into the neighboring Mediterranean Sea and Atlantic Ocean. The aim was to identify patterns of neutral and potential adaptive genetic variation by applying seascape genomic framework. Results from both genetic and otolith data suggested significant divergence among putative populations of common sole, confirming a clear separation between Western, Adriatic Sea and Eastern Mediterranean Sea. Evidence of fine-scale population structure in the Western Mediterranean Sea was observed at outlier loci level and in the Adriatic. Our study not only indicates that separation among Mediterranean sole population is led primarily by neutral processes, but it also suggests the presence of local adaptation influenced by environmental and spatial factors. The holistic approach by considering the spatio-temporal scales of variation confirmed that the same pattern of separation between these geographical sites is currently occurring and has occurred for many generations. Results showed the occurrence of population structure in Merluccius merluccius by detecting westward–eastward differentiation among populations and distinct subgroups at a fine geographical scale using outlier SNPs. These results enhance the knowledge of the population structure of commercially relevant species to support the application of spatial stock assessment models, including a redefinition of fishery management units.