Establishment and characterization of a bladder cancer cell line with enhanced doxorubicin resistance by mevalonate pathway activation.

Resistance to chemotherapy is a major problem in the treatment of urothelial bladder cancer. Several mechanisms have been identified in resistance to doxorubicin by analysis of resistant urothelial carcinoma (UC) cell lines, prominently activation of drug efflux pumps and diminished apoptosis. We have derived a new doxorubicin-resistant cell line from BFTC-905 UC cells, designated BFTC-905-DOXO-II. A doxorubicin-responsive green fluorescent protein (GFP) reporter assay indicated that resistance in BFTC-905-DOXO-II was not due to increased drug efflux pump activity, whereas caspase-3/7 activation was indeed diminished. Gene expression microarray analysis revealed changes in proapoptotic and antiapoptotic genes, but additionally induction of the mevalonate (cholesterol) biosynthetic pathway. Treatment with simvastatin restored sensitivity of BFTC-905-DOXO-II to doxorubicin to that of the parental cell line. Induction of the mevalonate pathway has been reported as a mechanism of chemoresistance in other cancers; this is the first observation in bladder cancer. Combinations of statins with doxorubicin-containing chemotherapy regimens may provide a therapeutic advantage in such cases.

Periostin secreted by mesenchymal stem cells supports tendon formation in an ectopic mouse model.

True tendon regeneration in human patients remains a vision of musculoskeletal therapies. In comparison to other mesenchymal lineages the biology of tenogenic differentiation is barely understood. Specifically, easy and efficient protocols are lacking that might enable tendon cell and tissue differentiation based on adult (stem) cell sources. In the murine mesenchymal progenitor cell line C3H10T½, overexpression of the growth factor bone morphogenetic protein 2 (BMP2) and a constitutively active transcription factor, Smad8 L+MH2, mediates tendon cell differentiation in vitro and the formation of tendon-like tissue in vivo. We hypothesized that during this differentiation secreted factors involved in extracellular matrix formation exert a major impact on tendon development. Gene expression analyses revealed four genes encoding secreted factors that are notably upregulated: periostin, C-type lectin domain family 3 (member b), RNase A4, and follistatin-like 1. These factors have not previously been implicated in tendon biology. Among these, periostin showed a specific expression in tenocytes of adult mouse Achilles tendon and in chondrocytes within the nonmineralized fibrocartilage zone of the enthesis with the calcaneus. Overexpression of periostin alone or in combination with constitutively active BMP receptor type in human mesenchymal stem cells and subsequent implantation into ectopic sites in mice demonstrated a reproducible moderate tenogenic capacity that has not been described before. Therefore, periostin may belong to the factors contributing to the development of tenogenic tissue.

The RNA binding protein Musashi1 regulates apoptosis, gene expression and stress granule formation in urothelial carcinoma cells.

The RNA-binding protein Musashi1 (MSI1) is a marker of progenitor cells in the nervous system functioning as a translational repressor. We detected MSI1 mRNA in several bladder carcinoma cell lines, but not in cultured normal uroepithelial cells, whereas the paralogous MSI2 gene was broadly expressed. Knockdown of MSI1 expression by siRNA induced apoptosis and a severe decline in cell numbers in 5637 bladder carcinoma cells. Microarray analysis of gene expression changes after MSI1 knockdown significantly up-regulated 735 genes, but down-regulated only 31. Up-regulated mRNAs contained a highly significantly greater number and density of Musashi binding sites. Therefore, a much larger set of mRNAs may be regulated by Musashi1, which may affect not only their translation, but also their turnover. The study confirmed p21(CIP1) and Numb proteins as targets of Musashi1, suggesting additionally p27(KIP1) in cell-cycle regulation and Jagged-1 in Notch signalling. A significant number of up-regulated genes encoded components of stress granules (SGs), an organelle involved in translational regulation and mRNA turnover, and impacting on apoptosis. Accordingly, heat shock induced SG formation was augmented by Musashi1 down-regulation. Our data show that ectopic MSI1 expression may contribute to tumorigenesis in selected bladder cancers through multiple mechanisms and reveal a previously unrecognized function of Musashi1 in the regulation of SG formation.

Prospectively isolated CD133/CD24-positive ependymal cells from the adult spinal cord and lateral ventricle wall differ in their long-term in vitro self-renewal and in vivo gene expression.

In contrast to ependymal cells located above the subventricular zone (SVZ) of the adult lateral ventricle wall (LVW), adult spinal cord (SC) ependymal cells possess certain neural stem cell characteristics. The molecular basis of this difference is unknown. In this study, antibodies against multiple cell surface markers were applied to isolate pure populations of SC and LVW ependymal cells, which allowed a direct comparison of their in vitro behavior and in vivo gene expression profile. Isolated CD133(+)/CD24(+)/CD45(-)/CD34(-) ependymal cells from the SC displayed in vitro self-renewal and differentiation capacity, whereas those from the LVW did not. SC ependymal cells showed a higher expression of several genes involved in cell division, cell cycle regulation, and chromosome stability, which is consistent with a long-term self-renewal capacity, and shared certain transcripts with neural stem cells of the embryonic forebrain. They also expressed several retinoic acid (RA)-regulated genes and responded to RA exposure. LVW ependymal cells showed higher transcript levels of many genes regulated by transforming growth factor-ß family members. Among them were Dlx2, Id2, Hey1, which together with Foxg1 could explain their potential to turn into neuroblasts under certain environmental conditions.

MethVisual - visualization and exploratory statistical analysis of DNA methylation profiles from bisulfite sequencing.

BACKGROUND: Exploration of DNA methylation and its impact on various regulatory mechanisms has become a very active field of research. Simultaneously there is an arising need for tools to process and analyse the data together with statistical investigation and visualisation. FINDINGS: MethVisual is a new application that enables exploratory analysis and intuitive visualization of DNA methylation data as is typically generated by bisulfite sequencing. The package allows the import of DNA methylation sequences, aligns them and performs quality control comparison. It comprises basic analysis steps as lollipop visualization, co-occurrence display of methylation of neighbouring and distant CpG sites, summary statistics on methylation status, clustering and correspondence analysis. The package has been developed for methylation data but can be also used for other data types for which binary coding can be inferred. The application of the package, as well as a comparison to existing DNA methylation analysis tools and its workflow based on two datasets is presented in this paper. CONCLUSIONS: The R package MethVisual offers various analysis procedures for data that can be binarized, in particular for bisulfite sequenced methylation data. R/Bioconductor has become one of the most important environments for statistical analysis of various types of biological and medical data. Therefore, any data analysis within R that allows the integration of various data types as provided from different technological platforms is convenient. It is the first and so far the only specific package for DNA methylation analysis, in particular for bisulfite sequenced data available in R/Bioconductor enviroment. The package is available for free at http://methvisual.molgen.mpg.de/ and from the Bioconductor Consortium http://www.bioconductor.org.

The genome-wide dynamics of the binding of Ldb1 complexes during erythroid differentiation.

One of the complexes formed by the hematopoietic transcription factor Gata1 is a complex with the Ldb1 (LIM domain-binding protein 1) and Tal1 proteins. It is known to be important for the development and differentiation of the erythroid cell lineage and is thought to be implicated in long-range interactions. Here, the dynamics of the composition of the complex-in particular, the binding of the negative regulators Eto2 and Mtgr1-are studied, in the context of their genome-wide targets. This shows that the complex acts almost exclusively as an activator, binding a very specific combination of sequences, with a positioning relative to transcription start site, depending on the type of the core promoter. The activation is accompanied by a net decrease in the relative binding of Eto2 and Mtgr1. A Chromosome Conformation Capture sequencing (3C-seq) assay also shows that the binding of the Ldb1 complex marks genomic interaction sites in vivo. This establishes the Ldb1 complex as a positive regulator of the final steps of erythroid differentiation that acts through the shedding of negative regulators and the active interaction between regulatory sequences.

Expression profile and transcription factor binding site exploration of imprinted genes in human and mouse.

BACKGROUND: In mammals, imprinted genes are regulated by an epigenetic mechanism that results in parental origin-specific expression. Though allele-specific regulation of imprinted genes has been studied for several individual genes in detail, little is known about their overall tissue-specific expression patterns and interspecies conservation of expression. RESULTS: We performed a computational analysis of microarray expression data of imprinted genes in human and mouse placentae and in a variety of adult tissues. For mouse, early embryonic stages were also included. The analysis reveals that imprinted genes are expressed in a broad spectrum of tissues for both species. Overall, the relative tissue-specific expression levels of orthologous imprinted genes in human and mouse are not highly correlated. However, in both species distinctive expression profiles are found in tissues of the endocrine pathways such as adrenal gland, pituitary, pancreas as well as placenta. In mouse, the placental and embryonic expression patterns of imprinted genes are highly similar. Transcription factor binding site (TFBS) prediction reveals correlation of tissue-specific expression patterns and the presence of distinct TFBS signatures in the upstream region of human imprinted genes. CONCLUSION: Imprinted genes are broadly expressed pre- and postnatally and do not exhibit a distinct overall expression pattern when compared to non-imprinted genes. The relative expression of most orthologous gene pairs varies significantly between human and mouse suggesting rapid species-specific changes in gene regulation. Distinct expression profiles of imprinted genes are confined to certain human and mouse hormone producing tissues, and placentae. In contrast to the overall variability, distinct expression profiles and enriched TFBS signatures are found in human and mouse endocrine tissues and placentae. This points towards an important role played by imprinted gene regulation in these tissues.

Normalization and quantification of differential expression in gene expression microarrays.

Array-based gene expression studies frequently serve to identify genes that are expressed differently under two or more conditions. The actual analysis of the data, however, may be hampered by a number of technical and statistical problems. Possible remedies on the level of computational analysis lie in appropriate preprocessing steps, proper normalization of the data and application of statistical testing procedures in the derivation of differentially expressed genes. This review summarizes methods that are available for these purposes and provides a brief overview of the available software tools.

Up-regulation of glucocorticoid-regulated genes in a mouse model of Rett syndrome.

Rett syndrome (RTT) is a severe form of mental retardation, which is caused by spontaneous mutations in the X-linked gene MECP2. How the loss of MeCP2 function leads to RTT is currently unknown. Mice lacking the Mecp2 gene initially show normal postnatal development but later acquire neurological phenotypes, including heightened anxiety, that resemble RTT. The MECP2 gene encodes a methyl-CpG-binding protein that can act as a transcriptional repressor. Using cDNA microarrays, we found that Mecp2-null animals differentially express several genes that are induced during the stress response by glucocorticoids. Increased levels of mRNAs for serum glucocorticoid-inducible kinase 1 (Sgk) and FK506-binding protein 51 (Fkbp5) were observed before and after onset of neurological symptoms, but plasma glucocorticoid was not significantly elevated in Mecp2-null mice. MeCP2 is bound to the Fkbp5 and Sgk genes in brain and may function as a modulator of glucocorticoid-inducible gene expression. Given the known deleterious effect of glucocorticoid exposure on brain development, our data raise the possibility that disruption of MeCP2-dependent regulation of stress-responsive genes contributes to the symptoms of RTT.

Long-term maintenance of hematopoietic stem cells does not require contact with embryo-derived stromal cells in cocultures.

We recently established that two midgestation-derived stromal clones--UG26-1B6, urogenital ridge-derived, and EL08-1D2, embryonic liver-derived--support the maintenance of murine adult bone marrow and human cord blood hematopoietic repopulating stem cells (HSCs). In this study, we investigate whether direct HSC-stroma contact is required for this stem cell maintenance. Adult bone marrow ckit+ Ly-6C- side population (K6-SP) cells and stromal cells were cocultured under contact or noncontact conditions. These experiments showed that HSCs were maintained for at least 4 weeks in culture and that direct contact between HSCs and stromal cells was not required. To find out which factors might be involved in HSC maintenance, we compared the gene expression profile of EL08-1D2 and UG26-1B6 with four HSC-nonsupportive clones. We found that EL08-1D2 and UG26-1B6 both expressed 21 genes at a higher level, including the putative secreted factors fibroblast growth factor-7, insulin-like growth factor-binding proteins 3 and 4, pleiotrophin, pentaxin-related, and thrombospondin 2, whereas 11 genes, including GPX-3 and HSP27, were expressed at a lower level. In summary, we show for the first time long-term maintenance of adult bone marrow HSCs in stroma noncontact cultures and identify some secreted molecules that may be involved in this support.

Analyzing gene expression time-courses.

Measuring gene expression over time can provide important insights into basic cellular processes. Identifying groups of genes with similar expression time-courses is a crucial first step in the analysis. As biologically relevant groups frequently overlap, due to genes having several distinct roles in those cellular processes, this is a difficult problem for classical clustering methods. We use a mixture model to circumvent this principal problem, with hidden Markov models (HMMs) as effective and flexible components. We show that the ensuing estimation problem can be addressed with additional labeled data-partially supervised learning of mixtures-through a modification of the Expectation-Maximization (EM) algorithm. Good starting points for the mixture estimation are obtained through a modification to Bayesian model merging, which allows us to learn a collection of initial HMMs. We infer groups from mixtures with a simple information-theoretic decoding heuristic, which quantifies the level of ambiguity in group assignment. The effectiveness is shown with high-quality annotation data. As the HMMs we propose capture asynchronous behavior by design, the groups we find are also asynchronous. Synchronous subgroups are obtained from a novel algorithm based on Viterbi paths. We show the suitability of our HMM mixture approach on biological and simulated data and through the favorable comparison with previous approaches. A software implementing the method is freely available under the GPL from http://ghmm.org/gql.

Robust inference of groups in gene expression time-courses using mixtures of HMMs.

MOTIVATION: Genetic regulation of cellular processes is frequently investigated using large-scale gene expression experiments to observe changes in expression over time. This temporal data poses a challenge to classical distance-based clustering methods due to its horizontal dependencies along the time-axis. We propose to use hidden Markov models (HMMs) to explicitly model these time-dependencies. The HMMs are used in a mixture approach that we show to be superior over clustering. Furthermore, mixtures are a more realistic model of the biological reality, as an unambiguous partitioning of genes into clusters of unique functional assignment is impossible. Use of the mixture increases robustness with respect to noise and allows an inference of groups at varying level of assignment ambiguity. A simple approach, partially supervised learning, allows to benefit from prior biological knowledge during the training. Our method allows simultaneous analysis of cyclic and non-cyclic genes and copes well with noise and missing values. RESULTS: We demonstrate biological relevance by detection of phase-specific groupings in HeLa time-course data. A benchmark using simulated data, derived using assumptions independent of those in our method, shows very favorable results compared to the baseline supplied by k-means and two prior approaches implementing model-based clustering. The results stress the benefits of incorporating prior knowledge, whenever available. AVAILABILITY: A software package implementing our method is freely available under the GNU general public license (GPL) at http://ghmm.org/gql

Gene expression changes in the course of neural progenitor cell differentiation.

The molecular changes underlying neural progenitor differentiation are essentially unknown. We applied cDNA microarrays with 13,627 clones to measure dynamic gene expression changes during the in vitro differentiation of neural progenitor cells that were isolated from the subventricular zone of postnatal day 7 mice and grown in vitro as neurospheres. In two experimental series in which we withdrew epidermal growth factor and added the neurotrophins Neurotrophin-4 or BDNF, four time points were investigated: undifferentiated cells grown as neurospheres, and cells 24, 48, and 96 hr after differentiation. Expression changes of selected genes were confirmed by semiquantitative RT-PCR. Ten different groups of gene expression dynamics obtained by cluster analysis are described. To correlate selected gene expression changes to the localization of respective proteins, we performed immunostainings of cultured neurospheres and of brain sections from adult mice. Our results provide new insights into the genetic program of neural progenitor differentiation and give strong hints to as yet unknown cellular communications within the adult subventricular zone stem cell niche.

Different molecular mechanisms underlie placental overgrowth phenotypes caused by interspecies hybridization, cloning, and Esx1 mutation.

To obtain a deeper insight into the genes and gene networks involved in the development of placentopathies, we have assessed global gene expression in three different models of placental hyperplasia caused by interspecies hybridization (IHPD), cloning by nuclear transfer, and mutation of the Esx1 gene, respectively. Comparison of gene expression profiles of approximately 13,000 expressed sequence tags (ESTs) identified specific subsets of genes with changed expression levels in IHPD, cloned, and Esx1 mutant placentas. Of interest, only one gene of known function and one EST of unknown function were found common to all three placentopathies; however, a significant number of ESTs were common to IHPD and cloned placentas. In contrast, only one gene was shared between IHPD and Esx1 mutant, and cloned and Esx1 mutant placentas, respectively. These genes common to different abnormal placental growth genotypes are likely to be important in the occurrence of placentopathy.

Decreased Fas expression in advanced-stage bladder cancer is not related to p53 status.

OBJECTIVES: The Fas-Fas ligand system is an important regulator of apoptosis and is involved in tumor development. Invasive cancers downregulate Fas expression to evade antitumor immune responses. Fas is a transcriptional target of p53, which is often mutated in bladder cancers. Therefore, Fas expression and its relation to p53 mutation was investigated. METHODS: Expression of Fas protein and p53 status was studied by immunohistochemistry in 83 bladder cancer specimens. In addition, mRNA levels for soluble (decoy) and membrane-bound forms of Fas were compared between 10 bladder cancer cell lines and primary uroepithelial cells by quantitative TaqMan polymerase chain reaction. Mutational analysis of the death domain of the Fas gene was performed in all cell lines. RESULTS: Organ-confined tumors maintained specific Fas staining at the cell membrane and often also in the cytoplasm. In higher stage carcinomas, Fas expression became restricted to a smaller fraction of cells or was lacking entirely. The correlation of Fas staining with tumor stage was highly significant but no correlation to tumor grade or survival was found. Furthermore, no statistically significant relationship was observed with either the presence or lack of mutated p53 accumulation. Membrane-bound Fas mRNA was decreased in most, and soluble Fas was increased in all transitional cell carcinoma lines compared with primary uroepithelial cells. No mutations in the death domain were detected. CONCLUSIONS: Fas downregulation occurring in advanced bladder cancer is unrelated to p53 mutations. The results of immunohistochemistry and mRNA studies of soluble and membrane-bound Fas in transitional cell carcinoma lines support the hypothesis of immune evasion in advanced bladder cancer.

Effect of routine repeat transurethral resection for superficial bladder cancer: a long-term observational study.

PURPOSE: We determined the long-term outcome in patients with superficial bladder cancer (Ta and T1) undergoing routine second transurethral bladder tumor resection (ReTURB) in regard to recurrence and progression. MATERIALS AND METHODS: We performed an inception cohort study of 124 consecutive patients with superficial bladder cancer undergoing transurethral resection and routine ReTURB (83) between November 1993 and October 1995 at a German university hospital. Immediately after transurethral resection all lesions were documented on a designed bladder map. ReTURB of the scar from initial resection and other suspicious lesions was performed at a mean of 7 weeks. Patients were followed until recurrence or death, or a minimum of 5 years. RESULTS: Residual tumor was found in 33% of all ReTURB cases, including 27% of Ta and 53% of T1 disease, and in 81% at the initial resection site. Five of the 83 patients underwent radical cystectomy due to ReTURB findings. The estimated risk of recurrence after years 1 to 3 was 18%, 29% and 32%, respectively. After 5 years 63% of the patients undergoing ReTURB were still disease-free (mean recurrence-free survival 62 months, median 87). Progression to muscle invasive disease was observed in only 2 patients (3%) after a mean observation of 61 months. CONCLUSIONS: These data suggest a favorable outcome regarding recurrence and progression in patients with superficial bladder cancer who undergo ReTURB. ReTURB is suggested at least in those at high risk when bladder preservation is intended.

Using hidden Markov models to analyze gene expression time course data.

MOTIVATION: Cellular processes cause changes over time. Observing and measuring those changes over time allows insights into the how and why of regulation. The experimental platform for doing the appropriate large-scale experiments to obtain time-courses of expression levels is provided by microarray technology. However, the proper way of analyzing the resulting time course data is still very much an issue under investigation. The inherent time dependencies in the data suggest that clustering techniques which reflect those dependencies yield improved performance. RESULTS: We propose to use Hidden Markov Models (HMMs) to account for the horizontal dependencies along the time axis in time course data and to cope with the prevalent errors and missing values. The HMMs are used within a model-based clustering framework. We are given a number of clusters, each represented by one Hidden Markov Model from a finite collection encompassing typical qualitative behavior. Then, our method finds in an iterative procedure cluster models and an assignment of data points to these models that maximizes the joint likelihood of clustering and models. Partially supervised learning--adding groups of labeled data to the initial collection of clusters--is supported. A graphical user interface allows querying an expression profile dataset for time course similar to a prototype graphically defined as a sequence of levels and durations. We also propose a heuristic approach to automate determination of the number of clusters. We evaluate the method on published yeast cell cycle and fibroblasts serum response datasets, and compare them, with favorable results, to the autoregressive curves method.

Gene expression profile of mouse bone marrow stromal cells determined by cDNA microarray analysis.

Bone marrow stromal cells (BMSC) have gained increased attention because of their multipotency and adult stem cell character. They have been shown to differentiate into other cell types of the mesenchymal lineage and also into non-mesenchymal cells. The exact identity of the original cells, which are isolated from bone marrow by their selective adherence to plastic, remains unknown to date. We have established and characterized mouse BMSC cultures and analyzed three independent samples by cDNA microarrays. The expression profile was compared with two previous expression studies of human BMSC and revealed a high degree of concordance between different techniques and species. To gain clues about the positional context and biology of the isolated cells within the bone marrow stroma, we searched our data for genes that encode proteins of the extracellular matrix, cell adhesion proteins, cytoskeletal proteins and cytokines/cytokine receptors. This analysis revealed a close association of BMSC with vascular cells and indicated that BMSC resemble pericytes.

Decrease of DNA methyltransferase 1 expression relative to cell proliferation in transitional cell carcinoma.

In many common cancers such as transitional cell carcinoma (TCC), specific genes are hypermethylated, whereas overall DNA methylation is diminished. Genome-wide DNA hypomethylation mostly affects repetitive sequences such as LINE-1 retrotransposons. Methylation of these sequences depends on adequate expression of DNA methyltransferase I (DNMT1) during DNA replication. Therefore, DNMT1 expression relative to proliferation was investigated in TCC cell lines and tissue as well as in renal carcinoma (RCC) cell lines, which also display hypomethylation, as indicated by decreased LINE-1 methylation. Cultured normal uroepithelial cells or normal bladder tissue served as controls. In all tumor cell lines, DNMT1 mRNA as well as protein was decreased relative to the DNA replication factor PCNA, and DNA hypomethylation was present. However, the extents of hypomethylation and DNMT1 downregulation did not correlate. Reporter gene assays showed that the differences in DNMT1 expression between normal and tumor cells were not established at the level of DNMT1 promoter regulation. Diminished DNMT1:PCNA mRNA ratios were also found in 28/45 TCC tissues but did not correlate with the extent of DNA hypomethylation. In addition, expression of the presumed de novo methyltransferases DNMT3A and DNMT3B mRNAs was investigated. DNMT3B overexpression was observed in about half of all high-stage TCC (DNMT3B vs. tumor stage, chi(2): p = 0.03), whereas overexpression of DNMT3A was rarer and less pronounced. Expression of DNMT3A and DNMT3B in most RCC lines was higher than in TCC lines. Our data indicate that DNMT1 expression does not increase adequately with cell proliferation in bladder cancer. This relative downregulation probably contributes to hypomethylation of repetitive DNA but does not determine its extent alone.

Activities of MAP-kinase pathways in normal uroepithelial cells and urothelial carcinoma cell lines.

It is often assumed that MAPK pathways drive proliferation of normal uroepithelial (UEC) and urothelial carcinoma (TCC) cells. To check this assumption, activities and inducibilities of promoters containing serum-response elements (SRE) or AP-1 binding sites were investigated in cultured UEC and seven TCC lines. Reporter plasmids dependent on SRE or AP-1 sites were highly active in UEC, but significantly less so in TCC lines. Reporter activity in TCC lines could be induced by constitutively active MEKK4 or TPA. Accordingly, phosphorylation of the MAPK pathway components MEK, ERK, and ELK1 was most pronounced in UEC and lower in TCC lines. MAPK-dependent promoter activities and bromodeoxyuridine incorporation decreased in UEC upon withdrawal of growth factors, but less so in TCC lines, in which serum diminution increased apoptosis. Likewise, E2F-dependent promoters responded to growth factors in UEC, but were more serum-independent in the TCC lines, which lack either RB1 or p16(INK4A). MEK inhibitors inhibited BrdU incorporation in UEC more strongly than in TCC lines. Thus, proliferation of normal uroepithelial cells is indeed associated with activation of MAPK pathways. However, autonomous proliferation of TCC lines--unexpectedly--appears much less dependent on MAPK activation and may rather be promoted by defects in cell cycle regulation.