Analysis of the differential expression profile of miRNAs in myocardial tissues of rats with burn injury

Figures & data

Validation of candidate miRNAs and correlation analysis

Figure 1. Effects of burn injury on the expression of C5b9, LDH, TnI, and CK-MB in rat myocardial tissues. (a). III burn rat model was validated by H&E staining of skin. (b). The expression levels of C5b9, LDH, TnI, and CK-MB in the burn injury group were significantly upregulated.*p-value < 0.05, **p-value < 0.01

Table 1. Differentially expressed miRNAs in the normal group compared with the burn injury group

miRNA	NC	MOD	Fold Change	log2(Fold Change)	Up-Down-Regulation	p-value
rno-miR-99a-3p	4.965	13.771	2.774	1.472	Up	0.014
rno-miR-582-5p	38.067	104.427	2.743	1.456	Up	0.011
rno-miR-551b-5p	5.259	10.967	2.085	1.06	Up	0.024
rno-miR-126b	10.85	3.06	0.282	−1.826	Down	0.034
rno-miR-344b-3p	1.84	5.358	2.912	1.542	Up	0.003
rno-miR-122-3p	10.996	22.97	2.089	1.063	Up	0.02
rno-miR-421-5p	5.341	11.733	2.197	1.135	Up	0.027
rno-miR-369-3p	480.761	3736.991	7.773	2.958	Up	0.001
rno-miR-181d-3p	3.774	9.695	2.569	1.361	Up	0.012
rno-miR-341	18.438	45.62	2.474	1.307	Up	0.046
rno-miR-154-3p	24	5.358	0.223	−2.163	Down	0.013
rno-miR-503-3p	34.748	176.664	5.084	2.346	Up	0.024
rno-miR-505-3p	0.92	4.332	4.709	2.235	Up	0.035
rno-miR-190b-3p	17.739	50.546	2.849	1.511	Up	0.003
rno-miR-450b-3p	17.478	42.315	2.421	1.276	Up	0.028
rno-miR-295-5p	43.715	16.816	0.385	−1.378	Down	0.032
rno-miR-667-5p	40.943	153.602	3.752	1.908	Up	0.004
rno-miR-1912-3p	13.692	6.37	0.465	−1.104	Down	0.036
rno-miR-19a-5p	4.432	1.783	0.402	−1.314	Down	0.027
rno-miR-3596b	2.864	12.229	4.269	2.094	Up	0.01
rno-miR-3551-5p	11.018	22.704	2.061	1.043	Up	0.032
rno-miR-879-5p	12.179	142.237	11.679	3.546	Up	0.001
rno-miR-9b-5p	4.286	1.527	0.356	−1.489	Down	0.021

Figure 2. miRNA profiles in myocardial tissues of the normal group and burn injury group. PCA and cluster analyses were used, and the MOD6 sample was removed. Differential expression analysis was then performed. The value of |log2(fold change)| ≥ 1 and p-value < 0.05 were used as standards. A total of 23 differentially expressed miRNA were screened, of which 6 were significantly downregulated, and 17 were significantly upregulated in the burn injury group. (a). PCA analysis. (b). Cluster analysis. c. Heatmap. d. Scatterplot

Figure 3. GO and KEGG enrichment analyses of differentially expressed miRNAs in the normal group compared with the burn injury group. Target mRNA of differentially expressed miRNAs were major enriched in BP category and CC category, in which metabolic process is the most important, and miRNA target genes are mainly enriched in PI3K/AKT, MAPK, Wnt, Rap1, and other signal pathways related to myocardial injury. (a). GO enrichment analysis. (b). KEGG enrichment analysis

Figure 4. miRNA–mRNA interaction network of differentially expressed miRNAs in the normal group compared with the burn injury group. The miRNA–mRNA interaction network was constructed by targeting genes in 15 signaling pathways mapped by KEGG, and only 18 differentially expressed miRNA were selected in the network regulatory map. Red indicated the up-regulated miRNAs green indicated the down-regulated miRNAs, pink indicated the target mRNA

Figure 5. Validation of candidate miRNAs and correlation analysis. qRT-PCR of differentially expressed miRNAs in the tissues of normal group compared with the tissue of burn injury group. Pearson’s analysis was used to identify the correlation between miRNAs and C5b9, LDH, TnI, and CK-MB. Finally, through ROC curve analysis, the sensitivity of miRNAs related to C5b9, TnI, and CK-MB was evaluated as indicators of burn injury. (a). qRT-PCR identified the differentially expressed miRNAs in the burn injury group compared with the normal group. (b). Pearson’s analysis. c. ROC curves

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.