Main

The response of a tissue to stress/injury can involve cell death and proliferation. However, it has become increasingly recognized that changes in cellular differentiation state can have an important role.1 In type I diabetes (T1D), the established view has been that the primary pathophysiological event is β-cell apoptosis due to a β-cell specific autoimmune response,2 leading to profound β-cell deficiency. Thus, there has been great interest in inducing β-cell neogenesis, but there has been controversy over how and even whether β-cell regeneration occurs.3

Activation of dedicated stem/progenitor cells within the pancreas and transdifferentiation of other differentiated cell types to β-cells are two potential mechanisms. In the past, the prevailing paradigm has been that neogenesis proceeds by the activation of facultative β-cell stem/progenitors within pancreatic ducts.4, 5, 6, 7 However, more recent studies have not found evidence of robust β-cell neogenesis from ducts.8, 9, 10, 11 β-cell neogenesis from other cell types within the pancreas, including acinar12 and centroacinar13 cells has also been reported.

Recently, we demonstrated robust β-cell neogenesis by transdifferentiation from preexisting α-cells in a model of T1D where severe β-cell deficiency was induced by high-dose alloxan.14, 15 In this model, β-cell neogenesis from α-cells was stimulated by pancreatic duct ligation (PDL).9, 15 Surgical reversal of PDL led to the recovery of β-cell mass and function by a combination of β-cell replication and β-cell neogenesis, demonstrating that β-cell regeneration by α- to β-cell neogenesis could be a robust approach to diabetes therapy,16 but PDL, which involves major surgery, is not a practical approach to therapy. Importantly, the relevance of α- to β-cell transdifferentiation to human biology remained unclear, as previous studies were performed in rodents.

Here, we report the occurrence of efficient islet cell transdifferentiation using an entirely pharmacologically based approach where the peptide caerulein,17, 18 substituting for PDL, stimulated β-cell transdifferentiation from α-cells in mice rendered severely β-cell deficient by alloxan. Following caerulein plus alloxan, many of the neogenic β-cells went on to form δ-cells. In murine and human T1D, a similar process appeared to occur, where α-cells transdifferentiated into β-cells, which went on to form δ-cells. This led to a marked δ-cell excess in both murine and human T1D.

The finding of endocrine cell transdifferentiation in T1D supports a new paradigm where β-cells, in addition to undergoing destruction by inflammatory mediators, undergo a dynamic process of neogenesis from α-cells and transdifferentiation to δ-cells. Controlling the neogenic process could lead to a new approach to diabetes therapy.

Results

Caerulein plus alloxan induced α-cells transdifferentiation into β-cells

As expected,17, 19 caerulein-induced severe disruption of the acinar pancreas (Supplementary Figures 1A–D) associated with almost the complete loss of the acinar cell marker amylase (compare Supplementary Figures 1E and G). This was transient, with recovery by 11 days post caerulein injection (Supplementary Figures 1A–D). Caerulin injection induced marked infiltration with CD3-positive lymphocytes, in contrast with alloxan injection that did not induce lymphocytic infiltration (Supplementary Figures 1E–H). Caerulein had no effect on islet morphology or hormone expression (Supplementary Figures 1A–E). Alloxan induced β-cell ablation, which we found previously to be required for endocrine cell transdifferentiation9 (Supplementary Figures 1G and H).

In our previous study of α- to β-cell transdifferentiation, employing PDL plus alloxan to induce transdifferentiation,9 the evidence for transdifferentiation as a cause of the increased β-cell mass was the finding of a high frequency of intermediates with glucagon–insulin coexpression, as well as the presence of little β-cell replication or neogenesis within ducts.9 The high mortality following PDL plus alloxan hampered the genetic lineage tracing, which is considered to be the gold standard for determining lineage relationships. However, the alloxan plus caerulein model had lower mortality, allowing us to perform such studies.

For genetic lineage tracing of α-cells, we used mice expressing membrane-localized enhanced green fluorescent protein (eGFP; membrane-localized enhanced green fluorescent protein (mG)) on their α-cells while expressing membrane-localized Tomato (mT) on all other cells (Glu-mT/mG; Figure 1a). Pancreases from Glu-mT/mG mice treated with alloxan alone or alloxan plus caerulein were analyzed 20 and 38 days post alloxan injection. In mice treated with alloxan alone, the few remaining β-cells expressed mT. Following alloxan plus caerulein, insulin-positive cells expressing mG, indicating an origin from recombined α-cells, were easily detected (Figures 1b–e, h, quantified in i, Supplementary Movie 1). At day 20, 5.5% of mG-positive cells were positive for insulin, which translated to 18% of all β-cells. All such cells were in islets. No mG-positive cells expressing glucagon or insulin were ever found in the ducts (Supplementary Figure 2). Following alloxan alone, mG-positive cells expressing insulin were extremely rare (1% at day 20). At day 38 post caerulein plus alloxan, 9% of β-cells were mG-positive, compared with 1% following alloxan alone (Figure 1i). Considering that the α-cell recombination efficiency in this experiment was 62% (Figure 1), approximately 29% of β-cells present at 20 days following alloxan plus caerulein arose from α-cells. Thus, caerulein stimulated a dramatic increase in α- to β-cell transdifferentiation.

Figure 1
figure 1

Genetic lineage tracing demonstrated transdifferentiation of α-cells to β- and δ-cells following alloxan plus caerulein. (a) Glu-mT/mG mice were injected with caerulein and alloxan as described in the Materials and Methods section. The recombination efficiency in α-cells was 0.62+/−0.047 (S.E.M.), calculated as the number of α-cells expressing eGFP/α-cells expressing eGFP(green) plus mT (red). (b and c) Representative low-power photomicrographs of islets from mice killed at days 20 and 38, respectively, showing transdifferentiated β- (blue) and δ-cells (white). (dh) High-power views of the regions from b and c showing lineage-marked cells expressing insulin (d and h), somatostatin (f and g) or both (e). (i) Quantification of β- and δ- cell transdifferentiation from α-cells. Note that no eGFP-positive cells expressing somatostatin were detected after alloxan alone. (j) Quantification of the frequency of islet cell types during alloxan alone or alloxan+caerulein at days 20 and 38. For the day 20 time point, 317 cells islet cells were quantified in the alloxan alone condition and 605 cells were quantified in the alloxan plus caerulein condition. For the day 38 time point, 309 cells islet cells were quantified in the alloxan alone condition and 478 cells were quantified in the alloxan plus caerulein condition. Error bars are S.E.M. Scale bar for cg=10 μm

Caerulein-induced transdifferentiation of α-cells into δ-cells through a β-cell intermediate

A comparison between days 20 and 38 of the frequency of insulin-positive and somatostatin-positive cells revealed a substantial increase over time in mG-positive cells expressing somatostatin and a corresponding decrease in cells expressing insulin (Figures 1f and g, quantified in i). Following alloxan plus caerulein, somatostatin-positive cells expressing mG, indicating an origin from recombined α-cells, rose from 2.1% of mG-positive cells at day 20 to 8.2% at day 38 (Figures 1b–e and h, quantified in i, Supplementary Movie 1), but were absent from samples from mice treated with alloxan alone. A cell coexpressing mG, insulin and somatostatin is shown in Figure 1e, consistent with transdifferentiated β-cells going on to form δ-cells. Following alloxan alone, no cells coexpressing mG and somatostatin were found (Figure 1i and Supplementary Movie 1).

Transdifferentiation of α- cells induced by alloxan plus caerulein proceeded sequentially, with β-cells being formed initially, followed by β- to δ-cell transdifferentiation

Transdifferentiation of α-cells to β- and δ-cells could have arisen by two possible routes; either as two independent transdifferentiation processes: α- to β-cells and α- to δ-cells in parallel, or as a linear process where α-cells convert to β- and then to δ-cells. To distinguish between those possibilities, we undertook a time course analysis to evaluate the frequency of intermediates in the process of transdifferentiation. This was done with C57Bl/6 mice, with a higher dose of caerulein (50 μg/kg/injection), increasing the number of transdifferentiating cells from 12+/−2.3 (Figure 1j) to 22+/−4 (Figure 2j).

Figure 2
figure 2

Time course analysis demonstrated the transdifferentiation from α- to β- to δ-cells following alloxan plus caerulein. Mice were killed and pancreas harvested at days 13, 20, 27 and 38 following the initiation of alloxan +/− caerulein treatment on day 1. (a) A typical islet of Langerhans from an untreated C57BL/6 mouse demonstrating a preponderance of β-cells (insulin-red) located internal to a ring of α-cells (glucagon-green) and δ-cells (somatostatin-white). DAPI (blue nuclear stain). (be) Islets from mice injected with alloxan alone demonstrated continued preponderance of α-cells, with few δ-cells or residual β-cells over time. (fi) Islets from mice injected with alloxan plus caerulein demonstrated a dramatic increase of β- and δ-cells over time. (j) Quantification of insulin-positive area over time, demonstrating a peak of β-cell formation at day 20 following alloxan plus caerulein (red line), with no significant change following alloxan alone (black line). (k) Quantification of the somatostatin-positive area normalized to the glucagon-positive area to control for effects of β-cell ablation on the somatostatin-positive area, demonstrating a continued increase of δ-cell formation over time following alloxan plus caerulein, with no significant change following alloxan alone. (l) Quantification of glucagon-positive area over time, demonstrating a decrease in the number of α-cells following alloxan plus caerulein (red line), with no significant change following alloxan alone. (m) Quantification of insulin–glucagon copositive area over time, demonstrating an increase of cells coexpressing insulin and glucagon following alloxan plus caerulein until day 38, when such cells were no longer present. There was no significant change in insulin–glucagon copositive cells following alloxan alone, never rising above the background. (n) Quantification of insulin–somatostatin copositive area over time, demonstrating an increase of cells coexpressing insulin and somatostatin following alloxan plus caerulein, peaking at day 27 and disappearing by day 38. There was no significant change in insulin–somatostatin copositive cells following alloxan alone, never rising above the background. Error bars are S.E.M. *P<0.05, **P<0.005, ***P<0.0005. Scale bar=75 μm

Mice were killed at days 13, 20, 27 and 38 following alloxan plus caerulein or alloxan alone. At day 13, islets in both groups consisted primarily of α-cells and δ-cells (Figures 2b and f, respectively-compare with normal islet in Figure 2a). However, at day 13, there were consistently more β-cells in islets of the alloxan plus caerulein group (7.8% of total endocrine cells) compared with the alloxan alone group (1.3% of total endocrine cells; Figure 2f versus b, quantified in j). This trend continued to day 20, when 22% of the total endocrine cells in the alloxan plus caerulein group expressed insulin (Figure 2g versus c, quantified in j), but the number of β-cells then declined at day 27 (Figure 2h versus d, quantified in j) and 3.7% at day 38 (Figure 2i versus e, quantified in j). Consistent with the lineage tracing study, the decreased number of β-cells occurred concomitantly with an increase in the number of δ-cells (Figures 2g–i, quantified in 2k). The number of α-cells did not change significantly (Figure 2l). This was expected, given our previous finding that α-cells replicate following PDL plus alloxan.14

The initial rise in β-cell number followed by a decline (Figure 2j), occurring concomitantly with a late rise in the number of δ-cells (Figure 2k), was consistent with the genetic lineage tracing experiment (Figure 1i). Thus, both the lineage tracing and time course experiments support a model where transdifferentiation occurred by a process where α-cells first converted into β-cells, with δ-cell formation occurring later.

Intermediate cells predicted to exist as a consequence of α- to β- to δ-cell transdifferentiation are present following alloxan plus caerulein

A sequential process of α- to β- to δ-cell transdifferentiation should lead to intermediates with glucagon–insulin or insulin–somatostatin but not glucagon–somatostatin coexpression, whereas if δ-cells formed directly from α-cells, intermediate cell types coexpressing glucagon and somatostatin should be found. As predicted by the α- to β- to δ-cell model, glucagon–insulin coexpression was high at day 13, the earliest time-point studied, declined at day 27 and persisted until day 38 (Figure 2m). Insulin–somatostatin coexpressing cells rose later, peaking at day 27 (Figure 2n). Glucagon–somatostatin coexpression was never found at any time-point.

The sequential model predicts that transcription factors expressed in β-cells should be found in cells with glucagon–insulin and insulin–somatostatin coexpression. Thus, we examined the expression of the transcription factors pancreas–duodenum homeobox 1 (PDX1) and NK homeobox 6.1 (NKX6.1). Both of the transcription factors are expressed in all β-cells, with NKX6.1 being highly β-cell specific,20, 21 while PDX1 is also expressed in a low percentage of δ-cells.22 Neither is detected in α-cells. Following alloxan plus caerulein, there was a large increase in islet cells coexpressing PDX1 and glucagon, reaching a peak at day 13 when almost 25% of glucagon-positive cells also expressed PDX1 (Figure 3a, quantified in e). Such cells were not found in untreated (Figure 3a, Day 0) or alloxan-injected mice (Supplementary Figure 3A, quantified in Figure 3e).

Figure 3
figure 3

PDX1 and NKX6.1 expression in somatostatin- and glucagon-positive cells following alloxan plus caerulein. (a) Cells coexpressing PDX1 (white) and glucagon (red) were found transiently following alloxan plus caerulein, peaking at day 13 (quantified in e). Coexpression of glucagon and NKX6.1 (green) was not found. (b) β-catenin (green), which localizes primarily to the cell membrane, was used to outline cells to confirm coexpression of PDX1 (white) and glucagon (red) in the same cell. PDX1-glucagon coexpression frequently occurred in the absence of insulin expression (Supplementary Figure 4), suggesting that PDX1 induction was an early event in the process of α- to β-cell transdifferentiation. (c) Cells coexpressing NKX6.1 (green) and somatostatin (red) were found transiently following alloxan plus caerulein, peaking at day 20, but not after alloxan alone (quantified in g; see also Supplementary Figure 3B). Note that in many of the cells coexpressing NKX6.1 and somatostatin at day 20, NKX6.1 was present in the cytoplasm (additional examples in Supplementary Figure 5). (d) β-catenin (white) was used to outline cells to confirm the coexpression of NKX6.1 (green) and somatostatin (red). No such cells were found in islets from mice injected with alloxan alone (Supplementary Figure 3B). Further examination of the cells coexpressing NKX6.1 and somatostatin revealed that in cells expressing somatostatin and NKX6.1 but not insulin, the NKX6.1 was invariably cytoplasmic (area outlined in red in Supplementary Figure 5E).39 In cells that coexpressed NKX6.1, somatostatin and insulin, NKX6.1 was usually nuclear (Supplementary Figures 5C–F), but in rare cells was cytoplasmic. No cells coexpressing insulin and somatostatin but not NKX6.1 were found. Cytoplasmic expression of NKX6.1 has been described, with the suggestion that it is controlled by the redox state of the cell and that its localization may be important in Type II diabetes.39 PDX1 also exhibited a transient increase in expression in cells expressing somatostatin, peaking at day 20 (quantified in f), which did not occur following alloxan alone (Supplementary Figure 3) In alloxan-injected mice, expression of PDX1 in δ-cells remained constant over time at 9.5% (Supplementary Figure 3B, quantified in Figure 3f), but rose to more than 40% at day 20 in alloxan plus caerulein treated mice (Figure 3c, quantified in f). Scale bar=75 μm

Because of its specificity for β-cells in the adult,20, 21 NKX6.1 was informative in revealing the sequence by which transdifferentiation occurred. Unlike PDX1, and in contrast to our previous study,9 no glucagon and NKX6.1 coexpression was found. With regard to δ-cell neogenesis, if δ-cells arose by a pathway involving a β-cell intermediate, there should be cells coexpressing somatostatin and NKX6.1, whereas if δ-cells arose directly from α-cells, NKX6.1 should never have been expressed and there should be no cells coexpressing somatostatin and NKX6.1. At day 20, 10% of the cells expressing somatostatin coexpressed NKX6.1 (Figure 3c, quantified in g). These observations support a model where δ-cell transdifferentiation from α-cells obligatorily passes through a β-cell intermediate.

Cellular replication is a central aspect of stem cell biology.23 In an independent experiment, mice were administered BrdU continuously following alloxan plus caerulein to mark cells that had replicated. Both BrdU-positive and negative cells that coexpressed insulin and somatostatin coexpression were found at day 20 (Supplementary Figure 6). This is consistent with our previous study of islet cell transdifferentiation9 and indicates that transdifferentiation occurred independently of replication.

In the time course experiment, six mice were maintained for at least 38 days following alloxan plus caerulein, from which two exhibited reduced blood glucose over time, whereas the others remained diabetic (Supplementary Figure 7A). The mice with lower blood glucose were killed at days 38 and 62. In contrast to mice treated with alloxan alone, where α-cells were the predominant islet cell type at day 38 (Supplementary Figure 7B), and the alloxan plus caerulein mice that remained diabetic at day 38 and where δ-cells predominated (Supplementary Figure 7C), the alloxan plus caerulein mice that became normoglycemic had many β-cells in their islets (Supplementary Figure 7D). The islets had an intermixing of α-, β- and δ-cells (Supplementary Figures 7D–H).

Induction of islet cell transdifferentiation by alloxan plus caerulein was islet autonomous

Previously, we and others9, 24 found that profound β-cell ablation was required to stimulate endocrine cell transdifferentiation. Moderate β-cell ablation (up to 80%) led to β-cell replication rather than transdifferentiation.25 This suggested that α-cells sensed the number of β-cells and only when that fell below a certain threshold did transdifferentiation occur. The sensing mechanism could be a local phenomenon where α-cells sensed the proximate β-cell number, while sensing of the total β-cell mass was also a possibility.

Because the extent of β-cell ablation following alloxan injection varied between islets, we were able to measure the extent of transdifferentiation in islets that differed in the number of preexisting β-cells that survived alloxan. Using genetic lineage marking to distinguish between preexisting β-cells and neogenic β-cells, we counted the number of preexisting (mT expressing) and neogenic (mG expressing) β-cells in each islet. In islets with efficient β-cell ablation, there were many neogenic β-cells (Figures 4a–c), whereas in islets with many β-cells that survived alloxan, there was little to no transdifferentiation, that is, the β-cells expressed mT on their membranes (Figures 4d–f). There was a threshold effect, with the vast majority of transdifferentiated β-cells being found in islets with 10 or fewer preexisting β-cells that survived alloxan treatment (Figure 4g, P=0.0016). Only two mG-positive β-cells were found in islets with greater than ten preexisting β-cells. The threshold effect is consistent with a locally acting sensing mechanism where each islet behaves autonomously with respect to β-cell neogenesis from α-cells.

Figure 4
figure 4

Endocrine cell transdifferentiation is islet autonomous. Photomicrographs of pancreas sections from lineage tracing mice killed 38 days following alloxan plus caerulein. (a) An islet with no preexisting β-cells (blue cytoplasmic insulin and red mT membrane fluorescence), demonstrating neogenic β-cells (blue cytoplasmic insulin and green eGFP membrane fluorescence; high-power views in b and c, neogenic cells marked with cyan arrowheads). (d) An islet with many preexisting β-cells (high-power views in e and f, preexisting cells marked with yellow arrowheads) demonstrating no neogenic β- or δ-cells. (g) Quantification of neogenic β- and δ- cells per islet. The number of preexisting β-cells in each islet (insulin+, mT+) was plotted against the number of neogenic β- plus δ-cells (insulin or somatostatin+, eGFP+) in that islet. The P-value for the difference between islets in which there were greater than 10 preexisting β-cells and at least one eGFP-positive cell expressing insulin or somatostatin (quadrant 2) versus islets in which there were ten or fewer β-cells expressing mT and at least one eGFP-positive cell expressing insulin or somatostatin (quadrant 4) was 0.0016, calculated using a two sample proportion test with continuity correction in R. Scale bar for bf=10 μm

Endocrine cell transdifferentiation in the NOD mouse model of T1D

A logical next step was to examine T1D, the prototypical disease that combines inflammation with β-cell deficiency, for evidence of islet cell transdifferentiation. As genetic lineage tracing is not available in NOD mice, hormone colocalization, which we and others have shown to indicate islet cell transdifferentiation,9, 26, 27 was used as a marker. As in the alloxan plus caerulein model, there was an increase in δ-cells in the islets of NOD mice as they progressed from prediabetic to diabetic (Figure 5a, quantified in b and Supplementary Text).

Figure 5
figure 5

Endocrine cell transdifferentiation leads to an increase in δ-cells during the progression to diabetes in NOD mice. (a) A low-power view (scale bar=200 μm) of a prediabetic NOD mouse pancreas, with higher digital magnification (scale bars=75 μm) of four islets that represent different stages in the progression of autoimmune diabetes with Islet 1 having an intact core of β-cells in the islet core, Islet 2 exhibiting islet disorganization, Islet 3 exhibiting islet disorganization, β-cell loss and δ-cell increase and Islet 4 being composed almost entirely of δ-cells. Note that Islet 4 has no insulitis, typical of the late-stage diabetic NOD mouse. Confocal high-power images of islets 2 and 3 are shown in Supplementary Figure 9. (b) Quantification of the frequency of islet cell types in control, prediabetic and late-stage NOD mice. (c) δ/α ratio between control, prediabetic and late-stage NOD mice. Areas within the islet in which there is coexpression of multiple hormones, consistent with transdifferentiation, are shown in dg (insulin and glucagon colocalization) and hk (insulin and somatostatin colocalization). Scale bars for dk=20 μm. (l) Quantification of insulin and glucagon colocalization. (m) Quantification of insulin and somatostatin colocalization. Quantitative analysis of hormone colocalization was done using the ImageJ Mander’s coefficient algorithm on confocal images of islets. Confocal images of islets 2 and 3 are shown in Supplementary Figure 9. (n) Quantification in each islet of the relationship between cells coexpressing insulin and glucagon or insulin and somatostatin. Islets with hormone coexpression were plotted by the extent of coexpression and further subdivided by islets with (red squares) and without (blue diamonds) insulitis. Error bars are S.E.M. *P<0.05. ***P<0.001

Endocrine cell transdifferentiation in human T1D

To study transdifferentiation in human T1D (Supplementary Table 1), we first asked whether there was a δ-cell excess. As expected, islets were composed primarily of α-cells early in the course of disease (compare normal human pancreas in Figure 6a and Supplementary Figure 10 with early stage T1D pancreases in Figure 6b and Supplementary Figure 11). As the number of years of diabetes increased, so did the frequency of δ-cells (compare Figure 6b and Supplementary Figure 11 with Figure 6c and Supplementary Figure 12). The δ- to α-cell ratio demonstrated a statistically significant increase as a direct function of the number of years of diabetes, which did not occur in non-diabetic or type II diabetic patients (Figures 6e and f). In fact, the δ-/α-cell ratio decreased in type II diabetes with duration of diabetes (Figure 6d and Supplementary Figure 13, quantified in Figure 6f), most likely due to the α-cell hyperplasia that occurs.28

Figure 6
figure 6

Pancreatic δ-cells increase over time in patients with T1D. (a) pancreas from non-diabetic donor (nPOD 6140), demonstrating normal islet morphology and hormone expression. (b) pancreas from nPOD Donor #6052 (T1D for 1 year), demonstrating the α-cell predominance found in most of the pancreas of this patient. (c) Pancreas from nPOD Donor #6159 (T1D for 44 years), demonstrating δ-cell predominance. (d) Pancreas from a patient with type II diabetes (T2D) for 20 years, demonstrating mild loss of β-cells and no increase in δ-cells. For ad, scale bar=280 μm for low-power images and=75 μm for high-power images. (e) Quantification of the somatostatin-positive area normalized to the glucagon-positive area, demonstrating an increase over time with T1D. Red diamonds mark the age of type I diabetics and red squares mark the duration of T1D, connected by a black line for each patient. Blue circles mark the age of non-diabetic patients. Pearson’s coefficient for T1D was 0.61 and for non-diabetic patients was −0.74. (f) Quantification of the somatostatin-positive area normalized to the glucagon-positive area in patients with T1D greater than 20 years, compared with non-diabetic and type II diabetic patients, demonstrating a statistically significant increase in the number of δ-cells in the type I diabetic patients with years of disease. (gr) Islet hormone coexpression in a non-diabetic patient (NPOD Donor #6178) (gj), nPOD Donors #6052 and #6038, with T1D for 1 year (kn) and 20 years (or), respectively. Coexpression was analyzed by the ImageJ colocalization analysis highlighter algorithm that analyzes two channels (red and green) and marks colocalized areas in white. Insulin and glucagon (red and green in g, k and o), insulin and somatostatin (red and green in h, l and p), and somatostatin and glucagon (red and green in i, m and q). Note that all of the colors from the Image J analysis are pseudocolors. True colors from the same islets are shown in j, n and r. Examples of insulin and glucagon coexpression are shown in high power in s and t (double positive areas are indicated by a dashed white line), and of insulin and somatostatin coexpression in ux (double positive areas are indicated by a dashed yellow line). Quantification of insulin and glucagon coexpression is presented in y and for insulin and somatostatin in z in nPOD Donors #6052, #6196, #6081 and #6038 (Supplementary Table 1). These patients were selected because their pancreas had islets with a significant number of β-cells and only islets with greater than 20 β-cells were subjected to quantitative analysis. No colocalization between glucagon and somatostatin was found and so it could not be quantified. Coexpression was indistinguishable from background in the normal patient, but there was extensive insulin–glucagon coexpression in nPOD patients 6052, 6196 and 6081. Donors 6081 and 6038 had a high level of insulin and somatostatin coexpression. Quantification was performed on confocal images (Supplementary Figure 14). Error bars are S.E.M. of individual islets. Islets from five non-diabetic patients were pooled. Scale bar for it=75 μm

For direct evidence of transdifferentiation, we examined for intermediates in the process, finding cells with glucagon–insulin or insulin–somatostatin coexpression in some patients (Figures 6k–x, quantified in y and z, Supplementary Figures 11 and 14), but not in all (Figure 6c, Supplementary Figure 12). Cells coexpressing islet hormones were not found in non-diabetic controls (Figures 6g–j, quantified in y and z, Supplementary Figures 10 and 14) or patients with type II diabetes (Figure 6d, Supplementary Figure 13). No glucagon–somatostatin coexpressing cells were found in any patient. As in the mouse, we found islet cells coexpressing insulin, somatostatin and NKX6.1 in T1D but not normal pancreases (Supplementary Figure 15, compare the non-diabetic patient in A–D where NKX6.1 was restricted to β-cells, with the T1D patient in E–G where there were cells that coexpressed insulin, somatostatin and NKX6.1).

Although the fact that some patients with T1D exhibited glucagon–insulin coexpression indicated that α-cells were transdifferentiating into β-cells, the possibility of direct conversion of β-cells into δ-cells remained. A 37-year-old female with T1D for 20 years (nPOD donor #6038) who had a substantial number of residual β-cells was informative had a significant number of β-cells in most islets (Figures 6o–r, Supplementary Figures 11D, 14D and E), and 33% coexpressed somatostatin (Figure 6z). Islets were found where essentially all of the insulin-positive cells coexpressed somatostatin, with no glucagon–insulin coexpressing cells (Figures 6o–r and v–x, Supplementary Figure 14D and E and J–M). This is most consistent with direct transdifferentiation from β- to δ-cells, but examination of additional patients will be required to determine how general this phenomenon is.

Discussion

Here, we demonstrate the existence of a novel pathway of islet cell neogenesis where α-cells transdifferentiate first to β-cells, with a subset going on to transdifferentiate into δ-cells. Transdifferentiation was induced by a single peptide, caerulein, which was sufficient in the context of severe β-cell deficiency to promote α- to β- to δ-cell transdifferentiation. A study of β-cell regeneration following profound β-cell ablation also found evidence of α- to β-cell transdifferentiation, but the efficiency of conversion was low.24 Ectopic expression of the transcription factor Pax4 under the control of the glucagon promoter led to conversion of duct cells into β-cells, further validating the ability of cells expressing glucagon to convert into β-cells.26, 29 However, an important difference is that we have not observed evidence of induction of transdifferentiating cells in the ducts, either here or in our previous study using PDL plus alloxan.9

Although the genetic lineage tracing proved conclusively the existence of transdifferentiation from α- to β- and α- to δ-cells, the time course analysis verified that the transdifferentiation sequence was from α- to β- to δ-cells. The temporal appearance of δ-cells and glucagon–insulin and insulin–somatostatin coexpressing intermediates, as well as the absence of glucagon–somatostatin coexpressing cells, is consistent with the transdifferentiation occurring from α- to β- to δ-cells, without direct α- to δ-cell transdifferentiation (Figure 7).

Figure 7
figure 7

Model of pancreatic islet cell transdifferentiation. Intermediates in the transdifferentiation process are documented as follows: intermediate 2 is shown in Supplementary Figures 4C and D (surrounded by a white line), intermediate 3 in Supplementary Figure 4C (surrounded by a yellow line) and intermediate 5 in Supplementary Figure 5 and Figures 3c and d

Our data strongly indicate that hormone colocalization reflects transdifferentiation, as we never observed colocalization in normal mice. It was extremely rare even in the condition of pure β-cell ablation, studied previously by others.24 We never observed glucagon–somatostatin coexpressing cells, which might have occurred if the hormone colocalization was an artifact of immunostaining. Similar to us, other groups have found hormone colocalization only in settings where transdifferentiation occurred as demonstrated by genetic lineage tracing.26, 27 Consistent with our previous study of transdifferentiation following PDL plus alloxan,9 we never observed hormone colocalization in ducts, nor did we ever observe mG-positive cells in ducts.

β-cell dedifferentiation has been reported as occurring in response to deletion of FoxO1 and in type II diabetes.30, 31 The finding of a high level of glucagon–insulin colocalization in the setting of profound β-cell ablation argues against β-cell dedifferentiation occurring following alloxan plus caerulein. First, there were simply not sufficient residual β-cells following alloxan to account for the large number of cells coexpressing glucagon and insulin. Also, dedifferentiation resulted in the loss of β-cell identity, that is, β-cell markers such as insulin, PDX1 and NKX6.1, before redifferentiation along a different lineage.30 The glucagon–insulin and insulin–somatostatin coexpression, as well as somatostatin and NKX6.1, is not consistent with that phenomenon30 and argues for transdifferentiation as opposed to dedifferentiation and redifferentiation as the mechanism responsible for α- to β- to δ-cell conversion.

The δ-cell excess in murine32 and human T1D, with precisely the same intermediates as found in the alloxan plus caerulein model, demonstrates the relevance of endocrine cell transdifferentiation to autoimmune diabetes. However, endocrine cell transdifferentiation following alloxan plus caerulein occurred only in islets where there were few β-cells, whereas in T1D, transdifferentiation occurred in islets with a substantial number of residual β-cells, most likely reflecting a difference in the stimulus for transdifferentiation in autoimmunity versus that induced by caerulein. The simplest model that is consistent with the data presented here, postulating the existence of two factors, one positive and one repressive, is shown in Supplementary Figure 16.

In some human patients with T1D, islets were found with numerous β-cells, all of which coexpressed somatostatin, and with no evidence of α- to β-cell transdifferentiation, suggesting that direct β- to δ-cell transdifferentiation was taking place without an α-cell origin. Such islets were fairly normal in morphology, with a large number of β-cells located more centrally, all of which coexpressed somatostatin but with no glucagon–insulin colocalization. Analysis of additional patient samples, which are rare, will be of value in determining the degree to which this occurs in the general population of type I diabetics. Of importance, islet cell transdifferentiation was found to occur regardless of age or duration of diabetes, suggesting that stimulating this process might be a viable approach to therapy in individuals with longstanding T1D. This is in contrast to β-cell replication, the other major approach to β-cell regeneration, which declines markedly with age,33, 34, 35 but is consistent with our previous finding in the PDL plus alloxan model, where α- to β-cell transdifferentiation occurred in old mice.9

The dramatic δ-cell excess in murine and human T1D, combined with the observation that β-cells appear to be able to transdifferentiate directly into δ-cells, indicates that the conventional paradigm of T1D where β-cells are destroyed by autoimmunity is incomplete. Rather, in addition to autoimmune destruction of β-cells, the view of islet cell kinetics in T1D should be modified to take into account the endocrine cell transdifferentiation from α- to β-cells and β- to δ-cells. If the process of β- to δ-cell transdifferentiation is reversible, that may provide another route to β-cell regeneration in T1D. Regardless, the finding of efficient α- to β-cell transdifferentiation in T1D, as well as the finding that a single peptide is sufficient to drive transdifferentiation, raises the prospect that β-cell regeneration by the mobilization of α- and possibly δ-cells functioning as facultative β-cell stem cells may be feasible.

Materials and Methods

Mice

Glu-mT/mG mice for α-cell lineage tracing were a generous gift from Dr. Mark Huising (The Salk Institute for Biological Studies, La Jolla, California (USA)). Gt(ROSA)26Sortm4(ACTB-tdTomato,-EGFP)Luo/J (The Jackson laboratory (Bar Harbor, Maine)–stock # 007576)36 were crossed with Glucagon–Cre mice, expressing cre recombinase under the control of the glucagon promoter.37 This transgene has been well characterized in multiple studies, exhibiting very low to undetectable background recombination.37, 38 F1 mice were genotyped and mice positive for both markers were selected for the study. The F1 mice (Glu-mT/mG) have the mT cassette deleted in cre-expressing cells, while a membrane-targeted eGFP (mG) cassette is activated. Thus, α-cells and cells derived from α-cells that underwent cre-mediated recombination have a green fluorescent membrane, whereas preexisting cells not arising from recombined α-cells have a red fluorescent membrane. To calculate the recombination efficiency in α-cells, 444 cells from 12 mice were counted.

C57BL/6 mice were purchased from Harlan Sprague Dawley, Inc. (Placentia, CA, USA). NOD mice were bred at the Sanford-Burnham Medical Research Institute, La Jolla, California. All animal experiments were approved by the Institutional Animal Care and Use Committee of the Sanford-Burnham Medical Research Institute in accordance with national regulations.

Alloxan plus caerulein

Mice were injected intravenously with 110 mg/kg alloxan (Sigma-Aldrich, St. Louis, MO, USA) in PBS (phosphate-buffered saline) at day 1. Mice with a blood glucose level greater than 500 mg/dl (OneTouch Ultra Mini, Johnson & Johnson, Milpitas, CA, USA) were randomized into two groups. One group was injected with caerulein (Sigma-Aldrich; six IP injections on days 1, 3 and 5). Lineage tracing studies were done with a caerulein dose of 7 mg/kg and the time course study with C57BL/6 mice used a dose of 10 mg/kg in PBS. Mice (alloxan alone and alloxan plus caerulein) were killed on days:

Day 13- Three C57BL/6 alloxan plus caerulein mice, three C57BL/6 alloxan alone mice.

Day 20- Three lineage tracing alloxan plus caerulein mice, three lineage tracing alloxan alone mice, seven C57BL/6 alloxan plus caerulein mice, seven C57BL/6 alloxan alone mice. Of the seven C57BL/6 alloxan plus caerulein mice and seven alloxan alone mice, one group from each (n=3) had BrdU (1 mg/ml, Sigma-Aldrich) added to the drinking water, whereas the second (n=4) had no BrdU added.

Day 27- Three C57BL/6 alloxan plus caerulein mice, three C57BL/6 alloxan alone mice.

Day 38- Three lineage tracing alloxan plus caerulein mice, three lineage tracing alloxan alone mice, four C57BL/6 alloxan plus caerulein and four C57BL/6 alloxan alone.

NOD mouse experiments

Twenty 8.5 month-old non-diabetic NOD mice (blood glucose levels 80–200 mg/dl) and five diabetic NOD mice (age 6–15 months-old) were killed. Of the 20 non-diabetic NOD mice, six had an islet morphology typical of early stages of diabetes–that is, some islets with insulitis, disorganized islet morphology and transdifferentiating cells. Fourteen non-diabetic mice had no evidence of islet autoimmunity and were not analyzed further.

Human pancreatic tissue

Human pancreatic tissue was obtained from the Network for Pancreatic Organ Donors with Diabetes (nPOD), a collaborative T1D research project sponsored by JDRF. Organ Procurement Organizations (OPO) partnering with nPOD to provide research resources are listed at http://www.jdrfnpod.org/our-partners.php. Use of that tissue was approved by the SBMRI IRB (IRB-2013-019-13). All images from nPOD samples are available from the nPOD at http://www.jdrfnpod.org/online-pathology.php. A summary of the human patients is in Supplementary Table 1.

Immunohistochemical staining

Tissue was fixed in 4% formaldehyde over night at 4 °C (USB; OH, USA), washed in PBS, followed by overnight in 30% sucrose at 4 °C, then embedded in OCT (Optimal Cutting Temperature) (Sakura Finetek; CA, USA) and frozen at 80 °C. Slides of 5-μm thickness were washed four times with PBS and treated with 0.3% Triton in PBS for 10 min. Antigen retrieval (when applicable) was done using with CitriSolvTM (Fisher Scientific, PA, USA) for 7 or 10 min in sub-boiling temperature. After washing with PBS for 10 min, slides were incubated in blocking solution with 5% normal donkey serum (Jackson Immuno Research, PA, USA) for 60 min at room temperature. Human slides were fixed in 4% formaldehyde for 20 min at 4 °C, whereas all cryosections of 5-μm thickness were incubated with antisera specific for insulin (1/200, guinea pig, USBIO, Swampscott, MA, USA), insulin (1/200, rabbit, Santa Cruz Biotechnology, Santa Cruz, CA, USA), glucagon (1/2000, mouse, Sigma-Aldrich), somatostatin (1/200, goat, Santa Cruz Biotechnology), PDX1 (1/2000, goat, Abcam, Cambridge, MA, USA), NKX6.1 (1/1000, mouse, Beta Cell Biology Consortium, Nashville, TN, USA), 5-bromodeoxyuridine (BrdU; mouse, GE healthcare, Piscataway, NJ, USA). Secondary antibodies for detection of guinea pig, rabbit, goat or mouse antibodies were labeled with: Alexa Fluor 488 (Invitrogen, Carlsbad, CA, USA), Alexa Fluor 649 (Jackson ImmunoResearch Laboratories, West Grove, PA, USA), Rhodamine Red (Jackson ImmunoResearch Laboratories), Cy5 (Jackson ImmunoResearch Laboratories) and Dylight 405 (Jackson ImmunoResearch Laboratories). Nuclei were visualized with DAPI (4',6-diamidino-2- phenylindole; Sigma-Aldrich).

Imaging and image analysis

All slides were scanned at a magnification of × 20 using the Aperio Scanscope FL system (Aperio Technologies Inc., Vista, CA, USA). The appropriate dyes were assigned and illumination levels were calibrated using a preset procedure, the parameters were saved and applied to all slides. The acquired digital images represent whole tissue sections. Each section was 6-μm thick and calculations per islet reflect a section through that islet. Sections were evaluated for image quality. All acquired images were subsequently placed in dedicated project folders, and stored on a designated local server. Slides were viewed and analyzed using the web-based Image Scope viewer. Slides were quantified using the ‘Area Quantification FL’ algorithm (version 11 Aperio Technologies Inc.). The algorithm was optimized using a preset procedure and the subsequent macro was saved and applied to all slides.

Confocal images of the lineage experiment were acquired using Zeiss LSM710 or 780 confocal microscopes (Zeiss, Jena, Germany) at the Waitt Advanced Biophotonics Center Core Facility of the Salk Institute for Biological Sciences. Other confocal images were acquired using Olympus Fluoview FV1000 (Olympus Optical, Tokyo, Japan) at the SBMRI Imaging Core Facility.

Images of islet cells from the lineage tracing experiment were scored as expressing eGFP or mT on their plasma membranes. Only cells that had continuous, unbroken staining for the membrane markers were included in the analysis of hormone expression. Quantification and determination of double- and triple-positive cells was done using the Aperio Scanscope FL system as described above. Four slides that were 100 μM apart were analyzed from each C57BL/6 or NOD mouse. Confocal images were analyzed with ImageJ (NIH, Bethesda, MD, USA) for visualized colocalization (colocalization highlighter algorithm) and quantified (Mander’s coefficient algorithm). Human samples were scanned (Aperio Scanscope FL system) and analyzed as described above. A complete set of the confocal images used for this analysis is in Supplementary Movie 1.

Statistics

Because an a priori hypothesis was being tested, P-values were calculated as one-tailed homoscedastic Student’s t-tests except in Figure 4g, where a two-sample proportion test was used.