Your new experience awaits. Try the new design now and help us make it even better

HYPOTHESIS AND THEORY article

Front. Endocrinol., 16 September 2025

Sec. Diabetes: Molecular Mechanisms

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1676627

A mathematical model of insulin action on acinar tissue linking histology and radiological imaging of the pancreas

  • 1Department of Diagnostic Medicine, Dell Medical School, University of Texas at Austin, Austin, TX, United States
  • 2Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, United States

The pancreas is smaller in individuals with type 1 and type 2 diabetes. The etiology of this reduced pancreatic volume is not fully understood, but it may be due to loss of insulin’s trophic influence on exocrine pancreatic tissue. Supporting this, histological studies have identified a zone of acinar cell hypertrophy and hyperplasia surrounding pancreatic islets, putatively due to insulin action on peri-islet acinar tissue. This study develops a mathematical model of pancreas size incorporating beta cell density, beta cell clustering, and the magnitude and spatial extent of acinar cell expansion to estimate the relationship between beta cell mass and pancreas size. This model indicates that growth of acinar tissue surrounding the beta cell is sufficient to account for the smaller pancreas volume observed in individuals with diabetes. Furthermore, single beta cells and smaller beta cell clusters have a greater influence on pancreas size on a per cell basis, as larger islets have greater overlap in the zone of insulin action. Thus, changes in pancreas volume may be more sensitive to loss of single beta cells or small islets than larger islets. The model provides a conceptual framework linking histological and radiological imaging to better understand the relationship between pancreas volume and beta cell mass.

Introduction

Pancreas size is smaller in individuals with type 1 (T1D) and, to a lesser extent, type 2 diabetes (T2D) (1). Moreover, pancreas size is dynamic, declining as individuals progress through presymptomatic stages of T1D (2) and increasing in individuals with T2D remission (3). Pancreas size also changes over the lifespan with growth during childhood and adolescence and decline in older age (4). Islet transplantation studies indicate that the size of the pancreas size is proportional to the number of isolated islets (5), although other studies have found mixed results (6). Collectively, these results suggest that pancreas size correlates with beta cell mass and/or function, although the precise nature of this correlation remains unknown. Furthermore, the extent of pancreas size reduction in T1D (30–40% smaller) (1) surpasses the small fraction of the pancreas made up of beta cells (1-2%) (7), implicating loss of acinar tissue.

While a smaller pancreas size in individuals with T1D was observed 8 decades prior (8), the implications of this smaller pancreas size have not been well characterized. Furthermore, unlike the pancreatic beta cells, which pose significant challenges for imaging due to their small size and scattered distribution throughout the pancreas, pancreas size can be readily assessed with magnetic resonance imaging (MRI) or computed tomography (CT). Radiological imaging of the pancreas may provide information about diabetes progression or therapeutic response not captured by blood tests (2). However, the relationship between pancreas size and beta cell mass are not well understood. For example, it is not known whether there is a linear relationship between beta cell mass and pancreas size or whether there is a threshold of beta cell loss needed to detect changes in pancreas size.

Two mechanisms have been postulated to link beta cell loss and acinar cell atrophy: loss of trophic effects of insulin and autoimmune destruction of acinar tissue. Supporting the former, pancreas size in individuals with insulin deficiency, but lacking autoimmunity, have similar pancreas size as T1D (9). Additionally, acinar cell hypertrophy is observed in peri-islet regions of db/db mice (10), while chemical destruction of the beta cell eliminates peri-islet acinar hypertrophy (11). Human acinar tissue exhibits limited hypertrophy surrounding islets, but increased proliferation, suggesting that hyperplasia may be more dominant in humans (12). Together, these findings suggest absence of insulin may be a primary driver of pancreas atrophy in diabetes.

There are no studies that connect structural insights from pancreas autopsy studies with imaging observations of the pancreas in individuals with diabetes. Imaging of the entire intact excised human pancreas can offer a critical bridge between high-resolution histology and whole-organ radiological imaging (13). This study builds a mathematical model of pancreas size incorporating beta cell density and clustering estimates from excised pancreas studies to evaluate the effect of acinar tissue expansion on radiological imaging. This model can be used to guide interpretation of pancreas size measures and estimate relative loss of beta cell mass in individuals with diabetes.

Research design and methods

All statistical analyses and mathematical modeling were performed using R (version 4.3.1) within the RStudio environment (Posit, Boston, MA). Custom scripts were developed to implement the modeling framework, and standard packages were used for data processing, visualization, and statistical inference.

Geometric representation of beta cell distribution and pancreas volume

We modeled each beta cell as a sphere with radius of 5 µm (diameter of 10 µm), in agreement with prior histology studies (14, 15). Equation 1 was used to calculate the volume of each beta cell:

Volumebeta Cell=43 πrbeta cell3(1)

This yields a volume for a single beta cell of 524 µm3, in agreement with prior reports (16).

Islets were modeled as clusters of a range of number of beta cells, reflecting growing realization that islet sizes are more diverse than previously appreciated as our tools for imaging these small cell clusters improves (13). We modeled clusters consisting of a single beta cell, 5 cells, 25 cells, 100 cells, 250 cells, and 1000 cells (17). Each cell cluster was modeled as a perfect sphere with the same volume as the sum of the individual cells packed together. This assumes no extracellular space of packing inefficiencies. In addition to analyses of single-sized cell clusters, we modeled a distribution of beta cell cluster sizes throughout the pancreas, approximating findings from autopsy studies (13).

For beta cell density, we modeled a maximum beta cell volume comprising 3% of the pancreas based on a recent study (13) updating prior estimates of 1-2% beta cell density. We assumed that individuals with longstanding T1D have a complete loss of beta cells, in accordance with a prior report of 99% beta cell loss (18), although we note other studies have found a small population of beta cells remaining in T1D (19). Accordingly, we note that all results display beta cell densities ranging from 0-3%, and thus can be used to estimate changes in beta cell mass spanning this range. Pancreas size is influenced by body size, age, and sex, leading to a range of differences in pancreas size between non-diabetic controls and individuals with T1D (1). For our model, we assumed a non-diabetic pancreas volume of 90ml and a pancreas volume of 60ml in individuals devoid of beta cells, derived by averaging published results from a meta-analysis reporting pancreas volume in controls and individuals with longstanding T1D (1).

Modeling insulin action

We modeled insulin action as a concentric sphere surrounding the beta cell (or cluster of beta cells) as displayed in Figure 1. In the first approach, termed the constant shell model, we assumed that the trophic effect of insulin acted homogenously within a sphere surrounding the beta cell cluster (Figure 1). The equation describing this shell is given in Equation 2:

Figure 1
Cross-section diagrams on the left depict concentric colored shells representing layers with a central cluster. A graph on the right compares relative insulin concentration against the distance from the cluster edge for different shell thicknesses (1, 5, 10, 20, and 30 micrometers) using colors matching the diagrams. Solid lines indicate the shell model, and dashed lines indicate the decay model.

Figure 1. Left) Cartoon depicting the influence of a beta cell on a concentric shell of acinar tissue surrounding a single beta cell (top) or cluster of beta cells (bottom). A shell of thickness 10 µm is shown for illustrative purposes. Right) Graph of insulin action on surrounding acinar tissue as a function of distance from the cell cluster edge. Colors display various distances from the beta cell cluster upon which insulin drives acinar cell growth and correspond with the colors in the cartoon. Two models of insulin action are displayed. The constant shell model (solid line) assumes insulin acts homogenously within a shell surrounding the beta cell cluster of prescribed thickness. The insulin gradient model (dashed line) assumes insulin action is strongest at the outer edge of the beta cell cluster and decays exponentially with increasing distance from the beta cell cluster. The decay length parameter models how quickly the concentration of insulin decreases as you move away from a beta cell cluster (lower means faster decay). Cartoon was created using https://BioRender.com.

VolumeConstant Shell=43 πrshell343 πrbeta cell3(2)

The thickness of this shell, rshell, was modeled as 1, 5, 10, 20, and 30 µm beyond the outer edge of the beta cell cluster, building upon a recent report using a 33 µm width for peri-islet insulin action (20).

In a second approach, termed the insulin gradient model, we assumed that insulin concentration was highest at the outer edge of the beta cell cluster and decayed exponentially with increasing distance from the beta cell cluster. This model includes a decay length parameter (λ) that represents how quickly the concentration of insulin decreases as you move away from a beta cell cluster. The volume of this insulin gradient model is given in Equation 3:

VolumeInsulin Gradient=rbeta cellrbeta cell+5λ4 πr2·e(rrbeta cell)/λdr(3)

The decay length, λ, was set at 1, 5, 10, 20, and 30 µm to determine the distance over which insulin has a significant effect on the surrounding tissue. A small decay length means insulin’s influence drops off rapidly and is limited to areas very close to the cluster, while a larger value means insulin spreads out further, affecting a wider region around the cluster. A plot of the effect of insulin concentration as a function of decay length and distance from the beta cell cluster edge is shown in Figure 1. For these initial analyses, we assumed that the acinar tissue volume doubled within the volume surrounding the beta cell cluster.

To simulate varying degrees of acinar cell expansion around beta cells, we adjusted the magnitude of acinar cell growth within the shell surrounding the beta cell cluster. The maximum expansion was set to 100%, corresponding to the doubling of acinar cell volume in the affected volume used in prior analyses. Expansion was also modeled as 50% and 20% increases to simulate smaller amounts of acinar tissue growth in response to insulin. These varying levels of expansion allowed us to examine how different degrees of acinar cell growth might influence pancreas size.

Results

The relationship between beta cell mass and pancreas size depended on the spatial extent of insulin action as well as the beta cell cluster size. In the constant shell model, a wider thickness of affected acinar tissue increased the slope of the line linking beta cell density and pancreas volume (Figure 2). Smaller clusters of beta cells resulted in higher pancreas volumes for a given beta cell density (Figure 2). Note the different y-axis in panels of Figure 2. This relationship between cell clustering and pancreas size results from greater overlap in the shells surrounding larger clusters of beta cells versus smaller clusters of individual beta cells. For example, a single beta cell with an influence radius of 20 µm leads to an increased acinar tissue volume of 64,973 µm3. 100 of these single beta cells with no overlap would lead to an increased volume of 6,497,300 µm3. In contrast, a cluster of 100 beta cells would only lead to an increased volume of 287,431 µm3, 22 times smaller, due to overlap of the sphere of influence of these 100 clustered cells. The dashed lines on Figure 2 display a pancreas volume of 90ml, the estimated volume of a pancreas with a full complement of beta cells. For single beta cell clusters, shell thicknesses of 10, 20, or 30 µm exceed the 90ml pancreas volume at 3% beta cell density. Clusters of 250 or 500 beta cells do not reach 90ml at any shell thickness modeled.

Figure 2
Six line graphs show the relationship between beta cell percentage and pancreas volume for different cell clusters and shell thicknesses. Each graph represents a different cell per cluster value: 1, 5, 25, 100, 250, and 500. Lines are color-coded for shell thickness from 1 to 30 micrometers. Pancreas volume increases with beta cell percentage and shell thickness, with notable differences across cluster sizes.

Figure 2. In the constant shell model, pancreas volume correlates with beta cell density. The slope of this relationship is influenced by both the beta cell cluster size (shown in different panels) and the thickness of acinar tissue expansion surrounding each cell cluster (displayed according to color). Smaller beta cell clusters lead to larger pancreas volume due to diminished overlap in the shells surrounding each cluster. The pancreas volume of a non-diabetic control (90ml) is shown as a dashed line. This model assumes acinar tissue doubles within the affected volume. Note the different y-axis maximum in each panel.

In the insulin gradient model (Figure 3), pancreas volume is higher than the corresponding shell model for a given beta cell cluster size. The correlation between beta cell density and pancreas volume is influenced by both the beta cell cluster size and the decay length of acinar tissue expansion surrounding each cell cluster. Again, note the different y-axis in panels of Figure 3. Increasing decay length corresponds to a larger volume of acinar tissue expansion. For single beta cell clusters, all decay lengths can exceed the 90ml pancreas volume at 3% beta cell density. For beta cell clusters of 250 only a 30 µm decay length exceeds the 90ml pancreas volume at 3% beta cell density. Clusters of 500 beta cells cannot reach a pancreas volume of 90ml at any decay length modeled.

Figure 3
Six line graphs display pancreas volume versus beta cell percentage for different cluster sizes: 1, 5, 25, 100, 250, and 500 cells per cluster. Each graph includes lines representing decay lengths of 1, 5, 10, 20, and 30 micrometers in blue, green, orange, red, and purple, respectively. Pancreas volume increases with both beta cell percentage and decay length, with steeper slopes for higher decay lengths.

Figure 3. In the exponential decay model, pancreas volume correlates with beta cell density and is higher than the corresponding shell model. The slope of this relationship is influenced by both the beta cell cluster size (shown in different panels) and the decay length of acinar tissue expansion surrounding each cell cluster (displayed according to color). Increasing decay length corresponds to a larger volume of acinar tissue expansion, and more strongly influences the decay model than the constant shell model. The pancreas volume of a non-diabetic control (90ml) is shown as a dashed line. This model assumes acinar tissue doubles within the affected volume. Note the different y-axis maximum in each panel.

Given the strong influence of beta cell cluster size on pancreas volume, and the existence of variable islet sizes within the pancreas, we modeled further results as a distribution of beta cell cluster sizes. Specifically, we modeled a distribution of 5% single beta cells, 10% clusters of 10 cells, 20% clusters of 25 cells, 40% clusters of 100 cells, 20% clusters of 250 cells, and 5% clusters of 1000 cells, approximating findings from a recent autopsy study (17). For this analysis, we modeled different magnitudes of acinar tissue expansion within the affected region, ranging from a 20% to 100% increase in volume (a doubling of acinar tissue volume within the region). We found that 100% expansion was able to recapitulate the 90ml pancreas volume at 3% beta cell density in both the shell model and decay model (Figure 4). A 50% expansion of acinar cell volume reached a 90ml pancreas volume for the 30 µm shell thickness as well as 20 and 30 µm decay length. The lowest expansion magnitude, a 20% increase in acinar volume, only reached a pancreas volume of 90ml at 3% beta cell density in the decay model with the largest (30 µm) decay length.

Figure 4
Six line graphs compare pancreas volume against beta cell percentage under two models: Constant Shell Model and Insulin Gradient Model. Each model displays graphs for 20%, 50%, and 100% acinar expansion. Lines represent different shell thicknesses or decay lengths of 1, 5, 10, 20, and 30 micrometers, with varying slopes, showing volume increase as beta cell percentage increases.

Figure 4. A higher magnitude of acinar tissue expansion surrounding beta cell clusters leads to increased pancreas volume. Graphs illustrate increases of 20%, 50%, and 100% (a doubling of acinar tissue volume within the volume surrounding the beta cell cluster). The top row displays results for the shell model while the bottom row shows results for the exponential decay model. This model assumes a heterogeneous distribution of 5% single beta cells, 10% clusters of 10 cells, 20% clusters of 25 cells, 40% clusters of 100 cells, 20% clusters of 250 cells, and 5% clusters of 1000 cells. Note the different y-axis maximum in each panel.

Discussion

This study models the relationship between peri-islet acinar tissue expansion seen on histology with reduced pancreas size observed in the diabetic pancreas. Using previously reported parameters for beta cell density and acinar cell expansion, we demonstrate that loss of acinar tissue surrounding islets is sufficient to explain the smaller pancreas observed in individuals with diabetes. Indeed, the model recapitulates autopsy studies demonstrating declines in pancreas size, beta cell number, and acinar cell number found in T1D (19). We further establish the range of islet sizes and extent of acinar tissue expansion that correspond with clinically observed changes on radiological imaging. Importantly, this model suggests that loss of individual beta cells and small beta cell clusters leads to larger declines in pancreas size than loss of large islets.

The model developed herein suggests that smaller clusters of beta cells have a greater influence on pancreas volume than larger clusters, due to reduced overlap in the volume of insulin action surround each cell. This finding has important implications for understanding changes in pancreas volume given the heterogenous range of islet sizes in the pancreas (13). Specifically, the loss of small islets may may cause more pancreas atrophy than large islets, on a per cell basis. This finding carries added significance given a number of recent reports demonstrating a preferential loss of small beta cell clusters in individuals with T1D, particularly in the early stages of the disease (17, 2124). Our mathematical model suggests that the declines in pancreas size found using imaging may be inherently sensitive to this process. Specifically, our model predicts a linear relationship between pancreas size and the number of beta cells, but only for a given shell thickness/decay length and beta cell cluster size. As the slope of this linear relationship is dependent on these variables, temporal heterogeneity in beta cell loss, such as preferential loss of small endocrine clusters in early T1D, would lead to a non-linear relationship between beta cell mass and pancreas size over time. In this case, the pancreas volume would decline more quickly per beta cell lost in early T1D, with a slower decline in pancreas volume in later T1D as larger islets are lost. Beyond T1D, this model may also inform our understanding of other alterations in beta cell mass and/or pancreas size. Individuals with type 2 diabetes have a 50% deficit in beta cell mass (25) and a pancreas size intermediate between nondiabetic controls and T1D (1), in agreement with our model. In addition to understanding how declines in beta cell mass affect pancreas size, this model can simulate increases in beta cell mass/pancreas size, such as those observed in T2D remission (3).

Statistician George Box is credited with coining the aphorism ‘All models are wrong, but some are useful.’ This model is proposed in this spirit, as a first step toward modeling insulin action on acinar tissue, rather than a comprehensive model incorporating all aspects of insulin signaling and diabetes pathophysiology. As such, this mathematical model is limited by several simplifying assumptions. Our model assumed that individuals with T1D had complete loss of beta cell mass, but note that residual beta cells persist, even in longstanding T1D (26). This study modeled beta cells and their effect on surrounding acinar cells as spherical, but note that studies suggest that beta cells are elongated (14). Our model did not include other islet cells, either as drivers of acinar tissue expansion or contributors to cell cluster volume. For simplicity, we also assume the amount of insulin secreted by individual beta cells is independent of the size of the beta cell cluster. Importantly, insulin distribution is currently modeled as passive diffusion away from the beta cell; the model does not explicitly account for insulin consumption or binding to acinar cells (27). Finally, our model is unable to distinguish acinar cell hypertrophy from hyperplasia, the latter of which may be more dominant in humans (12).

Further studies are needed to both refine this model and validate it with histological studies. Model refinement includes incorporating the effects of other endocrine cells and the extracellular matrix, incorporating stochastic effects, and more sophisticated modeling of the insulin gradients surrounding islets. Importantly, our model assumes isotropic passive diffusion of insulin away from the beta cell. The vasculature surrounding the islet has complex structure that controls insulin delivery to acinar cells and modulates crosstalk between endocrine and exocrine cells (28). Future work may incorporate islet vasculature structure similar to our work modeling chemotherapy delivery to solid tumors (29). Our model is currently applied to one single point in time. However, it can be adapted to model temporal changes in beta cell mass during development of T1D and the effect of interventions that slow this process to better understand longitudinal imaging results.

In conclusion, the reduced pancreas size observed in individuals with T1D likely reflects not only beta cell loss but also acinar tissue atrophy in response to loss of insulin signaling. Radiological imaging can quantify the size of the pancreas, offering potential insights into disease progression or therapeutic response that may not be captured by blood-based measures alone. However, the relationship between pancreas size and beta cell mass remains incompletely understood. By integrating histological data on beta cell density and clustering from excised human pancreata, we developed a mathematical model that links structural changes in acinar and beta cell compartments to imaging-based measures of pancreas size. This framework provides a quantitative tool to interpret pancreas imaging, estimate relative beta cell loss, and bridge the gap between high-resolution autopsy studies and whole-organ radiological observations, advancing our understanding of pancreatic remodeling in diabetes.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.

Author contributions

JV: Formal Analysis, Methodology, Visualization, Resources, Data curation, Validation, Software, Investigation, Funding acquisition, Project administration, Writing – review & editing, Conceptualization, Writing – original draft.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. Financial support was graciously provided by the National Institutes of Health (DK129979 and HD115556) and Breakthrough T1D (1-INO-2023-1340-A-N). We gratefully acknowledge philanthropic support from Thomas J. and Karen K. Gentry. The study sponsors were not involved in the design of the study and did not impose any restrictions regarding the publication of the report. The authors have no relevant conflicts of interest to disclose.

Conflict of interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that Generative AI was used in the creation of this manuscript. During the course of preparing this work, the author used ChatGPT and Microsoft Copilot for the purposes of rewriting code and select sentences for improved clarity and readability. Following the use of this tool/service, the authors formally reviewed the content for its accuracy and edited it as necessary. The authors take full responsibility for all the content of this publication.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Garcia TS, Rech TH, and Leitao CB. Pancreatic size and fat content in diabetes: A systematic review and meta-analysis of imaging studies. PLoS One. (2017) 12:e0180911. doi: 10.1371/journal.pone.0180911

PubMed Abstract | Crossref Full Text | Google Scholar

2. Virostko J, Wright JJ, Williams JM, Hilmes MA, Triolo TM, Broncucia H, et al. Longitudinal assessment of pancreas volume by MRI predicts progression to stage 3 type 1 diabetes. Diabetes Care. (2024) 47:393–400. doi: 10.2337/dc23-1681

PubMed Abstract | Crossref Full Text | Google Scholar

3. Al-Mrabeh A, Hollingsworth KG, Shaw JAM, McConnachie A, Sattar N, Lean MEJ, et al. 2-year remission of type 2 diabetes and pancreas morphology: a post-hoc analysis of the DiRECT open-label, cluster-randomized trial. Lancet Diabetes Endocrinol. (2020) 8:939–48. doi: 10.1016/S2213-8587(20)30303-X

PubMed Abstract | Crossref Full Text | Google Scholar

4. Saisho Y, Butler AE, Meier JJ, Monchamp T, Allen-Auerbach M, Rizza RA, et al. Pancreas volumes in humans from birth to age one hundred taking into account sex, obesity, and presence of type-2 diabetes. Clin Anat. (2007) 20:933–42. doi: 10.1002/ca.20543

PubMed Abstract | Crossref Full Text | Google Scholar

5. Nanno Y, Wilhelm JJ, Heller D, Schat R, Freeman ML, Trikudanathan G, et al. Combination of pancreas volume and HbA1c level predicts islet yield in patients undergoing total pancreatectomy and islet autotransplantation. Clin Transplant. (2020) 34:e14008. doi: 10.1111/ctr.14008

PubMed Abstract | Crossref Full Text | Google Scholar

6. Olehnik SK, Fowler JL, Avramovich G, and Hara M. Quantitative analysis of intra- and inter-individual variability of human beta-cell mass. Sci Rep. (2017) 7:16398. doi: 10.1038/s41598-017-16300-w

PubMed Abstract | Crossref Full Text | Google Scholar

7. Scharfmann R, Staels W, and Albagli O. The supply chain of human pancreatic beta cell lines. J Clin Invest. (2019) 129:3511–20. doi: 10.1172/JCI129484

PubMed Abstract | Crossref Full Text | Google Scholar

8. Gepts W. Pathologic anatomy of the pancreas in juvenile diabetes mellitus. Diabetes. (1965) 14:619–33. doi: 10.2337/diab.14.10.619

PubMed Abstract | Crossref Full Text | Google Scholar

9. Wright JJ, Williams JM, Letourneau-Freiberg LR, Kandasamy B, Reyes D, Kanegusuku AG, et al. Insulin deficiency from insulin gene mutation leads to smaller pancreas. Diabetes Care. (2023) 46:773–6. doi: 10.2337/dc22-2082

PubMed Abstract | Crossref Full Text | Google Scholar

10. Egozi A, Bahar Halpern K, Farack L, Rotem H, and Itzkovitz S. Zonation of pancreatic acinar cells in diabetic mice. Cell Rep. (2020) 32:108043. doi: 10.1016/j.celrep.2020.108043

PubMed Abstract | Crossref Full Text | Google Scholar

11. Kramer MF and Tan HT. The peri-insular acini of the pancreas of the rat. Z Zellforsch Mikrosk Anat. (1968) 86:163–70. doi: 10.1007/BF00348522

PubMed Abstract | Crossref Full Text | Google Scholar

12. Anzi S, Stolovich-Rain M, Klochendler A, Fridlich O, Helman A, Paz-Sonnenfeld A, et al. Postnatal exocrine pancreas growth by cellular hypertrophy correlates with a shorter lifespan in mammals. Dev Cell. (2018) 45:726–37.e723. doi: 10.1016/j.devcel.2018.05.024

PubMed Abstract | Crossref Full Text | Google Scholar

13. Lehrstrand J, Davies WIL, Hahn M, Korsgren O, Alanentalo T, and Ahlgren U. Illuminating the complete ss-cell mass of the human pancreas- signifying a new view on the islets of Langerhans. Nat Commun. (2024) 15:3318. doi: 10.1038/s41467-024-47686-7

PubMed Abstract | Crossref Full Text | Google Scholar

14. Cottle L, Gan WJ, Gilroy I, Samra JS, Gill AJ, Loudovaris T, et al. Structural and functional polarization of human pancreatic beta cells in islets from organ donors with and without type 2 diabetes. Diabetologia. (2021) 64:618–29. doi: 10.1007/s00125-020-05345-8

PubMed Abstract | Crossref Full Text | Google Scholar

15. Zimny ML and Blackard WG. The surface structure of isolated pancreatic islet cells. Cell Tissue Res. (1975) 164:467–71. doi: 10.1007/BF00219938

PubMed Abstract | Crossref Full Text | Google Scholar

16. Sasaki H, Saisho Y, Inaishi J, Watanabe Y, Tsuchiya T, Makio M, et al. Reduced beta cell number rather than size is a major contributor to beta cell loss in type 2 diabetes. Diabetologia. (2021) 64:1816–21. doi: 10.1007/s00125-021-05467-7

PubMed Abstract | Crossref Full Text | Google Scholar

17. Murrall K, Luckett T, Lekka C, Flaxman CS, Wyatt R, Akhbari P, et al. Small things matter: Lack of extra-islet beta cells in Type 1 diabetes. bioRxiv. (2025). doi: 10.1101/2025.04.11.648319

Crossref Full Text | Google Scholar

18. Butler AE, Galasso R, Meier JJ, Basu R, Rizza RA, and Butler PC. Modestly increased beta cell apoptosis but no increased beta cell replication in recent-onset type 1 diabetic patients who died of diabetic ketoacidosis. Diabetologia. (2007) 50:2323–31. doi: 10.1007/s00125-007-0794-x

PubMed Abstract | Crossref Full Text | Google Scholar

19. Wright JJ, Saunders DC, Dai C, Poffenberger G, Cairns B, Serreze DV, et al. Decreased pancreatic acinar cell number in type 1 diabetes. Diabetologia. (2020) 63:1418–23. doi: 10.1007/s00125-020-05155-y

PubMed Abstract | Crossref Full Text | Google Scholar

20. Kattner N, Hang Y, Krentz NAJ, Russell LA, Palmer M, Flaxman C, et al. Identification of a vimentin-expressing alpha-cell phenotype in CF and normal pancreas. J Endocrinol. (2025) 264(3):e240190. doi: 10.1530/JOE-24-0190

PubMed Abstract | Crossref Full Text | Google Scholar

21. Drotar DM, Vazquez Ramos GJA, Williams MD, David ST, Luce C, Smith JA, et al. Loss of insulin-positive cell clusters precedes the decrease in islet frequency and beta-cell area in type 1 diabetes. Diabetes. (2025). doi: 10.2337/db25-0326

PubMed Abstract | Crossref Full Text | Google Scholar

22. Granlund L and Lundberg M. Loss of insulin-expressing extra-islet cells in type 1 diabetes is accompanied with increased number of glucagon-expressing extra-islet cells. Virchows Archiv. (2025) 486:687–95. doi: 10.1007/s00428-024-03842-4

PubMed Abstract | Crossref Full Text | Google Scholar

23. Rippa A, Posgai AL, Currlin S, Brusko M, Williams MD, Kaddis JS, et al. 3D imaging of human pancreas suggests islet size and endocrine composition influence their loss in type 1 diabetes. bioRxiv. (2025). doi: 10.21203/rs.3.rs-6759572/v1

PubMed Abstract | Crossref Full Text | Google Scholar

24. van der Heide V, McArdle S, Nelson MS, Cerosaletti K, Gnjatic S, Mikulski Z, et al. Integrated histopathology of the human pancreas throughout stages of type 1 diabetes progression. bioRxiv. (2025). doi: 10.1101/2025.03.18.644000

PubMed Abstract | Crossref Full Text | Google Scholar

25. Butler AE, Janson J, Bonner-Weir S, Ritzel R, Rizza RA, and Butler PC. Beta-cell deficit and increased beta-cell apoptosis in humans with type 2 diabetes. Diabetes. (2003) 52:102–10. doi: 10.2337/diabetes.52.1.102

PubMed Abstract | Crossref Full Text | Google Scholar

26. Oram RA, Sims EK, and Evans-Molina C. Beta cells in type 1 diabetes: mass and function; sleeping or dead? Diabetologia. (2019) 62:567–77. doi: 10.1007/s00125-019-4822-4

PubMed Abstract | Crossref Full Text | Google Scholar

27. Korc M, Sankaran H, Wong KY, Williams JA, and Goldfine ID. Insulin receptors in isolated mouse pancreatic acini. Biochem Biophys Res Commun. (1978) 84:293–9. doi: 10.1016/0006-291X(78)90169-9

PubMed Abstract | Crossref Full Text | Google Scholar

28. Bonner-Weir S and Orci L. New perspectives on the microvasculature of the islets of Langerhans in the rat. Diabetes. (1982) 31:883–9. doi: 10.2337/diab.31.10.883

PubMed Abstract | Crossref Full Text | Google Scholar

29. Lorenzo G, Jarrett AM, Meyer CT, DiCarlo JC, Virostko J, Quaranta V, et al. A global sensitivity analysis of a mechanistic model of neoadjuvant chemotherapy for triple negative breast cancer constrained by in vitro and in vivo imaging data. Eng Comput. (2024) 40:1469–99. doi: 10.1007/s00366-023-01873-0

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: diabetes, type 1 diabetes, volume, peri-insular, islets, endocrine, exocrine

Citation: Virostko J (2025) A mathematical model of insulin action on acinar tissue linking histology and radiological imaging of the pancreas. Front. Endocrinol. 16:1676627. doi: 10.3389/fendo.2025.1676627

Received: 30 July 2025; Accepted: 03 September 2025;
Published: 16 September 2025.

Edited by:

Manami Hara, The University of Chicago, United States

Reviewed by:

Dirk Homann, Icahn School of Medicine at Mount Sinai, United States
Jaycob Dalton Warfel, University of Tennessee at Martin, United States

Copyright © 2025 Virostko. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: John Virostko, amFjay52aXJvc3Rrb0BhdXN0aW4udXRleGFzLmVkdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.