The Human Pancreatic Islet Transcriptome: Expression of Candidate Genes for Type 1 Diabetes and the Impact of Pro-Inflammatory Cytokines
Type 1 diabetes (T1D) is an autoimmune disease in which pancreatic beta cells are killed by infiltrating immune cells and by cytokines released by these cells. Signaling events occurring in the pancreatic beta cells are decisive for their survival or death in diabetes. We have used RNA sequencing (RNA–seq) to identify transcripts, including splice variants, expressed in human islets of Langerhans under control conditions or following exposure to the pro-inflammatory cytokines interleukin-1β (IL-1β) and interferon-γ (IFN-γ). Based on this unique dataset, we examined whether putative candidate genes for T1D, previously identified by GWAS, are expressed in human islets. A total of 29,776 transcripts were identified as expressed in human islets. Expression of around 20% of these transcripts was modified by pro-inflammatory cytokines, including apoptosis- and inflammation-related genes. Chemokines were among the transcripts most modified by cytokines, a finding confirmed at the protein level by ELISA. Interestingly, 35% of the genes expressed in human islets undergo alternative splicing as annotated in RefSeq, and cytokines caused substantial changes in spliced transcripts. Nova1, previously considered a brain-specific regulator of mRNA splicing, is expressed in islets and its knockdown modified splicing. 25/41 of the candidate genes for T1D are expressed in islets, and cytokines modified expression of several of these transcripts. The present study doubles the number of known genes expressed in human islets and shows that cytokines modify alternative splicing in human islet cells. Importantly, it indicates that more than half of the known T1D candidate genes are expressed in human islets. This, and the production of a large number of chemokines and cytokines by cytokine-exposed islets, reinforces the concept of a dialog between pancreatic islets and the immune system in T1D. This dialog is modulated by candidate genes for the disease at both the immune system and beta cell level.
Published in the journal:
. PLoS Genet 8(3): e32767. doi:10.1371/journal.pgen.1002552
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1002552
Summary
Type 1 diabetes (T1D) is an autoimmune disease in which pancreatic beta cells are killed by infiltrating immune cells and by cytokines released by these cells. Signaling events occurring in the pancreatic beta cells are decisive for their survival or death in diabetes. We have used RNA sequencing (RNA–seq) to identify transcripts, including splice variants, expressed in human islets of Langerhans under control conditions or following exposure to the pro-inflammatory cytokines interleukin-1β (IL-1β) and interferon-γ (IFN-γ). Based on this unique dataset, we examined whether putative candidate genes for T1D, previously identified by GWAS, are expressed in human islets. A total of 29,776 transcripts were identified as expressed in human islets. Expression of around 20% of these transcripts was modified by pro-inflammatory cytokines, including apoptosis- and inflammation-related genes. Chemokines were among the transcripts most modified by cytokines, a finding confirmed at the protein level by ELISA. Interestingly, 35% of the genes expressed in human islets undergo alternative splicing as annotated in RefSeq, and cytokines caused substantial changes in spliced transcripts. Nova1, previously considered a brain-specific regulator of mRNA splicing, is expressed in islets and its knockdown modified splicing. 25/41 of the candidate genes for T1D are expressed in islets, and cytokines modified expression of several of these transcripts. The present study doubles the number of known genes expressed in human islets and shows that cytokines modify alternative splicing in human islet cells. Importantly, it indicates that more than half of the known T1D candidate genes are expressed in human islets. This, and the production of a large number of chemokines and cytokines by cytokine-exposed islets, reinforces the concept of a dialog between pancreatic islets and the immune system in T1D. This dialog is modulated by candidate genes for the disease at both the immune system and beta cell level.
Introduction
Type 1 diabetes (T1D) is an autoimmune disease with a strong genetic component [1]. We have previously proposed that insulitis, the pancreatic islet inflammation present in T1D, results from a “dialog” between immune cells homing into the islets and the target beta cells. Beta cells contribute to this dialog by local release of cytokines and chemokines and by delivering immunogenic signals during the cell death process; this, together with signals generated by invading immune cells, contributes to trigger and amplify (or dampen) insulitis [2]. The amplification or resolution of insulitis, and its progression or not to disease, probably depends on an interplay between environmental triggers, such as dietary components or viral infections, and the patient's genetic background [2], [3], [4] acting at least in part at the pancreatic beta cell level [5], [6], [7]. It is thus important to identify the molecular mechanisms by which immune signals and genetic and/or environmental factors affect beta cell survival and the production of inflammatory mediators such as chemokines and cytokines.
Evaluation of the full transcriptome of beta cells exposed to pro-inflammatory cytokines such as interleukin-1β (IL-1β), tumor necrosis factor-α (TNF-α) and interferon-γ (IFN-γ) provides a snapshot of the responses of these cells under conditions that may prevail in early T1D [2]. Until recently, the only way to analyze large numbers of transcripts was via oligonucleotide array technology. By using this technology we have described expression of nearly 8,000 genes in rat and human islet cells, of which around 20% were modified by cytokines [8], [9], [10]. Arrays, however, can only identify known transcripts due to the need for complementary recognition of probes by the target mRNA. In recent years RNA-sequencing (RNA-seq) has emerged as a new and promising tool for transcriptomic studies. RNA-seq works in an unbiased way, without the need for a priori knowledge of the targets, and shows both high reproducibility and low frequency of false positives [11], [12]. Moreover, RNA-seq is able to identify between 25 and 75% more genes than cDNA microarrays, and allows identification of both whole genes and splice variants [12], [13], [14].
Transcripts of >90% of eukaryotic genes can undergo alternative splicing (AS), i.e. be spliced in more than one way [15]. AS is a basic mechanism for the generation of multiple structurally and functionally distinct mRNAs and protein isoforms from a single gene [15], [16], [17]. It varies in a tissue-specific manner, contributing to tissue specificity [18], [19], [20], and can be modulated by cellular signals such as those provided by pro-inflammatory cytokines [9]. The use of RNA-seq, coupled to dedicated bioinformatic tools, enables the identification of novel splice variants by transcripts with skipped exons, retained introns, alternative start sites, etc [16].
Against this background, we describe here the first RNA-seq analysis of human pancreatic islets. This was done by reverse transcribing and sequencing RNA from human islets obtained from five organ donors, exposed or not to the pro-inflammatory cytokines IL-1β and IFN-γ. The data showed very good internal consistency, and allowed us:
-
To describe the complete human islet cell transcriptome, including splice variants, which provides a novel and valuable resource for future genetic and functional studies;
-
To show that >60% of the candidate genes for T1D, previously believed to be mostly expressed in the immune system [21], are expressed in human islets, and that expression of many of these genes is modified by cytokines;
-
To characterize the impact of an inflammatory challenge, i.e. exposure to pro-inflammatory cytokines, on the human islet transcriptome;
-
To validate some of the key findings obtained by RNA-seq by other methods, e.g. real time RT-PCR, ELISA or histology, in independent samples of human islets and clonal or primary rat beta cells. For some of the novel genes, the use of specific siRNAs allowed clarification of their function in beta cells.
Methods
Ethics statement
Human islet collection and handling were approved by the local Ethical Committee in Pisa, Italy. Wistar rats were used according to the rules of the Belgian Regulations for Animal Care with approval of the Ethical Committee for Animal Experiments of the ULB.
Human islet isolation and culture and rat beta cell culture
Human islet preparations were obtained in collaboration with Pisa University [5], [22], [23], [24]. The donors, aged 68±3 (n = 15), were heart-beating organ donors with no medical history of diabetes or metabolic disorders. Donor information is summarized in Table 1. Preparations 1–5 were used for RNA-seq and preparations 6–15 for independent confirmation of key findings. Isolated islets were used for research when the pancreas was not suitable for clinical transplantation. The human islets were isolated using collagenase digestion and density gradient purification [25]. The islets were cultured in M199 culture medium containing 5.5 mM glucose and shipped within 1–5 days following isolation. Upon arrival, the human islet cells were cultured in Ham's F-10 medium containing 6.1 mM glucose, 10% fetal bovine serum (FBS), 2 mM GlutaMAX, 50 µM 3-isobutyl-1-methylxanthine, 1% BSA, 50 U/ml penicillin and 50 µg/ml streptomycin. The islets were exposed or not to cytokines in the same medium without FBS for 2 days [5], [26]. The following cytokine concentrations were used, based on previous dose-response experiments from our group [5], [27], [28]: recombinant human IL-1β (specific activity 1.8×107 U/mg; a kind gift from C.W. Reinolds, National Cancer Institute, Bethesda, MD, USA) at 50 U/ml; recombinant human IFN-γ (specific activity 2×107 U/mg; R&D Systems, Abingdon, UK) at 1000 U/ml. The evaluation of islet cell purity, i.e. the percentage of beta cells present in the preparations, was done by immunocytochemistry with an anti-insulin antibody (1/1000; Sigma, Bornem, Belgium) and donkey anti-mouse IgG rhodamine (1/200; Lucron Bioproducts, De Pinte, Belgium). Only preparations with more than 40% beta cells were used for the RNA-seq analyses; on average they contained 58% beta cells (Table 1), which is similar to the reported percentage of 54% in isolated human islets [29] and 55% in the human pancreas [30].
For confirmation and mechanistic studies of selected genes, we used the rat insulin-producing INS-1E cell line, kindly provided by C. Wollheim, University of Geneva, Geneva, Switzerland [31]. The cells were maintained in RPMI 1640 medium supplemented with 5% heat-inactivated FBS, 10 mM HEPES, 1 mM Na-pyruvate and 50 µM 2-mercaptoethanol [26]. Cells were exposed to 10 U/ml human IL-1β and 100 U/ml murine IFN-γ (R&D Systems). These cytokine concentrations were selected based on previous dose-response studies [28], [32]; lower cytokine concentrations and shorter time points were used for rodent experiments because rat beta cells are more sensitive than human islets to cytokine damage [33], [34]. Additional confirmation was done in autofluorescence-activated cell sorting (FACS)-purified primary rat beta cells. Pancreatic islets were isolated from adult male Wistar rats (Charles River Laboratories, Brussels, Belgium) and primary beta cells FACS-purified (FACSAria; BD Bioscience, San Jose, CA, USA) and cultured as described [35]. Primary beta cells were transfected with the synthetic double-stranded (ds) RNA polyinosinic-polycytidylic acid (PIC, InvivoGen) as described [6], [7].
RNA sequencing
Five human islet preparations were used for sequencing. Total RNA was isolated using the RNeasy Mini Kit (Qiagen, Venlo, The Netherlands) which favors purification of all RNA molecules longer than 200 nucleotides and sample preparation done as described by the manufacturer (Illumina, Eindhoven, The Netherlands). Briefly, mRNA was purified from two µg total RNA using oligo (dT) beads, before it was fragmented and randomly primed for reverse transcription followed by second-strand synthesis to create ds cDNA fragments. The generated cDNA had undergone paired-end repair to convert overhangs into blunt ends. After 3′-monoadenylation and adaptor ligation, cDNAs were purified on a 2% agarose gel and 200 basepair (bp) products were excised from the gel. Following gel digestion, purified cDNA was amplified by PCR using primers specific for the ligated adaptors. The generated libraries were submitted to quality control with the Agilent bioanalyzer 2100 (Agilent Technologies, Wokingham, UK) before sequencing. The RNA integrity number (RIN) values for all samples were 7.5 and above. 1 µL cDNA was loaded on an Agilent DNA chip (DNA-1000) to verify cDNA quality and quantity. Only libraries reaching satisfactory conditions were used for sequencing, on one sequencing lane of an Illumina Genome Analyzer II system (GAII, Illumina). The raw data generated during the sequencing procedure on the GAII will be deposited in Gene Expression Omnibus (GEO) under submission number GSE35296.
RNA–seq data analysis
Sequencing reads were mapped to the human genome (version GRCh37/hg19) using the program gem-mapper from the GEM suite (http://gemlibrary.sourceforge.net). The GEM mapper reports exhaustively all mappings and split-mappings up to a user-defined amount of mismatches (default 2 mismatches), disregarding presumptive base-calling errors as identified by low associated quality values. Mapped reads were used to quantify transcripts from the RefSeq reference database [36], using the Flux Capacitor approach that deconvolves reads mapping to exonic regions shared by multiple transcripts by optimizing a system of linear equations and thus obtains a number of reads specifically assigned to each alternative spliceform (http://flux.sammeth.net, see [37] for a short description). All genes and transcripts have been assigned a relative coverage rate as measured in RPKM units (“reads per kilobase per million mapped reads”) [38].
Lists of differentially expressed genes and transcripts were generated from the Flux Capacitor output using scripts in Perl or R (see legends to figures and tables).
To define genes up- or downregulated by cytokines, the log2 of the proportion between the sum of the RPKM for all gene transcripts under cytokine condition and the same sum in control condition was taken as measure of change in gene expression. The p-value was obtained by performing a Fisher exact test (number of reads mapped to the gene and number of reads mapped to all other genes in the cytokine condition versus the control condition) and corrected by the Benjamini-Hochberg method (taking for each gene the 5 samples as independent tests). A difference in gene expression was considered significant if the corrected p-value was <0.05. As additional criteria, a gene was considered to be “modified by cytokines” only if its expression changed significantly in one direction - i.e. “up” or “down” - across at least 4 out of 5 islet preparations and no significant change in the opposite direction was observed. In order to quantify cytokine-modified splicing, differences in so-called “splice indices” - the proportion between the RPKM for a transcript and the sum of the RPKM for all the transcripts from the same gene - under cytokine exposure were compared to the control condition. Additionally, a p-value on the significance of changes in splicing patterns was obtained by performing a Fisher exact test (number of reads assigned to a certain transcript after deconvolution versus the number of reads mapped to all other transcripts of the same gene, comparing cytokine with control condition) and was corrected by the Benjamini-Hochberg method (taking for each transcript the 5 samples as independent tests). A change in AS was considered significant if the corrected p-value was <0.05. Consistent with the study of altered gene expression, a transcript was considered as “modified by cytokines” only if its splicing changed significantly in one direction - “up” or “down” - in at least 4 out of 5 islet samples and no sample pair exhibited a significant change in the opposite direction.
Preferential association of the lists of up/downregulated genes/transcripts with molecular and cellular functions and canonical pathways was determined with Benjamini-Hochberg corrected Fisher tests using the Ingenuity Pathway Analysis (IPA, Ingenuity Systems, http://www.ingenuity.com) software. A similar analysis was performed using DAVID (Database for Annotation, Visualization and Integrated Discovery, http://david.abcc.ncifcrf.gov) [39]. While IPA is curated manually, DAVID is generated automatically from 3rd party databases. We used Gene Ontology Biological Process and Molecular Function, KEGG, InterPro and UCSC_TFBS for our DAVID analyses.
Networks of pairwise interactions between proteins, as described in the IntAct database, were obtained from the lists of up/downregulated genes using the PPI_spider [40] from the BioProfiling site (http://www.bioprofiling.de).
We employed an approach similar to the one used to define cytokine-modified genes to compare the untreated control islets to the adipose tissue, colon, kidney, liver and skeletal muscle tissue data available through the Illumina bodyMap2 project (accession number ERP000546 in the European Nucleotide Archive http://www.ebi.ac.uk/ena/data/view/ERP000546). A more detailed comparative analysis between pancreatic islets and other tissues, aiming to detect novel beta cell biomarkers, is under way and will be the subject of a future publication. The RPKM data and lists of cytokine-modified and human islet-specific genes are available in Dataset S1.
Human islet and rat beta cell RNA extraction, RT–PCR, and qRT–PCR
Human islet preparations for validation experiments were from donors other than those used for sequencing (Table 1). In some experiments confirmation was also done in clonal INS-1E and primary rat beta cells, to confirm that gene expression was indeed derived from beta cells (human islets contain different cell types, with beta cells constituting around 60% of the total population in the present samples; Table 1). Poly(A)+ mRNA was isolated using the Dynabeads mRNA DIRECT kit (Invitrogen, Paisley, UK) and reverse transcribed as previously described [26]. Quantitative PCR was performed using the iQ SYBR Green Supermix (BIO-RAD, Nazareth Eke, Belgium) on a LightCycler (Roche Diagnostics, Mannheim, Germany) or iCycler MyiQ Single Color (BIO-RAD) instrument [41], [42]. Data were expressed as number of copies using the standard curve method. Expression values were corrected for the housekeeping gene β-actin and/or glyceraldehyde-3-phosphate dehydrogenase (GAPDH). These housekeeping genes are not modified by pro-inflammatory stimuli under the present experimental conditions [43], [44], [45]. For the evaluation of splice variant expression, conventional PCR was done. Primers were designed for DnaJ homolog subfamily A member 3 (DNAJA3) on exon-spliced junctions between exon 9 and 11 to obtain a product of 267 bp for variant 1 (NM_005147.4) and 150 bp for variant 2 (NM_001135110.1). RNA from INS-1E cells, transfected with a control siRNA (siC) or a siRNA targeting Nova1, was retro-transcribed and the cDNAs amplified with gabrg2 primers. The samples were amplified using BioTAQ Red DNA Polymerase, 10× NH4 reaction buffer, 50 mM MgCl2 and 100 mM dNTP mix (BioLine, London, UK) in a Thermal Cycler (Applied Biosystems) using the following conditions: after 8 min of denaturation at 95°C, samples were run for 32–35 cycles consisting of 1 min at 95°C, 45 sec at 60°C and 1 min at 72°C. The final step was 5 min at 72°C. PCR products were visualized on 2.3% agarose gel, stained with SYBR Safe gel stain (Invitrogen). Primers used for qRT- and RT-PCR are listed in Table S1.
RNA interference
For RNA interference in rat beta cells, the following siRNAs were used: smart pool targeting MDA5 (reference 105259, Thermo Scientific), siBCL2A1 CAGGGAAGAUCUGGGAAAUGCUCUU, smart pool targeting BCL2A1 (reference 170929, Thermo Scientific), siNova1 stealth UUAGCAUGUCCUAAUAGCCCUGCGG (Invitrogen) and Allstars Negative Control siRNA (Qiagen, Venlo, the Netherlands). Cells were transfected with a mix of 30 nM of siRNA and Lipofectamine RNAiMAX (Invitrogen) diluted in Opti-MEM I (Invitrogen) as described [5]. The transfection efficiency was >90% [5], [46]. After overnight transfection the cells were cultured for 48 h before being retrieved for evaluation of RNA and protein expression.
Western blot and chemokine and cytokine ELISA
For Western blotting, equal amounts of proteins were loaded in 12% SDS-PAGE. Immunoblot analysis was performed using goat anti-Nova1 (0.03 µg/ml; Abcam, Cambridge, UK) and mouse anti-α-tubulin (1∶5000; Sigma) antibodies. The proteins were detected using horseradish peroxidase-conjugated secondary antibody (1∶5000; Santa Cruz Biotechnology) and chemiluminescence Supersignal (Pierce). Densitometric analysis was performed using analysis software Aida1D (Fujifilm, London, UK) and data were normalized for α-tubulin.
Release of the human chemokines CXCL1 (Gro-α), CXCL9 (Mig), CXCL10 (IP-10), CXCL11 (Itac), CCL2 (MCP-1), CCL3 (Mip-1-α), CCL5 (Rantes) and the cytokines IL-6 and IL-8 was measured in culture medium of control and cytokine-exposed human islets using a Custom Multi-Analyte ELISArray kit (SABiosciences, Frederick, MD, USA). Samples were processed following the manufacturer's instructions. This is a semi-quantitative assay that does not include a standard curve. Absorbance at 450 nm was measured, corrected by readings at 570 nm, normalized to the geometric mean of β-actin and GAPDH expression and expressed as arbitrary units.
Immunofluorescence
Human pancreatic tissue obtained from biopsies or organ donors were fixed in formaldehyde and embedded in paraffin. Sections were stained for double immunofluorescence with rabbit anti-Nova1 (1∶500; Merck-Millipore, Overijse, Belgium) and guinea pig anti-insulin (I2018, 1∶2000; Sigma) or mouse anti-glucagon antibodies using FITC and Cy3 as fluorochromes, respectively. The samples were analyzed by inverted fluorescence microscopy and images captured with Axiocam (Zeiss).
Assessment of apoptosis
The percentage of apoptotic cells was determined by two observers (one being blind to sample identity), after staining with the DNA-binding dyes propidium iodide and Hoechst 33342 (Sigma-Aldrich) as previously described [47]. At least 500 cells were counted per condition, with an agreement between findings obtained by the two observers of >90%.
Statistical analysis
Data for the confirmation experiments are presented as means ± SEM. Comparisons were performed by paired two-tailed Student's t-test or Mann Whitney test as indicated in the figure legends. A p-value≤0.05 was considered statistically significant. The statistical analysis of the RNA-seq data is described above.
Results
Sequencing of human islets and analysis of transcripts
RNA-seq data were obtained from 5 human islet preparations (Table 1) cultured under control condition or following a 48-h exposure to the cytokines IL-1β+IFN-γ. Each of these preparations was sequenced on a single lane of an Illumina GAII sequencer, with 10–51 million reads for control and 35–62 million reads for cytokine-treated islets. This provides sufficient sequencing depth to quantify gene expression and detect rare transcripts as previously shown [16]. The 51 nucleotide paired-end reads were mapped to the human genome (version hg19) using GEM software. Taking this approach, we were able to map on average 83% of the raw reads. GEM can report multiple mappings for a single read and we observed on average a redundancy (mappings to reads ratio) of 1.5 (Table S2).
Reads that align with exons or with overlapping exon junctions can be used to evaluate the levels of splicing. We used the Flux Capacitor software, which in brief takes as input a list of reads mapped to the genome and a list of transcript annotations, and subsequently produces a list of reads that are uniquely assigned to one of the transcripts. As reference transcript annotation, we employed the 34,102 annotated human mRNA and ncRNA sequences from RefSeq [48]. In a first step, the program interprets the mate information of mappings and filters off mappings that do not pair properly within the boundaries of annotated transcripts. For about half of the originally sequenced reads a mate in correct orientation and within exon boundaries of the annotated RefSeq transcripts could be identified, with only spurious redundancy (<1.01). The 34,102 transcripts from RefSeq correspond to 22,205 genes, and islets cultured under control condition were found to express a median of 17,787 genes, with numbers varying with sequencing depth (Table S2). Of these, 15,212 genes were expressed in all individuals while 3,841 genes were expressed in some but not all. 5,408 genes expressed in all individuals have AS annotated in RefSeq (see below).
Analyses of the qualitative agreement of expression levels between the individual islet preparations using Pearson correlation coefficients (PCC) indicated a high correlation (0.95) (Figure 1A). As gene expression follows Zipf's law [49], corresponding quantification values have been power-law normalized to meet the prerequisite for correlation studies [50]. For each sample-pair, the corresponding PCC provides a numerical condensation of the similarity between gene expression profiles: a PCC of 1.0 represents sample-pairs where all expression tuples fall along a line, whereas a PCC of 0 is assigned to sample pairs that do not exhibit linear correlation. For the purpose of comparison, 5 tissues from the Illumina Human Body Map project, i.e. colon, adipose, kidney, liver and skeletal muscle, have been subjected to an analogous procedure. The similarity in terms of correlation among gene expression levels was significantly higher between the islet samples (0.90–0.96) than in comparison with the 5 other tissues (0.53–0.88). In line with these observations, a heatmap with complete linkage as clustering function indicates that the 5 islet preparations clustered together, as compared to the other tissues (Figure 1B). It cannot be excluded that islet culture affects the human islet transcriptome, although in other studies differential gene expression between diabetic and non-diabetic individuals was maintained after culture [25], [51], [52].
For internal methodological validation, we selected 4 genes for confirmation by qRT-PCR in the same samples used for RNA-seq. The gene expression data using these two methods were essentially superposable (Figure S1).
The validation steps described above, including comparison between islet samples and against 5 other tissues, and the validation using qRT-PCR in the same samples, indicate that the RNA-seq of human islets provided reliable and reproducible data, as has been described for other tissues [11], [16], [38], [53], enabling us to proceed with the analyses described below.
Expression of candidate genes for type 1 diabetes in pancreatic islets
Based on the datasets above, we examined whether candidate genes for T1D, previously identified by genome-wide association studies (GWAS) [54], [55], are expressed in human islets. We considered genes as “expressed” with a median RPKM >1. Out of 41 candidate genes, 25 (i.e. 61%) were clearly expressed in human islets (Figure 2A and Table S3). We followed this up by functional studies in insulin-producing INS-1E cells and purified rat beta cells, to confirm gene expression and query the relevance of these genes at the beta cell level. We have previously shown that 2 of these genes, namely IFIH1/MDA5 and PTPN2, are expressed in pancreatic beta cells and regulate respectively local inflammation [6] and apoptosis [5], [56]. Pro-inflammatory cytokines and dsRNA, a by-product of viral infections, modulate expression of these 2 genes, indicating crosstalk between T1D candidate genes and environmental factors and local inflammation [5], [6], [56]. Indeed, knockdown of IFIH1/MDA5 in rat beta cells reduced the chemokine and cytokine expression induced by a 48-h exposure to PIC, a synthetic dsRNA (Figure S2). We now confirm in clonal INS-1E cells expression of an additional candidate gene, namely SH2B3 (Figure 2B), and its induction by the cytokines IL-1β+IFN-γ in a time-course study.
It is commonly thought that antioxidative defense mechanisms of pancreatic beta cells are very low, rendering the cells vulnerable to reactive oxygen species which contribute to the pathogenesis of diabetes (reviewed in [57]). This seems to be the case for rat beta cells [58], but we have previously shown that human beta cells are 5–10-fold more resistant than mouse or rat islets to oxidative stress generated by agents such as alloxan [33]. We compared expression of several free radical scavenging enzymes in human islets against 5 other tissues (Table S4). Human islets have robust expression of several of these enzymes, including a marked expression of catalase (median RPKM 26) and SOD2 (median RPKM 388). In line with these findings, we have previously observed that human islets have several-fold higher expression of antioxidant enzymes than rodent islets [59]. Expression levels in islets compared to liver were lower for 3 antioxidant enzymes, similar for 3 and significantly higher for 4 enzymes, suggesting that human islets, as opposed to rodent islets, may have a fair antioxidant capacity.
Analysis of cytokine-modified genes
From the 19,621 genes detected as “present” by the RNA-seq, a total of 3068 (16%) were significantly modified by a 48-h exposure to the pro-inflammatory cytokines IL-1β+IFN-γ. From these, 1416 and 1652 were respectively up- and downregulated. The complete list of cytokine-modulated genes is accessible at http://lmedex.ulb.ac.be/data.php; password will be provided upon request. These genes were manually curated (by DLE; see selected cytokine-modified genes in Table S5) or analyzed in a non-biased way using IPA (Figure 3). Table S5 indicates that many key beta cell functions were modified by cytokines, including glucose and lipid metabolism, protein synthesis and translation, kinases and phosphatases and transcription factors. The most important responses, however, were those related to inflammation, innate immune response and apoptosis. Thus, there was massive up-regulation of the expression of a large number of genes encoding chemokines and cytokines, of genes involved in IFN-γ signaling and NF-κB regulation, proteasome/antigen presentation, and other innate immune response/pro-inflammatory components. There was also up-regulation of many genes involved in apoptosis, free radical scavenging and DNA damage response (Table S5). These observations were supported by IPA, which showed that up-regulated genes belong prominently to the functions “Cell Death” and “Cellular Movement” (actually mainly chemokines) (Figure 3A). In the IPA “Diseases and Disorders” analysis (not shown) modified genes clustered in “Inflammatory Response”. As shown in Figure 3B, IPA canonical pathways indicated that the highest p-value was related to “Acute Phase Response Signaling”. Interestingly, among the top canonical pathways we also found several other inflammation-related headlines, such as “Role of macrophages…”, “Dendritic cell…”, “Altered T and B cell signaling”, “IL-17 signaling” and, reassuringly, “Type 1 diabetes mellitus signaling”. The fact that 6 of the top canonical pathways were related to IL-17 is of particular interest given that IL-17 signaling may play a direct role in beta cell apoptosis in human T1D [43].
The analysis using IPA was validated by a separate analysis using the public tool DAVID, which relies on copies of various public databases. A term enrichment analysis against Gene Ontology, KEGG (metabolic and regulatory pathways) and InterPro (protein conserved motifs) showed that the up-regulated genes were preferentially associated with immune response, apoptosis, cytokines and other terms related to inflammatory stress (Figure S3A–S3D). Noteworthy is that term enrichment analysis against UCSC_TFBS showed genes with potential binding sites for the transcription factors NF-κB, AP-1 (Jun) and BACH2 (not shown). Protein-protein interactions among the up-regulated genes were examined using the BioProfiling tool, which relies on the IntAct database (Figure S4). It shows several interactions related to inflammatory response and antigen processing and presentation, with a clear role for members of the NF-κB and STAT families. The observations by RNA-seq of cytokine-induced chemokines (Table S5) are in line with our previous observations using array analysis of human islets exposed to viral infection or pro-inflammatory cytokines [60] or qRT-PCR of human and mouse islets exposed to cytokines or isolated from pre-diabetic mice [42]. The RNA-seq findings were confirmed at the protein level by ELISA for nearly all chemokines studied (Figure 4), indicating that many of the observed gene expression changes are translated to functional proteins with potential relevance for the early pathogenesis of T1D.
“Molecular and Cellular Function” IPA showed that downregulated genes belonged to “Cell Morphology”, “Assembly and Organization”, “Growth and Proliferation” and “Movement” (Figure 3C) and amino acid metabolism in the IPA “Canonical Pathways” (Figure 3D). A DAVID term enrichment analysis produced similar results (Figure S3E–S3H).
We compared the presently observed cytokine-modified genes in human islets against our recently published array data in cytokine-exposed INS-1E cells [61]. For this comparison, we used genes with homology between human and rat, and probes present on the Affymetrix GeneChip Rat Genome 230 2.0 array. This selection encompassed 790 and 874 genes, considered as up- or down-regulated by RNA-seq, respectively. Of these, 53% and 50% were detected as respectively “up-” and “down-regulated” in cytokine-treated INS-1E cells (data not shown). When we focused on some of the most relevant cytokine-modulated genes (Table S5) the overlap was even higher between human islet and INS-1E genes, with respectively 76% and 63% of the NF-κB/other transcription factors and chemokines showing a similar variation. Considering the issues of species differences (human vs rat), differences between primary/clonal cells (islets vs INS-1E cells), methodological differences (RNA-seq vs array analysis to asses RNA expression) and timing of exposure to cytokines (48 h for human islets and 12–24 h for INS-1E cells), the observed correlation (50–53%) between genes expressed in human islets and INS-1E cells is reasonable, and suggest that many of the presently observed cytokine-modified genes are expressed in beta cells.
To further confirm expression of some of the cytokine-modified genes, we used independent samples for qRT-PCR evaluation. We selected genes potentially involved in apoptosis, namely Bcl-2 related protein A1 (BCL2A1) and Bcl-2 modifying factor (BMF) (Figure 5). In line with the RNA-seq data (Table S5), cytokines respectively increased and decreased expression of BCL2A1 and BMF (Figure 5A and 5B). In clonal INS-1E cells, BCL2A1 expression was induced by IL-1β+IFN-γ in a time-course study (Figure 5C). Efficient knockdown of BCL2A1 using two different siRNAs (Figure 5D) amplified cytokine-induced apoptosis (Figure 5E), demonstrating the anti-apoptotic role of BCL2A1 in beta cells. Another Bcl-2 family member that was up-regulated by cytokines in the RNA-seq is BBC3 (PUMA, Table S5). BBC3 was recently shown to be cytokine-induced at the mRNA and protein level and pro-apoptotic in beta cells [23].
Of interest, expression of several of the putative candidate genes for T1D (Figure 2) was modified by 48-h exposure to cytokines. Besides the ones already discussed above (PTPN2, IFIH1 and SH2B3), STAT-4, GLIS-3, CD55, RASGRP1, SKAP2 and a large number of HLA-related genes tended to increase expression following cytokine exposure (Table S5).
Splice variants in human islets and their regulation by Nova1 and pro-inflammatory cytokines
Within the 5 islet samples, we found evidence for 87.3% of the islet-expressed genes that have multiple RefSeq transcripts annotated to express more than one spliceform. The complete list of these transcripts is available online at http://lmedex.ulb.ac.be/data.php; password will be provided upon request. Since there is no available information on the regulation of splicing in human or rat islet cells, we examined the expression in human islets of 224 genes previously identified as splicing factors in other human tissues [62] and found that most of them are expressed in islets, and 69 significantly more than in at least 4 out of 5 selected background tissues (adipose tissue, colon, kidney, liver and skeletal muscle, data not shown). We detected expression of several so-called “neuron-specific” splicing factors, including Nova1. Nova1 participates in the splicing of several genes implicated in neuronal function and development [63], [64], and was previously detected by microarray profiling of human islets [65]. We confirmed by qRT-PCR that Nova1 is indeed well expressed in human islets, at levels comparable to brain and higher than in liver, spleen, colon and lung (Figure 6A). Expression of Nova1 at the protein level was confirmed in insulin-positive beta and glucagon-positive alpha cells in human pancreatic sections, while there was little or no staining in the exocrine pancreas (Figure 6B and 6C). To explore the splicing function of Nova1 in beta cells, the gene was knocked down by a specific siRNA in insulin-producing INS-1E cells, leading to a nearly 60% decrease in Nova1 mRNA and protein expression (Figure 6D and 6E). To test the functional impact of Nova1 knockdown, we evaluated the expression of splice variants of gamma-aminobutyric acid A receptor, gamma 2 (Gabrg2). Nova1 was previously shown to cause exon 9 inclusion in Gabrg2 transcripts in mouse brain [66]. Primers were designed on the flanking regions of this exon (Figure 6F) to differentiate between the long transcript variant with exon 9, the short variant without exon 9 and an intermediate undefined variant [67]. Knockdown of Nova1 modified the splicing pattern of the gabrg2 transcripts generating more of the short variant (Figure 6G), suggesting a functional role for this splicing factor in beta cells.
Exposure of human islets to the cytokines IL-1β+IFN-γ induced modifications in the splicing of 548 genes; of these 425 and 433 splice variants were respectively up- and downregulated by the cytokines, as evaluated by a conservative assessment (see Methods). IPA of transcripts exhibiting cytokine-modified AS indicates that a large number of transcripts were related to “Cell Death” or “Cellular Growth and Proliferation” (Figure 7A) and canonical pathways of T and B cells and PKA, calcium, AMPK and p53 signaling (Figure 7B). A DAVID term enrichment analysis yielded among the top terms “cell death” and “apoptosis” (not shown).
To validate the RNA-seq analysis, we selected DNAJA3 for PCR confirmation in independent samples. DNAJA3 is related to “Cell Death” and its variants 1 and 2 were respectively down- and up-regulated by cytokines in 5 out of 5 islet samples. By RT-PCR, the cytokines IL-1β and IFN-γ increased variant 2 expression in 3 independent human islet preparations (Figure S5).
Discussion
We presently describe the first global sequencing of RNAs expressed in human islets of Langerhans. The analysis identified 15,200 genes expressed in the five independent preparations, increasing by >2-fold the known expressed genes in human islets. There was a high correlation between the islet samples (0.90–0.96), clearly higher than the correlation observed between islets and five other tissues (0.53–0.88) used for external comparison. This, and the fact that around 20 genes identified as expressed and/or modified by cytokines in the present analysis were confirmed at the RNA and/or protein expression level by other methods, supports the reliability of the present observations. This is in line with previous studies in other tissues indicating that RNA-seq is a reliable and reproducible method to evaluate RNA expression [11], [16], [38], [53].
The human islets used in this analysis contained 58% beta cells on average (Table 1), and the transcriptome includes RNAs from non-beta endocrine cells, mostly alpha and delta cells [51], and ductal cells. The comparison against INS-1E cells suggests, nonetheless, that at least half of the presently identified cytokine-modified genes are expressed in beta cells.
Use of GWAS has revealed more than 40 loci containing putative genetic contributors to the pathogenesis of T1D [54], [55]; this number was further increased by a recent genome-wide meta-analysis of six diabetes cohorts [68]. While in T2D most candidate genes impact more on islet function than on insulin resistance and are hence considered to regulate beta cell function and development [69], [70], it is usually assumed that in T1D most if not all candidate genes modulate the immune system (reviewed in [21]). In this conventional view beta cells are regarded as “passive victims” of a process that starts and is regulated elsewhere. By using the presently generated datasets, we observed that 61% of the candidate genes for T1D are consistently expressed in human pancreatic islets. Furthermore, the present and previous observations [5], [6], [56] indicate that expression of many of these genes change following exposure to pro-inflammatory cytokines or dsRNA (a by-product of virus infection), agents that may contribute to triggering of T1D [2]. For at least two of these genes, namely IFIH1/MDA5 [6] (present data) and PTPN2 [5], [6], [56], there is experimental evidence pointing to their respective roles in production of chemokines/cytokines and beta cell apoptosis.
These observations are in line with the present analysis of gene expression in cytokine-treated human islets. Of note, only one time point (48 h cytokine exposure) was examined here, providing a snapshot of dynamic regulation of gene expression. It is conceivable that relevant cytokine-modulated genes at other time points were missed in the present analysis. Cytokines modified expression of 3,000 genes, mostly related to inflammation, innate immune response and apoptosis. Key chemokines and cytokines were among the most up-regulated genes in human islets, a finding confirmed at the protein level for CCL2, CCL5, CCL3, CXCL9, CXCL10, CXCL11, IL-6 and IL-8. This is in good agreement with findings in diabetes-prone NOD mice, where increased expression of CCL2, CXCL10 and other chemokines/cytokines are observed in the pre-diabetic period [42], [71], [72]. CCL2 and CXCL10 attract macrophages, and may contribute to the recruitment of immune cells during the early stages of insulitis, as suggested by the observation that transgenic expression of CCL2 in beta cells causes insulitis and diabetes [72]. Some of these observations have been recently confirmed in histological material from T1D patients. Thus, it was observed that pancreatic beta cells from islets affected by insulitis express CXCL10, while the infiltrating T cells express CXCR3, the receptor of CXCL10 [73], [74]. Islet cells themselves are probably an important source of chemokine production during inflammation, as suggested by the present findings. That chemokines are indeed produced by beta cells is supported by the observations that FACS-purified rat beta cells (>90% pure) or clonal rat beta cells (INS-1E cells) exposed to IL-1β+IFN-γ, or to dsRNA, show increased expression of mRNAs encoding CCL2, CXCL10, CCL20, CX3CL1 and IL-15, among others [9], [44], [61], [75]. This is confirmed by histology of pancreatic samples, showing expression of chemokines by beta cells [73], [74], [76].
The findings described above support the concept of a “dialogue” between beta cells and the invading macrophages and T cells in the course of insulitis, rather than a “monolog” where all action takes place at the level of the immune system and beta cells are no more than passive victims. Thus, activated mononuclear cells produce cytokines such as IFN-γ, IL-1β and TNF-α, triggering the release of chemokines and stimulatory cytokines by the beta cells. This, together with beta cell death and the putative presentation of neoantigens secondary to modified AS and up-regulation of the machinery for antigen presentation, will attract more mononuclear cells that also release multiple cytokines and chemokines, in a process modulated by candidate genes that are expressed and act at both the immune system and beta cell levels, as shown for MDA5 and PTPN2, among others.
One of the most deleterious consequences of islet inflammation is the progressive loss of pancreatic beta cells via apoptosis [2]. We presently observed modulation of the expression of several apoptosis-related genes in human islets exposed to cytokines. One of them, the anti-apoptotic Bcl-2 family member BCL2A1 [77], [78], was confirmed by qRT-PCR in both independent human islet preparations and in clonal rat insulin-producing INS-1E cells. Knock down of BCL2A1 by a specific siRNA augmented both basal and cytokine-induced apoptosis, confirming the relevant function of this protein in protecting beta cells against apoptosis (present data). Cytokine-induced expression of BCL2A1 in human islets has been previously observed by array analysis [60], [79], but the function of this gene in beta cells remained to be clarified. Of interest, BCL2A1 inhibits apoptosis induced by, among others, the BH3 only protein Bim [80], [81]. Bim was recently shown to be a crucial pro-apoptotic signal following inhibition of the candidate gene PTPN2 [56], a gene also detected in the present RNAseq.
We presently report another level of molecular regulation of beta cell function, namely AS. Interestingly, AS is modified by cytokine exposure as suggested by the present findings in human islets and previous observations from our group based on exon array analysis in rat beta cells [9]. Regulation of splicing in other tissues involves the cooperation between SR, hnRNPs proteins and several other tissue-specific regulators of splicing such as neuron-specific Nova or the neural/muscle-enriched Fox proteins [82], [83]. The well-characterized Nova proteins regulate numerous splicing events in the central nervous system [64], [84], and the present findings show that Nova1 is expressed in beta cells and affects splicing of at least one target gene, namely Gabrg2. Of interest, several of the known Nova target genes in brain are also expressed in beta cells, including neuroligin and neurexin family members, inhibitory synapse-associated neuroligin and neurexin binding partners [64], [85]. These findings are in line with previous observations that beta cells share expression of a large number of genes and proteins with the central nervous system [86], [87], [88]. This opens a new field of research, and new experiments are now required to determine how AS is regulated in beta cells, and how cytokines modify this process.
In conclusion, the present study identifies most of the transcripts present in human islets of Langerhans, providing a valuable dataset for future genetic and functional studies in pancreatic beta cells. It also shows that pro-inflammatory cytokines modify AS and the expression of nearly 20% of the genes expressed in human islet cells. Importantly, the present observations indicate that >60% of the known candidate genes for T1D are expressed in human islets. This, taken together with the cytokine-induced expression of a large number of chemokines and cytokines in human islets, reinforces the concept of a dialog between pancreatic islets and the immune system, which might be crucial for triggering insulitis and eventual progression to diabetes. The present study identifies a large number of the words used by pancreatic islets in this dialog, and points to candidate genes for T1D as one of the writers of the beta cell speeches.
Supporting Information
Zdroje
1. ToddJA 2010 Etiology of type 1 diabetes. Immunity 32 457 467
2. EizirikDLColliMLOrtisF 2009 The role of inflammation in insulitis and β-cell loss in type 1 diabetes. Nat Rev Endocrinol 5 219 226
3. NorrisJM 2010 Infant and childhood diet and type 1 diabetes risk: recent advances and prospects. Curr Diab Rep 10 345 349
4. BoettlerTvon HerrathM 2011 Protection against or triggering of Type 1 diabetes? Different roles for viral infections. Expert Rev Clin Immunol 7 45 53
5. MooreFColliMLCnopMEsteveMICardozoAK 2009 PTPN2, a candidate gene for type 1 diabetes, modulates interferon-γ-induced pancreatic β-cell apoptosis. Diabetes 58 1283 1291
6. ColliMLMooreFGurzovENOrtisFEizirikDL 2010 MDA5 and PTPN2, two candidate genes for type 1 diabetes, modify pancreatic β-cell responses to the viral by-product double-stranded RNA. Hum Mol Genet 19 135 146
7. ColliMLNogueiraTCAllagnatFCunhaDAGurzovEN 2011 Exposure to the viral by-product dsRNA or Coxsackievirus B5 triggers pancreatic β cell apoptosis via a Bim/Mcl-1 imbalance. PLoS Pathog 7 e1002267 doi:10.1371/journal.ppat.1002267
8. EizirikDLMooreFFlamezDOrtisF 2008 Use of a systems biology approach to understand pancreatic β-cell death in Type 1 diabetes. Biochem Soc Trans 36 321 327
9. OrtisFNaamaneNFlamezDLadriereLMooreF 2010 Cytokines interleukin-1β and tumor necrosis factor-α regulate different transcriptional and alternative splicing networks in primary β-cells. Diabetes 59 358 374
10. GurzovENBarthsonJMarhfourIOrtisFNaamaneN 2011 Pancreatic β-cells activate a JunB/ATF3-dependent survival pathway during inflammation. Oncogene Epub ahead of print
11. RichardHSchulzMHSultanMNurnbergerASchrinnerS 2010 Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments. Nucleic Acids Res 38 e112
12. SultanMSchulzMHRichardHMagenAKlingenhoffA 2008 A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science 321 956 960
13. TangFBarbacioruCNordmanELiBXuN 2010 RNA-Seq analysis to capture the transcriptome landscape of a single cell. Nat Protoc 5 516 535
14. van BakelHNislowCBlencoweBJHughesTR 2010 Most “dark matter” transcripts are associated with known genes. PLoS Biol 8 e1000371 doi:10.1371/journal.pbio.1000371
15. PanQShaiOLeeLJFreyBJBlencoweBJ 2008 Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet 40 1413 1415
16. WangZGersteinMSnyderM 2009 RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10 57 63
17. ModrekBLeeC 2002 A genomic view of alternative splicing. Nat Genet 30 13 19
18. SugnetCWSrinivasanKClarkTAO'BrienGClineMS 2006 Unusual intron conservation near tissue-regulated exons found by splicing microarrays. PLoS Comput Biol 2 e4 doi:10.1371/journal.pcbi.0020004
19. LeKMitsourasKRoyMWangQXuQ 2004 Detecting tissue-specific regulation of alternative splicing as a qualitative change in microarray data. Nucleic Acids Res 32 e180
20. CooperTAWanLDreyfussG 2009 RNA and disease. Cell 136 777 793
21. ConcannonPRichSSNepomGT 2009 Genetics of type 1A diabetes. N Engl J Med 360 1646 1654
22. GurzovENOrtisFCunhaDAGossetGLiM 2009 Signaling by IL-1β+IFN-γ and ER stress converge on DP5/Hrk activation: a novel mechanism for pancreatic β-cell apoptosis. Cell Death Differ 16 1539 1550
23. GurzovENGermanoCMCunhaDAOrtisFVanderwindenJM 2010 p53 up-regulated modulator of apoptosis (PUMA) activation contributes to pancreatic β-cell apoptosis induced by proinflammatory cytokines and endoplasmic reticulum stress. J Biol Chem 285 19910 19920
24. CunhaDALadrièreLOrtisFIgoillo-EsteveMGurzovEN 2009 Glucagon-like peptide-1 agonists protect pancreatic β-cells from lipotoxic endoplasmic reticulum stress through upregulation of BiP and JunB. Diabetes 58 2851 2862
25. MarchettiPBuglianiMLupiRMarselliLMasiniM 2007 The endoplasmic reticulum in pancreatic beta cells of type 2 diabetes patients. Diabetologia 50 2486 2494
26. CardozoAKOrtisFStorlingJFengYMRasschaertJ 2005 Cytokines downregulate the sarcoendoplasmic reticulum pump Ca2+ ATPase 2b and deplete endoplasmic reticulum Ca2+, leading to induction of endoplasmic reticulum stress in pancreatic β-cells. Diabetes 54 452 461
27. OrtisFCardozoAKCrispimDStorlingJMandrup-PoulsenT 2006 Cytokine-induced proapoptotic gene expression in insulin-producing cells is related to rapid, sustained, and nonoscillatory nuclear factor-κB activation. Mol Endocrinol 20 1867 1879
28. EizirikDLMandrup-PoulsenT 2001 A choice of death - the signal-transduction of immune-mediated beta cell apoptosis. Diabetologia 44 2115 2133
29. BrissovaMFowlerMJNicholsonWEChuAHirshbergB 2005 Assessment of human pancreatic islet architecture and composition by laser scanning confocal microscopy. J Histochem Cytochem 53 1087 1097
30. CabreraOBermanDMKenyonNSRicordiCBerggrenPO 2006 The unique cytoarchitecture of human pancreatic islets has implications for islet cell function. Proc Natl Acad Sci U S A 103 2334 2339
31. AsfariMJanjicDMedaPLiGHalbanPA 1992 Establishment of 2-mercaptoethanol-dependent differentiated insulin-secreting cell lines. Endocrinology 130 167 178
32. KutluBDarvilleMICardozoAKEizirikDL 2003 Molecular regulation of monocyte chemoattractant protein-1 expression in pancreatic β-cells. Diabetes 52 348 355
33. EizirikDLPipeleersDGLingZWelshNHellerstromC 1994 Major species differences between humans and rodents in the susceptibility to pancreatic β-cell injury. Proc Natl Acad Sci U S A 91 9253 9256
34. EizirikDLSandlerSWelshNCetkovic-CvrljeMNiemanA 1994 Cytokines suppress human islet function irrespective of their effects on nitric oxide generation. J Clin Invest 93 1968 1974
35. RasschaertJLadriereLUrbainMDogusanZKatabuaB 2005 Toll-like receptor 3 and STAT-1 contribute to double-stranded RNA+interferon-γ-induced apoptosis in primary pancreatic β-cells. J Biol Chem 280 33984 33991
36. PruittKDTatusovaTMaglottDR 2007 NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 35 D61 65
37. MontgomerySBSammethMGutierrez-ArcelusMLachRPIngleC 2010 Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 464 773 777
38. MortazaviAWilliamsBAMcCueKSchaefferLWoldB 2008 Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 5 621 628
39. DennisGJrShermanBTHosackDAYangJGaoW 2003 DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol 4 P3
40. AntonovAVDietmannSRodchenkovIMewesHW 2009 PPI spider: a tool for the interpretation of proteomics data in the context of protein-protein interaction networks. Proteomics 9 2740 2749
41. KharroubiIRasschaertJEizirikDLCnopM 2003 Expression of adiponectin receptors in pancreatic β cells. Biochem Biophys Res Commun 312 1118 1122
42. CardozoAKProostPGysemansCChenMCMathieuC 2003 IL-1β and IFN-γ induce the expression of diverse chemokines and IL-15 in human and rat pancreatic islet cells, and in islets from pre-diabetic NOD mice. Diabetologia 46 255 266
43. ArifSMooreFMarksKBouckenoogheTDayanCM 2011 Peripheral and islet interleukin-17 pathway activation characterizes human autoimmune diabetes and promotes cytokine-mediated beta-cell death. Diabetes 60 2112 2119
44. CardozoAKKruhofferMLeemanROrntoftTEizirikDL 2001 Identification of novel cytokine-induced genes in pancreatic β-cells by high-density oligonucleotide arrays. Diabetes 50 909 920
45. KutluBCardozoAKDarvilleMIKruhofferMMagnussonN 2003 Discovery of gene networks regulating cytokine-induced dysfunction and apoptosis in insulin-producing INS-1 cells. Diabetes 52 2701 2719
46. CunhaDAHekermanPLadrièreLBazarra-CastroAOrtisF 2008 Initiation and execution of lipotoxic ER stress in pancreatic β-cells. J Cell Sci 121 2308 2318
47. CnopMLadrièreLHekermanPOrtisFCardozoAK 2007 Selective inhibition of eukaryotic translation initiation factor 2α dephosphorylation potentiates fatty acid-induced endoplasmic reticulum stress and causes pancreatic β-cell dysfunction and apoptosis. J Biol Chem 282 3989 3997
48. PruittKDTatusovaTKlimkeWMaglottDR 2009 NCBI Reference Sequences: current status, policy and new initiatives. Nucleic Acids Res 37 D32 36
49. FurusawaCKanekoK 2003 Zipf's law in gene expression. Phys Rev Lett 90 088102
50. BengtssonMStahlbergARorsmanPKubistaM 2005 Gene expression profiling in single cells from the pancreatic islets of Langerhans reveals lognormal distribution of mRNA levels. Genome Res 15 1388 1392
51. MarchettiPDel GuerraSMarselliLLupiRMasiniM 2004 Pancreatic islets from type 2 diabetic patients have functional defects and increased apoptosis that are ameliorated by metformin. J Clin Endocrinol Metab 89 5535 5541
52. Del GuerraSLupiRMarselliLMasiniMBuglianiM 2005 Functional and molecular defects of pancreatic islets in human type 2 diabetes. Diabetes 54 727 735
53. NagalakshmiUWangZWaernKShouCRahaD 2008 The transcriptional landscape of the yeast genome defined by RNA sequencing. Science 320 1344 1349
54. BarrettJCClaytonDGConcannonPAkolkarBCooperJD 2009 Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet 41 703 707
55. ToddJAWalkerNMCooperJDSmythDJDownesK 2007 Robust associations of four new chromosome regions from genome-wide analyses of type 1 diabetes. Nat Genet 39 857 864
56. SantinIMooreFColliMLGurzovENMarselliL 2011 PTPN2, a candidate gene for type 1 diabetes, modulates pancreatic β-cell apoptosis via regulation of the BH3-only protein Bim. Diabetes 60 3279 3288
57. LenzenS 2008 Oxidative stress: the vulnerable β-cell. Biochem Soc Trans 36 343 347
58. TiedgeMLortzSDrinkgernJLenzenS 1997 Relation between antioxidant enzyme gene expression and antioxidative defense status of insulin-producing cells. Diabetes 46 1733 1742
59. WelshNMargulisBBorgLAWiklundHJSaldeenJ 1995 Differences in the expression of heat-shock proteins and antioxidant enzymes between human and rodent pancreatic islets: implications for the pathogenesis of insulin-dependent diabetes mellitus. Mol Med 1 806 820
60. YlipaastoPKutluBRasilainenSRasschaertJSalmelaK 2005 Global profiling of coxsackievirus- and cytokine-induced gene expression in human pancreatic islets. Diabetologia 48 1510 1522
61. MooreFNaamaneNColliMLBouckenoogheTOrtisF 2011 STAT1 is a master regulator of pancreatic β-cell apoptosis and islet inflammation. J Biol Chem 286 929 941
62. GrossoARGomesAQBarbosa-MoraisNLCaldeiraSThorneNP 2008 Tissue-specific splicing factor gene expression signatures. Nucleic Acids Res 36 4823 4832
63. ZhangCFriasMAMeleARuggiuMEomT 2010 Integrative modeling defines the Nova splicing-regulatory network and its combinatorial controls. Science 329 439 443
64. JensenKBDredgeBKStefaniGZhongRBuckanovichRJ 2000 Nova-1 regulates neuron-specific alternative splicing and is essential for neuronal viability. Neuron 25 359 371
65. KutluBBurdickDBaxterDRasschaertJFlamezD 2009 Detailed transcriptome atlas of the pancreatic β cell. BMC Med Genomics 2 3
66. DredgeBKDarnellRB 2003 Nova regulates GABA(A) receptor γ2 alternative splicing via a distal downstream UCAU-rich intronic splicing enhancer. Mol Cell Biol 23 4687 4700
67. JelenNUleJZivinMDarnellRB 2007 Evolution of Nova-dependent splicing regulation in the brain. PLoS Genet 3 e173 doi:10.1371/journal.pgen.0030173
68. BradfieldJPQuHQWangKZhangHSleimanPM 2011 A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci. PLoS Genet 7 e1002293 doi:10.1371/journal.pgen.1002293
69. FlorezJC 2008 Newly identified loci highlight beta cell dysfunction as a key cause of type 2 diabetes: Where are the insulin resistance genes? Diabetologia 51 1100 1110
70. McCarthyMI 2010 Genomics, type 2 diabetes, and obesity. N Engl J Med 363 2339 2350
71. ChenMCProostPGysemansCMathieuCEizirikDL 2001 Monocyte chemoattractant protein-1β is expressed in pancreatic islets from prediabetic NOD mice and in interleukin-1-exposed human and rat islet cells. Diabetologia 44 325 332
72. MartinAPGrisottoMGCanasto-ChibuqueCKunkelSLBrombergJS 2008 Islet expression of M3 uncovers a key role for chemokines in the development and recruitment of diabetogenic cells in NOD mice. Diabetes 57 387 394
73. RoepBOKleijwegtFSvan HalterenAGBonatoVBoggiU 2010 Islet inflammation and CXCL10 in recent-onset type 1 diabetes. Clin Exp Immunol 159 338 343
74. UnoSImagawaASaishoKOkitaKIwahashiH 2010 Expression of chemokines, CXC chemokine ligand 10 (CXCL10) and CXCR3 in the inflamed islets of patients with recent-onset autoimmune type 1 diabetes. Endocr J 57 991 996
75. LiuDCardozoAKDarvilleMIEizirikDL 2002 Double-stranded RNA cooperates with interferon-γ and IL-1β to induce both chemokine expression and nuclear factor-κB-dependent apoptosis in pancreatic β-cells: potential mechanisms for viral-induced insulitis and β-cell death in type 1 diabetes mellitus. Endocrinology 143 1225 1234
76. Igoillo-EsteveMMarselliLCunhaDALadriereLOrtisF 2010 Palmitate induces a pro-inflammatory response in human pancreatic islets that mimics CCL2 expression by beta cells in type 2 diabetes. Diabetologia 53 1395 1405
77. SimmonsMJFanGZongWXDegenhardtKWhiteE 2008 Bfl-1/A1 functions, similar to Mcl-1, as a selective tBid and Bak antagonist. Oncogene 27 1421 1428
78. KarsanAYeeEHarlanJM 1996 Endothelial cell death induced by tumor necrosis factor-α is inhibited by the Bcl-2 family member, A1. J Biol Chem 271 27201 27204
79. SarkarSAKutluBVelmuruganKKizaka-KondohSLeeCE 2009 Cytokine-mediated induction of anti-apoptotic genes that are linked to nuclear factor κB (NF-κB) signalling in human islets and in a mouse beta cell line. Diabetologia 52 1092 1101
80. ChenLWillisSNWeiASmithBJFletcherJI 2005 Differential targeting of prosurvival Bcl-2 proteins by their BH3-only ligands allows complementary apoptotic function. Mol Cell 17 393 403
81. HermanMDNymanTWelinMLehtioLFlodinS 2008 Completing the family portrait of the anti-apoptotic Bcl-2 proteins: crystal structure of human Bfl-1 in complex with Bim. FEBS Lett 582 3590 3594
82. ChenMManleyJL 2009 Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat Rev Mol Cell Biol 10 741 754
83. HartmannBValcarcelJ 2009 Decrypting the genome's alternative messages. Curr Opin Cell Biol 21 377 386
84. UleJJensenKBRuggiuMMeleAUleA 2003 CLIP identifies Nova-regulated RNA networks in the brain. Science 302 1212 1215
85. SuckowATComolettiDWaldropMAMosedaleMEgodageS 2008 Expression of neurexin, neuroligin, and their cytoplasmic binding partners in the pancreatic beta-cells and the involvement of neuroligin in insulin secretion. Endocrinology 149 6006 6017
86. AtoufFCzernichowPScharfmannR 1997 Expression of neuronal traits in pancreatic β cells. Implication of neuron-restrictive silencing factor/repressor element silencing transcription factor, a neuron-restrictive silencer. J Biol Chem 272 1929 1934
87. CardozoAKBerthouLKruhofferMOrntoftTNicollsMR 2003 Gene microarray study corroborates proteomic findings in rodent islet cells. J Proteome Res 2 553 555
88. MartensGAJiangLHellemansKHStangeGHeimbergH 2011 Clusters of conserved β cell marker genes for assessment of β cell phenotype. PLoS ONE 6 e24134 doi:10.1371/journal.pone.0024134
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2012 Číslo 3
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- PIF4–Mediated Activation of Expression Integrates Temperature into the Auxin Pathway in Regulating Hypocotyl Growth
- Metabolic Profiling of a Mapping Population Exposes New Insights in the Regulation of Seed Metabolism and Seed, Fruit, and Plant Relations
- A Splice Site Variant in the Bovine Gene Compromises Growth and Regulation of the Inflammatory Response
- Comprehensive Research Synopsis and Systematic Meta-Analyses in Parkinson's Disease Genetics: The PDGene Database