Construction and Analysis of ceRNA Networks Reveal the Key Genes Associated with Bovine Herpesvirus Type 1 Infection

Abstract

Background Virus infection can cause the changes of lncRNA expression levels to regulate the interaction between virus and host, but the relationship between BHV-1 infection and lncRNA has not been reported. Methods In this study, in order to reveal the molecular mechanism of RNA in BoHV-1 infection, the Madin-Darby bovine kidney (MDBK) cells were infected with BoHV-1, transcriptome sequencing were performed by next-generation sequencing at 18 h or 24 h or 33 h of viral infection and then based on the competitive endogenous RNA (ceRNA) theory, lncRNA-miRNA-mRNA networks were constructed using these high-throughput sequencing data. The network analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for functional annotation and exploration of ncRNA ceRNAs in BoHV-1 infection. Results The results showed that 48 lncRNAs, 123 mRNAs and 20 miRNAs as differentially expressed genes, and the mitogen activated protein kinase (MAPK) pathway and calcium signaling pathway were significantly enriched in the ceRNA network. Some differentially expressed lncRNA genes were randomly selected for verification by RT-qPCR, and the results showed that their expression trend was consistent with the results of transcriptome sequencing data. Conclusion This study revealed that BoHV-1 infection can affect the expression of RNAs in MDBK cells and the regulation of ceRNA network to carry out corresponding biological functions in the host, but further experimental studies are still necessary to prove the hub genes function in ceRNA network and the molecular mechanism in BoHV-1 infection.