0.1 Introduction

The bodymapRat package contains gene expression data on 652 RNA-Seq samples from a comprehensive rat transcriptomic BodyMap study. These samples include the sequence identifier information provided in the header of the FASTQ files which can be used as a surrogate for batch. These samples have not been normalized or pre-processed.

The data are provided in a SummarizedExperiment. The phenotypic information can be extracted using the colData() function and a description of the phenotypic data is listed in the table below:

Title Description
sraExperiment SRA Experiment ID
title Title of sample provided by the authors
geoAccession GEO Accession ID
BioSample BioSample ID
avgLength Average read length
instrument Machine identifier (from FASTQ header)
runID Run ID (from FASTQ header)
fcID Flow cell ID (from FASTQ header)
fcLane Flow cell lane (from FASTQ header)
tile Tile (from FASTQ header)
xtile xtile (from FASTQ header)
ytile ytile (from FASTQ header)
organ Body organ
sex Gender
stage Stage
techRep Technical replicate number
colOrgan Column of colors to help with plotting
rnaRIN RIN number
barcode barcode number

The data can be accessed as follows:

library(SummarizedExperiment) 
library(bodymapRat)

We use the bodymapRat() function to download the relevant files from Bioconductor’s ExperimentHub web resource. Running this function will download a SummarizedExperiment object, which contains read counts, as well as the metadata on the rows (genes) and columns (cells).

bm_rat <- bodymapRat()

# Get the expression data
counts = assay(bm_rat)
dim(counts)
## [1] 32637   652
counts[1:5, 1:5]
##                    SRR1169893 SRR1169894 SRR1169895 SRR1169896 SRR1169897
## ENSRNOG00000000001          1          0          0          1          4
## ENSRNOG00000000007          1          1          0          3          0
## ENSRNOG00000000008          7          4          2          3          7
## ENSRNOG00000000009          0          0          0          0          1
## ENSRNOG00000000010          0          1          0          0          0
# Get the meta data along columns
head(colData(bm_rat))
## DataFrame with 6 rows and 22 columns
##            sraExperiment      sraRun       title geoAccession sraSample
##              <character> <character> <character>     <factor>  <factor>
## SRR1169893     SRX471368  SRR1169893 Adr_F_002_1   GSM1328469 SRS558114
## SRR1169894     SRX471368  SRR1169894 Adr_F_002_1   GSM1328469 SRS558114
## SRR1169895     SRX471369  SRR1169895 Adr_F_002_2   GSM1328470 SRS558115
## SRR1169896     SRX471369  SRR1169896 Adr_F_002_2   GSM1328470 SRS558115
## SRR1169897     SRX471370  SRR1169897 Adr_F_002_3   GSM1328471 SRS558116
## SRR1169898     SRX471370  SRR1169898 Adr_F_002_3   GSM1328471 SRS558116
##               BioSample avgLength       organ         sex     stage   techRep
##                <factor> <integer> <character> <character> <numeric> <integer>
## SRR1169893 SAMN02642886        50     Adrenal           F         2         1
## SRR1169894 SAMN02642886        50     Adrenal           F         2         2
## SRR1169895 SAMN02642867        50     Adrenal           F         2         1
## SRR1169896 SAMN02642867        50     Adrenal           F         2         2
## SRR1169897 SAMN02642894        50     Adrenal           F         2         1
## SRR1169898 SAMN02642894        50     Adrenal           F         2         2
##               colOrgan         mix      rnaRIN     barcode  instrument
##            <character> <character> <character> <character> <character>
## SRR1169893       brown          M1         9.3          11   HWI-ST845
## SRR1169894       brown          M1         9.3          11   HWI-ST845
## SRR1169895       brown          M1         9.1           5   HWI-ST845
## SRR1169896       brown          M1         9.1           5   HWI-ST845
## SRR1169897       brown          M1         9.5           3  HWI-ST1131
## SRR1169898       brown          M1         9.5           3  HWI-ST1195
##                  runID        fcID      fcLane        tile       xtile
##            <character> <character> <character> <character> <character>
## SRR1169893      120326   D0VTJACXX           2        1101        1506
## SRR1169894      120525   D10G7ACXX           2        1101        1394
## SRR1169895      120326   D0VTJACXX           5        1101        1170
## SRR1169896      120525   D10G7ACXX           5        1101        1650
## SRR1169897      120424   C0P4UACXX           4        1101        1675
## SRR1169898      120525   C0TDUACXX           4        1101        1138
##                  ytile
##            <character>
## SRR1169893        2000
## SRR1169894        2133
## SRR1169895        2029
## SRR1169896        2126
## SRR1169897        2216
## SRR1169898        2067

The data in this package are used as an example data set in the qsmooth Bioconductor package.

1 References

Appendix

  1. Yu et al. (2013). A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages. Nature Communications 5:3230. PMID: 24510058. PMCID: PMC3926002.

A SessionInfo

sessionInfo()
## 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] bodymapRat_1.18.0           ExperimentHub_2.10.0       
##  [3] AnnotationHub_3.10.0        BiocFileCache_2.10.0       
##  [5] dbplyr_2.3.4                SummarizedExperiment_1.32.0
##  [7] Biobase_2.62.0              GenomicRanges_1.54.0       
##  [9] GenomeInfoDb_1.38.0         IRanges_2.36.0             
## [11] S4Vectors_0.40.0            BiocGenerics_0.48.0        
## [13] MatrixGenerics_1.14.0       matrixStats_1.0.0          
## [15] knitr_1.44                  BiocStyle_2.30.0           
## 
## loaded via a namespace (and not attached):
##  [1] KEGGREST_1.42.0               xfun_0.40                    
##  [3] bslib_0.5.1                   lattice_0.22-5               
##  [5] vctrs_0.6.4                   tools_4.3.1                  
##  [7] bitops_1.0-7                  generics_0.1.3               
##  [9] curl_5.1.0                    AnnotationDbi_1.64.0         
## [11] tibble_3.2.1                  fansi_1.0.5                  
## [13] RSQLite_2.3.1                 blob_1.2.4                   
## [15] pkgconfig_2.0.3               Matrix_1.6-1.1               
## [17] lifecycle_1.0.3               GenomeInfoDbData_1.2.11      
## [19] compiler_4.3.1                Biostrings_2.70.1            
## [21] httpuv_1.6.12                 htmltools_0.5.6.1            
## [23] sass_0.4.7                    RCurl_1.98-1.12              
## [25] yaml_2.3.7                    interactiveDisplayBase_1.40.0
## [27] later_1.3.1                   pillar_1.9.0                 
## [29] crayon_1.5.2                  jquerylib_0.1.4              
## [31] ellipsis_0.3.2                DelayedArray_0.28.0          
## [33] cachem_1.0.8                  abind_1.4-5                  
## [35] mime_0.12                     tidyselect_1.2.0             
## [37] digest_0.6.33                 purrr_1.0.2                  
## [39] dplyr_1.1.3                   bookdown_0.36                
## [41] BiocVersion_3.18.0            fastmap_1.1.1                
## [43] grid_4.3.1                    cli_3.6.1                    
## [45] SparseArray_1.2.0             magrittr_2.0.3               
## [47] S4Arrays_1.2.0                utf8_1.2.4                   
## [49] withr_2.5.1                   promises_1.2.1               
## [51] filelock_1.0.2                rappdirs_0.3.3               
## [53] bit64_4.0.5                   rmarkdown_2.25               
## [55] XVector_0.42.0                httr_1.4.7                   
## [57] bit_4.0.5                     png_0.1-8                    
## [59] shiny_1.7.5.1                 memoise_2.0.1                
## [61] evaluate_0.22                 rlang_1.1.1                  
## [63] Rcpp_1.0.11                   xtable_1.8-4                 
## [65] glue_1.6.2                    DBI_1.1.3                    
## [67] BiocManager_1.30.22           jsonlite_1.8.7               
## [69] R6_2.5.1                      zlibbioc_1.48.0