The CYP80A and CYP80G Are Involved in the Biosynthesis of Benzylisoquinoline Alkaloids in the Sacred Lotus (Nelumbo nucifera)

Bisbenzylisoquinoline and aporphine alkaloids are the two main pharmacological compounds in the ancient sacred lotus (Nelumbo nucifera). The biosynthesis of bisbenzylisoquinoline and aporphine alkaloids has attracted extensive attention because bisbenzylisoquinoline alkaloids have been reported as potential therapeutic agents for COVID-19. Our study showed that NnCYP80A can catalyze C-O coupling in both (R)-N-methylcoclaurine and (S)-N-methylcoclaurine to produce bisbenzylisoquinoline alkaloids with three different linkages. In addition, NnCYP80G catalyzed C-C coupling in aporphine alkaloids with extensive substrate selectivity, specifically using (R)-N-methylcoclaurine, (S)-N-methylcoclaurine, coclaurine and reticuline as substrates, but the synthesis of C-ring alkaloids without hydroxyl groups in the lotus remains to be elucidated. The key residues of NnCYP80G were also studied using the 3D structure of the protein predicted using Alphafold 2, and six key amino acids (G39, G69, A211, P288, R425 and C427) were identified. The R425A mutation significantly decreased the catalysis of (R)-N-methylcoclaurine and coclaurine inactivation, which might play important role in the biosynthesis of alkaloids with new configurations.

NnWOX1-1, NnWOX4-3, and NnWOX5-1 of lotus (Nelumbo nucifera Gaertn)promote root formation and enhance stress tolerance in transgenic Arabidopsis thaliana

BACKGROUND

Adventitious roots (ARs) represent an important organ system for water and nutrient uptake in lotus plants because of degeneration of the principal root. The WUSCHEL-related homeobox (WOX) gene regulates plant development and growth by affecting the expression of several other genes. In this study, three WOX genes, NnWOX1-1, NnWOX4-3, and NnWOX5-1, were isolated and their functions were assessed in Arabidopsis plants.

RESULTS

The full lengths of NnWOX1-1, NnWOX4-3, and NnWOX5-1 were 1038, 645, and 558 bp, encoding 362, 214, and 185 amino acid residues, respectively. Phylogenetic analysis classified NnWOX1-1 and NnWOX4-3 encoding proteins into one group, and NnWOX5-1 and MnWOX5 encoding proteins exhibited strong genetic relationships. The three genes were induced by sucrose and indoleacetic acid (IAA) and exhibited organ-specific expression characteristics. In addition to improving root growth and salt tolerance, NnWOX1-1 and NnWOX4-3 promoted stem development in transgenic Arabidopsis plants. A total of 751, 594, and 541 genes, including 19, 19, and 13 respective genes related to ethylene and IAA metabolism and responses, were enhanced in NnWOX1-1, NnWOX4-3, and NnWOX5-1 transgenic plants, respectively. Further analysis showed that ethylene production rates in transgenic plants increased, whereas IAA, peroxidase, and lignin content did not significantly change. Exogenous application of ethephon on lotus seedlings promoted AR formation and dramatically increased the fresh and dry weights of the plants.

CONCLUSIONS

NnWOX1-1, NnWOX4-3, and NnWOX5-1 influence root formation, stem development, and stress adaptation in transgenic Arabidopsis plants by affecting the transcription of multiple genes. Among these, changes in gene expression involving ethylene metabolism and responses likely critically affect the development of Arabidopsis plants. In addition, ethylene may represent an important factor affecting AR formation in lotus seedlings.

Transcription factor NnMYB5 controls petal color by regulating GLUTATHIONE S-TRANSFERASE2 in Nelumbo nucifera

Lotus (Nelumbo spp.) is an important aquatic ornamental genus in the family Nelumbonaceae comprising only 2 species: Nelumbo lutea with yellow flowers and Nelumbo nucifera with red or white flowers. The petal color variations between these 2 species have previously been associated with the potential activities of FLAVONOL SYNTHASE (FLS) and MYB5. However, the underlying genetic mechanisms of flower color divergence within the N. nucifera species remain unclear. Here, quantitative trait locus mapping led to the identification of MYB5, a candidate gene controlling petal color in N. nucifera. Genotyping of 213 natural lotus accessions revealed an 80 kb presence/absence variant (PAV) of the NnMYB5 gene that is associated with petal color variation. Transcriptome analysis, dual-luciferase, and yeast 1-hybrid assays showed that NnMYB5 could directly activate the anthocyanin transporter gene GLUTATHIONE S-TRANSFERASE2 (NnGST2). Heterologous expression of NnGST2 in Arabidopsis (Arabidopsis thaliana) and its overexpression in lotus petals induced anthocyanin accumulation. Deletion of the 80 kb PAV within NnMYB5 inactivated NnGST2 expression and blocked anthocyanin accumulation in white N. nucifera petals. In contrast, the anthocyanin deficiency of N. lutea occurred due to pseudogenized NlMYB5 alleles. Our results establish a regulatory link between NnMYB5 and NnGST2 in petal anthocyanin accumulation and demonstrate the independent mechanisms controlling flower coloration in Nelumbo.

Genome-Wide Analysis of WUSCHEL-Related Homeobox Gene Family in Sacred Lotus (Nelumbo nucifera)

WUSCHEL-related homeobox (WOX) is a plant-specific transcription factor (TF), which plays an essential role in the regulation of plant growth, development, and abiotic stress responses. However, little information is available on the specific roles of WOX TFs in sacred lotus (Nelumbo nucifera), which is a perennial aquatic plant with important edible, ornamental, and medicinal values. We identified 15 WOX TFs distributing on six chromosomes in the genome of N. nucifera. A total of 72 WOX genes from five species were divided into three clades and nine subclades based on the phylogenetic tree. NnWOXs in the same subclades had similar gene structures and conserved motifs. Cis-acting element analysis of the promoter regions of NnWOXs found many elements enriched in hormone induction, stress responses, and light responses, indicating their roles in growth and development. The Ka/Ks analysis showed that the WOX gene family had been intensely purified and selected in N. nucifera. The expression pattern analysis suggested that NnWOXs were involved in organ development and differentiation of N. nucifera. Furthermore, the protein-protein interaction analysis showed that NnWOXs might participate in the growth, development, and metabolic regulation of N. nucifera. Taken together, these findings laid a foundation for further analysis of NnWOX functions.

Combined analysis of microRNA and mRNA profiles provides insights into the pathogenic resistant mechanisms of the lotus rhizome rot

Lotus rhizome rot caused by Fusarium oxysporum is a common vascular fungal disease in plants that significantly impacts the yield. However, only a few studies have studied the mechanism of Nelumbo nucifera responding to lotus rhizome rot. Here, we investigated the pathogenic genes and miRNAs in lotus rhizome rot to uncover the pathogenic resistant mechanisms by transcriptome and small RNA sequencing of lotus roots after inoculation with Fusarium oxysporum. GO and KEGG functional enrichment analysis showed that differential miRNAs were mostly enriched in starch and sucrose metabolism, biosynthesis of secondary metabolites, glutathione metabolism, brassinosteroid biosynthesis and flavonoid biosynthesis pathways. Twenty-seven upregulated miRNAs, 19 downregulated miRNAs and their target genes were identified. Correlation analysis found that miRNAs negatively regulate target genes, which were also enriched in starch and sucrose metabolism and glutathione metabolism pathways. Their expression was measured by reverse transcription quantitative PCR (qRT-PCR), and the results were consistent with the transcriptome analysis, thus verifying the reliability of transcriptome data. We selected three miRNAs (miRNA858-y, miRNA171-z and a novel miRNA novel-m0005-5p) to test the relationship between miRNAs and their target genes. The activity of the GUS testing assay indicated that miRNA could decrease the GUS activity by inhibiting the expression of their target genes. Collectively, this study provides a comprehensive analysis of transcriptome and small RNA sequencing of lotus root after inoculation with Fusarium oxysporum, and we identified candidate miRNAs and their target genes for breeding strategies of Nelumbo nucifera.

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.

Pathway elucidation and microbial synthesis of proaporphine and bis-benzylisoquinoline alkaloids from sacred lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera) has been utilized as a food, medicine, and spiritual symbol for nearly 3000 years. The medicinal properties of lotus are largely attributed to its unique profile of benzylisoquinoline alkaloids (BIAs), which includes potential anti-cancer, anti-malarial and anti-arrhythmic compounds. BIA biosynthesis in sacred lotus differs markedly from that of opium poppy and other members of the Ranunculales, most notably in an abundance of BIAs possessing the (R)-stereochemical configuration and the absence of reticuline, a major branchpoint intermediate in most BIA producers. Owing to these unique metabolic features and the pharmacological potential of lotus, we set out to elucidate the BIA biosynthesis network in N. nucifera. Here we show that lotus CYP80G (NnCYP80G) and a superior ortholog from Peruvian nutmeg (Laurelia sempervirens; LsCYP80G) stereospecifically convert (R)-N-methylcoclaurine to the proaporphine alkaloid glaziovine, which is subsequently methylated to pronuciferine, the presumed precursor to nuciferine. While sacred lotus employs a dedicated (R)-route to aporphine alkaloids from (R)-norcoclaurine, we implemented an artificial stereochemical inversion approach to flip the stereochemistry of the core BIA pathway. Exploiting the unique substrate specificity of dehydroreticuline synthase from common poppy (Papaver rhoeas) and pairing it with dehydroreticuline reductase enabled de novo synthesis of (R)-N-methylcoclaurine from (S)-norcoclaurine and its subsequent conversion to pronuciferine. We leveraged our stereochemical inversion approach to also elucidate the role of NnCYP80A in sacred lotus metabolism, which we show catalyzes the stereospecific formation of the bis-BIA nelumboferine. Screening our collection of 66 plant O-methyltransferases enabled conversion of nelumboferine to liensinine, a potential anti-cancer bis-BIA from sacred lotus. Our work highlights the unique benzylisoquinoline metabolism of N. nucifera and enables the targeted overproduction of potential lotus pharmaceuticals using engineered microbial systems.

Evolution of isoform-level gene expression patterns across tissues during lotus species divergence

Both gene duplication and alternative splicing (AS) drive the functional diversity of gene products in plants, yet the relative contributions of the two key mechanisms to the evolution of gene function are largely unclear. Here, we studied AS in two closely related lotus plants, Nelumbo lutea and Nelumbo nucifera, and the outgroup Arabidopsis thaliana, for both single-copy and duplicated genes. We show that most splicing events evolved rapidly between orthologs and that the origin of lineage-specific splice variants or isoforms contributed to gene functional changes during species divergence within Nelumbo. Single-copy genes contain more isoforms, have more AS events conserved across species, and show more complex tissue-dependent expression patterns than their duplicated counterparts. This suggests that expression divergence through isoforms is a mechanism to extend the expression breadth of genes with low copy numbers. As compared to isoforms of local, small-scale duplicates, isoforms of whole-genome duplicates are less conserved and display a less conserved tissue bias, pointing towards their contribution to subfunctionalization. Through comparative analysis of isoform expression networks, we identified orthologous genes of which the expression of at least some of their isoforms displays a conserved tissue bias across species, indicating a strong selection pressure for maintaining a stable expression pattern of these isoforms. Overall, our study shows that both AS and gene duplication contributed to the diversity of gene function during the evolution of lotus.

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.

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.

Transcriptome Analyses Revealed the Key Metabolic Genes and Transcription Factors Involved in Terpenoid Biosynthesis in Sacred Lotus

As the largest group of structurally diverse metabolites, terpenoids are versatile natural compounds that act as metabolism mediators, plant volatiles, and ecological communicators. However, few terpenoid compounds have been identified in plant parts of sacred lotus (Nelumbo nucifera Gaertn.). To elucidate the molecular genetic basis of the terpene biosynthetic pathway, terpenes from different parts of the plant, including seeds (S), young leaves (YL), mature leaves (ML), white flowers (WF), yellow flowers (YF), and red flowers (RF), were identified by LC-MS/MS and the relative contents of the same terpenes in different parts were compared. The results indicate that all plant parts primarily consist of triterpenes, with only minor quantities of sesquiterpenes and diterpenes, and there were differences in the terpene content detected in different plant parts. To illustrate the biosynthesis of various terpenoids, RNA sequencing was performed to profile the transcriptomes of various plant parts, which generated a total of 126.95 GB clean data and assembled into 29,630 unigenes. Among these unigenes, 105 candidate unigenes are involved in the mevalonate (MVA) pathway, methyl-erythritol phosphate (MEP) pathway, terpenoid backbone biosynthesis pathway, and terpenoid synthases pathway. Moreover, the co-expression network between terpene synthase (TPS) and WRKY transcription factors provides new information for the terpene biosynthesis pathway.

Integrated mRNA and Small RNA Sequencing Reveals a microRNA Regulatory Network Associated with Starch Biosynthesis in Lotus (Nelumbo nucifera Gaertn.) Rhizomes

Internode starch biosynthesis is one of the most important traits in lotus rhizome because of its relation to crop productivity. Understanding the microRNA (miRNA) and mRNA expression profiles related to lotus internode starch biosynthesis would help develop molecular improvement strategies, but they are not yet well-investigated. To identify genes and miRNAs involved in internode starch biosynthesis, the cDNA and small RNA libraries of Z6-1, Z6-2, and Z6-3 were sequenced, and their expression were further studied. Through combined analyses of transcriptome data and small RNA sequencing data, a complex co-expression regulatory network was constructed, in which 20 miRNAs could modulate starch biosynthesis in different internodes by tuning the expression of 10 target genes. QRT-PCR analysis, transient co-expression experiment and dual luciferase assay comprehensively confirmed that NnumiR396a down-regulated the expression of NnSS2 and ultimately prevents the synthesis of amylopectin, and NnumiR396b down-regulated the expression of NnPGM2 and ultimately prevents the synthesis of total starch. Our results suggest that miRNAs play a critical role in starch biosynthesis in lotus rhizome, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing lotus rhizome.

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.

Color fading in lotus (Nelumbo nucifera) petals is manipulated both by anthocyanin biosynthesis reduction and active degradation

Flower color is a key trait that determines the ornamental quality of aquatic lotus (Nelumbo nucifera). Color fading significantly decreases the ornamental value of lotus flowers. However, the molecular mechanism underlying lotus petal discoloration remains largely unknown. Here, the anthocyanin content and global transcriptional profiling of lotus petals of cultivar 'Qiusanse' in four developmental stages were analyzed. Five anthocyanin components were detected, and the total anthocyanin content decreased as the petal color changed from red to nearly white. Moreover, the malondialdehyde (MDA) content and peroxidase (POD) activity increased during color fading. RNA-seq analysis revealed a total of 4,092 differentially expressed genes (DEGs) between petal developmental stages. Notably, oxidoreductase and hydrolase activity related genes were overrepresented in DEGs. The expression pattern of key anthocyanin biosynthesis genes including, CHS, F3H, ANS, UFGT, and transcription factor regulators, including MYBs, WRKYs and bHLHs were correlated with anthocyanin accumulation. Interestingly, DEGs associated with anthocyanin degradation and vacuolar pH regulation, including peroxidase, proton pumps regulators such as WRKY3 and MYB5-like, were significantly upregulated during the late stages of flowering. This study reveals for the first time the transcriptional dynamics during lotus petal discoloration. Our results suggest the involvement of anthocyanin biosynthesis repressors and degrading genes as well as pH regulators in controlling color fading of lotus petals. The study also provides valuable information and candidate genes for improving the lotus flower color.

Short-term exposure to silver nano-particles alters the physiology and induces stress-related gene expression in Nelumbo nucifera

Lotus (Nelumbo nucifera) was used as model plant in this study to explore its physiology and molecular response upon short-term exposure to silver nano-particles (AgNPs). Accumulation patterns demonstrated a potential uptake of AgNPs by roots and transport to the leaves as a likely key translocation route in lotus. AgNPs exposure was negatively correlated with lotus growth, including germination rate and petiole length in a concentration-dependent manner. Synthesis of chloroplast pigments in lotus leaves was enhanced by low AgNPs concentration, but were inhibited at high concentration. Hydrogen peroxide (H₂O₂) was detected in lotus leaves after AgNPs treatment. Proline accumulation in lotus leaves was induced with the increase in AgNPs concentration and exposure time. Antioxidant enzyme activities of superoxide dismutase (SOD), peroxidase (POD) as well as catalase (CAT) were enhanced after the first day of AgNPs exposure, but declined with increased exposure time, indicating a time-dependent toxicity of AgNPs. In addition, real-time PCR revealed that two detoxification-related genes, GSH1 and GST, could be activated on the first day of AgNPs exposure, but down-regulated with prolonged AgNPs treatment. Photosynthesis-related RbcS gene was up-regulated, however, no obvious difference in the expression of RbcL was observed after the first day of AgNPs exposure. Moreover, WRKY70a and WRKY70b transcription factors exhibited similar expression patterns, with the highest induction after a 5 mg/L AgNPs exposure on the first day, which decreased with prolonged exposure time. This study provides useful references for further evaluation of the toxic mechanism of AgNPs and their bio-effects on aquatic plants and ecosystems.

Lignin Synthesis, Affected by Sucrose in Lotus (Nelumbo nucifera) Seedlings, Was Involved in Regulation of Root Formation in the Arabidopsis thanliana

Adventitious roots (ARs) have an unmatched status in plant growth and metabolism due to the degeneration of primary roots in lotuses. In the present study, we sought to assess the effect of sucrose on ARs formation and observed that lignin synthesis was involved in ARs development. We found that the lignification degree of the ARs primordium was weaker in plants treated with 20 g/L sucrose than in 50 g/L sucrose treatment and control plants. The contents of lignin were lower in plants treated with 20 g/L sucrose and higher in plants treated with 50 g/L sucrose. The precursors of monomer lignin, including p-coumaric acid, caffeate, sinapinal aldehyde, and ferulic acid, were lower in the GL50 library than in the GL20 library. Further analysis revealed that the gene expression of these four metabolites had no novel difference in the GL50/GL20 libraries. However, a laccase17 gene (NnLAC17), involved in polymer lignin synthesis, had a higher expression in the GL50 library than in the GL20 library. Therefore, NnLAC17 was cloned and the overexpression of NnLAC17 was found to directly result in a decrease in the root number in transgenic Arabidopsis plants. These findings suggest that lignin synthesis is probably involved in ARs formation in lotus seedlings.

LAMP Assay for Distinguishing Fusarium oxysporum and Fusarium commune in Lotus (Nelumbo nucifera) Rhizomes

Yields of edible rhizome from cultivation of the perennial hydrophyte lotus (Nelumbo nucifera) can be severely reduced by rhizome rot disease caused by Fusarium species. There is a lack of rapid field-applicable methods for detection of these pathogens on lotus plants displaying symptoms of rhizome rot. Fusarium commune (91%) and Fusarium oxysporum (9%) were identified at different frequencies from lotus samples showing symptoms of rhizome rot. Because these two species can cause different severity of disease and their morphology is similar, molecular diagnostic-based methods to detect these two species were developed. Based on the comparison of the mitochondrial genome of the two species, three specific DNA loci targets were found. The designed primer sets for conventional PCR, quantitative PCR, and loop-mediated isothermal amplification (LAMP) precisely distinguished the above two species when isolated from lotus and other plants. The LAMP detection limits were 10 pg/μl and 1 pg/μl of total DNA for F. commune and F. oxysporum, respectively. We also carried out field-mimicked experiments on lotus seedlings and rhizomes (including inoculated samples and field-diseased samples), and the results indicated that the LAMP primer sets and the supporting portable methods are suitable for rapid diagnosis of the lotus disease in the field. The LAMP-based detection method will aid in the rapid identification of whether F. oxysporum or F. commune is infecting lotus plants with symptoms of rhizome rot and can facilitate efficient pesticide use and prevent disease spread through vegetative propagation of Fusarium-infected lotus rhizomes.

Turgor regulation defect 1 proteins play a conserved role in pollen tube reproductive innovation of the angiosperms

Sexual reproduction in angiosperms is siphonogamous, and the interaction between pollen tube and pistil is critical for successful fertilization. Our previous study demonstrated that mutation of the Arabidopsis turgor regulation defect 1 (TOD1) gene leads to reduced male fertility, a result of retarded pollen tube growth in the pistil. TOD1 encodes a Golgi-localized alkaline ceramidase, a key enzyme for the production of sphingosine-1-phosphate (S1P), which is involved in the regulation of turgor pressure in plant cells. However, whether TOD1s play a conserved role in the innovation of siphonogamy is largely unknown. In this study, we provide evidence that OsTOD1, which is similar to AtTOD1, is also preferentially expressed in rice pollen grains and pollen tubes. OsTOD1 knockout results in reduced pollen tube growth potential in rice pistil. Both the OsTOD1 genomic sequence with its own promoter and the coding sequence under the AtTOD1 promoter can partially rescue the attod1 mutant phenotype. Furthermore, TOD1s from other angiosperm species can partially rescue the attod1 mutant phenotype, while TOD1s from gymnosperm species are not able to complement the attod1 mutant phenotype. Our data suggest that TOD1 acts conservatively in angiosperms, and this opens up an opportunity to dissect the role of sphingolipids in pollen tube growth in angiosperms.

Biased allelic expression in tissues of F1 hybrids between tropical and temperate lotus (Nelumbo nuicfera)

The genome-wide allele-specific expression in F1 hybrids from the cross of tropical and temperate lotus unveils how cis-regulatory divergences affect genes in key pathways related to ecotypic divergence. Genetic variation, particularly cis-regulatory variation, plays a crucial role in phenotypic variation and adaptive evolution in plants. Temperate and tropical lotus, the two ecotypes of Nelumbo nucifera, show distinction in the degree of rhizome enlargement, which is associated with winter dormancy. To understand the roles of genome-wide cis-regulatory divergences on adaptive evolution of temperate and tropical lotus (Nelumbo nucifera), here we performed allele-specific expression (ASE) analyses on the tissues including flowers, leaves and rhizome from F1 hybrids of tropical and temperate lotus. For all investigated tissues in F1s, about 36% of genes showed ASE and about 3% of genes showed strong consistent ASE. Most of ASEs were biased towards the tropical parent in all surveyed samples, indicating that the tropical genome might be dominant over the temperate genome in gene expression of tissues from their F1 hybrids. We found that promoter sequences with similar allelic expression are more conserved than genes with significant or conditional ASE, suggesting the cis-regulatory sequence divergence underlie the allelic expression bias. We further uncovered biased genes being related to phenotypic differentiation between two lotus ecotypes, especially metabolic and phytohormone-related pathways in the rhizome. Overall, our study provides a global landscape of cis-regulatory variations between two lotus ecotypes and highlights their roles in rhizome growth variation for the climatic adaptation.

Enzymatic basis for stepwise C-glycosylation in the formation of flavonoid di-C-glycosides in sacred lotus (Nelumbo nucifera Gaertn.)

Lotus plumule, the embryo of the seed of the sacred lotus (Nelumbo nucifera), contains a high accumulation of secondary metabolites including flavonoids and possesses important pharmaceutical value. Flavonoid C-glycosides, which accumulate exclusively in lotus plumule, have attracted considerable attention in recent decades due to their unique chemical structure and special bioactivities. As well as mono-C-glycosides, lotus plumule also accumulates various kinds of di-C-glycosides by mechanisms which are as yet unclear. In this study we identified two C-glycosyltransferase (CGT) genes by mining sacred lotus genome data and provide in vitro and in planta evidence that these two enzymes (NnCGT1 and NnCGT2, also designated as UGT708N1 and UGT708N2, respectively) exhibit CGT activity. Recombinant UGT708N1 and UGT708N2 can C-glycosylate 2-hydroxyflavanones and 2-hydroxynaringenin C-glucoside, forming flavone mono-C-glycosides and di-C-glycosides, respectively, after dehydration. In addition, the above reactions were successfully catalysed by cell-free extracts from tobacco leaves transiently expressing NnCGT1 or NnCGT2. Finally, enzyme assays using cell-free extracts of lotus plumule suggested that flavone di-C-glycosides (vicenin-1, vicenin-3, schaftoside and isoschaftoside) are biosynthesized through sequentially C-glucosylating and C-arabinosylating/C-xylosylating 2-hydroxynaringenin. Taken together, our results provide novel insights into the biosynthesis of flavonoid di-C-glycosides by proposing a new biosynthetic pathway for flavone C-glycosides in N. nucifera and identifying a novel uridine diphosphate-glycosyltransferase (UGT708N2) that specifically catalyses the second glycsosylation, C-arabinosylating and C-xylosylating 2-hydroxynaringenin C-glucoside.

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 CONSTANS-LIKE gene of Nelumbo nucifera could promote potato tuberization

CONSTANS-LIKE 5 of Nelumbo nucifera is capable of promoting potato tuberization through CONSTANS-FLOWERING LOCUS T and gibberellin signaling pathways with a probable association with lotus rhizome enlargement. Lotus (Nelumbo nucifera) is an aquatic plant that is affiliated to the Nelumbonaceace family. It is widely used as an ornamental, vegetable, and medicinal herb with its rhizome being a popular vegetable. To explore the molecular mechanism underlying its rhizome enlargement, we conducted a systematic analysis on the CONSTANS-LIKE (COL) gene family, with the results, indicating that this gene plays a role in regulating potato tuber expansion. These analyses included phylogenetic relationships, gene structure, and expressional patterns of lotus COL family genes. Based on these analyses, NnCOL5 was selected for further study on its potential function in lotus rhizome formation. NnCOL5 was shown to be located in the nucleus, and its expression was positively associated with the enlargement of lotus rhizome. Besides, the overexpression of NnCOL5 in potato led to increased tuber weight and starch content under short-day conditions without changing the number of tubers. Further analysis suggested that the observed tuber changes might be mediated by affecting the expression of genes in CO-FT and GA signaling pathways. These results provide valuable insight in understanding the functions of COL gene as well as the enlargement of lotus rhizome.

The lotus NnFTIP1 and NnFT1 regulate flowering time in Arabidopsis

In higher plants, floral signals are mainly collected and transduced to FLOWERING LOCUS T (FT) in Arabidopsis and its orthologues. The movement of FT from leaves to the shoot apical meristem (SAM) is partially mediated by FT-INTERACTING PROTEIN1 (FTIP1). Although the functions of OsFTIP1 in rice and DOFTIP1 in orchid in FT transport have also been investigated, the FTIP1 homologue in lotus (Nelumbo nucifera Gaertn.), a type of horticultural plant with high economic and cultural value, has not been isolated, and the mechanism of NnFT1 transport has not been explored. Here, we revealed that NnFTIP1 mediates the transport of NnFT1 in ectopic transgenic lines in Arabidopsis. Overexpression of NnFTIP1 in the ftip1-1 background rescued the late flowering phenotype of ftip1-1, indicating that NnFTIP1 has a conserved function as FTIP1. NnFTIP1 and NnFT1 share similar tissue expression patterns and subcellular localization. NnFTIP1 and NnFT1 interact both in vitro and in vivo. In addition, NnFTIP1 affects NnFT1 transport from leaves to the SAM. Furthermore, we found that NnUOF8, a MYB-like transcription factor, directly regulates the expression of NnFTIP1. Our results suggest that the functions of FTIP1 and FT are conserved during evolution in flowering plants.

A bHLH gene NnTT8 of Nelumbo nucifera regulates anthocyanin biosynthesis

Lotus is an important aquatic ornamental plant, whose flower color is one of the key horticultural traits that determines its ornamental value. Previous studies revealed that anthocyanins largely determined the red color of lotus flower, which are also the main component that has beneficial effects on human health. However, the regulation mechanism of flower pigmentation in lotus flower remains unclear. In the present study, in order to further understand the regulatory mechanism underlying the anthocyanin biosynthesis, a bHLH gene NnTT8 was characterized to be phylogenetically close to AtTT8 and the bHLH proteins from other plant species that have been indicated to be involved in the positive regulation of anthocyanin biosynthesis. Complementation analysis in Arabidopsis tt8 mutant showed that NnTT8 could function similarly to AtTT8 in regulating anthocyanin and proanthocyanin biosynthesis. An MYB transcription factor capable of interacting with NnTT8 was also characterized from lotus. The identification of a bHLH transcription factor playing regulatory roles in anthocyanin biosynthesis is crucial, as it might help to obtain more in-depth insight into the coloration of lotus and help in breeding high anthocyanin content lotus variety that can be explored for lotus flower beverages.

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.

Genome-wide identification of MADS-box gene family in sacred lotus (Nelumbo nucifera) identifies a SEPALLATA homolog gene involved in floral development

BACKGROUND

Sacred lotus (Nelumbo nucifera) is a vital perennial aquatic ornamental plant. Its flower shape determines the horticultural and ornamental values. However, the mechanisms underlying lotus flower development are still elusive. MADS-box transcription factors are crucial in various features of plant development, especially in floral organogenesis and specification. It is still unknown how the MADS-box transcription factors regulate the floral organogenesis in lotus.

RESULTS

To obtain a comprehensive insight into the functions of MADS-box genes in sacred lotus flower development, we systematically characterized members of this gene family based on the available genome information. A total of 44 MADS-box genes were identified, of which 16 type I and 28 type II genes were categorized based on the phylogenetic analysis. Furthermore, the structure of MADS-box genes and their expressional patterns were also systematically analyzed. Additionally, subcellular localization analysis showed that they are mainly localized in the nucleus, of which a SEPALLATA3 (SEP3) homolog NnMADS14 was proven to be involved in the floral organogenesis.

CONCLUSION

These results provide some fundamental information about the MADS-box gene family and their functions, which might be helpful in not only understanding the mechanisms of floral organogenesis but also breeding of high ornamental value cultivars in lotus.

Gene expression profiling reveals the effects of light on adventitious root formation in lotus seedlings (Nelumbo nucifera Gaertn.)

BACKGROUND

Lotus is an aquatic horticultural crop that is widely cultivated in most regions of China and is used as an important off-season vegetable. The principal root of lotus is degenerated, and adventitious roots (ARs) are irreplaceable for plant growth. We found that no ARs formed under darkness and that exposure to high-intensity light significantly promoted the development of root primordia. Four differential expression libraries based on three light intensities were constructed to monitor metabolic changes, especially in indole-3-acetic acid (IAA) and sugar metabolism.

RESULTS

AR formation was significantly affected by light, and high light intensity accelerated AR development. Metabolic changes during AR formation under different light intensities were evaluated using gene expression profiling by high-throughput tag-sequencing. More than 2.2 × 104 genes were obtained in each library; the expression level of most genes was between 0.01 and 100 (FPKF value). Libraries constructed from plants grown under darkness (D/CK), under 5000 lx (E/CK), and under 20,000 lx (F/CK) contained 1739, 1683, and 1462 upregulated genes and 1533, 995, and 834 downregulated genes, respectively, when compared to those in the initial state (CK). Additionally, we found that 1454 and 478 genes had altered expression in a comparison of libraries D/CK and F/CK. Gene transcription between libraries D/F ranged from a 5-fold decrease to a 5-fold increase. Twenty differentially expressed genes (DEGs) were involved in the signal transduction pathway, 28 DEGs were related to the IAA response, and 35 DEGs were involved in sugar metabolism. We observed that the IAA content was enhanced after seed germination, even in darkness; this was responsible for AR formation. We also observed that sucrose could eliminate the negative effect of 150 μMol IAA during AR development.

CONCLUSIONS

AR formation was regulated by IAA, even in the dark, where induction and developmental processes could also be completed. In addition, 36 genes displayed altered expression in carbohydrate metabolism and ucrose metabolism was involved in AR development (expressed stage) according to gene expression and content change characteristics.

Comprehensive analysis of AGPase genes uncovers their potential roles in starch biosynthesis in lotus seed

BACKGROUND

Starch in the lotus seed contains a high proportion of amylose, which endows lotus seed a promising property in the development of hypoglycemic and low-glycemic index functional food. Currently, improving starch content is one of the major goals for seed-lotus breeding. ADP-glucose pyrophosphorylase (AGPase) plays an essential role in regulating starch biosynthesis in plants, but little is known about its characterization in lotus.

RESULTS

We describe the nutritional compositions of lotus seed among 30 varieties with starch as a major component. Comparative transcriptome analysis showed that AGPase genes were differentially expressed in two varieties (CA and JX) with significant different starch content. Seven putative AGPase genes were identified in the lotus genome (Nelumbo nucifera Gaertn.), which could be grouped into two subfamilies. Selective pressure analysis indicated that purifying selection acted as a vital force in the evolution of AGPase genes. Expression analysis revealed that lotus AGPase genes have varying expression patterns, with NnAGPL2a and NnAGPS1a as the most predominantly expressed, especially in seed and rhizome. NnAGPL2a and NnAGPS1a were co-expressed with a number of starch and sucrose metabolism pathway related genes, and their expressions were accompanied by increased AGPase activity and starch content in lotus seed.

CONCLUSIONS

Seven AGPase genes were characterized in lotus, with NnAGPL2a and NnAGPS1a, as the key genes involved in starch biosynthesis in lotus seed. These results considerably extend our understanding on lotus AGPase genes and provide theoretical basis for breeding new lotus varieties with high-starch content.

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.

Integrated Omics Analyses Identify Key Pathways Involved in Petiole Rigidity Formation in Sacred Lotus

Sacred lotus (Nelumbo nucifera Gaertn.) is a relic aquatic plant with two types of leaves, which have distinct rigidity of petioles. Here we assess the difference from anatomic structure to the expression of genes and proteins in two petioles types, and identify key pathways involved in petiole rigidity formation in sacred lotus. Anatomically, great variation between the petioles of floating and vertical leaves were observed. The number of collenchyma cells and thickness of xylem vessel cell wall was higher in the initial vertical leaves’ petiole (IVP) compared to the initial floating leaves’ petiole (IFP). Among quantified transcripts and proteins, 1021 and 401 transcripts presented 2-fold expression increment (named DEGs, genes differentially expressed between IFP and IVP) in IFP and IVP, 421 and 483 proteins exhibited 1.5-fold expression increment (named DEPs, proteins differentially expressed between IFP and IVP) in IFP and IVP, respectively. Gene function and pathway enrichment analysis displayed that DEGs and DEPs were significantly enriched in cell wall biosynthesis and lignin biosynthesis. In consistent with genes and proteins expressions in lignin biosynthesis, the contents of lignin monomers precursors were significantly different in IFP and IVP. These results enable us to understand lotus petioles rigidity formation better and provide valuable candidate genes information on further investigation.

Small RNA and Transcriptome Sequencing Reveals miRNA Regulation of Floral Thermogenesis in Nelumbo nucifera

The sacred lotus (Nelumbo nucifera Gaertn.) can produce heat autonomously and maintain a relatively stable floral chamber temperature for several days when blooming. Floral thermogenesis is critical for flower organ development and reproductive success. However, the regulatory role of microRNA (miRNA) underlying floral thermogenesis in N. nucifera remains unclear. To comprehensively understand the miRNA regulatory mechanism of thermogenesis, we performed small RNA sequencing and transcriptome sequencing on receptacles from five different developmental stages. In the present study, a total of 172 known miRNAs belonging to 39 miRNA families and 126 novel miRNAs were identified. Twenty-nine thermogenesis-related miRNAs and 3024 thermogenesis-related mRNAs were screened based on their expression patterns. Of those, seventeen differentially expressed miRNAs (DEMs) and 1765 differentially expressed genes (DEGs) had higher expression during thermogenic stages. The upregulated genes in the thermogenic stages were mainly associated with mitochondrial function, oxidoreductase activity, and the energy metabolism process. Further analysis showed that miR156_2, miR395a_5, miR481d, and miR319p may play an important role in heat-producing activity by regulating cellular respiration-related genes. This study provides comprehensive miRNA and mRNA expression profile of receptacle during thermogenesis in N. nucifera, which advances our understanding on the regulation of floral thermogenesis mediated by miRNA.

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.

Differences in flavonoid pathway metabolites and transcripts affect yellow petal colouration in the aquatic plant Nelumbo nucifera

BACKGROUND

The Asia lotus (Nelumbo nucifera Gaertn.) is an ornamental aquatic plant with high economic value. Flower colour is an important ornamental trait, with much of N. nucifera breeding focusing on its yellow flowers. To explore the yellow flower colouration mechanism in N. nucifera, we analysed its pigment constituents and content, as well as gene expression in the flavonoid pathway, in two N. nucifera cultivars.

RESULTS

We performed metabolomic and gene expression analyses in two N. nucifera cultivars with yellow and white flowers, Molinqiuse (MLQS) and Yeguangbei (YGB), respectively, at five stages of flower colouration. Based on phenotypic observation and metabolite analyses, the later stages of flower colouration (S3-S5) were determined to be key periods for differences between MLQS and YGB, with dihydroflavonols and flavonols differing significantly between cultivars. Dihydroquercetin, dihydrokaempferol, and isorhamnetin were significantly higher in MLQS than in YGB, whereas kaempferol was significantly higher in YGB. Most of the key homologous structural genes in the flavonoid pathway were significantly more active in MLQS than in YGB at stages S1-S4.

CONCLUSION

In this study, we performed the first analyses of primary and secondary N. nucifera metabolites during flower colouration, and found that isorhamnetin and kaempferol shunting resulted in petal colour differences between MLQS and YGB. Based on our data integration analyses of key enzyme expression in the putative flavonoid pathways of the two N. nucifera cultivars, NnFLS gene substrate specificity and differential expression of NnOMTs may be related to petal colour differences between MLQS and YGB. These results will contribute to determining the mechanism of yellow flower colouration in N. nucifera, and will improve yellow petal colour breeding in lotus species.

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.

Evolutionarily conserved function of the sacred lotus (Nelumbo nucifera Gaertn.) CER2-LIKE family in very-long-chain fatty acid elongation

Identification of NnCER2 and NnCER2-LIKE from Nelumbo nucifera, which are required for the very-long-chain fatty acid elongation, provides new evidence that CER2 proteins are evolutionarily conserved across the eudicots. CER2-LIKE family proteins have been described as core components of the fatty acid elongase complex in Arabidopsis, maize, and rice, having specific function in synthesis of the C₃₀ to C₃₄ fatty acyl-CoA precursors of cuticular waxes. Little is known about the functional conservation in this gene family across species. In this study, two CER2-LIKE family proteins, NnCER2 and NnCER2-LIKE, were characterized from sacred lotus (Nelumbo nucifera), which is an ancient basal eudicot. The transcriptional expression of NnCER2 and NnCER2-LIKE was found in floating leaf blades, emergent petioles and vertical leaves, petals, and anthers. The NnCER2 and NnCER2-LIKE proteins were localized to the endoplasmic reticulum and nucleus. Overexpressing NnCER2 and NnCER2-LIKE in Arabidopsis led to alteration of cuticle wax structure in inflorescence stems, and this was associated with elevated 30, 32, and 34 carbon length wax compounds, and their derivatives. The different substrate specificities of NnCER2 and NnCER2-LIKE were explored using co-expression with AtCER6 in yeast cells. These findings provide clear evidence that the function of CER2 family proteins in producing VLCFAs is highly conserved across the eudicots.

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

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

Whole genome re-sequencing reveals evolutionary patterns of sacred lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera or lotus) is an important aquatic plant in horticulture and ecosystems. As a foundation for exploring genomic variation and evolution among different germplasms, we re-sequenced 19 individuals from three cultivated temperate lotus subgroups (rhizome, seed and flower lotus), one wild temperate lotus subgroup (wild lotus), one tropical lotus group (Thai lotus) and an outgroup (Nelumbo lutea). Through genetic diversity and polymorphism analysis by non-missing SNP sites widely distributed in the whole genome, we confirmed that wild and Thai lotus exhibited greater differentiation with a higher genomic diversity compared to cultivated lotus. Rhizome lotus had the lowest genomic diversity and a closer relationship to wild lotus, whereas the genomes of seed and flower lotus were admixed. Genes in energy metabolism process and plant immunity evolved rapidly in lotus, reflecting local adaptation. We established that candidate genes in genomic regions with significant differentiation associated with temperate and tropical lotus divergence always exhibited highly divergent expression pattern. Together, this study comprehensive and credible interpretates important patterns of genetic diversity and relationships, gene evolution, and genomic signature from ecotypic differentiation of sacred lotus.

The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre-eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA-likes. Additionally, the MIRNAs participating in lotus pollen-pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.

The NnCenH3 protein and centromeric DNA sequence profiles of Nelumbo nucifera Gaertn. (sacred lotus) reveal the DNA structures and dynamics of centromeres in basal eudicots

Centromeres on eukaryotic chromosomes consist of large arrays of DNA repeats that undergo very rapid evolution. Nelumbo nucifera Gaertn. (sacred lotus) is a phylogenetic relict and an aquatic perennial basal eudicot. Studies concerning the centromeres of this basal eudicot species could provide ancient evolutionary perspectives. In this study, we characterized the centromeric marker protein NnCenH3 (sacred lotus centromere-specific histone H3 variant), and used a chromatin immunoprecipitation (ChIP)-based technique to recover the NnCenH3 nucleosome-associated sequences of sacred lotus. The properties of the centromere-binding protein and DNA sequences revealed notable divergence between sacred lotus and other flowering plants, including the following factors: (i) an NnCenH3 alternative splicing variant comprising only a partial centromere-targeting domain, (ii) active genes with low transcription levels in the NnCenH3 nucleosomal regions, and (iii) the prevalence of the Ty1/copia class of long terminal repeat (LTR) retrotransposons in the centromeres of sacred lotus chromosomes. In addition, the dynamic natures of the centromeric region showed that some of the centromeric repeat DNA sequences originated from telomeric repeats, and a pair of centromeres on the dicentric chromosome 1 was inactive in the metaphase cells of sacred lotus. Our characterization of the properties of centromeric DNA structure within the sacred lotus genome describes a centromeric profile in ancient basal eudicots and might provide evidence of the origins and evolution of centromeres. Furthermore, the identification of centromeric DNA sequences is of great significance for the assembly of the sacred lotus genome.

Small RNA and degradome profiling reveals miRNA regulation in the seed germination of ancient eudicot Nelumbo nucifera

BACKGROUND

MicroRNAs (miRNAs) play important roles in plant growth and development. MiRNAs and their targets have been widely studied in model plants, but limited knowledge is available concerning this small RNA population and their targets in sacred lotus (Nelumbo nucifera Gaertn.).

RESULTS

In this study, a total of 145 known miRNAs belonging to 47 families and 78 novel miRNAs were identified during seed germination using high-throughput small RNA sequencing. Furthermore, some miRNA families which have not yet been reported in monocot or eudicot species were detected in N. nucifera, indicating that these miRNAs was divergence from monocots and core eudicots during evolution. Using degradome sequencing, 2580 targets were detected for all the miRNAs. GO (Gene Ontology) and KEGG pathway analyses showed that many target genes enriched in "regulation of transcription" and involved in "carbohydrate", "amino acid and energy metabolism". Nine miRNAs and three corresponding targets of them were further validated by quantitative RT-PCR.

CONCLUSIONS

The results present here suggested that many miRNAs were involved in the regulation of seed germination of sacred lotus, providing a foundation for future studies of sacred lotus seed longevity. Comparative analysis of miRNAs from different plants also provided insight into the evolutionary gains and losses of miRNAs in plants.

Construction of a high-density, high-quality genetic map of cultivated lotus (Nelumbo nucifera) using next-generation sequencing

BACKGROUND

The sacred lotus (Nelumbo nucifera) is widely cultivated in China for its edible rhizomes and seeds. Traditional plant breeding methods have been used to breed cultivars with increased yields and quality of rhizomes and seeds with limited success. Currently, the available genetic maps and molecular markers in lotus are too limited to be useful for molecular genetics based breeding programs. However, the development of next-generation sequencing (NGS) technologies has enabled large-scale identification of single-nucleotide polymorphisms (SNPs) for genetic map construction. In this study, we constructed an SNP-based high-density genetic map for cultivated lotus using double digest restriction site-associated DNA sequencing (ddRADseq).

RESULTS

An F2 population of 96 individuals was derived from a cross between the rhizome lotus cultivar 'Juwuba' (male parent) and the seed lotus cultivar 'Mantianxing' (female parent). Genomic DNAs from this population were digested with the restriction enzymes EcoRI and MspI and then sequenced. In total, 133.65 Gb of raw data containing 1,088,935,610 pair-end reads were obtained. The coverage of reads on a reference genome was 7.2 % for the female parent, 6.56 % for the male parent, and 1.46 % for F2 individuals. From these reads, 10,753 valid SNP markers were used for genetic map construction. Finally, 791 bin markers (so-segregated adjacent SNPs treated as a bin marker), consisting of 8,971 SNP markers, were sorted into 8 linkage groups (LGs) that spanned 581.3 cM, with an average marker interval of 0.74 cM. A total of 809 genome sequence scaffolds, covering about 565.9 cM of the wild sacred lotus genome, were anchored on the genetic map, accounting for 70.6 % of the genome assembly.

CONCLUSIONS

This study reports the large-scale discovery of SNPs between cultivars of rhizome and seed lotus using a ddRADseq library combined with NGS. These SNPs have been used to construct the first high-density genetic map for cultivated lotus that can serve as a genomic reference and will facilitate genetic mapping of important traits in the parental cultivars.

Molecular and characterization of NnPPO cDNA from lotus (Nelumbo nucifera) in rhizome browning

The complete cDNA (NnPPO) of polyphenol oxidase in Nelumbo nucifera was successfully isolated, using Rapid amplification cDNA end (RACE) assays. The full-length cDNA of NnPPO was 2069 bp in size, containing a 1791 bp open reading frame coding 597 amino acids. The putative NnPPO possessed the conserved active sites and domains for PPO function. Phylogenetic analysis revealed that NnPPO shared high homology with PPO of high plants, and the homology modeling proved that NnPPO had the typical structure of PPO family. In order to characterize the role of NnPPO, Real-time PCR assay demonstrated that NnPPO mRNA was expressed in different tissues of N. nucifera including young leave, rhizome, flower, root and leafstalk, with the highest expression in rhizome. Patterns of NnPPO expression in rhizome illustrated its mRNA level was significantly elevated, which was consistent with the change of NnPPO activity during rhizome browning. Therefore, transcriptional activation of NnPPO was probably the main reason causing rhizome browning.

Identification and characterization of chalcone synthase cDNAs (NnCHS) from Nelumbo nucifera

Chalcone synthase (CHS) catalyzes the first committed step in flavonoids biosynthetic pathway. In this study, six full-length cDNAs (NnCHS) encoding CHS from Nelumbo nucifera were successfully isolated, using rapid amplification cDNA end (RACE) assay. The obtained cDNAs were 1426 bp in size, containing a 1167 bp open reading frame coding 389 amino acids. Exons-intron architecture of NnCHS gene was illustrated, consisting two exons inserted by a 426 bp intron. The putative NnCHS possessed all the conserved active sites for CHS function as well as the family signature. Phylogenetic analysis revealed that NnCHS shared high homology with CHS from high plants, and the homology-based structural modeling showed that NnCHS had the typical structure of CHS. Moreover, Real-time PCR assays demonstrated that NnCHS mRNAs were expressed in various tissues of N. nucifera, with the highest expression in red flower and lowest level in the leaves. Moreover, patterns of NnCHS expression illustrated short-time wounding or low temperature significantly induced the up-regulation of NnCHS mRNA.

Genome-Wide Identification of SSR and SNP Markers Based on Whole-Genome Re-Sequencing of a Thailand Wild Sacred Lotus (Nelumbo nucifera)

Genomic resources such as single nucleotide polymorphism (SNPs), insertions and deletions (InDels) and SSRs (simple sequence repeats) are essential for crop improvement and better utilization in genetic breeding. However, the resources for the sacred lotus (Nelumbo nucifera Gaertn.) are still limited. In the present study, to dissect large-scale genomic molecular marker resources for sacred lotus, we re-sequenced a Thailand sacred lotus cultivar ‘Chiang Mai wild lotus’ and compared with the reported lotus genome ‘Middle lake wild lotus’. A total of 3,180,059 SNPs, 328, 251 InDels and 14,191 SVs were found between the two genomes. The functional impact analyses of these SNPs indicated that they may be involved in metabolic processes, binding, catalytic activity, etc. Mining the genome sequences for SSRs showed that 191,657 SSRs were identified with a frequency of one SSR per 4.23 kb and 103,656 SSR primer pairs were designed. Furthermore, 14, 502 EST-SSRs were also indentified using the available RNA-seq data in the NCBI. A subset of 150 SSRs (genomic and EST-SSRs) was randomly selected for validation and genetic diversity analysis. The genotypes could be easily distinguished using these SSR markers and the ‘Chiang Mai wild lotus’ was obviously differentiated from the other Chinese accessions. This study provides considerable amounts of genomic resources and markers for the quantitative trait locus (QTL) identification and molecular selection of the species, which could have a potential role in various applications in sacred lotus breeding.

Molecular Cloning and Expression Analysis of a Catalase Gene (NnCAT) from Nelumbo nucifera

Rapid amplification cDNA end (RACE) assay was established to achieve the complete cDNA sequence of a catalase gene (NnCAT) from Nelumbo nucifera. The obtained full-length cDNA was 1666 bp in size and contained a 1476-bp open reading frame. The 3D structural model of NnCAT was constructed by homology modeling. The putative NnCAT possessed all the main characteristic amino acid residues and motifs of catalase (CAT) protein family, and the phylogenetic analysis revealed that NnCAT grouped together with high plants. Moreover, recombinant NnCAT showed the CAT activity (758 U/mg) at room temperature, holding high activity during temperature range of 20-50 °C, then the optimal pH of recombinant protein was assessed from pH 4 to pH 11. Additionally, real-time PCR assay demonstrated that NnCAT mRNA was expressed in various tissues of N. nucifera, with the highest expression in young leaf and lowest level in the root, and mRNA level of NnCAT was significantly augmented in response to short-time mechanical wounding. Different expression pattern of NnCAT gene suggested that NnCAT probably played a defensive role in the initial stages of oxidative stress, regulating the level of reactive oxygen species (ROS) by extracellular stimuli such as short-time mechanical wounding.

Proteomic and Epigenetic Analyses of Lotus (Nelumbo nucifera) Petals Between Red and White cultivars

Lotus is a vital aquatic ornamental plant with different flower colors. To explore the flower coloration mechanism in lotus, the constituents and contents of pigments in two lotus cultivars with red and white flowers were analyzed. Although flavones and flavonols were detected in both cultivars, anthocyanins could only be detected in the red cultivar. A comparative proteomics analysis on the flower petals between these two cultivars was conducted. A total of 88 differentially expressed proteins were identified with 36 more abundant and 52 less abundant in the red than in the white cultivar. Among them, four enzymes involved in the anthocyanin pathway were identified, i.e. flavanone 3-hydroxylase, anthocyanidin synthase, anthocyanidin 3-O-glucosyltransferase and glutathione S-transferase. Analysis of the expression patterns of anthocyanin biosynthetic genes indicated that the anthocyanindin synthase (ANS) gene might be the critical gene determining anthocyanin biosynthesis and accumulation in lotus flower. Further analysis showed that different methylation intensities on the promoter sequence of the ANS gene might result in the different flower coloration in the red and white cultivar. This study provides new insights into the mechanism of flower coloration in lotus, and may be helpful in its breeding and germplasm enhancement.

Cloning and expression of an APETALA1-like gene from Nelumbo nucifera

The objective of this study was to clone the full-length cDNA of the APETALA1 (AP1) gene from lotus and analyze its sequence and expression pattern. The full-length cDNA sequence of the NnAP1 gene was amplified from the petals of Nelumbo nucifera 'Hongxia' using RT-PCR and rapid amplification of cDNA ends. Bioinformatic methods were used to analyze the sequence characteristics of the gene. Quantitative real-time PCR methods were used to investigate the expression pattern of NnAP1 in various organs and during different developmental stages. The cloned full-length NnAP1 cDNA (GenBank accession No. KF361315) was 902 bp, containing a 795-bp open reading frame encoding 264 amino acids with a relative molecular mass of 30,288.4 and an isoelectric point of 9.13. NnAP1 had a MADS-box domain and a K-box domain, which is typical of the SQUA/AP1 gene family. A protein sequence identity search showed that NnAP1 was 75-96% similar to other plant AP1s. Phylogenetic tree analysis indicated that NnAP1 was very closely related to AP1 of Glycine max, suggesting that they shared the same protein ancestor. Quantitative real-time PCR analysis showed that NnAP1 was expressed in various organs during different developmental stages; it had the highest expression in blooming flowers and had trace expression in the young vegetative and flower senescence stages. Our analysis suggests that NnAP1 plays an important role in controlling floral meristem identity and floral organ formation.

RNA-Seq Uncovers SNPs and Alternative Splicing Events in Asian Lotus (Nelumbo nucifera)

RNA-Seq is an efficient way to comprehensively identify single nucleotide polymorphisms (SNPs) and alternative splicing (AS) events from the expressed genes. In this study, we conducted transcriptome sequencing of four Asian lotus (Nelumbo nucifera) cultivars using Illumina HiSeq2000 platform to identify SNPs and AS events in lotus. A total of 505 million pair-end RNA-Seq reads were generated from four cultivars, of which 86% were mapped to the lotus reference genome. Using the four sets of data together, a total of 357,689 putative SNPs were identified with an average density of one SNP per 2.2 kb. These SNPs were located in 1,253 scaffolds and 15,016 expressed genes. A/G and C/T were the two major types of SNPs in the Asian lotus transcriptome. In parallel, a total of 177,540 AS events were detected in the four cultivars and were distributed in 64% of the expressed genes of lotus. The predominant type of AS events was alternative 5’ first exon, which accounted for 41.2% of all the observed AS events, and exon skipping only accounted for 4.3% of all AS. Gene Ontology analysis was conducted to analyze the function of the genes containing SNPs and AS events. Validation of selected SNPs and AS events revealed that 74% of SNPs and 80% of AS events were reliable, which indicates that RNA-Seq is an efficient approach to uncover gene-associated SNPs and AS events. A large number of SNPs and AS events identified in our study will facilitate further genetic and functional genomics research in lotus.

Genome-wide discovery and analysis of phased small interfering RNAs in Chinese sacred lotus

Phased small interfering RNA (phasiRNA) generating loci (briefly as PHAS) in plants are a novel class of genes that are normally regulated by microRNAs (miRNAs). Similar to miRNAs, phasiRNAs encoded by PHAS play important regulatory roles by targeting protein coding transcripts in plant species. We performed a genome-wide discovery of PHAS loci in Chinese sacred lotus and identified a total of 106 PHAS loci. Of these, 47 loci generate 21 nucleotide (nt) phasiRNAs and 59 loci generate 24 nt phasiRNAs, respectively. We have also identified a new putative TAS3 and a putative TAS4 loci in the lotus genome. Our results show that some of the nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance proteins and MYB transcription factors potentially generate phasiRNAs. Furthermore, our results suggest that some large subunit (LSU) rRNAs can derive putative phasiRNAs, which is potentially resulted from crosstalk between small RNA biogenesis pathways that are employed to process rRNAs and PHAS loci, respectively. Some of the identified phasiRNAs have putative trans-targets with less than 4 mismatches, suggesting that the identified PHAS are involved in many different pathways. Finally, the discovery of 24 nt PHAS in lotus suggests that there are 24 nt PHAS in dicots.