Using the GEOfastq Package

library(GEOfastq)

Installation

GEOfastq can be installed from Bioconductor as follows:

if(!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("GEOfastq")

Overview of GEOfastq

The NCBI Gene Expression Omnibus (GEO) offers a convenient interface to explore high-throughput experimental data such as RNA-seq. GEO deposits RNA-seq data as sra files to the Sequence Read Archive (SRA) which can be converted to fastq files using fastq-dump. This conversion process can be quite slow and it is usually more convenient to download fastq files for a GEO accession generated by the European Nucleotide Archive (ENA). GEOfastq crawls GEO to retrieve metadata and ENA fastq urls, and then downloads them.

Getting Started using GEOfastq

To get fastq data for a GEO series, we first retrieve the metadata for a GEO accession:

gse_name <- 'GSE133758'
gse_text <- crawl_gse(gse_name)

Next, we extract the sample accessions for this study and retrieve the GEO metadata and ENA fastq url for an example:

gsm_names <- extract_gsms(gse_text)
gsm_name <- gsm_names[182]
srp_meta <- crawl_gsms(gsm_name)
#> 1 GSMs to process

Now that we have retrieved the necessary metadata, we are ready to download the fastq files for this sample:

data_dir <- tempdir()

# example using smaller file
srp_meta <- data.frame(
        run  = 'SRR014242',
        row.names = 'SRR014242',
        gsm_name = 'GSM315559',
        ebi_dir = get_dldir('SRR014242'), stringsAsFactors = FALSE)

res <- get_fastqs(srp_meta, data_dir)
#> Warning: Recv failure: Connection reset by peer

Session info

The following package and versions were used in the production of this vignette.

#> 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] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> other attached packages:
#> [1] GEOfastq_1.2.0
#> 
#> loaded via a namespace (and not attached):
#>  [1] Rcpp_1.0.7        knitr_1.36        magrittr_2.0.1    doParallel_1.0.16
#>  [5] R6_2.5.1          rlang_0.4.12      fastmap_1.1.0     foreach_1.5.1    
#>  [9] stringr_1.4.0     plyr_1.8.6        tools_4.1.1       parallel_4.1.1   
#> [13] xfun_0.27         jquerylib_0.1.4   htmltools_0.5.2   iterators_1.0.13 
#> [17] yaml_2.2.1        digest_0.6.28     sass_0.4.0        bitops_1.0-7     
#> [21] codetools_0.2-18  RCurl_1.98-1.5    evaluate_0.14     rmarkdown_2.11   
#> [25] stringi_1.7.5     compiler_4.1.1    bslib_0.3.1       jsonlite_1.7.2