ARID1A suppresses R-loop-mediated STING-type I interferon pathway activation of anti-tumor immunity.

Clinical trials have identified ARID1A mutations as enriched among patients who respond favorably to immune checkpoint blockade (ICB) in several solid tumor types independent of microsatellite instability. We show that ARID1A loss in murine models is sufficient to induce anti-tumor immune phenotypes observed in ARID1A mutant human cancers, including increased CD8+ T cell infiltration and cytolytic activity. ARID1A-deficient cancers upregulated an interferon (IFN) gene expression signature, the ARID1A-IFN signature, associated with increased R-loops and cytosolic single-stranded DNA (ssDNA). Overexpression of the R-loop resolving enzyme, RNASEH2B, or cytosolic DNase, TREX1, in ARID1A-deficient cells prevented cytosolic ssDNA accumulation and ARID1A-IFN gene upregulation. Further, the ARID1A-IFN signature and anti-tumor immunity were driven by STING-dependent type I IFN signaling, which was required for improved responsiveness of ARID1A mutant tumors to ICB treatment. These findings define a molecular mechanism underlying anti-tumor immunity in ARID1A mutant cancers.

Spatial CRISPR genomics identifies regulators of the tumor microenvironment.

While CRISPR screens are helping uncover genes regulating many cell-intrinsic processes, existing approaches are suboptimal for identifying extracellular gene functions, particularly in the tissue context. Here, we developed an approach for spatial functional genomics called Perturb-map. We applied Perturb-map to knock out dozens of genes in parallel in a mouse model of lung cancer and simultaneously assessed how each knockout influenced tumor growth, histopathology, and immune composition. Moreover, we paired Perturb-map and spatial transcriptomics for unbiased analysis of CRISPR-edited tumors. We found that in Tgfbr2 knockout tumors, the tumor microenvironment (TME) was converted to a fibro-mucinous state, and T cells excluded, concomitant with upregulated TGFß and TGFß-mediated fibroblast activation, indicating that TGFß-receptor loss on cancer cells increased TGFß bioavailability and its immunosuppressive effects on the TME. These studies establish Perturb-map for functional genomics within the tissue at single-cell resolution with spatial architecture preserved and provide insight into how TGFß responsiveness of cancer cells can affect the TME.

ZBTB46 coordinates angiogenesis and immunity to control tumor outcome.

Tumor angiogenesis and immunity show an inverse correlation in cancer progression and outcome1. Here, we report that ZBTB46, a repressive transcription factor and a widely accepted marker for classical dendritic cells (DCs)2,3, controls both tumor angiogenesis and immunity. Zbtb46 was downregulated in both DCs and endothelial cells by tumor-derived factors to facilitate robust tumor growth. Zbtb46 downregulation led to a hallmark pro-tumor microenvironment (TME), including dysfunctional vasculature and immunosuppressive conditions. Analysis of human cancer data revealed a similar association of low ZBTB46 expression with an immunosuppressive TME and a worse prognosis. In contrast, enforced Zbtb46 expression led to TME changes to restrict tumor growth. Mechanistically, Zbtb46-deficient endothelial cells were highly angiogenic, and Zbtb46-deficient bone marrow progenitors upregulated Cebpb and diverted the DC program to immunosuppressive myeloid lineage output, potentially explaining the myeloid lineage skewing phenomenon in cancer4. Conversely, enforced Zbtb46 expression normalized tumor vessels and, by suppressing Cebpb, skewed bone marrow precursors toward immunostimulatory myeloid lineage output, leading to an immune-hot TME. Remarkably, Zbtb46 mRNA treatment synergized with anti-PD1 immunotherapy to improve tumor management in preclinical models. These findings identify ZBTB46 as a critical factor for angiogenesis and for myeloid lineage skewing in cancer and suggest that maintaining its expression could have therapeutic benefits.

Conserved transcriptional connectivity of regulatory T cells in the tumor microenvironment informs new combination cancer therapy strategies.

While regulatory T (Treg) cells are traditionally viewed as professional suppressors of antigen presenting cells and effector T cells in both autoimmunity and cancer, recent findings of distinct Treg cell functions in tissue maintenance suggest that their regulatory purview extends to a wider range of cells and is broader than previously assumed. To elucidate tumoral Treg cell 'connectivity' to diverse tumor-supporting accessory cell types, we explored immediate early changes in their single-cell transcriptomes upon punctual Treg cell depletion in experimental lung cancer and injury-induced inflammation. Before any notable T cell activation and inflammation, fibroblasts, endothelial and myeloid cells exhibited pronounced changes in their gene expression in both cancer and injury settings. Factor analysis revealed shared Treg cell-dependent gene programs, foremost, prominent upregulation of VEGF and CCR2 signaling-related genes upon Treg cell deprivation in either setting, as well as in Treg cell-poor versus Treg cell-rich human lung adenocarcinomas. Accordingly, punctual Treg cell depletion combined with short-term VEGF blockade showed markedly improved control of PD-1 blockade-resistant lung adenocarcinoma progression in mice compared to the corresponding monotherapies, highlighting a promising factor-based querying approach to elucidating new rational combination treatments of solid organ cancers.

Identification of RIOK2 as a master regulator of human blood cell development.

Anemia is a major comorbidity in aging, chronic kidney and inflammatory diseases, and hematologic malignancies. However, the transcriptomic networks governing hematopoietic differentiation in blood cell development remain incompletely defined. Here we report that the atypical kinase RIOK2 (right open reading frame kinase 2) is a master transcription factor (TF) that not only drives erythroid differentiation, but also simultaneously suppresses megakaryopoiesis and myelopoiesis in primary human stem and progenitor cells. Our study reveals the previously uncharacterized winged helix-turn-helix DNA-binding domain and two transactivation domains of RIOK2 that are critical to regulate key hematopoietic TFs GATA1, GATA2, SPI1, RUNX3 and KLF1. This establishes RIOK2 as an integral component of the transcriptional regulatory network governing human hematopoietic differentiation. Importantly, RIOK2 mRNA expression significantly correlates with these TFs and other hematopoietic genes in myelodysplastic syndromes, acute myeloid leukemia and chronic kidney disease. Further investigation of RIOK2-mediated transcriptional pathways should yield therapeutic approaches to correct defective hematopoiesis in hematologic disorders.

PD-1 directed immunotherapy alters Tfh and humoral immune responses to seasonal influenza vaccine.

Anti-programmed death-1 (anti-PD-1) immunotherapy reinvigorates CD8 T cell responses in patients with cancer but PD-1 is also expressed by other immune cells, including follicular helper CD4 T cells (Tfh) which are involved in germinal centre responses. Little is known, however, about the effects of anti-PD-1 immunotherapy on noncancer immune responses in humans. To investigate this question, we examined the impact of anti-PD-1 immunotherapy on the Tfh-B cell axis responding to unrelated viral antigens. Following influenza vaccination, a subset of adults receiving anti-PD-1 had more robust circulating Tfh responses than adults not receiving immunotherapy. PD-1 pathway blockade resulted in transcriptional signatures of increased cellular proliferation in circulating Tfh and responding B cells compared with controls. These latter observations suggest an underlying change in the Tfh-B cell and germinal centre axis in a subset of immunotherapy patients. Together, these results demonstrate dynamic effects of anti-PD-1 therapy on influenza vaccine responses and highlight analytical vaccination as an approach that may reveal underlying immune predisposition to adverse events.

Dietary Intake Regulates the Circulating Inflammatory Monocyte Pool.

Caloric restriction is known to improve inflammatory and autoimmune diseases. However, the mechanisms by which reduced caloric intake modulates inflammation are poorly understood. Here we show that short-term fasting reduced monocyte metabolic and inflammatory activity and drastically reduced the number of circulating monocytes. Regulation of peripheral monocyte numbers was dependent on dietary glucose and protein levels. Specifically, we found that activation of the low-energy sensor 5'-AMP-activated protein kinase (AMPK) in hepatocytes and suppression of systemic CCL2 production by peroxisome proliferator-activator receptor alpha (PPARα) reduced monocyte mobilization from the bone marrow. Importantly, we show that fasting improves chronic inflammatory diseases without compromising monocyte emergency mobilization during acute infectious inflammation and tissue repair. These results reveal that caloric intake and liver energy sensors dictate the blood and tissue immune tone and link dietary habits to inflammatory disease outcome.

The cis-Regulatory Atlas of the Mouse Immune System.

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

Blockade of IL-22 signaling reverses erythroid dysfunction in stress-induced anemias.

Patients with myelodysplastic syndromes (MDSs) display severe anemia but the mechanisms underlying this phenotype are incompletely understood. Right open-reading-frame kinase 2 (RIOK2) encodes a protein kinase located at 5q15, a region frequently lost in patients with MDS del(5q). Here we show that hematopoietic cell-specific haploinsufficient deletion of Riok2 (Riok2f/+Vav1cre) led to reduced erythroid precursor frequency leading to anemia. Proteomic analysis of Riok2f/+Vav1cre erythroid precursors suggested immune system activation, and transcriptomic analysis revealed an increase in p53-dependent interleukin (IL)-22 in Riok2f/+Vav1cre CD4+ T cells (TH22). Further, we discovered that the IL-22 receptor, IL-22RA1, was unexpectedly present on erythroid precursors. Blockade of IL-22 signaling alleviated anemia not only in Riok2f/+Vav1cre mice but also in wild-type mice. Serum concentrations of IL-22 were increased in the subset of patients with del(5q) MDS as well as patients with anemia secondary to chronic kidney disease. This work reveals a possible therapeutic opportunity for reversing many stress-induced anemias by targeting IL-22 signaling.

A microRNA expression and regulatory element activity atlas of the mouse immune system.

To better define the control of immune system regulation, we generated an atlas of microRNA (miRNA) expression from 63 mouse immune cell populations and connected these signatures with assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation followed by sequencing (ChIP-seq) and nascent RNA profiles to establish a map of miRNA promoter and enhancer usage in immune cells. miRNA complexity was relatively low, with >90% of the miRNA compartment of each population comprising <75 miRNAs; however, each cell type had a unique miRNA signature. Integration of miRNA expression with chromatin accessibility revealed putative regulatory elements for differentially expressed miRNAs, including miR-21a, miR-146a and miR-223. The integrated maps suggest that many miRNAs utilize multiple promoters to reach high abundance and identified dominant and divergent miRNA regulatory elements between lineages and during development that may be used by clustered miRNAs, such as miR-99a/let-7c/miR-125b, to achieve distinct expression. These studies, with web-accessible data, help delineate the cis-regulatory elements controlling miRNA signatures of the immune system.

The Human Transcription Factors.

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

Protein Barcodes Enable High-Dimensional Single-Cell CRISPR Screens.

CRISPR pools are being widely employed to identify gene functions. However, current technology, which utilizes DNA as barcodes, permits limited phenotyping and bulk-cell resolution. To enable novel screening capabilities, we developed a barcoding system operating at the protein level. We synthesized modules encoding triplet combinations of linear epitopes to generate >100 unique protein barcodes (Pro-Codes). Pro-Code-expressing vectors were introduced into cells and analyzed by CyTOF mass cytometry. Using just 14 antibodies, we detected 364 Pro-Code populations; establishing the largest set of protein-based reporters. By pairing each Pro-Code with a different CRISPR, we simultaneously analyzed multiple phenotypic markers, including phospho-signaling, on dozens of knockouts. Pro-Code/CRISPR screens found two interferon-stimulated genes, the immunoproteasome component Psmb8 and a chaperone Rtp4, are important for antigen-dependent immune editing of cancer cells and identified Socs1 as a negative regulator of Pd-l1. The Pro-Code technology enables simultaneous high-dimensional protein-level phenotyping of 100s of genes with single-cell resolution.

A regulatory circuit controlled by extranuclear and nuclear retinoic acid receptor α determines T cell activation and function.

Ligation of retinoic acid receptor alpha (RARα) by RA promotes varied transcriptional programs associated with immune activation and tolerance, but genetic deletion approaches suggest the impact of RARα on TCR signaling. Here, we examined whether RARα would exert roles beyond transcriptional regulation. Specific deletion of the nuclear isoform of RARα revealed an RARα isoform in the cytoplasm of T cells. Extranuclear RARα was rapidly phosphorylated upon TCR stimulation and recruited to the TCR signalosome. RA interfered with extranuclear RARα signaling, causing suboptimal TCR activation while enhancing FOXP3+ regulatory T cell conversion. TCR activation induced the expression of CRABP2, which translocates RA to the nucleus. Deletion of Crabp2 led to increased RA in the cytoplasm and interfered with signalosome-RARα, resulting in impaired anti-pathogen immunity and suppressed autoimmune disease. Our findings underscore the significance of subcellular RA/RARα signaling in T cells and identify extranuclear RARα as a component of the TCR signalosome and a determinant of immune responses.

Human microglia maturation is underpinned by specific gene regulatory networks.

Microglia phenotypes are highly regulated by the brain environment, but the transcriptional networks that specify the maturation of human microglia are poorly understood. Here, we characterized stage-specific transcriptomes and epigenetic landscapes of fetal and postnatal human microglia and acquired corresponding data in induced pluripotent stem cell (iPSC)-derived microglia, in cerebral organoids, and following engraftment into humanized mice. Parallel development of computational approaches that considered transcription factor (TF) co-occurrence and enhancer activity allowed prediction of shared and state-specific gene regulatory networks associated with fetal and postnatal microglia. Additionally, many features of the human fetal-to-postnatal transition were recapitulated in a time-dependent manner following the engraftment of iPSC cells into humanized mice. These data and accompanying computational approaches will facilitate further efforts to elucidate mechanisms by which human microglia acquire stage- and disease-specific phenotypes.

Tridimensional Visualization and Analysis of Early Human Development.

Generating a precise cellular and molecular cartography of the human embryo is essential to our understanding of the mechanisms of organogenesis in normal and pathological conditions. Here, we have combined whole-mount immunostaining, 3DISCO clearing, and light-sheet imaging to start building a 3D cellular map of the human development during the first trimester of gestation. We provide high-resolution 3D images of the developing peripheral nervous, muscular, vascular, cardiopulmonary, and urogenital systems. We found that the adult-like pattern of skin innervation is established before the end of the first trimester, showing important intra- and inter-individual variations in nerve branches. We also present evidence for a differential vascularization of the male and female genital tracts concomitant with sex determination. This work paves the way for a cellular and molecular reference atlas of human cells, which will be of paramount importance to understanding human development in health and disease. PAPERCLIP.

WDR5 represents a therapeutically exploitable target for cancer stem cells in glioblastoma.

Glioblastomas (GBMs) are heterogeneous, treatment-resistant tumors driven by populations of cancer stem cells (CSCs). However, few molecular mechanisms critical for CSC population maintenance have been exploited for therapeutic development. We developed a spatially resolved loss-of-function screen in GBM patient-derived organoids to identify essential epigenetic regulators in the SOX2-enriched, therapy-resistant niche and identified WDR5 as indispensable for this population. WDR5 is a component of the WRAD complex, which promotes SET1 family-mediated Lys4 methylation of histone H3 (H3K4me), associated with positive regulation of transcription. In GBM CSCs, WDR5 inhibitors blocked WRAD complex assembly and reduced H3K4 trimethylation and expression of genes involved in CSC-relevant oncogenic pathways. H3K4me3 peaks lost with WDR5 inhibitor treatment occurred disproportionally on POU transcription factor motifs, including the POU5F1(OCT4)::SOX2 motif. Use of a SOX2/OCT4 reporter demonstrated that WDR5 inhibitor treatment diminished cells with high reporter activity. Furthermore, WDR5 inhibitor treatment and WDR5 knockdown altered the stem cell state, disrupting CSC in vitro growth and self-renewal, as well as in vivo tumor growth. These findings highlight the role of WDR5 and the WRAD complex in maintaining the CSC state and provide a rationale for therapeutic development of WDR5 inhibitors for GBM and other advanced cancers.

Mitochondrial protein import stress regulates the LC3 lipidation step of mitophagy through NLRX1 and RRBP1.

Protein import into mitochondria is a highly regulated process, yet how cells clear mitochondria undergoing dysfunctional protein import remains poorly characterized. Here we showed that mitochondrial protein import stress (MPIS) triggers localized LC3 lipidation. This arm of the mitophagy pathway occurs through the Nod-like receptor (NLR) protein NLRX1 while, surprisingly, without the engagement of the canonical mitophagy protein PINK1. Mitochondrial depolarization, which itself induces MPIS, also required NLRX1 for LC3 lipidation. While normally targeted to the mitochondrial matrix, cytosol-retained NLRX1 recruited RRBP1, a ribosome-binding transmembrane protein of the endoplasmic reticulum, which relocated to the mitochondrial vicinity during MPIS, and the NLRX1/RRBP1 complex in turn controlled the recruitment and lipidation of LC3. Furthermore, NLRX1 controlled skeletal muscle mitophagy in vivo and regulated endurance capacity during exercise. Thus, localization and lipidation of LC3 at the site of mitophagosome formation is a regulated step of mitophagy controlled by NLRX1/RRBP1 in response to MPIS.

Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

An atlas of genetic scores to predict multi-omic traits.

The use of omic modalities to dissect the molecular underpinnings of common diseases and traits is becoming increasingly common. But multi-omic traits can be genetically predicted, which enables highly cost-effective and powerful analyses for studies that do not have multi-omics1. Here we examine a large cohort (the INTERVAL study2; n = 50,000 participants) with extensive multi-omic data for plasma proteomics (SomaScan, n = 3,175; Olink, n = 4,822), plasma metabolomics (Metabolon HD4, n = 8,153), serum metabolomics (Nightingale, n = 37,359) and whole-blood Illumina RNA sequencing (n = 4,136), and use machine learning to train genetic scores for 17,227 molecular traits, including 10,521 that reach Bonferroni-adjusted significance. We evaluate the performance of genetic scores through external validation across cohorts of individuals of European, Asian and African American ancestries. In addition, we show the utility of these multi-omic genetic scores by quantifying the genetic control of biological pathways and by generating a synthetic multi-omic dataset of the UK Biobank3 to identify disease associations using a phenome-wide scan. We highlight a series of biological insights with regard to genetic mechanisms in metabolism and canonical pathway associations with disease; for example, JAK-STAT signalling and coronary atherosclerosis. Finally, we develop a portal ( https://www.omicspred.org/ ) to facilitate public access to all genetic scores and validation results, as well as to serve as a platform for future extensions and enhancements of multi-omic genetic scores.