Comparative transcriptome analysis of molecular mechanism underlying gray-to-red body color formation in red crucian carp (Carassius auratus, red var.)

Involvement of miR-495 in the skin pigmentation of rainbow trout (Oncorhynchus mykiss) through the regulation of mc1r

MicroRNAs (miRNAs) play crucial roles in skin pigmentation in animals. Rainbow trout (Oncorhynchus mykiss) is a key economic fish species worldwide, and skin color directly affects its economic value. However, the functions of miRNAs in rainbow trout skin pigmentation remain largely unknown. Herein, we overexpressed and silenced miR-495 in vitro and in vivo to investigate its functions. The analysis of spatial and temporal expression patterns suggested that miR-495 is a potential regulator during the process of skin pigmentation. In vitro, mc1r was validated as a direct target for miR-495 by dual-luciferase reporter assay, and overexpression of miR-495 significantly inhibited mc1r expression; in contrast, mc1r and its downstream gene mitf levels were markedly upregulated by decreased miR-495. In vivo, overexpressed miR-495 by injecting agomiR-495 led to a substantial decrease in the expression of mc1r and mitf in dorsal skin and liver, while the opposite results were obtained after miR-495 silencing by antagomiR-495. These findings suggested that miR-495 can target mc1r to regulate rainbow trout skin pigmentation, which provide a potential basis for using miRNAs as target drugs to treat pigmentation disorders and melanoma.

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.

Comparative Transcriptome Analysis of the Response to Vibrio parahaemolyticus and Low-Salinity Stress in the Swimming Crab Portunus trituberculatus

Vibrio parahaemolyticus is one of the main pathogenic bacteria of Portunus trituberculatus and causes mass mortality of P. trituberculatus in aquaculture. In addition, low-salinity stimulation makes P. trituberculatus more susceptible to V. parahaemolyticus infections. In order to elucidate the molecular mechanism of resistance to V. parahaemolyticus in P. trituberculatus, comparative transcriptomic analysis of blood cells stimulated by low salinity and V. parahaemolyticus was carried out in this study. Transcriptome sequencing of low-salinity stress and pathogen infection at different time points was completed using Illumina sequencing technology. A total of 5827, 6432, 5362 and 1784 differentially expressed genes (DEGs) involved in pathways related to ion transport and immunoregulation were found under low-salinity stress at 12, 24, 48 and 72 h compared with the control at 0 h. In contrast, 4854, 4814, 5535 and 6051 DEGs, which were significantly enriched in Toll and IMD signaling pathways, were found at 12, 24, 48 and 72 h compared with the control at 0 h under V. parahaemolyticus infection. Among them, 952 DEGs were shared in the two treatment groups, which were mainly involved in apoptosis and Hippo signaling pathway. Cluster analysis screened 103 genes that were differentially expressed in two factors that were negatively correlated, including immunoglobulin, leukocyte receptor cluster family, scavenger receptor, macroglobulin and other innate-immune-related genes. These results provide data support for the analysis of the mechanisms of immunity to V. parahaemolyticus under low-salinity stress in P. trituberculatus and help to elucidate the molecular mechanisms by which environmental factors affect immunity.

Genome-wide association study (GWAS) analysis of black color trait in the leopard coral grouper (Plectropomus leopardus) using whole genome resequencing

The leopard coral grouper (Plectropomus leopardus) is a coral reef fish species that exhibits rapid and diverse color variation. However, the presence of melanoma and the high proportion of individuals displaying black color in artificial breeding have led to reduced economic and ornamental value. To pinpoint single nucleotide polymorphisms (SNPs) and potential genes linked to the black pigmentation characteristic in this particular species, This study gathered a cohort of 360 specimens from diverse origins and conducted a comprehensive genome-wide association analysis (GWAS) employing whole-genome resequencing. As a result, 57 SNPs related to the black skin trait were identified, and a grand total of 158 genes were annotated within 50 kb of these SNPs. Subsequently, GWAS was applied to three populations (LED, QHH, and QHL), and the corresponding results were compared with the analysis results of the total population. The results of the four GWAS models showed significant enrichment in Rap1 signaling pathway, melanin biosynthesis, metabolic pathways, tyrosine metabolism, cAMP signaling pathway, AMPK signaling pathway, PI3K-Akt signaling pathway, EGFR tyrosine kinase inhibitor resistance, HIF-1 signaling pathway, Ras signaling pathway, MAPK signaling pathway, etc. (p < 0.05), which were mainly associated with eleven genes (POL4, MET, E2F2, COMT, ZBED1, TYRP2, FOXP2, THIKA, LORF2, MYH16 and SOX2). Significant differences (p < 0.05) were observed in the expression of all 11 genes in the dorsal skin tissue, in 10 genes except COMT in the ventral skin tissue, and in all 11 genes in the caudal fin tissue. These findings imply that the control of body color in the P. leopardus is the result of the joint action of multiple genes and signaling pathways. These findings will contribute to a more profound comprehension of the genetic attributes that underlie the development of black skin in the vibrant P. leopardus, thus furnishing a theoretical foundation for genetic enhancement.

Pigment Identification and Gene Expression Analysis during Erythrophore Development in Spotted Scat (Scatophagus argus) Larvae

Red coloration is considered an economically important trait in some fish species, including spotted scat, a marine aquaculture fish. Erythrophores are gradually covered by melanophores from the embryonic stage. Despite studies of black spot formation and melanophore coloration in the species, little is known about erythrophore development, which is responsible for red coloration. 1-phenyl 2-thiourea (PTU) is a tyrosinase inhibitor commonly used to inhibit melanogenesis and contribute to the visualization of embryonic development. In this study, spotted scat embryos were treated with 0.003% PTU from 0 to 72 h post fertilization (hpf) to inhibit melanin. Erythrophores were clearly observed during the embryonic stage from 14 to 72 hpf, showing an initial increase (14 to 36 hpf), followed by a gradual decrease (36 to 72 hpf). The number and size of erythrophores at 36 hpf were larger than those at 24 and 72 hpf. At 36 hpf, LC-MS and absorbance spectrophotometry revealed that the carotenoid content was eight times higher than the pteridine content, and β-carotene and lutein were the main pigments related to red coloration in spotted scat larvae. Compared with their expression in the normal hatching group, rlbp1b, rbp1.1, and rpe65a related to retinol metabolism and soat2 and apoa1 related to steroid hormone biosynthesis and steroid biosynthesis were significantly up-regulated in the PTU group, and rh2 associated with phototransduction was significantly down-regulated. By qRT-PCR, the expression levels of genes involved in carotenoid metabolism (scarb1, plin6, plin2, apoda, bco1, and rep65a), pteridine synthesis (gch2), and chromatophore differentiation (slc2a15b and csf1ra) were significantly higher at 36 hpf than at 24 hpf and 72 hpf, except for bco1. These gene expression profiles were consistent with the developmental changes of erythrophores. These findings provide insights into pigment cell differentiation and gene function in the regulation of red coloration and contribute to selective breeding programs for ornamental aquatic animals.

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.

Comparative transcriptome and miRNA analysis of skin pigmentation during embryonic development of Chinese soft-shelled turtle (Pelodiscus sinensis)

BACKGROUND

Aquatic animals show diverse body coloration, and the formation of animal body colour is a complicated process. Increasing evidence has shown that microRNAs (miRNAs) play important regulatory roles in many life processes. The role of miRNAs in pigmentation has been investigated in some species. However, the regulatory patterns of miRNAs in reptile pigmentation remain to be elucidated. In this study, we performed an integrated analysis of miRNA and mRNA expression profiles to explore corresponding regulatory patterns in embryonic body colour formation in the soft-shelled turtle Pelodiscus sinensis.

RESULTS

We identified 8 866 novel genes and 9 061 mature miRNAs in the skin of Chinese soft-shelled turtles in three embryonic stages (initial period: IP, middle period: MP, final period: FP). A total of 16 563 target genes of the miRNAs were identified. Furthermore, we identified 2 867, 1 840 and 4 290 different expression genes (DEGs) and 227, 158 and 678 different expression miRNAs (DEMs) in IP vs. MP, MP vs. FP, and IP vs. FP, respectively. Among which 72 genes and 25 miRNAs may be related to turtle pigmentation in embryonic development. Further analysis of the novel miRNA families revealed that some novel miRNAs related to pigmentation belong to the miR-7386, miR-138, miR-19 and miR-129 families. Novel_miR_2622 and novel_miR_2173 belong to the miR-19 family and target Kit and Gpnmb, respectively. The quantification of novel_miR_2622 and Kit revealed negative regulation, indicating that novel_miR_2622 may participate in embryonic pigmentation in P. sinensis by negatively regulating the expression of Kit.

CONCLUSIONS

miRNA act as master regulators of biological processes by controlling the expression of mRNAs. Considering their importance, the identified miRNAs and their target genes in Chinese soft-shelled turtle might be useful for investigating the molecular processes involved in pigmentation. All the results of this study may aid in the improvement of P. sinensis breeding traits for aquaculture.

MiR-382 Functions on the Regulation of Melanogenesis via Targeting dct in Rainbow Trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are a class of endogenous small noncoding RNAs that are involved in a diverse collection of biological processes as important post-transcriptional regulators. However, little is known about the molecular regulatory mechanism of miRNAs in fish skin pigmentation. In this study, we first confirmed that dopachrome tautomerase (dct), a key gene of melanogenesis, is a target of miR-382 in rainbow trout (Oncorhynchus mykiss) using luciferase reporter assay. The analysis of different developmental stages and tissue expression patterns between wild-type and yellow mutant rainbow trout suggested that miR-382 is a potential regulator during the process of skin pigmentation. In vitro, miR-382 mimics in rainbow trout primary liver cells significantly downregulated dct expression and resulted in decreased expression of key melanogenic genes including tyrosine-related protein 1 (tyrp1) and premelanosome protein (pmel), whereas the expression level of dct was markedly increased after transfected with miR-382 inhibitor. In vivo, overexpression of miR-382 by injection of miR-382 agomir significantly depressed the expression of dct in dorsal skin, tail fin, and liver and then reduced the expression levels of tyrp1 and pmel. Furthermore, transfection of miR-382 mimics inhibited cell proliferation and induced apoptosis. Taken together, our results identified a functional role of miR-382 in rainbow trout skin pigmentation through targeting dct, which facilitate understanding the regulatory mechanism of rainbow trout skin color at the post-transcriptional level and provide a theoretical basis for molecular breeding with skin color as the target trait.

Transcriptomic and proteomic analyses reveal the common and unique pathway(s) underlying different skin colors of leopard coral grouper (Plectropomus leopardus)

To gain a comprehensive and unbiased molecular understanding of the different skin colors of P. leopardus, we used Illumina HiSeq 2500 and TMT (Tandem Mass Tag) to compare transcription and protein levels between red and black skin of P. leopardus. We identified 797 upregulated and 314 downregulated genes (differentially expressed genes; DEGs) in red (RG) compared with black (BG) skin of P. leopardus. We also identified 377 differentially abundant proteins (DAPs), including 314 upregulated and 63 downregulated proteins. These DEGs and DAPs were significantly enriched in melanin synthesis (e.g., pyrimidine metabolism, Phenylalanine, tyrosine, and tryptophan biosynthesis, melanogenesis, phenylalanine metabolism, and tyrosine metabolism), oxidative phosphorylation (e.g., phosphonate and phosphinate metabolism, and oxidative phosphorylation), energy metabolism (e.g., HIF-1, glycolysis/gluconeogenesis, fatty acid biosynthesis, and fatty acid degradation), and signal transduction (e.g., Wnt, calcium, MAPK, and cGMP-PKG signaling pathways), etc. Further analysis of MAPKs showed that the activation levels of its main members JNK1 and ERK1/2 differed significantly between red and black skin colors. After RNAi was used to interfere with ERK1/2, it was found that the local skin of the tail of P. leopardus would turn black. Combined transcriptome and proteome analysis showed that most DEGs-DAPs in red skin were higher than in black skin (58 were upregulated, 1 was downregulated, and 4 were opposite). These DEGs-DAPs showed that the differences between red and black skin tissues of P. leopardus were related primarily to energy metabolism, signal transduction and cytoskeleton. These findings are not only conducive to understand the skin color regulation mechanism of P. leopardus and other coral reef fish, but also provide an important descriptive to the breeding of color strains. SIGNIFICANCE OF THE STUDY: The skin color of P. leopardus gradually darkens or blackens due to environmental factors such as changes in light intensity and human activities, and this directly affects its ornamental and economic value. In this study, RNAseq and TMT were used to conduct comparative quantitative transcriptomics and proteomics and analyze differences between red and black P. leopardus skin. The results showed that energy metabolism, signal transduction and cytoskeleton were the main metabolic pathways causing their skin color differences. These findings contribute to existing data describing fish skin color, and provide information about protein levels, which are of great significance to a deeper understanding of the skin color regulation mechanism in P. leopardus and other coral reef fishes.

Molecular and functional analysis of the microphthalmia-associated transcription factor (mitf) gene duplicates in red tilapia

In vertebrates, the microphthalmia-associated transcription factor (mitf) is at the hub of the melanin synthesis regulation network. However, little information is known about its molecular characterization, expression, location, or function in skin color differentiation and variation of red tilapia. The full-length cDNA sequences (1977 bp and 1999 bp) of mitfa and mitfb, encoding polypeptides of 491 and 514 amino acids, were effectively identified from red tilapia in this study. The Mitfa and Mitfb sequences of red tilapia clustered first with O. aureus, then with other teleost fish, according to phylogenetic analysis. Mitfa and mitfb mRNA were highly expressed in the brain, dorsal skin and eye tissues using quantitative real-time PCR. The mRNA expressions of mitfa and mitfb were the highest in the cleavage stage during the early development of red tilapia. Among three different colors of red tilapia, the expression levels of mitfa and mitfb were highest in the PB (pink with scattered black spots) dorsal skin. After overwintering, the mitfa and mitfb mRNA expressions were high in the dorsal skin of PB (color changed from pink to black). Mitfa and mitfb were mostly found in the epidermal layer of the dorsal skin, according to in situ hybridization (ISH) analysis. After injecting mitf-dsRNA duplicates along the tail vein of red tilapia, the activity of tyrosinase and the level of melanin in the dorsal skin both decreased significantly. The mRNA expressions of mitfa and its downstream genes (tyrb, tyrp1a and dct) decreased, whereas the mRNA expression of mitfb increased after mitfa-dsRNA injection. The mRNA expressions of mitfb, tyrb, tyrp1a and dct decreased, whereas the mRNA expression of mitfa increased after injecting mitfb-dsRNA. These findings suggest that mitf gene duplicates may play an important role in red tilapia skin color differentiation and variation via the melanogenesis pathway.

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.

CRISPR knockouts of pmela and pmelb engineered a golden tilapia by regulating relative pigment cell abundance

Premelanosome protein (pmel) is a key gene for melanogenesis. Mutations in this gene are responsible for white plumage in chicken, but its role in pigmentation of fish remains to be demonstrated. In this study we found that most fishes have two pmel genes arising from the teleost-specific whole genome duplication. Both pmela and pmelb were expressed at high levels in the eyes and skin of Nile tilapia. We mutated both genes in tilapia using CRISPR/Cas9. Homozygous mutation of pmela resulted in yellowish body color with weak vertical bars and a hypo-pigmented retinal pigment epithelium (RPE) due to significantly reduced number and size of melanophores. In contrast, we observed an increased number and size of xanthophores in mutants compared to wild-type fish. Homozygous mutation of pmelb resulted in a similar, but milder phenotype than pmela-/- mutants. Double mutation of pmela and pmelb resulted in loss of additional melanophores compared to the pmela-/- mutants, and also an increase in the number and size of xanthophores, producing a golden body color. The RPE pigmentation of pmela-/-;pmelb-/- was similar to pmela-/- mutants, with much less pigmentation than pmelb-/- mutants and wild-type fish. Taken together, our results indicate that, while both pmel genes are important for the formation of body color in tilapia, pmela plays a more important role than pmelb. To our knowledge, this is the first report on mutation of pmelb or both pmela;pmelb in fish. Studies on these mutants suggest new strategies for breeding golden tilapia, and also provide a new model for studies of pmel function in vertebrates.

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.

Water volume influences antibiotic resistomes and microbiomes during fish corpse decomposition

Corpse decomposition may cause serious pollution (e.g., releasing antibiotic resistance genes) to the water environment, thereby threatening public health. However, whether antibiotic resistance genes (ARGs) and microbiomes are affected by different water volumes during carcass decomposition remains unknown. Here, we investigated the effects of large/small water volumes on microbial communities and ARGs during fish cadaver decomposition by 16S rRNA high-throughput sequencing and high-throughput quantitative PCR. The results showed that the large water volume almost eliminated the effects of corpse decomposition on pH, total organic carbon (TOC), and total nitrogen (TN). When the water volume enlarged by 62.5 fold, the relative abundances of some ARGs resisting tetracycline and sulfonamide during carcass decomposition decreased by 217 fold on average, while there was also a mean 5267 fold increase of vancomycin resistance genes. Compared with the control group, the enriched types of ARGs varied between the large and small volume. Water volume, mobile genetic elements, and carcass decomposition were the most important factors affecting ARG profiles. Many opportunistic pathogens (like Bacteroides and Comamonas) were enriched in the corpse group. Bacteroides and Comamonas may be potential hosts of ARGs, indicating the potential for the spread of ARGs to humans by water pathogenic bacteria. This research highlights that the "dilution effect" can contribute to eliminating this adverse effect during corpse decomposition to a certain extent. It may provide references for environmental governance and public health.

Appropriate leucine supplementation promotes glucose metabolism and enhances energy homeostasis in juvenile crucian carp (Carassius auratus gibelio var. CAS III)

In order to characterize the molecular mechanisms by which leucine regulates carbohydrate metabolism and energy homeostasis, juvenile crucian carps (Carassius auratus gibelio var. CAS III) fed with a high carbohydrate diet were supplemented with different levels of dietary leucine: 0% (Leu0), 0.4% (Leu4), 0.8% (Leu8), 1.2% (Leu12), 1.6% (Leu16), 2.0% (Leu20), and 5.0% (Leu50). After 8 weeks, RNA sequencing was performed on samples collected from the Leu0, Leu8, Leu12 and Leu50 groups. Differentially expressed genes were then detected and analyzed. The results showed a total of 91.6 Gb of clean bases were generated. Moreover, a total of 1131, 5254, and 1539 DEGs were detected in Leu8, Leu12, and Leu50 compared with Leu0, respectively, encompassing 161 common DEGs. STEM analysis elucidated four significant clusters of DEGs that were associated with "glycerophospholipid metabolism," "glycerolipid metabolism," "PPAR signaling pathway," and "adipocytokine signaling pathway." Moreover, the mRNA expression levels of acyl-CoA synthetase long chain family member 5 (ACSL5), choline kinase beta (CHKB), cryptochrome-1 (CRY1), lon protease homolog 2, peroxisomal isoform X2 (LONP2), lipin 1 (LPIN1), membrane bound O-acyltransferase domain containing 2 (MBOAT2), phosphoenolpyruvate carboxykinase 1 (PEPCK), and uridine-cytidine kinase 2b (UCK2b) were then further investigated in all leucine treatment groups at starvation times of 0 h, 24 h, and 48 h. The results revealed that the expression levels of UCK2b and MBOAT2 were negatively correlated with the addition of leucine, whereas CHKB, LONP2, CRY1, PEPCK, and LPIN1 were positively correlated. In conclusion, dietary leucine supplementation below 1.2% enhanced carbohydrate metabolism in juvenile crucian carp fed with a high-carbohydrate diet, whereas concentrations above 2.0% is a better choice for energy homeostasis under starvation.

Brightened body coloration in female guppies (Poecilia reticulata) serves as an in vivo biomarker for environmental androgens: The example of 17β-trenbolone

In vivo testing systems for environmental androgens are scarce. The aim of this study was to evaluate the potential of male-specific brightened body coloration in female guppies (Poecilia reticulata) to serve as an in vivo biomarker of environmental androgens using 17β-trenbolone as an example. The high bioaccumulation of 17β-trenbolone in the skin of female guppies suggests that it is a potential target tissue of environmental androgens. The coloration index, pigment cell ultrastructure, pigment levels, sexual attractiveness, and reproductive capability of female guppies were analyzed following 28 days of exposure to 20 ng/L, 200 ng/L, and 2000 ng/L 17β-trenbolone. Increases in the coloration index caused by 17β-trenbolone exposure were attributable to increased pteridine and melanin levels. Decreases in the sexual attractiveness, number of offspring, and survival rate of offspring suggested that the changes in body coloration translated into adverse outcomes. Finally, mRNA sequencing indicated that 17β-trenbolone increased pteridine levels by activating genomic effects of androgen receptor on xanthine dehydrogenase and increased melanin levels by exerting non-genomic effects targeting microphthalmia-associated transcription factor, tyrosinase, and tyrosinase-related protein 1 that were mediated by mitogen-activated protein kinase and calcium signaling pathways. We have derived a robust adverse outcome pathway of environmental androgens, and our findings suggest that indicators at different biological levels related to brightened body coloration in female guppies can serve as less-invasive or noninvasive in vivo biomarkers of short-term exposure to environmental androgens.

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.

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.

Skin transcriptome, tissue distribution of mucin genes and discovery of simple sequence repeats in crucian carp (Carassius auratus)

Crucian carp (Carassius auratus) is one of the major freshwater species and important food fish in China. Fish skin acts as the first line of defence against pathogens, yet its molecular and immune mechanism remains unclear. In this study, a de novo transcriptome assembly of C. auratus skin was performed with the Illumina Hiseq 2000 platform. A total of 49,154,776 clean reads were assembled, among which 60,824 (46.86%), 37,103 (28.59%), 43,269 (33.33%) unigenes were annotated against National Center for Biotechnology Information, Gene Onotology and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. KEGG Orthology categories were significantly involved in immune system (20.50%), signal transduction (18.04%) and mucosal mucin genes (e.g., muc2, muc5AC, muc5B, muc17, muc18). The high expression of muc18 gene was observed in brain; that of muc2 in intestine; and that of muc5AC in skin, liver, spleen, intestine and muscle. Moreover, the potential 28,928 simple sequence repeats with the three most abundant dinucleotide repeat motifs (AC/GT, AG/CT, AT/AT) were detected in C. auratus. To authors' knowledge, this is the first report to describe the transcriptome analysis of C. auratus skin, and the outcome of this study contributed to the understanding of mucosal immune response of the skin and molecular markers in cyprinid species.

The Genetic Basis of Morphological Diversity in Domesticated Goldfish

Although domesticated goldfish strains exhibit highly diversified phenotypes in morphology, the genetic basis underlying these phenotypes is poorly understood. Here, based on analysis of transposable elements in the allotetraploid goldfish genome, we found that its two subgenomes have evolved asymmetrically since a whole-genome duplication event in the ancestor of goldfish and common carp. We conducted whole-genome sequencing of 27 domesticated goldfish strains and wild goldfish. We identified more than 60 million genetic variations and established a population genetic structure of major goldfish strains. Genome-wide association studies and analysis of strain-specific variants revealed genetic loci associated with several goldfish phenotypes, including dorsal fin loss, long-tail, telescope-eye, albinism, and heart-shaped tail. Our results suggest that accumulated mutations in the asymmetrically evolved subgenomes led to generation of diverse phenotypes in the goldfish domestication history. This study is a key resource for understanding the genetic basis of phenotypic diversity among goldfish strains.

Transcriptome analysis identifies immune-related genes and antimicrobial peptides in Siamese fighting fish (Betta splendens)

Siamese fighting fish (Betta splendens) is one of the most widely cultivated ornamental fish in global trade. However, transcriptomic data, which can reveal valuable genetic data for disease control and prevention, are extremely limited for this species. In this study, whole-body transcriptome sequencing of juvenile betta fish generated 4.457 GB of clean data and a total of 71,775 unigenes using the Illumina HiSeq4000 platform. These unigenes were functionally classified using 7 functional databases, yielding 45,316 NR (63.14%), 47,287 NT (65.88%), 39,105 Swiss-Prot (54.48%), 16,492 COG (22.98%), 37,694 KEGG (52.52%), 4,506 GO (6.28%), and 35,374 Interpro (49.28%) annotated unigenes. Furthermore, we also detected 13,834 SSRs distributed on 10,636 unigenes and 49,589 predicted CDSs. Based on KEGG analysis, five innate immune pathways (997 unigenes) were reported, including the NOD-like receptor signaling pathway, complement and coagulation cascades, toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway and cytosolic DNA-sensing pathway. Moreover, four antimicrobial peptide (AMP) families (hepcidin, piscidin, LEAP-2, and defensins) from the betta fish transcriptome were also identified. Additionally, cDNA and genomic DNA of two β-defensins was successfully isolated from four betta fish species. RT-PCR analysis showed that BsBD1 transcripts were most abundant in the muscle and kidney and BsBD2 transcripts were most abundant in the gill. The genomic organization showed that the BD1 and BD2 genes consisted of three exons and two introns according to the GT-AG rule. Most importantly, this is the first report of the betta fish whole-body transcriptome obtained by high-throughput sequencing. Our transcriptomic data and the discovery of betta fish AMPs should promote a better understanding of molecular immunology for disease prevention for further ornamental fish aquaculture.

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.

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.)

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

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.