Analysis artefacts of the INS-IGF2 fusion transcript
- Rasmus Wernersson1, 2Email author,
- Thomas Frogne†3,
- Claude Rescan†3,
- Lena Hansson†1,
- Christine Bruun3,
- Mads Grønborg3,
- Jan Nygaard Jensen3Email author and
- Ole Dragsbæk Madsen3Email author
© Wernersson et al. 2015
Received: 1 December 2014
Accepted: 13 July 2015
Published: 29 July 2015
In gene expression analysis, overlapping genes, splice variants, and fusion transcripts are potential sources of data analysis artefacts, depending on how the observed intensity is assigned to one, or more genes. We here exemplify this by an in-depth analysis of the INS-IGF2 fusion transcript, which has recently been reported to be among the highest expressed transcripts in human pancreatic beta cells and its protein indicated as a novel autoantigen in Type 1 Diabetes.
Through RNA sequencing and variant specific qPCR analyses we demonstrate that the true abundance of INS-IGF2 is >20,000 fold lower than INS in human beta cells, and we suggest an explanation to the nature of the artefacts which have previously led to overestimation of the gene expression level in selected studies. We reinvestigated the previous reported findings of detection of INS-IGF2 using antibodies both in Western blotting and immunohistochemistry. We found that the one available commercial antibody (BO1P) raised against recombinant INS-IGF2 show strong cross-reaction to native proinsulin, and we did not detect INS-IGF2 protein in the human beta cell line EndoC-βH1. Furthermore, using highly sensitive proteomics analysis we could not demonstrate INS-IGF2 protein in samples of human islets nor in EndoC-βH1.
Sequence features, such as fusion transcripts spanning multiple genes can lead to unexpected results in gene expression analysis, and care must be taken in generating and interpreting the results. For the specific case of INS-IGF2 we conclude that the abundance of the fusion transcript/protein is exceedingly lower than previously reported, and that current immuno-reagents available for detecting INS-IGF2 protein have a strong cross-reaction to native human proinsulin. Finally, we were unable to detect INS-IGF2 protein by proteomics analysis.
KeywordsInsulin INS-IGF2 Beta cell Fusion-transcript Gene expression Proteomics Antibody Bioinformatics Analysis artefact Pancreatic islets Diabetes
Gene expression analysis artefacts
As part of our research on pancreatic beta cell biology, our group has been investigating the transcripts originating from the insulin loci using both microarrays and RNA-sequencing platforms. On several occasions we have noticed issues with the way overlapping transcripts are handled in gene expression analysis pipelines, and this sparked the interest for looking further into the underlying reasons for the problems. Our initial observations were as follows:
For RNA sequencing analysis a related problem is encountered. The individual reads will typically be mapped to a set of transcripts from the organism in question, and will be assigned multiple times to transcripts sharing sequences. In most cases this is the desired behaviour, and will be a way to handle splicing variants of a given gene (and a more thorough analysis of exon/exon junctions can be investigated). However, in the case of transcripts spanning several genes, there is a risk that the signal from one gene bleeds into the overlapping transcript, and this may not be apparent when the results are reported.
With this in mind, we decided to re-investigate the observation by Nica et al.  that the INS-IGF2 fusion transcript was found to be the second most abundant mRNA in human beta cells, an observation that has generated much interest in the beta cell biology field since INS-IGF2 has previously been annotated to be translated  and moreover, recently proposed as a novel autoantigen in Type 1 Diabetes (T1D) .
Case background: beta cell biology
The pancreatic beta cell is a key component in regulating blood glucose homeostasis as it represents our sole source of insulin, which is required for peripheral tissues to internalize and utilize glucose as an energy source. Lack of adequate functional beta cell mass leading to relative insulin insufficiency is the common denominator of all forms of diabetes. T1D is characterized by an absolute loss of beta cells as a result of autoimmune destruction. Type 2 diabetes (T2D) is where a relative beta cell insufficiency results from failure of compensatory upregulation of insulin production capacity in response to increasing peripheral insulin resistance. Therefore much effort in, both preventive and curative, diabetes research has aimed at understanding the molecular networks, both at the level of transcriptional and translational regulation of genes controlling development and at the function and survival of the pancreatic beta cell. Recent advancements in the ability to purify single cell types from the Islet of Langerhans combined with omics’ approaches have led to detailed information of the transcriptomic environment of beta cells.
Gene-knockouts, and animal models, in particular targeting transcription factors combined with genome-wide association studies, have in the past decade provided us with a comprehensive map of genes and loci that may cause or contribute to diabetes, and it is evident that better understanding of molecular control(s) of functional beta cell mass will be key to design strategies to defeat diabetes.
All classical studies, not surprisingly, report a very high abundance of insulin mRNA, and insulin represents one of the most abundant cell specific mRNA’s known in our body. It was therefore surprising to realize that the second most abundant beta cell mRNA, recently reported in FACS purified human beta cells , was encoding INS-IGF2. Values comparable to glucagon mRNA levels in purified human alpha cells was reported for INS-IGF2 suggesting that a putative and highly abundant beta cell protein could have been overlooked. Moreover, recent studies suggested that INS-IGF2 could be yet another beta cell specific antigen recognized by autoantibodies in T1D .
In our reanalysis of the RNA sequencing data from Nica et al.  we did not map the reads to transcripts but rather investigated an unbiased mapping of the reads to the genomic region on chromosome 11 containing the INS, IGF2 and INS-IGF2 transcripts. By visual inspection of the distribution of the reads (Fig. 1), it is clearly observed that the regions specific to INS-IGF2 are lowly expressed, whereas the INS specific exons are (as expected in beta cells) highly expressed.
Both coding exons of INS being expressed at the same order of magnitude is a clear indication that (a) most reads mapped to the exon that is shared by both INS and INS-IGF2 belong to the INS-001 transcript, and (b) the perceived high expression of the INS-IGF2 transcript is due to the read count being assigned to both genes, where the difference in level is due to the different lengths of the two genes.
The data analysis artefact is especially pronounced in this case, since INS in itself constitutes a large part of the transcript pool in beta cells. In our own experiments using the human beta cell line EndoC-βH1 we observe the same trends, both with RNA sequencing and micro-array analysis (unpublished, data not shown).
qPCR analysis of INS-IGF2, INS and IGF2
cDNA input (ng)
GAPDH normalized exp. Rel. To INS
To finalise our conclusion that INS-IGF2 expression levels have been vastly overestimated, we investigated the abundance of the INS-IGF2 fusion protein using a strategy involving both immunochemical detection and mass spectrometry.
The control-antibody GS-9A8 (Madsen, Frank et al. ), a well characterized proinsulin-specific antibody binding a linear epitope spanning the B-C junction (Fig. 2) show full cross-reaction to recombinant INS-IGF2 (Fig. 3, panel a upper)—as would be expected with the intact proinsulin-B-C junction preserved in the INS-IGF2 protein (Fig. 2).
From band intensities in Western blotting experiments we conclude that the commercial polyclonal antibody BOA1 has higher affinity towards INS-IGF2 than for proinsulin and yet it fails to detect any INS-IGF2-like immunoreactivity in native human EndoC-βH1. In transduced EndoC-βH1 the expression of INS-IGF2 is much lower than proinsulin. Yet a weaker band is readily detected by BOA1 but completely absent in native human EndoC-βH1. Similarly, the GS-9A8 antibody having a presumed equal high affinity for both proinsulin and INS-IGF2 completely fails to detect an INS-IGF2 band in native human EndoC-βH1. In summary, these data are in line with the low levels of transcript we see in the human beta cell line, EndoC-βH1, which again raise the question do the protein exist?
We do observe that the antibody raised against INS-IGF2 (BO1P) indeed stain human beta cells on section of human pancreas (Fig. 3, panel b), but we cannot deduce if this is only due to the expected cross-reaction to proinsulin. In fact, the two antibodies, GS-9A8 and BO1P display very similar staining of proinsulin-like immuno-reactivity in human beta cells (Fig. 3, Panel b).
In order to further investigate the abundance of the INS-IGF2 protein level we conducted a comprehensive analysis of both EndoC-βH1 cells as well as isolated human islets using high-resolution mass spectrometry. Human islets were digested by trypsin and fractionated by either hydrophilic interaction chromatography (HILIC) or strong cation exchange (SCX) chromatography prior to reverse phase liquid chromatography tandem mass spectrometry (LC–MS). Our 2D fractionation method combined with high resolution MS allowed us to identify a total of 8,441 distinct human islet proteins. As expected, we were able to identify human insulin with five unique peptides and a high Mascot protein score (>500).
Moreover, we did identify proteins known to be expressed at a notoriously low level including transcription factors (e.g. FOS-B, NAC-2, SOX4 and SOX5), nuclear receptors (e.g. NR2F6 and NR5A2) and nuclear receptor coactivators (e.g. NCOA2, NCOA4 and NCOA6). However, even though we identified close to 8,500 proteins in our analysis we were not able to identify peptides specific to the INS-IGF2 protein in the samples. This indicates that INS-IGF2 is at most expressed in exceedingly low abundance in human islets.
In the proteomic analysis of the human beta cell line EndoC-βH1 we again identified a large number of distinct proteins (close to 7,000) including insulin and low abundant proteins but again, we were not able to identify peptides specific to the INS-IGF2 protein.
However, the caveat one has to keep in mind with MS-based proteomics experiments is the fact that absence of evidence is not evidence of absence. Therefore, in order to ensure that the INS-IFG2 protein could indeed be identified by our proteomics approach, HEK cells transduced with either INS-IGF2 or a control construct were subjected to the same analysis. Both samples were analysed by LC–MS and peptides specific for INS-IGF2 were readily observed only in the sample transduced with the INS-IGF2 construct. The missing evidence of INS-IGF2 protein in human islets and in EndoC-βH1 cells is therefore not related to our LC–MS approach.
In summary, based on sensitive proteomic and immunologic detection methods using adequate positive controls, we conclude that INS-IGF2 protein is expressed in exceedingly low abundance in normal human islet beta cells as well as in EndoC-βH1. This finding correlates well with the exceedingly low expression level of the INS-IGF2 mRNA.
To increase the awareness of the type of data interpretation artefacts fusion transcripts (or other overlapping sequence features) can inflict upon gene expression analysis for both RNA sequencing and micro-array experiments. In our opinion, the main problem is that the analysis artefacts can go unnoticed, which can lead to the publication of misleading results and wasted efforts in follow-up experimental work. In addition to the issue with the abundance of INS-IGF2 being over-reported, we have seen several published micro-array studies, where INS wrongly was reported to be missing on the array (with comments along the line of “…Ins gene was excluded from ranking as it was not annotated in our Affymetrix chip…” ).
The problem with over-reporting of overlapping transcripts may actually be more prevalent than immediately evident, as most cases are likely to be much less extreme than the case of insulin and INS-IGF2 in beta cells and thus go unnoticed. In fact, there are in total 459 human protein-coding transcripts that share exons in Ensembl 77, that could potentially be affected by the same misassignment as reported here.
On a more fundamental level, we recognize that biology is complicated and it is difficult to devise one true interpretation of the transcripts. Furthermore, a complicating factor is that as more and more transcripts variants are sequenced and deposited into the genomic databases (e.g. from large scale efforts such as the ENCODE project ), the risk of encountering this kind of problems is therefore likely to increase.
Part of this problem is, in essence, that these variants will carry the same weight in the analysis, even if some of them are exceedingly rare. On a more technical note, there is also the issue of what is classified as a gene in its own right as opposed to a transcript variant. In the case of INS-IGF2 part of the problem is that it’s listed as an independent gene and not as a splice variant of INS, as the analysis artefact could otherwise be side-stepped by aggregating the signal at the gene level.
While devising a “perfect” solution to the problem with overlapping genes is far from trivial, following an approach of listing which other transcripts may contribute to the signal of any given gene, could help to bring potential problems out in the open.
Lastly, it should be mentioned that several recent publications of islet transcriptomics data do not report on a faulty INS-IGF2 expression as they specifically handle this issue. Eizirik et al.  used the Flux Capacitor approach  to specifically handle the problem of reads mapping to exonic regions shared by multiple transcripts. Moran et al.  addressed the multiple transcript mapping issue by applying custom parameters in the RNA sequencing pipeline. The issue of handling exons from multiple transcripts is highly dependant on the software tool used, including its version and default parameters. For example, for users of the TopHat2 software package it should be noted, that the default settings for handling multiple transcript mapping was changed around September 2012 .
In the more specific case of INS-IGF2 it is important to avoid a chase for an elusive transcript for the many groups of scientists working with beta cells and islet biology. However, the big question remains on how important this transcript is? And does it indeed lead to the production of a novel beta cell specific protein? We show that the INS-IGF2 antibody as well as our proinsulin-specific antibody recognize an expected band of MW 22,000 in INS-IGF2 transduced cells (Hek293 and EndoC-βH1), as well as the native proinsulin band from the human beta cell line. The work of Kanatsuna et al.  concludes that INS-IGF2 represents a novel autoantigen in T1D, as they use in vitro translated/labeled INS-IGF2 to detect specific autoantibodies in T1D patients. However, “pre-proinsulin” autoantibodies would possibly also cross-react effectively in this assay, and such autoantibodies are well known in T1D. As an example, we show that proinsulin-B/C junction specific antibody fully cross-reacts to INS-IGF2. More specific reagents, such as antibodies raised against the unique c-terminal fragment of the predicted INS-IGF2 protein, are required in order to detect the protein specifically by immunocytochemical techniques.
Finally, our data failed to demonstrate the existence of the INS-IGF2 mRNA and protein in EndoC-βH1 cells. This was further supported by lack of detection of the protein by proteomics analysis, both in human islets and the beta cell line.
Overlapping sequence features are problematic to handle in expression analysis pipelines. This challenge potentially leads to situations were some variants appears to be much higher expressed than data can support, or in the extreme case to completely mask out other genes in the vicinity. A large part of the problem is the “black box” nature of the gene expression analysis pipeline in the eyes of researchers not experts in bioinformatics, and that such errors may therefore go unnoticed. While the underlying mapping issue is far from trivial to solve, we suggest that by bringing the problem out in the open (e.g. by listing other factors contributing to the signal of a given gene), it will be easier to flag potential problems with a given data set, and to help avoid over-, or underestimation of gene expression levels.
Our in-depth analysis concludes that INS-IGF2 abundance is >20,000 times lower or undetectable at the level of mRNA and protein, respectively, compared to previous reported findings in human beta cells.
Finally, we suggest that pre-proinsulin antibodies (T1D autoantibodies) will also cross-react to recombinant INS-IGF2 protein, as exemplified by a proinsulin-B-C junction specific monoclonal antibody, GS9A8.
In conclusion, we suggest that more specific antibody reagents (towards against INS-IGF2 unique epitopes) need to be developed to facilitate and support investigations into whether INS-IGF2 fusion protein actually might be induced to elicit an autoantigen response in T1D.
Reanalysis of RNA sequencing data
For the re-analysis of the human Islet of Langerhans RNA sequencing data from Nica et al. , data was downloaded from the European Genome-phenome Archive  (accession: EGAS00001000442) as processed BAM files .
The BAM files follow the original analysis performed by Nica et al. : briefly, the transcript-to-genome mapping was done using assembly version GRCh37/hg19 from February 2009 and GENCODE annotation v10 , which corresponds to Ensembl  v65. They further used BWA, Burrows-Wheeler Aligner , and filtered the alignments for correct orientation of the mapped mate pairs with a maximum insert size of 500 kb, and a minimum mapping quality score of 10.
In our re-analysis the processed BAM files were sorted and indexed and the relevant genomics region, chr11:2103768-2222439 (CGRh37/hg19 coordinates), was retrieved as Sequence Alignment/Map (SAM)-files via samtools . Finally the read depth were calculated via samtools (v 0.1.19) and plotted using R (v. 3.0.2) —we here observed that the read depth had been capped at 8,000. Exon/Intron coordinates for INS, IGF2 and INS-IGF2 transcripts were obtained from ENSEMBL  version 75 (also CGRh37/hg19 coordinates). See Additional file 1: Figure S1 for an in-depth overview of transcripts.
DNA sequences from the following transcripts were downloaded from ENSEMBL v. 75 : INS-IGF2-001 (Ensembl: ENST00000397270), INS-001 (Ensembl: ENST00000381330), IGF2-001 (Ensembl: ENST00000416167). For primers uniquely targeting INS and INS-IGF2 coding exons the primer placement strategy shown in Additional file 1: Figure S1(panel B) was used, and following the verification of the unique parts of the transcripts using global pairwise alignment the primer sets listed below were designed. For the IGF2 coding exons there are no overlaps with other genes (see Fig. 1 and Additional file 1: Figure S1 panel A), and the probes were trivial to design.
Primerset 1 (INS specific), Forward: GCAGGTGGAGCTGGGC, Reverse: GGAGGAGAGGGACAAAGCTG. Primerset 2 (INS-IGF2 specific), Forward: CTACCTAGTGTGCGGGGAA, Reverse: ATTGTTCCACAATGCCACGC. Primerset 3 (IGF2 specific), Forward: GATGCTGGTGCTTCTCACCT, Reverse: CAGACGAACTGGAGGGTGTC. For normalization of cDNA input we used GAPDH, Forward: AGGGCTGCTTTTAACTCTGGT, Reverse: CCCCACTTGATTTTGGAGGGA. Primers were synthesized by Eurofins.
Quantitative RT PCR was carried out by standard methods. Briefly, RNA was purified with RNeasy (Qiagen) and checked on a Nanodrop instrument. All OD260/280 was ~2. The cDNA was made from 0.5 µg RNA in 20 µl iSCRIPT (Roche) reaction. For qPCR 2.5 µl of 10× (EndoC-βH1 cells) or 1,000× (Human islets) diluted cDNA was analysed with specific primers (see Table 1) and Brilliant II SYBR green mastermix (Agilent) using a Stratagene 3000p instrument. All PCR products were verified by melting curve analysis.
The following cell lines were used in the study. EndoC-βH1 Human pancreatic beta cell line (EndoC-βH1) . The cells were grown as described in ref . HEK293 Human embryonic kidney cell line (ATCC CRL-1573). The HEK293 cells were grown in DMEM with 10% FBS and 1% P/S.
Human pancreatic islets
Human islets were obtained via Prodo Laboratories Inc., CA, USA. Human pancreases were procured from cadaveric donors after written informed research consent was provided by donor relatives. All experiments were performed in agreement with national legislation and institutional ethical rules. Upon arrival, the islets were cultured in CMRL1066 with 5.6 mM glucose supplemented with 10% FCS, 2 mM l-Glutamine, 100 units/ml Penicillin and 100 μg/ml streptomycin at 37°C, under 95% air/5% CO2.
Human islet donor data: 3 donors, 1 female, 2 males; age, mean 34 years (range 30–38); BMI, mean 31.5 (range 28.5–33.1); purity, mean 90% (range 85–95); viability, mean 94.7 (range 94–95%). None of the donors were diagnosed with diabetes.
The INS-IGF2 cDNA sequence was obtained from GenBank  (accession id: NM_001042376.2) and the corresponding DNA was synthesized and inserted into a lentivirus vector (Gentarget Inc.). Lentivirus was produced in HEK293 cells and purified into PBS by ultracentrifugation. The MOI was 10 for EndoC-βH1 cells and 2 for HEK293 cells. Selection was done with puromycin for 1–3 weeks. As control we used the same backbone vector with a GFP cDNA insert.
Human pancreas paraffin sections (Zyagen) were dewaxed and rehydrated to double-distilled H2O. Sections were microwave treated for 15 min in Tris-EGTA buffer (pH 9.0) and were allowed to cool for 30 min and rinsed in double-distilled H2O. Slides were blocked with TNB blocking buffer (Perkin Elmer) for 30 min, followed by incubation with the primary antibody guinea pig anti Glucagon (1:500; Linco) in combination with the mouse monoclonal anti proinsulin clone GS9A8 (1:100; Developmental Studies Hybridoma Bank), or the mouse polyclonal anti INS-IGF2 (1:100; BO1P Abnova) diluted in TNB o.n at room temperature. Three washes in PBS were then followed by incubation with fluorescent labelled secondary antibodies (Cy2 or Cy3 conjugated donkey anti guinea pig and anti rabbit IgG (1:400; Jackson ImmunoResearch)) and DAPI diluted in TNB and three washes in PBS before mounting. Images were acquired using confocal microscopy (Zeiss LSM 510).
Cultured HEK293 or EndoC-βH1 cells were harvested in Cell Extraction Buffer (Life Technologies) freshly supplemented with Complete Protease Inhibitor Cocktail (Roche). Samples were run on 12% NuPAGE SDS gels, under reducing conditions in MES-SDS buffer (Invitrogen) and then transferred on PVDF membrane for chemiluminescent Western blot. Anti-proinsulin and INS-IGF2 antibodies have been used at 1:1,000 dilution, monoclonal anti Beta actin-peroxidase (Sigma) at 1:25,000. Secondary antibody goat anti mouse IgG (Santa Cruz Biotechnology) was diluted 1:2,000.
Human islets or EndoC-βH1 cells were lysed in 0.5% RapiGest (Waters, Milford, MA, USA) dissolved in 50 mM ammonium bicarbonate. The islets were subsequently reduced in DTT (10 mM) and alkylated with 45 mM chloroacetamide. The sample was digested with trypsin (1:50 w/w) O.N. at 37 degrees. To precipitate the RapiGest (Waters, Milford, MA, USA) 10% TFA was added to a final concentration of 0.5%. The peptides were centrifuged for 20 min at room temperature and purified in a SepPak (50 mg) cartridge (Waters, Milford, MA, USA). The peptides were eluted with 50% acetonitrile and lyophilized. For HILIC fractionation, the peptides were dissolved in 4 ul buffer A (0.1% TFA) and 36 ul buffer B (98% acetonitrile, 0.1% TFA). The peptides were fractionated using a gradient of 55 min (95–60% B). The fractionation was done using an Agilent uHPLC 1290 system (Agilent Technologies, Santa Clara, USA) For SCX fractionation, the column (StageTips with 3 M Empore SCX disks, 3 M Empore, MN, USA) was conditioned with (1) 100% acetonitrile, (2) 50% acetonitrile pH4.0, (3) 50% acetonitrile pH 11, (4) 40% acetonitrile with 0.1% TFA. The peptides were loaded in 40% acetonitrile with 0.1% TFA and subsequently step eluted using a buffer containing 20 mM acetic acid, 20 mM boric acid, 20 mM phosphoric acid, 50% acetonitrile with pH adjusted to 4.0; 4.5; 5.0; 6.5; 7.5; 11.0, respectively. Peptides from HILIC and SCX fractionation were analysed by LC–MS using an Easy-nLC 1000 system coupled to a Q-Exactive Orbitrap mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). The peptides were separated on a 90 min gradient (8–40% acetonitrile in 0.1% formic acid) on a reverse-phase column (1.9 µm beads, Dr Maisch GmbH, Ammerbach, Germany). The raw data was searched using the Mascot search engine in Proteome Discoverer v1.4 (Thermo Scientific, Bremen, Germany) with a peptide false discovery rate of less than 1%, which was estimated on the number of reverse hits. Carbamidomethylation and oxidation of methionine were specified as variable modifications.
RNA sequencing analysis and initial discovery of the expression overestimation by RW and LH. TF did transductions, protein & RNA preparations and qPCR. CR did western blot and immunohistochemistry. CB cultured and validated the human islets. Mass spectrometry analysis by MG. Project lead by JNJ and ODM. Paper primarily drafted by RW, JNJ, LH and ODM with input from all authors. The final version of the manuscript was reviewed by all authors. All authors read and approved the final manuscript.
The authors would like to thank Lin Chen for technical assistance and Hanni Willenbrock Thomsen for in-house RNA sequencing analysis.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Nica AC, Ongen H, Irminger J-C, Bosco D, Berney T, Antonarakis SE et al (2013) Cell-type, allelic, and genetic signatures in the human pancreatic beta cell transcriptome. Genome Res 23:1554–1562PubMed CentralPubMedView ArticleGoogle Scholar
- Monk D, Sanches R, Arnaud P, Apostolidou S, Hills FA, Abu-Amero S et al (2006) Imprinting of IGF2 P0 transcript and novel alternatively spliced INS-IGF2 isoforms show differences between mouse and human. Hum Mol Genet 15:1259–1269PubMedView ArticleGoogle Scholar
- Kanatsuna N, Taneera J, Vaziri-Sani F, Wierup N, Larsson HE, Delli A et al (2013) Autoimmunity against INS-IGF2 protein expressed in human pancreatic islets. J Biol Chem 288:29013–29023PubMed CentralPubMedView ArticleGoogle Scholar
- Madsen OD, Frank BH, Steiner DF (1984) Human proinsulin-specific antigenic determinants identified by monoclonal antibodies. Diabetes 33:1012–1016PubMedView ArticleGoogle Scholar
- Taneera J, Lang S, Sharma A, Fadista J, Zhou Y, Ahlqvist E et al (2012) A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets. Cell Metab 16:122–134PubMedView ArticleGoogle Scholar
- Kellis M, Wold B, Snyder MP, Bernstein BE, Kundaje A, Marinov GK et al (2014) Defining functional DNA elements in the human genome. Proc Natl Acad Sci USA 111:6131–6138PubMed CentralPubMedView ArticleGoogle Scholar
- Eizirik DL, Sammeth M, Bouckenooghe T, Bottu G, Sisino G, Igoillo-Esteve M et al (2012) The human pancreatic islet transcriptome: expression of candidate genes for type 1 diabetes and the impact of pro-inflammatory cytokines. PLoS Genet 8:e1002552PubMed CentralPubMedView ArticleGoogle Scholar
- Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J et al (2010) Transcriptome genetics using second generation sequencing in a Caucasian population. Nature 464:773–777PubMedView ArticleGoogle Scholar
- Morán I, Akerman I, van de Bunt M, Xie R, Benazra M, Nammo T et al (2012) Human β cell transcriptome analysis uncovers lncRNAs that are tissue-specific, dynamically regulated, and abnormally expressed in type 2 diabetes. Cell Metab 16:435–448PubMed CentralPubMedView ArticleGoogle Scholar
- TopHat. http://ccb.jhu.edu/software/tophat/index.shtml. Accessed 1 Dec 2014
- Brooksbank C, Bergman MT, Apweiler R, Birney E, Thornton J (2014) The European Bioinformatics Institute’s data resources 2014. Nucleic Acids Res 42:D18–D25PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N et al (2009) 1000 Genome Project Data Processing Subgroup: the sequence alignment/map format and SAMtools. Bioinforma Oxf Engl 25:2078–2079View ArticleGoogle Scholar
- GENCODE–GENCODE data releases. http://www.gencodegenes.org/releases/. Accessed 1 Dec 2014
- Flicek P, Amode MR, Barrell D, Beal K, Billis K, Brent S et al (2014) Ensembl 2014. Nucleic Acids Res 42:D749–D755PubMed CentralPubMedView ArticleGoogle Scholar
- Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinforma Oxf Engl 25:1754–1760View ArticleGoogle Scholar
- R Core Team (2013) R: a language and environment for statistical computing. R Foundation for Statistical, Computing, ViennaGoogle Scholar
- Ravassard P, Hazhouz Y, Pechberty S, Bricout-Neveu E, Armanet M, Czernichow P et al (2011) A genetically engineered human pancreatic β cell line exhibiting glucose-inducible insulin secretion. J Clin Invest 121:3589–3597PubMed CentralPubMedView ArticleGoogle Scholar
- Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ et al (2013) GenBank. Nucleic Acids Res 41:36–42View ArticleGoogle Scholar