This pipeline performs Gene Set Enrichment Analysis (GSEA) using The Broad Institute GSEA tool with MSigDB - Class2: Canonical Pathways gene sets. For a given phenotype subtype, it shows what pathways are significantly enriched in each subtype by comparing gene expression profiles between subtypes. Here, the phenotype is mRNAseq_cNMF subtypes in HNSC-TP. This pipeline has the following features:
-
For each subtype, calculates enrichment scores (ES) using signal to noise (S2N) that checks similarity between subtypes in expression level then calculates p values through permutation test.
-
Lists pathways significantly enriched in each phenotype subtype and their enrichment scores (ES).
-
Lists top 20 core genes enriched in each significant gene set and their enrichment scores (ES).
-
Checks if the top core genes are up-regulated or down-regulated.
-
Checks if the top core genes are high expressed or low expressed.
-
Checks if the top core genes are significantly differently expressed genes.
Table 1. Get Full Table basic data info
basic data info |
---|
Number of Gene Sets: 192 |
Number of samples: 520 |
Original number of Gene Sets: 404 |
Maximum gene set size: 388 |
Table 2. Get Full Table pheno data info
phenotype info |
---|
pheno.type: 1 - 3 :[ clus1 ] 144 |
pheno.type: 2 - 3 :[ clus2 ] 214 |
pheno.type: 3 - 3 :[ clus3 ] 162 |
For the expression subtypes of 18378 genes in 521 samples, GSEA found enriched gene sets in each cluster using 520 gene sets in MSigDB canonical pathways. Top enriched gene sets are listed as below.
-
clus1
-
Top enriched gene sets are BIOCARTA NO1 PATHWAY, BIOCARTA AGR PATHWAY, BIOCARTA ALK PATHWAY, BIOCARTA AT1R PATHWAY, BIOCARTA BIOPEPTIDES PATHWAY, BIOCARTA CARM ER PATHWAY, BIOCARTA HDAC PATHWAY, BIOCARTA EGF PATHWAY, BIOCARTA ERK PATHWAY, BIOCARTA FMLP PATHWAY
-
And common core enriched genes are PPP1R12B, ACTA1, ARHGAP6, MYL2, MYLK, PIP5K1B, ACTN1, ACTN2, ACTN3, CAV1
-
clus2
-
Top enriched gene sets are BIOCARTA CHREBP2 PATHWAY, BIOCARTA WNT PATHWAY, KEGG FATTY ACID METABOLISM, KEGG GLYCINE SERINE AND THREONINE METABOLISM, KEGG VALINE LEUCINE AND ISOLEUCINE DEGRADATION, KEGG LYSINE DEGRADATION, KEGG SELENOAMINO ACID METABOLISM, KEGG STARCH AND SUCROSE METABOLISM, KEGG GLYCOSPHINGOLIPID BIOSYNTHESIS LACTO AND NEOLACTO SERIES, KEGG BUTANOATE METABOLISM
-
And common core enriched genes are UGT1A1, UGT1A10, UGT1A3, UGT1A4, UGT1A6, UGT1A7, UGT1A8, UGT1A9, UGT2A1, UGT2B15
-
clus3
-
Top enriched gene sets are BIOCARTA ERK PATHWAY, BIOCARTA RACCYCD PATHWAY, BIOCARTA PYK2 PATHWAY, BIOCARTA IL1R PATHWAY, KEGG PYRIMIDINE METABOLISM, KEGG AMINO SUGAR AND NUCLEOTIDE SUGAR METABOLISM, KEGG ETHER LIPID METABOLISM, KEGG AMINOACYL TRNA BIOSYNTHESIS, KEGG RIBOSOME, KEGG RNA POLYMERASE
-
And common core enriched genes are IFNG, POMP, PSMA1, PSMA2, PSMA3, PSMA4, PSMA5, PSMA6, PSMA7, PSMB1
Table 3. Get Full Table This table shows top 10 pathways which are significantly enriched in cluster clus1. It displays only significant gene sets satisfying nom.p.val.threshold (-1), fwer.p.val.threshold (-1) , fdr.q.val.threshold (0.25) and the default table is sorted by Normalized Enrichment Score (NES). Further details on NES statistics, please visit The Broad GSEA website.
GeneSet(GS) | Size(#genes) | genes.ES.table | ES | NES | NOM.p.val | FDR.q.val | FWER.p.val | Tag.. | Gene.. | Signal | FDR..median. | glob.p.val |
---|---|---|---|---|---|---|---|---|---|---|---|---|
BIOCARTA NO1 PATHWAY | 28 | genes.ES.table | 0.75 | 1.8 | 0 | 0.032 | 0.29 | 0.25 | 0.061 | 0.24 | 0 | 0 |
BIOCARTA AGR PATHWAY | 35 | genes.ES.table | 0.73 | 2 | 0 | 0.012 | 0.059 | 0.26 | 0.084 | 0.24 | 0 | 0 |
BIOCARTA ALK PATHWAY | 34 | genes.ES.table | 0.64 | 1.8 | 0 | 0.032 | 0.29 | 0.44 | 0.2 | 0.35 | 0 | 0 |
BIOCARTA AT1R PATHWAY | 32 | genes.ES.table | 0.58 | 2.1 | 0.0019 | 0.014 | 0.027 | 0.28 | 0.22 | 0.22 | 0 | 0.002 |
BIOCARTA BIOPEPTIDES PATHWAY | 40 | genes.ES.table | 0.63 | 2.1 | 0 | 0.0093 | 0.028 | 0.1 | 0.013 | 0.099 | 0 | 0 |
BIOCARTA CARM ER PATHWAY | 34 | genes.ES.table | 0.43 | 1.5 | 0.072 | 0.13 | 0.89 | 0.24 | 0.24 | 0.18 | 0.091 | 0 |
BIOCARTA HDAC PATHWAY | 27 | genes.ES.table | 0.64 | 1.8 | 0 | 0.033 | 0.28 | 0.52 | 0.27 | 0.38 | 0 | 0 |
BIOCARTA EGF PATHWAY | 30 | genes.ES.table | 0.57 | 1.8 | 0.0082 | 0.034 | 0.32 | 0.4 | 0.26 | 0.3 | 0 | 0 |
BIOCARTA ERK PATHWAY | 27 | genes.ES.table | 0.46 | 1.5 | 0.069 | 0.12 | 0.86 | 0.26 | 0.22 | 0.2 | 0.08 | 0 |
BIOCARTA FMLP PATHWAY | 35 | genes.ES.table | 0.4 | 1.3 | 0.18 | 0.24 | 0.98 | 0.29 | 0.24 | 0.22 | 0.2 | 0.002 |
Table S1. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | MYL2 | MYL2 | MYL2 | 27 | 1.5 | 0.44 | YES |
2 | ARHGAP6 | ARHGAP6 | ARHGAP6 | 413 | 0.43 | 0.54 | YES |
3 | PPP1R12B | PPP1R12B | PPP1R12B | 1332 | 0.23 | 0.55 | YES |
4 | PIP5K1B | PIP5K1B | PIP5K1B | 1692 | 0.2 | 0.59 | YES |
5 | MYLK | MYLK | MYLK | 1784 | 0.19 | 0.64 | YES |
6 | VCL | VCL | VCL | 3822 | 0.082 | 0.55 | NO |
7 | ARPC1B | ARPC1B | ARPC1B | 4005 | 0.076 | 0.56 | NO |
8 | ARHGAP1 | ARHGAP1 | ARHGAP1 | 4329 | 0.066 | 0.56 | NO |
9 | OPHN1 | OPHN1 | OPHN1 | 4937 | 0.052 | 0.55 | NO |
10 | LIMK1 | LIMK1 | LIMK1 | 5215 | 0.046 | 0.54 | NO |
11 | GSN | GSN | GSN | 5405 | 0.042 | 0.55 | NO |
12 | ROCK1 | ROCK1 | ROCK1 | 5448 | 0.042 | 0.56 | NO |
13 | SRC | SRC | SRC | 5662 | 0.038 | 0.55 | NO |
14 | PFN1 | PFN1 | PFN1 | 6036 | 0.033 | 0.54 | NO |
15 | CFL1 | CFL1 | CFL1 | 6213 | 0.03 | 0.54 | NO |
16 | ACTR2 | ACTR2 | ACTR2 | 6612 | 0.024 | 0.53 | NO |
17 | ARHGAP5 | ARHGAP5 | ARHGAP5 | 6772 | 0.022 | 0.52 | NO |
18 | DIAPH1 | DIAPH1 | DIAPH1 | 7177 | 0.018 | 0.51 | NO |
19 | RHOA | RHOA | RHOA | 7290 | 0.016 | 0.51 | NO |
20 | ARHGEF1 | ARHGEF1 | ARHGEF1 | 7380 | 0.016 | 0.51 | NO |
21 | ARPC2 | ARPC2 | ARPC2 | 7838 | 0.011 | 0.48 | NO |
22 | ACTR3 | ACTR3 | ACTR3 | 8431 | 0.0048 | 0.45 | NO |
23 | PIP5K1A | PIP5K1A | PIP5K1A | 8745 | 0.0017 | 0.44 | NO |
24 | ARHGEF11 | ARHGEF11 | ARHGEF11 | 9077 | -0.0013 | 0.42 | NO |
25 | ARPC4 | ARPC4 | ARPC4 | 9109 | -0.0017 | 0.42 | NO |
26 | ARPC1A | ARPC1A | ARPC1A | 10262 | -0.012 | 0.36 | NO |
27 | ARPC5 | ARPC5 | ARPC5 | 11350 | -0.022 | 0.31 | NO |
28 | ARPC3 | ARPC3 | ARPC3 | 11460 | -0.024 | 0.31 | NO |
29 | BAIAP2 | BAIAP2 | BAIAP2 | 13640 | -0.049 | 0.2 | NO |
30 | ARHGEF5 | ARHGEF5 | ARHGEF5 | 14645 | -0.065 | 0.17 | NO |
31 | ARHGAP4 | ARHGAP4 | ARHGAP4 | 16608 | -0.13 | 0.096 | NO |
Figure S1. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA NO1 PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_NO1_PATHWAY.heatmap_plot.clus1.28.png)
Figure S2. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA NO1 PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_RHO_PATHWAY.volcano_plot.clus1.1.png)
Table S2. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | ACTN2 | ACTN2 | ACTN2 | 21 | 1.5 | 0.26 | YES |
2 | ACTA1 | ACTA1 | ACTA1 | 23 | 1.5 | 0.52 | YES |
3 | ACTN3 | ACTN3 | ACTN3 | 433 | 0.42 | 0.57 | YES |
4 | TNS1 | TNS1 | TNS1 | 883 | 0.29 | 0.6 | YES |
5 | ITGA1 | ITGA1 | ITGA1 | 1307 | 0.23 | 0.61 | YES |
6 | PPP1R12B | PPP1R12B | PPP1R12B | 1332 | 0.23 | 0.65 | YES |
7 | CAV1 | CAV1 | CAV1 | 1791 | 0.19 | 0.66 | YES |
8 | ACTN1 | ACTN1 | ACTN1 | 2421 | 0.14 | 0.65 | YES |
9 | ITGB1 | ITGB1 | ITGB1 | 2836 | 0.12 | 0.65 | YES |
10 | RAPGEF1 | RAPGEF1 | RAPGEF1 | 3030 | 0.11 | 0.65 | YES |
11 | FYN | FYN | FYN | 3039 | 0.11 | 0.67 | YES |
12 | PXN | PXN | PXN | 3089 | 0.11 | 0.69 | YES |
13 | VCL | VCL | VCL | 3822 | 0.082 | 0.66 | NO |
14 | MAPK8 | MAPK8 | MAPK8 | 4036 | 0.075 | 0.66 | NO |
15 | BCAR1 | BCAR1 | BCAR1 | 4296 | 0.067 | 0.66 | NO |
16 | ZYX | ZYX | ZYX | 4367 | 0.065 | 0.67 | NO |
17 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.66 | NO |
18 | PTK2 | PTK2 | PTK2 | 5277 | 0.045 | 0.64 | NO |
19 | RAP1A | RAP1A | RAP1A | 5378 | 0.043 | 0.64 | NO |
20 | ROCK1 | ROCK1 | ROCK1 | 5448 | 0.042 | 0.64 | NO |
21 | SRC | SRC | SRC | 5662 | 0.038 | 0.64 | NO |
22 | JUN | JUN | JUN | 5847 | 0.035 | 0.63 | NO |
23 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.57 | NO |
24 | RHOA | RHOA | RHOA | 7290 | 0.016 | 0.56 | NO |
25 | SOS1 | SOS1 | SOS1 | 7663 | 0.013 | 0.54 | NO |
26 | CRKL | CRKL | CRKL | 7686 | 0.012 | 0.54 | NO |
27 | MAP2K1 | MAP2K1 | MAP2K1 | 8211 | 0.0071 | 0.51 | NO |
28 | CAPN1 | CAPN1 | CAPN1 | 8817 | 0.001 | 0.48 | NO |
29 | CAPNS1 | CAPNS1 | CAPNS1 | 8972 | -0.00034 | 0.47 | NO |
30 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.44 | NO |
31 | HRAS | HRAS | HRAS | 10029 | -0.0099 | 0.42 | NO |
32 | RAF1 | RAF1 | RAF1 | 10715 | -0.016 | 0.38 | NO |
33 | BCR | BCR | BCR | 11943 | -0.029 | 0.32 | NO |
34 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.31 | NO |
35 | MAP2K2 | MAP2K2 | MAP2K2 | 12299 | -0.033 | 0.31 | NO |
36 | CSK | CSK | CSK | 12816 | -0.038 | 0.29 | NO |
37 | CAPNS2 | CAPNS2 | CAPNS2 | 14593 | -0.064 | 0.21 | NO |
Figure S3. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA AGR PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_AGR_PATHWAY.heatmap_plot.clus1.13.png)
Figure S4. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA AGR PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_INTEGRIN_PATHWAY.volcano_plot.clus1.2.png)
Table S3. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | HSPB2 | HSPB2 | HSPB2 | 176 | 0.79 | 0.21 | YES |
2 | MEF2C | MEF2C | MEF2C | 207 | 0.7 | 0.41 | YES |
3 | TGFB3 | TGFB3 | TGFB3 | 743 | 0.32 | 0.47 | YES |
4 | TGFB2 | TGFB2 | TGFB2 | 1098 | 0.26 | 0.52 | YES |
5 | ATF2 | ATF2 | ATF2 | 3433 | 0.095 | 0.42 | NO |
6 | TGFBR1 | TGFBR1 | TGFBR1 | 3516 | 0.092 | 0.45 | NO |
7 | MEF2D | MEF2D | MEF2D | 3558 | 0.091 | 0.47 | NO |
8 | RAPGEF2 | RAPGEF2 | RAPGEF2 | 3679 | 0.087 | 0.49 | NO |
9 | MEF2A | MEF2A | MEF2A | 3907 | 0.079 | 0.5 | NO |
10 | TGFB1 | TGFB1 | TGFB1 | 4308 | 0.067 | 0.49 | NO |
11 | STAT1 | STAT1 | STAT1 | 4431 | 0.063 | 0.5 | NO |
12 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.5 | NO |
13 | MAP3K7 | MAP3K7 | MAP3K7 | 7673 | 0.012 | 0.35 | NO |
14 | ELK1 | ELK1 | ELK1 | 7699 | 0.012 | 0.35 | NO |
15 | CDC42 | CDC42 | CDC42 | 7701 | 0.012 | 0.35 | NO |
16 | CREB1 | CREB1 | CREB1 | 7720 | 0.012 | 0.36 | NO |
17 | MAP2K4 | MAP2K4 | MAP2K4 | 7787 | 0.011 | 0.36 | NO |
18 | MAPKAPK2 | MAPKAPK2 | MAPKAPK2 | 7978 | 0.0094 | 0.35 | NO |
19 | LOC729991-MEF2B | LOC729991-MEF2B | LOC729991-MEF2B | 8169 | 0.0076 | 0.34 | NO |
20 | MAP3K1 | MAP3K1 | MAP3K1 | 8199 | 0.0073 | 0.34 | NO |
21 | MAPK14 | MAPK14 | MAPK14 | 8829 | 0.001 | 0.31 | NO |
22 | RIPK1 | RIPK1 | RIPK1 | 8866 | 0.00062 | 0.3 | NO |
23 | MKNK1 | MKNK1 | MKNK1 | 9066 | -0.0013 | 0.29 | NO |
24 | DAXX | DAXX | DAXX | 9335 | -0.0036 | 0.28 | NO |
25 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.27 | NO |
26 | HRAS | HRAS | HRAS | 10029 | -0.0099 | 0.25 | NO |
27 | MAX | MAX | MAX | 10456 | -0.014 | 0.23 | NO |
28 | MYC | MYC | MYC | 10473 | -0.014 | 0.23 | NO |
29 | MAPKAPK5 | MAPKAPK5 | MAPKAPK5 | 11631 | -0.025 | 0.17 | NO |
30 | RAC1 | RAC1 | RAC1 | 11810 | -0.027 | 0.17 | NO |
31 | TRADD | TRADD | TRADD | 12692 | -0.037 | 0.13 | NO |
32 | DDIT3 | DDIT3 | DDIT3 | 12919 | -0.04 | 0.13 | NO |
33 | TRAF2 | TRAF2 | TRAF2 | 13831 | -0.051 | 0.098 | NO |
34 | MAP3K9 | MAP3K9 | MAP3K9 | 13990 | -0.053 | 0.1 | NO |
35 | MAP3K5 | MAP3K5 | MAP3K5 | 14205 | -0.057 | 0.11 | NO |
36 | RPS6KA5 | RPS6KA5 | RPS6KA5 | 14280 | -0.058 | 0.12 | NO |
37 | HSPB1 | HSPB1 | HSPB1 | 14376 | -0.06 | 0.13 | NO |
38 | MAP2K6 | MAP2K6 | MAP2K6 | 16679 | -0.13 | 0.045 | NO |
39 | PLA2G4A | PLA2G4A | PLA2G4A | 17276 | -0.17 | 0.059 | NO |
Figure S5. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA ALK PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_ALK_PATHWAY.heatmap_plot.clus1.27.png)
Figure S6. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA ALK PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_P38MAPK_PATHWAY.volcano_plot.clus1.3.png)
Table S4. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | ACTA1 | ACTA1 | ACTA1 | 23 | 1.5 | 0.48 | YES |
2 | PGR | PGR | PGR | 435 | 0.42 | 0.59 | YES |
3 | PRKACB | PRKACB | PRKACB | 3285 | 0.1 | 0.46 | NO |
4 | PAQR5 | PAQR5 | PAQR5 | 3910 | 0.079 | 0.45 | NO |
5 | PRKAR2A | PRKAR2A | PRKAR2A | 3922 | 0.079 | 0.48 | NO |
6 | ARPC1B | ARPC1B | ARPC1B | 4005 | 0.076 | 0.5 | NO |
7 | PRKAR2B | PRKAR2B | PRKAR2B | 4318 | 0.066 | 0.5 | NO |
8 | SRC | SRC | SRC | 5662 | 0.038 | 0.44 | NO |
9 | GNB1 | GNB1 | GNB1 | 5805 | 0.036 | 0.44 | NO |
10 | CAP1 | CAP1 | CAP1 | 5880 | 0.035 | 0.45 | NO |
11 | PRKAR1A | PRKAR1A | PRKAR1A | 5926 | 0.034 | 0.46 | NO |
12 | ACTR2 | ACTR2 | ACTR2 | 6612 | 0.024 | 0.43 | NO |
13 | PAQR7 | PAQR7 | PAQR7 | 6814 | 0.022 | 0.42 | NO |
14 | GNAS | GNAS | GNAS | 6865 | 0.021 | 0.43 | NO |
15 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.42 | NO |
16 | RPS6KA1 | RPS6KA1 | RPS6KA1 | 7131 | 0.018 | 0.42 | NO |
17 | GNAI1 | GNAI1 | GNAI1 | 7163 | 0.018 | 0.43 | NO |
18 | ARPC2 | ARPC2 | ARPC2 | 7838 | 0.011 | 0.4 | NO |
19 | ACTR3 | ACTR3 | ACTR3 | 8431 | 0.0048 | 0.36 | NO |
20 | ARPC4 | ARPC4 | ARPC4 | 9109 | -0.0017 | 0.33 | NO |
21 | HRAS | HRAS | HRAS | 10029 | -0.0099 | 0.28 | NO |
22 | ARPC1A | ARPC1A | ARPC1A | 10262 | -0.012 | 0.27 | NO |
23 | ARPC5 | ARPC5 | ARPC5 | 11350 | -0.022 | 0.22 | NO |
24 | ARPC3 | ARPC3 | ARPC3 | 11460 | -0.024 | 0.22 | NO |
25 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.19 | NO |
26 | PIN1 | PIN1 | PIN1 | 12372 | -0.033 | 0.19 | NO |
27 | PRKAR1B | PRKAR1B | PRKAR1B | 12892 | -0.039 | 0.18 | NO |
28 | MYT1 | MYT1 | MYT1 | 14377 | -0.06 | 0.11 | NO |
29 | GNGT1 | GNGT1 | GNGT1 | 14607 | -0.064 | 0.12 | NO |
30 | CCNB1 | CCNB1 | CCNB1 | 14962 | -0.072 | 0.12 | NO |
31 | CDK1 | CDK1 | CDK1 | 15682 | -0.092 | 0.12 | NO |
32 | CDC25C | CDC25C | CDC25C | 15913 | -0.1 | 0.13 | NO |
Figure S7. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA AT1R PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_AT1R_PATHWAY.heatmap_plot.clus1.5.png)
Figure S8. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA AT1R PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_MPR_PATHWAY.volcano_plot.clus1.4.png)
Table S5. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | AGT | AGT | AGT | 234 | 0.63 | 0.26 | YES |
2 | AGTR1 | AGTR1 | AGTR1 | 291 | 0.52 | 0.48 | YES |
3 | PRKCA | PRKCA | PRKCA | 1761 | 0.19 | 0.48 | YES |
4 | PRKCB | PRKCB | PRKCB | 3324 | 0.098 | 0.44 | YES |
5 | ATF2 | ATF2 | ATF2 | 3433 | 0.095 | 0.47 | YES |
6 | MEF2D | MEF2D | MEF2D | 3558 | 0.091 | 0.5 | YES |
7 | MEF2A | MEF2A | MEF2A | 3907 | 0.079 | 0.52 | YES |
8 | MAPK8 | MAPK8 | MAPK8 | 4036 | 0.075 | 0.55 | YES |
9 | EGFR | EGFR | EGFR | 4060 | 0.074 | 0.58 | YES |
10 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.56 | NO |
11 | PTK2 | PTK2 | PTK2 | 5277 | 0.045 | 0.55 | NO |
12 | SRC | SRC | SRC | 5662 | 0.038 | 0.55 | NO |
13 | JUN | JUN | JUN | 5847 | 0.035 | 0.55 | NO |
14 | CALM1 | CALM1 | CALM1 | 6049 | 0.032 | 0.56 | NO |
15 | PTK2B | PTK2B | PTK2B | 6162 | 0.031 | 0.56 | NO |
16 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.52 | NO |
17 | SOS1 | SOS1 | SOS1 | 7663 | 0.013 | 0.5 | NO |
18 | ELK1 | ELK1 | ELK1 | 7699 | 0.012 | 0.5 | NO |
19 | MAP2K4 | MAP2K4 | MAP2K4 | 7787 | 0.011 | 0.5 | NO |
20 | LOC729991-MEF2B | LOC729991-MEF2B | LOC729991-MEF2B | 8169 | 0.0076 | 0.48 | NO |
21 | MAP3K1 | MAP3K1 | MAP3K1 | 8199 | 0.0073 | 0.48 | NO |
22 | MAP2K1 | MAP2K1 | MAP2K1 | 8211 | 0.0071 | 0.48 | NO |
23 | CALM3 | CALM3 | CALM3 | 8268 | 0.0065 | 0.48 | NO |
24 | CALM2 | CALM2 | CALM2 | 8837 | 0.00092 | 0.45 | NO |
25 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.42 | NO |
26 | GNAQ | GNAQ | GNAQ | 9599 | -0.0062 | 0.42 | NO |
27 | HRAS | HRAS | HRAS | 10029 | -0.0099 | 0.4 | NO |
28 | RAF1 | RAF1 | RAF1 | 10715 | -0.016 | 0.37 | NO |
29 | PAK1 | PAK1 | PAK1 | 11209 | -0.021 | 0.35 | NO |
30 | RAC1 | RAC1 | RAC1 | 11810 | -0.027 | 0.33 | NO |
31 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.32 | NO |
32 | MAP2K2 | MAP2K2 | MAP2K2 | 12299 | -0.033 | 0.33 | NO |
Figure S9. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA BIOPEPTIDES PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_BIOPEPTIDES_PATHWAY.heatmap_plot.clus1.8.png)
Figure S10. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA BIOPEPTIDES PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_AT1R_PATHWAY.volcano_plot.clus1.5.png)
Table S6. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CAV3 | CAV3 | CAV3 | 13 | 1.6 | 0.046 | YES |
2 | MYLPF | MYLPF | MYLPF | 20 | 1.6 | 0.091 | YES |
3 | ACTN2 | ACTN2 | ACTN2 | 21 | 1.5 | 0.14 | YES |
4 | MYL2 | MYL2 | MYL2 | 27 | 1.5 | 0.18 | YES |
5 | MYLK2 | MYLK2 | MYLK2 | 127 | 0.97 | 0.2 | YES |
6 | FLNC | FLNC | FLNC | 148 | 0.87 | 0.23 | YES |
7 | THBS4 | THBS4 | THBS4 | 150 | 0.86 | 0.25 | YES |
8 | MYLK3 | MYLK3 | MYLK3 | 204 | 0.71 | 0.27 | YES |
9 | ITGA7 | ITGA7 | ITGA7 | 213 | 0.69 | 0.29 | YES |
10 | EGF | EGF | EGF | 350 | 0.46 | 0.3 | YES |
11 | LAMA2 | LAMA2 | LAMA2 | 355 | 0.46 | 0.31 | YES |
12 | IGF1 | IGF1 | IGF1 | 393 | 0.44 | 0.32 | YES |
13 | ACTN3 | ACTN3 | ACTN3 | 433 | 0.42 | 0.33 | YES |
14 | TNN | TNN | TNN | 460 | 0.41 | 0.34 | YES |
15 | RELN | RELN | RELN | 498 | 0.4 | 0.35 | YES |
16 | ITGB3 | ITGB3 | ITGB3 | 499 | 0.39 | 0.36 | YES |
17 | COL11A1 | COL11A1 | COL11A1 | 574 | 0.37 | 0.37 | YES |
18 | FN1 | FN1 | FN1 | 691 | 0.33 | 0.37 | YES |
19 | ITGA10 | ITGA10 | ITGA10 | 696 | 0.33 | 0.38 | YES |
20 | THBS2 | THBS2 | THBS2 | 728 | 0.32 | 0.39 | YES |
21 | ITGA9 | ITGA9 | ITGA9 | 760 | 0.32 | 0.4 | YES |
22 | ITGA11 | ITGA11 | ITGA11 | 776 | 0.31 | 0.4 | YES |
23 | TNXB | TNXB | TNXB | 788 | 0.31 | 0.41 | YES |
24 | COMP | COMP | COMP | 790 | 0.31 | 0.42 | YES |
25 | VEGFC | VEGFC | VEGFC | 809 | 0.3 | 0.43 | YES |
26 | COL5A1 | COL5A1 | COL5A1 | 849 | 0.3 | 0.44 | YES |
27 | COL4A6 | COL4A6 | COL4A6 | 852 | 0.3 | 0.44 | YES |
28 | COL5A2 | COL5A2 | COL5A2 | 868 | 0.3 | 0.45 | YES |
29 | RASGRF1 | RASGRF1 | RASGRF1 | 955 | 0.28 | 0.46 | YES |
30 | ITGA8 | ITGA8 | ITGA8 | 1080 | 0.26 | 0.46 | YES |
31 | COL3A1 | COL3A1 | COL3A1 | 1082 | 0.26 | 0.46 | YES |
32 | COL6A3 | COL6A3 | COL6A3 | 1090 | 0.26 | 0.47 | YES |
33 | FLT4 | FLT4 | FLT4 | 1120 | 0.26 | 0.48 | YES |
34 | COL5A3 | COL5A3 | COL5A3 | 1157 | 0.25 | 0.48 | YES |
35 | COL1A2 | COL1A2 | COL1A2 | 1162 | 0.25 | 0.49 | YES |
36 | TNC | TNC | TNC | 1206 | 0.24 | 0.49 | YES |
37 | COL1A1 | COL1A1 | COL1A1 | 1214 | 0.24 | 0.5 | YES |
38 | PDGFC | PDGFC | PDGFC | 1273 | 0.24 | 0.5 | YES |
39 | THBS1 | THBS1 | THBS1 | 1291 | 0.24 | 0.51 | YES |
40 | PDGFRA | PDGFRA | PDGFRA | 1300 | 0.23 | 0.52 | YES |
41 | ITGA1 | ITGA1 | ITGA1 | 1307 | 0.23 | 0.52 | YES |
42 | PDGFRB | PDGFRB | PDGFRB | 1311 | 0.23 | 0.53 | YES |
43 | MYL9 | MYL9 | MYL9 | 1337 | 0.23 | 0.54 | YES |
44 | COL4A4 | COL4A4 | COL4A4 | 1391 | 0.22 | 0.54 | YES |
45 | IBSP | IBSP | IBSP | 1411 | 0.22 | 0.54 | YES |
46 | COL6A1 | COL6A1 | COL6A1 | 1412 | 0.22 | 0.55 | YES |
47 | VTN | VTN | VTN | 1414 | 0.22 | 0.56 | YES |
48 | LAMA4 | LAMA4 | LAMA4 | 1415 | 0.22 | 0.56 | YES |
49 | COL6A2 | COL6A2 | COL6A2 | 1456 | 0.22 | 0.57 | YES |
50 | LAMA1 | LAMA1 | LAMA1 | 1507 | 0.21 | 0.57 | YES |
51 | ITGA5 | ITGA5 | ITGA5 | 1523 | 0.21 | 0.58 | YES |
52 | LAMA3 | LAMA3 | LAMA3 | 1549 | 0.21 | 0.58 | YES |
53 | ITGB6 | ITGB6 | ITGB6 | 1584 | 0.21 | 0.59 | YES |
54 | ITGA4 | ITGA4 | ITGA4 | 1679 | 0.2 | 0.59 | YES |
55 | PDGFB | PDGFB | PDGFB | 1705 | 0.19 | 0.59 | YES |
56 | SHC2 | SHC2 | SHC2 | 1722 | 0.19 | 0.6 | YES |
57 | LAMC2 | LAMC2 | LAMC2 | 1725 | 0.19 | 0.6 | YES |
58 | PRKCA | PRKCA | PRKCA | 1761 | 0.19 | 0.6 | YES |
59 | MYLK | MYLK | MYLK | 1784 | 0.19 | 0.61 | YES |
60 | CAV1 | CAV1 | CAV1 | 1791 | 0.19 | 0.61 | YES |
61 | TLN2 | TLN2 | TLN2 | 1804 | 0.18 | 0.62 | YES |
62 | KDR | KDR | KDR | 1913 | 0.18 | 0.62 | YES |
63 | HGF | HGF | HGF | 1924 | 0.18 | 0.62 | YES |
64 | FIGF | FIGF | FIGF | 2007 | 0.17 | 0.62 | YES |
65 | LAMC3 | LAMC3 | LAMC3 | 2023 | 0.17 | 0.63 | YES |
66 | PIK3CG | PIK3CG | PIK3CG | 2105 | 0.16 | 0.63 | YES |
67 | SHC3 | SHC3 | SHC3 | 2142 | 0.16 | 0.63 | YES |
68 | COL2A1 | COL2A1 | COL2A1 | 2157 | 0.16 | 0.63 | YES |
69 | AKT3 | AKT3 | AKT3 | 2213 | 0.16 | 0.64 | YES |
70 | PIK3R1 | PIK3R1 | PIK3R1 | 2283 | 0.15 | 0.64 | YES |
71 | LAMB2 | LAMB2 | LAMB2 | 2380 | 0.14 | 0.64 | YES |
72 | ACTN1 | ACTN1 | ACTN1 | 2421 | 0.14 | 0.64 | YES |
73 | LAMB1 | LAMB1 | LAMB1 | 2443 | 0.14 | 0.64 | YES |
74 | COL4A2 | COL4A2 | COL4A2 | 2459 | 0.14 | 0.64 | YES |
75 | PIK3CD | PIK3CD | PIK3CD | 2506 | 0.14 | 0.64 | YES |
76 | COL4A1 | COL4A1 | COL4A1 | 2515 | 0.14 | 0.65 | YES |
77 | ITGAV | ITGAV | ITGAV | 2581 | 0.13 | 0.65 | YES |
78 | THBS3 | THBS3 | THBS3 | 2610 | 0.13 | 0.65 | YES |
79 | PARVA | PARVA | PARVA | 2648 | 0.13 | 0.65 | YES |
80 | ITGA3 | ITGA3 | ITGA3 | 2672 | 0.13 | 0.66 | YES |
81 | PAK3 | PAK3 | PAK3 | 2767 | 0.12 | 0.65 | YES |
82 | LAMB4 | LAMB4 | LAMB4 | 2771 | 0.12 | 0.66 | YES |
83 | PARVB | PARVB | PARVB | 2804 | 0.12 | 0.66 | YES |
84 | ITGB1 | ITGB1 | ITGB1 | 2836 | 0.12 | 0.66 | YES |
85 | PDGFD | PDGFD | PDGFD | 2933 | 0.12 | 0.66 | YES |
86 | COL6A6 | COL6A6 | COL6A6 | 2945 | 0.12 | 0.66 | YES |
87 | LAMB3 | LAMB3 | LAMB3 | 2948 | 0.12 | 0.66 | YES |
88 | TNR | TNR | TNR | 2961 | 0.12 | 0.67 | YES |
89 | ROCK2 | ROCK2 | ROCK2 | 3027 | 0.11 | 0.67 | YES |
90 | ITGA6 | ITGA6 | ITGA6 | 3029 | 0.11 | 0.67 | YES |
91 | RAPGEF1 | RAPGEF1 | RAPGEF1 | 3030 | 0.11 | 0.67 | YES |
92 | FYN | FYN | FYN | 3039 | 0.11 | 0.68 | YES |
93 | ITGB5 | ITGB5 | ITGB5 | 3043 | 0.11 | 0.68 | YES |
94 | MYL7 | MYL7 | MYL7 | 3048 | 0.11 | 0.68 | YES |
95 | LAMC1 | LAMC1 | LAMC1 | 3058 | 0.11 | 0.68 | YES |
96 | PDGFA | PDGFA | PDGFA | 3072 | 0.11 | 0.69 | YES |
97 | PXN | PXN | PXN | 3089 | 0.11 | 0.69 | YES |
98 | PIK3R5 | PIK3R5 | PIK3R5 | 3156 | 0.1 | 0.69 | YES |
99 | MET | MET | MET | 3208 | 0.1 | 0.69 | YES |
100 | TLN1 | TLN1 | TLN1 | 3272 | 0.1 | 0.69 | YES |
101 | PRKCB | PRKCB | PRKCB | 3324 | 0.098 | 0.69 | YES |
102 | VAV2 | VAV2 | VAV2 | 3329 | 0.098 | 0.69 | YES |
103 | CAV2 | CAV2 | CAV2 | 3540 | 0.091 | 0.68 | NO |
104 | PARVG | PARVG | PARVG | 3584 | 0.09 | 0.68 | NO |
105 | CCND1 | CCND1 | CCND1 | 3631 | 0.088 | 0.68 | NO |
106 | VWF | VWF | VWF | 3747 | 0.085 | 0.68 | NO |
107 | FLNA | FLNA | FLNA | 3768 | 0.084 | 0.68 | NO |
108 | SPP1 | SPP1 | SPP1 | 3808 | 0.082 | 0.68 | NO |
109 | VCL | VCL | VCL | 3822 | 0.082 | 0.68 | NO |
110 | DOCK1 | DOCK1 | DOCK1 | 3825 | 0.082 | 0.68 | NO |
111 | SHC4 | SHC4 | SHC4 | 3854 | 0.081 | 0.69 | NO |
112 | ILK | ILK | ILK | 3889 | 0.08 | 0.69 | NO |
113 | IGF1R | IGF1R | IGF1R | 4006 | 0.075 | 0.68 | NO |
114 | MAPK8 | MAPK8 | MAPK8 | 4036 | 0.075 | 0.68 | NO |
115 | EGFR | EGFR | EGFR | 4060 | 0.074 | 0.68 | NO |
116 | BIRC2 | BIRC2 | BIRC2 | 4193 | 0.07 | 0.68 | NO |
117 | FLT1 | FLT1 | FLT1 | 4248 | 0.068 | 0.68 | NO |
118 | BCAR1 | BCAR1 | BCAR1 | 4296 | 0.067 | 0.68 | NO |
119 | ZYX | ZYX | ZYX | 4367 | 0.065 | 0.68 | NO |
120 | ITGB4 | ITGB4 | ITGB4 | 4485 | 0.061 | 0.67 | NO |
121 | ITGA2 | ITGA2 | ITGA2 | 4709 | 0.057 | 0.66 | NO |
122 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.66 | NO |
123 | PTEN | PTEN | PTEN | 5081 | 0.048 | 0.64 | NO |
124 | SOS2 | SOS2 | SOS2 | 5145 | 0.047 | 0.64 | NO |
125 | VAV3 | VAV3 | VAV3 | 5176 | 0.047 | 0.64 | NO |
126 | MYL12A | MYL12A | MYL12A | 5220 | 0.046 | 0.64 | NO |
127 | PTK2 | PTK2 | PTK2 | 5277 | 0.045 | 0.64 | NO |
128 | RAC2 | RAC2 | RAC2 | 5281 | 0.045 | 0.64 | NO |
129 | CRK | CRK | CRK | 5360 | 0.043 | 0.64 | NO |
130 | RAP1A | RAP1A | RAP1A | 5378 | 0.043 | 0.64 | NO |
131 | ROCK1 | ROCK1 | ROCK1 | 5448 | 0.042 | 0.63 | NO |
132 | PIP5K1C | PIP5K1C | PIP5K1C | 5587 | 0.039 | 0.63 | NO |
133 | SRC | SRC | SRC | 5662 | 0.038 | 0.62 | NO |
134 | VEGFB | VEGFB | VEGFB | 5713 | 0.037 | 0.62 | NO |
135 | JUN | JUN | JUN | 5847 | 0.035 | 0.62 | NO |
136 | CHAD | CHAD | CHAD | 5876 | 0.035 | 0.62 | NO |
137 | PPP1R12A | PPP1R12A | PPP1R12A | 6070 | 0.032 | 0.6 | NO |
138 | AKT1 | AKT1 | AKT1 | 6455 | 0.026 | 0.58 | NO |
139 | PIK3CA | PIK3CA | PIK3CA | 6594 | 0.025 | 0.58 | NO |
140 | CTNNB1 | CTNNB1 | CTNNB1 | 6626 | 0.024 | 0.58 | NO |
141 | GSK3B | GSK3B | GSK3B | 6665 | 0.024 | 0.58 | NO |
142 | ACTB | ACTB | ACTB | 6686 | 0.024 | 0.57 | NO |
143 | XIAP | XIAP | XIAP | 6705 | 0.023 | 0.57 | NO |
144 | ARHGAP5 | ARHGAP5 | ARHGAP5 | 6772 | 0.022 | 0.57 | NO |
145 | CCND2 | CCND2 | CCND2 | 6821 | 0.022 | 0.57 | NO |
146 | RAP1B | RAP1B | RAP1B | 6997 | 0.02 | 0.56 | NO |
147 | CCND3 | CCND3 | CCND3 | 7001 | 0.02 | 0.56 | NO |
148 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.56 | NO |
149 | DIAPH1 | DIAPH1 | DIAPH1 | 7177 | 0.018 | 0.55 | NO |
150 | CAPN2 | CAPN2 | CAPN2 | 7201 | 0.017 | 0.55 | NO |
151 | VASP | VASP | VASP | 7211 | 0.017 | 0.55 | NO |
152 | RHOA | RHOA | RHOA | 7290 | 0.016 | 0.55 | NO |
153 | PAK2 | PAK2 | PAK2 | 7510 | 0.014 | 0.54 | NO |
154 | PDPK1 | PDPK1 | PDPK1 | 7614 | 0.013 | 0.53 | NO |
155 | SOS1 | SOS1 | SOS1 | 7663 | 0.013 | 0.53 | NO |
156 | CRKL | CRKL | CRKL | 7686 | 0.012 | 0.53 | NO |
157 | ELK1 | ELK1 | ELK1 | 7699 | 0.012 | 0.53 | NO |
158 | CDC42 | CDC42 | CDC42 | 7701 | 0.012 | 0.53 | NO |
159 | LAMA5 | LAMA5 | LAMA5 | 7895 | 0.01 | 0.52 | NO |
160 | PPP1CA | PPP1CA | PPP1CA | 7921 | 0.0099 | 0.52 | NO |
161 | GRLF1 | GRLF1 | GRLF1 | 8099 | 0.0084 | 0.51 | NO |
162 | AKT2 | AKT2 | AKT2 | 8176 | 0.0075 | 0.5 | NO |
163 | FLNB | FLNB | FLNB | 8185 | 0.0074 | 0.5 | NO |
164 | MAP2K1 | MAP2K1 | MAP2K1 | 8211 | 0.0071 | 0.5 | NO |
165 | ACTN4 | ACTN4 | ACTN4 | 8225 | 0.007 | 0.5 | NO |
166 | VAV1 | VAV1 | VAV1 | 8577 | 0.0032 | 0.48 | NO |
167 | MAPK9 | MAPK9 | MAPK9 | 9328 | -0.0036 | 0.44 | NO |
168 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.43 | NO |
169 | PPP1CB | PPP1CB | PPP1CB | 9749 | -0.0075 | 0.42 | NO |
170 | MYL12B | MYL12B | MYL12B | 9755 | -0.0075 | 0.42 | NO |
171 | ACTG1 | ACTG1 | ACTG1 | 9993 | -0.0096 | 0.4 | NO |
172 | PIK3R2 | PIK3R2 | PIK3R2 | 10327 | -0.012 | 0.39 | NO |
173 | RAF1 | RAF1 | RAF1 | 10715 | -0.016 | 0.37 | NO |
174 | BAD | BAD | BAD | 10768 | -0.017 | 0.36 | NO |
175 | PPP1CC | PPP1CC | PPP1CC | 10860 | -0.018 | 0.36 | NO |
176 | PAK1 | PAK1 | PAK1 | 11209 | -0.021 | 0.34 | NO |
177 | BCL2 | BCL2 | BCL2 | 11595 | -0.025 | 0.32 | NO |
178 | RAC1 | RAC1 | RAC1 | 11810 | -0.027 | 0.31 | NO |
179 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.29 | NO |
180 | COL11A2 | COL11A2 | COL11A2 | 12443 | -0.034 | 0.28 | NO |
181 | ITGB8 | ITGB8 | ITGB8 | 12634 | -0.036 | 0.27 | NO |
182 | PIK3CB | PIK3CB | PIK3CB | 12946 | -0.04 | 0.25 | NO |
183 | PAK6 | PAK6 | PAK6 | 13455 | -0.046 | 0.22 | NO |
184 | PGF | PGF | PGF | 13942 | -0.053 | 0.2 | NO |
185 | PAK4 | PAK4 | PAK4 | 14104 | -0.055 | 0.19 | NO |
186 | ITGB7 | ITGB7 | ITGB7 | 14110 | -0.055 | 0.19 | NO |
187 | BIRC3 | BIRC3 | BIRC3 | 14392 | -0.06 | 0.18 | NO |
188 | RAC3 | RAC3 | RAC3 | 14667 | -0.066 | 0.17 | NO |
189 | BRAF | BRAF | BRAF | 14967 | -0.072 | 0.15 | NO |
190 | ERBB2 | ERBB2 | ERBB2 | 15187 | -0.078 | 0.14 | NO |
191 | PIK3R3 | PIK3R3 | PIK3R3 | 15707 | -0.092 | 0.12 | NO |
192 | VEGFA | VEGFA | VEGFA | 16216 | -0.11 | 0.092 | NO |
193 | PAK7 | PAK7 | PAK7 | 16361 | -0.12 | 0.088 | NO |
194 | MYL5 | MYL5 | MYL5 | 17135 | -0.16 | 0.05 | NO |
195 | PRKCG | PRKCG | PRKCG | 17294 | -0.17 | 0.046 | NO |
196 | MAPK10 | MAPK10 | MAPK10 | 17701 | -0.2 | 0.029 | NO |
197 | ITGA2B | ITGA2B | ITGA2B | 17971 | -0.24 | 0.022 | NO |
Figure S11. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA CARM ER PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_CARM_ER_PATHWAY.heatmap_plot.clus1.67.png)
Figure S12. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA CARM ER PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.KEGG_FOCAL_ADHESION.volcano_plot.clus1.6.png)
Table S7. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | ACTN2 | ACTN2 | ACTN2 | 21 | 1.5 | 0.21 | YES |
2 | CTNNA3 | CTNNA3 | CTNNA3 | 147 | 0.87 | 0.32 | YES |
3 | SORBS1 | SORBS1 | SORBS1 | 272 | 0.55 | 0.38 | YES |
4 | ACTN3 | ACTN3 | ACTN3 | 433 | 0.42 | 0.43 | YES |
5 | SNAI1 | SNAI1 | SNAI1 | 1059 | 0.26 | 0.43 | YES |
6 | PTPRM | PTPRM | PTPRM | 1226 | 0.24 | 0.46 | YES |
7 | PTPRB | PTPRB | PTPRB | 1641 | 0.2 | 0.46 | YES |
8 | TGFBR2 | TGFBR2 | TGFBR2 | 1838 | 0.18 | 0.48 | YES |
9 | ACTN1 | ACTN1 | ACTN1 | 2421 | 0.14 | 0.46 | YES |
10 | WASF3 | WASF3 | WASF3 | 2454 | 0.14 | 0.48 | YES |
11 | SNAI2 | SNAI2 | SNAI2 | 2510 | 0.14 | 0.5 | YES |
12 | FGFR1 | FGFR1 | FGFR1 | 2699 | 0.13 | 0.5 | YES |
13 | FYN | FYN | FYN | 3039 | 0.11 | 0.5 | YES |
14 | CTNNA2 | CTNNA2 | CTNNA2 | 3149 | 0.1 | 0.51 | YES |
15 | MET | MET | MET | 3208 | 0.1 | 0.52 | YES |
16 | TCF7 | TCF7 | TCF7 | 3418 | 0.095 | 0.52 | YES |
17 | TGFBR1 | TGFBR1 | TGFBR1 | 3516 | 0.092 | 0.52 | YES |
18 | VCL | VCL | VCL | 3822 | 0.082 | 0.52 | YES |
19 | PVRL3 | PVRL3 | PVRL3 | 3831 | 0.082 | 0.53 | YES |
20 | LEF1 | LEF1 | LEF1 | 3954 | 0.078 | 0.53 | YES |
21 | IGF1R | IGF1R | IGF1R | 4006 | 0.075 | 0.54 | YES |
22 | EGFR | EGFR | EGFR | 4060 | 0.074 | 0.55 | YES |
23 | WAS | WAS | WAS | 4432 | 0.063 | 0.54 | YES |
24 | PTPRJ | PTPRJ | PTPRJ | 4544 | 0.06 | 0.54 | YES |
25 | TCF7L1 | TCF7L1 | TCF7L1 | 4626 | 0.058 | 0.54 | YES |
26 | INSR | INSR | INSR | 4738 | 0.056 | 0.54 | YES |
27 | LMO7 | LMO7 | LMO7 | 4778 | 0.055 | 0.55 | YES |
28 | TCF7L2 | TCF7L2 | TCF7L2 | 5002 | 0.05 | 0.54 | NO |
29 | RAC2 | RAC2 | RAC2 | 5281 | 0.045 | 0.53 | NO |
30 | SRC | SRC | SRC | 5662 | 0.038 | 0.52 | NO |
31 | IQGAP1 | IQGAP1 | IQGAP1 | 5792 | 0.036 | 0.52 | NO |
32 | PTPN1 | PTPN1 | PTPN1 | 5846 | 0.035 | 0.52 | NO |
33 | CSNK2A2 | CSNK2A2 | CSNK2A2 | 5861 | 0.035 | 0.52 | NO |
34 | CTNNB1 | CTNNB1 | CTNNB1 | 6626 | 0.024 | 0.48 | NO |
35 | ACTB | ACTB | ACTB | 6686 | 0.024 | 0.48 | NO |
36 | EP300 | EP300 | EP300 | 6825 | 0.022 | 0.48 | NO |
37 | SMAD4 | SMAD4 | SMAD4 | 7043 | 0.02 | 0.47 | NO |
38 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.47 | NO |
39 | RHOA | RHOA | RHOA | 7290 | 0.016 | 0.46 | NO |
40 | CREBBP | CREBBP | CREBBP | 7325 | 0.016 | 0.46 | NO |
41 | MAP3K7 | MAP3K7 | MAP3K7 | 7673 | 0.012 | 0.44 | NO |
42 | CDC42 | CDC42 | CDC42 | 7701 | 0.012 | 0.44 | NO |
43 | ACTN4 | ACTN4 | ACTN4 | 8225 | 0.007 | 0.42 | NO |
44 | FER | FER | FER | 8694 | 0.0022 | 0.39 | NO |
45 | CSNK2A1 | CSNK2A1 | CSNK2A1 | 9020 | -0.0008 | 0.37 | NO |
46 | PTPRF | PTPRF | PTPRF | 9262 | -0.003 | 0.36 | NO |
47 | WASL | WASL | WASL | 9564 | -0.0059 | 0.34 | NO |
48 | SSX2IP | SSX2IP | SSX2IP | 9857 | -0.0084 | 0.33 | NO |
49 | PVRL2 | PVRL2 | PVRL2 | 9860 | -0.0084 | 0.33 | NO |
50 | CTNNA1 | CTNNA1 | CTNNA1 | 9907 | -0.0088 | 0.33 | NO |
51 | ACTG1 | ACTG1 | ACTG1 | 9993 | -0.0096 | 0.33 | NO |
52 | PTPN6 | PTPN6 | PTPN6 | 10178 | -0.011 | 0.32 | NO |
53 | YES1 | YES1 | YES1 | 10290 | -0.012 | 0.31 | NO |
54 | WASF2 | WASF2 | WASF2 | 10360 | -0.013 | 0.31 | NO |
55 | TJP1 | TJP1 | TJP1 | 10497 | -0.014 | 0.31 | NO |
56 | CSNK2B | CSNK2B | CSNK2B | 11656 | -0.026 | 0.25 | NO |
57 | RAC1 | RAC1 | RAC1 | 11810 | -0.027 | 0.24 | NO |
58 | CTNND1 | CTNND1 | CTNND1 | 11973 | -0.029 | 0.24 | NO |
59 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.23 | NO |
60 | SMAD2 | SMAD2 | SMAD2 | 12476 | -0.035 | 0.22 | NO |
61 | FARP2 | FARP2 | FARP2 | 12539 | -0.035 | 0.22 | NO |
62 | MLLT4 | MLLT4 | MLLT4 | 12683 | -0.037 | 0.22 | NO |
63 | CDH1 | CDH1 | CDH1 | 13261 | -0.044 | 0.19 | NO |
64 | PVRL1 | PVRL1 | PVRL1 | 13298 | -0.044 | 0.2 | NO |
65 | WASF1 | WASF1 | WASF1 | 13345 | -0.045 | 0.2 | NO |
66 | ACP1 | ACP1 | ACP1 | 13391 | -0.045 | 0.2 | NO |
67 | SMAD3 | SMAD3 | SMAD3 | 13421 | -0.046 | 0.21 | NO |
68 | BAIAP2 | BAIAP2 | BAIAP2 | 13640 | -0.049 | 0.2 | NO |
69 | NLK | NLK | NLK | 13667 | -0.049 | 0.21 | NO |
70 | PARD3 | PARD3 | PARD3 | 13808 | -0.051 | 0.21 | NO |
71 | RAC3 | RAC3 | RAC3 | 14667 | -0.066 | 0.17 | NO |
72 | ERBB2 | ERBB2 | ERBB2 | 15187 | -0.078 | 0.15 | NO |
73 | PVRL4 | PVRL4 | PVRL4 | 17430 | -0.18 | 0.051 | NO |
Figure S13. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA HDAC PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_HDAC_PATHWAY.heatmap_plot.clus1.26.png)
Figure S14. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA HDAC PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.KEGG_ADHERENS_JUNCTION.volcano_plot.clus1.7.png)
Table S8. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CAMK2A | CAMK2A | CAMK2A | 128 | 0.97 | 0.2 | YES |
2 | CAMK2B | CAMK2B | CAMK2B | 175 | 0.79 | 0.36 | YES |
3 | MAPT | MAPT | MAPT | 222 | 0.66 | 0.5 | YES |
4 | AGT | AGT | AGT | 234 | 0.63 | 0.63 | YES |
5 | PRKCA | PRKCA | PRKCA | 1761 | 0.19 | 0.58 | NO |
6 | MYLK | MYLK | MYLK | 1784 | 0.19 | 0.62 | NO |
7 | STAT4 | STAT4 | STAT4 | 2865 | 0.12 | 0.59 | NO |
8 | FYN | FYN | FYN | 3039 | 0.11 | 0.6 | NO |
9 | PRKCB | PRKCB | PRKCB | 3324 | 0.098 | 0.61 | NO |
10 | CAMK2D | CAMK2D | CAMK2D | 3828 | 0.082 | 0.6 | NO |
11 | MAPK8 | MAPK8 | MAPK8 | 4036 | 0.075 | 0.6 | NO |
12 | STAT2 | STAT2 | STAT2 | 4234 | 0.069 | 0.61 | NO |
13 | STAT1 | STAT1 | STAT1 | 4431 | 0.063 | 0.61 | NO |
14 | STAT5B | STAT5B | STAT5B | 4638 | 0.058 | 0.61 | NO |
15 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.62 | NO |
16 | GNA11 | GNA11 | GNA11 | 5103 | 0.048 | 0.61 | NO |
17 | KNG1 | KNG1 | KNG1 | 5273 | 0.045 | 0.61 | NO |
18 | PLCG1 | PLCG1 | PLCG1 | 5399 | 0.042 | 0.61 | NO |
19 | GNB1 | GNB1 | GNB1 | 5805 | 0.036 | 0.59 | NO |
20 | CALM1 | CALM1 | CALM1 | 6049 | 0.032 | 0.59 | NO |
21 | PTK2B | PTK2B | PTK2B | 6162 | 0.031 | 0.59 | NO |
22 | STAT5A | STAT5A | STAT5A | 6715 | 0.023 | 0.56 | NO |
23 | CAMK2G | CAMK2G | CAMK2G | 6721 | 0.023 | 0.57 | NO |
24 | JAK2 | JAK2 | JAK2 | 6768 | 0.023 | 0.57 | NO |
25 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.56 | NO |
26 | GNAI1 | GNAI1 | GNAI1 | 7163 | 0.018 | 0.56 | NO |
27 | SOS1 | SOS1 | SOS1 | 7663 | 0.013 | 0.53 | NO |
28 | MAP2K1 | MAP2K1 | MAP2K1 | 8211 | 0.0071 | 0.5 | NO |
29 | CALM3 | CALM3 | CALM3 | 8268 | 0.0065 | 0.5 | NO |
30 | MAPK14 | MAPK14 | MAPK14 | 8829 | 0.001 | 0.47 | NO |
31 | CALM2 | CALM2 | CALM2 | 8837 | 0.00092 | 0.47 | NO |
32 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.43 | NO |
33 | STAT6 | STAT6 | STAT6 | 9788 | -0.0078 | 0.42 | NO |
34 | RAF1 | RAF1 | RAF1 | 10715 | -0.016 | 0.37 | NO |
35 | CDK5 | CDK5 | CDK5 | 10935 | -0.018 | 0.37 | NO |
36 | STAT3 | STAT3 | STAT3 | 11471 | -0.024 | 0.34 | NO |
37 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.31 | NO |
38 | MAP2K2 | MAP2K2 | MAP2K2 | 12299 | -0.033 | 0.31 | NO |
39 | F2 | F2 | F2 | 12461 | -0.034 | 0.31 | NO |
40 | GNGT1 | GNGT1 | GNGT1 | 14607 | -0.064 | 0.2 | NO |
Figure S15. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA EGF PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_EGF_PATHWAY.heatmap_plot.clus1.29.png)
Figure S16. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA EGF PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_BIOPEPTIDES_PATHWAY.volcano_plot.clus1.8.png)
Table S9. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | PPP1R3A | PPP1R3A | PPP1R3A | 43 | 1.4 | 0.095 | YES |
2 | PRKAG3 | PRKAG3 | PRKAG3 | 48 | 1.4 | 0.19 | YES |
3 | PHKG1 | PHKG1 | PHKG1 | 141 | 0.92 | 0.25 | YES |
4 | PYGM | PYGM | PYGM | 143 | 0.9 | 0.31 | YES |
5 | SORBS1 | SORBS1 | SORBS1 | 272 | 0.55 | 0.34 | YES |
6 | PRKAA2 | PRKAA2 | PRKAA2 | 359 | 0.46 | 0.37 | YES |
7 | FBP2 | FBP2 | FBP2 | 367 | 0.45 | 0.4 | YES |
8 | PPARGC1A | PPARGC1A | PPARGC1A | 458 | 0.41 | 0.42 | YES |
9 | PPP1R3C | PPP1R3C | PPP1R3C | 546 | 0.38 | 0.44 | YES |
10 | SLC2A4 | SLC2A4 | SLC2A4 | 613 | 0.35 | 0.46 | YES |
11 | PDE3A | PDE3A | PDE3A | 732 | 0.32 | 0.48 | YES |
12 | ACACB | ACACB | ACACB | 1176 | 0.25 | 0.47 | YES |
13 | IRS1 | IRS1 | IRS1 | 1471 | 0.22 | 0.47 | YES |
14 | SHC2 | SHC2 | SHC2 | 1722 | 0.19 | 0.46 | YES |
15 | HK3 | HK3 | HK3 | 1898 | 0.18 | 0.47 | YES |
16 | PIK3CG | PIK3CG | PIK3CG | 2105 | 0.16 | 0.47 | YES |
17 | SHC3 | SHC3 | SHC3 | 2142 | 0.16 | 0.48 | YES |
18 | AKT3 | AKT3 | AKT3 | 2213 | 0.16 | 0.48 | YES |
19 | PIK3R1 | PIK3R1 | PIK3R1 | 2283 | 0.15 | 0.49 | YES |
20 | PRKAB2 | PRKAB2 | PRKAB2 | 2376 | 0.14 | 0.49 | YES |
21 | PIK3CD | PIK3CD | PIK3CD | 2506 | 0.14 | 0.5 | YES |
22 | CBLB | CBLB | CBLB | 2878 | 0.12 | 0.48 | YES |
23 | FOXO1 | FOXO1 | FOXO1 | 2882 | 0.12 | 0.49 | YES |
24 | IRS2 | IRS2 | IRS2 | 3005 | 0.11 | 0.49 | YES |
25 | RHOQ | RHOQ | RHOQ | 3028 | 0.11 | 0.5 | YES |
26 | RAPGEF1 | RAPGEF1 | RAPGEF1 | 3030 | 0.11 | 0.51 | YES |
27 | PIK3R5 | PIK3R5 | PIK3R5 | 3156 | 0.1 | 0.51 | YES |
28 | PRKACB | PRKACB | PRKACB | 3285 | 0.1 | 0.51 | NO |
29 | SHC4 | SHC4 | SHC4 | 3854 | 0.081 | 0.48 | NO |
30 | PRKAR2A | PRKAR2A | PRKAR2A | 3922 | 0.079 | 0.48 | NO |
31 | MAPK8 | MAPK8 | MAPK8 | 4036 | 0.075 | 0.48 | NO |
32 | PYGL | PYGL | PYGL | 4126 | 0.072 | 0.48 | NO |
33 | SH2B2 | SH2B2 | SH2B2 | 4146 | 0.071 | 0.49 | NO |
34 | PRKAR2B | PRKAR2B | PRKAR2B | 4318 | 0.066 | 0.48 | NO |
35 | INSR | INSR | INSR | 4738 | 0.056 | 0.46 | NO |
36 | SHC1 | SHC1 | SHC1 | 4745 | 0.056 | 0.46 | NO |
37 | CBL | CBL | CBL | 5031 | 0.05 | 0.45 | NO |
38 | CALML6 | CALML6 | CALML6 | 5048 | 0.049 | 0.46 | NO |
39 | SOS2 | SOS2 | SOS2 | 5145 | 0.047 | 0.45 | NO |
40 | PCK1 | PCK1 | PCK1 | 5301 | 0.044 | 0.45 | NO |
41 | CRK | CRK | CRK | 5360 | 0.043 | 0.45 | NO |
42 | PHKA1 | PHKA1 | PHKA1 | 5462 | 0.042 | 0.44 | NO |
43 | SOCS4 | SOCS4 | SOCS4 | 5509 | 0.041 | 0.44 | NO |
44 | PRKACA | PRKACA | PRKACA | 5610 | 0.039 | 0.44 | NO |
45 | INPP5K | INPP5K | INPP5K | 5659 | 0.038 | 0.44 | NO |
46 | SOCS3 | SOCS3 | SOCS3 | 5665 | 0.038 | 0.44 | NO |
47 | PRKAA1 | PRKAA1 | PRKAA1 | 5674 | 0.038 | 0.45 | NO |
48 | MTOR | MTOR | MTOR | 5786 | 0.036 | 0.44 | NO |
49 | GYS1 | GYS1 | GYS1 | 5804 | 0.036 | 0.44 | NO |
50 | PTPN1 | PTPN1 | PTPN1 | 5846 | 0.035 | 0.44 | NO |
51 | FLOT1 | FLOT1 | FLOT1 | 5922 | 0.034 | 0.44 | NO |
52 | PRKAR1A | PRKAR1A | PRKAR1A | 5926 | 0.034 | 0.44 | NO |
53 | CALM1 | CALM1 | CALM1 | 6049 | 0.032 | 0.44 | NO |
54 | PYGB | PYGB | PYGB | 6086 | 0.032 | 0.44 | NO |
55 | SOCS2 | SOCS2 | SOCS2 | 6131 | 0.031 | 0.44 | NO |
56 | ARAF | ARAF | ARAF | 6169 | 0.031 | 0.44 | NO |
57 | RPTOR | RPTOR | RPTOR | 6406 | 0.027 | 0.43 | NO |
58 | AKT1 | AKT1 | AKT1 | 6455 | 0.026 | 0.43 | NO |
59 | PIK3CA | PIK3CA | PIK3CA | 6594 | 0.025 | 0.42 | NO |
60 | GSK3B | GSK3B | GSK3B | 6665 | 0.024 | 0.42 | NO |
61 | NRAS | NRAS | NRAS | 6901 | 0.021 | 0.41 | NO |
62 | PHKB | PHKB | PHKB | 6948 | 0.02 | 0.41 | NO |
63 | MAPK1 | MAPK1 | MAPK1 | 7096 | 0.019 | 0.4 | NO |
64 | FBP1 | FBP1 | FBP1 | 7102 | 0.019 | 0.4 | NO |
65 | PPP1R3B | PPP1R3B | PPP1R3B | 7139 | 0.018 | 0.4 | NO |
66 | PDPK1 | PDPK1 | PDPK1 | 7614 | 0.013 | 0.38 | NO |
67 | SOS1 | SOS1 | SOS1 | 7663 | 0.013 | 0.37 | NO |
68 | CRKL | CRKL | CRKL | 7686 | 0.012 | 0.37 | NO |
69 | ELK1 | ELK1 | ELK1 | 7699 | 0.012 | 0.37 | NO |
70 | PRKAG1 | PRKAG1 | PRKAG1 | 7722 | 0.012 | 0.37 | NO |
71 | PPP1CA | PPP1CA | PPP1CA | 7921 | 0.0099 | 0.36 | NO |
72 | HK1 | HK1 | HK1 | 7970 | 0.0095 | 0.36 | NO |
73 | EIF4E2 | EIF4E2 | EIF4E2 | 8038 | 0.0089 | 0.36 | NO |
74 | AKT2 | AKT2 | AKT2 | 8176 | 0.0075 | 0.35 | NO |
75 | EXOC7 | EXOC7 | EXOC7 | 8190 | 0.0074 | 0.35 | NO |
76 | MAP2K1 | MAP2K1 | MAP2K1 | 8211 | 0.0071 | 0.35 | NO |
77 | CALM3 | CALM3 | CALM3 | 8268 | 0.0065 | 0.35 | NO |
78 | TSC1 | TSC1 | TSC1 | 8489 | 0.0041 | 0.34 | NO |
79 | CALM2 | CALM2 | CALM2 | 8837 | 0.00092 | 0.32 | NO |
80 | MKNK1 | MKNK1 | MKNK1 | 9066 | -0.0013 | 0.3 | NO |
81 | PTPRF | PTPRF | PTPRF | 9262 | -0.003 | 0.29 | NO |
82 | MAPK9 | MAPK9 | MAPK9 | 9328 | -0.0036 | 0.29 | NO |
83 | GRB2 | GRB2 | GRB2 | 9523 | -0.0055 | 0.28 | NO |
84 | PKLR | PKLR | PKLR | 9595 | -0.0062 | 0.28 | NO |
85 | PPP1CB | PPP1CB | PPP1CB | 9749 | -0.0075 | 0.27 | NO |
86 | RHEB | RHEB | RHEB | 9938 | -0.0091 | 0.26 | NO |
87 | TSC2 | TSC2 | TSC2 | 10157 | -0.011 | 0.25 | NO |
88 | PIK3R2 | PIK3R2 | PIK3R2 | 10327 | -0.012 | 0.24 | NO |
89 | RPS6KB2 | RPS6KB2 | RPS6KB2 | 10579 | -0.015 | 0.23 | NO |
90 | RAF1 | RAF1 | RAF1 | 10715 | -0.016 | 0.22 | NO |
91 | BAD | BAD | BAD | 10768 | -0.017 | 0.22 | NO |
92 | PPP1CC | PPP1CC | PPP1CC | 10860 | -0.018 | 0.22 | NO |
93 | INPP5D | INPP5D | INPP5D | 11478 | -0.024 | 0.18 | NO |
94 | PRKAG2 | PRKAG2 | PRKAG2 | 11496 | -0.024 | 0.18 | NO |
95 | KRAS | KRAS | KRAS | 11542 | -0.024 | 0.18 | NO |
96 | PRKCI | PRKCI | PRKCI | 11546 | -0.024 | 0.18 | NO |
97 | PRKAB1 | PRKAB1 | PRKAB1 | 11617 | -0.025 | 0.18 | NO |
98 | RPS6KB1 | RPS6KB1 | RPS6KB1 | 11685 | -0.026 | 0.18 | NO |
99 | IKBKB | IKBKB | IKBKB | 11926 | -0.028 | 0.17 | NO |
100 | TRIP10 | TRIP10 | TRIP10 | 12161 | -0.031 | 0.16 | NO |
101 | MAPK3 | MAPK3 | MAPK3 | 12238 | -0.032 | 0.16 | NO |
102 | MAP2K2 | MAP2K2 | MAP2K2 | 12299 | -0.033 | 0.16 | NO |
103 | GYS2 | GYS2 | GYS2 | 12456 | -0.034 | 0.15 | NO |
104 | EIF4E | EIF4E | EIF4E | 12656 | -0.037 | 0.14 | NO |
105 | PRKAR1B | PRKAR1B | PRKAR1B | 12892 | -0.039 | 0.13 | NO |
106 | ACACA | ACACA | ACACA | 12918 | -0.04 | 0.13 | NO |
107 | PIK3CB | PIK3CB | PIK3CB | 12946 | -0.04 | 0.13 | NO |
108 | FLOT2 | FLOT2 | FLOT2 | 13052 | -0.041 | 0.13 | NO |
109 | PHKA2 | PHKA2 | PHKA2 | 13081 | -0.041 | 0.13 | NO |
110 | RPS6 | RPS6 | RPS6 | 13209 | -0.043 | 0.13 | NO |
111 | CALML3 | CALML3 | CALML3 | 13378 | -0.045 | 0.12 | NO |
112 | PRKCZ | PRKCZ | PRKCZ | 13381 | -0.045 | 0.12 | NO |
113 | FASN | FASN | FASN | 13630 | -0.049 | 0.11 | NO |
114 | PCK2 | PCK2 | PCK2 | 14018 | -0.054 | 0.096 | NO |
115 | SOCS1 | SOCS1 | SOCS1 | 14625 | -0.064 | 0.067 | NO |
116 | HK2 | HK2 | HK2 | 14758 | -0.067 | 0.064 | NO |
117 | PHKG2 | PHKG2 | PHKG2 | 14764 | -0.068 | 0.069 | NO |
118 | BRAF | BRAF | BRAF | 14967 | -0.072 | 0.062 | NO |
119 | EIF4EBP1 | EIF4EBP1 | EIF4EBP1 | 15178 | -0.078 | 0.056 | NO |
120 | PIK3R3 | PIK3R3 | PIK3R3 | 15707 | -0.092 | 0.034 | NO |
121 | MKNK2 | MKNK2 | MKNK2 | 15823 | -0.097 | 0.034 | NO |
122 | PPP1R3D | PPP1R3D | PPP1R3D | 15831 | -0.097 | 0.04 | NO |
123 | SREBF1 | SREBF1 | SREBF1 | 15899 | -0.099 | 0.043 | NO |
124 | PDE3B | PDE3B | PDE3B | 15902 | -0.099 | 0.05 | NO |
125 | LIPE | LIPE | LIPE | 16653 | -0.13 | 0.017 | NO |
126 | CBLC | CBLC | CBLC | 17186 | -0.16 | -0.00089 | NO |
127 | GCK | GCK | GCK | 17235 | -0.16 | 0.0076 | NO |
128 | PRKX | PRKX | PRKX | 17497 | -0.18 | 0.0057 | NO |
129 | CALML5 | CALML5 | CALML5 | 17675 | -0.2 | 0.0094 | NO |
130 | MAPK10 | MAPK10 | MAPK10 | 17701 | -0.2 | 0.022 | NO |
131 | IRS4 | IRS4 | IRS4 | 17834 | -0.21 | 0.029 | NO |
Figure S17. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA ERK PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_ERK_PATHWAY.heatmap_plot.clus1.61.png)
Figure S18. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA ERK PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.KEGG_INSULIN_SIGNALING_PATHWAY.volcano_plot.clus1.9.png)
Table S10. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | MYH1 | MYH1 | MYH1 | 16 | 1.6 | 0.056 | YES |
2 | MYLPF | MYLPF | MYLPF | 20 | 1.6 | 0.11 | YES |
3 | ACTN2 | ACTN2 | ACTN2 | 21 | 1.5 | 0.17 | YES |
4 | MYH8 | MYH8 | MYH8 | 25 | 1.5 | 0.22 | YES |
5 | MYL2 | MYL2 | MYL2 | 27 | 1.5 | 0.28 | YES |
6 | MYH7 | MYH7 | MYH7 | 53 | 1.4 | 0.33 | YES |
7 | MYH2 | MYH2 | MYH2 | 68 | 1.3 | 0.37 | YES |
8 | MYH4 | MYH4 | MYH4 | 72 | 1.2 | 0.42 | YES |
9 | MYH6 | MYH6 | MYH6 | 82 | 1.2 | 0.46 | YES |
10 | MYH3 | MYH3 | MYH3 | 83 | 1.2 | 0.5 | YES |
11 | CTNNA3 | CTNNA3 | CTNNA3 | 147 | 0.87 | 0.53 | YES |
12 | MYH13 | MYH13 | MYH13 | 166 | 0.82 | 0.56 | YES |
13 | JAM2 | JAM2 | JAM2 | 368 | 0.45 | 0.57 | YES |
14 | MRAS | MRAS | MRAS | 387 | 0.44 | 0.58 | YES |
15 | ACTN3 | ACTN3 | ACTN3 | 433 | 0.42 | 0.6 | YES |
16 | PRKCQ | PRKCQ | PRKCQ | 582 | 0.37 | 0.6 | YES |
17 | MYH7B | MYH7B | MYH7B | 624 | 0.35 | 0.61 | YES |
18 | MYH15 | MYH15 | MYH15 | 936 | 0.28 | 0.6 | YES |
19 | MPDZ | MPDZ | MPDZ | 1050 | 0.27 | 0.61 | YES |
20 | JAM3 | JAM3 | JAM3 | 1085 | 0.26 | 0.62 | YES |
21 | MYL9 | MYL9 | MYL9 | 1337 | 0.23 | 0.61 | YES |
22 | MYH11 | MYH11 | MYH11 | 1344 | 0.23 | 0.62 | YES |
23 | CLDN5 | CLDN5 | CLDN5 | 1368 | 0.23 | 0.62 | YES |
24 | PRKCA | PRKCA | PRKCA | 1761 | 0.19 | 0.61 | NO |
25 | ZAK | ZAK | ZAK | 1978 | 0.17 | 0.6 | NO |
26 | CLDN14 | CLDN14 | CLDN14 | 2062 | 0.17 | 0.61 | NO |
27 | CLDN11 | CLDN11 | CLDN11 | 2074 | 0.16 | 0.61 | NO |
28 | AKT3 | AKT3 | AKT3 | 2213 | 0.16 | 0.61 | NO |
29 | ACTN1 | ACTN1 | ACTN1 | 2421 | 0.14 | 0.6 | NO |
30 | MYH10 | MYH10 | MYH10 | 2611 | 0.13 | 0.6 | NO |
31 | MYL7 | MYL7 | MYL7 | 3048 | 0.11 | 0.58 | NO |
32 | CTNNA2 | CTNNA2 | CTNNA2 | 3149 | 0.1 | 0.58 | NO |
33 | CLDN9 | CLDN9 | CLDN9 | 3233 | 0.1 | 0.57 | NO |
34 | PRKCB | PRKCB | PRKCB | 3324 | 0.098 | 0.57 | NO |
35 | RRAS | RRAS | RRAS | 3435 | 0.095 | 0.57 | NO |
36 | PPP2R2C | PPP2R2C | PPP2R2C | 3748 | 0.085 | 0.56 | NO |
37 | CLDN19 | CLDN19 | CLDN19 | 3873 | 0.08 | 0.55 | NO |
38 | AMOTL1 | AMOTL1 | AMOTL1 | 4350 | 0.066 | 0.53 | NO |
39 | EPB41L3 | EPB41L3 | EPB41L3 | 4471 | 0.062 | 0.52 | NO |
40 | GNAI2 | GNAI2 | GNAI2 | 4539 | 0.06 | 0.52 | NO |
41 | MYH9 | MYH9 | MYH9 | 4601 | 0.059 | 0.52 | NO |
42 | EPB41L2 | EPB41L2 | EPB41L2 | 4742 | 0.056 | 0.52 | NO |
43 | TJAP1 | TJAP1 | TJAP1 | 4792 | 0.055 | 0.52 | NO |
44 | HCLS1 | HCLS1 | HCLS1 | 5005 | 0.05 | 0.51 | NO |
45 | PTEN | PTEN | PTEN | 5081 | 0.048 | 0.5 | NO |
46 | MYL12A | MYL12A | MYL12A | 5220 | 0.046 | 0.5 | NO |
47 | SRC | SRC | SRC | 5662 | 0.038 | 0.48 | NO |
48 | CSDA | CSDA | CSDA | 5811 | 0.036 | 0.47 | NO |
49 | CSNK2A2 | CSNK2A2 | CSNK2A2 | 5861 | 0.035 | 0.47 | NO |
50 | LLGL1 | LLGL1 | LLGL1 | 5899 | 0.035 | 0.47 | NO |
51 | CTTN | CTTN | CTTN | 6435 | 0.027 | 0.44 | NO |
52 | AKT1 | AKT1 | AKT1 | 6455 | 0.026 | 0.44 | NO |
53 | RRAS2 | RRAS2 | RRAS2 | 6512 | 0.026 | 0.44 | NO |
54 | CTNNB1 | CTNNB1 | CTNNB1 | 6626 | 0.024 | 0.43 | NO |
55 | ACTB | ACTB | ACTB | 6686 | 0.024 | 0.43 | NO |
56 | SPTAN1 | SPTAN1 | SPTAN1 | 6811 | 0.022 | 0.42 | NO |
57 | NRAS | NRAS | NRAS | 6901 | 0.021 | 0.42 | NO |
58 | RAB3B | RAB3B | RAB3B | 7058 | 0.019 | 0.41 | NO |
59 | GNAI1 | GNAI1 | GNAI1 | 7163 | 0.018 | 0.4 | NO |
60 | RHOA | RHOA | RHOA | 7290 | 0.016 | 0.4 | NO |
61 | CLDN23 | CLDN23 | CLDN23 | 7569 | 0.014 | 0.38 | NO |
62 | CDC42 | CDC42 | CDC42 | 7701 | 0.012 | 0.38 | NO |
63 | MAGI2 | MAGI2 | MAGI2 | 7881 | 0.01 | 0.37 | NO |
64 | PPP2R1A | PPP2R1A | PPP2R1A | 7944 | 0.0098 | 0.36 | NO |
65 | AKT2 | AKT2 | AKT2 | 8176 | 0.0075 | 0.35 | NO |
66 | ACTN4 | ACTN4 | ACTN4 | 8225 | 0.007 | 0.35 | NO |
67 | ASH1L | ASH1L | ASH1L | 8289 | 0.0063 | 0.35 | NO |
68 | EXOC4 | EXOC4 | EXOC4 | 8415 | 0.005 | 0.34 | NO |
69 | RAB13 | RAB13 | RAB13 | 8775 | 0.0015 | 0.32 | NO |
70 | PPP2R2D | PPP2R2D | PPP2R2D | 8915 | 0.0002 | 0.31 | NO |
71 | CSNK2A1 | CSNK2A1 | CSNK2A1 | 9020 | -0.0008 | 0.31 | NO |
72 | MPP5 | MPP5 | MPP5 | 9086 | -0.0014 | 0.3 | NO |
73 | GNAI3 | GNAI3 | GNAI3 | 9176 | -0.0022 | 0.3 | NO |
74 | EXOC3 | EXOC3 | EXOC3 | 9512 | -0.0054 | 0.28 | NO |
75 | MYL12B | MYL12B | MYL12B | 9755 | -0.0075 | 0.27 | NO |
76 | CASK | CASK | CASK | 9906 | -0.0088 | 0.26 | NO |
77 | CTNNA1 | CTNNA1 | CTNNA1 | 9907 | -0.0088 | 0.26 | NO |
78 | ACTG1 | ACTG1 | ACTG1 | 9993 | -0.0096 | 0.26 | NO |
79 | PPP2R1B | PPP2R1B | PPP2R1B | 10012 | -0.0097 | 0.26 | NO |
80 | VAPA | VAPA | VAPA | 10027 | -0.0098 | 0.26 | NO |
81 | YES1 | YES1 | YES1 | 10290 | -0.012 | 0.24 | NO |
82 | SYMPK | SYMPK | SYMPK | 10326 | -0.012 | 0.24 | NO |
83 | TJP1 | TJP1 | TJP1 | 10497 | -0.014 | 0.23 | NO |
84 | PPP2CB | PPP2CB | PPP2CB | 10943 | -0.018 | 0.21 | NO |
85 | PPP2CA | PPP2CA | PPP2CA | 11335 | -0.022 | 0.19 | NO |
86 | KRAS | KRAS | KRAS | 11542 | -0.024 | 0.18 | NO |
87 | PRKCI | PRKCI | PRKCI | 11546 | -0.024 | 0.18 | NO |
88 | CSNK2B | CSNK2B | CSNK2B | 11656 | -0.026 | 0.17 | NO |
89 | CDK4 | CDK4 | CDK4 | 11701 | -0.026 | 0.17 | NO |
90 | MAGI1 | MAGI1 | MAGI1 | 12027 | -0.03 | 0.15 | NO |
91 | INADL | INADL | INADL | 12082 | -0.03 | 0.15 | NO |
92 | PARD6G | PARD6G | PARD6G | 12637 | -0.036 | 0.12 | NO |
93 | PRKCH | PRKCH | PRKCH | 12672 | -0.037 | 0.12 | NO |
94 | MLLT4 | MLLT4 | MLLT4 | 12683 | -0.037 | 0.12 | NO |
95 | F11R | F11R | F11R | 13022 | -0.041 | 0.11 | NO |
96 | PRKCD | PRKCD | PRKCD | 13341 | -0.045 | 0.09 | NO |
97 | PRKCZ | PRKCZ | PRKCZ | 13381 | -0.045 | 0.089 | NO |
98 | PRKCE | PRKCE | PRKCE | 13654 | -0.049 | 0.076 | NO |
99 | PARD3 | PARD3 | PARD3 | 13808 | -0.051 | 0.07 | NO |
100 | PPP2R2A | PPP2R2A | PPP2R2A | 14403 | -0.06 | 0.039 | NO |
101 | PPP2R2B | PPP2R2B | PPP2R2B | 14557 | -0.063 | 0.033 | NO |
102 | TJP2 | TJP2 | TJP2 | 14615 | -0.064 | 0.032 | NO |
103 | CLDN15 | CLDN15 | CLDN15 | 14724 | -0.067 | 0.029 | NO |
104 | CLDN2 | CLDN2 | CLDN2 | 14750 | -0.067 | 0.03 | NO |
105 | EPB41 | EPB41 | EPB41 | 14818 | -0.068 | 0.029 | NO |
106 | EPB41L1 | EPB41L1 | EPB41L1 | 15017 | -0.073 | 0.021 | NO |
107 | CLDN17 | CLDN17 | CLDN17 | 15175 | -0.078 | 0.015 | NO |
108 | MAGI3 | MAGI3 | MAGI3 | 15179 | -0.078 | 0.017 | NO |
109 | CLDN6 | CLDN6 | CLDN6 | 15409 | -0.084 | 0.0079 | NO |
110 | CLDN1 | CLDN1 | CLDN1 | 15476 | -0.086 | 0.0074 | NO |
111 | CLDN18 | CLDN18 | CLDN18 | 15737 | -0.094 | -0.0034 | NO |
112 | CLDN20 | CLDN20 | CLDN20 | 16347 | -0.12 | -0.033 | NO |
113 | MYH14 | MYH14 | MYH14 | 16416 | -0.12 | -0.032 | NO |
114 | IGSF5 | IGSF5 | IGSF5 | 16939 | -0.15 | -0.055 | NO |
115 | CLDN8 | CLDN8 | CLDN8 | 17124 | -0.16 | -0.06 | NO |
116 | MYL5 | MYL5 | MYL5 | 17135 | -0.16 | -0.054 | NO |
117 | PRKCG | PRKCG | PRKCG | 17294 | -0.17 | -0.057 | NO |
118 | PARD6A | PARD6A | PARD6A | 17312 | -0.17 | -0.052 | NO |
119 | LLGL2 | LLGL2 | LLGL2 | 17447 | -0.18 | -0.053 | NO |
120 | CLDN16 | CLDN16 | CLDN16 | 17502 | -0.18 | -0.049 | NO |
121 | PARD6B | PARD6B | PARD6B | 17874 | -0.22 | -0.061 | NO |
122 | CLDN10 | CLDN10 | CLDN10 | 18041 | -0.24 | -0.062 | NO |
123 | CRB3 | CRB3 | CRB3 | 18067 | -0.25 | -0.054 | NO |
124 | CLDN4 | CLDN4 | CLDN4 | 18101 | -0.26 | -0.046 | NO |
125 | CLDN7 | CLDN7 | CLDN7 | 18178 | -0.28 | -0.04 | NO |
126 | OCLN | OCLN | OCLN | 18285 | -0.32 | -0.035 | NO |
127 | CGN | CGN | CGN | 18301 | -0.32 | -0.024 | NO |
128 | CLDN3 | CLDN3 | CLDN3 | 18309 | -0.33 | -0.012 | NO |
129 | TJP3 | TJP3 | TJP3 | 18359 | -0.41 | 0.00011 | NO |
Figure S19. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA FMLP PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.BIOCARTA_FMLP_PATHWAY.heatmap_plot.clus1.88.png)
Figure S20. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA FMLP PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus1___Class2-CanonicalPathway.KEGG_TIGHT_JUNCTION.volcano_plot.clus1.10.png)
For the top enriched genes, if you want to check whether they are
-
up-regulated, please check the list of up-regulated genes
-
down-regulated, please check the list of down-regulated genes
For the top enriched genes, if you want to check whether they are
-
highly expressed genes, please check the list of high (top 30%) expressed genes
-
low expressed genes, please check the list of low (bottom 30%) expressed genes
An expression pattern of top(30%)/middle(30%)/low(30%) in this subtype against other subtypes is available in a heatmap
For the top enriched genes, if you want to check whether they are
-
significantly differently expressed genes by eBayes lm fit, please check the list of significant genes
Table 4. Get Full Table This table shows top 10 pathways which are significantly enriched in cluster clus2. It displays only significant gene sets satisfying nom.p.val.threshold (-1), fwer.p.val.threshold (-1) , fdr.q.val.threshold (0.25) and the default table is sorted by Normalized Enrichment Score (NES). Further details on NES statistics, please visit The Broad GSEA website.
GeneSet(GS) | Size(#genes) | genes.ES.table | ES | NES | NOM.p.val | FDR.q.val | FWER.p.val | Tag.. | Gene.. | Signal | FDR..median. | glob.p.val |
---|---|---|---|---|---|---|---|---|---|---|---|---|
BIOCARTA CHREBP2 PATHWAY | 40 | genes.ES.table | 0.32 | 1.3 | 0.18 | 0.84 | 0.99 | 0.12 | 0.12 | 0.11 | 0.73 | 0.3 |
BIOCARTA WNT PATHWAY | 25 | genes.ES.table | 0.4 | 1.5 | 0.088 | 1 | 0.92 | 0.16 | 0.15 | 0.14 | 1 | 0.69 |
KEGG FATTY ACID METABOLISM | 40 | genes.ES.table | 0.42 | 1.3 | 0.18 | 0.81 | 0.99 | 0.38 | 0.31 | 0.26 | 0.72 | 0.29 |
KEGG GLYCINE SERINE AND THREONINE METABOLISM | 29 | genes.ES.table | 0.59 | 1.5 | 0.038 | 1 | 0.93 | 0.55 | 0.25 | 0.42 | 1 | 0.62 |
KEGG VALINE LEUCINE AND ISOLEUCINE DEGRADATION | 43 | genes.ES.table | 0.41 | 1.4 | 0.11 | 1 | 0.95 | 0.56 | 0.41 | 0.33 | 1 | 0.51 |
KEGG LYSINE DEGRADATION | 43 | genes.ES.table | 0.32 | 1.3 | 0.16 | 0.88 | 0.98 | 0.46 | 0.38 | 0.29 | 0.75 | 0.34 |
KEGG SELENOAMINO ACID METABOLISM | 25 | genes.ES.table | 0.43 | 1.3 | 0.15 | 0.87 | 0.99 | 0.28 | 0.24 | 0.21 | 0.76 | 0.33 |
KEGG STARCH AND SUCROSE METABOLISM | 39 | genes.ES.table | 0.57 | 1.5 | 0.061 | 1 | 0.92 | 0.41 | 0.16 | 0.34 | 1 | 0.78 |
KEGG GLYCOSPHINGOLIPID BIOSYNTHESIS LACTO AND NEOLACTO SERIES | 25 | genes.ES.table | 0.52 | 1.4 | 0.093 | 0.92 | 0.97 | 0.28 | 0.13 | 0.24 | 0.78 | 0.36 |
KEGG BUTANOATE METABOLISM | 29 | genes.ES.table | 0.46 | 1.3 | 0.11 | 0.83 | 0.98 | 0.45 | 0.33 | 0.3 | 0.7 | 0.3 |
Table S11. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | UGT1A9 | UGT1A9 | UGT1A9 | 272 | 0.48 | 0.079 | YES |
2 | UGT1A6 | UGT1A6 | UGT1A6 | 273 | 0.48 | 0.17 | YES |
3 | UGT2A1 | UGT2A1 | UGT2A1 | 375 | 0.45 | 0.26 | YES |
4 | CP | CP | CP | 408 | 0.45 | 0.34 | YES |
5 | UGT1A8 | UGT1A8 | UGT1A8 | 623 | 0.4 | 0.41 | YES |
6 | UGT1A3 | UGT1A3 | UGT1A3 | 659 | 0.39 | 0.48 | YES |
7 | UGT1A7 | UGT1A7 | UGT1A7 | 1113 | 0.33 | 0.53 | YES |
8 | UGT1A4 | UGT1A4 | UGT1A4 | 1434 | 0.3 | 0.57 | YES |
9 | UGT1A10 | UGT1A10 | UGT1A10 | 2604 | 0.22 | 0.55 | YES |
10 | UGT1A1 | UGT1A1 | UGT1A1 | 2934 | 0.2 | 0.57 | YES |
11 | UGT2B15 | UGT2B15 | UGT2B15 | 2960 | 0.2 | 0.6 | YES |
12 | PPOX | PPOX | PPOX | 4022 | 0.15 | 0.57 | NO |
13 | UGT2B7 | UGT2B7 | UGT2B7 | 4986 | 0.12 | 0.55 | NO |
14 | FECH | FECH | FECH | 5219 | 0.11 | 0.56 | NO |
15 | MMAB | MMAB | MMAB | 5808 | 0.097 | 0.54 | NO |
16 | GUSB | GUSB | GUSB | 5912 | 0.094 | 0.55 | NO |
17 | CPOX | CPOX | CPOX | 6173 | 0.088 | 0.56 | NO |
18 | EARS2 | EARS2 | EARS2 | 8354 | 0.046 | 0.45 | NO |
19 | BLVRB | BLVRB | BLVRB | 9010 | 0.036 | 0.42 | NO |
20 | EPRS | EPRS | EPRS | 9338 | 0.032 | 0.41 | NO |
21 | COX15 | COX15 | COX15 | 9733 | 0.026 | 0.39 | NO |
22 | FTH1 | FTH1 | FTH1 | 9796 | 0.025 | 0.39 | NO |
23 | UROS | UROS | UROS | 10198 | 0.02 | 0.38 | NO |
24 | ALAD | ALAD | ALAD | 10526 | 0.016 | 0.36 | NO |
25 | COX10 | COX10 | COX10 | 10711 | 0.012 | 0.35 | NO |
26 | BLVRA | BLVRA | BLVRA | 12128 | -0.0068 | 0.28 | NO |
27 | ALAS1 | ALAS1 | ALAS1 | 12322 | -0.0095 | 0.27 | NO |
28 | ALAS2 | ALAS2 | ALAS2 | 13489 | -0.029 | 0.21 | NO |
29 | HMBS | HMBS | HMBS | 13524 | -0.03 | 0.21 | NO |
30 | UROD | UROD | UROD | 13680 | -0.032 | 0.21 | NO |
31 | HMOX1 | HMOX1 | HMOX1 | 15350 | -0.072 | 0.14 | NO |
32 | HMOX2 | HMOX2 | HMOX2 | 16857 | -0.15 | 0.082 | NO |
Figure S21. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA CHREBP2 PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.BIOCARTA_CHREBP2_PATHWAY.heatmap_plot.clus2.19.png)
Figure S22. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA CHREBP2 PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_PORPHYRIN_AND_CHLOROPHYLL_METABOLISM.volcano_plot.clus2.1.png)
Table S12. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | UGT1A9 | UGT1A9 | UGT1A9 | 272 | 0.48 | 0.049 | YES |
2 | UGT1A6 | UGT1A6 | UGT1A6 | 273 | 0.48 | 0.11 | YES |
3 | UGT2A1 | UGT2A1 | UGT2A1 | 375 | 0.45 | 0.17 | YES |
4 | GCK | GCK | GCK | 452 | 0.44 | 0.22 | YES |
5 | UGT1A8 | UGT1A8 | UGT1A8 | 623 | 0.4 | 0.26 | YES |
6 | UGT1A3 | UGT1A3 | UGT1A3 | 659 | 0.39 | 0.32 | YES |
7 | AMY2B | AMY2B | AMY2B | 757 | 0.38 | 0.36 | YES |
8 | UGT1A7 | UGT1A7 | UGT1A7 | 1113 | 0.33 | 0.39 | YES |
9 | AMY2A | AMY2A | AMY2A | 1175 | 0.32 | 0.42 | YES |
10 | UGT1A4 | UGT1A4 | UGT1A4 | 1434 | 0.3 | 0.45 | YES |
11 | AMY1A | AMY1A | AMY1A | 1456 | 0.3 | 0.49 | YES |
12 | ENPP3 | ENPP3 | ENPP3 | 1486 | 0.29 | 0.53 | YES |
13 | GBA3 | GBA3 | GBA3 | 1573 | 0.29 | 0.56 | YES |
14 | UGT1A10 | UGT1A10 | UGT1A10 | 2604 | 0.22 | 0.53 | YES |
15 | UGT1A1 | UGT1A1 | UGT1A1 | 2934 | 0.2 | 0.54 | YES |
16 | UGT2B15 | UGT2B15 | UGT2B15 | 2960 | 0.2 | 0.56 | YES |
17 | TREH | TREH | TREH | 4189 | 0.14 | 0.52 | NO |
18 | ENPP1 | ENPP1 | ENPP1 | 4811 | 0.12 | 0.5 | NO |
19 | UGT2B7 | UGT2B7 | UGT2B7 | 4986 | 0.12 | 0.51 | NO |
20 | UGDH | UGDH | UGDH | 5298 | 0.11 | 0.5 | NO |
21 | GUSB | GUSB | GUSB | 5912 | 0.094 | 0.48 | NO |
22 | GAA | GAA | GAA | 6333 | 0.084 | 0.47 | NO |
23 | AGL | AGL | AGL | 7136 | 0.068 | 0.44 | NO |
24 | GPI | GPI | GPI | 7781 | 0.055 | 0.41 | NO |
25 | PGM2L1 | PGM2L1 | PGM2L1 | 7858 | 0.054 | 0.41 | NO |
26 | HK1 | HK1 | HK1 | 12363 | -0.01 | 0.17 | NO |
27 | HK3 | HK3 | HK3 | 12618 | -0.014 | 0.16 | NO |
28 | GYS2 | GYS2 | GYS2 | 12834 | -0.017 | 0.15 | NO |
29 | GBE1 | GBE1 | GBE1 | 13337 | -0.026 | 0.12 | NO |
30 | UGP2 | UGP2 | UGP2 | 13712 | -0.033 | 0.11 | NO |
31 | PGM1 | PGM1 | PGM1 | 14113 | -0.041 | 0.09 | NO |
32 | PYGB | PYGB | PYGB | 14237 | -0.044 | 0.089 | NO |
33 | PGM2 | PGM2 | PGM2 | 14629 | -0.052 | 0.075 | NO |
34 | HK2 | HK2 | HK2 | 15422 | -0.074 | 0.041 | NO |
35 | GYS1 | GYS1 | GYS1 | 15719 | -0.085 | 0.036 | NO |
36 | MGAM | MGAM | MGAM | 15915 | -0.092 | 0.038 | NO |
37 | GANC | GANC | GANC | 16607 | -0.13 | 0.018 | NO |
38 | PYGL | PYGL | PYGL | 17774 | -0.29 | -0.0072 | NO |
39 | PYGM | PYGM | PYGM | 17782 | -0.29 | 0.032 | NO |
Figure S23. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA WNT PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.BIOCARTA_WNT_PATHWAY.heatmap_plot.clus2.3.png)
Figure S24. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA WNT PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_STARCH_AND_SUCROSE_METABOLISM.volcano_plot.clus2.2.png)
Table S13. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | WNT1 | WNT1 | WNT1 | 811 | 0.37 | 0.13 | YES |
2 | WIF1 | WIF1 | WIF1 | 1169 | 0.32 | 0.26 | YES |
3 | LEF1 | LEF1 | LEF1 | 1810 | 0.27 | 0.35 | YES |
4 | FRAT1 | FRAT1 | FRAT1 | 2763 | 0.21 | 0.4 | YES |
5 | SMAD4 | SMAD4 | SMAD4 | 5530 | 0.1 | 0.3 | NO |
6 | NLK | NLK | NLK | 5827 | 0.096 | 0.32 | NO |
7 | APC | APC | APC | 7490 | 0.061 | 0.26 | NO |
8 | TAB1 | TAB1 | TAB1 | 7535 | 0.06 | 0.29 | NO |
9 | FZD1 | FZD1 | FZD1 | 7710 | 0.057 | 0.3 | NO |
10 | CREBBP | CREBBP | CREBBP | 8581 | 0.043 | 0.28 | NO |
11 | CTNNB1 | CTNNB1 | CTNNB1 | 10078 | 0.022 | 0.21 | NO |
12 | HDAC1 | HDAC1 | HDAC1 | 10180 | 0.02 | 0.21 | NO |
13 | GSK3B | GSK3B | GSK3B | 10297 | 0.018 | 0.21 | NO |
14 | TLE1 | TLE1 | TLE1 | 11447 | 0.0029 | 0.15 | NO |
15 | CSNK2A1 | CSNK2A1 | CSNK2A1 | 12247 | -0.0083 | 0.11 | NO |
16 | CTBP1 | CTBP1 | CTBP1 | 12281 | -0.0088 | 0.11 | NO |
17 | AXIN1 | AXIN1 | AXIN1 | 12288 | -0.0089 | 0.12 | NO |
18 | MAP3K7 | MAP3K7 | MAP3K7 | 12293 | -0.009 | 0.12 | NO |
19 | CSNK1D | CSNK1D | CSNK1D | 12300 | -0.0092 | 0.12 | NO |
20 | PPP2CA | PPP2CA | PPP2CA | 13973 | -0.038 | 0.052 | NO |
21 | CCND1 | CCND1 | CCND1 | 14631 | -0.052 | 0.041 | NO |
22 | MYC | MYC | MYC | 14874 | -0.058 | 0.055 | NO |
23 | BTRC | BTRC | BTRC | 15352 | -0.072 | 0.063 | NO |
24 | CSNK1A1 | CSNK1A1 | CSNK1A1 | 15869 | -0.091 | 0.077 | NO |
25 | DVL1 | DVL1 | DVL1 | 16534 | -0.13 | 0.1 | NO |
Figure S25. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG FATTY ACID METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_FATTY_ACID_METABOLISM.heatmap_plot.clus2.20.png)
Figure S26. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG FATTY ACID METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.BIOCARTA_WNT_PATHWAY.volcano_plot.clus2.3.png)
Table S14. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CHDH | CHDH | CHDH | 24 | 0.67 | 0.12 | YES |
2 | CBS | CBS | CBS | 50 | 0.61 | 0.22 | YES |
3 | GNMT | GNMT | GNMT | 412 | 0.44 | 0.28 | YES |
4 | GATM | GATM | GATM | 589 | 0.41 | 0.34 | YES |
5 | PHGDH | PHGDH | PHGDH | 705 | 0.39 | 0.4 | YES |
6 | BHMT | BHMT | BHMT | 1098 | 0.33 | 0.44 | YES |
7 | PIPOX | PIPOX | PIPOX | 1219 | 0.32 | 0.49 | YES |
8 | MAOB | MAOB | MAOB | 1534 | 0.29 | 0.52 | YES |
9 | PSAT1 | PSAT1 | PSAT1 | 2911 | 0.2 | 0.48 | YES |
10 | MAOA | MAOA | MAOA | 3014 | 0.19 | 0.51 | YES |
11 | GAMT | GAMT | GAMT | 3304 | 0.18 | 0.52 | YES |
12 | DMGDH | DMGDH | DMGDH | 3534 | 0.17 | 0.54 | YES |
13 | GCAT | GCAT | GCAT | 4104 | 0.15 | 0.54 | YES |
14 | SARDH | SARDH | SARDH | 4111 | 0.15 | 0.56 | YES |
15 | AOC2 | AOC2 | AOC2 | 4327 | 0.14 | 0.57 | YES |
16 | PSPH | PSPH | PSPH | 4518 | 0.13 | 0.59 | YES |
17 | DAO | DAO | DAO | 5272 | 0.11 | 0.57 | NO |
18 | SHMT2 | SHMT2 | SHMT2 | 6301 | 0.085 | 0.52 | NO |
19 | SHMT1 | SHMT1 | SHMT1 | 6317 | 0.084 | 0.54 | NO |
20 | AGXT | AGXT | AGXT | 6789 | 0.074 | 0.53 | NO |
21 | SDS | SDS | SDS | 8397 | 0.045 | 0.45 | NO |
22 | DLD | DLD | DLD | 8562 | 0.043 | 0.44 | NO |
23 | CTH | CTH | CTH | 8604 | 0.042 | 0.45 | NO |
24 | ALAS1 | ALAS1 | ALAS1 | 12322 | -0.0095 | 0.25 | NO |
25 | SRR | SRR | SRR | 12940 | -0.019 | 0.22 | NO |
26 | ALAS2 | ALAS2 | ALAS2 | 13489 | -0.029 | 0.19 | NO |
27 | GLYCTK | GLYCTK | GLYCTK | 13525 | -0.03 | 0.2 | NO |
28 | AOC3 | AOC3 | AOC3 | 16991 | -0.16 | 0.037 | NO |
29 | GLDC | GLDC | GLDC | 17383 | -0.22 | 0.053 | NO |
Figure S27. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG GLYCINE SERINE AND THREONINE METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_GLYCINE_SERINE_AND_THREONINE_METABOLISM.heatmap_plot.clus2.4.png)
Figure S28. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG GLYCINE SERINE AND THREONINE METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_GLYCINE_SERINE_AND_THREONINE_METABOLISM.volcano_plot.clus2.4.png)
Table S15. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | NEIL1 | NEIL1 | NEIL1 | 438 | 0.44 | 0.13 | YES |
2 | LIG1 | LIG1 | LIG1 | 4045 | 0.15 | -0.009 | YES |
3 | MUTYH | MUTYH | MUTYH | 4139 | 0.14 | 0.038 | YES |
4 | UNG | UNG | UNG | 4789 | 0.12 | 0.048 | YES |
5 | NEIL3 | NEIL3 | NEIL3 | 4966 | 0.12 | 0.082 | YES |
6 | POLE2 | POLE2 | POLE2 | 5004 | 0.12 | 0.12 | YES |
7 | OGG1 | OGG1 | OGG1 | 5075 | 0.12 | 0.16 | YES |
8 | PARP1 | PARP1 | PARP1 | 5155 | 0.12 | 0.2 | YES |
9 | POLB | POLB | POLB | 5548 | 0.1 | 0.21 | YES |
10 | MBD4 | MBD4 | MBD4 | 5679 | 0.1 | 0.24 | YES |
11 | POLE | POLE | POLE | 5732 | 0.099 | 0.28 | YES |
12 | NTHL1 | NTHL1 | NTHL1 | 5860 | 0.095 | 0.3 | YES |
13 | XRCC1 | XRCC1 | XRCC1 | 5944 | 0.093 | 0.33 | YES |
14 | PCNA | PCNA | PCNA | 6208 | 0.087 | 0.35 | YES |
15 | LIG3 | LIG3 | LIG3 | 6461 | 0.081 | 0.36 | YES |
16 | FEN1 | FEN1 | FEN1 | 6756 | 0.075 | 0.38 | YES |
17 | POLD3 | POLD3 | POLD3 | 6900 | 0.072 | 0.39 | YES |
18 | PARP3 | PARP3 | PARP3 | 7000 | 0.07 | 0.41 | YES |
19 | PARP2 | PARP2 | PARP2 | 7036 | 0.069 | 0.44 | YES |
20 | PARP4 | PARP4 | PARP4 | 7558 | 0.059 | 0.43 | YES |
21 | POLD1 | POLD1 | POLD1 | 7625 | 0.058 | 0.45 | YES |
22 | NEIL2 | NEIL2 | NEIL2 | 8351 | 0.046 | 0.42 | NO |
23 | POLE3 | POLE3 | POLE3 | 8503 | 0.044 | 0.43 | NO |
24 | APEX1 | APEX1 | APEX1 | 8986 | 0.037 | 0.42 | NO |
25 | POLL | POLL | POLL | 9298 | 0.032 | 0.41 | NO |
26 | SMUG1 | SMUG1 | SMUG1 | 9880 | 0.024 | 0.39 | NO |
27 | APEX2 | APEX2 | APEX2 | 11364 | 0.0038 | 0.31 | NO |
28 | POLD2 | POLD2 | POLD2 | 12781 | -0.016 | 0.24 | NO |
29 | TDG | TDG | TDG | 12983 | -0.02 | 0.24 | NO |
30 | MPG | MPG | MPG | 13701 | -0.033 | 0.21 | NO |
31 | POLE4 | POLE4 | POLE4 | 14227 | -0.044 | 0.2 | NO |
32 | POLD4 | POLD4 | POLD4 | 15747 | -0.085 | 0.14 | NO |
Figure S29. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG VALINE LEUCINE AND ISOLEUCINE DEGRADATION.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_VALINE_LEUCINE_AND_ISOLEUCINE_DEGRADATION.heatmap_plot.clus2.7.png)
Figure S30. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG VALINE LEUCINE AND ISOLEUCINE DEGRADATION, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_BASE_EXCISION_REPAIR.volcano_plot.clus2.5.png)
Table S16. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | GSTA1 | GSTA1 | GSTA1 | 133 | 0.54 | 0.038 | YES |
2 | UGT1A9 | UGT1A9 | UGT1A9 | 272 | 0.48 | 0.07 | YES |
3 | UGT1A6 | UGT1A6 | UGT1A6 | 273 | 0.48 | 0.11 | YES |
4 | GSTA2 | GSTA2 | GSTA2 | 359 | 0.46 | 0.14 | YES |
5 | UGT2A1 | UGT2A1 | UGT2A1 | 375 | 0.45 | 0.18 | YES |
6 | ADH1C | ADH1C | ADH1C | 436 | 0.44 | 0.22 | YES |
7 | ADH6 | ADH6 | ADH6 | 496 | 0.43 | 0.25 | YES |
8 | CYP2F1 | CYP2F1 | CYP2F1 | 526 | 0.42 | 0.28 | YES |
9 | UGT1A8 | UGT1A8 | UGT1A8 | 623 | 0.4 | 0.31 | YES |
10 | UGT1A3 | UGT1A3 | UGT1A3 | 659 | 0.39 | 0.34 | YES |
11 | ADH7 | ADH7 | ADH7 | 689 | 0.39 | 0.37 | YES |
12 | ALDH3A1 | ALDH3A1 | ALDH3A1 | 702 | 0.39 | 0.4 | YES |
13 | UGT1A7 | UGT1A7 | UGT1A7 | 1113 | 0.33 | 0.41 | YES |
14 | GSTM4 | GSTM4 | GSTM4 | 1386 | 0.3 | 0.42 | YES |
15 | UGT1A4 | UGT1A4 | UGT1A4 | 1434 | 0.3 | 0.44 | YES |
16 | CYP2E1 | CYP2E1 | CYP2E1 | 1544 | 0.29 | 0.46 | YES |
17 | MGST1 | MGST1 | MGST1 | 1604 | 0.28 | 0.48 | YES |
18 | CYP2S1 | CYP2S1 | CYP2S1 | 1627 | 0.28 | 0.5 | YES |
19 | GSTM2 | GSTM2 | GSTM2 | 1670 | 0.28 | 0.52 | YES |
20 | GSTM3 | GSTM3 | GSTM3 | 1959 | 0.26 | 0.53 | YES |
21 | EPHX1 | EPHX1 | EPHX1 | 2157 | 0.24 | 0.54 | YES |
22 | CYP1B1 | CYP1B1 | CYP1B1 | 2288 | 0.24 | 0.55 | YES |
23 | CYP2C19 | CYP2C19 | CYP2C19 | 2389 | 0.23 | 0.56 | YES |
24 | UGT1A10 | UGT1A10 | UGT1A10 | 2604 | 0.22 | 0.57 | YES |
25 | CYP2C9 | CYP2C9 | CYP2C9 | 2900 | 0.2 | 0.57 | YES |
26 | UGT1A1 | UGT1A1 | UGT1A1 | 2934 | 0.2 | 0.59 | YES |
27 | UGT2B15 | UGT2B15 | UGT2B15 | 2960 | 0.2 | 0.6 | YES |
28 | DHDH | DHDH | DHDH | 3112 | 0.19 | 0.61 | YES |
29 | CYP3A5 | CYP3A5 | CYP3A5 | 3334 | 0.18 | 0.61 | YES |
30 | AKR1C1 | AKR1C1 | AKR1C1 | 3407 | 0.17 | 0.62 | YES |
31 | CYP3A7 | CYP3A7 | CYP3A7 | 3790 | 0.16 | 0.61 | YES |
32 | ADH1A | ADH1A | ADH1A | 4365 | 0.14 | 0.59 | YES |
33 | CYP1A1 | CYP1A1 | CYP1A1 | 4512 | 0.13 | 0.6 | YES |
34 | GSTA3 | GSTA3 | GSTA3 | 4563 | 0.13 | 0.6 | YES |
35 | CYP2B6 | CYP2B6 | CYP2B6 | 4747 | 0.13 | 0.61 | YES |
36 | AKR1C3 | AKR1C3 | AKR1C3 | 4804 | 0.12 | 0.61 | YES |
37 | GSTA4 | GSTA4 | GSTA4 | 4888 | 0.12 | 0.62 | YES |
38 | AKR1C2 | AKR1C2 | AKR1C2 | 4939 | 0.12 | 0.63 | YES |
39 | UGT2B7 | UGT2B7 | UGT2B7 | 4986 | 0.12 | 0.63 | YES |
40 | ALDH3B1 | ALDH3B1 | ALDH3B1 | 6234 | 0.086 | 0.57 | NO |
41 | MGST2 | MGST2 | MGST2 | 6540 | 0.08 | 0.56 | NO |
42 | GSTM1 | GSTM1 | GSTM1 | 6752 | 0.075 | 0.56 | NO |
43 | GSTT1 | GSTT1 | GSTT1 | 6832 | 0.074 | 0.56 | NO |
44 | CYP2C8 | CYP2C8 | CYP2C8 | 7657 | 0.058 | 0.52 | NO |
45 | CYP2C18 | CYP2C18 | CYP2C18 | 7769 | 0.056 | 0.52 | NO |
46 | GSTT2 | GSTT2 | GSTT2 | 7874 | 0.054 | 0.52 | NO |
47 | MGST3 | MGST3 | MGST3 | 8523 | 0.043 | 0.48 | NO |
48 | GSTM5 | GSTM5 | GSTM5 | 8648 | 0.042 | 0.48 | NO |
49 | ADH5 | ADH5 | ADH5 | 8935 | 0.037 | 0.47 | NO |
50 | ADH1B | ADH1B | ADH1B | 9668 | 0.027 | 0.43 | NO |
51 | ADH4 | ADH4 | ADH4 | 9841 | 0.025 | 0.42 | NO |
52 | GSTO2 | GSTO2 | GSTO2 | 9977 | 0.023 | 0.42 | NO |
53 | GSTZ1 | GSTZ1 | GSTZ1 | 11299 | 0.0047 | 0.35 | NO |
54 | GSTK1 | GSTK1 | GSTK1 | 14141 | -0.042 | 0.2 | NO |
55 | GSTO1 | GSTO1 | GSTO1 | 14222 | -0.044 | 0.19 | NO |
56 | GSTP1 | GSTP1 | GSTP1 | 14499 | -0.05 | 0.18 | NO |
57 | ALDH1A3 | ALDH1A3 | ALDH1A3 | 15336 | -0.072 | 0.14 | NO |
58 | ALDH3B2 | ALDH3B2 | ALDH3B2 | 17620 | -0.26 | 0.04 | NO |
Figure S31. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG LYSINE DEGRADATION.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_LYSINE_DEGRADATION.heatmap_plot.clus2.16.png)
Figure S32. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG LYSINE DEGRADATION, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_METABOLISM_OF_XENOBIOTICS_BY_CYTOCHROME_P450.volcano_plot.clus2.6.png)
Table S17. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | ACADSB | ACADSB | ACADSB | 1480 | 0.29 | -0.014 | YES |
2 | IVD | IVD | IVD | 1560 | 0.29 | 0.047 | YES |
3 | ALDH6A1 | ALDH6A1 | ALDH6A1 | 1714 | 0.28 | 0.1 | YES |
4 | BCKDHB | BCKDHB | BCKDHB | 1761 | 0.27 | 0.16 | YES |
5 | ABAT | ABAT | ABAT | 2066 | 0.25 | 0.2 | YES |
6 | MCCC1 | MCCC1 | MCCC1 | 2230 | 0.24 | 0.25 | YES |
7 | PCCA | PCCA | PCCA | 3147 | 0.19 | 0.24 | YES |
8 | ACAT2 | ACAT2 | ACAT2 | 3579 | 0.17 | 0.25 | YES |
9 | BCAT2 | BCAT2 | BCAT2 | 3620 | 0.16 | 0.29 | YES |
10 | PCCB | PCCB | PCCB | 4371 | 0.14 | 0.28 | YES |
11 | IL4I1 | IL4I1 | IL4I1 | 4635 | 0.13 | 0.29 | YES |
12 | HADH | HADH | HADH | 4726 | 0.13 | 0.32 | YES |
13 | ACADS | ACADS | ACADS | 4926 | 0.12 | 0.33 | YES |
14 | ALDH3A2 | ALDH3A2 | ALDH3A2 | 5247 | 0.11 | 0.34 | YES |
15 | EHHADH | EHHADH | EHHADH | 5589 | 0.1 | 0.35 | YES |
16 | OXCT2 | OXCT2 | OXCT2 | 5606 | 0.1 | 0.37 | YES |
17 | HMGCS2 | HMGCS2 | HMGCS2 | 5844 | 0.096 | 0.38 | YES |
18 | MCCC2 | MCCC2 | MCCC2 | 6068 | 0.09 | 0.39 | YES |
19 | ACADM | ACADM | ACADM | 6248 | 0.086 | 0.4 | YES |
20 | HIBADH | HIBADH | HIBADH | 6399 | 0.082 | 0.41 | YES |
21 | AUH | AUH | AUH | 6666 | 0.077 | 0.41 | YES |
22 | ALDH1B1 | ALDH1B1 | ALDH1B1 | 7044 | 0.069 | 0.4 | YES |
23 | HIBCH | HIBCH | HIBCH | 7330 | 0.064 | 0.4 | YES |
24 | DBT | DBT | DBT | 7461 | 0.061 | 0.41 | YES |
25 | ACAD8 | ACAD8 | ACAD8 | 7741 | 0.056 | 0.41 | NO |
26 | BCKDHA | BCKDHA | BCKDHA | 7977 | 0.052 | 0.41 | NO |
27 | DLD | DLD | DLD | 8562 | 0.043 | 0.38 | NO |
28 | OXCT1 | OXCT1 | OXCT1 | 8937 | 0.037 | 0.37 | NO |
29 | ALDH9A1 | ALDH9A1 | ALDH9A1 | 8995 | 0.036 | 0.38 | NO |
30 | ALDH2 | ALDH2 | ALDH2 | 9500 | 0.03 | 0.36 | NO |
31 | MUT | MUT | MUT | 9727 | 0.026 | 0.35 | NO |
32 | MCEE | MCEE | MCEE | 9789 | 0.026 | 0.35 | NO |
33 | HSD17B10 | HSD17B10 | HSD17B10 | 10360 | 0.017 | 0.33 | NO |
34 | ECHS1 | ECHS1 | ECHS1 | 10431 | 0.017 | 0.33 | NO |
35 | ACAT1 | ACAT1 | ACAT1 | 11395 | 0.0035 | 0.27 | NO |
36 | ACAA2 | ACAA2 | ACAA2 | 11529 | 0.0017 | 0.27 | NO |
37 | HMGCS1 | HMGCS1 | HMGCS1 | 11731 | -0.001 | 0.26 | NO |
38 | HADHB | HADHB | HADHB | 12856 | -0.018 | 0.2 | NO |
39 | HMGCL | HMGCL | HMGCL | 13123 | -0.022 | 0.19 | NO |
40 | HADHA | HADHA | HADHA | 13565 | -0.03 | 0.17 | NO |
41 | ACAA1 | ACAA1 | ACAA1 | 15851 | -0.09 | 0.068 | NO |
42 | AOX1 | AOX1 | AOX1 | 16187 | -0.1 | 0.074 | NO |
43 | ALDH7A1 | ALDH7A1 | ALDH7A1 | 17263 | -0.2 | 0.06 | NO |
Figure S33. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG SELENOAMINO ACID METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_SELENOAMINO_ACID_METABOLISM.heatmap_plot.clus2.17.png)
Figure S34. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG SELENOAMINO ACID METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_VALINE_LEUCINE_AND_ISOLEUCINE_DEGRADATION.volcano_plot.clus2.7.png)
Table S18. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | GSTA1 | GSTA1 | GSTA1 | 133 | 0.54 | 0.038 | YES |
2 | UGT1A9 | UGT1A9 | UGT1A9 | 272 | 0.48 | 0.07 | YES |
3 | UGT1A6 | UGT1A6 | UGT1A6 | 273 | 0.48 | 0.11 | YES |
4 | GSTA2 | GSTA2 | GSTA2 | 359 | 0.46 | 0.14 | YES |
5 | UGT2A1 | UGT2A1 | UGT2A1 | 375 | 0.45 | 0.18 | YES |
6 | ADH1C | ADH1C | ADH1C | 436 | 0.44 | 0.22 | YES |
7 | ADH6 | ADH6 | ADH6 | 496 | 0.43 | 0.25 | YES |
8 | UGT1A8 | UGT1A8 | UGT1A8 | 623 | 0.4 | 0.27 | YES |
9 | UGT1A3 | UGT1A3 | UGT1A3 | 659 | 0.39 | 0.31 | YES |
10 | ADH7 | ADH7 | ADH7 | 689 | 0.39 | 0.34 | YES |
11 | ALDH3A1 | ALDH3A1 | ALDH3A1 | 702 | 0.39 | 0.37 | YES |
12 | UGT1A7 | UGT1A7 | UGT1A7 | 1113 | 0.33 | 0.37 | YES |
13 | FMO3 | FMO3 | FMO3 | 1137 | 0.33 | 0.4 | YES |
14 | CYP2A6 | CYP2A6 | CYP2A6 | 1163 | 0.32 | 0.43 | YES |
15 | GSTM4 | GSTM4 | GSTM4 | 1386 | 0.3 | 0.44 | YES |
16 | UGT1A4 | UGT1A4 | UGT1A4 | 1434 | 0.3 | 0.46 | YES |
17 | CYP2D6 | CYP2D6 | CYP2D6 | 1517 | 0.29 | 0.48 | YES |
18 | MAOB | MAOB | MAOB | 1534 | 0.29 | 0.5 | YES |
19 | CYP2E1 | CYP2E1 | CYP2E1 | 1544 | 0.29 | 0.53 | YES |
20 | MGST1 | MGST1 | MGST1 | 1604 | 0.28 | 0.55 | YES |
21 | GSTM2 | GSTM2 | GSTM2 | 1670 | 0.28 | 0.57 | YES |
22 | GSTM3 | GSTM3 | GSTM3 | 1959 | 0.26 | 0.57 | YES |
23 | CYP2C19 | CYP2C19 | CYP2C19 | 2389 | 0.23 | 0.57 | YES |
24 | UGT1A10 | UGT1A10 | UGT1A10 | 2604 | 0.22 | 0.58 | YES |
25 | CYP2C9 | CYP2C9 | CYP2C9 | 2900 | 0.2 | 0.58 | YES |
26 | UGT1A1 | UGT1A1 | UGT1A1 | 2934 | 0.2 | 0.59 | YES |
27 | UGT2B15 | UGT2B15 | UGT2B15 | 2960 | 0.2 | 0.61 | YES |
28 | MAOA | MAOA | MAOA | 3014 | 0.19 | 0.62 | YES |
29 | CYP3A5 | CYP3A5 | CYP3A5 | 3334 | 0.18 | 0.62 | NO |
30 | CYP3A7 | CYP3A7 | CYP3A7 | 3790 | 0.16 | 0.61 | NO |
31 | ADH1A | ADH1A | ADH1A | 4365 | 0.14 | 0.59 | NO |
32 | GSTA3 | GSTA3 | GSTA3 | 4563 | 0.13 | 0.59 | NO |
33 | CYP2B6 | CYP2B6 | CYP2B6 | 4747 | 0.13 | 0.59 | NO |
34 | GSTA4 | GSTA4 | GSTA4 | 4888 | 0.12 | 0.59 | NO |
35 | UGT2B7 | UGT2B7 | UGT2B7 | 4986 | 0.12 | 0.59 | NO |
36 | FMO5 | FMO5 | FMO5 | 5448 | 0.11 | 0.58 | NO |
37 | FMO2 | FMO2 | FMO2 | 5643 | 0.1 | 0.58 | NO |
38 | ALDH3B1 | ALDH3B1 | ALDH3B1 | 6234 | 0.086 | 0.55 | NO |
39 | MGST2 | MGST2 | MGST2 | 6540 | 0.08 | 0.54 | NO |
40 | GSTM1 | GSTM1 | GSTM1 | 6752 | 0.075 | 0.54 | NO |
41 | GSTT1 | GSTT1 | GSTT1 | 6832 | 0.074 | 0.54 | NO |
42 | FMO4 | FMO4 | FMO4 | 7466 | 0.061 | 0.51 | NO |
43 | CYP2C8 | CYP2C8 | CYP2C8 | 7657 | 0.058 | 0.5 | NO |
44 | CYP2C18 | CYP2C18 | CYP2C18 | 7769 | 0.056 | 0.5 | NO |
45 | GSTT2 | GSTT2 | GSTT2 | 7874 | 0.054 | 0.5 | NO |
46 | MGST3 | MGST3 | MGST3 | 8523 | 0.043 | 0.47 | NO |
47 | GSTM5 | GSTM5 | GSTM5 | 8648 | 0.042 | 0.46 | NO |
48 | ADH5 | ADH5 | ADH5 | 8935 | 0.037 | 0.45 | NO |
49 | ADH1B | ADH1B | ADH1B | 9668 | 0.027 | 0.42 | NO |
50 | ADH4 | ADH4 | ADH4 | 9841 | 0.025 | 0.41 | NO |
51 | GSTO2 | GSTO2 | GSTO2 | 9977 | 0.023 | 0.4 | NO |
52 | GSTZ1 | GSTZ1 | GSTZ1 | 11299 | 0.0047 | 0.33 | NO |
53 | GSTK1 | GSTK1 | GSTK1 | 14141 | -0.042 | 0.18 | NO |
54 | GSTO1 | GSTO1 | GSTO1 | 14222 | -0.044 | 0.18 | NO |
55 | GSTP1 | GSTP1 | GSTP1 | 14499 | -0.05 | 0.17 | NO |
56 | ALDH1A3 | ALDH1A3 | ALDH1A3 | 15336 | -0.072 | 0.13 | NO |
57 | FMO1 | FMO1 | FMO1 | 15834 | -0.089 | 0.11 | NO |
58 | AOX1 | AOX1 | AOX1 | 16187 | -0.1 | 0.097 | NO |
59 | ALDH3B2 | ALDH3B2 | ALDH3B2 | 17620 | -0.26 | 0.04 | NO |
Figure S35. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG STARCH AND SUCROSE METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_STARCH_AND_SUCROSE_METABOLISM.heatmap_plot.clus2.2.png)
Figure S36. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG STARCH AND SUCROSE METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_DRUG_METABOLISM_CYTOCHROME_P450.volcano_plot.clus2.8.png)
Table S19. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CACNA1B | CACNA1B | CACNA1B | 65 | 0.59 | 0.095 | YES |
2 | GNB3 | GNB3 | GNB3 | 607 | 0.4 | 0.13 | YES |
3 | GNG3 | GNG3 | GNG3 | 1172 | 0.32 | 0.16 | YES |
4 | PRKX | PRKX | PRKX | 1245 | 0.32 | 0.2 | YES |
5 | TAS2R19 | TAS2R19 | TAS2R19 | 1615 | 0.28 | 0.23 | YES |
6 | TAS2R20 | TAS2R20 | TAS2R20 | 1640 | 0.28 | 0.28 | YES |
7 | TAS2R31 | TAS2R31 | TAS2R31 | 1654 | 0.28 | 0.32 | YES |
8 | TRPM5 | TRPM5 | TRPM5 | 1887 | 0.26 | 0.36 | YES |
9 | TAS1R3 | TAS1R3 | TAS1R3 | 2005 | 0.25 | 0.39 | YES |
10 | SCNN1A | SCNN1A | SCNN1A | 2218 | 0.24 | 0.42 | YES |
11 | TAS2R14 | TAS2R14 | TAS2R14 | 2223 | 0.24 | 0.46 | YES |
12 | ADCY6 | ADCY6 | ADCY6 | 2300 | 0.24 | 0.49 | YES |
13 | TAS2R43 | TAS2R43 | TAS2R43 | 2748 | 0.21 | 0.5 | YES |
14 | TAS2R4 | TAS2R4 | TAS2R4 | 2802 | 0.2 | 0.54 | YES |
15 | SCNN1B | SCNN1B | SCNN1B | 3268 | 0.18 | 0.54 | YES |
16 | PLCB2 | PLCB2 | PLCB2 | 3703 | 0.16 | 0.54 | YES |
17 | TAS1R1 | TAS1R1 | TAS1R1 | 3897 | 0.15 | 0.56 | YES |
18 | TAS2R10 | TAS2R10 | TAS2R10 | 4044 | 0.15 | 0.57 | YES |
19 | TAS2R5 | TAS2R5 | TAS2R5 | 4335 | 0.14 | 0.58 | YES |
20 | ADCY8 | ADCY8 | ADCY8 | 4769 | 0.13 | 0.58 | NO |
21 | SCNN1G | SCNN1G | SCNN1G | 6077 | 0.09 | 0.52 | NO |
22 | TAS2R13 | TAS2R13 | TAS2R13 | 6611 | 0.078 | 0.51 | NO |
23 | TAS2R38 | TAS2R38 | TAS2R38 | 8962 | 0.037 | 0.39 | NO |
24 | PDE1A | PDE1A | PDE1A | 9282 | 0.033 | 0.37 | NO |
25 | GNAS | GNAS | GNAS | 9615 | 0.028 | 0.36 | NO |
26 | PRKACA | PRKACA | PRKACA | 9901 | 0.024 | 0.35 | NO |
27 | CACNA1A | CACNA1A | CACNA1A | 10676 | 0.013 | 0.31 | NO |
28 | PRKACB | PRKACB | PRKACB | 12398 | -0.011 | 0.22 | NO |
29 | ITPR3 | ITPR3 | ITPR3 | 13558 | -0.03 | 0.16 | NO |
30 | GNB1 | GNB1 | GNB1 | 13847 | -0.036 | 0.15 | NO |
31 | GRM4 | GRM4 | GRM4 | 14111 | -0.041 | 0.14 | NO |
32 | ACCN1 | ACCN1 | ACCN1 | 14657 | -0.053 | 0.12 | NO |
33 | ADCY4 | ADCY4 | ADCY4 | 14723 | -0.055 | 0.13 | NO |
34 | KCNB1 | KCNB1 | KCNB1 | 17035 | -0.17 | 0.028 | NO |
35 | TAS2R60 | TAS2R60 | TAS2R60 | 17667 | -0.27 | 0.038 | NO |
Figure S37. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG GLYCOSPHINGOLIPID BIOSYNTHESIS LACTO AND NEOLACTO SERIES.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_GLYCOSPHINGOLIPID_BIOSYNTHESIS_LACTO_AND_NEOLACTO_SERIES.heatmap_plot.clus2.12.png)
Figure S38. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG GLYCOSPHINGOLIPID BIOSYNTHESIS LACTO AND NEOLACTO SERIES, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_TASTE_TRANSDUCTION.volcano_plot.clus2.9.png)
Table S20. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | SMC1B | SMC1B | SMC1B | 219 | 0.5 | 0.029 | YES |
2 | E2F5 | E2F5 | E2F5 | 238 | 0.49 | 0.068 | YES |
3 | CDKN2C | CDKN2C | CDKN2C | 1029 | 0.34 | 0.053 | YES |
4 | WEE2 | WEE2 | WEE2 | 1518 | 0.29 | 0.05 | YES |
5 | CDC14A | CDC14A | CDC14A | 1595 | 0.28 | 0.069 | YES |
6 | TFDP2 | TFDP2 | TFDP2 | 1865 | 0.26 | 0.076 | YES |
7 | CCNB3 | CCNB3 | CCNB3 | 1879 | 0.26 | 0.097 | YES |
8 | CDKN2A | CDKN2A | CDKN2A | 2313 | 0.23 | 0.092 | YES |
9 | CDC7 | CDC7 | CDC7 | 2720 | 0.21 | 0.087 | YES |
10 | E2F1 | E2F1 | E2F1 | 2800 | 0.2 | 0.099 | YES |
11 | CDC25C | CDC25C | CDC25C | 2816 | 0.2 | 0.12 | YES |
12 | TP53 | TP53 | TP53 | 2938 | 0.2 | 0.12 | YES |
13 | E2F2 | E2F2 | E2F2 | 2982 | 0.19 | 0.14 | YES |
14 | MCM2 | MCM2 | MCM2 | 3126 | 0.19 | 0.14 | YES |
15 | CHEK2 | CHEK2 | CHEK2 | 3153 | 0.19 | 0.16 | YES |
16 | ESPL1 | ESPL1 | ESPL1 | 3323 | 0.18 | 0.16 | YES |
17 | CCNE2 | CCNE2 | CCNE2 | 3474 | 0.17 | 0.17 | YES |
18 | CDC14B | CDC14B | CDC14B | 3512 | 0.17 | 0.18 | YES |
19 | PTTG2 | PTTG2 | PTTG2 | 3610 | 0.16 | 0.19 | YES |
20 | ORC1L | ORC1L | ORC1L | 3765 | 0.16 | 0.2 | YES |
21 | CDC45 | CDC45 | CDC45 | 3819 | 0.16 | 0.2 | YES |
22 | CDC25A | CDC25A | CDC25A | 3906 | 0.15 | 0.21 | YES |
23 | CDK1 | CDK1 | CDK1 | 4272 | 0.14 | 0.2 | YES |
24 | MCM6 | MCM6 | MCM6 | 4298 | 0.14 | 0.21 | YES |
25 | BUB1B | BUB1B | BUB1B | 4333 | 0.14 | 0.22 | YES |
26 | CDC6 | CDC6 | CDC6 | 4390 | 0.14 | 0.23 | YES |
27 | MCM5 | MCM5 | MCM5 | 4576 | 0.13 | 0.23 | YES |
28 | MAD2L1 | MAD2L1 | MAD2L1 | 4677 | 0.13 | 0.24 | YES |
29 | MCM7 | MCM7 | MCM7 | 4758 | 0.13 | 0.24 | YES |
30 | MCM3 | MCM3 | MCM3 | 5009 | 0.12 | 0.24 | YES |
31 | ATR | ATR | ATR | 5114 | 0.12 | 0.24 | YES |
32 | CDKN1B | CDKN1B | CDKN1B | 5142 | 0.12 | 0.25 | YES |
33 | TTK | TTK | TTK | 5189 | 0.11 | 0.26 | YES |
34 | STAG1 | STAG1 | STAG1 | 5198 | 0.11 | 0.27 | YES |
35 | ORC6L | ORC6L | ORC6L | 5264 | 0.11 | 0.27 | YES |
36 | RBL1 | RBL1 | RBL1 | 5280 | 0.11 | 0.28 | YES |
37 | MDM2 | MDM2 | MDM2 | 5379 | 0.11 | 0.28 | YES |
38 | CDK2 | CDK2 | CDK2 | 5512 | 0.1 | 0.28 | YES |
39 | SMAD4 | SMAD4 | SMAD4 | 5530 | 0.1 | 0.29 | YES |
40 | WEE1 | WEE1 | WEE1 | 5555 | 0.1 | 0.3 | YES |
41 | SMAD2 | SMAD2 | SMAD2 | 5785 | 0.098 | 0.3 | YES |
42 | GADD45G | GADD45G | GADD45G | 5882 | 0.095 | 0.3 | YES |
43 | ORC2L | ORC2L | ORC2L | 5914 | 0.094 | 0.3 | YES |
44 | ATM | ATM | ATM | 6023 | 0.091 | 0.3 | YES |
45 | CCNB2 | CCNB2 | CCNB2 | 6069 | 0.09 | 0.31 | YES |
46 | CDK4 | CDK4 | CDK4 | 6105 | 0.089 | 0.32 | YES |
47 | PCNA | PCNA | PCNA | 6208 | 0.087 | 0.32 | YES |
48 | BUB1 | BUB1 | BUB1 | 6221 | 0.087 | 0.32 | YES |
49 | CDKN1C | CDKN1C | CDKN1C | 6408 | 0.082 | 0.32 | YES |
50 | CCNA2 | CCNA2 | CCNA2 | 6507 | 0.08 | 0.32 | YES |
51 | CDKN2D | CDKN2D | CDKN2D | 6603 | 0.078 | 0.32 | YES |
52 | ORC4L | ORC4L | ORC4L | 6699 | 0.076 | 0.32 | YES |
53 | HDAC2 | HDAC2 | HDAC2 | 6706 | 0.076 | 0.33 | YES |
54 | SMAD3 | SMAD3 | SMAD3 | 6743 | 0.075 | 0.33 | YES |
55 | ORC3L | ORC3L | ORC3L | 6759 | 0.075 | 0.34 | YES |
56 | ANAPC7 | ANAPC7 | ANAPC7 | 6760 | 0.075 | 0.34 | YES |
57 | SMC3 | SMC3 | SMC3 | 7233 | 0.066 | 0.32 | NO |
58 | ANAPC10 | ANAPC10 | ANAPC10 | 7570 | 0.059 | 0.31 | NO |
59 | SMC1A | SMC1A | SMC1A | 7614 | 0.058 | 0.31 | NO |
60 | MCM4 | MCM4 | MCM4 | 7643 | 0.058 | 0.32 | NO |
61 | SKP2 | SKP2 | SKP2 | 7678 | 0.057 | 0.32 | NO |
62 | PTTG1 | PTTG1 | PTTG1 | 7762 | 0.056 | 0.32 | NO |
63 | CHEK1 | CHEK1 | CHEK1 | 7812 | 0.055 | 0.32 | NO |
64 | CDC26 | CDC26 | CDC26 | 7903 | 0.053 | 0.32 | NO |
65 | ANAPC1 | ANAPC1 | ANAPC1 | 7963 | 0.052 | 0.32 | NO |
66 | ORC5L | ORC5L | ORC5L | 8248 | 0.047 | 0.31 | NO |
67 | E2F3 | E2F3 | E2F3 | 8378 | 0.046 | 0.31 | NO |
68 | EP300 | EP300 | EP300 | 8457 | 0.044 | 0.31 | NO |
69 | STAG2 | STAG2 | STAG2 | 8501 | 0.044 | 0.31 | NO |
70 | CREBBP | CREBBP | CREBBP | 8581 | 0.043 | 0.31 | NO |
71 | RAD21 | RAD21 | RAD21 | 8991 | 0.037 | 0.29 | NO |
72 | PLK1 | PLK1 | PLK1 | 9036 | 0.036 | 0.29 | NO |
73 | CCNB1 | CCNB1 | CCNB1 | 9168 | 0.034 | 0.28 | NO |
74 | CDC23 | CDC23 | CDC23 | 9278 | 0.033 | 0.28 | NO |
75 | CCNE1 | CCNE1 | CCNE1 | 9669 | 0.027 | 0.26 | NO |
76 | ANAPC5 | ANAPC5 | ANAPC5 | 9747 | 0.026 | 0.26 | NO |
77 | CCNH | CCNH | CCNH | 9826 | 0.025 | 0.26 | NO |
78 | RBL2 | RBL2 | RBL2 | 9873 | 0.024 | 0.26 | NO |
79 | HDAC1 | HDAC1 | HDAC1 | 10180 | 0.02 | 0.24 | NO |
80 | GSK3B | GSK3B | GSK3B | 10297 | 0.018 | 0.24 | NO |
81 | FZR1 | FZR1 | FZR1 | 10370 | 0.017 | 0.23 | NO |
82 | PRKDC | PRKDC | PRKDC | 10417 | 0.017 | 0.23 | NO |
83 | ANAPC13 | ANAPC13 | ANAPC13 | 10459 | 0.016 | 0.23 | NO |
84 | CCND3 | CCND3 | CCND3 | 10476 | 0.016 | 0.23 | NO |
85 | ANAPC11 | ANAPC11 | ANAPC11 | 10601 | 0.014 | 0.23 | NO |
86 | BUB3 | BUB3 | BUB3 | 10638 | 0.014 | 0.22 | NO |
87 | ABL1 | ABL1 | ABL1 | 11281 | 0.005 | 0.19 | NO |
88 | YWHAH | YWHAH | YWHAH | 11375 | 0.0037 | 0.18 | NO |
89 | ANAPC2 | ANAPC2 | ANAPC2 | 11577 | 0.0011 | 0.17 | NO |
90 | CUL1 | CUL1 | CUL1 | 11611 | 0.0006 | 0.17 | NO |
91 | PKMYT1 | PKMYT1 | PKMYT1 | 11627 | 0.00031 | 0.17 | NO |
92 | SKP1 | SKP1 | SKP1 | 11726 | -0.00091 | 0.17 | NO |
93 | ANAPC4 | ANAPC4 | ANAPC4 | 12107 | -0.0066 | 0.15 | NO |
94 | YWHAE | YWHAE | YWHAE | 12294 | -0.0091 | 0.14 | NO |
95 | TFDP1 | TFDP1 | TFDP1 | 12352 | -0.01 | 0.14 | NO |
96 | MAD2L2 | MAD2L2 | MAD2L2 | 12442 | -0.011 | 0.13 | NO |
97 | CDC20 | CDC20 | CDC20 | 12639 | -0.014 | 0.12 | NO |
98 | RBX1 | RBX1 | RBX1 | 13095 | -0.022 | 0.098 | NO |
99 | YWHAG | YWHAG | YWHAG | 13564 | -0.03 | 0.075 | NO |
100 | CDKN2B | CDKN2B | CDKN2B | 13773 | -0.034 | 0.066 | NO |
101 | CDC27 | CDC27 | CDC27 | 14156 | -0.042 | 0.048 | NO |
102 | GADD45A | GADD45A | GADD45A | 14208 | -0.043 | 0.049 | NO |
103 | E2F4 | E2F4 | E2F4 | 14217 | -0.044 | 0.052 | NO |
104 | YWHAQ | YWHAQ | YWHAQ | 14344 | -0.046 | 0.049 | NO |
105 | CCND2 | CCND2 | CCND2 | 14468 | -0.049 | 0.046 | NO |
106 | CDC16 | CDC16 | CDC16 | 14606 | -0.052 | 0.043 | NO |
107 | RB1 | RB1 | RB1 | 14607 | -0.052 | 0.048 | NO |
108 | CCND1 | CCND1 | CCND1 | 14631 | -0.052 | 0.05 | NO |
109 | YWHAB | YWHAB | YWHAB | 14736 | -0.055 | 0.049 | NO |
110 | ZBTB17 | ZBTB17 | ZBTB17 | 14764 | -0.056 | 0.052 | NO |
111 | MYC | MYC | MYC | 14874 | -0.058 | 0.051 | NO |
112 | MAD1L1 | MAD1L1 | MAD1L1 | 14943 | -0.06 | 0.052 | NO |
113 | CDK7 | CDK7 | CDK7 | 15570 | -0.079 | 0.024 | NO |
114 | GADD45B | GADD45B | GADD45B | 15927 | -0.093 | 0.013 | NO |
115 | YWHAZ | YWHAZ | YWHAZ | 16009 | -0.096 | 0.016 | NO |
116 | CDC25B | CDC25B | CDC25B | 16117 | -0.1 | 0.018 | NO |
117 | TGFB3 | TGFB3 | TGFB3 | 16310 | -0.11 | 0.017 | NO |
118 | TGFB1 | TGFB1 | TGFB1 | 16576 | -0.13 | 0.013 | NO |
119 | CDKN1A | CDKN1A | CDKN1A | 16627 | -0.13 | 0.021 | NO |
120 | TGFB2 | TGFB2 | TGFB2 | 16937 | -0.16 | 0.017 | NO |
121 | CDK6 | CDK6 | CDK6 | 17185 | -0.19 | 0.019 | NO |
122 | SFN | SFN | SFN | 17381 | -0.21 | 0.026 | NO |
123 | CCNA1 | CCNA1 | CCNA1 | 17935 | -0.34 | 0.023 | NO |
Figure S39. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG BUTANOATE METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_BUTANOATE_METABOLISM.heatmap_plot.clus2.15.png)
Figure S40. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG BUTANOATE METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus2___Class2-CanonicalPathway.KEGG_CELL_CYCLE.volcano_plot.clus2.10.png)
For the top enriched genes, if you want to check whether they are
-
up-regulated, please check the list of up-regulated genes
-
down-regulated, please check the list of down-regulated genes
For the top enriched genes, if you want to check whether they are
-
highly expressed genes, please check the list of high (top 30%) expressed genes
-
low expressed genes, please check the list of low (bottom 30%) expressed genes
An expression pattern of top(30%)/middle(30%)/low(30%) in this subtype against other subtypes is available in a heatmap
For the top enriched genes, if you want to check whether they are
-
significantly differently expressed genes by eBayes lm fit, please check the list of significant genes
Table 5. Get Full Table This table shows top 10 pathways which are significantly enriched in cluster clus3. It displays only significant gene sets satisfying nom.p.val.threshold (-1), fwer.p.val.threshold (-1) , fdr.q.val.threshold (0.25) and the default table is sorted by Normalized Enrichment Score (NES). Further details on NES statistics, please visit The Broad GSEA website.
GeneSet(GS) | Size(#genes) | genes.ES.table | ES | NES | NOM.p.val | FDR.q.val | FWER.p.val | Tag.. | Gene.. | Signal | FDR..median. | glob.p.val |
---|---|---|---|---|---|---|---|---|---|---|---|---|
BIOCARTA ERK PATHWAY | 27 | genes.ES.table | 0.24 | 0.78 | 0.72 | 1 | 1 | 0.22 | 0.13 | 0.19 | 1 | 0.63 |
BIOCARTA RACCYCD PATHWAY | 25 | genes.ES.table | 0.28 | 1.2 | 0.27 | 0.9 | 1 | 0.2 | 0.1 | 0.18 | 0.82 | 0.35 |
BIOCARTA PYK2 PATHWAY | 27 | genes.ES.table | 0.23 | 0.89 | 0.58 | 0.98 | 1 | 0.3 | 0.16 | 0.25 | 0.97 | 0.45 |
BIOCARTA IL1R PATHWAY | 32 | genes.ES.table | 0.34 | 0.97 | 0.51 | 0.97 | 1 | 0.34 | 0.16 | 0.29 | 0.93 | 0.42 |
KEGG PYRIMIDINE METABOLISM | 96 | genes.ES.table | 0.21 | 0.94 | 0.5 | 0.98 | 1 | 0.083 | 0.056 | 0.079 | 0.95 | 0.44 |
KEGG AMINO SUGAR AND NUCLEOTIDE SUGAR METABOLISM | 42 | genes.ES.table | 0.21 | 0.83 | 0.68 | 1 | 1 | 0.21 | 0.1 | 0.19 | 1 | 0.59 |
KEGG ETHER LIPID METABOLISM | 28 | genes.ES.table | 0.32 | 0.94 | 0.54 | 0.93 | 1 | 0.18 | 0.044 | 0.17 | 0.9 | 0.36 |
KEGG AMINOACYL TRNA BIOSYNTHESIS | 40 | genes.ES.table | 0.35 | 1.4 | 0.13 | 0.74 | 0.97 | 1 | 0.65 | 0.35 | 0.57 | 0.26 |
KEGG RIBOSOME | 85 | genes.ES.table | 0.61 | 1.5 | 0.12 | 0.53 | 0.87 | 0.88 | 0.26 | 0.65 | 0.35 | 0.15 |
KEGG RNA POLYMERASE | 28 | genes.ES.table | 0.28 | 0.99 | 0.42 | 1 | 1 | 0.93 | 0.49 | 0.48 | 0.97 | 0.46 |
Table S21. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | PSMB8 | PSMB8 | PSMB8 | 2091 | 0.076 | -0.068 | YES |
2 | PSMD6 | PSMD6 | PSMD6 | 2137 | 0.074 | -0.025 | YES |
3 | PSME2 | PSME2 | PSME2 | 2249 | 0.07 | 0.011 | YES |
4 | PSMA5 | PSMA5 | PSMA5 | 2268 | 0.069 | 0.053 | YES |
5 | PSMB7 | PSMB7 | PSMB7 | 2306 | 0.068 | 0.092 | YES |
6 | POMP | POMP | POMP | 2338 | 0.067 | 0.13 | YES |
7 | PSMD13 | PSMD13 | PSMD13 | 2351 | 0.067 | 0.17 | YES |
8 | PSMA6 | PSMA6 | PSMA6 | 2698 | 0.056 | 0.19 | YES |
9 | PSMB5 | PSMB5 | PSMB5 | 2859 | 0.051 | 0.21 | YES |
10 | PSMC1 | PSMC1 | PSMC1 | 2867 | 0.051 | 0.24 | YES |
11 | PSMB2 | PSMB2 | PSMB2 | 3058 | 0.047 | 0.26 | YES |
12 | PSMB6 | PSMB6 | PSMB6 | 3119 | 0.045 | 0.28 | YES |
13 | PSME1 | PSME1 | PSME1 | 3153 | 0.044 | 0.31 | YES |
14 | PSMD12 | PSMD12 | PSMD12 | 3309 | 0.042 | 0.32 | YES |
15 | PSMA4 | PSMA4 | PSMA4 | 3381 | 0.04 | 0.34 | YES |
16 | PSMA3 | PSMA3 | PSMA3 | 3420 | 0.039 | 0.37 | YES |
17 | PSMB9 | PSMB9 | PSMB9 | 3422 | 0.039 | 0.39 | YES |
18 | PSMC3 | PSMC3 | PSMC3 | 3447 | 0.039 | 0.41 | YES |
19 | PSMA2 | PSMA2 | PSMA2 | 3506 | 0.037 | 0.43 | YES |
20 | PSMD7 | PSMD7 | PSMD7 | 3556 | 0.037 | 0.45 | YES |
21 | PSMC6 | PSMC6 | PSMC6 | 3690 | 0.034 | 0.47 | YES |
22 | PSMA7 | PSMA7 | PSMA7 | 3708 | 0.034 | 0.48 | YES |
23 | PSMB1 | PSMB1 | PSMB1 | 3744 | 0.033 | 0.5 | YES |
24 | PSMD8 | PSMD8 | PSMD8 | 3803 | 0.032 | 0.52 | YES |
25 | PSME3 | PSME3 | PSME3 | 3808 | 0.032 | 0.54 | YES |
26 | PSMD1 | PSMD1 | PSMD1 | 3936 | 0.03 | 0.55 | YES |
27 | PSMB3 | PSMB3 | PSMB3 | 3967 | 0.029 | 0.57 | YES |
28 | PSMF1 | PSMF1 | PSMF1 | 4007 | 0.029 | 0.58 | YES |
29 | PSMA1 | PSMA1 | PSMA1 | 4123 | 0.027 | 0.59 | YES |
30 | PSMC5 | PSMC5 | PSMC5 | 4556 | 0.019 | 0.58 | YES |
31 | SHFM1 | SHFM1 | SHFM1 | 4790 | 0.016 | 0.58 | YES |
32 | PSMD3 | PSMD3 | PSMD3 | 4845 | 0.015 | 0.58 | YES |
33 | PSMD14 | PSMD14 | PSMD14 | 4877 | 0.014 | 0.59 | YES |
34 | IFNG | IFNG | IFNG | 4938 | 0.013 | 0.6 | YES |
35 | PSMD4 | PSMD4 | PSMD4 | 5211 | 0.0085 | 0.58 | NO |
36 | PSMC4 | PSMC4 | PSMC4 | 5296 | 0.0071 | 0.58 | NO |
37 | PSMD11 | PSMD11 | PSMD11 | 5850 | -0.0014 | 0.56 | NO |
38 | PSMB4 | PSMB4 | PSMB4 | 6201 | -0.0067 | 0.54 | NO |
39 | PSMC2 | PSMC2 | PSMC2 | 6286 | -0.0078 | 0.54 | NO |
40 | PSMD2 | PSMD2 | PSMD2 | 6862 | -0.016 | 0.52 | NO |
41 | PSME4 | PSME4 | PSME4 | 7400 | -0.025 | 0.5 | NO |
42 | PSMA8 | PSMA8 | PSMA8 | 13189 | -0.15 | 0.28 | NO |
Figure S41. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA ERK PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.BIOCARTA_ERK_PATHWAY.heatmap_plot.clus3.20.png)
Figure S42. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA ERK PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_PROTEASOME.volcano_plot.clus3.1.png)
Table S22. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | IL1B | IL1B | IL1B | 320 | 0.31 | 0.074 | YES |
2 | IL18 | IL18 | IL18 | 334 | 0.31 | 0.16 | YES |
3 | POLR3G | POLR3G | POLR3G | 394 | 0.28 | 0.24 | YES |
4 | AIM2 | AIM2 | AIM2 | 604 | 0.22 | 0.29 | YES |
5 | CASP1 | CASP1 | CASP1 | 698 | 0.19 | 0.35 | YES |
6 | IRF7 | IRF7 | IRF7 | 711 | 0.19 | 0.4 | YES |
7 | PYCARD | PYCARD | PYCARD | 720 | 0.19 | 0.46 | YES |
8 | IFNA13 | IFNA13 | IFNA13 | 806 | 0.18 | 0.5 | YES |
9 | DDX58 | DDX58 | DDX58 | 890 | 0.16 | 0.55 | YES |
10 | IFNA1 | IFNA1 | IFNA1 | 975 | 0.15 | 0.59 | YES |
11 | CXCL10 | CXCL10 | CXCL10 | 1283 | 0.12 | 0.61 | YES |
12 | ZBP1 | ZBP1 | ZBP1 | 1680 | 0.094 | 0.61 | YES |
13 | RIPK3 | RIPK3 | RIPK3 | 2846 | 0.052 | 0.56 | NO |
14 | CCL5 | CCL5 | CCL5 | 2944 | 0.05 | 0.57 | NO |
15 | POLR1D | POLR1D | POLR1D | 3240 | 0.043 | 0.57 | NO |
16 | POLR3B | POLR3B | POLR3B | 3363 | 0.04 | 0.57 | NO |
17 | POLR3H | POLR3H | POLR3H | 3557 | 0.036 | 0.57 | NO |
18 | RELA | RELA | RELA | 4106 | 0.027 | 0.55 | NO |
19 | IRF3 | IRF3 | IRF3 | 4414 | 0.022 | 0.54 | NO |
20 | NFKBIB | NFKBIB | NFKBIB | 4421 | 0.021 | 0.55 | NO |
21 | POLR3A | POLR3A | POLR3A | 4660 | 0.018 | 0.54 | NO |
22 | RIPK1 | RIPK1 | RIPK1 | 5029 | 0.012 | 0.52 | NO |
23 | POLR1C | POLR1C | POLR1C | 5148 | 0.0095 | 0.52 | NO |
24 | IFNB1 | IFNB1 | IFNB1 | 5165 | 0.0092 | 0.52 | NO |
25 | NFKBIA | NFKBIA | NFKBIA | 5419 | 0.0052 | 0.51 | NO |
26 | NFKB1 | NFKB1 | NFKB1 | 5544 | 0.0032 | 0.5 | NO |
27 | IL33 | IL33 | IL33 | 5889 | -0.0019 | 0.49 | NO |
28 | IKBKG | IKBKG | IKBKG | 5923 | -0.0025 | 0.48 | NO |
29 | CHUK | CHUK | CHUK | 6098 | -0.0053 | 0.48 | NO |
30 | TMEM173 | TMEM173 | TMEM173 | 6311 | -0.0082 | 0.47 | NO |
31 | ADAR | ADAR | ADAR | 6661 | -0.013 | 0.45 | NO |
32 | POLR3F | POLR3F | POLR3F | 6687 | -0.014 | 0.46 | NO |
33 | CCL4 | CCL4 | CCL4 | 7340 | -0.024 | 0.43 | NO |
34 | TREX1 | TREX1 | TREX1 | 7833 | -0.031 | 0.41 | NO |
35 | IL6 | IL6 | IL6 | 7998 | -0.034 | 0.41 | NO |
36 | POLR3K | POLR3K | POLR3K | 8109 | -0.036 | 0.41 | NO |
37 | MAVS | MAVS | MAVS | 8501 | -0.042 | 0.4 | NO |
38 | TBK1 | TBK1 | TBK1 | 8519 | -0.042 | 0.42 | NO |
39 | IKBKE | IKBKE | IKBKE | 8737 | -0.046 | 0.42 | NO |
40 | POLR3C | POLR3C | POLR3C | 8821 | -0.047 | 0.43 | NO |
41 | POLR3D | POLR3D | POLR3D | 8956 | -0.05 | 0.43 | NO |
42 | IKBKB | IKBKB | IKBKB | 9933 | -0.068 | 0.4 | NO |
43 | CCL4L2 | CCL4L2 | CCL4L2 | 10514 | -0.08 | 0.39 | NO |
44 | POLR3GL | POLR3GL | POLR3GL | 12195 | -0.12 | 0.34 | NO |
Figure S43. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA RACCYCD PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.BIOCARTA_RACCYCD_PATHWAY.heatmap_plot.clus3.9.png)
Figure S44. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA RACCYCD PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_CYTOSOLIC_DNA_SENSING_PATHWAY.volcano_plot.clus3.2.png)
Table S23. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | RPL36A | RPL36A | RPL36A | 1134 | 0.13 | -0.024 | YES |
2 | RPL39 | RPL39 | RPL39 | 1913 | 0.082 | -0.044 | YES |
3 | RPL36AL | RPL36AL | RPL36AL | 1932 | 0.082 | -0.022 | YES |
4 | RPL9 | RPL9 | RPL9 | 2039 | 0.078 | -0.0057 | YES |
5 | RPLP2 | RPLP2 | RPLP2 | 2290 | 0.069 | 0.000021 | YES |
6 | RPL22L1 | RPL22L1 | RPL22L1 | 2305 | 0.068 | 0.018 | YES |
7 | RPS13 | RPS13 | RPS13 | 2372 | 0.066 | 0.034 | YES |
8 | RPL35 | RPL35 | RPL35 | 2380 | 0.066 | 0.052 | YES |
9 | RPL29 | RPL29 | RPL29 | 2397 | 0.065 | 0.069 | YES |
10 | RPL27A | RPL27A | RPL27A | 2447 | 0.064 | 0.084 | YES |
11 | RPL11 | RPL11 | RPL11 | 2579 | 0.059 | 0.094 | YES |
12 | RPL18A | RPL18A | RPL18A | 2585 | 0.059 | 0.11 | YES |
13 | RPS3 | RPS3 | RPS3 | 2621 | 0.058 | 0.12 | YES |
14 | RPL36 | RPL36 | RPL36 | 2653 | 0.057 | 0.14 | YES |
15 | RPL38 | RPL38 | RPL38 | 2684 | 0.056 | 0.15 | YES |
16 | RPL26L1 | RPL26L1 | RPL26L1 | 2807 | 0.053 | 0.16 | YES |
17 | RPS10 | RPS10 | RPS10 | 2812 | 0.052 | 0.18 | YES |
18 | RPS12 | RPS12 | RPS12 | 2833 | 0.052 | 0.19 | YES |
19 | RPS29 | RPS29 | RPS29 | 2848 | 0.052 | 0.2 | YES |
20 | RPS2 | RPS2 | RPS2 | 2902 | 0.05 | 0.21 | YES |
21 | MRPL13 | MRPL13 | MRPL13 | 2930 | 0.05 | 0.23 | YES |
22 | RPS15 | RPS15 | RPS15 | 2938 | 0.05 | 0.24 | YES |
23 | RPL34 | RPL34 | RPL34 | 2949 | 0.049 | 0.25 | YES |
24 | RPL27 | RPL27 | RPL27 | 2964 | 0.049 | 0.27 | YES |
25 | RPS16 | RPS16 | RPS16 | 2976 | 0.049 | 0.28 | YES |
26 | RPL14 | RPL14 | RPL14 | 3068 | 0.047 | 0.29 | YES |
27 | RPS25 | RPS25 | RPS25 | 3127 | 0.045 | 0.3 | YES |
28 | FAU | FAU | FAU | 3152 | 0.044 | 0.31 | YES |
29 | RPSA | RPSA | RPSA | 3205 | 0.044 | 0.32 | YES |
30 | RPL8 | RPL8 | RPL8 | 3241 | 0.043 | 0.33 | YES |
31 | UBA52 | UBA52 | UBA52 | 3275 | 0.042 | 0.34 | YES |
32 | RPS24 | RPS24 | RPS24 | 3290 | 0.042 | 0.35 | YES |
33 | RPL13 | RPL13 | RPL13 | 3305 | 0.042 | 0.36 | YES |
34 | RPS21 | RPS21 | RPS21 | 3336 | 0.041 | 0.37 | YES |
35 | RPL28 | RPL28 | RPL28 | 3408 | 0.04 | 0.38 | YES |
36 | RPL4 | RPL4 | RPL4 | 3483 | 0.038 | 0.38 | YES |
37 | RPS26 | RPS26 | RPS26 | 3600 | 0.036 | 0.39 | YES |
38 | RPS18 | RPS18 | RPS18 | 3602 | 0.036 | 0.4 | YES |
39 | RPL26 | RPL26 | RPL26 | 3615 | 0.036 | 0.41 | YES |
40 | RPS5 | RPS5 | RPS5 | 3719 | 0.034 | 0.41 | YES |
41 | RPL7A | RPL7A | RPL7A | 3729 | 0.034 | 0.42 | YES |
42 | RPS15A | RPS15A | RPS15A | 3741 | 0.033 | 0.43 | YES |
43 | RPS17 | RPS17 | RPS17 | 3759 | 0.033 | 0.44 | YES |
44 | RPL37A | RPL37A | RPL37A | 3787 | 0.033 | 0.44 | YES |
45 | RPL32 | RPL32 | RPL32 | 3823 | 0.032 | 0.45 | YES |
46 | RPL3L | RPL3L | RPL3L | 3839 | 0.032 | 0.46 | YES |
47 | RPS9 | RPS9 | RPS9 | 3843 | 0.032 | 0.47 | YES |
48 | RPL12 | RPL12 | RPL12 | 3851 | 0.032 | 0.48 | YES |
49 | RPL23A | RPL23A | RPL23A | 3867 | 0.031 | 0.48 | YES |
50 | RPS8 | RPS8 | RPS8 | 3894 | 0.031 | 0.49 | YES |
51 | RPL6 | RPL6 | RPL6 | 3895 | 0.031 | 0.5 | YES |
52 | RPS6 | RPS6 | RPS6 | 3916 | 0.03 | 0.51 | YES |
53 | RPL5 | RPL5 | RPL5 | 3929 | 0.03 | 0.52 | YES |
54 | RPL7 | RPL7 | RPL7 | 3937 | 0.03 | 0.52 | YES |
55 | RPL37 | RPL37 | RPL37 | 3948 | 0.03 | 0.53 | YES |
56 | RPL24 | RPL24 | RPL24 | 4023 | 0.028 | 0.54 | YES |
57 | RPS7 | RPS7 | RPS7 | 4087 | 0.027 | 0.54 | YES |
58 | RPS23 | RPS23 | RPS23 | 4103 | 0.027 | 0.55 | YES |
59 | RPL17 | RPL17 | RPL17 | 4122 | 0.027 | 0.55 | YES |
60 | RPL30 | RPL30 | RPL30 | 4198 | 0.026 | 0.56 | YES |
61 | RPS20 | RPS20 | RPS20 | 4242 | 0.025 | 0.56 | YES |
62 | RPS27A | RPS27A | RPS27A | 4244 | 0.025 | 0.57 | YES |
63 | RPL19 | RPL19 | RPL19 | 4285 | 0.024 | 0.57 | YES |
64 | RPL13A | RPL13A | RPL13A | 4302 | 0.024 | 0.58 | YES |
65 | RPS4X | RPS4X | RPS4X | 4303 | 0.024 | 0.58 | YES |
66 | RPL10A | RPL10A | RPL10A | 4340 | 0.023 | 0.59 | YES |
67 | RPS3A | RPS3A | RPS3A | 4366 | 0.022 | 0.59 | YES |
68 | RPL31 | RPL31 | RPL31 | 4430 | 0.021 | 0.6 | YES |
69 | RPS19 | RPS19 | RPS19 | 4468 | 0.021 | 0.6 | YES |
70 | RPL22 | RPL22 | RPL22 | 4519 | 0.02 | 0.6 | YES |
71 | RPL10 | RPL10 | RPL10 | 4531 | 0.02 | 0.61 | YES |
72 | RPL23 | RPL23 | RPL23 | 4607 | 0.018 | 0.61 | YES |
73 | RPS11 | RPS11 | RPS11 | 4759 | 0.016 | 0.6 | YES |
74 | RSL24D1 | RSL24D1 | RSL24D1 | 4796 | 0.016 | 0.61 | YES |
75 | RPLP0 | RPLP0 | RPLP0 | 4835 | 0.015 | 0.61 | YES |
76 | RPLP1 | RPLP1 | RPLP1 | 5277 | 0.0074 | 0.59 | NO |
77 | RPL3 | RPL3 | RPL3 | 5375 | 0.0059 | 0.58 | NO |
78 | RPS28 | RPS28 | RPS28 | 5537 | 0.0034 | 0.58 | NO |
79 | RPL41 | RPL41 | RPL41 | 5808 | -0.00081 | 0.56 | NO |
80 | RPL15 | RPL15 | RPL15 | 6554 | -0.012 | 0.52 | NO |
81 | RPS27L | RPS27L | RPS27L | 6664 | -0.013 | 0.52 | NO |
82 | RPL35A | RPL35A | RPL35A | 6830 | -0.016 | 0.52 | NO |
83 | RPS27 | RPS27 | RPS27 | 8072 | -0.035 | 0.46 | NO |
84 | RPS4Y1 | RPS4Y1 | RPS4Y1 | 12636 | -0.14 | 0.25 | NO |
85 | RPL21 | RPL21 | RPL21 | 15213 | -0.23 | 0.17 | NO |
Figure S45. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA PYK2 PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.BIOCARTA_PYK2_PATHWAY.heatmap_plot.clus3.17.png)
Figure S46. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA PYK2 PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_RIBOSOME.volcano_plot.clus3.3.png)
Table S24. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | WARS | WARS | WARS | 1138 | 0.13 | 0.017 | YES |
2 | RARS | RARS | RARS | 1598 | 0.098 | 0.05 | YES |
3 | GARS | GARS | GARS | 2141 | 0.074 | 0.064 | YES |
4 | FARS2 | FARS2 | FARS2 | 2267 | 0.069 | 0.098 | YES |
5 | SARS | SARS | SARS | 2272 | 0.069 | 0.14 | YES |
6 | HARS | HARS | HARS | 2471 | 0.063 | 0.16 | YES |
7 | CARS2 | CARS2 | CARS2 | 2568 | 0.06 | 0.19 | YES |
8 | VARS | VARS | VARS | 2642 | 0.057 | 0.22 | YES |
9 | KARS | KARS | KARS | 2710 | 0.056 | 0.25 | YES |
10 | YARS | YARS | YARS | 2721 | 0.055 | 0.28 | YES |
11 | FARSB | FARSB | FARSB | 3095 | 0.046 | 0.29 | YES |
12 | CARS | CARS | CARS | 3861 | 0.031 | 0.27 | YES |
13 | TARS | TARS | TARS | 4128 | 0.027 | 0.27 | YES |
14 | YARS2 | YARS2 | YARS2 | 4532 | 0.02 | 0.26 | YES |
15 | FARSA | FARSA | FARSA | 4795 | 0.016 | 0.25 | YES |
16 | QARS | QARS | QARS | 4918 | 0.014 | 0.26 | YES |
17 | HARS2 | HARS2 | HARS2 | 5275 | 0.0074 | 0.24 | YES |
18 | TARSL2 | TARSL2 | TARSL2 | 5346 | 0.0064 | 0.24 | YES |
19 | SARS2 | SARS2 | SARS2 | 5638 | 0.0019 | 0.22 | YES |
20 | DARS | DARS | DARS | 6345 | -0.0086 | 0.19 | YES |
21 | LARS2 | LARS2 | LARS2 | 6473 | -0.01 | 0.19 | YES |
22 | NARS2 | NARS2 | NARS2 | 6550 | -0.012 | 0.19 | YES |
23 | MARS | MARS | MARS | 6821 | -0.016 | 0.19 | YES |
24 | EARS2 | EARS2 | EARS2 | 6851 | -0.016 | 0.2 | YES |
25 | IARS2 | IARS2 | IARS2 | 6904 | -0.017 | 0.2 | YES |
26 | NARS | NARS | NARS | 6928 | -0.017 | 0.21 | YES |
27 | EPRS | EPRS | EPRS | 7039 | -0.019 | 0.22 | YES |
28 | RARS2 | RARS2 | RARS2 | 7073 | -0.02 | 0.23 | YES |
29 | IARS | IARS | IARS | 7219 | -0.022 | 0.23 | YES |
30 | MTFMT | MTFMT | MTFMT | 7302 | -0.023 | 0.24 | YES |
31 | LARS | LARS | LARS | 7523 | -0.026 | 0.24 | YES |
32 | WARS2 | WARS2 | WARS2 | 7756 | -0.03 | 0.25 | YES |
33 | AARS2 | AARS2 | AARS2 | 8380 | -0.04 | 0.24 | YES |
34 | VARS2 | VARS2 | VARS2 | 8470 | -0.041 | 0.26 | YES |
35 | AARS | AARS | AARS | 8787 | -0.047 | 0.27 | YES |
36 | DARS2 | DARS2 | DARS2 | 9110 | -0.053 | 0.28 | YES |
37 | PSTK | PSTK | PSTK | 10217 | -0.074 | 0.26 | YES |
38 | SEPSECS | SEPSECS | SEPSECS | 10868 | -0.087 | 0.28 | YES |
39 | PARS2 | PARS2 | PARS2 | 11445 | -0.1 | 0.31 | YES |
40 | MARS2 | MARS2 | MARS2 | 11996 | -0.12 | 0.35 | YES |
Figure S47. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: BIOCARTA IL1R PATHWAY.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.BIOCARTA_IL1R_PATHWAY.heatmap_plot.clus3.14.png)
Figure S48. Get High-res Image For the top 5 core enriched genes in the pathway: BIOCARTA IL1R PATHWAY, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_AMINOACYL_TRNA_BIOSYNTHESIS.volcano_plot.clus3.4.png)
Table S25. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | HSPA2 | HSPA2 | HSPA2 | 2097 | 0.076 | -0.092 | YES |
2 | EIF4A3 | EIF4A3 | EIF4A3 | 2244 | 0.07 | -0.078 | YES |
3 | CWC15 | CWC15 | CWC15 | 3077 | 0.047 | -0.11 | YES |
4 | CCDC12 | CCDC12 | CCDC12 | 3134 | 0.045 | -0.098 | YES |
5 | HSPA8 | HSPA8 | HSPA8 | 3166 | 0.044 | -0.086 | YES |
6 | SNRPB2 | SNRPB2 | SNRPB2 | 3183 | 0.044 | -0.074 | YES |
7 | CTNNBL1 | CTNNBL1 | CTNNBL1 | 3203 | 0.044 | -0.061 | YES |
8 | DHX38 | DHX38 | DHX38 | 3328 | 0.041 | -0.055 | YES |
9 | SF3B5 | SF3B5 | SF3B5 | 3329 | 0.041 | -0.042 | YES |
10 | PLRG1 | PLRG1 | PLRG1 | 3541 | 0.037 | -0.043 | YES |
11 | PPIL1 | PPIL1 | PPIL1 | 3618 | 0.036 | -0.036 | YES |
12 | SNRPA1 | SNRPA1 | SNRPA1 | 3665 | 0.034 | -0.028 | YES |
13 | SFRS5 | SFRS5 | SFRS5 | 3745 | 0.033 | -0.022 | YES |
14 | SNRPD1 | SNRPD1 | SNRPD1 | 3819 | 0.032 | -0.016 | YES |
15 | PUF60 | PUF60 | PUF60 | 4003 | 0.029 | -0.017 | YES |
16 | SNRPC | SNRPC | SNRPC | 4097 | 0.027 | -0.014 | YES |
17 | PRPF31 | PRPF31 | PRPF31 | 4134 | 0.026 | -0.0077 | YES |
18 | SNRPF | SNRPF | SNRPF | 4300 | 0.024 | -0.0094 | YES |
19 | TRA2A | TRA2A | TRA2A | 4395 | 0.022 | -0.0078 | YES |
20 | PPIE | PPIE | PPIE | 4416 | 0.022 | -0.0023 | YES |
21 | SF3B14 | SF3B14 | SF3B14 | 4527 | 0.02 | -0.0022 | YES |
22 | SNRPD2 | SNRPD2 | SNRPD2 | 4536 | 0.02 | 0.0034 | YES |
23 | ZMAT2 | ZMAT2 | ZMAT2 | 4610 | 0.018 | 0.0051 | YES |
24 | AQR | AQR | AQR | 4630 | 0.018 | 0.0097 | YES |
25 | SNRPB | SNRPB | SNRPB | 4678 | 0.018 | 0.013 | YES |
26 | PCBP1 | PCBP1 | PCBP1 | 4680 | 0.018 | 0.018 | YES |
27 | DHX15 | DHX15 | DHX15 | 4682 | 0.018 | 0.023 | YES |
28 | LSM7 | LSM7 | LSM7 | 4733 | 0.017 | 0.026 | YES |
29 | SNRPD3 | SNRPD3 | SNRPD3 | 4737 | 0.017 | 0.031 | YES |
30 | LSM5 | LSM5 | LSM5 | 4911 | 0.014 | 0.025 | YES |
31 | BUD31 | BUD31 | BUD31 | 4939 | 0.013 | 0.028 | YES |
32 | PPIH | PPIH | PPIH | 4994 | 0.012 | 0.029 | YES |
33 | HNRNPC | HNRNPC | HNRNPC | 5068 | 0.011 | 0.028 | YES |
34 | SFRS2B | SFRS2B | SFRS2B | 5093 | 0.01 | 0.03 | YES |
35 | SNW1 | SNW1 | SNW1 | 5121 | 0.0099 | 0.032 | YES |
36 | SNRNP70 | SNRNP70 | SNRNP70 | 5137 | 0.0097 | 0.034 | YES |
37 | EFTUD2 | EFTUD2 | EFTUD2 | 5158 | 0.0093 | 0.036 | YES |
38 | MAGOH | MAGOH | MAGOH | 5170 | 0.0091 | 0.038 | YES |
39 | SMNDC1 | SMNDC1 | SMNDC1 | 5251 | 0.0078 | 0.036 | YES |
40 | SNRPG | SNRPG | SNRPG | 5271 | 0.0075 | 0.037 | YES |
41 | SART1 | SART1 | SART1 | 5361 | 0.0061 | 0.034 | YES |
42 | PRPF6 | PRPF6 | PRPF6 | 5379 | 0.0058 | 0.035 | YES |
43 | PRPF38A | PRPF38A | PRPF38A | 5423 | 0.0052 | 0.034 | YES |
44 | SF3A2 | SF3A2 | SF3A2 | 5452 | 0.0048 | 0.034 | YES |
45 | RBM22 | RBM22 | RBM22 | 5507 | 0.0038 | 0.032 | YES |
46 | TXNL4A | TXNL4A | TXNL4A | 5510 | 0.0038 | 0.033 | YES |
47 | LSM3 | LSM3 | LSM3 | 5568 | 0.0029 | 0.031 | YES |
48 | PQBP1 | PQBP1 | PQBP1 | 5592 | 0.0027 | 0.031 | YES |
49 | PRPF19 | PRPF19 | PRPF19 | 5657 | 0.0017 | 0.028 | YES |
50 | BAT1 | BAT1 | BAT1 | 5783 | -0.00032 | 0.021 | YES |
51 | XAB2 | XAB2 | XAB2 | 5805 | -0.00074 | 0.02 | YES |
52 | SNRPA | SNRPA | SNRPA | 5811 | -0.00083 | 0.02 | YES |
53 | SFRS1 | SFRS1 | SFRS1 | 5916 | -0.0024 | 0.015 | YES |
54 | RBM25 | RBM25 | RBM25 | 6006 | -0.004 | 0.011 | YES |
55 | HNRNPM | HNRNPM | HNRNPM | 6063 | -0.0047 | 0.0097 | YES |
56 | PRPF40B | PRPF40B | PRPF40B | 6091 | -0.0052 | 0.0098 | YES |
57 | HNRNPK | HNRNPK | HNRNPK | 6102 | -0.0054 | 0.011 | YES |
58 | SYF2 | SYF2 | SYF2 | 6144 | -0.006 | 0.011 | YES |
59 | U2AF2 | U2AF2 | U2AF2 | 6171 | -0.0064 | 0.011 | YES |
60 | THOC4 | THOC4 | THOC4 | 6186 | -0.0065 | 0.012 | YES |
61 | SF3B2 | SF3B2 | SF3B2 | 6192 | -0.0066 | 0.014 | YES |
62 | NHP2L1 | NHP2L1 | NHP2L1 | 6250 | -0.0073 | 0.013 | YES |
63 | SF3B3 | SF3B3 | SF3B3 | 6258 | -0.0075 | 0.015 | YES |
64 | NCBP1 | NCBP1 | NCBP1 | 6277 | -0.0077 | 0.017 | YES |
65 | SFRS3 | SFRS3 | SFRS3 | 6339 | -0.0085 | 0.016 | YES |
66 | HNRNPA1 | HNRNPA1 | HNRNPA1 | 6347 | -0.0086 | 0.018 | YES |
67 | RBM17 | RBM17 | RBM17 | 6378 | -0.009 | 0.019 | YES |
68 | SFRS9 | SFRS9 | SFRS9 | 6386 | -0.0091 | 0.022 | YES |
69 | U2AF1 | U2AF1 | U2AF1 | 6400 | -0.0093 | 0.024 | YES |
70 | SF3A3 | SF3A3 | SF3A3 | 6496 | -0.011 | 0.022 | YES |
71 | DHX16 | DHX16 | DHX16 | 6516 | -0.011 | 0.024 | YES |
72 | TCERG1 | TCERG1 | TCERG1 | 6520 | -0.011 | 0.027 | YES |
73 | SLU7 | SLU7 | SLU7 | 6540 | -0.011 | 0.03 | YES |
74 | SFRS4 | SFRS4 | SFRS4 | 6544 | -0.011 | 0.033 | YES |
75 | DDX23 | DDX23 | DDX23 | 6645 | -0.013 | 0.032 | YES |
76 | SFRS2 | SFRS2 | SFRS2 | 6653 | -0.013 | 0.035 | YES |
77 | RBM8A | RBM8A | RBM8A | 6716 | -0.014 | 0.036 | YES |
78 | PRPF40A | PRPF40A | PRPF40A | 6765 | -0.015 | 0.038 | YES |
79 | HNRNPA1L2 | HNRNPA1L2 | HNRNPA1L2 | 6841 | -0.016 | 0.039 | YES |
80 | SFRS6 | SFRS6 | SFRS6 | 6873 | -0.016 | 0.042 | YES |
81 | ACIN1 | ACIN1 | ACIN1 | 6878 | -0.016 | 0.047 | YES |
82 | DHX8 | DHX8 | DHX8 | 6986 | -0.018 | 0.047 | YES |
83 | CRNKL1 | CRNKL1 | CRNKL1 | 7011 | -0.019 | 0.051 | YES |
84 | DDX5 | DDX5 | DDX5 | 7071 | -0.02 | 0.054 | YES |
85 | SNRNP40 | SNRNP40 | SNRNP40 | 7161 | -0.021 | 0.056 | YES |
86 | SF3B4 | SF3B4 | SF3B4 | 7205 | -0.022 | 0.06 | YES |
87 | SNRPE | SNRPE | SNRPE | 7231 | -0.022 | 0.065 | YES |
88 | LSM4 | LSM4 | LSM4 | 7258 | -0.022 | 0.071 | YES |
89 | LSM2 | LSM2 | LSM2 | 7329 | -0.024 | 0.074 | YES |
90 | PHF5A | PHF5A | PHF5A | 7436 | -0.025 | 0.076 | YES |
91 | PRPF18 | PRPF18 | PRPF18 | 7470 | -0.026 | 0.082 | YES |
92 | LSM6 | LSM6 | LSM6 | 7491 | -0.026 | 0.089 | YES |
93 | PRPF4 | PRPF4 | PRPF4 | 7529 | -0.026 | 0.095 | YES |
94 | HNRNPU | HNRNPU | HNRNPU | 7691 | -0.029 | 0.096 | YES |
95 | BCAS2 | BCAS2 | BCAS2 | 7702 | -0.029 | 0.1 | YES |
96 | THOC1 | THOC1 | THOC1 | 7789 | -0.03 | 0.11 | YES |
97 | DDX46 | DDX46 | DDX46 | 7817 | -0.031 | 0.12 | YES |
98 | SFRS7 | SFRS7 | SFRS7 | 7928 | -0.033 | 0.12 | YES |
99 | PRPF38B | PRPF38B | PRPF38B | 7967 | -0.034 | 0.13 | YES |
100 | CDC5L | CDC5L | CDC5L | 7989 | -0.034 | 0.14 | YES |
101 | SFRS13A | SFRS13A | SFRS13A | 8164 | -0.037 | 0.14 | YES |
102 | DDX42 | DDX42 | DDX42 | 8176 | -0.037 | 0.15 | YES |
103 | HNRNPA3 | HNRNPA3 | HNRNPA3 | 8257 | -0.038 | 0.16 | YES |
104 | HSPA6 | HSPA6 | HSPA6 | 8337 | -0.039 | 0.17 | YES |
105 | MAGOHB | MAGOHB | MAGOHB | 8349 | -0.039 | 0.18 | YES |
106 | PRPF8 | PRPF8 | PRPF8 | 8432 | -0.04 | 0.19 | YES |
107 | SF3B1 | SF3B1 | SF3B1 | 8477 | -0.041 | 0.2 | YES |
108 | PRPF3 | PRPF3 | PRPF3 | 8482 | -0.041 | 0.21 | YES |
109 | CDC40 | CDC40 | CDC40 | 8551 | -0.042 | 0.22 | YES |
110 | ISY1 | ISY1 | ISY1 | 8560 | -0.043 | 0.23 | YES |
111 | SNRNP27 | SNRNP27 | SNRNP27 | 8621 | -0.044 | 0.24 | YES |
112 | SR140 | SR140 | SR140 | 8626 | -0.044 | 0.25 | YES |
113 | HSPA1B | HSPA1B | HSPA1B | 8738 | -0.046 | 0.26 | YES |
114 | CHERP | CHERP | CHERP | 8742 | -0.046 | 0.28 | YES |
115 | TRA2B | TRA2B | TRA2B | 8743 | -0.046 | 0.29 | YES |
116 | THOC2 | THOC2 | THOC2 | 8769 | -0.046 | 0.3 | YES |
117 | SF3A1 | SF3A1 | SF3A1 | 8773 | -0.047 | 0.32 | YES |
118 | NAA38 | NAA38 | NAA38 | 9038 | -0.051 | 0.32 | YES |
119 | RBMX | RBMX | RBMX | 9062 | -0.052 | 0.33 | YES |
120 | SNRNP200 | SNRNP200 | SNRNP200 | 9180 | -0.054 | 0.34 | YES |
121 | THOC3 | THOC3 | THOC3 | 9405 | -0.058 | 0.35 | YES |
122 | WBP11 | WBP11 | WBP11 | 9915 | -0.067 | 0.34 | YES |
123 | HSPA1A | HSPA1A | HSPA1A | 10033 | -0.07 | 0.36 | YES |
124 | USP39 | USP39 | USP39 | 10601 | -0.082 | 0.35 | NO |
125 | HSPA1L | HSPA1L | HSPA1L | 15344 | -0.24 | 0.16 | NO |
Figure S49. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG PYRIMIDINE METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_PYRIMIDINE_METABOLISM.heatmap_plot.clus3.15.png)
Figure S50. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG PYRIMIDINE METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_SPLICEOSOME.volcano_plot.clus3.5.png)
Table S26. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | IFNE | IFNE | IFNE | 261 | 0.35 | 0.055 | YES |
2 | IFNK | IFNK | IFNK | 371 | 0.29 | 0.11 | YES |
3 | ISG15 | ISG15 | ISG15 | 488 | 0.25 | 0.15 | YES |
4 | NLRX1 | NLRX1 | NLRX1 | 599 | 0.22 | 0.19 | YES |
5 | IRF7 | IRF7 | IRF7 | 711 | 0.19 | 0.22 | YES |
6 | IFNA13 | IFNA13 | IFNA13 | 806 | 0.18 | 0.25 | YES |
7 | DDX58 | DDX58 | DDX58 | 890 | 0.16 | 0.28 | YES |
8 | IL8 | IL8 | IL8 | 971 | 0.15 | 0.3 | YES |
9 | IFNA1 | IFNA1 | IFNA1 | 975 | 0.15 | 0.33 | YES |
10 | IFIH1 | IFIH1 | IFIH1 | 1024 | 0.15 | 0.36 | YES |
11 | MAPK13 | MAPK13 | MAPK13 | 1233 | 0.13 | 0.37 | YES |
12 | CXCL10 | CXCL10 | CXCL10 | 1283 | 0.12 | 0.39 | YES |
13 | TRADD | TRADD | TRADD | 1625 | 0.097 | 0.39 | YES |
14 | CASP10 | CASP10 | CASP10 | 1802 | 0.087 | 0.4 | YES |
15 | DHX58 | DHX58 | DHX58 | 2646 | 0.057 | 0.37 | NO |
16 | MAPK14 | MAPK14 | MAPK14 | 2929 | 0.05 | 0.36 | NO |
17 | TNF | TNF | TNF | 2951 | 0.049 | 0.37 | NO |
18 | DDX3X | DDX3X | DDX3X | 3206 | 0.044 | 0.36 | NO |
19 | ATG12 | ATG12 | ATG12 | 3392 | 0.04 | 0.36 | NO |
20 | CYLD | CYLD | CYLD | 3558 | 0.036 | 0.36 | NO |
21 | FADD | FADD | FADD | 3795 | 0.032 | 0.35 | NO |
22 | RELA | RELA | RELA | 4106 | 0.027 | 0.34 | NO |
23 | CASP8 | CASP8 | CASP8 | 4141 | 0.026 | 0.34 | NO |
24 | TRIM25 | TRIM25 | TRIM25 | 4322 | 0.023 | 0.34 | NO |
25 | IRF3 | IRF3 | IRF3 | 4414 | 0.022 | 0.34 | NO |
26 | NFKBIB | NFKBIB | NFKBIB | 4421 | 0.021 | 0.34 | NO |
27 | AZI2 | AZI2 | AZI2 | 4469 | 0.021 | 0.34 | NO |
28 | RIPK1 | RIPK1 | RIPK1 | 5029 | 0.012 | 0.32 | NO |
29 | IFNB1 | IFNB1 | IFNB1 | 5165 | 0.0092 | 0.31 | NO |
30 | MAPK11 | MAPK11 | MAPK11 | 5322 | 0.0068 | 0.3 | NO |
31 | NFKBIA | NFKBIA | NFKBIA | 5419 | 0.0052 | 0.3 | NO |
32 | NFKB1 | NFKB1 | NFKB1 | 5544 | 0.0032 | 0.29 | NO |
33 | MAP3K7 | MAP3K7 | MAP3K7 | 5861 | -0.0015 | 0.28 | NO |
34 | IKBKG | IKBKG | IKBKG | 5923 | -0.0025 | 0.27 | NO |
35 | CHUK | CHUK | CHUK | 6098 | -0.0053 | 0.26 | NO |
36 | PIN1 | PIN1 | PIN1 | 6172 | -0.0064 | 0.26 | NO |
37 | TMEM173 | TMEM173 | TMEM173 | 6311 | -0.0082 | 0.26 | NO |
38 | TANK | TANK | TANK | 6314 | -0.0083 | 0.26 | NO |
39 | SIKE1 | SIKE1 | SIKE1 | 6374 | -0.0089 | 0.26 | NO |
40 | TRAF3 | TRAF3 | TRAF3 | 6810 | -0.015 | 0.24 | NO |
41 | MAP3K1 | MAP3K1 | MAP3K1 | 6843 | -0.016 | 0.24 | NO |
42 | OTUD5 | OTUD5 | OTUD5 | 6947 | -0.018 | 0.24 | NO |
43 | ATG5 | ATG5 | ATG5 | 7296 | -0.023 | 0.22 | NO |
44 | DAK | DAK | DAK | 7677 | -0.029 | 0.2 | NO |
45 | MAPK9 | MAPK9 | MAPK9 | 7940 | -0.033 | 0.2 | NO |
46 | TRAF6 | TRAF6 | TRAF6 | 8104 | -0.036 | 0.2 | NO |
47 | MAVS | MAVS | MAVS | 8501 | -0.042 | 0.18 | NO |
48 | TBK1 | TBK1 | TBK1 | 8519 | -0.042 | 0.19 | NO |
49 | IKBKE | IKBKE | IKBKE | 8737 | -0.046 | 0.19 | NO |
50 | TRAF2 | TRAF2 | TRAF2 | 9324 | -0.056 | 0.17 | NO |
51 | IKBKB | IKBKB | IKBKB | 9933 | -0.068 | 0.15 | NO |
52 | IL12B | IL12B | IL12B | 10880 | -0.087 | 0.11 | NO |
53 | TBKBP1 | TBKBP1 | TBKBP1 | 11247 | -0.096 | 0.11 | NO |
54 | MAPK8 | MAPK8 | MAPK8 | 13037 | -0.15 | 0.042 | NO |
55 | MAPK12 | MAPK12 | MAPK12 | 13833 | -0.18 | 0.033 | NO |
56 | DDX3Y | DDX3Y | DDX3Y | 13869 | -0.18 | 0.066 | NO |
57 | MAPK10 | MAPK10 | MAPK10 | 15010 | -0.22 | 0.048 | NO |
58 | RNF125 | RNF125 | RNF125 | 15876 | -0.27 | 0.054 | NO |
59 | IL12A | IL12A | IL12A | 17596 | -0.41 | 0.042 | NO |
Figure S51. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG AMINO SUGAR AND NUCLEOTIDE SUGAR METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_AMINO_SUGAR_AND_NUCLEOTIDE_SUGAR_METABOLISM.heatmap_plot.clus3.18.png)
Figure S52. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG AMINO SUGAR AND NUCLEOTIDE SUGAR METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_RIG_I_LIKE_RECEPTOR_SIGNALING_PATHWAY.volcano_plot.clus3.6.png)
Table S27. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CARD18 | CARD18 | CARD18 | 47 | 0.57 | 0.084 | YES |
2 | MEFV | MEFV | MEFV | 95 | 0.49 | 0.16 | YES |
3 | NOD2 | NOD2 | NOD2 | 248 | 0.36 | 0.2 | YES |
4 | IL1B | IL1B | IL1B | 320 | 0.31 | 0.25 | YES |
5 | IL18 | IL18 | IL18 | 334 | 0.31 | 0.29 | YES |
6 | CARD6 | CARD6 | CARD6 | 518 | 0.24 | 0.32 | YES |
7 | CASP1 | CASP1 | CASP1 | 698 | 0.19 | 0.34 | YES |
8 | PYCARD | PYCARD | PYCARD | 720 | 0.19 | 0.37 | YES |
9 | CASP5 | CASP5 | CASP5 | 730 | 0.19 | 0.4 | YES |
10 | CXCL1 | CXCL1 | CXCL1 | 744 | 0.19 | 0.42 | YES |
11 | IL8 | IL8 | IL8 | 971 | 0.15 | 0.44 | YES |
12 | MAPK13 | MAPK13 | MAPK13 | 1233 | 0.13 | 0.44 | YES |
13 | SUGT1 | SUGT1 | SUGT1 | 1305 | 0.12 | 0.46 | YES |
14 | MAPK3 | MAPK3 | MAPK3 | 1411 | 0.11 | 0.47 | YES |
15 | CXCL2 | CXCL2 | CXCL2 | 2131 | 0.074 | 0.44 | NO |
16 | CCL7 | CCL7 | CCL7 | 2686 | 0.056 | 0.42 | NO |
17 | NLRP1 | NLRP1 | NLRP1 | 2839 | 0.052 | 0.42 | NO |
18 | MAPK14 | MAPK14 | MAPK14 | 2929 | 0.05 | 0.42 | NO |
19 | CCL5 | CCL5 | CCL5 | 2944 | 0.05 | 0.43 | NO |
20 | TNF | TNF | TNF | 2951 | 0.049 | 0.43 | NO |
21 | TRIP6 | TRIP6 | TRIP6 | 3115 | 0.046 | 0.43 | NO |
22 | RIPK2 | RIPK2 | RIPK2 | 3417 | 0.039 | 0.42 | NO |
23 | CCL8 | CCL8 | CCL8 | 3699 | 0.034 | 0.41 | NO |
24 | RELA | RELA | RELA | 4106 | 0.027 | 0.39 | NO |
25 | TAB2 | TAB2 | TAB2 | 4107 | 0.027 | 0.4 | NO |
26 | CASP8 | CASP8 | CASP8 | 4141 | 0.026 | 0.4 | NO |
27 | NLRP3 | NLRP3 | NLRP3 | 4299 | 0.024 | 0.39 | NO |
28 | NOD1 | NOD1 | NOD1 | 4362 | 0.022 | 0.39 | NO |
29 | NFKBIB | NFKBIB | NFKBIB | 4421 | 0.021 | 0.4 | NO |
30 | MAPK11 | MAPK11 | MAPK11 | 5322 | 0.0068 | 0.35 | NO |
31 | NFKBIA | NFKBIA | NFKBIA | 5419 | 0.0052 | 0.34 | NO |
32 | NFKB1 | NFKB1 | NFKB1 | 5544 | 0.0032 | 0.34 | NO |
33 | MAP3K7 | MAP3K7 | MAP3K7 | 5861 | -0.0015 | 0.32 | NO |
34 | TNFAIP3 | TNFAIP3 | TNFAIP3 | 5887 | -0.0019 | 0.32 | NO |
35 | IKBKG | IKBKG | IKBKG | 5923 | -0.0025 | 0.32 | NO |
36 | HSP90AA1 | HSP90AA1 | HSP90AA1 | 6046 | -0.0046 | 0.31 | NO |
37 | BIRC2 | BIRC2 | BIRC2 | 6080 | -0.005 | 0.31 | NO |
38 | CHUK | CHUK | CHUK | 6098 | -0.0053 | 0.31 | NO |
39 | TAB3 | TAB3 | TAB3 | 6315 | -0.0083 | 0.3 | NO |
40 | MAPK1 | MAPK1 | MAPK1 | 6457 | -0.01 | 0.29 | NO |
41 | CCL13 | CCL13 | CCL13 | 6614 | -0.012 | 0.29 | NO |
42 | HSP90B1 | HSP90B1 | HSP90B1 | 7363 | -0.024 | 0.25 | NO |
43 | XIAP | XIAP | XIAP | 7743 | -0.03 | 0.23 | NO |
44 | MAPK9 | MAPK9 | MAPK9 | 7940 | -0.033 | 0.23 | NO |
45 | IL6 | IL6 | IL6 | 7998 | -0.034 | 0.23 | NO |
46 | TRAF6 | TRAF6 | TRAF6 | 8104 | -0.036 | 0.23 | NO |
47 | CCL11 | CCL11 | CCL11 | 9597 | -0.061 | 0.16 | NO |
48 | ERBB2IP | ERBB2IP | ERBB2IP | 9830 | -0.066 | 0.15 | NO |
49 | IKBKB | IKBKB | IKBKB | 9933 | -0.068 | 0.16 | NO |
50 | TAB1 | TAB1 | TAB1 | 9935 | -0.068 | 0.17 | NO |
51 | PSTPIP1 | PSTPIP1 | PSTPIP1 | 10323 | -0.076 | 0.16 | NO |
52 | PYDC1 | PYDC1 | PYDC1 | 12956 | -0.15 | 0.038 | NO |
53 | CCL2 | CCL2 | CCL2 | 13025 | -0.15 | 0.056 | NO |
54 | MAPK8 | MAPK8 | MAPK8 | 13037 | -0.15 | 0.079 | NO |
55 | CARD8 | CARD8 | CARD8 | 13333 | -0.16 | 0.087 | NO |
56 | BIRC3 | BIRC3 | BIRC3 | 13459 | -0.16 | 0.1 | NO |
57 | MAPK12 | MAPK12 | MAPK12 | 13833 | -0.18 | 0.11 | NO |
58 | CARD9 | CARD9 | CARD9 | 14318 | -0.19 | 0.11 | NO |
59 | NLRC4 | NLRC4 | NLRC4 | 14422 | -0.2 | 0.14 | NO |
60 | MAPK10 | MAPK10 | MAPK10 | 15010 | -0.22 | 0.14 | NO |
61 | NAIP | NAIP | NAIP | 15935 | -0.27 | 0.13 | NO |
Figure S53. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG ETHER LIPID METABOLISM.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_ETHER_LIPID_METABOLISM.heatmap_plot.clus3.16.png)
Figure S54. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG ETHER LIPID METABOLISM, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_NOD_LIKE_RECEPTOR_SIGNALING_PATHWAY.volcano_plot.clus3.7.png)
Table S28. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | VEGFC | VEGFC | VEGFC | 274 | 0.34 | 0.067 | YES |
2 | MMP1 | MMP1 | MMP1 | 580 | 0.22 | 0.11 | YES |
3 | TYMP | TYMP | TYMP | 708 | 0.19 | 0.15 | YES |
4 | PGF | PGF | PGF | 771 | 0.18 | 0.19 | YES |
5 | IL8 | IL8 | IL8 | 971 | 0.15 | 0.21 | YES |
6 | CDKN1A | CDKN1A | CDKN1A | 987 | 0.15 | 0.25 | YES |
7 | HRAS | HRAS | HRAS | 1050 | 0.14 | 0.28 | YES |
8 | MAPK3 | MAPK3 | MAPK3 | 1411 | 0.11 | 0.29 | YES |
9 | CDH1 | CDH1 | CDH1 | 1799 | 0.087 | 0.29 | YES |
10 | FGFR3 | FGFR3 | FGFR3 | 1852 | 0.085 | 0.31 | YES |
11 | MYC | MYC | MYC | 2028 | 0.078 | 0.32 | YES |
12 | DAPK3 | DAPK3 | DAPK3 | 2664 | 0.057 | 0.3 | NO |
13 | NRAS | NRAS | NRAS | 3390 | 0.04 | 0.27 | NO |
14 | RB1 | RB1 | RB1 | 3534 | 0.037 | 0.27 | NO |
15 | MAP2K2 | MAP2K2 | MAP2K2 | 4008 | 0.029 | 0.25 | NO |
16 | VEGFA | VEGFA | VEGFA | 4355 | 0.023 | 0.24 | NO |
17 | EGFR | EGFR | EGFR | 4701 | 0.017 | 0.22 | NO |
18 | MAP2K1 | MAP2K1 | MAP2K1 | 5031 | 0.012 | 0.21 | NO |
19 | ERBB2 | ERBB2 | ERBB2 | 5291 | 0.0073 | 0.2 | NO |
20 | THBS1 | THBS1 | THBS1 | 5679 | 0.0014 | 0.18 | NO |
21 | RASSF1 | RASSF1 | RASSF1 | 5735 | 0.00049 | 0.17 | NO |
22 | MAPK1 | MAPK1 | MAPK1 | 6457 | -0.01 | 0.14 | NO |
23 | VEGFB | VEGFB | VEGFB | 6522 | -0.011 | 0.13 | NO |
24 | CCND1 | CCND1 | CCND1 | 7375 | -0.024 | 0.093 | NO |
25 | RAF1 | RAF1 | RAF1 | 7446 | -0.025 | 0.096 | NO |
26 | ARAF | ARAF | ARAF | 7927 | -0.033 | 0.078 | NO |
27 | MMP9 | MMP9 | MMP9 | 8616 | -0.043 | 0.051 | NO |
28 | BRAF | BRAF | BRAF | 8631 | -0.044 | 0.061 | NO |
29 | KRAS | KRAS | KRAS | 8651 | -0.044 | 0.07 | NO |
30 | E2F2 | E2F2 | E2F2 | 9564 | -0.06 | 0.036 | NO |
31 | E2F3 | E2F3 | E2F3 | 9790 | -0.065 | 0.039 | NO |
32 | MDM2 | MDM2 | MDM2 | 9798 | -0.065 | 0.055 | NO |
33 | DAPK1 | DAPK1 | DAPK1 | 10195 | -0.073 | 0.051 | NO |
34 | CDK4 | CDK4 | CDK4 | 10389 | -0.077 | 0.06 | NO |
35 | MMP2 | MMP2 | MMP2 | 10530 | -0.08 | 0.072 | NO |
36 | DAPK2 | DAPK2 | DAPK2 | 11750 | -0.11 | 0.032 | NO |
37 | RPS6KA5 | RPS6KA5 | RPS6KA5 | 12198 | -0.12 | 0.037 | NO |
38 | TP53 | TP53 | TP53 | 13719 | -0.17 | -0.0033 | NO |
39 | CDKN2A | CDKN2A | CDKN2A | 13977 | -0.18 | 0.027 | NO |
40 | EGF | EGF | EGF | 17392 | -0.39 | -0.064 | NO |
41 | FIGF | FIGF | FIGF | 17999 | -0.47 | 0.02 | NO |
Figure S55. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG AMINOACYL TRNA BIOSYNTHESIS.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_AMINOACYL_TRNA_BIOSYNTHESIS.heatmap_plot.clus3.4.png)
Figure S56. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG AMINOACYL TRNA BIOSYNTHESIS, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_BLADDER_CANCER.volcano_plot.clus3.8.png)
Table S29. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | CDKN1A | CDKN1A | CDKN1A | 987 | 0.15 | 0.042 | YES |
2 | HRAS | HRAS | HRAS | 1050 | 0.14 | 0.13 | YES |
3 | CCNE1 | CCNE1 | CCNE1 | 1219 | 0.13 | 0.2 | YES |
4 | MAPK3 | MAPK3 | MAPK3 | 1411 | 0.11 | 0.26 | YES |
5 | CDK6 | CDK6 | CDK6 | 1884 | 0.084 | 0.28 | YES |
6 | RAC1 | RAC1 | RAC1 | 2855 | 0.052 | 0.26 | NO |
7 | RB1 | RB1 | RB1 | 3534 | 0.037 | 0.25 | NO |
8 | TFDP1 | TFDP1 | TFDP1 | 3646 | 0.035 | 0.27 | NO |
9 | RELA | RELA | RELA | 4106 | 0.027 | 0.26 | NO |
10 | RHOA | RHOA | RHOA | 4687 | 0.017 | 0.24 | NO |
11 | NFKBIA | NFKBIA | NFKBIA | 5419 | 0.0052 | 0.2 | NO |
12 | NFKB1 | NFKB1 | NFKB1 | 5544 | 0.0032 | 0.2 | NO |
13 | PAK1 | PAK1 | PAK1 | 5584 | 0.0028 | 0.2 | NO |
14 | IKBKG | IKBKG | IKBKG | 5923 | -0.0025 | 0.18 | NO |
15 | CHUK | CHUK | CHUK | 6098 | -0.0053 | 0.17 | NO |
16 | MAPK1 | MAPK1 | MAPK1 | 6457 | -0.01 | 0.16 | NO |
17 | AKT1 | AKT1 | AKT1 | 7134 | -0.02 | 0.14 | NO |
18 | CCND1 | CCND1 | CCND1 | 7375 | -0.024 | 0.14 | NO |
19 | RAF1 | RAF1 | RAF1 | 7446 | -0.025 | 0.15 | NO |
20 | IKBKB | IKBKB | IKBKB | 9933 | -0.068 | 0.056 | NO |
21 | CDK4 | CDK4 | CDK4 | 10389 | -0.077 | 0.08 | NO |
22 | CDK2 | CDK2 | CDK2 | 11468 | -0.1 | 0.085 | NO |
23 | PIK3CA | PIK3CA | PIK3CA | 12036 | -0.12 | 0.13 | NO |
24 | CDKN1B | CDKN1B | CDKN1B | 12302 | -0.12 | 0.19 | NO |
25 | PIK3R1 | PIK3R1 | PIK3R1 | 14952 | -0.22 | 0.19 | NO |
Figure S57. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG RIBOSOME.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_RIBOSOME.heatmap_plot.clus3.3.png)
Figure S58. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG RIBOSOME, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.BIOCARTA_RACCYCD_PATHWAY.volcano_plot.clus3.9.png)
Table S30. Get Full Table This table shows a Running Enrichment Score (RES) of each gene in this pathway, that is, the enrichment score at this point in the ranked list of genes. All genes are ranked by Signal-to-Noise (S2N), a measure of similarity as default and are used to obtain ES matrix of all genes. In this way, GSEA tool uses expression pattern of not only overlapped genes but also not-overlapped genes to produce ES matrix.
Rank | GENE | SYMBOL | DESC | LIST.LOC | S2N | RES | CORE_ENRICHMENT |
---|---|---|---|---|---|---|---|
1 | STX11 | STX11 | STX11 | 799 | 0.18 | 0.032 | YES |
2 | STX19 | STX19 | STX19 | 1357 | 0.12 | 0.051 | YES |
3 | VAMP8 | VAMP8 | VAMP8 | 1786 | 0.088 | 0.065 | YES |
4 | STX3 | STX3 | STX3 | 1935 | 0.082 | 0.091 | YES |
5 | VAMP3 | VAMP3 | VAMP3 | 2084 | 0.076 | 0.12 | YES |
6 | VTI1B | VTI1B | VTI1B | 2109 | 0.075 | 0.15 | YES |
7 | YKT6 | YKT6 | YKT6 | 2176 | 0.073 | 0.17 | YES |
8 | STX4 | STX4 | STX4 | 2516 | 0.061 | 0.18 | YES |
9 | STX10 | STX10 | STX10 | 2559 | 0.06 | 0.2 | YES |
10 | STX8 | STX8 | STX8 | 2874 | 0.051 | 0.21 | YES |
11 | STX6 | STX6 | STX6 | 3187 | 0.044 | 0.21 | YES |
12 | BNIP1 | BNIP1 | BNIP1 | 3359 | 0.041 | 0.22 | YES |
13 | STX17 | STX17 | STX17 | 3386 | 0.04 | 0.23 | YES |
14 | VTI1A | VTI1A | VTI1A | 3405 | 0.04 | 0.25 | YES |
15 | STX7 | STX7 | STX7 | 3528 | 0.037 | 0.26 | YES |
16 | STX12 | STX12 | STX12 | 3587 | 0.036 | 0.27 | YES |
17 | GOSR2 | GOSR2 | GOSR2 | 4102 | 0.027 | 0.25 | NO |
18 | STX5 | STX5 | STX5 | 4269 | 0.024 | 0.26 | NO |
19 | STX1A | STX1A | STX1A | 5359 | 0.0062 | 0.2 | NO |
20 | BET1L | BET1L | BET1L | 5493 | 0.004 | 0.19 | NO |
21 | STX16 | STX16 | STX16 | 5936 | -0.0028 | 0.17 | NO |
22 | STX18 | STX18 | STX18 | 6018 | -0.0042 | 0.17 | NO |
23 | VAMP7 | VAMP7 | VAMP7 | 6857 | -0.016 | 0.13 | NO |
24 | VAMP2 | VAMP2 | VAMP2 | 7309 | -0.023 | 0.11 | NO |
25 | GOSR1 | GOSR1 | GOSR1 | 7331 | -0.024 | 0.12 | NO |
26 | VAMP4 | VAMP4 | VAMP4 | 7345 | -0.024 | 0.13 | NO |
27 | USE1 | USE1 | USE1 | 7988 | -0.034 | 0.11 | NO |
28 | SEC22B | SEC22B | SEC22B | 8019 | -0.034 | 0.12 | NO |
29 | SNAP23 | SNAP23 | SNAP23 | 8277 | -0.038 | 0.13 | NO |
30 | BET1 | BET1 | BET1 | 8514 | -0.042 | 0.13 | NO |
31 | STX2 | STX2 | STX2 | 8764 | -0.046 | 0.14 | NO |
32 | STX1B | STX1B | STX1B | 10203 | -0.073 | 0.09 | NO |
33 | SNAP47 | SNAP47 | SNAP47 | 10865 | -0.087 | 0.091 | NO |
34 | VAMP5 | VAMP5 | VAMP5 | 10938 | -0.089 | 0.12 | NO |
35 | SNAP25 | SNAP25 | SNAP25 | 13765 | -0.17 | 0.044 | NO |
36 | TSNARE1 | TSNARE1 | TSNARE1 | 13909 | -0.18 | 0.11 | NO |
37 | VAMP1 | VAMP1 | VAMP1 | 16436 | -0.31 | 0.1 | NO |
Figure S59. Get High-res Image This plot shows mRNAseq_cNMF expression data heatmap (on the left) a RunningEnrichmentScore(RES) plot (on the top right) and a Signal2Noise(S2N) plot (on the bottom right) of genes in the pathway: KEGG RNA POLYMERASE.
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_RNA_POLYMERASE.heatmap_plot.clus3.13.png)
Figure S60. Get High-res Image For the top 5 core enriched genes in the pathway: KEGG RNA POLYMERASE, this volcano plot shows how much they are up/down-regulated and significant. The significance was calculated by empirical bayesian fit
![](HNSC-TP-mRNAseq_cNMF-clus3___Class2-CanonicalPathway.KEGG_SNARE_INTERACTIONS_IN_VESICULAR_TRANSPORT.volcano_plot.clus3.10.png)
For the top enriched genes, if you want to check whether they are
-
up-regulated, please check the list of up-regulated genes
-
down-regulated, please check the list of down-regulated genes
For the top enriched genes, if you want to check whether they are
-
highly expressed genes, please check the list of high (top 30%) expressed genes
-
low expressed genes, please check the list of low (bottom 30%) expressed genes
An expression pattern of top(30%)/middle(30%)/low(30%) in this subtype against other subtypes is available in a heatmap
For the top enriched genes, if you want to check whether they are
-
significantly differently expressed genes by eBayes lm fit, please check the list of significant genes
-
Gene set database = c2.cp.v3.0-2.symbols.gmt
-
Expression data file = HNSC-TP.uncv2.mRNAseq_RSEM_normalized_log2.txt
-
Phenotype data file = HNSC-TP.mergedcluster.txt
For the Gene Set Enrichment Analysis (GSEA), Broad GSEA-P-R.1.0 version is used with class2: canonical pathways geneses from MSigDB. Further details about statistics are available inThe Broad GSEA website.
In addition to the links below, the full results of the analysis summarized in this report can also be downloaded programmatically using firehose_get, or interactively from either the Broad GDAC website or TCGA Data Coordination Center Portal.