Mick at 11:48, 7 January 2026

2026-01-07T11:48:34Z

← Older revision		Revision as of 13:48, 7 January 2026
Line 82:		Line 82:
	<p>To start an R session in your terminal typing (or copy/pasting)</p>		<p>To start an R session in your terminal typing (or copy/pasting)</p>
	<pre>		<pre>
	R-4.2.2		/home/ctools/opt/R-4.4.2_22140/bin/R
	</pre>		</pre>
	<p>And load the library we need by typing</p>		<p>And load the library we need by typing</p>

WikiSysop: Created page with "

Overview and background

Groups

Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.

Assignment notes

While some question..."

2024-03-19T15:37:12Z

Created page with " <div class="page-content has-page-title"> <div id="overview-and-background" class="section level1"> <h1>Overview and background</h1> <div id="groups" class="section level2"> <h2>Groups</h2> <p>Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.</p> </div> <div id="assignment-notes" class="section level2"> <h2>Assignment notes</h2> <p>While some question..."

New page

<div class="page-content has-page-title">
<div id="overview-and-background" class="section level1">
<h1>Overview and background</h1>
<div id="groups" class="section level2">
<h2>Groups</h2>
<p>Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.</p>
</div>

<div id="assignment-notes" class="section level2">
<h2>Assignment notes</h2>
<p>While some questions might seem hard we naturally don’t ask questions/tasks which you have not been given the tools to solve in this assignment - so if you are stuck try thinking about what you have already learned before asking an instructor.</p>
</div>

<div id="assignment-overview" class="section level2">
<h2>Assignment overview</h2>
<p>In this assignment you are going to analyze RNA-sequencing data from real cancer patients to analyze the importance of alternative splicing in a clinical context</p>
</div>

<div id="biological-background" class="section level2">
<h2>Biological background</h2>
<p>Today you will be working with colorectal cancers - specifically Colon Adenocarcinoma (often abbreviated COAD). It is a cancer of the colon that is very frequent. The lifetime risk of developing
colorectal cancer is ~4% for both males and females. That means COAD represents ~10% of all cancers and results in the death of hundreds of thousands of people each year! (More info on COAD can be found on [https://en.wikipedia.org/wiki/Colorectal_cancer Wikipedia].</p>

<p>One important aspect of cancer is that tumors from different patients are extremely different even when they originate from the same tissue (more info on tumor heterogeneity [https://en.wikipedia.org/wiki/Tumour_heterogeneity here]). To improve treatment and prognosis we therefore try to classify COAD into cancer subtypes (a simple form of precision medicine). We currently think there are 5 subtypes (see [https://www.cell.com/cancer-cell/pdf/S1535-6108(18)30114-4.pdf Liu ''et al.'']) and today you will be working with CIN and GS. CIN is an abbreviation for Chromosomal INstable and GS means genome stable. More on that later.</p>

<p>To help us understand COAD subtypes you will today compare these to healthy adjacent tissue. For all samples a biopsy was taken and bulk RNA-seq performed. Low-quality samples have been removed.</p>

</div>
<div id="bioinformatic-background" class="section level2">
<h2>Bioinformatic background</h2>
<p>For background on transcriptomics and splicing please refer to today’s slides. The data you are working with is a randomly selected a subset of the TCGA COAD data (google TCGA if you want to know more). The data was quantified with Kallisto against the human transcriptome.</p>

<p>Today you will be using the 'pairedGSEA' R package we developed. This package is specifically designed to make it easy to do the following analysis:</p>

<ol style="list-style-type: decimal">
<li>Differential gene expression (aka DGE) via DESeq(2)</li>
<li>Differential gene usage (differential splicing) (aka DGU)</li>
<li>gene-set over-representation analysis (ORA) on DGU and DGE
results</li>
</ol>
<p>While at each step facilitating easy comparison of DGE and DGU.</p>
<hr />
</div>
</div>
<div id="assignment" class="section level1">
<h1>Assignment</h1>
<div id="step-1-determine-which-cancer-to-work-with" class="section level2">
<h2>Step 1: Determine which cancer to work with</h2>
<p>Determine which cancer type you will work with:</p>
<ul>
<li>If your birthday is within the first 6 months of the year (January-June) you will work with <strong>CIN</strong>.</li>
<li>If your birthday is within the last 6 months of the year (July-December) you will work with <strong>GS</strong>.</li>
</ul>
</div>
<div id="step-2-set-up-enviroment" class="section level2">
<h2>Step 2: Set up enviroment</h2>
<p>Log into the server as you usually do except this time you have to use the '-X' option. That means using:

<pre>
ssh -X username@pupil1.healthtech.dtu.dk</pre>.
</p>

<p>Make a directory for this exercise and move into it</p>
<pre>
mkdir transcriptomics_exercise
cd transcriptomics_exercise
</pre>

<p>Copy the exercise data of your cancer subtype to your folder</p>
<pre>
### for CIN subtype:
cp /home/projects/22126_NGS/exercises/transcriptomics/coad_iso_subset_cin.Rdata .

### For GS subtype:
cp /home/projects/22126_NGS/exercises/transcriptomics/coad_iso_subset_gs.Rdata .
</pre>

</div>
<div id="step-3-start-r-session-and-enviroment" class="section level2">
<h2>Step 3: Start R session and enviroment</h2>
<p>To start an R session in your terminal typing (or copy/pasting)</p>
<pre>
R-4.2.2
</pre>
<p>And load the library we need by typing</p>
<pre>
library(pairedGSEA)
</pre>

<p>This loads the functionality of the “pairedGSEA” R package.</p>
</div>
<div id="step-4-load-and-inspect-data" class="section level2">
<h2>Step 4: Load and inspect data</h2>
<p>Load the assignment data into your R session:</p>
<pre>
### for CIN subtype:
load('coad_iso_subset_cin.Rdata')

### For GS subtype:
load('coad_iso_subset_gs.Rdata')
</pre>
<p>This will give you two data objects in your R session:</p>
<ol style="list-style-type: decimal">
<li>A count matrix</li>
<li>A matrix with meta information about each sample in the count matrix.</li>
<li>A list of gene_sets that you should use for your ORA analysis (step 7).</li>
</ol>

<p>All objects can be directly used by the 'pairedGSEA'
package - no need to do any data modifications.</p>
<p><br></p>
<p>Use the following functions to take a look at the data:</p>
<pre>
### List objects in an R session
ls()

### Inspect the first lines of the object
head( <object_name> )
</pre>

<p><strong>Question</strong>: Which object contains what data?</p>

</div>
<div id="step-5-run-differential-analysis" class="section level2">
<h2>Step 5: Run differential analysis</h2>
<p>Next you will need to use the 'pairedGSEA' package and
here a bit of self-study is needed. <strong>Importantly</strong> you
should only run this analysis once per group - else we don’t have
enough computational power. You can download the
'pairedGSEA' vignette (short document showing how to use it)
<a href="https://www.dropbox.com/s/oalth29pxulffec/pairedGSEA.html?dl=1">here</a>.</p>
<p>Hints:</p>
<ol style="list-style-type: decimal">
<li>After reading the introduction you can skip to the
'3.3 Running the analysis' section.</li>
<li>For now you only need to use 'paired_diff()' as that
makes both differential analyses (both DGE and DGU).</li>
<li>There is no need to use the “store_results” option</li>
</ol>
<p><strong>Question</strong>: This will take a while to run (~10 min).
In the mean time take a closer look at the Liu <em>et al.</em> paper
(see above) and summarise what the difference between the CIN and GS
COAD subtypes are.</p>

</div>
<div id="step-6-inspect-diffrential-result" class="section level2">
<h2>Step 6: Inspect diffrential result</h2>
<p><strong>Question</strong>: Look at the first 10 lines of the result
file. Which gene is most significant (smallest p-value) for the DGE and
DGU analysis (respectively DESeq2 and DEXSeq)</p>

<p>The following code <em>example</em> counts how many significantly
differentially expressed genes are found:</p>
<pre>
sum( gi_diff_results$padj_deseq < 0.05, na.rm = T )
</pre>
<p><strong>Question</strong>: Modify the R code above to count how many
genes are DGE and DGU.</p>

<p><strong>Question</strong>: Use the 'nrow()' function to
calculate the fraction of genes that are DGE and DGU.</p>

<p>Now we are ready to do the gene-set enrichment analysis.</p>
</div>
<div id="step-7-run-gene-set-enrichment-analysis" class="section level2">
<h2>Step 7: Run Gene-Set Enrichment Analysis</h2>
<p>Use the vignette to help you use 'pairedGSEA' to run GSEA on both DGE and DGU results (see the vignette section 4: “Over-Representation Analysis”). You should use the 'gene_set_list' object you have already loaded into R instead of using the 'prepare_msigdb()' function.</p>

<p>Note: There is (again) no need to store the intermediary results.</p>

</div>
<div id="step-8-inspect-ora-result" class="section level2">
<h2>Step 8: Inspect ORA result</h2>
<p>What you have been analyzing so far is a subset of the entire dataset
(since the runtime else would have been 3-4x longer). To enable a more
realistic last step use <strong>one</strong> of these commands to load
the full results corresponding to what you have been working with.</p>
<pre>
### for CIN subtype:
load('/home/projects/22126_NGS/exercises/transcriptomics/03_coad_cin_ora.Rdata')
# loads the "cin_ora" object

### For GS subtype:
load('/home/projects/22126_NGS/exercises/transcriptomics/03_coad_gs_ora.Rdata')
# loads the gs_ora object
</pre>
<p>The following code <em>example</em> extract the ORA analysis of
either DGU and DGE and sorts it so the most significant gene-sets are at
the top.</p>

<pre>
### DGE:
dge_ora_sorted <- gi_paired_ora[
sort.list(gi_paired_ora$pval_deseq), # sort part
c('pathway','pval_deseq','enrichment_score_deseq') # select part
]

### DGU ORA:
dgu_ora_sorted <- gi_paired_ora[
sort.list(gi_paired_ora$pval_dexseq), # sort part
c('pathway','pval_dexseq','enrichment_score_dexseq') # select part
]
</pre>

<p><strong>Question</strong>: Look at the 10-15 most significant gene
sets from both analyses. What are the similarities and differences?</p>
</div>
<div id="step-9-visual-inspection-of-ora-result" class="section level2">
<h2>Step 9: Visual inspection of ORA result</h2>
<p><strong>Question</strong>: Based on your insights from step 8 use the 'plot_ora()' functionality to test if these are just examples or generalize to all the significant results. An example: If I from the 10-15 top gene-sets saw that only DGU had gene-sets covering “telomer” function I would use the 'plot_ora()' function to test this.</p>
<p><strong>Question</strong>: Try to make a hypothesis as to why this/these molecular functions might be important for cancer.</p>

</div>
<div id="step-10-critical-self-evaluation" class="section level2">

<h2>Step 10: Critical self evaluation</h2>
<p><strong>Question</strong>: Take a moment to think about what potential problems there could be with this assignment. Are there any obvious things we have not taken into consideration?</p>

</div>
<div id="step-11-repport-result" class="section level2">

<h2>Step 11: Report result</h2>
<p>Go to the blackboard and report one or more of the following:</p>
<ul>
<li>A keyword that showed a similar enrichment pattern in DGU and DGE</li>
<li>A keyword that showed preferential regulation through DGU or DGE</li>
</ul>
<hr/>
</div>
</div>
<div id="bonus-assignment" class="section level1">
<h1>Bonus Assignment</h1>
<p>Use 'pairedGSEA' to analyze the other COAD cancer subtype (the one you did not analyze). Are the gene-sets similar or different between the subtypes and analysis types?</p>
</div>
</div>

Rnaseq exercise - Revision history

Mick at 11:48, 7 January 2026

WikiSysop: Created page with " Overview and background Groups Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group. Assignment notes While some question..."

Overview and background

Groups

Assignment notes

WikiSysop: Created page with "

Overview and background

Groups

Please get into groups of 2-3. We don’t have enough computational power for all of you working alone. Please let the instructors know if you need help finding a group.

Assignment notes

While some question..."