Dataset Information

Reducing confounding and suppression effects in TCGA data: an integrated analysis of chemotherapy response in ovarian cancer.

ABSTRACT:

Background

Despite initial response in adjuvant chemotherapy, ovarian cancer patients treated with the combination of paclitaxel and carboplatin frequently suffer from recurrence after few cycles of treatment, and the underlying mechanisms causing the chemoresistance remain unclear. Recently, The Cancer Genome Atlas (TCGA) research network concluded an ovarian cancer study and released the dataset to the public. The TCGA dataset possesses large sample size, comprehensive molecular profiles, and clinical outcome information; however, because of the unknown molecular subtypes in ovarian cancer and the great diversity of adjuvant treatments TCGA patients went through, studying chemotherapeutic response using the TCGA data is difficult. Additionally, factors such as sample batches, patient ages, and tumor stages further confound or suppress the identification of relevant genes, and thus the biological functions and disease mechanisms.

Results

To address these issues, herein we propose an analysis procedure designed to reduce suppression effect by focusing on a specific chemotherapeutic treatment, and to remove confounding effects such as batch effect, patient's age, and tumor stages. The proposed procedure starts with a batch effect adjustment, followed by a rigorous sample selection process. Then, the gene expression, copy number, and methylation profiles from the TCGA ovarian cancer dataset are analyzed using a semi-supervised clustering method combined with a novel scoring function. As a result, two molecular classifications, one with poor copy number profiles and one with poor methylation profiles, enriched with unfavorable scores are identified. Compared with the samples enriched with favorable scores, these two classifications exhibit poor progression-free survival (PFS) and might be associated with poor chemotherapy response specifically to the combination of paclitaxel and carboplatin. Significant genes and biological processes are detected subsequently using classical statistical approaches and enrichment analysis.

Conclusions

The proposed procedure for the reduction of confounding and suppression effects and the semi-supervised clustering method are essential steps to identify genes associated with the chemotherapeutic response.

SUBMITTER: Hsu FH

PROVIDER: S-EPMC3481440 | biostudies-literature | 2012

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Reducing confounding and suppression effects in TCGA data: an integrated analysis of chemotherapy response in ovarian cancer.

Hsu Fang-Han FH Serpedin Erchin E Hsiao Tzu-Hung TH Bishop Alexander J R AJ Dougherty Edward R ER Chen Yidong Y

BMC genomics 20121026

<h4>Background</h4>Despite initial response in adjuvant chemotherapy, ovarian cancer patients treated with the combination of paclitaxel and carboplatin frequently suffer from recurrence after few cycles of treatment, and the underlying mechanisms causing the chemoresistance remain unclear. Recently, The Cancer Genome Atlas (TCGA) research network concluded an ovarian cancer study and released the dataset to the public. The TCGA dataset possesses large sample size, comprehensive molecular profil ...[more]

PMID: 23134756

Similar Datasets

Project description:BackgroundOvarian cancer is the most deadly gynecological cancer with a very poor prognosis. Xenograft mouse models have proven to be one very useful tool in testing candidate therapeutic agents and gene function in vivo. In this study we identify genes and gene networks important for the efficacy of a pre-clinical anti-tumor therapeutic, MT19c.MethodsIn order to understand how ovarian xenograft tumors may be growing and responding to anti-tumor therapeutics, we used genome-wide mRNA expression and DNA copy number measurements to identify key genes and pathways that may be critical for SKOV-3 xenograft tumor progression. We compared SKOV-3 xenografts treated with the ergocalciferol derived, MT19c, to untreated tumors collected at multiple time points. Cell viability assays were used to test the function of the PPARγ agonist, Rosiglitazone, on SKOV-3 cell growth.ResultsThese data indicate that a number of known survival and growth pathways including Notch signaling and general apoptosis factors are differentially expressed in treated vs. untreated xenografts. As tumors grow, cell cycle and DNA replication genes show increased expression, consistent with faster growth. The steroid nuclear receptor, PPARγ, was significantly up-regulated in MT19c treated xenografts. Surprisingly, stimulation of PPARγ with Rosiglitazone reduced the efficacy of MT19c and cisplatin suggesting that PPARγ is regulating a survival pathway in SKOV-3 cells. To identify which genes may be important for tumor growth and treatment response, we observed that MT19c down-regulates some high copy number genes and stimulates expression of some low copy number genes suggesting that these genes are particularly important for SKOV-3 xenograft growth and survival.ConclusionsWe have characterized the time dependent responses of ovarian xenograft tumors to the vitamin D analog, MT19c. Our results suggest that PPARγ promotes survival for some ovarian tumor cells. We propose that a combination of regulated expression and copy number can identify genes that are likely important for chemotherapy response. Our findings suggest a new approach to identify candidate genes that are critical for anti-tumor therapy.

Project description:BackgroundAbnormal expression of long non-coding RNAs (lncRNA) play a significant role in the incidence and progression of high-grade serous ovarian cancer (HGSOC), which is a leading cause of mortality among gynecologic malignant tumor patients. In this study, our aim is to identify lncRNA-associated competing endogenous RNA (ceRNA ) axes that could define more reliable prognostic parameters of HGSOC, and to investigate the lncRNAs' potential mechanism of in lymphocyte infiltration.MethodsThe RNA-seq and miRNA expression profiles were downloaded from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database; while for obtaining the differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs), we used edgeR, limma and DESeq2. After validating the RNA, miRNA and gene expressions, using integrated three RNA expression profiles (GSE18520, GSE27651, GSE54388) and miRNA profile (GSE47841) from the Gene Expression Omnibus (GEO) database, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analyses through ClusterProfiler. The prognostic value of these genes was determined with Kaplan-Meier survival analysis and Cox regression analysis. The ceRNA network was constructed using Cytoscape. The correlation between lncRNAs in ceRNA network and immune infiltrating cells was analyzed by using Tumor IMmune Estimation Resource (TIMER), and gene markers of tumor-infiltrating immune cells were identified using Spearman's correlation after removing the influence of tumor purity.ResultsA total of 33 DELs (25 upregulated and eight downregulated), 134 DEMs (76 upregulated and 58 downregulated), and 1,612 DEGs (949 upregulated and 663 downregulated) were detected that could be positively correlated with overall survival (OS) of HGSOC. With the 1,612 analyzed genes, we constructed a ceRNA network, which indicated a pre-dominant involvement of the immune-related pathways. Furthermore, our data revealed that LINC00665 influenced the infiltration level of macrophages and dendritic cells (DCs). On the other hand, FTX and LINC00665, which may play their possible roles through the ceRNA axis, demonstrated a potential to inhibit Tregs and prevent T-cell exhaustion.ConclusionWe defined several prognostic biomarkers for the incidence and progression of HGSOC and constructed a network for ceRNA axes; among which three were indicated to have a positive correlation with lymphocyte infiltration, namely: FTX-hsa-miR-150-5p-STK11, LINC00665-hsa-miR449b-5p-VAV3 and LINC00665-hsa-miR449b-5p-RRAGD.

Project description:BACKGROUND:Ovarian cancer remains a significant public health burden, with the highest mortality rate of all the gynecological cancers. This is attributable to the late stage at which the majority of ovarian cancers are diagnosed, coupled with the low and variable response of advanced tumors to standard chemotherapies. To date, clinically useful predictors of treatment response remain lacking. Identifying the genetic determinants of ovarian cancer survival and treatment response is crucial to the development of prognostic biomarkers and personalized therapies that may improve outcomes for the late-stage patients who comprise the majority of cases. METHODS:To identify constitutional genetic variations contributing to ovarian cancer mortality, we systematically investigated associations between germline polymorphisms and ovarian cancer survival using data from The Cancer Genome Atlas Project (TCGA). Using stage-stratified Cox proportional hazards regression, we examined >650,000 SNP loci for association with survival. We additionally examined whether the association of significant SNPs with survival was modified by somatic alterations. RESULTS:Germline polymorphisms at rs4934282 (AGAP11/C10orf116) and rs1857623 (DNAH14) were associated with stage-adjusted survival (p= 1.12e-07 and 1.80e-07, FDR q= 1.2e-04 and 2.4e-04, respectively). A third SNP, rs4869 (C10orf116), was additionally identified as significant in the exome sequencing data; it is in near-perfect LD with rs4934282. The associations with survival remained significant when somatic alterations. CONCLUSIONS:Discovery analysis of TCGA data reveals germline genetic variations that may play a role in ovarian cancer survival even among late-stage cases. The significant loci are located near genes previously reported as having a possible relationship to platinum and taxol response. Because the variant alleles at the significant loci are common (frequencies for rs4934282 A/C alleles = 0.54/0.46, respectively; rs1857623 A/G alleles = 0.55/0.45, respectively) and germline variants can be assayed noninvasively, our findings provide potential targets for further exploration as prognostic biomarkers and individualized therapies.

Dataset Information

Reducing confounding and suppression effects in TCGA data: an integrated analysis of chemotherapy response in ovarian cancer.

Background

Results

Conclusions

Publications

Reducing confounding and suppression effects in TCGA data: an integrated analysis of chemotherapy response in ovarian cancer.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets