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.

Age-Associated Changes in Endothelial Transcriptome and Epigenetic Landscapes Correlate With Elevated Risk of Cerebral Microbleeds.

Background Stroke is a leading global cause of human death and disability, with advanced aging associated with elevated incidences of stroke. Despite high mortality and morbidity of stroke, the mechanisms leading to blood-brain barrier dysfunction and development of stroke with age are poorly understood. In the vasculature of brain, endothelial cells (ECs) constitute the core component of the blood-brain barrier and provide a physical barrier composed of tight junctions, adherens junctions, and basement membrane. Methods and Results We show, in mice, the incidents of intracerebral bleeding increases with age. After isolating an enriched population of cerebral ECs from murine brains at 2, 6, 12, 18, and 24 months, we studied age-associated changes in gene expression. The study reveals age-dependent dysregulation of 1388 genes, including many involved in the maintenance of the blood-brain barrier and vascular integrity. We also investigated age-dependent changes on the levels of CpG methylation and accessible chromatin in cerebral ECs. Our study reveals correlations between age-dependent changes in chromatin structure and gene expression, whereas the dynamics of DNA methylation changes are different. Conclusions We find significant age-dependent downregulation of the Aplnr gene along with age-dependent reduction in chromatin accessibility of promoter region of the Aplnr gene in cerebral ECs. Aplnr is associated with positive regulation of vasodilation and is implicated in vascular health. Altogether, our data suggest a potential role of the apelinergic axis involving the ligand apelin and its receptor to be critical in maintenance of the blood-brain barrier and vascular integrity.

Ageing-associated small RNA cargo of extracellular vesicles.

Previous work on murine models and humans demonstrated global as well as tissue-specific molecular ageing trajectories of RNAs. Extracellular vesicles (EVs) are membrane vesicles mediating the horizontal transfer of genetic information between different tissues. We sequenced small regulatory RNAs (sncRNAs) in two mouse plasma fractions at five time points across the lifespan from 2-18 months: (1) sncRNAs that are free-circulating (fc-RNA) and (2) sncRNAs bound outside or inside EVs (EV-RNA). Different sncRNA classes exhibit unique ageing patterns that vary between the fcRNA and EV-RNA fractions. While tRNAs showed the highest correlation with ageing in both fractions, rRNAs exhibited inverse correlation trajectories between the EV- and fc-fractions. For miRNAs, the EV-RNA fraction was exceptionally strongly associated with ageing, especially the miR-29 family in adipose tissues. Sequencing of sncRNAs and coding genes in fat tissue of an independent cohort of aged mice up to 27 months highlighted the pivotal role of miR-29a-3p and miR-29b-3p in ageing-related gene regulation that we validated in a third cohort by RT-qPCR.

A Monoclonal Human Alveolar Epithelial Cell Line ("Arlo") with Pronounced Barrier Function for Studying Drug Permeability and Viral Infections.

In the development of orally inhaled drug products preclinical animal models regularly fail to predict pharmacological as well as toxicological responses in humans. Models based on human cells and tissues are potential alternatives to animal experimentation allowing for the isolation of essential processes of human biology and making them accessible in vitro. Here, the generation of a novel monoclonal cell line "Arlo," derived from the polyclonal human alveolar epithelium lentivirus immortalized cell line hAELVi via single-cell printing, and its characterization as a model for the human alveolar epithelium as well as a building block for future complex in vitro models is described. "Arlo" is systematically compared in vitro to primary human alveolar epithelial cells (hAEpCs) as well as to the polyclonal hAELVi cell line. "Arlo" cells show enhanced barrier properties with high transepithelial electrical resistance (TEER) of ≈3000 Ω cm2 and a potential difference (PD) of ≈30 mV under air-liquid interface (ALI) conditions, that can be modulated. The cells grow in a polarized monolayer and express genes relevant to barrier integrity as well as homeostasis as is observed in hAEpCs. Successful productive infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in a proof-of-principle study offers an additional, attractive application of "Arlo" beyond biopharmaceutical experimentation.

Transcriptional and Epigenetic Consequences of DMSO Treatment on HepaRG Cells.

Dimethyl sulfoxide (DMSO) is used to sustain or favor hepatocyte differentiation in vitro. Thus, DMSO is used in the differentiation protocol of the HepaRG cells that present the closest drug-metabolizing enzyme activities to primary human hepatocytes in culture. The aim of our study is to clarify its influence on liver-specific gene expression. For that purpose, we performed a large-scale analysis (gene expression and histone modification) to determine the global role of DMSO exposure during the differentiation process of the HepaRG cells. The addition of DMSO drives the upregulation of genes mainly regulated by PXR and PPARα whereas genes not affected by this addition are regulated by HNF1α, HNF4α, and PPARα. DMSO-differentiated-HepaRG cells show a differential expression for genes regulated by histone acetylation, while differentiated-HepaRG cells without DMSO show gene signatures associated with histone deacetylases. In addition, we observed an interplay between cytoskeleton organization and EMC remodeling with hepatocyte maturation.

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.

A comprehensive approach for genome-wide efficiency profiling of DNA modifying enzymes.

A precise understanding of DNA methylation dynamics is of great importance for a variety of biological processes including cellular reprogramming and differentiation. To date, complex integration of multiple and distinct genome-wide datasets is required to realize this task. We present GwEEP (genome-wide epigenetic efficiency profiling) a versatile approach to infer dynamic efficiencies of DNA modifying enzymes. GwEEP relies on genome-wide hairpin datasets, which are translated by a hidden Markov model into quantitative enzyme efficiencies with reported confidence around the estimates. GwEEP predicts de novo and maintenance methylation efficiencies of Dnmts and furthermore the hydroxylation efficiency of Tets. Its design also allows capturing further oxidation processes given available data. We show that GwEEP predicts accurately the epigenetic changes of ESCs following a Serum-to-2i shift and applied to Tet TKO cells confirms the hypothesized mutual interference between Dnmts and Tets.

Broad domains of histone marks in the highly compact Paramecium macronuclear genome.

The unicellular ciliate Paramecium contains a large vegetative macronucleus with several unusual characteristics, including an extremely high coding density and high polyploidy. As macronculear chromatin is devoid of heterochromatin, our study characterizes the functional epigenomic organization necessary for gene regulation and proper Pol II activity. Histone marks (H3K4me3, H3K9ac, H3K27me3) reveal no narrow peaks but broad domains along gene bodies, whereas intergenic regions are devoid of nucleosomes. Our data implicate H3K4me3 levels inside ORFs to be the main factor associated with gene expression, and H3K27me3 appears in association with H3K4me3 in plastic genes. Silent and lowly expressed genes show low nucleosome occupancy, suggesting that gene inactivation does not involve increased nucleosome occupancy and chromatin condensation. Because of a high occupancy of Pol II along highly expressed ORFs, transcriptional elongation appears to be quite different from that of other species. This is supported by missing heptameric repeats in the C-terminal domain of Pol II and a divergent elongation system. Our data imply that unoccupied DNA is the default state, whereas gene activation requires nucleosome recruitment together with broad domains of H3K4me3. In summary, gene activation and silencing in Paramecium run counter to the current understanding of chromatin biology.

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.

Strong Expansion of Human Regulatory T Cells for Adoptive Cell Therapy Results in Epigenetic Changes Which May Impact Their Survival and Function.

Adoptive transfer of regulatory T cells (Treg) is a promising new therapeutic option to treat detrimental inflammatory conditions after transplantation and during autoimmune disease. To reach sufficient cell yield for treatment, ex vivo isolated autologous or allogenic Tregs need to be expanded extensively in vitro during manufacturing of the Treg product. However, repetitive cycles of restimulation and prolonged culture have been shown to impact T cell phenotypes, functionality and fitness. It is therefore critical to scrutinize the molecular changes which occur during T cell product generation, and reexamine current manufacturing practices. We performed genome-wide DNA methylation profiling of cells throughout the manufacturing process of a polyclonal Treg product that has proven safety and hints of therapeutic efficacy in kidney transplant patients. We found progressive DNA methylation changes over the duration of culture, which were donor-independent and reproducible between manufacturing runs. Differentially methylated regions (DMRs) in the final products were significantly enriched at promoters and enhancers of genes implicated in T cell activation. Additionally, significant hypomethylation did also occur in promoters of genes implicated in functional exhaustion in conventional T cells, some of which, however, have been reported to strengthen immunosuppressive effector function in Tregs. At the same time, a set of reported Treg-specific demethylated regions increased methylation levels with culture, indicating a possible destabilization of Treg identity during manufacturing, which was independent of the purity of the starting material. Together, our results indicate that the repetitive TCR-mediated stimulation lead to epigenetic changes that might impact functionality of Treg products in multiple ways, by possibly shifting to an effector Treg phenotype with enhanced functional activity or by risking destabilization of Treg identity and impaired TCR activation. Our analyses also illustrate the value of epigenetic profiling for the evaluation of T cell product manufacturing pipelines, which might open new avenues for the improvement of current adoptive Treg therapies with relevance for conventional effector T cell products.

Rapid Base-Specific Calling of SARS-CoV-2 Variants of Concern Using Combined RT-PCR Melting Curve Screening and SIRPH Technology.

BACKGROUND: The emergence of novel variants of concern of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) demands fast and reliable detection of such variants in local populations. METHODS: Here we present a cost-efficient and fast workflow combining a prescreening of SARS-CoV-2-positive samples using reverse transcription polymerase chain reaction melting curve analysis with multiplexed IP-RP-HPLC-based single nucleotide primer extensions. RESULTS: The entire workflow from positive SARS-CoV-2 testing to base-specific identification of variants requires about 24 hours. CONCLUSIONS: We applied the sensitive method to monitor local variant of concern outbreaks in SARS-CoV-2-positive samples collected in a confined region of Germany.

The complete mitochondrial genome of the photosymbiotic sea slug Berghia stephanieae (Valdés, 2005) (Gastropoda, Nudibranchia).

Berghia stephanieae (Nudibranchia, Cladobranchia) is a photosymbiotic sea slug that feeds exclusively on sea anemones from the genus Exaiptasia. It then specifically incorporates dinoflagellates belonging to the Symbiodiniaceae obtained from their prey. Here, we present the complete mitochondrial genome sequence of B. stephanieae combining Oxford Nanopore long read and Illumina short-read sequencing data. The mitochondrial genome has a total length of 14,786 bp, it contains the 13 protein-encoding genes, 23 tRNAs, and two rRNAs and is similar to other nudibranchs except for the presence of a duplicated tRNA-Ser 1.

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.

Human IgA-Expressing Bone Marrow Plasma Cells Characteristically Upregulate Programmed Cell Death Protein-1 Upon B Cell Receptor Stimulation.

The functions of bone marrow plasma cells (BMPC) beyond antibody production are not fully elucidated and distinct subsets of BMPC suggest potential different functions. Phenotypic differences were identified for human BMPC depending on CD19 expression. Since CD19 is a co-stimulatory molecule of the B-cell-receptor (BCR), and IgA+ and IgM+ BMPC express the BCR on their surface, we here studied whether CD19 expression affects cellular responses, such as BCR signaling and the expression of checkpoint molecules. We analyzed 132 BM samples from individuals undergoing routine total hip arthroplasty. We found that both CD19+ and CD19- BMPC expressed BCR signaling molecules. Notably, the BCR-associated kinase spleen tyrosine kinase (SYK) including pSYK was higher expressed in CD19+ BMPC compared to CD19- BMPC. BCR stimulation also resulted in increased kinase phosphorylation downstream of the BCR while expression of CD19 remained stable afterwards. Interestingly, the BCR response was restricted to IgA+ BMPC independently of CD19 expression. With regard to the expression of checkpoint molecules, CD19- BMPC expressed higher levels of co-inhibitory molecule programmed cell death protein-1 (PD-1) than CD19+ BMPC. IgA+ BMPC characteristically upregulated PD-1 upon BCR stimulation in contrast to other PC subsets and inhibition of the kinase SYK abrogated PD-1 upregulation. In contrast, expression of PD-1 ligand, B and T lymphocyte attenuator (BTLA) and CD28 did not change upon BCR activation of IgA+ BMPC. Here, we identify a distinct characteristic of IgA+ BMPC that is independent of the phenotypic heterogeneity of the subsets according to their CD19 expression. The data suggest that IgA+ BMPC underlie different regulatory principles and/or exert distinct regulatory functions.

Targeted De-Methylation of the FOXP3-TSDR Is Sufficient to Induce Physiological FOXP3 Expression but Not a Functional Treg Phenotype.

CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

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.

DNA methylation signatures of monozygotic twins clinically discordant for multiple sclerosis.

Multiple sclerosis (MS) is an inflammatory, demyelinating disease of the central nervous system with a modest concordance rate in monozygotic twins, which strongly argues for involvement of epigenetic factors. We observe highly similar peripheral blood mononuclear cell-based methylomes in 45 MS-discordant monozygotic twins. Nevertheless, we identify seven MS-associated differentially methylated positions (DMPs) of which we validate two, including a region in the TMEM232 promoter and ZBTB16 enhancer. In CD4 + T cells we find an MS-associated differentially methylated region in FIRRE. Additionally, 45 regions show large methylation differences in individual pairs, but they do not clearly associate with MS. Furthermore, we present epigenetic biomarkers for current interferon-beta treatment, and extensive validation shows that the ZBTB16 DMP is a signature for prior glucocorticoid treatment. Taken together, this study represents an important reference for epigenomic MS studies, identifies new candidate epigenetic markers, and highlights treatment effects and genetic background as major confounders.

A comprehensive analysis of 195 DNA methylomes reveals shared and cell-specific features of partially methylated domains.

BACKGROUND: Partially methylated domains are extended regions in the genome exhibiting a reduced average DNA methylation level. They cover gene-poor and transcriptionally inactive regions and tend to be heterochromatic. We present a comprehensive comparative analysis of partially methylated domains in human and mouse cells, to identify structural and functional features associated with them. RESULTS: Partially methylated domains are present in up to 75% of the genome in human and mouse cells irrespective of their tissue or cell origin. Each cell type has a distinct set of partially methylated domains, and genes expressed in such domains show a strong cell type effect. The methylation level varies between cell types with a more pronounced effect in differentiating and replicating cells. The lowest level of methylation is observed in highly proliferating and immortal cancer cell lines. A decrease of DNA methylation within partially methylated domains tends to be linked to an increase in heterochromatic histone marks and a decrease of gene expression. Characteristic combinations of heterochromatic signatures in partially methylated domains are linked to domains of early and middle S-phase and late S-G2 phases of DNA replication. CONCLUSIONS: Partially methylated domains are prominent signatures of long-range epigenomic organization. Integrative analysis identifies them as important general, lineage- and cell type-specific topological features. Changes in partially methylated domains are hallmarks of cell differentiation, with decreased methylation levels and increased heterochromatic marks being linked to enhanced cell proliferation. In combination with broad histone marks, partially methylated domains demarcate distinct domains of late DNA replication.

Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.