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.

Macrophages derived from human induced pluripotent stem cells (iPSCs) serve as a high-fidelity cellular model for investigating HIV-1, dengue, and influenza viruses.

Macrophages are important target cells for diverse viruses and thus represent a valuable system for studying virus biology. Isolation of primary human macrophages is done by culture of dissociated tissues or from differentiated blood monocytes, but these methods are both time consuming and result in low numbers of recovered macrophages. Here, we explore whether macrophages derived from human induced pluripotent stem cells (iPSCs)-which proliferate indefinitely and potentially provide unlimited starting material-could serve as a faithful model system for studying virus biology. Human iPSC-derived monocytes were differentiated into macrophages and then infected with HIV-1, dengue virus, or influenza virus as model human viruses. We show that iPSC-derived macrophages support the replication of these viruses with kinetics and phenotypes similar to human blood monocyte-derived macrophages. These iPSC-derived macrophages were virtually indistinguishable from human blood monocyte-derived macrophages based on surface marker expression (flow cytometry), transcriptomics (RNA sequencing), and chromatin accessibility profiling. iPSC lines were additionally generated from non-human primate (chimpanzee) fibroblasts. When challenged with dengue virus, human and chimpanzee iPSC-derived macrophages show differential susceptibility to infection, thus providing a valuable resource for studying the species-tropism of viruses. We also show that blood- and iPSC-derived macrophages both restrict influenza virus at a late stage of the virus lifecycle. Collectively, our results substantiate iPSC-derived macrophages as an alternative to blood monocyte-derived macrophages for the study of virus biology. IMPORTANCE: Macrophages have complex relationships with viruses: while macrophages aid in the removal of pathogenic viruses from the body, macrophages are also manipulated by some viruses to serve as vessels for viral replication, dissemination, and long-term persistence. Here, we show that iPSC-derived macrophages are an excellent model that can be exploited in virology.

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.

Suppression of p53 response by targeting p53-Mediator binding with a stapled peptide.

DNA-binding transcription factors (TFs) remain challenging to target with molecular probes. Many TFs function in part through interaction with Mediator, a 26-subunit complex that controls RNA polymerase II activity genome-wide. We sought to block p53 function by disrupting the p53-Mediator interaction. Through rational design and activity-based screening, we characterize a stapled peptide, with functional mimics of both p53 activation domains, that blocks p53-Mediator binding and selectively inhibits p53-dependent transcription in human cells; importantly, this "bivalent" peptide has negligible impact, genome-wide, on non-p53 target genes. Our proof-of-concept strategy circumvents the TF entirely and targets the TF-Mediator interface instead, with desired functional outcomes (i.e., selective inhibition of p53 activation). Furthermore, these results demonstrate that TF activation domains represent viable starting points for Mediator-targeting molecular probes, as an alternative to large compound libraries. Different TFs bind Mediator through different subunits, suggesting this strategy could be broadly applied to selectively alter gene expression programs.

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.

Characterization of Global Gene Expression, Regulation of Metal Ions, and Infection Outcomes in Immune-Competent 129S6 Mouse Macrophages.

Nutritional immunity involves cellular and physiological responses to invading pathogens, such as limiting iron, increasing exposure to bactericidal copper, and altering zinc to restrict the growth of pathogens. Here, we examine infection of bone marrow-derived macrophages from 129S6/SvEvTac mice by Salmonella enterica serovar Typhimurium. The 129S6/SvEvTac mice possess a functional Slc11a1 (Nramp-1), a phagosomal transporter of divalent cations that plays an important role in modulating metal availability to the pathogen. We carried out global RNA sequencing upon treatment with live or heat-killed Salmonella at 2 h and 18 h postinfection and observed widespread changes in metal transport, metal-dependent genes, and metal homeostasis genes, suggesting significant remodeling of iron, copper, and zinc availability by host cells. Changes in host cell gene expression suggest infection increases cytosolic zinc while simultaneously limiting zinc within the phagosome. Using a genetically encoded sensor, we demonstrate that cytosolic labile zinc increases 45-fold at 12 h postinfection. Further, manipulation of zinc in the medium alters bacterial clearance and replication, with zinc depletion inhibiting both processes. Comparing the transcriptomic changes to published data on infection of C57BL/6 macrophages revealed notable differences in metal regulation and the global immune response. Our results reveal that 129S6 macrophages represent a distinct model system compared to C57BL/6 macrophages. Further, our results indicate that manipulation of zinc at the host-pathogen interface is more nuanced than that of iron or copper. The 129S6 macrophages leverage intricate means of manipulating zinc availability and distribution to limit the pathogen's access to zinc, while simultaneously ensuring sufficient zinc to support the immune response.

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

Detecting changes in the activity of a transcription factor (TF) in response to a perturbation provides insights into the underlying cellular process. Transcription Factor Enrichment Analysis (TFEA) is a robust and reliable computational method that detects positional motif enrichment associated with changes in transcription observed in response to a perturbation. TFEA detects positional motif enrichment within a list of ranked regions of interest (ROIs), typically sites of RNA polymerase initiation inferred from regulatory data such as nascent transcription. Therefore, we also introduce muMerge, a statistically principled method of generating a consensus list of ROIs from multiple replicates and conditions. TFEA is broadly applicable to data that informs on transcriptional regulation including nascent transcription (eg. PRO-Seq), CAGE, histone ChIP-Seq, and accessibility data (e.g., ATAC-Seq). TFEA not only identifies the key regulators responding to a perturbation, but also temporally unravels regulatory networks with time series data. Consequently, TFEA serves as a hypothesis-generating tool that provides an easy, rigorous, and cost-effective means to broadly assess TF activity yielding new biological insights.

Just 2% of SARS-CoV-2-positive individuals carry 90% of the virus circulating in communities.

We analyze data from the fall 2020 pandemic response efforts at the University of Colorado Boulder, where more than 72,500 saliva samples were tested for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using qRT-PCR. All samples were collected from individuals who reported no symptoms associated with COVID-19 on the day of collection. From these, 1,405 positive cases were identified. The distribution of viral loads within these asymptomatic individuals was indistinguishable from what has been previously observed in symptomatic individuals. Regardless of symptomatic status, ∼50% of individuals who test positive for SARS-CoV-2 seem to be in noninfectious phases of the disease, based on having low viral loads in a range from which live virus has rarely been isolated. We find that, at any given time, just 2% of individuals carry 90% of the virions circulating within communities, serving as viral "supercarriers" and possibly also superspreaders.

The MUC5B-associated variant rs35705950 resides within an enhancer subject to lineage- and disease-dependent epigenetic remodeling.

The G/T transversion rs35705950, located approximately 3 kb upstream of the MUC5B start site, is the cardinal risk factor for idiopathic pulmonary fibrosis (IPF). Here, we investigate the function and chromatin structure of this -3 kb region and provide evidence that it functions as a classically defined enhancer subject to epigenetic programming. We use nascent transcript analysis to show that RNA polymerase II loads within 10 bp of the G/T transversion site, definitively establishing enhancer function for the region. By integrating Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) analysis of fresh and cultured human airway epithelial cells with nuclease sensitivity data, we demonstrate that this region is in accessible chromatin that affects the expression of MUC5B. Through applying paired single-nucleus RNA- and ATAC-seq to frozen tissue from IPF lungs, we extend these findings directly to disease, with results indicating that epigenetic programming of the -3 kb enhancer in IPF occurs in both MUC5B-expressing and nonexpressing lineages. In aggregate, our results indicate that the MUC5B-associated variant rs35705950 resides within an enhancer that is subject to epigenetic remodeling and contributes to pathologic misexpression in IPF.

Competition between PRC2.1 and 2.2 subcomplexes regulates PRC2 chromatin occupancy in human stem cells.

Polycomb repressive complex 2 (PRC2) silences expression of developmental transcription factors in pluripotent stem cells by methylating lysine 27 on histone H3. Two mutually exclusive subcomplexes, PRC2.1 and PRC2.2, are defined by the set of accessory proteins bound to the core PRC2 subunits. Here we introduce separation-of-function mutations into the SUZ12 subunit of PRC2 to drive it into a PRC2.1 or 2.2 subcomplex in human induced pluripotent stem cells (iPSCs). We find that PRC2.2 occupies polycomb target genes at low levels and that homeobox transcription factors are upregulated when this complex is exclusively present. In contrast with previous studies, we find that chromatin occupancy of PRC2 increases drastically when it is forced to form PRC2.1. Additionally, several cancer-associated mutations also coerce formation of PRC2.1. We suggest that PRC2 chromatin occupancy can be altered in the context of disease or development by tuning the ratio of PRC2.1 to PRC2.2.

Combining signal and sequence to detect RNA polymerase initiation in ATAC-seq data.

The assay for transposase-accessible chromatin followed by sequencing (ATAC-seq) is an inexpensive protocol for measuring open chromatin regions. ATAC-seq is also relatively simple and requires fewer cells than many other high-throughput sequencing protocols. Therefore, it is tractable in numerous settings where other high throughput assays are challenging to impossible. Hence it is important to understand the limits of what can be inferred from ATAC-seq data. In this work, we leverage ATAC-seq to predict the presence of nascent transcription. Nascent transcription assays are the current gold standard for identifying regions of active transcription, including markers for functional transcription factor (TF) binding. We combine mapped short reads from ATAC-seq with the underlying peak sequence, to determine regions of active transcription genome-wide. We show that a hybrid signal/sequence representation classified using recurrent neural networks (RNNs) can identify these regions across different cell types.

Remodeling of Zn2+ homeostasis upon differentiation of mammary epithelial cells.

Zinc is the second most abundant transition metal in humans and an essential nutrient required for growth and development of newborns. During lactation, mammary epithelial cells differentiate into a secretory phenotype, uptake zinc from blood circulation, and export it into mother's milk. At the cellular level, many zinc-dependent cellular processes, such as transcription, metabolism of nutrients, and proliferation are involved in the differentiation of mammary epithelial cells. Using mouse mammary epithelial cells as a model system, we investigated the remodeling of zinc homeostasis during differentiation induced by treatment with the lactogenic hormones cortisol and prolactin. RNA-Seq at different stages of differentiation revealed changes in global gene expression, including genes encoding zinc-dependent proteins and regulators of zinc homeostasis. Increases in mRNA levels of three zinc homeostasis genes, Slc39a14 (ZIP14) and metallothioneins (MTs) I and II were induced by cortisol but not by prolactin. The cortisol-induced increase was partially mediated by the nuclear glucocorticoid receptor signaling pathway. An increase in the cytosolic labile Zn2+ pool was also detected in lactating mammary cells, consistent with upregulation of MTs. We found that the zinc transporter ZIP14 was important for the expression of a major milk protein, whey acid protein (WAP), as knockdown of ZIP14 dramatically decreased WAP mRNA levels. In summary, our study demonstrated remodeling of zinc homeostasis upon differentiation of mammary epithelial cells resulting in changes in cytosolic Zn2+ and differential expression of zinc homeostasis genes, and these changes are important for establishing the lactation phenotype.

Two-stage ML Classifier for Identifying Host Protein Targets of the Dengue Protease.

Flaviviruses such as dengue encode a protease that is essential for viral replication. The protease functions by cleaving well-conserved positions in the viral polyprotein. In addition to the viral polyprotein, the dengue protease cleaves at least one host protein involved in immune response. This raises the question, what other host proteins are targeted and cleaved? Here we present a new computational method for identifying putative host protein targets of the dengue virus protease. Our method relies on biochemical and secondary structure features at the known cleavage sites in the viral polyprotein in a two-stage classification process to identify putative cleavage targets. The accuracy of our predictions scaled inversely with evolutionary distance when we applied it to the known cleavage sites of several other flaviviruses-a good indication of the validity of our predictions. Ultimately, our classifier identified 257 human protein sites possessing both a similar target motif and accessible local structure. These proteins are promising candidates for further investigation. As the number of viral sequences expands, our method could be adopted to predict host targets of other flaviviruses.

Nascent transcript analysis of glucocorticoid crosstalk with TNF defines primary and cooperative inflammatory repression.

The glucocorticoid receptor (NR3C1, also known as GR) binds to specific DNA sequences and directly induces transcription of anti-inflammatory genes that contribute to cytokine repression, frequently in cooperation with NF-kB. Whether inflammatory repression also occurs through local interactions between GR and inflammatory gene regulatory elements has been controversial. Here, using global run-on sequencing (GRO-seq) in human airway epithelial cells, we show that glucocorticoid signaling represses transcription within 10 min. Many repressed regulatory regions reside within "hyper-ChIPable" genomic regions that are subject to dynamic, yet nonspecific, interactions with some antibodies. When this artifact was accounted for, we determined that transcriptional repression does not require local GR occupancy. Instead, widespread transcriptional induction through canonical GR binding sites is associated with reciprocal repression of distal TNF-regulated enhancers through a chromatin-dependent process, as evidenced by chromatin accessibility and motif displacement analysis. Simultaneously, transcriptional induction of key anti-inflammatory effectors is decoupled from primary repression through cooperation between GR and NF-kB at a subset of regulatory regions. Thus, glucocorticoids exert bimodal restraints on inflammation characterized by rapid primary transcriptional repression without local GR occupancy and secondary anti-inflammatory effects resulting from transcriptional cooperation between GR and NF-kB.

Transcriptional Responses to IFN-γ Require Mediator Kinase-Dependent Pause Release and Mechanistically Distinct CDK8 and CDK19 Functions.

Transcriptional responses to external stimuli remain poorly understood. Using global nuclear run-on followed by sequencing (GRO-seq) and precision nuclear run-on sequencing (PRO-seq), we show that CDK8 kinase activity promotes RNA polymerase II pause release in response to interferon-γ (IFN-γ), a universal cytokine involved in immunity and tumor surveillance. The Mediator kinase module contains CDK8 or CDK19, which are presumed to be functionally redundant. We implemented cortistatin A, chemical genetics, transcriptomics, and other methods to decouple their function while assessing enzymatic versus structural roles. Unexpectedly, CDK8 and CDK19 regulated different gene sets via distinct mechanisms. CDK8-dependent regulation required its kinase activity, whereas CDK19 governed IFN-γ responses through its scaffolding function (i.e., it was kinase independent). Accordingly, CDK8, not CDK19, phosphorylates the STAT1 transcription factor (TF) during IFN-γ stimulation, and CDK8 kinase inhibition blocked activation of JAK-STAT pathway TFs. Cytokines such as IFN-γ rapidly mobilize TFs to "reprogram" cellular transcription; our results implicate CDK8 and CDK19 as essential for this transcriptional reprogramming.

Identification and characterization of a novel anti-inflammatory lipid isolated from Mycobacterium vaccae, a soil-derived bacterium with immunoregulatory and stress resilience properties.

RATIONALE: Mycobacterium vaccae (NCTC 11659) is an environmental saprophytic bacterium with anti-inflammatory, immunoregulatory, and stress resilience properties. Previous studies have shown that whole, heat-killed preparations of M. vaccae prevent allergic airway inflammation in a murine model of allergic asthma. Recent studies also demonstrate that immunization with M. vaccae prevents stress-induced exaggeration of proinflammatory cytokine secretion from mesenteric lymph node cells stimulated ex vivo, prevents stress-induced exaggeration of chemically induced colitis in a model of inflammatory bowel disease, and prevents stress-induced anxiety-like defensive behavioral responses. Furthermore, immunization with M. vaccae induces anti-inflammatory responses in the brain and prevents stress-induced exaggeration of microglial priming. However, the molecular mechanisms underlying anti-inflammatory effects of M. vaccae are not known. OBJECTIVES: Our objective was to identify and characterize novel anti-inflammatory molecules from M. vaccae NCTC 11659. METHODS: We have purified and identified a unique anti-inflammatory triglyceride, 1,2,3-tri [Z-10-hexadecenoyl] glycerol, from M. vaccae and evaluated its effects in freshly isolated murine peritoneal macrophages. RESULTS: The free fatty acid form of 1,2,3-tri [Z-10-hexadecenoyl] glycerol, 10(Z)-hexadecenoic acid, decreased lipopolysaccharide-stimulated secretion of the proinflammatory cytokine IL-6 ex vivo. Meanwhile, next-generation RNA sequencing revealed that pretreatment with 10(Z)-hexadecenoic acid upregulated genes associated with peroxisome proliferator-activated receptor alpha (PPARα) signaling in lipopolysaccharide-stimulated macrophages, in association with a broad transcriptional repression of inflammatory markers. We confirmed using luciferase-based transfection assays that 10(Z)-hexadecenoic acid activated PPARα signaling, but not PPARγ, PPARδ, or retinoic acid receptor (RAR) α signaling. The effects of 10(Z)-hexadecenoic acid on lipopolysaccharide-stimulated secretion of IL-6 were prevented by PPARα antagonists and absent in PPARα-deficient mice. CONCLUSION: Future studies should evaluate the effects of 10(Z)-hexadecenoic acid on stress-induced exaggeration of peripheral inflammatory signaling, central neuroinflammatory signaling, and anxiety- and fear-related defensive behavioral responses.

Control of inducible gene expression links cohesin to hematopoietic progenitor self-renewal and differentiation.

Cohesin is important for 3D genome organization. Nevertheless, even the complete removal of cohesin has surprisingly little impact on steady-state gene transcription and enhancer activity. Here we show that cohesin is required for the core transcriptional response of primary macrophages to microbial signals, and for inducible enhancer activity that underpins inflammatory gene expression. Consistent with a role for inflammatory signals in promoting myeloid differentiation of hematopoietic stem and progenitor cells (HPSCs), cohesin mutations in HSPCs led to reduced inflammatory gene expression and increased resistance to differentiation-inducing inflammatory stimuli. These findings uncover an unexpected dependence of inducible gene expression on cohesin, link cohesin with myeloid differentiation, and may help explain the prevalence of cohesin mutations in human acute myeloid leukemia.

Detecting Differential Transcription Factor Activity from ATAC-Seq Data.

Transcription factors are managers of the cellular factory, and key components to many diseases. Many non-coding single nucleotide polymorphisms affect transcription factors, either by directly altering the protein or its functional activity at individual binding sites. Here we first briefly summarize high-throughput approaches to studying transcription factor activity. We then demonstrate, using published chromatin accessibility data (specifically ATAC-seq), that the genome-wide profile of TF recognition motifs relative to regions of open chromatin can determine the key transcription factor altered by a perturbation. Our method of determining which TFs are altered by a perturbation is simple, is quick to implement, and can be used when biological samples are limited. In the future, we envision that this method could be applied to determine which TFs show altered activity in response to a wide variety of drugs and diseases.

PI(3,5)P2 controls vacuole potassium transport to support cellular osmoregulation.

Lysosomes are dynamic organelles with critical roles in cellular physiology. The lysosomal signaling lipid phosphatidylinositol 3,5-bisphosphate (PI(3,5)P2) is a key regulator that has been implicated to control lysosome ion homeostasis, but the scope of ion transporters targeted by PI(3,5)P2 and the purpose of this regulation is not well understood. Through an unbiased screen in Saccharomyces cerevisiae, we identified loss-of-function mutations in the vacuolar H+-ATPase (V-ATPase) and in Vnx1, a vacuolar monovalent cation/proton antiporter, as suppressor mutations that relieve the growth defects and osmotic swelling of vacuoles (lysosomes) in yeast lacking PI(3,5)P2. We observed that depletion of PI(3,5)P2 synthesis in yeast causes a robust accumulation of multiple cations, most notably an ∼85 mM increase in the cellular concentration of potassium, a critical ion used by cells to regulate osmolarity. The accumulation of potassium and other cations in PI(3,5)P2-deficient yeast is relieved by mutations that inactivate Vnx1 or inactivate the V-ATPase and by mutations that increase the activity of a vacuolar cation export channel, Yvc1. Collectively, our data demonstrate that PI(3,5)P2 signaling orchestrates vacuole/lysosome cation transport to aid cellular osmoregulation.

Genome-wide dose-dependent inhibition of histone deacetylases studies reveal their roles in enhancer remodeling and suppression of oncogenic super-enhancers.

Histone deacetylase inhibitors (HDACIs) are known to alter gene expression by both up- and down-regulation of protein-coding genes in normal and cancer cells. However, the exact regulatory mechanisms of action remain uncharacterized. Here we investigated genome wide dose-dependent epigenetic and transcriptome changes in response to HDACI largazole in a transformed and a non-transformed cell line. Exposure to low nanomolar largazole concentrations (