MiRNA-seq analysis of spleen and head kidney tissue from aquacultured largemouth bass (Micropterus salmoides) in response to Aeromonas hydrophila infection.

Recently, the same fish diseases, which have been found in pond farming, have been found in the newly tested largemouth bass (Micropterus salmoides) system. Bacterial septicemia caused by Aeromonas hydrophila occurs frequently in largemouth bass culture leading to significant economic losses. To investigate the role miRNA in the largemouth bass disease resistance, twelve (2 tissues (spleen and head kidney) × 2 experimental groups (infected and control) × three biological replicates) small RNA libraries were constructed and sequenced with miRNA-seq. A total of 26 differentially expressed miRNAs, 8 upregulated and 18 downregulated, were identified in the spleen, and 19 differentially expressed miRNAs, 9 upregulated and 10 downregulated, were identified in head kidney (fold change ≥ 2 or ≤ 0.5 and P ≤ 0.05). The differentially expressed miRNAs with the largest fold change were selected for target gene prediction using GO and KEGG analysis. Six miRNAs in the spleen and 5 miRNAs in the head kidney were selected. Analysis showed that, of all the immune and metabolic pathways, the FoxO signaling pathway was enriched in both the spleen and head kidney groups. Common target genes of the pathway included AMP-activated catalytic subunit alpha 1 (prkaa1), phosphatidylinositol 3-kinase (pik3r3b), serine/threonine-protein kinase (plk2), and forkhead box protein G1 (foxg1a). MiRNAs (such as miR-126-5P, miR-126-3P) are involved in immune response and cell cycle functions as they regulate targeted genes in the FoxO pathway. These results will enhance our understanding of the molecular mechanisms underlying immune responses to bacterial septicemia and facilitate molecular-assisted selection of resistant strains of largemouth bass.

Differential transcriptomic analysis of crayfish (Procambarus clarkii) from a rice coculture system challenged by Vibrio parahaemolyticus.

Rice-crayfish (Procambarus clarkii) coculture is an effective farming mode and has been promoted in various regions of China. However, infection in crayfish can be a significant economic drain. We found crayfish infected with Vibrio parahemolyticus (VP), and to understand the molecular mechanisms of the immune responses of crayfish to VP infection, Illumina sequencing was employed to identify changes in the mRNA of hepatopancreatic tissue. A total of 47.30 and 43.01million high-quality transcriptome reads were generated from the hepatopancreatic samples of the experimental group (EG) and control group (CG), respectively. We found 5559 genes were significantly differentially expressed, including 2521 up-regulated genes (45.35%) and 3038 down-regulated genes (54.65%). These genes were enriched in 126 GO terms and 76 KEGG pathways (P ≤ 0.05), including the MAPK and PI3K-Akt signaling pathways and cell adhesion molecules, with 23 up-regulated genes and 3 down-regulated genes related to immune responses in the EG relative to the CG. Histopathological analysis revealed that the epithelial cells of the hepatopancreatic tubules in the EG were severely atrophic, necrotic, and exfoliated, resulting in thin and collapsing hepatopancreatic tubules. The expression patterns of 8 differentially expressed genes involved in immune responses were validated by quantitative real-time RT-PCR. These results provide a valuable basis for the immune responses of crayfish to acute hepatopancreatic necrosis disease at transcriptome level.

Integrated analysis of RNA-seq and microRNA-seq depicts miRNA-mRNA networks involved in stripe patterns of Botia superciliaris skin.

Botia superciliaris, an endemic cobitid fish in China, is widely accepted by Chinese consumers because its edibility. Recently, the black and yellow stripes of B. superciliaris skin have made this species increasingly popular as a novel ornamental fish. However, the genetic basis of the stripe patterns in B. superciliaris skin has not been extensively studied. In this study, Illumina sequencing was employed to identify the mRNAs and miRNAs involved in stripe pattern formation in B. superciliaris skin. A total of 147.25 and 155.15 million (M) high-quality transcriptome reads were generated from three black and yellow skin libraries respectively, which resulted in 159,327 unigenes that were used as reference sequences. A total of 3169 genes exhibited significantly differential expression patterns (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05), including 1891 upregulated genes (59.67%) and 1278 downregulated genes (40.33%) in black vs yellow skin. These genes were enriched in 50 GO terms and 10 KEGG pathways (q ≤ 0.05), including melanogenesis, with 21 upregulated genes and 5 downregulated genes in black vs yellow skin. Based on the zebrafish genome, miRNA-seq identified a total of 355 miRNAs, which included 38 novel miRNAs. Furthermore, 87 differentially expressed miRNAs including 50 upregulated and 37 downregulated miRNAs were identified in different color skin (fold-change ≥ 2 or ≤ 0.5 and q ≤ 0.05). Then, target prediction revealed a variety of putative target genes; differentially expressed mRNAs and miRNAs patterns of high-throughput sequencing were validated in 5 mRNAs and miR-217-5p by qRT-PCR. In vivo tests and dual-luciferase reporter assay revealed that overexpression of miR-217-5p can inhibit pheomelanin formation by targeting Zgc. In this study, a comparative analysis was conducted to profile the transcriptome of black and yellow skin for B. superciliaris, and we detected key genes and important miRNAs involved in the B. superciliaris skin pigmentation process. These results will enhance understanding of molecular mechanisms underlying skin pigmentation and facilitate molecular-assisted selection of highly valued skin colors.