Selection of suitable reference genes for quantitative RT-PCR normalization in the halophyte Halostachys caspica under salt and drought stress

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.

Halostachys caspica pathogenesis-related protein 10 acts as a cytokinin reservoir to regulate plant growth and development

Pathogenesis-related class 10 (PR-10) proteins play a role in plant growth and development, but the underlying molecular mechanisms are unclear. Here, we isolated a salt-induced PR-10 gene from the halophyte Halostachys caspica and named it HcPR10. HcPR10 was constitutively expressed during development and HcPR10 localized to the nucleus and cytoplasm. HcPR10-mediated phenotypes including bolting, earlier flowering, increased branch number and siliques per plant are highly correlated with increased cytokinin levels in transgenic Arabidopsis. Meanwhile, increased levels of cytokinin in plants is temporally correlated with HcPR10 expression patterns. Although the expression of cytokinin biosynthesis genes validated was not upregulated, cytokinin-related genes including chloroplast-related genes, cytokinin metabolism and cytokinin responses genes and flowering-related genes were significantly upregulated in the transgenic Arabidopsis compared to the wild type by transcriptome deep sequencing. Analysis of the crystal structure of HcPR10 revealed a trans-zeatin riboside (a type of cytokinin) located deep in its cavity, with a conserved conformation and protein-ligand interactions, supporting HcRP10 acts as a cytokinin reservoir. Moreover, HcPR10 in Halostachys caspica predominantly accumulated in vascular tissue, the site of long-distance translocation of plant hormones. Collectively, we draw that HcPR10 as a cytokinin reservoir induces cytokinin-related signal transduction in plants, thereby promoting plant growth and development. These findings could provide intriguing insights into the role of HcPR10 proteins in phytohormone regulation in plants and advance our understanding of cytokinin-mediated plant development and could facilitate the breeding of transgenic crops with earlier mature, higher yielding agronomic traits.

The miR169b/NFYA1 module from the halophyte Halostachys caspica endows salt and drought tolerance in Arabidopsis through multi-pathways

Salt and drought are the major abiotic stress factors plaguing plant growth, development and crop yields. Certain abiotic-stress tolerant plants have developed special mechanisms for adapting to adverse environments in the long process of evolution. Elucidating the molecular mechanisms by which they can exert resistance to abiotic stresses is beneficial for breeding new cultivars to guide agricultural production. Halostachys caspica, a perennial halophyte belonging to Halostachys in Amaranthaceae, is extremely tolerant to harsh environments, which is commonly grown in the saline-alkali arid desert area of Northwest, China. However, the molecular mechanism of stress tolerance is unclear. Nuclear Factor Y-A (NFYA) is a transcription factor that regulates the expression of downstream genes in plant response to adverse environments. It has also been reported that some members of the NFYA family are the main targets of miR169 in plants. In this study, we mainly focused on exploring the functions and preliminary mechanism of the miR169b/NFYA1 module from H. caspica to abiotic stress. The main results showed that RLM-RACE technology validated that HcNFYA1 was targeted by HcmiR169b, qRT-PCR revealed that HcmiR169b was repressed and HcNFYA1 was induced in the H. caspica branches under various abiotic stress as well ABA treatment and Arabidopsis stable transformation platform with molecular methods was applied to elucidate that the HcmiR169b/HcNFYA1 module conferred the salt and drought tolerance to plants by enhancing ABA synthesis and ABA signal transduction pathways, maintaining ROS homeostasis and the stability of cell membrane. HcNFYA1 is expected to be a candidate gene to improve plant resistance to salt and drought stresses.

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

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

Selection and Validation of Reference Genes for RT-qPCR Analysis in Aegilops tauschii (Coss.) under Different Abiotic Stresses

Aegilops tauschii (Coss.) is an aggressive and serious annual grass weed in China. Its DD genome is a rich source of genetic material and performs better under different abiotic stress conditions (salinity, drought, temperature, etc.). Reverse-transcribed quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for reference gene selection and validation. This work aimed to evaluate the stability of reference gene expression in Ae. tauschii under different abiotic stresses (salinity, drought, hot, and cold) and developmental stages (seedling and development). The results show that the ubiquitin-conjugating enzyme E2 36-like (UBC36) and protein microrchidia 2-like (HSP) are the most stable genes under control and salinity conditions, respectively. Under drought stress conditions, UBC36 is more stable as compared with others. Glyceraldehyde-3-phosphate dehydrogenase (GADPH) is the most stable reference gene during heat stress conditions and thioredoxin-like protein (YLS) under cold stress condition. Phosphate2A serine/threonine-protein phosphatase 2A (PP2A) and eukaryotic translation initiation factor 3 (ETIF3) are the most stable genes at seedling and developmental stages. Intracellular transport protein (CAC) is recommended as the most stable gene under different abiotic stresses and at developmental stages. Furthermore, the relative expression levels of NHX1 and DREB under different levels of salinity and drought stress conditions varied with the most (HSP and UBC36) and least (YLS and ACT) stable genes. This study provides reliable reference genes for understanding the tolerance mechanisms in Ae. tauschii under different abiotic stress conditions.

Identification of HvLRX, a new dehydration and light responsive gene in Tibetan hulless barley (Hordeum vulgare var. nudum)

BACKGROUND

Tibetan hulless barley (Hordeum vulgare var. nudum), adjusting to the harsh environment on Qinghai-Tibet Plateau, is a good subject for analyzing drought tolerance mechanism. Several unannotated differentially expressed genes (DEGs) were identified through our previous RNA-Seq study using two hulless barley accessions with contrasting drought tolerance. One of these DEGs, HVU010048.2, showed up-regulated pattern under dehydration stress in both drought tolerant (DT) and drought susceptible (DS) accessions, while its function in drought resistance remains unknown. This new gene was named as HvLRX (light responsive X), because its expression was induced under high light intensity while suppressed under dark.

OBJECTIVE

To provide preliminary bioinformatics prediction, expression pattern, and drought resistance function of this new gene.

METHODS

Bioinformatics analysis of HvLRX were conducted by MEGA, PlantCARE, ProtParam, CELLO et al. The expression pattern of HvLRX under different light intensity, dehydration shock, gradual drought stress, NaCl stress, polyethylene glycol (PEG) 6000 stress and abscisic acid (ABA) treatment was investigated by quantitative reverse transcription-polymerase chain reaction (RT-qPCR). The function of HvLRX was analyzed by virus induced gene silencing (VIGS) in hulless barley and by transgenic method in tobacco.

RESULTS

Full cDNAs of HvLRX were cloned and compared in three hulless barley accessions. Homologues of HvLRX protein in other plants were excavated and their phylogenetic relationship was analyzed. Several light responsive elements (ATC-motif, Box 4, G-box, Sp1, and chs-CMA1a) were identified in its promoter region. Its expression can be promoted under high light intensity, dehydration shock, gradual drought stress, PEG 6000, and NaCl stress, but was almost unchanged in ABA treatment. HvLRX-silenced plants had a higher leaf water loss rate (WLR) and a lower survival rate (SR) compared with controls under dehydration stress. The infected leaves of HvLRX-silenced plants lost their water content quickly and became withered at 10 dpi. The SR of HvLRX overexpressed transgenic tobacco plants was significantly higher than that of wild-type plants. These results indicated HvLRX play a role in drought resistance. Besides, retarded vegetative growth was detected in HvLRX-silenced hulless barley plants, which suggested that this gene is important for plant development.

CONCLUSIONS

This study provided data of bioinformatics, expression pattern, and function of HvLRX. To our knowledge, this is the first report of this new dehydration and light responsive gene.

The gut microbiota metabolite urolithin A inhibits NF-κB activation in LPS stimulated BMDMs

Inflammation is a natural defense process of the innate immune system, associated with the release of proinflammatory cytokines such as interleukin-1β, interleukin-6, interleukin-12 and TNFα; and enzymes including iNOS through the activation and nuclear translocation of NF-κB p65 due to the phosphorylation of IκBα. Regulation of intracellular Ca²⁺ is considered a promising strategy for the prevention of reactive oxygen species (ROS) production and accumulation of DNA double strand breaks (DSBs) that occurs in inflammatory-associated-diseases. Among the metabolites of ellagitannins that are produced in the gut microbiome, urolithin A (UA) has received an increasing attention as a novel candidate with anti-inflammatory and anti-oxidant effects. Here, we investigated the effect of UA on the suppression of pro-inflammatory molecules and NF-κB activation by targeting TLR4 signalling pathway. We also identified the influence of UA on Ca²⁺ entry, ROS production and DSBs availability in murine bone-marrow-derived macrophages challenged with lipopolysaccharides (LPS). We found that UA inhibits IκBα phosphorylation and supresses MAPK and PI3K activation. In addition, UA was able to reduce calcium entry, ROS production and DSBs availability. In conclusion, we suggest that urolithin A is a promising therapeutic agent for treating inflammatory diseases through suppression of NF-κB and preserving DNA through maintaining intracellular calcium and ROS homeostasis.

Identification of valid reference genes for quantitative RT-PCR in Caragana microphylla under salt and drought stresses

Caragana microphylla is a leguminosae plant and grows mainly in semi-arid areas of northwest China and Mongolia. However, the lack of studies on C. microphylla reference genes limits the accurate understanding of the molecular biology mechanisms in this crop under abiotic stresses. In this study, we selected nine candidate genes from salt-treated C. microphylla transcriptome data and evaluated their stability by using geNorm, NormFinder, BestKeeper, and RefFinder in salt and drought conditions. In addition, the relative expressions of Delta 1-pyrroline-5-carboxylate synthase 2 (P5CS2) and Catalase 2 (CAT2) were examined to confirm the stability of the candidate reference genes. As a results, glyceraldehyde-3-phosphate dehydrogenase C2 (GAPC2) and 26S proteasome regulatory subunit (RPN5) were the most stable in both salt and drought treatments. The relative expression of P5CS2 and CAT2 also showed more stable levels in normalization by GAPC2 and RPN5 than the most unstable gene, Ubiquitin 4 (UBQ4). Therefore, it is believed that these candidate reference genes selected and validated in our study could be used to study the molecular biological study of response to salt and drought stress in C. microphylla.

A 2-Cys peroxiredoxin gene from Tamarix hispida improved salt stress tolerance in plants

BACKGROUND

Peroxiredoxins (Prxs) are a large family of antioxidant enzymes that respond to biotic and abiotic stress by decomposing reactive oxygen species (ROS). In this study, the stress tolerance function of the Th2CysPrx gene was further analysed. It lays a foundation for further studies on the salt tolerance molecular mechanism of T. hispida and improved salt tolerance via transgenic plants.

RESULTS

In this study, the stress tolerance function of the Th2CysPrx gene was further analysed. The results of transgenic tobacco showed higher seed germination rates, root lengths, and fresh weight under salt stress than wild-type tobacco. Simultaneously, physiological indicators of transgenic tobacco and T. hispida showed that Th2CysPrx improved the activities of antioxidant enzymes and enhanced ROS removal ability to decrease cellular damage under salt stress. Moreover, Th2CysPrx improved the expression levels of four antioxidant genes (ThGSTZ1, ThGPX, ThSOD and ThPOD).

CONCLUSIONS

Overall, these results suggested that Th2CysPrx enhanced the salt tolerance of the transgenic plants. These findings lay a foundation for further studies on the salt tolerance molecular mechanism of T. hispida and improved salt tolerance via transgenic plants.

Identification of Appropriate Reference Genes for Normalizing miRNA Expression in Citrus Infected by Xanthomonas citri subsp. citri

MicroRNAs (miRNAs) are short noncoding RNA molecules that regulate gene expression at the posttranscriptional level. Reverse transcription-quantitative PCR (RT-qPCR) is one of the most common methods used for quantification of miRNA expression, and the levels of expression are normalized by comparing with reference genes. Thus, the selection of reference genes is critically important for accurate quantification. The present study was intended to identify appropriate miRNA reference genes for normalizing the level of miRNA expression in Citrus sinensis L. Osbeck and Citrus reticulata Blanco infected by Xanthomonas citri subsp. citri, which caused citrus canker disease. Five algorithms (Delta Ct, geNorm, NormFinder, BestKeeper and RefFinder) were used for screening reference genes, and two quantification approaches, poly(A) extension RT-qPCR and stem-loop RT-qPCR, were used to determine the most appropriate method for detecting expression patterns of miRNA. An overall comprehensive ranking output derived from the multi-algorithms showed that poly(A)-tailed miR162-3p/miR472 were the best reference gene combination for miRNA RT-qPCR normalization in citrus canker research. Candidate reference gene expression profiles determined by poly(A) RT-qPCR were more consistent in the two citrus species. To the best of our knowledge, this is the first systematic comparison of two miRNA quantification methods for evaluating reference genes. These results highlight the importance of rigorously assessing candidate reference genes and clarify some contradictory results in miRNA research on citrus.

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

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

Reference genes selection for Calotropis procera under different salt stress conditions

Calotropis procera is a perennial Asian shrub with significant adaptation to adverse climate conditions and poor soils. Given its increased salt and drought stress tolerance, C. procera stands out as a powerful candidate to provide alternative genetic resources for biotechnological approaches. The qPCR (real-time quantitative polymerase chain reaction), widely recognized among the most accurate methods for quantifying gene expression, demands suitable reference genes (RGs) to avoid over- or underestimations of the relative expression and incorrect interpretation. This study aimed at evaluating the stability of ten RGs for normalization of gene expression of root and leaf of C. procera under different salt stress conditions and different collection times. The selected RGs were used on expression analysis of three target genes. Three independent experiments were carried out in greenhouse with young plants: i) Leaf₁₀₀ = leaf samples collected 30 min, 2 h, 8 h and 45 days after NaCl-stress (100 mM NaCl); ii) Root₅₀ and iii) Root₂₀₀ = root samples collected 30 min, 2 h, 8 h and 1day after NaCl-stress (50 and 200 mM NaCl, respectively). Stability rank among the three algorithms used showed high agreement for the four most stable RGs. The four most stable RGs showed high congruence among all combination of collection time, for each software studied, with minor disagreements. CYP23 was the best RG (rank of top four) for all experimental conditions (Leaf₁₀₀, Root_50, and Root₂₀₀). Using appropriated RGs, we validated the relative expression level of three differentially expressed target genes (NAC78, CNBL4, and ND1) in Leaf₁₀₀ and Root₂₀₀ samples. This study provides the first selection of stable reference genes for C. procera under salinity. Our results emphasize the need for caution when evaluating the stability RGs under different amplitude of variable factors.

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

Carex rigescens is an ornamental turfgrass in northern China which has a relatively low maintenance cost and robust tolerance to many adverse environmental conditions, so it could be considered a new material for researching into plant stress resistance. However, suitable reference genes are vacant for obtaining reliable results in quantitative real-time PCR (qRT-PCR) analysis of C. rigescens in adversity research. In this study, we tested the expression stability of nine potential reference genes in leaves and roots under five different abiotic stress conditions, including cold, salt, heat, osmotic and cadmium (Cd). We then selected the best reference genes according to the analysis results calculated by three algorithmic programs (geNorm, NormFinder and BestKeeper) and used the RankAggreg package to merge the outputted data. The results showed that combinations of at least two reference genes should be used for reliable normalization except in heat-treated root samples, which require three reference genes. eIF-4α, GADPH, SAND and PEPKR1 and their combination were found to be the most stably expressed reference genes, while SAM, TUA-α and UPL7 were the three least stable reference genes among most of experimental samples. In addition, five stress-induced genes (Cu-Zn SOD, P5CS, LEA, GST, and APX) were chosen to verify the stability of the selected reference genes in various tissues and under various stress conditions. The results of this study will provide an important fundamental basis both for gene expression verification for transcriptomic and proteomic analyses and for gene expression analysis for future gene function research in C. rigescens.

Selection and validation of reference genes for normalization of qRT-PCR gene expression in wheat (Triticum durum L.) under drought and salt stresses

Eight candidate housekeeping genes were examined as internal controls for normalizing expression analysis of durum wheat (Triticum durum L.) under drought and salinity stress conditions. Quantitative real-time PCR was used to analyse gene expression of multiple stress levels, plant ages (24 and 50 days old), and plant tissues (leaf and root). The algorithms BestKeeper, NormFinder, GeNorm, the delta Ct method and the RefFinder were applied to determine the stability of candidate genes. Under drought stress, the most stable reference genes were glyceraldehyde-3 phosphate, ubiquitin and β-tubulin2, whereas under salinity stress conditions, eukaryotic elongation factor 1-α, glyceraldehyde-3 phosphate and actin were identified as the most stable reference genes. Validation with stress-responsive genes NAC29 and NAC6 demonstrated that the expression level of target genes could be determined reliably with combinations of up to three of the reference genes. This is the first report on reference genes appropriate for quantification of target gene expression in T. durum under drought and salt stresses. Results of this investigation may be applicable to other Triticum species.

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.

Validation of reference genes for accurate normalization of gene expression with quantitative real-time PCR in Haloxylon ammodendron under different abiotic stresses

Haloxylon ammodendron plays an important role in maintaining the structure and function of the entire ecosystem where it grows. No suitable reference genes have been reported in H. ammodendron plants to date. In this study, a total of 8 reference genes (18S, ACT1, ACT7, UBC18, TUA5, GAPDH, EF-1α and UBQ10) were selected from the available trancriptome database, and the expression stability of these 8 candidate genes was validated under different abiotic stress with three different statistical algorithms (geNorm, NormFinder and BestKeeper). The results produced from different models were in agreement with each other essentially: 18S and TUA5 were the most stable genes under drought stress, 18S, the most stable gene under heat stress and mechanical damage, ACT7 and UBC18, stable under salt stress while TUA5 and GAPDH expressed constantly under mechanical damage, and ACT1 expressed steadily under cold conditions. Expression profiles of several stress response genes, including FT-5, FT-9, DREB2A and DREB2C, were further confirmed with various candidate reference genes. None of the candidate genes showed a constant expression among all tested samples. Hence, it's essential to use more than one reference gene in order to guarantee the accuracy of quantitative real-time PCR. The results of this study will contribute to the accuracy and reliability in transcripts quantification, which is of significance to transcription-based studies and applications in this important shrub H. ammodendron.

Selection of relatively exact reference genes for gene expression studies in goosegrass (Eleusine indica) under herbicide stress

Goosegrass (Eleusine indica) is one of the most serious annual grassy weeds worldwide, and its evolved herbicide-resistant populations are more difficult to control. Quantitative real-time PCR (qPCR) is a common technique for investigating the resistance mechanism; however, there is as yet no report on the systematic selection of stable reference genes for goosegrass. This study proposed to test the expression stability of 9 candidate reference genes in goosegrass in different tissues and developmental stages and under stress from three types of herbicide. The results show that for different developmental stages and organs (control), eukaryotic initiation factor 4 A (eIF-4) is the most stable reference gene. Chloroplast acetolactate synthase (ALS) is the most stable reference gene under glyphosate stress. Under glufosinate stress, eIF-4 is the best reference gene. Ubiquitin-conjugating enzyme (UCE) is the most stable reference gene under quizalofop-p-ethyl stress. The gene eIF-4 is the recommended reference gene for goosegrass under the stress of all three herbicides. Moreover, pairwise analysis showed that seven reference genes were sufficient to normalize the gene expression data under three herbicides treatment. This study provides a list of reliable reference genes for transcript normalization in goosegrass, which will facilitate resistance mechanism studies in this weed species.