Elevated expression of the centromere protein-A(CENP-A)-encoding gene as a prognostic and predictive biomarker in human cancers.

Centromere protein-A (CENP-A), a histone-H3 variant, plays an essential role in cell division by ensuring proper formation and function of centromeres and kinetochores. Elevated CENP-A expression has been associated with cancer development. This study aimed to establish whether elevated CENP-A expression can be used as a prognostic and predictive cancer biomarker. Molecular profiling of CENP-A in human cancers was investigated using genomic, transcriptomic and patient information from databases, including COSMIC, Oncomine, Kaplan-Meier plotter and cBioPortal. A network of CENP-A co-expressed genes was derived from cBioPortal and analyzed using Ingenuity Pathway Analysis (IPA) and Oncomine protocols to explore the function of CENP-A and its predictive potential. Transcriptional and post-transcriptional regulation of CENP-A expression was analyzed in silico. It was found that CENP-A expression was elevated in 20 types of solid cancer compared with normal counterparts. Elevated CENP-A expression highly correlated with cancer progression and poor patient outcome. Genomic analysis indicated that the elevated CENP-A expression was not due to alterations in the sequence or copy number of the CENP-A gene. Furthermore, CENP-A can be regulated by key oncogenic proteins and tumor-suppressive microRNAs. CENP-A co-expression network analysis indicated that CENP-A function is associated with cell cycle progression. Oncomine analysis showed a strong correlation between elevated CENP-A expression and oncolytic response of breast cancer patients to taxane-based chemotherapy. In conclusion, elevated CENP-A expression is coupled to malignant progression of numerous types of cancer. It may be useful as a biomarker of poor patient prognosis and as a predictive biomarker for taxane-based chemotherapy.

Identification of the epigenetic reader CBX2 as a potential drug target in advanced prostate cancer.

BACKGROUND: While localized prostate cancer (PCa) can be effectively cured, metastatic disease inevitably progresses to a lethal state called castration-resistant prostate cancer (CRPC). Emerging evidence suggests that aberrant epigenetic repression by the polycomb group (PcG) complexes fuels PCa progression, providing novel therapeutic opportunities. RESULTS: In the search for potential epigenetic drivers of CRPC, we analyzed the molecular profile of PcG members in patient-derived xenografts and clinical samples. Overall, our results identify the PcG protein and methyl-lysine reader CBX2 as a potential therapeutic target in advanced PCa. We report that CBX2 was recurrently up-regulated in metastatic CRPC and that elevated CBX2 expression was correlated with poor clinical outcome in PCa cohorts. Furthermore, CBX2 depletion abrogated cell viability and induced caspase 3-mediated apoptosis in metastatic PCa cell lines. Mechanistically explaining this phenotype, microarray analysis in CBX2-depleted cells revealed that CBX2 controls the expression of many key regulators of cell proliferation and metastasis. CONCLUSIONS: Taken together, this study provides the first evidence that CBX2 inhibition induces cancer cell death, positioning CBX2 as an attractive drug target in lethal CRPC.

Integrated analysis of the prostate cancer small-nucleolar transcriptome reveals SNORA55 as a driver of prostate cancer progression.

Metastasis is the primary cause of death in prostate cancer (PCa) patients. Small nucleolar RNAs (snoRNAs) have long been considered "housekeeping" genes with no relevance for cancer biology. Emerging evidence has challenged this assumption, suggesting that snoRNA expression is frequently modulated during cancer progression. Despite this, no study has systematically addressed the prognostic and functional significance of snoRNAs in PCa. We performed RNA Sequencing on paired metastatic/non-metastatic PCa xenografts derived from clinical specimens. The clinical significance of differentially expressed snoRNAs was further investigated in two independent primary PCa cohorts (131 and 43 patients, respectively). The snoRNA demonstrating the strongest association with clinical outcome was quantified in PCa patient-derived serum samples and its functional relevance was investigated in PCa cells via gene expression profiling, pathway analysis and gene silencing. Our comparison revealed 21 differentially expressed snoRNAs in the metastatic vs. non-metastatic xenografts. Of those, 12 were represented in clinical databases and were further analyzed. SNORA55 emerged as a predictor of shorter relapse-free survival (results confirmed in two independent databases). SNORA55 was reproducibly detectable in serum samples from PCa patients. SNORA55 silencing in PCa cell lines significantly inhibited cell proliferation and migration. Pathway analysis revealed that SNORA55 expression is significantly associated with growth factor signaling and pro-inflammatory cytokine expression in PCa. Our results demonstrate that SNORA55 up-regulation predicts PCa progression and that silencing this non-coding gene affects PCa cell proliferation and metastatic potential, thus positioning it as both a novel biomarker and therapeutic target.

DNA methylation at enhancer regions: Novel avenues for epigenetic biomarker development.

Biomarkers are molecules or features which can provide clinically-relevant information about a particular disease state, thus providing useful tools for oncologists. Recently, a number of studies have demonstrated that DNA methylation holds great promise as a novel source of cancer biomarkers. Although promoter regions have been the focus of most investigations thus far, mounting evidence demonstrates that enhancer sequences also undergo extensive differential methylation in cancer cells. Moreover, enhancer methylation correlates with target gene expression better than promoter methylation, providing unexplored strategies for biomarker development. Here, we review important considerations associated with the clinical analysis of DNA methylation at distal regulatory regions. Notably, we highlight emerging literature addressing the methylation status of enhancers in development and cancer, and subsequently discuss how enhancer methylation can be exploited to guide disease management. While acknowledging current limitations, we propose that the methylation state of enhancer regions has the potential to headline the next generation of epigenetic biomarkers.

Polycomb genes are associated with response to imatinib in chronic myeloid leukemia.

AIM: Imatinib is a tyrosine kinase inhibitor that has revolutionized the treatment of chronic myeloid leukemia (CML). Despite its efficacy, about a third of patients discontinue the treatment due to therapy failure or intolerance. The rational identification of patients less likely to respond to imatinib would be of paramount clinical relevance. We have shown that transmembrane transporter hOCT1 genotyping predicts imatinib activity. In parallel, Polycomb group genes (PcGs) are epigenetic repressors implicated in CML progression and in therapy resistance. PATIENTS & METHODS: We measured the expression of eight PcGs in paired pre- and post-imatinib bone marrow samples from 30 CML patients. RESULTS: BMI1, PHC3, CBX6 and CBX7 expression was significantly increased during imatinib treatment. Post-treatment levels of CBX6 and CBX7 predicted 3-month response rate. Measurement of post-treatment BMI1 levels improved the predictive power of hOCT1 genotyping. CONCLUSION: These results suggest that the expression levels of PcGs might be useful for a more accurate risk stratification of CML patients.

Polycomb-mediated silencing in neuroendocrine prostate cancer.

BACKGROUND: Neuroendocrine prostate cancer (NEPC) is a highly aggressive subtype of prostate cancer (PCa) for which the median survival remains less than a year. Current treatments are only palliative in nature, and the lack of suitable pre-clinical models has hampered previous efforts to develop novel therapeutic strategies. Addressing this need, we have recently established the first in vivo model of complete neuroendocrine transdifferentiation using patient-derived xenografts. Few genetic differences were observed between parental PCa and relapsed NEPC, suggesting that NEPC likely results from alterations that are epigenetic in nature. Thus, we sought to identify targetable epigenetic regulators whose expression was elevated in NEPC using genome-wide profiling of patient-derived xenografts and clinical samples. RESULTS: Our data indicate that multiple members of the polycomb group (PcG) family of transcriptional repressors were selectively upregulated in NEPC. Notably, CBX2 and EZH2 were consistently the most highly overexpressed epigenetic regulators across multiple datasets from clinical and xenograft tumor tissues. Given the striking upregulation of PcG genes and other transcriptional repressors, we derived a 185-gene list termed 'neuroendocrine-associated repression signature' (NEARS) by overlapping transcripts downregulated across multiple in vivo NEPC models. In line with the striking upregulation of PcG family members, NEARS was preferentially enriched with PcG target genes, suggesting a driving role for PcG silencing in NEPC. Importantly, NEARS was significantly associated with high-grade tumors, metastatic progression, and poor outcome in multiple clinical datasets, consistent with extensive literature linking PcG genes and aggressive disease progression. CONCLUSIONS: We have explored the epigenetic landscape of NEPC and provided evidence of increased PcG-mediated silencing associated with aberrant transcriptional regulation of key differentiation genes. Our results position CBX2 and EZH2 as potential therapeutic targets in NEPC, providing opportunities to explore novel strategies aimed at reversing epigenetic alterations driving this lethal disease.

The long non-coding RNA PCGEM1 is regulated by androgen receptor activity in vivo.

BACKGROUND: Long non-coding RNAs (lncRNAs) can orchestrate oncogenic or tumor-suppressive functions in cancer biology. Accordingly, PCGEM1 and PRNCR1 were implicated in progression of prostate cancer (PCa) as transcriptional co-regulators of the androgen receptor (AR). However, these findings were recently refuted asserting that neither gene physically binds to the AR. Despite evidence for differing AR transcriptional programs in vivo and in vitro, studies investigating AR-regulation of these genes hitherto have only been conducted in vitro. Here, we further examine the relevance of PCGEM1 and PRNCR1 in PCa, and their relationship with AR signaling, using patient-derived xenograft models. FINDINGS: RNA sequencing of two distinct androgen-dependent models shows PCGEM1 to be considerably expressed, while PRNCR1 showed scant basal expression. PCGEM1 was sharply down-regulated following castration and up-regulated upon AR activation in vivo. However, we found no parallel evidence following AR stimulation in vitro. A PCGEM1-associated gene expression signature (PES) was significantly repressed in response to androgen ablation therapy and in hormone-refractory versus hormone-naïve PCa patients. Furthermore, we found PCGEM1 was uniformly distributed in PCa cell nucleus and cytoplasm which remained unaltered upon AR transcriptional activation. PCGEM1 was up-regulated in primary PCa but not in metastasized PCa. Accordingly, the PES was significantly down-regulated in advanced and higher grade PCa patients from multiple independent studies. CONCLUSION: Our results demonstrate PCGEM1 as an in vivo androgen-regulated transcript with potential nuclear and/or cytoplasmic function(s). Importantly, the clinical expression profile of PCGEM1 implicates it in the early stages of PCa warranting further research in this direction.

Identification of a long non-coding RNA as a novel biomarker and potential therapeutic target for metastatic prostate cancer.

Metastatic prostate cancer (PCa) is still an incurable disease. Long non-coding RNAs (lncRNAs) may be an overlooked source of cancer biomarkers and therapeutic targets. We therefore performed RNA sequencing on paired metastatic/non-metastatic PCa xenografts derived from clinical specimens. The most highly up-regulated transcript was LOC728606, a lncRNA now designated PCAT18. PCAT18 is specifically expressed in the prostate compared to 11 other normal tissues (p<0.05) and up-regulated in PCa compared to 15 other neoplasms (p<0.001). Cancer-specific up-regulation of PCAT18 was confirmed on an independent dataset of PCa and benign prostatic hyperplasia samples (p<0.001). PCAT18 was detectable in plasma samples and increased incrementally from healthy individuals to those with localized and metastatic PCa (p<0.01). We identified a PCAT18-associated expression signature (PES), which is highly PCa-specific and activated in metastatic vs. primary PCa samples (p<1E-4, odds ratio>2). The PES was significantly associated with androgen receptor (AR) signalling. Accordingly, AR activation dramatically up-regulated PCAT18 expression in vitro and in vivo. PCAT18 silencing significantly (p<0.001) inhibited PCa cell proliferation and triggered caspase 3/7 activation, with no effect on non-neoplastic cells. PCAT18 silencing also inhibited PCa cell migration (p<0.01) and invasion (p<0.01). These results position PCAT18 as a potential therapeutic target and biomarker for metastatic PCa.

The non-coding transcriptome as a dynamic regulator of cancer metastasis.

Since the discovery of microRNAs, non-coding RNAs (NC-RNAs) have increasingly attracted the attention of cancer investigators. Two classes of NC-RNAs are emerging as putative metastasis-related genes: long non-coding RNAs (lncRNAs) and small nucleolar RNAs (snoRNAs). LncRNAs orchestrate metastatic progression through several mechanisms, including the interaction with epigenetic effectors, splicing control and generation of microRNA-like molecules. In contrast, snoRNAs have been long considered "housekeeping" genes with no relevant function in cancer. However, recent evidence challenges this assumption, indicating that some snoRNAs are deregulated in cancer cells and may play a specific role in metastasis. Interestingly, snoRNAs and lncRNAs share several mechanisms of action, and might synergize with protein-coding genes to generate a specific cellular phenotype. This evidence suggests that the current paradigm of metastatic progression is incomplete. We propose that NC-RNAs are organized in complex interactive networks which orchestrate cellular phenotypic plasticity. Since plasticity is critical for cancer cell metastasis, we suggest that a molecular interactome composed by both NC-RNAs and proteins orchestrates cancer metastasis. Interestingly, expression of lncRNAs and snoRNAs can be detected in biological fluids, making them potentially useful biomarkers. NC-RNA expression profiles in human neoplasms have been associated with patients' prognosis. SnoRNA and lncRNA silencing in pre-clinical models leads to cancer cell death and/or metastasis prevention, suggesting they can be investigated as novel therapeutic targets. Based on the literature to date, we critically discuss how the NC-RNA interactome can be explored and manipulated to generate more effective diagnostic, prognostic, and therapeutic strategies for metastatic neoplasms.

Histone modifications, stem cells and prostate cancer.

Prostate cancer (PCa) is a very common neoplasm, which is generally treated by chemo-, radio-, and/or hormonal-therapy. After a variable time, PCa becomes resistant to conventional treatment, leading to patient death. Prostate tumor-initiating cells (TICs) and cancer repopulating cells (CRCs) are stem-like populations, driving respectively cancer initiation and progression. Histone modifiers (HMs) control gene expression in normal and cancer cells, thereby orchestrating key physiological and pathological processes. In particular, Polycomb group genes (PcGs) are a set of HMs crucial for lineage-specific gene silencing and stem cell self renewal. PcG products are organized into two main Polycomb Repressive Complexes (PRCs). At specific loci, PRC2 catalyzes histone H3 Lys27 trimethylation, which triggers gene silencing by recruiting PRC1, histone deacetylases and DNA methyl transferases. PRC1 catalyzes addition of the repressive mark histone H2A ubiquitination. Recently, the catalytic component of PRC1 (BMI1) was shown to play critical roles in prostate CRC self-renewal and resistance to chemotherapy, resulting in poorer prognosis. Similarly, pharmacological disruption of PRC2 by a small molecule inhibitor reduced the tumorigenicity and metastatic potential of prostate CRCs. Along with PcGs, some histone lysine demethylases (KDMs) are emerging as critical regulators of TIC/CRC biology. KDMs may be inhibited by specific small molecules, some of which display antitumor activity in PCa cells at micromolar concentrations. Since epigenetic gene regulation is crucial for stem cell biology, exploring the role of HMs in prostate cancer is a promising path that may lead to novel treatments.

Genistein versus ICI 182, 780: an ally or enemy in metastatic progression of prostate cancer.

BACKGROUND: Androgen signalling through the androgen receptor (AR) plays a critical role in prostate cancer (PCa) initiation and progression. Estrogen in synergy with androgen is essential for cell growth of the normal and malignant prostate. However, the exact role that estrogen and the estrogen receptor play in prostate carcinogenesis remains unclear. We have previously demonstrated the metastasis-promoting effect of an estrogen receptor beta (ERß) agonist (genistein) in a patient-derived PCa xenograft model mimicking localized and metastatic disease. METHODS: To test the hypothesis that the tumor-promoting activity of genistein was due to its estrogenic properties, we treated the xenograft-bearing mice with genistein and an anti-estrogen compound (ICI 182, 780) and compared the differential gene expression using microarrays. RESULTS: Using a second xenograft model which was derived from another patient, we showed that genistein promoted disease progression in vivo and ICI 182, 780 inhibited metastatic spread. The microarray analysis revealed that the metallothionein (MT) gene family was differentially expressed in tumors treated by these compounds. Using qRT-PCR, the differences in expression levels were validated in the metastatic and non-metastatic LTL313 PCa xenograft tumor lines, both of which were originally derived from the same PCa patient. CONCLUSIONS: Together our data provide evidence that genistein stimulates and ICI 182, 780 inhibits metastatic progression, suggesting that these effects may be mediated by ERß signalling.

Plasma miRNAs as biomarkers to identify patients with castration-resistant metastatic prostate cancer.

MicroRNAs (miRNAs) have emerged as key regulators of numerous biological processes, and increasing evidence suggests that circulating miRNAs may be useful biomarkers of clinical disease. In this study, we sought to identify plasma miRNAs that differentiate patients with metastatic castration resistant prostate cancer (mCRPC) from those with localized prostate cancer (PCa). Pooled plasma samples from patients with localized PCa or mCRPC (25 per group) were assayed using the Exiqon miRNA qPCR panel, and the differential expression of selected candidates was validated using qRT-PCR. We identified 63 miRNAs upregulated in mCRPC versus localized PCa, while only four were downregulated. Pearson's correlation analysis revealed two highly correlated groups: one consisting of miR-141, miR375 and miR-200c and the other including miR151-3p, miR423-3p, miR-126, miR152 and miR-21. A third group, containing miR-16 and miR-205, showed less correlation. One miRNA from each group (miR-141, miR151-3p and miR-16) was used for logistic regression analysis and proved to increase the sensitivity of the prostate-specific antigen (PSA) test alone. While no miRNA alone differentiated localized PCa and mCRPC, combinations had greater sensitivity and specificity. The expression of these 10 candidates was assayed for association with clinical parameters of disease progression through the cBio portal. Our results demonstrate that plasma levels of selected miRNAs are potential biomarkers to differentiate localized PCa and mCRPC.

The emerging role of histone lysine demethylases in prostate cancer.

Early prostate cancer (PCa) is generally treatable and associated with good prognosis. After a variable time, PCa evolves into a highly metastatic and treatment-refractory disease: castration-resistant PCa (CRPC). Currently, few prognostic factors are available to predict the emergence of CRPC, and no curative option is available. Epigenetic gene regulation has been shown to trigger PCa metastasis and androgen-independence. Most epigenetic studies have focused on DNA and histone methyltransferases. While DNA methylation leads to gene silencing, histone methylation can trigger gene activation or inactivation, depending on the target amino acid residues and the extent of methylation (me1, me2, or me3). Interestingly, some histone modifiers are essential for PCa tumor-initiating cell (TIC) self-renewal. TICs are considered the seeds responsible for metastatic spreading and androgen-independence. Histone Lysine Demethylases (KDMs) are a novel class of epigenetic enzymes which can remove both repressive and activating histone marks. KDMs are currently grouped into 7 major classes, each one targeting a specific methylation site. Since their discovery, KDM expression has been found to be deregulated in several neoplasms. In PCa, KDMs may act as either tumor suppressors or oncogenes, depending on their gene regulatory function. For example, KDM1A and KDM4C are essential for PCa androgen-dependent proliferation, while PHF8 is involved in PCa migration and invasion. Interestingly, the possibility of pharmacologically targeting KDMs has been demonstrated. In the present paper, we summarize the emerging role of KDMs in regulating the metastatic potential and androgen-dependence of PCa. In addition, we speculate on the possible interaction between KDMs and other epigenetic effectors relevant for PCa TICs. Finally, we explore the role of KDMs as novel prognostic factors and therapeutic targets. We believe that studies on histone demethylation may add a novel perspective in our efforts to prevent and cure advanced PCa.

Differentiation of mouse embryonic stem cells into endoderm without embryoid body formation.

Pluripotent embryonic stem cells hold a great promise as an unlimited source of tissue for treatment of chronic diseases such as Type 1 diabetes. Herein, we describe a protocol using all-trans-retinoic acid, basic fibroblast growth factor and dibutyryl cAMP (DBcAMP) in the absence of embryoid body formation, for differentiation of murine embryonic stem cells into definitive endoderm that may serve as pancreatic precursors. The produced cells were analyzed by quantitative PCR, immunohistochemistry and static insulin release assay for markers of trilaminar embryo, and pancreas. Differentiated cells displayed increased Sox17 and Foxa2 expression consistent with definitive endoderm production. There was minimal production of Sox7, an extraembryonic endoderm marker, and Oct4, a marker of pluripotency. There was minimal mesoderm or neuroectoderm formation based on expression levels of the markers brachyury and Sox1, respectively. Various assays revealed that the cell clusters generated by this protocol express markers of the pancreatic lineage including insulin I, insulin II, C-peptide, PDX-1, carboxypeptidase E, pan-cytokeratin, amylase, glucagon, PAX6, Ngn3 and Nkx6.1. This protocol using all-trans-retinoic acid, DBcAMP, in the absence of embryoid bodies, generated cells that have features of definitive endoderm that may serve as pancreatic endocrine precursors.

Locus co-occupancy, nucleosome positioning, and H3K4me1 regulate the functionality of FOXA2-, HNF4A-, and PDX1-bound loci in islets and liver.

The liver and pancreas share a common origin and coexpress several transcription factors. To gain insight into the transcriptional networks regulating the function of these tissues, we globally identify binding sites for FOXA2 in adult mouse islets and liver, PDX1 in islets, and HNF4A in liver. Because most eukaryotic transcription factors bind thousands of loci, many of which are thought to be inactive, methods that can discriminate functionally active binding events are essential for the interpretation of genome-wide transcription factor binding data. To develop such a method, we also generated genome-wide H3K4me1 and H3K4me3 localization data in these tissues. By analyzing our binding and histone methylation data in combination with comprehensive gene expression data, we show that H3K4me1 enrichment profiles discriminate transcription factor occupied loci into three classes: those that are functionally active, those that are poised for activation, and those that reflect pioneer-like transcription factor activity. Furthermore, we demonstrate that the regulated presence of H3K4me1-marked nucleosomes at transcription factor occupied promoters and enhancers controls their activity, implicating both tissue-specific transcription factor binding and nucleosome remodeling complex recruitment in determining tissue-specific gene expression. Finally, we apply these approaches to generate novel insights into how FOXA2, PDX1, and HNF4A cooperate to drive islet- and liver-specific gene expression.

SHIP is required for dendritic cell maturation.

Although several groups have investigated the role of SHIP in macrophage (M) development and function, SHIP's contribution to the generation, maturation, and innate immune activation of dendritic cells (DCs) is poorly understood. We show herein that SHIP negatively regulates the generation of DCs from bone marrow precursors in vitro and in vivo, as illustrated by the enhanced expansion of DCs from SHIP(-/-) GM-CSF cultures, as well as increased numbers of DCs in the spleens of SHIP-deficient mice. Interestingly, however, these SHIP(-/-) DCs display a relatively immature phenotype and secrete substantially lower levels of IL-12 after TLR ligand stimulation than wild type DCs. This, in turn, leads to a dramatically reduced stimulation of Ag-specific T cell proliferation and Th1 cell responses in vitro and in vivo. This immature phenotype of SHIP(-/-) DCs could be reversed with the PI3K inhibitors LY294002 and wortmannin, suggesting that SHIP promotes DC maturation by reducing the levels of the PI3K second messenger phosphatidylinositol-3,4,5-trisphosphate. These results are consistent with SHIP being a negative regulator of GM-CSF-derived DC generation but a positive regulator of GM-CSF-derived DC maturation and function.

Loss-of-function Additional sex combs like 1 mutations disrupt hematopoiesis but do not cause severe myelodysplasia or leukemia.

The Additional sex combs like 1 (Asxl1) gene is 1 of 3 mammalian homologs of the Additional sex combs (Asx) gene of Drosophila. Asx is unusual because it is required to maintain both activation and silencing of Hox genes in flies and mice. Asxl proteins are characterized by an amino terminal homology domain, by interaction domains for nuclear receptors, and by a C-terminal plant homeodomain protein-protein interaction domain. A recent study of patients with myelodysplastic syndrome (MDS) and chronic myelomonocytic leukemia (CMML) revealed a high incidence of truncation mutations that would delete the PHD domain of ASXL1. Here, we show that Asxl1 is expressed in all hematopoietic cell fractions analyzed. Asxl1 knockout mice exhibit defects in frequency of differentiation of lymphoid and myeloid progenitors, but not in multipotent progenitors. We do not detect effects on hematopoietic stem cells, or in peripheral blood. Notably, we do not detect severe myelodysplastic phenotypes or leukemia in this loss-of-function model. We conclude that Asxl1 is needed for normal hematopoiesis. The mild phenotypes observed may be because other Asxl genes have redundant function with Asxl1, or alternatively, MDS or oncogenic phenotypes may result from gain-of-function Asxl mutations caused by genomic amplification, gene fusion, or truncation of Asxl1.

A seriation approach for visualization-driven discovery of co-expression patterns in Serial Analysis of Gene Expression (SAGE) data.

BACKGROUND: Serial Analysis of Gene Expression (SAGE) is a DNA sequencing-based method for large-scale gene expression profiling that provides an alternative to microarray analysis. Most analyses of SAGE data aimed at identifying co-expressed genes have been accomplished using various versions of clustering approaches that often result in a number of false positives. PRINCIPAL FINDINGS: Here we explore the use of seriation, a statistical approach for ordering sets of objects based on their similarity, for large-scale expression pattern discovery in SAGE data. For this specific task we implement a seriation heuristic we term 'progressive construction of contigs' that constructs local chains of related elements by sequentially rearranging margins of the correlation matrix. We apply the heuristic to the analysis of simulated and experimental SAGE data and compare our results to those obtained with a clustering algorithm developed specifically for SAGE data. We show using simulations that the performance of seriation compares favorably to that of the clustering algorithm on noisy SAGE data. CONCLUSIONS: We explore the use of a seriation approach for visualization-based pattern discovery in SAGE data. Using both simulations and experimental data, we demonstrate that seriation is able to identify groups of co-expressed genes more accurately than a clustering algorithm developed specifically for SAGE data. Our results suggest that seriation is a useful method for the analysis of gene expression data whose applicability should be further pursued.

Expression of Groucho/TLE proteins during pancreas development.

BACKGROUND: The full-length mammalian homologs of groucho, Tle1, 2, 3, and 4, act as transcriptional corepressors and are recruited by transcription factors containing an eh1 or WRPW/Y domain. Many transcription factors critical to pancreas development contain a Gro/TLE interaction domain and several have been shown to require Gro/TLE interactions for proper function during neuronal development. However, a detailed analysis of the expression patterns of the Gro/TLE proteins in pancreas development has not been performed. Moreover, little is known about the ability of Gro/TLE proteins to interact with transcription factors in the pancreas. RESULTS: We describe the expression of Gro/TLE family members, and of 34 different transcription factors that contain a Gro/TLE interaction motif, in the pancreas utilizing nine SAGE libraries created from the developing and adult pancreas, as well as the GenePaint database. Next, we show the dynamic expression of Tle1, 2, 3, 4, 5 and 6 during pancreas development by qRT-PCR. To further define the cell-type specificity of the expression of these proteins we use immunofluorescence to co-localize them with Pdx1 at embryonic day 12.5 (E12.5), Ngn3 at E14.5, Pdx1, Nkx2-2, Insulin, Glucagon, Pancreatic polypeptide and Somatostatin at E18.5, as well as Insulin and Glucagon in the adult. We then show that Tle2 can interact with Nkx2-2, Hes1, Arx, and Nkx6-1 which are all critical factors in pancreas development. Finally, we demonstrate that Tle2 modulates the repressive abilities of Arx in a beta-cell line. CONCLUSION: Although Tle1, 2, 3, and 4 show overlapping expression in pancreatic progenitors and in the adult islet, the expression of these factors is restricted to different cell types during endocrine cell maturation. Of note, Tle2 and Tle3 are co-expressed with Gro/TLE interaction domain containing transcription factors that are essential for endocrine pancreas development. We further demonstrate that Tle2 can interact with several of these factors and that Tle2 modulate Arx's repressive activity. Taken together our studies suggest that Gro/TLE proteins play a role in the repression of target genes during endocrine cell specification.

Identification of transcripts with enriched expression in the developing and adult pancreas.

BACKGROUND: Despite recent advances, the transcriptional hierarchy driving pancreas organogenesis remains largely unknown, in part due to the paucity of comprehensive analyses. To address this deficit we generated ten SAGE libraries from the developing murine pancreas spanning Theiler stages 17-26, making use of available Pdx1 enhanced green fluorescent protein (EGFP) and Neurog3 EGFP reporter strains, as well as tissue from adult islets and ducts. RESULTS: We used a specificity metric to identify 2,536 tags with pancreas-enriched expression compared to 195 other mouse SAGE libraries. We subsequently grouped co-expressed transcripts with differential expression during pancreas development using K-means clustering. We validated the clusters first using quantitative real time PCR and then by analyzing the Theiler stage 22 pancreas in situ hybridization staining patterns of over 600 of the identified genes using the GenePaint database. These were then categorized into one of the five expression domains within the developing pancreas. Based on these results we identified a cascade of transcriptional regulators expressed in the endocrine pancreas lineage and, from this, we developed a predictive regulatory network describing beta-cell development. CONCLUSION: Taken together, this work provides evidence that the SAGE libraries generated here are a valuable resource for continuing to elucidate the molecular mechanisms regulating pancreas development. Furthermore, our studies provide a comprehensive analysis of pancreas development, and insights into the regulatory networks driving this process are revealed.