Dissecting the dynamic transcriptional landscape of early T helper cell differentiation into Th1, Th2, and Th1/2 hybrid cells.

Selective differentiation of CD4+ T helper (Th) cells into specialized subsets such as Th1 and Th2 cells is a key element of the adaptive immune system driving appropriate immune responses. Besides those canonical Th-cell lineages, hybrid phenotypes such as Th1/2 cells arise in vivo, and their generation could be reproduced in vitro. While master-regulator transcription factors like T-bet for Th1 and GATA-3 for Th2 cells drive and maintain differentiation into the canonical lineages, the transcriptional architecture of hybrid phenotypes is less well understood. In particular, it has remained unclear whether a hybrid phenotype implies a mixture of the effects of several canonical lineages for each gene, or rather a bimodal behavior across genes. Th-cell differentiation is a dynamic process in which the regulatory factors are modulated over time, but longitudinal studies of Th-cell differentiation are sparse. Here, we present a dynamic transcriptome analysis following Th-cell differentiation into Th1, Th2, and Th1/2 hybrid cells at 3-h time intervals in the first hours after stimulation. We identified an early bifurcation point in gene expression programs, and we found that only a minority of ~20% of Th cell-specific genes showed mixed effects from both Th1 and Th2 cells on Th1/2 hybrid cells. While most genes followed either Th1- or Th2-cell gene expression, another fraction of ~20% of genes followed a Th1 and Th2 cell-independent transcriptional program associated with the transcription factors STAT1 and STAT4. Overall, our results emphasize the key role of high-resolution longitudinal data for the characterization of cellular phenotypes.

Heritable changes in division speed accompany the diversification of single T cell fate.

Rapid clonal expansion of antigen-specific T cells is a fundamental feature of adaptive immune responses. It enables the outgrowth of an individual T cell into thousands of clonal descendants that diversify into short-lived effectors and long-lived memory cells. Clonal expansion is thought to be programmed upon priming of a single naive T cell and then executed by homogenously fast divisions of all of its descendants. However, the actual speed of cell divisions in such an emerging "T cell family" has never been measured with single-cell resolution. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. This comprehensive mapping of T cell family trees identifies a short burst phase, in which division speed is homogenously fast and maintained independent of external cytokine availability or continued T cell receptor stimulation. Thereafter, however, division speed diversifies, and model-based computational analysis using a Bayesian inference framework for tree-structured data reveals a segregation into heritably fast- and slow-dividing branches. This diversification of division speed is preceded already during the burst phase by variable expression of the interleukin-2 receptor alpha chain. Later it is accompanied by selective expression of memory marker CD62L in slower dividing branches. Taken together, these data demonstrate that T cell clonal expansion is structured into subsequent burst and diversification phases, the latter of which coincides with specification of memory versus effector fate.

Selective multi-kinase inhibition sensitizes mesenchymal pancreatic cancer to immune checkpoint blockade by remodeling the tumor microenvironment.

KRAS-mutant pancreatic ductal adenocarcinoma (PDAC) is highly immunosuppressive and resistant to targeted and immunotherapies. Among the different PDAC subtypes, basal-like mesenchymal PDAC, which is driven by allelic imbalance, increased gene dosage and subsequent high expression levels of oncogenic KRAS, shows the most aggressive phenotype and strongest therapy resistance. In the present study, we performed a systematic high-throughput combination drug screen and identified a synergistic interaction between the MEK inhibitor trametinib and the multi-kinase inhibitor nintedanib, which targets KRAS-directed oncogenic signaling in mesenchymal PDAC. This combination treatment induces cell-cycle arrest and cell death, and initiates a context-dependent remodeling of the immunosuppressive cancer cell secretome. Using a combination of single-cell RNA-sequencing, CRISPR screens and immunophenotyping, we show that this combination therapy promotes intratumor infiltration of cytotoxic and effector T cells, which sensitizes mesenchymal PDAC to PD-L1 immune checkpoint inhibition. Overall, our results open new avenues to target this aggressive and therapy-refractory mesenchymal PDAC subtype.

Brief homogeneous TCR signals instruct common iNKT progenitors whose effector diversification is characterized by subsequent cytokine signaling.

Innate-like T cell populations expressing conserved TCRs play critical roles in immunity through diverse developmentally acquired effector functions. Focusing on the prototypical lineage of invariant natural killer T (iNKT) cells, we sought to dissect the mechanisms and timing of fate decisions and functional effector differentiation. Utilizing induced expression of the semi-invariant NKT cell TCR on double positive thymocytes, an initially highly synchronous wave of iNKT cell development was triggered by brief homogeneous TCR signaling. After reaching a uniform progenitor state characterized by IL-4 production potential and proliferation, effector subsets emerged simultaneously, but then diverged toward different fates. While NKT17 specification was quickly completed, NKT1 cells slowly differentiated and expanded. NKT2 cells resembled maturing progenitors, which gradually diminished in numbers. Thus, iNKT subset diversification occurs in dividing progenitor cells without acute TCR input but utilizes multiple active cytokine signaling pathways. These data imply a two-step model of iNKT effector differentiation.

Fate mapping of single NK cells identifies a type 1 innate lymphoid-like lineage that bridges innate and adaptive recognition of viral infection.

Upon viral infection, natural killer (NK) cells expressing certain germline-encoded receptors are selected, expanded, and maintained in an adaptive-like manner. Currently, these are thought to differentiate along a common pathway. However, by fate mapping of single NK cells upon murine cytomegalovirus (MCMV) infection, we identified two distinct NK cell lineages that contributed to adaptive-like responses. One was equivalent to conventional NK (cNK) cells while the other was transcriptionally similar to type 1 innate lymphoid cells (ILC1s). ILC1-like NK cells showed splenic residency and strong cytokine production but also recognized and killed MCMV-infected cells, guided by activating receptor Ly49H. Moreover, they induced clustering of conventional type 1 dendritic cells and facilitated antigen-specific T cell priming early during MCMV infection, which depended on Ly49H and the NK cell-intrinsic expression of transcription factor Batf3. Thereby, ILC1-like NK cells bridge innate and adaptive viral recognition and unite critical features of cNK cells and ILC1s.

Early emergence of T central memory precursors programs clonal dominance during chronic viral infection.

Chronic cytomegalovirus (CMV) infection leads to long-term maintenance of extraordinarily large CMV-specific T cell populations. The magnitude of this so-called 'memory inflation' is thought to mainly depend on antigenic stimulation during the chronic phase of infection. However, by mapping the long-term development of CD8+ T cell families derived from single naive precursors, we find that fate decisions made during the acute phase of murine CMV infection can alter the level of memory inflation by more than 1,000-fold. Counterintuitively, a T cell family's capacity for memory inflation is not determined by its initial expansion. Instead, those rare T cell families that dominate the chronic phase of infection show an early transcriptomic signature akin to that of established T central memory cells. Accordingly, a T cell family's long-term dominance is best predicted by its early content of T central memory precursors, which later serve as a stem-cell-like source for memory inflation.

The CMV-Specific CD8+ T Cell Response Is Dominated by Supra-Public Clonotypes with High Generation Probabilities.

Evolutionary processes govern the selection of T cell clonotypes that are optimally suited to mediate efficient antigen-specific immune responses against pathogens and tumors. While the theoretical diversity of T cell receptor (TCR) sequences is vast, the antigen-specific TCR repertoire is restricted by its peptide epitope and the presenting major histocompatibility complex (pMHC). It remains unclear how many TCR sequences are recruited into an antigen-specific T cell response, both within and across different organisms, and which factors shape both of these distributions. Infection of mice with ovalbumin-expressing cytomegalovirus (IE2-OVA-mCMV) represents a well-studied model system to investigate T cell responses given their size and longevity. Here we investigated > 180,000 H2kb/SIINFEKL-recognizing TCR CDR3α or CDR3ß sequences from 25 individual mice spanning seven different time points during acute infection and memory inflation. In-depth repertoire analysis revealed that from a pool of highly diverse, but overall limited sequences, T cell responses were dominated by public clonotypes, partly with unexpectedly extreme degrees of sharedness between individual mice ("supra-public clonotypes"). Public clonotypes were found exclusively in a fraction of TCRs with a high generation probability. Generation probability and degree of sharedness select for highly functional TCRs, possibly mediated through elevating intraindividual precursor frequencies of clonotypes.

The Outcome of Ex Vivo TIL Expansion Is Highly Influenced by Spatial Heterogeneity of the Tumor T-Cell Repertoire and Differences in Intrinsic In Vitro Growth Capacity between T-Cell Clones.

PURPOSE: During our efforts to develop tumor-infiltrating lymphocyte (TIL) therapy to counter the devastating recurrence rate in patients with primary resectable pancreatic ductal adenocarcinoma (PDA), we found that PDA TILs can readily be expanded in vitro and that the majority of resulting TIL cultures show reactivity against the autologous tumor. However, the fraction of tumor-reactive T cells is low. We investigated to which extent this was related to the in vitro expansion. EXPERIMENTAL DESIGN: We compared the clonal composition of TIL preparations before and after in vitro expansion using T-cell receptor (TCR) deep sequencing. Our findings for PDA were benchmarked to experiments with melanoma TILs. RESULTS: We found that the TIL TCR repertoire changes dramatically during in vitro expansion, leading to loss of tumor- dominant T-cell clones and overgrowth by newly emerging T-cell clones that are barely detectable in the tumor. These changes are primarily driven by differences in the intrinsic in vitro expansion capacity of T-cell clones. Single-cell experiments showed an association between poor proliferative capacity and expression of markers related to antigen experience and dysfunction. Furthermore, we found that spatial heterogeneity of the TIL repertoire resulted in TCR repertoires that are greatly divergent between TIL cultures derived from distant tumor samples of the same patient. CONCLUSIONS: Culture-induced changes in clonal composition are likely to affect tumor reactivity of TIL preparations. TCR deep sequencing provides important insights into the factors that govern the outcome of in vitro TIL expansion and thereby a path toward optimization of the production of TIL preparations with high therapeutic efficacy.See related commentary by Lozano-Rabella and Gros, p. 4177.

Differential expansion of T central memory precursor and effector subsets is regulated by division speed.

While antigen-primed T cells proliferate at speeds close to the physiologic maximum of mammalian cells, T cell memory is maintained in the absence of antigen by rare cell divisions. The transition between these distinct proliferative programs has been difficult to resolve via population-based analyses. Here, we computationally reconstruct the proliferative history of single CD8+ T cells upon vaccination and measure the division speed of emerging T cell subsets in vivo. We find that slower cycling central memory precursors, characterized by an elongated G1 phase, segregate early from the bulk of rapidly dividing effector subsets, and further slow-down their cell cycle upon premature removal of antigenic stimuli. In contrast, curtailed availability of inflammatory stimuli selectively restrains effector T cell proliferation due to reduced receptivity for interleukin-2. In line with these findings, persistence of antigenic but not inflammatory stimuli throughout clonal expansion critically determines the later size of the memory compartment.

Concomitant Infection of S. mansoni and H. pylori Promotes Promiscuity of Antigen-Experienced Cells and Primes the Liver for a Lower Fibrotic Response.

Helicobacter pylori chronically colonizes the stomach and is strongly associated with gastric cancer. Its concomitant occurrence with helminths such as schistosomes has been linked to reduced cancer incidence, presumably due to suppression of H. pylori-associated pro-inflammatory responses. However, experimental evidence in support of such a causal link or the mutual interaction of both pathogens is lacking. We investigated the effects of co-infection during the different immune phases of S. mansoni infection. Surprisingly, co-infected mice had increased H. pylori gastric colonization during the interferon gamma (IFNγ) phase of schistosome infection but reduced infiltration of T cells in the stomach due to misdirection of antigen-experienced CXCR3+ T cells to the liver. Unexpectedly, H. pylori co-infection resulted in partial protection from schistosome-induced liver damage. Here, we demonstrate that an increase in fibrosis-protective IL-13Ra2 is associated with H. pylori infection. Thus, our study strongly points to an immunological interaction of anatomically isolated pathogens, eventually resulting in altered disease pathology.

Distinct Surface Expression of Activating Receptor Ly49H Drives Differential Expansion of NK Cell Clones upon Murine Cytomegalovirus Infection.

Natural killer (NK) cells show some features of adaptive immunity but have not been studied at the clonal level. Here, we used retrogenic color-barcoding and single-cell adoptive transfers to track clonal immune responses to murine cytomegalovirus (MCMV) infection, derived from individual NK cells expressing activating receptor Ly49H. Clonal expansion of single NK cells varied substantially, and this variation could not be attributed to the additional presence or absence of inhibitory Ly49 receptors. Instead, single-cell-derived variability correlated with distinct surface expression levels of Ly49H itself. Ly49Hhi NK cell clones maintained higher Ly49H expression and expanded more than their Ly49Hlo counterparts in response to MCMV. Thus, akin to adaptive processes shaping an antigen-specific T cell receptor (TCR) repertoire, the Ly49H+ NK cell population adapts to MCMV infection. This process relies on the clonal maintenance of distinct Ly49H expression levels, generating a repertoire of individual NK cells outfitted with distinct reactivity to MCMV.

Single-Cell Resolution of T Cell Immune Responses.

Single antigen-specific B or T lymphocytes are the smallest functional units, into which an adaptive immune response can be dissected. Today, novel high-throughput technologies are providing researches with increasingly complex information on the diverse phenotypic signatures of individual lymphocytes. With a focus on T cells, we summarize here, how computational approaches are becoming increasingly important to identify the relevant developmental boundaries and connections between these high-dimensional lymphocyte states. We then describe how these insights may be further expanded by novel experimental approaches that allow to map the fate of individual T cells and their progeny in vivo and in vitro. Finally, we highlight how these experiments have uncovered a probabilistic regulatory structure of T cell immune responses and briefly discuss, how two distinct theoretical frameworks used to describe this structure may be merged to best capture single T cell behavior in computational terms.

TCR Signal Quality Modulates Fate Decisions of Single CD4+ T Cells in a Probabilistic Manner.

To what extent the lineage decisions of activated CD4+ T cells are determined by the quality of T cell receptor (TCR) ligation is incompletely understood. Here, we show that individual T cells expressing identical TCRs take highly variable fate decisions despite binding the same ligand. We identify a mathematical model that correctly captures this probabilistic behavior and allows one to formalize changes in TCR signal quality-due to cognate versus altered peptide ligation-as changes of lineage-specific proliferation and differentiation rates. We show that recall responses also adhere to this probabilistic framework requiring recruitment of multiple memory clones to provide reliable differentiation patterns. By extending our framework to simulate hypothetical TCRs of distinct binding strength, we reconstruct primary and secondary response patterns emerging from a polyclonal TCR repertoire in silico. Collectively, these data suggest that individual T cells harboring distinct TCRs generate overlapping primary differentiation patterns that segregate only upon repetitive immunization.

Next-generation TCR sequencing - a tool to understand T-cell infiltration in human cancers.

Tumour-infiltrating lymphocytes (TILs) are known to mediate potent anti-tumour activity. As T-cell-based therapies start to reach clinical practice, it becomes increasingly important to understand what characterizes a successful anti-tumour T-cell response and to exploit this knowledge for patient stratification. Next-generation sequencing of T-cell receptors (TCRs) promises to provide insights into the complexity of the tumour T-cell infiltrate that go beyond the phenotypic level. A recent study by Chen et al made use of this novel technology to demonstrate that the TIL repertoire of oesophageal squamous cell carcinoma patients is distinct from that of non-tumour sites and is characterized by significant intratumoural heterogeneity. This study illustrates the great potential of the method and addresses several technical and biological hurdles that need to be considered. Careful sampling, normalization, and error correction will be required to optimally use TCR sequencing to answer biological questions and define predictive biomarkers, e.g. for cancer immunotherapy. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Stem-cell dynamics and lineage topology from in vivo fate mapping in the hematopoietic system.

In recent years, sophisticated fate-mapping tools have been developed to study the behavior of stem cells in the intact organism. These experimental approaches are beginning to yield a quantitative picture of how cell numbers are regulated during steady state and in response to challenges. Focusing on hematopoiesis and immune responses, we discuss how novel mathematical approaches driven by these fate-mapping data have provided insights into the dynamics and topology of cellular differentiation pathways in vivo. The combination of experiment and theory has allowed to quantify the degree of self-renewal in stem and progenitor cells, shown how native hematopoiesis differs fundamentally from post-transplantation hematopoiesis, and uncovered that the diversification of T lymphocytes during immune responses resembles tissue renewal driven by stem cells.

Identification of a tumor-reactive T-cell repertoire in the immune infiltrate of patients with resectable pancreatic ductal adenocarcinoma.

PURPOSE: The devastating prognosis of patients with resectable pancreatic ductal adenocarcinoma (PDA) presents an urgent need for the development of therapeutic strategies targeting disseminated tumor cells. Until now, T-cell therapy has been scarcely pursued in PDA, due to the prevailing view that it represents a poorly immunogenic tumor. EXPERIMENTAL DESIGN: We systematically analyzed T-cell infiltrates in tumor biopsies from 127 patients with resectable PDA by means of immunohistochemistry, flow cytometry, T-cell receptor (TCR) deep-sequencing and functional analysis of in vitro expanded T-cell cultures. Parallel studies were performed on tumor-infiltrating lymphocytes (TIL) from 44 patients with metastatic melanoma. RESULTS: Prominent T-cell infiltrates, as well as tertiary lymphoid structures harboring proliferating T-cells, were detected in the vast majority of biopsies from PDA patients. The notion that the tumor is a site of local T-cell expansion was strengthened by TCR deep-sequencing, revealing that the T-cell repertoire in the tumor is dominated by highly frequent CDR3 sequences that can be up to 10,000-fold enriched in tumor as compared to peripheral blood. In fact, TCR repertoire composition in PDA resembled that in melanoma. Moreover, in vitro expansion of TILs was equally efficient for PDA and melanoma, resulting in T-cell cultures displaying HLA class I-restricted reactivity against autologous tumor cells. CONCLUSIONS: The tumor-infiltrating T-cell response in PDA shows striking similarity to that in melanoma, where adoptive T-cell therapy has significant therapeutic impact. Our findings indicate that T-cell-based therapies may be used to counter disease recurrence in patients with resectable PDA.

Fundamental properties of unperturbed haematopoiesis from stem cells in vivo.

Haematopoietic stem cells (HSCs) are widely studied by HSC transplantation into immune- and blood-cell-depleted recipients. Single HSCs can rebuild the system after transplantation. Chromosomal marking, viral integration and barcoding of transplanted HSCs suggest that very low numbers of HSCs perpetuate a continuous stream of differentiating cells. However, the numbers of productive HSCs during normal haematopoiesis, and the flux of differentiating progeny remain unknown. Here we devise a mouse model allowing inducible genetic labelling of the most primitive Tie2(+) HSCs in bone marrow, and quantify label progression along haematopoietic development by limiting dilution analysis and data-driven modelling. During maintenance of the haematopoietic system, at least 30% or ∼5,000 HSCs are productive in the adult mouse after label induction. However, the time to approach equilibrium between labelled HSCs and their progeny is surprisingly long, a time scale that would exceed the mouse's life. Indeed, we find that adult haematopoiesis is largely sustained by previously designated 'short-term' stem cells downstream of HSCs that nearly fully self-renew, and receive rare but polyclonal HSC input. By contrast, in fetal and early postnatal life, HSCs are rapidly used to establish the immune and blood system. In the adult mouse, 5-fluoruracil-induced leukopenia enhances the output of HSCs and of downstream compartments, thus accelerating haematopoietic flux. Label tracing also identifies a strong lineage bias in adult mice, with several-hundred-fold larger myeloid than lymphoid output, which is only marginally accentuated with age. Finally, we show that transplantation imposes severe constraints on HSC engraftment, consistent with the previously observed oligoclonal HSC activity under these conditions. Thus, we uncover fundamental differences between the normal maintenance of the haematopoietic system, its regulation by challenge, and its re-establishment after transplantation. HSC fate mapping and its linked modelling provide a quantitative framework for studying in situ the regulation of haematopoiesis in health and disease.

Individual T helper cells have a quantitative cytokine memory.

The probabilistic expression of cytokine genes in differentiated T helper (Th) cell populations remains ill defined. By single-cell analyses and mathematical modeling, we show that one stimulation featured stable cytokine nonproducers as well as stable producers with wide cell-to-cell variability in the magnitude of expression. Focusing on interferon-γ (IFN-γ) expression by Th1 cells, mathematical modeling predicted that this behavior reflected different cell-intrinsic capacities and not mere gene-expression noise. In vivo, Th1 cells sort purified by secreted IFN-γ amounts preserved a quantitative memory for both probability and magnitude of IFN-γ re-expression for at least 1 month. Mechanistically, this memory resulted from quantitatively distinct transcription of individual alleles and was controlled by stable expression differences of the Th1 cell lineage-specifying transcription factor T-bet. Functionally, Th1 cells with graded IFN-γ production competence differentially activated Salmonella-infected macrophages for bacterial killing. Thus, individual Th cells commit to produce distinct amounts of a given cytokine, thereby generating functional intrapopulation heterogeneity.