Skip to content

Commit

Permalink
Restructure some figures for better pdf output
Browse files Browse the repository at this point in the history
  • Loading branch information
royfrancis committed Mar 21, 2024
1 parent 5a41d94 commit e99724e
Show file tree
Hide file tree
Showing 2 changed files with 59 additions and 27 deletions.
6 changes: 6 additions & 0 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -47,6 +47,12 @@ docker run --platform linux/amd64 --rm -u 1000:1000 -v ${PWD}:/qmd ghcr.io/nbisw

:exclamation: Output files are for local preview only. Do not push any rendered .html files or intermediates.

## Convert HTML slides to PDF

```
docker run --platform=linux/amd64 -v $PWD:/work astefanutti/decktape https://nbisweden.github.io/workshop-ngsintro/2403/topics/rnaseq/slide_rnaseq.html /work/slide_rnaseq.pdf
```

## Repo organisation

The source material is located on the *master* branch (default). The rendered material is located on the *gh-pages* branch. One only needs to update source materials in *master*. Changes pushed to the *master* branch is automatically rendered to the *gh-pages* branch using github actions.
Expand Down
80 changes: 53 additions & 27 deletions topics/rnaseq/slide_rnaseq.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -51,12 +51,12 @@ upid <- yaml::read_yaml(here("_quarto.yml"))$uppmax_project
## Workflow

::: {.columns}
::: {.column width="70%"}
::: {.column width="65%"}

![](assets/rnaseq_workflow.svg)

:::
::: {.column width="30%"}
::: {.column width="35%"}

![](assets/sequence.jpg)

Expand Down Expand Up @@ -205,28 +205,28 @@ Poor quality
## Read QC • PBSQ, PSQS

**Per base sequence quality**
![](assets/pbsq.jpg){width="60%"}
![](assets/pbsq.jpg){width="50%"}

**Per sequence quality scores**
![](assets/psqs.jpg){width="60%"}
![](assets/psqs.jpg){width="50%"}

::: {.notes}

## Read QC • PBSC, PSGC

**Per base sequence content**
![](assets/pbsc.jpg){width="60%"}
![](assets/pbsc.jpg){width="50%"}

**Per sequence GC content**
![](assets/psgc.jpg){width="60%"}
![](assets/psgc.jpg){width="50%"}

## Read QC • SDL, AC

**Sequence duplication level**
![](assets/sdl.jpg){width="60%"}
![](assets/sdl.jpg){width="50%"}

**Adapter content**
![](assets/ac.jpg){width="60%"}
![](assets/ac.jpg){width="50%"}

:::

Expand Down Expand Up @@ -303,7 +303,7 @@ Poor quality

:::

## Mapping
##

- Reads (FASTQ)

Expand All @@ -321,7 +321,6 @@ NCATCGTGGTATTTGCACATCTTTTCTTATCAAATAAAAAGTTTAACCTACTCAGTTATGCGCATACGTTTTTTGATGGC
```
>1 dna:chromosome chromosome:GRCz10:1:1:58871917:1 REF
GATCTTAAACATTTATTCCCCCTGCAAACATTTTCAATCATTACATTGTCATTTCCCCTC
CAAATTAAATTTAGCCAGAGGCGCACAACATACGACCTCTAAAAAAGGTGCTGTAACATG
```

- Annotation (GTF/GFF)
Expand Down Expand Up @@ -557,13 +556,16 @@ ENSG00000000460 33 27 23 19 27 42 69 44 40

![](assets/correlation.png){width="50%"}

- Pairwise correlation between samples must be high (>0.9)
Pairwise correlation between samples must be high (>0.9)

## MultiQC

[![](assets/multiqc.png)](https://multiqc.info/examples/rna-seq/multiqc_report.html)
[![](assets/multiqc.png){width="90%"}](https://multiqc.info/examples/rna-seq/multiqc_report.html)

## Normalisation
## Normalization

::: {.columns}
::: {.column width="70%"}

- Control for Sequencing depth, compositional bias and more
- Median of Ratios (DESeq2) and TMM (edgeR) perform the best
Expand All @@ -575,10 +577,16 @@ ENSG00000000460 33 27 23 19 27 42 69 44 40
- For own analysis, plots etc, use TPM
- Other solutions: spike-ins/house-keeping genes

<img src="assets/distribution.jpg" style="height:400px; position:fixed; right:0px; bottom: 0px; margin-right: 70px; margin-bottom: 140px;">

@dillies2013comprehensive, @evans2018selecting, @wagner2012measurement

:::
::: {.column width="30%"}

![](assets/distribution.jpg)

:::
:::

## Exploratory

```{r,echo=FALSE,eval=FALSE}
Expand Down Expand Up @@ -677,12 +685,15 @@ add_markers(p1, color = ~patient)
## Differential expression

:::: {.columns}
::: {.column width="70%"}
::: {.column width="65%"}

![](assets/dge.png){width="90%"}
![](assets/dge.png)

- Univariate testing gene-by-gene
- More descriptive, less predictive

:::
::: {.column width="30%"}
::: {.column width="35%"}

::: {.fragment}

Expand All @@ -704,15 +715,12 @@ ggplot(dfr,aes(gene,value,col=group))+
theme(legend.position="top")
```

![](assets/dge-scatter.png){width="90%"}
![](assets/dge-scatter.png){width="80%"}

:::
:::
:::

- Univariate testing gene-by-gene
- More descriptive, less predictive

::: {.notes}

## Differential expression
Expand Down Expand Up @@ -807,7 +815,7 @@ low counts [2] : 0, 0%
(mean count < 3)
```

## DGE
##

:::: {.columns}
::: {.column}
Expand Down Expand Up @@ -860,15 +868,24 @@ plotCounts(d1, "ENSG00000000003", intgroup = "type")
:::
::::

## Functional analysis • GO
## Functional analysis • Gene Ontology

::: {.columns}
::: {.column width="70%"}

- Gene set analysis (GSA)
- Gene set enrichment analysis (GSEA)
- Gene ontology / Reactome databases

![](assets/go.jpg){width="50%"}
![](assets/go.jpg){width="70%"}

<img src="assets/systembio.png" style="height:350px; position:fixed; right:0px; bottom:0px; margin-right: 70px; margin-bottom: 300px;">
:::
::: {.column width="30%"}

![](assets/systembio.png)

:::
:::

## Functional analysis • Kegg

Expand All @@ -892,12 +909,21 @@ plotCounts(d1, "ENSG00000000003", intgroup = "type")
## Further learning

::: {.columns}
::: {.column width="75%"}

- Griffith lab [RNA-Seq using HiSat & StringTie tutorial](https://github.com/griffithlab/rnaseq_tutorial/wiki)
- HBC Training [DGE using DeSeq2 tutorial](https://hbctraining.github.io/main/)
- [RNA-Seq Blog](https://www.rna-seqblog.com/)
- SciLifeLab [courses](https://www.scilifelab.se/education/courses%26training)

<img src="assets/help.png" style="height:250px; position:fixed; right:0px; bottom:0px; margin-right: 35px; margin-bottom: 40px;">
:::
::: {.column width="25%"}

![](assets/help.png)

:::
:::

## {background-image="/assets/images/cover.jpg"}

Expand Down

0 comments on commit e99724e

Please sign in to comment.