Search

Phenylethanoid glycosides (PhGs) are naturally occurring glycosides derived from plants with various biological activities. Glycosyltransferases catalyze the production of PhGs from phenylethanols via a transglycosylation reaction. The low activity and stability of glycosyltransferase limit its industrial application. An ancestral glycosyltransferase, UGTAn85, with heat resistance, alkali resistance, and high stability was resurrected using ancestral sequence reconstruction technology. This enzyme can efficiently convert phenylethanols to PhGs. The optimal reaction temperature and pH for UGTAn85 were found to be 70 °C and pH 10.0, respectively. This study employed a combination of structure-guided rational design and co-evolution analysis to enhance its catalytic activity. Potential mutation sites were identified through computer-aided design, including homology modeling, molecular docking, Rosetta dock design, molecular dynamics simulation, and co-evolution analysis. By targeted mutagenesis, the UGTAn85 mutant Q23E/N65D exhibited a 2.2-fold increase in enzyme activity (11.85 U/mg) and elevated affinity ( K m = 0.11 mM) for 2-phenylethanol compared to UGTAn85. Following a fed-batch reaction, 36.16 g/L 2-phenylethyl-β-d-glucopyranoside and 51.49 g/L salidroside could be produced within 24 h, respectively. The findings in this study provide a new perspective on enhancing the stability and activity of glycosyltransferases, as well as a potential biocatalyst for the industrial production of PhGs.

Rhodiola, belonging to the Crassulaceae family, is a perennial herbaceous plant genus. There are about 90 Rhodiola species worldwide, some of which have been reported to have medicinal properties. Rhodiola sachalinensis is a perennial medicinal herb within this genus and, in the present study, its chloroplast genome was sequenced, assembled, annotated and compared with 24 other Rhodiola species. The results obtained show that the chloroplast genome of R. sachalinensis is 151 595 bp long and has a CG content of 37.7%. The inverted repeats (IR) region of the Rhodiola chloroplast genome is the most conserved region, with the main differences being observed in the ycf1 and ndhF genes at the IRb-small single copy boundary, and rps19 and trnH genes at the IRa-large single copy boundary. Phylogenetic analysis showed that Rhodiola species form two major clades, and species with recorded medicinal properties, clustered together in one branch except for R. dumulosa. Within the genus, R. sachalinensis is most closely related to Rhodiola rosea, although comparative analyses showed that only R. sachalinensis and Rhodiola subopposita contained the psbZ gene, which encodes a highly conserved protein subunit of the Photosystem II core complex. Overall, the present study contributes to the understanding of the chloroplast genome of Rhodiola species, and provides a theoretical basis for the study of their genetic diversity and possible use as medicinal plants., (© 2024 The Author(s). FEBS Open Bio published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Rhodiola crenulata, a plant of great medicinal value found in cold high-altitude regions, has been excessively exploited due to the difficulty in cultivation. Understanding Rhodiola crenulata's adaptation mechanisms to cold environment can provide a theoretical basis for artificial breeding. Glutathione peroxidases (GPXs), critical enzymes found in plants, play essential roles in antioxidant defense through the ascorbate-glutathione cycle. However, it is unknown whether GPX5 contributes to Rhodiola crenulata's cold tolerance. In this study, we investigated the role of GPX5 in Rhodiola crenulata's cold tolerance mechanisms. By overexpressing Rhodiola crenulata GPX5 (RcGPX5) in yeast and Arabidopsis thaliana, we observed down-regulation of Arabidopsis thaliana GPX5 (AtGPX5) and increased cold tolerance in both organisms. Furthermore, the levels of antioxidants and enzyme activities in the ascorbate-glutathione cycle were elevated, and cold-responsive genes such as AtCBFs and AtCORs were induced. Additionally, RcGPX5 overexpressing lines showed insensitivity to exogenous abscisic acid (ABA), suggesting a negative regulation of the ABA pathway by RcGPX5. RcGPX5 also promoted the expression of several thioredoxin genes in Arabidopsis and interacted with two endogenous genes of Rhodiola crenulata, RcTrx2-3 and RcTrxo1, located in mitochondria and chloroplasts. These findings suggest a significantly different model in Rhodiola crenulata compared to Arabidopsis thaliana, highlighting a complex network involving the function of RcGPX5. Moreover, overexpressing RcGPX5 in Rhodiola crenulata hairy roots positively influenced the salidroside synthesis pathway, enhancing its pharmaceutical value for doxorubicin-induced cardiotoxicity. These results suggested that RcGPX5 might be a key component for Rhodiola crenulata to adapt to cold stress and overexpressing RcGPX5 could enhance the pharmaceutical value of the hairy roots., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Rosavin is the characteristic component of Rhodiola rosea L., an important medicinal plant used widely in the world that has been reported to possess multiple biological activities. However, the endangered status of wild Rhodiola has limited the supply of rosavin. In this work, we successfully engineered an Escherichia coli strain to efficiently produce rosavin as an alternative production method. Firstly, cinnamate: CoA ligase from Hypericum calycinum, cinnamoyl-CoA reductase from Lolium perenne, and uridine diphosphate (UDP)-glycosyltransferase (UGT) from Bacillus subtilis (Bs-YjiC) were selected to improve the titer of rosin in E. coli. Subsequently, four UGTs from the UGT91R subfamily were identified to catalyze the formation of rosavin from rosin, with SlUGT91R1 from Solanum lycopersicum showing the highest activity level. Secondly, production of rosavin was achieved for the first time in E. coli by incorporating the SlUGT91R1 and UDP-arabinose pathway, including UDP-glucose dehydrogenase, UDP-xylose synthase, and UDP-xylose 4-epimerase, into the rosin-producing stain, and the titer reached 430.5 ± 91.4 mg/L. Thirdly, a two-step pathway derived from L-arabinose, composed of L-arabinokinase and UDP-sugar pyrophosphorylase, was developed in E. coli to further optimize the supply of the precursor UDP-arabinose. Furthermore, 1203.7 ± 32.1 mg/L of rosavin was produced from D-glucose and L-arabinose using shake-flask fermentation. Finally, the production of rosavin reached 7539.1 ± 228.7 mg/L by fed-batch fermentation in a 5-L bioreactor. Thus, the microbe-based production of rosavin shows great potential for commercialization. This work provides an effective strategy for the biosynthesis of other valuable natural products with arabinose-containing units from D-glucose and L-arabinose., Competing Interests: Declaration of competing interest The authors declare the following competing financial interest(s): This work has been included in a patent application by the Tianjin Institute of Industrial Biotechnology., (Copyright © 2024 International Metabolic Engineering Society. Published by Elsevier Inc. All rights reserved.)

Rhodiola L. is a genus that has undergone rapid radiation in the mid-Miocene and may represent a typic case of adaptive radiation. Many species of Rhodiola have also been widely used as an important adaptogen in traditional medicines for centuries. However, a lack of high-quality chromosome-level genomes hinders in-depth study of its evolution and biosynthetic pathway of secondary metabolites. Here, we assembled two chromosome-level genomes for two Rhodiola species with different chromosome number and sexual system. The assembled genome size of R. chrysanthemifolia (2n = 14; hermaphrodite) and R. kirilowii (2n = 22; dioecious) were of 402.67 and 653.62 Mb, respectively, with approximately 57.60% and 69.22% of transposable elements (TEs). The size difference between the two genomes was mostly due to proliferation of long terminal repeat-retrotransposons (LTR-RTs) in the R. kirilowii genome. Comparative genomic analysis revealed possible gene families responsible for high-altitude adaptation of Rhodiola, including a homolog of plant cysteine oxidase 2 gene of Arabidopsis thaliana (AtPCO2), which is part of the core molecular reaction to hypoxia and contributes to the stability of Group VII ethylene response factors (ERF-VII). We found extensive chromosome fusion/fission events and structural variations between the two genomes, which might have facilitated the initial rapid radiation of Rhodiola. We also identified candidate genes in the biosynthetic pathway of salidroside. Overall, our results provide important insights into genome evolution in plant rapid radiations, and possible roles of chromosome fusion/fission and structure variation played in rapid speciation., (© 2023 Society for Experimental Biology and John Wiley & Sons Ltd.)

Rhodiola rosea L. is a vulnerable species in the Altai Republic (AR) and Russia in general. For the first time on the territory of AR, studies of the adaptive capabilities of the species and genetic differentiation using ISSR markers were carried out in seven cenopopulations (CP) of R. rosea in 2018 and 2020. The research was founded on the notion of conducting a comparative analysis of the morphogenetic structure of Rhodiola rosea populations in various ecological and geographical conditions of AR. The aim of this work is to evaluate the variability of morphometric traits of sexually mature living female R. rosea plants and to conduct a comparative analysis of genetic variability in cenopopulations (CP) both under undisturbed conditions and under stressful conditions of anthropogenic impact (grazing). Of the 8 primers used, HB12 turned out to be the most informative. The percentage of polymorphic loci in the populations between 0 and 88%. Two populations, located in favorable conditions at relatively low absolute altitudes (2000 m above sea level) (masl) in the undisturbed habitats of the Katun and Altai reserves of AR, were characterized by higher polymorphism. The share of polymorphic loci reached 80%. According to the analysis of statistical data, the highest values of morphometric parameters of the aerial parts of R. rosea plants and the highest potential seed productivity were also recorded in these habitats. Representatives of two high-mountain CPs (2400-2500 masl) in the Sailyugemsky National Park (SNP) were characterized by the lowest genetic polymorphism. Their genetic structure is the most homogeneous, since we have not found polymorphic loci. Due to spatial isolation, these individuals are reliably genetically differentiated. In addition, individuals of one type were subjected to stressful anthropogenic impact (grazing). Therefore, the smallest sizes and lowest potential seed productivity were recorded. Our research shows that alpine populations of R. rosea in AR, under conditions of anthropogenic stress, need protection for their gene pool.

The Tibetan Plateau and adjacent mountain regions (TP; including the Tibetan Plateau, Himalaya, Hengduan Mountains and Mountains of Central Asia) harbor great biodiversity, some lineages on which may have undergone rapid radiations. However, only a few studies have investigated the evolutionary pattern of such diversification in depth using genomic data. In this study, we reconstructed a robust phylogeny backbone of Rhodiola, a lineage that may have undergone rapid radiation in the TP, using Genotyping-by-sequencing data, and conducted a series of gene flow and diversification analyses. The concatenation and coalescent-based methods yield similar tree topologies, and five well-supported clades were revealed. Potential gene flow and introgression events were detected, both between species from different major clades and closely related species, suggesting pervasive hybridization and introgression. An initial rapid and later slowdown of the diversification rate was revealed, indicating niche filling. Molecular dating and correlation analyses showed that the uplift of TP and global cooling in the mid-Miocene might have played an important role in promoting the rapid radiation of Rhodiola. Our work demonstrates that gene flow and introgression might be an important contributor to rapid radiation possibly by quickly reassembling old genetic variation into new combinations., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Representatives of the Crassulaceae family's genus Rhodiola are succulents, making them distinctive in a changing environment. One of the most significant tools for analyzing plant resources, including numerous genetic processes in wild populations, is the analysis of molecular genetic polymorphism. This work aimed to look at the polymorphisms of allelic variations of the superoxide dismutase (SOD) and auxin response factor (ARF) gene families, as well as the genetic diversity of five Rhodiola species, using the retrotransposons-based fingerprinting approach. The multi-locus exon-primed intron-crossing (EPIC-PCR) profiling approach was used to examine allelic variations in the SOD and ARF gene families. We implemented the inter-primer binding site (iPBS) PCR amplification technique for genome profiling, which demonstrated a significant level of polymorphism in the Rhodiola samples studied. Natural populations of Rhodiola species have a great capacity for adaptation to unfavorable environmental influences. The genetic variety of wild populations of Rhodiola species leads to their improved tolerance of opposing environmental circumstances and species evolutionary divergence based on the diversity of reproductive systems.

Background: Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates., Results: The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP., Conclusion: The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species., (© 2023. The Author(s).)

Aims: Rhodiola was found to be a potential treatment for nonalcoholic fatty liver disease (NAFLD). The macrophage migration inhibitory factor (MIF)-regulated lipophagy and lipid metabolism might be the therapeutic targets of Rhodiola., Main Methods: A 16-week high-fat diet (HFD) was used to simulate a NAFLD mouse model. Rhodiola extract or normal saline were administrated to mice. Blood was collected to assess blood glucose and insulin, and livers were harvested to assess lipid accumulation and metabolism. In cell experiments, the active ingredient of Rhodiola, salidroside, and recombinant MIF protein (rMIF) were used to treat palmitate (PA)-incubated HepG2 cells, with MIF-siRNA or NC-siRNA transfection. Then, the level of lipophagy and lipid metabolism was examined., Key Findings: Rhodiola improved lipid accumulation and metabolism disorder of HFD mice. The oil red O staining of the liver showed that increased lipid droplets in the NAFLD liver could be relieved by Rhodiola; Rhodiola also alleviated the increasing body weight, liver weight, and HOMA-IR index of HFD mice. Results in cell experiments were consistent: salidroside relieved the lipid droplet accumulation and triglyceride release in PA cells, as well as reduced lipophagosome and lipid metabolism disorder in PA cells. However, all these effects of salidroside were partially blocked by MIF-siRNA transfection., Significance: Rhodiola reduces lipid accumulation in the liver of NAFLD by facilitating the MIF pathway and the downstream lipophagy and lipid metabolism. MIF may be an endogenous regulator of liver lipophagy and lipid metabolism and a potential therapeutic target for NAFLD., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Rhodiola imbricata (Crassulaceae) is a traditional trans-Himalayan endangered medicinal herb with immense therapeutic applications. Over the years, over-exploitation, un-managed harvesting, and lack of captive cultivation procedures persuaded threat to its wild habitat. Plant tissue culture and RNA-Seq-based molecular bioprospection of key regulatory genes aid the understanding of molecular dynamics involved in specialized metabolites (phenylethanoids and phenylpropanoids) biosynthesis and its sustainable production. Hence, comparative transcriptomic analysis was performed using leaf and root tissues from the wild and tissue-cultured plants, revealing tissue-specific production of salidroside and rosavin. The transcriptome profiling resulted in 345 million high-quality reads yielding 92,380 unique transcripts with an N 50 of 1260 bp. Tissue-specific gene expression analysis revealed that both phenylethanoids and phenylpropanoids biosynthesis are predominantly associated with the shikimate pathway. In addition to RNA-Seq data, the downstream biosynthesis pathways genes viz., phospho-2-dehydro-3-deoxyheptonate aldolase (DAHPS), 3-dehydroquinate synthase (DHQS), shikimate kinase (SK), chorismate mutase (CM), arogenate dehydrogenase (TYRAAT), aromatic-L-amino-acid decarboxylase (TDC), phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4-CL), cinnamoyl-CoA reductase (CCR), and cinnamyl alcohol dehydrogenase (CAD) showed higher expression pattern in wild plant tissues compared to tissue-cultured plants. The transcript fold expression determined by RT-qPCR results followed similar patterns as those observed in RNA-seq and targeted metabolite profiling data. Salidroside and rosavin content in wild plants exhibited 2.40 fold and 1.77 fold increase accumulation compared to the tissue-cultured plant. The present investigation explained the tissue and condition-specific significant differences between the expression of proposed biosynthetic pathway genes and salidroside and rosavin content. Additionally, NAC, bHLH, and ARF were the most abundant transcription factor families found in the transcriptomic analysis of R. imbricata. The generated transcriptome dataset provides a valuable gene(s)/transcription factors hub that can be used for the sustainable production of salidroside and rosavin in R. imbricata under tissue culture conditions., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Background: As a valuable medicinal plant, Rhodiola has a very long history of folk medicine used as an important adaptogen, tonic, and hemostatic. However, our knowledge of the chloroplast genome level of Rhodiola is limited. This drawback has limited studies on the identification, evolution, genetic diversity and other relevant studies on Rhodiola., Results: Six Rhodiola complete chloroplast genomes were determined and compared to another Rhodiola cp genome at the genome scale. The results revealed a cp genome with a typical quadripartite and circular structure that ranged in size from 150,771 to 151,891 base pairs. High similarity of genome organization, gene number, gene order, and GC content were found among the chloroplast genomes of Rhodiola. 186 (R. wallichiana) to 200 (R. gelida) SSRs and 144 pairs of repeats were detected in the 6 Rhodiola cp genomes. Thirteen mutational hotspots for genome divergence were determined and could be used as candidate markers for phylogenetic analyses and Rhodiola species identification. The phylogenetic relationships inferred by members of Rhodiola cluster into two clades: dioecious and hermaphrodite. Our findings are helpful for understanding Rhodiola's taxonomic, phylogenetic, and evolutionary relationships., Conclusions: Comparative analysis of chloroplast genomes of Rhodiola facilitates medicinal resource conservation, phylogenetic reconstruction and biogeographical research of Rhodiola., (© 2022. The Author(s).)

Background and Aims: Many plant taxa in the Qinghai-Tibetan Plateau (QTP) and the Hengduan Mountains (HM) radiated rapidly during the Quaternary but with frequent secondary contact between diverging populations. Incomplete lineage sorting and introgressive hybridization might be involved during the rapid radiation, but their effects on phylogeography have not been fully determined., Methods: We investigated the chloroplast DNA (cpDNA)/internal transcribed spacer (ITS) sequence variations of 611 samples of Rhodiola bupleuroides, R. discolor, R. fastigiata and R. chrysanthemifolia from the QTP and HM to compare the phylogeographic patterns between the four species with different evolutionary histories, geographic ranges and reproductive modes., Key Results: The divergence times of these species were consistent with the last peak of in situ speciation in the HM. While closely related species exhibited different phylogeographic patterns, they shared several ribotypes and haplotypes in sympatric populations, suggesting introgressive hybridization. A significant phylogenetic discordance between ribotypes and haplotypes was detected in three species, implying incomplete lineage sorting. Rhodiola discolor houses an extraordinary richness of cpDNA haplotypes, and this finding may be attributed to adaptive radiation., Conclusion: In addition to geographic isolation and climate oscillations during the Quaternary, both introgressive hybridization and incomplete lineage sorting play important roles in species that experienced rapid diversification in the QTP and HM., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Annals of Botany Company. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: The contents of some its crucial metabolites tend to decrease when Rhodiola crenulata is cultured at low altitude. Interestingly, it was found that an endophyte, Phialocephala fortinii, could alleviate this problem., Results: There were 16 151 differential genes including 14 706 up-regulated and 1445 down-regulated unigenes with significant differences (P < 0.05), and a total of 1432 metabolites exhibited statistically significant (P < 0.05) metabolic differences comprising 27 different marker metabolites which showed highly significant values of VIP > 5 and P < 0.01. Results highlight differential regulation of 20 enzymatic genes that are involved in the biosynthesis of five different marker metabolites including acetaldehyde, homocysteine, cyclopropylamine, 1-pyrrolinium and halistanol sulfate., Conclusions: The positive physiological effect of P. fortinii on R. crenulata encompasses differential regulation in carbohydrate metabolism, lipid metabolism and secondary metabolite synthesis. © 2020 Society of Chemical Industry., (© 2020 Society of Chemical Industry.)

The amount of plastome sequence data available has soared in the last decade, but the nature of plastome evolution during rapid radiations is largely unknown. Moreover, although there is increasing evidence showing that plastomes may have undergone adaptive evolution in order to allow adaptation to various environments, few studies have systematically investigated the role of the plastome in alpine adaptation. To address these questions, we sequenced and analyzed 12 representative species of Rhodiola, a genus which includes ca. 70 perennial herbs mainly growing in alpine habitats in the Qinghai-Tibet Plateau and the Hengduan Mountains. Rapid radiation in this genus was triggered by the uplift of the Qinghai-Tibet Plateau. We also included nine species of Crassulaceae as the outgroups. All plastomes were conserved with respect to size, structure, and gene content and order, with few variations: each contained 134 genes, including 85 protein-coding genes, 37 tRNAs, 8 rRNAs, and 4 potential pseudogenes. Four types of repeat sequence were detected. Slight contraction and expansion of the inverted repeats were also revealed. Both the genome-wide alignment and sequence polymorphism analyses showed that the inverted repeats and coding regions were more conserved than the single-copy regions and the non-coding regions. Positive selection analyses identified three genes containing sites of positive selection (rpl16, ndhA, ndhH), and one gene with a faster than average rate of evolution (psaA). The products of these genes may be involved in the adaptation of Rhodiola to alpine environments such as low CO 2 concentration and high-intensity light., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Roseroot (Rhodiola rosea L.) is a medicinal plant with adaptogenic properties and several pharmaceutically important metabolites. In this study, a full length cDNA encoding a UDPG gene of roseroot was identified, cloned and characterized. Its ORF (1425 bp) was transferred into E. coli, where the expression of the recombinant enzyme was confirmed. To monitor the enzyme activity, 3 precursors (tyramine, 4-hydroxyphenylpyruvate & tyrosol) of salidroside biosynthesis pathway were added to roseroot callus cultures and samples were harvested after 1, 6, 12, 24, 48 & 96 h. Along with the controls (without precursor feeding), each sample was subjected to HPLC and qRT-PCR for phytochemical and relative UDP-glycosyltransferase gene expression analysis, respectively. The HPLC analysis showed that the salidroside content significantly increased; reaching 0.5% of the callus dry weight (26-fold higher than the control) after 96 h when 2 mM tyrosol was given to the media. The expression of the UDP-glycosyltransferase increased significantly being the highest at 12 h after the feeding. The effect of tyramine and 4-hydroxyphenylpyruvate was not as pronounced as of tyrosol. Here, we introduce a R. rosea specific UDPG gene and its expression pattern after biotransformation of intermediate precursors in in vitro roseroot callus cultures., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Climate change profoundly influences species distributions. These effects are evident in poleward latitudinal range shifts for many taxa, and upward altitudinal range shifts for alpine species, that resulted from increased annual global temperatures since the Last Glacial Maximum (LGM, ca. 22,000 BP). For the latter, the ultimate consequence of upward shifts may be extinction as species in the highest alpine ecosystems can migrate no further, a phenomenon often characterized as "nowhere to go". To predict responses to climate change of the alpine plants on the Qinghai-Tibetan Plateau (QTP), we used ecological niche modelling (ENM) to estimate the range shifts of 14 Rhodiola species, beginning with the Last Interglacial (ca. 120,000-140,000 BP) through to 2050. Distributions of Rhodiola species appear to be shaped by temperature-related variables. The southeastern QTP, and especially the Hengduan Mountains, were the origin and center of distribution for Rhodiola, and also served as refugia during the LGM. Under future climate scenario in 2050, Rhodiola species might have to migrate upward and northward, but many species would expand their ranges contra the prediction of the "nowhere to go" hypothesis, caused by the appearance of additional potential habitat concomitant with the reduction of permafrost with climate warming.

Quaternary climatic oscillations have had tremendous effects on current distribution of species. Previous studies unraveled multiple microrefugia on the Qinghai-Tibetan Plateau (QTP) in two woody plants. Still we know little whether herbs growing in forests responded to climatic oscillations similarly. We herein conducted a phylogeographic study on Rhodiola sect. Trifida, an herbaceous group endemic to the QTP, which mainly growing on the forest floors, using plastid and ITS sequences as well as ecological niche modeling. The origin and divergence of major clades of sect. Trifida were in accordance with the last phase of the QTP uplifts. Mismatch distribution analysis indicated a range expansion dated to ca. 135 thousand years ago. A high frequency and an even distribution of private haplotypes in both plastid and ITS data sets throughout the distribution of sect. Trifida were detected. The ecological niche modeling results showed that there were suitable habitats on the QTP platform during the LGM. Our results found that multiple microrefugia existed on the QTP platform, supporting the hypothesis that species with similar geographic distribution and inhabiting the same community had similar responses to the Quaternary climatic oscillations. Furthermore, species delimitations in sect. Trifida need to be tested based on integrative evidence from morphological, ecological and genetic data., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Siberian ginseng (Eleutherococcus senticosus, Araliaceae) and roseroot (Rhodiola rosea, Rosaceae) are popular herbal supplements which have been shown to improve resilience to conditions such as stress and exhaustion. Using DNA barcoding methods we tested 25 Siberian ginseng and 14 roseroot products which are widely available to UK customers to test whether the herbal ingredient stated on the label is also in the product. All Siberian ginseng supplements contained E. senticosus, however, 36% also contained an Eleutherococcus species other than E. senticosus. In three out of the 13 roseroot products which produced amplifiable DNA, we could only retrieve sequences matching alfalfa (declared on the product label) and fenugreek (not declared). In the other 10 supplements Rhodiola was detected but only five matched the target species R. rosea. As DNA can get severely degraded during the manufacturing process we did not take the absence of Rhodiola DNA as proof for a compromised product. Contamination could explain the presence of non-target species such as fenugreek but is unlikely to be account for the detection of congeneric Rhodiola species in roseroot preparations. Our results therefore suggest that the substitution or mixing of the target medicinal ingredient in these two popular supplements with other species is common., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Salidroside is a bioactive tyrosine-derived phenolic natural product found in medicinal plants under the Rhodiola genus. In addition to their anti-fatigue and anti-anoxia roles in traditional medicine, Rhodiola total extract and salidroside have also displayed medicinal properties as anti-cardiovascular diseases and anti-cancer agents. The resulting surge in global demand of Rhodiola plants and salidroside has driven some species close to extinction. Here, we report the full elucidation of the Rhodiola salidroside biosynthetic pathway utilizing the first comprehensive transcriptomics and metabolomics datasets for Rhodiola rosea. Unlike the previously proposed pathway involving separate decarboxylation and deamination enzymatic steps from tyrosine to the key intermediate 4-hydroxyphenylacetaldehyde (4-HPAA), Rhodiola contains a pyridoxal phosphate-dependent 4-HPAA synthase that directly converts tyrosine to 4-HPAA. We further identified genes encoding the subsequent 4-HPAA reductase and tyrosol:UDP-glucose 8-O-glucosyltransferase, respectively, to complete salidroside biosynthesis in Rhodiola. We show that heterologous production of salidroside can be achieved in the yeast Saccharomyces cerevisiae as well as the plant Nicotiana benthamiana through transgenic expression of Rhodiola salidroside biosynthetic genes. This study provides new tools for engineering sustainable production of salidroside in heterologous hosts., (Copyright © 2017 The Author. Published by Elsevier Inc. All rights reserved.)

Endophyte is a factor that affects the physiology and metabolism of plant. However, limited information is available on the mechanism of interaction between endophyte and plant. To investigate the effects of endophytic fungus ZPRs-R11, that is, Trimmatostroma sp., on salidroside and tyrosol accumulations in Rhodiola crenulata, signal transduction, enzyme gene expression, and metabolic pathway were investigated. Results showed that hydrogen peroxide (H 2 O 2 ), nitric oxide (NO), and salicylic acid (SA) involved in fungus-induced salidroside and tyrosol accumulations. NO acted as an upstream signal of H 2 O 2 and SA. No up- or down-stream relationship was observed, but mutual coordination existed between H 2 O 2 and SA. Rate-limiting enzyme genes with the maximum expression activities were UDP-glucosyltransferase, tyrosine decarboxylase (TYDC), monoamine oxidase, phenylalanine ammonialyase (PAL), and cinnamic-4-hydroxylase sequentially. Nevertheless, the genes of tyrosine transaminase and pyruvate decarboxylase only indicated slightly higher activities than those in control. Thus, TYDC and PAL branches were the preferential pathways in ZPRs-R11-induced salidroside and tyrosol accumulation. Trimmatostroma sp. was a potential fungus for promoting salidroside and tyrosol accumulations. The present data also provided scientific basis for understanding complex interaction between endophytic fungus and R. crenulata.

The roots and rhizomes of Rhodiola crenulata and R. rosea have been used worldwide as adaptogens for hundreds of years. However, rapid growth in demand has resulted in merchants using other species of Rhodiola as adulterants. Here, we surveyed 518 individuals representing 47 of the 55 species in the genus, including 253 R. crenulata individuals from 16 populations and 98 R. rosea individuals from 11 populations, to evaluate the utility of the internal transcribed spacer 2 (ITS2) barcode for identification of Rhodiola species. We detected six haplotypes in R. crenulata and only one haplotype in R. rosea. An obvious overlap between intra- and inter-specific distance was detected, and the authentication efficacy of ITS2, which was assessed by BLAST1, a nearest distance method, and a tree test, was much lower than in other groups. However, R. crenulata and R. rosea could be exactly identified. Analysis showed that the secondary structure of ITS2 differs in R. crenulata and its closest relatives. Our results demonstrated that both a mini barcode from ITS2 and the structure of ITS2 are effective markers for the identification of R. crenulata and R. rosea. This study represents the most comprehensive database of ITS2 barcodes in Rhodiola to date and will be useful in Rhodiola species identification.

Three rosmarinic acid analogs produced by recombinant Escherichia coli, two xanthones from fungi and honokiol from plants, were explored as the substrates of E. coli harboring a glucosyltransferase mutant UGT73B6 FS to generate phenolic glucosides. Six new and two known compounds were isolated from the fermentation broth of the recombinant strain of the feeding experiments, and the compounds were identified by spectroscopy. The biotransformation of rosmarinic acid analogs and xanthones into corresponding glucosides was presented for the first time. This study not only demonstrated the substrate flexibility of the glucosyltransferase mutant UGT73B6 FS toward aromatic alcohols but also provided an effective and economical method to produce phenolic glucosides by fermentation circumventing the use of expensive precursor UDP-glucose.

Rhodiola crenulata, a well-known medicinal Tibetan herb, is mainly grown in high-altitude regions of the Tibet, Yunnan, and Sichuan provinces in China. In the past few years, increasing numbers of studies have been published on the potential pharmacological activities of R. crenulata, strengthening our understanding into its putitive active ingredient composition, pharmacological activity, and mechanism of action. These findings also provide strong evidence supporting the important medicinal and economical value of R. crenulata. Consequently, some Rhodiola species are becoming endangered because of overexploitation and environmental destruction. However, little is known about the genetic and genomic information of any Rhodiola species. Here we report the first draft assembly ofthe R. crenulata genome, which was 344.5 Mb (25.7 Mb Ns), accounting for 82% of the estimated genome size, with a scaffold N50 length of 144.7 kb and a contig N50 length of 25.4 kb. The R. crenulata genome is not only highly heterozygous but also highly repetitive, with ratios of 1.12% and 66.15%, respectively, based on the k-mer analysis. Furthermore, 226.6 Mb of transposable elements were detected, of which 77.03% were long terminal repeats. In total, 31 517 protein-coding genes were identified, capturing 86.72% of expected plant genes in BUSCO. Additionally, 79.73% of protein-coding genes were functionally annotated. R. crenulata is an important medicinal plant and also a potentially interesting model species for studying the adaptability of Rhodiola species to extreme environments. The genomic sequences of R. crenulata will be useful for understanding the evolutionary mechanism of the stress resistance gene and the biosynthesis pathways of the different medicinal ingredients, for example, salidroside in R. crenulata., (© The Authors 2017. Published by Oxford University Press.)

Polyphenols, which include phenolic acids, flavonoids, stilbenes, and phenylethanoids, are generally known as useful antioxidants. Tyrosol, hydroxytyrosol, and salidroside are typical phenylethanoids. Phenylethanoids are found in plants such as olive, green tea, and Rhodiola and have various biological activities, including the prevention of cardiovascular diseases, cancer, and brain damage. We used Escherichia coli to synthesize three phenylethanoids, tyrosol, hydroxytyrosol, and salidroside. To synthesize tyrosol, the aromatic aldehyde synthase (AAS) was expressed in E. coli. Hydroxytyrosol was synthesized using E. coli harboring AAS and HpaBC, which encodes hydroxylase. In order to synthesize salidroside, 12 uridine diphosphate-dependent glycosyltransferases (UGTs) were screened and UGT85A1 was found to convert tyrosol to salidroside. Using E. coli harboring AAS and UGT85A1, salidroside was synthesized. Through the optimization of these three E. coli strains, we were able to synthesize 531 mg/L tyrosol, 208 mg/L hydroxytyrosol, and 288 mg/L salidroside, respectively.

Rhodiola plants are a valuable resource in traditional Chinese medicine. The objective of this study was to evaluate the correlation between ribosomal DNA internal transcribed spacer (ITS) sequences and the three active components in Rhodiola plants. For this, we determined ITS sequence polymorphisms and the concentrations of active components salidroside, tyrosol, and gallic acid in different Rhodiola species from the Tibetan Plateau. In a total of 23 Rhodiola samples, 16 different haplotypes were defined based on their ITS sequences. Analysis of the active components in these same samples revealed that salidroside was not detected in species with haplotypes H 4 , H 5 , or H 10 , tyrosol was not detected with haplotypes H 3 , H 5 , H 7 , H 10 , H 14 , or H 15 , and gallic acid was detected in with all haplotypes except H 14 and H 15 . In addition, the concentrations of salidroside, tyrosol and gallic acid varied between samples with different haplotypes as well as those with the same haplotype, implying that no significant correlation exists between haplotype and salidroside, tyrosol or gallic acid concentrations. However, a statistically significant positive correlation was observed for among these three active components.

One strain of endophytic fungus ZPRa-R-1 was obtained for the capacity of promoting production of salidroside in Rhodiola crenulata. To explain the mechanism of salidroside biosynthesis in host plant, eight housekeeping genes were evaluated, and the evaluation method was created for the expression activities of four key enzyme genes PAL (phenylalanine ammonia-lyase), TyDC (tyrosine decarboxylase), TAT (tyrosine transaminase), UDPGT (UDP-glucosyltransferase) referenced double reference genes in biosynthesis pathway of salidroside in R. crenulata. Stabilities of housekeeping genes were confirmed by real-time fluorescent quantitative PCR technology and three softwares including geNorm, NormFinder and BestKeeper, then relative expressions of key enzyme genes were analysized by the 2-ΔΔCt method. The results showed that the most stable gene was GAPDH, followed by PCS, and the most appropriate reference of internal genes were combination with two genes in R. crenulata inoculated with endophytic fungus ZPRa-R-1. Under symbiosis conditions, regularity changes of key enzyme genes affected by endophytic fungus ZPRa-R-1 were as follows: the relative expression activity of PAL attached to peak value, which was 4.9 times as that of control group when inoculated ten days. The relative expression of TyDC reached the maximum value, which was 2.8 times of that control after inoculating 12 days. The relative expression of UDPGT actually reach 17.1 times than that of control after inoculating 8 days. However, the relative expression of TAT was not affected by this fungus. The changes of four key enzyme genes are positively correlated with the changes of salidroside content in R. crenulata.

Gastrodin, a phenolic glycoside, is the key ingredient of Gastrodia elata, a notable herbal plant that has been used to treat various conditions in oriental countries for centuries. Gastrodin is extensively used clinically for its sedative, hypnotic, anticonvulsive and neuroprotective properties in China. Gastrodin is usually produced by plant extraction or chemical synthesis, which has many disadvantages. Herein, we report unprecedented microbial synthesis of gastrodin via an artificial pathway. A Nocardia carboxylic acid reductase, endogenous alcohol dehydrogenases and a Rhodiola glycosyltransferase UGT73B6 transformed 4-hydroxybenzoic acid, an intermediate of ubiquinone biosynthesis, into gastrodin in Escherichia coli. Pathway genes were overexpressed to enhance metabolic flux toward precursor 4-hydroxybenzyl alcohol. Furthermore, the catalytic properties of the UGT73B6 toward phenolic alcohols were improved through directed evolution. The finally engineered strain produced 545mgl(-1) gastrodin in 48h. This work creates a new route to produce gastrodin, instead of plant extractions and chemical synthesis., (Copyright © 2016 International Metabolic Engineering Society. Published by Elsevier Inc. All rights reserved.)

Salidroside, a plant secondary metabolite in Rhodiola, has been demonstrated to have several adaptogenic properties as a medicinal herb. Due to the limitation of plant source, microbial production of salidroside by expression of plant uridine diphosphate glycosyltransferase (UGT) is promising. However, glycoside production usually remains hampered by poor expression of plant UGTs in microorganisms. Herein, we achieved salidroside production by expression of Rhodiola UGT72B14 in Escherichia coli (E. coli) and codon optimization was accordingly applied. UGT72B14 expression was optimized by changing 278 nucleotides and decreasing the G+C content to 51.05% without altering the amino acid sequence. The effect of codon optimization on UGT72B14 catalysis for salidroside production was assessed both in vitro and in vivo. In vitro, salidroside production by codon-optimized UGT72B14 is enhanced because of a significantly improved protein yield (increased by 4.8-fold) and an equivalently high activity as demonstrated by similar kinetic parameters (K M and V max), compared to that by wild-type protein. In vivo, both batch and fed-batch cultivation using the codon-optimized gene resulted in a significant increase in salidroside production, which was up to 6.7 mg/L increasing 3.2-fold over the wild-type UGT72B14.

Rhodiola alsia, which has been used widely in traditional Chinese medicine for a considerable time, grows on moist habitats at high altitude near the snow line. Microsatellite loci were developed for R. alsia to investigate its population genetics. In total, 17 polymorphic microsatellites were developed based on ESTs from the Illumina HiSeq(TM) 2000 platform. The microsatellite loci were checked for variability using 80 individuals of R. alsia sampled from four locations on the Qinghai-Tibet Plateau. The total number of alleles per locus ranged from 10 to 20, and the observed heterozygosity ranged from 0.000 to 1.000. The null allele frequency ranged from 0.000 to 0.324. These microsatellites are expected to be helpful in future studies of population genetics in R. alsia and related species.

DNA barcoding, the identification of species using one or a few short standardized DNA sequences, is an important complement to traditional taxonomy. However, there are particular challenges for barcoding plants, especially for species with complex evolutionary histories. We herein evaluated the utility of five candidate sequences - rbcL, matK, trnH-psbA, trnL-F and the internal transcribed spacer (ITS) - for barcoding Rhodiola species, a group of high-altitude plants frequently used as adaptogens, hemostatics and tonics in traditional Tibetan medicine. Rhodiola was suggested to have diversified rapidly recently. The genus is thus a good model for testing DNA barcoding strategies for recently diversified medicinal plants. This study analyzed 189 accessions, representing 47 of the 55 recognized Rhodiola species in the Flora of China treatment. Based on intraspecific and interspecific divergence and degree of monophyly statistics, ITS was the best single-locus barcode, resolving 66% of the Rhodiola species. The core combination rbcL+matK resolved only 40.4% of them. Unsurprisingly, the combined use of all five loci provided the highest discrimination power, resolving 80.9% of the species. However, this is weaker than the discrimination power generally reported in barcoding studies of other plant taxa. The observed complications may be due to the recent diversification, incomplete lineage sorting and reticulate evolution of the genus. These processes are common features of numerous plant groups in the high-altitude regions of the Qinghai-Tibetan Plateau.

The adulteration of herbal products is a threat to consumer safety. Here we surveyed the species composition of commercial Rhodiola products using DNA barcoding as a supervisory method. A Rhodiola dietary supplement DNA barcode database was successfully constructed using 82 voucher samples from 10 Rhodiola species. Based on the DNA barcoding standard operating procedure (SOP), we used this database to identify 100 Rhodiolae Crenulatae Radix et Rhizoma decoction piece samples that were purchased from drug stores and hospitals. The results showed that only 36 decoction piece sequences (40%) were authentic R. crenulata, which is recorded in Chinese Pharmacopeia, whereas the other samples were all adulterants and may indicate a potential safety issue. Among the adulterants, 35 sequences (38.9%) were authenticated as R. serrata, nine sequences (10%) were authenticated as R. rosea, which is documented in the United States Pharmacopeia, and the remaining samples were authenticated as other three Rhodiola species. This result indicates decoction pieces that are available in the market have complex origins and DNA barcoding is a convenient tool for market supervision.

Transcriptome sequencing is a powerful tool for the assessment of gene expression and the identification and characterization of molecular markers in non-model organisms. Rhodiola algida L. (Crassulaceae), endemic to the Qinghai-Tibetan Plateau, has long been used in traditional Chinese medicine to prevent altitude sickness and eliminate fatigue. Illumina-based high-throughput transcriptome sequencing of aboveground and underground tissues of R. algida respectively yielded 5.40 million and 5.18 million clean reads. A total of 82,664 unigenes averaging 577 bp in length were generated from the aboveground clean reads, with 86,237 unigenes of 502-bp mean length obtained from the underground tissues. Of 55,028 unigenes compared with sequences in UniProt databases, 20,413 unigenes had significant similarities with existing sequences in NR, NT, Swiss-Prot, GO, KEGG, and COG databases. Single nucleotide polymorphism (SNP) analysis identified 237,294 SNPs from 154,636 contigs of aboveground tissues and 197,540 SNPs from 144,963 underground-derived contigs. The information uncovered in this study should serve as a valuable resource for the characterization of important traits related to secondary metabolite formation and for the identification of associated molecular mechanisms., (Copyright © 2014 Elsevier B.V. All rights reserved.)

The evolution and current distribution of the Sino-Tibetan flora have been greatly affected by historical geological events, such as the uplift of the Qinghai-Tibetan Plateau (QTP), and Quaternary climatic oscillations. Rhodiola kirilowii, a perennial herb with its distribution ranging from the southeastern QTP and the Hengduan Mountains (HM) to adjacent northern China and central Asia, provides an excellent model to examine and disentangle the effect of both geological orogeny and climatic oscillation on the evolutionary history of species with such distribution patterns. We here conducted a phylogeographic study using sequences of two chloroplast fragments (trnL-F and trnS-G) and internal transcribed spacers in 29 populations of R. kirilowii. A total of 25 plastid haplotypes and 12 ITS ribotypes were found. Molecular clock estimation revealed deep divergence between the central Asian populations and other populations from the HM and northern China; this split occurred ca. 2.84 million year ago. The majority of populations from the mountains of northern China were dominated by a single haplotype or ribotype, while populations of the HM harbored both high genetic diversity and high haplotype diversity. This distribution pattern indicates that HM was either a diversification center or a refugium for R. kirilowii during the Quaternary climatic oscillations. The present distribution of this species on mountains in northern China may have resulted from a rapid glacial population expansion from the HM. This expansion was confirmed by the mismatch distribution analysis and negative Tajima's D and Fu's FS values, and was dated to ca. 168 thousand years ago. High genetic diversity and population differentiation in both plastid and ITS sequences were revealed; these imply restricted gene flow between populations. A distinct isolation-by-distance pattern was suggested by the Mantel test. Our results show that in old lineages, populations may harbour divergent genetic forms that are sufficient to maintain or even increase overall genetic diversity despite fragmentation and low within-population variation.

Rhodiola sachalinensis is an endangered species with important medicinal value. We used inter-simple sequence repeat (ISSR) and methylation-sensitive amplified polymorphism (MSAP) markers to analyze genetic and epigenetic differentiation in different populations of R. sachalinensis, including three natural populations and an ex situ population. Chromatographic fingerprint was used to reveal HPLC fingerprint differentiation. According to our results, the ex situ population of R. sachalinensis has higher level genetic diversity and greater HPLC fingerprint variation than natural populations, but shows lower epigenetic diversity. Most genetic variation (54.88%) was found to be distributed within populations, and epigenetic variation was primarily distributed among populations (63.87%). UPGMA cluster analysis of ISSR and MSAP data showed identical results, with individuals from each given population grouping together. The results of UPGMA cluster analysis of HPLC fingerprint patterns was significantly different from results obtained from ISSR and MSAP data. Correlation analysis revealed close relationships among altitude, genetic structure, epigenetic structure, and HPLC fingerprint patterns (R2 = 0.98 for genetic and epigenetic distance; R2 = 0.90 for DNA methylation level and altitude; R2 = -0.95 for HPLC fingerprint and altitude). Taken together, our results indicate that ex situ population of R. sachalinensis show significantly different genetic and epigenetic population structures and HPLC fingerprint patterns. Along with other potential explanations, these findings suggest that the ex situ environmental factors caused by different altitude play an important role in keeping hereditary characteristic of R. sachalinensis.

Rhodiola L. (Crassulaceae) is a mid-sized plant genus consisting of about 70 species, with most species distributed on the Qinghai-Tibetan Plateau (QTP) and the adjacent areas, and several species in north-east Asia, Europe, and North America. This study explored the origin and diversification history of Rhodiola and tested the biogeographic relationships between the QTP and other regions of the Northern Hemisphere. We sequenced the nuclear ribosomal internal transcribed spacers and eight plastid DNA fragments representing 55 species of Rhodiola, and reconstructed phylogenetic relationships with maximum parsimony, maximum likelihood and Bayesian inference. Several instances of incongruence between the nuclear and the plastid data sets were revealed, which can best be explained by reticulate evolution. Species of Rhodiola and Pseudosedum form a well-supported clade sister to Phedimus. Dating analysis suggested that the origin and diversification times of this group are largely correlated with the extensive uplifts of the Qinghai-Tibetan Plateau. Ancestral state reconstruction supports the hypothesis that Rhodiola originated on the QTP, and then dispersed to other regions of the Northern Hemisphere. Our findings highlight the importance of the uplifts of the Qinghai-Tibetan Plateau in promoting species diversification and the possible role of reticulate evolution in the diversification process. Our results also suggest the biogeographic significance of QTP as the source area in alpine plant evolution in the Northern Hemisphere., (Copyright © 2014 Elsevier Inc. All rights reserved.)

The phylogeographical patterns of Rhodiola dumulosa, an alpine plant species restrictedly growing in the crevices of rock piles, were investigated based on 4 fragments of the chloroplast genome. To cover the full distribution of R. dumulosa in China, 19 populations from 3 major disjunct distribution areas (northern, central, and northwestern China) were sampled. A total of 5881bp (after alignment) of chloroplast DNA (cpDNA) from 100 individuals were sequenced. The combined cpDNA data set yielded 36 haplotypes. The total genetic diversity of R. dumulosa was remarkably high (H(T) = 0.981). The interpopulation genetic differentiation was significantly large (F(ST) = 0.8537, P < 0.001), possibly due to the long-term isolation of the natural populations. N(ST) was significantly larger than G(ST) (P < 0.001), indicating the presence of phylogeographical structure among the R. dumulosa populations. We propose 2 migration steps to explain the current distribution of R. dumulosa in China. First, this species migrated from refugia in the Qinghai-Tibetan Plateau to northern areas via the intervening highlands when temperatures increased; second, the highland populations migrated toward the mountaintops when temperatures increased further because R. dumulosa is adapted to cold environments. During the second migration step, the common ancestral haplotypes may have been gradually lost.

Despite the strength of climatic variability at high latitudes and upper elevations, we still do not fully understand how plants in North America that are distributed between Arctic and alpine areas responded to the environmental changes of the Quaternary. To address this question, we set out to resolve the evolutionary history of the King's Crown, Rhodiola integrifolia using multi-locus population genetic and phylogenetic analyses in combination with ecological niche modeling. Our population genetic analyses of multiple anonymous nuclear loci revealed two major clades within R. integrifolia that diverged from each other ~ 700 kya: one occurring in Beringia to the north (including members of subspecies leedyi and part of subspecies integrifolia), and the other restricted to the Southern Rocky Mountain refugium in the south (including individuals of subspecies neomexicana and part of subspecies integrifolia). Ecological niche models corroborate our hypothesized locations of refugial areas inferred from our phylogeographic analyses and revealed some environmental differences between the regions inhabited by its two subclades. Our study underscores the role of geographic isolation in promoting genetic divergence and the evolution of endemic subspecies in R. integrifolia. Furthermore, our phylogenetic analyses of the plastid spacer region trnL-F demonstrate that among the native North American species, R. integrifolia and R. rhodantha are more closely related to one another than either is to R. rosea. An understanding of these historic processes lies at the heart of making informed management decisions regarding this and other Arctic-alpine species of concern in this increasingly threatened biome.

Tyrosine decarboxylase initializes salidroside biosynthesis. Metabolic characterization of tyrosine decarboxylase gene from Rhodiola crenulata (RcTYDC) revealed that it played an important role in salidroside biosynthesis. Recombinant 53 kDa RcTYDC converted tyrosine into tyramine. RcTYDC gene expression was induced coordinately with the expression of RcUDPGT (the last gene involved in salidroside biosynthesis) in SA/MeJA treatment; the expression of RcTYDC and RcUDPGT was dramatically upregulated by SA, respectively 49 folds and 36 folds compared with control. MeJA also significantly increased the expression of RcTYDC and RcUDPGT in hairy root cultures. The tissue profile of RcTYDC and RcUDPGT was highly similar: highest expression levels found in stems, higher expression levels in leaves than in flowers and roots. The gene expressing levels were consistent with the salidroside accumulation levels. This strongly suggested that RcTYDC played an important role in salidroside biosynthesis in R. crenulata. Finally, RcTYDC was used to engineering salidroside biosynthetic pathway in R. crenulata hairy roots via metabolic engineering strategy of overexpression. All the transgenic lines showed much higher expression levels of RcTYDC than non-transgenic one. The transgenic lines produced tyramine, tyrosol and salidroside at higher levels, which were respectively 3.21-6.84, 1.50-2.19 and 1.27-3.47 folds compared with the corresponding compound in non-transgenic lines. In conclusion, RcTYDC overexpression promoted tyramine biosynthesis that facilitated more metabolic flux flowing toward the downstream pathway and as a result, the intermediate tyrosol was accumulated more that led to the increased production of the end-product salidroside.

Salidroside, the 8-O-β-D-glucoside of tyrosol, is the main bioactive component of Rhodiola species and is found mainly in the plant roots. It is well known that glucosylation of tyrosol is the final step in the biosynthesis of salidroside; however, the biosynthetic pathway of tyrosol and its regulation are less well understood. A summary of the results of related studies revealed that the precursor of tyrosol might be tyramine, which is synthesized from tyrosine. In this study, a cDNA clone encoding tyrosine decarboxylase (TyrDC) was isolated from Rhodiola sachalinensis A. Bor using rapid amplification of cDNA ends. The resulting cDNA was designated RsTyrDC. RNA gel-blot analysis revealed that the predominant sites of expression in plants are the roots and high levels of transcripts are also found in callus tissue culture. Functional analysis revealed that tyrosine was best substrate of recombinant RsTyrDC. The over-expression of the sense-RsTyrDC resulted in a marked increase of tyrosol and salidroside content, but the levels of tyrosol and salidroside were 274 and 412%, respectively, lower in the antisense-RsTyrDC transformed lines than those in the controls. The data presented here provide in vitro and in vivo evidence that the RsTyrDC can regulate the tyrosol and salidroside biosynthesis, and the RsTyrDC is most likely to have an important function in the initial reaction of the salidroside biosynthesis pathway in R. sachalinensis.

Salidroside, the 8-O-β-D-glucoside of tyrosol, is a novel adaptogenic drug extracted from the medicinal plant Rhodiola sachalinensis A. Bor. Due to the scarcity of R. sachalinensis and its low yield of salidroside, there is great interest in enhancing production of salidroside by biotechnological manipulations. In this study, two putative UDP-glycosyltransferase (UGT) cDNAs, UGT72B14 and UGT74R1, were isolated from roots and cultured cells of methyl jasmonate (MeJA)-treated R. sachalinensis, respectively. The level of sequence identity between their deduced amino acid sequences was ca. 20%. RNA gel-blot analysis established that UGT72B14 transcripts were more abundant in roots, and UGT74R1 was highly expressed in the calli, but not in roots. Functional analysis indicated that recombinant UGT72B14 had the highest level of activity for salidroside production, and that the catalytic efficiency (Vmax/Km) of UGT72B14 was 620% higher than that of UGT74R1. The salidroside contents of the UGT72B14 and UGT74R1 transgenic hairy root lines of R. sachalinensis were also ∼420% and ∼50% higher than the controls, respectively. UGT72B14 transcripts were mainly detected in roots, and UGT72B14 had the highest level of activity for salidroside production in vitro and in vivo., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

Objective: To construct root-specific promoter of plant expression vector and transformate into Rhodiola hairy root., Methods: The expression vector pCA-Tob7: UGTR, driven by tobacco root-specific promoter TobRB7, was constructed from pCAM-BIA1301 by substituting the CaMV 35S promoter and GUS gene with TobRB7 and UGTR(a glycosyltransferase gene), respectively. The pCA-Tob7: UGTR vector was introduced into Agrobacterium tumefaciens strain, and the hairy root of Rhodiola were then transformed., Results: Some Kanamycin resistant plants were positive,which indicated that the expression vector had integrated into Rhodiola hairy root genome., Conclusion: Root-specific promoter of plant expression vector is constructed and the hairy root of Rhodiola are transformed for future research.

Aims: Rhodiola dumulosa (Crassulaceae) is a perennial diploid species found in high-montane areas. It is distributed in fragmented populations across northern, central and northwestern China. In this study, we aimed to (i) measure the genetic diversity of this species and that of its populations; (ii) describe the genetic structure of these populations across the entire distribution range in China; and (iii) evaluate the extent of gene flow among the naturally fragmented populations., Methods: Samples from 1089 individuals within 35 populations of R. dumulosa were collected, covering as much of the entire distribution range of this species within China as possible. Population genetic diversity and structure were analyzed using AFLP molecular markers. Gene flow among populations was estimated according to the level of population differentiation., Important Findings: The total genetic diversity of R. dumulosa was high but decreased with increasing altitude. Population-structure analysis indicated that the most closely related populations were geographically restricted and occurred in close proximity to each other. A significant isolation-by-distance pattern, caused by the naturally fragmented population distribution, was observed. At least two distinct gene pools were found in the 35 sampled populations, one composed of populations in northern China and the other composed of populations in central and northwestern China. The calculation of Nei's gene diversity index revealed that the genetic diversity in the northern China pool (0.1972) was lower than that in the central and northwestern China pool (0.2216). The populations were significantly isolated, and gene flow was restricted throughout the entire distribution. However, gene flow among populations on the same mountain appears to be unrestricted, as indicated by the weak genetic isolation among these populations.

Roseroot, Rhodiola rosea, is a perennial herbaceous plant of the family Crassulaceae. The rhizomes of 95 roseroot clones in the Norwegian germplasm collection were analysed and quantified for their content of the bioactive compounds rosavin, salidroside, rosin, cinnamyl alcohol and tyrosol using HPLC analysis. All five bioactive compounds were detected in all 95 roseroot clones but in highly variable quantities. The ranges observed for the different compounds were for rosavin 2.90-85.95 mg g(-1), salidroside 0.03-12.85 mg g(-1), rosin 0.08-4.75 mg g(-1), tyrosol 0.04-2.15 mg g(-1) and cinnamyl alcohol 0.02-1.18 mg g(-1). The frequency distribution of the chemical content of each clone did not reflect a certain geographic region of origin or the gender of the plant. Significant correlations were found for the contents of several of these bioactive compounds in individual roseroot clones. A low, but not significant correlation was found between AFLP markers previously used to study the genetic diversity of the roseroot clones and their production of the chemical compounds. The maximum level of rosavin, rosin and salidroside observed were higher than for any roseroot plant previously reported in literature, and the roseroot clones characterized in this study might therefore prove to be of high pharmacological value.

Objective: To acquire homozygous tetraploid germplasm of Rhodiola sachalinensis., Method: PEG-mediated protoplast fusions were conducted using callus of Rh. sachalinensis as materials. Protoplast fusion products were embedded and cultured in low-density, low-melting-point agar and marked according to the protoplast size, and single-celled sister lines were established to acquire genetically homozygous tetraploid germplasm., Result: R(D) and R(M) of newborn daughter cells or protoplasm, metaphase cells or protoplasm were approximately in line with the formula R(D) = 0.793 7R(M). The change range in diameter of the diploid cells without fusion, two protoplasts fusion product were: 16.7 microm < or = R < 21.3 microm, 21.0 microm < or = R' < 26.8 microm respectively. There is an overlap between the two diameter ranges. The protoplast inoculation density of 1 x 10(4) cells x mL(-1) was appropriate when protoplasts were anchored by low-intensity, low-melting-point agar. Under the conditions of this density, plating efficiency was high and single cell origin of the sister lines microclones grew rapidly, and it was easy to mark the single cell microclones, and separate from each other to subculture. The chromosome counts results showed that chromosome numbers of diploid and tetraploid of single cell lines were 26 and 52, respectively. The result from flow cytometry assay showed that there is no presence of chimerism in single-cell regeneration plantlets., Conclusion: The results of this study provide a scientific basis for polyploid breeding of Rh. sachalinensis.

Tyrosine decarboxylase (TyrDC) is an important enzyme in the secondary metabolism of several plant species, and was hypothesized to play a key role in the biosynthesis of salidroside, a pharmacologically valuable compound of roseroot. A 1520bp cDNA was cloned and sequenced, and turned out to contain an ORF of 963bp, which encodes a protein of 320 amino acids. The expression of the gene was studied by real-time PCR from leaves and roots of both high and low salidroside producer genotype of roseroot. The gene expression analysis showed the gene to be expressed in leaves as well as in roots; however, the expression was significantly higher in roots, which coincides with the fact that salidroside accumulates preferentially in the underground parts of the plant. The expression was also higher in the line accumulating high levels of salidroside, compared to the line with lower salidroside content. The difference in the expression intensity suggests a decisive role for this enzyme in the salidroside biosynthesis.

Objective: To research the genetic diversity of different Rhodiola rosea geographical populations in Tianshan Mountain, China;, Method: The genetic diversity of eighteen R. rosea geological populations from six niches was estimated using amplified fragment length polymorphism (AFLP) markers. The data of amplified bands were analyzed by the software POPGENE v1.31 (32-bit) and SPSS., Result: The nine primers employed produced a total of 238 discernable and reproducible amplified fragments. There were 228 polymorphic bands. The percentage of polymorphic bands with in different populations was 95.6%. Genetic diversity analysis showed that average number of alleles per loci was Na = 1.4883, effective number of alleles per loci Ne = 1.3907, Neis gene diversity index H = 0.2170, Shannon's information index I = 0.3108, the percentage of polymorphic loci P = 52.71, genetic differentiation among populations Gst = 0.364; UPGMA cluster analysis based on genetic distance data divided eighteen populations into two clusters: Cluster I composed of twelve populations and Cluster II 6 populations which distributed in attitude upper 3 175 m;, Conclusion: Our researches suggest that the best niche of R. rosea was at attitude between 3 150-3 250 m; this region is important for the conservation of R. rosea germplasm resource.

Salidroside, a novel effective adaptogenic drug extracted from the medicinal plant Rhodiola sachalinensis A. Bor, can be derived from phenylalanine or tyrosine. Due to the scarcity of R. sachalinensis and its low yield of salidroside, there is great interest in enhancing production of salidroside by the plant. In this study, a cDNA clone encoding phenylalanine ammonia-lyase (PAL) was isolated from R. sachalinensis using rapid amplification of cDNA ends. The resulting cDNA was designated PALrs1. It is 2407-bp long and encodes 710 deduced amino acid residues. Southern blot analysis of genomic DNA indicated that the PAL gene family is composed of three to five genes in the R. sachalinensis genome. Northern blot analysis revealed that transcripts of PALrs1 were present in calli, leaves and stems, but expression in roots was very low. The PALrs1 under the 35S promoter with double-enhancer sequences from CaMV-Omega and TMV-Omega fragments was transferred into R. sachalinensis via Agrobacterium tumefaciens. PCR and PCR-Southern blot confirmed that the PALrs1 gene had been integrated into the genome of transgenic plants. Northern blot analysis revealed that the PALrs1 gene had been expressed at the transcriptional level. High-performance liquid chromatography indicated that overexpression of the PALrs1 gene resulted in a 3.3-fold increase in p-coumaric acid content, as expected. In contrast, levels of tyrosol and salidroside were 4.7-fold and 7.7-fold, respectively, lower in PALrs1 transgenic plants than in controls. Furthermore, overexpression of the PALrs1 gene resulted in a 2.6-fold decrease in tyrosine content. These data suggest that overexpression of the PALrs1 gene and accumulation of p-coumaric acid did not facilitate tyrosol biosynthesis; tyrosol, as a phenylethanoid derivative, is not derived from phenylalanine; and reduced availability of tyrosine most likely resulted in a large reduction in tyrosol biosynthesis and accumulation of salidroside.

Isoenzyme markers and polyacrylamide gel electrophoresis have been used to study the genetic structure of populations of Rhodiola iremelica Boriss. (Grassulaceae), a Southern Ural endemic protected by the state and included in the Red Data Book of Bashkortostan Republic. A relatively large genetic variation at the species level has been found. The subdivision among populations (F(ST) = 0.115) is higher than in most cross-pollination angiosperms. No consistent pattern has been observed in the spatial distribution of its genetic variation. The relatively high differentiation among samples of R. iremelica characterized by small effective population sizes, may be accounted for by genetic drift, inbreeding, and a restricted gene flow. To preserve the population gene pool, in situ protection of the species in nature is insufficient. It seems advisable to create synthetic populations ex situ and reintroduce them into nature.

Salidroside is a novel effective adaptogenic drug extracted from the medicinal plant Rhodiola sachalinensis A. Bor. Because this plant is a rare resource and has low yield, there is great interest in enhancing the production of salidroside. In this study, a putative UDP-glucosyltransferase (UGT) cDNA, UGT73B6, was isolated from Rhodiola sachalinensis using a rapid amplification of cDNA ends (RACE) method. The cDNA was 1,598 bp in length encoding 480 deduced amino acid residues with a conserved UDP-glucose-binding domain (PSPG box). Southern blot analysis of genomic DNA indicated that UGT73B6 existed as a single copy gene in the R. sachalinensis genome. Northern blot analysis revealed that transcripts of UGT73B6 were present in roots, calli and stems, but not in leaves. The UGT73B6 under 35S promoter with double-enhancer sequences from CaMV-Omega and TMV-Omega fragments was transferred into R. sachalinensis via Agrobacterium tumefaciens. PCR, PCR-Southern and Southern blot analyses confirmed that the UGT73B6 gene had been integrated into the genome of transgenic calli and plants. Northern blot analysis revealed that the UGT73B6 gene had been expressed at the transcriptional level. High performance liquid chromatography (HPLC) analysis indicated that the overexpression of the UGT73B6 gene resulted in an evident increase of salidroside content. These data suggest that the cloned UGT73B6 can regulate the conversion of tyrosol aglycon to salidroside in R. sachalinensis. This is the first cloned glucosyltransferase gene involved in salidroside biosynthesis.