Secondary Cell Wall Biosynthesis Research Articles

Induction of secondary cell wall biosynthesis genes and their regulators by melatonin in Glycine max

Induction of secondary cell wall biosynthesis genes and their regulators by melatonin in Glycine max

Precocious cell differentiation occurs in proliferating cells in leaf primordia in Arabidopsis angustifolia3 mutant.

During leaf development, the timing of transition from cell proliferation to expansion is an important factor in determining the final organ size. However, the regulatory system involved in this transition remains less understood. To get an insight into this system, we investigated the compensation phenomenon, in which the cell number decreases while the cell size increases in organs with determinate growth. Compensation is observed in several plant species suggesting coordination between cell proliferation and expansion. In this study, we examined an Arabidopsis mutant of ANGUSTIFOLIA 3 (AN3)/GRF-INTERACTING FACTOR 1, a positive regulator of cell proliferation, which exhibits the compensation. Though the AN3 role has been extensively investigated, the mechanism underlying excess cell expansion in the an3 mutant remains unknown. Focusing on the early stage of leaf development, we performed kinematic, cytological, biochemical, and transcriptome analyses, and found that the cell size had already increased during the proliferation phase, with active cell proliferation in the an3 mutant. Moreover, at this stage, chloroplasts, vacuoles, and xylem cells developed earlier than in the wild-type cells. Transcriptome data showed that photosynthetic activity and secondary cell wall biosynthesis were activated in an3 proliferating cells. These results indicated that precocious cell differentiation occurs in an3 cells. Therefore, we suggest a novel AN3 role in the suppression of cell expansion/differentiation during the cell proliferation phase.

In silico characterization of a MYB TF (WsMYB46) from the medicinal plant Withania somnifera L. predicting its probable role in secondary cell wall biosynthesis

In silico characterization of a MYB TF (WsMYB46) from the medicinal plant Withania somnifera L. predicting its probable role in secondary cell wall biosynthesis

Reducing nitrogen application at high planting density enhances secondary cell wall formation and decreases stem lodging in rapeseed

Reducing the rate of nitrogen fertilizer (N rate) and achieving high plant density is an effective and economically beneficial approach to enhance rapeseed productivity. However, it is unclear how N rate affects stem lodging in rapeseed directly-sown at high densities. In this study, the effect of different N rates on stem lodging at high plant density was examined in field trial and pot experiments using two modern hybrid rapeseed varieties, Fengyou520 (FY520) and Huayouza62 (HYZ62). Increasing N fertilizer from low level to adequate level significantly increased yield by 39% in FY520 and 43% in HYZ62. However, with elevated N yield increases were limited by an increase in stem lodging. At T1 (the first inflorescence branch below the apex has the first flowers open) phloroglucinol-HCl staining for lignin was weaker under high N (N4: 4 g N pot−1) compared with adequate N (N2: 2 g N pot−1), and was correlated with weaker stem strength. Differences in stem lignin were also associated with variations in gene expression as revealed by transcriptomic profiling of middle stems at T1. Differentially expressed genes indicated significant enrichment for cell wall. Several genes encoding known transcriptional regulators of secondary cell walls (SCWs) biosynthesis and deposition were down-regulated with high N compared with adequate N. Many genes putatively involved in lignin biosynthesis, cellulose biosynthesis and assembly and xylan biosynthesis are down-regulated with high N, corresponding to lower accumulation of lignin, cellulose and hemicellulose in these plants. These findings not only provide us with gene candidates potentially affecting stem lodging in direct sown rapeseed, but also reveal that reducing N at high plant density enhances biosynthesis of stem cell walls to alleviate stem lodging.

Characterization of ZmCesAs for Secondary Cell Wall Biosynthesis in Maize

Characterization of ZmCesAs for Secondary Cell Wall Biosynthesis in Maize

GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, play a regulatory role in cotton fiber secondary cell wall biosynthesis.

GhHB14_D10 and GhREV_D5 regulated secondary cell wall formation and played an important role in fiber development. Cotton serves as an important source of natural fiber, and the biosynthesis of the secondary cell wall plays a pivotal role in determining cotton fiber quality. Nevertheless, the intricacies of this mechanism in cotton fiber remain insufficiently elucidated. This study investigates the functional roles of GhHB14_D10 and GhREV_D5, two HD-ZIP III transcription factors, in secondary cell wall biosynthesis in cotton fibers. Both GhHB14_D10 and GhREV_D5 were found to be localized in the nucleus with transcriptional activation activity. Ectopic overexpression of GhHB14_D10 and GhREV_D5 in Arabidopsis resulted in changed xylem differentiation, secondary cell wall deposition, and expression of genes related to the secondary cell wall. Silencing of GhHB14_D10 and GhREV_D5 in cotton led to enhanced fiber length, reduced cell wall thickness, cellulose contents and expression of secondary cell wall-related genes. Moreover, GhHB14_D10's direct interaction with GhREV_D5, and transcriptional regulation of cellulose biosynthesis genes GhCesA4-4 and GhCesA7-2 revealed their collaborative roles in secondary cell wall during cotton fiber development. Overall, these results shed light on the roles of GhHB14_D10 and GhREV_D5 in secondary cell wall biosynthesis, offering a strategy for the genetic improvement of cotton fiber quality.

Transcription factor PagMYB31 positively regulates cambium activity and negatively regulates xylem development in poplar.

Wood formation involves consecutive developmental steps, including cell division of vascular cambium, xylem cell expansion, secondary cell wall (SCW) deposition, and programmed cell death. In this study, we identified PagMYB31 as a coordinator regulating these processes in Populus alba × Populus glandulosa and built a PagMYB31-mediated transcriptional regulatory network. PagMYB31 mutation caused fewer layers of cambial cells, larger fusiform initials, ray initials, vessels, fiber and ray cells, and enhanced xylem cell SCW thickening, showing that PagMYB31 positively regulates cambial cell proliferation and negatively regulates xylem cell expansion and SCW biosynthesis. PagMYB31 repressed xylem cell expansion and SCW thickening through directly inhibiting wall-modifying enzyme genes and the transcription factor genes that activate the whole SCW biosynthetic program, respectively. In cambium, PagMYB31 could promote cambial activity through TRACHEARY ELEMENT DIFFERENTIATION INHIBITORY FACTOR (TDIF)/PHLOEM INTERCALATED WITH XYLEM (PXY) signaling by directly regulating CLAVATA3/ESR-RELATED (CLE) genes, and it could also directly activate WUSCHEL HOMEOBOX RELATED4 (PagWOX4), forming a feedforward regulation. We also observed that PagMYB31 could either promote cell proliferation through the MYB31-MYB72-WOX4 module or inhibit cambial activity through the MYB31-MYB72-VASCULAR CAMBIUM-RELATED MADS2 (VCM2)/PIN-FORMED5 (PIN5) modules, suggesting its role in maintaining the homeostasis of vascular cambium. PagMYB31 could be a potential target to manipulate different developmental stages of wood formation.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar.

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood. In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031. Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex. Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

Deciphering the intricate hierarchical gene regulatory network: unraveling multi-level regulation and modifications driving secondary cell wall formation.

Wood quality is predominantly determined by the amount and the composition of secondary cell walls (SCWs). Consequently, unraveling the molecular regulatory mechanisms governing SCW formation is of paramount importance for genetic engineering aimed at enhancing wood properties. Although SCW formation is known to be governed by a hierarchical gene regulatory network (HGRN), our understanding of how a HGRN operates and regulates the formation of heterogeneous SCWs for plant development and adaption to ever-changing environment remains limited. In this review, we examined the HGRNs governing SCW formation and highlighted the significant key differences between herbaceous Arabidopsis and woody plant poplar. We clarified many confusions in existing literatures regarding the HGRNs and their orthologous gene names and functions. Additionally, we revealed many network motifs including feed-forward loops, feed-back loops, and negative and positive autoregulation in the HGRNs. We also conducted a thorough review of post-transcriptional and post-translational aspects, protein-protein interactions, and epigenetic modifications of the HGRNs. Furthermore, we summarized how the HGRNs respond to environmental factors and cues, influencing SCW biosynthesis through regulatory cascades, including many regulatory chains, wiring regulations, and network motifs. Finally, we highlighted the future research directions for gaining a further understanding of molecular regulatory mechanisms underlying SCW formation.

Jasmonate activates secondary cell wall biosynthesis through MYC2-MYB46 module.

Formation of secondary cell wall (SCW) is tightly regulated spatiotemporally by various developmental and environmental signals. Successful fine-tuning of the trade-off between SCW biosynthesis and stress responses requires a better understanding of how plant growth is regulated under environmental stress conditions. However, the current understanding of the interplay between environmental signaling and SCW formation is limited. The lipid-derived plant hormone jasmonate (JA) and its derivatives are important signaling components involved in various physiological processes including plant growth, development, and abiotic/biotic stress responses. Recent studies suggest that JA is involved in SCW formation but the signaling pathway has not been studied for how JA regulates SCW formation. We tested this hypothesis using the transcription factor MYB46, a master switch for SCW biosynthesis, and JA treatments. Both the transcript and protein levels of MYB46, a master switch for SCW formation, were significantly increased by JA treatment, resulting in the upregulation of SCW biosynthesis. We then show that this JA-induced upregulation of MYB46 is mediated by MYC2, a central regulator of JA signaling, which binds to the promoter of MYB46. We conclude that this MYC2-MYB46 module is a key component of the plant response to JA in SCW formation.

Functional Diversification and the Plant Secondary Cell Wall.

Much evidence exists suggesting the presence of genetic functional diversification in plants, though literature associated with the role of functional diversification in the evolution of the plant secondary cell wall (SCW) has sparsely been compiled and reviewed in a recent context. This review aims to elucidate, through the examination of gene phylogenies associated with its biosynthesis and maintenance, the role of functional diversification in shaping the critical, dynamic, and characteristic organelle, the secondary cell wall. It will be asserted that gene families resulting from gene duplication and subsequent functional divergence are present and are heavily involved in SCW biosynthesis and maintenance. Furthermore, diversification will be presented as a significant driver behind the evolution of the many functional characteristics of the SCW. The structure and function of the plant cell wall and its constituents will first be explored, followed by a discussion on the phenomenon of gene duplication and the resulting genetic functional divergence that can emerge. Finally, the major constituents of the SCW and their individual relationships with duplication and divergence will be reviewed to the extent of current knowledge on the subject.

Cotton GhNAC4 promotes drought tolerance by regulating secondary cell wall biosynthesis and ribosomal protein homeostasis.

Drought has a severe impact on the quality and yield of cotton. Deciphering the key genes related to drought tolerance is important for understanding the regulation mechanism of drought stress and breeding drought-tolerant cotton cultivars. Several studies have demonstrated that NAC transcription factors are crucial in the regulation of drought stress, however, the related functional mechanisms are still largely unexplored. Here, we identified that NAC transcription factor GhNAC4 positively regulated drought stress tolerance in cotton. The expression of GhNAC4 was significantly induced by abiotic stress and plant hormones. Silencing of GhNAC4 distinctly impaired the resistance to drought stress and overexpressing GhNAC4 in cotton significantly enhanced the stress tolerance. RNA-seq analysis revealed that overexpression of GhNAC4 enriched the expression of genes associated with the biosynthesis of secondary cell walls and ribosomal proteins. We confirmed that GhNAC4 positively activated the expressions of GhNST1, a master regulator reported previously in secondary cell wall formation, and two ribosomal protein-encoding genes GhRPL12 and GhRPL18p, by directly binding to their promoter regions. Overexpression of GhNAC4 promoted the expression of downstream genes associated with the secondary wall biosynthesis, resulting in enhancing secondary wall deposition in the roots, and silencing of GhRPL12 and GhRPL18p significantly impaired the resistance to drought stress. Taken together, our study reveals a novel pathway mediated by GhNAC4 that promotes secondary cell wall biosynthesis to strengthen secondary wall development and regulates the expression of ribosomal protein-encoding genes to maintain translation stability, which ultimately enhances drought tolerance in cotton.

Effect of the HCT gene on lignin synthesis and fiber development in Gossypium barbadense

As one of the key enzymes in the metabolic pathway of phenylpropane, shikimate hydroxycinnamoyl transferase (HCT) is mainly involved in the biosynthesis of the plant secondary cell wall, which is closely related to cotton fiber quality. In this study, whole-genome identification and bioinformatics analysis of the HCT gene family were performed in G. barbadense. In the whole genome, we identified 136 GbHCT genes encoding 309–504 amino acids. Phylogenetic analysis divided the genome into 5 subfamilies, which were located on 25 chromosomes. Collinear analysis of polyploidization and tandem duplication events were the main driving forces for the rapid expansion and evolution of this family, and the genes underwent loose purifying selection constraints after duplication. Gene promoters identified a variety of cis-acting elements related to plant hormones and the stress response. Several members of the GbHCT family were highly expressed during the development of cotton fiber, and different members had different expression patterns in cotton fiber. After GbHCT114 gene silencing in cotton, the amount of stem surface trichomes and lignin content decreased, and the cell morphology and arrangement changed. After the GbHCT114 gene was overexpressed in Arabidopsis thaliana (L.) Heynh., the number of stem and leaf surface trichomes and the cross-sectional area of the secondary xylem duct cell wall increased. In addition, utilizing transcriptomic analysis, differentially expressed genes associated with lignin synthesis and fiber development were identified. Taken together, the results obtained in this study confirm that the GbHCT114 gene regulates plant trichome development, which lays a theoretical foundation for future studies on the function of GbHCT114 in cotton.

Tissue-Specific Transcriptomes in the Secondary Cell Wall Provide an Understanding of Stem Growth Enhancement in Solidago canadensis during Invasion.

Invasive plants generally present a significant enhancement in aboveground vegetative growth, which is mainly caused by variation in secondary cell wall (SCW) deposition and vascular tissue development. However, the coordination of the transcriptional regulators of SCW biosynthesis is complex, and a comprehensive regulation map has not yet been clarified at a transcriptional level to explain the invasive mechanism of S. canadensis. Here, RNA sequencing was performed in the phloem and xylem of two typical native (US01) and invasive (CN25) S. canadensis populations with different stem morphologies. A total of 296.14 million high-quality clean reads were generated; 438,605 transcripts and 156,968 unigenes were assembled; and 66,648 and 19,510 differential expression genes (DEGs) were identified in the phloem and xylem, respectively. Bioinformatics analysis indicated that the SCW transcriptional network was dramatically altered during the successful invasion of S.canadensis. Based on a comprehensive analysis of SCW deposition gene expression profiles, we revealed that the invasive population is dedicated to synthesizing cellulose and reducing lignification, leading to an SCW with high cellulose content and low lignin content. A hypothesis thus has been proposed to explain the enhanced stem growth of S. canadensis through the modification of the SCW composition.

Stone cell formation in the pedicel of pears and apples.

For the first time, stone cells in pear and apple pedicel were studied. The lignification of the pedicel outer part was correlated with flesh, and the secondary cell wall biosynthesis genes were activated. Fruit pedicels act as bridges between the fruit and the shoot. They have secondary thickened cell walls that presumably function in mechanical support, water and nutrient transport. Stone cells are cells with a secondary cell wall thickening. In pears, yet not in apples, the stone cells affect the flesh texture. There have been few reports on stone cell formation in pear and apple pedicels; therefore, we studied these cells for the first time. The apple pedicel had few stone cells in the cortex. The formation of stone cells in pear continued until seven weeks after flowering (WAF), and the density was significantly higher than in apple. The stone cell formation degree (SFD) of pear was 3.6-7.1 times higher than that of apple. Total lignin and lignin non-condensed structure (G and S units) content in the pear pedicle outer part was 1.5-2.7 times higher than that of the apple at harvest. The SFD of the pedicel outer part had a positive correlation with the G and S units content of the flesh. The total lignin and G and S units content between flesh and the pedicel outer part were positively correlated. Correlation analysis revealed a positive relationship between fruit and pedicel formation of the stone cells. The WGCNA showed that NST3 was linked to NAC028, MYB46, CESA, POD, LAC, and VSR6. These genes were highly expressed in the outer part of the pear pedicel, while they were suppressed in that issue of the apple at 4 WAF.

Association Mapping and Expression Analysis of the Genes Involved in the Wood Formation of Poplar.

Xylogenesis is a complex and sequential biosynthetic process controlled by polygenes. Deciphering the genetic architecture of this complex quantitative trait could provide valuable information for increasing wood biomass and improving its properties. Here, we performed genomic resequencing of 64 24-year-old trees (64 hybrids of section Aigeiros and their parents) grown in the same field and conducted full-sib family-based association analyses of two growth and six woody traits using GEMMA as a choice of association model selection. We identified 1342 significantly associated single nucleotide polymorphisms (SNPs), 673 located in the region upstream and downstream of 565 protein-encoding genes. The transcriptional regulation network of secondary cell wall (SCW) biosynthesis was further constructed based on the published data of poplar miRNA, transcriptome, and degradome. These provided a certain scientific basis for the in-depth understanding of the mechanism of poplar timber formation and the molecular-assisted breeding in the future.

The transcription factor ERF108 interacts with AUXIN RESPONSE FACTORs to mediate cotton fiber secondary cell wall biosynthesis.

Phytohormones play indispensable roles in plant growth and development. However, the molecular mechanisms underlying phytohormone-mediated regulation of fiber secondary cell wall (SCW) formation in cotton (Gossypium hirsutum) remain largely underexplored. Here, we provide mechanistic evidence for functional interplay between the AP2/ERF transcription factor GhERF108 and auxin response factors GhARF7-1 and GhARF7-2 in dictating the ethylene-auxin signaling crosstalk that regulates fiber SCW biosynthesis. Specifically, in vitro cotton ovule culture revealed that ethylene and auxin promote fiber SCW deposition. GhERF108 RNAi cotton displayed remarkably reduced cell wall thickness compared with controls. GhERF108 interacted with GhARF7-1 and GhARF7-2 to enhance the activation of the MYB transcription factor gene GhMYBL1 (MYB domain like protein 1) in fibers. GhARF7-1 and GhARF7-2 respond to auxin signals that promote fiber SCW thickening. GhMYBL1 RNAi and, GhARF7-1 and GhARF7-2 VIGS cotton displayed similar defects in fiber SCW formation as GhERF108 RNAi cotton. Moreover, the ethylene and auxin responses were reduced in GhMYBL1 RNAi plants. GhMYBL1 directly binds to the promoters of GhCesA4-1, GhCesA4-2, and GhCesA8-1 and activates their expression to promote cellulose biosynthesis, thereby boosting fiber SCW formation. Collectively, our findings demonstrate that the collaboration between GhERF108 and GhARF7-1 or GhARF7-2 establishes ethylene-auxin signaling crosstalk to activate GhMYBL1, ultimately leading to the activation of fiber SCW biosynthesis.

Two R2R3-MYB transcription factors from Chinese cedar (Cryptomeria fortunei Hooibrenk) are involved in the regulation of secondary cell wall formation

As the most abundant renewable energy source, wood comprises the secondary cell wall (SCW). SCW biosynthesis involves lignin and cellulose deposition. Increasing studies have illustrated that R2R3-MYB transcription factors (TFs) play pivotal roles in affecting lignin accumulation and SCW formation. Nevertheless, the regulatory roles of R2R3-MYBs are still unresolved in Cryptomeria fortunei Hooibrenk cambium and wood formation. To dissect the potentials of CfMYBs, we successfully cloned and intensively studied the functions of CfMYB4 and CfMYB5 in SCW formation and abiotic stress response. They both contained the conserved MYB domain capable of forming a special structure that could bind to the core motifs of downstream genes. The phylogenetic tree implied that two CfMYBs clustered into different evolutionary branches. They were predominantly expressed in the stem and were localized to the nucleus. Furthermore, CfMYB4 functioned as an activator to enhance lignin and cellulose accumulation, and increase the SCW thickness by elevating the expression levels of SCW-related genes. By contrast, CfMYB5 negatively regulated lignin and cellulose biosynthesis, and decreased SCW formation by reducing the expression of SCW biosynthetic genes. Our data not only highlight the regulatory functions of CfMYBs in lignin deposition but also provide critical insights into the development of strategies for the genetic improvement of Cryptomeria fortunei wood biomass.

Genome-wide identification and expression profiling of B3 transcription factor genes in Populus alba × Populus glandulosa.

B3-domain containing transcription factors (TFs) are well known to play important roles in various developmental processes, including embryogenesis, seed germination, etc. Characterizations and functional studies of the B3 TF superfamily in poplar are still limited, especially on their roles in wood formation. In this study, we conducted comprehensive bioinformatics and expression analysis of B3 TF genes in Populus alba × Populus glandulosa. A total of 160 B3 TF genes were identified in the genome of this hybrid poplar, and their chromosomal locations, syntenic relationships, gene structures, and promoter cis-acting elements were analyzed. Through domain structure and phylogenetic relationship analyses, these proteins were classified into four families LAV, RAV, ARF, and REM. Domain and conservation analyses revealed different gene numbers and different DNA-binding domains among families. Syntenic relationship analysis suggested that approximately 87% of the genes resulted from genome duplication (segmental or tandem), contributing to the expansion of the B3 family in P. alba × P. glandulosa. Phylogeny in seven species revealed the evolutionary relationship of B3 TF genes across different species. B3 domains among the eighteen proteins that were highly expressed in differentiating xylem had a high synteny, suggesting a common ancestor for these seven species. We performed co-expression analysis on the representative genes in two different ages of poplar, followed by pathways analysis. Among those genes co-expressed with four B3 genes, 14 were involved in lignin synthases and secondary cell walls biosynthesis, including PagCOMT2, PagCAD1, PagCCR2, PagCAD1, PagCCoAOMT1, PagSND2, and PagNST1. Our results provide valuable information for the B3 TF family in poplar and show the potential of B3 TF genes in engineering to improve wood properties.

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.