Novel ceRNA network construction associated with programmed cell death in acute rejection of heart allograft in mice

BACKGROUND: T cell-mediated acute rejection(AR) after heart transplantation(HT) ultimately results in graft failure and is a common indication for secondary transplantation. It’s a serious threat to heart transplant recipients. This study aimed to explore the novel lncRNA-miRNA-mRNA networks that contributed to AR in a mouse heart transplantation model. METHODS: The donor heart from Babl/C mice was transplanted to C57BL/6 mice with heterotopic implantation to the abdominal cavity. The control group was syngeneic heart transplantation with the same kind of mice donor. The whole-transcriptome sequencing was performed to obtain differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs) and lncRNAs (DElncRNAs) in mouse heart allograft. The biological functions of ceRNA networks was analyzed by GO and KEGG enrichment. Differentially expressed ceRNA involved in programmed cell death were further verified with qRT-PCR testing. RESULTS: Lots of DEmRNAs, DEmiRNAs and DElncRNAs were identified in acute rejection and control after heart transplantation, including up-regulated 4754 DEmRNAs, 1634 DElncRNAs, 182 DEmiRNAs, and down-regulated 4365 DEmRNAs, 1761 DElncRNAs, 132 DEmiRNAs. Based on the ceRNA theory, lncRNA-miRNA-mRNA regulatory networks were constructed in allograft acute rejection response. The functional enrichment analysis indicate that the down-regulated mRNAs are mainly involved in cardiac muscle cell contraction, potassium channel activity, etc. and the up-regulated mRNAs are mainly involved in T cell differentiation and mononuclear cell migration, etc. The KEGG pathway enrichment analysis showed that the down-regulated DEmRNAs were mainly enriched in adrenergic signaling, axon guidance, calcium signaling pathway, etc. The up-regulated DEmRNAs were enriched in the adhesion function, chemokine signaling pathway, apoptosis, etc. Four lncRNA-mediated ceRNA regulatory pathways, Pvt1/miR-30c-5p/Pdgfc, 1700071M16Rik/miR-145a-3p/Pdgfc, 1700071M16Rik/miR-145a-3p/Tox, 1700071M16Rik/miR-145a-3p/Themis2, were finally validated. In addition, increased expression of PVT1, 1700071M16Rik, Tox and Themis2 may be considered as potential diagnostic gene biomarkers in AR. CONCLUSION: We speculated that Pvt1/miR-30c-5p/Pdgfc, 1700071M16Rik/miR-145a-3p/Pdgfc, 1700071M16Rik/miR-145a-3p/Tox and 1700071M16Rik/miR-145a-3p/Themis2 interaction pairs may serve as potential biomarkers in AR after HT.