Integrin αvβ5 heterodimer is a specific marker of human pancreatic beta cells

The identification of cell surface markers specific to pancreatic beta cells is important for both the study of islet biology and for investigating the pathophysiology of diseases in which this cell type is lost or damaged. Following analysis of publicly available RNAseq data, we identified specific integrin subunits, integrin αv and integrin β5, that were expressed in beta cells. This finding was further elaborated using immunofluorescence analysis of histological sections derived from donor human pancreas. Despite the broad expression of specific integrin subunits, we found that expression of integrin αvβ5 heterodimers was restricted to beta cells and that this complex persisted in islet remnants of some type 1 diabetic individuals from which insulin expression had been lost. This study identifies αvβ5 heterodimers as a novel cell surface marker of human pancreatic beta cells, a finding that will aid in the identification and characterisation of this important cell type.

Integrin αvβ5 heterodimer is a specific marker of human pancreatic beta cells Jacqueline V. Schiesser 1,2 , Thomas Loudovaris 3 , Helen E. Thomas 3,4 , Andrew G. Elefanty 1,2,5,6 & Edouard G. Stanley 1,2,5,6* The identification of cell surface markers specific to pancreatic beta cells is important for both the study of islet biology and for investigating the pathophysiology of diseases in which this cell type is lost or damaged. Following analysis of publicly available RNAseq data, we identified specific integrin subunits, integrin αv and integrin β5, that were expressed in beta cells. This finding was further elaborated using immunofluorescence analysis of histological sections derived from donor human pancreas. Despite the broad expression of specific integrin subunits, we found that expression of integrin αvβ5 heterodimers was restricted to beta cells and that this complex persisted in islet remnants of some type 1 diabetic individuals from which insulin expression had been lost. This study identifies αvβ5 heterodimers as a novel cell surface marker of human pancreatic beta cells, a finding that will aid in the identification and characterisation of this important cell type.

T1D Type 1 diabetes
There is a dearth of cell surface markers available that specifically identify pancreatic beta cells-instead, enrichment strategies based on intracellular zinc granules (for example TSQ 1 or Newport Green 2 or negative selection are currently used. Given that integrins are often used to distinguish and isolate specific cell types from a variety of sources 3,4 , we sought to determine whether this group of proteins could identify specific sub-populations within human islets. Integrins are a large family of cell adhesion receptors that bind extracellular matrix (ECM) proteins and cell surface ligands, regulating a wide variety of processes such as tissue morphogenesis, cell adhesion and cytoskeletal organisation. The receptors themselves are heterodimeric in structure, consisting of non-covalently bound α and β subunits. Each α subunit preferentially associates with specific β subunits, a characteristic that confers an added level of cell type specificity. The αv-integrin subfamily comprises of five members (αvβ1, αvβ3, αvβ5, αvβ6, αvβ8), all of which recognise ECM and other peptides via a Arg-Gly-Asp (RGD) motif 5 . In the pancreas, interaction between the ECM and islet cells via adhesion-mediated integrin signaling has been shown to play a crucial role in maintenance of islet architecture and beta-cell function 6 .
In adults, the αv subunit was found to be expressed in the majority of the islet cells in intact tissue, but this expression was substantially decreased in isolated human islets 7 . Studies performed using mid-gestation (18)(19)(20) week) fetal pancreas found that the integrin heterodimers αvβ3 and αvβ5 were expressed in ductal cells and migrating endocrine cells. It was also observed that this expression decreased as the endocrine cells formed larger, more organised islet structures 8 . Further work demonstrated adult pancreatic cell preparations expressed the αv integrin subunit as well as both the β1 and β5 but not the β3 subunit 9 . This study suggested that αvβ1 was required for beta cell migration in the fetal pancreas and that the αvβ5 heterodimer had a significant role in adhesion to vitronectin 9 -an ECM component that has been previously identified in the developing and adult pancreas 9 .
In the present study, we show that the integrin heterodimer αvβ5 is specifically formed in beta cells of the adult human pancreas, in both control and type 1 diabetic subjects. Interestingly, analysis of pancreatic sections Islet isolation. Healthy human pancreata were obtained with informed consent from next of kin, from heart-beating, brain-dead donors, with research approval from the Human Research Ethics Committee at St Vincent's Hospital, Melbourne. Human islets were purified by intraductal perfusion and digestion of the pancreases with collagenase AF-1 (SERVA/Nordmark, Germany) 10 followed by purification using Ficoll density gradients 11 . Purified islets were cultured in Miami Media 1A (Mediatech/Corning 98-021, USA) supplemented with 2.5% human serum albumin (Australian Red Cross, Melbourne, VIC, Australia), in a 37 °C, 5% CO 2 incubator.
Immunofluorescence staining. Paraffin sections of donor human pancreas were obtained from the Tom Mandel Islet Isolation Program (St Vincent's Hospital, Victoria). Paraffin was removed using xylene, samples were rehydrated, and antigen retrieved using 10 mM citrate buffer. Samples were blocked for 1 h at room temperature in staining buffer (10% foetal calf serum (FCS) (Sigma-Aldrich; 12003C) in PBS) and 0.1% Triton-X (Sigma-Aldrich; T9284), stained overnight with primary antibodies at 4 °C, stained for 1 h at room temperature with secondary antibodies, and stained with DAPI (Sigma-Aldrich; D9542) for 5 min. Antibody details are provided in Table 3. Samples were mounted using Fluoromount-G (Southern Biotech; 0100-01) and imaged using a LSM780 inverted confocal microscope (Zeiss). Image analysis was performed using ImageJ (version 1.0).
Flow cytometry and sorting. Isolated human islets obtained from the Tom Mandel Islet Isolation Program were digested by resuspension in Accutase (Sigma-Aldrich; A6964) solution for 15 min at 37 °C. Following trituration, cells were washed in PBS and stained with primary antibody in FACS buffer (2% FCS in PBS) for 30 min on ice. Cells were then washed twice with FACS buffer and stained with the appropriate secondary antibody for 30 min on ice. Antibody details are provided in Table 3. Cells were then washed twice with FACS buffer and then resuspended in 1 μg/ml propidium iodide (Sigma-Aldrich; P4864) to exclude dead cells, prior to cell sorting. Flow sorting was performed on an BD Influx (BD Biosciences) or BD FACSAria Fusion (BD Biosciences). Isolated cells were resuspended in RLY buffer (Bioline) as a prelude to preparation of RNA.
Bulk RNAseq analysis. RNA 13 and edgeR 14,15 packages within R. Briefly, the counts per million (CPM) value was calculated using the cpm() function in edgeR, and genes expressed at low levels (defined as a CPM value below 0.5 in any sample) were filtered out. The filtered matrix was then used to create a DGEList object in edgeR, and the object TMM normalized using the calcNormFactors() function. Highly variable genes were identified by estimating the variance of each gene across samples, then sorting genes according to variance value. Unsupervised hierarchical clustering of samples was then performed using the top 500 variable genes. Heatmaps were created for the top 500 variable genes using the heatmap2() function in the gplots package ( Supplementary Figs. 1, 2).
Single cell RNAseq analysis. Processed RNA sequencing data was downloaded from GEO (GSE114412).
The count matrix was filtered to remove mitochondrially-encoded genes, genes with less than 1000 UMIFM counts and cells with greater than 25% mitochondrial-DNA content. Variation in the total counts of individual cells was then removed by normalizing the sum of counts for each individual cell to 10,000. The normalized counts were then used for dimensionality reduction and clustering for each dataset, which was performed using the Seurat package within R 16 . Briefly, highly variable genes were identified using the FindVariableFeatures() function within Seurat. Principle components were then computed and clustering was performed using Louvain community detection in the space of the first 30 principle components. UMAP projections were then computed using the first 30 principle components. Differentially expressed genes within clusters were then computed using the FindAllMarkers() function within Seurat and cluster identity assigned using the genes identified in Supplementary Fig. 1.

Results
In order to search for cell surface markers associated with specific cell types within human islets, we interrogated a previously published single cell transcriptional profiling (scRNAseq) data set, representing islets isolated from 4 independent non-diabetic donors 17,18 . Visualisation of this data set using uniform manifold approximation and projection (UMAP) analysis placed cells within distinct clusters (Fig. 1a). We assigned nominal identities to these clusters based on the expression of a key set of signature genes ( Supplementary Fig. 3a,b). This process enabled the identification of clusters corresponding to endocrine (alpha, beta, delta, gamma), exocrine (acinar, ductal) and immune (macrophage, mast cell, T-Cell, B-Cell) cells, as well as other cell types (endothelial, stellate, pericyte). We examined the gene expression profile associated with each cell cluster with a view to identifying potential cell surface markers that could be used to quantify islet composition using live cell techniques, such as flow cytometry. We paid particular attention to members of the integrin family, some of which have been previously associated with pancreatic endocrine development 8 and functionality of adult beta cells 19 . This analysis demonstrated that the integrin subunit αv was expressed in the endocrine cells of the pancreas. Subsequent analysis of publicly available RNAseq data showed that of the potential heterodimerisation partners, islets expressed beta  Fig. 3c,d). As expression of β1 has already been examined in the adult islet, we focussed on the next most highly expressed beta subunit, β5. In order to examine how the expression pattern of αv and β5 proteins related to the pattern suggested by RNAseq analysis, we used FACS to isolate specific islet subpopulations. Flow cytometry analysis showed that the vast majority of islet cells expressed some level of both the αv and β5 subunits. Islets also contained a smaller population that exclusively expressed the αv subunit (Fig. 1b). Using our analysis of single cell sequencing data as a guide, we purified populations expressing both subunits from isolated human islets and subjected these to bulk www.nature.com/scientificreports/ RNAseq analysis. As a point of reference, we also isolated a population of cells that only expressed the αv subunit, predicted to contain haematopoietic cells (See Supplementary Fig. 4).
Mapping the 20 most highly expressed genes in each isolated cell fraction to the gene expression profiles associated with specific cell clusters identified by scRNAseq analysis enabled us to estimate the cellular composition of the isolated cell populations (Fig. 1c-e). These associations can be visualised by means of a heat map ( Fig. 1c and Supplementary Fig. 4). This analysis suggested that the αv lo β5 lo population was enriched for endocrine cell types including alpha, delta and gamma cells (Fig. 1c-e) whilst the αv hi β5 hi population was enriched for pericytes and stellate cells (Supplementary Fig. 4b). Interestingly, the αv + β5 neg reference population appeared to be enriched for myeloid cells, most likely mast cells (Supplementary Fig. 4a).
To complement our RNAseq and flow cytometry analyses, we performed immunofluorescence staining on human pancreas from both control and T1D donors (Fig. 2, Supplementary Fig. 5). Consistent with our flow cytometry analysis, αv was found to be expressed throughout the islets; a pattern that was apparent in tissue from both control and T1D donors (Fig. 2a). Similarly, immunofluorescence analysis showed that β5 was expressed in endocrine cells (Fig. 2b), in addition to being present within endothelial and pericyte populations ( Supplementary Fig. 6). Immunofluorescence labelling also enabled us to examine the distribution of αvβ5 heterodimer using an antibody that specifically detects this complex. This analysis suggested that expression of the αvβ5 heterodimer was restricted to beta cells (Fig. 2c), as defined by the expression of INSULIN. In order to quantify the co-localisation of the αvβ5 heterodimer with either INSULIN or GLUCAGON in the islets of a control donor, we performed co-localisation analysis using the coloc2 plugin within ImageJ. This enabled us to calculate the Pearson's correlation coefficient for the proteins in question and to provide a quantitative measure of co-expression. From this analysis, we found that the mean Pearson's correlation coefficient for INSULIN and αvβ5 heterodimer was 1.00194 (Range 0.9434-1.0646, n = 5) indicating that these two proteins are indeed co-expressed. In contrast, we found that the mean Pearson's correlation coefficient for GLUCAGON and αvβ5 heterodimer was 0.0001 (Range − 0.286 to 0.198, n = 4) affirming that there is no co-expression of these two proteins in the islets of control donors. To enable quantification of the relationship between insulin and αvβ5 expression within islets, we examined four fields of view at 5× magnification from sections representing three control donors and 3 T1D donors that had been stained for insulin, glucagon and αvβ5. We counted the number of islets (defined by clusters containing cells that expressed insulin and glucagon) and scored them for the coexpression of αvβ5. This analysis showed that all islets expressed αvβ5 irrespective of whether they were from a control or T1D donor. Furthermore, we also identified clusters of αvβ5 expressing cells and scored these for the presence of insulin. As expected, for control donors, this analysis showed that there were no examples of αvβ5 expressing cell clusters that did not also express insulin. Conversely, in a subset of sections derived from T1D donors, we found rare instances of αvβ5 cell clusters that lacked insulin expression ( Supplementary Fig. 7). In these cases, clusters of insulin negative αvβ5 positive cells also expressed glucagon, consistent with the possibility of trans-differentiation of beta cells within the islets of individuals with type 1 diabetes 20 (Donor T1D-2, Fig. 2c, see Supplementary Fig. 3c for additional donor samples). We also examined expression of other αv heterodimers-αvβ1 and αvβ6. However, using validated commercially available antibodies, we were unable to detect specific expression of either heterodimer within the pancreas for either control or type 1 diabetic donors ( Supplementary Fig. 8).

Discussion
During a search for novel beta cell surface markers, we identified integrin subunits, αv and β5, whose expression overlapped in beta cell enriched islet sub-fractions. Using immunofluorescence analysis, we found that αvβ5 heterodimers were specifically formed in beta cells. Surprisingly, specificity in heterodimer formation was not simply a function of the overlap in expression of its component parts, with both the αv and β5 subunits showing much broader expression across all endocrine cells within the islet. This observation raises the possibility of cell intrinsic factors capable of regulating heterodimerisation.
Although antibodies directed against the integrin subunits αv and β5 could be used to enrich for endocrine cells from disaggregated islets, commercially available antibodies recognising the αvβ5 heterodimer were unsuitable for isolation of live cells using FACS. It is unclear if this limitation of the antibody was due to the position of epitope which it recognised (intra-versus extra-cellular) or a consequence of epitope destruction caused by enzymatic disaggregation. Nevertheless, our observation that heterodimer formation is restricted to beta cells suggests that this cell type has a unique relationship with its extracellular environment, potentially impacting its function and biology.

Data availability
Bulk RNA-seq data used in this study has been deposited in the Gene Expression Omnibus (GEO) data base and are available under the accession number GSE167589. All other data supporting the findings of this study are available from the corresponding author upon reasonable request.