Genome-wide identification and expression analysis of the CaNAC family members in chickpea during development, dehydration and ABA treatments

Characterization of NAC Gene Family in Ammopiptanthus mongolicus and Functional Analysis of AmNAC24, an Osmotic and Cold-Stress-Induced NAC Gene

The NAC family of transcription factors (TFs) is recognized as a significant group within the plant kingdom, contributing crucially to managing growth and development processes in plants, as well as to their response and adaptation to various environmental stressors. Ammopiptanthus mongolicus, a temperate evergreen shrub renowned for its remarkable resilience to low temperatures and drought stress, presents an ideal subject for investigating the potential involvement of NAC TFs in stress response mechanisms. Here, the structure, evolution, and expression profiles of NAC family TFs were analyzed systematically, and a cold and osmotic stress-induced member, AmNAC24, was selected and functionally characterized. A total of 86 NAC genes were identified in A. mongolicus, and these were divided into 15 groups. Up to 48 and 8 NAC genes were generated by segmental duplication and tandem duplication, respectively, indicating that segmental duplication is a predominant mechanism in the expansion of the NAC gene family in A. mongolicus. A considerable amount of NAC genes, including AmNAC24, exhibited upregulation in response to cold and osmotic stress. This observation is in line with the detection of numerous cis-acting elements linked to abiotic stress response in the promoters of A. mongolicus NAC genes. Subcellular localization revealed the nuclear residence of the AmNAC24 protein, coupled with demonstrable transcriptional activation activity. AmNAC24 overexpression enhanced the tolerance of cold and osmotic stresses in Arabidopsis thaliana, possibly by maintaining ROS homeostasis. The present study provided essential data for understanding the biological functions of NAC TFs in plants.

A Transcription Factor SlNAC4 Gene of Suaeda liaotungensis Enhances Salt and Drought Tolerance through Regulating ABA Synthesis

The NAC (NAM, ATAF1/2 and CUC2) transcription factors are ubiquitously distributed in plants and play critical roles in the construction of plant organs and abiotic stress response. In this study, we described the cloning of a Suaeda liaotungensis K. NAC transcription factor gene SlNAC4, which contained 1450 bp, coding a 331 amino acid. We found that SlNAC4 was highly expressed in stems of S. liaotungensis, and the expression of SlNAC4 was considerably up-regulated after salt, drought, and ABA treatments. Transcription analysis and subcellular localization demonstrated that the SlNAC4 protein was located both in the nucleus and cytoplasm, and contained a C-terminal transcriptional activator. The SlNAC4 overexpression Arabidopsis lines significantly enhanced the tolerance to salt and drought treatment and displayed obviously increased activity of antioxidant enzymes under salt and drought stress. Additionally, transgenic plants overexpressing SlNAC4 had a significantly higher level of physiological indices. Interestingly, SlNAC4 promoted the expression of ABA metabolism-related genes including AtABA1, AtABA3, AtNCED3, AtAAO3, but inhibited the expression of AtCYP707A3 in overexpression lines. Using a yeast one-hybrid (Y1H) assay, we identified that the SlNAC4 transcription factor could bind to the promoters of those ABA metabolism-related genes. These results indicate that overexpression of SlNAC4 in plants enhances the tolerance to salt and drought stress by regulating ABA metabolism.

Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Genome-wide analyses of the mung bean NAC gene family reveals orthologs, co-expression networking and expression profiling under abiotic and biotic stresses

BACKGROUND

Mung bean is a short-duration and essential food crop owing to its cash prominence in Asia. Mung bean seeds are rich in protein, fiber, antioxidants, and phytonutrients. The NAC transcription factors (TFs) family is a large plant-specific family, participating in tissue development regulation and abiotic and biotic stresses.

RESULTS

In this study, we perform genome-wide comparisons of VrNAC with their homologs from Arabidopsis. We identified 81 NAC transcription factors (TFs) in mung bean genome and named as per their chromosome location. A phylogenetic analysis revealed that VrNACs are broadly distributed in nine groups. Moreover, we identified 20 conserved motifs across the VrNACs highlighting their roles in different biological process. Based on the gene structure of the putative VrNAC and segmental duplication events might be playing a vital role in the expansion of mung bean genome. A comparative phylogenetic analysis of mung bean NAC together with homologs from Arabidopsis allowed us to classify NAC genes into 13 groups, each containing several orthologs and paralogs. Gene ontology (GO) analysis categorized the VrNACs into biological process, cellular components and molecular functions, explaining the functions in different plant physiology processes. A gene co-expression network analysis identified 173 genes involved in the transcriptional network of putative VrNAC genes. We also investigated how miRNAs potentially target VrNACs and shape their interactions with proteins. VrNAC1.4 (Vradi01g03390.1) was targeted by the Vra-miR165 family, including 9 miRNAs. Vra-miR165 contributes to leaf development and drought tolerance. We also performed qRT-PCR on 22 randomly selected VrNAC genes to assess their expression patterns in the NM-98 genotype, widely known for being tolerant to drought and bacterial leaf spot disease.

CONCLUSIONS

This genome-wide investigation of VrNACs provides a unique resource for further detailed investigations aimed at predicting orthologs functions and what role the play under abiotic and biotic stress, with the ultimate aim to improve mung bean production under diverse environmental conditions.

Genome-Wide Identification and Analysis of the NAC Transcription Factor Gene Family in Garden Asparagus (Asparagus officinalis)

As a large plant-specific gene family, the NAC (NAM, ATAF1/2, and CUC2) transcription factor is related to plant growth, development, and response to abiotic stresses. Although the draft genome of garden asparagus (Asparagus officinalis) has been released, the genome-wide investigation of the NAC gene family is still unavailable. In this study, a total of 85 A. officinalis NAC genes were identified, and a comprehensive analysis of the gene family was performed, including physicochemical properties, phylogenetic relationship, chromosome localization, gene structure, conserved motifs, intron/exon, cis-acting elements, gene duplication, syntenic analysis, and differential gene expression analysis. The phylogenetic analysis demonstrated that there were 14 subgroups in both A. officinalis and Arabidopsis thaliana, and the genes with a similar gene structure and motif distribution were clustered in the same group. The cis-acting regulatory analysis of AoNAC genes indicated four types of cis-acting elements were present in the promoter regions, including light-responsive, hormone-responsive, plant-growth-and-development-related, and stress-responsive elements. The chromosomal localization analysis found that 81 NAC genes in A. officinalis were unevenly distributed on nine chromosomes, and the gene duplication analysis showed three pairs of tandem duplicated genes and five pairs of segmental duplications, suggesting that gene duplication is possibly associated with the amplification of the A. officinalis NAC gene family. The differential gene expression analysis revealed one and three AoNAC genes that were upregulated and downregulated under different types of salinity stress, respectively. This study provides insight into the evolution, diversity, and characterization of NAC genes in garden asparagus and will be helpful for future understanding of their biological roles and molecular mechanisms in plants.

Insights into the gene and protein structures of the CaSWEET family members in chickpea (Cicer arietinum), and their gene expression patterns in different organs under various stress and abscisic acid treatments

'Sugars Will Eventually be Exported Transporters' (SWEETs) are a group of sugar transporters that play crucial roles in various biological processes, particularly plant stress responses. However, no information is available yet for the CaSWEET family in chickpea. Here, we identified all putative CaSWEET members in chickpea, and obtained their major characteristics, including physicochemical patterns, chromosomal distribution, subcellular localization, gene organization, conserved motifs and three-dimensional protein structures. Subsequently, we explored available transcriptome data to compare spatiotemporal transcript abundance of CaSWEET genes in various major organs. Finally, we studied the changes in their transcript levels in leaves and/or roots following dehydration and exogenous abscisic acid treatments using RT-qPCR to obtain valuable information underlying their potential roles in chickpea responses to water-stress conditions. Our results provide the first insights into the characteristics of the CaSWEET family members and a foundation for further functional characterizations of selected candidate genes for genetic engineering of chickpea.

Genome-wide identification and expression analysis of the GRAS gene family in response to drought stress in chickpea (Cicer arietinum L.)

The GRAS (gibberellic acid insensitive, repressor of GAI and scarecrow) transcription factors (TFs) regulate diverse biological processes involved in plant growth and development. These TFs are also known to regulate gene expression in response to various abiotic stress factors like cold, drought, etc. In chickpea one of the most devastating abiotic stress factors is terminal drought. The GRAS TF family has not been characterized in chickpea (Cicer arietinum L.) until now. In this study, we report 46 GRAS TF genes (CaGRAS genes) in the chickpea genome. The CaGRAS proteins were categorized into nine subfamilies based on their phylogenetic relationship with known GRAS members of Arabidopsis and soybean. The PAT subfamily was the largest consisting of ten CaGRAS members whereas the LAS subfamily was the smallest with only one member. Gene duplication analysis revealed that segmental duplication was the primary reason for the expansion of this gene family within the chickpea genome. The gene expression levels of CaGRAS genes were analysed using two different chickpea varieties contrasting for drought tolerance trait, i.e., ICC 4958 (drought tolerant) and ICC 1882 (drought sensitive). On exposure to drought stress, the two chickpea genotypes, exhibited differential drought response, which was quantified and estimated in terms of differences in leaf relative water content (RWC). The well-watered or control plants of the drought tolerant variety were able to maintain a higher leaf RWC by the end of the drought stress period, whereas the control plants of the drought sensitive variety continued to show a decline in leaf RWC. The two genotypes also differed in their root morphologies, under well-watered and drought stress conditions. The gene expression analysis revealed a potential role of PAT, SCR, SCL3 and SHR GRAS members in the regulation of differential response to drought, in the root tissues, for both the genotypes. CaGRAS 12 (SCR) was identified as a drought-responsive GRAS TF gene, which could serve as a potential candidate gene for utilization in developing chickpea varieties with improved drought tolerance. This study demonstrates the drought-responsive expression of CaGRAS genes in chickpea and also describes the morpho-physiological response of chickpea plants to drought stress conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03104-z.

Identification and functional characterization of Gh_D01G0514 (GhNAC072) transcription factor in response to drought stress tolerance in cotton

Cotton encounters long-term drought stress problems resulting in major yield losses. Transcription factors (TFs) plays an important role in response to biotic and abiotic stresses. The coexpression patterns of gene networks associated with drought stress tolerance were investigated using transcriptome profiles. Applying a weighted gene coexpression network analysis, we discovered a salmon module with 144 genes strongly linked to drought stress tolerance. Based on coexpression and RT-qPCR analysis GH_D01G0514 was selected as the candidate gene, as it was also identified as a hub gene in both roots and leaves with a consistent expression in response to drought stress in both tissues. For validation of GH_D01G0514, Virus Induced Gene Silencing was performed and VIGS plants showed significantly higher excised leaf water loss and ion leakage, while lower relative water and chlorophyll contents as compared to WT (Wild type) and positive control plants. Furthermore, the WT and positive control seedlings showed higher CAT and SOD activities, and lower activities of hydrogen peroxide and MDA enzymes as compared to the VIGS plants. Gh_D01G0514 (GhNAC072) was localized in the nucleus and cytoplasm. Y2H assay demonstrates that Gh_D01G0514 has a potential of auto activation. It was observed that the Gh_D01G0514 was highly upregulated in both tissues based on RNA Seq and RT-qPCR analysis. Thus, we inferred that, this candidate gene might be responsible for drought stress tolerance in cotton. This finding adds significantly to the existing knowledge of drought stress tolerance in cotton and deep molecular analysis are required to understand the molecular mechanisms underlying drought stress tolerance in cotton.

Genome-wide identification and characterization of NAC genes in Brassica juncea var. tumida

Background

NAC (NAM, ATAF1/2, and CUC2) transcription factors play an important role in plant growth and development. However, in tumorous stem mustard (Brassica juncea var. tumida), one of the economically important crops cultivated in southwest China and some southeast Asian countries, reports on the identification of NAC family genes are lacking. In this study, we conducted a genome-wide investigation of the NAC family genes in B. juncea var. tumida, based on its recently published genome sequence data.

Methods

The NAC genes were identified in B. juncea var. tumida using the bioinformatics approach on the whole genome level. Additionally, the expression of BjuNAC genes was analyzed under high- and low-temperature stresses by quantitative real-time PCR (qRT-PCR).

Results

A total of 300 BjuNAC genes were identified, of which 278 were mapped to specific chromosomes. Phylogenetic analysis of B. juncea var. tumida, Brassica rapa, Brassica nigra, rice and Arabidopsis thaliana NAC proteins revealed that all NAC genes were divided into 18 subgroups. Furthermore, gene structure analysis showed that most of the NAC genes contained two or three exons. Conserved motif analysis revealed that BjuNAC genes contain a conserved NAM domain. Additionally, qRT-PCR data indicated that thirteen BjuNAC genes with a varying degree of up-regulation during high-temperature stress. Conversely, four BjuNAC genes (BjuNAC006, BjuNAC083, BjuNAC170 and BjuNAC223) were up-regulated and two BjuNAC genes (BjuNAC074 and BjuNAC295) down-regulated under low temperature, respectively. Together, the results of this study provide a strong foundation for future investigation of the biological function of NAC genes in B. juncea var. tumida.

The biotechnological importance of the plant-specific NAC transcription factor family in crop improvement

Climate change, malnutrition, and food insecurity are the inevitable challenges being faced by the agriculture sector today. Plants are susceptible to extreme temperatures during the crucial phases of flowering and seed development, and elevated carbon levels also lead to yield losses. Productivity is also affected by floods and droughts. Therefore, increasing plant yield and stress tolerance are the priorities to be met through novel biotechnological interventions. The contributions of NAC genes towards enhancing plant survivability under stress is well known. Here we focus on the potential of NAC genes in the regulation of abiotic stress tolerance, secondary cell wall synthesis, lateral root development, yield potential, seed size and biomass, ROS signaling, leaf senescence, and programmed cell death. Once naturally tolerant candidate NAC genes have been identified, and the nature of their association with growth and fitness against multi-environmental stresses has been determined, they can be exploited for building inherent tolerance in future crops via transgenic technologies. An update on the latest developments is provided in this review, which summarizes the current understanding of the roles of NAC in the establishment of various stress-adaptive mechanisms in model and food crop plants.

Comprehensive analysis and identification of drought-responsive candidate NAC genes in three semi-arid tropics (SAT) legume crops

Background

Chickpea, pigeonpea, and groundnut are the primary legume crops of semi-arid tropics (SAT) and their global productivity is severely affected by drought stress. The plant-specific NAC (NAM - no apical meristem, ATAF - Arabidopsis transcription activation factor, and CUC - cup-shaped cotyledon) transcription factor family is known to be involved in majority of abiotic stresses, especially in the drought stress tolerance mechanism. Despite the knowledge available regarding NAC function, not much information is available on NAC genes in SAT legume crops.

Results

In this study, genome-wide NAC proteins – 72, 96, and 166 have been identified from the genomes of chickpea, pigeonpea, and groundnut, respectively, and later grouped into 10 clusters in chickpea and pigeonpea, while 12 clusters in groundnut. Phylogeny with well-known stress-responsive NACs in Arabidopsis thaliana, Oryza sativa (rice), Medicago truncatula, and Glycine max (soybean) enabled prediction of putative stress-responsive NACs in chickpea (22), pigeonpea (31), and groundnut (33). Transcriptome data revealed putative stress-responsive NACs at various developmental stages that showed differential expression patterns in the different tissues studied. Quantitative real-time PCR (qRT-PCR) was performed to validate the expression patterns of selected stress-responsive, Ca_NAC (Cicer arietinum - 14), Cc_NAC (Cajanus cajan - 15), and Ah_NAC (Arachis hypogaea - 14) genes using drought-stressed and well-watered root tissues from two contrasting drought-responsive genotypes of each of the three legumes. Based on expression analysis, Ca_06899, Ca_18090, Ca_22941, Ca_04337, Ca_04069, Ca_04233, Ca_12660, Ca_16379, Ca_16946, and Ca_21186; Cc_26125, Cc_43030, Cc_43785, Cc_43786, Cc_22429, and Cc_22430; Ah_ann1.G1V3KR.2, Ah_ann1.MI72XM.2, Ah_ann1.V0X4SV.1, Ah_ann1.FU1JML.2, and Ah_ann1.8AKD3R.1 were identified as potential drought stress-responsive candidate genes.

Conclusion

As NAC genes are known to play role in several physiological and biological activities, a more comprehensive study on genome-wide identification and expression analyses of the NAC proteins have been carried out in chickpea, pigeonpea and groundnut. We have identified a total of 21 potential drought-responsive NAC genes in these legumes. These genes displayed correlation between gene expression, transcriptional regulation, and better tolerance against drought. The identified candidate genes, after validation, may serve as a useful resource for molecular breeding for drought tolerance in the SAT legume crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07602-5.

Genome-wide identification and expression analysis of the NAC transcription factor family in tomato (Solanum lycopersicum) during aluminum stress

Background

The family of NAC proteins (NAM, ATAF1/2, and CUC2) represent a class of large plant-specific transcription factors. However, identification and functional surveys of NAC genes of tomato (Solanum lycopersicum) remain unstudied, despite the tomato genome being decoded for several years. This study aims to identify the NAC gene family and investigate their potential roles in responding to Al stress.

Results

Ninety-three NAC genes were identified and named in accordance with their chromosome location. Phylogenetic analysis found SlNACs are broadly distributed in 5 groups. Gene expression analysis showed that SlNACs had different expression levels in various tissues and at different fruit development stages. Cycloheximide treatment and qRT-PCR analysis indicated that SlNACs may aid regulation of tomato in response to Al stress, 19 of which were significantly up- or down-regulated in roots of tomato following Al stress.

Conclusion

This work establishes a knowledge base for further studies on biological functions of SlNACs in tomato and will aid in improving agricultural traits of tomato in the future.

The Soybean GmNAC019 Transcription Factor Mediates Drought Tolerance in Arabidopsis in an Abscisic Acid-Dependent Manner

Being master regulators of gene expression, transcription factors (TFs) play important roles in determining plant growth, development and reproduction. To date, many TFs have been shown to positively mediate plant responses to environmental stresses. In the current study, the biological functions of a stress-responsive NAC [NAM (No Apical Meristem), ATAF1/2 (Arabidopsis Transcription Activation Factor1/2), CUC2 (Cup-shaped Cotyledon2)]-TF encoding gene isolated from soybean (GmNAC019) in relation to plant drought tolerance and abscisic acid (ABA) responses were investigated. By using a heterologous transgenic system, we revealed that transgenic Arabidopsis plants constitutively expressing the GmNAC019 gene exhibited higher survival rates in a soil-drying assay, which was associated with lower water loss rate in detached leaves, lower cellular hydrogen peroxide content and stronger antioxidant defense under water-stressed conditions. Additionally, the exogenous treatment of transgenic plants with ABA showed their hypersensitivity to this phytohormone, exhibiting lower rates of seed germination and green cotyledons. Taken together, these findings demonstrated that GmNAC019 functions as a positive regulator of ABA-mediated plant response to drought, and thus, it has potential utility for improving plant tolerance through molecular biotechnology.

Genome-wide identification and expression profile analysis of nuclear factor Y family genes in Sorghum bicolor L. (Moench)

Members of the plant Heme Activator Protein (HAP) or NUCLEAR FACTOR Y (NF-Y) are trimeric transcription factor complexes composed of the NF-YA, NF-YB and NF-YC subfamilies. They bind to the CCAAT box in the promoter regions of the target genes and regulate gene expressions. Plant NF-Ys were reported to be involved in adaptation to several abiotic stresses as well as in development. In silico analysis of Sorghum bicolor genome resulted in the identification of a total of 42 NF-Y genes, among which 8 code for the SbNF-YA, 19 for SbNF-YB and 15 for the SbNF-YC subunits. Analysis was also performed to characterize gene structures, chromosomal distribution, duplication status, protein subcellular localizations, conserved motifs, ancestral protein sequences, miRNAs and phylogenetic tree construction. Phylogenetic relationships and ortholog predictions displayed that sorghum has additional NF-YB genes with unknown functions in comparison with Arabidopsis. Analysis of promoters revealed that they harbour many stress-related cis-elements like ABRE and HSE, but surprisingly, DRE and MYB elements were not detected in any of the subfamilies. SbNF-YA1, 2, and 6 were found upregulated under 200 mM salt and 200 mM mannitol stresses. While NF-YA7 appeared associated with high temperature (40°C) stress, NF-YA8 was triggered by both cold (4°C) and high temperature stresses. Among NF-YB genes, 7, 12, 15, and 16 were induced under multiple stress conditions such as salt, mannitol, ABA, cold and high temperatures. Likewise, NF-YC 6, 11, 12, 14, and 15 were enhanced significantly in a tissue specific manner under multiple abiotic stress conditions. Majority of the mannitol (drought)-inducible genes were also induced by salt, high temperature stresses and ABA. Few of the high temperature stress-induced genes are also induced by cold stress (NF-YA2, 4, 6, 8, NF-YB2, 7, 10, 11, 12, 14, 16, 17, NF-YC4, 6, 12, and 13) thus suggesting a cross talk among them. This work paves the way for investigating the roles of diverse sorghum NF-Y proteins during abiotic stress responses and provides an insight into the evolution of diverse NF-Y members.

A comprehensive analysis of the B3 superfamily identifies tissue-specific and stress-responsive genes in chickpea (Cicer arietinum L.)

The aim of this study was to provide a comprehensive analysis of the plant-specific B3 domain-containing transcription factors (TFs) in chickpea. Scanning of the chickpea genome resulted in the identification of 51 B3 domain-containing TFs that were located on seven out of eight chickpea chromosomes. Based on the presence of additional domains other than the B3 domain, the candidates were classified into four subfamilies, i.e., ARF (24), REM (19), LAV (6) and RAV (2). Phylogenetic analysis classified them into four groups in which members of the same group had similar intron-exon organization and motif composition. Genome duplication analysis of the candidate B3 genes revealed an event of segmental duplication that was instrumental in the expansion of the B3 gene family. Ka/Ks analysis showed that the B3 gene family was under purifying selection. Further, chickpea B3 genes showed maximum orthology with Medicago followed by soybean and Arabidopsis. Promoter analyses of the B3 genes led to the identification of several tissue-specific and stress-responsive cis-regulatory elements. Expression profiling of the candidate B3 genes using publicly available RNA-seq data of several chickpea tissues indicated their putative role in plant development and abiotic stress response. These findings were further validated by real-time expression analysis. Overall, this study provides a comprehensive analysis of the B3 domain-containing proteins in chickpea that would aid in devising strategies for crop manipulation in chickpea.

Systematic analysis of NAC transcription factors' gene family and identification of post-flowering drought stress responsive members in sorghum

SbNAC genes (131) encoding 183 proteins were identified from the sorghum genome and characterized. The expression patterns of SbSNACs were evaluated at three sampling time points under post-flowering drought stress. NAC proteins are specific transcription factors in plants, playing vital roles in development and response to various environmental stresses. Despite the fact that Sorghum bicolor is well-known for its drought-tolerance, it suffers from grain yield loss due to pre and post-flowering drought stress. In the present study, 131 SbNAC genes encoding 183 proteins were identified from the sorghum genome. The phylogenetic trees were constructed based on the NAC domains of sorghum, and also based on sorghum with Arabidopsis and 8 known NAC domains of other plants, which classified the family into 15 and 19 subfamilies, respectively. Based on the obtained results, 13 SbNAC proteins joined the SNAC subfamily, and these proteins are expected to be involved in response to abiotic stresses. Promoter analysis revealed that all SbNAC genes comprise different stress-associated cis-elements in their promoters. UTRs analysis indicated that 101 SbNAC transcripts had upstream open reading frames, while 39 of the transcripts had internal ribosome entry sites in their 5'UTR. Moreover, 298 miRNA target sites were predicted to exist in the UTRs of SbNAC transcripts. The expression patterns of SbSNACs were evaluated in three genotypes at three sampling time points under post-flowering drought stress. Based on the results, it could be suggested that some gene members are involved in response to drought stress at the post-flowering stage since they act as positive or negative transcriptional regulators. Following further functional analyses, some of these genes might be perceived to be promising candidates for breeding programs to enhance drought tolerance in crops.

Identification, Structural Characterization and Gene Expression Analysis of Members of the Nuclear Factor-Y Family in Chickpea (Cicer arietinum L.) under Dehydration and Abscisic Acid Treatments

In plants, the Nuclear Factor-Y (NF-Y) transcription factors (TFs), which include three distinct types of NF-YA, NF-YB, and NF-YC TFs, have been identified to play key roles in the regulation of various plant growth and developmental processes under both normal and environmental stress conditions. In this work, a total of 40 CaNF-Y-encoding genes, including eight CaNF-YAs, 21 CaNF-YBs, and 11 CaNF-YCs, were identified in chickpea, and their major gene and protein characteristics were subsequently obtained using various web-based tools. Of our interest, a phylogenetically-based analysis predicted 18 CaNF-Ys (eight CaNF-YAs, seven CaNF-YBs, and three CaNF-YCs) that potentially play roles in chickpea responses to dehydration according to their close relationship with the well-characterized GmNF-Ys in soybean. These results were in good agreement with the enrichment of drought-responsive cis-regulatory motifs and expression patterns obtained from in silico analyses using publically available transcriptome data. Most of the phylogenetically predicted drought-responsive CaNF-Y genes (15 of 18) were quantitatively validated to significantly respond to dehydration treatment in leaves and/or roots, further supporting the results of in silico analyses. Among these CaNF-Y genes, the transcript levels of CaNF-YA01 and CaNF-YC10 were the most highly accumulated in leaves (by approximately eight-fold) and roots (by approximately 18-fold), respectively, by dehydration. Furthermore, 12 of the 18 CaNF-Y genes were found to be responsive to the most well-known stress hormone, namely abscisic acid (ABA), in leaves and/or roots, suggesting that these genes may act in chickpea response to dehydration in ABA-dependent manner. Taken together, our study has provided a comprehensive and fundamental information for further functional analyses of selected CaNF-Y candidate genes, ultimately leading to the improvement of chickpea growth under water-limited conditions.

Current advances in chickpea genomics: applications and future perspectives

Chickpea genomics promises to illuminate our understanding of genome organization, structural variations, evolutionary and domestication-related insights and fundamental biology of legume crops. Unprecedented advancements of next generation sequencing (NGS) technologies have enabled in decoding of multiple chickpea genome sequences and generating huge genomic resources in chickpea both at functional and structural level. This review is aimed to update the current progress of chickpea genomics ranging from high density linkage map development, genome-wide association studies (GWAS), functional genomics resources for various traits, emerging role of abiotic stress responsive coding and non-coding RNAs after the completion of draft chickpea genome sequences. Additionally, the current efforts of whole genome re-sequencing (WGRS) approach of global chickpea germplasm to capture the global genetic diversity existing in the historically released varieties across the world and increasing the resolution of the previously identified candidate gene(s) of breeding importance have been discussed. Thus, the outcomes of these genomics resources will assist in genomics-assisted selection and facilitate breeding of climate-resilient chickpea cultivars for sustainable agriculture.

Genome-wide identification and expression profile analysis of the NAC transcription factor family during abiotic and biotic stress in woodland strawberry

The NAC transcription factors involved plant development and response to various stress stimuli. However, little information is available concerning the NAC family in the woodland strawberry. Herein, 37 NAC genes were identified from the woodland strawberry genome and were classified into 13 groups based on phylogenetic analysis. And further analyses of gene structure and conserved motifs showed closer relationship of them in every subgroup. Quantitative real-time PCR evaluation different tissues revealed distinct spatial expression profiles of the FvNAC genes. The comprehensive expression of FvNAC genes revealed under abiotic stress (cold, heat, drought, salt), signal molecule treatments (H2O2, ABA, melatonin, rapamycin), biotic stress (Colletotrichum gloeosporioides and Ralstonia solanacearum). Expression profiles derived from quantitative real-time PCR suggested that 5 FvNAC genes responded dramatically to the various abiotic and biotic stresses, indicating their contribution to abiotic and biotic stresses resistance in woodland strawberry. Interestingly, FvNAC genes showed greater extent responded to the cold treatment than other abiotic stress, and H2O2 exhibited a greater response than ABA, melatonin, and rapamycin. For biotic stresses, 3 FvNAC genes were up-regulated during infection with C. gloeosporioides, while 6 FvNAC genes were down-regulated during infection with R. solanacearum. In conclusion, this study identified candidate FvNAC genes to be used for the genetic improvement of abiotic and biotic stress tolerance in woodland strawberry.

Transcription Factors and Their Roles in Signal Transduction in Plants under Abiotic Stresses

In agricultural production, abiotic stresses are known as the main disturbance leading to negative impacts on crop performance. Research on elucidating plant defense mechanisms against the stresses at molecular level has been addressed for years in order to identify the major contributors in boosting the plant tolerance ability. From literature, numerous genes from different species, and from both functional and regulatory gene categories, have been suggested to be on the list of potential candidates for genetic engineering. Noticeably, enhancement of plant stress tolerance by manipulating expression of Transcription Factors (TFs) encoding genes has emerged as a popular approach since most of them are early stress-responsive genes and control the expression of a set of downstream target genes. Consequently, there is a higher chance to generate novel cultivars with better tolerance to either single or multiple stresses. Perhaps, the difficult task when deploying this approach is selecting appropriate gene(s) for manipulation. In this review, on the basis of the current findings from molecular and post-genomic studies, our interest is to highlight the current understanding of the roles of TFs in signal transduction and mediating plant responses towards abiotic stressors. Furthermore, interactions among TFs within the stress-responsive network will be discussed. The last section will be reserved for discussing the potential applications of TFs for stress tolerance improvement in plants.

Transcription Factors and Their Roles in Signal Transduction in Plants under Abiotic Stresses

In agricultural production, abiotic stresses are known as the main disturbance leading to negative impacts on crop performance. Research on elucidating plant defense mechanisms against the stresses at molecular level has been addressed for years in order to identify the major contributors in boosting the plant tolerance ability. From literature, numerous genes from different species, and from both functional and regulatory gene categories, have been suggested to be on the list of potential candidates for genetic engineering. Noticeably, enhancement of plant stress tolerance by manipulating expression of Transcription Factors (TFs) encoding genes has emerged as a popular approach since most of them are early stress-responsive genes and control the expression of a set of downstream target genes. Consequently, there is a higher chance to generate novel cultivars with better tolerance to either single or multiple stresses. Perhaps, the difficult task when deploying this approach is selecting appropriate gene(s) for manipulation. In this review, on the basis of the current findings from molecular and post-genomic studies, our interest is to highlight the current understanding of the roles of TFs in signal transduction and mediating plant responses towards abiotic stressors. Furthermore, interactions among TFs within the stress-responsive network will be discussed. The last section will be reserved for discussing the potential applications of TFs for stress tolerance improvement in plants.

Genetic Analysis of NBS-LRR Gene Family in Chickpea and Their Expression Profiles in Response to Ascochyta Blight Infection

Ascochyta blight is one of the major diseases of chickpea worldwide. The genetic resistance to ascochyta blight in chickpea is complex and governed by multiple QTLs. However, the molecular mechanism of quantitative disease resistance to ascochyta blight and the genes underlying these QTLs are still unknown. Most often disease resistance is determined by resistance (R) genes. The most predominant R-genes contain nucleotide binding site and leucine rich repeat (NBS-LRR) domains. A total of 121 NBS-LRR genes were identified in the chickpea genome. Ninety-eight of these genes contained all essential conserved domains while 23 genes were truncated. The NBS-LRR genes were grouped into eight distinct classes based on their domain architecture. Phylogenetic analysis grouped these genes into two major clusters based on their structural variation, the first cluster with toll or interleukin-1 like receptor (TIR) domain and the second cluster either with or without a coiled-coil domain. The NBS-LRR genes are distributed unevenly across the eight chickpea chromosomes and nearly 50% of the genes are present in clusters. Thirty of the NBS-LRR genes were co-localized with nine of the previously reported ascochyta blight QTLs and were tested as potential candidate genes for ascochyta blight resistance. Expression pattern of these genes was studied in two resistant (CDC Corinne and CDC Luna) and one susceptible (ICCV 96029) genotypes at different time points after ascochyta blight infection using real-time quantitative PCR. Twenty-seven NBS-LRR genes showed differential expression in response to ascochyta blight infection in at least one genotype at one time point. Among these 27 genes, the majority of the NBS-LRR genes showed differential expression after inoculation in both resistant and susceptible genotypes which indicates the involvement of these genes in response to ascochyta blight infection. Five NBS-LRR genes showed genotype specific expression. Our study provides a new insight of NBS-LRR gene family in chickpea and the potential involvement of NBS-LRR genes in response to ascochyta blight infection.

Transcriptome-Wide Identification and Expression Analysis of the NAC Gene Family in Tea Plant [Camellia sinensis (L.) O. Kuntze]

In plants, the NAC (NAM-ATAF1/2-CUC) family of proteins constitutes several transcription factors and plays vital roles in diverse biological processes, such as growth, development, and adaption to adverse factors. Tea, as a non-alcoholic drink, is known for its bioactive ingredients and health efficacy. Currently, knowledge about NAC gene family in tea plant remains very limited. In this study, a total of 45 CsNAC genes encoding NAC proteins including three membrane-bound members were identified in tea plant through transcriptome analysis. CsNAC factors and Arabidopsis counterparts were clustered into 17 subgroups after phylogenetic analysis. Conserved motif analysis revealed that CsNAC proteins with a close evolutionary relationship possessed uniform or similar motif compositions. The distribution of NAC family MTFs (membrane-associated transcription factors) among higher plants of whose genome-wide has been completed revealed that the existence of doubled TMs (transmembrane motifs) may be specific to fabids. Transcriptome analysis exhibited the expression profiles of CsNAC genes in different tea plant cultivars under non-stress conditions. Nine CsNAC genes, including the predicted stress-related and membrane-bound genes, were examined through qRT-PCR (quantitative real time polymerase chain reaction) in two tea plant cultivars, namely, 'Huangjinya' and 'Yingshuang'. The expression patterns of these genes were investigated in different tissues (root, stem, mature leaf, young leaf and bud) and under diverse environmental stresses (drought, salt, heat, cold and abscisic acid). Several CsNAC genes, including CsNAC17 and CsNAC30 that are highly orthologous to known stress-responsive ANAC072/RD26 were identified as highly responsive to abiotic stress. This study provides a global survey of tea plant NAC proteins, and would be helpful for the improvement of stress resistance in tea plant via genetic engineering.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

EcoTILLING-Based Association Mapping Efficiently Delineates Functionally Relevant Natural Allelic Variants of Candidate Genes Governing Agronomic Traits in Chickpea

The large-scale mining and high-throughput genotyping of novel gene-based allelic variants in natural mapping population are essential for association mapping to identify functionally relevant molecular tags governing useful agronomic traits in chickpea. The present study employs an alternative time-saving, non-laborious and economical pool-based EcoTILLING approach coupled with agarose gel detection assay to discover 1133 novel SNP allelic variants from diverse coding and regulatory sequence components of 1133 transcription factor (TF) genes by genotyping in 192 diverse desi and kabuli chickpea accessions constituting a seed weight association panel. Integrating these SNP genotyping data with seed weight field phenotypic information of 192 structured association panel identified eight SNP alleles in the eight TF genes regulating seed weight of chickpea. The associated individual and combination of all SNPs explained 10-15 and 31% phenotypic variation for seed weight, respectively. The EcoTILLING-based large-scale allele mining and genotyping strategy implemented for association mapping is found much effective for a diploid genome crop species like chickpea with narrow genetic base and low genetic polymorphism. This optimized approach thus can be deployed for various genomics-assisted breeding applications with optimal expense of resources in domesticated chickpea. The seed weight-associated natural allelic variants and candidate TF genes delineated have potential to accelerate marker-assisted genetic improvement of chickpea.

CarNAC4, a NAC-type chickpea transcription factor conferring enhanced drought and salt stress tolerances in Arabidopsis

CarNAC4 is a typical stress-responsive NAC transcription factor and enhances drought and salt stress tolerances in transgenic Arabidopsis. Chickpea (Cicer arietinum L.) is relatively vulnerable to abiotic stress conditions, but the tolerance mechanisms for such stresses in chickpea are largely unknown. To identify stress-related factors in chickpea, we previously constructed a cDNA library of chickpea leaves exposed to drought stress conditions. A cDNA encoding a putative NAC transcription factor (CarNAC4) was identified as a putative stress-responsive gene. Our study indicated that the transcript levels of CarNAC4 were enhanced in response to several abiotic stresses and phytohormones. Promoter analysis demonstrated that multiple stress-related cis-acting elements exist in promoter region of CarNAC4. CarNAC4 is localized in the nucleus and binds to the DNA sequence containing CGT[G/A], while the C-terminal region of CarNAC4 contains a transcriptional activation domain. Over-expression of CarNAC4 in Arabidopsis plants improved tolerance to drought and salt stresses. Transgenic plants exhibited greater reduced rates of water loss and more proline accumulation than Col-0 plants under drought stress and less MDA contents than Col-0 plants under salt stress. In addition, over-expression of CarNAC4 enhanced the expression of stress-responsive genes such as RD29A, ERD10, COR15A, COR47, KIN1 and DREB2A. These results indicated that CarNAC4 functions as a transcription factor involved in the regulation of drought and salt stress response.

NAC transcription factors in plant multiple abiotic stress responses: progress and prospects

Abiotic stresses adversely affect plant growth and agricultural productivity. According to the current climate prediction models, crop plants will face a greater number of environmental stresses, which are likely to occur simultaneously in the future. So it is very urgent to breed broad-spectrum tolerant crops in order to meet an increasing demand for food productivity due to global population increase. As one of the largest families of transcription factors (TFs) in plants, NAC TFs play vital roles in regulating plant growth and development processes including abiotic stress responses. Lots of studies indicated that many stress-responsive NAC TFs had been used to improve stress tolerance in crop plants by genetic engineering. In this review, the recent progress in NAC TFs was summarized, and the potential utilization of NAC TFs in breeding abiotic stress tolerant transgenic crops was also be discussed. In view of the complexity of field conditions and the specificity in multiple stress responses, we suggest that the NAC TFs commonly induced by multiple stresses should be promising candidates to produce plants with enhanced multiple stress tolerance. Furthermore, the field evaluation of transgenic crops harboring NAC genes, as well as the suitable promoters for minimizing the negative effects caused by over-expressing some NAC genes, should be considered.

Genome Wide Identification of LIM Genes in Cicer arietinum and Response of Ca-2LIMs in Development, Hormone and Pathogenic Stress

The eukaryotic lineage-specific LIM protein (LIN11, ISL1, and MEC3) family play pivotal role in modulation of actin dynamics and transcriptional regulation. The systematic investigation of this family has not been carried in detail and rare in legumes. Current study involves the mining of Cicer arietinum genome for the genes coding for LIM domain proteins and displayed significant homology with LIM genes of other species. The analysis led to the identification of 15 members, which were positioned on chickpea chromosomes. The phylogenetic and motif analysis suggested their categorization into two sub-families i.e., Ca-2LIMs and Ca-DA1/DAR, which comprised of nine and six candidates, respectively. Further sub-categories of Ca-2LIMs were recognised as αLIM, βLIM, δLIM and γLIM. The LIM genes within their sub-families displayed conserved genomic and motif organization. The expression pattern of Ca-2LIMs across developmental and reproductive tissues demonstrated strong correlation with established consensus. The Ca-2LIM belongs to PLIM and GLIM (XLIM) was found highly expressed in floral tissue. Others showed ubiquitous expression pattern with their dominance in stem. Under hormonal and pathogenic conditions these LIMs were found to up-regulate during salicylic acid, abscisic acid and Ascochyta rabiei treatment or infection; and down-regulated in response to jasmonic acid treatment. The findings of this work, particularly in terms of modulation of LIM genes under biotic stress will open up the way to further explore and establish the role of chickpea LIMs in plant defense response.

CTDB: An Integrated Chickpea Transcriptome Database for Functional and Applied Genomics

Chickpea is an important grain legume used as a rich source of protein in human diet. The narrow genetic diversity and limited availability of genomic resources are the major constraints in implementing breeding strategies and biotechnological interventions for genetic enhancement of chickpea. We developed an integrated Chickpea Transcriptome Database (CTDB), which provides the comprehensive web interface for visualization and easy retrieval of transcriptome data in chickpea. The database features many tools for similarity search, functional annotation (putative function, PFAM domain and gene ontology) search and comparative gene expression analysis. The current release of CTDB (v2.0) hosts transcriptome datasets with high quality functional annotation from cultivated (desi and kabuli types) and wild chickpea. A catalog of transcription factor families and their expression profiles in chickpea are available in the database. The gene expression data have been integrated to study the expression profiles of chickpea transcripts in major tissues/organs and various stages of flower development. The utilities, such as similarity search, ortholog identification and comparative gene expression have also been implemented in the database to facilitate comparative genomic studies among different legumes and Arabidopsis. Furthermore, the CTDB represents a resource for the discovery of functional molecular markers (microsatellites and single nucleotide polymorphisms) between different chickpea types. We anticipate that integrated information content of this database will accelerate the functional and applied genomic research for improvement of chickpea. The CTDB web service is freely available at http://nipgr.res.in/ctdb.html.

Genome-Wide Identification and Expression Analysis of the NAC Transcription Factor Family in Cassava

NAC [no apical meristem (NAM), Arabidopsis transcription activation factor [ATAF1/2] and cup-shaped cotyledon (CUC2)] proteins is one of the largest groups of plant specific transcription factors and plays a crucial role in plant growth, development, and adaption to the environment. Currently, no information is known about the NAC family in cassava. In this study, 96 NAC genes (MeNACs) were identified from the cassava genome. Phylogenetic analysis of the NACs from cassava and Arabidopsis showed that MeNAC proteins can be clustered into 16 subgroups. Gene structure analysis found that the number of introns of MeNAC genes varied from 0 to 5, with the majority of MeNAC genes containing two introns, indicating a small gene structure diversity of cassava NAC genes. Conserved motif analysis revealed that all of the identified MeNACs had the conserved NAC domain and/or NAM domain. Global expression analysis suggested that MeNAC genes exhibited different expression profiles in different tissues between wild subspecies and cultivated varieties, indicating their involvement in the functional diversity of different accessions. Transcriptome analysis demonstrated that MeNACs had a widely transcriptional response to drought stress and that they had differential expression profiles in different accessions, implying their contribution to drought stress resistance in cassava. Finally, the expression of twelve MeNAC genes was analyzed under osmotic, salt, cold, ABA, and H₂O₂ treatments, indicating that cassava NACs may represent convergence points of different signaling pathways. Taken together, this work found some excellent tissue-specific and abiotic stress-responsive candidate MeNAC genes, which would provide a solid foundation for functional investigation of the NAC family, crop improvement and improved understanding of signal transduction in plants. These data bring new insight on the complexity of the transcriptional control of MeNAC genes and support the hypothesis that NACs play an important role in plant growth, development, and adaption of environment.

Genome-wide identification and expression pattern of drought-responsive members of the NAC family in maize

NAC proteins are plant-specific transcription factors (TFs). Although they play a pivotal role in regulating distinct biological processes, TFs in maize are yet to be investigated comprehensively. Within the maize genome, we identified 152 putative NAC domain-encoding genes (ZmNACs), including eight membrane-bound members, by systematic sequence analysis and physically mapped them onto ten chromosomes of maize. In silico analysis of the ZmNACs and comparison with similar genes in other plants such as Arabidopsis, rice, and soybean, revealed a similar NAC sequence architecture. Phylogenetically, the ZmNACs were arranged into six distinct subgroups (I-VI) possessing conserved motifs. Phylogenetic analysis using stress-related NAC TFs from Arabidopsis, rice, and soybean as seeding sequences identified 24 of the 152 ZmNACs (all from Group II) as putative stress-responsive genes, including one dehydration-responsive ZmSNAC1 gene reported earlier. One drought-tolerant genotype (HKI577) and one susceptible genotype (PC13T-3) were used for studying the expression pattern of the NAC genes during drought stress. qRT-PCR based expression profiles of 11 genes predicted to be related to stress confirmed strong differential gene expression during drought stress. Phylogenetic analyses revealed that ZmNAC18, ZmNAC51, ZmNAC145, and ZmNAC72, which were up-regulated in the tolerant genotype and down-regulated in the susceptible genotype, belonged to the same group to which also belong other drought-responsive genes, namely SNAC1, OsNAC6, ANAC019, and ANAC055, which act as a transcriptional activator and are strongly induced under stress from various abiotic sources. Differentially expressed ZmNAC genes, alone or in combination with each other or with other type(s) of TFs, may control the general cellular machinery and regulate stress-responsive downstream genes. Alternatively, they may serve as a platform to regulate a broad set of genes, which are subsequently fine-tuned by specific regulators. This genome-wide identification and expression profiling opens new avenues for systematic functional analysis of new members of the NAC gene family, which may be exploited in developing lines that are better adapted to drought.