Performance and precision of double digestion RAD (ddRAD) genotyping in large multiplexed datasets of marine fish species

The development of Genotyping-By-Sequencing (GBS) technologies enables cost-effective analysis of large numbers of Single Nucleotide Polymorphisms (SNPs), especially in "non-model" species. Nevertheless, as such technologies enter a mature phase, biases and errors inherent to GBS are becoming evident. Here, we evaluated the performance of double digest Restriction enzyme Associated DNA (ddRAD) sequencing in SNP genotyping studies including high number of samples. Datasets of sequence data were generated from three marine teleost species (>5500 samples, >2.5 × 10¹² bases in total), using a standardized protocol. A common bioinformatics pipeline based on STACKS was established, with and without the use of a reference genome. We performed analyses throughout the production and analysis of ddRAD data in order to explore (i) the loss of information due to heterogeneous raw read number across samples; (ii) the discrepancy between expected and observed tag length and coverage; (iii) the performances of reference based vs. de novo approaches; (iv) the sources of potential genotyping errors of the library preparation/bioinformatics protocol, by comparing technical replicates. Our results showed use of a reference genome and a posteriori genotype correction improved genotyping precision. Individual read coverage was a key variable for reproducibility; variance in sequencing depth between loci in the same individual was also identified as an important factor and found to correlate to tag length. A comparison of downstream analysis carried out with ddRAD vs single SNP allele specific assay genotypes provided information about the levels of genotyping imprecision that can have a significant impact on allele frequency estimations and population assignment. The results and insights presented here will help to select and improve approaches to the analysis of large datasets based on RAD-like methodologies.

Sex change and effective population size: implications for population genetic studies in marine fish

Large variance in reproductive success is the primary factor that reduces effective population size (Ne) in natural populations. In sequentially hermaphroditic (sex-changing) fish, the sex ratio is typically skewed and biased towards the 'first' sex, while reproductive success increases considerably after sex change. Therefore, sex-changing fish populations are theoretically expected to have lower Ne than gonochorists (separate sexes), assuming all other parameters are essentially equal. In this study, we estimate Ne from genetic data collected from two ecologically similar species living along the eastern coast of South Africa: one gonochoristic, the 'santer' sea bream Cheimerius nufar, and one protogynous (female-first) sex changer, the 'slinger' sea bream Chrysoblephus puniceus. For both species, no evidence of genetic structuring, nor significant variation in genetic diversity, was found in the study area. Estimates of contemporary Ne were significantly lower in the protogynous species, but the same pattern was not apparent over historical timescales. Overall, our results show that sequential hermaphroditism may affect Ne differently over varying time frames, and that demographic signatures inferred from genetic markers with different inheritance modes also need to be interpreted cautiously, in relation to sex-changing life histories.

Exploring neutral and adaptive processes in expanding populations of gilthead sea bream, Sparus aurata L., in the North-East Atlantic

Recent studies in empirical population genetics have highlighted the importance of taking into account both neutral and adaptive genetic variation in characterizing microevolutionary dynamics. Here, we explore the genetic population structure and the footprints of selection in four populations of the warm-temperate coastal fish, the gilthead sea bream (Sparus aurata), whose recent northward expansion has been linked to climate change. Samples were collected at four Atlantic locations, including Spain, Portugal, France and the South of Ireland, and genetically assayed using a suite of species-specific markers, including 15 putatively neutral microsatellites and 23 expressed sequence tag-linked markers, as well as a portion of the mitochondrial DNA (mtDNA) control region. Two of the putatively neutral markers, Bld-10 and Ad-10, bore signatures of strong directional selection, particularly in the newly established Irish population, although the potential 'surfing effect' of rare alleles at the edge of the expansion front was also considered. Analyses after the removal of these loci suggest low but significant population structure likely affected by some degree of gene flow counteracting random genetic drift. No signal of historic divergence was detected at mtDNA. BLAST searches conducted with all 38 markers used failed to identify specific genomic regions associated to adaptive functions. However, the availability of genomic resources for this commercially valuable species is rapidly increasing, bringing us closer to the understanding of the interplay between selective and neutral evolutionary forces, shaping population divergence of an expanding species in a heterogeneous milieu.

Low cytochrome b variation in bream Abramis brama

Variability in cytochrome b (cytb) in European populations of bream Abramis brama was assessed. The cytb gene was found to be strongly conserved in A. brama relative to other cyprinid taxa. This limits the usefulness of this marker in examining geographical genetic structure in this species and raises interesting questions as to the recent evolutionary history of the species.

A highly permeable species boundary between two anadromous fishes

Meristic identification, mitochondrial DNA and a suite of microsatellite markers were employed to estimate the incidence of hybridization in wild populations of anadromous Allis shad Alosa alosa and twaite shad Alosa fallax in southern Irish riverine and estuarine waters. It was shown that 16% of the fishes examined were misclassified using meristic count of gill rakers. Next, a significant proportion of fishes that were robustly assigned to a species using nuclear markers were shown to possess the mtDNA of the other. The genomes of A. alosa and A. fallax in Ireland are extensively introgressed, which suggests a complex history of hybridization between these species, which can only partially be explained by recent man-made habitat changes.