If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Greehey Children's Cancer Research Institute, University of Texas Health Science Center, San Antonio, TexasDepartment of Cellular and Structural Biology, University of Texas Health Science Center, San Antonio, Texas
Division of Pediatric Neurooncology, German Cancer Research Center (DKFZ), Heidelberg, GermanyDepartment of Pediatric Oncology, Hematology and Immunology, Heidelberg University Hospital, Heidelberg, GermanyDivision of Neurosurgery, Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada
Department of Epidemiology and Biostatistics, University of Texas Health Science Center, San Antonio, TexasCancer Therapy and Research Center, University of Texas Health Science Center, San Antonio, Texas
CCU Neuropathology, German Cancer Research Center (DKFZ), Heidelberg, GermanyDepartment of Neuropathology, Heidelberg University Hospital, Heidelberg, Germany
Division of Pediatric Neurooncology, German Cancer Research Center (DKFZ), Heidelberg, GermanyDivision of Neurosurgery, Hospital for Sick Children, University of Toronto, Toronto, Ontario, Canada
Division of Pediatric Neurooncology, German Cancer Research Center (DKFZ), Heidelberg, GermanyDepartment of Pediatric Oncology, Hematology and Immunology, Heidelberg University Hospital, Heidelberg, Germany
Department of Molecular and Integrative Physiology, University of Kansas Medical Center, Kansas City, KansasUniversity of Kansas Cancer Center, Kansas City, Kansas
Address reprint requests to Luiz O.F. Penalva, Ph.D., University of Texas Health Science Center at San Antonio, 7703 Floyd Curl Dr., San Antonio, TX 78229-3900
Greehey Children's Cancer Research Institute, University of Texas Health Science Center, San Antonio, TexasDivision of Pediatric Neurooncology, German Cancer Research Center (DKFZ), Heidelberg, GermanyCancer Therapy and Research Center, University of Texas Health Science Center, San Antonio, Texas
Musashi1 (Msi1) is a highly conserved RNA-binding protein that is required during the development of the nervous system. Msi1 has been characterized as a stem cell marker, controlling the balance between self-renewal and differentiation, and has also been implicated in tumorigenesis, being highly expressed in multiple tumor types. We analyzed Msi1 expression in a large cohort of medulloblastoma samples and found that Msi1 is highly expressed in tumor tissue compared with normal cerebellum. Notably, high Msi1 expression levels proved to be a sign of poor prognosis. Msi1 expression was determined to be particularly high in molecular subgroups 3 and 4 of medulloblastoma. We determined that Msi1 is required for tumorigenesis because inhibition of Msi1 expression by small-interfering RNAs reduced the growth of Daoy medulloblastoma cells in xenografts. To characterize the participation of Msi1 in medulloblastoma, we conducted different high-throughput analyses. Ribonucleoprotein immunoprecipitation followed by microarray analysis (RIP-chip) was used to identify mRNA species preferentially associated with Msi1 protein in Daoy cells. We also used cluster analysis to identify genes with similar or opposite expression patterns to Msi1 in our medulloblastoma cohort. A network study identified RAC1, CTGF, SDCBP, SRC, PRL, and SHC1 as major nodes of an Msi1-associated network. Our results suggest that Msi1 functions as a regulator of multiple processes in medulloblastoma formation and could become an important therapeutic target.
RNA-binding proteins (RBPs) play a key role in gene regulation by modulating processes such as splicing, mRNA transport and stability, and translation. The human genome contains approximately 800 RBPs, many of which regulate specific groups of transcripts.
Some RBPs control mRNA subsets whose proteins are implicated in cancer-relevant processes, such as cell cycle, proliferation, apoptosis, and differentiation.
If the function or expression level of such RBPs is altered, the expression of multiple genes is affected, triggering or contributing to tumor formation.
An important example of such an RBP is Musashi1 (Msi1), a highly conserved protein with functions in nervous system development, stem cell self-renewal, and differentiation.
Msi1 can function either as a negative or positive regulator of translation by binding to specific motifs located in the 3′ untranslated regions (UTRs) of its target mRNAs.
Neural RNA-binding protein Musashi1 controls midline crossing of precerebellar neurons through posttranscriptional regulation of Robo3/Rig-1 expression.
Robo3 is a receptor involved in axonal guidance. Robo proteins direct the growing axons to the appropriate location and enable the crossing of the midline.
Expression of Msi1 was first observed in glioblastoma and medulloblastoma and has since been found in several cancers, such as hepatocellular carcinoma, cervical cancer, lung cancer, retinoblastoma, and breast, bladder, endometrial, and colon cancer.
Musashi1, an evolutionarily conserved neural RNA-binding protein, is a versatile marker of human glioma cells in determining their cellular origin, malignancy, and proliferative activity.
Enhanced therapeutic efficacy of oncolytic herpes vector G207 against human non-small cell lung cancer–expression of an RNA-binding protein Musashi1, as a marker for the tailored gene therapy.
Msi1 is particularly well characterized in breast and colon cancers. Approximately 40% of primary breast tumors and 100% of lymph node metastasis express increased levels of Msi1, which has been correlated to a poor prognosis.
The crypts of the intestines after exposure to toxic levels of 5-fluorouracil continued to express Msi1, suggesting that Msi1-positive cells are drug resistant.
Knockdown of Msi1 by RNA interference inhibited colon tumor cell growth as xenografts, whereas its ectopic expression in intestinal epithelium progenitor cells increased proliferation and tumorigenic capacity through the activation of Wnt and Notch pathways.
Musashi1, an evolutionarily conserved neural RNA-binding protein, is a versatile marker of human glioma cells in determining their cellular origin, malignancy, and proliferative activity.
In a previous study from our laboratory, we demonstrated that increased Msi1 expression is partially caused by loss or decreased expression of tumor suppressor microRNAs and via HuR, which increases mRNA stability and translation.
Musashi1, an evolutionarily conserved neural RNA-binding protein, is a versatile marker of human glioma cells in determining their cellular origin, malignancy, and proliferative activity.
Moreover, in an expression correlation study performed with The Cancer Genome Atlas glioblastoma data set, we identified Notch3 and Notch4 among the genes with the highest correlative expression with Msi1.
Knockdown of Msi1 in medulloblastoma cells affected their proliferation and growth as spheroids, suggesting a possible involvement with the so-called cancer stem cell population.
In this article, we report that Msi1 is a marker of poor prognosis in medulloblastoma. Msi1 seems to control tumor growth via a complex network of cancer-related genes. Altogether, our data demonstrate that Msi1 is a major player in medulloblastoma and may represent an important therapeutic target.
Materials and Methods
Cell Culture
The Daoy human medulloblastoma cell line and human embryonic kidney 293 (HEK293T) cells were obtained from ATCC (Manassas, VA). Cells were cultured using Dulbecco's modified Eagle medium (Hyclone, Logan, UT) supplemented with 10% fetal bovine serum (Atlanta Biologicals, Lawrenceville, GA) and penicillin and streptomycin (Gibco, Gaithersburg, MD) and incubated at 37°C in a humidified atmosphere of 5% carbon dioxide.
Tissue Arrays and Immunohistochemistry
Hematoxylin-eosin–stained sections from 259 paraffin blocks were prepared to define representative tumor regions. In addition, 10 samples of nontumor cerebellar tissues were included as a control. All tissue specimens were arrayed into a recipient block as previously described.
Evaluation of MSI1 (Abnova, Walnut, CA) immunostaining was performed in a semiquantitative manner (positive versus negative), and cases with clear evidence of cytoplasmic and/or nuclear staining were considered immunopositive by two investigators (A.K. and M.R.), who were blinded to clinical and molecular variables.
Xenografts
Xenograft models were established by using 5-week-old, male, athymic, nude mice. The mice were maintained with water and mouse chow ad libitum in a pathogen-free environment. To establish tumors, 5 × 107 Daoy cells were injected into the right and left flanks. Tumor size was measured by caliper, and volume was calculated as (length × width2) × 0.5. After 15 days, groups of five mice received three different treatments: Msi1 small-interfering (siRNA) injection, control siRNA injection, and nontreated control. siRNA was injected into the tumors every third day for a total of 5 doses. Msi1 siRNA (sequence targeting the coding region of Msi1; 5′-GGCUUUUGGAUUUGUGCAUdTdT-3′, where the lower case indicates deoxyribonucleotides) and a control siRNA were obtained from Ambion (Austin, TX). Efficacy of knockdown of Msi1 by siRNA was analyzed by real-time PCR and Western blot (see Supplemental Figure S1 at http://ajp.amjpathol.org). The siRNAs were mixed with 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) (Avanti Polar Lipids, Alabaster, AL) in the presence of excess t-butanol at a ratio of 1:10 (wt/wt) (siRNA/DOPC). Tween 20 (Sigma-Aldrich, St. Louis, MO) was added to the mixture in a ratio of 1:19 (Tween 20:siRNA/DOPC) and vortexed. The mixture was frozen in an acetone/dry ice bath and lyophilized. The siRNA preparation was reconstituted in 0.9% sterile saline before injecting at a dose of 5 μmol/L. At the end of treatment, animals were euthanized and the tumors removed and weighed. Size and weight between experiment and control tumors were compared using a paired two-tailed Student's t-test. P <0.05 was considered statistically significant. All values are expressed as mean ± SEM.
Immunohistochemistry of Xenografts
Tumors removed from xenografts were embedded in paraffin and cut into 4-μm sections and transferred to slides. The samples were deparaffinized, treated with a citrate buffer, and rehydrated in graded alcohol. Staining was completed by first blocking with avidin/biotin for 20 minutes and subsequent incubation with the appropriate antibody overnight at 4°C. Anti-Msi1 (Abcam, Cambridge, MA) and anti-Notch1 (Santa Cruz Biotechnology, Santa Cruz, CA) antibodies were used for immunohistochemistry. A secondary antibody conjugated to horseradish peroxidase was added to the slides, incubated for 30 minutes, and then developed with 3,3′-diaminobenzidine (Sigma-Aldrich). We used hematoxylin as a counterstain. Slides were dehydrated in graded alcohol (dH2O-95%), rinsed with xylene, and permanently mounted with Cryoseal Mounting Media 60 (Richard-Allen, Kalamazoo, MI).
Daoy cells were grown to 100% confluency in 150-mm dishes and then washed twice with ice-cold phosphate-buffered saline. Lysates were prepared with 1.5 volumes of polysomal lysis buffer [100 mmol/L KC, 25 mmol/L EDTA, 5 mmol/L MgCl2, 10 mmol/L HEPES, pH 7.0, 0.5% Nonidet P-40, 10% glycerol, 2 mmol/L dithiothreitol, 0.4 mmol/L vanadyl ribonucleoside complex, one tablet of complete protease inhibitor (Roche Applied Sciences, Indianapolis, IN), and RNaseOUT (Invitrogen, Carlsbad, CA)] and centrifuged, and the supernatant was removed and stored at −80°C. Protein A sepharose beads (GE Healthcare, Piscataway, NJ) were incubated overnight at 4°C with either anti-Msi1 polyclonal antibody or the preimmune serum.
The beads coated with antibody were resuspended in NT2 buffer supplemented with RNaseOUT, 0.2% vanadyl ribonucleoside complex, 2 mmol/L dithiothreitol, and 25 mmol/L EDTA. The cell lysate was added, and the mixture of beads and lysate was incubated at room temperature while tumbling end over end. After incubation, the beads were spun down and washed four times with ice-cold NT2 buffer. Subsequently, the material was digested with proteinase K, extracted with phenol/chloroform, and precipitated with isopropanol. RNA quality was analyzed by using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA) according to the manufacturer protocol.
Two hundred nanograms of RNA from each sample was labeled with Cy3 or Cy5 dye using an Agilent Low RNA Input Linear Amplification Kit (Agilent Technologies) according to manufacturer protocol. Dye-swap was performed on each sample to avoid dye-specific biases. After amplification, 750 ng of Cy3- and Cy5-labeled RNA were combined and hybridized to an Agilent Human Gene expression array (Agilent Technologies). After incubation, hybridized slides were washed according to manufacturer protocol. Microarray slides were then scanned on an Agilent Microarray Scanner G2565BA (Agilent Technologies), and data were extracted with Agilent Feature Extraction software version 9.1 (Agilent Technologies).
RIP-Chip Data Analysis
All analyses were performed within the R/Bioconductor software (R Foundation for Statistical Computing, Vienna, Austria). First, the quality of the microarrays was assessed by diagnostic plots (eg, a red/green density plot and an MA plot) and hierarchal clustering. The background was then corrected using a normexp + offset method.
Data were normalized within and between samples using the Loess and Aquantile methods, respectively. The LIMMA model was used to identify the mRNAs, which were preferentially associated with Msi1.
We used a conservative threshold of three-fold enrichment and an adjusted P<1 × 10−7. The microarray data were deposited in the NCBI Gene Expression Omnibus under the accession number GSE30904 (http://www.ncbi.nlm.nih.gov/geo).
Quantitative RT-PCR
We performed RIP-PCR to validate the RIP-chip data. Immunoprecipitated RNA was recovered and amplified using the NuGEN Ovation RNA Amplification kit (NuGEN Technologies, San Carlos, CA). Complementary DNA was synthesized using Applied Biosystems High Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Carlsbad, CA) using random priming. Quantitative PCR was performed using the TaqMan Gene Expression Master Mix (Applied Biosystems). All reactions were run on an ABI 7500 Real Time PCR machine (Applied Biosystems). The data were acquired using the ABI SDS software package version 2.0.1 (Applied Biosystems) and analyzed using the 2-ΔΔCt method.
Polysome Profiling
HEK293T cells were transfected in a 100-mm cell culture plate with an Msi1 expression construct and a green fluorescent protein (GFP) expression construct as negative control. Translation was arrested by adding 100 μg/mL of cycloheximide (Sigma-Aldrich) for 10 minutes for 72 hours after transfection. Cells were lysed in plate using a polysome lysis buffer (20 mmol/L Tris, pH 7.4, 250 mmol/L NaCl, 15 mmol/L MgCl2, 1 mmol/L dithiothreitol, 0.5% Triton X-100, 24 U/mL of Turbo DNase, and 100 μg/mL of cycloheximide). After scraping, the cells were triturated with a pipette. The lysate was incubated on ice for 10 minutes, clarified via centrifugation for 10 minutes at 11,592 × g at 4°C, and then overlaid on a 10% to 50% sucrose gradient in a buffer solution of 20 mmol/L Tris-HCl, pH 7.5, 100 mmol/L NaCl, 5 mmol/L MgCl2, and 0.1 mg/mL of cycloheximide. After ultracentrifugation at 113,422 × g at 4°C, 1-mL fractions were obtained on a Brandel density gradient fractionation system. A254 was measured during the entire fractionation process.
The monosomes and polysomes fractions were combined in separate tubes. RNA was extracted using the TRIzol reagent (Invitrogen) and resuspended in nuclease-free water (Ambion). Complementary DNA was synthesized using Applied Biosystems High Capacity cDNA Reverse Transcription Kit using random priming. Quantitative PCR was performed as described. Polysome association was calculated according to Peng et al,
Peng S, Chen LL, Lei XX, Yang L, Lin H, Carmichael GG, Huang Y: Genome-wide studies reveal that Lin28 enhances the translation of genes important for growth and survival of human embryonic stem cells. Stem Cells 29:496-504
with high polysome association taken as a proxy of high translation activity. In brief, the mRNA levels in the polysome and monosome fractions were calculated. Then, the percentage of mRNA in the polysome fraction relative to the total (polysome plus monosome fractions) was determined. The percentage values were normalized to the control GFP expression (setting the GFP value at 1).
Computational Data Analysis
To identify transcripts with similar expression as MSI1, we computed the Pearson correlation coefficients of the Msi1 expression profile across a medulloblastoma microarray study
to those of all other genes in the data set. Genes displaying strong positive correlation or anticorrelation (|R| > 0.5) with Msi1's profile were used in subsequent analysis.
To determine biological processes strongly associated with Msi1 and to explore the relationships of putative Msi1 targets, we analyzed all data sets with Pathway Studio 9 (Ariadne Genomics Inc., Rockville, MD), an automated text-mining tool based on published literature from various sources (eg, protein-protein interactions and gene co-expression).
Results
MSI1 Expression Is Elevated in Medulloblastoma and Is an Indicator of Poor Prognosis
We analyzed MSI1 expression in medulloblastoma versus normal cerebellar tissue using a large cohort of medulloblastoma samples (n = 103).
Overall, MSI1 expression levels in medulloblastoma are significantly higher than in adult cerebellum, where expression is restricted to the external granular cell layer (P < 0.001) (Figure 1A).
We also found that MSI1 expression is significantly elevated in medulloblastoma molecular subgroups 3 and 4, which are associated with a poor and intermediate prognosis, respectively, compared with the less severe WNT and SHH subgroups as defined by Taylor et al
Figure 1Distinct expression patterns of Msi1 and its prognostic relevance in medulloblastoma. A:MSI1 mRNA expression in the developing fetal cerebellum (FCB) and medulloblastoma (MB) compared with adult cerebellum (ACB) in the Toronto series (n = 103).
B:MSI1 mRNA expression is significantly elevated in both the group 3 and group 4 subgroups of medulloblastoma in the Toronto series (n = 103). MSI1 is elevated in the group 3 and group 4 subgroups of medulloblastoma compared with the WNT and SHH subgroups. C: Patterns of Msi1 protein expression: immunonegative (left panel) and immunopositive (right panel) examples on a medulloblastoma tissue microarray (n = 259). Kaplan-Meier plot of estimated overall survival (D) and progression-free survival (E) time distributions. The number of patients at risk is indicated for time increments of 12 months. Kaplan-Meier plots according to Msi1 expression in the entire cohort of 259 patients. *P ≤ 0.05, **P ≤ 0.01, and ***P ≤ 0.001.
Furthermore, Msi1 expression was evaluated by immunohistochemistry in 259 tumor samples on two medulloblastoma tissue microarrays. Among these samples, 59 (23%) were immunonegative, whereas the remaining 200 tumors (77%) demonstrated strong and homogeneous Msi1 expression (Figure 1C). Survival analysis by the Kaplan-Meier estimate and log-rank test revealed a significant association of Msi1 immunoreactivity with poor overall survival (P <0.0001) (Figure 1D) and progression-free survival (P <0.0001) (Figure 1E). We then delineated Msi1 expression in a subgroup-specific manner. For instance, we found that 98% of group 3 tumors have high Msi1 expression, whereas this number is only 53% in the case of WNT subgroup (see Supplemental Figure S3 at http://ajp.amjpathol.org). Moreover, when we focused on groups 3 and 4 tumors, we determined that high Msi1 expression denotes poor prognosis and overall survival compared with the Msi1-negative tumors (see Supplemental Figure S4 at http://ajp.amjpathol.org). These findings strongly support Msi1's role in tumorigenesis of high-risk medulloblastoma and, importantly, suggest its use as a prognostic marker.
Msi1 Contributes to Medulloblastoma Formation
We previously showed that Msi1 affects the viability of Daoy cells grown as tumor spheroids and cell proliferation by cell colony formation assay.
Building on these studies, we examined the role of Msi1 in medulloblastoma growth in xenografts of Daoy cells (Figure 2A). Fifteen days after cell implantation, tumors were treated with multiple injections of either a control (scrambled) or Msi1 siRNA. Tumor volume was measured at various time points, and tumor weight (Figure 2B) was determined on euthanasia. Tumors subjected to Msi1 siRNA treatment had a statistically significant decrease in tumor growth compared with the tumors injected with control siRNA and nontreated controls (P = 0.002) (Figure 2C). Tumors treated with control siRNA had a slight increase compared with the nontreated control, a similar effect that was previously observed.
Xenografts were immunohistochemically stained for Msi1 and Notch1 to corroborate the effects of the siRNA knockdown (Figure 2D). A decrease in Notch expression on Msi1 knockdown was expected based on the fact that Msi1 represses Numb, a Notch inhibitor.
Figure 2Msi1 silencing caused tumor regression of medulloblastoma cells grown as xenografts. Daoy medulloblastoma cells were implanted into the flank of nude mice and subjected to various treatments: nontreated control, scrambled siRNA injection, and Msi1 siRNA injection. Tumor size was measured at different time points, and tumor weight was calculated at the end of the experiment. A: Representative gross images of tumors from each of the treatment after extirpation. Scale bar = 1 cm. B: Graph representing the measurement of the tumor volume, in cubic millimeters, throughout the experiment. The data are shown as mean ± SE of mean (SEM) (n = 5). C: Tumor weights, in grams, at the final point of the experiment. Treatment with Msi1 siRNA resulted in reduction of tumor volume and weight. The data are shown as mean ± SEM (n = 5). Student's t-test was used to calculate P values. D: Immunohistochemical analysis of medulloblastoma xenograft tumors. Msi1 and Notch1 expression was assessed in each of the tumors. Representative images of each slide are shown. Notch1 expression is decreased when tumors are treated with Msi1 siRNA; thus, Notch1 is used as a readout for Notch signaling activity. Original magnification, ×40.
Genome-Wide Analyses Suggest that Msi1 Affects Tumorigenesis in Multiple Ways
Previous data from our laboratory indicated that through its binding to multiple RNA species, Msi1 controls a complex network of genes from multiple biological processes.
To investigate the impact of Msi1 in medulloblastoma, we conducted two genomic analyses. We started by mapping genes highly correlated (positively and negatively) with MSI1 expression by exploring the microarray data from our medulloblastoma cohort (n = 103, see Supplemental Table S1 at http://ajp.amjpathol.org). High correlation could indicate that the identified gene is a target or regulator of Msi1, the identified gene and Msi1 are co-regulated by the same factor, or the identified gene and Msi1 are functionally associated. Given the size of our data set, a Pearson correlation coefficient of ±0.5, corresponding to a P < 0.001, was used as cutoff for significant correlation or anticorrelation, resulting in 218 and 150 genes, respectively.
to identify mRNAs that are preferentially associated with Msi1 protein and to further determine how Msi1 influences tumorigenesis. Immunoprecipitations using Daoy cell extract were performed with anti-Msi1 antibodies and prebleed serum (see Supplemental Figure S5A at http://ajp.amjpathol.org).
Isolated mRNA populations were analyzed on microarrays. We used very stringent parameters (threefold minimum enrichment in experimental versus control samples and adjusted P ≤ 1 × 10−7) to filter the data. A summary plot presented in Supplemental Figure S5B (available at http://ajp.amjpathol.org) reveals the cutoff values and the genes with significant enrichment according to these criteria. The list includes 52 genes (Table 1; see also Supplemental Table S2 at http://ajp.amjpathol.org), 14 of which were previously identified in HEK293T cells.
We also noticed a substantial overlap (circa 60%) with mRNAs identified via RIP-chip in U251 glioblastoma cells (Vo et al, unpublished data). The presence of common targets across cell lines corroborates the efficiency of our approach to identify Msi1 target mRNAs. Data were validated by RIP-PCR (see Supplemental Figure S6 at http://ajp.amjpathol.org).
Table 1Genes Identified by RIP-chip are Associated with Cancer-Related Processes
For each target 3′ UTR, we counted the number of occurrences of the GUAG and UAG oligomers and normalized each count by the length of the UTR. We did the same for nontarget 3′ UTRs. On average, GUAG occurs 30% more frequently and UAG 50% more frequently in the target 3′ UTRs when compared with nontarget 3′ UTRs (P < 1.6 × 10−6 and < 9.2 × 10−20, respectively; Wilcoxon rank-sum test). Furthermore, we examined the most enriched targets in the RIP-chip data and observed the same distinct pattern of AUG and GUAG oligomers separated by a variable length linker region as was recently reported
(Figure 3). These results lend further support both for the putative Msi1 binding motif and for the accuracy of the identified targets.
Figure 3Identification of putative Msi1 binding sites. Highlighted in bold are the putative binding sites for Msi1. GUAG is the putative binding sequence for the RNA recognition motif 1 (RRM1) of Msi1, and UAG corresponds to the putative binding site for Msi1 RRM2.
Pathway Studio 9 determined significant correlation between identified Msi1 associated mRNAs and cancer relevant processes [apoptosis (∼44%), cell cycle (∼24%), cell differentiation (∼40%), cell proliferation (∼46%), cell survival (∼36%), and DNA repair (∼7%)]. All observed values showed enrichment over total transcriptome (P < 0.1–0.001) (Figure 4 and Table 1).
Figure 4Msi1 putative targets are enriched for cancer-related process. We used Pathway Studio 9 to analyze the biological processes connected with Msi1-associated mRNAs. The graph shows the percentage of genes associated with each biological process in the Msi1 subset (gray) and in total transcriptome (black). RIP-chip analysis of Msi1 results in an overall enrichment of mRNAs related to apoptosis, cell cycle, cell differentiation, cell proliferation, cell survival, and DNA repair versus the percentage of mRNAs that resides in whole transcriptome that are connected to cancer-related biological processes. Fisher's exact test was used to calculate P values.
To have a broader view of Msi1's impact on medulloblastoma and to identify the main mediators of its function, we conducted an in silico network analysis of the functional relationships between Msi1's associated mRNAs and correlated genes. Links between genes are drawn by Pathway Studio based on regulation, functional, and biological associations extracted from the literature (Figure 5; see also Supplemental Table S3 at http://ajp.amjpathol.org).
Figure 5Network analysis of Msi1. Blue ovals indicate Msi1 RIP-chip targets in Daoy cells (this study), and red ovals are genes whose expression is correlated to Msi1 (Pearson's R > 0.5). Links between nodes denote a functional relationship between the two genes as reported by literature, co-expression analysis, regulation, and protein-protein interaction analysis. The edges are directed to denote regulator-target relationships as determined by local connectivity. The figure illustrates the high connectivity of genes controlled by Msi1 expression.
Msi1 Potentially Regulates Translation of Target mRNAs
To determine changes in translation efficiency as a consequence of Msi1 binding, we performed a polysome association experiment. HEK293T cells were transfected with either a GFP (negative control) or a Msi1 expression construct.
Cellular lysate was prepared and the polysome fraction (which contains highly translated mRNAs) and monosome fraction (which contains poorly translated mRNAs) were separated by ultracentrifugation in a sucrose density gradient. Quantitative RT-PCR was then used to measure the impact of Msi1 transgenic expression on polysome association of a group mRNAs in experimental versus control samples (Figure 6). We focused on a set of mRNAs that had been identified as preferentially associated with Msi1 both in Daoy and HEK293T cells,
thus representing highly likely Msi1 targets. We observed that AASDHPPT, ATP6AP2, ATP6V1G1, and KIAA0101 have decreased polysome association on Msi1 expression, suggesting translational repression by Msi1 (as has been observed for other genes
Figure 6Msi1 affects ribosomal profile of associated mRNAs. 293T cells were transfected with expression constructs for Msi1 (pEF1-Msi1) or GFP (pEF1-GFP). After 72 hours of transfection, translation was arrested with cycloheximide, and cell lysates were prepared and analyzed on a sucrose density gradient. The polysome association for each selected mRNA was calculated according to Peng et al
Peng S, Chen LL, Lei XX, Yang L, Lin H, Carmichael GG, Huang Y: Genome-wide studies reveal that Lin28 enhances the translation of genes important for growth and survival of human embryonic stem cells. Stem Cells 29:496-504
and compared between Msi1 and GFP. GAPDH and ACTB, which were not identified to be associated with Msi1, were used as a negative control. Transgenic expression of Msi1 repressed polysome association of four of the five analyzed target mRNAs without affecting GAPDH and ACTB. Data were analyzed with Student's t-test and are presented as mean ± SEM. Experiments were performed in triplicate. *P < 0.01, ***P < 0.001.
The design of more efficient therapies requires a more elaborate study of this cell population, including the mapping and characterization of key regulators such as Msi1 that maintain the balance between self-renewal and differentiation under normal circumstances but can trigger or contribute to tumor formation when aberrantly expressed. We have previously determined that Msi1 expression is elevated in Daoy medulloblastoma cells grown as spheroids, which are enriched with cells with stem-like characteristics.
Our data suggest that Msi1 is a marker of poor prognosis of medulloblastoma because patients displaying tumors with high Msi1 expression had a worse overall survival and progression-free survival than patients with tumors that had lower Msi1 expression. In agreement with this observation, we demonstrated that Msi1 silencing in tumor xenografts caused an impact on tumor growth. Therefore, we conclude that Msi1 is not only a marker of poor prognosis in medulloblastoma but also is required for tumorigenesis.
Recent analyses of a large cohort of medulloblastoma gene expression data have defined four molecular subgroups of medulloblastoma with different prognostic outcomes (SHH, WNT, group 3, and group 4).
By examining the expression of MSI1 in this data set, we determined that high MSI1 expression is prevalent in subgroups 3 and 4, being particularly high in group 4. Groups 3 and 4 display a clear class separation from the other two subgroups of medulloblastoma, WNT and SHH. This stratification of MSI1 expression correlates with patients in groups 3 and 4 having a poor and intermediate prognosis, respectively, compared with the WNT and SHH subgroups, which have relatively good prognosis.
Msi1 seems to influence tumorigenesis in multiple ways as suggested by our RIP-chip studies. We focus the discussion on the major nodes identified by analysis of a functional network. We discuss selected genes whose degree in the network is higher than the median degree. Among the RIP-chip identified targets, we highlight RAC1, CTGF, and SDCBP. Rac1 is a member of the Ras superfamily of GTP-binding proteins that has many cellular roles, such as regulating gene expression, cell proliferation, cytoskeleton reorganization, and metastasis.
In many cancers, CTGF functions in biological processes, such as cell proliferation, migration, adhesion, extracellular matrix formation, and angiogenesis.
CTGF overexpression is found in many malignant tumors, such as melanoma, prostate cancer, pancreatic cancer, esophageal carcinoma, lymphoblastic leukemia, gastric cancer, and gliomas.
SDCBP, or syntenin, is an important molecule in connecting syndecan signaling to the cytoskeleton through the binding of its PDZ motif to syndecan proteoglycans.
This interaction allows cancer cells to acquire the capability to migrate and invade; in particular, this capability is garnered through the activation of c-Src, focal adhesion kinase, MAPK pathway, and NF-κB by syntenin.
Activation of the integrin effector kinase focal adhesion kinase in cancer cells is regulated by crosstalk between protein kinase C-alpha and the PDZ adapter protein mda-9/Syntenin.
In addition to its Msi1 interaction observed by RIP-chip analysis, SDCBP also correlates negatively with MSI1 gene expression (P < 0.0002).
Three major nodes come from the analysis of expression correlation: SRC, PRL, and SHC1. SRC (c-Src) is one of the first proto-oncogenes discovered and provided a launching point for the study of the molecular basis of cancer.
c-Src is one of the members of the SRC family of kinases and is implicated in the development and progression of various human cancers. Elevated expression of c-Src is observed in a number of different cancer types.
c-Src has a wide range of different kinase substrates, which includes proteins such as RasGAP SH3-domain binding protein, Ewing sarcoma breakpoint region 1, and heteronuclear ribonucleoprotein A2/B1.
The cellular effects of c-Src are evident and widespread. c-Src expression can transform normal fibroblast so that the cell is now able to proliferate faster, lose density-dependent growth inhibition, and form tumors in in vivo xenograft models. Like many other tumor types, c-Src is overexpressed in medulloblastoma.
Because c-Src expression is so potent and a potential therapeutic target, inhibition of c-Src kinase activity with novel inhibitors, such as pyrazolo-[3,4-day]-pyrimidine derivatives, can reduce the growth rate of medulloblastoma cells by arresting cell cycle progression at the G2/M transition and induce apoptosis.
PRL is a luteotropic hormone protein whose best characterized effect is the stimulation of mammary glands to produce milk. In the maternal central nervous system, PRL can stimulate the proliferation of oligodendrocyte precursor cells and enhance remyelination.
Evo-Devo and the evolution of cancer: a hypothesis for metamorphic therapies for the cancers of prolactin-influenced tumourigenesis: with special reference to glioblastoma multiforme (GBM).
SHC1 (also known as SHCA) is a signal transduction adaptor protein that relays extracellular signals downstream of receptor tyrosine kinases. SHC1 has a wide range of functions through its promotion of proliferation, angiogenesis, migration, and invasion.
In medulloblastoma, SHC1 expression is elevated in metastatic medulloblastoma samples compared with nonmetastatic samples, suggesting a potential role of SHC1 in mediating metastasis of medulloblastoma.
In summary, our study establishes a functional role of Msi1 in medulloblastoma through posttranscriptional control of a subset of cancer-related mRNAs. Studies such as this underline the importance of RBPs and posttranscriptional gene regulation in the molecular pathogenesis of cancer. We believe that further studies will likely demonstrate that more RBPs are directly involved in tumorigenesis. In fact, in a previous in silico study from our laboratory, we identified 53 RBPs aberrantly expressed in multiple tumor types.
A comprehensive in silico expression analysis of RNA binding proteins in normal and tumor tissue: identification of potential players in tumor formation.
Further systematic studies are required to understand the role of posttranscriptional gene regulation in tumorigenesis and its interaction with known cancer signaling pathways and to identify novel therapeutic targets for cancer treatment.
Assessing the efficacy of MSI1 knockdown by siRNA. A: Quantitative RT-PCR analysis of MSI1 knockdown using siRNAs. Daoy medulloblastoma cells were transfected with a control siRNA (at 50 nmol/L) or increasing concentrations of MSI1 siRNAs (10, 50, 100 nmol/L). Untreated Daoy cells are denoted as control. B: Western blot analysis of Msi1 knockdown using siRNAs. Control or MSI1 siRNAs were transfected into Daoy medulloblastoma cells and analyzed by Western blot. Concentrations of each siRNA are listed above each lane. Actin was used as a loading control.
MSI1 expression pattern in the Heidelberg medulloblastoma gene expression data set54 (n = 64). MSI1 is elevated in the group 3 and group 4 subtypes of medulloblastoma compared with the WNT and SHH subgroups.
Subgroup-specific delineation of Msi1 expression of tissue microarray data of medulloblastoma. Tumors were categorized into specific molecular subgroups of medulloblastoma and parsed into subcategories of high and low Msi1 expression. The percentage of each subcategory for each medulloblastoma subgroup is shown.
Kaplan-Meier analysis of overall survival of non-WNT/-SHH tumors and Msi1 expression. A total of 152 tumors of group 3 and group 4 were analyzed for Msi1 expression and correlated to Msi1 expression. Patients with tumors with positive Msi1 expression had a poorer overall survival compared with patients with tumors with negative Msi1 expression (P = 0.0187).
A: Immunoprecipitation-Western blotting demonstrates that the Msi1 antibody specifically immunoprecipitates Msi1 versus that of the prebleed serum. B: Summary plot shows the cutoff value for the RIP-chip data using an adjusted P ≤ 1 × 107 and a fold change of ≥3. Dots in red represent mRNAs that are preferentially associated with Msi1.
Validation of RIP-chip data using RIP-PCR. RIP was performed using a Msi1-specific antibody and IgG isotype antibody control. Bound RNA was extracted and amplified. Quantitative RT-PCR was used to analyze Msi1 target mRNAs, and the data were normalized to the IgG isotype immunoprecipitation.
Neural RNA-binding protein Musashi1 controls midline crossing of precerebellar neurons through posttranscriptional regulation of Robo3/Rig-1 expression.
Musashi1, an evolutionarily conserved neural RNA-binding protein, is a versatile marker of human glioma cells in determining their cellular origin, malignancy, and proliferative activity.
Peng S, Chen LL, Lei XX, Yang L, Lin H, Carmichael GG, Huang Y: Genome-wide studies reveal that Lin28 enhances the translation of genes important for growth and survival of human embryonic stem cells. Stem Cells 29:496-504
Activation of the integrin effector kinase focal adhesion kinase in cancer cells is regulated by crosstalk between protein kinase C-alpha and the PDZ adapter protein mda-9/Syntenin.
Evo-Devo and the evolution of cancer: a hypothesis for metamorphic therapies for the cancers of prolactin-influenced tumourigenesis: with special reference to glioblastoma multiforme (GBM).
A comprehensive in silico expression analysis of RNA binding proteins in normal and tumor tissue: identification of potential players in tumor formation.
Supported by grants from the Children's Brain Tumor Foundation, Association for Research of Childhood Cancer, The Max and Minnie Tomerlin Voelcker Fund, and NIH grant R01HG006015 (L.O.F.P.).
D.T.V., D.S., and M.R. contributed equally to this work.
Genes previously identified in a study with HEK293T cells.
00601-3&id=gr1.jpg){kind=link}
00601-3&id=gr2.jpg){kind=link}
00601-3&id=gr3.jpg){kind=link}
00601-3&id=gr4.jpg){kind=link}
00601-3&id=gr5.jpg){kind=link}
00601-3&id=gr6.jpg){kind=link}