Effect of Polyploidy Induction on Natural Metabolite Production in Medicinal Plants

Regulatory trends in engineering bioactive-phytocompounds

The secondary plant metabolites are of enormous importance because of their extensive medicinal, nutraceutical, and industrial applications. In plants, these secondary metabolites are often found in extremely small amounts, therefore, following the discovery of any prospective metabolite, the main constraining element is the ability to generate enough material for use in both industrial and therapeutic settings. In order to satisfy the rising demand for value-added metabolites, researchers prefer to use different molecular approaches for scalable and sustainable production of these phytocompounds. Here, we discuss the emerging regulatory trends in engineering these bioactive-phytocompounds and provide recommendation on successful employment of these state-of-the-art technologies for translation of these academic researches into novel process and products.

Promising approaches for simultaneous enhancement of medicinally significant benzylisoquinoline alkaloids in opium poppy

Benzylisoquinoline alkaloids (BIAs) produced in opium poppy have been evidenced to heal patients suffering from various diseases. They, therefore, hold an integral position in the herbal drug industry. Despite the adoption of several approaches for the large-scale production of BIAs, opium poppy remains the only platform in this purpose. The only disadvantage associated with producing BIAs in the plant is their small quantity. Thus, recruiting strategies that boost their levels is deemed necessary. All the methods which have been employed so far are just able to enhance a maximum of two BIAs. Thus, if these methods are utilized, a sizable amount of time and budget must be spent on the synthesis of all BIAs. Hence, the exploitation of strategies which increase the content of all BIAs at the same time is more commercially effective and time-saving, avoiding the laborious step of resolving the biosynthetic pathway of each compound. Exposure to biotic and abiotic elicitors, development of a synthetic auto-tetraploid, overexpression of a WRKY transcription factor, formation of an artificial metabolon, and suppression of a gene in the shikimate pathway and miRNA are strategies that turn opium poppy into a versatile bioreactor for the concurrent and massive production of BIAs. The last three strategies have never been applied for BIA biosynthetic pathways.

Biological Properties of Recently Described Wild Bramble Rubus oklejewiczii against the Species from Similar Niches

The aim of this study was to compare the biological properties, such as antiviral, antibacterial, and antioxidant, of recently described pentaploid species Rubus oklejewiczii with tetraploid taxa growing in similar habitats including R. plicatus, R. gracilis, and R. wimmerianus. The antiviral potential was analyzed against bacteriophages with different genetic material: phi6 (a surrogate for the SARS-CoV-2 virus), T7, phiX174, and MS2. Antibacterial properties of fruit and leaf extracts were determined against Staphylococcus aureus, Bacillus cereus, Escherichia coli, and Salmonella enterica. The total phenolic content, as well as anthocyanins, ascorbic acid, pH, and antioxidant properties (FRAP and DPPH) were determined. R. oklejewiczii leaf extract was characterized by the weakest antibacterial and antiviral properties, which was closely correlated with the lowest content of polyphenolic compounds and antioxidant properties. The strongest biological properties were observed for R. wimmerianus leaves. Fruit extracts were characterized by lower phenolic content and antioxidant activities than leaves, with the lowest values observed for R. oklejewiczii extract. The antibacterial properties of fruit extracts were strongest for R. gracilis. The strongest antiviral potential was observed for R. oklejewiczii and R. wimmerianus fruit extracts against the bacteriophage phi6, which correlated with the lowest pH and the highest ascorbic acid content. The positive effect of the higher ploidy of R. oklejewiczii for most of the analyzed biological properties was not observed except for the antiviral potential of fruit extract. Due to its large and tasty fruits, this species seems to be very promising for cultivation and attractive for consumers, even though most of its biological properties were not any better compared to other examined tetraploid species.

Comparative Analysis of Root Phenolic Profiles and Antioxidant Activity of Five Native and Invasive Solidago L. Species

The high environmental importance of invasive goldenrod has prompted research to find potential benefits that can be derived from these species. This study aimed to identify differences in root phenolic profiles among five Solidago species, thus providing valuable information on their potential applications and the botanical origin of the raw material. The roots of native S. virgaurea L., two alien species S. gigantea Aiton and S. canadensis L. and their hybrids S. ×niederederi Khek and S. ×snarskisii Gudž. & Žaln. were sampled from mixed-species stands in Lithuania. A complex of twelve phenolic acids and their derivatives was identified and quantified in methanol-water root extracts using the HPLC-PDA and LC/MS systems. The radical-scavenging capacities of the extracts were assessed by ABTS. The chemical content of the roots of S. virgaurea, S. gigantea and S. ×niederederi were statistically similar, while the roots of S. canadensis and S. ×snarskisii contained lower amounts of compounds than the other species. The PCA score-plot models of the phenolic profiles only partially confirmed the identification of S. ×niederederi and S. ×snarskisii as crosses between native and alien species. The findings from the phenolic profiles and the observed radical-scavenging activity of root extracts of Solidago species provide valuable insights into their potential applications in various fields.

Floral volatile composition of Jasminum sambac variants developed through colchicine

Jasmines are commercially grown for their fragrant flowers and essential oil. The present study investigates the composition of the volatile compounds from flowers of Jasminum sambac cv. Ramanthapuram Gundumalli and its variants that were evolved through colchicine. GC-MS analysis revealed that the flowers possessed major terpenes and sesquiterpenes such as Linalool, α-farnesene, germacrene-D, geranyl Linalool and D-Limonene as well as benzenoids (including benzyl acetate, benzyl alcohol and (Z)-Cinnamyl benzoate). The relative abundance of these volatile compounds in the variants have shown higher percentages than their wild-type (parent) which indicates that the variants possessed enhanced volatile composition. The new variations generated in floral volatile composition of J. sambac through polyploidisation are likely to have significant impact on the loose flower and perfume industries. Besides, the identified unique compounds can also be used as metabolic signatures to characterise the novel variants.

Transcriptomic and metabolomic analyses reveal the differential accumulation of phenylpropanoids and terpenoids in hemp autotetraploid and its diploid progenitor

BACKGROUND

Cannabis sativa, a dioecious plant that has been cultivated worldwide for thousands of years, is known for its secondary metabolites, especially cannabinoids, which possess several medicinal effects. In this study, we investigated the autopolyploidization effects on the biosynthesis and accumulation of these metabolites, transcriptomic and metabolomic analyses were performed to explore the gene expression and metabolic variations in industrial hemp autotetraploids and their diploid progenitors.

RESULTS

Through these analyses, we obtained 1,663 differentially expressed metabolites and 1,103 differentially expressed genes. Integrative analysis revealed that phenylpropanoid and terpenoid biosynthesis were regulated by polyploidization. No substantial differences were found in the cannabidiol or tetrahydrocannabinol content between tetraploids and diploids. Following polyploidization, some transcription factors, including nine bHLH and eight MYB transcription factors, affected the metabolic biosynthesis as regulators. Additionally, several pivotal catalytic genes, such as flavonol synthase/flavanone 3-hydroxylase, related to the phenylpropanoid metabolic pathway, were identified as being modulated by polyploidization.

CONCLUSIONS

This study enhances the overall understanding of the impact of autopolyploidization in C. sativa and the findings may encourage the application of polyploid breeding for increasing the content of important secondary metabolites in industrial hemp.

The effect of induced polyploidy on phytochemistry, cellular organelles and the expression of genes involved in thymol and carvacrol biosynthetic pathway in thyme (Thymus vulgaris)

Induced polyploidy usually results in larger vegetative and reproductive plant organs. In order to study the effect of chromosome doubling on Thymus vulgaris, three levels of colchicine concentration including 0.1, 0.3 and 0.5% (w/v) were applied for 6, 12 and 24 hours on apical meristem of 2- and 4-leaf seedlings. Ploidy level was evaluated by flow cytometry and microscopic chromosome counting. Chemical composition of essential oils extracted by hydro-distillation was analyzed by gas Chromatography/mass spectrometry (GC/MS) and gas Chromatography (GC). The application of 0.3% colchicine at 4-leaf seedling for 6 hours resulted in the highest survival rate and the highest number of tetraploid plants. Cytogenetic and flow cytometry analyses confirmed the increase of chromosome number from 2n=2x=30 in diploids to 2n=4x=60 in induced tetraploids. Tetraploid plants had larger leaves, taller and thicker stems, dense branching, longer trichome, larger stomata, larger guard cells, and decreased number of stomata. The number of chloroplasts and mitochondria increased significantly in tetraploid plants by 1.66 and 1.63 times, respectively. The expression of CYP71D178, CYP71D180 and CYP71D181 increased in tetraploids by 3.27, 7.39 and 2.15 times, respectively, probably resulting in higher essential oil compounds, as tetraploids outyielded the diploid plants by 64.7% in essential oil, 40.9% in thymol and 18.6% in carvacrol content.

Polyploidy as a strategy to increase taxane production in yew cell cultures: Obtaining and characterizing a Taxus baccata tetraploid cell line

Novel approaches to optimize the production of plant specialized metabolites are crucial to reach maximum productivity of plant biofactories. Plant polyploidization frequently enhances protein synthesis and thereby increases the biosynthesis of specialized metabolites. Paclitaxel is a valuable anticancer agent scarcely produced in nature. Therefore, plant biofactories represent a sustainable alternative source of this compound and related taxanes. With the aim of improving the productivity of Taxus spp. cell cultures, we induced polyploidy in vitro by treating immature embryos of Taxus baccata with colchicine. To obtain the polyploid cell lines, calli were induced from T. baccata plantlets previously treated with colchicine and ploidy levels were accurately identified using flow cytometry. In terms of cell morphology, tetraploid cells were about 3-fold bigger than the diploid cells. The expression of taxane pathway genes was higher in the tetraploid cell line compared to the diploid cells. Moreover, taxane production was 6.2-fold higher and the production peak was achieved 8 days earlier than in the diploid cell line, indicating a higher productivity. The obtained tetraploid cell line proved to be highly productive, constituting a step forward towards the development of a bio-sustainable production system for this chemotherapeutic drug.

Metabolic Effects of Elicitors on the Biosynthesis of Tropane Alkaloids in Medicinal Plants

Tropane alkaloids (TAs) are large secondary metabolite alkaloids that find extensive applications in the synthesis of antidotes, anesthetics, antiemetics, motion sickness drugs, and antispasmodics. The current production method primarily depends on extraction from medicinal plants of the Solanaceae family. Elicitation, as a highly effective biotechnological approach, offers significant advantages in augmenting the synthesis of secondary metabolites. The advantages include its simplicity of operation, low cost, and reduced risk of contamination. This review focuses on the impact of elicitation on the biosynthesis of TAs from three aspects: single-elicitor treatment, multiple-elicitor treatment, and the combination of elicitation strategy with other strategies. Some potential reasons are also proposed. Plant hormones and growth regulators, such as jasmonic acid (JA), salicylic acid (SA), and their derivatives, have been extensively employed in the separate elicitation processes. In recent years, novel elicitors represented by magnetic nanoparticles have emerged as significant factors in the investigation of yield enhancement in TAs. This approach shows promising potential for further development. The current utilization of multi-elicitor treatment is constrained, primarily relying on the combination of only two elicitors for induction. Some of these combinations have been found to exhibit synergistic amplification effects. However, the underlying molecular mechanism responsible for this phenomenon remains largely unknown. The literature concerning the integration of elicitation strategy with other strategies is limited, and several research gaps require further investigation. In conclusion, the impact of various elicitors on the accumulation of TAs is well-documented. However, further research is necessary to effectively implement elicitation strategies in commercial production. This includes the development of stable bioreactors, the elucidation of regulatory mechanisms, and the identification of more potent elicitors.

Cell Mutagenic Autopolyploidy Enhances Salinity Stress Tolerance in Leguminous Crops

Salinity stress affects plant growth and development by causing osmotic stress and nutrient imbalances through excess Na+, K+, and Cl- ion accumulations that induce toxic effects during germination, seedling development, vegetative growth, flowering, and fruit set. However, the effects of salt stress on growth and development processes, especially in polyploidized leguminous plants, remain unexplored and scantly reported compared to their diploid counterparts. This paper discusses the physiological and molecular response of legumes towards salinity stress-based osmotic and ionic imbalances in plant cells. A multigenic response involving various compatible solutes, osmolytes, ROS, polyamines, and antioxidant activity, together with genes encoding proteins involved in the signal transduction, regulation, and response mechanisms to this stress, were identified and discussed. This discussion reaffirms polyploidization as the driving force in plant evolution and adaptation to environmental stress constraints such as drought, feverish temperatures, and, in particular, salt stress. As a result, thorough physiological and molecular elucidation of the role of gene duplication through induced autopolyploidization and possible mechanisms regulating salinity stress tolerance in grain legumes must be further studied.

Cultivar-dependent phenotypic and chemotypic responses of drug-type Cannabis sativa L. to polyploidization

Cannabis sativa L. is a plant with a wide range of potential medicinal applications. In recent years, polyploidy has gained attention as a potential strategy for rapidly improving C. sativa, which, unlike other modern crops, has not yet benefitted from this established biotechnological application. Currently, no reports on high THCA and CBDA drug-type polyploid cultivars have been published. Moreover, it still needs to be clarified if different cultivars react similarly to polyploidization. For these reasons, we set out to evaluate and compare the phenotype and chemotype of three high Δ9-tetrahydrocannabinolic acid (THCA) and one high cannabidiolic acid (CBDA) drug-type cultivars in their diploid, triploid and tetraploid state through agronomic and metabolomic approaches. Our observations on plant morphology revealed a significant increase in plant height and leaf size with increasing ploidy levels in a cultivar-dependent manner. In contrast, cannabinoids were negatively affected by polyploidization, with the concentration of total cannabinoids, THCA, CBDA and cannabigerolic acid (CBGA) decreasing significantly in higher ploidy levels across all four cultivars. Headspace analysis of volatiles revealed that total volatile content decreased in triploids. On the other hand, tetraploids reacted differently depending on the cultivars. Two THCA dominant cultivars showed an increase in concentrations, while in the other two cultivars, concentrations decreased. Additionally, several rare compounds not present in diploids appeared in higher ploidy levels. Moreover, in one high THCA cultivar, a couple of elite tetraploid genotypes for cannabinoid and volatile production were identified, highlighting the role of cultivar and genotypic variability as an important factor in Cannabis sativa L. polyploids. Overall, our observations on plant morphology align with the giga phenotype observed in polyploids of other plant species. The decrease in cannabinoids and volatiles production in triploids have relevant implications regarding their commercial use. On the other hand, this study found that tetraploidization is a suitable approach to improve Cannabis sativa L. medicinal potential, although the response is cultivar and genotype-dependent. This work lays the ground for further improving, evaluating and harnessing Cannabis sativa L. chemical diversity by the breeding, biotechnological and pharmaceutical sectors.

Anti-AChE and Anti-BuChE Screening of the Fermentation Broth Extracts from Twelve Aspergillus Isolates and GC-MS and Molecular Docking Studies of the Most Active Extracts

Nowadays, the administration of cholinesterase enzyme (acetylcholinesterase: AChE and butyrylcholinesterase: BuChE) inhibitors is very common for the symptomatic treatment of Alzheimer's disease and the other forms of dementia and CNS disorders. In this paper, the anti-AChE and anti-BuChE activities of the fermentation broth ethyl acetate extracts from twelve Aspergillus isolates were evaluated by Ellman method. The results showed that A1 (Aspergillus flavus) and A5 (Aspergillus tubingensis, isolate 1) extracts with IC50 values of 46.77 μg/mL and 75.85 μg/mL possess the greatest ability to inhibit AChE and BuChE, respectively. GC-MS analysis of the extracts (A1 and A5) demonstrated that two alkaloids named 14-methyl-16-azabicyclo[10.3.1]hexadeca-1(15),12(16),13-triene (MAHT) and 6-chloro-2-methyl-7,8,9,10-tetrahydro-phenanthridine (CMTP) account for the highest percentage of A1 (26.95%) and A5 (25.5%) extracts, respectively. A 2-pyrazoline derivative, 5-hydroxy-3-(4-pyridinyl)-5-trifluoromethyl-1-(2,4,6-trimethylphenoxyacetyl)- (PHPTT), also constituted the high percentage (9.54%) of A5 extract. The anticholinesterase and neuroprotective effects of some 2-pyrazoline derivatives have been previously reported. The interaction study of MAHT with human AChE and CMTP and PHPTT with human BuChE using molecular docking indicated that these alkaloids bind to the active site gorge of the enzymes with high affinity. The best docking scores of MAHT, CMTP, and PHPTT were -7.1, -8.2, and -9.7 kcal/mol, respectively.

The known, unknown, and the intriguing about members of a critically endangered traditional medicinal plant genus Aconitum

Humanity will always be indebted to plants. In the ongoing scientific era, the 'Herbal Revolution' has helped discover several valuable medicinal plants and associated novel secondary metabolites from the diverse unexplored ecosystems, treating several diseases via phytotherapy. The Aconitum genus comprises several economically-important poisonous mountainous medicinal plant species whose unique biodiversity is on the verge of extinction due to illegal human intervention triggered habitat loss, over-harvesting, and unrestricted trading. Owing to its vast diversity of diterpene alkaloids, most species are extensively used to treat several ailments in rural parts of the world. Irrespective of this, many unexplored and intriguing prospects exist to understand and utilize this critical plant for human benefit. This systematic review tries to fill this gap by compiling information from the sporadically available literature known for ~300 Aconitum spp. regarding its nomenclature and classification, endangerment, plant morphology, ploidy, secondary metabolites, drug pharmacokinetics, conservation, and omics-based computational studies. We also depicted the disparity in the studied model organisms for this diverse genus. The absence of genomic/metagenomic data is becoming a limiting factor in understanding its plant physiology, metabolic pathways, and plant-microbes interactions, and therefore must be promoted. Additionally, government support and public participation are crucial in establishing conservation protocols to save this plant from endangerment.

Ecophysiological Differentiation among Two Resurrection Ferns and Their Allopolyploid Derivative

Theoretically, the coexistence of diploids and related polyploids is constrained by reproductive and competitive mechanisms. Although niche differentiation can explain the commonly observed co-occurrence of cytotypes, the underlying ecophysiological differentiation among cytotypes has hardly been studied. We compared the leaf functional traits of the allotetraploid resurrection fern Oeosporangium tinaei (HHPP) and its diploid parents, O. hispanicum (HH) and O. pteridioides (PP), coexisting in the same location. Our experimental results showed that all three species can recover physiological status after severe leaf dehydration, which confirms their 'resurrection' ability. However, compared with PP, HH had much higher investment per unit area of light-capturing surface, lower carbon assimilation rate per unit mass for the same midday water potential, higher non-enzymatic antioxidant capacity, higher carbon content, and lower contents of nitrogen, phosphorus, and other macronutrients. These traits allow HH to live in microhabitats with less availability of water and nutrients (rock crevices) and to have a greater capacity for resurrection. The higher assimilation capacity and lower antioxidant capacity of PP explain its more humid and nutrient-rich microhabitats (shallow soils). HHPP traits were mostly intermediate between those of HH and PP, and they allow the allotetraploid to occupy the free niche space left by the diploids.

Morphological, cariological, and phytochemical studies of diploid and autotetraploid Hippeastrum papilio plants

The polyploidization of Hippeastrum papilio influences its primary and secondary metabolism including the biosynthesis of bioactive alkaloids. Hippeastrum papilio is an ornamental plant that has advantages in comparison to the currently used plants for the extraction of galanthamine, a natural compound used for the cognitive treatment of Alzheimer's disease. In the present study, an autotetraploid line of H. papilio was induced for the first time, after treatment with 0.05% colchicine for 48 h. The chromosome number in diploids was found to be 2n = 2x = 22 and for autotetraploids 2n = 4x = 44. The flow cytometric analyses detected a DNA C-value of 14.88 ± 0.03 pg (1C) in diploids and 26.57 ± 0.12 pg in autotetraploids. The morphological, cytological, and phytochemical studies showed significant differences between diploids and autotetraploids. The length and width of stomata in autotetraploids were 22.47% and 17.94%, respectively, larger than those observed in the diploid leaves. The biomass of one-year-old autotetraploid H. papilio plants was reduced by 53.99% for plants' fresh weight, 56.53% for leaves' fresh weight, and 21.70% for bulb diameter. The GC-MS analysis of methanol extracts from one-year-old diploid and autotetraploid H. papilio plants revealed over 60 primary and secondary metabolites including alkaloids, phenolic acids, sterols, saccharides, and alcohols, among others. Principal component analysis of the metabolite profiles indicates a divergence of the metabolism between diploid and autotetraploid plants. The content of galanthamine and haemanthamine was found to be 49.73% and 80.10%, respectively, higher in the leaves of autotetraploids, compared to the diploid ones. The biosynthesis of the saccharides shows a tendency to be upregulated in tetraploid plants, while that of phenolic acids was downregulated. Polyploidization of H. papilio creates possibilities for further crop improvement aimed at high-galanthamine-producing genotypes.

Comparative transcriptome and metabolome analyses of four Panax species explore the dynamics of metabolite biosynthesis

Background

The genus Panax in the Araliaceae family has been used as traditional medicinal plants worldwide and is known to biosynthesize ginsenosides and phytosterols. However, genetic variation between Panax species has influenced their biosynthetic pathways is not fully understood.

Methods

Simultaneous analysis of transcriptomes and metabolomes obtained from adventitious roots of two tetraploid species (Panax ginseng and P. quinquefolius) and two diploid species (P. notoginseng and P. vietnamensis) revealed the diversity of their metabolites and related gene expression profiles.

Results

The transcriptome analysis showed that 2,3-OXIDOSQUALENE CYCLASEs (OSCs) involved in phytosterol biosynthesis are upregulated in the diploid species, while the expression of OSCs contributing to ginsenoside biosynthesis is higher in the tetraploid species. In agreement with these results, the contents of dammarenediol-type ginsenosides were higher in the tetraploid species relative to the diploid species.

Conclusion

These results suggest that a whole-genome duplication event has influenced the triterpene biosynthesis pathway in tetraploid Panax species during their evolution or ecological adaptation. This study provides a basis for further efforts to explore the genetic variation of the Panax genus.

Rhamnosyltransferases involved in the biosynthesis of flavone rutinosides in Chrysanthemum species

Linarin (acacetin-7-O-rutinoside), isorhoifolin (apigenin-7-O-rutinoside), and diosmin (diosmetin-7-O-rutinoside) are chemically and structurally similar flavone rutinoside (FR) compounds found in Chrysanthemum L. (Anthemideae, Asteraceae) plants. However, their biosynthetic pathways remain largely unknown. In this study, we cloned and compared FRs and genes encoding rhamnosyltransferases (RhaTs) among eight accessions of Chrysanthemum polyploids. We also biochemically characterized RhaTs of Chrysanthemum plants and Citrus (Citrus sinensis and Citrus maxima). RhaTs from these two genera are substrate-promiscuous enzymes catalyzing the rhamnosylation of flavones, flavanones, and flavonols. Substrate specificity analysis revealed that Chrysanthemum 1,6RhaTs preferred flavone glucosides (e.g. acacetin-7-O-glucoside), whereas Cs1,6RhaT preferred flavanone glucosides. The nonsynonymous substitutions of RhaTs found in some cytotypes of diploids resulted in the loss of catalytic function. Phylogenetic analysis and specialized pathways responsible for the biosynthesis of major flavonoids in Chrysanthemum and Citrus revealed that rhamnosylation activity might share a common evolutionary origin. Overexpression of RhaT in hairy roots resulted in 13-, 2-, and 5-fold increases in linarin, isorhoifolin, and diosmin contents, respectively, indicating that RhaT is mainly involved in the biosynthesis of linarin. Our findings not only suggest that the substrate promiscuity of RhaTs contributes to the diversity of FRs in Chrysanthemum species but also shed light on the evolution of flavone and flavanone rutinosides in distant taxa.

Validation of Reference Genes in a Population of Blueberry (Vaccinium corymbosum) Plants Regenerated in Colchicine

For the first time we report the validation of reference genes in plants from a population of blueberry (Vaccinium corymbosum) clones cultured in vitro on a colchicine-supplemented medium. Nodal segment explants of the cultivar Duke were regenerated by organogenesis under different periods of colchicine 1 mg/L exposure (1, 2, 3, 5, 30 days). The clones selected for the study showed variability for phenotypic traits after 2 years of adaptation to field conditions, compared to plants of the donor genotype that were regenerated on a medium without colchicine. Vaccinium myrtillus (GAPDH) and Vaccinium macrocarpon (ATP1, NADH, RPOB and COX2) were used as reference genomes for primer design. The results show that colchicine treatments can cause genomic changes in blueberry plants. At the molecular level, exposure of plants to colchicine in early periods could promote an increase in gene expression of specific genes such as ATP1, COX2, GAPDH, MATK, NADH and RPOB. However, prolonged exposure (30 days) could decrease gene expression of the genes studied. For qPCR assays, the primers designed for ATP1, COX2, GAPDH and MATK genes showed high efficiency. In addition, the GAPDH, ATP1, NADH and COX2 genes showed high stability and could be recommended as potential reference genes for gene expression assays in Vaccinium.