Quantification of Differential Transcription Factor Activity and Multiomics-Based Classification into Activators and Repressors: diffTF

Machine-guided cell-fate engineering

The creation of induced pluripotent stem cells (iPSCs) has enabled scientists to explore the function, mechanisms, and differentiation processes of many types of cells. One of the fastest and most efficient approaches is transcription factor (TF) over-expression. However, finding the right combination of TFs to over-express to differentiate iPSCs directly into other cell types is a difficult task. Here, we describe a machine-learning (ML) pipeline, called CellCartographer, that uses chromatin accessibility and transcriptomics data to design multiplex TF pooled-screening experiments for cell-type conversions that then may be iteratively refined. We validate this method by differentiating iPSCs into twelve cell types at low efficiency in preliminary screens and iteratively refine our TF combinations to achieve high-efficiency differentiation for six of these cell types in <6 days. Finally, we functionally characterize iPSC-derived cytotoxic T cells (iCytoTs), regulatory T cells (iTregs), type II astrocytes (iAstIIs), and hepatocytes (iHeps) to validate functionally accurate differentiation.

SWI/SNF-type complexes-transcription factor interplay: a key regulatory interaction

ATP-dependent switch/sucrose nonfermenting-type chromatin remodeling complexes (SWI/SNF CRCs) are multiprotein machineries altering chromatin structure, thus controlling the accessibility of genomic DNA to various regulatory proteins including transcription factors (TFs). SWI/SNF CRCs are highly evolutionarily conserved among eukaryotes. There are three main subtypes of SWI/SNF CRCs: canonical (cBAF), polybromo (pBAF), and noncanonical (ncBAF) in humans and their functional Arabidopsis counterparts SYD-associated SWI/SNF (SAS), MINU-associated SWI/SNF (MAS), and BRAHMA (BRM)-associated SWI/SNF (BAS). Here, we highlight the importance of interplay between SWI/SNF CRCs and TFs in human and Arabidopsis and summarize recent advances demonstrating their role in controlling important regulatory processes. We discuss possible mechanisms involved in TFs and SWI/SNF CRCs-dependent transcriptional control of gene expression. We indicate that Arabidopsis may serve as a valuable model for the identification of evolutionarily conserved SWI/SNF-TF interactions and postulate that further exploration of the TFs and SWI/SNF CRCs-interplay, especially in the context of the role of particular SWI/SNF CRC subtypes, TF type, as well as cell/tissue and conditions, among others, will help address important questions related to the specificity of SWI/SNF-TF interactions and the sequence of events occurring on their target genes.

Active repression of cell fate plasticity by PROX1 safeguards hepatocyte identity and prevents liver tumorigenesis

Cell fate plasticity enables development, yet unlocked plasticity is a cancer hallmark. While transcription master regulators induce lineage-specific genes to restrict plasticity, it remains unclear whether plasticity is actively suppressed by lineage-specific repressors. Here we computationally predict so-called safeguard repressors for 18 cell types that block phenotypic plasticity lifelong. We validated hepatocyte-specific candidates using reprogramming, revealing that prospero homeobox protein 1 (PROX1) enhanced hepatocyte identity by direct repression of alternative fate master regulators. In mice, Prox1 was required for efficient hepatocyte regeneration after injury and was sufficient to prevent liver tumorigenesis. In line with patient data, Prox1 depletion caused hepatocyte fate loss in vivo and enabled the transition of hepatocellular carcinoma to cholangiocarcinoma. Conversely, overexpression promoted cholangiocarcinoma to hepatocellular carcinoma transdifferentiation. Our findings provide evidence for PROX1 as a hepatocyte-specific safeguard and support a model where cell-type-specific repressors actively suppress plasticity throughout life to safeguard lineage identity and thus prevent disease.

Integrative multiomics reveals common endotypes across PSEN1, PSEN2, and APP mutations in familial Alzheimer's disease

BACKGROUND

PSEN1, PSEN2, and APP mutations cause Alzheimer's disease (AD) with an early age at onset (AAO) and progressive cognitive decline. PSEN1 mutations are more common and generally have an earlier AAO; however, certain PSEN1 mutations cause a later AAO, similar to those observed in PSEN2 and APP.

METHODS

We examined whether common disease endotypes exist across these mutations with a later AAO (~ 55 years) using hiPSC-derived neurons from familial Alzheimer's disease (FAD) patients harboring mutations in PSEN1A79V, PSEN2N141I, and APPV717I and mechanistically characterized by integrating RNA-seq and ATAC-seq.

RESULTS

We identified common disease endotypes, such as dedifferentiation, dysregulation of synaptic signaling, repression of mitochondrial function and metabolism, and inflammation. We ascertained the master transcriptional regulators associated with these endotypes, including REST, ASCL1, and ZIC family members (activation), and NRF1 (repression).

CONCLUSIONS

FAD mutations share common regulatory changes within endotypes with varying severity, resulting in reversion to a less-differentiated state. The regulatory mechanisms described offer potential targets for therapeutic interventions.

ChromatinHD connects single-cell DNA accessibility and conformation to gene expression through scale-adaptive machine learning

Gene regulation is inherently multiscale, but scale-adaptive machine learning methods that fully exploit this property in single-nucleus accessibility data are still lacking. Here, we develop ChromatinHD, a pair of scale-adaptive models that uses the raw accessibility data, without peak-calling or windows, to link regions to gene expression and determine differentially accessible chromatin. We show how ChromatinHD consistently outperforms existing peak and window-based approaches and find that this is due to a large number of uniquely captured, functional accessibility changes within and outside of putative cis-regulatory regions. Furthermore, ChromatinHD can delineate collaborating regulatory regions, including their preferential genomic conformations, that drive gene expression. Finally, our models also use changes in ATAC-seq fragment lengths to identify dense binding of transcription factors, a feature not captured by footprinting methods. Altogether, ChromatinHD, available at https://chromatinhd.org , is a suite of computational tools that enables a data-driven understanding of chromatin accessibility at various scales and how it relates to gene expression.

Optimized reporters for multiplexed detection of transcription factor activity

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 62 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways. A record of this paper's transparent peer review process is included in the supplemental information.

Integration of genetic and chromatin modification data pinpoints autoimmune-specific remodeling of enhancer landscape in CD4+ T cells

CD4+ T cells play a crucial role in adaptive immune responses and have been implicated in the pathogenesis of autoimmune diseases (ADs). Despite numerous studies, the molecular mechanisms underlying T cell dysregulation in ADs remain incompletely understood. Here, we used chromatin immunoprecipitation (ChIP)-sequencing of active chromatin and transcriptomic data from CD4+ T cells of healthy donors and patients with systemic lupus erythematosus (SLE), psoriasis, juvenile idiopathic arthritis (JIA), and Graves' disease to investigate the role of enhancers in AD pathogenesis. By generating enhancer-based gene regulatory networks (eGRNs), we identified disease-specific dysregulated pathways and potential downstream target genes of enhancers harboring AD-associated single-nucleotide polymorphisms (SNPs), which we also validated using chromatin-capture (HiC) data and CRISPR interference (CRISPRi) in primary CD4+ T cells. Our results suggest that alterations in the regulatory landscapes of CD4+ T cells, including enhancers, contribute to the development of ADs and provide a basis for developing new therapeutic approaches.

Decoding mutational hotspots in human disease through the gene modules governing thymic regulatory T cells

Computational strategies to extract meaningful biological information from multiomics data are in great demand for effective clinical use, particularly in complex immune-mediated disorders. Regulatory T cells (Tregs) are essential for immune homeostasis and self-tolerance, controlling inflammatory and autoimmune processes in many diseases with a multigenic basis. Here, we quantify the Transcription Factor (TF) differential occupancy landscape to uncover the Gene Regulatory Modules governing lineage-committed Tregs in the human thymus, and show that it can be used as a tool to prioritise variants in complex diseases. We combined RNA-seq and ATAC-seq and generated a matrix of differential TF binding to genes differentially expressed in Tregs, in contrast to their counterpart conventional CD4 single-positive thymocytes. The gene loci of both established and novel genetic interactions uncovered by the Gene Regulatory Modules were significantly enriched in rare variants carried by patients with common variable immunodeficiency, here used as a model of polygenic-based disease with severe inflammatory and autoimmune manifestations. The Gene Regulatory Modules controlling the Treg signature can, therefore, be a valuable resource for variant classification, and to uncover new therapeutic targets. Overall, our strategy can also be applied in other biological processes of interest to decipher mutational hotspots in individual genomes.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs-in our case NFkB-at the protein level.

Histone H3.3 lysine 9 and 27 control repressive chromatin at cryptic enhancers and bivalent promoters

Histone modifications are associated with distinct transcriptional states, but it is unclear whether they instruct gene expression. To investigate this, we mutate histone H3.3 K9 and K27 residues in mouse embryonic stem cells (mESCs). Here, we find that H3.3K9 is essential for controlling specific distal intergenic regions and for proper H3K27me3 deposition at promoters. The H3.3K9A mutation resulted in decreased H3K9me3 at regions encompassing endogenous retroviruses and induced a gain of H3K27ac and nascent transcription. These changes in the chromatin environment unleash cryptic enhancers, resulting in the activation of distinctive transcriptional programs and culminating in protein expression normally restricted to specialized immune cell types. The H3.3K27A mutant disrupts the deposition and spreading of the repressive H3K27me3 mark, particularly impacting bivalent genes with higher basal levels of H3.3 at promoters. Therefore, H3.3K9 and K27 crucially orchestrate repressive chromatin states at cis-regulatory elements and bivalent promoters, respectively, and instruct proper transcription in mESCs.

On the identification of differentially-active transcription factors from ATAC-seq data

ATAC-seq has emerged as a rich epigenome profiling technique, and is commonly used to identify Transcription Factors (TFs) underlying given phenomena. A number of methods can be used to identify differentially-active TFs through the accessibility of their DNA-binding motif, however little is known on the best approaches for doing so. Here we benchmark several such methods using a combination of curated datasets with various forms of short-term perturbations on known TFs, as well as semi-simulations. We include both methods specifically designed for this type of data as well as some that can be repurposed for it. We also investigate variations to these methods, and identify three particularly promising approaches (a chromVAR-limma workflow with critical adjustments, monaLisa and a combination of GC smooth quantile normalization and multivariate modeling). We further investigate the specific use of nucleosome-free fragments, the combination of top methods, and the impact of technical variation. Finally, we illustrate the use of the top methods on a novel dataset to characterize the impact on DNA accessibility of TRAnscription Factor TArgeting Chimeras (TRAFTAC), which can deplete TFs - in our case NFkB - at the protein level.

Optimized reporters for multiplexed detection of transcription factor activity

In any given cell type, dozens of transcription factors (TFs) act in concert to control the activity of the genome by binding to specific DNA sequences in regulatory elements. Despite their considerable importance in determining cell identity and their pivotal role in numerous disorders, we currently lack simple tools to directly measure the activity of many TFs in parallel. Massively parallel reporter assays (MPRAs) allow the detection of TF activities in a multiplexed fashion; however, we lack basic understanding to rationally design sensitive reporters for many TFs. Here, we use an MPRA to systematically optimize transcriptional reporters for 86 TFs and evaluate the specificity of all reporters across a wide array of TF perturbation conditions. We thus identified critical TF reporter design features and obtained highly sensitive and specific reporters for 60 TFs, many of which outperform available reporters. The resulting collection of "prime" TF reporters can be used to uncover TF regulatory networks and to illuminate signaling pathways.

Basis of gene-specific transcription regulation by the Integrator complex

The Integrator complex attenuates gene expression via the premature termination of RNA polymerase II (RNAP2) at promoter-proximal pausing sites. It is required for stimulus response, cell differentiation, and neurodevelopment, but how gene-specific and adaptive regulation by Integrator is achieved remains unclear. Here, we identify two sites on human Integrator subunits 13/14 that serve as binding hubs for sequence-specific transcription factors (TFs) and other transcription effector complexes. When Integrator is attached to paused RNAP2, these hubs are positioned upstream of the transcription bubble, consistent with simultaneous TF-promoter tethering. The TFs co-localize with Integrator genome-wide, increase Integrator abundance on target genes, and co-regulate responsive transcriptional programs. For instance, sensory cilia formation induced by glucose starvation depends on Integrator-TF contacts. Our data suggest TF-mediated promoter recruitment of Integrator as a widespread mechanism for targeted transcription regulation.

Gene Regulatory Programs that Specify Age-Related Differences during Thymocyte Development

T cell development is fundamental to immune system establishment, yet how this development changes with age remains poorly understood. Here, we construct a transcriptional and epigenetic atlas of T cell developmental programs in neonatal and adult mice, revealing the ontogeny of divergent gene regulatory programs and their link to age-related differences in phenotype and function. Specifically, we identify a gene module that diverges with age from the earliest stages of genesis and includes programs that govern effector response and cell cycle regulation. Moreover, we reveal that neonates possess more accessible chromatin during early thymocyte development, likely establishing poised gene expression programs that manifest later in thymocyte development. Finally, we leverage this atlas, employing a CRISPR-based perturbation approach coupled with single-cell RNA sequencing as a readout to uncover a conserved transcriptional regulator, Zbtb20, that contributes to age-dependent differences in T cell development. Altogether, our study defines transcriptional and epigenetic programs that regulate age-specific differences in T cell development.

Individual transcription factors modulate both the micromovement of chromatin and its long-range structure

The control of eukaryotic gene expression is intimately connected to highly dynamic chromatin structures. Gene regulation relies on activator and repressor transcription factors (TFs) that induce local chromatin opening and closing. However, it is unclear how nucleus-wide chromatin organization responds dynamically to the activity of specific TFs. Here, we examined how two TFs with opposite effects on local chromatin accessibility modulate chromatin dynamics nucleus-wide. We combine high-resolution diffusion mapping and dense flow reconstruction and correlation in living cells to obtain an imaging-based, nanometer-scale analysis of local diffusion processes and long-range coordinated movements of both chromatin and TFs. We show that the expression of either an individual transcriptional activator (CDX2) or repressor (SIX6) with large numbers of binding sites increases chromatin mobility nucleus-wide, yet they induce opposite coherent chromatin motions at the micron scale. Hi-C analysis of higher-order chromatin structures shows that induction of the pioneer factor CDX2 leads both to changes in local chromatin interactions and the distribution of A and B compartments, thus relating the micromovement of chromatin with changes in compartmental structures. Given that inhibition of transcription initiation and elongation by RNA Pol II has a partial impact on the global chromatin dynamics induced by CDX2, we suggest that CDX2 overexpression alters chromatin structure dynamics both dependently and independently of transcription. Our biophysical analysis shows that sequence-specific TFs can influence chromatin structure on multiple architectural levels, arguing that local chromatin changes brought by TFs alter long-range chromatin mobility and its organization.

Pharmacological induction of chromatin remodeling drives chemosensitization in triple-negative breast cancer

Targeted therapies have improved outcomes for certain cancer subtypes, but cytotoxic chemotherapy remains a mainstay for triple-negative breast cancer (TNBC). The epithelial-to-mesenchymal transition (EMT) is a developmental program co-opted by cancer cells that promotes metastasis and chemoresistance. There are no therapeutic strategies specifically targeting mesenchymal-like cancer cells. We report that the US Food and Drug Administration (FDA)-approved chemotherapeutic eribulin induces ZEB1-SWI/SNF-directed chromatin remodeling to reverse EMT that curtails the metastatic propensity of TNBC preclinical models. Eribulin induces mesenchymal-to-epithelial transition (MET) in primary TNBC in patients, but conventional chemotherapy does not. In the treatment-naive setting, but not after acquired resistance to other agents, eribulin sensitizes TNBC cells to subsequent treatment with other chemotherapeutics. These findings provide an epigenetic mechanism of action of eribulin, supporting its use early in the disease process for MET induction to prevent metastatic progression and chemoresistance. These findings warrant prospective clinical evaluation of the chemosensitizing effects of eribulin in the treatment-naive setting.

RFX4 is an intrinsic factor for neuronal differentiation through induction of proneural genes POU3F2 and NEUROD1

Proneural genes play a crucial role in neuronal differentiation. However, our understanding of the regulatory mechanisms governing proneural genes during neuronal differentiation remains limited. RFX4, identified as a candidate regulator of proneural genes, has been reported to be associated with the development of neuropsychiatric disorders. To uncover the regulatory relationship, we utilized a combination of multi-omics data, including ATAC-seq, ChIP-seq, Hi-C, and RNA-seq, to identify RFX4 as an upstream regulator of proneural genes. We further validated the role of RFX4 using an in vitro model of neuronal differentiation with RFX4 knock-in and a CRISPR-Cas9 knock-out system. As a result, we found that RFX4 directly interacts with the promoters of POU3F2 and NEUROD1. Transcriptomic analysis revealed a set of genes associated with neuronal development, which are highly implicated in the development of neuropsychiatric disorders, including schizophrenia. Notably, ectopic expression of RFX4 can drive human embryonic stem cells toward a neuronal fate. Our results strongly indicate that RFX4 serves as a direct upstream regulator of proneural genes, a role that is essential for normal neuronal development. Impairments in RFX4 function could potentially be related to the development of various neuropsychiatric disorders. However, understanding the precise mechanisms by which the RFX4 gene influences the onset of neuropsychiatric disorders requires further investigation through human genetic studies.

Cytotoxic CD8+ Temra cells show loss of chromatin accessibility at genes associated with T cell activation

As humans age, their memory T cell compartment expands due to the lifelong exposure to antigens. This expansion is characterized by terminally differentiated CD8+ T cells (Temra), which possess NK cell-like phenotype and are associated with chronic inflammatory conditions. Temra cells are predominantly driven by the sporadic reactivation of cytomegalovirus (CMV), yet their epigenomic patterns and cellular heterogeneity remain understudied. To address this gap, we correlated their gene expression profiles with chromatin openness and conducted single-cell transcriptome analysis, comparing them to other CD8+ subsets and CMV-responses. We confirmed that Temra cells exhibit high expression of genes associated with cytotoxicity and lower expression of costimulatory and chemokine genes. The data revealed that CMV-responsive CD8+ T cells (Tcmv) were predominantly derived from a mixed population of Temra and memory cells (Tcm/em) and shared their transcriptomic profiles. Using ATAC-seq analysis, we identified 1449 differentially accessible chromatin regions between CD8+ Temra and Tcm/em cells, of which only 127 sites gained chromatin accessibility in Temra cells. We further identified 51 gene loci, including costimulatory CD27, CD28, and ICOS genes, whose chromatin accessibility correlated with their gene expression. The differential chromatin regions Tcm/em cells were enriched in motifs that bind multiple transcriptional activators, such as Jun/Fos, NFkappaB, and STAT, whereas the open regions in Temra cells mainly contained binding sites of T-box transcription factors. Our single-cell analysis of CD8+CCR7loCD45RAhi sorted Temra population showed several subsets of Temra and NKT-like cells and CMC1+ Temra populations in older individuals that were shifted towards decreased cytotoxicity. Among CD8+CCR7loCD45RAhi sorted cells, we found a decreased proportion of IL7R+ Tcm/em-like and MAIT cells in individuals with high levels of CMV antibodies (CMVhi). These results shed new light on the molecular and cellular heterogeneity of CD8+ Temra cells and their relationship to aging and CMV infection.

Molecular Atlas of HER2+ Breast Cancer Cells Treated with Endogenous Ligands: Temporal Insights into Mechanisms of Trastuzumab Resistance

Trastuzumab therapy in HER2+ breast cancer patients has mixed success owing to acquired resistance to therapy. A detailed understanding of downstream molecular cascades resulting from trastuzumab resistance is yet to emerge. In this study, we investigate the cellular mechanisms underlying acquired resistance using trastuzumab-sensitive and -resistant cancer cells (BT474 and BT474R) treated with endogenous ligands EGF and HRG across time. We probe early receptor organization through microscopy and signaling events through multiomics measurements and assess the bioenergetic state through mitochondrial measurements. Integrative analyses of our measurements reveal significant alterations in EGF-treated BT474 HER2 membrane dynamics and robust downstream activation of PI3K/AKT/mTORC1 signaling. EGF-treated BT474R shows a sustained interferon-independent activation of the IRF1/STAT1 cascade, potentially contributing to trastuzumab resistance. Both cell lines exhibit temporally divergent metabolic demands and HIF1A-mediated stress responses. BT474R demonstrates inherently increased mitochondrial activity. HRG treatment in BT474R leads to a pronounced reduction in AR expression, affecting downstream lipid metabolism with implications for treatment response. Our results provide novel insights into mechanistic changes underlying ligand treatment in BT474 and BT474R and emphasize the pivotal role of endogenous ligands. These results can serve as a framework for furthering the understanding of trastuzumab resistance, with therapeutic implications for women with acquired resistance.

Stem cell-like reprogramming is required for leukemia-initiating activity in B-ALL

B cell acute lymphoblastic leukemia (B-ALL) is a multistep disease characterized by the hierarchical acquisition of genetic alterations. However, the question of how a primary oncogene reprograms stem cell-like properties in committed B cells and leads to a preneoplastic population remains unclear. Here, we used the PAX5::ELN oncogenic model to demonstrate a causal link between the differentiation blockade, the self-renewal, and the emergence of preleukemic stem cells (pre-LSCs). We show that PAX5::ELN disrupts the differentiation of preleukemic cells by enforcing the IL7r/JAK-STAT pathway. This disruption is associated with the induction of rare and quiescent pre-LSCs that sustain the leukemia-initiating activity, as assessed using the H2B-GFP model. Integration of transcriptomic and chromatin accessibility data reveals that those quiescent pre-LSCs lose B cell identity and reactivate an immature molecular program, reminiscent of human B-ALL chemo-resistant cells. Finally, our transcriptional regulatory network reveals the transcription factor EGR1 as a strong candidate to control quiescence/resistance of PAX5::ELN pre-LSCs as well as of blasts from human B-ALL.

Single-cell spatial multi-omics and deep learning dissect enhancer-driven gene regulatory networks in liver zonation

In the mammalian liver, hepatocytes exhibit diverse metabolic and functional profiles based on their location within the liver lobule. However, it is unclear whether this spatial variation, called zonation, is governed by a well-defined gene regulatory code. Here, using a combination of single-cell multiomics, spatial omics, massively parallel reporter assays and deep learning, we mapped enhancer-gene regulatory networks across mouse liver cell types. We found that zonation affects gene expression and chromatin accessibility in hepatocytes, among other cell types. These states are driven by the repressors TCF7L1 and TBX3, alongside other core hepatocyte transcription factors, such as HNF4A, CEBPA, FOXA1 and ONECUT1. To examine the architecture of the enhancers driving these cell states, we trained a hierarchical deep learning model called DeepLiver. Our study provides a multimodal understanding of the regulatory code underlying hepatocyte identity and their zonation state that can be used to engineer enhancers with specific activity levels and zonation patterns.

Computational tools for inferring transcription factor activity

Transcription factors (TFs) are essential players in orchestrating the regulatory landscape in cells. Still, their exact modes of action and dependencies on other regulatory aspects remain elusive. Since TFs act cell type-specific and each TF has its own characteristics, untangling their regulatory interactions from an experimental point of view is laborious and convoluted. Thus, there is an ongoing development of computational tools that estimate transcription factor activity (TFA) from a variety of data modalities, either based on a mapping of TFs to their putative target genes or in a genome-wide, gene-unspecific fashion. These tools can help to gain insights into TF regulation and to prioritize candidates for experimental validation. We want to give an overview of available computational tools that estimate TFA, illustrate examples of their application, debate common result validation strategies, and discuss assumptions and concomitant limitations.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

Multiomics dynamic learning enables personalized diagnosis and prognosis for pancancer and cancer subtypes

Artificial intelligence (AI) approaches in cancer analysis typically utilize a 'one-size-fits-all' methodology characterizing average patient responses. This manner neglects the diverse conditions in the pancancer and cancer subtypes of individual patients, resulting in suboptimal outcomes in diagnosis and treatment. To overcome this limitation, we shift from a blanket application of statistics to a focus on the explicit recognition of patient-specific abnormalities. Our objective is to use multiomics data to empower clinicians with personalized molecular descriptions that allow for customized diagnosis and interventions. Here, we propose a highly trustworthy multiomics learning (HTML) framework that employs multiomics self-adaptive dynamic learning to process each sample with data-dependent architectures and computational flows, ensuring personalized and trustworthy patient-centering of cancer diagnosis and prognosis. Extensive testing on a 33-type pancancer dataset and 12 cancer subtype datasets underscored the superior performance of HTML compared with static-architecture-based methods. Our findings also highlighting the potential of HTML in elucidating complex biological pathogenesis and paving the way for improved patient-specific care in cancer treatment.

Characterization of an eye field-like state during optic vesicle organoid development

Specification of the eye field (EF) within the neural plate marks the earliest detectable stage of eye development. Experimental evidence, primarily from non-mammalian model systems, indicates that the stable formation of this group of cells requires the activation of a set of key transcription factors. This crucial event is challenging to probe in mammals and, quantitatively, little is known regarding the regulation of the transition of cells to this ocular fate. Using optic vesicle organoids to model the onset of the EF, we generate time-course transcriptomic data allowing us to identify dynamic gene expression programmes that characterize this cellular-state transition. Integrating this with chromatin accessibility data suggests a direct role of canonical EF transcription factors in regulating these gene expression changes, and highlights candidate cis-regulatory elements through which these transcription factors act. Finally, we begin to test a subset of these candidate enhancer elements, within the organoid system, by perturbing the underlying DNA sequence and measuring transcriptomic changes during EF activation.

Chromatin Landscape Governing Murine Epidermal Differentiation

Chromatin landscape and regulatory networks are determinants in lineage specification and differentiation. To define the temporospatial differentiation axis in murine epidermal cells in vivo, we generated datasets profiling expression dynamics (RNA sequencing), chromatin accessibility (assay for transposase-accessible chromatin using sequencing), architecture (Hi-C), and histone modifications (chromatin immunoprecipitation followed by sequencing) in the epidermis. We show that many differentially regulated genes are suppressed during the differentiation process, with superenhancers controlling differentiation-specific epigenomic changes. Our data shows the relevance of the Dlx/Klf/Grhl combinatorial regulatory network in maintaining correct temporospatial gene expression during epidermal differentiation. We determined differential open compartments, topologically associating domain score, and looping in the basal cell and suprabasal cell epidermal fractions, with the evolutionarily conserved epidermal differentiation complex region showing distinct suprabasal cell-specific topologically associating domain and loop formation that coincided with superenhancer sites. Overall, our study provides a global genome-wide resource of chromatin dynamics that define unrecognized regulatory networks and the epigenetic control of Dlx3-bound superenhancer elements during epidermal differentiation.

GRaNIE and GRaNPA: inference and evaluation of enhancer-mediated gene regulatory networks

Enhancers play a vital role in gene regulation and are critical in mediating the impact of noncoding genetic variants associated with complex traits. Enhancer activity is a cell-type-specific process regulated by transcription factors (TFs), epigenetic mechanisms and genetic variants. Despite the strong mechanistic link between TFs and enhancers, we currently lack a framework for jointly analysing them in cell-type-specific gene regulatory networks (GRN). Equally important, we lack an unbiased way of assessing the biological significance of inferred GRNs since no complete ground truth exists. To address these gaps, we present GRaNIE (Gene Regulatory Network Inference including Enhancers) and GRaNPA (Gene Regulatory Network Performance Analysis). GRaNIE (https://git.embl.de/grp-zaugg/GRaNIE) builds enhancer-mediated GRNs based on covariation of chromatin accessibility and RNA-seq across samples (e.g. individuals), while GRaNPA (https://git.embl.de/grp-zaugg/GRaNPA) assesses the performance of GRNs for predicting cell-type-specific differential expression. We demonstrate their power by investigating gene regulatory mechanisms underlying the response of macrophages to infection, cancer and common genetic traits including autoimmune diseases. Finally, our methods identify the TF PURA as a putative regulator of pro-inflammatory macrophage polarisation.

BORIS/CTCFL-mediated chromatin accessibility alterations promote a pro-invasive transcriptional signature in melanoma cells

Melanoma is the deadliest form of skin cancer, due to its tendency to metastasize early. Brother of regulator of imprinted sites (BORIS), also known as CCCTC binding factor-like (CTCFL), is a transcription regulator that becomes ectopically expressed in melanoma. We recently showed that BORIS contributes to melanoma phenotype switching by altering the gene expression program of melanoma cells from an intermediate melanocytic state toward a more mesenchymal-like state. However, the mechanism underlying this transcriptional switch remains unclear. Here, ATAC-seq was used to study BORIS-mediated chromatin accessibility alterations in melanoma cells harboring an intermediate melanocytic state. The gene set that gained promoter accessibility, following ectopic BORIS expression, showed enrichment for biological processes associated with melanoma invasion, while promoters of genes associated with proliferation showed reduced accessibility. Integration of ATAC-seq and RNA-seq data demonstrated that increased chromatin accessibility was associated with transcriptional upregulation of genes involved in tumor progression processes, and the aberrant activation of oncogenic transcription factors, while reduced chromatin accessibility and downregulated genes were associated with repressed activity of tumor suppressors and proliferation factors. Together, these findings indicate that BORIS mediates transcriptional reprogramming in melanoma cells by altering chromatin accessibility and gene expression, shifting the cellular transcription landscape of melanoma cells toward a mesenchymal-like genetic signature.

Pharmacological Induction of mesenchymal-epithelial transition chemosensitizes breast cancer cells and prevents metastatic progression

The epithelial-mesenchymal transition (EMT) is a developmental program co-opted by tumor cells that aids the initiation of the metastatic cascade. Tumor cells that undergo EMT are relatively chemoresistant, and there are currently no therapeutic avenues specifically targeting cells that have acquired mesenchymal traits. We show that treatment of mesenchymal-like triple-negative breast cancer (TNBC) cells with the microtubule-destabilizing chemotherapeutic eribulin, which is FDA-approved for the treatment of advanced breast cancer, leads to a mesenchymal-epithelial transition (MET). This MET is accompanied by loss of metastatic propensity and sensitization to subsequent treatment with other FDA-approved chemotherapeutics. We uncover a novel epigenetic mechanism of action that supports eribulin pretreatment as a path to MET induction that curtails metastatic progression and the evolution of therapy resistance.

Proteomic profiling reveals distinct phases to the restoration of chromatin following DNA replication

Chromatin organization must be maintained during cell proliferation to preserve cellular identity and genome integrity. However, DNA replication results in transient displacement of DNA-bound proteins, and it is unclear how they regain access to newly replicated DNA. Using quantitative proteomics coupled to Nascent Chromatin Capture or isolation of Proteins on Nascent DNA, we provide time-resolved binding kinetics for thousands of proteins behind replisomes within euchromatin and heterochromatin in human cells. This shows that most proteins regain access within minutes to newly replicated DNA. In contrast, 25% of the identified proteins do not, and this delay cannot be inferred from their known function or nuclear abundance. Instead, chromatin organization and G1 phase entry affect their reassociation. Finally, DNA replication not only disrupts but also promotes recruitment of transcription factors and chromatin remodelers, providing a significant advance in understanding how DNA replication could contribute to programmed changes of cell memory.

ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

TF-Prioritizer: a Java pipeline to prioritize condition-specific transcription factors

BACKGROUND

Eukaryotic gene expression is controlled by cis-regulatory elements (CREs), including promoters and enhancers, which are bound by transcription factors (TFs). Differential expression of TFs and their binding affinity at putative CREs determine tissue- and developmental-specific transcriptional activity. Consolidating genomic datasets can offer further insights into the accessibility of CREs, TF activity, and, thus, gene regulation. However, the integration and analysis of multimodal datasets are hampered by considerable technical challenges. While methods for highlighting differential TF activity from combined chromatin state data (e.g., chromatin immunoprecipitation [ChIP], ATAC, or DNase sequencing) and RNA sequencing data exist, they do not offer convenient usability, have limited support for large-scale data processing, and provide only minimal functionality for visually interpreting results.

RESULTS

We developed TF-Prioritizer, an automated pipeline that prioritizes condition-specific TFs from multimodal data and generates an interactive web report. We demonstrated its potential by identifying known TFs along with their target genes, as well as previously unreported TFs active in lactating mouse mammary glands. Additionally, we studied a variety of ENCODE datasets for cell lines K562 and MCF-7, including 12 histone modification ChIP sequencing as well as ATAC and DNase sequencing datasets, where we observe and discuss assay-specific differences.

CONCLUSION

TF-Prioritizer accepts ATAC, DNase, or ChIP sequencing and RNA sequencing data as input and identifies TFs with differential activity, thus offering an understanding of genome-wide gene regulation, potential pathogenesis, and therapeutic targets in biomedical research.

Pluripotency factors are repurposed to shape the epigenomic landscape of neural crest cells

Yamanaka factors are essential for establishing pluripotency in embryonic stem cells, but their function in multipotent stem cell populations is poorly understood. Here, we show that OCT4 and SOX2 cooperate with tissue-specific transcription factors to promote neural crest formation. By assessing avian and human neural crest cells at distinct developmental stages, we characterized the epigenomic changes that occur during their specification, migration, and early differentiation. This analysis determined that the OCT4-SOX2 dimer is required to establish a neural crest epigenomic signature that is lost upon cell fate commitment. The OCT4-SOX2 genomic targets in the neural crest differ from those of embryonic stem cells, indicating the dimer displays context-specific functions. Binding of OCT4-SOX2 to neural crest enhancers requires pioneer factor TFAP2A, which physically interacts with the dimer to modify its genomic targets. Our results demonstrate how Yamanaka factors are repurposed in multipotent cells to control chromatin organization and define their developmental potential.

Identifying a novel role for the master regulator Tal1 in the Endothelial to Hematopoietic Transition

Progress in the generation of Hematopoietic Stem and Progenitor Cells (HSPCs) in vitro and ex vivo has been built on the knowledge of developmental hematopoiesis, underscoring the importance of understanding this process. HSPCs emerge within the embryonic vasculature through an Endothelial-to-Hematopoietic Transition (EHT). The transcriptional regulator Tal1 exerts essential functions in the earliest stages of blood development, but is considered dispensable for the EHT. Nevertheless, Tal1 is expressed with its binding partner Lmo2 and it homologous Lyl1 in endothelial and transitioning cells at the time of EHT. Here, we investigated the function of these genes using a mouse embryonic-stem cell (mESC)-based differentiation system to model hematopoietic development. We showed for the first time that the expression of TAL1 in endothelial cells is crucial to ensure the efficiency of the EHT process and a sustained hematopoietic output. Our findings uncover an important function of Tal1 during the EHT, thus filling the current gap in the knowledge of the role of this master gene throughout the whole process of hematopoietic development.

Phenotypic heterogeneity driven by plasticity of the intermediate EMT state governs disease progression and metastasis in breast cancer

The epithelial-to-mesenchymal transition (EMT) is frequently co-opted by cancer cells to enhance migratory and invasive cell traits. It is a key contributor to heterogeneity, chemoresistance, and metastasis in many carcinoma types, where the intermediate EMT state plays a critical tumor-initiating role. We isolate multiple distinct single-cell clones from the SUM149PT human breast cell line spanning the EMT spectrum having diverse migratory, tumor-initiating, and metastatic qualities, including three unique intermediates. Using a multiomics approach, we identify CBFβ as a key regulator of metastatic ability in the intermediate state. To quantify epithelial-mesenchymal heterogeneity within tumors, we develop an advanced multiplexed immunostaining approach using SUM149-derived orthotopic tumors and find that the EMT state and epithelial-mesenchymal heterogeneity are predictive of overall survival in a cohort of stage III breast cancer. Our model reveals previously unidentified insights into the complex EMT spectrum and its regulatory networks, as well as the contributions of epithelial-mesenchymal plasticity (EMP) in tumor heterogeneity in breast cancer.

Drug-microenvironment perturbations reveal resistance mechanisms and prognostic subgroups in CLL

The tumour microenvironment and genetic alterations collectively influence drug efficacy in cancer, but current evidence is limited and systematic analyses are lacking. Using chronic lymphocytic leukaemia (CLL) as a model disease, we investigated the influence of 17 microenvironmental stimuli on 12 drugs in 192 genetically characterised patient samples. Based on microenvironmental response, we identified four subgroups with distinct clinical outcomes beyond known prognostic markers. Response to multiple microenvironmental stimuli was amplified in trisomy 12 samples. Trisomy 12 was associated with a distinct epigenetic signature. Bromodomain inhibition reversed this epigenetic profile and could be used to target microenvironmental signalling in trisomy 12 CLL. We quantified the impact of microenvironmental stimuli on drug response and their dependence on genetic alterations, identifying interleukin 4 (IL4) and Toll-like receptor (TLR) stimulation as the strongest actuators of drug resistance. IL4 and TLR signalling activity was increased in CLL-infiltrated lymph nodes compared with healthy samples. High IL4 activity correlated with faster disease progression. The publicly available dataset can facilitate the investigation of cell-extrinsic mechanisms of drug resistance and disease progression.

Ranking reprogramming factors for cell differentiation

Transcription factor over-expression is a proven method for reprogramming cells to a desired cell type for regenerative medicine and therapeutic discovery. However, a general method for the identification of reprogramming factors to create an arbitrary cell type is an open problem. Here we examine the success rate of methods and data for differentiation by testing the ability of nine computational methods (CellNet, GarNet, EBseq, AME, DREME, HOMER, KMAC, diffTF and DeepAccess) to discover and rank candidate factors for eight target cell types with known reprogramming solutions. We compare methods that use gene expression, biological networks and chromatin accessibility data, and comprehensively test parameter and preprocessing of input data to optimize performance. We find the best factor identification methods can identify an average of 50-60% of reprogramming factors within the top ten candidates, and methods that use chromatin accessibility perform the best. Among the chromatin accessibility methods, complex methods DeepAccess and diffTF have higher correlation with the ranked significance of transcription factor candidates within reprogramming protocols for differentiation. We provide evidence that AME and diffTF are optimal methods for transcription factor recovery that will allow for systematic prioritization of transcription factor candidates to aid in the design of new reprogramming protocols.

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.

Epigenetic regulation of innate immune memory in microglia

Background

Microglia are the tissue-resident macrophages of the CNS. They originate in the yolk sac, colonize the CNS during embryonic development and form a self-sustaining population with limited turnover. A consequence of their relative slow turnover is that microglia can serve as a long-term memory for inflammatory or neurodegenerative events.

Methods

Using ATAC-, ChIP- and RNA-sequencing, we characterized the epigenomes and transcriptomes of FACS-purified microglia from mice exposed to different stimuli. A repeated endotoxin challenge (LPS) was used to induce tolerance in microglia, while genotoxic stress (DNA repair deficiency-induced accelerated aging through Ercc1 deficiency) resulted in primed (hypersensitive) microglia.

Results

Whereas the enrichment of permissive epigenetic marks at enhancer regions could explain training (hyper-responsiveness) of primed microglia to an LPS challenge, the tolerized response of microglia seems to be regulated by loss of permissive epigenetic marks. We identify that inflammatory stimuli and accelerated aging as a result of genotoxic stress activate distinct gene networks. These gene networks and associated biological processes are partially overlapping, which is likely driven by specific transcription factor networks, resulting in altered epigenetic signatures and distinct functional (desensitized vs. primed) microglia phenotypes.

Conclusion

This study provides insight into epigenetic profiles and transcription factor networks associated with transcriptional signatures of tolerized and trained microglia in vivo, leading to a better understanding of innate immune memory of microglia.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12974-022-02463-5.

Argonaute proteins regulate a specific network of genes through KLF4 in mouse embryonic stem cells

The Argonaute proteins (AGOs) are well known for their role in post-transcriptional gene silencing in the microRNA (miRNA) pathway. Here we show that in mouse embryonic stem cells, AGO1&2 serve additional functions that go beyond the miRNA pathway. Through the combined deletion of both Agos, we identified a specific set of genes that are uniquely regulated by AGOs but not by the other miRNA biogenesis factors. Deletion of Ago2&1 caused a global reduction of the repressive histone mark H3K27me3 due to downregulation at protein levels of Polycomb repressive complex 2 components. By integrating chromatin accessibility, prediction of transcription factor binding sites, and chromatin immunoprecipitation sequencing data, we identified the pluripotency factor KLF4 as a key modulator of AGO1&2-regulated genes. Our findings revealed a novel axis of gene regulation that is mediated by noncanonical functions of AGO proteins that affect chromatin states and gene expression using mechanisms outside the miRNA pathway.

•

AGO1&2 regulate a specific set of genes in mESCs, independently of the miRNA pathway

•

PRC2 proteins are downregulated in Ago2&1_KO mESCs, leading to H3K27me3 global loss

•

AGO1&2 regulate gene expression through the pluripotency factor KLF4

CDK7/12/13 inhibition targets an oscillating leukemia stem cell network and synergizes with venetoclax in acute myeloid leukemia

The heterogeneous response of acute myeloid leukemia (AML) to current anti‐leukemic therapies is only partially explained by mutational heterogeneity. We previously identified GPR56 as a surface marker associated with poor outcome across genetic groups, which characterizes two leukemia stem cell (LSC)‐enriched compartments with different self‐renewal capacities. How these compartments self‐renew remained unclear. Here, we show that GPR56⁺ LSC compartments are promoted in a complex network involving epithelial‐to‐mesenchymal transition (EMT) regulators besides Rho, Wnt, and Hedgehog (Hh) signaling. Unexpectedly, Wnt pathway inhibition increased the more immature, slowly cycling GPR56⁺CD34⁺ fraction and Hh/EMT gene expression, while Wnt activation caused opposite effects. Our data suggest that the crucial role of GPR56 lies in its ability to co‐activate these opposing signals, thus ensuring the constant supply of both LSC subsets. We show that CDK7 inhibitors suppress both LSC‐enriched subsets in vivo and synergize with the Bcl‐2 inhibitor venetoclax. Our data establish reciprocal transition between LSC compartments as a novel concept underlying the poor outcome in GPR56^high AML and propose combined CDK7 and Bcl‐2 inhibition as LSC‐directed therapy in this disease.

Bibliometric review of ATAC-Seq and its application in gene expression

With recent advances in high-throughput next-generation sequencing, it is possible to describe the regulation and expression of genes at multiple levels. An assay for transposase-accessible chromatin using sequencing (ATAC-seq), which uses Tn5 transposase to sequence protein-free binding regions of the genome, can be combined with chromatin immunoprecipitation coupled with deep sequencing (ChIP-seq) and ribonucleic acid sequencing (RNA-seq) to provide a detailed description of gene expression. Here, we reviewed the literature on ATAC-seq and described the characteristics of ATAC-seq publications. We then briefly introduced the principles of RNA-seq, ChIP-seq and ATAC-seq, focusing on the main features of the techniques. We built a phylogenetic tree from species that had been previously studied by using ATAC-seq. Studies of Mus musculus and Homo sapiens account for approximately 90% of the total ATAC-seq data, while other species are still in the process of accumulating data. We summarized the findings from human diseases and other species, illustrating the cutting-edge discoveries and the role of multi-omics data analysis in current research. Moreover, we collected and compared ATAC-seq analysis pipelines, which allowed biological researchers who lack programming skills to better analyze and explore ATAC-seq data. Through this review, it is clear that multi-omics analysis and single-cell sequencing technology will become the mainstream approach in future research.

Enhancer-silencer transitions in the human genome

Dual-function regulatory elements (REs), acting as enhancers in some cellular contexts and as silencers in others, have been reported to facilitate the precise gene regulatory response to developmental signals in Drosophila melanogaster. However, with few isolated examples detected, dual-function REs in mammals have yet to be systematically studied. We herein investigated this class of REs in the human genome and profiled their activity across multiple cell types. Focusing on enhancer–silencer transitions specific to the development of T cells, we built an accurate deep learning classifier of REs and identified about 12,000 silencers active in primary peripheral blood T cells that act as enhancers in embryonic stem cells. Compared with regular silencers, these dual-function REs are evolving under stronger purifying selection and are enriched for mutations associated with disease phenotypes and altered gene expression. In addition, they are enriched in the loci of transcriptional regulators, such as transcription factors (TFs) and chromatin remodeling genes. Dual-function REs consist of two intertwined but largely distinct sets of binding sites bound by either activating or repressing TFs, depending on the type of RE function in a given cell line. This indicates the recruitment of different TFs for different regulatory modes and a complex DNA sequence composition of these REs with dual activating and repressive encoding. With an estimated >6% of cell type–specific human silencers acting as dual-function REs, this overlooked class of REs requires a specific investigation on how their inherent functional plasticity might be a contributing factor to human diseases.

Unveiling E2F4, TEAD1 and AP-1 as regulatory transcription factors of the replicative senescence program by multi-omics analysis

Senescence, a stable state of growth arrest, affects many physiological and pathophysiological processes, especially aging. Previous work has indicated that transcription factors (TFs) play a role in regulating senescence. However, a systematic study of regulatory TFs during replicative senescence (RS) using multi-omics analysis is still lacking. Here, we generated time-resolved RNA-seq, reduced representation bisulfite sequencing (RRBS) and ATAC-seq datasets during RS of mouse skin fibroblasts, which demonstrated that an enhanced inflammatory response and reduced proliferative capacity were the main characteristics of RS in both the transcriptome and epigenome. Through integrative analysis and genetic manipulations, we found that transcription factors E2F4, TEAD1 and AP-1 are key regulators of RS. Overexpression of E2f4 improved cellular proliferative capacity, attenuated SA-β-Gal activity and changed RS-associated differentially methylated sites (DMSs). Moreover, knockdown of Tead1 attenuated SA-β-Gal activity and partially altered the RS-associated transcriptome. In addition, knockdown of Atf3, one member of AP-1 superfamily TFs, reduced Cdkn2a (p16) expression in pre-senescent fibroblasts. Taken together, the results of this study identified transcription factors regulating the senescence program through multi-omics analysis, providing potential therapeutic targets for anti-aging.

C/EBPB-dependent adaptation to palmitic acid promotes tumor formation in hormone receptor negative breast cancer

Epidemiological studies have established a positive association between obesity and the incidence of postmenopausal breast cancer. Moreover, it is known that obesity promotes stem cell-like properties of breast cancer cells. However, the cancer cell-autonomous mechanisms underlying this correlation are not well defined. Here we demonstrate that obesity-associated tumor formation is driven by cellular adaptation rather than expansion of pre-existing clones within the cancer cell population. While there is no correlation with specific mutations, cellular adaptation to obesity is governed by palmitic acid (PA) and leads to enhanced tumor formation capacity of breast cancer cells. This process is governed epigenetically through increased chromatin occupancy of the transcription factor CCAAT/enhancer-binding protein beta (C/EBPB). Obesity-induced epigenetic activation of C/EBPB regulates cancer stem-like properties by modulating the expression of key downstream regulators including CLDN1 and LCN2. Collectively, our findings demonstrate that obesity drives cellular adaptation to PA drives tumor initiation in the obese setting through activation of a C/EBPB dependent transcriptional network.

Obesity is linked to cancer risk in post-menopausal breast cancer. At the molecular level this is governed by cellular adaption to palmitic acid through epigenetic activation of a C/EBPB-dependent transcriptional network that drives tumor formation.

Deep Mining from Omics Data

Since the advent of high-throughput omics technologies, various molecular data such as genes, transcripts, proteins, and metabolites have been made widely available to researchers. This has afforded clinicians, bioinformaticians, statisticians, and data scientists the opportunity to apply their innovations in feature mining and predictive modeling to a rich data resource to develop a wide range of generalizable prediction models. What has become apparent over the last 10 years is that researchers have adopted deep neural networks (or "deep nets") as their preferred paradigm of choice for complex data modeling due to the superiority of performance over more traditional statistical machine learning approaches, such as support vector machines. A key stumbling block, however, is that deep nets inherently lack transparency and are considered to be a "black box" approach. This naturally makes it very difficult for clinicians and other stakeholders to trust their deep learning models even though the model predictions appear to be highly accurate. In this chapter, we therefore provide a detailed summary of the deep net architectures typically used in omics research, together with a comprehensive summary of the notable "deep feature mining" techniques researchers have applied to open up this black box and provide some insights into the salient input features and why these models behave as they do. We group these techniques into the following three categories: (a) hidden layer visualization and interpretation; (b) input feature importance and impact evaluation; and (c) output layer gradient analysis. While we find that omics researchers have made some considerable gains in opening up the black box through interpretation of the hidden layer weights and node activations to identify salient input features, we highlight other approaches for omics researchers, such as employing deconvolutional network-based approaches and development of bespoke attribute impact measures to enable researchers to better understand the relationships between the input data and hidden layer representations formed and thus the output behavior of their deep nets.

The Clock:Cycle complex is a major transcriptional regulator of Drosophila photoreceptors that protects the eye from retinal degeneration and oxidative stress

The aging eye experiences physiological changes that include decreased visual function and increased risk of retinal degeneration. Although there are transcriptomic signatures in the aging retina that correlate with these physiological changes, the gene regulatory mechanisms that contribute to cellular homeostasis during aging remain to be determined. Here, we integrated ATAC-seq and RNA-seq data to identify 57 transcription factors that showed differential activity in aging Drosophila photoreceptors. These 57 age-regulated transcription factors include two circadian regulators, Clock and Cycle, that showed sustained increased activity during aging. When we disrupted the Clock:Cycle complex by expressing a dominant negative version of Clock (ClkDN) in adult photoreceptors, we observed changes in expression of 15-20% of genes including key components of the phototransduction machinery and many eye-specific transcription factors. Using ATAC-seq, we showed that expression of ClkDN in photoreceptors leads to changes in activity of 37 transcription factors and causes a progressive decrease in global levels of chromatin accessibility in photoreceptors. Supporting a key role for Clock-dependent transcription in the eye, expression of ClkDN in photoreceptors also induced light-dependent retinal degeneration and increased oxidative stress, independent of light exposure. Together, our data suggests that the circadian regulators Clock and Cycle act as neuroprotective factors in the aging eye by directing gene regulatory networks that maintain expression of the phototransduction machinery and counteract oxidative stress.

ID3 promotes homologous recombination via non-transcriptional and transcriptional mechanisms and its loss confers sensitivity to PARP inhibition

The inhibitor of DNA-binding 3 (ID3) is a transcriptional regulator that limits interaction of basic helix-loop-helix transcription factors with their target DNA sequences. We previously reported that ID3 loss is associated with mutational signatures linked to DNA repair defects. Here we demonstrate that ID3 exhibits a dual role to promote DNA double-strand break (DSB) repair, particularly homologous recombination (HR). ID3 interacts with the MRN complex and RECQL helicase to activate DSB repair and it facilitates RAD51 loading and downstream steps of HR. In addition, ID3 promotes the expression of HR genes in response to ionizing radiation by regulating both chromatin accessibility and activity of the transcription factor E2F1. Consistently, analyses of TCGA cancer patient data demonstrate that low ID3 expression is associated with impaired HR. The loss of ID3 leads to sensitivity of tumor cells to PARP inhibition, offering new therapeutic opportunities in ID3-deficient tumors.

Epigenetic reprogramming of airway macrophages promotes polarization and inflammation in muco-obstructive lung disease

Lung diseases, such as cystic fibrosis and COPD, are characterized by mucus obstruction and chronic airway inflammation, but their mechanistic link remains poorly understood. Here, we focus on the function of the mucostatic airway microenvironment on epigenetic reprogramming of airway macrophages (AM) and resulting transcriptomic and phenotypical changes. Using a mouse model of muco-obstructive lung disease (Scnn1b-transgenic), we identify epigenetically controlled, differentially regulated pathways and transcription factors involved in inflammatory responses and macrophage polarization. Functionally, AMs from Scnn1b-transgenic mice have reduced efferocytosis and phagocytosis, and excessive inflammatory responses upon lipopolysaccharide challenge, mediated through enhanced Irf1 function and expression. Ex vivo stimulation of wild-type AMs with native mucus impairs efferocytosis and phagocytosis capacities. In addition, mucus induces gene expression changes, comparable with those observed in AMs from Scnn1b-transgenic mice. Our data show that mucostasis induces epigenetic reprogramming of AMs, leading to changes favoring tissue damage and disease progression. Targeting these altered AMs may support therapeutic approaches in patients with muco-obstructive lung diseases.

Differential analysis of chromatin accessibility and gene expression profiles identifies cis-regulatory elements in rat adipose and muscle

Chromatin accessibility is a key factor influencing gene expression. We optimized the Omni-ATAC-seq protocol and used it together with RNA-seq to investigate cis-regulatory elements in rat white adipose and skeletal muscle, two tissues with contrasting metabolic functions. While promoter accessibility correlated with RNA expression, integration of the two datasets identified tissue-specific differentially accessible regions (DARs) that predominantly localized in intergenic and intron regions. DARs were mapped to differentially expressed (DE) genes enriched in distinct biological processes in each tissue. Randomly selected DE genes were validated by qPCR. Top enriched motifs in DARs predicted binding sites for transcription factors (TFs) showing tissue-specific up-regulation. The correlation between differential chromatin accessibility at a given TF binding motif and differential expression of target genes further supported the functional relevance of that motif. Our study identified cis-regulatory regions that likely play a major role in the regulation of tissue-specific gene expression in adipose and muscle.