Background: Epimedium sagittatum (Sieb. Et Zucc.) Maxim, a traditional Chinese medicinal plant species, has been used extensively as genuine medicinal materials. Certain Epimedium species are endangered due to commercial overexploition, while sustainable application studies, conservation genetics, systematics, and marker-assisted selection (MAS) of Epimedium is less-studied due to the lack of molecular markers. Here, we report a set of expressed sequence tags (ESTs) and simple sequence repeats (SSRs) identified in these ESTs for E. sagittatum. Results: cDNAs of E. sagittatum are sequenced using 454 GS-FLX pyrosequencing technology. The raw reads are cleaned and assembled into a total of 76,459 consensus sequences comprising of 17,231 contigs and 59,228 singlets. About 38.5% (29,466) of the consensus sequences significantly match to the non-redundant protein database (E-value < 1e-10), 22,295 of which are further annotated using Gene Ontology (GO) terms. A total of 2,810 EST-SSRs is identified from the Epimedium EST dataset. Trinucleotide SSR is the dominant repeat type (55.2%) followed by dinucleotide (30.4%), tetranuleotide (7.3%), hexanucleotide (4.9%), and pentanucleotide (2.2%) SSR. The dominant repeat motif is AAG/CTT (23.6%) followed by AG/CT (19.3%), ACC/GGT (11.1%), AT/AT (7.5%), and AAC/GTT (5.9%). Thirty-two SSR-ESTs are randomly selected and primer pairs are synthesized for testing the transferability across 52 Epimedium species. Eighteen primer pairs (85.7%) could be successfully transferred to Epimedium species and sixteen of those show high genetic diversity with 0.35 of observed heterozygosity (Ho) and 0.65 of expected heterozygosity (He) and high number of alleles per locus (11.9). Conclusion: A large EST dataset with a total of 76,459 consensus sequences is generated, aiming to provide sequence information for deciphering secondary metabolism, especially for flavonoid pathway in Epimedium. A total of 2,810 EST-SSRs is identified from EST dataset and ~1580 EST-SSR markers are transferable. E. sagittatum EST-SSR transferability to the major Epimedium germplasm is up to 85.7%. Therefore, this EST dataset and EST-SSRs will be a powerful resource for further studies such as taxonomy, molecular breeding, genetics, genomics, and secondary metabolism in Epimedium species.

Tanshinones are the bioactivenor-diterpenoid constituents of the Chinese medicinal herb Danshen (Salvia miltiorrhiza). These groups of chemicals have the characteristic furan D-ring, which differentiates them from the phenolic abietane-type diterpenoids frequently found in the Lamiaceae family. However, how the 14,16-epoxy is formed has not been elucidated. Here, we report an improved genome assembly of Danshen using a highly homozygous genotype. We identify a cytochrome P450 (CYP71D) tandem gene array through gene expansion analysis. We show that CYP71D373 and CYP71D375 catalyze hydroxylation at carbon-16 (C16) and 14,16-ether (hetero)cyclization to form the D-ring, whereas CYP71D411 catalyzes upstream hydroxylation at C20. In addition, we discover a large biosynthetic gene cluster associated with tanshinone production. Collinearity analysis indicates a more specific origin of tanshinones inSalviagenus. It illustrates the evolutionary origin of abietane-type diterpenoids and those with a furan D-ring in Lamiaceae.

Secondary metabolites from plants play key roles in human medicine and chemical industries. Due to limited accumulation of secondary metabolites in plants and their important roles, characterization of key enzymes involved in biosynthetic pathway will enable metabolic engineering or synthetic biology to improve or produce the compounds in plants or microorganisms, which provides an alternative for production of these valuable compounds. Salvia miltiorrhiza, containing tanshinones and phenolic acids as its active compounds, has been widely used for the treatment of cardiovascular and cerebrovascular diseases. The biosynthetic analysis of secondary metabolites in S. miltiorrhiza has made great progress due to the successful genetic transformation system, simplified hairy roots system, and high-throughput sequencing. The cloned genes in S. miltiorrhiza had provided references for functional characterization of the post-modification steps involved in biosynthesis of tanshinones and phenolic acids, and further utilization of these steps in metabolic engineering. The strategies used in these studies could provide solid foundation for elucidation of biosynthetic pathways of diterpenoids and phenolic acids in other species. The present review systematically summarizes recent advances in biosynthetic pathway analysis of tanshinones and phenolic acids as well as synthetic biology and metabolic engineering applications of the rate-limiting genes involved in the secondary metabolism in S. miltiorrhiza.

Protoberberine alkaloids and benzophenanthridine alkaloids (BZDAs) are subgroups of benzylisoquinoline alkaloids (BIAs), which represent a diverse class of plant-specialized natural metabolites with many pharmacological properties. Microbial biosynthesis has been allowed for accessibility and scalable production of high-value BIAs. Here, we engineer Saccharomyces cerevisiae to de novo produce a series of protoberberines and BZDAs, including palmatine, berberine, chelerythrine, sanguinarine and chelirubine. An ER compartmentalization strategy is developed to improve vacuole protein berberine bridge enzyme (BBE) activity, resulting in >200% increase on the production of the key intermediate (S)-scoulerine. Another promiscuous vacuole protein dihydrobenzophenanthridine oxidase (DBOX) has been identified to catalyze two-electron oxidation on various tetrahydroprotoberberines at N7-C8 position and dihydrobenzophenanthridine alkaloids. Furthermore, cytosolically expressed DBOX can alleviate the limitation on BBE. This study highlights the potential of microbial cell factories for the biosynthesis of a diverse group of BIAs through engineering of heterologous plant enzymes.