
Figure 1
Heat map representation of differentially expressed genes belonging to the to “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. Arbitrary signal intensity acquired from microarray analysis is represented by colours (green, higher; red, lower expression). Log2 signal intensity values for any single gene were resized to Row Z-Score scale (from -2, the lowest expression to +2, the highest expression for single gene)
Table 1
Gene symbols, fold change in expression ratio, Entrez gene IDs, corrected p values and mean value of fold change ratio of studied genes
| GENE SYMBOL | RATIO D7/D1 | RATIO D15/D1 | RATIO D30/D1 | ADJUSTED P VALUE D7/D1 | ADJUSTED P VALUE D15/D1 | ADJUSTED P VALUE D30/D1 | ENTREZ GENE ID | MEAN RATIO |
|---|---|---|---|---|---|---|---|---|
| ACSL4 | 2,078475 | -1,049840 | -1,389865 | 0,000888 | 0,778847 | 0,027679 | 448980 | -0,120410 |
| HAAO | 2,260462 | 1,012229 | -1,484548 | 0,000495 | 0,948252 | 0,010911 | 100514283 | 0,596048 |
| PPARA | 2,156741 | 1,463697 | 1,431760 | 0,001697 | 0,036412 | 0,040484 | 397239 | 1,684066 |
| CPOX | 2,562460 | 1,513475 | 1,395239 | 0,000490 | 0,019946 | 0,042725 | 100511474 | 1,823725 |
| ALDOC | 3,100575 | 1,215495 | 1,317360 | 0,000301 | 0,274540 | 0,110445 | 100512013 | 1,877810 |
| PGM1 | 2,424519 | 1,680756 | 1,691567 | 0,000276 | 0,002924 | 0,002153 | 397566 | 1,932281 |
| TPI1 | 2,751177 | 1,737387 | 1,421065 | 0,000059 | 0,000779 | 0,006133 | 100157582 | 1,969876 |
| ME1 | 2,691975 | 1,528067 | 1,911722 | 0,000358 | 0,017206 | 0,001723 | 397538 | 2,043921 |
| ACSL3 | 2,049555 | 1,676912 | 2,729973 | 0,002097 | 0,008605 | 0,000218 | 100233169 | 2,152147 |
| FECH | 2,109918 | 2,722032 | 2,420194 | 0,000678 | 0,000117 | 0,000152 | 100322873 | 2,417382 |
| SLC25A1 | 3,272247 | 1,857686 | 2,239022 | 0,000043 | 0,000649 | 0,000117 | 100154310 | 2,456318 |
| OXSM | 2,234415 | 2,697009 | 2,514316 | 0,000560 | 0,000143 | 0,000144 | --- | 2,481914 |
| PFKFB4 | 2,978245 | 3,699473 | 2,573760 | 0,000132 | 0,000034 | 0,000122 | 100158056 | 3,083826 |
| PDK1 | 4,856476 | 1,981458 | 2,658224 | 0,000100 | 0,004418 | 0,000537 | 100286871 | 3,165386 |
| PCBD1 | 2,831266 | 3,234674 | 3,436204 | 0,000047 | 0,000015 | 0,000007 | 100155605 | 3,167382 |
| ACACA | 3,411749 | 3,730045 | 3,501388 | 0,000588 | 0,000307 | 0,000279 | 397324 | 3,547727 |
| TSPO | 3,123906 | 3,962945 | 3,669010 | 0,000048 | 0,000012 | 0,000009 | 396592 | 3,585287 |
| ACOT13 | 2,203769 | 4,264528 | 4,968930 | 0,000233 | 0,000008 | 0,000003 | 100514587 | 3,812409 |
| SCD | 5,037882 | 1,402767 | 5,796847 | 0,000092 | 0,087483 | 0,000025 | 396670 | 4,079166 |
| PDK3 | 2,773770 | 4,066151 | 5,486628 | 0,000182 | 0,000024 | 0,000006 | 100153858 | 4,108850 |
| AASS | 2,332931 | 4,119242 | 6,132447 | 0,007329 | 0,000476 | 0,000092 | 100513962 | 4,194873 |
| ACSS2 | 6,349718 | 2,503393 | 5,620930 | 0,000007 | 0,000095 | 0,000002 | 100153866 | 4,824680 |
| PANK1 | 7,855562 | 3,491674 | 7,158650 | 0,000007 | 0,000033 | 0,000002 | 100154650 | 6,168629 |
| P2RX7 | 6,165142 | 8,522042 | 5,795398 | 0,000027 | 0,000007 | 0,000010 | 497623 | 6,827528 |
| ALDH1L2 | 4,535793 | 8,626347 | 22,102755 | 0,000038 | 0,000004 | 0,000001 | 100151976 | 11,754965 |

Figure 2
The circle plot showing the differently expressed genes and z-score of the “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. The outer circle shows a scatter plot for each term of the fold change of the assigned genes. Green circles display up- regulation and red ones down- regulation. The inner circle shows the z-score of each GO BP term. The width of the each bar corresponds to the number of genes within GO BP term and the color corresponds to the z-score

Figure 3
The representation of the mutual relationship of differently expressed genes that belong to “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. The ribbons indicate which gene belongs to which categories. The middle circle represents logarithm from fold change (LogFC) between D7/D1, D15/D1 and D30/D1 respectively. The color of each block corresponds to the LogFC of each gene (green – up-regulated, red – down-regulated). The genes were sorted by logFC from most to least changed gene

Figure 4
Heatmap showing the gene occurrence between chosen differently expressed genes that belongs to “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. The yellow color is associated with gene occurrence in the GO Term. The intensity of the color is corresponding to amount of GO BP terms that each gene belongs to

Figure 5
STRING-generated interaction occurrence between chosen differently expressed genes that belongs to the to “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. The intensity of the edges reflects the strength of interaction score

Figure 6
Functional interaction (FI) between chosen differently expressed genes that belongs to the to “coenzyme metabolic process”, “cofactor biosynthetic process” and “cofactor metabolic process” GO BP terms. In following figure “->” stands for activating/catalyzing, “-|” for inhibition, “-” for FIs extracted from complexes or inputs, and “---” for predicted FIs