“…NextSeq2000 was used to generate the sequencing data of 36 million clusters (or 72 million F + R reads) per sample with the P3 flow cell at Marshal University Genomics Core (Huntington, WV, USA). RNA-Seq data analysis follows our previous work [ 34 , 35 ]: briefly, RNA-Seq read aligned to human genome by subread [ 36 ], read counts summarized based on RefSeq gene annotation by featurecount [ 37 ], expression level quantification by RPKM [ 38 ] with an in-house script, visualization of gene expression by MeV [ 39 ], prediction of differentially expressed (DE) genes by EdgeR (FDR = 0.05; |log2FC| > 1, and log2(Count Per million) > 0), and gene set enrichment analysis against hallmark gene sets by GSEA [ 40 , 41 ].…”