Single-cell chromatin accessibility and transposable element landscapes reveal shared features of tissue-residing immune cells.

Tissue adaptation is required for regulatory T (Treg) cell function within organs. Whether this program shares aspects with other tissue-localized immune populations is unclear. Here, we analyzed single-cell chromatin accessibility data, including the transposable element (TE) landscape of CD45+ immune cells from colon, skin, adipose tissue, and spleen. We identified features of organ-specific tissue adaptation across different immune cells. Focusing on tissue Treg cells, we found conservation of the Treg tissue adaptation program in other tissue-localized immune cells, such as amphiregulin-producing T helper (Th)17 cells. Accessible TEs can act as regulatory elements, but their contribution to tissue adaptation is not understood. TE landscape analysis revealed an enrichment of specific transcription factor binding motifs in TE regions within accessible chromatin peaks. TEs, specifically from the LTR family, were located in enhancer regions and associated with tissue adaptation. These findings broaden our understanding of immune tissue residency and provide an important step toward organ-specific immune interventions.

The effector program of human CD8 T cells supports tissue remodeling.

CD8 T lymphocytes are classically viewed as cytotoxic T cells. Whether human CD8 T cells can, in parallel, induce a tissue regeneration program is poorly understood. Here, antigen-specific assay systems revealed that human CD8 T cells not only mediated cytotoxicity but also promoted tissue remodeling. Activated CD8 T cells could produce the epidermal growth factor receptor (EGFR)-ligand amphiregulin (AREG) and sensitize epithelial cells for enhanced regeneration potential. Blocking the EGFR or the effector cytokines IFN-γ and TNF could inhibit tissue remodeling. This regenerative program enhanced tumor spheroid and stem cell-mediated organoid growth. Using single-cell gene expression analysis, we identified an AREG+, tissue-resident CD8 T cell population in skin and adipose tissue from patients undergoing abdominal wall or abdominoplasty surgery. These tissue-resident CD8 T cells showed a strong TCR clonal relation to blood PD1+TIGIT+ CD8 T cells with tissue remodeling abilities. These findings may help to understand the complex CD8 biology in tumors and could become relevant for the design of therapeutic T cell products.

MCT4 blockade increases the efficacy of immune checkpoint blockade.

BACKGROUND & AIMS: Intratumoral lactate accumulation and acidosis impair T-cell function and antitumor immunity. Interestingly, expression of the lactate transporter monocarboxylate transporter (MCT) 4, but not MCT1, turned out to be prognostic for the survival of patients with rectal cancer, indicating that single MCT4 blockade might be a promising strategy to overcome glycolysis-related therapy resistance. METHODS: To determine whether blockade of MCT4 alone is sufficient to improve the efficacy of immune checkpoint blockade (ICB) therapy, we examined the effects of the selective MCT1 inhibitor AZD3965 and a novel MCT4 inhibitor in a colorectal carcinoma (CRC) tumor spheroid model co-cultured with blood leukocytes in vitro and the MC38 murine CRC model in vivo in combination with an antibody against programmed cell death ligand-1(PD-L1). RESULTS: Inhibition of MCT4 was sufficient to reduce lactate efflux in three-dimensional (3D) CRC spheroids but not in two-dimensional cell-cultures. Co-administration of the MCT4 inhibitor and ICB augmented immune cell infiltration, T-cell function and decreased CRC spheroid viability in a 3D co-culture model of human CRC spheroids with blood leukocytes. Accordingly, combination of MCT4 and ICB increased intratumoral pH, improved leukocyte infiltration and T-cell activation, delayed tumor growth, and prolonged survival in vivo. MCT1 inhibition exerted no further beneficial impact. CONCLUSIONS: These findings demonstrate that single MCT4 inhibition represents a novel therapeutic approach to reverse lactic-acid driven immunosuppression and might be suitable to improve ICB efficacy.

Comparative analysis of genome-scale, base-resolution DNA methylation profiles across 580 animal species.

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.

Integrated single-cell profiling dissects cell-state-specific enhancer landscapes of human tumor-infiltrating CD8+ T cells.

Despite extensive studies on the chromatin landscape of exhausted T cells, the transcriptional wiring underlying the heterogeneous functional and dysfunctional states of human tumor-infiltrating lymphocytes (TILs) is incompletely understood. Here, we identify gene-regulatory landscapes in a wide breadth of functional and dysfunctional CD8+ TIL states covering four cancer entities using single-cell chromatin profiling. We map enhancer-promoter interactions in human TILs by integrating single-cell chromatin accessibility with single-cell RNA-seq data from tumor-entity-matching samples and prioritize cell-state-specific genes by super-enhancer analysis. Besides revealing entity-specific chromatin remodeling in exhausted TILs, our analyses identify a common chromatin trajectory to TIL dysfunction and determine key enhancers, transcriptional regulators, and deregulated genes involved in this process. Finally, we validate enhancer regulation at immunotherapeutically relevant loci by targeting non-coding regulatory elements with potent CRISPR activators and repressors. In summary, our study provides a framework for understanding and manipulating cell-state-specific gene-regulatory cues from human tumor-infiltrating lymphocytes.

Salt-inducible kinase 3 protects tumor cells from cytotoxic T-cell attack by promoting TNF-induced NF-κB activation.

BACKGROUND: Cancer immunotherapeutic strategies showed unprecedented results in the clinic. However, many patients do not respond to immuno-oncological treatments due to the occurrence of a plethora of immunological obstacles, including tumor intrinsic mechanisms of resistance to cytotoxic T-cell (TC) attack. Thus, a deeper understanding of these mechanisms is needed to develop successful immunotherapies. METHODS: To identify novel genes that protect tumor cells from effective TC-mediated cytotoxicity, we performed a genetic screening in pancreatic cancer cells challenged with tumor-infiltrating lymphocytes and antigen-specific TCs. RESULTS: The screening revealed 108 potential genes that protected tumor cells from TC attack. Among them, salt-inducible kinase 3 (SIK3) was one of the strongest hits identified in the screening. Both genetic and pharmacological inhibitions of SIK3 in tumor cells dramatically increased TC-mediated cytotoxicity in several in vitro coculture models, using different sources of tumor and TCs. Consistently, adoptive TC transfer of TILs led to tumor growth inhibition of SIK3-depleted cancer cells in vivo. Mechanistic analysis revealed that SIK3 rendered tumor cells susceptible to tumor necrosis factor (TNF) secreted by tumor-activated TCs. SIK3 promoted nuclear factor kappa B (NF-κB) nuclear translocation and inhibited caspase-8 and caspase-9 after TNF stimulation. Chromatin accessibility and transcriptome analyses showed that SIK3 knockdown profoundly impaired the expression of prosurvival genes under the TNF-NF-κB axis. TNF stimulation led to SIK3-dependent phosphorylation of the NF-κB upstream regulators inhibitory-κB kinase and NF-kappa-B inhibitor alpha on the one side, and to inhibition of histone deacetylase 4 on the other side, thus sustaining NF-κB activation and nuclear stabilization. A SIK3-dependent gene signature of TNF-mediated NF-κB activation was found in a majority of pancreatic cancers where it correlated with increased cytotoxic TC activity and poor prognosis. CONCLUSION: Our data reveal an abundant molecular mechanism that protects tumor cells from cytotoxic TC attack and demonstrate that pharmacological inhibition of this pathway is feasible.

Single-cell chromatin accessibility landscape identifies tissue repair program in human regulatory T cells.

Murine regulatory T (Treg) cells in tissues promote tissue homeostasis and regeneration. We sought to identify features that characterize human Treg cells with these functions in healthy tissues. Single-cell chromatin accessibility profiles of murine and human tissue Treg cells defined a conserved, microbiota-independent tissue-repair Treg signature with a prevailing footprint of the transcription factor BATF. This signature, combined with gene expression profiling and TCR fate mapping, identified a population of tissue-like Treg cells in human peripheral blood that expressed BATF, chemokine receptor CCR8 and HLA-DR. Human BATF+CCR8+ Treg cells from normal skin and adipose tissue shared features with nonlymphoid T follicular helper-like (Tfh-like) cells, and induction of a Tfh-like differentiation program in naive human Treg cells partially recapitulated tissue Treg regenerative characteristics, including wound healing potential. Human BATF+CCR8+ Treg cells from healthy tissue share features with tumor-resident Treg cells, highlighting the importance of understanding the context-specific functions of these cells.

Cancer-associated cells release citrate to support tumour metastatic progression.

Citrate is important for lipid synthesis and epigenetic regulation in addition to ATP production. We have previously reported that cancer cells import extracellular citrate via the pmCiC transporter to support their metabolism. Here, we show for the first time that citrate is supplied to cancer by cancer-associated stroma (CAS) and also that citrate synthesis and release is one of the latter's major metabolic tasks. Citrate release from CAS is controlled by cancer cells through cross-cellular communication. The availability of citrate from CAS regulated the cytokine profile, metabolism and features of cellular invasion. Moreover, citrate released by CAS is involved in inducing cancer progression especially enhancing invasiveness and organ colonisation. In line with the in vitro observations, we show that depriving cancer cells of citrate using gluconate, a specific inhibitor of pmCiC, significantly reduced the growth and metastatic spread of human pancreatic cancer cells in vivo and muted stromal activation and angiogenesis. We conclude that citrate is supplied to tumour cells by CAS and citrate uptake plays a significant role in cancer metastatic progression.

Targeting clinical epigenetic reprogramming for chemoprevention of metabolic and viral hepatocellular carcinoma.

OBJECTIVE: Hepatocellular carcinoma (HCC) is the fastest-growing cause of cancer-related mortality with chronic viral hepatitis and non-alcoholic steatohepatitis (NASH) as major aetiologies. Treatment options for HCC are unsatisfactory and chemopreventive approaches are absent. Chronic hepatitis C (CHC) results in epigenetic alterations driving HCC risk and persisting following cure. Here, we aimed to investigate epigenetic modifications as targets for liver cancer chemoprevention. DESIGN: Liver tissues from patients with NASH and CHC were analysed by ChIP-Seq (H3K27ac) and RNA-Seq. The liver disease-specific epigenetic and transcriptional reprogramming in patients was modelled in a liver cell culture system. Perturbation studies combined with a targeted small molecule screen followed by in vivo and ex vivo validation were used to identify chromatin modifiers and readers for HCC chemoprevention. RESULTS: In patients, CHC and NASH share similar epigenetic and transcriptomic modifications driving cancer risk. Using a cell-based system modelling epigenetic modifications in patients, we identified chromatin readers as targets to revert liver gene transcription driving clinical HCC risk. Proof-of-concept studies in a NASH-HCC mouse model showed that the pharmacological inhibition of chromatin reader bromodomain 4 inhibited liver disease progression and hepatocarcinogenesis by restoring transcriptional reprogramming of the genes that were epigenetically altered in patients. CONCLUSION: Our results unravel the functional relevance of metabolic and virus-induced epigenetic alterations for pathogenesis of HCC development and identify chromatin readers as targets for chemoprevention in patients with chronic liver diseases.

Structural cells are key regulators of organ-specific immune responses.

The mammalian immune system implements a remarkably effective set of mechanisms for fighting pathogens1. Its main components are haematopoietic immune cells, including myeloid cells that control innate immunity, and lymphoid cells that constitute adaptive immunity2. However, immune functions are not unique to haematopoietic cells, and many other cell types display basic mechanisms of pathogen defence3-5. To advance our understanding of immunology outside the haematopoietic system, here we systematically investigate the regulation of immune genes in the three major types of structural cells: epithelium, endothelium and fibroblasts. We characterize these cell types across twelve organs in mice, using cellular phenotyping, transcriptome sequencing, chromatin accessibility profiling and epigenome mapping. This comprehensive dataset revealed complex immune gene activity and regulation in structural cells. The observed patterns were highly organ-specific and seem to modulate the extensive interactions between structural cells and haematopoietic immune cells. Moreover, we identified an epigenetically encoded immune potential in structural cells under tissue homeostasis, which was triggered in response to systemic viral infection. This study highlights the prevalence and organ-specific complexity of immune gene activity in non-haematopoietic structural cells, and it provides a high-resolution, multi-omics atlas of the epigenetic and transcriptional networks that regulate structural cells in the mouse.

Genome-wide cooperation of EMT transcription factor ZEB1 with YAP and AP-1 in breast cancer.

Invasion, metastasis and therapy resistance are the major cause of cancer-associated deaths, and the EMT-inducing transcription factor ZEB1 is a crucial stimulator of these processes. While work on ZEB1 has mainly focused on its role as a transcriptional repressor, it can also act as a transcriptional activator. To further understand these two modes of action, we performed a genome-wide ZEB1 binding study in triple-negative breast cancer cells. We identified ZEB1 as a novel interactor of the AP-1 factors FOSL1 and JUN and show that, together with the Hippo pathway effector YAP, they form a transactivation complex, predominantly activating tumour-promoting genes, thereby synergising with its function as a repressor of epithelial genes. High expression of ZEB1, YAP, FOSL1 and JUN marks the aggressive claudin-low subtype of breast cancer, indicating the translational relevance of our findings. Thus, our results link critical tumour-promoting transcription factors: ZEB1, AP-1 and Hippo pathway factors. Disturbing their molecular interaction may provide a promising treatment option for aggressive cancer types.

Targeting Nuclear NOTCH2 by Gliotoxin Recovers a Tumor-Suppressor NOTCH3 Activity in CLL.

NOTCH signaling represents a promising therapeutic target in chronic lymphocytic leukemia (CLL). We compared the anti-neoplastic effects of the nuclear NOTCH2 inhibitor gliotoxin and the pan-NOTCH γ-secretase inhibitor RO4929097 in primary CLL cells with special emphasis on the individual roles of the different NOTCH receptors. Gliotoxin rapidly induced apoptosis in all CLL cases tested, whereas RO4929097 exerted a variable and delayed effect on CLL cell viability. Gliotoxin-induced apoptosis was associated with inhibition of the NOTCH2/FCER2 (CD23) axis together with concomitant upregulation of the NOTCH3/NR4A1 axis. In contrast, RO4929097 downregulated the NOTCH3/NR4A1 axis and counteracted the spontaneous and gliotoxin-induced apoptosis. On the cell surface, NOTCH3 and CD23 expression were mutually exclusive, suggesting that downregulation of NOTCH2 signaling is a prerequisite for NOTCH3 expression in CLL cells. ATAC-seq confirmed that gliotoxin targeted the canonical NOTCH signaling, as indicated by the loss of chromatin accessibility at the potential NOTCH/CSL site containing the gene regulatory elements. This was accompanied by a gain in accessibility at the NR4A1, NFκB, and ATF3 motifs close to the genes involved in B-cell activation, differentiation, and apoptosis. In summary, these data show that gliotoxin recovers a non-canonical tumor-suppressing NOTCH3 activity, indicating that nuclear NOTCH2 inhibitors might be beneficial compared to pan-NOTCH inhibitors in the treatment of CLL.

Author Correction: Mechanisms governing the pioneering and redistribution capabilities of the non-classical pioneer PU.1.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Chromatin mapping and single-cell immune profiling define the temporal dynamics of ibrutinib response in CLL.

The Bruton tyrosine kinase (BTK) inhibitor ibrutinib provides effective treatment for patients with chronic lymphocytic leukemia (CLL), despite extensive heterogeneity in this disease. To define the underlining regulatory dynamics, we analyze high-resolution time courses of ibrutinib treatment in patients with CLL, combining immune-phenotyping, single-cell transcriptome profiling, and chromatin mapping. We identify a consistent regulatory program starting with a sharp decrease of NF-κB binding in CLL cells, which is followed by reduced activity of lineage-defining transcription factors, erosion of CLL cell identity, and acquisition of a quiescence-like gene signature. We observe patient-to-patient variation in the speed of execution of this program, which we exploit to predict patient-specific dynamics in the response to ibrutinib based on the pre-treatment patient samples. In aggregate, our study describes time-dependent cellular, molecular, and regulatory effects for therapeutic inhibition of B cell receptor signaling in CLL, and it establishes a broadly applicable method for epigenome/transcriptome-based treatment monitoring.

Mechanisms governing the pioneering and redistribution capabilities of the non-classical pioneer PU.1.

Establishing gene regulatory networks during differentiation or reprogramming requires master or pioneer transcription factors (TFs) such as PU.1, a prototype master TF of hematopoietic lineage differentiation. To systematically determine molecular features that control its activity, here we analyze DNA-binding in vitro and genome-wide in vivo across different cell types with native or ectopic PU.1 expression. Although PU.1, in contrast to classical pioneer factors, is unable to access nucleosomal target sites in vitro, ectopic induction of PU.1 leads to the extensive remodeling of chromatin and redistribution of partner TFs. De novo chromatin access, stable binding, and redistribution of partner TFs both require PU.1's N-terminal acidic activation domain and its ability to recruit SWI/SNF remodeling complexes, suggesting that the latter may collect and distribute co-associated TFs in conjunction with the non-classical pioneer TF PU.1.

Precursors for Nonlymphoid-Tissue Treg Cells Reside in Secondary Lymphoid Organs and Are Programmed by the Transcription Factor BATF.

Specialized regulatory T (Treg) cells accumulate and perform homeostatic and regenerative functions in nonlymphoid tissues. Whether common precursors for nonlymphoid-tissue Treg cells exist and how they differentiate remain elusive. Using transcription factor nuclear factor, interleukin 3 regulated (Nfil3) reporter mice and single-cell RNA-sequencing (scRNA-seq), we identified two precursor stages of interleukin 33 (IL-33) receptor ST2-expressing nonlymphoid tissue Treg cells, which resided in the spleen and lymph nodes. Global chromatin profiling of nonlymphoid tissue Treg cells and the two precursor stages revealed a stepwise acquisition of chromatin accessibility and reprogramming toward the nonlymphoid-tissue Treg cell phenotype. Mechanistically, we identified and validated the transcription factor Batf as the driver of the molecular tissue program in the precursors. Understanding this tissue development program will help to harness regenerative properties of tissue Treg cells for therapy.

Systematic characterization of BAF mutations provides insights into intracomplex synthetic lethalities in human cancers.

Aberrations in genes coding for subunits of the BRG1/BRM associated factor (BAF) chromatin remodeling complexes are highly abundant in human cancers. Currently, it is not understood how these mostly loss-of-function mutations contribute to cancer development and how they can be targeted therapeutically. The cancer-type-specific occurrence patterns of certain subunit mutations suggest subunit-specific effects on BAF complex function, possibly by the formation of aberrant residual complexes. Here, we systematically characterize the effects of individual subunit loss on complex composition, chromatin accessibility and gene expression in a panel of knockout cell lines deficient for 22 BAF subunits. We observe strong, specific and sometimes discordant alterations dependent on the targeted subunit and show that these explain intracomplex codependencies, including the synthetic lethal interactions SMARCA4-ARID2, SMARCA4-ACTB and SMARCC1-SMARCC2. These data provide insights into the role of different BAF subcomplexes in genome-wide chromatin organization and suggest approaches to therapeutically target BAF-mutant cancers.

MTHFD1 interaction with BRD4 links folate metabolism to transcriptional regulation.

The histone acetyl reader bromodomain-containing protein 4 (BRD4) is an important regulator of chromatin structure and transcription, yet factors modulating its activity have remained elusive. Here we describe two complementary screens for genetic and physical interactors of BRD4, which converge on the folate pathway enzyme MTHFD1 (methylenetetrahydrofolate dehydrogenase, cyclohydrolase and formyltetrahydrofolate synthetase 1). We show that a fraction of MTHFD1 resides in the nucleus, where it is recruited to distinct genomic loci by direct interaction with BRD4. Inhibition of either BRD4 or MTHFD1 results in similar changes in nuclear metabolite composition and gene expression; pharmacological inhibitors of the two pathways synergize to impair cancer cell viability in vitro and in vivo. Our finding that MTHFD1 and other metabolic enzymes are chromatin associated suggests a direct role for nuclear metabolism in the control of gene expression.

Rbpj expression in regulatory T cells is critical for restraining TH2 responses.

The transcriptional regulator Rbpj is involved in T-helper (TH) subset polarization, but its function in Treg cells remains unclear. Here we show that Treg-specific Rbpj deletion leads to splenomegaly and lymphadenopathy despite increased numbers of Treg cells with a polyclonal TCR repertoire. A specific defect of Rbpj-deficient Treg cells in controlling TH2 polarization and B cell responses is observed, leading to the spontaneous formation of germinal centers and a TH2-associated immunoglobulin class switch. The observed phenotype is environment-dependent and can be induced by infection with parasitic nematodes. Rbpj-deficient Treg cells adopt open chromatin landscapes and gene expression profiles reminiscent of tissue-derived TH2-polarized Treg cells, with a prevailing signature of the transcription factor Gata-3. Taken together, our study suggests that Treg cells require Rbpj to specifically restrain TH2 responses, including their own excessive TH2-like differentiation potential.

HCV-Induced Epigenetic Changes Associated With Liver Cancer Risk Persist After Sustained Virologic Response.

BACKGROUND & AIMS: Chronic hepatitis C virus (HCV) infection is an important risk factor for hepatocellular carcinoma (HCC). Despite effective antiviral therapies, the risk for HCC is decreased but not eliminated after a sustained virologic response (SVR) to direct-acting antiviral (DAA) agents, and the risk is higher in patients with advanced fibrosis. We investigated HCV-induced epigenetic alterations that might affect risk for HCC after DAA treatment in patients and mice with humanized livers. METHODS: We performed genome-wide ChIPmentation-based ChIP-Seq and RNA-seq analyses of liver tissues from 6 patients without HCV infection (controls), 18 patients with chronic HCV infection, 8 patients with chronic HCV infection cured by DAA treatment, 13 patients with chronic HCV infection cured by interferon therapy, 4 patients with chronic hepatitis B virus infection, and 7 patients with nonalcoholic steatohepatitis in Europe and Japan. HCV-induced epigenetic modifications were mapped by comparative analyses with modifications associated with other liver disease etiologies. uPA/SCID mice were engrafted with human hepatocytes to create mice with humanized livers and given injections of HCV-infected serum samples from patients; mice were given DAAs to eradicate the virus. Pathways associated with HCC risk were identified by integrative pathway analyses and validated in analyses of paired HCC tissues from 8 patients with an SVR to DAA treatment of HCV infection. RESULTS: We found chronic HCV infection to induce specific genome-wide changes in H3K27ac, which correlated with changes in expression of mRNAs and proteins. These changes persisted after an SVR to DAAs or interferon-based therapies. Integrative pathway analyses of liver tissues from patients and mice with humanized livers demonstrated that HCV-induced epigenetic alterations were associated with liver cancer risk. Computational analyses associated increased expression of SPHK1 with HCC risk. We validated these findings in an independent cohort of patients with HCV-related cirrhosis (n = 216), a subset of which (n = 21) achieved viral clearance. CONCLUSIONS: In an analysis of liver tissues from patients with and without an SVR to DAA therapy, we identified epigenetic and gene expression alterations associated with risk for HCC. These alterations might be targeted to prevent liver cancer in patients treated for HCV infection.