Search
2024 Volume 4
Article Contents
ARTICLE   Open Access    

Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry

More Information
  • The rapid, reliable, and efficient characteristics of quantitative reverse transcription polymerase chain reaction (qRT-PCR) make it a highly advantageous method for assessing gene expression levels. The identification of stable reference genes is crucial for successful gene expression studies. Cultivated strawberry fruit has been extensively investigated as a model for studying the non-climacteric fruit ripening process. However, more research needs to be conducted on identifying suitable reference genes at different developmental stages of strawberry fruit. We selected the 'Yanli' and 'Chuliandeweidao' cultivars to screen potential reference genes in various tissues, organs, and developmental stages of strawberry fruit. Based on the analysis of high-quality haplotype-resolved genome and transcriptomic FPKM data, FaADPrf1 (ADP-ribosylation factor 1), FaGAPC2 (Glyceraldehyde-3-phosphate dehydrogenase), FaPPC1 (Peptidyl-prolyl cis-trans isomerase 1), and FaEF1-α (Elongation factor 1-alpha) were selected as candidate reference genes, along with the commonly used Fa26S rRNA, for normalization purposes. A qRT-PCR analysis showed 89.21% to 101.51% amplification efficiency for five candidate reference genes, with correlation coefficients (R2) exceeding 0.99. Reference genes' expression stability was assessed using GeNorm, NormFinder, BestKeeper, and Comparative delta-Ct method. RefFinder analysis determined that FaGAPC2 and FaADPrf1 were the most suitable reference genes, considering the results obtained from the abovementioned four methods. The calculated results were validated by studying the expression of FaMYB10, FaUGT1, and FaCHS in different developmental stages of 'Yanli' fruit. This validation confirmed that both FaGAPC2 and the combination of FaGAPC2 and FaADPrf1 could serve as suitable reference genes, with the combination of FaGAPC2 and FaADPrf1 being more suitable than the single FaGAPC2 in certain cases. In summary, we obtained suitable reference genes for research on cultivated strawberry fruit development, which will benefit further study on the ripening of non-climacteric fruits.
  • 加载中
  • Supplemental Table S1 FPKM values and other information of candidate reference genes.
    Supplemental Table S2 The FPKM data of FaMYB10, FaCHS and FaUGT1 in different development stages of 'Yanli' fruit.
    Supplemental Fig. S1 Nucleic acid sequence alignment between alleles of ADPrf1 gene.
    Supplemental Fig. S2 Nucleic acid sequence alignment between alleles of GAPC2 gene.
    Supplemental Fig. S3 Nucleic acid sequence alignment between alleles of PPC1 gene.
    Supplemental Fig. S4 Nucleic acid sequence alignment between alleles of EF1-α gene.
    Supplemental Fig. S5 Validation of the correlation coefficients (R2) for five candidate reference genes.
    Supplemental Fig. S6 Validation of melting curves for five candidate reference genes.
  • [1]

    Park SJ, Huh JW, Kim YH, Lee SR, Kim SH, et al. 2012. Selection of internal reference genes for normalization of quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis in the canine brain and other organs. Molecular Biotechnology 54:47−57

    doi: 10.1007/s12033-012-9543-6

    CrossRef   Google Scholar

    [2]

    Gomez C, Terrier N, Torregrosa L, Vialet S, Fournier-Level A, et al. 2009. Grapevine MATE-type proteins act as vacuolar H+-dependent acylated anthocyanin transporters. Plant Physiology 150:402−15

    doi: 10.1104/pp.109.135624

    CrossRef   Google Scholar

    [3]

    Li K, Xu N, Yang Y, Zhang J, Yin H. 2018. Identification and validation of reference genes for RT-qPCR normalization in Mythimna separata (Lepidoptera: Noctuidae). BioMed Research International 2018:1828253

    doi: 10.1155/2018/1828253

    CrossRef   Google Scholar

    [4]

    Kumar S, Bink MCAM, Volz RK, Bus VGM, Chagné D. 2012. Towards genomic selection in apple (Malus × domestica Borkh.) breeding programmes: prospects, challenges and strategies. Tree Genetics & Genomes 8:1−14

    doi: 10.1007/s11295-011-0425-z

    CrossRef   Google Scholar

    [5]

    Li C, Xu J, Deng Y, Sun H, Li Y. 2019. Selection of reference genes for normalization of cranberry (Vaccinium macrocarpon Ait.) gene expression under different experimental conditions. PLoS ONE 14:e0224798

    doi: 10.1371/journal.pone.0224798

    CrossRef   Google Scholar

    [6]

    Luo M, Gao Z, Li H, Li Q, Zhang C, et al. 2018. Selection of reference genes for miRNA qRT-PCR under abiotic stress in grapevine. Scientific Reports 8:4444

    doi: 10.1038/s41598-018-22743-6

    CrossRef   Google Scholar

    [7]

    Qu R, Miao Y, Cui Y, Cao Y, Zhou Y, et al. 2019. Selection of reference genes for the quantitative real-time PCR normalization of gene expression in Isatis indigotica fortune. BMC Molecular Biology 20:9

    doi: 10.1186/s12867-019-0126-y

    CrossRef   Google Scholar

    [8]

    Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology 3:research0034.1

    doi: 10.1186/gb-2002-3-7-research0034

    CrossRef   Google Scholar

    [9]

    Andersen CL, Jensen JL, Ørntoft TF. 2004. 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 Research 64:5245−50

    doi: 10.1158/0008-5472.CAN-04-0496

    CrossRef   Google Scholar

    [10]

    Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. 2004. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper – Excel-based tool using pair-wise correlations. Biotechnology Letters 26:509−15

    doi: 10.1023/B:BILE.0000019559.84305.47

    CrossRef   Google Scholar

    [11]

    Dong Z, Chen P, Zhang N, Huang S, Zhang H, et al. 2019. Evaluation of reference genes for quantitative real-time PCR analysis of gene expression in Hainan medaka (Oryzias curvinotus). Gene Reports 14:94−99

    doi: 10.1016/j.genrep.2018.11.008

    CrossRef   Google Scholar

    [12]

    Xie F, Xiao P, Chen D, Xu L, Zhang B. 2012. miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Molecular Biology 80:75−84

    doi: 10.1007/s11103-012-9885-2

    CrossRef   Google Scholar

    [13]

    Garg R, Sahoo A, Tyagi AK, Jain M. 2010. Validation of internal control genes for quantitative gene expression studies in chickpea (Cicer arietinum L.). Biochemical and Biophysical Research Communications 396:283−88

    doi: 10.1016/j.bbrc.2010.04.079

    CrossRef   Google Scholar

    [14]

    van de Moosdijk AAA, van Amerongen R. 2016. Identification of reliable reference genes for qRT-PCR studies of the developing mouse mammary gland. Scientific Reports 6:35595

    doi: 10.1038/srep35595

    CrossRef   Google Scholar

    [15]

    Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. 2005. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiology 139:5−17

    doi: 10.1104/pp.105.063743

    CrossRef   Google Scholar

    [16]

    Bai Y, Lv Y, Zeng M, Jia P, Lu H, et al. 2020. Selection of reference genes for normalization of gene expression in Thermobia domestica (Insecta: Zygentoma: Lepismatidae). Genes 12:21

    doi: 10.3390/genes12010021

    CrossRef   Google Scholar

    [17]

    Jain M, Nijhawan A, Tyagi AK, Khurana JP. 2006. Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochemical and Biophysical Research Communications 345:646−51

    doi: 10.1016/j.bbrc.2006.04.140

    CrossRef   Google Scholar

    [18]

    Expósito-Rodríguez M, Borges AA, Borges-Pérez A, Pérez JA. 2008. Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant Biology 8:131

    doi: 10.1186/1471-2229-8-131

    CrossRef   Google Scholar

    [19]

    Reid KE, Olsson N, Schlosser J, Peng F, Lund ST. 2006. An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biology 6:27

    doi: 10.1186/1471-2229-6-27

    CrossRef   Google Scholar

    [20]

    Wan H, Zhao Z, Qian C, Sui Y, Malik A, et al. 2010. Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber. Analytical Biochemistry 399:257−61

    doi: 10.1016/j.ab.2009.12.008

    CrossRef   Google Scholar

    [21]

    Remans T, Smeets K, Opdenakker K, Mathijsen D, Vangronsveld J, et al. 2008. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta 227:1343−49

    doi: 10.1007/s00425-008-0706-4

    CrossRef   Google Scholar

    [22]

    Gutierrez L, Mauriat M, Guénin S, Pelloux J, Lefebvre JF, et al. 2008. The lack of a systematic validation of reference genes: a serious pitfall undervalued in reverse transcription-polymerase chain reaction (RT-PCR) analysis in plants. Plant Biotechnology Journal 6:609−18

    doi: 10.1111/j.1467-7652.2008.00346.x

    CrossRef   Google Scholar

    [23]

    Li X, Gong P, Wang B, Wang C, Li M, et al. 2020. Selection and validation of experimental condition-specific reference genes for qRT-PCR in Metopolophium dirhodum (Walker) (Hemiptera: Aphididae). Scientific Reports 10:21951

    doi: 10.1038/s41598-020-78974-z

    CrossRef   Google Scholar

    [24]

    Erkan M, Wang SY, Wang CY. 2008. Effect of UV treatment on antioxidant capacity, antioxidant enzyme and decay in strawberries fruit. Postharvest Biology and Technology 48:163−71

    doi: 10.1016/j.postharvbio.2007.09.028

    CrossRef   Google Scholar

    [25]

    de L. de O. Pineli L, Moretti CL, dos Santos MS, Campos AB, Brasileiro AV, et al. 2011. Antioxidants and other chemical and physical characteristics of two strawberry cultivars at different ripeness stages. Journal of Food Composition and Analysis 24:11−16

    doi: 10.1016/j.jfca.2010.05.004

    CrossRef   Google Scholar

    [26]

    Zhang Y, Peng X, Liu Y, Li Y, Luo Y, et al. 2018. Evaluation of suitable reference genes for qRT-PCR normalization in strawberry (Fragaria × ananassa) under different experimental conditions. BMC Molecular Biology 19:8

    doi: 10.1186/s12867-018-0109-4

    CrossRef   Google Scholar

    [27]

    Amil-Ruiz F, Garrido-Gala J, Blanco-Portales R, Folta KM, et al. 2013. Identification and validation of reference genes for transcript normalization in strawberry (Fragaria × ananassa) defense responses. PLoS ONE 8:e70603

    doi: 10.1371/journal.pone.0070603

    CrossRef   Google Scholar

    [28]

    Galli V, Borowski JM, Perin EC, da Silva Messias R, Labonde J, et al. 2015. Validation of reference genes for accurate normalization of gene expression for real time-quantitative PCR in strawberry fruits using different cultivars and osmotic stresses. Gene 554:205−14

    doi: 10.1016/j.gene.2014.10.049

    CrossRef   Google Scholar

    [29]

    Jia H, Jiu S, Zhang C, Wang C, Tariq P, et al. 2016. Abscisic acid and sucrose regulate tomato and strawberry fruit ripening through the abscisic acid-stress-ripening transcription factor. Plant Biotechnology Journal 14:2045−65

    doi: 10.1111/pbi.12563

    CrossRef   Google Scholar

    [30]

    Wu J, Xu Z, Zhang Y, Chai L, Yi H, et al. 2014. An integrative analysis of the transcriptome and proteome of the pulp of a spontaneous late-ripening sweet orange mutant and its wild type improves our understanding of fruit ripening in citrus. Journal of Experimental Botany 65:1651−71

    doi: 10.1093/jxb/eru044

    CrossRef   Google Scholar

    [31]

    Wu L, Liu D, Wu J, Zhang R, Qin Z, et al. 2013. Regulation of FLOWERING LOCUS T by a MicroRNA in Brachypodium distachyon. The Plant Cell 25:4363−77

    doi: 10.1105/tpc.113.118620

    CrossRef   Google Scholar

    [32]

    Yi S, Mao J, Zhang X, Li X, Zhang Z, et al. 2022. FveARF2 negatively regulates fruit ripening and quality in strawberry. Frontiers in Plant Science 13:1023739

    doi: 10.3389/fpls.2022.1023739

    CrossRef   Google Scholar

    [33]

    López-Gómez C, Pino-Ángeles A, Órpez-Zafra T, Pinto-Medel MJ, Oliver-Martos B, et al. 2013. Candidate gene study of TRAIL and TRAIL Receptors: association with response to interferon beta therapy in multiple sclerosis patients. PLoS ONE 8:e62540

    doi: 10.1371/journal.pone.0062540

    CrossRef   Google Scholar

    [34]

    Xiao X, Ma J, Wang J, Wu X, Li P, et al. 2014. Validation of suitable reference genes for gene expression analysis in the halophyte Salicornia europaea by real-time quantitative PCR. Frontiers in Plant Science 5:788

    doi: 10.3389/fpls.2014.00788

    CrossRef   Google Scholar

    [35]

    Mao J, Wang Y, Wang B, Li J, Zhang C, et al. 2023. High-quality haplotype-resolved genome assembly of cultivated octoploid strawberry. Horticulture Research 10:uhad002

    doi: 10.1093/hr/uhad002

    CrossRef   Google Scholar

    [36]

    Chen J, Zhou J, Hong Y, Li Z, Chen X, et al. 2021. Genome-wide identification of ubiquitin proteasome subunits as superior reference genes for transcript normalization during receptacle development in strawberry cultivars. BMC Genomics 22:88

    doi: 10.1186/s12864-021-07393-9

    CrossRef   Google Scholar

    [37]

    Guo X, Xie Z, Zhang Y, Wang S. 2021. The FvCYP714C2 gene plays an important role in gibberellin synthesis in the woodland strawberry. Genes & Genomics 43:11−16

    doi: 10.1007/s13258-020-01011-w

    CrossRef   Google Scholar

    [38]

    Li X, Fan J, Gruber J, Guan R, Frentzen M, et al. 2013. Efficient selection and evaluation of transgenic lines of Crambe abyssinica. Frontiers in Plant Science 4:162

    doi: 10.3389/fpls.2013.00162

    CrossRef   Google Scholar

    [39]

    Chen C, Wu J, Hua Q, Tel-Zur N, Xie F, et al. 2019. Identification of reliable reference genes for quantitative real-time PCR normalization in pitaya. Plant Methods 15:70

    doi: 10.1186/s13007-019-0455-3

    CrossRef   Google Scholar

    [40]

    Hao X, Horvath DP, Chao WS, Yang Y, Wang X, et al. 2014. Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plant (Camellia sinensis (L.) O. Kuntze). International Journal of Molecular Sciences 15:22155−72

    doi: 10.3390/ijms151222155

    CrossRef   Google Scholar

    [41]

    Zeng W, Sun Z, Cai Z, Chen H, Lai Z, et al. 2017. Comparative transcriptome analysis of soybean response to bean pyralid larvae. BMC Genomics 18:871

    doi: 10.1186/s12864-017-4256-7

    CrossRef   Google Scholar

    [42]

    Yang H, Liu J, Huang S, Guo T, Deng L, et al. 2014. Selection and evaluation of novel reference genes for quantitative reverse transcription PCR (qRT-PCR) based on genome and transcriptome data in Brassica napus L. Gene 538:113−22

    doi: 10.1016/j.gene.2013.12.057

    CrossRef   Google Scholar

    [43]

    Wei L, Mao W, Jia M, Xing S, Ali U, et al. 2018. FaMYB44.2, a transcriptional repressor, negatively regulates sucrose accumulation in strawberry receptacles through interplay with FaMYB10. Journal of Experimental Botany 69:4805−20

    doi: 10.1093/jxb/ery249

    CrossRef   Google Scholar

    [44]

    Ric-Varas P, Barceló M, Rivera JA, Cerezo S, Matas AJ, et al. 2020. Exploring the use of fruit callus culture as a model system to study color development and cell wall remodeling during strawberry fruit ripening. Plants 9:805

    doi: 10.3390/plants9070805

    CrossRef   Google Scholar

    [45]

    Barry CS, Giovannoni JJ. 2007. Ethylene and fruit ripening. Journal of Plant Growth Regulation 26:143−59

    doi: 10.1007/s00344-007-9002-y

    CrossRef   Google Scholar

    [46]

    Zhang J, Lei Y, Wang B, Li S, Yu S, et al. 2020. The high-quality genome of diploid strawberry (Fragaria nilgerrensis) provides new insights into anthocyanin accumulation. Plant Biotechnology Journal 18:1908−24

    doi: 10.1111/pbi.13351

    CrossRef   Google Scholar

  • Cite this article

    Mao J, Li J, Wang Y, Zhang Z. 2024. Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry. Fruit Research 4: e010 doi: 10.48130/frures-0024-0003
    Mao J, Li J, Wang Y, Zhang Z. 2024. Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry. Fruit Research 4: e010 doi: 10.48130/frures-0024-0003

Figures(6)  /  Tables(2)

Article Metrics

Article views(2373) PDF downloads(289)

Other Articles By Authors

ARTICLE   Open Access    

Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry

Fruit Research  4 Article number: e010  (2024)  |  Cite this article

Abstract: The rapid, reliable, and efficient characteristics of quantitative reverse transcription polymerase chain reaction (qRT-PCR) make it a highly advantageous method for assessing gene expression levels. The identification of stable reference genes is crucial for successful gene expression studies. Cultivated strawberry fruit has been extensively investigated as a model for studying the non-climacteric fruit ripening process. However, more research needs to be conducted on identifying suitable reference genes at different developmental stages of strawberry fruit. We selected the 'Yanli' and 'Chuliandeweidao' cultivars to screen potential reference genes in various tissues, organs, and developmental stages of strawberry fruit. Based on the analysis of high-quality haplotype-resolved genome and transcriptomic FPKM data, FaADPrf1 (ADP-ribosylation factor 1), FaGAPC2 (Glyceraldehyde-3-phosphate dehydrogenase), FaPPC1 (Peptidyl-prolyl cis-trans isomerase 1), and FaEF1-α (Elongation factor 1-alpha) were selected as candidate reference genes, along with the commonly used Fa26S rRNA, for normalization purposes. A qRT-PCR analysis showed 89.21% to 101.51% amplification efficiency for five candidate reference genes, with correlation coefficients (R2) exceeding 0.99. Reference genes' expression stability was assessed using GeNorm, NormFinder, BestKeeper, and Comparative delta-Ct method. RefFinder analysis determined that FaGAPC2 and FaADPrf1 were the most suitable reference genes, considering the results obtained from the abovementioned four methods. The calculated results were validated by studying the expression of FaMYB10, FaUGT1, and FaCHS in different developmental stages of 'Yanli' fruit. This validation confirmed that both FaGAPC2 and the combination of FaGAPC2 and FaADPrf1 could serve as suitable reference genes, with the combination of FaGAPC2 and FaADPrf1 being more suitable than the single FaGAPC2 in certain cases. In summary, we obtained suitable reference genes for research on cultivated strawberry fruit development, which will benefit further study on the ripening of non-climacteric fruits.

    • Quantitative reverse transcription polymerase chain reaction (qRT-PCR) is a reliable and widely used technique in molecular biology research for quantifying gene expression levels[1]. It offers several advantages over conventional PCR, including high sensitivity, specificity, accuracy, and high throughput capacity[2, 3]. The precision of qRT-PCR outcomes relies on numerous elements, such as the stability of reference genes, the cDNA's quality, the efficiency of the cDNA polymerase, and the effectiveness of PCR amplification[4, 5]. The consistency of expression levels in reference genes is particularly critical for ensuring qRT-PCR outcomes' reliability[6, 7]. To identify appropriate reference genes, researchers have devised numerous techniques. GeNorm, NormFinder, and BestKeeper are software tools that were created utilizing Excel 2003, while the Comparative delta cycle threshold (ΔCt) method represents a conventional computational method[811]. Furthermore, RefFinder, an online platform integrating the aforementioned four methodologies, furnishes a ranking system for assessing the stability of candidate gene expression[12].

      Housekeeping genes, which play crucial roles in the fundamental life processes of cells, are commonly utilized as reference genes[13, 14]. These genes exhibit relatively stable expression levels across different tissues, organs, developmental stages, and under various biotic and abiotic stresses[15, 16]. Their expression products are generally essential for maintaining cellular life activities or the cytoskeleton of plant somatic cells[17]. Reference genes such as encoding α, β-tubulin (TUA and TUB), Actin (ACT), Histone B (H2B), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and elongation factor-1-alpha (EF1-α) are commonly selected, along with 18S ribosomal RNA (18S rRNA) and 28S ribosomal RNA (28S rRNA), for the normalization of qRT-PCR in plants[1721]. Nonetheless, prior research has provided evidence suggesting the nonexistence of optimal reference genes that demonstrate consistent expression levels in diverse temporal intervals, tissues, and in the presence of various stress-inducing stimuli. This can be attributed to the inherent relativity of gene expression stability, as the expression patterns of genes may exhibit species-specific or even tissue-specific characteristics, thereby challenging the identification of universally stable reference genes[22]. Exposito-Rodriguez et al. found that the expression of some housekeeping genes in tomato was not stable, even at the 6-leaf stage and 7-leaf stage, and from 1 mm to 8 mm in the bud[18]. Therefore, researchers must select the most suitable reference genes based on the specific characteristics of their research materials. Therefore, it is of utmost importance in scientific investigation to ascertain the stability of gene expression of internal reference genes under novel experimental circumstances and identify suitable reference genes[23].

      Strawberries are widely cultivated and highly favored by consumers worldwide due to their distinctive flavor and rich nutritional content, including minerals, vitamins, and microelements[24, 25]. Some research has been conducted on the selection of internal reference genes in strawberry. Zhang et al. analyzed of the expression stability of seven candidate reference genes in different tissues, various stages of fruit development, and under different light quality and low-temperature conditions in cultivated strawberries. The results revealed that the expression of DBP, HISTH4, ACTIN1, and GAPDH genes demonstrated greater stability[26]. Amil-Ruiz et al. tested the expression stability of 13 candidate reference genes in different cultivated strawberry varieties, ripening and senescent conditions, as well as under SA and JA treatments, revealing that FaRIB413, FaACTIN, FaEF1a, and FaGAPDH2 were suitable reference genes. Stress response constitutes a prominent investigation area within the strawberry research field[27]. Galli et al. investigated the expression stability of seven candidate reference genes under different abiotic stress conditions. The results revealed that the DBP gene exhibited the most stable expression under drought stress, while HISTH4 displayed the most stable expression pattern under osmotic and salt stress. On the other hand, GAPDH and 18S exhibited the least stable expression patterns across all conditions[28].

      Non-climacteric fruits, such as strawberries, represent an ideal model for studying the ripening process, which differs significantly from that of climacteric fruits like apples, peaches, bananas, and pears, in which ethylene plays a dominant role[24]. Previous studies suggested that abscisic acid (ABA) may influence strawberry fruit ripening by regulating cell wall degradation, anthocyanin biosynthesis, and growth, but the underlying mechanisms remain to be elucidated[29, 30]. In previous studies, 18S rRNA and 26S rRNA were used as internal reference genes in qRT-PCR to investigate gene expression during fruit ripening[31, 32]. Nevertheless, there is a lack of comprehensive experimental validation regarding the stability of gene expression levels during different developmental stages of strawberry fruit. The precise quantification of target gene expression in qRT-PCR experiments is contingent upon the careful selection of appropriate reference genes, as any variability in the expression levels of these reference genes may introduce inaccuracies that subsequently impact the dependability and authenticity of subsequent research outcomes[33]. Therefore, it is imperative to screen and identify suitable reference genes that exhibit stable expression across various tissues, organs, and developmental stages of cultivated octoploid strawberry. This will greatly facilitate future investigations in this field[34].

      Our objective was to identify a reference gene with relatively consistent expression throughout all sampled tissues, organs, and developmental stages of cultivated octoploid strawberry fruit. In a previous investigation, we successfully achieved a high-quality haplotype-resolved genome of cultivated octoploid strawberry and subsequently reannotated the second-generation transcriptome data based on this genome[35]. Based on the analysis of high-quality haplotype-resolved genome and the fragments per kilobase of exon model per million mapped fragments (FPKM) data of transcriptome, we selected candidate reference genes, along with the commonly used Fa26S rRNA. Primer3 was used to design primers for qPCR, and the analysis of the melting curve and standard curve showed that the primer design met the primary standard. We employed GeNorm, Normfinder, BestKeeper, and Comparative delta-Ct methods to analyze the Ct values of five candidate reference genes. RefFinder (http://blooge.cn/RefFinder/) calculated that FaGAPC2 and FaADPrf1 were suitable reference genes combined with the above four results. Validation using FaMYB10, FaUGT1 and FaCHS in different development stages of 'Yanli' fruit supported the calculation results of RefFinder. Further, it determined that FaGAPC2 and the combination of FaGAPC2 and FaADPrf1 could serve as suitable reference genes, with the combination of FaGAPC2 and FaADPrf1 being more suitable than the single FaGAPC2.

    • We are committed to finding a reference gene with relatively stable expression in different tissues, organs, and fruit development stages of cultivated strawberry. The candidate reference genes were screened using FPKM data from the next-generation transcriptome and the genes whose expression levels were higher than 1,000 respectively in root, shoot, leaf, and five stages of fruit development (SG = Small Green, BG = Big Green, W = White, TR = Turning Red, R = Red) were selected as candidate reference genes. ADP-ribosylation factor 1 (ADPrf1), glyceraldehyde-3-phosphate dehydrogenase (GAPC2), peptidyl-prolyl cis-trans isomerase 1 (PPC1), and elongation factor 1-alpha (EF1-α) were selected as candidate reference genes (Supplemental Table S1). We first verified the homology between the alleles of candidate reference genes, respectively. The results showed that each candidate reference gene had at least eight alleles, and the DNA sequence similarity between alleles was more than 95% (Supplemental Figs S1S4). In addition, previous studies on strawberry always used the Fa26S rRNA gene as a reference gene, so Fa26S rRNA was also selected as a candidate reference gene[3638].

      Primer3 version 4.1.0 was used to design primers of five candidate reference genes. The lengths of amplification fragments for five candidate reference genes ranged from 132 bp (Fa26S rRNA) to 258 bp (FaPPC1). In addition, except FaPPC1 (89.21%), all five candidate reference genes displayed an amplification efficiency exceeding 90%, while the correlation coefficients (R2) surpassed 0.99 (Table 1, Supplemental Fig. S5). The melting profiles of all potential reference genes displayed a singular peak, validating the specificity of the primer design and the existence of a sole PCR amplification product (Supplemental Fig. S6).

      Table 1.  Primers sequence and amplification characteristics of five candidate reference genes.

      Gene symbolGene namePrimer sequence (5'-3')Amplification length (bp)Amplification efficiency (%)Correlation coefficiency (R2)
      ADPrf1ADP–ribosylation factor 1F: 5'-TGCGAATTCTGATGGTCGGT-3'
      R: 5'-CTCCACAATGGACGGATCTT-3'
      144 bp95.43%0.9984
      GAPC2Glyceraldehyde-3-phosphate dehydrogenaseF: 5'-GAATCAACGGATTCGGAAGA-3'
      R: 5'-ACAATATCGGCACCAACTGA-3'
      231 bp101.51%0.9996
      EF1-αElongation factor 1 - alphaF: 5'-CACATCAACATTGTGGTCAT-3'
      R: 5'-GTCTCAAACTTCCACAAGGC-3'
      187 bp99.69%0.9984
      PPC1Peptidyl-prolyl cis-trans isomerase 1F: 5'-TACAAGGGATCGTCCTTCCA-3'
      R: 5'-ACCCAACCTTCTCGATGTTC-3'
      258 bp89.21%0.9955
      26S rRNA26S ribosomal RNAF: 5'-TAACCGCATCAGGTCTCCAA-3'
      R: 5'-CTCGAGCAGTTCTCCGACAG-3'
      132 bp95.69%0.9996
    • GeNorm software was utilized to calculate the M value, which indicates gene expression stability. Lower M values correspond to higher expression stability. The study showed that FaGAPC2 and FaEF1-α had the lowest M value (0.858), which represented the most stable expression, and Fa26S rRNA had the highest M value (1.272), which represented the least stable expression in 'Yanli'. Meanwhile, in 'Chuliandeweidao', FaGAPC2 and FaADPrf1 (1.612) had the most stable expression, and FaPPC1 (3.547) had the least stable expression (Fig. 1a). The pairwise variation (V) analysis indicated that increasing the number of genes increases the average stability of reference genes. Therefore, the qRT-PCR results obtained from the analysis of two reference genes will be more accurate in both Yanli and Chuliandeweidao (Fig.1b).

      Figure 1. 

      Expression stability value of five candidate reference genes calculated by geNorm in 'Yanli' and 'Chuliandeweidao'. (a) M value of five candidate reference genes. A lower M-value indicates more stable gene expression. (b) Pairwise variation (V) analysis of five candidate reference genes. A lower value indicates a more stable combination number of reference genes.

      NormFinder is a Visual Basic application that allows the calculation of reference gene stability, similar to GeNorm. NormFinder first calculates gene expression stability and then outputs specific numbers; the gene expression stability increases proportionally as the numerical value decreases. FaGAPC2 (0.425 for 'Yanli' and 0.806 for 'Chuliandeweidao') exhibited the lowest stable value, representing the highest expression stability in both varieties. However, the lowest expression stability gene in 'Yanli' and 'Chuliandeweidao' differed. Among the five candidate reference genes, Fa26S rRNA showed the lowest level of stability in terms of expression, with a value of 1.442 in 'Yanli'; however, FaPPC1 (3.455) was the lowest expression stability in 'Chuliandeweidao' (Fig. 2). Considering GeNorm and NormFinder, FaGAPC2 was a more suitable reference gene among the five candidates.

      Figure 2. 

      Expression stability value of five candidate reference genes calculated by NormFinder in 'Yanli' and 'Chuliandeweidao'. A lower value indicates more stable gene expression.

      BestKeeper calculates gene expression stability by combining the coefficient of variation (CV) and standard deviation (SD). The results showed that FaADPrf1 (0.86) had the highest expression stability due to its minimum SD value in 'Yanli', while FaEF1-α (1.52) showed the lowest expression stability (Fig. 3). However, there was little difference in SD value between FaADPrf1 and FaEF1-α. FaGAPC2 (1.00) had the highest expression stability y due to its minimum SD value, and FaPPC1 (3.59) displayed the lowest expression stability in 'Chuliandeweidao' (Fig. 3). The stability of expression determined by BestKeeper for the five candidate reference genes differed from that observed through GeNorm and NormFinder analyses.

      Figure 3. 

      Expression stability value of five candidate reference genes calculated by BestKeeper in 'Yanli' and 'Chuliandeweidao'. A lower value indicates more stable gene expression.

      Comparative ΔCt, another method employed in this study, calculates reference gene stability based on SD values. The lower SD value represented the high expression stability and vice versa. The results demonstrated that FaGAPC2 was the most stable reference gene with an SD value of 1.076 in 'Yanli' and 2.663 in 'Chuliandeweidao'. At the same time, FaPPC1 (1.610 in 'Yanli' and 4.145 in 'Chuliandeweidao') showed the least expression stability (Fig. 4).

      Figure 4. 

      Expression stability value of five candidate reference genes calculated by delta-CT method in 'Yanli' and 'Chuliandeweidao'. A lower value indicates more stable gene expression.

      The expression stability of the five candidate reference genes calculated by different methods was not the same, likely due to different statistical methods. Therefore, we used the online tool RefFinder to analyze all calculation results and get the most appropriate ranking. The lower value of FaGAPC2 (1.32 for 'Yanli' and 1.00 for 'Chuliandeweidao') and FaADPrf1 (1.86 for 'Yanli' and 1.68 for 'Chuliandeweidao') indicated that they are more suitable for reference genes. At the same time, the performance of FaPPC1 (4.00 for 'Yanli' and 5.00 for 'Chuliandeweidao') was poor according to all evaluation systems (Table 2). Based on the results from RefFinder and Pairwise variation, we considered the combination of FaGAPC2 and FaADPrf1 to be suitable as reference genes for qRT-PCR experiments.

      Table 2.  The comprehensive ranking of five candidate reference genes in 'Yanli' and 'Chuliandeweidao' analyzed by RefFinder.

      RankYanliChuliandeweidao
      Gene nameRanking valueGene nameRanking value
      1GAPC21.32GAPC21.00
      2ADPrf11.86ADPrf11.68
      3EF1-α2.5926S3.22
      426S3.98EF1-α3.72
      5PPC14.00PPC15.00
    • To further validate the stability of the selected reference genes, we proceeded to examine the expression of FaMYB10, FaUGT1, and FaCHS genes, which are known for their positive regulatory roles in anthocyanin synthesis (Supplemental Table S2). As the strawberry fruits developed, the expression levels of these three genes exhibited an increasing trend. We utilized the single reference gene FaGAPC2 and the combination of FaGAPC2 and FaADPrf1 to calculate the expression levels of the three genes at different developmental stages of the strawberry fruits. As shown in Fig. 5, regardless of whether a single gene or a combination of genes was used as the reference gene, the expression levels of FaMYB10, FaUGT1, and FaCHS showed an increasing trend during the development of the strawberry fruits. This result suggested that the selected reference genes were suitable for strawberry fruit ripening research.

      Figure 5. 

      Relative expression levels of FaMYB10, FaUGT1, and FaCHS in five different development stages of 'Yanli' fruit. The bars of different colors represent the relative expression levels of the validation genes calculated using different reference genes. (SG = Small Green, BG = Big Green, W = White, TR = Turning Red, R = Red).

    • qRT-PCR is widely regarded as the most robust and efficient technique for detecting gene expression levels and patterns in plants[39]. Researchers widely use it because of its specificity, high sensitivity, and simple operation. However, the accuracy of qRT-PCR is influenced by various factors, including RNA quality and integrity, instrument fluorescence detection sensitivity, and appropriate selection of reference genes[40, 41]. In this investigation, reference genes for different tissues, organs, and developmental stages of cultivated octoploid strawberry were identified using two cultivars, 'Yanli' and 'Chuliandeweidao', as experimental materials. Subsequently, stable reference gene expressions were obtained for subsequent qRT-PCR experiments.

      The amalgamation of transcriptome databases, facilitated by advancements in sequencing technology, has emerged as a mature approach for reference gene screening[42]. A previous study showed that the expression stability of reference genes UXS3, SAP5, and ARFA1E mined from transcriptome data was better than that of the traditional reference gene Actin7 in Brassica napus[43]. This study selected five housekeeping genes as candidate reference genes combined with the second-generation transcriptome data and high-quality haplotype-resolved genome. The expression stability of five candidates in different tissues and organs and five fruit development stages were evaluated by qRT-PCR and calculated by geNorm, NormFinder, BestKeeper, and Comparative delta-Ct. Owing to differences in operational logic and statistical methodologies, the ranking of candidate reference gene expression stability varied slightly across different analysis tools. For example, the analysis results of geNorm, NormFinder, and Comparative delta-Ct showed that the expression stability of FaGAPC2 was the best. In contrast, BestKeeper analysis showed that FaADPrf1 had the most stable expression level in 'Yanli'. This phenomenon also occurred in the study of reference genes in poplar (Populus deltoides), strawberry (Fragaria vesca), apple (Malus domestic), and other plants. RefFinder, a widely utilized program for comprehensive stability analysis of candidate reference genes, was employed for reference gene screening. To conduct a comprehensive evaluation of reference gene stability, RefFinder was utilized to thoroughly assess the stability of five selected reference genes. The expression stability ranking from high to low was FaGAPC2 (1.32) > FaADPrf1 (1.86) > FaEF1-α (2.59) > Fa26S rRNA (3.98) > FAPPC1 (4.00) in 'Yanli' and FaGAPC2 (1.00) > FaADPrf1 (1.68) > Fa26S rRNA (3.22) > FaEF1-α (3.72) > FaPPC1 (5.00) in 'Chuliandeweidao', which demonstrated that FaGAPC2 and FaADPrf1 were suitable reference gene for qRT-PCR in cultivated octoploid strawberry. Previous studies have investigated suitable reference genes under different stress conditions, but the ideal reference genes for distinct stress were not the same. While GAPDH demonstrates stable expression under low-temperature conditions, it is not a suitable reference gene under drought and salt stress. Additionally, the suitable reference genes under these two stress conditions are also different[28]. Therefore, further research is needed to select suitable reference genes under different stress conditions.

      FaMYB10, a member of the R2R3-MYB transcription factor family, exerts a pivotal role in anthocyanin biosynthesis in strawberry fruit[44, 45]. Previous studies and transcriptome data revealed that FaMYB10 predominantly exhibits expression during the TR and R stages of strawberry fruit development, thereby verifying the reliability of the selected reference genes. Our study demonstrated that the expression trend of FaMYB10 calculated using either FaGAPC2 or the combination of FaGAPC2 and FaADPrf1 as reference genes aligned with the FPKM values from the transcriptome. According to the FPKM values from the transcriptome, the expression level of FaUGT1 gradually increased during fruit development and reached its highest expression level at the R stage. However, when using FaGAPC2 as a single reference gene to calculate the expression trend of FaUGT1, the expression level was shown to be higher at the TR than at the R stage, which is inconsistent with the FPKM values. Nevertheless, when we used the combination of FaGAPC2 and FaADPrf1 as the reference genes to calculate the expression trend, it aligned with the FPKM values. This suggests that the combination of FaGAPC2 and FaADPrf1 was more suitable than FaGAPC2 alone when calculating the expression level of FaUGT1 gene. The Fa26S rRNA was always used to study the strawberry gene expression pattern. The Ct value analysis demonstrated that the expression level of the Fa26S rRNA gene was the highest among the five candidate genes in different tissues and organs and five development stages of cultivated strawberry. However, the expression pattern of Fa26S rRNA, with Ct values ranging from 7.99 to 13.28, displayed instability, highlighting its unsuitability as a reference gene for normalizing target gene expression across various tissues, organs, and fruit development stages (Fig. 6).

      Figure 6. 

      Boxplot analysis of the expression profiles of FaADPrf1, FaGAPC2, and Fa26S rRNA in root, crown stem, leaf, flower, and five fruit developmental stages. Solid dots represent the expression cycle threshold (CT) values of candidate reference genes in different organs. The line across the box represents the median. The boxes represent the 25/75 percentiles.

    • FaGAPC2 and FaADPrf1 are suitable reference genes for qRT-PCR.

    • The cultivated strawberry (Fragaria × ananassa) cultivars 'Yanli' and 'Chuliandeweidao' were grown in the solar greenhouse of Shenyang Agriculture University (China). Root, crown stem, leaf, flower, and fruits at five different development stages (Small Green = SG, Big Green = BG, White = W, Turning Red = TR, and Red = R) were collected in June 2023. All materials were stored at −80°C.

    • We used the CTAB method to extract the total RNA of the materials. The brief steps were: The samples were ground into powder form and transferred into a 1.5 ml RNase-free centrifuge tube. After adding 588 μl CTAB extraction solution and 12 μl β-mercaptoethanol, the samples were put into a 65 °C water bath for 30 min and shaken violently for 1 min every 5 min. Samples were shaken violently for 5 min after adding 600 μl chloroform/isoamyl alcohol (volume ratio = 24:1) and then centrifuged at 12,000 rpm at 4 °C for 10 min. Four hundred μl supernatant was extracted before adding 400 μl chloroform/isoamyl alcohol (volume ratio = 24:1), shaken violently for 5 min, then centrifuged at 12,000 rpm at 4 °C for 10 min. After absorbing 300 μl supernatant, adding 75 μl of 10 M LiCl2, and precipitating RNA overnight at −20 °C, the total RNA was subjected to two rounds of cleansing using absolute ethanol and subsequently solubilized in DEPC-treated water. NanoDrop 2000 and 1% agarose gel electrophoresis were utilized to assess the purity and integrity of the RNA samples.

      Complementary DNA (cDNA) was reverse transcribed using the PrimeScript™ RT reagent Kit with gDNA Eraser (TaKaRa, Dalian, China) following the manufacturer's protocol.

    • The qPCR was performed on QuantStudioTM 6 Flex Real-Time PCR System (Applied Biosystems). A total volume of 10 μl containing 0.5 μl cDNA, 1 μl gene-specific primers, 3.5 μl ddH2O, and 5 μl UltraSYBR Mixture was mixed and carried out qPCR using UltraSYBR Mixture (CWBio, Beijing, China). The PCR program was set based on the description provided by Zhang et al[46]. The relative mRNA levels were determined by employing the 2−ΔΔCᴛ approach. Each sample was examined in triplicate with three biological replicates. Primer sequences of candidate reference genes for qPCR are listed in Table 1 and FaMYB10 for qPCR are as follows: forward (5'-ACAGATGCAGGAAGAGCTGT-3') and reverse (5'-GTTCTTCCTGGCAATCGTCC-3'). Primer sequences of FaCHS are as follows: forward (5’-TCAACGGCCCAAACTATCCT-3’) and reverse (5’-TTAGCCTCAACCTGGTCCAG-3’). Primer sequences of FaUGT1 are as follows: forward (5’-CAGTAACAAGACCATCGCCG-3’) and reverse (5’-GAGTTCCAACCGCAATGTGT-3’). The design of all candidate genes and validation genes quantitative PCR primers was based on the high-quality haploid genome sequence of 'Yanli' (www.rosaceae.org/Analysis/14723107).

    • GeNorm, NormFinder, BestKeeper, and Comparative delta-Ct were used to calculate the expression stability of five candidate reference genes. The cycle threshold (Ct) value was first converted to the appropriate format and then used to analyze geNorm, Comparative delta-Ct, and NormFinder. BestKeeper is Excel spreadsheet software and can directly input Ct value for calculation. RefFinder is an online tool that can synthesize the above four results and rank candidate gene expression stability.

    • The authors confirm contribution to the paper as follows: carry out experiments and data analysis: Mao J; analyzed the transcriptome data analysis: Li J; helped ensure the completion of the experiment: Wang Y; draft manuscript preparation: Zhang Z. All authors reviewed the results and approved the final version of the manuscript.

    • The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

      • This work was supported by grants from the National Natural Science Foundation of China (31872072, 32130092).

      • The authors declare that they have no conflict of interest. Zhihong Zhang is the Editorial Board member of Fruit Research who was blinded from reviewing or making decisions on the manuscript. The article was subject to the journal's standard procedures, with peer-review handled independently of this Editorial Board member and the research groups.

      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (6)  Table (2) References (46)
  • About this article
    Cite this article
    Mao J, Li J, Wang Y, Zhang Z. 2024. Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry. Fruit Research 4: e010 doi: 10.48130/frures-0024-0003
    Mao J, Li J, Wang Y, Zhang Z. 2024. Selection and validation of reference genes for qRT-PCR in cultivated octoploid strawberry. Fruit Research 4: e010 doi: 10.48130/frures-0024-0003

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return