Mitch Report
Antony Kaspi & Mark Ziemann
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.
## x
## 1MMP -0.06107556
## 2-Cys 0.73109293
## 2A6 2.58510884
## 2MMP -2.37940783
## 3AT1 -0.29989906
## 3BETAHSD/D1 2.11566572| Profile metrics | |
|---|---|
| num_genesets | 4 |
| num_genes_in_profile | 21486 |
| duplicated_genes_present | 0 |
| num_profile_genes_in_sets | 156 |
| num_profile_genes_not_in_sets | 21330 |
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 | |
|---|---|
| num_genesets | 4 |
| num_genesets_excluded | 0 |
| num_genesets_included | 4 |
Differential pathway expression
Interactive enrichment scatterplot
Significance is calculated by -log10(p-value). All points shown are FDR<0.05.
Significance is calculated by -log10(p-value). Top N sets shown irrespective of FDR.
Results table
Top N= 4 gene sets| set | setSize | pANOVA | s.dist | p.adjustANOVA |
|---|---|---|---|---|
| MISC PEROXIDASES | 63 | 0.0014 | 0.2330 | 0.00559 |
| REDOX GLUTAREDOXINS | 29 | 0.0417 | 0.2190 | 0.05550 |
| MISC GLUTATHIONE S TRANSFERASES | 53 | 0.0131 | -0.1970 | 0.02610 |
| REDOX DISMUTASES AND CATALASES | 11 | 0.9100 | 0.0197 | 0.91000 |
Results (complete table)
Click HERE to show results for all gene sets
| set | setSize | pANOVA | s.dist | p.adjustANOVA |
|---|---|---|---|---|
| MISC PEROXIDASES | 63 | 0.0014 | 0.2330 | 0.00559 |
| REDOX GLUTAREDOXINS | 29 | 0.0417 | 0.2190 | 0.05550 |
| MISC GLUTATHIONE S TRANSFERASES | 53 | 0.0131 | -0.1970 | 0.02610 |
| REDOX DISMUTASES AND CATALASES | 11 | 0.9100 | 0.0197 | 0.91000 |
Detailed Gene set reports
MISC_PEROXIDASES
| 2 | |
|---|---|
| set | MISC_PEROXIDASES |
| setSize | 63 |
| pANOVA | 0.0014 |
| s.dist | 0.233 |
| p.adjustANOVA | 0.00559 |
Top enriched genes
| GeneID | Gene Rank |
|---|---|
| PER27 | 10114 |
| PER23 | 10099 |
| AT5G14130 | 10060 |
| PER2 | 10033 |
| PER50 | 10004 |
| PER71 | 9859 |
| AT4G26010 | 9818 |
| PER69 | 9755 |
| PER42 | 9597 |
| PER73 | 9580 |
| AT5G17820 | 9571 |
| PER10 | 9406 |
| PER39 | 9332 |
| PER32 | 9149 |
| PER37 | 9057 |
| PER30 | 8976 |
| PER51 | 8633 |
| PER20 | 8463 |
| PER63 | 8374 |
| PER62 | 8314 |
Click HERE to show all gene set members
| GeneID | Gene Rank |
|---|---|
| PER27 | 10114 |
| PER23 | 10099 |
| AT5G14130 | 10060 |
| PER2 | 10033 |
| PER50 | 10004 |
| PER71 | 9859 |
| AT4G26010 | 9818 |
| PER69 | 9755 |
| PER42 | 9597 |
| PER73 | 9580 |
| AT5G17820 | 9571 |
| PER10 | 9406 |
| PER39 | 9332 |
| PER32 | 9149 |
| PER37 | 9057 |
| PER30 | 8976 |
| PER51 | 8633 |
| PER20 | 8463 |
| PER63 | 8374 |
| PER62 | 8314 |
| PER16 | 7908 |
| PER7 | 7409 |
| PER61 | 7108 |
| PER12 | 7083 |
| PER24 | 7075 |
| PER64 | 6644 |
| PER60 | 6494 |
| PER19 | 5618 |
| PER8 | 5251 |
| PER35 | 4660 |
| PER56 | 4615 |
| PER59 | 4440 |
| PER9 | 4388 |
| PER17 | 4171 |
| PER43 | 3713 |
| PER11 | 3557 |
| PER72 | 3272 |
| PER49 | 2852 |
| PER31 | 2131 |
| PER45 | -226 |
| PER22 | -260 |
| PER26 | -1044 |
| PER70 | -1500 |
| PER53 | -2159 |
| PER33 | -4814 |
| PER48 | -4861 |
| APX6 | -5038 |
| PER52 | -5227 |
| PER40 | -5432 |
| PER66 | -6462 |
| PER46 | -7571 |
| PER14 | -8828 |
| PER21 | -9105 |
| PER67 | -9395 |
| PER15 | -9759 |
| PER65 | -10100 |
| PER28 | -10235 |
| PER58 | -10240 |
| PER38 | -10286 |
| PER25 | -10818 |
| AT4G33420 | -10899 |
| PER34 | -11198 |
| PER54 | -11265 |
REDOX_GLUTAREDOXINS
| 4 | |
|---|---|
| set | REDOX_GLUTAREDOXINS |
| setSize | 29 |
| pANOVA | 0.0417 |
| s.dist | 0.219 |
| p.adjustANOVA | 0.0555 |
Top enriched genes
| GeneID | Gene Rank |
|---|---|
| GRXS9 | 9750 |
| GRXC14 | 9734 |
| AT5G13810 | 9653 |
| GRXS11 | 9173 |
| GRXS6 | 8474 |
| GRXC10 | 8429 |
| AT5G58530 | 7440 |
| AT1G03850 | 7275 |
| GRXC8 | 7165 |
| AT3G28850 | 6703 |
| GRXS1 | 6521 |
| AT5G39865 | 6102 |
| AT2G41330 | 4390 |
| AT3G57070 | 3381 |
| GRXC6 | 1475 |
| AT5G01420 | 1188 |
| AT2G20270 | 588 |
| GRXC9 | -360 |
| AT4G10630 | -1239 |
| GRXC5 | -1629 |
Click HERE to show all gene set members
| GeneID | Gene Rank |
|---|---|
| GRXS9 | 9750 |
| GRXC14 | 9734 |
| AT5G13810 | 9653 |
| GRXS11 | 9173 |
| GRXS6 | 8474 |
| GRXC10 | 8429 |
| AT5G58530 | 7440 |
| AT1G03850 | 7275 |
| GRXC8 | 7165 |
| AT3G28850 | 6703 |
| GRXS1 | 6521 |
| AT5G39865 | 6102 |
| AT2G41330 | 4390 |
| AT3G57070 | 3381 |
| GRXC6 | 1475 |
| AT5G01420 | 1188 |
| AT2G20270 | 588 |
| GRXC9 | -360 |
| AT4G10630 | -1239 |
| GRXC5 | -1629 |
| AT1G77370 | -3026 |
| GRXC4 | -3078 |
| AT5G03870 | -3440 |
| AT1G32760 | -3954 |
| GRXS15 | -4477 |
| AT1G64500 | -7091 |
| GRXC1 | -7926 |
| AT5G06470 | -9605 |
| AT5G40370 | -10751 |
MISC_GLUTATHIONE_S_TRANSFERASES
| 1 | |
|---|---|
| set | MISC_GLUTATHIONE_S_TRANSFERASES |
| setSize | 53 |
| pANOVA | 0.0131 |
| s.dist | -0.197 |
| p.adjustANOVA | 0.0261 |
Top enriched genes
| GeneID | Gene Rank |
|---|---|
| GSTU10 | -11215 |
| GSTF7 | -11164 |
| GSTU4 | -10996 |
| GSTF3 | -10717 |
| GSTU18 | -10562 |
| ATGSTZ1 | -10552 |
| GSTF2 | -10487 |
| GSTU9 | -10469 |
| GSTU3 | -10313 |
| GSTL3 | -10298 |
| GSTU5 | -9565 |
| GSTU6 | -9341 |
| GSTL1 | -9261 |
| GSTU1 | -9151 |
| GSTF6 | -8805 |
| GSTU16 | -8703 |
| GSTU13 | -8171 |
| GSTT2 | -8151 |
| AT1G65820 | -7787 |
| GSTZ2 | -7481 |
Click HERE to show all gene set members
| GeneID | Gene Rank |
|---|---|
| GSTU10 | -11215 |
| GSTF7 | -11164 |
| GSTU4 | -10996 |
| GSTF3 | -10717 |
| GSTU18 | -10562 |
| ATGSTZ1 | -10552 |
| GSTF2 | -10487 |
| GSTU9 | -10469 |
| GSTU3 | -10313 |
| GSTL3 | -10298 |
| GSTU5 | -9565 |
| GSTU6 | -9341 |
| GSTL1 | -9261 |
| GSTU1 | -9151 |
| GSTF6 | -8805 |
| GSTU16 | -8703 |
| GSTU13 | -8171 |
| GSTT2 | -8151 |
| AT1G65820 | -7787 |
| GSTZ2 | -7481 |
| GSTU27 | -7263 |
| GSTU24 | -7022 |
| GSTU2 | -6915 |
| TCHQD | -6619 |
| ATGSTU21 | -6563 |
| GSTU19 | -6400 |
| GSTU26 | -5474 |
| GSTF12 | -4887 |
| GSTU7 | -4844 |
| GSTF13 | -4141 |
| GSTU23 | -1452 |
| GSTT3 | -457 |
| ATGSTF4 | -192 |
| GSTT1 | 101 |
| AT5G44000 | 760 |
| GSTU8 | 880 |
| GSTU25 | 2840 |
| GSTU11 | 3161 |
| GSTF10 | 3164 |
| GSTU17 | 3188 |
| GSTL2 | 4133 |
| GSTF9 | 5482 |
| GSTU14 | 5849 |
| GSTU12 | 6625 |
| GSTU28 | 6719 |
| GSTF14 | 7467 |
| AT4G19880 | 8035 |
| GSTU20 | 8064 |
| GSTU22 | 8398 |
| GSTF11 | 8608 |
| GSTF8 | 8724 |
| PER4 | 9967 |
| PER5 | 10025 |
REDOX_DISMUTASES_AND_CATALASES
| 3 | |
|---|---|
| set | REDOX_DISMUTASES_AND_CATALASES |
| setSize | 11 |
| pANOVA | 0.91 |
| s.dist | 0.0197 |
| p.adjustANOVA | 0.91 |
Top enriched genes
| GeneID | Gene Rank |
|---|---|
| CCS | 9855 |
| CAT3 | 9324 |
| CSD1 | 7556 |
| CSD2 | 7342 |
| CAT2 | 1505 |
| CSD3 | -3601 |
| FSD3 | -3701 |
| MSD1 | -6460 |
| CAT1 | -6928 |
| FSD2 | -8761 |
| FSD1 | -10303 |
Click HERE to show all gene set members
| GeneID | Gene Rank |
|---|---|
| CCS | 9855 |
| CAT3 | 9324 |
| CSD1 | 7556 |
| CSD2 | 7342 |
| CAT2 | 1505 |
| CSD3 | -3601 |
| FSD3 | -3701 |
| MSD1 | -6460 |
| CAT1 | -6928 |
| FSD2 | -8761 |
| FSD1 | -10303 |
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.3END of report