Comparative transcriptome analysis and HPLC reveal candidate genes associated with synthesis of bioactive constituents in Rheum palmatum complex

Content of bioactive constituents is one of the most important characteristics in Rheum palmatum complex. Increasing ingredient content through genetic breeding is an effective strategy to solve the contradiction between large market demand and resource depletion, but currently hampered by limited understanding of metabolite biosynthesis in rhubarb. In this study, deep transcriptome sequencing was performed to compare roots, stems, and leaves of two Rheum species (PL and ZK) that show different levels of anthraquinone contents. Approximately 0.52 billion clean reads were assembled into 58,782 unigenes, of which around 80% (46,550) were found to be functionally annotated in public databases. Expression patterns of differential unigenes between PL and ZK were thoroughly investigated in different tissues. This led to the identification of various differentially expressed genes (DEGs) involved in shikimate, MEP, MVA, and polyketide pathways, as well as those involved in catechin and gallic acid biosynthesis. Some structural enzyme genes were shown to be significantly up-regulated in roots of ZK with high anthraquinone content, implying potential central roles in anthraquinone synthesis. Taken together, our study provides insights for future functional studies to unravel the mechanisms underlying metabolite biosynthesis in rhubarb.