Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies

Reference Gene Validation for Quantitative PCR Under Various Biotic and Abiotic Stress Conditions in Toxoptera citricida (Hemiptera, Aphidiae)

The regulation of mRNA expression level is critical for gene expression studies. Currently, quantitative reverse transcription polymerase chain reaction (qRT-PCR) is commonly used to investigate mRNA expression level of genes under various experimental conditions. An important factor that determines the optimal quantification of qRT-PCR data is the choice of the reference gene for normalization. To advance gene expression studies in Toxoptera citricida (Kirkaldy), an important citrus pest and a main vector of the Citrus tristeza virus, we used five tools (GeNorm, NormFinder, BestKeeper, ΔCt methods, and RefFinder) to evaluate seven candidate reference genes (elongation factor-1 alpha [EF1α], beta tubulin [β-TUB], 18S ribosomal RNA [18S], RNA polymerase II large subunit (RNAP II), beta actin (β-ACT), alpha tubulin, and glyceraldhyde-3-phosphate dehydrogenase) under different biotic (developmental stages and wing dimorphism) and abiotic stress (thermal, starvation, and UV irradiation) conditions. The results showed that EF1α and 18S were the most stable genes under various biotic states, β-ACT and β-TUB during thermal stress, EF1α and RNAP II under starvation stress, and RNAP II, β-ACT, and EF1α under UV irradiation stress conditions. This study provides useful resources for the transcriptional profiling of genes in T. citricida and closely related aphid species.

Validation of reference genes for quantitative real-time PCR in Périgord black truffle (Tuber melanosporum) developmental stages

The symbiotic fungus Tuber melanosporum Vittad. (Périgord black truffle) belongs to the Ascomycota and forms mutualistic symbiosis with tree and shrub roots. This truffle has a high value in a global market and is cultivated in many countries of both hemispheres. The publication of the T. melanosporum genome has given researchers unique opportunities to learn more about the biology of the fungus. Real-time quantitative PCR (qRT-PCR) is a definitive technique for quantitating differences in transcriptional gene expression levels between samples. To facilitate gene expression studies and obtain more accurate qRT-PCR data, normalization relative to stable housekeeping genes is required. These housekeeping genes must show stable expression under given experimental conditions for the qRT-PCR results to be accurate. Unfortunately, there are no studies on the stability of housekeeping genes used in T. melanosporum development. In this study, we present a morphological and microscopical classification of the developmental stages of T. melanosporum fruit body, and investigate the expression levels of 12 candidate reference genes (18S rRNA; 5.8S rRNA; Elongation factor 1-alpha; Elongation factor 1-beta; α-tubulin; 60S ribosomal protein L29; β-tubulin; 40S ribosomal protein S1; 40S ribosomal protein S3; Glucose-6-phosphate dehydrogenase; β-actin; Ubiquitin-conjugating enzyme). To evaluate the suitability of these genes as endogenous controls, five software-based approaches and one web-based comprehensive tool (RefFinder) were used to analyze and rank the tested genes. We demonstrate here that the 18S rRNA gene shows the most stable expression during T. melanosporum development and that a set of three genes, 18S rRNA, Elongation factor 1-alpha and 40S ribosomal protein S3, is the most suitable to normalize qRT-PCR data from all the analyzed developmental stages; conversely, 18S rRNA, Glucose-6-phosphate dehydrogenase and Elongation factor 1-alpha are the most suitable genes for fruiting body developmental stages.

Lamina propria macrophage phenotypes in relation to Escherichia coli in Crohn's disease

Background

Abnormal handling of E. coli by lamina propria (LP) macrophages may contribute to Crohn’s disease (CD) pathogenesis. We aimed to determine LP macrophage phenotypes in CD, ulcerative colitis (UC) and healthy controls (HC), and in CD, to compare macrophage phenotypes according to E. coli carriage.

Methods

Mucosal biopsies were taken from 35 patients with CD, 9 with UC and 18 HCs. Laser capture microdissection was used to isolate E. coli-laden and unladen LP macrophages from ileal or colonic biopsies. From these macrophages, mRNA was extracted and cytokine and activation marker expression measured using RT-qPCR.

Results

E. coli-laden LP macrophages were identified commonly in mucosal biopsies from CD patients (25/35, 71 %), rarely in UC (1/9, 11 %) and not at all in healthy controls (0/18). LP macrophage cytokine mRNA expression was greater in CD and UC than healthy controls. In CD, E. coli-laden macrophages expressed high IL-10 & CD163 and lower TNFα, IL-23 & iNOS irrespective of macroscopic inflammation. In inflamed tissue, E. coli-unladen macrophages expressed high TNFα, IL-23 & iNOS and lower IL-10 & CD163. In uninflamed tissue, unladen macrophages had low cytokine mRNA expression, closer to that of healthy controls.

Conclusion

In CD, intra-macrophage E. coli are commonly found and LP macrophages express characteristic cytokine mRNA profiles according to E. coli carriage. Persistence of E. coli within LP macrophages may provide a stimulus for chronic inflammation.

β2-Microglobulin is an appropriate reference gene for RT-PCR-based gene expression analysis of hematopoietic stem cells

Real-time reverse transcription polymerase chain reaction (RT-PCR) is regarded as one of the most useful and powerful tools for characterizing hematopoietic stem cells (HSCs), because samples of extremely small cell numbers can be analyzed. The expression levels determined by RT-PCR are based on relative quantification; therefore, the selection of an appropriate reference gene with a relatively stable expression level under most conditions is crucial. Here, we determined that beta2-microglobulin (B2m) is an appropriate reference gene for analyzing mouse HSCs by a novel method using single-cell RT-PCR. Clonally sorted HSCs were subjected to RT reactions with exogenous RNA fragments and then to real-time PCR. Next, the relative gene expression levels of 4 well-known housekeeping genes were quantified in each single cell sample based on the threshold cycle of exogenous RNA. The analysis revealed that B2m expression was reproducibly detected in almost all HSCs and that B2m had the most stable expression level among the compared genes, even after the cells had been cultured under various conditions. Thus, our results indicate that B2m can reliably be used as a reference gene for the relative quantification of expression levels in HSCs across various conditions. Furthermore, our work proposes a novel method for the selection of appropriate reference genes.

Reference genes for real-time qPCR in leukocytes from asthmatic patients before and after anti-asthma treatment

The aim of this study was to develop a set of reference genes whose expression is stable and suitable for normalization of target gene expression measured in asthma patients during anti-asthmatic treatment. Real-time qPCR was used to determine expression of 7 candidate reference genes (18S rRNA, ACTB, B2M, GAPDH, POLR2A, RPL13A and RPL32) and 7 target genes in leukocytes from asthma patients before and after treatment with inhaled corticosteroids and leukotriene receptor antagonist. Variance of Cq values was analyzed and stability ranking was determined with geNorm. We further investigated how the different normalization strategies affected the consistency of conclusions if the specific investigated target gene is down-regulated or up-regulated after anti-asthmatic therapy. The top-ranking reference genes determined by geNorm, when samples before and after therapy were analyzed (ACTB, B2M and GAPDH) were different from those (POLR2A and B2M) when only samples before treatment were analyzed. Using only a single reference gene for normalization of 7 target gene expression compared to our strategy, there would be as low as 19% of consistency in conclusions. We suggest the use of the geometric mean of ACTB, B2M and GAPDH for normalization of qPCR data of target genes in pharmacogenomics studies in asthma patients before and after anti-asthmatic therapy, however if gene expression is measured only before anti-asthmatic treatment, we recommend the use of the geometric mean of POLR2A and B2M.

Evaluation of potential internal references for quantitative real-time RT-PCR normalization of gene expression in red drum (Sciaenops ocellatus)

Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) has been used extensively for studying gene expression in diverse organisms including fish. In this study, with an aim to identify reliable reference genes for qRT-PCR in red drum (Sciaenops ocellatus), an economic fish species, we determined the expression stability of seven housekeeping genes in healthy and bacterium-infected red drum. Each of the selected candidate genes was amplified by qRT-PCR from the brain, gill, heart, intestine, kidney, liver, muscle, and spleen of red drum infected with or without a bacterial pathogen for 12 and 48 h. The mRNA levels of the genes were analyzed with the geNorm and NormFinder algorithms. The results showed that in the absence of bacterial infection, translation initiation factor 3, NADH dehydrogenase 1, and QM-like protein may be used together as internal references across the eight examined tissues. Bacterial infection caused variations in the rankings of the most stable genes in a tissue-dependent manner. For all tissues, two genes sufficed for reliable normalization at both 12 and 48 h post-infection. However, the optimal gene pairs differed among tissues and, for four of the examined eight tissues, between infection points. These results indicate that when studying gene expression in red drum under conditions of bacterial infection, the optimal reference genes should be selected on the basis of tissue type and, for accurate normalization, infection stage.

Assays for the identification and prioritization of drug candidates for spinal muscular atrophy

Spinal muscular atrophy (SMA) is an autosomal recessive genetic disorder resulting in degeneration of α-motor neurons of the anterior horn and proximal muscle weakness. It is the leading cause of genetic mortality in children younger than 2 years. It affects ∼1 in 11,000 live births. In 95% of cases, SMA is caused by homozygous deletion of the SMN1 gene. In addition, all patients possess at least one copy of an almost identical gene called SMN2. A single point mutation in exon 7 of the SMN2 gene results in the production of low levels of full-length survival of motor neuron (SMN) protein at amounts insufficient to compensate for the loss of the SMN1 gene. Although no drug treatments are available for SMA, a number of drug discovery and development programs are ongoing, with several currently in clinical trials. This review describes the assays used to identify candidate drugs for SMA that modulate SMN2 gene expression by various means. Specifically, it discusses the use of high-throughput screening to identify candidate molecules from primary screens, as well as the technical aspects of a number of widely used secondary assays to assess SMN messenger ribonucleic acid (mRNA) and protein expression, localization, and function. Finally, it describes the process of iterative drug optimization utilized during preclinical SMA drug development to identify clinical candidates for testing in human clinical trials.

Evaluation and selection of reference genes for ecotoxicogenomic study of the green alga Closterium ehrenbergii using quantitative real-time PCR

The green alga Closterium ehrenbergii occurs in fresh water environments and has been suggested as a model for ecotoxicological assessment. Quantitative real-time PCR (qRT-PCR), with its high sensitivity and specificity, is a preferred method for reliable quantification of gene expression levels. qRT-PCR requires reference genes to normalize the transcription level of the target gene, and selection of appropriate references is crucial. Here, we evaluated nine housekeeping genes, that is, 18S rRNA, ACT, TUA, TUB, eIF, H4, UBQ, rps4, and GAPDH, using 34 RNA samples of C. ehrenbergii cultured in various environments (e.g. exposure to heat shock, UV, metals, and non-metallic chemicals). Each housekeeping gene tested displayed different ranges of C T values for each experimental condition. The gene stability was determined using the descriptive statistic software geNorm, which showed that ACT, H4, and TUA were the most suitable reference genes for all the conditions tested. In addition, at least three genes were required for proper normalization. With these references, we assessed the expression level of the heat shock protein 70 (HSP70) gene in C. ehrenbergii cells exposed to thermal and toxic contaminant stress and found that it was significantly up-regulated by these stressors. This study provides potential reference genes for gene expression studies on C. ehrenbergii with qRT-PCR.

Validation of housekeeping genes for gene expression analysis in glioblastoma using quantitative real-time polymerase chain reaction

Background

Quantitative real-time polymerase chain reaction (qPCR) is the most reliable tool for gene expression studies. Selection of housekeeping genes (HKGs) that are having most stable expression is critical to carry out accurate gene expression profiling. There is no 'universal' HKG having stable expression in all kinds of tissues under all experimental conditions.

Methods

The present study aims to identify most appropriate HKGs for gene expression analysis in glioblastoma (GBM) samples. Based on literature survey, six most commonly used HKGs that are invariant in GBM were chosen. We performed qPCR using RNA from formalin fixed paraffin embedded GBM samples and normal brain samples to investigate the expression pattern of HPRT, GAPDH, TBP, B2M, B2M, RPL13A, and RN18S1 with different abundance. A simple Δcycle threshold approach was employed to calculate the fold change.

Results

Our study shows that the expression of RPL13A and TBP were found to be most stable across all the samples and are thus suitable for gene expression analysis in human GBM. Except for TBP, none of the other conventionally used HKGs in GBM studies e.g., HPRT and GAPDH were found to be suitable as they showed variation in RNA expression.

Conclusion

Validation of HKGs is therefore immensely specific for a particular experimental setup and is crucial in assessing any new setup.

Relative expression of rRNA transcripts and 45S rDNA promoter methylation status are dysregulated in tumors in comparison with matched-normal tissues in breast cancer

Ribosomal RNA (rRNA) expression, one of the most important factors regulating ribosome production, is primarily controlled by a CG-rich 45 S rDNA promoter. However, the DNA methylation state of the 45 S rDNA promoter, as well as its effect on rRNA gene expression in types of human cancers is controversial. In the present study we analyzed the methylation status of the rDNA promoter (-380 to +53 bp) as well as associated rRNA expression levels in breast cancer cell lines and breast tumor-normal tissue pairs. We found that the aforementioned regulatory region was extensively methylated (74-96%) in all cell lines and in 68% (13/19 tumor-normal pairs) of the tumors. Expression levels of rRNA transcripts 18 S, 28 S, 5.8 S and 45 S external transcribed spacer (45 S ETS) greatly varied in the breast cancer cell lines regardless of their methylation status. Analyses of rRNA transcript expression levels in the breast tumor and normal matched tissues showed no significant difference when normalized with TBP. On the other hand, using the geometric mean of the rRNA expression values (GM-rRNA) as reference enabled us to identify significant changes in the relative expression of rRNAs in the tissue samples. We propose GM-rRNA normalization as a novel strategy to analyze expression differences between rRNA transcripts. Accordingly, the 18S rRNA/GM-rRNA ratio was significantly higher whereas the 5.8S rRNA/GM-rRNA ratio was significantly lower in breast tumor samples than this ratio in the matched normal samples. Moreover, the 18S rRNA/GM-rRNA ratio was negatively correlated with the 45 S rDNA promoter methylation level in the normal breast tissue samples, yet not in the breast tumors. Significant correlations observed between the expression levels of rRNA transcripts in the normal samples were lost in the tumor samples. We showed that the expression of rRNA transcripts may not be based solely on promoter methylation. Carcinogenesis may cause dysregulation of the correlation between spliced rRNA expression levels, possibly due to changes in rRNA processing, which requires further investigation.

Comparative proteomics: assessment of biological variability and dataset comparability

Background

Comparative proteomics in bacteria are often hampered by the differential nature of dataset quality and/or inherent biological deviations. Although common practice compensates by reproducing and normalizing datasets from a single sample, the degree of certainty is limited in comparison of multiple dataset. To surmount these limitations, we introduce a two-step assessment criterion using: (1) the relative number of total spectra (R_TS) to determine if two LC-MS/MS datasets are comparable and (2) nine glycolytic enzymes as internal standards for a more accurate calculation of relative amount of proteins. Lactococcus lactis HR279 and JHK24 strains expressing high or low levels (respectively) of green fluorescent protein (GFP) were used for the model system. GFP abundance was determined by spectral counting and direct fluorescence measurements. Statistical analysis determined relative GFP quantity obtained from our approach matched values obtained from fluorescence measurements.

Results

L. lactis HR279 and JHK24 demonstrates two datasets with an R_TS value less than 1.4 accurately reflects relative differences in GFP levels between high and low expression strains. Without prior consideration of R_TS and the use of internal standards, the relative increase in GFP calculated by spectral counting method was 3.92 ± 1.14 fold, which is not correlated with the value determined by the direct fluorescence measurement (2.86 ± 0.42 fold) with the p = 0.024. In contrast, 2.88 ± 0.92 fold was obtained by our approach showing a statistically insignificant difference (p = 0.95).

Conclusions

Our two-step assessment demonstrates a useful approach to: (1) validate the comparability of two mass spectrometric datasets and (2) accurately calculate the relative amount of proteins between proteomic datasets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0561-9) contains supplementary material, which is available to authorized users.

Validation of endogenous reference genes in Buglossoides arvensis for normalizing RT-qPCR-based gene expression data

Selection of a stably expressed reference gene (RG) is an important step for generating reliable and reproducible quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) gene expression data. We, in this study, have sought to validate RGs for Buglossoides arvensis, a high nutraceutical value plant whose refined seed oil is entering the market under the commercial trade name Ahiflower™. This weed plant has received attention for its natural ability to significantly accumulate the poly-unsaturated fatty acid (PUFA) stearidonic acid (SDA, C18:4n-3) in its seeds, which is uncommon for most plant species. Ten candidate RGs (β-Act, 18S rRNA, EF-1a, α-Tub, UBQ, α-actin, CAC, PP2a, RUBISCO, GAPDH) were isolated from B. arvensis and TaqMan™ compliant primers/probes were designed for RT-qPCR analysis. Abundance of these gene transcripts was analyzed across different tissues and growth regimes. Two of the most widely used algorithms, geNorm and NormFinder, showed variation in expression levels of these RGs. However, combinatorial analysis of the results clearly identified CAC and α-actin as the most stable and unstable RG candidates, respectively. This study has for the first time identified and validated RGs in the non-model system B. arvensis, a weed plant projected to become an important yet sustainable source of dietary omega-3 PUFA.

Selection of suitable reference genes for expression analysis in human glioma using RT-qPCR

In human glioma research, quantitative real-time reverse-transcription PCR is a frequently used tool. Considering the broad variation in the expression of candidate reference genes among tumor stages and normal brain, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. This study aimed at testing a panel of nine reference genes [beta-2-microglobulin, cytochrome c-1 (CYC1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hydroxymethylbilane synthase, hypoxanthine guanine phosphoribosyl transferase 1, ribosomal protein L13a (RPL13A), succinate dehydrogenase, TATA-box binding protein and 14-3-3 protein zeta] to identify and validate the most suitable reference genes for expression studies in human glioma of different grades (World Health Organization grades II-IV). After analysis of the stability values calculated using geNorm, NormFinder, and BestKeeper algorithms, GAPDH, RPL13A, and CYC1 can be indicated as reference genes applicable for accurate normalization of gene expression in glioma compared with normal brain and anaplastic astrocytoma or glioblastoma alone within this experimental setting. Generally, there are no differences in expression levels and variability of candidate genes in glioma tissue compared to normal brain. But stability analyses revealed just a small number of genes suitable for normalization in each of the tumor subgroups and across these groups. Nevertheless, our data show the importance of validation of adequate reference genes prior to every study.

Evaluation of housekeeping genes for normalizing real-time quantitative PCR assays in pig skeletal muscle at multiple developmental stages

Quantitative PCR (qPCR) requires a constantly expressed housekeeping gene as an internal control, the expression of which is similar in different biological samples. In the present study, we evaluated the applicability and compared the consistency of the gene expression of 16 reference genes, i.e., 10 common and 6 candidate genes, through qPCR assays in pig skeletal muscles at multiple developmental stages. The stability of these 16 potential reference genes was investigated using the geNorm and NormFinder methods. Our results indicated that DRAP1 and RNF7 were the most appropriate combination to normalize gene expression in the Yorkshire samples, the combination of DRAP1 and WSB2 were appropriate in the Tongcheng samples, H3F3A and DRAP1 in prenatal periods, DRAP1 and RNF7 in postnatal periods, and the combination of DRAP1 and WSB2 was most suitable for accurate normalization in whole samples. This result provides a strong basis for similar studies in pigs that explore the molecular mechanisms underlying skeletal muscle development.

An analysis of critical factors for quantitative immunoblotting

Immunoblotting (also known as Western blotting) combined with digital image analysis can be a reliable method for analyzing the abundance of proteins and protein modifications, but not every immunoblot-analysis combination produces an accurate result. I illustrate how sample preparation, protocol implementation, detection scheme, and normalization approach profoundly affect the quantitative performance of immunoblotting. This study implemented diagnostic experiments that assess an immunoblot-analysis workflow for accuracy and precision. The results showed that ignoring such diagnostics can lead to pseudoquantitative immunoblot data that markedly overestimate or underestimate true differences in protein abundance.

Classifying ten types of major cancers based on reverse phase protein array profiles

Gathering vast data sets of cancer genomes requires more efficient and autonomous procedures to classify cancer types and to discover a few essential genes to distinguish different cancers. Because protein expression is more stable than gene expression, we chose reverse phase protein array (RPPA) data, a powerful and robust antibody-based high-throughput approach for targeted proteomics, to perform our research. In this study, we proposed a computational framework to classify the patient samples into ten major cancer types based on the RPPA data using the SMO (Sequential minimal optimization) method. A careful feature selection procedure was employed to select 23 important proteins from the total of 187 proteins by mRMR (minimum Redundancy Maximum Relevance Feature Selection) and IFS (Incremental Feature Selection) on the training set. By using the 23 proteins, we successfully classified the ten cancer types with an MCC (Matthews Correlation Coefficient) of 0.904 on the training set, evaluated by 10-fold cross-validation, and an MCC of 0.936 on an independent test set. Further analysis of these 23 proteins was performed. Most of these proteins can present the hallmarks of cancer; Chk2, for example, plays an important role in the proliferation of cancer cells. Our analysis of these 23 proteins lends credence to the importance of these genes as indicators of cancer classification. We also believe our methods and findings may shed light on the discoveries of specific biomarkers of different types of cancers.

Identification of suitable endogenous controls for gene and miRNA expression studies in irradiated prostate cancer cells

This study aimed to to evaluate the stability of commonly used endogenous control genes for messenger RNA (mRNA) (N = 16) and miRNAs (N = 3) expression studies in prostate cell lines following irradiation. The stability of endogenous control genes expression in irradiated (6 Gy) versus unirradiated controls was quantified using NormFinder and coefficient of variation analyses. HPRT1 and 18S were identified as most and least stable endogenous controls, respectively, for mRNA expression studies in irradiated prostate cells. SNORD48 and miR16 miRNA endogenous controls tested were associated with low coefficient of variations following irradiation (6 Gy). This study highlights that commonly used endogenous controls can be responsive to radiation and validation is required prior to gene/miRNAs expression studies.

Gapdh gene expression is modulated by inflammatory arthritis and is not suitable for qPCR normalization

Gene expression studies are fundamental for the understanding of complex diseases, providing new insights into the pathogenic process and new tools for diagnostic and patient stratification. Gene profiling studies by real-time PCR require the use of reference genes for normalization and an appropriate validation is essential for accurate results. We performed a comprehensive assessment of six common housekeeping genes in the K/BxN serum-induced arthritis model in mice. Classical statistics and NormFinder analyses pointed out Gapdh as the less stable and therefore unsuitable as a reference control. Gapdh was considerably down-regulated in arthritic joints and therefore produced an overestimation of transcriptional changes. Hptr, B2m, and Rpl13a showed the most constant expression. Collectively our data advise against the use of Gapdh in gene expression studies in the acute phase of the K/BxN model and adds a cautionary note on the need to validate the reference genes for reliable, comparable, and reproducible results.

Exogenous reference gene normalization for real-time reverse transcription-polymerase chain reaction analysis under dynamic endogenous transcription

Quantitative real-time reverse transcription-polymerase chain reaction (qPCR) is widely used to investigate transcriptional changes following experimental manipulations to the nervous system. Despite the widespread utilization of qPCR, the interpretation of results is marred by the lack of a suitable reference gene due to the dynamic nature of endogenous transcription. To address this inherent deficiency, we investigated the use of an exogenous spike-in mRNA, luciferase, as an internal reference gene for the 2(-∆∆Ct) normalization method. To induce dynamic transcription, we systemically administered capsaicin, a neurotoxin selective for C-type sensory neurons expressing the TRPV-1 receptor, to adult male Sprague-Dawley rats. We later isolated nodose ganglia for qPCR analysis with the reference being either exogenous luciferase mRNA or the commonly used endogenous reference β-III tubulin. The exogenous luciferase mRNA reference clearly demonstrated the dynamic expression of the endogenous reference. Furthermore, variability of the endogenous reference would lead to misinterpretation of other genes of interest. In conclusion, traditional reference genes are often unstable under physiologically normal situations, and certainly unstable following the damage to the nervous system. The use of exogenous spike-in reference provides a consistent and easily implemented alternative for the analysis of qPCR data.

Selection of reference genes for qPCR- and ddPCR-based analyses of gene expression in Senescing Barley leaves

Leaf senescence is a tightly regulated developmental or stress-induced process. It is accompanied by dramatic changes in cell metabolism and structure, eventually leading to the disintegration of chloroplasts, the breakdown of leaf proteins, internucleosomal fragmentation of nuclear DNA and ultimately cell death. In light of the global and intense reorganization of the senescing leaf transcriptome, measuring time-course gene expression patterns in this model is challenging due to the evident problems associated with selecting stable reference genes. We have used oligonucleotide microarray data to identify 181 genes with stable expression in the course of dark-induced senescence of barley leaf. From those genes, we selected 5 candidates and confirmed their invariant expression by both reverse transcription quantitative PCR and droplet digital PCR (ddPCR). We used the selected reference genes to normalize the level of the expression of the following senescence-responsive genes in ddPCR assays: SAG12, ICL, AGXT, CS and RbcS. We were thereby able to achieve a substantial reduction in the data variability. Although the use of reference genes is not considered mandatory in ddPCR assays, our results show that it is advisable in special cases, specifically those that involve the following conditions: i) a low number of repeats, ii) the detection of low-fold changes in gene expression or iii) series data comparisons (such as time-course experiments) in which large sample variation greatly affects the overall gene expression profile and biological interpretation of the data.

Identification of Stably Expressed lncRNAs as Valid Endogenous Controls for Profiling of Human Glioma

Background: Recent research indicates that long non-coding RNAs (lncRNA) represent a new family of RNAs that is of fundamental importance for controlling transcription and translation. Thereby, there is increasing evidence that lncRNAs are also important in tumourigenesis. Thereby valid expression profiling using quantitative PCR requires suitable, stably expressed normalisers to achieve reliable and reproducible data. However, no systematic analysis of suitable references in lncRNA studies in human glioma has been performed yet.

Methods: In this study, we investigated 90 lncRNAs in 30 tissue specimen for the expression stability in human diffuse astrocytoma (WHO-Grade II), anaplastic astrocytoma (WHO-Grade III) and glioblastoma (WHO-Grade IV) both alone as well as in comparison with normal white matter. Our identification procedure included a rigorous bioinformatical selection process that resulted in the inclusion of only highly abundant, equally expressed lncRNAs for further analysis. Additionally, lncRNAs were classified according to their stability value using the NormFinder algorithm.

Results: We identified 24 appropriate normalisers suitable for studies in diffuse astrocytoma, 22 for studies in anaplastic astrocytoma and 12 for studies in glioblastoma. Comparing all three glioma entities 7 lncRNAs showed stable expression levels. Addition of normal brain tissue resulted in only 4 suitable lncRNAs.

Conclusions: Our findings indicate that 4 lncRNAs (HOXA6as, H19 upstream conserved 1 and 2, Zfhx2as and BC200) are suitable as normalisers in glioma and normal brain. These lncRNAs may thus be regarded as universal references being applicable for the accurate normalisation of lncRNA expression profiling in various glioma (WHO-Grades II-IV) alone and in combination with brain tissue. This enables to perform valid longitudinal studies, e.g. of glioma before and after malignisation to identify changes of lncRNA expressions probably driving malignant transformation.

Cytokine gene expression profiles in chicken spleen and intestinal tissues during Ascaridia galli infection

In the poultry production industry, chickens with access to outdoor areas are exposed to a wide range of parasites e.g. the helminth Ascaridia galli. By real-time quantitative RT-PCR, the relative gene expression of the T helper 1 (Th1) cytokine IFN-γ, the T helper 2 (Th2) cytokine IL-13, the anti-inflammatory cytokines IL-10 and TGF-β4 and the pro-inflammatory cytokine IL-17F were determined over a period of 3 weeks in A. galli and non-A. galli-infected chickens. A characteristic Th2 response was observed in the jejunum of A. galli-infected chickens with increased expression of IL-13 and decreased expression of IFN-γ from day 14 post infection. At the putative time of larvae invasion into the intestinal mucosa (day 7), an increased expression of IFN-γ, IL-10, and TGF-β4 was observed in the spleen. At the putative onset of the innate immune response (day 10), a decreased expression of jejunal IFN-γ and IL-13 was observed. Finally, at the expected period of an adaptive immune response (days 14-21) a general decreased expression of IFN-γ and TGF-β4 in spleen and IFN-γ in jejunum was followed by a decreased expression of IFN-γ and IL-10 at day 21 in caecal tonsils.

Identification of optimal reference genes for gene expression normalization in a wide cohort of endometrioid endometrial carcinoma tissues

Accurate normalization is a primary component of a reliable gene expression analysis based on qRT-PCR technique. While the use of one or more reference genes as internal controls is commonly accepted as the most appropriate normalization strategy, many qPCR-based published studies still contain data poorly normalized and reference genes arbitrarily chosen irrespective of the particular tissue and the specific experimental design. To date, no validated reference genes have been identified for endometrial cancer tissues. In this study, 10 normalization genes (GAPDH, B2M, ACTB, POLR2A, UBC, PPIA, HPRT1, GUSB, TBP, H3F3A) belonging to different functional and abundance classes in various tissues and used in different studies, were analyzed to determine their applicability. In total, 100 endometrioid endometrial cancer samples, which were carefully balanced according to their tumor grade, and 29 normal endometrial tissues were examined using SYBR Green Real-Time RT-PCR. The expression stability of candidate reference genes was determined and compared by means of geNorm and NormFinder softwares. Both algorithms were in agreement in identifying GAPDH, H3F3A, PPIA, and HPRT1 as the most stably expressed genes, only differing in their ranking order. Analysis performed on the expression levels of all candidate genes confirm HPRT1 and PPIA as the most stably expressed in the study groups regardless of sample type, to be used alone or better in combination. As the stable expression of HPRT1 and PPIA between normal and tumor endometrial samples fulfill the basic requirement of a reference gene to be used for normalization purposes, HPRT1 expression showed significant differences between samples from low-grade and high-grade tumors. In conclusion, our results recommend the use of PPIA as a single reference gene to be considered for improved reliability of normalization in gene expression studies involving endometrial tumor samples at different tumor degrees.

Quantifying mRNA and microRNA with qPCR in cervical carcinogenesis: a validation of reference genes to ensure accurate data

A number of recent studies have catalogued global gene expression patterns in a panel of normal, tumoral cervical tissues so that potential biomarkers can be identified. The qPCR has been one of the most widely used technologies for detecting these potential biomarkers. However, few studies have investigated a correct strategy for the normalization of data in qPCR assays for cervical tissues. The aim of this study was to validate reference genes in cervical tissues to ensure accurate quantification of mRNA and miRNA levels in cervical carcinogenesis. For this purpose, some issues for obtaining reliable qPCR data were evaluated such as the following: geNorm analysis with a set of samples which meet all of the cervical tissue conditions (Normal + CIN1 + CIN2 + CIN3 + Cancer); the use of individual Ct values versus pooled Ct values; and the use of a single (or multiple) reference genes to quantify mRNA and miRNA expression levels. Two different data sets were put on the geNorm to assess the expression stability of the candidate reference genes: the first dataset comprised the quantities of the individual Ct values; and the second dataset comprised the quantities of the pooled Ct values. Moreover, in this study, all the candidate reference genes were analyzed as a single “normalizer”. The normalization strategies were assessed by measuring p16^INK4a and miR-203 transcripts in qPCR assays. We found that the use of pooled Ct values, can lead to a misinterpretation of the results, which suggests that the maintenance of inter-individual variability is a key factor in ensuring the reliability of the qPCR data. In addition, it should be stressed that a proper validation of the suitability of the reference genes is required for each experimental setting, since the indiscriminate use of a reference gene can also lead to discrepant results.

Expressed repeat elements improve RT-qPCR normalization across a wide range of zebrafish gene expression studies

The selection and validation of stably expressed reference genes is a critical issue for proper RT-qPCR data normalization. In zebrafish expression studies, many commonly used reference genes are not generally applicable given their variability in expression levels under a variety of experimental conditions. Inappropriate use of these reference genes may lead to false interpretation of expression data and unreliable conclusions. In this study, we evaluated a novel normalization method in zebrafish using expressed repetitive elements (ERE) as reference targets, instead of specific protein coding mRNA targets. We assessed and compared the expression stability of a number of EREs to that of commonly used zebrafish reference genes in a diverse set of experimental conditions including a developmental time series, a set of different organs from adult fish and different treatments of zebrafish embryos including morpholino injections and administration of chemicals. Using geNorm and rank aggregation analysis we demonstrated that EREs have a higher overall expression stability compared to the commonly used reference genes. Moreover, we propose a limited set of ERE reference targets (hatn10, dna15ta1 and loopern4), that show stable expression throughout the wide range of experiments in this study, as strong candidates for inclusion as reference targets for qPCR normalization in future zebrafish expression studies. Our applied strategy to find and evaluate candidate expressed repeat elements for RT-qPCR data normalization has high potential to be used also for other species.

Validation of four reference genes for quantitative mRNA expression studies in a rat model of inflammatory injury

BACKGROUND

Real-time quantitative PCR (qPCR) is a technique frequently used to measure changes in mRNA expression. To ensure validity of experimental findings, it is important to normalize the qPCR data to reference genes that are stable and unaffected by the experimental treatment to correct for variability among samples. Unlike in some models of neuropathic pain, reference genes for models of inflammatory injury have not been validated. This study examined four candidate reference genes in an effort to identify and validate optimal genes for normalization of transcriptional changes occurring in the dorsal horn of the spinal cord and the rostral ventromedial medulla (RVM) following intraplantar injection of complete Freund's adjuvant (CFA).

RESULTS

The expression of hypoxanthine phosphoribosyltransferase 1 (Hprt1), beta-actin (Actb), mitogen-activated protein kinase 6 (Mapk6), and beta-2-microglobulin (B2m) was quantified in the dorsal horn and RVM of rats four days or two weeks after intraplantar injection of CFA or saline. The range of expression levels among these four genes differed by as much as 16-fold within the dorsal horn and the RVM. All four of these reference genes were stably expressed in both tissues and did not differ between saline and CFA-treated animals. Analyses using the statistical algorithms in geNorm and NormFinder programs determined that Mapk6 was the most stable gene and recommended the combination of Mapk6 and Actb, or Mapk6 and Hprt1, in such experimental conditions.

CONCLUSIONS

This study validated the four genes Hprt1, Actb, Mapk6 or B2m and showed that any one or combination of two of them are good reference genes for normalization of mRNA expression in qPCR experiments in the spinal cord and RVM in the CFA model of inflammatory injury.

Selection of reference genes for gene expression studies in Siberian Apricot (Prunus sibirica L.) Germplasm using quantitative real-time PCR

Quantitative real time reverse transcription polymerase chain reaction has been applied in a vast range of studies of gene expression analysis. However, real-time PCR data must be normalized with one or more reference genes. In this study, eleven putative consistently expressed genes (ACT, TUA, TUB, CYP, DNAj, ELFA, F-box27, RPL12, GAPDH, UBC and UBQ) in nine Siberian Apricot Germplasms (including much variability) were evaluated for their potential as references for the normalization of gene expression by NormFinder and geNorm programs. From our studies, ACT, UBC, CYP, UBQ and RPL12 as suitable for normalization were identified by geNorm, while UBC and CYP as the best pair by NormFinder. Moreover, UBC was selected as the most stably expressed gene by both algorithms in different Siberian Apricot seed samples. We also detected that a set of three genes (ACT, CYP and UBC) by geNorm as control for normalization could lead to accurate results. Furthermore, the expression levels of oleosin gene were analyzed to validate the suitability of the selected reference genes. These obtained experimental results could make an important contribution to normalize real-time PCR data for gene expression analysis in Siberian Apricot Germplasm.

Selection of reference genes as internal controls for gene expression in tissues of red abalone Haliotis rufescens (Mollusca, Vetigastropoda; Swainson, 1822)

The red abalone Haliotis rufescens is one of the most important species for aquaculture in Baja California, México, and despite this, few gene expression studies have been done in tissues such as gill, head and gonad. For this purpose, reverse transcription and quantitative real time PCR (RT-qPCR) is a powerful tool for gene expression evaluation. For a reliable analysis, however, it is necessary to select and validate housekeeping genes that allow proper transcription quantification. Stability of nine housekeeping genes (ACTB, BGLU, TUBB, CY, GAPDH, HPRTI, RPL5, SDHA and UBC) was evaluated in different tissues of red abalone (gill, head and gonad/digestive gland). Four-fold serial dilutions of cDNA (from 25 ngμL(-1) to 0.39 ngμL(-1)) were used to prepare the standard curve, and it showed gene efficiencies between 0.95 and 0.99, with R(2)=0.99. geNorm and NormFinder analysis showed that RPL5 and CY were the most stable genes considering all tissues, whereas in gill HPRTI and BGLU were most stable. In gonad/digestive gland, RPL5 and TUBB were the most stable genes with geNorm, while SDHA and HPRTI were the best using NormFinder. Similarly, in head the best genes were RPL5 and UBC with geNorm, and GAPDH and CY with NormFinder. The technical variability analysis with RPL5 and abalone gonad/digestive gland tissue indicated a high repeatability with a variation coefficient within groups ≤ 0.56% and between groups ≤ 1.89%. These results will help us for further research in reproduction, thermoregulation and endocrinology in red abalone.

Identification of suitable reference genes for the relative quantification of microRNAs in pleural effusion

Circulating cell-free microRNAs (miRNAs) are potential biomarkers of cancer. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is widely used in miRNA expression studies. The aim of this study was to identify suitable reference genes for RT-qPCR analyses of miRNA expression levels in pleural effusion. The expression levels of candidate reference miRNAs were investigated in 10 benign pleural effusion (BPE) and 10 lung adenocarcinoma-associated malignant pleural effusion (LA-MPE) samples using miRNA microarrays. The expression levels of candidate reference miRNAs, together with those of U6 small nuclear RNA (snRNA), RNU6B, RNU44 and RNU48 small RNAs, in 46 BPE and 45 LA-MPE samples were validated by RT-qPCR, and were analyzed using the NormFinder and BestKeeper algorithms. The impact of different normalization approaches on the detection of differential expression levels of miR-198 in BPE and LA-MPE samples was also assessed. As determined by the miRNA microarray data, five candidate reference miRNAs were identified. Following RT-qPCR validation, U6 snRNA, miR-192, miR-20a, miR-221, miR-222 and miR-16 were evaluated using the NormFinder and BestKeeper software programs. U6 snRNA and miR-192 were identified as single reference genes and the combination of these genes was preferred for the relative quantification of miRNA expression levels in pleural effusion. Normalization of miR-98 expression levels to those of U6 snRNA, miR-192 or a combination of these genes enabled the detection of a significant difference between BPE and LA-MPE samples. Therefore, U6 snRNA and miR-192 are recommended as reference genes for the relative quantification of miRNA expression levels in pleural effusion.

Reliable reference gene selection for Cordyceps militaris gene expression studies under different developmental stages and media

Cordyceps militaris is considered a model organism for the study of Cordyceps species, which are highly prized in traditional Chinese medicine. Gene expression analysis has become more popular and important in studies of this fungus. Reference gene validation under different experimental conditions is crucial for RT-qPCR analysis. In this study, eight candidate reference genes, actin, cox5, gpd, rpb1, tef1, try, tub, and ubi, were selected and their expression stability was evaluated in C. militaris samples using four algorithms, genorm, normfinder, bestkeeper, and the comparative ∆Ct method. Three sets of samples, five different developmental stages cultured in wheat medium and pupae, and all the samples pool were included. The results showed that rpb1 was the best reference gene during all developmental stages examined, while the most common reference genes, actin and tub, were not suitable internal controls. Cox5 also performed poorly and was less stable in our analysis. The ranks of ubi and gpd were inconsistent in different sample sets by different methods. Our results provide guidelines for reference gene selection at different developmental stages and also represent a foundation for more accurate and widespread use of RT-qPCR in C. militaris gene expression analysis.

Identification and validation of suitable reference genes for RT-qPCR analysis in mouse testis development

RT-qPCR is a commonly used method for evaluating gene expression; however, its accuracy and reliability are dependent upon the choice of appropriate reference gene(s), and there is limited information available on suitable reference gene(s) that can be used in mouse testis at different stages. In this study, using the RT-qPCR method, we investigated the expression variations of six reference genes representing different functional classes (Actb, Gapdh, Ppia, Tbp, Rps29, Hprt1) in mice testis during embryonic and postnatal development. The expression stabilities of putative reference genes were evaluated using five algorithms: geNorm, NormFinder, Bestkeeper, the comparative delta C(t) method and integrated tool RefFinder. Analysis of the results showed that Ppia, Gapdh and Actb were identified as the most stable genes and the geometric mean of Ppia, Gapdh and Actb constitutes an appropriate normalization factor for gene expression studies. The mRNA expression of AT1 as a test gene of interest varied depending upon which of the reference gene(s) was used as an internal control(s). This study suggested that Ppia, Gapdh and Actb are suitable reference genes among the six genes used for RT-qPCR normalization and provide crucial information for transcriptional analyses in future studies of gene expression in the developing mouse testis.

Reference genes for quantitative gene expression studies in multiple avian species

Quantitative real-time PCR (qPCR) rapidly and reliably quantifies gene expression levels across different experimental conditions. Selection of suitable reference genes is essential for meaningful normalization and thus correct interpretation of data. In recent years, an increasing number of avian species other than the chicken has been investigated molecularly, highlighting the need for an experimentally validated pan-avian primer set for reference genes. Here we report testing a set for 14 candidate reference genes (18S, ABL, GAPDH, GUSB, HMBS, HPRT, PGK1, RPL13, RPL19, RPS7, SDHA, TFRC, VIM, YWHAZ) on different tissues of the mallard (Anas platyrhynchos), domestic chicken (Gallus gallus domesticus), common crane (Grus grus), white-tailed eagle (Haliaeetus albicilla), domestic turkey (Meleagris gallopavo f. domestica), cockatiel (Nymphicus hollandicus), Humboldt penguin (Sphenicus humboldti), ostrich (Struthio camelus) and zebra finch (Taeniopygia guttata), spanning a broad range of the phylogenetic tree of birds. Primer pairs for six to 11 genes were successfully established for each of the nine species. As a proof of principle, we analyzed expression levels of 10 candidate reference genes as well as FOXP2 and the immediate early genes, EGR1 and CFOS, known to be rapidly induced by singing in the avian basal ganglia. We extracted RNA from microbiopsies of the striatal song nucleus Area X of adult male zebra finches after they had sang or remained silent. Using three different statistical algorithms, we identified five genes (18S, PGK1, RPS7, TFRC, YWHAZ) that were stably expressed within each group and also between the singing and silent conditions, establishing them as suitable reference genes. In conclusion, the newly developed pan-avian primer set allows accurate normalization and quantification of gene expression levels in multiple avian species.

Selection and assessment of reference genes for quantitative PCR normalization in migratory locust Locusta migratoria (Orthoptera: Acrididae)

Locusta migratoria is a classic hemimetamorphosis insect and has caused widespread economic damage to crops as a migratory pest. Researches on the expression pattern of functional genes in L. migratoria have drawn focus in recent years, especially with the release of genome information. Real-time quantitative PCR is the most reproducible and sensitive approach for detecting transcript expression levels of target genes, but optimal internal standards are key factors for its accuracy and reliability. Therefore, it's necessary to provide a systematic stability assessment of internal control for well-performed tests of target gene expression profile. In this study, twelve candidate genes (Ach, Act, Cht2, EF1α, RPL32, Hsp70, Tub, RP49, SDH, GAPDH, 18S, and His) were analyzed with four statistical methods: the delta Ct approach, geNorm, Bestkeeper and NormFinder. The results from these analyses aimed to choose the best suitable reference gene across different experimental situations for gene profile study in L. migratoria. The result demonstrated that for different developmental stages, EF1α, Hsp70 and RPL32 exhibited the most stable expression status for all samples; EF1α and RPL32 were selected as the best reference genes for studies involving embryo and larvae stages, while SDH and RP49 were identified for adult stage. The best-ranked reference genes across different tissues are RPL32, Hsp70 and RP49. For abiotic treatments, the most appropriate genes we identified were as follows: Act and SDH for larvae subjected to different insecticides; RPL32 and Ach for larvae exposed to different temperature treatments; and Act and Ach for larvae suffering from starvation. The present report should facilitate future researches on gene expression in L. migratoria with accessibly optimal reference genes under different experimental contexts.

Identification and validation of reference genes for quantitative real-time PCR normalization and its applications in lycium

Lycium barbarum and L. ruthenicum are extensively used as traditional Chinese medicinal plants. Next generation sequencing technology provides a powerful tool for analyzing transcriptomic profiles of gene expression in non-model species. Such gene expression can then be confirmed with quantitative real-time polymerase chain reaction (qRT-PCR). Therefore, use of systematically identified suitable reference genes is a prerequisite for obtaining reliable gene expression data. Here, we calculated the expression stability of 18 candidate reference genes across samples from different tissues and grown under salt stress using geNorm and NormFinder procedures. The geNorm-determined rank of reference genes was similar to those defined by NormFinder with some differences. Both procedures confirmed that the single most stable reference gene was ACNTIN1 for L. barbarum fruits, H2B1 for L. barbarum roots, and EF1α for L. ruthenicum fruits. PGK3, H2B2, and PGK3 were identified as the best stable reference genes for salt-treated L. ruthenicum leaves, roots, and stems, respectively. H2B1 and GAPDH1+PGK1 for L. ruthenicum and SAMDC2+H2B1 for L. barbarum were the best single and/or combined reference genes across all samples. Finally, expression of salt-responsive gene NAC, fruit ripening candidate gene LrPG, and anthocyanin genes were investigated to confirm the validity of the selected reference genes. Suitable reference genes identified in this study provide a foundation for accurately assessing gene expression and further better understanding of novel gene function to elucidate molecular mechanisms behind particular biological/physiological processes in Lycium.

The choice of reference gene affects statistical efficiency in quantitative PCR data analysis

Quantitative polymerase chain reaction (qPCR), a highly sensitive method of measuring gene expression, is widely used in biomedical research. To produce reliable results, it is essential to use stably expressed reference genes (RGs) for data normalization so that sample-to-sample variation can be controlled. In this study, we examine the effect of different RGs on statistical efficiency by analyzing a qPCR data set that contains 12 target genes and 3 RGs. Our results show that choosing the most stably expressed RG for data normalization does not guarantee reduced variance or improved statistical efficiency. We also provide a formula for determining when data normalization will improve statistical efficiency and hence increase the power of statistical tests in data analysis.

Selection of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in Brassica napus under various stress conditions

Data normalization is essential for reliable output of quantitative real-time reverse-transcription polymerase chain reaction (qRT-PCR) assays, as the unsuitable choice of reference gene(s), whose expression might be influenced by exogenous treatments in plant tissues, could cause misinterpretation of results. To date, no systematic studies on reference genes have been performed in stressed Brassica napus. In this study, we investigated the expression variations of nine candidate reference genes in 40 samples of B. napus leaves subjected to various exogenous treatments. Parallel analyses by geNorm and NormFinder revealed that optimal reference genes differed across the different sets of samples. The best-ranked reference genes were PP2A and TIP41 for salt stress, TIP41 and ACT7 for heavy metal (Cr(6+)) stress, PP2A and UBC21 for drought stress, F-box and SAND for cold stress, F-box and ZNF for salicylic acid stress, TIP41, ACT7, and PP2A for methyl jasmonate stress, TIP41 and ACT7 for abscisic acid stress, and TIP41, UBC21, and PP2A for Sclerotinia sclerotiorum stress. Two newly employed reference genes, TIP41 and PP2A, showed better performances, suggesting their suitability in multiple conditions. To further validate the suitability of the reference genes, the expression patterns of BnWRKY40 and BnMKS1 were studied in parallel. This study is the first systematic analysis of reference gene selection for qRT-PCR normalization in B. napus, an agriculturally important crop, under different stress conditions. The results will contribute toward more accurate and widespread use of qRT-PCR in gene analysis of the genus Brassica.

Validation of endogenous normalizing genes for expression analyses in adult human testis and germ cell neoplasms

The measurement of gene expression levels in cells and tissues typically depends on a suitable point of reference for inferring biological relevance. For quantitative (or real-time) RT-PCR assays, the method of choice is often to normalize gene expression data to an endogenous gene that is stably expressed across the samples analysed: a so-called normalizing or housekeeping gene. Although this is a valid strategy, the identification of stable normalizing genes has proved challenging and a gene showing stable expression across all cells or tissues is unlikely to exist. Therefore, it is necessary to define suitable normalizing genes for specific cells and tissues. Here, we report on the performance of a panel of nine commonly employed normalizing genes in adult human testis and testicular pathologies. Our analyses revealed significant variability in transcript abundance for commonly used normalizers, highlighting the importance of selecting appropriate normalizing genes as comparative measurements can yield variable results when different normalizing genes are employed. Based on our results, we recommend using RPS20, RPS29 or SRSF4 when analysing relative gene expression levels in human testis and associated testicular pathologies. OCT4 and SALL4 can be used with caution as second-tier normalizers when determining changes in gene expression in germ cells and germ cell tumour components, but the relative transcript abundance appears variable between different germ cell tumour types. We further recommend that such studies should be accompanied by additional assessment of histology and cellularity of each sample.

Quantitative real time RT-PCR study of pathogen-induced gene expression in rock bream (Oplegnathus fasciatus): internal controls for data normalization

Rock bream Oplegnathus fasciatus is an important economic fish species. In this study, we evaluated the appropriateness of six housekeeping genes as internal controls for quantitative real-time PCR (RT-qPCR) analysis of gene expression in rock bream before and after pathogen infection. The expression of the selected genes in eight tissues infected with Vibrio alginolyticus or megalocytivirus was determined by RT-qPCR, and the PCR data were analyzed with geNorm and NormFinder algorithms. The results showed that before pathogen infection, mediator of RNA polymerase II transcription subunit 8 and β-actin were ranked as the most stable genes across the examined tissues. After bacterial or viral infection, the stabilities of the housekeeping genes varied to significant extents in tissue-dependent manners, and no single pair of genes was identified as suitable references for all tissues for either of the pathogen stimuli. In addition, for the majority of tissues, the most stable genes during bacterial infection differed from those during viral infection. Nevertheless, optimum reference genes were identified for each tissue under different conditions. Taken together, these results indicate that tissue type and the nature of the infectious agent used in the study can all influence the choice of normalization factors, and that the optimum reference genes identified in this study will provide a useful guidance for the selection of internal controls in future RT-PCR study of gene expression in rock bream.

Validation of reference genes aiming accurate normalization of qRT-PCR data in Dendrocalamus latiflorus Munro

BACKGROUND

Dendrocalamus latiflorus Munro distributes widely in subtropical areas and plays vital roles as valuable natural resources. The transcriptome sequencing for D. latiflorus Munro has been performed and numerous genes especially those predicted to be unique to D. latiflorus Munro were revealed. qRT-PCR has become a feasible approach to uncover gene expression profiling, and the accuracy and reliability of the results obtained depends upon the proper selection of stable reference genes for accurate normalization. Therefore, a set of suitable internal controls should be validated for D. latiflorus Munro.

RESULTS

In this report, twelve candidate reference genes were selected and the assessment of gene expression stability was performed in ten tissue samples and four leaf samples from seedlings and anther-regenerated plants of different ploidy. The PCR amplification efficiency was estimated, and the candidate genes were ranked according to their expression stability using three software packages: geNorm, NormFinder and Bestkeeper. GAPDH and EF1α were characterized to be the most stable genes among different tissues or in all the sample pools, while CYP showed low expression stability. RPL3 had the optimal performance among four leaf samples. The application of verified reference genes was illustrated by analyzing ferritin and laccase expression profiles among different experimental sets. The analysis revealed the biological variation in ferritin and laccase transcript expression among the tissues studied and the individual plants.

CONCLUSIONS

geNorm, NormFinder, and BestKeeper analyses recommended different suitable reference gene(s) for normalization according to the experimental sets. GAPDH and EF1α had the highest expression stability across different tissues and RPL3 for the other sample set. This study emphasizes the importance of validating superior reference genes for qRT-PCR analysis to accurately normalize gene expression of D. latiflorus Munro.