For differential expression analysis of mRNA and microRNA was used edgeR package52 through TCGAbiolinks53. To identify differentially expressed lncRNAs between groups, it was used the paired/unpaired Student t test to assess the statistical difference of mean expression values between the two groups50. LncRNA with median value equal zero were excluded and fold change was calculated using median expression values. In all differential analysis, p-values were adjusted for False Discovery Rate (FDR) < 0.05 as multiple hypothesis test correction method.