Integrated mRNA–miRNA transcriptome profiling of blood immune responses potentially related to pulmonary fibrosis in forest musk deer

Abstract

BackgroundForest musk deer (FMD, Moschus Berezovskii) is a critically endangered species world-widely, the death of which can be caused by pulmonary disease in the farm. Pulmonary fibrosis (PF) was a huge threat to the health and survival of captive FMD. MicroRNAs (miRNAs) and messenger RNAs (mRNAs) have been involved in the regulation of immune genes and disease development. However, the regulatory profiles of mRNAs and miRNAs involved in immune regulation of FMD are unclear.MethodsIn this study, mRNA-seq and miRNA-seq in blood were performed to constructed coexpression regulatory networks between PF and healthy groups of FMD. The hub immune- and apoptosis-related genes in the PF blood of FMD were explored through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Further, protein–protein interaction (PPI) network of immune-associated and apoptosis-associated key signaling pathways were constructed based on mRNA-miRNA in the PF blood of the FMD. Immune hub DEGs and immune hub DEmiRNAs were selected for experimental verification using RT-qPCR.ResultsA total of 2744 differentially expressed genes (DEGs) and 356 differentially expressed miRNAs (DEmiRNAs) were identified in the PF blood group compared to the healthy blood group. Among them, 42 DEmiRNAs were negatively correlated with 20 immune DEGs from a total of 57 correlations. The DEGs were significantly associated with pathways related to CD molecules, immune disease, immune system, cytokine receptors, T cell receptor signaling pathway, Th1 and Th2 cell differentiation, cytokine-cytokine receptor interaction, intestinal immune network for IgA production, and NOD-like receptor signaling pathway. There were 240 immune-related DEGs, in which 186 immune-related DEGs were up-regulated and 54 immune-related DEGs were down-regulated. In the protein-protein interaction (PPI) analysis of immune-related signaling pathway, TYK2, TLR2, TLR4, IL18, CSF1, CXCL13, LCK, ITGB2, PIK3CB, HCK, CD40, CD86, CCL3, CCR7, IL2RA, TLR3, and IL4R were identified as the hub immune genes. The mRNA-miRNA coregulation analysis showed that let-7d, miR-324-3p, miR-760, miR-185, miR-149, miR-149-5p, and miR-1842-5p are key miRNAs that target DEGs involved in immune disease, immune system and immunoregulation.ConclusionThe development and occurrence of PF were significantly influenced by the immune-related and apoptosis-related genes present in PF blood. mRNAs and miRNAs associated with the development and occurrence of PF in the FMD.