A Single Shiny Solution from Sequencer to Pathway Analysis
Bioinformatics Unit, Core Facility,
Faculty of Medicine and Health Sciences
Linköping University
jyotirmoy.das@liu.se
Disclaimer - part of this presentation is under development and generated from a research project.
In the top figure, the unmethylated CpG target site matches with the U probe, enabling single-base extension and detection. It has a single-base mismatch to the M probe, which inhibits extension. If the CpG locus of interest is methylated (bottom figure), the reverse occurs.
Illumina arrays measure both methylated (C) and unmethylated (T) signal to get proportion of methylation at a CpG.
\[M = log_2 \frac{Meth}{Unmeth}\]
Better statistical properties, recommended for statistical testing
Remove bad samples and poor performing probes (CpGs)
Transform data to remove unwanted variation
Estimate means and variancs and borrow information acorss probes
Biological interpretation of results
e.g., gene expression
minfi, ChAMP, methylumi, limma, RnBeads, …
minfi, ChAMP, methylumi, limma, RnBeads, waterMelon, …
limma, bumhunter, DMRcate, …
Gviz, IGV, …
GenomicRanges
<a href="https://sourceforge.net/projects/methylr/" target="_blank">methylr</a>
What happens after you upload your data –
Please remember:
methylR Component tools
methylR 2.0 Component tools
[1] "R version 4.3.2 (2023-10-31)"
[2] "Platform: x86_64-pc-linux-gnu (64-bit)"
[3] "Running under: Ubuntu 22.04.4 LTS"
[4] ""
[5] "Matrix products: default"
[6] "BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3 "
[7] "LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.20.so; LAPACK version 3.10.0"
[8] ""
[9] "locale:"
[10] " [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C "
[11] " [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 "
[12] " [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 "
[13] " [7] LC_PAPER=en_US.UTF-8 LC_NAME=C "
[14] " [9] LC_ADDRESS=C LC_TELEPHONE=C "
[15] "[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C "
[16] ""
[17] "time zone: Etc/UTC"
[18] "tzcode source: system (glibc)"
[19] ""
[20] "attached base packages:"
[21] " [1] grid parallel stats4 stats graphics grDevices utils "
[22] " [8] datasets methods base "
[23] ""
[24] "other attached packages:"
[25] " [1] doRNG_1.8.6 "
[26] " [2] rngtools_1.5.2 "
[27] " [3] ChAMP_2.32.0 "
[28] " [4] RPMM_1.25 "
[29] " [5] cluster_2.1.4 "
[30] " [6] Illumina450ProbeVariants.db_1.38.0 "
[31] " [7] DMRcate_2.99.5 "
[32] " [8] dendextend_1.17.1 "
[33] " [9] gplots_3.1.3.1 "
[34] " [10] RColorBrewer_1.1-3 "
[35] " [11] tippy_0.1.0 "
[36] " [12] conumee_1.36.0 "
[37] " [13] IlluminaHumanMethylationEPICmanifest_1.0.0 "
[38] " [14] IlluminaHumanMethylationEPICanno.ilm10b2.hg19_0.6.0"
[39] " [15] IlluminaHumanMethylation450kmanifest_0.4.0 "
[40] " [16] IlluminaHumanMethylation450kanno.ilmn12.hg19_0.6.1 "
[41] " [17] dnaMethyAge_0.2.0 "
[42] " [18] pROC_1.18.5 "
[43] " [19] randomForest_4.7-1.1 "
[44] " [20] caret_6.0-94 "
[45] " [21] lattice_0.21-9 "
[46] " [22] lubridate_1.9.3 "
[47] " [23] forcats_1.0.0 "
[48] " [24] stringr_1.5.1 "
[49] " [25] purrr_1.0.2 "
[50] " [26] tibble_3.2.1 "
[51] " [27] tidyverse_2.0.0 "
[52] " [28] dlookr_0.6.3 "
[53] " [29] ReactomePA_1.46.0 "
[54] " [30] org.Hs.eg.db_3.18.0 "
[55] " [31] GOSemSim_2.28.1 "
[56] " [32] GO.db_3.18.0 "
[57] " [33] DOSE_3.28.2 "
[58] " [34] DO.db_2.9 "
[59] " [35] clusterProfiler_4.10.1 "
[60] " [36] zip_2.3.1 "
[61] " [37] webshot_0.5.5 "
[62] " [38] visNetwork_2.1.2 "
[63] " [39] viridis_0.6.5 "
[64] " [40] viridisLite_0.4.2 "
[65] " [41] Vennerable_3.1.0.9000 "
[66] " [42] UpSetR_1.4.0 "
[67] " [43] tinytex_0.50 "
[68] " [44] tidyr_1.3.1 "
[69] " [45] systemPipeShiny_1.12.0 "
[70] " [46] drawer_0.2.0.1 "
[71] " [47] stringdist_0.9.12 "
[72] " [48] spsUtil_0.2.2 "
[73] " [49] spsComps_0.3.3.0 "
[74] " [50] shinyWidgets_0.8.4 "
[75] " [51] shinyvalidate_0.1.3 "
[76] " [52] shinythemes_1.2.0 "
[77] " [53] shinymanager_1.0.410 "
[78] " [54] shinylogs_0.2.1 "
[79] " [55] shinyjs_2.1.0 "
[80] " [56] shinyfullscreen_1.1.0 "
[81] " [57] shinyFiles_0.9.3 "
[82] " [58] shinydashboardPlus_2.0.4 "
[83] " [59] shinydashboard_0.7.2 "
[84] " [60] shinyBS_0.61.1 "
[85] " [61] sf_1.0-16 "
[86] " [62] scales_1.3.0 "
[87] " [63] rmarkdown_2.26 "
[88] " [64] RecordLinkage_0.4-12.4 "
[89] " [65] ff_4.0.12 "
[90] " [66] bit_4.0.5 "
[91] " [67] RSQLite_2.3.6 "
[92] " [68] DBI_1.2.2 "
[93] " [69] readr_2.1.5 "
[94] " [70] r2symbols_1.4 "
[95] " [71] qpdf_1.3.3 "
[96] " [72] preprocessCore_1.64.0 "
[97] " [73] plyr_1.8.9 "
[98] " [74] plotrix_3.8-4 "
[99] " [75] plotly_4.10.4 "
[100] " [76] pathview_1.42.0 "
[101] " [77] manhattanly_0.3.0 "
[102] " [78] magrittr_2.0.3 "
[103] " [79] limma_3.58.1 "
[104] " [80] knitr_1.46 "
[105] " [81] JBrowseR_0.10.2 "
[106] " [82] illuminaio_0.44.0 "
[107] " [83] IlluminaHumanMethylationEPICv2manifest_1.0.0 "
[108] " [84] IlluminaHumanMethylationEPICv2anno.20a1.hg38_1.0.0 "
[109] " [85] minfi_1.48.0 "
[110] " [86] bumphunter_1.44.0 "
[111] " [87] locfit_1.5-9.9 "
[112] " [88] Biostrings_2.70.3 "
[113] " [89] XVector_0.42.0 "
[114] " [90] SummarizedExperiment_1.32.0 "
[115] " [91] MatrixGenerics_1.14.0 "
[116] " [92] matrixStats_1.3.0 "
[117] " [93] htmlwidgets_1.6.4 "
[118] " [94] gridExtra_2.3 "
[119] " [95] gprofiler2_0.2.3 "
[120] " [96] glue_1.7.0 "
[121] " [97] ggvenn_0.1.10 "
[122] " [98] ggrepel_0.9.5 "
[123] " [99] ggpubr_0.6.0 "
[124] "[100] ggfortify_0.4.16 "
[125] "[101] ggforce_0.4.2 "
[126] "[102] GenomicFeatures_1.54.4 "
[127] "[103] AnnotationDbi_1.64.1 "
[128] "[104] Biobase_2.62.0 "
[129] "[105] FactoMineR_2.10 "
[130] "[106] factoextra_1.0.7 "
[131] "[107] ggplot2_3.5.0 "
[132] "[108] explor_0.3.10 "
[133] "[109] DT_0.33 "
[134] "[110] dplyr_1.1.4 "
[135] "[111] doParallel_1.0.17 "
[136] "[112] iterators_1.0.14 "
[137] "[113] foreach_1.5.2 "
[138] "[114] data.table_1.15.4 "
[139] "[115] d3heatmap_0.9.0 "
[140] "[116] crosstalk_1.2.1 "
[141] "[117] corrplot_0.92 "
[142] "[118] colourpicker_1.3.0 "
[143] "[119] chromPlot_1.30.0 "
[144] "[120] biomaRt_2.58.2 "
[145] "[121] chromoMap_4.1.1 "
[146] "[122] ChAMPdata_2.34.0 "
[147] "[123] GenomicRanges_1.54.1 "
[148] "[124] GenomeInfoDb_1.38.8 "
[149] "[125] IRanges_2.36.0 "
[150] "[126] S4Vectors_0.40.2 "
[151] "[127] BiocGenerics_0.48.1 "
[152] "[128] Cairo_1.6-2 "
[153] "[129] BiocManager_1.30.22 "
[154] "[130] BBmisc_1.13 "
[155] "[131] shiny_1.8.1.1 "
[156] ""
[157] "loaded via a namespace (and not attached):"
[158] " [1] bslib_0.7.0 "
[159] " [2] httr_1.4.7 "
[160] " [3] gfonts_0.2.0 "
[161] " [4] biovizBase_1.50.0 "
[162] " [5] BiocParallel_1.36.0 "
[163] " [6] prettyunits_1.2.0 "
[164] " [7] ProtGenerics_1.34.0 "
[165] " [8] yulab.utils_0.1.4 "
[166] " [9] ggplotify_0.1.2 "
[167] " [10] lumi_2.54.0 "
[168] " [11] GenomicAlignments_1.38.2 "
[169] " [12] sparseMatrixStats_1.14.0 "
[170] " [13] pillar_1.9.0 "
[171] " [14] kableExtra_1.4.0 "
[172] " [15] Rgraphviz_2.46.0 "
[173] " [16] R6_2.5.1 "
[174] " [17] mime_0.12 "
[175] " [18] sysfonts_0.8.9 "
[176] " [19] edgeR_4.0.16 "
[177] " [20] Rttf2pt1_1.3.12 "
[178] " [21] reactable_0.4.4 "
[179] " [22] genefilter_1.84.0 "
[180] " [23] Rhdf5lib_1.24.2 "
[181] " [24] graphite_1.48.0 "
[182] " [25] Hmisc_5.1-2 "
[183] " [26] nor1mix_1.3-3 "
[184] " [27] parallelly_1.37.1 "
[185] " [28] kpmt_0.1.0 "
[186] " [29] caTools_1.18.2 "
[187] " [30] mgcv_1.9-0 "
[188] " [31] polyclip_1.10-6 "
[189] " [32] htmltools_0.5.8.1 "
[190] " [33] fansi_1.0.6 "
[191] " [34] showtext_0.9-7 "
[192] " [35] e1071_1.7-14 "
[193] " [36] fastICA_1.2-4 "
[194] " [37] anytime_0.3.9 "
[195] " [38] classInt_0.4-10 "
[196] " [39] car_3.1-2 "
[197] " [40] fgsea_1.28.0 "
[198] " [41] HDO.db_0.99.1 "
[199] " [42] scatterD3_1.0.1 "
[200] " [43] scatterplot3d_0.3-44 "
[201] " [44] clue_0.3-65 "
[202] " [45] rpart_4.1.21 "
[203] " [46] beanplot_1.3.1 "
[204] " [47] scatterpie_0.2.2 "
[205] " [48] tidyselect_1.2.1 "
[206] " [49] bsseq_1.38.0 "
[207] " [50] cowplot_1.1.3 "
[208] " [51] GenomeInfoDbData_1.2.11 "
[209] " [52] utf8_1.2.4 "
[210] " [53] fs_1.6.3 "
[211] " [54] future.apply_1.11.2 "
[212] " [55] graph_1.80.0 "
[213] " [56] R.oo_1.26.0 "
[214] " [57] ipred_0.9-14 "
[215] " [58] rtracklayer_1.62.0 "
[216] " [59] reactome.db_1.86.2 "
[217] " [60] JADE_2.0-4 "
[218] " [61] DelayedMatrixStats_1.24.0 "
[219] " [62] lazyeval_0.2.2 "
[220] " [63] sass_0.4.9 "
[221] " [64] carData_3.0-5 "
[222] " [65] munsell_0.5.1 "
[223] " [66] nanotime_0.3.7 "
[224] " [67] ROC_1.78.0 "
[225] " [68] ellipse_0.5.0 "
[226] " [69] lava_1.8.0 "
[227] " [70] treeio_1.24.3 "
[228] " [71] R.utils_2.12.3 "
[229] " [72] sva_3.50.0 "
[230] " [73] KEGGgraph_1.62.0 "
[231] " [74] bitops_1.0-7 "
[232] " [75] scrypt_0.1.6 "
[233] " [76] R.methodsS3_1.8.2 "
[234] " [77] labeling_0.4.3 "
[235] " [78] KEGGREST_1.42.0 "
[236] " [79] promises_1.3.0 "
[237] " [80] rhdf5filters_1.14.1 "
[238] " [81] zoo_1.8-12 "
[239] " [82] BiasedUrn_2.0.11 "
[240] " [83] globaltest_5.56.0 "
[241] " [84] DelayedArray_0.28.0 "
[242] " [85] marray_1.80.0 "
[243] " [86] assertthat_0.2.1 "
[244] " [87] tools_4.3.2 "
[245] " [88] ape_5.8 "
[246] " [89] ensembldb_2.26.0 "
[247] " [90] BiocFileCache_2.10.2 "
[248] " [91] rlang_1.1.3 "
[249] " [92] generics_0.1.3 "
[250] " [93] extrafont_0.19 "
[251] " [94] evaluate_0.23 "
[252] " [95] siggenes_1.76.0 "
[253] " [96] hrbrthemes_0.8.7 "
[254] " [97] BiocIO_1.12.0 "
[255] " [98] reshape2_1.4.4 "
[256] " [99] colorspace_2.1-0 "
[257] "[100] RBGL_1.78.0 "
[258] "[101] withr_3.0.0 "
[259] "[102] RCurl_1.98-1.14 "
[260] "[103] restfulr_0.0.15 "
[261] "[104] xtable_1.8-4 "
[262] "[105] aplot_0.2.2 "
[263] "[106] systemfonts_1.0.6 "
[264] "[107] mclust_6.1 "
[265] "[108] httpuv_1.6.15 "
[266] "[109] affyio_1.72.0 "
[267] "[110] latticeExtra_0.6-30 "
[268] "[111] units_0.8-5 "
[269] "[112] MASS_7.3-60 "
[270] "[113] reactR_0.5.0 "
[271] "[114] isva_1.9 "
[272] "[115] broom_1.0.5 "
[273] "[116] methylumi_2.48.0 "
[274] "[117] deldir_2.0-4 "
[275] "[118] rhdf5_2.46.1 "
[276] "[119] vctrs_0.6.5 "
[277] "[120] lifecycle_1.0.4 "
[278] "[121] ragg_1.3.0 "
[279] "[122] proxy_0.4-27 "
[280] "[123] codetools_0.2-19 "
[281] "[124] highr_0.10 "
[282] "[125] recipes_1.0.10 "
[283] "[126] nlme_3.1-163 "
[284] "[127] combinat_0.0-8 "
[285] "[128] future_1.33.2 "
[286] "[129] progress_1.2.3 "
[287] "[130] dbplyr_2.5.0 "
[288] "[131] jquerylib_0.1.4 "
[289] "[132] Rcpp_1.0.12 "
[290] "[133] rstudioapi_0.16.0 "
[291] "[134] patchwork_1.2.0 "
[292] "[135] stringi_1.8.3 "
[293] "[136] AnnotationFilter_1.26.0 "
[294] "[137] hms_1.1.3 "
[295] "[138] cachem_1.0.8 "
[296] "[139] multtest_2.58.0 "
[297] "[140] tidytree_0.4.6 "
[298] "[141] listenv_0.9.1 "
[299] "[142] geneLenDataBase_1.38.0 "
[300] "[143] jpeg_0.1-10 "
[301] "[144] ggtree_3.11.2 "
[302] "[145] enrichplot_1.22.0 "
[303] "[146] ExperimentHub_2.10.0 "
[304] "[147] ggfun_0.1.4 "
[305] "[148] HDF5Array_1.30.1 "
[306] "[149] SparseArray_1.2.4 "
[307] "[150] Formula_1.2-5 "
[308] "[151] interactiveDisplayBase_1.40.0 "
[309] "[152] leaps_3.1 "
[310] "[153] class_7.3-22 "
[311] "[154] memoise_2.0.1 "
[312] "[155] crayon_1.5.2 "
[313] "[156] gridGraphics_0.5-1 "
[314] "[157] rappdirs_0.3.3 "
[315] "[158] S4Arrays_1.2.1 "
[316] "[159] xml2_1.3.6 "
[317] "[160] filelock_1.0.3 "
[318] "[161] interp_1.1-6 "
[319] "[162] waiter_0.2.5 "
[320] "[163] png_0.1-8 "
[321] "[164] tzdb_0.4.0 "
[322] "[165] emmeans_1.10.1 "
[323] "[166] evd_2.3-6.1 "
[324] "[167] fastmap_1.1.1 "
[325] "[168] tidygraph_1.3.1 "
[326] "[169] flashClust_1.01-2 "
[327] "[170] pkgconfig_2.0.3 "
[328] "[171] cli_3.6.2 "
[329] "[172] affy_1.80.0 "
[330] "[173] RcppCCTZ_0.2.12 "
[331] "[174] prodlim_2023.08.28 "
[332] "[175] httpcode_0.3.0 "
[333] "[176] ggsignif_0.6.4 "
[334] "[177] nnet_7.3-19 "
[335] "[178] missMethyl_1.36.0 "
[336] "[179] BiocVersion_3.18.1 "
[337] "[180] textshaping_0.3.7 "
[338] "[181] multcompView_0.1-10 "
[339] "[182] gdtools_0.3.7 "
[340] "[183] timechange_0.3.0 "
[341] "[184] shinyAce_0.4.2 "
[342] "[185] ada_2.0-5 "
[343] "[186] foreign_0.8-85 "
[344] "[187] timeDate_4032.109 "
[345] "[188] splines_4.3.2 "
[346] "[189] askpass_1.2.0 "
[347] "[190] blob_1.2.4 "
[348] "[191] impute_1.76.0 "
[349] "[192] reshape_0.8.9 "
[350] "[193] annotate_1.80.0 "
[351] "[194] XML_3.99-0.16.1 "
[352] "[195] VariantAnnotation_1.48.1 "
[353] "[196] globals_0.16.3 "
[354] "[197] permute_0.9-7 "
[355] "[198] crul_1.4.2 "
[356] "[199] dichromat_2.0-0.1 "
[357] "[200] compiler_4.3.2 "
[358] "[201] rjson_0.2.21 "
[359] "[202] DNAcopy_1.76.0 "
[360] "[203] extrafontdb_1.0 "
[361] "[204] nleqslv_3.3.5 "
[362] "[205] AnnotationHub_3.10.1 "
[363] "[206] estimability_1.5 "
[364] "[207] formatR_1.14 "
[365] "[208] R.cache_0.16.0 "
[366] "[209] digest_0.6.35 "
[367] "[210] quadprog_1.5-8 "
[368] "[211] showtextdb_3.0 "
[369] "[212] graphlayouts_1.1.1 "
[370] "[213] base64_2.0.1 "
[371] "[214] wateRmelon_2.8.0 "
[372] "[215] fontLiberation_0.1.0 "
[373] "[216] vroom_1.6.5 "
[374] "[217] fontBitstreamVera_0.1.1 "
[375] "[218] zlibbioc_1.48.2 "
[376] "[219] tweenr_2.0.3 "
[377] "[220] ModelMetrics_1.2.2.2 "
[378] "[221] IlluminaHumanMethylationEPICanno.ilm10b4.hg19_0.6.0"
[379] "[222] fontawesome_0.5.2 "
[380] "[223] statmod_1.5.0 "
[381] "[224] ggraph_2.2.1 "
[382] "[225] openssl_2.1.1 "
[383] "[226] gson_0.1.0 "
[384] "[227] igraph_2.0.3 "
[385] "[228] mvtnorm_1.2-4 "
[386] "[229] yaml_2.3.8 "
[387] "[230] qvalue_2.34.0 "
[388] "[231] BSgenome_1.70.2 "
[389] "[232] later_1.3.2 "
[390] "[233] backports_1.4.1 "
[391] "[234] rstatix_0.7.2 "
[392] "[235] pagedown_0.20 "
[393] "[236] shadowtext_0.1.3 "
[394] "[237] Rsamtools_2.18.0 "
[395] "[238] fontquiver_0.2.1 "
[396] "[239] miniUI_0.1.1.1 "
[397] "[240] gtable_0.3.4 "
[398] "[241] abind_1.4-5 "
[399] "[242] xfun_0.43 "
[400] "[243] curl_5.2.1 "
[401] "[244] GEOquery_2.70.0 "
[402] "[245] KernSmooth_2.23-22 "
[403] "[246] styler_1.10.3 "
[404] "[247] survival_3.5-7 "
[405] "[248] goseq_1.54.0 "
[406] "[249] jsonlite_1.8.8 "
[407] "[250] shinyjqui_0.4.1 "
[408] "[251] svglite_2.1.3 "
[409] "[252] base64enc_0.1-3 "
[410] "[253] scrime_1.3.5 "
[411] "[254] billboarder_0.4.1 "
[412] "[255] Matrix_1.6-1.1 "
[413] "[256] Gviz_1.46.1 "
[414] "[257] shinytoastr_2.2.0 "
[415] "[258] fastmatch_1.1-4 "
[416] "[259] gower_1.0.1 "
[417] "[260] checkmate_2.3.1 "
[418] "[261] hardhat_1.3.1 "
[419] "[262] gtools_3.9.5 "
[420] "[263] htmlTable_2.4.2 "
[421] "[264] bit64_4.0.5 "
[422] "[265] farver_2.1.1 "
Other Information
Webservers:
@LiU
- https://methylr.research.liu.se
@SciLifeLab data center
- https://methylr.serve.scilifelab.se
Manual: https://methylr.netlify.app/
GitHub: https://github.com/JD2112/methylr
Docker: https://hub.docker.com/r/jd21/methylr
Google Group: https://groups.google.com/g/methylr
Test data: https://sourceforge.net/projects/methylr/
Contact: methylr@googlegroups.com
Acknowledgements
Linköping University IT Department
Core Facility, Linköping University
Clinical Genomics Linköping
NAISS and UPPMAX center, Sweden
SciLifeLab Data Center
© Jyotirmoy Das