A comparative whole genome sequencing analysis identified a candidate locus for lack of operculum in cultivated gilthead seabream (Sparus aurata).

The gilthead seabream (Sparus aurata, Sparidae family) is commonly used for aquaculture. Despite its great economic value, several problems in its cultivation remain. One of the major concerns is the high frequency of morphological abnormalities occurring during the early developmental stages. Partial and/or total lack of operculum is the most frequent anomaly affecting the fish cranial region. The existence of genetic factors that can at least partially determine this defect has been hypothesized. In this work, two DNA pools of highly related fry, one composed of normal-looking (control) fish and the other lacking an operculum (case), were constructed and whole-genome resequencing data produced from the two were compared. The analysis revealed a 1 Mb region on chromosome 2 with higher heterozygosity in the lack of operculum DNA pool than in the control DNA pool, consistent with the enrichment, in the first DNA pool, of one or more haplotypes causing or predisposing to the defect together with other normal haplotypes. A window-based FST analysis between the two DNA pools indicated that the same region had the highest divergence score. This region contained 2921 SNVs, 10 of which, with predicted high impacts (three splice donor and seven stop-gained variants), were detected in novel genes that are homologous to calcium-sensing receptor-like genes, probably involved in bone development. Other studies are needed to clarify the genetic mechanisms involved in predisposing fry to this deformity and then to identify associated markers that could be used in breeding programs to reduce the frequency of this defect in the broodstock.

Identification of a major locus determining a pigmentation defect in cultivated gilthead seabream (Sparus aurata).

The gilthead seabream (Sparus aurata) is an important cultivated species in the Mediterranean area. A major problem for the gilthead seabream aquaculture sector derives from the high frequency of phenotypic abnormalities, including discolorations. In this study, we applied a whole-genome resequencing approach to identify a genomic region affecting a pigmentation defect that occurred in a cultivated S. aurata population. Two equimolar DNA pools were constructed using DNA extracted from 30 normally coloured and 21 non-pigmented fish collected among the offspring of the same broodstock nucleus. Whole-genome resequencing reads from the two DNA pools were aligned to the S. aurata draft genome and variant calling was performed. A whole-genome heterozygosity scan from single pool sequencing data highlighted a peak of reduced heterozygosity of approximately 5 Mbp on chromosome 6 in the non-pigmented pool that was not present in the normally coloured pool. The comparison of the non-pigmented with the normally coloured fish using a whole-genome FST analysis detected three main regions within the coordinates previously detected with the heterozygosity analysis. The results support the presence of a major locus affecting this discoloration defect in this fish population. The results of this study have practical applications, including the possibility of eliminating this defect from the breeding stock, with direct economic advantages derived from the reduction of discarded fry. Other studies are needed to identify the candidate gene and the causative mutation, which could add information to understand the complex biology of fish pigmentation.

Genetic variability is unrelated to growth and parasite infestation in natural populations of the European eel (Anguilla anguilla).

Positive correlations between individual genetic heterozygosity and fitness-related traits (HFCs) have been observed in organisms as diverse as plants, marine bivalves, fish or mammals. HFCs are not universal and the strength and stability of HFCs seem to be variable across species, populations and ages. We analysed the relationship between individual genetic variability and two different estimators of fitness in natural samples of European eel, growth rate (using back-calculated length-at-age 1, 2 and 3) and parasite infestation by the swimbladder nematode Anguillicola crassus. Despite using a large data set of 22 expressed sequence tags-derived microsatellite loci and a large sample size of 346 individuals, no heterozygote advantage was observed in terms of growth rate or parasite load. The lack of association was evidenced by (i) nonsignificant global HFCs, (ii) a Multivariate General Linear Model showing no effect of heterozygosity on fitness components, (iii) single-locus analysis showing a lower number of significant tests than the expected false discovery rate, (iv) sign tests showing only a significant departure from expectations at one component, and, (v) a random distribution of significant single-locus HFCs that was not consistent across fitness components or sampling sites. This contrasts with the positive association observed in farmed eels in a previous study using allozymes, which can be explained by the nature of the markers used, with the allozyme study including many loci involved in metabolic energy pathways, while the expressed sequence tags-linked microsatellites might be located in genes or in the proximity of genes uncoupled with metabolism/growth.

Detecting genome-wide gene transcription profiles associated with high pollution burden in the critically endangered European eel.

The European eel illustrates an example of a critically endangered fish species strongly affected by human stressors throughout its life cycle, in which pollution is considered to be one of the factors responsible for the decline of the stock. The objective of our study was to better understand the transcriptional response of European eels chronically exposed to pollutants in their natural environment. A total of 42 pre-migrating (silver) female eels from lowly, highly and extremely polluted environments in Belgium and, for comparative purposes, a lowly polluted habitat in Italy were measured for polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and brominated flame retardants (BFRs). Multipollutant level of bioaccumulation was linked to their genome-wide gene transcription using an eel-specific array of 14,913 annotated cDNAs. Shared responses to pollutant exposure were observed when comparing the highly polluted site in Belgium with the relatively clean sites in Belgium and Italy. First, an altered pattern of transcription of genes was associated with detoxification, with a novel European eel CYP3A gene and gluthatione S-transferase transcriptionally up-regulated. Second, an altered pattern of transcription of genes associated with the oxidative phosphorylation pathway, with the following genes involved in the generation of ATP being transcriptionally down-regulated in individuals from the highly polluted site: NADH dehydrogenase, succinate dehydrogenase, ubiquinol-cytochrome c reductase, cytochrome c oxidase and ATP synthase. Although we did not measure metabolism directly, seeing that the transcription level of many genes encoding enzymes involved in the mitochondrial respiratory chain and oxidative phosphorylation were down-regulated in the highly polluted site suggests that pollutants may have a significant effect on energy metabolism in these fish.