Next Article in Journal
Genome-Wide Epistasis Study of Cerebrospinal Fluid Hyperphosphorylated Tau in ADNI Cohort
Previous Article in Journal
Upregulation of Heat-Shock Protein (hsp)-27 in a Patient with Heterozygous SPG11 c.1951C>T and SYNJ1 c.2614G>T Mutations Causing Clinical Spastic Paraplegia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selection of Reliable Reference Genes for Gene Expression Normalization in Sagittaria trifolia

1
College of Horticulture and Landscape, Yangzhou University, Yangzhou 225009, China
2
Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225009, China
3
Key Laboratory of Plant Functional Genomics of the Ministry of Education/Jiangsu Key Laboratory of Crop Genomics and Molecular Breeding, Yangzhou University, Yangzhou 225009, China
*
Author to whom correspondence should be addressed.
Genes 2023, 14(7), 1321; https://doi.org/10.3390/genes14071321
Submission received: 26 April 2023 / Revised: 15 June 2023 / Accepted: 21 June 2023 / Published: 23 June 2023
(This article belongs to the Special Issue Genetics and Genomics of Aquatic Vegetables)

Abstract

:
Real-time quantitative PCR (RT-qPCR) is a method with high sensitivity and convenience that has been extensively used to analyze the expression level of target genes. A reference gene with a highly stable expression is required to ensure the accuracy of experimental results. However, the report on appropriate reference genes in arrowheads (Sagittaria trifolia) is still limited. In this study, eight candidate reference genes (ACT5, UBQ, GAPDH, CYP, NAC, IDH, SLEEPER and PLA) were selected. The candidate genes were employed in a RT-qPCR assay in different tissues at different developmental stages of the same tissue (including corm, leaf and leafstalk) in arrowheads. Five statistical algorithms, GeNorm, NormFinder, BestKeeper, delta cycle threshold (ΔCt) and RefFinder, were used to evaluate the stability of these genes’ expressions in order to identify the appropriate reference genes. The results showed that UBQ was the optimum reference gene in leaf, leafstalk, root, stolon and corm, IDH exhibited the most stable expression during the expansion of corm, UBQ and PLA were the most stable reference genes in developmental stages of leaf and leafstalk, respectively. Finally, the reliability of reference genes was further confirmed by the normalization of PDS and EXP1 genes under different arrowhead tissues and developmental stages of corm, respectively. This study constitutes important guidance for the selection of reliable reference genes for analyzing the tissue- and developmental-stage-specific expression of genes in arrowheads.

1. Introduction

Real-time quantitative PCR (RT-qPCR) is an approach that uses fluorescence signals to monitor the total amount of products after each polymerase chain reaction (PCR) cycle [1]. Due to its high accuracy and strong specificity, it has been widely applied to analyze gene expression levels [2,3]. The gene expression analysis by RT-qPCR is affected by many factors, such as RNA quality, cDNA synthesis efficiency, primer specificity, and the stability of reference genes [4,5]. The expression level of the target gene can be calibrated by using stably expressed reference genes as standard, so as to reduce the difference between the target gene expression and the true value [6]. The stability of reference genes exerts a decisive effect; unstable reference genes will lead to erroneous results in gene expression quantification [7]. Selecting appropriate reference genes to standardize the expression of target genes is required to ensure the accuracy and reliability of RT-qPCR results [8,9].
Reference genes are a class of genes that is stably expressed in various tissues and cells, or stable in the same tissues and cells at all stages of plant growth. The expression of these reference genes will not be affected by the variation of external environmental factors and experimental conditions [10,11]. The commonly used reference genes include actin gene (ACT), glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH), tubulin gene (TUB), ubiquitin gene (UBQ), ribosomal RNA encoding gene (18S rRNA) and others [12,13]. Reference genes differ between species, and different experimental conditions also require different reference genes. Therefore, the study on the stability of reference gene expression in the materials from different groups is prerequisite and crucial for gene expression normalization.
Arrowhead (S. trifolia L.) is an herbaceous aquatic plant of the genus Sagittaria in the Alismaceae family. Arrowhead contains arrow-shaped leaves, white flowers, and expanded corms from an underground part. The corms show ovoid or spherical shape, and are the edible organ of arrowheads with high nutritional and medicinal values [14]. It is reported that arrowhead contains a variety of colchicine, thus showing several effects, including clearing heat, reducing swelling, preventing cancer and fighting cancer [15]. The nutrients in arrowhead corms are mainly composed of starch, protein, fat, carbohydrates and fiber [16]. Cultivated arrowheads are mainly distributed in the middle and lower reaches of the Yangtze River, South China and Southwest China, and are a kind of characteristic vegetable favored by consumers due to their abundant nutrition. They are also of great value in enriching the variety of vegetables on the table.
For arrowheads, the research on their molecular biology is lagging far behind other species. The analysis of gene expression is the most fundamental and frequently used technical tool in molecular biology research. However, there is still no consensus on which genes can be applied as reference genes in RT-qPCR experiments of arrowheads. In this study, eight potential genes (ACT5, UBQ, GAPDH, CYP, NAC, IDH, SLEEPER, and PLA) were selected as candidate genes. Then, we tested these genes expression levels in different tissues and three developmental stages of the same tissue (corm, leaf and leafstalk), and then performed the normalization assay and verification to identify stable reference genes for future RT-qPCR analysis.

2. Materials and Methods

2.1. Plant Materials

In this experiment, the S. trifolia. L. cv. ‘Baimati (BMT)’ was used as the material, which was planted in the aquatic vegetable experimental field of Yangzhou University. The arrowhead plants with 3–4 leaves were planted in the non-porous plant pot containing soil and maintained a 5 cm water layer on the surface of the soil. Early corm formation occurred at 6–7 weeks after planting. Here, the corm developmental stage was divided into three stages, including stage 1 (c1, 50 days after planting), stage 2 (c2, 70 days after planting) and stage 3 (c3, 90 days after planting). Leaf and leafstalk growth was divided into three stages, namely stage 1 (30 days after planting), stage 2 (50 days after planting) and stage 3 (70 days after planting). Different tissues (leaf, leafstalk, root, stolon, corm) were collected at stage 2 of corm development. The tissues, including corm, leaf and leafstalk, at three developmental stages were also harvested at three stages, respectively (Figure 1). Obtained samples were stored at −80 °C for subsequent experiment.

2.2. RNA Extraction and cDNA Synthesis

RNA was extracted using the Plant Total RNA Extraction Kit (Pudi, Shanghai, China). The concentration and quality of RNA was determined using a One-DropTM spectrophotometer (Thermo Fisher, Waltham, MA, USA). The OD260/280 of each RNA sample was 1.8~2.1. Purified RNA was stored in an ultra-low-temperature refrigerator at −80 °C. Then, 1000 ng RNA was used for reverse transcription to generate cDNA, and the obtained cDNA was stored at a −20 °C condition. All experimental procedures were conducted according to the manufacturer’s instructions. Three biological replicates were set for each sample.

2.3. Primer Design

The ORF (open reading frame) of eight reference genes (ACT5, UBQ, GAPDH, CYP, NAC, IDH, SLEEPER, and PLA) was retrieved from our transcriptome database and cloned from the ‘BMT’ cDNA. Specific primers of these genes were designed using Primer 6.0 software (Premier Biosoft, Palo Alto, CA, USA). The parameter for primer design was set as follows: melting temperatures (Tm) of 54–60 °C, GC content of 45~55%, primer length ranging between 18 and 25 bp, an amplicon length range of 100~250 bp. Designed primers for all target genes were synthesized by General Biology corporation (Chuzhou, Anhui, China). The sequences of all primers are listed in Table 1.

2.4. Primer’s Performance Analysis and RT-qPCR Test

The fragment of candidate reference genes was amplified using designed primers and a cDNA template from arrowhead corms. The amplification reaction system was 20 μL, including 10 μL of Taq PCR SuperMix (TransGen, Beijing, China), 1 μL of forward/reverse primers (10 μM), 2 μL of cDNA and 6 μL of ddH2O. The reaction conditions were as follows: pre-denaturation at 94 °C for 5 min, then denaturation at 94 °C for 30 s, annealing at 60 °C for 30 s, extension at 72 °C for 30 s, a total of 35 cycles, and finally 72 °C for 10 min. The amplified products were analyzed by 1% agarose gel electrophoresis to detect the uniqueness of the amplified bands. The amplification efficiency of primers was detected by a qPCR assay using a 10-fold diluted series (10×, 102×, 103×, 104× and 105×) of cDNA of arrowhead corms. The total reaction mixture for qPCR was 10 μL, including 4 μL of cDNA, 5 μL of iTap Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) and 0.5 μL of forward/reverse primer (10 μM). Reaction conditions involved pre-denaturation at 95 °C for 3 min, followed by 39 cycles at 95 °C/10 s and 60 °C/30 s. The amplicon was heated from 65 to 95 °C to obtain the dissociation curve. Three technical repetitions were set for each sample. The cDNA of all arrowhead tissues was used as the template for the analysis of the melting curve. The specificity of primer was further determined by melting curves. In addition, Cq values were used to draw standard curves, and the primer amplification efficiency (E) and correlation coefficient (R2) were analyzed. The computing formula of primer amplification efficiency was % E = (10[−1/slope] − 1) × 100%. The calculation of amplification efficiency requires that the standard curve has a satisfactory linear relationship (R2 > 0.99), and the amplification efficiency of all primers meeting the requirements of qPCR detection should be between 90% and 110%.
The expression levels of candidate reference genes in five plant tissues and three developmental stages of the same tissue (corm, leaf and leafstalk) were performed by RT-qPCR analysis, using the CFX96 PCR system (Bio-Rad, Hercules, CA, USA). Three biological replicates, each with three technical replicates, were set for the samples.

2.5. Statistical Analysis of Gene Expression Stability

In order to select the reference genes with high expression stability, four statistical algorithms were employed. Firstly, the M value (gene expression stability measurement) was computed using GeNorm to screen stable reference genes [17]. NormFinder software was used to reanalyze and calculate the S value (Normfinder stability) [18]. After evaluation by SD (standard deviation) and CV (coefficient of variation) using BestKeeper, the reference genes were ranked [19,20], and then the SD value was analyzed again by ΔCt [21]. Finally, RefFinder [22] (available online: http://blooge.cn/RefFinder/, accessed on 15 April 2023) was used to synthetically analyze data from all algorithms. The reliable reference genes in different tissues and developmental stages were screened. All algorithms were used according to the manufacturer’s instructions.

2.6. Expression Analysis of PDS and EXP1 Genes

PDS (encoding phytoene desaturase) and EXP1 (encoding one expansin) were chosen as target genes to verify the reliability of selected reference genes and their ORF sequences were listed in Figures S1 and S2. PDS was used to test the reliability of reference genes in different tissues, and EXP1 was used to verify the reference gene for corms at three developmental stages. The expression levels of two target genes, PDS and EXP1, were normalized using the top two stable and least stable reference genes. The reaction conditions of RT-qPCR consulting in the procedure are described above. Relative expression levels of PDS and EXP1 were calculated by the 2−ΔΔCt method. Three biological repetitions and technical repetitions were set for each reaction. The primers for gene expression verification of PDS and EXP1 are shown in Table 2.

2.7. Data Analysis

The significance of data was analyzed by one-way ANOVA using Duncan’s method. SPSS v20.0 software was used to conduct one-way ANOVA. Values of p < 0.05 were considered statistically significant.

3. Results

3.1. Specificity and Amplification Efficiency of Primers

PCR amplification was used to detect the specificity of primers, while electrophoresis results showed that the fragment size was correct and the band was single (Figure S3). The cDNA of all arrowhead tissues was used as the template, and melting curves of the selected reference genes were detected by qPCR. Generated melting curves showed a specific single peak, and no other miscellaneous peak was found, providing more evidence for the primer’s specificity (Figure 2). The standard curve method was employed to calculate the amplification efficiency (E) of candidate reference genes. The qPCR efficiency of eight reference genes was between 98% and 109%, and the linear correlation coefficient was above 0.99 (Table S1). The above results suggest that the specificity and amplification efficiency of the designed primers was great; the primers could be used in subsequent experiments.

3.2. Expression Levels of Candidate Reference Genes

The expression levels presented as the quantification cycle (Cq) values of eight candidate reference genes in the samples were obtained by RT-qPCR assay. The Cq values of all reference genes are shown in Figure 3 and Table S2. The Cq values of reference genes in samples from different tissues and different developmental stages of the same tissue (corm, leaf and leafstalk) ranged from 25.13 (GAPDH) to 35.26 (SLEEPER), 20.85 (UBQ) to 33.39 (NAC), 28.03 (UBQ) to 35.57 (NAC) and 26.72 (UBQ) to 34.55 (IDH), respectively. Low Cq values represented high expression levels [23]. In different arrowhead tissues, the expression of GAPDH was the highest, while SLEEPER exhibited the lowest expression (Figure 3a). The expression difference of PLA was the smallest, whereas that of GAPDH was the largest. In corms of three developmental stages, the expression levels of UBQ and GAPDH were significantly higher than those of other selected genes (Figure 3b). The expression difference of UBQ and CYP was the smallest and largest, respectively. During different developmental stages of the leaf and leafstalk, the expression of UBQ was significantly higher than that of other genes, and the expression of IDH showed significant differences among the detected samples (Figure 3c,d). However, it is not accurate to identify the stably expressed reference genes simply based on the evaluation of the original Cq value. Thus, we employed the four following statistical algorithms to screen the most stable reference genes.

3.3. GeNorm Analysis

The GeNorm algorithm was used to calculate the M value to evaluate the expression stability of reference genes [24]. The gene with the lowest M value is considered the most stable reference gene, while the gene with the highest M value is considered to be the least stable. Regrettably, this algorithm could not distinguish the M values of the top two stable genes, resulting in two most stable reference genes in the end. The M values of each gene in materials of different groups are shown in Table 3. The GeNorm algorithm required that the M value of the appropriate reference gene be less than 1.5, indicating the expression stability of each reference gene in the samples of different groups was satisfactory. UBQ and SLEEPER, with the lowest M value (0.37), were the most stable genes expressed in different tissues, followed by IDH (M value of 0.40), while GAPDH (M value of 1.46) was the least stable gene. In corms of different developmental stages, UBQ and PLA (M value of 0.19) were the most stable genes, followed by GAPDH (M value of 0.24), while CYP (M value of 0.56) was the least stable gene. Both SLEEPER and PLA (M value of 0.32 and 0.27) were the optimum reference genes in leaf and leafstalk developmental stages.

3.4. NormFinder Analysis

The NormFinder algorithm was applied to evaluate the expression stability of the eight selected reference genes, the S values of all the reference genes were calculated based on intra-group and inter-group variation [25]. The lower the S values, the more stable the gene expression. The S values of each gene in the materials from different groups are listed in Table 3. In different tissues of arrowheads, the expression of ACT5 was the most stable (S value of 0.22), followed by UBQ (S value of 0.25), and GAPDH was the least stable (S value of 2.95). In corms of different developmental stages, IDH, with the lowest S value (0.11), showed the highest expression stability, followed by ACT5 (S value of 0.13), and CYP was the worst stable reference gene with the highest S value (0.71). UBQ (S value of 0.27 and 0.26) was the optimum reference gene and IDH (S value of 1.04 and 0.74) was the more unstable reference gene in leaves and leafstalks at different developmental stages.

3.5. BestKeeper Analysis

The BestKeeper algorithm was used to calculate the standard deviation (SD) and coefficient of variation (CV) to rank the expression stability of eight reference genes [26]. The smaller the SD and CV values, the more stable the gene expression. The SD and CV values of each gene in the materials from different groups are shown in Table 3. BestKeeper algorithm believes that the SD value of the appropriate reference gene should be less than 1, so GAPDH and IDH were excluded from the analysis of different tissues and different developmental stages of leaf and leafstalk in arrowheads, respectively. In different tissues, NAC was the most stable gene (SD value of 0.72, CV value of 2.31), followed by UBQ (SD value of 0.76, CV value of 2.73), and GAPDH was the least stable gene (SD value of 2.99, CV value of 10.24). In corms of different developmental stages, SLEEPER (SD value of 0.29, CV value of 0.94) and UBQ (SD value of 0.30, CV value of 1.40) were judged to be the top two genes in terms of stable expression, and NAC (SD value of 0.67, CV value of 2.06) showed the least stable expression. During the different developmental stages of leaves and leafstalks, CYP (SD value of 0.21 and 0.33, CV value of 0.68 and 1.05) expression was the most stable, while IDH (SD value of 1.00 and 1.08, CV value of 3.03 and 3.24) was, to the contrary, not stable.

3.6. ΔCt Analysis

The expression stability of the selected reference genes was also evaluated using the ΔCt algorithm based on the SD. The gene with the lowest S value is thought to be the most stable [27]. As shown in Table 3, the S values of each gene in the arrowhead samples from different groups were assessed. In different tissues, UBQ showed the lowest S value (1.04) and it was the most stable reference gene, whereas GAPDH (S value of 3.03) was the least stable. The expression of IDH was the most stable (S value of 0.44) in arrowhead corms among the three developmental stages, followed by ACT5 (S value of 0.46), and CYP exhibited the least stable (S value of 0.78) expression. In leaves and leafstalks at different developmental stages, UBQ (S value of 0.69) and PLA (S value of 0.51) showed the most expression stability, while the expression of IDH (S value of 1.15 and 0.81) was the most unstable.

3.7. RefFinder Analysis

RefFinder is a comprehensive web-based tool that includes different algorithms from GeNorm, NormFinder, BestKeeper and ΔCt [22]. Based on the results of all statistic methods, the stability of the arrowhead candidate reference genes was comprehensively evaluated by RefFinder. The comprehensive analysis results are recorded in Table 3.
After comprehensive comparative analysis based on the above four algorithms, we found that UBQ was the most stable gene expressed in different arrowhead tissues, followed by SLEEPER, whereas GAPDH was the least stable. In different developmental stages of corms, IDH showed the most stable expression, followed by UBQ, and CYP obtained the lowest score in expression stability. The stability scores of UBQ and PLA were the highest during the developmental stages of leaves and leafstalks, respectively, whereas the stability score of IDH was the lowest.

3.8. Reference Gene Validation

Phytoene desaturase (PDS) is a key enzyme involved in carotenoid synthesis, which also has been proved to function in regulating chlorophyll deposition in plant tissues [28]. PDS is usually highly expressed in plant tissues that appear green. EXP1 (encoding one expansin) is the major factor regulating cell wall extension and plays an important role in plant organ expansion [29]. In order to verify the reliability and accuracy of the selected reference genes, we measured the relative expression of PDS in different tissues and EXP1 in corms from different developmental stages. The results show that the relative expression levels of PDS in green leaf and leafstalk are significantly higher than those in three other tissues when UBQ and SLEEPER were used as the reference genes. When the most unstable reference gene GAPDH was used as the internal standard, PDS was only highly expressed in the leaf (Figure 4a–c). During the corm expansion process, stably expressed IDH and UBQ were employed as the reference gene, and the relative expression level of EXP1 was up-regulated with the development of corms and peaked in stage 3 of corm development. When CYP was used as the reference gene, the EXP1 expression pattern was overestimated at stage 2 of corm development. (Figure 4d–f). The above results demonstrate that our screening for reference genes in varied arrowhead tissues and corms of different development stages was reliable.

4. Discussion

In recent years, RT-qPCR has become the preferred method in detecting gene expression [30]. Normalization of data by a reference gene with a highly stable expression is a prerequisite in ensuring the reliability of experimental results [31]. In this study, the expression levels of eight candidate reference genes (ACT5, UBQ, GAPDH, CYP, NAC, IDH, SLEEPER and PLA) were determined by using the samples from different tissues and different developmental stages of the same tissue (corm, leaf, leafstalk) in arrowheads. This provided a reference for studying the physiological and biochemical processes of arrowheads at the gene expression level.
Four statistical algorithms, GeNorm, NormFinder, BestKeeper and ΔCt, were used to evaluate the stability of these selected reference genes [32]. The most stable reference genes screened from several statistical algorithms were different, while the least stable reference gene was highly consistent. In different arrowhead tissues, the combination of UBQ and SLEEPER, ACT5, NAC and UBQ ranked first in GeNorm, NormFinder, BestKeeper and ∆Ct algorithms, respectively, and were considered to be the most stable reference genes. Whereas GAPDH was judged to be the least stable reference gene throughout the analysis of four algorithms. During the expansion of corms, both Normfinder and ∆Ct algorithms considered IDH to be the best reference gene, Genorm and BestKeeper algorithms indicated that the combination of UBQI, PLA and SLEEPER may be the most stable, respectively. CYP and NAC constituted the least stable reference genes based on the analysis of four algorithms. In three developmental stages of leaves and leafstalks, UBQ was considered as the most stably expressed gene by the NormFinder algorithm and ΔCt algorithm, and the expression of IDH showed the least stability. Since different statistical algorithms were applied, these differences in screening of optimal internal reference genes were acceptable. Finally, the RefFinder algorithm was used to comprehensively compare and rank these reference genes [33], suggesting that UBQ was the best reference gene in different tissues and growth of leaves, and IDH and PLA were the optimum references for gene expression normalization during the growth of corms and leafstalks, respectively. Based on the results that mentioned above, it is concluded that UBQ is the most suitable reference gene in arrowhead tissues.
UBQ and IDH have also been reported as reference genes in other plants. UBQ was verified to be a stably expressed reference gene in different parts and seeds of different developmental stages of Paeonia ostii [34], while a similar report has been found in Dimocarpus longan [35]. IDH has been used as the reference gene under specific conditions in other species, such as different developmental stages of flowers in Pinus massoniana, and various stress conditions in Eucalyptus [36,37]. The reliable reference genes evaluated by GeNorm, NormFinder and ∆Ct in our work were similar to previous findings.
In order to verify the accuracy of the selected optimal reference genes, we detected the expression levels of PDS and EXP1 in different tissues and corms of different developmental stages using different reference genes. PDS was closely linked to the accumulation of chlorophyll and was highly expressed in green plant tissues in general [38]. In Triticum aestivum, the expression of the PDS gene was highest in the leaf and petal, followed by the stem and inflorescence, and lowest in the decolorized root [39]. In Andrographis paniculata, PDS exhibited the highest transcription level in the leaf [38]. In Allium sativum, the PDS gene was highly expressed in green leaves compared to those in roots [40]. Considering our results in arrowheads, when UBQ and SLEEPER genes were applied to normalize gene expression, the PDS gene showed a significantly high expression in leaves and leafstalks, and a low expression in roots, stolons and corms. This result is consistent with the findings in previous research. However, the expression pattern of PDS normalized by the GAPDH gene presented an inconsistent result. EXP1 encoded a class of expansins, which contributed to cell expansion during the development of plant tissue [41]. In Psidium guajava, the expression of EXP1 increased with fruit ripening and reached the highest level at the ripening stage [42]. In Solanum tuberosum, EXP was closely related to the tuber expansion process, and the EXP gene exhibited high expression in growing tubers [43]. Under the normalization of IDH and UBQ, the expression level of EXP1 during corms development was consistent with that previously described. However, the expression level of EXP1 was overestimated at stage 2 of corm development when the unstable reference gene (CYP) was used to normalize gene expression. The assay of PDS and EXP1 gene expression further confirmed the accuracy of the selected reference genes in arrowheads.

5. Conclusions

In summary, five statistical algorithms were used to compare and screen suitable reference genes in the samples from different tissues and three developmental stages of the same tissue (corm, leaf and leafstalk) in arrowheads. The results concluded that UBQ was the most stable reference gene in different tissues and developmental stages of leaves, while IDH and PLA were the most stable reference genes during developmental stages of corms and leafstalks, respectively. The accuracy of the selected reference genes was verified by the expression normalization of the PDS gene in different tissues and the EXP1 gene in corms at three developmental stages. The current work will provide a reliable basis for the future research about the tissue- and developmental-stage-specific expression analysis of genes in arrowheads.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/genes14071321/s1. Figure S1: The ORF sequence of PDS gene in ‘BMT’ arrowhead. Figure S2: The ORF sequence of EXP1 gene in ‘BMT’ arrowhead. Figure S3: PCR amplification of specific primers for selected candidate reference genes. Table S1: Amplification efficiency of candidate reference genes primer pairs. Table S2: The quantification cycle (Cq) values of the candidate reference genes.

Author Contributions

Conceptualization, J.T. and M.M.; data curation, J.T., E.L. and J.L.; formal analysis, J.T., E.L. and J.L.; investigation, J.T., E.L. and J.L.; methodology, J.T., E.L. and J.L.; resources, Z.Z., B.H. and J.J.; supervision, Z.Z., B.H. and J.J.; writing—review and editing, J.T. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Jiangsu Modern Agricultural (Vegetable) Industrial Technology System (JATS[2022]495).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brattelid, T.; Levy, F.O. Quantification of GPCR mRNA using real-time RT-PCR. Methods Mol. Biol. 2011, 746, 165–193. [Google Scholar]
  2. Niu, L.; Tao, Y.; Chen, M.; Fu, Q.; Li, C.; Dong, Y.; Wang, X.; He, H.; Xu, Z. Selection of reliable reference genes for gene expression studies of a promising oilseed crop, Plukenetia volubilis, by real-time quantitative PCR. Int. J. Mol. Sci. 2015, 16, 12513–12530. [Google Scholar] [CrossRef] [Green Version]
  3. Sun, Z.; Li, S.; Sun, M. Selection of reliable reference genes for gene expression studies in Clonostachys rosea 67-1 under sclerotial induction. J. Microbiol. Methods 2015, 114, 62–65. [Google Scholar] [CrossRef] [PubMed]
  4. Gao, D.; Kong, F.; Sun, P.; Bi, G.; Mao, Y. Transcriptome-wide identification of optimal reference genes for expression analysis of Pyropia yezoensis responses to abiotic stress. BMC Genomics 2018, 19, 251. [Google Scholar] [CrossRef] [Green Version]
  5. Sinha, P.; Singh, V.K.; Suryanarayana, V.; Krishnamurthy, L.; Saxena, R.K.; Varshney, R.K. Evaluation and validation of housekeeping genes as reference for gene expression studies in pigeonpea (Cajanus cajan) under drought stress conditions. PLoS ONE 2015, 10, e0122847. [Google Scholar] [CrossRef] [PubMed]
  6. Dheda, K.; Huggett, J.F.; Chang, J.S.; Kim, L.U.; Bustin, S.A.; Johnson, M.A.; Rook, G.A.W.; Zumla, A. The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal. Biochem. 2005, 344, 141–143. [Google Scholar] [CrossRef]
  7. Wang, M.; Wang, Q.; Zhang, B. Evaluation and selection of reliable reference genes for gene expression under abiotic stress in cotton (Gossypium hirsutum L.). Gene 2013, 530, 44–50. [Google Scholar] [CrossRef]
  8. Andersen, C.L.; Jensen, J.L.; Orntoft, T.F. Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef] [Green Version]
  9. Zhu, X.; Wang, B.; Wang, X.; Wei, X. Screening of stable internal reference gene of Quinoa under hormone treatment and abiotic stress. Physiol. Mol. Biol. Plants 2021, 27, 2459–2470. [Google Scholar] [CrossRef] [PubMed]
  10. Popovici, V.; Goldstein, D.R.; Antonov, J.; Jaggi, R.; Delorenzi, M.; Wirapati, P. Selecting control genes for RT-QPCR using public microarray data. BMC. Bioinformatics 2009, 10, 42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wang, G.; Tian, C.; Wang, Y.; Wan, F.; Hu, L.; Xiong, A.; Tian, J. Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress. PeerJ 2019, 7, e7319. [Google Scholar] [CrossRef] [PubMed]
  12. Dheda, K.; Huggett, J.F.; Bustin, S.A.; Johnson, M.A.; Rook, G.; Zumla, A. Validation of housekeeping genes for normalizing RNA expression in real-time PCR. Biotechniques 2004, 37, 112–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gambarotta, G.; Ronchi, G.; Friard, O.; Galletta, P.; Perroteau, I.; Geuna, S. Identification and validation of suitable housekeeping genes for normalizing quantitative real-time PCR assays in injured peripheral nerves. PLoS ONE 2014, 9, e105601. [Google Scholar] [CrossRef] [PubMed]
  14. Itoh, K. Ecological studies on tubers of arrowhead (Sagittaria trifolia L.), a perennial weed in the paddy fields. Bull. Natl. Agric. Res. Cent. 1997, 34, 15–89. [Google Scholar]
  15. Zhang, Y.; Yang, G.; Wang, X.; Ni, G.; Cui, Z.; Yan, Z. Sagittaria trifolia tuber: Bioconstituents, processing, products, and health benefits. J. Sci. Food Agric. 2021, 101, 3085–3098. [Google Scholar] [CrossRef] [PubMed]
  16. Devi, R.; Kumar, S. Sagittaria trifolia L.: A potential nutraceutical of the Northeastern part of India. Biodivers. Conserv. 2019, 3, S1. [Google Scholar]
  17. Vandesompele, J.; De Preter, K.; Pattyn, F.; Poppe, B.; Van Roy, N.; De Paepe, A.; Speleman, F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3, research0034. [Google Scholar] [CrossRef] [Green Version]
  18. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. Normalization of real-time quantitative reverse transcription-PCR data: A mod-el-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar]
  19. Lin, L.; Han, X.; Chen, Y.; Wu, Q.; Wang, Y. Identification of appropriate reference genes for normalizing transcript expression by quantitative real-time PCR in Litsea cubeba. Mol. Genet. Genomics 2013, 288, 727–737. [Google Scholar] [CrossRef]
  20. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations. Biotechnol. Lett. 2004, 26, 509–515. [Google Scholar] [CrossRef] [PubMed]
  21. Silver, N.; Best, S.; Jiang, J.; Thein, S.L. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol. Biol. 2006, 7, 33. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Xie, F.; Wang, J.; Zhang, B. RefFinder: A web-based tool for comprehensively analyzing and identifying reference genes. Funct. Integr. Genomics 2023, 23, 125. [Google Scholar]
  23. Yang, Z.; Zhang, R.; Zhou, Z. Identification and validation of reference genes for gene expression analysis in Schima superba. Genes 2021, 12, 732. [Google Scholar] [CrossRef]
  24. Gantasala, N.P.; Papolu, P.K.; Thakur, P.K.; Kamaraju, D.; Sreevathsa, R.; Rao, U. Selection and validation of reference genes for quantitative gene expression studies by real-time PCR in eggplant (Solanum melongena L.). BMC Res. Notes 2013, 6, 312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, M.; Bhullar, N.K. Selection of suitable reference genes for qRT-PCR gene expression studies in rice. Methods Mol. Biol. 2021, 2238, 293–312. [Google Scholar] [PubMed]
  26. Brunner, A.M.; Yakovlev, I.A.; Strauss, S.H. Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol. 2004, 4, 14. [Google Scholar] [CrossRef] [Green Version]
  27. Aggarwal, A.; Jamwal, M.; Viswanathan, G.K.; Sharma, P.; Sachdeva, M.S.; Bansal, D.; Malhotra, P.; Das, R. Optimal reference gene selection for expression studies in human reticulocytes. J. Mol. Diagn. 2018, 20, 326–333. [Google Scholar] [CrossRef] [Green Version]
  28. Zhu, Y.; Jiang, J.; Yan, Y.; Chen, X. Isolation and characterization of phytoene desaturase cDNA involved in the beta-carotene biosynthetic pathway in Dunaliella salina. J. Agric. Food Chem. 2005, 53, 5593–5597. [Google Scholar] [CrossRef]
  29. Wu, Y.; Meeley, R.B.; Cosgrove, D.J. Analysis and expression of the alpha-expansin and beta-expansin gene families in maize. Plant Physiol. 2001, 126, 222–232. [Google Scholar] [CrossRef] [Green Version]
  30. Li, J.; Han, X.; Wang, C.; Qi, W.; Zhang, W.; Tang, L.; Zhao, X. Validation of suitable reference genes for RT-qPCR data in Achyranthes bidentata Blume under different experimental conditions. Front. Plant Sci. 2017, 8, 776. [Google Scholar] [CrossRef] [Green Version]
  31. Bu, J.; Zhao, J.; Liu, M. Expression stabilities of candidate reference genes for RT-qPCR in Chinese jujube (Ziziphus jujuba Mill.) under a variety of conditions. PLoS ONE 2016, 11, e0154212. [Google Scholar] [CrossRef] [Green Version]
  32. Yuan, Z.; Zhang, X.; Pang, Y.; Qi, Y.; Wang, Q.; Ren, S.; Hu, Y.; Zhao, Y.; Wang, T.; Huo, L. Screening of stably expressed internal reference genes for quantitative real-time PCR analysis in quail. Biol. Bull. 2022, 49, 418–427. [Google Scholar] [CrossRef]
  33. Chen, R.; Chen, W.; Tigabu, M.; Zhong, W.; Li, Y.; Ma, X.; Li, M. Screening and evaluation of stable reference genes for quantitative real-time polymerase chain reaction (qRT-PCR) analysis in Chinese fir roots under water, phosphorus, and nitrogen stresses. Forests 2019, 10, 1087. [Google Scholar] [CrossRef] [Green Version]
  34. Li, C.; Hu, L.; Wang, X.; Liu, H.; Tian, H.; Wang, J. Selection of reliable reference genes for gene expression analysis in seeds at different developmental stages and across various tissues in Paeonia ostii. Mol. Biol. Rep. 2019, 46, 6003–6011. [Google Scholar] [CrossRef] [PubMed]
  35. Lin, Y.; Lai, Z. Reference gene selection for qPCR analysis during somatic embryogenesis in longan tree. Plant Sci. 2010, 178, 359–365. [Google Scholar] [CrossRef]
  36. Moura, J.C.M.S.; Araujo, P.; Brito, M.d.S.; Souza, U.R.; Viana, J.d.O.F.; Mazzafera, P. Validation of reference genes from Eucalyptus spp. under different stress conditions. BMC Res. Notes 2012, 5, 634. [Google Scholar] [CrossRef] [Green Version]
  37. Chen, H.; Yang, Z.; Hu, Y.; Tan, J.; Jia, J.; Xu, H.; Chen, X. Reference genes selection for quantitative gene expression studies in Pinus massoniana L. Trees-Struct. Funct. 2016, 30, 685–696. [Google Scholar] [CrossRef]
  38. Choi, S.; Seo, Y.B.; Han Kyu, L.; Soo-Wan, N.; Gun-Do, K. Molecular cloning and overexpression of phytoene desaturase (CrtI) from Paracoccus haeundaensis. Microbiol. Biotechnol. Lett. 2018, 46, 145–153. [Google Scholar] [CrossRef] [Green Version]
  39. Cong, L.; Wang, C.; Li, Z.; Chen, L.; Yang, G.; Wang, Y.; He, G. cDNA cloning and expression analysis of wheat (Triticum aestivum L.) phytoene and zeta-carotene desaturase genes. Mol. Biol. Rep. 2010, 37, 3351–3361. [Google Scholar] [CrossRef]
  40. Tuan, P.A.; Kim, J.K.; Kim, H.H.; Lee, S.Y.; Park, N.I.; Park, S.U. Carotenoid accumulation and characterization of cDNAs encoding phytoene synthase and phytoene desaturase in Garlic (Allium sativum). J. Agric. Food Chem. 2011, 59, 5412–5417. [Google Scholar] [CrossRef]
  41. Zhou, Y.; Luo, S.; Hameed, S.; Xiao, D.; Zhan, J.; Wang, A.; He, L. Integrated mRNA and miRNA transcriptome analysis reveals a regulatory network for tuber expansion in Chinese yam (Dioscorea opposita). BMC Genomics 2020, 21, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Silva, A.I.R.; Palenius, H.G.N.; Guzman, G.H.; Solis, A.G.A.; Pina, C.G.; Dominguez, J.F.M. Ripening-related cDNAs in guava fruit (Psidium guajava L.). Characterization and expression analysis. Rev. Fitotec. Mex. 2013, 36, 117–125. [Google Scholar]
  43. Mao, Y. Research progress of tomato and potato extension protein. J. Heilongjiang Bayi Agric. Univ. 2013, 25, 5–8. [Google Scholar]
Figure 1. The definitions of tissue and corm developmental stages in arrowheads. (a) Five tissues of arrowheads. l: leaf, ls: leafstalk, st: stolon, r: root, c: corm. (b) Three stages of arrowheads corm development (c1–c3).
Figure 1. The definitions of tissue and corm developmental stages in arrowheads. (a) Five tissues of arrowheads. l: leaf, ls: leafstalk, st: stolon, r: root, c: corm. (b) Three stages of arrowheads corm development (c1–c3).
Genes 14 01321 g001
Figure 2. Melting curves of candidate reference genes.
Figure 2. Melting curves of candidate reference genes.
Genes 14 01321 g002
Figure 3. The quantification cycle (Cq) values of the candidate reference genes. (a) Cq of eight candidate references in different tissues. (b) Cq of eight candidate references in corms at different developmental stages. (c) Cq of eight candidate references in leaves at different developmental stages. (d) Cq of eight candidate references in leafstalks at different developmental stages. A line across the box-plot graph represents the median. Lower and upper boxes indicate the 25th percentile and 75th percentiles, respectively. The whisker caps represent the maximum and minimum values.
Figure 3. The quantification cycle (Cq) values of the candidate reference genes. (a) Cq of eight candidate references in different tissues. (b) Cq of eight candidate references in corms at different developmental stages. (c) Cq of eight candidate references in leaves at different developmental stages. (d) Cq of eight candidate references in leafstalks at different developmental stages. A line across the box-plot graph represents the median. Lower and upper boxes indicate the 25th percentile and 75th percentiles, respectively. The whisker caps represent the maximum and minimum values.
Genes 14 01321 g003
Figure 4. The relative expression levels of PDS and EXP1 gene normalized by different reference genes. (ac) Relative expression level of PDS in different arrowhead tissues. (df) Relative expression level of EXP1 in corms of different developmental stages. l: leaf, ls: leafstalk, st: stolon, r: root, c: corm, c1–c3: the three stages of arrowhead corm development. Error bars represent the standard deviation of three replicates. Different lowercase letters indicate significant differences at p  <  0.05.
Figure 4. The relative expression levels of PDS and EXP1 gene normalized by different reference genes. (ac) Relative expression level of PDS in different arrowhead tissues. (df) Relative expression level of EXP1 in corms of different developmental stages. l: leaf, ls: leafstalk, st: stolon, r: root, c: corm, c1–c3: the three stages of arrowhead corm development. Error bars represent the standard deviation of three replicates. Different lowercase letters indicate significant differences at p  <  0.05.
Genes 14 01321 g004
Table 1. Candidate reference genes and primer sequences.
Table 1. Candidate reference genes and primer sequences.
Gene SymbolGene NamePrimer Sequences (5′–3′)
ACT5Actin 5Forward
Reverse
CGGAGAGGCTGTCACGATT
CGAACTGATAGACGACATTGGAA
UBQUbiquitinForward
Reverse
AGGTGCTACTCTCCATCTGTTC
CATACTTGTTCCTGTCTGTCTTGT
GAPDHGlyceraldehyde-3-phospho
dehydrogense
Forward
Reverse
CTCCTCCTCGCAATACTCACA
CAGCCACAACTTCAACATCGT
CYPCyclophilinForward
Reverse
CTCCACCTTCCACCGTATCAT
GAACTTCAAGCCGTAGATAGACTC
NACNAC domain proteinForward
Reverse
ATTCCTCCAACTGTCTGATGAAC
ACTGCCTGATATGAGCCTGTT
IDHIsocitrate dehydrogenaseForward
Reverse
CTGAGAACAGCAGCGGTAAC
GGCAACGGTCGCATAAGAG
SLEEPERZinc finger BED domain-containing protein
DAYSLEEPER-like
Forward
Reverse
CCTGTATTGCTTGGAAGTGGTAT
CCGTCGTCTTGATTGAATGGAT
PLAPhospholipase A IForward
Reverse
GACCTTATCTGTGGCACTTCTAC
CTCTCCAAGTTGCTGCTTCATT
Table 2. Primer sequences of PDS and EXP1 used for RT-qPCR.
Table 2. Primer sequences of PDS and EXP1 used for RT-qPCR.
Gene SymbolGene NamePrimer Sequences (5′–3′)
PDSPhytoene desaturaseForward
Reverse
CAGTCTCCGTCTCCACCGTCAA
CCACCATCCGCAAGATACTTAGCC
EXP1ExpansinForward
Reverse
GTCAACATTGTGGTCACGGATTACG
TGTGCCATCAAGCCATCAGTCAG
Table 3. Expression stability analysis of candidate reference genes.
Table 3. Expression stability analysis of candidate reference genes.
MaterialsRankGeNormNormFinderBestKeeperDelta CTRefFinder
GeneStabilityGeneStabilityGeneSDCVGeneStabilityGeneValues
Different tissues1UBQ0.37ACT50.22NAC0.722.31UBQ1.04UBQ1.41
2SLEEPER0.37UBQ0.25UBQ0.762.73IDH1.07SLEEPER2.78
3IDH0.40IDH0.28SLEEPER0.852.62ACT51.08ACT53.03
4ACT50.46PLA0.39PLA0.862.63SLEEPER1.10IDH3.22
5PLA0.51SLEEPER0.51CYP0.943.23PLA1.12NAC3.83
6NAC0.70NAC1.00IDH0.962.99NAC1.45PLA4.47
7CYP0.94CYP1.52ACT50.983.00CYP1.81CYP6.44
8GAPDH1.46GAPDH2.95GAPDH2.9910.24GAPDH3.03GAPDH8.00
Corm developmental stages1UBQ0.19IDH0.11SLEEPER0.290.94IDH0.44IDH2.11
2PLA0.19ACT50.13UBQ0.301.40ACT50.46UBQ2.21
3GAPDH0.24SLEEPER0.31PLA0.321.05UBQ0.51PLA2.78
4IDH0.31UBQ0.37GAPDH0.411.85PLA0.52SLEEPER3.08
5ACT50.33PLA0.40IDH0.471.62SLEEPER0.54ACT53.31
6SLEEPER0.38GAPDH0.44ACT50.471.58GAPDH0.54GAPDH4.56
7NAC0.49NAC0.64CYP0.672.21NAC0.73NAC7.24
8CYP0.56CYP0.71NAC0.672.06CYP0.78CYP7.74
Leaf developmental stages1SLEEPER0.32UBQ0.27CYP0.210.68UBQ0.69UBQ1.68
2PLA0.32ACT50.33UBQ0.260.91ACT50.71ACT52.78
3ACT50.49GAPDH0.44GAPDH0.441.29SLEEPER0.76SLEEPER3.03
4UBQ0.59SLEEPER0.45NAC0.511.48GAPDH0.76CYP3.50
5GAPDH0.63CYP0.52ACT50.551.64CYP0.79GAPDH3.66
6CYP0.65PLA0.57PLA0.641.92PLA0.81PLA3.83
7NAC0.74NAC0.95SLEEPER0.651.96NAC1.06NAC6.09
8IDH0.84IDH1.04IDH1.003.03IDH1.15IDH8.00
Leafstalk developmental stages1SLEEPER0.27UBQ0.26CYP0.331.05PLA0.51PLA1.93
2PLA0.27PLA0.27GAPDH0.421.33UBQ0.51UBQ2.21
3UBQ0.35ACT50.33NAC0.451.33SLEEPER0.54SLEEPER2.91
4ACT50.43SLEEPER0.34UBQ0.562.01ACT50.54CYP3.50
5GAPDH0.47CYP0.41ACT50.651.97CYP0.56ACT53.94
6CYP0.48GAPDH0.44SLEEPER0.742.31GAPDH0.58GAPDH4.36
7NAC0.51NAC0.47PLA0.752.30NAC0.62NAC5.66
8IDH0.59IDH0.74IDH1.083.24IDH0.81IDH8.00
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tang, J.; Li, E.; Liu, J.; Zhang, Z.; Hua, B.; Jiang, J.; Miao, M. Selection of Reliable Reference Genes for Gene Expression Normalization in Sagittaria trifolia. Genes 2023, 14, 1321. https://doi.org/10.3390/genes14071321

AMA Style

Tang J, Li E, Liu J, Zhang Z, Hua B, Jiang J, Miao M. Selection of Reliable Reference Genes for Gene Expression Normalization in Sagittaria trifolia. Genes. 2023; 14(7):1321. https://doi.org/10.3390/genes14071321

Chicago/Turabian Style

Tang, Jing, Enjiao Li, Jiexia Liu, Zhiping Zhang, Bing Hua, Jiezeng Jiang, and Minmin Miao. 2023. "Selection of Reliable Reference Genes for Gene Expression Normalization in Sagittaria trifolia" Genes 14, no. 7: 1321. https://doi.org/10.3390/genes14071321

APA Style

Tang, J., Li, E., Liu, J., Zhang, Z., Hua, B., Jiang, J., & Miao, M. (2023). Selection of Reliable Reference Genes for Gene Expression Normalization in Sagittaria trifolia. Genes, 14(7), 1321. https://doi.org/10.3390/genes14071321

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop