In this tutorial we will analyse a public data set (Georg et al. 2021, https://doi.org/10.1016/j.cell.2021.12.040) of single cell phospho-proteomics measured with Cytometry by Time of Flight (CyTOF).
In this study, the authors performed CyTOF of whole blood samples from mild and severe COVID-19 patients during the acute and convalescent phase, and patients with other acute respiratory infections (Flu-like illness), as well as patients chronically infected by human immunodeficiency virus (HIV) or hepatitis B (HBV) and healthy controls. They analysed the T cell space and identified highly activated CD16+ T cells in severe COVID-19, which led the authors to hypothesise about the pathological role of these cytotoxic T cells. This hypothesis was then tested and confirmed with functional analyses, and found suitable mechanisms for their induction.
In this tutorial you will learn how to perform such computational analysis either in T cells, B cells, or monocytes! We invite you to give a look at the paper before the tutorial.
See you soon!
Please download the following data in advance!
https://box.hu-berlin.de/d/a40f4c6d19eb4be68d7e/
We will use so called .Rmd files which are a combination of R-code and nice formatting. The advantage is that the code and its results are directly connected.
In the last course, you got an introduction to R and its usage, but we took most of the installation-burden from you. After a year of R-usage you must be R-experts now. Therefore you have to install R locally now, making it easier to use for your own datasets.
For Windows, the installation steps are as follows. For Linux or Mac please reach out to us.
- Download R from https://cran.r-project.org/bin/windows/base/
- Install R.
- Download
rtools40-x86_64.exefrom https://cran.r-project.org/bin/windows/Rtools/rtools40.html - Install Rtools
- Download RStudio https://www.rstudio.com/products/rstudio/download/#download
- Install RStudio
There are many packages in R that hold much functionality you can directly use and do not have to code for yourself. We supply you with a list you should install.
For that:
- Open the installed Rstudio
- Create a new file (
File->New File->R Script) - Copy the following code inside that freshly generated script:
if (!require(pacman)) { install.packages("pacman") # If not already installed } pacman::p_load("tidyverse", install = TRUE) pacman::p_load("cowplot", install = TRUE) pacman::p_load("uwot", install = TRUE) pacman::p_load("data.table", install = TRUE) pacman::p_load("pheatmap", install = TRUE) # Todo pacman::p_load("needs", install = TRUE) pacman::p_load("knitr", install = TRUE) pacman::p_load("ggridges", install = TRUE) pacman::p_load("grDevices", install = TRUE) # Install BiocManager such that pacman can then find ComplexHeatmap by its own pacman::p_load("BiocManager", install = TRUE) pacman::p_load("ComplexHeatmap", install = TRUE) # Github packages have to be dealt in a special way # However you install them, you _need_ the "devtools" package pacman::p_load("devtools", install = TRUE) # Additionally, you might need "Rtools" pacman::p_load_gh("JinmiaoChenLab/Rphenograph", install = TRUE) pacman::p_load_gh("JinmiaoChenLab/cytofkit", install = TRUE)
- Either run each line of the code (
Ctrl+Enter per line) or all of the code at once (Ctrl+Shift+s) - There will be some red messages, but there should be no warnings nor errors.
Rstudio has a nice feature to organize projects.
- On the top right corner there is a button stating
Project: (None) - Click on that button and select
New Project - --> Select
New Directory - --> Select
New Project - --> Type in
InstructTutorial_IIas Directory name and select a folder where you want to place it inside.InstructTutorial_IImust not exist in the selected folder already.
After you set up this project, on the bottom right, find the button More and click Show Folder in New Window.
This opens the explorer in the freshly created directory. In this directory:
- Create a new directory called
codewhere you will later place all of the code. - Create a new directory called
data - Copy the
data_norm_sub5.csvfrom 2. (https://box.hu-berlin.de/d/a40f4c6d19eb4be68d7e/) into the newdatadirectory. - In 3.2 you generated a new code file but did not save it anywhere. Create the file again, now inside the project and save it in the
codedirectory under the filename01_setup.R
You can get all information, the code under https://github.com/gunthergl/instructtutorial as soon as we open it publicly. You do not necessarily go through all of that by yourself. Therefore we will open it rather late before the actual course.
Just let us know if you need help:
- Gunther Glehr (
prename.surname@ukr.de) - Rosario Astaburuaga Garcia (
prename.midname@charite.de) - Lev Petrov (
prename.surname@charite.de)