Melanoma progression requires deregulation of gene expression by currently uncharacterized epigenetic mechanisms. A mouse model based on changes in cell microenvironment was developed by our group to study melanocyte malignant transformation. Melanoma cell lines (4C11− and 4C11+) were obtained as result of 5 sequential anchorage blockades of non-tumorigenic melan-a melanocytes. Melan-a cells submitted to 4 de-adhesion cycles were also established (4C), are non-tumorigenic and represent an intermediary phase of tumor progression. The aim of this work was to identify factors contributing to epigenetic modifications in early and later phases of malignant transformation induced by anchorage impediment. Epigenetic alterations occur early in tumorigenesis; 4C cell line shows changes in global and gene-specific DNA methylation and histone marks. Many histone modifications differ between melan-a, 4C, 4C11− (non-metastatic melanoma cell line) and 4C11+ (metastatic melanoma cell line) which could be associated with changes in gene and microRNA expression. These epigenetic alterations seem to play a key role in malignant transformation since melanocytes treated with 5-Aza-2′-deoxycytidine before each anchorage blockade do not transform. Some epigenetic changes seem to be also responsible for the maintenance of malignant phenotype, since melanoma cell lines (4C11− and 4C11+) treated in vitro with 5-Aza-2′-deoxycytidine or Trichostatin A showed reduction of tumor growth in vivo. Changes in gene expression reflecting cell adaptation to new environment were also observed. We propose a model in which sustained microenvironmental stress in melanocytes results in epigenetic reprogramming. Thus, after adaptation, cells may acquire epigenetic marks that could contribute to the establishment of a malignant phenotype.

Chondrosarcomas are among the most malignant skeletal tumors. Dedifferentiated chondrosarcoma is a highly aggressive subtype of chondrosarcoma, with lung metastases developing within a few months of diagnosis in 90% of patients. In this paper we performed comparative analyses of the transcriptomes of five individual metastatic lung lesions that were surgically resected from a patient with dedifferentiated chondrosarcoma. We document for the first time a high heterogeneity of gene expression profiles among the individual lung metastases. Moreover, we reveal a signature of “multifunctional” genes that are expressed in all metastatic lung lesions. Also, for the first time, we document the occurrence of massive macrophage infiltration in dedifferentiated chondrosarcoma lung metastases.

Background

De novo retrotransposition of Alu elements has been recognized as a major driver for insertion polymorphisms in human populations. In this study, we exploited Alu-anchored bisulfite PCR libraries to identify evolutionarily recent Alu element insertions, and to investigate their genetic and epigenetic variation.

Results

A total of 327 putatively recent Alu insertions were identified, altogether represented by 1,762 sequence reads. Nearly all such de novo retrotransposition events (316/327) were novel. Forty-seven out of forty-nine randomly selected events, corresponding to nineteen genomic loci, were sequence-verified. Alu element insertions remained hemizygous in one or more individuals in sixteen of the nineteen genomic loci. The Alu elements were found to be enriched for young Alu families with characteristic sequence features, such as the presence of a longer poly(A) tail. In addition, we documented the occurrence of a duplication of the AT-rich target site in their immediate flanking sequences, a hallmark of retrotransposition. Furthermore, we found the sequence motif (TT/AAAA) that is recognized by the ORF2P protein encoded by LINE-1 in their 5'-flanking regions, consistent with the fact that Alu retrotransposition is facilitated by LINE-1 elements. While most of these Alu elements were heavily methylated, we identified an Alu localized 1.5 kb downstream of TOMM5 that exhibited a completely unmethylated left arm. Interestingly, we observed differential methylation of its immediate 5' and 3' flanking CpG dinucleotides, in concordance with the unmethylated and methylated statuses of its internal 5' and 3' sequences, respectively. Importantly, TOMM5's CpG island and the 3 Alu repeats and 1 MIR element localized upstream of this newly inserted Alu were also found to be unmethylated. Methylation analyses of two additional genomic loci revealed no methylation differences in CpG dinucleotides flanking the Alu insertion sites in the two homologous chromosomes, irrespective of the presence or absence of the insertion.

Conclusions

We anticipate that the combination of methodologies utilized in this study, which included repeat-anchored bisulfite PCR sequencing and the computational analysis pipeline herein reported, will prove invaluable for the generation of genetic and epigenetic variation maps.

Studies are beginning to emerge that demonstrate intriguing differences between human induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs). Here, we investigated the expression of key members of the Nodal embryonic signaling pathway, critical to the maintenance of pluripotency in hESCs. Western blot and Real-time RT-PCR analyses reveal slightly lower levels of Nodal (a TGF-β family member) and Cripto-1 (Nodal’s co-receptor) and a dramatic decrease in Lefty (Nodal’s inhibitor and TGF-β family member) in hiPSCs compared with hESCs. The noteworthy drop in hiPSC’s Lefty expression correlated with an increase in the methylation of Lefty B CpG island. Based on these findings, we addressed a more fundamental question related to the consequences of epigenetically reprogramming hiPSCs, especially with respect to maintaining a stable ESC phenotype. A global comparative analysis of 365 microRNAs (miRs) in two hiPSC vs. four hESC lines ultimately identified 10 highly expressed miRs in hiPCSs with >10-fold difference, which have been shown to be cancer related. These data demonstrate cancer hallmarks expressed by hiPSCs, which will require further assessment for their impact on future therapies.

The cystic fibrosis transmembrane conductance regulator (CFTR) gene is driven by a promoter that cannot alone account for the temporal and tissue-specific regulation of the gene. This has led to the search for additional regulatory elements that cooperate with the basal promoter to achieve coordinated expression. We previously identified two alternative upstream exons of the gene that were mutually exclusive of the first exon, and one of which showed temporal regulation in the human and sheep lung. We now demonstrate that this alternative splice product generates a stable protein, which initiates translation at an ATG in exon 4, and thus lacks the N terminus of CFTR. The other splice variant inhibits translation of the protein. In a search for the promoter used by the upstream exons, we identified a novel element that contributes to the activity of the basal CFTR promoter in airway epithelial cells, but does not function independently. Finally, we demonstrate that, in primary airway cells, skin fibroblasts, and both airway and intestinal cell lines, the CFTR promoter is unmethylated, irrespective of CFTR expression status. Thus, methylation is not the main cause of inactivation of CFTR transcription.

A limited number of reports have investigated the role of microRNAs in osteosarcoma. In this study, we performed miRNA expression profiling of osteosarcoma cell lines, tumor samples, and normal human osteoblasts. Twenty-two differentially expressed microRNAs were identified using high throughput real-time PCR analysis, and 4 (miR-135b, miR-150, miR-542-5p, and miR-652) were confirmed and validated in a different group of tumors. Both miR-135b and miR-150 have been previously shown to be important in cancer. We hypothesize that dysregulation of differentially expressed microRNAs may contribute to tumorigenesis. They might also represent molecular biomarkers or targets for drug development in osteosarcoma.

Context

Glioblastomas—uniformly fatal brain tumors—often have both monosomy of chromosome 10 and gains of the epidermal growth factor receptor (EGFR) gene locus on chromosome 7, an association for which the mechanism is poorly understood.

Objectives

To assess whether coselection of EGFR gains on 7p12 and monosomy 10 in glioblastomas promotes tumorigenic epidermal growth factor (EGF) signaling through loss of the annexin A7 (ANXA7) gene on 10q21.1–q21.2 and whether ANXA7 acts as a tumor suppressor gene by regulating EGFR in glioblastomas.

Design, Setting, and Patients

Multidimensional analysis of gene, coding sequence, promoter methylation, messenger RNA (mRNA) transcript, protein data for ANXA7 (and EGFR), and clinical patient data profiles of 543 high-grade gliomas from US medical centers and The Cancer Genome Atlas pilot project (made public 2006–2008; and unpublished, tumors collected 2001–2008). Functional analyses using LN229 and U87 glioblastoma cells.

Main Outcome Measures

Associations among ANXA7 gene dosage, coding sequence, promoter methylation, mRNA transcript, and protein expression. Effect of ANXA7 haploinsufficiency on EGFR signaling and patient survival. Joint effects of loss of ANXA7 and gain of EGFR expression on tumorigenesis.

Results

Heterozygous ANXA7 gene deletion is associated with significant loss of ANXA7 mRNA transcript expression (P=1×10−15; linear regression) and a reduction (mean [SEM]) of 91.5% (2.3%) of ANXA7 protein expression compared with ANXA7 wild-type glioblastomas (P=.004; unpaired t test). ANXA7 loss of function stabilizes the EGFR protein (72%–744% increase in EGFR protein abundance) and augments EGFR transforming signaling in glioblastoma cells. ANXA7 haploinsufficiency doubles tumorigenic potential of glioblastoma cells, and combined ANXA7 knockdown and EGFR overexpression promotes tumorigenicity synergistically. The heterozygous loss of ANXA7 in≈75% of glioblastomas in the The Cancer Genome Atlas plus infrequency of ANXA7 mutation (≈6% of tumors) indicates its role as a haploinsufficiency gene. ANXA7 mRNA transcript expression, dichotomized at the median, associates with patient survival in 191 glioblastomas (log-rank P=.008; hazard ratio [HR], 0.667; 95% confidence interval [CI], 0.493–0.902; 46.9 vs 74.8 deaths/100 person-years for high vs low ANXA7 mRNA expression) and with a separate group of 180 high-grade gliomas (log-rank P=.00003; HR, 0.476; 95% CI, 0.333–0.680; 21.8 vs 50.0 deaths/100 person-years for high vs low ANXA7 mRNA expression). Deletion of the ANXA7 gene associates with poor patient survival in 189 glioblastomas (log-rank P=.042; HR, 0.686; 95% CI, 0.476–0.989; 54.0 vs 80.1 deaths/100 person-years for wild-type ANXA7 vs ANXA7 deletion).

Conclusion

Haploinsufficiency of the tumor suppressor ANXA7 due to monosomy of chromosome 10 provides a clinically relevant mechanism to augment EGFR signaling in glioblastomas beyond that resulting from amplification of the EGFR gene.

Purpose

To understand the changes in gene expression in PV progenitor cells and their relationship to JAK2V617F

Experimental Design

mRNA isolated from CD34+ cells from 9 PV patients and normal controls was profiled using Affymetrix arrays. Gene expression change mediated by JAK2V617F was determined by profiling CD34+ cells transduced with the kinase and by analysis of leukemia cell lines harboring JAK2V617F, treated with an inhibitor.

Results

A PV expression signature was enriched for genes involved in hematopoietic development, inflammatory responses and cell proliferation. By quantitative RT-PCR, 23 genes were consistently deregulated in all patient samples. Several of these genes such as WT1 and KLF4 were regulated by JAK2, while others such as NFIB and EVI1 appeared to be deregulated in PV by a JAK2 independent mechanism. Using cell line models and comparing gene expression profiles of cell lines and PV CD34+ PV specimens, we have identified panels of JAK2 dependent 14 genes and JAK2 independent 12 genes. These two 14 and 12 gene sets could separate not only PV from normal CD34+ specimens, but also other MPN such as ET and MF from their normal counterparts.

Conclusions

A subset of the aberrant gene expression in PV progenitor cells can be attributed to the action of the mutant kinase, but there remain a significant number of genes characteristic of the disease but deregulated by as yet unknown mechanisms. Genes deregulated in PV as a result of the action of JAK2V617F or independent of the kinase may represent other targets for therapy.

Genomic DNA methylation contributes substantively to transcriptional regulations that underlie mammalian development and cellular differentiation. Much effort has been made to decipher the molecular mechanisms governing the establishment and maintenance of DNA methylation patterns. However, little is known about genome-wide variation of DNA methylation patterns. In this study, we introduced the concept of methylation entropy, a measure of the randomness of DNA methylation patterns in a cell population, and exploited it to assess the variability in DNA methylation patterns of Alu repeats and promoters. A few interesting observations were made: (i) within a cell population, methylation entropy varies among genomic loci; (ii) among cell populations, the methylation entropies of most genomic loci remain constant; (iii) compared to normal tissue controls, some tumors exhibit greater methylation entropies; (iv) Alu elements with high methylation entropy are associated with high GC content but depletion of CpG dinucleotides and (v) Alu elements in the intronic regions or far from CpG islands are associated with low methylation entropy. We further identified 12 putative allelic-specific methylated genomic loci, including four Alu elements and eight promoters. Lastly, using subcloned normal fibroblast cells, we demonstrated the highly variable methylation patterns are resulted from low fidelity of DNA methylation inheritance.

Trypanosoma congolense is one of the most economically important pathogens of livestock in Africa. Culture-derived parasites of each of the three main insect stages of the T. congolense life cycle, i.e., the procyclic, epimastigote and metacyclic stages, and bloodstream stage parasites isolated from infected mice, were used to construct stage-specific cDNA libraries and expressed sequence tags (ESTs or cDNA clones) in each library were sequenced. Thirteen EST clusters encoding different variant surface glycoproteins (VSGs) were detected in the metacyclic library and twenty-six VSG EST clusters were found in the bloodstream library, six of which are shared by the metacyclic library. Rare VSG ESTs are present in the epimastigote library, and none were detected in the procyclic library. ESTs encoding enzymes that catalyze oxidative phosphorylation and amino acid metabolism are about twice as abundant in the procyclic and epimastigote stages as in the metacyclic and bloodstream stages. In contrast, ESTs encoding enzymes involved in glycolysis, the citric acid cycle and nucleotide metabolism are about the same in all four developmental stages. Cysteine proteases, kinases and phosphatases are the most abundant enzyme groups represented by the ESTs. All four libraries contain T. congolense-specific expressed sequences not present in the T. brucei and T. cruzi genomes. Normalized cDNA libraries were constructed from the metacyclic and bloodstream stages, and found to be further enriched for T. congolense-specific ESTs. Given that cultured T. congolense offers an experimental advantage over other African trypanosome species, these ESTs provide a basis for further investigation of the molecular properties of these four developmental stages, especially the epimastigote and metacyclic stages for which it is difficult to obtain large quantities of organisms. The T. congolense EST databases are available at: http://www.sanger.ac.uk/Projects/T_congolense/EST_index.shtml.

We have previously shown that the microenvironment of human embryonic stem cells (hESCs) is able to change and reprogram aggressive cancer cells to a less aggressive state. Some mechanisms implicated in the phenotypic changes observed after this exposure are mainly associated with the Nodal signaling pathway, which plays a key role in tumor cell plasticity. However, several other molecular mechanisms might be related directly and/or indirectly to these changes, including microRNA (miRNA) regulation and DNA methylation.

Aim:

To further explore the epigenetic mechanisms potentially underlying the phenotypic changes that occur after exposing metastatic melanoma cells to a hESC microenvironment.

Materials & Methods:

A total of 365 miRNAs were screened using the TaqMan® Low Density Arrays. We also evaluated whether DNA methylation could be one of the factors regulating the expression of the inhibitor of Nodal, Lefty, in hESCs (where it is highly expressed) vs melanoma cells (where it is not expressed).

Results:

Using these experimental approaches, we identified miRNAs that are up- and down-regulated in melanoma cells exposed to a hESC microenvironment, such as miR-302a and miR-27b, respectively. We also demonstrate that Notch4 is one of the targets of miR-302a, which is upstream of Nodal. Additionally, one of the mechanisms that might explain the absence of the inhibitor of Nodal, Lefty, in cancer cells is silencing by DNA methylation, which provides new insights into the unregulated expression of Nodal in melanoma.

Conclusion:

These findings suggest that epigenetic changes such as DNA methylation and regulation by microRNAs might play a significant role in tumor cell plasticity and the metastatic phenotype.

Background

The medicinal leech, Hirudo medicinalis, is an important model system for the study of nervous system structure, function, development, regeneration and repair. It is also a unique species in being presently approved for use in medical procedures, such as clearing of pooled blood following certain surgical procedures. It is a current, and potentially also future, source of medically useful molecular factors, such as anticoagulants and antibacterial peptides, which may have evolved as a result of its parasitizing large mammals, including humans. Despite the broad focus of research on this system, little has been done at the genomic or transcriptomic levels and there is a paucity of openly available sequence data. To begin to address this problem, we constructed whole embryo and adult central nervous system (CNS) EST libraries and created a clustered sequence database of the Hirudo transcriptome that is available to the scientific community.

Results

A total of ~133,000 EST clones from two directionally-cloned cDNA libraries, one constructed from mRNA derived from whole embryos at several developmental stages and the other from adult CNS cords, were sequenced in one or both directions by three different groups: Genoscope (French National Sequencing Center), the University of Iowa Sequencing Facility and the DOE Joint Genome Institute. These were assembled using the phrap software package into 31,232 unique contigs and singletons, with an average length of 827 nt. The assembled transcripts were then translated in all six frames and compared to proteins in NCBI's non-redundant (NR) and to the Gene Ontology (GO) protein sequence databases, resulting in 15,565 matches to 11,236 proteins in NR and 13,935 matches to 8,073 proteins in GO. Searching the database for transcripts of genes homologous to those thought to be involved in the innate immune responses of vertebrates and other invertebrates yielded a set of nearly one hundred evolutionarily conserved sequences, representing all known pathways involved in these important functions.

Conclusions

The sequences obtained for Hirudo transcripts represent the first major database of genes expressed in this important model system. Comparison of translated open reading frames (ORFs) with the other openly available leech datasets, the genome and transcriptome of Helobdella robusta, shows an average identity at the amino acid level of 58% in matched sequences. Interestingly, comparison with other available Lophotrochozoans shows similar high levels of amino acid identity, where sequences match, for example, 64% with Capitella capitata (a polychaete) and 56% with Aplysia californica (a mollusk), as well as 58% with Schistosoma mansoni (a platyhelminth). Phylogenetic comparisons of putative Hirudo innate immune response genes present within the Hirudo transcriptome database herein described show a strong resemblance to the corresponding mammalian genes, indicating that this important physiological response may have older origins than what has been previously proposed.

Background

Blood feeding evolved independently in worms, arthropods and mammals. Among the adaptations to this peculiar diet, these animals developed an armament of salivary molecules that disarm their host's anti-bleeding defenses (hemostasis), inflammatory and immune reactions. Recent sialotranscriptome analyses (from the Greek sialo = saliva) of blood feeding insects and ticks have revealed that the saliva contains hundreds of polypeptides, many unique to their genus or family. Adult tsetse flies feed exclusively on vertebrate blood and are important vectors of human and animal diseases. Thus far, only limited information exists regarding the Glossina sialome, or any other fly belonging to the Hippoboscidae.

Results

As part of the effort to sequence the genome of Glossina morsitans morsitans, several organ specific, high quality normalized cDNA libraries have been constructed, from which over 20,000 ESTs from an adult salivary gland library were sequenced. These ESTs have been assembled using previously described ESTs from the fat body and midgut libraries of the same fly, thus totaling 62,251 ESTs, which have been assembled into 16,743 clusters (8,506 of which had one or more EST from the salivary gland library). Coding sequences were obtained for 2,509 novel proteins, 1,792 of which had at least one EST expressed in the salivary glands. Despite library normalization, 59 transcripts were overrepresented in the salivary library indicating high levels of expression. This work presents a detailed analysis of the salivary protein families identified. Protein expression was confirmed by 2D gel electrophoresis, enzymatic digestion and mass spectrometry. Concurrently, an initial attempt to determine the immunogenic properties of selected salivary proteins was undertaken.

Conclusions

The sialome of G. m. morsitans contains over 250 proteins that are possibly associated with blood feeding. This set includes alleles of previously described gene products, reveals new evidence that several salivary proteins are multigenic and identifies at least seven new polypeptide families unique to Glossina. Most of these proteins have no known function and thus, provide a discovery platform for the identification of novel pharmacologically active compounds, innovative vector-based vaccine targets, and immunological markers of vector exposure.

DNA methylation, the only known covalent modification of mammalian DNA, occurs primarily in CpG dinucleotides. 51% of CpGs in the human genome reside within repeats, and 25% within Alu elements. Despite that, no method has been reported for large-scale ascertainment of CpG methylation in repeats. Here we describe a sequencing-based strategy for parallel determination of the CpG-methylation status of thousands of Alu repeats, and a computation algorithm to design primers that enable their specific amplification from bisulfite converted genomic DNA. Using a single primer pair, we generated amplicons of high sequence complexity, and derived CpG-methylation data from 31 178 Alu elements and their 5′ flanking sequences, altogether representing over 4 Mb of a human cerebellum epigenome. The analysis of the Alu methylome revealed that the methylation level of Alu elements is high in the intronic and intergenic regions, but low in the regions close to transcription start sites. Several hypomethylated Alu elements were identified and their hypomethylated status verified by pyrosequencing. Interestingly, some Alu elements exhibited a strikingly tissue-specific pattern of methylation. We anticipate the amplicons herein described to prove invaluable as epigenome representations, to monitor epigenomic alterations during normal development, in aging and in diseases such as cancer.

Background

The medfly, Ceratitis capitata, is a highly invasive agricultural pest that has become a model insect for the development of biological control programs. Despite research into the behavior and classical and population genetics of this organism, the quantity of sequence data available is limited. We have utilized an expressed sequence tag (EST) approach to obtain detailed information on transcriptome signatures that relate to a variety of physiological systems in the medfly; this information emphasizes on reproduction, sex determination, and chemosensory perception, since the study was based on normalized cDNA libraries from embryos and adult heads.

Results

A total of 21,253 high-quality ESTs were obtained from the embryo and head libraries. Clustering analyses performed separately for each library resulted in 5201 embryo and 6684 head transcripts. Considering an estimated 19% overlap in the transcriptomes of the two libraries, they represent about 9614 unique transcripts involved in a wide range of biological processes and molecular functions. Of particular interest are the sequences that share homology with Drosophila genes involved in sex determination, olfaction, and reproductive behavior. The medfly transformer2 (tra2) homolog was identified among the embryonic sequences, and its genomic organization and expression were characterized.

Conclusion

The sequences obtained in this study represent the first major dataset of expressed genes in a tephritid species of agricultural importance. This resource provides essential information to support the investigation of numerous questions regarding the biology of the medfly and other related species and also constitutes an invaluable tool for the annotation of complete genome sequences. Our study has revealed intriguing findings regarding the transcript regulation of tra2 and other sex determination genes, as well as insights into the comparative genomics of genes implicated in chemosensory reception and reproduction.

Background

Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.

The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.

Results

We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.

Conclusion

In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Background

Dinoflagellates are important marine primary producers and grazers and cause toxic "red tides". These taxa are characterized by many unique features such as immense genomes, the absence of nucleosomes, and photosynthetic organelles (plastids) that have been gained and lost multiple times. We generated EST sequences from non-normalized and normalized cDNA libraries from a culture of the toxic species Alexandrium tamarense to elucidate dinoflagellate evolution. Previous analyses of these data have clarified plastid origin and here we study the gene content, annotate the ESTs, and analyze the genes that are putatively involved in DNA packaging.

Results

Approximately 20% of the 6,723 unique (11,171 total 3'-reads) ESTs data could be annotated using Blast searches against GenBank. Several putative dinoflagellate-specific mRNAs were identified, including one novel plastid protein. Dinoflagellate genes, similar to other eukaryotes, have a high GC-content that is reflected in the amino acid codon usage. Highly represented transcripts include histone-like (HLP) and luciferin binding proteins and several genes occur in families that encode nearly identical proteins. We also identified rare transcripts encoding a predicted protein highly similar to histone H2A.X. We speculate this histone may be retained for its role in DNA double-strand break repair.

Conclusion

This is the most extensive collection to date of ESTs from a toxic dinoflagellate. These data will be instrumental to future research to understand the unique and complex cell biology of these organisms and for potentially identifying the genes involved in toxin production.

Introduction

We have examined expression of microRNAs (miRNAs) in ependymomas to identify molecular markers of value for clinical management. miRNAs are non-coding RNAs that can block mRNA translation and affect mRNA stability. Changes in the expression of miRNAs have been correlated with many human cancers.

Materials and Methods

We have utilized TaqMan Low Density Arrays to evaluate the expression of 365 miRNAs in ependymomas and normal brain tissue. We first demonstrated the similarity of expression profiles of paired frozen tissue (FT) and paraffin-embedded specimens (FFPE). We compared the miRNA expression profiles of 34 FFPE ependymoma samples with 8 microdissected normal brain tissue specimens enriched for ependymal cells. miRNA expression profiles were then correlated with tumor location, histology and other clinicopathological features.

Results

We have identified miRNAs that are over-expressed in ependymomas, such as miR-135a and miR-17-5p, and down-regulated, such as miR-383 and miR-485-5p. We have also uncovered associations between expression of specific miRNAs which portend a worse prognosis. For example, we have identified a cluster of miRNAs on human chromosome 14q32 that is associated with time to relapse. We also found that miR-203 is an independent marker for relapse compared to the parameters that are currently used. Additionally, we have identified three miRNAs (let-7d, miR-596 and miR-367) that strongly correlate to overall survival.

Conclusion

We have identified miRNAs that are differentially expressed in ependymomas compared with normal ependymal tissue. We have also uncovered significant associations of miRNAs with clinical behavior. This is the first report of clinically relevant miRNAs in ependymomas.

Background

Chondrosarcomas are malignant cartilage tumors that do not respond to traditional chemotherapy or radiation. The 5-year survival rate of histologic grade III chondrosarcoma is less than 30%. An animal model of chondrosarcoma has been established - namely, the Swarm Rat Chondrosarcoma (SRC) - and shown to resemble the human disease. Previous studies with this model revealed that tumor microenvironment could significantly influence chondrosarcoma malignancy.

Methods

To examine the effect of the microenvironment, SRC tumors were initiated at different transplantation sites. Pyrosequencing assays were utilized to assess the DNA methylation of the tumors, and SAGE libraries were constructed and sequenced to determine the gene expression profiles of the tumors. Based on the gene expression analysis, subsequent functional assays were designed to determine the relevancy of the specific genes in the development and progression of the SRC.

Results

The site of transplantation had a significant impact on the epigenetic and gene expression profiles of SRC tumors. Our analyses revealed that SRC tumors were hypomethylated compared to control tissue, and that tumors at each transplantation site had a unique expression profile. Subsequent functional analysis of differentially expressed genes, albeit preliminary, provided some insight into the role that thymosin-β4, c-fos, and CTGF may play in chondrosarcoma development and progression.

Conclusion

This report describes the first global molecular characterization of the SRC model, and it demonstrates that the tumor microenvironment can induce epigenetic alterations and changes in gene expression in the SRC tumors. We documented changes in gene expression that accompany changes in tumor phenotype, and these gene expression changes provide insight into the pathways that may play a role in the development and progression of chondrosarcoma. Furthermore, specific functional analysis indicates that thymosin-β4 may have a role in chondrosarcoma metastasis.

Abnormal patterns of DNA methylation are observed in several types of human cancer. While localized DNA methylation of CpG islands has been associated with gene silencing, the effect that genome-wide loss of methylation has on tumorigenesis is not completely known. To examine its effect on tumorigenesis, we induced DNA demethylation in a rat model of human chondrosarcoma using 5-aza-2-deoxycytidine. Rat specific pyrosequencing assays were utilized to assess the methylation levels in both LINEs and satellite DNA sequences following 5-aza-2-deoxycytidine treatment. Loss of DNA methylation was accompanied by an increase in invasiveness of the rat chondrosarcoma cells, in vitro, as well as by an increase in tumor growth in vivo. Subsequent microarray analysis provided insight into the gene expression changes that result from 5-aza-2-deoxycytidine induced DNA demethylation. In particular, two genes that may function in tumorigenesis, sox-2 and midkine, were expressed at low levels in control cells but upon 5-aza-2-deoxycytidine treatment these genes became overexpressed. Promoter region DNA analysis revealed that these genes were methylated in control cells but became demethylated following 5-aza-2-deoxycytidine treatment. Following withdrawal of 5-aza-2-deoxycytidine, the rat chondrosarcoma cells reestablished global DNA methylation levels that were comparable to that of control cells. Concurrently, invasiveness of the rat chondrosarcoma cells, in vitro, decreased to a level indistinguishable to that of control cells. Taken together these experiments demonstrate that global DNA hypomethylation induced by 5-aza-2-deoxycytidine may promote specific aspects of tumorigenesis in rat chondrosarcoma cells.