Transcriptional pattern enriched for synaptic signaling is associated with shorter survival of patients with high-grade serous ovarian cancer
- Department of Medical Oncology, University Medical Center Groningen, University of Groningen, Netherlands
- Cancer Biology and Immunotherapies Group, Sanford Research, United States
- Penn Ovarian Cancer Research Center and Basser Center for BRCA, University of Pennsylvania, Perelman School of Medicine, United States
- Department of Obstetrics and Gynecology, University Medical Center Groningen, University of Groningen, Netherlands
Figures
Figure 1
Workflow indicates the data acquisition and relations between the methods.
Figure 2 with 3 supplements
Enrichment heatmap of hallmark gene sets in transcriptional components associated with patient overall survival.
Gene Set Enrichment Analysis for 14 transcriptional components (TCs) associated with overall survival (OS) identified through univariate or multivariate survival analyses are presented. Only Hallmark gene sets with significant enrichment (Bonferroni-corrected p-value) for at least one TC are shown. The heatmap displays Z-scores, which indicate the relative enrichment strength, with values truncated at a maximum of 4 for visualization purposes. The gene sets were clustered based on Pearson correlation using the Ward D2 method, providing insights into related biological processes captured by different TCs. In the right column, chromosomal locations of copy number alterations (CNAs) are shown, reflecting the downstream effects on gene expression that each TC captures. This integration of CNA information highlights the biological relevance of each TC and its contribution to gene expression variability and patient outcomes.
Figure 2—figure supplement 1
Association of OS-related TCs with clinicopathologic parameters.
The association between clinicopathologic parameters and the activity of the OS-associated 14 TCs was determined in the complete set of 1125 samples. Pearson correlation was used to calculate the association of each clinicopathologic parameter with each TC. The TCs were then ranked based on their association with OS. Abbreviations: TC = transcriptional component, OS = overall survival.
Figure 2—figure supplement 2
Enrichment heatmap for the KEGG gene set collection in OS-related TCs.
Gene set enrichment analysis (GSEA) results of 14 TCs associated with OS in univariate or multivariate survival analysis are presented, including KEGG gene sets that were included in the enrichment for at least one TC that passed the Bonferroni threshold for multiple testing correction. The gene sets were clustered using Pearson correlation and Ward D2, and the heatmap colors were based on Z-scores, truncated at a value of four. The right column shows the chromosomal location of a copy number alteration that the TC captures the downstream effects on gene expression levels. Abbreviations: TC = transcriptional component, OS = overall survival.
Figure 2—figure supplement 3
Enrichment heatmap for the REACTOME gene set collection in overall survival (OS)-related TCs.
Gene set enrichment analysis (GSEA) results of 14 TCs associated with OS in univariate or multivariate survival analysis are presented, including REACTOME gene sets that were included if the enrichment for at least one TC passed the Bonferroni threshold for multiple testing correction. The gene sets were clustered using Pearson correlation and Ward D2, and the heatmap colors were based on Z-scores, truncated at a value of four. The right column shows the chromosomal location of a copy number alteration that the TC captures the downstream effects on gene expression levels. Abbreviations: TC = transcriptional component.
Figure 3 with 3 supplements
Survival tree analysis of patients with platinum-treated HGSOC defines survival cohorts with distinct clinicopathologic and biological characteristics.
The results of the survival tree analysis of 294 patients with high-grade serous ovarian cancer (HGSOC) treated with platinum-based chemotherapy are presented. The analysis utilized 14 transcriptional components (TCs) associated with overall survival (OS), along with other clinicopathologic factors, including age, tumor stage, grade, and debulking status. The resulting tree identified nine distinct survival cohorts, each represented as a bar in the Sankey diagram, where the bar height corresponds to the number of patients in each cohort. Kaplan-Meier survival curves with accompanying number-at-risk tables are shown for each cohort, with survival data censored at 10 y. The names of the survival cohorts were based on enriched biological processes in the TCs, as determined by the chromosomal location of genes captured by a TC, GSEA, and co-functionality analysis of the top genes. The p-values in the Kaplan-Meier plots were derived from log-rank tests comparing survival distributions between groups. Abbreviations: TC = transcriptional component, ECM = extracellular matrix.
Figure 3—figure supplement 1
Survival tree analysis of patients with advanced-stage, high-grade serous ovarian cancer (HGSOC) defines survival cohorts with distinct clinicopathologic and biological characteristics.
Survival tree analysis of 265 patients with advanced-stage, platinum-treated HGSOC using 14 OS-associated TCs and other classifiers such as age, tumor stage, grade, and debulking status. The analysis resulted in nine survival cohorts, and the height of the bar in the Sankey diagram represents the number of patients in each cohort. The Kaplan-Meier plots and number-at-risk tables are presented with survival data censored at 10 y. The names of the survival cohorts were based on enriched biological processes in the TCs, as determined by the chromosomal location of genes captured by a TC, gene set enrichment analysis (GSEA), and co-functionality analysis of the top genes. The p-values in each panel show the p-value from the corresponding log-rank test between the two groups. Abbreviations: TC = transcriptional component, ECM = extracellular matrix.
Figure 3—figure supplement 2
The activity of TC121 in bulk transcriptomes of patients with different subtypes of ovarian cancer.
The boxplots display the activity scores of TC121 in different cancer subtypes, which are ordered based on their corresponding median activity scores.
Figure 3—figure supplement 3
Three overall survival (OS)-associated TCs capture the transcriptional effect of copy number alterations.
For each OS-related TC, the weight of each gene was plotted on its genomic location. Abbreviations: TC = transcriptional component.
Figure 4
Co-functionality network of top 100 absolute weighted genes in TC121.
Co-functionality network for the top 100 genes with the highest absolute weights in TC121 is presented. Genes were clustered based on predicted co-functionality (r>0.7) across datasets, with clusters identified using both Gene Ontology (GO) Biological Processes and Cellular Components databases. One primary cluster, containing more than five genes, exhibited strong enrichment for synaptic signaling in the GO Biological Processes database and for synaptic membranes in the GO Cellular Components database. This highlights the biological specificity of TC121 in regulating gene expression linked to synaptic functions.
Figure 5 with 3 supplements
The activity of TC121 in bulk transcriptomes of Cancer Cell Line Encyclopedia (CCLE), Genomics of Drug Sensitivity in Cancer (GDSC) cell lines, and Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) patient-derived samples.
(A) Cross-study TC projection of TC121 on CCLE and GDSC cell lines. The boxplots display the activity scores of TC121 in different tissue types, which are ordered based on their corresponding median activity scores. (B) Cross-study TC projection of TC121 on GEO and TCGA bulk transcriptomes resulted in the activity scores presented in the boxplots. Cancer types were ordered based on corresponding medians of TC121 activity scores. Abbreviations: TC = transcriptional component.
Figure 5—figure supplement 1
Activity of overall survival (OS)-associated TCs in bulk transcriptomes of Cancer Cell Line Encyclopedia (CCLE) cell lines.
The proportion of samples with outlier activity scores in individual tissue types for each of the 14 OS-associated TCs was obtained using an absolute cut-off value of 0.05. TCs were ordered based on their association with OS. Abbreviations: TC = transcriptional component.
Figure 5—figure supplement 2
The activity of overall survival (OS)-associated TCs in bulk transcriptomes of Genomics of Drug Sensitivity in Cancer (GDSC) cell lines.
The proportion of samples with outlier activity scores in individual tissue types for each of the 14 OS-associated TCs was obtained using an absolute cut-off value of 0.05. TCs were ordered based on their association with OS. Abbreviations: TC = transcriptional component.
Figure 5—figure supplement 3
Association heatmap for The Cancer Genome Atlas (TCGA) cluster sets with the activity scores of OS-related TCs.
This heatmap highlights the associations between each cluster set and the TCs, represented by log-transformed p-values from the corresponding Kruskal-Wallis test. A significant association between a TC and a cluster set indicates that at least one cluster within the cluster set exhibited significantly different activity scores for the corresponding TC compared to the other clusters. The heatmap colors were based on log-transformed p-values. Abbreviations: TC = transcriptional component. OS = overall survival.
Figure 6 with 2 supplements
Spatial transcriptomic profiles in ovarian cancer samples.
(A) We employed a permutation-based approach to pinpoint the areas of significant TC activity in spatial transcriptomic profiles. We ran 5000 permutations for each TC-profile combination, yielding a p-value that indicates the extent to which the TC activity in the corresponding profile differs from what would be expected by chance (the null distribution). We then transformed these p-values into logarithmic values and represented them using a heatmap. Heatmaps of activity scores of the TCs are presented in individual rows for the high-grade serous ovarian cancer (HGSOC), serous papillary, and endometrioid adenocarcinoma of ovary samples. The first column represents the stained images of the samples. The second to seventh columns show heatmaps corresponding to the mentioned TCs. (B) The heatmap illustrates the colocalization between two TC activities on spatial transcriptomic profiles from ovarian cancer samples. For each cell, the colocalization scores of the TCs at each of the three spatial transcriptomics samples OC 1, OC 2, and OC 3 are arranged in columns. A colocalization score of 4 between two TCs (red) indicates that the positively (+) and negatively (-) active regions of both TCs are perfectly colocalized. Conversely, a colocalization score of –4 between two TCs (blue) also indicates colocalization. Still, with inverse activity, i.e., the positively active regions of one TC are colocalized with the negatively active regions of the other TC or vice versa. A colocalization score close to 0 between two TCs (white) indicates that the activities of two TCs are spatially separated. The dashed and solid circles in the panel on the right side of the color bar represent two different TCs. Abbreviations: TC = transcriptional component.
Figure 6—figure supplement 1
Spatial transcriptomic profiles in eight ovarian cancer samples.
(A) We employed a permutation-based approach to pinpoint the areas of significant TC activity in spatial transcriptomic profiles. We ran 5000 permutations for each TC-profile combination, yielding a p-value that indicates the extent to which the TC activity in the corresponding profile differs from what would be expected by chance (the null distribution). We then transformed these p-values into logarithmic values and represented them using a heatmap. Heatmaps of activity scores of the TCs are presented in individual rows for the high-grade serous ovarian cancer (HGSOC) samples. The first column represents the stained images of the samples. The second to seventh columns show heatmaps corresponding to the mentioned TCs. (B) The heatmap illustrates the colocalization between two TC activities on spatial transcriptomic profiles from ovarian cancer samples. For each cell, the colocalization scores of the TCs at each of the eight spatial transcriptomics samples are arranged in columns. A colocalization score of 4 between two TCs (red) indicates that the positively (+) and negatively (-) active regions of both TCs are perfectly colocalized. Conversely, a colocalization score of –4 between two TCs (blue) also indicates colocalization. Still, with inverse activity, i.e., the positively active regions of one TC are colocalized with the negatively active regions of the other TC or vice versa. A colocalization score close to 0 between two TCs (white) indicates that the activities of two TCs are spatially separated. The dashed and solid circles in the panel on the right side of the color bar represent two different TCs. Abbreviations: TC = transcriptional component.
Figure 6—figure supplement 2
The activity scores of OS-associated TCs in different single-cell types from patients with HGSOC.
Abbreviations: TC = transcriptional component. OS = overall survival. HGSOC = high-grade serous ovarian cancer.
Additional files
-
Supplementary file 1
Patient characteristics.
Abbreviations: NA = not applicable, HGSOC = high-grade serous ovarian cancer, LGSOC = low-grade serous ovarian cancer, LMP = low malignant potential, OCCC = ovarian clear cell cancer.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp1-v1.docx
-
Supplementary file 2
Overview of the number of samples from each Gene Expression Omnibus (GEO) series and the corresponding study.
Abbreviations: PMID = PubMed identification number.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp2-v1.docx
-
Supplementary file 3
Multivariate permutation framework results from the univariate survival analysis.
Abbreviations: TC = transcriptional component. CI = confidence interval.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp3-v1.xlsx
-
Supplementary file 4
Permutation framework results from the multivariate survival analysis.
Abbreviations: TC = transcriptional component. CI = confidence interval.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp4-v1.xlsx
-
Supplementary file 5
Clinicopathological parameters per survival tree node in platinum-treated HGSOC.
Abbreviations: HGSOC = high-grade serous ovarian cancer, TC = transcriptional component, ECM = extracellular matrix.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp5-v1.docx
-
Supplementary file 6
Clinicopathological parameters per survival tree node in advanced-stage platinum-treated HGSOC.
Abbreviations: HGSOC = high-grade serous ovarian cancer, TC = transcriptional component, ECM = extracellular matrix.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp6-v1.docx
-
Supplementary file 7
GenetICA gene network analysis results.
Abbreviations: TC = transcriptional component.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp7-v1.xlsx
-
Supplementary file 8
Results from the univariate survival analysis on patients with ovarian clear cell carcinoma.
Abbreviations: TC = transcriptional component. CI = confidence interval.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp8-v1.xlsx
-
Supplementary file 9
Expression and function of the top 20 genes in TC121, according to literature.
Abbreviations: TC = transcriptional component.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp9-v1.docx
-
Supplementary file 10
Colocalization scores of every two TC combinations in spatial transcriptomic profiles from three ovarian cancer samples.
Abbreviations: TC = transcriptional component.
- https://cdn.elifesciences.org/articles/101369/elife-101369-supp10-v1.xlsx
-
MDAR checklist
- https://cdn.elifesciences.org/articles/101369/elife-101369-mdarchecklist1-v1.docx
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Transcriptional pattern enriched for synaptic signaling is associated with shorter survival of patients with high-grade serous ovarian cancer
eLife 13:RP101369.
https://doi.org/10.7554/eLife.101369.3
