Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Comparative Transcriptome Analysis Reveals the Genetic Basis of Skin Color Variation in Common Carp

  • Yanliang Jiang,

    Affiliation CAFS Key Laboratory of Aquatic Genomics and Beijing Key Laboratory of Fishery Biotechnology, Centre for Applied Aquatic Genomics, Chinese Academy of Fishery Sciences, Beijing, China

  • Songhao Zhang,

    Affiliation CAFS Key Laboratory of Aquatic Genomics and Beijing Key Laboratory of Fishery Biotechnology, Centre for Applied Aquatic Genomics, Chinese Academy of Fishery Sciences, Beijing, China

  • Jian Xu,

    Affiliation CAFS Key Laboratory of Aquatic Genomics and Beijing Key Laboratory of Fishery Biotechnology, Centre for Applied Aquatic Genomics, Chinese Academy of Fishery Sciences, Beijing, China

  • Jianxin Feng,

    Affiliation Henan Academy of Fishery Sciences, Zhengzhou, Henan, China

  • Shahid Mahboob,

    Affiliation Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia

  • Khalid A. Al-Ghanim,

    Affiliation Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia

  • Xiaowen Sun,

    Affiliation CAFS Key Laboratory of Aquatic Genomics and Beijing Key Laboratory of Fishery Biotechnology, Centre for Applied Aquatic Genomics, Chinese Academy of Fishery Sciences, Beijing, China

  • Peng Xu

    xupeng@cafs.ac.cn

    Affiliations CAFS Key Laboratory of Aquatic Genomics and Beijing Key Laboratory of Fishery Biotechnology, Centre for Applied Aquatic Genomics, Chinese Academy of Fishery Sciences, Beijing, China, Department of Zoology, College of Science, King Saud University, Riyadh, Saudi Arabia

Abstract

Background

The common carp is an important aquaculture species that is widely distributed across the world. During the long history of carp domestication, numerous carp strains with diverse skin colors have been established. Skin color is used as a visual criterion to determine the market value of carp. However, the genetic basis of common carp skin color has not been extensively studied.

Methodology/Principal Findings

In this study, we performed Illumina sequencing on two common carp strains: the reddish Xingguo red carp and the brownish-black Yellow River carp. A total of 435,348,868 reads were generated, resulting in 198,781 assembled contigs that were used as reference sequences. Comparisons of skin transcriptome files revealed 2,012 unigenes with significantly different expression in the two common carp strains, including 874 genes that were up-regulated in Xingguo red carp and 1,138 genes that were up-regulated in Yellow River carp. The expression patterns of 20 randomly selected differentially expressed genes were validated using quantitative RT-PCR. Gene pathway analysis of the differentially expressed genes indicated that melanin biosynthesis, along with the Wnt and MAPK signaling pathways, is highly likely to affect the skin pigmentation process. Several key genes involved in the skin pigmentation process, including TYRP1, SILV, ASIP and xCT, showed significant differences in their expression patterns between the two strains.

Conclusions

In this study, we conducted a comparative transcriptome analysis of Xingguo red carp and Yellow River carp skins, and we detected key genes involved in the common carp skin pigmentation process. We propose that common carp skin pigmentation depends upon at least three pathways. Understanding fish skin color genetics will facilitate future molecular selection of the fish skin colors with high market values.

Introduction

Coloration is one of the most diverse phenotypic traits in vertebrates, and it exerts multiple adaptive functions, such as species identification, camouflage, warning or threatening of predators, photoprotection, thermoregulation and photoreception [1]. Skin coloration is the result of diverse pigments synthesized by pigment cells or chromatophores, and it is affected by multiple factors, including environmental, nutritional, physiological, or genetic conditions. Among these factors, the most fundamental and important is the genetic basis of skin pigmentation: which genes are likely to be involved and how. Cellular, genetic and genomic approaches have been widely adopted to answer this question.

In humans and mammals, the pigment melanin is the primary determinant of skin color. The availability of the complete human genome sequence and of adequate genomic resources, along with genome-wide association studies (GWAS), has provided insight into the genetic basis of the pigmentation process. In contrast to mammals, which possess only one type of pigment cell (the melanocyte), and amphibians and reptiles, which possess xanthophores, erythrophores and reflecting iridophores, teleost fish possess up to six different pigment cells, including melanophores, xanthophores, erythrophores, iridophores, leucophores and cyanophores [2]. A diversity of pigment cells, associated with a series of cellular, physiological, genetic and environmental factors, makes fish skin pigmentation a complicated biological process. Extensive studies have been conducted on model fish species such as zebrafish and medaka in an effort to unravel the genetic basis of fish skin pigmentation, and dozens of genes have been reported to be involved in the pigmentation process, such as matp, oca4, sox10, kit, ednrb, slc24a5 and many others, through collecting and identifying the pigmentation mutations [3][5]. However, few genetic skin color studies have been conducted on non-model fish species.

The common carp (Cyprinus carpio), a freshwater fish that is especially widespread in Europe and Asia, was domesticated more than 2,000 years ago. Over its long history of domestication, the common carp has been introduced into various environments worldwide, resulting in hundreds of strains or varieties that display rich biodiversity, genetic polymorphisms and diverse skin colors. Skin color is an important economic trait for common carp, as it acts as an important criterion for visually determining quality and market value. Several studies have provided initial insights into the genetic basis of skin coloration in common carp. For instance, David et al. reported that the MC1R gene was associated with the development of black pigmentation in the ornamental Koi common carp [6]. A very recent study reported that 80% of identified zebrafish pigmentation genes were also present in the Oujiang color common carp strain [7]. In addition to the direct effects of genes, other regulators, such as microRNAs, play crucial roles in common carp skin coloration by regulating the expression of downstream pigmentation genes [8]. However, some details of the genetic mechanisms underlying common carp skin pigmentation are not well understood, such as the interactions among pigmentation genes and the genetic regulation underlying synthesis of different pigments.

The Xingguo red carp (Cyprinus carpio var. xingguonensis, XGC) and the Yellow River carp (Cyprinus carpio haematopterus Temminck et Schlegel, YRC) are two traditional domestic strains of common carp in China. XGC, a regional strain from Xingguo County in the Jiangxi Province in Southern China, is known for its red color and has approximately 1,300 years of cultural history [9]. YRC, a brownish-black common carp, was originally cultivated along the Yellow River basin thousands of years ago and is now a prominent strain in Northern China. It possesses advantageous traits such as strong cold tolerance, high efficiency of food conversion and good meat quality. The distinct skin colors of these two common carp strains make them suitable models for exploring the genetic basis of skin pigmentation. To better understand skin color genetics, we utilized the powerful approach of comparative transcriptome analysis using next-generation sequencing and examined transcript profiles from the skins of the XGC and YRC common carp strains. We obtained candidate genes that may be involved in the skin pigmentation process and identified gene pathways that may regulate the synthesis of different pigments. Understanding the molecular mechanisms of skin pigmentation in common carp will advance our knowledge of skin color genetics in fish and accelerate the molecular selection of fish species with consumer-favored skin colors.

Results and Discussion

Sequencing of short expressed reads from XGC and YRC

Diverse skin colors make fish good genetic models for understanding the skin pigmentation process. Various fish colorations are determined by the density and position of different pigment cells, which is believed to be primarily under genetic control. To better understand fish skin color genetics, we conducted a comparative transcriptome analysis between two common carp strains, XGC and YRC, using next-generation sequencing. First, we generated reference sequences for the subsequent analysis of skin genes that are differentially expressed between the two carp strains. Six tissues, including brain, blood, gill, head kidney, muscle and skin were collected and deep-sequenced using Illumina HiSeq 2000. As shown in Table 1, a total of 435 million paired-end reads were generated, of which 211 million were from YRC and 224 million were from XGC. The number of reads generated from each tissue ranged from 27.5 million to 36.9 million, with 2 outliers of 50.5 million and 46.3 million from XGC muscle and blood, respectively. After the removal of ambiguous nucleotides, low-quality sequences (Q<20) and short reads (length<15 bp), a total of 422 million clean reads (97%) were selected for further analysis.

thumbnail
Table 1. Summary of the raw sequencing data from YRC and XGC.

https://doi.org/10.1371/journal.pone.0108200.t001

Reference sequence assembly and annotation

All of the clean reads were pooled and de novo assembled to generate reference sequences using the Trinity assembler [10]. After the removal of sequence redundancy using CD-hit software [11], a total of 198,781 contigs, with a minimum length of 200 bp, a maximum length of 26,217 bp and an N50 of 1,970 bp, were generated as the reference sequences for subsequent analysis (Table 2). There were 43,310 contigs longer than 1,000 bp. To assess the quality of the sequencing and de novo assembly, all of the clean reads were mapped back to the assembled sequences. As shown in Table 1, the mapping ratio ranged from 77% to 90% with an average of 83%, indicating a high-quality sequence.

Gene prediction was performed on the assembled contigs using BlastX searches against three protein databases, including the NCBI non-redundant (nr) database, the UniProt database and the Ensembl zebrafish protein database, with an E-value cutoff of 1e−10. There were 62,343, 57,083, and 69,634 assembled contigs with significant hits against nr, UniProt and zebrafish, respectively (Table 2). Cumulatively, 73,928 assembled contigs had at least one significant hit against at least one of the three databases, allowing for the prediction of 20,028 unique genes.

Gene ontology (GO) annotation was then performed with these 73,928 unique gene-containing contigs using Blast2GO [12]. Of these, 58,841 contigs, corresponding to 16,849 unique genes, were assigned to at least one GO term (Table 2). As shown in Figure 1, a total of 38 GO terms were assigned, including 10 (26.3%) cellular component terms, 11 (28.9%) molecular function terms and 17 (60.7%) biological process terms. From the GO category of molecular function, binding was the most predominant term, accounting for 63.4% of the sequences annotated in that term, and it was followed by catalytic activity. In the biological processes category, cellular process was the most predominant term (20.2% of the sequences), and it was followed by single-organism process (16.6%) and metabolic process (14.9%).

thumbnail
Figure 1. Distribution of the most common GO term categories.

https://doi.org/10.1371/journal.pone.0108200.g001

Identification of differentially expressed genes in XGC skin compared with YRC skin

Based on the criteria that |fold-change|≧2 and p-value≦0.05, a total of 4,367 of the 198,781 assembled contigs showed significantly different expression in XGC skin compared with YRC skin. Differentially expressed contigs represented 2,012 unique genes, of which 874 genes were up-regulated in XGC and 1,138 genes were down-regulated (Figure 2, Table S1). To validate the differentially expressed genes identified by comparative transcriptome analysis, we randomly selected 20 representative genes for qRT-PCR confirmation of differential expression. The gene expression patterns in the result of qRT-PCR were compared with the data obtained from the comparative transcriptome analysis. As shown in Figure 3, the qRT-PCR expression patterns of 19 out of the 20 randomly selected differentially expressed genes were in agreement with the results of the comparative transcriptome analysis with only slight differences in expression levels, indicating that there was no consistent bias in the expression patterns (i.e., in the direction of the differential expression or in the degree of fold change) for either method. Melting-curve analysis showed that a single product was amplified for all tested genes, indicating that the reference assembly was largely accurate and that it did not contain a large number of chimeric transcripts.

thumbnail
Figure 2. Gene expression in the skin of two common carp strains.

The left panel shows an M-A plot, where the y-axis represents the logarithm of fold change and the x-axis represents the logarithm of transcript counts. The right panel shows a Venn diagram with the number of differentially expressed genes in the two carp strains.

https://doi.org/10.1371/journal.pone.0108200.g002

thumbnail
Figure 3. Comparison of gene expression patterns obtained using comparative transcriptome analysis and qRT-PCR.

Fold changes are expressed as the ratio of gene expression between XGC and YRC after normalization to β-actin.

https://doi.org/10.1371/journal.pone.0108200.g003

After carefully investigating the differentially expressed gene list from the skin tissues of XGC and YRC, we identified many putative pigmentation genes that may be associated with melanin synthesis, such as those encoding tyrosinase-related protein 1 (TYRP1), premelanosome protein (SILV/PMEL), cysteine/glutamate transporter (xCT/SLC7a11) and agouti signaling protein (ASIP). Melanin is the major pigment present in vertebrates [13]. There are two types of melanin, eumelanin and pheomelanin. Eumelanin corresponds to a brown/black color, while pheomelanin corresponds to a red/yellow color. As their names suggest, the skin color of XGC is red, while YRC is brownish-black, indicating that the amount and density of eumelanin are much lower in XGC skin than in YRC skin. Conversely, we hypothesized that the amount and density of pheomelanin might be much higher in XGC skin than in YRC skin. The observed differential expression patterns of melanin synthesis pathway genes endorsed this hypothesis. As shown in Table 3, the expression level of TYRP1 was 5.7-fold up-regulated in YRC skin compared with XGC, suggesting that TYRP1 is one key gene that contributes to brownish-black coloring in common carp by acting in the eumelanin synthesis pathway. TYRP1 is expressed specifically in melanocytes, and it plays a crucial role in pigmentation and color patterning by affecting melanin synthesis, stabilizing tyrosinase protein, modulating tyrosinase catalytic activity, maintaining melanosome structure, and affecting melanocyte proliferation and cell death [14][16]. TYRP1 was the first cloned pigmentation gene, and it was later mapped to the brown locus in mouse [17]. Mutations in TYRP1 result in failure to form a black coat color, as observed in mouse [18], dog [19], cat [20], cow [21], and sheep [22]. In human cells, TYRP1 is only detectable in cells containing eumelanin [23], and it is involved in the pigmentation differences among human populations [24]. These studies suggest essential roles for the TYRP1 gene in eumelanogenesis. Similar roles for TYRP1 have been reported in teleost fish. Knockdown of the TYRP1 gene in zebrafish and medaka revealed that black eumelanin formation essentially relies on the presence of TYRP1 [25]. Our results confirmed a role for TYRP1 in the skin coloration of the common carp.

thumbnail
Table 3. Detailed information about the differentially expressed genes involved in each of the three pigmentation-related pathways.

https://doi.org/10.1371/journal.pone.0108200.t003

Another significant eumelanin synthesis pathway gene we identified was SILV/PMEL. SILV encodes premelanosome protein, which is extensively expressed in pigment cells [26] and can catalyze the conversion of indole-5,6-quinone carboxylic acid into eumelanin [27]. Mutations in SILV causing pigmentation phenotypes have been reported in a number of vertebrate species. Inactivating SILV in mice led to a substantial reduction in eumelanin, indicating that SILV plays a critical role in maintaining efficient epidermal pigmentation [26]. In mice, SILV is required for normal melanosome development in skin melanocytes, choroid melanocytes and retinal pigment epithelium cells [26]. In the domesticated chicken, SILV mutations inhibited the production of all eumelanin in plumage and skin [28]. Similarly, in horses, SILV mutations caused a dilution of black pigment in the mane and tail [29]. Consistent with all these studies, our results showed that the expression level of SILV in YRC was significantly higher than in XGC, exhibiting a 5.2-fold up-regulation (Table 3). This result further confirmed the role of SILV in common carp eumelanogenesis.

Several other genes in our differentially expressed gene list are involved in the production of pheomelanin, the most significant of which are ASIP and xCT. In our results (Table 3), both ASIP and xCT were significantly up-regulated in XGC (by 3.2-fold and 3.5-fold, respectively), suggesting a significant effect on the production of pheomelanin in common carp. The functions of ASIP and xCT in pheomelanin synthesis have been widely studied in other species. ASIP can regulate pigmentation by antagonizing the binding of a-MSH to Mc1r, thus switching melanin synthesis from eumelanin to pheomelanin [30]. The ASIP protein also down-regulates genes such as TYRP1, DCT, and TYR to inhibit the synthesis of eumelanin [31], which is in agreement with our result that ASIP gene expression is up-regulated while TYRP1 gene expression is down-regulated in XGC skin. The important role of ASIP in determining skin or coat color has been demonstrated in human [32], mouse [33], fox [34], rat [35], pig [36] and other species. xCT (SLC7a11), a cysteine/glutamate exchanger, mediates the cellular uptake of cysteine [37]. Cysteine is an important component of pheomelanin; therefore, xCT directly affects the pheomelanin synthesis pathway [38]. Loss of xCT expression caused marked inhibition of pheomelanogenesis [38]. Conversely, high-level expression of xCT contributed to the production of pheomelanin, as indicated by our results.

Enrichment and pathway analysis

We attempted to categorize the 2,012 differentially expressed genes based on their likely functions by using gene annotation, GO term enrichment analysis, and KEGG pathway analysis. Using the Blast2GO software, all the differentially expressed genes were classified into different cellular, biological and functional gene ontologies. Together with the GO annotation of the reference sequences, the web-based program WEGO [39] revealed 17 significantly over-represented GO terms at the 2nd GO level that were enriched in XGC skin relative to YRC skin. The top 10 of these are shown in Table 4, and they include pigmentation (GO:0043473), binding (GO:0005488), enzyme regulator activity (GO:0030234), biological regulation (GO:0065007), multicellular organism process (GO:0032501), multi-organism process (GO:0051704), developmental process (GO:0032502), cellular process (GO:0009987), localization (GO:0051179), and extracellular region part (GO:0044421). The primary interest of this study was to detect expression signatures indicative of fish skin pigmentation; therefore, the 363 genes in the pigmentation GO category were considered to be most informative for further pathway analysis.

thumbnail
Table 4. The top ten GO terms in the GO enrichment results of genes with significantly different expression in XGC skin compared to YRC skin.

https://doi.org/10.1371/journal.pone.0108200.t004

As described in a pathway analysis of non-model fish species [40], downstream pathway analysis was conducted using KEGG pathway analysis combined with manual literature searches. We mainly focused on pigmentation-related pathways to reveal the mechanism of skin color variation in the common carp. We included 1) melanin biogenesis; 2) the Wnt signaling pathway; and 3) the MAPK signaling pathway. Table 3 lists the key differentially expressed genes involved in each pathway.

Melanin biogenesis is under complex controls by multiple agents via pathways. Slominski et al. [41] elucidated that myriad hormonal factors are involved in regulating melanogenesis, such as melanocortins (MSH), ß-endorphin, endothelins, histamines, eicosanoids, catecholamines, c-kit ligand estrogens, androgens, vitamin D, serotonin, melatonin, dopamine, acetylcholine, agouti proteins and their receptors. Among those hormonal regulators, MSH with its receptor Mc1r is the most important positive regulators, whereas agouti proteins stand out among the negative regulators [41]. Besides hormonal regulators, nutritional regulators such as aromatic amino acid L-tyrosine and L-DOPA play critical roles in melanogenesis as well [42]. The molecular and genetic basis of pigmentation including the structure of tyrosinase, the development and differentiation of pigment cell, is largely conserved between mammals and teleost [41], [43]. The putative gene pathways involved in the common carp skin pigmentation process are shown in Figure 4, indicating the critical roles of the hormonal and nutritional regulators in melanogenetic pathways in common carp (Figure 4). Briefly, MSH binds to the G-protein-coupled receptor Mc1r, resulting in up-regulated cAMP levels, which in turn trigger the eumelanin biosynthesis process. The synthesis of eumelanin is then catalyzed successively by Tyr, DCT, TYRP1 and SILV [44] (Figure 4). The switch from eumelanogenesis to pheomelanogenesis during melanin biosynthesis mainly depends on the presence of ASIP, an inverse agonist of Mc1r, which consequently make melanocytes start producing pheomelanin and stop producing eumelanin.

thumbnail
Figure 4. Diagram of putative gene pathways in the common carp skin pigmentation process.

represents the melanin biosynthesis pathway. represents the Wnt signaling pathway. represents the MAPK signaling pathway.

https://doi.org/10.1371/journal.pone.0108200.g004

Signaling pathway such as Wnt and MAPK is implicated in numerous development and physiological process. Several studies have been reported that both Wnt and MAPK signaling pathways take part in the development of pigment cells [45][47]. Intriguingly, our results indicated that both pathways are highly likely to be involved in melanin biosynthesis in common carp. A group of identified differentially expressed genes involved in the Wnt or MAPK signaling pathway have the potential affecting the melanin biosynthesis (Table 3). For instance, WNT5B, a member of the Wnt protein family, is expressed in pigment cells located at the extended posterior edges of the caudal fin in zebrafish [48]. In our results, WNT5B was 2.5-fold up-regulated in YRC, suggesting that WNT5B may be involved in the skin pigmentation process. SFRP5 (secreted frizzled-related protein 5), a member of the SFRP family that modulates Wnt signaling transduction, is highly expressed in the retinal pigment epithelium [49]. This finding, combined with our result that the expression level of SFRP5 is higher in the skin of brownish-black YRC than in red XGC (Figure 3), raises the likelihood that SFRP5 plays an important role in the Wnt signaling pathways involved in eumelanin synthesis.

The melanogenetic pathway in the skin could interact with other regulatory pathways such as nervous system, immune system, or circulatory system, via various organic compounds [50]. For instance, catecholamine, a compound produced in the skin, can be oxidized to form melanin, and also can affects the skin immune functions and systemic functions through communicating with specific cell surface receptors [50]. However, the details of interactions of regulatory pathways in common carp skin remain to be further investigated.

Materials and Methods

Ethics statement

This study was approved by the Animal Care and Use committee of the Centre for Applied Aquatic Genomics at the Chinese Academy of Fishery Sciences.

Fish sampling, RNA isolation and sequencing

As described in our previous study [51], the two strains of common carp used in this study were sampled from distinct breeding populations. XGC individuals were collected from the National Fish Hatchery of Xingguo red carp, Xingguo, Jiangxi province, China. YRC individuals were collected from the Henan Academy of Fishery Sciences, Henan province, China. Six tissues from 18 individuals of each strain were collected; the tissues included brain, blood, gill, head kidney, skin and muscle. All tissues were placed in 2 ml of RNAlater (Qiagen, Hilden, Germany), and kept at −20°C until RNA extraction. RNA isolation was performed using TRIzol (Invitrogen, Carlsbad, CA, USA) with DNase I in accordance with the manufacturer’s instructions. RNA quality was verified using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Equal amounts of high-quality RNA from each tissue were pooled and sent to HudsonAlpha Genomic Services Laboratory (Huntsville, AL, USA) for sequencing on an Illumina HiSeq2000 platform. All data generated were deposited into NIH Short Read Archive with accession number PRJNA254191.

De novo assembly of reference sequences

All raw sequencing reads were trimmed by removing adaptor sequences, ambiguous nucleotides, low quality sequences (Q<20) and short reads (length below 30 bp) using CLC Genomics Workbench (CLC bio, Aarhus, Denmark). The Trinity software was then used to assemble all of the cleaned reads using default parameters. The assembled sequences were then filtered using the CD-hit program to reduce redundancy. The resulting contigs that were larger than 200 bp were considered as the final non-redundant transcripts and were used as the reference sequences.

Transcriptome annotation and ontology

The assembled reference sequences were used as query sequences to search against the NCBI nr database, the UniProt database and the Ensembl zebrafish protein database. Searches were conducted using the BlastX program with an E-value cutoff of 1e−10. Gene ontology was performed by importing the zebrafish Blast results into Blast2GO software. GO terms were then assigned to each sequence automatically. The annotation output was categorized by cellular component, molecular function and biological process.

Differential gene expression analysis

Clean reads from the skin tissues of each common carp strain were first aligned to their respective reference sequences using Bowtie 2 [52]. As described previously [53], RSEM was then used to calculate and estimate gene or isoform abundances. The expression level of each transcript in each sample was then normalized using edgeR [54]. Transcripts with fold change values larger than 2 and p values lower than 0.05 were included in subsequent analyses as the differentially expressed genes.

qRT-PCR analysis

A total of 20 genes with significantly different expression in the two carp strains were selected for validation using qRT-PCR analysis. The house-keeping gene β-actin was used as an internal reference. All primer sequences are listed in Table S2. The first-strand cDNA was synthesized using the SuperScript III RT kit (Invitrogen) according to the manufacturer’s instructions. All cDNA samples were diluted to 100 ng/µl before use in qRT-PCR reactions on an ABI PRISM 7500 Real-Time Detection System (Life Technologies). The amplifications were performed in a total volume of 15 µl and included 7.5 µl of 2X SYBR Green Master Mix reagent, 1 µl of cDNA and 0.3 µl of each primer (10 µM). The thermal cycling profile consisted of an initial denaturation at 95°C for 10 min followed by 40 cycles of denaturation at 95°C for 15 s and annealing/extension at 60°C for 1 min. An additional temperature-ramping step from 95°C to 65°C was used to produce the melting curve. All reactions were conducted in triplicate and included negative controls with no template. Expression differences between XGC and YRC skin were assessed for statistical significance using a randomization test in REST software. The expression levels of genes were normalized to the levels of β-actin in the same sample.

GO enrichment analysis

The web-based program WEGO was used with default parameters for statistical analysis of GO term overrepresentation among the genes with differential expression in the skin of the two carp strains. The study set represented the frequency of GO terms in the differentially expressed genes, while the population set corresponded to the whole reference gene set.

Conclusions

To better understand the genetic mechanisms of skin color variation in common carp, we conducted a comparative analysis of the skin transcriptomes from two common carp strains with distinct skin colors, the red XGC and the brownish-black YRC. We detected 2,012 unique genes that were differentially expressed in the two strains. Further annotation, GO enrichment and pathway analysis indicated that key pigmentation-related genes were involved in at least three pathways that regulate pigment synthesis in the two strains. These results provide us with a valuable basis for understanding of the molecular mechanisms of skin pigmentation in teleost fish. Such understanding will facilitate the genetic selection and breeding of common carp with market-favored colors.

Supporting Information

Table S1.

Detailed information of the differentially expressed genes. The fold change cutoff is 2, and p-value cutoff is 0.05.

https://doi.org/10.1371/journal.pone.0108200.s001

(XLSX)

Table S2.

Primers used for qRT-PCR validation.

https://doi.org/10.1371/journal.pone.0108200.s002

(XLSX)

Author Contributions

Conceived and designed the experiments: PX. Performed the experiments: YJ SZ. Analyzed the data: YJ JX. Contributed reagents/materials/analysis tools: JF. Contributed to the writing of the manuscript: YJ SM KAG. Supervised the common carp genome project: XS.

References

  1. 1. Hubbard JK, Uy JA, Hauber ME, Hoekstra HE, Safran RJ (2010) Vertebrate pigmentation: from underlying genes to adaptive function. Trends Genet 26: 231–239.
  2. 2. Kelsh RN (2004) Genetics and evolution of pigment patterns in fish. Pigment Cell Res 17: 326–336.
  3. 3. Kelsh RN, Brand M, Jiang YJ, Heisenberg CP, Lin S, et al. (1996) Zebrafish pigmentation mutations and the processes of neural crest development. Development 123: 369–389.
  4. 4. Parichy DM (2006) Evolution of danio pigment pattern development. Heredity (Edinb) 97: 200–210.
  5. 5. Kelsh RN, Inoue C, Momoi A, Kondoh H, Furutani-Seiki M, et al. (2004) The Tomita collection of medaka pigmentation mutants as a resource for understanding neural crest cell development. Mech Dev 121: 841–859.
  6. 6. Bar I, Kaddar E, Velan A, David L (2013) Melanocortin receptor 1 and black pigmentation in the Japanese ornamental carp (Cyprinus carpio var. Koi). Front Genet 4: 6.
  7. 7. Wang C, Wachholtz M, Wang J, Liao X, Lu G (2014) Analysis of the skin transcriptome in two oujiang color varieties of common carp. PLoS One 9: e90074.
  8. 8. Yan B, Liu B, Zhu CD, Li KL, Yue LJ, et al. (2013) microRNA regulation of skin pigmentation in fish. J Cell Sci 126: 3401–3408.
  9. 9. Zhou J, Wu Q, Wang Z, Ye Y (2004) Genetic variation analysis within and among six varieties of common carp (Cyprinus carpio L.) in China using microsatellite markers. Genetika 40: 1389–1393.
  10. 10. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, et al. (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 29: 644–652.
  11. 11. Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658–1659.
  12. 12. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, et al. (2008) High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res 36: 3420–3435.
  13. 13. Riley PA (1997) Melanin. Int J Biochem Cell Biol 29: 1235–1239.
  14. 14. Rad HH, Yamashita T, Jin HY, Hirosaki K, Wakamatsu K, et al. (2004) Tyrosinase-related proteins suppress tyrosinase-mediated cell death of melanocytes and melanoma cells. Exp Cell Res 298: 317–328.
  15. 15. Sarangarajan R, Boissy RE (2001) Tyrp1 and oculocutaneous albinism type 3. Pigment Cell Res 14: 437–444.
  16. 16. Kobayashi T, Hearing VJ (2007) Direct interaction of tyrosinase with Tyrp1 to form heterodimeric complexes in vivo. J Cell Sci 120: 4261–4268.
  17. 17. Jackson IJ (1988) A cDNA encoding tyrosinase-related protein maps to the brown locus in mouse. Proc Natl Acad Sci U S A 85: 4392–4396.
  18. 18. Shibahara S, Tomita Y, Yoshizawa M, Shibata K, Tagami H (1992) Identification of mutations in the pigment cell-specific gene located at the brown locus in mouse. Pigment Cell Res Suppl 2: 90–95.
  19. 19. Schmutz SM, Berryere TG, Goldfinch AD (2002) TYRP1 and MC1R genotypes and their effects on coat color in dogs. Mamm Genome 13: 380–387.
  20. 20. Schmidt-Kuntzel A, Eizirik E, O’Brien SJ, Menotti-Raymond M (2005) Tyrosinase and tyrosinase related protein 1 alleles specify domestic cat coat color phenotypes of the albino and brown loci. J Hered 96: 289–301.
  21. 21. Berryere TG, Schmutz SM, Schimpf RJ, Cowan CM, Potter J (2003) TYRP1 is associated with dun coat colour in Dexter cattle or how now brown cow? Anim Genet 34: 169–175.
  22. 22. Gratten J, Beraldi D, Lowder BV, McRae AF, Visscher PM, et al. (2007) Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep. Proc Biol Sci 274: 619–626.
  23. 23. del Marmol V, Ito S, Jackson I, Vachtenheim J, Berr P, et al. (1993) TRP-1 expression correlates with eumelanogenesis in human pigment cells in culture. FEBS Lett 327: 307–310.
  24. 24. Alonso S, Izagirre N, Smith-Zubiaga I, Gardeazabal J, Diaz-Ramon JL, et al. (2008) Complex signatures of selection for the melanogenic loci TYR, TYRP1 and DCT in humans. BMC Evol Biol 8: 74.
  25. 25. Braasch I, Liedtke D, Volff JN, Schartl M (2009) Pigmentary function and evolution of tyrp1 gene duplicates in fish. Pigment Cell Melanoma Res 22: 839–850.
  26. 26. Hellstrom AR, Watt B, Fard SS, Tenza D, Mannstrom P, et al. (2011) Inactivation of Pmel alters melanosome shape but has only a subtle effect on visible pigmentation. PLoS Genet 7: e1002285.
  27. 27. Chakraborty AK, Platt JT, Kim KK, Kwon BS, Bennett DC, et al. (1996) Polymerization of 5,6-dihydroxyindole-2-carboxylic acid to melanin by the pmel 17/silver locus protein. Eur J Biochem 236: 180–188.
  28. 28. Karlsson AC, Kerje S, Hallbook F, Jensen P (2009) The Dominant white mutation in the PMEL17 gene does not cause visual impairment in chickens. Vet Ophthalmol 12: 292–298.
  29. 29. Brunberg E, Andersson L, Cothran G, Sandberg K, Mikko S, et al. (2006) A missense mutation in PMEL17 is associated with the Silver coat color in the horse. BMC Genet 7: 46.
  30. 30. Voisey J, Box NF, van Daal A (2001) A polymorphism study of the human Agouti gene and its association with MC1R. Pigment Cell Res 14: 264–267.
  31. 31. Sakai C, Ollmann M, Kobayashi T, Abdel-Malek Z, Muller J, et al. (1997) Modulation of murine melanocyte function in vitro by agouti signal protein. EMBO J 16: 3544–3552.
  32. 32. Wilson BD, Ollmann MM, Kang L, Stoffel M, Bell GI, et al. (1995) Structure and function of ASP, the human homolog of the mouse agouti gene. Hum Mol Genet 4: 223–230.
  33. 33. Bultman SJ, Michaud EJ, Woychik RP (1992) Molecular characterization of the mouse agouti locus. Cell 71: 1195–1204.
  34. 34. Vage DI, Lu D, Klungland H, Lien S, Adalsteinsson S, et al. (1997) A non-epistatic interaction of agouti and extension in the fox, Vulpes vulpes. Nat Genet 15: 311–315.
  35. 35. Kuramoto T, Nomoto T, Sugimura T, Ushijima T (2001) Cloning of the rat agouti gene and identification of the rat nonagouti mutation. Mamm Genome 12: 469–471.
  36. 36. Leeb T, Doppe A, Kriegesmann B, Brenig B (2000) Genomic structure and nucleotide polymorphisms of the porcine agouti signalling protein gene (ASIP). Anim Genet 31: 335–336.
  37. 37. Lim JC, Donaldson PJ (2011) Focus on molecules: the cystine/glutamate exchanger (System x(c)(−)). Exp Eye Res 92: 162–163.
  38. 38. Chintala S, Li W, Lamoreux ML, Ito S, Wakamatsu K, et al. (2005) Slc7a11 gene controls production of pheomelanin pigment and proliferation of cultured cells. Proc Natl Acad Sci U S A 102: 10964–10969.
  39. 39. Ye J, Fang L, Zheng H, Zhang Y, Chen J, et al. (2006) WEGO: a web tool for plotting GO annotations. Nucleic Acids Res 34: W293–297.
  40. 40. Li C, Zhang Y, Wang R, Lu J, Nandi S, et al. (2012) RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus. Fish Shellfish Immunol 32: 816–827.
  41. 41. Slominski A, Tobin DJ, Shibahara S, Wortsman J (2004) Melanin pigmentation in mammalian skin and its hormonal regulation. Physiol Rev 84: 1155–1228.
  42. 42. Slominski A, Zmijewski MA, Pawelek J (2012) L-tyrosine and L-DOPA as hormone-like regulators of melanocyte functions. Pigment Cell Melanoma Res 25: 14–27.
  43. 43. Braasch I, Schartl M, Volff JN (2007) Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol Biol 7: 74.
  44. 44. Kobayashi T, Urabe K, Winder A, Jimenez-Cervantes C, Imokawa G, et al. (1994) Tyrosinase related protein 1 (TRP1) functions as a DHICA oxidase in melanin biosynthesis. EMBO J 13: 5818–5825.
  45. 45. Stocker KM, Sherman L, Rees S, Ciment G (1991) Basic FGF and TGF-beta 1 influence commitment to melanogenesis in neural crest-derived cells of avian embryos. Development 111: 635–645.
  46. 46. Fujimura N, Taketo MM, Mori M, Korinek V, Kozmik Z (2009) Spatial and temporal regulation of Wnt/beta-catenin signaling is essential for development of the retinal pigment epithelium. Dev Biol 334: 31–45.
  47. 47. Squarzoni P, Parveen F, Zanetti L, Ristoratore F, Spagnuolo A (2011) FGF/MAPK/Ets signaling renders pigment cell precursors competent to respond to Wnt signal by directly controlling Ci-Tcf transcription. Development 138: 1421–1432.
  48. 48. Lin Sh (2010) Wnt5b signaling in zebrafish development and disease. PhD diss. University of Iowa. http://ir.uiowa.edu/etd/2740.
  49. 49. Chang JT, Esumi N, Moore K, Li Y, Zhang S, et al. (1999) Cloning and characterization of a secreted frizzled-related protein that is expressed by the retinal pigment epithelium. Hum Mol Genet 8: 575–583.
  50. 50. Slominski AT, Zmijewski MA, Skobowiat C, Zbytek B, Slominski RM, et al. (2012) Sensing the environment: regulation of local and global homeostasis by the skin’s neuroendocrine system. Adv Anat Embryol Cell Biol 212: v, vii, 1–115.
  51. 51. Xu J, Ji P, Zhao Z, Zhang Y, Feng J, et al. (2012) Genome-wide SNP discovery from transcriptome of four common carp strains. PLoS One 7: e48140.
  52. 52. Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nat Methods 9: 357–359.
  53. 53. Xu J, Li Q, Xu L, Wang S, Jiang Y, et al. (2013) Gene expression changes leading extreme alkaline tolerance in Amur ide (Leuciscus waleckii) inhabiting soda lake. BMC Genomics 14: 682.
  54. 54. Robinson MD, McCarthy DJ, Smyth GK (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140.