- Research article
- Open Access
Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology
BMC Genomics volume 22, Article number: 171 (2021)
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.
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.
In conclusion, we speculate that the alteration of the SNP in the NnADAP can affect the size and sugar content of the lotus rhizome.
Chinese lotus (Nelumbo nucifera) is popularly used for food and ornaments with agricultural and economic value. Based on their biological characteristics and observed ecological habits, field visits, and regional experiments, the lotus was documented to exist in two natural ecotypes, namely, temperate lotus and tropical lotus. These two ecotypes have distinctly different growth and development habits, especially in their underground stem development. Temperate lotus has an apparent annual growth and development cycle, withering in autumn while the underground stem expands into lotus rhizome and hibernates. Tropical lotus in a tropical climate will not have an enlarged rhizome but grows into a whip shaped rhizome during the entire growth period [1, 2]. In temperate lotus, the development stages of underground stem could be separated into four (S1 to S4) depending on the morphology [3, 4], S1 is described as the stolon stage, S2 as the middle stage where extension of underground stems occurs, in S3 thickening of underground stems occurs, and in S4 there is a continuous thickening of underground stems as well as simultaneous accumulation of starch and extension of internodes.
The AP2/ERF (APETALA2/Ethylene Responsive Factor) is a sizable transcription factors family in the plant kingdom. Their classification was widely studied in many model plants, including Arabidopsis and rice (Oryza sativa) , and many non-model plants such as the rubber tree (Hevea brasiliensis) , tea plant (Camellia sinensis) , and the Chinese cabbage (Brassica rapa ssp. pekinensis) . According to the domain number and gene sequence of the AP2 genes, the massive group has been divided into three subgroups , including the AP2 subfamily (containing two AP2/ERFs domain), the ERF subfamily (containing one AP2/ERF domain), and the RAV subfamily with one AP2/ERF domain. Besides, Sakuma classified the AP2/ERF family into five branches, including AP2, DREB, RAV, ERF, and Soloist .
The AP2/ERF superfamily members were demonstrated to be involved in plant growth [10,11,12,13], response to stress [14,15,16,17], and metabolism [18, 19]. The AP2 subfamily genes, WRINKLED (WRI), are involved in fatty acid production and biosynthesis pathway and have different biological functions in Arabidopsis [20, 21]. Besides, AP2/ERF proteins are involved in the carbohydrate biosynthesis and metabolism-regulation network in storage organs of Arabidopsis , rice (Oryza sativa) , and maize (Zea mays) [17, 24]. Nevertheless, previous studies also showed that the AP2 subfamily (L1140) play critical roles in the formation of lotus rhizome . Furthermore, two AP2-like transcription factors (NNU_12870 and NNU_17043) showed the same or an opposite trend with rhizome formation in tropical and temperate lotus . Here we have given an overview of this family and its expression in lotus rhizome, elucidating the relation between the AP2/ERF family and lotus rhizome development.
Quantitative trait loci (QTL) mapping studies provide numerous DNA markers and thus contributing to the improvement of plant breeding via marker-assisted selection (MAS) [25, 26]. The genetic linkage map of various plant types and floral organ characteristics in lotus has accelerated the process of molecular breeding and crop improvement [27, 28]. Although the genetic map with high-density of lotus has been constructed and QTL mapping of traits related to underground stem development has been obtained , the verification of specific genes needs to be studied further. Single nucleotide polymorphism (SNP) markers are vital for MAS, as even one SNP in 5′ regulatory region was capable of causing disability of seed shattering in rice  and SNP resulting in a truncated SHATTERING 3 (SH3)/SEED SHATTERING 4 (SH4) protein caused the seed shattering disability in Oryza glaberrima .
Previously, we detected several SNP markers in some genes in lotus, which showed a differential expression pattern between temperate and tropical lotus [4, 32]. To further confirm their potential relationship with the underground stem enlargement phenotypes, we selected one SNP in NnADAP (NNU_25830) belonging to the AP2 subfamily to show its preference in lotus varieties with various rhizome morphological characteristics. SNP genotyping showed that this SNP is associated with the agronomic traits of lotus rhizome. Besides, we analysed the expression of NnADAP and established a way that linked this marker with lotus rhizome traits.
Identification, distribution, and phylogenetic analysis of the AP2/ERF superfamily genes in lotus
The AP2 domain (PF00847) was blasted in the Lotus Database , where 121 predicted AP2/ERF proteins were identified. The length of these genes ranged from 278 to 8.01 k bp, and their corresponding polypeptide sequences ranged from 82 to 1.01 k aa (Additional file 1: Table S1). According to sequencing data, except for the soloist gene, the other 120 AP2/ERF family genes were mostly (90.8%) mapped to the ten megascaffolds 1–10 (Additional file 2: Fig. S1A). The rest (9.2%) were distributed in other small megascaffolds (Additional file 2: Fig. S1B). Megascaffold 1 and megascaffold 2 were the largest two megascaffolds (255 and 133 Mb), containing 55% of the total members (Additional file 2: Fig. S1A).
We constructed a phylogenetic analysis of AP2/ERF superfamily with Nelumbo nucifera and Arabidopsis by employing the maximum likelihood (ML) method. Sequences of 121 NnAP2/ERF genes and 143 AtAP2/ERF genes were retrieved from the Lotus Database  and the Arabidopsis information resource website, respectively . The phylogenetic tree in lotus showed that NnAP2/ERF genes were separated into five subfamilies according to the Arabidopsis classification method, including AP2, ERF, DREB, RAV, and Soloist (Fig. 1) [5, 37]. Eighteen genes belonging to the AP2 subfamily, seventeen of which had two AP2 domains and the one gene having only one AP2 domain. There were 55 genes in ERF subfamily and 42 genes belonging to the DREB subfamily, which had only one AP2 domain, and five genes in RAV subfamily with both an AP2 and B3 domain. There was one soloist gene (newGene_1464) sharing a similar protein sequence with Arabidopsis soloist gene At4g13040, which possessed one AP2 domain (Additional file 1: Table S1).
Structure analysis of lotus AP2/ERF superfamily
After classifying the AP2/ERF superfamily in lotus, we analysed its gene and protein structure. The gene structure analysis was performed through the GSDS (Gene Structure Display Server) 2.0 . The number of exon and intron differed in the four subgroups of the AP2/ERF family. The results showed that genes in AP2 subfamily possess more than ten exons individually. Genes in RAV subcategory and those in the soloist subfamily had two to five exons each, while 29.1% of members of the ERF subgroup had two exons with the remaining genes having only one exon. In the DREB subfamily, most of the genes had one exon with only four genes exceeding one exon (Fig. 2). To further study the characteristics of NnAP2 proteins, the motifs of 121 NnAP2 proteins were analysed using MEME (Multiple EM for Motif Elicitation, V5.0.1 ). Ten predicted motifs were set as default parameters and were listed in Additional file 1: Table S2. Members of the AP2 subgroup mainly own four to six motifs. All the AP2 subgroup proteins possessed motif 2 and motif 3 with the majority having motif 1, motif 6, and motif 3. The combination of motif 1, motif 2, motif 3, and motif 4 was similar to those in ERF subfamily. The comparison of DREB AND ERF had 37.8% of DREB subgroup proteins being similar with those of ERF subfamily, while 20% of the proteins had an additional motif (motif 8), and 48.9% had motif 7. All RAV subgroup proteins included motif 3, motif 1, motif 4, and motif 9, while the soloist protein had only motif 3. The structure analyses assembled using TBtools by Chen  were shown in Fig. 2. From these results, we deduced that most of the clustered genes in the phylogenetic tree share similar patterns of motif combination and may be involved in similar biological functions.
Gene expression analysis of the NnAP2/ERFs during lotus rhizome formation
To assess the prospective function of NnAP2/ERFs genes in the process of rhizome formation and enlargement, we explored gene expression based on previous RNA-seq data . The results showed that 63.6% of AP2/ERF superfamily members were detected during lotus rhizome formation, including eight AP2 subfamily members, 41 ERF genes, 23 DREB genes, four RAV members, and one soloist gene (Additional file 1: Table S3, Fig. 3). The overall expression of eight AP2 subgroup members was low with NNU_21756 and NNU_17043 showing a decreasing expression along with lotus rhizome enlargement. In contrast, the expression of NNU_12443 reached its peak during the middle stage of rhizome development. The other five genes showed different expressions between temperate and tropical lotus. Four RAV genes and one soloist gene showed no difference in their expression trend during the development process of tropical and temperate lotus underground stem. In DREB subgroup, 82.6 and 91.3% genes were down regulated in the rhizome of temperate lotus and tropical lotus, with a similar trend being observed in the ERF subfamily. More than half of the lotus AP2/ERF superfamily were expressed in lotus rhizome. Most of them showed a different expression pattern between temperate and tropical lotus, which indicates that this family may play a role in rhizome enlargement. We, therefore, focused on the genes that were up-regulated in temperate lotus and down-regulated in tropical lotus.
Acquisition and verification of target SNP in NnADAP
After the second-generation genome sequencing of several lotus varieties with different rhizome perimeters , we cloned an AP2 subfamily gene, NNU_25830 (herein known as NnADAP), which possesses an interesting SNP (megascaffold_20, 3,578,539) which differ with the change of rhizome perimeter. Then we verified the second-generation genome sequencing results by Sanger sequencing using genomic primers (NnADAP-G forward and reverse primers, Additional file 1: Table S4), and the sanger sequence results mostly kept congruent with high-throughput results (Additional file 1: Table S5). In tropical lotus, the SNP was “TT”, while in temperate lotus, it was “CC” or “C/T”, and the single basement change leads to variation in the amino acid sequence of proline to leucine.
To confirm the preference and effect of the SNP in different phenotypes of lotus rhizome, we collected a population of 95 lotus varieties for SNP genotype analysis (Additional file 1: Table S6). Among the population, the genome of 18 varieties have been sequenced , and this whole population included wild temperate lotus, cultivated temperate lotus, tropical lotus, N. lutea, and several hybrids, depending on the agricultural application and experience in cultivation [32, 41]. We conducted Kompetitive Allele-Specific PCR (KASP)  on the flanking sequences containing the target SNP with a pair of KASP primers (KASP forward and reverse primers, Additional file 1: Table S4). After KASP genotyping in the lotus collection, our target SNP had a clear preference in various lotus varieties (Fig. 4a; Additional file 1: Table S6). Thus, we divided the lotus collection into three groups named group I (Allele 1: CC), group II (Allele 2: TT), and group III (Heterozygote: CT) depending on the SNP genotype.
Association between the target SNP in NnADAP and rhizome biological indices in lotus
Lotus rhizome enlargement was accompanied by changes in morphological characteristics and sugar and starch content . In this study, we focused on the character of the three lotus groups to elucidate the function of the target SNP on NnADAP in traits of rhizome, including REI, perimeter, and soluble sugar content (Fig. 4b-d, Additional file 1: Table S6). The REI and perimeter in group I were significantly higher than the other two groups. At the same time, there was no difference between group II and group III, implying that lotus rhizome size in the SNP genotype CC plants may be larger. The soluble sugar content in group I showed a significant increase than group II and a significant decrease compared with group III. Group II had the lowest soluble sugar content with the SNP genotype TT, and the other two groups with CC and CT showed a higher sugar content, which probably revealed that SNP CC promotes soluble sugar accumulation.
The temporal-spatial expression of NnADAP in temperate lotus
Previous transcriptome data illuminates the differences in NnADAP expression in the lotus rhizome development process in temperate and tropical lotus . To further analyse the role played by NnADAP (with target SNP CC) in rhizome formation, we investigated the gene expression of critical plant parts during the temperate lotus formation process (Fig. 5). The expression of NnADAP was high in the petiole and rhizome, with the highest expression occurring at S2 stage of rhizome development. In the leaves, the expression of NnADAP was exceedingly significantly increased in S4. Together with the increasement of the starch and soluble sugar in temperate lotus rhizome [4, 43], we deduced that NnADAP (with target SNP CC) might be involved the accumulation of carbohydrate in leaves and carbohydrate partitioning.
The AP2/ERF family is widespread in plants and has many members that are involved in many phases of plant growth and development, including phytohormones biosynthesis [44, 45], floral organ development [46,47,48], and stress resistance [49,50,51,52]. As the lotus is becoming one of the model plants in the horticultural industry and offering immense contribution and allowing more in-depth studies focused on this species [53, 54], the analysis of lotus AP2/ERF superfamily lotus has a great value. Lotus rhizome has a crispy flavor when fresh and an agreeable taste after simmered in soup. Its richness in water, starch, and vitamin, lotus rhizome was used as a popular vegetable and snack in Asia. Several types of research have focused on food usage of lotus rhizome, including the changes that occur during cooking and some other processing method, including heat blanching [55, 56], the phenols and starches, and ADP-glucose pyrophosphorylase [57,58,59]. Further studies on lotus gene families and gene markers of economic traits are of great significance to molecular assisted lotus breeding.
Characteristics of AP2/ERF gene family in lotus
In the present study, we identified 121 NnAP2/ERF superfamily genes in lotus, which consisted of 42 DREB subgroup genes, 55 ERF genes, 18 AP2 members, five RAV subfamily genes, and one soloist gene (Fig. 1). The subfamilies ratio is similar to that in Arabidopsis and longan (Dimocarpus longan Lour.) , with the total numbers being less than that in Arabidopsis and several other species, including Chinese cabbage(Brassica rapa ssp. pekinensis) , grape (Vitis vinifera) , rice , and poplar (Populus trichocarpa) .
According to the gene structure and evolution rates, the AP2/ERF superfamily was separated into five subfamilies . The AP2 subfamily has the most introns and exons, and its protein structure is dominated by the combination of Motif 3, 2, 6, 5, 1, and 4 (Fig. 2). The difference between these members is the number of amino acids between different Motifs. In the DREB and ERF subfamilies, many genes lack introns. The ERF subfamily has 67.3% of genes without intron structure, while in the DREB the number is as high as 83.3%. It also has a similar phenomenon in tartary buckwheat (Fagopyum Tataricum)  and longan . The protein structures of the ERF subgroup and DREB subgroup has certain similarities. The DREB subfamily protein structure has more Motif 7 and Motif 8 than ERF. The protein structure of the RAV subfamily is the most conserved, with fixed combinations of Motif 3, 1, 4, and 9, while the protein structure of the soloist gene NewGene1464 only lacked Motif 3.
The expression of AP2/ERF family in lotus rhizome
In this study, eight AP2 subfamily genes expressed in the lotus underground stem with a low-level expression during lotus rhizome formation process, and most of them showed a different expression trend between temperate and tropical lotus. Flowering and tuberization are the representatives of reproductive and vegetative growth, respectively, which share an overlapping plant energy allocation. AP2 and AP2-like genes are involved in the initiation of flower growth [64,65,66,67], and some of them cause inhibition of flowering in various flowering plants, including rice, Arabidopsis, and Arabis alpina [48, 68, 69]. Flowering and tuber formation were triggered by certain overlapping signals in potato (Solanum tuberosum) , and the transition between flowering and tuberization was induced by ‘florigen’ and ‘tuberigen’, which share similar signaling pathway [71, 72]. In this study, the gene expression of NnADAP reached the highest in lotus rhizome and petiole in the initial of rhizome enlargement at S2 and reached the highest expression in leaf at S4. For the flowering time in temperate lotus is short (several days from S2 to S3), we speculated that NnADAP in lotus might also play an antagonistic role between lotus flowering and rhizome development.
DREB and ERF genes are widely researched in various plants focusing on stress resistance. Debbarma et al.  reviewed the ERF subfamily members in abiotic stresses tolerance, including drought, salinity, heat, cold, oxidation, and others in many plants, such as Arabidopsis , rice , tomato (Lycopersicon esculentum) , and cotton (Gossypium hirsutum) . There were 23 genes out of 42 DREB subfamily and 41 genes out of 55 ERF subfamily was expressed in the lotus rhizome (Fig. 3), which might indicate the different roles of AP2/ERF genes in lotus rhizome formation. The expression of DREB and ERF subfamily genes in lotus rhizome was high in the first two stages of lotus rhizome development, which may imply that the gene’s function focuses on responding to environmental stress in the early stage of lotus underground stem development. RAV subfamily members are also involved in stress resistance. Arabidopsis RAV1 gene adversely regulates drought and salt tolerance and seed development, and gene expression fluctuates because of plant hormones and abiotic stresses [78, 79]. Tomato SlRAV along with tomato SlERFs enhanced the tolerance to bacterial wilt . In soybean (Glycine max L.), RAV genes seem to increase the resistance of salt and drought and might be involved in varying signaling pathways . A total of four RAV subgroup genes and one soloist gene were detected in lotus rhizome and all of them exhibiting the same expression trend with the development of lotus rhizome in tropical lotus and temperate lotus. Among the four RAV genes, two of them showed an increase along with the rhizome formation while the other two genes showed a decreasing expression trend, indicating that lotus RAVs might play opposite roles in rhizome development.
SNP marker for measuring lotus rhizome production
The SNP marker has been researched in many economic traits of various species, such as maize , soybean , and longan [60, 84]. With the genetic and molecular studies in lotus, many molecular genetic markers have been researched and developed, including insertions and deletions (InDels), simple sequence repeats (SSRs), and abundant SNPs [32, 85, 86]. In lotus, the genome-wide analysis of SNP on AP2 domain demonstrated that, compared with other transcription factor gene families like MYB and RING finger, AP2/ERF gene family in lotus have a lower gene mutation of a loss of function , which implied that AP2 domain suffered a lower selection pressure. The SNPs in the protein-coding region of the AP2 gene are worth further study.
In this study, the target SNP in NnADAP showed a significant association between the SNP genotype and some rhizome phenotypes, including the perimeter and length of one rhizome internode, and sugar content, which constitutes the morphological characteristics and main nutritional characteristics of the rhizome. We collected 95 lotus germplasm resources as the SNP genotyping population, including wild temperate lotus, cultivated temperate lotus, and tropical lotus (Additional file 1: Table S6). Most of the flower characters have been collected, and the genome background of some lotus varieties was released [32, 87]. This study focused on the lotus rhizome traits of these lotus germplasm resources. Through the analysis of underground stem morphology data, the REI of this population ranged from 0.08 to 0.61, which achieves our research goal of the biological phenomenon of lotus rhizome enlarging with SNP alter. Perimeter and REI are important morphological indicators of lotus rhizome [88,89,90,91]. When the SNP is CC, the group has the maximum perimeter and REI, which shows that CC has a higher contribution to promoting lotus rhizome expansion than CT and TT. Also, the SNP genotyping result of the soluble sugar content of the CC-variety group was consistent with the perimeter and REI trends. This phenomenon reflects that there is a specific correlation between the SNP and the sugar content and the rhizome expansion . After preliminary genotyping and phenotyping, the gene markers are efficient for measuring lotus rhizome production and useful for rhizome lotus breeding.
In summary, we identified 121 AP2/ERF superfamily genes in lotus, classified them into five subgroups, analysed gene distribution and structure, and enumerated their expression patterns during the development of lotus rhizome. Additionally, we focused on an interesting SNP (C/T) on the NnADAP gene in the AP2 subfamily. Through verification and SNP genotyping experiments in the population, it was found that the different genotypes of this SNP correspond to the lotus rhizome phenotypes of the population groups. The lotus rhizome morphology data of group I (SNP: CC) was significantly larger than the other groups, and it had higher soluble sugar content. Together with the expression pattern of NnADAP, we concluded that the SNP alteration on this gene influences the phenotype of lotus rhizome.
A total of 95 lotus used in this study are planted and preserved in Wuhan Botanical Garden, Chinese academy of sciences, Wuhan, Hubei, China (N30°32′48.33″, E114°25′3.98″). Of the 95 accessions, including 21 Thai lotus, 67 temperate lotus, and seven wild lotus, and all the seeds and plants are from Wuhan Botanical Garden, Chinese academy of sciences (Additional file 1: Table S6). Wuhan Botanical Garden reserves the right of final interpretation. All the lotus germplasms are planted in the Wuhan Botanical Garden and some of them were analysed and sequenced in previous studies [4, 32]. Leaf samples were collected in June in 2017 and then stored at − 80 °C for DNA and RNA extraction. Rhizome samples were collected in April 2018.
Measurements of rhizome enlargement index, water, and soluble sugar of lotus rhizomes
After we measured the rhizome enlargement index, rhizome samples were dehydrated by using an oven heated at 100 °C for 1 h and kept in a 65 °C-oven to measure constant weight. The soluble sugar content was measured following the anthrone colorimetric method .
Identification and classification of the AP2/ERF genes in Nelumbo nucifera
From a previous study conducted by Huang et al. , numerous SNPs of lotus genomic variation were obtained followed by selection based on the transcriptome difference between temperate and tropical lotus. Several SNPs were selected including SNP (megascaffold_20:3578539, C/T) in NnADAP. Nelumbo nucifera genome and protein sequence were analysed from the Lotus database [25, 32, 93] and Arabidopsis AP2/ERF gene sequences were retrieved from the DATF database (http://datf.cbi. pku.edu.cn, last accessed date: April 20, 2016)  and TAIR website, (https://www.arabidopsis.org/, last accessed date: April 25, 2016). According to the description of Nakano et al. (2006), the AP2/ERF sequence was determined based on the presence of the AP2 domain, and all putative AP2/ERF proteins were compared with Arabidopsis AP2/ERF proteins to classify them into different groups.
Acquisition of target SNP and verification of the sequence of NnADAP
We developed specific primers based on the target SNP to verify the transcriptome study of the target SNP of NnADAP by employing the sanger sequence. Simply, different varieties’ DNA was extracted as templates to conduct polymerase chain reaction with Primer STAR Max DNA Polymerase (Takara), 10 uM forward and reverse primers each. The PCR products were then ligated on T-vector for sequencing.
Phylogenetic and conserved motif analysis and gene structure of the AP2/ERF genes
Multiple alignments of Arabidopsis and lotus AP2/ERF protein sequences were carried out using ClustalW with default parameters (Thompson et al., 2002). Unrooted phylogenetic tree of all AP2/ERF proteins was generated with MEGA (V10.0) using the maximum likelihood (ML) method with the following parameters: Poisson correction, pairwise deletion, and 1000 bootstrap replicate . Conserved motifs in Nelumbo nucifera AP2/ERF TFs were identified using the motif finding tool MEME (Multiple EM for Motif Elicitation, V5.0.1 ) with the following parameters: optimum motif width ≥ 10 and ≤ 200, the maximum motifs numbers 25, and the occurrences of a single motif distributed among the sequences with zero or one per sequence (−mod zoops). The gene structure data was download from LOTUS.DB website (http://lotus-db.wbgcas.cn/) [25, 32, 93] and analysed using TBtools .
RNA isolation and quantitative real-time PCR
The plant samples were collected in different stages of stem development for the qRT-PCR assay (Fig. 5B). Total RNA extraction and cDNA synthesis and qRT-PCR reactions were described as previous study . A pair of primers (NnADAP forward and reverse primers, Additional file 1: Table S4) were designed to investigate the expression profiles of NnADAP. The candidate gene was performed in three biological replicates. Relative gene expression was standardized by the expression of lotus β-actin NNU_24864 (Additional file 1: Table S4) and analysed using the 2-∆∆CT method.
KASP SNP genotyping
The platform of LGC Limited of the KASP (Kompetitive Allele Specific PCR) genotyping technology was employed in the present study. The complete coding sequences of different locus alleles used for KASP assays were generated in lotus dataset . The diagnostics SNP was discovered and KASP primers designed following KASP guidelines carrying FAM and HEX (Additional file 1: Table S4) tails at the 3′ end. The KASP experiment steps were carried out as described in previous study . After initial testing on real-time PCR (BioRad®, CFX-96) system, the NnADAP KASP assay was applied to varieties.
Availability of data and materials
Nelumbo nucifera genome and protein sequence analysed during the current study are downloaded from the Lotus database [25, 32, 93] and Arabidopsis AP2/ERF gene sequences are available in the DATF database (http://datf.cbi. pku.edu.cn)  and TAIR website (https://www.arabidopsis.org/). The information of lotus genome was downloaded from the NCBI Sequence Read Archive SRA with the accession number SRP021228 (http://trace.ncbi.nlm.nih.gov/Traces/sra/?study=SRP021228) and the transcriptome dataset analysed during the current study are available in NCBI SRA (http://www.ncbi.nlm.nih.gov/Traces/sra) with accession number SRA271278 from previous studies [4, 33].
APETALA2/Ethylene Responsive Factor
Quantitative trait loci
Single nucleotide polymorphism
Gene Structure Display Server
Multiple EM for Motif Elicitation
Kompetitive Allele-Specific PCR
Relative Fluorescence Units
Insertions and deletions
Simple sequence repeat
Zhang XY, Wang QC. Preliminary study of the eco-types of genetic resources of tropical Lotus. Chin Landscape Arch. 2006;(7):82-5.
Wang QC, Zhang XY. Lotus development prospect —from Chinese perspective. Chin Landscape Arch. 2011;27(1):50–3.
Cheng L, Li S, Yin J, Li L, Chen X. Genome-wide analysis of differentially expressed genes relevant to rhizome formation in Lotus root (Nelumbo nucifera Gaertn). PLoS One. 2013;8(6):e67116.
Yang M, Zhu L, Pan C, Xu L, Liu Y, Ke W, Yang P. Transcriptomic analysis of the regulation of rhizome formation in temperate and tropical lotus (Nelumbo nucifera). Sci Rep. 2015;5(4):13059.
Nakano T, Suzuki K, Fujimura T, Shinshi H. Genome-wide analysis of the ERF gene family in Arabidopsis and Rice. Plant Physiol. 2006;140(2):411–32.
Duan C, Argout X, Gebelin V, Summo M, Dufayard JF, Leclercq J, Kuswanhadi, Piyatrakul P, Pirrello J, Rio M, et al. Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing. BMC Genomics. 2013;14:30.
Wu ZJ, Li XH, Liu ZW, Li H, Wang YX, Zhuang J. Transcriptome-based discovery of AP2/ERF transcription factors related to temperature stress in tea plant (Camellia sinensis). Funct Integr Genomics. 2015;15(6):741–52.
Song X, Li Y, Hou X. Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genomics. 2013;14(1):573.
Sakuma Y, Liu Q, Dubouzet JG, Abe H, Shinozaki K, Yamaguchi-Shinozaki K. DNA-binding specificity of the ERF/AP2 domain of Arabidopsis DREBs, transcription factors involved in dehydration- and cold-inducible gene expression. Biochem Biophys Res Commun. 2002;290(3):998–1009.
Yamaguchi-Shinozaki K, Shinozaki K. A novel cis-acting element in an Arabidopsis gene is involved in responsiveness to drought, low-temperature, or high-salt stress. Plant Cell. 1994;6(2):251–64.
Thomashow MF. PLANT COLD ACCLIMATION: freezing tolerance genes and regulatory mechanisms. Annu Rev Plant Physiol Plant Mol Biol. 1999;50:571–99.
Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, et al. Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science (New York, NY). 2003;301(5633):653–7.
Hu YX, Wang YX, Liu XF, Li JY. Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development. Cell Res. 2004;14(1):8–15.
Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. AP2/ERF family transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819(2):86–96.
Licausi F, Ohme-Takagi M, Perata P. APETALA2/ethylene responsive factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 2013;199(3):639–49.
Hiss M, Laule O, Meskauskiene RM, Arif MA, Decker EL, Erxleben A, Frank W, Hanke ST, Lang D, Martin A, et al. Large-scale gene expression profiling data for the model moss Physcomitrella patens aid understanding of developmental progression, culture and stress conditions. Plant J. 2014;79(3):530–9.
Huang Y, Zhang B, Sun S, Xing G, Wang F, Li M, Tian Y, Xiong A. AP2/ERF transcription factors involved in response to tomato yellow leaf curly virus in tomato. Plant Genome. 2016;9(2):0.
Ji-Yu Z, Qing-Ju W, Zhong-Ren G. Progresses on plant AP2/ERF transcription factors. HEREDITAS. 2012;34(7):835–47.
Huang Q, Sun M, Yuan T, Wang Y, Shi M, Lu S, Tang B, Pan J, Wang Y, Kai G. The AP2/ERF transcription factor SmERF1L1 regulates the biosynthesis of tanshinones and phenolic acids in Salvia miltiorrhiza. Food Chem. 2019;274:368–75.
Baud S, Mendoza MS, To A, Harscoët E, Lepiniec L, Dubreucq B. WRINKLED1 specifies the regulatory action of LEAFY COTYLEDON2 towards fatty acid metabolism during seed maturation in Arabidopsis. Plant J. 2007;50(5):825–38.
To A, Joubès J, Barthole G, Lécureuil A, Scagnelli A, Jasinski S, Lepiniec L, Baud S. WRINKLED Transcription Factors Orchestrate Tissue-Specific Regulation of Fatty Acid Biosynthesis in Arabidopsis. Plant Cell. 2012;24(12):5007.
Cernac A, Benning C. WRINKLED1 encodes an AP2/EREB domain protein involved in the control of storage compound biosynthesis in Arabidopsis. Plant J. 2004;40(4):575–85.
Fu F-F, Xue H-W. Coexpression analysis identifies Rice starch Regulator1, a Rice AP2/EREBP family transcription factor, as a novel Rice starch biosynthesis regulator. Plant Physiol. 2010;154(2):927.
Li H, Xiao Q, Zhang C, Du J, Li X, Huang H, Wei B, Li Y, Yu G, Liu H, et al. Identification and characterization of transcription factor ZmEREB94 involved in starch synthesis in maize. J Plant Physiol. 2017;216:11–6.
Borgi W, Ghedira K, Chouchane N. Antiinflammatory and analgesic activities of Zizyphus lotus root barks. Fitoterapia. 2007;78(1):16–9.
Collard BCY, Jahufer MZZ, Brouwer JB, Pang ECK. An introduction to markers, quantitative trait loci (QTL) mapping and marker-assisted selection for crop improvement: the basic concepts. Euphytica. 2005;142:169–96.
Yang M, Han Y, VanBuren R, Ming R, Xu L, Han Y, Liu Y. Genetic linkage maps for Asian and American lotus constructed using novel SSR markers derived from the genome of sequenced cultivar. BMC Genomics. 2012;13:653.
Zhang Q, Li L, VanBuren R, Liu Y, Yang M, Xu L, Bowers JE, Zhong C, Han Y, Li S, et al. Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map. BMC Genomics. 2014;15(1):372.
Pan L, Xia Q, Quan Z, Liu H, Ke W, Ding Y. Development of novel EST-SSRs from sacred lotus (Nelumbo nucifera Gaertn) and their utilization for the genetic diversity analysis of N. nucifera. J Hered. 2010;101(1):71–82.
Konishi S, Izawa T, Lin SY, Ebana K, Fukuta Y, Sasaki T, Yano M. An SNP caused loss of seed shattering during rice domestication. Science (New York, NY). 2006;312(5778):1392–6.
Win KT, Yamagata Y, Doi K, Uyama K, Nagai Y, Toda Y, Kani T, Ashikari M, Yasui H, Yoshimura A. A single base change explains the independent origin of and selection for the nonshattering gene in African rice domestication. New Phytol. 2017;213(4):1925–35.
Huang L, Yang M, Li L, Li H, Yang D, Shi T, Yang P. Whole genome re-sequencing reveals evolutionary patterns of sacred lotus (Nelumbo nucifera). J Integr Plant Biol. 2017;60(1):2–15.
Wang K, Deng J, Damaris RN, Yang M, Xu L, Yang P. LOTUS-DB: an integrative and interactive database for Nelumbo nucifera study. Database (Oxford). 2015;2015:bav023.
Berardini TZ, Reiser L, Li D, Mezheritsky Y, Muller RJ, Strait E, Huala E. The arabidopsis information resource: making and mining the “gold standard” annotated reference plant genome. Genesis. 2015;53(8):474–85.
Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8(3):275–82.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547–9.
Hu L, Liu S. Genome-wide identification and phylogenetic analysis of the ERF gene family in cucumbers. Genet Mol Biol. 2011;34(4):624–34.
Guo A-Y, Hu B, Gao G, Zhang H, Luo J, Jin J. GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics. 2014;31(8):1296–7.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–8.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools - an integrative toolkit developed for interactive analyses of big biological data. Molecular Plant. 2020;13(8):1194-202.
Li J, Shi T, Huang L, He D, Nyong’A TM, Yang P. Systematic transcriptomic analysis provides insights into lotus (Nelumbo nucifera) seed development. Plant Growth Regul. 2018;86(3):339–50.
Rasheed A, Wen W, Gao F, Zhai S, Jin H, Liu J, Guo Q, Zhang Y, Dreisigacker S, Xia X, et al. Development and validation of KASP assays for genes underpinning key economic traits in bread wheat. Theor Appl Genet. 2016;129(10):1843–60.
Cao D, Damaris RN, Zhang Y, Liu M, Li M, Yang P. Proteomic analysis showing the signaling pathways involved in the rhizome enlargement process in Nelumbo nucifera. BMC Genomics. 2019;20(1):766.
Cai XT, Xu P, Zhao PX, Liu R, Yu LH, Xiang CB. Arabidopsis ERF109 mediates cross-talk between jasmonic acid and auxin biosynthesis during lateral root formation. Nat Commun. 2014;5:5833.
Xiao YY, Chen JY, Kuang JF, Shan W, Xie H, Jiang YM, Lu WJ. Banana ethylene response factors are involved in fruit ripening through their interactions with ethylene biosynthesis genes. J Exp Bot. 2013;64(8):2499–510.
Varkonyi-Gasic E, Lough RH, Moss SM, Wu R, Hellens RP. Kiwifruit floral gene APETALA2 is alternatively spliced and accumulates in aberrant indeterminate flowers in the absence of miR172. Plant Mol Biol. 2012;78(4–5):417–29.
Salemme M, Sica M, Iazzetti G, Gaudio L, Aceto S. The AP2-like gene OitaAP2 is alternatively spliced and differentially expressed in inflorescence and vegetative tissues of the orchid Orchis italica. PLoS One. 2013;8(10):e77454.
Lee Y-S, Lee D-Y, Cho L-H, An G. Rice miR172 induces flowering by suppressing OsIDS1 and SNB, two AP2 genes that negatively regulate expression of Ehd1 and florigens. Rice (New York, NY). 2014;7(1):31.
Ito Y, Katsura K, Maruyama K, Taji T, Kobayashi M, Seki M, Shinozaki K, Yamaguchi-Shinozaki K. Functional analysis of rice DREB1/CBF-type transcription factors involved in cold-responsive gene expression in transgenic rice. Plant Cell Physiol. 2006;47(1):141–53.
Qin F, Kakimoto M, Sakuma Y, Maruyama K, Osakabe Y, Tran LS, Shinozaki K, Yamaguchi-Shinozaki K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 2007;50(1):54–69.
Ramaiah M, Jain A, Raghothama KG. Ethylene response Factor070 regulates root development and phosphate starvation-mediated responses. Plant Physiol. 2014;164(3):1484–98.
Kosmacz M, Parlanti S, Schwarzländer M, Kragler F, Licausi F, Van Dongen JT. The stability and nuclear localization of the transcription factor RAP2.12 are dynamically regulated by oxygen concentration. Plant Cell Environ. 2015;38(6):1094–103.
Liu Z, Zhu H, Liu Y, Kuang J, Zhou K, Liang F, Liu Z, Wang D, Ke W. Construction of a high-density, high-quality genetic map of cultivated lotus (Nelumbo nucifera) using next-generation sequencing. BMC Genomics. 2016;17:466.
Lin Z, Zhang C, Cao D, Damaris RN, Yang P. The Latest Studies on Lotus (Nelumbo nucifera)-an Emerging Horticultural Model Plant. Int J Mol Sci. 2019;20(15):3680.
Ham YK, Hwang KE, Song DH, Kim YJ, Shin DJ, Kim KI, Lee HJ, Kim NR, Kim CJ. Lotus (Nelumbo nucifera) rhizome as an antioxidant dietary Fiber in cooked sausage: effects on physicochemical and sensory characteristics. Korean J Food Sci Anim Resour. 2017;37(2):219–27.
Zhao W, Xie W, Du S, Yan S, Li J, Wang Q. Changes in physicochemical properties related to the texture of lotus rhizomes subjected to heat blanching and calcium immersion. Food Chem. 2016;211:409–14.
Sun S, Zhang G, Ma C. Preparation, physicochemical characterization and application of acetylated lotus rhizome starches. Carbohydr Polym. 2016;135:10–7.
Cheng L, Liu X, Yin J, Yang J, Li Y, Hui L, Li S, Li L. Activity and expression of ADP-glucose pyrophosphorylase during rhizome formation in lotus (Nelumbo nucifera Gaertn.). Bot Stud. 2016;57(1):26.
Lin L, Huang J, Zhao L, Wang J, Wang Z, Wei C. Effect of granule size on the properties of lotus rhizome C-type starch. Carbohydr Polym. 2015;134:448–57.
Zhang S, Zhu C, Lyu Y, Chen Y, Zhang Z, Lai Z, Lin Y. Genome-wide identification, molecular evolution, and expression analysis provide new insights into the APETALA2/ethylene responsive factor (AP2/ERF) superfamily in Dimocarpus longan Lour. BMC Genomics. 2020;21(1):62.
Licausi F, Giorgi FM, Zenoni S, Osti F, Pezzotti M, Perata P. Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera. BMC Genomics. 2010;11:719.
Zhuang J, Cai B, Peng R-H, Zhu B, Jin X-F, Xue Y, Gao F, Fu X-Y, Tian Y-S, Zhao W, et al. Genome-wide analysis of the AP2/ERF gene family in Populus trichocarpa. Biochem Biophys Res Commun. 2008;371(3):468–74.
Liu M, Sun W, Ma Z, Zheng T, Huang L, Wu Q, Zhao G, Tang Z, Bu T, Li C, et al. Genome-wide investigation of the AP2/ERF gene family in tartary buckwheat (Fagopyum Tataricum). BMC Plant Biol. 2019;19(1):84.
Jofuku KD, den Boer BG, Van Montagu M, Okamuro JK. Control of Arabidopsis flower and seed development by the homeotic gene APETALA2. Plant Cell. 1994;6(9):1211–25.
Mathieu J, Yant LJ, Murdter F, Kuttner F, Schmid M. Repression of flowering by the miR172 target SMZ. PLoS Biol. 2009;7(7):e1000148.
Golembeski GS, Imaizumi T. Photoperiodic regulation of Florigen function in Arabidopsis thaliana. Arabidopsis Book. 2015;13:e0178.
Zheng XF, You YN, Diao Y, Zheng XW, Xie KQ, Zhou MQ, Hu ZL, Wang YW. Development and characterization of genic-SSR markers from different Asia lotus (Nelumbo nucifera) types by RNA-seq. Genet Mol Res. 2015;14(3):11171–84.
Yant L, Mathieu J, Dinh TT, Ott F, Lanz C, Wollmann H, Chen X, Schmid M. Orchestration of the floral transition and floral development in Arabidopsis by the bifunctional transcription factor APETALA2. Plant Cell. 2010;22(7):2156–70.
Bergonzi S, Albani MC, Ver Loren van Themaat E, Nordström KJ, Wang R, Schneeberger K, Moerland PD, Coupland G. Mechanisms of age-dependent response to winter temperature in perennial flowering of Arabis alpina. Science (New York, NY). 2013;340(6136):1094–7.
Chailakhyan MK, Yanina LI, Devedzhyan AG, Lotova GN. Photoperiodism and tuber formation in grafting of tobacco onto potato. Dokl Bot Sci. 1981;257:1276–80.
Jackson SD. Multiple signaling pathways control tuber induction in potato. Plant Physiol. 1999;119(1):1.
Turck F, Fornara F, Coupland G. Regulation and identity of Florigen: FLOWERING LOCUS T moves center stage. Annu Rev Plant Biol. 2008;59(1):573–94.
Debbarma J, Sarki YN, Saikia B, Boruah HPD, Singha DL, Chikkaputtaiah C. Ethylene response factor (ERF) family proteins in abiotic stresses and CRISPR–Cas9 genome editing of ERFs for multiple abiotic stress tolerance in crop plants: a review. Mol Biotechnol. 2019;61(2):153–72.
Gilmour SJ, Sebolt AM, Salazar MP, Everard JD, Thomashow MF. Overexpression of the Arabidopsis CBF3 transcriptional activator mimics multiple biochemical changes associated with cold acclimation. Plant Physiol. 2000;124(4):1854–65.
Abiri R, Shaharuddin NA, Maziah M, Yusof ZNB, Atabaki N, Sahebi M, Valdiani A, Kalhori N, Azizi P, Hanafi MM. Role of ethylene and the APETALA 2/ethylene response factor superfamily in rice under various abiotic and biotic stress conditions. Environ Exp Bot. 2017;134:33–44.
Hsieh TH, Lee JT, Yang PT, Chiu LH, Charng YY, Wang YC, Chan MT. Heterology expression of the Arabidopsis C-repeat/dehydration response element binding factor 1 gene confers elevated tolerance to chilling and oxidative stresses in transgenic tomato. Plant Physiol. 2002;129(3):1086–94.
Jin L-G, Liu J-Y. Molecular cloning, expression profile and promoter analysis of a novel ethylene responsive transcription factor gene GhERF4 from cotton (Gossypium hirstum). Plant Physiol Biochem. 2008;46(1):46–53.
Fu M, Kang HK, Son SH, Kim SK, Nam KH. A subset of Arabidopsis RAV transcription factors modulates drought and salt stress responses independent of ABA. Plant Cell Physiol. 2014;55(11):1892–904.
Shin HY, Nam KH. RAV1 negatively regulates seed development by directly repressing MINI3 and IKU2 in Arabidopsis. Mol Cell. 2018;41(12):1072–80.
Li CW, Su RC, Cheng CP, Sanjaya, You SJ, Hsieh TH, Chao TC, Chan MT. Tomato RAV transcription factor is a pivotal modulator involved in the AP2/EREBP-mediated defense pathway. Plant Physiol. 2011;156(1):213–27.
Zhao SP, Xu ZS, Zheng WJ, Zhao W, Wang YX, Yu TF, Chen M, Zhou YB, Min DH, Ma YZ, et al. Genome-wide analysis of the RAV family in soybean and functional identification of GmRAV-03 involvement in salt and drought stresses and exogenous ABA treatment. Front Plant Sci. 2017;8:905.
Yan J, Yang X, Shah T, Sánchez-Villeda H, Li J, Warburton M, Zhou Y, Crouch JH, Xu Y. High-throughput SNP genotyping with the GoldenGate assay in maize. Mol Breed. 2010;25(3):441–51.
Hyten DL, Cannon SB, Song Q, Weeks N, Fickus EW, Shoemaker RC, Specht JE, Farmer AD, May GD, Cregan PB. High-throughput SNP discovery through deep resequencing of a reduced representation library to anchor and orient scaffolds in the soybean whole genome sequence. BMC Genomics. 2010;11:38.
Wang B, Tan H-W, Fang W, Meinhardt LW, Mischke S, Matsumoto T, Zhang D. Developing single nucleotide polymorphism (SNP) markers from transcriptome sequences for identification of longan (Dimocarpus longan) germplasm. Horticulture Res. 2015;2:14065.
Hu J, Gui S, Zhu Z, Wang X, Ke W, Ding Y. Genome-wide identification of SSR and SNP markers based on whole-genome re-sequencing of a Thailand wild sacred Lotus (Nelumbo nucifera). PLoS One. 2015;10(11):e0143765.
Zhao M, Yang JX, Mao TY, Zhu HH, Xiang L, Zhang J, Chen LQ. Detection of Highly Differentiated Genomic Regions Between Lotus (Nelumbo nucifera Gaertn.) With Contrasting Plant Architecture and Their Functional Relevance to Plant Architecture. Front Plant Sci. 2018;9:1219.
Yang M, Fu J, Xiang QY, Liu YL. The Core-collection construction of flower Lotus based on AFLP molecular markers. Sci Agric Sin. 2011;44(15):3193–205.
Masuda J-i, Ozaki Y, Okubo H. Rhizome transition to storage organ is under phytochrome control in lotus (Nelumbo nucifera). Planta. 2007;226(4):909–15.
Masuda JI, Yoshimizu S, Ozaki Y, Okubo H. Rhythmic response of rhizome growth to light-break in Lotus (Nelumbo nucifera). J Fac Agric Kyushu Univ. 2007a;52(1):35–8.
Shen-Miller J. Sacred lotus, the long-living fruits of China antique. Seed Sci Res. 2002;12(3):131–43.
Chen WP. The structure,elongation and thickening of rhizome in Nelumbo nucifera Gaertn. J Integr Plant Biol. 1989;31(3):191–7.
McKenzie M, Corrigan V. Chapter 12 - Potato Flavor. In: Singh J, Kaur L, editors. Advances in Potato Chemistry and Technology (Second Edition). San Diego: Academic Press; 2016. p. 339–68.
Ming R, VanBuren R, Liu Y, Yang M, Han Y, Li L-T, Zhang Q, Kim M-J, Schatz MC, Campbell M, et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.). Genome Biol. 2013;14(5):R41.
Cao D, Lin Z, Huang L, Damaris RN, Li M, Yang P. A CONSTANS-LIKE gene of Nelumbo nucifera could promote potato tuberization. Planta. 2021;253(3):65.
Semagn K, Babu R, Hearne S, Olsen M. Single nucleotide polymorphism genotyping using Kompetitive allele specific PCR (KASP): overview of the technology and its application in crop improvement. Mol Breed. 2014;33(1):1–14.
The authors thank Wuhan Botanical Garden, Chinese Academy of Sciences for the help in germplasm sharing and Han Liu for the help in proof-reading and revising the manuscript.
This study was supported by the distinguished Professorship project to Pingfang Yang from Hubei University. The funding body has no roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests, but the corresponding author is on the editorial board of BMC Genomics.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The AP2/ERF family in Nelumbo nucfera. Table S2. The motif sequence of lotus AP2/ERF family. Table S3. The transcriptome data of AP2/ERF family in the temperate lotus and tropical lotus rhizome. Table S4. Primers used in this study. Table S5. Detection and the authentication of SNPs in genome sequencing results through first-generation sequencing. Table S6. The rhizome traits and SNP phenotype of 95 lotus germplasm resources
A, Mapping of AP2/ERF genes in lotus main megascaffolds. The unit of the length of each megascaffold is megabase. B, The numbers of these genes in each megascaffold
About this article
Cite this article
Cao, D., Lin, Z., Huang, L. et al. Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology. BMC Genomics 22, 171 (2021). https://doi.org/10.1186/s12864-021-07473-w
- Nelumbo nucifera