A functional genomic approach to identify reference genes for human pancreatic beta cell real-time quantitative RT-PCR analysis

Exposure of human pancreatic beta cells to pro-inflammatory cytokines or metabolic stressors is used to model events related to type 1 and type 2 diabetes, respectively. Quantitative real-time PCR is commonly used to quantify changes in gene expression. The selection of the most adequate reference gene(s) for gene expression normalization is an important pre-requisite to obtain accurate and reliable results. There are no universally applicable reference genes, and the human beta cell expression of commonly used reference genes can be altered by different stressors. Here we aimed to identify the most stably expressed genes in human beta cells to normalize quantitative real-time PCR gene expression. We used comprehensive RNA-sequencing data from the human pancreatic beta cell line EndoC-BH1, human islets exposed to cytokines or the free fatty acid palmitate in order to identify the most stably expressed genes. Genes were filtered based on their level of significance (adjusted P-value >0.05), fold-change (|fold-change| <1.5) and a coefficient of variation <10%. Candidate reference genes were validated by quantitative real-time PCR in independent samples. We identified a total of 264 genes stably expressed in EndoC-BH1 cells and human islets following cytokine- or palmitate-induced stress, displaying a low coefficient of variation. Validation by quantitative real-time PCR of the top five genes ARF1, CWC15, RAB7A, SIAH1 and VAPA corroborated their expression stability under most of the tested conditions. Further validation in independent samples indicated that the geometric mean of ACTB and VAPA expression can be used as a reliable normalizing factor in human beta cells.