Transcriptomic analysis identifies the key genes involved in stamen petaloid in lotus (Nelumbo nucifera)

Cell wall remodeling confers plant architecture with distinct wall structure in Nelumbo nucifera

Lotus (Nelumbo nucifera G.) is a perennial aquatic horticultural plant with diverse architectures. Distinct plant architecture (PA) has certain attractive and practical qualities, but its genetic morphogenesis in lotus remains elusive. In this study, we employ genome-wide association analysis (GWAS) for the seven traits of petiole length (PLL), leaf length (LL), leaf width (LW), peduncle length (PLF), flower diameter (FD), petal length (PeL), and petal width (PeW) in 301 lotus accessions. A total of 90 loci are identified to associate with these traits across 4 years of trials. Meanwhile, we perform RNA sequencing (RNA-seq) to analyze the differential expression of the gene (DEG) transcripts between large and small PA (LPA and SPA) of lotus stems (peduncles and petioles). As a result, eight key candidate genes are identified that are all primarily involved in plant cell wall remodeling significantly associated with PA traits by integrating the results of DEGs and GWAS. To verify this result, we compare the cell wall compositions and structures of LPA versus SPA in representative lotus germplasms. Intriguingly, compared with the SPA lotus, the LPA varieties have higher content of cellulose and hemicellulose, but less filling substrates of pectin and lignin. Additionally, we verified longer cellulose chains and higher cellulose crystallinity with less interference in LPA varieties. Taken together, our study illustrates how plant cell wall remodeling affects PA in lotus, shedding light on the genetic architecture of this significant ornamental trait and offering a priceless genetic resource for future genomic-enabled breeding.

Transcriptomics analyses reveal the key genes involved in stamen petaloid formation in Alcea rosea L

Alcea rosea L. is a traditional flower with a long cultivation history. It is extensively cultivated in China and is widely planted in green belt parks or used as cut flowers and potted ornamental because of its rich colors and flower shapes. Double-petal A. rosea flowers have a higher aesthetic value compared to single-petal flowers, a phenomenon determined by stamen petaloid. However, the underlying molecular mechanism of this phenomenon is still very unclear. In this study, an RNA-based comparative transcriptomic analysis was performed between the normal petal and stamen petaloid petal of A. rosea. A total of 3,212 differential expressed genes (DEGs), including 2,620 up-regulated DEGs and 592 down-regulated DEGs, were identified from 206,188 unigenes. Numerous DEGs associated with stamen petaloid were identified through GO and KEGG enrichment analysis. Notably, there were 63 DEGs involved in the plant hormone synthesis and signal transduction, including auxin, cytokinin, gibberellin, abscisic acid, ethylene, brassinosteroid, jasmonic acid, and salicylic acid signaling pathway and 56 key transcription factors (TFs), such as MADS-box, bHLH, GRAS, and HSF. The identification of these DEGs provides an important clue for studying the regulation pathway and mechanism of stamen petaloid formation in A. rosea and provides valuable information for molecular plant 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.

Combined Transcriptome and Proteome Analysis Provides Insights into Petaloidy in Pomegranate

Petaloidy leads to a plump floral pattern and increases the landscape value of ornamental pomegranates; however, research on the mechanism of petaloidy in ornamental pomegranates is limited. In this study, we aimed to screen candidate genes related to petaloidy. We performed transcriptomic and proteomic sequencing of the stamens and petals of single-petal and double-petal flowers of ornamental pomegranates. Briefly, 24,567 genes and 5865 proteins were identified, of which 5721 genes were quantified at both transcriptional and translational levels. In the petal and stamen comparison groups, the association between differentially abundant proteins (DAPs) and differentially expressed genes (DEGs) was higher than that between all genes and all proteins, indicating that petaloidy impacts the correlation between genes and proteins. The enrichment results of transcriptome, proteome, and correlation analyses showed that cell wall metabolism, jasmonic acid signal transduction, redox balance, and transmembrane transport affected petaloidy. Nine hormone-related DEGs/DAPs were selected, among which ARF, ILR1, LAX2, and JAR1 may promote petal doubling. Sixteen transcription factor DEGs/DAPs were selected, among which EREBP, LOB, MEF2, MYB, C3H, and trihelix may promote petal doubling. Our results provide transcriptomic and proteomic data on the formation mechanism of petaloidy and a theoretical basis for breeding new ornamental pomegranate varieties.

Plastid Phylogenomics of Paeonia and the Evolution of Ten Flower Types in Tree Peony

Paeonia suffruticosa Andr., a member of Paeoniaceae, is native to China. In its 1600 years' cultivation, more than 2000 cultivars for different purposes (ornamental, medicinal and oil use) have been inbred. However, there are still some controversies regarding the provenance of tree peony cultivars and the phylogenetic relationships between and within different cultivar groups. In this study, plastid genome sequencing was performed on 10 representative tree peony cultivars corresponding to 10 different flower types. Structure and comparative analyses of the plastid genomes showed that the total lengths of the chloroplast genome of the 10 cultivars ranged from 152,153 to 152,385 bp and encoded 84-88 protein-coding genes, 8 rRNAs and 31-40 tRNAs. The number of simple sequence repeats and interspersed repeat sequences of the 10 cultivars ranged from 65-68 and 40-42, respectively. Plastid phylogenetic relationships of Paeonia species/cultivars were reconstructed incorporating data from our newly sequenced plastid genomes and 15 published species, and results showed that subsect. Vaginatae was the closest relative to the central plains cultivar group with robust support, and that it may be involved in the formation of the group. Paeonia ostii was recovered as a successive sister group to this lineage. Additionally, eleven morphological characteristics of flowers were mapped to the phylogenetic skeleton to reconstruct the evolutionary trajectory of flower architecture in Paeoniaceae.

Genome-wide association study of traits in sacred lotus uncovers MITE-associated variants underlying stamen petaloid and petal number variations

Understanding the genetic variants responsible for floral trait diversity is important for the molecular breeding of ornamental flowers. Widely used in water gardening for thousands of years, the sacred lotus exhibits a wide range of diversity in floral organs. Nevertheless, the genetic variations underlying various morphological characteristics in lotus remain largely unclear. Here, we performed a genome-wide association study of sacred lotus for 12 well-recorded ornamental traits. Given a moderate linkage disequilibrium level of 32.9 kb, we successfully identified 149 candidate genes responsible for seven flower traits and plant size variations, including many pleiotropic genes affecting multiple floral-organ-related traits, such as NnKUP2. Notably, we found a 2.75-kb presence-and-absence genomic fragment significantly associated with stamen petaloid and petal number variations, which was further confirmed by re-examining another independent population dataset with petal number records. Intriguingly, this fragment carries MITE transposons bound by siRNAs and is related to the expression differentiation of a nearby candidate gene between few-petalled and double-petalled lotuses. Overall, these genetic variations and candidate genes responsible for diverse lotus traits revealed by our GWAS highlight the role of transposon variations, particularly MITEs, in shaping floral trait diversity.

Widely targeted metabolomics reveals stamen petaloid tissue of Paeonia lactiflora Pall. being a potential pharmacological resource

Paeonia lactiflflora Pall. has a long edible and medicinal history because of the very high content of biologically active compounds. However, little information is available about the metabolic basis of pharmacological values of P. lactiflora flowers. In this study, we investigated metabolites in the different parts of P. lactiflora flowers, including petal, stamen petaloid tissue and stamen, by widely targeted metabolomics approach. A total of 1102 metabolites were identified, among which 313 and 410 metabolites showed differential accumulation in comparison groups of petal vs. stamen petaloid tissue and stamen vs. stamen petaloid tissue. Differential accumulated metabolites analysis and KEGG pathway analysis showed that the flavonoids were the most critical differential metabolites. Furthermore, difference accumulation of flavonoids, phenolic acids, tannins and alkaloids might lead to the differences in antioxidant activities and tyrosinase inhibition effects. Indeed, stamen petaloid tissue displayed better antioxidant and anti-melanin production activities than petal and stamen through experimental verification. These results not only expand our understanding of metabolites in P. lactiflora flowers, but also reveal that the stamen petaloid tissues of P. lactiflora hold the great potential as promising ingredients for pharmaceuticals, functional foods and skincare products.

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.

Studies on Lotus Genomics and the Contribution to Its Breeding

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

Comparative transcriptomic analysis provides insight into carpel petaloidy in lotus (Nelumbo nucifera)

Lotus (Nelumbo nucifera) is a highly recognized flower with high ornamental value. Flower color and flower morphology are two main factors for flower lotus breeding. Petaloidy is a universal phenomenon in lotus flowers. However, the genetic regulation of floral organ petaloidy in lotus remains elusive. In this study, the transcriptomic analysis was performed among three organs, including petal, carpel petaloidy, and carpel in lotus. A total of 1,568 DEGs related to carpel petaloidy were identified. Our study identified one floral homeotic gene encoded by the MADS-box transcription factor, AGAMOUS (AG) as the candidate gene for petaloid in lotus. Meanwhile, a predicted labile boundary in floral organs of N. nucifera was hypothesized. In summary, our results explored the candidate genes related to carpel petaloidy, setting a theoretical basis for the molecular regulation of petaloid phenotype.

Development of the petaloid bracts of a paleoherb species, Saururus chinensis

Saururus chinensis is a core member of Saururaceae, an ancient, perianthless (lacking petals or sepals) family of the magnoliids in the Mesangiospermae, which is important for understanding the origin and evolution of early flowers due to its unusual floral composition and petaloid bracts. To compare their transcriptomes, RNA-seq abundance analysis identified 43,463 genes that were found to be differentially expressed in S. chinensis bracts. Of these, 5,797 showed significant differential expression, of which 1,770 were up-regulated and 4,027 down-regulated in green compared to white bracts. The expression profiles were also compared using cDNA microarrays, which identified 166 additional differentially expressed genes. Subsequently, qRT-PCR was used to verify and extend the cDNA microarray results, showing that the A and B class MADS-box genes were up-regulated in the white bracts. Phylogenetic analysis was performed on putative S. chinensis A and B-class of MADS-box genes to infer evolutionary relationships within the A and B-class of MADS-box gene family. In addition, nature selection and protein interactions of B class MADS-box proteins were inferred that B-class genes free from evolutionary pressures. The results indicate that petaloid bracts display anatomical and gene expression features normally associated with petals, as found in petaloid bracts of other species, and support an evolutionarily conserved developmental program for petaloid bracts.

Molecular Cloning and Expression Profile of Class E Genes Related to Sepal Development in Nelumbo nucifera

The lotus (Nelumbo Adans.) is an important aquatic plant with ornamental, medicinal and edible values and cultural connotations. It has single-, semi-double-, double- and thousand-petalled types of flower shape and is an ideal material for developmental research of flower doubling. The lotus is a basal eudicot species without a morphological difference between the sepals and petals and occupies a critical phylogenetic position in flowering plants. In order to investigate the genetic relationship between the sepals and petals in the lotus, the class E genes which affect sepal formation were focused on and analyzed. Here, SEPALLATA 1(NnSEP1) and its homologous genes AGAMOUS-LIKE MADS-BOXAGL9 (NnAGL9) and MADS-BOX TRANSCRIPTION FACTOR 6-like (NnMADS6-like) of the class E gene family were isolated from the flower buds of the Asian lotus (Nelumbo nucifera Gaertn.). The protein structure, subcellular localization and expression patterns of these three genes were investigated. All three genes were verified to locate in the nucleus and had typical MADS-box characteristics. NnSEP1 and NnMADS6-like were specifically expressed in the sepals, while NnAGL9 was highly expressed in the petals, suggesting that different developmental mechanisms exist in the formation of the sepals and petals in the lotus. The significant functional differences between NnSEP1, NnMADS6-like and NnAGL9 were also confirmed by a yeast two-hybrid assay. These results expand our knowledge on the class E gene family in sepal formation and will benefit fundamental research on the development of floral organs in Nelumbo.

Isolation and Functional Analysis of Genes Involved in Polyacylated Anthocyanin Biosynthesis in Blue Senecio cruentus

Polyacylated anthocyanins with multiple glycosyl and aromatic acyl groups tend to make flowers display bright and stable blue colours. However, there are few studies on the isolation and functional characterization of genes involved in the polyacylated anthocyanin biosynthesis mechanism, which limits the molecular breeding of truly blue flowers. Senecio cruentus is an important potted ornamental plant, and its blue flowers contain 3',7-polyacylated delphinidin-type anthocyanins that are not reported in any other plants, suggesting that it harbours abundant gene resources for the molecular breeding of blue flowers. In this study, using high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) analysis of blue, carmine and white colours of cineraria cultivars "Venezia" (named VeB, VeC, and VeW, respectively), we found that 3',7-polyacylated anthocyanin, cinerarin, was the main pigment component that determined the blue colour of ray florets of cineraria. Based on the transcriptome sequencing and differential gene expression (DEG) analysis combined with RT- and qRT-PCR, we found two genes encoding uridine diphosphate glycosyltransferase, named ScUGT1 and ScUGT4; two genes encoding acyl-glucoside-dependent glucosyltransferases which belong to glycoside hydrolase family 1 (GH1), named ScAGGT11 and ScAGGT12; one gene encoding serine carboxypeptidase-like acyltransferase ScSCPL2; and two MYB transcriptional factor genes ScMYB2 and ScMYB4, that were specifically highly expressed in the ray florets of VeB, which indicated that these genes may be involved in cinerarin biosynthesis. The function of ScSCPL2 was analysed by virus-induced gene silencing (VIGS) in cineraria leaves combined with HPLC-MS/MS. ScSCPL2 mainly participated in the 3' and 7-position acylation of cinerarin. These results will provide new insight into the molecular basis of the polyacylated anthocyanin biosynthesis mechanism in higher plants and are of great significance for blue flower molecular breeding of ornamental plants.

The complexity of alternative splicing and landscape of tissue-specific expression in lotus (Nelumbo nucifera) unveiled by Illumina- and single-molecule real-time-based RNA-sequencing

Alternative splicing (AS) plays a critical role in regulating different physiological and developmental processes in eukaryotes, by dramatically increasing the diversity of the transcriptome and the proteome. However, the saturation and complexity of AS remain unclear in lotus due to its limitation of rare obtainment of full-length multiple-splice isoforms. In this study, we apply a hybrid assembly strategy by combining single-molecule real-time sequencing and Illumina RNA-seq to get a comprehensive insight into the lotus transcriptomic landscape. We identified 211,802 high-quality full-length non-chimeric reads, with 192,690 non-redundant isoforms, and updated the lotus reference gene model. Moreover, our analysis identified a total of 104,288 AS events from 16,543 genes, with alternative 3ʹ splice-site being the predominant model, following by intron retention. By exploring tissue datasets, 370 tissue-specific AS events were identified among 12 tissues. Both the tissue-specific genes and isoforms might play important roles in tissue or organ development, and are suitable for ‘ABCE’ model partly in floral tissues. A large number of AS events and isoform variants identified in our study enhance the understanding of transcriptional diversity in lotus, and provide valuable resource for further functional genomic studies.

Differentially Expressed Genes between Carrot Petaloid Cytoplasmic Male Sterile and Maintainer during Floral Development

Petaloid cytoplasmic male sterility (CMS) is a maternally inherited loss of male fertility due to the complete conversion of stamens into petal-like organs, and CMS lines have been widely utilized in carrot breeding. Petaloid CMS is an ideal model not only for studying the mitochondrial-nuclear interaction but also for discovering genes that are essential for floral organ development. To investigate the comprehensive mechanism of CMS and homeotic organ alternation during carrot flower development, we conducted transcriptome analysis between the petaloid CMS line (P2S) and its maintainer line (P2M) at four flower developmental stages (T1-T4). A total of 2838 genes were found to be differentially expressed, among which 1495 genes were significantly downregulated and 1343 genes were significantly upregulated in the CMS line. Functional analysis showed that most of the differentially expressed genes (DEGs) were involved in protein processing in the endoplasmic reticulum, plant hormone signal transduction, and biosynthesis. A total of 16 MADS-box genes were grouped into class A, B, C, and E, but not class D, genes. Several key genes associated with oxidative phosphorylation showed continuously low expression from stage T2 in P2S, and the expression of DcPI and DcAG-like genes also greatly decreased at stage T2 in P2S. This indicated that energy deficiency might inhibit the expression of B- and C-class MADS-box genes resulting in the conversion of stamens into petals. Stamen petaloidy may act as an intrinsic stress, upregulating the expression of heat shock protein (HSP) genes and MADS-box genes at stages T3 and T4 in P2S, which results in some fertile revertants. This study will provide a better understanding of carrot petaloid CMS and floral development as a basis for further research.

Transcriptome analysis during floral organ development provides insights into stamen petaloidy in Lagerstroemia speciosa

As one of the most popular woody species that blooms in summer, Lagerstroemia speciosa has been used abundantly in urban landscape for its excellent floral beauty. For the first time, we discovered a double-flower variant with all petaloid stamens. To understand the molecular basis of this variation, we contrasted the transcriptomes of single- and double-flower buds at three stamen development stages. In total, 73,536 unigenes were mapped and 30,714 differently expressed genes (DEGs) were identified in the tissues. We focused on the DEGs expressing in both phenotypes and investigated the association of their expression profiles with their functions in transcription pathways. Furthermore, we performed WGCNA and identified co-expressed genes with four floral homeotic genes as hubs (MADS16, Unigene0026169; AP2, Unigene0042732; SOC1, Unigene0046314; AG, Unigene0056437). The expression of these hub genes has been conserved across the three developmental stages but significantly different between the two floral phenotypes. As a result, the robust transcriptional regulation of stamen petaloidy in double flowers was deduced. These findings will help to unravel the regulatory mechanisms of several specific genes, thereby providing a basis to study double-flower molecular breeding in L. speciosa.

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

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

Genome-wide analysis of spatiotemporal gene expression patterns during floral organ development in Brassica rapa

Flowering is a key agronomic trait that directly influences crop yield and quality and serves as a model system for elucidating the molecular basis that controls successful reproduction, adaptation, and diversification of flowering plants. Adequate knowledge of continuous series of expression data from the floral transition to maturation is lacking in Brassica rapa. To unravel the genome expression associated with the development of early small floral buds (< 2 mm; FB2), early large floral buds (2-4 mm; FB4), stamens (STs) and carpels (CPs), transcriptome profiling was carried out with a Br300K oligo microarray. The results showed that at least 6848 known nonredundant genes (30% of the genes of the Br300K) were differentially expressed during the floral transition from vegetative tissues to maturation. Functional annotation of the differentially expressed genes (DEGs) (fold change ≥ 5) by comparison with a close relative, Arabidopsis thaliana, revealed 6552 unigenes (4579 upregulated; 1973 downregulated), including 131 Brassica-specific and 116 functionally known floral Arabidopsis homologs. Additionally, 1723, 236 and 232 DEGs were preferentially expressed in the tissues of STs, FB2, and CPs. These DEGs also included 43 transcription factors, mainly AP2/ERF-ERF, NAC, MADS-MIKC, C2H2, bHLH, and WRKY members. The differential gene expression during flower development induced dramatic changes in activities related to metabolic processes (23.7%), cellular (22.7%) processes, responses to the stimuli (7.5%) and reproduction (1%). A relatively large number of DEGs were observed in STs and were overrepresented by photosynthesis-related activities. Subsequent analysis via semiquantitative RT-PCR, histological analysis performed with in situ hybridization of BrLTP1 and transgenic reporter lines (BrLTP promoter::GUS) of B. rapa ssp. pekinensis supported the spatiotemporal expression patterns. Together, these results suggest that a temporally and spatially regulated process of the selective expression of distinct fractions of the same genome leads to the development of floral organs. Interestingly, most of the differentially expressed floral transcripts were located on chromosomes 3 and 9. This study generated a genome expression atlas of the early floral transition to maturation that represented the flowering regulatory elements of Brassica rapa.

Transcription Profiles Reveal the Regulatory Synthesis of Phenols during the Development of Lotus Rhizome (Nelumbo nucifera Gaertn)

Lotus (Nelumbo nucifera Gaertn) is a wetland vegetable famous for its nutritional and medicinal value. Phenolic compounds are secondary metabolites that play important roles in the browning of fresh-cut fruits and vegetables, and chemical constituents are extracted from lotus for medicine due to their high antioxidant activity. Studies have explored in depth the changes in phenolic compounds during browning, while little is known about their synthesis during the formation of lotus rhizome. In this study, transcriptomic analyses of six samples were performed during lotus rhizome formation using a high-throughput tag sequencing technique. About 23 million high-quality reads were generated, and 92.14% of the data was mapped to the reference genome. The samples were divided into two stages, and we identified 23,475 genes in total, 689 of which were involved in the biosynthesis of secondary metabolites. A complex genetic crosstalk-regulated network involved in the biosynthesis of phenolic compounds was found during the development of lotus rhizome, and 25 genes in the phenylpropanoid biosynthesis pathway, 18 genes in the pentose phosphate pathway, and 30 genes in the flavonoid biosynthesis pathway were highly expressed. The expression patterns of key enzymes assigned to the synthesis of phenolic compounds were analyzed. Moreover, several differentially expressed genes required for phenolic compound biosynthesis detected by comparative transcriptomic analysis were verified through qRT-PCR. This work lays a foundation for future studies on the molecular mechanisms of phenolic compound biosynthesis during rhizome formation.

Genome-Wide DNA Methylation Profiling in the Lotus (Nelumbo nucifera) Flower Showing its Contribution to the Stamen Petaloid

DNA methylation is a vital epigenetic modification. Methylation has a significant effect on the gene expression influencing the regulation of different physiological processes. Current studies on DNA methylation have been conducted on model plants. Lotus (Nelumbo nucifera) is a basic eudicot exhibiting variations during development, especially in flower formation. DNA methylation profiling was conducted on different flower tissues of lotuses through whole genome bisulfite sequencing (WGBS) to investigate the effects of DNA methylation on its stamen petaloid. A map of methylated cytosines at the single base pair resolution for the lotus was constructed. When the stamen was compared with the stamen petaloid, the DNA methylation exhibited a global decrease. Genome-wide relationship analysis between DNA methylation and gene expression identified 31 different methylation region (DMR)-associated genes, which might play crucial roles in floral organ formation, especially in the stamen petaloid. One out of 31 DMR-associated genes, NNU_05638 was homolog with Plant U-box 33 (PUB33). The DNA methylation status of NNU_05638 promoter was distinct in three floral organs, which was confirmed by traditional bisulfite sequencing. These results provide further insights about the regulation of stamen petaloids at the epigenetic level in lotus.