Two chromosome-level genome assemblies of Rhodiola shed new light on genome evolution in rapid radiation and evolution of the biosynthetic pathway of salidroside.

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.

High-level production of Rhodiola rosea characteristic component rosavin from D-glucose and L-arabinose in engineered Escherichia coli.

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.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences.

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.

Genomic data provide a robust phylogeny backbone for Rhodiola L. (Crassulaceae) and reveal extensive reticulate evolution during its rapid radiation.

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.

Morphometric Characteristics and Genetic Issr Marker Variability in Rhodiola rosea L. (Crassulaceae) in Different Ecological and Geographic Conditions in the Altai Republic.

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.

Comparative analyses of chloroplast genomes from Six Rhodiola species: variable DNA markers identification and phylogenetic relationships within the genus.

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.

Comparative phylogeography study reveals introgression and incomplete lineage sorting during rapid diversification of Rhodiola.

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.

Conservation and innovation: Plastome evolution during rapid radiation of Rhodiola on the Qinghai-Tibetan Plateau.

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 CO2 concentration and high-intensity light.

Regulation by fungal endophyte of Rhodiola crenulata from enzyme genes to metabolites based on combination of transcriptome and metabolome.

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.

Gone with the trees: Phylogeography of Rhodiola sect. Trifida (Crassulaceae) reveals multiple refugia on the Qinghai-Tibetan Plateau.

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.

De novo biosynthesis of Gastrodin in Escherichia coli.

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.

Analysis of the relationship between ribosomal DNA ITS sequences and active components in Rhodiola plants.

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 H4, H5, or H10, tyrosol was not detected with haplotypes H3, H5, H7, H10, H14, or H15, and gallic acid was detected in with all haplotypes except H14 and H15. 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.

[Effect of endophytic fungus on expression of key enzyme genes in pathway of salidroside biosynthesis in Rhodiola crenulata].

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.

Isolation of microsatellite loci for Rhodiola alsia (Crassulaceae), an important ethno-medicinal herb endemic to the Qinghai-Tibetan plateau.

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.

Rapid radiation and dispersal out of the Qinghai-Tibetan Plateau of an alpine plant lineage Rhodiola (Crassulaceae).

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.

Phylogeographical patterns of an alpine plant, Rhodiola dumulosa (Crassulaceae), inferred from chloroplast DNA sequences.

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.

Overexpression of Rhodiola crenulata glutathione peroxidase 5 increases cold tolerance and enhances the pharmaceutical value of the hairy roots.

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.

Genetic Diversity in Natural Populations of Rhodiola Species of Different Adaptation Strategies.

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.

Comparative transcriptome and tissue-specific expression analysis of genes reveal tissue-cultured plants as an alternative source for phenylethanoids and phenylpropanoids in Rhodiola imbricata (Edgew.).

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 N50 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.

Rhodiola activates macrophage migration inhibitory factor to alleviate non-alcoholic fatty liver disease.

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.