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

Genome-wide investigation and analysis of NAC transcription factor family in Populus tomentosa and expression analysis under salt stress

The NAC transcription factor family is one of the largest families of TFs in plants, and members of NAC gene family play important roles in plant growth and stress response. Recent release of the haplotype-resolved genome assembly of P. tomentosa provide a platform for NAC protein genome-wide analysis. A total of 270 NAC genes were identified and a comprehensive overview of the PtoNAC gene family is presented, including gene promoter, structure and conserved motif analyses, chromosome localization and collinearity analysis, protein phylogeny, expression pattern, and interaction analysis. The results indicate that protein length, molecular weight, and theoretical isoelectric points of the NAC TF family vary, while gene structure and motif are relatively conserved. Chromosome mapping analysis showed that the P. tomentosa NAC genes are unevenly distributed on 19 chromosomes. The interchromosomal evolutionary results indicate 12 pairs of tandem and 280 segmental duplications. Segmental duplication is possibly related to amplification of P. tomentosa NAC gene family. Expression patterns of 35 PtoNAC genes from P. tomentosa subgroup were analysed under high salinity, and seven NAC genes were induced by this treatment. Promoter and protein interaction network analyses showed that PtoNAC genes are closely associated with growth, development, and abiotic and biotic stress, especially salt stress. These results provide a meaningful reference for follow-up studies of the functional characteristics of NAC genes in the mechanism of stress response and their potential roles in development of P. tomentosa.

Genome-wide analysis of AP2/ERF gene and functional analysis of CqERF24 gene in drought stress in quinoa

Quinoa is a crop with high nutritional value and strong stress resistance. AP2/ERF transcription factors play a key role in plant growth and development. In this study, 148 AP2/ERF genes were identified in quinoa, which were divided into 5 subfamilies, including ERF, AP2, DREB, RAV and Soloist. The results showed that the number of introns ranged from 0 to 11, and the Motif 1-Motif 4 was highly conserved in most CqAP2/ERF proteins. The 148 CqAP2/ERF genes were distributed on 19 chromosomes. There were 93 pairs of duplicating genes in this family, and gene duplication played a critical role in the expansion of this family. Protein-protein interaction indicated that the proteins in CqAP2/ERF subfamily exhibited complex interactions, and GO enrichment analysis indicated that 148 CqAP2/ERF proteins were involved in transcription factor activity. In addition, CqAP2/ERF gene contains a large number of elements related to hormones in promoter region (IAA, GA, SA, ABA and MeJA) and stresses (salt, drought, low temperature and anaerobic induction). Transcriptome analysis under drought stress indicated that most of the CqAP2/ERF genes were responsive to drought stress, and subcellular localization indicated that CqERF24 was location in the nucleus, qRT-PCR results also showed that most of the genes such as CqERF15, CqERF24, CqDREB03, CqDREB14, CqDREB37 and CqDREB43 also responded to drought stress in roots and leaves. Overexpression of CqERF24 in Arabidopsis thaliana enhanced drought resistance by increasing antioxidant enzyme activity and activation-related stress genes, and the gene is sensitive to ABA, while silencing CqERF24 in quinoa decreased drought tolerance. In addition, overexpression of CqERF24 in quinoa calli enhanced resistance to mannitol. These results lay a solid foundation for further study on the role of AP2/ERF family genes in quinoa under drought stress.

Molecular Characterization and Expression Analysis of YABBY Genes in Chenopodium quinoa

Plant-specific YABBY transcription factors play an important role in lateral organ development and abiotic stress responses. However, the functions of the YABBY genes in quinoa remain elusive. In this study, twelve YABBY (CqYAB) genes were identified in the quinoa genome, and they were distributed on nine chromosomes. They were classified into FIL/YAB3, YAB2, YAB5, INO, and CRC clades. All CqYAB genes consist of six or seven exons, and their proteins contain both N-terminal C2C2 zinc finger motifs and C-terminal YABBY domains. Ninety-three cis-regulatory elements were revealed in CqYAB gene promoters, and they were divided into six groups, such as cis-elements involved in light response, hormone response, development, and stress response. Six CqYAB genes were significantly upregulated by salt stress, while one was downregulated. Nine CqYAB genes were upregulated under drought stress, whereas six CqYAB genes were downregulated under cadmium treatment. Tissue expression profiles showed that nine CqYAB genes were expressed in seedlings, leaves, and flowers, seven in seeds, and two specifically in flowers, but no CqYAB expression was detected in roots. Furthermore, CqYAB4 could rescue the ino mutant phenotype in Arabidopsis but not CqYAB10, a paralog of CqYAB4, indicative of functional conservation and divergence among these YABBY genes. Taken together, these results lay a foundation for further functional analysis of CqYAB genes in quinoa growth, development, and abiotic stress responses.

Identification of the NAC Transcription Factor Family during Early Seed Development in Akebia trifoliata (Thunb.) Koidz

This study aimed to gain an understanding of the possible function of NACs by examining their physicochemical properties, structure, chromosomal location, and expression. Being a family of plant-specific transcription factors, NAC (petunia no apical meristem and Arabidopsis thaliana ATAF1, ATAF2, and CUC2) is involved in plant growth and development. None of the NAC genes has been reported in Akebia trifoliata (Thunb.) Koidz (A. trifoliata). In this study, we identified 101 NAC proteins (AktNACs) in the A. trifoliata genome by bioinformatic analysis. One hundred one AktNACs were classified into the following twelve categories based on the phylogenetic analysis of NAC protein: NAC-a, NAC-b, NAC-c, NAC-d, NAC-e, NAC-f, NAC-g, NAC-h, NAC-i, NAC-j, NAC-k, and NAC-l. The accuracy of the clustering results was demonstrated based on the gene structure and conserved motif analysis of AktNACs. In addition, we identified 44 pairs of duplication genes, confirming the importance of purifying selection in the evolution of AktNACs. The morphology and microstructure of early A. trifoliata seed development showed that it mainly underwent rapid cell division, seed enlargement, embryo formation and endosperm development. We constructed AktNACs co-expression network and metabolite correlation network based on transcriptomic and metabolomic data of A. trifoliata seeds. The results of the co-expression network showed that 25 AtNAC genes were co-expressed with 233 transcription factors. Metabolite correlation analysis showed that 23 AktNACs were highly correlated with 28 upregulated metabolites. Additionally, 25 AktNACs and 235 transcription factors formed co-expression networks with 141 metabolites, based on correlation analysis involving AktNACs, transcription factors, and metabolites. Notably, AktNAC095 participates in the synthesis of 35 distinct metabolites. Eight of these metabolites, strongly correlated with AktNAC095, were upregulated during early seed development. These studies may provide insight into the evolution, possible function, and expression of AktNACs genes.

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.

Genome-wide identification, molecular evolution and expression analysis of the non-specific lipid transfer protein (nsLTP) family in Setaria italica

BACKGROUND

Foxtail millet (Setaria italica L.) is a millet species with high tolerance to stressful environments. Plant non-specific lipid transfer proteins (nsLTPs) are a kind of small, basic proteins involved in many biological processes. So far, the genome of S. italica has been fully sequenced, and a comprehensive understanding of the evolution and expression of the nsLTP family is still lacking in foxtail millet.

RESULTS

Forty-five nsLTP genes were identified in S. italica and clustered into 5 subfamilies except three single genes (SinsLTP38, SinsLTP7, and SinsLTP44). The proportion of SinsLTPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures. Besides, 5 SinsLTP duplication events were investigated. Both tandem and segmental duplication contributed to nsLTP expansion in S. italica, and the duplicated SinsLTPs had mainly undergone purifying selection pressure, which suggested that the function of the duplicated SinsLTPs might not diverge much. Moreover, we identified the nsLTP members in 5 other monocots, and 41, 13, 10, 4, and 1 orthologous gene pairs were identified between S. italica and S. viridis, S. bicolor, Z. mays, O. sativa, and B. distachyon, respectively. The functional divergence within the nsLTP orthologous genes might be limited. In addition, the tissue-specific expression patterns of the SinsLTPs were investigated, and the expression profiles of the SinsLTPs in response to abiotic stress were analyzed, all the 10 selected SinsLTPs were responsive to drought, salt, and cold stress. Among the selected SinsLTPs, 2 paired duplicated genes shared almost equivalent expression profiles, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, conserved domain and gene structure, expansion pattern, and expression profile of the nsLTP family in S. italica. These findings could pave a way for further comparative genomic and evolution analysis of nsLTP family in foxtail millet and related monocots, and lay the foundation for the functional analysis of the nsLTPs in S. italica.

Comprehensive analysis of the NAC transcription factor gene family in Kandelia obovata reveals potential members related to chilling tolerance

BACKGROUND

Kandelia obovata is an important mangrove species extensively distributed in Eastern Asia that is susceptible to low-temperature stress. NAC (NAM, ATAF1/2 and CUC2) domain proteins are transcription factors (TFs) that play various roles in plant growth and development and in the plant response to environmental stresses. Nevertheless, genome-wide analyses of K. obovata NAC genes (KoNACs) and their responses to chilling stress have rarely been studied.

METHODS

The KoNAC gene family was identified and characterized using bioinformatic analysis, the subcellular location of some NAC proteins was confirmed using confocal microscopy analysis, and the KoNACs that responded to chilling stress were screened using RNA-seq and qRT-PCR analysis.

RESULTS

A total of 79 KoNACs were identified, and they were unequally distributed across all 18 chromosomes of K. obovata. The KoNAC proteins could be divided into 16 subgroups according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana. The KoNACs exhibited greater synteny with A. thaliana sequences than with Oryza sativa sequences, indicating that KoNACs underwent extensive evolution after the divergence of dicotyledons and monocotyledons. Segmental duplication was the main driving force of the expansions of KoNAC genes. Confocal microscopy analysis verified that the four randomly selected KoNACs localized to the nucleus, indicating the accuracy of the bioinformatic predictions. Tissue expression pattern analysis demonstrated that some KoNAC genes showed tissue-specific expression, suggesting that these KoNACs might be important for plant development and growth. Additionally, the expression levels of 19 KoNACs were significantly (15 positively and 4 negatively) induced by cold treatment, demonstrating that these KoNACs might play important roles during cold stress responses and might be candidate genes for the genetic engineering of K. obovata with enhanced chilling stress tolerance. Coexpression network analysis revealed that 381 coexpressed pairs (between 13 KoNACs and 284 other genes) were significantly correlated.

CONCLUSIONS

Seventy-nine KoNACs were identified in K. obovata, nineteen of which displayed chilling-induced expression patterns. These genes may serve as candidates for functional analyses of KoNACs engaged in chilling stress. Our results lay the foundation for evolutionary analyses of KoNACs and their molecular mechanisms in response to environmental stress.

Comprehensive Analysis of StSRO Gene Family and Its Expression in Response to Different Abiotic Stresses in Potato

As a highly conserved family of plant-specific proteins, SIMILAR-TO-RCD-ONE (SROs) play an essential role in plant growth, development and response to abiotic stresses. In this study, six StSRO genes were identified by searching the PARP, RST and WWE domains based on the genome-wide data of potato database DM v6.1, and they were named StSRO1–6 according to their locations on chromosomes. StSRO genes were comprehensively analyzed using bioinformatics methods. The results showed that six StSRO genes were irregularly distributed on five chromosomes. Phylogenetic analysis showed that 30 SRO genes of four species were distributed in three groups, while StSRO genes were distributed in groups II and III. The promoter sequence of StSRO genes contained many cis-acting elements related to hormones and stress responses. In addition, the expression level of StSRO genes in different tissues of doubled monoploid (DM) potato, as well as under salt, drought stresses and hormone treatments, was analyzed by RNA-seq data from the online database and quantitative real-time polymerase chain reaction (qRT-PCR) analysis. Furthermore, the expression level of StSRO genes was analyzed by transcriptome analysis under mild, moderate and severe salt stress. It was concluded that StSRO genes could respond to different abiotic conditions, but their expression level was significantly different. This study lays a foundation for further studies on the biological functions of the StSRO gene family.

Genome-Wide Identification and Characterization of SPL Family Genes in Chenopodium quinoa

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes encode a large family of plant-specific transcription factors that play important roles in plant growth, development, and stress responses. However, there is little information available on SPL genes in Chenopodiaceae. Here, 23 SPL genes were identified and characterized in the highly nutritious crop Chenopodium quinoa. Chromosome localization analysis indicated that the 23 CqSPL genes were unevenly distributed on 12 of 18 chromosomes. Two zinc finger-like structures and a nuclear location signal were present in the SBP domains of all CqSPLs, with the exception of CqSPL21/22. Phylogenetic analysis revealed that these genes were classified into eight groups (group I–VIII). The exon–intron structure and motif composition of the genes in each group were similar. Of the 23 CqSPLs, 13 were potential targets of miR156/7. In addition, 5 putative miR156-encoding loci and 13 putative miR157-encoding loci were predicted in the quinoa genome, and they were unevenly distributed on chromosome 1–4. The expression of several Cqu-MIR156/7 loci was confirmed by reverse transcription polymerase chain reaction in seedlings. Many putative cis-elements associated with light, stress, and phytohormone responses were identified in the promoter regions of CqSPLs, suggesting that CqSPL genes are likely involved in the regulation of key developmental processes and stress responses. Expression analysis revealed highly diverse expression patterns of CqSPLs among tissues. Many CqSPLs were highly expressed in leaves, flowers, and seeds, and their expression levels were low in the roots, suggesting that CqSPLs play distinct roles in the development and growth of quinoa. The expression of 13 of 23 CqSPL genes responded to salt treatment (11 up-regulated and 2 down-regulated). A total of 22 of 23 CqSPL genes responded to drought stress (21 up-regulated and 1 down-regulated). Moreover, the expression of 14 CqSPL genes was significantly altered following cadmium treatment (3 up-regulated and 11 down-regulated). CqSPL genes are thus involved in quinoa responses to salt/drought and cadmium stresses. These findings provide new insights that will aid future studies of the biological functions of CqSPLs in C. quinoa.

HuNAC20 and HuNAC25, Two Novel NAC Genes from Pitaya, Confer Cold Tolerance in Transgenic Arabidopsis

NAC transcription factors are one of the largest families of transcriptional regulators in plants, and members of the gene family play vital roles in regulating plant growth and development processes including biotic/abiotic stress responses. However, little information is available about the NAC family in pitaya. In this study, we conducted a genome-wide analysis and a total of 64 NACs (named HuNAC1-HuNAC64) were identified in pitaya (Hylocereus). These genes were grouped into fifteen subgroups with diversities in gene proportions, exon–intron structures, and conserved motifs. Genome mapping analysis revealed that HuNAC genes were unevenly scattered on all eleven chromosomes. Synteny analysis indicated that the segmental duplication events played key roles in the expansion of the pitaya NAC gene family. Expression levels of these HuNAC genes were analyzed under cold treatments using qRT-PCR. Four HuNAC genes, i.e., HuNAC7, HuNAC20, HuNAC25, and HuNAC30, were highly induced by cold stress. HuNAC7, HuNAC20, HuNAC25, and HuNAC30 were localized exclusively in the nucleus. HuNAC20, HuNAC25, and HuNAC30 were transcriptional activators while HuNAC7 was a transcriptional repressor. Overexpression of HuNAC20 and HuNAC25 in Arabidopsis thaliana significantly enhanced tolerance to cold stress through decreasing ion leakage, malondialdehyde (MDA), and H₂O₂ and O₂⁻ accumulation, accompanied by upregulating the expression of cold-responsive genes (AtRD29A, AtCOR15A, AtCOR47, and AtKIN1). This study presents comprehensive information on the understanding of the NAC gene family and provides candidate genes to breed new pitaya cultivars with tolerance to cold conditions through genetic transformation.

Genome-Wide Identification of NAC Transcription Factor Family in Juglans mandshurica and Their Expression Analysis during the Fruit Development and Ripening

The NAC (NAM, ATAF and CUC) gene family plays a crucial role in the transcriptional regulation of various biological processes and has been identified and characterized in multiple plant species. However, genome-wide identification of this gene family has not been implemented in Juglans mandshurica, and specific functions of these genes in the development of fruits remain unknown. In this study, we performed genome-wide identification and functional analysis of the NAC gene family during fruit development and identified a total of 114 JmNAC genes in the J. mandshurica genome. Chromosomal location analysis revealed that JmNAC genes were unevenly distributed in 16 chromosomes; the highest numbers were found in chromosomes 2 and 4. Furthermore, according to the homologues of JmNAC genes in Arabidopsis thaliana, a phylogenetic tree was constructed, and the results demonstrated 114 JmNAC genes, which were divided into eight subgroups. Four JmNAC gene pairs were identified as the result of tandem duplicates. Tissue-specific analysis of JmNAC genes during different developmental stages revealed that 39 and 25 JmNAC genes exhibited upregulation during the mature stage in walnut exocarp and embryos, indicating that they may serve key functions in fruit development. Furthermore, 12 upregulated JmNAC genes were common in fruit ripening stage in walnut exocarp and embryos, which demonstrated that these genes were positively correlated with fruit development in J. mandshurica. This study provides new insights into the regulatory functions of JmNAC genes during fruit development in J. mandshurica, thereby improving the understanding of characteristics and evolution of the JmNAC gene family.

Characterization of NAC family genes in Salvia miltiorrhiza and NAC2 potentially involved in the biosynthesis of tanshinones

The NAC (NAM, ATAF, and CUC) family members are specific transcription factors in plants. The large family is involved in many plant growth and developmental processes, as well as in abiotic/biotic stress responses. It has been well studied in the genomes of various plants, including Arabidopsis thaliana, tomato, and quinoa. However, identification and functional studies of NAC family members in medicinal Salvia miltiorrhiza are limited. Here, we systematically identified 84 NAC genes and named them according to their gene IDs in the recently sequenced genome. The phylogeny of NAC family protein sequences was analyzed using bioinformatics methods, which divided them into nine subfamilies. Then, their chromosomal locations, gene structures and conserved domains were analyzed comprehensively. To further investigate the regulatory functions of NACs in S. miltiorrhiza, we analyzed the response of 10 selected NAC genes to methyl jasmonate and used NAC2 for transgenic experiments. The overexpression of Sm-NAC2 decreased the tanshinone I and IIA contents by 56% and 62%, respectively. However, Sm-NAC2-RNAi promoted the accumulation of four tanshinones, tanshinone I, tanshinone IIA, cryptotanshinone, and dihydrotanshinone I, which increased 3.68-, 4.1-, 3.13- and 5.9- fold, respectively, compared with wild type. In the tanshinone biosynthetic pathways, the overexpression of Sm-NAC2 down-regulated CYP76AH1, and the silencing of Sm-NAC2 up-regulated the expression levels of HMGR1, DXS2, KSL2, and CYP76AH1. This study provides information on the evolution of Sm-NAC genes and their possible functions, and it lays a foundation for further research into the NAC family-associated regulation of tanshinone biosynthesis.

Comprehensive multiomics analysis reveals key roles of NACs in plant growth and development and its environmental adaption mechanism by regulating metabolite pathways

Abnormal environmental conditions induce polyploidization and exacerbate vulnerability to agricultural production. Polyploidization is a pivotal event for plant adaption to stress and the expansion of transcription factors. NACs play key roles in plant stress resistance and growth and development, but the adaptive mechanism of NACs during plant polyploidization remain to be explored. Here, we identified and analyzed NACs from 15 species and found that the expansion of NACs was contributed by polyploidization. The regulatory networks were systematically analyzed based on polyomics. NACs might influence plant phenotypes and were correlated with amino acids acting as nitrogen source, indicating that NACs play a vital role in plant development. More importantly, in quinoa and Arabidopsis thaliana, NACs enabled plants to resist stress by regulating flavonoid pathways, and the universality was further confirmed by the Arabidopsis population. Our study provides a cornerstone for future research into improvement of important agronomic traits by transcription factors in a changing global environment.

Genome-wide survey, characterization, and expression analysis of bZIP transcription factors in Chenopodium quinoa

BACKGROUND

Chenopodium quinoa Willd. (quinoa) is a pseudocereal crop of the Amaranthaceae family and represents a promising species with the nutritional content and high tolerance to stressful environments, such as soils affected by high salinity. The basic leucine zipper (bZIP) transcription factor represents exclusively in eukaryotes and can be related to many biological processes. So far, the genomes of quinoa and 3 other Amaranthaceae crops (Spinacia oleracea, Beta vulgaris, and Amaranthus hypochondriacus) have been fully sequenced. However, information about the bZIPs in these Amaranthaceae species is limited, and genome-wide analysis of the bZIP family is lacking in quinoa.

RESULTS

We identified 94 bZIPs in quinoa (named as CqbZIP1-CqbZIP94). All the CqbZIPs were phylogenetically splitted into 12 distinct subfamilies. The proportion of CqbZIPs was different in each subfamily, and members within the same subgroup shared conserved exon-intron structures and protein motifs. Besides, 32 duplicated CqbZIP gene pairs were investigated, and the duplicated CqbZIPs had mainly undergone purifying selection pressure, which suggested that the functions of the duplicated CqbZIPs might not diverge much. Moreover, we identified the bZIP members in 3 other Amaranthaceae species, and 41, 32, and 16 orthologous gene pairs were identified between quinoa and S. oleracea, B. vulgaris, and A. hypochondriacus, respectively. Among them, most were a single copy being present in S. oleracea, B. vulgaris, and A. hypochondriacus, and two copies being present in allotetraploid quinoa. The function divergence within the bZIP orthologous genes might be limited. Additionally, 11 selected CqbZIPs had specific spatial expression patterns, and 6 of 11 CqbZIPs were up-regulated in response to salt stress. Among the selected CqbZIPs, 3 of 4 duplicated gene pairs shared similar expression patterns, suggesting that these duplicated genes might retain some essential functions during subsequent evolution.

CONCLUSIONS

The present study provided the first systematic analysis for the phylogenetic classification, motif and gene structure, expansion pattern, and expression profile of the bZIP family in quinoa. Our results would lay an important foundation for functional and evolutionary analysis of CqbZIPs, and provide promising candidate genes for further investigation in tissue specificity and their functional involvement in quinoa's resistance to salt stress.