Comparative analysis between 59 pairs of matched normal and papillary thyroid cancer samples identified 455 differentially expressed lncRNAs (log2 fold change > 1 or < −1; adj. p-value < 1 × 10−7; Fig. 1A), being 71 lncRNAs upregulated and 384 lncRNAs downregulated in PTC (Supplemental Table S1). The same samples presented a total of 2016 mRNAs (log2 fold change > 1 or <−1; adj. p-value < 0.05; Supplemental Table S2) and 186 microRNAs (log2 fold change > 1 or < −1; adj. p-value < 0.05; Supplemental Table S3) differentially expressed.Figure 1Differentially expressed lncRNAs between N × T and WT × BRAFV600E were identified. (A) Volcano plot of DE lncRNAs between N × T (log2 fold change > 1 or < −1; adj. p-value < 1 × 10−7). (B) Volcano plot of DE lncRNAs between WT × BRAFV600E (log2 fold change > 1 or < −1; adj. p-value < 1 × 10−4). (C) Venn Diagram of common DE lncRNAs between N × T and WT × BRAFV600E. (D) Heatmap* of DE lncRNAs between N × T (log2 fold change > 3 or < −3; adj. p-value < 1 × 10−7). (E) Heatmap* of DE lncRNAs between WT × BRAFV600E (log2 fold change > 2.5 or < −2.5; adj. p-value < 1 × 10−4). *For hierarchical clustering, one minus Spearman rank correlation was performed.