The expression profile of the lncRNAs were analyzed by Deseq package (Affymetrix Inc., US). Samples were hybridized on the Human Clariom D (Thermo Fisher Scientific) gene chip. Background-adjustment, normalization, and log-transformation of signals intensity were performed with the Signal Space Transformation-Robust Multi-Array Average algorithm (RMA). Raw data were analyzed by the transcriptome analysis console (TAC) 4.0 software (Applied Biosystems, Foster City, CA, USA) awaiting further analysis [14 (link)]. The differentially expressed lncRNAs and mRNAs were screened according to the criteria of gene differential expression with |log2-fold change| (FC) more than 2 times and adjusted P<0.05. The differentially expressed lncRNAs were afterward clustered by a heatmap package while the hierarchical clustering diagram was drawn to show the results.