Abstract
Cabbage (Brassica oleracea) is one of the most important vegetable crops worldwide. qRT-PCR is a sensitive technique for gene expression studies and choosing the appropriate reference gene is essential to obtain reliable results. In the present work, 22 candidate reference genes were evaluated under various experimental conditions, including NaCl, drought stress treatment, temperature treatments (cold and heat) and a set of hormones stress (6-BA, NAA, and ABA) treatments, across a range of tissue types and cultivars. Gene expression data taken from 45 cabbage samples was analyzed using two algorithms, geNorm and NormFinder. Suitable combinations of reference genes for qRT-PCR normalization should be applied according to different experimental conditions. In this study, the genes EF1a, GAPC2 and SAND were verified as the suitable reference genes across all tested samples. Additionally, each experimental condition had a unique set of reference genes best suited to samples within the particular condition. To validate the suitability of the candidate reference genes, the gene expression of WSD1, a gene that may be involved in biosynthesis pathway of wax esters in cabbage, was measured across all 45 samples and normalized using the three best reference gene candidates. WSD1 displayed variation in gene expression across different tissues and cultivars, and exhibited diverse up- or down- regulated expression patterns under various treatments, which indicate that BoWSD1 may play an important role in the response to abiotic stresses in cabbage. Our results provide the foundation for gene expression analysis in Brassica oleracea and other species of Brassica vegetables.
Similar content being viewed by others
Abbreviations
- ABA:
-
Abscisic acid
- bp:
-
Base pair
- NAA:
-
1-Naphthaleneacetic acid
- PEG6000:
-
Polyethylene Glycol 6000
- qRT-PCR:
-
Quantitative real-time polymerase chain reaction
- 6-BA:
-
6-Benzylaminopurine
References
Andersen CL, Jensen JL, Orntoft 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 Res 64:5245–5250
Brulle F, Bernard F, Vandenbulcke F, Cuny D, Dumez S (2014) Identification of suitable qPCR reference genes in leaves of Brassica oleracea under abiotic stresses. Ecotoxicology 23:459–471
Bustin SA (2002) Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol 29:23–29
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55:611–622
Campos MD, Frederico AM, Nothnagel T, Arnholdt-Schmitt B, Cardoso H (2015) Selection of suitable reference genes for reverse transcription quantitative real-time PCR studies on different experimental systems from carrot (Daucus carota L.). Sci Hortic 186:115–123
Castro-Quezada P, Aarrouf J, Claverie M, Favery B, Mugniéry D, Lefebvre V, Caromel B (2013) Identification of reference genes for normalizing RNA expression in potato roots infected with cyst nematodes. Plant Mol Biol Rep 31:936–945
Cheng Y, Bian W, Pang X, Yu J, Ahammed GJ, Zhou G, Wang R, Ruan M, Li Z, Ye Q, Yao Z, Yang Y, Wan H (2017a) Genome-wide identification and evaluation of reference genes for quantitative RT-PCR analysis during tomato fruit development. Front Plant Sci 8:1440
Cheng Y, Pang X, Wan H, Ahammed GJ, Yu J, Yao Z, Ruan M, Ye Q, Li Z, Wang R, Yang Y, Zhou G (2017b) Identification of optimal reference genes for normalization of qPCR analysis during pepper fruit development. Front Plant Sci 8:1128
Delporte M, Legrand G, Hilbert JL, Gagneul D (2015) Selection and validation of reference genes for quantitative real-time PCR analysis of gene expression in Cichorium intybus. Front Plant Sci 6:651
Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GA, Zumla A (2005) The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization. Anal Biochem 344:141–143
Duan M, Wang J, Zhang X, Yang H, Wang H, Qiu Y, Song J, Guo Y, Li X (2017) Identification of optimal reference genes for expression analysis in radish (Raphanus sativus L.) and its relatives based on expression stability. Front Plant Sci 8:1605
Ferradás Y, Rey L, Martínez Ó, Rey M, González MV (2016) Identification and validation of reference genes for accurate normalization of real-time quantitative PCR data in kiwifruit. Plant Physiol Biochem 102:27–36
Fuentes A, Ortiz J, Saavedra N, Salazar LA, Meneses C, Arriagada C (2016) Reference gene selection for quantitative real-time PCR in Solanum lycopersicum L. inoculated with the mycorrhizal fungus Rhizophagus irregularis. Plant Physiol Biochem 101:124–131
Gutierrez L, Mauriat M, Pelloux J, Bellini C, Van Wuytswinkel O (2008) Towards a systematic validation of references in real-time RT-PCR. Plant Cell 20:1734–1735
He Y, Yan H, Hua W, Huang Y, Wang Z (2016) Selection and validation of reference genes for quantitative real-time PCR in Gentiana macrophylla. Front Plant Sci 7:945
Iskandar HM, Simpson RS, Casu RE, Bonnett GD, Maclean DJ, Manners JM (2004) Comparison of reference genes for quantitative real-time polymerase chain reaction analysis of gene expression in sugarcane. Plant Mol Biol Rep 22:325–337
Kim HA, Lim CJ, Kim S, Choe JK, Jo SH, Baek N, Kwon SY (2014) High-throughput sequencing and de novo assembly of Brassica oleracea var. Capitata L. for transcriptome analysis. PLoS ONE 9:92087
Kong Q, Gao L, Cao L, Liu Y, Saba H, Huang Y, Bie Z (2016) Assessment of suitable reference genes for quantitative gene expression studies in melon fruits. Front Plant Sci 7:1178
Lacerda A, Fonseca L, Blawid R, Boiteux L, Ribeiro S, Brasileiro A (2015) Reference gene selection for qPCR analysis in tomato-bipartite begomovirus interaction and validation in additional tomato-virus pathosystems. PLoS ONE 10:e0136820
Li F, Wu X, Lam P, Bird D, Zheng H, Samuels L, Jetter R, Kunst L (2008) Identification of the wax ester synthase/acyl-coenzyme A: diacylglycerol acyltransferase WSD1 required for stem wax ester biosynthesis in Arabidopsis. Plant Physiol 148:97–107
Li W, Qian YQ, Han L, Liu JX, Li ZJ, Ju GS, Sun ZY (2015) Validation of candidate reference genes for gene expression normalization in Buchloe dactyloides using quantitative real-time RT-PCR. Sci Hortic 197:99–106
Li J, Jia HX, Han XJ, Zhang J, Sun P, Lu MZ, Hu JJ (2016a) Selection of reliable reference genes for gene expression analysis under abiotic stresses in the desert biomass willow, Salix psammophila. Front Plant Sci 7:1505
Li MY, Song X, Wang F, Xiong AS (2016b) Suitable reference genes for accurate gene expression analysis in parsley (Petroselinum crispum) for abiotic stresses and hormone stimuli. Front Plant Sci 7:1481
Li MY, Wang F, Jiang Q, Wang GL, Tian C, Xiong AS (2016c) Validation and comparison of reference genes for qPCR normalization of celery (Apium graveolens) at different development stages. Front Plant Sci 7:313
Li J, Han X, Wang C, Qi W, Zhang W, Tang L, Zhao X (2017) Validation of suitable reference genes for qRT-PCR data in Achyranthes bidentata blume under different experimental conditions. Front Plant Sci 8:776
Liu S, Liu Y, Yang X, Tong C, Edwards D, Parkin IA, Zhao M, Ma J, Yu J, Huang S, Wang X, Wang J, Lu K, Fang Z, Bancroft I, Yang TJ, Hu Q, Wang X, Yue Z, Li H, Yang L, Wu J, Zhou Q, Wang W, King GJ, Pires JC, Lu C, Wu Z, Sampath P, Wang Z, Guo H, Pan S, Yang L, Min J, Zhang D, Jin D, Li W, Belcram H, Tu J, Guan M, Qi C, Du D, Li J, Jiang L, Batley J, Sharpe AG, Park BS, Ruperao P, Cheng F, Waminal NE, Huang Y, Dong C, Wang L, Li J, Hu Z, Zhuang M, Huang Y, Huang J, Shi J, Mei D, Liu J, Lee TH, Wang J, Jin H, Li Z, Li X, Zhang J, Xiao L, Zhou Y, Liu Z, Liu X, Qin R, Tang X, Liu W, Wang Y, Zhang Y, Lee J, Kim HH, Denoeud F, Xu X, Liang X, Hua W, Wang X, Wang J, Chalhoub B, Paterson AH (2014) The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat Commun 5:3930
Liu D, Tang J, Liu Z, Dong X, Zhuang M, Zhang Y, Lv H, Sun P, Liu Y, Li Z, Ye Z, Fang Z, Yang L (2017a) Cgl2 plays an essential role in cuticular wax biosynthesis in cabbage (Brassica oleracea L. var. capitata). BMC Plant Biol 17:223
Liu Z, Fang Z, Zhuang M, Zhang Y, Lv H, Liu Y, Li Z, Sun P, Tang J, Liu D, Zhang Z, Yang L (2017b) Fine-Mapping and analysis of cgl1, a gene conferring glossy trait in cabbage (Brassica oleracea L. var. capitata). Front Plant Sci 8:239
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25:402–408
Lv H, Fang Z, Yang L, Zhang Y, Wang Q, Liu Y, Zhuang M, Yang Y, Xie B, Liu B, Liu J, Kang J, Wang X (2014) Mapping and analysis of a novel candidate fusarium wilt resistance gene FOC1 in Brassica oleracea. BMC Genom 15:1094
Mamo S, Gal AB, Bodo S, Dinnyes A (2007) Quantitative evaluation and selection of reference genes in mouse oocytes and embryos cultured in vivo and in vitro. BMC Dev Biol 7:14
Medrano G, Guan P, Barlow-Anacker AJ, Gosain A (2017) Comprehensive selection of reference genes for quantitative RT-PCR analysis of murine extramedullary hematopoiesis during development. PLoS ONE 12:e0181881
Nazari F, Parham A, Maleki AF (2015) GAPDH, β-actin and β2-microglobulin, as three common reference genes, are not reliable for gene expression studies in equine adipose- and marrow-derived mesenchymal stem cells. J Anim Sci 57:1–8
Niu L, Tao YB, Chen MS, Fu Q, Li C, Dong Y, Wang X, He H, Xu ZF (2015a) Selection of reliable reference genes for gene expression studies of a promising oilseed crop, Plukenetia volubilis, by real-time quantitative PCR. Int J Mol Sci 16:12513–12530
Niu X, Qi J, Chen M, Zhang G, Tao A, Fang P, Xu J, Onyedinma SA, Su J (2015b) Reference genes selection for transcript normalization in kenaf (Hibiscus cannabinus L.) under salinity and drought stress. Peer J 3:e1347
Niu X, Qi J, Zhang G, Xu J, Tao A, Fang P, Su J (2015c) Selection of reliable reference genes for quantitative real-time PCR gene expression analysis in Jute (Corchorus capsularis) under stress treatments. Front Plant Sci 6:848
Obrero A, Die JV, Romen B, Gomez P, Nadal S, Gonzalez-Verdejo CL (2011) Selection of reference genes for gene expression studies in zucchini (Cucurbita pepo) using qPCR. J Agric Food Chem 59:402–5411
Park CS, Go YS, Suh MC (2016) Cuticular wax biosynthesis is positively regulated by WRINKLED4, an AP2/ERF-type transcription factor, in Arabidopsis stems. Plant J 88:257–270
Perini P, Pasquali G, Margis-Pinheiro M, de Oliviera PRD, Revers LF (2014) Reference genes for transcriptional analysis of flowering and fruit ripening stages in apple (Malus x domestica Borkh.). Mol Breed 34:829–842
Qi J, Yu S, Zhang F, Shen X, Zhao X, Yu Y, Zhang D (2010) Reference gene selection for real-time quantitative polymerase chain reaction of mrna transcript levels in chinese cabbage (Brassica rapa L. ssp. pekinensis). Plant Mol Biol Rep 28:597–604
Ramakers C, Ruijter JM, Deprez RHL, Moorman AF (2003) Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci Lett 339:62–66
Shivhare R, Lata C (2016) Selection of suitable reference genes for assessing gene expression in pearl millet under different abiotic stresses and their combinations. Sci Rep 6:230360
Silveira ÉD, Alves-Ferreira M, Guimarães LA, da Silva FR, de Campos Carneiro VT (2009) Selection of reference genes for quantitative real-time PCR expression studies in the apomictic and sexual grass Brachiaria brizantha. BMC Plant Biol 9:84
Vandesompele J, Preter KD, Pattyn F, Poppe B, Roy NV, Paepe AD, Speleman F (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3:research0034
Wan H, Zhao Z, Qian C, Sui Y, Malik AA, Chen J (2010) Selection of appropriate reference genes for gene expression studies by quantitative real-time polymerase chain reaction in cucumber. Anal Biochem 399:257–261
Wang Z, Chen Y, Fang H, Shi H, Chen K, Zhang Z, Tan X (2014) Selection of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in Brassica napus under various stress conditions. Mol Genet Genomics 289:1023–1035
Wang C, Cui HM, Huang TH, Liu TK, Hou XL, Li Y (2016) Identification and validation of reference genes for qRT-PCR analysis in non-heading chinese cabbage flowers. Front Plant Sci 7:811
Xiao X, Ma J, Wang J, Wu X, Li P, Yao Y (2014) Validation of suitable reference genes for gene expression analysis in the halophyte Salicornia europaea by real-time quantitative PCR. Front Plant Sci 5:788
Xu Y, Zhu X, Gong Y, Xu L, Wang Y, Liu L (2012) Evaluation of reference genes for gene expression studies in radish (Raphanus sativus L.) using quantitative real-time PCR. Biochem Biophys Res Commun 424:398–403
Xu H, Bao JD, Dai JS, Li Y, Zhu Y (2015a) Genomewide identification of new reference genes for qRT-PCR normalization under high temperature stress in rice endosperm. PLoS ONE 10:e0142015
Xu Y, Li H, Li X, Lin J, Wang Z, Yang Q, Chang Y (2015b) Systematic selection and validation of appropriate reference genes for gene expression studies by quantitative real-time PCR in pear. Acta Physiol Plant 37:40
Xu X, Liu X, Chen S, Li B, Wang X, Fan C, Wang G, Ni H (2016) Selection of relatively exact reference genes for gene expression studies in flixweed (Descurainia sophia) by quantitative real-time polymerase chain reaction. Pestic Biochem Physiol 127:59–66
Ye X, Zhang FM, Tao YH, Song SW, Fang JB (2015) Reference gene selection for quantitative real-time PCR normalization in different cherry genotypes, developmental stages and organs. Sci Hortic 181:182–188
You Y, Zhang L, Li P, Yang C, Ma F (2016) Selection of reliable reference genes for quantitative real-time PCR analysis in plum (Prunus salicina Lindl.) under different postharvest treatments. Sci Hortic 210:285–293
Zhang X, Liu Y, Fang Z, Li Z, Yang L, Zhuang M, Zhang Y, Lv H (2016) Comparative transcriptome analysis between broccoli (Brassica oleracea var. italica) and wild cabbage (Brassica macrocarpa Guss) in response to Plasmodiophora brassicae during different infection stages. Front Plant Sci 7:1929
Zhao X, Zhang X, Guo X, Li S, Han L, Song Z, Wang YY, Li JH, Li MJ (2016) Identification and validation of reference genes for qRT-PCR studies of gene expression in Dioscorea opposita. Biomed Res Int 2016:3089584
Zhuang HH, Fu YP, He W, Wang L, Wei YH (2015) Selection of appropriate reference genes for quantitative real-time PCR in Oxytropis ochrocephala Bunge using transcriptome datasets under abiotic stress treatments. Front Plant Sci 6:475
Acknowledgements
The study was partially supported by the Natural Science Foundation of Jiangsu Province (Nos. BK20191239 and BK20190262) and the National Key Research and Development Program of China (2017YFD0101804).
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
All authors declare that they have no potential conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Fig. S1
Dissociation curve analysis generated for all amplicons (DOC 816 kb)
Fig. S2
Expression levels of candidate reference genes across all samples. A line across the box depicts the median. The box indicates the 25% and 75% percentiles. Whiskers represent the maximum and minimum values (DOC 10980 kb)
Fig. S3
Average expression stability values (M) of the candidate reference genes. Average expression stability values (M) of the candidate reference genes were calculated by the geNorm software in cabbage samples under different experimental conditions, including different cultivars (a) and tissues (b), temperature treatment (c), NaCl treatment (d), PEG6000 treatment (e), hormonal treatment (f), and total (g). The lowest M-value indicates the most stable gene and vice versa (DOC 571 kb)
Fig. S4
Average expression stability values (M) of the candidate reference genes. Average expression stability values (M) of the candidate reference genes were calculated by the geNorm software in cabbage samples under different experimental conditions, including heat treatment (a), cold treatment (b), ABA treatment (c), 6-BA treatment (d), and NAA treatment (e). The lowest M-value indicates the most stable gene and vice versa (DOC 456 kb)
Fig. S5
Determination of the optimal number of reference genes. Every bar represents change in normalization accuracy upon stepwise addition of more reference genes according to the ranking in Fig. 4 and Fig. 5. The pairwise variation (Vn/Vn + 1) was calculated from the normalization factors NFn and NFn + 1, with a recommended cutoff threshold of 0.150 (DOC 64 kb)
Rights and permissions
About this article
Cite this article
Zeng, A., Xu, Y., Song, L. et al. Validation of suitable reference genes for qRT-PCR in cabbage (Brassica oleracea L.) under different abiotic stress experimental conditions. J. Plant Biochem. Biotechnol. 30, 184–195 (2021). https://doi.org/10.1007/s13562-020-00556-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13562-020-00556-x