Fluent genomics with plyranges and tximeta
Author(s)
Lee, S; Lawrence, M; Love, MI;
Journal Title
F1000Research
Publication Type
Journal Article
Abstract
We construct a simple workflow for fluent genomics data analysis using the R/Bioconductor ecosystem. This involves three core steps: import the data into an appropriate abstraction, model the data with respect to the biological questions of interest, and integrate the results with respect to their underlying genomic coordinates. Here we show how to implement these steps to integrate published RNA-seq and ATAC-seq experiments on macrophage cell lines. Using tximeta, we import RNA-seq transcript quantifications into an analysis-ready data structure, called the SummarizedExperiment, that contains the ranges of the reference transcripts and metadata on their provenance. Using SummarizedExperiments to represent the ATAC-seq and RNA-seq data, we model differentially accessible (DA) chromatin peaks and differentially expressed (DE) genes with existing Bioconductor packages. Using plyranges we then integrate the results to see if there is an enrichment of DA peaks near DE genes by finding overlaps and aggregating over log-fold change thresholds. The combination of these packages and their integration with the Bioconductor ecosystem provide a coherent framework for analysts to iteratively and reproducibly explore their biological data.
Publisher
F1000
WEHI Research Division(s)
Structural Biology
PubMed ID
32528659
Open Access at Publisher's Site
https://doi.org/10.12688/f1000research.22259.1
Rights Notice
Refer to copyright notice on published article.


Creation Date: 2020-06-15 10:33:57
Last Modified: 2020-06-15 10:39:46
An error has occurred. This application may no longer respond until reloaded. Reload 🗙