Combinatorial processing of bacterial and host-derived innate immune stimuli at the single-cell level

Macrophages maintain signaling fidelity in response to ligand mixtures

As immune sentinel cells, macrophages are required to respond specifically to diverse immune threats and initiate appropriate immune responses. This stimulus-response specificity (SRS) is in part encoded in the signaling dynamics of the NFκB transcription factor. While experimental stimulus-response studies have typically focused on single defined ligands, in physiological contexts cells are exposed to multi-ligand mixtures. It remains unclear how macrophages combine multi-ligand information and whether they are able to maintain SRS in such complex exposure conditions. Here, we leveraged an established mathematical model that captures the heterogeneous single-cell NFκB responses of macrophage populations to extend experimental studies with systematic simulations of complex mixtures containing up to five ligands. Live-cell microscopy experiments for some conditions validated model predictions but revealed a discrepancy when TLR3 and TLR9 are stimulated. Refining the model suggested that the observed but unexpected ligand antagonism arises from a limited capacity for endosomal transport which is required for responses to CpG and pIC. With the updated model, we systematically analyzed SRS across all combinatorial-ligand conditions and employed three ways of quantifying SRS involving trajectory decomposition into informative trajectory features or machine learning. Our findings show that macrophages most effectively distinguish single-ligand stimuli, and distinguishability declines as more ligands are combined. However, even in complex combinatorial conditions, macrophages still maintain statistically significant distinguishability. These results indicate a robustness of innate immune response specificity: even in the context of complex exposure conditions, the NFκB temporal signaling code of macrophages can still classify immune threats to direct an appropriate response.

Significance ≤120

Macrophages sense diverse pathogens within complex environments and respond appropriately. Experimental studies have found that the NFκB pathway responds with stimulus-specific dynamics when macrophages are exposed to single ligand stimuli. However, it remains unclear complex contexts might erode this stimulus-specificity. Here we systematically studies NFκB responses using a mathematical model that provides simulations of the heterogeneous population of single cell responses. We show that although the model is parameterized to single ligand data it can predict the responses to multi-ligand mixtures. Indeed, model validation uncovered signaling antagonism between two ligands and the underlying mechanism. Importantly, we found that NFκB signaling dynamics distinguish ligands within multi-ligand mixtures indicating a robustness of the NFκB temporal code that was not previously appreciated.

Quantifying dynamic pro-inflammatory gene expression and heterogeneity in single macrophage cells

Macrophages must respond appropriately to pathogens and other pro-inflammatory stimuli in order to perform their roles in fighting infection. One way in which inflammatory stimuli can vary is in their dynamics-that is, the amplitude and duration of stimulus experienced by the cell. In this study, we performed long-term live cell imaging in a microfluidic device to investigate how the pro-inflammatory genes IRF1, CXCL10, and CXCL9 respond to dynamic interferon-gamma (IFNγ) stimulation. We found that IRF1 responds to low concentration or short duration IFNγ stimulation, whereas CXCL10 and CXCL9 require longer or higherconcentration stimulation to be expressed. We also investigated the heterogeneity in the expression of each gene and found that CXCL10 and CXCL9 have substantial cell-to-cell variability. In particular, the expression of CXCL10 appears to be largely stochastic with a subpopulation of nonresponding cells across all the stimulation conditions tested. We developed both deterministic and stochastic models for the expression of each gene. Our modeling analysis revealed that the heterogeneity in CXCL10 can be attributed to a slow chromatin-opening step that is on a similar timescale to that of adaptation of the upstream signal. In this way, CXCL10 expression in individual cells can remain stochastic in response to each pulse of repeated stimulation, which we also validated by experiments. Together, we conclude that pro-inflammatory genes in the same signaling pathway can respond to dynamic IFNγ stimulus with very different response features and that upstream signal adaptation can contribute to shaping heterogeneous gene expression.

TLR2 is non-redundant in the population and subpopulation responses to Mycobacterium tuberculosis in macrophages and in vivo

Tuberculosis (TB), caused by the pathogenic bacterium Mycobacterium tuberculosis (Mtb), is a global health threat. Targeting host pathways that modulate protective or harmful components of inflammation has been proposed as a therapeutic strategy that could aid sterilization or mitigate TB-associated permanent tissue damage. In purified form, many Mtb components can activate innate immune pathways. However, knowledge of the pathways that contribute most to the observed response to live Mtb is incomplete, limiting the possibility of precise intervention. We took a systematic, unbiased approach to define the pathways that drive the earliest immune response to Mtb. Using a macrophage model of infection, we compared the bulk transcriptional response to infection with the response to a panel of Mtb-derived putative innate immune ligands. We identified two axes of response: an NF-kB-dependent response similarly elicited by all Mtb pathogen-associated molecular patterns (PAMPs) and a type I interferon axis unique to cells infected with live Mtb. Consistent with growing literature data pointing to TLR2 as a dominant Mtb-associated PAMP, the TLR2 ligand PIM6 most closely approximated the NF-kB-dependent response to the intact bacterium. Quantitatively, the macrophage response to Mtb was slower and weaker than the response to purified PIM6. On a subpopulation level, the TLR2-dependent response was heterogeneously induced, with only a subset of infected cells expressing key inflammatory genes known to contribute to the control of infection. Despite potential redundancies in Mtb ligand/innate immune receptor interactions during in vivo infection, loss of the TLR2/PIM6 interaction impacted the cellular composition of both the innate and adaptive compartments. IMPORTANCE Tuberculosis (TB) is a leading cause of death globally. Drug resistance is outpacing new antibiotic discovery, and even after successful treatment, individuals are often left with permanent lung damage from the negative consequences of inflammation. Targeting host inflammatory pathways has been proposed as an approach that could either improve sterilization or improve post-treatment lung health. However, our understanding of the inflammatory pathways triggered by Mycobacterium tuberculosis (Mtb) in infected cells and lungs is incomplete, in part because of the complex array of potential molecular interactions between bacterium and host. Here, we take an unbiased approach to identify the pathways most central to the host response to Mtb. We examine how individual pathways are triggered differently by purified Mtb products or infection with the live bacterium and consider how these pathways inform the emergence of subpopulation responses in cell culture and in infected mice. Understanding how individual interactions and immune pathways contribute to inflammation in TB opens the door to the possibility of developing precise therapeutic interventions.

Insights on the NF-κB System Using Live Cell Imaging: Recent Developments and Future Perspectives

The transcription factor family of nuclear factor kappa B (NF-κB) proteins is widely recognized as a key player in inflammation and the immune responses, where it plays a fundamental role in translating external inflammatory cues into precise transcriptional programs, including the timely expression of a wide variety of cytokines/chemokines. Live cell imaging in single cells showed approximately 15 years ago that the canonical activation of NF-κB upon stimulus is very dynamic, including oscillations of its nuclear localization with a period close to 1.5 hours. This observation has triggered a fruitful interdisciplinary research line that has provided novel insights on the NF-κB system: how its heterogeneous response differs between cell types but also within homogeneous populations; how NF-κB dynamics translate external cues into intracellular signals and how NF-κB dynamics affects gene expression. Here we review the main features of this live cell imaging approach to the study of NF-κB, highlighting the key findings, the existing gaps of knowledge and hinting towards some of the potential future steps of this thriving research field.

Stimulus-specific responses in innate immunity: Multilayered regulatory circuits

Immune sentinel cells initiate immune responses to pathogens and tissue injury and are capable of producing highly stimulus-specific responses. Insight into the mechanisms underlying such specificity has come from the identification of regulatory factors and biochemical pathways, as well as the definition of signaling circuits that enable combinatorial and temporal coding of information. Here, we review the multi-layered molecular mechanisms that underlie stimulus-specific gene expression in macrophages. We categorize components of inflammatory and anti-pathogenic signaling pathways into five layers of regulatory control and discuss unifying mechanisms determining signaling characteristics at each layer. In this context, we review mechanisms that enable combinatorial and temporal encoding of information, identify recurring regulatory motifs and principles, and present strategies for integrating experimental and computational approaches toward the understanding of signaling specificity in innate immunity.

Six distinct NFκB signaling codons convey discrete information to distinguish stimuli and enable appropriate macrophage responses

Macrophages initiate inflammatory responses via the transcription factor NFκB. The temporal pattern of NFκB activity determines which genes are expressed and thus, the type of response that ensues. Here, we examined how information about the stimulus is encoded in the dynamics of NFκB activity. We generated an mVenus-RelA reporter mouse line to enable high-throughput live-cell analysis of primary macrophages responding to host- and pathogen-derived stimuli. An information-theoretic workflow identified six dynamical features-termed signaling codons-that convey stimulus information to the nucleus. In particular, oscillatory trajectories were a hallmark of responses to cytokine but not pathogen-derived stimuli. Single-cell imaging and RNA sequencing of macrophages from a mouse model of Sjögren's syndrome revealed inappropriate responses to stimuli, suggestive of confusion of two NFκB signaling codons. Thus, the dynamics of NFκB signaling classify immune threats through six signaling codons, and signal confusion based on defective codon deployment may underlie the etiology of some inflammatory diseases.

Receptor-Ligand Kinetics Influence the Mechanism of Action of Covalently Linked TLR Ligands

We report a mechanistic study comparing the immune activation of conjugated Toll-like receptor (TLR) agonists and their unlinked mixtures. Herein, we synthesized a set of six linked dual agonists with different ligands, molecular structures, receptor locations, and biophysical characteristics. With these dimers, we ran a series of in vitro cell-based assays, comparing initial and overall NF-κB (nuclear factor kappa-light-chain-enhancer of activated B cells) activation, cytokine expression profiles, as well as time-resolved TNF-α (Tumor Necrosis Factor alpha) expression. We show that initial activation kinetics, ligand specificity, and the dose of the agonist influence the activity of these linked TLR systems. These results can help improve vaccine design by showing how linked TLR agonists can improve their potency with the appropriate selection of key criteria.

Co-stimulation with opposing macrophage polarization cues leads to orthogonal secretion programs in individual cells

Macrophages are innate immune cells that contribute to fighting infections, tissue repair, and maintaining tissue homeostasis. To enable such functional diversity, macrophages resolve potentially conflicting cues in the microenvironment via mechanisms that are unclear. Here, we use single-cell RNA sequencing to explore how individual macrophages respond when co-stimulated with inflammatory stimuli LPS and IFN-γ and the resolving cytokine IL-4. These co-stimulated macrophages display a distinct global transcriptional program. However, variable negative cross-regulation between some LPS + IFN-γ-specific and IL-4-specific genes results in cell-to-cell heterogeneity in transcription. Interestingly, negative cross-regulation leads to mutually exclusive expression of the T-cell-polarizing cytokine genes Il6 and Il12b versus the IL-4-associated factors Arg1 and Chil3 in single co-stimulated macrophages, and single-cell secretion measurements show that these specialized functions are maintained for at least 48 h. This study suggests that increasing functional diversity in the population is one strategy macrophages use to respond to conflicting environmental cues.

Integrative analysis suggests cell type-specific decoding of NF-κB dynamics

The complex signaling dynamics of transcription factors can encode both qualitative and quantitative information about the extracellular environment, which increases the information transfer capacity and potentially supports accurate cellular decision-making. An important question is how these signaling dynamics patterns are translated into functionally appropriate gene regulation programs. To address this question for transcription factors of the nuclear factor κB (NF-κB) family, we profiled the single-cell dynamics of two major NF-κB subunits, RelA and c-Rel, induced by a panel of pathogen-derived stimuli in immune and nonimmune cellular contexts. Diverse NF-κB-activating ligands produced different patterns of RelA and c-Rel signaling dynamic features, such as variations in duration or time-integrated activity. Analysis of nascent transcripts delineated putative direct targets of NF-κB as compared to genes controlled by other transcriptional and posttranscriptional mechanisms and showed that the transcription of more than half of the induced genes was tightly linked to specific dynamic features of NF-κB signaling. Fibroblast and macrophage cell lines shared a cluster of such "NF-κB dynamics-decoding" genes, as well as cell type-specific decoding genes. Dissecting the subunit specificity of dynamics-decoding genes suggested that target genes were most often linked to both RelA and c-Rel or to RelA alone. Thus, our analysis reveals the cell type-specific interpretation of pathogenic information through the signaling dynamics of NF-κB.

The versatile plasmacytoid dendritic cell: Function, heterogeneity, and plasticity

Since their identification as the natural interferon-producing cell two decades ago, plasmacytoid dendritic cells (pDCs) have been attributed diverse functions in the immune response. Their most well characterized function is innate, i.e., their rapid and robust production of type-I interferon (IFN-I) in response to viruses. However, pDCs have also been implicated in antigen presentation, activation of adaptive immune responses and immunoregulation. The mechanisms by which pDCs enact these diverse functions are poorly understood. One central debate is whether these functions are carried out by different pDC subpopulations or by plasticity in the pDC compartment. This chapter summarizes the latest reports regarding pDC function, heterogeneity, cell conversion and environmentally influenced plasticity, as well as the role of pDCs in infection, autoimmunity and cancer.

Techniques for Studying Decoding of Single Cell Dynamics

Cells must be able to interpret signals they encounter and reliably generate an appropriate response. It has long been known that the dynamics of transcription factor and kinase activation can play a crucial role in selecting an individual cell's response. The study of cellular dynamics has expanded dramatically in the last few years, with dynamics being discovered in novel pathways, new insights being revealed about the importance of dynamics, and technological improvements increasing the throughput and capabilities of single cell measurements. In this review, we highlight the important developments in this field, with a focus on the methods used to make new discoveries. We also include a discussion on improvements in methods for engineering and measuring single cell dynamics and responses. Finally, we will briefly highlight some of the many challenges and avenues of research that are still open.