Learn how to use the pick samples feature, how to view multiple genes in a single column, how to view a signature, and how to run a differential expression analysis
Description
This tutorial is made for those who have basic knowledge of how to use Xena. We will cover how to use the pick samples feature, how to view multiple genes in a single column, how to enter and view a signature, and how to run a differential expression analysis.
This tutorial assumes basic knowledge of how to build and read a Visual Spreadsheet. To get this, go through Basic Tutorial: Section 1. It also assumes basic knowledge of filtering. To get this, go through Basic Tutorial: Section 2.
Part A: 10 min
Part B: 5 min
Part C: 15 min
Part A
Create a visual spreadsheet with single column with multiple genes.
Filter to only Primary Tumor samples using the Pick Samples mode.
Remove nulls using the option in the filter menu
Part B
Enter and view a gene expression signature
Part C
Run a differential expression analysis.
We will investigate the PAM50 molecular subtypes in breast cancer. PAM50 is a 50-gene signature that classifies breast cancer into five molecular intrinsic subtypes: Luminal A, Luminal B, HER2-enriched, Basal-like, and Normal-like.
We will make a visual spreadsheet where we can explore the relationship between the PAM50 subtype call and the 50 genes that make up the PAM50 subtype call.
Start at https://xenabrowser.net/
Type 'TCGA Breast Cancer (BRCA)', select this study from the drop down menu, and click 'To first variable'.
Choose 'Phenotypic', select 'sample_type' from the dropdown menu, and click 'To second variable'.
Choose 'Phenotypic', click on 'advanced', type 'pam' into the search bar, select 'PAM50Call_RNAseq' from the dropdown menu, and click 'Done'. This will exit the wizard.
Click on 'Click to insert a column' after column C. Copy and paste the 50 genes, choose 'Gene Expression', and click 'Done'.
Click the handle in the lower right corner of column D, mutation. Move it to the right to make the column bigger.
List of 50 genes used to calculate the PAM50 subtype call:
UBE2T BIRC5 NUF2 CDC6 CCNB1 TYMS MYBL2 CEP55 MELK NDC80 RRM2 UBE2C CENPF PTTG1 EXO1 ORC6L ANLN CCNE1 CDC20 MKI67 KIF2C ACTR3B MYC EGFR KRT5 PHGDH CDH3 MIA KRT17 FOXC1 SFRP1 KRT14 ESR1 SLC39A6 BAG1 MAPT PGR CXXC5 MLPH BCL2 MDM2 NAT1 FOXA1 BLVRA MMP11 GPR160 FGFR4 GRB7 TMEM45B ERBB2
Click on the picker icon next to the filter menu to enter pick samples mode.
Click on the Primary Tumor samples.
Click the filter menu and select 'Keep samples'.
Exit pick samples mode by clicking on the picker icon again.
Click the filter menu and select 'Remove samples with nulls'.
More information:
We will now look at the TFAC30 gene signature and see how it relates to the PAM50 subtype calls. This gene expression signature over 30 genes predicts pathologic complete response (pCR) to preoperative weekly paclitaxel and fluorouracil-doxorubicin-cyclophosphamide (T/FAC) chemotherapy.
Click on 'Click to insert a column' after column D. Copy and paste the signature below, choose 'Gene Expression', and click 'Done'. Note you need to include the '=' as this tells Xena that you want the signature rather than to see all the genes individually.
TFAC30 gene expression signature:
=E2F3 + MELK + RRM2 + BTG3 - CTNND2 - GAMT - METRN - ERBB4 - ZNF552 - CA12 - KDM4B - NKAIN1 - SCUBE2 - KIAA1467 - MAPT - FLJ10916 - BECN1 - RAMP1 - GFRA1 - IGFBP4 - FGFR1OP - MDM2 - KIF3A - AMFR - MED13L - BBS4
We can now see that patient's samples that are labeled as 'Her2' and 'Basal' are predicted to be more likely to achieve pCR on TFAC chemotherapy.
More information
We will run a differential expression analysis comparing Basal samples to Luminal A and Luminal B samples.
Click the column menu for the PAM50 subtype call (column C) and chose 'Differential Expression'. This will open a new tab where we will run the analysis.
Choose the first subgroup to be 'Basal' and the second subgroup to be 'LumA' and 'LumB'. Hold the shift key while clicking to select multiple groups.
Click 'Submit'.
Note it can take a while for the analysis to run. Wait until it says 'Success' at the top.
More information