Cancer stem cells (CSCs) have been identified in a number of solid tumors, but not yet in rhabdomyosarcoma (RMS), the most frequently occurring soft tissue tumor in childhood. Hence, the aim of this study was to identify and characterize a CSC population in RMS using a functional approach. We found that embryonal rhabdomyosarcoma (eRMS) cell lines can form rhabdomyosarcoma spheres (short rhabdospheres) in stem cell medium containing defined growth factors over several passages. Using an orthotopic xenograft model, we demonstrate that a 100 fold less sphere cells result in faster tumor growth compared to the adherent population suggesting that CSCs were enriched in the sphere population. Furthermore, stem cell genes such as oct4, nanog, c-myc, pax3 and sox2 are significantly upregulated in rhabdospheres which can be differentiated into multiple lineages such as adipocytes, myocytes and neuronal cells. Surprisingly, gene expression profiles indicate that rhabdospheres show more similarities with neuronal than with hematopoietic or mesenchymal stem cells. Analysis of these profiles identified the known CSC marker CD133 as one of the genes upregulated in rhabdospheres, both on RNA and protein levels. CD133(+) sorted cells were subsequently shown to be more tumorigenic and more resistant to commonly used chemotherapeutics. Using a tissue microarray (TMA) of eRMS patients, we found that high expression of CD133 correlates with poor overall survival. Hence, CD133 could be a prognostic marker for eRMS. These experiments indicate that a CD133(+) CSC population can be enriched from eRMS which might help to develop novel targeted therapies against this pediatric tumor.

Salmonella Typhimurium (S. Tm) is a leading cause of diarrhea. The disease is triggered by pathogen invasion into the gut epithelium. Invasion is attributed to the SPI-1 type 3 secretion system (T1). T1 injects effector proteins into epithelial cells and thereby elicits rearrangements of the host cellular actin cytoskeleton and pathogen invasion. The T1 effector proteins SopE, SopB, SopE2 and SipA are contributing to this. However, the host cell factors contributing to invasion are still not completely understood. To address this question comprehensively, we used Hela tissue culture cells, a genome-wide siRNA library, a modified gentamicin protection assay and S. TmSipA, a sopBsopE2sopE mutant which strongly relies on the T1 effector protein SipA to invade host cells. We found that S. TmSipA invasion does not elicit membrane ruffles, nor promote the entry of non-invasive bacteria "in trans". However, SipA-mediated infection involved the SPIRE family of actin nucleators, besides well-established host cell factors (WRC, ARP2/3, RhoGTPases, COPI). Stage-specific follow-up assays and knockout fibroblasts indicated that SPIRE1 and SPIRE2 are involved in different steps of the S. Tm infection process. Whereas SPIRE1 interferes with bacterial binding, SPIRE2 influences intracellular replication of S. Tm. Hence, these two proteins might fulfill non-redundant functions in the pathogen-host interaction. The lack of co-localization hints to a short, direct interaction between S. Tm and SPIRE proteins or to an indirect effect.

Recent advancements of sequencing technology have opened up unprecedented opportunities in many application areas. Virus samples can now be sequenced efficiently with very deep coverage to infer the genetic diversity of the underlying virus populations. Several sequencing platforms with different underlying technologies and performance characteristics are available for viral diversity studies. Here, we investigate how the differences between two common platforms provided by 454/Roche and Illumina affect viral diversity estimation and the reconstruction of viral haplotypes. Using a mixture of ten HIV clones sequenced with both platforms and additional simulation experiments, we assessed the trade-off between sequencing coverage, read length, and error rate. For fixed costs, short Illumina reads can be generated at higher coverage and allow for detecting variants at lower frequencies. They can also be sufficient to assess the diversity of the sample if sequences are dissimilar enough, but, in general, assembly of full-length haplotypes is feasible only with the longer 454/Roche reads. The quantitative comparison highlights the advantages and disadvantages of both platforms and provides guidance for the design of viral diversity studies.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

MiRNAs and other small noncoding RNAs (sncRNAs) are key players in post-transcriptional gene regulation. HIV-1 derived small noncoding RNAs (sncRNAs) have been described in HIV-1 infected cells, but their biological functions still remain to be elucidated. Here, we approached the question whether viral sncRNAs may play a role in the RNA interference (RNAi) pathway or whether viral mRNAs are targeted by cellular miRNAs in human monocyte derived macrophages (MDM).

γ-radiation is an effective treatment for cancer. There is evidence that radiotherapy supports tumor-specific immunity. It was described that irradiation induces de novo protein synthesis and enhances antigen presentation, we therefore investigated whether γ-radiation results in increased expression of cancer-testis (CT) antigens and MHC-I, thus allowing efficient immunological control. This is relevant because the expression of CT-antigens and MHC-I on tumor cells is often heterogeneous. We found that the changes induced by γ-radiation promote the immunological recognition of the tumor, which is illustrated by the increased infiltration by lymphocytes after radiotherapy.

Recent molecular studies provide evidence for a significant transcriptional plasticity of tumor cell subpopulations that facilitate an active contribution to tumor vasculature. This feature is accompanied by morphological changes both in vitro and in vivo. Herein, we investigated the morphological plasticity of tumor cells with special focus on vasculogenic mimicry and neovascularisation in human melanoma and mouse xenografts of human melanoma cell lines. In melanoma xenograft experiments, different vessel markers and green fluorescent protein expression were used to show how melanoma cells contribute to neovascularization. Additionally, we analyzed neovascularization in 49 primary melanomas and 175 melanoma metastases using immunostaining for blood (CD34) and lymphatic (D2-40) vessel-specific markers. We found significantly more lymphatic vessels in primary melanomas than in melanoma metastases (p<0.0001). In contrast to the near absence of lymphatic vessels within metastases, we found extensive blood micro-neovascularization. Blood micro-neovascularization was absent in micro metastases (less than 2 mm). A significant inverse correlation between Glut-1 expression (implying local hypoxia) and the presence of microvessels indicates their functional activity as blood vessels (p<0.0001). We suggest that the hypoxic microenvironment in metastases contributes to a phenotype switch allowing melanoma cells to physically contribute to blood vessel formation.

To evaluate the frequency of MRE11/RAD50/NBS1 (MRN)-complex loss of protein expression in endometrial cancers (EC) and to determine whether loss of MRE11 renders the cancer cells sensitive to Poly(ADP-ribose) polymerase (PARP)-inhibitory treatment.

Identifying oral cancer lesions associated with high risk of relapse and predicting clinical outcome remain challenging questions in clinical practice. Genomic alterations may add prognostic information and indicate biological aggressiveness thereby emphasizing the need for genome-wide profiling of oral cancers. High-resolution array comparative genomic hybridization was performed to delineate the genomic alterations in clinically annotated primary gingivo-buccal complex and tongue cancers (n = 60). The specific genomic alterations so identified were evaluated for their potential clinical relevance. Copy-number changes were observed on chromosomal arms with most frequent gains on 3q (60%), 5p (50%), 7p (50%), 8q (73%), 11q13 (47%), 14q11.2 (47%), and 19p13.3 (58%) and losses on 3p14.2 (55%) and 8p (83%). Univariate statistical analysis with correction for multiple testing revealed chromosomal gain of region 11q22.1-q22.2 and losses of 17p13.3 and 11q23-q25 to be associated with loco-regional recurrence (P = 0.004, P = 0.003, and P = 0.0003) and shorter survival (P = 0.009, P = 0.003, and P 0.0001) respectively. The gain of 11q22 and loss of 11q23-q25 were validated by interphase fluorescent in situ hybridization (I-FISH). This study identifies a tractable number of genomic alterations with few underlying genes that may potentially be utilized as biological markers for prognosis and treatment decisions in oral cancers.

The chemoprotective properties of sulforaphane (SF), derived from cruciferous vegetables, are widely acknowledged to arise from its potent induction of xenobiotic-metabolizing and antioxidant enzymes. However, much less is known about the impact of SF on the efficacy of cancer therapy through the modulation of drug-metabolizing enzymes. To identify proteins modulated by a low concentration of SF, we treated HT29 colon cancer cells with 2.5 μM SF. Protein abundance changes were detected by stable isotope labeling of amino acids in cell culture. Among 18 proteins found to be significantly up-regulated, aldo-keto reductase 1C3 (AKR1C3), bioactivating the DNA cross-linking prodrug PR-104A, was further characterized. Preconditioning HT29 cells with SF reduced the EC50 of PR-104A 3.6-fold. The increase in PR-104A cytotoxicity was linked to AKR1C3 abundance and activity, both induced by SF in a dose-dependent manner. This effect was reproducible in a second colon cancer cell line, SW620, but not in other colon cancer cell lines where AKR1C3 abundance and activity were absent or barely detectable and could not be induced by SF. Interestingly, SF had no significant influence on PR-104A cytotoxicity in non-cancerous, immortalized human colonic epithelial cell lines expressing either low or high levels of AKR1C3. In conclusion, the enhanced response of PR-104A after preconditioning with SF was apparent only in cancer cells provided that AKR1C3 is expressed, while its expression in non-cancerous cells did not elicit such a response. Therefore, a subset of cancers may be susceptible to combined food-derived component and prodrug treatments with no harm to normal tissues.

The treatment of complex diseases often relies on combinatorial therapy, a strategy where drugs are used to target multiple genes simultaneously. Promising candidate genes for combinatorial perturbation often constitute epistatic genes, i.e., genes which contribute to a phenotype in a non-linear fashion. Experimental identification of the full landscape of genetic interactions by perturbing all gene combinations is prohibitive due to the exponential growth of testable hypotheses. Here we present a model for the inference of pairwise epistatic, including synthetic lethal, gene interactions from siRNA-based perturbation screens. The model exploits the combinatorial nature of siRNA-based screens resulting from the high numbers of sequence-dependent off-target effects, where each siRNA apart from its intended target knocks down hundreds of additional genes. We show that conditional and marginal epistasis can be estimated as interaction coefficients of regression models on perturbation data. We compare two methods, namely glinternet and xyz, for selecting non-zero effects in high dimensions as components of the model, and make recommendations for the appropriate use of each. For data simulated from real RNAi screening libraries, we show that glinternet successfully identifies epistatic gene pairs with high accuracy across a wide range of relevant parameters for the signal-to-noise ratio of observed phenotypes, the effect size of epistasis and the number of observations per double knockdown. xyz is also able to identify interactions from lower dimensional data sets (fewer genes), but is less accurate for many dimensions. Higher accuracy of glinternet, however, comes at the cost of longer running time compared to xyz. The general model is widely applicable and allows mining the wealth of publicly available RNAi screening data for the estimation of epistatic interactions between genes. As a proof of concept, we apply the model to search for interactions, and potential targets for treatment, among previously published sets of siRNA perturbation screens on various pathogens. The identified interactions include both known epistatic interactions as well as novel findings.

MAGE-C1/CT7 and MAGE-C2/CT10 are members of the large MAGE family of cancer-testis (CT) antigens. CT antigens are promising targets for immunotherapy in cancer because their expression is restricted to cancer and germ line cells and a proportion of cancer patients presents with immune responses against CT antigens, which clearly demonstrates their immunogenicity. This study investigates the expression of MAGE-C1/CT7 and MAGE-C2/CT10 in primary and metastatic melanoma. Immunohistochemical staining of tissue microarrays that consisted of 59 primary malignant melanomas of the skin, 163 lymph node and distant melanoma metastases and 68 melanoma cell lines was performed. We found MAGE-C1/CT7 expression in 15 out of 50 (24%) primary melanomas and 15 out of 50 (24%) cell lines, whereas MAGE-C2/CT10 was detected in 17 out of 51 (33%) primary melanomas and 14 out of 68 (17%) cell lines. MAGE-C1/CT7 and MAGE-C2/CT10 were both detected in 40% of melanoma metastases. Patients with MAGE-C1/CT7 or MAGE-C2/CT10 positive primary melanoma had significantly more lymph node metastases (p = 0.005 and p<0.001, resp.). Prediction of lymph node metastasis by MAGE-C1/CT7 and MAGE-C2/CT10 was independent of tumor cell proliferation rate (Ki67 labeling index) in a multivariate analysis (p = 0.01). Our results suggest that the expression of MAGE-C1/CT7 and MAGE-C2/CT10 in primary melanoma is a potent predictor of sentinel lymph node metastasis.

Ipilimumab and tremelimumab are human monoclonal antibodies (Abs) against cytotoxic T-lymphocyte antigen-4 (CTLA-4). Ipilimumab was the first agent to show a statistically significant benefit in overall survival in advanced melanoma patients. Currently, there is no proven association between the BRAFV600 mutation and the disease control rate in response to ipilimumab. This analysis was carried out to assess if BRAFV600 and NRAS mutation status affects the clinical outcome of anti-CTLA-4-treated melanoma patients. This is a retrospective multi-center analysis of 101 patients, with confirmed BRAF and NRAS mutation status, treated with anti-CTLA-4 antibodies from December 2006 until August 2012. The median overall survival, defined from the treatment start date with the anti-CTLA-4. Abs-treatment to death or till last follow up, of BRAFV600 or NRAS mutant patients (n = 62) was 10.12 months (95% CI 6.78-13.2) compared to 8.26 months (95% CI 6.02-19.9) in BRAFV600/NRASwt subpopulation (n = 39) (p = 0.67). The median OS of NRAS mutated patients (n = 24) was 12.1 months and although was prolonged compared to the median OS of BRAF mutated patients (n = 38, mOS = 8.03 months) or BRAFV600/NRASwt patients (n = 39, mOS = 8.26 months) the difference didn't reach statistical significance (p = 0.56). 69 patients were able to complete 4 cycles of anti-CTLA-4 treatment. Of the 24 patients treated with selected BRAF- or MEK-inhibitors, 16 patients received anti-CTLA 4 Abs following either a BRAF or MEK inhibitor with only 8 of them being able to finish 4 cycles of treatment. Based on our results, there is no difference in the median OS in patients treated with anti-CTLA-4 Abs implying that the BRAF/NRAS mutation status alone is not sufficient to predict the outcome of patients treated with anti-CTLA-4 Abs.