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

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

Identification and functional analysis of the DOF gene family in Populus simonii: implications for development and stress response

Background

Populus simonii, a notable native tree species in northern China, demonstrates impressive resistance to stress, broad adaptability, and exceptional hybridization potential. DOF family is a class of specific transcription factors that only exist in plants, widely participating in plant growth and development, and also playing an important role in abiotic stress response. To date, there have been no reported studies on the DOF gene family in P. simonii, and the expression levels of this gene family in different tissues of poplar, as well as its expression patterns under cold, heat, and other stress conditions, remain unclear.

Methods

In this study, DOF gene family were identified from the P. simonii genome, and various bioinformatics data on the DOF gene family, gene structure, gene distribution, promoters and regulatory networks were analyzed. Quantitative real time PCR technology was used to verify the expression patterns of the DOF gene family in different poplar tissues.

Results

This research initially pinpointed 41 PSDOF genes in P. simonii genome. The chromosomal localization results revealed that the PSDOF genes is unevenly distributed among 19 chromosomes, with the highest number of genes located on chromosomes 4, 5, and 11. A phylogenetic tree was constructed based on the homology between Arabidopsis thaliana and P. simonii, dividing the 41 PSDOF genes into seven subgroups. The expression patterns of PSDOF genes indicated that specific genes are consistently upregulated in various tissues and under different stress conditions, suggesting their pivotal involvement in both plant development and response to stress. Notably, PSDOF35 and PSDOF28 serve as pivotal hubs in the interaction network, playing a unique role in coordinating with other genes within the family.

Conclusion

The analysis enhances our comprehension of the functions of the DOF gene family in tissue development and stress responses within P. simonii. These findings provide a foundation for future exploration into the biological roles of DOF gene family.

Genome-wide identification and expression pattern analysis of the Aux/IAA (auxin/indole-3-acetic acid) gene family in alfalfa (Medicago sativa) and the potential functions under drought stress

BACKGROUND

Auxin/induced-3-acetic acid (Aux/IAA) is an important plant hormone that affects plant growth and resistance to abiotic stresses. Drought stress is a vital factor in reducing plant biomass yield and production quality. Alfalfa (Medicago sativa L.) is the most widely planted leguminous forage and one of the most economically valuable crops in the world. Aux/IAA is one of the early responsive gene families of auxin, playing a crucial role in response to drought stress. However, the characteristics of the Aux/IAA gene family in alfalfa and its potential function in response to drought stress are still unknown.

RESULT

A total of 41 Aux/IAA gene members were identified in alfalfa genome. The physicochemical, peptide structure, secondary and tertiary structure analysis of proteins encoded by these genes revealed functional diversity of the MsIAA gene. A phylogenetic analysis classified the MsIAA genes into I-X classes in two subgroups. And according to the gene domain structure, these genes were classified into typical MsIAA and atypical MsIAA. Gene structure analysis showed that the MsIAA genes contained 1-4 related motifs, and except for the third chromosome without MsIAAs, they were all located on 7 chromosomes. The gene duplication analysis revealed that segmental duplication and tandem duplication greatly affected the amplification of the MsIAA genes. Analysis of the Ka/Ks ratio of duplicated MsAux/IAA genes suggested purification selection pressure was high and functional differences were limited. In addition, identification and classification of promoter cis-elements elucidated that MsIAA genes contained numerous elements associated to phytohormone response and abiotic stress response. The prediction protein-protein interaction network showed that there was a complex interaction between the MsAux/IAA genes. Gene expression profiles were tissue-specific, and MsAux/IAA had a broad response to both common abiotic stress (ABA, salt, drought and cold) and heavy metal stress (Al and Pb). Furthermore, the expression patterns analysis of 41 Aux/IAA genes by the quantitative reverse transcription polymerase chain reaction (qRT-PCR) showed that Aux/IAA genes can act as positive or negative factors to regulate the drought resistance in alfalfa.

CONCLUSION

This study provides useful information for the alfalfa auxin signaling gene families and candidate evidence for further investigation on the role of Aux/IAA under drought stress. Future studies could further elucidate the functional mechanism of the MsIAA genes response to drought stress.

Widely Targeted Metabolomic Profiling Combined with Transcriptome Analysis Provides New Insights into Lipid Biosynthesis in Seed Kernels of Pinus koraiensis

Lipid-rich Pinus koraiensis seed kernels are highly regarded for their nutritional and health benefits. To ascertain the molecular mechanism of lipid synthesis, we conducted widely targeted metabolomic profiling together with a transcriptome analysis of the kernels in P. koraiensis cones at various developmental stages. The findings reveal that 148 different types of lipid metabolites, or 29.6% of total metabolites, are present in kernels. Among those metabolites, the concentrations of linoleic acid, palmitic acid, and α-linolenic acid were higher, and they steadily rose as the kernels developed. An additional 10 hub genes implicated in kernel lipid synthesis were discovered using weighted gene co-expression network analysis (WGCNA), gene interaction network analysis, oil body biosynthesis, and transcriptome analysis. This study used lipid metabolome and transcriptome analyses to investigate the mechanisms of key regulatory genes and lipid synthesis molecules during kernel development, which served as a solid foundation for future research on lipid metabolism and the creation of P. koraiensis kernel food.

Cloning and expression analysis of VrNAC13 gene in mung bean

To explore the role of NAC transcription factors in mung bean (Vigna ratiata), we here comprehensively analyzed VrNAC13 structure and expression patterns in the mung bean cultivar "Yulin No.1". The nucleotide sequence of VrNAC13 (GenBank accession number xp014518431.1) was determined by cloning and sequencing the gene. A predicted transcriptional activation domain in VrNAC13 was validated with a yeast one-hybrid assay. The composition and functional characteristics of VrNAC13 were analyzed using basic bioinformatics techniques, and the expression characteristics of VrNAC13 were analyzed via quantitative reverse transcription-PCR. The results showed that VrNAC13 was 1,068 bp in length and encoded a product of 355 amino acids. VrNAC13 was predicted to contain a NAM domain and to belong to the NAC transcription factor family. The protein was hydrophilic and contained several threonine phosphorylation sites. Phylogenetic analysis showed that VrNAC13 was highly similar in sequence to two Arabidopsis thaliana NAC proteins; we hypothesize that VrNAC13 may perform functions in mung bean similar to those of the two closely related proteins in Arabidopsis. Promoter analysis of VrNAC13 revealed cis-acting elements predicted to respond to abscisic acid (ABA), gibberellin, auxin, light, drought, low temperature, and other stressors. VrNAC13 was most highly expressed in the leaves and expressed at very low levels in the stem and root. It was experimentally determined to be induced by drought and ABA. Based on these results, VrNAC13 appears to regulate stress resistance in mung bean.

Genomic characterization of the NAC transcription factors, directed at understanding their functions involved in endocarp lignification of iron walnut (Juglans sigillata Dode)

The NAC (NAM, ATAF1/2, and CUC2) transcription factors (TF), one of the largest plant-specific gene families, play important roles in the regulation of plant growth and development, stress response and disease resistance. In particular, several NAC TFs have been identified as master regulators of secondary cell wall (SCW) biosynthesis. Iron walnut (Juglans sigillata Dode), an economically important nut and oilseed tree, has been widely planted in the southwest China. The thick and high lignified shell derived endocarp tissues, however, brings troubles in processing processes of products in industry. It is indispensable to dissect the molecular mechanism of thick endocarp formation for further genetic improvement of iron walnut. In the present study, based on genome reference of iron walnut, 117 NAC genes, in total, were identified and characterized in silico, which involves only computational analysis to provide insight into gene function and regulation. We found that the amino acids encoded by these NAC genes varied from 103 to 1,264 in length, and conserved motif numbers ranged from 2 to 10. The JsiNAC genes were unevenly distributed across the genome of 16 chromosomes, and 96 of these genes were identified as segmental duplication genes. Furthermore, 117 JsiNAC genes were divided into 14 subfamilies (A-N) according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana and common walnut (Juglans regia). Furthermore, tissue-specific expression pattern analysis demonstrated that a majority of NAC genes were constitutively expressed in five different tissues (bud, root, fruit, endocarp, and stem xylem), while a total of 19 genes were specifically expressed in endocarp, and most of them also showed high and specific expression levels in the middle and late stages during iron walnut endocarp development. Our result provided a new insight into the gene structure and function of JsiNACs in iron walnut, and identified key candidate JsiNAC genes involved in endocarp development, probably providing mechanistic insight into shell thickness formation across nut species.

Characterization of the NAC Transcription Factor in Passion Fruit (Passiflora edulis) and Functional Identification of PeNAC-19 in Cold Stress

The NAC (NAM, ATAF and CUC) gene family plays an important role in plant development and abiotic stress response. However, up to now, the identification and research of the NAC (PeNAC) family members of passion fruit are still lacking. In this study, 25 PeNACs were identified from the passion fruit genome, and their functions under abiotic stress and at different fruit-ripening stages were analyzed. Furthermore, we analyzed the transcriptome sequencing results of PeNACs under four various abiotic stresses (drought, salt, cold and high temperature) and three different fruit-ripening stages, and verified the expression results of some genes by qRT-PCR. Additionally, tissue-specific analysis showed that most PeNACs were mainly expressed in flowers. In particular, PeNAC-19 was induced by four various abiotic stresses. At present, low temperatures have seriously endangered the development of passion fruit cultivation. Therefore, PeNAC-19 was transformed into tobacco, yeast and Arabidopsis to study their function of resisting low temperature. The results show that PeNAC-19 responded to cold stress significantly in tobacco and Arabidopsis, and could improve the low temperature tolerance of yeast. This study not only improved the understanding of the PeNAC gene family characteristics and evolution, but also provided new insights into the regulation of the PeNAC gene at different stages of fruit maturation and abiotic stresses.

Identification and Analysis of the CBF Gene Family in Three Species of Acer under Cold Stress

The C-Repeat Binding Factor (CBF) gene family has been identified and characterized in multiple plant species, and it plays a crucial role in responding to low temperatures. Presently, only a few studies on tree species demonstrate the mechanisms and potential functions of CBFs associated with cold resistance, while our study is a novel report on the multi-aspect differences of CBFs among three tree species, compared to previous studies. In this study, genome-wide identification and analysis of the CBF gene family in Acer truncatum, Acer pseudosieboldianum, and Acer yangbiense were performed. The results revealed that 16 CBF genes (five ApseCBFs, four AcyanCBFs, and seven AtruCBFs) were unevenly distributed across the chromosomes, and most CBF genes were mapped on chromosome 2 (Chr2) and chromosome 11 (Chr11). The analysis of phylogenetic relationships, gene structure, and conserved motif showed that 16 CBF genes could be clustered into three subgroups; they all contained Motif 1 and Motif 5, and most of them only spanned one exon. The cis-acting elements analysis showed that some CBF genes might be involved in hormone and abiotic stress responsiveness. In addition, CBF genes exhibited tissue expression specificity. High expressions of ApseCBF1, ApseCBF3, AtruCBF1, AtruCBF4, AtruCBF6, AtruCBF7, and ApseCBF3, ApseCBF4, ApseCBF5 were detected on exposure to low temperature for 3 h and 24 h. Low expressions of AtruCBF2, AtruCBF6, AtruCBF7 were detected under cold stress for 24 h, and AtruCBF3 and AtruCBF5 were always down-regulated under cold conditions. Taken together, comprehensive analysis will enhance our understanding of the potential functions of the CBF genes on cold resistance, thereby providing a reference for the introduction of Acer species in our country.

Genome-wide identification and expression analysis of the NAC transcription factor family in Saccharum spontaneum under different stresses

The NAC (NAM, ATAF1/2, and CUC2) transcription factor family is one of the largest families unique to plants and is involved in plant growth and development, organs, morphogenesis, and stress responses. The NAC family has been identified in many plants. As the main source of resistance genes for sugarcane breeding, the NAC gene family in the wild species Saccharum spontaneum has not been systematically studied. In this study, 115 SsNAC genes were identified in the S. spontaneum genome, and these genes were heterogeneously distributed on 25 chromosomes. Phylogenetic analysis divided the SsNAC family members into 18 subgroups, and the gene structure and conserved motif analysis further supported the phylogenetic classification. Four groups of tandemly duplicated genes and nine pairs of segmentally duplicated genes were detected. The SsNAC gene has different expression patterns at different developmental stages of stems and leaves. Further qRT-PCR analysis showed that drought, low-temperature, salinity, pathogenic fungi, and other stresses as well as abscisic acid (ABA) and methyl jasmonate (MeJA) treatments significantly induced the expression of 12 SsNAC genes, indicating that these genes may play a key role in the resistance of S. spontaneum to biotic and abiotic stresses. In summary, the results from this study provide comprehensive information on the NAC transcription factor family, providing a reference for further functional studies of the SsNAC gene.

Genome-Wide Identification and Expression Analysis of the Aux/IAA Gene Family of the Drumstick Tree (Moringa oleifera Lam.) Reveals Regulatory Effects on Shoot Regeneration

Auxin plays a critical role in organogenesis in plants. The classical auxin signaling pathway holds that auxin initiates downstream signal transduction by degrading Aux/IAA transcription repressors that interact with ARF transcription factors. In this study, 23 MoIAA genes were identified in the drumstick tree genome. All MoIAA genes were located within five subfamilies based on phylogenetic evolution analysis; the gene characteristics and promoter cis-elements were also analyzed. The protein interaction network between the MoIAAs with MoARFs was complex. The MoIAA gene family responded positively to NAA treatment, exhibiting different patterns and degrees, notably for MoIAA1, MoIAA7 and MoIAA13. The three genes expressed and functioned in the nucleus; only the intact encoding protein of MoIAA13 exhibited transcriptional activation activity. The shoot regeneration capacity in the 35S::MoIAA13-OE transgenic line was considerably lower than in the wild type. These results establish a foundation for further research on MoIAA gene function and provide useful information for improved tissue culture efficiency and molecular breeding of M. oleifera.

Genome-Wide Identification of AP2/ERF Superfamily Genes in Juglans mandshurica and Expression Analysis under Cold Stress

Juglans mandshurica has strong freezing resistance, surviving temperatures as low as -40 °C, making it an important freeze tolerant germplasm resource of the genus Juglans. APETALA2/ethylene responsive factor (AP2/ERF) is a plant-specific superfamily of transcription factors that regulates plant development, growth, and the response to biotic and abiotic stress. In this study, phylogenetic analysis was used to identify 184 AP2/ERF genes in the J. mandshurica genome, which were classified into five subfamilies (JmAP2, JmRAV, JmSoloist, JmDREB, and JmERF). A significant amount of discordance was observed in the 184 AP2/ERF genes distribution of J. mandshurica throughout its 16 chromosomes. Duplication was found in 14 tandem and 122 segmental gene pairs, which indicated that duplications may be the main reason for JmAP2/ERF family expansion. Gene structural analysis revealed that 64 JmAP2/ERF genes contained introns. Gene evolution analysis among Juglandaceae revealed that J. mandshurica is separated by 14.23 and 15 Mya from Juglans regia and Carya cathayensis, respectively. Based on promoter analysis in J. mandshurica, many cis-acting elements were discovered that are related to light, hormones, tissues, and stress response processes. Proteins that may contribute to cold resistance were selected for further analysis and were used to construct a cold regulatory network based on GO annotation and JmAP2/ERF protein interaction network analysis. Expression profiling using qRT-PCR showed that 14 JmAP2/ERF genes were involved in cold resistance, and that seven and five genes were significantly upregulated under cold stress in female flower buds and phloem tissues, respectively. This study provides new light on the role of the JmAP2/ERF gene in cold stress response, paving the way for further functional validation of JmAP2/ERF TFs and their application in the genetic improvement of Juglans and other tree species.

Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant

The NAC (NAM, ATAF1/2, and CUC2) gene family, one of the largest transcription factor families in plants, acts as positive or negative regulators in plant response and adaption to various environmental stresses, including cold stress. Multiple reports on the functional characterization of NAC genes in Arabidopsis thaliana and other plants are available. However, the function of the NAC genes in the typical woody mangrove (Kandelia obovata) remains poorly understood. Here, a comprehensive analysis of NAC genes in K. obovata was performed with a pluri-disciplinary approach including bioinformatic and molecular analyses. We retrieved a contracted NAC family with 68 genes from the K. obovata genome, which were unevenly distributed in the chromosomes and classified into ten classes. These KoNAC genes were differentially and preferentially expressed in different organs, among which, twelve up-regulated and one down-regulated KoNAC genes were identified. Several stress-related cis-regulatory elements, such as LTR (low-temperature response), STRE (stress response element), ABRE (abscisic acid response element), and WUN (wound-responsive element), were identified in the promoter regions of these 13 KoNAC genes. The expression patterns of five selected KoNAC genes (KoNAC6, KoNAC15, KoNAC20, KoNAC38, and KoNAC51) were confirmed by qRT-PCR under cold treatment. These results strongly implied the putative important roles of KoNAC genes in response to chilling and other stresses. Collectively, our findings provide valuable information for further investigations on the function of KoNAC genes.

Comprehensive Analysis of NAC Genes Reveals Differential Expression Patterns in Response to Pst DC3000 and Their Overlapping Expression Pattern during PTI and ETI in Tomato

NAC (NAM/ATAF/CUC) transcription factors belong to a unique gene family in plants, which play vital roles in regulating diverse biological processes, including growth, development, senescence, and in response to biotic and abiotic stresses. Tomato (Solanum lycopersicum), as the most highly valued vegetable and fruit crop worldwide, is constantly attacked by Pseudomonas syringae pv. tomato DC3000 (Pst DC3000), causing huge losses in production. Thus, it is essential to conduct a comprehensive identification of the SlNAC genes involved in response to Pst DC3000 in tomato. In this study, a complete overview of this gene family in tomato is presented, including genome localization, protein domain architectures, physical and chemical features, and nuclear location score. Phylogenetic analysis identified 20 SlNAC genes as putative stress-responsive genes, named SSlNAC 1–20. Expression profiles analysis revealed that 18 of these 20 SSlNAC genes were significantly induced in defense response to Pst DC3000 stress. Furthermore, the RNA-seq data were mined and analyzed, and the results revealed the expression pattern of the 20 SSlNAC genes in response to Pst DC3000 during the PTI and ETI. Among them, SSlNAC3, SSlNAC4, SSlNAC7, SSlNAC8, SSlNAC12, SSlNAC17, and SSlNAC19 were up-regulated against Pst DC3000 during PTI and ETI, which suggested that these genes may participate in both the PTI and ETI pathway during the interaction between tomato and Pst DC3000. In addition, SSlNAC genes induced by exogenous hormones, including indole-3-acetic acid (IAA), abscisic acid (ABA), salicylic acid (SA), and methyl jasmonic acid (MeJA), were also recovered. These results implied that SSlNAC genes may participate in the Pst DC3000 stress response by multiple regulatory pathways of the phytohormones. In all, this study provides important clues for further functional analysis and of the regulatory mechanism of SSlNAC genes under Pst DC3000 stress.

Genome-wide identification and expression analysis of the MADS-box gene family during female and male flower development in Juglans mandshurica

The MADS-box gene family plays a crucial role in multiple developmental processes of plants, especially in floral organ specification and the regulation of fruit development and ripening. Juglans mandshurica is a precious fruit material whose quality and yield are determined by floral organ development. The molecular mechanism of J. mandshurica female and male flower development depending on MADS-box genes remains unclear. In our study, 67 JmMADS genes were identified and unevenly distributed on 15 of 16 J. mandshurica chromosomes. These genes were divided into two types [type I (Mα, Mγ, Mδ) and type II (MIKC)]. The gene structure and motif analyses showed that most genes belonging to the same type had similar gene structures and conserved motifs. The analysis of syntenic relationships showed that MADS-box genes in J. mandshurica, J. sigillata, and J. regia exhibited the highest homology and great collinearity. Analysis of cis-acting elements showed that JmMADS gene promoter regions contained light, stress and hormone response cis-acting elements. The gene expression patterns demonstrated that 30 and 26 JmMADS genes were specifically expressed in the female and male flowers, respectively. In addition, 12 selected genes common to J. mandshurica female and male flowers were significantly upregulated at the mature stage and were used to validate the reliability of the transcriptome data using quantitative real-time PCR. This comprehensive and systematic analysis of J. mandshurica MADS-box genes lays a foundation for future studies on MADS-box gene family functions.

Genome-wide identification and comprehensive analyses of NAC transcription factor gene family and expression analysis under Fusarium kyushuense and drought stress conditions in Passiflora edulis

The NAC gene family is one of the largest plant transcription factors (TFs) families and plays important roles in plant growth, development, metabolism, and biotic and abiotic stresses. However, NAC gene family has not been reported in passion fruit (Passiflora edulis). In this study, a total of 105 NAC genes were identified in the passion fruit genome and were unevenly distributed across all nine-passion fruit chromomere, with a maximum of 48 PeNAC genes on chromosome one. The physicochemical features of all 105 PeNAC genes varied including 120 to 3,052 amino acids, 3 to 8 conserved motifs, and 1 to 3 introns. The PeNAC genes were named (PeNAC001-PeNAC105) according to their chromosomal locations and phylogenetically grouped into 15 clades (NAC-a to NAC-o). Most PeNAC proteins were predicted to be localized in the nucleus. The cis-element analysis indicated the possible roles of PeNAC genes in plant growth, development, light, hormones, and stress responsiveness. Moreover, the PeNAC gene duplications including tandem (11 gene pairs) and segmental (12 gene pairs) were identified and subjected to purifying selection. All PeNAC proteins exhibited similar 3D structures, and a protein-protein interaction network analysis with known Arabidopsis proteins was predicted. Furthermore, 17 putative ped-miRNAs were identified to target 25 PeNAC genes. Potential TFs including ERF, BBR-BPC, Dof, and bZIP were identified in promoter region of all 105 PeNAC genes and visualized in a TF regulatory network. GO and KEGG annotation analysis exposed that PeNAC genes were related to different biological, molecular, and cellular terms. The qRT-PCR expression analysis discovered that most of the PeNAC genes including PeNAC001, PeNAC003, PeNAC008, PeNAC028, PeNAC033, PeNAC058, PeNAC063, and PeNAC077 were significantly upregulated under Fusarium kyushuense and drought stress conditions compared to controls. In conclusion, these findings lay the foundation for further functional studies of PeNAC genes to facilitate the genetic improvement of plants to stress resistance.

Characterization of Phytohormones and Transcriptomic Profiling of the Female and Male Inflorescence Development in Manchurian Walnut (Juglans mandshurica Maxim.)

Flowers are imperative reproductive organs and play a key role in the propagation of offspring, along with the generation of several metabolic products in flowering plants. In Juglans mandshurica, the number and development of flowers directly affect the fruit yield and subsequently its commercial value. However, owing to the lack of genetic information, there are few studies on the reproductive biology of Juglans mandshurica, and the molecular regulatory mechanisms underlying the development of female and male inflorescence remain unclear. In this study, phytohormones and transcriptomic sequencing analyses at the three stages of female and male inflorescence growth were performed to understand the regulatory functions underlying flower development. Gibberellin is the most dominant phytohormone that regulates flower development. In total, 14,579 and 7188 differentially expressed genes were identified after analyzing the development of male and female flowers, respectively, wherein, 3241 were commonly expressed. Enrichment analysis for significantly enriched pathways suggested the roles of MAPK signaling, phytohormone signal transduction, and sugar metabolism. Genes involved in floral organ transition and flowering were obtained and analyzed; these mainly belonged to the M-type MADS-box gene family. Three flowering-related genes (SOC1/AGL20, ANT, and SVP) strongly interacted with transcription factors in the co-expression network. Two key CO genes (CO3 and CO1) were identified in the photoperiod pathway. We also identified two GA20xs genes, one SVP gene, and five AGL genes (AGL8, AGL9, AGL15, AGL19, and AGL42) that contributed to flower development. The findings are expected to provide a genetic basis for the studies on the regulatory networks and reproductive biology in inflorescence development for J. mandshurica.

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.

Genome-Wide Identification and Analysis of the WRKY Gene Family and Cold Stress Response in Acer truncatum

WRKY transcription factors constitute one of the largest gene families in plants and are involved in many biological processes, including growth and development, physiological metabolism, and the stress response. In earlier studies, the WRKY gene family of proteins has been extensively studied and analyzed in many plant species. However, information on WRKY transcription factors in Acer truncatum has not been reported. In this study, we conducted genome-wide identification and analysis of the WRKY gene family in A. truncatum, 54 WRKY genes were unevenly located on all 13 chromosomes of A. truncatum, the highest number was found in chromosomes 5. Phylogenetic relationships, gene structure, and conserved motif identification were constructed, and the results affirmed 54 AtruWRKY genes were divided into nine subgroup groups. Tissue species analysis of AtruWRKY genes revealed which were differently exhibited upregulation in flower, leaf, root, seed and stem, and the upregulation number were 23, 14, 34, 18, and 8, respectively. In addition, the WRKY genes expression in leaf under cold stress showed that more genes were significantly expressed under 0, 6 and 12 h cold stress. The results of this study provide a new insight the regulatory function of WRKY genes under abiotic and biotic stresses.