A blood atlas of COVID-19 defines hallmarks of disease severity and specificity.

Treatment of severe COVID-19 is currently limited by clinical heterogeneity and incomplete description of specific immune biomarkers. We present here a comprehensive multi-omic blood atlas for patients with varying COVID-19 severity in an integrated comparison with influenza and sepsis patients versus healthy volunteers. We identify immune signatures and correlates of host response. Hallmarks of disease severity involved cells, their inflammatory mediators and networks, including progenitor cells and specific myeloid and lymphocyte subsets, features of the immune repertoire, acute phase response, metabolism, and coagulation. Persisting immune activation involving AP-1/p38MAPK was a specific feature of COVID-19. The plasma proteome enabled sub-phenotyping into patient clusters, predictive of severity and outcome. Systems-based integrative analyses including tensor and matrix decomposition of all modalities revealed feature groupings linked with severity and specificity compared to influenza and sepsis. Our approach and blood atlas will support future drug development, clinical trial design, and personalized medicine approaches for COVID-19.

IKAROS and AIOLOS directly regulate AP-1 transcriptional complexes and are essential for NK cell development.

Ikaros transcription factors are essential for adaptive lymphocyte function, yet their role in innate lymphopoiesis is unknown. Using conditional genetic inactivation, we show that Ikzf1/Ikaros is essential for normal natural killer (NK) cell lymphopoiesis and IKZF1 directly represses Cish, a negative regulator of interleukin-15 receptor resulting in impaired interleukin-15 receptor signaling. Both Bcl2l11 and BIM levels, and intrinsic apoptosis were increased in Ikzf1-null NK cells, which in part accounts for NK lymphopenia as both were restored to normal levels when Ikzf1 and Bcl2l11 were co-deleted. Ikzf1-null NK cells presented extensive transcriptional alterations with reduced AP-1 transcriptional complex expression and increased expression of Ikzf2/Helios and Ikzf3/Aiolos. IKZF1 and IKZF3 directly bound AP-1 family members and deletion of both Ikzf1 and Ikzf3 in NK cells resulted in further reductions in Jun/Fos expression and complete loss of peripheral NK cells. Collectively, we show that Ikaros family members are important regulators of apoptosis, cytokine responsiveness and AP-1 transcriptional activity.

Gut CD4+ T cell phenotypes are a continuum molded by microbes, not by TH archetypes.

CD4+ effector lymphocytes (Teff) are traditionally classified by the cytokines they produce. To determine the states that Teff cells actually adopt in frontline tissues in vivo, we applied single-cell transcriptome and chromatin analyses to colonic Teff cells in germ-free or conventional mice or in mice after challenge with a range of phenotypically biasing microbes. Unexpected subsets were marked by the expression of the interferon (IFN) signature or myeloid-specific transcripts, but transcriptome or chromatin structure could not resolve discrete clusters fitting classic helper T cell (TH) subsets. At baseline or at different times of infection, transcripts encoding cytokines or proteins commonly used as TH markers were distributed in a polarized continuum, which was functionally validated. Clones derived from single progenitors gave rise to both IFN-γ- and interleukin (IL)-17-producing cells. Most of the transcriptional variance was tied to the infecting agent, independent of the cytokines produced, and chromatin variance primarily reflected activities of activator protein (AP)-1 and IFN-regulatory factor (IRF) transcription factor (TF) families, not the canonical subset master regulators T-bet, GATA3 or RORγ.

Systematic dissection of sequence features affecting binding specificity of a pioneer factor reveals binding synergy between FOXA1 and AP-1.

Despite the unique ability of pioneer factors (PFs) to target nucleosomal sites in closed chromatin, they only bind a small fraction of their genomic motifs. The underlying mechanism of this selectivity is not well understood. Here, we design a high-throughput assay called chromatin immunoprecipitation with integrated synthetic oligonucleotides (ChIP-ISO) to systematically dissect sequence features affecting the binding specificity of a classic PF, FOXA1, in human A549 cells. Combining ChIP-ISO with in vitro and neural network analyses, we find that (1) FOXA1 binding is strongly affected by co-binding transcription factors (TFs) AP-1 and CEBPB; (2) FOXA1 and AP-1 show binding cooperativity in vitro; (3) FOXA1's binding is determined more by local sequences than chromatin context, including eu-/heterochromatin; and (4) AP-1 is partially responsible for differential binding of FOXA1 in different cell types. Our study presents a framework for elucidating genetic rules underlying PF binding specificity and reveals a mechanism for context-specific regulation of its binding.

Epigenetic programming underpins B cell dysfunction in human SLE.

Systemic lupus erythematosus (SLE) is characterized by the expansion of extrafollicular pathogenic B cells derived from newly activated naive cells. Although these cells express distinct markers, their epigenetic architecture and how it contributes to SLE remain poorly understood. To address this, we determined the DNA methylomes, chromatin accessibility profiles and transcriptomes from five human B cell subsets, including a newly defined effector B cell subset, from subjects with SLE and healthy controls. Our data define a differentiation hierarchy for the subsets and elucidate the epigenetic and transcriptional differences between effector and memory B cells. Importantly, an SLE molecular signature was already established in resting naive cells and was dominated by enrichment of accessible chromatin in motifs for AP-1 and EGR transcription factors. Together, these factors acted in synergy with T-BET to shape the epigenome of expanded SLE effector B cell subsets. Thus, our data define the molecular foundation of pathogenic B cell dysfunction in SLE.

Chromatin-bound RB targets promoters, enhancers, and CTCF-bound loci and is redistributed by cell-cycle progression.

The interaction of RB with chromatin is key to understanding its molecular functions. Here, for first time, we identify the full spectrum of chromatin-bound RB. Rather than exclusively binding promoters, as is often described, RB targets three fundamentally different types of loci (promoters, enhancers, and insulators), which are largely distinguishable by the mutually exclusive presence of E2F1, c-Jun, and CTCF. While E2F/DP facilitates RB association with promoters, AP-1 recruits RB to enhancers. Although phosphorylation in CDK sites is often portrayed as releasing RB from chromatin, we show that the cell cycle redistributes RB so that it enriches at promoters in G1 and at non-promoter sites in cycling cells. RB-bound promoters include the classic E2F-targets and are similar between lineages, but RB-bound enhancers associate with different categories of genes and vary between cell types. Thus, RB has a well-preserved role controlling E2F in G1, and it targets cell-type-specific enhancers and CTCF sites when cells enter S-phase.

A non-canonical repressor function of JUN restrains YAP activity and liver cancer growth.

Yes-associated protein (YAP) and its homolog, transcriptional coactivator with PDZ-binding motif (TAZ), are the main transcriptional downstream effectors of the Hippo pathway. Decreased Hippo pathway activity leads to nuclear translocation of YAP/TAZ where they interact with TEAD transcription factors to induce target gene expression. Unrestrained YAP/TAZ activity can lead to excessive growth and tumor formation in a short time, underscoring the evolutionary need for tight control of these two transcriptional coactivators. Here, we report that the AP-1 component JUN acts as specific repressor of YAP/TAZ at joint target sites to decrease YAP/TAZ activity. This function of JUN is independent of its heterodimeric AP-1 partner FOS and the canonical AP-1 function. Since expression of JUN is itself induced by YAP/TAZ, our work identifies a JUN-dependent negative feedback loop that buffers YAP/TAZ activity at joint genomic sites. This negative feedback loop gets disrupted in liver cancer to unlock the full oncogenic potential of YAP/TAZ. Our results thus demonstrate an additional layer of control for the interplay of YAP/TAZ and AP-1.

BACH2 regulates CD8(+) T cell differentiation by controlling access of AP-1 factors to enhancers.

T cell antigen receptor (TCR) signaling drives distinct responses depending on the differentiation state and context of CD8(+) T cells. We hypothesized that access of signal-dependent transcription factors (TFs) to enhancers is dynamically regulated to shape transcriptional responses to TCR signaling. We found that the TF BACH2 restrains terminal differentiation to enable generation of long-lived memory cells and protective immunity after viral infection. BACH2 was recruited to enhancers, where it limited expression of TCR-driven genes by attenuating the availability of activator protein-1 (AP-1) sites to Jun family signal-dependent TFs. In naive cells, this prevented TCR-driven induction of genes associated with terminal differentiation. Upon effector differentiation, reduced expression of BACH2 and its phosphorylation enabled unrestrained induction of TCR-driven effector programs.

Plasma cell differentiation is coupled to division-dependent DNA hypomethylation and gene regulation.

The epigenetic processes that regulate antibody-secreting plasma cells are not well understood. Here, analysis of plasma cell differentiation revealed DNA hypomethylation of 10% of CpG loci that were overrepresented at enhancers. Inhibition of DNA methylation enhanced plasma cell commitment in a cell-division-dependent manner. Analysis of B cells differentiating in vivo stratified by cell division revealed a fivefold increase in mRNA transcription coupled to DNA hypomethylation. Demethylation occurred first at binding motifs for the transcription factors NF-κB and AP-1 and later at those for the transcription factors IRF and Oct-2 and was coincident with activation and differentiation gene-expression programs in a cell-division-dependent manner. These data provide mechanistic insight into cell-division-coupled transcriptional and epigenetic reprogramming and suggest that DNA hypomethylation reflects the cis-regulatory history of plasma cell differentiation.

Induction of AP-1 by YAP/TAZ contributes to cell proliferation and organ growth.

Yes-associated protein (YAP) and its homolog transcriptional coactivator with PDZ-binding motif (TAZ) are key effectors of the Hippo pathway to control cell growth and organ size, of which dysregulation yields to tumorigenesis or hypertrophy. Upon activation, YAP/TAZ translocate into the nucleus and bind to TEAD transcription factors to promote transcriptional programs for proliferation or cell specification. Immediate early genes, represented by AP-1 complex, are rapidly induced and control later-phase transcriptional program to play key roles in tumorigenesis and organ maintenance. Here, we report that YAP/TAZ directly promote FOS transcription that in turn contributes to the biological function of YAP/TAZ. YAP/TAZ bind to the promoter region of FOS to stimulate its transcription. Deletion of YAP/TAZ blocks the induction of immediate early genes in response to mitogenic stimuli. FOS induction contributes to expression of YAP/TAZ downstream target genes. Genetic deletion or chemical inhibition of AP-1 suppresses growth of YAP-driven cancer cells, such as Lats1/2-deficient cancer cells as well as Gαq/11 mutated uveal melanoma. Furthermore, AP-1 inhibition almost completely abrogates the hepatomegaly induced by YAP overexpression. Our findings reveal a feed-forward interplay between immediate early transcription of AP-1 and Hippo pathway function.

Machine learning identifies activation of RUNX/AP-1 as drivers of mesenchymal and fibrotic regulatory programs in gastric cancer.

Gastric cancer (GC) is the fifth most common cancer worldwide and is a heterogeneous disease. Among GC subtypes, the mesenchymal phenotype (Mes-like) is more invasive than the epithelial phenotype (Epi-like). Although gene expression of the epithelial-to-mesenchymal transition (EMT) has been studied, the regulatory landscape shaping this process is not fully understood. Here we use ATAC-seq and RNA-seq data from a compendium of GC cell lines and primary tumors to detect drivers of regulatory state changes and their transcriptional responses. Using the ATAC-seq data, we developed a machine learning approach to determine the transcription factors (TFs) regulating the subtypes of GC. We identified TFs driving the mesenchymal (RUNX2, ZEB1, SNAI2, AP-1 dimer) and the epithelial (GATA4, GATA6, KLF5, HNF4A, FOXA2, GRHL2) states in GC. We identified DNA copy number alterations associated with dysregulation of these TFs, specifically deletion of GATA4 and amplification of MAPK9 Comparisons with bulk and single-cell RNA-seq data sets identified activation toward fibroblast-like epigenomic and expression signatures in Mes-like GC. The activation of this mesenchymal fibrotic program is associated with differentially accessible DNA cis-regulatory elements flanking upregulated mesenchymal genes. These findings establish a map of TF activity in GC and highlight the role of copy number driven alterations in shaping epigenomic regulatory programs as potential drivers of GC heterogeneity and progression.

Sox2 functions as a sequence-specific DNA sensor in neutrophils to initiate innate immunity against microbial infection.

Neutrophils express Toll-like receptors (TLRs) for the recognition of conserved bacterial elements to initiate antimicrobial responses. However, whether other cytosolic DNA sensors are expressed by neutrophils remains elusive. Here we found constitutive expression of the transcription factor Sox2 in the cytoplasm of mouse and human neutrophils. Neutrophil-specific Sox2 deficiency exacerbated bacterial infection. Sox2 directly recognized microbial DNA through its high-mobility-group (HMG) domain. Upon challenge with bacterial DNA, Sox2 dimerization was needed to activate a complex of the kinase TAK1 and its binding partner TAB2, which led to activation of the transcription factors NF-κB and AP-1 in neutrophils. Deficiency in TAK1 or TAB2 impaired Sox2-mediated antibacterial immunity. Overall, we reveal a previously unrecognized role for Sox2 as a cytosolic sequence-specific DNA sensor in neutrophils, which might provide potential therapeutic strategies for the treatment of infectious diseases.

Liver cancer development driven by the AP-1/c-Jun~Fra-2 dimer through c-Myc.

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death. HCC incidence is on the rise, while treatment options remain limited. Thus, a better understanding of the molecular pathways involved in HCC development has become a priority to guide future therapies. While previous studies implicated the Activator Protein-1 (AP-1) (Fos/Jun) transcription factor family members c-Fos and c-Jun in HCC formation, the contribution of Fos-related antigens (Fra-) 1 and 2 is unknown. Here, we show that hepatocyte-restricted expression of a single chain c-Jun~Fra-2 protein, which functionally mimics the c-Jun/Fra-2 AP-1 dimer, results in spontaneous HCC formation in c-Jun~Fra-2hep mice. Several hallmarks of human HCC, such as cell cycle dysregulation and the expression of HCC markers are observed in liver tumors arising in c-Jun~Fra-2hep mice. Tumorigenesis occurs in the context of mild inflammation, low-grade fibrosis, and Pparγ-driven dyslipidemia. Subsequent analyses revealed increased expression of c-Myc, evidently under direct regulation by AP-1 through a conserved distal 3' enhancer. Importantly, c-Jun~Fra-2-induced tumors revert upon switching off transgene expression, suggesting oncogene addiction to the c-Jun~Fra-2 transgene. Tumors escaping reversion maintained c-Myc and c-Myc target gene expression, likely due to increased c-Fos. Interfering with c-Myc in established tumors using the Bromodomain and Extra-Terminal motif inhibitor JQ-1 diminished liver tumor growth in c-Jun~Fra-2 mutant mice. Thus, our data establish c-Jun~Fra-2hep mice as a model to study liver tumorigenesis and identify the c-Jun/Fra-2-Myc interaction as a potential target to improve HCC patient stratification and/or therapy.

Anti-Inflammatory Chromatinscape Suggests Alternative Mechanisms of Glucocorticoid Receptor Action.

Despite the widespread use of glucocorticoids (GCs), their anti-inflammatory effects are not understood mechanistically. Numerous investigations have examined the effects of glucocorticoid receptor (GR) activation prior to inflammatory challenges. However, clinical situations are emulated by a GC intervention initiated in the midst of rampant inflammatory responses. To characterize the effects of a late GC treatment, we profiled macrophage transcriptional and chromatinscapes with Dexamethasone (Dex) treatment before or after stimulation by lipopolysaccharide (LPS). The late activation of GR had a similar gene-expression profile as from GR pre-activation, while ameliorating the disruption of metabolic genes. Chromatin occupancy of GR was not predictive of Dex-regulated gene expression, contradicting the "trans-repression by tethering" model. Rather, GR activation resulted in genome-wide blockade of NF-κB interaction with chromatin and directly induced inhibitors of NF-κB and AP-1. Our investigation using GC treatments with clinically relevant timing highlights mechanisms underlying GR actions for modulating the "inflamed epigenome."

The enhancer RNA, AANCR, regulates APOE expression in astrocytes and microglia.

Enhancers, critical regulatory elements within the human genome, are often transcribed into enhancer RNAs. The dysregulation of enhancers leads to diseases collectively termed enhanceropathies. While it is known that enhancers play a role in diseases by regulating gene expression, the specific mechanisms by which individual enhancers cause diseases are not well understood. Studies of individual enhancers are needed to fill this gap. This study delves into the role of APOE-activating noncoding RNA, AANCR, in the central nervous system, elucidating its function as a genetic modifier in Alzheimer's Disease. We employed RNA interference, RNaseH-mediated degradation, and single-molecule RNA fluorescence in situ hybridization to demonstrate that mere transcription of AANCR is insufficient; rather, its transcripts are crucial for promoting APOE expression. Our findings revealed that AANCR is induced by ATM-mediated ERK phosphorylation and subsequent AP-1 transcription factor activation. Once activated, AANCR enhances APOE expression, which in turn imparts an inflammatory phenotype to astrocytes. These findings demonstrate that AANCR is a key enhancer RNA in some cell types within the nervous system, pivotal for regulating APOE expression and influencing inflammatory responses, underscoring its potential as a therapeutic target in neurodegenerative diseases.

Functional analysis of protein interactions using coupled bi-fluorescence complementation/GFP nanobody techniques.

Transcription factors (TFs) form homo- or hetero-dimeric DNA binding complexes along with associated co-regulators that can have transcriptional repressor or activator functions. Defining the specific composition of the complexes is therefore key to understanding their biological role. Here, we utilized bimolecular fluorescence complementation (BiFC) to visualize the formation of defined TF dimers and associated co-regulators derived from the activator protein-1 (AP-1) and myocyte enhancer factor 2 (MEF2) families. Firstly, BiFC signals were observed in cells co-expressing TFs tagged with complimentary combinations of the split fluorescent protein, demonstrating the engineered formation of defined dimer complexes. Next, we applied this approach and determined that defined AP-1 dimers localized at discrete sub-nuclear locations. Subsequently, a combination of BiFC coupled with GFP binding peptide (GBP)-nanotrap allowed observation of protein-protein interactions between a co-regulator, HDAC4, and defined BiFC-MEF2 engineered dimers. To determine transactivation properties of defined TF dimers in a cellular system, the Gal4-DNA binding domain fused to GBP was utilized to assess the transcriptional properties of the BiFC-TF dimers using a generically applicable Gal4/UAS luciferase reporter gene assay system. Here, we report efficacy of a BiFC/GBP-nanobody approach that allows engineering, visualization, and functional analysis of defined TF dimers.

IRF2BP2 counteracts the ATF7/JDP2 AP-1 heterodimer to prevent inflammatory overactivation in acute myeloid leukemia (AML) cells.

Acute myeloid leukemia (AML) is a hematological malignancy characterized by abnormal proliferation and accumulation of immature myeloid cells in the bone marrow. Inflammation plays a crucial role in AML progression, but excessive activation of cell-intrinsic inflammatory pathways can also trigger cell death. IRF2BP2 is a chromatin regulator implicated in AML pathogenesis, although its precise role in this disease is not fully understood. In this study, we demonstrate that IRF2BP2 interacts with the AP-1 heterodimer ATF7/JDP2, which is involved in activating inflammatory pathways in AML cells. We show that IRF2BP2 is recruited by the ATF7/JDP2 dimer to chromatin and counteracts its gene-activating function. Loss of IRF2BP2 leads to overactivation of inflammatory pathways, resulting in strongly reduced proliferation. Our research indicates that a precise equilibrium between activating and repressive transcriptional mechanisms creates a pro-oncogenic inflammatory environment in AML cells. The ATF7/JDP2-IRF2BP2 regulatory axis is likely a key regulator of this process and may, therefore, represent a promising therapeutic vulnerability for AML. Thus, our study provides new insights into the molecular mechanisms underlying AML pathogenesis and identifies a potential therapeutic target for AML treatment.

Adenovirus small E1A directs activation of Alu transcription at YAP/TEAD- and AP-1-bound enhancers through interactions with the EP400 chromatin remodeler.

Alu retrotransposons, which form the largest family of mobile DNA elements in the human genome, have recently come to attention as a potential source of regulatory novelties, most notably by participating in enhancer function. Even though Alu transcription by RNA polymerase III is subjected to tight epigenetic silencing, their expression has long been known to increase in response to various types of stress, including viral infection. Here we show that, in primary human fibroblasts, adenovirus small e1a triggered derepression of hundreds of individual Alus by promoting TFIIIB recruitment by Alu-bound TFIIIC. Epigenome profiling revealed an e1a-induced decrease of H3K27 acetylation and increase of H3K4 monomethylation at derepressed Alus, making them resemble poised enhancers. The enhancer nature of e1a-targeted Alus was confirmed by the enrichment, in their upstream regions, of the EP300/CBP acetyltransferase, EP400 chromatin remodeler and YAP1 and FOS transcription factors. The physical interaction of e1a with EP400 was critical for Alu derepression, which was abrogated upon EP400 ablation. Our data suggest that e1a targets a subset of enhancer Alus whose transcriptional activation, which requires EP400 and is mediated by the e1a-EP400 interaction, may participate in the manipulation of enhancer activity by adenoviruses.

Assigning pathogenicity for TAB2 variants using a novel scalable functional assay and expanding TAB2 disease spectrum.

Haploinsufficiency of TGF-beta-activated kinase 1 (MAP3K7) binding protein 2 (TAB2) has been associated with congenital heart disease and more recently multiorgan structural abnormalities. Missense variant represents a major proportion of non-synonymous TAB2 variants reported in gnomAD (295/576) and Clinvar (16/73), most of which are variants of uncertain significance (VUSs). However, interpretation of TAB2 missense variants remains challenging because of lack of functional assays. To address this issue, we established a cell-based luciferase assay that enables high-throughput screening of TAB2 variants to assess the functional consequence for predicting variant pathogenicity. Using this platform, we screened 47 TAB2 variants including five pathogenic controls and one benign control, and the results showed that the transcriptional activity of activator protein 1 (AP-1) but not nuclear factor kappa B predicts the TAB2 variant pathogenicity. This assay provides accurate functional readout for both loss-of-function (LOF) and gain-of-function variants, which are associated with distinct phenotypes. In all, 22 out of 32 tested VUSs were reclassified. Genotype-Phenotype association showed that most patients with partial LOF variants do not exhibit congenital heart disease but high frequency of developmental delay, hypotonia and dysmorphic features, which suggests that genetic testing for TAB2 is needed for a broader spectrum of patients with more diverse phenotypes. Molecular modeling with Npl4 zinc finger (NZF) domain variants revealed that the stability of the NZF domain in TAB2 protein is crucial for AP-1 activation. In conclusion, we developed a highly effective functional assay for TAB2 variant prediction and interpretation.

AP-1 signaling modulates cardiac fibroblast stress responses.

Matrix remodeling outcomes largely dictate patient survival post myocardial infarction. Moreover, human-restricted noncoding regulatory elements have been shown to worsen fibrosis, but their mechanism of action remains elusive. Here, we demonstrate, using induced pluripotent stem cell-derived cardiac fibroblasts (iCFs), that inflammatory ligands abundant in the remodeling heart after infarction activate AP-1 transcription factor signaling pathways resulting in fibrotic responses. This observed signaling induces deposition of fibronectin matrix and is further capable of supporting immune cell adhesion; pathway inhibition blocks iCF matrix production and cell adhesion. Polymorphisms in the noncoding regulatory elements within the 9p21 locus (also referred to as ANRIL) redirect stress programs, and in iCFs, they transcriptionally silence the AP-1 inducible transcription factor GATA5. The presence of these polymorphisms modulate iCF matrix production and assembly and reduce cell-cell signaling. These data suggest that this signaling axis is a critical modulator of cardiac disease models and might be influenced by noncoding regulatory elements.