In this section, we will perform the same analysis performed using ELMER, but instead of doing it programmatically we will use TCGAbiolinksGUI (Silva et al. 2017).
First we will launch the TCGAbiolinksGUI.
Please download this two objects:
To create the MultiAssayExperiment object go to Integrative analysis/ELMER/Create input data
.
Select the DNA methylation object previously downloaded.
Select the gene expression object previously downloaded.
Fill the field Save as:
and click on Create MAE object.
The object will be created.
To perform ELMER analysis go to Integrative analysis/ELMER/Analysis
.
Select the MAE data created in the previous section.
Select the groups that will be analysed: Primary solid Tumor and Solid Tissue Normal.
We will identify probes that are hypomethylated in Primary solid Tumor compared to Solid Tissue Normal.
For the significant differently methylated probes identified before we will correlated with the 20 nearest genes. Change the value of the field Number of permutations
to 100
, Raw P-value cut-off
to 0.05
and Empirical P value cut-off
to 0.01
.
There will be no changes in the step 3.
There will be no changes in the step 4.
Click on Run the analysis
.
If the analysis identified significant regulatory TF the results will be saved into an R object.
To visualize the results go to Integrative analysis/ELMER/Visualize results
.
Click on Select results
and select the object created on the previous section.
Or the avarage DNA methylation levels of probes of a Motif vs the expression of a TF.
For each enriched motif you can verify the ranking of sigificances between the correlation of DNA methylation level on the significant paired probes with that motif vs the TF expression (for all human TF).
The enrichement of each motif can be visualized.
You can take a look for a gene which was the probe linked.
You can see the plot and its neraby genes.
It is possible to visualize the table with the significant differently methylated probes.
It is possible to visualize the table with the enriched motifs.
It is possible to visualize the table with the candidates regulatory TF.
## R version 4.1.1 (2021-08-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.14-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.14-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [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] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] MultiAssayExperiment_1.20.0 SummarizedExperiment_1.24.0
## [3] Biobase_2.54.0 GenomicRanges_1.46.0
## [5] GenomeInfoDb_1.30.0 IRanges_2.28.0
## [7] S4Vectors_0.32.0 MatrixGenerics_1.6.0
## [9] matrixStats_0.61.0 sesameData_1.11.13
## [11] rmarkdown_2.11 ExperimentHub_2.2.0
## [13] AnnotationHub_3.2.0 BiocFileCache_2.2.0
## [15] dbplyr_2.1.1 BiocGenerics_0.40.0
## [17] BiocStyle_2.22.0 dplyr_1.0.7
## [19] DT_0.19 ELMER_2.18.0
## [21] ELMER.data_2.17.1
##
## loaded via a namespace (and not attached):
## [1] utf8_1.2.2 R.utils_2.11.0
## [3] tidyselect_1.1.1 RSQLite_2.2.8
## [5] AnnotationDbi_1.56.0 htmlwidgets_1.5.4
## [7] grid_4.1.1 BiocParallel_1.28.0
## [9] munsell_0.5.0 ragg_1.1.3
## [11] codetools_0.2-18 withr_2.4.2
## [13] colorspace_2.0-2 filelock_1.0.2
## [15] highr_0.9 knitr_1.36
## [17] rstudioapi_0.13 ggsignif_0.6.3
## [19] labeling_0.4.2 GenomeInfoDbData_1.2.7
## [21] farver_2.1.0 bit64_4.0.5
## [23] downloader_0.4 vctrs_0.3.8
## [25] generics_0.1.1 xfun_0.27
## [27] biovizBase_1.42.0 R6_2.5.1
## [29] doParallel_1.0.16 clue_0.3-60
## [31] AnnotationFilter_1.18.0 bitops_1.0-7
## [33] cachem_1.0.6 reshape_0.8.8
## [35] DelayedArray_0.20.0 assertthat_0.2.1
## [37] vroom_1.5.5 promises_1.2.0.1
## [39] BiocIO_1.4.0 scales_1.1.1
## [41] nnet_7.3-16 gtable_0.3.0
## [43] ensembldb_2.18.0 rlang_0.4.12
## [45] systemfonts_1.0.3 GlobalOptions_0.1.2
## [47] splines_4.1.1 rtracklayer_1.54.0
## [49] rstatix_0.7.0 lazyeval_0.2.2
## [51] dichromat_2.0-0 broom_0.7.9
## [53] checkmate_2.0.0 BiocManager_1.30.16
## [55] yaml_2.2.1 reshape2_1.4.4
## [57] abind_1.4-5 crosstalk_1.1.1
## [59] GenomicFeatures_1.46.0 backports_1.2.1
## [61] httpuv_1.6.3 Hmisc_4.6-0
## [63] tools_4.1.1 ggplot2_3.3.5
## [65] ellipsis_0.3.2 jquerylib_0.1.4
## [67] RColorBrewer_1.1-2 Rcpp_1.0.7
## [69] plyr_1.8.6 base64enc_0.1-3
## [71] progress_1.2.2 zlibbioc_1.40.0
## [73] purrr_0.3.4 RCurl_1.98-1.5
## [75] prettyunits_1.1.1 ggpubr_0.4.0
## [77] rpart_4.1-15 GetoptLong_1.0.5
## [79] haven_2.4.3 ggrepel_0.9.1
## [81] cluster_2.1.2 magrittr_2.0.1
## [83] data.table_1.14.2 openxlsx_4.2.4
## [85] circlize_0.4.13 ProtGenerics_1.26.0
## [87] xtable_1.8-4 mime_0.12
## [89] hms_1.1.1 TCGAbiolinksGUI.data_1.13.1
## [91] evaluate_0.14 XML_3.99-0.8
## [93] rio_0.5.27 jpeg_0.1-9
## [95] readxl_1.3.1 gridExtra_2.3
## [97] shape_1.4.6 compiler_4.1.1
## [99] biomaRt_2.50.0 tibble_3.1.5
## [101] crayon_1.4.1 R.oo_1.24.0
## [103] htmltools_0.5.2 later_1.3.0
## [105] tzdb_0.1.2 Formula_1.2-4
## [107] tidyr_1.1.4 DBI_1.1.1
## [109] ComplexHeatmap_2.10.0 rappdirs_0.3.3
## [111] Matrix_1.3-4 car_3.0-11
## [113] readr_2.0.2 cli_3.0.1
## [115] R.methodsS3_1.8.1 parallel_4.1.1
## [117] Gviz_1.38.0 forcats_0.5.1
## [119] pkgconfig_2.0.3 GenomicAlignments_1.30.0
## [121] foreign_0.8-81 plotly_4.10.0
## [123] xml2_1.3.2 foreach_1.5.1
## [125] bslib_0.3.1 XVector_0.34.0
## [127] rvest_1.0.2 stringr_1.4.0
## [129] VariantAnnotation_1.40.0 digest_0.6.28
## [131] Biostrings_2.62.0 cellranger_1.1.0
## [133] TCGAbiolinks_2.22.0 htmlTable_2.3.0
## [135] restfulr_0.0.13 curl_4.3.2
## [137] shiny_1.7.1 Rsamtools_2.10.0
## [139] rjson_0.2.20 lifecycle_1.0.1
## [141] jsonlite_1.7.2 carData_3.0-4
## [143] viridisLite_0.4.0 BSgenome_1.62.0
## [145] fansi_0.5.0 pillar_1.6.4
## [147] lattice_0.20-45 KEGGREST_1.34.0
## [149] fastmap_1.1.0 httr_1.4.2
## [151] survival_3.2-13 interactiveDisplayBase_1.32.0
## [153] glue_1.4.2 zip_2.2.0
## [155] png_0.1-7 iterators_1.0.13
## [157] BiocVersion_3.14.0 bit_4.0.4
## [159] stringi_1.7.5 sass_0.4.0
## [161] blob_1.2.2 textshaping_0.3.6
## [163] latticeExtra_0.6-29 memoise_2.0.0
Silva, Tiago C., Antonio Colaprico, Catharina Olsen, Gianluca Bontempi, Michele Ceccarelli, Benjamin P. Berman, and Houtan Noushmehr. 2017. “TCGAbiolinksGUI: A Graphical User Interface to Analyze Cancer Molecular and Clinical Data.” bioRxiv. https://doi.org/10.1101/147496.