Evaluation And Validation of Reliable Reference Genes For Quantitative Real-Time PCR Analysis of The Gene Expression In Macadamia Integrifolia

Abstract

Background: Macadamia integrifolia, a new economically important crop, the kernel oil is rich in bioactive compound and monounsaturated fatty acid. Gene expression analysis of qRT-PCR is beneficial to understand the complex regulatory networks of macadamia.

Results: In this study, the expression stability of 11 traditional housekeeping genes including α-tubulin (TUBa), β–tubulin (TUBb), malate dehydrogenase (MDH), 18S ribosome RNA (18S), glyceraldehyde-3- phosphate dehydrogenase (GAPDH), α-elongation factor 1 (EF1a), β- elongation factor 1 (EF1b), ubiquitin (UBQ), ubiquitin-conjugating enzyme (UBC), cyclophilin (CYP) and actin (ACT) were accessed by qRT-PCR in macadamia seedlings under different experimental conditions and tissues. The expression stability of the 11 reference genes was evaluated by the online tool RefFinder, which include ΔCt, geNorm, NormFinder, BestKeeper four commonly software, and then determinated a comprehensive expression stability ranking by integrating above four ranking results based on the geometric mean. Our results show that ACT was the best stable genes for all samples, cold stress, NaCl sress, PEG stress, ABA treatment, MeJA treatment, stem and leaf tissue samples; EF1b is the most stable gene in GA treatment and heat stress samples; UBC and CYP were respectively ranked top in ethylene treatment and root tissue samples. Finally, the reliability of these results was further validated with a target gene SAD by qRT-PCR.

Conclusions: In summary, this study evaluated and validated the suitable reference genes for qRT-PCR under different experiment treatment and tissues, and will be useful for further gene expression studies on the molecular mechanisms in Macadamia.

Background

Macadamia is a genus that belongs to the Proteaceae family, is a long-lived evergreen tree of subtropical and tropical origin, and native to Australia. The commercial cultivars is composed of Macadamia integrifolia Maiden and Bethche, Macadamia tetraphylla and their hybrids[1]. Macadamia produce edible dried fruit[2], and also are consumed as an ingredient in making pastry, chocolate, oil, cosmetic and pharmaceutical industry. Therefore macadamia is the most economically important crop. Macadamias have also been expanded in South Africa, Kenya, China, the United State, Guatemala, Malawi, Brazil, Vietnam, Mexico and Thailand currently. In China the planting area reach to 295856 hm² until 2018, and account for more than 65% of the total area of the world grow, becoming the world’s largest country with planting area.

The shells of Macadamia nut is very hard to crack as the high strength and toughness, it is considered as an interesting research material to study for developing of impact- and puncture-resistant engineering material[3]. Moreover, tocotrienols[4] and squalene[5] these bioactive compounds were detected in macadamia oil, and these compounds with high antioxidant ability can be beneficial to health maintenance and disease prevention[6]. Diets containing macadamia nuts reduce plasma LDL cholesterol levels and decrease the risk of CVD[7]. The kernel of macadamia nut is rich of lipids. Macadamia oil contains higher monounsaturated fatty acid than other vegetable oil[5, 8], and hence oil content as a major indicator of macadamia kernel quality, but the mechanisms accumulating the unsaturated fatty acid are still unknown. Gene expression analysis of fatty acid synthesis is beneficial to understand the complex regulatory networks of macadamia oil accumulation. In macadamia, soluble acy-ACP desaturases[9–10] and thiesterases[11]which were responsible for accumulation of unsaturated fatty acids were cloned and characterized, but the transcriptional level regulation need to be further elucidate. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) is one of the most frequently-used technologies to study the expression pattern of genes[12–13]. And it was widely used in the analysis of gene expression levels with high sensitivity, specificity, reproducibility and accuracy[14]. The accuracy of qRT-PCR results depends largely on the selection of internal reference genes. Selecting suitable reference genes as internal control under different experimental conditions is very necessary because the starting amount of material, RNA integrity, reverse transcription efficiency, cDNA quality, and different tissue and cell type can significantly influence the accuracy of gene expression [15–16]. Studies have shown that the expression of housekeeping genes is not constant in different cell types and physiological states. The expression of housekeeping genes will change with the change of experimental conditions and organ specificity. If the change of internal reference gene expression is not considered, the accuracy of qRT-PCR results will be reduced or even the opposite or wrong conclusion will be reached. Therefore, in order to decipher the molecular mechanisms of UFAs biosynthesis (accumulation) promoted by SAD and FAT genes in Macadamia integrifolia kernel, the need for selecting suitable reference gene for qRT-PCR in macadamia is urgent .

In this study, we selected 11 traditional candidate reference genes, α-tubulin (TUBa), β–tubulin (TUBb), malate dehydrogenase (MDH), 18S ribosome RNA (18S), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), α-elongation factor 1 (EF1a), β- elongation factor 1 (EF1b), ubiquitin (UBQ), ubiquitin-conjugating enzyme (UBC), cyclophilin (CYP) and actin (ACT), tested its expression pattern under various treatments by qRT-PCR. The expression stability of these candidate reference genes were further validated and evaluated using statistical algorithms includingΔCt[17], geNrom[18], Normfinder[19], and BestKeeper[20]. The comprehensive stability ranking of these reference genes under each specific experimental condition were also performed. Additionally, one target gene, Δ9-Stearoyl-ACP desaturases (SAD) which first desaturation introduction in the synthesis of aliphatic chains in plant was used to validate the effectiveness of the selected reference genes. Finally, this work provides the basis for further research in exploring of gene expression profiling in macadamia oil accumulation.

Results

In order to select the best reference genes, RefFinder online tool was used to rank the expression stability values. This method include comparative ΔCt, geNorm, NormFinder and BestKeeper, four popular algorithms, determinates a comprehensive expression stability ranking combined with above four data generated from different algorithms(Fig. 2; Table 2).

ΔCt algorithm analysis

The comparative ΔCt algorithm ranks the potential candidate reference genes expression stabilities based on the standard deviation values(Supplemental Table 3) [17]. As determined by ΔCt, ACT was consistently identified as the most stable gene through the total samples(Supplemental Fig. 3). CYP was the most stable reference gene in macadamia root and leaf, while ACT was the most stable reference gene in the stem in the different organs dataset from macadamia. In the subset of cold stress, NaCl stress, PEG stress, ABA treatment and MeJA treatment, ACT was still considered as the best stable gene in the 11 candidate reference genes, and TUBa was the least stable gene under the cold stress, NaCl stress, PEG stress, ABA treatment and MeJA treatment. Under the heat stress EF1b ranked top than ACT gene, and was determined the best stable reference gene. UBQ was ranked as the stable reference gene, while 18S was the least stable reference gene under the GA treatment. In the ethylene treatment, UBC was the most highly ranked, and EF1a performanced the least stable gene(Fig. 2; Table 2, Supplemental Table 3).

geNorm Analysis

The geNorm program used the M-value to assess stability of the reference gene expression. The M-value was calculated at each step during stepwise exclusion of the least stable reference gene until two best gens were obtained[18]. The top two stable candidate reference gene were ACT / MDH for all samples(Supplemental Fig. 4), leaf and MeJA treatment, ACT and EF1b under cold stress, heat stress, NaCl stress and PEG stress, MDH / UBC in the GA treatment and ABA treatment, ACT and UBC for stem and ethylene treatment show the most stability. Additionally, as show in Fig. 2, Table 2 and Supplemental Table 4, the least stably expressed was TUBa under cold stress, NaCl stress, ABA treatment, MeJA treatment, stem and all samples, UBC under heat stress and leaf, TUBb under PEG stress, 18S under the GA treatment and root, EF1a under the ethylene treatment.

NormFinder Analysis 

The NormFinder algorithm ranking candidate reference gene stability based on a variance estimation approach[19]. The NormFinder ranking orders are similar to the geNorm analysis. According this analysis(Figure S4; Table 2; Supplemental Table 5), the most stable gene were ACT for cold stress, ABA treatment, MeJA treatment, stem and all samples, EF1b for heat stress, CYP for leaf, ethylene treatment and NaCl stress, UBC for PEG stress, UBQ for GA treatment, EF1a for root. (Supplemental Fig. 5)

BestKeeper Analysis

The BestKeeper algorithm determined the stabilities of candidate reference genes based on the CV ± SD values[20]. The stability rankings by BestKeeper was quite different from the ranking order determined by ΔCt, geNorm and NormFinder algorithms(Supplemental Fig. 6). ACT was considered as the best stable gene under all samples, PEG stress, and leaf. TUBb in heat stress, UBC for cold stress, MeJA treatment and stem subset, 18S for NaCl stress and root, MDH for ABA treatment and ethylene treatment, EF1b for GA treatment were all ranked in top stable genes in different situation. Additionally, 18S and TUBa were determined as the least stable reference gene in the most samples by BestKeeper algorithm(Fig. 2; Table 2; Supplemental Table 6).

Comprehensive Ranking

In order to obtain a consensus result of the best stable reference genes as recommended by the four above approached according to the RefFinder[21], the geometric mean of four algorithms corresponding rankings for each candidate gene were calculated(Table 2; Fig. 2; Supplemental Tables 7 and 8). ACT was ranked the first stable gene in all samples, cold stress, NaCl stress, PEG stress, ABA treatment, MeJA treatment, stem tissue and leaf, EF1b for heat stress and GA treatment, CYP for root, UBC for ethylene treatment (Fig. 2; Table 2).

4. Validation of the reference genes with the SAD gene

To further validate the selected reference genes, the relative expression level of a target gene SAD in macadamia under different experimental conditions were evaluated using qRT-PCR. It was normalized using the most stable reference genes and two least stable reference genes as the internal control both singly and combination under different treatment subsets (Fig. 3). The relative expression abundance of SAD showed substantial divergence when normalized to different kinds of reference genes. In this study, the water treatment was designed as a negative control, when the most stable genes ACT and CYP were employed respectively and together to normalized expression levels of SAD gene, the expression abundant was not significantly effect; when the least stable genes TUBa were used to normalized expression level of target gene, the expression abundant was up-regulation dramatically. When the EF1a was used as the reference, the expression pattern of SAD gene was similar with the most stable genes normalization results. The expression level of SAD gene gradually declined under the cold stress using ACT and EF1b as the internal control, while the expression level was increased dramatically at 6h and 24h when the least stable reference genes TUBa was used, and increased sharply at 24h for UBQ normalized as internal reference. Under heat stress, the expression pattern of SAD gene were similar when the most stable reference genes (ACT and EF1b) and the least stable reference genes( UBC and UBQ) were used as controls, but the highest expression level at 12h was 2.5-fold higher by using UBQ normalizations than using the most stable reference gene ACT normalization. In PEG stress, the expression level of SAD gene was up-regulated at 24h and then down-regulated at 48h when using ACT and UBC genes normalization, while the expression appeared overestimated when normalized using GAPDH and underestimated for TUBb. Similarly, in NaCl stress, the expression profiles of the SAD showed same trends when stable reference genes ACT and CYP were used as internal control, when the least stable reference gene TUBb and TUBa were used as reference control, the expression level of target gene was significantly overestimated and underestimated respectively. Under ET and ABA treatment, the relative expression trends of SAD gene showed a similarly trend, appeared slightly descent at 6h and then increased at 24h when the most stable genes were used in combination and independent for normalization, while appeared down-regulated at 12h and 24h when normalized by using the least stable reference genes. In MeJA stress, a similar expression pattern was observed when either ACT or MDH was used alone and combination for normalization. While the expression levels of SAD gene were overestimated when the least stable genes TUBa and EF1a were used for normalization. In GA stress, the relative expression of SAD gene increased at 6h when using EF1b and 18S as internal control, while declined slightly at 6h when normalized by using TUBa and UBQ as reference genes. Above all, the result showed that ACT and EF1b were more suitable for different treatment subset and TUBa failed standardized the expression data effectively.

In addition, to validate the selected reference genes, 4 genes involved fatty acid accumulation including FATA, FAD, Oleosin and SAD in qRT-PCR were compared with the expression pattern in RNA-seq results. In the qRT-PCR results, Actin was used as the normalizer to quantify the expression. As show in the Fig. 4, all the selected genes revealed similar expression file in the both methods. Therefore, our result provide the reliable reference gene in Macadamia.

Discussion

Gene expression quantification which as an important way for charactering gene function, has been accepted and applied widely in the field of genetic research[22]. qRT-PCR is commonly regarded as an practical method for accurate analysis of gene expression profiling with high sensitivity, specificity, reproducibility and less samples as MIQE(minimum information for publication of quantitative real-time PCR experiments, MIQE) a precondition[23]. There is no universally suitable reference genes under each experimental condition, and the reference genes expression is tissue-specific and various based on the physiological status of the organs or experimental condition[24–25]. Numerous studies have reported that the expression of housekeeping gene may be constant in some experiment condition, but varied considerably in other cases[26–28]. The nature protocols suggest that using a specifically and stably reference for normalization could ensure accurate results[29]. Therefore suitable reference genes selection and evaluation is critically mandatory for data normalization in gene expression analysis[30]. Suitable reference genes selection has been studied in the plant field recent years, including cucumber[31], Setaria viridis[32], sugarcane[33], wheat[34], soybean[35], Brassica napus[36], Taxus spp.[37], Gentiana macrophylla[38] Botrytis cinerea[39] and other various plants. To date, few studies have been carried out in macadamia involved in reference selecting. This has hindered the molecular functional studies of important genes involved in stress conditions and specific organisms. In present study, we systematically validated the stability of 11 traditional reference gene under the specific environmental condition or treatment as well as the in different tissues of Macadamia integrifolia seedlings. Traditional housekeeping genes involved in cytoskeleton structure(ACT, TUBa, TUBb), protein synthesis (EF1a, EF1b, 18S), biological metabolic process(GAPDH and UBQ), and multifunctional proteins (CYP and MDH), are selected and the open reading frame sequences of these genes were cloned from Macadamia integrifolia.

It is recommended that a comprehensive RefFinder algorithms could provide a reliable assessment of the expression stability of candidate reference gene under various experiment condition for real-time qRT-PCR analysis[40]. This statistical algorithms four methods, includingΔCt, geNorm, NormFinder and BestKeeper, that have been developed to assess the expression stability of candidate reference genes for accurate normalization in gene expression studies were used to rank the reference gene. The ranking results by different algorithms are relatively coincident mostly. In the heat treatment subset, EF1b and ACT are top two genes inΔCt, NormFinder and geNorm analysis, and ranked third and fourth in bestkeeper analysis. In ABA treatment, NaCl, stem, and all samples sets, TUBa was always ranked the least stable gene by four different algorithms. But sometimes the stability ranking showed significant differences due to the complementary different statistical program (Table 2). For example, under NaCl stress, ACT was ranked first by ΔCt and geNorm, while it was ranked sixth by BestKeeper. CYP was ranked first by NormFinder, and followed by ACT. The ΔCt method by comparing relative expression of pairs of genes within each sample to confidently identify stability of reference genes[17]. geNorm is one of the most commonly used method for systematic validation procedures of reference gene stability[18]. geNorm identifies two reference genes according to the similarity of expression profile in each sample and the lowest intra-group variation [41–42], and not suitable for distinguishing the gene of similar expression pattern. Andersen developed the NormFinder approach[19], which takes both the iner- and intra-group variations into account, and combines them into a stability value, and finally ranks the top genes with minimal inter- and intra-group variation. Most of time the results of NormFinder is similar to the geNorm results. For instance in the heat stress samples, EF1b and ACT were ranked top two stable reference genes, and TUBa and UBQ were the least stable reference genes, other genes ranked order have slightly change by geNorm and NormFinder. The BestKeeper software analysis is not restricted to the stability of reference genes but also target gene[20], it determines the stability ranking of the reference genes according to their SD and CV values[43–44]. Compare with geNorm and NormFinder, the BestKeeper results have significant difference owing to the different calculation strategy. For example, in the subset of all samples, ACT was ranked first genes by all the programs. Under cold stress, ACT was determined the most stable genes byΔCt, geNorm, NormFinder and comprehensive ranking, excepted BestKeeper. As these programs ranked different most stable reference genes [45–47], we selected the most reliable reference gene by using a comprehensive tool RefFinder, which that rank the reference genes based on the geometric mean of the weights of every gene calculated by each program[21]. The comprehensive results obtained on the basis of statistical analysis byΔCt, geNorm, NormFinder and BestKeeper showed mostly consistency in this study.

ACT is commonly known to be a reliable reference gene[48]. Although there are studies have provided evidence that ACT vary considerably in Zea mays [49], our finding support that ACT is the most stable reference gene for the relative quantification in Macadamia among all samples, cold stress, salt stress, PEG stress, ABA treatment and MeJA treatment samples. Elongation factors involved in translocation of aminoacyl-transfer RNA to the ribosome during protein synthesis. EF1a and EF1b are highly conserved in sequence and expression in eukaryotes, usually used as internal genes. EF1b was established as the most appropriate reference gene under GA treatment and heat treatment in this study, while EF1a is considered the most stable gene under complex developmental process in Populus [50]. UBC ranked first than ACT gene in ethylene treatment which consistent with the result in platycladus orientali [43]. CYP(Cyclophilin), a specific cytosolic binding protein, which is ranked the top stable gene in root in macadamia. TUBa mostly was not satisfactory for qRT-PCR analysis in present study, but it was sufficiently stable in Hererosigma akashiwo (Raphidophyceae)[51], which has been selected as reference gene in Ulva linza [52] and the diatom Pseudo-nitzschia multistriata[53].

Our previous studies found thatΔ⁹-Stearoyl-ACP desaturases(SAD) is the key of unsaturated fatty acid synthesis and accumulation in macadamia oil (no publish data). The expression profile of SAD and the fatty acid level were affected in Arabidopsis Crown Galls under drought and hypoxia stress conditions[54]. In this study, we quantified the SAD gene expression using two most stable reference genes both singly and combined as well as two least stable genes as an internal control in qRT-PCR analysis to validate the selected reference genes in macadamia. When the least stable genes TUBb was validated as a reference gene for normalization, the expression pattern of SAD gene was obviously overestimated, and underestimated using TUBa to reference under the NaCl treatment (Fig. 3). The SAD gene expression was induced using the stable reference gene, while down-regulation when normalization with the two least reference gene TUBa and EF1a after the 6h ET treatment. So these results demonstrated that the suitable selection of internal control is critically important for the normalization, and used reference without validation would be leading misinterpretation.

RNA-seq is an approach to analysis the transcriptome profiling of various species, and recently was used to search for candidate reference genes. In this study, the target gene expression profile consistent both in qRT-PCR results and RNA-seq results, and they supported each other. So, our finding indicated that the result of this experiment are credible. The reference genes evaluated in this study would be very useful for further gene expression analysis in molecular mechanism of fatty acid accumulation under different treatment in Macadamia.

Conclusions

The reference genes selected in current study will be helpful for accurate normalization of qRT-PCR data in M. integrifolia. In this study, we selected the most stable reference genes in different tissue of macadamia seedlings and under different treatment. We found that ACT shows a good stability in all samples. UBC could be used to normalization when treated by ethylene. In the GA treatment samples, EF1b is the best reference genes. The identification of the suitable reference genes in this study will facilitate the future work on gene expression studies in Macadamia to improve our standing on the molecular mechanisms of fatty acid synthesis and accumulation under different experiment condition.

Methods

Plant materials and stress treatments.

NanYa No1, a macadamia cultivar developed by the South Subtropical Crops Research Institute of China Academy of Tropical Agricultural Sciences(CATAS). The Seeds of cultivars NanYa No1 were provided by National Field Genebank for Tropical Fruits of the South Subtropical Crops Research Institute，Guangdong Province and Chinese Academy of Tropical Agricultural Sciences, Zhanjiang, China(110º27′E; 21º17′N). Macadamia seeds were sowed and maintained in a glass bottle with 1/2 SD medium(without hormone) kept at 20℃±2 . When the seedlings stem of macadamia sizes were 2.1-3.3 mm in diameter, different treatment could be conducted in this study.

For experiments, each treatment group was set up in triplicate. Control group was treated with water. The organ-specific series of samples (root, stem, leaf) were collected from the seedlings. For drought treatment, 20% PEG-6000 solution(w/v, polythyleneglycol) was applied to incubate the plants for 0, 12, 24, 48 and 72h. For cold and heat stress, plants in the pots were placed at chamber with the temperature of 4℃ and 42℃ respectively for 0, 2, 6, 12 and 24h. For salinity treatment, seedling were transferred to 1/2 SD medium with 150 mmol NaCl for 0, 12, 24, 48, 72h. For hormone treatment, plants were imposed in 100 μmol/L methyl jasmonate(MeJA), abscisic acid(ABA), gibberrellins(GA) or ethylene (ETH) for 0, 2, 6, 12 and 24h. After that, the roots, leaves and stems were sampled separately at different periods for each treatment used for expression analysis. The detailed information of samples collected from various tissues/experimental conditions were also listed in Table S1. All samples were frozen in liquid nitrogen and then stored at -80℃ prior to RNA isolation. 

Total RNA Isolation and cDNA Synthesis

The frozen samples were ground to a fine powder in liquid nitrogen using a pestle and mortar. The total RNA was extracted from the samples using Plant Total RNA Isolation Kit Plus (Cat. No.RE-05021/05022/05024, FOREGNE, China) following the manufacturer’s recommendation. RNA purity and concentration were measured by NanoDrop spectrophotometer (Thermo, Germany), and RNA integrity was verified by agrose gel electrophoresis with ethidium bromide staining. First-stand cDNAs were synthesized from 1μg total RNA with a final volume of 20 μL using the RTIII all-in-one Mix (Monad, China). The cDNA synthesis mixture was incubated at 37℃ for 2 minutes and then 55℃ for 15 minutes. This reaction was terminated by incubating it at 85℃ for 10 seconds. The cDNAs were diluted 1:50 with nuclease-free water before use in subsequent experiments. 

Selection of Candidate Reference Genes and Primer Design

Eleven classical housekeeping genes commonly used as internal control for expression studies. The open reading frame sequence of these genes were cloned from macadamia（GenBank accession number MN650750, MN627205, MN627206, MN627207, MN627208, MN627209, MN627210, MN627211, MN627212, MN627213, MN627214）. The candidate reference genes were listed in table 1 with their respective reference where they were first described. The amplification primers for real-time PCR were designed using Primer-BLAST in NCBI. The criteria for primer design were set as follows: primer lengths of 20-24 bp, GC contents of 45-55%, melting temperature(Tm) in a range of 55-60℃ and amplification lengths of 100-250 bp.  

Real-Time Quantitative Polymerase Chain Reaction

qRT-PCR was conducted in 96-well plates and performed on the LightCycler 480 (Roche Molecular Biochemicals, Mannheim, Germany). The reaction mix contained 1μL diluted cDNA，5μL Tap SYBR Green qPCR Premix (Nova, Jiangsu, China), 0.5ΜlμL the specific primer pairs, and 3.5μL ddH2O in a final volume of 10μL. Three biological replicates for all of the samples and three technical replicates of each biological replicate were carried out. The qRT-PCR protocol was follows: 95℃ for 3 min, 45 cycles of 95℃ 10 s, 60℃ for 30 s, 72℃ for 10s. To verified the specificity of each primer, a melting-curve analysis was include. The mean amplification efficiency of each primer pair was checked by the LinRegPCR program[55]. 

Analysis of real-time PCR data

The comparative ΔCt algorithm assesses the expression stability of reference genes by the mean standard deviation values from the test sample[17]. The ΔCt value indicated the variably transcription of candidates.

The geNorm VBA applet determines the best reference based on the geometric mean and pair wise variation of each gene from all the candidate reference genes in the total sample. Two parameters were used to evaluate the expression stability of the reference genes; the average expression stability value(M value), based on the pairwise variation between a particular gene compared to all others, and the pairwise variation(Vn/n+1), which determines the required number of genes to result in a more accurate normalization[18]. A cut-off value of Vn/n+1 <0.15 indicates that an additional reference gene makes no significant contribution to the normalization. 

While geNorm performs a stepwise exclusion of the least stably expressed genes, NormFinder uses a model-based approach, which calculates both inter- and intra-group variability to estimate the stability of gene expression. NormFinder identified the optimal normalization genes among a panel of candidates according to their expression stability in a given sample set or given experimental designs. This algorithm evaluates not only the overall expression variation of the candidate reference gens, but also the variation between subgroup of samples. In addition to data from real-time PCR , NormFinder can also analyze expression data obtained through other quantitative gene expression methods, such as microarray. Ranking the stability based on the stability value (SV) where lower stability value represents higher gene expression stability and vice versa[19]. 

Bestkeeper is used to determinate the stability of reference genes based on the coefficient of variance (CV) and the standard deviation(SD) of the average Ct values. The lower of CV±SD values, the more stable the candidate reference genes, and gene with SD >1 would be considered unacceptable and should be excluded[20].

RefFinder is a user-friendly web-based comprehensive tool developed for evaluating and screening reference genes from extensive experimental datasets. It integrates the currently available major computational programs (geNorm, Normfinder, BestKeeper, and the comparative ΔCt method) to compare and rank the tested candidate reference genes according to the geometric mean of four algorithms corresponding rankings. [21]. We analysis the qRT-PCR data on the RefFinder mirror site(http://blooge.cn/RefFinder/) . 

Evaluating Reference Genes Expression

The expression level of the target gene SAD was analyzed using the most stable and least stable reference genes after normalization across all the experimental sets. The amplification efficiencies of the target genes were also estimated by the LinRegPCR program. The average Ct value was calculated from three biological and technical replicates and used for relative expression analyses. The relative expression data were calculated according to the qRT-PCR 2^-ΔCt method and presented as fold change[56].

Abbreviations

SAD：Δ9-Stearoyl-ACP desaturases 

18S：18S ribosome RNA 

ACT：Actin 

CYP：Cyclophilin 

EF1a：Elongation factor 1-α 

EF1b：Elongation factor 1-β

GAPDH：Glyceraldehyde-3-phosphate dehydrogenase 

MDH：Malate dehydrogenease

UBC：Ubiquitin-conjugating enzyme 

TUBa：Tubulin-α 

TUBb：Tubulin -β

UBQ：Ubiquitin-A

FATA：fatty acid thioesterase

FAD：fatty acid desaturation

MeJA：jasmonate

ABA：abscisic acid 

GA ：gibberrellins 

ETH：ethylene

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

This manuscript has been reviewed and approved by all authors.

Availability of data and materials

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Competing interests

The authors declare that they have no competing interests.

Funding

This research was supported by Guangdong Basic and Applied Basic Research Foundation (No. 2021A1515012421), Central Public-interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences (No. 1630062021015, No.1630062020001), Hainan Provincial Natural Science Foundation of China (No. xxxxxxxx). The funding bodies played no role in study design, the collection and analysis of the data, data interpretation, and in writing the manuscript.

Authors’ contributions

Q.Y and Z.P.Y contributed equally to this work. Q.Y performed the experiments, prepared the figures and tables, and drafted the manuscript. Z.P.Y participated in design of the experiments and performed the statistical analysis. Y.L.Z prepared the plants and performed part experiments. H.Z and M.H.Z provided the Macadamia tree of NanYa No1 cultivars to collect Macadamia seeds. H.L revised the manuscript. All authors have read and approved the manuscript.

Acknowledgements

We thank National Tropical Plants Germplasm Resource Center and National Field Genebank for Tropical Fruits for providing macadamia seeds.

References

1. Hardner CM, Peace C, Lowe AJ, Neal J, Pisanu P, Powell M, Schmidt A, Spain C, Williams K. Genetic resources and domestication of macadamia. Horticultural Reviews. 2009; 35:1–125.
2. Cavaletto CG. Macadamia nuts. Handbook of tropical foods (Ed. HT Chan), New York: Marcel Dekker, Inc.1983; 361–397.
3. Schüler P, Speck T, Bührig-Polaczek A, Fleck C. Structure-function relationships in Macadamia integrifolia seed coats-fundamentals of the hierarchical microstructure. PLoS One. 2014; 9(8): e102913.
4. Maguire LS, O'Sullivan SM, Galvin K, O'Connor TP, O'Brien NM. Fatty acid profile, tocopherol, squalene and phytosterol content of walnuts, almonds, peanuts, hazelnuts and the macadamia nut. International Journal of Food Sciences and Nutrition. 2004; 55(3): 171–178.
5. Rengel A, Pérez E, Piombo G, Ricci J, Servent A, Tapia MS, Gibert O, Montet D. Lipid profile and antioxidant activity of macadamia nuts (Macadamia integrifolia) cultivated in Venezuela. Nature Science. 2015; 7: 535–547.
6. Vadivel V, Kunyanga CN, Biesalski HK. Health benefits of nut consumption with special reference to body weight control. Nutrition. 2012; 28:1089–1097.
7. Griel AE, Cao Y, Bagshaw DD, Cifelli AM, Holub B, Kris-Etherton PM. A macadamia nut-rich diet reduces total and LDL-cholesterol in mildly hypercholesterolemic men and women. The Journal of Nutrition. 2008; 138(4): 761e767.
8. Wang XF, Yang RN, Xue L, Zhang LX, Wang XP, Zhang Q, Li PW. Determination of fatty acid composition of 28 kinds of functional vegetable oil, Journal of Food Safety and Quality. 2017; 8(11):4336–4343.
9. Gummeson PO, Lenman M, Lee M, Singh S, Stymne S. Characterisation of acyl-ACP desaturases from Macadamia integrifolia Maiden & Betche and Nerium oleander L. Plant Science. 2000;154(1): 53–60.
10. Rodríguez MFR, Sánchez-García A, Salas JJ, Garcés R, Martínez-Forc E. Characterization of soluble acyl-ACP desaturases from Camelina sativa, Macadamia tetraphylla and Dolichandra unguis-cati. Journal of Plant Physiology. 2015; 178: 35–42.
11. Moreno-Pérez AJ, Sánchez-García A, Salas JJ, Garcés R, Martínez-Force E. Acyl-ACP thioesterases from macadamia (Macadamia tetraphylla) nuts:Cloning, characterization and their impact on oil composition. Plant Physiology and Biochemistry. 2011; 49: 82e87
12. Wong ML, Medrano JF. Real-time PCR for mRNA quantitation. BioTechniques. 2005; 39(1): 75–85.
13. Bustin SA, Benes V, Nolan T, Pfaffl MW. Quantitative real-time RT-PCR-a perspective. Journal of Molecular Endocrinology. 2005; 34: 597–601.
14. VanGuilder HD, Vrana KE, Freeman WM. Twenty-five years of quantitative PCR for gene expression analysis. BioTechniques. 2008; 44(5):619–626.
15. Bustin SA. Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): Trends and problems. Journal of Molecular Endocrinology. 2002; 29(1): 23–39.
16. Huggett J, Dheda K, Bustin S, Zumla A. Real-time RT-PCR normalization; strategies and considerations. Genes & Immunnity. 2005; 6(4):279–284.
17. Silver N, Best S, Jiang J, Thein S. Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Molecular Biology. 2006; 7(1):1–9.
18. Vandesompele J, Preter KD, Pattyn F, Poppe B, Roy NV, Paepe AD, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology. 2002; 3(7): 00341–003411.
19. Andersen CL, Jensen JL, Ørntoft TF. 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. 2004; 64(15):5245–5250.
20. Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnology letters. 2004; 26(6):509–515.
21. Xie F, Xiao P, Chen D, Xu L, Zhang B. miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs. Plant Molecular Biology. 2012; 80(1): 75–84.
22. Bustin SA. Real-time, fluorescence-based quantitative PCR: a snapshot of current procedures and preferences. Expert Review of Molecular Diagnostics. 2005; 5 (4): 493–498.
23. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clinical Chemistry. 2009; 55(4): 611–622.
24. Huang Y, Tan H, Yu J, Chen Y, Guo Z, Wang G, Zhang Q, Chen J, Zhang L, Diao Y. Stable internal reference genes for normalizing real-time quantitative PCR in Baphicacanthus cusia under Hormonal Stimuli and UV Irradiation, and in different plant organs. Frontiers in plant science. 2017; 8:668.
25. Ma R, Xu S, Zhao Y, Xia B, Wang R. Selection and validation of appropriate reference genes for quantitative real-Time PCR analysis of gene expression in Lycoris aurea. Frontier in Plant Science. 2016; 7:536.
26. Roman A. Volkov, Irina I. Panchuk, Fritz Schöffl. Heat-stress-dependency and developmental modulation of gene expression: The potential of house-keeping genes as internal standards in mRNA expression profiling using real-time RT-PCR. Journal of Experimental Botany. 2003; 54(391): 2343–2349.
27. Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiology. 2005; 139(1): 5–17.
28. Remans T, Smeets K, Opdenakker K, Mathijsen D, Vangronsveld J, Cuypers A. Normalisation of real-time RT-PCR gene expression measurements in Arabidopsis thaliana exposed to increased metal concentrations. Planta. 2008; 227(6): 1343–1349.
29. Nolan T, Hands RE, Bustin SA. Quantification of mRNA using real-time PCR. Nature Protocls. 2006;1: 1559–1582.
30. Brunner AM, Yakovlev IA, Strauss SH. Validating internal controls for quantitative plant gene expression studies. BMC Plant Biology. 2004; 4,14.
31. Migocka M, Papierniak A. Identification of suitable reference genes for studying gene expression in cucumber plants subjected to abiotic stress and growth regulators. Molecular Breeding. 2011; 28(3): 343–357.
32. Martins PK, Mafra V, Souza WRD, Ribeiro AP, Vinecky F, Basso MF, da Cunha BA, Kobayashi AK, Molinari HB. Selection of reliable reference genes for RT-qPCR analysis during developmental stages and abiotic stress in Setaria viridis. Scientific Reports. 2016, 6:28348.
33. Guo J, Ling H, Wu Q, Xu L, Que Y. The choice of reference genes for assessing gene expression in sugarcane under salinity and drought stresses. Scientific Reports. 2014; 4, 7042.
34. Long XY, Wang JR, Ouellet T, Rocheleau H, Wei YM, Pu ZE, Jiang QT, Lan XJ, Zheng YL. Genome-wide identification and evaluation of novel internal control genes for Q-PCR based transcript normalization in wheat. Plant Molecular Biology. 2010; 74(3):307–311.
35. Ma S, Niu H, Liu C, Zhang J, Hou C, Wang D. Expression stabilities of candidate reference genes for RT-qPCR under different stress conditions in soybean. PLoS One. 2013; 8(10): e75271.
36. Wang Z, Chen Y, Fang H, Shi H, Chen K, Zhang Z, Tan X. Selection of reference genes for quantitative reverse-transcription polymerase chain reaction normalization in Brassica napus under various stress conditions. Molecular Genetics and Genomics. 2014; 289(5): 1023–1035.
37. Zhang K, Fan W, Chen D, Jiang L, Li Y, Yao Z, Yang Y, Qiu D. Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp.. Scientific Reports. 2020; 10, 22205.
38. He Y, Yan H, Hua W, Huang Y, Wang Z. Selection and validation of reference genes for quantitative real-time PCR in Gentiana macrophylla. Frontier in Plant Science. 2016; 7:945.
39. Ren H, Wu X, Lyu Y, Zhou H, Xie X, Zhang X, Yang H. Selection of reliable reference genes for gene expression studies in Botrytis cinerea. Journal of Microbiological Methods. 2017; 142: 71–75.
40. Chen Y, Hu B, Tan Z, Liu J, Yang Z, Li Z, Huang B. Selection of reference genes for quantitative real-time PCR normalization in creeping bentgrass involved in four abiotic stresses. Plant Cell Reports. 2015; 34(10):1825–1834.
41. Jian B, Liu B, Bi Y, Hou W, Wu C, Han T. Validation of internal control for gene expression study in soybean by quantitative real-time PCR. BMC Molecular Biololgy. 2008; 9(59):1–14.
42. Cruz F, Kalaoun S, Nobile P, Colombo C, Almeida J, Barros LMG, Romano E, Grossi-de-Sá MF, Vaslin M, Alves-Ferreira M. Evaluation of coffee reference genes for relative expression studies by quantitative real-time RT-PCR. Molecular Breeding. 2009; 23(4): 607–616.
43. Chang E, Shi S, Liu J, Cheng T, Xue L, Yang X, Yang W, Lan Q, Jiang Z. Selection of reference genes for quantitative gene expression studies in Platycladus orientalis (Cupressaceae) using real-time PCR. PLoS ONE. 2012; 7(3):e33278.
44. Xiao X, Ma J, Wang J, Wu X, Li P, Yao Y. Validation of suitable reference genes for gene expression analysis in the halophyte Salicornia europaea by real-time quantitative PCR. Frontier in Plant Science. 2015; 5:788.
45. Jacob F, Guertler R, Naim S, Nixdorf S, Fedier A, Hacker NF, Heinzelmann-Schwarz V. Careful selection of reference genes is required for reliable performance of RT-qPCR in human normal and cancer cell lines. PLoS ONE. 2013;8(3): e59180.
46. Petriccione M, Mastrobuoni F, Zampella L, Scortichini M. Reference gene selection for normalization of RT-qPCR gene expression data from Actinidia deliciosa leaves infected with Pseudomonas syringae pv. actinidiae. Scienctific Reports. 2015; 5, 16961.
47. Rivera L, Lopez-Patino M, Milton D, Nieto T, Farto R. Effective qPCR methodology to quantify the expression of virulence genes in Aeromonas salmonicida subsp. salmonicida. Journal of Applied Microbiology. 2015; 118(4):792–802.
48. Pu Q, Li Z, Nie G, Zhou J, Liu L, Peng Y. Selection and validation of reference genes for quantitative real-time PCR in white clover (Trifolium repens L.) involved in five abiotic stresses. Plants. 2020; 9(8):996.
49. Díaz-Camino C, Conde R, Ovsenek N, Villanueva MA. Actin expression is induced and three isoforms are differentially expressed during germination in Zea mays. Journal of Experimental Botany. 2005; 56(412):557–565.
50. Xu M, Zhang B, Su X, Zhang S, Huang M. Reference gene selection for quantitative real-time polymerase chain reaction in Populus. Analytical biochemistry. 2011; 408(2): 337–339.
51. Ji N, Li L, Lin L, Lin S. Screening for suitable reference genes for quantitative real-Time PCR in Heterosigma akashiwo (Raphidophyceae). PLoS ONE. 2015; 10(7): e0132183.
52. Dong M, Zhang X, Chi X, Mou S, Xu J, Xu D, Wang W, Ye N. The validity of a reference gene is highly dependent on the experimental conditions in green alga Ulva linza. Current Genetics. 2012; 58(1): 13–20.
53. Adelfi MG, Borra M, Sanges R, Montresor M, Fontana A, Ferrante MI. Selection and validation of reference genes for qPCR analysis in the pennate diatoms Pseudo- nitzschia multistriata and P. arenysensis. Journal of Experimental Marine Biology and Ecology. 2014; 451: 74–81.
54. Klinkenberg J, Faist H, Saupe S, Lambertz S, Krischke M, Stingl N, Fekete A, Mueller MJ, Feussner I, Hedrich R, Deeken R. Two fatty acid desaturases, STEAROYL-ACYL CARRIER PROTEIN ∆9-DESATURASE6 and FATTY ACID DESATURASE3, are involved in drought and hypoxia stress signaling in Arabidopsis crown galls. Plant Physiology. 2014; 164 (2): 570–583.
55. Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, Moorman AF. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Research. 2009; 37(6), e45.
56. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 – ∆∆CT method. Methods. 2001; 25 (4) :402–408.