A co-ordinated transcriptional programme in the maternal liver supplies long chain polyunsaturated fatty acids to the conceptus using phospholipids.

The long and very long chain polyunsaturated fatty acids (LC-PUFAs) are preferentially transported by the mother to the fetus. Failure to supply LC-PUFAs is strongly linked with stillbirth, fetal growth restriction, and impaired neurodevelopmental outcomes. However, dietary supplementation during pregnancy is unable to simply reverse these outcomes, suggesting imperfectly understood interactions between dietary fatty acid intake and the molecular mechanisms of maternal supply. Here we employ a comprehensive approach combining untargeted and targeted lipidomics with transcriptional profiling of maternal and fetal tissues in mouse pregnancy. Comparison of wild-type mice with genetic models of impaired lipid metabolism allows us to describe maternal hepatic adaptations required to provide LC-PUFAs to the developing fetus. A late pregnancy-specific, selective activation of the Liver X Receptor signalling pathway dramatically increases maternal supply of LC-PUFAs within circulating phospholipids. Crucially, genetic ablation of this pathway in the mother reduces LC-PUFA accumulation by the fetus, specifically of docosahexaenoic acid (DHA), a critical nutrient for brain development.

Modeling methyl-sensitive transcription factor motifs with an expanded epigenetic alphabet.

BACKGROUND: Transcription factors bind DNA in specific sequence contexts. In addition to distinguishing one nucleobase from another, some transcription factors can distinguish between unmodified and modified bases. Current models of transcription factor binding tend not to take DNA modifications into account, while the recent few that do often have limitations. This makes a comprehensive and accurate profiling of transcription factor affinities difficult. RESULTS: Here, we develop methods to identify transcription factor binding sites in modified DNA. Our models expand the standard A/C/G/T DNA alphabet to include cytosine modifications. We develop Cytomod to create modified genomic sequences and we also enhance the MEME Suite, adding the capacity to handle custom alphabets. We adapt the well-established position weight matrix (PWM) model of transcription factor binding affinity to this expanded DNA alphabet. Using these methods, we identify modification-sensitive transcription factor binding motifs. We confirm established binding preferences, such as the preference of ZFP57 and C/EBPß for methylated motifs and the preference of c-Myc for unmethylated E-box motifs. CONCLUSIONS: Using known binding preferences to tune model parameters, we discover novel modified motifs for a wide array of transcription factors. Finally, we validate our binding preference predictions for OCT4 using cleavage under targets and release using nuclease (CUT&RUN) experiments across conventional, methylation-, and hydroxymethylation-enriched sequences. Our approach readily extends to other DNA modifications. As more genome-wide single-base resolution modification data becomes available, we expect that our method will yield insights into altered transcription factor binding affinities across many different modifications.

Epigenetic control and genomic imprinting dynamics of the Dlk1-Dio3 domain.

Genomic imprinting is an epigenetic process whereby genes are monoallelically expressed in a parent-of-origin-specific manner. Imprinted genes are frequently found clustered in the genome, likely illustrating their need for both shared regulatory control and functional inter-dependence. The Dlk1-Dio3 domain is one of the largest imprinted clusters. Genes in this region are involved in development, behavior, and postnatal metabolism: failure to correctly regulate the domain leads to Kagami-Ogata or Temple syndromes in humans. The region contains many of the hallmarks of other imprinted domains, such as long non-coding RNAs and parental origin-specific CTCF binding. Recent studies have shown that the Dlk1-Dio3 domain is exquisitely regulated via a bipartite imprinting control region (ICR) which functions differently on the two parental chromosomes to establish monoallelic expression. Furthermore, the Dlk1 gene displays a selective absence of imprinting in the neurogenic niche, illustrating the need for precise dosage modulation of this domain in different tissues. Here, we discuss the following: how differential epigenetic marks laid down in the gametes cause a cascade of events that leads to imprinting in the region, how this mechanism is selectively switched off in the neurogenic niche, and why studying this imprinted region has added a layer of sophistication to how we think about the hierarchical epigenetic control of genome function.

Epigenetic inheritance is unfaithful at intermediately methylated CpG sites.

DNA methylation at the CpG dinucleotide is considered a stable epigenetic mark due to its presumed long-term inheritance through clonal expansion. Here, we perform high-throughput bisulfite sequencing on clonally derived somatic cell lines to quantitatively measure methylation inheritance at the nucleotide level. We find that although DNA methylation is generally faithfully maintained at hypo- and hypermethylated sites, this is not the case at intermediately methylated CpGs. Low fidelity intermediate methylation is interspersed throughout the genome and within genes with no or low transcriptional activity, and is not coordinately maintained between neighbouring sites. We determine that the probabilistic changes that occur at intermediately methylated sites are likely due to DNMT1 rather than DNMT3A/3B activity. The observed lack of clonal inheritance at intermediately methylated sites challenges the current epigenetic inheritance model and has direct implications for both the functional relevance and general interpretability of DNA methylation as a stable epigenetic mark.

Mammary adipocyte flow cytometry as a tool to study mammary gland biology.

The mammary gland is a vital exocrine organ that has evolved in mammals to secrete milk and provide nutrition to ensure the growth and survival of the neonate The mouse mammary gland displays extraordinary plasticity each time the female undergoes pregnancy and lactation, including a sophisticated process of tertiary branching and alveologenesis to form a branched epithelial tree and subsequently milk-producing alveoli. Upon the cessation of lactation, the gland remodels back to a simple ductal architecture via highly regulated involution processes. At the cellular level, the plasticity is characterised by proliferation of mammary cell populations, differentiation and apoptosis, accompanied by major changes in cell function and morphology. The mammary epithelium requires a specific stromal environment to grow, known as the mammary fat pad. Mammary adipocytes are one of the most prominent cell types in the fat pad, but despite their vast proportion in the tissue and their crucial interaction with epithelial cells, their physiology remains largely unknown. Over the past decade, the need to understand the properties and contribution of mammary adipocytes has become more recognised. However, the development of adequate methods and protocols to study this cellular niche is still lagging, partially due to their fragile nature, the difficulty of isolating them, the lack of reliable cell surface markers and the heterogenous environment in this tissue, which differs from other adipocyte depots. Here, we describe a new rapid and simple flow cytometry protocol specifically designed for the analysis and isolation of mouse mammary adipocytes across mammary gland developmental stages.

Reassessment of weak parent-of-origin expression bias shows it rarely exists outside of known imprinted regions.

In mouse and human, genes subjected to genomic imprinting have been shown to function in development, behavior, and post-natal adaptations. Failure to correctly imprint genes in human is associated with developmental syndromes, adaptive, and metabolic disorders during life as well as numerous forms of cancer. In recent years researchers have turned to RNA-seq technologies applied to reciprocal hybrid strains of mice to identify novel imprinted genes, causing a threefold increase in genes reported as having a parental origin-specific expression bias. The functional relevance of parental origin-specific expression bias is not fully appreciated especially since many are reported with only minimal parental bias (e.g. 51:49). Here, we present an in-depth meta-analysis of previously generated RNA-seq data and show that the methods used to generate and analyze libraries greatly influence the calling of allele-specific expression. Validation experiments show that most novel genes called with parental-origin-specific allelic bias are artefactual, with the mouse strain contributing a larger effect on expression biases than parental origin. Of the weak novel genes that do validate, most are located at the periphery of known imprinted domains, suggesting they may be affected by local allele- and tissue-specific conformation. Together these findings highlight the need for robust tools, definitions, and validation of putative imprinted genes to provide meaningful information within imprinting databases and to understand the functional and mechanistic implications of the process.

Balanced gene dosage control rather than parental origin underpins genomic imprinting.

Mammalian parental imprinting represents an exquisite form of epigenetic control regulating the parent-specific monoallelic expression of genes in clusters. While imprinting perturbations are widely associated with developmental abnormalities, the intricate regional interplay between imprinted genes makes interpreting the contribution of gene dosage effects to phenotypes a challenging task. Using mouse models with distinct deletions in an intergenic region controlling imprinting across the Dlk1-Dio3 domain, we link changes in genetic and epigenetic states to allelic-expression and phenotypic outcome in vivo. This determined how hierarchical interactions between regulatory elements orchestrate robust parent-specific expression, with implications for non-imprinted gene regulation. Strikingly, flipping imprinting on the parental chromosomes by crossing genotypes of complete and partial intergenic element deletions rescues the lethality of each deletion on its own. Our work indicates that parental origin of an epigenetic state is irrelevant as long as appropriate balanced gene expression is established and maintained at imprinted loci.

Epigenetic changes induced by in utero dietary challenge result in phenotypic variability in successive generations of mice.

Transmission of epigenetic information between generations occurs in nematodes, flies and plants, mediated by specialised small RNA pathways, modified histones and DNA methylation. Similar processes in mammals can also affect phenotype through intergenerational or trans-generational mechanisms. Here we generate a luciferase knock-in reporter mouse for the imprinted Dlk1 locus to visualise and track epigenetic fidelity across generations. Exposure to high-fat diet in pregnancy provokes sustained re-expression of the normally silent maternal Dlk1 in offspring (loss of imprinting) and increased DNA methylation at the somatic differentially methylated region (sDMR). In the next generation heterogeneous Dlk1 mis-expression is seen exclusively among animals born to F1-exposed females. Oocytes from these females show altered gene and microRNA expression without changes in DNA methylation, and correct imprinting is restored in subsequent generations. Our results illustrate how diet impacts the foetal epigenome, disturbing canonical and non-canonical imprinting mechanisms to modulate the properties of successive generations of offspring.

Subnuclear localisation is associated with gene expression more than parental origin at the imprinted Dlk1-Dio3 locus.

At interphase, de-condensed chromosomes have a non-random three-dimensional architecture within the nucleus, however, little is known about the extent to which nuclear organisation might influence expression or vice versa. Here, using imprinting as a model, we use 3D RNA- and DNA-fluorescence-in-situ-hybridisation in normal and mutant mouse embryonic stem cell lines to assess the relationship between imprinting control, gene expression and allelic distance from the nuclear periphery. We compared the two parentally inherited imprinted domains at the Dlk1-Dio3 domain and find a small but reproducible trend for the maternally inherited domain to be further away from the periphery however we did not observe an enrichment of inactive alleles in the immediate vicinity of the nuclear envelope. Using Zfp57KO ES cells, which harbour a paternal to maternal epigenotype switch, we observe that expressed alleles are significantly further away from the nuclear periphery. However, within individual nuclei, alleles closer to the periphery are equally likely to be expressed as those further away. In other words, absolute position does not predict expression. Taken together, this suggests that whilst stochastic activation can cause subtle shifts in localisation for this locus, there is no dramatic relocation of alleles upon gene activation. Our results suggest that transcriptional activity, rather than the parent-of-origin, defines subnuclear localisation at an endogenous imprinted domain.

Preimplantation genetic testing for a chr14q32 microdeletion in a family with Kagami-Ogata syndrome and Temple syndrome.

INTRODUCTION: Kagami-Ogata syndrome (KOS14) and Temple syndrome (TS14) are two disorders associated with reciprocal alterations within the chr14q32 imprinted domain. Here, we present a work-up strategy for preimplantation genetic testing (PGT) to avoid the transmission of a causative micro-deletion. METHODS: We analysed DNA from the KOS14 index case and parents using methylation-sensitive ligation-mediated probe amplification and methylation pyrosequencing. The extent of the deletion was mapped using SNP arrays. PGT was performed in trophectoderm samples in order to identify unaffected embryos. Samples were amplified using multiple displacement amplification, followed by genome-wide SNP genotyping to determine the at-risk haplotype and next-generation sequencing to determine aneuploidies. RESULTS: A fully methylated pattern at the normally paternally methylated IG-DMR and MEG3 DMR in the KOS14 proband, accompanied by an unmethylated profile in the TS14 mother was consistent with maternal and paternal transmission of a deletion, respectively. Further analysis revealed a 108 kb deletion in both cases. The inheritance of the deletion on different parental alleles was consistent with the opposing phenotypes. In vitro fertilisation with intracytoplasmatic sperm injection and PGT were used to screen for deletion status and to transfer an unaffected embryo in this couple. A single euploid-unaffected embryo was identified resulting in a healthy baby born. DISCUSSION: We identify a microdeletion responsible for multigeneration KOS14 and TS14 within a single family where carriers have a 50% risk of transmitting the deletion to their offspring. We show that PGT can successfully be offered to couples with IDs caused by genetic anomalies.

Epigenetic Mechanisms of ART-Related Imprinting Disorders: Lessons From iPSC and Mouse Models.

The rising frequency of ART-conceived births is accompanied by the need for an improved understanding of the implications of ART on gametes and embryos. Increasing evidence from mouse models and human epidemiological data suggests that ART procedures may play a role in the pathophysiology of certain imprinting disorders (IDs), including Beckwith-Wiedemann syndrome, Silver-Russell syndrome, Prader-Willi syndrome, and Angelman syndrome. The underlying molecular basis of this association, however, requires further elucidation. In this review, we discuss the epigenetic and imprinting alterations of in vivo mouse models and human iPSC models of ART. Mouse models have demonstrated aberrant regulation of imprinted genes involved with ART-related IDs. In the past decade, iPSC technology has provided a platform for patient-specific cellular models of culture-associated perturbed imprinting. However, despite ongoing efforts, a deeper understanding of the susceptibility of iPSCs to epigenetic perturbation is required if they are to be reliably used for modelling ART-associated IDs. Comparing the patterns of susceptibility of imprinted genes in mouse models and IPSCs in culture improves the current understanding of the underlying mechanisms of ART-linked IDs with implications for our understanding of the influence of environmental factors such as culture and hormone treatments on epigenetically important regions of the genome such as imprints.

Variably methylated retrotransposons are refractory to a range of environmental perturbations.

The agouti viable yellow (Avy) allele is an insertional mutation in the mouse genome caused by a variably methylated intracisternal A particle (VM-IAP) retrotransposon. Avy expressivity is sensitive to a range of early-life chemical exposures and nutritional interventions, suggesting that environmental perturbations can have long-lasting effects on the methylome. However, the extent to which VM-IAP elements are environmentally labile with phenotypic implications is unknown. Using a recently identified repertoire of VM-IAPs, we assessed the epigenetic effects of different environmental contexts. A longitudinal aging analysis indicated that VM-IAPs are stable across the murine lifespan, with only small increases in DNA methylation detected for a subset of loci. No significant effects were observed after maternal exposure to the endocrine disruptor bisphenol A, an obesogenic diet or methyl donor supplementation. A genetic mouse model of abnormal folate metabolism exhibited shifted VM-IAP methylation levels and altered VM-IAP-associated gene expression, yet these effects are likely largely driven by differential targeting by polymorphic KRAB zinc finger proteins. We conclude that epigenetic variability at retrotransposons is not predictive of environmental susceptibility.

Defective folate metabolism causes germline epigenetic instability and distinguishes Hira as a phenotype inheritance biomarker.

The mechanism behind transgenerational epigenetic inheritance is unclear, particularly through the maternal grandparental line. We previously showed that disruption of folate metabolism in mice by the Mtrr hypomorphic mutation results in transgenerational epigenetic inheritance of congenital malformations. Either maternal grandparent can initiate this phenomenon, which persists for at least four wildtype generations. Here, we use genome-wide approaches to reveal genetic stability in the Mtrr model and genome-wide differential DNA methylation in the germline of Mtrr mutant maternal grandfathers. We observe that, while epigenetic reprogramming occurs, wildtype grandprogeny and great grandprogeny exhibit transcriptional changes that correlate with germline methylation defects. One region encompasses the Hira gene, which is misexpressed in embryos for at least three wildtype generations in a manner that distinguishes Hira transcript expression as a biomarker of maternal phenotypic inheritance.

Dlk1 dosage regulates hippocampal neurogenesis and cognition.

Neurogenesis in the adult brain gives rise to functional neurons, which integrate into neuronal circuits and modulate neural plasticity. Sustained neurogenesis throughout life occurs in the subgranular zone (SGZ) of the dentate gyrus in the hippocampus and is hypothesized to be involved in behavioral/cognitive processes such as memory and in diseases. Genomic imprinting is of critical importance to brain development and normal behavior, and exemplifies how epigenetic states regulate genome function and gene dosage. While most genes are expressed from both alleles, imprinted genes are usually expressed from either the maternally or the paternally inherited chromosome. Here, we show that in contrast to its canonical imprinting in nonneurogenic regions, Delta-like homolog 1 (Dlk1) is expressed biallelically in the SGZ, and both parental alleles are required for stem cell behavior and normal adult neurogenesis in the hippocampus. To evaluate the effects of maternally, paternally, and biallelically inherited mutations within the Dlk1 gene in specific behavioral domains, we subjected Dlk1-mutant mice to a battery of tests that dissociate and evaluate the effects of Dlk1 dosage on spatial learning ability and on anxiety traits. Importantly, reduction in Dlk1 levels triggers specific cognitive abnormalities that affect aspects of discriminating differences in environmental stimuli, emphasizing the importance of selective absence of imprinting in this neurogenic niche.

A spontaneous genetically induced epiallele at a retrotransposon shapes host genome function.

Intracisternal A-particles (IAPs) are endogenous retroviruses (ERVs) responsible for most insertional mutations in the mouse. Full-length IAPs harbour genes flanked by long terminal repeats (LTRs). Here, we identify a solo LTR IAP variant (Iap5-1solo) recently formed in the inbred C57BL/6J mouse strain. In contrast to the C57BL/6J full-length IAP at this locus (Iap5-1full), Iap5-1solo lacks DNA methylation and H3K9 trimethylation. The distinct DNA methylation levels between the two alleles are established during preimplantation development, likely due to loss of KRAB zinc finger protein binding at the Iap5-1solo variant. Iap5-1solo methylation increases and becomes more variable in a hybrid genetic background yet is unresponsive to maternal dietary methyl supplementation. Differential epigenetic modification of the two variants is associated with metabolic differences and tissue-specific changes in adjacent gene expression. Our characterisation of Iap5-1 as a genetically induced epiallele with functional consequences establishes a new model to study transposable element repression and host-element co-evolution.

Our genome provides a complete set of genetic instructions for life. It begins by directing the growth and development of the embryo, and subsequently supports all the cells of the adult body in their daily routines. Yet approximately 10% of the DNA in mammalian genomes is made up of sequences originating from past retroviral infections, leaving a calling card in our genetic code. While these segments of retroviral DNA can no longer produce new infectious viruses, some of them retain the ability to copy themselves and jump into new parts of the genome. This can be problematic if they jump into and disrupt an important piece of genetic code. To protect against this, our bodies have evolved the ability to chemically strap down retroviral sequences by adding methyl groups to them and by modifying the proteins they are wrapped around. However, some of these endogenous retroviruses can dodge such so-called epigenetic modifications and disrupt genome function as a result. Studying a population of widely used inbred laboratory mice, Bertozzi et al. have identified a retroviral element that evades these epigenetic restraints. They discovered that some mice carry a full-length retroviral sequence while others have a shortened version of the same element. The shorter sequence lacked the repressive epigenetic marks found on the longer version, and this affected the expression of nearby genes. Moreover, the repressive marks could be partially restored by breeding the short-version mice with a distantly related mouse strain. Bertozzi et al. highlight an important issue for research using mouse models. Inbred laboratory mouse strains are assumed to have a fixed genetic code which allows scientists to conclude that any observed differences in their experiments are not a product of background genetic variation. However, this study emphasizes that this assumption is not guaranteed, and that hidden genetic diversity may be present in ostensibly genetically identical mice, with important implications for experimental outcomes. In addition, Bertozzi et al. provide a new mouse model for researchers to study the evolution and regulation of retroviral sequences and the impact of these processes on cell function.

Genomic properties of variably methylated retrotransposons in mouse.

BACKGROUND: Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function. RESULTS: Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation. CONCLUSION: Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

Imprinting methylation predicts hippocampal volumes and hyperintensities and the change with age in later life.

Epigenetic imprinting is important for neurogenesis and brain function. Hippocampal volumes and brain hyperintensities in late life have been associated with early life circumstances. Epigenetic imprinting may underpin these associations. Methylation was measured at 982 sites in 13 imprinted locations in blood samples from a longitudinal cohort by bisulphite amplicon sequencing. Hippocampal volumes and hyperintensities were determined at age 64y and 72y using MRI. Hyperintensities were determined in white matter, grey matter and infratentorial regions. Permutation methods were used to adjust for multiple testing. At 64y, H19/IGF2 and NESPAS methylation predicted hippocampal volumes. PEG3 predicted hyperintensities in hippocampal grey matter, and white matter. GNASXL predicted grey matter hyperintensities. Changes with age were predicted for hippocampal volume (MEST1, KvDMR, L3MBTL, GNASXL), white matter (MEST1, PEG3) and hippocampal grey matter hyperintensities (MCTS2, GNASXL, NESPAS, L3MBTL, MCTS2, SNRPN, MEST1). Including childhood cognitive ability, years in education, or socioeconomic status as additional explanatory variables in regression analyses did not change the overall findings. Imprinting methylation in multiple genes predicts brain structures, and their change over time. These findings are potentially relevant to the development of novel tests of brain structure and function across the life-course, strategies to improve cognitive outcomes, and our understanding of early influences on brain development and function.

KRAB zinc finger protein diversification drives mammalian interindividual methylation variability.

Most transposable elements (TEs) in the mouse genome are heavily modified by DNA methylation and repressive histone modifications. However, a subset of TEs exhibit variable methylation levels in genetically identical individuals, and this is associated with epigenetically conferred phenotypic differences, environmental adaptability, and transgenerational epigenetic inheritance. The evolutionary origins and molecular mechanisms underlying interindividual epigenetic variability remain unknown. Using a repertoire of murine variably methylated intracisternal A-particle (VM-IAP) epialleles as a model, we demonstrate that variable DNA methylation states at TEs are highly susceptible to genetic background effects. Taking a classical genetics approach coupled with genome-wide analysis, we harness these effects and identify a cluster of KRAB zinc finger protein (KZFP) genes that modifies VM-IAPs in trans in a sequence-specific manner. Deletion of the cluster results in decreased DNA methylation levels and altered histone modifications at the targeted VM-IAPs. In some cases, these effects are accompanied by dysregulation of neighboring genes. We find that VM-IAPs cluster together phylogenetically and that this is linked to differential KZFP binding, suggestive of an ongoing evolutionary arms race between TEs and this large family of epigenetic regulators. These findings indicate that KZFP divergence and concomitant evolution of DNA binding capabilities are mechanistically linked to methylation variability in mammals, with implications for phenotypic variation and putative paradigms of mammalian epigenetic inheritance.