Transcription factor enrichment analysis (TFEA) quantifies the activity of multiple transcription factors from a single experiment

Comparative investigation of lung adenocarcinoma and squamous cell carcinoma transcriptome to reveal potential candidate biomarkers: An explainable AI approach

Patients with Non-Small Cell Lung Cancer (NSCLC) present a variety of clinical symptoms, such as dyspnea and chest pain, complicating accurate diagnosis. NSCLC includes subtypes distinguished by histological characteristics, specifically lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). This study aims to compare and identify abnormal gene expression patterns in LUAD and LUSC samples relative to adjacent healthy tissues using an explainable artificial intelligence (XAI) framework. The LASSO algorithm was employed to identify the top gene features in the LUAD and LUSC datasets. An ensemble-based extreme gradient boosting (XGBoost) machine learning (ML) algorithm was trained and interpreted using SHapley Additive exPlanations (SHAP), with top features undergoing biological annotation through survival and functional enrichment analyses. The XAI-based SHAP module addresses the opaque nature of ML models. Notably, 35 and 33 genes were identified for LUAD and LUSC, respectively, using the LASSO algorithm. Performance metrics such as average accuracy and Matthew's correlation coefficient were evaluated. The XGBoost model demonstrated an average accuracy of 99.1 % for LUAD and 98.6 % for LUSC. The SFTPC gene emerged as the most significant feature across both NSCLC subtypes. For LUAD, genes such as STX11, CLEC3B, EMP2, and LYVE1 significantly influenced the XAI-SHAP framework. Conversely, GKN2, OGN, SLC39A8, and MMRN1 were identified for LUSC. Survival analysis and functional validation of these genes highlighted the physiological functions observed to be dysregulated in the NSCLC subtypes. These identified genes have the potential to enhance current medical diagnostics and therapeutics.

TFIIH kinase CDK7 drives cell proliferation through a common core transcription factor network

How cyclin-dependent kinase 7 (CDK7) coordinately regulates the cell cycle and RNA polymerase II transcription remains unclear. Here, high-resolution cryo-electron microscopy revealed how two clinically relevant inhibitors block CDK7 function. In cells, CDK7 inhibition rapidly suppressed transcription, but constitutively active genes were disproportionately affected versus stimulus-responsive. Distinct transcription factors (TFs) regulate constitutive versus stimulus-responsive genes. Accordingly, stimulus-responsive TFs were refractory to CDK7 inhibition whereas constitutively active "core" TFs were repressed. Core TFs (n = 78) are predominantly promoter associated and control cell cycle and proliferative gene expression programs across cell types. Mechanistically, rapid suppression of core TF function can occur through CDK7-dependent phosphorylation changes in core TFs and RB1. Moreover, CDK7 inhibition depleted core TF protein levels within hours, consistent with durable target gene suppression. Thus, a major but unappreciated biological function for CDK7 is regulation of a TF cohort that drives proliferation, revealing an apparent universal mechanism by which CDK7 coordinates RNAPII transcription with cell cycle CDK regulation.

Mediator kinase inhibition suppresses hyperactive interferon signaling in Down syndrome

Hyperactive interferon (IFN) signaling is a hallmark of Down syndrome (DS), a condition caused by Trisomy 21 (T21); strategies that normalize IFN signaling could benefit this population. Mediator-associated kinases CDK8 and CDK19 drive inflammatory responses through incompletely understood mechanisms. Using sibling-matched cell lines with/without T21, we investigated Mediator kinase function in the context of hyperactive IFN in DS over a 75 min to 24 hr timeframe. Activation of IFN-response genes was suppressed in cells treated with the CDK8/CDK19 inhibitor cortistatin A (CA), via rapid suppression of IFN-responsive transcription factor (TF) activity. We also discovered that CDK8/CDK19 affect splicing, a novel means by which Mediator kinases control gene expression. To further probe Mediator kinase function, we completed cytokine screens and metabolomics experiments. Cytokines are master regulators of inflammatory responses; by screening 105 different cytokine proteins, we show that Mediator kinases help drive IFN-dependent cytokine responses at least in part through transcriptional regulation of cytokine genes and receptors. Metabolomics revealed that Mediator kinase inhibition altered core metabolic pathways in cell type-specific ways, and broad upregulation of anti-inflammatory lipid mediators occurred specifically in kinase-inhibited cells during hyperactive IFNγ signaling. A subset of these lipids (e.g. oleamide, desmosterol) serve as ligands for nuclear receptors PPAR and LXR, and activation of these receptors occurred specifically during hyperactive IFN signaling in CA-treated cells, revealing mechanistic links between Mediator kinases, lipid metabolism, and nuclear receptor function. Collectively, our results establish CDK8/CDK19 as context-specific metabolic regulators, and reveal that these kinases control gene expression not only via TFs, but also through metabolic changes and splicing. Moreover, we establish that Mediator kinase inhibition antagonizes IFN signaling through transcriptional, metabolic, and cytokine responses, with implications for DS and other chronic inflammatory conditions.

Histone Deacetylase Inhibitor Largazole Deactivates A Subset of Superenchancers and Causes Mitotic Chromosome Mis-alignment by Suppressing SP1 and BRD4

Histone deacetylase inhibitors have been investigated as potential therapeutic agents for cancer and other diseases. HDIs are known to promote histone acetylation, resulting in an open chromatin conformation and generally increased gene expression. In previous work, we reported that a subset of genes, particularly those regulated by superenhancers, can be suppressed by the HDAC inhibitor largazole. To elucidate the molecular mechanisms underlying gene repression by largazole, we conducted transposase-accessible chromatin sequencing, ChIP-seq, and RNA-seq studies. Our findings revealed that while largazole treatment generally enhances chromatin accessibility, it selectively decreases the accessibility of a subset of superenhancer regions. These genomic regions, showing the most significant changes in the presence of largazole, were enriched with transcription factor binding motifs for SP1, BRD4, CTCF, and YY1. ChIP-seq analysis confirmed reduced binding of BRD4 and SP1 at their respective sites on chromatin, particularly at superenhancers regulating genes such as ID1, c-Myc and MCMs. Largazole exerts its effects by inhibiting DNA replication, RNA processing, and cell cycle progression, partially through the suppression of SP1 expression. Depletion of SP1 by shRNA mimics several key biological effects of largazole and increases cellular sensitivity to the drug. Specific to cell cycle regulation, we demonstrated that largazole disrupts G/M transition by interfering with chromosome alignment during metaphase, a phenotype also observed with SP1 depletion. Our results suggest that largazole exerts its growth-inhibitory effect by suppressing BRD4 and SP1 at super-enhancers, leading to cytostatic responses and mitotic dysfunction.

Roles of Post-Translational Modifications of Transcription Factors Involved in Breast Cancer Hypoxia

BC is the most commonly diagnosed cancer and the second leading cause of cancer death among women worldwide. Cellular stress is a condition that leads to disrupted homeostasis by extrinsic and intrinsic factors. Among other stressors, hypoxia is a driving force for breast cancer (BC) progression and a general hallmark of solid tumors. Thus, intratumoral hypoxia is an important determinant of invasion, metastasis, treatment failure, prognosis, and patient mortality. Acquisition of the epithelial-mesenchymal transition (EMT) phenotype is also a consequence of tumor hypoxia. The cellular response to hypoxia is mainly regulated by the hypoxia signaling pathway, governed by hypoxia-inducible factors (HIFs), mainly HIF1α. HIFs are a family of transcription factors (TFs), which induce the expression of target genes involved in cell survival and proliferation, metabolic reprogramming, angiogenesis, resisting apoptosis, invasion, and metastasis. HIF1α cooperates with a large number of other TFs. In this review, we focused on the crosstalk and cooperation between HIF1α and other TFs involved in the cellular response to hypoxia in BC. We identified a cluster of TFs, proposed as the HIF1α-TF interactome, that orchestrates the transcription of target genes involved in hypoxia, due to their post-translational modifications (PTMs), including phosphorylation/dephosphorylation, ubiquitination/deubiquitination, SUMOylation, hydroxylation, acetylation, S-nitrosylation, and palmitoylation. PTMs of these HIF1α-related TFs drive their stability and activity, degradation and turnover, and the bidirectional translocation between the cytoplasm or plasma membrane and nucleus of BC cells, as well as the transcription/activation of proteins encoded by oncogenes or inactivation of tumor suppressor target genes. Consequently, PTMs of TFs in the HIF1α interactome are crucial regulatory mechanisms that drive the cellular response to oxygen deprivation in BC cells.

RNAseq of peripheral blood mononucleated cells exposed to platelet-rich fibrin and enamel matrix derivatives

Platelet-rich fibrin (PRF) and Enamel Matrix Derivatives (EMD) can support the local regenerative events in periodontal defects. There is reason to suggest that PRF and EMD exert part of their activity by targeting the blood-derived cells accumulating in the early wound healing blastema. However, the impact of PRF and EMD on blood cell response remains to be discovered. To this aim, we have exposed human peripheral blood mononucleated cells (PBMCs) to PRF lysates prepared by a swing-out rotor and EMD, followed by bulk RNA sequencing. A total of 111 and 8 genes are up- and down-regulated by PRF under the premise of an at least log2 two-fold change and a minus log10 significance level of two, respectively. Representative is a characteristic IFN response indicated by various human leukocyte antigens (HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQA2, HLA-DRA, HLA-DRB1, HLA-DRB5), gamma Fc receptors (FCGR1A, FCGR1B, FCGR3B), chemokines (CXCL9-11), and calprotectin (S100A8/9 and S100A12), complement (C1QA/B, C2) and interferon-induced guanylate-binding proteins (GBP1, GBP5). With EMD, 67 and 29 genes are up- and down-regulated, respectively. Characteristic of the upregulated genes are tensins (TNS1 and TNS3). Among the genes downregulated by EMD were epsilon Fc receptors (FCER1A; FCER2), Fc receptor-like proteins (FCRL1, FCRL3) and CX3CR1. Genes commonly upregulated by PRF and EMD were most noticeably NXPH4 and MN1, as well as FN1, MMP14, MERTK, and AXL. Our findings suggest that PRF provokes an inflammatory response, while EMD dampens IgE signaling in peripheral mononucleated blood cells.

Exploring the shared mechanism of fatigue between systemic lupus erythematosus and myalgic encephalomyelitis/chronic fatigue syndrome: monocytic dysregulation and drug repurposing

Background

SLE and ME/CFS both present significant fatigue and share immune dysregulation. The mechanisms underlying fatigue in these disorders remain unclear, and there are no standardized treatments. This study aims to explore shared mechanisms and predict potential therapeutic drugs for fatigue in SLE and ME/CFS.

Methods

Genes associated with SLE and ME/CFS were collected from disease target and clinical sample databases to identify overlapping genes. Bioinformatics analyses, including GO, KEGG, PPI network construction, and key target identification, were performed. ROC curve and correlation analysis of key targets, along with single-cell clustering, were conducted to validate their expression in different cell types. Additionally, an inflammation model was established using THP-1 cells to simulate monocyte activation in both diseases in vitro, and RT-qPCR was used to validate the expression of the key targets. A TF-mRNA-miRNA co-regulatory network was constructed, followed by drug prediction and molecular docking.

Results

Fifty-eight overlapping genes were identified, mainly involved in innate immunity and inflammation. Five key targets were identified (IL1β, CCL2, TLR2, STAT1, IFIH1). Single-cell sequencing revealed that monocytes are enriched with these targets. RT-qPCR confirmed significant upregulation of these targets in the model group. A co-regulatory network was constructed, and ten potential drugs, including suloctidil, N-Acetyl-L-cysteine, simvastatin, ACMC-20mvek, and camptothecin, were predicted. Simvastatin and camptothecin showed high affinity for the key targets.

Conclusion

SLE and ME/CFS share immune and inflammatory pathways. The identified key targets are predominantly enriched in monocytes at the single-cell level, suggesting that classical monocytes may be crucial in linking inflammation and fatigue. RT-qPCR confirmed upregulation in activated monocytes. The TF-mRNA-miRNA network provides a foundation for future research, and drug prediction suggests N-Acetyl-L-cysteine and camptothecin as potential therapies.

Evaluation of large language models for discovery of gene set function

Gene set enrichment is a mainstay of functional genomics, but it relies on gene function databases that are incomplete. Here we evaluate five large language models (LLMs) for their ability to discover the common functions represented by a gene set, supported by molecular rationale and a self-confidence assessment. For curated gene sets from Gene Ontology, GPT-4 suggests functions similar to the curated name in 73% of cases, with higher self-confidence predicting higher similarity. Conversely, random gene sets correctly yield zero confidence in 87% of cases. Other LLMs (GPT-3.5, Gemini Pro, Mixtral Instruct and Llama2 70b) vary in function recovery but are falsely confident for random sets. In gene clusters from omics data, GPT-4 identifies common functions for 45% of cases, fewer than functional enrichment but with higher specificity and gene coverage. Manual review of supporting rationale and citations finds these functions are largely verifiable. These results position LLMs as valuable omics assistants.

Iron-(Fe3+)-Dependent Reactivation of Telomerase Drives Colorectal Cancers

Over-consumption of iron-rich red meat and hereditary or genetic iron overload are associated with an increased risk of colorectal carcinogenesis, yet the mechanistic basis of how metal-mediated signaling leads to oncogenesis remains enigmatic. Using fresh colorectal cancer samples we identify Pirin, an iron sensor, that overcomes a rate-limiting step in oncogenesis, by reactivating the dormant human telomerase reverse transcriptase (hTERT) subunit of the telomerase holoenzyme in an iron-(Fe3+)-dependent manner and thereby drives colorectal cancers. Chemical genetic screens combined with isothermal dose-response fingerprinting and mass spectrometry identified a small molecule SP2509 that specifically inhibits Pirin-mediated hTERT reactivation in colorectal cancers by competing with iron-(Fe3+) binding. Our findings, first to document how metal ions reactivate telomerase, provide a molecular mechanism for the well-known association between red meat and increased incidence of colorectal cancers. Small molecules like SP2509 represent a novel modality to target telomerase that acts as a driver of 90% of human cancers and is yet to be targeted in clinic. Significance: We show how iron-(Fe3+) in collusion with genetic factors reactivates telomerase, providing a molecular mechanism for the association between iron overload and increased incidence of colorectal cancers. Although no enzymatic inhibitors of telomerase have entered the clinic, we identify SP2509, a small molecule that targets telomerase reactivation and function in colorectal cancers.

Cellular zinc status alters chromatin accessibility and binding of p53 to DNA

Zn2+ is an essential metal required by approximately 850 human transcription factors. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Using ATAC-seq to profile regions of accessible chromatin coupled with transcription factor enrichment analysis, we examined how increases and decreases in the labile zinc pool affect chromatin accessibility and transcription factor enrichment. We found 685 transcription factor motifs were differentially enriched, corresponding to 507 unique transcription factors. The pattern of perturbation and the types of transcription factors were notably different at promoters versus intergenic regions, with zinc-finger transcription factors strongly enriched in intergenic regions in elevated Zn2+ To test whether ATAC-seq and transcription factor enrichment analysis predictions correlate with changes in transcription factor binding, we used ChIP-qPCR to profile six p53 binding sites. We found that for five of the six targets, p53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc alter chromatin accessibility and transcription factor binding to DNA.

Transcription factor binding specificities of the oomycete Phytophthora infestans reflect conserved and divergent evolutionary patterns and predict function

BACKGROUND

Identifying the DNA-binding specificities of transcription factors (TF) is central to understanding gene networks that regulate growth and development. Such knowledge is lacking in oomycetes, a microbial eukaryotic lineage within the stramenopile group. Oomycetes include many important plant and animal pathogens such as the potato and tomato blight agent Phytophthora infestans, which is a tractable model for studying life-stage differentiation within the group.

RESULTS

Mining of the P. infestans genome identified 197 genes encoding proteins belonging to 22 TF families. Their chromosomal distribution was consistent with family expansions through unequal crossing-over, which were likely ancient since each family had similar sizes in most oomycetes. Most TFs exhibited dynamic changes in RNA levels through the P. infestans life cycle. The DNA-binding preferences of 123 proteins were assayed using protein-binding oligonucleotide microarrays, which succeeded with 73 proteins from 14 families. Binding sites predicted for representatives of the families were validated by electrophoretic mobility shift or chromatin immunoprecipitation assays. Consistent with the substantial evolutionary distance of oomycetes from traditional model organisms, only a subset of the DNA-binding preferences resembled those of human or plant orthologs. Phylogenetic analyses of the TF families within P. infestans often discriminated clades with canonical and novel DNA targets. Paralogs with similar binding preferences frequently had distinct patterns of expression suggestive of functional divergence. TFs were predicted to either drive life stage-specific expression or serve as general activators based on the representation of their binding sites within total or developmentally-regulated promoters. This projection was confirmed for one TF using synthetic and mutated promoters fused to reporter genes in vivo.

CONCLUSIONS

We established a large dataset of binding specificities for P. infestans TFs, representing the first in the stramenopile group. This resource provides a basis for understanding transcriptional regulation by linking TFs with their targets, which should help delineate the molecular components of processes such as sporulation and host infection. Our work also yielded insight into TF evolution during the eukaryotic radiation, revealing both functional conservation as well as diversification across kingdoms.

Cellular zinc status alters chromatin accessibility and binding of transcription factor p53 to genomic sites

Zinc (Zn2+) is an essential metal required by approximately 2500 proteins. Nearly half of these proteins act on DNA, including > 850 human transcription factors, polymerases, DNA damage response factors, and proteins involved in chromatin architecture. How these proteins acquire their essential Zn2+ cofactor and whether they are sensitive to changes in the labile Zn2+ pool in cells remain open questions. Here, we examine how changes in the labile Zn2+ pool affect chromatin accessibility and transcription factor binding to DNA. We observed both increases and decreases in accessibility in different chromatin regions via ATAC-seq upon treating MCF10A cells with elevated Zn2+ or the Zn2+-specific chelator tris(2-pyridylmethyl)amine (TPA). Transcription factor enrichment analysis was used to correlate changes in chromatin accessibility with transcription factor motifs, revealing 477 transcription factor motifs that were differentially enriched upon Zn2+ perturbation. 186 of these transcription factor motifs were enriched in Zn2+ and depleted in TPA, and the majority correspond to Zn2+ finger transcription factors. We selected TP53 as a candidate to examine how changes in motif enrichment correlate with changes in transcription factor occupancy by ChIP-qPCR. Using publicly available ChIP-seq and nascent transcription datasets, we narrowed the 50,000+ ATAC-seq peaks to 2164 TP53 targets and subsequently selected 6 high-probability TP53 binding sites for testing. ChIP-qPCR revealed that for 5 of the 6 targets, TP53 binding correlates with the local accessibility determined by ATAC-seq. These results demonstrate that changes in labile zinc directly alter chromatin accessibility and transcription factor binding to DNA.

Nascent transcription upon interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines

OBJECTIVES

The interferon-triggered innate immune response has been observed to be under strong diversifying selection to counteract the many pathogens hosts have to defend against. In particular, rewiring of gene transcription regulation allows organisms to rapidly acquire new phenotypes by removing and adding genes into the innate immune gene network. Dissecting the molecular processes by which this rewiring takes place, either by changing the DNA regulatory elements or by changing the activity of the regulators across species, is key to better understand this evolutionary process.

DATA DESCRIPTION

To better comprehend the evolutionary dynamics that have occurred in the initial transcriptional response to interferon in primates, we present Precision Run-On (PRO-seq) datasets made after 1 h of interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines. Further, we tested the difference between using either species' cognate interferon versus using the other orthologous interferon to account for any potential impacts in the interaction of the orthologous interferons with their cellular membrane receptors. This data provides insights into the regulatory mechanisms that drive species-specific responses to environmental perturbations, such as the one driven by the interactions of pathogens and their hosts.

Identification of miRNAs and target genes associated with lymph node metastasis in cervical cancer using bioinformatics analysis

This study was designed to identify the differentially expressed miRNAs (DEMs) and genes (DEGs) in metastatic cervical cancer using bioinformatic tools. In this study, fifty-seven DEMs (48 downregulated and 9 upregulated) were identified, among which miR-4459 and miR-3195 expression was negatively associated with overall survival of cervical cancer patients. Then, 476 target DEGs were determined, and protein-protein interaction (PPI) network was constructed. Seventeen hub genes (LONRF2, CCNE2, AURKA, SYT1, NEGR1, PPP1R12B, GABRP, RAD51, CDK1, FBLN5, PRKG1, CDC6, CACNA1C, MEOX2, ANLN, MYLK, and EDNRB) were finally selected to construct the miRNA-hub gene network. Overall, our study discovered the key miRNAs and mRNAs related to lymph node metastasis (LNM) in cervical cancer, which helps discover candidate therapeutic targets for cervical cancer.

Double DAP-seq uncovered synergistic DNA binding of interacting bZIP transcription factors

Many eukaryotic transcription factors (TF) form homodimer or heterodimer complexes to regulate gene expression. Dimerization of BASIC LEUCINE ZIPPER (bZIP) TFs are critical for their functions, but the molecular mechanism underlying the DNA binding and functional specificity of homo- versus heterodimers remains elusive. To address this gap, we present the double DNA Affinity Purification-sequencing (dDAP-seq) technique that maps heterodimer binding sites on endogenous genomic DNA. Using dDAP-seq we profile twenty pairs of C/S1 bZIP heterodimers and S1 homodimers in Arabidopsis and show that heterodimerization significantly expands the DNA binding preferences of these TFs. Analysis of dDAP-seq binding sites reveals the function of bZIP9 in abscisic acid response and the role of bZIP53 heterodimer-specific binding in seed maturation. The C/S1 heterodimers show distinct preferences for the ACGT elements recognized by plant bZIPs and motifs resembling the yeast GCN4 cis-elements. This study demonstrates the potential of dDAP-seq in deciphering the DNA binding specificities of interacting TFs that are key for combinatorial gene regulation.

Dynamic activity in cis-regulatory elements of leukocytes identifies transcription factor activation and stratifies COVID-19 severity in ICU patients

Transcription factor programs mediating the immune response to coronavirus disease 2019 (COVID-19) are not fully understood. Capturing active transcription initiation from cis-regulatory elements such as enhancers and promoters by capped small RNA sequencing (csRNA-seq), in contrast to capturing steady-state transcripts by conventional RNA-seq, allows unbiased identification of the underlying transcription factor activity and regulatory pathways. Here, we profile transcription initiation in critically ill COVID-19 patients, identifying transcription factor motifs that correlate with clinical lung injury and disease severity. Unbiased clustering reveals distinct subsets of cis-regulatory elements that delineate the cell type, pathway-specific, and combinatorial transcription factor activity. We find evidence of critical roles of regulatory networks, showing that STAT/BCL6 and E2F/MYB regulatory programs from myeloid cell populations are activated in patients with poor disease outcomes and associated with COVID-19 susceptibility genetic variants. More broadly, we demonstrate how capturing acute, disease-mediated changes in transcription initiation can provide insight into the underlying molecular mechanisms and stratify patient disease severity.

MEPP: more transparent motif enrichment by profiling positional correlations

Score-based motif enrichment analysis (MEA) is typically applied to regulatory DNA to infer transcription factors (TFs) that may modulate transcription and chromatin state in different conditions. Most MEA methods determine motif enrichment independent of motif position within a sequence, even when those sequences harbor anchor points that motifs and their bound TFs may functionally interact with in a distance-dependent fashion, such as other TF binding motifs, transcription start sites (TSS), sequencing assay cleavage sites, or other biologically meaningful features. We developed motif enrichment positional profiling (MEPP), a novel MEA method that outputs a positional enrichment profile of a given TF's binding motif relative to key anchor points (e.g. transcription start sites, or other motifs) within the analyzed sequences while accounting for lower-order nucleotide bias. Using transcription initiation and TF binding as test cases, we demonstrate MEPP's utility in determining the sequence positions where motif presence correlates with measures of biological activity, inferring positional dependencies of binding site function. We demonstrate how MEPP can be applied to interpretation and hypothesis generation from experiments that quantify transcription initiation, chromatin structure, or TF binding measurements. MEPP is available for download from https://github.com/npdeloss/mepp.

Single molecule microscopy to profile the effect of zinc status on transcription factor dynamics

The regulation of transcription is a complex process that involves binding of transcription factors (TFs) to specific sequences, recruitment of cofactors and chromatin remodelers, assembly of the pre-initiation complex and recruitment of RNA polymerase II. Increasing evidence suggests that TFs are highly dynamic and interact only transiently with DNA. Single molecule microscopy techniques are powerful approaches for tracking individual TF molecules as they diffuse in the nucleus and interact with DNA. Here we employ multifocus microscopy and highly inclined laminated optical sheet microscopy to track TF dynamics in response to perturbations in labile zinc inside cells. We sought to define whether zinc-dependent TFs sense changes in the labile zinc pool by determining whether their dynamics and DNA binding can be modulated by zinc. We used fluorescently tagged versions of the glucocorticoid receptor (GR), with two C4 zinc finger domains, and CCCTC-binding factor (CTCF), with eleven C2H2 zinc finger domains. We found that GR was largely insensitive to perturbations of zinc, whereas CTCF was significantly affected by zinc depletion and its dwell time was affected by zinc elevation. These results indicate that at least some transcription factors are sensitive to zinc dynamics, revealing a potential new layer of transcriptional regulation.

Integrated genomics approaches identify transcriptional mediators and epigenetic responses to Afghan desert particulate matter in small airway epithelial cells

Military Deployment to Southwest Asia and Afghanistan and exposure to toxic airborne particulates have been associated with an increased risk of developing respiratory disease, collectively termed deployment-related respiratory diseases (DRRDs). Our knowledge about how particulates mediate respiratory disease is limited, precluding the appropriate recognition or management. Central to this limitation is the lack of understanding of how exposures translate into dysregulated cell identity with dysregulated transcriptional programs. The small airway epithelium is involved in both the pathobiology of DRRD and fine particulate matter deposition. To characterize small airway epithelial cell epigenetic and transcriptional responses to Afghan desert particulate matter (APM) and investigate the functional interactions of transcription factors that mediate these responses, we applied two genomics assays, the assay for transposase accessible chromatin with sequencing (ATAC-seq) and Precision Run-on sequencing (PRO-seq). We identified activity changes in a series of transcriptional pathways as candidate regulators of susceptibility to subsequent insults, including signal-dependent pathways, such as loss of cytochrome P450 or P53/P63, and lineage-determining transcription factors, such as GRHL2 loss or TEAD3 activation. We further demonstrated that TEAD3 activation was unique to APM exposure despite similar inflammatory responses when compared with wood smoke particle exposure and that P53/P63 program loss was uniquely positioned at the intersection of signal-dependent and lineage-determining transcriptional programs. Our results establish the utility of an integrated genomics approach in characterizing responses to exposures and identifying genomic targets for the advanced investigation of the pathogenesis of DRRD.

Extracellular matrix stiffness controls cardiac valve myofibroblast activation through epigenetic remodeling

Aortic valve stenosis (AVS) is a progressive fibrotic disease that is caused by thickening and stiffening of valve leaflets. At the cellular level, quiescent valve interstitial cells (qVICs) activate to myofibroblasts (aVICs) that persist within the valve tissue. Given the persistence of myofibroblasts in AVS, epigenetic mechanisms have been implicated. Here, we studied changes that occur in VICs during myofibroblast activation by using a hydrogel matrix to recapitulate different stiffnesses in the valve leaflet during fibrosis. We first compared the chromatin landscape of qVICs cultured on soft hydrogels and aVICs cultured on stiff hydrogels, representing the native and diseased phenotypes respectively. Using assay for transposase-accessible chromatin sequencing (ATAC-Seq), we found that open chromatin regions in aVICs were enriched for transcription factor binding motifs associated with mechanosensing pathways compared to qVICs. Next, we used RNA-Seq to show that the open chromatin regions in aVICs correlated with pro-fibrotic gene expression, as aVICs expressed higher levels of contractile fiber genes, including myofibroblast markers such as alpha smooth muscle actin (αSMA), compared to qVICs. In contrast, chromatin remodeling genes were downregulated in aVICs compared to qVICs, indicating qVICs may be protected from myofibroblast activation through epigenetic mechanisms. Small molecule inhibition of one of these remodelers, CREB Binding Protein (CREBBP), prevented qVICs from activating to aVICs. Notably, CREBBP is more abundant in valves from healthy patients compared to fibrotic valves. Our findings reveal the role of mechanical regulation in chromatin remodeling during VIC activation and quiescence and highlight one potential therapeutic target for treating AVS.

Macrophage programming is regulated by a cooperative interaction between fatty acid binding protein 5 and peroxisome proliferator-activated receptor γ

Resolution of inflammation is an active process that is tightly regulated to achieve repair and tissue homeostasis. In the absence of resolution, persistent inflammation underlies the pathogenesis of chronic lung disease such as chronic obstructive pulmonary disease (COPD) with recurrent exacerbations. Over the course of inflammation, macrophage programming transitions from pro-inflammatory to pro-resolving, which is in part regulated by the nuclear receptor Peroxisome Proliferator-Activated Receptor γ (PPARγ). Our previous work demonstrated an association between Fatty Acid Binding Protein 5 (FABP5) expression and PPARγ activity in peripheral blood mononuclear cells of healthy and COPD patients. However, a role for FABP5 in macrophage programming has not been examined. Here, using a combination of in vitro and in vivo approaches, we demonstrate that FABP5 is necessary for PPARγ activation. In turn, PPARγ acts directly to increase FABP5 expression in primary human alveolar macrophages. We further illustrate that lack of FABP5 expression promotes a pro-inflammatory macrophage programming with increased secretion of pro-inflammatory cytokines and increased chromatin accessibility for pro-inflammatory transcription factors (e.g., NF-κB and MAPK). And finally, real-time cell metabolic analysis using the Seahorse technology shows an inhibition of oxidative phosphorylation in FABP5-deficient macrophages. Taken together, our data indicate that FABP5 and PPARγ reciprocally regulate each other's expression and function, consistent with a novel positive feedback loop between the two factors that mediates macrophage pro-resolving programming. Our studies highlight the importance of defining targets and regulatory mechanisms that control the resolution of inflammation and may serve to inform novel interventional strategies directed towards COPD.

Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries

Background

A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries.

Results

Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ^′ end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription.

Conclusions

Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08352-8).

Profiling Transcription Initiation in Peripheral Leukocytes Reveals Severity-Associated Cis-Regulatory Elements in Critical COVID-19

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the 'active cistrome,' offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.

Profiling Transcription Initiation in Peripheral Leukocytes Reveals Severity-Associated Cis-Regulatory Elements in Critical COVID-19

The contribution of transcription factors (TFs) and gene regulatory programs in the immune response to COVID-19 and their relationship to disease outcome is not fully understood. Analysis of genome-wide changes in transcription at both promoter-proximal and distal cis-regulatory DNA elements, collectively termed the 'active cistrome,' offers an unbiased assessment of TF activity identifying key pathways regulated in homeostasis or disease. Here, we profiled the active cistrome from peripheral leukocytes of critically ill COVID-19 patients to identify major regulatory programs and their dynamics during SARS-CoV-2 associated acute respiratory distress syndrome (ARDS). We identified TF motifs that track the severity of COVID- 19 lung injury, disease resolution, and outcome. We used unbiased clustering to reveal distinct cistrome subsets delineating the regulation of pathways, cell types, and the combinatorial activity of TFs. We found critical roles for regulatory networks driven by stimulus and lineage determining TFs, showing that STAT and E2F/MYB regulatory programs targeting myeloid cells are activated in patients with poor disease outcomes and associated with single nucleotide genetic variants implicated in COVID-19 susceptibility. Integration with single-cell RNA-seq found that STAT and E2F/MYB activation converged in specific neutrophils subset found in patients with severe disease. Collectively we demonstrate that cistrome analysis facilitates insight into disease mechanisms and provides an unbiased approach to evaluate global changes in transcription factor activity and stratify patient disease severity.

The Δ40p53 isoform inhibits p53-dependent eRNA transcription and enables regulation by signal-specific transcription factors during p53 activation

The naturally occurring Δ40p53 isoform heterotetramerizes with wild-type p53 (WTp53) to regulate development, aging, and stress responses. How Δ40p53 alters WTp53 function remains enigmatic because their co-expression causes tetramer heterogeneity. We circumvented this issue with a well-tested strategy that expressed Δ40p53:WTp53 as a single transcript, ensuring a 2:2 tetramer stoichiometry. Human MCF10A cell lines expressing Δ40p53:WTp53, WTp53, or WTp53:WTp53 (as controls) from the native TP53 locus were examined with transcriptomics (precision nuclear run-on sequencing [PRO-seq] and RNA sequencing [RNA-seq]), metabolomics, and other methods. Δ40p53:WTp53 was transcriptionally active, and, although phenotypically similar to WTp53 under normal conditions, it failed to induce growth arrest upon Nutlin-induced p53 activation. This occurred via Δ40p53:WTp53-dependent inhibition of enhancer RNA (eRNA) transcription and subsequent failure to induce mRNA biogenesis, despite similar genomic occupancy to WTp53. A different stimulus (5-fluorouracil [5FU]) also showed Δ40p53:WTp53-specific changes in mRNA induction; however, other transcription factors (TFs; e.g., E2F2) could then drive the response, yielding similar outcomes vs. WTp53. Our results establish that Δ40p53 tempers WTp53 function to enable compensatory responses by other stimulus-specific TFs. Such modulation of WTp53 activity may be an essential physiological function for Δ40p53. Moreover, Δ40p53:WTp53 functional distinctions uncovered herein suggest an eRNA requirement for mRNA biogenesis and that human p53 evolved as a tetramer to support eRNA transcription.

How does Δ40p53, a naturally occurring isoform of p53 that is linked to accelerated aging, alter WTp53 function? Using an innovative approach, this study reveals that Δ40p53 suppresses enhancer RNA transcription and allows other stimulus-specific transcription factors to modulate the p53 transcriptional response.