Comparison of whole-genome bisulfite sequencing library preparation strategies identifies sources of biases affecting DNA methylation data.

BACKGROUND: Whole-genome bisulfite sequencing (WGBS) is becoming an increasingly accessible technique, used widely for both fundamental and disease-oriented research. Library preparation methods benefit from a variety of available kits, polymerases and bisulfite conversion protocols. Although some steps in the procedure, such as PCR amplification, are known to introduce biases, a systematic evaluation of biases in WGBS strategies is missing. RESULTS: We perform a comparative analysis of several commonly used pre- and post-bisulfite WGBS library preparation protocols for their performance and quality of sequencing outputs. Our results show that bisulfite conversion per se is the main trigger of pronounced sequencing biases, and PCR amplification builds on these underlying artefacts. The majority of standard library preparation methods yield a significantly biased sequence output and overestimate global methylation. Importantly, both absolute and relative methylation levels at specific genomic regions vary substantially between methods, with clear implications for DNA methylation studies. CONCLUSIONS: We show that amplification-free library preparation is the least biased approach for WGBS. In protocols with amplification, the choice of bisulfite conversion protocol or polymerase can significantly minimize artefacts. To aid with the quality assessment of existing WGBS datasets, we have integrated a bias diagnostic tool in the Bismark package and offer several approaches for consideration during the preparation and analysis of WGBS datasets

TEX15 is an essential executor of MIWI2-directed transposon DNA methylation and silencing.

The PIWI protein MIWI2 and its associated PIWI-interacting RNAs (piRNAs) instruct DNA methylation of young active transposable elements (TEs) in the male germline. piRNAs are proposed to recruit MIWI2 to the transcriptionally active TE loci by base pairing to nascent transcripts, however the downstream mechanisms and effector proteins utilized by MIWI2 in directing de novo TE methylation remain incompletely understood. Here, we show that MIWI2 associates with TEX15 in foetal gonocytes. TEX15 is predominantly a nuclear protein that is not required for piRNA biogenesis but is essential for piRNA-directed TE de novo methylation and silencing. In summary, TEX15 is an essential executor of mammalian piRNA-directed DNA methylation

MERVL/Zscan4 Network Activation Results in Transient Genome-wide DNA Demethylation of mESCs.

Mouse embryonic stem cells are dynamic and heterogeneous. For example, rare cells cycle through a state characterized by decondensed chromatin and expression of transcripts, including the Zscan4 cluster and MERVL endogenous retrovirus, which are usually restricted to preimplantation embryos. Here, we further characterize the dynamics and consequences of this transient cell state. Single-cell transcriptomics identified the earliest upregulated transcripts as cells enter the MERVL/Zscan4 state. The MERVL/Zscan4 transcriptional network was also upregulated during induced pluripotent stem cell reprogramming. Genome-wide DNA methylation and chromatin analyses revealed global DNA hypomethylation accompanying increased chromatin accessibility. This transient DNA demethylation was driven by a loss of DNA methyltransferase proteins in the cells and occurred genome-wide. While methylation levels were restored once cells exit this state, genomic imprints remained hypomethylated, demonstrating a potential global and enduring influence of endogenous retroviral activation on the epigenome

DNA methylation repels binding of hypoxia-inducible transcription factors to maintain tumor immunotolerance.

BACKGROUND: Hypoxia is pervasive in cancer and other diseases. Cells sense and adapt to hypoxia by activating hypoxia-inducible transcription factors (HIFs), but it is still an outstanding question why cell types differ in their transcriptional response to hypoxia. RESULTS: We report that HIFs fail to bind CpG dinucleotides that are methylated in their consensus binding sequence, both in in vitro biochemical binding assays and in vivo studies of differentially methylated isogenic cell lines. Based on in silico structural modeling, we show that 5-methylcytosine indeed causes steric hindrance in the HIF binding pocket. A model wherein cell-type-specific methylation landscapes, as laid down by the differential expression and binding of other transcription factors under normoxia, control cell-type-specific hypoxia responses is observed. We also discover ectopic HIF binding sites in repeat regions which are normally methylated. Genetic and pharmacological DNA demethylation, but also cancer-associated DNA hypomethylation, expose these binding sites, inducing HIF-dependent expression of cryptic transcripts. In line with such cryptic transcripts being more prone to cause double-stranded RNA and viral mimicry, we observe low DNA methylation and high cryptic transcript expression in tumors with high immune checkpoint expression, but not in tumors with low immune checkpoint expression, where they would compromise tumor immunotolerance. In a low-immunogenic tumor model, DNA demethylation upregulates cryptic transcript expression in a HIF-dependent manner, causing immune activation and reducing tumor growth. CONCLUSIONS: Our data elucidate the mechanism underlying cell-type-specific responses to hypoxia and suggest DNA methylation and hypoxia to underlie tumor immunotolerance

A MILI-independent piRNA biogenesis pathway empowers partial germline reprogramming.

In mice, the pathway involving PIWI and PIWI-interacting RNA (PIWI-piRNA) is essential to re-establish transposon silencing during male-germline reprogramming. The cytoplasmic PIWI protein MILI mediates piRNA-guided transposon RNA cleavage as well as piRNA amplification. MIWI2's binding to piRNA and its nuclear localization are proposed to be dependent upon MILI function. Here, we demonstrate the existence of a piRNA biogenesis pathway that sustains partial MIWI2 function and reprogramming activity in the absence of MILI

Defective germline reprogramming rewires the spermatogonial transcriptome.

Defective germline reprogramming in Piwil4 (Miwi2)- and Dnmt3l-deficient mice results in the failure to reestablish transposon silencing, meiotic arrest and progressive loss of spermatogonia. Here we sought to understand the molecular basis for this spermatogonial dysfunction. Through a combination of imaging, conditional genetics and transcriptome analysis, we demonstrate that germ cell elimination in the respective mutants arises as a result of defective de novo genome methylation during reprogramming rather than because of a function for the respective factors within spermatogonia. In both Miwi2-/- and Dnmt3l-/- spermatogonia, the intracisternal-A particle (IAP) family of endogenous retroviruses is derepressed, but, in contrast to meiotic cells, DNA damage is not observed. Instead, we find that unmethylated IAP promoters rewire the spermatogonial transcriptome by driving expression of neighboring genes. Finally, spermatogonial numbers, proliferation and differentiation are altered in Miwi2-/- and Dnmt3l-/- mice. In summary, defective reprogramming deregulates the spermatogonial transcriptome and may underlie spermatogonial dysfunction

An endosiRNA-Based Repression Mechanism Counteracts Transposon Activation during Global DNA Demethylation in Embryonic Stem Cells

Erasure of DNA methylation and repressive chromatin marks in the mammalian germline leads to risk of transcriptional activation of transposable elements (TEs). Here, we used mouse embryonic stem cells (ESCs) to identify an endosiRNA-based mechanism involved in suppression of TE transcription. In ESCs with DNA demethylation induced by acute deletion of Dnmt1, we saw an increase in sense transcription at TEs, resulting in an abundance of sense/antisense transcripts leading to high levels of ARGONAUTE2 (AGO2)-bound small RNAs. Inhibition of Dicer or Ago2 expression revealed that small RNAs are involved in an immediate response to demethylation-induced transposon activation, while the deposition of repressive histone marks follows as a chronic response. In vivo, we also found TE-specific endosiRNAs present during primordial germ cell development. Our results suggest that antisense TE transcription is a “trap” that elicits an endosiRNA response to restrain acute transposon activity during epigenetic reprogramming in the mammalian germline.ISSN:1934-5909ISSN:1875-977

Role of small RNAs and chromatin in transposable element silencing during global demethylation

DNA methylation entails the addition of a methyl group to the 5-carbon of the cytosine base of the DNA. This modification is important during many biological processes such as imprinting, X-chromosome inactivation, cell differentiation as well as silencing of transposable elements (TEs). DNA methylation is dynamic during early mammalian development, despite being a more static mark in somatic cells. Global hypomethylation is a hallmark of epigenetic reprogramming in mammalian primordial germ cells (PGCs), the early embryo and in naïve embryonic stem cells (ESCs). Genome integrity is crucial during early development, as the germline DNA needs to be protected for future generations. Therefore, epigenetic reprogramming presents a critical phase for TE defence since presumably alternative silencing pathways need to be employed to limit their activity. In this thesis, I investigate the role of small RNAs to control TEs during global waves of DNA demethylation in cellular reprogramming, naïve pluripotency as well as early mammalian development. Following an introduction to the research questions, in chapter 3 I investigate the mechanism of TE regulation in an in vitro model of Dnmt1 deletion in mouse ES cells to recapitulate in vivo epigenetic reprogramming. I find that certain classes of TEs become transcriptionally upregulated and subsequently resilenced by a mechanism independent of DNA methylation. I identify ARGONAUTE 2 (AGO2) bound siRNAs as the prominent mechanism to control certain classes of TEs, while others appear to be regulated by redistribution of repressive histone modifications. In chapter 4, I construct Dicer constitutive and conditional KO ESCs in the background of the Dnmt1f l/f l ESCs using CRISPR-Cas9. I dissect the role of DNA methylation and of DICER dependent small RNAs on transcriptional changes of ESCs. Additionally, I find that DICER dependent small interfering RNAs (siRNAs) re-silence transcriptionally active TE classes. Finally, in chapter 5, I examine the role of small RNAs in TE silencing in different models of global hypomethylation in vivo and in vitro PGCs, during iPSC reprogramming and in a transition from serum to 2i culturing of mouse ESCs.Gates Cambridg

Transgenerational transmission of hedonic behaviors and metabolic phenotypes induced by maternal overnutrition

Maternal overnutrition has been associated with increased susceptibility to develop obesity and neurological disorders later in life. Most epidemiological as well as experimental studies have focused on the metabolic consequences across generations following an early developmental nutritional insult. Recently, it has been shown that maternal high-fat diet (HFD) affects third-generation female body mass via the paternal lineage. We showed here that the offspring born to HFD ancestors displayed addictive-like behaviors as well as obesity and insulin resistance up to the third generation in the absence of any further exposure to HFD. These findings, implicate that the male germ line is a major player in transferring phenotypic traits. These behavioral and physiological alterations were paralleled by reduced striatal dopamine levels and increased dopamine 2 receptor density. Interestingly, by the third generation a clear gender segregation emerged, where females showed addictive-like behaviors while male HFD offspring showed an obesogenic phenotype. However, methylome profiling of F1 and F2 sperm revealed no significant difference between the offspring groups, suggesting that the sperm methylome might not be the major carrier for the transmission of the phenotypes observed in our mouse model. Together, our study for the first time demonstrates that maternal HFD insult causes sustained alterations of the mesolimbic dopaminergic system suggestive of a predisposition to develop obesity and addictive-like behaviors across multiple generations.ISSN:2158-318

Comparison of whole-genome bisulfite sequencing library preparation strategies identifies sources of biases affecting DNA methylation data.

BACKGROUND: Whole-genome bisulfite sequencing (WGBS) is becoming an increasingly accessible technique, used widely for both fundamental and disease-oriented research. Library preparation methods benefit from a variety of available kits, polymerases and bisulfite conversion protocols. Although some steps in the procedure, such as PCR amplification, are known to introduce biases, a systematic evaluation of biases in WGBS strategies is missing. RESULTS: We perform a comparative analysis of several commonly used pre- and post-bisulfite WGBS library preparation protocols for their performance and quality of sequencing outputs. Our results show that bisulfite conversion per se is the main trigger of pronounced sequencing biases, and PCR amplification builds on these underlying artefacts. The majority of standard library preparation methods yield a significantly biased sequence output and overestimate global methylation. Importantly, both absolute and relative methylation levels at specific genomic regions vary substantially between methods, with clear implications for DNA methylation studies. CONCLUSIONS: We show that amplification-free library preparation is the least biased approach for WGBS. In protocols with amplification, the choice of bisulfite conversion protocol or polymerase can significantly minimize artefacts. To aid with the quality assessment of existing WGBS datasets, we have integrated a bias diagnostic tool in the Bismark package and offer several approaches for consideration during the preparation and analysis of WGBS datasets