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

Identification of pigmentation genes in skin, muscle and tail of a Thai-flag variety of Siamese fighting fish Betta splendens

Pigmentation genes expressed in skin, body muscle and tail of Thai-flag compared with Blue, White and Red varieties of Siamese fighting fish Betta splendens were identified. In total, 22,919 new unigenes were found. Pearson correlation and PCA analysis revealed that expression profiles of genes in muscle, skin and tail across solid color variety were similar. In contrast, those in skin and red tail part of Thai-flag were closely related but they showed different expression profiles with the white tail part. Moreover, 21,347-64,965 SNPs were identified in exonic regions of identified genes. In total, 28,899 genes were differentially expressed between paired comparisons of libraries where 13,907 genes (48.12 %) were upregulated and 14,992 genes (51.88 %) were downregulated. DEGs between paired libraries were 106-5775 genes relative to the compared libraries (56-2982 and 50-2782 for upregulated and downregulated DEGs). Interestingly, 432 pigmentation genes of B. splendens were found. Of these, 297 DEGs showed differential expression between varieties. Many DEGs in melanogenesis (Bsmcr1r, Bsmcr5r, and Bsslc2a15b), tyrosine metabolism (Bstyr, Bstyrp1b and Bsdct), stripe repressor (BsAsip1 and BsAsip2b), pteridine (Bsgch2) and carotenoid (BsBco2) biosynthesis were downregulated in the Thai-flag compared with solid color varieties. Expression of Bsbco1l, Bsfrem2b, Bskcnj13, Bszic2a and Bspah in skin, muscle and tail of Thai-flag, Blue, Red and White varieties was analyzed by qRT-PCR and revealed differential expression between fish varieties and showed anatomical tissue-preferred expression patterns in the same fish variety. The information could be applied to assist genetic-based development of new B. splendens varieties in the future.

Comparative Transcriptome Analysis of the Skin and the Peritoneal Wall Layer of Triplophysa stenura Distributed in High Elevations

A total of 81,868 All-Unigenes were sequenced and assembled by the transcriptome in the dorsal skin, the lateral skin, and the peritoneal wall layer of Triplophysa stenura with a total assembly length of 123,827,585 bp, and 68,750 unigenes were annotated to seven functional databases. A total of 588 DEGs were screened between the dorsal and lateral skin, 17,097 DEGs were screened between the dorsal skin and the peritoneal wall layer, and 16,598 DEGs were screened between the lateral skin and the peritoneal wall layer. Most of DEGs in three tissues were annotated to GO terms related to cellular structures, binding, cellular processes, and catalytic activity. They were also annotated to KEGG pathways such as the MAPK signaling pathway, PI3K-Akt signaling pathway, Wnt signaling pathway, melanogenesis, tyrosine metabolism, and cell cycle. A total of twenty-three DEGs were found to be enriched in the melanin synthesis pathway by a local Blast comparison, of which nine DEGs were significantly upregulated in the peritoneal wall layer and six DEGs were significantly upregulated in the dorsal and lateral skin. The results suggest that these genes may be associated with the molecular mechanism of melanin synthesis in T. stenura, and the differential regulation of genes may be related to the differences in UVR intensity and tissue sites of melanin synthesis. Further investigation is needed on how these genes specifically regulate melanin synthesis.

Whole Genome Resequencing Reveals Selection Signals Related to Wool Color in Sheep

Wool color is controlled by a variety of genes. Although the gene regulation of some wool colors has been studied in relative depth, there may still be unknown genetic variants and control genes for some colors or different breeds of wool that need to be identified and recognized by whole genome resequencing. Therefore, we used whole genome resequencing data to compare and analyze sheep populations of different breeds by population differentiation index and nucleotide diversity ratios (Fst and θπ ratio) as well as extended haplotype purity between populations (XP-EHH) to reveal selection signals related to wool coloration in sheep. Screening in the non-white wool color group (G1 vs. G2) yielded 365 candidate genes, among which PDE4B, GMDS, GATA1, RCOR1, MAPK4, SLC36A1, and PPP3CA were associated with the formation of non-white wool; an enrichment analysis of the candidate genes yielded 21 significant GO terms and 49 significant KEGG pathways (p < 0.05), among which 17 GO terms and 21 KEGG pathways were associated with the formation of non-white wool. Screening in the white wool color group (G2 vs. G1) yielded 214 candidate genes, including ABCD4, VSX2, ITCH, NNT, POLA1, IGF1R, HOXA10, and DAO, which were associated with the formation of white wool; an enrichment analysis of the candidate genes revealed 9 significant GO-enriched pathways and 19 significant KEGG pathways (p < 0.05), including 5 GO terms and 12 KEGG pathways associated with the formation of white wool. In addition to furthering our understanding of wool color genetics, this research is important for breeding purposes.

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 pigment cell composition, pigment content, tyrosinase content and activity of three kinds of loaches Misgurnus anguillicaudatus from Poyang Lake

Pigment cell composition, pigment content, tyrosinase content and activity analysis were investigated on three kinds of loaches Misgurnus anguillicaudatus: big blackspot loaches (BBL), small blackspot loaches (SBL) and non-blackspot loaches (NBL), from Poyang Lake. Results showed that there were three types of skin pigment cells, namely melanophores, xanthophores and iridophores. Melanophores in dorsum were more than those in abdomen. Melanophore cytosomes in BBL were larger than those in SBL and NBL, and melanosomes were the largest in stage four. The melanophores in dorsal skin of SBL or NBL were small cell bodies, spindle-like and in chain distribution. There was an extremely significant difference in melanin content in BBL between the dorsum and abdomen (P < 0.01). There were no significant differences in melanin abdominal content, lutein and carotenoid contents among three kinds of loaches (P > 0.05). In dorsal skin, tyrosinase content was the highest in BBL, and it was significantly lower in NBL than in BBL and SBL (P < 0.01). This study reveals the differences in pigment and tyrosinase content in three kinds of loaches and provides a theoretical basis for further study of the mechanism of black spot formation.

Nile Tilapia: A Model for Studying Teleost Color Patterns

The diverse color patterns of cichlid fishes play an important role in mate choice and speciation. Here we develop the Nile tilapia (Oreochromis niloticus) as a model system for studying the developmental genetics of cichlid color patterns. We identified 4 types of pigment cells: melanophores, xanthophores, iridophores and erythrophores, and characterized their first appearance in wild-type fish. We mutated 25 genes involved in melanogenesis, pteridine metabolism, and the carotenoid absorption and cleavage pathways. Among the 25 mutated genes, 13 genes had a phenotype in both the F0 and F2 generations. None of F1 heterozygotes had phenotype. By comparing the color pattern of our mutants with that of red tilapia (Oreochromis spp), a natural mutant produced during hybridization of tilapia species, we found that the pigmentation of the body and eye is controlled by different genes. Previously studied genes like mitf, kita/kitlga, pmel, tyrb, hps4, gch2, csf1ra, pax7b, and bco2b were proved to be of great significance for color patterning in tilapia. Our results suggested that tilapia, a fish with 4 types of pigment cells and a vertically barred wild-type color pattern, together with various natural and artificially induced color gene mutants, can serve as an excellent model system for study color patterning in vertebrates.

Insights Into the Albinism Mechanism for Two Distinct Color Morphs of Northern Snakehead, Channa argus Through Histological and Transcriptome Analyses

The great northern snakehead (Channa argus) is one of the most important economic and conservational fish in China. In this study, the melanocytes in the skin of two distinct color morphs C. argus were investigated and compared through employment of the microscopic analysis, hematoxylin and eosin (H&E) and Masson Fontana staining. Our results demonstrated the uneven distribution of melanocytes with extremely low density and most of them were in the state of aging or death. Meanwhile, there was no obvious pigment layer and melanocytes distribution pattern found in the albino-type (AT), while the melanocytes were evenly distributed with abundance in the bicolor-type (BT). The transcriptome analysis through Illumina HiSeq sequencing showed that a total of 34.93 Gb Clean Data was obtained, and Q30 base percentage reached 92.66%. The BT and AT northern snakeheads transcriptome data included a total of 56,039,701 and 60,410,063 clean reads (n = 3), respectively. In gene expression analyses, the sample correlation coefficients (r) were ranged between 0.92 and 1.00; the contribution of PC1 and PC2 were 50.25 and 13.73% by using PCA cluster analysis, the total number of DEGs were 1024 (559 up-regulated and 465 down-regulated), and the number of annotated DEGs was 767 (COG 172, KEGG 262, GO 288, SwissProt 548, Pfam 579 and NR 765). Additionally, 46,363 ± 873 and 44,947 ± 392 single nucleotide polymorphisms (SNPs) were compiled via genetic structure analysis, respectively. Ten key pigment-related genes were screened using qRT-PCR. And all of them revealed extremely higher expression levels in the skin of BT than those of AT. This is the first study to analyze the mechanism of albino characteristics of Channa via histology and transcriptomics, and also provide the oretical and practical support for the protection and development of germplasm resources for C. argus.