Overexpression of cotton GhNAC072 gene enhances drought and salt stress tolerance in transgenic Arabidopsis
BMC Genomics volume 23, Article number: 648 (2022)
Crops face several environmental stresses (biotic and abiotic), thus resulting in severe yield losses. Around the globe abiotic stresses are the main contributors of plant damages, primarily drought and salinity. Many genes and transcription factors are involved in abiotic and biotic stress responses. NAC TF (Transcription Factors) improves tolerance to stresses by controlling the physiological and enzyme activities of crops.
In current research, GhNAC072 a highly upregulated TF in RNA-Seq was identified as a hub gene in the co-expression network analysis (WGCNA). This gene was transformed to Arabidopsis thaliana to confirm its potential role in drought and salt stress tolerance. Significant variations were observed in the morpho-physiological traits with high relative leaf water contents, chlorophyll contents, higher germination and longer root lengths of the overexpressed lines and low excised leaf loss and ion leakage as compared to the wildtype plants. Besides, overexpressed lines have higher amounts of antioxidants and low oxidant enzyme activities than the wildtype during the period of stress exposure.
In summary, the above analysis showed that GhNAC072 might be the true candidate involved in boosting tolerance mechanisms under drought and salinity stress.
Crops play a vital part in our lives but, face several environmental stresses regularly. The plant perceives and transduces the stress signals, which leads to the expression of functional proteins that safeguard the plant. Genes that respond to stresses are primarily controlled by transcription factors, which are vital in plant resistance to the changing environmental conditions . A transcription factor (TF) is a type of important regulatory protein that regulates the expression level of a gene . It governs chromatin and transcription by recognizing certain DNA sequences, forming a complex mechanism that directs genome expression .
Abiotic stress is the harmful influence of abiotic factors on living things in a definite location. Drought, high temperature, cold, salt stress and other ecological extremes are among the stresses. Crop losses are primarily caused by salinity and drought around the globe . Plant cells produce a huge quantity of reactive oxygen species (ROS) when they are under stress. Under these circumstances, the activity of superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD) enzymes will rise to some extent to eliminate toxic compounds from the plant and preserve normal development .
Plant transcription factors include more than 80 different TF families. NAC, MYB, WRKY, bZIP, and ERF/DREB are only a few TFs that are involved in abiotic and biotic stress responses . The NAC TF is part of one of the biggest families of TFs involved in crops response to abiotic stresses . A NAC TF has a two way regulatory approach: The first is transcriptional regulation, which includes phosphorylation and ubiquitination miRNAs, which can control TF expression at the protein level. The alternative method is to attach to the target mRNA and perform posttranscriptional regulation. NAC TF can impact the final output by binding to sequences and regulating target genes . NAC TF improved tolerance to stress by controlling the physiological and enzyme activities of crops. So far, this family has been studied in a wide range of crops. In Malus baccata, the transformation of the MbNAC25 gene in Arabidopsis improves cold and salinity tolerance . Drought tolerance was dramatically increased and lowered in transgenic rice plants via overexpressing and RNA interfering in the OsDRAP1 transcription factor .
Transformation of ZmNAC33 transcription factor from maize was done in A. thaliana. The germination percentage under ABA and osmotic stress during the germination stage are higher than that of wildtype seedlings. Compared with the low wild-type drought stress, the overexpression lines show a higher survival rate and higher antioxidant enzyme activity . JUNGBRUNNEN1NAC TF improved drought tolerance in Solanum lycopersicum . The mOsNAC2 overexpressing lines were much more drought and salinity tolerant than the wildtype, with a higher survival rate thanks to higher ABA levels . Similar research finding in rice indicated that the ONAC022oe plants had lower rates of water loss and transpiration, a lower proportion of opening stomata collected fewer Na+ in roots, and had higher levels of proline and soluble carbohydrates than wildtype seedlings during drought and salinity stress . Under drought and salt stress, overexpressed lines had higher expressions of NtARF1, NtARF2, and NtARF8, implying that overexpression of SlNAC35 aided root growing by engaging auxin signaling and modulating NtARF expression in tomato . The ABA-dependent signaling system improves drought and cold stress tolerance in overexpressed lines by transforming the MLNAC5 gene from Miscanthus lutarioriparius as a stress responsive NAC TF .
Several experiments was done specifically on NAC072 TF in different crop species. A NAC TF from Oxytropis ochrocephala, OoNAC72 enhanced drought, salt stress tolerance, and ABA-dependent process regulation in Arabidopsis. OoNAC72 overexpression also increased the expression of stress-responsive genes like RD29A, RD29B, RD26, LEA14, ANACOR19, ZAT10, PP2CA, and NCED3 . In line with the above finding, overexpression of RcNAC72, a NAC TF from Rosa chinensis, improved drought tolerance and ABA sensitivity in Arabidopsis. Furthermore, in rose leaves, silencing RcNAC72 decreases tolerance to water loss stress and rehydration. The regulation mechanism of RcNAC72 via the ABA-dependent signaling pathway and the DRE/CBF-COR pathway in response to drought stress identified . There is an interaction between NAC072 and ABA-responsive element binding factor 3 (ABF3) in Arabidopsis, which acts as a positive regulator of ABA-responsive gene expression. ABF3 upregulates NAC072 expression in the ABA response, and NAC072 protein also interacts with ABF3 .
Transcriptome analysis in cotton species was very crucial for mining of genes involved in abiotic and biotic stresses. RNA-Seq analysis was done to identify differentially expressed genes in cadmium stress tolerance , under NaHCO3 alkaline stress , in cotton boll weevil pests , in Aspergillus flavus fungi responsible for Aflatoxin production , in fiber growth and development , in anther development and male sterility traits , and DEGs linked to Na2SO4 tolerance . To better understand the role of the NAC genes in drought and salinity stress, we cloned the NAC transcription factor gene GhNAC072 from G. hirsutum L. The function of this TF was investigated and recognized as a potential candidate for improved cotton cultivations under drought and salt stress.
Identification of NAC genes in G. hirsutum, G. arboreum and G. raimondii
We identified a total of 608 NAC TF in the three Gossypium species, 305 in G. hirsutum with 152 genes in GhAt subgenome, 144 genes in GhDt subgenome and 9 genes in scaffold region, 150 in G. arboreum and 153 in G. raimondii according to CottonFGD (www.cottonfgd.org) database. 295, 44, 123, 119 and 27 NAC genes were found in group I, Group II, Group III, Group IV and Group V respectively. The highest number of genes were found in group I with 153 NAC genes in G. hirsutum, 72 in G. arboreum and 70 in G. raimondii while the lowest scored in group V with 13, 7 and 7 NAC genes in G hirsutum, G. arboreum and G. raimondii respectively (Table S1), (Fig. 1).
Chromosomal mapping and gene ontology analysis
Genome wide analysis of the three Gossypium species revealed that, the NAC TF were distributed among all thirteen chromosomes. Both At05 and At11 chromosomes has the maximum number of genes (15 genes) from the At subgenome individually, whereas the lowest number of genes (6 genes) were located on chromosomes At04 and At06. Similarly, Dt11 and Dt05 had the most gene loci with 15 and 12 genes in the Dt Subgenome, whereas Dt03, Dt08, and Dt10 had 7, 7, and 5 genes, respectively. The remaining chromosomes had 6 and 15 genes on them. The gene distribution arrangement was similar among two diploid Gossypium species, G. arboreum and G. raimondii. The highest gene loci in G. arboreum were found on chromosomes A11 and A05, with 17 and 15 genes, respectively, whereas in G. raimondii, the highest gene loci were recorded on chromosomes D07 and D09, with 17 and 16 genes, respectively, while chromosomes D02, D03, and D13 had 3, 7 and 9 genes, respectively. The scaffold category harbored 33 NAC genes together (Fig S1).
There was an odd distribution of genes between the chromosomes of G. hirsutum and both diploid species, G. arboreum and G. raimondii. In G. hirsutum At subgenome and G. arboreum 5th and 9thchromosomes showed the same number of genes with 15 and 9 genes respectively. Similarly, chromosomes 2, 3, and 5 harbor the same number of genes in the GhDt subgenome and G. raimondii with 11, 7, and 14 genes respectively. Only six chromosomal pairs from the GhAt/A and GhDt/D contrasts were discovered to have a similar number of genes (Table 1). This depicts the rate of gene loss during evolution, as well as the role of progenitors in the development of tetraploid cotton species . The unequal distribution of genes on the cotton genome's chromosome could be due to gene replication or a partial segment of gene duplication that happened during the plant's long evolutionary history. With each duplication, cotton's whole gene sequence was doubled, and the extra genes were recombined or eliminated with time .
The Gene Ontology (GO) project was a large bioinformatics endeavor that aimed to create a computer description of how genes encode activities at the biological, molecular, and cellular system levels . In G. hirsutum, the biological function (GO:0,008,150) includes biological regulation, regulation of biological process, cellular process, and metabolic process that regulates transcription and DNA-dependent functions. Similarly, in molecular function—GO:0,003,674, the significant role was binding of nucleic acid and DNA functions (Fig S2).
In both diploid species the significant functions were similar to G. hirsutum, biological and molecular functions. In G. arboreum the biological process function—GO:0,008,150 deals with metabolic process, biological regulation, regulation of biological process, cellular process, developmental process, single-organism process, multicellular organismal process, reproduction, cellular component organization or biogenesis and response to stimulus functions. In molecular activity—GO:0,003,674, binding and nucleic acid binding transcription factor activity functions. For G. raimondii, the biological process—GO:0,008,150 mainly on metabolic process, cellular process, biological regulation, and regulation of biological process functions. In molecular function—GO:0,003,674, Binding which regulates organic cyclic compound binding and heterocyclic compound binding that leads to nucleic acid binding for final DNA binding processes. In all the three cotton species, the molecular function was similarly related to the binding of nucleic acids and there is no significant cellular function in the GO category.
Gene structure and motif identification analysis
In gene structure analysis of NAC TFs, high variability was found between all genes. In G. hirsutum, the highest exon–intron ratios were scored for Gh_D02G0122 with 9 exons and 8 introns, whereas just one exon was present in, Gh_A01G0806, Gh_A01G0131,Gh_D11G0608, Gh_A11G0524, Gh_Sca077544G01, Gh_Sca121194G01, Gh_Sca111532G01, Gh_Sca140422G01, Gh_Sca043542G01 and Gh_Sca083050G01 genes (Fig S1 A).
In G. arboreum, Ga03G0113 has the highest numbers of exons and introns (14 exons, 13 introns) followed by Ga12G0369 (9 exons, 8 introns) whereas, the lowest exon–intron ratio of 2:1 was found in Ga09G2767, Ga01G0440, Ga14G2861, Ga06G2355, and Ga13G0670 genes. Similarly, In G. raimondii Gorai.008G261400 (9 exons, 8 introns), possesses the highest number of exons and introns, followed by Gorai.005G013300, Gorai.007G114500, and Gorai.006G034600 with (7 exons, 8 introns), while a single exon was found in Gorai.002G019700 (Fig S3 C).
To better understand the structural evolution of NAC proteins, the patterns of motifs were investigated. In the three Gossypium species, MEME analysis (http://meme-suite.org/) showed a total of 20 unique motifs. Based on the known motifs, In G. hirsutum motifs 3, 4, 1, 7, 2, and 5 were conserved while motifs 6, 14, 15, and 17 were diversified once. In G. arboreum, motifs 2, 5, 1, 4, 3, 6, and 8 were conserved with 16, 18, 19, and 20 diversified. Similarly, in G. raimondii, motifs 3, 4, 1, 12, 2, 6, and 5 were conserved, however, 17, 18, 19, and 20 were found as diversified motifs. All three cotton species shared almost the same conserved and diversified motifs (Fig S3).
Co-expression network analysis of NAC TFs for hub gene identification
Co-expression network analysis was performed using 193 DEGs selected from the RNA-Seq of Gossypium hirsutum races related to drought stress tolerance. Correlation based relationships were plotted using a Pearson correlation coefficient greater than 0.99 (Fig. 2). The network has a total of 193 nodes connected in it, with 1768 edges in the leaves and 1801 in the roots. In case of leaves we found, 973 positive and 795 negative correlations, whereas in roots 993 were positive and 808 were negatively correlated. A threshold level of > 20 edges was considered as hub genes from the analysis. Owing to this, 14 genes in leaf tissues and 17 in the roots were considered as hub genes. Among the hub genes, highest correlation was recorded in Gh_A05G3322 with 27 and Gh_A03G0887 scored 29 correlations in leaf and root tissues respectively. The candidate gene Gh_D01G0514 (GhNAC072) was included in both correlations as a hub gene.
The GhNAC072 gene has a protein length of 346 aa and a CDS length of 1041 bp. There are 3 exons, 2 introns in GhNAC072 gene with a molecular weight of 38.4 KDa, and an isoelectric point of 9, with -0.573 grand average of hydropathy values (Table 2).
Expression of NAC genes in leaves and roots under drought and salinity stress
Already available RNA-Seq data (PRJNA663204) of three cotton lines from G. hirsutum was used to detect the expression variations of NAC genes in leaf and root tissues at different time intervals. The GhNAC genes exhibited differential expression patterns in both tissues. According to the expression analysis most of GhNAC genes were upregulated in leaf tissues than roots tissues collected from the three cotton lines. There were no significant variations observed from 24 to 48 h. GhNAC072which was also selected as a hub gene via WGCNA analysis showed an upregulated expression in leaf tissues of all Gossypium species (Fig. 3, Fig. 4B).
Similarly, the RNA-Seq analysis of NAC genes also showed consistent expressions in Mariegalante-85, Upland cotton, and Latifolium-40 under drought stress in roots tissues. The gene expression profile in both tissues leads to selection of Gh_D01G0514 as a candidate gene and increased our confidence (Fig. 4A).
Selection of overexpressed lines from T2 generation
After transformation, we found seven overexpressed positive lines at the T2 generation stage. We performed PCR and RT-qPCR analysis to check the expression levels of all the lines. OE-1, OE-3, and OE-5 overexpressed lines were selected for further experiments due to their highest expression levels (Fig. 5).
Evaluation of physiological parameters under drought and salt stress
A significant variation (p < 0.05) was observed among wildtype and overexpressed lines in excised leaf water loss. No significant variations were observed under control conditions, whereas, after drought and salt application, the wildtype seedlings had the highest leaf water loss and ion leakage as compared to the overexpressed lines. Significant differences (P < 0.05) were observed for relative leaf water contents and chlorophyll contents between the wildtype and overexpressed lines (Fig. 6).
Determination of oxidant and antioxidant enzymes activities
Our results showed that the enzyme activities have been significantly influenced by drought and salt treatments. The activities of CAT and SOD were observed highest in the overexpressed lines than the wildtype. Whereas, the levels of oxidant enzymes, MDA and H2O2 were lower in the overexpressed lines than the wildtype under drought and salinity stress. Current results confirmed that overexpressed lines face a minimal damage as compared to wildtype under stress conditions (Fig. 7).
GhNAC072 enhances germination rate and root length in overexpressed lines
Significant variations (P < 0.05) were observed for germination percentage in both, wildtype and overexpressed lines during drought and salinity stress conditions. The germination percentage of the overexpressed lines was higher than that of wildtype under 100 mM, 200 mM, and 300 mM mannitol concentrations. Under 300 mM mannitol concentration, the germination percentages of < 47% and > 67% were observed in case of wildtype and overexpressed lines respectively (Fig. 8). Similarly, under 200 mM NaCl concentration, the germination rate in wildtype was 40% but for overexpressed lines more than 72% germination rate was recorded in the current study.
Drought and salt stress treatments have a considerable impact on root length. There were no significant differences in root length between wildtype and overexpressed lines under normal circumstances. The overexpressed lines, on the other hand, had longer roots per line under both drought and salt treatments with 300 mM mannitol and 200 mM NaCl, respectively. This confirms that GhNAC072plays a major role in drought and salinity stress tolerance (Fig. 9).
Evaluation of abiotic stress responsive genes under drought and salt stress
The four responsive genes APF4, SOS1, RAB18, and RD22 were significantly upregulated during drought and salinity stress situations in the wildtype in the reverse transcription quantitative (RT-qPCR) analysis (Fig. 10). Current results indicate that the overexpression of GhNAC072 has a valuable effect on the upregulation of stress responsive genes expressions in A. thaliana, indicating that regulation of GhNAC072 might be crucial in the abiotic stress tolerance of crops.
Plants are vulnerable to several types of abiotic stresses, the most notable of which is drought stress. Drought has currently gotten a lot of attention due to its severe impact on crop production . Similarly, salinity leads to a significant loss in crop output every year. The lower osmotic potential of soil solutes in high salinity soils might result in a huge loss of crop production [29, 30].
The NAC gene family has 305 genes in G. hirsutum, 150 genes in G. arboreum, and 153 genes in G. raimondii. Both the exons and introns take part in the protein formation. Both follow the DNA packaging mechanism to fit inside a cell. Gene structure arrangements are thus made by comparing exons conferring to their size and reading frame . Exon–intron ratios in G. hirsutum and G. arboreum were lower, but there was upstream/downstream attachment at the beginning and end of the structures in G. raimondii. Many stress responsive genes that play a crucial part in boosting osmotic and salinity stress in cotton have low intron interruption . In the three Gossypium species, the distribution of genes was unequal, although six chromosomal pairs were discovered to possess a similar number of genes. Thus, showing the extent of gene loss during the evolution process by observing the contribution of progenitors in developing tetraploid Gossypium species throughout evolution . The unequal distribution of genes on the cotton genome's chromosomes could be because of gene repetition or a partial segment of gene duplication that happened during the plant's long evolutionary history. With each duplication, cotton's whole gene sequence was doubled, and the extra genes were recombined or eliminated throughout time .
NAC genes positively regulate abiotic stress related genes during harsh environmental conditions. Previous studies identified CaNAC46 in Capsicum annuum , ATNAC3 in Solanum lycopersicum , PgNACs in Pennisetum glaucum  improves drought and salinity tolerance, AmNAC11 from Ammopiptanthus mongolicus boosts drought and cold stress tolerance in A. thaliana, OsNAC2 and ONAC022in rice positively controls drought and salinity stress via ABA synthesis mechanisms [11, 12].
In case of physiological traits, the wildtype seedlings had the highest excised leaf water loss and ion leakage after drought and salinity stress treatments relative to the overexpressed lines. Whereas, relative leaf water content and chlorophyll content for the tested treatments were recorded significantly higher in the overexpressed lines. When the crops were subjected to numerous stress situations, the overexpressed lines outperformed the wildtype in terms of stress resistance. The GhDTX/MATE transgenic lines scored the highest amount of leaf water content. However, low in leaf water loss and cell membrane stability over the wildtype Arabidopsis seedlings during drought, salinity, and cold stress . Similar to this finding, leaf water content and chlorophyll content were increased, and the leaf water loss and ion leakage exhibited a substantial decrease in overexpressed lines compared with the WT lines were reported in GhMPK3 genes from cotton . The overexpressed lines had the highest CAT and SOD concentrations but lower levels of MDA and H2O2 oxidant enzymes in relative to the wildtype in drought and salinity stress environments. This showed that overexpressed lines face low oxidative stress as compared to the wildtype under stress conditions. Under drought and salinity situations, the level of oxidants (MDA and H2O2) were found to be considerably greater in the wildtype and significantly lower in the transgenic lines. Antioxidant levels in transgenic lines were found to be higher than in wildtype . Alike findings were stated in previous works from GhTOM transgenic lines in improving salinity stress tolerance in A. thaliana .
GhNAC072 enhances the germination rate and root length capacity of A. thaliana during drought and salinity stress treatments. The germination percentage and root length of GhNAC072 overexpressed lines were higher than that of wildtype under different mannitol and NaCl concentrations. This showed that the GhNAC072 gene enhances the tolerance ability in A. thaliana seedlings during stress conditions. Under long drought and extreme salt stress, the CaNAC46 transgenic lines improve root length and lateral branches than wildtype . GmNAC8 transgenic lines demonstrated greater tolerance and survival rates under drought stress, but GmNAC8 VIGS data suggested the contrary . The profiling of MeNAC genes was studied in cassava under drought, ABA, salinity, cold, and H2O2 conditions, revealing that NACs may serve as a site of junction for many signaling pathways . Similarly, functional validation using transformation and silencing of transgenic lines demonstrate that overexpression of ONAC066 enhanced while knockdown of ONAC066 diminished oxidative stress tolerance in rice . GhNAC072 expression profiling with the stress responsive genesAPF4, SOS1, RAB18, and RD22 upregulated during drought and salt stress conditions, signifying that GhNAC072 may have a crucial role in the abiotic stress tolerance. In line with these results, the accumulation of reactive oxygen species was influenced by CaNAC46. Moreover, during drought and salinity stress, CaNAC46 increased the profiling of responsive genes RD29B, RD20, LDB18, P5CS, and IAA4 . In line with this, in the mOsNAC2 overexpressing lines, reverse transcription quantitative PCR study revealed a considerable boost in the profiling of the ABA biosynthesis genes OsNCED1and OsNCED3, as well as the stress responsive genes OsP5CS1, OsLEA3, and OsRab16 .
Herewith, we confirmed on the basis of current results that the cotton GhNAC072 genes plays a significant role in improving drought and salinity stress tolerance. Physiological traits, biochemical parameters, germination rate, and root length evaluation, and stress responsive genes expression analysis showed that GhNAC072 is a potential involved in boosting tolerance mechanisms under drought and salinity stress. In summary, these results offer data that can be used for applications of GhNAC072 function in plants and can lead to genetic enhancement of crops in the field of molecular breeding and crop selection.
Materials and methods
Arabidopsis thaliana Col-0 was used to study overexpression. The seeds were vernalized in dark for two days at 4 °C before being moved to the growth chamber. The plants were grown in a growth chamber with a temperature of 22 °C, a photoperiod of 16 h of light/8 h of darkness, and relative humidity of 80%. Then, the A. thaliana seedlings were grown in pots with a 1:1 soil mixture of enriched soil and vermiculite. T0 – T3 seeds were screened, and T3 progeny was used for the final treatment .
NAC family and phylogenetic tree analysis
CottonFGD database (www.cottonfgd.org) was used to identify the NAC genes in cotton species . Using the identified genes from G. hirsutum species, the expression analysis of the all genes obtained from the RNA-Seq (NCBI Accession: PRJNA663204) was performed to identify the candidate gene “GhNAC072”. The amino acid sequence of the genes of NAC family members in G. hirsutum, G. arboreum, and G. raimondii, were used to develop the phylogenetic tree. The protein sequences of the selected species were aligned via ClustalX and the tree was built by MEGA 7.0 applying the neighbor-joining (NJ) technique under 1000 bootstrap replications .
Chromosome mapping and gene ontology analysis
The GFF3 file from CottonFGD (www.cottonfgd.org) and the gene IDs were used to determine the distribution of NAC TF genes across all chromosomes of A, D, and AD cotton genomes. Then, using TBtools, the gene locations were shown on the chromosomes using data from GTF/GFF3 file. The functional taxonomy of genes were established using an online program AgriGO (www.bioinfo.cau.edu.cn/agriGO), which was based on cellular composition, biological processes, and molecular function .
Gene structure analysis and motif identification
To arrange the gene structure of the Gossypium species G. hirsutum, G. arboreum, and G. raimondii, CDS sequences and genomic DNA sequences of the genes were uploaded to gene structure display server (http://gsds.gao-lab.org/) website . Motif identification analysis was performed by uploading the protein sequence of the cotton species to the Multiple Expectation Maximization for Motif Elicitation-MEME website (http://meme-suite.org/tools/meme), thus using the Mast-XML file and gene IDs to visualize the graph via TBtools software .
Co-expression network analysis of genes
We used already available RNA-Seq data of leaf and root tissues under drought stress to check the expression patterns of NAC gene family. The links between genes involved in drought stress tolerance were examined by coexpression network analysis. A coexpression regulation network of the genes linked to drought stress tolerance was created using the Cytoscape software (version 3.7.2) . The threshold for the coexpression network map was set as p ≥ 0.99. The topological coefficient of each node with a degree ≥ 20 was used to identify the network as hub genes.
RNA extraction, cDNA synthesis, and RT-qPCR analysis
Total RNA was extracted from drought and salt treated leaf and root tissues by using a Tiangen (www.tiangen.com) RNA extraction kit. RNA was subsequently reverse transcribed into cDNA using a TaKaRa reverse transcription kit. DNaseI treatment was used to eliminate genomic DNA contamination, following manufacturer procedures. Three biological repeats were maintained throughout the analysis . Real-time quantitative PCR was used to assess the transcript profiling of NAC in leaf and root samples under various stress conditions.The following settings were used to run the PCR: 95 °C for 10 min, which was followed by 45 cycles of 95 °C for 30 s and 60 °C for 10 s. The 2^^CT method of analysis was used and each sample was analyzed in triplicate . AtActin was used as the internal reference. The RT-qPCR analysis were carried out three times independently, with three technical repeats for each treatment time point for each experiment .
Plant transformation and screening
Overexpression was performed in Arabidopsis thaliana (Colombia-0). The seedlings were grown in a growth chamber with a temperature of 22 °C, a photoperiod of 16 h of light/8 h of darkness, and an 80% relative humidity. T0-T3 seeds were screened, and T3 lines were used for treatment testing. First generation seeds (T1) were collected from T0 generation Arabidopsis seedlings. T1 seeds were sown on media containing Antibiotics and select the segregant lines. Grow the segregant lines and harvest seeds from them for T3 generation. At this stage, we chose lines with 100% purity and do RT-qPCR and PCR for the final selection of the three best lines . From the seven transformed GhNAC072 lines, we chose OE-1, OE-3, and OE-5. Homozygous lines were used for stress treatment and phenotypic studies.
Determination of drought and salt tolerance in the transgenic lines
GhNAC072 overexpressed Arabidopsis lines, OE-1, OE-3, OE-5, and wildtype seeds were sterilized by soaking them in 15% Sodium hypochlorite for five minutes and then rinsing them three times in sterilized deionized water. The seeds were sterilized before being sown on ½ MS media and kept in a dark chamber at 4 °C for two days before being incubated at 22 °C with a 16 h light/8 h dark photoperiod. After 7 days active seedlings were transplanted into pots that have a 1:1 mix of vermiculite and humus. Three weeks later, the plants were exposed to abiotic stresses. Overexpressed lines OE-1, OE-3, OE-5, and WT plants were treated with 20% polyethylene glycol for drought stress treatment, while 200 mM NaCl was used for salt stress treatment. After eight days, physiological and phenotypic traits were measured. The samples were promptly freeze in liquid nitrogen and kept at a temperature of − 80 °C to be used for subsequent gene expression analysis. All sampling were replicated three times biologically and three times technically .
Germination rate and root elongation determination
Under drought and salt stress conditions, the germination rates of overexpressed lines and wild-types were evaluated. Overexpressed lines OE-1, OE-3, OE-5, and wild-type seedlings were sterilized and spread on plates at ½ MS media with mannitol supplementations of 0 mM, 100 mM, 200 mM, and 300 mM for drought simulation. Similarly, to apply salt stress, ½ MS was supplemented with 0, 100, 150, and 200 mM NaCl levels . The germination percentage was measured after 10 days. To simulate drought and salt stress, transgenic and wildtype seedlings were seeded on ½ MS media for 6 days before being transferred to ½ MS supplemented with different amounts of mannitol and NaCl levels. The seedlings were grown for 6 days before having their roots measured on the 7th day after treatment. All measurements were replicated three times biologically and three times technically .
Stress responsive genes expression analysis
The expression levels of four abiotic stress-responsive genes were explored for GhNAC072 by PCR using specific gene primers to further explore the function of GhNAC072 under conditions of exposure to drought and salinity stress in cotton. The four stress responsive genes include ABF4, SOS1, RAB18, and RD22 . Specific primers for these genes were designed using NCBI (https://www.ncbi.nlm.nih.gov/) database (Table S2).
The data were analyzed by one way analysis of variance (ANOVA), and significant differences among the individual means were determined using Least Significant Difference (P < 0.05 and P < 0.01) at the 5% and 1% confidence levels, respectively. In addition, the graphs were drawn using a PrismPad graph, and the gene family analysis was analyzed and organized using various web tools.
Availability of data and materials
The RNA-Seq data of the three G. hirsutum lines, Mariegalante-85, Upland cotton and Latifolium-40 under drought stress were found in NCBI (Accession: PRJNA663204). The Arabidopsis thaliana Col-0 planting materials used in the experiment came were obtained our laboratory. In addition, additional tables and figures were included as article supplementary files.
Cotton Functional Genomic Database
Least Significant Difference
- H2O2 :
Hydrogen per oxide
Han D, Du M, Zhou Z, Wang S, Li T, Han J, Xu T, Yang G. Overexpression of a Malus baccata NAC Transcription Factor Gene MbNAC25 Increases Cold and Salinity Tolerance in Arabidopsis. Int J Mol Sci. 2020;21(4):1198. https://doi.org/10.3390/ijms21041198.
Duan A, Yang XL, Feng K, Liu JX, Xu ZS, Xiong AS. Genome-wide analysis of NAC transcription factors and their response to abiotic stress in celery (Apium graveolens L.). Comp Biol Chem. 2020;84:107186–107186. https://doi.org/10.1016/j.compbiolchem.2019.107186.
Lambert SA, Jolma A, Campitelli LF, Das PK, Yin Y, Albu M, Chen X, Taipale J, Hughes TR, Weirauch MT. The Human Transcription Factors. Cell. 2018;172:650–65. https://doi.org/10.1016/j.cell.2018.01.029.
Mahajan S, Tuteja N. Cold, salinity and drought stresses: an overview. Arch Biochem Biophys. 2005;444:139–58. https://doi.org/10.1016/j.abb.2005.10.018.
Sewelam N, Kazan K, Schenk PM. Global Plant Stress Signaling: Reactive Oxygen Species at the Cross-Road. Front Plant Sci. 2016;7:187. https://doi.org/10.3389/fpls.2016.00187.
Baillo HE, Kimotho NR, Zhang Z, Xu P. Transcription Factors Associated with Abiotic and Biotic Stress Tolerance and Their Potential for Crops Improvement. Genes (Basel). 2019;10:771. https://doi.org/10.3390/genes10100771.
Yang C, Huang Y, Lv W, Zhang Y, Bhat JA, Kong J, Xing H, Zhao J, Zhao T. GmNAC8 acts as a positive regulator in soybean drought stress. Plant Sci. 2020;293:110442. https://doi.org/10.1016/j.plantsci.2020.110442.
Huang L, Wang Y, Wang W, Zhao X, Qin Q, Sun F, Hu F, Zhao Y, Li Z, Fu B, et al. Characterization of Transcription Factor Gene OsDRAP1 Conferring Drought Tolerance in Rice. Front Plant Sci. 2018;9:94. https://doi.org/10.3389/fpls.2018.00094.
Liu M, Ma Z, Sun W, Huang L, Wu Q, Tang Z, Bu T, Li C, Chen H. Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum). BMC Genom. 2019;20:1–16. https://doi.org/10.1186/s12864-019-5500-0.
Thirumalaikumar VP, Devkar V, Mehterov N, Ali S, Ozgur R, Turkan I, Mueller-roeber B, Balazadeh S. NAC transcription factor JUNGBRUNNEN1 enhances drought tolerance in tomato. Plant Biol J. 2018;1:354–66. https://doi.org/10.1111/pbi.12776.
Jiang D, Zhou L, Chen W, Ye N, Xia J, Zhuang C. Overexpression of a microRNA-targeted NAC transcription factor improves drought and salt tolerance in Rice via ABA-mediated pathways. Rice. 2019;12:76. https://doi.org/10.1186/s12284-019-0334-6.
Hong Y, Zhang H, Huang L, Li D, Song F. Overexpression of a stress-responsive NAC transcription factor gene ONAC022 improves drought and salt tolerance in rice. Front Plant Sci. 2016;7:1–19. https://doi.org/10.3389/fpls.2016.00004.
Wang G, Zhang S, Ma X, Wang Y, Kong F, Meng Q. A stress-associated NAC transcription factor (SlNAC35) from tomato plays a positive role in biotic and abiotic stresses. Physiol Plant. 2016;158:45–64. https://doi.org/10.1111/ppl.12444.
Yang H, Zhang D, Li H, Dong L, Lan H. Ectopic overexpression of the aldehyde dehydrogenase ALDH21 from Syntrichia caninervis in tobacco confers salt and drought stress tolerance. Plant Physiol Biochem. 2015;95:83–91. https://doi.org/10.1016/j.plaphy.2015.07.001.
Guan H, Liu X, Niu F, Zhao Q, Fan N, Cao D, Meng D, He W, Guo B, Wei Y. OoNAC72, a NAC-type Oxytropis ochrocephala transcription factor, conferring enhanced drought and salt stress tolerance in Arabidopsis. Front Plant Sci. 2019;10:890. https://doi.org/10.3389/fpls.2019.00890.
Jia X, Zeng Z, Lyu Y, Zhao S. Drought-Responsive NAC Transcription Factor RcNAC72 Is Recognized by RcABF4, Interacts with RcDREB2A to Enhance Drought Tolerance in Arabidopsis. Int J Mol Sci. 2022;23:1755. https://doi.org/10.3390/ijms23031755.
Li X, Li X, Li M, Yan Y, Liu X, Li L. Dual function of NAC072 in ABF3-mediated ABA-responsive gene regulation in Arabidopsis. Front Plant Sci. 2016;7:1075. https://doi.org/10.3389/fpls.2016.01075.
Han M, Lu X, Yu J, Chen X, Wang X, Malik WA, Wang J, Wang D, Wang S, Guo L, et al. Transcriptome Analysis Reveals Cotton (Gossypium hirsutum) Genes That Are Differentially Expressed in Cadmium Stress Tolerance. Int J Mol Sci. 2019;20(6):1479. https://doi.org/10.3390/ijms20061479.
Fan Y, Lu X, Chen X, Wang J, Wang D, Wang S, Guo L, Rui C, Zhang Y, Cui R, et al. Cotton transcriptome analysis reveals novel biological pathways that eliminate reactive oxygen species (ROS) under sodium bicarbonate (NaHCO3) alkaline stress. Genomics. 2021;113:1157–69. https://doi.org/10.1016/j.ygeno.2021.02.022.
Firmino AAP, Fonseca FCA, de Macedo LLP, Coelho RR, de Souza JDA, Togawa RC, Silva-Junior OB, Pappas-Jr GJ, Silva MCMD, Engler G. Transcriptome Analysis in Cotton Boll Weevil (Anthonomus grandis) and RNA Interference in Insect Pests. PLoS One. 2013;8(12):e85079. https://doi.org/10.1371/journal.pone.0085079.
Bedre R, Rajasekaran K, Mangu VR, Sanchez Timm LE, Bhatnagar D, Baisakh N. Genome-Wide Transcriptome Analysis of Cotton (Gossypium hirsutum L.) Identifies Candidate Gene Signatures in Response to Aflatoxin Producing Fungus Aspergillus flavus. PLos One. 2015;10(9):e0138025. https://doi.org/10.1371/journal.pone.0138025.
Li PT, Wang M, Lu QW, Ge Q, Rashid MHO, Liu AY, Gong JW, Shang HH, Gong WK, Li JW, et al. Comparative transcriptome analysis of cotton fiber development of Upland cotton (Gossypium hirsutum) and Chromosome Segment Substitution Lines from G. hirsutum x G. barbadense. BMC Genom. 2017;18:705. https://doi.org/10.1186/s12864-017-4077-8.
Cheng XQ, Zhang XY, Xue F, Zhu SH, Li YJ, Zhu QH, Liu F, Sun J. Characterization and transcriptome analysis of a dominant genic male sterile cotton mutant. BMC Plant Biol. 2020;20:312. https://doi.org/10.1186/s12870-020-02522-0.
Wang Q, Lu X, Chen X, Malik WA, Wang D, Zhao L, Wang J, Wang S, Guo L, Cui R, et al. Transcriptome analysis of upland cotton revealed novel pathways to scavenge reactive oxygen species (ROS) responding to Na2SO4 tolerance. Sci Rep. 2021;11:8670. https://doi.org/10.1038/s41598-021-87999-x.
Malik WA, Wang X, Wang X, Shu N, Cui R, Chen X, Wang D, Lu X, Yin Z, Wang J, et al. Genome-wide expression analysis suggests glutaredoxin genes response to various stresses in cotton. Int J Biol Macromol. 2020;153:470–91. https://doi.org/10.1016/j.ijbiomac.2020.03.021.
Li F, Fan G, Lu C, Xiao G, Zou C, Kohel RJ, Ma Z, Shang H, Ma X, et al. Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat Biotechnol. 2015;33:524–30. https://doi.org/10.1038/nbt.3208.
Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9. https://doi.org/10.1038/75556.
Zhu JK. Abiotic Stress Signaling and Responses in Plants. Cell. 2016;167:313–24. https://doi.org/10.1016/j.cell.2016.08.029.
Parihar P, Singh S, Singh R, Singh VP, Prasad SM. Effect of salinity stress on plants and its tolerance strategies: a review. Environ Sci Pollut Res Int. 2015;22:4056–75. https://doi.org/10.1007/s11356-014-3739-1.
Yang Y, Guo Y. Unraveling salt stress signaling in plants. J Integr Plant Biol. 2018;60:796–804. https://doi.org/10.1111/jipb.12689.
Pavesi G, Zambelli F, Caggese C, Pesole G. Exalign: a new method for comparative analysis of exon-intron gene structures. Nucleic Acids Res. 2008;36:e47. https://doi.org/10.1093/nar/gkn153.
He Q, Jones DC, Li W, Xie F, Ma J, Sun R, Wang Q, Zhu S, Zhang B. Genome-Wide Identification of R2R3-MYB Genes and Expression Analyses during Abiotic Stress in Gossypium raimondii. Sci Rep. 2016;6:1–14. https://doi.org/10.1038/srep22980.
Ma J, Wang LY, Dai JX, Wang Y, Lin D. The NAC-type transcription factor CaNAC46 regulates the salt and drought tolerance of transgenic Arabidopsis thaliana. BMC Plant Biol. 2021;21:11. https://doi.org/10.1186/s12870-020-02764-y.
Al-Abdallat AM, Ali-Sheikh-Omar MA, Alnemer LM. Overexpression of two ATNAC3-related genes improves drought and salt tolerance in tomato (Solanum lycopersicum L.). Plant Cell Tissue Organ Cult (PCTOC). 2014;120:989–1001. https://doi.org/10.1007/s11240-014-0652-8.
Dudhate A, Shinde H, Yu P, Tsugama D, Gupta SK, Liu S, Takano T. Comprehensive analysis of NAC transcription factor family uncovers drought and salinity stress response in pearl millet (Pennisetum glaucum). BMC Genom. 2021;22:70. https://doi.org/10.1186/s12864-021-07382-y.
Lu P, Magwanga RO, Kirungu JN, Hu Y, Dong Q, Cai X, Zhou Z, Wang X, Zhang Z, Hou Y, et al. Overexpression of cotton a DTX/MATE gene enhances drought, salt, and cold stress tolerance in transgenic arabidopsis. Front Plant Sci. 2019;10:299. https://doi.org/10.3389/fpls.2019.00299.
Sadau SB, Ahmad A, Tajo SM, Ibrahim S, Kazeem BB, Wei H, Yu S. Overexpression of GhMPK3 from Cotton Enhances Cold, Drought, and Salt Stress in Arabidopsis. Agronomy. 2021;11:1049. https://doi.org/10.3390/agronomy11061049.
Lu P, Magwanga RO, Lu H, Kirungu JN, Wei Y, Dong Q, Wang X, Cai X, Zhou Z, Wang K, et al. A novel G-protein-coupled receptors gene from upland cotton enhances salt stress tolerance in transgenic Arabidopsis. Genes. 2018;9:209. https://doi.org/10.3390/genes9040209.
Hu W, Wei Y, Xia Z, Yan Y, Hou X, Zou M, Lu C, Wang W, Peng M. Genome-wide identification and expression analysis of the NAC transcription factor family in cassava. PLoS ONE. 2015;10:1–25. https://doi.org/10.1371/journal.pone.0136993.
Yang X, Kirungu JN, Magwanga RO, Xu Y, Pu L, et al. Knockdown of GhIQD31 and GhIQD32 increases drought and salt stress sensitivity in Gossypium hirsutum. Plant Physiol Biochem. 2019;144:166–77. https://doi.org/10.1016/j.plaphy.2019.09.027.
Li Z, Li L, Zhou K, Zhang Y, Han X, Din Y, Ge X, Qin W, Wang P, et al. GhWRKY6 Acts as a Negative Regulator in Both Transgenic Arabidopsis and Cotton During Drought and Salt Stress. Front Genet. 2019;10:392. https://doi.org/10.3389/fgene.2019.00392.
Kumar S, Stecher G, Tamura K. MEGA7. molecular evolutionary genetics analysis version 70 for bigger datasets. Mol Biol Evol. 2016;33:1870–4. https://doi.org/10.1093/molbev/msw054.
Hu B, Jin J, Guo A, Zhang H, Luo J, Gao G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics. 2015;31:1296–7. https://doi.org/10.1093/bioinformatics/btu817.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Molecular plant. 2020;13:81194–202. https://doi.org/10.1016/j.molp.2020.06.009.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: A software Environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504. https://doi.org/10.1101/gr.1239303.
Mehari TG, Xu Y, Magwanga RO, Umer MJ, Shiraku ML, Hou Y, et al. Identification and functional characterization of Gh_D01G0514 (GhNAC072) transcription factor in response to drought stress tolerance in cotton. Plant Physiol Biochem. 2021;166:361–75. https://doi.org/10.1016/j.plaphy.2021.05.050.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8. https://doi.org/10.1038/nprot.2008.73.
Nieuwenhuizen NJ, Chen X, Wang MY, Matich AJ, Perez RL, Allan AC, Green SA, Atkinson RG. Natural variation in monoterpene synthesis in kiwifruit: Transcriptional regulation of terpene synthases by NAC and ETHYLENE-INSENSITIVE3-like transcription factors. Plant Physiol. 2015;167:1243–58. https://doi.org/10.1104/pp.114.254367.
Elasad M, Ahmad A, Wang H, Ma L, Yu S, Wei H. Overexpression of CDSP32 (GhTRX134) Cotton Gene Enhances Drought, Salt, and Oxidative Stress Tolerance in Arabidopsis. Plants (Basel). 2020;9:1388. https://doi.org/10.3390/plants9101388.
Li J, Yang Y, Sun K, Chen Y, Chen X, Li X. Exogenous melatonin enhances cold, salt and drought stress tolerance by improving antioxidant defense in tea plant (Camellia sinensis (L.) O. Kuntze). Molecules. 2019;24(9):1826. https://doi.org/10.3390/molecules24091826.
Shiraku ML, Magwanga RO, Cai X, Kirungu JN, Xu Y, Mehari TG, Hou Y, Wang Y, Agong SG, Peng R, et al. Functional Characterization of GhACX3 Gene Reveals Its Significant Role in Enhancing Drought and Salt Stress Tolerance in Cotton. Front Plant Sci. 2021;12:658755. https://doi.org/10.3389/fpls.2021.658755.
Mehari TG, Xu Y, Umer MJ, Shiraku ML, Hou Y, Wang Y, Yu S, Zhang X, Wang K, Cai X, et al. Multi-Omics-Based Identification and Functional Characterization of Gh_A06G1257 Proves Its Potential Role in Drought Stress Tolerance in Gossypium hirsutum. Front Plant Sci. 2021;12:746771. https://doi.org/10.3389/fpls.2021.746771.
We are grateful to all the colleagues in our laboratory that provided technical assistance.
This research was funded by the the National Key R&D Program of China ((2021YFE0101200), PSF/CRP/18thProtocol (07), National Natural Science Foundation of China (32072023, 32171994), Postgraduate Improvement Project of Henan Province (YJS2020JD47).
Ethics approval and consent to participate
The study complied with relevant institutional, national, and international guidelines and legislation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Figure S1. Chromosomal locations of NAC genes in three cotton species. The chromosomal position of cotton Species were mapped inline to their genome. A, G. hirsutum with At subgenome B, G. hirsutum with Dt subgenome C, G. arboreum D, G. raimondii E, Scaffold collections.
Supplementary Figure S2. Gene ontology (GO) annotation classification. A, Biological function for G. hirsutum B, Molecular function of G. hirsutum C, Biological function for G. arboreum D, Molecular function of G. arboreum E, Biological function for G. raimondii F, Molecular function of G. raimondii, there is no significant classification for cellular function, the logo stands for the level of significance.
Supplementary Figure S3. Phylogenetic relationships of gene structure analysis and Motif Identification in Gossypium species of NAC genes. A, G. hirsutum, B, G. arboreum, C, G. raimondii.
Supplementary Figure S4. PCR amplification and gel band formation of the 1038bp coding sequence gene Gh_D01G0514 (GhNAC072) using 5000bp marker.
Supplementary Table S1. Number of NAC genes and their assembly in G. hirsutum, G. arboreum, and G. raimondii.
Supplementary Table S2. List of primers for cloning, stress responsive genes and RT-qPCR analysis.
Supplementary Table S3. Data for germination rate, and root length of transgenic and wildtype during drought and salt stress treatment.
Supplementary Table S4. Physiological traits (ELWL, RLWC, CMS and Chlorophyll content) of transgenic and wildtype during drought and salt stress treatment.
Supplementary Table S5. Biochemical parameters (SOD, CAT, MDA and H2O2) of transgenic and wildtype during drought and salt stress treatment.
Supplementary Table S6. Data for relative expression of transgenic and wildtype during drought and salt stress treatment.
About this article
Cite this article
Mehari, T.G., Hou, Y., Xu, Y. et al. Overexpression of cotton GhNAC072 gene enhances drought and salt stress tolerance in transgenic Arabidopsis. BMC Genomics 23, 648 (2022). https://doi.org/10.1186/s12864-022-08876-z