Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

Single-cell (sc) RNA, ATAC and Multiome sequencing became powerful tools for uncovering biological and disease mechanisms. Unfortunately, manual analysis of sc data presents multiple challenges due to large data volumes and complexity of configuration parameters. This complexity, as well as not being able to reproduce a computational environment, affects the reproducibility of analysis results. The Scientific Data Analysis Platform (https://SciDAP.com) allows biologists without computational expertise to analyze sequencing-based data using portable and reproducible pipelines written in Common Workflow Language (CWL). Our suite of computational pipelines addresses the most common needs in scRNA-Seq, scATAC-Seq and scMultiome data analysis. When executed on SciDAP, it offers a user-friendly alternative to manual data processing, eliminating the need for coding expertise. In this protocol, we describe the use of SciDAP to analyze scMultiome data. Similar approaches can be used for analysis of scRNA-Seq, scATAC-Seq and scVDJ-Seq datasets.

Aryl hydrocarbon receptor and IL-13 signaling crosstalk in human keratinocytes and atopic dermatitis

Introduction

Atopic dermatitis (AD) is an allergic skin disease mediated by skin barrier impairment and IL-13-driven immune response. Activation of the aryl hydrocarbon receptor (AHR) has shown promise in early clinical trials for AD; however, the mechanism by which AHR partially ameliorates AD is not well known.

Methods

Gene expression data from human biopsies were analyzed, and compared to gene expression from RNA-sequencing in our in-vitro HaCaT cell model system. Western blot, ELISA qRT-PCR were used to further explore the relationship between AHR and IL-13 signaling in HaCaT cells.

Results

The AHR target gene CYP1A1 was decreased in lesional skin compared with healthy control skin (p = 4.30 × 10-9). Single-cell RNA sequencing (scRNAseq) demonstrated increased AHR expression (p < 1.0 × 10-4) and decreased CYP1A1 expression in lesional AD keratinocytes compared with healthy control keratinocytes (p < 0.001). Activation of AHR by AHR agonists in HaCaT cells reversed IL-13-dependent gene expression of several key genes in AD pathogenesis, most notably the eosinophil chemoattractant CCL26 (eotaxin-3). Differentially expressed genes in keratinocytes of patients with AD substantially overlapped with genes regulated by AHR agonists from HaCaT cells by RNAseq, but in reverse direction. Mechanistically, there was evidence for direct transcriptional effects of AHR; AHR binding motifs were identified in the differentially expressed genes from lesional AD keratinocytes compared to control keratinocytes, and AHR activation did not modify IL-13-dependent signal transducer and activator of transcription 6 (STAT6) translocation to the nucleus.

Discussion

Together, these data suggest that the AHR pathway is dysregulated in AD and that AHR modulates IL-13 downstream signaling in keratinocytes through genome-wide, transcriptional regulatory effects.

TET1 regulates responses to house dust mite by altering chromatin accessibility, DNA methylation, and gene expression in airway epithelial cells

Background

Previous studies have identified TET1 as a potential key regulator of genes linked to asthma. TET1 has been shown to transcriptionally respond to house dust mite extract, an allergen known to directly cause allergic asthma development, and regulate the expression of genes involved in asthma. How TET1 regulates expression of these genes, however, is unknown. TET1 is a DNA demethylase; therefore, most prior research on TET1-based gene regulation has focused on how TET1 affects methylation. However, TET1 can also interact directly with transcription factors and histone modifiers to regulate gene expression. Understanding how TET1 regulates expression to contribute to allergic responses and asthma development thus requires a comprehensive approach. To this end, we measured mRNA expression, DNA methylation, chromatin accessibility and histone modifications in control and TET1 knockdown human bronchial epithelial cells treated or untreated with house dust mite extract.

Results

Throughout our analyses, we detected strong similarities between the effects of TET1 knockdown alone and the effects of HDM treatment alone. One especially striking pattern was that both TET1 knockdown and HDM treatment generally led to decreased chromatin accessibility at largely the same genomic loci. Transcription factor enrichment analyses indicated that altered chromatin accessibility following the loss of TET1 may affect, or be affected by, CTCF and CEBP binding. TET1 loss also led to changes in DNA methylation, but these changes were generally in regions where accessibility was not changing.

Conclusions

TET1 regulates gene expression through different mechanisms (DNA methylation and chromatin accessibility) in different parts of the genome in the airway epithelial cells, which mediates inflammatory responses to allergen. Collectively, our data suggest novel molecular mechanisms through which TET1 regulates critical pathways following allergen challenges and contributes to the development of asthma.

TSLP shapes the pathogenic responses of memory CD4+ T cells in eosinophilic esophagitis

The cytokine thymic stromal lymphopoietin (TSLP) mediates type 2 immune responses, and treatments that interfere with TSLP activity are in clinical use for asthma. Here, we investigated whether TSLP contributes to allergic inflammation by directly stimulating human CD4+ T cells and whether this process is operational in eosinophilic esophagitis (EoE), a disease linked to variants in TSLP. We showed that about 10% of esophageal-derived memory CD4+ T cells from individuals with EoE and less than 3% of cells from control individuals expressed the receptor for TSLP and directly responded to TSLP, as determined by measuring the phosphorylation of STAT5, a transcription factor activated downstream of TSLP stimulation. Accordingly, increased numbers of TSLP-responsive memory CD4+ T cells were present in the circulation of individuals with EoE. TSLP increased the proliferation of CD4+ T cells, enhanced type 2 cytokine production, induced the increased abundance of its own receptor, and modified the expression of 212 genes. The epigenetic response to TSLP was associated with an enrichment in BATF and IRF4 chromatin-binding sites, and these transcription factors were induced by TSLP, providing a feed-forward loop. The numbers of circulating and esophageal CD4+ T cells responsive to TSLP correlated with the numbers of esophageal eosinophils, supporting a potential functional role for TSLP in driving the pathogenesis of EoE and providing the basis for a blood-based diagnostic test based on the extent of TSLP-induced STAT5 phosphorylation in circulating CD4+ T cells. These findings highlight the potential therapeutic value of TSLP inhibitors for the treatment of EoE.

PRC1 suppresses a female gene regulatory network to ensure testicular differentiation

Gonadal sex determination and differentiation are controlled by somatic support cells of testes (Sertoli cells) and ovaries (granulosa cells). In testes, the epigenetic mechanism that maintains chromatin states responsible for suppressing female sexual differentiation remains unclear. Here, we show that Polycomb repressive complex 1 (PRC1) suppresses a female gene regulatory network in postnatal Sertoli cells. We genetically disrupted PRC1 function in embryonic Sertoli cells after sex determination, and we found that PRC1-depleted postnatal Sertoli cells exhibited defective proliferation and cell death, leading to the degeneration of adult testes. In adult Sertoli cells, PRC1 suppressed specific genes required for granulosa cells, thereby inactivating the female gene regulatory network. Chromatin regions associated with female-specific genes were marked by Polycomb-mediated repressive modifications: PRC1-mediated H2AK119ub and PRC2-mediated H3K27me3. Taken together, this study identifies a critical Polycomb-based mechanism that suppresses ovarian differentiation and maintains Sertoli cell fate in adult testes.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

Vitamin D receptor and STAT6 interactome governs oesophageal epithelial barrier responses to IL-13 signalling

OBJECTIVE

The contribution of vitamin D (VD) deficiency to the pathogenesis of allergic diseases remains elusive. We aimed to define the impact of VD on oesophageal allergic inflammation.

DESIGN

We assessed the genomic distribution and function of VD receptor (VDR) and STAT6 using histology, molecular imaging, motif discovery and metagenomic analysis. We examined the role of VD supplementation in oesophageal epithelial cells, in a preclinical model of IL-13-induced oesophageal allergic inflammation and in human subjects with eosinophilic oesophagitis (EoE).

RESULTS

VDR response elements were enriched in oesophageal epithelium, suggesting enhanced VDR binding to functional gene enhancer and promoter regions. Metagenomic analysis showed that VD supplementation reversed dysregulation of up to 70% of the transcriptome and epigenetic modifications (H3K27Ac) induced by IL-13 in VD-deficient cells, including genes encoding the transcription factors HIF1A and SMAD3, endopeptidases (SERPINB3) and epithelial-mesenchymal transition mediators (TGFBR1, TIAM1, SRC, ROBO1, CDH1). Molecular imaging and chromatin immunoprecipitation showed VDR and STAT6 colocalisation within the regulatory regions of the affected genes, suggesting that VDR and STAT6 interactome governs epithelial tissue responses to IL-13 signalling. Indeed, VD supplementation reversed IL-13-induced epithelial hyperproliferation, reduced dilated intercellular spaces and barrier permeability, and improved differentiation marker expression (filaggrin, involucrin). In a preclinical model of IL-13-mediated oesophageal allergic inflammation and in human EoE, VD levels inversely associated with severity of oesophageal eosinophilia and epithelial histopathology.

CONCLUSIONS

Collectively, these findings identify VD as a natural IL-13 antagonist with capacity to regulate the oesophageal epithelial barrier functions, providing a novel therapeutic entry point for type 2 immunity-related diseases.

An atlas of gene regulatory networks for memory CD4 + T cells in youth and old age

Aging profoundly affects immune-system function, promoting susceptibility to pathogens, cancers and chronic inflammation. We previously identified a population of IL-10-producing, T follicular helper-like cells (" Tfh10 "), linked to suppressed vaccine responses in aged mice. Here, we integrate single-cell ( sc )RNA-seq, scATAC-seq and genome-scale modeling to characterize Tfh10 - and the full CD4 + memory T cell ( CD4 + TM ) compartment - in young and old mice. We identified 13 CD4 + TM populations, which we validated through cross-comparison to prior scRNA-seq studies. We built gene regulatory networks ( GRNs ) that predict transcription-factor control of gene expression in each T-cell population and how these circuits change with age. Through integration with pan-cell aging atlases, we identified intercellular-signaling networks driving age-dependent changes in CD4 + TM. Our atlas of finely resolved CD4 + TM subsets, GRNs and cell-cell communication networks is a comprehensive resource of predicted regulatory mechanisms operative in memory T cells, presenting new opportunities to improve immune responses in the elderly.

maxATAC: Genome-scale transcription-factor binding prediction from ATAC-seq with deep neural networks

Transcription factors read the genome, fundamentally connecting DNA sequence to gene expression across diverse cell types. Determining how, where, and when TFs bind chromatin will advance our understanding of gene regulatory networks and cellular behavior. The 2017 ENCODE-DREAM in vivo Transcription-Factor Binding Site (TFBS) Prediction Challenge highlighted the value of chromatin accessibility data to TFBS prediction, establishing state-of-the-art methods for TFBS prediction from DNase-seq. However, the more recent Assay-for-Transposase-Accessible-Chromatin (ATAC)-seq has surpassed DNase-seq as the most widely-used chromatin accessibility profiling method. Furthermore, ATAC-seq is the only such technique available at single-cell resolution from standard commercial platforms. While ATAC-seq datasets grow exponentially, suboptimal motif scanning is unfortunately the most common method for TFBS prediction from ATAC-seq. To enable community access to state-of-the-art TFBS prediction from ATAC-seq, we (1) curated an extensive benchmark dataset (127 TFs) for ATAC-seq model training and (2) built "maxATAC", a suite of user-friendly, deep neural network models for genome-wide TFBS prediction from ATAC-seq in any cell type. With models available for 127 human TFs, maxATAC is the largest collection of high-performance TFBS prediction models for ATAC-seq. maxATAC performance extends to primary cells and single-cell ATAC-seq, enabling improved TFBS prediction in vivo. We demonstrate maxATAC's capabilities by identifying TFBS associated with allele-dependent chromatin accessibility at atopic dermatitis genetic risk loci.

Role of transcription factor RUNX1 in the maintenance of T helper cell epigenome and rapid recall response

Memory T cells exhibit an accelerated response after a second encounter with the same antigen. In our previous studies, we showed that the ability to produce effector molecules within minutes of the secondary exposure to antigen (“rapid recall response”) is associated with the epigenetic state of the rapid recall genes. This memory epigenome is established during the primary T cell encounter with its cognate antigen. To maintain a “rapid recall response” even when antigen is eliminated, resting memory T cells would have to maintain open chromatin signature for years despite episodic homeostatic T cell proliferation. However, the mechanisms underlying the long-term maintenance of this epigenetic pattern remains unclear. Here, we used 10X Multiome (scRNA-seq and scATAC-seq obtained from the same cells) analysis to show that the association between the chromatin state of memory cells and their cytokine expression upon activation is Th lineage specific.

Bioinformatic analysis showed that open chromatin specific for memory but not naïve CD4 T cells is enriched with motifs of constitutively expressed transcription factors RUNX, SP1, and ETS. Our preliminary Cut&Tag data in total human resting CD4 T cells demonstrates that RUNX1 binds DNA close to T helper cytokines genes (IL3, IL4, IL9, IL10, IL13, IFNG) at “memory-specific” ATAC-seq regions, suggesting the role of RUNX1 in maintaining this “rapid recall” epigenome. Together with the dynamic model of transcriptome and epigenome rearrangements this data could help us to understand how memory T cells establish a platform that enables and maintains their unique ability to respond more rapidly to subsequent encounters with antigen.

Supported by R01AI153442-01A1 from NIAID, NIH

Abstract PO-068: Cholesterol auxotrophy promotes the expansion of centroacinar cells giving rise to the basal subtype of pancreatic adenocarcinoma

Gene expression analyses established at least two molecular subtypes of pancreatic adenocarcinoma (PDAC), the classical (or glandular), and the basal (or mesenchymal), each of which is associated with distinct prognoses and sensitivity to chemotherapy. It remains unclear, however, whether the basal carcinoma cells arise from a separate cell-of-origin, or are emerging from the pre-existing well-differentiated “classical” PDAC cells. To distinguish these alternatives, we conducted single-cell transcriptome analyses and virtual lineage tracing comparing cellular populations at pre-malignant stages in basal versus classical PDAC mouse models. We previously reported that chemical or genetic inhibition of the cholesterol biosynthetic pathway in KrasG12D; Trp53 (KPPC) mice predisposes to basal rather than glandular PDAC development because of the pancreas-specific increased sterol response element-binding protein 1 (SREBP1) activity and TGFβ signaling that induces cancer cell stemness and the EMT (PMID: 32976774). Pancreas-selective knockout of a conditional allele of cholesterol pathway gene, Nsdhl (NAD(P)-dependent steroid dehydrogenase-like), renders pancreatic epithelial cells cholesterol auxotrophs and drives basal PDAC in the majority of animals (KPPCN mice). At 5-6 weeks of age, grossly and microscopically tumor-free pancreatic tissues were selected for single-cell isolation and single-cell RNA sequencing (scRNA seq) using the 10X platform. After standard filtering and sample normalization procedures, downstream analyses included identification of relevant cell clusters using Seurat, lineage tracing algorithms, and in silico modeling of autocrine and paracrine signaling interactions between subsets of PDAC and non-malignant cells. Our key findings are as follows: 1) premalignant KPPCN pancreata exhibit a massive expansion of cancer-associated fibroblasts (CAFs) of predominantly inflammatory differentiation (iCAFs); 2) despite relatively fewer ADM and PanIN pre-malignant lesions in KPPCN compared to KPPC, scRNA seq identifies the significant expansion of epithelial cells with features of centroacinar and stem-like cells (increased expression of Aldh1a2, Nes, Sox9, Ly6a, Cxcl12, and Met); these centroacinar-like cells, while retaining epithelial identity (Epcam, Cdh1), also exhibit features of pluripotency by co-expression of Ins2 and other stem cell markers; 3) alignment with basal PDAC (KPPCN) and classical (KPPC) carcinoma cell populations strongly suggests the continuity of clonal evolution of the centroacinar-like cells towards the basal PDAC. While our genetic model does not recapitulate the multiple alternative pathways leading to basal PDAC development, cholesterol auxotrophy via SREBP1 may be a factor governing the expansion of undifferentiated precursors, which via interactions with cancer-promoting iCAFs, drive basal PDAC development.

Citation Format: Michael Kotliar, Aizhan Surumbayeva, Linara Gabitova, Suraj Peri, Diana Restifo, Kathy Q. Cai, Artem Barski, Igor Astsaturov. Cholesterol auxotrophy promotes the expansion of centroacinar cells giving rise to the basal subtype of pancreatic adenocarcinoma [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2021 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2021;81(22 Suppl):Abstract nr PO-068.

Aiolos regulates eosinophil migration into tissues

Expression of Ikaros family transcription factor IKZF3 (Aiolos) increases during murine eosinophil lineage commitment and maturation. Herein, we investigated Aiolos expression and function in mature human and murine eosinophils. Murine eosinophils deficient in Aiolos demonstrated gene expression changes in pathways associated with granulocyte-mediated immunity, chemotaxis, degranulation, ERK/MAPK signaling, and extracellular matrix organization; these genes had ATAC peaks within 1 kB of the TSS that were enriched for Aiolos-binding motifs. Global Aiolos deficiency reduced eosinophil frequency within peripheral tissues during homeostasis; a chimeric mouse model demonstrated dependence on intrinsic Aiolos expression by eosinophils. Aiolos deficiency reduced eosinophil CCR3 surface expression, intracellular ERK1/2 signaling, and CCL11-induced actin polymerization, emphasizing an impaired functional response. Aiolos-deficient eosinophils had reduced tissue accumulation in chemokine-, antigen-, and IL-13-driven inflammatory experimental models, all of which at least partially depend on CCR3 signaling. Human Aiolos expression was associated with active chromatin marks enriched for IKZF3, PU.1, and GATA-1-binding motifs within eosinophil-specific histone ChIP-seq peaks. Furthermore, treating the EOL-1 human eosinophilic cell line with lenalidomide yielded a dose-dependent decrease in Aiolos. These collective data indicate that eosinophil homing during homeostatic and inflammatory allergic states is Aiolos-dependent, identifying Aiolos as a potential therapeutic target for eosinophilic disease.

Runx Transcription Factors in T Cells-What Is Beyond Thymic Development?

Runx proteins (also known as Runt-domain transcription factors) have been studied for a long time as key regulators of cellular differentiation. RUNX2 has been described as essential for osteogenesis, whereas RUNX1 and RUNX3 are known to control blood cell development during different stages of cell lineage specification. However, recent studies show evidence of complex relationships between RUNX proteins, chromatin-modifying machinery, the cytoskeleton and different transcription factors in various non-embryonic contexts, including mature T cell homeostasis, inflammation and cancer. In this review, we discuss the diversity of Runx functions in mature T helper cells, such as production of cytokines and chemokines by different CD4 T cell populations; apoptosis; and immunologic memory acquisition. We then briefly cover recent findings about the contribution of RUNX1, RUNX2 and RUNX3 to various immunologic diseases. Finally, we discuss areas that require further study to better understand the role that Runx proteins play in inflammation and immunity.

Epigenetic Analysis of the Chromatin Landscape Identifies a Repertoire of Murine Eosinophil-Specific PU.1-Bound Enhancers

Eosinophils develop in the bone marrow from hematopoietic progenitors into mature cells capable of a plethora of immunomodulatory roles via the choreographed process of eosinophilopoiesis. However, the gene regulatory elements and transcription factors (TFs) orchestrating this process remain largely unknown. The potency and resulting diversity fundamental to an eosinophil's complex immunomodulatory functions and tissue specialization likely result from dynamic epigenetic regulation of the eosinophil genome, a dynamic eosinophil regulome. In this study, we applied a global approach using broad-range, next-generation sequencing to identify a repertoire of eosinophil-specific enhancers. We identified over 8200 active enhancers located within 1-20 kB of expressed eosinophil genes. TF binding motif analysis revealed PU.1 (Spi1) motif enrichment in eosinophil enhancers, and chromatin immunoprecipitation coupled with massively parallel sequencing confirmed PU.1 binding in likely enhancers of genes highly expressed in eosinophils. A substantial proportion (>25%) of these PU.1-bound enhancers were unique to murine, culture-derived eosinophils when compared among enhancers of highly expressed genes of three closely related myeloid cell subsets (macrophages, neutrophils, and immature granulocytes). Gene ontology analysis of eosinophil-specific, PU.1-bound enhancers revealed enrichment for genes involved in migration, proliferation, degranulation, and survival. Furthermore, eosinophil-specific superenhancers were enriched in genes whose homologs are associated with risk loci for eosinophilia and allergic diseases. Our collective data identify eosinophil-specific enhancers regulating key eosinophil genes through epigenetic mechanisms (H3K27 acetylation) and TF binding (PU.1).

Methylation quantitative trait locus analysis of chronic postsurgical pain uncovers epigenetic mediators of genetic risk

Background: Overlap of pathways enriched by single nucleotide polymorphisms and DNA-methylation underlying chronic postsurgical pain (CPSP), prompted pilot study of CPSP-associated methylation quantitative trait loci (meQTL). Materials & methods: Children undergoing spine-fusion were recruited prospectively. Logistic-regression for genome- and epigenome-wide CPSP association and DNA-methylation-single nucleotide polymorphism association/mediation analyses to identify meQTLs were followed by functional genomics analyses. Results: CPSP (n = 20/58) and non-CPSP groups differed in pain-measures. Of 2753 meQTLs, DNA-methylation at 127 cytosine-guanine dinucleotides mediated association of 470 meQTLs with CPSP (p < 0.05). At PARK16 locus, CPSP risk meQTLs were associated with decreased DNA-methylation at RAB7L1 and increased DNA-methylation at PM20D1. Corresponding RAB7L1/PM20D1 blood eQTLs (GTEx) and cytosine-guanine dinucleotide-loci enrichment for histone marks, transcription factor binding sites and ATAC-seq peaks suggest altered transcription factor-binding. Conclusion: CPSP-associated meQTLs indicate epigenetic mechanisms mediate genetic risk. Clinical trial registration: NCT01839461, NCT01731873 (ClinicalTrials.gov).

Super-enhancer switching drives a burst in gene expression at the mitosis-to-meiosis transition

Due to bursts in the expression of thousands of germline-specific genes, the testis has the most diverse and complex transcriptome of all organs. By analyzing the male germline of mice, we demonstrate that the genome-wide reorganization of super-enhancers (SEs) drives bursts in germline gene expression after the mitosis-to-meiosis transition. SE reorganization is regulated by two molecular events: the establishment of meiosis-specific SEs via A-MYB (MYBL1), a key transcription factor for germline genes, and the resolution of SEs in mitotically proliferating cells via SCML2, a germline-specific Polycomb protein required for spermatogenesis-specific gene expression. Prior to entry into meiosis, meiotic SEs are preprogrammed in mitotic spermatogonia to ensure the unidirectional differentiation of spermatogenesis. We identify key regulatory factors for both mitotic and meiotic enhancers, revealing a molecular logic for the concurrent activation of mitotic enhancers and suppression of meiotic enhancers in the somatic and/or mitotic proliferation phases.

AP-1 activity induced by co-stimulation is required for chromatin opening during T cell activation

Activation of T cells is dependent on the organized and timely opening and closing of chromatin. Herein, we identify AP-1 as the transcription factor that directs most of this remodeling. Chromatin accessibility profiling showed quick opening of closed chromatin in naive T cells within 5 h of activation. These newly opened regions were strongly enriched for the AP-1 motif, and indeed, ChIP-seq demonstrated AP-1 binding at >70% of them. Broad inhibition of AP-1 activity prevented chromatin opening at AP-1 sites and reduced the expression of nearby genes. Similarly, induction of anergy in the absence of co-stimulation during activation was associated with reduced induction of AP-1 and a failure of proper chromatin remodeling. The translational relevance of these findings was highlighted by the substantial overlap of AP-1-dependent elements with risk loci for multiple immune diseases, including multiple sclerosis, inflammatory bowel disease, and allergic disease. Our findings define AP-1 as the key link between T cell activation and chromatin remodeling.

Epigenomics of T Cell Activation and Memory

A characteristic feature of adaptive immunity is the rapid recall response of memory T cells, the ability of memory T cells to produce cytokines in response to pathogen much more rapidly than naïve T cells. We have recently found that this rapid recall ability is associated with epigenetic gene poising. Here, we aim to uncover the mechanism that underlie chromatin remodeling during initial activation of T cells and, thus, contribute to the establishment of memory phenotype. Upon encountering an antigen, naïve T helper cells become activated, differentiate into several lineages and contribute to immune response. Activation of T cells is dependent on organized and timely opening and closing of chromatin. Herein, we identify AP-1 as the transcription factor that directs most of this remodeling. Chromatin accessibility profiling showed quick opening of closed chromatin in naïve T cells within hours of activation. These newly open regions were strongly enriched for the AP-1 motif, and indeed, ChIP-Seq demonstrated AP-1 binding at more than 70% of them. Broad inhibition of AP-1 activity prevented chromatin opening at AP-1 sites and reduced expression of nearby genes. Similarly, induction of anergy in the absence of co-stimulation during activation, was associated with reduced induction of AP-1 and a failure of proper chromatin remodeling. The translational relevance of these findings was highlighted by the substantial overlap of AP-1–dependent elements with risk loci for multiple immune diseases, most notably multiple sclerosis, IBD and allergy. Our findings define AP-1 as the key link between T cell activation and chromatin remodeling.

KMT2C/D COMPASS complex-associated diseases [KCDCOM-ADs]: an emerging class of congenital regulopathies

The type 2 lysine methyltransferases KMT2C and KMT2D are large, enzymatically active scaffold proteins that form the core of nuclear regulatory structures known as KMT2C/D COMPASS complexes (complex of proteins associating with Set1). These evolutionarily conserved proteins regulate DNA promoter and enhancer elements, modulating the activity of diverse cell types critical for embryonic morphogenesis, central nervous system development, and post-natal survival. KMT2C/D COMPASS complexes and their binding partners enhance active gene expression of specific loci via the targeted modification of histone-3 tail residues, in general promoting active euchromatic conformations. Over the last 20 years, mutations in five key COMPASS complex genes have been linked to three human congenital syndromes: Kabuki syndrome (type 1 [KMT2D] and 2 [KDM6A]), Rubinstein-Taybi syndrome (type 1 [CBP] and 2 [EP300]), and Kleefstra syndrome type 2 (KMT2C). Here, we review the composition and biochemical function of the KMT2 complexes. The specific cellular and embryonic roles of the KMT2C/D COMPASS complex are highlight with a focus on clinically relevant mechanisms sensitive to haploinsufficiency. The phenotypic similarities and differences between the members of this new family of disorders are outlined and emerging therapeutic strategies are detailed.

Identification of anoctamin 1 (ANO1) as a key driver of esophageal epithelial proliferation in eosinophilic esophagitis

BACKGROUND

The pathology of eosinophilic esophagitis (EoE) is characterized by eosinophil-rich inflammation, basal zone hyperplasia (BZH), and dilated intercellular spaces, and the underlying processes that drive the pathologic manifestations of the disease remain largely unexplored.

OBJECTIVE

We sought to investigate the involvement of the calcium-activated chloride channel anoctamin 1 (ANO1) in esophageal proliferation and the histopathologic features of EoE.

METHODS

We examined mRNA and protein expression of ANO1 in esophageal biopsy samples from patients with EoE and in mice with EoE. We performed molecular and cellular analyses and ion transport assays on an in vitro esophageal epithelial 3-dimensional model system (EPC2-ALI) and murine models of EoE to define the relationship between expression and function of ANO1 and esophageal epithelial proliferation in patients with EoE.

RESULTS

We observed increased ANO1 expression in esophageal biopsy samples from patients with EoE and in mice with EoE. ANO1 was expressed within the esophageal basal zone, and expression correlated positively with disease severity (eosinophils/high-power field) and BZH. Using an in vitro esophageal epithelial 3-dimensional model system revealed that ANO1 undergoes chromatin modification and rapid upregulation of expression after IL-13 stimulation, that ANO1 is the primary apical IL-13-induced Cl- transport mechanism within the esophageal epithelium, and that loss of ANO1-dependent Cl- transport abrogated esophageal epithelial proliferation. Mechanistically, ANO1-dependent regulation of basal cell proliferation was associated with modulation of TP63 expression and phosphorylated cyclin-dependent kinase 2 levels.

CONCLUSIONS

These data identify a functional role for ANO1 in esophageal cell proliferation and BZH in patients with EoE and provide a rationale for pharmacologic intervention of ANO1 function in patients with EoE.

Dynamic reorganization of open chromatin underlies diverse transcriptomes during spermatogenesis

During spermatogenesis, germ cells undergo massive cellular reconstruction and dynamic chromatin remodeling to facilitate highly diverse transcriptomes, which are required for the production of functional sperm. However, it remains unknown how germline chromatin is organized to promote the dynamic, complex transcriptomes of spermatogenesis. Here, using ATAC-seq, we establish the varied landscape of open chromatin during spermatogenesis. We identify the reorganization of accessible chromatin in intergenic and intronic regions during the mitosis-to-meiosis transition. During the transition, mitotic-type open chromatin is closed while the de novo formation of meiotic-type open chromatin takes place. Contrastingly, differentiation processes such as spermatogonial differentiation and the meiosis-to-postmeiosis transition involve chromatin closure without the de novo formation of accessible chromatin. In spermiogenesis, the germline-specific Polycomb protein SCML2 promotes the closure of open chromatin at autosomes for gene suppression. Paradoxically, we identify the massive de novo formation of accessible chromatin when the sex chromosomes undergo meiotic sex chromosome inactivation, and this is also mediated by SCML2. These results reveal meiotic sex chromosome inactivation as an active process for chromatin organization. Together, our results unravel the genome-wide, dynamic reorganization of open chromatin and reveal mechanisms that underlie diverse transcriptomes during spermatogenesis.

CWL-Airflow: a lightweight pipeline manager supporting Common Workflow Language

BACKGROUND

Massive growth in the amount of research data and computational analysis has led to increased use of pipeline managers in biomedical computational research. However, each of the >100 such managers uses its own way to describe pipelines, leading to difficulty porting workflows to different environments and therefore poor reproducibility of computational studies. For this reason, the Common Workflow Language (CWL) was recently introduced as a specification for platform-independent workflow description, and work began to transition existing pipelines and workflow managers to CWL.

FINDINGS

Herein, we present CWL-Airflow, a package that adds support for CWL to the Apache Airflow pipeline manager. CWL-Airflow uses CWL version 1.0 specification and can run workflows on stand-alone MacOS/Linux servers, on clusters, or on a variety of cloud platforms. A sample CWL pipeline for processing of chromatin immunoprecipitation sequencing data is provided.

CONCLUSIONS

CWL-Airflow will provide users with the features of a fully fledged pipeline manager and the ability to execute CWL workflows anywhere Airflow can run-from a laptop to a cluster or cloud environment. CWL-Airflow is available under Apache License, version 2.0 (Apache-2.0), and can be downloaded from https://barski-lab.github.io/cwl-airflow, https://scicrunch.org/resolver/RRID:SCR_017196.

Co-Stimulation–Induced AP-1 Activity is Required for Chromatin Opening During T Cell Activation

Activation of T cells is dependent on organized and timely opening and closing of chromatin. Herein, we identify AP-1 as the transcription factor that directs most of this remodeling. Chromatin accessibility profiling showed quick opening of closed chromatin in naïve T cells within 5 hours of activation. These newly open regions were strongly enriched for the AP-1 motif, and indeed, ChIP-seq demonstrated AP-1 binding at more than 70% of them. Broad inhibition of AP-1 activity prevented chromatin opening at AP-1 sites and reduced expression of nearby genes. Similarly, induction of anergy in the absence of co-stimulation during activation, was associated with reduced induction of AP-1 and a failure of proper chromatin remodeling. The translational relevance of these findings was highlighted by the substantial overlap of AP-1–dependent elements with risk loci for multiple immune diseases, most notably multiple sclerosis. Our findings define AP-1 as the key link between T cell activation and chromatin remodeling.

Attenuated chromatin compartmentalization in meiosis and its maturation in sperm development

Germ cells manifest a unique gene expression program and regain totipotency in the zygote. Here, we perform Hi-C analysis to examine 3D chromatin organization in male germ cells during spermatogenesis. We show that the highly compartmentalized 3D chromatin organization characteristic of interphase nuclei is attenuated in meiotic prophase. Meiotic prophase is predominated by short-range intrachromosomal interactions that represent a condensed form akin to that of mitotic chromosomes. However, unlike mitotic chromosomes, meiotic chromosomes display weak genomic compartmentalization, weak topologically associating domains, and localized point interactions in prophase. In postmeiotic round spermatids, genomic compartmentalization increases and gives rise to the strong compartmentalization seen in mature sperm. The X chromosome lacks domain organization during meiotic sex-chromosome inactivation. We propose that male meiosis occurs amid global reprogramming of 3D chromatin organization and that strengthening of chromatin compartmentalization takes place in spermiogenesis to prepare the next generation of life.

A rapidly evolved domain, the SCML2 DNA-binding repeats, contributes to chromatin binding of mouse SCML2†

Genes involved in sexual reproduction diverge rapidly as a result of reproductive fitness. Here, we identify a novel protein domain in the germline-specific Polycomb protein SCML2 that is required for the establishment of unique gene expression programs after the mitosis-to-meiosis transition in spermatogenesis. We term this novel domain, which is comprised of rapidly evolved, DNA-binding repeat units of 28 amino acids, the SCML2 DNA-binding (SDB) repeats. These repeats are acquired in a specific subgroup of the rodent lineage, having been subjected to positive selection in the course of evolution. Mouse SCML2 has two DNA-binding domains: one is the SDB repeats and the other is an RNA-binding region, which is conserved in human SCML2. For the recruitment of SCML2 to target loci, the SDB repeats cooperate with the other functional domains of SCML2 to bind chromatin. The cooperative action of these domains enables SCML2 to sense DNA hypomethylation in an in vivo chromatin environment, thereby enabling SCML2 to bind to hypomethylated chromatin. We propose that the rapid evolution of SCML2 is due to reproductive adaptation, which has promoted species-specific gene expression programs in spermatogenesis.

Microbiota-sensitive epigenetic signature predicts inflammation in Crohn's disease

Altered response to the intestinal microbiota strongly associates with inflammatory bowel disease (IBD); however, how commensal microbial cues are integrated by the host during the pathogenesis of IBD is not understood. Epigenetics represents a potential mechanism that could enable intestinal microbes to modulate transcriptional output during the development of IBD. Here, we reveal a histone methylation signature of intestinal epithelial cells isolated from the terminal ilea of newly diagnosed pediatric IBD patients. Genes characterized by significant alterations in histone H3-lysine 4 trimethylation (H3K4me3) showed differential enrichment in pathways involving immunoregulation, cell survival and signaling, and metabolism. Interestingly, a large subset of these genes was epigenetically regulated by microbiota in mice and several microbiota-sensitive epigenetic targets demonstrated altered expression in IBD patients. Remarkably though, a substantial proportion of these genes exhibited H3K4me3 levels that correlated with the severity of intestinal inflammation in IBD, despite lacking significant differential expression. Collectively, these data uncover a previously unrecognized epigenetic profile of IBD that can be primed by commensal microbes and indicate sensitive targets in the epithelium that may underlie how microbiota predispose to subsequent intestinal inflammation and disease.

SCML2 promotes heterochromatin organization in late spermatogenesis

Spermatogenesis involves the progressive reorganization of heterochromatin. However, the mechanisms that underlie the dynamic remodeling of heterochromatin remain unknown. Here, we identify SCML2, a germline-specific Polycomb protein, as a critical regulator of heterochromatin organization in spermatogenesis. We show that SCML2 accumulates on pericentromeric heterochromatin (PCH) in male germ cells, where it suppresses PRC1-mediated monoubiquitylation of histone H2A at Lysine 119 (H2AK119ub) and promotes deposition of PRC2-mediated H3K27me3 during meiosis. In postmeiotic spermatids, SCML2 is required for heterochromatin organization, and the loss of SCML2 leads to the formation of ectopic patches of facultative heterochromatin. Our data suggest that, in the absence of SCML2, the ectopic expression of somatic lamins drives this process. Furthermore, the centromere protein CENP-V is a specific marker of PCH in postmeiotic spermatids, and SCML2 is required for CENP-V localization on PCH. Given the essential functions of PRC1 and PRC2 for genome-wide gene expression in spermatogenesis, our data suggest that heterochromatin organization and spermatogenesis-specific gene expression are functionally linked. We propose that SCML2 coordinates the organization of heterochromatin and gene expression through the regulation of Polycomb complexes.

Reuse of public, genome-wide, murine eosinophil expression data for hypotheses development

The eosinophil (Eos) surface phenotype and activation state is altered after recruitment into tissues and after exposure to pro-inflammatory cytokines. In addition, distinct Eos functional subsets have been described, suggesting that tissue-specific responses for Eos contribute to organ homeostasis. Understanding the mechanisms by which Eos subsets achieve their tissue-specific identity is currently an unmet goal for the eosinophil research community. Publicly archived expression data can be used to answer original questions, test and generate new hypotheses, and serve as a launching point for experimental design. With these goals in mind, we investigated the effect of genetic background, culture methods, and tissue residency on murine Eos gene expression using publicly available, genome-wide expression data. Eos differentiated from cultures have a gene expression profile that is distinct from that of native homeostatic Eos; thus, researchers can repurpose published expression data to aid in selecting the appropriate culture method to study their gene of interest. In addition, we identified Eos lung- and gastrointestinal-specific transcriptomes, highlighting the profound effect of local tissue environment on gene expression in a terminally differentiated granulocyte even at homeostasis. Expanding the "toolbox" of Eos researchers to include public-data reuse can reduce redundancy, increase research efficiency, and lead to new biological insights.

AP-1 transcription factor remodels chromatin during T cell activation

Upon encountering an antigen, naïve T helper cells activate, differentiate into several lineages and contribute to immune response. Epigenetic changes at cytokine and other effector genes during activation have been described, but the underlying mechanism is not understood. We profiled open chromatin during human T cell activation using ATAC-Seq and identified regions of epigenetic changes genome-wide. Open chromatin correlated with expression at nearby genes, which were enriched for functions related to T cell activation and immune response. As AP-1 and NFAT motifs were enriched in these newly open chromatin regions, we profiled their binding using ChIP-Seq. AP-1 was present at the majority of activation-specific open chromatin regions, which were often co-bound by both transcription factors. Notably, super enhancer formation was associated with AP-1 binding and chromatin opening during T cell activation. We also found that autoimmune disease-associated SNPs were significantly enriched in the newly open chromatin regions. Our data suggest that AP-1 activity contributes to the formation of the activated T cell epigenome and that mutations at the activation-specific open chromatin regions may contribute to autoimmune diseases.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. We show that transcription factors (TFs) occupy multiple loci of individual complex genetic disorders using novel computational methods. Application to 213 phenotypes and 1,544 TF binding datasets identifies 2,264 relationships between hundreds of TFs and 94 phenotypes, including AR in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of the systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many co-clustering human TFs, revealing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis, and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent upon EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new paradigms for disease origins.

TSLP signaling in CD4+ T cells programs a pathogenic T helper 2 cell state

Pathogenic T helper 2 (T_H2) cells, which produce increased amounts of the cytokines interleukin-5 (IL-5) and IL-13, promote allergic disorders, including asthma. Thymic stromal lymphopoietin (TSLP), a cytokine secreted by epithelial and innate immune cells, stimulates such pathogenic T_H2 cell responses. We found that TSLP signaling in mouse CD4⁺ T cells initiated transcriptional changes associated with T_H2 cell programming. IL-4 signaling amplified and stabilized the genomic response of T cells to TSLP, which increased the frequency of T cells producing IL-4, IL-5, and IL-13. Furthermore, the TSLP- and IL-4-programmed T_H2 cells had a pathogenic phenotype, producing greater amounts of IL-5 and IL-13 and other proinflammatory cytokines than did T_H2 cells stimulated with IL-4 alone. TSLP-mediated T_H2 cell induction involved distinct molecular pathways, including activation of the transcription factor STAT5 through the kinase JAK2 and repression of the transcription factor BCL6. Mice that received wild-type CD4⁺ T cells had exacerbated pathogenic T_H2 cell responses upon exposure to house dust mites compared to mice that received TSLP receptor-deficient CD4⁺ T cells. Transient TSLP signaling stably programmed pathogenic potential in memory T_H2 cells. In human CD4⁺ T cells, TSLP and IL-4 promoted the generation of T_H2 cells that produced greater amounts of IL-5 and IL-13. Compared to healthy controls, asthmatic children showed enhancement of such T cell responses in peripheral blood. Our data support a sequential cytokine model for pathogenic T_H2 cell differentiation and provide a mechanistic basis for the therapeutic targeting of TSLP signaling in human allergic diseases.

RNF8 and SCML2 cooperate to regulate ubiquitination and H3K27 acetylation for escape gene activation on the sex chromosomes

The sex chromosomes are enriched with germline genes that are activated during the late stages of spermatogenesis. Due to meiotic sex chromosome inactivation (MSCI), these sex chromosome-linked genes must escape silencing for activation in spermatids, thereby ensuring their functions for male reproduction. RNF8, a DNA damage response protein, and SCML2, a germline-specific Polycomb protein, are two major, known regulators of this process. Here, we show that RNF8 and SCML2 cooperate to regulate ubiquitination during meiosis, an early step to establish active histone modifications for subsequent gene activation. Double mutants of Rnf8 and Scml2 revealed that RNF8-dependent monoubiquitination of histone H2A at Lysine 119 (H2AK119ub) is deubiquitinated by SCML2, demonstrating interplay between RNF8 and SCML2 in ubiquitin regulation. Additionally, we identify distinct functions of RNF8 and SCML2 in the regulation of ubiquitination: SCML2 deubiquitinates RNF8-independent H2AK119ub but does not deubiquitinate RNF8-dependent polyubiquitination. RNF8-dependent polyubiquitination is required for the establishment of H3K27 acetylation, a marker of active enhancers, while persistent H2AK119ub inhibits establishment of H3K27 acetylation. Following the deposition of H3K27 acetylation, H3K4 dimethylation is established as an active mark on poised promoters. Together, we propose a model whereby regulation of ubiquitin leads to the organization of poised enhancers and promoters during meiosis, which induce subsequent gene activation from the otherwise silent sex chromosomes in postmeiotic spermatids.

Polycomb directs timely activation of germline genes in spermatogenesis

Maezawa et al. show that Polycomb-repressive complex 1 (PRC1) directs timely activation of germline genes during spermatogenesis.

During spermatogenesis, a large number of germline genes essential for male fertility are coordinately activated. However, it remains unknown how timely activation of this group of germline genes is accomplished. Here we show that Polycomb-repressive complex 1 (PRC1) directs timely activation of germline genes during spermatogenesis. Inactivation of PRC1 in male germ cells results in the gradual loss of a stem cell population and severe differentiation defects, leading to male infertility. In the stem cell population, RNF2, the dominant catalytic subunit of PRC1, activates transcription of Sall4, which codes for a transcription factor essential for subsequent spermatogenic differentiation. Furthermore, RNF2 and SALL4 together occupy transcription start sites of germline genes in the stem cell population. Once differentiation commences, these germline genes are activated to enable the progression of spermatogenesis. Our study identifies a novel mechanism by which Polycomb directs the developmental process by activating a group of lineage-specific genes.

Xist RNA repeat E is essential for ASH2L recruitment to the inactive X and regulates histone modifications and escape gene expression

Long non-coding RNA Xist plays a crucial role in establishing and maintaining X-chromosome inactivation (XCI) which is a paradigm of long non-coding RNA-mediated gene regulation. Xist has Xist-specific repeat elements A-F which are conserved among eutherian mammals, underscoring their functional importance. Here we report that Xist RNA repeat E, a conserved Xist repeat element in the Xist exon 7, interacts with ASH2L and contributes to maintenance of escape gene expression level on the inactive X-chromosome (Xi) during XCI. The Xist repeat E-deletion mutant female ES cells show the depletion of ASH2L from the Xi upon differentiation. Furthermore, a subset of escape genes exhibits unexpectedly higher expression in the repeat E mutant cells than the cells expressing wildtype Xist during X-inactivation, whereas the silencing of X-linked non-escape genes is not affected. We discuss the implications of these results to understand the role of ASH2L and Xist repeat E for histone modifications and escape gene regulation during random X-chromosome inactivation.

AP-1 transcription factor reprograms T cell epigenome during activation

Upon encountering an antigen, naïve T helper (Th) cells are activated, differentiate into several lineages and contribute to immune response. Epigenetic changes at several cytokine and other effector genes during activation has been described previously, but the mechanism behind these changes is not understood. Here, we have profiled open chromatin during human T cell activation using ATAC-Seq approach and identified the regions where epigenetic changes take place. T cell activation requires antigen signaling via T cell receptor and co-stimulation via CD28 which in turn result in nuclear translocation of NFAT and AP-1 transcription factors. For this reason, we profiled genome-wide distribution of these proteins during activation by ChIP-Seq. Our results show that 73% of activation-induced open chromatin loci contain binding sites for AP-1, whereas 39% of them bind NFAT. Over 90 % of the binding sites for NFAT were shared with AP-1. We and others have previously shown that the lack of co-stimulation results in reduction or delay of AP-1 translocation into the nucleus. In order to confirm the role of AP-1 in establishment of the open chromatin structure during activation we mapped open chromatin in the cells that were activated without co-stimulation. Our data showed that ATAC signal was reduced at the open chromatin sites bound by AP-1 in the absence of co-stimulation, suggesting that CD28 signaling and AP-1 transcription factor are involved in open chromatin formation during T cell activation.

Genomic integration of Wnt/β-catenin and BMP/Smad1 signaling coordinates foregut and hindgut transcriptional programs

Digestive system development is orchestrated by combinatorial signaling interactions between endoderm and mesoderm, but how these signals are interpreted in the genome is poorly understood. Here we identified the transcriptomes of Xenopus foregut and hindgut progenitors, which are conserved with mammals. Using RNA-seq and ChIP-seq we show that BMP/Smad1 regulates dorsal-ventral gene expression in both the endoderm and mesoderm, whereas Wnt/β-catenin acts as a genome-wide toggle between foregut and hindgut programs. Unexpectedly, β-catenin and Smad1 binding were associated with both transcriptional activation and repression, with Wnt-repressed genes often lacking canonical Tcf DNA binding motifs, suggesting a novel mode of direct repression. Combinatorial Wnt and BMP signaling was mediated by Smad1 and β-catenin co-occupying hundreds of cis-regulatory DNA elements, and by a crosstalk whereby Wnt negatively regulates BMP ligand expression in the foregut. These results extend our understanding of gastrointestinal organogenesis and of how Wnt and BMP might coordinate genomic responses in other contexts.

Rapid Recall Ability of Memory T cells is Encoded in their Epigenome

Even though T-cell receptor (TCR) stimulation together with co-stimulation is sufficient for the activation of both naïve and memory T cells, the memory cells are capable of producing lineage specific cytokines much more rapidly than the naïve cells. The mechanisms behind this rapid recall response of the memory cells are still not completely understood. Here, we performed epigenetic profiling of human resting naïve, central and effector memory T cells using ChIP-Seq and found that unlike the naïve cells, the regulatory elements of the cytokine genes in the memory T cells are marked by activating histone modifications even in the resting state. Therefore, the ability to induce expression of rapid recall genes upon activation is associated with the deposition of positive histone modifications during memory T cell differentiation. We propose a model of T cell memory, in which immunological memory state is encoded epigenetically, through poising and transcriptional memory.

Polycomb repressive complex 1 controls uterine decidualization

Uterine stromal cell decidualization is an essential part of the reproductive process. Decidual tissue development requires a highly regulated control of the extracellular tissue remodeling; however the mechanism of this regulation remains unknown. Through systematic expression studies, we detected that Cbx4/2, Rybp, and Ring1B [components of polycomb repressive complex 1 (PRC1)] are predominantly utilized in antimesometrial decidualization with polyploidy. Immunofluorescence analyses revealed that PRC1 members are co-localized with its functional histone modifier H2AK119ub1 (mono ubiquitination of histone-H2A at lysine-119) in polyploid cell. A potent small-molecule inhibitor of Ring1A/B E3-ubiquitin ligase or siRNA-mediated suppression of Cbx4 caused inhibition of H2AK119ub1, in conjunction with perturbation of decidualization and polyploidy development, suggesting a role for Cbx4/Ring1B-containing PRC1 in these processes. Analyses of genetic signatures by RNA-seq studies showed that the inhibition of PRC1 function affects 238 genes (154 up and 84 down) during decidualization. Functional enrichment analyses identified that about 38% genes primarily involved in extracellular processes are specifically targeted by PRC1. Furthermore, ~15% of upregulated genes exhibited a significant overlap with the upregulated Bmp2 null-induced genes in mice. Overall, Cbx4/Ring1B-containing PRC1 controls decidualization via regulation of extracellular gene remodeling functions and sheds new insights into underlying molecular mechanism(s) through transcriptional repression regulation.

BioWardrobe: an integrated platform for analysis of epigenomics and transcriptomics data

Development of next generation sequencing (NGS) has revolutionized molecular biology by enhancing our ability to perform genome-wide studies. However, due to the need for bioinformatics expertise and the sheer size of resulting datasets, use of these technologies is still beyond the capabilities of many laboratories. Here we present BioWardrobe, a platform that allows storing, visualizing and analyzing epigenomics and transcriptomics NGS data using a biologist friendly web-based graphical user interface without the need for programming expertize. BioWardrobe can be installed on consumer class hardware within an institutional local network. Analysis capabilities include predefined pipelines that allow the user to download data from either institutional core facility or public databases, perform quality control, map reads, and visualize data on a built-in mirror of the UCSC genome browser. Reads Per Kilobase per Million mapped (RPKMs) are calculated for RNA-Seq, and islands of enrichment are identified for ChIP-Seq and similar datasets. Advanced analysis capabilities include differential gene expression and binding analysis, and creation of average tag density profiles and heatmaps. BioWardrobe package and documentation is available at http://biowardrobe.com.

CD28 co-stimulatory signaling and AP-1 transcription factor are involved in the establishment of open chromatin during T cell activation

Upon encountering an antigen, naïve T helper (Th) cells are activated, differentiate into several lineages and contribute to immune response. Epigenetic changes at several cytokine and other effector genes during activation has been described previously, but the mechanism behind these changes is not understood. Here, we have profiled open chromatin during human T cell activation using ATAC-Seq approach and identified the regions where epigenetic changes take place genome-wide. Open chromatin formation correlated with induction of gene expression at nearby genes. This gene set was also enriched for functions related to T cell activation and immune response in GO analysis. T cell activation requires antigen signaling via T cell receptor and co-stimulation via CD28 which in turn result in nuclear translocation of NFAT and AP-1 transcription factors. For this reason, we profiled genome-wide distribution of these proteins during activation by ChIP-Seq. Our results show that 73% of activation-induced open chromatin loci contain binding sites for AP-1, whereas 39% of them bind NFAT. Over 90 % of the binding sites for NFAT were shared with AP-1 binding sites. We and others have previously shown that the lack of co-stimulation results in reduction or delay of AP-1 translocation into the nucleus. In order to confirm the role of AP-1 in establishment of the open chromatin structure during activation we mapped open chromatin in the cells that were activated without co-stimulation. Our data showed that ATAC signal was reduced at the open chromatin sites bound by AP-1 in the absence of co-stimulation, suggesting that CD28 signaling and AP-1 transcription factor are involved in open chromatin formation during T cell activation.

Nuclear Factor κB1/RelA Mediates Inflammation in Human Lung Epithelial Cells at Atmospheric Oxygen Levels

Oxygen levels range from 2% to 9% in vivo. Atmospheric O2 levels (21%) are known to induce cell proliferation defects and cellular senescence in primary cell cultures. However, the mechanistic basis of the deleterious effects of higher O2 levels is not fully understood. On the other hand, immortalized cells including cancer cell lines, which evade cellular senescence are normally cultured at 21% O2 and the effects of higher O2 on these cells are understudied. Here, we addressed this problem by culturing immortalized human bronchial epithelial (BEAS-2B) cells at ambient atmospheric, 21% O2 and lower, 10% O2. Our results show increased inflammatory response at 21% O2 but not at 10% O2. We found higher RelA binding at the NF-κB1/RelA target gene promoters as well as upregulation of several pro-inflammatory cytokines in cells cultured at 21% O2. RelA knockdown prevented the upregulation of the pro-inflammatory cytokines at 21% O2, suggesting NF-κB1/RelA as a major mediator of inflammatory response in cells cultured at 21% O2. Interestingly, unlike the 21% O2 cultured cells, exposure of 10% O2 cultured cells to H2O2 did not elicit inflammatory response, suggesting increased ability to tolerate oxidative stress in cells cultured at lower O2 levels.

Transcription Factor Repertoire of Homeostatic Eosinophilopoiesis

The production of mature eosinophils (Eos) is a tightly orchestrated process with the aim to sustain normal Eos levels in tissues while also maintaining low numbers of these complex and sensitive cells in the blood. To identify regulators of homeostatic eosinophilopoiesis in mice, we took a global approach to identify genome-wide transcriptome and epigenome changes that occur during homeostasis at critical developmental stages, including Eos-lineage commitment and lineage maturation. Our analyses revealed a markedly greater number of transcriptome alterations associated with Eos maturation (1199 genes) than with Eos-lineage commitment (490 genes), highlighting the greater transcriptional investment necessary for differentiation. Eos-lineage-committed progenitors (EoPs) were noted to express high levels of granule proteins and contain granules with an ultrastructure distinct from that of mature resting Eos. Our analyses also delineated a 976-gene Eos-lineage transcriptome that included a repertoire of 56 transcription factors, many of which have never previously been associated with Eos. EoPs and Eos, but not granulocyte-monocyte progenitors or neutrophils, expressed Helios and Aiolos, members of the Ikaros family of transcription factors, which regulate gene expression via modulation of chromatin structure and DNA accessibility. Epigenetic studies revealed a distinct distribution of active chromatin marks between genes induced with lineage commitment and genes induced with cell maturation during Eos development. In addition, Aiolos and Helios binding sites were significantly enriched in genes expressed by EoPs and Eos with active chromatin, highlighting a potential novel role for Helios and Aiolos in regulating gene expression during Eos development.

BioWardrobe: an integrated platform for analysis of epigenomics and transcriptomics data

High-throughput sequencing has revolutionized biology by enhancing our ability to perform genome-wide studies. However, due to lack of bioinformatics expertise, modern technologies are still beyond the capabilities of many laboratories. Herein, we present the BioWardrobe platform, which allows users to store, visualize and analyze epigenomics and transcriptomics data using a biologist-friendly web interface, without the need for programming expertise. Predefined pipelines allow users to download data, visualize results on a genome browser, calculate RPKMs (reads per kilobase per million) and identify peaks. Advanced capabilities include differential gene expression and binding analysis, and creation of average tag -density profiles and heatmaps. BioWardrobe can be found at http://biowardrobe.com.

Poised chromatin and bivalent domains facilitate the mitosis-to-meiosis transition in the male germline

Background

The male germline transcriptome changes dramatically during the mitosis-to-meiosis transition to activate late spermatogenesis genes and to transiently suppress genes commonly expressed in somatic lineages and spermatogenesis progenitor cells, termed somatic/progenitor genes.

Results

These changes reflect epigenetic regulation. Induction of late spermatogenesis genes during spermatogenesis is facilitated by poised chromatin established in the stem cell phases of spermatogonia, whereas silencing of somatic/progenitor genes during meiosis and postmeiosis is associated with formation of bivalent domains which also allows the recovery of the somatic/progenitor program after fertilization. Importantly, during spermatogenesis mechanisms of epigenetic regulation on sex chromosomes are different from autosomes: X-linked somatic/progenitor genes are suppressed by meiotic sex chromosome inactivation without deposition of H3K27me3.

Conclusions

Our results suggest that bivalent H3K27me3 and H3K4me2/3 domains are not limited to developmental promoters (which maintain bivalent domains that are silent throughout the reproductive cycle), but also underlie reversible silencing of somatic/progenitor genes during the mitosis-to-meiosis transition in late spermatogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0159-8) contains supplementary material, which is available to authorized users.

Neurotrophic tyrosine kinase receptor 1 is a direct transcriptional and epigenetic target of IL-13 involved in allergic inflammation

Although interleukin (IL)-13 and neurotrophins are functionally important for the pathogenesis of immune responses, the interaction of these pathways has not been explored. Herein, by interrogating IL-13-induced responses in human epithelial cells we show that neurotrophic tyrosine kinase receptor, type 1 (NTRK1), a cognate, high-affinity receptor for nerve growth factor (NGF), is an early transcriptional IL-13 target. Induction of NTRK1 was accompanied by accumulation of activating epigenetic marks in the promoter; transcriptional and epigenetic changes were signal transducer and activator of transcription 6 dependent. Using eosinophilic esophagitis as a model for human allergic inflammation, we found that NTRK1 was increased in inflamed tissue and dynamically expressed as a function of disease activity and that the downstream mediator of NTRK1 signaling early growth response 1 protein was elevated in allergic inflammatory tissue compared with control tissue. Unlike NTRK1, its ligand NGF was constitutively expressed in control and disease states, indicating that IL-13-stimulated NTRK1 induction is a limiting factor in pathway activation. In epithelial cells, NGF and IL-13 synergistically induced several target genes, including chemokine (C-C motif) ligand 26 (eotaxin-3). In summary, we have demonstrated that IL-13 confers epithelial cell responsiveness to NGF by regulating NTRK1 levels by a transcriptional and epigenetic mechanism and that this process likely contributes to allergic inflammation.

Functional characterization of human T cell hyporesponsiveness induced by CTLA4-Ig

During activation, T cells integrate multiple signals from APCs and cytokine milieu. The blockade of these signals can have clinical benefits as exemplified by CTLA4-Ig, which blocks interaction of B7 co-stimulatory molecules on APCs with CD28 on T cells. Variants of CTLA4-Ig, abatacept and belatacept are FDA approved as immunosuppressive agents in arthritis and transplantation, yet murine studies suggested that CTLA4-Ig could be beneficial in a number of other diseases. However, detailed analysis of human CD4 cell hyporesponsivness induced by CTLA4-Ig has not been performed. Herein, we established a model to study the effect of CTLA4-Ig on the activation of human naïve T cells in a human mixed lymphocytes system. Comparison of human CD4 cells activated in the presence or absence of CTLA4-Ig showed that co-stimulation blockade during TCR activation does not affect NFAT signaling but results in decreased activation of NF-κB and AP-1 transcription factors followed by a profound decrease in proliferation and cytokine production. The resulting T cells become hyporesponsive to secondary activation and, although capable of receiving TCR signals, fail to proliferate or produce cytokines, demonstrating properties of anergic cells. However, unlike some models of T cell anergy, these cells did not possess increased levels of the TCR signaling inhibitor CBLB. Rather, the CTLA4-Ig-induced hyporesponsiveness was associated with an elevated level of p27kip1 cyclin-dependent kinase inhibitor.

SCML2 establishes the male germline epigenome through regulation of histone H2A ubiquitination

Gametogenesis is dependent on the expression of germline-specific genes. However, it remains unknown how the germline epigenome is distinctly established from that of somatic lineages. Here we show that genes commonly expressed in somatic lineages and spermatogenesis-progenitor cells undergo repression in a genome-wide manner in late stages of the male germline and identify underlying mechanisms. SCML2, a germline-specific subunit of a Polycomb repressive complex 1 (PRC1), establishes the unique epigenome of the male germline through two distinct antithetical mechanisms. SCML2 works with PRC1 and promotes RNF2-dependent ubiquitination of H2A, thereby marking somatic/progenitor genes on autosomes for repression. Paradoxically, SCML2 also prevents RNF2-dependent ubiquitination of H2A on sex chromosomes during meiosis, thereby enabling unique epigenetic programming of sex chromosomes for male reproduction. Our results reveal divergent mechanisms involving a shared regulator by which the male germline epigenome is distinguished from that of the soma and progenitor cells.