Effects of toxic Microcystis aeruginosa on the silver carp Hypophthalmichtys molitrix revealed by hepatic RNA-seq and miRNA-seq

High-throughput sequencing was applied to analyze the effects of toxic Microcystis aeruginosa on the silver carp Hypophthalmichthys molitrix. Silver carps were exposed to two cyanobacteria species (toxic and non-toxic) for RNA-seq and miRNA-seq analysis. RNA-seq revealed that the liver tissue contained 105,379 unigenes. Of these genes, 143 were significantly differentiated, 82 were markedly up-regulated, and 61 were remarkably down-regulated. GO term enrichment analysis indicated that 35 of the 154 enriched GO terms were significantly enriched. KEGG pathway enrichment analysis demonstrated that 17 of the 118 enriched KEGG pathways were significantly enriched. A considerable number of disease/immune-associated GO terms and significantly enriched KEGG pathways were also observed. The sequence length determined by miRNA-seq was mainly distributed in 20-23 bp and composed of 882,620 unique small RNAs, and 53% of these RNAs were annotated to miRNAs. As confirmed, 272 known miRNAs were differentially expressed, 453 novel miRNAs were predicted, 112 miRNAs were well matched with 7,623 target genes, and 203 novel miRNAs were matched with 15,453 target genes. qPCR also indicated that Steap4, Cyp7a1, CABZ01088134.1, and PPP1R3G were significantly differentially expressed and might play major roles in the toxic, detoxifying, and antitoxic mechanisms of microcystin in fish.