A Latent Pro-survival Function for the Mir-290-295 Cluster in Mouse Embryonic Stem Cells

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of thousands of distinct mRNAs. While some regulatory interactions help to maintain basal cellular functions, others are likely relevant in more specific settings, such as response to stress. Here we describe such a role for the mir-290-295 cluster, the dominant miRNA cluster in mouse embryonic stem cells (mESCs). Examination of a target list generated from bioinformatic prediction, as well as expression data following miRNA loss, revealed strong enrichment for apoptotic regulators, two of which we validated directly: Caspase 2, the most highly conserved mammalian caspase, and Ei24, a p53 transcriptional target. Consistent with these predictions, mESCs lacking miRNAs were more likely to initiate apoptosis following genotoxic exposure to gamma irradiation or doxorubicin. Knockdown of either candidate partially rescued this pro-apoptotic phenotype, as did transfection of members of the mir-290-295 cluster. These findings were recapitulated in a specific mir-290-295 deletion line, confirming that they reflect miRNA functions at physiological levels. In contrast to the basal regulatory roles previously identified, the pro-survival phenotype shown here may be most relevant to stressful gestations, where pro-oxidant metabolic states induce DNA damage. Similarly, this cluster may mediate chemotherapeutic resistance in a neoplastic context, making it a useful clinical target.National Institutes of Health (U.S.) (NIH grant RO1-GM34277)National Cancer Institute (U.S.) (NCI grant PO1-CA42063)National Cancer Institute (U.S.) (NCI Cancer Center Support (core) grant P30-CA14051

Transcription forms and remodels supercoiling domains unfolding large-scale chromatin structures

DNA supercoiling is an inherent consequence of twisting DNA and is critical for regulating gene expression and DNA replication. However, DNA supercoiling at a genomic scale in human cells is uncharacterized. To map supercoiling we used biotinylated-trimethylpsoralen as a DNA structure probe to show the genome is organized into supercoiling domains. Domains are formed and remodeled by RNA polymerase and topoisomerase activities and are flanked by GC-AT boundaries and CTCF binding sites. Under-wound domains are transcriptionally active, enriched in topoisomerase I, “open” chromatin fibers and DNaseI sites, but are depleted of topoisomerase II. Furthermore DNA supercoiling impacts on additional levels of chromatin compaction as under-wound domains are cytologically decondensed, topologically constrained, and decompacted by transcription of short RNAs. We suggest that supercoiling domains create a topological environment that facilitates gene activation providing an evolutionary purpose for clustering genes along chromosomes

Extensive and coordinated transcription of noncoding RNAs within cell-cycle promoters

Transcription of long noncoding RNAs (lncRNAs) within gene regulatory elements can modulate gene activity in response to external stimuli, but the scope and functions of such activity are not known. Here we use an ultrahigh-density array that tiles the promoters of 56 cell-cycle genes to interrogate 108 samples representing diverse perturbations. We identify 216 transcribed regions that encode putative lncRNAs, many with RT-PCR–validated periodic expression during the cell cycle, show altered expression in human cancers and are regulated in expression by specific oncogenic stimuli, stem cell differentiation or DNA damage. DNA damage induces five lncRNAs from the CDKN1A promoter, and one such lncRNA, named PANDA, is induced in a p53-dependent manner. PANDA interacts with the transcription factor NF-YA to limit expression of pro-apoptotic genes; PANDA depletion markedly sensitized human fibroblasts to apoptosis by doxorubicin. These findings suggest potentially widespread roles for promoter lncRNAs in cell-growth control.National Institutes of Health (U.S.)National Institute of Arthritis and Musculoskeletal and Skin Diseases (U.S.) (NIAMS) (K08-AR054615))National Cancer Institute (U.S.) (NIH/(NCI) (R01-CA118750))National Cancer Institute (U.S.) (NIH/(NCI) R01-CA130795))Juvenile Diabetes Research Foundation InternationalAmerican Cancer SocietyHoward Hughes Medical Institute (Early career scientist)Stanford University (Graduate Fellowship)National Science Foundation (U.S.) (Graduate Research Fellowship)United States. Dept. of Defense (National Defense Science and Engineering Graduate Fellowship

Tumor Suppressor RASSF1A Promoter: p53 Binding and Methylation

Oncogenes and tumor suppressors work in concert to regulate cell growth or death, which is a pair of antagonist factors for regulation of tumorigenesis. Here we show promoter characteristic of tumor suppressor RASSF1A, which revealed a p53 binding site in the distal and a GC-rich region in the proximal promoter region of RASSF1A, in despite of TATA box-less. The GC-rich region, which is ∼300 bp upstream from the RASSF1A ATG, showed the strongest promoter activity in an assay of RASSF1A-driving GFP expression. Methylation analysis of the CpG island showed that 78.57% of the GC sties were methylated in testis tumor samples compared with methylation-less in normal testis. Hypermethylation of the GC-rich region is associated with RASSF1A silencing in human testis tumors. In addition, electrophoretic mobility shift assay indicated that p53 protein bound to the RASSF1A promoter. Further chromatin immunoprecipitation confirmed p53 binding to the RASSF1A. Moreover, p53 binding to the promoter down-regulated RASSF1A expression. These results suggest that p53 protein specifically binds to the RASSF1A promoter and inhibits its expression. Our results provide new insight into the mechanism of action of tumor suppressors and may be a starting point for development of new approaches to cancer treatment

Integrating 5-Hydroxymethylcytosine into the Epigenomic Landscape of Human Embryonic Stem Cells

Covalent modification of DNA distinguishes cellular identities and is crucial for regulating the pluripotency and differentiation of embryonic stem (ES) cells. The recent demonstration that 5-methylcytosine (5-mC) may be further modified to 5-hydroxymethylcytosine (5-hmC) in ES cells has revealed a novel regulatory paradigm to modulate the epigenetic landscape of pluripotency. To understand the role of 5-hmC in the epigenomic landscape of pluripotent cells, here we profile the genome-wide 5-hmC distribution and correlate it with the genomic profiles of 11 diverse histone modifications and six transcription factors in human ES cells. By integrating genomic 5-hmC signals with maps of histone enrichment, we link particular pluripotency-associated chromatin contexts with 5-hmC. Intriguingly, through additional correlations with defined chromatin signatures at promoter and enhancer subtypes, we show distinct enrichment of 5-hmC at enhancers marked with H3K4me1 and H3K27ac. These results suggest potential role(s) for 5-hmC in the regulation of specific promoters and enhancers. In addition, our results provide a detailed epigenomic map of 5-hmC from which to pursue future functional studies on the diverse regulatory roles associated with 5-hmC

Evidence That Gene Activation and Silencing during Stem Cell Differentiation Requires a Transcriptionally Paused Intermediate State

A surprising portion of both mammalian and Drosophila genomes are transcriptionally paused, undergoing initiation without elongation. We tested the hypothesis that transcriptional pausing is an obligate transition state between definitive activation and silencing as human embryonic stem cells (hESCs) change state from pluripotency to mesoderm. Chromatin immunoprecipitation for trimethyl lysine 4 on histone H3 (ChIP-Chip) was used to analyze transcriptional initiation, and 3′ transcript arrays were used to determine transcript elongation. Pluripotent and mesodermal cells had equivalent fractions of the genome in active and paused transcriptional states (∼48% each), with ∼4% definitively silenced (neither initiation nor elongation). Differentiation to mesoderm changed the transcriptional state of 12% of the genome, with roughly equal numbers of genes moving toward activation or silencing. Interestingly, almost all loci (98–99%) changing transcriptional state do so either by entering or exiting the paused state. A majority of these transitions involve either loss of initiation, as genes specifying alternate lineages are archived, or gain of initiation, in anticipation of future full-length expression. The addition of chromatin dynamics permitted much earlier predictions of final cell fate compared to sole use of conventional transcript arrays. These findings indicate that the paused state may be the major transition state for genes changing expression during differentiation, and implicate control of transcriptional elongation as a key checkpoint in lineage specification

The Prevalence and Regulation of Antisense Transcripts in Schizosaccharomyces pombe

A strand-specific transcriptome sequencing strategy, directional ligation sequencing or DeLi-seq, was employed to profile antisense transcriptome of Schizosaccharomyces pombe. Under both normal and heat shock conditions, we found that polyadenylated antisense transcripts are broadly expressed while distinct expression patterns were observed for protein-coding and non-coding loci. Dominant antisense expression is enriched in protein-coding genes involved in meiosis or stress response pathways. Detailed analyses further suggest that antisense transcripts are independently regulated with respect to their sense transcripts, and diverse mechanisms might be potentially involved in the biogenesis and degradation of antisense RNAs. Taken together, antisense transcription may have profound impacts on global gene regulation in S. pombe

Detection and Removal of Biases in the Analysis of Next-Generation Sequencing Reads

Since the emergence of next-generation sequencing (NGS) technologies, great effort has been put into the development of tools for analysis of the short reads. In parallel, knowledge is increasing regarding biases inherent in these technologies. Here we discuss four different biases we encountered while analyzing various Illumina datasets. These biases are due to both biological and statistical effects that in particular affect comparisons between different genomic regions. Specifically, we encountered biases pertaining to the distributions of nucleotides across sequencing cycles, to mappability, to contamination of pre-mRNA with mRNA, and to non-uniform hydrolysis of RNA. Most of these biases are not specific to one analyzed dataset, but are present across a variety of datasets and within a variety of genomic contexts. Importantly, some of these biases correlated in a highly significant manner with biological features, including transcript length, gene expression levels, conservation levels, and exon-intron architecture, misleadingly increasing the credibility of results due to them. We also demonstrate the relevance of these biases in the context of analyzing an NGS dataset mapping transcriptionally engaged RNA polymerase II (RNAPII) in the context of exon-intron architecture, and show that elimination of these biases is crucial for avoiding erroneous interpretation of the data. Collectively, our results highlight several important pitfalls, challenges and approaches in the analysis of NGS reads