Fig. 4

Functional distribution of the DEGs in candidate pathways in response to shading at time points of 4h and 8h. The enrichment score indicates intensiveness of DEGs (fold change ≥2) in a certain pathway (Enrichment score = (m/n): (M/N); m, the number of DEGs mapped to a certain pathway; n, the number of all DEGs annotated in transcriptome datasets; M, the number of unigenes mapped to a certain pathway; N, the number of all unigenes annotated in transcriptome datasets). A large enrichment score denotes a high degree of intensiveness. The p value (ranging from 0~1) was calculated using hypergeometric test through Bonferroni Correction and less p value means greater intensiveness. Gene number means number of DEGs mapped to a certain pathway according to KEGG database. S 4h and S 8h indicate time points at 4h and 8h throughout shading period