Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 7-Day Trial for You or Your Team.

Learn More →

Comprehensive Transcriptome Analysis of Different Skin Colors to Evaluate Genes Related to the Production of Pigment in Celestial Goldfish

Comprehensive Transcriptome Analysis of Different Skin Colors to Evaluate Genes Related to the... Article Comprehensive Transcriptome Analysis of Different Skin Colors to Evaluate Genes Related to the Production of Pigment in Celestial Goldfish Rongni Li, Yansheng Sun, Ran Cui * and Xin Zhang * Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China * Correspondence: ranc20220918@163.com (R.C.); 13901156204@163.com (X.Z.) Simple Summary: Skin color is the result of various pigments synthesized by dyeing groups, while goldfish have a greater variety of pigment cells and a wider range of skin colors. Up to now, several pigment cells have been reported, including melanophores, erythrocytes, xantho-phores, irido- phores, leucophores and cyanophores. Therefore, the key genes associated with sieve pigments will provide clues for the production of goldfish with a diversity of skin tones, benefiting the goldfish industry. In this study, high-throughput sequencing was performed on the back skin tissue of ce- lestial goldfish with different skin colors (white, yellow, brown), and it was found that melanin production, tyrosine metabolism pathway, Wnt signaling pathway, and MAPK signaling pathway may play an important role in pigmentation. The study provides valuable data for screening candi- date genes involved in pigment formation, and the results of this study will lay a foundation for further study on the expression characteristics and gene network analysis of pigment related genes. Abstract: Skin color is an important phenotypic feature of vertebrate fitness under natural condi- tions. Celestial goldfish, a common goldfish breed in China, mainly shows three kinds of skin colors including white, yellow and brown. However, the molecular genetic basis of this phenotype is still unclear. In this study, high-throughput sequencing was carried out on the back skin tissues of ce- lestial goldfish with different skin colors. About 58.46 Gb of original data were generated, filtered and blasted, and 74,297 mRNAs were obtained according to the reference transcriptome. A total of 4653 differentially expressed genes were screened out among the brown, yellow and white groups, Citation: Li, R.; Sun, Y.; Cui, R.; and the expression of melanogenesis related genes in brown goldfish was significantly higher than Zhang, X. Comprehensive the other two groups. There are 19 common differentially expressed genes among three groups, of Transcriptome Analysis of Different which eight genes are related to pigment production, including tyrp1a, slc2a11b, mlana, gch2, Skin Colors to Evaluate Genes loc113060382, loc113079820, loc113068772 and loc113059134. RT-qPCR verified that the expression Related to the Production of Pigment in Celestial Goldfish. patterns of randomly selected differentially expressed transcripts were highly consistent with those Biology 2023, 12, 7. https://doi.org/ obtained by RNA sequencing. GO and KEGG annotation revealed that these differentially ex- 10.3390/biology12010007 pressed genes were mostly enriched in pathways of the production of pigment, including melano- genesis, tyrosine metabolism, Wnt signaling pathway, MAPK signaling pathway etc. These results Academic Editors: Mingxing Chu, indicated that the external characteristics of goldfish are consistent with the analysis results at tran- Ran Di and Yufang Liu scriptome level. The results of this study will lay a foundation for further study on the expression Received: 3 December 2022 characteristics and gene network analysis of pigment related genes. Revised: 17 December 2022 Accepted: 19 December 2022 Keywords: celestial goldfish; skin color; transcriptome; pigment production Published: 21 December 2022 Copyright: © 2022 by the authors. Li- censee MDPI, Basel, Switzerland. 1. Introduction This article is an open access article Skin color is the result of diverse pigments synthesized by chromatophores; different distributed under the terms and con- ditions of the Creative Commons At- kinds of pigment cells are found in various vertebrates, while there are more kinds of tribution (CC BY) license (https://cre- pigment cells and a wider range of skin colors in goldfish. Up to now, several pigment ativecommons.org/licenses/by/4.0/). Biology 2023, 12, 7. https://doi.org/10.3390/biology12010007 www.mdpi.com/journal/biology Biology 2023, 12, 7 2 of 12 cells have been reported, including melanophores, erythrocytes, xanthophores, irido- phores, leucophores and cyanophores [1,2]. Although the production of skin color in fish is complex and varies depending on the breed, feed, and water environment, genes are the most basic and essential element that influences skin pigmentation. RNA-seq is a study of gene structure and function at the overall transcriptome level, based on high-throughput sequencing technology, and reveals the molecular mechanism in specific biological processes [3]. Many studies have been devoted to explaining the ge- netic basis of pigment production by RNA-seq in recent years, and a certain number of important genes and biological pathways have been found, including melanogenesis and tyrosine metabolism. Jiang et al. [4] performed transcriptome analysis of skin color during and after the overwintering of red tilapia, showing that pigmentation-related genes in- clude tyrp1, tyr, pmel, MITF and apoptosis and autophagy-related genes possibly regulate the molecular mechanisms of color change. Gan et al. [5] performed comparative tran- scriptome analysis on four skin colors of goldfish and concluded that the melanogenesis and tyrosine metabolism pathway may affect the skin pigment agglutination process. Luo et al. [6] found that biological pathway was mainly enriched in tyrosine metabolism and melanogenesis through transcriptome analysis of koi fish with white, brown and red skin colors. Chen et al. [7] reported that, consistent with goldfish, the differential expressed genes of various skin colors in Silurus meridionalis is also mainly concentrated in mel- anogenesis, melanin metabolism and melanin agglutination hormone activity pathway. In addition, tyrosinase (TYR) [8–10], solute carrier family 24 member 5 [11,12], GTP_cy- clohydroI (gch) [13,14] and microphthalmia associated transcription factor (MITF) [15] have been reported as important genes involved in melanogenesis. Celestial goldfish, as a unique variety in China, is deeply loved by people with eyes turned upward and color- ful body colors. However, the key genes and biological processes, which are involved in the production of pigment in celestial goldfish, are still unclear. The molecular process of the production of pigment in celestial goldfish has not yet undergone a comprehensive examination. Therefore, clarifying the roles of the differential expression genes can contribute to understanding the changes associated with the pro- duction of pigment in different skin color groups. The purpose of this study is to obtain an overview of transcriptome profiles in three kinds of skin colors including white, yellow and brown in skin tissue. In addition, we identified differentially expressed genes that were possibly involved in the production of pigment and screened out common differen- tial expression genes among three groups. In-depth knowledge of the molecular processes underlying the production of pigment in celestial goldfish is provided in this report. 2. Material and Method 2.1. Animal Sample Collection The celestial goldfish used in the experiment was taken from Xiaotangshan Aquacul- ture Base, Fisheries Research Institute, Beijing Academy of Agricultural and Forestry Sci- ences (Beijing, China). The experimental group used three brown and three yellow gold- fish, and the control group used three white goldfish. All three groups of goldfish came from half-sibs, were raised under the same conditions until they were 4 months old, and then were slaughtered on the same day. The fish body length is 5–6 cm, and the goldfish samples collected in the observation experiment were faded and the body color was sta- bilized. The skin tissues from the backs of nine goldfish were extracted and stored in liq- uid nitrogen or at −80 °C for later use. 2.2. RNA Extraction and Sequencing According to the instructions, skin tissue samples in three groups were taken, and the total RNA of each sample was extracted by RNA extraction kit (Invitrogen, Carlsbad, CA, USA). To analyze the integrity of the extracted total RNA and whether there was DNA contamination, 10 g/L agarose gel was used. The purity of RNA was detected by Biology 2023, 12, 7 3 of 12 nucleic acid protein determination spectrophotometer (IMPLEN, Munich, Bavaria, Ger- many). The Agilent 2100 bioanalyzer was used to accurately detect the integrity and total amount of RNA. According to the instructions, NEB-Next Ultra TM RNA Library Prep Kit (Illumina) was used to build the library. The cDNA was synthesized by magnetic bead enrichment method, and the cDNA end was repaired, amplified and purified to obtain the cDNA library of the skin tissue of celestial goldfish. 2.3. Bioinformatic Analysis The original data were filtered by removing N-containing, jointed and low-quality reads, and the GC content, Q20 and Q30 of clean data were calculated. Carassius auratus reference genome and gene model annotation files were downloaded from the NCBI web- site (https://www.ncbi.nlm.nih.gov/assembly/GCA_003368295.1 (accessed on 1 January 2022), and the parametric analysis of transcriptome sequencing data was performed. HISAT2 [16] was used to align paired-end clean reads to the reference genome. The mapped readings of each skin sample were put together using StringTie [16] with a refer- ence-based method, and then new transcripts were annotated by Pfam [17], GO [18], KEGG [19] and other databases, and the new genes were predicted. Using the Feature- Counts program, the readings mapped to each gene were calculated [20], and the FPKM of each gene (the expected number of every thousand base segments of transcript se- quence segments sequenced per million base pairs) was calculated based on the gene length. Differential expression analysis was performed using DESeq [21], with a default significant difference of two fold changes and a q value of 0.05. 2.4. Gene ontology and kyoto encyclopedia of genes and genomes enrichment analysis Cluster-Profile software was used to analyze the enrichment of differentially ex- pressed genes by GO function and KEGG pathway, and the significant enrichment thresh- old was p value < 0.05. The main biochemical metabolic pathways and signal transduction pathways involved in color change were screened out. 2.5. Validation of Differentially Expressed Transcripts RT-qPCR verified the expression pattern of 10 differentially expressed transcripts randomly selected from the skin tissue of the back. Total RNA of the sample was reverse transcribed into cDNA by prime script TM 1st stand cDNA synthesis kit. β-actin was se- lected as endogenous control gene (Table S7). Three biological replicates and triple reac- tions for each sample were used for all qPCR validations. The relative expression of each gene was calculated by 2-△△Ct method, and the log2 value of multiple expression be- tween the two groups was calculated. The column chart was drawn together with the RNA-seq sequencing results, and the trend was compared. 2.6. Analytical Statistics The means and standard error of the means (SEM) are used to express data. In order to control the error rate, Hochberg and Benjamini methods are used to adjust the p value; differences were deemed statistically significant at p < 0.05. 3. Results 3.1. Blast Analysis of Transcriptome Sequencing The transcriptome profiling of the skin tissue from goldfish of different colors were obtained to evaluate the genes involved in the production of melanin. A total of 58.46 Gb of filtered base was obtained, and the base amount of each sample reached more than 5.96 Gb, and 87.72~92.96% of these readings were specifically compared with the reference ge- nome (https://www.ncbi.nlm.nih.gov/assembly/GCA_003368295.1, accessed on 1 January 2022). The base error rate of each sample sequencing is less than 0. 03%, and Q30 is greater than 93.90%, which indicates that the sequencing data quality is reliable (Table S1). Biology 2023, 12, 7 4 of 12 Celestial goldfish is a non-model species, and its genetic annotation in production is usually not perfect. After the alignment was completed, the annotated transcription re- gions of the genome were searched, and the aligned sequences were assembled, so as to mine new transcripts or new genes in the species. In this study, new transcripts were as- sembled and annotated in Pfam, GO and KEGG databases and 2205 new genes were pre- dicted and annotated (Table S2). 3.2. Differentially Expressed Genes in Goldfish with Different Skin Colors Through the comparison of different skin color groups, the differentially expressed genes are shown (Figure 1a,b). When compared between brown and yellow groups, a total of 1000 differentially expressed genes were screened out, of which 651 were up-regulated and 349 were down-regulated. When comparing between white and brown groups, we screened out 3033 differentially expressed genes, of which 1196 genes were up-regulated and 1837 genes were down-regulated. Consistent with the visual contrast between the brown and white groups, the number of differentially expressed genes in this group is more than other groups through RNA-seq. When comparing between white and yellow groups, a total of 1942 differentially expressed genes were screened out, of which 761 were up-regulated and 1181 were down-regulated. Some DE mRNAs were specifically ex- pressed in the brown group, including: loc113077208, loc113119661, loc113085971, loc113066815; Some DE mRNA was specifically expressed in the white group, including: clic5a, loc113104870, loc113118795; There are also some specific expression of DE mRNAs in the yellow group, including: slc2a11b, asz1, loc113092873. These genes may regulate the production of pigment in goldfish (Table S3). The differentially expressed genes of three groups of goldfish were clustered and drawn into a heat map (Figure 1c). It can be seen that the expression patterns of the brown group and the yellow group are more similar, while the white group and the other two groups are at a far distance. Biology 2023, 12, 7 5 of 12 Figure 1. (a). The different skin colors in the celestial goldfish. Left: white (W); center：yellow (Y); right: brown (B). (b) The bar graph of number of DEGs among three different skin colors. (c) Clus- tering heat map of DEGs in the three skin color groups. (d) DEGs number and Venn diagram of the overlap of the different skin color groups. W: white; Y: yellow; B: brown. 3.3. Functional Analysis of Differentially Expressed Genes GO enrichment analysis was conducted between brown, yellow and white group ac- cording to three categories: biological process, cell components and molecular functions. The most significant 30 entries were selected from the GO enrichment analysis results to draw scatter plots (Figure 2, Table S4). GO enrichment entries in the white and yellow groups are significant in the biological process (BP), including peptide cross-linking (GO:0018149), small molecule biosynthetic process (GO:0044283), protein tyrosine phos- phatase activity (GO:0004725), etc. The GO entries of the brown and yellow groups, white and brown groups were both significantly enriched in the cellular component (CC), mainly including melanosome transport-related: troponin complex (GO:0005861), myofil- ament (GO:0036379), contractile fiber (GO:0043292), etc. GO enrichment concluded that differentially expressed genes were mainly enriched in the pathways related to tyrsine metabolism and melanosome transport. Biology 2023, 12, 7 6 of 12 Figure 2. (a). The scatter figure of GO enrichment in white vs. yellow. (b) The scatter figure of GO enrichment of DEGs in brown vs. yellow. (c) The scatter figure of GO enrichment of DEGs in white vs. brown. Count: number of genes in each pathway. GeneRatio: ratio of the number of target genes divided by the total number of genes in each pathway. The differentially expressed genes of three groups in goldfish were compared to KEGG database for pathway enrichment analysis, and the most significant 20 KEGG path- ways were selected to draw columns (Figure 3, Table S5). According to the results of KEGG enrichment, KEGG pathways among three groups are all enriched to pigment pro- duction related, including Wnt signaling pathway (dre04310), tyrosine metabolism (dre00350), melanogenesis (dre04916) and MAPK signaling pathway (dre04010). These re- sults reflect that the difference of skin color of goldfish may be the result of the regulation of multiple signal pathways, and the different skin color characteristics in this study are closely related to pigment production. Biology 2023, 12, 7 7 of 12 Figure 3. (a). The bar graph of KEGG pathway enrichment of DEGs in white vs. yellow. (b) The bar graph of KEGG pathway enrichment of DEGs in brown vs. yellow. (c) The bar graph of KEGG path- way enrichment of DEGs in white vs. brown. 3.4. Common Expression Gene Analysis According to the overlapping of differential genes among different skin color groups, the results showed that 19 differential genes belonged to the common differential genes among three groups (Figure 1d, Tables 1 and S6), among which eight genes were reported to be related to pigment production and regulation [8,14,22–26], including: tyrp1a, slc2a11b, mlana, gch2, loc113060382, loc113079820, loc113068772 and loc113059134. The re- maining genes mainly regulate cell proliferation, differentiation and apoptosis, including transcription elongation factor a protein 3, translation initiation factor if-2, RAS and EF- hand domain containing, translocation protein Sec62, et al. Table 1. A total of 19 common differentially expressed genes among three groups. FPKM Gene ID Description White Yellow Brown 113116617 0.053316 14.85341 4.939858 uncharacterized 113059134 0.101642 11.21777 3.472121 perilipin-3-like novel.91 0.024732 41.3923 11.47289 -- 113094292 0 22.96619 9.53373 putative defense protein 3 113054403 0 7.702924 1.333329 diacylglycerol O-acyltransferase 2 113060382 0.021323 4.63168 1.063391 retinol dehydrogenase 7-like solute carrier family facilitated glucose 113072589 0 3.77947 0.555086 transporter member 11-like melanoma antigen recognized by T-cells 1- 113081155 0.065602 12.08961 3.294872 like novel.919 0.052276 2.696407 9.566884 -- 113056131 0.017281 4.126237 0.910635 GTP cyclohydrolase 1-like melanoma antigen recognized by T-cells 1- 113079820 0.023412 8.795228 2.98262 like 113116699 18.29243 526.6035 80.16275 translation initiation factor IF-2-like Biology 2023, 12, 7 8 of 12 113044679 0.872271 0.010594 0.200235 uncharacterized 113107421 0.154702 2.39033 0.707736 RAS and EF-hand domain containing transcription elongation factor A protein 3- 113059190 1.056176 0.085419 5.333497 like 113068772 0.148587 0.020619 1.319104 myosin-6-like 113119606 3.652859 1.230851 0.193719 translocation protein SEC62-like C6-dihydroxyindole-2-carboxylic acid oxi- 113106535 0.017673 0.195103 1.180688 dase-like novel.58 0.154461 0.976559 4.569909 -- 3.5. Validation of Differentially Expressed Genes Expression Eight genes were randomly selected from the differentially expressed genes for qPCR verification, and the results showed that the gene expression trend obtained by qPCR was basically consistent with the results of RNA-seq (Figure 4). It indicates that the sequencing results of transcriptome in this study are reliable, and the differentially expressed genes among different groups can be accurately detected. Figure 4. The expression levels of mRNA verified by qPCR. LOC113045861-LOC113043321: left: yel- low group, right: brown group; LOC113075300: left: white group, right: brown group; LOC113040942-LOC113066051: left: yellow group, right: white group. * indicates p < 0.05. 4. Discussion 4.1. Composition Analysis of Pigment Cells in Different Skin Colors Skin color plays an important role in animal courtship, escape from enemies and sig- nal communication, and is one of the most diverse forms of kinvertebrates. Due to diverse color pattern, goldfish could be a perfect model to study the genetic mechanism of pig- mentation. The research on the formation mechanism of fish skin color provides an im- portant theoretical reference for the study of physiology, pathology and genetics of hu- man and other biological pigment cells. Although the production of skin color in fish is complex and varies depending on the breed, feed, and water environment [27,28], genes are the most basic and essential element that influences skin pigmentation. In this study, GO enrichment reflected the different components of pigment cells in different skin col- ors. The skin cells of white goldfish are mainly composed of white pigment cell, while those of yellow goldfish are mainly yellow pigment cells. Studies have reported that that characteristics and differentiation processes of white pigment and yellow pigment cells Biology 2023, 12, 7 9 of 12 are similar in medaka [23], which also explains why there is less difference in cellular composition between white and yellow groups than in comparison with the brown group. Inconsistent with white and yellow groups, there were dendritic melanocytes in the skin of brown groups, and GO enrichment enriched the related pathways with protein tyrosine phosphatase activity. Tyrosinase (TYR) is the rate limiting enzyme for melanin synthesis, and its gene mutation will cause albinism in animals. Since the cloning of human tyrosi- nase gene in 1987 [29], much research has been carried out on the relationship between tyrosine gene and skin color (mainly albinism) in humans and animals [30,31]. Melanin is a polymer synthesized by tyrosine hydroxylation by tyrosinase catalyzed in vivo, which involves three enzymes in the tyrosinase gene family, tyrosinase (TYR), tyrosinase-related protein 1 (TRP1) and tyrosinase-related protein-2 (TRP2). After melanin is produced, it travels from the dendritic tip of melanocytes to keratinocytes. Furthermore, GO enrich- ment also enriched the related pathways with the structure of iridescent cells, actin and enzyme activity. It has been reported that the directional movement of pigment particles is closely related to actin. 4.2. Regulated Genes in Brown Skin Indicated the Molecular Mechanism of the Melanogenesis Pathway At present, most studies on the skin color of goldfish are focused on pigment cell metabolism, pigment deposition and production, especially melanin production and transport, which is consistent with the results in this study. For example, Jiang et al. [32] performed comparative transcriptome analysis on the skin colors of Xingguo Red Carp and Yellow River Carp and also concluded that melanogenesis pathway affected the skin pigment agglutination process in common carp. According to our research, the differen- tially expressed genes were not only enriched in the melanogenesis pathways mentioned above in the existing studies, but MAPK signaling pathway and Wnt signaling pathway were also enriched. Wnt signaling pathways, particularly Wnt1, Wnt3a and Wnt5a, have been reported to play an important role in the development of melanocytes [33–36]. Wnt3a has been reported to reduce the number of neuronal cells and glial cells in quail, while significantly increasing the number of melanocytes [34]. Another report on zebrafish sug- gests that after injection of domain-inactivated Wnt1, the Wnt signaling pathway is blocked, the number of pigment cells is reduced, and neural crest cells are easily differen- tiated into neuronal cells and glial cells [35]. Vertebrate melanocytes almost always origi- nate from the neural crest, and the cells derived from the neural crest undergo obvious differentiation until they are dispersed, and produce peripheral neurons, glial cells and melanocytes at precise locations [36]. Furthermore, phosphorylation of ERK in MAPK sig- naling pathway inhibits the production of melanin in zebrafish, while inhibition of phos- phorylation of ERK can promote the production of melanin in zebrafish by increasing the expression of MITF genes, thereby upregulating downstream genes TYR, TRP1 and TRP2 [37]. The expression of MITF was also able to differentiate pluripotent stem cells into mel- anocytes in medaka [38]. Therefore, MAPK signaling pathway is involved in the activation of melanocyte receptors and the complex mechanism of ligand activation up-regulates MITF by binding to the extracellular domain of the receptor [38,39], which can effectively regulate the development of melanocytes. 4.3. Identification of Candidate Genes Related to the Production of Pigment According to the analysis of sequencing data, there are 19 common differentially ex- pressed genes in three groups, most of which are related to pigment production. The ex- pression levels of 12 genes were the highest in the yellow group, but the lowest in the white group, most of which have the function of promoting the aggregation and diffusion of yellow pigment. For example, loc113059134 involves perilipoprotein, which is reported to promote pigment aggregation and is highly expressed in yellow pigment cells that me- diate the aggregation and transproduction of red/yellow pigments [22]. The gene Slc2a11b plays an important role in the growth and differentiation of yellow and white pigment Biology 2023, 12, 7 10 of 12 cells, which is consistent with previous studies on medaka [23]. GTP cyclohydrolase (Gch) is an important cofactor of nitric oxide synthase, a rate-limiting enzyme in pteridine syn- thesis, and catalyzes the de novo synthesis of tetrahydrobioptrexate (H4 bioptera) by GTP5. The expression of gch2 is the initial step of melanocyte and yellow pigment cells differentiation [1], and may play an important role in goldfish pigment aggregation. Fur- thermore, mlana and loc113079820 are essential pigment agglutination genes in the process of melanosome production and maintenance [24]. Among these 19 genes, the expression level of five genes is the highest in the brown group, and two genes are related to the synthesis and regulation of melanin. The gene loc113068772 involves myosin, which is related to melanosome transport and melanin aggregation [25,26]; tyrp1a involves in di- hydroxyindole carboxylic acid oxidase and participates in the production of melanin [8]. The expression levels of the remaining two genes were the highest in the white group, one was not annotated, and the other was related to Sec62. Consistent with the above func- tions, transcriptome sequencing showed that tyrp1 and mlana were significantly upregu- lated in the black group, while slc2a11b and gch2 were significantly down regulated. In a word, eight genes related to the production of pigment could be selected from common differentially expressed genes as candidate genes, which are worthy of future analysis. 5. Conclusions In conclusion, the current study paints thorough transcriptome profiles of different skin color groups (white, yellow and brown) in celestial goldfish, and eight candidate genes were selected as important functional genes involved in the color variation, includ- ing tyrp1a, slc2a11b, mlana, gch2, loc113060382, loc113079820, loc113068772 and loc113059134. Many differentially expressed genes are also enriched in important biolog- ical pathways related to pigment production, including MAPK signaling pathway, Wnt signaling pathway, tyrosine metabolism, melanogenesis, etc. Our findings lay the ground- work for future research on the molecular mechanisms behind skin color by identifying possible regulators that may be associated with pigment production in goldfish. Supplementary Materials: The following supporting materials can be downloaded at: https://www.mdpi.com/article/10.3390/biology12010007/s1. Table S1. Overview of the data for RNA sequencing. Table S2. Identified novel genes annotate by Pfam and GO function. Table S3. Differen- tially expressed mRNAs in skin tissues among three groups. Table S4. GO enrichment analysis for DE mRNAs in skin tissues among three groups. Table S5. KEGG enrichment analysis for DE mRNAs in skin tissues among three groups. Table S6. A total of 19 common differentially expressed genes among three groups. Table S7. Primers sequences for qPCR. Author Contributions: Conceptualization and supervision, X.Z.; data analysis, R.L. and Y.S.; writ- ing—original draft preparation, R.C.; writing—review and editing, R.L. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded as an independent project in the Beijing Academy of Agriculture and Forestry Sciences (Research and Application of Key Technologies for Beijing Characteristic Fish Breeding and Healthy Aquaculture) and National Freshwater Aquatic Germplasm Resource Bank (FGRC: 18537). Institutional Review Board Statement: All animal handling and manipulation procedures were based on the standards of the Chinese Council on Animal Care. Before sampling, all samples were anesthetized with 400 ppm of MS222. Data Availability Statement: The datasets provided in this study have been submitted to the online NCBI database (GSE220201). Conflicts of interest: The authors declare no conflicts of interest. Biology 2023, 12, 7 11 of 12 Reference 1. Kelsh, R.N. Genetics and Evolution of Pigment Patterns in Fish. Pigment Cell Res. 2010, 17, 326–336. 2. Barsh, G.S. The genetics of pigmentation: From fancy genes to complex traits. Trends Genet. 1996, 12, 299–305. 3. Mortazavi, A.; Williams, B.A.; Mccue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA- Seq. Nat. Methods 2008, 5, 621–628. 4. Jiang, B.; Wang, L.; Luo, M.; Fu, J.; Dong, Z. Transcriptome Analysis of Skin Color Variation During and after Overwintering of Malaysian Red Tilapia. Fish Physiol. Biochem. 2022, 48, 669–682. 5. Gan, W.; Chung, D.; Yu, W.; Chen, Z.; Song, S.; Ren, J. Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish. Sci. Rep. 2021, 11, 1815. 6. Luo, M.; Wang, L.; Yin, H.; Zhu, W.; Dong, Z. Integrated analysis of long non-coding RNA and mRNA expression in different colored skin of koi carp. BMC Genom. 2019, 20, 515. 7. Chen, Y.; Gong, Q.; Lai, J.; Song, M.; Long, Z. Transcriptome analysis identifies candidate genes associated with skin color variation in Triplophysa siluroides. Comp. Biochem. Physiol. Part D Genom. Proteom. 2020, 35, 100682. 8. Hearing, V.J.; Jiménez, M. Mammalian tyrosinase—The critical regulatory control point in melanocyte pigmentation. Int. J. Biochem. 1987, 19, 1141–1147. 9. Ferrini, U.; Mileo, A.M.; Hearing, V.J. Microheterogeneity of melanosome-bound tyrosinase from the harding-passey murine melanoma. Int. J. Biochem. 1987, 19, 227–234. 10. Kowichi, J.; Hiroyuki, H.; Thuraiayah, V.; Dong, L.; Jamal, D. Molecular Control of Melanogenesis in Malignant Melanoma: Functional Assessment of Tyrosinase and Lamp Gene Families by UV Exposure and Gene Co-Transfection, and Cloning of a cDNA Encoding Calnexin, A Possible Melanogenesis “Chaperone”. J. Dermatol. 1994, 21, 894–906. 11. Ginger, R.S.; Askew, S.E.; Ogborne, R.M.; Wilson, S.; Green, M.R. SLC24A5 Encodes a trans-Golgi Network Protein with Potassium-dependent Sodium-Calcium Exchange Activity That Regulates Human Epidermal Melanogenesis. J. Biol. Chem. 2008, 283, 5486. 12. Lamason, R.L.; Mohideen, M.A.P.; Mest, J.R.; Wong, A.C.; Norton, H.L.; Aros, M.C.; Jurynec, M.J.; Mao, X.; Humphreville, V.R.; Humbert, J.E.; et al. SLC24A5, a Putative Cation Exchanger, Affects Pigmentation in Zebrafish and Humans. Science 2005, 16, 13. Braasch, I.; Schartl, M.; Volff, J. Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol. Biol. 2007, 7, 74. 14. Pelletier, I.; Boyle, P.; Mcdonald, T.; Hesslinger, C.; Ziegler, I. The Pteridine Pathway in the Zebrafish, Danio Rerio: Development in Neural Crest-Derived Cells and Its Control by GTP Cyclohydrolase I. In Chemistry and Biology of Pteridines and Folates; Springer: Boston, MA, USA, 2002. 15. Wu, L.C.; Lin, Y.Y.; Yang, S.Y.; Weng, Y.T.; Tsai, Y.T. Antimelanogenic effect of c-phycocyanin through modulation of tyrosinase expression by upregulation of ERK and downregulation of p38 MAPK signaling pathways. J. Biomed. Sci. 2011, 18, 74. 16. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. 17. Liang, S.; Haitao, L.; Dechao, B.; Guoguang, Z.; Kuntao, Y.; Changhai, Z.; Yuanning, L.; Runsheng, C.; Yi, Z. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. 18. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. 19. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. 20. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general-purpose read summarization program. Bioinformatics 2013, 30, 923–930. 21. Likun, W.; Zhixing, F.; Xi, W.; Xiaowo, W.; Xuegong, Z. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics (Oxf. Engl.) 2010, 26, 136–138. 22. Granneman, J.G.; Kimler, V.A.; Zhang, H.; Ye, X.; Luo, X.; Postlethwait, J.H.; Thummel, R. Lipid droplet biology and evolution illuminated by the characterization of a novel perilipin in teleost fish. eLife 2017, 6, e21771. 23. Kimura, T.; Nagao, Y.; Hashimoto, H.; Yamamoto-Shiraishi, Y.I.; Yamamoto, S.; Yabe, T.; Takada, S.; Kinoshita, M.; Kuroiwa, A.; Naruse, K. Leucophores are similar to xanthophores in their specification and differentiation processes in medaka. Proc. Natl. Acad. Sci. USA 2014, 111, 7343–7348. 24. Teulings, H.E.; Willemsen, K.J.; Glykofridis, I.; Krebbers, G.; Komen, L.; Kemp, E.H.; Wolkerstorfer, A.; van der Veen, J.P.W.; Luiten, R.M.; Tjin, E. P. The antibody response against MART-1 differs in patients with melanoma-associated leucoderma and vitiligo. Pigment Cell Melanoma Res. 2015, 27, 1086–1096. 25. Heimsath, E.G.; Yim, Y.I.; Mustapha, M.; Hammer, J.A.; Cheney, R.E. Myosin-X knockout is semi-lethal and demonstrates that myosin-X functions in neural tube closure, pigmentation, hyaloid vasculature regression, and filopodia formation. Sci. Rep. 2017, 7, 17354. 26. Funkenstein, B.; Rebhan, Y.; Balas, V.; Pliatner, A.; Skopal, T.; Nadjar-Boger, E.; Du, S.J. New insights into the regulation of expression and biological activity of myostatin in fish: Lessons from sea bream. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2008, 151, S11–S12. 27. Wallat, G.K.; Lazur, A.M.; Chapman, F.A. Carotenoids of Different Types and Concentrations in Commercial Formulated Fish Diets Affect Color and Its Development in the Skin of the Red Oranda Variety of Goldfish. N. Am. J. Aquac. 2005, 67, 42–51. Biology 2023, 12, 7 12 of 12 28. Cunha, L.D.; Besen, K.P.; Ha, N.; Uczay, J.; Fabregat, T. Biofloc technology (BFT) improves skin pigmentation of goldfish (Carassius auratus). Aquaculture 2020, 522, 735132. 29. Kwon, B.S.; Haq, A.K.; Pomerantz, S.H.; Halaban, R. Isolation and sequence of a cDNA clone for human tyrosinase that maps at the mouse c-albino locus. Proc. Natl. Acad. Sci. USA 1987, 84, 7473–7477. 30. Giebel, L.B.; Strunk, K.M.; King, R.A.; Hanifin, J.M.; Spritz, R.A. A frequent tyrosinase gene mutation in classic, tyrosinase- negative (type IA) oculocutaneous albinism. Proc. Natl. Acad. Sci. USA 1990, 87, 3255–3258. 31. Goto, M.; Sato-Matsumura, K.C.; Sawamura, D.; Yokota, K.; Nakamura, H.; Shimizu, H. Tyrosinase gene analysis in Japanese patients with oculocutaneous albinism. J. Dermatol. Sci. 2004, 35, 215–220. 32. Jiang, Y.; Zhang, S.; Xu, J.; Feng, J.; Mahboob, S.; Al-Ghanim, K.A.; Sun, X.; Xu, P. Comparative Transcriptome Analysis Reveals the Genetic Basis of Skin Color Variation in Common Carp. PLoS ONE 2014, 9, e108200. 33. Hari, L.; Brault, V.; Kleber, M.; Lee, H.Y.; Ille, F.; Leimeroth, R.; Paratore, C.; Suter, U.; Kemler, R.; Sommer, L. Lineage-specific requirements of beta-catenin in neural crest development. J. Cell Biol. 2002, 159, 867–880. 34. Jin, E.J.; Erickson, C.A.; Takada, S.; Burrus, L.W. Wnt and BMP Signaling Govern Lineage Segregation of Melanocytes in the Avian Embryo. Dev. Biol. 2001, 233, 22–37. 35. Dorsky, R.I.; Randall, T.; Raible, D. Control of neural crest cell fate by the Wnt signalling pathway. Nature 1998, 396, 370–373. 36. Del Bino, S.; Duval, C.; Bernerd, F. Clinical and Biological Characterization of Skin Pigmentation Diversity and Its Consequences on UV Impact. Int. J. Mol. Sci. 2018, 19, 2668. 37. Kim, Y.M.; Lee, E.C.; Lim, H.M.; Seo, Y.K. Rice Bran Ash Mineral Extract Increases Pigmentation through the p-ERK Pathway in Zebrafish (Danio rerio). Int. J. Mol. Sci. 2019, 20, 2172. 38. Béjar, J.; Hong, Y.; Schartl, M. Mitf expression is sufficient to direct differentiation of medaka blastula derived stem cells to melanocytes. Development 2003, 130, 6545–6553. 39. Chae, K.J.; Subedi, L.; Jeong, M.; Park, Y.U.; Kim, Y.C. Gomisin N Inhibits Melanogenesis through Regulating the PI3K/Akt and MAPK/ERK Signaling Pathways in Melanocytes. Int. J. Mol. Sci. 2017, 18, 471. Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim re- sponsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Biology Multidisciplinary Digital Publishing Institute

Comprehensive Transcriptome Analysis of Different Skin Colors to Evaluate Genes Related to the Production of Pigment in Celestial Goldfish

Li, Rongni; Sun, Yansheng; Cui, Ran; Zhang, Xin
Biology , Volume 12 (1) – Dec 21, 2022
Download PDF
12 pages

Loading...

Page 2

Loading...

Page 3

Loading...

Page 4

Loading...

Page 5

Loading...

Page 6

Loading...

Page 7

Loading...

Page 8

Loading...

Page 9

Loading...

Page 10

Loading...

Page 11

Loading...

Page 12

 
/lp/multidisciplinary-digital-publishing-institute/comprehensive-transcriptome-analysis-of-different-skin-colors-to-9UbXGkGaU2

References (42)

Publisher
Multidisciplinary Digital Publishing Institute
Copyright
© 1996-2022 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. Terms and Conditions Privacy Policy
ISSN
2079-7737
DOI
10.3390/biology12010007
Publisher site
See Article on Publisher Site

Abstract

Article Comprehensive Transcriptome Analysis of Different Skin Colors to Evaluate Genes Related to the Production of Pigment in Celestial Goldfish Rongni Li, Yansheng Sun, Ran Cui * and Xin Zhang * Fisheries Science Institute, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China * Correspondence: ranc20220918@163.com (R.C.); 13901156204@163.com (X.Z.) Simple Summary: Skin color is the result of various pigments synthesized by dyeing groups, while goldfish have a greater variety of pigment cells and a wider range of skin colors. Up to now, several pigment cells have been reported, including melanophores, erythrocytes, xantho-phores, irido- phores, leucophores and cyanophores. Therefore, the key genes associated with sieve pigments will provide clues for the production of goldfish with a diversity of skin tones, benefiting the goldfish industry. In this study, high-throughput sequencing was performed on the back skin tissue of ce- lestial goldfish with different skin colors (white, yellow, brown), and it was found that melanin production, tyrosine metabolism pathway, Wnt signaling pathway, and MAPK signaling pathway may play an important role in pigmentation. The study provides valuable data for screening candi- date genes involved in pigment formation, and the results of this study will lay a foundation for further study on the expression characteristics and gene network analysis of pigment related genes. Abstract: Skin color is an important phenotypic feature of vertebrate fitness under natural condi- tions. Celestial goldfish, a common goldfish breed in China, mainly shows three kinds of skin colors including white, yellow and brown. However, the molecular genetic basis of this phenotype is still unclear. In this study, high-throughput sequencing was carried out on the back skin tissues of ce- lestial goldfish with different skin colors. About 58.46 Gb of original data were generated, filtered and blasted, and 74,297 mRNAs were obtained according to the reference transcriptome. A total of 4653 differentially expressed genes were screened out among the brown, yellow and white groups, Citation: Li, R.; Sun, Y.; Cui, R.; and the expression of melanogenesis related genes in brown goldfish was significantly higher than Zhang, X. Comprehensive the other two groups. There are 19 common differentially expressed genes among three groups, of Transcriptome Analysis of Different which eight genes are related to pigment production, including tyrp1a, slc2a11b, mlana, gch2, Skin Colors to Evaluate Genes loc113060382, loc113079820, loc113068772 and loc113059134. RT-qPCR verified that the expression Related to the Production of Pigment in Celestial Goldfish. patterns of randomly selected differentially expressed transcripts were highly consistent with those Biology 2023, 12, 7. https://doi.org/ obtained by RNA sequencing. GO and KEGG annotation revealed that these differentially ex- 10.3390/biology12010007 pressed genes were mostly enriched in pathways of the production of pigment, including melano- genesis, tyrosine metabolism, Wnt signaling pathway, MAPK signaling pathway etc. These results Academic Editors: Mingxing Chu, indicated that the external characteristics of goldfish are consistent with the analysis results at tran- Ran Di and Yufang Liu scriptome level. The results of this study will lay a foundation for further study on the expression Received: 3 December 2022 characteristics and gene network analysis of pigment related genes. Revised: 17 December 2022 Accepted: 19 December 2022 Keywords: celestial goldfish; skin color; transcriptome; pigment production Published: 21 December 2022 Copyright: © 2022 by the authors. Li- censee MDPI, Basel, Switzerland. 1. Introduction This article is an open access article Skin color is the result of diverse pigments synthesized by chromatophores; different distributed under the terms and con- ditions of the Creative Commons At- kinds of pigment cells are found in various vertebrates, while there are more kinds of tribution (CC BY) license (https://cre- pigment cells and a wider range of skin colors in goldfish. Up to now, several pigment ativecommons.org/licenses/by/4.0/). Biology 2023, 12, 7. https://doi.org/10.3390/biology12010007 www.mdpi.com/journal/biology Biology 2023, 12, 7 2 of 12 cells have been reported, including melanophores, erythrocytes, xanthophores, irido- phores, leucophores and cyanophores [1,2]. Although the production of skin color in fish is complex and varies depending on the breed, feed, and water environment, genes are the most basic and essential element that influences skin pigmentation. RNA-seq is a study of gene structure and function at the overall transcriptome level, based on high-throughput sequencing technology, and reveals the molecular mechanism in specific biological processes [3]. Many studies have been devoted to explaining the ge- netic basis of pigment production by RNA-seq in recent years, and a certain number of important genes and biological pathways have been found, including melanogenesis and tyrosine metabolism. Jiang et al. [4] performed transcriptome analysis of skin color during and after the overwintering of red tilapia, showing that pigmentation-related genes in- clude tyrp1, tyr, pmel, MITF and apoptosis and autophagy-related genes possibly regulate the molecular mechanisms of color change. Gan et al. [5] performed comparative tran- scriptome analysis on four skin colors of goldfish and concluded that the melanogenesis and tyrosine metabolism pathway may affect the skin pigment agglutination process. Luo et al. [6] found that biological pathway was mainly enriched in tyrosine metabolism and melanogenesis through transcriptome analysis of koi fish with white, brown and red skin colors. Chen et al. [7] reported that, consistent with goldfish, the differential expressed genes of various skin colors in Silurus meridionalis is also mainly concentrated in mel- anogenesis, melanin metabolism and melanin agglutination hormone activity pathway. In addition, tyrosinase (TYR) [8–10], solute carrier family 24 member 5 [11,12], GTP_cy- clohydroI (gch) [13,14] and microphthalmia associated transcription factor (MITF) [15] have been reported as important genes involved in melanogenesis. Celestial goldfish, as a unique variety in China, is deeply loved by people with eyes turned upward and color- ful body colors. However, the key genes and biological processes, which are involved in the production of pigment in celestial goldfish, are still unclear. The molecular process of the production of pigment in celestial goldfish has not yet undergone a comprehensive examination. Therefore, clarifying the roles of the differential expression genes can contribute to understanding the changes associated with the pro- duction of pigment in different skin color groups. The purpose of this study is to obtain an overview of transcriptome profiles in three kinds of skin colors including white, yellow and brown in skin tissue. In addition, we identified differentially expressed genes that were possibly involved in the production of pigment and screened out common differen- tial expression genes among three groups. In-depth knowledge of the molecular processes underlying the production of pigment in celestial goldfish is provided in this report. 2. Material and Method 2.1. Animal Sample Collection The celestial goldfish used in the experiment was taken from Xiaotangshan Aquacul- ture Base, Fisheries Research Institute, Beijing Academy of Agricultural and Forestry Sci- ences (Beijing, China). The experimental group used three brown and three yellow gold- fish, and the control group used three white goldfish. All three groups of goldfish came from half-sibs, were raised under the same conditions until they were 4 months old, and then were slaughtered on the same day. The fish body length is 5–6 cm, and the goldfish samples collected in the observation experiment were faded and the body color was sta- bilized. The skin tissues from the backs of nine goldfish were extracted and stored in liq- uid nitrogen or at −80 °C for later use. 2.2. RNA Extraction and Sequencing According to the instructions, skin tissue samples in three groups were taken, and the total RNA of each sample was extracted by RNA extraction kit (Invitrogen, Carlsbad, CA, USA). To analyze the integrity of the extracted total RNA and whether there was DNA contamination, 10 g/L agarose gel was used. The purity of RNA was detected by Biology 2023, 12, 7 3 of 12 nucleic acid protein determination spectrophotometer (IMPLEN, Munich, Bavaria, Ger- many). The Agilent 2100 bioanalyzer was used to accurately detect the integrity and total amount of RNA. According to the instructions, NEB-Next Ultra TM RNA Library Prep Kit (Illumina) was used to build the library. The cDNA was synthesized by magnetic bead enrichment method, and the cDNA end was repaired, amplified and purified to obtain the cDNA library of the skin tissue of celestial goldfish. 2.3. Bioinformatic Analysis The original data were filtered by removing N-containing, jointed and low-quality reads, and the GC content, Q20 and Q30 of clean data were calculated. Carassius auratus reference genome and gene model annotation files were downloaded from the NCBI web- site (https://www.ncbi.nlm.nih.gov/assembly/GCA_003368295.1 (accessed on 1 January 2022), and the parametric analysis of transcriptome sequencing data was performed. HISAT2 [16] was used to align paired-end clean reads to the reference genome. The mapped readings of each skin sample were put together using StringTie [16] with a refer- ence-based method, and then new transcripts were annotated by Pfam [17], GO [18], KEGG [19] and other databases, and the new genes were predicted. Using the Feature- Counts program, the readings mapped to each gene were calculated [20], and the FPKM of each gene (the expected number of every thousand base segments of transcript se- quence segments sequenced per million base pairs) was calculated based on the gene length. Differential expression analysis was performed using DESeq [21], with a default significant difference of two fold changes and a q value of 0.05. 2.4. Gene ontology and kyoto encyclopedia of genes and genomes enrichment analysis Cluster-Profile software was used to analyze the enrichment of differentially ex- pressed genes by GO function and KEGG pathway, and the significant enrichment thresh- old was p value < 0.05. The main biochemical metabolic pathways and signal transduction pathways involved in color change were screened out. 2.5. Validation of Differentially Expressed Transcripts RT-qPCR verified the expression pattern of 10 differentially expressed transcripts randomly selected from the skin tissue of the back. Total RNA of the sample was reverse transcribed into cDNA by prime script TM 1st stand cDNA synthesis kit. β-actin was se- lected as endogenous control gene (Table S7). Three biological replicates and triple reac- tions for each sample were used for all qPCR validations. The relative expression of each gene was calculated by 2-△△Ct method, and the log2 value of multiple expression be- tween the two groups was calculated. The column chart was drawn together with the RNA-seq sequencing results, and the trend was compared. 2.6. Analytical Statistics The means and standard error of the means (SEM) are used to express data. In order to control the error rate, Hochberg and Benjamini methods are used to adjust the p value; differences were deemed statistically significant at p < 0.05. 3. Results 3.1. Blast Analysis of Transcriptome Sequencing The transcriptome profiling of the skin tissue from goldfish of different colors were obtained to evaluate the genes involved in the production of melanin. A total of 58.46 Gb of filtered base was obtained, and the base amount of each sample reached more than 5.96 Gb, and 87.72~92.96% of these readings were specifically compared with the reference ge- nome (https://www.ncbi.nlm.nih.gov/assembly/GCA_003368295.1, accessed on 1 January 2022). The base error rate of each sample sequencing is less than 0. 03%, and Q30 is greater than 93.90%, which indicates that the sequencing data quality is reliable (Table S1). Biology 2023, 12, 7 4 of 12 Celestial goldfish is a non-model species, and its genetic annotation in production is usually not perfect. After the alignment was completed, the annotated transcription re- gions of the genome were searched, and the aligned sequences were assembled, so as to mine new transcripts or new genes in the species. In this study, new transcripts were as- sembled and annotated in Pfam, GO and KEGG databases and 2205 new genes were pre- dicted and annotated (Table S2). 3.2. Differentially Expressed Genes in Goldfish with Different Skin Colors Through the comparison of different skin color groups, the differentially expressed genes are shown (Figure 1a,b). When compared between brown and yellow groups, a total of 1000 differentially expressed genes were screened out, of which 651 were up-regulated and 349 were down-regulated. When comparing between white and brown groups, we screened out 3033 differentially expressed genes, of which 1196 genes were up-regulated and 1837 genes were down-regulated. Consistent with the visual contrast between the brown and white groups, the number of differentially expressed genes in this group is more than other groups through RNA-seq. When comparing between white and yellow groups, a total of 1942 differentially expressed genes were screened out, of which 761 were up-regulated and 1181 were down-regulated. Some DE mRNAs were specifically ex- pressed in the brown group, including: loc113077208, loc113119661, loc113085971, loc113066815; Some DE mRNA was specifically expressed in the white group, including: clic5a, loc113104870, loc113118795; There are also some specific expression of DE mRNAs in the yellow group, including: slc2a11b, asz1, loc113092873. These genes may regulate the production of pigment in goldfish (Table S3). The differentially expressed genes of three groups of goldfish were clustered and drawn into a heat map (Figure 1c). It can be seen that the expression patterns of the brown group and the yellow group are more similar, while the white group and the other two groups are at a far distance. Biology 2023, 12, 7 5 of 12 Figure 1. (a). The different skin colors in the celestial goldfish. Left: white (W); center：yellow (Y); right: brown (B). (b) The bar graph of number of DEGs among three different skin colors. (c) Clus- tering heat map of DEGs in the three skin color groups. (d) DEGs number and Venn diagram of the overlap of the different skin color groups. W: white; Y: yellow; B: brown. 3.3. Functional Analysis of Differentially Expressed Genes GO enrichment analysis was conducted between brown, yellow and white group ac- cording to three categories: biological process, cell components and molecular functions. The most significant 30 entries were selected from the GO enrichment analysis results to draw scatter plots (Figure 2, Table S4). GO enrichment entries in the white and yellow groups are significant in the biological process (BP), including peptide cross-linking (GO:0018149), small molecule biosynthetic process (GO:0044283), protein tyrosine phos- phatase activity (GO:0004725), etc. The GO entries of the brown and yellow groups, white and brown groups were both significantly enriched in the cellular component (CC), mainly including melanosome transport-related: troponin complex (GO:0005861), myofil- ament (GO:0036379), contractile fiber (GO:0043292), etc. GO enrichment concluded that differentially expressed genes were mainly enriched in the pathways related to tyrsine metabolism and melanosome transport. Biology 2023, 12, 7 6 of 12 Figure 2. (a). The scatter figure of GO enrichment in white vs. yellow. (b) The scatter figure of GO enrichment of DEGs in brown vs. yellow. (c) The scatter figure of GO enrichment of DEGs in white vs. brown. Count: number of genes in each pathway. GeneRatio: ratio of the number of target genes divided by the total number of genes in each pathway. The differentially expressed genes of three groups in goldfish were compared to KEGG database for pathway enrichment analysis, and the most significant 20 KEGG path- ways were selected to draw columns (Figure 3, Table S5). According to the results of KEGG enrichment, KEGG pathways among three groups are all enriched to pigment pro- duction related, including Wnt signaling pathway (dre04310), tyrosine metabolism (dre00350), melanogenesis (dre04916) and MAPK signaling pathway (dre04010). These re- sults reflect that the difference of skin color of goldfish may be the result of the regulation of multiple signal pathways, and the different skin color characteristics in this study are closely related to pigment production. Biology 2023, 12, 7 7 of 12 Figure 3. (a). The bar graph of KEGG pathway enrichment of DEGs in white vs. yellow. (b) The bar graph of KEGG pathway enrichment of DEGs in brown vs. yellow. (c) The bar graph of KEGG path- way enrichment of DEGs in white vs. brown. 3.4. Common Expression Gene Analysis According to the overlapping of differential genes among different skin color groups, the results showed that 19 differential genes belonged to the common differential genes among three groups (Figure 1d, Tables 1 and S6), among which eight genes were reported to be related to pigment production and regulation [8,14,22–26], including: tyrp1a, slc2a11b, mlana, gch2, loc113060382, loc113079820, loc113068772 and loc113059134. The re- maining genes mainly regulate cell proliferation, differentiation and apoptosis, including transcription elongation factor a protein 3, translation initiation factor if-2, RAS and EF- hand domain containing, translocation protein Sec62, et al. Table 1. A total of 19 common differentially expressed genes among three groups. FPKM Gene ID Description White Yellow Brown 113116617 0.053316 14.85341 4.939858 uncharacterized 113059134 0.101642 11.21777 3.472121 perilipin-3-like novel.91 0.024732 41.3923 11.47289 -- 113094292 0 22.96619 9.53373 putative defense protein 3 113054403 0 7.702924 1.333329 diacylglycerol O-acyltransferase 2 113060382 0.021323 4.63168 1.063391 retinol dehydrogenase 7-like solute carrier family facilitated glucose 113072589 0 3.77947 0.555086 transporter member 11-like melanoma antigen recognized by T-cells 1- 113081155 0.065602 12.08961 3.294872 like novel.919 0.052276 2.696407 9.566884 -- 113056131 0.017281 4.126237 0.910635 GTP cyclohydrolase 1-like melanoma antigen recognized by T-cells 1- 113079820 0.023412 8.795228 2.98262 like 113116699 18.29243 526.6035 80.16275 translation initiation factor IF-2-like Biology 2023, 12, 7 8 of 12 113044679 0.872271 0.010594 0.200235 uncharacterized 113107421 0.154702 2.39033 0.707736 RAS and EF-hand domain containing transcription elongation factor A protein 3- 113059190 1.056176 0.085419 5.333497 like 113068772 0.148587 0.020619 1.319104 myosin-6-like 113119606 3.652859 1.230851 0.193719 translocation protein SEC62-like C6-dihydroxyindole-2-carboxylic acid oxi- 113106535 0.017673 0.195103 1.180688 dase-like novel.58 0.154461 0.976559 4.569909 -- 3.5. Validation of Differentially Expressed Genes Expression Eight genes were randomly selected from the differentially expressed genes for qPCR verification, and the results showed that the gene expression trend obtained by qPCR was basically consistent with the results of RNA-seq (Figure 4). It indicates that the sequencing results of transcriptome in this study are reliable, and the differentially expressed genes among different groups can be accurately detected. Figure 4. The expression levels of mRNA verified by qPCR. LOC113045861-LOC113043321: left: yel- low group, right: brown group; LOC113075300: left: white group, right: brown group; LOC113040942-LOC113066051: left: yellow group, right: white group. * indicates p < 0.05. 4. Discussion 4.1. Composition Analysis of Pigment Cells in Different Skin Colors Skin color plays an important role in animal courtship, escape from enemies and sig- nal communication, and is one of the most diverse forms of kinvertebrates. Due to diverse color pattern, goldfish could be a perfect model to study the genetic mechanism of pig- mentation. The research on the formation mechanism of fish skin color provides an im- portant theoretical reference for the study of physiology, pathology and genetics of hu- man and other biological pigment cells. Although the production of skin color in fish is complex and varies depending on the breed, feed, and water environment [27,28], genes are the most basic and essential element that influences skin pigmentation. In this study, GO enrichment reflected the different components of pigment cells in different skin col- ors. The skin cells of white goldfish are mainly composed of white pigment cell, while those of yellow goldfish are mainly yellow pigment cells. Studies have reported that that characteristics and differentiation processes of white pigment and yellow pigment cells Biology 2023, 12, 7 9 of 12 are similar in medaka [23], which also explains why there is less difference in cellular composition between white and yellow groups than in comparison with the brown group. Inconsistent with white and yellow groups, there were dendritic melanocytes in the skin of brown groups, and GO enrichment enriched the related pathways with protein tyrosine phosphatase activity. Tyrosinase (TYR) is the rate limiting enzyme for melanin synthesis, and its gene mutation will cause albinism in animals. Since the cloning of human tyrosi- nase gene in 1987 [29], much research has been carried out on the relationship between tyrosine gene and skin color (mainly albinism) in humans and animals [30,31]. Melanin is a polymer synthesized by tyrosine hydroxylation by tyrosinase catalyzed in vivo, which involves three enzymes in the tyrosinase gene family, tyrosinase (TYR), tyrosinase-related protein 1 (TRP1) and tyrosinase-related protein-2 (TRP2). After melanin is produced, it travels from the dendritic tip of melanocytes to keratinocytes. Furthermore, GO enrich- ment also enriched the related pathways with the structure of iridescent cells, actin and enzyme activity. It has been reported that the directional movement of pigment particles is closely related to actin. 4.2. Regulated Genes in Brown Skin Indicated the Molecular Mechanism of the Melanogenesis Pathway At present, most studies on the skin color of goldfish are focused on pigment cell metabolism, pigment deposition and production, especially melanin production and transport, which is consistent with the results in this study. For example, Jiang et al. [32] performed comparative transcriptome analysis on the skin colors of Xingguo Red Carp and Yellow River Carp and also concluded that melanogenesis pathway affected the skin pigment agglutination process in common carp. According to our research, the differen- tially expressed genes were not only enriched in the melanogenesis pathways mentioned above in the existing studies, but MAPK signaling pathway and Wnt signaling pathway were also enriched. Wnt signaling pathways, particularly Wnt1, Wnt3a and Wnt5a, have been reported to play an important role in the development of melanocytes [33–36]. Wnt3a has been reported to reduce the number of neuronal cells and glial cells in quail, while significantly increasing the number of melanocytes [34]. Another report on zebrafish sug- gests that after injection of domain-inactivated Wnt1, the Wnt signaling pathway is blocked, the number of pigment cells is reduced, and neural crest cells are easily differen- tiated into neuronal cells and glial cells [35]. Vertebrate melanocytes almost always origi- nate from the neural crest, and the cells derived from the neural crest undergo obvious differentiation until they are dispersed, and produce peripheral neurons, glial cells and melanocytes at precise locations [36]. Furthermore, phosphorylation of ERK in MAPK sig- naling pathway inhibits the production of melanin in zebrafish, while inhibition of phos- phorylation of ERK can promote the production of melanin in zebrafish by increasing the expression of MITF genes, thereby upregulating downstream genes TYR, TRP1 and TRP2 [37]. The expression of MITF was also able to differentiate pluripotent stem cells into mel- anocytes in medaka [38]. Therefore, MAPK signaling pathway is involved in the activation of melanocyte receptors and the complex mechanism of ligand activation up-regulates MITF by binding to the extracellular domain of the receptor [38,39], which can effectively regulate the development of melanocytes. 4.3. Identification of Candidate Genes Related to the Production of Pigment According to the analysis of sequencing data, there are 19 common differentially ex- pressed genes in three groups, most of which are related to pigment production. The ex- pression levels of 12 genes were the highest in the yellow group, but the lowest in the white group, most of which have the function of promoting the aggregation and diffusion of yellow pigment. For example, loc113059134 involves perilipoprotein, which is reported to promote pigment aggregation and is highly expressed in yellow pigment cells that me- diate the aggregation and transproduction of red/yellow pigments [22]. The gene Slc2a11b plays an important role in the growth and differentiation of yellow and white pigment Biology 2023, 12, 7 10 of 12 cells, which is consistent with previous studies on medaka [23]. GTP cyclohydrolase (Gch) is an important cofactor of nitric oxide synthase, a rate-limiting enzyme in pteridine syn- thesis, and catalyzes the de novo synthesis of tetrahydrobioptrexate (H4 bioptera) by GTP5. The expression of gch2 is the initial step of melanocyte and yellow pigment cells differentiation [1], and may play an important role in goldfish pigment aggregation. Fur- thermore, mlana and loc113079820 are essential pigment agglutination genes in the process of melanosome production and maintenance [24]. Among these 19 genes, the expression level of five genes is the highest in the brown group, and two genes are related to the synthesis and regulation of melanin. The gene loc113068772 involves myosin, which is related to melanosome transport and melanin aggregation [25,26]; tyrp1a involves in di- hydroxyindole carboxylic acid oxidase and participates in the production of melanin [8]. The expression levels of the remaining two genes were the highest in the white group, one was not annotated, and the other was related to Sec62. Consistent with the above func- tions, transcriptome sequencing showed that tyrp1 and mlana were significantly upregu- lated in the black group, while slc2a11b and gch2 were significantly down regulated. In a word, eight genes related to the production of pigment could be selected from common differentially expressed genes as candidate genes, which are worthy of future analysis. 5. Conclusions In conclusion, the current study paints thorough transcriptome profiles of different skin color groups (white, yellow and brown) in celestial goldfish, and eight candidate genes were selected as important functional genes involved in the color variation, includ- ing tyrp1a, slc2a11b, mlana, gch2, loc113060382, loc113079820, loc113068772 and loc113059134. Many differentially expressed genes are also enriched in important biolog- ical pathways related to pigment production, including MAPK signaling pathway, Wnt signaling pathway, tyrosine metabolism, melanogenesis, etc. Our findings lay the ground- work for future research on the molecular mechanisms behind skin color by identifying possible regulators that may be associated with pigment production in goldfish. Supplementary Materials: The following supporting materials can be downloaded at: https://www.mdpi.com/article/10.3390/biology12010007/s1. Table S1. Overview of the data for RNA sequencing. Table S2. Identified novel genes annotate by Pfam and GO function. Table S3. Differen- tially expressed mRNAs in skin tissues among three groups. Table S4. GO enrichment analysis for DE mRNAs in skin tissues among three groups. Table S5. KEGG enrichment analysis for DE mRNAs in skin tissues among three groups. Table S6. A total of 19 common differentially expressed genes among three groups. Table S7. Primers sequences for qPCR. Author Contributions: Conceptualization and supervision, X.Z.; data analysis, R.L. and Y.S.; writ- ing—original draft preparation, R.C.; writing—review and editing, R.L. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded as an independent project in the Beijing Academy of Agriculture and Forestry Sciences (Research and Application of Key Technologies for Beijing Characteristic Fish Breeding and Healthy Aquaculture) and National Freshwater Aquatic Germplasm Resource Bank (FGRC: 18537). Institutional Review Board Statement: All animal handling and manipulation procedures were based on the standards of the Chinese Council on Animal Care. Before sampling, all samples were anesthetized with 400 ppm of MS222. Data Availability Statement: The datasets provided in this study have been submitted to the online NCBI database (GSE220201). Conflicts of interest: The authors declare no conflicts of interest. Biology 2023, 12, 7 11 of 12 Reference 1. Kelsh, R.N. Genetics and Evolution of Pigment Patterns in Fish. Pigment Cell Res. 2010, 17, 326–336. 2. Barsh, G.S. The genetics of pigmentation: From fancy genes to complex traits. Trends Genet. 1996, 12, 299–305. 3. Mortazavi, A.; Williams, B.A.; Mccue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA- Seq. Nat. Methods 2008, 5, 621–628. 4. Jiang, B.; Wang, L.; Luo, M.; Fu, J.; Dong, Z. Transcriptome Analysis of Skin Color Variation During and after Overwintering of Malaysian Red Tilapia. Fish Physiol. Biochem. 2022, 48, 669–682. 5. Gan, W.; Chung, D.; Yu, W.; Chen, Z.; Song, S.; Ren, J. Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish. Sci. Rep. 2021, 11, 1815. 6. Luo, M.; Wang, L.; Yin, H.; Zhu, W.; Dong, Z. Integrated analysis of long non-coding RNA and mRNA expression in different colored skin of koi carp. BMC Genom. 2019, 20, 515. 7. Chen, Y.; Gong, Q.; Lai, J.; Song, M.; Long, Z. Transcriptome analysis identifies candidate genes associated with skin color variation in Triplophysa siluroides. Comp. Biochem. Physiol. Part D Genom. Proteom. 2020, 35, 100682. 8. Hearing, V.J.; Jiménez, M. Mammalian tyrosinase—The critical regulatory control point in melanocyte pigmentation. Int. J. Biochem. 1987, 19, 1141–1147. 9. Ferrini, U.; Mileo, A.M.; Hearing, V.J. Microheterogeneity of melanosome-bound tyrosinase from the harding-passey murine melanoma. Int. J. Biochem. 1987, 19, 227–234. 10. Kowichi, J.; Hiroyuki, H.; Thuraiayah, V.; Dong, L.; Jamal, D. Molecular Control of Melanogenesis in Malignant Melanoma: Functional Assessment of Tyrosinase and Lamp Gene Families by UV Exposure and Gene Co-Transfection, and Cloning of a cDNA Encoding Calnexin, A Possible Melanogenesis “Chaperone”. J. Dermatol. 1994, 21, 894–906. 11. Ginger, R.S.; Askew, S.E.; Ogborne, R.M.; Wilson, S.; Green, M.R. SLC24A5 Encodes a trans-Golgi Network Protein with Potassium-dependent Sodium-Calcium Exchange Activity That Regulates Human Epidermal Melanogenesis. J. Biol. Chem. 2008, 283, 5486. 12. Lamason, R.L.; Mohideen, M.A.P.; Mest, J.R.; Wong, A.C.; Norton, H.L.; Aros, M.C.; Jurynec, M.J.; Mao, X.; Humphreville, V.R.; Humbert, J.E.; et al. SLC24A5, a Putative Cation Exchanger, Affects Pigmentation in Zebrafish and Humans. Science 2005, 16, 13. Braasch, I.; Schartl, M.; Volff, J. Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol. Biol. 2007, 7, 74. 14. Pelletier, I.; Boyle, P.; Mcdonald, T.; Hesslinger, C.; Ziegler, I. The Pteridine Pathway in the Zebrafish, Danio Rerio: Development in Neural Crest-Derived Cells and Its Control by GTP Cyclohydrolase I. In Chemistry and Biology of Pteridines and Folates; Springer: Boston, MA, USA, 2002. 15. Wu, L.C.; Lin, Y.Y.; Yang, S.Y.; Weng, Y.T.; Tsai, Y.T. Antimelanogenic effect of c-phycocyanin through modulation of tyrosinase expression by upregulation of ERK and downregulation of p38 MAPK signaling pathways. J. Biomed. Sci. 2011, 18, 74. 16. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. 17. Liang, S.; Haitao, L.; Dechao, B.; Guoguang, Z.; Kuntao, Y.; Changhai, Z.; Yuanning, L.; Runsheng, C.; Yi, Z. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013, 41, e166. 18. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. 19. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. 20. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general-purpose read summarization program. Bioinformatics 2013, 30, 923–930. 21. Likun, W.; Zhixing, F.; Xi, W.; Xiaowo, W.; Xuegong, Z. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics (Oxf. Engl.) 2010, 26, 136–138. 22. Granneman, J.G.; Kimler, V.A.; Zhang, H.; Ye, X.; Luo, X.; Postlethwait, J.H.; Thummel, R. Lipid droplet biology and evolution illuminated by the characterization of a novel perilipin in teleost fish. eLife 2017, 6, e21771. 23. Kimura, T.; Nagao, Y.; Hashimoto, H.; Yamamoto-Shiraishi, Y.I.; Yamamoto, S.; Yabe, T.; Takada, S.; Kinoshita, M.; Kuroiwa, A.; Naruse, K. Leucophores are similar to xanthophores in their specification and differentiation processes in medaka. Proc. Natl. Acad. Sci. USA 2014, 111, 7343–7348. 24. Teulings, H.E.; Willemsen, K.J.; Glykofridis, I.; Krebbers, G.; Komen, L.; Kemp, E.H.; Wolkerstorfer, A.; van der Veen, J.P.W.; Luiten, R.M.; Tjin, E. P. The antibody response against MART-1 differs in patients with melanoma-associated leucoderma and vitiligo. Pigment Cell Melanoma Res. 2015, 27, 1086–1096. 25. Heimsath, E.G.; Yim, Y.I.; Mustapha, M.; Hammer, J.A.; Cheney, R.E. Myosin-X knockout is semi-lethal and demonstrates that myosin-X functions in neural tube closure, pigmentation, hyaloid vasculature regression, and filopodia formation. Sci. Rep. 2017, 7, 17354. 26. Funkenstein, B.; Rebhan, Y.; Balas, V.; Pliatner, A.; Skopal, T.; Nadjar-Boger, E.; Du, S.J. New insights into the regulation of expression and biological activity of myostatin in fish: Lessons from sea bream. Comp. Biochem. Physiol. Part A Mol. Integr. Physiol. 2008, 151, S11–S12. 27. Wallat, G.K.; Lazur, A.M.; Chapman, F.A. Carotenoids of Different Types and Concentrations in Commercial Formulated Fish Diets Affect Color and Its Development in the Skin of the Red Oranda Variety of Goldfish. N. Am. J. Aquac. 2005, 67, 42–51. Biology 2023, 12, 7 12 of 12 28. Cunha, L.D.; Besen, K.P.; Ha, N.; Uczay, J.; Fabregat, T. Biofloc technology (BFT) improves skin pigmentation of goldfish (Carassius auratus). Aquaculture 2020, 522, 735132. 29. Kwon, B.S.; Haq, A.K.; Pomerantz, S.H.; Halaban, R. Isolation and sequence of a cDNA clone for human tyrosinase that maps at the mouse c-albino locus. Proc. Natl. Acad. Sci. USA 1987, 84, 7473–7477. 30. Giebel, L.B.; Strunk, K.M.; King, R.A.; Hanifin, J.M.; Spritz, R.A. A frequent tyrosinase gene mutation in classic, tyrosinase- negative (type IA) oculocutaneous albinism. Proc. Natl. Acad. Sci. USA 1990, 87, 3255–3258. 31. Goto, M.; Sato-Matsumura, K.C.; Sawamura, D.; Yokota, K.; Nakamura, H.; Shimizu, H. Tyrosinase gene analysis in Japanese patients with oculocutaneous albinism. J. Dermatol. Sci. 2004, 35, 215–220. 32. Jiang, Y.; Zhang, S.; Xu, J.; Feng, J.; Mahboob, S.; Al-Ghanim, K.A.; Sun, X.; Xu, P. Comparative Transcriptome Analysis Reveals the Genetic Basis of Skin Color Variation in Common Carp. PLoS ONE 2014, 9, e108200. 33. Hari, L.; Brault, V.; Kleber, M.; Lee, H.Y.; Ille, F.; Leimeroth, R.; Paratore, C.; Suter, U.; Kemler, R.; Sommer, L. Lineage-specific requirements of beta-catenin in neural crest development. J. Cell Biol. 2002, 159, 867–880. 34. Jin, E.J.; Erickson, C.A.; Takada, S.; Burrus, L.W. Wnt and BMP Signaling Govern Lineage Segregation of Melanocytes in the Avian Embryo. Dev. Biol. 2001, 233, 22–37. 35. Dorsky, R.I.; Randall, T.; Raible, D. Control of neural crest cell fate by the Wnt signalling pathway. Nature 1998, 396, 370–373. 36. Del Bino, S.; Duval, C.; Bernerd, F. Clinical and Biological Characterization of Skin Pigmentation Diversity and Its Consequences on UV Impact. Int. J. Mol. Sci. 2018, 19, 2668. 37. Kim, Y.M.; Lee, E.C.; Lim, H.M.; Seo, Y.K. Rice Bran Ash Mineral Extract Increases Pigmentation through the p-ERK Pathway in Zebrafish (Danio rerio). Int. J. Mol. Sci. 2019, 20, 2172. 38. Béjar, J.; Hong, Y.; Schartl, M. Mitf expression is sufficient to direct differentiation of medaka blastula derived stem cells to melanocytes. Development 2003, 130, 6545–6553. 39. Chae, K.J.; Subedi, L.; Jeong, M.; Park, Y.U.; Kim, Y.C. Gomisin N Inhibits Melanogenesis through Regulating the PI3K/Akt and MAPK/ERK Signaling Pathways in Melanocytes. Int. J. Mol. Sci. 2017, 18, 471. Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim re- sponsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Journal

Biology – Multidisciplinary Digital Publishing Institute

Published: Dec 21, 2022

Keywords: celestial goldfish; skin color; transcriptome; pigment production

There are no references for this article.