Heyun Song1,2
Yanling Liu1,3
Gangqiang Dong4Minghua Zhang1,2Yuxin Wang1,2Jia Xin1,2Yanyan Su4Heng Sun1,3*
Mei Yang1,3*- 1Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, China
- 2University of Chinese Academy of Sciences, Beijing, China
- 3Aquatic Plant Research Center, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan, China
- 4Amway (China) Botanical R&D Centre, Wuxi, China
NAC (NAM, ATAF, and CUC) is a ubiquitously expressed plant-specific transcription factor (TF) family which is involved in the regulation of various biological processes. However, a systematic characterization of NAC gene family is yet to be reported in lotus. Here, 82 NnNAC genes which included five predicted membrane-bound NAC proteins were identified in the lotus genome. Phylogenetic analysis revealed seven-subfamily clusters (I–VII) of NnNAC proteins, with homologous gene pairs displaying similar conserved motifs and gene structure characteristics. Transactivation assay of NnNAC proteins revealed an extensive transcriptional activation capacity which is mediated by the highly divergent C-terminal activation domain (AD). Expression analysis of NnNAC genes in lotus tissues showed high transcript levels in root, stamen, petal and seed coat. In addition, 30 and 29 differentially expressed NnNAC candidate genes putatively involved in lotus seed development and response to complete submergence stress, respectively, were identified. Overall, our study provides potentially useful candidate gene resources for future molecular breeding of lotus varieties with novel agronomic traits.
Introduction
Transcription factors (TFs) are master regulators of gene expression. NAC is one of the largest and most characterized plant-specific TF superfamily with 117 and 152 members in the dicotyledonous Arabidopsis thaliana and Glycine max, respectively, and 151 members in the monocotyledonous Oryza sativa (Nuruzzaman et al., 2010; Dung et al., 2011; Yan et al., 2017; Li et al., 2018a). NAC protein generally harbors a highly-conserved N-terminal DNA-binding domain (BD) and a highly variable C-terminal transcriptional regulatory region (TRR) (Puranik et al., 2012). The DNA binding of NAC TFs is determined by the presence of NAC domains, which are classified into five subdomains designated, A to E. The highly conserved C and D subdomains are necessary for DNA interaction, subdomain A is necessary for protein dimerization or heterodimerization, whereas the highly variable subdomains B and E determine the functional diversity of NAC proteins (Puranik et al., 2012). The highly divergent C-terminal TRRs are responsible for transcriptional regulation (Christiansen and Gregersen, 2014; Lindemose et al., 2014). Most NAC TFs are located in the cell nucleus; however, some membrane-associated NAC proteins show subcellular localization at the endoplasmic reticulum or plasma membrane, which could be attributed to the α-helical transmembrane (TM) motif in the C-terminus (Seo et al., 2010; Li et al., 2016). Membrane-bound NAC genes have previously been reported in other species, such as A. thaliana (Arabidopsis), O. sativa (rice), and S. lycopersicum (tomato) with 13, 5, and 13 gene members, respectively (Kim et al., 2010a; Bhattacharjee et al., 2017). Membrane-bound NAC genes have important regulatory functions in stress responses and are strictly post-translationally regulated under specific conditions (Dung et al., 2011).
The first reported NAC gene developed no shoot apical meristem (SAM) in Petunia embryos and showed an equiformal phenotype in SAM formation in Arabidopsis mutants (Souer et al., 1996; Takada et al., 2001). The NAC TF family has been associated with various biological processes, including tissue development, hormone response, organogenesis, secondary cell wall (SCW) biosynthesis, and stress response (Zhong et al., 2006; Yamaguchi et al., 2010; Nakano et al., 2015; Yuan et al., 2020). In addition to regulating these processes, numerous NAC TFs have been shown to modulate seed development in some plant species. For example, three seed-size related NAC genes, ONAC020, ONAC023, and ONAC026 with varied expression levels among rice accessions were highly expressed during seed development (Mathew et al., 2016). Similarly, ONAC025, ONAC127, and ONAC129 were found to be involved in seed development by modulating rice grain filling during reproductive period (Mathew et al., 2020; Ren et al., 2021). Plants constantly encounter diverse biotic and abiotic stresses throughout their lives, which cause adverse effects on their growth and development. Increasing evidence suggests that NAC proteins play pivotal roles in abiotic stress tolerance (Puranik et al., 2012; Shao et al., 2015; Yuan et al., 2020; Singh et al., 2021). In Arabidopsis, ANAC019, ANAC055, and ANAC072 are transcriptional activators of drought stress response, and their overexpression improved drought tolerance in transgenic plants (Tran et al., 2004). Similarly, a NAC SlJUB1 gene, was shown to enhance tomato drought tolerance by activating the expression of several stress related genes, such as SlDREB1, SlDREB, and SlDELLA (Thirumalaikumar et al., 2018). Due to their roles in regulating important biological processes, comprehensive identification of stress related NAC genes is thus crucial for studying stress response mechanisms in plants.
Lotus (Nelumbo nucifera Gaertn.) is an old domesticated perennial wetland plant in the family Nelumbonaceae, which contains a single genus, Nelumbo, with two extant species: N. nucifera Gaertn. and N. lutea Pers. (Ming et al., 2013; Sun et al., 2021). The lotus seeds, rhizomes, and flowers have versatile uses, including as popular vegetable and ornamental plant with diverse medicinal properties (Sun et al., 2020). Lotus is popularly cultivated in Asia, and it is mainly categorized as seed-, rhizome-, and flower-lotus based on agronomic traits and distinct uses (Yang et al., 2015). Seed-lotus is predominantly cultivated for its edible lotus seeds, which are rich in proteins, vitamins, minerals, essential amino acids, and other numerous health-promoting bioactive components (Li et al., 2018b; Sun et al., 2021). Rhizome-lotus is produced for its edible underground stems, while flower-lotus is popular for its diversely colored and shaped flowers (Yang et al., 2015). TFs, including MYB, WRKY, bHLH, ERF, and bZIP have previously been reported to play crucial roles in lotus growth and development (Chen et al., 2013; Deng et al., 2016; Li et al., 2018b; Li et al., 2019; Sun et al., 2020). For example, Deng et al. (2016) demonstrated the NnMYB5-mediated transcriptional activation of anthocyanin synthesis, and its overexpression induced higher anthocyanin accumulation in immature transgenic Arabidopsis seeds and flower stalks. In addition, NnANT and NnAP2 have been reported as potential negative regulators of lotus seed size and development (Li et al., 2018b). Moreover, the lotus NAC genes have been implicated in root browning under anaerobic stress and the formation of adventitious roots (Cheng et al., 2019; Min et al., 2019). However, despite previous reports, comprehensive mapping of genome-wide NAC TF family, which could potentially improve our knowledge on the gene family function, distribution, and evolution in sacred lotus is still lacking.
The recently improved sacred lotus genome assembly by Shi et al. (2020) provides an excellent opportunity to investigate the genome-wide distribution and evolution of NAC gene family in lotus. This study aimed to conduct a comprehensive analysis to clarify sequence features, phylogenesis, genome synteny, expression patterns, subcellular localization, and transcriptional activation capacity of NAC gene family members in sacred lotus. The study also integrated the analysis of possible roles of NnNAC genes in lotus seed development and response to complete submergence stress. Our results will not only improve the understanding on NnNAC gene family function, genome distribution, and evolution, but also provides vital candidate genes for future molecular breeding in lotus.
Materials and Methods
Identification and Phylogenetic Analysis of NAC Genes in the Sacred Lotus Genome
To identify NnNAC genes, we retrieved the gene annotation gff3 file containing all the predicted protein sequences from lotus reference genome (N. Gaertn.) (Li et al., 2021a). The NnNAC genes were predicted using PlantTFDB 4.0 (http://planttfdb.cbi.pku.edu.cn/), and their complete amino acid sequences were confirmed with Clustalx v.1.81 and MEME software (http://meme-suite.org/index.html). Predicted molecular weights (Mw) and isoelectric points (pI) of NnNAC proteins were calculated using ExPASy software (http://web.expasy.org/protparam/). Membrane-bound NnNAC members were predicted using the TMHMM server v.2.0 (http://www.cbs.dtu.dk/services/TMHMM/).
Multiple sequence alignments of the identified 82 NnNAC proteins were carried out using Clustalx v.1.81 with default settings. Phylogenetic tree was constructed in MEGA7 using the Neighbor-Joining algorithm, the evolutionary distances were computed with the Poisson correction method, and tree nodes were evaluated with 1,000 bootstrap replications (Kumar et al., 2016).
Gene Structure, Conserved Motifs, Promoter, Chromosomal Location and Genome Scale Syntenic Analysis of NnNAC Genes
To obtain the intron and exon structures, the coding and genomic sequences of NnNAC genes were analyzed with the Gene Structure Display Server (GSDS v.2.0) online program (http://gsds.cbi.pku.edu.cn/). Conserved motifs were identified using the MEME online tool (http://meme-suite.org/index.html) with default settings. For promoter sequence analysis, a 2-kb length sequence upstream of the start codon of each gene was downloaded for cis-elements prediction using PlantCARE program (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (Lescot et al., 2002). Distribution of NnNAC genes within the lotus chromosomes was visualized using TBtools software (Chen et al., 2020). NAC genes were sequentially named from chromosome (Chr) 1 to 8 with numerical digits in the order of their physical location from top to bottom of the chromosomes.
Genome sequences of A. thaliana and O. sativa were downloaded from Phytozome v13 (https://phytozome-next.jgi.doe.gov/). Genome scale syntenic analysis between the three plant species was conducted by MCScan X, and TBtools software was used to visualize their relationship (Wang et al., 2012; Chen et al., 2020).
Expression Profiling of NnNAC Genes
To investigate the expression patterns of NnNAC genes in lotus tissues, the transcriptome data corresponding to gene expression abundance in various tissues were obtained from Nelumbo Genome Database (http://nelumbo.biocloud.net/nelumbo/home). Public transcriptome datasets on seed development and response to complete submergence stress were retrieved from the National Center for Biotechnology Information (NCBI) with accession numbers, SRP127765 and PRJNA723672, respectively (Li et al., 2018b; Deng et al., 2022). We performed quality control of the downloaded transcriptome data sets with FastQC and Trimmomatics programs (Bolger et al., 2014; Wingett and Andrews, 2018). Subsequently, gene expression levels were quantified by fragments per kilobase of transcript per million mapped fragments (FPKM) using StringTie software (Pertea et al., 2015). Gene expression was visualized using TBtools software. DESeq2 R package v.1.10.1 was used for pairwise comparisons to identify differentially expressed genes (DEGs) with filter the criteria, Fold Change (FC) ≥ 2 and False Discovery Rate (FDR) < 0.05 (Love et al., 2014).
qRT-PCR Analysis
High-quality RNAs were reverse transcribed to cDNA using TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (Lot#M31212, Beijing TransGene Biotech, Beijing, China). qRT-PCR analysis was carried out on a StepOnePlus™ Real-time PCR System (Applied Biosystems, United States). The relative gene expression level was calculated and normalized using NnACTIN (Gene ID NNU_24864) used as the internal standard. All primers used for qRT-PCR were listed in Supplementary Table S1.
Subcellular Localization of NnNAC Genes
ProtComp v9.0 (http://linux1.softberry.com/berry.phtml) was used to obtain the subcellular localization predictions of NnNAC proteins. Coding regions of NnNAC genes were cloned into entry vectors (pDONR221 Zeo) using high fidelity primers (Supplementary Table S1) following the BP-clonase kit instruction manual (Lot#2335893, Invitrogen by Thermo Fisher Scientific, United States). Transformed plasmids were then cloned into PMDC43 vectors with LR-clonase according to the manufacturer’s instructions. These plasmids were subsequently inserted into the Agrobacterium tumefaciens cv. GV3101 using electric shock method. Young leaves of tobacco plants (4–6 weeks-old) were selected for transformation assay according to the methods of Kokkirala et al. (2010). A confocal laser scanning microscope (Leica TCS SP2; Leica microsystems, Wetzlar, Germany) was used to photograph the agroinfiltrated leaves 48 h after infiltration.
Transcriptional Activation Analysis
Transactivation capacity of NnNAC proteins was assessed using yeast transactivation assay (Sun et al., 2018). The full-length coding sequences of NnNAC genes and the truncated NnNAC-AD and NnNAC-BD forms were independently cloned in pGBKT7 vectors. The resulting plasmids were subsequently transformed into Y2H gold strains according to manufacturer’s instructions (Cat.No.630489, Clontech) then plated on SD growth mediums without Trp, but with X-α-Gal (SD-Trp + X-α-Gal) to identify the transactivation activity of NnNAC proteins. SD mediums were photographed after 36 h incubation. Primers used for this assay were listed in Supplementary Table S1.
Results
Genome-Wide Identification of NAC Genes in the Lotus Genome and Their Gene Structural and Conserved Motifs Analyses
A total of 82 complete non-redundant NnNAC genes were identified in the lotus genome (N. nucifera Gaertn.) from PlantTFDB (Supplementary Table S2). Multiple sequence alignment revealed the presence of conserved NAC domain and the five (A–E) subdomains among all predicted NnNAC proteins (Figures 1A, 2B). Five putative membrane-bound NAC proteins each harboring a single predicted TM at the C-terminal were identified, including NnNAC6, NnNAC12, NnNAC23, NnNAC27, and NnNAC44 (Figure 1A; Supplementary Figure S1; Supplementary Table S2). Analysis of NnNAC protein features revealed that their amino acid residue, molecular weight (Mw), and protein isoelectric point (pI) ranged from 152 to 683, 7.63 to 77.94 KDa, and 4.48 to 9.92, respectively (Figure 1B; Supplementary Table S3).
FIGURE 1. NAC domains, protein sequence features and phylogenetic analysis in lotus. (A) Schematic representation of the highly conserved NnNAC key domain and subdomains A–E (purple squares) at the N-terminal. The blue square represents the transmembrane motif (TM). (B) Box plots showing statistical distribution of amino acid sequence length, molecular weight (Mw) and isoelectric point (pI) of NnNACs. (C) Phylogenetic clustering of 82 identified lotus NAC proteins. The purple triangle represents NnNAC protein harboring the TM motif.
FIGURE 2. Conserved motifs, cis-elements and NnNAC gene structure analyses. (A) Clusters and conserved motifs of NnNAC proteins. (B) Screening of conserved NAC domain in lotus. (C) Identification of cis-regulatory elements in the 2-kb upstream promoter regions of NnNAC genes. (D) The lotus NAC gene structures.
To understand the evolutionary relationships among the NnNAC gene family, a phylogenetic tree was constructed using the 82 identified protein sequences. All NnNAC proteins were clustered into seven subfamilies (I–VII) (Figure 1C). Notably, the five predicted membrane-bound proteins were exclusively present in I and II subfamily clusters (Figure 1C). Conserved motif analysis identified five conserved protein motifs, with motif 2 and motif 5 as the most conserved (Figure 2A). Interestingly, NnNACs within the same phylogenetic cluster tended to harbor similar conserved motifs, while paired NnNACs, such as NnNAC1 and NnNAC14, NnNAC30 and NnNAC38, NnNAC46 and NnNAC47, and NnNAC67 and NnNAC75 shared equiformal conserved motifs (Figure 2A).
Gene structure variability due to different exon and intron combinations is useful for understanding the diverse gene functions and genome evolution. To explore the structural diversity of NnNAC genes, the exon/intron arrangements were analyzed. The numbers of exons among NnNACs genes varied from one to seven, with most genes (64.2%) having three exons (Supplementary Figure S2). NnNAC38 had seven exons which was the highest number observed, while NnNAC55 and NnNAC81 only harbored one exon each (Figure 2D). Notably, genes sharing the same phylogenetic cluster showed similar conserved structural characteristics (Figure 2).
Chromosomal Location and Syntenic Analysis of NAC Genes in Lotus, Arabidopsis and Rice
Chromosomal localization revealed uneven distribution of NnNAC genes within the lotus genome with each of the eight lotus chromosomes harboring at least one gene. The highest gene density was observed in Chr 1 with 20 NnNACs, representing 24% of total genes identified, while only seven NnNAC genes were anchored on Chr 7. In addition, 8, 8, 12, 9, 10, and 8 NnNAC genes were anchored on Chr 2, 3, 4, 5, 6, and 8, respectively (Figure 3A; Supplementary Figure S3). Interestingly, some NnNACs indicated potential evidence of gene duplication events in the lotus genome (Supplementary Figure S3). For example, NnNAC18 and its homolog NnNAC19 located within a 0.77-Mb interval on Chr 1 shared about 99.8% identity (Supplementary Figure S4A). Similarly, four NnNAC69, NnNAC70, NnNAC71, and NnNAC72 anchored within a 0.41-Mb chromosomal region on Chr 7 shared about 98.8% amino acid sequence identity (Supplementary Figure S4B). In addition, 16 pairs of segmentally duplicated NnNAC genes were identified by intragenomic synteny analysis. Highly conserved amino acid sequences were observed among some collinear gene pairs, such as NnNAC54 and its paralog NnNAC60 located on Chr 6 which shared about 87.5% amino acid sequence identity, and NnNAC67 on Chr 6 which shared about 88.5% identity with NnNAC75 anchored on Chr 8 (Supplementary Figures S4C,D).
FIGURE 3. Genome scale synteny analysis of lotus, Arabidopsis, and rice NAC genes. (A) Syntenic relationship of NnNAC genes. All 82 NnNAC genes are marked according to their chromosome distribution in the lotus genome and syntenic gene pairs are connected with red lines. (B) Syntenic pairs of NAC genes between Arabidopsis and lotus. (C) Syntenic pairs of NAC genes between lotus and rice. The lotus, Arabidopsis, and rice chromosomes are shown in orange, green, and purple bars, respectively. Syntenic gene pairs are linked with red lines.
Genome synteny analysis is an effective method for inferring genome evolutionary history. The evolution and collinearity of NAC gene family was explored using the complete genome sequences of sacred lotus, Arabidopsis and rice. As a result, a total of 59 orthologous gene pairs sharing high similarities were identified between Arabidopsis to lotus (Figure 3B). In contrast, only 26 orthologous gene pairs were identified between rice to lotus genomes (Figure 3C). These results suggest that lotus, an ancient true dicot, is more likely to be evolutionarily closer to dicots than monocots.
Screening of NAC Gene Cis-Elements in Lotus
To study the transcriptional regulation of NnNAC genes, cis-elements were identified using the PLANTCARE database. As a result, numerous cis-regulatory elements linked to stress responses, development, and phytohormone responses were detected in the promoter regions of NnNAC genes. Of the 82 gene promoters screened, ARE element, CGTCA-motif, TGACG-motif, G-box, ABRE, and W box were the six most commonly identified cis-elements with percentage frequencies of 89, 84.1, 82.9, 76.8, 71.9, and 65.9, respectively (Figure 2C; Supplementary Figure S5). The W box and ABRE are both associated with drought and salt stress responses, while TGACG-motif and CGTCA-motif are involved in phytohormone response and development. These results suggest the likely roles of NnNAC genes as crucial stress response and development regulators in lotus.
Subcellular Localization and Transcriptional Activation Assay of NnNAC Proteins
Bioinformatic studies have predicted that about 74 (∼90%) of the identified NnNAC proteins were nuclear while others were extracellularly located (Supplementary Table S3). Three predicted nuclear located NnNAC genes, including NnNAC20, NnNAC23, and NnNAC40 and a membrane located NnNAC12 gene were selected for subcellular localization to test the prediction accuracy results. Microscopic observations showed nuclear localization of NnNAC20, NnNAC23 and NnNAC40, while the GFP signals of NnNAC12 were observed in plasma membranes (Figure 4A).
FIGURE 4. Subcellular localization and transcriptional activation assay of NnNAC proteins. (A) Subcellular localization of four NnNAC proteins in tobacco leaves. (B) Transactivation assay of NnNAC proteins in the Y2H system. I: NnNAC12; II: NnNAC23; III: NnNAC20; IV: NnNAC52; V: NnNAC40; AD and BD represent NnNAC protein with AD and BD domains, respectively; P: positive control (pGBKT7-53); N: negative control. (C) Schematic representation of NnNAC protein truncations.
Transcriptional activation assay is useful for determining the regulatory capacity of TFs. Transactivation capacities of NnNAC proteins were determined using the yeast two-hybrid system (Y2H). Full-length cDNA of five selected NAC proteins, including NnNAC12, NnNAC20, NnNAC23, NnNAC40, and NnNAC52 were individually cloned into GAL4DB vector. As a result, blue colonies indicating transcriptional activation activity were observed in all tested NnNAC proteins except NnNAC52 (Figure 4B). To investigate the transactivation domain, we further truncated the NnNAC proteins into conserved N-terminal DNA-binding domains (BD) and high divergent C-terminal activation domain (AD) regions (Figure 4C). Expectedly, transactivation activity was observed for the AD domain and not for BD domain (Figure 4B).
Expression Profiling of NnNAC Genes in Lotus Tissues
Published transcriptome data of leaf, petiole, petal, receptacle (immature and mature), stamen (immature and mature), carpel (unpollinated and pollinated), seed coat (6, 12, and 18 days after pollination, DAP), root and rhizome were obtained and used to investigate the tissue expression patterns of NnNAC genes (Figure 5A). Overall, 64 (∼78%) of the 82 NnNAC genes were expressed in at least one lotus tissue, whereas 14 (∼17%) were expressed in all tested tissues (FPKM ≥ 1). Low transcript abundance of FPKM < 5 in all tested tissues was observed in 42 (∼51%) NnNAC genes. Conversely, the expressions of NnNAC25, NnNAC43, NnNAC46, NnNAC47, NnNAC71, and NnNAC72 were not detected in all analyzed tissues (Supplementary Table S4). Notably, highly expressed NnNAC genes with FPKM ≥ 20 were identified in root, stamen, petal, and seed coat (Figure 5B; Supplementary Table S4). NnNAC genes, including NnNAC24, NnNAC31, NnNAC41, NnNAC48, and NnNAC50 exhibited tissue-specific expression in roots, while NnNAC27 and NnNAC42 were only expressed in root and stamen. In contrast, NnNAC20, NnNAC29, NnNAC37, NnNAC52, NnNAC53, and NnNAC66 were expressed in all tested tissues. Interestingly, phylogenetically close NnNAC genes, including NnNAC24/64, NnNAC31/39, and NnNAC44/45 exhibited comparable expression patterns. In addition, duplicated gene pairs, such as NnNAC18/19 and NnNAC69/70/71/72 also shared similar expression patterns.
FIGURE 5. Expression profiling of NnNAC genes in different lotus tissues. (A) Lotus tissues used to analyze the expression of NnNAC genes. (B) Heat-map clustering based on FPKM expression of 82 NnNAC genes in different lotus tissues.
Identification of NnNAC Genes Involved in Lotus Seed Development and Complete Submergence Stress Response
Multiple complex pathways respond to different plant signals in a coordinated manner that involves the recruitment of various transcription factors (Li et al., 2018b; Sun et al., 2020). The published transcriptome data of seed cotyledons at 9, 12 and 15 DAP in “China Antique” (CA) and “Jianxuan 17” (JX) were used to investigate the possible roles of NnNACs in lotus seed development. A total of 10 and 23 NnNAC DEGs were subsequently identified in CA and JX, respectively, while 20 NnNAC DEGs were identified between CA and JX. Overall, 30 NnNAC DEGs were identified during lotus seed development (Supplementary Table S5). Of these, 14 (46.7%) NnNACs were expressed in the cotyledon in at least one development stage (FPKM ≥ 1) (Figure 6A). Highly expressed NnNAC genes, including NnNAC12, NnNAC20, NnNAC37, NnNAC38, NnNAC48, and NnNAC52 with FPKM abundance ≥10, and with varied expression patterns during lotus seed development were identified, suggesting their likely roles in the regulation of lotus seed development. The expression of NnNAC38 showed a continuous increase from 9 to 15 DAP in CA, while it was downregulated in JX. In contrast, the expression profile of NnNAC37 showed a continuous increase from 9 to 15 DAP in both CA and JX. Notably, preferential expression of NnNAC20 and NnNAC52 gene pairs over other NnNACs was observed, with the latter showing a continuous decrease in expression from 9 to 15 DAP in CA. In contrast, the expression of NnNAC20 showed a continuous downregulation in CA, whereas it was upregulated in JX from 9 to 12 DAP. Four NnNAC genes, including NnNAC23, NnNAC37, NnNAC38, and NnNAC40 were randomly selected to evaluate the reliability of RNA-Seq by qRT-PCR. The results showed comparable expression changes with of transcriptome data (Supplementary Figure S6). For example, qRT-PCR and FPKM expression results showed that the levels of NnNAC37 were persistently elevated from 9 to 15 DAP in both CA and JX. NAC genes have previously been implicated in the regulation of various stress signaling pathways. As a perennial aquatic plant, lotus experience frequent submergence stress during its life cycle due to flooding (Deng et al., 2022). To determine potential submergence responsive NnNACs, the published transcriptome of completely submerged “Qiu Xing” at different time intervals (0, 3, 6, and 24 h and 5 days) was analyzed. As a result, a total of 29 NnNAC DEGs were identified (Figure 6B; Supplementary Table S6). Seventeen highly expressed NnNAC genes with FPKM ≥ 10, and exhibiting two different expression patterns were identified. Of these genes, nine (∼52.9%), including NnNAC12, NnNAC28, NnNAC32, NnNAC37, NnNAC40, NnNAC41, NnNAC42, NnNAC52, and NnNAC73 were all significantly upregulated with transcript fold increase ranging from 2 to 165.5 after 3 h complete submergence treatment. For example, a 165.5-fold increase in the expression level of NnNAC32 compared to that of control was observed after 3 h complete submergence, suggesting that it could be a crucial complete submergence stress responsive regulator in lotus. On the contrary, four NnNAC genes, NnNAC20, NnNAC35, NnNAC49, and NnNAC56 showed decreased expression profiles after complete submergence treatment, suggesting their unlikely involvement in the lotus submergence stress response. Taken together, these results provide important possible candidate NAC regulators for studying the mechanisms of seed development and response to complete submergence in lotus.
FIGURE 6. Expression analysis of NnNAC genes during lotus seed development and response to complete submergence. (A) Heat-map showing the expression patterns of 30 NnNAC genes involved in the lotus seed development (FPKM values). (B) Heat-map showing the expression patterns of 29 complete submergence stress responsive NnNAC genes in lotus (FPKM values). “SY” represent complete submergence stress treatment.
Discussion
TFs function as important switches of transcriptional networks that accurately regulate gene expression (Jiang et al., 2016; Droge-Laser et al., 2018; Kumar et al., 2021). NAC TF family has been identified as one of the largest plant-specific regulators involved in diverse biological processes of numerous plant species, such as Arabidopsis, rice, tomato, and maize (Christiansen and Gregersen, 2014; Su et al., 2015; Singh et al., 2021). However, despite their crucial roles and the availability of complete genome sequence, a comprehensive characterization of NAC gene family has not been conducted in the aquatic lotus plant.
Identification and Sequence Characteristics of NAC Gene Family in Lotus
Previous studies have reported that members of the NAC TF family are widely distributed in lower and higher plants, such as moss with 31 members, and Arabidopsis, rice, soybean, wheat, tomato, grape, and cucumber with 117, 151, 269, 263, 101, 70, and 83 members, respectively (Nuruzzaman et al., 2010; Xu et al., 2014; Ma et al., 2018; Singh et al., 2021). Our study identified for the first time 82 NnNAC genes in the published genome of sacred lotus, and this number is very close to that of cucumber. As a basal eudicot, lotus experienced a whole genome duplication (WGD) event about 18 million years ago (Ming et al., 2013; Wang et al., 2013). Consequently, genes associated with signal transduction might have been preferentially retained during WGD to form the present NAC TF family (De Smet et al., 2013). Despite wide variations in gene length, predicted protein molecular weight, isoelectric point, and exons/introns organization among the identified NnNAC genes, they remain relatively conserved with about 64% of the genes having three exons (Figure 2D; Supplementary Figure S2). Similar observations were also made in Arabidopsis, rice, tomato, tobacco, and maize (Nuruzzaman et al., 2010; Peng et al., 2015; Su et al., 2015; Li et al., 2018a; Li et al., 2021b; Singh et al., 2021). Homologous NnNAC gene pairs predominantly shared higher degree of similarities in predicted protein features, gene structures, and conserved protein motifs, which is consistent with previous findings that duplicate genes derived from a common ancestor evolve independently at the same rate with few variations (Zhang and Wang, 2005; Wang et al., 2010; Rinerson et al., 2015). Membrane-bound NAC TFs have been identified in various plant species, providing direct evidence of their involvement in many biological processes (De Clercq et al., 2013; Li et al., 2016; Bhattacharjee et al., 2017). Here, we identified five membrane-bound NAC TFs which showed over-representation in subfamilies I and II in the lotus genome (Figure 1C; Supplementary Figure S1). All the five genes belonged to the classical membrane-bound NAC TFs, harboring a single α-helical TM located at the C-terminal, which is consistent with the membrane-bound NACs identified in Arabidopsis and rice (Kim et al., 2007; Ng et al., 2013). Interestingly, a divergent form of NAC protein with two TM located in front of the NAC domain has been reported in tomato and cotton, indicating evolution and functional differentiation (Olsen et al., 2005; Bhattacharjee et al., 2017; Sun et al., 2018). Studies have shown that membrane-bound NAC TFs could be activated by specific signal on the cell membrane (Kim et al., 2007; Seo et al., 2008). In this study, a predicted NAC protein with TM, NnNAC23, was localized in the nucleus, suggesting that this gene might be regulated by post-translation and being cut off the transmembrane domain.
Functional and Evolutionary Relationship Assessment of NnNACs Using Gene Expression Patterns
Plenty of studies have reported the ubiquitous distribution of NAC TFs in the plant kingdom and their association with various biological processes (Kim et al., 2010b; Nakano et al., 2015; Yuan et al., 2020). The presence of cis-elements is crucial for the expression of NAC TF genes and their functions. Here, we comprehensively identified numerous cis-elements with various functions in the promoter region of NnNAC genes, suggesting their likely involvement in various biological processes, such as stress responses and regulation of lotus development. For example, classical stress-responsive cis-elements, including LTREs, MYB and DREB (Tran et al., 2004; Nakashima et al., 2012; Yan et al., 2017; Yuan et al., 2020), were identified in the promoter region of NnNAC genes. These results provide a basis for future studies on the transcriptional regulation mechanisms and functions of NnNAC genes.
The expression patterns of a gene family can be used to predict their functions and evolutionary relationships (Wang et al., 2013). Gene expression pattern showed multiple expression of NnNAC genes in different lotus tissues. For example, tissue-specific expression in roots was observed in five NnNAC genes, including NnNAC24, NnNAC31, NnNAC41, NnNAC48, and NnNAC50, suggesting their potential association with lotus root development. Conversely, several NnNAC genes, such as NnNAC20, NnNAC29, NnNAC37, NnNAC52, NnNAC53, and NnNAC66 were expressed in all tested tissues, thus they could be essential for various developmental processes in the aquatic lotus (Figure 5B). Overall, members of the NnNAC gene family displayed evidence of functional redundancy and evolutionary diversity in the lotus genome. Although varied expression patterns were observed among NnNACs, most homologous gene pairs shared similar expression patterns, for example, NnNAC31 and NnNAC39 were uniquely expressed in the root, NnNAC29 and NnNAC37 were constitutively expressed in all tested tissues, while NnNAC18/19 and NnNAC25/26 were not expressed in any of the tissues tested (Figure 5B). Interestingly, similar observations were made for Arabidopsis and cotton NAC genes (Jensen et al., 2010; Sun et al., 2018). Taken together, our results indicated that NAC gene family might have undergone expansion through gene duplication during evolution, resulting in functional redundancy of gene members.
Identification of Putative Candidate NnNAC Genes Involved in Lotus Seed Development and Complete Submergence Stress Response
Lotus seed has important commercial, nutritional and medicinal values, and improving its yield and quality is one of the breeding targets today (Sun et al., 2021). Despite previous studies focused on lotus seed development, little is still known about whether NnNAC genes are associated with this process (Wang et al., 2016; Li et al., 2018b). In rice, ONAC020, ONAC023, and ONAC026 showed high expression level during seed development, and their sequence variations were correlated with seed size (Mathew et al., 2016). Herein, 30 NnNAC DEGs during lotus seed development were identified for the first time, with six genes, including NnNAC12, NnNAC20, NnNAC37, NnNAC38, NnNAC48, and NnNAC52 being preferentially expressed, suggesting their potential involvement in the regulation of lotus seed development (Figure 6A).
Efforts have been made to clarify the role of NAC TFs in responses to diverse stress, such as drought, salinity and flooding (Puranik et al., 2012; Shao et al., 2015; Yuan et al., 2020; Singh et al., 2021). Lotus is an aquatic plant that experiences frequent submergence stress due to sudden frequent flooding events. Twenty-nine NnNAC DEGs associated with submergence stress response were identified in the lotus genome, of which 17 genes showed high expression levels. For example, NnNAC32, an orthologous of Arabidopsis stress response AtNAP (At1g69490) gene (Sakuraba et al., 2015; Seok et al., 2017), was significantly upregulated under complete submergence treatment, which suggest its possible regulatory role in lotus submergence stress response (Supplementary Figure S7). This result might further confirm the functional differentiation of NAC gene family. In addition, two lotus homologs, NnNAC20 and NnNAC52, sharing ∼46% amino acid sequence identity, and orthologous to the Arabidopsis stress related ANAC2 (At1g01720) NAC TF gene, were identified (Lahiri et al., 2021). However, NnNAC20 and NnNAC52 showed different expression patterns under complete submergence, suggesting they could be diverging functionally. Notably, thirteen common NnNAC genes were differentially expressed during seed development and stress response, indicating they might play multiple biological functions in lotus. This systematic characterization of lotus NAC genes in will pave way for future functional analysis to clarify the roles of these candidate genes in lotus.
Conclusion
Our study identified 82 NnNAC genes for the first time in the sacred lotus genome, and conducted a comprehensive sequence characterization, gene structure, expression profiling and transcriptional activation. In addition, NnNAC genes involved in lotus seed development and response to complete submergence stress were highlighted. In future, functional characterization of candidate NnNAC genes identified in this study is warranted. Our results provide potentially valuable NnNAC gene resources for further lotus genetic improvement.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
HnS conceived, designed, and analyzed the data. YL, MZ, and GD performed experiments. HnS and HgS wrote the manuscript. HgS and MY supervised the research. YW, JX, and YS revised the manuscript. All authors have read and approved the manuscript.
Funding
This project was supported by the National Natural Science Foundation of China (Grant No. 31872136), the Biological Resources Program CAS (Grant No. KFJ-BRP-007-009), and the Hubei Provincial Natural Science Foundation of China (Grant No. 2020CFB484).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.901838/full#supplementary-material
Supplementary Figure S1 | Structures of membrane-bound NnNAC proteins in lotus.
Supplementary Figure S2 | Statistical summary of exon numbers in NnNAC genes.
Supplementary Figure S3 | Chromosome localization of NnNAC genes in lotus.
Supplementary Figure S4 | Comparison of NAC amino acid sequence similarities.
Supplementary Figure S5 | Cis-element analysis in the promoters of NAC genes in lotus.
Supplementary Figure S6 | Validation of gene expression through qRT-PCR.
Supplementary Figure S7 | Phylogenetic relationship between NnNAC proteins and AtNAC proteins in Arabidopsis.
Supplementary Table S1 | Primer sequences used for vector construction and qRT-PCR validation.
Supplementary Table S2 | Protein sequences of NnNAC genes.
Supplementary Table S3 | Statistical analyses of sequence length, molecular weight, isoelectric point, and predicted subcellar localizations for NAC proteins in lotus.
Supplementary Table S4 | The expression of NnNAC genes in different organs and different tissues in Lotus (FPKM values).
Supplementary Table S5 | The expression of NnNAC DEGs involved in lotus seed development (FPKM values).
Supplementary Table S6 | The expression of NnNAC DEGs involved in completely submerged stress response (FPKM values).
References
Bhattacharjee, P., Das, R., Mandal, A., and Kundu, P. (2017). Functional Characterization of Tomato Membrane-Bound NAC Transcription Factors. Plant Mol. Biol. 93, 511–532. doi:10.1007/s11103-016-0579-z
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 13, 1194–1202. doi:10.1016/j.molp.2020.06.009
Chen, S., Xiang, Y., Deng, J., Liu, Y., and Li, S. (2013). Simultaneous Analysis of Anthocyanin and Non-Anthocyanin Flavonoid in Various Tissues of Different lotus (Nelumbo) Cultivars by HPLC-DAD-ESI-MSn. PLoS One 8, e62291. doi:10.1371/journal.pone.0062291
Cheng, L. B., Han, Y. Y., Liu, H. Y., Jiang, R. Z., and Li, S. Y. (2019). Transcriptomic Analysis Reveals Ethylene's Regulation Involved in Adventitious Roots Formation in lotus (Nelumbo nucifera Gaertn.). Acta Physiol. Plant 41, 97. doi:10.1007/s11738-019-2895-9
Christiansen, M. W., and Gregersen, P. L. (2014). Members of the Barley NAC Transcription Factor Gene Family Show Differential Co-Regulation with Senescence-Associated Genes during Senescence of Flag Leaves. J. Exp. Bot. 65, 4009–4022. doi:10.1093/jxb/eru046
De Clercq, I., Vermeirssen, V., Van Aken, O., Vandepoele, K., Murcha, M. W., Law, S. R., et al. (2013). The Membrane-Bound NAC Transcription Factor ANAC013 Functions in Mitochondrial Retrograde Regulation of the Oxidative Stress Response in Arabidopsis. Plant Cell 25, 3472–3490. doi:10.1105/tpc.113.117168
De Smet, R., Adams, K. L., Vandepoele, K., Van Montagu, M. C. E., Maere, S., and Van de Peer, Y. (2013). Convergent Gene Loss Following Gene and Genome Duplications Creates Single-Copy Families in Flowering Plants. Proc. Natl. Acad. Sci. U.S.A. 110, 2898–2903. doi:10.1073/pnas.1300127110
Deng, J., Li, M., Huang, L., Yang, M., and Yang, P. (2016). Genome-Wide Analysis of the R2R3 MYB Subfamily Genes in Lotus (Nelumbo Nucifera). Plant Mol. Biol. Rep. 34, 1016–1026. doi:10.1007/s11105-016-0981-3
Deng, X., Yang, D., Sun, H., Liu, J., Song, H., Xiong, Y., et al. (2022). Time-course Analysis and Transcriptomic Identification of Key Response Strategies to Complete Submergence in Nelumbo Nucifera. Hortic. Res. 9, uhac001. doi:10.1093/hr/uhac001
Droge-Laser, W., Snoek, B. L., Snel, B., and Weiste, C. (2018). The Arabidopsis bZIP Transcription Factor Family-An Update. Curr. Opin. Plant Biol. 45, 36–49. doi:10.1016/j.pbi.2018.05.001
Jensen, M. K., Kjaersgaard, T., Nielsen, M. M., Galberg, P., Petersen, K., O'Shea, C., et al. (2010). The Arabidopsis Thaliana NAC Transcription Factor Family: Structure-Function Relationships and Determinants of ANAC019 Stress Signaling. Biochem. J. 426, 183–196. doi:10.1042/BJ20091234
Jiang, J., Ma, S., Ye, N., Jiang, M., Cao, J., and Zhang, J. (2016). WRKY Transcription Factors in Plant Responses to Stresses. J. Integr. Plant Biol. 59, 86–101. doi:10.1111/jipb.12513
Kim, S.-G., Lee, S., Ryu, J., and Park, C.-M. (2010a). Probing Protein Structural Requirements for Activation of Membrane-Bound NAC Transcription Factors in Arabidopsis and Rice. Plant Sci. 178, 239–244. doi:10.1016/j.plantsci.2009.12.007
Kim, S.-G., Lee, S., Seo, P. J., Kim, S.-K., Kim, J.-K., and Park, C.-M. (2010b). Genome-Scale Screening and Molecular Characterization of Membrane-Bound Transcription Factors in Arabidopsis and Rice. Genomics 95, 56–65. doi:10.1016/j.ygeno.2009.09.003
Kim, S.-Y., Kim, S.-G., Kim, Y.-S., Seo, P. J., Bae, M., Yoon, H.-K., et al. (2007). Exploring Membrane-Associated NAC Transcription Factors in Arabidopsis: Implications for Membrane Biology in Genome Regulation. Nucleic Acids Res. 35, 203–213. doi:10.1093/nar/gkl1068
Kokkirala, V. R., Yonggang, P., Abbagani, S., Zhu, Z., and Umate, P. (2010). Subcellular Localization of Proteins of Oryza Sativa L. In the Model Tobacco and Tomato Plants. Plant Signal Behav. 5, 1336–1341. doi:10.4161/psb.5.11.13318
Kumar, R., Das, S., Mishra, M., Choudhury, D. R., Sharma, K., Kumari, A., et al. (2021). Emerging Roles of NAC Transcription Factor in Medicinal Plants: Progress and Prospects. 3 Biotech. 11, 1–14. doi:10.1007/s13205-021-02970-x
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 33, 1870–1874. doi:10.1093/molbev/msw054
Lahiri, A., Zhou, L., He, P., and Datta, A. (2021). Detecting Drought Regulators Using Stochastic Inference in Bayesian Networks. PLos One 16, e0255486. doi:10.1371/journal.pone.0255486
Le, D. T., Nishiyama, R., Watanabe, Y., Mochida, K., Yamaguchi-Shinozaki, K., Shinozaki, K., et al. (2011). Genome-Wide Survey and Expression Analysis of the Plant-Specific NAC Transcription Factor Family in Soybean during Development and Dehydration Stress. DNA Res. 18, 263–276. doi:10.1093/dnares/dsr015
Lescot, M., Dehais, P., Thijs, G., Marchal, K., Moreau, Y., Peer, Y. V., et al. (2002). PlantCARE, a Database of Plant Cis-Acting Regulatory Elements and a Portal to Tools for In Silico Analysis of Promoter Sequences. Nucleic Acids Res. 30, 325–327. doi:10.1093/nar/30.1.325
Li, H., Yang, X., Zhang, Y., Gao, Z., Liang, Y., Chen, J., et al. (2021a). Nelumbo Genome Database, an Integrative Resource for Gene Expression and Variants of Nelumbo N. Sci. Data 8, 38. doi:10.1038/s41597-021-00828-8
Li, J., Shi, T., Huang, L., He, D., Nyong’A, T. M., and Yang, P. (2018b). Systematic Transcriptomic Analysis Provides Insights into lotus (Nelumbo Nucifera) Seed Development. Plant Growth Regul. 86, 339–350. doi:10.1007/s10725-018-0433-1
Li, J., Xiong, Y., Li, Y., Ye, S., Yin, Q., Gao, S., et al. (2019). Comprehensive Analysis and Functional Studies of WRKY Transcription Factors in Nelumbo Nucifera. Int. J. Mol. Sci. 20, 5006. doi:10.3390/ijms20205006
Li, S., Wang, N., Ji, D., Xue, Z., Yu, Y., Jiang, Y., et al. (2016). Evolutionary and Functional Analysis of Membrane-Bound NAC Transcription Factor Genes in Soybean. Plant Physiol. 172, 1804–1820. doi:10.1104/pp.17.0001810.1104/pp.16.01132
Li, W., Li, X., Chao, J., Zhang, Z., Wang, W., and Guo, Y. (2018a). NAC Family Transcription Factors in Tobacco and Their Potential Role in Regulating Leaf Senescence. Front. Plant Sci. 9, 1900. doi:10.3389/fpls.2018.01900
Li, W., Zeng, Y., Yin, F., Wei, R., and Mao, X. (2021b). Genome-Wide Identification and Comprehensive Analysis of the NAC Transcription Factor Family in Sunflower during Salt and Drought Stress. Sci. Rep. 11, 19865. doi:10.1038/s41598-021-98107-4
Lindemose, S., Jensen, M. K., de Velde, J. V., O'Shea, C., Heyndrickx, K. S., Workman, C. T., et al. (2014). A DNA-Binding-Site Landscape and Regulatory Network Analysis for NAC Transcription Factors in Arabidopsis Thaliana. Nucleic Acids Res. 42, 7681–7693. doi:10.1093/nar/gku502
Love, M. I., Huber, W., and Anders, S. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Ma, X., Zhang, Y., Turečková, V., Xue, G.-P., Fernie, A. R., Mueller-Roeber, B., et al. (2018). The NAC Transcription Factor SlNAP2 Regulates Leaf Senescence and Fruit Yield in Tomato. Plant Physiol. 177, 1286–1302. doi:10.1104/pp.18.00292
Mathew, I. E., Das, S., Mahto, A., and Agarwal, P. (2016). Three Rice NAC Transcription Factors Heteromerize and Are Associated with Seed Size. Front. Plant Sci. 7, 1638. doi:10.3389/fpls.2016.01638
Mathew, I. E., Priyadarshini, R., Mahto, A., Jaiswal, P., Parida, S. K., and Agarwal, P. (2020). SUPER STARCHY1/ONAC025 Participates in Rice Grain Filling. Plant Direct 4, e00249. doi:10.1002/pld3.249
Min, T., Niu, L.-F., Xie, J., Yi, Y., Wang, L.-M., Ai, Y.-W., et al. (2019). Effects of Vacuum Packaging on NAC Gene Expression in Fresh-Cut Lotus Root. J. Amer. Soc. Hort. Sci. 145, 36–44. doi:10.21273/JASHS04806-19
Ming, R., Vanburen, R., Liu, Y., Yang, M., Han, Y., Li, L. T., et al. (2013). Genome of the Long-Living Sacred lotus (Nelumbo Nucifera Gaertn.). Genome Biol. 14, R41. doi:10.1186/gb-2013-14-5-r41
Nakano, Y., Yamaguchi, M., Endo, H., Rejab, N. A., and Ohtani, M. (2015). NAC-MYB-Based Transcriptional Regulation of Secondary Cell Wall Biosynthesis in Land Plants. Front. Plant Sci. 6, 288. doi:10.3389/fpls.2015.00288
Nakashima, K., Takasaki, H., Mizoi, J., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2012). NAC Transcription Factors in Plant Abiotic Stress Responses. Biochem. Biophys. Acta 1819, 97–103. doi:10.1016/j.bbagrm.2011.10.005
Ng, S., Ivanova, A., Duncan, O., Law, S. R., Van Aken, O., De Clercq, I., et al. (2013). A Membrane-Bound NAC Transcription Factor, ANAC017, Mediates Mitochondrial Retrograde Signaling in Arabidopsis. Plant Cell 25, 3450–3471. doi:10.1105/tpc.113.113985
Nuruzzaman, M., Manimekalai, R., Sharoni, A. M., Satoh, K., Kondoh, H., Ooka, H., et al. (2010). Genome-Wide Analysis of NAC Transcription Factor Family in Rice. Gene 465, 30–44. doi:10.1016/j.gene.2010.06.008
Olsen, A. N., Ernst, H. A., Leggio, L. L., and Skriver, K. (2005). NAC Transcription Factors: Structurally Distinct, Functionally Diverse. Trends Plant Sci. 10, 79–87. doi:10.1016/j.tplants.2004.12.010
Peng, X., Zhao, Y., Li, X., Wu, M., Chai, W., Sheng, L., et al. (2015). Genome-Wide Identification, Classification and Analysis of NAC Type Gene Family in Maize. J. Genet. 94, 377–390. doi:10.1007/s12041-015-0526-9
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T.-C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 33, 290–295. doi:10.1038/nbt.3122
Puranik, S., Sahu, P. P., Srivastava, P. S., and Prasad, M. (2012). NAC Proteins: Regulation and Role in Stress Tolerance. Trends Plant Sci. 17, 369–381. doi:10.1016/j.tplants.2012.02.004
Ren, Y., Huang, Z., Jiang, H., Wang, Z., Wu, F., Xiong, Y., et al. (2021). A Heat Stress Responsive NAC Transcription Factor Heterodimer Plays Key Roles in Rice Grain Filling. J. Exp. Bot. 72, 2947–2964. doi:10.1093/jxb/erab027
Rinerson, C. I., Rabara, R. C., Tripathi, P., Shen, Q. J., and Rushton, P. J. (2015). The Evolution of WRKY Transcription Factors. BMC Plant Biol. 15, 66. doi:10.1186/s12870-015-0456-y
Sakuraba, Y., Kim, Y.-S., Han, S.-H., Lee, B.-D., and Paek, N.-C. (2015). The Arabidopsis Transcription Factor NAC016 Promotes Drought Stress Responses by Repressing AREB1 Transcription through a Trifurcate Feed-Forward Regulatory Loop Involving NAP. Plant Cell 27, 1771–1787. doi:10.1105/tpc.15.00222
Seo, P. J., Kim, M. J., Park, J.-Y., Kim, S.-Y., Jeon, J., Lee, Y.-H., et al. (2010). Cold Activation of a Plasma Membrane-Tethered NAC Transcription Factor Induces a Pathogen Resistance Response in Arabidopsis. Plant J. 61, 661–671. doi:10.1111/j.1365-313X.2009.04091.x
Seo, P. J., Kim, S.-G., and Park, C.-M. (2008). Membrane-Bound Transcription Factors in Plants. Trends Plant Sci. 13, 550–556. doi:10.1016/j.tplants.2008.06.008
Seok, H.-Y., Woo, D.-H., Nguyen, L. V., Tran, H. T., Tarte, V. N., Mehdi, S. M. M., et al. (2017). Arabidopsis AtNAP Functions as a Negative Regulator via Repression of AREB1 in Salt Stress Response. Planta 245, 329–341. doi:10.1007/s00425-016-2609-0
Shao, H., Wang, H., and Tang, X. (2015). NAC Transcription Factors in Plant Multiple Abiotic Stress Responses: Progress and Prospects. Front. Plant Sci. 6, 902. doi:10.3389/fpls.2015.00902
Shi, T., Rahmani, R. S., Gugger, P. F., Wang, M., Li, H., Zhang, Y., et al. (2020). Distinct Expression and Methylation Patterns for Genes with Different Fates Following a Single Whole-Genome Duplication in Flowering Plant. Mol. Biol. Evol. 37, 2394–2413. doi:10.1093/molbev/msaa105
Singh, S., Koyama, H., Bhati, K. K. K., and Alok, A. (2021). Correction to: The Biotechnological Importance of the Plant-Specific NAC Transcription Factor Family in Crop Improvement. J. Plant Res. 134, 643. doi:10.1007/s10265-021-01281-9
Souer, E., van Houwelingen, A., Kloos, D., Mol, J., and Koes, R. (1996). The No Apical Meristem Gene of Petunia Is Required for Pattern Formation in Embryos and Flowers and Is Expressed at Meristem and Primordia Boundaries. Cell 85, 159–170. doi:10.1016/S0092-8674(00)81093-4
Su, H., Zhang, S., Yin, Y., Zhu, D., and Han, L. (2015). Genome-Wide Analysis of NAM-ATAF1,2-CUC2 Transcription Factor Family in Solanum Lycopersicum. J. Plant Biochem. Biotechnol. 24, 176–183. doi:10.1007/s13562-014-0255-9
Sun, H., Hu, M., Li, J., Chen, L., Li, M., Zhang, S., et al. (2018). Comprehensive Analysis of NAC Transcription Factors Uncovers Their Roles during Fiber Development and Stress Response in Cotton. BMC Plant Biol. 18, 150. doi:10.1186/s12870-018-1367-5
Sun, H., Li, J., Song, H., Yang, D., Deng, X., Liu, J., et al. (2020). Comprehensive Analysis of AGPase Genes Uncovers Their Potential Roles in Starch Biosynthesis in Lotus Seed. BMC Plant Biol. 20, 457. doi:10.1186/s12870-020-02666-z
Sun, H., Liu, Y., Ma, J., Wang, Y., Song, H., Li, J., et al. (2021). Transcriptome Analysis Provides Strategies for Postharvest lotus Seeds Preservation. Postharvest Biol. Technol. 179, 111583. doi:10.1016/j.postharvbio.2021.111583
Takada, S., Hibara, K., Ishida, T., and Tasaka, M. (2001). The CUP-SHAPED COTYLEDON1 Gene of Arabidopsis Regulates Shoot Apical Meristem Formation. Development 128, 1127–1135. doi:10.1242/dev.128.7.1127
Thirumalaikumar, V. P., Devkar, V., Mehterov, N., Ali, S., Ozgur, R., Turkan, I., et al. (2018). NAC Transcription Factor JUNGBRUNNEN1 Enhances Drought Tolerance in Tomato. Plant Biotechnol. J. 16, 354–366. doi:10.1111/pbi.12776
Tran, L.-S. P., Nakashima, K., Sakuma, Y., Simpson, S. D., Fujita, Y., Maruyama, K., et al. (2004). Isolation and Functional Analysis of Arabidopsis Stress-Inducible NAC Transcription Factors that Bind to a Drought-Responsive Cis-Element in the Early Responsive to Dehydration Stress 1 Promoter. Plant Cell 16, 2481–2498. doi:10.1105/tpc.104.022699
Wang, K., Guo, W., Yang, Z., Hu, Y., Zhang, W., Zhou, B., et al. (2010). Structure and Size Variations between 12 A and 12 D Homoeologous Chromosomes Based on High-Resolution Cytogenetic Map in Allotetraploid Cotton. Chromosoma 119, 255–266. doi:10.1007/s00412-009-0254-0
Wang, L., Fu, J., Li, M., Fragner, L., Weckwerth, W., and Yang, P. (2016). Metabolomic and Proteomic Profiles Reveal the Dynamics of Primary Metabolism during Seed Development of lotus (Nelumbo nucifera). Front. Plant Sci. 7, 750. doi:10.3389/fpls.2016.00750
Wang, Y., Fan, G., Liu, Y., Sun, F., Shi, C., Liu, X., et al. (2013). The Sacred lotus Genome Provides Insights into the Evolution of Flowering Plants. Plant J. 76, 557–567. doi:10.1111/tpj.12313
Wang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: A Toolkit for Detection and Evolutionary Analysis of Gene Synteny and Collinearity. Nucleic Acids Res. 40, e49. doi:10.1093/nar/gkr1293
Wingett, S. W., and Andrews, S. (2018). FastQ Screen: A Tool for Multi-Genome Mapping and Quality Control. F1000Research 7, 1338. doi:10.12688/f1000research.15931.2
Xu, B., Ohtani, M., Yamaguchi, M., Toyooka, K., Wakazaki, M., Sato, M., et al. (2014). Contribution of NAC Transcription Factors to Plant Adaptation to Land. Science 343, 1505–1508. doi:10.1126/science.1248417
Yamaguchi, M., Ohtani, M., Mitsuda, N., Kubo, M., Ohme-Takagi, M., Fukuda, H., et al. (2010). VND-INTERACTING2, a NAC Domain Transcription Factor, Negatively Regulates Xylem Vessel Formation in Arabidopsis. Plant Cell 22, 1249–1263. doi:10.1105/tpc.108.064048
Yan, H., Zhang, A., Ye, Y., Xu, B., Chen, J., He, X., et al. (2017). Genome-Wide Survey of Switchgrass NACs Family Provides New Insights into Motif and Structure Arrangements and Reveals Stress-Related and Tissue-Specific NACs. Sci. Rep. 7, 3056. doi:10.1038/s41598-017-03435-z
Yang, M., Zhu, L., Pan, C., Xu, L., Liu, Y., Ke, W., et al. (2015). Transcriptomic Analysis of the Regulation of Rhizome Formation in Temperate and Tropical Lotus (Nelumbo Nucifera). Sci. Rep. 5, 13059. doi:10.1038/srep13059
Yuan, C., Li, C., Lu, X., Zhao, X., Yan, C., Wang, J., et al. (2020). Comprehensive Genomic Characterization of NAC Transcription Factor Family and Their Response to Salt and Drought Stress in Peanut. BMC Plant Biol. 20, 454. doi:10.1186/s12870-020-02678-9
Zhang, Y., and Wang, L. (2005). The WRKY Transcription Factor Superfamily: Its Origin in Eukaryotes and Expansion in Plants. BMC Evol. Biol. 5, 1. doi:10.1186/1471-2148-5-1
Keywords: lotus, NAC transcription factor, transactivation, development, stress
Citation: Song H, Liu Y, Dong G, Zhang M, Wang Y, Xin J, Su Y, Sun H and Yang M (2022) Genome-Wide Characterization and Comprehensive Analysis of NAC Transcription Factor Family in Nelumbo nucifera. Front. Genet. 13:901838. doi: 10.3389/fgene.2022.901838
Received: 22 March 2022; Accepted: 17 May 2022;
Published: 08 June 2022.
Edited by:
Turgay Unver, FicusBio, TurkeyReviewed by:
Ming Li, Hubei University, ChinaPeitao Lü, Fujian Agriculture and Forestry University, China
Copyright © 2022 Song, Liu, Dong, Zhang, Wang, Xin, Su, Sun and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Heng Sun, sunheng@wbgcas.cn; Mei Yang, yangmei815815@wbgcas.cn