Abstract: Anthraquinones are not only the main active constituents but also the index components for the quality control of Rhei Radix et Rhizoma. To study the anthraquinone biosynthesis, Rheum palmatum L. seedlings were subjected to a high-throughput transcriptomic sequencing analysis by Illumina HiSeq^TM 2000 150PE. The Illumina sequencing generated a total of 11.04 G clean data resulting in 736 309 74 clean reads, deposited in the sequence read archive (SRA accession SRP160030). Trinity do novo assembly yielded 93 646 unigenes, with an average of 1 108 nt. Functional annotation revealed that all unigenes were successfully annotated in the NR, NT, Swiss-port, PFAM, and KOG databases. GO enrichments showed that 57 subgroups were involved in biological process, cellular component, and molecular function. KEGG analysis indicated that 1 107 unigenes were implicated in 19 standard secondary metabolic pathways. 172 unigenes were analyzed to encode 28 key enzymes during the MVA, MEP, shikimic acid, and polyketide pathways related to anthraquinone biosynthesis. 125 CYP450 and 73 UGTs unigenes were related the modification of secondary metabolites in R. palmatum L. Furthermore, seven unigenes with full length cDNAs were successfully verified by RT-PCR and sequencing analyses. Then, MISA prediction produced a number of 18 885 simple sequence repeats (SSRs). Herein, the transcriptomic gene expression profiles of R. palmatum L. and candidate genes during the anthraquinone biosynthesis pathway were obtained for the first time. The results provided basic information for subsequent gene function characterization, secondary metabolic pathway analysis, and anthraquinone biosynthesis and regulation elucidation in R. palmatum L.