The aim of our scientific work was to investigate the aberrant expression of miRNAs that correlate with recurrent and invasive nature of NFPAs. We assessed 84 cancer–specific miRNAs in fresh tumor tissues from Bulgarian patients. Additionally, we wanted to analyze and select a specific miRNA profile that could be used as a prognostic recurrence and size-related biomarker.
Invasion-related Mirnas
Expression analysis showed deregulated miRNAs when comparing invasive NFPA (n = 8) and non-invasive NFPA (n = 12) (Table 2). The most significantly down-regulated were miR-210-3p (p < 0.001), miR-149 (p < 0.05) and miR-29b (p < 0.001), miR-103а (p < 0.001), and miR-107 (p < 0.001) in invasive NFPA with fold regulation (FR) = -6.12; FR = -8.4; FR = -4.03, FR = -2.71, FR = -2.71, respectively. We found miR-100 slightly up-regulated, FR = 3.43 (Fig. 1).
Table 2
Deregulated miRNA in invasive NFPA compared to non-invasive NFPA
miRNA name | FR | p-value |
miR-149-3p | -8.4 | 0.01 |
miR-210-3p | -6.12 | (p < 0.001) |
miR-29b-3p | -4.03 | (p < 0.001) |
miR-103a-3p | -2.71 | (p < 0.001) |
miR-107 | -2.71 | (p < 0.001) |
miR-100 | 3.43 | 0.04 |
FR - fold regulation |
Additionally, we observed aberrant expression of some short-noncoding RNAs, such as U6, SNORD49A, SNORD38B that were included in the panel as potential reference controls.
Recurrence-related Mirnas
Furthermore, we divided invasive NFPA into two groups. Recurrent tumors were compared to newly diagnosed tumors (Table 3). Significantly up-regulated were miR-20а, miR-19a, miR-16, miR-17, while miR-106a was slightly up-regulated (Fig. 2).
Table 3
Deregulated miRs in recurrent NFPAs compared to NFPAs without recurrence
miR name | FR | p-value |
miR-20а | 2.35 | 0.005 |
miR-16 | 2.04 | 0.009 |
miR-19a | 2.07 | 0.007 |
miR-17 | 2.12 | 0.012 |
miR-106a | 1.88 | 0.02 |
FR - fold regulation |
Size-related Mirnas
Additionally, we compared miRNA expression in giant NFPAs (n = 5) and other macroadenomas (n = 15). We found that miR-186 (FR = -24.08, p = 0.001), miR-17 (FR = -24.88, p = 0.006) and miR-210 (FR = -22.77, p = 0.012) were the most significantly down-regulated, while miR-9 (FR = 51.94) (p = 0.0015) was the most highly expressed and up-regulated. MiR-9 is a brain specific miRNA that showed significantly increased expression in our patients with giant NFPAs (Table 4).
Table 4
Deregulated miRs in giant NFPAs compared to other macroadenomas
miRNA name | FR | p-value |
Up-regulated |
miR-9 | 51.94 | 0.0015 |
Down-regulated |
miR-17 | -24.88 | 0.006 |
miR-186 | -24.08 | 0.001 |
miR-210-3p | -22.77 | 0.012 |
FR - fold regulation |
The correaltion anaylsis between clinocopathological characteristics and miRNA expession showed a possitive correlation between miR-17, miR-210, miR-186, miR-29b, miR-20, miR-149, miR-106a (r > 0.7, p > 0.01) and invasiveness of NFPA.
Roc Curve Analysis
ROC (Receiver Оperating Characteristic) curves and AUC (Area Under the Curve) were used to calculate the biomarker power of the selected miRNAs and their diagnostic potential.
ROC analysis of deregulated miRNA showed that miR-210-3p (AUC = 0.903, p = 0.003, 95% Confidence Interval (CI) = 0.721-1), miR-149 (AUC = 1, p < 0.001, 95% CI = 1), and miR-29b (AUC = 0.898, p = 0.004, 95% CI = 1) could be used as diagnostic biomarkers to differentiate between invasive NFPA and non-invasive NFPA (Fig. 3).
ROC analysis of miRNA potential to differentiate recurrent NFPA showed that the combination of miR-17, miR-20 and miR-106a is a suitable biomarker (AUC = 1, p = 0.025, 95% CI = 1) (Fig. 4).
Target Analysis
The cell-signaling pathways regulated by some of the most disturbed miRNAs were analyzed by GeneGlobe tool. Some of the target genes of our disturbed miRNAs were important tumor-suppressors and cell-signaling regulators (Table 5).
Table 5
Selected putative miRNA targets of our selected deregulated miRNAs
miRNA | Some of the putative mRNA targets | Function of the targets |
miR-186 | PIK3R3, BCL2L11, DICER1, MAPK1, VEGFA, CDK5R1, WNT5A, PIK3- AKT | angiogenesis, cell proliferation and invasion |
miR-20а | TP53INP1, RB1, CASP2, CASP7 | tumor-suppressors, apoptosis |
miR-17 | TP53INP1, RB1, CASP7, VEGFA | tumor-suppressors, angiogenesis, apoptosis |
miR-16 | DICER1, TSC1, BCL2, FGFR1, NF1, TSC1 | oncogenes, tumour-suppressors, gene expression |
miR-106а | RB1, CASP2, CASP7 | tumor-suppressors, apoptosis |
miR-210-3p | ARMC1 BDNF B4GALT5 EFNA3, ISCU MID1IP1, NEUROD2, NDUFA4, SCN9A, FGFR | neural development and erythropoiesis, apoptosis, migration and cell division |
miR-149-3p | MAPK, TGF | cell growth, proliferation |
miR-29b-3p | TET3, COL1A1, DNMTs, | DNA methylation, cell connection |
miR-9 | PTCH1, HOXA11, MAP3K1,PIK3C2A, RAB11FIP4 SCN2B, RABs, EFNA1, IGF2BP2, KLF13 | apoptosis, DNA reparation, development, cell growth, intracellular vesicular trafficking |