CUT&RUNTools 2.0: a pipeline for single-cell and bulk-level CUT&RUN and CUT&Tag data analysis

Auricular malformations are driven by copy number variations in a hierarchical enhancer cluster and a dominant enhancer recapitulates human pathogenesis

Enhancers, through the combinatorial action of transcription factors (TFs), dictate both the spatial specificity and the levels of gene expression, and their aberrations can result in diseases. While a HMX1 downstream enhancer is associated with ear malformations, the mechanisms underlying bilateral constricted ear (BCE) remain unclear. Here, we identify a copy number variation (CNV) containing three enhancers-collectively termed the positional identity hierarchical enhancer cluster (PI-HEC)-that drives BCE by coordinately regulating HMX1 expression. Each enhancer exhibits distinct activity-location-structure features, and the dominant enhancer with high mobility group (HMG)-box combined with Coordinator and homeodomain TF motifs modulating its activity and specificity, respectively. Mouse models demonstrate that neural crest-derived fibroblasts with aberrant Hmx1 expression in the basal pinna, along with ectopic distal pinna expression, disrupt outer ear development, affecting cartilage, muscle, and epidermis. Our findings elucidate mammalian ear morphogenesis and underscore the complexity of synergistic regulation among enhancers and between enhancers and transcription factors.

C/EBPα activates Irf8 expression in myeloid progenitors at the +56-kb enhancer to initiate cDC1 development

Development of type 1 conventional dendritic cells (cDC1s) underlies the capacity to generate antiviral and antitumor immune responses. Here, we identify the basis for cDC1 development from its earliest progenitors, determining the hierarchy of several required transcription factors and uncovering a series of mandatory cis interactions between constituent enhancers within the Irf8 superenhancer. We produced in vivo mutations of two C/EBPα binding sites that comprise the Irf8 +56-kilobase (kb) enhancer that markedly reduced IRF8 expression in all myeloid progenitors and impaired cDC1 development. These sites did not bind RUNX1 or RUNX3, and C/EBPα expression was instead regulated by their action at the Cebpa +37-kb enhancer, placing RUNX factors upstream of Cebpa in regulating Irf8. Last, we demonstrate that cis interactions between the +56-kb Irf8 enhancer and the previously reported +41- and +32-kb Irf8 enhancers are mandatory in the sequential progression of these stage-specific constituent elements.

IL-10 targets IRF transcription factors to suppress IFN and inflammatory response genes by epigenetic mechanisms

Interleukin-10 (IL-10) is pivotal in suppressing innate immune activation, in large part by suppressing induction of inflammatory genes. Despite decades of research, the molecular mechanisms underlying this inhibition have not been resolved. Here we utilized an integrated epigenomic analysis to investigate IL-10-mediated suppression of LPS and TNF responses in primary human monocytes. Instead of inhibiting core TLR4-activated pathways such as NF-κB, MAPK-AP-1 and TBK1-IRF3 signaling, IL-10 targeted IRF transcription factor activity and DNA binding, particularly IRF5 and an IRF1-mediated amplification loop. This resulted in suppression of inflammatory NF-κB target genes and near-complete suppression of interferon-stimulated genes. Mechanisms of gene inhibition included downregulation of chromatin accessibility, de novo enhancer formation and IRF1-associated H3K27ac activating histone marks. These results provide a mechanism by which IL-10 suppresses inflammatory NF-κB target genes, highlight the role of IRF1 in inflammatory gene expression and describe the suppression of IFN responses by epigenetic mechanisms.

EpiMapper: A new tool for analyzing high-throughput sequencing from CUT&Tag

Since the invention of next-generation sequencing, new methods have been developed to understand the regulation of gene expression through epigenetic markers. Among these, CUT&Tag (Cleavage Under Targets and Tagmentation) analysis has emerged as an efficient epigenomic profiling technique with low input requirements, high sensitivity, and low background signals. Although wet-lab techniques are available, data analysis remains challenging for scientists without expert-level computational skills. Therefore, we developed EpiMapper, a new Python package that simplifies the data analysis of CUT&Tag sequencing and similar techniques, such as ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) or ChIP-seq (chromatin immunoprecipitation [ChIP] assays with sequencing), and allows biomedical scientists to easily interpret the results. This new package includes every necessary step, from quality control to annotation and differential peak analysis. In particular, EpiMapper has improved functionality (e.g., reproducibility assessment) compared to previous analysis protocols and visualization plots and provides new features, such as genome annotation and differential peak analysis. Using three case studies, two on CUT&Tag and one on ATAC-seq data, the EpiMapper Python package successfully reproduced previous results.

Activation of the imprinted Prader-Willi syndrome locus by CRISPR-based epigenome editing

Epigenome editing with DNA-targeting technologies such as CRISPR-dCas9 can be used to dissect gene regulatory mechanisms and potentially treat associated disorders. For example, Prader-Willi syndrome (PWS) results from loss of paternally expressed imprinted genes on chromosome 15q11.2-q13.3, although the maternal allele is intact but epigenetically silenced. Using CRISPR repression and activation screens in human induced pluripotent stem cells (iPSCs), we identified genomic elements that control the expression of the PWS gene SNRPN from the paternal and maternal chromosomes. We showed that either targeted transcriptional activation or DNA demethylation can activate the silenced maternal SNRPN and downstream PWS transcripts. However, these two approaches function at unique regions, preferentially activating different transcript variants and involving distinct epigenetic reprogramming mechanisms. Remarkably, transient expression of the targeted demethylase leads to stable, long-term maternal SNRPN expression in PWS iPSCs. This work uncovers targeted epigenetic manipulations to reprogram a disease-associated imprinted locus and suggests possible therapeutic interventions.

Orthopedia regulates melanocortin 4 receptor transcription and energy homeostasis

Disruption of hypothalamic melanocortin 4 receptors (MC4Rs) causes obesity in mice and humans. Here, we investigated the transcriptional regulation of MC4R in the hypothalamus. In mice, we show that the homeodomain transcription factor Orthopedia (OTP) is enriched in MC4R neurons in the paraventricular nucleus (PVN) of the hypothalamus and directly regulates Mc4r transcription. Deletion of Otp in PVN neurons during development or adulthood reduced Mc4r expression, causing increased food intake and obesity. In humans, four of the five carriers of rare predicted functional OTP variants in UK Biobank had obesity. To explore a causal role for human OTP variants, we generated mice with a loss-of-function OTP mutation identified in a child with severe obesity. Heterozygous knock-in mice exhibited hyperphagia and obesity, reversed by treatment with an MC4R agonist. Our findings demonstrate that OTP regulates mammalian energy homeostasis and enable the diagnosis and treatment of individuals with obesity due to OTP deficiency.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programmes and repeats in pluripotent cells

H3K9me3 heterochromatin, established by lysine methyltransferases (KMTs) and compacted by heterochromatin protein 1 (HP1) isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3 heterochromatin stability is presently limited to individual domains and DNA repeats. Here we engineered Suv39h2-knockout mouse embryonic stem cells to degrade remaining two H3K9me3 KMTs within 1 hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A 'binary switch' governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMT depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening and exit from pluripotency within 12 h. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3 heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Integrated multi-omics analysis of zinc-finger proteins uncovers roles in RNA regulation

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation.

Hyperactivating EZH2 to augment H3K27me3 levels in regulatory T cells enhances immune suppression by driving early effector differentiation

The immunosuppressive function of regulatory T (Treg) cells is essential for maintaining immune homeostasis. Enhancer of zeste homolog 2 (EZH2), a histone H3 lysine 27 (H3K27) methyltransferase, plays a key role in maintaining Treg cell function upon CD28 co-stimulation, and Ezh2 deletion in Treg cells causes autoimmunity. Here, we assess whether increasing H3K27me3 levels, by using an Ezh2Y641F gain-of-function mutation, will improve Treg cell function. We find that Treg cells expressing Ezh2Y641F display an effector Treg phenotype, are poised for improved homing to organ tissues, and can accelerate remission from autoimmunity. The H3K27me3 landscape and transcriptome of naive Ezh2Y641F Treg cells exhibit a redistribution of H3K27me3 modifications that recapitulates the gene expression profile of activated Ezh2WT Treg cells after CD28 co-stimulation. Altogether, increased H3K27me3 levels promote the differentiation of effector Treg cells that can better suppress autoimmunity.

Distinct H3K9me3 heterochromatin maintenance dynamics govern different gene programs and repeats in pluripotent cells

H3K9me3-heterochromatin, established by lysine methyltransferases (KMTs) and compacted by HP1 isoforms, represses alternative lineage genes and DNA repeats. Our understanding of H3K9me3-heterochromatin stability is presently limited to individual domains and DNA repeats. We engineered Suv39h2 KO mouse embryonic stem cells to degrade remaining two H3K9me3-KMTs within one hour and found that both passive dilution and active removal contribute to H3K9me3 decay within 12-24 hours. We discovered four different H3K9me3 decay rates across the genome and chromatin features and transcription factor binding patterns that predict the stability classes. A "binary switch" governs heterochromatin compaction, with HP1 rapidly dissociating from heterochromatin upon KMTs' depletion and a particular threshold level of HP1 limiting pioneer factor binding, chromatin opening, and exit from pluripotency within 12 hr. Unexpectedly, receding H3K9me3 domains unearth residual HP1β peaks enriched with heterochromatin-inducing proteins. Our findings reveal distinct H3K9me3-heterochromatin maintenance dynamics governing gene networks and repeats that together safeguard pluripotency.

Review and Evaluate the Bioinformatics Analysis Strategies of ATAC-seq and CUT&Tag Data

Efficient and reliable profiling methods are essential to study epigenetics. Tn5, one of the first identified prokaryotic transposases with high DNA-binding and tagmentation efficiency, is widely adopted in different genomic and epigenomic protocols for high-throughputly exploring the genome and epigenome. Based on Tn5, the Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and the Cleavage Under Targets and Tagmentation (CUT&Tag) were developed to measure chromatin accessibility and detect DNA-protein interactions. These methodologies can be applied to large amounts of biological samples with low-input levels, such as rare tissues, embryos, and sorted single cells. However, fast and proper processing of these epigenomic data has become a bottleneck because massive data production continues to increase quickly. Furthermore, inappropriate data analysis can generate biased or misleading conclusions. Therefore, it is essential to evaluate the performance of Tn5-based ATAC-seq and CUT&Tag data processing bioinformatics tools, many of which were developed mostly for analyzing chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here, we conducted a comprehensive benchmarking analysis to evaluate the performance of eight popular software for processing ATAC-seq and CUT&Tag data. We compared the sensitivity, specificity, and peak width distribution for both narrow-type and broad-type peak calling. We also tested the influence of the availability of control IgG input in CUT&Tag data analysis. Finally, we evaluated the differential analysis strategies commonly used for analyzing the CUT&Tag data. Our study provided comprehensive guidance for selecting bioinformatics tools and recommended analysis strategies, which were implemented into Docker/Singularity images for streamlined data analysis.

Single-cell transcriptomics reveals tumor microenvironment remodeling in hepatocellular carcinoma with varying tumor subclonal complexity

Introduction

The complexity of tumor cell subclonal structure has been extensively investigated in hepatocellular carcinoma. However, the role of subclonal complexity in reshaping the tumor microenvironment (TME) remains poorly understood.

Methods

We integrated single-cell transcriptome sequencing data from four independent HCC cohorts, involving 30 samples, to decode the associations between tumor subclonal complexity and the TME. We proposed a robust metric to accurately quantify the degree of subclonal complexity for each sample based on discrete copy number variations (CNVs) profiles.

Results

We found that tumor cells in the high-complexity group originated from the cell lineage with FGB overexpression and exhibited high levels of transcription factors associated with poor survival. In contrast, tumor cells in low-complexity patients showed activation of more hallmark signaling pathways, more active cell-cell communications within the TME and a higher immune activation status. Additionally, cytokines signaling activity analysis suggested a link between HMGB1 expressed by a specific endothelial subtype and T cell proliferation.

Discussion

Our study sheds light on the intricate relationship between the complexity of subclonal structure and the TME, offering novel insights into potential therapeutic targets for HCC.

An autoimmune transcriptional circuit drives FOXP3+ regulatory T cell dysfunction

Autoimmune diseases, among the most common disorders of young adults, are mediated by genetic and environmental factors. Although CD4+FOXP3+ regulatory T cells (Tregs) play a central role in preventing autoimmunity, the molecular mechanism underlying their dysfunction is unknown. Here, we performed comprehensive transcriptomic and epigenomic profiling of Tregs in the autoimmune disease multiple sclerosis (MS) to identify critical transcriptional programs regulating human autoimmunity. We found that up-regulation of a primate-specific short isoform of PR domain zinc finger protein 1 (PRDM1-S) induces expression of serum and glucocorticoid-regulated kinase 1 (SGK1) independent from the evolutionarily conserved long PRDM1, which led to destabilization of forkhead box P3 (FOXP3) and Treg dysfunction. This aberrant PRDM1-S/SGK1 axis is shared among other autoimmune diseases. Furthermore, the chromatin landscape profiling in Tregs from individuals with MS revealed enriched activating protein-1 (AP-1)/interferon regulatory factor (IRF) transcription factor binding as candidate upstream regulators of PRDM1-S expression and Treg dysfunction. Our study uncovers a mechanistic model where the evolutionary emergence of PRDM1-S and epigenetic priming of AP-1/IRF may be key drivers of dysfunctional Tregs in autoimmune diseases.

Innate immune memory after brain injury drives inflammatory cardiac dysfunction

The medical burden of stroke extends beyond the brain injury itself and is largely determined by chronic comorbidities that develop secondarily. We hypothesized that these comorbidities might share a common immunological cause, yet chronic effects post-stroke on systemic immunity are underexplored. Here, we identify myeloid innate immune memory as a cause of remote organ dysfunction after stroke. Single-cell sequencing revealed persistent pro-inflammatory changes in monocytes/macrophages in multiple organs up to 3 months after brain injury, notably in the heart, leading to cardiac fibrosis and dysfunction in both mice and stroke patients. IL-1β was identified as a key driver of epigenetic changes in innate immune memory. These changes could be transplanted to naive mice, inducing cardiac dysfunction. By neutralizing post-stroke IL-1β or blocking pro-inflammatory monocyte trafficking with a CCR2/5 inhibitor, we prevented post-stroke cardiac dysfunction. Such immune-targeted therapies could potentially prevent various IL-1β-mediated comorbidities, offering a framework for secondary prevention immunotherapy.

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Direct neuronal reprogramming of mouse astrocytes is associated with multiscale epigenome remodeling and requires Yy1

Direct neuronal reprogramming is a promising approach to regenerate neurons from local glial cells. However, mechanisms of epigenome remodeling and co-factors facilitating this process are unclear. In this study, we combined single-cell multiomics with genome-wide profiling of three-dimensional nuclear architecture and DNA methylation in mouse astrocyte-to-neuron reprogramming mediated by Neurogenin2 (Ngn2) and its phosphorylation-resistant form (PmutNgn2), respectively. We show that Ngn2 drives multilayered chromatin remodeling at dynamic enhancer-gene interaction sites. PmutNgn2 leads to higher reprogramming efficiency and enhances epigenetic remodeling associated with neuronal maturation. However, the differences in binding sites or downstream gene activation cannot fully explain this effect. Instead, we identified Yy1, a transcriptional co-factor recruited by direct interaction with Ngn2 to its target sites. Upon deletion of Yy1, activation of neuronal enhancers, genes and ultimately reprogramming are impaired without affecting Ngn2 binding. Thus, our work highlights the key role of interactors of proneural factors in direct neuronal reprogramming.

Nuclear ATP-citrate lyase regulates chromatin-dependent activation and maintenance of the myofibroblast gene program

Differentiation of cardiac fibroblasts to myofibroblasts is necessary for matrix remodeling and fibrosis in heart failure. We previously reported that mitochondrial calcium signaling drives α-ketoglutarate-dependent histone demethylation, promoting myofibroblast formation. Here we investigate the role of ATP-citrate lyase (ACLY), a key enzyme for acetyl-CoA biosynthesis, in histone acetylation regulating myofibroblast fate and persistence in cardiac fibrosis. We show that inactivation of ACLY prevents myofibroblast differentiation and reverses myofibroblasts towards quiescence. Genetic deletion of Acly in post-activated myofibroblasts prevents fibrosis and preserves cardiac function in pressure-overload heart failure. TGFβ stimulation enhances ACLY nuclear localization and ACLY-SMAD2/3 interaction, and increases H3K27ac at fibrotic gene loci. Pharmacological inhibition of ACLY or forced nuclear expression of a dominant-negative ACLY mutant prevents myofibroblast formation and H3K27ac. Our data indicate that nuclear ACLY activity is necessary for myofibroblast differentiation and persistence by maintaining histone acetylation at TGFβ-induced myofibroblast genes. These findings provide targets to prevent and reverse pathological fibrosis.

greenPipes: an integrated data analysis pipeline for greenCUT&RUN and CUT&RUN genome-localization datasets

MOTIVATION

To study gene regulation through transcription factors and chromatin modifiers, a variety of genome-wide techniques are used. Recently, CUT&RUN-based technologies have become popular, but a pipeline for the comprehensive analysis of CUT&RUN datasets is currently lacking. Here, we present the "greenPipes" package, which includes fine-tuned parameters specifically for bioinformatic analyses of greenCUT&RUN and CUT&RUN datasets. greenPipes provides additional functionalities for data analysis and data integration with other -omics technologies, which are either not available in other pipelines developed for CUT&RUN datasets or scattered in the literature as individual packages.

AVAILABILITY AND IMPLEMENTATION

Source code and a manual of the greenPipes are freely available on GitHub website (https://github.com/snizam001/greenPipes). The test datasets, comprehensive annotation files, and other datasets are available at https://osf.io/ruhj9/.

CONTACT

n.sheikh@dkfz-heidelberg.de or m.timmers@dkfz-heidelberg.de.

SUPPLEMENTARY INFORMATION

The handbook of greenPipes is available online at Bioinformatics as Supplementary text.

Increased EZH2 function in regulatory T cells promotes their capacity to suppress autoimmunity by driving effector differentiation prior to activation

The immunosuppressive function of regulatory T (Treg) cells is essential for maintaining immune homeostasis. Enhancer of zeste homolog 2 (EZH2), a histone H3 lysine 27 (H3K27) methyltransferase, plays a key role in maintaining Treg cell function upon CD28 co-stimulation, and Ezh2 deletion in Treg cells causes autoimmunity. Here we assessed whether increased EZH2 activity in Treg cells would improve Treg cell function. Using an Ezh2 gain-of-function mutation, Ezh2 Y641F , we found that Treg cells expressing Ezh2 Y641F displayed an increased effector Treg phenotype and were poised for improved homing to organ tissues. Expression of Ezh2 Y641F in Treg cells led to more rapid remission from autoimmunity. H3K27me3 profiling and transcriptomic analysis revealed a redistribution of H3K27me3, which prompted a gene expression profile in naïve Ezh2 Y641F Treg cells that recapitulated aspects of CD28-activated Ezh2 WT Treg cells. Altogether, increased EZH2 activity promotes the differentiation of effector Treg cells that can better suppress autoimmunity.

Highlights

EZH2 function promotes effector differentiation of Treg cells.EZH2 function promotes Treg cell migration to organ tissues.EZH2 function in Treg cells improves remission from autoimmunity.EZH2 function poises naïve Treg cells to adopt a CD28-activated phenotype.

A phenotypic screening platform for identifying chemical modulators of astrocyte reactivity

Disease, injury and aging induce pathological reactive astrocyte states that contribute to neurodegeneration. Modulating reactive astrocytes therefore represent an attractive therapeutic strategy. Here we describe the development of an astrocyte phenotypic screening platform for identifying chemical modulators of astrocyte reactivity. Leveraging this platform for chemical screening, we identify histone deacetylase 3 (HDAC3) inhibitors as effective suppressors of pathological astrocyte reactivity. We demonstrate that HDAC3 inhibition reduces molecular and functional characteristics of reactive astrocytes in vitro. Transcriptional and chromatin mapping studies show that HDAC3 inhibition disarms pathological astrocyte gene expression and function while promoting the expression of genes associated with beneficial astrocytes. Administration of RGFP966, a small molecule HDAC3 inhibitor, blocks reactive astrocyte formation and promotes neuroprotection in vivo in mice. Collectively, these results establish a platform for discovering modulators of reactive astrocyte states, inform the mechanisms that control astrocyte reactivity and demonstrate the therapeutic benefits of modulating astrocyte reactivity for neurodegenerative diseases.

Immune-related transcriptomic and epigenetic reconfiguration in BV2 cells after lipopolysaccharide exposure: an in vitro omics integrative study

BACKGROUND

Molecular alterations affecting microglia have been consistently associated with the inflammatory response. These cells can have pro- or anti-inflammatory activity, phenotypes thought to be regulated by epigenetic mechanisms. Still, little is known about the details on how epigenetic marks regulate the expression of genes in the context of an inflammatory response.

METHODS

Through CUT&RUN, we profiled four genome-wide histone marks (HM) (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in lipopolysaccharide-exposed cells and compared their distributions to control cells. Transcriptomic profiles were determined through RNA-seq and differentially expressed genes were identified and contrasted with the epigenetic landscapes. Other downstream analyses were also included in this study.

RESULTS

Our results illustrate an effectively induced M1 phenotype in microglial cells derived from LPS exposure. We observed differential bound regions associated with the genes classically involved in the inflammatory response in the expected direction according to each histone modification. Consistently, our transcriptomic analysis yielded a conspicuous illustration of the LPS-induced immune activity showing the up-regulation of Nf-κB-induced mRNAs (TNF-α, nfκbiz, nfκbia) and other important genes (Marco, Il-6, etc.). Furthermore, we integrated both omics profiles and identified an important reconfiguration of the genome induced by LPS. The latter was depicted by 8 different chromatin states that changed between conditions and that associated with unique clusters of differentially expressed genes, which not only represented regulatory elements, but also underlined distinct biological functions (inhibition of morphogenesis; changes in metabolism, homeostasis, and cytokine regulation; activation of the inflammatory response).

CONCLUSION

This study exhibits important differences in the distribution of histone modifications in treated and control BV2 cells, constituting an epigenetic reconfiguration that leads to the inflammatory response. Also, it highlights the importance of these marks' regulatory role in gene expression and provides possible targets for further studies in the context of inflammation.

The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

Polygenic regression uncovers trait-relevant cellular contexts through pathway activation transformation of single-cell RNA sequencing data

Advances in single-cell RNA sequencing (scRNA-seq) techniques have accelerated functional interpretation of disease-associated variants discovered from genome-wide association studies (GWASs). However, identification of trait-relevant cell populations is often impeded by inherent technical noise and high sparsity in scRNA-seq data. Here, we developed scPagwas, a computational approach that uncovers trait-relevant cellular context by integrating pathway activation transformation of scRNA-seq data and GWAS summary statistics. scPagwas effectively prioritizes trait-relevant genes, which facilitates identification of trait-relevant cell types/populations with high accuracy in extensive simulated and real datasets. Cellular-level association results identified a novel subpopulation of naive CD8+ T cells related to COVID-19 severity and oligodendrocyte progenitor cell and microglia subsets with critical pathways by which genetic variants influence Alzheimer's disease. Overall, our approach provides new insights for the discovery of trait-relevant cell types and improves the mechanistic understanding of disease variants from a pathway perspective.

Genome-wide mapping of i-motifs reveals their association with transcription regulation in live human cells

i-Motifs (iMs) are four-stranded DNA structures that form at cytosine (C)-rich sequences in acidic conditions in vitro. Their formation in cells is still under debate. We performed CUT&Tag sequencing using the anti-iM antibody iMab and showed that iMs form within the human genome in live cells. We mapped iMs in two human cell lines and recovered C-rich sequences that were confirmed to fold into iMs in vitro. We found that iMs in cells are mainly present at actively transcribing gene promoters, in open chromatin regions, they overlap with R-loops, and their abundance and distribution are specific to each cell type. iMs with both long and short C-tracts were recovered, further extending the relevance of iMs. By simultaneously mapping G-quadruplexes (G4s), which form at guanine-rich regions, and comparing the results with iMs, we proved that the two structures can form in independent regions; however, when both iMs and G4s are present in the same genomic tract, their formation is enhanced. iMs and G4s were mainly found at genes with low and high transcription rates, respectively. Our findings support the in vivo formation of iM structures and provide new insights into their interplay with G4s as new regulatory elements in the human genome.

Bioinformatics for wet-lab scientists: practical application in sequencing analysis

BACKGROUND

Genomics data is available to the scientific community after publication of research projects and can be investigated for a multitude of research questions. However, in many cases deposited data is only assessed and used for the initial publication, resulting in valuable resources not being exploited to their full depth. MAIN: A likely reason for this is that many wetlab-based researchers are not formally trained to apply bioinformatic tools and may therefore assume that they lack the necessary experience to do so themselves. In this article, we present a series of freely available, predominantly web-based platforms and bioinformatic tools that can be combined in analysis pipelines to interrogate different types of next-generation sequencing data. Additionally to the presented exemplary route, we also list a number of alternative tools that can be combined in a mix-and-match fashion. We place special emphasis on tools that can be followed and used correctly without extensive prior knowledge in programming. Such analysis pipelines can be applied to existing data downloaded from the public domain or be compared to the results of own experiments.

CONCLUSION

Integrating transcription factor binding to chromatin (ChIP-seq) with transcriptional output (RNA-seq) and chromatin accessibility (ATAC-seq) can not only assist to form a deeper understanding of the molecular interactions underlying transcriptional regulation but will also help establishing new hypotheses and pre-testing them in silico.

Low-level repressive histone marks fine-tune gene transcription in neural stem cells

Coordinated regulation of gene activity by transcriptional and translational mechanisms poise stem cells for a timely cell-state transition during differentiation. Although important for all stemness-to-differentiation transitions, mechanistic understanding of the fine-tuning of gene transcription is lacking due to the compensatory effect of translational control. We used intermediate neural progenitor (INP) identity commitment to define the mechanisms that fine-tune stemness gene transcription in fly neural stem cells (neuroblasts). We demonstrate that the transcription factor FruitlessC (FruC) binds cis-regulatory elements of most genes uniquely transcribed in neuroblasts. Loss of fruC function alone has no effect on INP commitment but drives INP dedifferentiation when translational control is reduced. FruC negatively regulates gene expression by promoting low-level enrichment of the repressive histone mark H3K27me3 in gene cis-regulatory regions. Identical to fruC loss-of-function, reducing Polycomb Repressive Complex 2 activity increases stemness gene activity. We propose low-level H3K27me3 enrichment fine-tunes gene transcription in stem cells, a mechanism likely conserved from flies to humans.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

The ETS transcription factor ETV6 constrains the transcriptional activity of EWS-FLI to promote Ewing sarcoma

Transcription factors (TFs) are frequently mutated in cancer. Paediatric cancers exhibit few mutations genome-wide but frequently harbour sentinel mutations that affect TFs, which provides a context to precisely study the transcriptional circuits that support mutant TF-driven oncogenesis. A broadly relevant mechanism that has garnered intense focus involves the ability of mutant TFs to hijack wild-type lineage-specific TFs in self-reinforcing transcriptional circuits. However, it is not known whether this specific type of circuitry is equally crucial in all mutant TF-driven cancers. Here we describe an alternative yet central transcriptional mechanism that promotes Ewing sarcoma, wherein constraint, rather than reinforcement, of the activity of the fusion TF EWS-FLI supports cancer growth. We discover that ETV6 is a crucial TF dependency that is specific to this disease because it, counter-intuitively, represses the transcriptional output of EWS-FLI. This work discovers a previously undescribed transcriptional mechanism that promotes cancer.

YAP-mediated mechanotransduction in urinary bladder remodeling: Based on RNA-seq and CUT&Tag

Yes-associated protein (YAP) is an important transcriptional coactivator binding to transcriptional factors that engage in many downstream gene transcription. Partial bladder outlet obstruction (pBOO) causes a massive burden to patients and finally leads to bladder fibrosis. Several cell types engage in the pBOO pathological process, including urothelial cells, smooth muscle cells, and fibroblasts. To clarify the function of YAP in bladder fibrosis, we performed the RNA-seq and CUT&Tag of the bladder smooth muscle cell to analyze the YAP ablation of human bladder smooth muscle cells (hBdSMCs) and immunoprecipitation of YAP. 141 differentially expressed genes (DEGs) were identified through RNA-seq between YAP-knockdown and nature control. After matching with the results of CUT&Tag, 36 genes were regulated directly by YAP. Then we identified the hub genes in the DEGs, including CDCA5, CENPA, DTL, NCAPH, and NEIL3, that contribute to cell proliferation. Thus, our study provides a regulatory network of YAP in smooth muscle proliferation. The possible effects of YAP on hBdSMC might be a vital target for pBOO-associated bladder fibrosis.

A genetic disorder reveals a hematopoietic stem cell regulatory network co-opted in leukemia

The molecular regulation of human hematopoietic stem cell (HSC) maintenance is therapeutically important, but limitations in experimental systems and interspecies variation have constrained our knowledge of this process. Here, we have studied a rare genetic disorder due to MECOM haploinsufficiency, characterized by an early-onset absence of HSCs in vivo. By generating a faithful model of this disorder in primary human HSCs and coupling functional studies with integrative single-cell genomic analyses, we uncover a key transcriptional network involving hundreds of genes that is required for HSC maintenance. Through our analyses, we nominate cooperating transcriptional regulators and identify how MECOM prevents the CTCF-dependent genome reorganization that occurs as HSCs differentiate. We show that this transcriptional network is co-opted in high-risk leukemias, thereby enabling these cancers to acquire stem cell properties. Collectively, we illuminate a regulatory network necessary for HSC self-renewal through the study of a rare experiment of nature.

Molecular mechanisms of ferroptosis and the potential therapeutic targets of ferroptosis signaling pathways for glioblastoma

Ferroptosis is a newly identified form of cell death that differs from autophagy, apoptosis and necrosis, and its molecular characteristics include iron-dependent lipid reactive oxygen species accumulation, mitochondrial morphology changes, and membrane permeability damage. These characteristics are closely related to various human diseases, especially tumors of the nervous system. Glioblastoma is the most common primary malignant tumor of the adult central nervous system, and the 5-year survival rate is only 4%-5%. This study reviewed the role and mechanism of ferroptosis in glioblastoma and the research status and progress on ferroptosis as a potential therapeutic target. The mechanism of ferroptosis is related to the intracellular iron metabolism level, lipid peroxide content and glutathione peroxidase 4 activity. It is worth exploring how ferroptosis can be applied in disease treatment; however, the relation between ferroptosis and other apoptosis methods is poorly understood and methods of applying ferroptosis to drug-resistant tumors are insufficient. Ferroptosis is a promising therapeutic target for glioblastoma. In-depth studies of its mechanism of action in glioblastoma and applications for clinical treatment are expected to provide insights for glioblastoma patients.

Ferroptosis: A new therapeutic target for bladder cancer

Bladder cancer (BC) is the most frequent type of urinary system cancer. The prognosis of BC is poor due to high metastasis rates and multidrug resistance. Hence, development of novel therapies targeting BC cell death is urgently needed. As a novel cell death type with strong antitumor potential, ferroptosis has been investigated by many groups for its potential in BC treatment. As an iron-dependent cell death process, ferroptosis is characterized by excessive oxidative phospholipids. The molecular mechanisms of ferroptosis include iron overload and the system Xc-GSH-GPX4 signaling pathway. A recent study revealed that ferroptosis is involved in the metastasis, treatment, and prognosis of BC. Herein, in this review, we comprehensively summarize the mechanism of ferroptosis, address newly identified targets involved in ferroptosis, and discuss the potential of new clinical therapies targeting ferroptosis in BC.

The therapeutic potential of natural products for treating pancreatic cancer

Pancreatic cancer is one of the most malignant tumors of the digestive tract, with the poor prognosis and low 5-year survival rate less than 10%. Although surgical resection and chemotherapy as gemcitabine (first-line treatment) has been applied to the pancreatic cancer patients, the overall survival rates of pancreatic cancer are quite low due to drug resistance. Therefore, it is of urgent need to develop alternative strategies for its treatment. In this review, we summarized the major herbal drugs and metabolites, including curcumin, triptolide, Panax Notoginseng Saponins and their metabolites etc. These compounds with antioxidant, anti-angiogenic and anti-metastatic activities can inhibit the progression and metastasis of pancreatic cancer. Expecting to provide comprehensive information of potential natural products, our review provides valuable information and strategies for pancreatic cancer treatment.

In situ tools for chromatin structural epigenomics

Technological progress over the past 15 years has fueled an explosion in genome-wide chromatin profiling tools that take advantage of low-cost short-read sequencing technologies to map particular chromatin features. Here, we survey the recent development of epigenomic tools that provide precise positions of chromatin proteins genome-wide in intact cells or nuclei. Some profiling tools are based on tethering Micrococcal Nuclease to chromatin proteins of interest in situ, whereas others similarly tether Tn5 transposase to integrate DNA sequencing adapters (tagmentation) and so eliminate the need for library preparation. These in situ cleavage and tagmentation tools have gained in popularity over the past few years, with many protocol enhancements and adaptations for single-cell and spatial chromatin profiling. The application of experimental and computational tools to address problems in gene regulation, eukaryotic development, and human disease are helping to define the emerging field of chromatin structural epigenomics.

Exploring epigenetic reprogramming during central nervous system infection

Epigenetics involves the study of various modes of adaptable transcriptional regulation, contributing to cell identity, characteristics, and function. During central nervous system (CNS) infection, epigenetic mechanisms can exert pronounced control over the maturation and antimicrobial properties of nearly every immune cell type. Epigenetics is a relatively new field, with the first mention of these marks proposed only a half-century ago and a substantial body of immunological epigenetic research emerging only in the last few decades. Here, we review the best-characterized epigenetic marks and their functions as well as illustrate how various immune cell populations responding to CNS infection utilize these marks to organize their activation state and inflammatory processes. We also discuss the metabolic and clinical implications of epigenetic marks and the rapidly expanding set of tools available to researchers that are enabling elucidation of increasingly detailed genetic regulatory pathways. These considerations paint an intricate picture of inflammatory regulation, where epigenetic marks influence genetic, signaling, and environmental elements to orchestrate a tailored immunological response to the threat at hand, cementing epigenetics as an important player in immunity.

Advancements in Genomic and Behavioral Neuroscience Analysis for the Study of Normal and Pathological Brain Function

Psychiatric and neurological disorders are influenced by an undetermined number of genes and molecular pathways that may differ among afflicted individuals. Functionally testing and characterizing biological systems is essential to discovering the interrelationship among candidate genes and understanding the neurobiology of behavior. Recent advancements in genetic, genomic, and behavioral approaches are revolutionizing modern neuroscience. Although these tools are often used separately for independent experiments, combining these areas of research will provide a viable avenue for multidimensional studies on the brain. Herein we will briefly review some of the available tools that have been developed for characterizing novel cellular and animal models of human disease. A major challenge will be openly sharing resources and datasets to effectively integrate seemingly disparate types of information and how these systems impact human disorders. However, as these emerging technologies continue to be developed and adopted by the scientific community, they will bring about unprecedented opportunities in our understanding of molecular neuroscience and behavior.

ePeak: from replicated chromatin profiling data to epigenomic dynamics

We present ePeak, a Snakemake-based pipeline for the identification and quantification of reproducible peaks from raw ChIP-seq, CUT&RUN and CUT&Tag epigenomic profiling techniques. It also includes a statistical module to perform tailored differential marking and binding analysis with state of the art methods. ePeak streamlines critical steps like the quality assessment of the immunoprecipitation, spike-in calibration and the selection of reproducible peaks between replicates for both narrow and broad peaks. It generates complete reports for data quality control assessment and optimal interpretation of the results. We advocate for a differential analysis that accounts for the biological dynamics of each chromatin factor. Thus, ePeak provides linear and nonlinear methods for normalisation as well as conservative and stringent models for variance estimation and significance testing of the observed marking/binding differences. Using a published ChIP-seq dataset, we show that distinct populations of differentially marked/bound peaks can be identified. We study their dynamics in terms of read coverage and summit position, as well as the expression of the neighbouring genes. We propose that ePeak can be used to measure the richness of the epigenomic landscape underlying a biological process by identifying diverse regulatory regimes.

Dual function NFI factors control fetal hemoglobin silencing in adult erythroid cells

The mechanisms by which the fetal-type β-globin-like genes HBG1 and HBG2 are silenced in adult erythroid precursor cells remain a fundamental question in human biology and have therapeutic relevance to sickle cell disease (SCD) and β-thalassemia. Here, we identify via a CRISPR-Cas9 genetic screen two members of the NFI transcription factor family – NFIA and NFIX – as HBG1/2 repressors. NFIA and NFIX are expressed at elevated levels in adult erythroid cells compared to fetal cells, and function cooperatively to repress HBG1/2 in cultured cells and in human-to-mouse xenotransplants. Genomic profiling, genome editing, and DNA binding assays demonstrate that the potent concerted activity of NFIA and NFIX is explained in part by their ability to stimulate the expression of BCL11A, a known silencer of the HBG1/2 genes, and in part by directly repressing the HBG1/2 genes. Thus, NFI factors emerge as versatile regulators of the fetal-to-adult switch in β-globin production.

AUTS2 Syndrome: Molecular Mechanisms and Model Systems

AUTS2 syndrome is a genetic disorder that causes intellectual disability, microcephaly, and other phenotypes. Syndrome severity is worse when mutations involve 3' regions (exons 9-19) of the AUTS2 gene. Human AUTS2 protein has two major isoforms, full-length (1259 aa) and C-terminal (711 aa), the latter produced from an alternative transcription start site in exon 9. Structurally, AUTS2 contains the putative "AUTS2 domain" (∼200 aa) conserved among AUTS2 and its ohnologs, fibrosin, and fibrosin-like-1. Also, AUTS2 contains extensive low-complexity sequences and intrinsically disordered regions, features typical of RNA-binding proteins. During development, AUTS2 is expressed by specific progenitor cell and neuron types, including pyramidal neurons and Purkinje cells. AUTS2 localizes mainly in cell nuclei, where it regulates transcription and RNA metabolism. Some studies have detected AUTS2 in neurites, where it may regulate cytoskeletal dynamics. Neurodevelopmental functions of AUTS2 have been studied in diverse model systems. In zebrafish, auts2a morphants displayed microcephaly. In mice, excision of different Auts2 exons (7, 8, or 15) caused distinct phenotypes, variously including neonatal breathing abnormalities, cerebellar hypoplasia, dentate gyrus hypoplasia, EEG abnormalities, and behavioral changes. In mouse embryonic stem cells, AUTS2 could promote or delay neuronal differentiation. Cerebral organoids, derived from an AUTS2 syndrome patient containing a pathogenic missense variant in exon 9, exhibited neocortical growth defects. Emerging technologies for analysis of human cerebral organoids will be increasingly useful for understanding mechanisms underlying AUTS2 syndrome. Questions for future research include whether AUTS2 binds RNA directly, how AUTS2 regulates neurogenesis, and how AUTS2 modulates neural circuit formation.

Variant to function mapping at single-cell resolution through network propagation

Genome-wide association studies in combination with single-cell genomic atlases can provide insights into the mechanisms of disease-causal genetic variation. However, identification of disease-relevant or trait-relevant cell types, states and trajectories is often hampered by sparsity and noise, particularly in the analysis of single-cell epigenomic data. To overcome these challenges, we present SCAVENGE, a computational algorithm that uses network propagation to map causal variants to their relevant cellular context at single-cell resolution. We demonstrate how SCAVENGE can help identify key biological mechanisms underlying human genetic variation, applying the method to blood traits at distinct stages of human hematopoiesis, to monocyte subsets that increase the risk for severe Coronavirus Disease 2019 (COVID-19) and to intermediate lymphocyte developmental states that predispose to acute leukemia. Our approach not only provides a framework for enabling variant-to-function insights at single-cell resolution but also suggests a more general strategy for maximizing the inferences that can be made using single-cell genomic data.

Characterizing the Copy Number Variation of Non-Coding RNAs Reveals Potential Therapeutic Targets and Prognostic Markers of LUSC

Lung squamous cell carcinoma (LUSC) has a poor clinical prognosis and a lack of available targeted therapies. Therefore, there is an urgent need to identify novel prognostic markers and therapeutic targets to assist in the diagnosis and treatment of LUSC. With the development of high-throughput sequencing technology, integrated analysis of multi-omics data will provide annotation of pathogenic non-coding variants and the role of non-coding sequence variants in cancers. Here, we integrated RNA-seq profiles and copy number variation (CNV) data to study the effects of non-coding variations on gene regulatory network. Furthermore, the 372 long non-coding RNAs (lncRNA) regulated by CNV were used as candidate genes, which could be used as biomarkers for clinical application. Nine lncRNAs including LINC00896, MCM8-AS1, LINC01251, LNX1-AS1, GPRC5D-AS1, CTD-2350J17.1, LINC01133, LINC01121, and AC073130.1 were recognized as prognostic markers for LUSC. By exploring the association of the prognosis-related lncRNAs (pr-lncRNAs) with immune cell infiltration, GPRC5D-AS1 and LINC01133 were highlighted as markers of the immunosuppressive microenvironment. Additionally, the cascade response of pr-lncRNA-CNV-mRNA-physiological functions was revealed. Taken together, the identification of prognostic markers and carcinogenic regulatory mechanisms will contribute to the individualized treatment for LUSC and promote the development of precision medicine.