Background: Highly pathogenic avian influenza virus (HPAIV) is considered a global threat to both human health and the poultry industry. MicroRNAs (miRNA) can modulate the immune system by affecting gene expression patterns in HPAIV-infected chickens. Objectives: To gain further insights into the role of miRNAs in immune responses against H5N1 infection, as well as the development of strategies for breeding disease-resistant chickens, we characterized miRNA expression patterns in tracheal tissues from H5N1-infected Ri chickens. Methods: miRNAs expression was analyzed from two H5N1-infected Ri chicken lines using small RNA sequencing. The target genes of differentially expressed (DE) miRNAs were predicted using miRDB. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were then conducted. Furthermore, using quantitative real-time polymerase chain reaction, we validated the expression levels of DE miRNAs (miR-22-3p, miR-146b-3p, miR-27b-3p, miR-128-3p, miR-2188-5p, miR-451, miR-205a, miR-203a, miR-21-3p, and miR-200a-3p) from all comparisons and their immune-related target genes. Results: A total of 53 miRNAs were significantly expressed in the infection samples of the resistant compared to the susceptible line. Network analyses between the DE miRNAs and target genes revealed that DE miRNAs may regulate the expression of target genes involved in the transforming growth factor-beta, mitogen-activated protein kinase, and Toll-like receptor signaling pathways, all of which are related to influenza A virus progression. Conclusions Collectively, our results provided novel insights into the miRNA expression patterns of tracheal tissues from H5N1-infected Ri chickens. More importantly, our findings offer insights into the relationship between miRNA and immune-related target genes and the role of miRNA in HPAIV infections in chickens.

Background: Highly pathogenic avian influenza viruses (HPAIVs) is an extremely contagious and high mortality rates in chickens resulting in substantial economic impact on the poultry sector. Therefore, it is necessary to elucidate the pathogenic mechanism of HPAIV for infection control. Objective: Gene set enrichment analysis (GSEA) can effectively avoid the limitations of subjective screening for differential gene expression. Therefore, we performed GSEA to compare HPAI-infected resistant and susceptible Ri chicken lines. Methods: The Ri chickens Mx(A)/BF2(B21) were chosen as resistant, and the chickens Mx(G)/ BF2(B13) were selected as susceptible by genotyping the Mx and BF2 genes. The tracheal tissues of HPAIV H5N1 infected chickens were collected for RNA sequencing followed by GSEA analysis to define gene subsets to elucidate the sequencing results. Results: We identified four differentially expressed pathways, which were immune-related pathways with a total of 78 genes. The expression levels of cytokines (IL-1β, IL-6, IL-12), chemokines (CCL4 and CCL5), type interferons and their receptors (IFN-β, IFNAR1, IFNAR2, and IFNGR1), Jak-STAT signaling pathway genes (STAT1, STAT2, and JAK1), MHC class I and II and their co-stimulatory molecules (CD80, CD86, CD40, DMB2, BLB2, and B2M), and interferon stimulated genes (EIF2AK2 and EIF2AK1) in resistant chickens were higher than those in susceptible chickens. Conclusions Resistant Ri chickens exhibit a stronger antiviral response to HPAIV H5N1 compared with susceptible chickens. Our findings provide insights into the immune responses of genetically disparate chickens against HPAIV.