Comparative genome/transcriptome analysis probes Boraginales' phylogenetic position, WGDs in Boraginales, and key enzyme genes in the alkannin/shikonin core pathway

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Genome Assembly of Cordia subcordata, a Coastal Protection Species in Tropical Coral Islands

Cordia subcordata trees or shrubs, belonging to the Boraginaceae family, have strong resistance and have adapted to their habitat on a tropical coral island in China, but the lack of genome information regarding its genetic background is unclear. In this study, the genome was assembled using both short/long whole genome sequencing reads and Hi-C reads. The assembled genome was 475.3 Mb, with 468.7 Mb (99.22%) of the sequences assembled into 16 chromosomes. Repeat sequences accounted for 54.41% of the assembled genome. A total of 26,615 genes were predicted, and 25,730 genes were functionally annotated using different annotation databases. Based on its genome and the other 17 species, phylogenetic analysis using 336 single-copy genes obtained from ortholog analysis showed that C. subcordata was a sister to Coffea eugenioides, and the divergence time was estimated to be 77 MYA between the two species. Gene family evolution analysis indicated that the significantly expanded gene families were functionally related to chemical defenses against diseases. These results can provide a reference to a deeper understanding of the genetic background of C. subcordata and can be helpful in exploring its adaptation mechanism on tropical coral islands in the future.

Genome-Wide Comparison and Functional Characterization of HMGR Gene Family Associated with Shikonin Biosynthesis in Lithospermum erythrorhizon

3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), as the rate-limiting enzyme in the mevalonate pathway, is essential for the biosynthesis of shikonin in Lithospermum erythrorhizon. However, in the absence of sufficient data, the principles of a genome-wide in-depth evolutionary exploration of HMGR family members in plants, as well as key members related to shikonin biosynthesis, remain unidentified. In this study, 124 HMGRs were identified and characterized from 36 representative plants, including L. erythrorhizon. Vascular plants were found to have more HMGR family genes than nonvascular plants. The phylogenetic tree revealed that during lineage and species diversification, the HMGRs evolved independently and intronless LerHMGRs emerged from multi-intron HMGR in land plants. Among them, Pinus tabuliformis and L. erythrorhizon had the most HMGR gene duplications, with 11 LerHMGRs most likely expanded through WGD/segmental and tandem duplications. In seedling roots and M9 cultured cells/hairy roots, where shikonin biosynthesis occurs, LerHMGR1 and LerHMGR2 were expressed significantly more than other genes. The enzymatic activities of LerHMGR1 and LerHMGR2 further supported their roles in catalyzing the conversion of HMG-CoA to mevalonate. Our findings provide insight into the molecular evolutionary properties and function of the HMGR family in plants and a basis for the genetic improvement of efficiently produced secondary metabolites in L. erythrorhizon.

The discovery and characterization of AeHGO in the branching route from shikonin biosynthesis to shikonofuran in Arnebia euchroma

Shikonin derivatives are natural naphthoquinone compounds and the main bioactive components produced by several boraginaceous plants, such as Lithospermum erythrorhizon and Arnebia euchroma. Phytochemical studies utilizing both L. erythrorhizon and A. euchroma cultured cells indicate the existence of a competing route branching out from the shikonin biosynthetic pathway to shikonofuran. A previous study has shown that the branch point is the transformation from (Z)-3''-hydroxy-geranylhydroquinone to an aldehyde intermediate (E)-3''-oxo-geranylhydroquinone. However, the gene encoding the oxidoreductase that catalyzes the branch reaction remains unidentified. In this study, we discovered a candidate gene belonging to the cinnamyl alcohol dehydrogenase family, AeHGO, through coexpression analysis of transcriptome data sets of shikonin-proficient and shikonin-deficient cell lines of A. euchroma. In biochemical assays, purified AeHGO protein reversibly oxidized (Z)-3''-hydroxy-geranylhydroquinone to produce (E)-3''-oxo-geranylhydroquinone followed by reversibly reducing (E)-3''-oxo-geranylhydroquinone to (E)-3''-hydroxy-geranylhydroquinone, resulting in an equilibrium mixture of the three compounds. Time course analysis and kinetic parameters showed that the reduction of (E)-3''-oxo-geranylhydroquinone was stereoselective and efficient in presence of NADPH, which determined that the overall reaction proceeded from (Z)-3''-hydroxy-geranylhydroquinone to (E)-3''-hydroxy-geranylhydroquinone. Considering that there is a competition between the accumulation of shikonin and shikonofuran derivatives in cultured plant cells, AeHGO is supposed to play an important role in the metabolic regulation of the shikonin biosynthetic pathway. Characterization of AeHGO should help expedite the development of metabolic engineering and synthetic biology toward production of shikonin derivatives.

Excretion of triacylglycerol as a matrix lipid facilitating apoplastic accumulation of a lipophilic metabolite shikonin

Plants produce a large variety of lipophilic metabolites, many of which are secreted by cells and accumulated in apoplasts. These compounds often play a role to protect plants from environmental stresses. However, little is known about how these lipophilic compounds are secreted into apoplastic spaces. In this study, we used shikonin-producing cultured cells of Lithospermum erythrorhizon as an experimental model system to analyze the secretion of lipophilic metabolites, taking advantage of its high production rate and the clear inducibility in culture. Shikonin derivatives are lipophilic red naphthoquinone compounds that accumulate exclusively in apoplastic spaces of these cells and also in the root epidermis of intact plants. Microscopic analysis showed that shikonin is accumulated in the form of numerous particles on the cell wall. Lipidomic analysis showed that L. erythrorhizon cultured cells secrete an appreciable portion of triacylglycerol (24-38% of total triacylglycerol), composed predominantly of saturated fatty acids. Moreover, in vitro reconstitution assay showed that triacylglycerol encapsulates shikonin derivatives with phospholipids to form lipid droplet-like structures. These findings suggest a novel role for triacylglycerol as a matrix lipid, a molecular component involved in the secretion of specialized lipophilic metabolites.

Comparison of Boraginales Plastomes: Insights into Codon Usage Bias, Adaptive Evolution, and Phylogenetic Relationships

The Boraginales (Boraginaceae a.l.) comprise more than 2450 species worldwide. However, little knowledge exists of the characteristics of the complete plastid genome. In this study, three new sequences representing the first pt genome of Heliotropiaceae and Cordiaceae were assembled and compared with other Boraginales species. The pt genome sizes of Cordia dichotoma, Heliotropium arborescens, and Tournefortia montana were 151,990 bp, 156,243 bp, and 155,891 bp, respectively. Multiple optimal codons were identified, which may provide meaningful information for enhancing the gene expression of Boraginales species. Furthermore, codon usage bias analyses revealed that natural selection and other factors may dominate codon usage patterns in the Boraginales species. The boundaries of the IR/LSC and IR/SSC regions were significantly different, and we also found a signal of obvious IR region expansion in the pt genome of Nonea vesicaria and Arnebia euchroma. Genes with high nucleic acid diversity (pi) values were also calculated, which may be used as potential DNA barcodes to investigate the phylogenetic relationships in Boraginales. psaI, rpl33, rpl36, and rps19 were found to be under positive selection, and these genes play an important role in our understanding of the adaptive evolution of the Boraginales species. Phylogenetic analyses implied that Boraginales can be divided into two groups. The existence of two subfamilies (Lithospermeae and Boragineae) in Boraginaceae is also strongly supported. Our study provides valuable information on pt genome evolution and phylogenetic relationships in the Boraginales species.

Root-associated bacteria modulate the specialised metabolome of Lithospermum officinale L

Bacteria influence plant growth and development and therefore are attractive resources for applications in agriculture. However, little is known about the impact of these microorganisms on secondary metabolite (SM) production by medicinal plants. Here we assessed, for the first time, the effects of bacteria on the modulation of SM production in the medicinal plant Lithospermum officinale (Boraginaceae family) with a focus on the naphthoquinones alkannin/shikonin and their derivatives (A/Sd). The study was conducted in an in vitro cultivation system developed for that purpose, as well as in a greenhouse. Targeted and non-targeted metabolomics were performed, and expression of the gene PGT encoding for a key enzyme in the A/S biosynthesis pathway was evaluated with qPCR. Three strains, Chitinophaga sp. R-73072, Xanthomonas sp. R-73098 and Pseudomonas sp. R-71838 induced a significant increase of A/Sd in L. officinale in both systems, demonstrating the strength of our approach for screening A/Sd-inducing bacteria. The bacterial treatments altered other plant metabolites derived from the shikimate pathway as well. Our results demonstrate that bacteria influence the biosynthesis of A/Sd and interact with different metabolic pathways. This work highlights the potential of bacteria to increase the production of SM in medicinal plants and reveals new patterns in the metabolome regulation of L. officinale.

Transcriptional dynamics of Chitinophaga sp. strain R-73072-mediated alkannin/shikonin biosynthesis in Lithospermum officinale

Plants are colonized by a wide range of bacteria, several of which are known to confer benefits to their hosts such as enhancing plant growth and the biosynthesis of secondary metabolites (SMs). Recently, it has been shown that Chitinophaga sp. strain R-73072 enhances the production of alkannin/shikonin, SMs of pharmaceutical and ecological importance. However, the mechanisms by which this bacterial strain increases these SMs in plants are not yet understood. To gain insight into these mechanisms, we analyzed the molecular responses of Lithospermum officinale, an alkannin/shikonin producing member of Boraginaceae, to inoculation with R-73072 in a gnotobiotic system using comparative transcriptomics and targeted metabolite profiling of root samples. We found that R-73072 modulated the expression of 1,328 genes, of which the majority appeared to be involved in plant defense and SMs biosynthesis including alkannin/shikonin derivatives. Importantly, bacterial inoculation induced the expression of genes that predominately participate in jasmonate and ethylene biosynthesis and signaling, suggesting an important role of these phytohormones in R-73072-mediated alkannin/shikonin biosynthesis. A detached leaf bioassay further showed that R-73072 confers systemic protection against Botrytis cinerea. Finally, R-73072-mediated coregulation of genes involved in plant defense and the enhanced production of alkannin/shikonin esters further suggest that these SMs could be important components of the plant defense machinery in alkannin/shikonin producing species.

Inventory of ATP-binding cassette proteins in Lithospermum erythrorhizon as a model plant producing divergent secondary metabolites

ATP-binding cassette (ABC) proteins are the largest membrane transporter family in plants. In addition to transporting organic substances, these proteins function as ion channels and molecular switches. The development of multiple genes encoding ABC proteins has been associated with their various biological roles. Plants utilize many secondary metabolites to adapt to environmental stresses and to communicate with other organisms, with many ABC proteins thought to be involved in metabolite transport. Lithospermum erythrorhizon is regarded as a model plant for studying secondary metabolism, as cells in culture yielded high concentrations of meroterpenes and phenylpropanoids. Analysis of the genome and transcriptomes of L. erythrorhizon showed expression of genes encoding 118 ABC proteins, similar to other plant species. The number of expressed proteins in the half-size ABCA and full-size ABCB subfamilies was ca. 50% lower in L. erythrorhizon than in Arabidopsis, whereas there was no significant difference in the numbers of other expressed ABC proteins. Because many ABCG proteins are involved in the export of organic substances, members of this subfamily may play important roles in the transport of secondary metabolites that are secreted into apoplasts.

Comparative analysis of the main medicinal substances and applications of Echium vulgare L. and Echium plantagineum L.: A review

ETHNOPHARMACOLOGICAL RELEVANCE

Echium vulgare L. and Echium plantagineum L. originated in the Mediterranean, and were later domesticated in Africa, America, Asia, Europe and Oceania, where they were widely used to treat many diseases including cough, urinary tract infection, fever, inflammation and muscle strain.

AIM OF THE STUDY

The purpose of this review is to provide scientific literature on the traditional uses, bioactive chemical components and pharmacological activities of two species of Echium, and to critically analyze the information provided, so as to understand the current work on these two species and explore the possible prospect of this plant in pharmaceutical research.

METHODS

Systematic review and meta-analysis were conducted according to Prisma guidelines, and the related literatures searched on Google Academic, Science Direct, Baidu Scholars and China National Knowledge Infrastructure (CNKI) up to June 2021 were reviewed. The key words used are: Echium, E.vulgare, E.plantagineum, plant components, chemical components, pharmacological activities, pharmaceutical products and applications. Thereafter all eligible studies are analyzed and summarized in this review. The selection of manuscripts is based on the following inclusion criteria: the article has years of research or publication, is published in English, Portuguese or Spanish and Chinese, and there are keywords in the title, abstract, keywords or full text of the article. For the selection of manuscripts, first, select articles according to titles, then summarize them, and finally, analyze the full text of the publication. Elimination criteria: 1. Duplicate reports; 2. There are research design defects and poor quality; 3. Incomplete data and unclear ending effect; 4. The statistical method is wrong and cannot be corrected.

RESULTS

The pharmacological characteristics of E.vulgare and E.plantagineum can basically support their traditional use, but the medicinal substances contained in them are quite different in composition and content, and the development and application of corresponding products are also different.

CONCLUSIONS

At present, there is little clinical data about drugs related to the two species, and more research is needed in the future, especially human experiments and clinical trials, to evaluate the cellular and molecular mechanisms based on pharmacological, biological activity and safety studies, and to provide more powerful scientific basis for their traditional medicinal properties. In addition, the further application and development of the medicinal products of E.vulgare and E.plantagineum still need to be precise and identified, so as to give full play to their medicinal potential.

Integrative analysis of the shikonin metabolic network identifies new gene connections and reveals evolutionary insight into shikonin biosynthesis

Plant specialized 1,4-naphthoquinones present a remarkable case of convergent evolution. Species across multiple discrete orders of vascular plants produce diverse 1,4-naphthoquinones via one of several pathways using different metabolic precursors. Evolution of these pathways was preceded by events of metabolic innovation and many appear to share connections with biosynthesis of photosynthetic or respiratory quinones. Here, we sought to shed light on the metabolic connections linking shikonin biosynthesis with its precursor pathways and on the origins of shiknoin metabolic genes. Downregulation of Lithospermum erythrorhizon geranyl diphosphate synthase (LeGPPS), recently shown to have been recruited from a cytoplasmic farnesyl diphosphate synthase (FPPS), resulted in reduced shikonin production and a decrease in expression of mevalonic acid and phenylpropanoid pathway genes. Next, we used LeGPPS and other known shikonin pathway genes to build a coexpression network model for identifying new gene connections to shikonin metabolism. Integrative in silico analyses of network genes revealed candidates for biochemical steps in the shikonin pathway arising from Boraginales-specific gene family expansion. Multiple genes in the shikonin coexpression network were also discovered to have originated from duplication of ubiquinone pathway genes. Taken together, our study provides evidence for transcriptional crosstalk between shikonin biosynthesis and its precursor pathways, identifies several shikonin pathway gene candidates and their evolutionary histories, and establishes additional evolutionary links between shikonin and ubiquinone metabolism. Moreover, we demonstrate that global coexpression analysis using limited transcriptomic data obtained from targeted experiments is effective for identifying gene connections within a defined metabolic network.

Shikonin, a naphthalene ingredient: Therapeutic actions, pharmacokinetics, toxicology, clinical trials and pharmaceutical researches

BACKGROUND

Shikonin is one of the major phytochemical components of Lithospermum erythrorhizon (Purple Cromwell), which is a type of medicinal herb broadly utilized in traditional Chinese medicine. It is well established that shikonin possesses remarkable therapeutic actions on various diseases, with the underlying mechanisms, pharmacokinetics and toxicological effects elusive. Also, the clinical trial and pharmaceutical study of shikonin remain to be comprehensively delineated.

PURPOSE

The present review aimed to systematically summarize the updated knowledge regarding the therapeutic actions, pharmacokinetics, toxicological effects, clinical trial and pharmaceutical study of shikonin.

METHODS

The information contained in this review article were retrieved from some authoritative databases including Web of Science, PubMed, Google scholar, Chinese National Knowledge Infrastructure (CNKI), Wanfang Database and so on, till August 2021.

RESULTS

Shikonin exerts multiple therapeutic efficacies, such as anti-inflammation, anti-cancer, cardiovascular protection, anti-microbiomes, analgesia, anti-obesity, brain protection, and so on, mainly by regulating the NF-κB, PI3K/Akt/MAPKs, Akt/mTOR, TGF-β, GSK3β, TLR4/Akt signaling pathways, NLRP3 inflammasome, reactive oxygen stress, Bax/Bcl-2, etc. In terms of pharmacokinetics, shikonin has an unfavorable oral bioavailability, 64.6% of the binding rate of plasma protein, and enhances some metabolic enzymes, particularly including cytochrome P450. In regard to the toxicological effects, shikonin may potentially cause nephrotoxicity and skin allergy. The above pharmacodynamics and pharmacokinetics of shikonin have been validated by few clinical trials. In addition, pharmaceutical innovation of shikonin with novel drug delivery system such as nanoparticles, liposomes, microemulsions, nanogel, cyclodextrin complexes, micelles and polymers are beneficial to the development of shikonin-based drugs.

CONCLUSIONS

Shikonin is a promising phytochemical for drug candidates. Extensive and intensive explorations on shikonin are warranted to expedite the utilization of shikonin-based drugs in the clinical setting.

Exploring the evolutionary process of alkannin/shikonin O-acyltransferases by a reliable Lithospermum erythrorhizon genome

Increasing genome data are coming out. Genome size estimation plays an essential role in guiding genome assembly. Several months ago, other researchers were the first to publish a draft genome of the red gromwell (i.e. Lithospermum erythrorhizon). However, we considered that the genome size they estimated and assembled was incorrect. This study meticulously estimated the L. erythrorhizon genome size to should be ∼708.74 Mb and further provided a reliable genome version (size ≈ 693.34 Mb; contigN50 length ≈ 238.08 Kb) to support our objection. Furthermore, according to our genome, we identified a gene family of the alkannin/shikonin O-acyltransferases (i.e. AAT/SAT) that catalysed enantiomer-specific acylations in the alkannin/shikonin biosynthesis (a characteristic metabolic pathway in L. erythrorhizon's roots) and further explored its evolutionary process. The results indicated that the existing AAT/SAT were not generated from only one round of gene duplication but three rounds; after different rounds of gene duplication, the existing AAT/SAT and their recent ancestors were under positive selection at different amino acid sites. These suggested that a combined power from gene duplication plus positive selection plausibly propelled AAT/SAT's functional differentiation in evolution.

Asterid Phylogenomics/Phylotranscriptomics Uncover Morphological Evolutionary Histories and Support Phylogenetic Placement for Numerous Whole-Genome Duplications

Asterids are one of the most successful angiosperm lineages, exhibiting extensive morphological diversity and including a number of important crops. Despite their biological prominence and value to humans, the deep asterid phylogeny has not been fully resolved, and the evolutionary landscape underlying their radiation remains unknown. To resolve the asterid phylogeny, we sequenced 213 transcriptomes/genomes and combined them with other data sets, representing all accepted orders and nearly all families of asterids. We show fully supported monophyly of asterids, Berberidopsidales as sister to asterids, monophyly of all orders except Icacinales, Aquifoliales, and Bruniales, and monophyly of all families except Icacinaceae and Ehretiaceae. Novel taxon placements benefited from the expanded sampling with living collections from botanical gardens, resolving hitherto uncertain relationships. The remaining ambiguous placements here are likely due to limited sampling and could be addressed in the future with relevant additional taxa. Using our well-resolved phylogeny as reference, divergence time estimates support an Aptian (Early Cretaceous) origin of asterids and the origin of all orders before the Cretaceous-Paleogene boundary. Ancestral state reconstruction at the family level suggests that the asterid ancestor was a woody terrestrial plant with simple leaves, bisexual, and actinomorphic flowers with free petals and free anthers, a superior ovary with a style, and drupaceous fruits. Whole-genome duplication (WGD) analyses provide strong evidence for 33 WGDs in asterids and one in Berberidopsidales, including four suprafamilial and seven familial/subfamilial WGDs. Our results advance the understanding of asterid phylogeny and provide numerous novel evolutionary insights into their diversification and morphological evolution.

CYP82AR Subfamily Proteins Catalyze C-1' Hydroxylations of Deoxyshikonin in the Biosynthesis of Shikonin and Alkannin

Shikonin and S-enantiomer alkannin are important naphthoquinone derivatives present in many Boraginaceae species. We report that cytochrome P450 monooxygenases (CYPs) from a new CYP82AR subfamily catalyzed hydroxylations of deoxyshikonin at C-1' position of isoprenoid side chain. Two homologues were discovered from each species of the four Boraginaceae plants. One CYP preferred converting deoxyshikonin into shikonin, and the other stereoselectively hydroxylated deoxyshikonin into alkannin. The discovery might be a general feature of shikonin/alkannin-producing Boraginaceae plants.

Convergent evolution of plant specialized 1,4-naphthoquinones: metabolism, trafficking, and resistance to their allelopathic effects

Plant 1,4-naphthoquinones encompass a class of specialized metabolites known to mediate numerous plant-biotic interactions. This class of compounds also presents a remarkable case of convergent evolution. The 1,4-naphthoquinones are synthesized by species belonging to nearly 20 disparate orders spread throughout vascular plants, and their production occurs via one of four known biochemically distinct pathways. Recent developments from large-scale biology and genetic studies corroborate the existence of multiple pathways to synthesize plant 1,4-naphthoquinones and indicate that extraordinary events of metabolic innovation and links to respiratory and photosynthetic quinone metabolism probably contributed to their independent evolution. Moreover, because many 1,4-naphthoquinones are excreted into the rhizosphere and they are highly reactive in biological systems, plants that synthesize these compounds also needed to independently evolve strategies to deploy them and to resist their effects. In this review, we highlight new progress made in understanding specialized 1,4-naphthoquinone biosynthesis and trafficking with a focus on how these discoveries have shed light on the convergent evolution and diversification of this class of compounds in plants. We also discuss how emerging themes in metabolism-based herbicide resistance may provide clues to mechanisms plants employ to tolerate allelopathic 1,4-naphthoquinones.

Alcohol Dehydrogenase Activity Converts 3″-Hydroxy-geranylhydroquinone to an Aldehyde Intermediate for Shikonin and Benzoquinone Derivatives in Lithospermum erythrorhizon

Shikonin derivatives are red naphthoquinone pigments produced by several boraginaceous plants, such as Lithospermum erythrorhizon. These compounds are biosynthesized from p-hydroxybenzoic acid and geranyl diphosphate. The coupling reaction that yields m-geranyl-p-hydroxybenzoic acid has been actively characterized, but little is known about later biosynthetic reactions. Although 3″-hydroxy-geranylhydroquinone produced from geranylhydroquinone by CYP76B74 has been regarded as an intermediate of shikonin derivatives, the next intermediate has not yet been identified. This study describes a novel alcohol dehydrogenase activity in L. erythrorhizon cell cultures. This enzyme was shown to oxidize the 3″-alcoholic group of (Z)-3″-hydroxy-geranylhydroquinone to an aldehyde moiety concomitant with the isomerization at the C2'-C3' double bond from the Z-form to the E-form. An enzyme oxidizing this substrate was not detected in other plant cell cultures, suggesting that this enzyme is specific to L. erythrorhizon. The reaction product, (E)-3″-oxo-geranylhydroquinone, was further converted to deoxyshikonofuran, another meroterpenoid metabolite produced in L. erythrorhizon cells. Although nonenzymatic cyclization occurred slowly, it was more efficient in the presence of crude enzymes of L. erythrorhizon cells. This activity was detected in both shikonin-producing and nonproducing cells, suggesting that the aldehyde intermediate at the biosynthetic branch point between naphthalene and benzo/hydroquinone ring formation likely constitutes a key common intermediate in the synthesis of shikonin and benzoquinone products, respectively.

Hybrid de novo genome assembly of red gromwell (Lithospermum erythrorhizon) reveals evolutionary insight into shikonin biosynthesis

Lithospermum erythrorhizon (red gromwell; zicao) is a medicinal and economically valuable plant belonging to the Boraginaceae family. Roots from L. erythrorhizon have been used for centuries based on the antiviral and wound-healing properties produced from the bioactive compound shikonin and its derivatives. More recently, shikonin, its enantiomer alkannin, and several other shikonin/alkannin derivatives have collectively emerged as valuable natural colorants and as novel drug scaffolds. Despite several transcriptomes and proteomes having been generated from L. erythrorhizon, a reference genome is still unavailable. This has limited investigations into elucidating the shikonin/alkannin pathway and understanding its evolutionary and ecological significance. In this study, we obtained a de novo genome assembly for L. erythrorhizon using a combination of Oxford Nanopore long-read and Illumina short-read sequencing technologies. The resulting genome is ∼367.41 Mb long, with a contig N50 size of 314.31 kb and 27,720 predicted protein-coding genes. Using the L. erythrorhizon genome, we identified several additional p-hydroxybenzoate:geranyltransferase (PGT) homologs and provide insight into their evolutionary history. Phylogenetic analysis of prenyltransferases suggests that PGTs originated in a common ancestor of modern shikonin/alkannin-producing Boraginaceous species, likely from a retrotransposition-derived duplication event of an ancestral prenyltransferase gene. Furthermore, knocking down expression of LePGT1 in L. erythrorhizon hairy root lines revealed that LePGT1 is predominantly responsible for shikonin production early in culture establishment. Taken together, the reference genome reported in this study and the provided analysis on the evolutionary origin of shikonin/alkannin biosynthesis will guide elucidation of the remainder of the pathway.