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

Transcriptional Control of Seed Life: New Insights into the Role of the NAC Family

Transcription factors (TFs) regulate gene expression by binding to specific sequences on DNA through their DNA-binding domain (DBD), a universal process. This update conveys information about the diverse roles of TFs, focusing on the NACs (NAM-ATAF-CUC), in regulating target-gene expression and influencing various aspects of plant biology. NAC TFs appeared before the emergence of land plants. The NAC family constitutes a diverse group of plant-specific TFs found in mosses, conifers, monocots, and eudicots. This update discusses the evolutionary origins of plant NAC genes/proteins from green algae to their crucial roles in plant development and stress response across various plant species. From mosses and lycophytes to various angiosperms, the number of NAC proteins increases significantly, suggesting a gradual evolution from basal streptophytic green algae. NAC TFs play a critical role in enhancing abiotic stress tolerance, with their function conserved in angiosperms. Furthermore, the modular organization of NACs, their dimeric function, and their localization within cellular compartments contribute to their functional versatility and complexity. While most NAC TFs are nuclear-localized and active, a subset is found in other cellular compartments, indicating inactive forms until specific cues trigger their translocation to the nucleus. Additionally, it highlights their involvement in endoplasmic reticulum (ER) stress-induced programmed cell death (PCD) by activating the vacuolar processing enzyme (VPE) gene. Moreover, this update provides a comprehensive overview of the diverse roles of NAC TFs in plants, including their participation in ER stress responses, leaf senescence (LS), and growth and development. Notably, NACs exhibit correlations with various phytohormones (i.e., ABA, GAs, CK, IAA, JA, and SA), and several NAC genes are inducible by them, influencing a broad spectrum of biological processes. The study of the spatiotemporal expression patterns provides insights into when and where specific NAC genes are active, shedding light on their metabolic contributions. Likewise, this review emphasizes the significance of NAC TFs in transcriptional modules, seed reserve accumulation, and regulation of seed dormancy and germination. Overall, it effectively communicates the intricate and essential functions of NAC TFs in plant biology. Finally, from an evolutionary standpoint, a phylogenetic analysis suggests that it is highly probable that the WRKY family is evolutionarily older than the NAC family.

Genome-wide identification of NAC transcription factors and regulation of monoterpenoid indole alkaloid biosynthesis in Catharanthus roseus

NAC transcription factors (TFs) are crucial to growth and defense responses in plants. Though NACs have been characterized for their role in several plants, comprehensive information regarding their role in Catharanthus roseus, a perennial ornamental plant, is lacking. Homology modelling was employed to identify and characterize NACs in C. roseus. In-vitro propagation of C. roseus plants was carried out using cell suspension and nodal culture and were elicited with two auxin-antagonists, 5-fluoro Indole Acetic Acid (5-F-IAA) and α-(phenyl ethyl-2-oxo)-Indole-Acetic-Acid (PEO-IAA) for the enhanced production of monoterpenoid indole alkaloids (MIAs) namely catharanthine, vindoline, and vinblastine. Analyses revealed the presence of 47 putative CrNAC genes in the C. roseus genome, primarily localized in the nucleus. Phylogenetic analysis categorized these CrNACs into eight clusters, demonstrating the highest synteny with corresponding genes in Camptotheca acuminata. Additionally, at least one defense or hormone-responsive cis-acting element was identified in the promoter region of all the putative CrNACs. Of the two elicitors, 5-F-IAA was effective at 200 µM to elicit a 3.07-fold increase in catharanthine, 2.76-fold in vindoline, and 2.4-fold in vinblastine production in nodal culture. While a relatively lower increase in MIAs was recorded in suspension culture. Validation of RNA-Seq by qRT-PCR showed upregulated expression of stress-related genes (CrNAC-07 and CrNAC-24), and downregulated expression of growth-related gene (CrNAC-25) in elicited nodal culture of C. roseus. Additionally, the expression of genes involved in the biosynthesis of MIAs was significantly upregulated upon elicitation. The current study provides the first report on the role of CrNACs in regulating the biosynthesis of MIAs.

Comprehensive Analysis of NAC Transcription Factors Reveals Their Evolution in Malvales and Functional Characterization of AsNAC019 and AsNAC098 in Aquilaria sinensis

NAC is a class of plant-specific transcription factors that are widely involved in the growth, development and (a)biotic stress response of plants. However, their molecular evolution has not been extensively studied in Malvales, especially in Aquilaria sinensis, a commercial and horticultural crop that produces an aromatic resin named agarwood. In this study, 1502 members of the NAC gene family were identified from the genomes of nine species from Malvales and three model plants. The macroevolutionary analysis revealed that whole genome duplication (WGD) and dispersed duplication (DSD) have shaped the current architectural structure of NAC gene families in Malvales plants. Then, 111 NAC genes were systemically characterized in A. sinensis. The phylogenetic analysis suggests that NAC genes in A. sinensis can be classified into 16 known clusters and four new subfamilies, with each subfamily presenting similar gene structures and conserved motifs. RNA-seq analysis showed that AsNACs presents a broad transcriptional response to the agarwood inducer. The expression patterns of 15 AsNACs in A. sinensis after injury treatment indicated that AsNAC019 and AsNAC098 were positively correlated with the expression patterns of four polyketide synthase (PKS) genes. Additionally, AsNAC019 and AsNAC098 were also found to bind with the AsPKS07 promoter and activate its transcription. This comprehensive analysis provides valuable insights into the molecular evolution of the NAC gene family in Malvales plants and highlights the potential mechanisms of AsNACs for regulating secondary metabolite biosynthesis in A. sinensis, especially for the biosynthesis of 2-(2-phenyl) chromones in agarwood.

Gene Cloning and Characterization of Transcription Factor FtNAC10 in Tartary Buckwheat (Fagopyrum tataricum (L.) Gaertn.)

NAC transcription factors play a significant role in plant stress responses. In this study, an NAC transcription factor, with a CDS of 792 bp encoding 263 amino acids, was cloned from Fagopyrum tataricum (L.) Gaertn. (F. tataricum), a minor cereal crop, which is rich in flavonoids and highly stress resistant. The transcription factor was named FtNAC10 (NCBI accession number: MK614506.1) and characterized as a member of the NAP subgroup of NAC transcriptions factors. The gene exhibited a highly conserved N-terminal, encoding about 150 amino acids, and a highly specific C-terminal. The resulting protein was revealed to be hydrophilic, with strong transcriptional activation activity. FtNAC10 expression occurred in various F. tataricum tissues, most noticeably in the root, and was regulated differently under various stress treatments. The over-expression of FtNAC10 in transgenic Arabidopsis thaliana (A. thaliana) seeds inhibited germination, and the presence of FtNAC10 enhanced root elongation under saline and drought stress. According to phylogenetic analysis and previous reports, our experiments indicate that FtNAC10 may regulate the stress response or development of F. tataricum through ABA-signaling pathway, although the mechanism is not yet known. This study provides a reference for further analysis of the regulatory function of FtNAC10 and the mechanism that underlies stress responses in Tartary buckwheat.

Cloning and functional analysis of the BrCUC2 gene in Brassica rapa L

The CUP-SHAPED COTYLEDON2 (CUC2) gene plays an important role in the formation of apical meristem and organ edges in plants. The apical meristematic tissue of Brassica rapa (B. rapa) is associated with cold resistance, however, the role of the CUC2 gene in cold resistance of B.rapa is unclear. In this study, we used bioinformatics software to analyze the structure of BrCUC2 gene, real-time fluorescence quantitative PCR to detect the expression level of BrCUC2, constructed transgenic Arabidopsis thaliana by the flower dipping method and subcellular localization for functional validation. The results showed that, we isolated a 1104 bp open reading frame of BrCUC2 from the winter B. rapa cultivar 'Longyou 7'. The BrCUC2 contains a highly conserved domain belonging to the NAM superfamily. Its homologus CUC genes contain similar conserved motifs and are closely related to Brassica oleracea (B.oleracea), and the N-terminal of amino acid sequence contains NAC domain. BrCUC2 protein was localized in the nucleus and self-activation tests showed that pGBKT7-BrCUC2 had self-activation. Tissue-specific expression analysis and promoter β-Glucuronidase (GUS) activity showed that BrCUC2 had high expression levels in B. rapa growth points and A. thaliana leaf edges, stems and growth points. After low-temperature stress, BrCUC2 showed greater expression in 'Longyou 7,' which presents strong cold resistance and concave growth points, than in 'Longyou 99,' which presents weak cold resistance and protruding growth points. BrCUC2 promoter contains multiple elements related to stress responses. BrCUC2 overexpression revealed that the phenotype did not differ from that of the wild type during the seedling stage but showed weak growth and a dwarf phenotype during the flowering and mature stages. After low-temperature treatment, the physiological indexes and survival rate of BrCUC2-overexpression lines of Arabidopsis thaliana (A. thaliana) were better than those of the wild type within 12 h, although differences were not observed after 24 h. These results showed that BrCUC2 improved the low-temperature tolerance of transgenic A. thaliana within a short time. It can provide a foundation for the study of cold resistance in winter B. rapa.

Superoxide dismutase gene family in cassava revealed their involvement in environmental stress via genome-wide analysis

Superoxide dismutase (SOD) is a crucial metal-containing enzyme that plays a vital role in catalyzing the dismutation of superoxide anions, converting them into molecular oxygen and hydrogen peroxide, essential for enhancing plant stress tolerance. We identified 8 SOD genes (4 CSODs, 2 FSODs, and 2 MSODs) in cassava. Bioinformatics analyses provided insights into chromosomal location, phylogenetic relationships, gene structure, conserved motifs, and gene ontology annotations. MeSOD genes were classified into two groups through phylogenetic analysis, revealing evolutionary connections. Promoters of these genes harbored stress-related cis-elements. Duplication analysis indicated the functional significance of MeCSOD2/MeCSOD4 and MeMSOD1/MeMSOD2. Through qRT-PCR, MeCSOD2 responded to salt stress, MeMSOD2 to drought, and cassava bacterial blight. Silencing MeMSOD2 increased XpmCHN11 virulence, indicating MeMSOD2 is essential for cassava's defense against XpmCHN11 infection. These findings enhance our understanding of the SOD gene family's role in cassava and contribute to strategies for stress tolerance improvement.

Genome-Wide Analysis of WRKY and NAC Transcription Factors in Carica papaya L. and Their Possible Role in the Loss of Drought Tolerance by Recent Cultivars through the Domestication of Their Wild Ancestors

A genome-wide analysis for two families of key transcription factors (TF; WRKY and NAC) involved in drought response revealed 46 WRKY and 66 NAC members of the Carica papaya genome. A phylogenetic analysis grouped the CpWRKY proteins into three groups (I, II a, b, c, d, e and III), while the CpNAC proteins were clustered into 15 groups. The conserved domains, chromosomal localization and promoter cis-acting elements were also analyzed. In addition, from a previous transcriptome study of two contrasting genotypes in response to 14 days of water deficit stress (WDS), we found that 29 of the 46 CpWRKYs genes and 25 of the 66 CpNACs genes were differentially expressed in response to the WDS. In the present paper, the native wild genotype (WG) (collected in its center of origin) consistently showed a higher expression (transcripts per million; TPM and fold change; FC) than the commercial genotype (CG) in almost all the members of the CpWRKY and CpNAC gene families. To corroborate this, we selected CpWRKY50 and CpNAC83.1 for further evaluation by RT-qPCR. Consistently, the WG showed higher relative expression levels (REL) after 14 days of WDS than the CG, in both the leaves and roots. The results suggest that the CpWRKY and CpNAC TF families are important for drought tolerance in this species. The results may also suggest that, during the domestication process, the ability of the native (wild) C. papaya genotypes to respond to drought (including the overexpression of the CpWRKY and CpNAC genes) was somehow reduced in the current commercial genotypes.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Gene Family in Litchi chinensis

NAC proteins play an essential role in the growth and development of litchi, especially during reproductive development. However, a comprehensive analysis of the litchi NAC gene family is currently absent. Based on information from the litchi genome, we found that the 112 NAC genes of litchi show an uneven distribution on the chromosomes. Phylogenetic and conserved structural domain analyses indicated that different types of variability were exhibited in the family of litchi NACs (LcNACs). Gene covariance analysis showed that the LcNACs showed better similarity in the same genus than with Arabidopsis. We further investigated the differential expression patterns of LcNACs in buds and rudimentary leaves of litchi. qRT-PCR results implied that they were involved in the process. Profiling of LcNAC promoter elements in litchi showed that they were extensively involved in light response, phytohormone regulation, abiotic stress response, and plant growth and development processes. This study provides new insights into the identification, structural characterization, tissue-specific expression analysis, and promoter response elements of LcNACs. It reveals the characteristics of the LcNACs and lays the foundation for the subsequent understanding of its biological functions and molecular regulatory mechanisms.

Genome-Wide Identification of NAC Gene Family and Expression Analysis under Abiotic Stresses in Avena sativa

In this study, a total of 177 NAC members were identified in Avena sativa, located on 21 chromosomes. Phylogenetic analysis showed that AsNAC proteins could be divided into seven subfamilies (I-VII), and that proteins in the same subfamily have similar protein motifs. Gene structure analysis found that NAC introns ranged from 1 to 17. Cis-element analysis of the promoter indicated that the gene family may have stress-related elements and growth regulation elements. Through qRT-PCR experiments, we speculated that AsNACs genes can respond to abiotic stresses such as cold, freezing, salt, and saline alkali. This study provides a theoretical basis for further exploring the function of the NAC gene family in A. sativa.

Genome-wide identification and expression analysis of NAC family genes in Ginkgo biloba L

NAC (NAM, ATAF, CUC2) transcription factors constitute one of the largest families of plant-specific transcription factors with important roles in plant growth and development and in biotic and abiotic stresses. The physicochemical properties, gene structure, cis-acting elements and expression patterns of NAC transcription factors in Ginkgo biloba were analysed using bioinformatics, and expression of this gene family was analysed via quantitative reverse transcription PCR. The family of G. biloba NAC transcription factors had 50 members, distributed on 12 chromosomes and divided into 11 groups. Members in the same group share a similar gene structure and motif distribution. Transcriptome data analysis of G. biloba showed that 35 genes were expressed in eight tissues. Correlation analysis suggested that GbNAC007 and GNAC008 might be involved in flavonoid biosynthesis. Expression levels of 12 GbNACs under cold, het, and salt stresses were analysed. Results indicate that NAC transcription factors play an important role in response to abiotic stresses. This study provides a reference for the functional analysis of the G. biloba family of NAC transcription factors, as well as a resource for studies on the involvement of this family in responses to abiotic stresses and flavonoid biosynthesis.

Analysis of proteomic changes in cassava cv. Kasetsart 50 caused by Sri Lankan cassava mosaic virus infection

BACKGROUND

Sri Lankan cassava mosaic virus (SLCMV) is a plant virus causing significant economic losses throughout Southeast Asia. While proteomics has the potential to identify molecular markers that could assist the breeding of virus resistant cultivars, the effects of SLCMV infection in cassava have not been previously explored in detail.

RESULTS

Liquid Chromatography-Tandem Mass Spectrometry (LC/MS-MS) was used to identify differentially expressed proteins in SLCMV infected leaves, and qPCR was used to confirm changes at mRNA levels. LC/MS-MS identified 1,813 proteins, including 479 and 408 proteins that were upregulated in SLCMV-infected and healthy cassava plants respectively, while 109 proteins were detected in both samples. Most of the identified proteins were involved in biosynthetic processes (29.8%), cellular processes (20.9%), and metabolism (18.4%). Transport proteins, stress response molecules, and proteins involved in signal transduction, plant defense responses, photosynthesis, and cellular respiration, although present, only represented a relatively small subset of the detected differences. RT-qPCR confirmed the upregulation of WRKY 77 (A0A140H8T1), WRKY 83 (A0A140H8T7), NAC 6 (A0A0M4G3M4), NAC 35 (A0A0M5JAB4), NAC 22 (A0A0M5J8Q6), NAC 54 (A0A0M4FSG8), NAC 70 (A0A0M4FEU9), MYB (A0A2C9VER9 and A0A2C9VME6), bHLH (A0A2C9UNL9 and A0A2C9WBZ1) transcription factors. Additional upregulated transcripts included receptors, such as receptor-like serine/threonine-protein kinase (RSTK) (A0A2C9UPE4), Toll/interleukin-1 receptor (TIR) (A0A2C9V5Q3), leucine rich repeat N-terminal domain (LRRNT_2) (A0A2C9VHG8), and cupin (A0A199UBY6). These molecules participate in innate immunity, plant defense mechanisms, and responses to biotic stress and to phytohormones.

CONCLUSIONS

We detected 1,813 differentially expressed proteins infected cassava plants, of which 479 were selectively upregulated. These could be classified into three main biological functional groups, with roles in gene regulation, plant defense mechanisms, and stress responses. These results will help identify key proteins affected by SLCMV infection in cassava plants.

Overexpression of MePMEI1 in Arabidopsis enhances Pb tolerance

Pb is one of the most ubiquitously distributed heavy metal pollutants in soils and has serious negative effects on plant growth, food safety, and public health. Pectin methylesterase inhibitors (PMEIs) play a pivotal role in regulating the integrity of plant cell walls; however, the molecular basis by which PMEIs promote plant resistance to abiotic stress remains poorly understood. In this study, we identified a novel PMEI gene, MePMEI1, from Manihot esculenta, and determined its role in plant resistance to Pb stress. The expression of MePMEI1 was remarkably upregulated in the roots, stems, and leaves of cassava plants following exposure to Pb stress. An analysis of subcellular localization revealed that the MePMEI1 protein was localized in the cell wall. MePMEI1 inhibited commercial orange peel pectin methyltransferase (PME), and the expression of MePMEI1 in Arabidopsis decreased the PME activity, indicating that MePMEI1 can inhibit PME activity in the cell wall. Additionally, the overexpression of MePMEI1 in Arabidopsis reduced oxidative damage and induced the thickening of cell walls, thus contributing to Pb tolerance. Altogether, the study reports a novel mechanism by which the MePMEI1 gene, which encodes the PMEI protein in cassava, plays an essential role in promoting tolerance to Pb toxicity by regulating the thickness of cell walls. These results provide a theoretical basis for the MePMEI1-mediated plant breeding for increasing heavy metal tolerance and provide insights into controlling Pb pollution in soils through phytoremediation in future studies.

Overexpression of cotton GhNAC072 gene enhances drought and salt stress tolerance in transgenic Arabidopsis

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

In summary, the above analysis showed that GhNAC072 might be the true candidate involved in boosting tolerance mechanisms under drought and salinity stress.

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.

Current status and impending progress for cassava structural genomics

We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Genome-Wide Identification of the NAC Gene Family in Zanthoxylum bungeanum and Their Transcriptional Responses to Drought Stress

NAC (NAM, ATAF1/2, and CUC2) transcription factors (TFs) are one of the largest plant-specific TF families and play a pivotal role in adaptation to abiotic stresses. The genome-wide analysis of NAC TFs is still absent in Zanthoxylum bungeanum. Here, 109 ZbNAC proteins were identified from the Z. bungeanum genome and were classified into four groups with Arabidopsis NAC proteins. The 109 ZbNAC genes were unevenly distributed on 46 chromosomes and included 4 tandem duplication events and 17 segmental duplication events. Synteny analysis of six species pairs revealed the closely phylogenetic relationship between Z. bungeanum and C. sinensis. Twenty-four types of cis-elements were identified in the ZbNAC promoters and were classified into three types: abiotic stress, plant growth and development, and response to phytohormones. Co-expression network analysis of the ZbNACs revealed 10 hub genes, and their expression levels were validated by real-time quantitative polymerase chain reaction (qRT-PCR). Finally, ZbNAC007, ZbNAC018, ZbNAC047, ZbNAC072, and ZbNAC079 were considered the pivotal NAC genes for drought tolerance in Z. bungeanum. This study represented the first genome-wide analysis of the NAC family in Z. bungeanum, improving our understanding of NAC proteins and providing useful information for molecular breeding of Z. bungeanum.

The future of cassava in the era of biotechnology in Southern Africa

Cassava (Manihot esculenta) is a major staple food and the world's fourth source of calories. Biotechnological contributions to enhancing this crop, its advances, and present issues must be assessed regularly. Functional genomics, genomic-assisted breeding, molecular tools, and genome editing technologies, among other biotechnological approaches, have helped improve the potential of economically important crops like cassava by addressing some of its significant constraints, such as nutrient deficiency, toxicity, poor starch quality, disease susceptibility, low yield capacity, and postharvest deterioration. However, the development, improvement, and subsequent acceptance of the improved cultivars have been challenging and have required holistic approaches to solving them. This article provides an update of trends and gaps in cassava biotechnology, reviewing the relevant strategies used to improve cassava crops and highlighting the potential risk and acceptability of improved cultivars in Southern Africa.

A population based expression atlas provides insights into disease resistance and other physiological traits in cassava (Manihot esculenta Crantz)

Cassava, a food security crop in Africa, is grown throughout the tropics and subtropics. Although cassava can provide high productivity in suboptimal conditions, the yield in Africa is substantially lower than in other geographies. The yield gap is attributable to many challenges faced by cassava in Africa, including susceptibility to diseases and poor soil conditions. In this study, we carried out 3’RNA sequencing on 150 accessions from the National Crops Resources Research Institute, Uganda for 5 tissue types, providing population-based transcriptomics resources to the research community in a web-based queryable cassava expression atlas. Differential expression and weighted gene co-expression network analysis were performed to detect 8820 significantly differentially expressed genes (DEGs), revealing similarity in expression patterns between tissue types and the clustering of detected DEGs into 18 gene modules. As a confirmation of data quality, differential expression and pathway analysis targeting cassava mosaic disease (CMD) identified 27 genes observed in the plant–pathogen interaction pathway, several previously identified CMD resistance genes, and two peroxidase family proteins different from the CMD2 gene. Present research work represents a novel resource towards understanding complex traits at expression and molecular levels for the development of resistant and high-yielding cassava varieties, as exemplified with CMD.

AvNAC030, a NAC Domain Transcription Factor, Enhances Salt Stress Tolerance in Kiwifruit

Kiwifruit (Actinidia chinensis Planch) is suitable for neutral acid soil. However, soil salinization is increasing in kiwifruit production areas, which has adverse effects on the growth and development of plants, leading to declining yields and quality. Therefore, analyzing the salt tolerance regulation mechanism can provide a theoretical basis for the industrial application and germplasm improvement of kiwifruit. We identified 120 NAC members and divided them into 13 subfamilies according to phylogenetic analysis. Subsequently, we conducted a comprehensive and systematic analysis based on the conserved motifs, key amino acid residues in the NAC domain, expression patterns, and protein interaction network predictions and screened the candidate gene AvNAC030. In order to study its function, we adopted the method of heterologous expression in Arabidopsis. Compared with the control, the overexpression plants had higher osmotic adjustment ability and improved antioxidant defense mechanism. These results suggest that AvNAC030 plays a positive role in the salt tolerance regulation mechanism in kiwifruit.

Identification and analysis of Chrysanthemum nankingense NAC transcription factors and an expression analysis of OsNAC7 subfamily members

NAC (NAM, ATAF1-2, and CUC2) transcription factors (TFs) play a vital role in plant growth and development, as well as in plant response to biotic and abiotic stressors (Duan et al., 2019; Guerin et al., 2019). Chrysanthemum is a plant with strong stress resistance and adaptability; therefore, a systematic study of NAC TFs in chrysanthemum is of great significance for plant breeding. In this study, 153 putative NAC TFs were identified based on the Chrysanthemum nankingense genome. According to the NAC family in Arabidopsis and rice, a rootless phylogenetic tree was constructed, in which the 153 CnNAC TFs were divided into two groups and 19 subfamilies. Moreover, the expression levels of 12 CnNAC TFs belonging to the OsNAC7 subfamily were analyzed in C. nankingense under osmotic and salt stresses, and different tissues were tested during different growth periods. The results showed that these 12 OsNAC7 subfamily members were involved in the regulation of root and stem growth, as well as in the regulation of drought and salt stresses. Finally, we investigated the function of the CHR00069684 gene, and the results showed that CHR00069684 could confer improved salt and low temperature resistance, enhance ABA sensitivity, and lead to early flowering in tobacco. It was proved that members of the OsNAC7 subfamily have dual functions including the regulation of resistance and the mediation of plant growth and development. This study provides comprehensive information on analyzing the function of CnNAC TFs, and also reveals the important role of OsNAC7 subfamily genes in response to abiotic stress and the regulation of plant growth. These results provide new ideas for plant breeding to control stress resistance and growth simultaneously.

Genome-wide identification, characterization, and expression analysis of the NAC transcription factor family in orchardgrass (Dactylis glomerata L.)

Background

Orchardgrass (Dactylis glomerata L.) is one of the most important cool-season perennial forage grasses that is widely cultivated in the world and is highly tolerant to stressful conditions. However, little is known about the mechanisms underlying this tolerance. The NAC (NAM, ATAF1/2, and CUC2) transcription factor family is a large plant-specific gene family that actively participates in plant growth, development, and response to abiotic stress. At present, owing to the absence of genomic information, NAC genes have not been systematically studied in orchardgrass. The recent release of the complete genome sequence of orchardgrass provided a basic platform for the investigation of DgNAC proteins.

Results

Using the recently released orchardgrass genome database, a total of 108 NAC (DgNAC) genes were identified in the orchardgrass genome database and named based on their chromosomal location. Phylogenetic analysis showed that the DgNAC proteins were distributed in 14 subgroups based on homology with NAC proteins in Arabidopsis, including the orchardgrass-specific subgroup Dg_NAC. Gene structure analysis suggested that the number of exons varied from 1 to 15, and multitudinous DgNAC genes contained three exons. Chromosomal mapping analysis found that the DgNAC genes were unevenly distributed on seven orchardgrass chromosomes. For the gene expression analysis, the expression levels of DgNAC genes in different tissues and floral bud developmental stages were quite different. Quantitative real-time PCR analysis showed distinct expression patterns of 12 DgNAC genes in response to different abiotic stresses. The results from the RNA-seq data revealed that orchardgrass-specific NAC exhibited expression preference or specificity in diverse abiotic stress responses, and the results indicated that these genes may play an important role in the adaptation of orchardgrass under different environments.

Conclusions

In the current study, a comprehensive and systematic genome-wide analysis of the NAC gene family in orchardgrass was first performed. A total of 108 NAC genes were identified in orchardgrass, and the expression of NAC genes during plant growth and floral bud development and response to various abiotic stresses were investigated. These results will be helpful for further functional characteristic descriptions of DgNAC genes and the improvement of orchardgrass in breeding programs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07485-6.

Comprehensive analysis of NAC transcription factor family uncovers drought and salinity stress response in pearl millet (Pennisetum glaucum)

BACKGROUND

Pearl millet (Pennisetum glaucum) is a cereal crop that possesses the ability to withstand drought, salinity and high temperature stresses. The NAC [NAM (No Apical Meristem), ATAF1 (Arabidopsis thaliana Activation Factor 1), and CUC2 (Cup-shaped Cotyledon)] transcription factor family is one of the largest transcription factor families in plants. NAC family members are known to regulate plant growth and abiotic stress response. Currently, no reports are available on the functions of the NAC family in pearl millet.

RESULTS

Our genome-wide analysis found 151 NAC transcription factor genes (PgNACs) in the pearl millet genome. Thirty-eight and 76 PgNACs were found to be segmental and dispersed duplicated respectively. Phylogenetic analysis divided these NAC transcription factors into 11 groups (A-K). Three PgNACs (- 073, - 29, and - 151) were found to be membrane-associated transcription factors. Seventeen other conserved motifs were found in PgNACs. Based on the similarity of PgNACs to NAC proteins in other species, the functions of PgNACs were predicted. In total, 88 microRNA target sites were predicted in 59 PgNACs. A previously performed transcriptome analysis suggests that the expression of 30 and 42 PgNACs are affected by salinity stress and drought stress, respectively. The expression of 36 randomly selected PgNACs were examined by quantitative reverse transcription-PCR. Many of these genes showed diverse salt- and drought-responsive expression patterns in roots and leaves. These results confirm that PgNACs are potentially involved in regulating abiotic stress tolerance in pearl millet.

CONCLUSION

The pearl millet genome contains 151 NAC transcription factor genes that can be classified into 11 groups. Many of these genes are either upregulated or downregulated by either salinity or drought stress and may therefore contribute to establishing stress tolerance in pearl millet.

An OsNAM gene plays important role in root rhizobacteria interaction in transgenic Arabidopsis through abiotic stress and phytohormone crosstalk

Overexpression of Bacillus amyloliquefaciens SN13-responsive OsNAM gene in Arabidopsis reveals its important role in beneficial plant and plant growth promoting rhizobacteria interaction by conferring stress tolerance and phytohormone modulation. Salinity is one of the major constraints that affect crop development and yield. Plants respond and adapt to salt stress via complex mechanisms that involve morpho-physiological, biochemical, and molecular changes. The expression of numerous genes is known to alter during various abiotic stresses and impart stress tolerance. Recently, some known rhizospheric microbes have also been used to mitigate the effects of abiotic stresses; however, the molecular basis of such interactions remains elusive. Therefore, the present investigation was aimed to elucidate the plant growth-promoting rhizobacteria (PGPR; Bacillus amyloliquefaciens-SN13) -induced crosstalk among salinity and phytohormones in OsNAM-overexpressed Arabidopsis plants. Transgenic plants showed increased germination percentage compared to wild-type (WT) seeds under 100 mM of NaCl. Phenotypic data showed increased root length, rosette diameter, leaf size, and biomass in transgenics than WT plants. Transgenic plants can also better maintain membrane integrity and osmolyte concentration under salinity as compared to WT. Further, gene expression analysis of AP2/ERF, GST, ERD4, and ARF2 genes showed differential expression and their positive modulation in transgenic Arabidopsis exposed to salt stress in the presence of SN13 as compared to uninoculated WT. Modulation in IAA, ABA, and GA content in inoculated plants showed the more pronounced positive effects of SN13 on transgenic plants that supported our findings on Arabidopsis-SN13 interaction. Overall, the study concludes that SN13 positively modulated expression of stress-responsive genes under salinity and alter phytohormones levels in OsNAM-overexpressed plants suggesting its extensive role in cross-talk among salinity and phytohormones in response to PGPR.

Genome-Wide Identification and Gene Expression Analysis of Acyl-Activating Enzymes Superfamily in Tomato (Solanum lycopersicum) Under Aluminum Stress

In response to changing environments, plants regulate gene expression and subsequent metabolism to acclimate and survive. A superfamily of acyl-activating enzymes (AAEs) has been observed in every class of creatures on planet. Some of plant AAE genes have been identified and functionally characterized to be involved in growth, development, biotic, and abiotic stresses via mediating diverse metabolic pathways. However, less information is available about AAEs superfamily in tomato (Solanum lycopersicum), the highest value fruit and vegetable crop globally. In this study, we aimed to identify tomato AAEs superfamily and investigate potential functions with respect to aluminum (Al) stress that represents one of the major factors limiting crop productivity on acid soils worldwide. Fifty-three AAE genes of tomato were identified and named on the basis of phylogenetic relationships between Arabidopsis and tomato. The phylogenetic analysis showed that AAEs could be classified into six clades; however, clade III contains no AAE genes of tomato. Synteny analyses revealed tomato vegetable paralogs and Arabidopsis orthologs. The RNA-seq and quantitative reverse-transcriptase PCR (qRT-PCR) analysis indicated that 9 out of 53 AAEs genes were significantly up- or downregulated by Al stress. Numerous cis-acting elements implicated in biotic and abiotic stresses were detected in the promoter regions of SlAAEs. As the most abundantly expressed gene in root apex and highly induced by Al, there are many potential STOP1 cis-acting elements present in the promoter of SlAAE3-1, and its expression in root apex was specific to Al. Finally, transgenic tobacco lines overexpressing SlAAE3-1 displayed increased tolerance to Al. Altogether, our results pave the way for further studies on the functional characterization of SlAAE genes in tomato with a wish of improvement in tomato crop in the future.

Genomic analyses of heat stress transcription factors (HSFs) in simulated drought stress response and storage root deterioration after harvest in cassava

Heat shock factors (HSFs) play crucial roles in various plant stress responses. However, the current knowledge about HSFs in cassava, an important crop, is still insufficient. In this research, we identified 32 cassava HSF genes (MeHSFs) and clustered them into three groups (A, B, C) based on phylogenetic analysis and structural characteristics. Conserved motif analyses showed that MeHSFs display domains characteristic to HSF transcription factors. Gene structure analyses suggested that 29 MeHSFs contained only two exons. All identified 32 cassava MeHSFs were distributed on 13 chromosomes. Their expression profiles revealed that the different MeHSFs were expressed differentially in different tissues, most high expression genes belonged to group A. The similar MeHSFs were up-regulated after treatment with both PEG and abscisic acid (ABA), which implied that these MeHSFs may participate in resistance to simulated drought stress associated with the ABA signaling pathway. In addition, several MeHSFs were induced during postharvest physiological deterioration (PPD) in cassava. Our results provided basic but important knowledge for future gene function analysis of MeHSFs toward efforts in improving tolerance to abiotic stress and PPD in cassava.

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.

Genome-wide analysis of the NAC transcription factor family in broomcorn millet (Panicum miliaceum L.) and expression analysis under drought stress

BACKGROUND

Broomcorn millet is a drought-tolerant cereal that is widely cultivated in the semiarid regions of Asia, Europe, and other continents; however, the mechanisms underlying its drought-tolerance are poorly understood. The NAM, ATAF1/2, and CUC2 (NAC) transcription factors form a large plant-specific gene family that is involved in the regulation of tissue development and abiotic stress. To date, NAC transcription factors have not been systematically researched in broomcorn millet.

RESULTS

In the present study, a total of 180 NAC (PmNAC) genes were identified from the broomcorn millet genome and named uniformly according to their chromosomal distribution. Phylogenetic analysis demonstrated that the PmNACs clustered into 12 subgroups, including the broomcorn millet-specific subgroup Pm_NAC. Gene structure and protein motif analyses indicated that closely clustered PmNAC genes were relatively conserved within each subgroup, while genome mapping analysis revealed that the PmNAC genes were unevenly distributed on broomcorn millet chromosomes. Transcriptome analysis revealed that the PmNAC genes differed greatly in expression in various tissues and under different drought stress durations. The expression of 10 selected genes under drought stress was analyzed using quantitative real-time PCR.

CONCLUSION

In this study, 180 NAC genes were identified in broomcorn millet, and their phylogenetic relationships, gene structures, protein motifs, chromosomal distribution, duplication, expression patterns in different tissues, and responses to drought stress were studied. These results will be useful for the further study of the functional characteristics of PmNAC genes, particularly with regards to drought resistance.

Genome-Wide Identification and Characterization of the NAC Transcription Factor Family in Musa Acuminata and Expression Analysis during Fruit Ripening

Banana (Musa acuminata, AAA group) is a representative climacteric fruit with essential nutrients and pleasant flavors. Control of its ripening determines both the fruit quality and the shelf life. NAC (NAM, ATAF, CUC2) proteins, as one of the largest superfamilies of transcription factors, play crucial roles in various functions, especially developmental processes. Thus, it is important to conduct a comprehensive identification and characterization of the NAC transcription factor family at the genomic level in M. acuminata. In this article, a total of 181 banana NAC genes were identified. Phylogenetic analysis indicated that NAC genes in M. acuminata, Arabidopsis, and rice were clustered into 18 groups (S1–S18), and MCScanX analysis disclosed that the evolution of MaNAC genes was promoted by segmental duplication events. Expression patterns of NAC genes during banana fruit ripening induced by ethylene were investigated using RNA-Seq data, and 10 MaNAC genes were identified as related to fruit ripening. A subcellular localization assay of selected MaNACs revealed that they were all localized to the nucleus. These results lay a good foundation for the investigation of NAC genes in banana toward the biological functions and evolution.

Identification and Expression Analysis of the NAC Gene Family in Coffea canephora

The NAC gene family is one of the largest families of transcriptional regulators in plants, and it plays important roles in the regulation of growth and development as well as in stress responses. Genome-wide analyses have been performed in diverse plant species, but there is still no systematic analysis of the NAC genes of Coffea canephora Pierre ex A. Froehner. In this study, we identified 63 NAC genes from the genome of C. canephora. The basic features and comparison analysis indicated that the NAC gene members increased via duplication events during the evolution of the plant. Phylogenetic analysis divided the NAC proteins from C. canephora, Arabidopsis and rice into 16 subgroups. Analysis of the expression patterns of CocNACs under cold stress and coffee bean development indicated that 38 CocNACs were differentially expressed under cold stress; six genes may play important roles in the process of cold acclimation, and four genes among 54 CocNACs showing a variety of expression patterns during different developmental stages of coffee beans may be positively related to the bean development. This study can expand our understanding of the functions of the CocNAC gene family in cold responses and bean development, thereby potentially intensifying the molecular breeding programs of Coffea spp. plants.

Genome-wide analysis of the NF-Y gene family in peach (Prunus persica L.)

Background

Nuclear Factor Y (NF-Y) is a heterotrimeric complex composed of three unique subunits: NF-YA, NF-YB, and NF-YC. The NF-Y transcription factor complex binds to the CCAAT box of eukaryotic promoters, playing a vital role in various biological processes in plants. However, the NF-Y gene family has not yet been reported from the peach genome. The current study identified and classified candidate peach NF-Y genes for further functional analysis of this family.

Results

The current study identified 24 Nuclear Factor Y (NF-Y) transcription factor subunits (6 NF-YA, 12 NF-YB, and 6 NF-YC subunits) in peach. These NF-Y subunits were described with respect to basic physicochemical characteristics, chromosome locations, gene structures, and conserved domains. Based on an analysis of the phylogenetic relationships among peach NF-Ys, six pairs of paralogous NF-Ys were detected. The expansion of the peach NF-Y family occurred by segmental and tandem duplication. Phylogenetic gene synteny of NF-Y proteins was observed between peach and Arabidopsis, and five pairs of paralogous NF-Y proteins from peach and Arabidopsis were identified. Twenty-four peach NF-Ys displayed a diversity of tissue expression patterns. In addition, drought-responsive cis-elements were observed in peach NF-Y promoters, and 9 peach NF-Y genes were shown to distinctly increase their transcript abundances under drought stress.

Conclusions

This study identified 24 NF-Y genes in the peach genome and analysed their properties at different levels, providing a foundation for researchers to understand this gene family in peach. The up-regulation of 9 NF-Y genes under drought stress indicates that they can serve as candidate functional genes to further study drought resistance in peach.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5968-7) contains supplementary material, which is available to authorized users.

Global Dynamics in Protein Disorder during Maize Seed Development

Intrinsic protein disorder is a physicochemical attribute of some proteins lacking tridimensional structure and is collectively known as intrinsically disordered proteins (IDPs). Interestingly, several IDPs have been associated with protective functions in plants and with their response to external stimuli. To correlate the modulation of the IDPs content with the developmental progression in seed, we describe the expression of transcripts according to the disorder content of the proteins that they codify during seed development, from the early embryogenesis to the beginning of the desiccation tolerance acquisition stage. We found that the total expression profile of transcripts encoding for structured proteins is highly increased during middle phase. However, the relative content of protein disorder is increased as seed development progresses. We identified several intrinsically disordered transcription factors that seem to play important roles throughout seed development. On the other hand, we detected a gene cluster encoding for IDPs at the end of the late phase, which coincides with the beginning of the acquisition of desiccation tolerance. In conclusion, the expression pattern of IDPs is highly dependent on the developmental stage, and there is a general reduction in the expression of transcripts encoding for structured proteins as seed development progresses. We proposed maize seeds as a model to study the regulation of protein disorder in plant development and its involvement in the acquisition of desiccation tolerance in plants.

Genome-Wide Identification, Characterization, and Expression Analysis of the NAC Transcription Factor in Chenopodium quinoa

The NAC (NAM, ATAF, and CUC) family is one of the largest families of plant-specific transcription factors. It is involved in many plant growth and development processes, as well as abiotic/biotic stress responses. So far, little is known about the NAC family in Chenopodium quinoa. In the present study, a total of 90 NACs were identified in quinoa (named as CqNAC1-CqNAC90) and phylogenetically divided into 14 distinct subfamilies. Different subfamilies showed diversities in gene proportions, exon-intron structures, and motif compositions. In addition, 28 CqNAC duplication events were investigated, and a strong subfamily preference was found during the NAC expansion in quinoa, indicating that the duplication event was not random across NAC subfamilies during quinoa evolution. Moreover, the analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios suggested that the duplicated CqNACs might have mainly experienced purifying selection pressure with limited functional divergence. Additionally, 11 selected CqNACs showed significant tissue-specific expression patterns, and all the CqNACs were positively regulated in response to salt stress. The result provided evidence for selecting candidate genes for further characterization in tissue/organ specificity and their functional involvement in quinoa's strong salinity tolerance.

The Class III Peroxidase (POD) Gene Family in Cassava: Identification, Phylogeny, Duplication, and Expression

The class III peroxidase (POD) enzymes participate in plant development, hormone signaling, and stress responses. However, little is known about the POD family in cassava. Here, we identified 91 cassava POD genes (MePODs) and classified them into six subgroups using phylogenetic analysis. Conserved motif analysis demonstrated that all MePOD proteins have typical peroxidase domains, and gene structure analysis showed that MePOD genes have between one and nine exons. Duplication pattern analysis suggests that tandem duplication has played a role in MePOD gene expansion. Comprehensive transcriptomic analysis revealed that MePOD genes in cassava are involved in the drought response and postharvest physiological deterioration. Several MePODs underwent transcriptional changes after various stresses and related signaling treatments were applied. In sum, we characterized the POD family in cassava and uncovered the transcriptional control of POD genes in response to various stresses and postharvest physiological deterioration conditions. These results can be used to identify potential target genes for improving the stress tolerance of cassava crops.

Differential expression of the Hsf family in cassava under biotic and abiotic stresses

Heat shock transcription factors (Hsfs) are important regulators of biotic and abiotic stress responses in plants. Currently, the Hsf gene family is not well understood in cassava, an important tropical crop. In the present study, 32 MeHsf genes were identified from the cassava genome database, which were divided into three types based on functional domain and motif distribution analyses. Analysis of the differential expression of the genes belonging to the Hsf family in cassava was carried out based on published cassava transcriptome data from tissues/organs (leaf blade, leaf midvein, lateral buds, organized embryogenic structures, friable embryogenic callus, fibrous roots, storage roots, stem, petiole, shoot apical meristem, and root apical meristem) under abiotic stress (cold, drought) or biotic stress (mealybugs. cassava brown streak disease, cassava bacterial blight). The results show the expression diversity of cassava Hsfs genes in various tissues/organs. The transcript levels of MeHsfB3a, MeHsfA6a, MeHsfA2a, and MeHsfA9b were upregulated by abiotic and biotic stresses, such as cold, drought, cassava bacterial blight, cassava brown streak disease, and mealybugs, indicating their potential roles in mediating the response of cassava plants to environment stresses. Further interaction network and co-expression analyses suggests that Hsf genes may interact with Hsp70 family members to resist environmental stresses in cassava. These results provide valuable information for future studies of the functional characterization of the MeHsf gene family.

Genome-wide analysis of the NAC transcription factor family in Tartary buckwheat (Fagopyrum tataricum)

Background

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family represents a group of large plant-specific transcriptional regulators, participating in plant development and response to external stress. However, there is no comprehensive study on the NAC genes of Tartary buckwheat (Fagopyrum tataricum), a large group of extensively cultivated medicinal and edible plants. The recently published Tartary buckwheat genome permits us to explore all the FtNAC genes on a genome-wide basis.

Results

In the present study, 80 NAC (FtNAC) genes of Tartary buckwheat were obtained and named uniformly according to their distribution on chromosomes. Phylogenetic analysis of NAC proteins in both Tartary buckwheat and Arabidopsis showed that the FtNAC proteins are widely distributed in 15 subgroups with one subgroup unclassified. Gene structure analysis found that multitudinous FtNAC genes contained three exons, indicating that the structural diversity in Tartary buckwheat NAC genes is relatively low. Some duplication genes of FtNAC have a conserved structure that was different from others, indicating that these genes may have a variety of functions. By observing gene expression, we found that FtNAC genes showed abundant differences in expression levels in various tissues and at different stages of fruit development.

Conclusions

In this research, 80 NAC genes were identified in Tartary buckwheat, and their phylogenetic relationships, gene structures, duplication, global expression and potential roles in Tartary buckwheat development were studied. Comprehensive analysis will be useful for a follow-up study of functional characteristics of FtNAC genes and for the development of high-quality Tartary buckwheat varieties.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5500-0) contains supplementary material, which is available to authorized users.

Identification and characterization of drought-responsive CC-type glutaredoxins from cassava cultivars reveals their involvement in ABA signalling

BACKGROUND

CC-type glutaredoxins (GRXs) are plant-specific glutaredoxin, play regulatory roles in response of biotic and abiotic stress. However, it is not clear whether the CC-type GRXs are involve in drought response in cassava (Manihot esculenta), an important tropical tuber root crop.

RESULTS

Herein, genome-wide analysis identified 18 CC-type GRXs in the cassava genome, of which six (namely MeGRXC3, C4, C7, C14, C15, and C18) were induced by drought stress in leaves of two cassava cultivars Argentina 7 (Arg7) and South China 124 (SC124). Exogenous abscisic acid (ABA) application induced the expression of all the six CC-type GRXs in leaves of both Arg7 and SC124 plants. Overexpression of MeGRXC15 in Arabidopsis (Col-0) increases tolerance of ABA on the sealed agar plates, but results in drought hypersensitivity in soil-grown plants. The results of microarray assays show that MeGRXC15 overexpression affected the expression of a set of transcription factors which involve in stress response, ABA, and JA/ET signalling pathway. The results of protein interaction analysis show that MeGRXC15 can interact with TGA5 from Arabidopsis and MeTGA074 from cassava.

CONCLUSIONS

CC-type glutaredoxins play regulatory roles in cassava response to drought possibly through ABA signalling pathway.

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.

Functional analysis of MeCIPK23 and MeCBL1/9 in cassava defense response against Xanthomonas axonopodis pv. manihotis

KEY MESSAGE

MeCIPK23 interacts with MeCBL1/9, and they confer improved defense response, providing potential genes for further genetic breeding in cassava. Cassava (Manihot esculenta) is an important food crop in tropical area, but its production is largely affected by cassava bacterial blight. However, the information of defense-related genes in cassava is very limited. Calcium ions play essential roles in plant development and stress signaling pathways. Calcineurin B-like proteins (CBLs) and CBL-interacting protein kinases (CIPKs) are crucial components of calcium signals. In this study, systematic expression profile of 25MeCIPKs in response to Xanthomonas axonopodis pv. manihotis (Xam) infection was examined, by which seven candidate MeCIPKs were chosen for functional investigation. Through transient expression in Nicotiana benthamiana leaves, we found that six MeCIPKs (MeCIPK5, MeCIPK8, MeCIPK12, MeCIPK22, MeCIPK23 and MeCIPK24) conferred improved defense response, via regulating the transcripts of several defense-related genes. Notably, we found that MeCIPK23 interacted with MeCBL1 and MeCBL9, and overexpression of these genes conferred improved defense response. On the contrary, virus-induced gene silencing of either MeCIPK23 or MeCBL1/9 or both genes resulted in disease sensitive in cassava. To our knowledge, this is the first study identifying MeCIPK23 as well as MeCBL1 and MeCBL9 that confer enhanced defense response against Xam.

The Late Embryogenesis Abundant Protein Family in Cassava (Manihot esculenta Crantz): Genome-Wide Characterization and Expression during Abiotic Stress

Late embryogenesis abundant (LEA) proteins, as a highly diverse group of polypeptides, play an important role in plant adaptation to abiotic stress; however, LEAs from cassava have not been studied in cassava. In this study, 26 LEA members were genome-wide identified from cassava, which were clustered into seven subfamily according to evolutionary relationship, protein motif, and gene structure analyses. Chromosomal location and duplication event analyses suggested that 26 MeLEAs distributed in 10 chromosomes and 11 MeLEA paralogues were subjected to purifying selection. Transcriptomic analysis showed the expression profiles of MeLEAs in different tissues of stem, leaves, and storage roots of three accessions. Comparative transcriptomic analysis revealed that the function of MeLEAs in response to drought may be differentiated in different accessions. Compared with the wild subspecies W14, more MeLEA genes were activated in cultivated varieties Arg7 and SC124 after drought treatment. Several MeLEA genes showed induction under various stresses and related signaling treatments. Taken together, this study demonstrates the transcriptional control of MeLEAs in tissue development and the responses to abiotic stress in cassava and identifies candidate genes for improving crop resistance to abiotic stress.

Genome-wide analysis of the NAC transcription factor family and their expression during the development and ripening of the Fragaria × ananassa fruits

NAC proteins are a family of transcription factors which have a variety of important regulatory roles in plants. They present a very well conserved group of NAC subdomains in the N-terminal region and a highly variable domain at the C-terminus. Currently, knowledge concerning NAC family in the strawberry plant remains very limited. In this work, we analyzed the NAC family of Fragaria vesca, and a total of 112 NAC proteins were identified after we curated the annotations from the version 4.0.a1 genome. They were placed into the ligation groups (pseudo-chromosomes) and described its physicochemical and genetic features. A microarray transcriptomic analysis showed six of them expressed during the development and ripening of the Fragaria x ananassa fruit. Their expression patterns were studied in fruit (receptacle and achenes) in different stages of development and in vegetative tissues. Also, the expression level under different hormonal treatments (auxins, ABA) and drought stress was investigated. In addition, they were clustered with other NAC transcription factor with known function related to growth and development, senescence, fruit ripening, stress response, and secondary cell wall and vascular development. Our results indicate that these six strawberry NAC proteins could play different important regulatory roles in the process of development and ripening of the fruit, providing the basis for further functional studies and the selection for NAC candidates suitable for biotechnological applications.

Genome-Wide Analyses of the NAC Transcription Factor Gene Family in Pepper (Capsicum annuum L.): Chromosome Location, Phylogeny, Structure, Expression Patterns, Cis-Elements in the Promoter, and Interaction Network

The NAM, ATAF1/2, and CUC2 (NAC) transcription factors form a large plant-specific gene family, which is involved in the regulation of tissue development in response to biotic and abiotic stress. To date, there have been no comprehensive studies investigating chromosomal location, gene structure, gene phylogeny, conserved motifs, or gene expression of NAC in pepper (Capsicum annuum L.). The recent release of the complete genome sequence of pepper allowed us to perform a genome-wide investigation of Capsicum annuum L. NAC (CaNAC) proteins. In the present study, a comprehensive analysis of the CaNAC gene family in pepper was performed, and a total of 104 CaNAC genes were identified. Genome mapping analysis revealed that CaNAC genes were enriched on four chromosomes (chromosomes 1, 2, 3, and 6). In addition, phylogenetic analysis of the NAC domains from pepper, potato, Arabidopsis, and rice showed that CaNAC genes could be clustered into three groups (I, II, and III). Group III, which contained 24 CaNAC genes, was exclusive to the Solanaceae plant family. Gene structure and protein motif analyses showed that these genes were relatively conserved within each subgroup. The number of introns in CaNAC genes varied from 0 to 8, with 83 (78.9%) of CaNAC genes containing two or less introns. Promoter analysis confirmed that CaNAC genes are involved in pepper growth, development, and biotic or abiotic stress responses. Further, the expression of 22 selected CaNAC genes in response to seven different biotic and abiotic stresses [salt, heat shock, drought, Phytophthora capsici, abscisic acid, salicylic acid (SA), and methyl jasmonate (MeJA)] was evaluated by quantitative RT-PCR to determine their stress-related expression patterns. Several putative stress-responsive CaNAC genes, including CaNAC72 and CaNAC27, which are orthologs of the known stress-responsive Arabidopsis gene ANAC055 and potato gene StNAC30, respectively, were highly regulated by treatment with different types of stress. Our results also showed that CaNAC36 plays an important role in the interaction network, interacting with 48 genes. Most of these genes are in the mitogen-activated protein kinase (MAPK) family. Taken together, our results provide a platform for further studies to identify the biological functions of CaNAC genes.

Comprehensive analysis of NAC transcription factors and their expression during fruit spine development in cucumber (Cucumis sativus L.)

The cucumber (Cucumis sativus L.) is an important vegetable crop worldwide, and fruit trichomes or spines are an important trait for external fruit quality. The mechanisms underlying spine formation are not well understood, but the plant-specific NAC family of transcription factors may play important roles in fruit spine initiation and development. In this study, we conducted a genome-wide survey and identified 91 NAC gene homologs in the cucumber genome. Clustering analysis classified these genes into six subfamilies; each contained a varying number of NAC family members with a similar intron-exon structure and conserved motifs. Quantitative real-time PCR analysis revealed tissue-specific expression patterns of these genes, including 10 and 12 that exhibited preferential expression in the stem and fruit, respectively. Thirteen of the 91 NAC genes showed higher expression in the wild-type plant than in its near-isogenic trichome mutant, suggesting their important roles in fruit spine development. Exogenous application of four plant hormones promoted spine formation and increased spine density on the cucumber fruits; several NAC genes showed differential expression over time in response to phytohormone treatments on cucumber fruit, implying their essential roles in fruit-trichome development. Among the NAC genes identified, 12 were found to be targets of 13 known cucumber micro-RNAs. Collectively, these findings provide a useful resource for further analysis of the interactions between NAC genes and genes underlying trichome organogenesis and development during fruit spine development in cucumber.

NAC Family Transcription Factors in Tobacco and Their Potential Role in Regulating Leaf Senescence

The NAC family is one of the largest families of plant-specific transcription factors (TFs) and NAC proteins play important regulatory roles in a variety of developmental and stress response processes in plants. Members of the NAC family TFs have been shown to be important regulators of leaf senescence in a number of plant species. Here we report the identification of the NAC family in tobacco (Nicotiana tabacum) and characterization of the potential role of some of the tobacco NAC TFs in regulating leaf senescence. A total of 154 NAC genes (NtNACs) were identified and clustered together with the Arabidopsis NAC family into fifteen groups (a-o). Transcriptome data analysis followed by qRT-PCR validation showed that the majority of the senescence-up-regulated NtNACs fall into subgroups NAC-b and f. A number of known senescence regulators from Arabidopsis also belong to these two subgroups. Among these senescence-up-regulated NtNACs, NtNAC080, a close homolog of AtNAP, is a positive regulator of leaf senescence. Overexpression of NtNAC080 caused early senescence in Arabidopsis leaves and NtNAC080 mutation induced by Cas9/gRNA in tobacco led to delayed leaf senescence.

Comparative transcriptome analysis of axillary buds in response to the shoot branching regulators gibberellin A3 and 6-benzyladenine in Jatropha curcas

Cytokinin (CK) is the primary hormone that positively regulates axillary bud outgrowth. However, in many woody plants, such as Jatropha curcas, gibberellin (GA) also promotes shoot branching. The molecular mechanisms underlying GA and CK interaction in the regulation of bud outgrowth in Jatropha remain unclear. To determine how young axillary buds respond to GA₃ and 6-benzyladenine (BA), we performed a comparative transcriptome analysis of the young axillary buds of Jatropha seedlings treated with GA₃ or BA. Two hundred and fifty genes were identified to be co-regulated in response to GA₃ or BA. Seven NAC family members were down-regulated after treatment with both GA₃ and BA, whereas these genes were up-regulated after treatment with the shoot branching inhibitor strigolactone. The expressions of the cell cycle genes CDC6, CDC45 and GRF5 were up-regulated after treatment with both GA₃ and BA, suggesting they may promote bud outgrowth via regulation of the cell cycle machinery. In the axillary buds, BA significantly increased the expression of GA biosynthesis genes JcGA20oxs and JcGA3ox1, and down-regulated the expression of GA degradation genes JcGA2oxs. Overall, the comprehensive transcriptome data set provides novel insight into the responses of young axillary buds to GA and CK.

The essence of NAC gene family to the cultivation of drought-resistant soybean (Glycine max L. Merr.) cultivars

BACKGROUND

The NAC gene family is notable due to its large size, as well as its relevance in crop cultivation - particularly in terms of enhancing stress tolerance of plants. These plant-specific proteins contain NAC domain(s) that are named after Petunia NAM and Arabidopsis ATAF1/2 and CUC2 transcription factors based on the consensus sequence they have. Despite the knowledge available regarding NAC protein function, an extensive study on the possible use of GmNACs in developing soybean cultivars with superior drought tolerance is yet to be done.

RESULTS

In response to this, our study was carried out, mainly through means of phylogenetic analysis (rice and Arabidopsis NAC genes served as seeding sequences). Through this, 139 GmNAC genes were identified and later grouped into 17 clusters. Furthermore, real-time quantitative PCR was carried out on drought-stressed and unstressed leaf tissues of both sensitive (B217 and H228) and tolerant (Jindou 74 and 78) cultivars. This was done to analyze the gene expression of 28 dehydration-responsive GmNAC genes. Upon completing the analysis, it was found that GmNAC gene expression is actually dependent on genotype. Eight of the 28 selected genes (GmNAC004, GmNAC021, GmNAC065, GmNAC066, GmNAC073, GmNAC082, GmNAC083 and GmNAC087) were discovered to have high expression levels in the drought-resistant soybean varieties tested. This holds true for both extreme and standard drought conditions. Alternatively, the drought-sensitive cultivars exhibited lower GmNAC expression levels in comparison to their tolerant counterparts.

CONCLUSION

The study allowed for the identification of eight GmNAC genes that could be focused upon in future attempts to develop superior soybean varieties, particularly in terms of drought resistance. This study revealed that there were more dehydration-responsive GmNAC genes as (GmNAC004, GmNAC005, GmNAC020 and GmNAC021) in addition to what were reported in earlier inquiries. It is important to note though, that discovering such notable genes is not the only goal of the study. It managed to put emphasis on the significance of further understanding the potential of soybean GmNAC genes, for the purpose of enhancing tolerance towards abiotic stress in general. This scientific inquiry has also revealed that cultivar genotypes tend to differ in their drought-induced gene expression.

MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network

Recent achievements in plant microRNA (miRNA), a large class of small and non-coding RNAs, are very exciting. A wide array of techniques involving forward genetic, molecular cloning, bioinformatic analysis, and the latest technology, deep sequencing have greatly advanced miRNA discovery. A tiny miRNA sequence has the ability to target single/multiple mRNA targets. Most of the miRNA targets are transcription factors (TFs) which have paramount importance in regulating the plant growth and development. Various families of TFs, which have regulated a range of regulatory networks, may assist plants to grow under normal and stress environmental conditions. This present review focuses on the regulatory relationships between miRNAs and different families of TFs like; NF-Y, MYB, AP2, TCP, WRKY, NAC, GRF, and SPL. For instance NF-Y play important role during drought tolerance and flower development, MYB are involved in signal transduction and biosynthesis of secondary metabolites, AP2 regulate the floral development and nodule formation, TCP direct leaf development and growth hormones signaling. WRKY have known roles in multiple stress tolerances, NAC regulate lateral root formation, GRF are involved in root growth, flower, and seed development, and SPL regulate plant transition from juvenile to adult. We also studied the relation between miRNAs and TFs by consolidating the research findings from different plant species which will help plant scientists in understanding the mechanism of action and interaction between these regulators in the plant growth and development under normal and stress environmental conditions.

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 characterization and analysis of bZIP transcription factor gene family related to abiotic stress in cassava

The basic leucine zipper (bZIP) transcription factor family plays crucial roles in various aspects of biological processes. Currently, no information is available regarding the bZIP family in the important tropical crop cassava. Herein, 77 bZIP genes were identified from cassava. Evolutionary analysis indicated that MebZIPs could be divided into 10 subfamilies, which was further supported by conserved motif and gene structure analyses. Global expression analysis suggested that MebZIPs showed similar or distinct expression patterns in different tissues between cultivated variety and wild subspecies. Transcriptome analysis of three cassava genotypes revealed that many MebZIP genes were activated by drought in the root of W14 subspecies, indicating the involvement of these genes in the strong resistance of cassava to drought. Expression analysis of selected MebZIP genes in response to osmotic, salt, cold, ABA, and H₂O₂ suggested that they might participate in distinct signaling pathways. Our systematic analysis of MebZIPs reveals constitutive, tissue-specific and abiotic stress-responsive candidate MebZIP genes for further functional characterization in planta, yields new insights into transcriptional regulation of MebZIP genes, and lays a foundation for understanding of bZIP-mediated abiotic stress response.