Abstract
This study aimed to detect changes in the complete transcriptome of MDCK cells after infection with the H5N1 and H3N2 canine influenza viruses using high-throughput sequencing, search for differentially expressed RNAs in the transcriptome of MDCK cells infected with H5N1 and H3N2 using comparative analysis, and explain the differences in the pathogenicity of H5N1 and H3N2 at the transcriptome level. Based on the results of our comparative analysis, significantly different levels of expression were found for 2,464 mRNAs, 16 miRNAs, 181 lncRNAs, and 262 circRNAs in the H3N2 infection group and 448 mRNAs, 12 miRNAs, 77 lncRNAs, and 189 circRNAs in the H5N1 infection group. Potential functions were predicted by performing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the target genes of miRNAs, lncRNAs and circRNAs, and the ncRNA-mRNA regulatory network was constructed based on differentially expressed RNAs. A greater number of pathways regulating immune metabolism were altered in the H3N2 infection group than in the H5N1 infection group, which may be one reason why the H3N2 virus is less pathogenic than is the H5N1 virus. This study provides detailed data on the production of ncRNAs during infection of MDCK cells by the canine influenza viruses H3N2 and H5N1, analyzed differences in the total transcriptomes between H3N2- and H5N1-infected MDCK cells, and explained these differences with regard to the pathogenicity of H3N2 and H5N1 at the transcriptional level.