Fig. 2

The expression profiles of mRNAs and lncRNAs. a, c The expression profiles of mRNAs (upper panel) and lncRNAs (lower panel). The heat map shows the expression profile of mRNAs and lncRNAs. The top and left panel is the sample and gene tree constructed by Pearson Correlation; the value represents the log2 transformed values of (FPKM+ 1). b, d Dynamic changes in expression profile of mRNAs and lncRNAs. The value represents the pairwise Pearson Correlation. The correlation between every two samples was calculated by log2-transformed (FPKM+ 1) gene expression values. e The distributions of Shannon entropy-based temporal specificity scores that were calculated for distinct classes of mRNAs and lncRNAs. f Age-dependent clusters with significantly enriched GO terms and pathways map. Hierarchical clustering analysis grouped the age-dependent genes into 7 clusters. g Temporal profiles of age-dependent DEGs (upper panel) and DE lncRNAs (lower panel). Modules in different colors represent different temporal expression patterns which were significantly enriched by STEM analysis (non-significant modules were not shown). Numbers in the top indicates the module number and the lower numbers indicate numbers of DEGs and DE lncRNAs in each module