- Open Access
Over-expression of a NAC 67 transcription factor from finger millet (Eleusine coracana L.) confers tolerance against salinity and drought stress in rice
BMC Biotechnologyvolume 16, Article number: 35 (2016)
NAC proteins (NAM (No apical meristem), ATAF (Arabidopsis transcription activation factor) and CUC (cup-shaped cotyledon)) are plant-specific transcription factors reported to be involved in regulating growth, development and stress responses. Salinity responsive transcriptome profiling in a set of contrasting finger millet genotypes through RNA-sequencing resulted in the identification of a NAC homolog (EcNAC 67) exhibiting differential salinity responsive expression pattern.
Full length cDNA of EcNAC67 was isolated, characterized and validated for its role in abiotic stress tolerance through agrobacterium mediated genetic transformation in a rice cultivar ASD16.
Bioinformatics analysis of putative NAC transcription factor (TF) isolated from a salinity tolerant finger millet showed its genetic relatedness to NAC67 family TFs in related cereals. Putative transgenic lines of rice over-expressing EcNAC67 were generated through Agrobacterium mediated transformation and presence/integration of transgene was confirmed through PCR and southern hybridization analysis. Transgenic rice plants harboring EcNAC67 showed enhanced tolerance against drought and salinity under greenhouse conditions. Transgenic rice plants were found to possess higher root and shoot biomass during stress and showed better revival ability upon relief from salinity stress. Upon drought stress, transgenic lines were found to maintain higher relative water content and lesser reduction in grain yield when compared to non-transgenic ASD16 plants. Drought induced spikelet sterility was found to be much lower in the transgenic lines than the non-transgenic ASD16.
Results revealed the significant role of EcNAC67 in modulating responses against dehydration stress in rice. No detectable abnormalities in the phenotypic traits were observed in the transgenic plants under normal growth conditions. Results indicate that EcNAC67 can be used as a novel source for engineering tolerance against drought and salinity stress in rice and other crop plants.
One of the major challenges in future agriculture is to sustain food grain production under changing climate and limited natural resources such as water and nutrients. Abiotic stresses viz., drought, salinity, flooding and temperature extremes are becoming major threats to increased agricultural productivity under fragile environments. Among these, drought and salinity remain at the top in affecting productivity of major food grains such as rice, wheat, maize etc., and predicted climate change may cause serious salinization of more than 50 % of arable lands by 2050 . This necessitates genetic manipulation of tolerance against these two major abiotic stresses in major food crops which remains very difficult through conventional breeding methods due to complexity of tolerance mechanisms [2, 3]. In this context, advancements in the fields of molecular breeding and genetic engineering offer us a powerful tool for genetic manipulation of these traits . During the recent past, genetic engineering has been successful in developing transgenic crop plants engineered for their tolerance against biotic/abiotic stresses, enhanced nutritional quality and various agronomic traits [4–6]. Genetic engineering strategies for improving salinity tolerance in crop plants includes manipulation of various metabolic pathways viz., accumulation of osmolytes, antioxidant enzymes, regulating the uptake/compartmentalization of salts, transcriptional factors and various signalling pathway components etc., . In this context, identification and validation of novel genes associated with various component traits controlling salinity tolerance in candidate crops/organisms is an important step which will allow us to manipulate salinity tolerance in any agriculturally important crop.
Recently, numerous studies have shown that transcription factors (TFs) play an important role in regulating responses against various stresses in plants and some of them have been shown to be essential for stress tolerance . Genetic manipulation of transcription factors have been demonstrated to be an effective strategy in manipulating complex traits like drought/salinity tolerance rather than modification of individual genes involved in key metabolic pathways due to the ability of TFs in modulating expression of hundreds of downstream genes involved in various metabolic pathways associated with stress tolerance in plants [8–11]. Several studies have clearly demonstrated that numerous TFs, such as DREB, bZIP, zinc-finger, MYB, WRKY and NAC, directly or indirectly regulate plant responses under abiotic stress conditions [11–14]. The NAC protein (NAM, No apical meristem; ATAF, Arabidopsis transcription activation factor and CUC, Cup-shaped cotyledon) super-family is one of the largest plant-specific TF families containing a highly conserved N-terminal DNA-binding domain, a nuclear localization signal sequence and a variable C-terminal domain [15–18]. Around 138, 158, 149 and 289 NAC family members have been reported from Arabidopsis, Rice, Setaria and Populus trichocarpa, respectively . NAC TFs play an important role in growth, development including pattern formation of embryos and flowers , secondary wall formation [21, 22], leaf senescence  and root development  in plants. Besides being involved in plant growth and development; NAC TFs were reported to be involved in modulation of responses against various biotic and abiotic stresses [10, 24–26] showing their potential to improve biotic and abiotic stress tolerance through genetic engineering .
Several genomics and bioinformatics studies have led to identification of number of drought/salinity responsive NAC TFs in Arabidopsis [28, 29], rice  and soy bean . Results of above studies suggest that stress-responsive NAC TFs may have important roles in providing tolerance against abiotic stresses and their over-expression can improve stress tolerance in crop plants. Transgenic rice plants engineered with a NAC TF (OsNAC6) were found to exhibit enhanced tolerance against various abiotic (drought and salinity) and biotic stresses . Similarly, transgenic cotton plants engineered with a stress responsive NAC TF namely SNAC1 in rice were found to exhibit enhanced tolerance against drought and salinity stresses. The transgenic cotton plants were found to possess enhanced root development and reduced transpiration rate . Similar observations were reported when rice plants were engineered with OsNAC10  and SNAC1 . In another study, TaNAC67 from wheat was found to improve tolerance against drought, salinity and freezing stresses in Arabidopsis .
It has been hypothesized that exploitation of highly saline tolerant “halophytes” or wild germplasm may serve as an excellent strategy for understanding physiological/molecular mechanisms underlying salinity tolerance, and thereby, leading to identification of novel candidate gene(s) for engineering salinity tolerance in agriculturally important crop plants [36–38]. Finger millet (Eleusine coracana L.) is one of the resilient cereal crops belonging to the family, Poaceae and genetically close to rice  which is known for its high degree of tolerance against drought, salinity and blast disease [40, 41]. Our earlier studies on transcriptome profiling of salinity responsiveness in a set of contrasting finger millet genotypes resulted in the identification of several novel putative candidate genes for functional validation . In the present study, efforts were made to isolate and validate the function of a novel NAC transcription factor namely EcNAC67 exhibiting contrasting salinity responsive expression pattern between the susceptible and tolerant finger millet genotypes. Full length gene encoding NAC transcription factor i.e. EcNAC67 was isolated from a salinity tolerant finger millet genotype, Trichy 1, and transgenic plants of a rice variety ASD16 over-expressing EcNAC67 were developed and evaluated for their responses against drought and salinity stresses.
Genetic material and stress treatments
Based on the results of our earlier study , a putative candidate gene namely a NAC transcription factor 67 (EcNAC67) exhibiting contrasting salinity response between susceptible (CO 12) and tolerant (Trichy 1) finger millet genotypes was selected for functional validation. Saline tolerant finger millet genotype “Trichy 1” was grown under normal greenhouse conditions up to 21 days (when plants were 4–5 leaf stage) and salinity stress was imposed by adding 300 mM of NaCl by maintaining suitable control plants irrigated with normal water. Leaf (top 3 leaves), root and shoot tissues were collected from both control and salinity stressed plants (20 days after stress) of Trichy 1 and used for expression analysis.
Expression analysis of EcNAC67
Tissue samples collected from control and salinity stressed plants of Trichy 1 were frozen in liquid nitrogen and used for total RNA extraction as per manufacturer’s protocol (Biobasic Inc., Canada). Equal amount of DNAse treated total RNA (about 1 μg) was converted to sscDNA using Transcriptor High Fidelity cDNA Synthesis Kit (Roche, Germany) and used for qRT-PCR analysis (StepOne Plus, Applied Biosystems, USA) by following default cycling conditions (10 min 95 °C, 40 cycles of 15 s at 95 °C and 60 s at 60 °C). The reaction mixture contained SYBRGreen Master mix (Roche Diagnostics) 300 nM of gene specific primers (NAC RT-F 5’-TCAGCAGCAGATGATGGTG-3’ and NAC RT-R 5’-CGGATCAGGTTCAGGTTCTTCG-3’) (see Additional file 1) and 2 μl of cDNA in each 15 μl reaction. “No template controls” (NTC) containing all of the RT-PCR reagents except the cDNA template were also maintained to rule out cross contamination. Abundance (relative quantity) of mRNAs was calculated using the comparative Ct (ΔΔCt method; ). qRT-PCR analysis was repeated using samples collected from three biological replications including two technical replications per biological replication and Actin was used as an endogenous reference gene for the normalization of Ct values (see Additional file 1).
Isolation and characterization of cDNA encoding EcNAC67
Gene specific primers were designed for isolating full length cDNA encoding for the candidate NAC transcription factor (EcNAC67) based on the alignment of finger millet transcript reads against its rice homologue at both 5’ and 3’ UTRs . Primers were designed with the flanking restriction sites viz., BamHI in the forward primer (5’-cgc ggatcc CAG GAG GGA GAG AGG AAA GAG-3’) and KpnI site in the reverse primer (5’- cgcggtacc C GGA TCA GGT TCA GGT TCT TCG-3’). Full length cDNA encoding EcNAC67 was PCR amplified (150 ng of each primer, 200 mM dNTPs, 2.5 U XT5 DNA polymerase in a 50 mL reaction, with 94 °C, 5 min for 1 cycle, 94 °C, 1 min, 60 °C, 1 min, and 72 °C, 1.5 min for 30 cycles; and 72 °C, 7 min for 1 cycle) from the sscDNA synthesized from salinity tolerant finger millet genotype Trichy1. Amplified PCR products were purified and cloned in pTZ57RT vector as per manufacturer’s protocol using InsTAclone PCR Cloning Kit (ThermoScientific, USA) and sequenced (SciGenom Labs, India).
In silico analysis
Nucleotide and translated amino acid sequence analysis was performed using BLASTn/BLASTp search against RNA/cDNA/protein sequences in the NCBI database (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Multiple sequence alignment of deduced amino acid sequences of EcNAC67 against other known NAC sequences from related crop species was carried out using CLUSTALW tool in BioEdit software and used for phylogenetic analysis (MEGA6 software using the maximum likelihood method with 1,000 bootstrap replications). Various other properties viz., secondary structure (PSIPHRED; http://bioinf.cs.ucl.ac.uk/psipred/); pI/Mw (Compute pI/Mw; http://web.expasy.org/compute_pi/); functional region (PROSITE; http://www.expasy.org/) and subcellular localization (ProtCompv9.0; http://www.softberry.com/berry.phtml?topic = protcomppl) were also analyzed.
Functional validation of EcNAC67
Construction of plant transformation vector harboring EcNAC67
Full length cDNA encoding EcNAC67 was released from pTZ57RT through BamHI/KpnI restriction digestion and ligated in a plant transformation vector pCAMBIA1300 under the control of RD29 promoter and nos terminator. Putative recombinant clones were selected and used for confirmation of the presence and orientation of the transgene through PCR analysis (using M13F and M13R primers) and sequencing. Then, pCAMBIA1300 harboring the transgene was mobilized into the Agrobacterium strain LBA4404 through freeze thaw method .
Genetic transformation of rice variety ASD16 using pCAMBIA1300 harboring EcNAC67
Immature embryos of rice variety ASD16 were co-cultivated with Agrobacterium strain LBA4404 harboring pCAMBIA1300 + EcNAC67 as suggested by Hiei and Komari . Transformed calli were subjected to selection in a medium containing 50 mg/l hygromycin and putative transgenic calli were regenerated in the presence of 1 mg/ml of NAA and 3 mg/ml 6-BA and rooted on half MS media containing 50 mg/lhygromycin . Putative transgenic plants were transferred to transgenic green house for hardening and establishment.
Molecular characterization of putative transgenic plants (T0)
Putative transgenic plants of ASD16 engineered with EcNAC67 were subjected to PCR analysis for confirming the presence of selectable marker gene (hygromycin), transgene (EcNAC67) and vir gene using gene specific primers (Additional file 1) and southern hybridization analysis. Genomic DNA isolated from the putative transgenic plants was digested with BamHI, electrophoresed on 1 % agarose gel and blotted onto a positively charged nylon membrane along with suitable non-transgenic ASD16 and hybridized using P32 labeled fragments of hygromycin (hpt) marker gene . After hybridization, membranes were washed, dried and exposed to X-ray film (Kodak Photo Film) overnight and developed.
Evaluation of transgenic rice plants engineered with EcNAC67 against salinity and drought
Transgenic ASD16 rice plants (T0) confirmed through PCR and southern hybridization analysis were selfed and forwarded to T1 generation. About 50 plants (T1) from each event were raised and subjected to PCR analysis to identify transgene positive and negative plants. Positive transgenic plants in all the events were allowed to grow till maturity in soil filled pots and observations on morphological characters were recorded. Seeds were collected from each plant individually and used for evaluation of tolerance against salinity and drought. Transgenic lines (T2) were evaluated for their level of salinity tolerance at both germination and vegetative stage along with non-transgenic (NT) controls. At germination stage, seeds of both transgenic and non-transgenic ASD16 were germinated in petri plates containing different concentrations of NaCl (0 mM, 75 mM, 100 mM and 150 mM). After every 24 h, fresh NaCl solutions of respective concentration were used for replacing existing solutions. Germination percentage, root length (cm) and shoot length (cm) were recorded on 10th day. For vegetative stage screening, seeds of both transgenic and non-transgenic ASD16 were germinated in petri plates (upto 7 days) and then transferred to hydroponics system in trays filled with Yoshida solution (grown up to 30 days). Presence of transgene was again confirmed through PCR analysis and salinity stress was imposed by adding 100 mM NaCl to the Yoshida solution. Effect of salinity stress on both transgenic and non-transgenic plants was assessed based on the development of wilting and drying of leaves. After 35 days of 100 mM NaCl stress, plants were transferred back to normal condition (in pots filled with soil) and allowed to grow till maturity to assess the recovery ability.
Transgenic plants were also evaluated for their tolerance against drought stress. Seeds of different transgenic lines were germinated in por trays and PCR analysis was carried out to identify transgene positive plants. Positive plants were transplanted in pots containing 2.5 kg of soil and kept in trays filled with water along with non-transgenic controls. All the plants were maintained under greenhouse conditions (70 % humidity, 28 °C, 12 h light/12 h dark cycle). Drought stress was imposed by withholding water at 54 days after sowing (DAS) and pots were saturated with water for overnight before imposing the stress. Observations on soil moisture content(%) and leaf relative water content (%) were recorded at 8 days after stress when non-transgenic ASD16 started showing severe stress symptoms. RWC was calculated by using the formula,
Where, FW is fresh weight of leaf, DW is dry weight of leaf and TW is turgid weight of leaf.
Abiotic stress responsiveness of EcNAC67 in rice
Leaf samples were collected from control and drought/salinity stressed plants of all transgenic lines along with non-transgenic ASD16 plants and used for analyzing the expression of transgene through qRT-PCR as described in previous section (Expression analysis of EcNAC67). Salinity responsive expression pattern of all the transgenic lines was compared against non-transgenic ASD16 lines. Ubiquitin was used as an endogenous reference gene for the normalization of Ct values (see Additional file 1).
Isolation of a salinity responsive NAC transcription factor from finger millet
In our earlier study , RNA-seq analysis was carried out in a set of contrasting finger millet genotypes to monitor the salinity responsive changes at transcript level. A putative candidate NAC domain-containing protein homologous to a rice NAC67 protein (LOC_Os03g60080) was found to be significantly up-regulated in the tolerant finger millet genotype Trichy 1 when compared to the susceptible genotype CO 12. Based on the sequence information of finger millet contigs mapped at 5’ and 3’ UTRs of OsNAC67, gene specific primers were designed to isolate full length cDNA of EcNAC67 from Trichy 1 genotype. Cloning and sequencing data revealed that EcNAC67 was found to be 1178 bp in size including a 36 bp 5’UTR, 969 bp open reading frame and 173 bp 3’UTR.
Validation of differential regulation of EcNAC67 during salinity stress
qRT-PCR analysis of salinity responsive expression pattern in control and salinity stressed leaves, shoots and roots of Trichy 1 revealed that, EcNAC67 transcript was found to be abundant in leaves, followed by roots and minimum in shoots (Fig. 1). Transcript abundance was found to be highly up-regulated in all the three tissues during salinity (300 mM) and the level of up-regulation was found to be the highest in shoots followed by leaves and roots (Fig. 1).
Cloning, sequencing and in silico characterization of EcNAC67
Full length cDNA encoding putative NAC transcription factor was isolated from the salinity tolerant finger millet genotype Trichy 1 using gene specific primers designed based on the alignment of finger millet transcript contigs against rice homologue (LOC_Os03g60080) through RT-PCR. Cloning and sequence analysis revealed that it belongs to NAC 67 family (based on the structural domain) of transcription factors. EcNAC67 shared significant homology with various NAC 67 family members viz., 87 % with S. italica (XP_004981295), 79 % with O. brachyantha (XP_006650786.1), 77 % with Z. mays (XP_008644394.1) and 75 % with B. distachyon (XP_003563852.1). It was predicted that it encodes for a protein of 322 amino acids with a molecular mass of 35.98 kDa and pI of 6.34. Conserved domain analysis revealed that EcNAC67 contains two functional domains namely, N-terminal NAC domain (24–176 AA) which was found to be highly conserved (Fig. 2a) and a C-terminal domain which was exhibiting greater degree of variation across different members. Prediction of putative structure of EcNAC67 protein revealed that EcNAC67 contains 9 coils, 3 alpha helices and 5 beta-strands (see Additional file 2). Prediction of sub-cellular localization using ProtComp v9.0 online tool suggested that EcNAC67 was a typical nucleus localization protein with a certainty score of 0.9 to1. Phylogenetic analysis of EcNAC67 with other NAC family members of different species showed that EcNAC67 formed a separate lineage from already reported E. coracana NAC as well as other NAC transcription factors reported from other species (Fig. 2b). Evolutionary divergence between sequences using Poisson correction model showed that EcNAC67 was having least divergence with S.italica NAC67 and S.bicolor SNAC1 (see Additional file 3).
Generation and characterization of transgenic rice plants over-expressing EcNAC67
Ninety two putative T0 transgenic plants were screened for transgene integration by PCR using primers specific to selectable marker gene hygromycin (hpt), transgene (EcNAC67) and also vir gene of Agrobacterium. All transgenic plants were found to be positive for hpt and transgene EcNAC67 but were negative for vir gene of Agrobacterium showing that PCR amplification is not due to contaminating Agrobacterium. Two PCR positive transgenic plants (T0) were selected from each transgenic event and used for confirmation of transgene integration by Southern blot analysis. Results of southern blot analysis indicated that three transgenic lines namely, EcNAC67-E1, EcNAC67-E4 and EcNAC67-E6 were found to contain single-copy of the transgene; two transgenic lines namely, EcNAC67-E2 and EcNAC67-E3 were found to harbor two copies of the transgene and one transgenic line EcNAC67-E5 was found to possess 3 copies of the transgene (Fig. 3). One of the plants (Plant # A) from a putative transgenic line EcNAC67-E3 was found to be positive for the presence of transgene in the PCR analysis, but did not show hybridization signals in Southern blot analysis (Fig. 3). Southern confirmed and PCR positive transgenic lines were selfed and harvested individually.
Fifty seeds from each T0 lines were germinated in soil filled portrays and allowed for germination. PCR analysis using gene specific primers of both transgene and hpt revealed that 35–38 plants out of 50 plants were positive in single copy insertion lines (i.e. EcNAC67-E1 EcNAC67-E4 and EcNAC67-E6); 42 – 47 plants out of 50 were positive in two copy insertion lines (i.e., EcNAC67-E2 and EcNAC67-E3) showing a normal 3:1 Mendelian segregation ratio. PCR positive T1 plants were selfed and T2 transgenic lines were evaluated for their tolerance against abiotic stresses viz., drought and salinity.
Over-expression of EcNAC67 confers tolerance against salinity in rice
Non-transgenic and transgenic ASD16 lines were allowed to germinate in petri plates containing 75, 100 and 150 mM NaCl solutions and effect of salinity on the development of shoot and roots was measured on 10th day. Salinity stress had significant effect on the growth of both shoot and root in non-transgenic ASD16 than the transgenic lines. All the transgenic rice lines were found to possess relatively less reduction in their shoot and root length when compared to non-transgenic ASD16 at 75 mM, 100 mM and 150 mM NaCl stress (Fig. 5). Non-transgenic ASD16 plants were found to have very small root at 150 mM NaCl stress but transgenic lines were having 3–4 times longer roots suggesting that the over-expression of EcNAC67 in rice confers enhanced level of tolerance against salinity stress.
At 100 mM NaCl stress during vegetative stage, non-transgenic ASD16 plants exhibited growth retardation and wilting of terminal leaves at 11 days after stress at which all the transgenic lines were found to be healthy (Fig. 6a). At 33 days after stress, non-transgenic ASD16 plants were found to be severely affected and most of the leaves were found to be dried (Fig. 6b). All the transgenic lines retained greenness in leaves and had higher root/shoot biomass (Fig. 6d). Non-transgenic ASD 16 plants showed 60 % reduction in the total biomass during salinity stress where as in the transgenic lines it ranged between 31 – 44 % (Fig. 6d). After 33 days of 100 mM NaCl stress, all the plants were allowed for recovery by transferring to pots filled with soil and irrigated using normal water. All the transgenic lines were able to recover completely within 15 days after revival and were able to reach maturity and set seeds, whereas non-transgenic ASD16 plants were not able to revive from salinity injury and died (Fig. 6c and e).
Over-expression of EcNAC67 improves drought tolerance in rice
To understand the effect of EcNAC67 over-expression on drought tolerance, T2 transgenic plants were grown in pots and severe drought stress was imposed. Both the non-transgenic and transgenic ASD 16 lines were subjected to approximately equal intensity of stress by allowing the soil moisture in pots to reach around 15 – 16 %. Upon drought, transgenic plants showed much delayed leaf-rolling symptom when compared to non-transgenic ASD16 plants (Fig. 7). All the transgenic plants were able to maintain 20 % (approx.) higher relative water content in the leaves (Fig. 8).
Upon re-watering, transgenic lines were able to revive and resumed growth rapidly when compared to non-transgenic plants (Fig. 7). Non-transgenic plants showed 47 % reduction in final plant height whereas transgenic lines showed height reduction between 15-21 % (Table 1). Further, reduction in grain yield/plant was 69.6 % in non-transgenic ASD 16 whereas reduction among transgenic lines ranged from 23.3 % to 29 % with exception to EcNAC67-E2 transgenic line where yield reduction was around 50 % (Table 1). Further, drought induced spikelet sterility was found to be higher in the non-transgenic ASD16 (40.82 %) than in the transgenic lines where it ranged between 11.70 and 21.84 % (Table 1).
Drought/salinity responsive expression pattern of transgene
qRT-PCR analysis of transgene expression in leaves of control, salinity and drought stressed plants of both non-transgenic and transgenic ASD16 showed that expression of transgene was absent in the non-transgenic ASD16 rice plants whereas all the transgenic lines were found to possess significantly increased level of transgene expression under stress conditions (Fig. 4). Among the transgenic lines, single copy lines (EcNAC67-E1 and EcNAC67-E4) were found to have maximum induction of transgene expression under both drought and salinity treatments when compared to the lines having multiple integrations of transgene (Fig. 4).
Drought and salinity are becoming major abiotic stresses limiting agricultural productivity worldwide. Predicted climate change is expected to increase the frequency of occurrence of these stresses and posing serious threat to global food security. Developing drought and salinity tolerant crop varieties will help in sustaining increased productivity under agricultural areas that are prone to such stresses. Conventional breeding efforts are resulting in a slow progress in achieving this goal due to complexity of mechanisms controlling tolerance against these stresses and lack of reliable high throughput phenotyping. In this context, use of biotechnological tools viz., marker assisted breeding and genetic engineering offers us a powerful tool for genetic manipulation of these traits. Under genetic engineering, one of the promising strategies is to modulate the expression levels of stress responsive transcription factors that might regulate wide array of downstream genes/pathways and thus bringing the desired levels of tolerance to plants. Among the TFs, NAC transcription factors were shown to provide enhanced abiotic stress tolerance by regulating a wide array of stress related genes [27, 32, 47–49].
Recently, through various genome-wide sequencing and gene expression profiling experiments, several NAC TF family members have been identified and characterized in Arabidopsis, rice, wheat and other plants [50–57]. It has been demonstrated that over-expression of stress-responsive NAC TFs can significantly improve abiotic stress tolerance in plants [18, 26, 34, 35, 58, 59]. In this context, identification and characterization of novel stress responsive NAC TFs from resilient crop species like finger millet can provide greater insight into this unique group of transcription factors. In our previous study, a salinity responsive NAC67 transcription factor (homologous to rice NAC TF; LOC_Os03g60080) was identified through RNA-sequencing in a set of contrasting finger millet genotype differing for their degree of salinity tolerance . In the present study, full length cDNA encoding stress inducible NAC67 TF of finger millet was isolated, cloned and characterized. Results of this study confirmed that EcNAC67 might act as a key TF in imparting abiotic stress tolerance.
Cloning, sequencing and analysis of deduced amino acids sequence of EcNAC67 suggested that it shared significant similarity with already reported NAC67 family of monocots i.e. S. italica, S. bicolor and Z. mays. Phylogenetic analysis showed genetic relatedness of EcNAC67 with NAC67 of S. italica (84 %) and S. bicolor SNAC1 (81.7 %). EcNAC67 was found to possess highly conserved N-terminal DNA binding domain when compared with other stress responsive NAC TFs and highly variable C-terminal transcriptional activation domain.
Among several NAC TFs reported in rice, OsSNAC1 , SNAC2  and Arabidopsis RD26 , ANAC019, ANAC055, and ANAC072 , very few were found to be stress-responsive. In this study, expression analysis of EcNAC67 revealed its salinity stress responsiveness where EcNAC67 transcripts were highly up-regulated under long term high salinity stress in leaves, roots and shoots of salinity tolerant finger millet genotype Trichy 1, indicating its involvement in salinity stress tolerance.
Ultimate aim of this study was to validate the role of salinity responsive EcNAC67 from finger millet for abiotic stress tolerance in rice. As there were several reports indicating that constitutive over-expression of transcription factors resulted in growth retardation in transgenic plants, EcNAC67 was cloned under the control of a stress inducible (RD29) promoter and used for transforming an elite rice cultivar ASD16. Putative transgenic lines were identified through PCR analysis and confirmed through southern hybridization (Fig. 3). In T1 generation, positive plants were identified through PCR analysis and the results showed the segregation of transgene in 3:1, following typical Mendelian segregation pattern. PCR positive T1 transgenic lines were forwarded to T2 generation and evaluated for their tolerance against salinity at both germination and vegetative stage along with non-transgenic controls. At germination stage, progenies of transgenic ASD16 lines over-expressing EcNAC67 showed better germination and shoot/root length under salinity stress (Fig. 5).
Transgenic lines exhibited better growth rate and found to possess lesser number of dried leaves when compared to the non-transgenic plants during salinity. All the transgenic lines were found to exhibit enhanced level of salinity tolerance when compared to non-transgenic ASD16 plants by recording lesser reduction in their root, shoot and total biomass during salinity. Similar observations were made by Liu et al. , where transgenic cotton plants over-expressed with OSNAC1 were found to accumulate more biomass when compared to the non-transgenic controls. There was no significant difference between the transgenic and non-transgenic ASD16 plants in their root length/growth under normal condition(s), but all the transgenic lines were found to possess significantly longer roots than non-transgenic ASD16 plants during salinity stress. Furthermore, all the transgenic ASD16 plants were found to recover rapidly from salinity induced injuries when compared to the non-transgenic plants when both were allowed to grow under normal conditions after 33 days of stress (Fig. 6c).
All the transgenic rice lines (T2) were evaluated for their drought tolerance ability in comparison with non-transgenic ASD16 plants. Non-transgenic ASD16 plants started showing drought induced leaf rolling and drying symptoms on 6th day after stress whereas transgenic lines showed much delayed leaf-rolling symptoms (Fig. 7). RWC (%) is a widely used parameter to determine the plant’s internal water status and is a direct reflection of the ability of the genotypes to retain turgidity during water deficit conditions . In this study, all the transgenic rice lines engineered with EcNAC67 were found to retain 17-22 % higher RWC than the non-transgenic ASD16 plants at more or less equal intensity of drought (Fig. 8). This enhanced water retention capacity of EcNAC67 engineered transgenic rice lines might be due to better stomatal regulation during dehydration similar to the previous reports on SNAC1  and TaNAC67 . Further experiments are needed in this aspect to understand the precise role of EcNAC67 in providing enhanced tolerance against drought stress. Further, all the transgenic lines were found to show better and rapid recovery upon rehydration when compared to non-transgenic ASD16 suggesting multiple roles of EcNAC67 in protecting cellular proteins/enzymes/organelles during dehydration apart from regulation of water loss through stomata. Grain yield of both non-transgenic and transgenic ASD16 lines were found to be significantly affected by intermittent drought but the % reduction in yield was much higher in non-transgenic ASD16 than most of the transgenic lines. Lesser reduction in the grain yield of transgenic rice plants may be due to relatively lesser reduction in spikelet fertility (11.7 % to 21.84 %) than that of their non-transgenic counterpart where it showed around 40.8 % spikelet sterility. Drought stress negatively affects successful pollination and fertilization by decreasing the amount of viable pollen thus leading to spikelet sterility  indicating that non-transgenic ASD16 had undergone severe drought stress during panicle/pollen development as compared to transgenic lines at equal intensity of soil moisture stress.
The qRT-PCR analysis of the EcNAC67 transgene did not show any detectable levels under control condition as it’s was driven by a stress inducible promoter, RD29 (Data not shown). qRT-PCR analysis of transgene expression in the transgenic lines revealed that two lines possessing single copy insertions (EcNAC67-E1 and EcNAC67-E4) were found to exhibit maximum level of expression under both salinity and drought stressed conditions (Fig. 4) as compared to the lines having multiple copies of transgene which may be attributed to co-suppression of transgene expression in case of multiple copy events as reported earlier . Abundance of ECNAC67 transcripts was found to be highly correlated with root proliferation and development in the transgenic plants suggesting the probable role of EcNAC67 in root development as reported in few other studies viz., EcNAC1 in tobacco , AtNAC2 and SNAC1 in rice [33, 47].
Results of this study indicated that EcNAC67 can serve as a novel source for engineering salinity/drought tolerance in crop plants. Even though growth retardation is a common adverse physiological disturbance reported in transgenic plants over-expressing TFs, morphological and agronomical characters of EcNAC67 engineered rice plants were comparable to non-transgenic ASD16 plants under normal well-watered conditions. No adverse effects were noticed in terms of plant height, panicle length and grain yield/plant. This shows the practical applicability of EcNAC67 for genetic improvement of abiotic stress tolerance in rice and other agricultural crops of economic importance.
Availability of data and materials
Nucleotide sequence of the reported candidate gene EcNAC67 from finger millet (Eleusine coracana L.) is available in the NCBI-GenBank database (Accession # KU500625). All other supporting data are included as additional files.
Wang W, Vinocur B, Altman A. Plant responses to drought, salinity and extreme temperatures: towards genetic engineering for stress tolerance. Planta. 2003;218(1):1–14.
Gepts P. A Comparison between crop domestication, classical plant breeding, and genetic engineering. Crop Sci. 2002;42(6):1780–90.
Athar HR, Ashraf M. Strategies for crop improvement against salinity and drought stress: An overview. In: Ashraf M, Ozturk M, Athar HR, editors. Salinity and Water Stress, vol. 44. Netherlands: Springer; 2009. p. 1–16.
Mittler R, Blumwald E. Genetic engineering for modern agriculture: challenges and perspectives. Ann Rev Plant Biol. 2010;61:443–62.
Singla-Pareek SL, Reddy M, Sopory SK. Transgenic approach towards developing abiotic stress tolerance in plants. Proc Indian Nat Sci Acad. 2001;B67(5):265–84.
Yadav RC, Solanke AU, Kumar P, Pattanayak D, Yadav NR, Kumar PA. Genetic engineering for tolerance to climate change-related traits. In: Kole C, editor. Genomics and breeding for climate-resilient crops. Berlin Heidelberg: Springer-Verlagr; 2013. p. 285–330.
Turan S, Cornish K, Kumar S. Salinity tolerance in plants: breeding and genetic engineering. Aus J Crop Sci. 2012;6(9):1337–48.
Lata C, Yadav A, Prasad M. Role of plant transcription factors in abiotic stress tolerance. In: Shanker A, VenkateswarluB,editors. Abiotic stress response in plants - Physiological, biochemical and genetic perspective. InTech. 2011. doi:10.5772/23172 .
Khong G, Richaud F, Coudert Y, Pati PK, Santi C, Périn C, et al. Modulating rice stress tolerance by transcription factors. Biotechnol Genet Eng Rev. 2008;25:381–404.
Nakashima K, Takasaki H, Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K. NAC transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012;1819(2):97–103.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53:247–73.
Yadav NR, Taunk J, Rani A, Aneja B, Yadav RC. Role of transcription factors in abiotic stress tolerance in crop plants. In: Tuteja N, Gill SS, editors. Climate change and plant abiotic stress tolerance. Germany: Wiley-VCH Verlag GmbH & Co; 2013. p. 605–40.
Agarwal P, Agarwal P, Reddy MK, Sopory S. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Rep. 2006;25(12):1263–74.
Shen H, Liu C, Zhang Y, Meng X, Zhou X, Chu C, et al. OsWRKY30 is activated by MAP kinases to confer drought tolerance in rice. Plant Mol Biol. 2012;80(3):241–53.
Xie Q. Arabidopsis NAC1 transduces auxin signal downstream of TIR1 to promote lateral root development. Genes Dev. 2000;14(23):3024–36.
Duval M, Hsieh TF, Kim SY, Thomas TL. Molecular characterization of AtNAM: a member of the Arabidopsis NAC domain superfamily. Plant Mol Biol. 2002;50(2):237–48.
Ernst HA, Nina Olsen A, Skriver K, Larsen S, Lo-Leggio L. Structure of the conserved domain of ANAC, a member of the NAC family of transcription factors. EMBO Rep. 2004;5(3):297–303.
Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, Xiong L. Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Nat Acad Sci USA. 2006;103(35):12987–92.
Jin J, Zhang H, Kong L, Gao G, Luo J. PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 2014;42(Database issue):D1182–7.
Souer E, van Houwelingen A, Kloos D, Mol J, Koes R. 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. 1996;85(2):159–70.
Zhong R, Demura T, Ye ZH. SND1, a NAC domain transcription factor, is a key regulator of secondary wall synthesis in fibers of Arabidopsis. Plant Cell. 2006;18(11):3158–70.
Zhong R, Richardson E, Ye ZH. Two NAC domain transcription factors, SND1 and NST1, function redundantly in regulation of secondary wall synthesis in fibers of Arabidopsis. Planta. 2007;225(6):1603–11.
Yang SD, Seo PJ, Yoon HK, Park CM. The Arabidopsis NAC transcription factor VNI2 integrates abscisic acid signals into leaf senescence via the COR/RD genes. Plant Cell. 2011;23(6):2155–68.
Sun L, Li D, Zhang H, Song F. Functions of NAC transcription factors in biotic and abiotic stress responses in plants. Hereditas(Beijing). 2012;34(8):993–1002.
Nuruzzaman M, Sharoni AM, Kikuchi S. Roles of NAC transcription factors in the regulation of biotic and abiotic stress responses in plants. Front microbiol. 2013;4:1–16.
Ramegowda V, Senthil-Kumar M, Nataraja KN, Reddy MK, Mysore KS, Udayakumar M. Expression of a finger millet transcription factor, EcNAC1, in tobacco confers abiotic stress-tolerance. PloS one. 2012;7(7), e40397.
Tran LS, Nishiyama R, Yamaguchi-Shinozaki K, Shinozaki K. Potential utilization of NAC transcription factors to enhance abiotic stress tolerance in plants by biotechnological approach. GM crops. 2010;1(1):32–9.
Fujita M, Fujita Y, Maruyama K, Seki M, Hiratsu K, Ohme-Takagi M, et al. A dehydration-induced NAC protein, RD26, is involved in a novel ABA-dependent stress-signaling pathway. Plant J. 2004;39(6):863–76.
Jiang Y, Deyholos MK. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006;6(1):25.
Fang Y, You J, Xie K, Xie W, Xiong L. Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol Genet Genomics. 2008;280(6):547–63.
Mochida K, Yoshida T, Sakurai T, Yamaguchi-Shinozaki K, Shinozaki K, Tran LS. In silico analysis of transcription factor repertoire and prediction of stress responsive transcription factors in soybean. DNA Res. 2009;16(6):353–69.
Nakashima K, Tran LS, Van-Nguyen D, Fujita M, Maruyama K, Todaka D, et al. Functional analysis of a NAC‐type transcription factor OsNAC6 involved in abiotic and biotic stress‐responsive gene expression in rice. Plant J. 2007;51(4):617–30.
Liu G, Li X, Jin S, Liu X, Zhu L, Nie Y, et al. Overexpression of Rice NAC gene SNAC1 improves drought and salt tolerance by enhancing root development and reducing transpiration rate in transgenic Cotton. PLoS ONE. 2014;9(1):1–10.
Jeong JS, Kim YS, Baek KH, Jung H, Ha SH, Do-Choi Y, Kim M, Reuzeau C, Kim JK. Root-specific expression of OsNAC10 improves drought tolerance and grain yield in Rice under field drought conditions. Plant Physiol. 2010;153(1):185–97.
Mao X, Chen S, Li A, Zhai C, Jing R. Novel NAC Transcription Factor TaNAC67 confers enhanced multi-abiotic stress tolerances in Arabidopsis. PLoS ONE. 2014;9(1):1–10.
Mehta PA, Sivaprakash K, Parani M, Venkataraman G, Parida AK. Generation and analysis of expressed sequence tags from the salt-tolerant mangrove species Avicennia marina (Forsk) Vierh. Theor Appl Genet. 2005;110(3):416–24.
Mishra RN, Reddy PS, Nair S, Markandeya G, Reddy AR, Sopory SK, et al. Isolation and characterization of expressed sequence tags (ESTs) from subtracted cDNA libraries of Pennisetum glaucum seedlings. Plant Mol Biol. 2007;64(6):713–32.
Garg R, Verma M, Agrawal S, Shankar R, Majee M, Jain M. Deep transcriptome sequencing of wild halophyte Rice, Porteresia coarctata, provides novel insights into the salinity and submergence tolerance factors. DNA Res. 2014;21(1):69–84.
Devos K, Gale M. Comparative genetics in the grasses. Plant Mol Biol. 1997;35(1–2):3–15.
Shailaja H, Thirumeni S. Evaluation of salt-tolerance in finger millet (Eleusine coracana) genotypes at seedling stage. Indian J Agric Sci. 2007;77(10):672–4.
Agarwal S, Kumar A, Singh PK, Singh A. Responses of some genotypes of finger millet (Eleusine coracana Gaertn.) for their salt tolerance. Int J Curr Res. 2011;3:45–50.
Rahman H, Jagadeeshselvam N, Valarmathi R, Sachin B, Sasikala R, Senthil N, Sudhakar D, Robin S, Muthurajan R. Transcriptome analysis of salinity responsiveness in contrasting genotypes of finger millet (Eleusine coracana L.) through RNA-sequencing. Plant Mol Biol. 2014;85(4–5):485–503.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔC T method. Methods. 2001;25(4):402–8.
Weigel D, Glazebrook J. Transformation of Agrobacterium using the freeze-thaw method. Cold Spring Harbor Protocols. 2006 (7): doi:10.1101/pdb.prot4666.
Hiei Y, Komari T. Agrobacterium-mediated transformation of rice using immature embryos or calli induced from mature seed. Nat Protoc. 2008;3(5):824–34.
Southern E. Southern blotting. Nat Protoc. 2006;1(2):518–25.
He XJ, Mu RL, Cao WH, Zhang ZG, Zhang JS, Chen SY. AtNAC2, a transcription factor downstream of ethylene and auxin signaling pathways, is involved in salt stress response and lateral root development. Plant J. 2005;44(6):903–16.
Hu H, You J, Fang Y, Zhu X, Qi Z, Xiong L. Characterization of transcription factor gene SNAC2 conferring cold and salt tolerance in rice. Plant Mol Biol. 2008;67(1–2):169–81.
Puranik S, Sahu PP, Srivastava PS, Prasad M. NAC proteins: regulation and role in stress tolerance. Trends Plant Sci. 2012;17(6):369–81.
Nuruzzaman M, Manimekalai R, Sharoni AM, Satoh K, Kondoh H, Ooka H, et al. Genome-wide analysis of NAC transcription factor family in rice. Gene. 2010;465(1–2):30–44.
Puranik S, Sahu PP, Mandal SN, Suresh BV, Parida SK, Prasad M. Comprehensive genome-wide Survey, genomic constitution and expression profiling of the NAC transcription factor family in Foxtail Millet (Setaria italica L.). PLoS ONE. 2013;8(5):1–16.
Singh AK, Sharma V, Pal AK, Acharya V, Ahuja PS. Genome-wide organization and expression profiling of the NAC transcription factor family in Potato (Solanum tuberosum L.). DNA Res. 2013;20(4):403–23.
Chen Z, Guo X, Chen Z, Chen W, Liu D, Zheng Y, et al. Genome-wide characterization of developmental stage-and tissue-specific transcription factors in wheat. BMC Genomics. 2015;16(1):1–15.
Huang H, Wang Y, Wang S, Wu X, Yang K, Niu Y, et al. Transcriptome-wide survey and expression analysis of stress-responsive NAC genes in Chrysanthemumlavandulifolium. Plant Sci. 2012;193–194:18–27.
Ma J, Wang F, Li M, Jiang Q, Tan G, Xiong A. Genome wide analysis of the NAC transcription factor family in Chinese cabbage to elucidate responses to temperature stress. Sci Hortic. 2014;165:82–90.
Pinheiro GL, Marques CS, Costa MDBL, Reis PAB, Alves MS, Carvalho CM, et al. Complete inventory of soybean NAC transcription factors: Sequence conservation and expression analysis uncover their distinct roles in stress response. Gene. 2009;444(1–2):10–23.
Shiriga K, Sharma R, Kumar K, Yadav SK, Hossain F, Thirunavukkarasu N. Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize. Meta Gene. 2014;2:407–17.
You J, Zhang L, Song B, Qi X, Chan Z. Systematic analysis and identification of stress-responsive genes of the NAC gene family in Brachypodium distachyon. PLoS ONE. 2015;10(3), e0122027.
Zheng X, Chen B, Lu G, Han B. Overexpression of a NAC transcription factor enhances rice drought and salt tolerance. Biochem Biophys Res Commun. 2009;379(4):985–9.
Tran LS, Nakashima K, Sakuma Y, Simpson SD, Fujita Y, Maruyama K, et al. 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. 2004;16(9):2481–98.
Alizadeh A. Soil, water, plants Relationship. 3rd ed. Iran: Emam Reza University Press; 2002.
Alqudah AM, Samarah NH, Mullen RE. Drought stress effect on crop pollination, seed set, yield and quality. In: Lichtfocus E, editor. Alternative Farming Systems, Biotechnology, Drought Stress and Ecological Fertilisation. Netherlands: Springer; 2011. p. 193–213.
Daniel RG. Controlling gene expression in transgenics. Curr Opin Plant Biol. 1998;1(2):166–72.
The financial support of Department of Biotechnology, Government of India, New Delhi (Grant Number BT/PR-10482/AGR/02/564/2008) is greatly acknowledged. The authors are grateful to University Grants Commission, Govt. of India, New Delhi for providing fellowship to the first author. Technical help rendered by Dr. M. Parani, Professor and Head, Dept. of Genetic Engineering, SRM University, Chennai, India by providing the plant transformation vector pCAMBIA1300 with stress inducible RD29 promoter and Dr. S. Robin, Professor and Head, Dept. of Rice, Tamil Nadu Agricultural University, Coimbatore, India in reviewing the manuscript is greatly acknowledged.
The study was funded by Department of Biotechnology, Government of India, New Delhi (Grant Number BT/PR-10482/AGR/02/564/2008). Publication charges for this article were met from author’s institutional resources.
This article has been published as part of BMC Biotechnology Volume 16 Supplement 1, 2016: Proceedings of the Indian Genetics Congress 2015: Biotechnology. The full contents of the supplement are available online at http://bmcbiotechnol.biomedcentral.com/articles/supplements/volume-16-supplement-1.
The authors declare that they have no competing interests.
HR and VR were involved in generation and evaluation of transgenic lines; JN did the bioinformatic analysis of EcNAC67; SD and RM were involved in designing, planning and execution of the work and manuscript preparation. All authors read and approved the final manuscript.
Table showing the list of primers used in the study. (PDF 42 kb)
Secondary structure of deduced EcNAC67 protein sequence predicted by PSIPRED protein structure prediction server. (PDF 283 kb)
Table showing the estimates of Evolutionary Divergence among NAC homologs. Analyses were conducted using the Poisson correction model. All positions containing gaps and missing data were eliminated. There were a total of 241 positions in the final dataset. Evolutionary analyses were conducted in MEGA6. (PDF 191 kb)