Figure 4. Enrichment analysis of Gene Ontology (GO) terms and
KEGG pathways
(a ) The number of significantly enriched GO terms obtained
against the top 500 DEG sequences. The treatment contrasts are
represented with similar colors. The GO terms are categorized into
Biological process (BP), Molecular function (MF) and Cellular process
(CC). (b ) Number of enriched GO terms obtained from direct
count of sequences in Biological process (BP) category. (c )The
KEGG metabolic pathways significantly enriched in red vs control and
(d ) blue vs control. corresponding circle sizes represent the
counts of sequences with p-value <0.01.
Figure 5. DEGs from light signaling and circadian
rhythm. Schematic representation of red and blue light signal
perception by higher plants which are likely to be involved in light
regulated anthocyanin biosynthesis. The associated DEGs from red and
blue light vs control contrasts were represented in two color code boxes
based on log2 fold changes. Red light treatment is shown left and blue
light treatment on the right side of the box.Abbreviations and gene definitions: PHYA/B, Phytochrome;
CRY1/2, Cryptochrome; COP1, constitutive photomorphogenic; Elf3, early
flowering 3; FT, flowering locus T; PRR5, pseudo response regulator 5;
CK2α, casein kinase II subunit alpha; Cab, chlorophyll a/b binding
protein; Exp, expansin; HY5, elongated hypocotyl; FKF1- flavin-binding
kelch domain F box protein; SPA, suppressor of PhyA; HY5, elongated
hypocotyl 5; CDF1, cycling DOF factor 1; Co, Constans.
Figure 6. DEGs from anthocyanin, carotenoid and ABA
biosynthesis. Schematic representation of bilberry anthocyanin
biosynthetic pathway branching from phenylpropanoid biosynthesis (top)
and representation of carotenoid biosynthetic pathway leading towards
abscisic acid (ABA) biosynthesis and catabolism. DEGs from flavonoid,
carotenoid-ABA pathway genes and selected TFs visualized as heatmap
based on log2 fold changes obtained from light treatments against the
control samples. Enzyme abbreviations: PAL,
phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase; 4CL,
4-coumarate:CoA ligase; CCoAOMT, caffeoyl-CoA O-methyltransferase;
Prx,cPrx- peroxidases, cationic peroxidases; CHS, chalcone synthase;
CHI, chalcone isomerase; F3H, flavanone 3-hydroxylase; F3’H, flavonoid
3’ hydroxylase; F3’5’H, flavonoid 3’5’ hydroxylase; DFR, dihydroflavonol
4-reductase; ANS, anthocyanidin synthase; UFGT, UDP-glucose flavonoid
3-O -glucosyltransferase; PSY, phytoene synthase; ZDS,
zeta-carotene desaturase; BCH, beta-carotene hydroxylase; CYP 450-BCH,
carotenoid β -ring hydroxylase of cytochrome P450 family; LUT1,
lutein deficient 1; NCED, 9-cis -epoxycarotenoid dioxygenase; ABA
8’ hyd, Abscisic acid 8’ hydroxylase, PYR/PYL, Pyrabactin-resistance
like.
Figure 7. DEGs from SNARE mediated vesicular trafficking.Schematic representation of vesicular transport of anthocyanins mediated
by SNARE proteins interaction through transmembrane - Endoplasmic
reticulum (ER) - golgi network. The associated transporter genes from
red and blue light vs control contrasts DEGs were represented in color
code boxes based on log2 fold changes. Red light treatment is shown left
and blue light treatment on the right side of the box.Abbreviations and gene definitions: Stx, syntaxin-like;
Bos1, Gos, Qb type golgi SNAP receptor complex; Ykt6, Sec22, vesicle
transporter/VAMP like protein; SNAP, soluble NSF attachment protein;
SNARE, ‘SNAP REceptor’.