Miguel Fribourg

Cytokine Response Is Determined by Duration of Receptor and Signal Transducers and Activators of Transcription 3 (STAT3) Activation

Background: Interleukin-6 and interleukin-10 both activate the same signaling mediator, STAT3, yet generate nearly opposing responses. Results: Interleukin-6 and interleukin-10 signaling lead to different durations of STAT3 activation and consequently distinct responses. Conclusion: The duration of receptor and STAT3 activation determines the specific cytokine response. Significance: This work reveals a signaling coding mechanism that relieves a cellular information bottleneck. Paradoxically, the pro-inflammatory cytokine IL-6 and the anti-inflammatory cytokine IL-10 both activate STAT3, yet generate nearly opposing cellular responses. Here, we show that the temporal pattern of STAT3 activation codes for the specific cytokine response. A computational model of IL-6 and IL-10 signaling predicted that IL-6 stimulation results in transient activation of STAT3, with a rapid decline in phosphorylation and nuclear localization. In contrast, simulated IL-10 signaling resulted in sustained STAT3 activation. The predicted STAT3 patterns produced by each cytokine were confirmed experimentally in human dendritic cells. Time course microarray studies further showed that the dynamic genome-wide transcriptional responses were nearly identical at early time points following stimulation (when STAT3 is active in response to both IL-6 and IL-10) but divergent at later times (when STAT3 is active only in response to IL-10). Truncating STAT3 activation after IL-10 stimulation caused IL-10 to elicit an IL-6–like transcriptional and secretory response. That the duration of IL-10 receptor and STAT3 activation can direct distinct responses reveals a complex cellular information-coding mechanism that may be relevant to improving the prediction of the effects of drug candidates using this mechanism.

Single-cell analysis shows that paracrine signaling by first responder cells shapes the interferon-β response to viral infection

Signals produced by early responding cells to viral infection calibrate the immune response of the cell population. Taking the time to respond One of the earliest responses of cells to viral infection is expression of the gene encoding interferon-β (Ifnb1). Within a population of cells exposed to virus, the overall immune response is shaped by the percentage of cells that become infected, the proportion of those infected cells that initially express Ifnb1, and the extent of Ifnb1 expression per cell. Patil et al. used single-cell imaging to quantify the amounts of Ifnb1 and viral mRNAs in infected human dendritic cells over time. Mathematical simulations of viral infection indicated that differences in the timing of Ifnb1 induction in individual cells contributed to the dynamics of the population. Indeed, further experiments showed that the earliest responders released paracrine signals that controlled the timing of the antiviral response of the remaining cells within the population. These results have implications for understanding how individual immune cells coordinate the extent of the immune response to viral infection. Immune responses to viral infection are stochastic processes, which initiate in a limited number of cells that then propagate the response. A key component of the response to viral infection entails the synthesis and secretion of type I interferons (IFNs), including the early induction of the gene encoding IFN-β (Ifnb1). With single-cell analysis and mathematical modeling, we investigated the mechanisms underlying how increases in the amount of Ifnb1 mRNA per cell and in the numbers of cells expressing Ifnb1 calibrate the response to viral infection. We used single-cell, single-molecule assays to quantify the early induction of Ifnb1 expression (the Ifnb1 response) in human monocyte-derived dendritic cells infected with Newcastle disease virus, thus retaining the physiological stoichiometry of transcriptional regulators to both alleles of the Ifnb1 gene. We applied computational methods to extract the stochastic features that underlie the cell-to-cell variations in gene expression over time. Integration of simulations and experiments identified the role of paracrine signaling in increasing the number of cells that express Ifnb1 over time and in calibrating the immune response to viral infection.

Interactive Big Data Resource to Elucidate Human Immune Pathways and Diseases

Many functionally important interactions between genes and proteins involved in immunological diseases and processes are unknown. The exponential growth in public high-throughput data offers an opportunity to expand this knowledge. To unlock human-immunology-relevant insight contained in the global biomedical research effort, including all public high-throughput datasets, we performed immunological-pathway-focused Bayesian integration of a comprehensive, heterogeneous compendium comprising 38,088 genome-scale experiments. The distillation of this knowledge into immunological networks of functional relationships between molecular entities (ImmuNet), and tools to mine this resource, are accessible to the public at http://immunet.princeton.edu. The predictive capacity of ImmuNet, established by rigorous statistical validation, is easily accessed by experimentalists to generate data-driven hypotheses. We demonstrate the power of this approach through the identification of unique host-virus interaction responses, and we show how ImmuNet complements genetic studies by predicting disease-associated genes. ImmuNet should be widely beneficial for investigating the mechanisms of the human immune system and immunological diseases.

Human Dendritic Cell Response Signatures Distinguish 1918, Pandemic, and Seasonal H1N1 Influenza Viruses.

ABSTRACT Influenza viruses continue to present global threats to human health. Antigenic drift and shift, genetic reassortment, and cross-species transmission generate new strains with differences in epidemiology and clinical severity. We compared the temporal transcriptional responses of human dendritic cells (DC) to infection with two pandemic (A/Brevig Mission/1/1918, A/California/4/2009) and two seasonal (A/New Caledonia/20/1999, A/Texas/36/1991) H1N1 influenza viruses. Strain-specific response differences included stronger activation of NF-κB following infection with A/New Caledonia/20/1999 and a unique cluster of genes expressed following infection with A/Brevig Mission/1/1918. A common antiviral program showing strain-specific timing was identified in the early DC response and found to correspond with reported transcript changes in blood during symptomatic human influenza virus infection. Comparison of the global responses to the seasonal and pandemic strains showed that a dramatic divergence occurred after 4 h, with only the seasonal strains inducing widespread mRNA loss. IMPORTANCE Continuously evolving influenza viruses present a global threat to human health; however, these host responses display strain-dependent differences that are incompletely understood. Thus, we conducted a detailed comparative study assessing the immune responses of human DC to infection with two pandemic and two seasonal H1N1 influenza strains. We identified in the immune response to viral infection both common and strain-specific features. Among the stain-specific elements were a time shift of the interferon-stimulated gene response, selective induction of NF-κB signaling by one of the seasonal strains, and massive RNA degradation as early as 4 h postinfection by the seasonal, but not the pandemic, viruses. These findings illuminate new aspects of the distinct differences in the immune responses to pandemic and seasonal influenza viruses.

Model of influenza A virus infection: dynamics of viral antagonism and innate immune response.

Viral antagonism of host responses is an essential component of virus pathogenicity. The study of the interplay between immune response and viral antagonism is challenging due to the involvement of many processes acting at multiple time scales. Here we develop an ordinary differential equation model to investigate the early, experimentally measured, responses of human monocyte-derived dendritic cells to infection by two H1N1 influenza A viruses of different clinical outcomes: pandemic A/California/4/2009 and seasonal A/New Caledonia/20/1999. Our results reveal how the strength of virus antagonism, and the time scale over which it acts to thwart the innate immune response, differs significantly between the two viruses, as is made clear by their impact on the temporal behavior of a number of measured genes. The model thus sheds light on the mechanisms that underlie the variability of innate immune responses to different H1N1 viruses.

G Protein-Coupled Receptor Signaling to Kir Channels in Xenopus Oocytes

Kir3 (or GIRK) channels have been known for nearly three decades to be activated by direct interactions with the βγ subunits of heterotrimeric G (Gαβγ) proteins in a membrane-delimited manner. Gα also interacts with GIRK channels and since PTX-sensitive Gα subunits show higher affinity of interaction they confer signaling specificity to G Protein- Coupled Receptors (GPCRs) that normally couple to these G protein subunits. In heterologous systems, overexpression of non PTX-sensitive Gα subunits scavenges the available Gβγ and biases GIRK activation through GPCRs that couple to these Gα subunits. Moreover, all Kir channels rely on their direct interactions with the phospholipid PIP2 to maintain their activity. Thus, signals that activate phospholipase C (e.g. through Gq signaling) to hydrolyze PIP2 result in inhibition of Kir channel activity. In this review, we illustrate with experiments performed in Xenopus oocytes that Kir channels can be used efficiently as reporters of GPCR function through Gi, Gs or Gq signaling. The membrane-delimited nature of this expression system makes it highly efficient for constructing dose-response curves yielding highly reproducible apparent affinities of different ligands for each GPCR tested.

Erythropoietin Receptor-Mediated Molecular Crosstalk Promotes T Cell Immunoregulation and Transplant Survival

Although spontaneous kidney transplant acceptance/tolerance occurs in mice and occasionally in humans, mechanisms remain unclear. Herein we test the hypothesis that EPO, a hormone predominantly produced by the adult kidney, has immunomodulating properties that are required for spontaneous kidney graft acceptance. In vitro, in a manner dependent on the EPO receptor and CD131 on antigen-presenting cells, EPO induced the secretion of active TGFβ by antigen-presenting cells, which in turn converted naïve CD4+ T cells into functional Foxp3+ regulatory T cells (Treg). In murine transplant models, pharmacologic downregulation of kidney-derived EPO prevented spontaneous Treg generation. In a controlled, prospective cohort clinical study, EPO administration at doses used to correct anemia augmented the frequency of peripheral CD4+CD25+CD127lo T cells in humans with CKD. Furthermore, EPO directly inhibited conventional T cell proliferation in vitro via tyrosine phosphatase SHP-1-dependent uncoupling of IL-2Rβ signaling. Conversely, EPO-initiated signals facilitated Treg proliferation by augmenting IL-2Rγ signaling and maintaining constitutively quenched IL-2Rβ signaling. In additional murine transplant models, recombinant EPO administration prolonged heart allograft survival, whereas pharmacologic downregulation of kidney-derived EPO reduced the expression of TGFβ mRNA and abrogated kidney allograft acceptance. Together, our findings delineate the protolerogenic properties of EPO in inhibiting conventional T cells while simultaneously promoting Treg induction, and suggest that manipulating the EPO/EPO receptor signaling axis could be exploited to prevent and/or treat T cell-mediated pathologies, including transplant rejection.

Allospecific Memory B Cell Responses are Dependent on Autophagy

Long‐lived, donor‐reactive memory B cells (Bmems) can produce alloantibodies that mediate transplant injury. Autophagy, an intrinsic mechanism of cell organelle/component recycling, is required for Bmem survival in infectious and model antigen systems, but whether autophagy affects alloreactive Bmem is unknown. We studied mice with an inducible yellow fluorescent protein (YFP) reporter expressed under the activation‐induced cytidine deaminase (AID) promoter active in B cells undergoing germinal center reactions. Up to 12 months after allogeneic sensitization, splenic YFP+ B cells were predominantly IgD–IgM–IgG+ and expressed CD73, CD80, and PD‐L2, consistent with Bmems. Labeled cells contained significantly more cells with autophagosomes and more autophagosomes per cell than unlabeled, naïve B cells. To test for a functional link, we quantified alloantibody formation in mice with B cells conditionally deficient in the requisite autophagy gene ATG7. These experiments revealed absent B cell ATG7 (1) prevented B cell autophagy, (2) inhibited secondary alloantibody responses without altering primary alloantibody formation, and (3) diminished frequencies of alloreactive Bmems. Pharmacological autophagy inhibition with 3‐methyladenine had similar effects on wild‐type mice. Together with new documentation of increased autophagosomes within human Bmems, our data indicate that targeting autophagy has potential for eliminating donor‐reactive Bmems in transplant recipients.

Regulatory architecture of the LβT2 gonadotrope cell underlying the response to gonadotropin-releasing hormone

The LβT2 mouse pituitary cell line has many characteristics of a mature gonadotrope and is a widely used model system for studying the developmental processes and the response to gonadotropin-releasing hormone (GnRH). The global epigenetic landscape, which contributes to cell-specific gene regulatory mechanisms, and the single-cell transcriptome response variation of LβT2 cells have not been previously investigated. Here, we integrate the transcriptome and genome-wide chromatin accessibility state of LβT2 cells during GnRH stimulation. In addition, we examine cell-to-cell variability in the transcriptional response to GnRH using Gel bead-in-Emulsion Drop-seq technology. Analysis of a bulk RNA-seq data set obtained 45 min after exposure to either GnRH or vehicle identified 112 transcripts that were regulated >4-fold by GnRH (FDR < 0.05). The top regulated transcripts constitute, as determined by Bayesian massive public data integration analysis, a human pituitary-relevant coordinated gene program. Chromatin accessibility [assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq)] data sets generated from GnRH-treated LβT2 cells identified more than 58,000 open chromatin regions, some containing notches consistent with bound transcription factor footprints. The study of the most prominent open regions showed that 75% were in transcriptionally active promoters or introns, supporting their involvement in active transcription. Lhb, Cga, and Egr1 showed significantly open chromatin over their promoters. While Fshb was closed over its promoter, several discrete significantly open regions were found at −40 to −90 kb, which may represent novel upstream enhancers. Chromatin accessibility determined by ATAC-seq was associated with high levels of gene expression determined by RNA-seq. We obtained high-quality single-cell Gel bead-in-Emulsion Drop-seq transcriptome data, with an average of >4,000 expressed genes/cell, from 1,992 vehicle- and 1,889 GnRH-treated cells. While the individual cell expression patterns showed high cell-to-cell variation, representing both biological and measurement variation, the average expression patterns correlated well with bulk RNA-seq data. Computational assignment of each cell to its precise cell cycle phase showed that the response to GnRH was unaffected by cell cycle. To our knowledge, this study represents the first genome-wide epigenetic and single-cell transcriptomic characterization of this important gonadotrope model. The data have been deposited publicly and should provide a resource for hypothesis generation and further study.

Elucidation of molecular kinetic schemes from macroscopic traces using system identification

Overall cellular responses to biologically-relevant stimuli are mediated by networks of simpler lower-level processes. Although information about some of these processes can now be obtained by visualizing and recording events at the molecular level, this is still possible only in especially favorable cases. Therefore the development of methods to extract the dynamics and relationships between the different lower-level (microscopic) processes from the overall (macroscopic) response remains a crucial challenge in the understanding of many aspects of physiology. Here we have devised a hybrid computational-analytical method to accomplish this task, the SYStems-based MOLecular kinetic scheme Extractor (SYSMOLE). SYSMOLE utilizes system-identification input-output analysis to obtain a transfer function between the stimulus and the overall cellular response in the Laplace-transformed domain. It then derives a Markov-chain state molecular kinetic scheme uniquely associated with the transfer function by means of a classification procedure and an analytical step that imposes general biological constraints. We first tested SYSMOLE with synthetic data and evaluated its performance in terms of its rate of convergence to the correct molecular kinetic scheme and its robustness to noise. We then examined its performance on real experimental traces by analyzing macroscopic calcium-current traces elicited by membrane depolarization. SYSMOLE derived the correct, previously known molecular kinetic scheme describing the activation and inactivation of the underlying calcium channels and correctly identified the accepted mechanism of action of nifedipine, a calcium-channel blocker clinically used in patients with cardiovascular disease. Finally, we applied SYSMOLE to study the pharmacology of a new class of glutamate antipsychotic drugs and their crosstalk mechanism through a heteromeric complex of G protein-coupled receptors. Our results indicate that our methodology can be successfully applied to accurately derive molecular kinetic schemes from experimental macroscopic traces, and we anticipate that it may be useful in the study of a wide variety of biological systems.