Fausto Tinti

Identifying fishes through DNA barcodes and microarrays.

Background International fish trade reached an import value of 62.8 billion Euro in 2006, of which 44.6% are covered by the European Union. Species identification is a key problem throughout the life cycle of fishes: from eggs and larvae to adults in fisheries research and control, as well as processed fish products in consumer protection. Methodology/Principal Findings This study aims to evaluate the applicability of the three mitochondrial genes 16S rRNA (16S), cytochrome b (cyt b), and cytochrome oxidase subunit I (COI) for the identification of 50 European marine fish species by combining techniques of “DNA barcoding” and microarrays. In a DNA barcoding approach, neighbour Joining (NJ) phylogenetic trees of 369 16S, 212 cyt b, and 447 COI sequences indicated that cyt b and COI are suitable for unambiguous identification, whereas 16S failed to discriminate closely related flatfish and gurnard species. In course of probe design for DNA microarray development, each of the markers yielded a high number of potentially species-specific probes in silico, although many of them were rejected based on microarray hybridisation experiments. None of the markers provided probes to discriminate the sibling flatfish and gurnard species. However, since 16S-probes were less negatively influenced by the “position of label” effect and showed the lowest rejection rate and the highest mean signal intensity, 16S is more suitable for DNA microarray probe design than cty b and COI. The large portion of rejected COI-probes after hybridisation experiments (>90%) renders the DNA barcoding marker as rather unsuitable for this high-throughput technology. Conclusions/Significance Based on these data, a DNA microarray containing 64 functional oligonucleotide probes for the identification of 30 out of the 50 fish species investigated was developed. It represents the next step towards an automated and easy-to-handle method to identify fish, ichthyoplankton, and fish products.

Outlier SNP markers reveal fine-scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (FCT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (FCT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.

DNA Microarrays for Identifying Fishes

In many cases marine organisms and especially their diverse developmental stages are difficult to identify by morphological characters. DNA-based identification methods offer an analytically powerful addition or even an alternative. In this study, a DNA microarray has been developed to be able to investigate its potential as a tool for the identification of fish species from European seas based on mitochondrial 16S rDNA sequences. Eleven commercially important fish species were selected for a first prototype. Oligonucleotide probes were designed based on the 16S rDNA sequences obtained from 230 individuals of 27 fish species. In addition, more than 1200 sequences of 380 species served as sequence background against which the specificity of the probes was tested in silico. Single target hybridisations with Cy5-labelled, PCR-amplified 16S rDNA fragments from each of the 11 species on microarrays containing the complete set of probes confirmed their suitability. True-positive, fluorescence signals obtained were at least one order of magnitude stronger than false-positive cross-hybridisations. Single nontarget hybridisations resulted in cross-hybridisation signals at approximately 27% of the cases tested, but all of them were at least one order of magnitude lower than true-positive signals. This study demonstrates that the 16S rDNA gene is suitable for designing oligonucleotide probes, which can be used to differentiate 11 fish species. These data are a solid basis for the second step to create a “Fish Chip” for approximately 50 fish species relevant in marine environmental and fisheries research, as well as control of fisheries products.

Spatio-temporal population structuring and genetic diversity retention in depleted Atlantic Bluefin tuna of the Mediterranean Sea

Fishery genetics have greatly changed our understanding of population dynamics and structuring in marine fish. In this study, we show that the Atlantic Bluefin tuna (ABFT, Thunnus thynnus), an oceanic predatory species exhibiting highly migratory behavior, large population size, and high potential for dispersal during early life stages, displays significant genetic differences over space and time, both at the fine and large scales of variation. We compared microsatellite variation of contemporary (n = 256) and historical (n = 99) biological samples of ABFTs of the central-western Mediterranean Sea, the latter dating back to the early 20th century. Measures of genetic differentiation and a general heterozygote deficit suggest that differences exist among population samples, both now and 96–80 years ago. Thus, ABFTs do not represent a single panmictic population in the Mediterranean Sea. Statistics designed to infer changes in population size, both from current and past genetic variation, suggest that some Mediterranean ABFT populations, although still not severely reduced in their genetic potential, might have suffered from demographic declines. The short-term estimates of effective population size are straddled on the minimum threshold (effective population size = 500) indicated to maintain genetic diversity and evolutionary potential across several generations in natural populations.

Novel Tools for Conservation Genomics: Comparing Two High-Throughput Approaches for SNP Discovery in the Transcriptome of the European Hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Comparative analysis of AFLPs and SSRs efficiency in resolving population genetic structure of Mediterranean Solea vulgaris

The performance of different molecular markers in the assessment of population structure was tested using samples of Solea vulgaris collected in the Mediterranean within and outside the hypothetical dispersal ability of the species. A total of 172 individuals belonging to four population samples were analysed using 15 microsatellites [simple sequence repeats (SSRs)] and 153 amplified fragment length polymorphisms (AFLPs). Considering the global qualitative patterns, we found a correlation between SSRs and AFLPs in detecting genetic differentiation among samples. However, on a small geographical scale, AFLPs were able to discriminate individuals from neighbouring populations whereas SSRs were not, and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with SSRs. The high number of loci analysed with the AFLP technique could increase the probability to include outlier loci in the analysis; however, the neutrality test performed on our data set did not show evidence of selection acting on the S. vulgaris samples. Even if the choice of the molecular marker depends mainly on the biological question to be addressed, the higher power of discrimination and the comparative technical ease of obtaining data from AFLPs with respect to SSRs suggest the use of AFLPs for many population genetics studies.

Molecular Organization of 5S rDNAs in Rajidae (Chondrichthyes): Structural Features and Evolution of Piscine 5S rRNA Genes and Nontranscribed Intergenic Spacers

The genomic and gene organisation of 5S rDNA clusters have been extensively characterized in bony fish and eukaryotes, providing general issues for understanding the molecular evolution of this multigene DNA family. By contrast, the 5S rDNA features have been rarely investigated in cartilaginous fish (only three species). Here, we provide evidence for a dual 5S rDNA gene system in the Rajidae by sequence analysis of the coding region (5S) and adjacent nontranscribed spacer (NTS) in five Mediterranean species of rays (Rajidae), and in a large number of piscine taxa including lampreys and bony fish. As documented in several bony fish, two functional 5S rDNA types were found here also in the rajid genome: a short one (I) and a long one (II), distinguished by distinct 5S and NTS sequences. That the ancestral piscine genome had these two 5S rDNA loci might be argued from the occurrence of homologous dual gene systems that exist in several fish taxa and from 5S phylogenetic relationships. An extensive analysis of NTS-II sequences of Rajidae and Dasyatidae revealed the occurrence of large simple sequence repeat (SSR) regions that are formed by microsatellite arrays. The localization and organization of SSR within the NTS-II are conserved in Rajiformes since the Upper Cretaceous. The direct correlation between the SSRs extension and the NTS length indicated that they might play a role in the maintenance of the larger 5S rDNA clusters in rays. The phylogenetic analysis indicated that NTS-II is a valuable systematic tool limited to distantly related taxa of Rajiformes.

Mitochondrial DNA Sequence Variation Suggests the Lack of Genetic Heterogeneity in the Adriatic and Ionian Stocks of Sardina pilchardus

Abstract: A genetic stock structure analysis of 11 sardine samples from the Adriatic Sea and Ionian neighboring area was carried out through sequence variation analysis of a 307-bp cytochrome b gene fragment in order to identify self-recruiting units in the Adriatic Sardina pilchardus stock. The overall lack of genetic subdivision among samples detected by analysis of molecular variance, pairwise Φst values, and the exact test of population differentiation indicates this sardine stock is part of a larger self-recruiting population whose boundaries are larger than the investigated area. This conclusion is in agreement with preliminary allozymic and mitochondrial DNA restriction fragment length polymorphism data, but contradicts the previous identification of 2 subpopulations of sardines in the Adriatic Sea argued on morphologic differences, which could be rather attributed to different hydrographic or ecologic conditions occurring in different areas of the Adriatic Sea. The reduced gene flow observed between Adriatic-Ionian and Spanish sardine geographic samples (P < 0.001) suggests that reproductively isolated populations of sardines may occur in the Mediterranean Sea.

Rapid miocene-pliocene dispersal and evolution of Mediterranean rajid fauna as inferred by mitochondrial gene variation.

Rajidae (colloquially known as skates and rays) experienced multiple and parallel adaptive radiations allowing high species diversity and great differences of species composition between regional faunas. Nevertheless, they show considerable conservation of bio‐ecological, morphological and reproductive traits. The evolutionary history and dispersal of North‐east Atlantic and Mediterranean rajid fauna were investigated throughout the sequence analysis of the control region and 16S rDNA mitochondrial genes. Molecular estimates of divergence times indicated recent origin and rapid dispersal of the present species. Compared with the ancient origin of the family (Late Cretaceous), the present species diversity arose in a relatively narrow time‐window (12 Myr) from Middle Miocene to Early Pleistocene, likely by speciation processes related to dramatic geological and climatic events in the Mediterranean. Nucleotide substitution rates and phylogenetic relationships indicated Mediterranean endemic skates derived from sister species with wider distribution during Late Pliocene–Pleistocene. Skate phylogeny and systematics obtained using mitochondrial gene variation were largely consistent with those based on morpho‐anatomical data.

Microsatellite DNA variation reveals high gene flow and panmictic populations in the Adriatic shared stocks of the European squid and cuttlefish (Cephalopoda)

In the semienclosed Adriatic Sea, the shared stocks of the cephalopods Loligo vulgaris and Sepia officinalis represent important marine fisheries resources exploited by all coastal countries. The improving of knowledge on the demographic features of these shared stocks is internationally relevant for adopting responsible management and conservation of these marine resources. Analyses of microsatellite variation in geographical samples collected from all parts of the Adriatic Sea were performed using arrays of species-specific di-nucleotide and tri-nucleotide loci. In L. vulgaris the level of genetic variability was consistent with that observed in other loliginid species, whereas the S. officinalis stock showed a microsatellite variation markedly lower than that estimated for the Atlantic and Mediterranean populations collected around the Iberian peninsula. The weak spatial genetic differentiation, the discordant results of the genetic divergence estimators and the lack of any geographical cline in the spatial genetic differences suggest the occurrence of single genetically homogeneous populations within the Adriatic stocks of both species, recommending a coordinated management of the squid and cuttlefish by the Adriatic fishing countries. On the contrary, significant differences detected in temporal replicates of S. officinalis might suggest that allelic frequency can change relating to reproductive behaviour.