Competitive binding of STATs to receptor phospho-Tyr motifs accounts for altered cytokine responses
Abstract
Cytokines elicit pleiotropic and non-redundant activities despite strong overlap in their usage of receptors, JAKs and STATs molecules. We use IL-6 and IL-27 to ask how two cytokines activating the same signaling pathway have different biological roles. We found that IL-27 induces more sustained STAT1 phosphorylation than IL-6, with the two cytokines inducing comparable levels of STAT3 phosphorylation. Mathematical and statistical modelling of IL-6 and IL-27 signaling identified STAT3 binding to GP130, and STAT1 binding to IL-27Rα, as the main dynamical processes contributing to sustained pSTAT1 levels by IL-27. Mutation of Tyr613 on IL-27Rα decreased IL-27-induced STAT1 phosphorylation by 80% but had limited effect on STAT3 phosphorylation. Strong receptor/STAT coupling by IL-27 initiated a unique gene expression program, which required sustained STAT1 phosphorylation and IRF1 expression and was enriched in classical Interferon Stimulated Genes. Interestingly, the STAT/receptor coupling exhibited by IL-6/IL-27 was altered in patients with systemic lupus erythematosus (SLE). IL-6/IL-27 induced a more potent STAT1 activation in SLE patients than in healthy controls, which correlated with higher STAT1 expression in these patients. Partial inhibition of JAK activation by sub-saturating doses of Tofacitinib specifically lowered the levels of STAT1 activation by IL-6. Our data show that receptor and STATs concentrations critically contribute to shape cytokine responses and generate functional pleiotropy in health and disease.
Data availability
Python (version 3.7) codes for the ABC-SMC model selection and parameter inference can be found in the public repository "https://github.com/PollyJeffrey/Cytokine_modelling", along with the results of the analysis. Phospho-proteomic and proteomic datasets were uploaded to the Proteome Exchange platform with accession numbers PXD024657 and PXD024188 respectively. RNA-seq dataset was uploaded in the GSE database with accession number GSE164479.
Article and author information
Author details
Funding
Horizon 2020 Framework programme (714680)
- Stephan Wilmes
- Jonathan Martinez-Fabregas
- Paul K Fyfe
- Ignacio Moraga Gonzalez
Wellcome Trust (202323/Z/16/Z)
- Stephan Wilmes
- Ignacio Moraga Gonzalez
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Frederik Graw, Heidelberg University, Germany
Ethics
Human subjects: This study was authorized by the French Competent Authority dealing with Research on Human Biological Samples namely the French Ministry of Research. The Authorization number is ECH 19/04. all patients gave their written informed consent
Version history
- Received: December 22, 2020
- Accepted: April 18, 2021
- Accepted Manuscript published: April 19, 2021 (version 1)
- Version of Record published: May 5, 2021 (version 2)
Copyright
© 2021, Wilmes et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,377
- views
-
- 377
- downloads
-
- 21
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
The reconstruction of complete microbial metabolic pathways using ‘omics data from environmental samples remains challenging. Computational pipelines for pathway reconstruction that utilize machine learning methods to predict the presence or absence of KEGG modules in incomplete genomes are lacking. Here, we present MetaPathPredict, a software tool that incorporates machine learning models to predict the presence of complete KEGG modules within bacterial genomic datasets. Using gene annotation data and information from the KEGG module database, MetaPathPredict employs deep learning models to predict the presence of KEGG modules in a genome. MetaPathPredict can be used as a command line tool or as a Python module, and both options are designed to be run locally or on a compute cluster. Benchmarks show that MetaPathPredict makes robust predictions of KEGG module presence within highly incomplete genomes.
-
- Computational and Systems Biology
- Evolutionary Biology
As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.