ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Note

The importance of selecting the appropriate reference genes for quantitative real time PCR as illustrated using colon cancer cells and tissue

[version 1; peer review: 2 approved]
* Equal contributors
PUBLISHED 22 Jan 2016
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Quantitative real-time reverse-transcription polymerase chain reaction (RT-qPCR) remains the most sensitive technique for nucleic acid quantification. Its popularity is reflected in the remarkable number of publications reporting RT-qPCR data. Careful normalisation within RT-qPCR studies is imperative to ensure accurate quantification of mRNA levels. This is commonly achieved through the use of reference genes as an internal control to normalise the mRNA levels between different samples. The selection of appropriate reference genes can be a challenge as transcript levels vary with physiology, pathology and development, making the information within the transcriptome flexible and variable. In this study, we examined the variation in expression of a panel of nine candidate reference genes in HCT116 and HT29 2-dimensional and 3-dimensional cultures, as well as in normal and cancerous colon tissue. Using normfinder we identified the top three most stable genes for all conditions. Further to this we compared the change in expression of a selection of PKC coding genes when the data was normalised to one reference gene and three reference genes. Here we demonstrated that there is a variation in the fold changes obtained dependent on the number of reference genes used. As well as this, we highlight important considerations namely; assay efficiency tests, inhibition tests and RNA assessment which should also be implemented into all RT-qPCR studies. All this data combined demonstrates the need for careful experimental design in RT-qPCR studies to help eliminate false interpretation and reporting of results.

Keywords

Quantitative real-time PCR, Normalisation, Reference Genes, NormFinder, Colon Cancer

Introduction

Gene expression analysis is a critical and important tool in molecular diagnostics and medicine14. Quantification of RNA transcripts is carried out using one of four common methods; reverse transcription polymerase chain reaction (RT-PCR)5, RNase protection assays6, northern blotting and in situ hybridisation7, and less commonly now using cDNA arrays8. At present, the most popular and widely used method for gene expression is fluorescence based quantitative real time PCR (RT-qPCR)9. It is the most sensitive and flexible of the quantitative methods with a capacity to detect and measure minute amounts of nucleic acids10,11. There are two types of quantitative methods that can be applied within RT-qPCR; absolute quantification and relative quantification. Absolute quantification relates the PCR signal to a standard curve to determine the input copy number of the gene of interest. In contrast, relative quantification evaluates the change in expression of a target gene relative to a reference group, for example an untreated control12.

When employing RT-qPCR to compare mRNA levels between two different test conditions, it is imperative that reference genes are utilised carefully9,10. Normalisation of the data with these reference genes is essential for correcting results of different amounts of input RNA, uneven loading, reverse-transcription yield, efficiency of amplification and variation within experimental conditions9,13. The mRNA of reference genes should be stably expressed and their expression should not be affected by experimental condition or by any human disease14. Numerous studies have demonstrated that common reference genes, such as β-Actin and GAPDH, which are largely accepted as being stably expressed within cells, can in fact show large variations in expression1518. Despite the awareness that validation of the stability of reference genes is an essential component for accurate RT-qPCR analysis, this consideration is still largely disregarded1921.

Further to this, it is reported that over 90% of gene expression analysis published in high impact journals used only one reference gene2224. It has since been widely documented that normalisation of data with a single reference gene can lead to inaccurate interruption of results10,25,26. Taken together, this highlights the importance of selecting the optimal number and type of reference genes for any RT-qPCR study. Other essential considerations such as; analysis of assay efficiency, testing for inhibition with biological samples and reporting the quality and integrity of input RNA are all highlighted in the ‘MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments’10.

In this study, we sought to highlight the importance of carefully-designed RT-qPCR studies in order to avoid the reporting of inaccurate and misleading information. We test a panel of nine candidate reference genes and report their stability between 2-dimensional and 3-dimensional HCT116 and HT29 colon cancer cell lines, as well as between normal and cancerous tissue from colon cancer patients. We also demonstrate useful tests that should be implemented within RT-qPCR studies to ensure that studies comply with the MIQE guidelines.

Methods

Cell culture

HCT116 (ATCC® CCL-247™) and HT29 (ATCC® HTB-38™) cell lines were obtained from ATCC. These cell lines were cultured in complete Dulbecco's modified essential medium (DMEM) supplemented with 10% of foetal bovine serum, 1% of penicillin/streptomycin and 1% of L-glutamine. All cells were incubated at 37°C in a humidified 95% air/5% CO2 environment. Cellular suspensions were obtained by adding 0.5% trypsin to the cultures and incubating at 37°C at 5% CO2.

3-dimensional cell cultures

Individual wells of a 6-well plate were coated with MatrigelTM (BD Biosciences) and placed in an incubator at 37°C for 30 min. Cell lines were trypsinized and counted. 50,000 cells/ml were resuspended in DMEM supplemented with 2% MatrigelTM. Cells were placed in MatrigelTM coated wells for 30 min at 37°C, after which DMEM supplemented with 2% MatrigelTM was added to the cultures. Cells were maintained in culture for 6 days in an incubator at 37°C, 5% CO2 with fresh medium added every 2 days. On day 6, cultures were harvested using EDTA/PBS and either fixed with paraformaldehyde (PFA) for confocal analysis (Zeiss LSM 710) or used for RNA extraction.

Clinical samples

Following ethical approval from the University Hospital Limerick’s Ethics Committee (ethical approval number 73/11), tissue samples measuring approximately 0.5cm in diameter were collected from patients undergoing surgery in University Hospital Limerick. Normal tissue from the patients was also collected approximately 10 cm away from the cancer tissue. Specimens were immediately placed in Allprotect tissue reagent (Qiagen) and stored at -80°C.

RNA extraction and cDNA synthesis

2-dimensional and 3-dimensional cell cultures were trypsinised as described above and frozen tissue was immersed in liquid nitrogen and ground into powder. Lysis buffer was added to the cells and tissue and the samples transferred to tubes using a 21-gauge needle. Total RNA was extracted as per Qiagen RNeasy Mini Kit instructions. RNA was quantified using a Nanodrop Spectrophotometer (Thermo Scientific) and stored at -80 degrees. RNA purity was evaluated by the ratio of absorbance at 260/280 nm and RNA quality was evaluated through visualization of the 28S:18S ribosomal RNA ratio on a 1% agarose gel. Total RNA (1 μg) was synthesised into cDNA using Vilo cDNA synthesis kit (Invitrogen) and stored at -20 degrees.

Real-time PCR

Real-time PCR was conducted using the ABI 7900 HT instrument (Applied Biosystems) following supplier instructions. Taqman® Gene Expression Assay Kits (Applied Biosystems) were used to analyse the gene expression of PKC coding genes. Data was normalised to either one reference gene or three reference genes (see below).

Assay efficiency test

The efficiency of each assay was determined by means of a calibration curve with the logarithm of the initial template concentration plotted on the x axis and the Cq plotted on the y axis. The slope of the graph was obtained and the PCR efficiency was calculated using the equation: 10-1/slope-1.

Inhibition test

Real-time PCR was conducted on corn DNA using a corn gene assay with a known Cq value of 24–26. Samples of cDNA from 2D and 3D HCT116 and HT29 cultures and from patient tissue was added to the reaction to test for an inhibitory components that may be present in these biological samples.

Selection of reference genes

All nine reference genes (Table 1) were purchased as pre-designed Taqman® Gene Expression Assays. The Cq value of each reference gene was determined for all biological samples. Normfinder was used to determine the most stable reference genes between 2D and 3D cell cultures as well as between normal and cancer tissue. Differences in gene expression levels of the PKC coding genes was determined using Pair Wise Fixed Reallocation Randomisation Test© as per REST© software. Within the software data was normalised to either the top reference gene or the top three reference genes.

Table 1. Description of the nine candidate housekeeper genes used in the study.

The accession numbers for each gene are taken from the National Center for Biotechnology Information.

SymbolNameAccession NumberFunction
B2MBeta 2 MicroglobulinNM_004048Important cell surface structure
PMM1Phosphomannomutase 1NM_002676Synthesis of the GDP-mannose and
dolichol-phosphate-mannose
TBPTATA Box Binding ProteinNM_001172085Transcription factor
RPLPOLarge Ribosomal ProteinNM_053275Ribosomal Protein
GUSBBeta GlucuronidaseNM_000181Glycoprotein
PGK1Phosphoglycerate Kinase 1NM_000291Glycolytic enzyme
ACTBBeta ActinNM_001101Cytoskeleton Protein
PPIAPeptidylprolyl Isomerase A
(Cyclophilin A)
NM_021130Catalyses the cis-trans isomerization
of proline imidic peptide bonds

Results

HCT116HCT116HCT116HCT116HCT116HCT116
2D2D2D3D3D3D
Gene Name Cq ValueCq ValueCq ValueCq ValueCq ValueCq Value
B2M22.14023422.1482922.12171222.45304322.41338322.370667
PMM129.40187529.27669328.97587629.9371729.79297429.651934
TBP26.8657426.9227326.81259728.27432428.13855228.062511
HRPT125.18523625.14627825.17779726.21830226.3284526.191511
RPLPO19.65694219.512119.68759220.15293720.1549720.111223
GUSB25.55094125.29034225.57396326.34415426.32348626.43148
PGK123.2743223.2743222.97286824.04533824.09898824.020983
ACTB27.57358727.3792527.50420426.93467126.87247826.94638
Dataset 1.Cq Values for reference genes in HCT116 cell lines.
The three Cq values for each reference gene is displayed for the 2-dimensional and 3-dimensional HCT116 cell cultures.
PRKCAPRKCEPRKCHPRKCQPRKCDPRKCZPRKCIB2MPMM1RPLPO
Cq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq Value
2D25.53479630.02182625.59707531.21085425.40270825.95780227.50833322.14023429.40187519.656942
2D25.65810830.02231225.60774231.49434325.49303425.92390427.67855522.1482929.27669319.5121
2D25.60383629.93523425.51254331.40276325.43472925.87727227.75883922.12171228.97587619.687592
3D24.35798829.41300824.76040530.88133424.92872225.32808326.55817422.45304329.9371720.152937
3D24.50706329.65078424.61309830.6015124.6246425.26624726.58215922.41338329.79297420.15497
Dataset 2.Cq Values for PKC coding genes and reference genes in HCT116 cell lines.
The three Cq values for each PKC coding gene and the appropriate reference genes is displayed for the 2-dimensional and 3-dimensional HCT116 cell cultures.
HT29HT29HT29HT29HT29HT29
2D2D2D3D3D3D
Gene Name Cq ValueCq ValueCq ValueCq ValueCq ValueCq Value
B2M19.99526419.94171319.93212920.515220.5898420.54583
PMM128.31492428.35563728.46154828.18984428.3085328.233744
TBP28.27079428.46041128.6335828.22400728.18035928.147839
HRPT124.83818424.83040224.85071824.75309424.73955724.66608
RPLPO20.07800120.14605920.15740419.74885219.75024419.656834
GUSB32.1632.17555232.37540429.42348525.49056825.269575
PGK123.16693122.69402522.74411623.6407423.66552524.008387
ACTB27.91863428.49796128.5322328.00906628.46990228.509188
Dataset 3.Cq Values for reference genes in HT29 cell lines.
The three Cq values for each reference gene is displayed for the 2-dimensional and 3-dimensional HT29 cell cultures.
PRKCAPRKCEPRKCHPRKCQPRKCDPRKCZPRKCIPMM1HRPTIPP1A
Cq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq Value
2D24.91495729.658425.35773331.51306524.59426925.44904523.83582328.31492424.83818419.944347
2D24.90939729.64091925.75363531.36996824.60879525.36551323.68431728.35563724.83040219.939817
2D24.92087729.64402825.77827631.86542324.68634225.35390323.75362628.46154824.85071819.802155
3D24.55074728.83278525.57593331.53737424.41540725.28512823.69891428.18984424.75309419.77429
3D24.42678528.67113725.60593831.45032124.38903825.17842123.62211628.3085324.73955719.735796
Dataset 4.Cq Values for PKC coding genes and reference genes in HT29 cell lines.
The three Cq values for each PKC coding gene and the appropriate reference genes is displayed for the 2-dimensional and 3-dimensional HT29 cell cultures.
Normal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueNormal TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer TissueCancer Tissue
Gene Name Cq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq Value
B2M21.9811321.85053421.9943218.03864917.94903617.99037714.78949514.82860214.79241217.3292117.27062417.32316817.38188617.44112217.46734423.72337723.69744523.65869118.58782818.59538518.7098916.55572116.73210316.58836718.19460318.2984818.21300918.57953818.69130118.687683
PMM128.5305628.78611828.69026626.07149126.19139526.26872624.51562324.54615624.63321726.23112926.4382426.2020225.67713425.93589825.80242733.5853732.02775232.4075426.65495126.91845726.79593825.84567625.7682125.64930226.39778726.51436426.4529226.56239326.5393326.572369
TBP33.22138633.38522733.0471529.04721628.97449728.82614527.27225527.27140627.31162631.89127332.3082832.37838431.38556930.94126331.0381232.9373332.2694131.95796828.78349528.76446528.78997827.3289327.3316427.53794929.04220829.09940729.25805329.88672430.01182229.895754
HRPT129.5743729.28316329.18941925.20540225.12086125.05607224.16598324.15162524.09254526.26012626.16370426.16266425.66966825.6037425.7452529.46576129.23300429.3155224.50010724.47905324.41435623.95137423.97200424.03919823.91689324.03495423.96445824.36298224.30720124.345879
RPLPO26.82353226.97188826.78510721.59473821.43688821.25205419.9963520.24705120.0668522.01290122.1831122.0771821.2820921.2125721.09284226.39476226.64626126.67112220.48825520.34034320.29869319.11118519.33826819.0291719.58546319.41458319.3966620.48320.42498820.36848
GUSB29.28672628.98339528.70597624.79539124.74218224.7705723.58712823.52171723.38170424.84998124.91282724.91384524.23187824.312924.16479930.02648429.7829229.93482824.2144424.08232123.97618522.4580922.49598522.5758322.65715222.60894622.4592924.32337824.35683324.128782
PGK126.73733926.6777426.5610523.00561522.95558222.96566420.96109820.88859221.68743522.79748322.79748322.89833822.3363622.276622.22342327.09194227.0527227.1946322.63693822.62550722.64176221.0975521.05491321.04413621.14839721.15241221.15013922.34392422.23067322.338123
ACTB29.7849729.44600929.52177824.96075825.03824624.80567721.59635221.51061822.02836422.7577322.68566722.7798922.12115922.58889421.97669233.66527631.91437534.123624.71927624.86236624.52747321.23394621.20209321.17825122.16221222.34902822.07873321.41009321.47676521.247252
Dataset 5.Cq Values for reference genes in normal and colon cancer tissue.
The three Cq values for each reference gene is displayed for the normal and colon cancer tissue.
PRKCAPRKCEPRKCHPRKCQPRKCDPRKCZPRKCIPGK1GUSBPP1A
Cq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq ValueCq Value
Normal Tissue 26.82627331.69177627.90734930.44792625.6178527.59838526.76390626.73733929.28672625.230883
Normal Tissue 27.03364231.49269927.91168630.50201425.78627627.59964226.91211526.6777428.98339525.126457
Normal Tissue 26.92827831.4889127.8960230.3921525.61997227.53412426.79927326.5610528.70597625.074203
Normal Tissue 31.60501735.6781333.95111535.1304930.41161533.09014530.91241623.00561524.79539120.968544
Normal Tissue 31.97254835.42518232.95326636.2203630.27338232.4345830.96775222.95558224.74218221.130848
Normal Tissue 31.85316735.6673633.55610734.57665330.04595432.86409830.82991222.96566424.7705721.146055
Normal Tissue 24.90762130.07382825.89400729.9265824.0594725.3273425.62509720.96109823.58712819.425198
Normal Tissue 25.10222629.86619425.88662529.9630523.96940225.50418925.61238520.88859223.52171719.42026
Normal Tissue 25.14757229.91539425.94743330.13907423.99979425.57900825.7332221.68743523.38170419.48943
Normal Tissue 25.1244228.16028426.39134828.69866223.86183225.89435825.01956222.79748324.84998121.279922
Normal Tissue 25.16895928.17232126.2754728.8834123.80961225.91702325.07737722.79748324.91282721.340395
Normal Tissue 24.94547728.15164826.34275228.77118523.8312425.77544624.96575422.89833824.91384521.253801
Normal Tissue 25.5252929.78129427.4886531.04027226.21488428.04591226.8676422.3363624.23187820.708403
Normal Tissue 25.58450529.91686427.5566630.542326.24748827.98285926.7994222.276624.312920.777203
Normal Tissue 25.5981529.72649227.5189430.48505626.186727.90770326.61946322.22342324.16479920.506758
Cancer Tissue 29.16719832.9432629.60906432.87066727.35580828.98259428.31248727.09194230.02648424.951777
Cancer Tissue 29.212432.83979829.62227232.36624527.446929.26442128.40253827.0527229.7829224.978256
Cancer Tissue 29.18106532.77680229.79466432.45764527.46087629.3086128.53318627.1946329.93482824.989716
Cancer Tissue 31.9304337.24955733.91776736.74299631.92115236.23898732.19390522.63693824.2144420.300858
Cancer Tissue 32.8954937.00724435.3647536.74299632.13334734.36691332.56575422.62550724.08232120.273832
Cancer Tissue 31.99256935.32929634.20422435.8620332.1998534.9576232.27038222.64176223.97618520.300154
Cancer Tissue 26.091730.75228926.90740430.53380225.36768526.2622226.96924621.0975522.4580919.597858
Cancer Tissue 26.32310130.50395626.89798230.31556125.22425326.2472626.8980521.05491322.49598519.55419
Cancer Tissue 24.38890829.34485226.36338429.21767424.9808226.12120425.9455921.04413622.5758319.532415
Cancer Tissue 24.49708429.55182326.25680729.3860925.26531825.95678325.771321.14839722.65715219.893276
Cancer Tissue 24.48128529.31637626.30442829.28660825.26972426.04345525.80748721.15241222.60894619.778872
Cancer Tissue 26.15909430.33476826.39735828.69730425.99310727.34946326.22256921.15013922.4592919.908897
Cancer Tissue 26.21276930.22098726.7500228.88694625.94893627.90910526.15552922.34392424.32337819.399368
Cancer Tissue 26.15151630.18375226.42757428.87178625.83584427.3647726.26631222.23067324.35683319.487078
Dataset 6.Cq Values for PKC coding genes and reference genes in normal and colon cancer tissue.
The three Cq values for each PKC coding gene and the appropriate reference genes is displayed for the normal and colon cancer tissue.
Dilution FactorCq ValueAverage Cq Value
1021.45344721.529772
21.615067
21.520802
124.03341523.90450733
23.832119
23.847988
0.127.25300427.26691967
27.365187
27.182568
0.0131.27787231.537185
32.356636
30.977047
0.00134.5810734.429821
Dataset 7.Cq values for sample assay efficiency test.
The Cq values obtained when HT29 cDNA was serial diluted and the gene PRKCA was amplified.
Sample TypeCq Value
No sample25.21
2D HCT11625.2
3D HCT11625.25
2D HT2925.21
3D HT2925.15
Normal Tissue Sample25.2
Dataset 8.Cq values for inhibition assay test.
The Cq values obtained for the corn assay when each of the sample types indicated are added to the RT-qPCR.

Comparison of reference genes in HCT116 2-dimensional and 3-dimensional cultures

In this study, we wanted to compare and validate the stability of reference genes used in quantitative real time PCR (RT-qPCR). To do this, HCT116 cells were grown in 2-dimensional and 3-dimensional cultures (Figure 1A). Following this, RNA was extracted from the cultures and cDNA was synthesised. Quantitative real time PCR was utilised to measure the variability in RNA transcript levels of 9 reference genes (RG) (Table 1) in the 2-dimensional and 3-dimensional cultures. The expression levels of the candidate reference genes were determined using the raw Cq values and NormFinder was then utilised to verify the stability of the genes. Normfinder ranks the RGs according to their stability values under the tested conditions. The top three stable genes when comparing 2-Dimensional and 3-dimensional HCT116 cultures were B2M, PMM1 and RPLPO, with B2M and PMM1 showing identical stability levels (Figure 1B, Dataset 1). Next, we wanted to elucidate the benefit of normalising data to more than one RG. To do this, we compared the expression of seven PKC coding genes in 3-Dimensional HCT116 cultures compared to 2-dimensional HCT116 cultures. The data was normalised to either one RG, B2M, or normalised to three RGs, B2M, PMM1 and RPLPO (Figure 1C, Dataset 2). Results indicate that using one RG gives fold changes that are greater than the fold changes obtained using three RGs.

6033185f-aa3c-4bb4-bce6-efe516249d61_figure1.gif

Figure 1. Reference genes in 2-dimensional and 3-dimensional HCT116 cultures.

The stability of the nine candidate reference genes between 2D and 3D HCT116 cultures was analysed using NormFinder. (A) Immunofluorescence images of HCT116 cells in 2D (100X) (left panel) and 3D (right panel) cell cultures (63X). (B) Table displaying the stability levels of the nine candidate reference genes between the 2D and 3D cultures. (C) Graph representing the fold change of PKC coding genes in 3D cultures compared to 2D cultures when using one reference gene (B2M) versus three reference genes (B2M, PMM1 and RPLPO).

Comparison of reference genes in HT29 2-dimensional and 3-dimensional cultures

Next, we compared the stability of the same 9 candidate reference genes in HT29 cultures. The cells were grown in 2-dimensional and 3-dimensional cultures (Figure 2A) before using RT-qPCR to determine the stability of the RGs between the two conditions. Normfinder revealed the most stable RGs were PMM1, HRPTI, PP1A and TBP (Figure 2B, Dataset 3) with PMM1 and HRPTI having a value of 0.001 and PP1A and TBP having a value of 0.002. Again, we examined the expression of the PKC coding genes in the cultures and normalised the data to one RG, PMM1, or three RGs, PMM1, HRPTI and PP1A (Figure 2C, Dataset 4). Our results indicate that there is variation in the fold changes obtained when using one RG versus three RGs. In some instances, genes that are found to be down-regulated when normalising with one RG are in fact up-regulated when normalising with three RGs.

6033185f-aa3c-4bb4-bce6-efe516249d61_figure2.gif

Figure 2. Reference genes in 2-dimensional and 3-dimensional HT29 cultures.

The stability of the nine candidate reference genes between 2D and 3D HT29 cultures was analysed using NormFinder. (A) Immunofluorescence images of HT29 cells in 2D (100X) (left panel) and 3D (right panel) cell cultures (63X). (B) Table displaying the stability levels of the nine candidate reference genes between the 2D and 3D cultures. (C) Graph representing the fold change of PKC coding genes in 3D cultures compared to 2D cultures when using one reference gene (PMM1) versus three reference genes (PMM1, HRPT1 and PP1A).

Comparison of reference genes in normal colon tissue versus colon cancer tissue

Following this, we wanted to examine the stability of the nine candidate RGs in normal and colon cancer tissue. We used fresh tissue samples that were excised from both the cancer tissue and normal distant tissue of individual patients (Figure 3A). As above, the expression levels of the nine candidate RGs were determined and Normfinder was used to establish the stability of the genes. PGK1, GUSB and PP1A were ranked as the most stable genes between normal and cancerous tissue (Figure 3B, Dataset 5). Next, we examined the change in PKC coding genes in colon cancer tissue when the data was normalised to one RG, PGK1, and normalised to three RGs, PGK1, GUSB and PP1A (Figure 3C, Dataset 6). The results demonstrate that using one RG can present fold changes that are up to 2-fold greater than when using three RGs.

6033185f-aa3c-4bb4-bce6-efe516249d61_figure3.gif

Figure 3. Reference genes in normal and colon cancer tissue.

The stability of the nine candidate reference genes between normal and cancer tissue was analysed using NormFinder. (A) Surgical image of specimen resected from a colon cancer patient. (B) Table displaying the stability levels of the nine candidate reference genes between the normal and cancer tissue. (C) Graph representing the fold change of PKC coding genes in cancer tissue compared to normal tissue (n=21) when using one reference gene (PGK1) versus three reference genes (PGK1, GUSB and PP1A).

Considerations when conducting RT-qPCR

Taken together, the results indicate that variations in fold changes can occur depending on the RG used to normalise data; making the selection of the correct RGs an imperative part of RT-qPCR studies. Further to this, the testing and reporting of assay efficiency is also essential to prevent the reporting of misinformation. Taking this into consideration, we examined the efficiency of all the RG assays and PKC coding genes assays (Figure 4A, Dataset 7). This information was inputted into the REST© software when establishing changes in gene expression between tested conditions. Another important consideration when designing RT-qPCR studies is the testing of your cDNA for any contaminants which could lead to the inhibition of the RT-qPCR reaction. For this reason, we added our samples to a standard RT-qPCR reaction using corn DNA and a gene that is known to have a Cq value of 24–26. If there were contaminants present in our cDNA samples this would inhibit the reaction resulting in a reduction in the Cq values10. However, we found no change in the Cq values for the reactions with the cDNA added, indicating the samples do not have any contaminates that will affect the amplification of our genes (Figure 4B, Dataset 8). It is also essential to report the quality assessment of the RNA templates, such as the RNA quantity, quality and integrity. We evaluated the RNA purity by the ratio of absorbance at 260/280 nm and RNA quality was assessed through visualization of the 28S:18S ribosomal RNA ratio on a 1% agarose gel (Figure 4C,D).

6033185f-aa3c-4bb4-bce6-efe516249d61_figure4.gif

Figure 4. Considerations to comply with during RT-qPCR.

(A) Representative graph of assay efficiency check. (B) Graph representing the inhibition test for all biological samples. (C) Representative graph from Nanodrop Spectrophotometer displaying the quantity and purity of the RNA. (D) Representative image of agarose gel displaying the 28S:18S ribosomal RNA ratio for RNA samples.

Discussion

The first publications using fluorescence-based quantitative real time PCR (RT-qPCR)2730 emerged almost a decade ago and since this time it has become the leading technique for gene expression analysis31,32. While RT-qPCR remains the most sensitive method for the detection of RNA transcripts33 there are also many challenges associated with the technique34,35. One of the major difficulties is the selection of appropriate reference genes for the normalisation of data. Hence the purpose of this study was to evaluate the stability in expression of nine candidate reference genes in two colon cancer cell lines as well as in normal and cancerous tissue from colon cancer patients. To help find the most suitable reference genes we selected genes which display a variation of functions within cells (Table 1).

Firstly, we examined the stability of the nine candidate reference genes between 2-dimensional and 3-dimensional HCT116 and HT29 cultures (Figure 1A, Figure 2A). The use of 3-dimensional cell cultures as cancer models is becoming increasingly popular3638; making the availability of appropriate reference genes important to help reduce the reporting of misinformation. When we examined the variation in expression between 2-dimensional and 3-dimensional HCT116 cells we found B2M, RPLPO and PMM1 to be the most stable genes between these two conditions (Figure 1B). Many publications have highlighted the problems associated with normalisation of data using only one reference gene22,23, for this reason we wanted to investigate differences in fold changes associated with normalising data to one reference gene compared to three reference genes. To do this, we investigated the change in expression in a selection of protein kinase c (PKC) coding genes between 2-dimensional and 3-dimensional HCT116 cultures. We examined PKCs as they are a group of proteins that are extensively studied for their role in oncogenic signalling39. Interestingly, when normalising the data to the reference gene B2M alone we found the change in expression of PKC coding genes was greater compared to normalisation with the reference genes, B2M, RPLPO and PMM1 together (Figure 1C). This finding highlights the need for normalisation with more than one reference gene to help eliminate the misinterpretation of fold changes in target genes.

Next we wanted to establish the stability of these reference genes in 2-dimensional and 3-dimensional HT29 cultures (Figure 2A). Normfinder ranked PMM1, HRPTI, PP1A and TBP as the most stable genes between these cultures (Figure 2B). It is important to note that despite the fact the treatments here were the same; there was a difference in the selected reference genes for HCT116 and HT29 cultures. This again emphasises the need to conduct stability tests on a panel of reference genes prior to all RT-qPCR studies to ensure data is normalised correctly. Again, we examined the difference in fold changes of PKC coding genes when normalising with varying numbers of reference genes. Importantly, we found that some target genes showing a down regulation when normalised with PMM1 showed no change when normalised to PMM1, HRPT1 and PP1A (Figure 2C).

RT-qPCR is the most common method used for the quantification of individual genetic differences in normal versus cancerous tissue9,34. Recent publications demonstrated that 97% of RT-qPCR studies contained on colorectal cancer contained information that was unreliable21. Thus, when examining difference in mRNA levels between normal and diseased tissue it is imperative the correct reference genes are used to normalise the data to prevent the presence of misleading information in the literature. Using normal and cancer tissue from CRC patients (Figure 3A) we examined the stability of the nine candidate reference genes, finding PGK1, GUSB and PP1A to be the most stably expressed (Figure 3B). As before, we compared the expression of PKC coding genes in normal and cancer tissue with the data normalised to either PGK1 alone or PGK1, GUSB and PP1A together. Strikingly we found that using only one reference gene results in a fold change that is up to 2 fold greater than when using three reference genes. This is a very important observation as it clearly displays that the misuse of reference genes could lead to the incorrect reporting of a dysregulated genes in cancerous tissue.

Although the selection of the correct reference genes is a key challenge when conducting RT-qPCR studies there are other aspects of experimental design that also need to be considered10. In this study, we highlighted appropriate tests to comply with necessary measures for RT-qPCR studies (Table 2). When utilising relative quantification it is essential that the gene assay of the reference gene and the target gene are amplified with comparable efficiencies34. For this reason, we examined the efficiency of all gene assays using a calibration curve (Figure 4A) and we used this value when evaluating the fold change between conditions. Another important consideration in experimental design is establishing the presence or absence of biological contaminants in samples which may inhibit the RT-qPCR reaction40. We designed an inhibition assay test and displayed that there was no inhibitors present in any of the samples (Figure 4B). Finally, the documenting of the quality assessment of RNA templates is critical within RT-qPCR studies as it has been observed that there is a difference in gene expression stability between intact and degraded RNA samples from the same tissue and higher gene-specific variation in degraded samples34,41. In this study, we documented the RNA purity by the ratio of absorbance at 260/280 nm and RNA quality through visualization of the 28S:18S ribosomal RNA ratio on a 1% agarose gel (Figure 4C,D).

Table 2. Checklist of tests to conduct when designing RT-qPCR studies.

ChecklistSuggested Test
Correct reference genesTest a panel of candidate reference genes using Normfinder
Efficiency of primer
assays
Conduct a calibration curve and use the slope of the graph to
calculate PCR efficiency with the following equation: 10-1/slope-1
Inhibition within samplesAdd samples to a standard RT-qPCR reaction and look for
changes in the Cq values
RNA purityMeasure the ratio of absorbance at 260/280 nm
RNA integrityVisualization of the 28S:18S ribosomal RNA ratio on a 1%
agarose gel

Our data clearly demonstrates that the variability in the expression of reference genes can lead to false interpretation of results; making the selection of the correct genes essential when normalizing RNA concentrations in RT-qPCR analyses. Further to this we have demonstrated appropriate tests to create studies which comply with the MIQE guidelines. The implementation of these guidelines10,42 should be employed by all reviewers when accepting gene expression studies for publication as it will help eliminate the reporting of inaccurate and misleading information.

Data availability

F1000Research: Dataset 1. Cq Values for reference genes in HCT116 cell lines, 10.5256/f1000research.7656.d11180343

F1000Research: Dataset 2. Cq Values for PKC coding genes and reference genes in HCT116 cell lines, 10.5256/f1000research.7656.d11180444

F1000Research: Dataset 3. Cq Values for reference genes in HT29 cell lines, 10.5256/f1000research.7656.d11180545

F1000Research: Dataset 4. Cq Values for PKC coding genes and reference genes in HT29 cell lines, 10.5256/f1000research.7656.d11180746

F1000Research: Dataset 5. Cq Values for reference genes in normal and colon cancer tissue, 10.5256/f1000research.7656.d11180847

F1000Research: Dataset 6. Cq Values for PKC coding genes and reference genes in normal and colon cancer tissue, 10.5256/f1000research.7656.d11180948

F1000Research: Dataset 7. Cq values for sample assay efficiency test, 10.5256/f1000research.7656.d11181049

F1000Research: Dataset 8. Cq values for inhibition assay test, 10.5256/f1000research.7656.d11181150

Consent

Written informed consent for publication of their clinical details and clinical images was obtained from the patients.

(Ethical approval number 73/11, University Hospital Limerick, Limerick, Ireland).

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 22 Jan 2016
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Dowling CM, Walsh D, Coffey JC and Kiely PA. The importance of selecting the appropriate reference genes for quantitative real time PCR as illustrated using colon cancer cells and tissue [version 1; peer review: 2 approved] F1000Research 2016, 5:99 (https://doi.org/10.12688/f1000research.7656.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 22 Jan 2016
Views
28
Cite
Reviewer Report 24 Feb 2016
Verónica Ayllón Cases, Gene Regulation, Stem Cells and Development Group, Department of Genomic Oncology, Centre for Genomics and Oncological Research (Genyo), Pfizer-University of Granada-Regional Government of Andalusia, Granada, Spain 
Approved
VIEWS 28
This article by Dowling et al. demonstrates the importance of choosing the right reference genes (RG) when performing RT-qPCR experiments. They have compared the effect of using a single RG versus three RGs on the gene expression values of the ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Ayllón Cases V. Reviewer Report For: The importance of selecting the appropriate reference genes for quantitative real time PCR as illustrated using colon cancer cells and tissue [version 1; peer review: 2 approved]. F1000Research 2016, 5:99 (https://doi.org/10.5256/f1000research.8245.r12476)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 09 Mar 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    09 Mar 2016
    Author Response
    We thank the reviewer for the suggestions. We have adjusted the manuscript accordingly and we agree that it is better presented now.

    The values given on the datasets corresponds to three ... Continue reading
  • Author Response 09 Mar 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    09 Mar 2016
    Author Response
    We thank the reviewer for the suggestions. We have adjusted the manuscript accordingly and we agree that it is better presented now.

    The values given on the datasets corresponds to three ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 09 Mar 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    09 Mar 2016
    Author Response
    We thank the reviewer for the suggestions. We have adjusted the manuscript accordingly and we agree that it is better presented now.

    The values given on the datasets corresponds to three ... Continue reading
  • Author Response 09 Mar 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    09 Mar 2016
    Author Response
    We thank the reviewer for the suggestions. We have adjusted the manuscript accordingly and we agree that it is better presented now.

    The values given on the datasets corresponds to three ... Continue reading
Views
30
Cite
Reviewer Report 02 Feb 2016
Gary Loughran, School of Biochemistry and Cell Biology, University College Cork, Cork, Ireland 
Approved
VIEWS 30
Dowling et al. reinforce the necessity of using more than one reference gene (RG) for qRT-PCR. They show that not only should more than one RG be used for normalisation but that a panel of RGs should be tested at ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Loughran G. Reviewer Report For: The importance of selecting the appropriate reference genes for quantitative real time PCR as illustrated using colon cancer cells and tissue [version 1; peer review: 2 approved]. F1000Research 2016, 5:99 (https://doi.org/10.5256/f1000research.8245.r12053)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 05 Feb 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    05 Feb 2016
    Author Response
    We thank the reviewer for the suggestion.

    The sequences of the primers and probes we used in our assays are pre-designed and are the proprietary of Life Technologies who are unable ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 05 Feb 2016
    Catríona Dowling, Graduate Entry Medical School, University of Limerick, Limerick, Ireland
    05 Feb 2016
    Author Response
    We thank the reviewer for the suggestion.

    The sequences of the primers and probes we used in our assays are pre-designed and are the proprietary of Life Technologies who are unable ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 22 Jan 2016
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.