Next Article in Journal
Variants in KMT2A in Three Individuals with Previous Suspicion of 22q11.2 Deletion Syndrome
Previous Article in Journal
The Impact of Cancer-Associated Fibroblasts on the Biology and Progression of Colorectal Carcinomas
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 15
Issue 2
10.3390/genes15020210
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR

by
Fangyan Wang
1,†,
Peng Li
1,†,
Qiuxu Liu
2,
Gang Nie
1,
Yongqun Zhu
2,* and
Xinquan Zhang
1,*
1
College of Grassland Science and Technology, Sichuan Agricultural University, Chengdu 611130, China
2
Institute of Agricultural Resources and Environment, Sichuan Academy of Agricultural Sciences, Chengdu 610066, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Genes 2024, 15(2), 210; https://doi.org/10.3390/genes15020210
Submission received: 8 January 2024 / Revised: 28 January 2024 / Accepted: 3 February 2024 / Published: 6 February 2024
(This article belongs to the Section Plant Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Quantitative reverse transcription PCR (qRT-PCR) can screen applicable reference genes of species, and reference genes can be used to reduce experimental errors. Sudan grass (Sorghum sudanense (Piper) Stapf) is a high-yield, abiotic-tolerant annual high-quality forage with a wide range of uses. However, no studies have reported reference genes suitable for Sudan grass. Therefore, we found eight candidate reference genes, including UBQ10, HIS3, UBQ9, Isoform0012931, PP2A, ACP2, eIF4α, and Actin, under salt stress (NaCl), drought stress (DR), acid aluminum stress (AlCl3), and methyl jasmonate treatment (MeJA). By using geNorm, NormFinder, BestKeeper, and RefFinder, we ranked eight reference genes on the basis of their expression stabilities. The results indicated that the best reference gene was PP2A under all treatments. eIF4α can be used in CK, MeJA, NaCl, and DR. HIS3 can serve as the best reference gene in AlCl3. Two target genes (Isoform0007606 and Isoform0002387) belong to drought-stress-response genes, and they are highly expressed in Sudan grass according to transcriptome data. They were used to verify eight candidate reference genes under drought stress. The expression trends of the two most stable reference genes were similar, but the trend in expression for Actin showed a significant difference. The reference genes we screened provided valuable guidance for future research on Sudan grass.

1. Introduction

Sudan grass (Sorghum sudanense (Piper) Stapf) belongs to the sorghum family, native to Africa. It has many tillers, good palatability, and high nutritional value while suppressing weeds and increasing soil organic matter [1,2,3]. It is a favorite food for cattle, sheep, and fish. After consuming Sudan grass, fresh milk production can be improved, which can bring high economic benefits [4]. Therefore, Sudan grass is an important feed crop for the development of animal husbandry in many countries [5]. Marginal land is a kind of land with a poor climate, poor physical properties, and difficult to cultivate. It is widely distributed throughout the world. If marginal land can be effectively used, it can significantly fill the gap in the demand for high-quality forage grass [6]. However, abiotic stresses on marginal land, such as salinity, acid aluminum, and drought, threaten the production of Sudan grass [7,8,9]. Therefore, it is urgent to research the mechanism of resistance of Sudan grass to abiotic stress. To study gene expression in more detail and depth, it is important to pinpoint appropriate reference genes for Sudan grass.
Quantitative reverse transcription PCR (qRT-PCR) is often used in the field of molecular biology because of its strong specificity and high sensitivity [10]. In order to achieve precise test results, it is essential to employ suitable reference genes to normalize and standardize the expression levels of target genes [11]. At present, the ordinarily used evaluation method of reference genes is based on Ct values of qRT-PCR; widely used software includes geNorm ((ver.3.5), NormFinder (ver.0.953), and BestKeeper (ver.1.0) [12].
It is widely recognized that the optimal reference genes can vary depending on the species and specific situations [13]. TUB-α and PGK are optimal reference genes in Toona ciliate, but PP2C59 and UBC5B are the most stable reference genes in leaves, and HIS1 and ACT7 can serve as reference genes in young stems under Hypsipyla robusta treatment [14]. GADPH is not optimal as a reference gene for ABA treatment in Isodon rubescens, while it is suitable to be used in the roots of Avena sativa L. under drought stress and the leaves under salt stress [15,16]. Under low-temperature pressure, the optimal reference genes in leaves of Caucasian clover (Trifolium ambiguum M. Bieb.) are APA and EFTu-GTP, and RCD1 and NLI2 are the most stable in the roots [17]. Meanwhile, different analysis methods will lead to different results due to their different algorithms [18]. GAPC and TUB2 can serve as reference genes in various developmental stages of purple- and white-flower genotype (DPW) samples of Allium wallichii by GeNorm, TUB2 and EF1α are optimal in DPW by NormFinder, and UBC is the best in DPW by BestKeeper [19]. Hence, identifying the most suitable reference genes among various plant varieties, growth stages, experimental treatments, and other conditions is important.
Usually, when studying gene expression, we will choose more than one gene to serve as reference genes for result correction to minimize experimental errors [20]. However, Web of Science (www.webofscience.com) did not include any information about the reference genes of Sudan grass. Only sorghum, which is in the same genus as Sudan grass, had reports on the reference genes. Therefore, we will fill the gap of the reference genes of Sudan grass through this experiment. In this study, eight reference genes were assessed through geNorm, NormFinder, BestKeeper, and RefFinder, namely, Eukaryotic initiation factor 4 α (eIF4α), Histone H3 (HIS3), Polyubiquitin 9 (UBQ9), Polyubiquitin 10 (UBQ10), Protein phosphatase 2A (PP2A), Actin gene (Actin), Acyl carrier protein 2 (ACP2), and Isoform0012931, in salt stress (NaCl), drought stress (DR), acid aluminum stress (AlCl3), and methyl jasmonate treatment (MeJA), to facilitate future research on Sudan grass.

2. Materials and Methods

2.1. Plant Material

Sudan grass cv. Chuansu No. 1, which was selected from the Sichuan Academy of Agricultural Sciences, was used in our research. Sudan grass seeds were sterilized by immersing them in 10% sodium hypochlorite solution for 5 min. Afterward, all of the Sudan grass seeds were rinsed three times with sterile distilled water to ensure thorough cleaning. Then, the seeds were sown in quartz-filled plastic containers (30 cm × 20 cm) for germination. Distilled water was poured into the containers until just over the quartz, and the containers were covered with plastic wrap to prevent the water from evaporating quickly. After germination, the distilled water was poured away. Then, the Sudan grass was watered with full-strength Hoagland’s solution. The plastic containers were kept in an illumination incubator, the temperature was 25 °C, the light intensity was 100 µmol/(m·s), and the day/night cycle was 16 h and 8 h at Sichuan Agricultural University.

2.2. Material Treatments

When the number of leaves reached 4–6, the materials were subjected to stress treatments. For acid aluminum stress, seedlings were cultivated with 0.1 mol/L AlCl3 for 5 days. For methyl jasmonate treatment, seedlings were grown with 200 µmol/L methyl jasmonate for 5 days. To induce salt stress, the seedlings underwent a 5-day cultivation period with 200 mmol/L NaCl. To induce drought stress, seedlings underwent 15% (w/v) PEG 6000 for 5 days. Meanwhile, the control group was irrigated with Hoagland’s solution. Each treatment was tested three times independently. Samples were harvested at 0, 1, 2, 3, and 5 days under the five treatments. Then, all the samples were stored at −80 °C.

2.3. Total RNA Extraction and First-Strand cDNA Synthesis

Total RNA was extracted from the Sudan grass material using the HiPure Plant RNA Mini Kit (Magen Biotechnology Co., Ltd., Guangzhou, China). Afterward, the purity and concentration of total RNA were measured by a spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Total RNA (1 µg) was reverse-transcribed using an M5 Super plus qPCR RT kit with gDNA remover (Mei5 Biotechnology Co., Ltd., Beijing, China) into cDNA. Then, the cDNA samples were stored at −20 °C.

2.4. Screening Candidate Reference Genes and PCR Primer Design

Through preliminary experiments, seven common reference genes (eIF4α, PP2A, HIS3, UBQ9, UBQ10, Actin, and ACP2) and one unigene obtained from the Sudan grass transcriptome data (Isoform0012931) were selected. Primer Premier version 5.00 was used to design the primers for the qRT-PCR (Table 1).

2.5. qRT-PCR Amplification Procedure

The qRT-PCR was carried out with a Bio-Rad CFX96 quantitative PCR instrument (BIO-RAD, Hercules, CA, USA). The 10 μL reaction mixture contained 1 μL of quintuple-diluted cDNA, 5 μL of 2× M5 HiPer SYBR Premix EsTaq (Mei5 Biotechnology Co., Ltd., Beijing, China), 0.5 µL of reverse primer and forward primer (10 µmol/L) (total 1 µL), and 3 μL of DNase/RNase-free water. The amplification procedure was 95 °C for 30 s, 40 cycles of 95 °C for 3 s, and 60 °C for 30 s. In order to confirm the specificity of the primers, the melting curve was included. The technical and biological samples were each run in triplicate for all qRT-PCRs. Each gene was able to have a standard curve established through the process of quadruple continuous dilution of the original cDNA, and 4−1, 4−2, 4−3, 4−4, 4−5 were used as the template. The formula for the amplification efficiency (E) was E = 10−1/slope.

2.6. Data Analysis and Validation of Selected Candidate Reference Genes

The stability of the eight reference genes was assessed through geNorm (https://seqyuan.shinyapps.io/seqyuan_prosper/, accessed on 12 October 2023), NormFinder [21], and BestKeeper [12]. Then, RefFinder (http://blooge.cn/RefFinder/?type=reference, accessed on 10 October 2023) was used to evaluate the stability comprehensively. The two top-ranked genes (eIF4α and PP2A) for drought stress and one low-ranked gene (Actin) were used to calculate target genes’ (Isoform0007606 and Isoform0002387) expression levels to verify the reliability of the reference genes.

3. Results

3.1. Verification of Primer Specificity and Effectiveness

We used agarose gel (1%) electrophoresis to identify the specificity of the primers. The results indicated that eight candidate reference genes (eIF4α, PP2A, HIS3, UBQ9, UBQ10, Actin, ACP2, and Isoform0012931) were amplified with a single band (Figure 1), and no primer dimer existed.
Moreover, qRT-PCR was carried out using the cDNA template of Sudan grass under different stresses, and the melting curves were analyzed. It was found that each melting curve indicated a single peak with good repeatability (Figure 2). As shown in Table 1, E ranged from 88.2% to 116.7%, and the correlation coefficient (R2) ranged between 0.972 and 0.996. Based on these results, it appears that the qRT-PCR primers demonstrated satisfactory specificity and efficiency, which means they were suitable for use in subsequent experiments.

3.2. Expression Distribution of Eight Candidate Reference Genes

According to the findings, Ct values were used to indicate the expression levels of the reference genes, with higher gene expression levels corresponding to lower Ct values. After analyzing the Ct values we obtained from Sudan grass, it was observed that the range of Ct values was between 19.73 and 36.03. The ranking of the expression levels was UBQ10 > HIS3 > UBQ9 > Isoform0012931 > PP2A > ACP2 > eIF4α > Actin (Figure 3). Among the eight reference genes, UBQ10 had the highest mean Ct value (23.66), while Actin had a low expression level in comparison with the other genes (27.27).

3.3. Expression Stability Analysis of Candidate Reference Genes under Five Treatments

We employed four analytical tools to evaluate the eight reference genes in this research. Comprehensive rankings of the eight reference genes are listed in Table 2. We found that different reference genes exhibited different expression stability under different treatments. Meanwhile, using different analytics software can result in different rankings.

3.3.1. GeNorm Analysis

GeNorm was created by Vandesompele et al. in 2002 to assist with screening reference genes, determining the ideal quantity of reference genes, and assessing reference genes by calculating their M value [22]. The reference gene’s stability improves as the M value decreases. In this study, PP2A and eIF4α were the optimal genes with the same M values of 0.188, 0.331, 0.438, and 0.284 in CK, MeJA, NaCl, and DR, while Actin was the least stable gene, with M values of 0.476, 0.639, and 0.674 in CK, MeJA, and DR. UBQ10 had the highest M value of 1.217 in NaCl. For acid aluminum stress, UBQ10 and Isoform0012931 had the highest stability, whereas ACP2 had the lowest. Out of all the treatments, PP2A and eIF4α exhibited higher stability, whereas Actin had the lowest (Figure 4A–E).
Additionally, the paired variant V value of the standardized factor can be calculated by the software after introducing a new reference gene. If Vn/Vn+1 < 0.15, then the ideal number of reference genes would be n. On the other hand, if Vn/Vn+1 > 0.15, then the number would be n + 1. Based on the data in Figure 5, it appears that using just two reference genes is enough to normalize target genes across various treatments, as all pairwise variations were below 0.15.

3.3.2. NormFinder Analysis

NormFinder was designed by Andersen et al. in 2004 to screen reference genes for stability [21]. It first evaluates the variance of reference genes within and between groups. Then, the reference gene with the lowest expression stability value is the most suitable [23]. The stability values of the five treatments are shown in Figure 6A. eIF4α was the top rank in CK, NaCl, and DR with stability values of 0.17, 0.181, and 0.098; HIS3 was the top-ranked gene in MeJA with a value of 0.177, and PP2A was the top in AlCl3 with a value of 0.258. Actin had the lowest values in CK, MeJA, and DR, which is the same as geNorm.

3.3.3. BestKeeper Analysis

BestKeeper is a program designed by Michael et al. to rank reference genes and target genes based on their stabilities. Through the calculations of this program, we can obtain the correlation coefficient (r), standard deviation (SD), and coefficient of variation (CV) of a pairing. A high value of r corresponds to a low SD and CV, indicating greater stability of the reference gene [12]. Figure 6B displays the findings. eIF4α was the top rank in CK. PP2A was the top rank in MeJA and DR. ACP2 was the top rank in AlCl3 and NaCl. However, Actin was the last one in all treatments except AlCl3. And for AlCl3, UBQ10 was the most unstable gene.

3.3.4. Comparative ΔCt Analysis

To determine the expression stability differences between samples, the mean standard deviation (SD) was ranked using the ΔCt method, the stability increasing as the SD decreased [24]. eIF4α had the highest stability in CK, NaCl, and DR with SD of 0.38, 0.9, and 0.52. HIS3 was the top rank in MeJA, with an SD of 0.52. PP2A was at the top of the ranking in AlCl3 with an SD of 0.54 (Figure 6C). Actin had the lowest values in CK, MeJA, and DR, which is the same as for geNorm and NormFinder.

3.3.5. RefFinder Analysis

RefFinder is an innovative analysis software that utilizes a combination of four computational methods and calculates the geometric mean of gene weights to provide a comprehensive ranking. Stable gene expression is typically associated with a lower geometric mean [25]. In CK, NaCl, and DR, eIF4α was the top-ranked gene. PP2A was the top rank in MeJA and AlCl3. Nevertheless, Actin was the most unstable gene in four stresses, CK, MeJA, Al, and DR (Figure 6D). And for NaCl, UBQ10 was the most unstable.

3.4. Validation of Candidate Reference Genes

To ensure the dependability of the rankings for the reference genes, the relative expression trend in the target genes was analyzed using the two most stable reference genes (eIF4α, PP2A), one combined reference gene (eIF4α + PP2A), and the most unstable reference gene (Actin) under drought stress. Two target genes (Isoform0007606 and Isoform0002387) that belong to the drought-stress-response genes are highly expressed according to transcriptomic data of Sudan grass under drought stress [26]. They belong to COR (Isoform0007606) and MAPK (Isoform0002387). The polypeptide encoded by the COR gene can form a hydrophilic lipid α-helix, which can stabilize the structure of the cell membrane and make the plant resistant when the plant is dehydrated by low temperature, drought, or osmotic stress [27]. Under drought stress, cotton activates the MAPK cascade GhMAP3K62-GhMKK16-GhMPK32 pathway and phosphorylates downstream target gene GhEDT1, thereby activating the expression of ABA synthesizing gene GhNCED3, promoting ABA accumulation, and thus, participating in the molecular mechanism of drought response [28]. Therefore, they are appropriate for the target genes.
In Figure 7, we used the 2−ΔΔCt method to normalize the target genes. When Actin served as a reference gene, we discovered that the relative expression values of the two target genes were higher than those of other groups. When Isoform0007606 was the target gene, the relative expression values using eIF4α, PP2A, eIF4α + PP2A, and Actin were higher on the first day compared to the beginning. Then, the values of eIF4α, PP2A, and eIF4α + PP2A started a steady decline while the values of Actin started to rise (Figure 7A). As shown in Figure 7B, when using Actin, the expression trend of Isoform0002387 continued to rise from the beginning, and by the third day, it was 1.7 times that of beginning Nevertheless, the values of eIF4α, PP2A, and eIF4α + PP2A started a steady decline from the beginning to the second day. There was a slight rise on the third day.

4. Discussion

Ideally, reference genes remain consistent across all experimental conditions and are stable in all tissues and growth stages, but it has been shown in many studies that stable expression genes will change across diverse species and conditions [14,29]. An effective way to analyze gene expression is through the use of qRT-PCR, known for its high efficiency. However, choosing appropriate reference genes is crucial when analyzing gene expression levels, as the accuracy of the measurements will be impacted by the RNA quality, efficiency of reverse transcription, or cDNA quality [30].
Sudan grass is an inexpensive, fast-growing plant that has great forage value, but there are not studies on reference gene screening of Sudan grass under different stresses [9]. Therefore, we screened eight reference genes from transcriptome data and common reference genes using five algorithms. As can be seen from Table 2, most of the rankings were different between the five algorithms. However, comparative ΔCt and NormFinder had the same order for CK, methyl jasmonate treatment, and salt stress. Meanwhile, comparative ΔCt and RefFinder had the same ranking for DR. Under aluminum stress, UBQ9 ranked second in BestKeeper, but in the other four algorithms it ranks very low. In comparative ΔCt, NormFinder, and RefFinder, HIS3 was the top rank in MeJA, but it was an unstable reference gene in BestKeeper. So, the ranking of BestKeeper is always significantly different from those of the other four algorithms.
Due to drought and the application of large amounts of fertilizer or other factors, global soil salinization is becoming increasingly serious, and there are about 10 million km2 of saline–alkali land in the world. Therefore, more and more attention has been paid to the research of saline soils [31]. It can be seen from the result of the ranking that eIF4α and PP2A are the most stable under CK, MeJA, NaCl, and DR eIF4α belongs to the eukaryotic initiation factor 4α family. It is required for the binding of mRNA to 40S ribosomal subunits. It can also unwind double-stranded RNA [32]. PP2A belongs to the phosphoprotein phosphatase family of non-metal dependent serine/threonine protein phosphatases. They are universal in eukaryotes. PP2A regulates plant vesicle traffic and various processes of plant development and stress responses [33]. The results are consistent with the sequencing results of reference genes of sorghum. For example, Ma et al. found that GAPDH and eIF4α had the highest expression stability in leaves of sorghum under low nitrogen treatment, while Actin had the lowest [34]. Palakolanu et al. found that PP2A and CYP were the reference genes with the best performance under abiotic stress in sorghum [35]. María et al. found that ARI 8 and PP2A were the top ranks under glyphosate stress in sorghum [36]. Otherwise, eIF4α and PP2A always express stably in Larix olgensis, Vaccinium corymbosum × angustifolium, Chenopodium quinoa Willd, oat, Suaeda glauca L., and so on, especially in salt stress and drought stress [37,38,39,40,41].
However, in our study, Actin is the most unstable gene in most treatments. This result is similar to the research of Zhang et al. in Papaver somniferum L. but contrary to the study of You et al. in Prunus persica [42,43]. Studies on the reference genes in various plants at different growth stages showed that Actin and UBQ were the most stable [44]. But our study researched different abiotic stresses, which is probably the reason why Actin and UBQ express unstably. Additionally, the reference gene we screened from the transcriptome data (Isoform0012931) usually had a poor performance, always being ranked sixth or seventh. Nevertheless, its performance was better in aluminum stress and it could be ranked in the top three.
In conclusion, we identified eight reference genes in Sudan grass under four abiotic stresses and used five arithmetics for evaluation. Because of the use of different algorithms, the ranking results are also different. The best reference gene was PP2A under all treatments. eIF4α can be used as a reference gene in CK, MeJA, NaCl, and DR. HIS3 can be used as an optimal reference gene in AlCl3. Actin was the most unstable in most treatments. Our research can provide data support for the next step of functional analysis of Sudan grass. However, it still has room for improvement; we will keep investigating more reference genes in various tissues, under various stresses, and in various varieties of Sudan grass.

Author Contributions

Conceptualization, Q.L. and F.W.; data curation, Q.L., P.L. and F.W.; writing—original draft preparation, F.W. and P.L.; writing—review and editing, G.N.; funding acquisition, Y.Z. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the China Agriculture Research System of MOF and MARA (CARS-34), the Sichuan Province Breeding Research grant (2021YFYZ0013), Sichuan Forage Innovation Team of the National Modern Agricultural Industry Technology System (SCCXTD-2020-16), the China Agriculture Research System of MOF and MARA (CARS-34-51), and Achievement Transformation Project of Sichuan Science and Technology Department (2021NZZJ0004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Liu, M.X.; Wang, Y.W.; Han, J.G.; Mao, P.S. Phenolic Compounds from Chinese Sudangrass, Sorghum, Sorghum-Sudangrass Hybrid, and Their Antioxidant Properties. Crop Sci. 2011, 51, 247–258. [Google Scholar] [CrossRef]
  2. Rangaswami Ayyangar, G.N.; Ponnaiya, B.W.N. Studies in Sorghum sudanense, Stapf-The Sudan grass. Proc. Natl. Acad. Sci. India 1939, 10, 237–254. [Google Scholar] [CrossRef]
  3. Zamfir, M.C.; Schitea, M.; Zamfir, I. The variability study of some quantitative traits in Sudan grass Sorghum sudanense Piper. (Staph.). Rom. Agric. Res. 2001, 1, 23–29. [Google Scholar]
  4. Chen, Q.; Wei, S.S.; Deng, Z.R.; Yin, L.P.; He, B.; Kong, X.L. Optimization of DNA Extraction from Seeds of Sorghum sudanense (Piper) Stapf. Not. Bot. Horti Agrobot. Cluj-Napoca 2009, 37, 256–260. [Google Scholar]
  5. Han, P.A.; Lu, X.P.; Mi, F.G.; Dong, J.; Xue, C.L.; Li, J.K.; Han, B.; Zhang, X.Y. Proteomic analysis of heterosis in the leaves of sorghum-sudangrass hybrids. Acta Biochim. Biophys. Sin. 2016, 48, 161–173. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, Y.; Guo, X.; Liu, J.; Zhang, Z. Assessment of marginal land potential for energy plants in China. Land Dev. Eng. Res. 2017, 2, 1–7. [Google Scholar]
  7. Bibi, A.; Sadaqat, H.A.; Akram, H.M.; Khan, T.M.; Usman, B.F. Physiological and agronomic responses of sudangrass to water stress. J. Agric. Res. 2010, 48, 369–380. [Google Scholar]
  8. Buchanan, C.D.; Lim, S.Y.; Salzman, R.A.; Kagiampakis, L.; Morishige, D.T.; Weers, B.D.; Klein, R.R.; Pratt, L.H.; Cordonnier-Pratt, M.M.; Klein, P.E.; et al. Sorghum bicolor’s transeriptome response to dehydration, high salinity and ABA. Plant Mol. Biol. 2005, 58, 699–720. [Google Scholar] [CrossRef] [PubMed]
  9. Saraeian, A.; Hadi, A.; Raji, F.; Ghassemi, A.; Johnson, M. Cadmium removal from aqueous solution by low-cost native and surface modified Sorghum × drummondii (Sudangrass). J. Environ. Chem. Eng. 2018, 6, 3322–3331. [Google Scholar] [CrossRef]
  10. Huggett, J.; Dheda, K.; Bustin, S.; Zumla, A. Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6, 279–284. [Google Scholar] [CrossRef]
  11. Yuan, W.; Wan, H.; Yang, Y. Characterization and Selection of Reference Genes for Real-time Quantitative RT-PCR of Plants. Chin. Bull. Bot. 2012, 47, 427–436. [Google Scholar]
  12. 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]
  13. Zhang, K.; Li, M.N.; Cao, S.H.; Sun, Y.; Long, R.C.; Kang, J.M.; Yan, L.; Cui, H.T. Selection and validation of reference genes for target gene analysis with quantitative real-time PCR in the leaves and roots of Carex rigescens under abiotic stress. Ecotoxicol. Environ. Saf. 2019, 168, 127–137. [Google Scholar] [CrossRef]
  14. Song, H.Y.; Mao, W.M.; Duan, Z.H.; Que, Q.M.; Zhou, W.; Chen, X.Y.; Li, P. Selection and validation of reference genes for measuring gene expression in Toona ciliata under different experimental conditions by quantitative real-time PCR analysis. BMC Plant Biol. 2020, 20, 450. [Google Scholar] [CrossRef]
  15. Lian, C.L.; Zhang, B.; Yang, J.F.; Lan, J.X.; Yang, H.; Guo, K.H.; Li, J.J.; Chen, S.Q. Validation of suitable reference genes by various algorithms for gene expression analysis in Isodon rubescens under different abiotic stresses. Sci. Rep. 2022, 12, 19599. [Google Scholar] [CrossRef]
  16. Tajti, J.; Pal, M.; Janda, T. Validation of Reference Genes for Studying Different Abiotic Stresses in Oat (Avena sativa L.) by RT-qPCR. Plants 2021, 10, 1272. [Google Scholar] [CrossRef] [PubMed]
  17. Yin, X.; He, T.; Yi, K.; Zhao, Y.; Hu, Y.; Liu, J.; Zhang, X.; Meng, L.; Wang, L.; Liu, H.; et al. Comprehensive evaluation of candidate reference genes for quantitative real-time PCR-based analysis in Caucasian clover. Sci. Rep. 2021, 11, 3269. [Google Scholar] [CrossRef]
  18. Liu, Q.; Qi, X.; Yan, H.; Huang, L.; Nie, G.; Zhang, X. Reference Gene Selection for Quantitative Real-Time Reverse-Transcriptase PCR in Annual Ryegrass (Lolium multiflorum) Subjected to Various Abiotic Stresses. Molecules 2018, 23, 172. [Google Scholar] [CrossRef]
  19. Lin, Y.; Liu, G.; Rao, Y.; Wang, B.; Tian, R.; Tan, Y.; Peng, T. Identification and validation of reference genes for qRT-PCR analyses under different experimental conditions in Allium wallichii. J. Plant Physiol. 2023, 281, 153925. [Google Scholar] [CrossRef]
  20. Bai, S.; Wang, X.; Guo, M.; Cheng, G.; Khan, A.; Yao, W.; Gao, Y.; Li, J. Selection and Evaluation of Reference Genes for Quantitative Real-Time PCR in Tomato (Solanum lycopersicum L.) Inoculated with Oidium neolycopersici. Agronomy 2022, 12, 3171. [Google Scholar] [CrossRef]
  21. Andersen, C.L.; Jensen, J.L.; Ørntoft, 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] [PubMed]
  22. 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.1. [Google Scholar] [CrossRef] [PubMed]
  23. Hu, A.; Yang, X.; Zhu, J.; Wang, X.; Liu, J.; Wang, J.; Wu, H.; Zhang, H.; Zhang, H. Selection and validation of appropriate reference genes for RT-qPCR analysis of Nitraria sibirica under various abiotic stresses. BMC Plant Biol. 2022, 22, 592. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, Y.; Zhu, L.; Xue, J.; Yang, J.; Hu, H.; Cui, J.; Xu, J. Selection and verification of appropriate reference genes for expression normalization in Cryptomeria fortunei under abiotic stress and hormone treatments. Genes 2021, 12, 791. [Google Scholar] [CrossRef] [PubMed]
  25. Xie, F.; Xiao, P.; Chen, D.; Xu, L.; Zhang, B. miRDeepFinder: A miRNA analysis tool for deep sequencing of plant small RNAs. Plant Mol. Biol. 2012, 80, 75–84. [Google Scholar] [CrossRef]
  26. Liu, Q.; Wang, F.; Xu, Y.; Lin, C.; Li, X.; Xu, W.; Wang, H.; Zhu, Y. Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing. Plants 2023, 12, 2624. [Google Scholar] [CrossRef]
  27. Huang, Y.; Liu, Y.; Zhu, Y.; Liu, Z. Effects of COR gene on plant cold resistance engineering. Guizhou Agric. Sci. 2014, 42, 37–42. [Google Scholar]
  28. Chen, L.; Zhang, B.; Xia, L.; Yue, D.; Han, B.; Sun, W.; Wang, F.; Lindsey, K.; Zhang, X.; Yang, X. The GhMAP3K62-GhMKK16-GhMPK32 kinase cascade regulates drought tolerance by activating GhEDT1-mediated ABA accumulation in cotton. J. Adv. Res. 2023, 51, 13–25. [Google Scholar] [CrossRef]
  29. Razavi, S.A.; Afsharpad, M.; Modarressi, M.H.; Zarkesh, M.; Yaghmaei, P.; Nasiri, S.; Tavangar, S.M.; Gholami, H.; Daneshafrooz, A.; Hedayati, M. Validation of reference genes for normalization of relative qRT-PCR studies in papillary thyroid carcinoma. Sci. Rep. 2019, 9, 15241. [Google Scholar] [CrossRef]
  30. Liu, Y.P.; Zhang, Y.; Liu, F.; Liu, T.; Chen, J.Y.; Fu, G.; Zheng, C.Y.; Su, D.D.; Wang, Y.N.; Zhou, H.K.; et al. Establishment of reference (housekeeping) genes via quantitative real-time PCR for investigation of the genomic basis of abiotic stress resistance in Psammochloa villosa (Poaceae). J. Plant Physiol. 2022, 268, 153575. [Google Scholar] [CrossRef]
  31. Zhang, H.; Yu, F.; Xie, P.; Sun, S.; Qiao, X.; Tang, S.; Chen, C.; Yang, S.; Mei, C.; Yang, D.; et al. A Gγ protein regulates alkaline sensitivity in crops. Science 2023, 379, eade8416. [Google Scholar] [CrossRef]
  32. Rogers, G.W.; Komar, A.A.; Merrick, W.C. eIF4A: The godfather of the DEAD box helicases. In Progress in Nucleic Acid Research and Molecular Biology; Academic Press: Cambridge, MA, USA, 2002; Volume 72, pp. 307–331. [Google Scholar]
  33. Máthé, C.; M-Hamvas, M.; Freytag, C.; Garda, T. The Protein Phosphatase PP2A Plays Multiple Roles in Plant Development by Regulation of Vesicle Traffic—Facts and Questions. Int. J. Mol. Sci. 2021, 22, 975. [Google Scholar] [CrossRef]
  34. Ma, J.; Sun, Y.; Wang, Y.; Duan, Y. Screening of reference genes for qRT-PCR analysis in sorghum (Sorghum bicolor) under low nitrogen stress. J. Agric. Biotechnol. 2017, 25, 805–812. [Google Scholar]
  35. Sudhakar Reddy, P.; Srinivas Reddy, D.; Sivasakthi, K.; Bhatnagar-Mathur, P.; Vadez, V.; Sharma, K.K. Evaluation of sorghum [Sorghum bicolor (L.)] reference genes in various tissues and under abiotic stress conditions for quantitative real-time PCR data normalization. Front. Plant Sci. 2016, 7, 529. [Google Scholar] [CrossRef]
  36. Ulrich, M.N.; Muñiz-Padilla, E.; Corach, A.; Hopp, E.; Tosto, D. Validation of reference genes for quantitative PCR in Johnsongrass (Sorghum halepense L.) under glyphosate stress. Plants 2021, 10, 1555. [Google Scholar] [CrossRef]
  37. Li, D.; Yu, S.; Zeng, M.; Liu, X.; Yang, J.; Li, C. Selection and Validation of Appropriate Reference Genes for Real-Time Quantitative PCR Analysis in Needles of Larix olgensis under Abiotic Stresses. Forests 2020, 11, 193. [Google Scholar] [CrossRef]
  38. Deng, Y.; Li, Y.; Sun, H. Selection of reference genes for RT-qPCR normalization in blueberry (Vaccinium corymbosum × angustifolium) under various abiotic stresses. FEBS Open Bio 2020, 10, 1418–1435. [Google Scholar] [CrossRef] [PubMed]
  39. 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]
  40. Duan, Z.L.; Han, W.H.; Yan, L.; Wu, B. Reference gene selections for real time quantitative PCR analysis of gene expression in different oat tissues and under salt stress. Biol. Plant. 2020, 64, 838–844. [Google Scholar] [CrossRef]
  41. Wang, M.; Ren, T.; Marowa, P.; Du, H.; Xu, Z. Identification and selection of reference genes for gene expression analysis by quantitative real-time PCR in Suaeda glauca’s response to salinity. Sci. Rep. 2021, 11, 8569. [Google Scholar] [CrossRef] [PubMed]
  42. You, S.; Cao, K.; Chen, C.; Li, Y.; Wu, J.; Zhu, G.; Fang, W.; Wang, X.; Wang, L. Selection and validation reference genes for qRT-PCR normalization in different cultivars during fruit ripening and softening of peach (Prunus persica). Sci. Rep. 2021, 11, 7302. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, Z.; Li, C.; Zhang, J.; Chen, F.; Gong, Y.; Li, Y.; Su, Y.; Wei, Y.; Zhao, Y. Selection of the reference gene for expression normalization in Papaver somniferum L. under abiotic stress and hormone treatment. Genes 2020, 11, 124. [Google Scholar] [CrossRef] [PubMed]
  44. Yakupjan, H.; Guan, B.; Zhang, F. Research Progress in Plant Reference Genes. Biotechnol. Bull. 2011, 7, 7–11. [Google Scholar]
Genes 15 00210 g001
Figure 1. PCR products of 8 reference genes. M, DNA marker; 1, eIF4α; 2, PP2A; 3, HIS3; 4, UBQ9; 5, UBQ10; 6, Isoform0012931; 7, ACP2; 8, Actin.
Figure 1. PCR products of 8 reference genes. M, DNA marker; 1, eIF4α; 2, PP2A; 3, HIS3; 4, UBQ9; 5, UBQ10; 6, Isoform0012931; 7, ACP2; 8, Actin.
Genes 15 00210 g001
Genes 15 00210 g002
Figure 2. Melting curves for eight genes of CK. Melting temperatures were visualized by plotting the negative first derivative of fluorescence relative to the temperature in Celsius (–(d/dT)).
Figure 2. Melting curves for eight genes of CK. Melting temperatures were visualized by plotting the negative first derivative of fluorescence relative to the temperature in Celsius (–(d/dT)).
Genes 15 00210 g002
Genes 15 00210 g003
Figure 3. Ct values of 8 reference genes in Sudan grass under five treatments. Different colors represent different reference genes. The reference gene represented by each color corresponds to the name of the reference gene in the vertical coordinate below.
Figure 3. Ct values of 8 reference genes in Sudan grass under five treatments. Different colors represent different reference genes. The reference gene represented by each color corresponds to the name of the reference gene in the vertical coordinate below.
Genes 15 00210 g003
Genes 15 00210 g004
Figure 4. Expression stability values of 8 candidate reference genes validated by geNorm analysis. (A) CK, (B) acid aluminum stress, (C) methyl jasmonate treatment, (D) salt stress, and (E) drought stress.
Figure 4. Expression stability values of 8 candidate reference genes validated by geNorm analysis. (A) CK, (B) acid aluminum stress, (C) methyl jasmonate treatment, (D) salt stress, and (E) drought stress.
Genes 15 00210 g004
Genes 15 00210 g005
Figure 5. Pairwise variation (V) analysis of candidate reference genes.
Figure 5. Pairwise variation (V) analysis of candidate reference genes.
Genes 15 00210 g005
Genes 15 00210 g006
Figure 6. Stability analysis of 8 candidate reference genes using four methods. The color within a pane becomes whiter as the value decreases. (A) Analysis results from NormFinder (B) Analysis results from BestKeeper (C) Analysis results from Delta CT (D) Analysis results from RefFinder.
Figure 6. Stability analysis of 8 candidate reference genes using four methods. The color within a pane becomes whiter as the value decreases. (A) Analysis results from NormFinder (B) Analysis results from BestKeeper (C) Analysis results from Delta CT (D) Analysis results from RefFinder.
Genes 15 00210 g006
Genes 15 00210 g007
Figure 7. Relative expression trends of two target genes under drought stress. The target genes we used were screened from Sudan grass transcriptome data. (A) COR (Isoform0007606), (B) MAPK (Isoform0002387).
Figure 7. Relative expression trends of two target genes under drought stress. The target genes we used were screened from Sudan grass transcriptome data. (A) COR (Isoform0007606), (B) MAPK (Isoform0002387).
Genes 15 00210 g007
Table 1. Primer sequences (forward and reverse 5′-3′) of eight candidate reference genes and two target genes.
Table 1. Primer sequences (forward and reverse 5′-3′) of eight candidate reference genes and two target genes.
Gene SymbolGene DescriptionForward/Reverse Primer (5′-3′)Primers Tm (°C)Amplicon Length (bp)RT-qPCR Efficiency (%)R2
eIF4αEukaryotic initiation factor 4 αCAACTTTGTCACCCGCGATGA
TCCAGAAACCTTAGCAGCCCA		57.57
57.57		17899.50.993
PP2AProtein phosphatase 2AAACCCGCAAAACCCCAGACTA
TACAGGTCGGGCTCATGGAAC		57.57
59.52		176111.30.980
HIS3Histone H3ACTTCAAGACTGATCTGCGTTTCC
CATGGATGGCACAAAGGTTGG		57.86
57.57		146102.20.988
UBQ9Polyubiquitin 9TACAGTTCTACAAGGTGGACGAC
GCAGTAGTGGCGGTCGAAGT		57.77
59.50		15892.70.992
UBQ10Polyubiquitin 10CCGTGGTGGCCAGTAAGTTC
GGACTCAACATGGGCTCTGC		59.50
59.50		20488.20.987
Isoform0012931 ATGGCCAACCGCTGGGTC
CTTGGGGTCAGTGAAGAACTTGT		59.46
57.77		286116.70.972
ACP2Acyl carrier protein 2ACGAACTTGTTGCGGCAGAAG
GAACAAGAAGGGATGCGCTGG		57.57
59.52		128116.30.988
ActinActin geneAAGTGCGACGTGGATATTAGGA
TCTTGGGCGGAAAGAATTAGA		55.81
53.66		344101.50.996
Isoform0007606Cold response geneTTCGGCACTTCCTTCCTCA
GAATAGACGCAGAATAACAGCAATA		58.70
58.40		204
Isoform0002387Mitogen-activated protein kinaseCCGAGCAATTTGTTCCTAA
GCATCATCTGGTGAGCCTAT		53.80
54.70		284
Table 2. Comprehensive rankings of 8 candidate reference genes.
Table 2. Comprehensive rankings of 8 candidate reference genes.
Rank12345678
(A) CK
Delta CTeIF4αPP2AUBQ9HIS3UBQ10ACP2Isoform0012931Actin
BestKeepereIF4αUBQ10PP2AHIS3UBQ9Isoform0012931ACP2Actin
NormfindereIF4αPP2AUBQ9HIS3UBQ10ACP2Isoform0012931Actin
GeNormPP2A|eIF4αUBQ9UBQ10HIS3Isoform0012931ACP2Actin
RefFindereIF4αPP2AUBQ9UBQ10HIS3ACP2Isoform0012931Actin
(B) Salt stress
Delta CTeIF4αPP2AHIS3ACP2UBQ9Isoform0012931ActinUBQ10
BestKeeperACP2UBQ9eIF4αPP2AHIS3Isoform0012931UBQ10Actin
NormfindereIF4αPP2AHIS3ACP2UBQ9Isoform0012931ActinUBQ10
GeNormPP2A|eIF4αHIS3UBQ9ACP2Isoform0012931ActinUBQ10
RefFindereIF4αPP2AACP2HIS3UBQ9Isoform0012931ActinUBQ10
(C) Drought stress
Delta CTeIF4αPP2AUBQ10ACP2HIS3UBQ9Isoform0012931Actin
BestKeeperPP2AUBQ10HIS3eIF4αIsoform0012931ACP2UBQ9Actin
NormfindereIF4αPP2AACP2UBQ10HIS3UBQ9Isoform0012931Actin
GeNormPP2A|eIF4αACP2UBQ10HIS3UBQ9Isoform0012931Actin
RefFindereIF4αPP2AUBQ10ACP2HIS3UBQ9Isoform0012931Actin
(D) Aluminum stress
Delta CTPP2AHIS3Isoform0012931UBQ10eIF4αActinUBQ9ACP2
BestKeeperACP2UBQ9eIF4αPP2AHIS3ActinIsoform0012931UBQ10
NormfinderPP2AHIS3Isoform0012931eIF4αUBQ10ActinACP2UBQ9
GeNormIsoform0012931|UBQ10HIS3PP2AActineIF4αUBQ9ACP2
RefFinderPP2AHIS3Isoform0012931UBQ10eIF4αACP2UBQ9Actin
(E) Methyl jasmonate treatment
Delta CTHIS3eIF4αUBQ10PP2AUBQ9ACP2Isoform0012931Actin
BestKeeperPP2AUBQ10ACP2UBQ9eIF4αHIS3Isoform0012931Actin
NormfinderHIS3eIF4αUBQ10PP2AUBQ9ACP2Isoform0012931Actin
GeNormPP2A|eIF4αUBQ10HIS3UBQ9ACP2Isoform0012931Actin
RefFinderPP2AeIF4αHIS3UBQ10UBQ9ACP2Isoform0012931Actin
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

Wang, F.; Li, P.; Liu, Q.; Nie, G.; Zhu, Y.; Zhang, X. Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR. Genes 2024, 15, 210. https://doi.org/10.3390/genes15020210

AMA Style

Wang F, Li P, Liu Q, Nie G, Zhu Y, Zhang X. Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR. Genes. 2024; 15(2):210. https://doi.org/10.3390/genes15020210

Chicago/Turabian Style

Wang, Fangyan, Peng Li, Qiuxu Liu, Gang Nie, Yongqun Zhu, and Xinquan Zhang. 2024. "Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR" Genes 15, no. 2: 210. https://doi.org/10.3390/genes15020210

APA Style

Wang, F., Li, P., Liu, Q., Nie, G., Zhu, Y., & Zhang, X. (2024). Selection and Validation of Reference Genes in Sudan Grass (Sorghum sudanense (Piper) Stapf) under Various Abiotic Stresses by qRT-PCR. Genes, 15(2), 210. https://doi.org/10.3390/genes15020210

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