The analysis of CAGE (Cap Analysis of Gene Expression) time-course has been proposed by the FANTOM5 Consortium to extend the understanding of the sequence of events facilitating cell state transition at the level of promoter regulation. To identify the most prominent transcriptional regulations induced by growth factors in human breast cancer, we apply here the Complexity Invariant Dynamic Time Warping motif EnRichment (CIDER) analysis approach to the CAGE time-course datasets of MCF-7 cells stimulated by epidermal growth factor (EGF) or heregulin (HRG). We identify a multi-level cascade of regulations rooted by the Serum Response Factor (SRF) transcription factor, connecting the MAPK-mediated transduction of the HRG stimulus to the negative regulation of the MAPK pathway by the members of the DUSP family phosphatases. The finding confirms the known primary role of FOS and FOSL1, members of AP-1 family, in shaping gene expression in response to HRG induction. Moreover, we identify a new potential regulation of DUSP5 and RARA (known to antagonize the transcriptional regulation induced by the estrogen receptors) by the activity of the AP-1 complex, specific to HRG response. The results indicate that a divergence in AP-1 regulation determines cellular changes of breast cancer cells stimulated by ErbB receptors.

Although it is generally accepted that cellular differentiation requires changes to transcriptional networks, dynamic regulation of promoters and enhancers at specific sets of genes has not been previously studied en masse. Exploiting the fact that active promoters and enhancers are transcribed, we simultaneously measured their activity in 19 human and 14 mouse time courses covering a wide range of cell types and biological stimuli. Enhancer RNAs, then messenger RNAs encoding transcription factors, dominated the earliest responses. Binding sites for key lineage transcription factors were simultaneously overrepresented in enhancers and promoters active in each cellular system. Our data support a highly generalizable model in which enhancer transcription is the earliest event in successive waves of transcriptional change during cellular differentiation or activation.

MicroRNAs are small non-coding RNAs that inhibit the translation of target mRNAs. In humans, most microRNAs are transcribed by RNA polymerase II as long primary transcripts and processed by sequential cleavage of the two RNase III enzymes, DROSHA and DICER, into precursor and mature microRNAs, respectively. Although the fundamental functions of microRNAs in RNA silencing have been gradually uncovered, less is known about the regulatory mechanisms of microRNA expression. Here, we report that telomerase reverse transcriptase (TERT) extensively affects the expression levels of mature microRNAs. Deep sequencing-based screens of short RNA populations revealed that the suppression of TERT resulted in the downregulation of microRNAs expressed in THP-1 cells and HeLa cells. Primary and precursor microRNA levels were also reduced under the suppression of TERT. Similar results were obtained with the suppression of either BRG1 (also called SMARCA4) or nucleostemin, which are proteins interacting with TERT and functioning beyond telomeres. These results suggest that TERT regulates microRNAs at the very early phases in their biogenesis, presumably through non-telomerase mechanism(s).

Deciphering the most common modes by which chromatin regulates transcription, and how this is related to cellular status and processes is an important task for improving our understanding of human cellular biology. The FANTOM5 and ENCODE projects represent two independent large scale efforts to map regulatory and transcriptional features to the human genome. Here we investigate chromatin features around a comprehensive set of transcription start sites in four cell lines by integrating data from these two projects.

VEGF-C/VEGFR-3 signaling plays a central role in lymphatic development, regulating the budding of lymphatic progenitor cells from embryonic veins and maintaining the expression of PROX1 during later developmental stages. However, how VEGFR-3 activation translates into target gene expression is still not completely understood. We used cap analysis of gene expression (CAGE) RNA sequencing to characterize the transcriptional changes invoked by VEGF-C in LECs and to identify the transcription factors (TFs) involved. We found that MAFB, a TF involved in differentiation of various cell types, is rapidly induced and activated by VEGF-C. MAFB induced expression of PROX1 as well as other TFs and markers of differentiated LECs, indicating a role in the maintenance of the mature LEC phenotype. Correspondingly, E14.5 Mafb(-/-) embryos showed impaired lymphatic patterning in the skin. This suggests that MAFB is an important TF involved in lymphangiogenesis.

The Pax6 transcription factor is expressed in cerebellar granule cells and when mutated, as in the Sey/Sey mouse, produces granule cells with disturbed survival and migration and with defects in neurite extension. The impact of Pax6 on other genes in the context of cerebellar development has not been identified. In this study, we performed transcriptome comparisons between wildtype and Pax6-null whole cerebellar tissue at embryonic day (E) 13.5, 15.5 and 18.5 using Affymetrix arrays (U74Av2). Statistical analyses identified 136 differentially regulated transcripts (FDR 0.05, 1.2-fold change cutoff) over time in Pax6-null cerebellar tissue. In parallel we examined the Math1-null granuloprival cerebellum and identified 228 down-regulated transcripts (FDR 0.05, 1.2-fold change cutoff). The intersection of these two microarray datasets produced a total of 21 differentially regulated transcripts. For a subset of the identified transcripts, we used qRT-PCR to validate the microarray data and demonstrated the expression in the rhombic lip lineage and differential expression in Pax6-null cerebellum with in situ hybridisation analysis. The candidate genes identified in this way represent direct or indirect Pax6-downstream genes involved in cerebellar development.

The number of neural progenitor cells (NPCs) decreases with advancing age, and the mechanisms responsible for this decline is unclear. Here, we demonstrate the importance of genetics as a modulator for the age-related decline in NPCs. We systematically quantified the number of proliferating NPCs in the rostral migratory stream, the rostral extension of the subventricular zone, in 9 inbred mouse strains from 3 to 18 months of age. A striking negative impact of age and significant interstrain differences in the number of NPCs was detected at 3 and 12 months of age. Extended proliferative profiles of C57BL/6J and DBA/2J from 3 to 24 months of age revealed differential dynamics of the age-related decline in NPCs. Statistically significant interaction effects for age and strain were detected over the 3- to 7-month period. Strain differences were mapped to several genetic loci suggesting complex genetic control of NPC proliferation at different ages. Furthermore, correlational analyses revealed the differential regulation of cell proliferation and genes that may underlie the proliferative deficits of NPCs in the aging brain.

Evolutionarily conserved RFX transcription factors (TFs) regulate their target genes through a DNA sequence motif called the X-box. Thereby they regulate cellular specialization and terminal differentiation. Here, we provide a comprehensive analysis of all the eight human RFX genes (RFX1-8), their spatial and temporal expression profiles, potential upstream regulators and target genes.

Endomyocardial biopsy is a valuable tool in cardiac diagnostics but is limited by low diagnostic yield and significant complication risks. Meanwhile, recent developments in transcriptomic and proteomic technologies promise a wealth of biological data from minimal tissue samples. To take advantage of the minimal tissue amount needed for molecular analyses, we have developed a sub-millimeter endovascular biopsy device, considerably smaller than current clinical equipment, and devised a low-input RNA-sequencing protocol for analyzing small tissue samples. In in vivo evaluation in swine, 81% of biopsy attempts (n = 157) were successful. High quality RNA-sequencing data was generated from 91% of the sequenced cardiac micro-biopsy samples (n = 32). Gene expression signatures of samples taken with the novel device were comparable with a conventional device. No major complications were detected either during procedures or during 7 days' follow-up, despite acquiring a relatively large number of biopsies (median 30) in each animal. In conclusion, the novel device coupled with RNA-sequencing provides a feasible method to obtain molecular data from the myocardium. The method is less traumatic and has a higher flexibility compared to conventional methods, enabling safer and more targeted sampling from different parts of the myocardium.

DEAD-Box Helicase 3 X-Linked (DDX3X) is frequently mutated in the Wingless (WNT) and Sonic hedghog (SHH) subtypes of medulloblastoma-the commonest malignant childhood brain tumor, but whether DDX3X functions as a medulloblastoma oncogene or tumor suppressor gene is not known. Here, we show that Ddx3x regulates hindbrain patterning and development by controlling Hox gene expression and cell stress signaling. In mice predisposed to Wnt- or Shh medulloblastoma, Ddx3x sensed oncogenic stress and suppressed tumor formation. WNT and SHH medulloblastomas normally arise only in the lower and upper rhombic lips, respectively. Deletion of Ddx3x removed this lineage restriction, enabling both medulloblastoma subtypes to arise in either germinal zone. Thus, DDX3X is a medulloblastoma tumor suppressor that regulates hindbrain development and restricts the competence of cell lineages to form medulloblastoma subtypes.

Zebrafish, a popular organism for studying embryonic development and for modeling human diseases, has so far lacked a systematic functional annotation program akin to those in other animal models. To address this, we formed the international DANIO-CODE consortium and created a central repository to store and process zebrafish developmental functional genomic data. Our data coordination center ( https://danio-code.zfin.org ) combines a total of 1,802 sets of unpublished and re-analyzed published genomic data, which we used to improve existing annotations and show its utility in experimental design. We identified over 140,000 cis-regulatory elements throughout development, including classes with distinct features dependent on their activity in time and space. We delineated the distinct distance topology and chromatin features between regulatory elements active during zygotic genome activation and those active during organogenesis. Finally, we matched regulatory elements and epigenomic landscapes between zebrafish and mouse and predicted functional relationships between them beyond sequence similarity, thus extending the utility of zebrafish developmental genomics to mammals.

The cellular localization and development of neuronal intranuclear inclusions (NIIs) in cortex and striatum of R6/2 HD transgenic mice were studied to ascertain the relationship of NIIs to symptom formation in these mice and gain clues regarding the possible relationship of NII formation to neuropathology in Huntington's disease (HD). All NIIs observed in R6/2 mice were ubiquitinated, and no evidence was observed for a contribution to them from wild-type huntingtin; they were first observed in cortex and striatum at 3.5 weeks of age. In cortex, NIIs increased rapidly in size and prevalence after their appearance. Generally, cortical projection neurons developed NIIs more rapidly than cortical interneurons containing calbindin or parvalbumin. Few cortical somatostatinergic interneurons, however, formed NIIs. In striatum, calbindinergic projection neurons and parvalbuminergic interneurons rapidly formed NIIs, but they formed more gradually in cholinergic interneurons, and few somatostatinergic interneurons developed NIIs. Striatal NIIs tended to be smaller than those in cortex. The early accumulation of NIIs in cortex and striatum in R6/2 mice is consistent with the early appearance of motor and learning abnormalities in these mice, and the eventual pervasiveness of NIIs at ages at which severe abnormalities are evident is consistent with their contribution to a neuronal dysfunction underlying the abnormalities. That cortex develops larger NIIs than striatum, however, is inconsistent with the preferential loss of striatal neurons in HD but is consistent with recent evidence of early morphological abnormalities in cortical neurons in HD. That calbindinergic and parvalbuminergic striatal neurons develop large NIIs is consistent with a contribution of nuclear aggregate formation to their high degree of vulnerability in HD.

The pathway for RNA interference is widespread in metazoans and participates in numerous cellular tasks, from gene silencing to chromatin remodeling and protection against retrotransposition. The unicellular eukaryote Trypanosoma cruzi is missing the canonical RNAi pathway and is unable to induce RNAi-related processes. To further understand alternative RNA pathways operating in this organism, we have performed deep sequencing and genome-wide analyses of a size-fractioned cDNA library (16-61 nt) from the epimastigote life stage. Deep sequencing generated 582,243 short sequences of which 91% could be aligned with the genome sequence. About 95-98% of the aligned data (depending on the haplotype) corresponded to small RNAs derived from tRNAs, rRNAs, snRNAs and snoRNAs. The largest class consisted of tRNA-derived small RNAs which primarily originated from the 3' end of tRNAs, followed by small RNAs derived from rRNA. The remaining sequences revealed the presence of 92 novel transcribed loci, of which 79 did not show homology to known RNA classes.

PROM1 is the gene encoding prominin-1 or CD133, an important cell surface marker for the isolation of both normal and cancer stem cells. PROM1 transcripts initiate at a range of transcription start sites (TSS) associated with distinct tissue and cancer expression profiles. Using high resolution Cap Analysis of Gene Expression (CAGE) sequencing we characterize TSS utilization across a broad range of normal and developmental tissues. We identify a novel proximal promoter (P6) within CD133(+) melanoma cell lines and stem cells. Additional exon array sampling finds P6 to be active in populations enriched for mesenchyme, neural stem cells and within CD133(+) enriched Ewing sarcomas. The P6 promoter is enriched with respect to previously characterized PROM1 promoters for a HMGI/Y (HMGA1) family transcription factor binding site motif and exhibits different epigenetic modifications relative to the canonical promoter region of PROM1.

The cerebellum assists coordination of somatomotor, respiratory, and autonomic actions. Purkinje cell alterations or loss appear in sudden infant death and sudden death in epilepsy victims, possibly contributing to the fatal event. We evaluated breathing patterns in 12 wild-type (WT) and Lurcher mutant mice with 100% developmental cerebellar Purkinje cell loss under baseline (room air), and recovery from hypercapnia, a concern in sudden death events. Six mutant and six WT mice were exposed to 4-min blocks of increasing CO2 (2, 4, 6, and 8%), separated by 4-min recovery intervals in room air. Breath-by-breath patterns, including depth of breathing and end-expiratory pause (EEP) durations during recovery, were recorded. No baseline genotypic differences emerged. However, during recovery, EEP durations significantly lengthened in mutants, compared to WT mice, following the relatively low levels of CO2 exposure. Additionally, mutant mice exhibited signs of post-sigh disordered breathing during recovery following each exposure. Developmental cerebellar Purkinje cell loss significantly affects compensatory breathing patterns following mild CO2 exposure, possibly by inhibiting recovery from elevated CO2. These data implicate cerebellar Purkinje cells in the ability to recover from hypercarbia, suggesting that neuropathologic changes or loss of these cells contribute to inadequate ventilatory recovery to increased environmental CO2. Multiple disorders, including sudden infant death syndrome (SIDS) and sudden unexpected death in epilepsy (SUDEP), appear to involve both cardiorespiratory failure and loss or injury to cerebellar Purkinje cells; the findings support the concept that such neuropathology may precede and exert a prominent role in these fatal events.

Math1 is the defining molecule of the cerebellar rhombic lip and Pax6 is downstream in the Math1 pathway. In the present study, we discover that Wntless (Wls) is a novel molecular marker of the cells in the interior face of the rhombic lip throughout normal mouse cerebellar development. Wls expression is found complementary to the expression of Math1 and Pax6, which are localized to the exterior face of the rhombic lip. To determine the interaction between these molecules, we examine the loss-of-Math1 or loss-of-Pax6 in the cerebellum, i.e., the Math1-null and Pax6-null (Sey) mutant cerebella. The presence of Wls-positive cells in the Math1-null rhombic lip indicates that Wls expression is independent of Math1. In the Sey mutant cerebellum, there is an expansion of Wls-expressing cells into regions that are normally colonized by Pax6-expressing cells. The ectopic expression of Wls in the Pax6-null cerebellum suggests a negative interaction between Wls-expressing cells and Pax6-positive cells. These findings suggest that the rhombic lip is dynamically patterned by the expression of Wls, Math1, and Pax6. We also examine five rhombic lip cell markers (Wls, Math1, Pax6, Lmx1a, and Tbr2) to identify four molecularly distinct compartments in the rhombic lip during cerebellar development. The existence of spatial compartmentation in the rhombic lip and the interplay between Wls, Math1, and Pax6 in the rhombic lip provides novel views of early cerebellar development.

Choroid plexus carcinomas (CPCs) are poorly understood and frequently lethal brain tumors with few treatment options. Using a mouse model of the disease and a large cohort of human CPCs, we performed a cross-species, genome-wide search for oncogenes within syntenic regions of chromosome gain. TAF12, NFYC, and RAD54L co-located on human chromosome 1p32-35.3 and mouse chromosome 4qD1-D3 were identified as oncogenes that are gained in tumors in both species and required for disease initiation and progression. TAF12 and NFYC are transcription factors that regulate the epigenome, whereas RAD54L plays a central role in DNA repair. Our data identify a group of concurrently gained oncogenes that cooperate in the formation of CPC and reveal potential avenues for therapy.

Wntless (Wls) is implicated in the Wnt signaling pathway by regulating the secretion of Wnt molecules. During brain development, Wls is expressed in the isthmic organizer (ISO) and rhombic lip (RL). Wls regulates Wnt1 secretion at the ISO which is required to induce midbrain-hindbrain structures. However, Wls function in the RL is not known. Here, we employed Nestin-cre to delete Wls specifically in the RL during mid-gestation. The loss-of-Wls leads to an abnormal RL during development and cerebellar vermis hypoplasia at birth. The Wls conditional knockout (cKO) has rudimentary foliation with an absence of Bergmann glia fibers in the external germinal layer (EGL). The Wls-cKO cerebellum also exhibits ectopia of several cell types and aberrations in granule cell organization. Finally, there is a loss of 85% of unipolar brush cells. From these findings, Wls-expressing cells in the rhombic lip are implicated in the orchestration of cerebellar development.

The promoters of immediate early genes (IEGs) are rapidly activated in response to an external stimulus. These genes, also known as primary response genes, have been identified in a range of cell types, under diverse extracellular signals and using varying experimental protocols. Whereas genomic dissection on a case-by-case basis has not resulted in a comprehensive catalogue of IEGs, a rigorous meta-analysis of eight genome-wide FANTOM5 CAGE (cap analysis of gene expression) time course datasets reveals successive waves of promoter activation in IEGs, recapitulating known relationships between cell types and stimuli: we obtain a set of 57 (42 protein-coding) candidate IEGs possessing promoters that consistently drive a rapid but transient increase in expression over time. These genes show significant enrichment for known IEGs reported previously, pathways associated with the immediate early response, and include a number of non-coding RNAs with roles in proliferation and differentiation. Surprisingly, we also find strong conservation of the ordering of activation for these genes, such that 77 pairwise promoter activation orderings are conserved. Using the leverage of comprehensive CAGE time series data across cell types, we also document the extensive alternative promoter usage by such genes, which is likely to have been a barrier to their discovery until now. The common activation ordering of the core set of early-responding genes we identify may indicate conserved underlying regulatory mechanisms. By contrast, the considerably larger number of transiently activated genes that are specific to each cell type and stimulus illustrates the breadth of the primary response.

In obesity, white adipose tissue (WAT) inflammation is linked to insulin resistance. Increased adipocyte chemokine (C-C motif) ligand 2 (CCL2) secretion may initiate adipose inflammation by attracting the migration of inflammatory cells into the tissue. Using an unbiased approach, we identified adipose microRNAs (miRNAs) that are dysregulated in human obesity and assessed their possible role in controlling CCL2 production. In subcutaneous WAT obtained from 56 subjects, 11 miRNAs were present in all subjects and downregulated in obesity. Of these, 10 affected adipocyte CCL2 secretion in vitro and for 2 miRNAs (miR-126 and miR-193b), regulatory circuits were defined. While miR-126 bound directly to the 3'-untranslated region of CCL2 mRNA, miR-193b regulated CCL2 production indirectly through a network of transcription factors, many of which have been identified in other inflammatory conditions. In addition, overexpression of miR-193b and miR-126 in a human monocyte/macrophage cell line attenuated CCL2 production. The levels of the two miRNAs in subcutaneous WAT were significantly associated with CCL2 secretion (miR-193b) and expression of integrin, α-X, an inflammatory macrophage marker (miR-193b and miR-126). Taken together, our data suggest that miRNAs may be important regulators of adipose inflammation through their effects on CCL2 release from human adipocytes and macrophages.