Selection of reliable reference genes for quantitative RT-PCR in garlic under salt stress

Pungency related gene network in Allium sativum L., response to sulfur treatments

Pungency of garlic (Allium sativum L.) is generated from breakdown of the alk(en)yl cysteine sulphoxide (CSO), alliin and its subsequent breakdown to allicin under the activity of alliinase (All). Based on recent evidence, two other important genes including Sulfite reductase (SiR) and Superoxide dismutase (SOD) are thought to be related to sulfur metabolism. These three gene functions are in sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. In this research, the order and priority of simultaneous expression of three genes including All, SiR and SOD were measured on some garlic ecotypes of Iran, collected from Zanjan, Hamedan and Gilan, provinces under sulfur concentrations (0, 6, 12, 24 and 60 g/ per experimental unit: pot) using real-time quantitative PCR (RT-qPCR) analysis. For understanding the network interactions between studied genes and other related genes, in silico gene network analysis was constructed to investigate various mechanisms underlying stimulation of A. sativum L. to cope with imposed sulfur. Complicated network including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-networks to have a deeper insight. Analysis of q-RT-PCR data revealed the highest expression in All and SiR genes respectively. To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. Complex sub-network exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed their involvement in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic. The findings of this research can be used to extend further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency for health benefits and industry demands.

Real-time expression and in silico characterization of pea genes involved in salt and water-deficit stress

BACKGROUND

To tolerate salt and water-deficit stress, the plant adapts to the adverse environment by regulating its metabolism and expressing certain stress-induced metabolic pathways. This research analyzed the relative expression of four pea genes (P5CR, PAL1, SOD, and POX) in three pea varieties (Climax, Green grass, and Meteor) under different levels of salt and water-deficit stress.

METHODS AND RESULTS

The experiments on salt stress and water-deficit stress were carried out within greenhouse settings under controlled environment. The saturation percentage was employed to create artificial salinity conditions: Control without NaCl treatment, Treatment 1: 50 mM NaCl treatment, Treatment 2: 75 mM NaCl treatment, and Treatment 3: 100 mM NaCl treatment. Field capacity (FC) was used for the development of artificial water-deficit treatments in the pots, i.e., Treatment 1 (Control; water application 100% of FC), Treatment 2 (water application 75% of FC), and Treatment 3 (water application 50% of FC). Pea genes involved in biosynthetic pathways of proline, flavonoids, and enzymatic antioxidant enzymes including P5CR, PAL1, SOD, and POX were selected based on literature. Quantitative real-time PCR using cDNA as a template was used to analyze the gene expression. Pea genes were analyzed for phylogenetic analysis in closely related crops having similarity percent identity 80% and above. In silico characterization of selected proteins including the family classification was done by the NCBI CDD and INTERPRO online servers. Results from RT-qPCR analysis showed increased expression of P5CR, PAL1, and POX genes, while SOD gene expression decreased under both stresses. Climax exhibited superior stress tolerance with elevated expression of P5CR and PAL1, while Meteor showed better tolerance through increased POX expression. Phylogenetic analysis revealed common ancestry with other species like chickpea, red clover, mung bean, and barrel clover, suggesting the cross relationship among these plant species. Conserved domain analysis of respective proteins revealed that these proteins contain PLNO 2688, PLN02457, Cu-Zn Superoxide dismutase, and secretory peroxidase conserved domains. Furthermore, protein family classification indicated that the oxidation-reduction process is the most common chemical process involved in these stresses given to pea plant which validates the relationship of these proteins.

CONCLUSIONS

Salt and water-deficit stresses trigger distinct metabolic pathways, leading to the up-regulation of specific genes and the synthesis of corresponding proteins. These findings further emphasize the conservation of stress-tolerance-related genes and proteins across various plant species. This knowledge enhances our understanding of plant adaptation to stress and offers opportunities for developing strategies to improve stress resilience in crops, thereby addressing global food security challenges.

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

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

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

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

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.

Screening and Verification of Reference Genes for Analysis of Gene Expression in Garlic (Allium sativum L.) under Cold and Drought Stress

The principal objective of this study was to screen and verify reference genes appropriate for gene expression evaluation during plant growth and development under distinct growth conditions. Nine candidate reference genes were screened based on garlic transcriptome sequence data. RT-qPCR was used to detect the expression levels of the aforementioned reference genes in specific tissues under drought and cold stress. Then, geNorm, NormFinder, BestKeeper, and ReFinder were used to consider the consistency of the expression levels of candidate reference genes. Finally, the stress-responsive gene expression of ascorbate peroxidase (APX) was quantitatively evaluated to confirm the chosen reference genes. Our results indicated that there were variations in the abundance and stability of nine reference gene transcripts underneath cold and drought stress, among which ACT and UBC-E2 had the highest transcript abundance, and 18S rRNA and HIS3 had the lowest transcript abundance. UBC and UBC-E2 were the most stably expressed genes throughout all samples; UBC and UBC-E2 were the most stably expressed genes during cold stress, and ACT and UBC were the most stably expressed genes under drought stress. The most stably expressed genes in roots, pseudostems, leaves, and cloves were EF1, ACT, HIS3, UBC, and UBC-E2, respectively, while GAPDH was the most unstable gene during drought and cold stress conditions and in exclusive tissues. Taking the steady reference genes UBC-E2, UBC, and ACT as references during drought and cold stress, the reliability of the expression levels was further demonstrated by detecting the expression of AsAPX. Our work thereby offers a theoretical reference for the evaluation of gene expression in garlic in various tissues and under stress conditions.

Transcriptome-wide identification of NAC (no apical meristem/Arabidopsis transcription activation factor/cup-shaped cotyledon) transcription factors potentially involved in salt stress response in garlic

Soil salinity has been an increasing problem worldwide endangering crop production and human food security. It is an ideal strategy to excavate stress resistant genes and develop salt tolerant crops. NAC (no apical meristem/Arabidopsis transcription activation factor/cup-shaped cotyledon) transcription factors have been demonstrated to be involved in salt stress response. However, relevant studies have not been observed in garlic, an important vegetable consumed in the world. In this study, a total of 46 AsNAC genes encoding NAC proteins were identified in garlic plant by transcriptome data. Phylogenetic analysis showed that the examined AsNAC proteins were clustered into 14 subgroups. Motif discovery revealed that the conserved domain region was mainly composed of five conserved subdomains. Most of the genes selected could be induced by salt stress in different tissues, indicating a potential role in salt stress response. Further studies may focus on the molecular mechanisms of the AsNAC genes in salt stress response. The results of the current work provided valuable resources for researchers aimed at developing salt tolerant crops.

A grobacterium-mediated genetic transformation and cloning of candidate reference genes in suspension cells of Artemisia pallens Wall. ex DC

A reliable and stable Agrobacterium-mediated genetic transformation system for Artemisia pallens has been developed using cell suspension cultures derived from cotyledon explants. Cotyledon, attached cotyledon, and compound leaves were found to be suitable for the induction of callus among five different types of explants tested. The yellow friable callus derived from attached cotyledon was used to initiate suspension cultures in Suspension Culture Medium (SCM) which was supplemented with 2,4-dichlorophenoxyacetic acid (2,4-D) at 2.0 mg L-1 and in combination with different concentrations of Zeatin (ZEA) at 0.25 mg L-1. Two different shock treatments, cold shock (at 4 ℃) for 20 min and heat shock (at 45 ℃) treatment for 5 min, heat shock treatment increased the transformation efficiency. The supplementation of Pluronic F-68 (0.05%) significantly enhanced the transformation efficiency of suspension cultures, whereas Silwet L-77 (0.05%) leads to more browning of the cells and reduced the transformation efficiency. The maximum GUS intensity was recorded with an optimal intensity of blue spots in the transformed cells. The highest GUS fluorometric activity measured was 879.4 ± 113.7 nmol 4MU/mg/min in transformed cell suspension cultures. The hygromycin-resistant calli showed intense blue color in GUS histochemical assay. The transgene integration into the plant genome was confirmed by polymerase chain reaction (PCR) using uidA specific primers in six hygromycin-resistant cell lines. The partial coding sequence of three candidate reference genes, i.e., ADP-ribosylation factor (Arf), β-actin (Act), and ubiquitin (Ubi), and carotenoid biosynthesis pathway gene, i.e., Phytoene desaturase (Pds) were cloned, sequenced, and submitted to NCBI for the first time. The quantitative mRNA expression of the transgene (uidA) and internal ApPds gene were evaluated in transgenic callus lines. The present Agrobacterium-mediated genetic transformation protocol could help in better understanding of the metabolic pathways of this medicinally important plant and its genetic improvement.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03251-x.

Identification of Early Salt-Stress-Responsive Proteins in In Vitro Prunus Cultured Excised Roots

Fruit-tree rootstock selection is a challenge under a scenario of growing environmental stresses in which the soil and climate are greatly affected. Salinization is an increasing global process that severely affects soil fertility. The selection of rootstocks with the ability to tolerate salt stress is essential. Excised root cultures may be an excellent experimental approach to study stress physiology and a predictive tool to assess possible tolerance. In this study, we show how protein changes in response to salt stress evaluated in excised root cultures of Prunus cerasus (moderate salt-sensitive cultivar) could be representative of these changes in the roots of whole plants. The 2D electrophoresis of root extracts and subsequent spot identification by MALDI-TOF/TOF-MS show 16 relevant proteins differentially expressed in roots as a response to 60 mM NaCl. Cytoplasmic isozyme fructose 1,6-bisphosphate aldolase shows relevant changes in its relative presence of isoforms as a response to saline stress, while the total level of enzymes remains similar. Ferredoxin-NADP⁺ reductase increases as a response to salinity, even though the measured activity is not significantly different. The observed changes are congruent with previous proteomic studies on the roots of whole plants that are involved in protection mechanisms against salt stress.

Transcriptome Characterization of the Roles of Abscisic Acid and Calcium Signaling during Water Deficit in Garlic

Garlic (Allium sativum L.) is one of the most important vegetable crops, and breeding drought-tolerant varieties is a vital research goal. However, the underlying molecular mechanisms in response to drought stress in garlic are still limited. In this study, garlic seedlings were subjected to 15% PEG6000 for 0, 1, 4, and 12 h, respectively, to simulate drought stress. Changes of transcriptomes as a result of drought stress in garlic leaves were determined by de novo assembly using the Illumina platform. In total, 96,712 unigenes and 11,936 differentially expressed genes (DEGs) were identified in the presence of drought conditions. Transcriptome profiling revealed that the DEGs were mainly enriched in the biosynthesis of secondary metabolites, MAPK signaling pathway, starch and sucrose metabolism, phenylpropanoid biosynthesis, and plant hormone signal transduction. Genes involved in abscisic acid and calcium signaling were further investigated and discussed. Our results indicated that a coordinated interplay between abscisic acid and calcium is required for drought-induced response in garlic.

Determination of reliable reference genes for gene expression studies in Chinese chive (Allium tuberosum) based on the transcriptome profiling

Chinese chive (Allium tuberosum) is widely cultivated around the world for its unique flavor, nutrient, and medicinal values, yet its molecular mechanism on flavor formation and other metabolic pathways remains intangible. The elucidation of these complex processes begins with investigating the expression of the genes of interest, however the appropriate reference genes (RGs) for normalizing the gene expression are still unavailable in A. tuberosum. To fill this lacuna, transcriptome-wide screening was undertaken to identify the most stable genes according to the analysis of their FPKM values. The expression stability of the RGs was further evaluated using geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The comprehensive analysis showed that GLY1 and SKP1, instead of two traditionally used RGs (eIF1α and ACT2), were the most stable genes across diverse A. tuberosum tissues, indicating the necessity to carefully validate the stability of RGs prior to their use for normalizations. As indicated by geNorm, the normalizations with at least two RGs could give more accurate results. qRT-PCR experiments were conducted with randomly selected genes, demonstrating that normalization with a combination of GLY1 and SKP1 resulted in reliable normalization results. Our finding represents the first attempt toward establishing a standardized qRT-PCR analysis in this economically important vegetable.

Validation of Reference Genes for Quantitative PCR in Johnsongrass (Sorghum halepense L.) under Glyphosate Stress

Weeds are one of the main causes of the decrease in crop yields, with Johnsongrass (Sorghum halepense L.) being one of the most significant. Weeds can be controlled by herbicides, but some have developed resistance. Quantitative PCR is the technique of choice for studying gene expression related to herbicide resistance because of its high sensitivity and specificity, although its quantitative accuracy is highly dependent on the stability of the reference genes. Thus, in this study we evaluated the stability of different reference genes of glyphosate-resistant S. halepense. Nine genes frequently used as reference genes were selected: MDH, ADP, PP2A, EIF4α, ACT, ARI8, DnaJ, Hsp70, and ALS1, and their expression analyzed in susceptible and resistant biotypes at 0, 24 and 72 h post-application of glyphosate. The stability was analyzed with the geNorm, NormFinder, and BestKeeper software programs and using the ΔCt method. RefFinder was used to generate a comprehensive stability ranking. The results showed that PP2A and ARI8 were the most stable genes under the test conditions. EPSPS expression was also verified against the best two and the worst two reference genes. This study provides useful information for gene expression analysis under glyphosate stress and will facilitate resistance mechanism studies in this weed species.

Current Practices for Reference Gene Selection in RT-qPCR of Aspergillus: Outlook and Recommendations for the Future

Aspergillus is a genus of filamentous fungi with vast geographic and ecological distributions. Species within this genus are clinically, agriculturally and biotechnologically relevant, leading to increasing interest in elucidating gene expression dynamics of key metabolic and physiological processes. Reverse-transcription quantitative Polymerase Chain Reaction (RT-qPCR) is a sensitive and specific method of quantifying gene expression. A crucial step for comparing RT-qPCR results between strains and experimental conditions is normalisation to experimentally validated reference gene(s). In this review, we provide a critical analysis of current reference gene selection and validation practices for RT-qPCR gene expression analyses of Aspergillus. Of 90 primary research articles obtained through our PubMed query, 17 experimentally validated the reference gene(s) used. Twenty reference genes were used across the 90 studies, with beta-tubulin being the most used reference gene, followed by actin, 18S rRNA and glyceraldehyde 3-phosphate dehydrogenase. Sixteen of the 90 studies used multiple reference genes for normalisation. Failing to experimentally validate the stability of reference genes can lead to conflicting results, as was the case for four studies. Overall, our review highlights the need to experimentally validate reference genes in RT-qPCR studies of Aspergillus.

Molecular Cloning and Functional Characterization of Heat Stress-Responsive Superoxide Dismutases in Garlic (Allium sativum L.)

Superoxide dismutases (SODs) are key antioxidant enzymes that can detoxify the superoxide radicals generated by various stresses. Although various plant SODs have been suggested to improve stress tolerance, SODs in garlic, an economically important vegetable grown worldwide, remain relatively unknown. In this study, we found that heat stress strongly induced the activities of Cu/ZnSODs, FeSODs, and MnSODs in garlic leaves. In addition, we cloned four garlic SODs (AsSODs) and suggest that heat stress-increased SOD activity was reflected at least by the induction of these AsSODs. The results of the agro-infiltration assay suggested that the cloned AsSODs encoded functional SOD enzymes belonging to the Cu/ZnSOD and MnSOD families. As a first step toward understanding the enzymatic antioxidant system in garlic plants, our results provide a solid foundation for an in-depth analysis of the physiological functions of the AsSOD family.

Selection and validation of appropriate reference genes for RT-qPCR analysis of flowering stages and different genotypes of Iris germanica L

Iris germanica L. is a perennial herbaceous plant that has been widely cultivated worldwide and is popular for its elegant and vibrantly colorful flowers. Selection of appropriate reference genes is the prerequisite for accurate normalization of target gene expression by quantitative real-time PCR. However, to date, the most suitable reference genes for flowering stages have not been elucidated in I. germanica. In this study, eight candidate reference genes were examined for the normalization of RT-qPCR in three I. germanica cultivars, and their stability were evaluated by four different algorithms (GeNorm, NormFinder, BestKeeper, and Ref-finder). The results revealed that IgUBC and IgGAPDH were the most stable reference genes in '00246' and 'Elizabeth', and IgTUB and IgUBC showed stable expression in '2010200'. IgUBC and IgGAPDH were the most stable in all samples, while IgUBQ showed the least stability. Finally, to validate the reliability of the selected reference genes, the expression patterns of IgFT (Flowering Locus T gene) was analyzed and emphasized the importance of appropriate reference gene selection. This work presented the first systematic study of reference genes selection during flower bud development and provided guidance to research of the molecular mechanisms of flowering stages in I. germanica.

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

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

Identification and evaluation of reference genes for quantitative real-time PCR analysis in Passiflora edulis under stem rot condition

Passion fruit (Passiflora edulis), an important tropical and subtropical fruit, has a high edible and medicinal value. Stem rot disease is one of the most important diseases of passion fruit. An effective way for control and prevention of this disease is to identify the genes associated with resistance to this disease. Quantitative real-time PCR (RT-qPCR) has mainly been widely applied to detect gene expression because of its simplicity, fastness, low cost and high sensitivity. One of the requirements for RT-qPCR is the availability of suitable reference genes for normalization of gene expression. However, currently, no Passiflora edulis reference genes have been identified andthus it has hindered the gene expression studies in this plant. The present study aimed to address this issue. We analyzed sixteen candidate reference genes, including nine common (GAPDH, UBQ, ACT1, ACT2, EF-1α-1, EF-1α-2, TUA, NADP, and GBP) and seven novel genes (C13615, C24590, C27182, C10445, C21209, C22199, and C22526), in different tissues (stem, leaf, flower and fruit) of two accessions under stem rot condition. We calculated the expression stability in twenty-four samples using the ΔCt, GeNorm, NormFinder, BestKeeper and RefFinder. The results showed that both C21209 and EF-1α-2 were sufficient to normalize gene expression under stem rot, whereas the commonly used reference genes, GAPDH and UBQ, were the least stable ones. The expression patterns of PeUFC under stem rot condition normalized by stable and unstable reference genes indicated the suitability of using the optimal reference genes. To our knowledge, this is the first systematic study of reference genes in Passiflora edulis, which identified a number of reliable reference genes suitable for gene expression studies in Passiflora edulis by RT-qPCR.