Genome-wide association analyses of ovarian cancer patients undergoing primary debulking surgery identify candidate genes for residual disease.

Survival from ovarian cancer depends on the resection status after primary surgery. We performed genome-wide association analyses for resection status of 7705 ovarian cancer patients, including 4954 with high-grade serous carcinoma (HGSOC), to identify variants associated with residual disease. The most significant association with resection status was observed for rs72845444, upstream of MGMT, in HGSOC (p = 3.9 × 10-8). In gene-based analyses, PPP2R5C was the most strongly associated gene in HGSOC after stage adjustment. In an independent set of 378 ovarian tumours from the AGO-OVAR 11 study, variants near MGMT and PPP2R5C correlated with methylation and transcript levels, and PPP2R5C mRNA levels predicted progression-free survival in patients with residual disease. MGMT encodes a DNA repair enzyme, and PPP2R5C encodes the B56γ subunit of the PP2A tumour suppressor. Our results link heritable variation at these two loci with resection status in HGSOC.

Mitochondrial DNA mutations drive aerobic glycolysis to enhance checkpoint blockade response in melanoma.

The mitochondrial genome (mtDNA) encodes essential machinery for oxidative phosphorylation and metabolic homeostasis. Tumor mtDNA is among the most somatically mutated regions of the cancer genome, but whether these mutations impact tumor biology is debated. We engineered truncating mutations of the mtDNA-encoded complex I gene, Mt-Nd5, into several murine models of melanoma. These mutations promoted a Warburg-like metabolic shift that reshaped tumor microenvironments in both mice and humans, consistently eliciting an anti-tumor immune response characterized by loss of resident neutrophils. Tumors bearing mtDNA mutations were sensitized to checkpoint blockade in a neutrophil-dependent manner, with induction of redox imbalance being sufficient to induce this effect in mtDNA wild-type tumors. Patient lesions bearing >50% mtDNA mutation heteroplasmy demonstrated a response rate to checkpoint blockade that was improved by ~2.5-fold over mtDNA wild-type cancer. These data nominate mtDNA mutations as functional regulators of cancer metabolism and tumor biology, with potential for therapeutic exploitation and treatment stratification.

PBRM1 and KDM5C cooperate to define high-angiogenesis tumors and increased antiangiogenic response in renal cancer.

Vascular endothelial growth factor receptor tyrosine kinase inhibitors (VEGFR-TKIs) are key antiangiogenic drugs for renal cancer treatment. While Von Hippel-Lindau dysfunction constitutes the base for VEGFR-TKIs sensitivity, the role for individual and concurrent mutations in the genes encoding for the chromatin remodelers Polybromo-1 (PBRM1) and Lysine Demethylase 5C (KDM5C) is poorly understood. Here, we analyzed the tumor mutational and expression profiles of 155 unselected clear cell RCC (ccRCC) cases treated with first-line VEGFR-TKIs and the ccRCC cases of IMmotion151 trial were used for validation. We found that concurrent PBRM1 and KDM5C (PBRM1&KDM5C) mutations occurred in 4-9% of cases and were enriched in Memorial Sloan Kettering Cancer Center favorable-risk patients. In our cohort, tumors only mutated in PBRM1 or concurrently mutated in PBRM1 and KDM5C had increased angiogenesis (P=0.0068 and 0.039; respectively), and tumors only mutated in KDM5C showed a similar trend. Best response to VEGFR-TKIs corresponded to PBRM1&KDM5C mutated cases, followed by those mutated only in KDM5C or only in PBRM1 (P=0.050, 0.040 and 0.027 versus non-mutated cases, respectively), with a trend for longer progression free survival (PFS) in the group with only PBRM1 mutated (HR=0.64; P=0.059). Validation in the IMmotion151 trial revealed a similar correlation with increased angiogenesis and the PFS of patients in the VEGFR-TKI-arm was the longest in PBRM1&KDM5C mutated cases, intermediate for only PBRM1 or only KDM5C mutated patients and the shortest in non-mutated cases (P=0.009 and 0.025, for PBRM1&KDM5C and PBRM1 versus non-mutated cases). In conclusion, somatic PBRM1 and KDM5C mutations are common in patients with metastatic ccRCC and likely cooperate increasing tumor angiogenesis and VEGFR-TKI-based antiangiogenic therapy benefit.

Deubiquitinase USP9X loss sensitizes renal cancer cells to mTOR inhibition.

Mammalian target of rapamycin (mTOR) is a central regulator of mammalian metabolism and physiology. Aberrant hyperactivation of the mTOR pathway promotes tumor growth and metastasis, and can also promote tumor resistance to chemotherapy and cancer drugs; this makes mTOR an attractive cancer therapeutic target. mTOR inhibitors have been approved to treat cancer; however, the mechanisms underlying drug sensitivity remain poorly understood. Here, whole exome sequencing of three chromophobe renal cell carcinoma (chRCC) patients with exceptional mTOR inhibitor sensitivity revealed that all three patients shared somatic mutations in the deubiquitinase gene USP9X. The clonal characteristics of the mutations, which were amassed by studying multiple patients' primary and metastatic samples from various years, together with the low USP9X mutation rate in unselected chRCC series, reinforced a causal link between USP9X and mTOR inhibitor sensitivity. Rapamycin treatment of USP9X-depleted HeLa and renal cancer 786-O cells, along with the pharmacological inhibition of USP9X, confirmed that this protein plays a role in patients' sensitivity to mTOR inhibitors. USP9X was not found to exert a direct effect on mTORC1, but subsequent ubiquitylome analyses identified p62 as a direct USP9X target. Increased p62 ubiquitination and the augmented rapamycin effect upon bortezomib treatment, together with the results of p62 and LC3 immunofluorescence assays, suggested that dysregulated autophagy in USP9X-depleted cells can have a synergistic effect with mTOR inhibitors. In summary, we show that USP9X constitutes a potential novel marker of sensitivity to mTOR inhibitors in chRCC patients, and represents a clinical strategy for increasing the sensitivity to these drugs.

Clinical parameters affecting survival outcomes in patients with low-grade serous ovarian carcinoma: an international multicentre analysis.

BACKGROUND: Women with low-grade ovarian serous carcinoma (LGSC) benefit from surgical treatment; however, the role of chemotherapy is controversial. We examined an international database through the Ovarian Cancer Association Consortium to identify factors that affect survival in LGSC. METHODS: We performed a retrospective cohort analysis of patients with LGSC who had had primary surgery and had overall survival data available. We performed univariate and multivariate analyses of progression-free survival and overall survival, and generated Kaplan-Meier survival curves. RESULTS: Of the 707 patients with LGSC, 680 (96.2%) had available overall survival data. The patients' median age overall was 54 years. Of the 659 patients with International Federation of Obstetrics and Gynecology stage data, 156 (23.7%) had stage I disease, 64 (9.7%) had stage II, 395 (59.9%) had stage III, and 44 (6.7%) had stage IV. Of the 377 patients with surgical data, 200 (53.0%) had no visible residual disease. Of the 361 patients with chemotherapy data, 330 (91.4%) received first-line platinum-based chemotherapy. The median follow-up duration was 5.0 years. The median progression-free survival and overall survival were 43.2 months and 110.4 months, respectively. Multivariate analysis indicated a statistically significant impact of stage and residual disease on progression-free survival and overall survival. Platinum-based chemotherapy was not associated with a survival advantage. CONCLUSION: This multicentre analysis indicates that complete surgical cytoreduction to no visible residual disease has the most impact on improved survival in LGSC. This finding could immediately inform and change practice.

Tumour mitochondrial DNA mutations drive aerobic glycolysis to enhance checkpoint blockade.

The mitochondrial genome encodes essential machinery for respiration and metabolic homeostasis but is paradoxically among the most common targets of somatic mutation in the cancer genome, with truncating mutations in respiratory complex I genes being most over-represented1. While mitochondrial DNA (mtDNA) mutations have been associated with both improved and worsened prognoses in several tumour lineages1-3, whether these mutations are drivers or exert any functional effect on tumour biology remains controversial. Here we discovered that complex I-encoding mtDNA mutations are sufficient to remodel the tumour immune landscape and therapeutic resistance to immune checkpoint blockade. Using mtDNA base editing technology4 we engineered recurrent truncating mutations in the mtDNA-encoded complex I gene, Mt-Nd5, into murine models of melanoma. Mechanistically, these mutations promoted utilisation of pyruvate as a terminal electron acceptor and increased glycolytic flux without major effects on oxygen consumption, driven by an over-reduced NAD pool and NADH shuttling between GAPDH and MDH1, mediating a Warburg-like metabolic shift. In turn, without modifying tumour growth, this altered cancer cell-intrinsic metabolism reshaped the tumour microenvironment in both mice and humans, promoting an anti-tumour immune response characterised by loss of resident neutrophils. This subsequently sensitised tumours bearing high mtDNA mutant heteroplasmy to immune checkpoint blockade, with phenocopy of key metabolic changes being sufficient to mediate this effect. Strikingly, patient lesions bearing >50% mtDNA mutation heteroplasmy also demonstrated a >2.5-fold improved response rate to checkpoint inhibitor blockade. Taken together these data nominate mtDNA mutations as functional regulators of cancer metabolism and tumour biology, with potential for therapeutic exploitation and treatment stratification.

Genomic and immune landscape Of metastatic pheochromocytoma and paraganglioma.

The mechanisms triggering metastasis in pheochromocytoma/paraganglioma are unknown, hindering therapeutic options for patients with metastatic tumors (mPPGL). Herein we show by genomic profiling of a large cohort of mPPGLs that high mutational load, microsatellite instability and somatic copy-number alteration burden are associated with ATRX/TERT alterations and are suitable prognostic markers. Transcriptomic analysis defines the signaling networks involved in the acquisition of metastatic competence and establishes a gene signature related to mPPGLs, highlighting CDK1 as an additional mPPGL marker. Immunogenomics accompanied by immunohistochemistry identifies a heterogeneous ecosystem at the tumor microenvironment level, linked to the genomic subtype and tumor behavior. Specifically, we define a general immunosuppressive microenvironment in mPPGLs, the exception being PD-L1 expressing MAML3-related tumors. Our study reveals canonical markers for risk of metastasis, and suggests the usefulness of including immune parameters in clinical management for PPGL prognostication and identification of patients who might benefit from immunotherapy.

Comprehensive molecular analysis of immortalization hallmarks in thyroid cancer reveals new prognostic markers.

BACKGROUND: Comprehensive molecular studies on tumours are needed to delineate immortalization process steps and identify sensitive prognostic biomarkers in thyroid cancer. METHODS AND RESULTS: In this study, we extensively characterize telomere-related alterations in a series of 106 thyroid tumours with heterogeneous clinical outcomes. Using a custom-designed RNA-seq panel, we identified five telomerase holoenzyme-complex genes upregulated in clinically aggressive tumours compared to tumours from long-term disease-free patients, being TERT and TERC denoted as independent prognostic markers by multivariate regression model analysis. Characterization of alterations related to TERT re-expression revealed that promoter mutations, methylation and/or copy gains exclusively co-occurred in clinically aggressive tumours. Quantitative-FISH (fluorescence in situ hybridization) analysis of telomere lengths showed a significant shortening in these carcinomas, which matched with a high proliferative rate measured by Ki-67 immunohistochemistry. RNA-seq data analysis indicated that short-telomere tumours exhibit an increased transcriptional activity in the 5-Mb-subtelomeric regions, site of several telomerase-complex genes. Gene upregulation enrichment was significant for specific chromosome-ends such as the 5p, where TERT is located. Co-FISH analysis of 5p-end and TERT loci showed a more relaxed chromatin configuration in short telomere-length tumours compared to normal telomere-length tumours. CONCLUSIONS: Overall, our findings support that telomere shortening leads to a 5p subtelomeric region reorganization, facilitating the transcription and accumulation of alterations at TERT-locus.

Cabozantinib combination therapy for the treatment of solid tumors: a systematic review.

Background: Cabozantinib monotherapy is approved for the treatment of several types of solid tumors. Investigation into the use of cabozantinib combined with other therapies is increasing. To understand the evidence in this area, we performed a systematic review of cabozantinib combination therapy for the treatment of solid tumors in adults. Methods: This study was designed in accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses, and the protocol was registered with PROSPERO (CRD42020144680). On 9 October 2020, we searched for clinical trials and observational studies of cabozantinib as part of a combination therapy for solid tumors using Embase, MEDLINE, and Cochrane databases, and by screening relevant congress abstracts. Eligible studies reported clinical or safety outcomes, or biomarker data. Randomized and observational studies with a sample size of fewer than 25 and studies of cabozantinib monotherapy were excluded. For each study, quality was assessed using National Institute for Health and Care Excellence methodology, and the study characteristics were described qualitatively. This study was funded by Ipsen. Results: Of 2421 citations identified, 32 articles were included (6 with results from randomized studies, 24 with results from non-randomized phase I or II studies, and 2 with results from both). The most commonly studied tumor types were metastatic urothelial carcinoma/genitourinary tumors and castration-resistant prostate cancer (CRPC). Findings from randomized studies suggested that cabozantinib combined with other therapies may lead to better progression-free survival than some current standards of care in renal cell carcinoma, CRPC, and non-small-cell lung cancer. The most common adverse events were hypertension, diarrhea, and fatigue. Conclusion: This review demonstrates the promising efficacy outcomes of cabozantinib combined with other therapies, and a safety profile similar to cabozantinib alone. However, the findings are limited by the fact that most of the identified studies were reported as congress abstracts only. More evidence from randomized trials is needed to explore cabozantinib as a combination therapy further.

Cabozantinib for the treatment of solid tumors: a systematic review.

Background: Cabozantinib is approved, in various settings, for the treatment of renal cell carcinoma, medullary thyroid cancer, and hepatocellular carcinoma, and it has been investigated for the treatment of other cancers. With the available evidence and the real-world performance of cabozantinib compared with clinical trial data, we performed a systematic review of cabozantinib monotherapy as treatment for solid tumors in adults. Methods: This study was designed in accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses and registered with PROSPERO (CRD42020144680). We searched for clinical and observational studies of cabozantinib monotherapy for solid tumors using Embase, MEDLINE, and Cochrane databases (October 2020), and screened relevant congress abstracts. Eligible studies reported clinical or safety outcomes, or biomarker data. Small studies (n < 25) and studies of cabozantinib combination therapies were excluded. Quality was assessed using National Institute for Health and Care Excellence methodology, and study characteristics were described qualitatively. Results: Of 2888 citations, 114 were included (52 randomized studies, 29 observational studies, 32 nonrandomized phase I or II studies or pilot trials, and 1 analysis of data from a randomized study and a nonrandomized study). Beyond approved indications, other tumors studied were castration-resistant prostate cancer, urothelial carcinoma, Ewing sarcoma, osteosarcoma, uveal melanoma, non-small-cell lung cancer, Merkel cell carcinoma, glioblastoma, pheochromocytomas and paragangliomas, cholangiocarcinoma, gastrointestinal stromal tumor, colorectal cancer, salivary gland cancer, carcinoid and pancreatic neuroendocrine tumors, and breast, endometrial and ovarian cancers. The most common adverse events were hypertension, diarrhea, and fatigue. Conclusion: The identified evidence demonstrates the positive efficacy/effectiveness of cabozantinib monotherapy in various solid tumor types, with safety findings being consistent with those observed with other VEGFR-targeting tyrosine kinase inhibitors. When available, real-world findings were consistent with the data reported from clinical trials. A limitation of this review is the high proportion of abstracts; however, this allowed us to capture the most up-to-date findings.

Single-nucleotide polymorphism associations with efficacy and toxicity in metastatic castration-resistant prostate cancer treated with cabazitaxel.

Background: The aim of this study was to evaluate the impact of certain single-nucleotide polymorphisms (SNPs) in cabazitaxel activity and toxicity in patients with metastatic castration-resistant prostate cancer (mCRPC). Patients & methods: 56 SNPs in five genes (CYP3A4, CYP3A5, ABCB1, TUBB1 and CYP2C8) were genotyped in 67 mCRPC patients and their correlation with outcomes analyzed. Results:TUBB1-rs151352 (hazard ratio: 0.52) and CYP2C8-rs1341164 (hazard ratio: 0.53) were associated with better overall survival, and CYP2C8-rs1058932 with biochemical progression (odds ratio: 6.60) in multivariate analysis. ABCB1-rs17327624 correlated with severe toxicity ≥grade 3 (odds ratio: 8.56) and CYP2C8-rs11572093 with asthenia (odds ratio: 8.12). Conclusion: Genetic variants in mCRPC patients could explain different outcomes with cabazitaxel. Nonetheless, the small sample size and the high number of SNPs analyzed mean that the results are only hypothesis-generating and require further validation.

Genome-Wide Meta-Analysis Identifies Variants in DSCAM and PDLIM3 That Correlate with Efficacy Outcomes in Metastatic Renal Cell Carcinoma Patients Treated with Sunitinib.

Individual response to sunitinib in metastatic renal cell carcinoma (mRCC) patients is highly variable. Earlier, sunitinib outcome was related to single nucleotide polymorphisms (SNPs) in CYP3A5 and ABCB1. Our aim is to provide novel insights into biological mechanisms underlying sunitinib action. We included mRCC patients from the European EuroTARGET consortium (n = 550) and the RIKEN cohort in Japan (n = 204) which were analysed separately and in a meta-analysis of genome-wide association studies (GWAS). SNPs were tested for association with progression-free survival (PFS) and overall survival (OS) using Cox regression. Summary statistics were combined using a fixed effect meta-analysis. SNP rs28520013 in PDLIM3 and the correlated SNPs rs2205096 and rs111356738 both in DSCAM, showed genome-wide significance (p < 5 × 10−8) with PFS and OS in the meta-analysis. The variant T-allele of rs28520013 associated with an inferior PFS of 5.1 months compared to 12.5 months in non-carriers (p = 4.02 × 10−10, HR = 7.26). T-allele carriers of rs28520013 showed an inferior OS of 6.9 months versus 30.2 months in non-carriers (p = 1.62 × 10−8, HR = 5.96). In this GWAS we identified novel genetic variants in PDLIM3 and DSCAM that impact PFS and OS in mRCC patients receiving sunitinib. The underlying link between the identified genes and the molecular mechanisms of sunitinib action needs to be elucidated.

Considerations on diagnosis and surveillance measures of PTEN hamartoma tumor syndrome: clinical and genetic study in a series of Spanish patients.

BACKGROUND: The limited knowledge about the PTEN hamartoma tumor syndrome (PHTS) makes its diagnosis a challenging task. We aimed to define the clinical and genetic characteristics of this syndrome in the Spanish population and to identify new genes potentially associated with the disease. RESULTS: We reviewed the clinical data collected through a specific questionnaire in a series of 145 Spanish patients with a phenotypic features compatible with PHTS and performed molecular characterization through several approaches including next generation sequencing and whole exome sequencing (WES). Macrocephaly, mucocutaneous lesions, gastrointestinal polyposis and obesity are prevalent phenotypic features in PHTS and help predict the presence of a PTEN germline variant in our population. We also find that PHTS patients are at risk to develop cancer in childhood or adolescence. Furthermore, we observe a high frequency of variants in exon 1 of PTEN, which are associated with renal cancer and overexpression of KLLN and PTEN. Moreover, WES revealed variants in genes like NEDD4 that merit further research. CONCLUSIONS: This study expands previously reported findings in other PHTS population studies and makes new contributions regarding clinical and molecular aspects of PHTS, which are useful for translation to the clinic and for new research lines.

Clinical pharmacogenetic analysis in 5,001 individuals with diagnostic Exome Sequencing data.

Exome sequencing is utilized in routine clinical genetic diagnosis. The technical robustness of repurposing large-scale next-generation sequencing data for pharmacogenetics has been demonstrated, supporting the implementation of preemptive pharmacogenetic strategies based on adding clinical pharmacogenetic interpretation to exomes. However, a comprehensive study analyzing all actionable pharmacogenetic alleles contained in international guidelines and applied to diagnostic exome data has not been performed. Here, we carried out a systematic analysis based on 5001 Spanish or Latin American individuals with diagnostic exome data, either Whole Exome Sequencing (80%), or the so-called Clinical Exome Sequencing (20%) (60 Mb and 17 Mb, respectively), to provide with global and gene-specific clinical pharmacogenetic utility data. 788 pharmacogenetic alleles, distributed through 19 genes included in Clinical Pharmacogenetics Implementation Consortium guidelines were analyzed. We established that Whole Exome and Clinical Exome Sequencing performed similarly, and 280 alleles in 11 genes (CACNA1S, CYP2B6, CYP2C9, CYP4F2, DPYD, G6PD, NUDT15, RYR1, SLCO1B1, TPMT, and UGT1A1) could be used to inform of pharmacogenetic phenotypes that change drug prescription. Each individual carried in average 2.2 alleles and overall 95% (n = 4646) of the cohort could be informed of at least one actionable pharmacogenetic phenotype. Differences in variant allele frequency were observed among the populations studied and the corresponding gnomAD population for 7.9% of the variants. In addition, in the 11 selected genes we uncovered 197 novel variants, among which 27 were loss-of-function. In conclusion, we provide with the landscape of actionable pharmacogenetic information contained in diagnostic exomes, that can be used preemptively in the clinics.

PharmVar GeneFocus: CYP3A5.

The Pharmacogene Variation Consortium (PharmVar) catalogs star (*) allele nomenclature for the polymorphic human CYP3A5 gene. Genetic variation within the CYP3A5 gene locus impacts the metabolism of several clinically important drugs, including the immunosuppressants tacrolimus, sirolimus, cyclosporine, and the benzodiazepine midazolam. Variable CYP3A5 activity is of clinical importance regarding tacrolimus metabolism. This GeneFocus provides a CYP3A5 gene summary with a focus on aspects regarding standardized nomenclature. In addition, this review also summarizes recent changes and updates, including the retirement of several allelic variants and provides an overview of how PharmVar CYP3A5 star allele nomenclature is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Clinical Pharmacogenetics Implementation Consortium Guideline for the Use of Aminoglycosides Based on MT-RNR1 Genotype.

Aminoglycosides are widely used antibiotics with notable side effects, such as nephrotoxicity, vestibulotoxicity, and sensorineural hearing loss (cochleotoxicity). MT-RNR1 is a gene that encodes the 12s rRNA subunit and is the mitochondrial homologue of the prokaryotic 16s rRNA. Some MT-RNR1 variants (i.e., m.1095T>C; m.1494C>T; m.1555A>G) more closely resemble the bacterial 16s rRNA subunit and result in increased risk of aminoglycoside-induced hearing loss. Use of aminoglycosides should be avoided in individuals with an MT-RNR1 variant associated with an increased risk of aminoglycoside-induced hearing loss unless the high risk of permanent hearing loss is outweighed by the severity of infection and safe or effective alternative therapies are not available. We summarize evidence from the literature supporting this association and provide therapeutic recommendations for the use of aminoglycosides based on MT-RNR1 genotype (updates at https://cpicpgx.org/guidelines/ and www.pharmgkb.org).

Analysis of Telomere Maintenance Related Genes Reveals NOP10 as a New Metastatic-Risk Marker in Pheochromocytoma/Paraganglioma.

One of the main problems we face with PPGL is the lack of molecular markers capable of predicting the development of metastases in patients. Telomere-related genes, such as TERT and ATRX, have been recently described in PPGL, supporting the association between the activation of immortalization mechanisms and disease progression. However, the contribution of other genes involving telomere preservation machinery has not been previously investigated. In this work, we aimed to analyze the prognostic value of a comprehensive set of genes involved in telomere maintenance. For this study, we collected 165 PPGL samples (97 non-metastatic/63 metastatic), genetically characterized, in which the expression of 29 genes of interest was studied by NGS. Three of the 29 genes studied, TERT, ATRX and NOP10, showed differential expression between metastatic and non-metastatic cases, and alterations in these genes were associated with a shorter time to progression, independent of SDHB-status. We studied telomere length by Q-FISH in patient samples and in an in vitro model. NOP10 overexpressing tumors displayed an intermediate-length telomere phenotype without ALT, and in vitro results suggest that NOP10 has a role in telomerase-dependent telomere maintenance. We also propose the implementation of NOP10 IHC to better stratify PPGL patients.

Circulating Levels of the Interferon-γ-Regulated Chemokines CXCL10/CXCL11, IL-6 and HGF Predict Outcome in Metastatic Renal Cell Carcinoma Patients Treated with Antiangiogenic Therapy.

Sunitinib and pazopanib are standard first-line treatments for patients with metastatic renal cell carcinoma (mRCC). Nonetheless, as the number of treatment options increases, there is a need to identify biomarkers that can predict drug efficacy and toxicity. In this prospective study we evaluated a set of biomarkers that had been previously identified within a secretory signature in mRCC patients. This set includes tumor expression of c-Met and serum levels of HGF, IL-6, IL-8, CXCL9, CXCL10 and CXCL11. Our cohort included 60 patients with mRCC from 10 different Spanish hospitals who received sunitinib (n = 51), pazopanib (n = 4) or both (n = 5). Levels of biomarkers were studied in relation to response rate, progression-free survival (PFS) and overall survival (OS). High tumor expression of c-Met and high basal serum levels of HGF, IL-6, CXCL11 and CXCL10 were significantly associated with reduced PFS and/or OS. In multivariable Cox regression analysis, CXCL11 was identified as an independent biomarker predictive of shorter PFS and OS, and HGF was an independent predictor of reduced PFS. Correlation analyses using our cohort of patients and patients from TCGA showed that HGF levels were significantly correlated with those of IL-6, CXCL11 and CXCL10. Bioinformatic protein-protein network analysis revealed a significant interaction between these proteins, all this suggesting a coordinated expression and secretion. We also developed a prognostic index that considers this group of biomarkers, where high values in mRCC patients can predict higher risk of relapse (HR 5.28 [2.32-12.0], p < 0.0001). In conclusion, high plasma HGF, CXCL11, CXCL10 and IL-6 levels are associated with worse outcome in mRCC patients treated with sunitinib or pazopanib. Our findings also suggest that these factors may constitute a secretory cluster that acts coordinately to promote tumor growth and resistance to antiangiogenic therapy.