Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genome evolution of the ancient hexaploid Platanus × acerifolia (London planetree)

Whole-genome duplication (WGD; i.e., polyploidy) and chromosomal rearrangement (i.e., genome shuffling) significantly influence genome structure and organization. Many polyploids show extensive genome shuffling relative to their pre-WGD ancestors. No reference genome is currently available for Platanaceae (Proteales), one of the sister groups to the core eudicots. Moreover, Platanus × acerifolia (London planetree; Platanaceae) is a widely used street tree. Given the pivotal phylogenetic position of Platanus and its 2-y flowering transition, understanding its flowering-time regulatory mechanism has significant evolutionary implications; however, the impact of Platanus genome evolution on flowering-time genes remains unknown. Here, we assembled a high-quality, chromosome-level reference genome for P. × acerifolia using a phylogeny-based subgenome phasing method. Comparative genomic analyses revealed that P. × acerifolia (2n = 42) is an ancient hexaploid with three subgenomes resulting from two sequential WGD events; Platanus does not seem to share any WGD with other Proteales or with core eudicots. Each P. × acerifolia subgenome is highly similar in structure and content to the reconstructed pre-WGD ancestral eudicot genome without chromosomal rearrangements. The P. × acerifolia genome exhibits karyotypic stasis and gene sub-/neo-functionalization and lacks subgenome dominance. The copy number of flowering-time genes in P. × acerifolia has undergone an expansion compared to other noncore eudicots, mainly via the WGD events. Sub-/neo-functionalization of duplicated genes provided the genetic basis underlying the unique flowering-time regulation in P. × acerifolia. The P. × acerifolia reference genome will greatly expand understanding of the evolution of genome organization, genetic diversity, and flowering-time regulation in angiosperms.

An Updated Review on Nelumbo Nucifera Gaertn: Chemical Composition, Nutritional Value and Pharmacological Activities

Nelumbo nucifera Gaertn is a recognised herbal plant in ancient medical sciences. Each portion of the plant leaf, flower, seed and rhizome is utilised for nutritional and medicinal purposes. The chemical compositions like phenol, alkaloids, glycoside, terpenoids and steroids have been isolated. The plant contains various nutritional values like lipids, proteins, amino acids, minerals, carbohydrates, and fatty acids. Traditional medicine confirms that the phytochemicals of plants give significant benefits to the treatment of various diseases such as leukoderma, smallpox, dysentery, haematemesis, coughing, haemorrhage, metrorrhagia, haematuria, fever, hyperlipidaemia, cholera, hepatopathy and hyperdipsia. To verify the traditional claims, researchers have conducted scientific biological in vivo and in vitro screenings, which have exhibited that the plant keeps various notable pharmacological activities such as anticancer, hepatoprotective, antioxidant, antiviral, hypolipidemic, anti-obesity, antipyretic, hypoglycaemic, antifungal, anti-inflammatory and antibacterial activities. This review, summaries the nutritional composition, chemical constituents and biological activities substantiated by the researchers done in vivo and in vitro.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

Metabolite profiling and screening of callus browning-related genes in lotus (Nelumbo nucifera)

Callus browning is a major drawback to lotus callus proliferation and regeneration. However, the underlying mechanism of its formation remains largely unknown. Herein, we aimed to explore the metabolic and molecular basis of lotus callus browning by combining histological staining, high-throughput metabolomics, and transcriptomic assays for lotus callus at three browning stages. Histological stained brown callus cross sections displayed severe cell death symptoms, accompanied by an obvious accumulation of polyphenols and lignified materials. Widely targeted metabolomics revealed extensively decreased accumulation of most detected flavonoids and benzylisoquinoline alkaloids (BIAs), as well as a few phenolic acids, amino acids and their derivatives in callus with browning symptoms. Conversely, the contents of most detected tannins were significantly increased. Subsequent comparative transcriptomics identified a set of differentially expressed genes (DEGs) associated with the biosynthesis and regulation of flavonoids and BIAs in lotus. Notably, callus browning was coupled with significantly up-regulated expression of two polyphenol oxidase (PPO) and 17 peroxidase (POD) encoding genes, while the expression of ethylene associated genes remained at marginal levels. These results suggest that lotus callus browning is primarily controlled at the level of metabolism, wherein the oxidation of flavonoids and BIAs is crucially decisive.

Identification of QTLs and Putative Candidate Genes for Plant Architecture of Lotus Revealed by Regional Association Mapping

The lotus (Nelumbo Adans.) is one of the most economically relevant ornamental aquatic plants. Plant architecture (PA) is an important trait for lotus classification, cultivation, breeding, and applications. However, the underlying genetic and molecular basis controlling PA remains poorly understood. In this study, an association study for PA-related traits was performed with 93 genome-wide microsatellite markers (simple sequence repeat, SSR) and 51 insertion-deletion (InDel) markers derived from the candidate regions using a panel of 293 lotus accessions. Phenotypic data analysis of the five PA-related traits revealed a wide normal distribution and high heritability from 2013 to 2016, which indicated that lotus PA-related traits are highly polygenic traits. The population structure (Q-matrix) and the relative kinships (K-matrix) of the association panels were analyzed using 93 SSR markers. The mixed linear model (MLM) taking Q-matrix and K-matrix into account was used to estimate the association between markers and the traits. A total of 26 markers and 65 marker-trait associations were identified by considering associations with p < 0.001 and Q < 0.05. Based on the significant markers, two QTLs on Chromosome 1 were identified, and two candidate genes were preliminarily determined. The results of our study provided useful information for the lotus breeding aiming at different PA phenotypes using a molecular-assisted selection (MAS) method and also laid the foundation for the illustration of the molecular mechanism underlying the major QTL and key markers associated with lotus PA.

Elucidation of the (R)-enantiospecific benzylisoquinoline alkaloid biosynthetic pathways in sacred lotus (Nelumbo nucifera)

Benzylisoquinoline alkaloids (BIAs) are a structurally diverse group of plant specialized metabolites found mainly in members of the order Ranunculales, including opium poppy (Papaver somniferum), for which BIA biosynthetic pathways leading to the critical drugs morphine, noscapine, and sanguinarine have been elucidated. Sacred lotus (Nelumbo nucifera), in the order Proteales, accumulates medicinal BIAs in the proaporphine, aporphine, and bisbenzylisoquinoline structural subgroups with a prevalence of R enantiomers, opposed to the dominant S configuration occurring in the Ranunculales. Nevertheless, distinctive BIA biosynthetic routes in sacred lotus have not been explored. In planta labeling experiments and in vitro assays with recombinant enzymes and plant protein extracts showed that dopamine and 4-hydroxyphenylacetaldehyde derived from L-tyrosine serve as precursors for the formation of (R,S)-norcoclaurine in sacred lotus, whereas only (R)-norcoclaurine byproducts are favored in the plant by action of R-enantiospecific methyltransferases and cytochrome P450 oxidoreductases (CYPs). Enzymes responsible for the R-enantiospecific formation of proaporphine (NnCYP80Q1) and bisbenzylisoquinoline (NnCYP80Q2) scaffolds, and a methylenedioxy bridge introduction on aporphine substrates (NnCYP719A22) were identified, whereas additional aspects of the biosynthetic pathways leading to the distinctive alkaloid profile are discussed. This work expands the availability of molecular tools that can be deployed in synthetic biology platforms for the production of high-value alkaloids.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

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.

Tree2GD: a phylogenomic method to detect large-scale gene duplication events

MOTIVATION

Whole-genome duplication events have long been discovered throughout the evolution of eukaryotes, contributing to genome complexity and biodiversity and leaving traces in the descending organisms. Therefore, an accurate and rapid phylogenomic method is needed to identify the retained duplicated genes on various lineages across the target taxonomy.

RESULTS

Here, we present Tree2GD, an integrated method to identify large-scale gene duplication events by automatically perform multiple procedures, including sequence alignment, recognition of homolog, gene tree/species tree reconciliation, Ks distribution of gene duplicates and synteny analyses. Application of Tree2GD on 2 datasets, 12 metazoan genomes and 68 angiosperms, successfully identifies all reported whole-genome duplication events exhibited by these species, showing effectiveness and efficiency of Tree2GD on phylogenomic analyses of large-scale gene duplications.

AVAILABILITY AND IMPLEMENTATION

Tree2GD is written in Python and C++ and is available at https://github.com/Dee-chen/Tree2gd.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

Background

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

Results

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

Conclusions

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01420-1.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

Evaluation of Intracellular Gene Transfers from Plastome to Nuclear Genome across Progressively Improved Assemblies for Arabidopsis thaliana and Oryza sativa

DNA originating from organellar genomes are regularly discovered in nuclear sequences during genome assembly. Nevertheless, such insertions are sometimes omitted during the process of nuclear genome assembly because the inserted DNA is assigned to organellar genomes, leading to a systematic underestimation of their frequency. With the rapid development of high-throughput sequencing technology, more inserted fragments from organelle genomes can now be detected. Therefore, it is necessary to be aware of the insertion events from organellar genomes during nuclear genome assembly to properly attribute the impact and rate of such insertions in the evolution of nuclear genomes. Here, we investigated the impact of intracellular gene transfer (IGT) from the plastome to the nuclear genome using genome assemblies that were refined through time with technological improvements from two model species, Arabidopsis thaliana and Oryza sativa. We found that IGT from the plastome to the nuclear genome is a dynamic and ongoing process in both A. thaliana and O. sativa, and mostly occurred recently, as the majority of transferred sequences showed over 95% sequence similarity with plastome sequences of origin. Differences in the plastome-to-nuclear genome IGT between A. thaliana and O. sativa varied among the different assembly versions and were associated with the quality of the nuclear genome assembly. IGTs from the plastome to nuclear genome occurred more frequently in intergenic regions, which were often associated with transposable elements (TEs). This study provides new insights into intracellular genome evolution and nuclear genome assembly by characterizing and comparing IGT from the plastome into the nuclear genome for two model plant species.

Identification of QTLs and a putative candidate gene involved in rhizome enlargement of Asian lotus (Nelumbo nucifera)

QTL mapping studies identified three reliable QTLs of rhizome enlargement in lotus. NnBEL6 located within the confidence interval of the major QTL cqREI-LG2 is a key candidate gene enhancing rhizome enlargement. Lotus (Nelumbo) is perennial aquatic plant with nutritional, pharmacological, and ornamental significance. Rhizome is an underground lotus stem that acts as a storage organ and as a reproductive tissue for asexual production. The enlargement of lotus rhizome is an important adaptive strategy for surviving the cold winter. The aims of this study were to identify quantitative trait loci (QTLs) for rhizome enlargement traits including rhizome enlargement index (REI) and number of enlarged rhizome (NER), and to uncover their associated candidate genes. A high-density genetic linkage map was constructed, consisting of 2935 markers binned from 236,840 SNPs. A total of 14 significant QTLs were detected for REI and NER, which explained 6.7-22.3% of trait variance. Three QTL regions were repeatedly identified in at least 2 years, and a major QTL, designated cqREI-LG2, with a rhizome-enlargement effect and about 20% of the phenotypic contribution was identified across the 3 climatic years. A candidate NnBEL6 gene located within the confidence interval of cqREI-LG2 was considered to be putatively involved in lotus rhizome enlargement. The expression of NnBEL6 was exclusively induced by rhizome swelling. Sequence comparison of NnBEL6 among lotus cultivars revealed a functional Indel site in its promoter that likely initiates the rhizome enlargement process. Transgenic potato assay was used to confirm the role of NnBEL6 in inducing tuberization. The successful identification QTLs and functional validation of NnBEL6 gene reported in this study will enrich our knowledge on the genetic basis of rhizome enlargement in lotus.

The slow-evolving Acorus tatarinowii genome sheds light on ancestral monocot evolution

Monocots are one of the most diverse groups of flowering plants, and tracing the evolution of their ancestral genome into modern species is essential for understanding their evolutionary success. Here, we report a high-quality assembly of the Acorus tatarinowii genome, a species that diverged early from all the other monocots. Genome-wide comparisons with a range of representative monocots characterized Acorus as a slowly evolved genome with one whole-genome duplication. Our inference of the ancestral monocot karyotypes provides new insights into the chromosomal evolutionary history assigned to modern species and reveals the probable molecular functions and processes related to the early adaptation of monocots to wetland or aquatic habitats (that is, low levels of inorganic phosphate, parallel leaf venation and ephemeral primary roots). The evolution of ancestral gene order in monocots is constrained by gene structural and functional features. The newly obtained Acorus genome offers crucial evidence for delineating the origin and diversification of monocots, including grasses.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Comparative analyses of American and Asian lotus genomes reveal insights into petal color, carpel thermogenesis and domestication

Nelumbo lutea (American lotus), which differs from Nelumbo nucifera (Asian lotus) morphologically, is one of the two remaining species in the basal eudicot family Nelumbonaceae. Here, we assembled the 843-Mb genome of American lotus into eight pseudochromosomes containing 31 382 protein-coding genes. Comparative analyses revealed conserved synteny without large chromosomal rearrangements between the genomes of American and Asian lotus and identified 29 533 structural variants (SVs). Carotenoid and anthocyanin pigments determine the yellow and red petal colors of American and Asian lotus, respectively. The structural genes encoding enzymes of the carotenoid and anthocyanin biosynthesis pathways were conserved between two species but differed in expression. We detected SVs caused by repetitive sequence expansion or contraction among the anthocyanin biosynthesis regulatory MYB genes. Further transient overexpression of candidate NnMYB5 induced anthocyanin accumulation in lotus petals. Alternative oxidase (AOX), uncoupling proteins (UCPs), and sugar metabolism and transportation contributed to carpel thermogenesis. Carpels produce heat with sugars transported from leaves as the main substrates, because there was weak tonoplast sugar transporter (TST) activity, and with SWEETs were highly expressed during thermogenesis. Cell proliferation-related activities were particularly enhanced in the warmer carpels compared with stamens during the cold night before blooming, which suggested that thermogenesis plays an important role in flower protogyny. Population genomic analyses revealed deep divergence between American and Asian lotus, and independent domestication affecting seed, rhizome, and flower traits. Our findings provide a high-quality reference genome of American lotus for exploring the genetic divergence and variation between two species and revealed possible genomic bases for petal color, carpel thermogenesis and domestication in lotus.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.

The adaptive evolution of Euryale ferox to the aquatic environment through paleo-hexaploidization

Occupation of living space is one of the main driving forces of adaptive evolution, especially for aquatic plants whose leaves float on the water surface and thus have limited living space. Euryale ferox, from the angiosperm basal family Nymphaeaceae, develops large, rapidly expanding leaves to compete for space on the water surface. Microscopic observation found that the cell proliferation of leaves is almost completed underwater, while the cell expansion occurs rapidly after they grow above water. To explore the mechanism underlying the specific development of leaves, we performed sequences assembly and analyzed the genome and transcriptome dynamics of E. ferox. Through reconstruction of the three sub-genomes generated from the paleo-hexaploidization event in E. ferox, we revealed that one sub-genome was phylogenetically closer to Victoria cruziana, which also exhibits gigantic floating leaves. Further analysis revealed that while all three sub-genomes promoted the evolution of the specific leaf development in E. ferox, the genes from the sub-genome closer to V. cruziana contributed more to this adaptive evolution. Moreover, we found that genes involved in cell proliferation and expansion, photosynthesis, and energy transportation were over-retained and showed strong expression association with the leaf development stages, such as the expression divergence of SWEET orthologs as energy uploaders and unloaders in the sink and source leaf organs of E. ferox. These findings provide novel insights into the genome evolution through polyploidization, as well as the adaptive evolution regarding the leaf development accomplished through biased gene retention and expression sub-functionalization of multi-copy genes in E. ferox.

Twenty years of plant genome sequencing: achievements and challenges

Publication of the complete genome sequence of Arabidopsis thaliana, the first plant reference genome, in December 2000 heralded the beginning of the plant genome era. Over the past 20 years reference genomes have been generated for hundreds of plant species, spanning non-vascular to flowering plants. Releasing these plant genomes has dramatically advanced studies in all disciplines of plant biology. Importantly, multiple reference-level genomes have been generated for the major crops and their progenitors, enabling the creation of pan-genomes and exploration of domestication history and natural variations that can be adopted by modern crop breeding. We summarize the progress of plant genome sequencing and the challenges of sequencing more complex plant genomes and generating pan-genomes.

Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans)

Lotus (Nelumbo Adans.), a relict plant, is the testimony of long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remains unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus, respectively. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10-9 base/generation in lotus shoot stem cells. The past history of lotus revealed that the ancestors of lotus in eastern and southern Asia could be traced back ~20 million years ago (Mya) and experienced twice significant bottlenecks and population splits. We further identified the selected genes among three lotus groups in different habitats, suggesting that 453 genes between tropical and temperate group and 410 genes between two subgroups from Northeastern China and the Yangtze River - Yellow River Basin might play important roles in natural selection in lotus's adaptation and resilience. Our findings not only improve an understanding of the lotus evolutionary history and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for the relict plants.

Jasmonate-Responsive Transcription Factors NnWRKY70a and NnWRKY70b Positively Regulate Benzylisoquinoline Alkaloid Biosynthesis in Lotus (Nelumbo nucifera)

Lotus (Nelumbo nucifera) is a large aquatic plant that accumulates pharmacologically significant benzylisoquinoline alkaloids (BIAs). However, little is known about their biosynthesis and regulation. Here, we show that the two group III WRKY transcription factors (TFs), NnWRKY70a and NnWRKY70b, positively regulate the BIA biosynthesis in lotus. Both NnWRKY70s are jasmonic acid (JA) responsive, with their expression profiles highly correlated to the BIA concentration and BIA pathway gene expression. A dual-luciferase assay showed that NnWRKY70a could transactivate the NnTYDC promoter, whereas NnWRKY70b could activate promoters of the three BIA structural genes, including NnTYDC, NnCYP80G, and Nn7OMT. In addition, the transient overexpression of NnWRKY70a and NnWRKY70b in lotus petals significantly elevated the BIA alkaloid concentrations. Notably, NnWRKY70b seems to be a stronger BIA biosynthesis regulator, because it dramatically induced more BIA structural gene expressions and BIA accumulation than NnWRKY70a. A yeast two-hybrid assay further revealed that NnWRKY70b physically interacted with NnJAZ1 and two other group III WRKY TFs (NnWRKY53b and NnWRKY70a), suggesting that it may cooperate with the other group III WRKYs to adjust the lotus BIA biosynthesis via the JA-signaling pathway. To illustrate the mechanism underlying NnWRKY70b-mediated BIA regulation in the lotus, a simplified model is proposed. Our study provides useful insights into the regulatory roles of WRKY TFs in the biosynthesis of secondary metabolites.

Characterization of Genomic Variation from Lotus (Nelumbo Adans.) Mutants with Wide and Narrow Tepals

Compared with rose, chrysanthemum, and water lily, the absence of short-wide and long-narrow tepals of ornamental lotus (Nelumbo Adans.) limits the commercial value of flowers. In this study, the genomes of two groups of lotus mutants with wide-short and narrow-long tepals were resequenced to uncover the genomic variation and candidate genes associated with tepal shape. In group NL (short for N. lutea, containing two mutants and one control of N. lutea), 716,656 single nucleotide polymorphisms (SNPs) and 221,688 insertion-deletion mutations (Indels) were obtained, while 639,953 SNPs and 134,6118 Indels were obtained in group WSH (short for ‘Weishan Hong’, containing one mutant and two controls of N. nucifera ‘Weishan Hong’). Only a small proportion of these SNPs and Indels was mapped to exonic regions of genome: 1.92% and 0.47%, respectively, in the NL group, and 1.66% and 0.48%, respectively, in the WSH group. Gene Ontology (GO) analysis showed that out of 4890 (NL group) and 1272 (WSH group) annotated variant genes, 125 and 62 genes were enriched (Q < 0.05), respectively. Additionally, in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 104 genes (NL group) and 35 genes (WSH group) were selected (p < 0.05). Finally, there were 306 candidate genes that were sieved to determine the development of tepal shape in lotus plants. It will be an essential reference for future identification of tepal-shaped control genes in lotus plants. This is the first comprehensive report of genomic variation controlling tepal shape in lotus, and the mutants in this study are promising materials for breeding novel lotus cultivars with special tepals.

BREC: an R package/Shiny app for automatically identifying heterochromatin boundaries and estimating local recombination rates along chromosomes

Background

Meiotic recombination is a vital biological process playing an essential role in genome's structural and functional dynamics. Genomes exhibit highly various recombination profiles along chromosomes associated with several chromatin states. However, eu-heterochromatin boundaries are not available nor easily provided for non-model organisms, especially for newly sequenced ones. Hence, we miss accurate local recombination rates necessary to address evolutionary questions.

Results

Here, we propose an automated computational tool, based on the Marey maps method, allowing to identify heterochromatin boundaries along chromosomes and estimating local recombination rates. Our method, called BREC (heterochromatin Boundaries and RECombination rate estimates) is non-genome-specific, running even on non-model genomes as long as genetic and physical maps are available. BREC is based on pure statistics and is data-driven, implying that good input data quality remains a strong requirement. Therefore, a data pre-processing module (data quality control and cleaning) is provided. Experiments show that BREC handles different markers' density and distribution issues.

Conclusions

BREC's heterochromatin boundaries have been validated with cytological equivalents experimentally generated on the fruit fly Drosophila melanogaster genome, for which BREC returns congruent corresponding values. Also, BREC's recombination rates have been compared with previously reported estimates. Based on the promising results, we believe our tool has the potential to help bring data science into the service of genome biology and evolution. We introduce BREC within an R-package and a Shiny web-based user-friendly application yielding a fast, easy-to-use, and broadly accessible resource. The BREC R-package is available at the GitHub repository https://github.com/GenomeStructureOrganization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04233-1.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

Re-examination of two diatom reference genomes using long-read sequencing

Background

The marine diatoms Thalassiosira pseudonana and Phaeodactylum tricornutum are valuable model organisms for exploring the evolution, diversity and ecology of this important algal group. Their reference genomes, published in 2004 and 2008, respectively, were the product of traditional Sanger sequencing. In the case of T. pseudonana, optical restriction site mapping was employed to further clarify and contextualize chromosome-level scaffolds. While both genomes are considered highly accurate and reasonably contiguous, they still contain many unresolved regions and unordered/unlinked scaffolds.

Results

We have used Oxford Nanopore Technologies long-read sequencing to update and validate the quality and contiguity of the T. pseudonana and P. tricornutum genomes. Fine-scale assessment of our long-read derived genome assemblies allowed us to resolve previously uncertain genomic regions, further characterize complex structural variation, and re-evaluate the repetitive DNA content of both genomes. We also identified 1862 previously undescribed genes in T. pseudonana. In P. tricornutum, we used transposable element detection software to identify 33 novel copia-type LTR-RT insertions, indicating ongoing activity and rapid expansion of this superfamily as the organism continues to be maintained in culture. Finally, Bionano optical mapping of P. tricornutum chromosomes was combined with long-read sequence data to explore the potential of long-read sequencing and optical mapping for resolving haplotypes.

Conclusion

Despite its potential to yield highly contiguous scaffolds, long-read sequencing is not a panacea. Even for relatively small nuclear genomes such as those investigated herein, repetitive DNA sequences cause problems for current genome assembly algorithms. Determining whether a long-read derived genomic assembly is ‘better’ than one produced using traditional sequence data is not straightforward. Our revised reference genomes for P. tricornutum and T. pseudonana nevertheless provide additional insight into the structure and evolution of both genomes, thereby providing a more robust foundation for future diatom research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07666-3.

Research advances in and prospects of ornamental plant genomics

The term 'ornamental plant' refers to all plants with ornamental value, which generally have beautiful flowers or special plant architectures. China is rich in ornamental plant resources and known as the "mother of gardens". Genomics is the science of studying genomes and is useful for carrying out research on genome evolution, genomic variations, gene regulation, and important biological mechanisms based on detailed genome sequence information. Due to the diversity of ornamental plants and high sequencing costs, the progress of genome research on ornamental plants has been slow for a long time. With the emergence of new sequencing technologies and a reduction in costs since the whole-genome sequencing of the first ornamental plant (Prunus mume) was completed in 2012, whole-genome sequencing of more than 69 ornamental plants has been completed in <10 years. In this review, whole-genome sequencing and resequencing of ornamental plants will be discussed. We provide analysis with regard to basic data from whole-genome studies of important ornamental plants, the regulation of important ornamental traits, and application prospects.

Genetic dissection of rhizome yield-related traits in Nelumbo nucifera through genetic linkage map construction and QTL mapping

Lotus (Nelumbo nucifera) is a perennial aquatic plant with great value in ornamentation, nutrition, and medicine. Being a storage organ, lotus rhizome is not only used for vegetative reproduction, but also as a popular vegetable in Southeast Asia. Rhizome development, especially enlargement, largely determines its yield and hence becomes one of the major concerns in rhizome lotus breeding and cultivation. To obtain the genetic characteristic of this trait, and discover markers or genes associated with this trait, an F₂ population was generated by crossing between temperate and tropical cultivars with contrasting rhizome enlargement. Based on this F₂ population and Genotyping-by-Sequencing (GBS) technique, a genetic map was constructed with 1475 bin markers containing 12,113 SNP markers. Six traits associated with rhizome yield were observed over 3 years. Quantitative trait locus (QTL) mapping analysis identified 22 QTLs that are associated with at least one of these traits, among which 9 were linked with 3 different intervals. Comparison of the genes located in these three intervals with our previous transcriptomic data showed that light and phytohormone signaling might contribute to the development and enlargement of lotus rhizome. The QTLs obtained here could also be used for marker-assisted breeding of rhizome lotus.

A chromosome-level genome assembly for the tertiary relict plant Tetracentron sinense oliv. (trochodendraceae)

Tetracentron sinense and Trochodendron aralioides are two Tertiary relict species of large trees in the family Trochodendraceae with narrow distributions on the mainland and islands of eastern Asia. They belong to the order Trochodendrales, which is one of the four early-diverged eudicot lineages. These two relict species provide a good system in which to examine genomic changes that occurred as they survived during repeated climatic oscillations in the Quaternary. We sequenced the genome of Te. sinense and compared it with that of Tr. aralioides. We found that Te. sinense has a smaller genome size (986.3 Mb) than that of Tr. aralioides (1610 Mb). Repetitive elements made the major contribution to the contrasting genome sizes in the two species, with most bursts of repeats occurring within the past four million years when the climate oscillated greatly. These species share two rounds of whole-genome duplications. The mainland species Te. sinense had a larger effective population size than the island species Tr. aralioides after the largest glaciation during the Quaternary climatic oscillation. However, soon after this recovery stage, the effective population sizes of both species continued to decrease, although the current effective population size of Te. sinense is still larger than that of Tr. aralioides. We recovered three distinctly diverged clades through resequencing the genomes of 50 individuals across the distributional range of Te. sinense in China. Our results provide an important genomic resource with which to examine early trait evolution in the core eudicots and assist efforts to conserve this relict tree species.

The Tetracentron genome provides insight into the early evolution of eudicots and the formation of vessel elements

Background

Tetracentron sinense is an endemic and endangered deciduous tree. It belongs to the Trochodendrales, one of four early diverging lineages of eudicots known for having vesselless secondary wood. Sequencing and resequencing of the T. sinense genome will help us understand eudicot evolution, the genetic basis of tracheary element development, and the genetic diversity of this relict species.

Results

Here, we report a chromosome-scale assembly of the T. sinense genome. We assemble the 1.07 Gb genome sequence into 24 chromosomes and annotate 32,690 protein-coding genes. Phylogenomic analyses verify that the Trochodendrales and core eudicots are sister lineages and showed that two whole-genome duplications occurred in the Trochodendrales approximately 82 and 59 million years ago. Synteny analyses suggest that the γ event, resulting in paleohexaploidy, may have only happened in core eudicots. Interestingly, we find that vessel elements are present in T. sinense, which has two orthologs of AtVND7, the master regulator of vessel formation. T. sinense also has several key genes regulated by or regulating TsVND7.2 and their regulatory relationship resembles that in Arabidopsis thaliana. Resequencing and population genomics reveals high levels of genetic diversity of T. sinense and identifies four refugia in China.

Conclusions

The T. sinense genome provides a unique reference for inferring the early evolution of eudicots and the mechanisms underlying vessel element formation. Population genomics analysis of T. sinense reveals its genetic diversity and geographic structure with implications for conservation.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-020-02198-7.

Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history

Lotus (family: Nelumbonaceae) are perennial aquatic plants that represent one of the most ancient basal dicots. In the present study, we resequenced 296 lotus accessions from various geographical locations and germplasms to explore their genomic diversity and population structure. This germplasm set consisted of four accessions of American wild lotus and 292 accessions of Asian lotus, which were divided into four subgroups: wild, rhizome, flower and seed. Total single nucleotide polymorphisms (SNPs) suggested that the wild lotus had the highest variant number (7 191 010). Population structure and genome diversity analysis indicated that the American wild lotus demonstrated a distant genetic relationship with the Asian lotus. Furthermore, the seed and rhizome lotus groups had not originated from a single source but rather had a more complex multisource origin. Besides that, the seed lotus showed higher genetic diversity, which might have been due to the gene flow from the flower lotus to seed lotus by artificial crossing, and the rhizome lotus showed a much lower genetic diversity than the other groups. The present study provides SNP markers for lotus genomic diversity analysis, which will be useful for guiding lotus breeding.

The Establishment of an Efficient Callus Induction System for Lotus (Nelumbo nucifera)

The lotus (Nelumbo nucifera) is one of the most popular aquatic plants in Asia, and has emerged as a novel model for studying flower and rhizome development, and primary and secondary metabolite accumulation. Here, we developed a highly efficient callus induction system for the lotus by optimizing a series of key factors that affect callus formation. The highest efficient callus production was induced on immature cotyledon and embryo explants grown on Murashige and Skoog (MS) basal medium containing an optimized combination of 3 mg/L 2,4-dichlorophenoxyacetic acid (2,4-D) and 0.5 mg/L 6-benzylaminopurine (6-BA). In addition, lotus callus induction was proven to be influenced by lotus genotypes, light conditions, the developmental stages of explants and the time of explant sampling. Collecting immature cotyledons from seeds of the genotype “Shilihe 1”, at 9 days post pollination, and to culture the explants in darkness, are proposed as the optimum conditions for lotus callus induction. Interestingly, highly efficient callus induction was also observed in explants of immature embryo derived aseptic seedlings; and a small amount of lotus benzylisoquinoline alkaloid (BIA) and obvious expression of BIA biosynthetic genes were detected in lotus callus.

Distinct Expression and Methylation Patterns for Genes with Different Fates following a Single Whole-Genome Duplication in Flowering Plants

For most sequenced flowering plants, multiple whole-genome duplications (WGDs) are found. Duplicated genes following WGD often have different fates that can quickly disappear again, be retained for long(er) periods, or subsequently undergo small-scale duplications. However, how different expression, epigenetic regulation, and functional constraints are associated with these different gene fates following a WGD still requires further investigation due to successive WGDs in angiosperms complicating the gene trajectories. In this study, we investigate lotus (Nelumbo nucifera), an angiosperm with a single WGD during the K-pg boundary. Based on improved intraspecific-synteny identification by a chromosome-level assembly, transcriptome, and bisulfite sequencing, we explore not only the fundamental distinctions in genomic features, expression, and methylation patterns of genes with different fates after a WGD but also the factors that shape post-WGD expression divergence and expression bias between duplicates. We found that after a WGD genes that returned to single copies show the highest levels and breadth of expression, gene body methylation, and intron numbers, whereas the long-retained duplicates exhibit the highest degrees of protein-protein interactions and protein lengths and the lowest methylation in gene flanking regions. For those long-retained duplicate pairs, the degree of expression divergence correlates with their sequence divergence, degree in protein-protein interactions, and expression level, whereas their biases in expression level reflecting subgenome dominance are associated with the bias of subgenome fractionation. Overall, our study on the paleopolyploid nature of lotus highlights the impact of different functional constraints on gene fate and duplicate divergence following a single WGD in plant.

Over 100 Million Years of Enzyme Evolution Underpinning the Production of Morphine in the Papaveraceae Family of Flowering Plants

Phylogenomic analysis of whole genome sequences of five benzylisoquinoline alkaloid (BIA)-producing species from the Ranunculales and Proteales orders of flowering plants revealed the sequence and timing of evolutionary events leading to the diversification of these compounds. (S)-Reticuline is a pivotal intermediate in the synthesis of many BIAs and our analyses revealed parallel evolution between the two orders, which diverged ∼122 million years ago (MYA). Berberine is present in species across the entire Ranunculales, and we found co-evolution of genes essential for production of the protoberberine class. The benzophenanthridine class, which includes the antimicrobial compound sanguinarine, is specific to the Papaveraceae family of Ranunculales, and biosynthetic genes emerged after the split with the Ranunculaceae family ∼110 MYA but before the split of the three Papaveraceae species used in this study at ∼77 MYA. The phthalideisoquinoline noscapine and morphinan class of BIAs are exclusive to the opium poppy lineage. K_s estimation of paralogous pairs indicates that morphine biosynthesis evolved more recently than 18 MYA in the Papaver genus. In the preceding 100 million years gene duplication, neofunctionalization and recruitment of additional enzyme classes, combined with gene clustering, gene fusion, and gene amplification, resulted in emergence of medicinally valuable BIAs including morphine and noscapine.

Isolation and characterization of two O-methyltransferases involved in benzylisoquinoline alkaloid biosynthesis in sacred lotus (Nelumbo nucifera)

Benzylisoquinoline alkaloids (BIAs) are a major class of plant metabolites with many pharmacological benefits. Sacred lotus (Nelumbo nucifera) is an ancient aquatic plant of medicinal value because of antiviral and immunomodulatory activities linked to its constituent BIAs. Although more than 30 BIAs belonging to the 1-benzylisoquinoline, aporphine, and bisbenzylisoquinoline structural subclasses and displaying a predominant R-enantiomeric conformation have been isolated from N. nucifera, its BIA biosynthetic genes and enzymes remain unknown. Herein, we report the isolation and biochemical characterization of two O-methyltransferases (OMTs) involved in BIA biosynthesis in sacred lotus. Five homologous genes, designated NnOMT1-5 and encoding polypeptides sharing >40% amino acid sequence identity, were expressed in Escherichia coli Functional characterization of the purified recombinant proteins revealed that NnOMT1 is a regiospecific 1-benzylisoquinoline 6-O-methyltransferase (6OMT) accepting both R- and S-substrates, whereas NnOMT5 is mainly a 7-O-methyltransferase (7OMT), with relatively minor 6OMT activity and a strong stereospecific preference for S-enantiomers. Available aporphines were not accepted as substrates by either enzyme, suggesting that O-methylation precedes BIA formation from 1-benzylisoquinoline intermediates. Km values for NnOMT1 and NnOMT5 were 20 and 13 μm for (R,S)-norcoclaurine and (S)-N-methylcoclaurine, respectively, similar to those for OMTs from other BIA-producing plants. Organ-based correlations of alkaloid content, OMT activity in crude extracts, and OMT gene expression supported physiological roles for NnOMT1 and NnOMT5 in BIA metabolism, occurring primarily in young leaves and embryos of sacred lotus. In summary, our work identifies two OMTs involved in BIA metabolism in the medicinal plant N. nucifera.

Comprehensive Analysis and Functional Studies of WRKY Transcription Factors in Nelumbo nucifera

The WRKY family is one of the largest transcription factor (TF) families in plants and plays central roles in modulating plant stress responses and developmental processes, as well as secondary metabolic regulations. Lotus (Nelumbo nucifera) is an aquatic crop that has significant food, ornamental and pharmacological values. Here, we performed an overview analysis of WRKY TF family members in lotus, and studied their functions in environmental adaptation and regulation of lotus benzylisoquinoline alkaloid (BIA) biosynthesis. A total of 65 WRKY genes were identified in the lotus genome and they were well clustered in a similar pattern with their Arabidopsis homologs in seven groups (designated I, IIa-IIe, and III), although no lotus WRKY was clustered in the group IIIa. Most lotus WRKYs were functionally paired, which was attributed to the recently occurred whole genome duplication in lotus. In addition, lotus WRKYs were regulated dramatically by salicilic acid (SA), jasmonic acid (JA), and submergence treatments, and two lotus WRKYs, NnWRKY40a and NnWRKY40b, were significantly induced by JA and promoted lotus BIA biosynthesis through activating BIA biosynthetic genes. The investigation of WRKY TFs for this basal eudicot reveals new insights into the evolution of the WRKY family, and provides fundamental information for their functional studies and lotus breeding.

The Latest Studies on Lotus (Nelumbo nucifera)-an Emerging Horticultural Model Plant

Lotus (Nelumbo nucifera) is a perennial aquatic basal eudicot belonging to a small family Nelumbonaceace, which contains only one genus with two species. It is an important horticultural plant, with its uses ranging from ornamental, nutritional to medicinal values, and has been widely used, especially in Southeast Asia. Recently, the lotus obtained a lot of attention from the scientific community. An increasing number of research papers focusing on it have been published, which have shed light on the mysteries of this species. Here, we comprehensively reviewed the latest advancement of studies on the lotus, including phylogeny, genomics and the molecular mechanisms underlying its unique properties, its economic important traits, and so on. Meanwhile, current limitations in the research of the lotus were addressed, and the potential prospective were proposed as well. We believe that the lotus will be an important model plant in horticulture with the generation of germplasm suitable for laboratory operation and the establishment of a regeneration and transformation system.

Benzylisoquinoline Alkaloids Biosynthesis in Sacred Lotus

Sacred lotus (Nelumbo nucifera Gaertn.) is an ancient aquatic plant used throughout Asia for its nutritional and medicinal properties. Benzylisoquinoline alkaloids (BIAs), mostly within the aporphine and bisbenzylisoquinoline structural categories, are among the main bioactive constituents in the plant. The alkaloids of sacred lotus exhibit promising anti-cancer, anti-arrhythmic, anti-HIV, and anti-malarial properties. Despite their pharmacological significance, BIA metabolism in this non-model plant has not been extensively investigated. In this review, we examine the diversity of BIAs in sacred lotus, with an emphasis on the distinctive stereochemistry of alkaloids found in this species. Additionally, we discuss our current understanding of the biosynthetic genes and enzymes involved in the formation of 1-benzylisoquinoline, aporphine, and bisbenzylisoquinoline alkaloids in the plant. We conclude that a comprehensive functional characterization of alkaloid biosynthetic enzymes using both in vitro and in vivo methods is required to advance our limited knowledge of BIA metabolism in the sacred lotus.