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

Transcriptome analysis reveals hypoxic response key genes and modules as well as adaptive mechanism of crucian carp (Carassius auratus) gill under hypoxic stress

Fish gill tissue is a primary organ responsive to acute oxygen deprivation or dissolved oxygen (DO) fluctuations in aquatic environments. However, the adaptive mechanism of crucian carp to hypoxic stress remains largely unknown. Here, we investigated gill physiological and transcriptomic changes of crucian carp exposed to hypoxic conditions (dissolved oxygen concentration of 0.6 ± 0.3 mg/L) for different durations (0 d, 1 d, 2d, 3d, 4 d, and 5d). Transcriptomic analysis revealed that the hypoxia group (0.6 ± 0.3 mg/L DO) exhibited a reduction in interlamellar cell mass (ILCM) on the gill filaments, compared with the control group (6.6 ± 0.3 mg/L DO). With prolonged hypoxia stress, the epithelial cells in the gill lamellae became sparse at 3 d to 5 d, and gill vacuoles were increased. A total of 3,502 differentially expressed genes (DEGs) were identified, and 3 hypoxia-specific modules were screened through differential expression analysis, weighted gene co-expression network analysis (WGCNA), and Bayesian network analysis. The apoptosis, necroptosis, efferocytosis and FoxO signaling pathways were significantly enriched based on the KEGG enrichment pathway analysis. The VEGF pathway genes are significantly expressed, enhancing the generation of microvessels in the gill filaments, and improving the capacity to carry oxygen, thus enabling the crucian carp to adapt to hypoxia stress. Hypoxia activated glycolysis, enhanced anaerobic metabolism, promoted β-oxidation of fatty acids, providing energy and maintaining normal physiological metabolism, eventually improving antioxidant and immune capabilities in crucian carp. In summary, this study reveals the molecular mechanism by which crucian carp adapt to hypoxic stress. Our findings provide valuable references for promoting the healthy aquaculture of hypoxic-sensitive fish and breeding hypoxia-tolerant fish varieties.

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.

miR-330 targeting BCO2 is involved in carotenoid metabolism to regulate skin pigmentation in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

MicroRNAs (miRNAs) play a critical role in regulating skin pigmentation. As a key economic trait, skin color directly affects the market value of rainbow trout (Oncorhynchus mykiss), however, the regulatory mechanism of most miRNAs in fish skin color is still unclear.

RESULTS

In this study, the full-length cDNA sequence of β-carotene oxygenase 2 (BCO2, a key regulator of carotenoid metabolism) from the rainbow trout was obtained using rapid-amplification of cDNA ends (RACE) technology, and qRT-PCR was used to investigate the differential expression of miR-330 and BCO2 in 14 developmental stages and 13 tissues between wild-type rainbow trout (WTrt) and yellow mutant rainbow trout (YMrt). Additionally, the function of miR-330 was verified by overexpression and silencing in vitro and in vivo. The results showed that the complete cDNA sequence of BCO2 was 2057 bp with a 1707 bp ORF, encoding a 568 amino acid protein having a molecular weight of 64.07 kD. Sequence alignment revealed that higher conservation of BCO2 protein amongst fishes than amongst other vertebrates, which was further confirmed by phylogenetic analysis. The analysis of spatial and temporal expression patterns suggested that BCO2 and miR-330 were abundantly expressed from fertilized-stage to multi-cell as well as in the dorsal and ventral skin of WTrt and YMrt, and their expression patterns were opposite in most of the same periods and tissues. In vitro, luciferase reporter assay confirmed that BCO2 was a direct target of miR-330, and transfection of miR-330 mimics into rainbow trout liver cells resulted in a decrease in the expression of BCO2; conversely, miR-330 inhibitor had the opposite effect to the miR-330 mimics. In vivo, miR-330 agomir significantly decreased BCO2 expression in dorsal skin, tail fin, and liver. Furthermore, overexpression of miR-330 could suppress cell proliferation and induce apoptosis.

CONCLUSION

Our results showed that miR-330 is involved in the regulation of skin pigmentation in rainbow trout by targeting BCO2 and shows its promise as a potential molecular target to assist the selection of rainbow trout with better skin color patterns.

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.