Estimating enrichment of repetitive elements from high-throughput sequence data

Foxg1 bimodally tunes L1-mRNA and -DNA dynamics in the developing murine neocortex

Foxg1 masters telencephalic development via a pleiotropic control over its progression. Expressed within the central nervous system (CNS), L1 retrotransposons are implicated in progression of its histogenesis and tuning of its genomic plasticity. Foxg1 represses gene transcription, and L1 elements share putative Foxg1-binding motifs, suggesting the former might limit telencephalic expression (and activity) of the latter. We tested such a prediction, in vivo as well as in engineered primary neural cultures, using loss- and gain-of-function approaches. We found that Foxg1-dependent, transcriptional L1 repression specifically occurs in neopallial neuronogenic progenitors and post-mitotic neurons, where it is supported by specific changes in the L1 epigenetic landscape. Unexpectedly, we discovered that Foxg1 physically interacts with L1-mRNA and positively regulates neonatal neopallium L1-DNA content, antagonizing the retrotranscription-suppressing activity exerted by Mov10 and Ddx39a helicases. To the best of our knowledge, Foxg1 represents the first CNS patterning gene acting as a bimodal retrotransposon modulator, limiting transcription of L1 elements and promoting their amplification, within a specific domain of the developing mouse brain.

RepEnTools: an automated repeat enrichment analysis package for ChIP-seq data reveals hUHRF1 Tandem-Tudor domain enrichment in young repeats

BACKGROUND

Repeat elements (REs) play important roles for cell function in health and disease. However, RE enrichment analysis in short-read high-throughput sequencing (HTS) data, such as ChIP-seq, is a challenging task.

RESULTS

Here, we present RepEnTools, a software package for genome-wide RE enrichment analysis of ChIP-seq and similar chromatin pulldown experiments. Our analysis package bundles together various software with carefully chosen and validated settings to provide a complete solution for RE analysis, starting from raw input files to tabular and graphical outputs. RepEnTools implementations are easily accessible even with minimal IT skills (Galaxy/UNIX). To demonstrate the performance of RepEnTools, we analysed chromatin pulldown data by the human UHRF1 TTD protein domain and discovered enrichment of TTD binding on young primate and hominid specific polymorphic repeats (SVA, L1PA1/L1HS) overlapping known enhancers and decorated with H3K4me1-K9me2/3 modifications. We corroborated these new bioinformatic findings with experimental data by qPCR assays using newly developed primate and hominid specific qPCR assays which complement similar research tools. Finally, we analysed mouse UHRF1 ChIP-seq data with RepEnTools and showed that the endogenous mUHRF1 protein colocalizes with H3K4me1-H3K9me3 on promoters of REs which were silenced by UHRF1. These new data suggest a functional role for UHRF1 in silencing of REs that is mediated by TTD binding to the H3K4me1-K9me3 double mark and conserved in two mammalian species.

CONCLUSIONS

RepEnTools improves the previously available programmes for RE enrichment analysis in chromatin pulldown studies by leveraging new tools, enhancing accessibility and adding some key functions. RepEnTools can analyse RE enrichment rapidly, efficiently, and accurately, providing the community with an up-to-date, reliable and accessible tool for this important type of analysis.

EZH2 inhibition stimulates repetitive element expression and viral mimicry in resting splenic B cells

Mammalian cells repress expression of repetitive genomic sequences by forming heterochromatin. However, the consequences of ectopic repeat expression remain unclear. Here we demonstrate that inhibitors of EZH2, the catalytic subunit of the Polycomb repressive complex 2 (PRC2), stimulate repeat misexpression and cell death in resting splenic B cells. B cells are uniquely sensitive to these agents because they exhibit high levels of histone H3 lysine 27 trimethylation (H3K27me3) and correspondingly low DNA methylation at repeat elements. We generated a pattern recognition receptor loss-of-function mouse model, called RIC, with mutations in Rigi (encoding for RIG-I), Ifih1 (MDA5), and Cgas. In both wildtype and RIC mutant B cells, EZH2 inhibition caused loss of H3K27me3 at repetitive elements and upregulated their expression. However, NF-κB-dependent expression of inflammatory chemokines and subsequent cell death was suppressed by the RIC mutations. We further show that inhibition of EZH2 in cancer cells requires the same pattern recognition receptors to activate an interferon response. Together, the results reveal chemokine expression induced by EZH2 inhibitors in B cells as a novel inflammatory response to genomic repeat expression. Given the overlap of genes induced by EZH2 inhibitors and Epstein-Barr virus infection, this response can be described as a form of viral mimicry.

Pediatric glioma histone H3.3 K27M/G34R mutations drive abnormalities in PML nuclear bodies

BACKGROUND

Point mutations in histone variant H3.3 (H3.3K27M, H3.3G34R) and the H3.3-specific ATRX/DAXX chaperone complex are frequent events in pediatric gliomas. These H3.3 point mutations affect many chromatin modifications but the exact oncogenic mechanisms are currently unclear. Histone H3.3 is known to localize to nuclear compartments known as promyelocytic leukemia (PML) nuclear bodies, which are frequently mutated and confirmed as oncogenic drivers in acute promyelocytic leukemia.

RESULTS

We find that the pediatric glioma-associated H3.3 point mutations disrupt the formation of PML nuclear bodies and this prevents differentiation down glial lineages. Similar to leukemias driven by PML mutations, H3.3-mutated glioma cells are sensitive to drugs that target PML bodies. We also find that point mutations in IDH1/2-which are common events in adult gliomas and myeloid leukemias-also disrupt the formation of PML bodies.

CONCLUSIONS

We identify PML as a contributor to oncogenesis in a subset of gliomas and show that targeting PML bodies is effective in treating these H3.3-mutated pediatric gliomas.

Genome sequencing and resequencing identified three horizontal gene transfers and uncovered the genetic mechanism on the intraspecies adaptive evolution of Gastrodia elata Blume

Horizontal gene transfer is a rare and useful genetic mechanism in higher plants. Gastrodia elata Blume (GE) (Orchidaceae), well known as traditional medicinal material in East Asia, adopts a heterotrophic lifestyle, thus being considered to be more prone to horizontal gene transfer (HGT). GE is a "polytypic species" that currently comprised of five recognized forms according to the plant morphology. G. elata Blume forma elata (GEE) and G. elata Bl.f.glauca (GEG) are two common forms that naturally grow in different habitats with difference in altitude and latitude. G. elata Bl.f.viridis (GEV) often occurs sporadically in cultivated populations of GEE and GEG. However, the genetic relationships and genetic mechanism underpinned the divergent ecological adaptations of GEE and GEG have not been revealed. Here, we assembled a chromosome-level draft genome of GEE with 1.04 Gb. Among predicted 17,895 protein coding genes, we identified three HGTs. Meanwhile, we resequenced 10 GEE accessions, nine GEG accessions, and 10 GEV accessions, and identified two independent genetic lineages: GEG_pedigree (GEG individuals and GEV individuals collected from GEG populations) and GEE_pedigree (GEE individuals and GEV individuals collected from GEE populations), which strongly support the taxonomic status of GEE and GEG as subspecies, not as different forms. In highly differentiated genomic regions of GEE_pedigree and GEG_pedigree, three chalcone synthase-encoding genes and one Phox/Bem1p (PB1) domain of encoding Auxin (AUX)/Indoleacetic acid (IAA) were identified in selection sweeping genome regions, which suggested that differentiation between GEE_pedigree and GEG_pedigree was promoted by the selection of genes related to photoresponse and growth and development. Overall, this new genome would be helpful for breeding and utilization of GE and the new findings would deepen the understanding about ecological adaptation and evolution of GE.

Interplay between activation of endogenous retroviruses and inflammation as common pathogenic mechanism in neurological and psychiatric disorders

Human endogenous retroviruses (ERVs) are ancestorial retroviral elements that were integrated into our genome through germline infections and insertions during evolution. They have repeatedly been implicated in the aetiology and pathophysiology of numerous human disorders, particularly in those that affect the central nervous system. In addition to the known association of ERVs with multiple sclerosis and amyotrophic lateral sclerosis, a growing number of studies links the induction and expression of these retroviral elements with the onset and severity of neurodevelopmental and psychiatric disorders. Although these disorders differ in terms of overall disease pathology and causalities, a certain degree of (subclinical) chronic inflammation can be identified in all of them. Based on these commonalities, we discuss the bidirectional relationship between ERV expression and inflammation and highlight that numerous entry points to this reciprocal sequence of events exist, including initial infections with ERV-activating pathogens, exposure to non-infectious inflammatory stimuli, and conditions in which epigenetic silencing of ERV elements is disrupted.

TERRA regulates DNA G-quadruplex formation and ATRX recruitment to chromatin

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.

The chromatin remodeller ATRX facilitates diverse nuclear processes, in a stochastic manner, in both heterochromatin and euchromatin

The chromatin remodeller ATRX interacts with the histone chaperone DAXX to deposit the histone variant H3.3 at sites of nucleosome turnover. ATRX is known to bind repetitive, heterochromatic regions of the genome including telomeres, ribosomal DNA and pericentric repeats, many of which are putative G-quadruplex forming sequences (PQS). At these sites ATRX plays an ancillary role in a wide range of nuclear processes facilitating replication, chromatin modification and transcription. Here, using an improved protocol for chromatin immunoprecipitation, we show that ATRX also binds active regulatory elements in euchromatin. Mutations in ATRX lead to perturbation of gene expression associated with a reduction in chromatin accessibility, histone modification, transcription factor binding and deposition of H3.3 at the sequences to which it normally binds. In erythroid cells where downregulation of α-globin expression is a hallmark of ATR-X syndrome, perturbation of chromatin accessibility and gene expression occurs in only a subset of cells. The stochastic nature of this process suggests that ATRX acts as a general facilitator of cell specific transcriptional and epigenetic programmes, both in heterochromatin and euchromatin.

Recent Bioinformatic Progress to Identify Epigenetic Changes Associated to Transposable Elements

Transposable elements (TEs) are recognized for their great impact on the functioning and evolution of their host genomes. They are associated to various deleterious effects, which has led to the evolution of regulatory epigenetic mechanisms to control their activity. Despite these negative effects, TEs are also important actors in the evolution of genomes by promoting genetic diversity and new regulatory elements. Consequently, it is important to study the epigenetic modifications associated to TEs especially at a locus-specific level to determine their individual influence on gene functioning. To this aim, this short review presents the current bioinformatic tools to achieve this task.

Histone H3.3 phosphorylation promotes heterochromatin formation by inhibiting H3K9/K36 histone demethylase

Histone H3.3 is an H3 variant which differs from the canonical H3.1/2 at four residues, including a serine residue at position 31 which is evolutionarily conserved. The H3.3 S31 residue is phosphorylated (H3.3 S31Ph) at heterochromatin regions including telomeres and pericentric repeats. However, the role of H3.3 S31Ph in these regions remains unknown. In this study, we find that H3.3 S31Ph regulates heterochromatin accessibility at telomeres during replication through regulation of H3K9/K36 histone demethylase KDM4B. In mouse embryonic stem (ES) cells, substitution of S31 with an alanine residue (H3.3 A31 -phosphorylation null mutant) results in increased KDM4B activity that removes H3K9me3 from telomeres. In contrast, substitution with a glutamic acid (H3.3 E31, mimics S31 phosphorylation) inhibits KDM4B, leading to increased H3K9me3 and DNA damage at telomeres. H3.3 E31 expression also increases damage at other heterochromatin regions including the pericentric heterochromatin and Y chromosome-specific satellite DNA repeats. We propose that H3.3 S31Ph regulation of KDM4B is required to control heterochromatin accessibility of repetitive DNA and preserve chromatin integrity.

CHAF1A/B mediate silencing of unintegrated HIV-1 DNAs early in infection

Early events of the retroviral life cycle are the targets of many host restriction factors that have evolved to prevent establishment of infection. Incoming retroviral DNAs are transcriptionally silenced before integration in most cell types, and efficient viral gene expression occurs only after formation of the provirus. The molecular machinery for silencing unintegrated retroviral DNAs of HIV-1 remains poorly characterized. Here, we identified the histone chaperones CHAF1A and CHAF1B as essential factors for silencing of unintegrated HIV-1 DNAs. Using RNAi-mediated knockdown (KD) of multiple histone chaperones, we found that KD of CHAF1A or CHAF1B resulted in a pronounced increase in expression of incoming viral DNAs. The function of these two proteins in silencing was independent of their interaction partner RBBP4. Viral DNA levels accumulated to significantly higher levels in CHAF1A KD cells over controls, suggesting enhanced stabilization of actively transcribed DNAs. Chromatin immunoprecipitation assays revealed no major changes in histone loading onto viral DNAs in the absence of CHAF1A, but levels of the H3K9 trimethylation silencing mark were reduced. KD of the H3K9me3-binding protein HP1γ accelerated the expression of unintegrated HIV-1 DNAs. While CHAF1A was critical for silencing HIV-1 DNAs, it showed no role in silencing of unintegrated retroviral DNAs of mouse leukemia virus. Our study identifies CHAF1A and CHAF1B as factors involved specifically in silencing of HIV-1 DNAs early in infection. The results suggest that these factors act by noncanonical pathways, distinct from their histone loading activities, to mediate silencing of newly synthesized HIV-1 DNAs.

Differential enrichment of H3K9me3 at annotated satellite DNA repeats in human cell lines and during fetal development in mouse

BACKGROUND

Trimethylation of histone H3 on lysine 9 (H3K9me3) at satellite DNA sequences has been primarily studied at (peri)centromeric regions, where its level shows differences associated with various processes such as development and malignant transformation. However, the dynamics of H3K9me3 at distal satellite DNA repeats has not been thoroughly investigated.

RESULTS

We exploit the sets of publicly available data derived from chromatin immunoprecipitation combined with massively parallel DNA sequencing (ChIP-Seq), produced by the The Encyclopedia of DNA Elements (ENCODE) project, to analyze H3K9me3 at assembled satellite DNA repeats in genomes of human cell lines and during mouse fetal development. We show that annotated satellite elements are generally enriched for H3K9me3, but its level in cancer cell lines is on average lower than in normal cell lines. We find 407 satellite DNA instances with differential H3K9me3 enrichment between cancer and normal cells including a large 115-kb cluster of GSATII elements on chromosome 12. Differentially enriched regions are not limited to satellite DNA instances, but instead encompass a wider region of flanking sequences. We found no correlation between the levels of H3K9me3 and noncoding RNA at corresponding satellite DNA loci. The analysis of data derived from multiple tissues identified 864 instances of satellite DNA sequences in the mouse reference genome that are differentially enriched between fetal developmental stages.

CONCLUSIONS

Our study reveals significant differences in H3K9me3 level at a subset of satellite repeats between biological states and as such contributes to understanding of the role of satellite DNA repeats in epigenetic regulation during development and carcinogenesis.

UHRF1 regulates alternative splicing by binding to splicing factors and U snRNAs

The well-established functions of UHRF1 converge to DNA biological processes, as exemplified by DNA methylation maintenance and DNA damage repair during cell cycles. However, the potential effect of UHRF1 on RNA metabolism is largely unexplored. Here, we revealed that UHRF1 serves as a novel alternative RNA splicing regulator. The protein interactome of UHRF1 identified various splicing factors. Among them, SF3B3 could interact with UHRF1 directly and participate in UHRF1-regulated alternative splicing events. Furthermore, we interrogated the RNA interactome of UHRF1, and surprisingly, we identified U snRNAs, the canonical spliceosome component, in the UHRF1-associated RNA components. Unexpectedly, we found H3R2 methylation status determines the binding preference of U snRNA, especially U2 snRNA. The involvement of U snRNAs in UHRF1-containing complex and their binding preference to specific chromatin configuration imply a finely orchestrated mechanism at play. Our results provided the resources and pinpointed the molecular basis of UHRF1-mediated alternative RNA splicing, which will help us better our understanding of the physiological and pathological roles of UHRF1 in disease development.

Mutations inhibiting KDM4B drive ALT activation in ATRX-mutated glioblastomas

Alternative Lengthening of Telomeres (ALT) is a telomere maintenance pathway utilised in 15% of cancers. ALT cancers are strongly associated with inactivating mutations in ATRX; yet loss of ATRX alone is insufficient to trigger ALT, suggesting that additional cooperating factors are involved. We identify H3.3^G34R and IDH1/2 mutations as two such factors in ATRX-mutated glioblastomas. Both mutations are capable of inactivating histone demethylases, and we identify KDM4B as the key demethylase inactivated in ALT. Mouse embryonic stem cells inactivated for ATRX, TP53, TERT and KDM4B (KDM4B knockout or H3.3^G34R) show characteristic features of ALT. Conversely, KDM4B over-expression in ALT cancer cells abrogates ALT-associated features. In this work, we demonstrate that inactivation of KDM4B, through H3.3^G34R or IDH1/2 mutations, acts in tandem with ATRX mutations to promote ALT in glioblastomas.

Alternative Lengthening of Telomeres (ALT) is a telomere maintenance pathway utilised in 15% of cancers that have been associated with mutations in ATRX. Here the authors reveal a functional role of histone demethylases KDM4B in regulating ALT activation.

Innovation of heterochromatin functions drives rapid evolution of essential ZAD-ZNF genes in Drosophila

Contrary to dogma, evolutionarily young and dynamic genes can encode essential functions. We find that evolutionarily dynamic ZAD-ZNF genes, which encode the most abundant class of insect transcription factors, are more likely to encode essential functions in Drosophila melanogaster than ancient, conserved ZAD-ZNF genes. We focus on the Nicknack ZAD-ZNF gene, which is evolutionarily young, poorly retained in Drosophila species, and evolves under strong positive selection. Yet we find that it is necessary for larval development in D. melanogaster. We show that Nicknack encodes a heterochromatin-localizing protein like its paralog Oddjob, also an evolutionarily dynamic yet essential ZAD-ZNF gene. We find that the divergent D. simulans Nicknack protein can still localize to D. melanogaster heterochromatin and rescue viability of female but not male Nicknack-null D. melanogaster. Our findings suggest that innovation for rapidly changing heterochromatin functions might generally explain the essentiality of many evolutionarily dynamic ZAD-ZNF genes in insects.

Retinoblastoma Tumor Suppressor Protein Roles in Epigenetic Regulation

Simple Summary

Loss of function of the retinoblastoma gene (RB1) is the rate-limiting step in the initiation of both the hereditary and sporadic forms of retinoblastoma tumor. Furthermore, loss of function of the retinoblastoma tumor suppressor protein (pRB) is frequently found in most human cancers. In retinoblastoma, tumor progression is driven by epigenetic changes following pRB loss. This review focuses on the diverse functions of pRB in epigenetic regulation.

Abstract

Mutations that result in the loss of function of pRB were first identified in retinoblastoma and since then have been associated with the propagation of various forms of cancer. pRB is best known for its key role as a transcriptional regulator during cell cycle exit. Beyond the ability of pRB to regulate transcription of cell cycle progression genes, pRB can remodel chromatin to exert several of its other biological roles. In this review, we discuss the diverse functions of pRB in epigenetic regulation including nucleosome mobilization, histone modifications, DNA methylation and non-coding RNAs.

Silencing and Transcriptional Regulation of Endogenous Retroviruses: An Overview

Almost half of the human genome is made up of transposable elements (TEs), and about 8% consists of endogenous retroviruses (ERVs). ERVs are remnants of ancient exogenous retrovirus infections of the germ line. Most TEs are inactive and not detrimental to the host. They are tightly regulated to ensure genomic stability of the host and avoid deregulation of nearby gene loci. Histone-based posttranslational modifications such as H3K9 trimethylation are one of the main silencing mechanisms. Trim28 is one of the identified master regulators of silencing, which recruits most prominently the H3K9 methyltransferase Setdb1, among other factors. Sumoylation and ATP-dependent chromatin remodeling factors seem to contribute to proper localization of Trim28 to ERV sequences and promote Trim28 interaction with Setdb1. Additionally, DNA methylation as well as RNA-mediated targeting of TEs such as piRNA-based silencing play important roles in ERV regulation. Despite the involvement of ERV overexpression in several cancer types, autoimmune diseases, and viral pathologies, ERVs are now also appreciated for their potential positive role in evolution. ERVs can provide new regulatory gene elements or novel binding sites for transcription factors, and ERV gene products can even be repurposed for the benefit of the host.

An atlas of dynamic chromatin landscapes in mouse fetal development

The Encyclopedia of DNA Elements (ENCODE) project has established a genomic resource for mammalian development, profiling a diverse panel of mouse tissues at 8 developmental stages from 10.5 days after conception until birth, including transcriptomes, methylomes and chromatin states. Here we systematically examined the state and accessibility of chromatin in the developing mouse fetus. In total we performed 1,128 chromatin immunoprecipitation with sequencing (ChIP-seq) assays for histone modifications and 132 assay for transposase-accessible chromatin using sequencing (ATAC-seq) assays for chromatin accessibility across 72 distinct tissue-stages. We used integrative analysis to develop a unified set of chromatin state annotations, infer the identities of dynamic enhancers and key transcriptional regulators, and characterize the relationship between chromatin state and accessibility during developmental gene regulation. We also leveraged these data to link enhancers to putative target genes and demonstrate tissue-specific enrichments of sequence variants associated with disease in humans. The mouse ENCODE data sets provide a compendium of resources for biomedical researchers and achieve, to our knowledge, the most comprehensive view of chromatin dynamics during mammalian fetal development to date.

A role for the rare endogenous retrovirus β4 in development of Japanese fancy mice

Two coat-color mutations, nonagouti, which changes coat color from wild-type agouti to black, and piebald, which induces irregular white spotting, are the characteristics of Japanese fancy mouse strain JF1/Ms. In our Communications Biology article, we reported that insertion of a rare type of endogenous retrovirus β4 has caused both coat color mutations. Although there are some reports on the roles of β4 in the mouse genome, further studies on β4 will uncover new features of endogenous retrovirus sequences.

In light of their recent finding that insertion of a rare endogenous retrovirus, β4, is the cause of the characteristic coat coloring in agouti and piebald mice, Akira Tanave and Tsuyoshi Koide now discuss the origin and expansion of this element as well as potential roles of β4 in the mouse genome.

SIRT7 mediates L1 elements transcriptional repression and their association with the nuclear lamina

Long interspersed elements-1 (LINE-1, L1) are retrotransposons that hold the capacity of self-propagation in the genome with potential mutagenic outcomes. How somatic cells restrict L1 activity and how this process becomes dysfunctional during aging and in cancer cells is poorly understood. L1s are enriched at lamin-associated domains, heterochromatic regions of the nuclear periphery. Whether this association is necessary for their repression has been elusive. Here we show that the sirtuin family member SIRT7 participates in the epigenetic transcriptional repression of L1 genome-wide in both mouse and human cells. SIRT7 depletion leads to increased L1 expression and retrotransposition. Mechanistically, we identify a novel interplay between SIRT7 and Lamin A/C in L1 repression. Our results demonstrate that SIRT7-mediated H3K18 deacetylation regulates L1 expression and promotes L1 association with elements of the nuclear lamina. The failure of such activity might contribute to the observed genome instability and compromised viability in SIRT7 knockout mice. Overall, our results reveal a novel function of SIRT7 on chromatin organization by mediating the anchoring of L1 to the nuclear envelope, and a new functional link of the nuclear lamina with transcriptional repression.

Telescope: Characterization of the retrotranscriptome by accurate estimation of transposable element expression

Characterization of Human Endogenous Retrovirus (HERV) expression within the transcriptomic landscape using RNA-seq is complicated by uncertainty in fragment assignment because of sequence similarity. We present Telescope, a computational software tool that provides accurate estimation of transposable element expression (retrotranscriptome) resolved to specific genomic locations. Telescope directly addresses uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We demonstrate the utility of our approach through single locus analysis of HERV expression in 13 ENCODE cell types. When examined at this resolution, we find that the magnitude and breadth of the retrotranscriptome can be vastly different among cell types. Furthermore, our approach is robust to differences in sequencing technology and demonstrates that the retrotranscriptome has potential to be used for cell type identification. We compared our tool with other approaches for quantifying transposable element (TE) expression, and found that Telescope has the greatest resolution, as it estimates expression at specific TE insertions rather than at the TE subfamily level. Telescope performs highly accurate quantification of the retrotranscriptomic landscape in RNA-seq experiments, revealing a differential complexity in the transposable element biology of complex systems not previously observed. Telescope is available at https://github.com/mlbendall/telescope.

Almost half of the human genome is composed of transposable elements (TEs), but their contribution to the transcriptome, their cell-type specific expression patterns, and their role in disease remains poorly understood. Recent studies have found many elements to be actively expressed and involved in key cellular processes. For example, human endogenous retroviruses (HERVs) are reported to be involved in human embryonic stem cell differentiation. Discovering which exact HERVs are differentially expressed in RNA-seq data would be a major advance in understanding such processes. However, because HERVs have a high level of sequence similarity it is hard to identify which exact HERV is differentially expressed. To solve this problem, we developed a computer program which addressed uncertainty in fragment assignment by reassigning ambiguously mapped fragments to the most probable source transcript as determined within a Bayesian statistical model. We call this program, “Telescope”. We then used Telescope to identify HERV expression in 13 well-studied cell types from the ENCODE consortium and found that different cell types could be characterized by enrichment for different HERV families, and for locus specific expression. We also showed that Telescope performed better than other methods currently used to determine TE expression. The use of this computational tool to examine new and existing RNA-seq data sets may lead to new understanding of the roles of TEs in health and disease.

RB1 Deletion in Retinoblastoma Protein Pathway-Disrupted Cells Results in DNA Damage and Cancer Progression

Proliferative control in cancer cells is frequently disrupted by mutations in the retinoblastoma protein (RB) pathway. Intriguingly, RB1 mutations can arise late in tumorigenesis in cancer cells whose RB pathway is already compromised by another mutation. In this study, we present evidence for increased DNA damage and instability in cancer cells with RB pathway defects when RB1 mutations are induced. We generated isogenic RB1 mutant genotypes with CRISPR/Cas9 in a number of cell lines. Cells with even one mutant copy of RB1 have increased basal levels of DNA damage and increased mitotic errors. Elevated levels of reactive oxygen species as well as impaired homologous recombination repair underlie this DNA damage. When xenografted into immunocompromised mice, RB1 mutant cells exhibit an elevated propensity to seed new tumors in recipient lungs. This study offers evidence that late-arising RB1 mutations can facilitate genome instability and cancer progression that are beyond the preexisting proliferative control deficit.

Paired-end mappability of transposable elements in the human genome

Though transposable elements make up around half of the human genome, the repetitive nature of their sequences makes it difficult to accurately align conventional sequencing reads. However, in light of new advances in sequencing technology, such as increased read length and paired-end libraries, these repetitive regions are now becoming easier to align to. This study investigates the mappability of transposable elements with 50 bp, 76 bp and 100 bp paired-end read libraries. With respect to those read lengths and allowing for 3 mismatches during alignment, over 68, 85, and 88% of all transposable elements in the RepeatMasker database are uniquely mappable, suggesting that accurate locus-specific mapping of older transposable elements is well within reach.

Transposable Elements and Their Epigenetic Regulation in Mental Disorders: Current Evidence in the Field

Transposable elements (TEs) are highly repetitive DNA sequences in the human genome that are the relics of previous retrotransposition events. Although the majority of TEs are transcriptionally inactive due to acquired mutations or epigenetic processes, around 8% of TEs exert transcriptional activity. It has been found that TEs contribute to somatic mosaicism that accounts for functional specification of various brain cells. Indeed, autonomous retrotransposition of long interspersed element-1 (LINE-1) sequences has been reported in the neural rat progenitor cells from the hippocampus, the human fetal brain and the human embryonic stem cells. Moreover, expression of TEs has been found to regulate immune-inflammatory responses, conditioning immunity against exogenous infections. Therefore, aberrant epigenetic regulation and expression of TEs emerged as a potential mechanism underlying the development of various mental disorders, including autism spectrum disorders (ASD), schizophrenia, bipolar disorder, major depression, and Alzheimer's disease (AD). Consequently, some studies revealed that expression of some sequences of human endogenous retroviruses (HERVs) appears only in a certain group of patients with mental disorders (especially those with schizophrenia, bipolar disorder, and ASD) but not in healthy controls. In addition, it has been found that expression of HERVs might be related to subclinical inflammation observed in mental disorders. In this article, we provide an overview of detrimental effects of transposition on the brain development and immune mechanisms with relevance to mental disorders. We show that transposition is not the only mechanism, explaining the way TEs might shape the phenotype of mental disorders. Other mechanisms include the regulation of gene expression and the impact on genomic stability. Next, we review current evidence from studies investigating expression and epigenetic regulation of specific TEs in various mental disorders. Most consistently, these studies indicate altered expression of HERVs and methylation of LINE-1 sequences in patients with ASD, schizophrenia, and mood disorders. However, the contribution of TEs to the etiology of AD is poorly documented. Future studies should further investigate the mechanisms linking epigenetic processes, specific TEs and the phenotype of mental disorders to disentangle causal associations.

Early life exposures, neurodevelopmental disorders, and transposable elements

Transposable elements make up a much larger portion of the genome than protein-coding genes, yet we know relatively little about their function in the human genome. However, we are beginning to more fully understand their role in brain development, neuroinflammation, and adaptation to environmental insults such as stress. For instance, glucocorticoid receptor activation regulates transposable elements in the brain following acute stress. Early life is a period of substantial brain development during which transposable elements play a role. Environmental exposures and experiences during early life that promote abnormal regulation of transposable elements may lead to a cascade of events that ultimately increase susceptibility to disorders later in life. Recent attention to transposable elements in psychiatric illness has begun to clarify associations indicative of dysregulation of different classes of transposable elements in stress-related and neurodevelopmental illness. Though individual susceptibility or resiliency to psychiatric illness has not been explained by traditional genetic studies, the wide inter-individual variability in transposable element composition in the human genome make TEs attractive candidates to elucidate this differential susceptibility. In this review, we discuss evidence that regulation of transposable elements in the brain are stage-specific, sensitive to environmental factors, and may be impacted by early life perturbations. We further present evidence of associations with stress-related and neurodevelopmental psychiatric illness from a developmental perspective.

On the Importance to Acknowledge Transposable Elements in Epigenomic Analyses

Eukaryotic genomes comprise a large proportion of repeated sequences, an important fraction of which are transposable elements (TEs). TEs are mobile elements that have a significant impact on genome evolution and on gene functioning. Although some TE insertions could provide adaptive advantages to species, transposition is a highly mutagenic event that has to be tightly controlled to ensure its viability. Genomes have evolved sophisticated mechanisms to control TE activity, the most important being epigenetic silencing. However, the epigenetic control of TEs can also affect genes located nearby that can become epigenetically regulated. It has been proposed that the combination of TE mobilization and the induced changes in the epigenetic landscape could allow a rapid phenotypic adaptation to global environmental changes. In this review, we argue the crucial need to take into account the repeated part of genomes when studying the global impact of epigenetic modifications on an organism. We emphasize more particularly why it is important to carefully consider TEs and what bioinformatic tools can be used to do so.

PML modulates H3.3 targeting to telomeric and centromeric repeats in mouse fibroblasts

Targeted deposition of histone variant H3.3 into chromatin is paramount for proper regulation of chromatin integrity, particularly in heterochromatic regions including repeats. We have recently shown that the promyelocytic leukemia (PML) protein prevents H3.3 from being deposited in large heterochromatic PML-associated domains (PADs). However, to what extent PML modulates H3.3 loading on chromatin in other areas of the genome remains unexplored. Here, we examined the impact of PML on targeting of H3.3 to genes and repeat regions that reside outside PADs. We show that loss of PML increases H3.3 deposition in subtelomeric, telomeric, pericentric and centromeric repeats in mouse embryonic fibroblasts, while other repeat classes are not affected. Expression of major satellite, minor satellite and telomeric non-coding transcripts is altered in Pml-null cells. In particular, telomeric Terra transcripts are strongly upregulated, in concordance with a marked reduction in H4K20me3 at these sites. Lastly, for most genes H3.3 enrichment or gene expression outcomes are independent of PML. Our data argue towards the importance of a PML-H3.3 axis in preserving a heterochromatin state at centromeres and telomeres.

Transposable elements are regulated by context-specific patterns of chromatin marks in mouse embryonic stem cells

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs.

Transposable elements (TEs) fulfill essential but poorly understood roles in genome organization and gene expression control. Here the authors show that the regulation of TEs occurs through overlapping epigenetic mechanisms that control the expression and chromatin signatures at TEs.

Screening and Identification of PLK1-Polo Box Binding Peptides by High-Throughput Sequencing of Phage-Selected Libraries

BACKGROUND

Human proteome contains a plethora of short linear peptide motifs that is crucial for signaling and other cellular processes. These motifs are difficult to identify due to lack of systematic approach for their detection.

OBJECTIVES

Here we demonstrate the use of peptide phage display in combination with high throughput next generation sequencing to identify enriched peptide sequences through biopanning process against polo box domain (PBD) of mitotic polo like kinase 1 (Plk1).

METHODS

Purified recombinant Plk1 and two unrelated controls namely B-lymphocyte antigen (CD20) and fluorescent protein (mCherry) were subjected to peptide phage display analysis. Bacterially-propagated phage DNA was amplified by PCR using triplet bar coded primers to tag the pool from each amplicon.

RESULTS

Proteomic peptide phage display along with next generation sequencing and Bioinformatics analysis demonstrated several known and putative novel interactions which were potentially related to Plk1-PBD. With our strategy, we were able to identify and characterize several Plk1-PBD binding peptides, as well as define more precisely, consensus sequences.

CONCLUSION

We believe that this information could provide valuable tools for exploring novel interaction involved in Plk1 signaling as well as to choose peptides for Plk1 specific drug development.

Widespread roles of enhancer-like transposable elements in cell identity and long-range genomic interactions

A few families of transposable elements (TEs) have been shown to evolve into cis-regulatory elements (CREs). Here, to extend these studies to all classes of TEs in the human genome, we identified widespread enhancer-like repeats (ELRs) and find that ELRs reliably mark cell identities, are enriched for lineage-specific master transcription factor binding sites, and are mostly primate-specific. In particular, elements of MIR and L2 TE families whose abundance co-evolved across chordate genomes, are found as ELRs in most human cell types examined. MIR and L2 elements frequently share long-range intra-chromosomal interactions and binding of physically interacting transcription factors. We validated that eight L2 and nine MIR elements function as enhancers in reporter assays, and among 20 MIR-L2 pairings, one MIR repressed and one boosted the enhancer activity of L2 elements. Our results reveal a previously unappreciated co-evolution and interaction between two TE families in shaping regulatory networks.

Ribosomal DNA copy loss and repeat instability in ATRX-mutated cancers

ATRX (alpha thalassemia/mental retardation X-linked) complexes with DAXX to deposit histone variant H3.3 into repetitive heterochromatin. Recent genome sequencing studies in cancers have revealed mutations in ATRX and their association with ALT (alternative lengthening of telomeres) activation. Here we report depletion of ATRX in mouse ES cells leads to selective loss in ribosomal RNA gene (rDNA) copy number. Supporting this, ATRX-mutated human ALT-positive tumors also show a substantially lower rDNA copy than ALT-negative tumors. Further investigation shows that the rDNA copy loss and repeat instability are caused by a disruption in H3.3 deposition and thus a failure in heterochromatin formation at rDNA repeats in the absence of ATRX. We also find that ATRX-depleted cells are reduced in ribosomal RNA transcription output and show increased sensitivity to RNA polymerase I (Pol I) transcription inhibitor CX5461. In addition, human ALT-positive cancer cell lines are also more sensitive to CX5461 treatment. Our study provides insights into the contribution of ATRX loss of function to tumorigenesis through the loss of rDNA stability and suggests the therapeutic potential of targeting Pol I transcription in ALT cancers.

Silencing of transposable elements may not be a major driver of regulatory evolution in primate iPSCs

Transposable elements (TEs) comprise almost half of primate genomes and their aberrant regulation can result in deleterious effects. In pluripotent stem cells, rapidly evolving KRAB-ZNF genes target TEs for silencing by H3K9me3. To investigate the evolution of TE silencing, we performed H3K9me3 ChIP-seq experiments in induced pluripotent stem cells from 10 human and 7 chimpanzee individuals. We identified four million orthologous TEs and found the SVA and ERV families to be marked most frequently by H3K9me3. We found little evidence of inter-species differences in TE silencing, with as many as 82% of putatively silenced TEs marked at similar levels in humans and chimpanzees. TEs that are preferentially silenced in one species are a similar age to those silenced in both species and are not more likely to be associated with expression divergence of nearby orthologous genes. Our data suggest limited species-specificity of TE silencing across 6 million years of primate evolution.

The epigenetic alterations of endogenous retroelements in aging

Endogenous retroelements, transposons that mobilize through RNA intermediates, include some of the most abundant repetitive sequences of the human genome, such as Alu and LINE-1 sequences, and human endogenous retroviruses. Recent discoveries demonstrate that these mobile genetic elements not only act as intragenomic parasites, but also exert regulatory roles in living cells. The risk of genomic instability represented by endogenous retroelements is normally counteracted by a series of epigenetic control mechanisms which include, among the most important, CpG DNA methylation. Indeed, most of the genomic CpG sites subjected to DNA methylation in the nuclear DNA are carried by these repetitive elements. As other parts of the genome, endogenous retroelements and other transposable elements are subjected to deep epigenetic alterations during aging, repeatedly observed in the context of organismal and cellular senescence, in human and other species. This review summarizes the current status of knowledge about the epigenetic alterations occurring in this large, non-genic portion of the genome in aging and age-related conditions, with a focus on the causes and the possible functional consequences of these alterations.

Comparison of ChIP-Seq Data and a Reference Motif Set for Human KRAB C2H2 Zinc Finger Proteins

KRAB C2H2 zinc finger proteins (KZNFs) are the largest and most diverse family of human transcription factors, likely due to diversifying selection driven by novel endogenous retroelements (EREs), but the vast majority lack binding motifs or functional data. Two recent studies analyzed a majority of the human KZNFs using either ChIP-seq (60 proteins) or ChIP-exo (221 proteins) in the same cell type (HEK293). The ChIP-exo paper did not describe binding motifs, however. Thirty-nine proteins are represented in both studies, enabling the systematic comparison of the data sets presented here. Typically, only a minority of peaks overlap, but the two studies nonetheless display significant similarity in ERE binding for 32/39, and yield highly similar DNA binding motifs for 23 and related motifs for 34 (MoSBAT similarity score >0.5 and >0.2, respectively). Thus, there is overall (albeit imperfect) agreement between the two studies. For the 242 proteins represented in at least one study, we selected a highest-confidence motif for each protein, utilizing several motif-derivation approaches, and evaluating motifs within and across data sets. Peaks for the majority (158) are enriched (96% with AUC >0.6 predicting peak vs. nonpeak) for a motif that is supported by the C2H2 “recognition code,” consistent with intrinsic sequence specificity driving DNA binding in cells. An additional 63 yield motifs enriched in peaks, but not supported by the recognition code, which could reflect indirect binding. Altogether, these analyses validate both data sets, and provide a reference motif set with associated quality metrics.

How to Tackle Challenging ChIP-Seq, with Long-Range Cross-Linking, Using ATRX as an Example

Chromatin immunoprecipitation coupled with high-throughput, next-generation DNA sequencing (ChIP-seq) has enabled researchers to establish the genome-wide patterns of chromatin modifications and binding of chromatin-associated proteins. Well-established protocols produce robust ChIP-seq data for many proteins by sequencing the DNA obtained following immunoprecipitation of fragmented chromatin using a wide range of specific antibodies. In general, the quality of these data mainly depends on the specificity and avidity of the antibody used. However, even using optimal antibodies, ChIP-seq can become more challenging when the protein associates with chromatin via protein-protein interactions rather than directly binding DNA. An example of such a protein is the alpha-thalassaemia mental retardation X-linked (ATRX) protein; a chromatin remodeler that associates with the histone chaperone DAXX, in the deposition of the replication-independent histone variant H3.3 and plays an important role in maintaining chromatin integrity. Inherited mutations of ATRX cause syndromal mental retardation (ATR-X Syndrome) whereas acquired mutations are associated with myelodysplasia, acute myeloid leukemia (ATMDS syndrome), and a range of solid tumors. Therefore, high quality ChIP-seq data have been needed to analyze the genome-wide distribution of ATRX, to advance our understanding of its normal role and to comprehend how mutations contribute to human disease. Here, we describe an optimized ChIP-seq protocol for ATRX which can also be used to produce high quality data sets for other challenging proteins which are indirectly associated with DNA and complement the ChIP-seq toolkit for genome-wide analyses of histone chaperon complexes and associated chromatin remodelers. Although not a focus of this chapter, we will also provide some insight for the analysis of the large dataset generated by ChIP-seq. Even though this protocol has been fully optimized for ATRX, it should also provide guidance for efficient ChIP-seq analysis, using the appropriate antibodies, for other proteins interacting indirectly with DNA.

Analysis of RNA-Seq Data Using TEtranscripts

Transposable elements (TE) are mobile genetic elements that can readily change their genomic position. When not properly silenced, TEs can contribute a substantial portion to the cell's transcriptome, but are typically ignored in most RNA-seq data analyses. One reason for leaving TE-derived reads out of RNA-seq analyses is the complexities involved in properly aligning short sequencing reads to these highly repetitive regions. Here we describe a method for including TE-derived reads in RNA-seq differential expression analysis using an open source software package called TEtranscripts. TEtranscripts is designed to assign both uniquely and ambiguously mapped reads to all possible gene and TE-derived transcripts in order to statistically infer the correct gene/TE abundances. Here, we provide a detailed tutorial of TEtranscripts using a published qPCR validated dataset.Barbara McClintock laid the foundation for TE research with her discoveries in maize of mobile genetic elements capable of inserting into novel locations in the genome, altering the expression of nearby genes [1]. Since then, our appreciation of the contribution of repetitive TE-derived sequences to eukaryotic genomes has vastly increased. With the publication of the first human genome draft by the Human Genome Project, it was determined that nearly half of the human genome is derived from TE sequences [2, 3], with varying levels of repetitive DNA present in most plant and animal species. More recent studies looking at distantly related TE-like sequences have estimated that up to two thirds of the human genome might be repeat-derived [4], with the vast majority of these sequences attributed to retrotransposons that require transcription as part of the mobilization process, as discussed below.

Recent advances in ChIP-seq analysis: from quality management to whole-genome annotation

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis can detect protein/DNA-binding and histone-modification sites across an entire genome. Recent advances in sequencing technologies and analyses enable us to compare hundreds of samples simultaneously; such large-scale analysis has potential to reveal the high-dimensional interrelationship level for regulatory elements and annotate novel functional genomic regions de novo. Because many experimental considerations are relevant to the choice of a method in a ChIP-seq analysis, the overall design and quality management of the experiment are of critical importance. This review offers guiding principles of computation and sample preparation for ChIP-seq analyses, highlighting the validity and limitations of the state-of-the-art procedures at each step. We also discuss the latest challenges of single-cell analysis that will encourage a new era in this field.

PML protein organizes heterochromatin domains where it regulates histone H3.3 deposition by ATRX/DAXX

Maintenance of chromatin homeostasis involves proper delivery of histone variants to the genome. The interplay between different chaperones regulating the supply of histone variants to distinct chromatin domains remains largely undeciphered. We report a role of promyelocytic leukemia (PML) protein in the routing of histone variant H3.3 to chromatin and in the organization of megabase-size heterochromatic PML-associated domains that we call PADs. Loss of PML alters the heterochromatic state of PADs by shifting the histone H3 methylation balance from K9me3 to K27me3. Loss of PML impairs deposition of H3.3 by ATRX and DAXX in PADs but preserves the H3.3 loading function of HIRA in these regions. Our results unveil an unappreciated role of PML in the large-scale organization of chromatin and demonstrate a PML-dependent role of ATRX/DAXX in the deposition of H3.3 in PADs. Our data suggest that H3.3 loading by HIRA and ATRX-dependent H3K27 trimethylation constitute mechanisms ensuring maintenance of heterochromatin when the integrity of these domains is compromised.

Polycomb repressive complex 1 provides a molecular explanation for repeat copy number dependency in FSHD muscular dystrophy

Repression of repetitive elements is crucial to preserve genome integrity and has been traditionally ascribed to constitutive heterochromatin pathways. FacioScapuloHumeral Muscular Dystrophy (FSHD), one of the most common myopathies, is characterized by a complex interplay of genetic and epigenetic events. The main FSHD form is linked to a reduced copy number of the D4Z4 macrosatellite repeat on 4q35, causing loss of silencing and aberrant expression of the D4Z4-embedded DUX4 gene leading to disease. By an unknown mechanism, D4Z4 copy-number correlates with FSHD phenotype. Here we show that the DUX4 proximal promoter (DUX4p) is sufficient to nucleate the enrichment of both constitutive and facultative heterochromatin components and to mediate a copy-number dependent gene silencing. We found that both the CpG/GC dense DNA content and the repetitive nature of DUX4p arrays are important for their repressive ability. We showed that DUX4p mediates a copy number-dependent Polycomb Repressive Complex 1 (PRC1) recruitment, which is responsible for the copy-number dependent gene repression. Overall, we directly link genetic and epigenetic defects in FSHD by proposing a novel molecular explanation for the copy number-dependency in FSHD pathogenesis, and offer insight into the molecular functions of repeats in chromatin regulation.

The molecular basis of the organization of repetitive DNA-containing constitutive heterochromatin in mammals

Constitutive heterochromatin is composed mainly of repetitive elements and represents the typical inert chromatin structure in eukaryotic cells. Approximately half of the mammalian genome is made of repeat sequences, such as satellite DNA, telomeric DNA, and transposable elements. As essential genes are not present in these regions, most of these repeat sequences were considered as junk DNA in the past. However, it is now clear that these regions are essential for chromosome stability and the silencing of neighboring genes. Genetic and biochemical studies have revealed that histone methylation at H3K9 and its recognition by heterochromatin protein 1 represent the fundamental mechanism by which heterochromatin forms. Although this molecular mechanism is highly conserved from yeast to human cells, its detailed epigenetic regulation is more complex and dynamic for each distinct constitutive heterochromatin structure in higher eukaryotes. It can also vary according to the developmental stage. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) analysis is a powerful tool to investigate the epigenetic regulation of eukaryote genomes, but non-unique reads are usually discarded during standard ChIP-seq data alignment to reference genome databases. Therefore, specific methods to obtain global epigenetic information concerning repetitive elements are needed. In this review, we focus on such approaches and we summarize the latest molecular models for distinct constitutive heterochromatin types in mammals.

An RB-EZH2 Complex Mediates Silencing of Repetitive DNA Sequences

Repetitive genomic regions include tandem sequence repeats and interspersed repeats, such as endogenous retroviruses and LINE-1 elements. Repressive heterochromatin domains silence expression of these sequences through mechanisms that remain poorly understood. Here, we present evidence that the retinoblastoma protein (pRB) utilizes a cell-cycle-independent interaction with E2F1 to recruit enhancer of zeste homolog 2 (EZH2) to diverse repeat sequences. These include simple repeats, satellites, LINEs, and endogenous retroviruses as well as transposon fragments. We generated a mutant mouse strain carrying an F832A mutation in Rb1 that is defective for recruitment to repetitive sequences. Loss of pRB-EZH2 complexes from repeats disperses H3K27me3 from these genomic locations and permits repeat expression. Consistent with maintenance of H3K27me3 at the Hox clusters, these mice are developmentally normal. However, susceptibility to lymphoma suggests that pRB-EZH2 recruitment to repetitive elements may be cancer relevant.

Quantitative analysis of ChIP-seq data uncovers dynamic and sustained H3K4me3 and H3K27me3 modulation in cancer cells under hypoxia

Background

A comprehensive assessment of the epigenetic dynamics in cancer cells is the key to understanding the molecular mechanisms underlying cancer and to improving cancer diagnostics, prognostics and treatment. By combining genome-wide ChIP-seq epigenomics and microarray transcriptomics, we studied the effects of oxygen deprivation and subsequent reoxygenation on histone 3 trimethylation of lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in a breast cancer cell line, serving as a model for abnormal oxygenation in solid tumors. A priori, epigenetic markings and gene expression levels not only are expected to vary greatly between hypoxic and normoxic conditions, but also display a large degree of heterogeneity across the cell population. Where traditionally ChIP-seq data are often treated as dichotomous data, the model and experiment here necessitate a quantitative, data-driven analysis of both datasets.

Results

We first identified genomic regions with sustained epigenetic markings, which provided a sample-specific reference enabling quantitative ChIP-seq data analysis. Sustained H3K27me3 marking was located around centromeres and intergenic regions, while sustained H3K4me3 marking is associated with genes involved in RNA binding, translation and protein transport and localization. Dynamic marking with both H3K4me3 and H3K27me3 (hypoxia-induced bivalency) was found in CpG-rich regions at loci encoding factors that control developmental processes, congruent with observations in embryonic stem cells.

Conclusions

In silico-identified epigenetically sustained and dynamic genomic regions were confirmed through ChIP-PCR in vitro, and obtained results are corroborated by published data and current insights regarding epigenetic regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0090-4) contains supplementary material, which is available to authorized users.

Setdb1 maintains hematopoietic stem and progenitor cells by restricting the ectopic activation of nonhematopoietic genes

Setdb1, an H3K9 histone methyltransferase, is essential for the maintenance of HSPCs. Setdb1 restricts the activation of nonhematopoietic genes, such as gluconeogenic pathway genes, to maintain HSPCs.

Detection of differential DNA methylation in repetitive DNA of mice and humans perinatally exposed to bisphenol A

Developmental exposure to bisphenol A (BPA) has been shown to induce changes in DNA methylation in both mouse and human genic regions; however, the response in repetitive elements and transposons has not been explored. Here we present novel methodology to combine genomic DNA enrichment with RepeatMasker analysis on next-generation sequencing data to determine the effect of perinatal BPA exposure on repetitive DNA at the class, family, subfamily, and individual insertion level in both mouse and human samples. Mice were treated during gestation and lactation to BPA in chow at 0, 50, or 50,000 ng/g levels and total BPA was measured in stratified human fetal liver tissue samples as low (non-detect to 0.83 ng/g), medium (3.5 to 5.79 ng/g), or high (35.44 to 96.76 ng/g). Transposon methylation changes were evident in human classes, families, and subfamilies, with the medium group exhibiting hypomethylation compared to both high and low BPA groups. Mouse repeat classes, families, and subfamilies did not respond to BPA with significantly detectable differential DNA methylation. In human samples, 1251 individual transposon loci were detected as differentially methylated by BPA exposure, but only 19 were detected in mice. Of note, this approach recapitulated the discovery of a previously known mouse environmentally labile metastable epiallele, Cabp^IAP. Thus, by querying repetitive DNA in both mouse and humans, we report the first known transposons in humans that respond to perinatal BPA exposure.

Integration and Fixation Preferences of Human and Mouse Endogenous Retroviruses Uncovered with Functional Data Analysis

Endogenous retroviruses (ERVs), the remnants of retroviral infections in the germ line, occupy ~8% and ~10% of the human and mouse genomes, respectively, and affect their structure, evolution, and function. Yet we still have a limited understanding of how the genomic landscape influences integration and fixation of ERVs. Here we conducted a genome-wide study of the most recently active ERVs in the human and mouse genome. We investigated 826 fixed and 1,065 in vitro HERV-Ks in human, and 1,624 fixed and 242 polymorphic ETns, as well as 3,964 fixed and 1,986 polymorphic IAPs, in mouse. We quantitated >40 human and mouse genomic features (e.g., non-B DNA structure, recombination rates, and histone modifications) in ±32 kb of these ERVs' integration sites and in control regions, and analyzed them using Functional Data Analysis (FDA) methodology. In one of the first applications of FDA in genomics, we identified genomic scales and locations at which these features display their influence, and how they work in concert, to provide signals essential for integration and fixation of ERVs. The investigation of ERVs of different evolutionary ages (young in vitro and polymorphic ERVs, older fixed ERVs) allowed us to disentangle integration vs. fixation preferences. As a result of these analyses, we built a comprehensive model explaining the uneven distribution of ERVs along the genome. We found that ERVs integrate in late-replicating AT-rich regions with abundant microsatellites, mirror repeats, and repressive histone marks. Regions favoring fixation are depleted of genes and evolutionarily conserved elements, and have low recombination rates, reflecting the effects of purifying selection and ectopic recombination removing ERVs from the genome. In addition to providing these biological insights, our study demonstrates the power of exploiting multiple scales and localization with FDA. These powerful techniques are expected to be applicable to many other genomic investigations.

The dynamic genome: transposons and environmental adaptation in the nervous system

Classically thought as genomic clutter, the functional significance of transposable elements (TEs) has only recently become a focus of attention in neuroscience. Increasingly, studies have demonstrated that the brain seems to have more retrotransposition and TE transcription relative to other somatic tissues, suggesting a unique role for TEs in the central nervous system. TE expression and transposition also appear to vary by brain region and change in response to environmental stimuli such as stress. TEs appear to serve a number of adaptive roles in the nervous system. The regulation of TE expression by steroid, epigenetic and other mechanisms in interplay with the environment represents a significant and novel avenue to understanding both normal brain function and disease.

New players in heterochromatin silencing: histone variant H3.3 and the ATRX/DAXX chaperone

A number of studies have demonstrated that various components of the ATRX/DAXX/Histone H3.3 complex are important for heterochromatin silencing at multiple genomic regions. We provide an overview of the individual components (ATRX, DAXX and/or H3.3) tested in each study and propose a model where the ATRX/DAXX chaperone complex deposits H3.3 to maintain the H3K9me3 modification at heterochromatin throughout the genome.

Differential DNA Methylation Analysis without a Reference Genome

Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for DNA methylation analysis in non-model organisms, we developed an integrated approach for studying DNA methylation differences independent of a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, carp, sea bass, and zebrafish). Bioinformatically, we developed the RefFreeDMA software to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions between samples or groups of individuals (http://RefFreeDMA.computational-epigenetics.org). The identified regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell-type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.

Formation of Extrachromosomal Circular DNA from Long Terminal Repeats of Retrotransposons in Saccharomyces cerevisiae

Extrachromosomal circular DNA (eccDNA) derived from chromosomal Ty retrotransposons in yeast can be generated in multiple ways. Ty eccDNA can arise from the circularization of extrachromosomal linear DNA during the transpositional life cycle of retrotransposons, or from circularization of genomic Ty DNA. Circularization may happen through nonhomologous end-joining (NHEJ) of long terminal repeats (LTRs) flanking Ty elements, by Ty autointegration, or by LTR–LTR recombination. By performing an in-depth investigation of sequence reads stemming from Ty eccDNAs obtained from populations of Saccharomyces cerevisiae S288c, we find that eccDNAs predominantly correspond to full-length Ty1 elements. Analyses of sequence junctions reveal no signs of NHEJ or autointegration events. We detect recombination junctions that are consistent with yeast Ty eccDNAs being generated through recombination events within the genome. This opens the possibility that retrotransposable elements could move around in the genome without an RNA intermediate directly through DNA circularization.