RNA sequencing and short-read assembly were utilized to produce a transcriptome of livers from loaches (Misgurnus anguillicaudatus) fed with three different diets respectively containing fresh fish oil (FO group), medium oxidized fish oil (MO group) and high oxidized fish oil (HO group). A total of 60,663 unigenes were obtained in this study, with mean length 848.74 bp. 50,814, 49,584 and 49,814 unigenes were respectively obtained from FO, MO and HO groups. There were 2,343 differentially expressed genes between FO and MO, with 855 down- and 1,488 up-regulated genes in the MO group. 2,813 genes were differentially expressed between FO and HO, including 1,256 down- and 1,552 up-regulated genes in the HO group. 2,075 differentially expressed genes were found in the comparison of MO and HO, including 1,074 up- and 1,001 down-regulated genes in the MO group. Some differentially expressed genes, such as fatty acid transport protein (fatp), fatty acid binding protein (fabp), apolipoprotein (apo), peroxisome proliferator activated receptor-gamma (ppar-γ), acetyl-CoA synthetase (acs) and arachidonate 5-lipoxygenase (alox5), were involved in lipid metabolism, suggesting these genes in the loach were responsive to dietary oxidized fish oil. Results of transcriptome profilings here were validated using quantitative real time PCR in fourteen randomly selected unigenes. The present study provides insights into hepatic transcriptome profile of the loach, which is a valuable resource for studies of loach genomics. More importantly, this study identifies some important genes responsible for dietary oxidized fish oil, which will benefit researches of lipid metabolism in fish.

Elongation of very long-chain fatty acid (Elovl) proteins are key enzymes that catalyze the rate-limiting step in the fatty acid elongation pathway. The most recently discovered member of the Elovl family, Elovl8, has been proposed to be a fish-specific elongase with two gene paralogs described in teleosts. However, the biological functions of Elovl8 are still to be elucidated. In this study, we showed that in contrast to previous findings, elovl8 is not unique to teleosts, but displays a rather unique and ample phylogenetic distribution. For functional determination, we generated elovl8a (elovl8a^−/−) and elovl8b (elovl8b^−/−) zebrafish using CRISPR/Cas9 technology. Fatty acid composition in vivo and zebrafish liver cell experiments suggest that the substrate preference of Elovl8 overlapped with other existing Elovl enzymes. Zebrafish Elovl8a could elongate the polyunsaturated fatty acids (PUFAs) C18:2n-6 and C18:3n-3 to C20:2n-6 and C20:3n-3, respectively. Along with PUFA, zebrafish Elovl8b also showed the capacity to elongate C18:0 and C20:1. Gene expression quantification suggests that Elovl8a and Elovl8b may play a potentially important role in fatty acid biosynthesis. Overall, our results provide novel insights into the function of Elovl8a and Elovl8b, representing additional fatty acid elongases not previously described in chordates.

Dojo loach, Misgurnus anguillicaudatus is a freshwater fish species of the loach family Cobitidae, using its posterior intestine as an accessory air-breathing organ. Little is known about the molecular regulatory mechanisms in the formation of intestinal air-breathing function of M. anguillicaudatus. Here high-throughput sequencing of mRNAs was performed from six developmental stages of posterior intestine of M. anguillicaudatus: 4-Dph (days post hatch) group, 8-Dph group, 12-Dph group, 20-Dph group, 40-Dph group and Oyd (one-year-old) group. These six libraries were assembled into 81300 unigenes. Totally 40757 unigenes were annotated. Subsequently, 35291 differentially expressed genes (DEGs) were scanned among different developmental stages and clustered into 20 gene expression profiles. Finally, 15 key pathways and 25 key genes were mined, providing potential targets for candidate gene selection involved in formation of intestinal air-breathing function in M. anguillicaudatus. This is the first report of developmental transcriptome of posterior intestine in M. anguillicaudatus, offering a substantial contribution to the sequence resources for this species and providing a deep insight into the formation mechanism of its intestinal air-breathing function. This report demonstrates that M. anguillicaudatus is a good model for studies to identify and characterize the molecular basis of accessory air-breathing organ development in fish.

Elongation of very long-chain fatty acids protein 6 (Elovl6) has been reported to be associated with clinical treatments of a variety of metabolic diseases. However, there is no systematic and comprehensive study to reveal the regulatory role of Elovl6 in mRNA, protein and phosphorylation levels. We established the first knock-out (KO), elovl6^−/−, in zebrafish. Compared with wild type (WT) zebrafish, KO presented significant higher whole-body lipid content and lower content of fasting blood glucose. We utilized RNA-Seq, tandem mass tag (TMT) labeling-based quantitative technology and liquid chromatography-tandem mass spectrometry (LC-MS/MS) to perform the transcriptomic, proteomic and phosphoproteomic analyses of livers from WT and elovl6^−/− zebrafish. There were 734 differentially expressed genes (DEG) and 559 differentially expressed proteins (DEP) between elovl6^−/− and WT zebrafish, identified out of quantifiable 47251 transcripts and 5525 proteins. Meanwhile, 680 differentially expressed phosphoproteins (DEPP) with 1054 sites were found out of quantifiable 1230 proteins with 3604 sites. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) analysis of the transcriptomic and proteomic data further suggested that the abnormal lipid metabolism and glucose metabolism in KO were mainly related to fatty acid degradation and biosynthesis, glycolysis/gluconeogenesis and PPAR signaling pathway. Based on phosphoproteomic analyses, some kinases critical for lipid metabolism and glucose metabolism, including ribosomal protein S6 kinase (Rps6kb), mitogen-activated protein kinase14 (Mapk14) and V-akt murine thymoma viral oncogene homolog 2-like (Akt2l), were identified. These results allowed us to catch on the regulatory networks of elovl6 on lipid and glucose metabolism in zebrafish. To our knowledge, this is the first multi-omic study of zebrafish lacking elovl6, which provides strong datasets to better understand many lipid/glucose metabolic risks posed to human health.