Colour and pigmentation variation attracted fish biologists for a while, but high-throughput genomic studies investigating the molecular basis of body pigmentation remain still limited to few species and conservation biology issues ignored. Using 75,684 SNPs, we investigated the genomic basis of pigmentation pattern variation among individuals of the Atlantic and Mediterranean clades of the brown trout (Salmo trutta), a polytypic species in which Atlantic hatchery individuals are commonly used to supplement local wild populations. Using redundancy analyses and genome-wide association studies, a set of 384 “colour patterning loci” (CPL) significantly correlated with pigmentation traits such as the number of red and black spots on flanks, but also the presence of a large black stain on the opercular bone was identified. CPLs map onto 35 out of 40 brown trout linkage groups indicating a polygenic basis to pigmentation patterns. They are mostly located in coding regions (43.4%) of 223 candidate genes, and correspond to GO-terms known to be involved in pigmentation (e.g. calcium and ion-binding, cell adhesion). Annotated genes especially include 24 candidates with known pigmentation effects (e.g. SOX10, PEML, SLC45A2), but also the Gap-junction Δ2 (GJD2) gene that was previously showed be differentially expressed in trout skin. Patterns of admixture were found significantly distinct when using either the full SNP data set or the set of CPLs, indicating that pigmentation patterns accessible to practitioners are not a reliable proxy of genome-wide admixture. Consequences for management are discussed.

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single-and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single-and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium-and ionbinding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.