Multi-omics analysis Identifies PTTG1 as a prognostic biomarker associated with immunotherapy and chemotherapy resistance

doi:10.21203/rs.3.rs-4923978/v1

Download PDF

Research Article

Multi-omics analysis Identifies PTTG1 as a prognostic biomarker associated with immunotherapy and chemotherapy resistance

https://doi.org/10.21203/rs.3.rs-4923978/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 25 Oct, 2024

Read the published version in BMC Cancer →

You are reading this latest preprint version

Background

Pituitary tumor-transforming gene 1 (PTTG1) is an important gene in tumour development. However, the relevance of PTTG1 in tumour prognosis, immunotherapy response, and medication sensitivity in human pan-cancer has to be determined.

Methods

TIMER, GEPIA, the human protein atlas, GEPIA, TISCH2, and cBioportal examined the gene expression, protein expression, prognostic value, and genetic modification landscape of PTTG1 in 33 malignancies based on the TCGA cohort. The association between PTTG1 and tumour immunity, tumour microenvironment, immunotherapy response, and anticancer drug sensitivity was investigated using GSCA, TIDE, and CellMiner CDB. Molecular docking was used to validate the possible chemotherapeutic medicines for PTTG1. Additionally, siRNA-mediated knockdown was employed to confirm the probable role of PTTG1 in paclitaxel-resistant cells.

Results

PTTG1 is overexpressed and associated with poor survival in most tumors. Functional enrichment study revealed that PTTG1 is involved in the cell cycle and DNA replication. A substantial connection between PTTG1 expression and immune cell infiltration points to PTTG1's possible role in the tumour microenvironment. High PTTG1 expression is associated with tumour immunotherapy resistance. The process could be connected to PTTG1, which mediates T cell exhaustion and promotes cytotoxic T lymphocyte malfunction. Furthermore, PTTG1 was found to be substantially linked with sensitivity to several anticancer medications. Suppressing PTTG1 with siRNA reduced clone formation and migration, implying that PTTG1 may play a role in paclitaxel resistance.

Conclusion

PTTG1 shows potential as a cancer diagnostic, prognostic, and chemosensitivity marker. Increased PTTG1 expression is linked to resistance to cancer treatment. The mechanism could be linked to PTTG1's role in promoting cytotoxic T lymphocyte dysfunction and mediating T cell exhaustion. It is feasible to consider PTTG1, which is expressed on Treg and Tprolif cells, as a new therapeutic target for overcoming immunotherapy resistance.

PTTG1

Tumor immunity

Immunotherapy

Drug sensitivity

Pan-cancer

The pituitary tumor-transforming gene 1 (PTTG1), also known as SECURIN was isolated from pituitary tumors and identified as a gene with tumorigenic effect in vivo. It encodes a homolog of the yeast Securin protein and acts as a sister chromatid regulator that inhibits separase function and thus prevents sister chromatid separation.^1,2. As a transcriptional activator, PTTG1 can activate C-Myc and up-regulate the expression of cyclinB1 and CDK1 ³. PTTG1 plays a role in chromosomal stability maintenance, orderly cell cycle monitoring, genetic stability and mitotic fidelity assurance, malignant transformation, tumour induction, invasion, and metastasis^4,5.

Research studies indicate that overexpression of PTTG1 is associated with unfavorable tumor phenotype and adverse prognosis. This shows that PTTG1 plays a critical role in the occurrence and progression of tumours. For instance, hepatocellular carcinoma has an elevated overexpression of PTTG1, which enhances the synthesis of asparagine to activate mTOR and advance the tumour⁶. Furthermore, it was shown that a higher level of PTTG1 was associated with poorer prognosis for LUAD patients⁷, and targeted PTTG1 inhibition improves radiation-induced antitumour immunity in lung adenocarcinoma⁸. Additionally, the predictive values of PTTG1 for immunotherapy response⁸ and chemosensitivity ⁹ were also revealed. Although the importance of PTTG1 in tumors is becoming more widely known, there have been no comprehensive studies of its function in tumor occurrence, recurrence, immunotherapy or chemotherapy.

In this work, we evaluated the gene expression, prognostic value, and relate to immunotherapy response of PTTG1 in many well-established databases. We discovered that PTTG1 is overexpressed in several tumours and is related with poor survival in a variety of malignancies. PTTG1 promotes cytotoxic T cell malfunction, which contributes to immunotherapy resistance. PTTG1 was found to be significantly related with sensitivity to several anticancer medicines, particularly MET inhibitors. The current findings highlighted the critical function of PTTG1 in carcinogenesis, immunotherapy, and chemotherapy response, identifying a new target for cancer treatment.

Gene and protein expression analysis

The differential expression of PTTG1 in tumor and neighboring normal tissues across all TCGA tumors was analyzed by TIMER (Tumour Immune Estimation Resource, http://timer.cistrome.org/)¹⁰. The Wilcoxon test was used to determine the statistical significance of differential expression. The PTTG1 gene expression profile and pathological stage plot across all tumor samples and paired normal tissues were analyzed by GEPIA2 (Gene Expression Profiling Interactive Analysis 2, http://gepia2.cancer-pku.cn/#index)¹¹. The immunohistochemistry-based protein expression profiles of PTTG1 were analyzed by the Human Protein Atlas (hpa, https://www.proteinatlas.org/)¹².

Kaplan–Meier Test

Survival analysis of PTTG1 in tumor patients was obtained from the GEPIA2 based on the TCGA datasets. Kaplan-Meier Plotter ¹³ (http://kmplot.com/analysis/) was used to validate the prognostic value of PTTG1 in NSCLC patients. The Kaplan-Meier method was performed to plot the overall survival (OS) or disease-free survival (DFS, also called relapse-free survival and RFS) curves. Survival differences were assessed by the log-rank test using the median of PTTG1 as a cutoff value. The cox proportional hazard ratio (HR) with 95% confidence intervals and log‑rank P‑values were calculated.

Genetic alteration analysis

The PTTG1 alteration landscape in human cancers was depicted based on TCGA pan-cancer atlas studies by cBioPortal (https://www.cbioportal.org/) ¹⁴. To determine the correlation with copy number variation and survival, patients were divided into delete, amplification, and wild-type groups based on PTTG1 mutation type. The Kaplan-Meier method was performed to plot the overall survival curves in TCGA datasets by GSCA (Gene Set Cancer Analysis, https://guolab.wchscu.cn/GSCA/#/) ¹⁵.

Similar Genes Detection and PPI network construction

To identify genes with a similar expression pattern to PTTG1 in pan-cancer, the Pearson correlation coefficients were analyzed using the GEPIA2 database. The similar genes were analyzed using the STRING database (https://string-db.org/) to generate the PPI network. PPI pairs with a combined score > 0.4 were extracted. The PPI network was visualized using Cytoscape 3.7.2 software ¹⁶, and the most important module was performed using the MCODE plug-in in Cytoscape software. The Metascape web-based tool (https://metascape.org/gp/index.html) ¹⁷ was used for functional enrichment analysis (GO analysis and KEGG pathway analysis) of these genes.

PTTG1 and tumor immune microenvironment

The correlation of PTTG1 expression with immune infiltration level in diverse cancer types was analyzed by TIMER (Tumor IMmune Estimation Resource, https://cistrome.shinyapps.io/timer/)¹⁰. The TCGA pan-cancer data set was retrieved from the UCSC database. The PTTG1 expression data in individual samples was retrieved, and a log2(x + 1) transformation was applied. The immune, stromal, and ESTIMATE scores for each sample were determined using the R package ESTIMATE (https://bioinformatics.mdanderson.org/estimate/rpackage.html). Pearson correlation analysis was used to identify significantly linked immune infiltration scores.

PTTG1 and immunotherapy response

TIDE (Tumor Immune Dysfunction and Exclusion, http://tide.dfci.harvard.edu/) is a computational framework created to assess the potential of tumour immune escape from the gene expression profiles of cancer samples¹⁸. The 'Biomarker Evaluation' module was used to compare PTTG1 to other reported biomarkers in terms of their predictive ability for response outcome and overall survival. The 'Query Gene' module calculated the PTTG1 gene signature in T-cell dysfunction using data from the TCGA, PRECOG, and METABRIC databases. The T cell dysfunction score of PTTG1 is defined as the Wald test z score ¹⁸. The two-sided Wald test was used to calculate the relationship between cytotoxic T lymphocyte CTL levels and overall survival. The Kaplan-Meier plot divides tumours into two groups: 'High CTL' has above-average CTL values across all samples, and 'Low CTL' has values below average. To demonstrate the relationship between CTL levels and survival outcomes, samples were divided based on their PTTG1 expression levels.

Expression level of PTTG1 at the single‑cell level

TISCH2 (Tumor Immune Single-cell Hub 2, http://tisch.comp-genomics.org/) is a scRNA-seq database that focusses on the tumour microenvironment. It provides precise cell-type annotation at the single-cell level, allowing the analysis of TME across various cancer types¹⁹. The expression level of PTTG1 at the single-cell level was visualized by TISCH2 in the NSCLC_GSE139555, KIRC_GSE111360, and LIHC_GSE140228 datasets. The PTTG1 expression between immunotherapy responders and non-responders was investigated in the LIHC_GSE125449_aPDL1aCTLA4 dataset.

PTTG1 and chemo response

CellMiner Cross-Database (CellMinerCDB, discover.nci.nih.gov/cellminercdb)²⁰ enables the analysis of molecular and pharmacological data across cancer cell line databases to identify medications that match genomic determinants of response. CellMinerCDB investigated the relationship between PTTG1 expression and chemotherapy response using the GDSC (Sanger/Massachusetts General Hospital Genomics of Drug Sensitivity in Cancer) dataset.

Moleular Docking

The X-ray crystal structure of PTTG1 (7nj0) was obtained from the RCSB PDB protein databank (http://www.rcsb.org/). The 3D formats of ABT-263 (ZINC150338726) and NSC-207895 (ZINC5180959) were obtained from ZINC15 drug database (http://zinc15.docking.org/)²¹. Discovery Studio 4.5 Client deleted the water and ligand molecules from the crystal structures before starting the docking simulation. PyRx software was used to simulate drug docking and PTTG1 interactions. We calculated interaction energies to predict docking positions and pick the binding pose with the lowest binding energy (kcal mol^− 1). The results were visualised and analysed with Discovery Studio 4.5 Client.

Cell lines

Human lung cancer cell lines A549, H358 and the paclitaxel resistant cell lines A549-TXR, H358-TXR were kindly provided by Professor Wang Luo (University of Michigan, USA). Cells were maintained in RPMI 1640 or F12K supplemented with 10% FBS and 1% antibiotic-antimycotic. All cell lines were cultured at 37°C in a 5% CO2 cell culture incubator. Mycoplasma contamination was excluded in these cell lines.

Real-Time Quantitative PCR

Total RNA was reverse transcribed into cDNA using Reverse Transcription Kit (Thermo Fisher Scientific). qRT-PCR was performed using Power SYBR™ Green (Thermo Fisher Scientific) and ABI 7300 detection system (Applied Biosystems). qRT-PCR data were normalized to the expression of housekeeping gene GAPDH, and relative expressions were calculated using the 2^−ΔΔCt method. The oligonucleotide primers were shown as follows: PTTG1 (forward: 5’- ATGAATGCGGCTGTTAAGACCTG-3’, reverse: 5’- TCCCATCTAAGGCTTTGATTGAAGG-3’). All tests were performed in triplicate, and the data were presented as mean ± SD.

siRNA-mediated knockdown

Cells were seeded at the desired concentration in 60 mm plates and then transfected with 10 nM experimental siRNA oligonucleotides or non-targeting controls 24 hours later. Lipofectamine® RNAiMax Reagent (Invitrogen, USA) was used to knockdown cells in OptiMEM medium, following the manufacturer's recommendations. The knockdown efficiency was measured using qPCR. PTTG1 siRNA sequences for knockdown tests are listed below: sense: 5'-UGUGGUUGCUAAGGAUGGGCUTT-3'; antisense: 5'- AGCCCAUCCUUAGCAACCACATT-3'.

Clone Formation Assay

A549-TXR and H358-TXR cells (250 per well) were planted in 12-well plates and grown overnight. Cells were transfected with 10 nM experimental siRNA oligonucleotides or non-targeting controls 24 hours after plating. After two weeks, the cells were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet. Each well's clone count was recorded.

Basement membrane migration assays

Cells were treated with the siRNAs listed above to perform invasion tests. After 48 hours of transfection, cells were trypsinised, and diluted to the desired concentration. Cells were seeded into basement membrane matrix Boyden chambers (8-mm pore size, BD) located in the insert of a 24-well culture plate. The lower compartment received 20% FBS as a chemoattractant. After 48 hours, the non-migrating cells were gently removed using a cotton swab. Cells on the chamber's lower side were stained with Diff-QuikTM Stain Set (SIEMENS), then air dried and photographed.

Statistical analysis

ANOVA was used to compare PTTG1 expression levels in tumour and normal samples from the TCGA and GTEx datasets. Pearson's correlation coefficient was used to analyse the correlation between Immuno Score, Stromal Score, and ESTIMATE Score, as well as immunological checkpoints. The P-value < 0.05 was judged significant.

PTTG1 is overexpressed in pan cancer and associated with tumor stages

We performed a gene expression examination for PTTG1 on the gene expression matrix derived from the TCGA datasets using TIMER. PTTG1 mRNA expression was considerably higher in BLCA, BRCA, CHOL, ESCA, HNSC, KICH, KIRC, KIRP, LIHC, LUAD, LUSC, PRAD, READ, STAD, and UCEC, but significantly lower in THCA. The remaining cancer types exhibited no significant variation in PTTG1 mRNA levels between normal and tumor cells (Fig. 1A). The PTTG1 mRNA expression in pan-cancer was then analyzed using RNA sequencing datasets from the TCGA and GTEx projects by GEPIA. In addition to the 14 cancer types listed above, PTTG1 was significantly overexpressed in ACC, CESC, DLBC, GBM, OV, PAAD, SKCM, THYM, and UCS while being underexpressed in LAML and TGCT (Fig. 1B). To validate the diagnostic and prognostic relevance of PTTG1 in pan cancer, we used GEPIA to investigate its expression in tumor pathological stages. The findings revealed that PTTG1 was overexpressed in samples from advanced cancer patients with KIRC, KIRP, ACC, LUAD, LIHC, and BRCA (Fig. 1C). We then obtained the immunohistochemistry staining of PTTG1 from the HPA database. The PTTG1 protein stained darker in LUAD, BRCA, LIHC, and KIRC samples (Fig. 1D). The findings indicate that the PTTG1 protein may be significantly expressed in these tumors.

PTTG1 is a predictor of poor patient survival

To investigate the prognostic benefit of PTTG1, we compared survival (OS and DFS) differences between high and low PTTG1 expression groups in pan-cancer. The high PTTG1 expression group had substantially lower overall survival (P < 0.05) compared to the low expression group in ACC, KIRC, KIRP, LGG, LIHC, LUAD, MESO, PAAD, THCA, and UVM (Fig. 2A). Increased PTTG1 expression was linked to decreased disease-free survival in ACC, KIRC, KIRP, LGG, LIHC, MESO, PAAD, PRAD, SARC, and UVM (Fig. 2B). These findings indicate that PTTG1 is overexpressed and linked with poor outcomes in ACC, KIRC, KIRP, LIHC, LUAD, and PAAD, and has the potential to be a predictive biomarker.

Associations of PTTG1 expression with genomic alterations in pan-cancer

cBioPortal was used to analyse genetic variations (mutation, structural variant, amplification, and deep deletion) in PTTG1 across malignancies. In the majority of cancer types, PTTG1 gene amplification was obvious. High amplification was observed in KIRC, CHOL, and UCS (Supplementary Figs. 1A, B). The GSCA database was then used to analyse the frequency of change in single nucleotide variation (SNV) and copy number variation (CNV) in the PTTG1 gene. The findings revealed that SNV of PTTG1 is more frequently found in UCEC (0.88%, Fig. 3B) and STAD (0.68%, Fig. 3B), with low frequency in the other cancer types (Fig. 3A). The most common detrimental mutation (Missense_Mutation, Splice_Site, Frame_Shift_Del, etc.) in the PTTG1 gene is missense mutation. Point mutation analysis revealed that the SNV of PTTG1 had C > T and C > A transversions (Fig. 3C). Next, we looked for CNV of the PTTG1 gene across malignancies. The findings revealed that CNV of PTTG1 was common in KIRC, ACC, CHOL, BLCA, TGCT, and LUSC, but uncommon in THCA, LAML, and PRAD (Fig. 3D). We investigated at the Kaplan-Meier survival curves for PTTG1 CNV (homo amplification and homo deletion) and broad type (Fig. 3E). KIRP patients with homo deletion and homo amplification of PTTG1 genes exhibited a poor prognosis compared to the wild type group (P < 0.05, Fig. 3F).

Protein-protein interactions of PTTG1 and similar genes in pan-cancer

In order to investigate the probable mechanism of PTTG1 in carcinogenesis, we identified the top 100 genes with similar expression patterns to PTTG1 in the TCGA cancer cohort using GEPIA2. The top five genes with similar expression patterns ordered by Pearson correlation coefficient were AURKB (R = 0.69, P < 0.0001), CDC20 (R = 0.70, P < 0.0001), CCNB1 (R = 0.69, P < 0.0001), KIF2C (R = 0.66, P < 0.0001), and RAD54L (R = 0.64, P < 0.0001) (Fig. 4A). Metascape was used to visualise and analyse the interactome network composed of the top 100 PTTG1 related genes. The PPI network has 90 nodes and 599 edges, indicating a wide range of interactions between these proteins (Fig. 4B). The MCODE algorithm was then applied to this network to find sites where proteins are highly linked. Following that, each MCODE network underwent GO enrichment analysis. The results revealed that five MCODE complexes were discovered in the PPI network, with 'cell cycle' being the most common biological meaning (Fig. 4C). To further study the functions of the top 100 PTTG1 similar genes, GO enrichment analysis and KEGG pathway analysis were used. These proteins have a significant role in various cellular processes, including nuclear division, organelle fission, microtubule binding, and tubulin binding (P < 0.0001) (Fig. 4D). Pathway enrichment revealed that PTTG1-related genes function in multiple pathways, including cell cycle and DNA replication (P < 0.0001) (Fig. 4D). These findings indicated that cancer-related genes and pathways were common in PTTG1 similar genes.

Correlations between PTTG1 expression and immune cell infiltration in pan-cancer

We investigated the relationship between immune cell infiltration and PTTG1 gene expression in 33 malignancies using GSCA. ImmuCellAI was performed to get the infiltrates of immune cells in each TCGA sample. In most types of cancer, PTTG1 correlated positively with T cell exhaustion, Th1, B cells, DC, effector memory T cells, Infiltration Score, CD8 naive T cells, CD8 T cells, cytotoxic T cells, and nTreg, but negatively with CD4_T, CD4_naive, central memory T cells, MAIT, neutrophil, NKT, and Th17 cells (Fig. 5A, Supplementary table 1). The ESTIMATE package in R was used to determine the connection between ESTIMATEScore, ImmunoScore, and StromalScore and PTTG1 expression in TCGA pan-cancer datasets. PTTG1 expression showed a positive correlation with ImmunoScore, StromalScore, and ESTIMATEScore in GBMLGG and KIPAN (Pearson R > 0.2 and Pearson P value < 0.05). PTTG1 expression had a negative connection with ImmunoScore, StromalScore, and ESTIMATEScore in GBM, TGCT, READ, STES, and LUSC (Pearson R<-0.2 and Pearson P value < 0.05, Fig. 5B, Supplementary table 2). The connection of PTTG1 with immune cell infiltration may help to explain the poor prognosis caused by its high expression.

The expression level of PTTG1 is associated with response to immunotherapy treatment

To investigate the possibility of PTTG1 as an immune checkpoint blockade response biomarker, the relationship between PTTG1 and immunotherapy response was examined using the TIDE database. The results showed that high expression of PTTG1 significantly affected the efficacy of immune checkpoint blockade (anti-PD1) and adoptive T cell therapy (ACT), reducing the OS of patients in the Braun2020_PD1_Kidney_Clear cohort, Lauss2017_ACT_Melanoma cohort, and Riaz017_PD1_Melanoma_lpi. Prog cohort (Fig. 6A.B). To assess the accuracy of PTTG1 as an immunotherapy response biomarker, we examined it with other biomarkers previously associated with tumour immune evasion by TIDE. The area under the receiver operating characteristic curve (AUC) was used to assess the predictive ability of these biomarkers. PTTG1 produced an AUC greater than or equal to 0.5 in 11 of the 16 immune checkpoint blockade sub-cohorts (Fig. 6C). PTTG1 outperformed previously published biomarkers in predicting immunotherapy prognosis for melanoma patients (Nathanson2017_CTLA4_Melanoma_Post cohort, AUC = 0.8636 and Gide2019_PD1 + CTLA4_Melanoma cohort, AUC = 0.7000) (Fig. 6C). These findings revealed that PTTG1 played a role in antitumor immune response and promoted immunotherapy resistance.

PTTG1 exacerbates CTL dysfunction and contributes to resistance against immunotherapy interventions

The relationship between PTTG1 and Cytotoxic T lymphocyte (CTL) dysfunction was evaluated by TIDE. The results showed that a higher CTL level indicates better patient survival, but only when PTTG1 has a low expression level in glioma (OS, z = 2.81, P = 0.00501), myeloma (OS, z = 3.66, P = 0.00025), LUAD (OS, z = 2.25, P = 0.0241), LUSC (OS, z = 2.02, P = 0.0432), COAD (OS, z = 3.08, P = 0.0021), and KIRC (PFS, z = 2.75, P = 0.00604) (Fig. 7A and Supplementary table 3). This finding suggests that a higher PTTG1 level in tumours reduces the positive connection between CTL and survival. Gene co-expression study revealed that the expression of PTTG1 was positively correlated with immune inhibitory marker such as EDNRB and C10orf54. And negatively correlated with immune stimulatory marker such as HMGB1, CD70, TNFSF9 and TNFRSF18 (Fig. 7B). Tumour mutation burden (TMB) and homologous recombination deficit (HRD) are two clinically important immunological markers that are closely related to tumour immunotherapy. The link between PTTG1 expression and TMB or HRD scores was examined utilising the Sangerbox platform. We discovered a strong positive connection between PTTG1 mRNA expression and TMB or HRD scores in the majority of malignant tumours, particularly KICH (Fig. 7C and D). These results suggested that PTTG1 was related to promoting tumor immune escape and immunotherapy resistance.

PTTG1 mainly expressed in regulatory T cells and proliferative T cells

To investigate the probable ways by which PTTG1 affects theted tumor immune microenvironment, public single-cell RNA-seq (scRNA) datasets were used to determine the expression level of PTTG1 on various immune cells. We discovered that PTTG1 was extensively expressed in immune cells, but it was particularly high in regulatory T cells (Treg) and proliferative T cells (Tprolif) in NSCLC, KIRC, and LIHC cohots (Fig. 8A, B, C). Interestingly, immunotherapy responders' hepatic progenitor cells and malignant cells expressed PTTG1 more than non-responders in LIHC patients, indicating that PTTG1 was related with immune suppressive cells and immunotherapy resistance (Fig. 8D).

Cellular clusters and single-cell PTTG1 expression profile in different cellular types of NSCLC (A), KIRC (B), and LIHC (C). (D) Single-cell PTTG1 expression profile of non-responder and responder in LIHC patients with PD-L1/CTLA-4 treatment.

Correlation of PTTG1 and drug response

The CellMinerCDB was used to look into the relationship between PTTG1 expression and drug sensitivity (IC50) in the GDSC dataset. Pearson's correlation analysis demonstrated a negative connection between high PTTG1 expression and the -log10(IC50) values of 22 drugs (Fig. 9A, Supplementary table 4), particularly MEK1/MEK2 inhibitors such as RDEA119, trametinib, CI-1040, selumetinib, and PD-0325901 (Fig. 9B). These data showed that increasing PTTG1 could reduce drug sensitivity for these small compounds. PTTG1 showed a substantial positive correlation with the -log10(IC50) values of 13 drugs (Fig. 9A, Supplementary table 4), including NSC-207895 (Fig. 9C) and ABT-263 (Navitoclax) (Fig. 9E). The results suggested that these 13 small compounds have the potential to be PTTG1 target medicines. A molecular docking model was developed to investigate the mechanism of interaction between PTTG1 and prospective target medicines. We discovered that NSC-207895 could bind to PTTG1 and remain in the binding pocket surrounded by critical residues (Lys1189, Pro1144, Leu1140, His1142, Cys1160, and Arg1197) (Fig. 9D). ABT-263 could attach to PTTG1 and remain in the binding pocket surrounded by critical residues (Ala1671, Arg1675, Cys928, Glu925, Ala929, Leu945, Arg1638, Lys1645, Asp941, Gly931, Gln933, and Gln1674) (Fig. 9F). Our findings suggested that PTTG1 could be a therapeutic target for MEK1/MEK2 inhibitor resistance. ABT-263 and NSC-207895 may be effective anticancer drugs that target PTTG1.

Validation of PTTG1 expression and function in NSCLC

We further confirmed the expression of PTTG1 in tumours by utilising the GEO database. The findings showed that NSCLC tissues had considerably greater levels of PTTG1 expression than did normal lung tissues (P < 0.01, Fig. 10A). Additionally, we used the online Kaplan Meier plotter database to determine the predictive significance of PTTG1 in 1161 patients with LUAD. Poor patient survival was observed to be substantially correlated with high expression of PTTG1 (low vs. high: 110.27 vs. 60.73, P < 0.01) (Fig. 10B). In the multi-drug resistance lung cancer dataset GSE77209, which included the parental cell lines H1299, H1355, and paclitaxel-carboplatin-resistant cells, we examined the expression of PTTG1. In H1299 and H1355, resistant cells expressed more PTTG1 than the parental cells did (Fig. 10C). The expression of PTTG1 in chemoresistance cells was confirmed by qRT-PCR. The expression levels of PTTG1 in the paclitaxel-resistant cell lines A549-TXR and H358-TXR were much higher than those in the parental cells, A549 and H358 (Fig. 10D). We employed short interfering RNAs to knock down PTTG1 in A549- TXR and H358- TXR in order to ascertain the biological roles of PTTG1 in paclitaxel-resistant cells (Fig. 10D). The findings showed that PTTG1 knockdown reduced A549-TXR and H358-TXR's capacity to generate clones (Fig. 10F). The transwell experiment verifies that PTTG1 knockdown prevents A549-TXR and H358-TXR cells from migrating (Fig. 10G). These results imply that PTTG1 is involved in both chemoresistance and carcinogenesis.

Human life is seriously threatened by cancer. It is useful to employ biomarkers that are often expressed in many tumour types as targets for treatment and as tumour biomarkers. PTTG1 is involved in the growth of tumours and may function as an oncogene in the development and growth of tumours ²². We analyzed the role of PTTG1 in pan-cancer and confirmed that PTTG1 was overexpressed in 15 of 33 tumor types of TCGA datasets. The expression of PTTG1 was related to tumor classification in several cancers like ACC, LUAD, LIHC, BRCA, KIRC and KIRP. Moreover, high PTTG1 expression was related to poor patient survival, consistent with the previously reported role of PTTG1 in kidney renal clear cell carcinoma ²³ and hepatocellular carcinoma ⁶. An increasing number of research have investigated the relationship between genetic changes and cancer progression. We discovered that PTTG1 gene amplification resulted in significant protein expression in most tumours. Pathway enrichment revealed that PTTG1 comparable genes were involved in a variety of pathways, including cancer-related pathways such as 'cell cycle' and 'DNA replication'. These findings imply that PTTG1 has the potential to serve as a tumour diagnostic and prognostic marker. Furthermore, inhibiting aberrant PTTG1 expression at the genetic and protein levels could be a possible therapeutic method for reversing carcinogenesis.

Tumor-produced chemicals and non-cancerous elements make up the tumour microenvironment. It is essential for the initiation, growth, metastasis, and reaction to treatment of tumours. One of the main mechanisms of tumour cell immune escape that leads to immunological dysfunction in cancer patients is T cell depletion^24,25. Recent research has shown that T cell exhaustion is considered a mechanism of resistance for cellular immunotherapies ²⁶ and a significant marker of poor outcomes in many cancers including breast cancer ²⁷ and renal cell carcinoma ²⁸. Revitalization of exhausted T cells might improve immunity. In this context, recovered fatigued T cells may be a useful source of predictive biomarkers for a possible target. revitalising exhausted T cells might improve immunity. In this context, recovered fatigued T cells may be a useful source of predictive biomarkers for a possible target. In this study, we discovered that PTTG1 is positively connected with T cell exhaustion but negatively correlated with the ImmunoScore, StromalScore, and ESTIMATEScore in most types of tumours. Furthermore, PTTG1 was shown to be strongly expressed in immune-modulating cells such as Treg and Tprolif cells from NSCLC, KIRC, and LIHC patients. Treg cells produce an immuno-suppressive environment and are less sensitive to immune checkpoint inhibitors^29,30. It was noteworthy that Tprolif cells were able to exhibit high expression levels of immunological exhaustion markers, including PDCD1, HAVCR2, CTLA4, LAG3, and TIGIT³¹. In both KIRC and LIHC, our results demonstrated a strong positive correlation between PTTG1 and the immune-suppressive genes HAVCR2, CTLA4, LAG3, TIGIT, and PDCD1. These findings imply that the immune-suppressive microenvironment and T cell exhaustion may be connected to the elevated expression of PTTG1.

The immuno-suppressive tumour microenvironment makes immune checkpoint drugs significantly less effective than anticipated in the treatment of cancer³². To investigate PTTG1's involvement in immunotherapy resistance, the link between PTTG1 and immunotherapy response was examined. According to the results, patients with KIRC and melanoma undergoing immune checkpoint inhibition (anti-PD1) and adoptive T cell treatment (ACT) had a worse prognosis when their expression of PTTG1 was elevated. An important part of the antitumor immune system is played by cytotoxic T lymphocytes (CTLs). Immunotherapy resistance and tumour immune evasion are significantly aided by CTL malfunction. We discovered that higher CTL levels predicted longer patient life only when PTTG1 expression was low in glioma, myeloma, NSCLC, COAD, and KIRC. As a result, increased PTTG1 levels in tumours will reduce the positive correlation between CTL and survival. TMB and HRD are prognostic indicators for immune checkpoint inhibitor clinical effectiveness³³. Particularly in KICH, we discovered a significant positive relationship between PTTG1 and TMB or HRD scores in most malignant tumours. These findings imply a relationship between tumour immunotherapy resistance and PTTG1 overexpression. The process might be connected to PTTG1's role in increasing cytotoxic T lymphocyte malfunction and mediating T cell exhaustion. It is possible to view PTTG1, which is expressed on Treg and Tprolif cells, as a novel therapeutic target to overcome immunotherapy resistance.

A increasing amount of research indicates that PTTG1 plays a crucial role in chemosensitivity³⁴. Pancreatic ductal adenocarcinoma patients with decreased PTTG1 have good treatment response with great sensitivity and selectivity³⁵. In ovarian cancer cell lines, downregulating PTTG1 improved saracatinib sensitivity³⁶. PTTG1 gene suppression reduces prostate cancer cell sensitivity to paclitaxel-induced apoptosis³⁷. Consistent with earlier research, we found that PTTG1 was linked to chemosensitivity, particularly in MEK1/MEK2 inhibitors such as PD-0325901, RDEA119, trametinib, CI-1040, and selumetinib. ABT-263 and NSC-207895 may be employed as an efficient anticancer medication by targeting PTTG1, according to the molecular docking studies. Furthermore, we confirmed the expression and function of PTTG1 in NSCLC. We discovered that PTTG1 was considerably greater in NSCLC tissues compared to non-tumor lung tissues. Furthermore, elevated PTTG1 expression was found to be substantially associated with poor patient survival in NSCLC. PTTG1 expression levels in paclitaxel-resistant cell lines A549-PTX and H358-PTX were much higher than in the parental cells A549 and H358, which is consistent with recent studies in LUAD³⁸. PTTG1 knockdown inhibited clone formation and migration ability of A549- TXR and H358- TXR cells.

To summarise, our investigations indicate that PTTG1 may have potential as a tumour diagnostic, prognostic, and chemosensitivity marker. Reversing carcinogenesis and chemoresistance may be achieved by blocking aberrant PTTG1 expression at the genetic and protein levels. Increased PTTG1 expression is associated with resistance to tumour treatment. The process might be connected to PTTG1's role in increasing cytotoxic T lymphocyte malfunction and mediating T cell exhaustion. It is possible to view PTTG1, which is expressed on Treg and Tprolif cells, as a novel therapeutic target to overcome immunotherapy resistance.

ACC

Adrenocortical carcinoma

BLCA

Bladder Urothelial Carcinoma BRCA,Breast invasive carcinoma

CESC

Cervical squamous cell carcinoma and endocervical adenocarcinoma

CHOL

Cholangio carcinoma

COAD

Colon adenocarcinoma

DLBC

Lymphoid Neoplasm Diffuse Large B-cell Lymphoma

ESCA

Esophageal carcinoma

GBM

Glioblastoma multiforme

HNSC

Head and Neck squamous cell carcinoma

KICH

Kidney Chromophobe

KIRC

Kidney renal clear cell carcinoma

KIRP

Kidney renal papillary cell carcinoma

LAML

Acute Myeloid Leukemia

LGG

Brain Lower Grade Glioma

LIHC

Liver hepatocellular carcinoma

LUAD

Lung adenocarcinoma

LUSC

Lung squamous cell carcinoma

MESO

Mesothelioma

Ovarian serous cystadenocarcinoma

PAAD

Pancreatic adenocarcinoma

PCPG

Pheochromocytoma and Paraganglioma

PRAD

Prostate adenocarcinoma

READ

Rectum adenocarcinoma

SARC

Sarcoma

SKCM

Skin Cutaneous Melanoma

STAD

Stomach adenocarcinoma

TGCT

Testicular Germ Cell Tumors

THCA

Thyroid carcinoma

THYM

Thymoma

UCEC

Uterine Corpus Endometrial Carcinoma

UCS

Uterine Carcinosarcoma

UVM

Uveal Melanoma.

Author contributions

Lihui Wang and Chunlin Zou designed the study. Handong Wei, Yashu Ma and Shuxing Chen collected the data. All Authors analyzed the data and wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Joint Project on Regional High-Incidence Diseases Research of Guangxi Natural Science Foundation under Grant (No. 2023GXNSFAA026338), the National Natural Science Foundation of China (81803564), and the National College Students Innovation and Entrepreneurship Training Program (No. 202210598040；No. 202210598038).

CONFLICT OF INTEREST

The authors declare that there are no conflicts of interest.

Ethics, Consent to Participate, and Consent to Publish declarations

Not applicable.

Yanagida, M. (2005). Basic mechanism of eukaryotic chromosome segregation. Philos Trans R Soc Lond B Biol Sci 360, 609-621. 10.1098/rstb.2004.1615.
Wirth, K.G., Wutz, G., Kudo, N.R., Desdouets, C., Zetterberg, A., Taghybeeglu, S., Seznec, J., Ducos, G.M., Ricci, R., Firnberg, N., et al. (2006). Separase: a universal trigger for sister chromatid disjunction but not chromosome cycle progression. J Cell Biol 172, 847-860. 10.1083/jcb.200506119.
Marangos, P., and Carroll, J. (2008). Securin regulates entry into M-phase by modulating the stability of cyclin B. Nat Cell Biol 10, 445-451. 10.1038/ncb1707.
Bernal, J.A., Luna, R., Espina, A., Lazaro, I., Ramos-Morales, F., Romero, F., Arias, C., Silva, A., Tortolero, M., and Pintor-Toro, J.A. (2002). Human securin interacts with p53 and modulates p53-mediated transcriptional activity and apoptosis. Nat Genet 32, 306-311. 10.1038/ng997.
Moreno-Mateos, M.A., Espina, A.G., Torres, B., Gamez del Estal, M.M., Romero-Franco, A., Rios, R.M., and Pintor-Toro, J.A. (2011). PTTG1/securin modulates microtubule nucleation and cell migration. Mol Biol Cell 22, 4302-4311. 10.1091/mbc.E10-10-0838.
Zhou, Q., Li, L., Sha, F., Lei, Y., Tian, X., Chen, L., Chen, Y., Liu, H., and Guo, Y. (2023). PTTG1 Reprograms Asparagine Metabolism to Promote Hepatocellular Carcinoma Progression. Cancer Res 83, 2372-2386. 10.1158/0008-5472.CAN-22-3561.
Bai, L., Li, L.H., Liang, J., and Li, E.X. (2022). Prognostic Significance of PTTG1 and Its Methylation in Lung Adenocarcinoma. J Oncol 2022, 3507436. 10.1155/2022/3507436.
Chen, Z., Cao, K., Hou, Y., Lu, F., Li, L., Wang, L., Xia, Y., Zhang, L., Chen, H., Li, R., et al. (2021). PTTG1 knockdown enhances radiation-induced antitumour immunity in lung adenocarcinoma. Life Sci 277, 119594. 10.1016/j.lfs.2021.119594.
Hong, H., Jin, Z., Qian, T., Xu, X., Zhu, X., Fei, Q., Yang, J., Sui, C., and Xu, M. (2021). Falcarindiol Enhances Cisplatin Chemosensitivity of Hepatocellular Carcinoma via Down-Regulating the STAT3-Modulated PTTG1 Pathway. Front Pharmacol 12, 656697. 10.3389/fphar.2021.656697.
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J.S., Li, B., and Liu, X.S. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res 77, e108-e110. 10.1158/0008-5472.CAN-17-0307.
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res 45, W98-W102. 10.1093/nar/gkx247.
Uhlen, M., Fagerberg, L., Hallstrom, B.M., Lindskog, C., Oksvold, P., Mardinoglu, A., Sivertsson, A., Kampf, C., Sjostedt, E., Asplund, A., et al. (2015). Proteomics. Tissue-based map of the human proteome. Science 347, 1260419. 10.1126/science.1260419.
Gyorffy, B. (2024). Transcriptome-level discovery of survival-associated biomarkers and therapy targets in non-small-cell lung cancer. Br J Pharmacol 181, 362-374. 10.1111/bph.16257.
de Bruijn, I., Kundra, R., Mastrogiacomo, B., Tran, T.N., Sikina, L., Mazor, T., Li, X., Ochoa, A., Zhao, G., Lai, B., et al. (2023). Analysis and Visualization of Longitudinal Genomic and Clinical Data from the AACR Project GENIE Biopharma Collaborative in cBioPortal. Cancer Res 83, 3861-3867. 10.1158/0008-5472.CAN-23-0816.
Liu, C.J., Hu, F.F., Xie, G.Y., Miao, Y.R., Li, X.W., Zeng, Y., and Guo, A.Y. (2023). GSCA: an integrated platform for gene set cancer analysis at genomic, pharmacogenomic and immunogenomic levels. Brief Bioinform 24. 10.1093/bib/bbac558.
Shannon, P., Markiel, A., Ozier, O., Baliga, N.S., Wang, J.T., Ramage, D., Amin, N., Schwikowski, B., and Ideker, T. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 13, 2498-2504. 10.1101/gr.1239303.
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A.H., Tanaseichuk, O., Benner, C., and Chanda, S.K. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 10, 1523. 10.1038/s41467-019-09234-6.
Jiang, P., Gu, S., Pan, D., Fu, J., Sahu, A., Hu, X., Li, Z., Traugh, N., Bu, X., Li, B., et al. (2018). Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med 24, 1550-1558. 10.1038/s41591-018-0136-1.
Sun, D., Wang, J., Han, Y., Dong, X., Ge, J., Zheng, R., Shi, X., Wang, B., Li, Z., Ren, P., et al. (2021). TISCH: a comprehensive web resource enabling interactive single-cell transcriptome visualization of tumor microenvironment. Nucleic Acids Res 49, D1420-D1430. 10.1093/nar/gkaa1020.
Luna, A., Elloumi, F., Varma, S., Wang, Y., Rajapakse, V.N., Aladjem, M.I., Robert, J., Sander, C., Pommier, Y., and Reinhold, W.C. (2021). CellMiner Cross-Database (CellMinerCDB) version 1.2: Exploration of patient-derived cancer cell line pharmacogenomics. Nucleic Acids Res 49, D1083-D1093. 10.1093/nar/gkaa968.
Sterling, T., and Irwin, J.J. (2015). ZINC 15--Ligand Discovery for Everyone. J Chem Inf Model 55, 2324-2337. 10.1021/acs.jcim.5b00559.
Wang, L., and Liu, X. (2024). Multi-Omics Analysis of the Oncogenic Value of Pituitary Tumor-Transforming Gene 1 (PTTG1) in Human Cancers. Front Biosci (Landmark Ed) 29, 87. 10.31083/j.fbl2902087.
Zhang, X., Ji, H., Huang, Y., Zhu, B., and Xing, Q. (2023). Elevated PTTG1 predicts poor prognosis in kidney renal clear cell carcinoma and correlates with immunity. Heliyon 9, e13201. 10.1016/j.heliyon.2023.e13201.
Wherry, E.J., and Kurachi, M. (2015). Molecular and cellular insights into T cell exhaustion. Nat Rev Immunol 15, 486-499. 10.1038/nri3862.
Pauken, K.E., and Wherry, E.J. (2015). Overcoming T cell exhaustion in infection and cancer. Trends Immunol 36, 265-276. 10.1016/j.it.2015.02.008.
Chow, A., Perica, K., Klebanoff, C.A., and Wolchok, J.D. (2022). Clinical implications of T cell exhaustion for cancer immunotherapy. Nat Rev Clin Oncol 19, 775-790. 10.1038/s41571-022-00689-z.
Terranova-Barberio, M., Pawlowska, N., Dhawan, M., Moasser, M., Chien, A.J., Melisko, M.E., Rugo, H., Rahimi, R., Deal, T., Daud, A., et al. (2020). Exhausted T cell signature predicts immunotherapy response in ER-positive breast cancer. Nat Commun 11, 3584. 10.1038/s41467-020-17414-y.
Ficial, M., Jegede, O.A., Sant'Angelo, M., Hou, Y., Flaifel, A., Pignon, J.C., Braun, D.A., Wind-Rotolo, M., Sticco-Ivins, M.A., Catalano, P.J., et al. (2021). Expression of T-Cell Exhaustion Molecules and Human Endogenous Retroviruses as Predictive Biomarkers for Response to Nivolumab in Metastatic Clear Cell Renal Cell Carcinoma. Clin Cancer Res 27, 1371-1380. 10.1158/1078-0432.CCR-20-3084.
Principe, D.R., Chiec, L., Mohindra, N.A., and Munshi, H.G. (2021). Regulatory T-Cells as an Emerging Barrier to Immune Checkpoint Inhibition in Lung Cancer. Front Oncol 11, 684098. 10.3389/fonc.2021.684098.
Kumagai, S., Togashi, Y., Kamada, T., Sugiyama, E., Nishinakamura, H., Takeuchi, Y., Vitaly, K., Itahashi, K., Maeda, Y., Matsui, S., et al. (2020). The PD-1 expression balance between effector and regulatory T cells predicts the clinical efficacy of PD-1 blockade therapies. Nat Immunol 21, 1346-1358. 10.1038/s41590-020-0769-3.
Yang, C., Guo, Y., Wu, Z., Huang, J., and Xiang, B. (2022). Comprehensive Analysis of Cuproptosis-Related Genes in Prognosis and Immune Infiltration of Hepatocellular Carcinoma Based on Bulk and Single-Cell RNA Sequencing Data. Cancers (Basel) 14. 10.3390/cancers14225713.
O'Donnell, J.S., Teng, M.W.L., and Smyth, M.J. (2019). Cancer immunoediting and resistance to T cell-based immunotherapy. Nat Rev Clin Oncol 16, 151-167. 10.1038/s41571-018-0142-8.
Wang, Y., Tong, Z., Zhang, W., Zhang, W., Buzdin, A., Mu, X., Yan, Q., Zhao, X., Chang, H.H., Duhon, M., et al. (2021). FDA-Approved and Emerging Next Generation Predictive Biomarkers for Immune Checkpoint Inhibitors in Cancer Patients. Front Oncol 11, 683419. 10.3389/fonc.2021.683419.
Vlotides, G., Eigler, T., and Melmed, S. (2007). Pituitary tumor-transforming gene: physiology and implications for tumorigenesis. Endocr Rev 28, 165-186. 10.1210/er.2006-0042.
Sahni, S., Nahm, C., Ahadi, M.S., Sioson, L., Byeon, S., Chou, A., Maloney, S., Moon, E., Pavlakis, N., Gill, A.J., et al. (2023). Gene expression profiling of pancreatic ductal adenocarcinomas in response to neoadjuvant chemotherapy. Cancer Med 12, 18050-18061. 10.1002/cam4.6411.
Nakachi, I., Helfrich, B.A., Spillman, M.A., Mickler, E.A., Olson, C.J., Rice, J.L., Coldren, C.D., Heasley, L.E., Geraci, M.W., and Stearman, R.S. (2016). PTTG1 Levels Are Predictive of Saracatinib Sensitivity in Ovarian Cancer Cell Lines. Clin Transl Sci 9, 293-301. 10.1111/cts.12413.
Castilla, C., Flores, M.L., Medina, R., Perez-Valderrama, B., Romero, F., Tortolero, M., Japon, M.A., and Saez, C. (2014). Prostate cancer cell response to paclitaxel is affected by abnormally expressed securin PTTG1. Mol Cancer Ther 13, 2372-2383. 10.1158/1535-7163.MCT-13-0405.
Dalvi, M.P., Wang, L., Zhong, R., Kollipara, R.K., Park, H., Bayo, J., Yenerall, P., Zhou, Y., Timmons, B.C., Rodriguez-Canales, J., et al. (2017). Taxane-Platin-Resistant Lung Cancers Co-develop Hypersensitivity to JumonjiC Demethylase Inhibitors. Cell Rep 19, 1669-1684. 10.1016/j.celrep.2017.04.077.

No competing interests reported.

Download PDF

Journal Publication

published 25 Oct, 2024

Read the published version in BMC Cancer →

Editorial decision: Revision requested
19 Aug, 2024
Editor assigned by journal
19 Aug, 2024
Submission checks completed at journal
19 Aug, 2024
First submitted to journal
16 Aug, 2024

You are reading this latest preprint version

Multi-omics analysis Identifies PTTG1 as a prognostic biomarker associated with immunotherapy and chemotherapy resistance

Status:

Journal Publication

Version 1

Abstract

Background

Methods

Results

Conclusion

Figures

Introduction

Materials and methods

Gene and protein expression analysis

Kaplan–Meier Test

Genetic alteration analysis

Similar Genes Detection and PPI network construction

PTTG1 and tumor immune microenvironment

PTTG1 and immunotherapy response

Expression level of PTTG1 at the single‑cell level

PTTG1 and chemo response

Moleular Docking

Cell lines

Real-Time Quantitative PCR

siRNA-mediated knockdown

Clone Formation Assay

Basement membrane migration assays

Statistical analysis

Results

PTTG1 is overexpressed in pan cancer and associated with tumor stages

PTTG1 is a predictor of poor patient survival

Associations of PTTG1 expression with genomic alterations in pan-cancer

Protein-protein interactions of PTTG1 and similar genes in pan-cancer

Correlations between PTTG1 expression and immune cell infiltration in pan-cancer

The expression level of PTTG1 is associated with response to immunotherapy treatment

PTTG1 exacerbates CTL dysfunction and contributes to resistance against immunotherapy interventions

PTTG1 mainly expressed in regulatory T cells and proliferative T cells

Correlation of PTTG1 and drug response

Validation of PTTG1 expression and function in NSCLC

Discussion

Abbreviations

Declarations

References

Additional Declarations

Supplementary Files

Status:

Journal Publication

Version 1