Identification and modelling of human breast cancer subtypes

Résumé Le cancer du sein est le cancer le plus commun chez les femmes et est responsable de presque 30% de tous les nouveaux cas de cancer en Europe. On estime le nombre de décès liés au cancer du sein en Europe est à plus de 130.000 par an. Ces chiffres expliquent l'impact social considérable de cette maladie. Les objectifs de cette thèse étaient: (1) d'identifier les prédispositions et les mécanismes biologiques responsables de l'établissement des sous-types spécifiques de cancer du sein; (2) les valider dans un modèle ín vivo "humain-dans-souris"; et (3) de développer des traitements spécifiques à chaque sous-type de cancer du sein identifiés. Le premier objectif a été atteint par l'intermédiaire de l'analyse des données d'expression de gènes des tumeurs, produite dans notre laboratoire. Les données obtenues par puces à ADN ont été produites à partir de 49 biopsies des tumeurs du sein provenant des patientes participant dans l'essai clinique EORTC 10994/BIG00-01. Les données étaient très riches en information et m'ont permis de valider des données précédentes des autres études d'expression des gènes dans des tumeurs du sein. De plus, cette analyse m'a permis d'identifier un nouveau sous-type biologique de cancer du sein. Dans la première partie de la thèse, je décris I identification des tumeurs apocrines du sein par l'analyse des puces à ADN et les implications potentielles de cette découverte pour les applications cliniques. Le deuxième objectif a été atteint par l'établissement d'un modèle de cancer du sein humain, basé sur des cellules épithéliales mammaires humaines primaires (HMECs) dérivées de réductions mammaires. J'ai choisi d'adapter un système de culture des cellules en suspension basé sur des mammosphères précédemment décrit et pat décidé d'exprimer des gènes en utilisant des lentivirus. Dans la deuxième partie de ma thèse je décris l'établissement d'un système de culture cellulaire qui permet la transformation quantitative des HMECs. Par la suite, j'ai établi un modèle de xénogreffe dans les souris immunodéficientes NOD/SCID, qui permet de modéliser la maladie humaine chez la souris. Dans la troisième partie de ma thèse je décris et je discute les résultats que j'ai obtenus en établissant un modèle estrogène-dépendant de cancer du sein par transformation quantitative des HMECs avec des gènes définis, identifiés par analyse de données d'expression des gènes dans le cancer du sein. Les cellules transformées dans notre modèle étaient estrogène-dépendantes pour la croissance, diploïdes et génétiquement normales même après la culture cellulaire in vitro prolongée. Les cellules formaient des tumeurs dans notre modèle de xénogreffe et constituaient des métastases péritonéales disséminées et du foie. Afin d'atteindre le troisième objectif de ma thèse, j'ai défini et examiné des stratégies de traitement qui permettent réduire les tumeurs et les métastases. J'ai produit un modèle de cancer du sein génétiquement défini et positif pour le récepteur de l'estrogène qui permet de modéliser le cancer du sein estrogène-dépendant humain chez la souris. Ce modèle permet l'étude des mécanismes impliqués dans la formation des tumeurs et des métastases. Abstract Breast cancer is the most common cancer in women and accounts for nearly 30% of all new cancer cases in Europe. The number of deaths from breast cancer in Europe is estimated to be over 130,000 each year, implying the social impact of the disease. The goals of this thesis were first, to identify biological features and mechanisms --responsible for the establishment of specific breast cancer subtypes, second to validate them in a human-in-mouse in vivo model and third to develop specific treatments for identified breast cancer subtypes. The first objective was achieved via the analysis of tumour gene expression data produced in our lab. The microarray data were generated from 49 breast tumour biopsies that were collected from patients enrolled in the clinical trial EORTC 10994/BIG00-01. The data set was very rich in information and allowed me to validate data of previous breast cancer gene expression studies and to identify biological features of a novel breast cancer subtype. In the first part of the thesis I focus on the identification of molecular apacrine breast tumours by microarray analysis and the potential imptìcation of this finding for the clinics. The second objective was attained by the production of a human breast cancer model system based on primary human mammary epithelial cells {HMECs) derived from reduction mammoplasties. I have chosen to adopt a previously described suspension culture system based on mammospheres and expressed selected target genes using lentiviral expression constructs. In the second part of my thesis I mainly focus on the establishment of a cell culture system allowing for quantitative transformation of HMECs. I then established a xenograft model in immunodeficient NOD/SCID mice, allowing to model human disease in a mouse. In the third part of my thesis I describe and discuss the results that I obtained while establishing an oestrogen-dependent model of breast cancer by quantitative transformation of HMECs with defined genes identified after breast cancer gene expression data analysis. The transformed cells in our model are oestrogen-dependent for growth; remain diploid and genetically normal even after prolonged cell culture in vitro. The cells farm tumours and form disseminated peritoneal and liver metastases in our xenograft model. Along the lines of the third objective of my thesis I defined and tested treatment schemes allowing reducing tumours and metastases. I have generated a genetically defined model of oestrogen receptor alpha positive human breast cancer that allows to model human oestrogen-dependent breast cancer in a mouse and enables the study of mechanisms involved in tumorigenesis and metastasis.

Identification, using a siRNA screen approach, and molecular characterization of a new Fas pathway modulator: STK11/

Abstract : Apoptosis is an evolutionarily conserved cellular suicide mechanism that can be triggered by activation of various pathways, such as the Fas-Pathway. Upon stimulation by its specific ligand (FasL), present at the surface of Cytotoxic Τ lymphocytes, the death receptor Fas initiates a signaling cascade culminating in the activation of cellular caspases, leading thus to cell death of the target cell (e.g. transformed cell). Dysregulation of apoptosis in general, and of Fas pathway in particular, was shown to contribute to pathogenesis of cancers and many human diseases. Even though, during the last decades the molecular mechanisms of apoptosis have been widely studied, it is important to better understand the mechanisms leading to apoptosis, to improve our understanding of pathological processes, and generate more subtle apoptosis-modulating therapies to fight cancer and other diseases. In order to identify new components of the Fas signaling pathway, a screen based on the mechanism of RNA interference was undertaken. After a first and a second manual whole-kinome screen, we identified several strong positive hits that showed a protection against Fas ligand-induced apoptosis with distinct siRNAs, notably STK11, an interesting tumor suppressor mutated in several sporadic and inherited cancers. The STK11 functional characterization reveals that this kinase represents an apically acting general pro-apoptotic modulator of the extrinsic pathway (FasL, TRAIL, TNF-induced apoptosis), but not of the intrinsic apoptotic pathway. The STK11 action on the Fas pathway was shown to be dependent on its kinase activity, but independent of AMPK, a well-characterized STK11 downstream substrate. Furthermore, STK11 was shown to interact with caspase-8, a major mediator of the extrinsic pathway, and modulate its activity through an unclear mechanism that may involve an STK11-dependant caspase-8 phosphorylation. This modification may allow a proper caspase-8 polyubiquitination and activation in p62 sequestosmes aggregates, but may also increase the activation of caspase-8 at the DISC level. In addition, we observed that STK11 modulate not only the apoptotic pathway induced by Fas engagement, but also FasL-induced JNK and NF- KB, sustaining an upstream role of this kinase in the pathway. In conclusion, our report reveals that STK11 is an important pro-apoptotic modulator of the Fas pathway in particular, and extrinsic pathway in general. Our finding could explain, at least partially, why inactivating mutations of the kinase leads to cancer, by allowing resistance to apoptosis and accordingly evasion of immune surveillance. Résumé : L'apoptose est un mécanisme de suicide cellulaire, conservé dans diverses espèces, et qui au niveau moléculaire est déclenché par différentes voies de signalisation, comme par exemple lors de l'activation du récepteur Fas. La liaison du ligand FasL au récepteur de la mort Fas, induit une cascade de signalisation qui conduit à l'activation des caspases. Les lymphocytes Τ cytotoxiques peuvent utiliser la voie Fas pour induire la mort et se débarrasser de cellules dangereuses pour le reste de l'organisme, tel que les cellules transformées. La dysrégulation de l'apoptose en général, et de la voie Fas en particulier, peut contribuer à diverses maladies telles que le cancer. Même si ces dernières décennies, les mécanismes moléculaires conduisant à l'apoptose ont été extensivement étudiés, il reste néanmoins important de mieux comprendre le phénomène d'apoptose, pour améliorer notre compréhension des processus pathologiques, mais surtout dans le but de développer de nouvelles thérapies ciblant l'apoptose contre le cancer et d'autres pathologies. Pour identifier de nouveau constituants de la voie Fas, un criblage génétique basé sur l'interférence à l'ARN a été entrepris. Après un premier et un deuxième criblage d'une librairie du kinome, nous avons identifié différentes protéines qui pourraient jouer un rôle positif dans la voie Fas, et en particulier la protéine suppresseur de tumeur STK11, qui est fréquemment mutée dans divers cancers sporadiques et héréditaires. La caractérisation fonctionnelle de STK11 a révélé que cette kinase était un modulateur apical de la voie extrinsèque de l'apoptose en général (Fas, TNF, TRAIL), mais pas de la voie intrinsèque. L'action de STK11 sur la voie Fas est dépendante de sa fonction kinase, mais indépendante de l'AMPK, un substrat bien caractérisé de STK11. De plus, STK11 interagît avec la caspase-8, un constituant majeur de la voie Fas, et module son activité, par un mécanisme encore peu clair qui pourrait impliquer une phosphorylation de la caspase-8 par STK11. Cette modification pourrait permettre une activation optimale de la caspase-8 en jouant un rôle dans le processus de polyubiquitination de la caspase-8, phénomène qui semble être important pour l'activation de la caspase-8 dans des agrégats protéiques avec p62, mais qui pourrait aussi augmenter son activation au niveau du DISC. Finalement, nous avons observé que STK11 modulait non seulement la voie apoptotique déclenchée par l'activation de Fas, mais aussi les voies non-apoptotiques de Fas, comme JNK et NF-KB. En conclusion notre étude, révèle que STK11 est un important modulateur pro- apoptotique de la voie Fas, et de la voie extrinsèque en général. Cette découverte pourrait expliquer, du moins partiellement, pourquoi les mutations inactivatrices de STK11 conduisent au cancer, par une augmentation de la résistance à l'apoptose et donc par l'évasion de la surveillance immunitaire.

Performance of blood pressure-to-height ratio at a single screening visit for the identification of hypertension in children.

BACKGROUND: The diagnosis of hypertension in children is difficult because of the multiple sex-, age-, and height-specific thresholds to define elevated blood pressure (BP). Blood pressure-to-height ratio (BPHR) has been proposed to facilitate the identification of elevated BP in children. OBJECTIVE: We assessed the performance of BPHR at a single screening visit to identify children with hypertension that is sustained elevated BP. METHOD: In a school-based study conducted in Switzerland, BP was measured at up to three visits in 5207 children. Children had hypertension if BP was elevated at the three visits. Sensitivity, specificity, negative predictive value (NPV), and positive predictive value (PPV) for the identification of hypertension were assessed for different thresholds of BPHR. The ability of BPHR at a single screening visit to discriminate children with and without hypertension was evaluated with receiver operating characteristic (ROC) curve analyses. RESULTS: The prevalence of systolic/diastolic hypertension was 2.2%. Systolic BPHR had a better performance to identify hypertension compared with diastolic BPHR (area under the ROC curve: 0.95 vs. 0.84). The highest performance was obtained with a systolic BPHR threshold set at 0.80 mmHg/cm (sensitivity: 98%; specificity: 85%; PPV: 12%; and NPV: 100%) and a diastolic BPHR threshold set at 0.45 mmHg/cm (sensitivity: 79%; specificity: 70%; PPV: 5%; and NPV: 99%). The PPV was higher among tall or overweight children. CONCLUSION: BPHR at a single screening visit had a high performance to identify hypertension in children, although the low prevalence of hypertension led to a low PPV.

Identification by a monoclonal antibody of a glycolipid highly expressed by cells from the human myeloid lineage.

A monoclonal antibody (MAb) HL-C5, which bound selectively to cells of the myeloid lineage tested, was derived from a fusion between P3/NS2/1-AG8 myeloma cells and splenocytes from a mouse immunized with cells of the promyelocytic leukemia line HL-60. Among a panel of 29 human cell lines derived from either hematopoietic or solid tumors, MAb HL-C5 was found to react exclusively with cells from the five differentiated acute myeloid leukemia lines, HL-60, ML1, ML2, ML3, KG-1B and not with the less differentiated myeloid lines. Fluorescence-activated cell sorter analysis of normal bone marrow samples confirmed that the reactivity of MAb HL-C5 was limited to myeloid cells, from the promyelocytic stage of differentiation to the mature granulocytes. Indirect immunoperoxidase staining of cytocentrifuge preparations of normal bone marrow and peripheral blood leukocytes confirmed these results and showed that MAb HL-C5 stained neutrophils but not eosinophils or basophils. The antigen recognized by HL-C5 was recovered in the upper phase of chloroform-methanol-water lipid extracts prepared from HL-60 cells. By competitive binding experiments, it was found that MAb HL-C5 recognizes the same antigenic determinant as MAb WGHS 29-1, which has been reported to react with glycolipids containing the sugar sequence lacto-N-fucopentaose 111. Autoradiographs of thin layer chromatograms of HL-60 glycolipid extracts which were revealed by incubation with MAb HL-C5 or WGHS 29-1 followed by the addition of 125I-labelled rabbit anti-mouse immunoglobulin antibody confirmed that the two MAbs reacted with the same or structurally very similar glycolipids.

Systematic Identification of High Crash Locations Final Report, 2001

Federal and state policy makers increasingly emphasize the need to reduce highway crash rates. This emphasis is demonstrated in Iowa’s recently released draft Iowa Strategic Highway Safety Plan and by the U.S. Department of Transportation’s placement of “improved transportation safety” at the top of its list of strategic goals. Thus, finding improved methods to enhance highway safety has become a top priority at highway agencies. The objective of this project is to develop tools and procedures by which Iowa engineers can identify potentially hazardous roadway locations and designs, and to demonstrate the utility of these tools by developing candidate lists of high crash locations in the State. An initial task, building an integrated database to facilitate the tools and procedures, is an important product, in and of itself. Accordingly, the Iowa Department of Transportation (Iowa DOT) Geographic Information Management System (GIMS) and Geographic Information System Accident Analysis and Location System (GIS-ALAS) databases were integrated with available digital imagery. (The GIMS database contains roadway characteristics, e.g., lane width, surface and shoulder type, and traffic volume, for all public roadways. GIS-ALAS records include data, e.g., vehicles, drivers, roadway conditions, and the crash severity, for crashes occurring on public roadways during then past 10 years.)

Right coronary artery flow velocity and volume assessment with spiral K-space sampled breathhold velocity-encoded MRI at 3 tesla: accuracy and reproducibility.

PURPOSE: To evaluate accuracy and reproducibility of flow velocity and volume measurements in a phantom and in human coronary arteries using breathhold velocity-encoded (VE) MRI with spiral k-space sampling at 3 Tesla. MATERIALS AND METHODS: Flow velocity assessment was performed using VE MRI with spiral k-space sampling. Accuracy of VE MRI was tested in vitro at five constant flow rates. Reproducibility was investigated in 19 healthy subjects (mean age 25.4 +/- 1.2 years, 11 men) by repeated acquisition in the right coronary artery (RCA). RESULTS: MRI-measured flow rates correlated strongly with volumetric collection (Pearson correlation r = 0.99; P < 0.01). Due to limited sample resolution, VE MRI overestimated the flow rate by 47% on average when nonconstricted region-of-interest segmentation was used. Using constricted region-of-interest segmentation with lumen size equal to ground-truth luminal size, less than 13% error in flow rate was found. In vivo RCA flow velocity assessment was successful in 82% of the applied studies. High interscan, intra- and inter-observer agreement was found for almost all indices describing coronary flow velocity. Reproducibility for repeated acquisitions varied by less than 16% for peak velocity values and by less than 24% for flow volumes. CONCLUSION: 3T breathhold VE MRI with spiral k-space sampling enables accurate and reproducible assessment of RCA flow velocity.

Identification of molecular apocrine breast tumours by microarray analysis.

Previous microarray studies on breast cancer identified multiple tumour classes, of which the most prominent, named luminal and basal, differ in expression of the oestrogen receptor alpha gene (ER). We report here the identification of a group of breast tumours with increased androgen signalling and a 'molecular apocrine' gene expression profile. Tumour samples from 49 patients with large operable or locally advanced breast cancers were tested on Affymetrix U133A gene expression microarrays. Principal components analysis and hierarchical clustering split the tumours into three groups: basal, luminal and a group we call molecular apocrine. All of the molecular apocrine tumours have strong apocrine features on histological examination (P=0.0002). The molecular apocrine group is androgen receptor (AR) positive and contains all of the ER-negative tumours outside the basal group. Kolmogorov-Smirnov testing indicates that oestrogen signalling is most active in the luminal group, and androgen signalling is most active in the molecular apocrine group. ERBB2 amplification is commoner in the molecular apocrine than the other groups. Genes that best split the three groups were identified by Wilcoxon test. Correlation of the average expression profile of these genes in our data with the expression profile of individual tumours in four published breast cancer studies suggest that molecular apocrine tumours represent 8-14% of tumours in these studies. Our data show that it is possible with microarray data to divide mammary tumour cells into three groups based on steroid receptor activity: luminal (ER+ AR+), basal (ER- AR-) and molecular apocrine (ER- AR+).

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

Chroma binary similarity and local alignment applied to cover song identification

We present a new technique for audio signal comparison based on tonal subsequence alignment and its application to detect cover versions (i.e., different performances of the same underlying musical piece). Cover song identification is a task whose popularity has increased in the Music Information Retrieval (MIR) community along in the past, as it provides a direct and objective way to evaluate music similarity algorithms.This article first presents a series of experiments carried outwith two state-of-the-art methods for cover song identification.We have studied several components of these (such as chroma resolution and similarity, transposition, beat tracking or Dynamic Time Warping constraints), in order to discover which characteristics would be desirable for a competitive cover song identifier. After analyzing many cross-validated results, the importance of these characteristics is discussed, and the best-performing ones are finally applied to the newly proposed method. Multipleevaluations of this one confirm a large increase in identificationaccuracy when comparing it with alternative state-of-the-artapproaches.

Identification of the Best Practices for Design, Construction, and Repair of Bridge Approches, 2005

Bridge approach settlement and the formation of the bump is a common problem in Iowa that draws upon considerable resources for maintenance and creates a negative perception in the minds of transportation users. This research study was undertaken to investigate bridge approach problems and develop new concepts for design, construction, and maintenance that will reduce this costly problem. As a result of the research described in this report, the following changes are suggested for implementation on a pilot test basis: • Use porous backfill behind the abutment and/or geocomposite drainage systems to improve drainage capacity and reduce erosion around the abutment. • On a pilot basis, connect the approach slab to the bridge abutment. Change the expansion joint at the bridge to a construction joint of 2 inch. Use a more effective joint sealing system at the CF joint. Change the abutment wall rebar from #5 to #7 for non-integral abutments. • For bridges with soft foundation or embankment soils, implement practices of better compaction, preloading, ground improvement, soil removal and replacement, or soil reinforcement that reduce time-dependent post construction settlements.

Multiuser detection in a dynamic environment — Part II: Joint user identification and parameter estimation

The problem of jointly estimating the number, the identities, and the data of active users in a time-varying multiuser environment was examined in a companion paper (IEEE Trans. Information Theory, vol. 53, no. 9, September 2007), at whose core was the use of the theory of finite random sets on countable spaces. Here we extend that theory to encompass the more general problem of estimating unknown continuous parameters of the active-user signals. This problem is solved here by applying the theory of random finite sets constructed on hybrid spaces. We doso deriving Bayesian recursions that describe the evolution withtime of a posteriori densities of the unknown parameters and data.Unlike in the above cited paper, wherein one could evaluate theexact multiuser set posterior density, here the continuous-parameter Bayesian recursions do not admit closed-form expressions. To circumvent this difficulty, we develop numerical approximationsfor the receivers that are based on Sequential Monte Carlo (SMC)methods (“particle filtering”). Simulation results, referring to acode-divisin multiple-access (CDMA) system, are presented toillustrate the theory.

Radio Frequency Identification Tagging as a Mechanism of Creating a Viable Producer’s Brand in the Cattle Industry, May 2005

This manuscript reports on a project to examine the feasibility of extensive radio frequency identification (RFID) tagging to determine product provenance in the meat production industry. The investigators examined existing technologies and meat production processes as well as emerging technologies in RFID tagging to assess the potential of RFID technologies for provenance assurance. While RFID technologies hold tremendous promise for traceability, the current state of the technology and production process creates challenges for effectively creating full traceability. However, RFID holds tremendous potential for improving processing throughput, which will help make RFIDbased traceability more attractive for adoption by meat processors.

Applause identification and its relevance to archival of Carnatic music

A Carnatic music concert is made up of a sequence of pieces, where each piece corresponds to a particular genre and ra¯aga (melody). Unlike a western music concert, the artist may be applauded intra-performance inter-performance. Most Carnatic music that is archived today correspond to a single audio recordings of entire concerts.The purpose of this paper is to segment single audio recordings into a sequence of pieces using thecharacteristic features of applause and music. Spectral flux, spectral entropy change quite significantly from music to applause and vice-versa. The characteristics of these features for a subset of concerts was studied. A threshold based approach was used to segment the pieces into music fragments and applauses. Preliminary resultson recordings 19 concerts from matched microphones show that the EER is about 17% for a resolution of 0.25 seconds. Further, a parameter called CUSUM is estimatedfor the applause regions. The CUSUM values determine the strength of the applause. The CUSUM is used to characterise the highlights of a concert.

Identification of the Best Practices for Design, Construction, and Repair of Bridge Approaches, January 2005

Bridge approach settlement and the formation of the bump is a common problem in Iowa that draws upon considerable resources for maintenance and creates a negative perception in the minds of transportation users. This research study was undertaken to investigate bridge approach problems and develop new concepts for design, construction, and maintenance that will reduce this costly problem. As a result of the research described in this report, the following changes are suggested for implementation on a pilot test basis: • Use porous backfill behind the abutment and/or geocomposite drainage systems to improve drainage capacity and reduce erosion around the abutment. • On a pilot basis, connect the approach slab to the bridge abutment. Change the expansion joint at the bridge to a construction joint of 2 inch. Use a more effective joint sealing system at the CF joint. Change the abutment wall rebar from #5 to #7 for non-integral abutments. • For bridges with soft foundation or embankment soils, implement practices of better compaction, preloading, ground improvement, soil removal and replacement, or soil reinforcement that reduce time-dependent post construction settlements.

Identification of cancer stem cells in Ewing's sarcoma.

Cancer stem cells that display tumor-initiating properties have recently been identified in several distinct types of malignancies, holding promise for more effective therapeutic strategies. However, evidence of such cells in sarcomas, which include some of the most aggressive and therapy-resistant tumors, has not been shown to date. Here, we identify and characterize cancer stem cells in Ewing's sarcoma family tumors (ESFT), a highly aggressive pediatric malignancy believed to be of mesenchymal stem cell (MSC) origin. Using magnetic bead cell separation of primary ESFT, we have isolated a subpopulation of CD133+ tumor cells that display the capacity to initiate and sustain tumor growth through serial transplantation in nonobese diabetic/severe combined immunodeficiency mice, re-establishing at each in vivo passage the parental tumor phenotype and hierarchical cell organization. Consistent with the plasticity of MSCs, in vitro differentiation assays showed that the CD133+ cell population retained the ability to differentiate along adipogenic, osteogenic, and chondrogenic lineages. Quantitative real-time PCR analysis of genes implicated in stem cell maintenance revealed that CD133+ ESFT cells express significantly higher levels of OCT4 and NANOG than their CD133- counterparts. Taken together, our observations provide the first identification of ESFT cancer stem cells and demonstration of their MSC properties, a critical step towards a better biological understanding and rational therapeutic targeting of these tumors.