Identification of the novel FOXP3-dependent Treg cell transcription factor MEOX1 by high-dimensional analysis of human CD4+ T cells.

CD4+ T cells play a central role in the adaptive immune response through their capacity to activate, support and control other immune cells. Although these cells have become the focus of intense research, a comprehensive understanding of the underlying regulatory networks that orchestrate CD4+ T cell function and activation is still incomplete. Here, we analyzed a large transcriptomic dataset consisting of 48 different human CD4+ T cell conditions. By performing reverse network engineering, we identified six common denominators of CD4+ T cell functionality (CREB1, E2F3, AHR, STAT1, NFAT5 and NFATC3). Moreover, we also analyzed condition-specific genes which led us to the identification of the transcription factor MEOX1 in Treg cells. Expression of MEOX1 was comparable to FOXP3 in Treg cells and can be upregulated by IL-2. Epigenetic analyses revealed a permissive epigenetic landscape for MEOX1 solely in Treg cells. Knockdown of MEOX1 in Treg cells revealed a profound impact on downstream gene expression programs and Treg cell suppressive capacity. These findings in the context of CD4+ T cells contribute to a better understanding of the transcriptional networks and biological mechanisms controlling CD4+ T cell functionality, which opens new avenues for future therapeutic strategies.

Parallel recovery of chromatin accessibility and gene expression dynamics from frozen human regulatory T cells.

Epigenetic features such as DNA accessibility dictate transcriptional regulation in a cell type- and cell state- specific manner, and mapping this in health vs. disease in clinically relevant material is opening the door to new mechanistic insights and new targets for therapy. Assay for Transposase Accessible Chromatin Sequencing (ATAC-seq) allows chromatin accessibility profiling from low cell input, making it tractable on rare cell populations, such as regulatory T (Treg) cells. However, little is known about the compatibility of the assay with cryopreserved rare cell populations. Here we demonstrate the robustness of an ATAC-seq protocol comparing primary Treg cells recovered from fresh or cryopreserved PBMC samples, in the steady state and in response to stimulation. We extend this method to explore the feasibility of conducting simultaneous quantitation of chromatin accessibility and transcriptome from a single aliquot of 50,000 cryopreserved Treg cells. Profiling of chromatin accessibility and gene expression in parallel within the same pool of cells controls for cellular heterogeneity and is particularly beneficial when constrained by limited input material. Overall, we observed a high correlation of accessibility patterns and transcription factor dynamics between fresh and cryopreserved samples. Furthermore, highly similar transcriptomic profiles were obtained from whole cells and from the supernatants recovered from ATAC-seq reactions. We highlight the feasibility of applying these techniques to profile the epigenomic landscape of cells recovered from cryopreservation biorepositories.

Single cell sequencing reveals endothelial plasticity with transient mesenchymal activation after myocardial infarction.

Endothelial cells play a critical role in the adaptation of tissues to injury. Tissue ischemia induced by infarction leads to profound changes in endothelial cell functions and can induce transition to a mesenchymal state. Here we explore the kinetics and individual cellular responses of endothelial cells after myocardial infarction by using single cell RNA sequencing. This study demonstrates a time dependent switch in endothelial cell proliferation and inflammation associated with transient changes in metabolic gene signatures. Trajectory analysis reveals that the majority of endothelial cells 3 to 7 days after myocardial infarction acquire a transient state, characterized by mesenchymal gene expression, which returns to baseline 14 days after injury. Lineage tracing, using the Cdh5-CreERT2;mT/mG mice followed by single cell RNA sequencing, confirms the transient mesenchymal transition and reveals additional hypoxic and inflammatory signatures of endothelial cells during early and late states after injury. These data suggest that endothelial cells undergo a transient mes-enchymal activation concomitant with a metabolic adaptation within the first days after myocardial infarction but do not acquire a long-term mesenchymal fate. This mesenchymal activation may facilitate endothelial cell migration and clonal expansion to regenerate the vascular network.

Two populations of self-maintaining monocyte-independent macrophages exist in adult epididymis and testis.

Macrophages are the principal immune cells of the epididymis and testis, but their origins, heterogeneity, development, and maintenance are not well understood. Here, we describe distinct populations of epididymal and testicular macrophages that display an organ-specific cellular identity. Combining in vivo fate-mapping, chimeric and parabiotic mouse models with in-depth cellular analyses, we found that CD64hiMHCIIlo and CD64loMHCIIhi macrophage populations of epididymis and testis arise sequentially from yolk sac erythro-myeloid progenitors, embryonic hematopoiesis, and nascent neonatal monocytes. While monocytes were the major developmental source of both epididymal and testicular macrophages, both populations self-maintain in the steady-state independent of bone marrow hematopoietic precursors. However, after radiation-induced macrophage ablation or during infection, bone marrow-derived circulating monocytes are recruited to the epididymis and testis, giving rise to inflammatory macrophages that promote tissue damage. These results define the layered ontogeny, maintenance and inflammatory response of macrophage populations in the male reproductive organs.

ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation.

While footprinting analysis of ATAC-seq data can theoretically enable investigation of transcription factor (TF) binding, the lack of a computational tool able to conduct different levels of footprinting analysis has so-far hindered the widespread application of this method. Here we present TOBIAS, a comprehensive, accurate, and fast footprinting framework enabling genome-wide investigation of TF binding dynamics for hundreds of TFs simultaneously. We validate TOBIAS using paired ATAC-seq and ChIP-seq data, and find that TOBIAS outperforms existing methods for bias correction and footprinting. As a proof-of-concept, we illustrate how TOBIAS can unveil complex TF dynamics during zygotic genome activation in both humans and mice, and propose how zygotic Dux activates cascades of TFs, binds to repeat elements and induces expression of novel genetic elements.

BCG Vaccination in Humans Elicits Trained Immunity via the Hematopoietic Progenitor Compartment.

Induction of trained immunity by Bacille-Calmette-Guérin (BCG) vaccination mediates beneficial heterologous effects, but the mechanisms underlying its persistence and magnitude remain elusive. In this study, we show that BCG vaccination in healthy human volunteers induces a persistent transcriptional program connected to myeloid cell development and function within the hematopoietic stem and progenitor cell (HSPC) compartment in the bone marrow. We identify hepatic nuclear factor (HNF) family members 1a and b as crucial regulators of this transcriptional shift. These findings are corroborated by higher granulocyte numbers in BCG-vaccinated infants, HNF1 SNP variants that correlate with trained immunity, and elevated serum concentrations of the HNF1 target alpha-1 antitrypsin. Additionally, transcriptomic HSPC remodeling was epigenetically conveyed to peripheral CD14+ monocytes, displaying an activated transcriptional signature three months after BCG vaccination. Taken together, transcriptomic, epigenomic, and functional reprogramming of HSPCs and peripheral monocytes is a hallmark of BCG-induced trained immunity in humans.

Cxcr4 distinguishes HSC-derived monocytes from microglia and reveals monocyte immune responses to experimental stroke.

Monocyte-derived and tissue-resident macrophages are ontogenetically distinct components of the innate immune system. Assessment of their respective functions in pathology is complicated by changes to the macrophage phenotype during inflammation. Here we find that Cxcr4-CreER enables permanent genetic labeling of hematopoietic stem cells (HSCs) and distinguishes HSC-derived monocytes from microglia and other tissue-resident macrophages. By combining Cxcr4-CreER-mediated lineage tracing with Cxcr4 inhibition or conditional Cxcr4 ablation in photothrombotic stroke, we find that Cxcr4 promotes initial monocyte infiltration and subsequent territorial restriction of monocyte-derived macrophages to infarct tissue. After transient focal ischemia, Cxcr4 deficiency reduces monocyte infiltration and blunts the expression of pattern recognition and defense response genes in monocyte-derived macrophages. This is associated with an altered microglial response and deteriorated outcomes. Thus, Cxcr4 is essential for an innate-immune-system-mediated defense response after cerebral ischemia. We further propose Cxcr4-CreER as a universal tool to study functions of HSC-derived cells.

Scalable Prediction of Acute Myeloid Leukemia Using High-Dimensional Machine Learning and Blood Transcriptomics.

Acute myeloid leukemia (AML) is a severe, mostly fatal hematopoietic malignancy. We were interested in whether transcriptomic-based machine learning could predict AML status without requiring expert input. Using 12,029 samples from 105 different studies, we present a large-scale study of machine learning-based prediction of AML in which we address key questions relating to the combination of machine learning and transcriptomics and their practical use. We find data-driven, high-dimensional approaches-in which multivariate signatures are learned directly from genome-wide data with no prior knowledge-to be accurate and robust. Importantly, these approaches are highly scalable with low marginal cost, essentially matching human expert annotation in a near-automated workflow. Our results support the notion that transcriptomics combined with machine learning could be used as part of an integrated -omics approach wherein risk prediction, differential diagnosis, and subclassification of AML are achieved by genomics while diagnosis could be assisted by transcriptomic-based machine learning.

Transcriptional Signature Derived from Murine Tumor-Associated Macrophages Correlates with Poor Outcome in Breast Cancer Patients.

Tumor-associated macrophages (TAMs) are frequently the most abundant immune cells in cancers and are associated with poor survival. Here, we generated TAM molecular signatures from K14cre;Cdh1flox/flox;Trp53flox/flox (KEP) and MMTV-NeuT (NeuT) transgenic mice that resemble human invasive lobular carcinoma (ILC) and HER2+ tumors, respectively. Determination of TAM-specific signatures requires comparison with healthy mammary tissue macrophages to avoid overestimation of gene expression differences. TAMs from the two models feature a distinct transcriptomic profile, suggesting that the cancer subtype dictates their phenotype. The KEP-derived signature reliably correlates with poor overall survival in ILC but not in triple-negative breast cancer patients, indicating that translation of murine TAM signatures to patients is cancer subtype dependent. Collectively, we show that a transgenic mouse tumor model can yield a TAM signature relevant for human breast cancer outcome prognosis and provide a generalizable strategy for determining and applying immune cell signatures provided the murine model reflects the human disease.

Enzymatic Activity of HPGD in Treg Cells Suppresses Tconv Cells to Maintain Adipose Tissue Homeostasis and Prevent Metabolic Dysfunction.

Regulatory T cells (Treg cells) are important for preventing autoimmunity and maintaining tissue homeostasis, but whether Treg cells can adopt tissue- or immune-context-specific suppressive mechanisms is unclear. Here, we found that the enzyme hydroxyprostaglandin dehydrogenase (HPGD), which catabolizes prostaglandin E2 (PGE2) into the metabolite 15-keto PGE2, was highly expressed in Treg cells, particularly those in visceral adipose tissue (VAT). Nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ)-induced HPGD expression in VAT Treg cells, and consequential Treg-cell-mediated generation of 15-keto PGE2 suppressed conventional T cell activation and proliferation. Conditional deletion of Hpgd in mouse Treg cells resulted in the accumulation of functionally impaired Treg cells specifically in VAT, causing local inflammation and systemic insulin resistance. Consistent with this mechanism, humans with type 2 diabetes showed decreased HPGD expression in Treg cells. These data indicate that HPGD-mediated suppression is a tissue- and context-dependent suppressive mechanism used by Treg cells to maintain adipose tissue homeostasis.

Ceramide Synthase Schlank Is a Transcriptional Regulator Adapting Gene Expression to Energy Requirements.

Maintenance of metabolic homeostasis requires adaption of gene regulation to the cellular energy state via transcriptional regulators. Here, we identify a role of ceramide synthase (CerS) Schlank, a multiple transmembrane protein containing a catalytic lag1p motif and a homeodomain, which is poorly studied in CerSs, as a transcriptional regulator. ChIP experiments show that it binds promoter regions of lipases lipase3 and magro via its homeodomain. Mutation of nuclear localization site 2 (NLS2) within the homeodomain leads to loss of DNA binding and deregulated gene expression, and NLS2 mutants can no longer adjust the transcriptional response to changing lipid levels. This mechanism is conserved in mammalian CerS2 and emphasizes the importance of the CerS protein rather than ceramide synthesis. This study demonstrates a double role of CerS Schlank as an enzyme and a transcriptional regulator, sensing lipid levels and transducing the information to the level of gene expression.

Western Diet Triggers NLRP3-Dependent Innate Immune Reprogramming.

Long-term epigenetic reprogramming of innate immune cells in response to microbes, also termed "trained immunity," causes prolonged altered cellular functionality to protect from secondary infections. Here, we investigated whether sterile triggers of inflammation induce trained immunity and thereby influence innate immune responses. Western diet (WD) feeding of Ldlr-/- mice induced systemic inflammation, which was undetectable in serum soon after mice were shifted back to a chow diet (CD). In contrast, myeloid cell responses toward innate stimuli remained broadly augmented. WD-induced transcriptomic and epigenomic reprogramming of myeloid progenitor cells led to increased proliferation and enhanced innate immune responses. Quantitative trait locus (QTL) analysis in human monocytes trained with oxidized low-density lipoprotein (oxLDL) and stimulated with lipopolysaccharide (LPS) suggested inflammasome-mediated trained immunity. Consistently, Nlrp3-/-/Ldlr-/- mice lacked WD-induced systemic inflammation, myeloid progenitor proliferation, and reprogramming. Hence, NLRP3 mediates trained immunity following WD and could thereby mediate the potentially deleterious effects of trained immunity in inflammatory diseases.

Cellular Differentiation of Human Monocytes Is Regulated by Time-Dependent Interleukin-4 Signaling and the Transcriptional Regulator NCOR2.

Human in vitro generated monocyte-derived dendritic cells (moDCs) and macrophages are used clinically, e.g., to induce immunity against cancer. However, their physiological counterparts, ontogeny, transcriptional regulation, and heterogeneity remains largely unknown, hampering their clinical use. High-dimensional techniques were used to elucidate transcriptional, phenotypic, and functional differences between human in vivo and in vitro generated mononuclear phagocytes to facilitate their full potential in the clinic. We demonstrate that monocytes differentiated by macrophage colony-stimulating factor (M-CSF) or granulocyte macrophage colony-stimulating factor (GM-CSF) resembled in vivo inflammatory macrophages, while moDCs resembled in vivo inflammatory DCs. Moreover, differentiated monocytes presented with profound transcriptomic, phenotypic, and functional differences. Monocytes integrated GM-CSF and IL-4 stimulation combinatorically and temporally, resulting in a mode- and time-dependent differentiation relying on NCOR2. Finally, moDCs are phenotypically heterogeneous and therefore necessitate the use of high-dimensional phenotyping to open new possibilities for better clinical tailoring of these cellular therapies.

Generation and functional characterization of MDSC-like cells.

Myeloid-derived suppressor cells (MDSC) are critical in regulating immune responses by suppressing antigen presenting cells (APC) and T cells. We previously observed that incubation of peripheral blood monocytes with interleukin (IL)-10 during their differentiation to monocyte-derived dendritic cells (moDCs) results in the generation of an APC population with a CD14+HLA-DRlowphenotype (IL-10-APC) with reduced stimulatory capacity similar to human MDSC. Co-incubation experiments now revealed that the addition of IL-10-APC to moDC caused a reduction of DC-induced T-cell proliferation, of the expression of maturation markers, and of secreted cytokines and chemokines such as TNF-α, IL-6, MIP-1α and Rantes. Addition of IL-10-APC increased the immunosuppressive molecule osteoactivin and its corresponding receptor syndecan-4 on moDC. Moreover, CD14+HLA-DRlow MDSC isolated from healthy donors expressed high levels of osteoactivin, which was even further upregulated by the auxiliary addition of IL-10. Using transcriptome analysis, we identified a set of molecules and pathways mediating these effects. In addition, we found that IL-10-APC as well as human isolated MDSC expressed higher levels of programmed death (PD)-1, PD-ligand-1 (PD-L1), glucocorticoid-induced-tumor-necrosis-factor-receptor-related-protein (GITR) and GITR-ligand. Inhibition of osteoactivin, syndecan-4, PD-1 or PD-L1 on MDSC by using blocking antibodies restored the stimulatory capacity of DC in co-incubation experiments. Activation of MDSC with Dectin-1 ligand curdlan reduced the expression of osteoactivin and PD-L1. Our results demonstrate that osteoactivin/syndecan-4 and PD-/PD-L1 are key molecules that are profoundly involved in the inhibitory effects of MDSC on DC function and might be promising tools for clinical application.

The transcriptional regulator network of human inflammatory macrophages is defined by open chromatin.

Differentiation of inflammatory macrophages from monocytes is characterized by an orderly integration of epigenetic and transcriptional regulatory mechanisms guided by lineage-determining transcription factors such as PU.1. Further activation of macrophages leads to a stimulus- or microenvironment-specific signal integration with subsequent transcriptional control established by the action of tissue- or signal-associated transcription factors. Here, we assess four histone modifications during human macrophage activation and integrate this information with the gene expression data from 28 different macrophage activation conditions in combination with GM-CSF. Bioinformatically, for inflammatory macrophages we define a unique network of transcriptional and epigenetic regulators (TRs), which was characterized by accessible promoters independent of the activation signal. In contrast to the general accessibility of promoters of TRs, mRNA expression of central TRs belonging to the TR network displayed stimulus-specific expression patterns, indicating a second level of transcriptional regulation beyond epigenetic chromatin changes. In contrast, stringent integration of epigenetic and transcriptional regulation was observed in networks of TRs established from somatic tissues and tissue macrophages. In these networks, clusters of TRs with permissive histone marks were associated with high gene expression whereas clusters with repressive chromatin marks were associated with absent gene expression. Collectively, these results support that macrophage activation during inflammation in contrast to lineage determination is mainly regulated transcriptionally by a pre-defined TR network.

ATF3 Is a Key Regulator of Macrophage IFN Responses.

Cytokines and IFNs downstream of innate immune pathways are critical for mounting an appropriate immune response to microbial infection. However, the expression of these inflammatory mediators is tightly regulated, as uncontrolled production can result in tissue damage and lead to chronic inflammatory conditions and autoimmune diseases. Activating transcription factor 3 (ATF3) is an important transcriptional modulator that limits the inflammatory response by controlling the expression of a number of cytokines and chemokines. However, its role in modulating IFN responses remains poorly defined. In this study, we demonstrate that ATF3 expression in macrophages is necessary for governing basal IFN-ß expression, as well as the magnitude of IFN-ß cytokine production following activation of innate immune receptors. We found that ATF3 acted as a transcriptional repressor and regulated IFN-ß via direct binding to a previously unidentified specific regulatory site distal to the Ifnb1 promoter. Additionally, we observed that ATF3 itself is a type I IFN-inducible gene, and that ATF3 further modulates the expression of a subset of inflammatory genes downstream of IFN signaling, suggesting it constitutes a key component of an IFN negative feedback loop. Consistent with this, macrophages deficient in Atf3 showed enhanced viral clearance in lymphocytic choriomeningitis virus and vesicular stomatitis virus infection models. Our study therefore demonstrates an important role for ATF3 in modulating IFN responses in macrophages by controlling basal and inducible levels of IFNß, as well as the expression of genes downstream of IFN signaling.

Co-existence of intact stemness and priming of neural differentiation programs in mES cells lacking Trim71.

Regulatory networks for differentiation and pluripotency in embryonic stem (ES) cells have long been suggested to be mutually exclusive. However, with the identification of many new components of these networks ranging from epigenetic, transcriptional, and translational to even post-translational mechanisms, the cellular states of pluripotency and early differentiation might not be strictly bi-modal, but differentiating stem cells appear to go through phases of simultaneous expression of stemness and differentiation genes. Translational regulators such as RNA binding proteins (RBPs) and micro RNAs (miRNAs) might be prime candidates for guiding a cell from pluripotency to differentiation. Using Trim71, one of two members of the Tripartite motif (Trim) protein family with RNA binding activity expressed in murine ES cells, we demonstrate that Trim71 is not involved in regulatory networks of pluripotency but regulates neural differentiation. Loss of Trim71 in mES cells leaves stemness and self-maintenance of these cells intact, but many genes required for neural development are up-regulated at the same time. Concordantly, Trim71(-/-) mES show increased neural marker expression following treatment with retinoic acid. Our findings strongly suggest that Trim71 keeps priming steps of differentiation in check, which do not pre-require a loss of the pluripotency network in ES cells.

Genome of the haloarchaeon Natronomonas moolapensis, a neutrophilic member of a previously haloalkaliphilic genus.

The genus Natronomonas contains two species, one haloalkaliphile (N. pharaonis) and one neutrophile (N. moolapensis). Here, we report the genome sequence of N. moolapensis strain 8.8.11. The overall genome properties are similar for the two species. Only the neutrophile contains bacteriorhodopsin and a membrane glycolipid.

The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions.

Watermelon, Citrullus lanatus, is an important cucurbit crop grown throughout the world. Here we report a high-quality draft genome sequence of the east Asia watermelon cultivar 97103 (2n = 2× = 22) containing 23,440 predicted protein-coding genes. Comparative genomics analysis provided an evolutionary scenario for the origin of the 11 watermelon chromosomes derived from a 7-chromosome paleohexaploid eudicot ancestor. Resequencing of 20 watermelon accessions representing three different C. lanatus subspecies produced numerous haplotypes and identified the extent of genetic diversity and population structure of watermelon germplasm. Genomic regions that were preferentially selected during domestication were identified. Many disease-resistance genes were also found to be lost during domestication. In addition, integrative genomic and transcriptomic analyses yielded important insights into aspects of phloem-based vascular signaling in common between watermelon and cucumber and identified genes crucial to valuable fruit-quality traits, including sugar accumulation and citrulline metabolism.

Haloquadratum walsbyi: limited diversity in a global pond.

BACKGROUND: Haloquadratum walsbyi commonly dominates the microbial flora of hypersaline waters. Its cells are extremely fragile squares requiring >14%(w/v) salt for growth, properties that should limit its dispersal and promote geographical isolation and divergence. To assess this, the genome sequences of two isolates recovered from sites at near maximum distance on Earth, were compared. PRINCIPAL FINDINGS: Both chromosomes are 3.1 MB in size, and 84% of each sequence was highly similar to the other (98.6% identity), comprising the core sequence. ORFs of this shared sequence were completely synteneic (conserved in genomic orientation and order), without inversion or rearrangement. Strain-specific insertions/deletions could be precisely mapped, often allowing the genetic events to be inferred. Many inferred deletions were associated with short direct repeats (4-20 bp). Deletion-coupled insertions are frequent, producing different sequences at identical positions. In cases where the inserted and deleted sequences are homologous, this leads to variant genes in a common synteneic background (as already described by others). Cas/CRISPR systems are present in C23(T) but have been lost in HBSQ001 except for a few spacer remnants. Numerous types of mobile genetic elements occur in both strains, most of which appear to be active, and with some specifically targetting others. Strain C23(T) carries two ∼6 kb plasmids that show similarity to halovirus His1 and to sequences nearby halovirus/plasmid gene clusters commonly found in haloarchaea. CONCLUSIONS: Deletion-coupled insertions show that Hqr. walsbyi evolves by uptake and precise integration of foreign DNA, probably originating from close relatives. Change is also driven by mobile genetic elements but these do not by themselves explain the atypically low gene coding density found in this species. The remarkable genome conservation despite the presence of active systems for genome rearrangement implies both an efficient global dispersal system, and a high selective fitness for this species.