Skip to main content

Comparative transcriptome analysis of trout skin pigment cells



Enormous variability in skin colour and patterning is a characteristic of teleost fish, including Salmonidae fishes, which present themselves as a suitable model for studying mechanisms of pigment patterning. In order to screen for candidate genes potentially involved in the specific skin pigment pattern in marble trout (labyrinthine skin pattern) and brown trout (spotted skin pattern), we conducted comparative transcriptome analysis between differently pigmented dermis sections of the adult skin of the two species.


Differentially expressed genes (DEGs) possibly associated with skin pigment pattern were identified. The expression profile of 27 DEGs was further tested with quantitative real-time PCR on a larger number of samples. Expression of a subset of ten of these genes was analysed in hybrid (marble x brown) trout individuals and compared with the complexity of their skin pigment pattern. A correlation between the phenotype and the expression profile assessed for hybrid individuals was detected for four (gja5, clcn2, cdkn1a and tjp1) of the ten candidate genes tested. The potential role of these genes in skin pigment pattern maintenance is discussed.


Our results indicate that the maintenance of different pigment patterns in trout is dependent upon specific communication—involving gap junctions, tight junctions and ion channels—between chromatophores present in differentially pigmented skin regions.


Skin pigmentation and pigment patterning play an important role in the survival of animals. Pigments can influence thermoregulation, social communication, mate choice and predator avoidance, among others, and provide protection against the harmful effect of solar radiation [1, 2]. Even closely related animals exhibit extremely divergent patterns indicative of great plasticity in skin patterning [3]. Enormous variability in skin colour and pattern is characteristic of teleost fish, in which there are six known pigment cell types: melanophores, iridophores, xanthophores, erythrophores, leucophores and cyanophores [4, 5]. These chromatophores synthesise and retain their specific pigments and inner structures intracellularly [6]. Skin pigmentation and pattern thus reflect the numbers and arrangements of the chromatophores themselves. Although only one pigment cell type is known in mammals, there is much conservation of pigmentation-related genes in teleosts and mammals [7]. Therefore, given also the great variability created by numerous pigment cell types, teleosts are important models for studying the cellular, physiological, genetic and environmental factors involved in pigment pattern formation and maintenance [2, 8, 9].

The genetic basis of pigment production and subsequent skin pigmentation is well understood in model species such as zebrafish [10,11,12], and some mechanisms regulating pigment patterning in zebrafish have been recently proposed (e.g. [13,14,15]). Many studies conducted on zebrafish show a requirement for interactions among three types of pigment cell—melanophore, xanthophore and iridophore—for stripe formation in adult fish [16,17,18,19] and the importance of cell–cell communication in the formation and maintenance of different pigment patterns. Through induction of single gene mutations, important roles in colour patterning have been identified in ion channels [20], tight junctions [21] and gap junctions [16, 22, 23]. The results of those studies show that change in bioelectrical communication between pigment cells can greatly affect skin pigment pattern [24]. Along with chemical diffusion bioelectric signalling mechanisms are fundamental properties for pattern self-generation, as envisaged in the Turing theoretical model for morphogenesis in living systems [24,25,26].

Much research on zebrafish has been undertaken on artificially produced mutants. However, much less is known about the genetics of skin pigment patterns in natural populations. Salmonidae fishes are characterised by a huge variety of naturally present pigment patterns and thus present themselves as a suitable model for studying the mechanism of pigment pattern formation and maintenance.

One factor that may cause great variety in salmonids could be the relatively recent whole genome duplication (WGD) event (salmonid-specific 4th WGD, or Ss4R [27],) that took place some 25–100 million years ago. Analysis of the rainbow trout genome [28] has shown that half (48%) of the protein-coding genes retain both ohnologues (paralogues formed by a WGD event). Several instances of gene duplication in pigment pathways have also been observed [7, 29, 30] and may have led to diversity of skin colour and patterns in this group of fishes.

In the present study, we focused on a very closely related [31], yet phenotypically divergent, sister taxon pair in the genus Salmo: marble trout (S. marmoratus), with labyrinthine skin pattern, and brown trout (S. trutta), with spotted skin pattern (Fig. 1). Changes in skin pigmentation and pattern can occur over surprisingly short evolutionary timescales, as in these Salmo species, while the reoccurrence of similar colour patterns across large phylogenetic distances is common. A labyrinthine pattern is found in pufferfish (Takifugu exascurus), some zebrafish mutants (Cx41.8 M7 [3]), partially in jaguar/obelix [20], and also in an isolated population of brown trout inhabiting the River Otra, Norway [32], as well as in some hybrids between brown trout and brook trout (S. trutta x Salvelinus fontinalis [33]).

Fig. 1

Marble and brown trout and their typical skin pigment pattern: a, Labyrinthine pattern on the skin of marble trout; b, The part of the dermis used for preparing cryosections; Scale bar: 50 μm; c, Spot pattern on the skin of brown trout. Circles in a and c highlight the differently pigmented parts dissected with biopsy punches

Morphological differences in chromatophore ultrastructure, and position or arrangement, or both, in the dermis of marble and brown trout have been revealed using transmission electron microscopy [34]. In zebrafish, position, interaction and presence or absence of different pigment cells provide the basis for different skin pigmentation and pattern. In these two trout species, however, the position or arrangement of the chromatophores in differently pigmented skin regions is much more complex, with a higher level of organisation of chromatophores observed in the skin of brown trout than that of marble trout. Meanwhile, a new pigment cell type, erythrophore type 2, and its ultrastructure in brown trout have been described for the first time in salmonids [34].

Recently, microarray analysis has been undertaken to compare gene expression profiles of whole skin samples of marble and brown trout, with dermis and epidermis, and all differently pigmented regions, considered as a single sample [35]. Four candidate genes for labyrinthine skin pattern have been described—hdac1, vps18, dct and scg2a—and it has been proposed that the formation of the labyrinthine pattern at least partially depends upon the Wnt signalling pathway and is based upon a reaction–diffusion mechanism. Only genes involved in the biological process ‘pigmentation during development’ were considered as informative [35]. In the present study, to further explore the pathways and genes expressed in the skin and potentially causative to skin pigment pattern in Salmo, we employed high-throughput mRNA sequencing technology on differentially pigmented dermis sections of marble trout and brown trout. These first RNA-seq data on differently pigmented dermis regions for these two species enabled (i) an overview of the transcriptome in dark, white and red coloured skin regions, and (ii) identification of differentially expressed genes (DEGs) possibly associated with skin pigment pattern. The expression profile of 27 DEGs was tested further with quantitative real-time PCR (qPCR) on a larger number of samples, while expression of a subset of ten of these genes was analysed in hybrid trout individuals (marble x brown trout) and compared with the complexity of their skin pigment pattern.


Illumina sequencing, data assembly and annotation

Sequencing of RNA isolated from trout dermis resulted in a total of 812,144,586 pair-end reads of 100 bp each (414,612,488 in marble trout and 397,532,098 in brown trout) that were mapped to the gene model of Atlantic salmon (S. salar) genome (INSDC: AGKD00000000.4) [36]. All transcripts of the gene model were annotated with Swissprot to obtain the Gene symbols for Gene Ontology analysis. The assembled transcriptomes were submitted to the SRA database under accession number SRP157513. Results of transcriptome sequencing and mapping to the salmon genome are reported in Additional file 1.

Genes highly expressed in trout skin

The top 50 genes most highly expressed in trout skin included genes of the keratin family and collagen proteins, with keratin, type II cytoskeletal cochlear-like the most highly expressed. The other genes highly expressed in marble and brown trout skin are reported in Additional files 2 and 3.

Differential gene expression

Using NOISeq numerous DEGs between marble and brown trout pigment cell genes were identified. A total of 52,012 predicted genes were identified as expressed differentially between marble and brown trout skin. Of these, 3715 genes were up-regulated (≥ 2 fold) and 4327 down-regulated (≤ 2 fold) in the skin of marble trout compared to that of brown trout. Among the up- and down-regulated genes several are known to encode membrane proteins involved in, for example, signalling, cell junctions, or ion permeability. The 100 most up- and down-regulated genes that differed between marble and brown trout samples are reported in Additional files 4 and 5. Several key genes involved in skin pigmentation (e.g. dct, tyr, tyrp1, pmel, ednrb) showed differences in expression profiles between the two species, as did some genes involved in cell–cell junctions (e.g. gja1, gja5, gjb6, gjd2, mpp3, tjp1) and the establishment of membrane and action potentials (e.g. clcn2, ctr2, kcnk3). Among these groups of differentially expressed genes, 27 were selected for additional testing by qPCR.

K-mean clustering was used to group genes with a similar pattern of expression among samples. Many genes were overexpressed specifically in the red spot of brown trout compared to other skin regions of both species, indicating that there are specific pathways involved in either carotenoid metabolism or communication between these specific pigment cells and their surrounding cells, or both, resulting in the formation of a red spot.

Functional classification of DEGs

DEGs were annotated through GO classification analysis and grouped into three categories (biological process, cellular component and molecular function) based on their putative functions. Salmon genes were annotated with BLAST [37, 38] against Swissprot and Swissprot Gene Ontology annotations at

For the GO analysis, 18,087, 61,751 and 87,737 DEGs between marble and brown trout were grouped by molecular function, cellular component and biological process categories, respectively. Altogether, DEGs were mapped to 61 categories, as presented in Fig. 2. The GO terms containing the largest number of DEGs were cell, cell part, cellular process, single-organism process, binding and membrane.

Fig. 2

GO functional classification of transcripts according to their assigned molecular function, cellular component and biological process (up-regulated, MT > BT; down-regulated, BT > MT)

KEGG pathway analysis

KEGG pathway analysis was carried out to categorise and annotate the DEGs, which were mapped and visualised in a way similar to the KEGG Mapper tool. The pathways with the largest number of genes expressed differentially between marble and brown trout were Neuroactive ligand-receptor interaction, NOD-like receptor signalling pathway, PI3K-Akt signalling pathway, Pathways in cancer, Focal adhesion, and MAPK signalling pathway (Fig. 3a). Other interesting differentially expressed pathways included: Cell adhesion molecules (CAMs), Calcium signalling pathway, Tight junction, Adherens junction, Gap junction, Cell cycle, p53 signalling pathway, Wnt signalling pathway, Melanogenesis and Tyrosine metabolism, all of which could affect and explain the difference observed between the two species in their skin pigment cell organisation. Expression profiles of DEGs from Melanogenesis and Gap junction are presented in Fig. 3b and c. The most represented pathways included, interestingly, the Thyroid hormone signalling pathway, Insulin signalling pathway and Insulin secretion. Many genes involved in these pathways were detected also during analysis of the first 100 up- and down-regulated transcripts with the highest fold change.

Fig. 3

KEGG pathways: a, KEGG pathway classification of genes (green, up-regulated, MT > BT; red, down-regulated, BT > MT); b, Heat map visualisation of the expression of genes involved in melanogenesis; c, Heat map visualisation of the expression of genes involved in gap junctions. The columns display the samples and the rows the genes (red, down-regulated; green, up-regulated)

When comparing only the dark pigmented regions of the two species, the pathways with the largest number of DEGs were those describing neurodegenerative disorders, such as Huntington’s, Alzheimer’s and Parkinson’s diseases. When comparing only the light regions, the highest represented pathways were those associated with cancer, PI3K-Akt signalling, Focal adhesion, Huntington’s disease and Alzheimer’s disease.

The highest absolute number of DEGs was observed when comparing red skin with other skin regions in brown or marble trout. Focal adhesion, ECM receptor interaction, Cell adhesion molecules and MAPK signalling pathway were among the pathways with the largest number of genes represented and with a large differential expression profile.

Quantitative real-time PCR

In order to validate the transcriptome sequencing results, we selected 27 DEGs for qPCR to determine their relative expression in marble and brown trout skin (Table 1). These genes could be classified into two groups based on their known functions in model organisms, including zebrafish: (1) genes involved in pigmentation pathways—tyr, tyrp1, dct (tyrp2), mitf, pmel, mc1r, ednrb, sox10, aox1, cdkn1a, wnt10a and wnt9b—, and (2) genes with a molecular function in cell binding known to be involved in ion exchange, chromatophore interactions and pigment pattern: gja1, gja4, gja5, gja9, gjb1, gjb3, gjb4, gjd2, tjp1, tspan3, slc7a2, slc01c1, kcnc2, clcn2 and pcdhac2. qPCR (Fig. 4) showed that six (dct, wnt10a, gja9, gjd2, tjp1, pcdhac2) of the 27 selected genes had a significantly higher level of expression in brown trout skin than in marble trout skin, four others (mitf, mc1r, sox10, ednrb) were expressed significantly more, and one (tyr) significantly less, in the red spot of brown trout, while five genes (cdkn1a, gjb1, slc7a2, kcnc2, clcn2) had significantly higher expression levels in the skin of marble trout than in brown trout. Gene gja5 exhibited an interesting pattern of expression, significantly higher in (black and red) spots of brown trout than in other skin parts of either species. The qPCR results were mostly consistent with the data from transcriptome sequencing, with 23 out of 27 genes having the same expression profile.

Table 1 Primers used for the qPCR analysis
Fig. 4

Gene expression patterns obtained with qPCR. Log2 fold changes are expressed as the ratio of gene expression after normalisation to reference genes. Bars represent the means ± SD. Grey, brown trout dark spot; light grey, brown trout light region; red, brown trout red spot; blue, marble trout dark region; light blue, marble trout light region. Expression of each gene was compared with expression in the black spot region of brown trout and statistically evaluated by using unpaired Student’s t-test. A p-value of < 0.05 was considered statistically significant (*p < 0.05; **p < 0.01; ***p < 0.001)

Correlation between expression of candidate genes and phenotype of hybrids

From analysis of the qPCR products, ten candidate genes (dct, cdkn1a, gja5, tjp1, wnt10a, gjd2, slc7a2, kcnc2, clcn2 and pcdhac2) with the most differential expression between the two species were selected and their expression further tested in differently pigmented skin regions of hybrid samples.

The skin pigment pattern variation of hybrid animals was evaluated and compared to those of brown trout and marble trout. The complexity of colour pattern in each individual fish, with a lower score representing greater complexity, was plotted against the overall colour tone of the pattern, defined as the proportion of white area [33] (Fig. 5). The complexity of pattern variation of hybrids was then compared with the expression of candidate genes. The expression of clcn2, cdkn1a and tjp1 was partially consistent with pattern complexity, while clcn2 and cdkn1a demonstrated a higher level of expression in hybrids closer to the marble phenotype than to the brown trout phenotype, and tjp1 had a higher expression in hybrids closer to the brown trout phenotype (Fig. 6). The overall expression of gja5 was partly correlated with pattern complexity, though the red spot samples showed a strict overexpression in hybrids with a phenotype closer to marble trout (Fig. 6). Figures 5 and 6 show the average expression value for each hybrid, and expression values for every differently pigmented region of their skin is shown in Additional file 6. For other genes no correlation between expression profile and pattern complexity was observed.

Fig. 5

Skin pigment pattern variation: a, Skin pigment pattern variation in genetically pure marble and brown trout, and in hybrid animals; b, quantification of pigment patterns in hybrids (pink dots) displaying different patterns, ranging from brown trout phenotype (grey dots) to marble trout phenotype (blue dots)

Fig. 6

Gene expression patterns of four candidate genes in hybrids: a, expression of cdkn1a gene averaged over skin regions in hybrids; b, expression of tjp1 gene averaged over skin regions in hybrids; c, expression of clcn2 gene averaged over skin regions in hybrids; d, expression of gja5 gene in the red spots of hybrids


Various mechanisms are thought to be responsible for the great variety in skin pigment patterns observed in animals. Recently, the role of homotypic and heterotypic intercellular interactions has been demonstrated in pigment progenitor cell migration, and in the formation and maintenance of a specific pigment pattern in zebrafish [10, 17,18,19]. In order to screen for genes potentially involved in maintenance of a specific pigment pattern in marble and brown trout, we conducted comparative transcriptome analysis of differently pigmented regions of the adult skin of the two species. The results indicate that the maintenance of different skin pigment patterns in trout is dependent upon specific communication—involving gap junctions, tight junctions and ion channels—between chromatophores present in differentially pigmented skin regions. Some candidate genes are involved in the distribution and differentiation of pigment progenitor cells [3, 39].

DEGs and genes potentially correlated with skin pigment pattern

KEGG pathway analysis showed that a substantial number of DEGs between the two species are associated with pigmentation-related pathways—e.g. MAPK signalling pathway, Wnt signalling pathway, Melanogenesis, Tyrosine metabolism—and cell–cell communication, such as Focal adhesion, Cell adhesion molecules (CAMs), Calcium signalling pathway, Tight junction, Adherens junction and Gap junctions.

Differential expressions were also observed when differently coloured regions of skin from the same individual fish were compared, most probably associated with specific expression profiles of individual pigment cell types predominating in a particularly coloured skin region [34]. When comparing only dark regions or only light regions between the two species, several DEGs were found. These genes are involved in Ion binding and exchange, Ion flow, Cell–cell junctions, Cellular signal transductions (Ras subfamily of proteins) and Actin regulation (Rho GTPases). Genes involved in Melanogenesis (pmel, tyrp1, dct) were overexpressed in the black spot of brown trout when compared to other pigmented regions of brown trout skin.

The results from the present study indicate that DEGs with a possible role in trout skin pigment pattern maintenance can be classified into two major groups: (1) genes involved in pathways that have a known role in pigmentation (e.g. mc1r, tyrp1, dct, wnt6), and (2) genes with a molecular function of binding, including those coding for proteins involved in cell–cell communication and ion channels (e.g. gjd2, clcn2, kcnk3, kcj12), known to act as factors that can affect the organization of pigment cells in the dermis [3, 23]. We selected 27 genes from both groups for qPCR confirmation of differential expression. RNA isolated from differently pigmented regions of ten brown and ten marble trout samples were used for individual qPCR validation of the sequencing results. The expression profile of most of the DEGs were in agreement with the results of the comparative transcriptome analysis, indicating that there was no consistent bias in the expression profile.

In addition, expression of ten candidate genes, with consistent differences in expression between marble and brown trout skin, was analysed in hybrid individuals and correlation with phenotype complexity was assessed. The phenotype complexity was comparable with various levels of brown trout introgression (Additional file 7), estimated from the genotyping of 36 SNPs.

Of four differentially expressed informative genes—hdac1, vps18, dct and scg2a—previously proposed to be involved in skin pigment pattern maintenance in marble and brown trout [35], only dct demonstrated a large difference in the present study, and was expressed more than twice as much in the marble trout skin than in the brown trout skin. As explained below, different expression of this gene does not correlate with the phenotype of hybrid individuals, nor do the remaining genes that are involved in pathways that have a known role in pigmentation. These other genes are important for pigment synthesis; their differential expression thus results in different pigmentation intensity and, most probably, not in a different pigment pattern.

Candidate genes determined from correlation between gene expression and phenotype

A positive correlation between phenotype and the expression profiles assessed in hybrid individuals with a variety of different skin pigment patterns was observed for four of the tested candidate genes: gja5, clcn2, cdkn1a and tjp1. Each of these genes is included in the above-mentioned second group of candidate genes, with a molecular function of binding and cell–cell communication. (However, their expression was not entirely consistent with the assessed complexity—or simplicity—of the pattern, perhaps from an accumulated bias from the steps of pattern evaluation: e.g. quality of photographs taken in the field; the dark parts of the hybrid skin not being separated into spotted and non-spotted regions; parr marks present in some hybrids affecting evaluation of the pattern.)

Gene cdkn1a (also referred to in the literature as p21) encodes a cyclin dependent kinase inhibitor 1a that was up-regulated in skin with labyrinthine pattern (marble trout). This gene functions as a regulator of cell-cycle progression in the phases G1 and S. Furthermore, it has been shown that cdkn1a is involved in patterning of eye-spots in the pigmented epithelium of the retina [40]. Another interesting role of this gene involves the regulation of cell migration and actin stress fibre formation through regulation of the Rho pathway [41], which could have a role also in skin pigment patterning. Finally, cdkn1a is known to be important in the development of melanocytes, promoting their differentiation and melanogenesis [42].

Šestáková et al. [43] reported a higher level of P21 protein in lightly pigmented, faster proliferating melanocytes than in dark-pigmented, slowly proliferating melanocytes. In a previous study, we observed a different ultrastructure of the melanophores present in the dark spot than in the light region of brown trout skin [34]. In addition, we noticed two types of melanophores—stagnate and active—present in the primary culture of pigment cells isolated from brown trout skin, with a slightly different colour and different speed of movement (data not shown, available from authors upon request). The presence of only one type of melanophore was noticed in the skin of marble trout, reflecting mostly the characteristics of the active type of melanophore. We therefore suggest that there may be up to two types of melanophore present in trout skin, with a higher share of slow, almost stagnate, melanophores in brown trout compared to a higher share of fast active melanophores in marble trout. A large proportion of the stagnate type melanophore in dark spots of brown trout skin indicates it could potentially be involved in spot maintenance. Also, the movement of the two types of melanophore may be differentially regulated through cdkn1a and the Rho-pathway and thus could affect the final skin pigment pattern. To substantiate our hypothesis KEGG analysis of the light region versus the black spot was conducted. Pathways with the largest numbers of DEGs included PI3K-Akt signalling pathway, Focal adhesion, ECM-receptor interaction, Cell adhesion molecules, Regulation of actin cytoskeleton and Axon guidance. These results show that different cell movement and attachments could indeed affect the perpetuation of spots.

A second candidate gene, clcn2, encodes a plasma membrane voltage-gated chloride channel Clcn2 (also known as ClC-2). As with other chloride channels, this is a transmembrane protein involved in the transport of chloride ions and, through that, the maintenance of chloride ion homeostasis and the transmission of electrical signals between cells. The functions of chloride channels include membrane potential stabilization, signal transduction, regulation of cell volume and of a tight junction, and transepithelial transport [44]. To the best of our knowledge, the clcn2 gene has not been described previously to be involved in pigment patterning. Meanwhile, another chloride channel—CIC-7—has been shown to have an important role in pigmentation in mice [45], while in zebrafish a potassium channel—Kir7.1—has been connected to specific pigment cell interactions and depolarizations, affecting pigment patterning [20].

The other two candidate genes, tjp1 and gja5, are involved in cell–cell junctions. The gene tjp1 (tight junction protein 1, also known as zo-1) encodes a protein involved in tight junctions and was up-regulated in spotted skin (red and black spots in brown trout). It has already been demonstrated that the expression of this gene can affect pigment patterning. In a previous study, a zebrafish with a mutation that disrupts the functioning of the tjp1a protein was examined [21]. The mutant had a spotted pattern rather than the usual striped pattern. It was shown that tjp1a was expressed in dense iridophores but its expression was down-regulated in loose-form iridophores and in melanophores and xanthophores. In the mutant, the dense iridophores interrupted the dark stripes, which reorganized to form spots. Iridophores in trout skin were not recognized as loose or dense forms, but rather were classified based upon the size of their inclusions, into large (L) and small (S) types [34]. Nevertheless, we hypothesize that tjp1 might play a distinct role in the functioning of these iridophores. Namely, type L iridophores were present exclusively in the light area of brown trout skin [34]. With the overexpression of tjp1, these cells could prevent melanophores from forming or from maintaining the labyrinthine pattern, or both, but rather organizing into spots. Another, very likely, role of the tjp1 gene in skin pigment patterning is connected to the observation that Tjp1 protein is a scaffolding protein that links tight junction proteins (such as claudins and occludins) to the actin cytoskeleton and could, as a consequence, affect the rearrangement of the cytoskeleton and the movement of the cells [46]. Together with cdkn1a, which has a similar role, the expression profile of these two genes could explain the diverse movements of two types of melanophores and their communication with other pigment cells resulting in different pigment patterns [47, 48].

The gene gja5 (gap junction protein alpha 5, also known as cx40) was overexpressed in dark and red spots in brown trout. This gene codes for a protein that is a member of the connexin family, which forms another type of cell junction: gap junctions. These are connections between cells, have a role in cell–cell communication and are involved in the maintenance of ionic and metabolic homeostasis and transmission of various ions and electrical signals [49]. Dark and red spots in brown trout are characterized by both melanophores (stagnate, see above) and erythrophores [34]. Thus, it is hypothesized that gja5 might play an important role in communication (heterotypic interaction) between these two cell types or between erythrophores (homotypic interaction), resulting in a spotted pattern. The role of gap junction proteins in pigment pattern formation and maintenance has been confirmed several times in zebrafish, mainly through studying a range of mutants with various pigment patterns as a consequence of the modified cx41.8 (also known as gja5b) gene (missense point mutations, deletions) [3, 22, 23, 50], a homolog of gja5. Nevertheless, to the best of our knowledge, this is the first study demonstrating a connection between gja5 expression level and skin pigment pattern. Differential expression of another connexin—gja1—between light and dark pigmented regions has been shown previously in sticklebacks [51], where a higher level of expression was detected in tissue from light bars. The high expression level of gja5 in red spots of hybrids showed a direct correlation with phenotype, with an overexpression in red spots present in skin with a general marble trout phenotype, or labyrinthine pattern. We additionally examined the correlation between gja5 expression in red spots and overall colour tone that is directly connected to the labyrinthine pattern. The hybrids formed two groups, where higher gja5 expression was connected to a darker overall colour tone (Additional file 8). A possible explanation for this, at first sight paradoxical, result involves cell interactions and the specific cell environment [52,53,54], with possible indirect mechanisms involved in pigment pattern maintenance and chromatophore interactions [13,14,15,16, 18]. Erythrophores forming red spots in brown trout and some hybrid individuals are otherwise not present in marble trout skin (labyrinthine pattern). We hypothesise that more intense specific homotypic (between erythrophores) or heterotypic (erythrophore–melanophore), or both, interaction and molecular or signal transductions are required in hybrids in order to maintain the red spot pattern on a labyrinthine pattern background. Since the gap junction works as a channel, connecting the cytoplasm of two cells and permitting the passage of various ions and molecules, and is affected by depolarisation, the roles of all four candidate genes for pigment pattern maintenance in Salmo could be tightly entwined. It has been shown that melanophores present in zebrafish jaguar mutants, with a mutation in the potassium channel Kir7.1, were constantly depolarised, and not only when in contact with xanthophores, as observed in wild-type melanophores [55]. The consequence of this constant depolarisation is a variant pattern in jaguar mutant skin, where usual zebrafish stripes are interrupted. On the other hand, depolarisation could also affect the actin cytoskeleton rearrangement [56]. The involvement of all four candidate genes suggests the importance of bioelectrical signalling, depolarisation, cytoskeleton rearrangement and cellular movement in skin pigment pattern maintenance in trout. To confirm this hypothesis further studies need to be conducted.


In summary, we sequenced transcripts from the skin of marble trout and brown trout for the first time and identified more than 8000 DEGs (with a fold change ≥2) in the skin of the two species. These DEG data provide a reference for screening genes for their effect on pigment patterning and other processes taking place in trout skin. The results significantly enhance our understanding of the composition of the trout skin transcriptome and the potential differences in gene expression associated with skin pigment pattern maintenance, and provide a foundation for future studies. We proved the validity of the mRNA-seq results using qPCR, where 27 candidate genes were tested on a large number of samples. Our results suggest an important role of cell–cell interactions, junctions and cytoskeleton rearrangements in trout skin pigment pattern maintenance as a consequence of differential gene expression.


Experimental animals

Adult marble trout individuals (of age > 2 years) were collected from Tolminka fish farm in Tolmin, Slovenia, and adult brown trout (of age > 2 years) from Bled fish farm, Bled, Slovenia. With an age of at least two years these trout had reached sexual maturity and already developed their adult skin pigment pattern. Their body length was approximately 30 cm. All fish were first-generation offspring reared in their respective fish farm from wild-caught parents originating from Zadlaščica stream (marble trout) and Malešnica stream (brown trout), Slovenia. Individuals had been fed with Biomar INICIO Plus (Denmark) as a starter feed and Biomar EFICO Enviro 920 when adult. Prior to skin sample collection fish were sedated in anaesthetic Tricaine-S (MS-222, Western chemical, Ferndale, USA) and killed by a blow to the neck.

In order to verify any correlation between the expression of ten candidate genes and the skin pigment pattern, hybrids between the two species with a variety of pigment patterns were sampled from Volarja stream (Soča river basin, Slovenia) where stocking with brown trout has been intensive for several decades and where ‘pure’ marble trout could not be found due to hybridisation. Hybrids with various levels of brown trout introgression are a result of hybridisation between marble and introduced brown trout, and backcrosses. Level of introgression was assessed by genotyping a panel of 36 diagnostic nuclear SNPs as reported previously [57]. The introgression levels for hybrid individuals are reported in Additional file 7. Skin from thirteen hybrid trout was selectively dissected with biopsy punches from differentially pigmented parts (light and dark region, red spots where present). In total, 37 different samples were collected.

All methods described were carried out in accordance with relevant guidelines and regulations. Fish skin sampling was approved by the Ministry of Agriculture and Environment, Slovenia, under decision letter number U34401–60/2013/4.

Sample preparation

Small pieces of skin from the lateral part of the trunk of the body (Fig. 1) were obtained with a 2 mm biopsy punch (Kai Group, Japan) to selectively dissect differently pigmented regions of skin in marble trout (dark (MTD) and light (MTL) region) and brown trout (light region (BTL), black (BTB) and red (BTR) spots). Skin pieces were immediately immersed in liquid nitrogen and subsequently used for preparation of cryosections. Only sections of the dermis where pigment cells are located (see [34] for more details about pigment cell location) were used for total RNA extraction (Fig. 1).

Total RNA extraction and transcriptome sequencing

Total RNA was isolated from cryosections of skin using a PicoPure RNA Isolation Kit (Thermo Fisher Scientific, Waltham, MA USA) according to manufacturer’s instructions. All samples were treated with DNA-free DNA Removal Kit (Thermo Fisher Scientific) to remove traces of genomic DNA. Concentration of total RNA and RNA integrity value (RIN) were checked using an RNA 6000 Pico LabChip of Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). RNA extracts with the highest RNA integrity of identically pigmented skin samples from three individuals per species were pooled. All RIN values were > 8.5. In total, five libraries from 1 μg of total RNA isolated and pooled from six individual fish (three per species), labelled as MTD, MTL, BTB, BTL and BTR, were prepared using random hexamer primers. cDNA libraries were sequenced by Beijing Genomics Institute (Hong Kong, China) using 2x100bp pair-end high-throughput mRNA sequencing (RNA-seq) on two lanes (1st lane with two marble trout samples, 2nd lane with three brown trout samples) of Illumina Hiseq 2000 (Illumina, San Diego, CA).

Bioinformatical analysis of RNA-seq

Paired end sequences from the mRNA-seq were mapped with Bowtie2 [58] against the transcripts (gene model) of Atlantic salmon genome (ICSASG_v2) applying standard parameters, and a raw count table was extracted using SAMtools and the command ‘Samtools idxstats’ (SAMtools version 0.1.19-96b5f2294a). Reads were mapped directly to the gene model of Atlantic salmon, as this proved to be less error prone than mapping to the genome due to sequence differences between Atlantic salmon and marble and brown trout and difficult prediction of exon/intron boundaries. The read counts were normalized to total reads in the sample with the smallest number of reads. Genes for which the read count was less than ten in the highest expressed sample were discarded and differentially expressed genes (DEG) were computed with NOISeq [59].

Differential gene expression

Differential expression analysis was performed using the R/Bioconductor package NOISeq (1-PNOI < 0.05, equivalent to false discovery rate [FDR] adjusted P-value [60, 61]). For NOISeq, the probability of a gene being differentially expressed that was provided by the method was used for the threshold (0.8). Since there were no true biological replicates available, technical replicates were simulated by NOISeq-sim, where single samples were used in NOISeq and additional samples from the same species run as replicates in making comparisons between the species. As NOISes-sim only simulates technical replicates, we accounted for some of the biological variability by using RNA from three individual fish per sample for transcriptome sequencing. The differential expression of selected genes was validated by running qPCR on multiple samples. qPCR also provided information on the biological variation of the genes of interest. We ran NOISeq for comparison on pooled libraries, using the library-normalised raw read count for each gene. K-mean clustering was used to group genes with similar patterns of expression. Gene Ontology [62, 63] was used to pick out genes that were annotated with ECM. Gene model transcripts of Atlantic salmon were annotated by BLAST+ against Swissprot to obtain the gene symbol for each reference gene. Gene Ontology annotation terms for Swissprot were downloaded from The output ECM gene list could be cross-checked for DEGs from NOISeq. KEGG pathway analysis [64] was performed by mapping the KEGG annotated DEGs from NOISeq to KEGG pathways as described in the KEGG Mapper tool. Both raw expressions of genes and DEGs as fold change were plotted in pathways, and the ratio of number of up-regulated genes to number of down-regulated genes, or vice versa, were calculated as a means of ranking up- or down-regulated pathways.

Quantitative real-time PCR

Samples for RNA isolation were prepared using biopsy punches to selectively dissect differently pigmented parts of the skin. Skin pieces were then homogenised with Percellys24 homogenizer (Bertin corp.) with the use of zirconium oxide beads (Bertin corp). Total RNA was isolated using RNeasy Plus Universal Mini Kit according to manufacturer’s instructions and subsequently treated with RNase-free DNase I (Thermo Fischer Scientific). Only RNA with good purity, checked with a NanoVue spectrophotometer (GE Healthcare, Little Chalfont, Buckinghamshire, UK), with values A260/280 ~ 2.0 was used for cDNA synthesis, performed with High-Capacity cDNA Reverse Transcription Kit (Thermo Fischer Scientific).

qPCR was performed using SYBR Green PCR Master Mix (Thermo Fischer Scientific, MA USA) and conducted on a Viia7 Real-Time PCR System (Applied Biosystem, Thermo Fisher Scientific). All primers for selected candidate genes (Table 1) were designed based on transcriptome sequences (SRP157513) obtained using Primer-BLAST [65] with an amplicon size of approximately 100–200 bp. Ten potential reference genes—rps20, hprt1, RNA polII, ppia1, pgk, sdha, b2m, ef1a, β-actin, 18 s—were analysed to identify two that could be used as reference genes: rps20 and pgk1 were chosen as having a stable expression, i.e. at the same level in all tested samples. The expression of candidate genes was analysed with ten biological replicates of marble (dark, MTD; light region, MTL) and ten of brown trout (light region, BTL; black spot, BTB; red spot, BTR), each performed in triplicate. A no-template-control (NTC) was used to check for potential external contamination. Standard curve analysis was performed for reference and target genes to assess amplification efficiency, which was comparable between genes. The conditions for all reactions were 50 °C for 2 min, 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min. At the end of every run, melting curve analysis was performed in order to confirm a unique amplicon reaction.

Expression of ten candidate genes with the highest differential expression between the two species was additionally tested in differently pigmented skin regions of hybrid samples using the protocol described above.

Differential gene expression results were calculated using the Pfaffl method [66], normalised to the geometric mean of the above-mentioned reference genes and with the use of amplicon-specific efficiency. Student’s t-test was performed to assess the statistical significance of differential expression between samples, and results with p < 0.05 were considered statistically significant.

Quantification of pigment patterns

Pigment patterns of hybrid animals and genetically pure brown and marble trout were quantified and analysed as described previously [33]. Briefly, colour pattern complexity, overall tone and the proportion of red spot regions were quantified from images using in-house Python scripts and the OpenCV library. Pattern simplicity score (PSS) is defined as the area-weighted mean isoperimetric quotient of the contours extracted from each image:

$$ PSS=\sum \limits_i{w}_i\ {Q}_i, $$

where \( {Q}_i=\frac{4\uppi {S}_i}{L_i^2} \) is the isoperimetric quotient (or circularity) of each contour, \( {w}_i=\frac{S_i}{\sum \limits_i{S}_i} \) is the area weight, and Si and Li are, respectively, the surface area and the perimeter of each contour. The overall colour tone of a pattern was defined and calculated as the proportion of non-pigmented area (white pixels in a binarized image). The proportion of red spot regions was calculated from the average proportion of reddish pixels in each image.



brown trout black spot skin region


brown trout light skin region


brown trout red spot skin region


Cell adhesion molecules


complementary DNA


differentially expressed genes


Deoxyribonucleic acid


gene ontology


Kyoto Encyclopedia of Genes and Genomes


messenger ribonucleic acid


marble trout dark skin region


marble trout light skin region




polymerase chain reaction


Pattern simplicity score


quantitative real-time PCR


RNA integrity value


ribonucleic acid


ribonucleic acid sequencing


single nucleotid polymorphisms

type L iridophores:

large iridophores

type S iridophores:

small iridophores


whole genome duplication


  1. 1.

    Hubbard JK, Uy JAC, Hauber ME, Hoekstra HE, Safran RJ. Vertebrate pigmentation: from underlying genes to adaptive function. Trends Genet. 2010;26(5):231–9.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Kelsh RN. Genetics and evolution of pigment patterns in fish. Pigment Cell Res. 2004;17(4):326–36.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Watanabe M, Kondo S. Changing clothes easily: connexin41.8 regulates skin pattern variation. Pigment Cell Melanoma Res. 2012;25(3):326–30.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Fujii R. The regulation of motile activity in fish chromatophores. Pigment Cell Res. 2000;13(5):300–19.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Burton D. The physiology of flatfish chromatophores. Microsc Res Tech. 2002;58(6):481–7.

    PubMed  Article  Google Scholar 

  6. 6.

    Bagnara TJ, Matsumoto J. Comparative anatomy and physiology of pigment cells in nonmammalian tissues. In: J. J. Nordlund, R. E. Boissy, V. J. Hearing, R. A. King WSO and JO, editor. The pigmentary system. Oxford: Wiley-Blackwell; 2007. p. 11–59.

  7. 7.

    Braasch I, Schartl M, Volff JN. Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol Biol. 2007;7:1–18.

    Article  CAS  Google Scholar 

  8. 8.

    Ceinos RM, Guillot R, Kelsh RN, Cerdá-Reverter JM, Rotllant J. Pigment patterns in adult fish result from superimposition of two largely independent pigmentation mechanisms. Pigment Cell Melanoma Res. 2015;28(2):196–209.

    CAS  PubMed  Article  Google Scholar 

  9. 9.

    Kelsh RN, Schmid B, Eisen JS. Genetic analysis melanophore development in zebrafish embryos. Dev Biol. 2000;225(2):277–93.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Irion U, Singh AP, Nuesslein-Volhard C. The Developmental Genetics of Vertebrate Color Pattern Formation: Lessons from Zebrafish. In: Wassarman P, editor. Essays on developmental Biology, PT B. San Diego: Elsevier Academic Press Inc; 2016. p. 141. (Current Topics in Developmental Biology; vol. 117).

    Google Scholar 

  11. 11.

    Kelsh RN, Brand M, Jiang YJ, Heisenberg CP, Lin S, Haffter P, et al. Zebrafish pigmentation mutations and the processes of neural crest development. Development. 1996;123:369–89.

    CAS  PubMed  Google Scholar 

  12. 12.

    Parichy DM. Homology and the evolution of novelty during Danio adult pigment pattern development. J Exp Zool Part B Mol Dev Evol. 2007;308B(5):578–90.

    Article  Google Scholar 

  13. 13.

    Nüsslein-Volhard C, Singh AP. How fish color their skin: a paradigm for development and evolution of adult patterns: multipotency, plasticity, and cell competition regulate proliferation and spreading of pigment cells in zebrafish coloration. BioEssays. 2017;39(3):1–11.

    Article  Google Scholar 

  14. 14.

    Patterson LB, Bain EJ, Parichy DM. Pigment cell interactions and differential xanthophore recruitment underlying zebrafish stripe reiteration and Danio pattern evolution. Nat Commun. 2014;5:5299.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Singh AP, Nüsslein-Volhard C. Zebrafish stripes as a model for vertebrate colour pattern formation. Curr Biol. 2015;25(2):R81–92.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Maderspacher F, Nüsslein-Volhard C. Formation of the adult pigment pattern in zebrafish requires leopard and obelix dependent cell interactions. Development. 2003;130(15):3447–57.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Frohnhofer HG, Krauss J, Maischein H-M, Nusslein-Volhard C. Iridophores and their interactions with other chromatophores are required for stripe formation in zebrafish. Development. 2013;140(14):2997–3007.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Patterson LB, Parichy DM. Interactions with iridophores and the tissue environment required for patterning melanophores and xanthophores during zebrafish adult pigment stripe formation. PLoS Genet. 2013;9(5):1–14.

    Article  CAS  Google Scholar 

  19. 19.

    Nakamasu A, Takahashi G, Kanbe A, Kondo S. Interactions between zebrafish pigment cells responsible for the generation of Turing patterns. Proc Natl Acad Sci U S A. 2009;106(21):8429–34.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Iwashita M, Watanabe M, Ishii M, Chen T, Johnson SL, Kurachi Y, et al. Pigment pattern in jaguar/obelix zebrafish is caused by a Kir7.1 mutation: implications for the regulation of melanosome movement. PLoS Genet. 2006;2(11):1861–70.

    CAS  Article  Google Scholar 

  21. 21.

    Fadeev A, Krauss J, Frohnhöfer HG, Irion U, Nüsslein-Volhard C. Tight junction protein 1a regulates pigment cell organisation during zebrafish colour patterning. Elife. 2015;4:1–25.

    Article  CAS  Google Scholar 

  22. 22.

    Watanabe M, Iwashita M, Ishii M, Kurachi Y, Kawakami A, Kondo S, et al. Spot pattern of leopard Danio is caused by mutation in the zebrafish connexin41.8 gene. EMBO Rep. 2006;7(9):893–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Irion U, Frohnhöfer HG, Krauss J, Çolak Champollion T, Maischein H-M, Geiger-Rudolph S, et al. Gap junctions composed of connexins 41.8 and 39.4 are essential for colour pattern formation in zebrafish. Elife. 2014;3:1–16.

    Article  Google Scholar 

  24. 24.

    Adams DS, Levin M. Endogenous voltage gradients as mediators of cell-cell communication: strategies for investigating bioelectrical signals during pattern formation. Cell Tissue Res. 2013;352(1):95–122.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Boettiger A, Ermentrout B, Oster G. The neural origins of shell structure and pattern in aquatic mollusks. Proc Natl Acad Sci. 2009;106(16):6837–42.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Oster F. Lateral inhibition models of developmental processes. Math Comput Model. 1988;90(9):265–86.

    Google Scholar 

  27. 27.

    Allendorf FW, Thorgaard GH. Tetraploidy and the evolution of salmonid fishes. In: Turner BJ, editor. Evolutionary genetics of fishes. Boston: Springer US; 1984. p. 1–53.

    Google Scholar 

  28. 28.

    Berthelot C, Brunet F, Chalopin D, Juanchich A, Bernard M, Noël B, et al. The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. Nat Commun. 2014;5:1–10.

    Article  Google Scholar 

  29. 29.

    Braasch I, Brunet F, Volff JN, Schartl M. Pigmentation pathway evolution after whole-genome duplication in fish. Genome Biol Evol. 2009;1(0):479–93.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. 30.

    Braasch I, Liedtke D, Volff JN, Schartl M. Pigmentary function and evolution of tyrp1 gene duplicates in fish. Pigment Cell Melanoma Res. 2009;22(6):839–50.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Pustovrh G, Snoj A, Bajec SS. Molecular phylogeny of Salmo of the western Balkans, based upon multiple nuclear loci. Genet Sel Evol. 2014;46:1–12.

    Article  Google Scholar 

  32. 32.

    Skaala O, Solberg G. Biochemical genetic variability and taxonomy of a marmorated salmonid in river Otra, Norway. Nord J Freshw Res. 1997;73:3–12.

    Google Scholar 

  33. 33.

    Miyazawa S, Okamoto M, Kondo S. Blending of animal colour patterns by hybridization. Nat Commun. 2010;1(5):1–6.

    PubMed Central  Google Scholar 

  34. 34.

    Djurdjevič I, Kreft ME, Sušnik Bajec S. Comparison of pigment cell ultrastructure and organisation in the dermis of marble trout and brown trout, and first description of erythrophore ultrastructure in salmonids. J Anat. 2015;227(5):583–95.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  35. 35.

    Sivka U, Snoj A, Palandačić A, Sušnik Bajec S. Identification of candidate genes involved in marble color pattern formation in genus Salmo. Comp Biochem Physiol - Part D Genomics Proteomics. 2013;8(3):244–9.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Lien S, Koop BF, Sandve SR, Miller JR, Kent MP, Nome T, et al. The Atlantic salmon genome provides insights into rediploidization. Nature. 2016;533(7602):200–5.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10(1):421.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  38. 38.

    Consortium TU. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2018;47(D1):D506–15.

    Article  CAS  Google Scholar 

  39. 39.

    Hou L, Arnheiter H, Pavan WJ. Interspecies difference in the regulation of melanocyte development by SOX10 and MITF. Proc Natl Acad Sci U S A. 2006;103(24):9081–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Bishop AJR, Kosaras B, Hollander MC, Fornace A, Sidman RL, Schiestl RH. p21 controls patterning but not homologous recombination in RPE development. DNA Repair (Amst). 2006;5(1):111–20.

    CAS  Article  Google Scholar 

  41. 41.

    Besson A, Gurian-West M, Schmidt A, Hall A, Roberts JM. Activation p27 Kip1 modulates cell migration through the regulation of RhoA activation. Genes Dev. 2004;18(8):862–76.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Eller MS, Ostrom K, Gilchrest BA. DNA damage enhances melanogenesis. Proc Natl Acad Sci U S A. 1996;93(3):1087–92.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Šestáková B, Ondrušová L, Vachtenheim J. Cell cycle inhibitor p21/ WAF1/ CIP1 as a cofactor of MITF expression in melanoma cells. Pigment Cell Melanoma Res. 2010;23(2):238–51.

    PubMed  Article  CAS  Google Scholar 

  44. 44.

    Bi MM, Hong S, Zhou HY, Wang HW, Wang LN, Zheng YJ. Chloride channelopathies of CLC-2. Int J Mol Sci. 2013;15(1):218–49.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  45. 45.

    Bellono NW, Oancea EV. Ion transport in pigmentation. Arch Biochem Biophys. 2014;563:35–41.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Miyoshi J, Takai Y. Structural and functional associations of apical junctions with cytoskeleton. Biochim Biophys Acta Biomembr. 2008;1778(3):670–91.

    CAS  Article  Google Scholar 

  47. 47.

    Bendris N, Lemmers B, Blanchard JM. Cell cycle, cytoskeleton dynamics and beyond: the many functions of cyclins and CDK inhibitors. Cell Cycle. 2015;14(12):1786–98.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Sells MA, Knaus UG, Bagrodia S, Ambrose DM, Bokoch GM, Chernoff J. Human p21-activated kinase (Pak1) regulates actin organization in mammalian cells. Curr Biol. 1997;7(3):202–10.

    CAS  PubMed  Article  Google Scholar 

  49. 49.

    Goodenough DA, Paul DL. Gap junctions. Cold Spring Harb Perspect Biol. 2009;1(1):–a002576.

  50. 50.

    Watanabe M, Watanabe D, Kondo S. Polyamine sensitivity of gap junctions is required for skin pattern formation in zebrafish. Sci Rep. 2012;2:1–5.

    Google Scholar 

  51. 51.

    Greenwood AK, Cech JN, Peichel CL. Molecular and developmental contributions to divergent pigment patterns in marine and freshwater sticklebacks. Evol Dev. 2012;14(4):351–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Brooks AN, Turkarslan S, Beer KD, Lo FY, Baliga NS. Adaptation of cells to new environments. Wiley Interdiscip Rev Syst Biol Med. 2011;3(5):491–618.

    Article  CAS  Google Scholar 

  53. 53.

    Isakson BE, Seedorf GJ, Lubman R, Evans W, Boitano S. Cell-cell communication in heterocellular cultures of alveolar epithelial cells. Am J Respir Cell Mol Biol. 2003;29:552–61.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Martins AJ, Narayanan M, Prüstel T, Fixsen B, Park K, Gottschalk RA, et al. Environment tunes propagation of cell-to-cell variation in the human macrophage gene network. Cell Syst. 2017;4(4):379–392.e12.

    CAS  PubMed  Article  Google Scholar 

  55. 55.

    Inaba M, Yamanaka H, Kondo S. Pigment pattern formation by contact-dependent depolarization. Science. 2012;335(6069):677.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Chifflet S, Hernández JA. The plasma membrane potential and the organization of the actin cytoskeleton of epithelial cells. Int J Cell Biol. 2012;2012:1–13.

    Article  CAS  Google Scholar 

  57. 57.

    Sušnik Bajec S, Pustovrh G, Jesenšek D, Snoj A. Population genetic SNP analysis of marble and brown trout in a hybridization zone of the Adriatic watershed in Slovenia. Biol Conserv. 2015;184:239–50.

    Article  Google Scholar 

  58. 58.

    Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  60. 60.

    Tarazona S, Furió-Tarí P, Turrà D, Di PA, Nueda MJ, Ferrer A, et al. Data quality aware analysis of differential expression in RNA-seq with NOISeq R/bioc package. Nucleic Acids Res. 2015;43(21):e140.

    PubMed  PubMed Central  Google Scholar 

  61. 61.

    Tarazona S, García-Alcalde F, Dopazo J, Ferrer A, Conesa A. Differential expression in RNA-seq: a matter of depth. Genome Res. 2011;21(12):2213–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    The Gene Ontology Consortium. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2018;47(D1):D330–8.

    PubMed Central  Article  CAS  Google Scholar 

  64. 64.

    Aoki-Kinoshita KF, Kanehisa M. Gene annotation and pathway mapping in KEGG. Methods Mol Biol. 2007;396:71–91.

  65. 65.

    Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):45e–45.

    Article  Google Scholar 

Download references


The authors would like to thank Tolmin Angling Club and Bled Fishing Club for providing samples. Many thanks go to Iain F. Wilson for English revision of the manuscript.


The authors acknowledge the financial support of the Slovenian Research Agency (Research core funding No. P4–0220). This research formed part of ID’s PhD, financed (Postgraduate research funding) by the Slovenian Research Agency. The funding agency had no role in study design, data collection and analysis, in the interpretation of data, or in writing the manuscript.

Availability of data and materials

The datasets generated and analysed during the current study are available in the SRA repository, within the article’s additional files or available from the corresponding author (see the text for more information).

Author information




ID and SSB conceived the project. ID prepared the samples for RNA sequencing and performed qPCR experiments. TF performed the bioinformatics analysis of the acquired NGS data. SM analysed the pattern variation of hybrids. ID and SSB wrote the initial draft of the manuscript, which was critically revised and approved by all co-authors.

Corresponding author

Correspondence to Simona Sušnik Bajec.

Ethics declarations

Ethics approval and consent to participate

All methods described were carried out in accordance with relevant guidelines and regulations. Fish skin sampling were approved by the Ministry of Agriculture and Environment, Slovenia, under decision letter number U34401–60/2013/4.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Summary results of transcriptome sequencing and mapping. (PDF 21 kb)

Additional file 2:

Table S2. 50 most expressed transcripts in marble trout. (PDF 51 kb)

Additional file 3:

Table S3. 50 most expressed transcripts in brown trout. (PDF 61 kb)

Additional file 4:

Table S4. 100 transcripts with the highest fold change between marble and brown trout samples. (PDF 108 kb)

Additional file 5:

Table S5. 100 transcripts with the highest fold change between brown and marble trout samples. (PDF 111 kb)

Additional file 6:

Figure S6. Gene expression patterns of four candidate genes in differently pigmented regions in the skin of hybrids. Red, red spot; dark brown, dark region; light brown, light region. (PDF 345 kb)

Additional file 7:

Table S7.. Level of brown trout introgression in hybrid individuals. (PDF 417 kb)

Additional file 8:

Figure S8. Correlation between the expression of gja5 gene in the red spots of hybrids and overall colour tone of the skin in hybrids. (PDF 315 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Djurdjevič, I., Furmanek, T., Miyazawa, S. et al. Comparative transcriptome analysis of trout skin pigment cells. BMC Genomics 20, 359 (2019).

Download citation


  • Salmo marmoratus
  • Salmo trutta
  • Pigment cell
  • DEG
  • RNA-seq
  • qPCR