A Chromosome-Level Genome Assembly of Garlic (Allium sativum) Provides Insights into Genome Evolution and Allicin Biosynthesis

Incomplete lineage sorting and gene flow within Allium (Amayllidaceae)

The phylogeny and systematics of the genus Allium have been studied with a variety of diverse data types, including an increasing amount of molecular data. However, strong phylogenetic discordance and high levels of uncertainty have prevented the identification of a consistent phylogeny. The difficulty in establishing phylogenetic consensus and evidence for genealogical discordance make Allium a compelling test case to assess the relative contribution of incomplete lineage sorting (ILS), gene flow and gene tree estimation error on phylogenetic reconstruction. In this study, we obtained 75 transcriptomes of 38 Allium species across 10 subgenera. Whole plastid genome, single copy genes and consensus CDS were generated to estimate phylogenetic trees both using coalescence and concatenation methods. Multiple approaches including coalescence simulation, quartet sampling, reticulate network inference, sequence simulation, theta of ILS and reticulation index were carried out across the CDS gene trees to investigate the degrees of ILS, gene flow and gene tree estimation error. Afterward, a regression analysis was used to test the relative contributions of each of these forms of uncertainty to the final phylogeny. Despite extensive topological discordance among gene trees, we found a fully supported species tree that agrees with the most of well-accepted relationships and establishes monophyly of the genus Allium. We presented clear evidence for substantial ILS across the phylogeny of Allium. Further, we identified two ancient hybridization events for the formation of the second evolutionary line and subg. Butomissa as well as several introgression events between recently diverged species. Our regression analysis revealed that gene tree inference error and gene flow were the two most dominant factors explaining for the overall gene tree variation, with the difficulty in disentangling the effects of ILS and gene tree estimation error due to a positive correlation between them. Based on our efforts to mitigate the methodological errors in reconstructing trees, we believed ILS and gene flow are two principal reasons for the oft-reported phylogenetic heterogeneity of Allium. This study presents a strongly-supported and well-resolved phylogenetic backbone for the sampled Allium species, and exemplifies how to untangle heterogeneity in phylogenetic signal and reconstruct the true evolutionary history of the target taxa.

Identification of male-fertility gene AsaNRF1 and molecular marker development in cultivated garlic (Allium sativum L.)

Garlic cultivars are predominantly characterized by their sterility and reliance on asexual reproduction, which have traditionally prevented the use of hybrid breeding for cultivar improvement in garlic. Our investigation has revealed a notable exception in the garlic line G398, which demonstrates the ability to produce fertile pollen. Notably, at the seventh stage of anther development, callose degradation in the sterile line G390 was impeded, while G398 exhibited normal callose degradation. Transcriptome profiling revealed an enhanced expression of the callose-degrading gene, AsaNRF1, in the mature flower buds of the fertile line G398 compared to the sterile line G390. An insertion in the promoter of AsaNRF1 in G390 was identified, which led to its reduced expression at the tetrad stage and consequently delayed callose degradation, potentially resulting in the male sterility of G390. A discriminatory marker was developed to distinguish between fertile G398 and sterile G390, facilitating the assessment of male fertility in garlic germplasm resources. This study introduces a practical approach to harnessing garlic hybridization, which can further facilitate the breeding of new cultivars and the creation of novel male-fertile garlic germplasm using modern molecular biology methods.

Genome-wide identification, classification and expression analysis of the heat shock transcription factor family in Garlic (Allium sativum L.)

BACKGROUND

The heat shock transcription factor (HSF) plays a crucial role in the regulatory network by coordinating responses to heat stress as well as other stress signaling pathways. Despite extensive studies on HSF functions in various plant species, our understanding of this gene family in garlic, an important crop with nutritional and medicinal value, remains limited. In this study, we conducted a comprehensive investigation of the entire garlic genome to elucidate the characteristics of the AsHSF gene family.

RESULTS

In this study, we identified a total of 17 AsHSF transcription factors. Phylogenetic analysis classified these transcription factors into three subfamilies: Class A (9 members), Class B (6 members), and Class C (2 members). Each subfamily was characterized by shared gene structures and conserved motifs. The evolutionary features of the AsHSF genes were investigated through a comprehensive analysis of chromosome location, conserved protein motifs, and gene duplication events. These findings suggested that the evolution of AsHSF genes is likely driven by both tandem and segmental duplication events. Moreover, the nucleotide diversity of the AsHSF genes decreased by only 0.0002% from wild garlic to local garlic, indicating a slight genetic bottleneck experienced by this gene family during domestication. Furthermore, the analysis of cis-acting elements in the promoters of AsHSF genes indicated their crucial roles in plant growth, development, and stress responses. qRT-PCR analysis, co-expression analysis, and protein interaction prediction collectively highlighted the significance of Asa6G04911. Subsequent experimental investigations using yeast two-hybridization and yeast induction experiments confirmed its interaction with HSP70/90, reinforcing its significance in heat stress.

CONCLUSIONS

This study is the first to unravel and analyze the AsHSF genes in garlic, thereby opening up new avenues for understanding their functions. The insights gained from this research provide a valuable resource for future investigations, particularly in the functional analysis of AsHSF genes.

Nano-CeO2 activates physical and chemical defenses of garlic (Allium sativum L.) for reducing antibiotic resistance genes in plant endosphere

The transmission of manure- and wastewater-borne antibiotic-resistant bacteria (ARB) to plants contributes to the proliferation of antimicrobial resistance in agriculture, necessitating effective strategies for preventing the spread of antibiotic resistance genes (ARGs) from ARB in the environment to humans. Nanomaterials are potential candidates for efficiently controlling the dissemination of ARGs. The present study investigated the abundance of ARGs in hydroponically grown garlic (Allium sativum L.) following nano-CeO2 (nCeO2) application. Specifically, root exposure to nCeO2 (1, 2.5, 5, 10 mg L-1, 18 days) reduced ARG abundance in the endosphere of bulbs and leaves. The accumulation of ARGs (cat, tet, and aph(3')-Ia) in garlic bulbs decreased by 24.2-32.5 % after nCeO2 exposure at 10 mg L-1. Notably, the lignification extent of garlic stem-disc was enhanced by 10 mg L-1 nCeO2, thereby accelerating the formation of an apoplastic barrier to impede the upward transfer of ARG-harboring bacteria to garlic bulbs. Besides, nCeO2 upregulated the gene expression related to alliin biosynthesis and increased allicin content by 15.9-16.2 %, promoting a potent antimicrobial defense for reducing ARG-harboring bacteria. The potential exposure risks associated with ARGs and Ce were evaluated according to the estimated daily intake (EDI). The EDI of ARGs exhibited a decrease exceeding 95 %, while the EDI of Ce remained below the estimated oral reference dose. Consequently, through stimulating physical and chemical defenses, nCeO2 contributed to a reduced EDI of ARGs and Ce, highlighting its potential for controlling ARGs in plant endosphere within the framework of nano-enabled agrotechnology.

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae)

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

Crocus genome reveals the evolutionary origin of crocin biosynthesis

Crocus sativus (saffron) is a globally autumn-flowering plant, and its stigmas are the most expensive spice and valuable herb medicine. Crocus specialized metabolites, crocins, are biosynthesized in distant species, Gardenia (eudicot) and Crocus (monocot), and the evolution of crocin biosynthesis remains poorly understood. With the chromosome-level Crocus genome assembly, we revealed that two rounds of lineage-specific whole genome triplication occurred, contributing important roles in the production of carotenoids and apocarotenoids. According to the kingdom-wide identification, phylogenetic analysis, and functional assays of carotenoid cleavage dioxygenases (CCDs), we deduced that the duplication, site positive selection, and neofunctionalization of Crocus-specific CCD2 from CCD1 members are responsible for the crocin biosynthesis. In addition, site mutation of CsCCD2 revealed the key amino acids, including I143, L146, R161, E181, T259, and S292 related to the catalytic activity of zeaxanthin cleavage. Our study provides important insights into the origin and evolution of plant specialized metabolites, which are derived by duplication events of biosynthetic genes.

Genes encoding γ‑glutamyl‑transpeptidases in the allicin biosynthetic pathway in garlic (Allium sativum)

Allicin is a thiosulphate molecule produced in garlic (Allium sativum) and has a wide range of biological actions and pharmaceutical applications. Its precursor molecule is the non-proteinogenic amino acid alliin (S-allylcysteine sulphoxide). The alliin biosynthetic pathway in garlic involves a group of enzymes, members of which are the γ-glutamyl-transpeptidase isoenzymes, Allium sativum γ-glutamyl-transpeptidase AsGGT1, AsGGT2 and AsGGT3, which catalyze the removal of the γ-glutamyl group from γ-glutamyl-S-allyl-L-cysteine to produce S-allyl-L-cysteine. This removal is followed by an S-oxygenation, which leads to the biosynthesis of alliin. The aim of the present study is to annotate previously discovered genes of garlic γ-glutamyl-transpeptidases, as well as a fourth candidate gene (AsGGT4) that has yet not been described. The annotation includes identifying the loci of the genes in the garlic genome, revealing the overall structure and conserved regions of these genes, and elucidating the evolutionary history of these enzymes through their phylogenetic analysis. The genomic structure of γ-glutamyl-transpeptidase genes is conserved; each gene consists of seven exons, and these genes are located on different chromosomes. AsGGT3 and AsGGT4 enzymes contain a signal peptide. To that end, the AsGGT3 protein sequence was corrected; four indel events occurring in AsGGT3 coding regions suggested that at least in the garlic variety Ershuizao, AsGGT3 may be a pseudogene. Finally, the use of protein structure prediction tools allowed the visualization of the tertiary structure of the candidate peptide.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

AlliumDB: a central portal for comparative and functional genomics in Allium

The genus Allium belongs to the botanical family Amaryllidaceae and includes economically important crops such as onion, garlic, bunching onion, and leek, used as vegetables, spices, and traditional medicines. The large sizes of Allium genomes hamper the genetic dissection of agronomically important traits and molecular breeding. With the growing accumulation of genomic, resequencing, transcriptome, and phenotypic data, the demand for an integrative Allium database is increasing. Here we present a user-friendly database, AlliumDB (https://allium.qau.edu.cn), as a functional genomics hub integrating public and in-house data. The database contains all currently available nuclear and organelle genomes for Allium species, with genes comprehensively annotated based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, orthology, gene families, protein families (Pfam), and non-coding RNA families (Rfam). Transcriptome and variation profiles are integrated into dynamic visualization tools. We took phenotypic photographs and generated trait records for hundreds of Allium germplasms collected worldwide, which are included in the database. We incorporated JBrowse for the visualization of gene structures, RNA sequencing data, and variation data. Analysis tools such as the basic local alignment search tool (BLAST), sequence fetch, enrichment, and motif analyses are available to explore potential gene functions. This database incorporates comprehensive Allium genotypic and phenotypic datasets. As the community assembles new genomes and generates resequencing data for Allium germplasms, the database will be improved and continuously updated with these multi-omics data and comparative genomic studies. We expect the AlliumDB database to become a key resource for the study of Allium crops.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae)

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

High-quality assembly and methylome of a Tibetan wild tree peony genome (Paeonia ludlowii) reveal the evolution of giant genome architecture

Tree peony belongs to one of the Saxifragales families, Paeoniaceae. It is one of the most famous ornamental plants, and is also a promising woody oil plant. Although two Paeoniaceae genomes have been released, their assembly qualities are still to be improved. Additionally, more genomes from wild peonies are needed to accelerate genomic-assisted breeding. Here we assemble a high-quality and chromosome-scale 10.3-Gb genome of a wild Tibetan tree peony, Paeonia ludlowii, which features substantial sequence divergence, including around 75% specific sequences and gene-level differentials compared with other peony genomes. Our phylogenetic analyses suggest that Saxifragales and Vitales are sister taxa and, together with rosids, they are the sister taxon to asterids. The P. ludlowii genome is characterized by frequent chromosome reductions, centromere rearrangements, broadly distributed heterochromatin, and recent continuous bursts of transposable element (TE) movement in peony, although it lacks recent whole-genome duplication. These recent TE bursts appeared during the uplift and glacial period of the Qinghai-Tibet Plateau, perhaps contributing to adaptation to rapid climate changes. Further integrated analyses with methylome data revealed that genome expansion in peony might be dynamically affected by complex interactions among TE proliferation, TE removal, and DNA methylation silencing. Such interactions also impact numerous recently duplicated genes, particularly those related to oil biosynthesis and flower traits. This genome resource will not only provide the genomic basis for tree peony breeding but also shed light on the study of the evolution of huge genome structures as well as their protein-coding genes.

Deprivation of Sexual Reproduction during Garlic Domestication and Crop Evolution

Garlic, originating in the mountains of Central Asia, has undergone domestication and subsequent widespread introduction to diverse regions. Human selection for adaptation to various climates has resulted in the development of numerous garlic varieties, each characterized by specific morphological and physiological traits. However, this process has led to a loss of fertility and seed production in garlic crops. In this study, we conducted morpho-physiological and transcriptome analyses, along with whole-genome resequencing of 41 garlic accessions from different regions, in order to assess the variations in reproductive traits among garlic populations. Our findings indicate that the evolution of garlic crops was associated with mutations in genes related to vernalization and the circadian clock. The decline in sexual reproduction is not solely attributed to a few mutations in specific genes, but is correlated with extensive alterations in the genetic regulation of the annual cycle, stress adaptations, and environmental requirements. The regulation of flowering ability, stress response, and metabolism occurs at both the genetic and transcriptional levels. We conclude that the migration and evolution of garlic crops involve substantial and diverse changes across the entire genome landscape. The construction of a garlic pan-genome, encompassing genetic diversity from various garlic populations, will provide further insights for research into and the improvement of garlic crops.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

Genome resequencing reveals the evolutionary history of garlic reproduction traits

The propagation of cultivated garlic relies on vegetative cloves, thus flowers become non-essential for reproduction in this species, driving the evolution of reproductive feature-derived traits. To obtain insights into the evolutionary alteration of reproductive traits in the clonally propagated garlic, the evolutionary histories of two main reproduction-related traits, bolting and flower differentiation, were explored by genome analyses using 134 accessions displaying wide diversity in these two traits. Resequencing identified 272.8 million variations in the garlic genome, 198.0 million of which represent novel variants. Population analysis identified five garlic groups that have evolved into two clades. Gene expression, single-cell transcriptome sequencing, and genome-wide trait association analyses have identified numerous candidates that correlate with reproductive transition and flower development, some of which display distinct selection signatures. Selective forces acting on the B-box zinc finger protein-encoding Asa2G00291.1, the global transcription factor group E protein-encoding Asa5G01527.1, and VERNALIZATION INSENSITIVE 3-like Asa3G03399.1 appear to be representative of the evolution of garlic bolting. Plenty of novel genomic variations and trait-related candidates represent valuable resources for biological studies of garlic. Numerous selective signatures from genes associated with the two chosen reproductive traits provide important insights into the evolutionary history of reproduction in this clonally propagated crop.

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.

Host-specific adaptation drove the coevolution of leek yellow stripe virus and Allium plants

Host adaptation plays a crucial role in virus evolution and is a consequence of long-term interactions between virus and host in a complex arms race between host RNA silencing and viral RNA silencing suppressor (RSS) as counterdefense. Leek yellow stripe virus (LYSV), a potyvirus causing yield loss of garlic, infects several species of Allium plants. The unexpected discovery of an interspecific hybrid of garlic, leek, and great-headed (GH) garlic motivated us to explore the host-adaptive evolution of LYSV. Here, using Bayesian phylogenetic comparative methods and a functional assay of viral RSS activity, we show that the evolutionary context of LYSV has been shaped by the host adaptation of the virus during its coevolution with Allium plants. Our phylogenetic analysis revealed that LYSV isolates from leek and their taxonomic relatives (Allium ampeloprasum complex; AAC) formed a distinct monophyletic clade separate from garlic isolates and are likely to be uniquely adapted to AAC. Our comparative studies on viral accumulation indicated that LYSV accumulated at a low level in leek, whereas LYSVs were abundant in other Allium species such as garlic and its relatives. When RSS activity of the viral P1 and HC-Pro of leek LYSV isolate was analyzed, significant synergism in RSS activity between the two proteins was observed in leek but not in other species, suggesting that viral RSS activity may be important for the viral host-specific adaptation. We thus consider that LYSV may have undergone host-specific evolution at least in leek, which must be driven by speciation of its Allium hosts. IMPORTANCE Potyviruses are the most abundant plant RNA viruses and are extremely diversified in terms of their wide host range. Due to frequent host switching during their evolution, host-specific adaptation of potyviruses may have been shaped by numerous host factors. However, any critical determinants for viral host range remain largely unknown, possibly because of the repeated gain and loss of virus infectivity of plants. Leek yellow stripe virus (LYSV) is a species of the genus Potyvirus, which has a relatively narrow host range, generally limited to hosts in the genus Allium. Our investigations on leek and leek relatives (Allium ampeloprasum complex), which must have been generated through interspecies hybridization, revealed that LYSV accumulation remained low in leek as a result of viral host adaptation in competition with host resistance such as RNA silencing. This study presents LYSV as an ideal model to study the process of host-adaptive evolution and virus-host coevolution.

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.

Lineage-specific amplification and epigenetic regulation of LTR-retrotransposons contribute to the structure, evolution, and function of Fabaceae species

BACKGROUND

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) are ubiquitous and make up the majority of nearly all sequenced plant genomes, whereas their pivotal roles in genome evolution, gene expression regulation as well as their epigenetic regulation are still not well understood, especially in a large number of closely related species.

RESULTS

Here, we analyzed the abundance and dynamic evolution of LTR-RTs in 54 species from an economically and agronomically important family, Fabaceae, and also selected two representative species for further analysis in expression of associated genes, transcriptional activity and DNA methylation patterns of LTR-RTs. Annotation results revealed highly varied proportions of LTR-RTs in these genomes (5.1%~68.4%) and their correlation with genome size was highly positive, and they were significantly contributed to the variance in genome size through species-specific unique amplifications. Almost all of the intact LTR-RTs were inserted into the genomes 4 Mya (million years ago), and more than 50% of them were inserted in the last 0.5 million years, suggesting that recent amplifications of LTR-RTs were an important force driving genome evolution. In addition, expression levels of genes with intronic, promoter, and downstream LTR-RT insertions of Glycine max and Vigna radiata, two agronomically important crops in Fabaceae, showed that the LTR-RTs located in promoter or downstream regions suppressed associated gene expression. However, the LTR-RTs within introns promoted gene expression or had no contribution to gene expression. Additionally, shorter and younger LTR-RTs maintained higher mobility and transpositional potential. Compared with the transcriptionally silent LTR-RTs, the active elements showed significantly lower DNA methylation levels in all three contexts. The distributions of transcriptionally active and silent LTR-RT methylation varied across different lineages due to the position of LTR-RTs located or potentially epigenetic regulation.

CONCLUSION

Lineage-specific amplification patterns were observed and higher methylation level may repress the activity of LTR-RTs, further influence evolution in Fabaceae species. This study offers valuable clues into the evolution, function, transcriptional activity and epigenetic regulation of LTR-RTs in Fabaceae genomes.

DREB1 and DREB2 Genes in Garlic (Allium sativum L.): Genome-Wide Identification, Characterization, and Stress Response

Dehydration-responsive element-binding (DREB) transcription factors (TFs) of the A1 and A2 subfamilies involved in plant stress responses have not yet been reported in Allium species. In this study, we used bioinformatics and comparative transcriptomics to identify and characterize DREB A1 and A2 genes redundant in garlic (Allium sativum L.) and analyze their expression in A. sativum cultivars differing in the sensitivity to cold and Fusarium infection. Eight A1 (AsaDREB1.1-1.8) and eight A2 (AsaDREB2.1-2.8) genes were identified. AsaDREB1.1-1.8 genes located in tandem on chromosome 1 had similar expression patterns, suggesting functional redundancy. AsaDREB2.1-2.8 were scattered on different chromosomes and had organ- and genotype-specific expressions. AsaDREB1 and AsaDREB2 promoters contained 7 and 9 hormone- and stress-responsive cis-regulatory elements, respectively, and 13 sites associated with TF binding and plant development. In both Fusarium-resistant and -sensitive cultivars, fungal infection upregulated the AsaDREB1.1-1.5, 1.8, 2.2, 2.6, and 2.8 genes and downregulated AsaDREB2.5, but the magnitude of response depended on the infection susceptibility of the cultivar. Cold exposure strongly upregulated the AsaDREB1 genes, but downregulated most AsaDREB2 genes. Our results provide the foundation for further functional analysis of the DREB TFs in Allium crops and could contribute to the breeding of stress-tolerant varieties.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

The Sapria himalayana genome provides new insights into the lifestyle of endoparasitic plants

BACKGROUND

Sapria himalayana (Rafflesiaceae) is an endoparasitic plant characterized by a greatly reduced vegetative body and giant flowers; however, the mechanisms underlying its special lifestyle and greatly altered plant form remain unknown. To illustrate the evolution and adaptation of S. himalayasna, we report its de novo assembled genome and key insights into the molecular basis of its floral development, flowering time, fatty acid biosynthesis, and defense responses.

RESULTS

The genome of S. himalayana is ~ 1.92 Gb with 13,670 protein-coding genes, indicating remarkable gene loss (~ 54%), especially genes involved in photosynthesis, plant body, nutrients, and defense response. Genes specifying floral organ identity and controlling organ size were identified in S. himalayana and Rafflesia cantleyi, and showed analogous spatiotemporal expression patterns in both plant species. Although the plastid genome had been lost, plastids likely biosynthesize essential fatty acids and amino acids (aromatic amino acids and lysine). A set of credible and functional horizontal gene transfer (HGT) events (involving genes and mRNAs) were identified in the nuclear and mitochondrial genomes of S. himalayana, most of which were under purifying selection. Convergent HGTs in Cuscuta, Orobanchaceae, and S. himalayana were mainly expressed at the parasite-host interface. Together, these results suggest that HGTs act as a bridge between the parasite and host, assisting the parasite in acquiring nutrients from the host.

CONCLUSIONS

Our results provide new insights into the flower development process and endoparasitic lifestyle of Rafflesiaceae plants. The amount of gene loss in S. himalayana is consistent with the degree of reduction in its body plan. HGT events are common among endoparasites and play an important role in their lifestyle adaptation.

A developmental gradient reveals biosynthetic pathways to eukaryotic toxins in monocot geophytes

Numerous eukaryotic toxins that accumulate in geophytic plants are valuable in the clinic, yet their biosynthetic pathways have remained elusive. A lead example is the >150 Amaryllidaceae alkaloids (AmAs) including galantamine, an FDA-approved treatment for Alzheimer's disease. We show that while AmAs accumulate to high levels in many tissues in daffodils, biosynthesis is localized to nascent, growing tissue at the base of leaves. A similar trend is found for the production of steroidal alkaloids (e.g. cyclopamine) in corn lily. This model of active biosynthesis enabled elucidation of a complete set of biosynthetic genes for the production of AmAs. Taken together, our work sheds light on the developmental and enzymatic logic of diverse alkaloid biosynthesis in daffodil. More broadly, it suggests a paradigm for biosynthesis regulation in monocot geophytes where plants are protected from herbivory through active charging of newly formed cells with eukaryotic toxins that persist as aboveground tissue develops.

Genome-Wide Identification, Expression, and Response to Fusarium Infection of the SWEET Gene Family in Garlic (Allium sativum L.)

Proteins of the SWEET (Sugar Will Eventually be Exported Transporters) family play an important role in plant development, adaptation, and stress response by functioning as transmembrane uniporters of soluble sugars. However, the information on the SWEET family in the plants of the Allium genus, which includes many crop species, is lacking. In this study, we performed a genome-wide analysis of garlic (Allium sativum L.) and identified 27 genes putatively encoding clade I-IV SWEET proteins. The promoters of the A. sativum (As) SWEET genes contained hormone- and stress-sensitive elements associated with plant response to phytopathogens. AsSWEET genes had distinct expression patterns in garlic organs. The expression levels and dynamics of clade III AsSWEET3, AsSWEET9, and AsSWEET11 genes significantly differed between Fusarium-resistant and -susceptible garlic cultivars subjected to F. proliferatum infection, suggesting the role of these genes in the garlic defense against the pathogen. Our results provide insights into the role of SWEET sugar uniporters in A. sativum and may be useful for breeding Fusarium-resistant Allium cultivars.

Identification and characterization of aldehyde dehydrogenase (ALDH) gene superfamily in garlic and expression profiling in response to drought, salinity, and ABA

In response to biotic and abiotic stressors, aldehydes are detoxified and converted to carboxylic acids by aldehyde dehydrogenases (ALDHs), which are enzymes that use NAD+/NADP+ as cofactors. Garlic (Allium sativum L.) has not yet undergone a systematic examination of the ALDH superfamily, despite the genome sequence having been made public. In this investigation, we identified, characterized, and profiled the expression of the garlic ALDH gene family over the entire genome. The ALDH Gene Nomenclature Committee (AGNC) classification was used to classify and name the 34 ALDH genes that were discovered. Except for chromosome 8, all AsALDH genes were dispersed across the chromosomes. AsALDH genes have various localizations, according to predictions about subcellular localization. The AsALDH proteins are more varied and closely related to rice than to Arabidopsis, according to a study of conserved motifs and phylogenetic relationships. The presence of stress modulation pathways is indicated by the abundance of stress-related cis-elements in the AsALDH genes' promoter regions. Analysis of the RNA-seq data showed that AsALDHs expressed differently in various tissues and at various developmental stages. Nine AsALDHs were chosen for study using RT-qPCR, and the results revealed that the majority of the genes were upregulated in response to ABA and downregulated in response to salinity and drought. The results of this study improved our knowledge of the traits, evolutionary background, and biological functions of AsALDHs genes in growth and development.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

IGTCM: An integrative genome database of traditional Chinese medicine plants

Fully understanding traditional Chinese medicines (TCMs) is still challenging because of the extreme complexity of their chemical components and mechanisms of action. The TCM Plant Genome Project aimed to obtain genetic information, determine gene functions, discover regulatory networks of herbal species, and elucidate the molecular mechanisms involved in the disease prevention and treatment, thereby accelerating the modernization of TCMs. A comprehensive database that contains TCM-related information will provide a vital resource. Here, we present an integrative genome database of TCM plants (IGTCM) that contains 14,711,220 records of 83 annotated TCM-related herb genomes, including 3,610,350 genes, 3,534,314 proteins and corresponding coding sequences, and 4,032,242 RNAs, as well as 1033 non-redundant component records for 68 herbs, downloaded and integrated from the GenBank and RefSeq databases. For minimal interconnectivity, each gene, protein, and component was annotated using the eggNOG-mapper tool and Kyoto Encyclopedia of Genes and Genomes database to acquire pathway information and enzyme classifications. These features can be linked across several species and different components. The IGTCM database also provides visualization and sequence similarity search tools for data analyses. These annotated herb genome sequences in IGTCM database are a necessary resource for systematically exploring genes related to the biosynthesis of compounds that have significant medicinal activities and excellent agronomic traits that can be used to improve TCM-related varieties through molecular breeding. It also provides valuable data and tools for future research on drug discovery and the protection and rational use of TCM plant resources. The IGTCM database is freely available at http://yeyn.group:96/.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

Turning Garlic into a Modern Crop: State of the Art and Perspectives

Garlic is cultivated worldwide for the value of its bulbs, but its cultivation is challenged by the infertility of commercial cultivars and the accumulation of pathogens over time, which occurs as a consequence of vegetative (clonal) propagation. In this review, we summarize the state of the art of garlic genetics and genomics, highlighting recent developments that will lead to its development as a modern crop, including the restoration of sexual reproduction in some garlic strains. The set of tools available to the breeder currently includes a chromosome-scale assembly of the garlic genome and multiple transcriptome assemblies that are furthering our understanding of the molecular processes underlying important traits like the infertility, the induction of flowering and bulbing, the organoleptic properties and resistance to various pathogens.

Comparative analysis of two kinds of garlic seedings: qualities and transcriptional landscape

BACKGROUND

Facility cultivation is widely applied to meet the increasing demand for high yield and quality, with light intensity and light quality being major limiting factors. However, how changes in the light environment affect development and quality are unclear in garlic. When garlic seedlings are grown, they can also be exposed to blanching culture conditions of darkness or low-light intensity to ameliorate their appearance and modify their bioactive compounds and flavor.

RESULTS

In this study, we determined the quality and transcriptomes of 14-day-old garlic and blanched garlic seedlings (green seedlings and blanched seedlings) to explore the mechanisms by which seedlings integrate light signals. The findings revealed that blanched garlic seedlings were taller and heavier in fresh weight compared to green garlic seedlings. In addition, the contents of allicin, cellulose, and soluble sugars were higher in the green seedlings. We also identified 3,872 differentially expressed genes between green and blanched garlic seedlings. The Kyoto Encyclopedia of Genes and Genomes analysis suggested enrichment for plant-pathogen interactions, phytohormone signaling, mitogen-activated protein kinase signaling, and other metabolic processes. In functional annotations, pathways related to the growth and formation of the main compounds included phytohormone signaling, cell wall metabolism, allicin biosynthesis, secondary metabolism and MAPK signaling. Accordingly, we identified multiple types of transcription factor genes involved in plant-pathogen interactions, plant phytohormone signaling, and biosynthesis of secondary metabolites among the differentially expressed genes between green and blanched garlic seedlings.

CONCLUSIONS

Blanching culture is one facility cultivation mode that promotes chlorophyll degradation, thus changing the outward appearance of crops, and improves their flavor. The large number of DEGs identified confirmed the difference of the regulatory machinery under two culture system. This study increases our understanding of the regulatory network integrating light and darkness signals in garlic seedlings and provides a useful resource for the genetic manipulation and cultivation of blanched garlic seedlings.

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.

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.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

Integrated Genomic and Transcriptomic Elucidation of Flowering in Garlic

Commercial cultivars of garlic are sterile, and therefore efficient breeding of this crop is impossible. Recent restoration of garlic fertility has opened new options for seed production and hybridization. Transcriptome catalogs were employed as a basis for garlic genetic studies, and in 2020 the huge genome of garlic was fully sequenced. We provide conjoint genomic and transcriptome analysis of the regulatory network in flowering garlic genotypes. The genome analysis revealed phosphatidylethanolamine-binding proteins (PEBP) and LEAFY (LFY) genes that were not found at the transcriptome level. Functions of TFL-like genes were reduced and replaced by FT-like homologs, whereas homologs of MFT-like genes were not found. The discovery of three sequences of LFY-like genes in the garlic genome and confirmation of their alternative splicing suggest their role in garlic florogenesis. It is not yet clear whether AsLFY1 acts alone as the "pioneer transcription factor" or AsLFY2 also provides these functions. The presence of several orthologs of flowering genes that differ in their expression and co-expression network advocates ongoing evolution in the garlic genome and diversification of gene functions. We propose that the process of fertility deprivation in garlic cultivars is based on the loss of transcriptional functions of the specific genes.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

Achievements and perspectives of synthetic biology in botanical insecticides

Botanical insecticides are the origin of all insecticidal compounds. They have been widely used to control pests in crops for a long time. Currently, the commercial production of botanical insecticides extracted from plants is limited because of insufficient raw material supply. Synthetic biology is a promising and effective approach for addressing the current problems of the production of botanical insecticides. It is an emerging biological research hotspot in the field of botanical insecticides. However, the biosynthetic pathways of many botanical insecticides are not completely elucidated. On the other hand, the cytotoxicity of botanical pesticides and low efficiency of these biosynthetic enzymes in new hosts make it still challenging for their heterologous production. In the present review, we summarized the recent developments in the heterologous production of botanical insecticides, analyzed the current challenges, and discussed the feasible production strategies, focusing on elucidating biosynthetic pathways, enzyme engineering, host engineering, and cytotoxicity engineering. Looking to the future, synthetic biology promises to further advance heterologous production of more botanical pesticides.

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.

The genome and gene editing system of sea barleygrass provide a novel platform for cereal domestication and stress tolerance studies

The tribe Triticeae provides important staple cereal crops and contains elite wild species with wide genetic diversity and high tolerance to abiotic stresses. Sea barleygrass (Hordeum marinum Huds.), a wild Triticeae species, thrives in saline marshlands and is well known for its high tolerance to salinity and waterlogging. Here, a 3.82-Gb high-quality reference genome of sea barleygrass is assembled de novo, with 3.69 Gb (96.8%) of its sequences anchored onto seven chromosomes. In total, 41 045 high-confidence (HC) genes are annotated by homology, de novo prediction, and transcriptome analysis. Phylogenetics, non-synonymous/synonymous mutation ratios (Ka/Ks), and transcriptomic and functional analyses provide genetic evidence for the divergence in morphology and salt tolerance among sea barleygrass, barley, and wheat. The large variation in post-domestication genes (e.g. IPA1 and MOC1) may cause interspecies differences in plant morphology. The extremely high salt tolerance of sea barleygrass is mainly attributed to low Na⁺ uptake and root-to-shoot translocation, which are mainly controlled by SOS1, HKT, and NHX transporters. Agrobacterium-mediated transformation and CRISPR/Cas9-mediated gene editing systems were developed for sea barleygrass to promote its utilization for exploration and functional studies of hub genes and for the genetic improvement of cereal crops.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

Genomic insights into the evolutionary history and diversification of bulb traits in garlic

Background

Garlic is an entirely sterile crop with important value as a vegetable, condiment, and medicine. However, the evolutionary history of garlic remains largely unknown.

Results

Here we report a comprehensive map of garlic genomic variation, consisting of amazingly 129.4 million variations. Evolutionary analysis indicates that the garlic population diverged at least 100,000 years ago, and the two groups cultivated in China were domesticated from two independent routes. Consequently, 15.0 and 17.5% of genes underwent an expression change in two cultivated groups, causing a reshaping of their transcriptomic architecture. Furthermore, we find independent domestication leads to few overlaps of deleterious substitutions in these two groups due to separate accumulation and selection-based removal. By analysis of selective sweeps, genome-wide trait associations and associated transcriptomic analysis, we uncover differential selections for the bulb traits in these two garlic groups during their domestication.

Conclusions

This study provides valuable resources for garlic genomics-based breeding, and comprehensive insights into the evolutionary history of this clonal-propagated crop.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02756-1.

WUSCHEL-related homeobox genes cooperate with cytokinin to promote bulbil formation in Lilium lancifolium

The bulbil is an important vegetative reproductive organ in triploid tiger lily (Lilium lancifolium). Based on our previously obtained transcriptome data, we screened two WUSCHEL-related homeobox (WOX) genes closely related to bulbil formation, LlWOX9 and LlWOX11. However, the biological functions and regulatory mechanisms of LlWOX9 and LlWOX11 are unclear. In this study, we cloned the full-length coding sequences of LlWOX9 and LlWOX11. Transgenic Arabidopsis (Arabidopsis thaliana) showed increased branch numbers, and the overexpression of LlWOX9 and LlWOX11 in stem segments promoted bulbil formation, while the silencing of LlWOX9 and LlWOX11 inhibited bulbil formation, indicating that LlWOX9 and LlWOX11 are positive regulators of bulbil formation. Cytokinin type-B response regulators could bind to the promoters of LlWOX9 and LlWOX11 and promote their transcription. LlWOX11 could enhance cytokinin pathway signaling by inhibiting the transcription of type-A LlRR9. Our study enriches the understanding of the regulation of plant development by the WOX gene family and lays a foundation for further research on the molecular mechanism of bulbil formation in lily.

The chromosome-level genome for Toxicodendron vernicifluum provides crucial insights into Anacardiaceae evolution and urushiol biosynthesis

The lacquer tree (Toxicodendron vernicifluum (Stokes) F.A. Barkley) is an important tree with economic, industrial, and medicinal values. Here, we generated the reference genome of T. vernicifluum at the chromosome level with 491.93 Mb in size, in which 98.26% of the assembled contigs were anchored onto 15 pseudochromosomes with the scaffold N50 of 32.97 Mb. Comparative genomic analysis revealed the gene families related to urushiol biosynthesis were expanded, contributing to the ecological fitness and biological adaptability of the lacquer tree. We combined multi-omics data to identify genes that encode key enzymes in the T. vernicifluum urushiol and lignin biosynthetic pathways. Furthermore, the unique active metabolites, such as butin and fisetin, in cultivar lacquers were identified by metabolism profiling. Our work would provide crucial insights into metabolite synthesis such as urushiol and lignin, meanwhile offer a basis for further exploration of the cultivation and breeding of T. vernicifluum and other Anacardiaceae members.

•

We provide the first chromosome-level genome for Toxicodendron vernicifluum

•

We explore the phylogenetic position of lacquer tree

•

We identify genes involved in the urushiol and lignin biosynthetic pathways

•

Analysis of the unique active metabolites in cultivar lacquers by LC-MS/MS is reported

Research Progress on the Synthetic Biology of Botanical Biopesticides

The production and large-scale application of traditional chemical pesticides will bring environmental pollution and food safety problems. With the advantages of high safety and environmental friendliness, botanical biopesticides are in line with the development trend of modern agriculture and have gradually become the mainstream of modern pesticide development. However, the traditional production of botanical biopesticides has long been faced with prominent problems, such as limited source and supply, complicated production processes, and excessive consumption of resources. In recent years, the rapid development of synthetic biology will break through these bottlenecks, and many botanical biopesticides are produced using synthetic biology, such as emodin, celangulin, etc. This paper reviews the latest progress and application prospect of synthetic biology in the development of botanical pesticides so as to provide new ideas for the analysis of synthetic pathways and heterologous and efficient production of botanical biopesticides and accelerate the research process of synthetic biology of natural products.

Single-molecule real-time sequencing of the full-length transcriptome of purple garlic (Allium sativum L. cv. Leduzipi) and identification of serine O-acetyltransferase family proteins involved in cysteine biosynthesis

BACKGROUND

Garlic (Allium sativum L.), whose bioactive components are mainly organosulfur compounds (OSCs), is a herbaceous perennial widely consumed as a green vegetable and a condiment. Yet, the metabolic enzymes involved in the biosynthesis of OSCs are not identified in garlic.

RESULTS

Here, a full-length transcriptome of purple garlic was generated via PacBio and Illumina sequencing, to characterize the garlic transcriptome and identify key proteins mediating the biosynthesis of OSCs. Overall, 22.56 Gb of clean data were generated, resulting in 454 698 circular consensus sequence (CCS) reads, of which 83.4% (379 206) were identified as being full-length non-chimeric reads - their further transcript clustering facilitated identification of 36 571 high-quality consensus reads. Once corrected, their genome-wide mapping revealed that 6140 reads were novel isoforms of known genes, and 2186 reads were novel isoforms from novel genes. We detected 1677 alternative splicing events, finding 2902 genes possessing either two or more poly(A) sites. Given the importance of serine O-acetyltransferase (SERAT) in cysteine biosynthesis, we investigated the five SERAT homologs in garlic. Phylogenetic analysis revealed a three-tier classification of SERAT proteins, each featuring a serine acetyltransferase domain (N-terminal) and one or two hexapeptide transferase motifs. Template-based modeling showed that garlic SERATs shared a common homo-trimeric structure with homologs from bacteria and other plants. The residues responsible for substrate recognition and catalysis were highly conserved, implying a similar reaction mechanism. In profiling the five SERAT genes' transcript levels, their expression pattern varied significantly among different tissues.

CONCLUSION

This study's findings deepen our knowledge of SERAT proteins, and provide timely genetic resources that could advance future exploration into garlic's genetic improvement and breeding. © 2021 Society of Chemical Industry.

Pretreatment with Chitosan Prevents Fusarium Infection and Induces the Expression of Chitinases and β-1,3-Glucanases in Garlic (Allium sativum L.)

Fusarium infection decreases the yield of garlic (Allium sativum L.); however, the knowledge about garlic response to fungal attack is limited. Chitosan induces plant defense response to stress conditions. Here, we analyzed the effects of chitosan with low (Ch1, 39 kDa) and medium (Ch2, 135 kDa) molecular weight on Fusarium infection in garlic. Ch1 and Ch2 at concentrations 0.125–0.400 mg/mL suppressed the growth of Fusarium proliferatum cultures in vitro. Pretreatment of garlic bulbs with Ch1 or Ch2 prevented disease symptoms after F. proliferatum inoculation, while exerting early inhibitory and late stimulatory effects on chitinase and β-1,3-glucanase activities. Ch1/Ch2 treatment of garlic already infected with F. proliferatum caused transcriptional upregulation of chitinases and β-1,3-glucanases at the early stage, which was maintained at the late stage in Ch2-treated samples, but not in Ch1-treated samples, where transcriptional inhibition was observed. The stimulatory effect of Ch2 pretreatment on the expression of chitinase and endo-β-1,3-glucanase genes was stronger than that of Ch1 pretreatment, suggesting that Ch2 could be more effective than Ch1 in pre-sowing treatment of garlic bulbs. Our results provide insights into the effects of chitosan on the garlic response to Fusarium, suggesting a novel strategy to protect garlic crop against fungal infection.

Genome-wide Identification and Characterization of the GRAS Transcription Factors in Garlic (Allium sativum L.)

GRAS transcription factors play crucial roles in plant growth and development and have been widely explored in many plant species. Garlic (Allium sativum L.) is an important crop owing to its edible and medicinal properties. However, no GRAS transcription factors have been identified in this crop. In this study, 46 garlic GRAS genes were identified and assigned to 16 subfamilies using the GRAS members of Arabidopsis thaliana, Oryza sativa, and Amborella trichopoda as reference queries. Expression analysis revealed that garlic GRAS genes showed distinct differences in various garlic tissues, as well as during different growth stages of the bulbs. Five of these 46 genes were identified as DELLA-like protein-encoding genes and three of which, Asa2G00237.1/Asa2G00240.1 and Asa4G02090.1, responded to exogenous GA3 treatment, and showed a significant association between their transcription abundance and bulb traits in 102 garlic accessions, thereby indicating their role in regulating the growth of garlic bulbs. These results will lay a useful foundation for further investigation of the biological functions of GRAS genes and guiding the genetic breeding of garlic in the future.