Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong.

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Structure basis of the caffeic acid O-methyltransferase from Ligusiticum chuanxiong to understand its selective mechanism.

Caffeic acid O-methyltransferase from Ligusticum chuanxiong (LcCOMT) showed strict regiospecificity despite a relative degree of preference. Compared with caffeic acid, methyl caffeate was the preferential substrate by its low Km and high Kcat. In this study, we obtained the SAM binary (1.80 Å) and SAH binary (1.95 Å) complex LcCOMT crystal structures, and established the ternary complex structure with methyl caffeate by molecular docking. The active site of LcCOMT included phenolic substrate pocket, SAM/SAH ligand pocket and conserved catalytic residues as well. The regiospecificity of LcCOMT that permitted only 3-hydroxyl group to be methylated arise from the interactions between the active site and the phenyl ring. However, the propanoid tail governed the relative preference of LcCOMT. The ester group in methyl caffeate stabilized the anionic intermediate caused by His268-Asp269 pair, whereas caffeic acid was unable to stabilize the anionic intermediate due to the adjacent carboxylate anion in the propanoid tail. Ser183 residue formed an additional hydrogen bond with SAH and its role was identified by S183A mutation. Ile318 residue might be a potential site for determination of substrate preference, and its mutation led to the change of tertiary conformation. The results supported the selective mechanism of LcCOMT.

Uncovering dynamic evolution in the plastid genome of seven Ligusticum species provides insights into species discrimination and phylogenetic implications.

Ligusticum L., one of the largest members in Apiaceae, encompasses medicinally important plants, the taxonomic statuses of which have been proved to be difficult to resolve. In the current study, the complete chloroplast genomes of seven crucial plants of the best-known herbs in Ligusticum were presented. The seven genomes ranged from 148,275 to 148,564 bp in length with a highly conserved gene content, gene order and genomic arrangement. A shared dramatic decrease in genome size resulted from a lineage-specific inverted repeat (IR) contraction, which could potentially be a promising diagnostic character for taxonomic investigation of Ligusticum, was discovered, without affecting the synonymous rate. Although a higher variability was uncovered in hotspot divergence regions that were unevenly distributed across the chloroplast genome, a concatenated strategy for rapid species identification was proposed because separate fragments inadequately provided variation for fine resolution. Phylogenetic inference using plastid genome-scale data produced a concordant topology receiving a robust support value, which revealed that L. chuanxiong had a closer relationship with L. jeholense than L. sinense, and L. sinense cv. Fuxiong had a closer relationship to L. sinense than L. chuanxiong, for the first time. Our results not only furnish concrete evidence for clarifying Ligusticum taxonomy but also provide a solid foundation for further pharmaphylogenetic investigation.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships.

BACKGROUND: The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach. RESULTS: Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses. CONCLUSION: The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Sequencing and Comparative Analysis of the Chloroplast Genome of Angelica polymorpha and the Development of a Novel Indel Marker for Species Identification.

The genus Angelica (Apiaceae) comprises valuable herbal medicines. In this study, we determined the complete chloroplast (CP) genome sequence of A. polymorpha and compared it with that of Ligusticum officinale (GenBank accession no. NC039760). The CP genomes of A. polymorpha and L. officinale were 148,430 and 147,127 bp in length, respectively, with 37.6% GC content. Both CP genomes harbored 113 unique functional genes, including 79 protein-coding, four rRNA, and 30 tRNA genes. Comparative analysis of the two CP genomes revealed conserved genome structure, gene content, and gene order. However, highly variable regions, sufficient to distinguish between A. polymorpha and L. officinale, were identified in hypothetical chloroplast open reading frame1 (ycf1) and ycf2 genic regions. Nucleotide diversity (Pi) analysis indicated that ycf4⁻chloroplast envelope membrane protein (cemA) intergenic region was highly variable between the two species. Phylogenetic analysis revealed that A. polymorpha and L. officinale were well clustered at family Apiaceae. The ycf4-cemA intergenic region in A. polymorpha carried a 418 bp deletion compared with L. officinale. This region was used for the development of a novel indel marker, LYCE, which successfully discriminated between A. polymorpha and L. officinale accessions. Our results provide important taxonomic and phylogenetic information on herbal medicines and facilitate their authentication using the indel marker.

Identification of Species in the Aromatic Spice Family Apiaceae Using DNA Mini-barcodes.

The species of the aromatic plant family Apiaceae are mainly used as spices and foods, but the family also includes medicinal and some poisonous plant species. Due to the similar chemical compounds or aroma and morphology, the poisonous species are often mistaken for the edible aromatic species. It is therefore imperative to correctly identify the species present at the initial raw stage samples to ensure product safety and efficacy. At the molecular level, plant species can be identified using DNA loci either from nuclear or plastid genome with easily available universal oligonucleotides, a technique called DNA barcoding. However, this is possible when single-species plant material is present but may not work on a mixture of plants species. Another disadvantage is that using universal oligonucleotides is of limited help, especially if the adulterating material is present in low quantities. On the other hand, if using the species-specific oligonucleotides, only single specific adulterating plant material could be detected and, consequently, the unexpected adulterants may go undetected. Therefore, in the current work, four degenerated oligonucleotides from ITS1 and ITS2 regions of the nuclear genome were designed that can bind to a variety of Apiaceae genera only and not to other genera belonging to different plant families. These family-specific oligonucleotides were able to amplify a diagnostic PCR product from 16 Apiaceae species that, upon sequencing, revealed the identity of the plant it was derived from. The size of these products is around 140 bp for ITS1 and approximately 80 bp for the ITS2 region. The size range of the amplified products falls in the category of a desired mini-barcode size to be used for damaged/fragmented DNA and next generation sequencing.

[Analysis of endophytic fungi community of Ligusticum chuanxiong using PCR-DGGE].

OBJECTIVE: To research the diversity of endophytic fungal communities among Ligusticum chuanxiong growing at 5 areas in Sichuan province, and illuminate the developing mechanism of geoherbs from the microecological perspective. METHOD: The PCR-DGGE and DNA sequencing techniques were used to analyze the endophytic fungi community of L. chuanxiong. RESULT: The community of endophytic fungi present difference among different growing areas. Though minor difference were found among individuals at the same area, similarity among individuals from the same growing areas were higher significantly than those from different growing areas. Compared with the other 4 growing areas, L. chuanxiong from Shiyang town, Dujiangyan city had more abundant endophytic fungi and low similarity to others, and which probably had special types of fungi. CONCLUSION: The abundant and stable endophytic fungal community is an important factor for the development of geoherb L. chuanxiong at Shiyang town, Dujiangyan city.

Metabolic profiling and identification of the genetic varieties and agricultural origin of Cnidium officinale and Ligusticum chuanxiong.

Quality control methods for Cnidium officinale and Ligusticum chuanxiong are lacking because their quality is influenced by multiple factors. Thus, there is a need to develop a multifactorial method for measuring quality that is both standardized and practical. Here, we report a profiling method based on gas chromatography-mass spectrometry (GC-MS) to discriminate among the genetic varieties and agricultural origins of C. officinale and L. chuanxiong. Our metabolome analysis identified 68 metabolites, 13 of which were newly identified in our samples. The S-plot of the OPLS discriminant analysis enabled us to determine significant biomarkers. Using only double-compound biomarkers, the samples were successfully classified into distinct groups defined by genetic variety and cultivation origin. This method will simplify the process of searching for quality control markers that can be used to determine genetic variety and agricultural origin.

Authentication of the oriental medicinal plant Ligusticum tenuissimum (Nakai) Kitagawa (Korean Go-Bon) by multiplex PCR.

The oriental medicinal plant Ligusticum tenuissimum (Korean name, Go-Bon) is widely used in Korea and China. L. tenuissimum (Go-Bon) has been employed in the treatment of headache and common cold, and as a fever remedy. The internal transcribed spacer (ITS) region was sequenced from thirty-four Go-Bon samples collected from botanical gardens and markets in Korea and China to identify and authenticate L. tenuissimum. Based on the ITS sequences, the thirty-four Go-Bon samples were classified into three groups: L. tenuissimum (Korean Go-Bon), L. jeholense (Chinese Go-Bon), and unknown Chinese Ligusticum species. Three specific primers were designed to identify the three groups of Ligusticum species using multiplex PCR. The established multiplex-PCR was proved to be effective for the differentiation of L. tenuissimum in commercial plant materials.

Molecular identification of "Chuanxiong" by nucleotide sequence and multiplex single base extension analysis on chloroplast trnK gene.

Chloroplast trnK gene sequences of Cnidium officinale and Ligusticum chuanxiong were determined to establish an effective method for identifying Japanese Senkyu and Chinese Chuanxiong, the two which have the same drug name in Chinese characters, similar external feature, but different botanical origins. Three sites of nucleotide differences were found between these 2 species at positions 767,924 and 964 from upstream in trnK gene sequence, allowing molecular identification of the two plants and crude drugs. Further, three kinds of specific primers of 14 mer, 23 mer and 30 mer long were designed to detect these 3 sites of marker nucleotides. By using multiplex single base extension (MSBE) analysis with the 3 specific primers, C. officinale and L. chuanxiong could be distinguished clearly by the electrophoretograms, where 3 peaks with different color of ddTMP, ddCMP and ddTMP were observed in case of C. officinale and those of ddGMP, ddAMP and ddGMP in L. chuanxiong. Moreover, trnK gene sequence of "Dongxiong," a kind of Chuanxiong cultivated in Northeast China, suggested that its botanical origin was C. officinale.

[matK and its nucleotide sequencing of crude drug chuanxiong and phylogenetic relationship between their species from China and Japan].

AIM: To provide more molecular evidences for species relationship between Chuanxiong (Ligusticum chuanxiong Hort.) from China and Japanese Chuanxiong (Senkyu in Japanese) (Cnidium officinale Makino). METHODS: To sequence such two genes as internal transcribed spacer (ITS) from nuclear rDNA and maturase for lysine (matK) in tRNA(lys) (UUU) intron from chloroplast DNA of both Ligusticum chuanxiong and Cnidium officinale using PCR direct sequencing and to analyze the sequence variation of two genes between these two species. RESULTS: The matK gene sequence of Ligusticum chuanxiong and Cnidium officinale is 1268 bp in length, coding 422 amino acids of maturase protein. ITS gene sequence 699 bp, consisting of 54 bp of 18S rRNA-3', 215 bp of ITS1, 162 bp of 5.8S rRNA, 222 bp of ITS2, 46 bp of 26S rRNA-5'. Multiple sequence alignment shows that the sequence of two genes between dried crude drug and fresh voucher material of Ligusticum chuanxiong and Cnidium officinale, there is 1 variable site (T-->C) in matK (upstream at 595 nt) and ITS (ITS1 at 54 nt) between Ligusticum chuanxiong and Cnidium officinale. CONCLUSION: Based on homology analysis of two genes plastid matK and nuclear ITS, the origin of Chuanxiong from China and Japan ought to be identical, the scientific name Cnidium officinale of Japanese Chuanxiong should be changed to Ligusticum chuanxiong.