A genomic data archive from the Network for Pancreatic Organ donors with Diabetes.

The Network for Pancreatic Organ donors with Diabetes (nPOD) is the largest biorepository of human pancreata and associated immune organs from donors with type 1 diabetes (T1D), maturity-onset diabetes of the young (MODY), cystic fibrosis-related diabetes (CFRD), type 2 diabetes (T2D), gestational diabetes, islet autoantibody positivity (AAb+), and without diabetes. nPOD recovers, processes, analyzes, and distributes high-quality biospecimens, collected using optimized standard operating procedures, and associated de-identified data/metadata to researchers around the world. Herein describes the release of high-parameter genotyping data from this collection. 372 donors were genotyped using a custom precision medicine single nucleotide polymorphism (SNP) microarray. Data were technically validated using published algorithms to evaluate donor relatedness, ancestry, imputed HLA, and T1D genetic risk score. Additionally, 207 donors were assessed for rare known and novel coding region variants via whole exome sequencing (WES). These data are publicly-available to enable genotype-specific sample requests and the study of novel genotype:phenotype associations, aiding in the mission of nPOD to enhance understanding of diabetes pathogenesis to promote the development of novel therapies.

Milder loss of insulin-containing islets in individuals with type 1 diabetes and type 2 diabetes-associated TCF7L2 genetic variants.

AIMS/HYPOTHESIS: TCF7L2 variants are the strongest genetic risk factor for type 2 diabetes. In individuals with type 1 diabetes, these variants are associated with a higher C-peptide AUC, a lower glucose AUC during an OGTT, single autoantibody positivity near diagnosis, particularly in individuals older than 12 years of age, and a lower frequency of type 1 diabetes-associated HLA genotypes. Based on initial observations from clinical cohorts, we tested the hypothesis that type 2 diabetes-predisposing TCF7L2 genetic variants are associated with a higher percentage of residual insulin-containing cells (ICI%) in pancreases of donors with type 1 diabetes, by examining genomic data and pancreatic tissue samples from the Network for Pancreatic Organ donors with Diabetes (nPOD) programme. METHODS: We analysed nPOD donors with type 1 diabetes (n=110; mean±SD age at type 1 diabetes onset 12.2±7.9 years, mean±SD diabetes duration 15.3±13.7 years, 53% male, 80% non-Hispanic White, 12.7% African American, 7.3% Hispanic) using data pertaining to residual beta cell number; quantified islets containing insulin-positive beta cells in pancreatic tissue sections; and expressed these values as a percentage of the total number of islets from each donor (mean ± SD ICI% 9.8±21.5, range 0-92.2). RESULTS: Donors with a high ICI% (≥5) (n=30; 27%) vs a low ICI% (<5) (n=80; 73%) were older at onset (15.3±6.9 vs 11.1±8 years, p=0.013), had a shorter diabetes duration at donor tissue procurement (7.0±7.4 vs 18.5±14.3 years, p<0.001), a higher African ancestry score (0.2±0.3 vs 0.1±0.2, p=0.043) and a lower European ancestry score (0.7±0.3 vs 0.9±0.3, p=0.023). After adjustment for age of onset (p=0.105), diabetes duration (p<0.001), BMI z score (p=0.145), sex (p=0.351) and African American race (p=0.053), donors with the TCF7L2 rs7903146 T allele (TC or TT, 45.5%) were 2.93 times (95% CI 1.02, 8.47) more likely to have a high ICI% than those without it (CC) (p=0.047). CONCLUSIONS/INTERPRETATION: Overall, these data support the presence of a type 1 diabetes endotype associated with a genetic factor that predisposes to type 2 diabetes, with donors in this category exhibiting less severe beta cell loss. It is possible that in these individuals the disease pathogenesis may include mechanisms associated with type 2 diabetes and thus this may provide an explanation for the poor response to immunotherapies to prevent type 1 diabetes or its progression in a subset of individuals. If so, strategies that target both type 1 diabetes and type 2 diabetes-associated factors when they are present may increase the success of prevention and treatment in these individuals.

Exocrine and Endocrine Inflammation Increases Cellular Replication in the Pancreatic Duct Compartment in Type 1 Diabetes.

Context: We recently demonstrated increased cellular proliferation in the pancreatic ductal gland (PDG) compartment of organ donors with type 1 diabetes, suggesting that PDGs may harbor progenitor cells capable of pancreatic regeneration. Objective: We evaluated the impact of diabetes and pancreatic inflammation on PDG and interlobular duct (ILD) cellular proliferation and profiles. Methods: Endocrine hormone expression (insulin, glucagon, somatostatin, pancreatic polypeptide) and proliferating Ki67+ cells were localized within the PDG and ILD compartments by multicolor immunohistochemistry in cross-sections from the head, body, and tail regions of pancreata from those with (n = 31) or without type 1 diabetes (n = 43). Whole-slide scanned images were analyzed using digital pathology. Results: Type 1 diabetes donors with insulitis or histologically identified pancreatitis had increased cellular replication in the ILD and PDG compartments. Interestingly, while cellular proliferation within the pancreatic ductal tree was significantly increased in type 1 diabetes (PDG mean = 3.36%, SEM = 1.06; ILD mean = 2.78%, SEM = 0.97) vs nondiabetes(ND) subjects without pancreatic inflammation (PDG mean = 1.18%, SEM = 0.42; ILD mean = 0.74%, SEM = 0.15, P < 0.05), robust replication was also observed in ND donors with pancreatitis (PDG mean = 3.52%, SEM = 1.33; ILD mean = 2.18%, SEM = 0.54, P < 0.05). Few polyhormonal cells were present in the ILD (type 1 diabetes = 0.04 ± 0.02%; ND = 0.08 ± 0.03%, P = 0.40) or PDG compartment (type 1 diabetes = 0.02 ± 0.01%; ND = 0.08 ± 0.13%, P = 0.63). Conclusion: These data suggest that increased pancreatic ductal cell replication is associated with sustained pancreatic inflammation; however, as replicating cells were hormone-negative, PDGs do not appear to represent a compelling endogenous source of hormone-positive endocrine cells.

Limited extent and consequences of pancreatic SARS-CoV-2 infection.

Concerns that infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of coronavirus disease 2019 (COVID-19), may cause new-onset diabetes persist in an evolving research landscape, and precise risk assessment is hampered by, at times, conflicting evidence. Here, leveraging comprehensive single-cell analyses of in vitro SARS-CoV-2-infected human pancreatic islets, we demonstrate that productive infection is strictly dependent on the SARS-CoV-2 entry receptor ACE2 and targets practically all pancreatic cell types. Importantly, the infection remains highly circumscribed and largely non-cytopathic and, despite a high viral burden in infected subsets, promotes only modest cellular perturbations and inflammatory responses. Similar experimental outcomes are also observed after islet infection with endemic coronaviruses. Thus, the limits of pancreatic SARS-CoV-2 infection, even under in vitro conditions of enhanced virus exposure, challenge the proposition that in vivo targeting of ß cells by SARS-CoV-2 precipitates new-onset diabetes. Whether restricted pancreatic damage and immunological alterations accrued by COVID-19 increase cumulative diabetes risk, however, remains to be evaluated.

Single-cell analysis of the human pancreas in type 2 diabetes using multi-spectral imaging mass cytometry.

Type 2 diabetes mellitus (T2D) is a chronic age-related disorder characterized by hyperglycemia due to the failure of pancreatic beta cells to compensate for increased insulin demand. Despite decades of research, the pathogenic mechanisms underlying T2D remain poorly defined. Here, we use imaging mass cytometry (IMC) with a panel of 34 antibodies to simultaneously quantify markers of pancreatic exocrine, islet, and immune cells and stromal components. We analyze over 2 million cells from 16 pancreata obtained from donors with T2D and 13 pancreata from age-similar non-diabetic controls. In the T2D pancreata, we observe significant alterations in islet architecture, endocrine cell composition, and immune cell constituents. Thus, both HLA-DR-positive CD8 T cells and macrophages are enriched intra-islet in the T2D pancreas. These efforts demonstrate the utility of IMC for investigating complex events at the cellular level in order to provide insights into the pathophysiology of T2D.

Altered cellular localisation and expression, together with unconventional protein trafficking, of prion protein, PrPC, in type 1 diabetes.

AIMS/HYPOTHESIS: Normal cellular prion protein (PrPC) is a conserved mammalian glycoprotein found on the outer plasma membrane leaflet through a glycophosphatidylinositol anchor. Although PrPC is expressed by a wide range of tissues throughout the body, the complete repertoire of its functions has not been fully determined. The misfolded pathogenic isoform PrPSc (the scrapie form of PrP) is a causative agent of neurodegenerative prion diseases. The aim of this study is to evaluate PrPC localisation, expression and trafficking in pancreases from organ donors with and without type 1 diabetes and to infer PrPC function through studies on interacting protein partners. METHODS: In order to evaluate localisation and trafficking of PrPC in the human pancreas, 12 non-diabetic, 12 type 1 diabetic and 12 autoantibody-positive organ donor tissue samples were analysed using immunofluorescence analysis. Furthermore, total RNA was isolated from 29 non-diabetic, 29 type 1 diabetic and 24 autoantibody-positive donors to estimate PrPC expression in the human pancreas. Additionally, we performed PrPC-specific immunoblot analysis on total pancreatic protein from non-diabetic and type 1 diabetic organ donors to test whether changes in PrPC mRNA levels leads to a concomitant increase in PrPC protein levels in human pancreases. RESULTS: In non-diabetic and type 1 diabetic pancreases (the latter displaying both insulin-positive [INS(+)] and -negative [INS(-)] islets), we found PrPC in islets co-registering with beta cells in all INS(+) islets and, strikingly, unexpected activation of PrPC in alpha cells within diabetic INS(-) islets. We found PrPC localised to the plasma membrane and endoplasmic reticulum (ER) but not the Golgi, defining two cellular pools and an unconventional protein trafficking mechanism bypassing the Golgi. We demonstrate PrPC co-registration with established protein partners, neural cell adhesion molecule 1 (NCAM1) and stress-inducible phosphoprotein 1 (STI1; encoded by STIP1) on the plasma membrane and ER, respectively, linking PrPC function with cyto-protection, signalling, differentiation and morphogenesis. We demonstrate that both PRNP (encoding PrPC) and STIP1 gene expression are dramatically altered in type 1 diabetic and autoantibody-positive pancreases. CONCLUSIONS/INTERPRETATION: As the first study to address PrPC expression in non-diabetic and type 1 diabetic human pancreas, we provide new insights for PrPC in the pathogenesis of type 1 diabetes. We evaluated the cell-type specific expression of PrPC in the human pancreas and discovered possible connections with potential interacting proteins that we speculate might address mechanisms relevant to the role of PrPC in the human pancreas.

ACE2 chromogenic immunostaining protocol optimized for formalin-fixed paraffin-embedded human tissue sections.

Angiotensin-converting enzyme 2 (ACE2) is a key cellular entry factor for severe acute respiratory syndrome coronavirus 2. Hence, identifying cell types that express ACE2 is important for understanding the pathophysiology of coronavirus disease 2019. We performed extensive testing of multiple primary antibodies across various human tissue types. Here, we describe an optimized protocol for immunostaining of ACE2 in formalin-fixed paraffin-embedded human pancreas, small intestine, and kidney tissue sections obtained from organ donors and autopsies. For complete details on the use and execution of this protocol, please refer to Kusmartseva et al. (2020).

Observing Islet Function and Islet-Immune Cell Interactions in Live Pancreatic Tissue Slices.

Live pancreatic tissue slices allow for the study of islet physiology and function in the context of an intact islet microenvironment. Slices are prepared from live human and mouse pancreatic tissue embedded in agarose and cut using a vibratome. This method allows for the tissue to maintain viability and function in addition to preserving underlying pathologies such as type 1 (T1D) and type 2 diabetes (T2D). The slice method enables new directions in the study of the pancreas through the maintenance of the complex structures and various intercellular interactions that comprise the endocrine and exocrine tissues of the pancreas. This protocol demonstrates how to perform staining and time-lapse microscopy of live endogenous immune cells within pancreatic slices along with assessments of islet physiology. Further, this approach can be refined to discern immune cell populations specific for islet cell antigens using major histocompatibility complex-multimer reagents.

Image-Based Machine Learning Algorithms for Disease Characterization in the Human Type 1 Diabetes Pancreas.

Emerging data suggest that type 1 diabetes affects not only the ß-cell-containing islets of Langerhans, but also the surrounding exocrine compartment. Using digital pathology, machine learning algorithms were applied to high-resolution, whole-slide images of human pancreata to determine whether the tissue composition in individuals with or at risk for type 1 diabetes differs from those without diabetes. Transplant-grade pancreata from organ donors were evaluated from 16 nondiabetic autoantibody-negative controls, 8 nondiabetic autoantibody-positive subjects with increased type 1 diabetes risk, and 19 persons with type 1 diabetes (0 to 12 years' duration). HALO image analysis algorithms were implemented to compare architecture of the main pancreatic duct as well as cell size, density, and area of acinar, endocrine, ductal, and other nonendocrine, nonexocrine tissues. Type 1 diabetes was found to affect exocrine area, acinar cell density, and size, whereas the type of difference correlated with the presence or absence of insulin-positive cells remaining in the pancreas. These changes were not observed before disease onset, as indicated by modeling cross-sectional data from pancreata of autoantibody-positive subjects and those diagnosed with type 1 diabetes. These data provide novel insights into anatomic differences in type 1 diabetes pancreata and demonstrate that machine learning can be adapted for the evaluation of disease processes from cross-sectional data sets.

Pancreatlas: Applying an Adaptable Framework to Map the Human Pancreas in Health and Disease.

Human tissue phenotyping generates complex spatial information from numerous imaging modalities, yet images typically become static figures for publication, and original data and metadata are rarely available. While comprehensive image maps exist for some organs, most resources have limited support for multiplexed imaging or have non-intuitive user interfaces. Therefore, we built a Pancreatlas resource that integrates several technologies into a unique interface, allowing users to access richly annotated web pages, drill down to individual images, and deeply explore data online. The current version of Pancreatlas contains over 800 unique images acquired by whole-slide scanning, confocal microscopy, and imaging mass cytometry, and is available at https://www.pancreatlas.org. To create this human pancreas-specific biological imaging resource, we developed a React-based web application and Python-based application programming interface, collectively called Flexible Framework for Integrating and Navigating Data (FFIND), which can be adapted beyond Pancreatlas to meet countless imaging or other structured data-management needs.

Expression of SARS-CoV-2 Entry Factors in the Pancreas of Normal Organ Donors and Individuals with COVID-19.

Diabetes is associated with increased mortality from severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Given literature suggesting a potential association between SARS-CoV-2 infection and diabetes induction, we examined pancreatic expression of angiotensin-converting enzyme 2 (ACE2), the key entry factor for SARS-CoV-2 infection. Specifically, we analyzed five public scRNA-seq pancreas datasets and performed fluorescence in situ hybridization, western blotting, and immunolocalization for ACE2 with extensive reagent validation on normal human pancreatic tissues across the lifespan, as well as those from coronavirus disease 2019 (COVID-19) cases. These in silico and ex vivo analyses demonstrated prominent expression of ACE2 in pancreatic ductal epithelium and microvasculature, but we found rare endocrine cell expression at the mRNA level. Pancreata from individuals with COVID-19 demonstrated multiple thrombotic lesions with SARS-CoV-2 nucleocapsid protein expression that was primarily limited to ducts. These results suggest SARS-CoV-2 infection of pancreatic endocrine cells, via ACE2, is an unlikely central pathogenic feature of COVID-19-related diabetes.

Organisation of the human pancreas in health and in diabetes.

For much of the last century, our knowledge regarding the pancreas in type 1 and type 2 diabetes was largely derived from autopsy studies of individuals with these disorders or investigations utilising rodent models of either disease. While many important insights emanated from these efforts, the mode for investigation has increasingly seen change due to the availability of transplant-quality organ-donor tissues, improvements in pancreatic imaging, advances in metabolic assessments of living patients, genetic analyses, technological advances for laboratory investigation and more. As a result, many long-standing notions regarding the role for and the changes that occur in the pancreas in individuals with these disorders have come under question, while, at the same time, new issues (e.g., beta cell persistence, disease heterogeneity, exocrine contributions) have arisen. In this article, we will consider the vital role of the pancreas in human health and physiology, including discussion of its anatomical features and dual (exocrine and endocrine) functions. Specifically, we convey changes that occur in the pancreas of those with either type 1 or type 2 diabetes, with careful attention to the facets that may contribute to the pathogenesis of either disorder. Finally, we discuss the emerging unknowns with the belief that understanding the role of the pancreas in type 1 and type 2 diabetes will lead to improvements in disease diagnosis, understanding of disease heterogeneity and optimisation of treatments at a personalised level. Graphical abstract.

Removing Formaldehyde-Induced Peptidyl Crosslinks Enables Mass Spectrometry Imaging of Peptide Hormone Distributions from Formalin-Fixed Paraffin-Embedded Tissues.

Linking molecular and chemical changes to human disease states depends on the availability of appropriate clinical samples, mostly preserved as formalin-fixed paraffin-embedded (FFPE) specimens stored in tissue banks. Mass spectrometry imaging (MSI) enables the visualization of the spatiotemporal distribution of molecules in biological samples. However, MSI is not effective for imaging FFPE tissues because of the chemical modifications of analytes, including complex crosslinking between nucleophilic moieties. Here we used an MS-compatible inorganic nucleophile, hydroxylamine hydrochloride, to chemically reverse inter- and intra-crosslinks from endogenous molecules. The analyte restoration appears specific for formaldehyde-reactive amino acids. This approach enabled the MSI-assisted localization of pancreatic peptides expressed in the alpha, beta, and gamma cells. Pancreatic islet-like distributions of islet hormones were observed in human FFPE tissues preserved for more than five years, demonstrating that samples from biobanks can effectively be investigated with MSI.

Long-term culture of human pancreatic slices as a model to study real-time islet regeneration.

The culture of live pancreatic tissue slices is a powerful tool for the interrogation of physiology and pathology in an in vitro setting that retains near-intact cytoarchitecture. However, current culture conditions for human pancreatic slices (HPSs) have only been tested for short-term applications, which are not permissive for the long-term, longitudinal study of pancreatic endocrine regeneration. Using a culture system designed to mimic the physiological oxygenation of the pancreas, we demonstrate high viability and preserved endocrine and exocrine function in HPS for at least 10 days after sectioning. This extended lifespan allowed us to dynamically lineage trace and quantify the formation of insulin-producing cells in HPS from both non-diabetic and type 2 diabetic donors. This technology is expected to be of great impact for the conduct of real-time regeneration/developmental studies in the human pancreas.

Publisher Correction: Long-term culture of human pancreatic slices as a model to study real-time islet regeneration.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Large-scale electron microscopy database for human type 1 diabetes.

Autoimmune ß-cell destruction leads to type 1 diabetes, but the pathophysiological mechanisms remain unclear. To help address this void, we created an open-access online repository, unprecedented in its size, composed of large-scale electron microscopy images ('nanotomy') of human pancreas tissue obtained from the Network for Pancreatic Organ donors with Diabetes (nPOD; www.nanotomy.org). Nanotomy allows analyses of complete donor islets with up to macromolecular resolution. Anomalies we found in type 1 diabetes included (i) an increase of 'intermediate cells' containing granules resembling those of exocrine zymogen and endocrine hormone secreting cells; and (ii) elevated presence of innate immune cells. These are our first results of mining the database and support recent findings that suggest that type 1 diabetes includes abnormalities in the exocrine pancreas that may induce endocrine cellular stress as a trigger for autoimmunity.

Pancreas tissue slices from organ donors enable in situ analysis of type 1 diabetes pathogenesis.

In type 1 diabetes (T1D), autoimmune destruction of pancreatic ß cells leads to insulin deficiency and loss of glycemic control. However, knowledge about human pancreas pathophysiology in T1D remains incomplete. To address this limitation, we established a pancreas tissue slice platform of donor organs with and without diabetes, facilitating the first live cell studies of human pancreas in T1D pathogenesis to our knowledge. We show that pancreas tissue slices from organ donors allow thorough assessment of processes critical for disease development, including insulin secretion, ß cell physiology, endocrine cell morphology, and immune infiltration within the same donor organ. Using this approach, we compared detailed pathophysiological profiles for 4 pancreata from donors with T1D with 19 nondiabetic control donors. We demonstrate that ß cell loss, ß cell dysfunction, alterations of ß cell physiology, and islet infiltration contributed differently to individual cases of T1D, allowing insight into pathophysiology and heterogeneity of T1D pathogenesis. Thus, our study demonstrates that organ donor pancreas tissue slices represent a promising and potentially novel approach in the search for successful prevention and reversal strategies of T1D.

Temporal Analysis of Amylase Expression in Control, Autoantibody-Positive, and Type 1 Diabetes Pancreatic Tissues.

Within the human pancreas, exocrine and endocrine cells control secretion of digestive enzymes and production of hormones to maintain metabolic homeostasis, respectively. While the vast majority of type 1 diabetes research efforts have focused on endocrine function and autoimmunity, recent studies identified a series of unique features (e.g., reduced weight and volume, increased density of leukocytes) within the exocrine pancreas in this disease, but the mechanisms underlying these aberrancies are unknown. Therefore, we histologically assessed amylase, insulin, glucagon, lipase, and/or trypsinogen in 78 organ donor pancreata from birth through adulthood in control subjects and those at various stages of type 1 diabetes. While amylase-positive (AMY+) acinar cells were detectable in pancreata from all study groups, tissues from individuals >2 years of age contained clusters of acinar cells devoid of amylase (AMY-). A majority of these AMY- cell clusters localized proximal to islets (i.e., peri-islet). Additionally, most AMY- clusters were positive for the exocrine enzymes lipase and trypsinogen. Interestingly, type 1 diabetes pancreata displayed significant reductions in the frequency of these AMY- cell clusters. These results support a contribution of the islet-acinar axis in pancreatic development and underscore a potential role for the exocrine pancreas in the pathogenesis of type 1 diabetes.

A Map of Human Type 1 Diabetes Progression by Imaging Mass Cytometry.

Type 1 diabetes (T1D) results from the autoimmune destruction of insulin-producing ß cells. A comprehensive picture of the changes during T1D development is lacking due to limited sample availability, inability to sample longitudinally, and the paucity of technologies enabling comprehensive tissue profiling. Here, we analyzed 1,581 islets from 12 human donors, including eight with T1D, using imaging mass cytometry (IMC). IMC enabled simultaneous measurement of 35 biomarkers with single-cell and spatial resolution. We performed pseudotime analysis of islets through T1D progression from snapshot data to reconstruct the evolution of ß cell loss and insulitis. Our analyses revealed that ß cell destruction is preceded by a ß cell marker loss and by recruitment of cytotoxic and helper T cells. The approaches described herein demonstrate the value of IMC for improving our understanding of T1D pathogenesis, and our data lay the foundation for hypothesis generation and follow-on experiments.

Hospital time prior to death and pancreas histopathology: implications for future studies.

AIMS/HYPOTHESIS: Diabetes research studies routinely rely upon the use of tissue samples from human organ donors. It remains unclear whether the length of hospital stay prior to organ donation affects the presence of cells infiltrating the pancreas or the frequency of replicating beta cells. METHODS: To address this, 39 organ donors without diabetes were matched for age, sex, BMI and ethnicity in groups of three. Within each group, donors varied by length of hospital stay immediately prior to organ donation (<3 days, 3 to <6 days, or ≥6 days). Serial sections from tissue blocks in the pancreas head, body and tail regions were immunohistochemically double stained for insulin and CD45, CD68, or Ki67. Slides were electronically scanned and quantitatively analysed for cell positivity. RESULTS: No differences in CD45+, CD68+, insulin+, Ki67+ or Ki67+/insulin+ cell frequencies were found when donors were grouped according to duration of hospital stay. Likewise, no interactions were observed between hospitalisation group and pancreas region, age, or both; however, with Ki67 staining, cell frequencies were greater in the body vs the tail region of the pancreas (∆ 0.65 [unadjusted 95% CI 0.25, 1.04]; p = 0.002) from donors <12 year of age. Interestingly, frequencies were less in the body vs tail region of the pancreas for both CD45+ cells (∆ -0.91 [95% CI -1.71, -0.10]; p = 0.024) and insulin+ cells (∆ -0.72 [95% CI -1.10, -0.34]; p < 0.001). CONCLUSIONS/INTERPRETATION: This study suggests that immune or replicating beta cell frequencies are not affected by the length of hospital stay prior to donor death in pancreases used for research. DATA AVAILABILITY: All referenced macros (adopted and developed), calculations, programming code and numerical dataset files (including individual-level donor data) are freely available on GitHub through Zenodo at https://doi.org/10.5281/zenodo.1034422.