Full transcriptome analysis of rhabdomyosarcoma, normal and fetal skeletal muscle: statistical comparison of multiple SAGE libraries

Rhabdomyosarcoma (RMS) is the most frequent soft tissue sarcoma in children. Improved treatment strategies have increased overall survival, but the response of approximately one‐third of the patients is still poor. To increase the knowledge of RMS pathogenesis, we performed the first full transcriptome analysis of RMS using serial analysis of gene expression (SAGE). With a G‐test for the simultaneous comparison of subsets of SAGE libraries of normal skeletal muscle, embryonal (ERMS) and alveolar (ARMS) RMS, we identified 251 differentially expressed genes. A literature‐mining procedure demonstrated that 158 of these genes have not previously been associated with RMS or normal muscle. Gene Ontology (GO) analysis assigned 198 of the 251 genes to muscle‐specific classes, including those involved in normal myogenic development, as well as tumor‐related classes. Prominent GO classes were those associated with proliferation and actin reorganization, which are processes that play roles during early muscle development, muscle function, and tumor progression. Using custom microarrays, we confirmed the (up‐ or down‐) regulation of 80% of 98 differentially expressed genes. Another SAGE library of 19‐ to 22‐week‐old fetal skeletal muscle was compared with the RMS and normal muscle transcriptomes. Cluster analysis showed that the RMS and fetal muscle SAGE libraries formed one cluster distinct from normal muscle samples. Moreover, the expression profile of 86% of the differentially expressed genes between normal muscle and RMS was highly similar in fetal muscle and RMS. In conclusion, the G‐test is a robust tool for analyzing groups of SAGE libraries and correctly identifies genes marking the difference between fully differentiated skeletal muscle and RMS. This study not only substantiates the close association between embryonic myogenesis and RMS development but also provides a rich source of candidate genes to further elucidate the etiology of RMS or to identify diagnostic and/or prognostic markers.

[1]  Francesco Hofmann,et al.  Overexpression, genomic amplification and therapeutic potential of inhibiting the UbcH10 ubiquitin conjugase in human carcinomas of diverse anatomic origin , 2004, Oncogene.

[2]  Yvonne Braun,et al.  Profiling and functional annotation of mRNA gene expression in pediatric rhabdomyosarcoma and Ewing's sarcoma , 2004, International journal of cancer.

[3]  M. Tsao,et al.  Skp2 Gene Copy Number Aberrations Are Common in Non-Small Cell Lung Carcinoma, and Its Overexpression in Tumors with ras Mutation Is a Poor Prognostic Marker , 2004, Clinical Cancer Research.

[4]  J. Cameron,et al.  Analysis of novel tumor markers in pancreatic and biliary carcinomas using tissue microarrays. , 2004, Human pathology.

[5]  Javed Khan,et al.  Expression profiling identifies the cytoskeletal organizer ezrin and the developmental homeoprotein Six-1 as key metastatic regulators , 2004, Nature Medicine.

[6]  F. Baas,et al.  Expression of complement components in the peripheral nervous system. , 2003, Human molecular genetics.

[7]  Gayle M. Smythe,et al.  Notch-Mediated Restoration of Regenerative Potential to Aged Muscle , 2003, Science.

[8]  Y. Miyagi,et al.  Forced expression of antisense 14-3-3beta RNA suppresses tumor cell growth in vitro and in vivo. , 2003, Carcinogenesis.

[9]  Li Deng,et al.  Differential expression in SAGE: accounting for normal between-library variation , 2003, Bioinform..

[10]  Purvesh Khatri,et al.  Onto-Tools, the toolkit of the modern biologist: Onto-Express, Onto-Compare, Onto-Design and Onto-Translate , 2003, Nucleic Acids Res..

[11]  U. Heinzmann,et al.  Molecular characterization of Patched‐associated rhabdomyosarcoma , 2003, The Journal of pathology.

[12]  D. Monté,et al.  ICBP90 belongs to a new family of proteins with an expression that is deregulated in cancer cells , 2003, British Journal of Cancer.

[13]  M. Daly,et al.  PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes , 2003, Nature Genetics.

[14]  J. Ruijter,et al.  Statistical evaluation of SAGE libraries: consequences for experimental design. , 2002, Physiological genomics.

[15]  C. Marcelle,et al.  FGFR4 signaling is a necessary step in limb muscle differentiation. , 2002, Development.

[16]  J. Rowley,et al.  Identifying novel transcripts and novel genes in the human genome by using novel SAGE tags , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[17]  Oana A. Tomescu,et al.  Identification of target genes regulated by PAX3 and PAX3-FKHR in embryogenesis and alveolar rhabdomyosarcoma. , 2002, Genomics.

[18]  Ji Huang,et al.  [Serial analysis of gene expression]. , 2002, Yi chuan = Hereditas.

[19]  B. Smeets,et al.  Changes in mRNA expression profile underlie phenotypic adaptations in creatine kinase‐deficient muscles , 2001, FEBS letters.

[20]  P. Lollini,et al.  Identification of new genes related to the myogenic differentiation arrest of human rhabdomyosarcoma cells. , 2001, Gene.

[21]  M. Ringnér,et al.  Classification and diagnostic prediction of cancers using gene expression profiling and artificial neural networks , 2001, Nature Medicine.

[22]  D. Stekel,et al.  The comparison of gene expression from multiple cDNA libraries. , 2000, Genome research.

[23]  Yixin Wang,et al.  POWER_SAGE: comparing statistical tests for SAGE experiments , 2000, Bioinform..

[24]  G. Prendergast,et al.  The c-Myc-interacting adaptor protein Bin1 activates a caspase-independent cell death program , 2000, Oncogene.

[25]  J. Rowley,et al.  Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[26]  Michael L. Bittner,et al.  cDNA microarrays detect activation of a myogenic transcription program by the PAX3-FKHR fusion oncogene , 1999, Nature Genetics.

[27]  P M Bossuyt,et al.  Genes differentially expressed in medulloblastoma and fetal brain. , 1999, Physiological genomics.

[28]  J. Squire,et al.  Application of comparative genomic hybridization, spectral karyotyping, and microarray analysis in the identification of subtype-specific patterns of genomic changes in rhabdomyosarcoma. , 1999, Neoplasia.

[29]  M. G. Koerkamp,et al.  Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. , 1999, Molecular biology of the cell.

[30]  S. Welle,et al.  Inventory of high-abundance mRNAs in skeletal muscle of normal men. , 1999, Genome research.

[31]  M. Bittner,et al.  Gene expression profiling of alveolar rhabdomyosarcoma with cDNA microarrays. , 1998, Cancer research.

[32]  F. Baas,et al.  Allelotype of pediatric rhabdomyosarcoma , 1997, Oncogene.

[33]  J. Hsieh,et al.  Epstein-Barr virus immortalization: Notch2 interacts with CBF1 and blocks differentiation , 1997, Journal of virology.

[34]  F. Barr,et al.  The PAX3-FKHR fusion protein created by the t(2;13) translocation in alveolar rhabdomyosarcomas is a more potent transcriptional activator than PAX3 , 1995, Molecular and cellular biology.

[35]  G. Hannon,et al.  Correlation of terminal cell cycle arrest of skeletal muscle with induction of p21 by MyoD , 1995, Science.

[36]  T. Triche,et al.  Cytogenetic studies in subgroups of rhabdomyosarcoma , 1992, Genes, chromosomes & cancer.

[37]  M. Oshimura,et al.  Human chromosome 11 contains two different growth suppressor genes for embryonal rhabdomyosarcoma. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[38]  P. Tonin,et al.  Muscle-specific gene expression in rhabdomyosarcomas and stages of human fetal skeletal muscle development. , 1991, Cancer research.

[39]  F. James Rohlf,et al.  Biometry: The Principles and Practice of Statistics in Biological Research , 1969 .