Explorative data analysis of MCL reveals gene expression networks implicated in survival and prognosis supported by explorative CGH analysis

BackgroundMantle cell lymphoma (MCL) is an incurable B cell lymphoma and accounts for 6% of all non-Hodgkin's lymphomas. On the genetic level, MCL is characterized by the hallmark translocation t(11;14) that is present in most cases with few exceptions. Both gene expression and comparative genomic hybridization (CGH) data vary considerably between patients with implications for their prognosis.MethodsWe compare patients over and below the median of survival. Exploratory principal component analysis of gene expression data showed that the second principal component correlates well with patient survival. Explorative analysis of CGH data shows the same correlation.ResultsOn chromosome 7 and 9 specific genes and bands are delineated which improve prognosis prediction independent of the previously described proliferation signature. We identify a compact survival predictor of seven genes for MCL patients. After extensive re-annotation using GEPAT, we established protein networks correlating with prognosis. Well known genes (CDC2, CCND1) and further proliferation markers (WEE1, CDC25, aurora kinases, BUB1, PCNA, E2F1) form a tight interaction network, but also non-proliferative genes (SOCS1, TUBA1B CEBPB) are shown to be associated with prognosis. Furthermore we show that aggressive MCL implicates a gene network shift to higher expressed genes in late cell cycle states and refine the set of non-proliferative genes implicated with bad prognosis in MCL.ConclusionThe results from explorative data analysis of gene expression and CGH data are complementary to each other. Including further tests such as Wilcoxon rank test we point both to proliferative and non-proliferative gene networks implicated in inferior prognosis of MCL and identify suitable markers both in gene expression and CGH data.

[1]  Hanno Steen,et al.  Development of human protein reference database as an initial platform for approaching systems biology in humans. , 2003, Genome research.

[2]  R. Gascoyne,et al.  Mantle cell lymphoma: a clinicopathologic study of 80 cases. , 1997, Blood.

[3]  D. Alberts,et al.  Variation in cyclooxygenase expression levels within the colorectum , 2003, Molecular carcinogenesis.

[4]  P. Grambsch,et al.  Martingale-based residuals for survival models , 1990 .

[5]  Paul Nurse,et al.  Cyclins and cell cycle control , 1991, Current Biology.

[6]  J. Davis.  Bioinformatics and Computational Biology Solutions Using R and Bioconductor , 2007 .

[7]  Jean YH Yang,et al.  Bioconductor: open software development for computational biology and bioinformatics , 2004, Genome Biology.

[8]  F. Nielsen,et al.  A Family of Insulin-Like Growth Factor II mRNA-Binding Proteins Represses Translation in Late Development , 1999, Molecular and Cellular Biology.

[9]  H. Herschman,et al.  Regulation of prostaglandin synthase-1 and prostaglandin synthase-2 , 1994, Cancer and Metastasis Reviews.

[10]  D. Koshland,et al.  The CDC20 gene product of Saccharomyces cerevisiae, a beta-transducin homolog, is required for a subset of microtubule-dependent cellular processes , 1991, Molecular and cellular biology.

[11]  Amy S. Lee,et al.  Stress induction of GRP78/BiP and its role in cancer. , 2006, Current molecular medicine.

[12]  H. Bourbon,et al.  Nucleolin, the major nucleolar protein of growing eukaryotic cells: an unusual protein structure revealed by the nucleotide sequence. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[14]  D. Pinkel,et al.  Mantle-cell lymphoma genotypes identified with CGH to BAC microarrays define a leukemic subgroup of disease and predict patient outcome. , 2005, Blood.

[15]  R. Eisenman,et al.  Deconstructing myc. , 2001, Genes & development.

[16]  R. Watson,et al.  B‐Myb protein in cellular proliferation, transcription control, and cancer: Latest developments , 1999, Journal of cellular physiology.

[17]  P. Nurse,et al.  Animal cell cycles and their control. , 1992, Annual review of biochemistry.

[18]  D. Treré,et al.  The quantity of nucleolar proteins nucleolin and protein B23 is related to cell doubling time in human cancer cells. , 1995, Laboratory investigation; a journal of technical methods and pathology.

[19]  B. Schaar,et al.  CENP-E is a putative kinetochore motor that accumulates just before mitosis , 1992, Nature.

[20]  Karin Milde-Langosch,et al.  The Fos family of transcription factors and their role in tumourigenesis. , 2005, European journal of cancer.

[21]  M. Cleary,et al.  PBX2 and PBX3, new homeobox genes with extensive homology to the human proto-oncogene PBX1 , 1991, Molecular and cellular biology.

[22]  P. Kaldis,et al.  Cdc2–cyclin E complexes regulate the G1/S phase transition , 2005, Nature Cell Biology.

[23]  A. Ristimäki,et al.  Cyclooxygenase Gene Expression in Inflammation and Angiogenesis a , 1993, Annals of the New York Academy of Sciences.

[24]  Rafael A. Irizarry,et al.  Bioinformatics and Computational Biology Solutions using R and Bioconductor , 2005 .

[25]  A. Patel,et al.  myc function and regulation. , 1992, Annual review of biochemistry.

[26]  Walter Krämer,et al.  Review of Modern applied statistics with S, 4th ed. by W.N. Venables and B.D. Ripley. Springer-Verlag 2002 , 2003 .

[27]  A. Bader,et al.  Molecular targets of the oncogenic transcription factor jun. , 2003, Current cancer drug targets.

[28]  Gordon K. Smyth,et al.  limma: Linear Models for Microarray Data , 2005 .

[29]  Gordon K Smyth,et al.  Statistical Applications in Genetics and Molecular Biology Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments , 2011 .

[30]  B. Calabretta,et al.  Cyclin D1-dependent regulation of B-myb activity in early stages of neuroblastoma differentiation , 2002, Cell Death and Differentiation.

[31]  Christian von Mering,et al.  STRING: known and predicted protein–protein associations, integrated and transferred across organisms , 2004, Nucleic Acids Res..

[32]  A. Rosenwald,et al.  The Ki67 proliferation index is a quantitative indicator of clinical risk in mantle cell lymphoma. , 2006, Blood.

[33]  Brian D. Ripley,et al.  Modern Applied Statistics with S Fourth edition , 2002 .

[34]  T. Kapoor,et al.  Correcting improper chromosome–spindle attachments during cell division , 2004, Nature Cell Biology.

[35]  F. Westermann,et al.  Neuroblastoma-derived sulfhydryl oxidase, a new member of the sulfhydryl oxidase/Quiescin6 family, regulates sensitization to interferon gamma-induced cell death in human neuroblastoma cells. , 2003, Cancer research.

[36]  C. Braak.  Canonical Correspondence Analysis: A New Eigenvector Technique for Multivariate Direct Gradient Analysis , 1986 .

[37]  K. Klempnauer,et al.  Regulation of B-Myb activity by cyclin D1 , 2000, Oncogene.

[38]  Ash A. Alizadeh,et al.  The lymphochip: a specialized cDNA microarray for the genomic-scale analysis of gene expression in normal and malignant lymphocytes. , 1999, Cold Spring Harbor symposia on quantitative biology.

[39]  A. Almasan,et al.  Opposing roles of E2Fs in cell proliferation and death , 2004, Cancer biology & therapy.

[40]  S. Akira,et al.  A nuclear factor for IL‐6 expression (NF‐IL6) is a member of a C/EBP family. , 1990, The EMBO journal.

[41]  F. Wilcoxon.  Individual Comparisons by Ranking Methods , 1945 .

[42]  E. Noordijk,et al.  Mantle-cell lymphoma: a population-based clinical study. , 1996, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[43]  D. Dewitt.  Prostaglandin endoperoxide synthase: regulation of enzyme expression. , 1991, Biochimica et biophysica acta.

[44]  L. Staudt,et al.  Specific secondary genetic alterations in mantle cell lymphoma provide prognostic information independent of the gene expression-based proliferation signature. , 2007, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[45]  Jörg Schultz,et al.  Genome Expression Pathway Analysis Tool – Analysis and visualization of microarray gene expression data under genomic, proteomic and metabolic context , 2007, BMC Bioinformatics.

[46]  K. Franssila,et al.  Ki‐67 expression level, histological subtype, and the International Prognostic Index as outcome predictors in mantle cell lymphoma , 2002, European journal of haematology.

[47]  C. Weiss,et al.  Deregulated Repression of c-Jun Provides a Potential Link to its Role in Tumorigenesis , 2004, Cell cycle.

[48]  Alexander F. Markham,et al.  ASPM is a major determinant of cerebral cortical size , 2002, Nature Genetics.

[49]  M. Barbacid,et al.  Cell cycle kinases in cancer. , 2007, Current opinion in genetics & development.

[50]  M. Introna,et al.  Independent regulation of c-myc, B-myb, and c-myb gene expression by inducers and inhibitors of proliferation in human B lymphocytes. , 1992, Journal of immunology.

[51]  G. Brittinger,et al.  Repp86: a new prognostic marker in mantle cell lymphoma , 2005, European journal of haematology.

[52]  D. Altieri,et al.  A novel anti-apoptosis gene, survivin, expressed in cancer and lymphoma , 1997, Nature Medicine.

[53]  C. Peschel,et al.  Targeting the proteasome in mantle cell lymphoma: A promising therapeutic approach , 2006, Leukemia & lymphoma.

[54]  Warren S. Alexander,et al.  A family of cytokine-inducible inhibitors of signalling , 1997, Nature.

[55]  L. Staudt,et al.  The proliferation gene expression signature is a quantitative integrator of oncogenic events that predicts survival in mantle cell lymphoma. , 2003, Cancer cell.

[56]  R. Prentice,et al.  Commentary on Andersen and Gill's "Cox's Regression Model for Counting Processes: A Large Sample Study" , 1982 .

[57]  G. Maga,et al.  Proliferating cell nuclear antigen (PCNA): a dancer with many partners , 2003, Journal of Cell Science.

[58]  M. Preece,et al.  Characterisation of the growth regulating gene IMP3, a candidate for Silver-Russell syndrome , 2002, Journal of medical genetics.

[59]  P. Knoepfler,et al.  HoxB8 requires its Pbx-interaction motif to block differentiation of primary myeloid progenitors and of most cell line models of myeloid differentiation , 2001, Oncogene.

[60]  D. Brautigan,et al.  Protein phosphatase 6 subunit with conserved Sit4-associated protein domain targets IkappaBepsilon. , 2006, The Journal of biological chemistry.

[61]  G. Gorbsky,et al.  Survivin dynamics increases at centromeres during G2/M phase transition and is regulated by microtubule-attachment and Aurora B kinase activity , 2004, Journal of Cell Science.

[62]  She Chen,et al.  Phosphorylation of Cdc20 by Bub1 provides a catalytic mechanism for APC/C inhibition by the spindle checkpoint. , 2004, Molecular cell.

[63]  K A Schafer,et al.  The Cell Cycle: A Review , 1998, Veterinary pathology.

[64]  A. Mulichak,et al.  The structure of mammalian cyclooxygenases. , 2003, Annual review of biophysics and biomolecular structure.

[65]  H. Pollard,et al.  Molecular dissection of nucleolin's role in growth and cell proliferation: new insights , 1999, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[66]  Stella Pelengaris,et al.  c-MYC: more than just a matter of life and death , 2002, Nature Reviews Cancer.

[67]  H. Friess,et al.  Cloning of a gene highly overexpressed in cancer coding for a novel KH-domain containing protein , 1997, Oncogene.

[68]  C. Meijer,et al.  Cell Cycle-controlled Interaction of Nucleolin with the Retinoblastoma Protein and Cancerous Cell Transformation* , 2006, Journal of Biological Chemistry.

[69]  D. Brautigan,et al.  Protein Phosphatase 6 Subunit with Conserved Sit4-associated Protein Domain Targets IκBϵ* , 2006, Journal of Biological Chemistry.