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

Genome-wide identification of NAC transcription factors in Acer paxii, and their expression dynamics during leaf aging

BACKGROUND

The NAC family (consisting of NAM, ATAF1/2, and CUC2) represents a crucial plant-specific transcription factor family, contributing significantly to various aspects of plant growth, development, and reactions to abiotic stresses. Yet, the underlying mechanism of NAC regulation in Acer paxii has not been discussed.

OBJECTIVES

Identification of NAC genes (ApNACs) in the genome of Acer paxii and exploring the regulatory network of NACs in mediating leaf senescence.

METHODS

A thorough genome-wide analysis of the NAC gene family in the Acer paxii genome was performed.

RESULTS

We identified 117 ApNACs from the Acer paxii genome database, which were irregularly distributed across 13 chromosomes. Phylogenetic analysis revealed that ApNAC genes were partitioned into 16 subgroups. There are four kinds of cis-regulatory elements in the promoter region of the ApNACs gene. We compared the expression levels of ApNAC genes at different times using transcriptome data and selected six ApNAC genes for qRT-PCR, which the results showed basic consistency with the transcriptome results. Six ApNACs and 187 TFs from different families were selected, and it was found that the TF families with the highest correlation were WRKY, MYB, bZIP, and ERF, and these TF families were reported to participate in the regulation function in senescence.

CONCLUSION

This study provides important data support for identifying the NAC gene family of Acer paxii and the regulatory function of the ApNAC genes on plant senescence, which will help to understand the NAC-mediated regulatory network in Acer paxii.

In silico genome-wide analysis of the growth-regulating factor gene family and their expression profiling in Vitis vinifera under biotic stress

Growth regulatory factors (GRFs) are transcription factors that encode the proteins involved in plant growth and development. However, no comprehensive analysis of Vitis vinifera GRF genes has yet been conducted. In the current study, we performed a genome-wide analysis of the GRF gene family to explore the VvGRF gene's role in Vitis vinifera. We identified 30 VvGRF genes in the Vitis vinifera genome, localized over 20 chromosomes. Based on evolutionary analysis, 49 GRF genes (nine AtGRF, ten FvGRF, and 30 VvGRF) were clustered into six groups. Many cis-elements involved in light control, defense, and plant growth have been identified in the promoter region of VvGRF genes, and multiple miRNAs have been predicted to be involved in regulating VvGRF gene expression. Protein-protein interaction analysis showed that nine VvGRF proteins formed a complex protein interaction network. Furthermore, the gene expression analysis of VvGRF revealed that VvGRF-5 and VvGRF-6 were highly upregulated suggesting that these genes are involved in biotic responses. This study provides comprehensive insights into the functional characteristics and occurrence of the VvGRF gene family in Vitis vinifera, which may be applied in breeding programs to enhance the growth of Vitis vinifera varieties under stress and growth changes.

Integration of Genetic and Imaging Data to Detect QTL for Root Traits in Interspecific Soybean Populations

Wild soybean, which has many desirable traits, such as adaptability to climate change-related stresses, is a valuable resource for expanding the narrow genetic diversity of cultivated soybeans. Plants require roots to adapt to different environments and optimize water and nutrient uptake to support growth and facilitate the storage of metabolites; however, it is challenging and costly to evaluate root traits under field conditions. Previous studies of quantitative trait loci (QTL) have been mainly based on cultivated soybean populations. In this study, an interspecific mapping population from a cross between wild soybean 'PI483463' and cultivar 'Hutcheson' was used to investigate QTLs associated with root traits using image data. Our results showed that 39 putative QTLs were distributed across 10 chromosomes (chr.). Seventeen of these were clustered in regions on chr. 8, 14, 15, 16, and 17, accounting for 19.92% of the phenotypic variation. We identified five significant QTL clusters influencing root-related traits, such as total root length, surface area, lateral total length, and number of tips, across five chr., with favorable alleles from both wild and cultivated soybeans. Furthermore, we identified eight candidate genes controlling these traits based on functional annotation. These genes were highly expressed in root tissues and directly or indirectly affected soybean root growth, development, and stress responses. Our results provide valuable insights for breeders aiming to optimize soybean root traits and leveraging genetic diversity from wild soybean species to develop varieties with improved root morphological traits, ultimately enhancing overall plant growth, productivity, and resilience.

Genome-wide analyses of the NAC transcription factor gene family in Acer palmatum provide valuable insights into the natural process of leaf senescence

Acer palmatum is a deciduous shrub or small tree. It is a popular ornamental plant because of its beautiful leaves, which change colour in autumn. This study revealed 116 ApNAC genes within the genome of A. palmatum. These genes are unevenly distributed on the 13 chromosomes of A. palmatum. An analysis of the phylogenetic tree of Arabidopsis thaliana NAC family members revealed that ApNAC proteins could be divided into 16 subgroups. A comparison of ApNAC proteins with NAC genes from other species suggested their potential involvement in evolutionary processes. Studies suggest that tandem and segmental duplications may be key drivers of the expansion of the ApNAC gene family. Analysis of the transcriptomic data and qRT‒PCR results revealed significant upregulation of most ApNAC genes during autumn leaf senescence compared with their expression levels in summer leaves. Coexpression network analysis revealed that the expression profiles of 10 ApNAC genes were significantly correlated with those of 200 other genes, most of which are involved in plant senescence processes. In conclusion, this study contributes to elucidating the theoretical foundation of the ApNAC gene family and provides a valuable basis for future investigations into the role of NAC genes in regulating leaf senescence in woody ornamental plants.

Genome-Wide Identification of the WRKY Gene Family in Four Cotton Varieties and the Positive Role of GhWRKY31 in Response to Salt and Drought Stress

The WRKY gene family is ubiquitously distributed in plants, serving crucial functions in stress responses. Nevertheless, the structural organization and evolutionary dynamics of WRKY genes in cotton have not been fully elucidated. In this study, a total of 112, 119, 217, and 222 WRKY genes were identified in Gossypium arboreum, Gossypium raimondii, Gossypium hirsutum, and Gossypium barbadense, respectively. These 670 WRKY genes were categorized into seven distinct subgroups and unequally distributed across chromosomes. Examination of conserved motifs, domains, cis-acting elements, and gene architecture collectively highlighted the evolutionary conservation and divergence within the WRKY gene family in cotton. Analysis of synteny and collinearity further confirmed instances of expansion, duplication, and loss events among WRKY genes during cotton evolution. Furthermore, GhWRKY31 transgenic Arabidopsis exhibited heightened germination rates and longer root lengths under drought and salt stress. Silencing GhWRKY31 in cotton led to reduced levels of ABA, proline, POD, and SOD, along with downregulated expression of stress-responsive genes. Yeast one-hybrid and molecular docking assays confirmed the binding capacity of GhWRKY31 to the W box of GhABF1, GhDREB2, and GhRD29. The findings collectively offer a systematic and comprehensive insight into the evolutionary patterns of cotton WRKYs, proposing a suitable regulatory framework for developing cotton cultivars with enhanced resilience to drought and salinity stress.

Genome-wide identification, molecular evolution and expression analysis of the B-box gene family in mung bean (Vigna radiata L.)

BACKGROUND

Mung bean (Vigna radiata L.) is an important warm-season grain legume. Adaptation to extreme environmental conditions, supported by evolution, makes mung bean a rich gene pool for stress tolerance traits. The exploration of resistance genes will provide important genetic resources and a theoretical basis for strengthening mung bean breeding. B-box (BBX) proteins play a major role in developmental processes and stress responses. However, the identification and analysis of the mung bean BBX gene family are still lacking.

RESULTS

In this study, 23 VrBBX genes were identified through comprehensive bioinformatics analysis and named based on their physical locations on chromosomes. All the VrBBXs were divided into five groups based on their phylogenetic relationships, the number of B-box they contained and whether there was an additional CONSTANS, CO-like and TOC1 (CCT) domain. Homology and collinearity analysis indicated that the BBX genes in mung bean and other species had undergone a relatively conservative evolution. Gene duplication analysis showed that only chromosomal segmental duplication contributed to the expansion of VrBBX genes and that most of the duplicated gene pairs experienced purifying selection pressure during evolution. Gene structure and motif analysis revealed that VrBBX genes clustered in the same group shared similar structural characteristics. An analysis of cis-acting elements indicated that elements related to stress and hormone responses were prevalent in the promoters of most VrBBXs. The RNA-seq data analysis and qRT-PCR of nine VrBBX genes demonstrated that VrBBX genes may play a role in response to environmental stress. Moreover, VrBBX5, VrBBX10 and VrBBX12 are important candidate genes for plant stress response.

CONCLUSIONS

In this study, we systematically analyzed the genomic characteristics and expression patterns of the BBX gene family under ABA, PEG and NaCl treatments. The results will help us better understand the complexity of the BBX gene family and provide valuable information for future functional characteristics of specific genes in this family.

Genome-Wide Identification of NAC Family Genes in Oat and Functional Characterization of AsNAC109 in Abiotic Stress Tolerance

The plant-specific NAC gene family is one of the largest transcription factor families, participating in plant growth regulation and stress response. Despite extensive characterization in various plants, our knowledge of the NAC family in oat is lacking. Herein, we identified 333 NAC genes from the latest release of the common oat genome. We provide a comprehensive overview of the oat NAC gene family, covering gene structure, chromosomal localization, phylogenetic characteristics, conserved motif compositions, and gene duplications. AsNAC gene expression in different tissues and the response to various abiotic stresses were characterized using RT-qPCR. The main driver of oat NAC gene family expansion was identified as segmental duplication using collinearity analysis. In addition, the functions of AsNAC109 in regulating abiotic stress tolerance in Arabidopsis were clarified. This is the first genome-wide investigation of the NAC gene family in cultivated oat, which provided a unique resource for subsequent research to elucidate the mechanisms responsible for oat stress tolerance and provides valuable clues for the improvement of stress resistance in cultivated oat.

Novel loci associated with resistance to downy and powdery mildew in grapevine

Among the main challenges in current viticulture, there is the increasing demand for sustainability in the protection from fungal diseases, such as downy mildew (DM) and powdery mildew (PM). Breeding disease-resistant grapevine varieties is a key strategy for better managing fungicide inputs. This study explores the diversity of grapevine germplasm (cultivated and wild) from Caucasus and neighboring areas to identify genotypes resistant to DM and PM, based on 13 Simple Sequence Repeat (SSR) loci and phenotypical (artificial pathogen inoculation) analysis, and to identify loci associated with DM and PM resistance, via Genome-Wide Association Analysis (GWAS) on Single Nucleotide Polymorphism (SNP) profiles. SSR analysis revealed resistant alleles for 16 out of 88 genotypes. Phenotypic data identified seven DM and 31 PM resistant genotypes. GWAS identified two new loci associated with DM resistance, located on chromosome 15 and 16 (designated as Rpv36 and Rpv37), and two with PM resistance, located on chromosome 6 and 17 (designated as Ren14 and Ren15). The four novel loci identified genomic regions rich in genes related to biotic stress response, such as genes involved in pathogen recognition, signal transduction and resistance response. This study highlights potential candidate genes associated with resistance to DM and PM, providing valuable insights for breeding programs for resistant varieties. To optimize their utilization, further functional characterization studies are recommended.

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

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

Genome-wide identification and characterization of NAC transcription factor family members in Trifolium pratense and expression analysis under lead stress

BACKGROUND

The NAC TF family is widely involved in plant responses to various types of stress. Red clover (Trifolium pratense) is a high-quality legume, and the study of NAC genes in red clover has not been comprehensive. The aim of this study was to analyze the NAC gene family of red clover at the whole-genome level and explore its potential role in the Pb stress response.

RESULTS

In this study, 72 TpNAC genes were identified from red clover; collinearity analysis showed that there were 5 pairs of large fragment replicators of TpNAC genes, and red clover was found to be closely related to Medicago truncatula. Interestingly, the TpNAC genes have more homologs in Arabidopsis thaliana than in soybean (Glycine max). There are many elements in the TpNAC genes promoters that respond to stress. Gene expression analysis showed that all the TpNAC genes responded to Pb stress. qRT-PCR showed that the expression levels of TpNAC29 and TpNAC42 were significantly decreased after Pb stress. Protein interaction network analysis showed that 21 TpNACs and 23 other genes participated in the interaction. In addition, the TpNAC proteins had three possible 3D structures, and the secondary structure of these proteins were mainly of other types. These results indicated that most TpNAC members were involved in the regulation of Pb stress in red clover.

CONCLUSION

These results suggest that most TpNAC members are involved in the regulation of Pb stress in red clover. TpNAC members play an important role in the response of red clover to Pb stress.

New strategies to address world food security and elimination of malnutrition: future role of coarse cereals in human health

As we face increasing challenges of world food security and malnutrition, coarse cereals are coming into favor as an important supplement to human staple foods due to their high nutritional value. In addition, their functional components, such as flavonoids and polyphenols, make them an important food source for healthy diets. However, we lack a systematic understanding of the importance of coarse cereals for world food security and nutritional goals. This review summarizes the worldwide cultivation and distribution of coarse cereals, indicating that the global area for coarse cereal cultivation is steadily increasing. This paper also focuses on the special adaptive mechanisms of coarse cereals to drought and discusses the strategies to improve coarse cereal crop yields from the perspective of agricultural production systems. The future possibilities, challenges, and opportunities for coarse cereal production are summarized in the face of food security challenges, and new ideas for world coarse cereal production are suggested.

Genome-Wide Analysis and Identification of 1-Aminocyclopropane-1-Carboxylate Synthase (ACS) Gene Family in Wheat (Triticum aestivum L.)

Ethylene has an important role in regulating plant growth and development as well as responding to adversity stresses. The 1-aminocyclopropane-1-carboxylate synthase (ACS) is the rate-limiting enzyme for ethylene biosynthesis. However, the role of the ACS gene family in wheat has not been examined. In this study, we identified 12 ACS members in wheat. According to their position on the chromosome, we named them TaACS1-TaACS12, which were divided into four subfamilies, and members of the same subfamilies had similar gene structures and protein-conserved motifs. Evolutionary analysis showed that fragment replication was the main reason for the expansion of the TaACS gene family. The spatiotemporal expression specificity showed that most of the members had the highest expression in roots, and all ACS genes contained W box elements that were related to root development, which suggested that the ACS gene family might play an important role in root development. The results of the gene expression profile analysis under stress showed that ACS members could respond to a variety of stresses. Protein interaction prediction showed that there were four types of proteins that could interact with TaACS. We also obtained the targeting relationship between TaACS family members and miRNA. These results provided valuable information for determining the function of the wheat ACS gene, especially under stress.

Comprehensive analysis of NAC transcription factors in Scutellaria baicalensis and their response to exogenous ABA and GA3

The NAC is a plant-specific family of transcription factor that plays important roles in various biological processes. Scutellaria baicalensis Georgi, belongs to the Lamiaceae family and has been widely used as a traditional herb with a wide range of pharmacological activities, including antitumor, heat-clearing, and detoxifying functions. However, no study on the NAC family in S. baicalensis has been conducted to date. In the present study, we identified 56 SbNAC genes using genomic and transcriptome analyses. These 56 SbNACs were unevenly distributed across nine chromosomes and were phylogenetically divided into six clusters. Cis-element analysis identified plant growth and development-, phytohormone-, light-, and stress-responsive elements were present in SbNAC genes promoter regions. Protein-protein interaction analysis was performed using Arabidopsis homologous proteins. Potential transcription factors, including bHLH, ERF, MYB, WRKY, and bZIP, were identified and constructed a regulatory network with SbNAC genes. The expression of 12 flavonoid biosynthetic genes was significantly upregulated with abscisic acid (ABA) and gibberellin (GA₃) treatments. Eight SbNAC genes (SbNAC9/32/33/40/42/43/48/50) also exhibited notable variation with two phytohormone treatments, among which SbNAC9 and SbNAC43 showed the most significant variation and deserved further study. Additionally, SbNAC44 displayed a positive correlation with C4H3, PAL5, OMT3, and OMT6, while SbNAC25 had negatively correlated with OMT2, CHI, F6H2, and FNSII-2. This study constitutes the first analysis of SbNAC genes and lays the basis foundation for further functional studies of SbNAC genes family members, while it may also facilitate the genetic improvement of plants and breeding of elite S. baicalensis varieties.

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 Investigation of the NAC Transcription Factor Family in Apocynum venetum Revealed Their Synergistic Roles in Abiotic Stress Response and Trehalose Metabolism

NAC (NAM, ATAF1/2, and CUC2) transcription factors (TFs) are one of the most prominent plant-specific TF families and play essential roles in plant growth, development and adaptation to abiotic stress. Although the NAC gene family has been extensively characterized in many species, systematic analysis is still relatively lacking in Apocynum venetum (A. venetum). In this study, 74 AvNAC proteins were identified from the A. venetum genome and were classified into 16 subgroups. This classification was consistently supported by their gene structures, conserved motifs and subcellular localizations. Nucleotide substitution analysis (Ka/Ks) showed the AvNACs to be under the influence of strong purifying selection, and segmental duplication events were found to play the dominant roles in the AvNAC TF family expansion. Cis-elements analysis demonstrated that the light-, stress-, and phytohormone-responsive elements being dominant in the AvNAC promoters, and potential TFs including Dof, BBR-BPC, ERF and MIKC_MADS were visualized in the TF regulatory network. Among these AvNACs, AvNAC58 and AvNAC69 exhibited significant differential expression in response to drought and salt stresses. The protein interaction prediction further confirmed their potential roles in the trehalose metabolism pathway with respect to drought and salt resistance. This study provides a reference for further understanding the functional characteristics of NAC genes in the stress-response mechanism and development of A. venetum.

Transcriptomic and biochemical analyses of drought response mechanism in mung bean (Vignaradiata (L.) Wilczek) leaves

Drought is a major factor that limiting mung bean development. To clarify the molecular mechanism of mung bean in response to drought stress, 2 mung bean groups were established, the experimental group (drought-treated) and the control group (normal water management). With prominent difference of 2 groups in stomatal conductance, relative water content and phenotype, leaf samples were collected at 4 stages, and the physiological index of MDA, POD, chlorophyll, and soluble proteins were estimated. RNA-seq was used to obtain high quality data of samples, and differentially expressed genes were identified by DESeq2. With GO and KEGG analysis, DEGs were enriched into different classifications and pathways. WGCNA was used to detect the relationship between physiological traits and genes, and qPCR was performed to confirm the accuracy of the data. We obtained 169.49 Gb of clean data from 24 samples, and the Q30 of each date all exceeded 94%. In total, 8963 DEGs were identified at 4 stages between the control and treated samples, and the DEGs were involved in most biological processes. 1270 TFs screened from DEGs were clustered into 158 TF families, such as AP2, RLK-Pelle-DLSVA, and NAC TF families. Genes related to physiological traits were closely related to plant hormone signaling, carotenoid biosynthesis, chlorophyll metabolism, and protein processing. This paper provides a large amount of data for drought research in mung bean.

Transcriptome-wide identification of NAC (no apical meristem/Arabidopsis transcription activation factor/cup-shaped cotyledon) transcription factors potentially involved in salt stress response in garlic

Soil salinity has been an increasing problem worldwide endangering crop production and human food security. It is an ideal strategy to excavate stress resistant genes and develop salt tolerant crops. NAC (no apical meristem/Arabidopsis transcription activation factor/cup-shaped cotyledon) transcription factors have been demonstrated to be involved in salt stress response. However, relevant studies have not been observed in garlic, an important vegetable consumed in the world. In this study, a total of 46 AsNAC genes encoding NAC proteins were identified in garlic plant by transcriptome data. Phylogenetic analysis showed that the examined AsNAC proteins were clustered into 14 subgroups. Motif discovery revealed that the conserved domain region was mainly composed of five conserved subdomains. Most of the genes selected could be induced by salt stress in different tissues, indicating a potential role in salt stress response. Further studies may focus on the molecular mechanisms of the AsNAC genes in salt stress response. The results of the current work provided valuable resources for researchers aimed at developing salt tolerant crops.

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.