- Research article
- Open Access
The chemotaxonomic classification of Rhodiola plants and its correlation with morphological characteristics and genetic taxonomy
- Zhenli Liu†1,
- Yuanyan Liu†2,
- Chunsheng Liu2,
- Zhiqian Song1,
- Qing Li1,
- Qinglin Zha3,
- Cheng Lu3,
- Chun Wang1,
- Zhangchi Ning2,
- Yuxin Zhang2,
- Cheng Tian2 and
- Aiping Lu3, 4Email author
© Liu et al.; licensee Chemistry Central Ltd. 2013
- Received: 18 March 2013
- Accepted: 10 July 2013
- Published: 12 July 2013
Rhodiola plants are used as a natural remedy in the western world and as a traditional herbal medicine in China, and are valued for their ability to enhance human resistance to stress or fatigue and to promote longevity. Due to the morphological similarities among different species, the identification of the genus remains somewhat controversial, which may affect their safety and effectiveness in clinical use.
In this paper, 47 Rhodiola samples of seven species were collected from thirteen local provinces of China. They were identified by their morphological characteristics and genetic and phytochemical taxonomies. Eight bioactive chemotaxonomic markers from four chemical classes (phenylpropanoids, phenylethanol derivatives, flavonoids and phenolic acids) were determined to evaluate and distinguish the chemotaxonomy of Rhodiola samples using an HPLC-DAD/UV method. Hierarchical cluster analysis (HCA) and principal component analysis (PCA) were applied to compare the two classification methods between genetic and phytochemical taxonomy.
The established chemotaxonomic classification could be effectively used for Rhodiola species identification.
- Rhodiola plants
- Morphological characteristic
- Genetic taxonomy
- Phytochemical taxonomy
The genus Rhodiola L. (Crassulaceae) comprises approximately 96 species found in the alpine regions of Asia and Europe. A total of 73 species, 2 subspecies and 7 varieties are found in China [1, 2]. Rhodiola species, historically used as adaptogens in Russia and northern Europe and as a traditional herbal medicine in China, are valued for their ability to enhance human resistance to stress or fatigue and to promote longevity [3–5]. Rhodiola plants are mainly distributed in southwest and northwest of China, with most species located in Tibet and in Sichuan province. In China, the Rhodiola species called Hongjingtian have been used as an important adaptogen, hemostatic, and tonic in traditional Tibetan medicines for thousands of years . The phytochemical extracts of Rhodiola plants are widely used throughout Europe, Asia and the United States, with biological activities including anti-allergenic and anti-inflammatory effects and enhanced mental alertness, as well as a variety of other therapeutic applications . Because of their commercial utility, Rhodiola plants are now cultivated in many locations in Europe and Asia. Most notably, the roots and rhizomes of R. crenulata (RC) have high activities and have been accepted by the Pharmacopoeia of China . In addition, many Rhodiola plants, such as R. sachalinensis (RS), R. himalensis (D. Dons) S. H. Fu (RH), R. serrata H. Ohba (RSE), R. rosea L. (RR), R. kirilowii (Regel) Maxim (RK) and R. fastigiata (HK. F. et Thoma) S. H. Fu (RF), etc., are also used as Hongjingtian in China. However, the identification of the closely related species of Rhodiola plants is often difficult due to their generally similar morphology.
Phytochemical investigations show that there are six important classes of constituents in Rhodiola rhizomes, including phenylpropanoids, phenylethanol derivatives, flavonoids, monoterpernes, triterpenes and phenolic acids [8–10]. Using animal models, bioassay-guided fractionation of various extracts of plant adaptogens have shown that the active components are mainly phenylpropanoids and phenylethanol derivatives, including salidroside, rosavin and tyrosol [4, 11–15]. In the Chinese Pharmacopoeia, salidroside is chosen as a marker compound for quality control . Phenylpropanoids, such as rosarin, rosavin and rosin, are not only typical for Rhodiola rhizomes but are also pharmacologically active as antioxidants and neuro-stimulants [16–18]. Compounds such as tyrosol and gallic acid had been proven to be good radical scavengers [19, 20]. A recent study revealed that rhodionin and salidroside might have anti-tumor effects [21, 22]. Additionally, rhodionin is recognized to be involved in learning and memory [23, 24]. The above 8 compounds from four chemical classes were selected as chemotaxonomic markers in the present paper. Because they exhibit variety bioactivities, meanwhile, phenylpropanoids and phenylethanol derivatives are characteristic in Rhodiola plants. Therefore, it is worthwhile to study the variety of co-existing phytochemical constituents in the plant, which may be responsible for its unique pharmacological activity.
The current taxonomical status of the genus Rhodiola has become quite complex. The rationale and defining criteria for the boundaries of the genus remain somewhat controversial . The morphologies of different species of commercial Rhodiola plants are too similar to distinguish visually. With the development of DNA sequencing methods and the discovery of the polymerase chain reaction (PCR) for DNA amplification, biological systematic analysis has increasingly been based on DNA sequence analysis. The genotypes identified by PCR amplification suffice to predict the species of plants. In addition to genetic taxonomy and other classical morphological and non-morphological methods, phytochemical taxonomy can also provide supplementary information in species identification . The chemotype of a plant species has traditionally been defined as by profile of natural products, and the genotype has been defined as its genetic constitution or DNA sequence.
Here, we collected 47 Rhodiola samples used clinically in different provinces of China. They were identified through morphological characterization and genotyping. In addition, eight bioactive compounds from four chemical classes (phenylpropanoids, phenylethanol derivatives, flavonoids and phenolic acids) were used as chemotaxonomic markers to elucidate the phytochemical taxonomy by an HPLC-DAD/UV method.
Chemicals and materials
The origins of the 47 Rhodiola samples
The data were obtained using an Agilent 1100 Series HPLC with DAD. The analytical conditions for recording chromatograms of the marker compounds in Rhodiola samples were as follows. A Zobax SB-C18 column (4.6 mm× 150 mm, 5 μm; Agilent Technologies) was used. The mobile phase consisted of MeCN (A) and 0.2% HAc (B) with a linear gradient elution at a flow rate of 1.0 mL/min. The gradient program (A/B, v/v) was as follows: 5:95 (t = 0 min), 8:92 (t = 10 min), 18:82 (t = 43 min) and 37:63 (t = 60 min). The detection wavelength program was 275 nm (t = 0 min), 250 nm (t = 30 min) and 332 nm (t = 45 min). The column temperature was set to 40°C. The detection wavelength was selected by DAD according to max UV absorption of each reference.
The samples were pulverized, and the powder (1.0 g) was accurately weighed and extracted with 25 mL of methanol by ultrasonication for 30 min. After cooling, the solution was filtered through a 0.2 μm membrane filter and stored at 4°C until analysis. A 5-μL aliquot solution was injected for HPLC analysis. Each sample was prepared in triplicate and relative standard deviation (RSD) was calculated for all the samples.
The eight standards were prepared and serially diluted with methanol to obtain seven different concentrations used for plotting standard curves, respectively. Method precision was determined by injecting one Hongjingtian sample solution six times consecutively. Reproducibility was studied through six independently prepared samples from a single batch of Hongjingtian. The stability test was performed by successively injecting the same sample solution over 24 hours. The limit of detection (LOD) and limit of quantity (LOQ) were determined at a signal-to-noise ratio (S/N) of 3 and 10, respectively. Standard solutions were diluted to series of appropriate concentrations with methanol and 5μL aliquots of the diluted solutions were injected into the HPLC for analysis.
Nucleic acids were extracted and purified from deep-frozen plant materials. Sample vouchers were deposited in the collection at Beijing University of Chinese Medicine. Genomic DNA was extracted from silica gel-dried root material using a plant DNA extraction kit (Tiangen, Beijing, China) according to the manufacturer’s protocol with some modifications. The quality of the isolated DNA was verified from absorbance measurements at wavelengths 230, 260 and 280 nm and on a 1% (w/v) ethidium bromide-stained agarose gel.
PCR, cloning and sequencing
The PCR procedure was designed according to the instruction manual of a GeneRacer Kit (Invitrogen, Carlsbad, CA, USA). To achieve the 5′-end cDNA sequence, two rounds of thermal asymmetric interlaced PCR were performed as described in Liu and Chen et al. . All of the obtained fragments were sequenced in both directions by ligating into pGEM&nonBR;T vector and using an ABI 3730XL Genetic Analyzer (Applied Biosystems). PCR was performed using a touchdown strategy: 94°C for 4 min, followed by 10 cycles of 94°C for 75s, 53°C for 5 min, 0.2°C/s to 41°C, and 72°C for 5 min, followed by 35cycles of 94°C for 1 min, 45°C for 2 min, and 72°C for 5 min. The PCR products were run on a 1% (w/v) ethidium bromide-stained agarose gel with a 6 × orange loading buffer (Fermentas, Vilnius, Lithuania). The expected size band (680 bp) was excised from the gel and eluted using a Qiaquick Gel Extraction kit from Qiagen. The eluted PCR product was cloned into pGEM-T Easy Vector (Promega Corporation, Madison, USA) and sequenced using the BigDye Terminator Cycle Sequencing Kit (PE Applied Biosystems, Warrington, UK).The full-length deduced amino acid sequence was aligned with the publicly available HQT groups using ClustalX and MEGA version 4.0 software. A neighbor-joining (NJ) tree was constructed based on standard parameters with bootstrap testing of 1000 replicates. All the DNA sequences obtained were submitted to GenBank-NCBI for comparison with the deposited sequences using the tool BLAST .
Statistical and multivariate analysis
The statistical analysis was performed using the SAS 9.1.3 statistical package (order no. 195557) for PCA and HCA. PCA and HCA were used to show the unsupervised clustering pattern of the Rhodiola species. PCA and HCA were used to observe the natural interrelationships among the chemical components for each of the Rhodiola samples. The critical p value for all analyses in this study was set to 0.05.
Morphological characteristic of the collected samples
The 47 collected Rhodiola samples were identified according to their morphological characteristics, and their collection locations are listed in Table 1. According to Flora of China , the morphological characteristic of Rhodiola plants is as following: stems dimorphic with usually very stout caudex or rhizome, usually with brown or blackish, membranous, scalelike leaves, sharply differentiated from much more slender, erect or ascending, leafy flowering stems. The roots and rhizomes of Rhodiola are used as the medicinal parts of the plants. The supplementary characteristic information from the aerial parts of the collected samples could not be obtained in this experiment. The morphologies of some Rhodiola samples are too similar to distinguish visually. In addition, the current taxonomical status of the genus Rhodiola has become quite complex. Accordingly, the species of some samples were tentatively identified. Among the 47 Rhodiola samples, RC-1 ~ RC-19 were identified as R. crenulata; RS-20 ~ RS-21 were identified as R. sachalinensis; RH-22 was identified as R. himalensis (D. Dons) S. H. Fu; RSE-23 ~ RSE-37 were identified as R. serrata H. Ohba; RF-38 was identified as R. fastigiata (HK. F. et Thoma) S. H. Fu; RR-39 was identified as R. rosea L.; RK-40 ~ RK-47 were identified as R. kirilowii (Regel) Maxim, as shown in Table 1.
Genetic taxonomic of the collected samples
To verify the accuracy of the identity of the species of the 47 Rhodiola samples, a phylogenetic tree was constructed based on the DNA sequences obtained from GenBank. For some samples, it was difficult to extract the exact DNA when the roots had been stored too long and were very dry, even when the roots were contaminated by microbes or the DNA was degraded. Therefore, only 34 Rhodiola samples were identified according to their DNA sequences. Among them, 23 Rhodiola samples were unambiguously identified with the similarities higher than 98%. The other 11 Rhodiola samples were tentatively identified. All sequences were submitted to GenBank (http://www.ncbi.nlm.nih.gov/genbank/) and their IDs were shown in the Additional file 1.
Phytochemical taxonomy of the collected samples
Based on the morphological characteristic and genetic taxonomy, both of them have some limitations in species identification. The classification of plants based on chemotypes can be used as a powerful chemotaxonomic tool that provides a detailed view of the differences and similarities between species. The 8 pure bioactive compounds classified into four types were used as chemotaxonomic markers to distinguish among the different Rhodiola samples.
To gain high sensitivity and good peak capacity, the chromatographic conditions were optimized, as described in the HPLC-DAD/UV analysis section. MeCN and 0.2%HAc were used as the mobile phase to improve the retention behavior of the constituents on the HPLC column. The wavelength for the detection of compounds was selected by DAD. The chromatograms at 275 nm could provide maximum absorption of gallic acid, tyrosol, salidroside and (+) catechin. The wave length for the detection of rosarin, rosavin and rosin was at 250 nm, and 332 nm was used for the detection of rhodionin.
Linear regression, LODs and LOQs, precisions, reproducibility, stability and recovery for eight compounds
Recover (n = 6)
(n = 2)
(n = 6) RSD (%)
(n = 6) RSD (%)
Y = 218.51X-0.4
Y = 384.13X-3.9
Y = 2923.6X-13.8
Y = 339.77X + 3.7
Y = 1954.5X-2.7
Y = 2140.4X + 8.4
Y = 3076.8X-4.2
Y = 166.15X + 1.3
Accordingly, the 23 Rhodiola samples classified by genetic taxonomy were analyzed by chemotaxonomic classification using the four types of bioactive compound as reference markers. HCA showed that Rhodiola samples were divided into five branches according to their chemotaxonomy (Figure 2b). Of the two classification methods for classifying the genus Rhodiola samples, HCA results showed considerably comparable results for both the genotype- and chemotype-based classification methods. This observation may be because different genotypes caused different chemotypes due to the genotype-dependent production of metabolites.
In this study, the 47 Rhodiola samples used commercially in China were identified by their morphological characteristics and genetic and phytochemical taxonomies. In the morphological characteristics, there exists variation between populations of the same species at different life stages and from different environments . If the samples collected are not intact, the accuracy of the identification will be affected. The morphologies of some Rhodiola samples are too similar to distinguish visually. Genetic taxonomy can provide exact classification of species submitted to GenBank. However, there may be uncertainty in the extraction of pure DNA from every sample, and the procedure of analyses is time-consuming. Here, eight bioactive compounds from four chemical classes (phenylpropanoids, phenylethanol derivatives, flavonoids and phenolic acids) were used as chemotaxonomic markers to evaluate and distinguish the chemotypes of 47 Rhodiola samples by an HPLC-DAD/UV method. First, 23 Rhodiola samples classified by genetic taxonomy and morphological characteristics were analyzed by chemotaxonomic classification, which showed considerably comparable results. This analysis indicated that different genotypes caused different chemotypes due to the genotype-dependent production of metabolites. Next, all the 47 Rhodiola samples were analyzed by PCA and HCA based on the content of the eight bioactive references. All the samples were divided into four clusters according to the established phytochemical taxonomic method. Consequently, chemotyping became useful for distinguishing morphologically similar species, by identifying variants of the chemotaxonomic markers. However, certain limitations exist in the present studies. The number of collected samples from RS, RH, RF and RR species is low, which affected the accuracy of this phytochemical taxonomic classification. Moreover, the clusters of genes for biosynthesis regulating the metabolite profiling need to be investigated in future.
This study was financially supported by the National Science Foundation of China (Project No. 30772726, No. 30825047, No. 30902000 and No. 81001623).
- Li T, Zhang H: Identification and comparative determination of rhodionin in traditional tibetan medicinal plants of fourteen Rhodiola species by high-performance liquid chromatography-photodiode array detection and electrospray ionization-mass spectrometry. Chem Pharm Bull(Tokyo). 2008, 56: 807-814. 10.1248/cpb.56.807.View ArticleGoogle Scholar
- Fu SX: Flora of China. 1984, Beijing: Science PressGoogle Scholar
- Mattioli L, Titomanlio F, Perfumi M: Effects of a Rhodiola rosea L. extract on the acquisition, expression, extinction, and reinstatement of morphine-induced conditioned place preference in mice. Psychopharmacology. 2012, 221: 183-193. 10.1007/s00213-012-2686-0.View ArticleGoogle Scholar
- Mattioli L, Funari C, Perfumi M: Effects of Rhodiola rosea L. extract on behavioural and physiological alterations induced by chronic mild stress in female rats. J Psychopharmacol. 2009, 23: 130-142.View ArticleGoogle Scholar
- Tolonen A, Hohtola A, Jalonen J: Comparison of electrospray ionization and atmospheric pressure chemical ionization techniques in the analysis of the main constituents from Rhodiola rosea extracts by liquid chromatography/mass spectrometry. J Mass Spectrom. 2003, 38: 845-853. 10.1002/jms.497.View ArticleGoogle Scholar
- Yang YCHTN, Lu SL, Hung RF, Wang ZX: Zang Yao Zhi. 1991, Xining: Qinghai People’s Publishing HouseGoogle Scholar
- Pharmacopoeia: Pharmacopoeia of the People’s Republic of China. 2010, Beijing: Chemical Industry Press, 144Google Scholar
- Ali Z, Fronczek FR, Khan IA: Phenylalkanoids and monoterpene analogues from the roots of Rhodiola rosea. Planta Med. 2008, 74: 178-181. 10.1055/s-2008-1034288.View ArticleGoogle Scholar
- Panossian AWG, Sarris J: Rosenroot (Rhodiola rosea): traditional use, chemical composition, pharmacology and clinical efficacy. Phytomedicine. 2010, 7: 481-493.View ArticleGoogle Scholar
- Hohtola A: Bioactive compounds from northern plants. Adv Exp Med Biol. 2010, 698: 99-109. 10.1007/978-1-4419-7347-4_8.View ArticleGoogle Scholar
- Kurkin VA, Zapesochnaya GG: Chemical composition and pharmacological properties of Rhodiola rosea L. Chemical and Pharmaceutical Journal. 1986, 20: 1231-1244.Google Scholar
- Zapesochnaya GG, Kurkin VA, Boyko VP, Kolkhir VK: Phenylpropanoids promising biologically active compounds of medicinal plants. Khim Farm Zh. 1995, 29: 47-50.Google Scholar
- Peschel W, Prieto JM, Karkour C, Williamson EM: Effect of provenance, plant part and processing on extract profiles from cultivated European Rhodiola rosea L. for medicinal use. Phytochemistry. 2012, 86: 92-102.View ArticleGoogle Scholar
- Cifani C, Micioni Di BM, Vitale G, Ruggieri V, Ciccocioppo R, Massi M: Effect of salidroside, active principle of Rhodiola rosea extract, on binge eating. Physiol Behav. 2010, 101: 555-562. 10.1016/j.physbeh.2010.09.006.View ArticleGoogle Scholar
- Wiedenfeld H, Dumaa M, Malinowski M, Furmanowa M, Narantuya S: Phytochemical and analytical studies of extracts from Rhodiola rosea and Rhodiola quadrifida. Pharmazie. 2007, 62: 308-311.Google Scholar
- Huang SC, Lee FT, Kuo TY, Yang JH, Chien CT: Attenuation of long-term Rhodiola rosea supplementation on exhaustive swimming-evoked oxidative stress in the rat. Chin J Physiol. 2009, 52: 316-324. 10.4077/CJP.2009.AMH029.View ArticleGoogle Scholar
- Schriner SE, Abrahamyan A, Avanessian A, Bussel I, Maler S, Gazarian M, Holmbeck MA, Jafari M: Decreased mitochondrial superoxide levels and enhanced protection against paraquat in Drosophila melanogaster supplemented with Rhodiola rosea. Free Radic Res. 2009, 43: 836-843. 10.1080/10715760903089724.View ArticleGoogle Scholar
- Kucinskaite A, Briedis V, Savickas A: Experimental analysis of therapeutic properties of Rhodiola rosea L. and its possible application in medicine. Medicina (Kaunas). 2004, 40: 614-619.Google Scholar
- Sun L, Isaak CK, Zhou Y, Petkau JC OK, Liu Y, Siow YL: Salidroside and tyrosol from Rhodiola protect H9c2 cells from ischemia/reperfusion-induced apoptosis. Life Sci. 2012, 91: 151-158. 10.1016/j.lfs.2012.06.026.View ArticleGoogle Scholar
- Lee MW, Lee YA, Park HM, Toh SH, Lee EJ, Jang HD, Kim YH: Antioxidative phenolic compounds from the roots of Rhodiola sachalinensis A. Bor. Arch Pharm Res. 2000, 23: 455-458. 10.1007/BF02976571.View ArticleGoogle Scholar
- Sun C, Wang Z, Zheng Q, Zhang H: Salidroside inhibits migration and invasion of human fibrosarcoma HT1080 cells. Phytomedicine. 2012, 19: 355-363. 10.1016/j.phymed.2011.09.070.View ArticleGoogle Scholar
- Hu X, Lin S, Yu D, Qiu S, Zhang X, Mei R: A preliminary study: the anti-proliferation effect of salidroside on different human cancer cell lines. Cell Biol Toxicol. 2010, 26: 499-507. 10.1007/s10565-010-9159-1.View ArticleGoogle Scholar
- Choe KI, Kwon JH, Park KH, Oh MH, Kim MH, Kim HH, Cho SH, Chung EK, Ha SY, Lee MW: The antioxidant and anti-inflammatory effects of phenolic compounds isolated from the root of rhodiola sachalinensis a. BOR. Molecules. 2012, 17: 11484-11494. 10.3390/molecules171011484.View ArticleGoogle Scholar
- Kobayashi K, Yamada K, Murata T, Hasegawa T, Takano F, Koga K, Fushiya S, Batkhuu J, Yoshizaki F: Constituents of Rhodiola rosea showing inhibitory effect on lipase activity in mouse plasma and alimentary canal. Planta Med. 2008, 74: 1716-1719. 10.1055/s-0028-1088318.View ArticleGoogle Scholar
- de Kuppler AL, Steiner U, Sulyok M, Krska R, Oerke EC: Genotyping and phenotyping of Fusarium graminearum isolates from Germany related to their mycotoxin biosynthesis. Int J Food Microbiol. 2011, 151: 78-86. 10.1016/j.ijfoodmicro.2011.08.006.View ArticleGoogle Scholar
- Desjardins AE: Natural product chemistry meets genetics: when is a genotype a chemotype. J Agric Food Chem. 2008, 56: 7587-7592. 10.1021/jf801239j.View ArticleGoogle Scholar
- Liu YG, Chen Y: High-efficiency thermal asymmetric interlaced PCR for amplification of unknown flanking sequences. Biotechniques. 2007, 43: 649-650. 10.2144/000112601. 652, 654 passimView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.View ArticleGoogle Scholar
- Sandasi M, Kamatou GP, Viljoen AM: An untargeted metabolomic approach in the chemotaxonomic assessment of two Salvia species as a potential source of alpha-bisabolol. Phytochemistry. 2012, 84: 94-101.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.