Background Repeat elements constitute a large proportion of the human genome and recent evidence indicates that repeat element expression has functional roles in both physiological and pathological states. Specifically for cancer, transcription of endogenous retrotransposons is often suppressed to attenuate an anti-tumor immune response, whereas aberrant expression of heterochromatin-derived satellite RNA has been identified as a tumor driver. These insights demonstrate separate functions for the dysregulation of distinct repeat subclasses in either the attenuation or progression of human solid tumors. For hematopoietic malignancies, such as Acute Myeloid Leukemia (AML), only very few studies on the expression/dysregulation of repeat elements were done. Methods To study the expression of repeat elements in AML, we performed total-RNA sequencing of healthy CD34 + cells and of leukemic blast cells from primary AML patient material. We also developed an integrative bioinformatic approach that can quantify the expression of repeat transcripts from all repeat subclasses (SINE/ALU, LINE, ERV and satellites) in relation to the expression of gene and other non-repeat transcripts (i.e. R/G ratio). This novel approach can be used as an instructive signature for repeat element expression and has been extended to the analysis of poly(A)-RNA sequencing datasets from Blueprint and TCGA consortia that together comprise 120 AML patient samples. Results We identified that repeat element expression is generally down-regulated during hematopoietic differentiation and that relative changes in repeat to gene expression can stratify risk prediction of AML patients and correlate with overall survival probabilities. A high R/G ratio identifies AML patient subgroups with a favorable prognosis, whereas a low R/G ratio is prevalent in AML patient subgroups with a poor prognosis. Conclusions We developed an integrative bioinformatic approach that defines a general model for the analysis of repeat element dysregulation in physiological and pathological development. We find that changes in repeat to gene expression (i.e. R/G ratios) correlate with hematopoietic differentiation and can sub-stratify AML patients into low-risk and high-risk subgroups. Thus, the definition of a R/G ratio can serve as a valuable biomarker for AML and could also provide insights into differential patient response to epigenetic drug treatment.
Histone H3 lysine 9 (H3K9) methylation is a central epigenetic modification that defines heterochromatin from unicellular to multicellular organisms. In mammalian cells, H3K9 methylation can be catalyzed by at least six distinct SET domain enzymes: Suv39h1/Suv39h2, Eset1/Eset2 and G9a/Glp. We used mouse embryonic fibroblasts (MEFs) with a conditional mutation for Eset1 and introduced progressive deletions for the other SET domain genes by CRISPR/Cas9 technology. Compound mutant MEFs for all six SET domain lysine methyltransferase (KMT) genes lack all H3K9 methylation states, derepress nearly all families of repeat elements and display genomic instabilities. Strikingly, the 6KO H3K9 KMT MEF cells no longer maintain heterochromatin organization and have lost electron-dense heterochromatin. This is a compelling analysis of H3K9 methylation-deficient mammalian chromatin and reveals a definitive function for H3K9 methylation in protecting heterochromatin organization and genome integrity.
DNA methylation is an important component in vital biological functions such as embryonic development, carcinogenesis, and heritable regulation. Accurate methods to assess genomic methylation status are crucial to its effective use in many scenarios, especially in the detection and diagnosis of disease. Methylation aligners, such as Bismark and bwa-meth, frequently assign significantly higher MapQ values than can be supported by the uniqueness of the region reads are mapped to. These incorrectly high MapQs result in inappropriate methylation calling in repetitive regions. We observe reads that should map to separate locations (possibly having different methylation states) actually end up mapping to the same locus, causing apparent mixed methylation at such loci. Methylation calling can be improved by using Bismap mappability data to filter out insufficiently unique reads. However, simply filtering out Cs in insufficiently unique regions is not adequate as it is prone to over-filtering Cs in small mappability dips. These Cs can in fact often be called using reads anchored in a nearby mappable region. We have created a new feature for the MethylDackel methylation caller to perform read-based filtering. This new methylation calling method resolves some of the apparent mixed methylation to either 0% or 100% methylation and removes many unsupportable methylation calls. We examined methylation calls with and without read-based filtering in or near the 7830 genes containing ClinVar variants in a methylation sequencing data set from the NA12878 cell line. Use of this improved method corrected 41,143 mixed methylation Cs to 0% methylation, and 22,345 to 100% methylation throughout the genome.
Repeat element transcription plays a vital role in early embryonic development. Expression of repeats such as MERVL characterises mouse embryos at the 2-cell stage, and defines a 2-cell-like cell (2CLC) population in a mouse embryonic stem cell culture. Repeat element sequences contain binding sites for numerous transcription factors. We identify the forkhead domain transcription factor FOXD3 as a regulator of repeat element transcription in mouse embryonic stem cells. FOXD3 binds to and recruits the histone methyltransferase SUV39H1 to MERVL and major satellite repeats, consequentially repressing the transcription of these repeats by the establishment of the H3K9me3 heterochromatin modification. Notably, depletion of FOXD3 leads to the de-repression of MERVL and major satellite repeats as well as a subset of genes expressed in the 2-cell state, shifting the balance between the stem cell and 2-cell like population in culture. Thus, FOXD3 acts as a negative regulator of repeat transcription, ascribing a novel function to this transcription factor.
Repeat element transcription plays a vital role in early embryonic development. The expression of repeats such as MERVL characterises mouse embryos at the 2-cell stage and defines a 2-cell-like cell (2CLC) population in a mouse embryonic stem cell culture. Repeat element sequences contain binding sites for numerous transcription factors. We identify the forkhead domain transcription factor FOXD3 as a regulator of major satellite repeats and MERVL transcription in mouse embryonic stem cells. FOXD3 binds to and recruits the histone methyltransferase SUV39H1 to MERVL and major satellite repeats, consequentially repressing the transcription of these repeats by the establishment of the H3K9me3 heterochromatin modification. Notably, depletion of FOXD3 leads to the de-repression of MERVL and major satellite repeats as well as a subset of genes expressed in the 2-cell state, shifting the balance between the stem cell and 2-cell-like population in culture. Thus, FOXD3 acts as a negative regulator of repeat transcription, ascribing a novel function to this transcription factor.
SUMMARY:Due to the rapidly increasing scale and diversity of epigenomic data, modular and scalable analysis workflows are of wide interest. Here we present snakePipes, a workflow package for processing and downstream analysis of data from common epigenomic assays: ChIP-seq, RNA-seq, Bisulfite-seq, ATAC-seq, Hi-C and single-cell RNA-seq. snakePipes enables users to assemble variants of each workflow and to easily install and upgrade the underlying tools, via its simple command-line wrappers and yaml files.AVAILABILITY AND IMPLEMENTATION:snakePipes can be installed via conda: `conda install -c mpi-ie -c bioconda -c conda-forge snakePipes'. Source code (https://github.com/maxplanck-ie/snakepipes) and documentation (https://snakepipes.readthedocs.io/en/latest/) are available online.SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
MOTIVATION:This paper presents Parkour, a software package for sample processing and quality management of next generation sequencing data and samples.RESULTS:Starting with user requests, Parkour allows tracking and assessing samples based on predefined quality criteria through different stages of the sample preparation workflow. Ideally suited for academic core laboratories, the software aims to maximize efficiency and reduce turnaround time by intelligent sample grouping and a clear assignment of staff to work units. Tools for automated invoicing, interactive statistics on facility usage and simple report generation minimize administrative tasks. Provided as a web application, Parkour is a convenient tool for both deep sequencing service users and laboratory personal. A set of web APIs allow coordinated information sharing with local and remote bioinformaticians. The flexible structure allows workflow customization and simple addition of new features as well as the expansion to other domains.AVAILABILITY AND IMPLEMENTATION:The code and documentation are available at https://github.com/maxplanck-ie/parkour.SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
Advanced age is not only a major risk factor for a range of disorders within an aging individual but may also enhance susceptibility for disease in the next generation. In humans, advanced paternal age has been associated with increased risk for a number of diseases. Experiments in rodent models have provided initial evidence that paternal age can influence behavioral traits in offspring animals, but the overall scope and extent of paternal age effects on health and disease across the life span remain underexplored. Here, we report that old father offspring mice showed a reduced life span and an exacerbated development of aging traits compared with young father offspring mice. Genome-wide epigenetic analyses of sperm from aging males and old father offspring tissue identified differentially methylated promoters, enriched for genes involved in the regulation of evolutionarily conserved longevity pathways. Gene expression analyses, biochemical experiments, and functional studies revealed evidence for an overactive mTORC1 signaling pathway in old father offspring mice. Pharmacological mTOR inhibition during the course of normal aging ameliorated many of the aging traits that were exacerbated in old father offspring mice. These findings raise the possibility that inherited alterations in longevity pathways contribute to intergenerational effects of aging in old father offspring mice.
Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) is an invaluable tool for mapping chromatin-associated proteins. However, sample preparation is still a largely individual and labor-intensive process that hinders assay throughput and comparability. Here, we present a novel method for ultra-parallelized high-throughput ChIP-seq that addresses the aforementioned problems. The method, called RELACS (Restriction Enzyme-based Labeling of Chromatin in Situ ), employs barcoding of chromatin within intact nuclei extracted from different sources (e.g. tissues, treatments, time points). Barcoded nuclei are pooled and processed within the same ChIP, for maximal comparability and significant workload reduction. The choice of user-friendly, straightforward, enzymatic steps for chromatin fragmentation and barcoding makes RELACS particularly suitable for implementation large-scale clinical studies and scarce samples. RELACS can generate ChIP-seq libraries from hundreds of samples within three days and with less than 1000 cells per sample.
Chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) is an invaluable tool for mapping chromatin-associated proteins. Current barcoding strategies aim to improve assay throughput and scalability but intense sample handling and lack of standardization over cell types, cell numbers and epitopes hinder wide-spread use in the field. Here, we present a barcoding method to enable high-throughput ChIP-seq using common molecular biology techniques. The method, called RELACS (restriction enzyme-based labeling of chromatin in situ) relies on standardized nuclei extraction from any source and employs chromatin cutting and barcoding within intact nuclei. Barcoded nuclei are pooled and processed within the same ChIP reaction, for maximal comparability and workload reduction. The innovative barcoding concept is particularly user-friendly and suitable for implementation to standardized large-scale clinical studies and scarce samples. Aiming to maximize universality and scalability, RELACS can generate ChIP-seq libraries for transcription factors and histone modifications from hundreds of samples within three days.
Maternal folic acid (FA) supplementation prior to and during gestation is recommended for the prevention of neural tube closure defects in the developing embryo. Prior studies, however, suggested that excessive FA supplementation during gestation can be associated with toxic effects on the developing organism. Here, we address whether maternal dietary folic acid supplementation at 40 mg/kg chow (FD), restricted to a period prior to conception, affects neurobehavioural development in the offspring generation. Detailed behavioural analyses showed reversal learning impairments in the Morris water maze in offspring derived from dams exposed to FD prior to conceiving. Furthermore, offspring of FD dams showed minor and transient gene expression differences relative to controls. Our data suggest that temporary exposure of female germ cells to FD is sufficient to cause impaired cognitive flexibility in the subsequent generation.
This release brings the following new features: Python 3.6 support cooler file format support PCA computation for A / B compartment visualizations of viewpoints many bug fixes change license to GPLv3
Dietary intake of methyl donors, such as folic acid and methionine, shows considerable intra-individual variation in human populations. While it is recognized that maternal departures from the optimum of dietary methyl donor intake can increase the risk for mental health issues and neurological disorders in offspring, it has not been explored whether paternal dietary methyl donor intake influences behavioral and cognitive functions in the next generation. Here, we report that elevated paternal dietary methyl donor intake in a mouse model, transiently applied prior to mating, resulted in offspring animals (methyl donor-rich diet (MD) F1 mice) with deficits in hippocampus-dependent learning and memory, impaired hippocampal synaptic plasticity and reduced hippocampal theta oscillations. Gene expression analyses revealed altered expression of the methionine adenosyltransferase Mat2a and BK channel subunit Kcnmb2, which was associated with changes in Kcnmb2 promoter methylation in MD F1 mice. Hippocampal overexpression of Kcnmb2 in MD F1 mice ameliorated altered spatial learning and memory, supporting a role of this BK channel subunit in the MD F1 behavioral phenotype. Behavioral and gene expression changes did not extend into the F2 offspring generation. Together, our data indicate that paternal dietary factors influence cognitive and neural functions in the offspring generation.
The primary problem with the explosion of biomedical datasets is not the data, not computational resources, and not the required storage space, but the general lack of trained and skilled researchers to manipulate and analyze these data. Eliminating this problem requires development of comprehensive educational resources. Here we present a community-driven framework that enables modern, interactive teaching of data analytics in life sciences and facilitates the development of training materials. The key feature of our system is that it is not a static but a continuously improved collection of tutorials. By coupling tutorials with a web-based analysis framework, biomedical researchers can learn by performing computation themselves through a web browser without the need to install software or search for example datasets. Our ultimate goal is to expand the breadth of training materials to include fundamental statistical and data science topics and to precipitate a complete re-engineering of undergraduate and graduate curricula in life sciences. This project is accessible at https://training.galaxyproject.org.
Dietary restriction regimes extend lifespan in various animal models. Here we show that longevity in male C57BL/6J mice subjected to every-other-day feeding is associated with a delayed onset of neoplastic disease that naturally limits lifespan in these animals. We compare more than 200 phenotypes in over 20 tissues in aged animals fed with a lifelong every-other-day feeding or ad libitum access to food diet to determine whether molecular, cellular, physiological and histopathological aging features develop more slowly in every-other-day feeding mice than in controls. We also analyze the effects of every-other-day feeding on young mice on shorter-term every-other-day feeding or ad libitum to account for possible aging-independent restriction effects. Our large-scale analysis reveals overall only limited evidence for a retardation of the aging rate in every-other-day feeding mice. The data indicate that every-other-day feeding-induced longevity is sufficiently explained by delays in life-limiting neoplastic disorders and is not associated with a more general slowing of the aging process in mice.Dietary restriction can extend the life of various model organisms. Here, Xie et al. show that intermittent periods of fasting achieved through every-other-day feeding protect mice against neoplastic disease but do not broadly delay organismal aging in animals.