Transmission of trained immunity and heterologous resistance to infections across generations.

Intergenerational inheritance of immune traits linked to epigenetic modifications has been demonstrated in plants and invertebrates. Here we provide evidence for transmission of trained immunity across generations to murine progeny that survived a sublethal systemic infection with Candida albicans or a zymosan challenge. The progeny of trained mice exhibited cellular, developmental, transcriptional and epigenetic changes associated with the bone marrow-resident myeloid effector and progenitor cell compartment. Moreover, the progeny of trained mice showed enhanced responsiveness to endotoxin challenge, alongside improved protection against systemic heterologous Escherichia coli and Listeria monocytogenes infections. Sperm DNA of parental male mice intravenously infected with the fungus C. albicans showed DNA methylation differences linked to immune gene loci. These results provide evidence for inheritance of trained immunity in mammals, enhancing protection against infections.

Peptide-specific recognition of human cytomegalovirus strains controls adaptive natural killer cells.

Natural killer (NK) cells are innate lymphocytes that lack antigen-specific rearranged receptors, a hallmark of adaptive lymphocytes. In some people infected with human cytomegalovirus (HCMV), an NK cell subset expressing the activating receptor NKG2C undergoes clonal-like expansion that partially resembles anti-viral adaptive responses. However, the viral ligand that drives the activation and differentiation of adaptive NKG2C+ NK cells has remained unclear. Here we found that adaptive NKG2C+ NK cells differentially recognized distinct HCMV strains encoding variable UL40 peptides that, in combination with pro-inflammatory signals, controlled the population expansion and differentiation of adaptive NKG2C+ NK cells. Thus, we propose that polymorphic HCMV peptides contribute to shaping of the heterogeneity of adaptive NKG2C+ NK cell populations among HCMV-seropositive people.

Dynamic readers for 5-(hydroxy)methylcytosine and its oxidized derivatives.

Tet proteins oxidize 5-methylcytosine (mC) to generate 5-hydroxymethyl (hmC), 5-formyl (fC), and 5-carboxylcytosine (caC). The exact function of these oxidative cytosine bases remains elusive. We applied quantitative mass-spectrometry-based proteomics to identify readers for mC and hmC in mouse embryonic stem cells (mESC), neuronal progenitor cells (NPC), and adult mouse brain tissue. Readers for these modifications are only partially overlapping, and some readers, such as Rfx proteins, display strong specificity. Interactions are dynamic during differentiation, as for example evidenced by the mESC-specific binding of Klf4 to mC and the NPC-specific binding of Uhrf2 to hmC, suggesting specific biological roles for mC and hmC. Oxidized derivatives of mC recruit distinct transcription regulators as well as a large number of DNA repair proteins in mouse ES cells, implicating the DNA damage response as a major player in active DNA demethylation.

LAG-3 Inhibitory Receptor Expression Identifies Immunosuppressive Natural Regulatory Plasma Cells.

B lymphocytes can suppress immunity through interleukin (IL)-10 production in infectious, autoimmune, and malignant diseases. Here, we have identified a natural plasma cell subset that distinctively expresses the inhibitory receptor LAG-3 and mediates this function in vivo. These plasma cells also express the inhibitory receptors CD200, PD-L1, and PD-L2. They develop from various B cell subsets in a B cell receptor (BCR)-dependent manner independently of microbiota in naive mice. After challenge they upregulate IL-10 expression via a Toll-like receptor-driven mechanism within hours and without proliferating. This function is associated with a unique transcriptome and epigenome, including the lowest amount of DNA methylation at the Il10 locus compared to other B cell subsets. Their augmented accumulation in naive mutant mice with increased BCR signaling correlates with the inhibition of memory T cell formation and vaccine efficacy after challenge. These natural regulatory plasma cells may be of broad relevance for disease intervention.

EpiSegMix: a flexible distribution hidden Markov model with duration modeling for chromatin state discovery.

MOTIVATION: Automated chromatin segmentation based on ChIP-seq (chromatin immunoprecipitation followed by sequencing) data reveals insights into the epigenetic regulation of chromatin accessibility. Existing segmentation methods are constrained by simplifying modeling assumptions, which may have a negative impact on the segmentation quality. RESULTS: We introduce EpiSegMix, a novel segmentation method based on a hidden Markov model with flexible read count distribution types and state duration modeling, allowing for a more flexible modeling of both histone signals and segment lengths. In a comparison with existing tools, ChromHMM, Segway, and EpiCSeg, we show that EpiSegMix is more predictive of cell biology, such as gene expression. Its flexible framework enables it to fit an accurate probabilistic model, which has the potential to increase the biological interpretability of chromatin states. AVAILABILITY AND IMPLEMENTATION: Source code: https://gitlab.com/rahmannlab/episegmix.

Mutual Regulation of Transcriptomes between Murine Pneumocytes and Fibroblasts Mediates Alveolar Regeneration in Air-Liquid Interface Cultures.

Alveolar type 2 and club cells are part of the stem cell niche of the lung and their differentiation is required for pulmonary homeostasis and tissue regeneration. A disturbed crosstalk between fibroblasts and epithelial cells contributes to the loss of lung structure in chronic lung diseases. Therefore, it is important to understand how fibroblasts and lung epithelial cells interact during regeneration. Here, we analyzed the interaction of fibroblasts and the alveolar epithelium modeled in air-liquid interface cultures. Single-cell transcriptomics showed that cocultivation with fibroblasts leads to increased expression of type 2 markers in pneumocytes, activation of regulons associated with the maintenance of alveolar type 2 cells (e.g., Etv5), and transdifferentiation of club cells toward pneumocytes. This was accompanied by an intensified transepithelial barrier. Vice versa, the activation of NF-κB pathways and the CEBPB regulon and the expression of IL-6 and other differentiation factors (e.g., fibroblast growth factors) were increased in fibroblasts cocultured with epithelial cells. Recombinant IL-6 enhanced epithelial barrier formation. Therefore, in our coculture model, regulatory loops were identified by which lung epithelial cells mediate regeneration and differentiation of the alveolar epithelium in a cooperative manner with the mesenchymal compartment.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states.

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Resident memory CD4+ T lymphocytes mobilize from bone marrow to contribute to a systemic secondary immune reaction.

Resident memory T lymphocytes (TRM ) of epithelial tissues and the Bm protect their host tissue. To what extent these cells are mobilized and contribute to systemic immune reactions is less clear. Here, we show that in secondary immune reactions to the measles-mumps-rubella (MMR) vaccine, CD4+ TRM are mobilized into the blood within 16 to 48 h after immunization in humans. This mobilization of TRM is cognate: TRM recognizing other antigens are not mobilized, unless they cross-react with the vaccine. We also demonstrate through methylome analyses that TRM are mobilized from the Bm. These mobilized cells make significant contribution to the systemic immune reaction, as evidenced by their T-cell receptor Vß clonotypes represented among the newly generated circulating memory T-cells, 14 days after vaccination. Thus, TRM of the Bm confer not only local, but also systemic immune memory.

DNA methylation biomarkers in colorectal cancer: Clinical applications for precision medicine.

Colorectal cancer (CRC) is the second leading cause of cancer death worldwide that is attributed to gradual long-term accumulation of both genetic and epigenetic changes. To reduce the mortality rate of CRC and to improve treatment efficacy, it will be important to develop accurate noninvasive diagnostic tests for screening, acute and personalized diagnosis. Epigenetic changes such as DNA methylation play an important role in the development and progression of CRC. Over the last decade, a panel of DNA methylation markers has been reported showing a high accuracy and reproducibility in various semi-invasive or noninvasive biosamples. Research to obtain comprehensive panels of markers allowing a highly sensitive and differentiating diagnosis of CRC is ongoing. Moreover, the epigenetic alterations for cancer therapy, as a precision medicine strategy will increase their therapeutic potential over time. Here, we discuss the current state of DNA methylation-based biomarkers and their impact on CRC diagnosis. We emphasize the need to further identify and stratify methylation-biomarkers and to develop robust and effective detection methods that are applicable for a routine clinical setting of CRC diagnostics particularly at the early stage of the disease.

Identification of an FXR-modulated liver-intestine hybrid state in iPSC-derived hepatocyte-like cells.

BACKGROUND & AIMS: Pluripotent stem cell (PSC)-derived hepatocyte-like cells (HLC) have enormous potential as a replacement for primary hepatocytes in drug screening, toxicology and cell replacement therapy, but their genome-wide expression patterns differ strongly from primary human hepatocytes (PHH). METHODS: We differentiated human induced pluripotent stem cells (hiPSC) via definitive endoderm to HLC and characterized the cells by single-cell and bulk RNA-seq, with complementary epigenetic analyses. We then compared HLC to PHH and publicly available data on human fetal hepatocytes (FH) ex vivo; we performed bioinformatics-guided interventions to improve HLC differentiation via lentiviral transduction of the nuclear receptor FXR and agonist exposure. RESULTS: Single-cell RNA-seq revealed that transcriptomes of individual HLC display a hybrid state, where hepatocyte-associated genes are expressed in concert with genes that are not expressed in PHH - mostly intestinal genes - within the same cell. Bulk-level overrepresentation analysis, as well as regulon analysis at the single-cell level, identified sets of regulatory factors discriminating HLC, FH, and PHH, hinting at a central role for the nuclear receptor FXR in the functional maturation of HLC. Combined FXR expression plus agonist exposure enhanced the expression of hepatocyte-associated genes and increased the ability of bile canalicular secretion as well as lipid droplet formation, thereby increasing HLCs' similarity to PHH. The undesired non-liver gene expression was reproducibly decreased, although only by a moderate degree. CONCLUSION: In contrast to physiological hepatocyte precursor cells and mature hepatocytes, HLC co-express liver and hybrid genes in the same cell. Targeted modification of the FXR gene regulatory network improves their differentiation by suppressing intestinal traits whilst inducing hepatocyte features. LAY SUMMARY: Generation of human hepatocytes from stem cells represents an active research field but its success is hampered by the fact that the stem cell-derived 'hepatocytes' still show major differences to hepatocytes obtained from a liver. Here, we identified an important reason for the difference, specifically that the stem cell-derived 'hepatocyte' represents a hybrid cell with features of hepatocytes and intestinal cells. We show that a specific protein (FXR) suppresses intestinal and induces liver features, thus bringing the stem cell-derived cells closer to hepatocytes derived from human livers.

Suv39h-dependent H3K9me3 marks intact retrotransposons and silences LINE elements in mouse embryonic stem cells.

Heterochromatin is required to restrict aberrant expression of retrotransposons, but it remains poorly defined due to the underlying repeat-rich sequences. We dissected Suv39h-dependent histone H3 lysine 9 trimethylation (H3K9me3) by genome-wide ChIP sequencing in mouse embryonic stem cells (ESCs). Refined bioinformatic analyses of repeat subfamilies indicated selective accumulation of Suv39h-dependent H3K9me3 at interspersed repetitive elements that cover ∼5% of the ESC epigenome. The majority of the ∼8,150 intact long interspersed nuclear elements (LINEs) and endogenous retroviruses (ERVs), but only a minor fraction of the >1.8 million degenerate and truncated LINEs/ERVs, are enriched for Suv39h-dependent H3K9me3. Transcriptional repression of intact LINEs and ERVs is differentially regulated by Suv39h and other chromatin modifiers in ESCs but governed by DNA methylation in committed cells. These data provide a function for Suv39h-dependent H3K9me3 chromatin to specifically repress intact LINE elements in the ESC epigenome.

Quantitative comparison of within-sample heterogeneity scores for DNA methylation data.

DNA methylation is an epigenetic mark with important regulatory roles in cellular identity and can be quantified at base resolution using bisulfite sequencing. Most studies are limited to the average DNA methylation levels of individual CpGs and thus neglect heterogeneity within the profiled cell populations. To assess this within-sample heterogeneity (WSH) several window-based scores that quantify variability in DNA methylation in sequencing reads have been proposed. We performed the first systematic comparison of four published WSH scores based on simulated and publicly available datasets. Moreover, we propose two new scores and provide guidelines for selecting appropriate scores to address cell-type heterogeneity, cellular contamination and allele-specific methylation. Most of the measures were sensitive in detecting DNA methylation heterogeneity in these scenarios, while we detected differences in susceptibility to technical bias. Using recently published DNA methylation profiles of Ewing sarcoma samples, we show that DNA methylation heterogeneity provides information complementary to the DNA methylation level. WSH scores are powerful tools for estimating variance in DNA methylation patterns and have the potential for detecting novel disease-associated genomic loci not captured by established statistics. We provide an R-package implementing the WSH scores for integration into analysis workflows.

GeneTrail 3: advanced high-throughput enrichment analysis.

We present GeneTrail 3, a major extension of our web service GeneTrail that offers rich functionality for the identification, analysis, and visualization of deregulated biological processes. Our web service provides a comprehensive collection of biological processes and signaling pathways for 12 model organisms that can be analyzed with a powerful framework for enrichment and network analysis of transcriptomic, miRNomic, proteomic, and genomic data sets. Moreover, GeneTrail offers novel workflows for the analysis of epigenetic marks, time series experiments, and single cell data. We demonstrate the capabilities of our web service in two case-studies, which highlight that GeneTrail is well equipped for uncovering complex molecular mechanisms. GeneTrail is freely accessible at: http://genetrail.bioinf.uni-sb.de.

Epigenetic control of region-specific transcriptional programs in mouse cerebellar and cortical astrocytes.

Astrocytes from the cerebral cortex (CTX) and cerebellum (CB) share basic molecular programs, but also form distinct spatial and functional subtypes. The regulatory epigenetic layers controlling such regional diversity have not been comprehensively investigated so far. Here, we present an integrated epigenome analysis of methylomes, open chromatin, and transcriptomes of astroglia populations isolated from the cortex or cerebellum of young adult mice. Besides a basic overall similarity in their epigenomic programs, cortical astrocytes and cerebellar astrocytes exhibit substantial differences in their overall open chromatin structure and in gene-specific DNA methylation. Regional epigenetic differences are linked to differences in transcriptional programs encompassing genes of region-specific transcription factor networks centered around Lhx2/Foxg1 in CTX astrocytes and the Zic/Irx families in CB astrocytes. The distinct epigenetic signatures around these transcription factor networks point to a complex interconnected and combinatorial regulation of region-specific transcriptomes. These findings suggest that key transcription factors, previously linked to temporal, regional, and spatial control of neurogenesis, also form combinatorial networks important for astrocytes. Our study provides a valuable resource for the molecular basis of regional astrocyte identity and physiology.

Neuronal deficiency of p38α-MAPK ameliorates symptoms and pathology of APP or Tau-transgenic Alzheimer's mouse models.

Alzheimer's disease (AD) is the leading cause of dementia with very limited therapeutic options. Amyloid ß (Aß) and phosphorylated Tau (p-Tau) are key pathogenic molecules in AD. P38α-MAPK is specifically activated in AD lesion sites. However, its effects on AD pathogenesis, especially on p-Tau-associated brain pathology, and the underlying molecular mechanisms remain unclear. We mated human APP-transgenic mice and human P301S Tau-transgenic mice with mapk14-floxed and neuron-specific Cre-knock-in mice. We observed that deletion of p38α-MAPK specifically in neurons improves the cognitive function of both 9-month-old APP and Tau-transgenic AD mice, which is associated with decreased Aß and p-Tau load in the brain. We further used next-generation sequencing to analyze the gene transcription in brains of p38α-MAPK deficient and wild-type APP-transgenic mice, which indicated that deletion of p38α-MAPK regulates the transcription of calcium homeostasis-related genes, especially downregulates the expression of grin2a, a gene encoding NMDAR subunit NR2A. Cell culture experiments further verified that deletion of p38α-MAPK inhibits NMDA-triggered calcium influx and neuronal apoptosis. Our systemic studies of AD pathogenic mechanisms using both APP- and Tau-transgenic mice suggested that deletion of neuronal p38α-MAPK attenuates AD-associated brain pathology and protects neurons in AD pathogenesis. This study supports p38α-MAPK as a novel target for AD therapy.

Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data.

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-seq, ATAC-seq and NOMe-seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However, we also find a high proportion of assay specific NDRs that are often 'called' by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-seq and DNase I-seq provide a superb high NDR calling rate for relatively low sequencing costs in comparison to NOMe-seq, NOMe-seq singles out for its genome-wide coverage allowing to not only detect NDRs but also endogenous DNA methylation and as we show here genome wide segmentation into heterochromatic B domains and local phasing of nucleosomes outside of NDRs. In summary, our comparisons strongly suggest to consider assay specific differences for the experimental design and for generalized and comparative functional interpretations.

The dynamics of genome-wide DNA methylation reprogramming in mouse primordial germ cells.

Genome-wide DNA methylation reprogramming occurs in mouse primordial germ cells (PGCs) and preimplantation embryos, but the precise dynamics and biological outcomes are largely unknown. We have carried out whole-genome bisulfite sequencing (BS-Seq) and RNA-Seq across key stages from E6.5 epiblast to E16.5 PGCs. Global loss of methylation takes place during PGC expansion and migration with evidence for passive demethylation, but sequences that carry long-term epigenetic memory (imprints, CpG islands on the X chromosome, germline-specific genes) only become demethylated upon entry of PGCs into the gonads. The transcriptional profile of PGCs is tightly controlled despite global hypomethylation, with transient expression of the pluripotency network, suggesting that reprogramming and pluripotency are inextricably linked. Our results provide a framework for the understanding of the epigenetic ground state of pluripotency in the germline.

Two are better than one: HPoxBS - hairpin oxidative bisulfite sequencing.

The controlled and stepwise oxidation of 5mC to 5hmC, 5fC and 5caC by Tet enzymes is influencing the chemical and biological properties of cytosine. Besides direct effects on gene regulation, oxidised forms influence the dynamics of demethylation and re-methylation processes. So far, no combined methods exist which allow to precisely determine the strand specific localisation of cytosine modifications along with their CpG symmetric distribution. Here we describe a comprehensive protocol combining conventional hairpin bisulfite with oxidative bisulfite sequencing (HPoxBS) to determine the strand specific distribution of 5mC and 5hmC at base resolution. We apply this method to analyse the contribution of local oxidative effects on DNA demethylation in mouse ES cells. Our method includes the HPoxBS workflow and subsequent data analysis using our developed software tools. Besides a precise estimation and display of strand specific 5mC and 5hmC levels at base resolution we apply the data to predict region specific activities of Dnmt and Tet enzymes. Our experimental and computational workflow provides a precise double strand display of 5mC and 5hmC modifications at single base resolution. Based on our data we predict region specific Tet and Dnmt enzyme efficiencies shaping the distinct locus levels and patterns of 5hmC and 5mC.