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.

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.

Elaboración artesanal de mermeladas de tres ecotipos de tuna (Opuntia ficus indica f. inerme) roja, anaranjada y verde

La tuna es una fruta presente principalmente en zonas áridas y semiáridas. Su aporte nutricional podría aprovecharse consumiéndola en forma directa, como así también a través de productos elaborados en forma artesanal. Dichos productos podrían incorporarse a la economía del mercado de consumo interno y eventualmente externo. Por tal motivo, se evalúa la posibilidad de desarrollar mermeladas con tres ecotipos de tuna de la misma variedad, considerando la aceptación del producto. Se trabajó con tunas (Opuntia ficus indica f. inerme), cultivadas en el Centro Regional de Investigaciones Científicas y Tecnológicas (CRICYT) ubicado en el Parque General San Martín, Capital, Mendoza. El objetivo general del trabajo estuvo dirigido a la elaboración de mermeladas de tres ecotipos de tuna: tuna roja o morada, tuna anaranjada o amarilla y tuna verde o blanca acorde a la legislación vigente. Los ítems específicos estuvieron conducidos hacia el desarrollo de una mermelada utilizando una fruta no convencional, para determinar el valor nutricional del mismo, y sus características fisicoquímicas, ampliar la producción de subproductos de tuna en la región árida y semiárida, para así finalmente evaluar el grado de aceptabilidad sensorial y comercial del producto. Se caracterizaron los tres ecotipos de tuna y sus respectivas mermeladas mediante determinaciones fisicoquímicas, se verificó el contenido nutricional y se realizaron los análisis microbiológicos de cada una de ellas. Dentro de las determinaciones fisicoquímicas de las mermeladas se incluyó K, P, Ca y vitamina C, teniendo en cuenta que se parte de una materia prima con buenas cantidades de estos minerales y vitamina C con respecto a otras frutas. Se llevaron a cabo pruebas de evaluación sensorial del producto terminado y se consultó el precio que estarían dispuestos a pagar los potenciales consumidores. Los jueces entrenados contribuyeron al diseño de un perfil sensorial para cada una de las mermeladas y evaluaron en forma global al producto obtenido. Los resultados obtenidos en cuanto a la caracterización de la materia prima arrojan un buen contenido en sólidos solubles, vitamina C y en minerales como el fósforo, calcio, potasio y un contenido bajo en sodio con respecto a otras frutas que se consumen habitualmente. Los análisis microbiológicos realizados para evaluar la conservación permiten expresar que se trata de un producto de bajo riesgo microbiológico, bajo las condiciones de conservación que se obtuvieron. En base a lo expuesto, se puede concluir que es técnicamente posible elaborar artesanalmente mermeladas de tres ecotipos de tuna, las cuales poseen un buen aporte nutricional, son aptas para el consumo humano desde el punto de vista bromatológico y aceptadas sensorialmente por los consumidores.

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1956-1972, compiled from statistics about ICCAT fishery region L1

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1999-2007, compiled from statistics about ICCAT fishery region T11

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Compilation of North Atlantic albacore tuna (Thunnus alalunga) fork length frequencies in 2 centimeter intervals from catches in the North Atlantic for the period 1956-2010

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1987-2009, compiled from statistics about ICCAT fishery region L3

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Compilation of Thunnus thynnus (bluefin tuna) catches in the Northeast Atlantic Ocean 1950-2010

The data have been extracted and compiled from various sources but mainly from the ICES data base. The ICES data are from catch databases downloaded from the ICES website on 2014-01-14. These data are resolved by ICES area, country and year. During inspection of these data, it was noted that Norwegian data for years before 1950 had not been entered into the catch database on the ICES website. ICES has been notified of this omission by B. R. MacKenzie. The Norwegian data from ICES Bulletins. Statistiques has been added. Additional historical bluefin tuna catch data from other fishery reports and sources have been included in the data file for years preceding those when countries started reported their landings officially to ICES. These additional data have been reported in the literature previously (MacKenzie and Myers 2007, Fisheries Research).

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1973-2011, compiled from statistics about ICCAT fishery region L2

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Spatially explicit estimates of North Atlantic albacore tuna (Thunnus alalunga) biomass in the North Atlantic for the period 1956-2010 - Gridded data product (NetCDF)

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species, together with the model estimates achieved from these data, allowing models inter-comparison and evaluation of model skills. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. Length frequencies of catch were also extracted according to the definition of fisheries for the period 1956-2010. Using these data, an application of the spatial ecosystem and population dynamics model (SEAPODYM) was developed for the North Atlantic albacore population and fisheries and provided the first spatially explicit estimate of albacore density in the North Atlantic by life stage. These densities by life stage (larval recruits, young immature fish adult mature fish and total biomass) are provided in gridded file (Netcdf) at resolution of 2° x 2° x month.

Spatially explicit estimates of stock size, structure and biomass of North Atlantic albacore tuna (Thunnus alalunga) in the North Atlantic for the period 1987-2011, compiled from statistics about ICCAT fishery region L7

The development of the ecosystem approach and models for the management of ocean marine resources requires easy access to standard validated datasets of historical catch data for the main exploited species. They are used to measure the impact of biomass removal by fisheries and to evaluate the models skills, while the use of standard dataset facilitates models inter-comparison. North Atlantic albacore tuna is exploited all year round by longline and in summer and autumn by surface fisheries and fishery statistics compiled by the International Commission for the Conservation of Atlantic Tunas (ICCAT). Catch and effort with geographical coordinates at monthly spatial resolution of 1° or 5° squares were extracted for this species with a careful definition of fisheries and data screening. In total, thirteen fisheries were defined for the period 1956-2010, with fishing gears longline, troll, mid-water trawl and bait fishing. However, the spatialized catch effort data available in ICCAT database represent a fraction of the entire total catch. Length frequencies of catch were also extracted according to the definition of fisheries above for the period 1956-2010 with a quarterly temporal resolution and spatial resolutions varying from 1°x 1° to 10°x 20°. The resolution used to measure the fish also varies with size-bins of 1, 2 or 5 cm (Fork Length). The screening of data allowed detecting inconsistencies with a relatively large number of samples larger than 150 cm while all studies on the growth of albacore suggest that fish rarely grow up over 130 cm. Therefore, a threshold value of 130 cm has been arbitrarily fixed and all length frequency data above this value removed from the original data set.

Compilation of Thunnus thynnus (bluefin tuna) catche in the Northeast Atlantic Ocean 1906-1949

The data have been extracted and compiled from various sources but mainly from the ICES data base. The ICES data are from catch databases downloaded from the ICES website on 2014-01-14. These data are resolved by ICES area, country and year. During inspection of these data, it was noted that Norwegian data for years before 1950 had not been entered into the catch database on the ICES website. ICES has been notified of this omission by B. R. MacKenzie. The Norwegian data from ICES Bulletins. Statistiques has been added. Additional historical bluefin tuna catch data from other fishery reports and sources have been included in the data file for years preceding those when countries started reported their landings officially to ICES. These additional data have been reported in the literature previously (MacKenzie and Myers 2007, Fisheries Research).

Stomach content records of albacore (Thunnus alalunga) and bluefin tuna (Thunnus thynnus) in the North Atlantic Drift Region focusing on the Bay of Biscay between 2004 and 2011

The tuna stomach database from AZTI-Tecnalia corresponds to 7 years of sampling from 2004 to 2011. Due to the absence of continuity in the different projects dealing with the feeding ecology of tunas, the sampling could not be performed every year for both species, and no sample was collected in 2008. However, the fish stomach content record contents composition - by prey weight - of 1525 albacore caught in the Bay of Biscay and surrounding waters of the North Atlantic Drift Region in 2005 (n=397), 2006 (n=196), 2007 (n=37), 2009 (n=95), 2010 (n=566) and 2011 (n=234) ; and of 686 bluefin tunas caught in the Southeastern Bay of Biscay in 2004 (n=32), 2005 (n=36), 2006 (n=3), 2009 (n=257), 2010 (n=233) and 2011 (n=125). Samples have been obtained from scientific research surveys (using a variety of different fishing gears), from commercial fisheries catches, from individual fish voluntarily sampled by recreational fishermen and from fish accidentally stranded on coastlines. Each predator is identified by an ID and its length and wet weight are given. In case the wet weight could not be measured, it was estimated through a length-weight relationship equation and is indicated in the comment for the Predator mass column. The total weight of each prey is given, as well as the weight of each prey taxonomic group in each stomach.

Otolith microstructure and trace elemental concentrations in juvenile Pacific bluefin tuna (Thunnus orientalis) caught off of Southern California in 2012

We combined longitudinal analyses of otolith microstructure and trace elemental composition in ~ age 1-2 Pacific bluefin tuna (PBT, n = 24) for inferring the arrival of individuals in the California Current Large Marine Ecosystem (CCLME). Element:Ca ratios in transverse otolith sections (9-12 rows, triplicate ablations from coreprimordium to edge, ø50 µm) were quantified for eight elements: Li, Mg, Mn, Co, Cu, Zn, Sr, and Ba, which was followed by microstructure analysis to provide age estimates corresponding to each ablation spot. Age estimates from otoliths ranged from 328 to 498 days post hatch. The combined elemental signatures of four elements (Ba, Mg, Co, Cu) showed a significant increase at the otolith edge in approximately half of the individuals (30-60 days prior to catch). Given the different oceanographic properties of oligotrophic open Pacific vs. high nutrient, upwelling CCLME waters, this signal is consistent with the entry of the fish into the CCLME, which was estimated to occur primarily in July after a transoceanic migration of ~1.5-2.0 months.