The NAC (NAM, ATAF1/2, and CUC2) family of proteins is one of the largest plant-specific transcription factor (TF) families and its members play varied roles in plant growth, development, and stress responses.
Crops are constantly challenged by a variety of abiotic and biotic factors that have negative impacts on their growth, development and yields [1]. Biotic stressors, mainly bacteria, fungi and viruses, can cause diseases and impact agricultural crops, forestry plantations, and native plant communities [2]. To resist or tolerate these adverse external stresses, plants have evolved many effective stress response strategies, including changes in morphology, establishing defensive systems, and constructing physiological, biochemical and molecular regulatory networks.
Plants have two main branches of defense against various pathogens: PAMP (pathogen-associated molecular patterns)-triggered immunity (PTI) and effector-triggered immunity (ETI) [3]. PTI is a basal defense mechanism in plants that is triggered by the recognition of PAMPs through PRRs (pattern recognition receptors) at the plant cell surface [4]. ETI is an accelerated and amplified PTI response, which arises from the interactions between plant resistance proteins and pathogen effector proteins [5,6][5][6]. Upon being attacked by pathogens, plants will activate a series of complex molecular regulatory networks, such as reactive oxygen species (ROS) signaling [7], phytohormone signaling [8], changes in redox status [9], and inorganic ion fluxes [10], to prevent further pathogen invasion.
Much of the plant response to pathogens involves transcriptional reprogramming. Plants have established high-efficiency gene expression networks to regulate multiple specific stress responsive genes in a coordinated manner. Such regulation of the large-scale expression of genes requires a concerted function of different types of transcription factors (TFs). TFs specifically bind to cis-elements and/or trans-acting factors in the promoters of target genes and act as transcriptional activators or repressors. Several TF families, such as NAC (NAM, ATAF and CUC), MYB (myeloblastosis-related proteins), WRKY (WRKYGQK), bZIP (basic leucine zipper domain), bHLH (basic helixloop-helix), CAMTA (CaM-binding transcription activator) and ERF/AP2 (ethylene responsive factor/apetala2), have crucial roles in abiotic and biotic stress responses [11]. So far, a large number of studies have proved that NAC TFs play roles in plant growth, development, abiotic stress response, and disease resistance [11,12,13,14,15,16,17][11][12][13][14][15][16][17].
Pathogens directly affect plant growth and development and decrease the quality and yields of the crop. Three major cereal crops, namely rice (Oryza sativa), maize (Zea mays), and wheat (Triticum aestivum), provide two-thirds of the food consumed all over the world [18]. Barley (Hordeum vulgare) and soybean (Glycine max) are also important crops, in the grain and legume families, while tomato (Solanum lycopersicum) is a fruit crop. NAC TFs comprise a gene family with 151 members in rice [19[19][20][21],20,21], 157 in maize [22], 559 in wheat [23], 167 in barley [24], 152 in soybean [25], and 93 in tomato [26]. Recent studies demonstrate that NAC TFs can be induced by pathogen infection, regulate the expression of downstream genes, and integrate hormone and other signals to offer plants resistance against pathogens.
NAC TFs constitute one of the largest groups of plant-specific TFs, with more than 100 members in most plants, including model plants like Arabidopsis or most crops, and the highest reported number of 559 in wheat [23]. The NAC acronym is derived from three reported proteins that contain a highly conserved domain in their N-terminal region (the NAC domain): (i) NAM (No Apical Meristem), (ii) ATAF1/2 (Arabidopsis thaliana Transcription Activator Factor 1/2) and (iii) CUC2 (Cup-shaped Cotyledon 2) [20,21][20][21]. NAC TFs play essential roles in diverse biological processes, such as growth, development, senescence and morphogenesis, and are widely involved in signaling pathways in response to different phytohormones and multiple abiotic and biotic stress [27,28,29][27][28][29].
NAC proteins are widespread in land plants, from “simple” bryophytes to “complex” angiosperms, and also in freshwater green algae (Charophytes) [30[30][31][32],31,32], a sister clade to land plants and likely their evolutionary root [33]. NAC proteins originated about 725–1200 million years ago (Mya), during the diversification of the charophytes/embryophytes (the Streptophytes), experienced a first expansion in bryophytes about 470 Mya, and a second expansion in angiosperms during the early Cretaceous period [27]. Along with the divergence of vascular plants, NAC proteins extensively expanded [30] [30] shortly before the origination and radiation of angiosperms through different (segmental or tandem) duplications [31,34,35][31][34][35]. In angiosperms, there are two detectable evolutionary pathways of the NAC proteins: ancient duplications that occurred before the divergence of dicots and monocots and recent duplications in a particular lineage of dicots or monocots [30,36][30][36]. As a result of duplication, some family members possess redundant functions [37,38][37][38].
According to the clustering analyses by Zhu et al., NAC proteins consist of 21 subfamilies among bryophytes and vascular plants [30]. Among these subfamilies, 15 were found only in angiosperms, while the other six occur in both flowering plants and earlier diverged lycophytes [30,33][30][33]. In 2015, another team recategorized NAC proteins into six major orthologous groups (Group I-VI) using more than 2000 non-redundant sequences from 24 different species of green plants by comparative genomic and gene functional analyses [31]. NAC group I is considered the basal NAC group, because its members are involved in secondary wall and wood formation, which might have been essential to water conduction or support during the adaptation of plants to land environments [39,40][39][40]. Group II NACs function in specific developmental processes among different organs, as well as other processes such as ethylene-auxin pathways. NAC group III is named the TMM (transmembrane motifs) Group, as 142 of the 164 NAC proteins in this group contained a TMM in their C-terminal regions that can anchor them in biomembranes (especially on the endoplasmic reticulum), and can be hydrolyzed by proteolytic enzymes activated by signaling [41,42,43][41][42][43]. The members in NAC group IV have diverse functions, such as ANAC009, which controls the reorientation and timing of cell division, and ANAC042 which regulates longevity of Arabidopsis. NAC group V is proposed to be the Stress Group, as most of its members are involved in stress responses. Finally, NAC group VI contains many species-specific sequences, which experienced whole-genome duplication events through their evolutionary history [31].
NAC proteins usually contain two relatively independent domains: one well-conserved N-terminal NAC domain of 151–159 amino acids and a relatively divergent C-terminal Transcriptional Activation Region (TAR) [20,44,45][20][44][45]. There are a few kinds of atypical NAC proteins, which have variations in either their NAC domains or TAR [46]. The NAC domain exists in all NAC proteins, and the most conserved consensus sequences are D-D/E-L-I/V, E-W-Y-F-F, G-Y-W-K, and M-H-E-Y [46]. The NAC domain contains nuclear localization signals (NLSs) and its main function is translocating the protein from the cytoplasmic matrix to the nucleus [44,47][44][47] and forming homo- or heterodimers, that the state in which these proteins bind DNA [48,49,50][48][49][50]. Some proteins also contain nuclear export signals (NESs) in the NAC domain, which might mediate their export out of the nucleus for degradation after their missions are completed [32]. Crystallographic structures of ANAC019 from Arabidopsis and SNAC1 (STRESS RESPONSIVE NAC1) from rice revealed that the NAC domain contains mainly twisted antiparallel β-sheet(s) flanked by a few α-helices. The central β-sheet is responsible for dimerization of the proteins, by forming stabilized salt bridges between conserved amino acids in two subunits, and also for interacting with the major groove of DNA in the dimer form [49,50,51][49][50][51]. The NAC domain can be further divided into five subdomains, designated A-E, distinguished by blocks of heterogeneous amino acids or gaps, among which subdomains A, C, and D are highly conserved compared to B and E. It is assumed that subdomain A promotes functional dimerization, subdomains C and D play parts in DNA binding, due to their content of basic amino acids, and in nuclear translocation of the protein with the help of the NLSs inside the subdomains, while the variable subdomains B and E diversify the roles of the NAC proteins [44,48][44][48].
The TAR region, also known as TRR (transcriptional regulatory region), is related to transcriptional regulation and corresponds to the distinct functions of NAC proteins [46[46][48],48], acting as both activator and/or repressor to downstream target genes together with the NAC domain [27]. Bioinformatic analysis revealed common motifs in the C-terminal regions of some NAC subfamilies [48,52][48][52]. In addition, an α-helical TMM could be found at the C-terminal end of many NAC proteins [43] [43] that contain a conserved N-terminal NAC domain, a variable middle TAR region, and a C-terminal TMM motif [42]. With the help of the TMM, the nascent NAC protein can be anchored to the endoplasmic reticulum membrane or plasma membrane, retaining it outside of the nucleus in a dormant state. Subsequently, the NAC TF can be cleaved through either protease- or ubiquitin proteasome-mediated proteolytic events upon stimulation by internal and/or environmental signals, which would activate the mature nuclear-associated form to enter the nucleus and regulate downstream genes, which is an adaptive strategy to environment [41,42,43][41][42][43]. Recently, transmembrane domains were found in the N-terminal ends of a small number of NAC proteins [32].
In both model and crop plants, numerous NAC genes are induced in response to pathogen infection [46]. Over-expression or silencing of certain NAC genes results in enhanced or reduced resistance to pathogens [53,54,55,56][53][54][55][56], suggesting that NAC TFs could positively or negatively regulate plant defense response. Here, we introduce the roles of NAC TFs during the response to pathogens in several important crops, including rice, wheat, barley, potato, soybean, maize, tomato and lettuce.
Rice is one of the most important food crops in the world, as it is consumed by more than 3 billion people [57]. Rice blast caused by the fungal pathogens Magnaporthe oryzae and M. grisea is one of the most devastating diseases threatening rice production. Several NACs have been found to respond to rice blast, therefore, it is important to understand the role of these TFs in response to both this disease and other pathogens in rice.
The rice genome is predicted to contain 151 NAC genes [19], but only a few have been characterized. OsNAC6 was the first NAC TF found to be involved in disease resistance in rice [58]. OsNAC6 can be induced by both abiotic and biotic stress, including blast disease. Transgenic rice plants over-expressing OsNAC6 had improved tolerance to rice blast, dehydration and high salt stress, but exhibited negative effects on growth and yields [58]. Similarly, resistance was obtained from the over-expression of OsNAC111 or OsNAC58. Over-expression of OsNAC111 or OsNAC58 in rice increased resistance to M. oryzae or the bacterial blight pathogen Xoo (Xanthomonas oryzae pv. oryzae), respectively [59,60][59][60]. In addition, inoculation with M. oryzae can induce OsNAC111 expression, which then activates several downstream defense genes, such as PR1 (PATHOGENESIS-RELATED1) and PR8, suggesting that the OsNAC111 TF positively regulates the expression of a specific set of PR genes in the rice blast response [59]. ONAC122 and ONAC131 also play positive regulatory roles during rice blast resistance and can be induced by infection with M. grisea and treatment with exogenous defense signaling molecules, such as salicylic acid (SA), methyl jasmonate (MeJA), or 1-aminocyclopropane-1-carboxylic acid (a precursor of ethylene) [61]. Metabolomics analysis was recently used, for the first time, to explore the molecular mechanism of NAC-mediated disease resistance. Liu et al. found that ONAC066 promotes the resistance against fungal blast and bacterial blight in rice by regulating the accumulation of soluble sugars and amino acids as well as the up-regulation of the PR gene [62].
The hypersensitive response (HR) of plants against disease is a type of programmed cell death. Kaneda et al. found that over-expression of OsNAC4 in rice could lead to HR cell death, accompanied with the loss of plasma membrane integrity, fragmentation of nuclear DNA, and typical morphological changes. HR cell death is noticeably decreased in OsNAC4 knock-down lines after induction by an avirulent pathogen N1141 [63], suggesting that OsNAC4 is a key positive regulator of plant hypersensitive cell death. In addition, the osnac60 mutant had increased susceptibility to M. oryzae, while the over-expressing plant had an enhanced defense response, including increased programmed cell death, ROS accumulation, and callous deposition and up-regulation of defense-related genes [64].
On the contrary, Yoshii et al. found that mutation of RICE DWARF VIRUS MULTIPLCATION 1 (RIM1), which encodes a NAC TF, resulted in a loss of susceptibility to Rice Dwarf Virus (RDV) but not to Rice Transitory Yellowing Virus (RTYV) and Rice Stripe Virus (RSV) [65]. The accumulation of virus capsid protein was significantly reduced in rim1 mutant plants inoculated with RDV, which impairs the multiplication of the virus, while proliferation of the virus was stimulated with over-expression of RIM1. Therefore, it was proposed that RIM1 negatively regulates rice resistance to RDV, by acting as a host factor that is required for multiplication of the virus [65].
This summary shows that 9 NAC TFs (OsNAC4, OsNAC6, OsNAC58, OsNAC60, ONAC066, OsNAC111, ONAC122, ONAC131, and RIM1) have been validated to take part in defense responses against pathogen attack. OsNAC6, OsNAC58, OsNAC60, ONAC066, OsNAC111, ONAC122 and ONAC131 positively regulate rice resistance to fungus Magnaporthe, while RIM1 negatively regulates rice resistance to the virus RDV.
Wheat is the most widely distributed food crop and occupies the largest planting area in the world [66]. The planting area and total output account for one third of all food crops, and about one-third of the population of the world uses wheat as its main food [67]. Thus, it is quite urgent to develop novel wheat varieties that have improved yield potential and increased tolerances to biotic and abiotic stresses. There are 559 TaNAC genes in the wheat genome, 186 of which were recently identified during infection of the resistant hexaploid wheat line N9134 with the fungal agents that cause rust and powdery mildew [23].
Powdery mildew, caused by the biotrophic Blumeria graminis f. sp. tritici (Bgt), is a widespread disease in wheat [68]. Zhou et al. isolated and characterized three wheat TaNAC6s, TaNAC6-A, TaNAC6-B and TaNAC6-D. Analysis of the expression patterns of the TaNAC6s showed that TaNAC6-A and TaNAC6-D were up-regulated early after Bgt inoculation. Over-expression of the TaNAC6s enhanced the resistance to Bgt, while silencing them reduced the resistance in wheat, indicating that the TaNAC6s play positive roles in resistance against Bgt [69]. Recently, wheat TaNACL-D1 was demonstrated to be involved in the disease response to Fusarium head blight (FHB) [70], a devastating disease of wheat and barley in humid and semi-humid regions of the world. TaNACL-D1 was recently shown to interact with TRITICUM AESTIVUM FUSARIUM RESISTANCE ORPHAN GENE (TaFROG), which enhances wheat resistance against FHB, and lines over-expressing TaNACL-D1 were more resistant to FHB disease [70,71][70][71].
Previous studies have shown that wheat TaNAC4 and TaNAC8 shared high homology with rice OsNAC4 and OsNAC8. TaNAC4 and TaNAC8 expression is induced by infection with Puccinia striiformis f. sp. tritici (Pst, which causes stripe rust) as well as by treatments with MeJA or ethylene (ET), but not by treatments with abscisic acid (ABA) or SA [72,73][72][73]. This indicates that TaNAC4 and TaNAC8 of wheat might be two important components in defense-signaling pathway and play essential roles in resistance to pathogen.
TaNAC1 is strongly expressed in wheat roots and is involved in responses to Pst infection and treatments with defense-related hormones. Knockdown of TaNAC1 enhances resistance to Puccinia stripe rust in wheat, while over-expression of TaNAC1 in Arabidopsis enhances susceptibility and reduces systemic-acquired resistance to Pseudomonas syringae pv tomato DC3000 (Pst DC3000). Similarly, knockdown of TaNAC21/22 in wheat enhanced resistance against Puccinia stripe rust [74,75][74][75]. Silencing of TaNAC2 or TaNAC30 enhances the resistance against stripe rust by significantly increasing H2O2 generation and decreasing hyphal growth at the early stage of the interaction between Pst and wheat [76,77][76][77]. These data indicate that TaNAC1, TaNAC2, TaNAC21/22, and TaNAC30 negatively regulate stripe rust resistance in wheat.
Barley is the fourth most important grain crop in the world [78] [78] and, like the other cereal crops, is threatened by various pests and diseases [79]. A search of the public barley sequence database identified 48 NAC genes (HvNACs), while the expression profiles of 46 HvNACs were investigated in various tissues with and without ABA or MeJA treatment. The HvNAC proteins have conserved functions in secondary cell wall biosynthesis, leaf senescence, root development, seed development and hormone-regulated stress response [80].
Barley HvNAC6 has a high similarity to the rice OsNAC6 in pathogen resistance. Transient over-expression of the gene in barley increased the penetration resistance of epidermal cells to the powdery mildew pathogen Blumeria graminis f. sp. hordei (Bgh), while silencing of the gene, using RNA interference (RNAi), enhanced the sensitivity to Bgh. Silencing of HvNAC6 also changed the accumulation of ABA, which was not affected by Bgh inoculation, indicating that HvNAC6 acts as an ABA-mediated defense response regulator to maintain basal resistance against Bgh [81,82][81][82]. In Arabidopsis, the expression of the HvNAC6 homologue ATAF1 was induced by Bgh, and the ataf1-1 mutant line displayed reduced penetration resistance to Bgh. HvNAC6 and ATAF1 play conserved positive roles in penetration resistance in both monocots and dicots, respectively [82].
Another NAC TF in barley is HvSNAC1, which promotes resistance to barley ramularia leaf spot (RLS) disease [83]. RLS is a new emerged barley disease, caused by the ascomycete fungus Ramularia collo-cygni and which broke out in Europe a decade ago [84]. Over-expression of HvSNAC1 in barley significantly reduced the severity of RLS, but had no effects on other pathogenic diseases, such as eyespot, powdery mildew, or blast. Further analysis showed that dark-induced leaf senescence is delayed in HvSNAC1 over-expression lines, indicating that HvSNAC1 may inhibit plant senescence [83].
Tomato is the highest value vegetable and fruit crop worldwide, at an annual production of 100 million tons, and makes a huge nutritional contribution to the human diet [85,86,87][85][86][87]. At the same time, tomato is constantly attacked by various pathogens, causing huge losses in production [88,89,90][88][89][90]. So far, 93 putative NAC proteins were identified in tomato [26]. Similar to other crops, tomato NAC TFs play roles in abiotic and biotic stress responses, as well as the development of the plant [91,92,93,94,95][91][92][93][94][95].
Infection by pathogens induces the expression of SlNAC1 in tomato, but plays dual functions in resistance to different pathogens. SlNAC1 expression is specifically induced in tomato by the replication enhancer (REn) of Tomato leaf curl virus (TLCV) [96], and its over-expression increases the accumulation of viral DNA in infected cells, indicating that SlNAC1 play negative roles in resistance against TLCV. SlNAC1 is also induced by Pseudomonas infection, but plays reversed roles in defense signaling [56]. In Pseudomonas-infected plants, expression of SlNAC1 increased rapidly while degradation of the SlNAC1 protein was suppressed. Further research proved that SlNAC1 could be ubiquitinated by SINA3, a ubiquitin ligase, but the expression of SINA3 was decreased in infected plants. Thus, pathogen infection counteracts the degradation of the SlNAC1 protein. These data suggest that SlNAC1 plays a positive role in resistance to Pseudomonas infection [56,97][56][97].
The NAC protein, Solanum lycopersicum Stress-related NAC1 (SlSRN1), was identified in tomato by virus-induced gene silencing technology [98]. The expression of SlSRN1 can be significantly induced by infection with Botrytis cinerea and Pst DC3000, while silencing of SlSRN1 leads to increased severity of the diseases. Silencing of SlSRN1 accelerates accumulation of ROS but reduces expression of defense genes after infection by B. cinerea. These results demonstrate that SlSRN1 is a positive regulator of the defense response against B. cinerea and Pst DC3000 in tomato [98]. Recently, six NAC TFs (SlNAC24, SlNAC20, SlNAC39, SlNAC47, SlNAC61 and SlNAC69) were studied in response to Tomato yellow leaf curl virus (TYLCV) infection in tomato. Four NAC genes (SlNAC20, SlNAC24, SlNAC47, and SlNAC61) were induced after TYLCV infection in resistant plants, and SlNAC61 played positive roles in response to TYLCV infection, according to Virus-induced gene silencing analysis. Furthermore, the six NAC TFs could interact with protein phosphatase 2C (PP2C), mitogen-activated protein kinase 3 (MPK3), and some defense response TFs, such as WRKY, MYB, and even NAC, by binding the promoters of these genes, indicating that NAC TFs have an complex response mechanism during TYLCV infection [94]. Recently, it was found that SlNAC082, a ribosomal stress mediator, was involved in the process of infection by citrus exocortis viroid (CEVd) in tomato. A higher expression level of SlNAC082 was detected in the CEVd-infected tomato leaves. CEVd and its derived viroid small RNAs were found to co-sediment with tomato ribosomes in vivo and caused alterations in ribosome biogenesis in the infected tomato plants. The alterations in both the rRNA processing and the induction of SlNAC082 were correlated with the degree of viroid symptomology [95].
Stomata play an active part in the plant innate immune response, and serve as an entrance for pathogen into plant cells [99]. The genes JA2 (Jasmonic Acid 2) and JA2L (JA2-like) both encode two NAC TFs that are closely related to ANAC019/ANAC055/ANAC072. These NACs were preferentially expressed in guard cells of tomato leaves [100]. In JA2-SRDX (SUPERMAN REPRESSION DOMAIN X) plants, Pst DC3000-induced stomatal closure was impaired at 1 h post infection (hpi), and pathogen-triggered stomatal reopening remained normal at 4 hpi, indicating that JA2 is required for Pst DC3000-induced stomatal closure but not stomatal reopening. By contrast, the Pst DC3000-induced stomatal closure was largely normal at 1 hpi, but the pathogen-triggered stomatal reopening was substantially impaired at 4 hpi in JA2L-AS plants expressing an antisense version of the JA2L cDNA, indicating that JA2L is required for pathogen-regulated stomatal reopening [100].
Potato is one of the four major food crops around the world. A NAC TF, StNACb4 from potato, was identified and characterised. StNACb4 has been shown to promote resistance to bacterial wilt caused by Ralstonia solanacearum [101]. Transgenic tobacco plants were generated in which the expression of StNACb4 was constitutively up-regulated or suppressed using RNAi. StNACb4 was found specifically in the phloem of the vascular system of the stems and leaves, and up-regulated upon infection with R. solanacearum or by treatment with SA, ABA and MeJA in transgenic tobacco. Silencing StNACb4 reduced the tolerance of tobacco to R. solanacearum, and over-expression of the gene enhanced the tolerance to this pathogen [101]. These results are consistent with findings on StNAC43, another potato NAC TF, which can increase the deposition of resistance-related metabolites to reinforce the secondary cell wall and improve resistance to late blight disease [102]. These data demonstrate that both StNACb4 and StNAC43 are positive regulators of disease resistance of potato.
Maize is not only an important and widely distributed cereal crop, but also a model plant for genetic research [103]. Plant diseases induced by pathogens cause huge yield losses, up to 41.1% every year [104]. Lu et al. identified 157 non-redundant maize NAC genes, which were unevenly distributed on 10 maize chromosomes [22]. Further sequence and evolutionary relationship analysis showed that 19 maize NAC genes were related to stress responses [105]. ZmNAC41 and ZmNAC100 were transcriptionally induced during infection by Colletotrichum graminicola and defense signals, and were also expressed during leaf senescence in maize. In addition, ZmNAC41 was up-regulated in response to the fungal biotroph Ustilago maydis. Interestingly, the transcripts of ZmNAC41 and ZmNAC100 are induced by JA and SA, respectively, suggesting that ZmNAC41 and ZmNAC100 could function in the defense response [106]. When the upstream promoters of maize NAC genes were analyzed, a MYC binding site was detected in ZmNAC15, ZmNAC38 and ZmNAC41, while a WRKY-binding motif was detected in ZmNAC15, ZmNAC36, ZmNAC41, and ZmNAC100. In short, the ZmNAC15, ZmNAC36, ZmNAC38, ZmNAC41, and ZmNAC100 genes all contained potential binding elements for TFs known to be involved in the plant defense network [106].
Soybean is a main source of high-quality proteins and a vegetable oil that provide nutrition for animals and humans [107]. NAC TFs are believed to play vital roles in soybean development and disease resistance. Six NAC-like genes, designated GmNAC1–GmNAC6, were cloned and characterized from soybean a decade ago. These genes had similar genomic organization and high sequence similarity, especially in the NAC domains, but exhibited different expression patterns during seed development [108]. Subsequently, more NAC proteins were identified in soybean [109], most of which are involved in development and abiotic stress, such as GmNAC30 [110] [110], GmNAC81 [110[110][111],111], GmNAC109 [112] [112] and GmNAC8 [113], while up to now only GmNAC42 is reported to be involved in plant disease resistance [114,115,116][114][115][116]. Soybean GmNAC42-1 is a homolog of the Arabidopsis ANAC042-1, which is an indole alkaloid plant antitoxin regulator. Over-expression of GmNAC42-1 in elicited hairy roots significantly increases the amount of glyceollin in soybean, suggesting this protein is an essential and positive regulator of glyceollin biosynthesis. GmNAC42 is annotated as a systemic acquired resistance (SAR) gene and functions in soybean disease resistance because glyceollins are defensive metabolites (phytoalexins) derived from isoflavones in soybean [114]. MYB TFs are also involved in the glyceollin gene regulatory network. GmMYB29A1 and GmMYB29A2 were up-regulated in hairy roots treated with a wall glucan elicitor from P. sojae, and the expression of GmNAC42-1 and GmMYB29A1 were increased with the over-expression of GmMYB29A2, indicating that GmNAC42-1 was also regulated by GmMYB29A2 during glyceollin biosynthesis [116].
In lettuce, LsNAC069, a NAC TF with a C-terminal TMM motif, is a target of the RxLR-like effectors of the fungus Bremia lactucae [117]. RxLR effectors are characterized by a conserved RxLR (Arg-x-Leu-Arg) motif in the N-terminal domain, and B. lactucae secretes potential RxLR effectors during the infection process. LsNAC069 silencing increases resistance to Pseudomonas cichorii bacteria. LsNAC069 is relocalized from the ER to the nucleus when wild-type plants are treated with Phytophthora capsici culture filtrate, but this process could be prevented by the protease inhibitor TPCK (N-tosyl-L-phenylalanine chloromethyl ketone), indicating that the LsNAC069 needs proteolytic cleavage to be untethered from the ER and relocalized to the nucleus. However, the susceptibility to B. lactucae was not significantly altered in LsNAC069 silenced lettuce lines, and the process of LsNAC069 relocalization was inhibited upon the expression of B. lactucae effectors. Moreover, both co-localization and yeast two-hybrid experiments demonstrated that LsNAC069 could interact with B. lactucae effectors. Together these data demonstrate that B. lactucae can cause disease in lettuce through its RxLR effectors inhibiting the hydrolysis and relocalization of LsNAC069 from the ER to the nucleus, which suppresses the activation of genes downstream of LsNAC069 [117].
As mentioned above, we constructed a phylogenetic tree of NAC TFs cited in the review except LsNAC069 in lettuce, which cannot be searched out, and constructed a diagram of those genes (Figure 1). Meanwhile, we also made a diagram of the NAC domain structure of most NAC protein sequences mentioned here (Figure 2).
Figure 1. Phylogenetic and gene structure analysis of NAC TFs mentioned in this review. The left is phylogenetic analysis of some NAC TFs from rice, wheat, barley, maize, soybean, tomato and potato. The right is gene maps of the same NAC TFs, showing the numbers of exon and intron, the absence or presence of upstream/downstream sequences and their relative length.
Figure 2. Diagram of the NAC proteins structure and sequence alignment of NAC domains. (a) Schematic representation showing a typical NAC protein with a highly conserved NAC domain at the N-terminal which is further divided into five conserved subdomains (A–E). The NAC domain contains nuclear localization signals (NLSs) and DNA-binding domain (DBD). The C-terminal region is a relatively divergent transcriptional activation region (TAR). In some cases, the C-terminal may have a transmembrane motif (TMM). (b) Sequence alignment of NAC domains for most of crops cited in the review. Subdomains A to E are shown by horizontal lines above the sequence. The red color represents amino acids with more than 75% common consensus sequences, and the cyan color represents amino acids with more than 50% common consensus sequences.