Metabolic modeling of single Th17 cells reveals regulators of autoimmunity.

Metabolism is a major regulator of immune cell function, but it remains difficult to study the metabolic status of individual cells. Here, we present Compass, an algorithm to characterize cellular metabolic states based on single-cell RNA sequencing and flux balance analysis. We applied Compass to associate metabolic states with T helper 17 (Th17) functional variability (pathogenic potential) and recovered a metabolic switch between glycolysis and fatty acid oxidation, akin to known Th17/regulatory T cell (Treg) differences, which we validated by metabolic assays. Compass also predicted that Th17 pathogenicity was associated with arginine and downstream polyamine metabolism. Indeed, polyamine-related enzyme expression was enhanced in pathogenic Th17 and suppressed in Treg cells. Chemical and genetic perturbation of polyamine metabolism inhibited Th17 cytokines, promoted Foxp3 expression, and remodeled the transcriptome and epigenome of Th17 cells toward a Treg-like state. In vivo perturbations of the polyamine pathway altered the phenotype of encephalitogenic T cells and attenuated tissue inflammation in CNS autoimmunity.

B-cell-directed CAR T-cell therapy activates CD8+ cytotoxic CARneg bystander T cells in patients and nonhuman primates.

ABSTRACT: Chimeric antigen receptor (CAR) T cells hold promise as a therapy for B-cell-derived malignancies, and despite their impressive initial response rates, a significant proportion of patients ultimately experience relapse. Although recent studies have explored the mechanisms of in vivo CAR T-cell function, little is understood about the activation of surrounding CARneg bystander T cells and their potential to enhance tumor responses. We performed single-cell RNA sequencing on nonhuman primate (NHP) and patient-derived T cells to identify the phenotypic and transcriptomic hallmarks of bystander activation of CARneg T cells following B-cell-targeted CAR T-cell therapy. Using a highly translatable CD20 CAR NHP model, we observed a distinct population of activated CD8+ CARneg T cells emerging during CAR T-cell expansion. These bystander CD8+ CARneg T cells exhibited a unique transcriptional signature with upregulation of natural killer-cell markers (KIR3DL2, CD160, and KLRD1), chemokines, and chemokine receptors (CCL5, XCL1, and CCR9), and downregulation of naïve T-cell-associated genes (SELL and CD28). A transcriptionally similar population was identified in patients after a tisagenlecleucel infusion. Mechanistic studies revealed that interleukin-2 (IL-2) and IL-15 exposure induced bystander-like CD8+ T cells in a dose-dependent manner. In vitro activated and patient-derived T cells with a bystander phenotype efficiently killed leukemic cells through a T-cell receptor-independent mechanism. Collectively, to our knowledge, these data provide the first comprehensive identification and profiling of CARneg bystander CD8+ T cells following B-cell-targeting CAR T-cell therapy and suggest a novel mechanism through which CAR T-cell infusion might trigger enhanced antileukemic responses. Patient samples were obtained from the trial #NCT03369353, registered at www.ClinicalTrials.gov.

High-Specificity Targeted Functional Profiling in Microbial Communities with ShortBRED.

Profiling microbial community function from metagenomic sequencing data remains a computationally challenging problem. Mapping millions of DNA reads from such samples to reference protein databases requires long run-times, and short read lengths can result in spurious hits to unrelated proteins (loss of specificity). We developed ShortBRED (Short, Better Representative Extract Dataset) to address these challenges, facilitating fast, accurate functional profiling of metagenomic samples. ShortBRED consists of two components: (i) a method that reduces reference proteins of interest to short, highly representative amino acid sequences ("markers") and (ii) a search step that maps reads to these markers to quantify the relative abundance of their associated proteins. After evaluating ShortBRED on synthetic data, we applied it to profile antibiotic resistance protein families in the gut microbiomes of individuals from the United States, China, Malawi, and Venezuela. Our results support antibiotic resistance as a core function in the human gut microbiome, with tetracycline-resistant ribosomal protection proteins and Class A beta-lactamases being the most widely distributed resistance mechanisms worldwide. ShortBRED markers are applicable to other homology-based search tasks, which we demonstrate here by identifying phylogenetic signatures of antibiotic resistance across more than 3,000 microbial isolate genomes. ShortBRED can be applied to profile a wide variety of protein families of interest; the software, source code, and documentation are available for download at http://huttenhower.sph.harvard.edu/shortbred.

CXCR3 regulates stem and proliferative CD8+ T cells during chronic infection by promoting interactions with DCs in splenic bridging channels.

Production of effector CD8+ T cells during persistent infection requires a stable pool of stem-like cells that can give rise to effector cells via a proliferative intermediate population. In infection models marked by T cell exhaustion, this process can be transiently induced by checkpoint blockade but occurs spontaneously in mice chronically infected with the protozoan intracellular parasite Toxoplasma gondii. We observe distinct locations for parasite-specific T cell subsets, implying a link between differentiation and anatomical niches in the spleen. Loss of the chemokine receptor CXCR3 on T cells does not prevent white pulp-to-red pulp migration but reduces interactions with CXCR3 ligand-producing dendritic cells (DCs) and impairs memory-to-intermediate transition, leading to a buildup of memory T cells in the red pulp. Thus, CXCR3 increases T cell exposure to differentiation-inducing signals during red pulp migration, providing a dynamic mechanism for modulating effector differentiation in response to environmental signals.

Batf-mediated epigenetic control of effector CD8+ T cell differentiation.

The response of naive CD8+ T cells to their cognate antigen involves rapid and broad changes to gene expression that are coupled with extensive chromatin remodeling, but the mechanisms governing these changes are not fully understood. Here, we investigated how these changes depend on the basic leucine zipper ATF-like transcription factor Batf, which is essential for the early phases of the process. Through genome scale profiling, we characterized the role of Batf in chromatin organization at several levels, including the accessibility of key regulatory regions, the expression of their nearby genes, and the interactions that these regions form with each other and with key transcription factors. We identified a core network of transcription factors that cooperated with Batf, including Irf4, Runx3, and T-bet, as indicated by their colocalization with Batf and their binding in regions whose accessibility, interactions, and expression of nearby genes depend on Batf. We demonstrated the synergistic activity of this network by overexpressing the different combinations of these genes in fibroblasts. Batf and Irf4, but not Batf alone, were sufficient to increase accessibility and transcription of key loci, normally associated with T cell function. Addition of Runx3 and T-bet further contributed to fine-tuning of these changes and was essential for establishing chromatin loops characteristic of T cells. These data provide a resource for studying the epigenomic and transcriptomic landscape of effector differentiation of cytotoxic T cells and for investigating the interdependency between transcription factors and its effects on the epigenome and transcriptome of primary cells.

Identification and Tracking of Alloreactive T Cell Clones in Rhesus Macaques Through the RM-scTCR-Seq Platform.

T cell receptor (TCR) clonotype tracking is a powerful tool for interrogating T cell mediated immune processes. New methods to pair a single cell's transcriptional program with its TCR identity allow monitoring of T cell clonotype-specific transcriptional dynamics. While these technologies have been available for human and mouse T cells studies, they have not been developed for Rhesus Macaques (RM), a critical translational organism for autoimmune diseases, vaccine development and transplantation. We describe a new pipeline, 'RM-scTCR-Seq', which, for the first time, enables RM specific single cell TCR amplification, reconstruction and pairing of RM TCR's with their transcriptional profiles. We apply this method to a RM model of GVHD, and identify and track in vitro detected alloreactive clonotypes in GVHD target organs and explore their GVHD driven cytotoxic T cell signature. This novel, state-of-the-art platform fundamentally advances the utility of RM to study protective and pathogenic T cell responses.

T cell self-reactivity during thymic development dictates the timing of positive selection.

Functional tuning of T cells based on their degree of self-reactivity is established during positive selection in the thymus, although how positive selection differs for thymocytes with relatively low versus high self-reactivity is unclear. In addition, preselection thymocytes are highly sensitive to low-affinity ligands, but the mechanism underlying their enhanced T cell receptor (TCR) sensitivity is not fully understood. Here we show that murine thymocytes with low self-reactivity experience briefer TCR signals and complete positive selection more slowly than those with high self-reactivity. Additionally, we provide evidence that cells with low self-reactivity retain a preselection gene expression signature as they mature, including genes previously implicated in modulating TCR sensitivity and a novel group of ion channel genes. Our results imply that thymocytes with low self-reactivity downregulate TCR sensitivity more slowly during positive selection, and associate membrane ion channel expression with thymocyte self-reactivity and progress through positive selection.

Spatiotemporal single-cell profiling reveals that invasive and tissue-resident memory donor CD8+ T cells drive gastrointestinal acute graft-versus-host disease.

Organ infiltration by donor T cells is critical to the development of acute graft-versus-host disease (aGVHD) in recipients after allogeneic hematopoietic stem cell transplant (allo-HCT). However, deconvoluting the transcriptional programs of newly recruited donor T cells from those of tissue-resident T cells in aGVHD target organs remains a challenge. Here, we combined the serial intravascular staining technique with single-cell RNA sequencing to dissect the tightly connected processes by which donor T cells initially infiltrate tissues and then establish a pathogenic tissue residency program in a rhesus macaque allo-HCT model that develops aGVHD. Our results enabled creation of a spatiotemporal map of the transcriptional programs controlling donor CD8+ T cell infiltration into the primary aGVHD target organ, the gastrointestinal (GI) tract. We identified the large and small intestines as the only two sites demonstrating allo-specific, rather than lymphodepletion-driven, T cell infiltration. GI-infiltrating donor CD8+ T cells demonstrated a highly activated, cytotoxic phenotype while simultaneously developing a canonical tissue-resident memory T cell (TRM) transcriptional signature driven by interleukin-15 (IL-15)/IL-21 signaling. We found expression of a cluster of genes directly associated with tissue invasiveness, including those encoding adhesion molecules (ITGB2), specific chemokines (CCL3 and CCL4L1) and chemokine receptors (CD74), as well as multiple cytoskeletal proteins. This tissue invasion transcriptional signature was validated by its ability to discriminate the CD8+ T cell transcriptome of patients with GI aGVHD from those of GVHD-free patients. These results provide insights into the mechanisms controlling tissue occupancy of target organs by pathogenic donor CD8+ TRM cells during aGVHD in primate transplant recipients.

The epigenetic landscape of T cell exhaustion.

Exhausted T cells in cancer and chronic viral infection express distinctive patterns of genes, including sustained expression of programmed cell death protein 1 (PD-1). However, the regulation of gene expression in exhausted T cells is poorly understood. Here, we define the accessible chromatin landscape in exhausted CD8+ T cells and show that it is distinct from functional memory CD8+ T cells. Exhausted CD8+ T cells in humans and a mouse model of chronic viral infection acquire a state-specific epigenetic landscape organized into functional modules of enhancers. Genome editing shows that PD-1 expression is regulated in part by an exhaustion-specific enhancer that contains essential RAR, T-bet, and Sox3 motifs. Functional enhancer maps may offer targets for genome editing that alter gene expression preferentially in exhausted CD8+ T cells.

Improving the oral efficacy of CNS drug candidates: discovery of highly orally efficacious piperidinyl piperidine M2 muscarinic receptor antagonists.

In search of a backup M(2) muscarinic receptor antagonist to the previously reported compound 1, we discovered compound (+)-14, which showed superior oral efficacy in animal models. The improvement of oral efficacy was achieved by modulating both the molecular weight and lipophilicity of the lead compounds.

Assessing the burden of diabetes mellitus in emergency departments in the United States: the National Hospital Ambulatory Medical Care Survey (NHAMCS).

OBJECTIVE: To evaluate the performance of three alternative methods to identify diabetes in patients visiting Emergency Departments (EDs), and to describe the characteristics of patients with diabetes who are not identified when the alternative methods are used. RESEARCH DESIGN AND METHODS: We used data from the National Hospital Ambulatory Medical Care Survey (NHAMCS) 2009 and 2010. We assessed the sensitivity and specificity of using providers' diagnoses and diabetes medications (both excluding and including biguanides) to identify diabetes compared to using the checkbox for diabetes as the gold standard. We examined the characteristics of patients whose diabetes was missed using multivariate Poisson regression models. RESULTS: The checkbox identified 5,567 ED visits by adult patients with diabetes. Compared to the checkbox, the sensitivity was 12.5% for providers' diagnoses alone, 20.5% for providers' diagnoses and diabetes medications excluding biguanides, and 21.5% for providers' diagnoses and diabetes medications including biguanides. The specificity of all three of the alternative methods was >99%. Older patients were more likely to have diabetes not identified. Patients with self-payment, those who had glucose measured or received IV fluids in the ED, and those with more diagnosis codes and medications, were more likely to have diabetes identified. CONCLUSIONS: NHAMCS's providers' diagnosis codes and medication lists do not identify the majority of patients with diabetes visiting EDs. The newly introduced checkbox is helpful in measuring ED resource utilization by patients with diabetes.

The active site of a zinc-dependent metalloproteinase influences the computed pK(a) of ligands coordinated to the catalytic zinc ion.

TNF-alpha converting enzyme (TACE) is a multidomain, membrane-anchored protein that includes a Zn-dependent protease domain. It releases the soluble form of cytokine tumor necrosis factor-alpha (TNF-alpha) from its membrane-bound precursor. TACE is a metalloprotease containing a catalytic glutamic acid, Glu-406, and a Zn(2+) ion ligated to three imidazoles. The protonation states of the active site glutamic acid and inhibitors are important factors in understanding the potency of inhibitors with acidic zinc-ligating groups such as hydroxamic and carboxylic acids. Density functional methods were utilized to compute pK(a) values using a model of the catalytic site of TACE and to predict a concomitant mechanism of binding, consistent with lowering the pK(a) of the bound ligand and raising the pK(a) of the active site Glu-406. Weak acids, such as hydroxamic acids, bind in their neutral form and then transfer an acidic proton to Glu-406. Stronger acids, such as carboxylic acids, bind in their anionic form and require preprotonation of Glu-406. Similar binding events would be expected for other zinc-dependent proteases.