Original Research ARTICLE
Analysis of Purified Pancreatic Islet Beta and Alpha Cell Transcriptomes Reveals 11β-Hydroxysteroid Dehydrogenase (Hsd11b1) as a Novel Disallowed Gene
- 1Section of Cell Biology and Functional Genomics, Department of Medicine, Imperial College London, London, UK
- 2Department of Neurobiology, Physiology, and Behavior, College of Biological Sciences, University of California, Davis, Davis, CA, USA
We and others have previously identified a group of genes, dubbed “disallowed,” whose expression is markedly lower in pancreatic islets than in other mammalian cell types. Forced mis-expression of several members of this family leads to defective insulin secretion, demonstrating the likely importance of disallowance for normal beta cell function. Up to now, transcriptomic comparisons have been based solely on data from whole islets. This raises the possibilities that (a) there may be important differences in the degree of disallowance of family members between beta and other either neuroendocrine cells; (b) beta (or alpha) cell disallowed genes may have gone undetected. To address this issue, we survey here recent massive parallel sequencing (RNA-Seq) datasets from purified mouse and human islet cells. Our analysis reveals that the most strongly disallowed genes are similar in beta and alpha cells, with 11β-hydroxysteroid dehydrogenase (Hsd11b1) mRNA being essentially undetectable in both cell types. The analysis also reveals that several genes involved in cellular proliferation, including Yap1 and Igfbp4, and previously assumed to be disallowed in both beta and alpha cells, are selectively repressed only in the beta cell. The latter finding supports the view that beta cell growth is selectively restricted in adults, providing a mechanism to avoid excessive insulin production and the risk of hypoglycaemia. Approaches which increase the expression or activity of selected disallowed genes in the beta cell may provide the basis for novel regenerative therapies in type 2 diabetes.
Pancreatic beta cells detect small fluctuations in circulating glucose levels by coupling oxidative metabolism to the regulation of ATP-sensitive K+ channels, and consequently to Ca2+ influx (Rutter et al., 2015). The loss or dysfunction of beta cells contributes to all forms of diabetes mellitus, a disease affecting more than 400 m individuals worldwide1.
Studies in the early 1990s (Sekine et al., 1994) revealed that two “housekeeping” genes, LDHA and SLC16A1 (Monocarboxylate transporter-1, MCT-1), which are abundant in essentially all mammalian cell types and permit vigorous glycolytic flux during anaerobosis, are expressed at vanishingly low levels in beta cells. Subsequent studies by ourselves (Pullen et al., 2010; Pullen and Rutter, 2013) and others (Thorrez et al., 2011; Lemaire et al., 2016) have provided a list of ∼60 genes which are selectively disallowed in these cells, of which there is general consensus on a list of ∼11 genes (Pullen and Rutter, 2013). Re-expression of Slc16a1 or Ldha (Zhao and Rutter, 1998; Ishihara et al., 1999; Ainscow et al., 2000; Pullen et al., 2012) as well as the acyl-CoA thioesterase, Acot7 (Martinez-Sanchez et al., 2016) in the beta cell leads to defects in insulin secretion, suggesting that the silencing of these genes in beta cells is likely to be functionally relevant.
Previous studies to identify islet disallowed genes have, however, analyzed whole islet transcriptome data (Pullen et al., 2010; Thorrez et al., 2011). Because islets are composed of multiple cell types (Elayat et al., 1995), this has not given a clear picture for any one cell type: the possibility consequently exists that certain genes may be less “disallowed” in the less abundant islet endocrine cells (notably alpha and delta) than in beta cells.
It has therefore been of interest to explore this question using datasets recently made available from highly purified islet cell types (Benner et al., 2014; Adriaenssens et al., 2016; DiGruccio et al., 2016), as well as our own, previously unpublished data. With this goal in mind, we have used a similar strategy to previous analyses but taking advantage of the increased dynamic range of RNA-Seq and the purified cell type datasets to reveal a more detailed insight of genes disallowed alpha and beta cells.
While we confirm that many previously identified islet disallowed genes are indeed disallowed in both alpha and beta cells, we also reveal a number of genes which are expressed at a far lower level in beta cells and whole islets. Strikingly, 11β-hydroxysteroid dehydrogenase (Hsd11b1) a critical enzyme for the conversion of the inactive precursor of corticosterone, 11-dehydrocorticosterone (11-DHC) (Seckl and Walker, 2001), is found to be remarkably weakly expressed in both beta and alpha cells, as well as in delta cells. This suggests that previous findings demonstrating actions of 11-DHC on insulin secretion (Davani et al., 2000) may reflect a requirement for non-beta cells in the activation of this molecule.
We also note that genes involved in cell proliferation (e.g., Yap1, Igfbp4) are more selectively inactivated in beta than alpha cells. This is consistent with an evolutionarily driven mechanism to suppress the proliferation of beta cells in adults and thus a risk of life-threatening hypoglycaemia.
Materials and Methods
Identification of Disallowed Genes in Adult Mouse Tissues and Purified Cells
RNA-Seq datasets for a range of normal mouse tissues were obtained from public repositories and our own data. These included: 10 datasets for FACS-purified alpha cells from three studies (Benner et al., 2014; Adriaenssens et al., 2016; DiGruccio et al., 2016); 14 FACS-purified beta cell datasets from four studies (Benner et al., 2014; Adriaenssens et al., 2016; DiGruccio et al., 2016; E-MTAB-2266); and 15 islet datasets from our own and other studies (Kone et al., 2014; GSE90531). Non-islet tissue datasets covering brain, heart, kidney, liver, lung, spleen, and thymus were included for comparison. Dispersed, FACS-purified cells contain many differences to dissected tissues including lack of blood and endothelial cells. To control for these differences, eight datasets of FACS-purified hepatocytes were included. Islet, alpha, and beta cell data were derived from data series with GEO or ArrayExpress accession numbers: GSE54973, E-MTAB-2266, GSE80673, GSE76017, GSE90531, and E-MTAB-2791. Non-islet tissue and cell data were derived from data series with accession numbers: GSE36025, GSE74747, GSE65207, and GSE68806.
Where data were in paired-end format, the second read was discarded for each fragment to allow a more direct comparison with the single-ended data. Reads were mapped to the mouse genome (GRCm38) using HiSat2 (Kim et al., 2015) and annotated transcripts quantified with featureCounts (Liao et al., 2014). Differential expression analysis was performed using DESeq2 (Love et al., 2014). Alpha cell, beta cell, and islet data were then subjected to pairwise comparisons against each non-islet tissue and cell type. There is an acknowledged problem that detecting differential expression by significance alone can result in the selection of genes with small but consistent changes which are of questionable biological significance while excluding larger but more variable changes. In an attempt to overcome this, the fold-change and significance were combined into a single metric, the π-value (Xiao et al., 2014). The π-values for each pairwise comparison were combined where the direction of change was consistent (e.g., genes whose expression was higher or lower in beta cells than all other non-islet datasets) by taking the π-value closest to zero (i.e., least combined significance and fold-change) for each gene. Genes were ranked by π-value allowing the most consistently under- and over-expressed genes to be identified at both ends of the resulting list. Since the number of datasets included for alpha cells, beta cells, and islets was not consistent, this is likely to have affected the number of differentially expressed genes detected in each case. For this reason, we have concentrated on comparing the ‘most disallowed’ genes at the top of the ranked gene lists in these cell types rather than using a significance/fold-change threshold to define disallowance.
The top 50 disallowed genes from alpha and beta cells were combined and hierarchical clustering performed on the basis of fold-change in gene expression in each tissue/cell type relative to islets. Weighted gene co-expression network analysis (WGCNA) was also performed on variance-stabilizing transformed count data exported from DESeq2 (Langfelder and Horvath, 2008). Functional characterization of gene lists using Gene Ontology (GO) Biological Process and PANTHER protein class was performed through the PANTHER website (Mi et al., 2017).
Developmental Regulation of Disallowed Genes
Transcriptomes of FACS-purified beta cells at different peri- and post-natal maturation stages were generated using the mIns1-H2b-mCherry reporter line (Benner et al., 2014) (Jax # 028589). Breeders homozygous for mIns1-H2b-mCherry were crossed to wild-type C57bl6 mates to ensure offspring were uniformly hemizygous for the mCherry reporter. Islets from single litters (a mix of male and female is expected) were pooled for each sample to obtain sufficient material. Pooled islets were dissociated, sorted and collect in Trizol for RNA isolation and library construction. Each time point was collected in the morning and FACS-sorted at the conclusion of the islet prep the same day. Each time point was done in duplicate or triplicate. Data are deposited under GEO accession number GSE88779.
Gene Expression Association with Type 2 Diabetes Status in Human Islets
To investigate whether the expression of the human orthologues of these genes are altered during the progression of type 2 diabetes, we exploited RNA-Seq data on pancreatic islets isolated from 89 human donors (Fadista et al., 2014). Gene expression was quantified using the pipeline described above using the GRCh37 genome, and differential expression relative to diabetes status determined by DESeq2. The top 20 beta cell disallowed genes were investigated for a significant association between expression and diabetes status.
Once the datasets were assembled, the analysis yielded the genes which were most enriched and most disallowed in alpha and beta cells (Figure 1A). The results were ranked by the π-value which combined both the log fold change and adjusted p-value into a single metric. It should be noted that our analysis is based on pairwise comparisons between the islet cell type and each of the non-islet tissues, and the smallest π-value from each of these comparisons is used for each gene. Since this analysis is based on the smallest π-value from comparisons with each other tissue it identifies genes which are consistently over- or under-represented in a particular cell type, rather than those with the largest median difference.
FIGURE 1. (A) Ranked list of the top 20 most disallowed (selectively down-regulated) genes in pancreatic alpha and beta cells, with overlap between the two lists highlighted by lines. (B) Functional annotation of the top 50 disallowed gene in each cell type on the basis of Gene Ontology (GO) Biological Process.
The similarity between the 20 genes most consistently repressed in alpha and beta cells, relative to non-islet tissues, is notable. Thus, 11 of these genes are common between the two lists. Functional classification of the top 50 disallowed genes in each cell type again revealed the similarity between the two lists, with ‘cellular process,’ ‘metabolic process’ and ‘biological regulation’ being the top three GO Biological Process terms in both lists (Figure 1B).
To study the similarities between the lists for the two cell types, the top 50 disallowed genes for each were merged and clustered on the basis of differential expression across all cell/tissue types (Figure 2A). While 27 of these genes were shared between the two cell types, the differences between the lists were equally interesting. One cluster contained five genes which were disallowed selectively in beta cells, but not alpha cells (Fam13a, Acot1, Ablim3, Chn2, Rxrg). The genes in this cluster exhibit lower expression in purified beta cells than in whole islets, and higher expression in alpha cells. A smaller cluster of three genes (Pcsk5, Yap1, Colex12) has a similar expression pattern, whereas the majority of the other beta cell disallowed genes also had low expression in alpha cells.
FIGURE 2. (A) Hierarchical cluster analysis of a combined list of the top 50 disallowed genes from both alpha and beta cells. Membership of the top 50 alpha and/or beta cell list is indicated by the gray bars on the right. Heatmap shading denotes fold-change in gene expression relative to islets, with blue representing up- and red down-regulation. A cluster of 50 genes selectively repressed in beta but not alpha cells was identified by weighted gene co-expression network analysis (WGCNA), and genes were functionally annotated on the basis of Biological Process (B) and Protein Class (C).
We sought to further investigate the differences in selectively repressed genes between alpha and beta cells using WGCNA of normalized count data from all the tissue/cell types. This approach identified a single module, containing 50 genes, which was significantly down-regulated in beta cells (p = 3 × 10-17) and islets (p = 2 × 10-12) but not alpha cells (p = 0.3; Supplementary Table S1). This module included most of the genes in the clusters described above. Searching for enrichment of GO terms revealed the enzyme-linked receptor signaling pathway (p = 0.023). This observation provides insights into possible differences in the proliferative capacity of alpha and beta cells.
Functional classification of the genes within this module showed that many were associated with metabolic processes (Figure 2B). A preponderance of nucleic acid binding, transcription factor and signaling molecules among the protein classes (Figure 2C) also indicates that selective silencing of this module in beta cells may contribute to the regulation of beta cell identity.
Figure 3 shows the intersection of data between previous analyses and the current analysis of islet disallowed genes (A) and between the different cell types and islets (B) and reveals that while there is considerable overlap between these datasets, we also noted genes not previously classed as disallowed.
FIGURE 3. Comparison of disallowed gene expression in isolated mouse islet cells versus intact islets. Venn diagram showing the overlap between the top 50 disallowed islet genes from this study (Yellow) with lists from previous studies by Pullen et al. (2010; Red) and Thorrez et al. (2011; Green) (A). The overlap between the top 50 disallowed genes from islets (Green), alpha (Yellow), and beta cells (Red) is also shown (B).
We next compared the levels of expression of five of the genes disallowed in alpha and/or beta cells (Figure 4). Of these, the most dramatically disallowed is Hsd11b1 with nearly a 1000-fold gradient existing between expression in brain versus purified alpha or beta cells, where mRNA levels were at or below the level of detection (<10 normalized counts). Relatively abundant expression in other islet cell types presumably explains its considerably higher expression in isolated islets.
FIGURE 4. Expression of selected disallowed genes across different mouse tissues. Gene expression is presented in normalized read counts per gene along with mean expression level across all tissues are shown for some of the genes highlighted in this study.
Interestingly, both Igfbp4 and Yap1 are approximately one order of magnitude more highly expressed in alpha than beta cells, indicating that the disallowance of these genes is likely to be more relevant to the beta cell phenotype than that of alpha cells.
Next, we sought to determine the time points at which members of the family of disallowed genes become inactivated in beta cells through analyzing RNA-Seq data from purified beta cells at various developmental stages (Figure 5). Strikingly, of those sampled, the majority were already very weakly expressed by embryonic day 18 (E18). Gas1 provided an exception in that mRNA levels encoding this enzyme declined (albeit from low levels) between E18 and P7.
FIGURE 5. Changes in disallowed gene expression during post-natal development. Expression of disallowed genes in purified insulin-positive cells from pre- and post-natal mice is shown along with the expression of glucokinase (Gck) for comparison. The six genes exhibiting the highest expression during this period are shown (A), and Gas1 is also shown as one of the few genes exhibiting a consistent down-regulation during this period (B).
Finally, we investigated whether increased expression of any of these disallowed genes is associated with type 2 diabetes (Table 1). Of the top 20 beta cell disallowed genes, increased expression of six was associated with diabetes with an adjusted p-value < 0.1, with IGF1 showing the most significant association.
TABLE 1. The association of expression of beta cell disallowed genes in human islets with type 2 diabetes status.
Comparison of Current Methodology to Previous Studies
The approach taken in this analysis differs from previous analyses in a number of ways, in addition to the use of data from purified cell types rather than isolated pancreatic islets. The largest difference is the use of RNA-Seq rather than microarray data. The principle advantages of RNA-Seq are the higher dynamic range (Zhao et al., 2014) and the fact that quantification is based on the whole length of the gene rather than the limited coverage of probesets. However, in addition to using a different platform, the analysis is also different in the current study.
Previous analyses have concentrated either on the statistical significance of the difference (Thorrez et al., 2011), or on a fold-change threshold (Pullen et al., 2010). The present analysis combined both measures into a π-value with the aim of capturing both highly significant (consistent) changes of a smaller magnitude and less significant (more variable) changes with a greater magnitude. Furthermore, ranking genes by this combined score is arguably more biologically meaningful than ranking by either significance or fold-change alone. We note, however, that demonstration of the relative strengths of this approach versus others will require future biological assessments, e.g., overexpression of candidate genes (Pullen et al., 2012; Martinez-Sanchez et al., 2016).
It is possible to assess the combined effect of these changes by comparing the current analysis of whole islets with previous analyses (Figure 3). The overlap between the top 50 islet disallowed genes from this study and the lists of 39 and 40 genes from the Pullen et al. (2010) and Thorrez et al. (2011) papers, respectively, is just two genes: Slc16a1 and Igfbp4. This is reassuring since there is strong evidence that Slc16a1 is specifically silenced in islets (Zhao et al., 2001), and that mis-expression of this gene in human (Otonkoski et al., 2007) and mouse (Pullen et al., 2012) beta cells causes inappropriate stimulation of insulin secretion in response to pyruvate.
The current list also contains six genes which overlap with the Pullen list but none with the Thorrez list. It is interesting to note that five genes common to both the Pullen and Thorrez lists are not included in the current islet list: Ldha, Cd302, Pdgfra, Cxcl12, and Oat. Of these, both Ldha and Oat are within the top 100 of the current list, so while their rank may be altered slightly in the current analysis, there is still strong evidence for their categorization as islet disallowed genes.
Pdgfra is an interesting disallowed gene because there is evidence that a gradual lowering in its expression plays a role in the age-dependent decline in beta cell proliferative capacity (Chen H. et al., 2011). However, Pdgfra also displays relatively low expression in liver and strikingly low expression in purified hepatocytes. While this does not detract from the important role of Pdgfra repression in beta cells, it does exclude it from our strict definition of islet disallowed genes.
Alpha and Beta Cell Disallowed Genes
The main aim of this study was to identify novel disallowed genes in alpha and beta cells, particularly those whose low expression may have been masked by expression in other islet cell types. When searching for genes specifically expressed in a particular cell type, low levels of contamination with other cell types have relatively little impact because they are unlikely to lead to false negatives and will only occasionally lead to false positives if the level of contamination is high. However, contamination is likely to be more of an issue when searching for cells specifically repressed in one cell type, since the expression in the contaminating cells will be easily detectable over the low background.
A striking feature of the lists of alpha and beta cell disallowed genes is their remarkable similarity. In contrast to enriched genes, which contain few similarities between cell types (Xin et al., 2016), the disallowed genes were notably similar (Figure 1). Looking at this aspect of cell identity suggests that alpha and beta cells are more similar than the impression given by looking at highly expressed genes, and is presumably an important aspect of reprogramming from an alpha cell to a beta fate (and vice versa) in some circumstances (Collombat et al., 2009; Thorel et al., 2010).
There is also considerable overlap between alpha cell, beta cell, and islet disallowed gene lists (Figure 3B), which includes several previously identified islet disallowed genes: Slc16a1, Fxyd1, Ndfg2, Gas1. However, it is particularly interesting to consider a few examples of disallowed genes which are either newly identified or ranked more highly in the current analysis of purified cells.
Hsd11b1 is the highest ranked beta cell disallowed gene, and second on the alpha cell list. While Hsd11b1 is also included on the islet disallowed list, the reason for its high rank in the purified cell types becomes apparent from a plot of its expression profile (Figure 4). Although the islet expression level is over an order of magnitude lower than any non-islet tissue, expression in alpha and beta cells is a further order of magnitude below this. Indeed, since the DESeq2 tool used in this analysis adds 0.5 to all normalized counts it can be seen from Figure 4 that no reads mapped to Hsd11b1 in two alpha cell, and three beta cell samples, and fewer than 10 normalized counts were detected in any alpha or beta cell sample. This level would normally be considered below the threshold of detectable expression, i.e., the gene is essentially silent in both cell types. Of note, although expression of Hsd11b1 is still very low in islets, the levels detected in the intact micro-organ are much higher than in purified alpha or beta cells. We find no evidence of expression in delta cells (DiGruccio et al., 2016) so this most likely reflects significant Hsd11b1 expression in other non-endocrine islet cell types (endothelial, pericytes, etc.).
HSD11B1 encodes an enzyme which converts inert cortisone into active cortisol (or 11-dehydroxycorticosterone to corticosterone in rodents), thus regulating local glucocorticoid action in metabolically active tissues (Seckl and Walker, 2001). Glucocorticoids have long been associated with adverse metabolic effects, including increased insulin resistance. Indeed, steroid-induced diabetes is a significant complication of glucocorticoid therapy (Hwang and Weiss, 2014). Hsd11b1 is reportedly the major regulator of tissue-specific glucocorticoid effects (Morgan et al., 2014). Overexpression of Hsd11b1 in either liver (Paterson et al., 2004) or adipose tissue (Masuzaki et al., 2003) renders mice insulin resistant, conversely knockout of this gene protects mice from glucocorticoid-induced insulin resistance (Morgan et al., 2014). In humans HSD11B1 hyperexpression has been associated with abnormal glucose metabolism and obesity in several studies (Nascimento et al., 2015). Whereas the diabetogenic effects of glucocorticoids have mainly been attributed to increased insulin resistance, there is some evidence that they involve beta cells directly.
Many reports of the effects of glucocorticoids on beta cells have highlighted inhibitory roles such as suppression of both insulin secretion (Chen F. et al., 2011) and cellular proliferation (Colvin et al., 2013) by the synthetic glucocorticoid, dexamethasone. Indeed, it would be adaptive for endocrine cells not to generate cortisol/corticosterone locally as this would preserve their ability to respond to circulating cortisol [derived from hypothalamus pituitary adrenal (HPA) axis activity] without any interference of local 11 HSD activity. After all, GCs are the functional antagonist of insulin and the two are cross-regulatory.
We have previously shown (Huising et al., 2011) that acute dexamethasone stimulation inhibits beta cell GLP1 receptor expression in a glucocorticoid receptor-dependent manner, confirming earlier findings. It is interesting to note that acute restraint stress elicited similar effects in mouse islets after 3 h, and these had normalized after 12 h. This suggests that acute stress transiently attenuates the sensitivity of beta cells to incretins by suppressing the expression of its receptors in order to facilitate the hyperglycemia induced by glucocorticoids. This acute, transient repression of insulin release by glucocorticoids is, in the above context, adaptive. However, in metabolic syndrome, the same mechanisms may rapidly become maladaptive.
However, other studies suggest that glucocorticoids may stimulate insulin secretion and/or beta cell survival in certain circumstances. For example, in vitro pre-treatment with low concentrations (200 nM) of corticosterone for 18 h increased glucose-stimulated insulin secretion from isolated islets (Hult et al., 2009). Dex also reportedly increased glucose-stimulated insulin secretion in islets isolated from rats that were treated in vivo with Dex (Rafacho et al., 2010), although it is difficult to dissociate effects on the islet from other systemic actions. Glucocorticoid pre-treatment also improved the function of transplanted human islets (Lund et al., 2008). Finally, overexpression of Hsd11b1 is reported to offer some protection against high-fat diet-induced glucose intolerance, through increasing both islet mass and function (Turban et al., 2012). Together, these findings suggest that the anti-inflammatory action of glucocorticoids may enhance beta cell function under certain conditions, although prolonged and/or high level treatment may interfere with both beta cell function and proliferation.
While we demonstrate here strikingly low expression of Hsd11b1 in both alpha and beta cells, there are apparently contradictory reports in the literature concerning the expression of this gene in beta cells and islets. Schmid et al. (2011) reported HSD11B1 expression at the protein level throughout rat islets and in INS-1 insulinoma cells, although this finding is sharply refuted in a commentary by Liu et al. (2011) who question the specificity of the antibody used. These authors also refer to their previous paper (Chowdhury et al., 2015) reporting HSD11B1 expression in glucagon-positive cells of wild-type mouse islets by immunofluorescence, and in islets by Western (immuno-) blotting. There are further reports of HSD11B1 detection in mouse islets after treatment with high-fat diet (Turban et al., 2012) or in ob/ob mice (Davani et al., 2000). It is difficult to reconcile these reports with the very low levels detected here in both purified alpha and beta cells by RNA-Seq. Firstly, there are clearly questions over the specificity of some of the antibodies used in the previous studies. Secondly, we do detect substantially higher expression in islets than in purified alpha or beta cells, and previous use of qPCR and Western blotting were mainly performed on isolated islets rather than purified cells. Our results suggest that the Hsd11b1 expression detected in whole islets comes almost exclusively from islet cell types other than alpha and beta cells. Of course, since RNA-Seq measures the transcript level it is possible that is does not accurately reflect the protein level. However, it is difficult to envisage how such a large discrepancy could occur. Finally, the cells used for the RNA-Seq analysis were all from mice on normal diet, and it is possible that Hsd11b1 expression is induced under certain conditions such as high-fat diet.
Interestingly, HSD11B1 levels were undetectable in 4/7 purified human beta cell preparations sampled by Blodgett et al. (2015) and at <1 RPKM in the remaining three samples. HSD11B1 mRNA was also absent in 2/7 purified human alpha cell samples, and at around 1 RPKM in the other five. The latter data indicate that the enzyme is also likely to be disallowed in both human beta and alpha cells, as well as in rodent cells. In each species this presumably provides a mechanism for “protecting” these cell types from the effects of 11-DHC, or ensuring that the actions of the precursor on the secretory function of these cells (Davani et al., 2000) is via the presence of this enzyme in neighboring islet cell types, as discussed above.
Mgll encodes monoglyceride lipase which is responsible for the unregulated hydrolysis of monoacylglycerides (MAG) in adipocytes and other tissues (Karlsson et al., 1997). Two other lipases, Adipose Triglyceride Lipase (ATGL) and Hormone Sensitive Lipase (HSL/LIPE) are both expressed in beta cells and have been proposed to stimulate insulin secretion through the production of lipid-derived coupling factors including MAG (Fex et al., 2009; Attané et al., 2016). By degrading this MAG, it is likely that expression of MGLL would remove this coupling factor and decrease insulin secretion. Interestingly, both Mgll and Lipe are dysregulated in beta cells null for PGC-1α/β. Knockout of these PPARγ co-activators impairs fatty-acid potentiated insulin secretion and was accompanied by down-regulation of Lipe and up-regulation of Mgll. It is tempting to speculate that maintained expression of Atgl and Lipe, along with suppression of Mgll are necessary to allow MAG to act as a coupling factor. It is interesting to note that overexpression of Acot7, another disallowed gene involved in lipid metabolism, impairs glucose-stimulated insulin secretion, although this appeared to be through increased ATP consumption rather than decreased MAG levels (Martinez-Sanchez et al., 2016).
Expression of PGC-1α is reduced in the islets from subjects with type 2 diabetes (Ling et al., 2008), and this is paralleled by a significant increase in Mgll expression in islets from diabetic subjects (Table 1). Although it is unclear whether there is a direct link between PGC-1α and Mgll expression, the dysregulation of lipid metabolism associated with these changes could contribute to the impaired beta cell function in type 2 diabetes.
Numerous studies have reported the pro-proliferative role of IGF-1 on beta cells (reviewed in Stewart et al., 2015), and exogenous IGF-1 treatment has been shown to increase beta cell proliferation (Sieradzki et al., 1988). However, the effects of locally produced IGF-1 may differ from those stimulated by raised circulating levels of the growth factor. Beta cell-specific overexpression of IGF-1 increased the survival and proliferation of beta cells after challenge with multiple low dose streptozotocin-treatment (George et al., 2002). Since this overexpression had little effect on systemic serum levels of IGF-1 it presumably acted in an autocrine/paracrine manner.
The action of IGF-1 is modulated by a family of IGF binding proteins which regulate its bioavailability. Igfbp4 is strongly disallowed in beta cells, and Igfbp2 is also within the top 50 disallowed genes for both alpha and beta cells. IGFBP4 is generally considered to have a negative regulatory role on IGF-1, for example, both IGFBP4 and IGFBP2 were found to be down-regulated through promoter hypermethylation in lung carcinomas (Sato et al., 2006). In addition to binding to IGF-1, thus preventing it binding to IGF receptors, IGFBP4 has also been reported to act independently of IGF (reviewed in Durai et al., 2006). However, there is growing evidence that IGFBP4 may differentially regulate proliferation in different contexts. Transfecting primary renal cancer cell lines with IGFBP4 increased proliferation, invasion, and motility, whereas knocking down the gene in metastatic renal cancer cell lines decreased proliferation (Ueno et al., 2011).
The complexity of its varying roles in other tissues makes it difficult to predict what the effect of Igfbp4 overexpression in islets would be. It also makes it difficult to suggest an adaptive reason for its relative down-regulation in beta cells, and how this might interact with the disallowance of Igf1. It is notable that increased expression of both IGF1 and IGFBP4 in islets show some association with type 2 diabetes (Table 1). Of the top 20 beta cell disallowed genes, IGF1 showed the strongest association with diabetes status (p = 0.0185), whereas the association with increased IGFBP4 was less strong (p = 0.0786). Future studies, in which Igfbp1 is selectively over-expressed in beta cells, will be needed to solve this conundrum.
Yap1 is a potent driver of cell growth via the Hippo pathway (Hansen et al., 2015). Interestingly, a survey of human adult alpha and beta cell expression data reveals levels of ∼1 RPKM for YAP1 in both cell types, albeit with slightly lower levels in beta cells (Blodgett et al., 2015). Nonetheless, in mouse cells, we found that Yap1 and Pdgfra were selectively repressed in beta cells, consistent with a lower proliferative capacity of these cells relative to alpha cells.
We therefore suspect that its suppression in beta cells provides a mechanism to inhibit beta cell expansion, and hence fatal hypoglycaemia. The absence of Yap1, Igf1 and potentially Igfbp4 may also contribute to the low proliferative index of beta cells in rodents but, more particularly in man, when cell division is barely detectable after the age of ∼21 (Cnop et al., 2010). Moreover, this, alongside the absence of Pdgfra (see above) may equally restrict compensatory beta cell growth in times of metabolic need including insulin resistance provoked by diet, aging or pregnancy. Thus, mechanisms designed to increase the expression or these receptors, or to engage their downstream signaling pathways, may provide new approaches toward regeneration of beta cell mass in type 2 diabetes.
Developmental Regulation of Expression
It is interesting to investigate at what stage in the developing islet expression of these genes is down-regulated relative to other tissues. It has previously been reported in developing rat islets that Ldha and Slc16a1 expression peaked at post-natal day 1 before reducing to near adult levels by P21 (Thorrez et al., 2011). The current investigation of purified insulin-positive cells from mice at stages E18 – P21 did not reveal consistent down-regulation of disallowed genes during this period (Figure 5A). Most of the beta cell disallowed genes exhibited low expression (<5 RPKM) by E18 and showed little clear change during the post-natal period. One exception was Gas1 which did show considerable down-regulation between E18 and P7. However, this pattern did not extend across the disallowed genes investigated. One possible explanation for these differing results is that Thorrez et al. (2011) measured expression in whole islets versus the purified insulin-positive cells reported here. It is therefore possible that the down-regulation observed in whole islets represents increases in the proportion of beta (and other endocrine) cells in the developing islet, whilst the expression of disallowed genes within these remains unchanged during this period, as demonstrated experimentally here.
The present report provides several novel insights. Notably, we identify a novel member of the beta cell disallowed group, Hsd11b1. Secondly, we reveal that assumptions as to the expression of established members of the family, Yap1 and Igf4bp1, as being disallowed in the alpha as well as the beta cell, were misplaced. We further identify a module of 50 genes, including Yap1, which are selectively down-regulated in beta but not alpha cells. Our findings may provide the basis of new approaches toward improving beta cell identity and function in type 2 diabetes.
TP and GR conceived and planned the study. TP and MH selected and analyzed data. TP, GR, and MH wrote and edited the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supported by a Wellcome Trust Senior Investigator Award (WT098424AIA), MRC Programmes (MR/J0003042/1; MR/L020149/1, “DIVA”), Biotechnology and Biological Sciences Research Council (BB/J015873/1), Diabetes UK (BDA/11/0004210, BDA/15/0005275) and Imperial Confidence in Concept (ICiC) grants, and a Royal Society Wolfson Research Merit Award to GR and a Career Development Award from the Juvenile Diabetes Research Foundation to MH (2-2013-54) TP was a DRWF Post-doctoral Fellow.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fgene.2017.00041/full#supplementary-material
TABLE S1 | Selectively down-regulated beta cell gene module.
Adriaenssens, A. E., Svendsen, B., Lam, B. Y. H., Yeo, G. S. H., Holst, J. J., Reimann, F., et al. (2016). Transcriptomic profiling of pancreatic alpha, beta and delta cell populations identifies delta cells as a principal target for ghrelin in mouse islets. Diabetologia 59, 2156–2165. doi: 10.1007/s00125-016-4033-1
Ainscow, E. K., Zhao, C., and Rutter, G. A. (2000). Acute overexpression of lactate dehydrogenase-A perturbs beta-cell mitochondrial metabolism and insulin secretion. Diabetes Metab. Res. Rev. 49, 1149–1155. doi: 10.2337/diabetes.49.7.1149
Attané, C., Peyot, M.-L., Lussier, R., Poursharifi, P., Zhao, S., Zhang, D., et al. (2016). A beta cell ATGL-lipolysis/adipose tissue axis controls energy homeostasis and body weight via insulin secretion in mice. Diabetologia 59, 2654–2663. doi: 10.1007/s00125-016-4105-2
Benner, C., van der Meulen, T., Cacéres, E., Tigyi, K., Donaldson, C. J., and Huising, M. O. (2014). The transcriptional landscape of mouse beta cells compared to human beta cells reveals notable species differences in long non-coding RNA and protein-coding gene expression. BMC Genomics 15:620. doi: 10.1186/1471-2164-15-620
Blodgett, D. M., Nowosielska, A., Afik, S., Pechhold, S., Cura, A. J., Kennedy, N. J., et al. (2015). Novel observations from next generation RNA sequencing of highly purified human adult and fetal islet cell subsets. Diabetes Metab. Res. Rev. 64, 3172–3181. doi: 10.2337/db15-0039
Chen, F., Zhu, Y., Tang, X., Sun, Y., Jia, W., Sun, Y., et al. (2011). Dynamic regulation of PDX-1 and FoxO1 Expression by FoxA2 in dexamethasone-induced pancreatic β-cells dysfunction. Endocrinology 152, 1779–1788. doi: 10.1210/en.2010-1048
Chowdhury, S., Grimm, L., Gong, Y. J. K. Y., Wang, B., Li, B., Srikant, C. C. B., et al. (2015). Decreased 11β-hydroxysteroid dehydrogenase 1 level and activity in murine pancreatic islets caused by insulin-like growth factor I overexpression. PLoS ONE 10:e0136656. doi: 10.1371/journal.pone.0136656
Cnop, M., Hughes, S. J., Igoillo-Esteve, M., Hoppa, M. B., Sayyed, F., van de Laar, L., et al. (2010). The long lifespan and low turnover of human islet beta cells estimated by mathematical modelling of lipofuscin accumulation. Diabetologia 53, 321–330. doi: 10.1007/s00125-009-1562-x
Collombat, P., Xu, X., Ravassard, P., Sosa-Pineda, B., Dussaud, S., Billestrup, N., et al. (2009). The ectopic expression of Pax4 in the mouse pancreas converts progenitor cells into α and subsequently β cells. Cell 138, 449–462. doi: 10.1016/j.cell.2009.05.035
Colvin, E. S., Ma, H.-Y., Chen, Y.-C., Hernandez, A. M., and Fueger, P. T. (2013). Glucocorticoid-induced suppression of β-cell proliferation is mediated by Mig6. Endocrinology 154, 1039–1046. doi: 10.1210/en.2012-1923
Davani, B., Khan, A., Hult, M., Martensson, E., Okret, S., Efendic, S., et al. (2000). Type 1 11 -hydroxysteroid dehydrogenase mediates glucocorticoid activation and insulin release in pancreatic islets. J. Biol. Chem. 275, 34841–34844. doi: 10.1074/jbc.C000600200
DiGruccio, M. R., Mawla, A. M., Donaldson, C. J., Noguchi, G. M., Vaughan, J., Cowing-Zitron, C., et al. (2016). Comprehensive alpha, beta and delta cell transcriptomes reveal that ghrelin selectively activates delta cells and promotes somatostatin release from pancreatic islets. Mol. Metab. 5, 449–458. doi: 10.1016/j.molmet.2016.04.007
Durai, R., Davies, M., Yang, W., Yang, S., Seifalian, A., Goldspink, G., et al. (2006). Biology of insulin-like growth factor binding protein-4 and its role in cancer (review). Int. J. Oncol. 28, 1317–1325. doi: 10.3892/ijo.28.6.1317
Fadista, J., Vikman, P., Laakso, E. O., Mollet, I. G., Esguerra, J., Lou Taneera, J., et al. (2014). Global genomic and transcriptomic analysis of human pancreatic islets reveals novel genes influencing glucose metabolism. Proc. Natl. Acad. Sci. U.S.A. 111, 13924–13929. doi: 10.1073/pnas.1402665111
Fex, M., Haemmerle, G., Wierup, N., Dekker-Nitert, M., Rehn, M., Ristow, M., et al. (2009). A beta cell-specific knockout of hormone-sensitive lipase in mice results in hyperglycaemia and disruption of exocytosis. Diabetologia 52, 271–280. doi: 10.1007/s00125-008-1191-9
George, M., Ayuso, E., Casellas, A., Costa, C., Devedjian, J. C., Bosch, F., et al. (2002). Beta cell expression of IGF-I leads to recovery from type 1 diabetes. J. Clin. Invest. 109, 1153–1163. doi: 10.1172/JCI12969
Huising, M. O., Pilbrow, A. P., Matsumoto, M., van der Meulen, T., Park, H., Vaughan, J. M., et al. (2011). Glucocorticoids differentially regulate the expression of CRFR1 and CRFR2α in MIN6 insulinoma cells and rodent islets. Endocrinology 152, 138–150. doi: 10.1210/en.2010-0791
Hult, M., Ortsäter, H., Schuster, G., Graedler, F., Beckers, J., Ploner, A., et al. (2009). Short-term glucocorticoid treatment increases insulin secretion in islets derived from lean mice through multiple pathways and mechanisms. Mol. Cell. Endocrinol. 301, 109–116. doi: 10.1016/j.mce.2008.09.038
Ishihara, H., Wang, H., Drewes, L. R., Wollheim, C. B., DeFronzo, R., Virkarmaki, A., et al. (1999). Overexpression of monocarboxylate transporter and lactate dehydrogenase alters insulin secretory responses to pyruvate and lactate in beta cells. J. Clin. Invest. 104, 1621–1629. doi: 10.1172/JCI7515
Karlsson, M., Contreras, J. A., Hellman, U., Tornqvist, H., and Holm, C. (1997). cDNA cloning, tissue distribution, and identification of the catalytic triad of monoglyceride lipase: evolutionary relationship to esterases, lysophospholipases, and haloperoxidases. J. Biol. Chem. 272, 27218–27223. doi: 10.1074/jbc.272.43.27218
Kone, M., Pullen, T. J., Sun, G., Ibberson, M., Martinez-Sanchez, A., Sayers, S., et al. (2014). LKB1 and AMPK differentially regulate pancreatic β-cell identity. FASEB J. 28, 4972–4985. doi: 10.1096/fj.14-257667
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi: 10.1093/bioinformatics/btt656
Ling, C., Del Guerra, S., Lupi, R., Rönn, T., Granhall, C., Luthman, H., et al. (2008). Epigenetic regulation of PPARGC1A in human type 2 diabetic islets and effect on insulin secretion. Diabetologia 51, 615–622. doi: 10.1007/s00125-007-0916-5
Liu, J.-L., Srikant, C. B., and Chowdhury, S. (2011). Is 11β-HSD1 expressed in islet β-cells and regulated by corticotropin-releasing hormone? Proc. Natl. Acad. Sci. U.S.A. 108, E1390. doi: 10.1073/pnas.1116146109
Lund, T., Fosby, B., Korsgren, O., Scholz, H., and Foss, A. (2008). Glucocorticoids reduce pro-inflammatory cytokines and tissue factor in vitro and improve function of transplanted human islets in vivo. Transpl. Int. 21, 669–678. doi: 10.1111/j.1432-2277.2008.00664.x
Martinez-Sanchez, A., Pullen, T. J., Chabosseau, P., Zhang, Q., Haythorne, E., Cane, M. C., et al. (2016). Disallowance of Acot7 in β-cells is required for normal glucose tolerance and insulin secretion. Diabetes Metab. Res. Rev. 65, 1268–1282. doi: 10.2337/db15-1240
Masuzaki, H., Yamamoto, H., Kenyon, C. J., Elmquist, J. K., Morton, N. M., Paterson, J. M., et al. (2003). Transgenic amplification of glucocorticoid action in adipose tissue causes high blood pressure in mice. J. Clin. Invest. 112, 83–90. doi: 10.1172/JCI17845
Mi, H., Huang, X., Muruganujan, A., Tang, H., Mills, C., Kang, D., et al. (2017). PANTHER version 11: expanded annotation data from gene ontology and reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 45, D183–D189. doi: 10.1093/nar/gkw1138
Morgan, S. A., McCabe, E. L., Gathercole, L. L., Hassan-Smith, Z. K., Larner, D. P., Bujalska, I. J., et al. (2014). 11β-HSD1 is the major regulator of the tissue-specific effects of circulating glucocorticoid excess. Proc. Natl. Acad. Sci. U.S.A. 111, E2482–E2491. doi: 10.1073/pnas.1323681111
Nascimento, F. V., do Piccoli, V., Beer, M. A., Frankenberg, A. D., von, Crispim, D., et al. (2015). Association of HSD11B1 polymorphic variants and adipose tissue gene expression with metabolic syndrome, obesity and type 2 diabetes mellitus: a systematic review. Diabetol. Metab. Syndr. 71, 2711–2716. doi: 10.1186/S13098-015-0036-1
Otonkoski, T., Jiao, H., Kaminen-Ahola, N., Tapia-Paez, I., Ullah, M. S., Parton, L. E., et al. (2007). Physical exercise-induced hypoglycemia caused by failed silencing of monocarboxylate transporter 1 in pancreatic beta cells. Am. J. Hum. Genet. 81, 467–474. doi: 10.1086/520960
Paterson, J. M., Morton, N. M., Fievet, C., Kenyon, C. J., Holmes, M. C., Staels, B., et al. (2004). Metabolic syndrome without obesity: hepatic overexpression of 11beta-hydroxysteroid dehydrogenase type 1 in transgenic mice. Proc. Natl. Acad. Sci. U.S.A. 101, 7088–7093. doi: 10.1073/pnas.0305524101
Pullen, T. J., Khan, A. M., Barton, G., Butcher, S. A., Sun, G., and Rutter, G. A. (2010). Identification of genes selectively disallowed in the pancreatic islet. Islets 2, 89–95. doi: 10.4161/isl.2.2.11025
Pullen, T. J., Sylow, L., Sun, G., Halestrap, A. P., Richter, E. A., and Rutter, G. A. (2012). Overexpression of monocarboxylate transporter-1 (Slc16a1) in mouse pancreatic β-cells leads to relative hyperinsulinism during exercise. Diabetes Metab. Res. Rev. 61, 1719–1725. doi: 10.2337/db11-1531
Rafacho, A., Marroquí, L., Taboga, S. R., Abrantes, J. L. F., Silveira, L. R., Boschero, A. C., et al. (2010). Glucocorticoids in Vivo induce both insulin hypersecretion and enhanced glucose sensitivity of stimulus-secretion coupling in isolated rat islets. Endocrinology 151, 85–95. doi: 10.1210/en.2009-0704
Rutter, G. A., Pullen, T. J., Hodson, D. J., and Martinez-Sanchez, A. (2015). Pancreatic β-cell identity, glucose sensing and the control of insulin secretion. Biochem. J. 466, 203–218. doi: 10.1042/BJ20141384
Sato, H., Yazawa, T., Suzuki, T., Shimoyamada, H., Okudela, K., Ikeda, M., et al. (2006). Growth regulation via insulin-like growth factor binding protein-4 and -2 in association with mutant K-ras in lung epithelia. Am. J. Pathol. 169, 1550–1566. doi: 10.2353/ajpath.2006.051068
Schmid, J., Ludwig, B., Schally, A. V., Steffen, A., Ziegler, C. G., Block, N. L., et al. (2011). Modulation of pancreatic islets-stress axis by hypothalamic releasing hormones and 11beta-hydroxysteroid dehydrogenase. Proc. Natl. Acad. Sci. U.S.A. 108, 13722–13727. doi: 10.1073/pnas.1110965108
Seckl, J. R., and Walker, B. R. (2001). Minireview: 11β-hydroxysteroid dehydrogenase type 1— a tissue-specific amplifier of glucocorticoid action 1. Endocrinology 142, 1371–1376. doi: 10.1210/endo.142.4.8114
Sekine, N., Cirulli, V., Regazzi, R., Brown, L., Gine, E., Tamarit-Rodriguez, J., et al. (1994). Low lactate dehydrogenase and high mitochondrial glycerol phosphate dehydrogenase in pancreatic beta-cells. Potential role in nutrient sensing. J. Biol. Chem. 269, 4895–4902.
Sieradzki, J., Fleck, H., Chatterjee, A. K., and Schatz, H. (1988). Stimulatory effect of insulin-like growth factor-I on [3H]thymidine incorporation, DNA content and insulin biosynthesis and secretion of isolated pancreatic rat islets. J. Endocrinol. 117, 59–62. doi: 10.1677/JOE.0.1170059
Stewart, A. F., Hussain, M. A., García-Ocaña, A., Vasavada, R. C., Bhushan, A., Bernal-Mizrachi, E., et al. (2015). Human β-cell proliferation and intracellular signaling: part 3. Diabetes 64, 1872–1885. doi: 10.2337/db14-1843
Thorel, F., Népote, V., Avril, I., Kohno, K., Desgraz, R., Chera, S., et al. (2010). Conversion of adult pancreatic α-cells to β-cells after extreme β-cell loss. Nature 464, 1149–1154. doi: 10.1038/nature08894
Thorrez, L., Laudadio, I., Van Deun, K., Quintens, R., Hendrickx, N., Granvik, M., et al. (2011). Tissue-specific disallowance of housekeeping genes: the other face of cell differentiation. Genome Res. 21, 95–105. doi: 10.1101/gr.109173.110
Turban, S., Liu, X., Ramage, L., Webster, S. P., Walker, B. R., Dunbar, D. R., et al. (2012). Optimal elevation of β-cell 11β-hydroxysteroid dehydrogenase type 1 is a compensatory mechanism that prevents high-fat diet–induced β-cell failure. Diabetes Metab. Res. Rev. 61, 642–652.
Ueno, K., Hirata, H., Majid, S., Tabatabai, Z. L., Hinoda, Y., and Dahiya, R. (2011). IGFBP-4 activates the Wnt/beta-catenin signaling pathway and induces M-CAM expression in human renal cell carcinoma. Int. J. Cancer 129, 2360–2369. doi: 10.1002/ijc.25899
Xiao, Y., Hsiao, T.-H., Suresh, U., Chen, H.-I. H., Wu, X., Wolf, S. E., et al. (2014). A novel significance score for gene selection and ranking. Bioinformatics 30, 801–807. doi: 10.1093/bioinformatics/btr671
Xin, Y., Kim, J., Ni, M., Wei, Y., Okamoto, H., Lee, J., et al. (2016). Use of the fluidigm C1 platform for RNA sequencing of single mouse pancreatic islet cells. Proc. Natl. Acad. Sci. U.S.A. 113, 3293–3298. doi: 10.1073/pnas.1602306113
Zhao, C., and Rutter, G. A. (1998). Overexpression of lactate dehydrogenase A attenuates glucose-induced insulin secretion in stable MIN-6 β-cell lines. FEBS Lett. 430, 213–216. doi: 10.1016/S0014-5793(98)00600-0
Zhao, C., Wilson, M. C., Schuit, F., Halestrap, A. P., and Rutter, G. A. (2001). Expression and distribution of lactate/monocarboxylate transporter isoforms in pancreatic islets and the exocrine pancreas. Diabetes Metab. Res. Rev. 50, 361–366. doi: 10.2337/diabetes.50.2.361
Keywords: Hsd11b1, glucocorticoids, islets of langerhans, beta cells, alpha cells, diabetes mellitus, disallowed genes
Citation: Pullen TJ, Huising MO and Rutter GA (2017) Analysis of Purified Pancreatic Islet Beta and Alpha Cell Transcriptomes Reveals 11β-Hydroxysteroid Dehydrogenase (Hsd11b1) as a Novel Disallowed Gene. Front. Genet. 8:41. doi: 10.3389/fgene.2017.00041
Received: 16 December 2016; Accepted: 23 March 2017;
Published: 10 April 2017.
Edited by:Jeff Schwartz, Griffith University, Australia
Reviewed by:Zong Wei, Salk Institute for Biological Studies, USA
Stefan Bonn, German Center for Neurodegenerative Diseases (HZ), Germany
Copyright © 2017 Pullen, Huising and Rutter. 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) or licensor 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.