Background: Sacred lotus is a basal eudicot with agricultural, medicinal, cultural and religious importance. It was domesticated in Asia about 7,000 years ago, and cultivated for its rhizomes and seeds as a food crop. It is particularly noted for its 1,300-year seed longevity and exceptional water repellency, known as the lotus effect. The latter property is due to the nanoscopic closely packed protuberances of its self-cleaning leaf surface, which have been adapted for the manufacture of a self-cleaning industrial paint, Lotusan. Results: The genome of the China Antique variety of the sacred lotus was sequenced with Illumina and 454 technologies, at respective depths of 101× and 5.2×. The final assembly has a contig N50 of 38.8 kbp and a scaffold N50 of 3.4 Mbp, and covers 86.5% of the estimated 929 Mbp total genome size. The genome notably lacks the paleo-triplication observed in other eudicots, but reveals a lineage-specific duplication. The genome has evidence of slow evolution, with a 30% slower nucleotide mutation rate than observed in grape. Comparisons of the available sequenced genomes suggest a minimum gene set for vascular plants of 4,223 genes. Strikingly, the sacred lotus has 16 COG2132 multi-copper oxidase family proteins with root-specific expression; these are involved in root meristem phosphate starvation, reflecting adaptation to limited nutrient availability in an aquatic environment. Conclusions: The slow nucleotide substitution rate makes the sacred lotus a better resource than the current standard, grape, for reconstructing the pan-eudicot genome, and should therefore accelerate comparative analysis between eudicots and monocots.

Background: The genus Nelumbo Adans. comprises two living species, N. nucifera Gaertan. (Asian lotus) and N. lutea Pers. (American lotus). A genetic linkage map is an essential resource for plant genetic studies and crop improvement but has not been generated for Nelumbo. We aimed to develop genomic simple sequence repeat (SSR) markers from the genome sequence and construct two genetic maps for Nelumbo to assist genome assembly and integration of a genetic map with the genome sequence. Results: A total of 86,089 SSR motifs were identified from the genome sequences. Di- and tri-nucleotide repeat motifs were the most abundant, and accounted for 60.73% and 31.66% of all SSRs, respectively. AG/GA repeats constituted 51.17% of dinucleotide repeat motifs, followed by AT/TA (44.29%). Of 500 SSR primers tested, 386 (77.20%) produced scorable alleles with an average of 2.59 per primer, and 185 (37.00%) showed polymorphism among two parental genotypes, N. nucifera ‘Chinese Antique’ and N. lutea ‘AL1’, and six progenies of their F₁ population. The normally segregating markers, which comprised 268 newly developed SSRs, 37 previously published SSRs and 53 sequence-related amplified polymorphism markers, were used for genetic map construction. The map for Asian lotus was 365.67 cM with 47 markers distributed in seven linkage groups. The map for American lotus was 524.51 cM, and contained 177 markers distributed in 11 genetic linkage groups. The number of markers per linkage group ranged from three to 34 with an average genetic distance of 3.97 cM between adjacent markers. Moreover, 171 SSR markers contained in linkage groups were anchored to 97 genomic DNA sequence contigs of ‘Chinese Antique’. The 97 contigs were merged into 60 scaffolds. Conclusion: Genetic mapping of SSR markers derived from sequenced contigs in Nelumbo enabled the associated contigs to be anchored in the linkage map and facilitated assembly of the genome sequences of ‘Chinese Antique’. The present study reports the first construction of genetic linkage maps for Nelumbo, which can serve as reference linkage maps to accelerate characterization germplasm, genetic mapping for traits of economic interest, and molecular breeding with marker-assisted selection.

Wild rice (Zizania latifolia) is an important genetic resource and it has important ecological functions in aquatic ecosystems as well. Here, we used landscape genetics to investigate how the landscape features, including rivers, mountains and habitat fragmentation, affect the genetic connectivity or create dispersal barriers for Z. latifolia. In this report, seventeen populations from the Sanjiang Plain and its surrounding areas were genotyped by ten microsatellite markers. A high genetic differentiation and genetic discontinuity were found among populations within each river investigated, suggesting that the rivers be not acting as corridors for dispersal. Meanwhile, genetic discontinuity was detected from different sides of the Lesser Khingan and its branch, the Qinghei Mountains, demonstrating that gene flow was blocked by the complex topography of mountains. The results that historical gene flow was much higher than the contemporary gene flow might infer that recent habitat fragmentation resulted in decreased gene flow. For all sampled Z. latifolia populations, the result of low level of genetic variation (na = 2.3, HE = 0.328) and high genetic divergence (FST = 0.405, Dest = 0.414 and Øpt = 0.424) was consistent with the decreased gene flow, an inbreeding system and repeated genetic bottlenecks. A conservation strategy for protecting as many populations as possible to maximize genomic representation of the species is proposed. In addition, dredging the watercourses should be carried out to improve habitat stability and to facilitate connectivity among populations.

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.

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.

The predominant alkaloids in lotus leaves are aporphine alkaloids. These are the most important active components and have many pharmacological properties, but little is known about their biosynthesis. We used digital gene expression technology to identify differentially-expressed genes (DEGs) between two lotus cultivars with different alkaloid contents at four leaf development stages. We also predicted potential genes involved in aporphine alkaloid biosynthesis by weighted gene co-expression network analysis (WGCNA). Approximately 335 billion nucleotides were generated; and 94% of which were aligned against the reference genome. Of 22 thousand expressed genes, 19,000 were differentially expressed between the two cultivars at the four stages. Gene Ontology (GO) enrichment analysis revealed that catalytic activity and oxidoreductase activity were enriched significantly in most pairwise comparisons. In Kyoto Encyclopaedia of Genes and Genomes (KEGG) analysis, dozens of DEGs were assigned to the categories of biosynthesis of secondary metabolites, isoquinoline alkaloid biosynthesis, and flavonoid biosynthesis. The genes encoding norcoclaurine synthase (NCS), norcoclaurine 6-O-methyltransferase (6OMT), coclaurine N-methyltransferase (CNMT), N-methylcoclaurine 3′-hydroxylase (NMCH), and 3′-hydroxy-N-methylcoclaurine 4′-O-methyltransferase (4′OMT) in the common pathways of benzylisoquinoline alkaloid biosynthesis and the ones encoding corytuberine synthase (CTS) in aporphine alkaloid biosynthetic pathway, which have been characterized in other plants, were identified in lotus. These genes had positive effects on alkaloid content, albeit with phenotypic lag. The WGCNA of DEGs revealed that one network module was associated with the dynamic change of alkaloid content. Eleven genes encoding proteins with methyltransferase, oxidoreductase and CYP450 activities were identified. These were surmised to be genes involved in aporphine alkaloid biosynthesis. This transcriptomic database provides new directions for future studies on clarifying the aporphine alkaloid pathway.