Addiction to Runx1 is partially attenuated by loss of p53 in the Eµ-Myc lymphoma model

The Runx genes function as dominant oncogenes that collaborate potently with Myc or loss of p53 to induce lymphoma when over-expressed. Here we examined the requirement for basal Runx1 activity for tumor maintenance in the Eµ-Myc model of Burkitt's lymphoma. While normal Runx1fl/fl lymphoid cells permit mono-allelic deletion, primary Eµ-Myc lymphomas showed selection for retention of both alleles and attempts to enforce deletion in vivo led to compensatory expansion of p53null blasts retaining Runx1. Surprisingly, Runx1 could be excised completely from established Eµ-Myc lymphoma cell lines in vitro without obvious effects on cell phenotype. Established lines lacked functional p53, and were sensitive to death induced by introduction of a temperature-sensitive p53 (Val135) allele. Transcriptome analysis of Runx1-deleted cells revealed a gene signature associated with lymphoid proliferation, survival and differentiation, and included strong de-repression of recombination-activating (Rag) genes, an observation that was mirrored in a panel of human acute leukemias where RUNX1 and RAG1,2 mRNA expression were negatively correlated. Notably, despite their continued growth and tumorigenic potential, Runx1null lymphoma cells displayed impaired proliferation and markedly increased sensitivity to DNA damage and dexamethasone-induced apoptosis, validating Runx1 function as a potential therapeutic target in Myc-driven lymphomas regardless of their p53 status.

Using genomics for drug discovery in medulloblastoma

Tese de doutoramento, Medicina (Neurocirurgia), Universidade de Lisboa, Faculdade de Medicina, 2014

Application of Machine Learning Techniques to Acute Myeloid Leukemia

Thesis (Master's)--University of Washington, 2016-03

Identification de déterminants impliqués dans la différenciation des cellules souches embryonnaires

Les cellules souches ont attiré l’attention du public ces dernières années, grâce non-seulement à leur utilisation comme thérapies visant à s’attaquer à certains types de cancers, mais aussi en relation avec leur potentiel dans le domaine de la médecine regénérative. Il est établi que le destin cellulaire des cellules souches embryonnaires (ESC) est régulé de façon intensive par un groupe de facteur clés agissant sur leur pluripotence. Il est néanmoins envisageable que certains déterminants influençant l’auto-renouvellement et la différenciation de ces cellules soient toujours inconnus. Afin de tester cette hypothèse, nous avons généré, en utilisant une méthode par infections virales, une collection de ESC contenant des délétions chromosomales chevauchantes que nous avons baptisée DelES (Deletion in ES cells). Cette librairie contient plus de 1000 clones indépendants dont les régions délétées couvrent environ 25% du génome murin. À l’aide de cette ressource, nous avons conduit un criblage de formation de corps embryoïdes (EB), démontrant que plusieurs clones délétés avaient un phénotype de différenciation anormal. Nos études de complémentation sur un groupe de clones ont par la suite permis l’identification de Rps14 - un gène codant pour une protéine ribosomale (RP) comme étant haploinsuffisant pour la formation de EB. Dans un deuxième temps, l’analyse approfondie des résultats de notre crible a permis d’identifier un groupe de gènes codants pour des RP qui semblent essentiels pour la différenciation des ESC, mais dispensables pour leur auto-renouvellement. De manière intéressante, les phénotypes anormaux de formation en EB les plus marqués sont associés à des délétions de RP qui se retrouvent au site de sortie des ARN messagers (ARNm) du ribosome, soit Rps5, Rps14 et Rps28. Étonnament, alors qu’un débalancement des RP conduit généralement à une réponse de type p53, l’haploinsuffisance de ces trois gènes ne peut être renversée par une simple réduction des niveaux d’expression de ce gène suppresseur de tumeurs. Finalement, nos études de profilage polysomal et de séquençage à haut-débit montrent une signature spécifique de gènes liés au mésoderme chez un clone hétérozygote pour Rps5, suggérant ainsi une explication au phénotype de différenciation p53-indépendant identifié chez ces ESC. Nos travaux rapportent donc la création d’une ressource intéressante de génomique fonctionnelle qui a permis de mettre à jour le rôle essentiel que jouent les RP dans le processus de formation de EB. Nos résultats permettent aussi de documenter une réponse p53-indépendante suite à un débalancement de RP dans un contexte opposant l’auto-renouvellement et la différenciation des ESC.

Identificação de genes diferencialmente expressos em câncer de laringe

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Negative regulation of NF-κB by the ING4 tumor suppressor in breast cancer

Nuclear Factor kappa B (NF-κB) is a key mediator of normal immune response but contributes to aggressive cancer cell phenotypes when aberrantly activated. Here we present evidence that the Inhibitor of Growth 4 (ING4) tumor suppressor negatively regulates NF-κB in breast cancer. We surveyed primary breast tumor samples for ING4 protein expression using tissue microarrays and a newly generated antibody. We found that 34% of tumors expressed undetectable to low levels of the ING4 protein (n = 227). Tumors with low ING4 expression were frequently large in size, high grade, and lymph node positive, suggesting that down-regulation of ING4 may contribute to breast cancer progression. In the same tumor set, we found that low ING4 expression correlated with high levels of nuclear phosphorylated p65/RelA (p-p65), an activated form of NF-κB (p = 0.018). Fifty seven percent of ING4-low/p-p65-high tumors were lymph node-positive, indicating a high metastatic tendency of these tumors. Conversely, ectopic expression of ING4 inhibited p65/RelA phosphorylation in T47D and MCF7 breast cancer cells. In addition, ING4 suppressed PMA-induced cell invasion and NF-κB-target gene expression in T47D cells, indicating that ING4 inhibited NF-κB activity in breast cancer cells. Supportive of the ING4 function in the regulation of NF-κB-target gene expression, we found that ING4 expression levels inversely correlated with the expression of NF-κB-target genes in primary breast tumors by analyzing public gene expression datasets. Moreover, low ING4 expression or high expression of the gene signature composed of a subset of ING4-repressed NF-κB-target genes was associated with reduced disease-free survival in breast cancer patients. Taken together, we conclude that ING4 negatively regulates NF-κB in breast cancer. Consequently, down-regulation of ING4 leads to activation of NF-κB, contributing to tumor progression and reduced disease-free patient survival in breast cancer.

Inflammasome-dependent IFN-γ drives pathogenesis in Streptococcus pneumoniae meningitis

The pathology associated with Streptococcus pneumoniae meningitis results largely from activation of immune-associated pathways. We systematically investigated the production of IFN subtypes, as well as their influence on pathology, in a mouse model of S. pneumoniae meningitis. Despite the occurrence of a mixed IFN type I/II gene signature, no evidence for production or involvement of type I IFNs in disease progression was found. In contrast, type II IFN (IFN-γ) was strongly induced, and IFN-γ(-/-) mice were significantly protected from severe disease. Using intracellular cytokine staining and targeted cell-depletion approaches, NK cells were found to be the dominant source of IFN-γ. Furthermore, production of IFN-γ was found to be dependent upon ASC and IL-18, indicating that an ASC-dependent inflammasome pathway was responsible for mediating IFN-γ induction. The influence of IFN-γ gene deletion on a range of processes known to be involved in bacterial meningitis pathogenesis was examined. Although neutrophil numbers in the brain were similar in infected wild-type and IFN-γ(-/-) mice, both monocyte recruitment and CCL2 production were less in infected IFN-γ(-/-) mice compared with infected wild-type controls. Additionally, gene expression of NO synthase was strongly diminished in infected IFN-γ(-/-) mice compared with infected controls. Finally, bacterial clearance was enhanced in IFN-γ(-/-) mice, although the underlying mechanism remains unclear. Together, these data suggest that inflammasome-dependent IFN-γ contributes via multiple pathways to pathology during S. pneumoniae meningitis.

Circulating microRNA profiling in patients with advanced non-squamous NSCLC receiving bevacizumab/erlotinib followed by platinum-based chemotherapy at progression (SAKK 19/05).

OBJECTIVES Molecular subclassification of non small-cell lung cancer (NSCLC) is essential to improve clinical outcome. This study assessed the prognostic and predictive value of circulating micro-RNA (miRNA) in patients with non-squamous NSCLC enrolled in the phase II SAKK (Swiss Group for Clinical Cancer Research) trial 19/05, receiving uniform treatment with first-line bevacizumab and erlotinib followed by platinum-based chemotherapy at progression. MATERIALS AND METHODS Fifty patients with baseline and 24 h blood samples were included from SAKK 19/05. The primary study endpoint was to identify prognostic (overall survival, OS) miRNA's. Patient samples were analyzed with Agilent human miRNA 8x60K microarrays, each glass slide formatted with eight high-definition 60K arrays. Each array contained 40 probes targeting each of the 1347 miRNA. Data preprocessing included quantile normalization using robust multi-array average (RMA) algorithm. Prognostic and predictive miRNA expression profiles were identified by Spearman's rank correlation test (percentage tumor shrinkage) or log-rank testing (for time-to-event endpoints). RESULTS Data preprocessing kept 49 patients and 424 miRNA for further analysis. Ten miRNA's were significantly associated with OS, with hsa-miR-29a being the strongest prognostic marker (HR=6.44, 95%-CI 2.39-17.33). Patients with high has-miR-29a expression had a significantly lower survival at 10 months compared to patients with a low expression (54% versus 83%). Six out of the 10 miRNA's (hsa-miRN-29a, hsa-miR-542-5p, hsa-miR-502-3p, hsa-miR-376a, hsa-miR-500a, hsa-miR-424) were insensitive to perturbations according to jackknife cross-validation on their HR for OS. The respective principal component analysis (PCA) defined a meta-miRNA signature including the same 6 miRNA's, resulting in a HR of 0.66 (95%-CI 0.53-0.82). CONCLUSION Cell-free circulating miRNA-profiling successfully identified a highly prognostic 6-gene signature in patients with advanced non-squamous NSCLC. Circulating miRNA profiling should further be validated in external cohorts for the selection and monitoring of systemic treatment in patients with advanced NSCLC.

Ribosomal and Immune Transcripts Associate with Relapse in Acquired ADAMTS13-Deficient Thrombotic Thrombocytopenic Purpura.

Approximately 40% of patients who survive acute episodes of thrombotic thrombocytopenic purpura (TTP) associated with severe acquired ADAMTS13 deficiency experience one or more relapses. Risk factors for relapse other than severe ADAMTS13 deficiency and ADAMTS13 autoantibodies are unknown. ADAMTS13 autoantibodies, TTP episodes following infection or type I interferon treatment and reported ensuing systemic lupus erythematosus in some patients suggest immune dysregulation. This cross-sectional study asked whether autoantibodies against RNA-binding proteins or peripheral blood gene expression profiles measured during remission are associated with history of prior relapse in acquired ADAMTS13-deficient TTP. Peripheral blood from 38 well-characterized patients with autoimmune ADAMTS13-deficient TTP in remission was examined for autoantibodies and global gene expression. A subset of TTP patients (9 patients, 24%) exhibited a peripheral blood gene signature composed of elevated ribosomal transcripts that associated with prior relapse. A non-overlapping subset of TTP patients (9 patients, 24%) displayed a peripheral blood type I interferon gene signature that associated with autoantibodies to RNA-binding proteins but not with history of relapse. Patients who had relapsed bimodally expressed higher HLA transcript levels independently of ribosomal transcripts. Presence of any one potential risk factor (ribosomal gene signature, elevated HLA-DRB1, elevated HLA-DRB5) associated with relapse (OR = 38.4; p = 0.0002) more closely than any factor alone or all factors together. Levels of immune transcripts typical of natural killer (NK) and T lymphocytes positively correlated with ribosomal gene expression and number of prior episodes but not with time since the most recent episode. Flow cytometry confirmed elevated expression of cell surface markers encoded by these transcripts on T and/or NK cell subsets of patients who had relapsed. These data associate elevated ribosomal and immune transcripts with relapse history in acquired, ADAMTS13-deficient TTP.

Blocking the epithelial-to-mesenchymal transition pathway abrogates resistance to anti-folate chemotherapy in lung cancer.

Anticancer therapies currently used in the clinic often can neither eradicate the tumor nor prevent disease recurrence due to tumor resistance. In this study, we showed that chemoresistance to pemetrexed, a multi-target anti-folate (MTA) chemotherapeutic agent for non-small cell lung cancer (NSCLC), is associated with a stem cell-like phenotype characterized by an enriched stem cell gene signature, augmented aldehyde dehydrogenase activity and greater clonogenic potential. Mechanistically, chemoresistance to MTA requires activation of epithelial-to-mesenchymal transition (EMT) pathway in that an experimentally induced EMT per se promotes chemoresistance in NSCLC and inhibition of EMT signaling by kaempferol renders the otherwise chemoresistant cancer cells susceptible to MTA. Relevant to the clinical setting, human primary NSCLC cells with an elevated EMT signaling feature a significantly enhanced potential to resist MTA, whereas concomitant administration of kaempferol abrogates MTA chemoresistance, regardless of whether it is due to an intrinsic or induced activation of the EMT pathway. Collectively, our findings reveal that a bona fide activation of EMT pathway is required and sufficient for chemoresistance to MTA and that kaempferol potently regresses this chemotherapy refractory phenotype, highlighting the potential of EMT pathway inhibition to enhance chemotherapeutic response of lung cancer.

TAZ AS A REGULATOR OF MESENCHYMAL TRANSFORMATION AND CLINICAL AGGRESSIVENESS IN GLIOMAS

Glioblastoma multiforme (GBM) is an aggressive, high grade brain tumor. Microarray studies have shown a subset of GBMs with a mesenchymal gene signature. This subset is associated with poor clinical outcome and resistance to treatment. To establish the molecular drivers of this mesenchymal transition, we correlated transcription factor expression to the mesenchymal signature and identified transcriptional co-activator with PDZ-binding motif (TAZ) to be highly associated with the mesenchymal shift. High TAZ expression correlated with worse clinical outcome and higher grade. These data led to the hypothesis that TAZ is critical to the mesenchymal transition and aggressive clinical behavior seen in GBM. We investigated the expression of TAZ, its binding partner TEAD, and the mesenchymal marker FN1 in human gliomas. Western analyses demonstrated increased expression of TAZ, TEAD4, and FN1 in GBM relative to lower grade gliomas. We also identified CpG islands in the TAZ promoter that are methylated in most lower grade gliomas, but not in GBMs. TAZ-methylated glioma stem cell (GSC) lines treated with a demethylation agent showed an increase in mRNA and protein TAZ expression; therefore, methylation may be another novel way TAZ is regulated since TAZ is epigenetically silenced in tumors with a better clinical outcome. To further characterize the role of TAZ in gliomagenesis, we stably silenced or over-expressed TAZ in GSCs. Silencing of TAZ decreased invasion, self-renewal, mesenchymal protein expression, and tumor-initiating capacity. Over-expression of TAZ led to an increase in invasion, mesenchymal protein expression, mesenchymal differentiation, and tumor-initiating ability. These actions are dependent on TAZ interacting with TEAD since all these effects were abrogated with TAZ could not bind to TEAD. We also show that TAZ and TEAD directly bind to mesenchymal gene promoters. Thus, TAZ-TEAD interaction is critically important in the mesenchymal shift and in the aggressive clinical behavior of GBM. We identified TAZ as a regulator of the mesenchymal transition in gliomas. TAZ could be used as a biomarker to both estimate prognosis and stratify patients into clinically relevant subgroups. Since mesenchymal transition is correlated to tumor aggressiveness, strategies to target and inhibit TAZ-TEAD and the downstream gene targets may be warranted in alternative treatment.

Microarray analysis of autoimmune diseases by machine learning procedures

—Microarray-based global gene expression profiling, with the use of sophisticated statistical algorithms is providing new insights into the pathogenesis of autoimmune diseases. We have applied a novel statistical technique for gene selection based on machine learning approaches to analyze microarray expression data gathered from patients with systemic lupus erythematosus (SLE) and primary antiphospholipid syndrome (PAPS), two autoimmune diseases of unknown genetic origin that share many common features. The methodology included a combination of three data discretization policies, a consensus gene selection method, and a multivariate correlation measurement. A set of 150 genes was found to discriminate SLE and PAPS patients from healthy individuals. Statistical validations demonstrate the relevance of this gene set from an univariate and multivariate perspective. Moreover, functional characterization of these genes identified an interferon-regulated gene signature, consistent with previous reports. It also revealed the existence of other regulatory pathways, including those regulated by PTEN, TNF, and BCL-2, which are altered in SLE and PAPS. Remarkably, a significant number of these genes carry E2F binding motifs in their promoters, projecting a role for E2F in the regulation of autoimmunity.

Mouse p53-deficient cancer models as platforms for obtaining genomic predictors of human cancer clinical outcomes

Mutations in the TP53 gene are very common in human cancers, and are associated with poor clinical outcome. Transgenic mouse models lacking the Trp53 gene or that express mutant Trp53 transgenes produce tumours with malignant features in many organs. We previously showed the transcriptome of a p53-deficient mouse skin carcinoma model to be similar to those of human cancers with TP53 mutations and associated with poor clinical outcomes. This report shows that much of the 682-gene signature of this murine skin carcinoma transcriptome is also present in breast and lung cancer mouse models in which p53 is inhibited. Further, we report validated gene-expression-based tests for predicting the clinical outcome of human breast and lung adenocarcinoma. It was found that human patients with cancer could be stratified based on the similarity of their transcriptome with the mouse skin carcinoma 682-gene signature. The results also provide new targets for the treatment of p53-defective tumours.

Prior knowledge transfer across transcriptional data sets and technologies using compositional statistics yields new mislabelled ovarian cell line

Here, we describe gene expression compositional assignment (GECA), a powerful, yet simple method based on compositional statistics that can validate the transfer of prior knowledge, such as gene lists, into independent data sets, platforms and technologies. Transcriptional profiling has been used to derive gene lists that stratify patients into prognostic molecular subgroups and assess biomarker performance in the pre-clinical setting. Archived public data sets are an invaluable resource for subsequent in silico validation, though their use can lead to data integration issues. We show that GECA can be used without the need for normalising expression levels between data sets and can outperform rank-based correlation methods. To validate GECA, we demonstrate its success in the cross-platform transfer of gene lists in different domains including: bladder cancer staging, tumour site of origin and mislabelled cell lines. We also show its effectiveness in transferring an epithelial ovarian cancer prognostic gene signature across technologies, from a microarray to a next-generation sequencing setting. In a final case study, we predict the tumour site of origin and histopathology of epithelial ovarian cancer cell lines. In particular, we identify and validate the commonly-used cell line OVCAR-5 as non-ovarian, being gastrointestinal in origin. GECA is available as an open-source R package.

A global expression-based analysis of the consequences of the t(4;14) translocation in myeloma

Purpose: Our purpose in this report was to define genes and pathways dysregulated as a consequence of the t(4;14) in myeloma, and to gain insight into the downstream functional effects that may explain the different prognosis of this subgroup.Experimental Design: Fibroblast growth factor receptor 3 (FGFR3) overexpression, the presence of immunoglobulin heavy chain-multiple myeloma SET domain (IgH-MMSET) fusion products and the identification of t(4;14) breakpoints were determined in a series of myeloma cases. Differentially expressed genes were identified between cases with (n = 55) and without (n = 24) a t(4;14) by using global gene expression analysis.Results: Cases with a t(4;14) have a distinct expression pattern compared with other cases of myeloma. A total of 127 genes were identified as being differentially expressed including MMSET and cyclin D2, which have been previously reported as being associated with this translocation. Other important functional classes of genes include cell signaling, apoptosis and related genes, oncogenes, chromatin structure, and DNA repair genes. Interestingly, 25% of myeloma cases lacking evidence of this translocation had up-regulation of the MMSET transcript to the same level as cases with a translocation.Conclusions: t(4;14) cases form a distinct subgroup of myeloma cases with a unique gene signature that may account for their poor prognosis. A number of non-t(4;14) cases also express MMSET consistent with this gene playing a role in myeloma pathogenesis.