modMine: flexible access to modENCODE data.

In an effort to comprehensively characterize the functional elements within the genomes of the important model organisms Drosophila melanogaster and Caenorhabditis elegans, the NHGRI model organism Encyclopaedia of DNA Elements (modENCODE) consortium has generated an enormous library of genomic data along with detailed, structured information on all aspects of the experiments. The modMine database (http://intermine.modencode.org) described here has been built by the modENCODE Data Coordination Center to allow the broader research community to (i) search for and download data sets of interest among the thousands generated by modENCODE; (ii) access the data in an integrated form together with non-modENCODE data sets; and (iii) facilitate fine-grained analysis of the above data. The sophisticated search features are possible because of the collection of extensive experimental metadata by the consortium. Interfaces are provided to allow both biologists and bioinformaticians to exploit these rich modENCODE data sets now available via modMine.

Oncogenic ETS fusions promote DNA damage and proinflammatory responses via pericentromeric RNAs in extracellular vesicles.

Aberrant expression of the E26 transformation-specific (ETS) transcription factors characterizes numerous human malignancies. Many of these proteins, including EWS:FLI1 and EWS:ERG fusions in Ewing sarcoma (EwS) and TMPRSS2:ERG in prostate cancer (PCa), drive oncogenic programs via binding to GGAA repeats. We report here that both EWS:FLI1 and ERG bind and transcriptionally activate GGAA-rich pericentromeric heterochromatin. The respective pathogen-like HSAT2 and HSAT3 RNAs, together with LINE, SINE, ERV, and other repeat transcripts, are expressed in EwS and PCa tumors, secreted in extracellular vesicles (EVs), and are highly elevated in plasma of patients with EwS with metastatic disease. High human satellite 2 and 3 (HSAT2,3) levels in EWS:FLI1- or ERG-expressing cells and tumors were associated with induction of G2/M checkpoint, mitotic spindle, and DNA damage programs. These programs were also activated in EwS EV-treated fibroblasts, coincident with accumulation of HSAT2,3 RNAs, proinflammatory responses, mitotic defects, and senescence. Mechanistically, HSAT2,3-enriched cancer EVs induced cGAS-TBK1 innate immune signaling and formation of cytosolic granules positive for double-strand RNAs, RNA-DNA, and cGAS. Hence, aberrantly expressed ETS proteins derepress pericentromeric heterochromatin, yielding pathogenic RNAs that transmit genotoxic stress and inflammation to local and distant sites. Monitoring HSAT2,3 plasma levels and preventing their dissemination may thus improve therapeutic strategies and blood-based diagnostics.

Cloud-based uniform ChIP-Seq processing tools for modENCODE and ENCODE.

BACKGROUND: Funded by the National Institutes of Health (NIH), the aim of the Model Organism ENCyclopedia of DNA Elements (modENCODE) project is to provide the biological research community with a comprehensive encyclopedia of functional genomic elements for both model organisms C. elegans (worm) and D. melanogaster (fly). With a total size of just under 10 terabytes of data collected and released to the public, one of the challenges faced by researchers is to extract biologically meaningful knowledge from this large data set. While the basic quality control, pre-processing, and analysis of the data has already been performed by members of the modENCODE consortium, many researchers will wish to reinterpret the data set using modifications and enhancements of the original protocols, or combine modENCODE data with other data sets. Unfortunately this can be a time consuming and logistically challenging proposition. RESULTS: In recognition of this challenge, the modENCODE DCC has released uniform computing resources for analyzing modENCODE data on Galaxy (https://github.com/modENCODE-DCC/Galaxy), on the public Amazon Cloud (http://aws.amazon.com), and on the private Bionimbus Cloud for genomic research (http://www.bionimbus.org). In particular, we have released Galaxy workflows for interpreting ChIP-seq data which use the same quality control (QC) and peak calling standards adopted by the modENCODE and ENCODE communities. For convenience of use, we have created Amazon and Bionimbus Cloud machine images containing Galaxy along with all the modENCODE data, software and other dependencies. CONCLUSIONS: Using these resources provides a framework for running consistent and reproducible analyses on modENCODE data, ultimately allowing researchers to use more of their time using modENCODE data, and less time moving it around.

Deep SAGE analysis of the Caenorhabditis elegans transcriptome.

We employed the Tag-seq technique to generate global transcription profiles for different strains and life stages of the nematode C. elegans. Tag-seq generates cDNA tags as does Serial Analysis of Gene Expression (SAGE), but the method yields a much larger number of tags, generating much larger data sets than SAGE. We examined differences in the performance of SAGE and Tag-seq by comparing gene expression data for 13 pairs of libraries. We identified genes for which expression was consistently changed in long-lived worms. Additional genes emerged in the deeper Tag-seq profiles, including several 'signature' genes found among those zup-regulated in long-lived dauer larvae (cki-1, aak-2 and daf-16). Fifty to sixty percent of the genes differentially expressed in daf-2(-) versus daf-2(+) adults had fragmentary or no functional annotation, suggesting the involvement of as yet unstudied pathways in aging. We were able to distinguish between changes in gene expression associated with altered genotype or altered growth conditions. We found 62 cases of possible mRNA isoform switching in the 13 Tag-seq libraries, whereas the 13 SAGE libraries allowed detection of only 15 such occurrences. We observed strong expression of anti-sense transcripts for several mitochondrial genes, but nuclear anti-sense transcripts were neither abundant nor consistently expressed among the libraries.

Ewing Sarcoma-Derived Extracellular Vesicles Impair Dendritic Cell Maturation and Function.

Ewing sarcoma (EwS) is an aggressive pediatric cancer of bone and soft tissues characterized by scant T cell infiltration and predominance of immunosuppressive myeloid cells. Given the important roles of extracellular vesicles (EVs) in cancer-host crosstalk, we hypothesized that EVs secreted by EwS tumors target myeloid cells and promote immunosuppressive phenotypes. Here, EVs were purified from EwS and fibroblast cell lines and exhibited characteristics of small EVs, including size (100-170 nm) and exosome markers CD63, CD81, and TSG101. Treatment of healthy donor-derived CD33+ and CD14+ myeloid cells with EwS EVs but not with fibroblast EVs induced pro-inflammatory cytokine release, including IL-6, IL-8, and TNF. Furthermore, EwS EVs impaired differentiation of these cells towards monocytic-derived dendritic cells (moDCs), as evidenced by reduced expression of co-stimulatory molecules CD80, CD86 and HLA-DR. Whole transcriptome analysis revealed activation of gene expression programs associated with immunosuppressive phenotypes and pro-inflammatory responses. Functionally, moDCs differentiated in the presence of EwS EVs inhibited CD4+ and CD8+ T cell proliferation as well as IFNγ release, while inducing secretion of IL-10 and IL-6. Therefore, EwS EVs may promote a local and systemic pro-inflammatory environment and weaken adaptive immunity by impairing the differentiation and function of antigen-presenting cells.

ELT-2 is the predominant transcription factor controlling differentiation and function of the C. elegans intestine, from embryo to adult.

Starting with SAGE-libraries prepared from C. elegans FAC-sorted embryonic intestine cells (8E-16E cell stage), from total embryos and from purified oocytes, and taking advantage of the NextDB in situ hybridization data base, we define sets of genes highly expressed from the zygotic genome, and expressed either exclusively or preferentially in the embryonic intestine or in the intestine of newly hatched larvae; we had previously defined a similarly expressed set of genes from the adult intestine. We show that an extended TGATAA-like sequence is essentially the only candidate for a cis-acting regulatory motif common to intestine genes expressed at all stages. This sequence is a strong ELT-2 binding site and matches the sequence of GATA-like sites found to be important for the expression of every intestinal gene so far analyzed experimentally. We show that the majority of these three sets of highly expressed intestinal-specific/intestinal-enriched genes respond strongly to ectopic expression of ELT-2 within the embryo. By flow-sorting elt-2(null) larvae from elt-2(+) larvae and then preparing Solexa/Illumina-SAGE libraries, we show that the majority of these genes also respond strongly to loss-of-function of ELT-2. To test the consequences of loss of other transcription factors identified in the embryonic intestine, we develop a strain of worms that is RNAi-sensitive only in the intestine; however, we are unable (with one possible exception) to identify any other transcription factor whose intestinal loss-of-function causes a phenotype of comparable severity to the phenotype caused by loss of ELT-2. Overall, our results support a model in which ELT-2 is the predominant transcription factor in the post-specification C. elegans intestine and participates directly in the transcriptional regulation of the majority (>80%) of intestinal genes. We present evidence that ELT-2 plays a central role in most aspects of C. elegans intestinal physiology: establishing the structure of the enterocyte, regulating enzymes and transporters involved in digestion and nutrition, responding to environmental toxins and pathogenic infections, and regulating the downstream intestinal components of the daf-2/daf-16 pathway influencing aging and longevity.

Akt-mediated YB-1 phosphorylation activates translation of silent mRNA species.

YB-1 is a broad-specificity RNA-binding protein that is involved in regulation of mRNA transcription, splicing, translation, and stability. In both germinal and somatic cells, YB-1 and related proteins are major components of translationally inactive messenger ribonucleoprotein particles (mRNPs) and are mainly responsible for storage of mRNAs in a silent state. However, mechanisms regulating the repressor activity of YB-1 are not well understood. Here we demonstrate that association of YB-1 with the capped 5' terminus of the mRNA is regulated via phosphorylation by the serine/threonine protein kinase Akt. In contrast to its nonphosphorylated form, phosphorylated YB-1 fails to inhibit cap-dependent but not internal ribosome entry site-dependent translation of a reporter mRNA in vitro. We also show that similar to YB-1, Akt is associated with inactive mRNPs and that activated Akt may relieve translational repression of the YB-1-bound mRNAs. Using Affymetrix microarrays, we found that many of the YB-1-associated messages encode stress- and growth-related proteins, raising the intriguing possibility that Akt-mediated YB-1 phosphorylation could, in part, increase production of proteins regulating cell proliferation, oncogenic transformation, and stress response.

Discovery of novel alternatively spliced C. elegans transcripts by computational analysis of SAGE data.

BACKGROUND: Alternative RNA splicing allows cells to produce multiple protein isoforms from one gene. These isoforms may have specialized functions, and may be tissue- or stage-specific. Our aim was to use computational analysis of SAGE and genomic data to predict alternatively spliced transcripts expressed in C. elegans. RESULTS: We predicted novel alternatively spliced variants and confirmed five of eighteen candidates selected for experimental validation by RT-PCR tests and DNA sequencing. CONCLUSION: We show that SAGE data can be efficiently used to discover alternative mRNA isoforms, including those with skipped exons or retained introns. Our results also imply that C. elegans may produce a larger number of alternatively spliced transcripts than initially estimated.

Genes that may modulate longevity in C. elegans in both dauer larvae and long-lived daf-2 adults.

We used Serial Analysis of Gene Expression (SAGE) to compare the global transcription profiles of long-lived mutant daf-2 adults and dauer larvae, aiming to identify aging-related genes based on similarity of expression patterns. Genes that are expressed similarly in both long-lived types potentially define a common life-extending program. Comparison of eight SAGE libraries yielded a set of 120 genes, the expression of which was significantly different in long-lived worms vs. normal adults. The gene annotations indicate a strong link between oxidative stress and life span, further supporting the hypothesis that metabolic activity is a major determinant in longevity. The SAGE data show changes in mRNA levels for electron transport chain components, elevated expression of glyoxylate shunt enzymes and significantly reduced expression for components of the TCA cycle in longer-lived nematodes. We propose a model for enhanced longevity through a cytochrome c oxidase-mediated reduction in reactive oxygen species commonly held to be a major contributor to aging.

An interactive tool for visualization of relationships between gene expression profiles.

BACKGROUND: Application of phenetic methods to gene expression analysis proved to be a successful approach. Visualizing the results in a 3-dimentional space may further enhance these techniques. RESULTS: We designed and built TreeBuilder3D, an interactive viewer for visualizing the hierarchical relationships between expression profiles such as SAGE libraries or microarrays. The program allows loading expression data as plain text files and visualizing the relative differences of the analyzed datasets in 3-dimensional space using various distance metrics. CONCLUSION: TreeBuilder3D provides a simple interface and has a small size. Written in Java, TreeBuilder3D is a platform-independent, open source application, which may be useful in analysis of large-scale gene expression data.

The modENCODE Data Coordination Center: lessons in harvesting comprehensive experimental details.

The model organism Encyclopedia of DNA Elements (modENCODE) project is a National Human Genome Research Institute (NHGRI) initiative designed to characterize the genomes of Drosophila melanogaster and Caenorhabditis elegans. A Data Coordination Center (DCC) was created to collect, store and catalog modENCODE data. An effective DCC must gather, organize and provide all primary, interpreted and analyzed data, and ensure the community is supplied with the knowledge of the experimental conditions, protocols and verification checks used to generate each primary data set. We present here the design principles of the modENCODE DCC, and describe the ramifications of collecting thorough and deep metadata for describing experiments, including the use of a wiki for capturing protocol and reagent information, and the BIR-TAB specification for linking biological samples to experimental results. modENCODE data can be found at http://www.modencode.org.

Translational activation of snail1 and other developmentally regulated transcription factors by YB-1 promotes an epithelial-mesenchymal transition.

Increased expression of the transcription/translation regulatory protein Y-box binding protein-1 (YB-1) is associated with cancer aggressiveness, particularly in breast carcinoma. Here we establish that YB-1 levels are elevated in invasive breast cancer cells and correlate with reduced expression of E-cadherin and poor patient survival. Enforced expression of YB-1 in noninvasive breast epithelial cells induced an epithelial-mesenchymal transition (EMT) accompanied by enhanced metastatic potential and reduced proliferation rates. YB-1 directly activates cap-independent translation of messenger RNAs encoding Snail1 and other transcription factors implicated in downregulation of epithelial and growth-related genes and activation of mesenchymal genes. Hence, translational regulation by YB-1 is a restriction point enabling coordinated expression of a network of EMT-inducing transcription factors, likely acting together to promote metastatic spread.