Survival Analysis Tutorial

Warmup Questions

If the video was too fast for you, or if you want to dive deeper into the concepts covered, we have also prepared a supplemental notebook. The topic of the Data & Survival Curve video is covered in Chapter 1. Please note that this is a complementary resource. To answer the warmup questions, going through the introductory video should suffice.

Survival Data And Estimation Of The Survival Curve

Consider the following data:

In our graph, the subjects (A, N, F, ...) are in rows. The start of the blue line tells us when we started observing the subject; this could be the time of the intervention, for example. The end of the blue line tells us when we last checked or observed an event. Observed events are marked with an X, all other data are censored. For example, we started observing subject K in week 5, and the event occurred in week 8. The total survival time for K is therefore 3 weeks. The survival time for N is 4 weeks, and at the last visit the event has not yet occurred.

Use the graph above and construct a data table ready for survival analysis. Use Excel or a similar spreadsheet program. Feed the data into Orange and answer the following questions:

By definition, median survival is the time at which half of the participants in a study population are expected to experience the event of interest. In Orange, you can view its value in the Kaplan-Meier widget.

Congratulations, this concludes the quiz for this section!

Before proceeding to the next section, please take a minute to give us some feedback on this introductory section of the tutorial. While answering these questions is optional, your answers will help us improve the software and the tutorial.

Exploring Survival Differences Between Cohorts Of Patients

Here we take a look at the U.S. Veteran's Administration Lung Cancer Trial (from Kalbfleisch D and Prentice RL, 1980), in which male patients with advanced, inoperable lung cancer received either standard therapy or investigational chemotherapy. The data set includes 137 patients, 9 of whom left the study before death. The study was designed to assess the benefit of test chemotherapy and analyze the effects of other covariates.

Your first task is to load the data into Orange and view the survival curve for the entire cohort of patients.

You can easily access this study's data in Orange using the Datasets widget. Search for "Veterans" and the dataset should appear at the top. This data already includes the information on time and event features (time since diagnosis [months], survival event), and you can feed the data to the Kaplan-Meier plot directly from the Datasets widget; you don't need to use the As Survival widget here, although it doesn't hurt.

The next task is to split the data into two cohorts. Our lung cancer trial dataset contains a variable called "Treatment" that indicates whether patients received standard chemotherapy or an alternative, test chemotherapy. We can use this variable to split the patient into two cohorts, the control and the test group. Use the Kaplan-Meier widget to answer the following questions:

While you can calculate the desired summary statistics by combining some of the other widgets in Orange, you can get the answer to the questions here directly with the Kaplan-Meier widget by grouping the data in this widget accordingly. The graph legend reports the counts, where "N" in the legend refers to the size of the cohort, and "n" refers to the number of uncensored (surviving) patients.

Using Different Features for Splitting Patients into Cohorts

Let us analyze the effect of other variables in this data set. How does age affect patient survival? How is survival affected by the Karnofsky performance score?

The Karnofsky Performance Score is a scale used to quantify a patient's general well-being and ability to carry out daily activities. Its typical values are 10 to 30 for fully hospitalized patients, 40 to 60 partially hospitalized patients, and 70 to 90 for patients able to care for themselves. Observe the differences between survival for two cohorts of patients split by the median value for age, and similar for two cohorts based on Karnofsky performance score. You can use the Discretize widget and "Equal frequency" discretization for grouping by the median value.

Congratulations, this concludes the quiz for this section!

Before proceeding to the next section, we again ask you to give us some feedback, this time on this current section on exploring of survival features. Thank you!

Expression Data In Survival Analysis

Let’s look at the cervical cancer data from The Cancer Genome Atlas (TCGA). The dataset is available in the Datasets widget under “TCGA-CESC”. It includes survival data and gene expression values for 306 patients with cervical cancer.

In the literature, genes associated with regulating the UV response have been researched as potential therapeutic targets for cervical cancer treatment (e.g., Gu et al., 2019).

Your task is to identify a handful of potential marker genes associated with survival. For simplicity, you will look only at genes associated with a down-regulated ultraviolet radiation (UV) response. Oncogenes from the human papillomavirus (HPV), the leading cause of cervical cancer, have a complicated relationship with cellular response to UV.

The gene set associated with a down-regulated UV response can be found in the Gene Sets widget under the Molecular Signatures Database (MSigDB) among the Hallmark Gene Sets under the name "HALLMARK_UV_RESPONSE_DN".

Gene Ranking

Do not forget to install the Bioinformatics add-on first, and use Gene Sets widget from this add-on to answer the question.

Now load the TCGA-CESC data and consider only the genes that make up the down-regulated UV response hallmark gene set.

Gene Expression-Based Cohort Analysis

Create two patient groups based on the expression of the highest-ranked gene from the set of downregulated UV response genes. Then, chart the corresponding survival curves.

To answer the following quiz questions, you will have to assemble the longest data analysis pipeline yet, using widgets such as Genes, Gene Sets, Rank Survival Features, Discretize, and Kaplan-Meier. You will become a survival analysis expert!

We again assume that, given a gene, cohorts are created by splitting the patients according to the median gene expression.

Congratulations, this concludes the quiz for this section!

Please give us your feedback on the Gene Ranking chapter of our tutorial. Thank you!

Gene Sets Ranking

In 2019, Xia et al. published a paper identifying a new set of proteins that cause endogenous DNA damage when overproduced. The researchers discovered these DNA damage-up proteins (DDPs) in E. coli and identified 284 human homologs that are overrepresented among known cancer drivers.

They constructed three sets of genes: a set of all 284 candidate DDPs, a set of DDPs excluding known cancer drivers, and a set of DDPs excluding validated DDPs. Known cancer drivers are those whose gain- or loss-of-function in driving cancer has been established in the literature. Validated DDPs, on the other hand, refer to a subset of DDPs identified as actual DNA damage initiators in human cells by Xia et al. They evaluated the association of the three DDP gene sets with overall survival by calculating a gene set enrichment score for each sample and comparing two cohorts, one above and one below the top tertiles. Using any of the three DDP gene sets resulted in significant survival differences between the formed cohorts. These results indicate that there are genes among the discovered human DDP candidates that were previously unknown to drive cancer.

Check whether the three DDP gene sets (all, known excluded, validated excluded) are associated with decreased overall survival in the BRCA dataset, even when you split the patients into cohorts by the median of the enrichment scores (Xia and co-authors used tertiles for splitting, see above).

The BRCA dataset is available in the Datasets widget under the name “TCGA-BRCA”. Loading the dataset might take a minute.

The gene sets are available here. Use the worflow with File, Genes, and Gene Sets widget (see the image below) to load the gene sets into Orange. Once you load custom gene sets into the Gene sets widget, continue using the widget like you usually would. Alternatively, you can download a pre-constructed workflow and open it with Orange.

Remember to pass the TCGA-BRCA data through the As Survival Data widget. You will need to use the Single Sample Scoring widget to score the gene set. After scoring, remove the gene expression information using Select Columns and move the gene set scores to the Features section. You can score each gene set individually using Discretize and Kaplan-Meier combination, or just use the Rank Survival Features widget to calculate the p-values for the differences in the survival curves.

We can use the same procedure to check the utility of other gene sets. Your next task is to rank all fifty gene sets from the Hallmark Gene Set collection with regards to overall survival in the BRCA dataset (order by p-value).

Your workflow from the previous question is ready to answer this question as well; just switch to the Hallmark Gene Sets in the Gene Sets widget. Remember to seected all of the listed gene sets, as the Gene Sets widgets output only the setected sets. You will also need to move the computed scores in the Select Columns widget from the Metas panel to the Features.

The 3rd and 4th ranking gene sets Hallmark Gene Set collection are associated with regulating DNA damage. Form groups based on these two gene sets and plot the corresponding survival curves.

Here, you could stay with the workflow from the previous two questions, select the 3rd and 4th gene sets in the Rank Survival Features widget, and than pass it to Discretize to create an indicator variable for the patient cohorts. Make sure you discretize only the two gene sets, and not the Time variable.

Congratulations, this completes the quiz for this chapter and all the topic quizzes in this tutorial!

Again, we ask you to answer the following questions to help us improve our tutorial in the future. Thank you!