Quantitative-Proteomic Comparison of Alpha and Beta Cells to Uncover Novel Targets for Lineage Reprogramming

Type-1 diabetes (T1D) is an autoimmune disease in which insulin-secreting pancreatic beta cells are destroyed by the immune system. An emerging strategy to regenerate beta-cell mass is through transdifferentiation of pancreatic alpha cells to beta cells. We previously reported two small molecules, BRD7389 and GW8510, that induce insulin expression in a mouse alpha cell line and provide a glimpse into potential intermediate cell states in beta-cell reprogramming from alpha cells. These small-molecule studies suggested that inhibition of kinases in particular may induce the expression of several beta-cell markers in alpha cells. To identify potential lineage reprogramming protein targets, we compared the transcriptome, proteome, and phosphoproteome of alpha cells, beta cells, and compound-treated alpha cells. Our phosphoproteomic analysis indicated that two kinases, BRSK1 and CAMKK2, exhibit decreased phosphorylation in beta cells compared to alpha cells, and in compound-treated alpha cells compared to DMSO-treated alpha cells. Knock-down of these kinases in alpha cells resulted in expression of key beta-cell markers. These results provide evidence that perturbation of the kinome may be important for lineage reprogramming of alpha cells to beta cells.


Introduction
Type-1 diabetes (T1D) is a chronic autoimmune disease affecting 35 million patients worldwide. In T1D, insulin-secreting pancreatic beta cells are destroyed by autoreactive immune cells [1,2]._ENREF_1 The most common treatment for T1D is daily injection of insulin; however, this treatment cannot always ensure optimal glucose homeostasis, leading to complications such as blindness, heart disease, limb amputation, and ultimately death [3]. Another therapeutic strategy involves transplantation of pancreatic islets [4] via infusion through the portal vein into the liver, but high cost [5], limited donor availability, and beta-cell toxicity [6,7] of immunosuppressive drugs severely restrict the use of this treatment protocol.
Since T1D is characterized by extreme loss of beta-cell mass, replenishing the beta-cell population by converting other pancreatic cell types, such as alpha cells, to beta cells may be a viable therapeutic strategy [8,9]. For example, lineage reprogramming of pancreatic alpha cells to beta cells by ectopic expression of transcription factor, Pax4, has been shown to restore normoglycemia in mouse models of diabetes [10]. An alternative, and potentially more therapeutically promising, approach for lineage reprogramming involves the use of small molecules. This approach, which we term chemical transdifferentiation [11], is potentially safer, as it does not require viral delivery, and offers greater dosage and temporal control. However, identification of proteins or pathways that can be targeted for chemical transdifferentiation requires knowledge of the similarities and differences in the molecular and physiological architecture of the initial state and the final, transdifferentiated state.
To gain such knowledge, we compared the transcriptomes and proteomes of alpha and beta cell lines. In addition to analyzing the basal state in cell culture, we also used two small molecules, BRD7389 and GW8510, that we previously discovered to induce insulin expression in alpha cells ( Figure S1 in File S1) [12,13]. BRD7389 and GW8510 inhibit several kinases, such as the RSK family and CDK2, suggesting that lineage reprogramming of alpha cells to insulin-producing cells may involve downregulation of specific phosphorylated proteins in alpha cells. We reasoned that a signature of such phosphorylated proteins would be its presence at a higher level in alpha cells compared to beta cells; further, BRD7389 and GW8510 treatment of alpha cells should lower the levels of such phosphorylated proteins, as the alpha cells progress towards an intermediary cell state. To identify phosphoproteins with these signatures, we compared the phosphoproteomes of alpha and beta cell lines, and also examined the alterations in the phosphoproteomes of alpha cell line upon treatment with either BRD7389 or GW8510.
We identified two proteins, BR serine/threonine kinase 1 (Brsk1) and calcium/calmodulin-dependent protein kinase kinase 2 (Camkk2), which were each phosphorylated at higher levels for specific phosphosites in alpha cells as compared to those in beta cells. Moreover, both BRD7389 and GW8510 treatment of alpha cells led to decreased phosphorylation at these phosphosites for both proteins. Knock-down of either kinase was sufficient to induce mRNA and protein expression of beta-cell markers, including insulin, in alpha cells. Furthermore, inhibitors of Camkk2 pathway induced beta-cell markers in alpha cells. These studies highlight the importance of fully characterizing cell states using both genetic and proteomic technologies in developing a cellular reprogramming strategy.

Cell culture and reagents
Mouse pancreatic cell lines aTC1 and bTC3 were obtained from ATCC and cultured at 37uC and 5% CO 2 in low-glucose (1 g/L) DMEM supplemented with 10% fetal bovine serum (Hyclone), 50 U/mL penicillin and 50 mg/mL streptomycin. Media was changed every 3-4 days. KN62 and KN93 were obtained from Sigma-Aldrich.

Microarray analyses
GenePattern was used to perform Gene Set Enrichment Analysis (GSEA) on previously reported data [14], using the gene set database c2.all.v3.0.symbols.gmt, chip platform mouse 430a_2.chip, and 1000 permutations. GENE-E [15] was used for comparative marker selection, and to generate the heat map. The data are available as GEO accession number GSE36379.
Cell mitochondrial respiration and glycolytic activity aTC1 cells were treated with 0.1% DMSO or 0.85 mM BRD7389 for 3 days. Cells were seeded at 40,000/well into polyornithine pre-coated 96-well cell culture plates from Seahorse Bioscience. Seahorse Bioscience XF-96 extracellular flux analyzer was used to measure OCR and ECAR over time and analyte addition. Cell glycolytic activity was assessed by adding 10 mM glucose at 34 min, while the maximum glycolytic activity was induced with 5 mM oligomycin at 55 min. Glycolytic reserve was determined by adding the glucose analogue 2-deoxyglucose (100 mM) into the cell media at 76 min. Mitochondrial respiration was assessed by adding: 1 mM oligomycin into the media at 34 min, 1 mM of the ionophore CCCP, and a combination of 5 mM rotenone and 5 mM antimycin A at 76 min. aTC1 cells media contained 0.1% DMSO.
Non-phosphorylated and phosphorylated peptides were prepared as previously described [16]. Briefly, peptides were fractionated offline by strong cation exchange (SCX) chromatography on an Akta Purifier 10 system (GE Healthcare) using a polysulfoethyl A strong cation exchange (SCX) column from PolyLC (25069.4 mm, 5 mm particle size, 200 A pore size). A 160 min gradient was utilized for the fractionation with using 7 mM KH 2 PO 4 , pH 2.65, 30% MeCN as solvent A and 7 mM KH 2 PO 4 , pH 2.65, 350 mM KCl, 30% MeCN as Solvent B. For proteome analyses 5% of each SCX fraction was aliquoted, and fractions were combined into 24 proteome samples that were subsequently desalted using StageTips. [17] For phosphoproteome analysis, 95% of each SCX fraction was used, and the SCX fractions were combined into 12 phosphopeptide samples. All SCX fractions were subsequently desalted with reversed phase tC18 SepPak columns.
Phosphopeptide enrichment was completed using immobilized metal affinity chromatography (IMAC) as previously described [16]. Briefly, peptides were reconstituted in 200 ml of 40% MeCN, 0.1% formic acid and incubated for with 10 ml of packed Phos-Select beads (Sigma) for 1 hr. Subsequently, IMAC beads were loaded on C18 StageTips, washed twice with 40% MeCN, 0.1% formic acid and once with 50 ml 1% formic acid. Phosphorylated peptides were then transeluted to the C18 portion of the StageTip by washing with 70 ml of 500 mM K 2 HPO 4 (pH 7.0) three times. StageTips were then washed with 50 ml of 1% formic acid and peptides were eluted with 80 ml of 50% MeCN/0.1% formic acid and dried to completeness.
Peptide samples were reconstituted in 3% MeCN/1% formic acid analyzed by LC-MS/MS using an Agilent 1200 LC coupled online to an LTQ Orbitrap Velos mass spectrometer (Thermo Fisher Scientific). Peptides were loaded onto a fused-silica capillary column (New Objective) packed in-house with 14 cm of C18 reversed phase media (3 mm ReproSil-Pur C18-AQ media, Dr. Maisch GmbH) and eluted into the mass spectrometer using a 70 min linear gradient (,0.29%B/min) from 10% solvent A (0.1% formic acid) to 30% solvent B (0.1% formic acid/90% acetonitrile). The mass spectrometer was operated in the data dependent mode where an MS1 scan (R = 60 K) was acquired with the Orbitrap analyzer followed by acquisition of CID MS/ MS scans with the iontrap analyzer on the top 16 most abundant ions. An MS1 ion target of 1610 6 ions and an MS2 target of 1610 4 ions were utilized for acquisition. The maximum ion time was set to 1 s for MS1 scans and set to 100 ms for MS/MS scans. The dynamic exclusion time was 120 s, the repeat count was set to 2, and the repeat duration was set to 20 s. Monoisotopic precursor selection and non-peptide monoisotopic recognition functions were enabled.
All proteome and phosphoproteome data were processed together using the MaxQuant software package v1.0.13.13 and searched against the International Protein Index protein sequence database (IPI version 3.70, mouse). For searching, the enzyme specificity was set to trypsin, the number of missed cleavages was set to 2, and the MS/MS tolerance was set to 0.5 Da. Oxidation of methionines, N-terminal protein acetylation, and phosphorylation of serines, threonines, and tyrosines were searched as variable modifications, while carbamidomethylation of cysteines was searched as a fixed modification. Peptides harboring oxidized methionines or N-terminal acetylations were used for protein quantification. Log 2 SILAC ratios of bTC3 versus aTC1 cells followed a normal distribution that was fitted using least-squares regression. Mean and standard deviation values derived from the gaussian fit were used to calculate p-values.
The peptide and protein false discovery rate (FDR) was set to 0.01 and the site FDR was set to 1. Phosphosite localization was completed as described in Mertins et al. [18].

siRNA experiments
Mouse Brsk1 and Camkk2 Gene Solution siRNAs (Qiagen) were used to perform the gene knock-down experiments. Quantitect primer sets for each gene were used to determine gene expression. Mouse aTC1 cells were plated in 96-well Corning (3340) plates at 40,000 cells/well in 100 mL DMEM. 0.3 mL/well Lipofectami-neTM RNAiMAX (Invitrogen) and Opti-MEM Media was used to transfect the mix of 4 different siRNA constructs into cells. Cells were incubated in transfection mix for 6 h at 37uC, before changing to fresh media. Cells were cultured for 3 days, followed by qPCR or immunohistochemistry. Statistical significance was determined using t-test.

shRNA experiments
The different shRNA lentiviruses were obtained from the RNAi Consortium (TRC; http://www.broadinstitute.org/rnai/public/) in 96-well format, with approximately 108 viral particles/mL per well. Hairpin identities are listed in Table S1 in File S1. Mouse aTC1 cells were plated in 96-well Corning (3340) plates at 40,000 cells/well in 100 mL standard culture media. The next day, polybrene was added to each well (6 mg/mL), and cells were spininfected with 1 mL virus at 2,250 rpm for 30 min at 30uC. Media was changed 4 h later to fresh media. The next day, media containing 1 mg/mL puromycin was added. Puromycin-supplemented media was changed every 3 days. After 10 days in culture, cells were lysed and mRNA extracted using Qiagen RNeasy 96 Kit. Statistical significance was determined using t-test.

Immunocytochemistry
Cells were fixed with 4% PFA for 15 min, followed by a 0.2% Triton-X-100 permeabilization for 20 min, and blocking with PBS supplemented with 2% BSA for 2 h. Fixed cells were then incubated with a mix of 1:200 rabbit anti-Pdx1 (Abcam) and 1:500 Endocrine-Cell Phosphoproteomics PLOS ONE | www.plosone.org guinea pig anti-insulin (Sigma) overnight at 4uC. As secondary antibodies, AlexaFluor594 anti-rabbit and AlexaFluor488 (Invitrogen) anti-guinea pig were used. Images were acquired on an ImageXpress Micro automated microscope (Molecular Devices). Exposure settings: 600 ms for Pdx1 (TxRed), 200 ms for insulin (FITC), and 8 ms for DAPI. Statistical significance was determined using t-test.

Gene expression
Total mRNA from siRNA and compound-treated cells were extracted with Qiagen RNeasy Plus Mini Kit. qPCRs were performed with Power SybrGreen PCR Master Mix and an Applied Biosystems 7900HT plate reader. Mouse primers were obtained from IDT and are listed in Table S2 in File S1.

Experimental strategy and rationale
We were interested in comparing the genetic and proteomic levels of various transcription factors in alpha and beta cells, as such proteins are linked with lineage reprogramming. Thus, we compared the gene-expression profiles of the mouse alpha-cell (aTC1) and beta-cell (bTC3) lines ( Figure 1A). In parallel, we compared the proteomes and phosphoproteomes of alpha cells and beta cells using stable isotope labeling by amino acids in cell culture (SILAC; Figure 1B) [19][20][21]. We used a three-state labeling strategy, with two states corresponding to the individual cell populations, and the third state corresponding to a mixed proteome from alpha and beta cells. This approach enabled four protein intensity measurements for alpha cells and beta cells using only two replicates ( Figure S2 and S3 in File S1). We analyzed the phosphoproteomes of alpha cells and beta cells by enrichment of phosphopeptides using immobilized metal-affinity chromatography ( Figure 1B), and also determined the changes in the phosphoproteome of alpha cells upon treatment with BRD7389 (1 h, 120 h) or GW8510 (120 h) ( Figure 1C).

Transcriptome analysis and phenotypic studies
We used previously described gene-expression profiling data [14], where the transcript levels of ,14,000 genes of alpha and beta cells were determined using GeneChip Mouse Genome 430A arrays (Affymetrix). A similar comparison of alpha and beta cell transcriptomes has been reported, using cell sorting of primary human islets [22]. To identify the genes most differentially expressed between alpha cells and beta cells, we used comparative marker selection [23], which employs statistical approaches such as false-discovery rate (FDR) and family-wise error rate (FWER) to account for multiple-hypothesis testing. From these top genes, we identified the transcription factors that were differentially expressed in alpha cells and beta cells ( Figure S4 in File S1). A number of transcription factors known to be involved in alpha-and beta-cell specification were identified, including Pdx1, Mnx1, and Irx1. To examine the pathways involved with differentially expressed genes, we used gene-set enrichment analysis (GSEA) [24,25], which determines whether a pre-defined set of genes shows statistically significant, concordant differences between two biological states (Figure 2A, B). We found that, among other gene sets, genes corresponding to oxidative phosphorylation, hypoxia, branched-chain amino acid catabolism, and electron transport were enriched in alpha cells (Figure 2A, Table S3 in File S1). In beta cells, gene sets corresponding to membrane trafficking and antigen presentation and processing were enriched ( Figure 2B, Table S4 in File S1).
Since GSEA suggested that the alpha cell line may have a higher metabolic rate than beta cells, we compared the mitochondrial respiration and glycolytic activity in aTC1 cells, BRD7389-treated aTC1 cells, and two beta cell lines: bTC3 and INS-1E cells. We observed higher glycolysis rate, maximum glycolytic capacity, and glycolytic reserve in aTC1 cells, as measured by extracellular acidification rate ( Figure 2C, S5 in File S1). We also determined mitochondrial respiration in cell populations by quantifying extracellular oxygen consumption rates (OCR) in the presence of various mitochondrial inhibitors [26]. aTC1 cells had elevated basal respiration rates, energy coupling efficiency, and maximum respiration rate compared to the beta cell lines ( Figure 2D). Further, treatment of aTC1 cells with BRD7389 led to a decrease in the overall glycolytic activity and mitochondrial respiration rates, a phenotype reminiscent of the beta cell lines ( Figure 2D) [27][28][29][30], further supporting the hypothesis that BRD7389 directs alpha cells towards an intermediate cell state.

Proteome and PTM coverage and analysis
We performed quantitative proteomics on the alpha cells, compound-treated alpha cells, and beta cells using SILAC [19][20][21]. Briefly, cells were cultured in media with amino acids containing light, medium, or heavy isotopes of arginine and lysine. Cell lysates from these three SILAC states were mixed and proteins were reduced, alkylated, and digested with trypsin. The SILAC media (e.g., light, medium, or heavy) from which each observed tryptic peptide originates determines the isotopic label state of the lysine or arginine residue in each peptide. The peptides were subsequently fractionated using strong cation exchange chromatography and phosphopeptides were enriched using immobilized metal-ion chromatography, where phosphopeptides chelate to a bead-immobilized iron (III) complex. Protein ratios were determined from the Protein Groups table generated by the MaxQuant [31,32] software package v1.0.13.13. Two independent SILAC triple-labeling experiments were completed for both the proteome and phosphoproteome analyses of alpha cells and beta cells, with Pearson correlation coefficient r = 0.9 between replicates ( Figure S3 in File S1). In total, 3,250 proteins and 9,978 phosphopeptides were quantified across both biological replicates. Approximately 10% of proteins had greater than a three-fold difference in abundance between alpha and beta cells (Table 1), including known hormones (e.g., insulin, glucagon) and kinases ( Table 2). Exploration of role of these proteins may provide further insight into the similarities and differences between alpha and beta cells. Neuropeptide Y and somatostatin, which are known to inhibit insulin secretion, were expressed at higher levels in beta cells compared to that in alpha cells [33][34][35]. Pdx1 and amylin, which play a key role in several beta-cell functions, were also expressed at higher levels in beta cells. A comparison of the phosphoproteomes of alpha and beta cells revealed interesting differences. For example, while glucagon is expressed at a higher level in alpha cells, the phosphorylation levels at serines 150 and 152 are higher in beta cells. Further, while insulin is expressed at higher levels in beta cells, the level of phosphorylation at threonine 51 is higher in alpha cells. We were able to detect known p38 MAP kinase-mediated phosphorylation of transcription factors MafA and Pdx1 at serine342 and serine 269, respectively [36][37][38]. The phosphoproteome of alpha cells was not significantly perturbed in the first hour of BRD7389 treatment, as only 2% of phosphopeptides showed alteration in levels. However, after five days of BRD7389 treatment, 9% of phosphopeptides were upregulated and 5% were downregulated. BRD7389 and GW8510 inhibit RSK and CDK2 kinases, respectively. Hence, we expected a decrease in phosphorylation of protein substrates of Rsk and Cdk2 upon BRD7389 or GW8510 treatment of alpha cells. Accordingly, we observed a decrease in S6 and EF2 phosphorylation, known substrates of RSK and CDK2, respectively (Table S5 in File S1).

Suppression of Brsk1 and Camkk2 elicits beta cell features in aTC1 cells
Although several kinases were differentially expressed in alpha and beta cells, we focused our attention on BRSK1 and CAMKK2, as they were more highly phosphorylated in alpha cells compared to beta cells, and in untreated alpha cells compared to compound-treated cells (Table S6 in File S1). We were unable to quantify absolute protein levels of these kinases in our proteomic analyses, but gene expression analysis suggested higher mRNA levels of Camkk2 in beta cells. Since the role of these proteins in maintenance of alpha cell state was unknown, we initiated our investigation by performing siRNA-mediated knock-down of Brsk1 or Camkk2 in aTC1 cells, and determined the expression of alpha and beta cell-specific genes. Brsk1 and Camkk2 mRNA levels were reduced by 64% and 85%, respectively, after three-day siRNA transfection ( Figure S6 in File S1). Expression of Camkk2 was also reduced in siBrsk1-treated cells (53%), suggesting that Brsk1 may be an upstream activator of Camkk2. Brsk1 knock-down resulted in the upregulation of the mRNA levels of genes important to beta-cell identity, including Pdx1, Ngn3, Nkx2.2, and glucokinase (Gck) (Figure 3A), while Camkk2 down-regulation elicited the expression of Nkx2.2, Pax4, MafB, and the hormones glucagon (Gcg) and insulin (Ins2) ( Figure 3B). Modest levels of Pdx1 protein expression were detected by immunofluorescence in siBrsk1-and siCamkk2-treated cells ( Figure S6 in File S1). Consistent with these observations, STO-609, a Camkk antagonist, and KN62 and KN93, two selective Ca 2+ /calmodulin-dependent protein kinase II (CaMKII) inhibitors, also increased Pdx1 protein expression in aTC1 cells ( Figure 3C).
Since PDX1 protein levels after 3 days of gene knock-down were detectable but faint, we performed a 10-day knock-down of Brsk1, Camkk2, and Stk11, an upstream modulator of Brsk1, using lentiviral shRNAs (Figure 4). When using hairpins targeting Brsk1 and Camkk2, a decrease in these kinases corresponded to a concomitant increase in Pdx1 mRNA expression ( Figure 4A-B). Hairpins targeting Stk11 (two out of three) also increased Pdx1 expression ( Figure 4C). The fact that not all hairpins had the same effect might reflect off-target silencing effects of these hairpins, which are part of the limitations of using small-interfering RNAs Figure 3. Knock-down of Brsk1 and Camkk2 elicit the expression of key beta cell-specific genes in aTC1 cells. mRNA expression of beta and alpha cell-specific genes was measured using qPCR after knock-down of (A) Brsk1 or (B) Camkk2. Fold changes were calculated compared to untransfected control aTC1 cells, and normalized to Gapdh expression. Significance was determined by t-test. *p,0.05, ***p,0.005. (C) aTC1 cells transfected with the indicated siRNA, or treated with the indicated compounds, were analyzed for Pdx1 protein expression by immunofluorescence. The percentage of Pdx1 + cells was calculated for each treatment. The significance was determined by ttest. *p,0.05, **p,0.01, and ***p,0.001. doi:10.1371/journal.pone.0095194.g003 [39]. Moreover, individual hairpins to Camkk2, Brsk1, and Stk11 elicited the expression of Pdx1 protein, as detected by immunofluorescence ( Figure 4D-S). PDX1-positive cells also co-stained for insulin when Brsk1 or Camkk2 were knocked down, and more modestly with knock-down of Stk11 (Figure 4 I,M,Q). Taken together, our results suggest a role for CAMK family proteins in maintaining an alpha-cell phenotype.

Discussion
In this study, we compared the transcriptomes, proteomes, and phosphoproteomes of a pancreatic alpha and beta cell line, with the aim of identifying new targets whose perturbation may induce expression of beta-cell markers in alpha cells. Gene-expression analysis suggested enhanced mitochondrial respiration and glycolysis in alpha cells compared to beta cells, which we confirmed experimentally. Furthermore, BRD7389 treatment of aTC1 lowered mitochondrial respiration and glycolytic activity as these cells progressed towards a state intermediate between alpha and beta cells. Beta cell lines, in particular INS-1E, have been reported to have high proton leak through the mitochondrial membrane [27], and to have strongly reduced lactate production [28][29][30], which is consistent with our observed low spare respiration and glycolytic capacity. The observed differences in the mitochondrial respiration and glycolytic activity could also arise from differences in the growth rate of the alpha and beta cell lines. Indeed, the beta cell lines used herein have a lower growth rate (doubling time: 2-3 days) compared to that of aTC1 cells (doubling time: 1-2 days).
Most lineage-reprogramming efforts involve overexpression of transcription factors that occur at higher levels in the target cell type. We have identified several transcription factors that are enriched in alpha cells and beta cells in our microarray data. Expected genes were readily identified; for example, Pdx1 plays an important role in beta cell maturation and function [40]. Interestingly, one of the identified transcription factors, motor neuron and pancreas homeobox 1 (Mnx1), was recently reported to play an important role in cell fate determination in zebrafish [41]. Physiological differences between alpha and beta cells are also reflected in the annotated gene sets enriched in these cells. Genes involved in respiration and branched-chain amino acid were enriched in alpha cells, while the genes associated with membrane trafficking and lysosomes were enriched in beta cells, pointing to the relative importance of these biological processes in the two cell types.
Our prior studies suggested that the inhibition of specific kinases in alpha cells can lead to expression of key markers of beta cells [12,13]. Here, we identified two kinases whose knock-down induced beta-cell markers in alpha cells. BRSK1, also known as SAD-B, is a serine/threonine kinase strongly expressed in mammalian forebrain to determine neuronal polarization [42]. BRSK1 is activated both in vitro and in vivo by serine/threonineprotein kinase 11 (STK11), also known as liver kinase B1 (LKB1), which plays an important role in the polarization of cortical axons [43][44][45]. STK11 is in turn a tumor suppressor and regulator of energy homeostasis through the activation of AMP-activated protein kinase (AMPK) [46]. A beta cell-specific deletion of Stk11 causes an increase in beta cell-mass and insulin secretion, effects mediated by an activation of the mTOR pathway [47]. In combination with STK11, CAMKK2 (Camkkb) phosphorylates and activates AMPK [48,49], while in combination with CAMKK1 (Camkka), it phosphorylates and activates CAMK1 and CAMK4 [50]. Camkk2 null mice are protected from obesity, insulin resistance, and glucose intolerance [51]. It is noteworthy that downregulation of Camkk2 in alpha cells increases the expression of both insulin and glucagon. Directed differentiation of embryonic stem cells into mature hormone secreting endocrine cells similarly involve polyhormonal cellular states, often described as a hallmark of non-fully differentiated cells. Interestingly, Camkk2 has been reported to play a role in neuronal differentiation [52] and bone remodeling by stimulating osteoblast formation [53] further supporting the role of this kinase to impact cell fate decision.
This chemical-genetic and phosphoproteomic study provides evidence that perturbing the alpha-cell kinome may be an attractive approach for lineage reprogramming, and points to a potential role for BRSK1, the CAMK kinase family, and STK11 in determining endocrine-cell identity. Gene silencing and compound treatments induced mRNA and protein expression of key beta cell-specific genes, such as Pdx1 and insulin; the LKB1/ BRSK1 and CAMKK2 pathways might converge at the level of AMPK, as suggested by previous work [46,51]. The presence of detectable levels of both glucagon and insulin in both aTC1 and bTC3 cells suggests that these cell lines are not perfect models of alpha and beta cells. However, we have previously shown that the protein targets identified from these cell lines using small molecules, BRD7389 and GW8510, translate in native pancreatic cells of humans. [12,13] These findings support the notion that these cell lines can lead to the discovery of novel protein targets relevant for the primary cells. A potential drawback of our study is the assumption that phosphorylation of kinases is correlated with kinase activation. Many examples exist where kinases are negatively regulated by phosphorylation (e.g., glycogen synthase kinase 3b). Further studies are required to dissect the opposing actions of kinases and reprogramming. This may include delineation of the signaling pathways of these kinases in endocrine cells, and understanding the signaling role of their phosphorylated forms. We anticipate that such studies may spur development of small molecules for lineage reprogramming of alpha cells to insulin-producing cells.