Total RNA was extracted from OM-MSCs using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. RNA integrity was analyzed with an Agilent RNA 6000 Chip in an Agilent 2100 Bioanalyser (Agilent Technologies). A microRNA library was then constructed according to the instruction manual, and the final library was qualified, quantified, and sequenced with paired-end reads on the HiSeq X-ten platform (BGI-Shenzhen, China).
Hierarchical cluster heatmaps were generated with pheatmap (v1.0.8) to show the gene expression levels in each sample. Differential expression analysis was performed using DESeq2 (v1.4.5). Differentially expressed genes between hypoxia and normoxia cultured groups were identified as those with |log2 (fold change) | > 1; the Q value was obtained after the P value was corrected using the false discovery rate, and a Q value ≤ 0.001 was considered statistically significant. All clean data have been submitted to the SRA database, and the SRA accession is PRJNA624725.
Hierarchical cluster heatmaps were generated with pheatmap (v1.0.8) to show the gene expression levels in each sample. Differential expression analysis was performed using DESeq2 (v1.4.5). Differentially expressed genes between hypoxia and normoxia cultured groups were identified as those with |log2 (fold change) | > 1; the Q value was obtained after the P value was corrected using the false discovery rate, and a Q value ≤ 0.001 was considered statistically significant. All clean data have been submitted to the SRA database, and the SRA accession is PRJNA624725.