Single‐cell transcriptomes reveal characteristic features of human pancreatic islet cell types
2015; Springer Nature; Volume: 17; Issue: 2 Linguagem: Inglês
10.15252/embr.201540946
ISSN1469-3178
AutoresJin Li, Johanna Klughammer, Matthias Farlik, Thomas Penz, Andreas Spittler, Charlotte Barbieux, Ekaterine Berishvili, Christoph Bock, Stefan Kubicek,
Tópico(s)Single-cell and spatial transcriptomics
ResumoScientific Report21 December 2015Open Access Single-cell transcriptomes reveal characteristic features of human pancreatic islet cell types Jin Li Jin Li CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Johanna Klughammer Johanna Klughammer CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Matthias Farlik Matthias Farlik CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Thomas Penz Thomas Penz CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Andreas Spittler Andreas Spittler Medical University of Vienna, Anna Spiegel Forschungsgebäude, Vienna, Austria Search for more papers by this author Charlotte Barbieux Charlotte Barbieux Department of Surgery, Cell Isolation and Transplantation Center, Geneva University Hospitals, University of Geneva, Geneva, Switzerland Search for more papers by this author Ekaterine Berishvili Ekaterine Berishvili Department of Surgery, Cell Isolation and Transplantation Center, Geneva University Hospitals, University of Geneva, Geneva, Switzerland Search for more papers by this author Christoph Bock Corresponding Author Christoph Bock CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Department of Laboratory Medicine, Medical University of Vienna, Vienna, Austria Max Planck Institute for Informatics, Saarbrücken, Germany Search for more papers by this author Stefan Kubicek Corresponding Author Stefan Kubicek CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Christian Doppler Laboratory for Chemical Epigenetics and Antiinfectives, CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Jin Li Jin Li CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Johanna Klughammer Johanna Klughammer CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Matthias Farlik Matthias Farlik CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Thomas Penz Thomas Penz CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Andreas Spittler Andreas Spittler Medical University of Vienna, Anna Spiegel Forschungsgebäude, Vienna, Austria Search for more papers by this author Charlotte Barbieux Charlotte Barbieux Department of Surgery, Cell Isolation and Transplantation Center, Geneva University Hospitals, University of Geneva, Geneva, Switzerland Search for more papers by this author Ekaterine Berishvili Ekaterine Berishvili Department of Surgery, Cell Isolation and Transplantation Center, Geneva University Hospitals, University of Geneva, Geneva, Switzerland Search for more papers by this author Christoph Bock Corresponding Author Christoph Bock CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Department of Laboratory Medicine, Medical University of Vienna, Vienna, Austria Max Planck Institute for Informatics, Saarbrücken, Germany Search for more papers by this author Stefan Kubicek Corresponding Author Stefan Kubicek CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Christian Doppler Laboratory for Chemical Epigenetics and Antiinfectives, CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria Search for more papers by this author Author Information Jin Li1,‡, Johanna Klughammer1,‡, Matthias Farlik1,‡, Thomas Penz1,‡, Andreas Spittler2, Charlotte Barbieux3, Ekaterine Berishvili3, Christoph Bock 1,4,5 and Stefan Kubicek 1,6 1CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria 2Medical University of Vienna, Anna Spiegel Forschungsgebäude, Vienna, Austria 3Department of Surgery, Cell Isolation and Transplantation Center, Geneva University Hospitals, University of Geneva, Geneva, Switzerland 4Department of Laboratory Medicine, Medical University of Vienna, Vienna, Austria 5Max Planck Institute for Informatics, Saarbrücken, Germany 6Christian Doppler Laboratory for Chemical Epigenetics and Antiinfectives, CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, Vienna, Austria ‡These authors contributed equally to this work *Corresponding author. Tel: +43 1 40160 70070; Fax: +43 1 40160 970 000; E-mail: [email protected] *Corresponding author. Tel: +43 1 40160 70036; Fax: +43 1 40160 970 000; E-mail: [email protected] EMBO Reports (2016)17:178-187https://doi.org/10.15252/embr.201540946 PDFDownload PDF of article text and main figures. Peer ReviewDownload a summary of the editorial decision process including editorial decision letters, reviewer comments and author responses to feedback. ToolsAdd to favoritesDownload CitationsTrack CitationsPermissions ShareFacebookTwitterLinked InMendeleyWechatReddit Figures & Info Abstract Pancreatic islets of Langerhans contain several specialized endocrine cell types, which are commonly identified by the expression of single marker genes. However, the established marker genes cannot capture the complete spectrum of cellular heterogeneity in human pancreatic islets, and existing bulk transcriptome datasets provide averages across several cell populations. To dissect the cellular composition of the human pancreatic islet and to establish transcriptomes for all major cell types, we performed single-cell RNA sequencing on 70 cells sorted from human primary tissue. We used this dataset to validate previously described marker genes at the single-cell level and to identify specifically expressed transcription factors for all islet cell subtypes. All data are available for browsing and download, thus establishing a useful resource of single-cell expression profiles for endocrine cells in human pancreatic islets. Synopsis We provide a resource of single-cell transcriptomes from primary human pancreatic islets and describe the detection of specific marker genes, cellular heterogeneity, and mouse–human species differences. We performed single-cell RNA-seq on 70 cells sorted from primary human pancreatic islets. Bioinformatic analysis of single-cell transcriptomes rediscovered all the major pancreatic cell types including new cell type-specific marker genes. All the data are provided as an online resource for browsing and download. Introduction Located within the pancreas, the islets of Langerhans are composed of endocrine cells expressing glucagon (alpha cells), insulin (beta cells), somatostatin (delta cells), pancreatic polypeptide (PP cells), and ghrelin (epsilon cells). Furthermore, they are heavily vascularized and innervated, and in contact with the surrounding acinar and ductal cells of the exocrine pancreas. Pancreatic islets function as highly specialized micro-organs that monitor and maintain blood glucose homeostasis. While damage to beta cells causes diabetes, the other pancreatic cell types may also contribute to pathogenesis in ways that are not well understood. Recent studies showed that both alpha 1 and delta cells 2 have the potential to replenish beta cell mass in animal models. Development of diabetes correlates with global changes in the transcriptome of pancreatic islets 3. These gene expression changes could reflect alterations in the cell subtype composition of the islet and/or changes in the transcriptomes of beta cells or other individual cell types. Analyzing islet cell-specific gene expression changes has the potential to shed light on the etiology of diabetes. Recently, alpha and beta cell purification protocols from human 456 and mouse islets 78 have yielded initial maps of cell type-specific transcriptomes. The available transcriptome datasets further comprise primary mouse and human alpha cells, beta cells, and delta cells, a number of rodent alpha and beta cell lines, and one human beta cell line 49101112. Despite the rapid progress in this field, a comprehensive transcriptome database for individual human islet cell types is still missing, and no transcriptome data are currently available for PP cells. Recent advances in next-generation sequencing and library preparation enabled for the first time the transcriptome characterization of single cells from primary tissue. For example, this approach was successfully used to establish transcriptome profiles and dissect cell type heterogeneity for primary tissue obtained from the lung 13, the spleen, and the brain 1415. Here, we used single-cell RNA-seq to establish a comprehensive transcriptome database for the cell types that are present in primary human pancreatic islets. Principal component analysis in combination with visualization as biplots identified alpha cells, beta cells, delta cells, PP cells, acinar cells, and pancreatic duct cells directly from the single-cell transcriptome profiles. We illustrate the utility of this resource by discovering novel cell type-specific marker genes, and we identified human-specific expression patterns in alpha and beta cells. All data are readily available for user-friendly online browsing and download to foster research on pancreatic islet biology and diabetes-related mechanisms in human. Results and Discussion Single-cell transcriptomes recapitulate pancreatic endocrine cell types Primary human pancreatic islets of Langerhans were disassociated into single cells, and these cells were sorted into individual wells of a 96-well plate by FACS 16. The Smart-Seq2 protocol 17 was then applied to obtain single-cell transcriptomes. Following the generation and amplification of cDNA, we determined the levels of beta-actin expression by qRT–PCR and selected all cell-containing wells for library preparation and next-generation sequencing (Fig 1A). Seventy cells were sequenced in total, of which 64 cells passed quality control (see 3) and were included in the analysis (Fig EV1A and B, and Dataset EV1). We obtained an average of 12.7 million high-quality reads per single cell, of which 62.9% aligned to the human reference genome. RNA expression levels were calculated using the BitSeq software which uses RPKM normalization and corrects for non-uniform read distribution along the transcripts (e.g., 3-prime bias) 18. Data quality was validated by assessing the relation between expression level and transcript length in native RNA (Fig EV1C) as well as ERCC spike-in controls (Fig EV1D). While transcript length and expression level were not correlated in the ERCC spike-in controls, we detected a negative correlation (r = −0.405) in the native RNA which was in the range of what had been previously reported as biologically significant finding 19. However, a potential bias due to transcript length normalization cannot be completely excluded; therefore, comparing expression levels of different transcripts/genes should be performed with caution. To define global similarities among the single cells and the marker genes that drive these similarities, we performed principal component analysis (PCA) on the transcriptome dataset and displayed the results as biplots. PCA on the full dataset separates a group of 18 cells based on high glucagon (GCG) and transthyretin (TTR) expression and a group of 9 cells expressing pancreatic polypeptide (PPY) from a heterogeneous group of 37 cells (Fig 1B). In a second PCA on the 37 yet undefined cells, we identified a group of 12 cells with high insulin (INS) expression, a group of 11 cells characterized by PRSS2, CTRB2, REG3A, REG1A, and REG1B and a group of two somatostatin (SST)-expressing cells. In a third PCA on the remaining 12 undefined cells, a group of 8 cells was characterized by keratin18 (KRT18) and keratin8 (KRT8). Based on the expression profiles of the identified marker genes, we were able to uniquely assign 60 out of 64 single-cell transcriptomes to the alpha, beta, delta, PP, acinar, or ductal cell type (Fig 1C). Figure 1. Single-cell transcriptomes recapitulate the major pancreatic cell types Workflow for obtaining and analyzing single-cell RNA-seq data from human pancreatic islets. Iterative PCA/biplot-based approach for the identification of cell types and cell type-defining transcripts from single-cell RNA-seq data. Expression (scaled RPKM values) of cell type-defining genes as identified in (B) across all single cells. Transcripts and single cells are grouped by cell type as identified in (B). Display of transcriptional similarity between all single cells by MDS. The coloring scheme is based on the cell types as identified in (B). Relative expression (scaled RPKM value) of canonical marker genes for the 6 identified pancreatic islet cell populations represented by bubble size and projected onto the MDS profile. Download figure Download PowerPoint Click here to expand this figure. Figure EV1. Statistical analysis of single-cell RNA-seq data Number of detected transcripts and total aligned reads for each single cell. The red line denotes 500 transcripts, below which samples were excluded from the analysis. Scatter plots displaying the correlation between the number of input ERCC RNA molecules and measured RPKM values in four representative single cells. The Pearson correlation (r) is noted in the upper left corner. Scatter plots correlating raw read counts (left) and RPKM normalized expression values (right) with transcript length. All adequately covered transcripts (> 25 reads) of all 64 samples were included in this analysis. The observed negative correlation after RPKM normalization lies in the range of what had been reported previously 19. Raw counts and RPKM normalized expression values for 10 groups of ERCC spike-in controls. The amount of molecules spiked into the sequencing reaction is constant within one group, whereas the length of the transcripts varies considerably. Within each group, equal expression values across different transcript lengths thereby confirm that RPKM normalization does not systematically penalize longer transcripts in the assessed ERCC transcript length range of ˜200 to 2,000 bp. ERCC spike-in controls covered by more than 25 reads are indicated in blue and those with ≤ 25 reads in red. MDS displaying transcriptional similarity between corresponding cell types of a published dataset 20 (prefix: ext) and the current dataset. Download figure Download PowerPoint As an additional validation of our cell type classification, we visualized the global transcriptional similarity of individual pancreatic cells by multidimensional scaling (MDS), where each single-cell transcriptome was colored by the cell type derived from PCA (Fig 1D). When mapped upon the MDS plot, the known cell type-specific marker genes INS, GCG, PPY, SST, REG1A, and KRT8 show the expected expression patterns, with different amounts of variability within the subgroups (Fig 1E). The validity of our single-cell RNA-seq dataset was further confirmed in direct comparison to an external dataset consisting of bulk RNA-seq data for whole islet, beta, and acinar cells 20. Using MDS, we show high transcriptional similarity between the corresponding cell types of both datasets (Fig EV1E). The expression information of individual cells and merged expression values for each cell type is available in Dataset EV2. To rule out technical reasons as a major source of gene expression variability, we identified presumably pure alpha and beta cells among the assessed single cells (Fig EV2A). Their transcription profiles were used to simulate transcriptomes with defined percentages of alpha and beta cell contribution (Fig EV2B). Individual alpha and beta cells were then compared to these virtual transcriptomes to estimate upper limits for potential cross-contamination (Fig EV2C–E). All beta cell transcriptomes were found to be free from any alpha cell contribution, whereas beta cell profiles could explain a small proportion (< 3%) of the variance observed in 8 of the 18 alpha cells studied. However, given that these alpha cells further show higher unexplained variance, it is likely that they are characterized by high inherent variability rather than cross-contamination from beta cells. We conclude that the differences between alpha and beta cell heterogeneity are in line with biological rather than technical effects which supports the hypothesis that alpha cells might be more plastic than beta cells 4. Click here to expand this figure. Figure EV2. Assessing cross-contamination between alpha and beta cells Scatter plot displaying single alpha and beta cells, 500-cell islet samples, as well as bulk islet and beta cell samples from published datasets according to their weighted mean of scaled expression values in alpha and beta cell-specific profile genes. The three selected profile cells for each cell type are indicated by their sample ID. Pure and mixed expression profiles consisting of 233 alpha cell-specific genes and 252 beta cell-specific genes. Alpha and beta cell-specific profiles are calculated based on the expression values of the three selected profile cells only, while profile genes were selected based on all single cells classified as alpha or beta cells, which is why the expression gradients in the mix profiles do not always follow the same direction. Profile correlation curves for each individual sample. The maximum of each curve defines the maximum variance that can be explained (y-axis) by the corresponding mix profile (x-axis) providing a measure for the composition of the respective sample. Diagram explaining the transition from profile correlation curves to sample composition estimates. The profile composition that explains most variance is linearly scaled to the maximum variance explained. Sample composition estimates for each assessed sample. The differences between the 500-cell islet samples and bulk samples might be explained by technical effects that enrich for alpha cells during islet cultivation, disassociation, and FACS purification. Download figure Download PowerPoint The heterogeneity within the different cell types was further explored by separate PCAs for each cell type (Appendix Fig S1). Particularly for endocrine cells, heterogeneity was mainly driven by expression differences of marker genes as identified in the initial cell type classification by PCA, suggesting that these cell types are characterized by a spectrum of marker gene expression levels. While this analysis provides evidence for transcriptional heterogeneity, more cells are needed to thoroughly characterize subgroups within the different cell types. A transcriptome resource to reveal marker genes of human pancreatic cell types To maximize the utility of our dataset for the identification of cell type-specific expression patterns, we generated a resource of genome browser tracks of all individual cells as well as cumulative tracks for the cell type clusters identified by PCA (http://islet-transcriptome.computational-epigenetics.org/). One interesting use of this resource is the analysis of master regulatory transcription factors, which are fundamental for the development and the maintenance of different pancreatic cell types based on animal models and human genetics. The genome browser tracks illustrate the beta cell-specific expression of PDX1, a master regulatory transcription factor directly controlling insulin expression. In contrast, the alpha lineage-defining factor ARX is expressed in both alpha and some PP cells (Fig 2A and Appendix Fig S2). Other transcription factors that are important for pancreas development have different degrees of cell type-specific expression in mature human islets, including pan-endocrine (PAX6), beta cell-specific (PAX4), and duct/delta (HHEX) patterns (Fig 2B). While MAFA is transcribed in beta cells specifically, we observed robust MAFB expression in alpha, beta, and delta cells. Half of the beta cells studied expressed MAFA and MAFB concomitantly. In addition to these previously described factors, we also observed cell type-specific expression for transcriptional regulators, which have not yet been extensively characterized in the endocrine pancreas. For example, MORF4L1 shows a similar pan-endocrine pattern to the canonical islet cell marker NEUROD1 (Fig 2C). A subset of alpha cells express IRX2 (Fig 2D), some beta cells show high expression of polycomb ring finger oncogene (BMI1) (Fig 2E), and PP cells can be characterized by the transcription factors ETV1 and MEIS1 (Fig 2F). Figure 2. Expression of cell type-specific transcription factors at single-cell resolution A. Merged UCSC Genome Browser tracks for the PDX1 and ARX loci. The respective tracks for all single cells are presented in Appendix Fig S2. B. Relative expression (scaled RPKM value) of important transcription factors represented by bubble size and projected onto the MDS profile. C–F. Cell type-specific expression of pan-endocrine (C), alpha cell (D), beta cell (E), and PP cell (F) transcription factors (red bar: mean expression). The statistical significance of the differential gene expression is presented in Appendix Fig S6. Download figure Download PowerPoint We further performed pairwise correlation analysis on transcript level to identify genes, of which the expression profiles correlate highly (r > 0.9) with those of the endocrine marker genes INS, GCG, SST, and PPY (Fig EV3). While several highly correlated genes could be identified for INS and SST (e.g., zinc transporter SLC39A4 and Notch pathway component DLK1 for INS and transcription factors NKX6-3, ZNF430 for SST), the expression profiles of GCG and PPY did not show high correlation with any other genes. Click here to expand this figure. Figure EV3. Assessing pairwise correlation of endocrine marker genesCorrelation matrix displaying all genes (y-axis) that are highly correlated (r > 0.9) with at least one of the endocrine marker genes (x-axis). Different transcripts of the same gene are indicated by “Tx”. Download figure Download PowerPoint To extend our analysis beyond transcription factors and known marker genes, we performed pairwise comparisons of cell type-specific transcriptomes by gene set enrichment analysis (Dataset EV3). Interestingly, we observed strong enrichment of a gene set containing the REST-binding motif in all endocrine cell types compared to acinar and ductal cells (Fig 3A). Most genes that contain the REST motif in their promoters are expressed in alpha, beta, delta, and PP cells, whereas they are repressed in ductal and acinar cells (Fig 3B). The transcriptional repressor REST targets the REST-binding motif. In line with the target gene expression pattern, REST is specifically expressed in ductal and acinar cells (Fig 3C and Appendix Fig S3). Figure 3. Single-cell transcriptomes reveal unique features of human islets A. Heatmap displaying the P-values obtained by pairwise Gene Set Enrichment Analysis (GSEA) for the REST-binding motif. B. Relative expression (scaled RPKM value) of genes contained in the REST-binding motif gene set. C. Merged UCSC Genome Browser tracks for REST. The respective tracks for all single cells are presented in Appendix Fig S3. D. Expression of DLK1 and GC in human islet cell types (red bar: mean expression). The statistical significance of the differential gene expression is presented in Appendix Fig S6. E–H. Co-staining of DLK1 (E, F) or GC (G, H) with insulin and with glucagon in representative human (E, G) and mouse (F, H) islets. Download figure Download PowerPoint Finally, based on pairwise differential expression analysis between the pancreatic cell types, genes with highly specific expression patterns were identified (Fig EV4 and Appendix Fig S4, Datasets EV4 and EV5). We then used these data to assess islet cell type-specific expression in two areas of high relevance for diabetes research–diabetes risk genes and mouse–human species differences. Click here to expand this figure. Figure EV4. Specific expression of selected marker genesRelative expression (scaled RPKM value) of interesting genes across all single cells represented by bubble size and projected onto the MDS profile as displayed in Fig 1D. Download figure Download PowerPoint Genomewide association studies (GWAS) have identified genomic loci conferring increased risk for the development of diabetes. We examined whether any of the diabetes-related genes predicted by GWAS were specifically expressed in one of the pancreatic islet cell types and genes differentially expressed between the endocrine and exocrine cell types (Fig EV5A). For both type 1 and type 2 diabetes, we identified GWAS genes with beta cell- and endocrine-specific expression. Other genes show broader expression patterns, emphasizing the complexity of functional annotation of diabetes GWAS results. Furthermore, key MODY (Mature Onset of Diabetes in Young) 21 genes PDX1, PAX4, INS, HNF1A, GCK are predominantly specific to beta cells (Fig EV5B). Click here to expand this figure. Figure EV5. Assessing cell type specificity of genes identified in diabetes-related GWAS Cell type specificity for genes reported in diabetes-related GWAS. Each gene reported in a diabetes-related GWAS (search for “Diabetes” on GWAS Catalog) was assigned to the pancreatic cell type in which it was found to be most specifically expressed. Ranking was performed as described in Appendix Fig S7B and Dataset EV6. Heatmap showing mean expression values for the most cell type-specific diabetes-associated GWAS genes in the different here identified human islet cell types. Specifically, only genes with a specificity rank lower than 500 (dashed line in panel A) are listed, and genes with equal expression in multiple cell types are not shown. The numbers in the colored boxes indicate the number of studies in which the respective gene has been reported. Heatmaps are colored by mean ln(RPKM), the mean of the natural logarithm of the RPKM values across all cells of the respective cell type. Download figure Download PowerPoint To investigate species-specific differences of alpha and beta cell transcriptomes, we assessed the degree to which the previously identified differentially expressed mouse genes 79 are also differentially expressed in human islets and vice versa (Appendix Fig S5). We found that the human alpha cell-specific gene group-specific component (vitamin D binding protein) GC and the human beta cell-specific gene DLK1 (Fig 3D) displayed opposite expression patterns as to what had been reported in mouse islet cells. To confirm the cell type-specific expression of DLK1 and GC, we performed immunofluorescence staining on both human and mouse pancreatic tissue sections. In human islets, DLK1 was specifically expressed only in insulin-positive cells (Fig 3E), whereas this protein was observed in glucagon-positive cells in mouse tissue (Fig 3F). Similarly, GC expression showed alpha cell specificity in human tissues (Fig 3G), whereas it was co-expressed with insulin in mouse tissues (Fig 3H). These results suggest that two of the most differentially expressed cell type-specific marker genes for human alpha and beta cells have opposite expression patterns in mouse islets. Pancreatic islets comprise different cell types with characteristic transcriptomes, which confounds transcriptome studies that focus on whole pancreatic islets in physiological and pathological conditions. Lineage-labeled transgenic mice have made it possible to obtain transcriptomes for highly pure alpha and beta cell populations in mouse. For human islets, however, cell type-specific enrichment strategies depend on the availability of specific antibodies. Efforts have been made to measure the transcription of individual genes in single human islet cells by qRT–PCR 22, but our dataset is the first to provide genomewide transcriptional information of human islets at single-cell resolution. Using single-cell data, we also for the first time defined the transcriptomes of human delta cells and PP cells, thereby providing reference transcriptomes for all major endocrine cell types in human pancreatic islets. We illustrated the practical utility of our resource and dataset by three case studies. First, after confirming the cell type specificity of the major transcription factors involved in pancreatic endocrine lineage determination, we identified transcripts encoding transcription factors expressed in islet cells. These include the pan-endocrine marker MORF4L1, alpha cell-specific IRX2, beta cell-specific BMI1, and PP cell-specific MEIS1 and ETV1. These data can provide the basis for future functional studies in the roles of these transcription factors in the pancreas and in diabetes. In a second example, we analyzed cell type-specific enrichment of previously characterized gene sets. The specific expression of REST-motif-containing genes in the endocrine cell types led us to identify the specific expression of the transcriptional repressor REST in the exocrine pancreas. REST recruits a large complex of chromatin regulators, including many factors that allow pharmacological modulation like histone deacetylases and the histone demethylase LSD1. REST repression in non-endocrine cells activates the promoters of important beta cell transcription factors
Referência(s)