AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs for atlas-scale datasets like Human Pancreas Analysis Program (HPAP), we develop AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX shows the higher performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulates known islet pathobiology and shows differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

Computational workflow and interactive analysis of single-cell expression profiling of islets generated by the Human Pancreas Analysis Program.

Type 1 and Type 2 diabetes are distinct genetic diseases of the pancreas which are defined by the abnormal level of blood glucose. Understanding the initial molecular perturbations that occur during the pathogenesis of diabetes is of critical importance in understanding these disorders. The inability to biopsy the human pancreas of living donors hampers insights into early detection, as the majority of diabetes studies have been performed on peripheral leukocytes from the blood, which is not the site of pathogenesis. Therefore, efforts have been made by various teams including the Human Pancreas Analysis Program (HPAP) to collect pancreatic tissues from deceased organ donors with different clinical phenotypes. HPAP is designed to define the molecular pathogenesis of islet dysfunction by generating detailed datasets of functional, cellular, and molecular information in pancreatic tissues of clinically well-defined organ donors with Type 1 and Type 2 diabetes. Moreover, data generated by HPAP continously become available through a centralized database, PANC-DB, thus enabling the diabetes research community to access these multi-dimensional data prepublication. Here, we present the computational workflow for single-cell RNA-seq data analysis of 258,379 high-quality cells from the pancreatic islets of 67 human donors generated by HPAP, the largest existing scRNA-seq dataset of human pancreatic tissues. We report various computational steps including preprocessing, doublet removal, clustering and cell type annotation across single-cell RNA-seq data from islets of four distintct classes of organ donors, i.e. non-diabetic control, autoantibody positive but normoglycemic, Type 1 diabetic, and Type 2 diabetic individuals. Moreover, we present an interactive tool, called CellxGene developed by the Chan Zuckerberg initiative, to navigate these high-dimensional datasets. Our data and interactive tools provide a reliable reference for singlecell pancreatic islet biology studies, especially diabetes-related conditions.

AnnoSpat annotates cell types and quantifies cellular arrangements from spatial proteomics.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs, we developed AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX show the superior performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulated known islet pathobiology and showed differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8+ T cells infiltration in islets during type 1 diabetes progression.

BET bromodomain protein inhibition reverses chimeric antigen receptor extinction and reinvigorates exhausted T cells in chronic lymphocytic leukemia.

Chimeric antigen receptor (CAR) T cells have induced remarkable antitumor responses in B cell malignancies. Some patients do not respond because of T cell deficiencies that hamper the expansion, persistence, and effector function of these cells. We used longitudinal immune profiling to identify phenotypic and pharmacodynamic changes in CD19-directed CAR T cells in patients with chronic lymphocytic leukemia (CLL). CAR expression maintenance was also investigated because this can affect response durability. CAR T cell failure was accompanied by preexisting T cell-intrinsic defects or dysfunction acquired after infusion. In a small subset of patients, CAR silencing was observed coincident with leukemia relapse. Using a small molecule inhibitor, we demonstrated that the bromodomain and extra-terminal (BET) family of chromatin adapters plays a role in downregulating CAR expression. BET protein blockade also ameliorated CAR T cell exhaustion as manifested by inhibitory receptor reduction, enhanced metabolic fitness, increased proliferative capacity, and enriched transcriptomic signatures of T cell reinvigoration. BET inhibition decreased levels of the TET2 methylcytosine dioxygenase, and forced expression of the TET2 catalytic domain eliminated the potency-enhancing effects of BET protein targeting in CAR T cells, providing a mechanism linking BET proteins and T cell dysfunction. Thus, modulating BET epigenetic readers may improve the efficacy of cell-based immunotherapies.

Epigenetic scarring of exhausted T cells hinders memory differentiation upon eliminating chronic antigenic stimulation.

Exhausted CD8 T cells (TEX) are a distinct state of T cell differentiation associated with failure to clear chronic viruses and cancer. Immunotherapies such as PD-1 blockade can reinvigorate TEX cells, but reinvigoration is not durable. A major unanswered question is whether TEX cells differentiate into functional durable memory T cells (TMEM) upon antigen clearance. Here, using a mouse model, we found that upon eliminating chronic antigenic stimulation, TEX cells partially (re)acquire phenotypic and transcriptional features of TMEM cells. These 'recovering' TEX cells originated from the T cell factor (TCF-1+) TEX progenitor subset. Nevertheless, the recall capacity of these recovering TEX cells remained compromised as compared to TMEM cells. Chromatin-accessibility profiling revealed a failure to recover core memory epigenetic circuits and maintenance of a largely exhausted open chromatin landscape. Thus, despite some phenotypic and transcriptional recovery upon antigen clearance, exhaustion leaves durable epigenetic scars constraining future immune responses. These results support epigenetic remodeling interventions for TEX cell-targeted immunotherapies.

The Identity of Human Tissue-Emigrant CD8+ T Cells.

Lymphocyte migration is essential for adaptive immune surveillance. However, our current understanding of this process is rudimentary, because most human studies have been restricted to immunological analyses of blood and various tissues. To address this knowledge gap, we used an integrated approach to characterize tissue-emigrant lineages in thoracic duct lymph (TDL). The most prevalent immune cells in human and non-human primate efferent lymph were T cells. Cytolytic CD8+ T cell subsets with effector-like epigenetic and transcriptional signatures were clonotypically skewed and selectively confined to the intravascular circulation, whereas non-cytolytic CD8+ T cell subsets with stem-like epigenetic and transcriptional signatures predominated in tissues and TDL. Moreover, these anatomically distinct gene expression profiles were recapitulated within individual clonotypes, suggesting parallel differentiation programs independent of the expressed antigen receptor. Our collective dataset provides an atlas of the migratory immune system and defines the nature of tissue-emigrant CD8+ T cells that recirculate via TDL.

Joint profiling of chromatin accessibility and CAR-T integration site analysis at population and single-cell levels.

Chimeric antigen receptor (CAR)-T immunotherapy has yielded impressive results in several B cell malignancies, establishing itself as a powerful means to redirect the natural properties of T lymphocytes. In this strategy, the T cell genome is modified by the integration of lentiviral vectors encoding CAR that direct tumor cell killing. However, this therapeutic approach is often limited by the extent of CAR-T cell expansion in vivo. A major outstanding question is whether or not CAR-T integration itself enhances the proliferative competence of individual T cells by rewiring their regulatory landscape. To address this question, it is critical to define the identity of an individual CAR-T cell and simultaneously chart where the CAR-T vector integrates into the genome. Here, we report the development of a method called EpiVIA (https://github.com/VahediLab/epiVIA) for the joint profiling of the chromatin accessibility and lentiviral integration site analysis at the population and single-cell levels. We validate our technique in clonal cells with previously defined integration sites and further demonstrate the ability to measure lentiviral integration sites and chromatin accessibility of host and viral genomes at the single-cell resolution in CAR-T cells. We anticipate that EpiVIA will enable the single-cell deconstruction of gene regulation during CAR-T therapy, leading to the discovery of cellular factors associated with durable treatment.

TCF-1-Centered Transcriptional Network Drives an Effector versus Exhausted CD8 T Cell-Fate Decision.

TCF-1 is a key transcription factor in progenitor exhausted CD8 T cells (Tex). Moreover, this Tex cell subset mediates responses to PD-1 checkpoint pathway blockade. However, the role of the transcription factor TCF-1 in early fate decisions and initial generation of Tex cells is unclear. Single-cell RNA sequencing (scRNA-seq) and lineage tracing identified a TCF-1+Ly108+PD-1+ CD8 T cell population that seeds development of mature Tex cells early during chronic infection. TCF-1 mediated the bifurcation between divergent fates, repressing development of terminal KLRG1Hi effectors while fostering KLRG1Lo Tex precursor cells, and PD-1 stabilized this TCF-1+ Tex precursor cell pool. TCF-1 mediated a T-bet-to-Eomes transcription factor transition in Tex precursors by promoting Eomes expression and drove c-Myb expression that controlled Bcl-2 and survival. These data define a role for TCF-1 in early-fate-bifurcation-driving Tex precursor cells and also identify PD-1 as a protector of this early TCF-1 subset.

A mass spectrometry-based assay using metabolic labeling to rapidly monitor chromatin accessibility of modified histone proteins.

Histone post-translational modifications (PTMs) contribute to chromatin accessibility due to their chemical properties and their ability to recruit enzymes responsible for DNA readout and chromatin remodeling. To date, more than 400 different histone PTMs and thousands of combinations of PTMs have been identified, the vast majority with still unknown biological function. Identification and quantification of histone PTMs has become routine in mass spectrometry (MS) but, since raising antibodies for each PTM in a study can be prohibitive, lots of potential is lost from MS datasets when uncharacterized PTMs are found to be significantly regulated. We developed an assay that uses metabolic labeling and MS to associate chromatin accessibility with histone PTMs and their combinations. The labeling is achieved by spiking in the cell media a 5x concentration of stable isotope labeled arginine and allow cells to grow for at least one cell cycle. We quantified the labeling incorporation of about 200 histone peptides with a proteomics workflow, and we confirmed that peptides carrying PTMs with extensively characterized roles in active transcription or gene silencing were in highly or poorly labeled forms, respectively. Data were further validated using next-generation sequencing to assess the transcription rate of chromatin regions modified with five selected PTMs. Furthermore, we quantified the labeling rate of peptides carrying co-existing PTMs, proving that this method is suitable for combinatorial PTMs. We focus on the abundant bivalent mark H3K27me3K36me2, showing that H3K27me3 dominantly represses histone swapping rate even in the presence of the more permissive PTM H3K36me2. Together, we envision this method will help to generate hypotheses regarding histone PTM functions and, potentially, elucidate the role of combinatorial histone codes.

A Thpok-Directed Transcriptional Circuitry Promotes Bcl6 and Maf Expression to Orchestrate T Follicular Helper Differentiation.

The generation of high-affinity neutralizing antibodies, the objective of most vaccine strategies, occurs in B cells within germinal centers (GCs) and requires rate-limiting "help" from follicular helper CD4+ T (Tfh) cells. Although Tfh differentiation is an attribute of MHC II-restricted CD4+ T cells, the transcription factors driving Tfh differentiation, notably Bcl6, are not restricted to CD4+ T cells. Here, we identified a requirement for the CD4+-specific transcription factor Thpok during Tfh cell differentiation, GC formation, and antibody maturation. Thpok promoted Bcl6 expression and bound to a Thpok-responsive region in the first intron of Bcl6. Thpok also promoted the expression of Bcl6-independent genes, including the transcription factor Maf, which cooperated with Bcl6 to mediate the effect of Thpok on Tfh cell differentiation. Our findings identify a transcriptional program that links the CD4+ lineage with Tfh differentiation, a limiting factor for efficient B cell responses, and suggest avenues to optimize vaccine generation.

Oncogenic Notch Promotes Long-Range Regulatory Interactions within Hyperconnected 3D Cliques.

Chromatin loops enable transcription-factor-bound distal enhancers to interact with their target promoters to regulate transcriptional programs. Although developmental transcription factors such as active forms of Notch can directly stimulate transcription by activating enhancers, the effect of their oncogenic subversion on the 3D organization of cancer genomes is largely undetermined. By mapping chromatin looping genome-wide in Notch-dependent triple-negative breast cancer and B cell lymphoma, we show that beyond the well-characterized role of Notch as an activator of distal enhancers, Notch regulates its direct target genes by instructing enhancer repositioning. Moreover, a large fraction of Notch-instructed regulatory loops form highly interacting enhancer and promoter spatial clusters termed "3D cliques." Loss- and gain-of-function experiments show that Notch preferentially targets hyperconnected 3D cliques that regulate the expression of crucial proto-oncogenes. Our observations suggest that oncogenic hijacking of developmental transcription factors can dysregulate transcription through widespread effects on the spatial organization of cancer genomes.

Epigenetic stability of exhausted T cells limits durability of reinvigoration by PD-1 blockade.

Blocking Programmed Death-1 (PD-1) can reinvigorate exhausted CD8 T cells (TEX) and improve control of chronic infections and cancer. However, whether blocking PD-1 can reprogram TEX into durable memory T cells (TMEM) is unclear. We found that reinvigoration of TEX in mice by PD-L1 blockade caused minimal memory development. After blockade, reinvigorated TEX became reexhausted if antigen concentration remained high and failed to become TMEM upon antigen clearance. TEX acquired an epigenetic profile distinct from that of effector T cells (TEFF) and TMEM cells that was minimally remodeled after PD-L1 blockade. This finding suggests that TEX are a distinct lineage of CD8 T cells. Nevertheless, PD-1 pathway blockade resulted in transcriptional rewiring and reengagement of effector circuitry in the TEX epigenetic landscape. These data indicate that epigenetic fate inflexibility may limit current immunotherapies.

Targeted genomic analysis reveals widespread autoimmune disease association with regulatory variants in the TNF superfamily cytokine signalling network.

BACKGROUND: Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored. METHODS: We took a hypothesis-driven approach using available genome-wide datasets to identify genetic variants regulating gene expression in the TNF superfamily cytokine signalling network and the association of these variants with autoimmune and autoinflammatory disease. Using paired gene expression and genetic data, we identified genetic variants associated with gene expression, expression quantitative trait loci (eQTLs), in four peripheral blood cell subsets. We then examined whether eQTLs were dependent on gene expression level or the presence of active enhancer chromatin marks. Using these eQTLs as genetic markers of the TNF superfamily signalling network, we performed targeted gene set association analysis in eight autoimmune and autoinflammatory disease genome-wide association studies. RESULTS: Comparison of TNF superfamily network gene expression and regulatory variants across four leucocyte subsets revealed patterns that differed between cell types. eQTLs for genes in this network were not dependent on absolute gene expression levels and were not enriched for chromatin marks of active enhancers. By examining autoimmune disease risk variants among our eQTLs, we found that risk alleles can be associated with either increased or decreased expression of co-stimulatory TNF superfamily cytokines, receptors or downstream signalling molecules. Gene set disease association analysis revealed that eQTLs for genes in the TNF superfamily pathway were associated with six of the eight autoimmune and autoinflammatory diseases examined, demonstrating associations beyond single genome-wide significant hits. CONCLUSIONS: This systematic analysis of the influence of regulatory genetic variants in the TNF superfamily network reveals widespread and diverse roles for these cytokines in susceptibility to a number of immune-mediated diseases.

BRD4 assists elongation of both coding and enhancer RNAs by interacting with acetylated histones.

Small-molecule BET inhibitors interfere with the epigenetic interactions between acetylated histones and the bromodomains of the BET family proteins, including BRD4, and they potently inhibit growth of malignant cells by targeting cancer-promoting genes. BRD4 interacts with the pause-release factor P-TEFb and has been proposed to release RNA polymerase II (Pol II) from promoter-proximal pausing. We show that BRD4 occupies widespread genomic regions in mouse cells and directly stimulates elongation of both protein-coding transcripts and noncoding enhancer RNAs (eRNAs), in a manner dependent on bromodomain function. BRD4 interacts with elongating Pol II complexes and assists Pol II in progression through hyperacetylated nucleosomes by interacting with acetylated histones via bromodomains. On active enhancers, the BET inhibitor JQ1 antagonizes BRD4-associated eRNA synthesis. Thus, BRD4 is involved in multiple steps of the transcription hierarchy, primarily by facilitating transcript elongation both at enhancers and on gene bodies independently of P-TEFb.

A mouse model of HIES reveals pro- and anti-inflammatory functions of STAT3.

Mutations of STAT3 underlie the autosomal dominant form of hyperimmunoglobulin E syndrome (HIES). STAT3 has critical roles in immune cells and thus, hematopoietic stem cell transplantation (HSCT), might be a reasonable therapeutic strategy in this disease. However, STAT3 also has critical functions in nonhematopoietic cells and dissecting the protean roles of STAT3 is limited by the lethality associated with germline deletion of Stat3. Thus, predicting the efficacy of HSCT for HIES is difficult. To begin to dissect the importance of STAT3 in hematopoietic and nonhematopoietic cells as it relates to HIES, we generated a mouse model of this disease. We found that these transgenic mice recapitulate multiple aspects of HIES, including elevated serum IgE and failure to generate Th17 cells. We found that these mice were susceptible to bacterial infection that was partially corrected by HSCT using wild-type bone marrow, emphasizing the role played by the epithelium in the pathophysiology of HIES.

Type I IFN induces binding of STAT1 to Bcl6: divergent roles of STAT family transcription factors in the T follicular helper cell genetic program.

CD4(+) T follicular helper cells (TFH) are critical for the formation and function of B cell responses to infection or immunization, but also play an important role in autoimmunity. The factors that contribute to the differentiation of this helper cell subset are incompletely understood, although several cytokines including IL-6, IL-21, and IL-12 can promote TFH cell formation. Yet, none of these factors, nor their downstream cognate STATs, have emerged as nonredundant, essential drivers of TFH cells. This suggests a model in which multiple factors can contribute to the phenotypic characteristics of TFH cells. Because type I IFNs are often generated in immune responses, we set out to investigate whether these factors are relevant to TFH cell differentiation. Type I IFNs promote Th1 responses, thus one possibility was these factors antagonized TFH-expressed genes. However, we show that type I IFNs (IFN-α/ß) induced B cell lymphoma 6 (Bcl6) expression, the master regulator transcription factor for TFH cells, and CXCR5 and programmed cell death-1 (encoded by Pdcd1), key surface molecules expressed by TFH cells. In contrast, type I IFNs failed to induce IL-21, the signature cytokine for TFH cells. The induction of Bcl6 was regulated directly by STAT1, which bound to the Bcl6, Cxcr5, and Pdcd1 loci. These data suggest that type I IFNs (IFN-α/ß) and STAT1 can contribute to some features of TFH cells but are inadequate in inducing complete programming of this subset.

Interleukin-27 priming of T cells controls IL-17 production in trans via induction of the ligand PD-L1.

Interleukin-27 (IL-27) is a key immunosuppressive cytokine that counters T helper 17 (Th17) cell-mediated pathology. To identify mechanisms by which IL-27 might exert its immunosuppressive effect, we analyzed genes in T cells rapidly induced by IL-27. We found that IL-27 priming of naive T cells upregulated expression of programmed death ligand 1 (PD-L1) in a signal transducer and activator of transcription 1 (STAT1)-dependent manner. When cocultured with naive CD4(+) T cells, IL-27-primed T cells inhibited the differentiation of Th17 cells in trans through a PD-1-PD-L1 interaction. In vivo, coadministration of naive TCR transgenic T cells (2D2 T cells) with IL-27-primed T cells expressing PD-L1 inhibited the development of Th17 cells and protected from severe autoimmune encephalomyelitis. Thus, these data identify a suppressive activity of IL-27, by which CD4(+) T cells can restrict differentiation of Th17 cells in trans.

Early Th1 cell differentiation is marked by a Tfh cell-like transition.

Follicular helper T (Tfh) cells comprise an important subset of helper T cells; however, their relationship with other helper lineages is incompletely understood. Herein, we showed interleukin-12 acting via the transcription factor STAT4 induced both Il21 and Bcl6 genes, generating cells with features of both Tfh and Th1 cells. However, STAT4 also induced the transcription factor T-bet. With ChIP-seq, we defined the genome-wide targets of T-bet and found that it repressed Bcl6 and other markers of Tfh cells, thereby attenuating the nascent Tfh cell-like phenotype in the late phase of Th1 cell specification. Tfh-like cells were rapidly generated after Toxoplasma gondii infection in mice, but T-bet constrained Tfh cell expansion and consequent germinal center formation and antibody production. Our data argue that Tfh and Th1 cells share a transitional stage through the signal mediated by STAT4, which promotes both phenotypes. However, T-bet represses Tfh cell functionalities, promoting full Th1 cell differentiation.

Optimal intervention strategies for cyclic therapeutic methods.

External control of a genetic regulatory network is used for the purpose of avoiding undesirable states such as those associated with a disease. Certain types of cancer therapies, such as chemotherapy, are given in cycles with each treatment being followed by a recovery period. During the recovery period, the side effects tend to gradually subside. In this paper, it is shown how an optimal cyclic intervention strategy can be devised for any Markovian genetic regulatory network. The effectiveness of optimal cyclic therapies is demonstrated through numerical studies for random networks. Furthermore, an optimal cyclic policy is derived to control the behavior of a regulatory model of the mammalian cell-cycle network.

Intervention in gene regulatory networks via a stationary mean-first-passage-time control policy.

A prime objective of modeling genetic regulatory networks is the identification of potential targets for therapeutic intervention. To date, optimal stochastic intervention has been studied in the context of probabilistic Boolean networks, with the control policy based on the transition probability matrix of the associated Markov chain and dynamic programming used to find optimal control policies. Dynamical programming algorithms are problematic owing to their high computational complexity. Two additional computationally burdensome issues that arise are the potential for controlling the network and identifying the best gene for intervention. This paper proposes an algorithm based on mean first-passage time that assigns a stationary control policy for each gene candidate. It serves as an approximation to an optimal control policy and, owing to its reduced computational complexity, can be used to predict the best control gene. Once the best control gene is identified, one can derive an optimal policy or simply utilize the approximate policy for this gene when the network size precludes a direct application of dynamic programming algorithms. A salient point is that the proposed algorithm can be model-free. It can be directly designed from time-course data without having to infer the transition probability matrix of the network.