Co-expression analysis and prediction of lncRNA function

The Pearson correlation coefficient (r) of each lncRNA and mRNA were calculated, and the correlation tests were performed. Differential lncRNA and mRNA relationship pairs were screened; r>0.85 and the adjusted P<0.05 following correction of the false discovery rate (FDR) was regarded as differential. The function of each lncRNA was analyzed by the relationship pairs in the regulatory lncRNA-mRNA network. The clusterProfiler package in R was used to analyze the pathways and BPs in which the target genes of these lncRNAs were involved. Pathways and BPs with the number of parameter enrichment genes count ≥2 and P<0.05 were regarded as significant enrichment results.