Peptide-specific recognition of human cytomegalovirus strains controls adaptive natural killer cells.

Natural killer (NK) cells are innate lymphocytes that lack antigen-specific rearranged receptors, a hallmark of adaptive lymphocytes. In some people infected with human cytomegalovirus (HCMV), an NK cell subset expressing the activating receptor NKG2C undergoes clonal-like expansion that partially resembles anti-viral adaptive responses. However, the viral ligand that drives the activation and differentiation of adaptive NKG2C+ NK cells has remained unclear. Here we found that adaptive NKG2C+ NK cells differentially recognized distinct HCMV strains encoding variable UL40 peptides that, in combination with pro-inflammatory signals, controlled the population expansion and differentiation of adaptive NKG2C+ NK cells. Thus, we propose that polymorphic HCMV peptides contribute to shaping of the heterogeneity of adaptive NKG2C+ NK cell populations among HCMV-seropositive people.

Resident memory CD4+ T lymphocytes mobilize from bone marrow to contribute to a systemic secondary immune reaction.

Resident memory T lymphocytes (TRM ) of epithelial tissues and the Bm protect their host tissue. To what extent these cells are mobilized and contribute to systemic immune reactions is less clear. Here, we show that in secondary immune reactions to the measles-mumps-rubella (MMR) vaccine, CD4+ TRM are mobilized into the blood within 16 to 48 h after immunization in humans. This mobilization of TRM is cognate: TRM recognizing other antigens are not mobilized, unless they cross-react with the vaccine. We also demonstrate through methylome analyses that TRM are mobilized from the Bm. These mobilized cells make significant contribution to the systemic immune reaction, as evidenced by their T-cell receptor Vß clonotypes represented among the newly generated circulating memory T-cells, 14 days after vaccination. Thus, TRM of the Bm confer not only local, but also systemic immune memory.

Principles for circadian orchestration of metabolic pathways.

Circadian rhythms govern multiple aspects of animal metabolism. Transcriptome-, proteome- and metabolome-wide measurements have revealed widespread circadian rhythms in metabolism governed by a cellular genetic oscillator, the circadian core clock. However, it remains unclear if and under which conditions transcriptional rhythms cause rhythms in particular metabolites and metabolic fluxes. Here, we analyzed the circadian orchestration of metabolic pathways by direct measurement of enzyme activities, analysis of transcriptome data, and developing a theoretical method called circadian response analysis. Contrary to a common assumption, we found that pronounced rhythms in metabolic pathways are often favored by separation rather than alignment in the times of peak activity of key enzymes. This property holds true for a set of metabolic pathway motifs (e.g., linear chains and branching points) and also under the conditions of fast kinetics typical for metabolic reactions. By circadian response analysis of pathway motifs, we determined exact timing separation constraints on rhythmic enzyme activities that allow for substantial rhythms in pathway flux and metabolite concentrations. Direct measurements of circadian enzyme activities in mouse skeletal muscle confirmed that such timing separation occurs in vivo.

Membrane Tension Acts Through PLD2 and mTORC2 to Limit Actin Network Assembly During Neutrophil Migration.

For efficient polarity and migration, cells need to regulate the magnitude and spatial distribution of actin assembly. This process is coordinated by reciprocal interactions between the actin cytoskeleton and mechanical forces. Actin polymerization-based protrusion increases tension in the plasma membrane, which in turn acts as a long-range inhibitor of actin assembly. These interactions form a negative feedback circuit that limits the magnitude of membrane tension in neutrophils and prevents expansion of the existing front and the formation of secondary fronts. It has been suggested that the plasma membrane directly inhibits actin assembly by serving as a physical barrier that opposes protrusion. Here we show that efficient control of actin polymerization-based protrusion requires an additional mechanosensory feedback cascade that indirectly links membrane tension with actin assembly. Specifically, elevated membrane tension acts through phospholipase D2 (PLD2) and the mammalian target of rapamycin complex 2 (mTORC2) to limit actin nucleation. In the absence of this pathway, neutrophils exhibit larger leading edges, higher membrane tension, and profoundly defective chemotaxis. Mathematical modeling suggests roles for both the direct (mechanical) and indirect (biochemical via PLD2 and mTORC2) feedback loops in organizing cell polarity and motility-the indirect loop is better suited to enable competition between fronts, whereas the direct loop helps spatially organize actin nucleation for efficient leading edge formation and cell movement. This circuit is essential for polarity, motility, and the control of membrane tension.

Three-Dimensional Gradients of Cytokine Signaling between T Cells.

Immune responses are regulated by diffusible mediators, the cytokines, which act at sub-nanomolar concentrations. The spatial range of cytokine communication is a crucial, yet poorly understood, functional property. Both containment of cytokine action in narrow junctions between immune cells (immunological synapses) and global signaling throughout entire lymph nodes have been proposed, but the conditions under which they might occur are not clear. Here we analyze spatially three-dimensional reaction-diffusion models for the dynamics of cytokine signaling at two successive scales: in immunological synapses and in dense multicellular environments. For realistic parameter values, we observe local spatial gradients, with the cytokine concentration around secreting cells decaying sharply across only a few cell diameters. Focusing on the well-characterized T-cell cytokine interleukin-2, we show how cytokine secretion and competitive uptake determine this signaling range. Uptake is shaped locally by the geometry of the immunological synapse. However, even for narrow synapses, which favor intrasynaptic cytokine consumption, escape fluxes into the extrasynaptic space are expected to be substantial (≥20% of secretion). Hence paracrine signaling will generally extend beyond the synapse but can be limited to cellular microenvironments through uptake by target cells or strong competitors, such as regulatory T cells. By contrast, long-range cytokine signaling requires a high density of cytokine producers or weak consumption (e.g., by sparsely distributed target cells). Thus in a physiological setting, cytokine gradients between cells, and not bulk-phase concentrations, are crucial for cell-to-cell communication, emphasizing the need for spatially resolved data on cytokine signaling.

Derivation of Ca2+ signals from puff properties reveals that pathway function is robust against cell variability but sensitive for control.

Ca(2+) is a universal second messenger in eukaryotic cells transmitting information through sequences of concentration spikes. A prominent mechanism to generate these spikes involves Ca(2+) release from the endoplasmic reticulum Ca(2+) store via inositol 1,4,5-trisphosphate (IP(3))-sensitive channels. Puffs are elemental events of IP(3)-induced Ca(2+) release through single clusters of channels. Intracellular Ca(2+) dynamics are a stochastic system, but a complete stochastic theory has not been developed yet. We formulate the theory in terms of interpuff interval and puff duration distributions because, unlike the properties of individual channels, they can be measured in vivo. Our theory reproduces the typical spectrum of Ca(2+) signals like puffs, spiking, and bursting in analytically treatable test cases as well as in more realistic simulations. We find conditions for spiking and calculate interspike interval (ISI) distributions. Signal form, average ISI and ISI distributions depend sensitively on the details of cluster properties and their spatial arrangement. In contrast to that, the relation between the average and the standard deviation of ISIs does not depend on cluster properties and cluster arrangement and is robust with respect to cell variability. It is controlled by the global feedback processes in the Ca(2+) signaling pathway (e.g., via IP(3)-3-kinase or endoplasmic reticulum depletion). That relation is essential for pathway function because it ensures frequency encoding despite the randomness of ISIs and determines the maximal spike train information content. Hence, we find a division of tasks between global feedbacks and local cluster properties that guarantees robustness of function while maintaining sensitivity of control of the average ISI.

Fundamental properties of Ca2+ signals.

BACKGROUND: Ca2+ is a ubiquitous and versatile second messenger that transmits information through changes of the cytosolic Ca2+ concentration. Recent investigations changed basic ideas on the dynamic character of Ca2+ signals and challenge traditional ideas on information transmission. SCOPE OF REVIEW: We present recent findings on key characteristics of the cytosolic Ca2+ dynamics and theoretical concepts that explain the wide range of experimentally observed Ca2+ signals. Further, we relate properties of the dynamical regulation of the cytosolic Ca2+ concentration to ideas about information transmission by stochastic signals. MAJOR CONCLUSIONS: We demonstrate the importance of the hierarchal arrangement of Ca2+ release sites on the emergence of cellular Ca2+ spikes. Stochastic Ca2+ signals are functionally robust and adaptive to changing environmental conditions. Fluctuations of interspike intervals (ISIs) and the moment relation derived from ISI distributions contain information on the channel cluster open probability and on pathway properties. GENERAL SIGNIFICANCE: Robust and reliable signal transduction pathways that entail Ca2+ dynamics are essential for eukaryotic organisms. Moreover, we expect that the design of a stochastic mechanism which provides robustness and adaptivity will be found also in other biological systems. Ca2+ dynamics demonstrate that the fluctuations of cellular signals contain information on molecular behavior. This article is part of a Special Issue entitled Biochemical, biophysical and genetic approaches to intracellular calcium signaling.

Competing feedback loops shape IL-2 signaling between helper and regulatory T lymphocytes in cellular microenvironments.

Cytokines are pleiotropic and readily diffusible messenger molecules, raising the question of how their action can be confined to specific target cells. The T cell cytokine interleukin-2 (IL-2) is essential for the homeostasis of regulatory T (Treg) cells that suppress (auto)immunity and stimulates immune responses mediated by conventional T cells. We combined mathematical modeling and experiments to dissect the dynamics of the IL-2 signaling network that links the prototypical IL-2 producers, conventional T helper (Th) cells, and Treg cells. We show how the IL-2-induced upregulation of high-affinity IL-2 receptors (IL-2R) establishes a positive feedback loop of IL-2 signaling. This feedback mediates a digital switch for the proliferation of Th cells and functions as an analog amplifier for the IL-2 uptake capacity of Treg cells. Unlike other positive feedbacks in cell signaling that augment signal propagation, the IL-2/IL-2R loop enhances the capture of the signal molecule and its degradation. Thus Treg and Th cells can compete for IL-2 and restrict its range of action through efficient cellular uptake. Depending on activation status and spatial localization of the cells, IL-2 may be consumed exclusively by Treg or Th cells, or be shared between them. In particular, a Treg cell can deprive a stimulated Th cell of its IL-2, but only when the cells are located in close proximity, within a few tens of micrometers. The present findings explain how IL-2 can play two distinct roles in immune regulation and point to a hitherto largely unexplored spatiotemporal complexity of cytokine signaling.

Distribution modeling quantifies collective TH cell decision circuits in chronic inflammation.

Immune responses are tightly regulated by a diverse set of interacting immune cell populations. Alongside decision-making processes such as differentiation into specific effector cell types, immune cells initiate proliferation at the beginning of an inflammation, forming two layers of complexity. Here, we developed a general mathematical framework for the data-driven analysis of collective immune cell dynamics. We identified qualitative and quantitative properties of generic network motifs, and we specified differentiation dynamics by analysis of kinetic transcriptome data. Furthermore, we derived a specific, data-driven mathematical model for T helper 1 versus T follicular helper cell-fate decision dynamics in acute and chronic lymphocytic choriomeningitis virus infections in mice. The model recapitulates important dynamical properties without model fitting and solely by using measured response-time distributions. Model simulations predict different windows of opportunity for perturbation in acute and chronic infection scenarios, with potential implications for optimization of targeted immunotherapy.

Transcriptional profiling upon T cell stimulation reveals down-regulation of inflammatory pathways in T and B cells in SLE versus Sjögren's syndrome.

Systemic lupus erythematosus (SLE) and primary Sjögren's syndrome (pSS) share clinical as well as pathogenic similarities. Although previous studies suggest various abnormalities in different immune cell compartments, dedicated cell-type specific transcriptomic signatures are often masked by patient heterogeneity. Here, we performed transcriptional profiling of isolated CD4, CD8, CD16 and CD19 lymphocytes from pSS and SLE patients upon T cell stimulation, in addition to a steady-state condition directly after blood drawing, in total comprising 581 sequencing samples. T cell stimulation, which induced a pronounced inflammatory response in all four cell types, gave rise to substantial re-modulation of lymphocyte subsets in the two autoimmune diseases compared to healthy controls, far exceeding the transcriptomic differences detected at steady-state. In particular, we detected cell-type and disease-specific down-regulation of a range of pro-inflammatory cytokine and chemokine pathways. Such differences between SLE and pSS patients are instrumental for selective immune targeting by future therapies.

Calcium signaling: from single channels to pathways.

Ca(2+) is not only one of the most versatile and ubiquitous second messengers but also a well-established representative example of cell signaling. The identification of most key elements involved in Ca(2+) signaling enables a mechanistic and quantitative understanding of this particular pathway. Cellular behavior relies in general on the orchestration of molecular behavior leading to reliable cellular responses that allow for regulation and adaptation. Ca(2+) signaling uses a hierarchical organization to transform single molecule behavior into cell wide signals. We have recently shown experimentally that this organization carries single channel signatures onto the whole cell level and renders Ca(2+) oscillations stochastic. Here, we briefly review the co-evolution of experimental and theoretical studies in Ca(2+) -signaling and show how dynamic bottom-up modeling can be used to address -biological questions and illuminate biological principles of cell signaling.

Data-Driven Mathematical Model of Apoptosis Regulation in Memory Plasma Cells.

Memory plasma cells constitutively produce copious amounts of antibodies, imposing a critical risk factor for autoimmune disease. We previously found that plasma cell survival requires secreted factors such as APRIL and direct contact to stromal cells, which act in concert to activate NF-κB- and PI3K-dependent signaling pathways to prevent cell death. However, the regulatory properties of the underlying biochemical network are confounded by the complexity of potential interaction and cross-regulation pathways. Here, based on flow-cytometric quantification of key signaling proteins in the presence or absence of the survival signals APRIL and contact to the stromal cell line ST2, we generated a quantitative model of plasma cell survival. Our model emphasizes the non-redundant nature of the two plasma cell survival signals APRIL and stromal cell contact, and highlights a requirement for differential regulation of individual caspases. The modeling approach allowed us to unify distinct data sets and derive a consistent picture of the intertwined signaling and apoptosis pathways regulating plasma cell survival.

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.

Timescales of IP(3)-evoked Ca(2+) spikes emerge from Ca(2+) puffs only at the cellular level.

The behavior of biological systems is determined by the properties of their component molecules, but the interactions are usually too complex to understand fully how molecular behavior generates cellular behavior. Ca(2+) signaling by inositol trisphosphate receptors (IP(3)R) offers an opportunity to understand this relationship because the cellular behavior is defined largely by Ca(2+)-mediated interactions between IP(3)R. Ca(2+) released by a cluster of IP(3)R (giving a local Ca(2+) puff) diffuses and ignites the behavior of neighboring clusters (to give repetitive global Ca(2+) spikes). We use total internal reflection fluorescence microscopy of two mammalian cell lines to define the temporal relationships between Ca(2+) puffs (interpuff intervals, IPI) and Ca(2+) spikes (interspike intervals) evoked by flash photolysis of caged IP(3). We find that IPI are much shorter than interspike intervals, that puff activity is stochastic with a recovery time that is much shorter than the refractory period of the cell, and that IPI are not periodic. We conclude that Ca(2+) spikes do not arise from oscillatory dynamics of IP(3)R clusters, but that repetitive Ca(2+) spiking with its longer timescales is an emergent property of the dynamics of the whole cluster array.

Mitochondrial energetic metabolism: a simplified model of TCA cycle with ATP production.

Mitochondria play a central role in cellular energetic metabolism. The essential parts of this metabolism are the tricarboxylic acid (TCA) cycle, the respiratory chain and the adenosine triphosphate (ATP) synthesis machinery. Here a simplified model of these three metabolic components with a limited set of differential equations is presented. The existence of a steady state is demonstrated and results of numerical simulations are presented. The relevance of a simple model to represent actual in vivo behavior is discussed.

Modeling Cell-to-Cell Communication Networks Using Response-Time Distributions.

Cell-to-cell communication networks have critical roles in coordinating diverse organismal processes, such as tissue development or immune cell response. However, compared with intracellular signal transduction networks, the function and engineering principles of cell-to-cell communication networks are far less understood. Major complications include: cells are themselves regulated by complex intracellular signaling networks; individual cells are heterogeneous; and output of any one cell can recursively become an additional input signal to other cells. Here, we make use of a framework that treats intracellular signal transduction networks as "black boxes" with characterized input-to-output response relationships. We study simple cell-to-cell communication circuit motifs and find conditions that generate bimodal responses in time, as well as mechanisms for independently controlling synchronization and delay of cell-population responses. We apply our modeling approach to explain otherwise puzzling data on cytokine secretion onset times in T cells. Our approach can be used to predict communication network structure using experimentally accessible input-to-output measurements and without detailed knowledge of intermediate steps.

Diverse drug-resistance mechanisms can emerge from drug-tolerant cancer persister cells.

Cancer therapy has traditionally focused on eliminating fast-growing populations of cells. Yet, an increasing body of evidence suggests that small subpopulations of cancer cells can evade strong selective drug pressure by entering a 'persister' state of negligible growth. This drug-tolerant state has been hypothesized to be part of an initial strategy towards eventual acquisition of bona fide drug-resistance mechanisms. However, the diversity of drug-resistance mechanisms that can expand from a persister bottleneck is unknown. Here we compare persister-derived, erlotinib-resistant colonies that arose from a single, EGFR-addicted lung cancer cell. We find, using a combination of large-scale drug screening and whole-exome sequencing, that our erlotinib-resistant colonies acquired diverse resistance mechanisms, including the most commonly observed clinical resistance mechanisms. Thus, the drug-tolerant persister state does not limit--and may even provide a latent reservoir of cells for--the emergence of heterogeneous drug-resistance mechanisms.

Rhythmic degradation explains and unifies circadian transcriptome and proteome data.

The rich mammalian cellular circadian output affects thousands of genes in many cell types and has been the subject of genome-wide transcriptome and proteome studies. The results have been enigmatic because transcript peak abundances do not always follow the peaks of gene-expression activity in time. We posited that circadian degradation of mRNAs and proteins plays a pivotal role in setting their peak times. To establish guiding principles, we derived a theoretical framework that fully describes the amplitudes and phases of biomolecules with circadian half-lives. We were able to explain the circadian transcriptome and proteome studies with the same unifying theory, including cases in which transcripts or proteins appeared before the onset of increased production rates. Furthermore, we estimate that 30% of the circadian transcripts in mouse liver and Drosophila heads are affected by rhythmic posttranscriptional regulation.

Reliable encoding of stimulus intensities within random sequences of intracellular Ca2+ spikes.

Ca(2+) is a ubiquitous intracellular messenger that regulates diverse cellular activities. Extracellular stimuli often evoke sequences of intracellular Ca(2+) spikes, and spike frequency may encode stimulus intensity. However, the timing of spikes within a cell is random because each interspike interval has a large stochastic component. In human embryonic kidney (HEK) 293 cells and rat primary hepatocytes, we found that the average interspike interval also varied between individual cells. To evaluate how individual cells reliably encoded stimuli when Ca(2+) spikes exhibited such unpredictability, we combined Ca(2+) imaging of single cells with mathematical analyses of the Ca(2+) spikes evoked by receptors that stimulate formation of inositol 1,4,5-trisphosphate (IP3). This analysis revealed that signal-to-noise ratios were improved by slow recovery from feedback inhibition of Ca(2+) spiking operating at the whole-cell level and that they were robust against perturbations of the signaling pathway. Despite variability in the frequency of Ca(2+) spikes between cells, steps in stimulus intensity caused the stochastic period of the interspike interval to change by the same factor in all cells. These fold changes reliably encoded changes in stimulus intensity, and they resulted in an exponential dependence of average interspike interval on stimulation strength. We conclude that Ca(2+) spikes enable reliable signaling in a cell population despite randomness and cell-to-cell variability, because global feedback reduces noise, and changes in stimulus intensity are represented by fold changes in the stochastic period of the interspike interval.