Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Methods & Resources
    • Archive
    • Subjects
  • Collections
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Research Article
Transparent Process
Open Access

Single-cell RNA sequencing of mouse islets exposed to proinflammatory cytokines

View ORCID ProfileJennifer S Stancill  Correspondence email, View ORCID ProfileMoujtaba Y Kasmani, Achia Khatun, Weiguo Cui, View ORCID ProfileJohn A Corbett  Correspondence email
Jennifer S Stancill
1Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jennifer S Stancill
  • For correspondence: jstancill@mcw.edu
Moujtaba Y Kasmani
2Department of Microbiology and Immunology, Medical College of Wisconsin, Milwaukee, WI, USA
3Blood Research Institute, Versiti, Milwaukee, WI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Moujtaba Y Kasmani
Achia Khatun
2Department of Microbiology and Immunology, Medical College of Wisconsin, Milwaukee, WI, USA
3Blood Research Institute, Versiti, Milwaukee, WI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Weiguo Cui
2Department of Microbiology and Immunology, Medical College of Wisconsin, Milwaukee, WI, USA
3Blood Research Institute, Versiti, Milwaukee, WI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John A Corbett
1Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for John A Corbett
  • For correspondence: jcorbett@mcw.edu
Published 21 April 2021. DOI: 10.26508/lsa.202000949
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Abstract

Exposure to proinflammatory cytokines is believed to contribute to pancreatic β-cell damage during diabetes development. Although some cytokine-mediated changes in islet gene expression are known, the heterogeneity of the response is not well-understood. After 6-h treatment with IL-1β and IFN-γ alone or together, mouse islets were subjected to single-cell RNA sequencing. Treatment with both cytokines together led to expression of inducible nitric oxide synthase mRNA (Nos2) and antiviral and immune-associated genes in a subset of β-cells. Interestingly, IL-1β alone activated antiviral genes. Subsets of δ- and α-cells expressed Nos2 and exhibited similar gene expression changes as β-cells, including increased expression of antiviral genes and repression of identity genes. Finally, cytokine responsiveness was inversely correlated with expression of genes encoding heat shock proteins. Our findings show that all islet endocrine cell types respond to cytokines, IL-1β induces the expression of protective genes, and cellular stress gene expression is associated with inhibition of cytokine signaling.

Introduction

Pancreatic islets of Langerhans are highly vascularized micro-organs containing heterogeneous populations of hormone-secreting endocrine cells that are responsible for controlling whole-body glucose homeostasis. The most abundant endocrine cell type is the insulin-secreting β-cell, which makes up about 70% of the rodent islet. Other endocrine cells include glucagon-secreting α-cells and somatostatin-secreting δ-cells, which comprise ∼15% and 5% of the islet, respectively. There are also non-endocrine cells in the islets, including endothelial cells and resident macrophages (Bonner-Weir & Orci, 1982; Setum et al, 1991). Type 1 diabetes is an autoimmune disease caused by selective destruction of β-cells (Gepts, 1965) that leads to insulin deficiency and severe hyperglycemia if not treated by daily insulin injections. Whereas the killing of these cells is primarily mediated by T-cell–dependent mechanisms (Like et al, 1985; Wicker et al, 1986; Bergman & Haskins, 1997), inflammatory cytokines such as IL-1β and IFN-γ produced by macrophages and T lymphocytes (Mandrup-Poulsen et al, 1985; Wright et al, 1988; Wright & Lacy, 1989) are also believed to contribute to β-cell damage and disease development. The role of cytokines in the disease process is based on the ability of these soluble mediators to impair β-cell function and induce β-cell death following ex vivo treatment. For example, IL-1β inhibits oxidative phosphorylation and insulin secretion while also inducing ER stress and DNA damage (Padgett et al, 2013; Eizirik et al, 2020). Nitric oxide, produced in β-cells after the expression of inducible nitric oxide synthase (iNOS) in response to IL-1β exposure, mediates these damaging effects (Southern et al, 1990; Corbett et al, 1991, 1992a, 1993; Welsh et al, 1991; Corbett & McDaniel, 1992). However, whereas IL-1β and nitric oxide are primarily viewed as damaging to the islet, we have recently shown that nitric oxide protects β-cells against DNA damage–induced apoptosis and against viral infection (Oleson et al, 2016; Stafford et al, 2020).

Cytokines have been shown to activate multiple pathways in β-cells. IL-1β stimulates nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) and MAPK signaling (Eizirik et al, 1996; Larsen et al, 1998). NF-κB is required for the expression of inflammatory genes such as Nos2 (Saldeen & Welsh, 1994; Kwon et al, 1995; Flodstrom et al, 1996; Darville & Eizirik, 1998). IFN-γ primarily activates JAK-STAT pathways as well as interferon regulatory factor (IRF)-1 expression (Mamane et al, 1999; Schneider et al, 2014). Interestingly, many inflammatory genes, including Nos2, contain IFN response elements in their promoters (Xie et al, 1993), perhaps explaining why IFN-γ potentiates their expression in β-cells (Heitmeier et al, 1997). IFN-γ is also responsible for the induction of many antiviral response genes that function to limit viral replication (Samuel, 2001).

In addition to Nos2, cytokines are known to induce robust islet gene expression changes, many of which are nitric oxide dependent (Cardozo et al, 2001b; Ylipaasto et al, 2005; Ortis et al, 2010; Hughes et al, 2011; Meares et al, 2013). However, previous studies have been limited by the use of whole islets, which contain several different cell types (Ylipaasto et al, 2005; Meares et al, 2013). Although β-cells make up most of the islet, it may not be accurate to assume that expression changes observed in islets can be attributed solely to changes in β-cells. This limitation can partially be overcome by using FACS-purified β-cells, which has been performed by several groups (Eizirik et al, 1992; Corbett et al, 1992b; Strandell et al, 1995; Cardozo et al, 2001a, 2001b), but this purification still ignores the heterogeneity that exists within the β-cell population (Pipeleers et al, 1994; Bader et al, 2016; Dorrell et al, 2016). In addition, the responses of non–β-cells within the islet to cytokines have largely been ignored, potentially because of the naturally low cell numbers in which these populations exist.

Here, we used single-cell RNA sequencing (scRNA-seq) to address these limitations. After exposing isolated mouse islets for 6 h to IL-1β or IFN-γ alone or in combination, we captured and sequenced transcriptomes from more than 7,700 single cells. Our analysis revealed heterogeneity in the β-cell population, with only 70% of β-cells increasing Nos2 in response to treatment with both cytokines. Cytokine nonresponsive cells were enriched for heat shock proteins and other chaperones, suggesting induction of a stress response. To our surprise, we observed that cytokines increase Nos2 expression in most of the δ-cells and in some of the α-cells, and that the response to cytokines is largely the same among all islet endocrine cell types. These responses include increase in antiviral and other immune response genes as well as repression of endocrine cell identity genes. Importantly, expression of selected genes was confirmed by qRT-PCR to be nitric oxide-independent. Our results suggest that islet responses to cytokines are not unique to β-cells but are similar throughout the entire endocrine population of the islet and are characterized by the stimulation of protective gene expression. Taken together with observations that nitric oxide protects β-cells from apoptosis and viral infection, the studies described here support a model in which the primary functions of cytokines are to protect islet endocrine cells from damage.

Results

scRNA-seq of mouse islets after 6-h cytokine exposure reveals cellular heterogeneity

Although some cytokine-mediated changes in islet gene expression are known, and several of the genes induced are protective, we do not have a complete picture of the cell type intrinsic changes stimulated by cytokines. To address this question, scRNA-seq was performed on untreated mouse islets (Sample 1) or mouse islets treated for 6 h with 10 U/ml IL-1β alone (Sample 2), 150 U/ml IFN-γ alone (Sample 3), or the combination of the two cytokines (Sample 4) (Fig 1A). After 6-h exposure, the levels of nitric oxide produced are low, such that observed gene expression changes are primarily nitric oxide-independent (Corbett et al, 1992a). After quality control, we recovered about 2,000 cells per sample. Cells from all samples were visualized using Uniform Manifold Approximation and Projection (UMAP), an algorithm that unbiasedly grouped the cells into 17 clusters based on similarity of gene expression (Fig 1B and Tables S1 and S2). We assigned endocrine cell identities (β-, α-, δ-, and PP-cells) based on enrichment of the four primary islet hormones (insulin, glucagon, somatostatin, and pancreatic polypeptide, respectively) (Figs 1B and S1). β-cells (Ins1, Ins2) comprised 70% (5,436 cells) of our dataset, α-cells (Gcg) 12% (921 cells), δ-cells (Sst) 8% (651 cells), and PP-cells (Ppy) 2% (159 cells). Our analysis also identified two clusters (6 and 11, made up of β- and non-β endocrine cells, respectively) significantly enriched for geminin (Gmnn), an inhibitor of DNA replication that accumulates during the S phase of mitosis, suggesting that they may represent proliferative cells. Using characteristic gene expression, we identified the cell types of the non-endocrine clusters, which each make up less than 2% of our dataset (Figs 1B and S1): endothelial cells (Cd34), macrophages (Ccr5), mesenchymal cells (Col1a1), ductal cells (Krt17), and acinar cells (Cpa1). Notably, Cluster 7, which was primarily comprised of cells originating from samples treated with IL-1β (Fig S2), was highly enriched for Nos2 (iNOS) mRNA, and is mostly made up of β-cells, based on hormone gene expression (Fig 1C). Indeed, Cluster 7 is localized with the rest of the β-cell clusters in the UMAP projection (Fig 1B).

Table S1 Enriched genes in each cluster shown in Fig 1B.

Table S2 Mean genes and UMIs (Unique Molecular Identifiers) per sample per cluster.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1. Single-cell RNA sequencing of mouse islet cells after 6 h cytokine exposure.

(A) Schematic of experimental design. (B) Uniform Manifold Approximation and Projection plot depicting clusters of cells from all four samples. Clusters are numbered and color-coded according to the legend shown at the right. Cell identity was assigned based on enrichment for genes indicated. (C) Dot plot indicating expression levels and percent of cells expressing Ins1, Ins2, Gcg, Sst, Ppy, and Nos2 mRNA in each of the 17 clusters. (D) Significantly enriched genes in cluster 7 were subjected to functional annotation clustering using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Select enriched categories of genes are shown. (E, F, G, H, I) Selected representative genes from the categories shown in (D). Fold change refers to changes in Cluster 7 compared with all other clusters.

Figure S1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S1. Heat map of top five enriched genes in each cluster.

The top five genes with the greatest positive ln (Fold change) in each cluster compared with all other clusters are shown. Red indicates higher scaled expression, and blue indicates lower scaled expression.

Figure S2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S2. Sample origins of each cluster.

Stacked bar plot showing the percentage of each cluster (shown in Fig 1B) that originates from each sample.

We used the Database for Annotation, Visualization and Integrated Discovery (DAVID) (Dennis et al, 2003) to identify significantly enriched categories among the genes up-regulated in Cluster 7 (Fig 1D). By far the most enriched category was “Immune response,” with an enrichment score of nine. This group includes genes such as guanylate binding protein 5 (Gbp5), β defensin 1 (Defb1), CD40 antigen (Cd40), aconitate decarboxylase 1 (Acod1), Toll-like receptor 2 (Tlr2), an apolipoprotein B editing complex member (Apobec3), and a subunit of the transporter associated with antigen processing (Tap2) (Fig 1E). A number of guanylate binding proteins, a family of GTPases with antiviral activities (Anderson et al, 1999), were also strongly up-regulated in Cluster 7 cells (Fig 1F). Several genes encoding proteins with antiviral roles are also enriched in Cluster 7, including viperin (Rsad2); several IFN-induced genes that help recognize and degrade viral RNA (Ifit1, Ifit3, and Isg20); an IFN-stimulated gene that inhibits viral replication (Isg15); interferon regulatory factor 1 (Irf1), a master regulator of the antiviral defense; and a member of the 2′–5′-oligoadenylate synthetase family (Oasl2), which facilitates RNase L-mediated degradation of viral RNA (Fig 1G) (Schneider et al, 2014; Zhu et al, 2014). As expected, in addition to Nos2, we observed enrichment in several other genes involved in NF-κB signaling, including Icam1, encoding inflammatory intercellular adhesion molecule-1; Sod2, encoding manganese-dependent superoxide dismutase; Nfkbia, encoding an NF-κB inhibitor; and Nfkb1, encoding subunit 1 of NF-κB (Fig 1H). Finally, one of the most highly enriched genes in Cluster 7, Cxcl10, is a proinflammatory chemokine (Cardozo et al, 2003), and is up-regulated along with several other members of the chemokine family (Fig 1I).

Distinct β-cell clusters are seen after 6-h cytokine exposure

Based on our initial clustering analysis (Fig 1B), there appeared to be several distinct groups of β-cells. To further characterize this heterogeneity, we performed a separate clustering analysis using only Ins1- and Ins2-expressing cells in our dataset (Clusters 0, 1, 4–8). Computationally isolating and clustering β-cells revealed 7 clusters (Figs 2A and S3A and Table S3). One (β-cell Cluster 4) appeared to only be present in cells treated with IL-1β, whereas another (β-cell Cluster 6) was only present in cells treated with IFN-γ (Fig 2B), and both clusters had high expression of Nos2 mRNA (Fig 2C). To determine the distinguishing characteristics of these two clusters of iNOS-expressing β-cells, we again used DAVID to identify enriched categories of up-regulated genes. There were several categories enriched in both clusters (Immune response, Guanylate binding, and Herpes simplex infection), but β-cell Cluster 6, the smaller of the two, had a higher correlation to these categories (Fig 2D and E). Closer examination revealed that IL-1β-regulated genes, including Defb1, Icam1, Gbp5, and Cd40 were primarily expressed in β-cell Cluster 4, while IFN-induced genes, including Cxcl10, Oasl2, Ifit3, Rsad2 (viperin), and Ddx58 (Rig-I) were enriched in β-cell Cluster 6 (Fig 2C). Differential expression analysis between these two clusters corroborates these conclusions (Fig S3B and C), as several of the top 10 genes increased in Cluster 4 compared to Cluster 6 are related to IL-1 signaling (Nfkb1 and Icam1), whereas genes increased in Cluster 6 include those regulated by IFNs (Rsad2, Ifit3, and Oasl2). In addition to the Nos2-expressing clusters, we also identified a cluster (β-cell Cluster 2) enriched for ribosomal proteins, chaperones, RNA-binding proteins, and mitochondrial proteins (Fig 2F). Genes enriched in these cells included those encoding subunits of heat shock protein 70 (Hspa1a and Hspa1b), eukaryotic translation initiation factor 1 (Eif1), heat shock protein 40 (Dnajb1), and ribosomal protein L9 (Rpl9) (Fig 2C). Enrichment of genes involved in these stress responses suggests that cells in β-cell Cluster 2 had been stressed before treatment.

Table S3 Enriched genes in each β-cell cluster shown in Fig 2A.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2. β-cell clusters after 6-h cytokine exposure.

(A) Uniform Manifold Approximation and Projection plot depicting clusters of β-cells from all four samples. Clusters are numbered and color-coded according to the legend shown at the right. (B) Uniform Manifold Approximation and Projection plots segregated by sample. (C) Dot plot indicating expression levels and percent of cells expressing selected genes in each of the seven β-cell clusters. (D, E, F) Significantly enriched genes in cluster 4 (D), cluster 6 (E), or cluster 2 (F) were subjected to functional annotation clustering using DAVID. Select enriched categories of genes are shown for each cluster.

Figure S3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S3. Enriched genes in each β-cell cluster.

(A) The top 10 genes with the greatest positive ln (Fold change) in each cluster compared with all other clusters are shown in a heat map. Red indicates higher scaled expression, and blue indicates lower scaled expression. (B, C) Summary of differential expression analysis performed between β-cell Clusters 4 and 6. The top 10 genes (based on the ln fold change) increased in Cluster 4 compared with Cluster 6 (B) and in Cluster 6 compared with Cluster 4 (C) are shown.

Cytokine requirements for differential gene expression in β-cells

To determine how the cytokines individually and together influenced gene expression, we performed differential expression analysis comparing the β-cells in Samples 2, 3, and 4 to those in Sample 1 (Table S4). We used a threshold of 0.5 ln(fold change) and an adjusted P-value of 0.01 to determine if expression of a gene was significantly changed. In response to: (1) IL-1β alone, 231 genes (218 up, 13 down); (2) IFN-γ alone, 65 genes (58 up, 7 down); and (3) both cytokines, 225 genes (199 up, 26 down) were modified (Fig 3A). Of the genes up-regulated in response to both IL-1β and IFN-γ (Sample 4), 92 (46.23%) were IL-1β–dependent (also up-regulated in Sample 2), 40 (20.10%) were IFN-γ–dependent (also up-regulated in Sample 3), 13 (6.53%) were induced by either cytokine (up-regulated in Samples 2, 3, and 4), and 54 (27.14%) required both cytokines for significant induction (they were only up-regulated in Sample 4) (Fig 3B). Examples of up-regulated genes from scRNA-seq data are shown (Fig 3C and D), and selected genes were validated by qRT-PCR using mouse islets (Fig 3E). Nos2 and Sod2 (manganese-dependent superoxide dismutase) were among the IL-1β–dependent genes, as expected, as were Defb1 and Icam1 (Fig 3C–E). Among the IFN-γ–dependent genes we found Irf1, as anticipated, as well as Cxcl10 (Fig 3C–E). The latter gene was significantly induced in response to IFN-γ by scRNA-seq, but, interestingly, was about 13-fold higher in response to both cytokines together (Fig 3C and E), suggesting that IL-1β dramatically augments its expression. Because of this dramatic increase in expression with exposure to both cytokines, the increase in expression by IFN-γ alone did not reach statistical significance by qRT-PCR (Fig 3E). Very few genes were induced in all three cytokine-treated samples, but an antiviral guanylate-binding protein (Gbp2) was among them (Fig 3C). Finally, among the genes requiring both cytokines for significant induction were Gbp5, Rsad2, and Cd40 (Fig 3C–E). Gbp5 was not significantly up-regulated by IL-1β alone via scRNA-seq analysis, but was by qRT-PCR (Fig 3D and E). As expected, these genes are not regulated by nitric oxide, as addition of NOS inhibitor NG-monomethyl-ʟ-arginine (NMMA) does not prevent the up-regulation by cytokines (Fig 3E). We were surprised to observe induction of antiviral genes by either IL-1β (Defb1, Gbp2, and Gbp5) or IFN-γ (Irf1 and Gbp2), as antiviral function is typically attributed to IFN signaling, but not to IL-1 signaling (Samuel, 2001). We did not observe induction of pro-apoptotic genes in any of our samples, as has been suggested (Fig S4) (Grunnet et al, 2009). Together, these results suggest that (1) in β-cells, IL-1β has a greater influence on gene expression than IFN-γ, but a number of genes require both cytokines to be significantly induced; and (2) IL-1β and IFN-γ both individually up-regulate antiviral genes in β-cells.

Table S4 Differential expression analysis of β-cells in response to each treatment.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3. Cytokine requirements for different genes in β-cells.

(A) Stacked bar plot showing the number of significantly up-regulated or down-regulated genes in β-cells in S2 (IL-1β), S3 (IFN-γ), or S4 (IL-1β + IFN-γ) compared with S1 (Untreated). ln (Fold change) > 0.5 and an adjusted P-value < 0.01 were used as cutoffs to determine significance. (B) Genes significantly up-regulated in β-cells treated with IL-1β + IFN-γ (S4) were compared with genes significantly up-regulated in β-cells treated with IL-1β (S2) or with IFN-γ (S3). Pie chart shows the percentage of up-regulated genes in S4 that are IL-1β–dependent (also up-regulated in S2), IFN-γ–dependent (also up-regulated in S3), IL-1β- or IFN-γ–dependent (up-regulated in S2, S3, and S4), or IL-1β– and IFN-γ–dependent (only up-regulated in S4). (C) ln (fold change) from single-cell RNA sequencing data of selected representative genes in the categories shown in (B). Green bars: IL-1β–dependent. Teal bars: IFN-γ–dependent. Yellow bars: IL-1β– or IFN-γ–dependent. Purple bars: IL-1β and IFN-γ–dependent. (D) Violin plots showing the expression level of selected genes in β-cells from each of the four samples. (E) Expression of selected genes (compared with Gapdh levels) determined by qRT-PCR using mouse islets treated with the indicated treatments for 6 h. Concentrations are as follows: 10 U/ml IL-1β, 150 U/ml IFN-γ, and 2 mM NMMA. Error bars represent SEM. qRT-PCR results are from three independent experiments with statistical significance indicated. *P < 0.05 (versus untreated).

Figure S4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S4. Expression of pro-apoptotic genes in β-cells.

Violin plots showing the expression levels of selected pro-apoptotic genes in β-cells in each of the four samples.

Nos2 mRNA is increased in a subset of δ- and α-cells

Our initial clustering analysis revealed that most “cytokine-responsive” cells (those that express Nos2 mRNA in response to IL-1β and IFN-γ exposure) were found in Cluster 7 of our overall dataset, a β-cell–enriched cluster (Fig 1C). Indeed, comparison of the expression of Nos2 and Ins1 in Sample 4 showed enhanced accumulation of Nos2 mRNA in 71% of β-cells (Fig 4A and B). However, closer examination of the α- and δ-cell clusters (Clusters 2 and 3, respectively) revealed that a subset of these non–β-cells also express Nos2 mRNA in response to treatment with both cytokines (Fig 4C and E). 69% of δ-cells and 39% of α-cells express Nos2 in Sample 4, suggesting that subsets of these cell types are “cytokine responsive” (Fig 4D and F and Tables S5 and S6). Nos2 mRNA was up-regulated in β-cells and δ-cells in response to IL-1β alone or in combination with IFN-γ (Fig 4G and H), whereas α-cells required both cytokines (Fig 4I). In agreement, mouse αTC1 cells increase Nos2 mRNA in response to 6 h exposure to 10 U/ml IL-1β and 150 U/ml IFN-γ (Fig S5A). Fifty percent of PP-cells expressed Nos2 in response to both cytokines, but this increase was not statistically significant (Fig S6C and Table S7). These findings are unexpected given that β-cells are thought to be the only islet endocrine cell type capable of generating nitric oxide upon IL-1β stimulation (Corbett et al, 1992b; Corbett & McDaniel, 1995; Arnush et al, 1998).

Table S5 Differential expression analysis of δ-cells in response to each treatment.

Table S6 Differential expression analysis of α-cells in response to each treatment.

Table S7 Differential expression analysis of PP-cells in response to each treatment.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4. Populations of α- and δ-cells express Nos2 mRNA in response to IL-1β.

(A, C, E) Uniform Manifold Approximation and Projection plots showing the overlap in expression of Nos2, shown in green, and Ins1 (A), Sst (C), or Gcg (E) shown in red. Intermediary colors denote cells that co-express Nos2 with each respective hormone. Expression levels of each gene are indicated to the right of each panel. (B, D, F) Pie charts showing the percentage of cytokine-responsive (Nos2-expressing) and nonresponsive (Nos2 non-expressing) β-cells (B), δ-cells (D), or α-cells (F) in Sample 4. (G, H, I) Violin plots showing the expression levels of Nos2 mRNA in β-cells (G), δ-cells (H), or α-cells (I) in each sample.

Figure S5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S5. Expression of selected genes in αTC1 cells in response to cytokines.

(A, B, C, D, E) Relative expression of Nos2 (A), Ptgs2 (B), Icam1 (C), Cxcl10 (D), and Gbp5 (E) (all compared with Gapdh levels) determined by qRT-PCR using αTC1 Clone 6 cells treated with the indicated treatments for 6 h. Concentrations are as follows: 10 U/ml IL-1β, 150 U/ml IFN-γ, and 2 mM NMMA. Error bars represent SEM. Results are from 3 to 4 independent experiments with statistical significance indicated. *P < 0.05 (versus untreated).

Figure S6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S6. Cytokine-responsive genes in PP-cells.

(A) Genes significantly up-regulated in PP-cells treated with IL-1β + IFN-γ (S4) were compared with genes significantly up-regulated in response to IL-1β (S2) or to IFN-γ (S3). Pie chart shows the percentage of up-regulated genes in S4 PP-cells that are IL-1β–dependent (also up-regulated in S2), IFN-γ–dependent (also up-regulated in S3), or IL-1β– and IFN-γ–dependent (only up-regulated in S4). (B) ln (fold change) of selected representative genes in the categories shown in (A) in S4 PP-cells. Green bars: IL-1β–dependent. Teal bars: IFN-γ–dependent. Purple bars: IL-1β– and IFN-γ–dependent. (C) Violin plots showing the expression levels of selected genes in PP-cells from each of the four samples.

IL-1β and IFN-γ induce many of the same genes in β-, δ-, and α-cells

Because β-, δ-, and α-cells all respond to cytokine stimulation by expressing iNOS, we hypothesized that they may display similar changes in gene expression in addition to Nos2. To test this, we compared the expression of the top up-regulated genes in β-cells to their expression levels in δ- and α-cells. These genes were significantly induced in all three cell types, with comparable fold changes (Fig 5A and Tables S4–S6). In addition to fold change, we also compared the percentage of each cell type with detectable transcripts for each gene. For most genes, fewer than 5% of β-cells expressed the transcript in untreated cells, but this increased to at least 70% in response to both cytokines (Fig 5B). For a few genes, however, including Mt2 and Sod2, a similar percentage of cells expressed the transcript in both untreated and cytokine-treated samples (Fig 5B). Importantly, the patterns of gene expression were comparable in both δ- and α-cells (Fig 5C and D), suggesting that β-cells and non–β endocrine cells respond to initial cytokine exposure in a similar manner.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5. Top up-regulated genes in endocrine cells in response to IL-1β and IFN-γ.

(A) Histogram showing ln (fold change) of selected genes in β-cells (turquoise bars), δ-cells (magenta bars), and α-cells (yellow bars) from Sample 4 compared with Sample 1. The top 11 genes induced in β-cells are shown. (B, C, D) Histograms showing the percent of β-cells (B), δ-cells (C), and α-cells (D) in the untreated sample (hatched bars) and in the sample treated with IL-1β + IFN-γ (solid bars) expressing transcripts of selected genes shown in (A).

Cytokine-responsive genes in non–β-endocrine cells

Like our analysis of the β-cells (Fig 3), we performed differential expression analysis comparing the δ-, α-, and PP-cells in Samples 2, 3, and 4 to those in Sample 1 (Tables S5–S7). In contrast to the β-cells, in which 46% of genes were induced by IL-1β, (Fig 3A), only 25% and 32% of genes were induced by IL-1β in δ- and α-cells, respectively (Fig 6A and C). Accordingly, δ- and α-cells required stimulation by both cytokines to induce a greater percentage of cytokine-dependent genes than did β-cells (Fig 6A and C), suggesting that β-cells may be more sensitive to IL-1β stimulation than other islet endocrine cells. In support of this conclusion, 35% and 38% of the genes that were significantly induced by IL-1β alone in β-cells required both IL-1β and IFN-γ for induction in δ- and α-cells, respectively. Defb1, Sod2, and Icam1 are induced in response to IL-1β, whereas Irf1 is dependent on IFN-γ in all three cell types (Figs 3C–E and 6B and D–F). There are some genes, such as Cxcl10, which are induced by IFN-γ alone in β-cells but required stimulation by both cytokines in δ- and α-cells (Figs 3C–E and 6B and D–F). In addition, Cd40 and Rsad2 are induced by both cytokines in β- and δ-cells but are not significantly induced in α-cells (Figs 3C–E and 6B and D–F). In accordance with our scRNA-seq data, αTC1 cells treated with IL-1β and IFN-γ for 6 h have increased mRNA accumulation of Icam1, Cxcl10, and Gbp5 in a nitric oxide–independent manner (Fig S5C–E). Although not detected by scRNA-seq, Ptgs2 mRNA levels (encoding cyclooxygenase 2 or Cox2) were increased by 6-h IL-1β exposure in αTC1 cells, further suggesting the induction of NF-κB driven gene expression in non-β endocrine cells in response to cytokines (Fig S5B). In general, the expression of fewer genes was increased in PP-cells in response to IL-1β when compared to the other islet endocrine cells, but Cxcl10, Icam1, Defb1, and Gbp5 expression levels were increased after stimulation by both cytokines (Fig S6A–C).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6. Cytokine-responsive genes in δ- and α-cells.

(A, C) Genes significantly up-regulated in δ- and α-cells treated with IL-1β + IFN-γ (S4) were compared with genes significantly up-regulated in δ- and α-cells with IL-1β (S2) or with IFN-γ (S3). Pie charts show the percentage of up-regulated genes in S4 δ-cells (A) and α-cells (C) that are IL-1β–dependent (also up-regulated in S2), IFN-γ–dependent (also up-regulated in S3), IL-1β– or IFN-γ–dependent (up-regulated in S2, S3, and S4), or IL-1β– and IFN-γ–dependent (only up-regulated in S4). (B, D) ln (fold change) of selected representative genes in the categories shown in (A) and (C) in S4 δ-cells (B) and S4 α-cells (D). Green bars: IL-1β–dependent. Teal bars: IFN-γ–dependent. Yellow bars: IL-1β– or IFN-γ–dependent. Purple bars: IL-1β– and IFN-γ–dependent. (E, F) Violin plots showing the expression levels of selected genes in δ-cells (E) and α-cells (F) from each of the four samples.

These data suggest that although many of the same genes are induced in all islet endocrine cells in response to cytokines, β-cells may be more sensitive to stimulation by individual cytokines than are δ-, α-, and PP-cells.

Repression of endocrine cell identity genes in an IL-1β–dependent manner

Others have reported repression of β-cell identity factors, including Pdx1, Mafa, and Slc2a2, in response to prolonged IL-1β treatment (of at least 24 h) (Cardozo et al, 2001a; Nordmann et al, 2017). Therefore, we were surprised to observe repression of these and other β-cell identity factors (Ucn3, Nkx6-1, and Hadh) after only 6 h of treatment (Fig 7A and E). Repression is IL-1β–dependent, as it was observed only in Samples 2 (IL-1β alone) and 4 (IL-1β + IFN-γ), but not in Sample 3 (IFN-γ alone) (Fig 7A and E). Pdx1, Mafa, and Nkx6-1 are transcription factors that maintain β-cell identity and promote insulin transcription and secretion (Matsuoka et al, 2004; Henseleit et al, 2005; Gao et al, 2014); Ucn3 encodes a maturation factor (Blum et al, 2012); Slc2a2 encodes the glucose transporter Glut2 that has a glucose sensing function in β-cells (Guillam et al, 1997); and Hadh encodes an enzyme that participates in β-oxidation of fatty acids (Hardy et al, 2007). Repression of identity factors did not coincide with increased expression of β-cell disallowed genes, such as Ldha (lactate dehydrogenase A), Slc16a1 (monocarboxylate transporter 1), and Hk1 (hexokinase) (Fig S7A) nor with genes associated with dedifferentiation, such as Gast (gastrin), Neurog3, and Aldh1a3 (Fig S7B) (Talchai et al, 2012; Kim-Muller et al, 2016; Dahan et al, 2017). Furthermore, we observed repression of δ-cell identity factor Hhex, which promotes δ-cell specification and somatostatin transcription, in an IL-1β–dependent manner (Zhang et al, 2014) (Fig 7B and E). Repression of α-cell identity factor Arx, which promotes α-cell specification and maturity (Wilcox et al, 2013), did not achieve statistical significance, but is observed when evaluated by both scRNA-seq and qRT-PCR (Fig 7C and E). Not only was the expression level of each factor reduced, but the percentage of cells with detectable levels of the gene was also dramatically reduced by cytokines (Fig 7D). Importantly, repression of these factors was not affected by addition of NMMA (Fig 7E). These data suggest that IL-1β represses not only β-cell identity factors, but also δ- and α-cell identity factors in a nitric oxide–independent manner.

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7. IL-1β represses endocrine cell identity genes.

(A, B, C) Violin plots showing the expression level of selected identity genes in β-cells (A), δ-cells (B), and α-cells (C) from each of the four samples. (D) Histogram showing the percent of untreated cells and cells treated with IL-1β + IFN-γ expressing detectable transcripts of selected genes. (E) Expression of selected genes (compared with Gapdh levels) determined by qRT-PCR using mouse islets treated with the indicated treatments for 6 h. Concentrations are as follows: 10 U/ml IL-1β, 150 U/ml IFN-γ, and 2 mM NMMA. Error bars represent SEM. qRT-PCR results are from three independent experiments with statistical significance indicated. *P < 0.05, ‡P = 0.05 (versus untreated).

Figure S7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S7. Expression of genes associated with dedifferentiation in β-cells.

(A, B) Violin plots showing the expression level of disallowed genes (A) and of dedifferentiation genes (B) in β-cells in each of the four samples.

Islet stress is negatively correlated with “cytokine responsiveness”

Whereas most β-cells treated with IL-1β and IFN-γ (Sample 4) expressed Nos2 (71%), nearly one third (29%) failed to increase Nos2 expression. We and others have shown that islet stress such as heat shock or ER stress attenuates the ability of cytokines to stimulate iNOS expression in β-cells (Bellmann et al, 1995; Scarim et al, 1998; Burkart et al, 2000; Weber et al, 2004b). Islet stress, as identified by increased expression of heat shock protein (Hsp)70, is associated with inhibition of signaling such that IL-1–induced NF-κB activation and MAPK signaling, and IFN-γ–stimulated STAT phosphorylation, are attenuated (Bellmann et al, 1995, Scarim et al, 1998; Weber et al, 2003, 2004a, 2004b). Indeed, comparison of the expression patterns of Nos2 and Hspa1a, encoding a subunit of Hsp70, suggested a negative correlation between the two genes in β-cell clusters: cells expressing high levels of Hspa1a do not express Nos2 in response to IL-1β and IFN-γ (Figs 8A and B and S8A). These cytokine nonresponsive cells appeared to correlate by UMAP visualization to β-cell Cluster 2, which was enriched in genes encoding ribosomal proteins and chaperones (Fig 2F). Indeed, differential expression analysis using Nos2-expressing versus non-expressing β-cells in Sample 4 revealed significant enrichment of genes encoding heat shock proteins including Hspa1a, Hspa1b, Dnajc24, and Dnajc3 (Fig 8E and Table S8).

Table S8 Differential expression analysis of cytokine-responsive (expressing Nos2 in response to IL-1β + IFN-γ) versus nonresponsive (not expressing Nos2 in response to IL-1β + IFN-γ) β-cells.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8. Heat shock protein gene expression is negatively correlated with Nos2 expression and identity gene repression.

(A, B, C, D) Uniform Manifold Approximation and Projection plots segregated by sample and colored to indicate expression level of Hspa1a (A), Nos2 (B), Gbp5 (C), and Hadh (D) in β-cell clusters. Arrowhead points out the population of cells with high expression of Hspa1a that fails to induce Nos2 and Gbp5 or to repress Hadh in response to IL-1β and IFN-γ. Panel (A) is repeated in Fig S8A. (E, F) Differential expression analysis comparing cytokine-responsive (Nos2-expressing) and nonresponsive (Nos2 non-expressing) β-cells treated with IL-1β and IFN-γ. Histograms show the expression level (ln (fold change)) of selected genes that are significantly up-regulated (E) or down-regulated (F) in cytokine nonresponsive β-cells.

Figure S8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure S8. Additional genes with negative correlations with Hspa1a expression.

(A, B, C, D, E) Uniform Manifold Approximation and Projection plots segregated by sample and colored to indicate expression level of Hspa1a (A), Defb1 (B), Icam1 (C) Cxcl10 (D), and Mafa (A) in β-cell clusters. Arrowhead points out the population of cells with high expression of Hspa1a in Sample 4. Panel (A) is repeated in Fig 8A.

In addition to the negative correlation between Nos2 and heat shock gene expression, cells expressing high levels of Hspa1a also did not express antiviral genes, including Gbp5, Gbp2, Gbp4, Defb1, Rsad2, and Irf1 in response to cytokine treatment (Figs 8C and F and S8A and B). Other cytokine-induced genes follow a similar pattern of expression in cells expressing high levels of Hspa1a, including Icam1 and Cxcl10 (Fig S8A, C, and D). Finally, β-cells with no or low expression of Hspa1a had reduced expression of identity gene Hadh in response to IL-1β treatment, whereas cells with high expression of Hspa1a maintained Hadh expression (Fig 8A and D). A similar pattern of expression was observed between Hspa1a and Mafa (Fig S8A and E). In accordance, differential expression analysis revealed that cells failing to increase Nos2 had significantly higher expression of identity genes, including Mafa, Hadh, Ucn3, Slc2a2, and Pdx1 in response to cytokine treatment (Fig 8E). These observations suggest that activation of the heat shock response (or stress response) not only prevents cytokine-induced Nos2 expression, but also prevents other aspects of cytokine signaling, including induction of antiviral genes and identity gene repression.

Discussion

Autoimmune, or type 1 diabetes, is characterized by selective destruction of insulin-producing pancreatic β-cells that is preceded by an inflammatory reaction in and around the islets of Langerhans (Gepts, 1965; Like et al, 1985; Wicker et al, 1986; Bergman & Haskins, 1997). T lymphocytes and macrophages that infiltrate the islet participate in this inflammatory reaction, releasing cytokines such as IFN-γ and IL-1β (Wright et al, 1988; Wright & Lacy, 1989). IL-1β is known to inhibit β-cell function and cause islet damage in a nitric oxide–dependent manner (Mandrup-Poulsen et al, 1985; Bendtzen et al, 1986). Nitric oxide is produced by β-cells after IL-1–induced iNOS expression and attenuates insulin secretion by inhibiting oxidative metabolism and decreasing ATP levels in β-cells (Corbett et al, 1993; Corbett & McDaniel, 1995).

Over the years, studies have partially elucidated how cytokines and nitric oxide affect islet or β-cell gene expression (Cardozo et al, 2001a, 2001b; Hughes et al, 2011; Meares et al, 2013). However, these have been limited by the use of intact islets, which represents a heterogeneous cell population, or purified β-cells, which ignores the potential influences of population heterogeneity. In addition, few studies have explored the effects of cytokines on islet non-β-cells. Finally, many previous studies were performed after prolonged treatment periods (24 h or greater), which miss immediate or early effects of cytokines and may be confounded by secondary and tertiary changes due to the actions of the cytokines. Here, we chose a short treatment period of 6 h, or an early treatment, that allowed us to determine how islet cells initially respond to cytokine stimulation, and diminish the effects of secondary responses, such as nitric oxide, that may modify gene expression profiles (Corbett et al, 1992a).

The initial clustering analysis of our scRNA-seq data revealed 17 distinct clusters, delineating several different cell types, most of which were β-cells (Figs 1 and S1). The primary Nos2-expressing cluster, Cluster 7, was enriched for genes involved in the immune response, guanylate binding, antiviral defense, NF-κB signaling, and chemokine signaling, and primarily consisted of β-cells (Fig 1). When we computationally isolated and re-clustered only the β-cells (Clusters 0, 1, 4–8), two clusters were enriched for Nos2 (Figs 2 and S3). One (β-cell Cluster 4) appeared to be primarily enriched for IL-1β–induced genes, whereas the other (β-cell Cluster 6) was enriched for IFN-γ–induced genes. Intriguingly, we also identified a cluster of β-cells with high expression of ribosomal proteins, heat shock proteins, and other chaperones. This cluster persisted across the samples, suggesting that these cells were stressed before cytokine exposure.

Although not obvious in our initial clustering analysis, we found that subpopulations of δ- and α-cells express Nos2 mRNA (Fig 4) in response to cytokines. Like β-cells, δ-cells increased Nos2 in response to IL-1β alone, but α-cells required both IL-1β and IFN-γ to significantly increase Nos2. These findings were unexpected, given that previous studies using rodent and human islets have suggested that β-cells are the primary endocrine cellular source of iNOS-derived nitric oxide in the islet (Corbett et al, 1992b; Corbett & McDaniel, 1995; Arnush et al, 1998). There are many potential reasons for this discrepancy that include species-dependent effects (rat or human versus mouse) and differences in cytokines used (IL-1β alone versus IL-1β and IFN-γ). However, the most likely reason for the differences is that most colocalization studies did not specifically focus on iNOS expression in α- or δ-cells but focused on iNOS expression in insulin-positive cells. Thus, iNOS expression in non–β-endocrine cells was likely missed or attributed to non-endocrine cells (macrophages or endothelial cells) known to express iNOS and to be present in islets. With these caveats in mind, our studies indicate that all endocrine cell types in islets are responsive to cytokines, including enhanced expression of iNOS. Whether or not iNOS induction results in nitric oxide production in these non–β-cells, and how nitric oxide affects δ- and α-cell function and survival in the context of cytokine exposure has yet to be determined.

The similarities in response to cytokines among β-, δ-, and α-cells do not end with Nos2 induction. Our analysis revealed that the top up-regulated genes in β-cells in response to IL-1β and IFN-γ (Sample 4) are similarly induced in δ- and α-cells, both in fold change and in percentage of cells expressing the transcript (Fig 5). In fact, of the top 100 genes induced in Sample 4 β-cells, 61 are also induced in δ-cells, and 66 are induced in α-cells. Similarities in the expression of genes that were repressed by cytokines were also observed. IL-1β has previously been shown to repress genes critical for maintaining β-cell function and identity, including Pdx1, Mafa, and Slc2a2 (Cardozo et al, 2001a; Nordmann et al, 2017). However, these studies observe this repression after prolonged cytokine exposure of 24 h or longer (Cardozo et al, 2001a; Nordmann et al, 2017). We show a significant repression of these same genes after only 6 h of IL-1β treatment, suggesting that this is an immediate response to cytokine exposure rather than an indication of dedifferentiation after prolonged exposure (Fig 7). Furthermore, we found repression of genes that maintain δ- and α-cell identity, such as the transcription factors Hhex and Arx, in the δ- and α-cell populations, respectively (Fig 7). Like the repression of β-cell identity genes, repression of Hhex and Arx only occurs upon exposure to IL-1β, but not to IFN-γ alone. These exciting results suggest that: (1) IL-1β is a primary mediator of the decrease in expression of genes associated with endocrine cell identity; (2) this action is not modified by IFN-γ; and (3) the loss of identity is not unique to β-cells but occurs in all islet endocrine cells.

Although IFN-γ is a known activator of antiviral responses (Samuel, 2001), our analysis revealed that IL-1β induces several protective genes in islet endocrine cells, a role for this cytokine that is not well-characterized. Defb1 (β-defensin 1) is an antimicrobial peptide (Semple & Dorin, 2012) and is robustly increased in β-, δ-, and α-cells in a IL-1β–dependent manner (Figs 3 and 6). Expression of this gene was confirmed by qRT-PCR analysis and is nitric oxide independent (Fig 3E). Guanylate-binding proteins are known to be IFN-regulated and have broad antiviral and antimicrobial activities (Anderson et al, 1999), but here, we found Gbp2 to be regulated by either IL-1β or IFN-γ alone (Figs 3 and 6). In addition, Gbp5 was significantly induced by IL-1β in islets by qRT-PCR analysis and was trending to be induced by IL-1β by scRNA-seq analysis in β-cells (Fig 3), but not in δ- or α-cells (Fig 6). Recently, both Gbp2 and Gbp5 have been shown to inhibit viral glycoprotein maturation (Braun et al, 2019). These findings demonstrate that both IL-1β and IFN-γ individually induce antiviral genes in islet endocrine cells.

There is a difference between β-cells and non–β-endocrine cells when comparing the changes in gene expression in response to the individual cytokines. In β-cells, IL-1β had the greatest effect on gene expression, significantly increasing more than 200 genes, whereas IFN-γ significantly increased about 65 genes (Fig 3). This trend held true when assessing the genes induced in cells treated with both cytokines, in which about 46% were IL-1β–dependent and 20% were IFN-γ–dependent, whereas 27% were only induced under conditions in which both cytokines were present (Fig 3). δ-, α-, and PP-cells had a greater percentage of genes (53%, 48%, and 69%, respectively) that required the presence of both cytokines for significant up-regulation (Figs 6 and S6). The difference in cytokine requirement may suggest that endocrine non–β-cells are less sensitive to the effects of either cytokine alone, or may be due to limitations in our experimental approach, as the α-, δ-, and PP-cell populations were at least five times smaller than the β-cell population captured (5,436 cells versus 921, 651, and 159 cells). Because of lower cell numbers, statistical significance was more difficult to achieve in α-, δ-, and PP-cells. Independent of the individual cytokine requirements, our results show that all islet endocrine cells respond to cytokine exposure in a similar manner: by rapidly up-regulating iNOS, antiviral genes, and additional genes involved in the immune response, and by repressing transcription factors and genes critical for maintenance of cell identity.

Our analyses suggest that all islet endocrine cells respond to cytokines in a similar manner, including up-regulation of Nos2 and antiviral pathways and repression of islet identity genes (Fig 9). Generally, cytokines and nitric oxide are viewed as damaging to β-cells based on observations that exogenous IL-1β treatment impairs oxidative phosphorylation and insulin secretion and eventually causes β-cell death with prolonged exposure in a nitric oxide–dependent manner (Mandrup-Poulsen et al, 1985; Bendtzen et al, 1986; Southern et al, 1990; Welsh et al, 1991; Corbett & McDaniel, 1992, 1995; Corbett et al, 1992a, 1993; Cunningham & Green, 1994; Mandrup-Poulsen, 1996). However, our recent studies suggest that the ability of nitric oxide to inhibit mitochondrial oxidation and ATP generation protects β-cells from apoptosis (Oleson et al, 2016) and viral infection (Stafford et al, 2020). Our observations here that IL-1β induces antiviral genes (Fig 3) and fails to induce apoptotic mediators (Fig S4) dovetail with these recent findings and further support the idea that the primary functions of IL-1β and nitric oxide are to protect β-cells from damage.

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9. Model.

Upon stimulation by IL-1β and IFN-γ, islet endocrine cells (β-, α-, and δ-cells) respond by up-regulating inducible nitric oxide synthase (Nos2) and a subset of antiviral genes, while simultaneously down-regulating identity factors, such as Pdx1, Mafa, Hhex, and Arx. However, induction of cellular stress, such as heat shock or ER stress, attenuates these aspects of cytokine signaling.

Finally, we observed that cytokine responsiveness in β-, α-, and δ-cells was heterogeneous. For example, Nos2 was induced in 71% of β-cells, 69% of δ-cells, and 39% of α-cells in response to IL-1β and IFN-γ (Fig 4). The population of cells that failed to express Nos2 under these conditions corresponded with β-cell Cluster 2. The identifying feature of this cluster was elevated expression of genes encoding heat shock proteins and chaperones (Figs 2 and 8). Indeed, we observed that cells with high Hspa1a expression failed to express Nos2 (Fig 8), or in other words, the expression of heat shock genes is enhanced in endocrine cells that fail to respond to cytokines. This finding is in agreement with previous studies demonstrating that in response to heat shock (Scarim et al, 1998) or ER stress (Weber et al, 2004a) cytokines fail to stimulate iNOS expression in rodent and human islets (Steer et al, 2006). Furthermore, this stress response is typical of human islet preparations (Welsh et al, 1995) and is likely responsible for the conclusions that nitric oxide does not mediate the damaging effects of cytokines on human islets (Eizirik et al, 1994; Burkart et al, 2000). More interesting are the findings that cells with high expression of Hspa1a also maintained expression of identity genes and did not express antiviral genes after cytokine treatment, suggesting a negative correlation between heat shock and cytokine signaling more broadly than just iNOS expression (Figs 8 and S8). It is likely that the stress signature observed in this cluster is an artifact of the islet isolation process and does not represent an endogenous cell population. Together with the observation that cytokine-induced nitric oxide promotes heat shock protein expression in rat islets (Helqvist et al, 1991), our results suggest that the heat shock response may provide an internal “off switch” to prevent cytokine-mediated islet damage after protective signaling has been initiated. In other words, nitric oxide may not only protect β-cells through inhibition of mitochondrial oxidation (Oleson et al, 2016; Stafford et al, 2020), but may also provide a means to shut down prolonged cytokine signaling, which is known to be damaging.

Materials and Methods

Materials, cells, and animals

Male C57BL6/J mice were purchased from Jackson Laboratories and housed in the MCW Biomedical Resource Center. All animal use and experimental procedures were approved by the Institutional Animal Care and Use Committees at the Medical College of Wisconsin. αTC1 Clone 6 cells (CRL-2934) were purchased from American Type Culture Collection. Connaught Medical Research Laboratories 1066 medium, DMEM, HBSS, ι-glutamine, sodium pyruvate, Hepes, penicillin, and streptomycin were purchased from Thermo Fisher Scientific. FBS was purchased from HyClone. Trypsin (0.05% in 0.53 mM EDTA) was purchased from Corning. Human recombinant IL-1β and mouse IFN-γ were purchased from PeproTech. NG-Monomethyl-ʟ-arginine (NMMA) is from Enzo Life Sciences.

Islet isolation, islet and cell culture, and treatment

Islets were isolated from adult male C57BL6/J mice by collagenase digestion and were cultured at 37°C and 5% CO2 in Connaught Medical Research Laboratories supplemented with 10% heat-inactivated FBS and containing 5.5 mM glucose as previously described (Kelly et al, 2003). Islets from eight mice were pooled before separation into four samples for cytokine treatment. αTC1 Clone 6 cells were maintained at 37°C and 5% CO2, and were cultured in DMEM supplemented with 10% heat-inactivated FBS and containing 16.7 mM glucose, as previously described (Waas et al, 2020). Intact islets or adherent αTC1 cells were untreated or were treated for 6 h with 10 U/ml IL-1β alone, 150 U/ml IFN-γ alone, the combination of the two cytokines with or without 2 mM NMMA, or NMMA alone before dissociation and preparation for sequencing or before RNA isolation for qRT-PCR.

scRNA-seq

After cytokine treatment, islets were dissociated into single-cell suspensions by incubation in 0.48 mM EDTA in phosphate buffered saline followed by agitation in 1 mg/ml trypsin in Ca2+/Mg2+-free HBSS. Cells were filtered and resuspended in HBSS + 0.04% BSA. Single-cell suspensions were loaded into the Chromium Controller (10x Genomics). scRNA-seq libraries were prepared using the Chromium Single Cell 3′ v3 Reagent Kit (10x Genomics) according to manufacturer’s protocol. Libraries were loaded onto a NextSeq 500/550 High Output Kit v2.5 flow cell (150 cycles; Illumina) with the following conditions: 26 cycles for read 1, 98 cycles for read 2, and 8 cycles for the i7 index read. CellRanger (10x Genomics) functions “mkfastq” and “count” were used to demultiplex the sequencing data and generate gene-barcode matrices. Sequencing depth was as follows: Sample 1 (108,744,355), Sample 2 (42,013,628), Sample 3 (70,207,661), and Sample 4 (54,588,606). All scRNA-seq analysis was performed in R (version 3.6.1) using the package Seurat (version 3.1.0) (Butler et al, 2018). Number of genes detected per cell and percent of mitochondrial genes were plotted, and outlier cells were removed (number of genes less than 200 or greater than 4,000, or percent mitochondrial genes over 12%) to filter out doublets and cells with low read quality, leaving 7,750 of the original 12,571 cells. Cell cycle genes were regressed. Principal component analysis was performed, and the top 30 principal components were used for UMAP analysis, with clustering performed using the Louvain algorithm. All samples were normalized using Seurat’s default normalization settings. Briefly, reads in each cell for each gene were divided by the total number of reads within that cell, multiplied by a factor of 10,000, and transformed using the natural logarithm. Fold change values are expressed as natural log (ln) throughout the article, as this is the default method used by Seurat.

qRT-PCR

Total RNA was isolated from mouse islets or αTC1 cells using the RNeasy kit (QIAGEN). First-strand synthesis was performed using oligo(dT)s and Thermo Fisher Scientific Maxima H Minus reverse transcriptase per the manufacturer’s instructions. Quantitative PCR was performed using the SsoFast EvaGreen supermix (Bio-Rad) and a Bio-Rad CFX96 Real-Time system. Primers were purchased from Integrated DNA Technologies and are listed in Table S9. For relative quantification, gene expression was normalized to Gapdh using the comparative ΔCt method (Nolan et al, 2006).

Table S9 Primer sequences.

Statistical analysis

For differential expression analysis, P-values were calculated using the Wilcoxon test and adjusted to avoid false positives using Bonferroni correction. The threshold level of significance was set to adjusted P-value < 0.01 and ln(fold-change) > 0.5. For qRT-PCR analysis, statistical comparisons were made between groups using one-way ANOVA with Tukey post hoc test in GraphPad Prism software, and the threshold level of significance was set to P < 0.05.

Data Availability

Sequencing data from this publication have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus database under accession number GSE156175.

Acknowledgements

The authors thank Dr. Polly Hansen and Joshua Stafford (Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI) for technical assistance and Dr. Polly Hansen, Aaron Naatz, Joshua Stafford, and Chay Teng Yeo (Department of Biochemistry, Medical College of Wisconsin, Milwaukee, WI) for helpful discussions related to this project and for proofreading the manuscript. This research was completed in part with computational resources and technical support provided by the Research Computing Center at MCW. Funding: This work was supported by the National Institute of Diabetes and Digestive and Kidney Diseases grant DK-052194 (to JA Corbett), the National Institute of Allergy and Infectious Diseases grant AI-044458 (to JA Corbett) and AI-125741 and AI-148403 (to W Cui), the Juvenile Diabetes Research Foundation grant SRA-2019-829-S-B (to JA Corbett), and a grant from the Medical College of Wisconsin Cancer Center (to JA Corbett). JS Stancill was supported by the National Heart, Lung, and Blood grant T32-HL134643. MY Kasmani was supported by the National Institute of Diabetes and Digestive and Kidney Diseases grant DK-127526 and is a member of the Medical Scientist Training Program at the Medical College of Wisconsin, which is partially supported by a training grant from the National Institute of General Medical Sciences (T32-GM080202). This work was also supported by gifts from the Scott Tilton Foundation and from the Forest County Potawatomi Foundation.

Author Contributions

  • JS Stancill: conceptualization, data curation, formal analysis, supervision, validation, investigation, visualization, project administration, and writing—original draft, review, and editing.

  • MY Kasmani: software, formal analysis, visualization, methodology, and writing—original draft, review, and editing.

  • A Khatun: data curation, investigation, methodology, and writing—review and editing.

  • W Cui: resources, supervision, funding acquisition, methodology, project administration, and writing—original draft, review, and editing.

  • JA Corbett: conceptualization, resources, supervision, funding acquisition, project administration, and writing—original draft, review, and editing.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received October 26, 2020.
  • Revision received April 1, 2021.
  • Accepted April 7, 2021.
  • © 2021 Stancill et al.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Anderson SL,
    2. Carton JM,
    3. Lou J,
    4. Xing L,
    5. Rubin BY
    (1999) Interferon-induced guanylate binding protein-1 (gbp-1) mediates an antiviral effect against vesicular stomatitis virus and encephalomyocarditis virus. Virology 256: 8–14. doi:10.1006/viro.1999.9614
    OpenUrlCrossRefPubMed
  2. ↵
    1. Arnush M,
    2. Heitmeier MR,
    3. Scarim AL,
    4. Marino MH,
    5. Manning PT,
    6. Corbett JA
    (1998) Il-1 produced and released endogenously within human islets inhibits beta cell function. J Clin Invest 102: 516–526. doi:10.1172/JCI844
    OpenUrlCrossRefPubMed
  3. ↵
    1. Bader E,
    2. Migliorini A,
    3. Gegg M,
    4. Moruzzi N,
    5. Gerdes J,
    6. Roscioni SS,
    7. Bakhti M,
    8. Brandl E,
    9. Irmler M,
    10. Beckers J, et al.
    (2016) Identification of proliferative and mature beta-cells in the islets of langerhans. Nature 535: 430–434. doi:10.1038/nature18624
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bellmann K,
    2. Wenz A,
    3. Radons J,
    4. Burkart V,
    5. Kleemann R,
    6. Kolb H
    (1995) Heat shock induces resistance in rat pancreatic islet cells against nitric oxide, oxygen radicals and streptozotocin toxicity in vitro. J Clin Invest 95: 2840–2845. doi:10.1172/JCI117989
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bendtzen K,
    2. Mandrup-Poulsen T,
    3. Nerup J,
    4. Nielsen JH,
    5. Dinarello CA,
    6. Svenson M
    (1986) Cytotoxicity of human pi 7 interleukin-1 for pancreatic islets of langerhans. Science 232: 1545–1547. doi:10.1126/science.3086977
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Bergman B,
    2. Haskins K
    (1997) Autoreactive t-cell clones from the nonobese diabetic mouse. Proc Soc Exp Biol Med 214: 41–48. doi:10.3181/00379727-214-44067
    OpenUrlCrossRefPubMed
  7. ↵
    1. Blum B,
    2. Hrvatin S,
    3. Schuetz C,
    4. Bonal C,
    5. Rezania A,
    6. Melton DA
    (2012) Functional beta-cell maturation is marked by an increased glucose threshold and by expression of urocortin 3. Nat Biotechnol 30: 261–264. doi:10.1038/nbt.2141
    OpenUrlCrossRefPubMed
  8. ↵
    1. Bonner-Weir S,
    2. Orci L
    (1982) New perspectives on the microvasculature of the islets of langerhans in the rat. Diabetes 31: 883–889. doi:10.2337/diab.31.10.883
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Braun E,
    2. Hotter D,
    3. Koepke L,
    4. Zech F,
    5. Gross R,
    6. Sparrer KMJ,
    7. Muller JA,
    8. Pfaller CK,
    9. Heusinger E,
    10. Wombacher R, et al.
    (2019) Guanylate-binding proteins 2 and 5 exert broad antiviral activity by inhibiting furin-mediated processing of viral envelope proteins. Cell Rep 27: 2092–2104.e10. doi:10.1016/j.celrep.2019.04.063
    OpenUrlCrossRefPubMed
  10. ↵
    1. Burkart V,
    2. Liu H,
    3. Bellmann K,
    4. Wissing D,
    5. Jaattela M,
    6. Cavallo MG,
    7. Pozzilli P,
    8. Briviba K,
    9. Kolb H
    (2000) Natural resistance of human beta cells toward nitric oxide is mediated by heat shock protein 70. J Biol Chem 275: 19521–19528. doi:10.1074/jbc.M002265200
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Butler A,
    2. Hoffman P,
    3. Smibert P,
    4. Papalexi E,
    5. Satija R
    (2018) Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol 36: 411–420. doi:10.1038/nbt.4096
    OpenUrlCrossRefPubMed
  12. ↵
    1. Cardozo AK,
    2. Heimberg H,
    3. Heremans Y,
    4. Leeman R,
    5. Kutlu B,
    6. Kruhoffer M,
    7. Orntoft T,
    8. Eizirik DL
    (2001a) A comprehensive analysis of cytokine-induced and nuclear factor-kappa b-dependent genes in primary rat pancreatic beta-cells. J Biol Chem 276: 48879–48886. doi:10.1074/jbc.M108658200
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Cardozo AK,
    2. Kruhoffer M,
    3. Leeman R,
    4. Orntoft T,
    5. Eizirik DL
    (2001b) Identification of novel cytokine-induced genes in pancreatic beta-cells by high-density oligonucleotide arrays. Diabetes 50: 909–920. doi:10.2337/diabetes.50.5.909
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Cardozo AK,
    2. Proost P,
    3. Gysemans C,
    4. Chen MC,
    5. Mathieu C,
    6. Eizirik DL
    (2003) Il-1beta and ifn-gamma induce the expression of diverse chemokines and il-15 in human and rat pancreatic islet cells, and in islets from pre-diabetic nod mice. Diabetologia 46: 255–266. doi:10.1007/s00125-002-1017-0
    OpenUrlCrossRefPubMed
  15. ↵
    1. Corbett JA,
    2. Lancaster JR Jr..,
    3. Sweetland MA,
    4. McDaniel ML
    (1991) Interleukin-1 beta-induced formation of epr-detectable iron-nitrosyl complexes in islets of langerhans. Role of nitric oxide in interleukin-1 beta-induced inhibition of insulin secretion. J Biol Chem 266: 21351–21354. doi:10.1016/s0021-9258(18)54642-1
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Corbett JA,
    2. McDaniel ML
    (1992) Does nitric oxide mediate autoimmune destruction of beta-cells? Possible therapeutic interventions in iddm. Diabetes 41: 897–903. doi:10.2337/diab.41.8.897
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Corbett JA,
    2. McDaniel ML
    (1995) Intraislet release of interleukin 1 inhibits beta cell function by inducing beta cell expression of inducible nitric oxide synthase. J Exp Med 181: 559–568. doi:10.1084/jem.181.2.559
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Corbett JA,
    2. Sweetland MA,
    3. Wang JL,
    4. Lancaster JR Jr..,
    5. McDaniel ML
    (1993) Nitric oxide mediates cytokine-induced inhibition of insulin secretion by human islets of langerhans. Proc Natl Acad Sci U S A 90: 1731–1735. doi:10.1073/pnas.90.5.1731
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Corbett JA,
    2. Wang JL,
    3. Hughes JH,
    4. Wolf BA,
    5. Sweetland MA,
    6. Lancaster JR Jr..,
    7. McDaniel ML
    (1992a) Nitric oxide and cyclic gmp formation induced by interleukin 1 beta in islets of langerhans. Evidence for an effector role of nitric oxide in islet dysfunction. Biochem J 287: 229–235. doi:10.1042/bj2870229
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Corbett JA,
    2. Wang JL,
    3. Sweetland MA,
    4. Lancaster JR Jr..,
    5. McDaniel ML
    (1992b) Interleukin 1 beta induces the formation of nitric oxide by beta-cells purified from rodent islets of langerhans. Evidence for the beta-cell as a source and site of action of nitric oxide. J Clin Invest 90: 2384–2391. doi:10.1172/JCI116129
    OpenUrlCrossRefPubMed
  21. ↵
    1. Cunningham JM,
    2. Green IC
    (1994) Cytokines, nitric oxide and insulin secreting cells. Growth Regul 4: 173–180.
    OpenUrlPubMed
  22. ↵
    1. Dahan T,
    2. Ziv O,
    3. Horwitz E,
    4. Zemmour H,
    5. Lavi J,
    6. Swisa A,
    7. Leibowitz G,
    8. Ashcroft FM,
    9. In’t Veld P,
    10. Glaser B, et al.
    (2017) Pancreatic beta-cells express the fetal islet hormone gastrin in rodent and human diabetes. Diabetes 66: 426–436. doi:10.2337/db16-0641
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Darville MI,
    2. Eizirik DL
    (1998) Regulation by cytokines of the inducible nitric oxide synthase promoter in insulin-producing cells. Diabetologia 41: 1101–1108. doi:10.1007/s001250051036
    OpenUrlCrossRefPubMed
  24. ↵
    1. Dennis G Jr..,
    2. Sherman BT,
    3. Hosack DA,
    4. Yang J,
    5. Gao W,
    6. Lane HC,
    7. Lempicki RA
    (2003) David: Database for annotation, visualization, and integrated discovery. Genome Biol 4: P3. doi:10.1186/gb-2003-4-5-p3
    OpenUrlCrossRefPubMed
  25. ↵
    1. Dorrell C,
    2. Schug J,
    3. Canaday PS,
    4. Russ HA,
    5. Tarlow BD,
    6. Grompe MT,
    7. Horton T,
    8. Hebrok M,
    9. Streeter PR,
    10. Kaestner KH, et al.
    (2016) Human islets contain four distinct subtypes of beta cells. Nat Commun 7: 11756. doi:10.1038/ncomms11756
    OpenUrlCrossRefPubMed
  26. ↵
    1. Eizirik DL,
    2. Cagliero E,
    3. Bjorklund A,
    4. Welsh N
    (1992) Interleukin-1 beta induces the expression of an isoform of nitric oxide synthase in insulin-producing cells, which is similar to that observed in activated macrophages. FEBS Lett 308: 249–252. doi:10.1016/0014-5793(92)81285-t
    OpenUrlCrossRefPubMed
  27. ↵
    1. Eizirik DL,
    2. Flodstrom M,
    3. Karlsen AE,
    4. Welsh N
    (1996) The harmony of the spheres: Inducible nitric oxide synthase and related genes in pancreatic beta cells. Diabetologia 39: 875–890. doi:10.1007/BF00403906
    OpenUrlCrossRefPubMed
  28. ↵
    1. Eizirik DL,
    2. Pasquali L,
    3. Cnop M
    (2020) Pancreatic beta-cells in type 1 and type 2 diabetes mellitus: Different pathways to failure. Nat Rev Endocrinol 16: 349–362. doi:10.1038/s41574-020-0355-7
    OpenUrlCrossRef
  29. ↵
    1. Eizirik DL,
    2. Sandler S,
    3. Welsh N,
    4. Cetkovic-Cvrlje M,
    5. Nieman A,
    6. Geller DA,
    7. Pipeleers DG,
    8. Bendtzen K,
    9. Hellerstrom C
    (1994) Cytokines suppress human islet function irrespective of their effects on nitric oxide generation. J Clin Invest 93: 1968–1974. doi:10.1172/JCI117188
    OpenUrlCrossRefPubMed
  30. ↵
    1. Flodstrom M,
    2. Welsh N,
    3. Eizirik DL
    (1996) Cytokines activate the nuclear factor kappa b (nf-kappa b) and induce nitric oxide production in human pancreatic islets. FEBS Lett 385: 4–6. doi:10.1016/0014-5793(96)00337-7
    OpenUrlCrossRefPubMed
  31. ↵
    1. Gao T,
    2. McKenna B,
    3. Li C,
    4. Reichert M,
    5. Nguyen J,
    6. Singh T,
    7. Yang C,
    8. Pannikar A,
    9. Doliba N,
    10. Zhang T, et al.
    (2014) Pdx1 maintains beta cell identity and function by repressing an alpha cell program. Cell Metab 19: 259–271. doi:10.1016/j.cmet.2013.12.002
    OpenUrlCrossRefPubMed
  32. ↵
    1. Gepts W
    (1965) Pathologic anatomy of the pancreas in juvenile diabetes mellitus. Diabetes 14: 619–633. doi:10.2337/diab.14.10.619
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Grunnet LG,
    2. Aikin R,
    3. Tonnesen MF,
    4. Paraskevas S,
    5. Blaabjerg L,
    6. Storling J,
    7. Rosenberg L,
    8. Billestrup N,
    9. Maysinger D,
    10. Mandrup-Poulsen T
    (2009) Proinflammatory cytokines activate the intrinsic apoptotic pathway in beta-cells. Diabetes 58: 1807–1815. doi:10.2337/db08-0178
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Guillam MT,
    2. Hummler E,
    3. Schaerer E,
    4. Yeh JI,
    5. Birnbaum MJ,
    6. Beermann F,
    7. Schmidt A,
    8. Deriaz N,
    9. Thorens B
    (1997) Early diabetes and abnormal postnatal pancreatic islet development in mice lacking glut-2. Nat Genet 17: 327–330. doi:10.1038/ng1197-327
    OpenUrlCrossRefPubMed
  35. ↵
    1. Hardy OT,
    2. Hohmeier HE,
    3. Becker TC,
    4. Manduchi E,
    5. Doliba NM,
    6. Gupta RK,
    7. White P,
    8. Stoeckert CJ Jr..,
    9. Matschinsky FM,
    10. Newgard CB, et al.
    (2007) Functional genomics of the beta-cell: Short-chain 3-hydroxyacyl-coenzyme a dehydrogenase regulates insulin secretion independent of k+ currents. Mol Endocrinol 21: 765–773. doi:10.1210/me.2006-0411
    OpenUrlCrossRefPubMed
  36. ↵
    1. Heitmeier MR,
    2. Scarim AL,
    3. Corbett JA
    (1997) Interferon-gamma increases the sensitivity of islets of langerhans for inducible nitric-oxide synthase expression induced by interleukin 1. J Biol Chem 272: 13697–13704. doi:10.1074/jbc.272.21.13697
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Helqvist S,
    2. Polla BS,
    3. Johannesen J,
    4. Nerup J
    (1991) Heat shock protein induction in rat pancreatic islets by recombinant human interleukin 1 beta. Diabetologia 34: 150–156. doi:10.1007/BF00418268
    OpenUrlCrossRefPubMed
  38. ↵
    1. Henseleit KD,
    2. Nelson SB,
    3. Kuhlbrodt K,
    4. Hennings JC,
    5. Ericson J,
    6. Sander M
    (2005) Nkx6 transcription factor activity is required for alpha- and beta-cell development in the pancreas. Development 132: 3139–3149. doi:10.1242/dev.01875
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Hughes KJ,
    2. Meares GP,
    3. Hansen PA,
    4. Corbett JA
    (2011) Foxo1 and sirt1 regulate beta-cell responses to nitric oxide. J Biol Chem 286: 8338–8348. doi:10.1074/jbc.M110.204768
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Kelly CB,
    2. Blair LA,
    3. Corbett JA,
    4. Scarim AL
    (2003) Isolation of islets of langerhans from rodent pancreas. Methods Mol Med 83: 3–14. doi:10.1385/1-59259-377-1:003
    OpenUrlCrossRefPubMed
  41. ↵
    1. Kim-Muller JY,
    2. Fan J,
    3. Kim YJ,
    4. Lee SA,
    5. Ishida E,
    6. Blaner WS,
    7. Accili D
    (2016) Aldehyde dehydrogenase 1a3 defines a subset of failing pancreatic beta cells in diabetic mice. Nat Commun 7: 12631. doi:10.1038/ncomms12631
    OpenUrlCrossRefPubMed
  42. ↵
    1. Kwon G,
    2. Corbett JA,
    3. Rodi CP,
    4. Sullivan P,
    5. McDaniel ML
    (1995) Interleukin-1 beta-induced nitric oxide synthase expression by rat pancreatic beta-cells: Evidence for the involvement of nuclear factor kappa b in the signaling mechanism. Endocrinology 136: 4790–4795. doi:10.1210/endo.136.11.7588208
    OpenUrlCrossRefPubMed
  43. ↵
    1. Larsen CM,
    2. Wadt KA,
    3. Juhl LF,
    4. Andersen HU,
    5. Karlsen AE,
    6. Su MS,
    7. Seedorf K,
    8. Shapiro L,
    9. Dinarello CA,
    10. Mandrup-Poulsen T
    (1998) Interleukin-1beta-induced rat pancreatic islet nitric oxide synthesis requires both the p38 and extracellular signal-regulated kinase 1/2 mitogen-activated protein kinases. J Biol Chem 273: 15294–15300. doi:10.1074/jbc.273.24.15294
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Like AA,
    2. Weringer EJ,
    3. Holdash A,
    4. McGill P,
    5. Atkinson D,
    6. Rossini AA
    (1985) Adoptive transfer of autoimmune diabetes mellitus in biobreeding/worcester (bb/w) inbred and hybrid rats. J Immunol 134: 1583–1587.
    OpenUrlAbstract/FREE Full Text
  45. ↵
    1. Mamane Y,
    2. Heylbroeck C,
    3. Genin P,
    4. Algarte M,
    5. Servant MJ,
    6. LePage C,
    7. DeLuca C,
    8. Kwon H,
    9. Lin R,
    10. Hiscott J
    (1999) Interferon regulatory factors: The next generation. Gene 237: 1–14. doi:10.1016/s0378-1119(99)00262-0
    OpenUrlCrossRefPubMed
  46. ↵
    1. Mandrup-Poulsen T
    (1996) The role of interleukin-1 in the pathogenesis of iddm. Diabetologia 39: 1005–1029. doi:10.1007/BF00400649
    OpenUrlCrossRefPubMed
  47. ↵
    1. Mandrup-Poulsen T,
    2. Bendtzen K,
    3. Nielsen JH,
    4. Bendixen G,
    5. Nerup J
    (1985) Cytokines cause functional and structural damage to isolated islets of langerhans. Allergy 40: 424–429. doi:10.1111/j.1398-9995.1985.tb02681.x
    OpenUrlCrossRefPubMed
  48. ↵
    1. Matsuoka TA,
    2. Artner I,
    3. Henderson E,
    4. Means A,
    5. Sander M,
    6. Stein R
    (2004) The mafa transcription factor appears to be responsible for tissue-specific expression of insulin. Proc Natl Acad Sci U S A 101: 2930–2933. doi:10.1073/pnas.0306233101
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Meares GP,
    2. Fontanilla D,
    3. Broniowska KA,
    4. Andreone T,
    5. Lancaster JR Jr..,
    6. Corbett JA
    (2013) Differential responses of pancreatic beta-cells to ros and rns. Am J Physiol Endocrinol Metab 304: E614–E622. doi:10.1152/ajpendo.00424.2012
    OpenUrlCrossRefPubMed
  50. ↵
    1. Nolan T,
    2. Hands RE,
    3. Bustin SA
    (2006) Quantification of mrna using real-time rt-pcr. Nat Protoc 1: 1559–1582. doi:10.1038/nprot.2006.236
    OpenUrlCrossRefPubMed
  51. ↵
    1. Nordmann TM,
    2. Dror E,
    3. Schulze F,
    4. Traub S,
    5. Berishvili E,
    6. Barbieux C,
    7. Boni-Schnetzler M,
    8. Donath MY
    (2017) The role of inflammation in beta-cell dedifferentiation. Sci Rep 7: 6285. doi:10.1038/s41598-017-06731-w
    OpenUrlCrossRef
  52. ↵
    1. Oleson BJ,
    2. Broniowska KA,
    3. Naatz A,
    4. Hogg N,
    5. Tarakanova VL,
    6. Corbett JA
    (2016) Nitric oxide suppresses beta-cell apoptosis by inhibiting the DNA damage response. Mol Cell Biol 36: 2067–2077. doi:10.1128/MCB.00262-16
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Ortis F,
    2. Naamane N,
    3. Flamez D,
    4. Ladriere L,
    5. Moore F,
    6. Cunha DA,
    7. Colli ML,
    8. Thykjaer T,
    9. Thorsen K,
    10. Orntoft TF, et al.
    (2010) Cytokines interleukin-1beta and tumor necrosis factor-alpha regulate different transcriptional and alternative splicing networks in primary beta-cells. Diabetes 59: 358–374. doi:10.2337/db09-1159
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Padgett LE,
    2. Broniowska KA,
    3. Hansen PA,
    4. Corbett JA,
    5. Tse HM
    (2013) The role of reactive oxygen species and proinflammatory cytokines in type 1 diabetes pathogenesis. Ann N Y Acad Sci 1281: 16–35. doi:10.1111/j.1749-6632.2012.06826.x
    OpenUrlCrossRefPubMed
  55. ↵
    1. Pipeleers D,
    2. Kiekens R,
    3. Ling Z,
    4. Wilikens A,
    5. Schuit F
    (1994) Physiologic relevance of heterogeneity in the pancreatic beta-cell population. Diabetologia 37: S57–S64. doi:10.1007/BF00400827
    OpenUrlCrossRefPubMed
  56. ↵
    1. Saldeen J,
    2. Welsh N
    (1994) Interleukin-1 beta induced activation of nf-kappa b in insulin producing rinm5f cells is prevented by the protease inhibitor n alpha-p-tosyl-l-lysine chloromethylketone. Biochem Biophys Res Commun 203: 149–155. doi:10.1006/bbrc.1994.2161
    OpenUrlCrossRefPubMed
  57. ↵
    1. Samuel CE
    (2001) Antiviral actions of interferons. Clin Microbiol Rev 14: 778–809. doi:10.1128/CMR.14.4.778-809.2001
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Scarim AL,
    2. Heitmeier MR,
    3. Corbett JA
    (1998) Heat shock inhibits cytokine-induced nitric oxide synthase expression by rat and human islets. Endocrinology 139: 5050–5057. doi:10.1210/endo.139.12.6366
    OpenUrlCrossRefPubMed
  59. ↵
    1. Schneider WM,
    2. Chevillotte MD,
    3. Rice CM
    (2014) Interferon-stimulated genes: A complex web of host defenses. Annu Rev Immunol 32: 513–545. doi:10.1146/annurev-immunol-032713-120231
    OpenUrlCrossRefPubMed
  60. ↵
    1. Semple F,
    2. Dorin JR
    (2012) Beta-defensins: Multifunctional modulators of infection, inflammation and more? J Innate Immun 4: 337–348. doi:10.1159/000336619
    OpenUrlCrossRefPubMed
  61. ↵
    1. Setum CM,
    2. Serie JR,
    3. Hegre OD
    (1991) Confocal microscopic analysis of the nonendocrine cellular component of isolated adult rat islets of langerhans. Transplantation 51: 1131–1133. doi:10.1097/00007890-199105000-00043
    OpenUrlCrossRefPubMed
  62. ↵
    1. Southern C,
    2. Schulster D,
    3. Green IC
    (1990) Inhibition of insulin secretion by interleukin-1 beta and tumour necrosis factor-alpha via an l-arginine-dependent nitric oxide generating mechanism. FEBS Lett 276: 42–44. doi:10.1016/0014-5793(90)80502-a
    OpenUrlCrossRefPubMed
  63. ↵
    1. Stafford JD,
    2. Shaheen ZR,
    3. Yeo CT,
    4. Corbett JA
    (2020) Inhibition of mitochondrial oxidative metabolism attenuates emcv replication and protects beta-cells from virally mediated lysis. J Biol Chem 295: 16655–16664. doi:10.1074/jbc.RA120.014851
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Steer SA,
    2. Scarim AL,
    3. Chambers KT,
    4. Corbett JA
    (2006) Interleukin-1 stimulates beta-cell necrosis and release of the immunological adjuvant hmgb1. PLoS Med 3: e17. doi:10.1371/journal.pmed.0030017
    OpenUrlCrossRefPubMed
  65. ↵
    1. Strandell E,
    2. Buschard K,
    3. Saldeen J,
    4. Welsh N
    (1995) Interleukin-1 beta induces the expression of hsp70, heme oxygenase and mn-sod in facs-purified rat islet beta-cells, but not in alpha-cells. Immunol Lett 48: 145–148. doi:10.1016/0165-2478(95)02459-x
    OpenUrlCrossRefPubMed
  66. ↵
    1. Talchai C,
    2. Xuan S,
    3. Lin HV,
    4. Sussel L,
    5. Accili D
    (2012) Pancreatic beta cell dedifferentiation as a mechanism of diabetic beta cell failure. Cell 150: 1223–1234. doi:10.1016/j.cell.2012.07.029
    OpenUrlCrossRefPubMed
  67. ↵
    1. Waas M,
    2. Snarrenberg ST,
    3. Littrell J,
    4. Lipinski RAJ,
    5. Hansen PA,
    6. Corbett JA,
    7. Gundry RL
    (2020) Surfacegenie: A web-based application for prioritizing cell-type-specific marker candidates. Bioinformatics 36: 3447–3456. doi:10.1093/bioinformatics/btaa092
    OpenUrlCrossRef
  68. ↵
    1. Weber SM,
    2. Chambers KT,
    3. Bensch KG,
    4. Scarim AL,
    5. Corbett JA
    (2004a) Ppargamma ligands induce er stress in pancreatic beta-cells: Er stress activation results in attenuation of cytokine signaling. Am J Physiol Endocrinol Metab 287: E1171–E1177. doi:10.1152/ajpendo.00331.2004
    OpenUrlCrossRefPubMed
  69. ↵
    1. Weber SM,
    2. Scarim AL,
    3. Corbett JA
    (2003) Inhibition of ifn-gamma -induced stat1 activation by 15- deoxy-delta 12,14-prostaglandin j2. Am J Physiol Endocrinol Metab 284: E883–E891. doi:10.1152/ajpendo.00515.2002
    OpenUrlCrossRefPubMed
  70. ↵
    1. Weber SM,
    2. Scarim AL,
    3. Corbett JA
    (2004b) Ppargamma is not required for the inhibitory actions of pgj2 on cytokine signaling in pancreatic beta-cells. Am J Physiol Endocrinol Metab 286: E329–E336. doi:10.1152/ajpendo.00392.2003
    OpenUrlCrossRefPubMed
  71. ↵
    1. Welsh N,
    2. Eizirik DL,
    3. Bendtzen K,
    4. Sandler S
    (1991) Interleukin-1 beta-induced nitric oxide production in isolated rat pancreatic islets requires gene transcription and may lead to inhibition of the krebs cycle enzyme aconitase. Endocrinology 129: 3167–3173. doi:10.1210/endo-129-6-3167
    OpenUrlCrossRefPubMed
  72. ↵
    1. Welsh N,
    2. Margulis B,
    3. Borg LA,
    4. Wiklund HJ,
    5. Saldeen J,
    6. Flodstrom M,
    7. Mello MA,
    8. Andersson A,
    9. Pipeleers DG,
    10. Hellerstrom C, et al.
    (1995) Differences in the expression of heat-shock proteins and antioxidant enzymes between human and rodent pancreatic islets: Implications for the pathogenesis of insulin-dependent diabetes mellitus. Mol Med 1: 806–820. doi:10.1007/bf03401895
    OpenUrlCrossRefPubMed
  73. ↵
    1. Wicker LS,
    2. Miller BJ,
    3. Mullen Y
    (1986) Transfer of autoimmune diabetes mellitus with splenocytes from nonobese diabetic (nod) mice. Diabetes 35: 855–860. doi:10.2337/diab.35.8.855
    OpenUrlAbstract/FREE Full Text
  74. ↵
    1. Wilcox CL,
    2. Terry NA,
    3. Walp ER,
    4. Lee RA,
    5. May CL
    (2013) Pancreatic alpha-cell specific deletion of mouse arx leads to alpha-cell identity loss. PLoS One 8: e66214. doi:10.1371/journal.pone.0066214
    OpenUrlCrossRefPubMed
  75. ↵
    1. Wright JR Jr..,
    2. Lacy PE
    (1989) Silica prevents the induction of diabetes with complete freund’s adjuvant and low-dose streptozotocin in rats. Diabetes Res 11: 51–54.
    OpenUrlPubMed
  76. ↵
    1. Wright JR Jr..,
    2. Lefkowith JB,
    3. Schreiner G,
    4. Lacy PE
    (1988) Essential fatty acid deficiency prevents multiple low-dose streptozotocin-induced diabetes in cd-1 mice. Proc Natl Acad Sci U S A 85: 6137–6141. doi:10.1073/pnas.85.16.6137
    OpenUrlAbstract/FREE Full Text
  77. ↵
    1. Xie QW,
    2. Whisnant R,
    3. Nathan C
    (1993) Promoter of the mouse gene encoding calcium-independent nitric oxide synthase confers inducibility by interferon gamma and bacterial lipopolysaccharide. J Exp Med 177: 1779–1784. doi:10.1084/jem.177.6.1779
    OpenUrlAbstract/FREE Full Text
  78. ↵
    1. Ylipaasto P,
    2. Kutlu B,
    3. Rasilainen S,
    4. Rasschaert J,
    5. Salmela K,
    6. Teerijoki H,
    7. Korsgren O,
    8. Lahesmaa R,
    9. Hovi T,
    10. Eizirik DL, et al.
    (2005) Global profiling of coxsackievirus- and cytokine-induced gene expression in human pancreatic islets. Diabetologia 48: 1510–1522. doi:10.1007/s00125-005-1839-7
    OpenUrlCrossRefPubMed
  79. ↵
    1. Zhang J,
    2. McKenna LB,
    3. Bogue CW,
    4. Kaestner KH
    (2014) The diabetes gene hhex maintains delta-cell differentiation and islet function. Genes Dev 28: 829–834. doi:10.1101/gad.235499.113
    OpenUrlAbstract/FREE Full Text
  80. ↵
    1. Zhu J,
    2. Zhang Y,
    3. Ghosh A,
    4. Cuevas RA,
    5. Forero A,
    6. Dhar J,
    7. Ibsen MS,
    8. Schmid-Burgk JL,
    9. Schmidt T,
    10. Ganapathiraju MK, et al.
    (2014) Antiviral activity of human oasl protein is mediated by enhancing signaling of the rig-i rna sensor. Immunity 40: 936–948. doi:10.1016/j.immuni.2014.05.007
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Single-cell RNA sequencing of mouse islets exposed to proinflammatory cytokines
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Single islet cell responses to cytokines
Jennifer S Stancill, Moujtaba Y Kasmani, Achia Khatun, Weiguo Cui, John A Corbett
Life Science Alliance Apr 2021, 4 (6) e202000949; DOI: 10.26508/lsa.202000949

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Single islet cell responses to cytokines
Jennifer S Stancill, Moujtaba Y Kasmani, Achia Khatun, Weiguo Cui, John A Corbett
Life Science Alliance Apr 2021, 4 (6) e202000949; DOI: 10.26508/lsa.202000949
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 4, No. 6
June 2021
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Data Availability
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Genomics & Functional Genomics
  • Immunology
  • Molecular Biology

Related Articles

  • No related articles found.

Cited By...

  • No citing articles found.
  • Google Scholar

More in this TOC Section

  • NLRP3 controls ATM activation
  • RIF1 stabilizes replication domains
  • Chromatin context-dependent interferon response
Show more Research Article

Similar Articles

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2023 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.