A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum.

BACKGROUND: Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism. RESULTS: We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening. CONCLUSION: Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

Transcriptomic divergence of the Rheum palmatum complex derived from top-geoherb and non-geoherb areas provides the insights into geoherbalism properties of rhubarb.

Geoherb usually represents high-quality medicinal herbs with better clinical therapeutic effects, and elucidating the geoherbalism is essential for the quality improvement of traditional Chinese Medicine. However, few researches were conducted to clarify the geoherbalism based on a large scale of transcriptomics. In the present study, we compared the transcriptomes of Rheum palmatum complex derived from top-geoherb and non-geoherb areas to show the geoherbalism properties of rhubarb. A total of 412.32 Gb clean reads were obtained with unigene numbers of 100,615 after assembly. Based on the obtained transcriptome datasets, key enzyme-encoding genes involved in the anthraquinones biosynthesis were also obtained. We also found that 21 anthraquinone-related unigenes were differentially expressed between two different groups, and some of these DEGs were correlated to the content accumulation of five free anthraquinones, indicating that the gene expression profiles may promote the geoherbalism formation of rhubarb. In addition, the selective pressure analyses indicated that most paired orthologous genes between these two groups were subject to negative selection, and only a low proportion of orthologs under positive selection were detected. Functional annotation analyses indicated that these positive-selected genes related to the functions such as gene expression, substance transport, stress response and metabolism, indicating that discrepant environment also enhanced the formation of geoherbalism. Our study not only provided insights for the genetic mechanism of geoherbalism of rhubarb, but also laid more genetic cues for the future rhubarb germplasms improvement and utilization.

Phylogenomics, phylogeography and germplasms authentication of the Rheum palmatum complex based on complete chloroplast genomes.

As a traditional Chinese medicine, rhubarb is used to treat several diseases such as severe acute pancreatitis, sepsis and chronic renal failure. However, few studies focused on the authentication of germplasm for the Rheum palmatum complex, and no studies have been conducted to elucidate the evolutionary history of the R. palmatum complex using plastome datasets. Hence, we aim to develop the potential molecular markers to identify the elite germplasms of rhubarb and explore the divergence and biogeographic history of the R. palmatum complex based on the newly sequenced chloroplast genome datasets. Chloroplast genomes of thirty-five the R. palmatum complex germplasms were sequenced, and the length ranged from 160,858 to 161,204 bp. The structure, gene content and gene order were highly conserved across all genomes. Eight InDels and sixty-one SNPs loci could be used to authenticate the high-quality germplasms of rhubarb in specific areas. Phylogenetic analysis revealed that all rhubarb germplasms were clustered in the same clade with high bootstrap support values and Bayesian posterior probabilities. According to the molecular dating result, the intraspecific divergence of the complex occurred in the Quaternary, which might be affected by climatic fluctuation. The biogeography reconstruction indicated that the ancestor of the R. palmatum complex might originate from the Himalaya-Hengduan Mountains or/and Bashan-Qinling Mountains, and then spread to surrounding areas. Several useful molecular markers were developed to identify rhubarb germplasms, and our study will provide further understanding on speciation, divergence and biogeography of the R. palmatum complex.

Comparative Transcriptome Analyses of Different Rheum officinale Tissues Reveal Differentially Expressed Genes Associated with Anthraquinone, Catechin, and Gallic Acid Biosynthesis.

Rheum officinale Baill. is an important traditional Chinese medicinal herb, its dried roots and rhizomes being widely utilized to cure diverse diseases. However, previous studies mainly focused on the active compounds and their pharmacological effects, and the molecular mechanism underlying the biosynthesis of these ingredients in R. officinale is still elusive. Here, we performed comparative transcriptome analyses to elucidate the differentially expressed genes (DEGs) in the root, stem, and leaf of R. officinale. A total of 236,031 unigenes with N50 of 769 bp was generated, 136,329 (57.76%) of which were annotated. A total of 5884 DEGs was identified after the comparative analyses of different tissues; 175 and 126 key enzyme genes with tissue-specific expression were found in the anthraquinone, catechin/gallic acid biosynthetic pathway, respectively, and some of these key enzyme genes were verified by qRT-PCR. The phylogeny of the PKS III family in Polygonaceae indicated that probably only PL_741 PKSIII1, PL_11549 PKSIII5, and PL_101745 PKSIII6 encoded PKSIII in the polyketide pathway. These results will shed light on the molecular basis of the tissue-specific accumulation and regulation of secondary metabolites in R. officinale, and lay a foundation for the future genetic diversity, molecular assisted breeding, and germplasm resource improvement of this essential medicinal plant.