Proteomic analysis showing the signaling pathways involved in the rhizome enlargement process in Nelumbo nucifera

Spatiotemporal transcriptomic atlas of rhizome formation in Oryza longistaminata

Rhizomes are modified stems that grow underground and produce new individuals genetically identical to the mother plant. Recently, a breakthrough has been made in efforts to convert annual grains into perennial ones by utilizing wild rhizomatous species as donors, yet the developmental biology of this organ is rarely studied. Oryza longistaminata, a wild rice species featuring strong rhizomes, provides a valuable model for exploration of rhizome development. Here, we first assembled a double-haplotype genome of O. longistaminata, which displays a 48-fold improvement in contiguity compared to the previously published assembly. Furthermore, spatiotemporal transcriptomics was performed to obtain the expression profiles of different tissues in O. longistaminata rhizomes and tillers. Two spatially reciprocal cell clusters, the vascular bundle 2 cluster and the parenchyma 2 cluster, were determined to be the primary distinctions between the rhizomes and tillers. We also captured meristem initiation cells in the sunken area of parenchyma located at the base of internodes, which is the starting point for rhizome initiation. Trajectory analysis further indicated that the rhizome is regenerated through de novo generation. Collectively, these analyses revealed a spatiotemporal transcriptional transition underlying the rhizome initiation, providing a valuable resource for future perennial crop breeding.

Genome-Wide Identification and Co-Expression Networks of WOX Gene Family in Nelumbo nucifera

WUSCHEL-related homeobox (WOX) genes are a class of plant-specific transcription factors, regulating the development of multiple tissues. However, the genomic characterizations and expression patterns of WOX genes have not been analyzed in lotus. In this study, 15 NnWOX genes were identified based on the well-annotated reference genome of lotus. According to the phylogenetic analysis, the NnWOX genes were clustered into three clades, i.e., ancient clade, intermediate clade, and WUS clade. Except for the conserved homeobox motif, we further found specific motifs of NnWOX genes in different clades and divergence gene structures, suggesting their distinct functions. In addition, two NnWOX genes in the ancient clade have conserved expression patterns and other NnWOX genes exhibit different expression patterns in lotus tissues, suggesting a low level of functional redundancy in lotus WOX genes. Furthermore, we constructed the gene co-expression networks for each NnWOX gene. Based on weighted gene co-expression network analysis (WGCNA), ten NnWOX genes and their co-expressed genes were assigned to the modules that were significantly related to the cotyledon and seed coat. We further performed RT-qPCR experiments, validating the expression levels of ten NnWOX genes in the co-expression networks. Our study reveals comprehensive genomic features of NnWOX genes in lotus, providing a solid basis for further function studies.

Screening of Early Flowering Lotus (Nelumbo nucifera Gaertn.) Cultivars and Effects of Different Cultivars on Flowering Period

Flowering time is an important trait that determines the breeding process of ornamental plants. The flowering period of lotus (Nelumbo nucifera Gaertn.) is mainly concentrated in June-August. During this period, the weather is hot and there are few tourists, which made many lotus scenic spots difficult to operate. People have a strong demand for early flowering lotus cultivars. In this paper, 30 lotus cultivars with high ornamental value were selected as materials and their phenological periods were observed for two consecutive years in 2019 and 2020. A number of cultivars with early flowering potential and stable flowering periods, such as 'Fenyanzi', 'Chengshanqiuyue', 'Xianghumingyue' and 'Wuzhilian', were screened by K-Means clustering method. The relationship between accumulated temperature and flowering time of 19 lotus cultivars at different growth stages was analyzed. It was found that lotus cultivars with early flowering traits could adapt well to the changes of early environmental temperature and were not affected by low temperature. On the other hand, by analyzing the relationship between different traits and flowering time of three typical cultivars, such as rhizome weight, phenological period, etc., it shows that the nutrient content of the rhizome and the early morphology of plants will affect the flowering time. These results provide a reference for the formation of a systematic lotus early flowering cultivar breeding mechanism and the establishment of a perfect flowering regulation technology system, which can further improve the ornamental value of lotus and promote industrial development.

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.

Expression Patterns of Key Genes in the Photoperiod and Vernalization Flowering Pathways in Lilium longiflorum with Different Bulb Sizes

Lilium longiflorum is a wild Lilium, and its flowering transition requires a long period of cold exposure to meet the demand of vernalization. The responses of different sized bulbs to cold exposure and photoperiod are different, and the floral transition pathways of small and large bulbs are different. In this study, small and large bulbs were placed in cold storage for different weeks and then cultured at a constant ambient temperature of 25 °C under long day (LD) and short day (SD) conditions. Then, the flowering characteristics and expression patterns of key genes related to the vernalization and photoperiod pathways in different groups were calculated and analyzed. The results showed that the floral transition of Lilium longiflorum was influenced by both vernalization and photoperiod, that vernalization and LD conditions can significantly improve the flowering rate of Lilium longiflorum, and that the time from planting to visible flowering buds’ appearance was decreased. The flowering time and rate of large bulbs were greatly influenced by cold exposure, and the vernalization pathway acted more actively at the floral transition stage. The floral transition of small bulbs was affected more by the photoperiod pathway. Moreover, it was speculated that cold exposure may promote greater sensitivity of the small bulbs to LD conditions. In addition, the expression of LlVRN1, LlFKF1, LlGI, LlCO5, LlCO7, LlCO16, LlFT1, LlFT3 and LlSOC1 was high during the process of floral transition, and LlCO13, LlCO14 and LlCO15 were highly expressed in the vegetative stage. The expression of LlCO13 and LlCO14 was different under different lighting conditions, and the flowering induction function of LlCO9 and LlFT3 was related to vernalization. Moreover, LlFKF1, LlGI, LlCO5, LlCO16, LlSOC1 and LlFT2 were involved in the entire growth process of plants, while LlCO6, LlCO16 and LlFT1 are involved in the differentiation and formation of small bulblets of plants after the inflorescence stage, and this process is also closely related to LD conditions. This study has great significance for understanding the molecular mechanisms of the vernalization and photoperiod flowering pathways of Lilium longiflorum.

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.

Studies on Lotus Genomics and the Contribution to Its Breeding

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

Expression rewiring and methylation of non-coding RNAs involved in rhizome phenotypic variations of lotus ecotypes

Non-coding RNAs (ncRNAs), including miRNAs, lncRNAs, and circRNAs, emerge as crucial components for gene regulation. Nelumbo nucifera (lotus), a horticulturally important plant, differentiates into a temperate ecotype of enlarged rhizomes and a tropical ecotype of thin rhizomes. Nevertheless, whether and how ncRNAs can be rewired in expression and differentially methylated contributing to adaptive divergence of this storage organ in lotus ecotypes is unclear. Herein, we study the expression behaviors and DNA methylation patterns of ncRNAs in temperate and tropical lotus rhizomes. By whole transcriptome sequencing, we found both mRNAs and lncRNAs have divergent expression patterns between ecotypes, whereas miRNAs and circRNAs tended to be accession-specific or noisier in expression. The differentially expressed ncRNAs are involved in phenotypic differentiation of lotus rhizome between ecotypes, as the genes that interacted with them in the competing endogenous RNA network are enriched in functions including carbohydrate metabolism and plant hormone signaling, being critical to rhizome enlargement. Intriguingly, ncRNA-targeted genes are less prone to show positive selection or differential expression during ecotypic divergence due to constraints from ncRNA-mRNA interactions. The methylation levels of ncRNAs generally tend to be higher in temperate lotus than in tropical lotus, and differential methylation of lncRNAs also tends to have expression changes. Overall, our study of ncRNAs and their targets highlights the role of ncRNAs in rhizome growth variation between lotus ecotypes through expression rewiring and methylation modification.

Distinct methylome patterns contribute to ecotypic differentiation in the growth of the storage organ of a flowering plant (sacred lotus)

DNA methylation is an epigenetic modification involved in phenotypic diversity, plant development, and environmental responses. However, the mechanisms of DNA methylation underpinning the adaption of lotus (Nelumbo nucifera) ecotypes to high and low latitudes remain unsolved, especially adaptive evolution of their storage organs. Tropical and temperate lotus ecotypes have thin and enlarged rhizomes which are adapted to low and high latitudes, respectively. Here, we investigated the DNA methylomes and transcriptomes of rhizomes of the temperate and tropical lotus to address this issue. Compared with that of the tropical lotus, the DNA of the temperate lotus was significantly more hypermethylated, indicating an increase in global DNA methylation in the lotus, with rhizome enlargement. Meanwhile, genes associated with differentially methylated regions in their promoters tended to be differentially expressed between the two ecotypes. Interestingly, the genes with their expression negatively correlated with methylation levels in their promoters and genomic regions displayed significantly higher transposon coverage, while the genes showing a significant positive correlation between expression and methylation showed lesser transposon coverage. Further, we identified that DNA methylation, especially in the promoter region, was significantly correlated with the expression of many starch-biosynthetic, gibberellin-, and brassinosteroid-signalling genes associated with rhizome differentiation. Overall, our study unveiled that distinct global and local methylation patterns between the two lotus ecotypes contribute to their expression differences and adaptive phenotypic divergence of their storage organs, highlighting the role of DNA methylation in shaping the ecotypic differentiation of lotus.

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.

Proteome analysis of bermudagrass stolons and rhizomes provides new insights into the adaptation of plant stems to aboveground and underground growth

As an important perennial warm-season turfgrass species, bermudagrass (Cynodon dactylon L.) forms underground-growing rhizomes and aboveground-growing stolons simultaneously, making it a fast propagating clonal plant with strong regeneration ability. In the current study, we compared the internode proteomes of rhizomes and stolons at the same developmental stage in the bermudagrass cultivar Yangjiang using iTRAQ. The results indicated that 228 protein species were differentially accumulated in the two specialized stems. In agreement with the different contents of starch, chlorophyll, anthocyanin and H₂O₂ in the two types of stems, photosynthesis and flavonoid biosynthesis were enriched with differentially accumulated protein species (DAPs) in stolons, whereas starch and sucrose metabolism, glycolysis, and H₂O₂ metabolism were enriched with DAPs in rhizomes. Burying stolons in the soil resulted in the gradual degradation of chlorophyll and anthocyanin, accumulation of starch, and increment of H₂O₂, which is similar to the physiological characteristics of rhizomes. These results collectively revealed that stolons and rhizomes of bermudagrass have significant differences at the proteome level and light might play important regulatory roles in the discrepancy of the proteome profiles and specialization of the two stems, providing new insights into the adaptation of plant stems to aboveground and underground growth. BIOLOGICAL SIGNIFICANCE: As two types of specialized stems that grow underground and aboveground respectively, rhizomes and stolons play important roles in overwintering and ecological invasion of many perennial and clonal plants. However, because rhizomes and stolons rarely coexist in single plant species, the differences between the two stems remain unclear at the molecular level. In this study, through an iTRAQ comparative proteomic analysis, we reported the identification of 228 differentially accumulated protein species (DAPs) in rhizomes and stolons of bermudagrass for the first time. We found that the 228 DAPs were interconnected to form protein networks in regulating diverse cellular activities and biochemical reactions. We also observed that stolons growing underground showed similar physiological activities and DAP expression as those of underground-growing rhizomes, suggesting that light might play important regulatory roles in the specialization of stolons and rhizomes. These results expanded our understanding of the mysterious adaption of plant stems to different growth conditions.

Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology

Background

The AP2/ERF family is widely present in plants and plays a crucial regulatory role in plant growth and development. As an essential aquatic horticultural model plant, lotus has an increasingly prominent economic and research value.

Results

We have identified and analysed the AP2/ERF gene family in the lotus. Initially, 121 AP2/ERF family genes were identified. By analysing their gene distribution and protein structure, and their expression patterns during the development of lotus rhizome, combined with previous studies, we obtained an SNP (megascaffold_20:3578539) associated with lotus rhizome phenotype. This SNP was in the NnADAP gene of the AP2 subfamily, and the changes in SNP (C/T) caused amino acid conversion (proline/leucine). We constructed a population of 95 lotus varieties for SNP verification. Through population typing experiments, we found that the group with SNP CC had significantly larger lotus rhizome and higher soluble sugar content among the population.

Conclusions

In conclusion, we speculate that the alteration of the SNP in the NnADAP can affect the size and sugar content of the lotus rhizome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07473-w.

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.

Developmental regulation of stolon and rhizome

Stolons and rhizomes are modified stems for vegetative reproduction. While stolons grow above the ground, rhizomes grow beneath the ground. Stolons and rhizomes maintain the genotypes of hybrids and hence are invaluable for agricultural propagation. Diploid strawberry is a model for studying stolon development. At the axillary meristems, gibberellins and MADS box gene SOC1 promote stolon formation, while the DELLA repressor inhibits stolon development. Photoperiod regulates stolon formation through regulating GA biosynthesis or balancing asexual with sexual mode of reproduction in the axillary meristems. In rhizomatous wild rice, the BLADE-ON-PETIOLE gene promotes sheath-to-blade ratio to confer rhizome tip stiffness and support underground growth. Together, this review aims to encourage further investigations into stolon and rhizome to benefit agriculture and environment.

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.