GC skew is a conserved property of unmethylated CpG island promoters across vertebrates

RNA biogenesis and RNA metabolism factors as R-loop suppressors: a hidden role in genome integrity

Genome integrity relies on the accuracy of DNA metabolism, but as appreciated for more than four decades, transcription enhances mutation and recombination frequencies. More recent research provided evidence for a previously unforeseen link between RNA and DNA metabolism, which is often related to the accumulation of DNA-RNA hybrids and R-loops. In addition to physiological roles, R-loops interfere with DNA replication and repair, providing a molecular scenario for the origin of genome instability. Here, we review current knowledge on the multiple RNA factors that prevent or resolve R-loops and consequent transcription-replication conflicts and thus act as modulators of genome dynamics.

p53-dependent R-loop formation and HPV pathogenesis

R-loops are trimeric RNA: DNA hybrids that are important physiological regulators of transcription; however, their aberrant formation or turnover leads to genomic instability and DNA breaks. High-risk human papillomaviruses (HPV) are the causative agents of genital as well as oropharyngeal cancers and exhibit enhanced amounts of DNA breaks. The levels of R-loops were found to be increased up to 50-fold in cells that maintain high-risk HPV genomes and were readily detected in squamous cell cervical carcinomas in vivo but not in normal cells. The high levels of R-loops in HPV-positive cells were present on both viral and cellular sites together with RNase H1, an enzyme that controls their resolution. Depletion of RNase H1 in HPV-positive cells further increased R-loop levels, resulting in impaired viral transcription and replication along with reduced expression of the DNA repair genes such as FANCD2 and ATR, both of which are necessary for viral functions. Overexpression of RNase H1 decreased total R-loop levels, resulting in a reduction of DNA breaks by over 50%. Furthermore, increased RNase H1 expression blocked viral transcription and replication while enhancing the expression of factors in the innate immune regulatory pathway. This suggests that maintaining elevated R-loop levels is important for the HPV life cycle. The E6 viral oncoprotein was found to be responsible for inducing high levels of R-loops by inhibiting p53's transcriptional activity. Our studies indicate that high R-loop levels are critical for HPV pathogenesis and that this depends on suppressing the p53 pathway.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

Integrated Analyses of DNA Methylation and Gene Expression of Rainbow Trout Muscle under Variable Ploidy and Muscle Atrophy Conditions

Rainbow trout, Oncorhynchus mykiss, is an important cool, freshwater aquaculture species used as a model for biological research. However, its genome reference has not been annotated for epigenetic markers affecting various biological processes, including muscle growth/atrophy. Increased energetic demands during gonadogenesis/reproduction provoke muscle atrophy in rainbow trout. We described DNA methylation and its associated gene expression in atrophying muscle by comparing gravid, diploid females to sterile, triploid females. Methyl Mini-seq and RNA-Seq were simultaneously used to characterize genome-wide DNA methylation and its association with gene expression in rainbow trout muscle. Genome-wide enrichment in the number of CpGs, accompanied by depleted methylation levels, was noticed around the gene transcription start site (TSS). Hypermethylation of CpG sites within ±1 kb on both sides of TSS (promoter and gene body) was weakly/moderately associated with reduced gene expression. Conversely, hypermethylation of the CpG sites in downstream regions of the gene body +2 to +10 kb was weakly associated with increased gene expression. Unlike mammalian genomes, rainbow trout gene promotors are poor in CpG islands, at <1% compared to 60%. No signs of genome-wide, differentially methylated (DM) CpGs were observed due to the polyploidy effect; only 1206 CpGs (0.03%) were differentially methylated, and these were primarily associated with muscle atrophy. Twenty-eight genes exhibited differential gene expression consistent with methylation levels of 31 DM CpGs. These 31 DM CpGs represent potential epigenetic markers of muscle atrophy in rainbow trout. The DM CpG-harboring genes are involved in apoptosis, epigenetic regulation, autophagy, collagen metabolism, cell membrane functions, and Homeobox proteins. Our study also identified genes explaining higher water content and modulated glycolysis previously shown as characteristic biochemical signs of rainbow trout muscle atrophy associated with sexual maturation. This study characterized DNA methylation in the rainbow trout genome and its correlation with gene expression. This work also identified novel epigenetic markers associated with muscle atrophy in fish/lower vertebrates.

A CpG island promoter drives the CXXC5 gene expression

CXXC5 is a member of the zinc-finger CXXC family that binds to unmethylated CpG dinucleotides. CXXC5 modulates gene expressions resulting in diverse cellular events mediated by distinct signaling pathways. However, the mechanism responsible for CXXC5 expression remains largely unknown. We found here that of the 14 annotated CXXC5 transcripts with distinct 5' untranslated regions encoding the same protein, transcript variant 2 with the highest expression level among variants represents the main transcript in cell models. The DNA segment in and at the immediate 5'-sequences of the first exon of variant 2 contains a core promoter within which multiple transcription start sites are present. Residing in a region with high G-C nucleotide content and CpG repeats, the core promoter is unmethylated, deficient in nucleosomes, and associated with active RNA polymerase-II. These findings suggest that a CpG island promoter drives CXXC5 expression. Promoter pull-down revealed the association of various transcription factors (TFs) and transcription co-regulatory proteins, as well as proteins involved in histone/chromatin, DNA, and RNA processing with the core promoter. Of the TFs, we verified that ELF1 and MAZ contribute to CXXC5 expression. Moreover, the first exon of variant 2 may contain a G-quadruplex forming region that could modulate CXXC5 expression.

A functional LSD1 coregulator screen reveals a novel transcriptional regulatory cascade connecting R-loop homeostasis with epigenetic regulation

The lysine specific demethylase 1 (LSD1) plays a pivotal role in cellular differentiation by regulating the expression of key developmental genes in concert with different coregulatory proteins. This process is impaired in different cancer types and incompletely understood. To comprehensively identify functional coregulators of LSD1, we established a novel tractable fluorescent reporter system to monitor LSD1 activity in living cells. Combining this reporter system with a state-of-the-art multiplexed RNAi screen, we identify the DEAD-box helicase 19A (DDX19A) as a novel coregulator and demonstrate that suppression of Ddx19a results in an increase of R-loops and reduced LSD1-mediated gene silencing. We further show that DDX19A binds to tri-methylated lysine 27 of histone 3 (H3K27me3) and it regulates gene expression through the removal of transcription promoting R-loops. Our results uncover a novel transcriptional regulatory cascade where the downregulation of genes is dependent on the LSD1 mediated demethylation of histone H3 lysine 4 (H3K4). This allows the polycomb repressive complex 2 (PRC2) to methylate H3K27, which serves as a binding site for DDX19A. Finally, the binding of DDX19A leads to the efficient removal of R-loops at active promoters, which further de-represses LSD1 and PRC2, establishing a positive feedback loop leading to a robust repression of the target gene.

Best practices for the visualization, mapping, and manipulation of R-loops

R-loops represent an abundant class of large non-B DNA structures in genomes. Even though they form transiently and at modest frequencies, interfering with R-loop formation or dissolution has significant impacts on genome stability. Addressing the mechanism(s) of R-loop-mediated genome destabilization requires a precise characterization of their distribution in genomes. A number of independent methods have been developed to visualize and map R-loops, but their results are at times discordant, leading to confusion. Here, we review the main existing methodologies for R-loop mapping and assess their limitations as well as the robustness of existing datasets. We offer a set of best practices to improve the reproducibility of maps, hoping that such guidelines could be useful for authors and referees alike. Finally, we propose a possible resolution for the apparent contradictions in R-loop mapping outcomes between antibody-based and RNase H1-based mapping approaches.

G-quadruplex-R-loop interactions and the mechanism of anticancer G-quadruplex binders

Genomic DNA and cellular RNAs can form a variety of non-B secondary structures, including G-quadruplex (G4) and R-loops. G4s are constituted by stacked guanine tetrads held together by Hoogsteen hydrogen bonds and can form at key regulatory sites of eukaryote genomes and transcripts, including gene promoters, untranslated exon regions and telomeres. R-loops are 3-stranded structures wherein the two strands of a DNA duplex are melted and one of them is annealed to an RNA. Specific G4 binders are intensively investigated to discover new effective anticancer drugs based on a common rationale, i.e.: the selective inhibition of oncogene expression or specific impairment of telomere maintenance. However, despite the high number of known G4 binders, such a selective molecular activity has not been fully established and several published data point to a different mode of action. We will review published data that address the close structural interplay between G4s and R-loops in vitro and in vivo, and how these interactions can have functional consequences in relation to G4 binder activity. We propose that R-loops can play a previously-underestimated role in G4 binder action, in relation to DNA damage induction, telomere maintenance, genome and epigenome instability and alterations of gene expression programs.

Emerging roles for R-loop structures in the management of topological stress

R-loop structures are a prevalent class of alternative non-B DNA structures that form during transcription upon invasion of the DNA template by the nascent RNA. R-loops form universally in the genomes of organisms ranging from bacteriophages, bacteria, and yeasts to plants and animals, including mammals. A growing body of work has linked these structures to both physiological and pathological processes, in particular to genome instability. The rising interest in R-loops is placing new emphasis on understanding the fundamental physicochemical forces driving their formation and stability. Pioneering work in Escherichia coli revealed that DNA topology, in particular negative DNA superhelicity, plays a key role in driving R-loops. A clear role for DNA sequence was later uncovered. Here, we review and synthesize available evidence on the roles of DNA sequence and DNA topology in controlling R-loop formation and stability. Factoring in recent developments in R-loop modeling and single-molecule profiling, we propose a coherent model accounting for the interplay between DNA sequence and DNA topology in driving R-loop structure formation. This model reveals R-loops in a new light as powerful and reversible topological stress relievers, an insight that significantly expands the repertoire of R-loops' potential biological roles under both normal and aberrant conditions.

Ultra-deep Coverage Single-molecule R-loop Footprinting Reveals Principles of R-loop Formation

R-loops are a prevalent class of non-B DNA structures that have been associated with both positive and negative cellular outcomes. DNA:RNA immunoprecipitation (DRIP) approaches based on the anti-DNA:RNA hybrid S9.6 antibody revealed that R-loops form dynamically over conserved genic hotspots. We have developed an orthogonal approach that queries R-loops via the presence of long stretches of single-stranded DNA on their looped-out strand. Nondenaturing sodium bisulfite treatment catalyzes the conversion of unpaired cytosines to uracils, creating permanent genetic tags for the position of an R-loop. Long-read, single-molecule PacBio sequencing allows the identification of R-loop 'footprints' at near nucleotide resolution in a strand-specific manner on long single DNA molecules and at ultra-deep coverage. Single-molecule R-loop footprinting coupled with PacBio sequencing (SMRF-seq) revealed a strong agreement between S9.6-based and bisulfite-based R-loop mapping and confirmed that R-loops form over genic hotspots, including gene bodies and terminal gene regions. Based on the largest single-molecule R-loop dataset to date, we show that individual R-loops form nonrandomly, defining discrete sets of overlapping molecular clusters that pileup through larger R-loop zones. R-loops most often map to intronic regions and their individual start and stop positions do not match with intron-exon boundaries, reinforcing the model that they form cotranscriptionally from unspliced transcripts. SMRF-seq further established that R-loop distribution patterns are not simply driven by intrinsic DNA sequence features but most likely also reflect DNA topological constraints. Overall, DRIP-based and SMRF-based approaches independently provide a complementary and congruent view of R-loop distribution, consolidating our understanding of the principles underlying R-loop formation.

The evolution of CpG density and lifespan in conserved primate and mammalian promoters

Gene promoters are evolutionarily conserved across holozoans and enriched in CpG sites, the target for DNA methylation. As animals age, the epigenetic pattern of DNA methylation degrades, with highly methylated CpG sites gradually becoming demethylated while CpG islands increase in methylation. Across vertebrates, aging is a trait that varies among species. We used this variation to determine whether promoter CpG density correlates with species’ maximum lifespan. Human promoter sequences were used to identify conserved regions in 131 mammals and a subset of 28 primate genomes. We identified approximately 1000 gene promoters (5% of the total), that significantly correlated CpG density with lifespan. The correlations were performed via the phylogenetic least squares method to account for trait similarity by common descent using phylogenetic branch lengths. Gene set enrichment analysis revealed no significantly enriched pathways or processes, consistent with the hypothesis that aging is not under positive selection. However, within both mammals and primates, 95% of the promoters showed a positive correlation between increasing CpG density and species lifespan, and two thirds were shared between the primate subset and mammalian datasets. Thus, these genes may require greater buffering capacity against age-related dysregulation of DNA methylation in longer-lived species.

Non-canonical DNA/RNA structures during Transcription-Coupled Double-Strand Break Repair: Roadblocks or Bona fide repair intermediates?

Although long overlooked, it is now well understood that DNA does not systematically assemble into a canonical double helix, known as B-DNA, throughout the entire genome but can also accommodate other structures including DNA hairpins, G-quadruplexes and RNA:DNA hybrids. Notably, these non-canonical DNA structures form preferentially at transcriptionally active loci. Acting as replication roadblocks and being targeted by multiple machineries, these structures weaken the genome and render it prone to damage, including DNA double-strand breaks (DSB). In addition, secondary structures also further accumulate upon DSB formation. Here we discuss the potential functions of pre-existing or de novo formed nucleic acid structures, as bona fide repair intermediates or repair roadblocks, especially during Transcription-Coupled DNA Double-Strand Break repair (TC-DSBR), and provide an update on the specialized protein complexes displaying the ability to remove these structures to safeguard genome integrity.

Comparative DNA methylomic analyses reveal potential origins of novel epigenetic biomarkers of insulin resistance in monocytes from virally suppressed HIV-infected adults

BACKGROUND

Compared to healthy individuals, those with stably repressed HIV experience a higher risk of developing insulin resistance, a hallmark of pre-diabetes and a major determinant for cardiometabolic diseases. Although epigenetic processes, including in particular DNA methylation, appear to be dysregulated in individuals with insulin resistance, little is known about where these occur in the genomes of immune cells and the origins of these alterations in HIV-infected individuals. Here, we examined the genome-wide DNA methylation states of monocytes in HIV-infected individuals (n = 37) with varying levels of insulin sensitivity measured by the homeostatic model assessment of insulin resistance (HOMA-IR).

RESULTS

By profiling DNA methylation at single-nucleotide resolution using the Illumina Infinium HumanMethylation450 BeadChip in monocytes from insulin-resistant (IR; HOMA-IR ≥ 2.0; n = 14) and insulin-sensitive (IS; HOMA-IR < 2.0; n = 23) individuals, we identified 123 CpGs with significantly different DNA methylation levels. These CpGs were enriched at genes involved in pathways relating to glucose metabolism, immune activation, and insulin-relevant signaling, with the majority (86.2%) being hypomethylated in IR relative to IS individuals. Using a stepwise multiple logistic regression analysis, we observed 4 CpGs (cg27655935, cg02000426, cg10184328, and cg23085143) whose methylation levels independently predicted the insulin-resistant state at a higher confidence than that of clinical risk factors typically associated with insulin resistance (i.e., fasting glucose, 120-min oral glucose tolerance test, Framingham Risk Score, and Total to HDL cholesterol ratio). Interestingly, 79 of the 123 CpGs (64%) exhibited remarkably similar levels of methylation as that of hematopoietic stem cells (HSC) in monocytes from IR individuals, implicating epigenetic defects in myeloid differentiation as a possible origin for the methylation landscape underlying the insulin resistance phenotype. In support of this, gene ontology analysis of these 79 CpGs revealed overrepresentation of these CpGs at genes relevant to HSC function, including involvement in stem cell pluripotency, differentiation, and Wnt signaling pathways.

CONCLUSION

Altogether, our data suggests a possible role for DNA methylation in regulating monocyte activity that may associate with the insulin-resistant phenotype. The methylomic landscape of insulin resistance in monocytes could originate from epigenetic dysregulation during HSC differentiation through the myeloid lineage. Understanding the factors involved with changes in the myeloid trajectory may provide further insight into the development of insulin resistance. Furthermore, regulation of specific genes that were implicated in our analysis reveal possible targets for modulating immune activity to ameliorate insulin resistance.

High-resolution, strand-specific R-loop mapping via S9.6-based DNA-RNA immunoprecipitation and high-throughput sequencing

R-loops are prevalent three-stranded non-B DNA structures composed of an RNA-DNA hybrid and a single strand of DNA. R-loops are implicated in various basic nuclear processes, such as class-switch recombination, transcription termination and chromatin patterning. Perturbations in R-loop metabolism have been linked to genomic instability and have been implicated in human disorders, including cancer. As a consequence, the accurate mapping of these structures has been of increasing interest in recent years. Here, we describe two related immunoprecipitation-based methods for mapping R-loop structures: basic DRIP-seq (DNA-RNA immunoprecipitation followed by high-throughput DNA sequencing), an easy, robust, but resolution-limited technique; and DRIPc-seq (DNA-RNA immunoprecipitation followed by cDNA conversion coupled to high-throughput sequencing), a high-resolution and strand-specific iteration of the method that permits accurate R-loop mapping genome wide. Briefly, after gentle DNA extraction and restriction digestion with a cocktail of enzymes, R-loop structures are immunoprecipitated with the anti-RNA-DNA hybrid S9.6 antibody. Compared with DRIP-seq, in which the immunoprecipitated DNA is directly sequenced, DRIPc-seq permits the recovery of the RNA moiety of R-loops, and these RNA strands are subjected to strand-specific RNA sequencing (RNA-seq) analysis. DRIPc-seq can be performed in 5 d and can be applied to any cell type, provided sufficient starting material can be collected. Accurately mapping R-loop distribution in various cell lines and under varied conditions is essential to understanding the formation, roles and dynamic resolution of these important structures.

Interplay between DNA sequence and negative superhelicity drives R-loop structures

Three-stranded R-loop structures form during transcription when the nascent RNA transcript rehybridizes to the template DNA strand. This creates an RNA:DNA hybrid and forces the nontemplate DNA strand into a single-stranded, looped-out state. R-loops form universally over conserved hotspot regions. To date, the physicochemical bases underlying R-loop formation remain unclear. Using a “first-principle” mathematical approach backed by experimental validation, we elucidated the relative contributions of DNA sequence and DNA topology to R-loop formation. Our work provides a quantitative assessment of the energies underlying R-loop formation and of their interplay. It further reveals these structures as important regulators of the DNA topological state.

R-loops are abundant three-stranded nucleic-acid structures that form in cis during transcription. Experimental evidence suggests that R-loop formation is affected by DNA sequence and topology. However, the exact manner by which these factors interact to determine R-loop susceptibility is unclear. To investigate this, we developed a statistical mechanical equilibrium model of R-loop formation in superhelical DNA. In this model, the energy involved in forming an R-loop includes four terms—junctional and base-pairing energies and energies associated with superhelicity and with the torsional winding of the displaced DNA single strand around the RNA:DNA hybrid. This model shows that the significant energy barrier imposed by the formation of junctions can be overcome in two ways. First, base-pairing energy can favor RNA:DNA over DNA:DNA duplexes in favorable sequences. Second, R-loops, by absorbing negative superhelicity, partially or fully relax the rest of the DNA domain, thereby returning it to a lower energy state. In vitro transcription assays confirmed that R-loops cause plasmid relaxation and that negative superhelicity is required for R-loops to form, even in a favorable region. Single-molecule R-loop footprinting following in vitro transcription showed a strong agreement between theoretical predictions and experimental mapping of stable R-loop positions and further revealed the impact of DNA topology on the R-loop distribution landscape. Our results clarify the interplay between base sequence and DNA superhelicity in controlling R-loop stability. They also reveal R-loops as powerful and reversible topology sinks that cells may use to nonenzymatically relieve superhelical stress during transcription.

On the length, weight and GC content of the human genome

Objective

Basic parameters commonly used to describe genomes including length, weight and relative guanine-cytosine (GC) content are widely cited in absence of a primary source. By using updated data and original software we determined these values to the best of our knowledge as standard reference for the whole human nuclear genome, for each chromosome and for mitochondrial DNA. We also devised a method to calculate the relative GC content in the whole messenger RNA sequence set and in transcriptomes by multiplying the GC content of each gene by its mean expression level.

Results

The male nuclear diploid genome extends for 6.27 Gigabase pairs (Gbp), is 205.00 cm (cm) long and weighs 6.41 picograms (pg). Female values are 6.37 Gbp, 208.23 cm, 6.51 pg. The individual variability and the implication for the DNA informational density in terms of bits/volume were discussed. The genomic GC content is 40.9%. Following analysis in different transcriptomes and species, we showed that the greatest deviation was observed in the pathological condition analysed (trisomy 21 leukaemic cells) and in Caenorhabditis elegans. Our results may represent a solid basis for further investigation on human structural and functional genomics while also providing a framework for other genome comparative analysis.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4137-z) contains supplementary material, which is available to authorized users.

Walking along chromosomes with super-resolution imaging, contact maps, and integrative modeling

Chromosome organization is crucial for genome function. Here, we present a method for visualizing chromosomal DNA at super-resolution and then integrating Hi-C data to produce three-dimensional models of chromosome organization. Using the super-resolution microscopy methods of OligoSTORM and OligoDNA-PAINT, we trace 8 megabases of human chromosome 19, visualizing structures ranging in size from a few kilobases to over a megabase. Focusing on chromosomal regions that contribute to compartments, we discover distinct structures that, in spite of considerable variability, can predict whether such regions correspond to active (A-type) or inactive (B-type) compartments. Imaging through the depths of entire nuclei, we capture pairs of homologous regions in diploid cells, obtaining evidence that maternal and paternal homologous regions can be differentially organized. Finally, using restraint-based modeling to integrate imaging and Hi-C data, we implement a method-integrative modeling of genomic regions (IMGR)-to increase the genomic resolution of our traces to 10 kb.

DNA Topoisomerase I differentially modulates R-loops across the human genome

BACKGROUND

Co-transcriptional R-loops are abundant non-B DNA structures in mammalian genomes. DNA Topoisomerase I (Top1) is often thought to regulate R-loop formation owing to its ability to resolve both positive and negative supercoils. How Top1 regulates R-loop structures at a global level is unknown.

RESULTS

Here, we perform high-resolution strand-specific R-loop mapping in human cells depleted for Top1 and find that Top1 depletion results in both R-loop gains and losses at thousands of transcribed loci, delineating two distinct gene classes. R-loop gains are characteristic for long, highly transcribed, genes located in gene-poor regions anchored to Lamin B1 domains and in proximity to H3K9me3-marked heterochromatic patches. R-loop losses, by contrast, occur in gene-rich regions overlapping H3K27me3-marked active replication initiation regions. Interestingly, Top1 depletion coincides with a block of the cell cycle in G0/G1 phase and a trend towards replication delay.

CONCLUSIONS

Our findings reveal new properties of Top1 in regulating R-loop homeostasis in a context-dependent manner and suggest a potential role for Top1 in modulating the replication process via R-loop formation.

Dynamic bimodal changes in CpG and non-CpG methylation genome-wide upon CGGBP1 loss-of-function

OBJECTIVES

Although CpG methylation is well studied, mechanisms of non-CpG methylation in mammals remains elusive. Studying proteins with non-CpG cytosine methylation-sensitive DNA-binding, such as human CGGBP1, can unveil cytosine methylation regulatory mechanisms. Here we have resequenced a published genome-wide bisulfite sequencing library and analyzed it at base level resolution. CpG, CHG and CHH (where H is any nucleotide other than G) methylation states in non-targeting or CGGBP1-targeting shmiR lentivirus-transduced cells have been analyzed to identify how CGGBP1 regulates CpG and non-CpG methylation.

RESULTS

We report that CGGBP1 acts as a dynamic bimodal balancer of methylation. Both gain and loss of methylation observed upon CGGBP1 depletion were spatially overlapping at annotated functional regions and not identifiable with any sequence motifs but clearly associated with GC-skew. CGGBP1 depletion caused clustered methylation changes in cis, upstream of R-loop forming promoters. This was complemented by clustered occurrences of methylation changes in proximity of transcription start sites of known cytosine methylation regulatory genes, altered expression of which can regulate cytosine methylation in trans. Despite low coverage, our data provide reliable estimates of the spectrum of methylation changes regulated by CGGBP1 in all cytosine contexts genome-wide through a combination of cis and trans-acting mechanisms.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

Unlinking the methylome pattern from nucleotide sequence, revealed by large-scale in vivo genome engineering and methylome editing in medaka fish

The heavily methylated vertebrate genomes are punctuated by stretches of poorly methylated DNA sequences that usually mark gene regulatory regions. It is known that the methylation state of these regions confers transcriptional control over their associated genes. Given its governance on the transcriptome, cellular functions and identity, genome-wide DNA methylation pattern is tightly regulated and evidently predefined. However, how is the methylation pattern determined in vivo remains enigmatic. Based on in silico and in vitro evidence, recent studies proposed that the regional hypomethylated state is primarily determined by local DNA sequence, e.g., high CpG density and presence of specific transcription factor binding sites. Nonetheless, the dependency of DNA methylation on nucleotide sequence has not been carefully validated in vertebrates in vivo. Herein, with the use of medaka (Oryzias latipes) as a model, the sequence dependency of DNA methylation was intensively tested in vivo. Our statistical modeling confirmed the strong statistical association between nucleotide sequence pattern and methylation state in the medaka genome. However, by manipulating the methylation state of a number of genomic sequences and reintegrating them into medaka embryos, we demonstrated that artificially conferred DNA methylation states were predominantly and robustly maintained in vivo, regardless of their sequences and endogenous states. This feature was also observed in the medaka transgene that had passed across generations. Thus, despite the observed statistical association, nucleotide sequence was unable to autonomously determine its own methylation state in medaka in vivo. Our results apparently argue against the notion of the governance on the DNA methylation by nucleotide sequence, but instead suggest the involvement of other epigenetic factors in defining and maintaining the DNA methylation landscape. Further investigation in other vertebrate models in vivo will be needed for the generalization of our observations made in medaka.

The genomes of vertebrate animals are naturally and extensively modified by methylation. The DNA methylation is essential to normal functions of cells, hence the whole animal, since it governs gene expression. Defects in the establishment and maintenance of proper methylation pattern are commonly associated with various developmental abnormalities and diseases. How exactly is the normal pattern defined in vertebrate animals is not fully understood, but recent researches with computational analyses and cultured cells suggested that DNA sequence is a primary determinant of the methylation pattern. This study encompasses the first experiments that rigorously test this notion in whole animal (medaka fish). In statistical sense, we observed the very strong correlation between DNA sequence and methylation state. However, by introducing unmethylated and artificially methylated native genomic DNA sequences into the genome, we demonstrated that the artificially conferred methylation states were robustly maintained in the animal, independent of the sequence and native state. Our results thus demonstrate that genome-wide DNA methylation pattern is not autonomously determined by the DNA sequence, which underpins the vital role of DNA methylation pattern as a core epigenetic element.

Ssb1 and Ssb2 cooperate to regulate mouse hematopoietic stem and progenitor cells by resolving replicative stress

Hematopoietic stem and progenitor cells (HSPCs) are vulnerable to endogenous damage and defects in DNA repair can limit their function. The 2 single-stranded DNA (ssDNA) binding proteins SSB1 and SSB2 are crucial regulators of the DNA damage response; however, their overlapping roles during normal physiology are incompletely understood. We generated mice in which both Ssb1 and Ssb2 were constitutively or conditionally deleted. Constitutive Ssb1/Ssb2 double knockout (DKO) caused early embryonic lethality, whereas conditional Ssb1/Ssb2 double knockout (cDKO) in adult mice resulted in acute lethality due to bone marrow failure and intestinal atrophy featuring stem and progenitor cell depletion, a phenotype unexpected from the previously reported single knockout models of Ssb1 or Ssb2 Mechanistically, cDKO HSPCs showed altered replication fork dynamics, massive accumulation of DNA damage, genome-wide double-strand breaks enriched at Ssb-binding regions and CpG islands, together with the accumulation of R-loops and cytosolic ssDNA. Transcriptional profiling of cDKO HSPCs revealed the activation of p53 and interferon (IFN) pathways, which enforced cell cycling in quiescent HSPCs, resulting in their apoptotic death. The rapid cell death phenotype was reproducible in in vitro cultured cDKO-hematopoietic stem cells, which were significantly rescued by nucleotide supplementation or after depletion of p53. Collectively, Ssb1 and Ssb2 control crucial aspects of HSPC function, including proliferation and survival in vivo by resolving replicative stress to maintain genomic stability.

Telomeres in ICF syndrome cells are vulnerable to DNA damage due to elevated DNA:RNA hybrids

DNA:RNA hybrids, nucleic acid structures with diverse physiological functions, can disrupt genome integrity when dysregulated. Human telomeres were shown to form hybrids with the lncRNA TERRA, yet the formation and distribution of these hybrids among telomeres, their regulation and their cellular effects remain elusive. Here we predict and confirm in several human cell types that DNA:RNA hybrids form at many subtelomeric and telomeric regions. We demonstrate that ICF syndrome cells, which exhibit short telomeres and elevated TERRA levels, are enriched for hybrids at telomeric regions throughout the cell cycle. Telomeric hybrids are associated with high levels of DNA damage at chromosome ends in ICF cells, which are significantly reduced with overexpression of RNase H1. Our findings suggest that abnormally high TERRA levels in ICF syndrome lead to accumulation of telomeric hybrids that, in turn, can result in telomeric dysfunction.

ICF syndrome cells exhibit shortened telomeres and elevated levels of the noncoding RNA TERRA. Here the authors show this is associated with high levels of DNA damage, suggesting an increase in telomere dysfunction due to the formation of DNA: RNA hybrids

R-loopDB: a database for R-loop forming sequences (RLFS) and R-loops

R-loopDB (http://rloop.bii.a-star.edu.sg) was originally constructed as a collection of computationally predicted R-loop forming sequences (RLFSs) in the human genic regions. The renewed R-loopDB provides updates, improvements and new options, including access to recent experimental data. It includes genome-scale prediction of RLFSs for humans, six other animals and yeast. Using the extended quantitative model of RLFSs (QmRLFS), we significantly increased the number of RLFSs predicted in the human genes and identified RLFSs in other organism genomes. R-loopDB allows searching of RLFSs in the genes and in the 2 kb upstream and downstream flanking sequences of any gene. R-loopDB exploits the Ensembl gene annotation system, providing users with chromosome coordinates, sequences, gene and genomic data of the 1 565 795 RLFSs distributed in 121 056 genic or proximal gene regions of the covered organisms. It provides a comprehensive annotation of Ensembl RLFS-positive genes including 93 454 protein coding genes, 12 480 long non-coding RNA and 7 568 small non-coding RNA genes and 7 554 pseudogenes. Using new interface and genome viewers of R-loopDB, users can search the gene(s) in multiple species with keywords in a single query. R-loopDB provides tools to carry out comparative evolution and genome-scale analyses in R-loop biology.

Nascent Connections: R-Loops and Chromatin Patterning

RNA molecules, such as long noncoding RNAs (lncRNAs), have critical roles in regulating gene expression, chromosome architecture, and the modification states of chromatin. Recent developments suggest that RNA also influences gene expression and chromatin patterns through the interaction of nascent transcripts with their DNA template via the formation of co-transcriptional R-loop structures. R-loop formation over specific, conserved, hotspots occurs at thousands of genes in mammalian genomes and represents an important and dynamic feature of mammalian chromatin. Here, focusing primarily on mammalian systems, I describe the accumulating connections and possible mechanisms linking R-loop formation and chromatin patterning. The possible contribution of aberrant R-loops to pathological conditions is also discussed.

R-loop: an emerging regulator of chromatin dynamics

The dynamic structure of chromatin, which exists in two conformational states: heterochromatin and euchromatin, alters the accessibility of the DNA to regulatory factors during transcription, replication, recombination, and DNA damage repair. Chemical modifications of histones and DNA, as well as adenosine triphospahate-dependent nucleosome remodeling, have been the major focus of research on chromatin dynamics over the past two decades. However, recent studies using a DNA-RNA hybrid-specific antibody and next-generation sequencing approaches have revealed that the formation of R-loops, one of the most common non-canonical DNA structures, is an emerging regulator of chromatin states. This review focuses on recent insights into the interplay between R-loop formation and the epigenetic modifications of chromatin in normal and disease states.

Regulation of chromatin structure and cell fate by R-loops

Hybridization of RNA to its template DNA strand during transcription induces formation of R-loops-RNA:DNA hybrids with unpaired non-template DNA strands. Although unresolved R-loops can be detrimental, some R-loops contribute to regulation of chromatin structure. Consequently, R-loops help regulate gene expression and play important roles in numerous cellular processes.

Distribution of CpG Motifs in Upstream Gene Domains in a Reef Coral and Sea Anemone: Implications for Epigenetics in Cnidarians

Coral reefs are under assault from stressors including global warming, ocean acidification, and urbanization. Knowing how these factors impact the future fate of reefs requires delineating stress responses across ecological, organismal and cellular scales. Recent advances in coral reef biology have integrated molecular processes with ecological fitness and have identified putative suites of temperature acclimation genes in a Scleractinian coral Acropora hyacinthus. We wondered what unique characteristics of these genes determined their coordinate expression in response to temperature acclimation, and whether or not other corals and cnidarians would likewise possess these features. Here, we focus on cytosine methylation as an epigenetic DNA modification that is responsive to environmental stressors. We identify common conserved patterns of cytosine-guanosine dinucleotide (CpG) motif frequencies in upstream promoter domains of different functional gene groups in two cnidarian genomes: a coral (Acropora digitifera) and an anemone (Nematostella vectensis). Our analyses show that CpG motif frequencies are prominent in the promoter domains of functional genes associated with environmental adaptation, particularly those identified in A. hyacinthus. Densities of CpG sites in upstream promoter domains near the transcriptional start site (TSS) are 1.38x higher than genomic background levels upstream of -2000 bp from the TSS. The increase in CpG usage suggests selection to allow for DNA methylation events to occur more frequently within 1 kb of the TSS. In addition, observed shifts in CpG densities among functional groups of genes suggests a potential role for epigenetic DNA methylation within promoter domains to impact functional gene expression responses in A. digitifera and N. vectensis. Identifying promoter epigenetic sequence motifs among genes within specific functional groups establishes an approach to describe integrated cellular responses to environmental stress in reef corals and potential roles of epigenetics on survival and fitness in the face of global climate change.