Facing the diversity of omics data and the difficulty of selecting one result over all those produced by several methods, consensus strategies have the potential to reconcile multiple inputs and to ...produce robust results. Here, we introduce ClustOmics, a generic consensus clustering tool that we use in the context of cancer subtyping. ClustOmics relies on a non-relational graph database, which allows for the simultaneous integration of both multiple omics data and results from various clustering methods. This new tool conciliates input clusterings, regardless of their origin, their number, their size or their shape. ClustOmics implements an intuitive and flexible strategy, based upon the idea of evidence accumulation clustering. ClustOmics computes co-occurrences of pairs of samples in input clusters and uses this score as a similarity measure to reorganize data into consensus clusters. We applied ClustOmics to multi-omics disease subtyping on real TCGA cancer data from ten different cancer types. We showed that ClustOmics is robust to heterogeneous qualities of input partitions, smoothing and reconciling preliminary predictions into high-quality consensus clusters, both from a computational and a biological point of view. The comparison to a state-of-the-art consensus-based integration tool, COCA, further corroborated this statement. However, the main interest of ClustOmics is not to compete with other tools, but rather to make profit from their various predictions when no gold-standard metric is available to assess their significance.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We undertook a systematic study focused on the matricellular protein Thrombospondin-1 (THBS1) to uncover molecular mechanisms underlying the role of THBS1 in glioblastoma (GBM) development. THBS1 was ...found to be increased with glioma grades. Mechanistically, we show that the TGFβ canonical pathway transcriptionally regulates THBS1, through SMAD3 binding to the THBS1 gene promoter. THBS1 silencing inhibits tumour cell invasion and growth, alone and in combination with anti-angiogenic therapy. Specific inhibition of the THBS1/CD47 interaction using an antagonist peptide decreases cell invasion. This is confirmed by CD47 knock-down experiments. RNA sequencing of patient-derived xenograft tissue from laser capture micro-dissected peripheral and central tumour areas demonstrates that THBS1 is one of the gene with the highest connectivity at the tumour borders. All in all, these data show that TGFβ1 induces THBS1 expression via Smad3 which contributes to the invasive behaviour during GBM expansion. Furthermore, tumour cell-bound CD47 is implicated in this process.
The structure of messenger RNA is important for post-transcriptional regulation, mainly because it affects binding of trans-acting factors. However, little is known about the in vivo structure of ...full-length mRNAs. Here we present hiCLIP, a biochemical technique for transcriptome-wide identification of RNA secondary structures interacting with RNA-binding proteins (RBPs). Using this technique to investigate RNA structures bound by Staufen 1 (STAU1) in human cells, we uncover a dominance of intra-molecular RNA duplexes, a depletion of duplexes from coding regions of highly translated mRNAs, an unexpected prevalence of long-range duplexes in 3' untranslated regions (UTRs), and a decreased incidence of single nucleotide polymorphisms in duplex-forming regions. We also discover a duplex spanning 858 nucleotides in the 3' UTR of the X-box binding protein 1 (XBP1) mRNA that regulates its cytoplasmic splicing and stability. Our study reveals the fundamental role of mRNA secondary structures in gene expression and introduces hiCLIP as a widely applicable method for discovering new, especially long-range, RNA duplexes.
Full text
Available for:
DOBA, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, SBMB, SIK, UILJ, UKNU, UL, UM, UPUK
Hi-C is one of the main methods for investigating spatial co-localisation of DNA in the nucleus. However, the raw sequencing data obtained from Hi-C experiments suffer from large biases and spurious ...contacts, making it difficult to identify true interactions. Existing methods use complex models to account for biases and do not provide a significance threshold for detecting interactions. Here we introduce a simple binomial probabilistic model that resolves complex biases and distinguishes between true and false interactions. The model corrects biases of known and unknown origin and yields a p-value for each interaction, providing a reliable threshold based on significance. We demonstrate this experimentally by testing the method against a random ligation dataset. Our method outperforms previous methods and provides a statistical framework for further data analysis, such as comparisons of Hi-C interactions between different conditions. GOTHiC is available as a BioConductor package (http://www.bioconductor.org/packages/release/bioc/html/GOTHiC.html).
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract
G-quadruplex ligands exert their antiproliferative effects through telomere-dependent and telomere-independent mechanisms, but the inter-relationships among autophagy, cell growth arrest and ...cell death induced by these ligands remain largely unexplored. Here, we demonstrate that the G-quadruplex ligand 20A causes growth arrest of cancer cells in culture and in a HeLa cell xenografted mouse model. This response is associated with the induction of senescence and apoptosis. Transcriptomic analysis of 20A treated cells reveals a significant functional enrichment of biological pathways related to growth arrest, DNA damage response and the lysosomal pathway. 20A elicits global DNA damage but not telomeric damage and activates the ATM and autophagy pathways. Loss of ATM following 20A treatment inhibits both autophagy and senescence and sensitizes cells to death. Moreover, disruption of autophagy by deletion of two essential autophagy genes ATG5 and ATG7 leads to failure of CHK1 activation by 20A and subsequently increased cell death. Our results, therefore, identify the activation of ATM by 20A as a critical player in the balance between senescence and apoptosis and autophagy as one of the key mediators of such regulation. Thus, targeting the ATM/autophagy pathway might be a promising strategy to achieve the maximal anticancer effect of this compound.
The establishment of the embryonic and trophoblast lineages is a developmental decision underpinned by dramatic differences in the epigenetic landscape of the two compartments. However, it remains ...unknown how epigenetic information and transcription factor networks map to the 3D arrangement of the genome, which in turn may mediate transcriptional divergence between the two cell lineages. Here, we perform promoter capture Hi-C experiments in mouse trophoblast (TSC) and embryonic (ESC) stem cells to understand how chromatin conformation relates to cell-specific transcriptional programmes. We find that key TSC genes that are kept repressed in ESCs exhibit interactions between H3K27me3-marked regions in ESCs that depend on Polycomb repressive complex 1. Interactions that are prominent in TSCs are enriched for enhancer-gene contacts involving key TSC transcription factors, as well as TET1, which helps to maintain the expression of TSC-relevant genes. Our work shows that the first developmental cell fate decision results in distinct chromatin conformation patterns establishing lineage-specific contexts involving both repressive and active interactions.
This protocol explains how to use the online integrated pipeline 'peak-motifs' (http://rsat.ulb.ac.be/rsat/) to predict motifs and binding sites in full-size peak sets obtained by chromatin ...immunoprecipitation-sequencing (ChIP-seq) or related technologies. The workflow combines four time- and memory-efficient motif discovery algorithms to extract significant motifs from the sequences. Discovered motifs are compared with databases of known motifs to identify potentially bound transcription factors. Sequences are scanned to predict transcription factor binding sites and analyze their enrichment and positional distribution relative to peak centers. Peaks and binding sites are exported as BED tracks that can be uploaded into the University of California Santa Cruz (UCSC) genome browser for visualization in the genomic context. This protocol is illustrated with the analysis of a set of 6,000 peaks (8 Mb in total) bound by the Drosophila transcription factor Krüppel. The complete workflow is achieved in about 25 min of computational time on the Regulatory Sequence Analysis Tools (RSAT) Web server. This protocol can be followed in about 1 h.
The discovery of functionally relevant KRAS effectors in lung and pancreatic ductal adenocarcinoma (LUAD and PDAC) may yield novel molecular targets or mechanisms amenable to inhibition strategies. ...Phospholipids availability has been appreciated as a mechanism to modulate KRAS oncogenic potential. Thus, phospholipid transporters may play a functional role in KRAS-driven oncogenesis. Here, we identified and systematically studied the phospholipid transporter PITPNC1 and its controlled network in LUAD and PDAC.
Genetic modulation of KRAS expression as well as pharmacological inhibition of canonical effectors was completed. PITPNC1 genetic depletion was performed in in vitro and in vivo LUAD and PDAC models. PITPNC1-deficient cells were RNA sequenced, and Gene Ontology and enrichment analyses were applied to the output data. Protein-based biochemical and subcellular localization assays were run to investigate PITPNC1-regulated pathways. A drug repurposing approach was used to predict surrogate PITPNC1 inhibitors that were tested in combination with KRASG12C inhibitors in 2D, 3D, and in vivo models.
PITPNC1 was increased in human LUAD and PDAC, and associated with poor patients' survival. PITPNC1 was regulated by KRAS through MEK1/2 and JNK1/2. Functional experiments showed PITPNC1 requirement for cell proliferation, cell cycle progression and tumour growth. Furthermore, PITPNC1 overexpression enhanced lung colonization and liver metastasis. PITPNC1 regulated a transcriptional signature which highly overlapped with that of KRAS, and controlled mTOR localization via enhanced MYC protein stability to prevent autophagy. JAK2 inhibitors were predicted as putative PITPNC1 inhibitors with antiproliferative effect and their combination with KRASG12C inhibitors elicited a substantial anti-tumour effect in LUAD and PDAC.
Our data highlight the functional and clinical relevance of PITPNC1 in LUAD and PDAC. Moreover, PITPNC1 constitutes a new mechanism linking KRAS to MYC, and controls a druggable transcriptional network for combinatorial treatments.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Lysosomes play a key role in regulating cell death in response to cancer therapies, yet little is known on the possible role of lysosomes in the therapeutic efficacy of G-quadruplex DNA ligands (G4L) ...in cancer cells. Here, we investigate the relationship between the modulation of lysosomal membrane damage and the degree to which cancer cells respond to the cytotoxic effects of G-quadruplex ligands belonging to the triarylpyridine family. Our results reveal that the lead compound of this family, 20A promotes the enlargement of the lysosome compartment as well as the induction of lysosome-relevant mRNAs. Interestingly, the combination of 20A and chloroquine (an inhibitor of lysosomal functions) led to a significant induction of lysosomal membrane permeabilization coupled to massive cell death. Similar effects were observed when chloroquine was added to three new triarylpyridine derivatives. Our findings thus uncover the lysosomal effects of triarylpyridines compounds and delineate a rationale for combining these compounds with chloroquine to increase their anticancer effects.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK