Selection of Suitable Reference Genes for Quantitative Real-time PCR in Sapium sebiferum

Metabolome and Transcriptome Reveal Chlorophyll, Carotenoid, and Anthocyanin Jointly Regulate the Color Formation of Triadica sebifera

The Chinese tallow tree (Triadica sebifera) is an economically important plant on account of its ornamental value and oil-producing seeds. Leaf colour is a key characteristic of T. sebifera, with yellow-, red- and purple-leaved varieties providing visually impressive displays during autumn. In this study, we performed metabolomic and transcriptomic analyses to gain a better understanding of the mechanisms underlying leaf colour development in purple-leaved T. sebifera at three stages during the autumnal colour transition, namely, green, hemi-purple, and purple leaves. We accordingly detected 370 flavonoid metabolites and 10 anthocyanins, among the latter of which, cyanidin-3-xyloside and peonidin-3-O-glucoside were identified as the predominant compounds in hemi-purple and purple leaves. Transcriptomic analysis revealed that structural genes associated with the anthocyanin biosynthetic pathway, chlorophyll synthesis pathway and carotenoid synthesis pathway were significantly differential expressed at the three assessed colour stages. Additionally, transcription factors associated with the MYB-bHLH-WD40 complex, including 22 R2R3-MYBs, 79 bHLHs and 44 WD40 genes, were identified as candidate regulators of the anthocyanin biosynthetic pathway. Moreover, on the basis of the identified differentially accumulated anthocyanins and key genes, we generated genetic and metabolic regulatory networks for anthocyanin biosynthesis in T. sebifera. These findings provide comprehensive information on the leaf transcriptome and three pigments of T. sebifera, thereby shedding new light on the mechanisms underlying the autumnal colouring of the leaves of this tree.

Screening of Reference Genes under Biotic Stress and Hormone Treatment of Mung Bean (Vigna radiata) by Quantitative Real-Time PCR

Mung bean (Vigna radiata) production has been greatly threatened by numerous diseases. Infection with these pathogens causes extensive changes in gene expression and the activation of hormone signal transduction. Quantitative real-time PCR (qRT-PCR) is the most common technique used for gene expression validation. Screening proper reference genes for mung bean under pathogen infection and hormone treatment is a prerequisite for ensuring the accuracy of qRT-PCR data in mung bean disease-resistance research. In this study, six candidate reference genes (Cons4, ACT, TUA, TUB, GAPDH, and EF1α) were selected to evaluate the expression stability under four soil-borne disease pathogens (Pythium myriotylum, Pythium aphanidermatum, Fusarium oxysporum, and Rhizoctonia solani) and five hormone treatments (SA, MeJA, ETH, ABA, and GA3). In the samples from different treatments, the Ct value distribution of the six candidate reference genes was different. Under the condition of hormone treatment, the Ct value ranged from a minimum of 17.87 for EF1α to a maximum of 29.63 for GAPDH. Under the condition of pathogen infection, the Ct value ranged from a minimum of 19.43 for EF1α to a maximum of 31.82 for GAPDH. After primer specificity analysis, it was found that GAPDH was not specific, so the five reference genes Cons4, ACT, TUA, TUB, and EF1α were used in subsequent experiments. The software products GeNorm, NormFinder, BestKeeper and RefFinder were used for qRT-PCR data analysis. In general, the best candidates reference genes were: TUA for SA, ABA, GA3, and Pythium myriotylum treatment; TUB for ETH treatment; ACT for MeJA and Fusarium oxysporum treatment; and EF1α for Pythium aphanidermatum and Rhizoctonia solani treatment. The most stably expressed genes in all samples were TUA, while Cons4 was the least stable reference gene. Finally, the reliability of the reference gene was further validated by analysis of the expression profiles of four mung bean genes (Vradi0146s00260, Vradi0158s00480, Vradi07g23860, and Vradi11g03350) selected from transcriptome data. Our results provide more accurate information for the normalization of qRT-PCR data in mung bean response to pathogen interaction.

Validation of stable reference genes in peripheral blood mononuclear cells for expression studies involving vector-borne haemoparasitic diseases in bovines

Normalization of gene expression data using appropriate reference genes is critical to diminish any technical bias in an experiment involving quantitative real-time PCR (qPCR). To the best of our knowledge, this is the first report offering a systematic assessment of 14 potential reference genes (RPLP0, ACTB, RPS28, YWHAZ, SDHA, PPIA, RPS9, RPS15, UXT, GAPDH, B2M, BACH1, HMBS, and PPIB) for the identification of the most stable normalizers for qPCR of target genes in peripheral blood mononuclear cells (PBMCs) of bovines for vector-borne haemoparasitic diseases such as anaplasmosis, babesiosis, theileriosis, and trypanosomiasis. A total of 38 blood samples were collected from healthy as well as diseased cattle and buffaloes representing different haemoparasitic diseases. RNA isolated from the PBMCs was subjected to qPCR for the 14 prospective internal control genes. The comprehensive ranking of the genes was accomplished by the RefFinder tool that integrates the results of three algorithms (geNorm, NormFinder, and BestKeeper) and the comparative C_T method. RPS15, B2M, and GAPDH were ranked to be the most stable genes, whereas, PPIA and HMBS emerged to be the least suitable genes. Validation of the selected reference genes by the qPCR analysis of two immunity genes, ISG15 and GPX7 was congruent with the observations of this study. We recommend that a panel of three reference genes including RPS15, B2M, and GAPDH could prove useful in delineating the transcriptional landscape of PBMCs for vector-borne haemoparasitic diseases in bovines.

I Choose You: Selecting Accurate Reference Genes for qPCR Expression Analysis in Reproductive Tissues in Arabidopsis thaliana

Quantitative real-time polymerase chain reaction (qPCR) is a widely used method to analyse the gene expression pattern in the reproductive tissues along with detecting gene levels in mutant backgrounds. This technique requires stable reference genes to normalise the expression level of target genes. Nonetheless, a considerable number of publications continue to present qPCR results normalised to a single reference gene and, to our knowledge, no comparative evaluation of multiple reference genes has been carried out in specific reproductive tissues of Arabidopsis thaliana. Herein, we assessed the expression stability levels of ten candidate reference genes (UBC9, ACT7, GAPC-2, RCE1, PP2AA3, TUA2, SAC52, YLS8, SAMDC and HIS3.3) in two conditional sets: one across flower development and the other using inflorescences from different genotypes. The stability analysis was performed using the RefFinder tool, which combines four statistical algorithms (geNorm, NormFinder, BestKeeper and the comparative ΔCt method). Our results showed that RCE1, SAC52 and TUA2 had the most stable expression in different flower developmental stages while YLS8, HIS3.3 and ACT7 were the top-ranking reference genes for normalisation in mutant studies. Furthermore, we validated our results by analysing the expression pattern of genes involved in reproduction and examining the expression of these genes in published mutant backgrounds. Overall, we provided a pool of appropriate reference genes for expression studies in reproductive tissues of A. thaliana, which will facilitate further gene expression studies in this context. More importantly, we presented a framework that will promote a consistent and accurate analysis of gene expression in any scientific field. Simultaneously, we highlighted the relevance of clearly defining and describing the experimental conditions associated with qPCR to improve scientific reproducibility.

Stepwise Optimization of Real-Time RT-PCR Analysis

Quantitative real-time reverse transcription PCR (qRT-PCR) analysis has been used routinely to quantify gene expression levels. Primer design and the optimization of qRT-PCR parameters are critical for the accuracy and reproducibility of qRT-PCR analysis. Computational tool-assisted primer design often overlooks the presence of homologous sequences of the gene of interest and the sequence similarities between homologous genes in a plant genome. This sometimes results in skipping the optimization of qRT-PCR parameters due to the false confidence in the quality of the designed primers. Here we present a stepwise optimization protocol for single nucleotide polymorphisms (SNPs)-based sequence-specific primer design and sequential optimization of primer sequences, annealing temperatures, primer concentrations, and cDNA concentration range for each reference and target gene. The goal of this optimization protocol is to achieve a standard cDNA concentration curve with an R² ≥ 0.9999 and efficiency (E) = 100 ± 5% for the best primer pair of each gene, which serves as the prerequisite for using the 2^-ΔΔCT method for data analysis.

Selection and validation of reference genes for quantitative real-time PCR of Quercus mongolica Fisch. ex Ledeb under abiotic stresses

Quercus mongolica Fisch. ex Ledeb is the main species of coniferous and broadleaved mixed forests in northeast and north China, which has high ornamental, economic, and ecological value. The appropriate reference genes must be selected for quantitative real-time PCR to reveal the molecular mechanisms of stress responses and their contribution to breeding of Q. mongolica. In the present study, we chose 11 candidate reference genes (TUA, CYP18, HIS4, RPS13, ACT97, TUB1, UBQ10, UBC5, SAND, PP2A, and SAMDC) and used four programs (GeNorm, NormFinder, BestKeeper, and RefFinder) to assess the expression stability of the above genes in roots, stems, and leaves under five abiotic stress factors (cold, salt, drought, weak light, and heavy metal). The findings revealed that under various experimental environments, the most stable genes were different; CYP18, ACT97, and RPS13 ranked the highest under most experimental environments. Moreover, two genes induced by stress, CMO and P5CS2, were chosen to demonstrate the reliability of the selected reference genes in various tissues under various stress conditions. Our research provides a significant basis for subsequent gene function studies of Q. mongolica.

Increased Diversity of Rhizosphere Bacterial Community Confers Adaptability to Coastal Environment for Sapium sebiferum Trees

Sapium sebiferum (L.) Roxb. is an economically important tree in eastern Asia, and it exhibits many traits associated with good forestation species in coastal land. However, scarce research has been conducted to elucidate the effects of rhizosphere bacterial diversity on the adaptability and viability of S. sebiferum trees grown in the coastal environment. Field trials were conducted, and rhizosphere soil samples were collected from typical coastal and forestry nursery environments. Rhizosphere bacterial communities were evaluated using 16S rRNA pyrosequencing. A total of 43 bacterial phyla were detected in all the coastal and nursery rhizospheric soil samples. Relatively higher rhizosphere community diversity was found in coastal field-grown trees. Proteobacteria, Acidobacteriota, Bacteroidota, Chloroflex, and Gemmatimonadota were dominant bacterial phyla in rhizosphere communities of tallow trees. However, the rare groups in the coastal rhizosphere soils, with a relative abundance lower than 1%, including Latescibacterota, Methylomirabilota, NB1-j, and Nitrospirota, were largely absent in the nursery field-grown tree’s rhizosphere soils. LEfSe analysis identified a total of 43 bacterial groups that were more significantly abundant in the coastal rhizosphere environment than in that of forestry nursery grown trees. Further, our cladogram analysis identified Nitrospirota, Methylomirabilota, NB1-j, and Latescibacterota as biomarkers for the coastal environment at the phylum taxonomic level. These results suggested that the adaptability of S. sebiferum trees in coastal environment might be promoted by rhizosphere microbial interactions. Complex tree–microbe interactions might enhance the resistance of the trees to coastal environment, partially by recruiting certain bacterial microbiome species, which is of high saline-alkali resistance.

A Selection of Reliable Reference Genes for Gene Expression Analysis in the Female and Male Flowers of Salix suchowensis

Salix is a dioecious plant. Research on the molecular regulation mechanism of male and female inflorescence differentiation and development is necessary to analyze sex differentiation in the willow and the underlying mechanisms of unisexual flower development. However, at present, there are no reference genes suitable for stable expression in the process of willow inflorescence development. In this study, Salix suchowensis was used as the research material, nine candidate reference genes (α-TUB1, α-TUB2, ACT, H2A, DnaJ, CDC2, GAPDH, TIP41, β-TUB) were selected, and qRT-PCR technology was used to detect the expression of each candidate reference gene in female and male flowers at different developmental stages and using five algorithms (geNorm, Normfinder, Delta Ct, BestKeeper, and RefFinder) to comprehensively evaluate the stability of candidate reference genes. The results showed that ACT and DnaJ were stably expressed in all samples and could be used as reference genes. In addition, the reliability of the screening results was further verified via an expression pattern analysis of the CFS gene that encodes flower specific transcription factor in different samples. The stable reference genes selected in this study provide the basis for future research on the expression analysis of functional genes related to the development of male and female flowers of S. suchowensis.

Reference gene selection for real-time quantitative PCR assays in different tissues of Huperzia serrata based on full-length transcriptome sequencing

Huperzia serrata (H. serrata) produces various types of effective lycopodium alkaloids, especially Huperzine A (HupA), which is a promising drug for the treatment of Alzheimer's disease. Numerous studies focused on the chemistry, bioactivities, toxicology, and clinical trials of HupA; however, the public genomic and transcriptomic resources are very limited for H. serrata research, especially for the selection of optimum reference genes. Based on the full-length transcriptome datasets and previous studies, 10 traditional and three new candidate reference genes were selected in different tissue of H. serrata. Then, two optimal reference genes GAPDHB and HisH2A were confirmed by four analysis methods. In order to further verify the accuracy of the two reference genes, they were used to analyze the expression patterns of four HupA-biosynthetic genes (lysine decarboxylas, RS-norcoclaurine 6-O-methyltransferase, cytochrome P45072A1, and copper amine oxidase). The data suggested that the expression pattern of HupA-biosynthetic genes was consistent with them in transcriptome sequencing in different tissue of H. serrata. This study identified that GAPDHB and HisH2A provides the reliable normalization for analyzing the HupA biosynthetic gene expression in different tissues of H. serrata on the transcriptional level.

An optimized protocol for stepwise optimization of real-time RT-PCR analysis

Computational tool-assisted primer design for real-time reverse transcription (RT) PCR (qPCR) analysis largely ignores the sequence similarities between sequences of homologous genes in a plant genome. It can lead to false confidence in the quality of the designed primers, which sometimes results in skipping the optimization steps for qPCR. However, the optimization of qPCR parameters plays an essential role in the efficiency, specificity, and sensitivity of each gene's primers. Here, we proposed an optimized approach to sequentially optimizing primer sequences, annealing temperatures, primer concentrations, and cDNA concentration range for each reference (and target) gene. Our approach started with a sequence-specific primer design that should be based on the single-nucleotide polymorphisms (SNPs) present in all the homologous sequences for each of the reference (and target) genes under study. By combining the efficiency calibrated and standard curve methods with the 2^-ΔΔCt method, the standard cDNA concentration curve with a logarithmic scale was obtained for each primer pair for each gene. As a result, an R² ≥ 0.9999 and the efficiency (E) = 100 ± 5% should be achieved for the best primer pair of each gene, which serve as the prerequisite for using the 2^-ΔΔCt method for data analysis. We applied our newly developed approach to identify the best reference genes in different tissues and at various inflorescence developmental stages of Tripidium ravennae, an ornamental and biomass grass, and validated their utility under varying abiotic stress conditions. We also applied this approach to test the expression stability of six reference genes in soybean under biotic stress treatment with Xanthomonas axonopodis pv. glycines (Xag). Thus, these case studies demonstrated the effectiveness of our optimized protocol for qPCR analysis.

Selection of suitable reference genes for qRT-PCR expression analysis of Codonopsis pilosula under different experimental conditions

Codonopsis pilosula is a well-known medicinal plant. Although its transcriptome sequence has been published, suitable reference genes have not been systematically identified for conducting expression analyses via quantitative real-time polymerase chain reaction (qRT-PCR). To screen appropriate genes for use with this species, we applied four different methods-GeNorm, NormFinder, BestKeeper, and RefFinder-to evaluate the stability of 13 candidates: CpiEF1Bb, CpiCACS, CpiF-Box, Cpiβ-Tubulin, CpiGAPDH, CpiActin2, CpiAPT1, CpiActin7, CpiActin8, CpiRPL6, CpiHAF1, CpiTubulin6, and CpiUBQ12. Expression was examined by qRT-PCR for various tissue types, chemical treatments, and developmental stages. For all tested samples, CpiGAPDH proved to be the most stable. Comprehensive analysis indicated that the most stable internal reference genes were CpiF-Box and CpiCACS in different tissues and at different developmental stages, respectively. Under NaCl stress, CpiAPT1 was the best internal reference gene. For methyl jasmonate and abscisic acid treatments, CpiGAPDH and CpiF-Box, respectively, presented the highest degree of expression stability. Based on these findings, we chose CpiSPL9 as the target gene for validating the suitability of these selected reference genes. All of these results provide a foundation for accurate quantification of expression levels by genes of interest in C. pilosula.

Identification and validation of superior housekeeping gene(s) for qRT-PCR data normalization in Agave sisalana (a CAM-plant) under abiotic stresses

The adaptive mechanisms in Agave species enable them to survive and exhibit remarkable tolerance to abiotic stresses. Quantitative real-time PCR is a highly reliable approach for validation of targeted differential gene expression. However, stable housekeeping gene(s) is prerequisite for accurate normalization of expression data by qRT-PCR. Till date, no systematic validation study for candidate housekeeping gene identification or evaluation has been carried-out in Agave species. A total of 17 candidate housekeeping genes were identified from the de novo assembled transcriptomic data of A. sisalana and rigorously analyzed for expression stability assessment under drought, heat, cold and NaCl stress. Different statistical algorithms like geNorm, BestKeeper, NormFinder, and RefFinder on expression data determined the superior housekeeping gene(s) for accurate normalization of the gene of interest (GOI). The comprehensive evaluation revealed the β-Tub 4, WIN-1 and CYC-A as the most stable, while EEF1α, GAPDH, and UBE2 were ranked as the least stable genes in pooled samples. Pairwise combination by geNorm showed that up to two housekeeping genes would be adequate to normalize the GOI expression data precisely. Validation of identified most and least stable housekeeping genes was carried-out by normalizing the expression data of AsHSP20 under abiotic stress conditions. Copy number of AsHSP20 gene supports the reliability of the genes used for normalization. This is the first report on the screening and validation of the housekeeping genes under abiotic stress condition in A. sisalana that would assist to understand the stress tolerance mechanisms by novel gene identification and accurate validation.

Selection and Validation of Appropriate Reference Genes for Quantitative RT-PCR Analysis in Rubia yunnanensis Diels Based on Transcriptome Data

Real-time quantitative polymerase chain reaction (RT-qPCR) has been widely applied in gene expression and transcription abundance analysis because of its high sensitivity, good repeatability, and strong specificity. Selection of relatively stable reference genes is a precondition in order to obtain the reliable analysis results. However, little is known about evaluation of a set of reference genes through scientific experiments in Rubia plants. Here, 15 candidate reference genes were selected from R. yunnanensis transcriptome database and analyzed under abiotic stresses, hormone treatments, and different tissues. Among these 15 candidate reference genes, heterogeneous nuclear ribonucleoprotein (hnRNP), TATA binding protein (TBP), ribosomal protein L5 (RPL5), malate dehydrogenase (MDH), and elongation factor 1-alpha (EF-1α) were indicated as the five most stable reference genes by four statistical programs (geNorm, NormFinder, BestKeeper, and RefFinder). Ultimately, the validity of reference genes was confirmed by normalizing the expression of o-succinylbenzoate-CoA ligase (OSBL) and isochorismate synthase (ICS) involved in the anthraquinone biosynthesis pathway in different tissues and hormone treatments. Meanwhile, four other putative genes involved in the anthraquinone biosynthesis pathway were also normalized with the selected reference genes, which showed similar expression levels with those given by transcriptome data. This work is the first research that aims at a systematic validation on the stability of reference genes selected from R. yunnanensis transcriptome data and will be conducive to analyze gene expression in R. yunnanensis or other Rubia species.

Selection of reliable reference genes for gene expression analysis in seeds at different developmental stages and across various tissues in Paeonia ostii

Paeonia ostii seeds have recently been identified as a new source of α-linolenic acid in China. Studying the gene expression patterns of unsaturated fatty acid-related genes would be helpful for understanding the mechanism of α-linolenic acid accumulation. Quantitative real-time polymerase chain reaction (qRT-PCR) is a useful method for reliably evaluating gene expression, and it is necessary to select reliable reference genes for data normalization in qRT-PCR analysis. In this study, we evaluated the expression stability of 12 candidate reference genes using four mathematical algorithms (∆C_t, BestKeeper, NormFinder, and geNorm). The web-based tool RefFinder was used to integrate the results and to provide a comprehensive ranking order. The expression stability ranking orders of reference genes were different caculated by these four algorithms, and the ranking order analyzed by the RefFinder was UBQ > Tip41 > UCE > EF-1α > α-TUB > PP2A > ACT > GAPDH > SAM > CYP > β-TUB > 18S at the different seed development stages, and UBQ > Tip41 > EF-1α > α-TUB > PP2A > UCE > GAPDH > SAM > ACT > CYP > 18S > β-TUB in P. ostii tissues. UBQ and Tip41 are the two most stable whereas 18S and β-TUB are the two least stable reference genes for gene expression in various tissues and seeds at different developmental stages in P. ostii.

Screening the Reference Genes for Quantitative Gene Expression by RT-qPCR During SE Initial Dedifferentiation in Four Gossypium hirsutum Cultivars that Have Different SE Capability

RNA sequencing (RNA-Seq)-based gene expression analysis is applicable to a wide range of biological purposes in various species. Reverse transcription quantitative PCR (RT-qPCR) is also used to assess target gene expression utilizing stably expressed reference genes as internal control under a given set of conditions. However, investigations of the reference genes for RT-qPCR normalization in the process of somatic embryogenesis (SE) initial dedifferentiation in Gossypium hirsutum are rarely reported. In this study, on the basis of our previous transcriptome data of three different induction stages during SE initial dedifferentiation process in four G. hirsutum cultivars that have different SE capability, 15 candidate genes were selected during SE initial dedifferentiation process, and their expression stability was evaluated by geNorm, NormFinder, and BestKeeper. The results indicated that the two genes of endonuclease 4 (ENDO4) and 18S ribosomal RNA (18S rRNA) showed stable expression in the four different G. hirsutum cultivars, endowing them to be appropriate reference genes during three induction stages in the four cotton cultivars. In addition, the stability and reliability of the two reference genes of ENDO4 and 18S rRNA were further verified by comparing the expressions of auxin-responsive protein 22 (AUX22) and ethylene-responsive transcription factor 17 (ERF17) between RT-qPCR results and the RNA-seq data, which showed strong positive correlation coefficient (R² = 0.8396–0.9984), validating again the steady expression of ENDO4 and 18S rRNA as the reliable reference genes. Our results provide effective reference genes for RT-qPCR normalization during SE process in different G. hirsutum cultivars.

Selection and Validation of Reference Genes for Gene Expression Studies by RT-PCR in Dalbergia odorifera

Perennial tree Dalbergia odorifera T. Chen could form the precious heartwood used to produce chinese traditional medicine, rosewood furniture and fragrances. However the formation of heartwood is time-consuming and low efficient, leading to the severe destruction of its wild resources. Thus, it is urgent to study the molecular mechanism of heartwood formation in D. odorifera. But till now, there is no report about the reference gene selection in this species. In this study, the expression stability of nine candidate reference genes were evaluated across different tissues and stems treated by wound and chemical stimulators. Four algorithms were applied to obtain the robust genes. The results support HIS2, GAPDH, and CYP to be the most stable reference genes in samples under different wound treatments while DNAj was the least stable. In different tissues, HIS2, UBQ, and RPL were the most stable reference genes while DNAj was the least stable. The selected reference genes were validated through the normalization of the qRT-PCR data of six heartwood related genes in terpene biosynthesis pathway and ethylene signal pathway. The results showed that their expression levels were accurate when they were normalized by the most stable reference gene HIS2, or by the combination of the two or three most stable reference genes. These results demonstrated that these selected reference genes are reliable.

Selection and evaluation of reference genes for qRT-PCR analysis in Euscaphis konishii Hayata based on transcriptome data

BACKGROUND

Quantitative real-time reverse transcription-polymerase chain reaction has been widely used in gene expression analysis, however, to have reliable and accurate results, reference genes are necessary to normalize gene expression under different experimental conditions. Several reliable reference genes have been reported in plants of Traditional Chinese Medicine, but none have been identified for Euscaphis konishii Hayata.

RESULTS

In this study, 12 candidate reference genes, including 3 common housekeeping genes and 9 novel genes based on E. konishii Hayata transcriptome data were selected and analyzed in different tissues (root, branch, leaf, capsule and seed), capsule and seed development stages. Expression stability was calculated using geNorm and NormFinder, the minimal number of reference genes required for accurate normalization was calculated by Vn/Vn + 1 using geNorm. EkEEF-5A-1 and EkADF2 were the two most stable reference genes for all samples, while EkGSTU1 and EkGAPDH were the most stable reference genes for tissue samples. For seed development stages, EkGAPDH and EkEEF-5A-1 were the most stable genes, whereas EkGSTU1 and EkGAPDH were identified as the two most stable genes in the capsule development stages. Two reference genes were sufficient to normalize gene expression across all sample sets.

CONCLUSION

Results of this study revealed that suitable reference genes should be selected for different experimental samples, and not all the common reference genes are suitable for different tissue samples and/or experimental conditions. In this study, we present the first data of reference genes selection for E. konishii Hayata based on transcriptome data, our data will facilitate further studies in molecular biology and gene function on E. konishii Hayata and other closely related species.

Evaluation of reference genes for reverse transcription quantitative real-time PCR (RT-qPCR) studies in Silene vulgaris considering the method of cDNA preparation

Accurate gene expression measurements are essential in studies of both crop and wild plants. Reverse transcription quantitative real-time PCR (RT-qPCR) has become a preferred tool for gene expression estimation. A selection of suitable reference genes for the normalization of transcript levels is an essential prerequisite of accurate RT-qPCR results. We evaluated the expression stability of eight candidate reference genes across roots, leaves, flower buds and pollen of Silene vulgaris (bladder campion), a model plant for the study of gynodioecy. As random priming of cDNA is recommended for the study of organellar transcripts and poly(A) selection is indicated for nuclear transcripts, we estimated gene expression with both random-primed and oligo(dT)-primed cDNA. Accordingly, we determined reference genes that perform well with oligo(dT)- and random-primed cDNA, making it possible to estimate levels of nucleus-derived transcripts in the same cDNA samples as used for organellar transcripts, a key benefit in studies of cyto-nuclear interactions. Gene expression variance was estimated by RefFinder, which integrates four different analytical tools. The SvACT and SvGAPDH genes were the most stable candidates across various organs of S. vulgaris, regardless of whether pollen was included or not.