Adverse stem cell clones within a single patient's tumor predict clinical outcome in AML patients.

Acute myeloid leukemia (AML) patients suffer dismal prognosis upon treatment resistance. To study functional heterogeneity of resistance, we generated serially transplantable patient-derived xenograft (PDX) models from one patient with AML and twelve clones thereof, each derived from a single stem cell, as proven by genetic barcoding. Transcriptome and exome sequencing segregated clones according to their origin from relapse one or two. Undetectable for sequencing, multiplex fluorochrome-guided competitive in vivo treatment trials identified a subset of relapse two clones as uniquely resistant to cytarabine treatment. Transcriptional and proteomic profiles obtained from resistant PDX clones and refractory AML patients defined a 16-gene score that was predictive of clinical outcome in a large independent patient cohort. Thus, we identified novel genes related to cytarabine resistance and provide proof of concept that intra-tumor heterogeneity reflects inter-tumor heterogeneity in AML.

Systems-level effects of ectopic galectin-7 reconstitution in cervical cancer and its microenvironment.

BACKGROUND: Galectin-7 (Gal-7) is negatively regulated in cervical cancer, and appears to be a link between the apoptotic response triggered by cancer and the anti-tumoral activity of the immune system. Our understanding of how cervical cancer cells and their molecular networks adapt in response to the expression of Gal-7 remains limited. METHODS: Meta-analysis of Gal-7 expression was conducted in three cervical cancer cohort studies and TCGA. In silico prediction and bisulfite sequencing were performed to inquire epigenetic alterations. To study the effect of Gal-7 on cervical cancer, we ectopically re-expressed it in the HeLa and SiHa cervical cancer cell lines, and analyzed their transcriptome and SILAC-based proteome. We also examined the tumor and microenvironment host cell transcriptomes after xenotransplantation into immunocompromised mice. Differences between samples were assessed with the Kruskall-Wallis, Dunn's Multiple Comparison and T tests. Kaplan-Meier and log-rank tests were used to determine overall survival. RESULTS: Gal-7 was constantly downregulated in our meta-analysis (p < 0.0001). Tumors with combined high Gal-7 and low galectin-1 expression (p = 0.0001) presented significantly better prognoses (p = 0.005). In silico and bisulfite sequencing assays showed de novo methylation in the Gal-7 promoter and first intron. Cells re-expressing Gal-7 showed a high apoptosis ratio (p < 0.05) and their xenografts displayed strong growth retardation (p < 0.001). Multiple gene modules and transcriptional regulators were modulated in response to Gal-7 reconstitution, both in cervical cancer cells and their microenvironments (FDR < 0.05 %). Most of these genes and modules were associated with tissue morphogenesis, metabolism, transport, chemokine activity, and immune response. These functional modules could exert the same effects in vitro and in vivo, even despite different compositions between HeLa and SiHa samples. CONCLUSIONS: Gal-7 re-expression affects the regulation of molecular networks in cervical cancer that are involved in diverse cancer hallmarks, such as metabolism, growth control, invasion and evasion of apoptosis. The effect of Gal-7 extends to the microenvironment, where networks involved in its configuration and in immune surveillance are particularly affected.

Mechismo: predicting the mechanistic impact of mutations and modifications on molecular interactions.

Systematic interrogation of mutation or protein modification data is important to identify sites with functional consequences and to deduce global consequences from large data sets. Mechismo (mechismo.russellab.org) enables simultaneous consideration of thousands of 3D structures and biomolecular interactions to predict rapidly mechanistic consequences for mutations and modifications. As useful functional information often only comes from homologous proteins, we benchmarked the accuracy of predictions as a function of protein/structure sequence similarity, which permits the use of relatively weak sequence similarities with an appropriate confidence measure. For protein-protein, protein-nucleic acid and a subset of protein-chemical interactions, we also developed and benchmarked a measure of whether modifications are likely to enhance or diminish the interactions, which can assist the detection of modifications with specific effects. Analysis of high-throughput sequencing data shows that the approach can identify interesting differences between cancers, and application to proteomics data finds potential mechanistic insights for how post-translational modifications can alter biomolecular interactions.

Analysis and prediction of pathways in HeLa cells by integrating biological levels of organization with systems-biology approaches.

It has recently begun to be considered that cancer is a systemic disease and that it must be studied at every level of complexity using many of the currently available approaches, including high-throughput technologies and bioinformatics. To achieve such understanding in cervical cancer, we collected information on gene, protein and phosphoprotein expression of the HeLa cell line and performed a comprehensive analysis of the different signaling pathways, transcription networks and metabolic events in which they participate. A total expression analysis by RNA-Seq of the HeLa cell line showed that 19,974 genes were transcribed. Of these, 3,360 were over-expressed, and 2,129 under-expressed when compared to the NHEK cell line. A protein-protein interaction network was derived from the over-expressed genes and used to identify central elements and, together with the analysis of over-represented transcription factor motifs, to predict active signaling and regulatory pathways. This was further validated by Metal-Oxide Affinity Chromatography (MOAC) and Tandem Mass Spectrometry (MS/MS) assays which retrieved phosphorylated proteins. The 14-3-3 family members emerge as important regulators in carcinogenesis and as possible clinical targets. We observed that the different over- and under-regulated pathways in cervical cancer could be interrelated through elements that participate in crosstalks, therefore belong to what we term "meta-pathways". Additionally, we highlighted the relations of each one of the differentially represented pathways to one or more of the ten hallmarks of cancer. These features could be maintained in many other types of cancer, regardless of mutations or genomic rearrangements, and favor their robustness, adaptations and the evasion of tissue control. Probably, this could explain why cancer cells are not eliminated by selective pressure and why therapy trials directed against molecular targets are not as effective as expected.

Prokaryotic regulatory systems biology: Common principles governing the functional architectures of Bacillus subtilis and Escherichia coli unveiled by the natural decomposition approach.

Escherichia coli and Bacillus subtilis are two of the best-studied prokaryotic model organisms. Previous analyses of their transcriptional regulatory networks have shown that they exhibit high plasticity during evolution and suggested that both converge to scale-free-like structures. Nevertheless, beyond this suggestion, no analyses have been carried out to identify the common systems-level components and principles governing these organisms. Here we show that these two phylogenetically distant organisms follow a set of common novel biologically consistent systems principles revealed by the mathematically and biologically founded natural decomposition approach. The discovered common functional architecture is a diamond-shaped, matryoshka-like, three-layer (coordination, processing, and integration) hierarchy exhibiting feedback, which is shaped by four systems-level components: global transcription factors (global TFs), locally autonomous modules, basal machinery and intermodular genes. The first mathematical criterion to identify global TFs, the κ-value, was reassessed on B. subtilis and confirmed its high predictive power by identifying all the previously reported, plus three potential, master regulators and eight sigma factors. The functionally conserved cores of modules, basal cell machinery, and a set of non-orthologous common physiological global responses were identified via both orthologous genes and non-orthologous conserved functions. This study reveals novel common systems principles maintained between two phylogenetically distant organisms and provides a comparison of their lifestyle adaptations. Our results shed new light on the systems-level principles and the fundamental functions required by bacteria to sustain life.