RTCGA package

The Cancer Genome Atlas (TCGA) Data Portal provides a platform for researchers to search, download, and analyze data sets generated by TCGA. It contains clinical information, genomic characterization data, and high level sequence analysis of the tumor genomes. The key is to understand genomics to improve cancer care.

RTCGA package offers download and integration of the variety and volume of TCGA data using patient barcode key, what enables easier data possession. This may have a benefcial infuence on development of science and improvement of patients’ treatment. RTCGA is an open-source R package, available to download from Bioconductor

source("http://bioconductor.org/biocLite.R")
biocLite("RTCGA")

or from GitHub

biocLite("RTCGA/RTCGA")

Furthermore, RTCGA package transforms TCGA data into form which is convenient to use in R statistical package. Those data transformations can be a part of statistical analysis pipeline which can be more reproducible with RTCGA.

Use cases and examples are shown in RTCGA packages vignettes:

browseVignettes("RTCGA")

How to download methylation data to gain the same datasets as in RTCGA.methylation package?

There are many available date times of TCGA data releases. To see them all just type:

library(RTCGA)
checkTCGA('Dates')

Version 1.0 of RTCGA.methylation package contains methylation datasets which were released 2015-11-01. They were downloaded in the following way (which is mainly copied from http://rtcga.github.io/RTCGA/:

Available cohorts

All cohort names can be checked using:

(cohorts <- infoTCGA() %>% 
   rownames() %>% 
   sub("-counts", "", x=.))

For all cohorts the following code downloads the mRNA data.

Downloading tarred files

#dir.create( "data2" )
releaseDate <- "2015-11-01"
sapply(cohorts,
             function(element){
                try({
downloadTCGA( cancerTypes = element, 
              dataSet = "Merge_methylation__humanmethylation27",
              destDir = "data2", 
              date = releaseDate )
                })
})

Reading downloaded methylation dataset

Shortening paths and directories

list.files( "data2") %>% 
   file.path( "data2", .) %>%
   file.rename( to = substr(.,start=1,stop=50))

Removing `NA` files from data2

The existance of NA files mean that there were no methylation data for these cohorts.

list.files( "data2") %>%
   file.path( "data2", .) %>%
   sapply(function(x){
      if (x == "data2/NA")
         file.remove(x)      
   })

Paths to methylation data

Below is the code that removes unneeded “MANIFEST.txt” file from each methylation cohort folder.

list.files( "data2") %>% 
   file.path( "data2", .) %>%
   sapply(function(x){
      file.path(x, list.files(x)) %>%
         grep(pattern = "MANIFEST.txt", x = ., value=TRUE) %>%
         file.remove()
      })

Below is the code that automatically gives the path to files for all available methylation cohorts types downloaded to data2 folder.

list.files("data2") %>%
   file.path("data2", .) %>%
   sapply(function(y){
      file.path(y, list.files(y)) %>%
      assign(value = .,
              x = paste0(list.files(y) %>%
                           gsub(x = .,
                                pattern = "\\..*",
                                replacement = "") %>%
                            gsub(x=.,
                                     pattern="-",
                                     replacement = "_"),
                         ".methylation.path"),
              envir = .GlobalEnv)
   })

Reading methylation data using `readTCGA`

Because of the fact that methylation data are transposed in downloaded files, there has been prepared special function readTCGA to read and transpose data automatically. Code is below

ls() %>%
   grep("methylation\\.path", x = ., value = TRUE) %>% 
   sapply(function(element){
      try({
         readTCGA(get(element, envir = .GlobalEnv),
               dataType = "methylation") %>%
         assign(value = .,
                x = sub("\\.path", "", x = element),
                envir = .GlobalEnv )
      }) 
     invisible(NULL)
    })

Using `RTCGA` package to download methylation data that are included in `RTCGA.methylation` package

Date of datasets release: 2015-11-01

Marcin Kosiński

2017-10-31

RTCGA package

How to download methylation data to gain the same datasets as in RTCGA.methylation package?

Available cohorts

Downloading tarred files

Reading downloaded methylation dataset

Shortening paths and directories

Removing `NA` files from data2

Paths to methylation data

Reading methylation data using `readTCGA`

Saving methylation data to `RTCGA.methylation` package

Using RTCGA package to download methylation data that are included in RTCGA.methylation package

Date of datasets release: 2015-11-01

Marcin Kosiński

2017-10-31

RTCGA package

How to download methylation data to gain the same datasets as in RTCGA.methylation package?

Available cohorts

Downloading tarred files

Reading downloaded methylation dataset

Shortening paths and directories

Removing NA files from data2

Paths to methylation data

Reading methylation data using readTCGA

Saving methylation data to RTCGA.methylation package

Using `RTCGA` package to download methylation data that are included in `RTCGA.methylation` package

Removing `NA` files from data2

Reading methylation data using `readTCGA`

Saving methylation data to `RTCGA.methylation` package