Transcriptomics Using the Enriched Arabidopsis Shoot Apex Reveals Developmental Priming Genes Involved in Plastic Plant Growth under Salt Stress Conditions

In the shoot apical meristem (SAM), the homeostasis of the stem cell population supplying new cells for organ formation is likely a key mechanism of multicellular plant growth and development. As plants are sessile organisms and constantly encounter environmental abiotic stresses, postembryonic development from the shoot stem cell population must be considered with surrounding abiotic stresses for plant adaptation. However, the underlying molecular mechanisms for plant adaptation remain unclear. Previous studies found that the stem-cell-related mutant clv3-2 has the property of salt tolerance without the differential response of typical stress-responsive genes compared to those in WT Ler. Based on these facts, we hypothesized that shoot meristems contain developmental priming genes having comprehensively converged functions involved in abiotic stress response and development. To better understand the biological process of developmental priming genes in the SAM, we performed RNA sequencing (RNA-seq) and transcriptome analysis through comparing genome-wide gene expression profiles between enriched shoot apex and leaf tissues. As a result, 121 putative developmental priming genes differentially expressed in the shoot apex compared to the leaf were identified under normal and salt stress conditions. RNA-seq experiments also revealed the shoot apex-specific responsive genes for salt stress conditions. Based on combinatorial comparisons, 19 developmental priming genes were finally identified, including developmental genes related to cell division and abiotic/biotic-stress-responsive genes. Moreover, some priming genes showed CLV3-dependent responses under salt stress conditions in the clv3-2. These results presumably provide insight into how shoot meristem tissues have relatively high viability against stressful environmental conditions for the developmental plasticity of plants.


Introduction
In plants, the shoot apical meristem (SAM) activities, including cell division and differentiation, are regulated mainly by the CLAVATA3 peptide (CLV3p)-CLV receptor complex mediated signaling pathway, negatively inhibiting the WUSCHEL (WUS) transcription factor (TF) [1]. WUS proteins expressed in the organizing center cells move toward upper stem cells via plasmodesmata and activate CLV3 expression directly in stem cells to promote the proliferation of undifferentiated stem cells [1]. The processed CLV3p are then secreted from stem cells to the neighboring cells and CLV receptor complexes, including CLV1 and BARELY ANY MERISTEM 1 (BAM1), recognize CLV3p ligands to induce downstream MAPK signaling to inhibit WUS expression [2]. These negative feedback loops between CLV3 and WUS eventually regulate the homeostasis of the stem cell population to control proper growth and postembryonic development of the aboveground lateral organs.
In nature, plants respond continuously to various environmental challenges during their lifespan because of a sessile property, and those abiotic/biotic stimuli usually affect plant growth and development. Therefore, plants must adapt to their growth, including endogenous developmental processes and morphological changes via their plastic properties we found 121 differentially expressed genes (DEGs) that displayed various criteria of gene ontology (GO) terms, such as development, abiotic/biotic stress response, phytohormone, and metabolism via GO enrichment analysis. RNA-seq experiments also revealed the shoot apex-specific responsive genes for salt stress conditions. Based on combinatorial comparisons, 19 developmental priming genes were finally identified, including the developmental genes related to cell division and abiotic/biotic-stress-responsive genes. Moreover, some priming genes showed CLV3-dependent responses under salt stress conditions in the clv3-2. These results show how shoot meristem tissues have a relatively high tolerance to stress conditions for the developmental plasticity of plants.

Developmental Plasticity Involved in the Shoot Growth and Development against Salt Stress Conditions
Soil-grown mutant plants harboring deficiencies in shoot meristem development showed the salt-tolerant phenotype [15]. Consistent with this, seedling growth of shoot meristem-related mutants, such as clv1 bam1 and clv3-2 displayed enhanced survival rates under 100 and 150 mM NaCl conditions compared to the corresponding WT Columbia-0 (Col-0) and Landsberg erecta (Ler), respectively ( Figure 1). Because leucine-rich repeatreceptor-like kinases (LRR-RLKs) CLV1 and BAM1 act as receptors for CLV3p in the CLV-WUS negative feedback pathway regulating stem cell homeostasis in the SAM [1,2], these mutants usually display increased phenotypes of shoot meristem size and floral organ numbers, which are not related to stress tolerance. Therefore, this finding suggests that stem cell signaling consisting of peptide signals and receptors may function in abiotic resistance and maintenance in the SAM. In addition, WT Ler showed more sensitive responses to shoot growth than another WT Col-0 under salt stress conditions ( Figure 1). Although Col-0 seedlings showed reduced growth under 100 mM NaCl conditions compared to control conditions ( Figure 1A), the survival rate was little affected ( Figure 1B). In contrast, the survival rate, including defects of germination rate and shoot growth, was reduced significantly under 150 mM NaCl conditions ( Figure 1B). Unlike Col-0, those defects of Ler seedlings appeared under 100 mM NaCl conditions ( Figure 1B). Consistently, it has been known that rosette size and ion leakage in Ler were reduced and increased compared to those in Col-0, respectively, in the analysis using 160 Arabidopsis accessions under salt stress conditions [16]. These results suggest that the genetic information depending on Arabidopsis accessions is also essential for developmental plasticity on shoot growth against unfavorable salt stress conditions.

Identification of Developmental Priming Genes in the Shoot Apex
Because previous studies reported that canonical stress-responsive genes, including KIN1, RD29A, RAB18, and DREB2A, were not distinguishable between WT and clv3-2 [15], we investigated the biological role of developmentally primed (intrinsically accumulated) genes acting as stress regulators in the shoot meristem compared to other leaf tissues by performing RNA sequencing (RNA-seq) to understand how shoot-stem-cell-related mutants show the salt-tolerant phenotype. For this, shoot apices and leaves of more saltsensitive WT Ler seedlings grown for nine days and treated without or with 200 mM NaCl for one day were harvested separately. A recent study showed that RNA-seq results using enriched shoot apices could represent SAM-specific tendency compared to those using functional zone/layer-specific cells [14,17].
Since heat-shock priming genes expressed in the SAM were previously found under moderate and severe stress conditions in RNA-seq analysis [14], we are also hypothesized that developmental priming genes are significantly expressed in the SAM than other tissues under both control (but primed, see the Introduction) and salt stress conditions. Therefore, we first isolated differentially expressed genes (DEGs) in the shoot apices compared to those in leaves under non-stress conditions in all three repetitive RNA-seq experiments (Figures 2A,B and S1). Four hundred and seventy-three DEGs in the shoot apex were identified; 216 genes were upregulated, and 257 genes were downregulated (Figure 2A ,B and  Table S1). Second, 484 DEGs in the shoot apex were isolated; 148 genes were upregulated, and 336 genes were downregulated under salt stress conditions (Figures 2C,D and S2 and  Table S2). Overall, 121 developmental priming genes were expressed significantly higher or lower in the shoot apex under both non-stress control and stress conditions ( Figure 2E,F and Table S3).
Gene Ontology (GO) enrichment analysis was performed using DAVID bioinformatics resources to determine the biological function of developmental priming genes in the SAM against salt stress conditions. Among the 121 DEGs, 70 and 51 genes were upregulated and downregulated, respectively, under control and stress conditions ( Figure 2E). Based on the criteria of enrichment, 14 GO terms in biological process (BP) and 11 GO terms in the molecular function (MF) were identified in upregulated 70 DEGs ( Figure 3A). An analysis of GO terms in biological processes included the following: auxin biosynthetic process ( (Table S4). Consistent with the BP category, the highest value of GO enrichment in the MF category also revealed auxin-related IAA-amido synthetase activity (GO:0010279) with 106.7-fold enrichment, and transcription factor activity (GO:0003700) included various types of transcription factors involved in cell proliferation and auxin regulation (Table 1). Interestingly, several genes related to the regulation of phytohormones, such as auxin and brassinosteroid, were in- Gene Ontology (GO) enrichment analysis was performed using DAVID bioinformatics resources to determine the biological function of developmental priming genes in the SAM against salt stress conditions. Among the 121 DEGs, 70 and 51 genes were upregulated and downregulated, respectively, under control and stress conditions ( Figure 2E). Based on the criteria of enrichment, 14 GO terms in biological process (BP) and 11 GO terms in the molecular function (MF) were identified in upregulated 70 DEGs ( Figure 3A). the BP category, the highest value of GO enrichment in the MF category also revealed auxin-related IAA-amido synthetase activity (GO:0010279) with 106.7-fold enrichment, and transcription factor activity (GO:0003700) included various types of transcription factors involved in cell proliferation and auxin regulation (Table 1). Interestingly, several genes related to the regulation of phytohormones, such as auxin and brassinosteroid, were involved in plant growth and cell elongation/division by a KEGG enrichment analysis of upregulated DEGs in the shoot apex ( Figure S3). In the downregulated DEGs, 11 GO terms in BP and 9 GO terms in MF were identified ( Figure 3B and Table S5). In the BP category, enriched GO terms included the following: regulation of the defense response by callose (GO:2000071) with an enrichment value of 177.5-fold; response to oomycetes (GO:0002239) with 55.5-fold enrichment; SA mediated signaling pathway (GO:0009863) with 32.9-fold enrichment; cellular response to hypoxia (GO:0071456) with 7.3-fold enrichment; response to oxidative stress (GO:0006979) with 5.7-fold enrichment; defense response to fungus (GO:0050832) with 5.6-fold enrichment; defense response (GO:0006952) with 5.4-fold enrichment; signal transduction (GO:0007165) with 4.7-fold enrichment; defense response to bacteria (GO:0042742) with 4.4-fold enrichment; response to abscisic acid (GO:0009737) with 4.1-fold enrichment; and  In the downregulated DEGs, 11 GO terms in BP and 9 GO terms in MF were identified ( Figure 3B and Table S5). In the BP category, enriched GO terms included the following: regulation of the defense response by callose (GO:2000071) with an enrichment value of 177.5-fold; response to oomycetes (GO:0002239) with 55.5-fold enrichment; SA mediated signaling pathway (GO:0009863) with 32.9-fold enrichment; cellular response to hypoxia (GO:0071456) with 7.3-fold enrichment; response to oxidative stress (GO:0006979) with 5.7-fold enrichment; defense response to fungus (GO:0050832) with 5.6-fold enrichment; defense response (GO:0006952) with 5.4-fold enrichment; signal transduction (GO:0007165) with 4.7-fold enrichment; defense response to bacteria (GO:0042742) with 4.4-fold enrichment; response to abscisic acid (GO:0009737) with 4.1-fold enrichment; and protein phosphorylation (GO:0006468) with 2.9-fold enrichment ( Table S5). The GO terms in BP were significantly related to biotic and abiotic stress responses. Analysis through the PANTHER classification system (http://www.pantherdb.org accessed on 22 September 2022) showed that even those genes included in the GO terms, such as signal transduction and protein phosphorylation, are primarily involved in the disease responses. Consistent with the function based on GO terms in the BP category, the GO terms of downregulated priming genes in the MF category mainly displayed the catalytic enzyme activities acting in the abiotic/biotic stress responses.

Clustering Analysis of Developmental Priming Genes in the Shoot Apex
K-Means/Medians Clustering analysis with average fold changes of 107 genes from 121 DEGs, expressed differentially in the shoot apex under both control and stress conditions, was performed to investigate the developmental priming gene more precisely. In this analysis, 14 genes showing no read (zero) value in the RNA-seq results were removed because of the unavailable fold ratio (Table S6). Twelve of them were classified as an upregulated state in developmental priming genes. Seven encoded several types of transcription factors, including KNOX ARABIDOPSIS THALIANA MEINOX (KNATM), SHOOT MERISTEMLESS (STM), and MYB117, involved in regulating shoot and axillary meristem development [31,37,38].
As a result, two clustering groups showed increased or decreased patterns between control and NaCl conditions via the Pearson Correlation distance metric ( Figure 4 and Table S7). In the top 10 genes of increased DEGs in clustering analysis, most genes showed that increased patterns were changed in downregulated priming state. Those genes are involved in plant growth and biotic/abiotic stress responses (Table 2). In the top 10 genes of decreased DEGs in clustering analysis, the most dramatic change in the GH3.3 and ARA-BINOGALACTAN PROTEIN 1 (AGP1) genes acting in the regulation of salt tolerance and reproductive development occurred between the different primed states from upregulated to downregulated states ( Table 2).

Identification of Salt-Inducible Genes Only in the Shoot Apex
The control and NaCl-treated RNA-seq results of shoot apex or leaf tissues were compared to find salt-inducible genes specifically enriched in the shoot apex. DEGs in the shoot apices treated with NaCl compared to those under control conditions in all three repetitive RNA-seq experiments were isolated ( Figures 5A,B and S4). Four hundred and ninety-nine differentially expressed genes were isolated in the shoot apex; 327 genes were upregulated by the NaCl treatment, and 172 genes were downregulated ( Figure 5A,B and Table S8). Four hundred and seventy-six DEGs were identified in the leaf tissue compared to those under control conditions; 377 genes were upregulated under salt stress conditions, and 99 genes were downregulated ( Figures 5C,D and S5 and Table S9). By comparing salt-inducible DEGs between the shoot apex and leaf, 379 preferentially enriched DEGs were found only in the shoot apex, 356 DEGs only in the leaf, and 120 DEGs in both the shoot apex and leaf tissues ( Figure 5E). Interestingly, the GO terms related to development were revealed to the largest portion of preferentially enriched in the shoot apex (40%) but not those in the RNA-seq results of overlap (3%) and leaf (4%). In contrast, GO terms related to abiotic/biotic response showed the largest portion in overlap (58%) and leaf (44%) than the shoot apex (22%) ( Figure 5F). The GO terms related to development in the shoot apex mainly included developmental processes involving cell division. Although all three categories contained GO terms related to leaf senescence (GO:0010150), those genes in the three categories were not overlapped, suggesting that they may function differently in developmental plasticity (Table S10).

Identification of Developmental Priming Genes Related to Salt in the Shoot Apex
Based on the previous comparisons, Venn diagram analysis using developmental priming genes and salt-inducible DEGs revealed 117 developmental priming genes expressed more specifically in the shoot apex, of which 19 genes also showed salt-inducible changes only in the shoot apex ( Figure 6A and Table S11). This suggests that these 19

Identification of Developmental Priming Genes Related to Salt in the Shoot Apex
Based on the previous comparisons, Venn diagram analysis using developmental priming genes and salt-inducible DEGs revealed 117 developmental priming genes expressed more specifically in the shoot apex, of which 19 genes also showed salt-inducible changes only in the shoot apex ( Figure 6A and Table S11). This suggests that these 19 genes presumably play essential roles in developmental plasticity under salt stress conditions. According to GO enrichment analysis using DAVID, 98 genes displayed GO terms related to the biotic response, hormone regulation, organ development, and metabolic response, while 19 genes mainly showed GO terms related to the regulation of cell cycle/division ( Figure 6B,C). Although 10 of the 19 genes were not included in GO enrichment analysis using DAVID, individual analysis of those genes revealed GO terms related to abiotic/biotic stress responses (Table S12). Under control conditions, the genes associated with cell division and microtubule-based movement were included as upregulated priming genes in the shoot apex compared to other leaf tissue. By contrast, the salt treatment reduced the expression change of these genes. On the other hand, the expression of abiotic/bioticstress-responsive genes associated with downregulated priming genes in the shoot apex under normal conditions was induced by salt stress conditions. These results revealed the functional behavior of developmental priming genes depending on salt stresses in the WT SAM.
Because previous studies showed that the clv3-2 mutant exhibiting the increase in cell division or the decrease in cell differentiation has the property of salt resistance compared to WT Ler via the non-canonical responses [15], qRT-PCR analysis was performed to verify the role of developmental priming genes in salt tolerance shown in clv3-2. As a result, some genes involved in cell division showed similar accumulation patterns in the shoot apex of Ler and clv3-2, respectively ( Figure 7A). In contrast, some genes classified as abiotic/biotic stress responses displayed abnormal responses only in clv3-2 under salt stress conditions ( Figure 7B). These results suggest that such genes involved in cell division may function in the recovery step of cell division after removing the stress stimuli instead of salt tolerance for plastic growth and development under stress conditions. In addition, they indicate that such genes involved in stress responses may play an important role in salt tolerance shown in clv3-2.
with cell division and microtubule-based movement were included as upregulated priming genes in the shoot apex compared to other leaf tissue. By contrast, the salt treatment reduced the expression change of these genes. On the other hand, the expression of abiotic/biotic-stress-responsive genes associated with downregulated priming genes in the shoot apex under normal conditions was induced by salt stress conditions. These results revealed the functional behavior of developmental priming genes depending on salt stresses in the WT SAM. Because previous studies showed that the clv3-2 mutant exhibiting the increase in cell division or the decrease in cell differentiation has the property of salt resistance compared to WT Ler via the non-canonical responses [15], qRT-PCR analysis was performed to verify the role of developmental priming genes in salt tolerance shown in clv3-2. As a result, some genes involved in cell division showed similar accumulation patterns in the shoot apex of Ler and clv3-2, respectively ( Figure 7A). In contrast, some genes classified as abiotic/biotic stress responses displayed abnormal responses only in clv3-2 under salt stress conditions ( Figure 7B). These results suggest that such genes involved in cell division may function in the recovery step of cell division after removing the stress stimuli instead of salt tolerance for plastic growth and development under stress conditions. In addition, they indicate that such genes involved in stress responses may play an important role in salt tolerance shown in clv3-2.

Discussion
In this study, we provided novel information on developmental priming genes associated with growth plasticity under unfavorable stress conditions using RNA-seq experiments and analysis. Developmental priming genes differentially expressed in the shoot apex under both control and salt stress conditions compared to leaf tissue contained ap- Figure 7. Expression patterns of developmental priming genes in WT Ler and clv3-2 shoot apices without (−) or with (+) salt stress conditions. Developmental priming genes related to cell division (A) and abiotic/biotic stress responses (B) were examined by qRT-PCR and ACT2 genes were used to normalize each gene expression. The data are presented as the mean ± SD (n = 3~5). The different letters indicate a significant difference (p < 0.05) according to the one-way ANOVA with post-hoc Turkey's HSD test.

Discussion
In this study, we provided novel information on developmental priming genes associated with growth plasticity under unfavorable stress conditions using RNA-seq experiments and analysis. Developmental priming genes differentially expressed in the shoot apex under both control and salt stress conditions compared to leaf tissue contained approximately 10% TFs, including two bHLHs (EGL3, At1g49830), three homeodomaincontaining proteins (STM, ATHB13, ATHB33), two C2H2 zinc finger proteins (SGR5, ZFP2), two MYB domain proteins (AS1/ATPHAN, MYB117), two WRKYs (WRKY6, WRKY35), an AP2 protein (ANT), and a basic leucine zipper (ATbZIP42) ( Table 1). These types of TFs, including MYB, HD-ZIP, AP2, and WRKY, are known to mediate the transcriptional regulations in response to salt [53]. By clustering analysis of these TF genes, ARABIDOP-SIS THALIANA HOMEOBOX PROTEIN 33 (ATHB33)/ZINC-FIGNER HOMEODOMAIN 5 (ZHD5), STM, ATHB13, ASYMMETRIC LEAVES 1 (AS1)/PHANTASTICA (ATPHAN), and AINTEGUMENTA (ANT) genes were relatively related ( Figure S6). The ATHB33 gene has been reported to function in seed germination and root growth via negative regulation by ABA-induced AUXIN RESPONSE FACTOR 2 (ARF2), and the overexpression of ATHB33/ZHD5 displayed accelerated plant growth, including floral architecture and leaf [25,26]. A recent study showed that STM expression is downregulated by the repressive class B ARFs, such as ETTIN (ETT/ARF3) and ARF4, to promote the initiation of flowers [21]. In addition to the developmental function of STM, recent reports showed that transgenic plants overexpressing STM adapt to drought stress via ABA-inducible MYB96, and the mechanical stress caused by curvature in the boundary domain of the shoot meristem induces STM expression correlated with auxin depletion [22,23]. ATHB13 is also known to function in plant development of stem elongation, silique, seed production, and ovules via the regulation of cell division [18]. Because the activation of ATHB13 showed disease resistance to powdery mildew species and tolerance to drought and salinity [19,20], these results may provide the role of ATHB13 in the crosstalk among development and abiotic/biotic stress resistance. Interestingly, SGR5 encoding zinc finger protein displayed the plastic adaptation role via the negative gravitropism of inflorescence stem under heat stress conditions [27]. In addition, AS1/ATPHAN, ANT, and WRKY6 also showed multi-functional roles in development and abiotic/biotic stress responses (Table 1). Therefore, these various previous reports indicate the role of phytohormones, such as auxin and ABA for plant adaptation in the SAM, and suggest that the developmental priming genes, including TFs derived from RNA-seq experiments, may be reasonable candidates for the development and stress responses at the SAM.
After sensing salt stresses, plants can be placed in early and late responses. Immediate early responses are composed of calcium burst, ROS production, and activities of various pump/ion transporters, which play an essential role in ion homeostasis and sodium sequestering [19]. After early responses, the growth rate of plants changes toward the quiescent phase via the regulation of ABA levels and signaling, and the growth is then recovered through the raised levels of jasmonate (JA), brassinosteroid (BR), and gibberellic acid (GA) for homeostasis [54,55]. By salt application, microtubules are broken down in the early response phase, and the microtubule network is rebuilt dynamically during the quiescent phase [54,56]. Interestingly, in addition to TFs, other candidates as developmental priming genes were isolated in this study, and those putative functions based on GO enrichment analysis showed development-related biosynthesis and regulation of auxin, BR, and nitrile, as well as cell division and microtubule-based movement in the upregulated priming genes at the shoot apex ( Figure 3A and Table S4), whereas most downregulated priming genes contained ROS-related catalytic enzymes based on GO terms of molecular function ( Figure 3B and Table S5). Overall, these suggest that the putative developmental priming genes isolated from RNA-seq experiments probably represent early and late responsive genes after salt application.
A previous study of cell division in shoot meristems using EdU staining showed a decreased rate of cell division by a salt treatment and recovered cell division after removing salt stresses [15]. Because salt stress triggers defective cell division, we could also find various cell-division-related genes in the group of priming and salt-inducible DEGs in the shoot apex ( Figure 5C and Table S12). Interestingly, most genes classified as cell division in Table S12 function in the process between metaphase and anaphase via chromosome segregation during cell division. TPX2 (targeting protein for Xklp2) has been known to act as an activator of Aurora A kinase, which is one component of the spindle assembly checkpoint (SAC) complex delaying chromosome segregation in animals by sequestering CELL DIVISION CYCLE 20 (CDC20) [57,58]. In addition, TPX2 colocalized with AtAurora1 and γ-tubulin on spindle microtubules to control the meristematic cell proliferation in Arabidopsis [59]. CDC20.1, as a cofactor of the anaphase-promoting complex (APC), regulates Aurora localization and meiotic chromosome segregation, and the cdc20.1 mutant showed defects of meiosis and male fertility [60]. In another salt priming candidate gene CYC1BAT, B-type cyclin and cdk1 form the mitosis promoting factor (MPF) complex phosphorylating TXP2 to activate Aurora kinase [61]. TOPII (topoisomerase II) is also important for chromosome stability, including chromosome segregation for DNA repair and mitosis/meiosis in eukaryotes [62,63]. However, since the expression patterns of these genes involved in the process of cell division under control and salt stress conditions in WT Ler and clv3-2 were almost similar ( Figure 7A), the developmental priming genes related to cell division are probably not critical for the direct response against salt stimuli in clv3-2 displaying salt tolerance compared to WT [15]. Instead, it seems to be important for resilience in plastic plant growth after salt application. Indeed, previous studies showed that the reduced rate of cell division under salt stress conditions in the SAM was almost recovered after removing the salt treatment [15]. Future studies will be needed to determine the biological function of these priming genes in regulating shoot development and stress responses in planta.

Plant Materials and Growth Conditions
The Arabidopsis thaliana accessions such as Columbia-0 (Col-0) and Landsberg erecta (Ler) were used in this study as the wild-type. The genetic crossing for the clv1 bam1 double mutant was reported previously [15]. For the analysis of survival rate under salt stresses, Ler, clv3-2, Col-0, and clv1 bam1 seeds were sowed on half MS plates containing 0.5 × MS, 0.5% Sucrose and 0.8% phytoagar without or with 100 or 150 mM NaCl, and seedlings in half MS plates were grown at 22-23 • C with 75 µmol m −2 s −1 light intensity under 16 h light/8 h dark photoperiod condition for 3 weeks. For RNA-seq and qRT-PCR experiments, Ler and clv3-2 seedlings were grown vertically on half MS plates for 9 days. Then, seedlings showing similar growth status were transferred to new half MS plates without or with 200 mM NaCl and grown for 1 day without any morphological change.

RNA Extraction and RNA Sequencing
Seedlings treated without or with NaCl treatment were dissected into shoot apex and leaves. Three biological replicates at each tissue with corresponding condition were used for RNA extraction, which was carried out with TRI Reagent (Molecular Research Center, Cincinnati, OH, USA) according to the manufacturer's instructions. The libraries were prepared for 150 bp paired-end sequencing using TruSeq stranded mRNA Sample Preparation Kit (Illumina, CA, USA). mRNA molecules were purified and fragmented from 1 µg of total RNA using oligo (dT) magnetic beads. The fragmented mRNAs were synthesized as single-stranded cDNAs through random hexamer priming. By applying this as a template for second strand synthesis, double-stranded cDNA was prepared. After sequential process of end repair, A-tailing, and adapter ligation, cDNA libraries were amplified with PCR (Polymerase Chain Reaction). The quality of these cDNA libraries was evaluated with the Agilent 2100 BioAnalyzer (Agilent, CA, USA). They were quantified with the KAPA library quantification kit (Kapa Biosystems, MA, USA) according to the manufacturer's library quantification protocol. Following cluster amplification of denatured templates, sequencing was progressed as paired-end (2 × 150 bp) using Illumina NovaSeq6000 (Illumina, San Diego, CA, USA). Raw data of RNA-seq reads obtained in this study were deposited with the accession nos. SRR20998325−SRR20998336 (BioProject: PRJNA868122) in the NCBI Short Read Archive (SRA) database.

Transcriptome Data Analysis
The adapter sequences and the ends of the reads less than Phred quality score 20 were trimmed, and simultaneously, reads shorter than 50 bp were removed by using cutadapt v.2.8 [64]. Filtered reads were mapped to the reference genome related to the species using the aligner STAR v.2.7.1a [65] following ENCODE standard options (refer to "Alignment" of "Help" section in the html report) with "-quantMode TranscriptomeSAM" option for estimation of transcriptome expression level. Gene expression estimation was performed by RSEM v.1.3.1 [66] considering the direction of the reads which are corresponding to the library protocol using option "-strandedness". To improve the accuracy of the measurement, "-estimate-rspd" option was applied. All other options were set to default values. To normalize sequencing depth among samples, FPKM (fragments per kilobase of transcript per million) and TPM (transcripts per millions) values were calculated. Based on the estimated read counts in the previous step, DEGs were identified using the R package called TCC v.1.26.0 [67]. TCC package applies robust normalization strategies to compare tag count data. Normalization factors were calculated using the iterative DESeq2/edgeR method [68,69]. Q-value was calculated based on the p-value using the p.adjust function of R package with default parameter settings. The DEGs were identified based on the q value threshold less than 0.05 for correcting errors caused by multiple-testing [70]. DEGs were visualized by volcano, scatter, and MA plots. Heatmap and clustering analysis using the value of FPKM or the log2 ratio of fold change were visualized by the MeV 4.9.0 software (https://mev.tm4.org/#/about, accessed on 22 September 2022). The analysis of Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) for developmental priming genes or shoot apex-specific differentially expressed genes was performed using the DAVID software (https://david.ncifcrf.gov/, accessed on 22 September 2022), which is based on the database of Arabidopsis TAIR10 release.

cDNA Synthesis and qRT-PCR
First-stranded cDNAs from 1 µg of total RNA were synthesized using the ImProm-II TM reverse transcriptase (Promega, Madison, WI, USA) according to the manufacturer's instructions. Quantitative reverse transcription-polymerization chain reaction (qRT-qPCR) assays were performed with primers of developmental priming genes identified from RNA-seq experiments (Table S13). Each reaction was carried out using MIC qPCR cycler (Bio Molecular Systems, Upper Coomera, QLD, Australia) with SYBR Green Realtime PCR Master Mix (TOYOBO, Osaka, Japan).

Statistical Analyses
All quantitative results were analyzed by one-way analysis of variation (ANOVA) using SPSS (ver. 27) statistic software. Significant differences and different letters were determined by one-way ANOVA with post-hoc Turkey's HSD multiple comparison tests.

Conclusions
This study provides a fundamental resource for understanding the developmental priming effect in shoot meristems orchestrating the aboveground growth. Under normal (non-stress) conditions, the putative priming genes regulating phytohormones, development, and abiotic/biotic stress responses are expressed differently in the shoot apex than in leaf tissues. Because many stress-related genes are highly or lowly pre-accumulated in the shoot apex, the transcriptome data also support the idea of a primed state in the SAM to respond efficiently to incoming severe stress conditions [12,13]. In addition, some cell division and stress-responsive genes are probably associated with a wide range of salt stress responses for plastic growth and development under salt stress conditions. Because the defective cell division within the shoot meristems affects overall growth and development, these developmental priming candidate genes isolated from RNA-seq experiments using enriched shoot apex may provide new insight into understanding plastic growth adaptation under salt stress conditions. Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11192546/s1, Figure S1: Quality validation of RNAseq experiments using leaf and shoot apex tissues under control conditions; Figure S2: Quality validation of RNA-seq experiments using leaf and shoot apex tissues under salt stress conditions; Figure S3: Upregulated developmental priming genes in a KEGG pathway; Figure S4: Quality validation of RNA-seq experiments using shoot apex tissues under control or NaCl conditions; Figure  S5: Quality validation of RNA-seq experiments using leaf tissues under control or NaCl conditions; Figure S6: Clustering analysis of TF genes classified as developmental priming genes; Table S1: Developmental priming genes in the shoot apex compared to leaf under the control condition; Table  S2: Developmental priming genes in the shoot apex compared to leaf under the NaCl condition; Table S3: Developmental priming genes in the shoot apex under both control and NaCl conditions; Table S4: GO enrichment analysis of upregulated priming genes; Table S5: GO enrichment analysis of downregulated priming genes; Table S6: List of removed 14 genes from clustering analysis; Table S7: Fold change values of developmental priming genes in the shoot apex used in clustering analysis; Table S8: Salt-inducible DEGs in the shoot apex; Table S9: Salt-inducible DEGs in the leaf; Table  S10: GO enrichment (BP, biological process) analysis of salt-inducible DEGs in the shoot apex, leaf, and both tissues; Table S11: FRKM values of 117 developmental priming genes in the shoot apex;   Data Availability Statement: RNA-seq data analyzed in this study are available in the NCBI SRA database (SRR20998325−SRR20998336).

Conflicts of Interest:
The authors declare no conflict of interest.