Where R is the number of reads in a sample assigned to a gene, N is the total number of mapped reads in the sample, and L is the length of the gene [16 (link)].
Pairwise analyses of differential gene expression were carried out with biological replicates using DEGseq version 1.16. The gene expression levels were determined using a negative binomial distribution model, and genes with q ≤ 0.05 and |log2ratio| ≥ 1 were identified as differentially expressed genes (DEGs) [17 (link)]. The expression levels of all divergent genes were log2-transformed, the Euclidean distance was calculated, and clustering was performed using Hierarchical Cluster methods in R software (version 3.2.5). Quantitative reverse transcription PCR (qRT-PCR) was performed using a Power SYBRGreen RT-PCR Reagents Kit (Applied Biosystems, Carlsbad, CA, USA) in triplicate for each sample and reference gene. Actin was used as a reference gene. Quantification of all chosen lncRNA and mRNA expressions was conducted using the comparative CT ratio between EXOP, StAR with Actin(the primers data showed in