Assay For Transposase-accessible Chromatin Sequencing Research Articles

Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets.

In castration-resistant prostate cancer (CRPC), the loss of androgen receptor (AR) dependence leads to clinically aggressive tumors with few therapeutic options. We used ATAC-seq (assay for transposase-accessible chromatin sequencing), RNA-seq, and DNA sequencing to investigate 22 organoids, six patient-derived xenografts, and 12 cell lines. We identified the well-characterized AR-dependent and neuroendocrine subtypes, as well as two AR-negative/low groups: a Wnt-dependent subtype, and a stem cell-like (SCL) subtype driven by activator protein-1 (AP-1) transcription factors. We used transcriptomic signatures to classify 366 patients, which showed that SCL is the second most common subtype of CRPC after AR-dependent. Our data suggest that AP-1 interacts with the YAP/TAZ and TEAD proteins to maintain subtype-specific chromatin accessibility and transcriptomic landscapes in this group. Together, this molecular classification reveals drug targets and can potentially guide therapeutic decisions.

Single-Cell RNA and ATAC Sequencing Reveal Hemodialysis-Related Immune Dysregulation of Circulating Immune Cell Subpopulations.

BackgroundAn increased risk of infection, malignancy, and cardiovascular diseases in maintenance hemodialysis patients is associated with hemodialysis-related immunity disturbances. Although defects in T-lymphocyte-dependent immune responses and preactivation of antigen-presenting cells have been documented in hemodialysis patients, the effects of long-term hemodialysis on the transcriptional program and chromosomal accessibility of circulating immune cell subpopulations remain poorly defined.MethodsWe integrated single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to characterize the transcriptome profiles of peripheral mononuclear cells (PBMCs) from healthy controls and maintenance hemodialysis patients. Validation of differentially expressed genes in CD4+ T cells and monocytes were performed by magnetic bead separation and quantitative real-time PCR.ResultsWe identified 16 and 15 PBMC subgroups in scRNA-seq and scATAC-seq datasets, respectively. Hemodialysis significantly suppressed the expression levels of T cell receptor (TCR) genes in CD4+ T cell subsets (e.g., TRAV4, CD45, CD3G, CD3D, CD3E) and major histocompatibility complex II (MHC-II) pathway-related genes in monocytes (HLA-DRB1, HLA-DQA2, HLA-DQA1, HLA-DPB1). Downstream pathways of TCR signaling, including PI3K-Akt-mTOR, MAPK, TNF, and NF-κB pathways, were also inhibited in CD4+ T cell subpopulations during the hemodialysis procedure. Hemodialysis altered cellular communication patterns between PBMC subgroups, particularly TGF-TGFBR, HVEM-BTLA, and IL16-CD4 signalings between CD4+ T cells and monocytes. Additionally, we found that hemodialysis inhibited the expression of AP-1 family transcription factors (JUN, JUND, FOS, FOSB) by interfering with the chromatin accessibility profile.ConclusionsOur study provides a valuable framework for future investigations of hemodialysis-related immune dysregulation and identifies potential therapeutic targets for reconstituting the circulating immune system in maintenance hemodialysis patients.

Twist2-driven chromatin remodeling governs the postnatal maturation of dermal fibroblasts.

Dermal fibroblasts lose stem cell potency after birth, which prevents regenerative healing. However, the underlying intracellular mechanisms are largely unknown. We uncover the postnatal maturation of papillary fibroblasts (PFs) driven by the extensive Twist2-mediated remodeling of chromatin accessibility. A loss of the regenerative ability of postnatal PFs occurs with decreased H3K27ac levels. Single-cell transcriptomics, assay for transposase-accessible chromatin sequencing (ATAC-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) reveal the postnatal maturation trajectory associated with the loss of the regenerative trajectory in PFs, which is characterized by a marked decrease in chromatin accessibility and H3K27ac modifications. Histone deacetylase inhibition delays spontaneous chromatin remodeling, thus maintaining the regenerative ability of postnatal PFs. Genomic analysis identifies Twist2 as a major regulator within chromatin regions with decreased accessibility during the postnatal period. When Twist2 is genetically deleted in dermal fibroblasts, the intracellular cascade of postnatal maturation is significantly delayed. Our findings reveal the comprehensive intracellular mechanisms underlying intrinsic postnatal changes in dermal fibroblasts.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Mesenchymal Stem/Stromal Cells Derived from Human Placenta.

Mesenchymal stem/stromal cells derived from placenta (PMSCs) are an attractive source for regenerative medicine because of their multidifferentiation potential and immunomodulatory capabilities. However, the cellular and molecular heterogeneity of PMSCs has not been fully characterized. Here, we applied single-cell RNA sequencing (scRNA-seq) and assay for transposase-accessible chromatin sequencing (scATAC-seq) techniques to cultured PMSCs from human full-term placenta. Based on the inferred characteristics of cell clusters, we identify several distinct subsets of PMSCs with specific characteristics, including immunomodulatory-potential and highly proliferative cell states. Furthermore, integrative analysis of gene expression and chromatin accessibility showed a clearer chromatin accessibility signature than those at the transcriptional level on immunomodulatory-related genes. Cell cycle gene-related heterogeneity can be more easily distinguished at the transcriptional than the chromatin accessibility level in PMSCs. We further reveal putative subset-specific cis-regulatory elements regulating the expression of immunomodulatory- and proliferation-related genes in the immunomodulatory-potential and proliferative subpopulations, respectively. Moreover, we infer a novel transcription factor PRDM1, which might play a crucial role in maintaining immunomodulatory capability by activating PRDM1-regulon loop. Collectively, our study first provides a comprehensive and integrative view of the transcriptomic and epigenomic features of PMSCs, which paves the way for a deeper understanding of cellular heterogeneity and offers fundamental biological insight of PMSC subset-based cell therapy.

Transcriptional states and chromatin accessibility during bovine myoblasts proliferation and myogenic differentiation.

ObjectivesAlthough major advances have been made in bovine epigenome study, the epigenetic basis for fetal skeletal muscle development still remains poorly understood. The aim is to recapitulated the time course of fetal skeletal muscle development in vitro, and explore the dynamic changes of chromatin accessibility and gene expression during bovine myoblasts proliferation and differentiation.MethodsPDGFR‐ cells were isolated from bovine fetal skeletal muscle, then cultured and induced myogenic differentiation in vitro in a time‐course study (P, D0, D2,and D4). The assay for transposase‐accessible chromatin sequencing (ATAC‐seq) and RNA sequencing (RNA‐seq) were performed.ResultsAmong the enriched transcriptional factors with high variability, we determined the effects of MAFF, ZNF384, and KLF6 in myogenesis using RNA interference (RNAi). In addition, we identified both stage‐specific genes and chromatin accessibility regions to reveal the sequential order of gene expression, transcriptional regulatory, and signal pathways involved in bovine skeletal muscle development. Further investigation integrating chromatin accessibility and transcriptome data was conducted to explore cis‐regulatory regions in line with gene expression. Moreover, we combined bovine GWAS results of growth traits with regulatory regions defined by chromatin accessibility, providing a suggestive means for a more precise annotation of genetic variants of bovine growth traits.ConclusionOverall, these findings provide valuable information for understanding the stepwise regulatory mechanisms in skeletal muscle development and conducting beef cattle genetic improvement programs.

A combinatorial indexing strategy for low-cost epigenomic profiling of plant single cells

Understanding how cis-regulatory elements facilitate gene expression is a key question in biology. Recent advances in single-cell genomics have led to the discovery of cell-specific chromatin landscapes that underlie transcription programs in animal models. However, the high equipment and reagent costs of commercial systems limit their applications for many laboratories. In this study, we developed a combinatorial index and dual PCR barcode strategy to profile the Arabidopsis thaliana root single-cell epigenome without any specialized equipment. We generated chromatin accessibility profiles for 13 576 root nuclei with an average of 12 784 unique Tn5 integrations per cell. Integration of the single-cell assay for transposase-accessible chromatin sequencing and RNA sequencing data sets enabled the identification of 24 cell clusters with unique transcription, chromatin, and cis-regulatory signatures. Comparison with single-cell data generated using the commercial microfluidic platform from 10X Genomics revealed that this low-cost combinatorial index method is capable of unbiased identification of cell-type-specific chromatin accessibility. We anticipate that, by removing cost, instrumentation, and other technical obstacles, this method will be a valuable tool for routine investigation of single-cell epigenomes and provide new insights into plant growth and development and plant interactions with the environment.

Cell-specific cis-regulatory elements and mechanisms of non-coding genetic disease in human retina and retinal organoids.

Cis-regulatory elements (CREs) play a critical role in the development and disease-states of all human cell types. In the retina, CREs have been implicated in several inherited disorders. To better characterize human retinal CREs, we performed single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) and single-nucleus RNA sequencing (snRNA-seq) on the developing and adult human retina and on induced pluripotent stem cell (iPSC)-derived retinal organoids. These analyses identified developmentally dynamic, cell-class-specific CREs, enriched transcription-factor-binding motifs, and putative target genes. CREs in the retina and organoids are highly correlated at the single-cell level, and this supports the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we disrupted a disease-associated CRE at 5q14.3, confirming its principal target gene as the miR-9-2 primary transcript and demonstrating its role in neurogenesis and gene regulation in mature glia. This study provides a resource for characterizing human retinal CREs and showcases organoids as a model to study the function of CREs that influence development and disease.

DCas9 fusion to computer-designed PRC2 inhibitor reveals functional TATA box in distal promoter region.

SUMMARYBifurcation of cellular fates, a critical process in development, requires histone 3 lysine 27 methylation (H3K27me3) marks propagated by the polycomb repressive complex 2 (PRC2). However, precise chromatin loci of functional H3K27me3 marks are not yet known. Here, we identify critical PRC2 functional sites at high resolution. We fused a computationally designed protein, EED binder (EB), which competes with EZH2 and thereby inhibits PRC2 function, to dCas9 (EBdCas9) to allow for PRC2 inhibition at a precise locus using gRNA. Targeting EBdCas9 to four different genes (TBX18, p16, CDX2, and GATA3) results in precise H3K27me3 and EZH2 reduction, gene activation, and functional outcomes in the cell cycle (p16) or trophoblast transdifferentiation (CDX2 and GATA3). In the case of TBX18, we identify a PRC2-controlled, functional TATA box >500 bp upstream of the TBX18 transcription start site (TSS) using EBdCas9. Deletion of this TATA box eliminates EBdCas9-dependent TATA binding protein (TBP) recruitment and transcriptional activation. EBdCas9 technology may provide a broadly applicable tool for epigenomic control of gene regulation.

Cytotoxicity-Related Gene Expression and Chromatin Accessibility Define a Subset of CD4+ T Cells That Mark Progression to Type 1 Diabetes.

Type 1 diabetes in children is heralded by a preclinical phase defined by circulating autoantibodies to pancreatic islet antigens. How islet autoimmunity is initiated and then progresses to clinical diabetes remains poorly understood. Only one study has reported gene expression in specific immune cells of children at risk associated with progression to islet autoimmunity. We analyzed gene expression with RNA sequencing in CD4+ and CD8+ T cells, natural killer (NK) cells, and B cells, and chromatin accessibility by assay for transposase-accessible chromatin sequencing (ATAC-seq) in CD4+ T cells, in five genetically at risk children with islet autoantibodies who progressed to diabetes over a median of 3 years ("progressors") compared with five children matched for sex, age, and HLA-DR who had not progressed ("nonprogressors"). In progressors, differentially expressed genes (DEGs) were largely confined to CD4+ T cells and enriched for cytotoxicity-related genes/pathways. Several top-ranked DEGs were validated in a semi-independent cohort of 13 progressors and 11 nonprogressors. Flow cytometry confirmed that progression was associated with expansion of CD4+ cells with a cytotoxic phenotype. By ATAC-seq, progression was associated with reconfiguration of regulatory chromatin regions in CD4+ cells, some linked to differentially expressed cytotoxicity-related genes. Our findings suggest that cytotoxic CD4+ T cells play a role in promoting progression to type 1 diabetes.

Genome-wide DNA methylation analysis in pediatric acute myeloid leukemia

AbstractWe investigated genome-wide DNA methylation patterns in 64 pediatric patients with acute myeloid leukemia (AML). Based on unsupervised clustering with the 567 most variably methylated cytosine guanine dinucleotide (CpG) sites, patients were categorized into 4 clusters associated with genetic alterations. Clusters 1 and 3 were characterized by the presence of known favorable prognostic factors, such as RUNX1-RUNX1T1 fusion and KMT2A rearrangement with low MECOM expression, and biallelic CEBPA mutations (all 8 patients), respectively. Clusters 2 and 4 comprised patients exhibiting molecular features associated with adverse outcomes, namely internal tandem duplication of FLT3 (FLT3-ITD), partial tandem duplication of KMT2A, and high PRDM16 expression. Depending on the methylation values of the 1243 CpG sites that were significantly different between FLT3-ITD+ and FLT3-ITD− AML, patients were categorized into 3 clusters: A, B, and C. The STAT5-binding motif was most frequently found close to the 1243 CpG sites. All 8 patients with FLT3-ITD in cluster A harbored high PRDM16 expression and experienced adverse events, whereas only 1 of 7 patients with FLT3-ITD in the other clusters experienced adverse events. PRDM16 expression levels were also related to DNA methylation patterns, which were drastically changed at the cutoff value of PRDM16/ABL1 = 0.10. The assay for transposase-accessible chromatin sequencing of AMLs supported enhanced chromatin accessibility around genomic regions, such as HOXB cluster genes, SCHIP1, and PRDM16, which were associated with DNA methylation changes in AMLs with FLT3-ITD and high PRDM16 expression. Our results suggest that DNA methylation levels at specific CpG sites are useful to support genetic alterations and gene expression patterns of patients with pediatric AML.

Baf45a Mediated Chromatin Remodeling Promotes Transcriptional Activation for Osteogenesis and Odontogenesis

Chromatin remodeling, specifically the tissue-specific regulation in mineralized tissues, is an understudied avenue of gene regulation. Here we show that Baf45a and Baf45d, two Baf45 homologs belong to ATPase-dependent SWI/SNF chromatin remodeling complex, preferentially expressed in osteoblasts and odontoblasts compared to Baf45b and Baf45c. Recently, biochemical studies revealed that BAF45A associates with Polybromo-associated BAF (PBAF) complex. However, the BAF45D subunit belongs to the polymorphic canonical BRG1-associated factor (cBAF) complex. Protein profiles of osteoblast and odontoblast differentiation uncovered a significant increase of BAF45A and PBAF subunits during early osteoblast and odontoblast maturation. Chromatin immunoprecipitation sequencing (ChIP-seq) during the bone marrow stromal cells (BMSCs) differentiation showed higher histone H3K9 and H3K27 acetylation modifications in the promoter of Baf45a and Baf45d and increased binding of bone and tooth specific transcription factor RUNX2. Overexpression of Baf45a in osteoblasts activates genes essential for the progression of osteoblast maturation and mineralization. Furthermore, shRNA-mediated knockdown of Baf45a in odontoblasts leads to markedly altered genes responsible for the proliferation, apoptosis, DNA repair, and modest decrease in dentinogenic marker gene expression. Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) assay in Baf45a knockout osteoblasts revealed a noticeable reduction in chromatin accessibility of osteoblast and odontoblast specific genes, along with transcription factor Atf4 and Klf4. Craniofacial mesenchyme-specific loss of Baf45a modestly reduced the mineralization of the tooth and mandibular bone. These findings indicated that BAF45A-dependent mineralized tissue-specific chromatin remodeling through PBAF-RUNX2 crosstalk results in transcriptional activation is critical for early differentiation and matrix maturation of mineralized tissues.

Single-cell ATAC-seq of fetal human retina and stem-cell-derived retinal organoids shows changing chromatin landscapes during cell fate acquisition.

We previously used single-cell transcriptomic analysis to characterize human fetal retinal development and assessed the degree to which retinal organoids recapitulate normal development. We now extend the transcriptomic analyses to incorporate single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), a powerful method used to characterize potential gene regulatory networks through the changes in accessible chromatin that accompany cell-state changes. The combination of scATAC-seq and single-cell RNA sequencing (scRNA-seq) provides a view of developing human retina at an unprecedented resolution. We identify key transcription factors relevant to specific fates and the order of the transcription factor cascades that define each of the major retinal cell types. The changing chromatin landscape is largely recapitulated in retinal organoids; however, there are differences in Notch signaling and amacrine cell gene regulation. The datasets we generated constitute an excellent resource for the continued improvement of retinal organoid technology and have the potential to inform and accelerate regenerative medicine approaches to retinal diseases.

An Improved EMSA-based Method to Prioritize Candidate cis-REs for Further Functional Validation.

Cells are the complex product of gene expression programs that involve the coordinated transcription of thousands of genes controlled by cis-regulatory elements (cis-REs). Therefore, identification of cis-REs is the key to decipher the mechanisms underlying the regulation of gene expression. Here, we describe a simple and time-effective protocol of fine-mapping cis-REs by using an electrophoresis mobility shift assay (EMSA)-based, in vitro, high-throughput (HTP) technique called regulatory element-sequencing (Reel-seq). Reel-seq can be applied to identify cis-REs at a high resolution and sensitivity over large genome regions, in a systematic and continuous manner. It can be used to prioritize candidate cis-REs as a complement to the existing approaches, such as massive parallel reporter assay (MPRA), chromatin immunoprecipitation DNA-sequencing (ChIP-seq), and the assay for transposase-accessible chromatin sequencing (ATAC-seq). Graphical abstract: Generation of the Reel-seq Library 1 and 2 (A) and identification of cis-REs by an electrophoresis mobility shift assay (EMSA)-based Reel-seq screen (B). NE: nuclear extract; NGS: next generation sequencing.

Parallel Single-Cell Multiomics Analysis of Neonatal Skin Reveals the Transitional Fibroblast States that Restrict Differentiation into Distinct Fates

One of the keys to achieving skin regeneration lies within understanding the heterogeneity of neonatal fibroblasts, which support skin regeneration. However, the molecular underpinnings regulating the cellular states and fates of these cells are not fully understood. To investigate this, we performed a parallel multiomics analysis by processing neonatal murine skin for single-cell Assay for Transposase-Accessible Chromatin sequencing and single-cell RNA sequencing separately. Our approach revealed that fibroblast clusters could be sorted into papillary and reticular lineages on the basis of transcriptome profiling, as previously reported. However, single-cell Assay for Transposase-Accessible Chromatin sequencing analysis of neonatal fibroblast lineage markers, such as Dpp4/Cd26, Corin, and Dlk1 along with markers of myofibroblasts, revealed accessible chromatin in all fibroblast populations despite their lineage-specific transcriptome profiles. These results suggest that accessible chromatin does not always translate to gene expression and that many fibroblast lineage markers reflect a fibroblast state, which includes neonatal papillary fibroblasts, reticular fibroblasts, and myofibroblasts. This analysis also provides a possible explanation as to why these marker genes can be promiscuously expressed in different fibroblast populations under different conditions. Our single-cell Assay for Transposase-Accessible Chromatin sequencing analysis also revealed that the functional lineage restriction between dermal papilla and adipocyte fates is regulated by distinct chromatin landscapes. Finally, we have developed a webtool for our multiomics analysis: https://skinregeneration.org/scatacseq-and-scrnaseq-data-from-thompson-et-al-2021-2/.

Single-Cell RNA Sequencing and Assay for Transposase-Accessible Chromatin Using Sequencing Reveals Cellular and Molecular Dynamics of Aortic Aging in Mice.

The impact of vascular aging on cardiovascular diseases has been extensively studied; however, little is known regarding the cellular and molecular mechanisms underlying age-related vascular aging in aortic cellular subpopulations. Approach and Results: Transcriptomes and transposase-accessible chromatin profiles from the aortas of 4-, 26-, and 86-week-old C57/BL6J mice were analyzed using single-cell RNA sequencing and assay for transposase-accessible chromatin sequencing. By integrating the heterogeneous transcriptome and chromatin accessibility data, we identified cell-specific TF (transcription factor) regulatory networks and open chromatin states. We also determined that aortic aging affects cell interactions, inflammation, cell type composition, dysregulation of transcriptional control, and chromatin accessibility. Endothelial cells 1 have higher gene set activity related to cellular senescence and aging than do endothelial cells 2. Moreover, construction of senescence trajectories shows that endothelial cell 1 and fibroblast senescence is associated with distinct TF open chromatin states and an mRNA expression model. Our data provide a system-wide model for transcriptional and epigenetic regulation during aortic aging at single-cell resolution.

Chromatin accessibility analysis identifies the transcription factor ETV5 as a suppressor of adipose tissue macrophage activation in obesity

Activation of adipose tissue macrophages (ATMs) contributes to chronic inflammation and insulin resistance in obesity. However, the transcriptional regulatory machinery involved in ATM activation during the development of obesity is not fully understood. Here, we profiled the chromatin accessibility of blood monocytes and ATMs from obese and lean mice using assay for transposase-accessible chromatin sequencing (ATAC-seq). We found that monocytes and ATMs from obese and lean mice exhibited distinct chromatin accessibility status. There are distinct regulatory elements that are specifically associated with monocyte or ATM activation in obesity. We also discovered several transcription factors that may regulate monocyte and ATM activation in obese mice, specifically a predicted transcription factor named ETS translocation variant 5 (ETV5). The expression of ETV5 was significantly decreased in ATMs from obese mice and its downregulation was mediated by palmitate stimulation. The decrease in ETV5 expression resulted in macrophage activation. Our results also indicate that ETV5 suppresses endoplasmic reticulum (ER) stress and Il6 expression in macrophages. Our work delineates the changes in chromatin accessibility in monocytes and ATMs during obesity, and identifies ETV5 as a critical transcription factor suppressing ATM activation, suggesting its potential use as a therapeutic target in obesity-related chronic inflammation.

Advanced Genomics-Based Approaches for Defining Allograft Rejection With Single Cell Resolution.

Solid organ transplant recipients require long-term immunosuppression for prevention of rejection. Calcineurin inhibitor (CNI)-based immunosuppressive regimens have remained the primary means for immunosuppression for four decades now, yet little is known about their effects on graft resident and infiltrating immune cell populations. Similarly, the understanding of rejection biology under specific types of immunosuppression remains to be defined. Furthermore, development of innovative, rationally designed targeted therapeutics for mitigating or preventing rejection requires a fundamental understanding of the immunobiology that underlies the rejection process. The established use of microarray technologies in transplantation has provided great insight into gene transcripts associated with allograft rejection but does not characterize rejection on a single cell level. Therefore, the development of novel genomics tools, such as single cell sequencing techniques, combined with powerful bioinformatics approaches, has enabled characterization of immune processes at the single cell level. This can provide profound insights into the rejection process, including identification of resident and infiltrating cell transcriptomes, cell-cell interactions, and T cell receptor α/β repertoires. In this review, we discuss genomic analysis techniques, including microarray, bulk RNAseq (bulkSeq), single-cell RNAseq (scRNAseq), and spatial transcriptomic (ST) techniques, including considerations of their benefits and limitations. Further, other techniques, such as chromatin analysis via assay for transposase-accessible chromatin sequencing (ATACseq), bioinformatic regulatory network analyses, and protein-based approaches are also examined. Application of these tools will play a crucial role in redefining transplant rejection with single cell resolution and likely aid in the development of future immunomodulatory therapies in solid organ transplantation.

Deep characterization of paired chromatin and transcriptomes in four immune cell types from multiple sclerosis patients.

Background: The putative involvement of chromatin states in multiple sclerosis (MS)is thus far unclear. Here we determined the association of chromatin-accessibility with concurrent genetic, epigeneticand transcriptional events. Material & methods: We generated paired assay for transposase-accessible chromatin sequencing and RNA-sequencing profiles from sorted blood immune CD4+ and CD8+ T cells, CD14+ monocytesand CD19+ B cells from healthy controls (HCs) and MS patients. Results: We identified differentially accessible regions between MS patients and HCs, primarily in CD4+ and CD19+. CD4+ regions were enriched for MS-associated single nucleotide polymorphisms and differentially methylated loci. In the vicinity of differentially accessible regions of CD4+ cells, 42 differentially expressed genes were identified. The top two dysregulated genes identified in this multilayer analysis were CCDC114 and SERTAD1. Conclusion: These findings provide new insight into the primary role of CD4+ and CD19+cells in MS.

ShinyArchR.UiO: user-friendly,integrative and open-source tool for visualization of single-cell ATAC-seq data using ArchR.

ABSTRACT MotivationMapping of chromatin accessibility landscapes in single-cells and the integration with gene expression enables a better understanding of gene regulatory mechanisms defining cell identities and cell-fate determination in development and disease. Generally, raw data generated from single-cell Assay for Transposase-Accessible Chromatin sequencing (scATAC-seq) are deposited in repositories that are generally inaccessible due to lack of in-depth knowledge of computational programming.ResultsWe have developed ShinyArchR.UiO, an R-based shiny app, that facilitates scATAC-seq data accessibility and visualization in a user-friendly, interactive and open-source web interface. ShinyArchR.UiO is an application that can streamline collaborative efforts for interpretation of massive chromatin accessibility datasets and allow for open access data sharing for wider audiences.Availability and implementation https://Github.com/EskelandLab/ShinyArchRUiO and a demo server with a hematopoietic tutorial dataset https://cancell.medisin.uio.no/ShinyArchR.UiOSupplementary information Supplementary data are available at Bioinformatics online.

Deconvolution of multiplexed transcriptional responses to wood smoke particles defines rapid aryl hydrocarbon receptor signaling dynamics

The heterogeneity of respirable particulates and compounds complicates our understanding of transcriptional responses to air pollution. Here, we address this by applying precision nuclear run-on sequencing and the assay for transposase-accessible chromatin sequencing to measure nascent transcription and chromatin accessibility in airway epithelial cells after wood smoke particle (WSP) exposure. We used transcription factor enrichment analysis to identify temporally distinct roles for ternary response factor–serum response factor complexes, the aryl hydrocarbon receptor (AHR), and NFκB in regulating transcriptional changes induced by WSP. Transcription of canonical targets of the AHR, such as CYP1A1 and AHRR, was robustly increased after just 30 min of WSP exposure, and we discovered novel AHR-regulated pathways and targets including the DNA methyltransferase, DNMT3L. Transcription of these genes and associated enhancers rapidly returned to near baseline by 120 min after exposure. The kinetics of AHR- and NFκB-regulated responses to WSP were distinguishable based on the timing of both transcriptional responses and chromatin remodeling, with induction of several cytokines implicated in maintaining NFκB-mediated responses through 120 min of exposure. In aggregate, our data establish a direct and primary role for AHR in mediating airway epithelial responses to WSP and identify crosstalk between AHR and NFκB signaling in controlling proinflammatory gene expression. This work also defines an integrated genomics-based strategy for deconvoluting multiplexed transcriptional responses to heterogeneous environmental exposures.