Transcriptome analysis identifies candidate genes associated with skin color variation in Triplophysa siluroides

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.

An integrative investigation of developmental toxicities induced by triphenyltin in a larval coral reef fish, Amphiprion ocellaris

Triphenyltin (TPT) is widely distributed on coastlines, which makes coral reef fish a potential target of TPT pollution. However, the negative effects of TPT on coral reef fish remain poorly understood. Therefore, in the present study, the larval coral reef fish Amphiprion ocellaris was used to investigate the developmental toxicities of TPT at environmentally relevant concentrations (0, 1, 10 and 100 ng/L). After TPT exposure for 14 d, the cumulative mortality increased, and growth was suppressed. In addition, TPT exposure inhibited the development of melanophores and xanthophores and delayed white strip formation, which might be responsible for the disruption of the genes (erbb3b, mitfa, kit, xdh, tyr, oca2, itk and trim33) related to pigmentation. TPT exposure also attenuated ossification of head skeletal elements and the vertebral column and inhibited the expression of genes (bmp2, bmp4 and sp7) related to skeletal development. The observed developmental toxicities on growth, pigmentation and skeleton development might be associated with the disruption of thyroid hormones and the genes related to thyroid hormone regulation (tshβ, thrα, thrβ, tg, tpo, dio2, and ttr). In addition, TPT exposure interfered with locomotor and shoaling behavior, and the related genes dbh, avp and avpr1aa. Taken together, our results suggest that TPT pollution might threaten the development of one of the most iconic coral reef fish, which might produce disastrous consequences on the health of coral reef ecosystems.

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.

Distinct skin morphological and transcriptomic profiles between wild and albino Oscar Astronotus ocellatus

Oscar Astronotus ocellatus is an important ornamental fish, including albino and wild varieties. Albino individuals attract aquarium hobbyists due to their unique body color, but studies on the species' albinism mechanism are currently scarce. Here, we investigated the morphological and transcriptomic profiles of the skin of albino and wild Oscar. The results showed that the albino type had fewer oval-shaped melanophores and immature melanosomes but that the wild type contained more stellate-shaped melanophores and mature melanosomes. Albino Oscar had a degenerative pigment layer without obvious melanin deposition and content, while the wild type contained more concentrated melanin within the pigment layer. A total of 272,392 unigenes were detected, 109 of which were identified as differentially expressed genes (DEGs) between albino and wild Oscar. Pathways of DEGs, including those involved in complement and coagulation cascades, novobiocin biosynthesis, Th1 and Th2 cell differentiation, and tropane, piperidine and pyridine alkaloid biosynthesis, were significantly enriched. DEGs, including upregulated Sfrp5 and Tat, and downregulated Wnt-10a, Ppp3c, Notch1 and Trim27 involved in the Wnt signaling pathway, Notch signaling pathway, tyrosine metabolism, MAPK signaling pathway and melanogenesis, might be associated with the albinism of Oscar. This study characterized the difference in melanophore morphology between wild and albino Oscar and identified some albinism-related candidate genes and signaling pathways, helping to understand the genetic mechanism of fish albinism.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

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.

Integrative mRNA-miRNA interaction analysis reveals the molecular mechanism of skin color variation between wild-type and yellow mutant rainbow trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss) is an important economic fish in China. Skin color affects the economic value of trout. However, the molecular mechanism of the skin color variation between wild-type (WR) and yellow mutant rainbow trout (YR) is unclear. We sequenced mRNAs and miRNAs of dorsal skin to identify key color variation-associated mRNAs and miRNAs between WR and YR. Overall, 2060 out of 3625 differentially expressed genes were upregulated in YR, and 196 out of 275 differentially expressed miRNAs were downregulated in WR. We identified three key YR-upregulated genes related to the formation of xanthophores (GCH1, SLC2A11, and SOX10). Interestingly, several genes related to melanogenesis (TYR, TYRP1, TYRP2, MC1R, MITF, PMEL, SLC45A2, and OCA2) were downregulated in WR. Integrated analysis identified five miRNAs that target at least two skin color-related genes (miR-495-y, miR-543-y, miR-665-z, miR-433-y, and miR-382-x). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses of target genes identified noncoding RNA metabolic process as the most significantly enriched GO term, and several metabolic pathways associated with skin color were enriched significantly, such as tyrosine metabolism, histidine metabolism, and vitamin B6 metabolism. Quantitative real-time PCR of selected mRNAs and miRNAs validated the reliability of the integrated analysis. This study provides in-depth insights into the molecular mechanism of skin color variation between WR and YR, which will accelerate the genetic selection and breeding of rainbow trout with consumer-favored traits.

Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish

Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.