RESUMEN
Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.
RESUMEN
The gilthead sea bream (Sparus aurata L.) is a marine fish of great importance for fisheries and aquaculture. It has also a peculiar sex-determination system, being a protandrous hermaphrodite. Here we report the construction of a first-generation genetic linkage map for S. aurata, based on 204 microsatellite markers. Twenty-six linkage groups (LG) were found. The total map length was 1241.9 cM. The ratio between sex-specific map lengths was 1:1.2 (male:female). Comparison with a preliminary radiation hybrid (RH) map reveals a good concordance, as all markers located in a single LG are located in a single RH group, except for Ad-25 and CId-31. Comparison with the Tetraodon nigroviridis genome revealed a considerable number of evolutionary conserved regions (ECRs) between the two species. The mean size of ECRs was 182 bp (sequence identity 60-90%). Forty-one ECRs have a known chromosomal location in the pufferfish genome. Despite the limited number of anchoring points, significant syntenic relationships were found. The linkage map presented here provides a robust comparative framework for QTL analysis in S. aurata and is a step toward the identification of genetic loci involved both in the determination of economically important traits and in the individual timing of sex reversal.