Ets transcription factor GABP controls T cell homeostasis and immunity

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Histone Lactylation Drives CD8 T Cell Metabolism and Function

The activation and functional differentiation of CD8 T cells are linked to metabolic pathways that result in the production of lactate. Lactylation is a lactate-derived histone post-translational modification (hPTM); however, the relevance of histone lactylation in the context of CD8 T cell activation and function is not known. Here, we show the enrichment of H3K18-lactylation (H3K18la) and H3K9-lactylation (H3K9la) in human and murine CD8 T cells which act as transcription initiators of key genes regulating CD8 T cell phenotype and function. Further, we note distinct impacts of H3K18la and H3K9la on CD8 T cell subsets linked to their specific metabolic profiles. Importantly, we demonstrate that modulation of H3K18la and H3K9la by targeting metabolic and epigenetic pathways regulates CD8 T cell effector function including anti-tumor immunity in preclinical models. Overall, our study uncovers the unique contributions of H3K18la and H3K9la in modulating CD8 T cell phenotype and function intricately associated with metabolic state.

Gfi1 controls the formation of effector CD8 T cells during chronic infection and cancer

During chronic infections and tumor progression, CD8 T cells gradually lose their effector functions and become exhausted. These exhausted CD8 T cells are heterogeneous and comprised of different subsets, including self-renewing progenitors that give rise to Ly108 - CX3CR1 + effector-like cells. Generation of these effector-like cells is essential for the control of chronic infections and tumors, albeit limited. However, the precise cues and mechanisms directing the formation and maintenance of exhausted effector-like are incompletely understood. Using genetic mouse models challenged with LCMV Clone 13 or syngeneic tumors, we show that the expression of a transcriptional repressor, growth factor independent 1 (Gfi1) is dynamically regulated in exhausted CD8 T cells, which in turn regulates the formation of exhausted effector-like cells. Gfi1 deletion in T cells dysregulates the chromatin accessibility and transcriptomic programs associated with the differentiation of LCMV Clone 13-specific CD8 T cell exhaustion, preventing the formation of effector-like and terminally exhausted cells while maintaining progenitors and a newly identified Ly108 + CX3CR1 + state. These Ly108 + CX3CR1 + cells have a distinct chromatin profile and may represent an alternative target for therapeutic interventions to combat chronic infections and cancer. In sum, we show that Gfi1 is a critical regulator of the formation of exhausted effector-like cells.

Functional and transcriptional heterogeneity within the massively expanding HLADR+CD38+ CD8 T cell population in acute febrile dengue patients

IMPORTANCE

CD8 T cells play a crucial role in protecting against intracellular pathogens such as viruses by eliminating infected cells and releasing anti-viral cytokines such as interferon gamma (IFNγ). Consequently, there is significant interest in comprehensively characterizing CD8 T cell responses in acute dengue febrile patients. Previous studies, including our own, have demonstrated that a discrete population of CD8 T cells with HLADR+ CD38+ phenotype undergoes massive expansion during the acute febrile phase of natural dengue virus infection. Although about a third of these massively expanding HLADR+ CD38+ CD8 T cells were also CD69high when examined ex vivo, only a small fraction of them produced IFNγ upon in vitro peptide stimulation. Therefore, to better understand such functional diversity of CD8 T cells responding to dengue virus infection, it is important to know the cytokines/chemokines expressed by these peptide-stimulated HLADR+CD38+ CD8 T cells and the transcriptional profiles that distinguish the CD69+IFNγ+, CD69+IFNγ-, and CD69-IFNγ- subsets.

NRF1-mediated mitochondrial biogenesis antagonizes innate antiviral immunity

Mitochondrial biogenesis is the process of generating new mitochondria to maintain cellular homeostasis. Here, we report that viruses exploit mitochondrial biogenesis to antagonize innate antiviral immunity. We found that nuclear respiratory factor-1 (NRF1), a vital transcriptional factor involved in nuclear-mitochondrial interactions, is essential for RNA (VSV) or DNA (HSV-1) virus-induced mitochondrial biogenesis. NRF1 deficiency resulted in enhanced innate immunity, a diminished viral load, and morbidity in mice. Mechanistically, the inhibition of NRF1-mediated mitochondrial biogenesis aggravated virus-induced mitochondrial damage, promoted the release of mitochondrial DNA (mtDNA), increased the production of mitochondrial reactive oxygen species (mtROS), and activated the innate immune response. Notably, virus-activated kinase TBK1 phosphorylated NRF1 at Ser318 and thereby triggered the inactivation of the NRF1-TFAM axis during HSV-1 infection. A knock-in (KI) strategy that mimicked TBK1-NRF1 signaling revealed that interrupting the TBK1-NRF1 connection ablated mtDNA release and thereby attenuated the HSV-1-induced innate antiviral response. Our study reveals a previously unidentified antiviral mechanism that utilizes a NRF1-mediated negative feedback loop to modulate mitochondrial biogenesis and antagonize innate immune response.

Tox4 regulates transcriptional elongation and reinitiation during murine T cell development

HMG protein Tox4 is a regulator of PP1 phosphatases with unknown function in development. Here we show that Tox4 conditional knockout in mice reduces thymic cellularity, partially blocks T cell development, and decreases ratio of CD8 to CD4 through decreasing proliferation and increasing apoptosis of CD8 cells. In addition, single-cell RNA-seq discovered that Tox4 loss also impairs proliferation of the fast-proliferating double positive (DP) blast population within DP cells in part due to downregulation of genes critical for proliferation, notably Cdk1. Moreover, genes with high and low expression level are more dependent on Tox4 than genes with medium expression level. Mechanistically, Tox4 may facilitate transcriptional reinitiation and restrict elongation in a dephosphorylation-dependent manner, a mechanism that is conserved between mouse and human. These results provide insights into the role of TOX4 in development and establish it as an evolutionarily conserved regulator of transcriptional elongation and reinitiation.

Multi-omic approach associates blood methylome with bronchodilator drug response in pediatric asthma

BACKGROUND

Albuterol is the drug most widely used as asthma treatment among African Americans despite having a lower bronchodilator drug response (BDR) than other populations. Although BDR is affected by gene and environmental factors, the influence of DNA methylation is unknown.

OBJECTIVE

This study aimed to identify epigenetic markers in whole blood associated with BDR, study their functional consequences by multi-omic integration, and assess their clinical applicability in admixed populations with a high asthma burden.

METHODS

We studied 414 children and young adults (8-21 years old) with asthma in a discovery and replication design. We performed an epigenome-wide association study on 221 African Americans and replicated the results on 193 Latinos. Functional consequences were assessed by integrating epigenomics with genomics, transcriptomics, and environmental exposure data. Machine learning was used to develop a panel of epigenetic markers to classify treatment response.

RESULTS

We identified 5 differentially methylated regions and 2 CpGs genome-wide significantly associated with BDR in African Americans located in FGL2 (cg08241295, P = 6.8 × 10-9) and DNASE2 (cg15341340, P = 7.8 × 10-8), which were regulated by genetic variation and/or associated with gene expression of nearby genes (false discovery rate < 0.05). The CpG cg15341340 was replicated in Latinos (P = 3.5 × 10-3). Moreover, a panel of 70 CpGs showed good classification for those with response and nonresponse to albuterol therapy in African American and Latino children (area under the receiver operating characteristic curve for training, 0.99; for validation, 0.70-0.71). The DNA methylation model showed similar discrimination as clinical predictors (P > .05).

CONCLUSIONS

We report novel associations of epigenetic markers with BDR in pediatric asthma and demonstrate for the first time the applicability of pharmacoepigenetics in precision medicine of respiratory diseases.

NetBID2 provides comprehensive hidden driver analysis

Many signaling and other genes known as "hidden" drivers may not be genetically or epigenetically altered or differentially expressed at the mRNA or protein levels, but, rather, drive a phenotype such as tumorigenesis via post-translational modification or other mechanisms. However, conventional approaches based on genomics or differential expression are limited in exposing such hidden drivers. Here, we present a comprehensive algorithm and toolkit NetBID2 (data-driven network-based Bayesian inference of drivers, version 2), which reverse-engineers context-specific interactomes and integrates network activity inferred from large-scale multi-omics data, empowering the identification of hidden drivers that could not be detected by traditional analyses. NetBID2 has substantially re-engineered the previous prototype version by providing versatile data visualization and sophisticated statistical analyses, which strongly facilitate researchers for result interpretation through end-to-end multi-omics data analysis. We demonstrate the power of NetBID2 using three hidden driver examples. We deploy NetBID2 Viewer, Runner, and Cloud apps with 145 context-specific gene regulatory and signaling networks across normal tissues and paediatric and adult cancers to facilitate end-to-end analysis, real-time interactive visualization and cloud-based data sharing. NetBID2 is freely available at https://jyyulab.github.io/NetBID .

NRF-2α and mitophagy underlie enhanced mitochondrial functions and biogenesis induced by T-2 toxin in GH3 cells

T-2 toxin is a natural contaminant in grain cereals produced by species of Fusarium. Studies indicate that T-2 toxin can positively affect mitochondrial function, but the underlying mechanism is unclear. In this study, we examined the role of nuclear respiratory factor 2α (NRF-2α) in T-2 toxin-activated mitochondrial biogenesis and the direct target genes of NRF-2α. Furthermore, we investigated T-2 toxin-induced autophagy and mitophagy, and the role of mitophagy in changes in mitochondrial function and apoptosis. It was found that T-2 toxin significantly increased NRF-2α levels and nuclear localization of NRF-2α was induced. NRF-2α deletion significantly increased the production of reactive oxygen species (ROS), abrogated T-2 toxin-induced increases in ATP and mitochondrial complex I activity, and inhibited the mitochondrial DNA copy number. Meanwhile, With chromatin immunoprecipitation sequencing (ChIP-Seq), various novel NRF-2α target genes were identified, such as mitochondrial iron-sulphur subunits (Ndufs 3,7) and mitochondrial transcription factors (Tfam, Tfb1m, and Tfb2m). Some target genes were also involved in mitochondrial fusion and fission (Drp1), mitochondrial translation (Yars2) and splicing (Ddx55), and mitophagy. Further studies showed that T-2 toxin induced Atg5 dependent autophagy and Atg5/PINK1-dependent mitophagy. In addition, mitophagy defects increase ROS production, inhibit ATP levels and the expression of genes related to mitochondrial dynamics, and promote apoptosis in the presence of T-2 toxins. Altogether, these results suggest that NRF-2α plays a critical role in promoting mitochondrial function and biogenesis through regulation of mitochondrial genes, and, interestingly, mitophagy caused by T-2 toxin positively affected mitochondrial function and protected cell survival against T-2 toxin.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Cluster-independent marker feature identification from single-cell omics data using SEMITONES

Identification of cell identity markers is an essential step in single-cell omics data analysis. Current marker identification strategies typically rely on cluster assignments of cells. However, cluster assignment, particularly for developmental data, is nontrivial, potentially arbitrary, and commonly relies on prior knowledge. In response, we present SEMITONES, a principled method for cluster-free marker identification. We showcase and evaluate its application for marker gene and regulatory region identification from single-cell data of the human haematopoietic system. Additionally, we illustrate its application to spatial transcriptomics data and show how SEMITONES can be used for the annotation of cells given known marker genes. Using several simulated and curated data sets, we demonstrate that SEMITONES qualitatively and quantitatively outperforms existing methods for the retrieval of cell identity markers from single-cell omics data.

Hierarchical regulation of the resting and activated T cell epigenome by major transcription factor families

T cell activation, a key early event in the adaptive immune response, is subject to elaborate transcriptional control. Here, we examined how the activities of eight major transcription factor (TF) families are integrated to shape the epigenome of naïve and activated CD4 and CD8 T cells. By leveraging extensive polymorphisms in evolutionarily divergent mice, we identified the “heavy lifters” positively influencing chromatin accessibility. Members of Ets, Runx, and TCF/Lef TF families occupied the vast majority of accessible chromatin regions, acting as “housekeepers”, “universal amplifiers”, and “placeholders”, respectively, at sites that maintained or gained accessibility upon T cell activation. Additionally, a small subset of strongly induced immune response genes displayed a non-canonical TF recruitment pattern. Our study provides a key resource and foundation for the understanding of transcriptional and epigenetic regulation in T cells and offers a new perspective on the hierarchical interactions between critical TFs.

Whole-genome analysis of TET dioxygenase function in regulatory T cells

TET methylcytosine dioxygenases are essential for the stability and function of regulatory T cells (Treg cells), which maintain immune homeostasis and self‐tolerance and express the lineage‐determining transcription factor Foxp3. Here, we use whole‐genome analyses to show that the transcriptional program and epigenetic features (DNA modification, chromatin accessibility) of Treg cells are attenuated in the absence of Tet2 and Tet3. Conversely, the addition of the TET activator vitamin C during TGFβ‐induced iTreg cell differentiation in vitro potentiates the expression of Treg signature genes and alters the epigenetic landscape to better resemble that of Treg cells generated in vivo. Vitamin C enhances IL‐2 responsiveness in iTreg cells by increasing IL2Rα expression, STAT5 phosphorylation, and STAT5 binding, mimicking the IL‐2/STAT5 dependence of Treg cells generated in vivo. In summary, TET proteins play essential roles in maintaining Treg molecular features and promoting their dependence on IL‐2. TET activity during endogenous Treg development and potentiation of TET activity by vitamin C during iTreg differentiation are necessary to maintain the transcriptional and epigenetic features of Treg cells.

Resident memory T cells in tumor-distant tissues fortify against metastasis formation

As a critical machinery for rapid pathogen removal, resident memory T cells (T_RMs) are locally generated after the initial encounter. However, their development accompanying tumorigenesis remains elusive. Using a murine breast cancer model, we show that T_RMs develop in the tumor, the contralateral mammary mucosa, and the pre-metastatic lung. Single-cell RNA sequencing of T_RMs reveals two phenotypically distinct populations representing their active versus quiescent phases. These T_RMs in different tissue compartments share the same TCR clonotypes and transcriptomes with a subset of intratumoral effector/effector memory T cells (T_Eff/EMs), indicating their developmental ontogeny. Furthermore, CXCL16 is highly produced by tumor cells and CXCR6^- T_Eff/EMs are the major subset preferentially egressing the tumor to form distant T_RMs. Functionally, releasing CXCR6 retention in the primary tumor amplifies tumor-derived T_RMs in the lung and leads to superior protection against metastases. This immunologic fortification suggests a potential strategy to prevent metastasis in clinical oncology.

GABPA is a master regulator of luminal identity and restrains aggressive diseases in bladder cancer

TERT promoter mutations occur in the majority of glioblastoma, bladder cancer (BC), and other malignancies while the ETS family transcription factors GABPA and its partner GABPB1 activate the mutant TERT promoter and telomerase in these tumors. GABPA depletion or the disruption of the GABPA/GABPB1 complex by knocking down GABPB1 was shown to inhibit telomerase, thereby eliminating the tumorigenic potential of glioblastoma cells. GABPA/B1 is thus suggested as a cancer therapeutic target. However, it is unclear about its role in BC. Here we unexpectedly observed that GABPA ablation inhibited TERT expression, but robustly increased proliferation, stem, and invasive phenotypes and cisplatin resistance in BC cells, while its overexpression exhibited opposite effects, and inhibited in vivo metastasizing in a xenograft transplant model. Mechanistically, GABPA directly activates the transcription of FoxA1 and GATA3, key transcription factors driving luminal differentiation of urothelial cells. Consistently, TCGA/GEO dataset analyses show that GABPA expression is correlated positively with luminal while negatively with basal signatures. Luminal tumors express higher GABPA than do basal ones. Lower GABPA expression is associated with the GABPA gene methylation or deletion (especially in basal subtype of BC tumors), and predicted significantly shorter patient survival based on TCGA and our cohort of BC patient analyses. Taken together, GABPA dictates luminal identity of BC cells and inhibits aggressive diseases in BC by promoting cellular differentiation despite its stimulatory effect on telomerase/TERT activation. Given these biological functions and its frequent methylation and/or deletion, GABPA serves as a tumor suppressor rather than oncogenic factor in BC. The GABPA effect on oncogenesis is context-dependent and its targeting for telomerase inhibition in BC may promote disease metastasizing.

Specific subfamilies of transposable elements contribute to different domains of T lymphocyte enhancers

Transposable elements (TEs) compose nearly half of mammalian genomes and provide building blocks for cis-regulatory elements. Using high-throughput sequencing, we show that 84 TE subfamilies are overrepresented, and distributed in a lineage-specific fashion in core and boundary domains of CD8⁺ T cell enhancers. Endogenous retroviruses are most significantly enriched in core domains with accessible chromatin, and bear recognition motifs for immune-related transcription factors. In contrast, short interspersed elements (SINEs) are preferentially overrepresented in nucleosome-containing boundaries. A substantial proportion of these SINEs harbor a high density of the enhancer-specific histone mark H3K4me1 and carry sequences that match enhancer boundary nucleotide composition. Motifs with regulatory features are better preserved within enhancer-enriched TE copies compared to their subfamily equivalents located in gene deserts. TE-rich and TE-poor enhancers associate with both shared and unique gene groups and are enriched in overlapping functions related to lymphocyte and leukocyte biology. The majority of T cell enhancers are shared with other immune lineages and are accessible in common hematopoietic progenitors. A higher proportion of immune tissue-specific enhancers are TE-rich compared to enhancers specific to other tissues, correlating with higher TE occurrence in immune gene-associated genomic regions. Our results suggest that during evolution, TEs abundant in these regions and carrying motifs potentially beneficial for enhancer architecture and immune functions were particularly frequently incorporated by evolving enhancers. Their putative selection and regulatory cooption may have accelerated the evolution of immune regulatory networks.

Survival of Naïve T Cells Requires the Expression of Let-7 miRNAs

Maintaining the diversity and constant numbers of naïve T cells throughout the organism's lifetime is necessary for efficient immune responses. Naïve T cell homeostasis, which consists of prolonged survival, occasional proliferation and enforcement of quiescence, is tightly regulated by multiple signaling pathways which are in turn controlled by various transcription factors. However, full understanding of the molecular mechanisms underlying the maintenance of the peripheral T cell pool has not been achieved. In the present study, we demonstrate that T cell-specific deficiency in let-7 miRNAs results in peripheral T cell lymphopenia resembling that of Dicer1 knockout mice. Deletion of let-7 leads to profound T cell apoptosis while overexpression prevents it. We further show that in the absence of let-7, T cells cannot sustain optimal levels of the pro-survival factor Bcl2 in spite of the intact IL-7 signaling, and re-expression of Bcl2 in let-7 deficient T cells completely rescues the survival defect. Thus, we have uncovered a novel let-7-dependent mechanism of post-transcriptional regulation of naïve T cell survival in vivo.