library(TumourMethData)
#> Warning: replacing previous import 'HDF5Array::h5ls' by 'rhdf5::h5ls' when
#> loading 'TumourMethData'DNA methylation is a repressive epigenetic modification involving the addition of methyl groups to DNA and occurs almost exclusively at CpG dinucleotides in mammals. Altered DNA methylation plays a profound role in the development and progression of cancer. However, much of our knowledge of DNA methylation in cancer has been garnered from methylation microarrays which measure methylation at only a small subset (generally <1%) of the almost 30 million CpG sites in humans, mostly those located close to gene promoters. Thus, whole genome bisulfite sequenicng (WGBS) studies in tumours which measure DNA methylation across the entire genome provide an invaluable resource for gaining a comprehensive understanding of DNA methylation changes in cancer, especially at regulatory regions located far from genes.
While packages such as curatedTCGAData provide DNA methylation data generated with microarrays for a range of different cancer types,
TumourMethData provides a collection of whole genome DNA methylation datasets for several different cancers (primary prostate cancer, prostate cancer metastases, esophageal cancer and rhabdoid tumour at present) as well as matching normal samples where available.
These whole genome methylation datasets are provided as RangedSummarizedExperiments, facilitating easy download of the data and extraction of methylation values for regions of interest.
Furthermore, RNA-seq transcripts counts are also provided for several of the datasets, enabling thorough analysis of how DNA methylation is associated with transcription and how this relationship is perturbed in cancer.
We can view the available datasets with TumourMethDatasets.
# Show available methylation datasets
data("TumourMethDatasets", package = "TumourMethData")
print(TumourMethDatasets)
#> dataset_name cancer_type technology genome_build
#> 1 cpgea_wgbs_hg38 prostate WGBS hg38
#> 2 tcga_wgbs_hg38 various WGBS hg38
#> 3 mcrpc_wgbs_hg38 prostate WGBS hg38
#> 4 mcrpc_wgbs_hg38_chr11 prostate WGBS hg38
#> 5 cao_esophageal_wgbs_hg19 esophageal WGBS hg19
#> 6 target_rhabdoid_wgbs_hg19 rhabdoid WGBS hg19
#> number_tumour_samples number_normal_samples wgbs_coverage_available
#> 1 187 187 FALSE
#> 2 39 8 FALSE
#> 3 100 0 TRUE
#> 4 100 0 TRUE
#> 5 10 9 FALSE
#> 6 69 0 FALSE
#> dataset_size_gb transcript_counts_available
#> 1 40.00 TRUE
#> 2 5.40 TRUE
#> 3 16.00 TRUE
#> 4 0.76 TRUE
#> 5 2.00 TRUE
#> 6 4.50 TRUE
#> notes
#> 1
#> 2
#> 3
#> 4 This dataset is a subset of the data in mcrpc_wgbs_hg38 for example purposes
#> 5
#> 6 Methylation values are not as precise as in other datasets. The original \n methylation values were integers between 0 and 10 with separate values for the C and G positions of each CpG site.\n The mean of these values was divided by 10 to produce the methylation values here, \n with CpG sites missing methylation values for either to C or G given an NA value
#> original_publication
#> 1 A genomic and epigenomic atlas of prostate cancer in Asian populations; Nature; 2020
#> 2 DNA methylation loss in late-replicating domains is linked to mitotic cell division; Nature genetics; 2018
#> 3 The DNA methylation landscape of advanced prostate cancer; Nature genetics; 2020
#> 4 The DNA methylation landscape of advanced prostate cancer; Nature genetics; 2020
#> 5 Multi-faceted epigenetic dysregulation of gene expression promotes esophageal squamous cell carcinoma; Nature communications; 2020
#> 6 Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways; Cancer Cell; 2016We use download_meth_dataset to download the methylation dataset we are interested in using mcrpc_wgbs_hg38_chr11 as an example.
# Download esophageal WGBS data
mcrpc_wgbs_hg38_chr11 = download_meth_dataset(dataset = "mcrpc_wgbs_hg38_chr11")
#> see ?TumourMethData and browseVignettes('TumourMethData') for documentation
#> loading from cache
#> require("rhdf5")
print(mcrpc_wgbs_hg38_chr11)
#> class: RangedSummarizedExperiment
#> dim: 1333114 100
#> metadata(5): genome is_h5 ref_CpG chrom_sizes descriptive_stats
#> assays(2): beta cov
#> rownames: NULL
#> rowData names(0):
#> colnames(100): DTB_003 DTB_005 ... DTB_265 DTB_266
#> colData names(4): metastatis_site subtype age sexsessionInfo()
#> R version 4.3.1 (2023-06-16)
#> Platform: x86_64-pc-linux-gnu (64-bit)
#> Running under: Ubuntu 22.04.3 LTS
#>
#> Matrix products: default
#> BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
#> LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
#>
#> 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
#>
#> time zone: America/New_York
#> tzcode source: system (glibc)
#>
#> attached base packages:
#> [1] stats4 stats graphics grDevices utils datasets methods
#> [8] base
#>
#> other attached packages:
#> [1] HDF5Array_1.30.0 rhdf5_2.46.0
#> [3] TumourMethData_1.0.0 SummarizedExperiment_1.32.0
#> [5] Biobase_2.62.0 GenomicRanges_1.54.0
#> [7] GenomeInfoDb_1.38.0 IRanges_2.36.0
#> [9] S4Vectors_0.40.0 BiocGenerics_0.48.0
#> [11] MatrixGenerics_1.14.0 matrixStats_1.0.0
#>
#> loaded via a namespace (and not attached):
#> [1] tidyselect_1.2.0 dplyr_1.1.3
#> [3] blob_1.2.4 R.utils_2.12.2
#> [5] filelock_1.0.2 Biostrings_2.70.1
#> [7] bitops_1.0-7 fastmap_1.1.1
#> [9] RCurl_1.98-1.12 BiocFileCache_2.10.0
#> [11] promises_1.2.1 digest_0.6.33
#> [13] mime_0.12 lifecycle_1.0.3
#> [15] ellipsis_0.3.2 KEGGREST_1.42.0
#> [17] interactiveDisplayBase_1.40.0 RSQLite_2.3.1
#> [19] magrittr_2.0.3 compiler_4.3.1
#> [21] rlang_1.1.1 sass_0.4.7
#> [23] tools_4.3.1 utf8_1.2.4
#> [25] yaml_2.3.7 knitr_1.44
#> [27] S4Arrays_1.2.0 bit_4.0.5
#> [29] curl_5.1.0 DelayedArray_0.28.0
#> [31] abind_1.4-5 withr_2.5.1
#> [33] purrr_1.0.2 R.oo_1.25.0
#> [35] grid_4.3.1 fansi_1.0.5
#> [37] ExperimentHub_2.10.0 xtable_1.8-4
#> [39] Rhdf5lib_1.24.0 cli_3.6.1
#> [41] rmarkdown_2.25 crayon_1.5.2
#> [43] generics_0.1.3 httr_1.4.7
#> [45] DBI_1.1.3 cachem_1.0.8
#> [47] zlibbioc_1.48.0 AnnotationDbi_1.64.0
#> [49] BiocManager_1.30.22 XVector_0.42.0
#> [51] vctrs_0.6.4 Matrix_1.6-1.1
#> [53] jsonlite_1.8.7 bit64_4.0.5
#> [55] jquerylib_0.1.4 glue_1.6.2
#> [57] BiocVersion_3.18.0 later_1.3.1
#> [59] tibble_3.2.1 pillar_1.9.0
#> [61] rappdirs_0.3.3 htmltools_0.5.6.1
#> [63] rhdf5filters_1.14.0 GenomeInfoDbData_1.2.11
#> [65] R6_2.5.1 dbplyr_2.3.4
#> [67] evaluate_0.22 shiny_1.7.5.1
#> [69] lattice_0.22-5 AnnotationHub_3.10.0
#> [71] R.methodsS3_1.8.2 png_0.1-8
#> [73] memoise_2.0.1 httpuv_1.6.12
#> [75] bslib_0.5.1 Rcpp_1.0.11
#> [77] SparseArray_1.2.0 xfun_0.40
#> [79] pkgconfig_2.0.3