Deciphering the high-quality genome sequence of coriander that causes controversial feelings

The gradual establishment of complex coumarin biosynthetic pathway in Apiaceae

Complex coumarins (CCs) represent characteristic metabolites found in Apiaceae plants, possessing significant medical value. Their essential functional role is likely as protectants against pathogens and regulators responding to environmental stimuli. Utilizing genomes and transcriptomes from 34 Apiaceae plants, including our recently sequenced Peucedanum praeruptorum, we conduct comprehensive phylogenetic analyses to reconstruct the detailed evolutionary process of the CC biosynthetic pathway in Apiaceae. Our results show that three key enzymes - p-coumaroyl CoA 2'-hydroxylase (C2'H), C-prenyltransferase (C-PT), and cyclase - originated successively at different evolutionary nodes within Apiaceae through various means of gene duplications: ectopic and tandem duplications. Neofunctionalization endows these enzymes with novel functions necessary for CC biosynthesis, thus completing the pathway. Candidate genes are cloned for heterologous expression and subjected to in vitro enzymatic assays to test our hypothesis regarding the origins of the key enzymes, and the results precisely validate our evolutionary inferences. Among the three enzymes, C-PTs are likely the primary determinant of the structural diversity of CCs (linear/angular), due to divergent activities evolved to target different positions (C-6 or C-8) of umbelliferone. A key amino acid variation (Ala161/Thr161) is identified and proven to play a crucial role in the alteration of enzymatic activity, possibly resulting in distinct binding forms between enzymes and substrates, thereby leading to different products. In conclusion, this study provides a detailed trajectory for the establishment and evolution of the CC biosynthetic pathway in Apiaceae. It explains why only a portion, not all, of Apiaceae plants can produce CCs and reveals the mechanisms of CC structural diversity among different Apiaceae plants.

Fusarium Wilt of Coriander: Root Cause Analysis and Varietal Tolerance Development

Since 2012, growers of coriander, Coriandrum sativum L., in Israel have been suffering from summer wilting that can result in entire fields collapsing. The current study aimed to determine the cause of the phenomenon and find a genetic solution to the problem. The disease was reproduced in a growth chamber using naturally-infested soil from a commercial field. Wilt became apparent within two weeks, and after ten weeks, all plants died compared to plants in sterilized soil from the same source. Fusarium oxysporum was isolated from infected plants, and Koch's postulates were completed. Sequence analysis of the Elongation Factor (EF1α) encoding gene of the pathogen had a 99.54% match to F. oxysporum f. sp. coriandrii. Several coriander varieties were screened for resistance or tolerance to the disease. In four independent experiments, only the cultivar 'Smadi' showed high tolerance, while other genotypes were susceptible. In a trial in a naturally infested field, the cultivar 'Smadi' outperformed the commercial cultivar 'Blair'. 'Smadi' provides a cropping solution to many Israeli farmers, yet this winter cultivar bolts early in the summer. There is a further need to characterize the tolerance mechanism and inheritance for informed breeding of late-bolting Fusarium-resistant coriander.

Genomic, transcriptomic, and metabolomic analyses provide insights into the evolution and development of a medicinal plant Saposhnikovia divaricata (Apiaceae)

Saposhnikovia divaricata, 2n = 2x = 16, as a perennial species, is widely distributed in China, Mongolia, Russia, etc. It is a traditional Chinese herb used to treat tetanus, rubella pruritus, rheumatic arthralgia, and other diseases. Here, we assembled a 2.07 Gb and N50 scaffold length of 227.67 Mb high-quality chromosome-level genome of S. divaricata based on the PacBio Sequel II sequencing platform. The total number of genes identified was 42 948, and 42 456 of them were functionally annotated. A total of 85.07% of the genome was composed of repeat sequences, comprised mainly of long terminal repeats (LTRs) which represented 73.7% of the genome sequence. The genome size may have been affected by a recent whole-genome duplication event. Transcriptional and metabolic analyses revealed bolting and non-bolting S. divaricata differed in flavonoids, plant hormones, and some pharmacologically active components. The analysis of its genome, transcriptome, and metabolome helped to provide insights into the evolution of bolting and non-bolting phenotypes in wild and cultivated S. divaricata and lays the basis for genetic improvement of the species.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, Ks, and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis

BACKGROUND

Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies.

RESULTS

A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins.

CONCLUSIONS

A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds.

The genome of the toxic invasive species Heracleum sosnowskyi carries an increased number of genes despite absence of recent whole-genome duplications

Heracleum sosnowskyi, belonging to a group of giant hogweeds, is a plant with large effects on ecosystems and human health. It is an invasive species that contributes to the deterioration of grassland ecosystems. The ability of H. sosnowskyi to produce linear furanocoumarins (FCs), photosensitizing compounds, makes it very dangerous. At the same time, linear FCs are compounds with high pharmaceutical value used in skin disease therapies. Despite this high importance, it has not been the focus of genetic and genomic studies. Here, we report a chromosome-scale assembly of Sosnowsky's hogweed genome. Genomic analysis revealed an unusually high number of genes (55106) in the hogweed genome, in contrast to the 25-35 thousand found in most plants. However, we did not find any traces of recent whole-genome duplications not shared with its confamiliar, Daucus carota (carrot), which has approximately thirty thousand genes. The analysis of the genomic proximity of duplicated genes indicates on tandem duplications as a main reason for this increase. We performed a genome-wide search of the genes of the FC biosynthesis pathway and surveyed their expression in aboveground plant parts. Using a combination of expression data and phylogenetic analysis, we found candidate genes for psoralen synthase and experimentally showed the activity of one of them using a heterologous yeast expression system. These findings expand our knowledge on the evolution of gene space in plants and lay a foundation for further analysis of hogweed as an invasive plant and as a source of FCs.

Genome-wide mining and characterization of MATE transporters in Coriandrum sativum L

Multidrug and Toxic Compound Extrusion (MATE) proteins are responsible for the transport of a wide range of metabolites out of plant cells. This helps to protect the cells from toxins and other harmful compounds. MATE proteins also play a role in plant development, by regulating the transport of hormones and other signalling molecules. They transport a wide variety of substances, including organic acids, plant hormones, flavonoids, alkaloids, terpenes and other secondary metabolites. MATE proteins are thought to play similar roles in Coriander, in addition to stress responses. The MATE genes in the coriander genome have been identified and characterized. Detailed genome homology search and domain identification analysis have identified 91 MATE proteins in the genome assembly of coriander. A phylogenetic analysis of the identified proteins divided them into five major clades. The functions of the transporters in each cluster were predicted based on the clustering pattern of the functionally characterized proteins. The amino acid sequences, exon-intron structures and motif details of all the 91 proteins are identified and described. This is the first work on the MATE transporters in coriander and the results deliver clues for the molecular mechanisms behind the stress responses and secondary metabolite transport in coriander.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11 Mb, with an N50 scaffold of 46.20 Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40 Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62 Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

A comprehensive review on genomic resources in medicinally and industrially important major spices for future breeding programs: Status, utility and challenges

In the global market, spices possess a high-value but low-volume commodities of commerce. The food industry depends largely on spices for taste, flavor, and therapeutic properties in replacement of cheap synthetic ones. The estimated growth rate for spices demand in the world is ∼3.19%. Since spices grow in limited geographical regions, India is one of the leading producer of spices, contributing 25-30 percent of total world trade. Hitherto, there has been no comprehensive review of the genomic resources of industrially important major medicinal spices to overcome major impediments in varietal improvement and management. This review focuses on currently available genomic resources of 24 commercially significant spices, namely, Ajwain, Allspice, Asafoetida, Black pepper, Cardamom large, Cardamom small, Celery, Chillies, Cinnamon, Clove, Coriander, Cumin, Curry leaf, Dill seed, Fennel, Fenugreek, Garlic, Ginger, Mint, Nutmeg, Saffron, Tamarind, Turmeric and Vanilla. The advent of low-cost sequencing machines has contributed immensely to the voluminous data generation of these spices, cracking the complex genomic architecture, marker discovery, and understanding comparative and functional genomics. This review of spice genomics resources concludes the perspective and way forward to provide footprints by uncovering genome assemblies, sequencing and re-sequencing projects, transcriptome-based studies, non-coding RNA-mediated regulation, organelles-based resources, developed molecular markers, web resources, databases and AI-directed resources in candidate spices for enhanced breeding potential in them. Further, their integration with molecular breeding could be of immense use in formulating a strategy to protect and expand the production of the spices due to increased global demand.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Comparative-genomic analysis reveals dynamic NLR gene loss and gain across Apiaceae species

Introduction: Nucleotide-binding leucine-rich repeat (NLR) genes play a crucial role in green plants' responding to various pathogens. Genome-scale evolutionary studies of NLR genes are important for discovering and applying functional NLR genes. However, little is known about the evolution of NLR genes in the Apiaceae family including agricultural and medical plants. Methods: In this study, comparative genomic analysis was performed in four Apiaceae species to trace the dynamic evolutionary patterns of NLR genes during speciation in this family. Results: The results revealed different number of NLR genes in these four Apiaceae species, namely, Angelica sinensis (95), Coriandrum sativum (183), Apium graveolens (153) and Daucus carota (149). Phylogenetic analysis demonstrated that NLR genes in these four species were derived from 183 ancestral NLR lineages and experienced different levels of gene-loss and gain events. The contraction pattern of the ancestral NLR lineages was discovered during the evolution of D. carota, whereas a different pattern of contraction after first expansion of NLR genes was observed for A. sinensis, C. sativum and A. graveolens. Discussion: Taken together, rapid and dynamic gene content variation has shaped evolutionary history of NLR genes in Apiaceae species.

PlantNLRatlas: a comprehensive dataset of full- and partial-length NLR resistance genes across 100 chromosome-level plant genomes

Plants have evolved two layers of protection against biotic stress: PAMP-triggered immunity (PTI) and effector-triggered immunity (ETI). The primary mechanism of ETI involves nucleotide-binding leucine-rich repeat immune receptors (NLRs). Although NLR genes have been studied in several plant species, a comprehensive database of NLRs across a diverse array of species is still lacking. Here, we present a thorough analysis of NLR genes across 100 high-quality plant genomes (PlantNLRatlas). The PlantNLRatlas includes a total of 68,452 NLRs, of which 3,689 are full-length and 64,763 are partial-length NLRs. The majority of NLR groups were phyletically clustered. In addition, the domain sequences were found to be highly conserved within each NLR group. Our PlantNLRatlas dataset is complementary to RefPlantNLR, a collection of NLR genes which have been experimentally confirmed. The PlantNLRatlas should prove helpful for comparative investigations of NLRs across a range of plant groups, including understudied taxa. Finally, the PlantNLRatlas resource is intended to help the field move past a monolithic understanding of NLR structure and function.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

The chromosome-level genome of female ginseng (Angelica sinensis) provides insights into molecular mechanisms and evolution of coumarin biosynthesis

Coumarins are natural products with important medicinal values, and include simple coumarins, furanocoumarins and pyranocoumarins. Female ginseng (Angelica sinensis) is a renowned herb with abundant coumarins, originated in China and known for the treatment of female ailments for thousands of years. The molecular basis of simple coumarin biosynthesis in A. sinensis and the evolutionary history of the genes involved in furanocoumarin biosynthesis are largely unknown. Here, we generated the first chromosome-scale genome of A. sinensis. It has a genome size of 2.37 Gb, which was generated by combining PacBio and Hi-C sequencing technologies. The genome was predicted to contain 43 202 protein-coding genes dispersed mainly on 11 pseudochromosomes. We not only provided evidence for whole-genome duplication (WGD) specifically occurring in the Apioideae subfamily, but also demonstrated the vital role of tandem duplication for phenylpropanoid biosynthesis in A. sinensis. Combined analyses of transcriptomic and metabolomic data revealed key genes and candidate transcription factors regulating simple coumarin biosynthesis. Furthermore, phylogenomic synteny network analyses suggested prenyltransferase genes involved in furanocoumarin biosynthesis evolved independently in the Moraceae, Fabaceae, Rutaceae and Apiaceae after ζ and ε WGD. Our work sheds light on coumarin biosynthesis, and provides a benchmark for accelerating genetic research and molecular breeding in A. sinensis.

Genome-wide identification and characterization profile of phosphatidy ethanolamine-binding protein family genes in carrot

Members of the family of Phosphatidy Ethanolamine-Binding Protein (PEBP) have been shown to be key regulators of the transition of plants from vegetative to reproductive phases. Here, a total of 12 PEBP proteins were identified in the carrot (Daucus carota L.) genome and classified into FT-like (4), TFL1-like (6), and MFT-like 2) subfamilies, that had different lengths (110-267 aa) and were distributed unevenly across seven chromosomes. Moreover, 13 and 31 PEBP proteins were identified in other two Apiaceae species, celery (Apium graveolens L.) and coriander (Coriandrum sativum L.). The phylogenetic and evolutionary results of these PEBP family proteins were obtained based on the protein sequences. In the three Apiaceae species, purifying selection was the main evolutionary force, and WGD, segmental duplication, and dispersed duplication have played key roles in the PEBP family expansion. The expression analysis showed that carrot PEBP genes exhibited relatively broad expression patterns across various tissues. In the period of bolting to flowering, the carrot FT-like subfamily genes were upregulated as positive regulators, and TFL1-like subfamily genes remained at lower expression levels as inhibitors. More interestingly, the members of carrot FT-like genes had different temporal-spatial expression characteristics, suggesting that they have different regulatory functions in the carrot reproductive phase. In summary, this study contributes to our understanding of the PEBP family proteins and provides a foundation for exploring the mechanism of carrot bolting and flowering for the breeding of cultivars with bolting resistance.

Genome-wide identification and analysis of terpene synthase (TPS) genes in celery reveals their regulatory roles in terpenoid biosynthesis

Terpenes are an important class of secondary metabolites in celery, which determine its flavor. Terpene synthase (TPS) has been established as a key enzyme in the biosynthesis of terpenes. This study systematically analyzed all members of the TPS gene family of celery (Apium graveolens) based on whole genome data. A total of 39 celery TPS genes were identified, among which TPS-a and TPS-b represented the two largest subfamilies. 77 cis-element types were screened in the promoter regions of AgTPS genes, suggesting the functional diversity of members of this family. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that AgTPS genes were enriched in multiple terpenoid biosynthesis pathways. Transcript abundance analysis and qRT-PCR showed that most AgTPS genes were differentially expressed in different tissues and colors of celery, with AgTPS 6, 9, and 11 expressed differentially in tissues, while AgTPS31, 32, and 38 are expressed differently in colors. More than 70% of the celery volatile compounds identified by HS-SPME-GC/MS were terpene, and the most critical compounds were β-Myrcene, D-Limonene, β-Ocimene and γ-Terpinene. Principal component analysis (PCA) showed that compounds (E)-β-Ocimene, D-Limonene, β-Myrcene and γ-Terpinene predominantly accounted for the variation. Further correlation analysis between gene expression and terpenoid accumulation showed that the four genes AgTPS9, 25, 31 and 38 genes may have positive regulatory effects on the synthesis of D-Limonene and β-Myrcene in celery. Overall, this study identified key candidate genes that regulate the biosynthesis of volatile compounds and provide the foothold for the development and utilization of terpenoids in celery.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Differential hydroxylation efficiency of the two non-heme carotene hydroxylases: DcBCH1, rather than DcBCH2, plays a major role in carrot taproot

Carotene hydroxylase plays an important role in catalyzing the hydroxylation of carotene to xanthopylls, including two types: non-heme carotene hydroxylase (BCH type) and heme-containing cytochrome P450 hydroxylase (P450 type). Two BCH-encoding genes were annotated in the carrot genome. However, the role of BCHs and whether there are functional interactions between the duplicated BCHs in carrot remains unclear. In this study, two BCH encoding genes, DcBCH1 and DcBCH2, were cloned from carrot. The relative expression level of DcBCH1 was much higher than that of DcBCH2 in carrot taproots with different carotene accumulation levels. Overexpression of DcBCH1 in 'KRD' (high carotene accumulated) carrot changed the taproot color from orange to yellow, accompanied by substantial reductions in α-carotene and β-carotene. There was no obvious change in taproot color between transgenic 'KRD' carrot overexpressing DcBCH2 and control carrot. Simultaneously, the content of α-carotene in the taproot of DcBCH2-overexpressing carrot decreased, but the content of β-carotene did not change significantly in comparison with control carrot. Using the CRISPR/Cas9 system to knock out DcBCH1 in 'KRD' carrot lightened the taproot color from orange to pink-orange; the content of α-carotene in the taproot increased slightly, while the β-carotene content was still significantly decreased, compared with control carrot. In DcBCH1-knockout carrot, the transcript level of DcBCH2 was significantly increased. These results indicated that in carrot taproot, DcBCH1 played the main function of BCH enzyme, which could hydroxylate α-carotene and β-carotene; DcBCH1 and DcBCH2 had functional redundancy, and these two DcBCHs could partially compensate for each other.

GGDB: A Grameneae genome alignment database of homologous genes hierarchically related to evolutionary events

The genomes of Gramineae plants have been preferentially sequenced owing to their economic value. These genomes are often quite complex, for example harboring many duplicated genes, and are the main source of genetic innovation and often the result of recurrent polyploidization. Deciphering these complex genome structures and linking duplicated genes to specific polyploidization events are important for understanding the biology and evolution of plants. However, efforts have been hampered by the complexity of analyzing these genomes. Here, we analyzed 29 well-assembled and up-to-date Gramineae genome sequences by hierarchically relating duplicated genes in collinear regions to specific polyploidization or speciation events. We separated duplicated genes produced by each event, established lists of paralogous and orthologous genes, and ultimately constructed an online database, GGDB (http://www.grassgenome.com/). Homologous gene lists from each plant and between plants can be displayed, searched, and downloaded from the database. Interactive comparison tools are deployed to demonstrate homology among user-selected plants and to draw genome-scale or local alignment figures and gene-based phylogenetic trees corrected by exploiting gene collinearity. Using these tools and figures, users can easily detect structural changes in genomes and explore the effects of paleo-polyploidy on crop genome structure and function. The GGDB will provide a useful platform for improving our understanding of genome changes and functional innovation in Gramineae plants.

Transcriptomics and Genomics Analysis Uncover the Differentially Expressed Chlorophyll and Carotenoid-Related Genes in Celery

Celery (Apium graveolens L.), a plant from Apiaceae, is one of the most important vegetables and is grown worldwide. Carotenoids can capture light energy and transfer it to chlorophyll, which plays a central role in photosynthesis. Here, by performing transcriptomics and genomics analysis, we identified and conducted a comprehensive analysis of chlorophyll and carotenoid-related genes in celery and six representative species. Significantly, different contents and gene expression patterns were found among three celery varieties. In total, 237 and 290 chlorophyll and carotenoid-related genes were identified in seven species. No notable gene expansion of chlorophyll biosynthesis was detected in examined species. However, the gene encoding ζ-carotene desaturase (ZDS) enzyme in carotenoid was expanded in celery. Comparative genomics and RNA-seq analyses revealed 16 and 5 key genes, respectively, regulating chlorophyll and carotenoid. An intriguing finding is that chlorophyll and carotenoid-related genes were coordinately regulated by transcriptional factors, which could be distinctively classified into positive- and negative-regulation groups. Six CONSTANS (CO)-like transcription factors co-regulated chlorophyll and carotenoid-related genes were identified in celery. In conclusion, this study provides new insights into the regulation of chlorophyll and carotenoid by transcription factors.

Origin, evolution, breeding, and omics of Apiaceae: a family of vegetables and medicinal plants

Many of the world's most important vegetables and medicinal crops, including carrot, celery, coriander, fennel, and cumin, belong to the Apiaceae family. In this review, we summarize the complex origins of Apiaceae and the current state of research on the family, including traditional and molecular breeding practices, bioactive compounds, medicinal applications, nanotechnology, and omics research. Numerous molecular markers, regulatory factors, and functional genes have been discovered, studied, and applied to improve vegetable and medicinal crops in Apiaceae. In addition, current trends in Apiaceae application and research are also briefly described, including mining new functional genes and metabolites using omics research, identifying new genetic variants associated with important agronomic traits by population genetics analysis and GWAS, applying genetic transformation, the CRISPR-Cas9 gene editing system, and nanotechnology. This review provides a reference for basic and applied research on Apiaceae vegetable and medicinal plants.

Recent applications of metabolomics in plant breeding

Metabolites play a central role in maintaining organismal life and in defining crop phenotypes, such as nutritional value, fragrance, color, and stress resistance. Among the 'omes' in biology, the metabolome is the closest to the phenotype. Consequently, metabolomics has been applied to crop improvement as method for monitoring changes in chemical compositions, clarifying the mechanisms underlying cellular functions, discovering markers and diagnostics, and phenotyping for mQTL, mGWAS, and metabolite-genome predictions. In this review, 359 reports of the most recent applications of metabolomics to plant breeding-related studies were examined. In addition to the major crops, more than 160 other crops including rare medicinal plants were considered. One bottleneck associated with using metabolomics is the wide array of instruments that are used to obtain data and the ambiguity associated with metabolite identification and quantification. To further the application of metabolomics to plant breeding, the features and perspectives of the technology are discussed.

Agri-biotechnology of coriander (Coriandrum sativum L.): an inclusive appraisal

Bountiful expression of bioactivity of phytochemicals obtained from spice crops like coriander gifts them the label of being natural antioxidants. It is well-accepted and time-tested towards contributing to human wellbeing. The accomplishment of coriander production is fundamentally influenced by genetic, agroclimatic, and agronomic factors. Despite the fact that there are very restricted options to manage the first two factors, the third one is apparently imperative to arbitrate as far as the elevated yield and enhanced quality are concerned. On the other hand, an indomitable, object-oriented, controlled agrotechnological and biotechnological intervention can also contribute towards better yield and quality of coriander. There are several accounts of the successful use of such technologies in order to genetically improve the qualitative and quantitative indicators of coriander. However, often these areas are not comprehensively explored and utilized. In that context, the present review highlights the botanical features, origin and distribution, multi-dimensional importance, pre- and post-harvest crop management, phytochemical production, and germplasm conservation, including the in vitro-based regeneration methods along with molecular marker-based biotechnological and omics approaches attempted in coriander until date. In addition, the possibility of the yet-to-be-explored agri-biotechnological methods and their potential for genetic improvement of this crop has also been reviewed in this appraisal. KEY POINTS: • Coriander, used both as an herb and spice, is popular in the pharmaceutical and culinary industries. • The current review provides insight into agrotechnological and biotechnological interventions for better yield and quality. • Provides novel ideas to harness the comprehensive qualitative and quantitative genetic improvement based on the potential use of promising biotechnological tools and techniques.

A chromosome-scale genome assembly for the holly (Ilex polyneura) provides insights into genomic adaptations to elevation in Southwest China

Southwest China is a plant diversity hotspot. The near-cosmopolitan genus Ilex (c. 664 spp., Aquifoliaceae) reaches its maximum diversity in this region, with many narrow-range and a few widespread species. Divergent selection on widespread species leads to local adaptation, with consequences for both conservation and utilization, but is counteracted by geneflow. Many Ilex species are utilized as teas, medicines, ornamentals, honey plants, and timber, but variation below the species level is largely uninvestigated. We therefore studied the widespread Ilex polyneura, which occupies most of the elevational range available and is cultivated for its decorative leafless branches with persistent red fruits. We assembled a chromosome-scale genome using approximately 100x whole genome long-read and short-read sequencing combined with Hi-C sequencing. The genome is approximately 727.1 Mb, with a contig N50 size of 5 124 369 bp and a scaffold N50 size of 36 593 620 bp, for which the BUSCO score was 97.6%, and 98.9% of the assembly was anchored to 20 pseudochromosomes. Out of 32 838 genes predicted, 96.9% were assigned functions. Two whole genome duplication events were identified. Using this genome as a reference, we conducted a population genomics study of 112 individuals from 21 populations across the elevation range using restriction site-associated DNA sequencing (RADseq). Most populations clustered into four clades separated by distance and elevation. Selective sweep analyses identified 34 candidate genes potentially under selection at different elevations, with functions related to responses to abiotic and biotic stresses. This first high-quality genome in the Aquifoliales will facilitate the further domestication of the genus.

Chromosome-level genome assembly of Amomum tsao-ko provides insights into the biosynthesis of flavor compounds

Amomum tsao-ko is an economically important spice plant in the ginger family (Zingiberaceae). The dried ripe fruit has been widely used as spice and medicine in Southeast Asia due to its distinct flavor metabolites. However, there is little genomic information available to understand the biosynthesis of its characteristic flavor compounds. Here, we present a high-quality chromosome-level genome of A. tsao-ko with a total length of 2.08 Gb assembled into 24 chromosomes. Potential relationships between genetic variation and chemical constituents were analyzed by a genome-wide association study of 119 representative A. tsao-ko specimens in China. Metabolome and transcriptome correlation analysis of different plant organs and fruit developmental stages revealed the proposed biosynthesis of the characteristic bicyclononane aldehydes and aromatic metabolites in A. tsao-ko fruit. Transcription factors of 20 families may be involved in the regulatory network of terpenoids. This study provides genomic and chemical insights into the biosynthesis of characteristic aroma and flavor constituents, which can be used to improve the quality of A. tsao-ko as food and medicine.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

Integrative genome, transcriptome, microRNA, and degradome analysis of water dropwort (Oenanthe javanica) in response to water stress

Water dropwort (Liyang Baiqin, Oenanthe javanica (BI.) DC.) is an aquatic perennial plant from the Apiaceae family with abundant protein, dietary fiber, vitamins, and minerals. It usually grows in wet soils and can even grow in water. Here, whole-genome sequencing of O. javanica via HiSeq 2000 sequencing technology was reported for the first time. The genome size was 1.28 Gb, including 42,270 genes, of which 93.92% could be functionally annotated. An online database of the whole-genome sequences of water dropwort, Water dropwortDB, was established to share the results and facilitate further research on O. javanica (database homepage: http://apiaceae.njau.edu.cn/waterdropwortdb ). Water dropwortDB offers whole-genome and transcriptome sequences and a Basic Local Alignment Search Tool. Comparative analysis with other species showed that the evolutionary relationship between O. javanica and Daucus carota was the closest. Twenty-five gene families of O. javanica were found to be expanded, and some genetic factors (such as genes and miRNAs) related to phenotypic and anatomic differentiation in O. javanica under different water conditions were further investigated. Two miRNA and target gene pairs (miR408 and Oja15472, miR171 and Oja47040) were remarkably regulated by water stress. The obtained reference genome of O. javanica provides important information for future work, thus making in-depth genetic breeding and gene editing possible. The present study also provides a foundation for the understanding of the O. javanica response to water stress, including morphological, anatomical, and genetic differentiation.

Impact of Genomic and Transcriptomic Resources on Apiaceae Crop Breeding Strategies

The Apiaceae taxon is one of the most important families of flowering plants and includes thousands of species used for food, flavoring, fragrance, medical and industrial purposes. This study had the specific intent of reviewing the main genomics and transcriptomic data available for this family and their use for the constitution of new varieties. This was achieved starting from the description of the main reproductive systems and barriers, with particular reference to cytoplasmic (CMS) and nuclear (NMS) male sterility. We found that CMS and NMS systems have been discovered and successfully exploited for the development of varieties only in Foeniculum vulgare, Daucus carota, Apium graveolens and Pastinaca sativa; whereas, strategies to limit self-pollination have been poorly considered. Since the constitution of new varieties benefits from the synergistic use of marker-assisted breeding in combination with conventional breeding schemes, we also analyzed and discussed the available SNP and SSR marker datasets (20 species) and genomes (8 species). Furthermore, the RNA-seq studies aimed at elucidating key pathways in stress tolerance or biosynthesis of the metabolites of interest were limited and proportional to the economic weight of each species. Finally, by aligning 53 plastid genomes from as many species as possible, we demonstrated the precision offered by the super barcoding approach to reconstruct the phylogenetic relationships of Apiaceae species. Overall, despite the impressive size of this family, we documented an evident lack of molecular data, especially because genomic and transcriptomic resources are circumscribed to a small number of species. We believe that our contribution can help future studies aimed at developing molecular tools for boosting breeding programs in crop plants of the Apiaceae family.

Haplotype-resolved genome of diploid ginger (Zingiber officinale) and its unique gingerol biosynthetic pathway

Ginger (Zingiber officinale), the type species of Zingiberaceae, is one of the most widespread medicinal plants and spices. Here, we report a high-quality, chromosome-scale reference genome of ginger 'Zhugen', a traditionally cultivated ginger in Southwest China used as a fresh vegetable, assembled from PacBio long reads, Illumina short reads, and high-throughput chromosome conformation capture (Hi-C) reads. The ginger genome was phased into two haplotypes, haplotype 1 (1.53 Gb with a contig N50 of 4.68 M) and haplotype 0 (1.51 Gb with a contig N50 of 5.28 M). Homologous ginger chromosomes maintained excellent gene pair collinearity. In 17,226 pairs of allelic genes, 11.9% exhibited differential expression between alleles. Based on the results of ginger genome sequencing, transcriptome analysis, and metabolomic analysis, we proposed a backbone biosynthetic pathway of gingerol analogs, which consists of 12 enzymatic gene families, PAL, C4H, 4CL, CST, C3'H, C3OMT, CCOMT, CSE, PKS, AOR, DHN, and DHT. These analyses also identified the likely transcription factor networks that regulate the synthesis of gingerol analogs. Overall, this study serves as an excellent resource for further research on ginger biology and breeding, lays a foundation for a better understanding of ginger evolution, and presents an intact biosynthetic pathway for species-specific gingerol biosynthesis.

Development and characterization of genic SSR-FDM for stem gall disease resistance in coriander (Coriandrum sativum L.) and its cross species transferability

Coriander (Coriandrum sativum L.) is well known vegetable and spice crop grown globally for its leaves and seeds. Stem gall (Protomyces macrosporus L.) is a fungal disease affecting its quality and yield. However, no information is available on SSR markers linked to disease resistance in coriander. Hence, development of co-dominant genetic markers is prerequisite for disease investigations in coriander. In-house stem gall resistance and susceptible cultivars transcriptome data were utilized. Totally, 59,933 and 56,861 transcripts were examined, 9141 and 8346 Simple Sequence Repeats (SSR) were identified and the most abundant type was the tri, followed by di, tetra, penta and hexa nucleotide repeats. A total of ten selected SSR-Functional Domain Markers (FDM) were developed based on functional annotation terms associated with pathogen response and validated among ten coriander cultivars and their transferability was examined in five fennel (Foeniculum vulgare L.) cultivars. Nine primer pairs resulted from amplified bands. Marker ACorSGD-1 shown monomorphic bands among coriander genotypes except Acr-1 showed heterologouse and multiple bands in fennel cultivars. Markers ACorSGD-4, 5, 7 and 9 shown presence in resistant cultivars and absence of bands among susceptible cultivars of coriander and thus, considered to be the candidate markers for disease screening. Marker ACorSGD-6 shown monomorphic bands among coriander. Markers ACorSGD-1, 2, 3, and 5 shown transferability among fennel cultivars. A total of 136 alleles in coriander and fennel were produced. Using UPGMA clustering method a dendrogram was generated and cultivars were grouped into two separate clusters with coriander and fennel. Identified and developed SSR-FDM markers are useful for linkage mapping for disease resistant in coriander.

De novo chromosome-level assembly of the Centella asiatica genome

Centella asiatica is a herbaceous, perennial species indigenous to India and Southeast Asia. C. asiatica possesses several medicinal properties: anti-aging, anti-inflammatory, wound healing and memory enhancing. The lack of available genomics resources significantly impedes the improvement of C. asiatica varieties through molecular breeding. Here, we combined the 10× Genomics linked-read technology and the long-range HiC technique to obtain the genome assembly. The final assembly contained nine pseudomolecules, corresponding to the haploid chromosome number in C. asiatica. These nine chromosomes covered 402,536,584 bases or 93.6% of the 430-Mb assembly. Comparative genomics analyses based on single-copy orthologous genes showed that C. asiatica and the common ancestor of Coriandrum sativum (coriander) and Daucus carota (carrot) diverged about 48 million years ago. This assembly provides a valuable reference genome for future molecular studies, varietal development through marker-assisted breeding and comparative genomics studies in C. asiatica.

The celery genome sequence reveals sequential paleo-polyploidizations, karyotype evolution and resistance gene reduction in apiales

Celery (Apium graveolens L. 2n = 2x = 22), a member of the Apiaceae family, is among the most important and globally grown vegetables. Here, we report a high-quality genome sequence assembly, anchored to 11 chromosomes, with total length of 3.33 Gb and N50 scaffold length of 289.78 Mb. Most (92.91%) of the genome is composed of repetitive sequences, with 62.12% of 31 326 annotated genes confined to the terminal 20% of chromosomes. Simultaneous bursts of shared long-terminal repeats (LTRs) in different Apiaceae plants suggest inter-specific exchanges. Two ancestral polyploidizations were inferred, one shared by Apiales taxa and the other confined to Apiaceae. We reconstructed 8 Apiales proto-chromosomes, inferring their evolutionary trajectories from the eudicot common ancestor to extant plants. Transcriptome sequencing in three tissues (roots, leaves and petioles), and varieties with different-coloured petioles, revealed 4 and 2 key genes in pathways regulating anthocyanin and coumarin biosynthesis, respectively. A remarkable paucity of NBS disease-resistant genes in celery (62) and other Apiales was explained by extensive loss and limited production of these genes during the last ~10 million years, raising questions about their biotic defence mechanisms and motivating research into effects of chemicals, for example coumarins, that give off distinctive odours. Celery genome sequencing and annotation facilitates further research into important gene functions and breeding, and comparative genomic analyses in Apiales.

Genome-Wide Identification and Comparative Analysis of ARF Family Genes in Three Apiaceae Species

The family Apiaceae includes many important vegetables and medicinal plants. Auxin response factors (ARFs) play critical roles in regulating plant growth and development. Here, we performed a comprehensive analysis of the ARF gene family in three Apiaceae species, celery, coriander, and carrot, and compared the results with the ARF gene family of lettuce, Arabidopsis, and grape. We identified 156 ARF genes in all six species and 89 genes in the three Apiaceae species, including 28, 34, and 27 in celery, coriander, and carrot, respectively. The paralogous gene number in coriander was far greater than that in carrot and celery. Our analysis revealed that ARF genes of the three Apiaceae species in 34 branches of the phylogenetic tree underwent significant positive selection. Additionally, our findings indicated that whole-genome duplication played an important role in ARF gene family expansion. Coriander contained a greater number of ARF genes than celery and carrot because of more gene duplications and less gene losses. We also analyzed the expression of ARF genes in three tissues by RNA-seq and verified the results by quantitative real-time PCR. Furthermore, we found that several paralogous genes exhibited divergent expression patterns. Overall, this study provides a valuable resource for exploring how ARF family genes regulate plant growth and development in other plants. Since this is the first report of the ARF gene family in Apiaceae, our results will serve as a guide for comparative and functional analyses of ARF and other gene families in Apiaceae.

Chromosome-scale genome assembly provides insights into the evolution and flavor synthesis of passion fruit (Passiflora edulis Sims)

Passion fruit (Passiflora edulis Sims) is an economically valuable fruit that is cultivated in tropical and subtropical regions of the world. Here, we report an ~1341.7 Mb chromosome-scale genome assembly of passion fruit, with 98.91% (~1327.18 Mb) of the assembly assigned to nine pseudochromosomes. The genome includes 23,171 protein-coding genes, and most of the assembled sequences are repetitive sequences, with long-terminal repeats (LTRs) being the most abundant. Phylogenetic analysis revealed that passion fruit diverged after Brassicaceae and before Euphorbiaceae. Ks analysis showed that two whole-genome duplication events occurred in passion fruit at 65 MYA and 12 MYA, which may have contributed to its large genome size. An integrated analysis of genomic, transcriptomic, and metabolomic data showed that 'alpha-linolenic acid metabolism', 'metabolic pathways', and 'secondary metabolic pathways' were the main pathways involved in the synthesis of important volatile organic compounds (VOCs) in passion fruit, and this analysis identified some candidate genes, including GDP-fucose Transporter 1-like, Tetratricopeptide repeat protein 33, protein NETWORKED 4B isoform X1, and Golgin Subfamily A member 6-like protein 22. In addition, we identified 13 important gene families in fatty acid pathways and eight important gene families in terpene pathways. Gene family analysis showed that the ACX, ADH, ALDH, and HPL gene families, especially ACX13/14/15/20, ADH13/26/33, ALDH1/4/21, and HPL4/6, were the key genes for ester synthesis, while the TPS gene family, especially PeTPS2/3/4/24, was the key gene family for terpene synthesis. This work provides insights into genome evolution and flavor trait biology and offers valuable resources for the improved cultivation of passion fruit.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

On genetic diversity in caraway: Genotyping of a large germplasm collection

Caraway (Carum carvi) is a widespread and frequently used spice and medicinal plant with a long history of cultivation. However, due to ongoing climatic changes, the cultivation is becoming increasingly risky. To secure caraway cultivation in future, timely breeding efforts to develop adapted material are necessary. Analysis of genetic diversity can accompany this process, for instance, by revealing untapped gene pools. Here, we analyzed 137 accessions using genotyping by sequencing (GBS). Hence, we can report a broad overview of population structure and genetic diversity of caraway. Population structure was determined using a principal coordinate analysis, a Bayesian clustering analysis, phylogenetic trees and a neighbor network based on 13,155 SNPs. Genotypic data indicate a clear separation of accessions into two subpopulations, which correlates with the flowering type (annual vs. biennial). Four winter-annual accessions were closer related to biennial accessions. In an analysis of molecular variance, genetic variation between the two subpopulations was 7.84%. In addition, we estimated the genome size for 35 accessions by flow cytometry. An average genome size of 4.282 pg/2C (± 0.0096 S.E.) was estimated. Therefore, we suggest a significantly smaller genome size than stated in literature.

Identification and expression analysis of candidate genes associated with stem gall disease in Coriander (Coriandrum sativum L.) cultivars

Coriander (Coriandrum sativum L.) is a well-known spice and aromatic crop cultivated globally. Stem gall disease is one of the major constraints for its leaf and seed quality used for consumption and also affecting the yield. The identification of resistance genes and further characterization of such genes could help to understand the molecular basis of resistance and lay a solid ground for cloning of stem gall resistance genes in coriander. To evaluate the genetic expression of disease resistance-relevant genes in popularly grown coriander cultivars in India such as Pant Haritma, Hisar Sugandh, Hisar Surabhi, Hisar Anand, Rajendra Swathi, ACr-1, ACr-2, AgCr-1, CO-2 and CS-6 were used for LRR, GDSL, USP, ANK and PDR gene expression using Real Time PCR along with 18S housekeeping gene as internal control for the normalization. Result revealed the different expression pattern of genes among the cultivars tested. Highest expression was shown in cultivar AgCr-1 followed by Pant Haritma, Hisar Sugandh and ACr-1, and least expression in Hisar Anand, ACr-2, CO-2, Rajendra Swathi and CS-6. Domain analysis revealed the conserved domain relevance of the genes. This is the first report on stem gall resistance gene expression in coriander. The identified genes have a potential role in coriander and further utilize in crop improvement program. We hypothesize that contrasting cultivars can be a good source for candidate gene evaluation and further to use them as potential markers and used in hybridization program focus on incorporating and develop durable disease-resistance into the adapted cultivars of the region.

Transcriptome Analysis Reveals Candidate Genes for Petroselinic Acid Biosynthesis in Fruits of Coriandrum sativum L

Petroselinic acid (18:1Δ⁶), a monounsaturated cis Δ-6 fatty acid, has many prospective applications in functional foods and for the nutraceutical and pharmaceutical industries. Up to 80% of petroselinic acid has been found in the oil from fruits of coriander (Coriandrum sativum L.), which make it an ideal source for investigating the biosynthesis of petroselinic acid. A coriander acyl-acyl carrier protein desaturase was identified to be involved in its biosynthesis more than two decades ago, but since then little further progress in this area has been reported. In this study, the fatty acid profiles of coriander fruits at six developmental stages were analyzed. Fruit samples from three developmental stages with rapid accumulation of petroselinic acid were used for RNA sequencing using the Illumina Hiseq4000 platform. The transcriptome analysis presented 93 323 nonredundant unigenes and 8545 differentially expressed genes. Functional annotation and combined gene expression data revealed candidate genes potentially involved in petroselinic acid biosynthesis and its incorporation into triacylglycerols. Tissue-specific examination of q-PCR validation further suggested that ACPD1/3, KAS I-1, FATB-1/3, and DGAT2 may be highly involved. Bioinformatic analysis of CsFATB and CsDGAT2 identified their putative key amino acids or functional motifs. These results provide a molecular foundation for petroselinic acid biosynthesis in coriander fruit and facilitate its genetic engineering in other hosts.