Validation of reference genes for quantitative RT-PCR normalization in Suaeda aralocaspica, an annual halophyte with heteromorphism and C4 pathway without Kranz anatomy

Identification and selection of reference genes for analysis of gene expression by quantitative real-time PCR in the euhalophyte Suaeda altissima (L.) Pall

Сoding sequences of seven housekeeping genes: actin SaACT7, ubiquitin-conjugating protein SaUBC10, glyceraldehyde-3-phosphate dehydrogenase SaGAPDH, protein of the large subunit of ribosomes SaL2, α-tubulin SaTUA, translation elongation factor SaeEF1α, and protein phosphatase SaPP2A were identified as candidate reference genes for expression analysis of target genes in the extremely salt tolerant plant Suaeda altissima (L.) Pall. The expression profiles of the genes differed. SaACT7 and SaeEF1α demonstrated the highest expression levels, while the lowest expression levels were found for SaPP2A and SaTUA. SaPP2A and SaeEF1α genes were the most stably expressed at different steady-state salinity levels and different nitrate concentrations in nutrient solutions (NSs). SaL2, SaPP2A, and SaeEF1α genes showed the greatest stability of expression when nitrate was added to nutrient solution of plants grown under conditions of nitrate deficiency. Less constant expression was demonstrated in this experiment by SaACT7 and SaTUA. SaL2, SaACT7, SaeEF1α, and SaUBC10 genes showed the smallest expression changes under salt shock. To validate the use of the most stably expressed genes for normalization of gene expression, we checked them as reference genes to study the expression of the nitrate transporter gene SaNPF6.3 in S. altissima roots under conditions of different salinity and different nitrate supply.

Selection of Reference Genes for Expression Normalization by RT-qPCR in Dracocephalum moldavica L

Dracocephalum moldavica is widely used as an ornamental, medicine, and perfume in industry. Real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) is widely and accurately utilized for gene expression evaluations. Selecting optimal reference genes is essential for normalizing RT-qPCR results. However, the identification of suitable reference genes in D. moldavica has not been documented. A total of 12 reference genes in D. moldavica were identified by PEG6000 (15%) treatment under hypertonia conditions in different tissues (roots, stem, leaves, flower, seeds and sepal) and during three stages of flower development, then used to validate the expression stability. There were four algorithms (delta Ct, geNorm, NormFinder, and BestKeeper) used to analyze the stability. Finally, the RefFinder program was employed to evaluate the candidate reference genes' stability. The results showed that ACTIN, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and EF1α (elongation factor-1α) were stable reference genes under the PEG6000 treatment. Heat shock protein 70 (HSP70) was the most stable gene across different flower development stages. ADP-ribosylation factor (ARF) was the most stable gene in different tissues and total samples. This study provides reliable gene expression studies for future research in D. moldavica.

Identification of Reliable Reference Genes under Different Stresses and in Different Tissues of Toxicodendron succedaneum

Toxicodendron succedaneum (L.) Kuntze (T. succedaneum) is an economic tree species that produces urushiol and urushi wax, and it is of great value in industry and medicine. However, the stability of reference genes (RGs) has not been systematically reported in T. succedaneum to date. In this study, the expression of 10 candidate RGs was analyzed by RT-qPCR in different tissues (roots, stems, leaves), stress treatments (high/low temperature, drought), and hormone stimulation (jasmonic acid, JA). Then, the stability ranking of 10 candidate genes was evaluated by ∆Ct analysis and three software programs: geNorm, NormFinder, and BestKeeper. Finally, RefFinder was used to comprehensively analyze the expression stability of 10 candidate genes. The comprehensive analysis showed that TsRG05/06, TsRG01/06, and TsRG03/ACT were stable under high/low-temperature stress, drought stress, and JA treatment, respectively. TsRG03 and ACT had stable expression in different tissues. While the TsRG03 and ACT were recommended as the suitable RGs for T. succedaneum in all samples. Meanwhile, UBQ was the least suitable as a reference gene for T. succedaneum. In addition, the results of geNorm showed that the combination of two stable RGs could make the results of gene expression more accurate. These results provide alternative RGs for the study of gene function, correction, and normalization of target gene expression and directed molecular breeding in T. succedaneum.

Genome-Wide Characterization and Analysis of the bHLH Transcription Factor Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cell C4 Anatomy

Basic helix-loop-helix (bHLH) transcription factors play important roles in plant growth, development, metabolism, hormone signaling pathways, and responses to abiotic stresses. However, comprehensive genomic and functional analyses of bHLH genes have not yet been reported in desert euhalophytes. Suaeda aralocaspica, an annual C4 halophyte without Kranz anatomy, presents high photosynthetic efficiency in harsh natural habitats and is an ideal plant for identifying transcription factors involved in stress resistance. In this study, 83 bHLH genes in S. aralocaspica were identified and categorized into 21 subfamilies based on conserved motifs, gene structures, and phylogenetic analysis. Functional annotation enrichment revealed that the majority of SabHLHs were enriched in Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways involved in the response to stress conditions, as transcription factors. A number of cis-acting elements related to plant hormones and stress responses were also predicted in the promoter regions of SabHLHs, which were confirmed by expression analysis under various abiotic stress conditions (NaCl, mannitol, low temperature, ABA, GA3, MeJA, and SA); most were involved in tolerance to drought and salinity. SabHLH169 (076) protein localized in the nucleus was involved in transcriptional activity, and gene expression could be affected by different light qualities. This study is the first comprehensive analysis of the bHLH gene family in S. aralocaspica. These data will facilitate further characterization of their molecular functions in the adaptation of desert plants to abiotic stress.

The Developmental Delay of Seedlings With Cotyledons Only Confers Stress Tolerance to Suaeda aralocaspica (Chenopodiaceae) by Unique Performance on Morphology, Physiology, and Gene Expression

Cotyledons play an important role in seedling establishment, although they may just exist for a short time and become senescent upon the emergence of euphylla. So far, the detailed function of cotyledons has not been well understood. Suaeda aralocaspica is an annual halophyte distributed in cold deserts; its cotyledons could exist for a longer time, even last until maturity, and they must exert a unique function in seedling development. Therefore, in this study, we conducted a series of experiments to investigate the morphological and physiological performances of cotyledons under salt stress at different developmental stages. The results showed that the cotyledons kept growing slowly to maintain the normal physiological activities of seedlings by balancing phytohormone levels, accumulating osmoprotectants and antioxidants, and scavenging reactive oxygen species (ROS). Salt stress activated the expression of osmoprotectant-related genes and enhanced the accumulation of related primary metabolites. Furthermore, differentially expressed transcriptional profiles of the cotyledons were also analyzed by cDNA-AFLP to gain an understanding of cotyledons in response to development and salt stress, and the results revealed a progressive increase in the expression level of development-related genes, which accounted for a majority of the total tested TDFs. Meanwhile, key photosynthetic and important salt stress-related genes also actively responded. All these performances suggest that "big cotyledons" are experiencing a delayed but active developmental process, by which S. aralocaspica may survive the harsh condition of the seedling stage.

Optimum Sterilization Method for In Vitro Cultivation of Dimorphic Seeds of the Succulent Halophyte Suaeda aralocaspica

Suaeda aralocaspica is an annual halophyte in the Amaranthaceae in the saline deserts of central Asia. This plant has succulent leaves and grape-like fruits and is a potential horticultural plant. To obtain the efficient sterilization method and optimal culture conditions, two types of seeds produced from a single plant of S. aralocaspica were treated with 75% ethanol for different time durations first, and then sodium hypochlorite (NaClO) or mercury chloride (HgCl2), with five different timing treatments were used for second seed surface sterilization. Sterilized seeds were germinated on a Murashige and Skoog (MS) medium at different potential hydrogenation (pH) levels, to examine germination and seedling performance. The results showed that the highest germination percentage of brown seeds was 100% and that of black seeds was 17%. Thus, brown seeds were more suitable for further culture experiments than black seeds. For brown seeds, the sterilization effect of NaClO was better than that of HgCl2, based on the results of seed germination, contamination, and seedling survival. Rinsing with 75% ethanol for 60 s, sterilizing with NaClO for 8 min, and cultivating at pH 8.0 MS for 7 days was the best of all sterilization procedures and cultivation methods tested, which has been successfully applied to S. aralocaspica in vitro culture. The optimized protocol described here can be used as the reference for the Suaeda genus.

Genome-Wide Identification and Analysis of the Phosphoenolpyruvate Carboxylase Gene Family in Suaeda aralocaspica, an Annual Halophyte With Single-Cellular C4 Anatomy

Phosphoenolpyruvate carboxylase (PEPC) plays pivotal roles in the carbon fixation of photosynthesis and a variety of metabolic and stress pathways. Suaeda aralocaspica belongs to a single-cellular C₄ species and carries out a photosynthetic pathway in an unusually elongated chlorenchyma cell, which is expected to have PEPCs with different characteristics. To identify the different isoforms of PEPC genes in S. aralocaspica and comparatively analyze their expression and regulation patterns as well as the biochemical and enzymatic properties in this study, we characterized a bacterial-type PEPC (BTPC; SaPEPC-4) in addition to the two plant-type PEPCs (PTPCs; SaPEPC-1 and SaPEPC-2) using a genome-wide identification. SaPEPC-4 presented a lower expression level in all test combinations with an unknown function; two SaPTPCs showed distinct subcellular localizations and different spatiotemporal expression patterns but positively responded to abiotic stresses. Compared to SaPEPC-2, the expression of SaPEPC-1 specifically in chlorenchyma cell tissues was much more active with the progression of development and under various stresses, particularly sensitive to light, implying the involvement of SaPEPC-1 in a C₄ photosynthetic pathway. In contrast, SaPEPC-2 was more like a non-photosynthetic PEPC. The expression trends of two SaPTPCs in response to light, development, and abiotic stresses were also matched with the changes in PEPC activity in vivo (native) or in vitro (recombinant), and the biochemical properties of the two recombinant SaPTPCs were similar in response to various effectors while the catalytic efficiency, substrate affinity, and enzyme activity of SaPEPC-2 were higher than that of SaPEPC-1 in vitro. All the different properties between these two SaPTPCs might be involved in transcriptional (e.g., specific cis-elements), posttranscriptional [e.g., 5'-untranslated region (5'-UTR) secondary structure], or translational (e.g., PEPC phosphorylation/dephosphorylation) regulatory events. The comparative studies on the different isoforms of the PEPC gene family in S. aralocaspica may help to decipher their exact role in C₄ photosynthesis, plant growth/development, and stress resistance.

Identification and selection of reference genes for gene expression analysis by quantitative real-time PCR in Suaeda glauca's response to salinity

Quantitative real-time polymerase chain reaction (qPCR) using a stable reference gene is widely used for gene expression research. Suaeda glauca L. is a succulent halophyte and medicinal plant that is extensively used for phytoremediation and extraction of medicinal compounds. It thrives under high-salt conditions, which promote the accumulation of high-value secondary metabolites. However, a suitable reference gene has not been identified for gene expression standardization in S. glauca under saline conditions. Here, 10 candidate reference genes, ACT7, ACT11, CCD1, TUA5, UPL1, PP2A, DREB1D, V-H⁺-ATPase, MPK6, and PHT4;5, were selected from S. glauca transcriptome data. Five statistical algorithms (ΔCq, geNorm, NormFinder, BestKeeper, and RefFinder) were applied to determine the expression stabilities of these genes in 72 samples at different salt concentrations in different tissues. PP2A and TUA5 were the most stable reference genes in different tissues and salt treatments, whereas DREB1D was the least stable. The two reference genes were sufficient to normalize gene expression across all sample sets. The suitability of identified reference genes was validated with MYB and AP2 in germinating seeds of S. glauca exposed to different NaCl concentrations. Our study provides a foundational framework for standardizing qPCR analyses, enabling accurate gene expression profiling in S. glauca.

Effects of genetic and environmental factors on variations of seed heteromorphism in Suaeda aralocaspica

Seed heteromorphism is an adaptive strategy towards adversity in many halophytes. However, the underlying mechanisms and ecological significance of seed heteromorphism have not been deeply explored. Using Suaeda aralocaspica, a typical C₄ annual halophyte without Kranz anatomy, we studied seed morphology, differentiation of morphs and fruit-setting patterns, and correlated these traits with germination responses, seed characteristics and heteromorphic seed ratio. To elucidate the genetic basis of seed heteromorphism, we analysed correlated patterns of gene expression for seed development-related genes as well. We observed that S. aralocaspica produced three types of seed morph: brown, large black and small black with differences in colour, size, mass and germination behaviour; the latter two were further distinguished by their origin in female or bisexual flowers, respectively. Further analysis revealed that seed heteromorphism was associated with genetic aspects including seed positioning, seed coat differentiation and seed developmental gene expression, while variations in seed heteromorphism may be associated with environmental conditions, e.g. annual precipitation, temperature, daylight and their monthly distribution in different calendar years. Seed heteromorphism and its variations in S. aralocaspica show multilevel regulation of the bet-hedging strategy that influences phenotypic plasticity, which is a consequence of internal genetic and external environmental factor interaction. Our findings contribute to the understanding of seed heteromorphism as a potential adaptive trait of desert plant species.

Selection and Validation of Reference Genes for the qRT-PCR Assays of Populus ussuriensis Gene Expression under Abiotic Stresses and Related ABA Treatment

Populus ussuriensis Kom. is one of the most important tree species for forest renewal in the eastern mountainous areas of Northeast China due to its fast growth, high yield, and significant commercial and ecological value. The selection of optimal reference genes for the normalization of qRT-PCR data is essential for the analysis of relative gene expression. In this study, fourteen genes were selected and assessed for their expression stability during abiotic stress (drought, high salinity, and cold stress) and after the treatment with the drought-related hormone ABA. Three algorithms were used, geNorm, NormFinder, and BestKeeper, and a comprehensive ranking of candidate reference genes was produced based on their output. The most appropriate reference genes were UBQ10 and RPL24 for drought and ABA treatment, UBQ10 and TUB3 for cold stress, and UBQ10 and 60S rRNA for high salinity. Overall, UBQ10 was the most stable reference gene for use as an internal control, whereas PP2A was the least stable. The expression of two target genes (P5CS2 and GI) was used to further verify that the selected reference genes were suitable for gene expression normalization. This work comprehensively assesses the stability of reference genes in Populus ussuriensis and identifies suitable reference genes for normalization during qRT-PCR analysis.

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.

Selection of reference genes for quantitative real-time PCR analysis in halophytic plant Rhizophora apiculata

Rhizophora apiculata is a halophytic, small mangrove tree distributed along the coastal regions of the tropical and subtropical areas of the world. They are natural genetic reservoirs of salt adaptation genes and offer a unique system to explore adaptive mechanisms under salinity stress. However, there are no reliable studies available on selection and validation of reference genes for quantitative real-time polymerase chain reaction (qRT-PCR) in R. apiculata physiological tissues and in salt stress conditions. The selection of appropriate candidate reference gene for normalization of qRT-PCR data is a crucial step towards relative analysis of gene expression. In the current study, seven genes such as elongation factor 1α (EF1α), Ubiquitin (UBQ), β-tubulin (β-TUB), Actin (ACT), Ribulose1,5-bisphosphate carboxylase/oxygenase (rbcL), Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), and 18S rRNA (18S) were selected and analyzed for their expression stability. Physiological tissues such as leaf, root, stem, and flower along with salt stress leaf samples were used for selection of candidate reference genes. The high-quality expression data was obtained from biological replicates and further analyzed using five different programs such as geNorm, NormFinder, BestKeeper, Delta Ct and RefFinder. All algorithms comprehensively ranked EF1α followed by ACT as the most stable candidate reference genes in R. apiculata physiological tissues. Moreover, β-TUB and 18S were ranked as moderately stable candidate reference genes, while GAPDH and rbcL were least stable reference genes. Under salt stress, EF1α was comprehensively recommended top-ranked candidate reference gene followed by ACT and 18S. In order to validate the identified most stable candidate reference genes, EF1α, ACT, 18S and UBQ were used for relative gene expression level of sodium/proton antiporter (NHX) gene under salt stress. The expression level of NHX varied according to the internal control which showed the importance of selection of appropriate reference gene. Taken together, this is the first ever systematic attempt of selection and validation of reference gene for qRT-PCR in R. apiculata physiological tissues and in salt stress. This study would promote gene expression profiling of salt stress tolerance related genes in R. apiculata.

The m6A Reader ECT2 Controls Trichome Morphology by Affecting mRNA Stability in Arabidopsis

The epitranscriptomic mark N⁶-methyladenosine (m⁶A) can be written, read, and erased via the action of a complex network of proteins. m⁶A binding proteins read m⁶A marks and transduce their downstream regulatory effects by altering RNA metabolic processes. The characterization of m⁶A readers is an essential prerequisite for understanding the roles of m⁶A in plants, but the identities of m⁶A readers have been unclear. Here, we characterized the YTH-domain family protein ECT2 as an Arabidopsis thaliana m⁶A reader whose m⁶A binding function is required for normal trichome morphology. We developed the formaldehyde cross-linking and immunoprecipitation method to identify ECT2-RNA interaction sites at the transcriptome-wide level. This analysis demonstrated that ECT2 binding sites are strongly enriched in the 3' untranslated regions (3' UTRs) of target genes and led to the identification of a plant-specific m⁶A motif. Sequencing analysis suggested that ECT2 plays dual roles in regulating 3' UTR processing in the nucleus and facilitating mRNA stability in the cytoplasm. Disruption of ECT2 accelerated the degradation of three ECT2 binding transcripts related to trichome morphogenesis, thereby affecting trichome branching. The results shed light on the underlying mechanisms of the roles of m⁶A in RNA metabolism, as well as plant development and physiology.

Mechanism of Oxymatrine-induced Human Esophageal Cancer Cell Apoptosis by the Endoplasmic Reticulum Stress Pathway

Endoplasmic reticulum stress is one of the mechanisms of cell apoptosis. In this study, the mechanism of oxymatrine-induced human esophageal cancer Eca-109 cell apoptosis by the endoplasmic reticulum stress pathway was investigated. Eca-109 cells were cultured in vitro with different doses of oxymatrine (0.5, 1, 2 μg/mL) for 48 h. The cell viability and proliferation inhibition rate were examined by MTT assay and cell cycle assay. The apoptosis rate was examined by Annexin V-FITC/propidium iodide assay. The expression of endoplasmic reticulum stress markers, including binding immunoglobulin protein and CCAAT-enhancer-binding protein homologous protein, were determined by real-time quantitative polymerase chain reaction and western blotting, respectively. MTT data showed that oxymatrine significantly inhibited the proliferation of Eca-109 cells. The cell apoptosis rate was quantified by flow cytometry. The expression of binding immunoglobulin protein was markedly downregulated in oxymatrine-treated Eca-109 cells while that of CCAAT-enhancer-binding protein homologous protein was upregulated. Oxymatrine inhibited proliferation and induced apoptosis of human esophageal carcinoma Eca-109 cells. Thus, oxymatrine may be a potential agent for treating human esophageal cancer.

Reference gene selection for qRT-PCR assays in Stellera chamaejasme subjected to abiotic stresses and hormone treatments based on transcriptome datasets

Background

Stellera chamaejasme Linn, an important poisonous plant of the China grassland, is toxic to humans and livestock. The rapid expansion of S. chamaejasme has greatly damaged the grassland ecology and, consequently, seriously endangered the development of animal husbandry. To draft efficient prevention and control measures, it has become more urgent to carry out research on its adaptive and expansion mechanisms in different unfavorable habitats at the genetic level. Quantitative real-time polymerase chain reaction (qRT-PCR) is a widely used technique for studying gene expression at the transcript level; however, qRT-PCR requires reference genes (RGs) as endogenous controls for data normalization and only through appropriate RG selection and qRT-PCR can we guarantee the reliability and robustness of expression studies and RNA-seq data analysis. Unfortunately, little research on the selection of RGs for gene expression data normalization in S. chamaejasme has been reported.

Method

In this study, 10 candidate RGs namely, 18S, 60S, CYP, GAPCP1, GAPDH2, EF1B, MDH, SAND, TUA1, and TUA6, were singled out from the transcriptome database of S. chamaejasme, and their expression stability under three abiotic stresses (drought, cold, and salt) and three hormone treatments (abscisic acid, ABA; gibberellin, GA; ethephon, ETH) were estimated with the programs geNorm, NormFinder, and BestKeeper.

Result

Our results showed that GAPCP1 and EF1B were the best combination for the three abiotic stresses, whereas TUA6 and SAND, TUA1 and CYP, GAPDH2 and 60S were the best choices for ABA, GA, and ETH treatment, respectively. Moreover, GAPCP1 and 60S were assessed to be the best combination for all samples, and 18S was the least stable RG for use as an internal control in all of the experimental subsets. The expression patterns of two target genes (P5CS2 and GI) further verified that the RGs that we selected were suitable for gene expression normalization.

Discussion

This work is the first attempt to comprehensively estimate the stability of RGs in S. chamaejasme. Our results provide suitable RGs for high-precision normalization in qRT-PCR analysis, thereby making it more convenient to analyze gene expression under these experimental conditions.

Reference gene selection for molecular studies of dormancy in wild oat (Avena fatua L.) caryopses by RT-qPCR method

Molecular studies of primary and secondary dormancy in Avena fatua L., a serious weed of cereal and other crops, are intended to reveal the species-specific details of underlying molecular mechanisms which in turn may be useable in weed management. Among others, quantitative real-time PCR (RT-qPCR) data of comparative gene expression analysis may give some insight into the involvement of particular wild oat genes in dormancy release, maintenance or induction by unfavorable conditions. To assure obtaining biologically significant results using this method, the expression stability of selected candidate reference genes in different data subsets was evaluated using four statistical algorithms i.e. geNorm, NormFinder, Best Keeper and ΔCt method. Although some discrepancies in their ranking outputs were noticed, evidently two ubiquitin-conjugating enzyme homologs, AfUBC1 and AfUBC2, as well as one homolog of glyceraldehyde 3-phosphate dehydrogenase AfGAPDH1 and TATA-binding protein AfTBP2 appeared as more stably expressed than AfEF1a (translation elongation factor 1α), AfGAPDH2 or the least stable α-tubulin homolog AfTUA1 in caryopses and seedlings of A. fatua. Gene expression analysis of a dormancy-related wild oat transcription factor VIVIPAROUS1 (AfVP1) allowed for a validation of candidate reference genes performance. Based on the obtained results it can be recommended that the normalization factor calculated as a geometric mean of Cq values of AfUBC1, AfUBC2 and AfGAPDH1 would be optimal for RT-qPCR results normalization in the experiments comprising A. fatua caryopses of different dormancy status.

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.