date generated: 2021-03-25

Background

Mitch performs unidimensional and multidimensional gene set enrichment analysis. The concept behind this dates to work by Cox and Mann (https://doi.org/10.1186/1471-2105-13-S16-S12). This implementation is suited to R based workflows of multi-omics datasets. This software was developed by Antony Kaspi and Mark Ziemann. Learn more about Mitch at the website: https://github.com/markziemann/Mitch

Input profiles

Here is the first few lines of the input profile.

##                   del3        del5        der3        der5
## 1MMP        -2.1945558  -0.2573510 -0.06107556 -0.06943941
## 2-Cys        1.6562904   2.8083698  0.73109293  1.11074220
## 2A6          1.7809241   2.9365255  2.58510884  0.87763168
## 2MMP        -4.3494649 -10.4988977 -2.37940783 -6.21455109
## 3AT1        -0.6118777   1.3787620 -0.29989906  0.03357734
## 3BETAHSD/D1  0.9976927   0.3206492  2.11566572  2.94991600
Here are some metrics about the input data profile:
Profile metrics
Profile metrics
num_genes_in_profile 19846
duplicated_genes_present 0
num_profile_genes_in_sets 143
num_profile_genes_not_in_sets 19703

Here is a plot of the input profiles. Note the dynamic ranges.

Here is the contour plot of the profile including all detected genes.

Input genesets

Here are some metrics about the gene sets used: GMT file of genesets:
Gene sets metrics
Gene sets metrics
num_genesets 4
num_genesets_excluded 0
num_genesets_included 4
Genes by sector
del3 del5 der3 der5 Count
1 -1 -1 -1 -1 3244
2 1 -1 -1 -1 870
3 -1 1 -1 -1 1413
4 1 1 -1 -1 1058
5 -1 -1 1 -1 1028
6 1 -1 1 -1 515
7 1 0 1 -1 1
8 -1 1 1 -1 728
9 0 1 1 -1 1
10 1 1 1 -1 813
11 -1 -1 1 0 1
12 -1 -1 -1 1 1036
13 1 -1 -1 1 1006
14 -1 1 -1 1 674
15 1 1 -1 1 1347
16 -1 -1 0 1 1
17 -1 -1 1 1 993
18 1 -1 1 1 1283
19 -1 1 1 1 863
20 1 1 1 1 2971
Number of significant gene sets (FDR<0.05)= 4

Gene sets by sector


Gene sets by sector
s.del3 s.del5 s.der3 s.der5 Count
1 -1 1 -1 -1 1
2 -1 1 1 -1 1
3 1 1 1 -1 1
4 -1 -1 1 1 1

Interactive enrichment scatterplot

All sets with FDR<0.05. Try hovering over the points.

Top N sets irrespective of FDR. Try hovering over the points.

A heatmap of S values for top results


A plot of effect size versus significance


Significance is the -log2(p.adjustMANOVA) and effect size is the s.dist which is the hypotenuse of the s scores.

Results table

Top N= 4 gene sets
set setSize pMANOVA p.adjustMANOVA s.dist s.del3 s.del5 s.der3 s.der5 p.del3 p.del5 p.der3 p.der5
REDOX DISMUTASES AND CATALASES 11 6.73e-03 6.73e-03 0.561 0.2670 0.3830 0.0284 -0.3080 0.1250 0.0277 0.87000 7.66e-02
MISC GLUTATHIONE S TRANSFERASES 49 5.57e-06 2.23e-05 0.442 -0.0191 0.0802 -0.1910 -0.3900 0.8170 0.3310 0.02060 2.26e-06
REDOX GLUTAREDOXINS 28 4.81e-03 6.41e-03 0.326 -0.1490 0.1000 0.2130 -0.1690 0.1710 0.3590 0.05160 1.22e-01
MISC PEROXIDASES 55 9.57e-05 1.91e-04 0.321 -0.1590 -0.0030 0.2790 0.0165 0.0418 0.9690 0.00035 8.32e-01


Results (complete table)


Click HERE to show results for all gene sets

Complete results
set setSize pMANOVA p.adjustMANOVA s.dist s.del3 s.del5 s.der3 s.der5 p.del3 p.del5 p.der3 p.der5
REDOX DISMUTASES AND CATALASES 11 6.73e-03 6.73e-03 0.561 0.2670 0.3830 0.0284 -0.3080 0.1250 0.0277 0.87000 7.66e-02
MISC GLUTATHIONE S TRANSFERASES 49 5.57e-06 2.23e-05 0.442 -0.0191 0.0802 -0.1910 -0.3900 0.8170 0.3310 0.02060 2.26e-06
REDOX GLUTAREDOXINS 28 4.81e-03 6.41e-03 0.326 -0.1490 0.1000 0.2130 -0.1690 0.1710 0.3590 0.05160 1.22e-01
MISC PEROXIDASES 55 9.57e-05 1.91e-04 0.321 -0.1590 -0.0030 0.2790 0.0165 0.0418 0.9690 0.00035 8.32e-01



Detailed Gene set reports


REDOX_DISMUTASES_AND_CATALASES

REDOX_DISMUTASES_AND_CATALASES
metric value
setSize 11
pMANOVA 0.00673
p.adjustMANOVA 0.00673
s.dist 0.561
s.del3 0.267
s.del5 0.383
s.der3 0.0284
s.der5 -0.308
p.del3 0.125
p.del5 0.0277
p.der3 0.87
p.der5 0.0766




Top 20 genes
Gene del5 der5
MSD1 7064 -7092
CAT1 6670 -7183
FSD1 3273 -8793
CSD1 3378 -7850
FSD2 9196 -2797
FSD3 2576 -5061

Click HERE to show all gene set members

All member genes
del3 del5 der3 der5
CAT1 8934 6670 -6527 -7183
CAT2 -8778 -5571 1409 -9321
CAT3 9270 9789 8505 5140
CCS 8616 4946 8953 8679
CSD1 4813 3378 6942 -7850
CSD2 6885 7242 6754 2279
CSD3 2848 -7354 -3392 1128
FSD1 -5776 3273 -9691 -8793
FSD2 9383 9196 -8258 -2797
FSD3 497 2576 -3484 -5061
MSD1 -8179 7064 -6089 -7092





MISC_GLUTATHIONE_S_TRANSFERASES

MISC_GLUTATHIONE_S_TRANSFERASES
metric value
setSize 49
pMANOVA 5.57e-06
p.adjustMANOVA 2.23e-05
s.dist 0.442
s.del3 -0.0191
s.del5 0.0802
s.der3 -0.191
s.der5 -0.39
p.del3 0.817
p.del5 0.331
p.der3 0.0206
p.der5 2.26e-06




Top 20 genes
Gene der5 der3
GSTU10 -9625 -10538
GSTF7 -9654 -10490
GSTU4 -9278 -10335
GSTF2 -9461 -9865
GSTU9 -9462 -9849
GSTF3 -9146 -10077
ATGSTZ1 -8920 -9921
GSTU18 -8596 -9931
GSTL1 -9565 -8726
GSTU3 -8316 -9700
GSTF6 -9580 -8299
GSTL3 -8187 -9686
GSTU5 -8432 -9003
GSTU1 -8649 -8622
AT1G65820 -9375 -7350
GSTU13 -8243 -7706
GSTU27 -9003 -6851
GSTU19 -8307 -6033
GSTU6 -5638 -8802
GSTU24 -6761 -6617

Click HERE to show all gene set members

All member genes
del3 del5 der3 der5
AT1G65820 -2238 -4736 -7350 -9375
AT4G19880 9060 9087 7370 7403
AT5G44000 -2848 4702 712 -6781
ATGSTF4 -9587 -3752 -177 -8822
ATGSTZ1 8312 2168 -9921 -8920
GSTF10 -5867 1802 2931 -6903
GSTF11 9591 9841 7877 8156
GSTF12 96 5834 -4594 4265
GSTF14 -7603 -5748 6866 -508
GSTF2 4208 -1602 -9865 -9461
GSTF3 4163 -2118 -10077 -9146
GSTF6 948 324 -8299 -9580
GSTF7 -6108 -3332 -10490 -9654
GSTF8 963 285 7979 9623
GSTF9 8250 9616 5049 -359
GSTL1 -8096 -9764 -8726 -9565
GSTL2 8003 6667 3813 1538
GSTL3 -9565 -1859 -9686 -8187
GSTT1 6718 6644 95 -5105
GSTT2 2809 -7491 -7687 2357
GSTT3 2494 -7272 -428 5928
GSTU1 8158 9425 -8622 -8649
GSTU10 3641 -7843 -10538 -9625
GSTU11 7116 3573 2928 -6383
GSTU12 -5593 -60 6105 -4913
GSTU13 6517 7143 -7706 -8243
GSTU14 3927 5944 5389 5334
GSTU17 -5467 2248 2951 1114
GSTU18 384 6758 -9931 -8596
GSTU19 -7894 7384 -6033 -8307
GSTU2 8882 44 -6514 -5210
GSTU20 8926 9847 7398 9762
GSTU22 -6639 -2331 7697 4887
GSTU23 -7431 1522 -1359 -5569
GSTU24 3807 8750 -6617 -6761
GSTU25 -9949 -8720 2641 6076
GSTU26 -8809 -2322 -5153 -7774
GSTU27 -5996 -2337 -6851 -9003
GSTU28 -1809 5001 6189 7003
GSTU3 8115 -1190 -9700 -8316
GSTU4 1925 3608 -10335 -9278
GSTU5 -5318 2818 -9003 -8432
GSTU6 -1954 -2286 -8802 -5638
GSTU7 1618 8492 -4554 -5605
GSTU8 -9326 -9239 820 -3741
GSTU9 -9855 -6311 -9849 -9462
GSTZ2 7439 -4008 -7058 1243
PER4 -8140 -734 9050 -4250
TCHQD -2105 -8237 -6239 -5590





REDOX_GLUTAREDOXINS

REDOX_GLUTAREDOXINS
metric value
setSize 28
pMANOVA 0.00481
p.adjustMANOVA 0.00641
s.dist 0.326
s.del3 -0.149
s.del5 0.1
s.der3 0.213
s.der5 -0.169
p.del3 0.171
p.del5 0.359
p.der3 0.0516
p.der5 0.122




Top 20 genes
Gene der3 der5
GRXS9 8866 -6192
GRXS11 8374 -5362
AT5G13810 8782 -4631
AT3G28850 6174 -6552
AT2G20270 552 -8112
GRXS1 6007 -306

Click HERE to show all gene set members

All member genes
del3 del5 der3 der5
AT1G03850 -6017 -1740 6694 218
AT1G32760 524 3790 -3723 -2773
AT1G64500 -7958 3058 -6685 -3126
AT1G77370 258 -1847 -2850 -6576
AT2G20270 -7072 4280 552 -8112
AT2G41330 3480 4370 4053 3074
AT3G28850 1657 -2924 6174 -6552
AT3G57070 4967 4978 3125 344
AT4G10630 -2298 6279 -1160 221
AT5G01420 -571 -3586 1108 747
AT5G03870 -3381 3718 -3241 7308
AT5G06470 -9159 -4960 -9040 -8038
AT5G13810 642 6100 8782 -4631
AT5G39865 6147 5866 5620 6000
AT5G40370 -7759 -1777 -10108 -8735
AT5G58530 -7329 -3705 6841 6336
GRXC1 -9317 -8469 -7478 -8739
GRXC10 -2434 6586 7723 7883
GRXC14 9614 9658 8851 9433
GRXC4 -3570 429 -2900 -6914
GRXC5 923 4007 -1522 -3630
GRXC6 -4846 3862 1380 841
GRXC9 -7893 -3858 -334 -7088
GRXS1 -5735 -8068 6007 -306
GRXS11 8915 8770 8374 -5362
GRXS15 2384 -7900 -4207 -3643
GRXS6 -1416 -2669 7761 8208
GRXS9 4168 2017 8866 -6192





MISC_PEROXIDASES

MISC_PEROXIDASES
metric value
setSize 55
pMANOVA 9.57e-05
p.adjustMANOVA 0.000191
s.dist 0.321
s.del3 -0.159
s.del5 -0.003
s.der3 0.279
s.der5 0.0165
p.del3 0.0418
p.del5 0.969
p.der3 0.00035
p.der5 0.832




Top 20 genes
Gene der3 del3
AT5G14130 9126 -9550
PER20 7751 -9096
PER39 8513 -7355
AT4G26010 8923 -5771
PER62 7621 -5806
PER7 6814 -6216
PER35 4298 -9207
AT5G17820 8713 -4379
PER59 4098 -8296
PER10 8575 -3812
PER11 3290 -9230
PER9 4051 -6251
PER64 6121 -4092
PER71 8957 -2721
PER17 3849 -5498
PER72 3024 -6523
PER56 4255 -3266
PER69 8870 -1492
PER73 8721 -1051
PER8 4837 -1845

Click HERE to show all gene set members

All member genes
del3 del5 der3 der5
APX6 -3040 -6228 -4738 839
AT4G26010 -5771 2535 8923 7616
AT4G33420 -49 -3288 -10241 -9590
AT5G14130 -9550 -7444 9126 7914
AT5G17820 -4379 4737 8713 5345
PER10 -3812 3347 8575 3015
PER11 -9230 -9931 3290 -7904
PER12 4541 1599 6519 5161
PER15 -5003 -5936 -9187 -9372
PER16 3187 6403 7257 9667
PER17 -5498 -6361 3849 -4504
PER19 4718 4008 5175 9018
PER2 1337 4960 9106 9996
PER20 -9096 -8725 7751 -1195
PER21 -6049 -3545 -8581 1611
PER22 -2283 -7557 -240 -9478
PER23 9094 7485 9155 9686
PER25 -9323 -8948 -10167 -8725
PER27 4301 6454 9168 10121
PER28 -9013 -8107 -9624 -9132
PER30 588 8156 8201 -2498
PER31 490 -1375 1990 1921
PER32 2743 3213 8354 3596
PER34 4683 -3724 -10522 -9665
PER35 -9207 -1072 4298 3425
PER37 424 3260 8273 4691
PER38 -3590 -7155 -9674 -1673
PER39 -7355 6197 8513 9165
PER40 8249 4518 -5113 -2204
PER42 8370 9802 8737 7857
PER43 5150 -3243 3434 3153
PER45 -6025 -452 -211 1989
PER49 -1376 -5201 2652 -9148
PER50 6296 5594 9082 9951
PER51 8467 3515 7900 9043
PER52 -6456 -5826 -4915 -9512
PER53 3679 -1569 -2028 -8548
PER54 -6490 -9122 -10583 -9603
PER56 -3266 3005 4255 981
PER58 8295 -1795 -9629 -7977
PER59 -8296 794 4098 -966
PER62 -5806 -20 7621 -1735
PER63 5743 6413 7676 7878
PER64 -4092 3332 6121 6102
PER65 -7782 -9637 -9498 -7892
PER66 1138 4139 -6091 -8431
PER67 198 1078 -8850 66
PER69 -1492 5890 8870 6640
PER7 -6216 3613 6814 -1010
PER70 -3340 -2333 -1406 -2818
PER71 -2721 4833 8957 8767
PER72 -6523 -991 3024 -1457
PER73 -1051 2609 8721 8907
PER8 -1845 1767 4837 -6790
PER9 -6251 1701 4051 514





Here is the session info with all the versions of packages used.

sessionInfo()
## R version 4.0.3 (2020-10-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.2 LTS
## 
## Matrix products: default
## BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.9.0
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.9.0
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
##  [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
##  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] parallel  stats4    stats     graphics  grDevices utils     datasets 
## [8] methods   base     
## 
## other attached packages:
##  [1] pkgload_1.2.0               GGally_2.1.1               
##  [3] ggplot2_3.3.3               beeswarm_0.3.1             
##  [5] gtools_3.8.2                tibble_3.1.0               
##  [7] dplyr_1.0.5                 echarts4r_0.4.0            
##  [9] vioplot_0.3.5               zoo_1.8-9                  
## [11] sm_2.2-5.6                  UpSetR_1.4.0               
## [13] eulerr_6.1.0                kableExtra_1.3.4           
## [15] mitch_1.2.2                 DESeq2_1.30.0              
## [17] SummarizedExperiment_1.20.0 Biobase_2.50.0             
## [19] MatrixGenerics_1.2.0        matrixStats_0.58.0         
## [21] GenomicRanges_1.42.0        GenomeInfoDb_1.26.0        
## [23] IRanges_2.24.0              S4Vectors_0.28.0           
## [25] BiocGenerics_0.36.0         gplots_3.1.1               
## [27] reshape2_1.4.4             
## 
## loaded via a namespace (and not attached):
##  [1] colorspace_2.0-0       ellipsis_0.3.1         rprojroot_2.0.2       
##  [4] XVector_0.30.0         rstudioapi_0.13        farver_2.1.0          
##  [7] bit64_4.0.5            AnnotationDbi_1.52.0   fansi_0.4.2           
## [10] xml2_1.3.2             splines_4.0.3          cachem_1.0.4          
## [13] geneplotter_1.68.0     knitr_1.31             polyclip_1.10-0       
## [16] jsonlite_1.7.2         annotate_1.68.0        shiny_1.6.0           
## [19] compiler_4.0.3         httr_1.4.2             assertthat_0.2.1      
## [22] Matrix_1.3-2           fastmap_1.1.0          later_1.1.0.1         
## [25] prettyunits_1.1.1      htmltools_0.5.1.1      tools_4.0.3           
## [28] gtable_0.3.0           glue_1.4.2             GenomeInfoDbData_1.2.4
## [31] Rcpp_1.0.6             jquerylib_0.1.3        vctrs_0.3.6           
## [34] svglite_2.0.0          polylabelr_0.2.0       xfun_0.22             
## [37] stringr_1.4.0          testthat_3.0.2         rvest_1.0.0           
## [40] mime_0.10              lifecycle_1.0.0        XML_3.99-0.6          
## [43] zlibbioc_1.36.0        MASS_7.3-53.1          scales_1.1.1          
## [46] hms_1.0.0              promises_1.2.0.1       RColorBrewer_1.1-2    
## [49] yaml_2.2.1             memoise_2.0.0          gridExtra_2.3         
## [52] sass_0.3.1             reshape_0.8.8          stringi_1.5.3         
## [55] RSQLite_2.2.4          highr_0.8              genefilter_1.72.0     
## [58] desc_1.3.0             caTools_1.18.1         BiocParallel_1.24.1   
## [61] rlang_0.4.10           pkgconfig_2.0.3        systemfonts_1.0.1     
## [64] bitops_1.0-6           evaluate_0.14          lattice_0.20-41       
## [67] purrr_0.3.4            labeling_0.4.2         htmlwidgets_1.5.3     
## [70] bit_4.0.4              tidyselect_1.1.0       plyr_1.8.6            
## [73] magrittr_2.0.1         R6_2.5.0               generics_0.1.0        
## [76] DelayedArray_0.16.0    DBI_1.1.1              withr_2.4.1           
## [79] pillar_1.5.1           survival_3.2-10        RCurl_1.98-1.3        
## [82] crayon_1.4.1           KernSmooth_2.23-18     utf8_1.2.1            
## [85] rmarkdown_2.7          progress_1.2.2         locfit_1.5-9.4        
## [88] grid_4.0.3             blob_1.2.1             digest_0.6.27         
## [91] webshot_0.5.2          xtable_1.8-4           httpuv_1.5.5          
## [94] munsell_0.5.0          viridisLite_0.3.0      bslib_0.2.4           
## [97] tcltk_4.0.3

END of report