Candidate Reference Genes Research Articles

Unbiased RNA-Seq-driven identification and validation of reference genes for quantitative RT-PCR analyses of pooled cancer exosomes

BackgroundExosomes are extracellular vesicles (EVs) derived from endocytic compartments of eukaryotic cells which contain various biomolecules like mRNAs or miRNAs. Exosomes influence the biologic behaviour and progression of malignancies and are promising candidates as non-invasive diagnostic biomarkers or as targets for therapeutic interventions. Usually, quantitative real-time polymerase chain reaction (qRT-PCR) is used to assess gene expression in cancer exosomes, however, the ideal reference genes for normalization yet remain to be identified.ResultsIn this study, we performed an unbiased analysis of high-throughput mRNA and miRNA-sequencing data from exosomes of patients with various cancer types and identify candidate reference genes and miRNAs in cancer exosomes. The expression stability of these candidate reference genes was evaluated by the coefficient of variation “CV” and the average expression stability value “M”. We subsequently validated these candidate reference genes in exosomes from an independent cohort of ovarian cancer patients and healthy control individuals by qRT-PCR.ConclusionsOur study identifies OAZ1 and hsa-miR-6835-3p as the most reliable individual reference genes for mRNA and miRNA quantification, respectively. For superior accuracy, we recommend the use of a combination of reference genes - OAZ1/SERF2/MPP1 for mRNA and hsa-miR-6835-3p/hsa-miR-4468-3p for miRNA analyses.

Reference genes for quantitative real-time PCR normalization of Cenostigma pyramidale roots under salt stress and mycorrhizal association.

Cenostigma pyramidale is a native legume of the Brazilian semiarid region which performs symbiotic association with arbuscular mycorrhizal fungi (AMF), being an excellent model for studying genes associated with tolerance against abiotic and biotic stresses. In RT-qPCR approach, the use of reference genes is mandatory to avoid incorrect interpretation of the relative expression. This study evaluated the stability of ten candidate reference genes (CRGs) from C. pyramidale root tissues under salt stress (three collection times) and associated with AMF (three different times of salinity). The de novo transcriptome was obtained via RNA-Seq sequencing. Three algorithms were used to calculate the stability of CRGs under different conditions: (i) global (Salt, Salt+AMF, AMF and Control, and collection times), (ii) only non-inoculated plants, and (iii) AMF (only inoculated plants). HAG2, SAC1, aRP3 were the most stable CRGs for global and AMF assays, whereas HAG2, SAC1, RHS1 were the best for salt stress assay. This CRGs were used to validate the relative expression of two up-regulated transcripts in Salt2h (RAP2-3 and PIN8). Our study provides the first set of reference genes for C. pyramidale under salinity and AMF, supporting future researches on gene expression with this species.

Selection and validation of suitable reference genes in adipose tissue of Jianzhou Da'er Goat (Capra hircus)

ABSTRACT Reverse transcription quantitative real-time PCR (RT-qPCR) is a common and high-efficiency technique for detecting the mRNA levels of genes where suitable reference genes were introduced to normalize the data of RT-qPCR. However, little is known about the suitable reference genes in the adipose tissue of Jianzhou da’er goat. Here, six housekeeping genes, including the ACTB, ALAS1, EIF3K, HSP90, RPLP0, and UXT, were selected as candidate reference genes, while the FASN, HSL, LPL, and PID1 were used as the target genes for the validation of the suitability of identified reference genes. After detecting the expression levels of selected genes, the geNorm, NormFinder, and RefFinder softwares were applied to analyze the obtained data of candidate reference genes. The suitability analysis showed that the HSP90 was the most stable candidate reference gene, followed by the ALAS1 and RPLP0 gene. The validation experiment not only verified that the HSP90 was more stable than ACTB as the reference gene, but also confirmed that the HSP90 alone and a combination of HSP90, ALAS1, and RPLP0 could play the same roles in the normalization of target genes. Therefore, we strongly suggested that the HSP90 alone and the combination of three genes (HSP90, ALAS1, and RPLP0) could be used as the suitable reference for normalizing gene expression data in adipose tissue of Jianzhou da’er goat.

Reference gene screening for Real-time quantitative PCR in Polygonum multiflorum

To select suitable references gene of Polygonum multiflorum for gene expression analysis in different tissues, five candidate reference genes like Actin,GAPDH,SAND,PP2A,TIP41 were selected from the transcriptome data of P. multiflorum, then the specific primers were designed. The expression stability of the five reference genes in different tissues of P. multiflorum was analyzed by Real-time quantitative PCR through avilable analysis methods such as geNorm, NormFinder, BestKeeper, Delta CT and RefFinder, to ensure the reliability of the analysis results. The results showed that there were significant differences in the expression levels and stability of candidate genes in different tissues of P. multiflorum. Ct distribution analysis of the expression levels of candidate genes showed that the expression levels of Actin and GAPDH genes were relatively high in different tissues, while the expression levels of SAND, PP2A and TIP41 were lower. The stability of each candidate gene was analyzed by different methods. The results of geNorm analysis showed that the expression of PP2A and GAPDH was the most stable, the expression stability of SAND was the worst, the stability of PP2A was the highest in both NormFinder and Delta CT, the stability of SAND was the lowest, and the stability of Actin was the most stable in BestKeeper analysis. Through the comprehensive evaluation and analysis of the stability of candidate genes by RefFinder, it is concluded that the stability of PP2A gene is the highest, followed by GAPDH, Actin, TIP41, SAND, and SAND gene is the worst. Therefore, the PP2A gene is an ideal reference gene for the analysis of gene expression in different tissues of P. multiflorum.

Selection of reference genes for quantitative PCR analysis in poultry red mite (Dermanyssus gallinae).

Poultry red mites (PRMs, Dermanyssus gallinae) are harmful ectoparasites that affect farmed chickens and cause serious economic losses in the poultry industry worldwide. Acaricides are used for PRM control; however, some PRMs have developed acaricide-resistant properties, which have indicated the need for different approaches for PRM control. Therefore, it is necessary to elucidate the biological status of PRMs to develop alternative PRM control strategies. Quantitative polymerase chain reaction (qPCR) allows analysis of the biological status at the transcript level. However, reference genes are preferable for accurate comparison of expression level changes given the large variation in the quality of the PRM samples collected in each farm. This study aimed to identify candidate reference genes with stable expression levels in the different blood feeding states and life stages of PRMs. First, we selected candidates based on the following criteria: sufficient expression intensity and no significant expression difference between fed and starved states. We selected and characterized seven candidate reference genes. Among them, we evaluated the gene expression stability between the starved and fed states using RefFinder; moreover, we compared their expression levels in each life-stage and identified two reference genes, Elongation factor 1-alpha (ELF1A)-like and apolipophorins-like. Finally, we evaluated the utility of the candidates as reference genes, and the use of ELF1A-like and apolipophorins-like successfully normalized ATP synthase subunit g -like gene expression. Thus, ELF1A-like and apolipophorins-like could be suitable reference genes in PRMs.

RpoB and efp are stable candidate reference genes for quantitative real-time PCR analysis in Saccharopolyspora spinosa

Spinosad (spinosyn A and spinosyn D), the secondary metabolite produced by Saccharopolyspora spinosa, is a potent insecticide with low effects on the environment and mammals. Strategies such as metabolic engineering, mutagenesis and fermentation process optimization have been employed for its production enhancement. Quantitative real-time polymerase chain reaction(qRT-PCR) is one of the preferred methods for the evaluation of gene transcript levels, whose accuracy and sensitivity depend on the normalization using optimal reference genes. However, no single reference gene is universally appropriate for all strains under various conditions. In this study, the transcriptional stability of 35 candidate reference genes, including 23 traditionally used reference genes and 12 novel reference genes identified in S. spinosa homologous species, was analysed in S. spinosa ATCC 49460 and spinosad high-yield strain S. spinosa S3-3 under three fermentation phases. The transcriptional stability of these genes was assessed by three statistical algorithms, geNorm, NormFinder and BestKeeper. The overall rankings suggested that rpoB and efp were the most stable candidate reference genes, and they may be the most promising reference genes for further study on the measurement of expression levels of target genes involved in the biosynthetic process of spinosad.

Identification of quantitative polymerase chain reaction reference genes suitable for normalising gene expression in the brain of normal and dystrophic mice and dogs.

Background: In addition to progressive, debilitating muscle degeneration, ~50% of patients with Duchenne muscular dystrophy (DMD) have associated cognitive and behavioural disorders secondary to deficiency of dystrophin protein in the brain. The brain expresses a variety of dystrophin isoforms (Dp427, Dp140 and Dp71) whose functions remain to be fully elucidated. Detailed comparative analysis of gene expression in healthy and dystrophin-deficient brain is fundamental to understanding the functions of each isoform, and the consequences of their deficiency, with animal models representing a key tool in this endeavour. Reverse transcription quantitative real-time PCR (RT-qPCR) is a widely used method to study gene expression. However, accurate quantitative assessment requires normalisation of expression data using validated reference genes. The aim of this study was to identify a panel of suitable reference genes that can be used to normalise gene expression in the brain of healthy and dystrophic dogs and mice. Methods: Using the DE50-MD dog and mdx mouse models of DMD we performed RT-qPCR from fresh frozen brain tissue and employed the geNorm, BestKeeper and Normfinder algorithms to determine the stability of expression of a panel of candidate reference genes across healthy and dystrophic animals, and across different brain regions. Results: We show that SDHA, UBC and YWHAZ are suitable reference genes for normalising gene expression in healthy and dystrophic canine brain, and GAPDH, RPL13A and CYC1 in healthy and dystrophic murine brain. Notably, there was no overlap in the highest performing reference genes between the two species. Conclusions: Our findings suggest that gene expression normalisation is possible across six regions of the canine brain, and three regions of the murine brain. Our results should facilitate future work to study gene expression in the brains of normal and dystrophic dogs and mice and thus decipher the transcriptional consequences of dystrophin deficiency in the brain.

Reference gene screening for analyzing gene expression in the heart, liver, spleen, lung and kidney of forest musk deer.

The screening of reference genes for real-time quantitative PCR (qPCR) in forest musk deer (FMD) tissue is of great significance to the basic research on FMD. However, there are few reports on the stability analysis of FMD reference genes so far. In this study, We used qPCR to detect the expression levels of 11 reference gene candidates (18S rRNA, beta-actin [ACTB], glyceraldehyde-3-phosphate dehydrogenase [GAPDH], TATA box-binding protein [TBP], hypoxanthine phosphoribosyltransferase 1 [HPRT1], tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein zeta polypeptide [YWHAZ], hydroxymethylbilane synthase [HMBS], eukaryotic translation elongation factor 1 alpha 1 [EEF1A1], succinate dehydrogenase complex flavoprotein subunit A [SDHA], peptidylprolyl isomerase B [PPIB], and ubiquitin C [UBC]) in heart, liver, spleen, lung and kidney of FMD. After removing 18S rRNA on account of its high expression level, geNorm, NormFinder, BestKeeper and ΔCt algorithms were used to evaluate the expression stability of the remaining genes in the five organs, and further comprehensive ranking was calculated by RefFinder. According to the results, the selected reference genes with the most stable expression in the heart of FMD are SDHA and YWHAZ, while in the liver are ACTB and SDHA; in the spleen and lung are YWHAZ and HPRT1; in the kidney are YWHAZ and PPIB. The use of common reference genes in all five organs is not recommended. The analyses showed that tissue is an important variability factor in genes expression stability. Meanwhile, the result can be used as a reference for the selection of reference genes for qPCR in further study.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Expression Analysis of Candidate Genes in Carposina sasakii (Lepidoptera: Carposinidae)

Carposina sasakii is one of the most important pests on the quality of stone and pome fruits. Investigation of a gene expression level in the species is hampered because of the gap of validated reference genes. The expression variation in the transcription levels of eight candidate reference genes, Actin (ACT), Tubulinbeta-1 (TUB), Ribosomal protein 49 (RP49), Elongation factor1-alpha (EF-1a), Elongation factor1-b (EF-1b), Elongation factor1-d (EF-1d), Ribosomal proteinL13 (RPL13) and Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), were analyzed by quantitative real-time PCR (qPCR). The stability and ranking of these gene expression profiles in three organ types (head, thorax and abdomen), three developmental stages (larva, pupa and moth), and five diapause states (non-diapause, pre-diapause, diapause 0 d, diapause 20 d and diapause 60 d) were assessed using two algorithm-based methods, geNorm and NormFinder. EF-1a, ACT and GAPDH were evaluated to be the three stable reference genes based on the important observations and comprehensive analysis, whereas TUB and EF-1b showed low expression stability. Best gene combinations for different qPCR analysis in C. sasakii could be chosen from the three stable reference genes, the using of two reference genes is sufficient to effectively normalize qPCR data in C. sasakii. The study laid the foundation for gene expression analysis in C. sasakii and provided new information for the selection of reference genes.

Reference gene selection for miRNA and mRNA normalization in potato in response to potato virus Y

This was the first report on evaluating candidate reference genes for quantifying the expression profiles of both coding (e.g., mRNA) and non-coding (e.g., miRNA) genes in potato response to potato virus Y (PVY) inoculation. The reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) method was employed to quantify the expression profiles of eight selected candidate reference genes; their expression stability was analyzed by four statistical algorithms, i.e., geNorm, BestKeeper, NormFinder and RefFinder. The most stable reference genes were sEF1a, sTUBb and seIF5 with a high stability. The least stable ones were sPP2A, sSUI1 and sGAPDH. The same reference gene allows for normalization of both miRNA and mRNA levels from a single RNA sample using cDNAs synthesized in a single RT reaction, in which a stem-loop primer was used for miRNAs and the oligo (dT) for mRNAs.

Reference gene selections for real time quantitative PCR analysis of gene expression in different oat tissues and under salt stress

Appropriate choice of reference genes for data normalization is of critical importance for accurate real time reverse transcription quantitative PCR (RT-qPCR) analysis of gene expression. Oat is an agriculturally important crop cultivated widely around the world for grain or forage, and appropriate reference genes for reliable gene expression analysis remain to be identified. In this study, we selected nine candidate reference genes based on available oat RNA-seq data. We then conducted a systematic evaluation of the relative expression stability of these genes in different tissues and under salt stress using statistical algorithms, geNorm, NormFinder, and BestKeeper. Our findings reveal that the highest-ranked reference genes for accurate data normalization should be selected according to specific sample subsets. For different tissues, the combination of two reference genes [elongation factor 1-alpha (EF1α) and serine/threonine protein phosphatase 2] was sufficient for accurate normalization. For salt stress treatment, the combination of two reference genes (EF1α and TATA-binding protein) was sufficient for accurate normalization. Moreover, the commonly used reference genes, actin and glyceraldehyde 3-phosphate dehydrogenase, were least suitable for data normalization of oat samples. Expression of a salt stress-inducible transcription factor Avena sativa WRKY2 was investigated to validate the reference genes identified in this study. This is the first systematic study of reference gene selection in cultivated oat and provides guidelines to obtain more accurate RT-qPCR results in this species.

Validation of superior reference genes for qRT-PCR and Western blot analyses in marine Emiliania huxleyi-virus model system.

To search for a set of reference genes for reliable gene expression analysis in the globally important marine coccolithophore Emiliania huxleyi-virus model system. Fifteen housekeeping genes (CDKA, CYP15, EFG3, POLAI, RPL30, RPL13, SAMS, COX1, GPB1-2, HSP90, TUA, TUB, UBA1, CAM3 and GAPDH) were evaluated for their stability as potential reference genes for qRT-PCR using ΔCt, geNorm, NormFinder, Bestkeeper and RefFinder software. CDKA, TUA and TUB genes were tested as loading controls for Western blot in the same sample panel. Additionally, target genes associated with cell apoptosis, that is metacaspase genes, were applied to validate the selection of reference genes. The analysis results demonstrated that putative housekeeping genes exhibited significant variations in both mRNA and protein content during virus infection. After a comprehensive analysis with all the algorithms, CDKA and GAPDH were recommended as the most stable reference genes for E huxleyi virus (EhV) infection treatments. For Western blot, significant variation was seen for TUA and TUB, whereas CDKA was stably expressed, consistent with the results of qRT-PCR. CDKA and GAPDH are the best choice for gene and protein expression analysis than the other candidate reference genes under EhV infection conditions. The stable internal control genes identified in this work will help to improve the accuracy and reliability of gene expression analysis and gain insight into complex E. huxleyi-EhV interaction regulatory networks.

Transcriptome-based screening and the optimal reference genes for real-time quantitative PCR in Rehmannia chingii and R. henryi

Real time quantitative PCR (qPCR) is a powerful tool for studying the expression of specific genes. The accuracy and reliability of qPCR analysis data require the selection of reference genes with stable expression. However, the reference genes that can be used for qPCR of Rehmannia chingii and R. henryi have not yet been identified. In this study, based on the transcriptome data of R. chingii and R. henryi, we initially selected genes with relatively stable expression in different samples. We screened six candidate reference genes in R. chingii and R. henryi and calculated their expression abundance by real time qPCR. Their expression stability was evaluated by three algorithms geNorm, NormFinder, and BestKeeper. Although the results obtained by different algorithms were not completely consistent, R. chingii type 2A phosphatase activator TIP41 and R. chingii 18S ribosomal RNA had the highest expression stability in six different samples of R. chingii, and R. henryi 18S ribosomal RNA and R. henryi actin showed the most stable expression in different samples of R. henryi. In addition, based on transcriptome data, four genes were screened in R. chingii and R. henryi, and the expression stability of the selected reference genes was further verified. This study laid the foundation for further analysis and verification of the functions of important genes in R. chingii and R. henryi.

Identification of suitable mRNAs and microRNAs as reference genes for expression analyses in skin cells under sex hormone exposure.

Quantitative RT-PCR is the most accurate technique for the study of gene expression profiles, however, to ensure the accuracy of qPCR results, suitable reference genes are necessary for data normalization. Hormones influence the development and function of skin cells, regulating the expression of genes and miRNAs. Nevertheless, the stability of reference genes after sex hormone treatment has not been thoroughly investigated. In this study, we evaluated the expression of a set of candidate mRNAs and microRNsA (miRNA) as reference genes in keratinocytes (HaCaT cells), primary human fibroblasts and a melanoma cell line (LM-36 cells) under testosterone or 17β-estradiol treatment. Two algorithms, namely geNorm, Best-Keeper, and the comparative ΔCt method were used to evaluate the expression stability of the candidate reference genes. The comprehensive ranking showed that TBP and miR-191-5p are the most stable expressed genes across all cultured cells under hormone treatment. Furthermore, we observed that GAPDH, HPRT1 and U6 snRNA expression may be altered by hormone exposure, thus, these genes are not recommended as reference genes. In conclusion, the present study provides, to the best of our knowledge, the first evaluation of expressed mRNA(s) and miRNA(s) as reference genes in three different types of skin cells under the stimulation of sex hormones.

Selection and validation of reference genes for qPCR in the human dermatophyte Trichophyton rubrum exposed to different carbon sources which promote adhesion-inducing conditions.

The present study aimed to identify reference genes for qPCR analysis of T.rubrum growth in culture media which promote adhesion-inducing conditions to the host tissue. We investigated the suitability of six candidate reference genes: β-act, β-tub, ef1-α, gapdh, sdha and rpl2 in reference strain of Trichophyton rubrum in response to different environmental stimuli. The stability of these genes was determined by NormFinder, geNorm and BestKeeper software. Our data obtained from the three algorithms revealed that mRNA expression levels of two candidate reference genes, ef1-α and β-tub, remained the most stable in response to different carbon sources, while different sample sets had their own most stable reference genes, highlighting the importance of the choice of internal controls in qPCR experiments. We then checked the stability of ef1-α and β-tub reference genes expression in different T.rubrum strains, suggesting that these two genes are reliable for normalisation of qPCR. Finally, we validated the suitability of selected reference genes as internal controls for target gene (SUB3) using the 2-ΔΔCt method. The best result indicating an increase of SUB3 transcript of T.rubrum was found when the two the most stable reference (ef1-α and β-tub ) genes were used, as revealed by all three algorithms. We recommend the use of ef1-α and β-tub as reference genes for qPCR analysis of target gene expression in T.rubrum exposed to different carbon sources which promote adhesion-inducing conditions.

Selection and validation of experimental condition-specific reference genes for qRT-PCR in Metopolophium dirhodum (Walker) (Hemiptera: Aphididae)

Metopolophium dirhodum (Walker) (Hemiptera: Aphididae) is one of the most common aphid pests of winter cereals. To facilitate accurate gene expression analyses with qRT-PCR assays, the expression stability of candidate reference genes under specific experimental conditions must be verified before they can be used to normalize target gene expression levels. In this study, 10 candidate reference genes in M. dirhodum were analyzed by qRT-PCR under various experimental conditions. Their expression stability was evaluated with delta Ct, BestKeeper, geNorm, and NormFinder methods, and the final stability ranking was determined with RefFinder. The results indicate that the most appropriate sets of internal controls were SDHB and RPL8 across geographic population; RPL8, Actin, and GAPDH across developmental stage; SDHB and NADH across body part; RPL8 and Actin across wing dimorphism and temperature; RPL4 and EF1A across starvation stress; AK and RPL4 across insecticide treatments; RPL8 and NADH across antibiotic treatments; RPL8, RPL4, Actin, and NADH across all samples. The results of this study provide useful insights for establishing a standardized qRT-PCR procedure for M. dirhodum and may be relevant for identifying appropriate reference genes for molecular analyses of related insects.

Reference Gene Selection for qRT-PCR Normalization in Iris germanica L.

Quantitative real-time PCR (qPCR) is an effective and widely used method to analyze expression patterns of target genes. Selection of stable reference genes is a prerequisite for accurate normalization of target gene expression by qRT-PCR. In Iris germanica L., no studies have yet been published regarding the evaluation of potential reference genes. In this study, nine candidate reference genes were assessed at different flower developmental stages and in different tissues by four different algorithms (GeNorm, NormFinder, BestKeeper, and RefFinder). The results revealed that ACT11 (Actin 11) and EF1α (Elongation factor 1 alpha) were the most stable reference genes in different tissues, whereas TUA (Tubulin alpha) and UBC9 (Ubiquitin-protein ligase 9) were the most stable ones in different flower developmental stages. UBC9 and ACT11 were the most stable reference genes in all of the tested samples, while the SAMDC (S-Adenosylmethionine decarboxylase) showed the least stability. Finally, to validate the suitability of the selected reference genes, the relative expression level of IgTPS (beta-caryophyllene synthase) was assessed and highlighted the importance of suitable reference gene selection. This work constitutes the first systematic evaluation of potential reference genes in I. germanica and provides guidelines for future research on gene function and molecular mechanisms on I. germanica and related species.

Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp.

Quantitative real-time PCR (qRT-PCR) is commonly used to measure gene expression to further explore gene function, while suitable reference genes must be stably expressed under different experimental conditions to obtain accurate and reproducible data for relative quantification. Taxol or paclitaxel is an important anticancer compound mainly identified in Taxus spp. The molecular mechanism of the regulation of taxol biosynthesis is current research goal. However, in the case of Taxus spp., few reports were published on screening suitable reference genes as internal controls for qRT-PCR. Here, eight reference genes were selected as candidate reference genes for further study. Common statistical algorithms geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to analyze the data from samples collected from a cell line of Taxus × media under various experimental conditions and from tissues of Taxus chinensis var. mairei. The expression patterns of TcMYC under salicylic acid treatment differed significantly, with the best and worst reference genes in the cell line. This study screened out suitable reference genes (GAPDH1 and SAND) under different treatments and tissues for the accurate and reliable normalization of the qRT-PCR expression data of Taxus spp. At the same time, this study will aid future research on taxol biosynthesis-related genes expression in Taxus spp., and can also be directly used to other related species.

Identification of reference genes for real-time polymerase chain reaction gene expression studies in Nile rats fed Water-Soluble Palm Fruit Extract.

The Nile rat (Arvicanthis niloticus) is a novel diurnal carbohydrate-sensitive rodent useful for studies on type 2 diabetes mellitus (T2DM) and the metabolic syndrome. Hepatic responses to T2DM and any interventions thereof can be evaluated via transcriptomic gene expression analysis. However, the study of gene expression via real-time reverse transcription quantitative polymerase chain reaction (RT-qPCR) requires identification of stably expressed reference genes for accurate normalisation. This study describes the evaluation and identification of stable reference genes in the livers from Control Nile rats as well as those supplemented with Water-Soluble Palm Fruit Extract, which has been previously shown to attenuate T2DM in this animal model. Seven genes identified as having stable expression in RNA-Sequencing transcriptome analysis were chosen for verification using real-time RT-qPCR. Six commonly used reference genes from previous literature and two genes from a previous microarray gene expression study in Nile rats were also evaluated. The expression data of these 15 candidate reference genes were analysed using the RefFinder software which incorporated analyses performed by various algorithms. The Hpd, Pnpla6 and Vpp2 genes were identified as the most stable across the 36 samples tested. Their applicability was demonstrated through the normalisation of the gene expression profiles of two target genes, Cela1 and Lepr.In conclusion, three novel reference genes which can be used for robust normalisation of real-time RT-qPCR data were identified, thereby facilitating future hepatic gene expression studies in the Nile rat.

Screening and validation of reference genes for qPCR analysis in gonads and embryos of Takifugu bimaculatus

Suitable reference genes are one of the necessary conditions for obtaining reliable results by real-time fluorescence quantitative PCR (qPCR). In this study, the expression of the 10 common candidate reference genes (18s rRNA, rps27, cnbp, rpl7, ube2, hsp-at, gapdh, β-actin, rpl13a, 1-ef1a) at different developmental stages of gonad and embryo of Takifugu bimaculatus were analyzed by qPCR. And the expression stability of these reference genes is analyzed by GeNorm, NormFinder and Bestkeeper softwares. The results showed that the expression stability of 1-ef1a was the highest (p < 0.05) combined with the results of three analysis softwares, and the optimal number of reference genes was 3, which were 1-ef1a, hsp-at and rps27. Based on this result, we analyzed the expression of objective genes (foxl2 and dmrt1) at different developmental stages (early, medium-term and late) of gonads and embryos (eye vesicle stage, motility stage, heart anlage stage, retinal pigmentation stage, newly hatched larva) of T. bimaculatus by qPCR. Among which the highest expression level of foxl2 was in the middle stage of ovary (p < 0.05), and dmrt1 was expressed in the early stage of testis with the highest level (p < 0.05). The expression levels of foxl2 and dmrt1 in the early stages of embryonic development were low or even almost non-existent, and had a slight increase in the late embryonic development (newly hatched larva). The expression profile of objective genes are consistent with their roles in the development of gonads and embryos, which indicates that 1-ef1a, hsp-at and rps27 are the suitable reference gene combination for detecting the expression of objective genes by qPCR in gonads and embryos of T. bimaculatus. This research also provides the reliability data for screening the reference genes by qPCR in fish.