Extensive miRNA expression analysis in craniopharyngiomas.

PURPOSE: Craniopharyngiomas are benign tumors of the sellar or parasellar regions. They arise from the remnants of Rathke's pouch and are considered a "developmental disease." microRNAs are short non-coding RNAs that play a key regulatory role in the control of expression of entire gene networks. We performed an extensive analysis of miRNAs in craniopharyngiomas aiming to identify a miRNA expression signature that might aid in the prognosis of disease progression and outcome. METHODS: Thirty-seven craniopharyngioma samples from twenty-three patients, ten age-matched controls from autopsy, and ten infant controls from the developing pituitary from autopsy were evaluated for the expression of 754 miRNAs using TaqMan® Low Density Arrays (TLDAs) v2.0 (Applied Biosystems, Foster City, CA). RESULTS: Among the most differentially expressed miRNAs, downregulation of miR-132 appears to be a marker of aggressiveness and also plays a role in epithelial-mesenchymal transition. CONCLUSIONS: This is the first time that an extensive study of miRNA expression has been performed in craniopharyngiomas. Further research needs to be performed to investigate the potential role of miR-132 in the development and progression of craniopharyngiomas, and its value as a prognostic marker of aggressiveness.

A gene signature for a long-term survivor of an atypical teratoid/rhabdoid tumor.

Atypical teratoid/rhabdoid tumors (AT/RTs) are aggressive brain tumors that are commonly associated with a dismal prognosis. However, there have been isolated reports of long-term survival that was not necessarily correlated with other prognostic factors such as age, clinical stage, or extent of surgical resection. Here, we report the case of a 6-year-old boy with AT/RT who remained disease-free for 8 years after undergoing subtotal surgical resection followed only by radiation therapy. On recurrence, the tumor rapidly progressed, leading to the patient's death a short time later. To further characterize this case and learn more about the tumor biology of long-term survivors, we compared the gene expression (GE) profiles from representative samples obtained from primary, recurrent, and progressive disease tumors of the above-mentioned patient along with a cohort of primary untreated AT/RT samples using cDNA microarrays. Global GE analysis and unsupervised hierarchical clustering showed the three events clustered together and distinctly apart from the rest of the samples. This indicates a GE background that is maintained throughout the course of the disease. This unique case suggests that there may be specific clinical characteristics associated with distinctive molecular subtypes of AT/RT. The identification and characterization of AT/RT subtypes could lead to advances in both prognosis and treatment of these tumors.

Study of the gene expression and microRNA expression profiles of malignant rhabdoid tumors originated in the brain (AT/RT) and in the kidney (RTK).

PURPOSE: Malignant rhabdoid tumors (MRT) can occur in a variety of anatomical sites. The most frequent locations are the brain, where they are named atypical teratoid/rhabdoid tumors (AT/RT), and the kidney, where they are named rhabdoid tumors of the kidney (RTK). MRTs at all sites are recognized as the same entity due to their similar morphology, aggressive behavior, and a common genetic abnormality, an inactivating mutation of the SMARCB1/INI-1/hSNF5/BAF47 gene. We aim to investigate potential molecular differences between AT/RT and RTK. METHODS: We analyzed the microRNA (miRNA) and gene expression (GE) profiles of 10 RTK, 13 AT/RT, and 2 human MRT cell lines (G401-RTK and MON-AT/RT). Illumina V2 MicroRNA Chips (Illumina, Inc., CA, USA) were used for miRNA analysis, and Illumina HT-12 whole genome expression arrays were used for GE analysis. RESULTS: The distribution of p values from GE showed a significant difference between RTK and AT/RT, with 20 % of the genes having p values ≤0.05 and the principal component analysis of the GE data showed separation between RTK and AT/RT. However, the miRNA expression failed to identify the different tumor groups. Among the 122 genes significantly differentially expressed between AT/RT and RTK, we found both genes related to brain development (i.e., FABP7, 22-fold increase in AT/RT) and genes related to kidney development (i.e., TCF21, sixfold increase in RTK). CONCLUSION: Based on our results, we hypothesized that although MRT are indeed the same tumor, independent of the site of origin, the GE differences reflect the influence of microenvironment over tumor development.

Histone deacetylases expression in atypical teratoid rhabdoid tumors.

PURPOSE: Atypical teratoid rhabdoid tumors (ATRTs) are rare, highly malignant central nervous system tumors that occur during infancy and early childhood. Their poor outcome and resistance to conventional chemotherapies and radiotherapy, urges the development of new therapies. Recent studies have evaluated the effects of histone deacetylase inhibitors (HDACi) as a new potential treatment for ATRTs. However, most HDACi act unselectively against all, or at least several, histone deacetylase (HDAC) family members. We hypothesized that specific HDAC family members are deregulated in ATRT and therefore a more selective class of HDACi would be beneficial to patients with ATRT. METHODS: To test our hypothesis, we evaluated the expression level of different HDAC family members in ATRTs. Eight ATRTs were compared to six medulloblastoma samples in regards to the level of expression of the 18 HDAC family members as determined by microarray gene expression profiling. RESULTS: HDAC1 was the only member of the HDAC family to be significantly differentially expressed in ATRTs (FC = 4.728; p value = 0.00003). CONCLUSIONS: A class of HDACi specifically targeting HDAC1 may allow for the desired therapeutic benefits with fewer side effects for children with ATRT.

Putative multifunctional signature of lung metastases in dedifferentiated chondrosarcoma.

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.

Identification of microRNAs as potential prognostic markers in ependymoma.

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.

BTECH: a platform to integrate genomic, transcriptomic and epigenomic alterations in brain tumors.

The identification of molecular signatures predictive of clinical behavior and outcome in brain tumors has been the focus of many studies in the recent years. Despite the wealth of data that are available in the public domain on alterations in the genome, epigenome and transcriptome of brain tumors, the underlying molecular mechanisms leading to tumor initiation and progression remain largely unknown. Unfortunately, most of these data are scattered in multiple databases and supplementary materials of publications, thus making their retrieval, evaluation, comparison and visualization a rather arduous task. Here we report the development and implementation of an open access database (BTECH), a community resource for the deposition of a wide range of molecular data derived from brain tumor studies. This comprehensive database integrates multiple datasets, including transcript profiles, epigenomic CpG methylation data, DNA copy number alterations and structural chromosomal rearrangements, tumor-associated gene lists, SNPs, genomic features concerning Alu repeats and general genomic annotations. A genome browser has also been developed that allows for the simultaneous visualization of the different datasets and the various annotated features. Besides enabling an integrative view of diverse datasets through the genome browser, we also provide links to the original references for users to have a more accurate understanding of each specific dataset. This integrated platform will facilitate uncovering interactions among genetic and epigenetic factors associated with brain tumor development. BTECH is freely available at http://cmbteg.childrensmemorial.org/.

Microenvironment alters epigenetic and gene expression profiles in Swarm rat chondrosarcoma tumors.

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.

Transcriptional patterns in both host and bacterium underlie a daily rhythm of anatomical and metabolic change in a beneficial symbiosis.

Mechanisms for controlling symbiont populations are critical for maintaining the associations that exist between a host and its microbial partners. We describe here the transcriptional, metabolic, and ultrastructural characteristics of a diel rhythm that occurs in the symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri. The rhythm is driven by the host's expulsion from its light-emitting organ of most of the symbiont population each day at dawn. The transcriptomes of both the host epithelium that supports the symbionts and the symbiont population itself were characterized and compared at four times over this daily cycle. The greatest fluctuation in gene expression of both partners occurred as the day began. Most notable was an up-regulation in the host of >50 cytoskeleton-related genes just before dawn and their subsequent down-regulation within 6 h. Examination of the epithelium by TEM revealed a corresponding restructuring, characterized by effacement and blebbing of its apical surface. After the dawn expulsion, the epithelium reestablished its polarity, and the residual symbionts began growing, repopulating the light organ. Analysis of the symbiont transcriptome suggested that the bacteria respond to the effacement by up-regulating genes associated with anaerobic respiration of glycerol; supporting this finding, lipid analysis of the symbionts' membranes indicated a direct incorporation of host-derived fatty acids. After 12 h, the metabolic signature of the symbiont population shifted to one characteristic of chitin fermentation, which continued until the following dawn. Thus, the persistent maintenance of the squid-vibrio symbiosis is tied to a dynamic diel rhythm that involves both partners.

Upregulation of mir-221 and mir-222 in atypical teratoid/rhabdoid tumors: potential therapeutic targets.

PURPOSE: The aim of this study is to search for new therapeutic targets for atypical teratoid-rhabdoid tumors (ATRT). METHODS: To achieve this, we compared the expression of 365 microRNAs among ATRT, medulloblastomas, and normal brain. RESULTS: MiR-221 and miR-222 were within the top differentially expressed microRNAs. The deregulated expression of miR221/222 was demonstrated to inhibit the expression of the tumor suppressor and inhibitor of cell cycle p27(Kip1). Here, we demonstrated the negative regulation of p27(Kip1) by miR-221/222 in ATRT using microarray, real-time reverse transcriptase polymerase chain reaction, and immunohistochemistry. CONCLUSION: As anti-miR therapy was recently proposed as an alternative treatment for cancer, these findings suggest that anti-miR-221/222 therapy might have therapeutic potential in ATRT.

SNP-based prediction of the human germ cell methylation landscape.

Base substitution occurs at a high rate at CpG dinucleotides due to the frequent methylation of CpG and the deamination of methylated cytosine to thymine. If these substitutions occur in germ cells, they constitute a heritable mutation that may eventually rise to polymorphic frequencies, hence resulting in a SNP that is methylation associated. In this study, we sought to identify clusters of methylation associated SNPs as a basis for prediction of methylation landscapes of germ cell genomes. Genomic regions enriched with methylation associated SNPs, namely "methylation associated SNP clusters", were identified with an agglomerative hierarchical clustering algorithm. Repetitive elements, segmental duplications, and syntenic tandem DNA repeats were enriched in methylation associated SNP clusters. The frequency of methylation associated SNPs in Alu Y/S elements exhibited a gradient pattern suggestive of linear spreading, being higher in proximity to methylation associated SNP clusters and lower closer to CpG islands. Interestingly, methylation associated SNP clusters were over-represented near the transcriptional initiation sites of immune response genes. We propose a de novo DNA methylation model during germ cell development whereby a pattern is established by long-range chromatic interactions through syntenic repeats combined with regional methylation spreading from methylation associated SNP clusters.

Epigenetically reprogramming metastatic tumor cells with an embryonic microenvironment.

UNLABELLED: 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.

Effects of colonization, luminescence, and autoinducer on host transcription during development of the squid-vibrio association.

The light-organ symbiosis between the squid Euprymna scolopes and the luminous bacterium Vibrio fischeri offers the opportunity to decipher the hour-by-hour events that occur during the natural colonization of an animal's epithelial surface by its microbial partners. To determine the genetic basis of these events, a glass-slide microarray was used to characterize the light-organ transcriptome of juvenile squid in response to the initiation of symbiosis. Patterns of gene expression were compared between animals not exposed to the symbiont, exposed to the wild-type symbiont, or exposed to a mutant symbiont defective in either of two key characters of this association: bacterial luminescence or autoinducer (AI) production. Hundreds of genes were differentially regulated as a result of symbiosis initiation, and a hierarchy existed in the magnitude of the host's response to three symbiont features: bacterial presence > luminescence > AI production. Putative host receptors for bacterial surface molecules known to induce squid development are up-regulated by symbiont light production, suggesting that bioluminescence plays a key role in preparing the host for bacteria-induced development. Further, because the transcriptional response of tissues exposed to AI in the natural context (i.e., with the symbionts) differed from that to AI alone, the presence of the bacteria potentiates the role of quorum signals in symbiosis. Comparison of these microarray data with those from other symbioses, such as germ-free/conventionalized mice and zebrafish, revealed a set of shared genes that may represent a core set of ancient host responses conserved throughout animal evolution.

An annotated cDNA library of juvenile Euprymna scolopes with and without colonization by the symbiont Vibrio fischeri.

BACKGROUND: Biologists are becoming increasingly aware that the interaction of animals, including humans, with their coevolved bacterial partners is essential for health. This growing awareness has been a driving force for the development of models for the study of beneficial animal-bacterial interactions. In the squid-vibrio model, symbiotic Vibrio fischeri induce dramatic developmental changes in the light organ of host Euprymna scolopes over the first hours to days of their partnership. We report here the creation of a juvenile light-organ specific EST database. RESULTS: We generated eleven cDNA libraries from the light organ of E. scolopes at developmentally significant time points with and without colonization by V. fischeri. Single pass 3' sequencing efforts generated 42,564 expressed sequence tags (ESTs) of which 35,421 passed our quality criteria and were then clustered via the UIcluster program into 13,962 nonredundant sequences. The cDNA clones representing these nonredundant sequences were sequenced from the 5' end of the vector and 58% of these resulting sequences overlapped significantly with the associated 3' sequence to generate 8,067 contigs with an average sequence length of 1,065 bp. All sequences were annotated with BLASTX (E-value < -03) and Gene Ontology (GO). CONCLUSION: Both the number of ESTs generated from each library and GO categorizations are reflective of the activity state of the light organ during these early stages of symbiosis. Future analyses of the sequences identified in these libraries promise to provide valuable information not only about pathways involved in colonization and early development of the squid light organ, but also about pathways conserved in response to bacterial colonization across the animal kingdom.

Pax9 and Jagged1 act downstream of Gli3 in vertebrate limb development.

From early in limb development the transcription factor Gli3 acts to define boundaries of gene expression along the anterior-posterior (AP) axis, establishing asymmetric patterns required to provide positional information. As limb development proceeds, posterior mesenchyme expression of Sonic hedgehog (Shh) regulates Gli3 transcription and post-translational processing to specify digit number and identity. The molecular cascades dependent on Gli3 at later stages of limb development, which link early patterning events with final digit morphogenesis, remain poorly characterised. By analysing the transcriptional consequences of loss of Gli3 in the anterior margin of the E11.5 and E12.5 limb bud in the polydactylous mouse mutant extra-toes (Gli3(Xt/Xt)), we have identified a number of known and novel transcripts dependent on Gli3 in the limb. In particular, we demonstrated that the genes encoding the paired box transcription factor Pax9, the Notch ligand Jagged1 and the cell surface receptor Cdo are dependent on Gli3 for correct expression in the anterior limb mesenchyme. Analysis of expression in compound Shh;Gli3 mutant mouse embryos and in both in vitro and in vivo Shh signaling assays, further defined the importance of Shh regulated processing of Gli3 in controlling gene expression. In particular Pax9 regulation by Shh and Gli3 was shown to be context dependent, with major differences between the limb and somite revealed by Shh bead implantation experiments in the chick. Jagged1 was shown to be induced by Shh in the chick limb and in a C3H10T1/2 cell based signaling assay, with Shh;Gli3 mutant analysis indicating that expression is dependent on Gli3 derepression. Our data have also revealed that perturbation of early patterning events within the Gli3(Xt/Xt) limb culminates in a specific delay of anterior chondrogenesis which is subsequently realised as extra digits.

Integrative annotation of 21,037 human genes validated by full-length cDNA clones.

The human genome sequence defines our inherent biological potential; the realization of the biology encoded therein requires knowledge of the function of each gene. Currently, our knowledge in this area is still limited. Several lines of investigation have been used to elucidate the structure and function of the genes in the human genome. Even so, gene prediction remains a difficult task, as the varieties of transcripts of a gene may vary to a great extent. We thus performed an exhaustive integrative characterization of 41,118 full-length cDNAs that capture the gene transcripts as complete functional cassettes, providing an unequivocal report of structural and functional diversity at the gene level. Our international collaboration has validated 21,037 human gene candidates by analysis of high-quality full-length cDNA clones through curation using unified criteria. This led to the identification of 5,155 new gene candidates. It also manifested the most reliable way to control the quality of the cDNA clones. We have developed a human gene database, called the H-Invitational Database (H-InvDB; http://www.h-invitational.jp/). It provides the following: integrative annotation of human genes, description of gene structures, details of novel alternative splicing isoforms, non-protein-coding RNAs, functional domains, subcellular localizations, metabolic pathways, predictions of protein three-dimensional structure, mapping of known single nucleotide polymorphisms (SNPs), identification of polymorphic microsatellite repeats within human genes, and comparative results with mouse full-length cDNAs. The H-InvDB analysis has shown that up to 4% of the human genome sequence (National Center for Biotechnology Information build 34 assembly) may contain misassembled or missing regions. We found that 6.5% of the human gene candidates (1,377 loci) did not have a good protein-coding open reading frame, of which 296 loci are strong candidates for non-protein-coding RNA genes. In addition, among 72,027 uniquely mapped SNPs and insertions/deletions localized within human genes, 13,215 nonsynonymous SNPs, 315 nonsense SNPs, and 452 indels occurred in coding regions. Together with 25 polymorphic microsatellite repeats present in coding regions, they may alter protein structure, causing phenotypic effects or resulting in disease. The H-InvDB platform represents a substantial contribution to resources needed for the exploration of human biology and pathology.

A comprehensive nonredundant expressed sequence tag collection for the developing Rattus norvegicus heart.

Congenital heart defects affect approximately 1,000,000 people in the United States, with 40,000 new births contributing to that number every year. A large percentage of these defects can be attributed to septal defects. We assembled a nonredundant collection of over 12,000 expressed sequence tags (ESTs) from a total of 30,000 ESTs, with the ultimate goal of identifying spatially and/or temporally regulated genes during heart septation. These ESTs were compiled from nonnormalized, normalized, and serially subtracted cDNA libraries derived from two sets of tissue samples. The first includes microdissected rat hearts from embryonic (E) days E13, E15, and E16.5-E18.5 and adult heart. The second includes hearts from embryonic days E17, E19, and E21 and postnatal (P) days P1, P12, P74, and P200. Over 6,000 novel ESTs were identified in the libraries derived from these two sets of tissues, all of which have been contributed to the NCBI rat UniGene collection. It is anticipated that such EST and cDNA clone resources will prove invaluable to gene expression studies aimed at the understanding of the molecular mechanisms underlying heart septation defects.