Gene Expression Variations of Red-White Skin Coloration in Common Carp (Cyprinus carpio)

Barramundi (Lates calcarifer) rare coloration patterns: a multiomics approach to understand the "panda" phenotype

The barramundi (Lates calcarifer), a significant aquaculture species, typically displays silver to bronze coloration. However, attention is now drawn to rare variants like the "panda" phenotype, characterized by blotch-like patterns of black (PB) and golden (PG) patches. This phenotype presents an opportunity to explore the molecular mechanisms underlying color variations in teleosts. Unlike stable color patterns in many fish, the "panda" variant demonstrates phenotypic plasticity, responding dynamically to unknown cues. We propose a complex interplay of genetic factors and epigenetic modifications, focusing on DNA methylation. Through a multiomics approach, we analyze transcriptomic and methylation patterns between PB and PG patches. Our study reveals differential gene expression related to melanosome trafficking and chromatophore differentiation. Although the specific gene responsible for the PB-PG difference remains elusive, candidate genes like asip1, asip2, mlph, and mreg have been identified. Methylation emerges as a potential contributor to the "panda" phenotype, with changes in gene promoters like hand2 and dynamin possibly influencing coloration. This research lays the groundwork for further exploration into rare barramundi color patterns, enhancing our understanding of color diversity in teleosts. Additionally, it underscores the "panda" phenotype's potential as a model for studying adult skin coloration.

Analysis of the Skin and Brain Transcriptome of Normally Pigmented and Pseudo-Albino Southern Flounder (Paralichthys lethostigma) Juveniles to Study the Molecular Mechanisms of Hypopigmentation and Its Implications for Species Survival in the Natural Environment

Southern flounder skin pigmentation is a critical phenotypic characteristic for this species' survival in the natural environment. Normal pigmentation allows rapid changes of color for concealment to capture prey and UV light protection. In contrast, highly visible hypopigmented pseudo-albinos exhibit a compromised immune system and are vulnerable to predation, sensitive to UV exposure, and likely have poor survival in the wild. Skin and brain tissue samples from normally pigmented and hypopigmented individuals were analyzed with next-generation RNA sequencing. A total of 1,589,613 transcripts were used to identify 952,825 genes to assemble a de novo transcriptome, with 99.43% of genes mapped to the assembly. Differential gene expression and gene enrichment analysis of contrasting tissues and phenotypes revealed that pseudo-albino individuals appeared more susceptible to environmental stress, UV light exposure, hypoxia, and osmotic stress. The pseudo-albinos' restricted immune response showed upregulated genes linked to cancer development, signaling and response, skin tissue formation, regeneration, and healing. The data indicate that a modified skin collagen structure likely affects melanocyte differentiation and distribution, generating the pseudo-albino phenotype. In addition, the comparison of the brain transcriptome revealed changes in myelination and melanocyte stem cell activity, which may indicate modified brain function, reduced melanocyte migration, and impaired vision.

Genetic mapping and molecular mechanism behind color variation in the Asian vine snake

BACKGROUND

Reptiles exhibit a wide variety of skin colors, which serve essential roles in survival and reproduction. However, the molecular basis of these conspicuous colors remains unresolved.

RESULTS

We investigate color morph-enriched Asian vine snakes (Ahaetulla prasina), to explore the mechanism underpinning color variations. Transmission electron microscopy imaging and metabolomics analysis indicates that chromatophore morphology (mainly iridophores) is the main basis for differences in skin color. Additionally, we assemble a 1.77-Gb high-quality chromosome-anchored genome of the snake. Genome-wide association study and RNA sequencing reveal a conservative amino acid substitution (p.P20S) in SMARCE1, which may be involved in the regulation of chromatophore development initiated from neural crest cells. SMARCE1 knockdown in zebrafish and immunofluorescence verify the interactions among SMARCE1, iridophores, and tfec, which may determine color variations in the Asian vine snake.

CONCLUSIONS

This study reveals the genetic associations of color variation in Asian vine snakes, providing insights and important resources for a deeper understanding of the molecular and genetic mechanisms related to reptilian coloration.

Comparative transcriptome analysis of skin color-associated genes in leopard coral grouper (Plectropomus leopardus)

BACKGROUND

The leopard coral grouper (Plectropomus leopardus) is an important economic species in East Asia-Pacific countries. To meet the market demand, leopard coral grouper is facing overfishing and their population is rapidly declining. With the improvement of the artificial propagation technique, the leopard coral grouper has been successfully cultured by Fisheries Research Institute in Taiwan. However, the skin color of farmed individuals is often lacking bright redness. As such, the market price of farmed individuals is lower than wild-type.

RESULTS

To understand the genetic mechanisms of skin coloration in leopard coral grouper, we compared leopard coral grouper with different skin colors through transcriptome analysis. Six cDNA libraries generated from wild-caught leopard coral grouper with different skin colors were characterized by using the Illumina platform. Reference-guided de novo transcriptome data of leopard coral grouper obtained 24,700 transcripts, and 1,089 differentially expressed genes (DEGs) were found between red and brown skin color individuals. The results showed that nine candidate DEGs (epha2, sema6d, acsl4, slc7a5, hipk1, nol6, timp2, slc25a42, and kdf1) significantly associated with skin color were detected by using comparative transcriptome analysis and quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The findings may provide genetic information for further skin color research, and to boost the market price of farmed leopard coral grouper by selective breeding.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

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

DNA Methylation Analyses Unveil a Regulatory Landscape in the Formation of Nacre Color in Pearl Oyster Pinctada fucata martensii

Pearl color is regulated by genetics, biological pigments, and organic matrices and an important factor that influences the pearl economic value. The epigenetic regulation mechanism underlying pearl pigmentation remains poorly understood. In this study, we collected the mantle pallial (MP) and mantle central (MC) of the golden-lipped strain, and MP of the silver-lipped strain of pearl oyster Pinctada fucata martensii. The whole-genome bisulfite sequencing (WGBS) technology was employed to investigate the possible implication of epigenetic factors regulating nacre color variation. Our results revealed approximately 2.5% of the cytosines in the genome of the P. fucata martensii were methylated, with the CG methylation type was in most abundance. Overall, we identified 12, 621 differentially methylated regions (DMRs) corresponding to 3,471 DMR-associated genes (DMGs) between the two comparison groups. These DMGs were principally enriched into KEGG metabolic pathways including ABC transporters, Terpenoid backbone biosynthesis, and fatty acid degradation. In addition, integrating information about DMGs, DEGs, and function annotation indicated eight genes LDLR, NinaB, RDH, CYP, FADS, fn3, PU-1, KRMP as the candidate genes related to pigmentation of nacre color. A further study proved that the pigment in nacre is violaxanthin. The results of our study provide the support that there is an association between nacre color formation and DNA methylation profiles and will help to reveal the epigenetic regulation of nacre pigmentation formation in pearl oyster P. fucata martensii.

Transcriptome analysis of skin color variation during and after overwintering of Malaysian red tilapia

The commercial value of red tilapia is hampered by variations in skin color during overwintering. In this study, three types of skin of red tilapia, including the skin remained pink color during and after overwintering (P), the skin changed from pink color to black color during overwintering and remained black color after overwintering (P-B), and the skin changed from pink color to black color during overwintering but recovered to pink color when the temperature rose after overwintering (P-B-P), were used to analyze their molecular mechanisms of color variation. The transcriptome results revealed that the P, P-B, and P-B-P libraries had 43, 42, and 43 million clean reads, respectively. The top 10 abundance mRNAs and specific mRNAs (specificity measure SPM > 0.9) were screened. After comparing intergroup gene expression levels, there were 2528, 1924, and 1939 differentially expressed genes (DEGs) between P-B-P and P-B, P-B-P and P, and P-B and P, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of color-related mRNAs showed that a number of DEGs, including tyrp1, tyr, pmel, mitf, mc1r, asip, tat, hpdb, and foxd3, might play a potential role in pigmentation. Additionally, the co-expression patterns of genes were detected within the pigment-related pathways by the PPI network from P-B vs. P group. Furthermore, DEGs from the apoptosis and autophagy pathways, such as baxα, beclin1, and atg7, might be involved in the fading of red tilapia melanocytes. The findings will aid in understanding the molecular mechanism underlying skin color variation in red tilapia during and after overwintering as well as lay a foundation for future research aimed at improving red tilapia skin color characteristics.

Analysis of yellow mutant rainbow trout transcriptomes at different developmental stages reveals dynamic regulation of skin pigmentation genes

Yellow mutant rainbow trout (YR), an economically important aquaculture species, is popular among consumers due to its excellent meat quality and attractive appearance. Skin color is a key economic trait for YR, but little is known about the molecular mechanism of skin color development. In this study, YR skin transcriptomes were analyzed to explore temporal expression patterns of pigmentation-related genes in three different stages of skin color development. In total, 16,590, 16,682, and 5619 genes were differentially expressed between fish at 1 day post-hatching (YR1d) and YR45d, YR1d and YR90d, and YR45d and YR90d. Numerous differentially expressed genes (DEGs) associated with pigmentation were identified, and almost all of them involved in pteridine and carotenoid synthesis were significantly upregulated in YR45d and YR90d compared to YR1d, including GCH1, PTS, QDPR, CSFIR1, SLC2A11, SCARB1, DGAT2, PNPLA2, APOD, and BCO2. Interestingly, many DEGs enriched in melanin synthesis pathways were also significantly upregulated, including melanogenesis (MITF, MC1R, SLC45A2, OCA2, and GPR143), tyrosine metabolism (TYR, TYRP1, and DCT), and MAPK signaling (KITA) pathways. Using short time-series expression miner, we identified eight differential gene expression pattern profiles, and DEGs in profile 7 were associated with skin pigmentation. Protein–protein interaction network analysis showed that two modules were related to xanthophores and melanophores. In addition, 1,812,329 simple sequence repeats and 2,011,334 single-nucleotide polymorphisms were discovered. The results enhance our understanding of the molecular mechanism underlying skin pigmentation in YR, and could accelerate the molecular breeding of fish species with valuable skin color traits and will likely be highly informative for developing new therapeutic approaches to treat pigmentation disorders and melanoma.

Pterin-based pigmentation in animals

Pterins are one of the major sources of bright coloration in animals. They are produced endogenously, participate in vital physiological processes and serve a variety of signalling functions. Despite their ubiquity in nature, pterin-based pigmentation has received little attention when compared to other major pigment classes. Here, we summarize major aspects relating to pterin pigmentation in animals, from its long history of research to recent genomic studies on the molecular mechanisms underlying its evolution. We argue that pterins have intermediate characteristics (endogenously produced, typically bright) between two well-studied pigment types, melanins (endogenously produced, typically cryptic) and carotenoids (dietary uptake, typically bright), providing unique opportunities to address general questions about the biology of coloration, from the mechanisms that determine how different types of pigmentation evolve to discussions on honest signalling hypotheses. Crucial gaps persist in our knowledge on the molecular basis underlying the production and deposition of pterins. We thus highlight the need for functional studies on systems amenable for laboratory manipulation, but also on systems that exhibit natural variation in pterin pigmentation. The wealth of potential model species, coupled with recent technological and analytical advances, make this a promising time to advance research on pterin-based pigmentation in animals.

Transcriptome sequencing and analysis for the pigmentation of scale and skin in common carp (Cyprinus carpio)

BACKGROUND

Teleost scale not only provides a protective layer resisting penetration and pathogens but also participate in coloration. It is interesting to study the mechanism of teleost scale formation. Furthermore, whether there existed consensus genes between scale coloration and skin coloration has not been examined yet.

METHODS AND RESULTS

We analyzed the transcriptome profiles of red scale, white scale, red skin, and white skin of common carp (Cyprinus carpio). Pair-wise comparison identified 3391 differentially expressed genes (DEGs) between scale and skin, respectively. The 1765 up-regulated genes (UEGs) in scale, as the down-regulated genes in skin, preferred mineralization and other scale development-related processes. The 1626 skin UEGs were enriched in the morphogenesis of skin and appendages. We also identified 195 UEGs in white scale and 223 UEGs in red scale. The white scale UEGs primarily participated in regulation of growth and cell migration. The UEGs in red scale preferred pigment cell differentiation and retinoid metabolic process. A total of 22 DEGs had consensus expression patterns in skin and scale of the same coloration. The expression levels of these DEGs clearly grouped skin and scale of the same coloration together with principle component analysis and correlation analysis. Eleven consensus DEGs were homologous to the orthologs of Poropuntius huangchuchieni, 82% of which were under strong purifying selection. Eight processes including lipid storage and lipid catabolism were shared in both scale pigmentation and skin pigmentation.

CONCLUSIONS

We identified consensus DEGs and biological processes in scale and skin pigmentation. Our transcriptome analysis will contribute to further elucidation of mechanisms of teleost scale formation and coloration.

Differences in DNA Methylation Between Disease-Resistant and Disease-Susceptible Chinese Tongue Sole (Cynoglossus semilaevis) Families

DNA methylation, the most widely studied and most well-understood epigenetic modification, has been reported to play crucial roles in diverse processes. Although it has been found that DNA methylation can modulate the expression of immune-related genes in teleosts, a systemic analysis of epigenetic regulation on teleost immunity has rarely been performed. In this research, we employed whole-genome bisulfite sequencing to investigate the genome-wide DNA methylation profiles in select disease-resistant Cynoglossus semilaevis (DR-CS, family 14L006) and disease-susceptible C. semilaevis (DS-CS, family 14L104) against Vibrio harveyi infection. The results showed that following selective breeding, DR-CS had higher DNA methylation levels and different DNA methylation patterns, with 3,311 differentially methylated regions and 6,456 differentially methylated genes. Combining these data with the corresponding transcriptome data, we identified several immune-related genes that exhibited differential expression levels that were modulated by DNA methylation. Specifically, DNA methylation of tumor necrosis factor–like and lipopolysaccharide-binding protein-like was significantly correlated with their expression and significantly contributed to the disease resistance of the selected C. semilaevis family. In conclusion, we suggest that artificial selection for disease resistance in Chinese tongue sole causes changes in DNA methylation levels in important immune-related genes and that these epigenetic changes are potentially involved in multiple immune responses in Chinese tongue sole.

Morphological Characteristics and Comparative Transcriptome Analysis of Three Different Phenotypes of Pristella maxillaris

Pristella maxillaris is known as the X-ray fish based on its translucent body. However, the morphological characteristics and the molecular regulatory mechanisms of these translucent bodies are still unknown. In this study, the following three phenotypes, a black-and-gray body color or wild-type (WT), a silvery-white body color defined as mutant I (MU1), and a fully transparent body with a visible visceral mass named as mutant II (MU2), were investigated to analyze their chromatophores and molecular mechanisms. The variety and distribution of pigment cells in the three phenotypes of P. maxillaris significantly differed by histological assessment. Three types of chromatophores (melanophores, iridophores, and xanthophores) were observed in the WT, whereas MU1 fish were deficient in melanophores, and MU2 fish lacked melanophores and iridophores. Transcriptome sequencing of the skin and peritoneal tissues of P. maxillaris identified a total of 166,089 unigenes. After comparing intergroup gene expression levels, more than 3,000 unigenes with significantly differential expression levels were identified among three strains. Functional annotation and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of the differentially expressed genes (DEGs) identified a number of candidates melanophores and iridophores genes that influence body color. Some DEGs that were identified using transcriptome analysis were confirmed by quantitative real-time PCR. This study serves as a global survey of the morphological characteristics and molecular mechanism of different body colors observed in P. maxillaris and thus provides a valuable theoretical foundation for the molecular regulation of the transparent phenotype.

Integrated analysis of long non-coding RNA and mRNA expression in different colored skin of koi carp

BACKGROUND

Long non-coding RNAs (lncRNAs) perform crucial roles in biological process involving complex mechanisms. However, information regarding their abundance, characteristics and potential functions linked to fish skin color is limited. Herein, Illumina sequencing and bioinformatics were conducted on black, white, and red skin of Koi carp (Cyprinus carpio L.).

RESULTS

A total of 590,415,050 clean reads, 446,614 putative transcripts, 4252 known and 72,907 novel lncRNAs were simultaneously obtained, including 92 significant differentially expressed lncRNAs and 722 mRNAs. Ccr_lnc5622441 and Ccr_lnc765201 were up-regulated in black and red skin, Ccr_lnc14074601 and Ccr_lnc2382951 were up-regulated in white skin, and premelanosome protein a (Pmela), Pmelb and tyrosinase (Tyr) were up-regulated in black skin. The expression patterns of 18 randomly selected differentially expressed genes were validated using the quantitative real-time PCR method. Moreover, 70 lncRNAs acting on 107 target mRNAs in cis and 79 lncRNAs acting on 41,625 target mRNAs in trans were investigated. The resulting co-expression networks revealed that a single lncRNA can connect with numerous mRNAs, and vice versa. To further reveal their potential functions, Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, and membrane, pigment cell development, cAMP signaling, melanogenesis and tyrosine metabolism appear to affect skin pigmentation. Additionally, three lncRNAs (Ccr_lnc142711, Ccr_lnc17214525 and Ccr_lnc14830101) and three mRNAs (Asip, Mitf and Tyr) involved in the melanogenesis pathway were investigated in terms of potential functions in embryogenesis and different tissues.

CONCLUSIONS

The findings broaden our understanding of lncRNAs and skin color genetics, and provide new insight into the mechanisms underlying lncRNA-mediated pigmentation and differentiation in Koi carp.

Identification and characterization of skin color microRNAs in Koi carp (Cyprinus carpio L.) by Illumina sequencing

Background

MicroRNAs (miRNAs) are endogenous, small (21–25 nucleotide), non-coding RNAs that play important roles in numerous biological processes. Koi carp exhibit diverse color patterns, making it an ideal subject for studying the genetics of pigmentation. However, the influence of miRNAs on skin color regulation and variation in Koi carp is poorly understood.

Results

Herein, we performed small RNA (sRNA) analysis of the three main skin colors in Koi carp by Illumina sequencing. The results revealed 330, 397, and 335 conserved miRNAs (belonging to 81 families) and 340, 353, and 351 candidate miRNAs in black, red, and white libraries, respectively. A total of 164 differentially expressed miRNAs (DEMs) and 14 overlapping DEMs were identified, including miR-196a, miR-125b, miR-202, miR-205-5p, miR-200b, and etc. Target prediction and functional analysis of color-related miRNAs such as miR-200b, miR-206, and miR-196a highlighted putative target genes, including Mitf, Mc1r, Foxd3, and Sox10 that are potentially related to pigmentation. Determination of reference miRNAs for relative quantification showed that let-7a was the most abundant single reference gene, and let-7a and miR-26b was the most abundant combination.

Conclusions

The findings provide novel insight into the molecular mechanisms determining skin color differentiation in Koi carp, and serve as a valuable reference for future studies on tissue-specific miRNA abundance in Koi carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5189-5) contains supplementary material, which is available to authorized users.

Dynamic regulation of mRNA and miRNA associated with the developmental stages of skin pigmentation in Japanese ornamental carp

The Japanese ornamental carp (Cyprinus carpio var. Koi) is famous for multifarious colors and patterns, making it commonly culture and trade across the world. Although functional genes and inheritance of color traits have been commonly studied, seldom attentions were focused on the genetic regulation during the developmental process of pigmentation. To better understand the mechanism of skin color development, we observed the morphogenesis of pigment cells during the post-embryonic stages and analysed the temporal expression pattern of mRNAs/miRNAs profiles in four distinct developmental stages. 59 and 103 differentially expressed genes/miRNAs (DEGs/DEMs) associated with pigmentation and skin were identified, including pax7, mitf, tyr, tyrp1, etc., and the highest DEGs were detected at 11 days post hatching (dph). In addition, the functional characteristics of mRNAs/miRNAs associated with pteridine and carotenoid pathway were also examined. Furthermore, 65 miRNA-mRNA interaction pairs related to pigmentation, pteridines and carotenoids metabolism were detected between different stages. Interestingly, the largest pairs appeared in the transition from 11 dph to 48 dph, which had the similar trend with DEGs further manifesting the importance of 11 dph. This study produced a comprehensive programme of DEGs/DEMs during color development, which will provide resources to understand the regulation mechanism in color formation. The understanding of genetic basis in color formation might promote the production and breeding of the Koi carp.

Comparative Transcriptome and DNA methylation analyses of the molecular mechanisms underlying skin color variations in Crucian carp (Carassius carassius L.)

Background

Crucian carp is a popular ornamental strain in Asia with variants in body color. To further explore the genetic mechanisms underlying gray and red body color formation in crucian carp, the skin transcriptomes and partial DNA methylation sites were obtained from red crucian carp (RCC) and white crucian carp (WCC). Here, we show significant differences in mRNA expression and DNA methylation sites between skin tissues of RCC and WCC.

Results

Totals of 3434 and 3683 unigenes had significantly lower and higher expression in WCC, respectively, compared with unigenes expressed in RCC. Some potential genes for body color development were further identified by quantitative polymerase chain reaction, such as mitfa, tyr, tyrp1, and dct, which were down-regulated, and foxd3, hpda, ptps, and gch1, which were up-regulated. A KEGG pathway analysis indicated that the differentially expressed genes were mainly related to mitogen activated protein kinase (MAPK), Wnt, cell cycle, and endocytosis signaling pathways, as well as variations in melanogenesis in crucian carp. In addition, some differentially expressed DNA methylation site genes were related to pigmentation, including mitfa, tyr, dct, foxd3, and hpda. The differentially expressed DNA methylation sites were mainly involved in signaling pathways, including MAPK, cAMP, endocytosis, melanogenesis, and Hippo.

Conclusions

Our study provides the results of comparative transcriptome and DNA methylation analyses between RCC and WCC skin tissues and reveals that the molecular mechanism of body color variation in crucian carp is strongly related to disruptions in gene expression and DNA methylation during pigmentation.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0564-9) contains supplementary material, which is available to authorized users.

Comparative transcriptome and histological analyses provide insights into the prenatal skin pigmentation in goat (Capra hircus)

The Youzhou dark goat is a natural mutant with dark skin over the whole body including the visible mucous membranes. In the present study, we characterized 100-day-old fetal skin at the histomorphological and transcriptomic levels in dark-skinned (Youzhou dark goat) and white-skinned (Yudong white goat) goats with deep RNA sequencing, quantitative PCR, and histological methods. Histological analysis indicated that there were marked differences in both melanin distribution and epidermal ultrastructure between the hyperpigmented and normal skin in two breeds of goat. Subsequent analyses suggested that a presumed structure variation (duplication or insertion) in ASIP might be responsible for its lower expression in the hyperpigmented skin (Youzhou dark goat) by determining the distribution of melanocytes across the body at early development stage. Analyses for genes with differential expression between the dark-skinned and white-skinned goats indicated the network composed of ASIP-MC1R, ECM-receptor interaction, and MAPK signaling might play crucial roles in the determination of skin pigmentation in fetal goats. Moreover, we also identified 1,616 novel transcripts in goat skin by RNA sequencing, which may represent two distinct groups of transcript based on their characteristics. Our findings contribute to the understanding of the characteristics of global gene expression in early-stage skin pigmentation and development and describe an animal model for human diseases associated with pigmentation.

Morphological Characters and Transcriptome Profiles Associated with Black Skin and Red Skin in Crimson Snapper (Lutjanus erythropterus)

In this study, morphology observation and illumina sequencing were performed on two different coloration skins of crimson snapper (Lutjanus erythropterus), the black zone and the red zone. Three types of chromatophores, melanophores, iridophores and xanthophores, were organized in the skins. The main differences between the two colorations were in the amount and distribution of the three chromatophores. After comparing the two transcriptomes, 9200 unigenes with significantly different expressions (ratio change ≥ 2 and q-value ≤ 0.05) were found, of which 5972 were up-regulated in black skin and 3228 were up-regulated in red skin. Through the function annotation, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the differentially transcribed genes, we excavated a number of uncharacterized candidate pigment genes as well as found the conserved genes affecting pigmentation in crimson snapper. The patterns of expression of 14 pigment genes were confirmed by the Quantitative real-time PCR analysis between the two color skins. Overall, this study shows a global survey of the morphological characters and transcriptome analysis of the different coloration skins in crimson snapper, and provides valuable cellular and genetic information to uncover the mechanism of the formation of pigment patterns in snappers.