Identification of the NAC Transcription Factors and Their Function in ABA and Salinity Response in Nelumbo nucifera

Nelumbo nucifera Gaertn. is an important perennial aquatic herb that has high ornamental, edible, medicinal, and economic value, being widely distributed and used in China. The NAC superfamily (NAM, ATAF1/2, CUC2) plays critical roles in plant growth, development, and response to abiotic and biotic stresses. Though there have been a few reports about NAC genes in lotus, systematic analysis is still relatively lacking. The present study aimed to characterize all the NAC genes in the lotus and obtain better insights on the NnNACs in response to salt stress by depending on ABA signaling. Here, 97 NAC genes were identified by searching the whole lotus genome based on the raw HMM models of the conserved NAM domain and NAC domain. They were characterized by bioinformatics analysis and divided into 18 subgroups based on the phylogenetic tree. Cis-element analysis demonstrated that NAC genes are responsive to biotic and abiotic stresses, light, low temperature, and plant hormones. Meanwhile, NAC genes had tissue expression specificity. qRT-PCR analysis indicated that NAC genes could be upregulated or downregulated by NaCl treatment, ABA, and fluoridone. In addition, NAC016, NAC025, and NAC070, whose encoding genes were significantly induced by NaCl and ABA, were located in the nucleus. Further analysis showed the three NAC proteins had transcriptional activation capabilities. The co-expression network analysis reflected that NAC proteins may form complexes with other proteins to play a role together. Our study provides a theoretical basis for further research to be conducted on the regulatory mechanisms of salinity resistance in the lotus.


Introduction
Plant transcription factors (TFs) play crucial roles in plant growth, development, and abiotic and biotic responses [1]. Considering the importance of TFs in responding to various stimuli, analyzing their function mechanism is becoming indispensable for vegetable crop breeding. With more plant genomes sequenced and transgenic technology developed, a large number of TFs have been identified, and their functions have been determined. NAC (non-apical meristem (NAM), Arabidopsis transcription activation factor (ATAF), and calathiform cotyledon (CUC)) transcription factors characterized by the presence of a highly conserved N-terminal NAC domain and diverse C-terminal domains have been widely studied [2,3]. The N-terminal NAC domain consists of about 160 amino acids that could bind to the cis-elements of their target genes, and C-terminal domains are a fundamental area to interact with other proteins and regulate other gene expression [4,5]. The structure of NAC proteins being well-identified provided a framework for further understanding the functions of NAC proteins at the transcriptional level. Subsequently, research on the functions of NAC transcription factors mushroomed in a variety of species.
Gene family analysis as a basic method plays an important role in investigating the functions of plant transcription factors. The members of the NAC family have been identified and reported in many species by gene family analysis. For instance, 108 members have been identified in Arabidopsis, 151 members in rice, 150 members in soybean, 261 members in M. sinensis, 132 members in peanut, 151 members in sunflower, 87 members in maize, 111 members in celery, 104 members in pepper, and 110 members in potato [6][7][8][9][10][11][12]. NAC proteins have not only been identified, but their functions have also been well studied. NAC TFs play important roles in many vital biological processes and plant stress responses [13][14][15][16][17][18].
Considering that abiotic stresses could adversely affect crop growth and productivity, it is very important to study the functions of NAC TFs in improving crop resistance to abiotic stresses, which have been well elucidated in many species in last decade. In grain crop rice, ONAC066, ONAC095, OsNAC2, OsNAC006, OsNAC67, OsNAC016, and OsNAP function as positive or negative regulators of the drought stress response [15,[19][20][21][22][23][24]. In vegetable crops, NAC TF genes, such as JUB1, SlVOZ1, SlNAC35, StNAC053, and CaNAC46, also play an important role in plant responses to drought stresses [25][26][27][28][29]. In addition, NAC TFs have been widely focused on regarding plant resistance to salinity, which reflects the fact that one NAC TF can have multiple functions. Potato StNAC053, grapevine VvNAC17, pepper CaNAC46, rice ONAC022, and wheat TaNAC29 enhance not only plant drought tolerance but also salt tolerance [25,26,[30][31][32]. These research works provide us with a clear research direction and a very important basis from which to elucidate the regulation mechanism of NAC TFs.
NAC TFs can regulate several plant signal networks and be regulated by environmental stimuli, plant hormones, and miRNAs. Environmental stimuli, such as light, temperature, drought, salinity, water, etc., could induce the expression of NAC TFs [33][34][35]. Furthermore, plant hormones are the essential regulators of plant growth and development. The interplay between plant hormones and NAC TFs has become an interesting research topic. Wu et al. discussed interaction between NAC TFs and ethylene and abscisic acid (ABA) pathways in controlling tomato fruit development and ripening [36]. OsNAC3, OsNAC45, and VvNAC17 improved transgenic plant ABA sensitivity and salt tolerance [31,37,38]. Further study demonstrated that CrNAC036 can directly bind and negatively regulate CrNCED5 expression to hinder ABA biosynthesis in citrus [39]. Interestingly, recent studies showed microRNAs (miRNAs) comprising about 22 nt endogenous noncoding RNAs could take part in plant abiotic responses by targeting NAC genes. In other words, NAC genes were regulated by plant miRNAs. Peu-miR164 regulated its target genes PeNAC070 and PeNAC081, which are involved in the poplar response to abiotic stress [40]. osa-miR164c negatively regulated OsCUC1 and OsCUC3, which are involved in rice meristem/organ boundary specification [41]. In addition, miR1514a modulates NAC700 transcription to trigger phasiRNA formation in response to water deficit in legumes [42]. These reports indicate that the NAC TF-related networks are complex and important for plant growth and development.
Lotus (Nelumbo nucifera Gaertn.), popular ornamental, edible, and medicinal plant which belongs to the small family Nelumbonaceae, will be an important model plant in horticulture [43]. With three lotus genomes sequenced, the lotus TFs, such as AP2/ERF, WRKY, bHLH, and the MADS gene family, have been widely studied [44][45][46][47][48][49][50]. However, the NAC gene family has not yet been comprehensively studied in lotus. Though 82 NnNAC genes have been found using the PlantTFDB 4.0, the NAC gene family still needs to be identified by searching the whole lotus genome. Here, all NAC members of lotus were identified, and their distribution, phylogenetic relationship, gene structure, and synteny were analyzed. In addition, certain candidate NAC TFs which had a higher expression in lotus leaves were further studied though subcellular localization and expression level analysis by qRT-PCR. These results will further deepen the understanding of NAC TFs and provide several candidate genes related to salinity and abscisic acid (ABA) response in lotus.

Identification of NAC Transcription Factors in Nelumbo nucifera
To identify all NAC TFs in Nelumbo nucifera, the conserved NAM domain (PF02365) and NAC domain (PF01849) were used to seek the respective raw HMM models. Based on these two HMM models, 107 and 106 candidate members were identified, and the redundant forms of the same genes were removed. In total, 97 NAC TF members were identified and named NAC001-NAC097 according to their chromosome location (Table 1 and Figure 1). The basic information, including genomic length, amino acid residues, pI, and molecular weight, was calculated. For example, the polypeptide lengths of the predicted NACs varied widely, ranging from 76 to 684. In addition, the pI data showed that NACs have both acidic and alkaline members. The information also contained the members' gene IDs, which could help researchers obtain recent news on NACs (Table 1).

Phylogenetic Analysis and Classification of NAC Genes
Eighteen NAC subgroups have been classified in rice and Arabidopsis [51]. In our results, 18 new subgroups were redefined in which no lotus NAC members belong to ONAC001 and OsNAC3. However, the lotus NAC members family had two specific subgroups named NNACO23 and NNAC014, which indicated lotus is species-specific compared to Arabidopsis and rice ( Figure 2).

Gene Structure and Motif Composition of the NAC Gene Family
The sequence features including introns, exons, and motifs were visualized by the TBtools software [52]. According to the frequency of occurrence, motifs 1-5 were the five most frequently presented motifs, covering almost all NAC members. Furthermore, motifs 8-10 only existed in the members NnNAC082, 083, 084, 085, and 031, which belong to subgroup NNAC023. Interestingly, NnNAC009 and NnNAC017 possessed two copies of motif 3. In addition, the motifs in subgroup ONAC003 showed a character of random arrangement. The exhibition of exon-intron structure showed that there is evolutionary diversity. Most of the members had more than three exons, and only four members had one exon. NnNAC069, 27, 28, and 79 possessed relatively longer intron lengths (Figure 3).

Cis-Acting Elements and Synteny Analysis of NAC Gene Promotors
To investigate the functions of NAC genes, the cis-acting elements of NAC gene promotors were analyzed by the online tool PlantCARE; the data were then visualized by TBtools. All NAC genes had the "light responsiveness" element, suggesting that the "light responsiveness" element is fundamental for plant growth and development. In addition, plant hormones were the important factors for plant growth, development, and response to the environment. Furthermore, the plant hormone elements "abscisic acid responsiveness", "gibberellin responsiveness", "MeJA responsiveness", "salicylic acid responsiveness", and "auxin responsiveness" were widely presented in NAC gene promotors (Figure 4). The "defense and stress responsiveness" element, which weas also widely presented in NAC gene promotors, played important roles in plant responses to abiotic and biotic stresses. The "low-temperature responsiveness" element was also important because lotus root is occasionally affected by low temperatures. The above elements suggest that NAC genes in lotus root may function in diverse aspects of lotus root growth and development. To compare the evolutionary events between lotus and Arabidopsis, synteny analysis was carried out using the orthologous NAC gene pairs. A total of 69 collinear genes were retrieved from the two species ( Figure S1).

Expression Patterns of NAC Genes in Different Plant Tissues
To investigate the tissue-specific expression of NAC genes, we obtained data from the "Nelumbo Genome Database" which contained the expression levels of NAC genes in the petiole, leaf, root, and cotyledon (Table S2). The results showed that NAC genes had tissue expression specificity, indicating that their functions were diverse. For example, NAC006 was only expressed in the leaf, root, and petiole. Unlike NAC006, three genes, NAC020, 012, and 071, were only expressed in the cotyledon. The leaf was the most important organ, and its phenotypes were easily observed after stress treatment. We selected the 10 most highly expressed genes-NAC003, 016, 025, 034, 043, 045, 060, 065, 070, and 078-for further research.

qRT-PCR Analysis of the NAC Genes
To investigate the functions of the NAC genes, the above selected genes were analyzed by qRT-PCR under NaCl, ABA, and ABA + FL treatment ( Figure 5). Under NaCl treatment, the transcription levels of nine NAC genes increased gradually, with NnNAC078 being the sole exception. Under ABA treatment, the 10 NAC genes were all induced, and the transcription level of NnNAC003 fluctuated. Fluoridone, the inhibitor of ABA, also induced gene transcription. In most cases, FL affects the action of ABA. In addition, the interaction of ABA and FL had positive effects on the expression levels of NAC genes including NnNAC003, NnNAC025, NnNAC045, and NnNAC065. The expression levels of NnNAC16, NnNAC25, and NnNAC70 were significantly increased under ABA and NaCl treatment, while their expression levels were significantly reduced under ABA and FL treatment; they were selected for further research.

Subcellular Localization of NAC Proteins
Subcellular localization could indicate the role and position of NAC-encoded proteins. In order to understand the molecular characteristics of hormone-and stress-responsive lotus root NACs, we selected three lotus root NAC genes for subcellular localization. Green fluorescence of the NnNAC16, NnNAC25, and NnNAC70 proteins all appeared in the nucleus, while green fluorescence of the control protein appeared on the cell membrane, in the cytoplasm, and in the nucleus, which indicated that the NnNAC16, NnNAC25, and NnNAC70 proteins function in the nucleus ( Figure 6). Figure 6. Subcellular localization of the soybean NAC proteins. The first row of pictures above is a control, and those below it are the NnNAC016, NnNAC025, and NnNAC070 proteins. Scale bars = 25 µm.

The Transcriptional Activation Capabilities Analysis of NAC Proteins
Transcription factors generally have self-activation ability. The full-length cDNAs of NnNAC16, NnNAC25, and NnNAC70 were cloned into a PGBKT7 vector, and the selfactivation ability of NnNAC16, NnNAC25, and NnNAC70 was verified by yeast experimentation. The results showed that NnNAC16, NnNAC25, NnNAC70, and the negative control pGBKT7 could all grow colonies on an SD/-Trp medium, while only the negative control pGBKT7 failed to grow colonies on an SD/-Trp/-Ade/His medium. This indicated that NnNAC16, NnNAC25, and NnNAC70 all had self-activation ability (Figure 7). To further investigate the functions of the three NnNAC genes, their homologs in Arabidopsis were identified. Then, their co-expression networks were drawn using the website STRING. The results indicated that the functions of NAC genes were diverse and closely connected with NaCl, drought, and the plant hormone ABA ( Figure S2). Coexpression networks allowed us to better study NAC genes in lotus root.

Discussion
With the completion of the genome sequencing of lotus root, several gene families were identified [44,47,53]. NAC proteins have diverse functions and play vital roles in plant growth and development, including responses to biotic and abiotic stresses [14,17,[54][55][56][57][58][59][60]. However, functional research of lotus NAC genes has rarely been reported. Eighty-two NAC genes were obtained using the website PlantTFDB, and we studied their potential functions in lotus seeds and response to complete submergence [61]. Compared to NAC genes of other plants, the number NAC genes in lotus root still needs to be found. Here, we identified 97 NAC gene members by searching the whole lotus genome based on the raw HMM models for the conserved NAM domain (PF02365) and NAC domain (PF01849), which could look for as many genes as possible. Then, the 97 members were divided into 18 subgroups based on the classification method used in other plants [62,63]. The promotor, as an important structure, plays vital roles in plant growth and development. All NAC genes had the "light responsiveness" element, which is the fundamental element for plant growth and development. ABA is an important phytohormone regulating plant growth, development, and stress responses [64]. The "abscisic acid responsiveness" element in the promotor of NAC genes indicated that NAC genes may be regulated by the ABA-related transcription factors. Meanwhile, the qRT-PCR analysis also showed that NAC genes were induced under ABA treatment (Figure 4). In addition, the co-expression network of Arabidopsis homologs of three NnNAC genes also indicated the functions of NAC genes and ABA signaling are closely connected (Figure 7). For example, PwNAC11 activated ERD1 by interaction with ABA-related proteins, ABF3 and DREB2A, to enhance drought tolerance in transgenic Arabidopsis [65]. OsNAC45 was involved in ABA response and salt tolerance in rice [38]. In addition, "gibberellin responsiveness", "MeJA responsiveness", "salicylic acid responsiveness", and "auxin responsiveness" were also widely presented in NAC gene promotors, which indicates NAC genes are multifunctional, interacting with many plant hormones' signaling. These could provide us with a new avenue by which to study the function of NAC genes. Previous research showed that MeJA enhanced ethylene synthesis by inducing NAC genes in kiwifruit [66]. In addition, NAC genes could function in response to biotic stresses based on the "defense and stress responsiveness" element in their promotors. In total, the functions of NAC genes still need to be further studied.
The lotus is an important ornamental, edible, and medicinal plant originating from India and China. It is famous for its products such as fresh rhizomes, lotus root starch, drinks, teas, and lotus seeds [67]. However, the plant's yield and product quality are affected by the salinity of soil [68]. Therefore, to excavate the key resistance genes is an important work for breeding superior lotus varieties. Furthermore, our results showed that NAC genes could be candidate genes to improve lotus resistance to slat stress. For instance, NnNAC003, NnNAC016, NnNAC043, NnNAC045, NnNAC060, and NnNAC070 were significantly induced under NaCl treatment (Figure 4). Further results indicated that their functions in response to salt stress were dependent on ABA signaling (Figures 3 and 4). In order to understand where NAC proteins function, we conducted a subcellular localization assay which demonstrated that NAC proteins function in the nucleus, as do most transcription factors ( Figure 5). Like most transcription factors, the three NAC proteins, NnNAC16, NnNAC25, and NnNAC70, have transcriptional activation, which could bind them to the related cis-elements to regulate the expression levels of target genes. In addition, NAC proteins may form complexes with other proteins to play a role together (Figure 7).
In total, the 97 NAC genes may be involved in responses to biotic and abiotic stresses by depending on ABA or other plant hormone signaling and would be the key candidate genes for further functional research by modern genetic and biological techniques.

The Identification of NAC Transcription Factor Members
The Nelumbo nucifera Geartn. genome (China Antique v2.0) was obtained from the Nelumbo Genome Database [48]. The raw HMM models for the conserved NAM domain (PF02365) and NAC domain (PF01849) were downloaded from the Pfam database (https: //pfam.xfam.org/ (accessed on 27 May 2022)). the raw HMM models were then used to search the Nelumbo nucifera Geartn. genome with the software TBtools [52] to obtain the candidate NAC transcription factor members. Then, certain NAC transcription factor members were determined by their sequence search in the Pfam database with the threshold of an e-value <e−5. Next, to exclude repeated sequences, all NAC sequences were aligned using DANMAN Version 7 and checked manually.

Basic Information and Phylogenetic Analyses
The NAC genes were visualized on the chromosomes using the software TBtools [52]. Their CDS length, molecular weight, and isoelectric point were predicted by the ExPASy tool (https://www.expasy.org/ (accessed on 5 June 2022)). The amino acid sequences of Arabidopsis and Nelumbo nucifera Geartn. were download from TAIR (https://www. arabidopsis.org/ (accessed on 5 June 2022)) and the Nelumbo Genome Database [48], respectively. Based on their amino acid sequences, the phylogenetic tree was constructed in MEGA7.0 by the neighbor-joining (NJ) algorithm; 1,000 replicates were used to evaluate the significance of nodes [69]. The annotation on the phylogenetic tree was created in the Adobe Illustrator SC5 software.

Synteny and Gene Structure Analyses
The applet MCScanX (default parameters) in TBtools was used to examine duplicate genes. Then, the homology of the NAC gene between Nelumbo nucifera Geartn. and Arabidopsis was visualized using the applet Dual Synteny Plotter in TBtools [52]. The conserved motifs of NAC proteins were identified using the online tool MEME Suite. The cis-element of the NAC gene promotor was analyzed via the website PlantCARE [70]. The motifs and cis-element were exhibited in one figure by TBtools.

Microarray Analysis of NAC Gene Expression Patterns
The tissue expression data of NAC genes was downloaded from the Nelumbo Genome Database. The analysis was carried out, which included the 97 NAC genes in the different tissues and development stages, and then it was visualized by TBtools.

Subcellular Localization Analysis
The selected proteins were fused with the expression vector pCAMBIA1300, containing a GFP tag for subcellular localization analysis, as described previously [71]. A diluted agrobacterium solution (NAC-pCAMBIA1300 or pCAMBIA1300) was injected into the leaves of tobacco. After 24 h of incubation, fluorescence images were captured by a confocal laser scanning microscope (TCS SP8; Leica). The primers used in this study can be found in Table S1.

Plant Materials and Treatments
Lotus seeds were germinated in soil in a greenhouse. Two-week-old seedlings were subjected to exogenous ABA and salinity treatments. For these treatments, seedlings were transferred to solutions containing 200 mM NaCl or 100 µM ABA. Seedlings were sampled at 0, 1, 2, 4, 8, 12, and 24 h after treatment. The samples were then dropped immediately into liquid nitrogen and stored at −80 • C.

Quantitative Real-Time PCR (qRT-PCR)
The samples were ground in liquid nitrogen, and the powder was used to extract the RNA. The cDNA synthesis was conducted following the manufacturer's instructions using a PrimeScript™ 1st Strand cDNA Synthesis Kit (TaKaRa, Beijing China). Then, the expression patterns were analyzed with a BIO-RAD CFX MaestroTM system. In this assay, three biological and technical replicates were made. The histogram was made with the GraphPad Prism 8.0 software, and the standard error was calculated with the SPSS software. The primers used in this study can be found in Table S1.

The Transcriptional Activation Assays
The transactivation ability of NAC TFs was tested in yeast cells. The fusion plasmids NnNACs:pGBKT7 and the control plasmid pGBKT7 were transformed into Y2H Gold yeast cells. Transformed yeast competent cells were spread on plates containing SD/-Trp or SD/-Trp/-His/-Ade.

Conflicts of Interest:
The authors declare no conflict of interest.