CoRegNet is an R/Bioconductor package to analyze large-scale transcriptomic data by highlighting sets of co-regulators. Based on a transcriptomic dataset, CoRegNet can be used to: reconstruct a ...large-scale co-regulatory network, integrate regulation evidences such as transcription factor binding sites and ChIP data, estimate sample-specific regulator activity, identify cooperative transcription factors and analyze the sample-specific combinations of active regulators through an interactive visualization tool. In this study CoRegNet was used to identify driver regulators of bladder cancer.
CoRegNet is available at http://bioconductor.org/packages/CoRegNet
remy.nicolle@issb.genopole.fr or mohamed.elati@issb.genopole.fr
Supplementary data are available at Bioinformatics online.
Relationships between genetic alterations, such as co-occurrence or mutual exclusivity, are often observed in cancer, where their understanding may provide new insights into etiology and clinical ...management. In this study, we combined statistical analyses and computational modeling to explain patterns of genetic alterations seen in 178 patients with bladder tumors (either muscle-invasive or non-muscle-invasive). A statistical analysis on frequently altered genes identified pair associations, including co-occurrence or mutual exclusivity. Focusing on genetic alterations of protein-coding genes involved in growth factor receptor signaling, cell cycle, and apoptosis entry, we complemented this analysis with a literature search to focus on nine pairs of genetic alterations of our dataset, with subsequent verification in three other datasets available publicly. To understand the reasons and contexts of these patterns of associations while accounting for the dynamics of associated signaling pathways, we built a logical model. This model was validated first on published mutant mice data, then used to study patterns and to draw conclusions on counter-intuitive observations, allowing one to formulate predictions about conditions where combining genetic alterations benefits tumorigenesis. For example, while CDKN2A homozygous deletions occur in a context of FGFR3-activating mutations, our model suggests that additional PIK3CA mutation or p21CIP deletion would greatly favor invasiveness. Furthermore, the model sheds light on the temporal orders of gene alterations, for example, showing how mutual exclusivity of FGFR3 and TP53 mutations is interpretable if FGFR3 is mutated first. Overall, our work shows how to predict combinations of the major gene alterations leading to invasiveness through two main progression pathways in bladder cancer.
Immune checkpoint inhibitors (ICIs) have proved to be an effective treatment for up to 40% of muscle-invasive bladder cancer (MIBC), but there is still a need for better performing biomarkers ...allowing to improve prediction of response to ICI. Response to immunotherapy in soft-tissue sarcoma, melanoma and renal cell carcinoma have been recently linked to the presence of tertiary lymphoid structures (TLS) in the tumour. TLS are organised aggregates of T, B and dendritic cells, participating in adaptive antitumor immune response. The chemokine CXCL13 is involved in the formation of TLS, and is reported as a reliable transcriptomic marker of TLS.
In this study, we sought to assess whether CXCL13 transcript expression can be a prognostic biomarker for ICI-treated MIBC patients and also investigated whether it can serve a biomarker of TLS in MIBC.
We analysed transcriptomic data from three publicly available MIBC cohorts and evaluated pathological slides from the TCGA-BLCA cohort for TLS presence and stage of maturation.
We showed that CXCL13 was independently associated with both prolonged survival (HR = 0.8, 95% CI 0.68–0.94) and objective response (p < 0.0001) in patients treated with ICI, at the difference of others immunological signatures. However, it was not a predictor for non–ICI-treated MIBC, suggesting a predictive effect of ICI efficacy. Finally, we validated that CXCL13 expression was correlated with tumour TLS in TCGA data set (p < 0.001), and can serve as a marker of TLS in bladder cancer.
These results support that CXCL13 expression, as a surrogate for tumour TLS, is a relevant candidate predictive biomarker of response to ICI for patients with advanced-stage bladder cancer.
Display omitted
•CXCL13 expression predicts survival for advanced-stage MIBC treated with ICI.•CXCL13 expression correlates with the presence of TLS in TCGA tumours.•Ba/Sq and stroma-rich subtypes have higher CXCL13 expression and harbour more TLS.•We suggest CXCL13 as a predictive biomarker for ICI response in bladder cancer.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
The Mitogen-Activated Protein Kinase (MAPK) network consists of tightly interconnected signalling pathways involved in diverse cellular processes, such as cell cycle, survival, apoptosis and ...differentiation. Although several studies reported the involvement of these signalling cascades in cancer deregulations, the precise mechanisms underlying their influence on the balance between cell proliferation and cell death (cell fate decision) in pathological circumstances remain elusive. Based on an extensive analysis of published data, we have built a comprehensive and generic reaction map for the MAPK signalling network, using CellDesigner software. In order to explore the MAPK responses to different stimuli and better understand their contributions to cell fate decision, we have considered the most crucial components and interactions and encoded them into a logical model, using the software GINsim. Our logical model analysis particularly focuses on urinary bladder cancer, where MAPK network deregulations have often been associated with specific phenotypes. To cope with the combinatorial explosion of the number of states, we have applied novel algorithms for model reduction and for the compression of state transition graphs, both implemented into the software GINsim. The results of systematic simulations for different signal combinations and network perturbations were found globally coherent with published data. In silico experiments further enabled us to delineate the roles of specific components, cross-talks and regulatory feedbacks in cell fate decision. Finally, tentative proliferative or anti-proliferative mechanisms can be connected with established bladder cancer deregulations, namely Epidermal Growth Factor Receptor (EGFR) over-expression and Fibroblast Growth Factor Receptor 3 (FGFR3) activating mutations.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Motivation: Genomic DNA regions are frequently lost or gained during tumor progression. Array Comparative Genomic Hybridization (array CGH) technology makes it possible to assess these changes in DNA ...in cancers, by comparison with a normal reference. The identification of systematically deleted or amplified genomic regions in a set of tumors enables biologists to identify genes involved in cancer progression because tumor suppressor genes are thought to be located in lost genomic regions and oncogenes, in gained regions. Array CGH profiles should also improve the classification of tumors. The achievement of these goals requires a methodology for detecting the breakpoints delimiting altered regions in genomic patterns and assigning a status (normal, gained or lost) to each chromosomal region. Results: We have developed a methodology for the automatic detection of breakpoints from array CGH profile, and the assignment of a status to each chromosomal region. The breakpoint detection step is based on the Adaptive Weights Smoothing (AWS) procedure and provides highly convincing results: our algorithm detects 97, 100 and 94% of breakpoints in simulated data, karyotyping results and manually analyzed profiles, respectively. The percentage of correctly assigned statuses ranges from 98.9 to 99.8% for simulated data and is 100% for karyotyping results. Our algorithm also outperforms other solutions on a public reference dataset. Availability: The R package GLAD (Gain and Loss Analysis of DNA) is available upon request
Molecular motors play important roles in force generation, migration, and intracellular trafficking. Changes in specific motor activities are altered in numerous diseases. KIF20A, a motor protein of ...the kinesin-6 family, is overexpressed in bladder cancer, and KIF20A levels correlate negatively with clinical outcomes. We report here a new role for the KIF20A kinesin motor protein in intracellular mechanics. Using optical tweezers to probe intracellular mechanics and surface AFM to probe cortical mechanics, we first confirm that bladder urothelial cells soften with an increasing cancer grade. We then show that inhibiting KIF20A makes the intracellular environment softer for both high- and low-grade bladder cancer cells. Upon inhibition of KIF20A, cortical stiffness also decreases in lower grade cells, while it surprisingly increases in higher grade malignant cells. Changes in cortical stiffness correlate with the interaction of KIF20A with myosin IIA. Moreover, KIF20A inhibition negatively regulates bladder cancer cell motility irrespective of the underlying substrate stiffness. Our results reveal a central role for a microtubule motor in cell mechanics and migration in the context of bladder cancer.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
In HPV-associated carcinomas, some examples of cancer-related genes altered by viral insertion and corresponding to potential therapeutic targets have been described, but no quantitative assessment ...of these events, including poorly recurrent targets, has been reported to date. To document these occurrences, we built and analyzed a database comprised of 1455 cases, including HPV genotypes and tumor localizations. Host DNA sequences targeted by viral integration were classified as "non-recurrent" (one single reported case; 838 loci), "weakly recurrent" (two reported cases; 82 loci), and highly recurrent (≥3 cases; 43 loci). Whereas the overall rate of cancer-related target genes was 3.3% in the Gencode database, this rate increased to 6.5% in "non-recurrent", 11.4% in "weakly recurrent", and 40.1% in "highly recurrent" genes targeted by integration (
= 4.9 × 10
). This rate was also significantly higher in tumors associated with high-risk HPV16/18/45 than other genotypes. Among the genes targeted by HPV insertion, 30.2% corresponded to direct or indirect druggable targets, a rate rising to 50% in "highly recurrent" targets. Using data from the literature and the DepMap 23Q4 release database, we found that genes targeted by viral insertion could be new candidates potentially involved in HPV-associated oncogenesis. A more systematic characterization of HPV/host fusion DNA sequences in HPV-associated cancers should provide a better knowledge of HPV-driven carcinogenesis and favor the development of personalize patient treatments.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
Purpose: To gain insight into genomic and trancriptomic subtypes of ductal carcinomas in situ of the breast (DCIS).
Experimental Design: We did a combined phenotypic and genomic analysis of a series ...of 57 DCIS integrated with gene expression profile analysis
for 26 of the 57 cases.
Results: Thirty-two DCIS exhibited a luminal phenotype; 21 were ERBB2 positive, and 4 were ERBB2/estrogen receptor (ER) negative with
1 harboring a bona fide basal-like phenotype. Based on a CGH analysis, genomic types were identified in this series of DCIS with the 1q gain/16q
loss combination observed in 3 luminal DCIS, the mixed amplifier pattern including all ERBB2, 12 luminal and 2 ERBB2 - /ER - DCIS, and the complex copy number alteration profile encompassing 14 luminal and 1 ERBB2 - /ER - DCIS. Eight cases (8 of 57; 14%) presented a TP53 mutation, all being amplifiers. Unsupervised analysis of gene expression profiles of 26 of the 57 DCIS showed that luminal
and ERBB2 -amplified, ER-negative cases clustered separately. We further investigated the effect of high and low copy number changes
on gene expression. Strikingly, amplicons but also low copy number changes especially on 1q, 8q, and 16q in DCIS regulated
the expression of a subset of genes in a very similar way to that recently described in invasive ductal carcinomas.
Conclusions: These combined approaches show that the molecular heterogeneity of breast ductal carcinomas exists already in in situ lesions and further indicate that DCIS and invasive ductal carcinomas share genomic alterations with a similar effect on
gene expression profile.
Muscle-invasive bladder cancer (MIBC) is an aggressive neoplasm with poor prognosis, lacking effective therapeutic targets. Oncogenic dependency on members of the TAM tyrosine kinase receptor family ...(TYRO3, AXL, MERTK) has been reported in several cancer types, but their role in bladder cancer has never been explored.
TAM receptor expression was evaluated in two series of human bladder tumours by gene expression (TCGA and CIT series), immunohistochemistry and western blotting analyses (CIT series). The role of the different TAM receptors was assessed by loss-of-function experiments and pharmaceutical inhibition in vitro and in vivo.
We reported a significantly higher expression of TYRO3, but not AXL or MERTK, in both non-MIBCs and MIBCs, compared to normal urothelium. Loss-of-function experiments identified a TYRO3-dependency of bladder carcinoma-derived cells both in vitro and in a mouse xenograft model, whereas AXL and MERTK depletion had only a minor impact on cell viability. Accordingly, TYRO3-dependent bladder tumour cells were sensitive to pharmacological treatment with two pan-TAM inhibitors. Finally, growth inhibition upon TYRO3 depletion relies on cell cycle inhibition and apoptosis associated with induction of tumour-suppressive signals.
Our results provide a preclinical proof of concept for TYRO3 as a potential therapeutic target in bladder cancer.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Background APOBEC-driven mutagenesis and functional positive selection of mutated genes may synergistically drive the higher frequency of some hotspot driver mutations compared to other mutations ...within the same gene, as we reported for FGFR3 S249C. Only a few APOBEC-associated driver hotspot mutations have been identified in bladder cancer (BCa). Here, we systematically looked for and characterised APOBEC-associated hotspots in BCa. Methods We analysed 602 published exome-sequenced BCas, for part of which gene expression data were also available. APOBEC-associated hotspots were identified by motif-mapping, mutation signature fitting and APOBEC-mediated mutagenesis comparison. Joint analysis of DNA hairpin stability and gene expression was performed to predict driver or passenger hotspots. Aryl hydrocarbon receptor (AhR) activity was calculated based on its target genes expression. Effects of AhR knockout/inhibition on BCa cell viability were analysed. Results We established a panel of 44 APOBEC-associated hotspot mutations in BCa, which accounted for about half of the hotspot mutations. Fourteen of them overlapped with the hotspots found in other cancer types with high APOBEC activity. They mostly occurred in the DNA lagging-strand templates and the loop of DNA hairpins. APOBEC-associated hotspots presented systematically a higher prevalence than the other mutations within each APOBEC-target gene, independently of their functional impact. A combined analysis of DNA loop stability and gene expression allowed to distinguish known passenger from known driver hotspot mutations in BCa, including loss-of-function mutations affecting tumour suppressor genes, and to predict new candidate drivers, such as AHR Q383H. We further characterised AHR Q383H as an activating driver mutation associated with high AhR activity in luminal tumours. High AhR activity was also found in tumours presenting amplifications of AHR and its co-receptor ARNT. We finally showed that BCa cells presenting those different genetic alterations were sensitive to AhR inhibition. Conclusions Our study identified novel potential drivers within APOBEC-associated hotspot mutations in BCa reinforcing the importance of APOBEC mutagenesis in BCa. It could allow a better understanding of BCa biology and aetiology and have clinical implications such as AhR as a potential therapeutic target. Our results also challenge the dogma that all hotspot mutations are drivers and mostly gain-of-function mutations affecting oncogenes. Keywords: Bladder cancer, APOBEC, Mutagenesis, Stem-loop, Driver mutation, Passenger mutation, Oncogene, Tumour suppressor gene, Aryl hydrocarbon receptor
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK