Integrated analysis of mRNA-miRNA expression in Tilapia infected with Tilapia lake virus (TiLV) and identifies primarily immuneresponse genes

We investigated differential gene expression in Tilapia infected with the Tilapia Lake virus (TiLV).We used highthroughput sequencing to identify mRNAs and miRNAs involved in TiLV infection progression We identified 25,359 differentially expressed genes that included 863 new genes. We identified 1770, 4142 and 4947 differently expressed genes comparing non-infected controls with 24 and 120 h infections and between the infected groups, respectively. These genes were enriched to 291 GO terms and 62 KEGG pathways and included immune system progress and virion genes. High-throughput miRNA sequencing identified 316 conserved miRNAs, 525 known miRNAs and 592 novel miRNAs. Furthermore, 138, 198 and 153 differently expressed miRNAs were found between the 3 groups listed above, respectively. Target prediction revealed numerous genes including erythropoietin isoform X2, double-stranded RNA-specific adenosine deaminase isoform Xl, bone morphogenetic protein 4 and tapasin-related protein that are involved in immune responsiveness. Moreover, these target genes overlapped with differentially expressed mRNAs obtained from RNA-seq. These target genes were significantly enriched to GO terms and KEGG pathways including immune system progress, virion and Wnt signaling pathways. Expression patterns of differentially expressed mRNA and miRNAs were validated in 20 mRNA and 19 miRNAs by qRT-PCR. We also were able to construct a miRNA-mRNA target network that can further understand the molecular mechanisms on the pathogenesis of TiLV and guide future research in developing effective agents and strategies to combat TiLV infections in Tilapia.