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

GC skew is a measure of the strand asymmetry in the distribution of guanines and cytosines. GC skew favors R-loops, a type of three stranded nucleic acid structures that form upon annealing of an RNA strand to one strand of DNA, creating a persistent RNA:DNA hybrid. Previous studies show that GC skew is prevalent at thousands of human CpG island (CGI) promoters and transcription termination regions, which correspond to hotspots of R-loop formation. Here, we investigated the conservation of GC skew patterns in 60 sequenced chordates genomes. We report that GC skew is a conserved sequence characteristic of the CGI promoter class in vertebrates. Furthermore, we reveal that promoter GC skew peaks at the exon 1/ intron1 junction and that it is highly correlated with gene age and CGI promoter strength. Our data also show that GC skew is predictive of unmethylated CGI promoters in a range of vertebrate species and that it imparts significant DNA hypomethylation for promoters with intermediate CpG densities. Finally, we observed that terminal GC skew is conserved for a subset of vertebrate genes that tend to be located significantly closer to their downstream neighbors, consistent with a role for R-loop formation in transcription termination

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

Prevalent, Dynamic, and Conserved R-Loop Structures Associate with Specific Epigenomic Signatures in Mammals.

R-loops are three-stranded nucleic acid structures formed upon annealing of an RNA strand to one strand of duplex DNA. We profiled R-loops using a high-resolution, strand-specific methodology in human and mouse cell types. R-loops are prevalent, collectively occupying up to 5% of mammalian genomes. R-loop formation occurs over conserved genic hotspots such as promoter and terminator regions of poly(A)-dependent genes. In most cases, R-loops occur co-transcriptionally and undergo dynamic turnover. Detailed epigenomic profiling revealed that R-loops associate with specific chromatin signatures. At promoters, R-loops associate with a hyper-accessible state characteristic of unmethylated CpG island promoters. By contrast, terminal R-loops associate with an enhancer- and insulator-like state and define a broad class of transcription terminators. Together, this suggests that the retention of nascent RNA transcripts at their site of expression represents an abundant, dynamic, and programmed component of the mammalian chromatin that affects chromatin patterning and the control of gene expression

Sequencing the transcriptome of milk production: milk trumps mammary tissue

Background: Studies of normal human mammary gland development and function have mostly relied on cell culture, limited surgical specimens, and rodent models. Although RNA extracted from human milk has been used to assay the mammary transcriptome non-invasively, this assay has not been adequately validated in primates. Thus, the objectives of the current study were to assess the suitability of lactating rhesus macaques as a model for lactating humans and to determine whether RNA extracted from milk fractions is representative of RNA extracted from mammary tissue for the purpose of studying the transcriptome of milk-producing cells. Results: We confirmed that macaque milk contains cytoplasmic crescents and that ample high-quality RNA can be obtained for sequencing. Using RNA sequencing, RNA extracted from macaque milk fat and milk cell fractions more accurately represented RNA from mammary epithelial cells (cells that produce milk) than did RNA from whole mammary tissue. Mammary epithelium-specific transcripts were more abundant in macaque milk fat, whereas adipose or stroma-specific transcripts were more abundant in mammary tissue. Functional analyses confirmed the validity of milk as a source of RNA from milk-producing mammary epithelial cells. Conclusions: RNA extracted from the milk fat during lactation accurately portrayed the RNA profile of milk-producing mammary epithelial cells in a non-human primate. However, this sample type clearly requires protocols that minimize RNA degradation. Overall, we validated the use of RNA extracted from human and macaque milk and provided evidence to support the use of lactating macaques as a model for human lactation

Inactive DNMT3B splice variants modulate de novo DNA methylation.

Inactive DNA methyltransferase (DNMT) 3B splice isoforms are associated with changes in DNA methylation, yet the mechanisms by which they act remain largely unknown. Using biochemical and cell culture assays, we show here that the inactive DNMT3B3 and DNMT3B4 isoforms bind to and regulate the activity of catalytically competent DNMT3A or DNMT3B molecules. DNMT3B3 modestly stimulated the de novo methylation activity of DNMT3A and also counteracted the stimulatory effects of DNMT3L, therefore leading to subtle and contrasting effects on activity. DNMT3B4, by contrast, significantly inhibited de novo DNA methylation by active DNMT3 molecules, most likely due to its ability to reduce the DNA binding affinity of co-complexes, thereby sequestering them away from their substrate. Immunocytochemistry experiments revealed that in addition to their effects on the intrinsic catalytic function of active DNMT3 enzymes, DNMT3B3 and DNMT34 drive distinct types of chromatin compaction and patterns of histone 3 lysine 9 tri-methylation (H3K9me3) deposition. Our findings suggest that regulation of active DNMT3 members through the formation of co-complexes with inactive DNMT3 variants is a general mechanism by which DNMT3 variants function. This may account for some of the changes in DNA methylation patterns observed during development and disease