PEPATAC: an optimized pipeline for ATAC-seq data analysis with serial alignments

Methods for constructing and evaluating consensus genomic interval sets

The amount of genomic region data continues to increase. Integrating across diverse genomic region sets requires consensus regions, which enable comparing regions across experiments, but also by necessity lose precision in region definitions. We require methods to assess this loss of precision and build optimal consensus region sets. Here, we introduce the concept of flexible intervals and propose three novel methods for building consensus region sets, or universes: a coverage cutoff method, a likelihood method, and a Hidden Markov Model. We then propose three novel measures for evaluating how well a proposed universe fits a collection of region sets: a base-level overlap score, a region boundary distance score, and a likelihood score. We apply our methods and evaluation approaches to several collections of region sets and show how these methods can be used to evaluate fit of universes and build optimal universes. We describe scenarios where the common approach of merging regions to create consensus leads to undesirable outcomes and provide principled alternatives that provide interoperability of interval data while minimizing loss of resolution.

Predictive Prioritization of Enhancers Associated with Pancreas Disease Risk

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function

Kidney failure, the decrease of kidney function below a threshold necessary to support life, is a major cause of morbidity and mortality. We performed a genome-wide association study (GWAS) of 406,504 individuals in the UK Biobank, identifying 430 loci affecting kidney function in middle-aged adults. To investigate the cell types affected by these loci, we integrated the GWAS with human kidney candidate cis-regulatory elements (cCREs) identified using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq). Overall, 56% of kidney function heritability localized to kidney tubule epithelial cCREs and an additional 7% to kidney podocyte cCREs. Thus, most heritable differences in adult kidney function are a result of altered gene expression in these two cell types. Using enhancer assays, allele-specific scATAC-seq and machine learning, we found that many kidney function variants alter tubule epithelial cCRE chromatin accessibility and function. Using CRISPRi, we determined which genes some of these cCREs regulate, implicating NDRG1, CCNB1 and STC1 in human kidney function.

Inhibition of Renin Expression Is Regulated by an Epigenetic Switch From an Active to a Poised State

BACKGROUND

Renin-expressing cells are myoendocrine cells crucial for the maintenance of homeostasis. Renin is regulated by cAMP, p300 (histone acetyltransferase p300)/CBP (CREB-binding protein), and Brd4 (bromodomain-containing protein 4) proteins and associated pathways. However, the specific regulatory changes that occur following inhibition of these pathways are not clear.

METHODS

We treated As4.1 cells (tumoral cells derived from mouse juxtaglomerular cells that constitutively express renin) with 3 inhibitors that target different factors required for renin transcription: H-89-dihydrochloride, PKA (protein kinase A) inhibitor; JQ1, Brd4 bromodomain inhibitor; and A-485, p300/CBP inhibitor. We performed assay for transposase-accessible chromatin with sequencing (ATAC-seq), single-cell RNA sequencing, cleavage under targets and tagmentation (CUT&Tag), and chromatin immunoprecipitation sequencing for H3K27ac (acetylation of lysine 27 of the histone H3 protein) and p300 binding on biological replicates of treated and control As4.1 cells.

RESULTS

In response to each inhibitor, Ren1 expression was significantly reduced and reversible upon washout. Chromatin accessibility at the Ren1 locus did not markedly change but was globally reduced at distal elements. Inhibition of PKA led to significant reductions in H3K27ac and p300 binding specifically within the Ren1 super-enhancer region. Further, we identified enriched TF (transcription factor) motifs shared across each inhibitory treatment. Finally, we identified a set of 9 genes with putative roles across each of the 3 renin regulatory pathways and observed that each displayed differentially accessible chromatin, gene expression, H3K27ac, and p300 binding at their respective loci.

CONCLUSIONS

Inhibition of renin expression in cells that constitutively synthesize and release renin is regulated by an epigenetic switch from an active to poised state associated with decreased cell-cell communication and an epithelial-mesenchymal transition. This work highlights and helps define the factors necessary for renin cells to alternate between myoendocrine and contractile phenotypes.

ATRX restricts Human Cytomegalovirus (HCMV) viral DNA replication through heterochromatinization and minimizes unpackaged viral genomes

ATRX limits the accumulation of human cytomegalovirus (HCMV) Immediate Early (IE) proteins at the start of productive, lytic infections, and thus is a part of the cell-intrinsic defenses against infecting viruses. ATRX is a chromatin remodeler and a component of a histone chaperone complex. Therefore, we hypothesized ATRX would inhibit the transcription of HCMV IE genes by increasing viral genome heterochromatinization and decreasing its accessibility. To test this hypothesis, we quantitated viral transcription and genome structure in cells replete with or depleted of ATRX. We found ATRX did indeed limit viral IE transcription, increase viral genome chromatinization, and decrease viral genome accessibility. The inhibitory effects of ATRX extended to Early (E) and Late (L) viral protein accumulation, viral DNA replication, and progeny virion output. However, we found the negative effects of ATRX on HCMV viral DNA replication were independent of its effects on viral IE and E protein accumulation but correlated with viral genome heterochromatinization. Interestingly, the increased number of viral genomes synthesized in ATRX-depleted cells were not efficiently packaged, indicating the ATRX-mediated restriction to HCMV viral DNA replication may benefit productive infection by increasing viral fitness. Our work mechanistically describes the antiviral function of ATRX and introduces a novel, pro-viral role for this protein, perhaps explaining why, unlike during infections with other herpesviruses, it is not directly targeted by a viral countermeasure in HCMV infected cells.

Cactus: A user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

The ever decreasing cost of Next-Generation Sequencing coupled with the emergence of efficient and reproducible analysis pipelines has rendered genomic methods more accessible. However, downstream analyses are basic or missing in most workflows, creating a significant barrier for non-bioinformaticians. To help close this gap, we developed Cactus, an end-to-end pipeline for analyzing ATAC-Seq and mRNA-Seq data, either separately or jointly. Its Nextflow-, container-, and virtual environment-based architecture ensures efficient and reproducible analyses. Cactus preprocesses raw reads, conducts differential analyses between conditions, and performs enrichment analyses in various databases, including DNA-binding motifs, ChIP-Seq binding sites, chromatin states, and ontologies. We demonstrate the utility of Cactus in a multi-modal and multi-species case study as well as by showcasing its unique capabilities as compared to other ATAC-Seq pipelines. In conclusion, Cactus can assist researchers in gaining comprehensive insights from chromatin accessibility and gene expression data in a quick, user-friendly, and reproducible manner.

Spinal motor neuron development and metabolism are transcriptionally regulated by Nuclear Factor IA

Neural circuits governing all motor behaviors in vertebrates rely on the proper development of motor neurons and their precise targeting of limb muscles. Transcription factors are essential for motor neuron development, regulating their specification, migration, and axonal targeting. While transcriptional regulation of the early stages of motor neuron specification is well-established, much less is known about the role of transcription factors in the later stages of maturation and terminal arborization. Defining the molecular mechanisms of these later stages is critical for elucidating how motor circuits are constructed. Here, we demonstrate that the transcription factor Nuclear Factor-IA (NFIA) is required for motor neuron positioning, axonal branching, and neuromuscular junction formation. Moreover, we find that NFIA is required for proper mitochondrial function and ATP production, providing a new and important link between transcription factors and metabolism during motor neuron development. Together, these findings underscore the critical role of NFIA in instructing the assembly of spinal circuits for movement.

Variants in tubule epithelial regulatory elements mediate most heritable differences in human kidney function

Kidney disease is highly heritable; however, the causal genetic variants, the cell types in which these variants function, and the molecular mechanisms underlying kidney disease remain largely unknown. To identify genetic loci affecting kidney function, we performed a GWAS using multiple kidney function biomarkers and identified 462 loci. To begin to investigate how these loci affect kidney function, we generated single-cell chromatin accessibility (scATAC-seq) maps of the human kidney and identified candidate cis-regulatory elements (cCREs) for kidney podocytes, tubule epithelial cells, and kidney endothelial, stromal, and immune cells. Kidney tubule epithelial cCREs explained 58% of kidney function SNP-heritability and kidney podocyte cCREs explained an additional 6.5% of SNP-heritability. In contrast, little kidney function heritability was explained by kidney endothelial, stromal, or immune cell-specific cCREs. Through functionally informed fine-mapping, we identified putative causal kidney function variants and their corresponding cCREs. Using kidney scATAC-seq data, we created a deep learning model (which we named ChromKid) to predict kidney cell type-specific chromatin accessibility from sequence. ChromKid and allele specific kidney scATAC-seq revealed that many fine-mapped kidney function variants locally change chromatin accessibility in tubule epithelial cells. Enhancer assays confirmed that fine-mapped kidney function variants alter tubule epithelial regulatory element function. To map the genes which these regulatory elements control, we used CRISPR interference (CRISPRi) to target these regulatory elements in tubule epithelial cells and assessed changes in gene expression. CRISPRi of enhancers harboring kidney function variants regulated NDRG1 and RBPMS expression. Thus, inherited differences in tubule epithelial NDRG1 and RBPMS expression may predispose to kidney disease in humans. We conclude that genetic variants affecting tubule epithelial regulatory element function account for most SNP-heritability of human kidney function. This work provides an experimental approach to identify the variants, regulatory elements, and genes involved in polygenic disease.

Vaginal microbes alter epithelial transcriptome and induce epigenomic modifications providing insight into mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health, with microbial communities dominated by Lactobacillus species considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes, including Gardnerella vaginalis, have been associated with adverse reproductive outcomes. However, how host-microbial interactions alter specific molecular pathways and impact cervical and vaginal epithelial function remains unclear. Using RNA-sequencing, we characterized the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulates genes associated with an activated innate immune response. Unexpectedly, G. vaginalis specifically induced inflammasome pathways through activation of NLRP3-mediated increases in caspase-1, IL-1β and cell death, while live L. crispatus had minimal transcriptomic changes on epithelial cells. L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells that was confirmed by ATAC-sequencing showing reduced chromatin accessibility. This study reveals new insights into host-microbe interactions in the lower reproductive tract and suggests potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Nuclear actin structure regulates chromatin accessibility

Polymerized β-actin may provide a structural basis for chromatin accessibility and actin transport into the nucleus can guide mesenchymal stem cell (MSC) differentiation. Using MSC, we show that using CK666 to inhibit Arp2/3 directed secondary actin branching results in decreased nuclear actin structure, and significantly alters chromatin access measured with ATACseq at 24 h. The ATAC-seq results due to CK666 are distinct from those caused by cytochalasin D (CytoD), which enhances nuclear actin structure. In addition, nuclear visualization shows Arp2/3 inhibition decreases pericentric H3K9me3 marks. CytoD, alternatively, induces redistribution of H3K27me3 marks centrally. Such alterations in chromatin landscape are consistent with differential gene expression associated with distinctive differentiation patterns. Further, knockdown of the non-enzymatic monomeric actin binding protein, Arp4, leads to extensive chromatin unpacking, but only a modest increase in transcription, indicating an active role for actin-Arp4 in transcription. These data indicate that dynamic actin remodeling can regulate chromatin interactions.

A human neural crest model reveals the developmental impact of neuroblastoma-associated chromosomal aberrations

Early childhood tumours arise from transformed embryonic cells, which often carry large copy number alterations (CNA). However, it remains unclear how CNAs contribute to embryonic tumourigenesis due to a lack of suitable models. Here we employ female human embryonic stem cell (hESC) differentiation and single-cell transcriptome and epigenome analysis to assess the effects of chromosome 17q/1q gains, which are prevalent in the embryonal tumour neuroblastoma (NB). We show that CNAs impair the specification of trunk neural crest (NC) cells and their sympathoadrenal derivatives, the putative cells-of-origin of NB. This effect is exacerbated upon overexpression of MYCN, whose amplification co-occurs with CNAs in NB. Moreover, CNAs potentiate the pro-tumourigenic effects of MYCN and mutant NC cells resemble NB cells in tumours. These changes correlate with a stepwise aberration of developmental transcription factor networks. Together, our results sketch a mechanistic framework for the CNA-driven initiation of embryonal tumours.

The Methyltransferases METTL7A and METTL7B Confer Resistance to Thiol-Based Histone Deacetylase Inhibitors

Histone deacetylase inhibitors (HDACi) are part of a growing class of epigenetic therapies used for the treatment of cancer. Although HDACis are effective in the treatment of T-cell lymphomas, treatment of solid tumors with this class of drugs has not been successful. Overexpression of the multidrug resistance protein P-glycoprotein (P-gp), encoded by ABCB1, is known to confer resistance to the HDACi romidepsin in vitro, yet increased ABCB1 expression has not been associated with resistance in patients, suggesting that other mechanisms of resistance arise in the clinic. To identify alternative mechanisms of resistance to romidepsin, we selected MCF-7 breast cancer cells with romidepsin in the presence of the P-gp inhibitor verapamil to reduce the likelihood of P-gp-mediated resistance. The resulting cell line, MCF-7 DpVp300, does not express P-gp and was found to be selectively resistant to romidepsin but not to other HDACis such as belinostat, panobinostat, or vorinostat. RNA-sequencing analysis revealed upregulation of the mRNA coding for the putative methyltransferase, METTL7A, whose paralog, METTL7B, was previously shown to methylate thiol groups on hydrogen sulfide and captopril. As romidepsin has a thiol as the zinc-binding moiety, we hypothesized that METTL7A could inactivate romidepsin and other thiol-based HDACis via methylation of the thiol group. We demonstrate that expression of METTL7A or METTL7B confers resistance to thiol-based HDACis and that both enzymes are capable of methylating thiol-containing HDACis. We thus propose that METTL7A and METTL7B confer resistance to thiol-based HDACis by methylating and inactivating the zinc-binding thiol.

The ncBAF Complex Regulates Transcription in AML Through H3K27ac Sensing by BRD9

The non-canonical BAF complex (ncBAF) subunit BRD9 is essential for acute myeloid leukemia (AML) cell viability but has an unclear role in leukemogenesis. Because BRD9 is required for ncBAF complex assembly through its DUF3512 domain, precise bromodomain inhibition is necessary to parse the role of BRD9 as a transcriptional regulator from that of a scaffolding protein. To understand the role of BRD9 bromodomain function in regulating AML, we selected a panel of five AML cell lines with distinct driver mutations, disease classifications, and genomic aberrations and subjected these cells to short-term BRD9 bromodomain inhibition. We examined the bromodomain-dependent growth of these cell lines, identifying a dependency in AML cell lines but not HEK293T cells. To define a mechanism through which BRD9 maintains AML cell survival, we examined nascent transcription, chromatin accessibility, and ncBAF complex binding genome-wide after bromodomain inhibition. We identified extensive regulation of transcription by BRD9 bromodomain activity, including repression of myeloid maturation factors and tumor suppressor genes, while standard AML chemotherapy targets were repressed by inhibition of the BRD9 bromodomain. BRD9 bromodomain activity maintained accessible chromatin at both gene promoters and gene-distal putative enhancer regions, in a manner that qualitatively correlated with enrichment of BRD9 binding. Furthermore, we identified reduced chromatin accessibility at GATA, ETS, and AP-1 motifs and increased chromatin accessibility at SNAIL-, HIC-, and TP53-recognized motifs after BRD9 inhibition. These data suggest a role for BRD9 in regulating AML cell differentiation through modulation of accessibility at hematopoietic transcription factor binding sites.

SIGNIFICANCE

The bromodomain-containing protein BRD9 is essential for AML cell viability, but it is unclear whether this requirement is due to the protein's role as an epigenetic reader. We inhibited this activity and identified altered gene-distal chromatin regulation and transcription consistent with a more mature myeloid cell state.

ROCCO: a robust method for detection of open chromatin via convex optimization

MOTIVATION

Analysis of open chromatin regions across multiple samples from two or more distinct conditions can determine altered gene regulatory patterns associated with biological phenotypes and complex traits. The ATAC-seq assay allows for tractable genome-wide open chromatin profiling of large numbers of samples. Stable, broadly applicable genomic annotations of open chromatin regions are not available. Thus, most studies first identify open regions using peak calling methods for each sample independently. These are then heuristically combined to obtain a consensus peak set. Reconciling sample-specific peak results post hoc from larger cohorts is particularly challenging, and informative spatial features specific to open chromatin signals are not leveraged effectively.

RESULTS

We propose a novel method, ROCCO, that determines consensus open chromatin regions across multiple samples simultaneously. ROCCO employs robust summary statistics and solves a constrained optimization problem formulated to account for both enrichment and spatial dependence of open chromatin signal data. We show this formulation admits attractive theoretical and conceptual properties as well as superior empirical performance compared to current methodology.

AVAILABILITY AND IMPLEMENTATION

Source code, documentation, and usage demos for ROCCO are available on GitHub at: https://github.com/nolan-h-hamilton/ROCCO. ROCCO can also be installed as a stand-alone binary utility using pip/PyPI.

β1 integrin signaling governs necroptosis via the chromatin-remodeling factor CHD4

Fibrosis, characterized by sustained activation of myofibroblasts and excessive extracellular matrix (ECM) deposition, is known to be associated with chronic inflammation. Receptor-interacting protein kinase 3 (RIPK3), the central kinase of necroptosis signaling, is upregulated in fibrosis and contributes to tumor necrosis factor (TNF)-mediated inflammation. In bile-duct-ligation-induced liver fibrosis, we found that myofibroblasts are the major cell type expressing RIPK3. Genetic ablation of β1 integrin, the major profibrotic ECM receptor in fibroblasts, not only abolished ECM fibrillogenesis but also blunted RIPK3 expression via a mechanism mediated by the chromatin-remodeling factor chromodomain helicase DNA-binding protein 4 (CHD4). While the function of CHD4 has been conventionally linked to the nucleosome-remodeling deacetylase (NuRD) and CHD4-ADNP-HP1(ChAHP) complexes, we found that CHD4 potently repressed a set of genes, including Ripk3, with high locus specificity but independent of either the NuRD or the ChAHP complex. Thus, our data uncover that β1 integrin intrinsically links fibrotic signaling to RIPK3-driven inflammation via a novel mode of action of CHD4.

Vaginal microbes alter epithelial transcriptomic and epigenomic modifications providing insight into the molecular mechanisms for susceptibility to adverse reproductive outcomes

The cervicovaginal microbiome is highly associated with women's health with microbial communities dominated by Lactobacillus spp. being considered optimal. Conversely, a lack of lactobacilli and a high abundance of strict and facultative anaerobes including Gardnerella vaginalis , have been associated with adverse reproductive outcomes. However, the molecular pathways modulated by microbe interactions with the cervicovaginal epithelia remain unclear. Using RNA-sequencing, we characterize the in vitro cervicovaginal epithelial transcriptional response to different vaginal bacteria and their culture supernatants. We showed that G. vaginalis upregulated genes were associated with an activated innate immune response including anti-microbial peptides and inflammasome pathways, represented by NLRP3-mediated increases in caspase-1, IL-1β and cell death. Cervicovaginal epithelial cells exposed to L. crispatus showed limited transcriptomic changes, while exposure to L. crispatus culture supernatants resulted in a shift in the epigenomic landscape of cervical epithelial cells. ATAC-sequencing confirmed epigenetic changes with reduced chromatin accessibility. This study reveals new insight into host-microbe interactions in the lower reproductive tract and suggest potential therapeutic strategies leveraging the vaginal microbiome to improve reproductive health.

Seq2science: an end-to-end workflow for functional genomics analysis

Sequencing databases contain enormous amounts of functional genomics data, making them an extensive resource for genome-scale analysis. Reanalyzing publicly available data, and integrating it with new, project-specific data sets, can be invaluable. With current technologies, genomic experiments have become feasible for virtually any species of interest. However, using and integrating this data comes with its challenges, such as standardized and reproducible analysis. Seq2science is a multi-purpose workflow that covers preprocessing, quality control, visualization, and analysis of functional genomics sequencing data. It facilitates the downloading of sequencing data from all major databases, including NCBI SRA, EBI ENA, DDBJ, GSA, and ENCODE. Furthermore, it automates the retrieval of any genome assembly available from Ensembl, NCBI, and UCSC. It has been tested on a variety of species, and includes diverse workflows such as ATAC-, RNA-, and ChIP-seq. It consists of both generic as well as advanced steps, such as differential gene expression or peak accessibility analysis and differential motif analysis. Seq2science is built on the Snakemake workflow language and thus can be run on a range of computing infrastructures. It is available at https://github.com/vanheeringen-lab/seq2science.

cisDynet: An integrated platform for modeling gene-regulatory dynamics and networks

Chromatin accessibility sequencing has been widely used for uncovering genetic regulatory mechanisms and inferring gene regulatory networks. However, effectively integrating large-scale chromatin accessibility datasets has posed a significant challenge. This is due to the lack of a comprehensive end-to-end solution, as many existing tools primarily emphasize data preprocessing and overlook downstream analyses. To bridge this gap, we have introduced cisDynet, a holistic solution that combines streamlined data preprocessing using Snakemake and R functions with advanced downstream analysis capabilities. cisDynet excels in conventional data analyses, encompassing peak statistics, peak annotation, differential analysis, motif enrichment analysis, and more. Additionally, it allows to perform sophisticated data exploration, such as tissue-specific peak identification, time course data modeling, integration of RNA-seq data to establish peak-to-gene associations, constructing regulatory networks, and conducting enrichment analysis of genome-wide association study (GWAS) variants. As a proof of concept, we applied cisDynet to reanalyze comprehensive ATAC-seq datasets across various tissues from the Encyclopedia of DNA Elements (ENCODE) project. The analysis successfully delineated tissue-specific open chromatin regions (OCRs), established connections between OCRs and target genes, and effectively linked these discoveries with 1861 GWAS variants. Furthermore, cisDynet was instrumental in dissecting the time course open chromatin data of mouse embryonic development, revealing the dynamic behavior of OCRs over developmental stages and identifying key transcription factors governing differentiation trajectories. In summary, cisDynet offers researchers a user-friendly solution that minimizes the need for extensive coding, ensures the reproducibility of results, and greatly simplifies the exploration of epigenomic data.

Dosage amplification dictates oncogenic regulation by the NKX2-1 lineage factor in lung adenocarcinoma

Amplified oncogene expression is a critical and widespread driver event in cancer, yet our understanding of how amplification-mediated elevated dosage mediates oncogenic regulation is limited. Here, we find that the most significant focal amplification event in lung adenocarcinoma (LUAD) targets a lineage super-enhancer near the NKX2-1 lineage transcription factor. The NKX2-1 super-enhancer is targeted by focal and co-amplification with NKX2-1, and activation or repression controls NKX2-1 expression. We find that NKX2-1 is a widespread dependency in LUAD cell lines, where NKX2-1 pioneers enhancer accessibility to drive a lineage addicted state in LUAD, and NKX2-1 confers persistence to EGFR inhibitors. Notably, we find that oncogenic NKX2-1 regulation requires expression above a minimum dosage threshold-NKX2-1 dosage below this threshold is insufficient for cell viability, enhancer remodeling, and TKI persistence. Our data suggest that copy-number amplification can be a gain-of-function alteration, wherein amplification elevates oncogene expression above a critical dosage required for oncogenic regulation and cancer cell survival.

Investigation of the effects of phthalates on in vitro thyroid models with RNA-Seq and ATAC-Seq

Introduction

Phthalates are a class of endocrine-disrupting chemicals that have been shown to negatively correlate with thyroid hormone serum levels in humans and to cause a state of hyperactivity in the thyroid. However, their mechanism of action is not well described at the molecular level.

Methods

We analyzed the response of mouse thyroid organoids to the exposure to a biologically relevant dose range of the phthalates bis(2-ethylhexyl) phthalate (DEHP), di-iso-decylphthalate (DIDP), di-iso-nonylphthalate (DINP), and di-n-octylphthalate (DnOP) for 24 h and simultaneously analyzed mRNA and miRNA expression via RNA sequencing. Using the expression data, we performed differential expression and gene set enrichment analysis. We also exposed the human thyroid follicular epithelial cell line Nthy-ori 3-1 to 1 µM of DEHP or DINP for 5 days and analyzed changes in chromatin accessibility via ATAC-Seq.

Results

Dose-series analysis showed how the expression of several genes increased or decreased at the highest dose tested. As expected with the low dosing scheme, the compounds induced a modest response on the transcriptome, as we identified changes in only mmu-miR-143-3p after DINP treatment and very few differentially expressed genes. No effect was observed on thyroid markers. Ing5, a component of histones H3 and H4 acetylation complexes, was consistently upregulated in three out of four conditions compared to control, and we observed a partial overlap among the genes differentially expressed by the treatments. Gene set enrichment analysis showed enrichment in the treatment samples of the fatty acid metabolism pathway and in the control of pathways related to the receptor signalling and extracellular matrix organization. ATAC-Seq analysis showed a general increase in accessibility compared to the control, however we did not identify significant changes in accessibility in the identified regions.

Discussion

With this work, we showed that despite having only a few differentially expressed genes, diverse analysis methods could be applied to retrieve relevant information on phthalates, showing the potential of in vitro thyroid-relevant systems for the analysis of endocrine disruptors.

Ribosome biogenesis disruption mediated chromatin structure changes revealed by SRAtac, a customizable end to end analysis pipeline for ATAC-seq

The nucleolus is a large nuclear body that serves as the primary site for ribosome biogenesis. Recent studies have suggested that it also plays an important role in organizing chromatin architecture. However, to establish a causal relationship between nucleolar ribosome assembly and chromatin architecture, genetic tools are required to disrupt nucleolar ribosome biogenesis. In this study, we used ATAC-seq to investigate changes in chromatin accessibility upon specific depletion of two ribosome biogenesis components, RPOA-2 and GRWD-1, in the model organism Caenorhabditis elegans. To facilitate the analysis of ATAC-seq data, we introduced two tools: SRAlign, an extensible NGS data processing workflow, and SRAtac, a customizable end-to-end ATAC-seq analysis pipeline. Our results revealed highly comparable changes in chromatin accessibility following both RPOA-2 and GRWD-1 perturbations. However, we observed a weak correlation between changes in chromatin accessibility and gene expression. While our findings corroborate the idea of a feedback mechanism between ribosomal RNA synthesis, nucleolar ribosome large subunit biogenesis, and chromatin structure during the L1 stage of C. elegans development, they also prompt questions regarding the functional impact of these alterations on gene expression.

Unravelling human hematopoietic progenitor cell diversity through association with intrinsic regulatory factors

Hematopoietic stem and progenitor cell (HSPC) transplantation is an essential therapy for hematological conditions, but finer definitions of human HSPC subsets with associated function could enable better tuning of grafts and more routine, lower-risk application. To deeply phenotype HSPCs, following a screen of 328 antigens, we quantified 41 surface proteins and functional regulators on millions of CD34+ and CD34- cells, spanning four primary human hematopoietic tissues: bone marrow, mobilized peripheral blood, cord blood, and fetal liver. We propose more granular definitions of HSPC subsets and provide new, detailed differentiation trajectories of erythroid and myeloid lineages. These aspects of our revised human hematopoietic model were validated with corresponding epigenetic analysis and in vitro clonal differentiation assays. Overall, we demonstrate the utility of using molecular regulators as surrogates for cellular identity and functional potential, providing a framework for description, prospective isolation, and cross-tissue comparison of HSPCs in humans.

FACT regulates pluripotency through proximal and distal regulation of gene expression in murine embryonic stem cells

BACKGROUND

The FACT complex is a conserved histone chaperone with critical roles in transcription and histone deposition. FACT is essential in pluripotent and cancer cells, but otherwise dispensable for most mammalian cell types. FACT deletion or inhibition can block induction of pluripotent stem cells, yet the mechanism through which FACT regulates cell fate decisions remains unclear.

RESULTS

To explore the mechanism for FACT function, we generated AID-tagged murine embryonic cell lines for FACT subunit SPT16 and paired depletion with nascent transcription and chromatin accessibility analyses. We also analyzed SPT16 occupancy using CUT&RUN and found that SPT16 localizes to both promoter and enhancer elements, with a strong overlap in binding with OCT4, SOX2, and NANOG. Over a timecourse of SPT16 depletion, nucleosomes invade new loci, including promoters, regions bound by SPT16, OCT4, SOX2, and NANOG, and TSS-distal DNaseI hypersensitive sites. Simultaneously, transcription of Pou5f1 (encoding OCT4), Sox2, Nanog, and enhancer RNAs produced from these genes' associated enhancers are downregulated.

CONCLUSIONS

We propose that FACT maintains cellular pluripotency through a precise nucleosome-based regulatory mechanism for appropriate expression of both coding and non-coding transcripts associated with pluripotency.

The pioneer factor SOX9 competes for epigenetic factors to switch stem cell fates

During development, progenitors simultaneously activate one lineage while silencing another, a feature highly regulated in adult stem cells but derailed in cancers. Equipped to bind cognate motifs in closed chromatin, pioneer factors operate at these crossroads, but how they perform fate switching remains elusive. Here we tackle this question with SOX9, a master regulator that diverts embryonic epidermal stem cells (EpdSCs) into becoming hair follicle stem cells. By engineering mice to re-activate SOX9 in adult EpdSCs, we trigger fate switching. Combining epigenetic, proteomic and functional analyses, we interrogate the ensuing chromatin and transcriptional dynamics, slowed temporally by the mature EpdSC niche microenvironment. We show that as SOX9 binds and opens key hair follicle enhancers de novo in EpdSCs, it simultaneously recruits co-factors away from epidermal enhancers, which are silenced. Unhinged from its normal regulation, sustained SOX9 subsequently activates oncogenic transcriptional regulators that chart the path to cancers typified by constitutive SOX9 expression.

Correction of transposase sequence bias in ATAC-seq data with rule ensemble modeling

Chromatin accessibility assays have revolutionized the field of transcription regulation by providing single-nucleotide resolution measurements of regulatory features such as promoters and transcription factor binding sites. ATAC-seq directly measures how well the Tn5 transposase accesses chromatinized DNA. Tn5 has a complex sequence bias that is not effectively scaled with traditional bias-correction methods. We model this complex bias using a rule ensemble machine learning approach that integrates information from many input k-mers proximal to the ATAC sequence reads. We effectively characterize and correct single-nucleotide sequence biases and regional sequence biases of the Tn5 enzyme. Correction of enzymatic sequence bias is an important step in interpreting chromatin accessibility assays that aim to infer transcription factor binding and regulatory activity of elements in the genome.

Maternal diet alters long-term innate immune cell memory in fetal and juvenile hematopoietic stem and progenitor cells in nonhuman primate offspring

Maternal overnutrition increases inflammatory and metabolic disease risk in postnatal offspring. This constitutes a major public health concern due to increasing prevalence of these diseases, yet mechanisms remain unclear. Here, using nonhuman primate models, we show that maternal Western-style diet (mWSD) exposure is associated with persistent pro-inflammatory phenotypes at the transcriptional, metabolic, and functional levels in bone marrow-derived macrophages (BMDMs) from 3-year-old juvenile offspring and in hematopoietic stem and progenitor cells (HSPCs) from fetal and juvenile bone marrow and fetal liver. mWSD exposure is also associated with increased oleic acid in fetal and juvenile bone marrow and fetal liver. Assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling of HSPCs and BMDMs from mWSD-exposed juveniles supports a model in which HSPCs transmit pro-inflammatory memory to myeloid cells beginning in utero. These findings show that maternal diet alters long-term immune cell developmental programming in HSPCs with proposed consequences for chronic diseases featuring altered immune/inflammatory activation across the lifespan.

Epiblast-like stem cells established by Wnt/β-catenin signaling manifest distinct features of formative pluripotency and germline competence

Different formative pluripotent stem cells harboring similar functional properties have been recently established to be lineage neutral and germline competent yet have distinct molecular identities. Here, we show that WNT/β-catenin signaling activation sustains transient mouse epiblast-like cells as epiblast-like stem cells (EpiLSCs). EpiLSCs display metastable formative pluripotency with bivalent cellular energy metabolism and unique transcriptomic features and chromatin accessibility. We develop single-cell stage label transfer (scSTALT) to study the formative pluripotency continuum and reveal that EpiLSCs recapitulate a unique developmental period in vivo, filling the gap of the formative pluripotency continuum between other published formative stem cells. WNT/β-catenin signaling activation counteracts differentiation effects of activin A and bFGF by preventing complete dissolution of naive pluripotency regulatory network. Moreover, EpiLSCs have direct competence toward germline specification, which is further matured by an FGF receptor inhibitor. Our EpiLSCs can serve as an in vitro model for mimicking and studying early post-implantation development and pluripotency transition.

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

Single-cell transcriptome and accessible chromatin dynamics during endocrine pancreas development

Delineating gene regulatory networks that orchestrate cell-type specification is a continuing challenge for developmental biologists. Single-cell analyses offer opportunities to address these challenges and accelerate discovery of rare cell lineage relationships and mechanisms underlying hierarchical lineage decisions. Here, we describe the molecular analysis of mouse pancreatic endocrine cell differentiation using single-cell transcriptomics, chromatin accessibility assays coupled to genetic labeling, and cytometry-based cell purification. We uncover transcription factor networks that delineate β-, α-, and δ-cell lineages. Through genomic footprint analysis, we identify transcription factor-regulatory DNA interactions governing pancreatic cell development at unprecedented resolution. Our analysis suggests that the transcription factor Neurog3 may act as a pioneer transcription factor to specify the pancreatic endocrine lineage. These findings could improve protocols to generate replacement endocrine cells from renewable sources, like stem cells, for diabetes therapy.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

HNF4A modulates glucocorticoid action in the liver

The glucocorticoid receptor (GR) is a nuclear receptor critical to the regulation of energy metabolism and inflammation. The actions of GR are dependent on cell type and context. Here, we demonstrate the role of liver lineage-determining factor hepatocyte nuclear factor 4A (HNF4A) in defining liver specificity of GR action. In mouse liver, the HNF4A motif lies adjacent to the glucocorticoid response element (GRE) at GR binding sites within regions of open chromatin. In the absence of HNF4A, the liver GR cistrome is remodeled, with loss and gain of GR recruitment evident. Loss of chromatin accessibility at HNF4A-marked sites associates with loss of GR binding at weak GRE motifs. GR binding and chromatin accessibility are gained at sites characterized by strong GRE motifs, which show GR recruitment in non-liver tissues. The functional importance of these HNF4A-regulated GR sites is indicated by an altered transcriptional response to glucocorticoid treatment in the Hnf4a-null liver.

BET Protein Inhibition Regulates Macrophage Chromatin Accessibility and Microbiota-Dependent Colitis

Introduction

In colitis, macrophage functionality is altered compared to normal homeostatic conditions. Loss of IL-10 signaling results in an inappropriate chronic inflammatory response to bacterial stimulation. It remains unknown if inhibition of bromodomain and extra-terminal domain (BET) proteins alters usage of DNA regulatory elements responsible for driving inflammatory gene expression. We determined if the BET inhibitor, (+)-JQ1, could suppress inflammatory activation of macrophages in Il10-/- mice.

Methods

We performed ATAC-seq and RNA-seq on Il10-/- bone marrow-derived macrophages (BMDMs) cultured in the presence and absence of lipopolysaccharide (LPS) with and without treatment with (+)-JQ1 and evaluated changes in chromatin accessibility and gene expression. Germ-free Il10-/- mice were treated with (+)-JQ1, colonized with fecal slurries and underwent histological and molecular evaluation 14-days post colonization.

Results

Treatment with (+)-JQ1 suppressed LPS-induced changes in chromatin at distal regulatory elements associated with inflammatory genes, particularly in regions that contain motifs for AP-1 and IRF transcription factors. This resulted in attenuation of inflammatory gene expression. Treatment with (+)-JQ1 in vivo resulted in a mild reduction in colitis severity as compared with vehicle-treated mice.

Conclusion

We identified the mechanism of action associated with a new class of compounds that may mitigate aberrant macrophage responses to bacteria in colitis.

Linking big biomedical datasets to modular analysis with Portable Encapsulated Projects

Background

Organizing and annotating biological sample data is critical in data-intensive bioinformatics. Unfortunately, metadata formats from a data provider are often incompatible with requirements of a processing tool. There is no broadly accepted standard to organize metadata across biological projects and bioinformatics tools, restricting the portability and reusability of both annotated datasets and analysis software.

Results

To address this, we present the Portable Encapsulated Project (PEP) specification, a formal specification for biological sample metadata structure. The PEP specification accommodates typical features of data-intensive bioinformatics projects with many biological samples. In addition to standardization, the PEP specification provides descriptors and modifiers for project-level and sample-level metadata, which improve portability across both computing environments and data processing tools. PEPs include a schema validator framework, allowing formal definition of required metadata attributes for data analysis broadly. We have implemented packages for reading PEPs in both Python and R to provide a language-agnostic interface for organizing project metadata.

Conclusions

The PEP specification is an important step toward unifying data annotation and processing tools in data-intensive biological research projects. Links to tools and documentation are available at http://pep.databio.org/.