Despite growing consensus that long intergenic non-coding ribonucleic acids (lincRNAs) are modulators of cancer, the knowledge about the deoxyribonucleic acid (DNA) methylation patterns of lincRNAs in cancers remains limited. In this study, we constructed DNA methylation profiles for 4629 tumors and 705 normal tissue samples from 20 different types of human cancer by reannotating data of DNA methylation arrays. We found that lincRNAs had different promoter methylation patterns in cancers. We classified 2461 lincRNAs into two categories and three subcategories, according to their promoter methylation patterns in tumors. LincRNAs with resistant methylation patterns in tumors had conserved transcriptional regulation regions and were ubiquitously expressed across normal tissues. By integrating cancer subtype data and patient clinical information, we identified lincRNAs with promoter methylation patterns that were associated with cancer status, subtype or prognosis for several cancers. Network analysis of aberrantly methylated lincRNAs in cancers showed that lincRNAs with aberrant methylation patterns might be involved in cancer development and progression. The methylated and demethylated lincRNAs identified in this study provide novel insights for developing cancer biomarkers and potential therapeutic targets.

Tauopathies such as Alzheimer's Disease (AD) are neurodegenerative disorders for which there is presently no cure. They are named after the abnormal oligomerization/aggregation of the neuronal microtubule-associated Tau protein. Besides its role as a microtubule-associated protein, a DNA-binding capacity and a nuclear localization for Tau protein has been described in neurons. While questioning the potential role of Tau-DNA binding in the development of tauopathies, we have carried out a large-scale analysis of the interaction of Tau protein with the neuronal genome under physiological and heat stress conditions using the ChIP-on-chip technique that combines Chromatin ImmunoPrecipitation (ChIP) with DNA microarray (chip). Our findings show that Tau protein specifically interacts with genic and intergenic DNA sequences of primary culture of neurons with a preference for DNA regions positioned beyond the ±5000 bp range from transcription start site. An AG-rich DNA motif was found recurrently present within Tau-interacting regions and 30% of Tau-interacting regions overlapped DNA sequences coding for lncRNAs. Neurological processes affected in AD were enriched among Tau-interacting regions with in vivo gene expression assays being indicative of a transcriptional repressor role for Tau protein, which was exacerbated in neurons displaying nuclear pathological oligomerized forms of Tau protein.

The 3' untranslated regions (3' UTRs) of transcripts serve as important hubs for posttranscriptional gene expression regulation. Here, we find that the exonisation of intergenic Alu elements introduced new terminal exons and polyadenylation sites during human genome evolution. While Alu exonisation from introns has been described previously, we shed light on a novel mechanism to create alternative 3' UTRs, thereby opening opportunities for differential posttranscriptional regulation. On the mechanistic level, we show that intergenic Alu exonisation can compete both with alternative splicing and polyadenylation in the upstream gene. Notably, the Alu-derived isoforms are often expressed in a tissue-specific manner, and the Alu-derived 3' UTRs can alter mRNA stability. In summary, we demonstrate that intergenic elements can affect processing of preceding genes, and elucidate how intergenic Alu exonisation can contribute to tissue-specific posttranscriptional regulation by expanding the repertoire of 3' UTRs.

Thousands of disease-associated SNPs (daSNPs) are located in intergenic regions (IGR), making it difficult to understand their association with disease phenotypes. Recent analysis found that non-coding daSNPs were frequently located in or approximate to regulatory elements, inspiring us to try to explain the disease phenotypes of IGR daSNPs through nearby regulatory sequences. Hence, after locating the nearest distal regulatory element (DRE) to a given IGR daSNP, we applied a computational method named INTREPID to predict the target genes regulated by the DRE, and then investigated their functional relevance to the IGR daSNP's disease phenotypes. 36.8% of all IGR daSNP-disease phenotype associations investigated were possibly explainable through the predicted target genes, which were enriched with, were functionally relevant to, or consisted of the corresponding disease genes. This proportion could be further increased to 60.5% if the LD SNPs of daSNPs were also considered. Furthermore, the predicted SNP-target gene pairs were enriched with known eQTL/mQTL SNP-gene relationships. Overall, it's likely that IGR daSNPs may contribute to disease phenotypes by interfering with the regulatory function of their nearby DREs and causing abnormal expression of disease genes.

The epigenetic landscape of cancer includes both focal hypermethylation and broader hypomethylation in a genome-wide manner. By means of a comprehensive genomic analysis on 6637 tissues of 21 tumor types, we here show that the degrees of overall methylation in CpG island (CGI) and demethylation in intergenic regions, defined as 'backbone', largely vary among different tumors. Depending on tumor type, both CGI methylation and backbone demethylation are often associated with clinical, epidemiological and biological features such as age, sex, smoking history, anatomic location, histological type and grade, stage, molecular subtype and biological pathways. We found connections between CGI methylation and hypermutability, microsatellite instability, IDH1 mutation, 19p gain and polycomb features, and backbone demethylation with chromosomal instability, NSD1 and TP53 mutations, 5q and 19p loss and long repressive domains. These broad epigenetic patterns add a new dimension to our understanding of tumor biology and its clinical implications.

How many intergenically encoded non-coding RNAs (ncRNAs) are expressed during various developmental stages in Drosophila? Previous analyses in one or a few developmental stages indicated abundant expression of intergenic ncRNAs. However, some reported that ncRNAs have been recently falsified, and, in general, the false positive rate for ncRNA detection is unknown. In this report, we used reverse transcription-PCR (RT-PCR), a more robust method, to detect ncRNAs from the intergenic regions that are expressed in four major developmental stages (6-8 h embryo, 20-22 h embryo, larvae and adult). We tested 1027 regions, approximately 10% of all intergenic regions, and detected transcription by RT-PCR. The results from 18 342 RT-PCR experiments revealed evidence for transcription in 72.7% of intergenic regions in the developmental process. The early developmental stage appears to be associated with more abundant ncRNAs than later developmental stages. In the early stage, we detected 43.6% of intergenic regions that encode transcripts in the triplicate RT-PCR experiments, yielding an estimate of 5006 intergenic regions in the entire genome likely encoding ncRNAs. We compared the RT-PCR-related approach with previous tiling array-based approach and observed that the latter method is insensitive to short ncRNAs, especially the molecules less than 120 bp. We measured false positive rates for the analyzed genomic approaches including the RT-PCR and tiling array method.

The discovery of structured non-coding RNAs (ncRNAs) in bacteria can reveal new facets of biology and biochemistry. Comparative genomics analyses executed by powerful computer algorithms have successfully been used to uncover many novel bacterial ncRNA classes in recent years. However, this general search strategy favors the discovery of more common ncRNA classes, whereas progressively rarer classes are correspondingly more difficult to identify. In the current study, we confront this problem by devising several methods to select subsets of intergenic regions that can concentrate these rare RNA classes, thereby increasing the probability that comparative sequence analysis approaches will reveal their existence. By implementing these methods, we discovered 224 novel ncRNA classes, which include ROOL RNA, an RNA class averaging 581 nt and present in multiple phyla, several highly conserved and widespread ncRNA classes with properties that suggest sophisticated biochemical functions and a multitude of putative cis-regulatory RNA classes involved in a variety of biological processes. We expect that further research on these newly found RNA classes will reveal additional aspects of novel biology, and allow for greater insights into the biochemistry performed by ncRNAs.

The Escherichia coli RutR protein is the master regulator of genes involved in pyrimidine catabolism. Here we have used chromatin immunoprecipitation in combination with DNA microarrays to measure the binding of RutR across the chromosome of exponentially growing E. coli cells. Twenty RutR-binding targets were identified and analysis of these targets generated a DNA consensus logo for RutR binding. Complementary in vitro binding assays showed high-affinity RutR binding to 16 of the 20 targets, with the four low-affinity RutR targets lacking predicted key binding determinants. Surprisingly, most of the DNA targets for RutR are located within coding segments of the genome and appear to have little or no effect on transcript levels in the conditions tested. This contrasts sharply with other E. coli transcription factors whose binding sites are primarily located in intergenic regions. We suggest that either RutR has yet undiscovered function or that evolution has been slow to eliminate non-functional DNA sites for RutR because they do not have an adverse effect on cell fitness.

Missense point mutations in the TP53 gene are frequent genetic alterations in human tumor tissue and cell lines derived thereof. Mutant p53 (mutp53) proteins have lost sequence-specific DNA binding, but have retained the ability to interact in a structure-selective manner with non-B DNA and to act as regulators of transcription. To identify functional binding sites of mutp53, we established a small library of genomic sequences bound by p53(R273H) in U251 human glioblastoma cells using chromatin immunoprecipitation (ChIP). Mutp53 binding to isolated DNA fragments confirmed the specificity of the ChIP. The mutp53 bound DNA sequences are rich in repetitive DNA elements, which are dispersed over non-coding DNA regions. Stable down-regulation of mutp53 expression strongly suggested that mutp53 binding to genomic DNA is functional. We identified the PPARGC1A and FRMD5 genes as p53(R273H) targets regulated by binding to intronic and intra-genic sequences. We propose a model that attributes the oncogenic functions of mutp53 to its ability to interact with intronic and intergenic non-B DNA sequences and modulate gene transcription via re-organization of chromatin.

Replication of eukaryotic chromosomes initiates at multiple sites called replication origins. Replication origins are best understood in the budding yeast Saccharomyces cerevisiae, where several complementary studies have mapped their locations genome-wide. We have collated these datasets, taking account of the resolution of each study, to generate a single list of distinct origin sites. OriDB provides a web-based catalogue of these confirmed and predicted S.cerevisiae DNA replication origin sites. Each proposed or confirmed origin site appears as a record in OriDB, with each record comprising seven pages. These pages provide, in text and graphical formats, the following information: genomic location and chromosome context of the origin site; time of origin replication; DNA sequence of proposed or experimentally confirmed origin elements; free energy required to open the DNA duplex (stress-induced DNA duplex destabilization or SIDD); and phylogenetic conservation of sequence elements. In addition, OriDB encourages community submission of additional information for each origin site through a User Notes facility. Origin sites are linked to several external resources, including the Saccharomyces Genome Database (SGD) and relevant publications at PubMed. Finally, a Chromosome Viewer utility allows users to interactively generate graphical representations of DNA replication data genome-wide. OriDB is available at www.oridb.org.

The transcription of ribosomal RNA (rRNA) is critical to life. Despite its importance, ribosomal DNA (rDNA) is not included in current genome assemblies and, consequently, genomic analyses to date have excluded rDNA. Here, we show that short sequence reads can be aligned to a genome assembly containing a single rDNA repeat. Integrated analysis of ChIP-seq, DNase-seq, MNase-seq and RNA-seq data reveals several novel findings. First, the coding region of active rDNA is contained within nucleosome-depleted open chromatin that is highly transcriptionally active. Second, histone modifications are located not only at the rDNA promoter but also at novel sites within the intergenic spacer. Third, the distributions of active modifications are more similar within and between different cell types than repressive modifications. Fourth, UBF, a positive regulator of rRNA transcription, binds to sites throughout the genome. Lastly, the insulator binding protein CTCF associates with the spacer promoter of rDNA, suggesting that transcriptional insulation plays a role in regulating the transcription of rRNA. Taken together, these analyses confirm and expand the results of previous ChIP studies of rDNA and provide novel avenues for exploration of chromatin-mediated regulation of rDNA.

A common practice in computational genomic analysis is to use a set of 'background' sequences as negative controls for evaluating the false-positive rates of prediction tools, such as gene identification programs and algorithms for detection of cis-regulatory elements. Such 'background' sequences are generally taken from regions of the genome presumed to be intergenic, or generated synthetically by 'shuffling' real sequences. This last method can lead to underestimation of false-positive rates. We developed a new method for generating artificial sequences that are modeled after real intergenic sequences in terms of composition, complexity and interspersed repeat content. These artificial sequences can serve as an inexhaustible source of high-quality negative controls. We used artificial sequences to evaluate the false-positive rates of a set of programs for detecting interspersed repeats, ab initio prediction of coding genes, transcribed regions and non-coding genes. We found that RepeatMasker is more accurate than PClouds, Augustus has the lowest false-positive rate of the coding gene prediction programs tested, and Infernal has a low false-positive rate for non-coding gene detection. A web service, source code and the models for human and many other species are freely available at http://repeatmasker.org/garlic/.

The genome consists of non-B-DNA structures such as G-quadruplexes (G4) that are involved in the regulation of genome stability and transcription. Telomeric-repeat containing RNA (TERRA) is capable of folding into G-quadruplex and interacting with chromatin remodeler ATRX. Here we show that TERRA modulates ATRX occupancy on repetitive sequences and over genes, and maintains DNA G-quadruplex structures at TERRA target and non-target sites in mouse embryonic stem cells. TERRA prevents ATRX from binding to subtelomeric regions and represses H3K9me3 formation. G4 ChIP-seq reveals that G4 abundance decreases at accessible chromatin regions, particularly at transcription start sites (TSS) after TERRA depletion; such G4 reduction at TSS is associated with elevated ATRX occupancy and differentially expressed genes. Loss of ATRX alleviates the effect of gene repression caused by TERRA depletion. Immunostaining analyses demonstrate that knockdown of TERRA diminishes DNA G4 signals, whereas silencing ATRX elevates G4 formation. Our results uncover an epigenetic regulation by TERRA that sequesters ATRX and preserves DNA G4 structures.

Genomic DNA (gDNA) contamination is an inherent problem during RNA purification that can lead to non-specific amplification and aberrant results in reverse transcription quantitative PCR (RT-qPCR). Currently, there is no alternative to RT(-) controls to evaluate the impact of the gDNA background on RT-PCR data. We propose a novel method (ValidPrime) that is more accurate than traditional RT(-) controls to test qPCR assays with respect to their sensitivity toward gDNA. ValidPrime measures the gDNA contribution using an optimized gDNA-specific ValidPrime assay (VPA) and gDNA reference sample(s). The VPA, targeting a non-transcribed locus, is used to measure the gDNA contents in RT(+) samples and the gDNA reference is used to normalize for GOI-specific differences in gDNA sensitivity. We demonstrate that the RNA-derived component of the signal can be accurately estimated and deduced from the total signal. ValidPrime corrects with high precision for both exogenous (spiked) and endogenous gDNA, contributing ∼60% of the total signal, whereas substantially reducing the number of required qPCR control reactions. In conclusion, ValidPrime offers a cost-efficient alternative to RT(-) controls and accurately corrects for signals derived from gDNA in RT-qPCR.

In Paramecium, the regeneration of a functional somatic genome at each sexual event relies on the elimination of thousands of germline DNA sequences, known as Internal Eliminated Sequences (IESs), from the zygotic nuclear DNA. Here, we provide evidence that IESs' length and sub-terminal bases jointly modulate IES excision by affecting DNA conformation in P. tetraurelia. Our study reveals an excess of complementary base pairing between IESs' sub-terminal and contiguous sites, suggesting that IESs may form DNA loops prior to cleavage. The degree of complementary base pairing between IESs' sub-terminal sites (termed Cin-score) is positively associated with IES length and is shaped by natural selection. Moreover, it escalates abruptly when IES length exceeds 45 nucleotides (nt), indicating that only sufficiently large IESs may form loops. Finally, we find that IESs smaller than 46 nt are favored targets of the cellular surveillance systems, presumably because of their relatively inefficient excision. Our findings extend the repertoire of cis-acting determinants for IES recognition/excision and provide unprecedented insights into the distinct selective pressures that operate on IESs and somatic DNA regions. This information potentially moves current models of IES evolution and of mechanisms of IES recognition/excision forward.

TET1 is a 5-methylcytosine dioxygenase and its DNA demethylating activity has been implicated in pluripotency and reprogramming. However, the precise role of TET1 in DNA methylation regulation outside of developmental reprogramming is still unclear. Here, we show that overexpression of the TET1 catalytic domain but not full length TET1 (TET1-FL) induces massive global DNA demethylation in differentiated cells. Genome-wide mapping reveals that 5-hydroxymethylcytosine production by TET1-FL is inhibited as DNA methylation increases, which can be explained by the preferential binding of TET1-FL to unmethylated CpG islands (CGIs) through its CXXC domain. TET1-FL specifically accumulates 5-hydroxymethylcytosine at the edges of hypomethylated CGIs, while knockdown of endogenous TET1 induces methylation spreading from methylated edges into hypomethylated CGIs. We also found that gene expression changes after TET1-FL overexpression are relatively small and independent of its dioxygenase function. Thus, our results identify TET1 as a maintenance DNA demethylase that does not purposely decrease methylation levels, but specifically prevents aberrant methylation spreading into CGIs in differentiated cells.

Pathogenic Vibrio species account for 3-5 million annual life-threatening human infections. Virulence is driven by bacterial hemolysin and toxin gene expression often positively regulated by the winged helix-turn-helix (wHTH) HlyU transcriptional regulator family and silenced by histone-like nucleoid structural protein (H-NS). In the case of Vibrio parahaemolyticus, HlyU is required for virulence gene expression associated with type 3 Secretion System-1 (T3SS1) although its mechanism of action is not understood. Here, we provide evidence for DNA cruciform attenuation mediated by HlyU binding to support concomitant virulence gene expression. Genetic and biochemical experiments revealed that upon HlyU mediated DNA cruciform attenuation, an intergenic cryptic promoter became accessible allowing for exsA mRNA expression and initiation of an ExsA autoactivation feedback loop at a separate ExsA-dependent promoter. Using a heterologous E. coli expression system, we reconstituted the dual promoter elements which revealed that HlyU binding and DNA cruciform attenuation were strictly required to initiate the ExsA autoactivation loop. The data indicate that HlyU acts to attenuate a transcriptional repressive DNA cruciform to support T3SS1 virulence gene expression and reveals a non-canonical extricating gene regulation mechanism in pathogenic Vibrio species.

Recent studies suggest that transcription takes place at DNA double-strand breaks (DSBs), that transcripts at DSBs are processed by Drosha and Dicer into damage-induced small RNAs (diRNAs), and that diRNAs are required for DNA repair. However, diRNAs have been mostly detected in reporter constructs or repetitive sequences, and their existence at endogenous loci has been questioned by recent reports. Using the homing endonuclease I-PpoI, we have investigated diRNA production in genetically unperturbed human and mouse cells. I-PpoI is an ideal tool to clarify the requirements for diRNA production because it induces DSBs in different types of loci: the repetitive 28S locus, unique genes and intergenic loci. We show by extensive sequencing that the rDNA locus produces substantial levels of diRNAs, whereas unique genic and intergenic loci do not. Further characterization of diRNAs emerging from the 28S locus reveals the existence of two diRNA subtypes. Surprisingly, Drosha and its partner DGCR8 are dispensable for diRNA production and only one diRNAs subtype depends on Dicer processing. Furthermore, we provide evidence that diRNAs are incorporated into Argonaute. Our findings provide direct evidence for diRNA production at endogenous loci in mammalian cells and give insights into RNA processing at DSBs.

The autonomous transcription of integrated retroviruses strongly depends on genetic and epigenetic effects of the chromatin at the site of integration. These effects are mostly suppressive and proviral activity can be finally silenced by mechanisms, such as DNA methylation and histone modifications. To address the role of the integration site at the whole-genome-scale, we performed clonal analysis of provirus silencing with an avian leucosis/sarcoma virus-based reporter vector and correlated the transcriptional silencing with the epigenomic landscape of respective integrations. We demonstrate efficient provirus silencing in human HCT116 cell line, which is strongly but not absolutely dependent on the de novo DNA methyltransferase activity, particularly of Dnmt3b. Proviruses integrated close to the transcription start sites of active genes into the regions enriched in H3K4 trimethylation display long-term stability of expression and are resistant to the transcriptional silencing after over-expression of Dnmt3a or Dnmt3b. In contrast, proviruses in the intergenic regions tend to spontaneous transcriptional silencing even in Dnmt3a(-/-) Dnmt3b(-/-) cells. The silencing of proviruses within genes is accompanied with DNA methylation of long terminal repeats, whereas silencing in intergenic regions is DNA methylation-independent. These findings indicate that the epigenomic features of integration sites are crucial for their permissivity to the proviral expression.

Post-translational modifications of core histones affect various cellular processes, primarily through transcription. However, their relationship with the termination of transcription has remained largely unknown. In this study, we show that DNA damage-activated AKT phosphorylates threonine 45 of core histone H3 (H3-T45). By genome-wide chromatin immunoprecipitation sequencing (ChIP-seq) analysis, H3-T45 phosphorylation was distributed throughout DNA damage-responsive gene loci, particularly immediately after the transcription termination site. H3-T45 phosphorylation pattern showed close-resemblance to that of RNA polymerase II C-terminal domain (CTD) serine 2 phosphorylation, which establishes the transcription termination signal. AKT1 was more effective than AKT2 in phosphorylating H3-T45. Blocking H3-T45 phosphorylation by inhibiting AKT or through amino acid substitution limited RNA decay downstream of mRNA cleavage sites and decreased RNA polymerase II release from chromatin. Our findings suggest that AKT-mediated phosphorylation of H3-T45 regulates the processing of the 3' end of DNA damage-activated genes to facilitate transcriptional termination.