Upgrade to Pro — share decks privately, control downloads, hide ads and more …

recount-webinar

 recount-webinar

recount workflow: Accessing over 70,000 human RNA-seq samples with Bioconductor

7382f7fe30561274624635116513ca37?s=128

Leonardo Collado-Torres

April 17, 2018
Tweet

Transcript

  1. 11 recount workflow: Accessing over 70,000 human RNA-seq samples with

    Bioconductor Leonardo Collado-Torres @fellgernon sites.google.com/view/lalresearchgroup/our-weekly-webinar speakerdeck.com/lcolladotor/recount-webinar April 17, 2018
  2. Reference genome Reads f1000research.com/articles/6-1558/v1

  3. f1000research.com/articles/6-1558/v1

  4. GTEx TCGA slide adapted from Shannon Ellis

  5. SRA

  6. http://rail.bio/ Slide adapted from Ben Langmead

  7. http://blogs.citrix.com/2012/10/17/announcing-general-availability-of-sharefile-with-storagezones/

  8. https://jhubiostatistics.shinyapps.io/recount/

  9. jx 1 jx 2 jx 3 jx 4 jx 5

    jx 6 Coverage Reads Gene Isoform 1 Isoform 2 Potential isoform 3 exon 1 exon 2 exon 3 exon 4 Expressed region 1: potential exon 5 f1000research.com/articles/6-1558/v1
  10. f1000research.com/articles/6-1558/v1

  11. exon 1 exon 2 exon 3 f1000research.com/articles/6-1558/v1

  12. disjoint exon 1 disjoint exon 2 disjoint exon 3 f1000research.com/articles/6-1558/v1

  13. f1000research.com/articles/6-1558/v1

  14. f1000research.com/articles/6-1558/v1

  15. > library('recount') > download_study( 'ERP001942', type='rse-gene') > load(file.path('ERP001942 ', 'rse_gene.Rdata'))

    > rse <- scale_counts(rse_gene) github.com/leekgroup/recount-analyses/
  16. >library('recount') > download_study('SRP029880', type='rse-gene') > download_study('SRP059039', type='rse-gene') > load(file.path('SRP029880 ',

    'rse_gene.Rdata')) > load(file.path('SRP059039', 'rse_gene.Rdata')) > mdat <- do.call(cbind, dat) github.com/leekgroup/recount-analyses/
  17. Collado Torres et al. Nat. Biotech 2017

  18. slide adapted from Kai Kammers Can combine with genotype data

    to identify eQTLs
  19. biorxiv.org/content/early/2018/01/12/247346

  20. http://www.rna-seqblog.com/ Can we use expression data to predict tissue? slide

    adapted from Shannon Ellis
  21. Number of Regions 589 589 589 589 Number of Samples

    (N) 4,769 4,769 7,193 8,951 97.3% 96.5% 71.9% 50.6% Tissue prediction is accurate across data sets doi.org/10.1093/nar/gky102
  22. Number of Regions 589 589 589 589 589 Number of

    Samples (N) 4,769 4,769 613 6,579 8,951 97.3% 96.5% 91.0% 70.2% Prediction is more accurate in healthy tissue 50.6% doi.org/10.1093/nar/gky102
  23. > library('recount') > download_study( 'ERP001942', type='rse-gene') > load(file.path('ERP001942 ', 'rse_gene.Rdata'))

    > rse <- scale_counts(rse_gene) > rse_with_pred <- add_predictions(rse_gene) github.com/leekgroup/recount-analyses/
  24. Ashkaun Razmara, in prep.

  25. expression data for ~70,000 human samples samples phenotypes ? GTEx

    N=9,962 TCGA N=11,284 SRA N=49,848 samples expression estimates gene exon junctions ERs Answer meaningful questions about human biology and expression sex tissue M Blood F Heart F Liver slide adapted from Shannon Ellis
  26. Code Example: research.libd.org/recount-brain/example_PMI/example_PMI.html research.libd.org/recount-brain/example_PMI/example_PMI.Rmd Replicates part of the GTEx PMI

    paper by Ferreira et al. doi.org/10.1038/s41467-017-02772-x Ashkaun Razmara, in prep.
  27. slide adapted from Jeff Leek

  28. The recount2 team Hopkins Kai Kammers Shannon Ellis Margaret Taub

    Kasper Hansen Jeff Leek Ben Langmead OHSU Abhinav Nellore LIBD Leonardo Collado-Torres Andrew Jaffe recount-brain Ashkaun Razmara Funding and hosting NIH R01 GM105705 NIH 1R21MH109956 CONACyT 351535 AWS in Education Seven Bridges IDIES SciServer
  29. expression data for ~70,000 human samples (Multiple) Postdoc positions available

    to - develop methods to process and analyze data from recount2 - use recount2 to address specific biological questions This project involves the Hansen, Leek, Langmead and Battle labs at JHU Contact: Kasper D. Hansen (khansen@jhsph.edu | www.hansenlab.org) Code for making a 5 min video on recount2: github.com/lcolladotor/biopeerprize2018