Comprehensive Genome-Wide Analysis and Expression Pattern Profiling of PLATZ Gene Family Members in Solanum Lycopersicum L. under Multiple Abiotic Stresses

PLATZ (plant AT-rich sequence and zinc-binding) family proteins with two conserved zinc-dependent DNA-binding motifs are transcription factors specific to the plant kingdom. The functions of PLATZ proteins in growth, development, and adaptation to multiple abiotic stresses have been investigated in various plant species, but their role in tomato has not been explored yet. In the present work, 20 non-redundant Solanum lycopersicum PLATZ (SlPLATZ) genes with three segmentally duplicated gene pairs and four tandemly duplicated gene pairs were identified on eight tomato chromosomes. The comparative modeling and gene ontology (GO) annotations of tomato PLATZ proteins indicated their probable roles in defense response, transcriptional regulation, and protein metabolic processes as well as their binding affinity for various ligands, including nucleic acids, peptides, and zinc. SlPLATZ10 and SlPLATZ17 were only expressed in 1 cm fruits and flowers, respectively, indicating their preferential involvement in the development of these organs. The expression of SlPLATZ1, SlPLATZ12, and SlPLATZ19 was up- or down-regulated following exposure to various abiotic stresses, whereas that of SlPLATZ11 was induced under temperature stresses (i.e., cold and heat stress), revealing their probable function in the abiotic stress tolerance of tomato. Weighted gene co-expression network analysis corroborated the aforementioned findings by spotlighting the co-expression of several stress-associated genes with SlPLATZ genes. Confocal fluorescence microscopy revealed the localization of SlPLATZ–GFP fusion proteins in the nucleus, hinting at their functions as transcription factors. These findings provide a foundation for a better understanding of the structure and function of PLATZ genes and should assist in the selection of potential candidate genes involved in the development and abiotic stress adaptation in tomato.


Introduction
A transcription factor (TF) is a regulatory protein with one or more DNA-binding domains that can control the expression of target genes by binding directly to their promoter regions. From a combinatorial gene regulation perspective, TFs can also physically interact with each other to control diverse cellular processes [1]. The evolution of the TF families does not always demonstrate a robust corresponding correlation among plants, of Arabidopsis by repressing the expression of CBL4/SOS3 and CBL10/SCaBP8 [27]. The overexpression of cotton PLATZ1 in Arabidopsis enhances the osmotic and salt stress tolerance of the resulting transgenic lines at the germination and seedling stages [28], while that of soybean PLATZ1 in Arabidopsis shows delayed germination under osmotic stress [29]. Furthermore, several reports have demonstrated that PLATZ genes in multiple plant species are responsive to various hormones and abiotic stresses [28,[30][31][32].
Tomato (Solanum lycopersicum L.) is of commercial value due to its widespread production and consumption around the globe. As a result of the abundance of several nutrients, such as lycopene, lutein, Zeaxanthin, potassium, ascorbic acid, and β-carotene in tomato, it is also the second most significant vegetable crop in the world [33]. With the availability of whole genome sequences of tomato cultivars as well as with its numerous interesting features, such as compound leaves, a sympodial shoot, and fleshy fruit, which cannot be analyzed in other model plants (i.e., Arabidopsis and rice), tomato has also been extensively studied as a model organism for fleshy-fruited plants [34]. Genome-wide analyses utilize various bioinformatics tools to explore the characteristics of gene family members, gene structures, phylogenetic relationships, expression profiles, gene ontologies, putative gene functions, and so forth. The data extracted from genome-wide identification studies provide researchers with valuable information and an initial framework for further functional verification of potential candidate genes that can be exploited for plant genetic improvement [35]. A genome-wide characterization of the PLATZ gene family in diverse plant species has been conducted [14,32,[36][37][38][39][40], but there has been no systematic investigation of their phylogenetic relationships or characteristics in solanaceous crops. The productivity and fruit quality of tomato are severely affected by various adverse environmental stresses, and hence it is essential to identify the potential novel genes implicated in the abiotic stress responses of tomato for the utilization of these genes in the genetic improvement of tomato [41]. Therefore, we conducted expression profiling of tomato PLATZ genes at various developmental stages and on exposure to diverse abiotic stresses. In addition to the genome-wide characterization of SlPLATZ proteins, their subcellular localization assays and gene co-expression network analysis using RNA sequencing data were carried out to shed light on their probable functional role in tomato. This study aimed to provide a solid foundation for further functional verification of potential PLATZ transcription factors for the improvement of abiotic stress adaptation in tomato.

Identification and Sequence Analysis of PLATZ Family Members in Tomato
A total of 20 non-redundant PLATZ genes were identified in the tomato genome and designated as SlPLATZ1-SlPLATZ20 in accordance with their chromosomal positions. In silico analyses revealed that the lengths of predicted PLATZ proteins varied from 137 (SlPLATZ6) to 255 (SlPLATZ14) amino acids (aa), and their molecular weight ranged between 15.94 and 29.16 kDa for SlPLATZ6 and SlPLATZ14, respectively. The theoretical pI value of only one protein (SlPLATZ6) was less than 7, whereas those of the other 19 proteins were greater than 7, revealing that all tomato PLATZ family proteins were basic proteins, except SlPLATZ6. The GRAVY of SlPLATZ proteins ranged between −0.217 and −0.709, showing that all these proteins were hydrophilic. Further information on SlPLATZ family proteins, including locus IDs, is shown in Table 1.

Multiple Protein Sequence Alignment and Phylogenetic Analysis
Multiple sequence alignment of tomato PLATZ proteins with a PLATZ protein from the model plant Arabidopsis underlined the evolutionary conservation of PLATZ proteins from tomato and Arabidopsis (Figure 1), particularly the presence of the two conserved regions, the first region harboring five cysteine and three histidine residues [C-x 2 -H-x 11 -C-x 2 -C-x( 4-5 )-C-x 2 -C-x( 3-5 )-H-x 2 -H] and the second region with four cysteine residues [C-x 2 -C-x( 10-12 )-C-x 3 -C], crucial for zinc-dependent DNA binding [38,39]. However, it is interesting to note that SlPLATZ6 is the only protein devoid of the second conserved region. To investigate the evolutionary relationship of the SlPLATZ gene family, a phylogenetic tree was constructed using PLATZ proteins from 10 different plant species. (Figure 2). PLATZ family members from various plant species were divided into two primary clades, which can be further categorized into seven groups based on their phylogenetic affinity. The evolutionary tree constructed indicated a major grouping of homologous proteins from monocots (maize, rice, and sorghum) and dicots (Brassica rapa, tomato, Arabidopsis, and potato) together with lower plants. Intriguingly, PLATZ proteins from the green alga Chlamydomonas reinhardtii were distributed specifically in groups 2 and 3 of clade I and clustered with the homologs from monocots, dicots, and the basal angiosperm Amborella trichopoda. The paralogous proteins from the moss Physcomitrella patens were clustered together in groups 1 and 3, but PpPLATZ13 was paired with the homolog from A. trichopoda in group 6. The basal angiosperm homologs, distributed in groups 1, 3, 5, 6, and 7, were clustered not only with those from dicots and monocots but also with those from the lower plants. The 20 tomato PLATZ proteins were allocated to four of the seven groups, except groups 3, 5, and 6. The greatest number (12) of SlPLATZs was arranged in group 2, while the smallest number (1) was in group 4. Groups 1 and 7 included four and three SlPLATZ proteins, respectively. PLATZ proteins from tomato were grouped specifically with those of their evolutionarily closely related species, potato.  (Figure 3). Among the three proteins having an additional domain, SlPLATZ4 was the only member harboring a TM domain in the rear of the PLATZ domain, while a BBOX domain positioned before the PLATZ domain was predicted in SlPLATZ15 and SlPLATZ16. Analysis of the exon-intron distribution of the SlPLATZ genes indicated that there was no gene lacking an intron in their coding sequences ( Figure S1). Except for SlPLATZ5, which contained only one intron, all SlPLATZ genes have between two and three introns. Intriguingly, the paralogous gene pairs in the same clades, such as SlPLATZ15 and SlPLATZ16, shared the identical gene structure in terms of intron number and exon length. in group 6. The basal angiosperm homologs, distributed in groups 1, 3, 5, 6, and 7, were clustered not only with those from dicots and monocots but also with those from the lower plants. The 20 tomato PLATZ proteins were allocated to four of the seven groups, except groups 3, 5, and 6. The greatest number (12) of SlPLATZs was arranged in group 2, while the smallest number (1) was in group 4. Groups 1 and 7 included four and three SlPLATZ proteins, respectively. PLATZ proteins from tomato were grouped specifically with those of their evolutionarily closely related species, potato.
-H] and the blue box marks the second conserved region [C-x2-C-x(10-12)-C-x3-C], which are essential for zinc-dependent DNA binding. The yellow underline indicates the location of the PLATZ domain. The elements of the secondary structure determined using the ESPript 3.0 web server are represented by gray bars (α-helices) and blue-green arrows (β-strands) above the alignment.  (3)(4)(5) -H-x 2 -H] and the blue box marks the second conserved region [C-x 2 -C-x (10-12) -C-x 3 -C], which are essential for zinc-dependent DNA binding. The yellow underline indicates the location of the PLATZ domain. The elements of the secondary structure determined using the ESPript 3.0 web server are represented by gray bars (α-helices) and blue-green arrows (β-strands) above the alignment.  (Figure 3). Among the three proteins having an additional do-  To further analyze the structural diversity of the PLATZ gene family, we explored the composition of conserved motifs in SlPLATZ proteins with respect to the PLATZ homologous proteins from Arabidopsis and rice ( Figures S2 and S3). The paralogous gene pairs and orthologous gene pairs from the three plant species harbored almost identical motif compositions. The motif organization of most SlPLATZ proteins belonging to phylogenetic group 2 is almost identical, with the exception of SlPLATZ5, which contained only motifs 1, 3, and 6. Motifs 1 and 5, which resided in the PLATZ domain region, were the most prevalent motifs in nearly all the PLATZ homologs of three plant species, except for AtPLATZ9 and OsPLATZ15, which were devoid of motif 1 and 5, respectively. Intriguingly, motif 7 was uniquely present in most SlPLATZ protein members belonging to group 2, but not identified in any PLATZ homologs from rice and Arabidopsis. In addition, motif 9 was specific to certain PLATZ proteins from all plant species belonging to group 1.

Chromosomal Distribution, Gene Duplication, and Microsynteny Analysis of SlPLATZ Genes
SlPLATZ genes were unequally dispersed on eight of the twelve chromosomes (Chr), except for Chr03, Chr05, Chr09, and Chr12 ( Figure S4). The majority of tomato chromosomes harbored a small number of PLATZ genes (1-3) residing near the distal end of the chromosomes, while Chr02 possessed 10 genes located in the central region of the chromosome. Of the chromosomes with a lesser number of PLATZ genes, Chr07 had three To further analyze the structural diversity of the PLATZ gene family, we explored the composition of conserved motifs in SlPLATZ proteins with respect to the PLATZ homologous proteins from Arabidopsis and rice ( Figures S2 and S3). The paralogous gene pairs and orthologous gene pairs from the three plant species harbored almost identical motif compositions. The motif organization of most SlPLATZ proteins belonging to phylogenetic group 2 is almost identical, with the exception of SlPLATZ5, which contained only motifs 1, 3, and 6. Motifs 1 and 5, which resided in the PLATZ domain region, were the most prevalent motifs in nearly all the PLATZ homologs of three plant species, except for AtPLATZ9 and OsPLATZ15, which were devoid of motif 1 and 5, respectively. Intriguingly, motif 7 was uniquely present in most SlPLATZ protein members belonging to group 2, but not identified in any PLATZ homologs from rice and Arabidopsis. In addition, motif 9 was specific to certain PLATZ proteins from all plant species belonging to group 1.

Chromosomal Distribution, Gene Duplication, and Microsynteny Analysis of SlPLATZ Genes
SlPLATZ genes were unequally dispersed on eight of the twelve chromosomes (Chr), except for Chr03, Chr05, Chr09, and Chr12 ( Figure S4). The majority of tomato chromosomes harbored a small number of PLATZ genes (1-3) residing near the distal end of the chromosomes, while Chr02 possessed 10 genes located in the central region of the chromosome. Of the chromosomes with a lesser number of PLATZ genes, Chr07 had three genes, which was followed by Chr08 with two genes. The remaining five chromosomes (Chr01, Chr04, Chr06, Chr10, and Chr11) carried only one gene each.
A total of seven duplicated gene pairs were identified in the tomato PLATZ gene family ( Table 2, Figures S4 and S5). Among the 20 SlPLATZ genes, four gene pairs (SlPLATZ2/SlPLATZ3, SlPLATZ3/SlPLATZ4, SlPLATZ6/SlPLATZ7, and SlPLATZ15/SlPLATZ16) were predicted to be tandemly duplicated, while three gene pairs (SlPLATZ13/SlPLATZ18, SlPLATZ14/SlPLATZ20, and SlPLATZ17/SlPLATZ18) were detected as segmental duplicates. The duplicated genes of each pair were located in the same phylogenetic group. The mode of selection pressure imposed on the SlPLATZ gene duplicates was identified based on Ka/Ks values. The Ka/Ks ratios of all SlPLATZ gene duplicates were lower than 1, revealing that these genes had undergone strong purifying/negative selection with little variation after duplication. The computation of the Ka/Ks ratios of gene duplicates indicated that the duplication events were predicted to have occurred between 1.05 and 31.43 million years ago (MYA). These duplicated gene pairs may have diverged in the last 1-31 million years. Furthermore, a comparative microsyntenic analysis was performed to investigate the phylogenetic relationship between tomato, Arabidopsis, and rice. No orthologous gene pairs were predicted between rice and Arabidopsis, whereas we identified seven orthologs between tomato and Arabidopsis and two orthologs between tomato and rice, respectively ( Figure 4). SlPLATZ18 was orthologous to three PLATZ genes (AtPLATZ3, AtPLATZ11, and AtPLATZ12) from Arabidopsis, while SlPLATZ12 and SlPLATZ17 were orthologous to AtPLATZ6 and AtPLATZ11, respectively. Intriguingly, SlPLATZ1 and SlPLATZ11 were collinear with the PLATZ genes from both Arabidopsis and rice. SlPLATZ1 and SlPLATZ11 were orthologous to AtPLATZ10 and OsPLATZ6, and OsPLATZ11 and SlPLATZ2, respectively.

Analysis of Cis-Acting Elements and miRNA Target Sites in SlPLATZ Genes
We analyzed the cis-regulatory elements in the upstream regions of SlPLATZ gene family members to explore their probable involvement in the transcriptional regulation on exposure to diverse abiotic stresses. Different numbers of cis-acting elements associated with phytohormones and abiotic stress adaptation were predicted in the promoter sequences of tomato PLATZ genes ( Figure S6 and Table S4). The identified stress-responsive cis-elements were LTR (related to low temperature and hypersalinity stress and defense), TC-rich repeats (involved in defense and stress responses), drought-responsive MYBbinding site (MBS), and WUN motifs (implicated in wounding responses), whereas those associated with hormones were TGA elements and AuxRR-core (related to the auxin response), ABA-responsive ABRE elements, TGA elements and AuxRR-core (with roles in the auxin response), TCA-elements (implicated in the SA response), the TGACG-and CGTCA-motifs (involved in the jasmonic acid [JA] response), and gibberellic acid (GA) response-related TATC-and P-boxes.

Analysis of Cis-Acting Elements and miRNA Target Sites in SlPLATZ Genes
We analyzed the cis-regulatory elements in the upstream regions of SlPLATZ gene family members to explore their probable involvement in the transcriptional regulation on exposure to diverse abiotic stresses. Different numbers of cis-acting elements associated with phytohormones and abiotic stress adaptation were predicted in the promoter sequences of tomato PLATZ genes ( Figure S6 and Table S4). The identified stress-responsive cis-elements were LTR (related to low temperature and hypersalinity stress and defense), TC-rich repeats (involved in defense and stress responses), drought-responsive MYB-binding site (MBS), and WUN motifs (implicated in wounding responses), whereas those associated with hormones were TGA elements and AuxRR-core (related to the auxin response), ABA-responsive ABRE elements, TGA elements and AuxRR-core (with roles in the auxin response), TCA-elements (implicated in the SA response), the TGACG-and Transcription factors (TFs) and microRNAs (miRNAs) are crucial gene regulatory factors, and miRNA-mediated regulation of transcription factors plays a vital role in plant development and stress tolerance [42]. Thus, we examined the miRNA target sequences in SlPLATZ genes to shed light on the relationship between miRNA and tomato PLATZ transcription factors in environmental stress response. A total of 12 out of 20 tomato PLATZ genes were predicted to be the target genes of diverse tomato miRNAs, including sly-miR159, sly-miR162, sly-miR171a, sly-miR319a, sly-miR319b, sly-miR319c-3p, sly-miR396a-5p, sly-miR396b, sly-miR403-5p, sly-miR482a, sly-miR5303, sly-miR6024, sly-miR9469-3p, sly-miR9472-3p, and sly-miR9478-3p, which are involved in the development and stress response of tomato (Table S5).

Analysis of Three-Dimensional Structures of Predicted SlPLATZ Proteins
The 3D model structures of tomato PLATZ proteins predicted through the I-TASSER server showed creditable values of parameters used in comparative modeling ( Figure 5), including the percentage of sequence identity, coverage, and Z-scores (Tables S6 and S7). This indicated the reliability of the predicted SlPLATZ models. The binding sites in the resulting 3D models were determined with Discovery Studio v.21.1, and the secondary structural components predicted in the models ranged from 2-8 for α-helices, 5-8 for β-strands, and 8-15 for coils, respectively. The majority of duplicated gene pairs in the SlPLATZ gene family displayed similar secondary structural components, but SlPLATZ3 and SlPLATZ4 had variable numbers of β-strands and coils, which were 5-8 and 10-12, respectively (Table S8). The prediction of the binding residues based on an alignment between the template and the generated SlPLATZ models revealed that five SlPLATZ proteins (SlPLATZ1, SlPLATZ2, SlPLATZ5, SlPLATZ13, and SlPLATZ17) showed binding affinity to nucleic acids, while eight proteins (SlPLATZ3, SlPLATZ4, SlPLATZ6, SlPLATZ7, SlPLATZ9, SlPLATZ15, SlPLATZ16, and SlPLATZ19) displayed binding affinity to zinc. SlPLATZ18 The prediction of the binding residues based on an alignment between the template and the generated SlPLATZ models revealed that five SlPLATZ proteins (SlPLATZ1, SlPLATZ2, SlPLATZ5, SlPLATZ13, and SlPLATZ17) showed binding affinity to nucleic acids, while eight proteins (SlPLATZ3, SlPLATZ4, SlPLATZ6, SlPLATZ7, SlPLATZ9, SlPLATZ15, SlPLATZ16, and SlPLATZ19) displayed binding affinity to zinc. SlPLATZ18 was the only protein that showed an affinity to both nucleic acids and zinc. The remaining six proteins that lacked binding affinity to both nucleic acids and zinc (SlPLATZ8, SlPLATZ10, SlPLATZ11, SlPLATZ12, SlPLATZ14, and SlPLATZ20) were predicted to bind to other diverse ligands such as peptide, mannose, copper, calcium, magnesium, cyanide, phosphate, glutathione, chlorophyll-a, and glycine (Table S9).
Gene ontology (GO) terms predicted from the modeling of SlPLATZ proteins by the I-TASSER web tool suggested that the modeled SlPLATZ proteins might function in diverse biological processes, such as cellular defense response, transcriptional regulation, gene expression, translation, protein maturation, and protein metabolic processes (Table S10).

Expression Profiling of Tomato PLATZ Genes in Various Organs
To further explore the possible functional role of SlPLATZ genes in developmental processes, we studied their expression patterns in different tomato tissues (leaves, stems, roots, flowers, 1 cm fruits, IM fruits, MG fruits, B fruits, and B5 fruits) through RT−qPCR analysis in the tomato cultivar Ailsa Craig. The higher homology in the cDNA sequences of most SlPLATZ genes inhibits the availability of the sequence regions to design the primers for the expression analysis ( Figure S7). Thus, nine tomato PLATZ genes, viz., SlPLATZ1, SlPLATZ5, SlPLATZ10, SlPLATZ11, SlPLATZ12, SlPLATZ14, SlPLATZ17, SlPLATZ19, and SlPLATZ20, were studied for expression profiling in various tomato organs. The expression of SlPLATZ5, SlPLATZ14, and SlPLATZ20 was not determined in any of the organs studied. Among the remaining six SlPLATZ genes, two genes (SlPLATZ10 and SlPLATZ17) displayed tissue-specific expression, while the other genes (i.e., SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19) exhibited varied transcript levels across the different organs tested ( Figure 6).
The transcripts of SlPLATZ10 were detected solely in 1 cm fruits but not in any other organs examined. Likewise, the expression of SlPLATZ17 was determined in no organ, except in flowers. SlPlATZ1 was expressed in all tomato organs, except in breaker 5 fruits. Relative to the control (leaves), SlPlATZ1 showed peak expression in 1 cm fruits (~12-fold), which was followed by flowers (~10-fold) and immature green fruits (~6-fold), respectively. It is noteworthy that the expression of SlPlATZ1 was higher in the early developmental stages of tomato fruits, and then gradually diminished at the late developmental stages. The transcript levels of SlPlATZ11 were higher (>6-fold) in flowers, roots, and stems compared with the control. However, its expression was downregulated (2.7-fold relative to the control) in 1 cm fruits and was not determined at all in the remaining developmental stages of tomato fruits.
The expression of SlPlATZ12 was detected in various organs tested, and it showed the highest expression (>15-fold) in stems, followed by roots (>7-fold), and 1 cm fruits (2.6-fold), respectively, compared to the control. However, in contrast to the control, the transcript levels of SlPlATZ12 were downregulated in the other developmental stages of fruits (i.e., IM fruits (>2-fold), MG fruits (>1.6-fold), breaker fruits (>3.6-fold), and breaker 5 fruits (>8.8-fold)). SlPlATZ19 showed the highest transcription levels in flowers (>1.5-fold over the control). Its expression is also higher in the vegetative organs (i.e., leaves, stems, and roots) compared with the five fruit developmental stages. The transcripts of SlPlATZ19 were reduced in 1 cm fruits (2.5-fold), IM fruits (2.7-fold), MG fruits (4.4-fold), breaker fruits (~2-fold) and breaker 5 fruits (4.5-fold) in comparison with the control. Plants 2022, 11, x 13 of 30 Error bars indicate the standard deviations of the means of three independent biological replicates. *, **, and *** represent the significant difference between the control samples (leaves) and the samples from the other organs according to Student's t−test, at p−values <0.05, <0.01, and <0.001, respectively.

Expression Analysis of SlPLATZ Genes in Response to Abiotic Stresses and Phytohormone Treatment
To explicate the probable functions of SlPLATZ genes in response to various abiotic stresses and the stress hormone ABA, the expression profiles of SlPLATZ genes

Expression Analysis of SlPLATZ Genes in Response to Abiotic Stresses and Phytohormone Treatment
To explicate the probable functions of SlPLATZ genes in response to various abiotic stresses and the stress hormone ABA, the expression profiles of SlPLATZ genes (SlPLATZ1, SlPLATZ15, SlPLATZ10, SlPLATZ11, SlPLATZ12, SlPLATZ14, SlPLATZ17, SlPLATZ19, and SlPLATZ20) in the leaves sampled before and after exposure to different stress treatments were analyzed via qRT−PCR assay. SlPLATZ5, SlPLATZ10, SlPLATZ14, SlPLATZ17, and SlPLATZ20 did not display detectable expression in all leaves harvested from control, abiotic stress-, or ABA-treated plants, while the other tomato PLATZ genes (SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19) were differentially expressed under these treatments ( Figure 7A-E).  Figure 7A-E). SlPLATZ1, SlPLATZ12, and SlPLATZ19 displayed responses to drought treatment, while the transcript levels of SlPLATZ11 were not significantly altered at all time points of the stress period ( Figure 7A). In comparison with the control (0 h), the transcript levels of SlPLATZ1 were significantly induced by >2-fold at 24 h after drought treatment and peaked (>3-fold) at 48 h in response to drought stress. Although the expression of SlPLATZ1 was unchanged at 60 h upon exposure to drought conditions, it was upregulated (>2-fold vs. the control) again at the last stage of the drought period. Similarly, SlPLATZ19 was highly up-regulated (2.9-to 8.2-fold) at all stages of drought treatment SlPLATZ1, SlPLATZ12, and SlPLATZ19 displayed responses to drought treatment, while the transcript levels of SlPLATZ11 were not significantly altered at all time points of the stress period ( Figure 7A). In comparison with the control (0 h), the transcript levels of SlPLATZ1 were significantly induced by >2-fold at 24 h after drought treatment and peaked (>3-fold) at 48 h in response to drought stress. Although the expression of SlPLATZ1 was unchanged at 60 h upon exposure to drought conditions, it was upregulated (>2-fold vs. the control) again at the last stage of the drought period. Similarly, SlPLATZ19 was highly up-regulated (2.9-to 8.2-fold) at all stages of drought treatment relative to the control. By contrast, the mRNA transcript levels of SlPLATZ12 were considerably downregulated (>2-to >55-fold relative to the control) throughout the drought stress period.
Under salt stress, the transcripts of SlPLATZ1 were initially repressed by~2-fold at 1 h after salt treatment compared with the control, but significantly upregulated (2.3-fold) at the last time point of the stress period. SlPLATZ12 expression was induced (>1.6-fold over the control) at 3 h of exposure to saline conditions but decreased from 1.5-to 1.8-fold at the later stages of salt treatment. Salinity stress attenuated the expression of SlPLATZ11 and SlPLATZ19 from 2.4 to 3-fold at 1 h and 9 h of the stress period, respectively ( Figure 7B).
In response to low-temperature conditions, SlPLATZ11 showed upregulation from 1.6-to 2-fold relative to the control at all stages of cold treatment. Likewise, the expression levels of SlPLATZ12 were induced by 1.6-to 2-fold throughout the cold stress period. By contrast, the expression of SlPLATZ1 was downregulated (1.9-to 2.6-fold over the control) at all time points following exposure to cold treatment. The transcripts of SlPLATZ19 were initially unchanged at 1 h following cold exposure but decreased significantly, 1.3-to 4-fold, at the later time points of low-temperature stress ( Figure 7C).
Under heat stress, the expression of SlPLATZ1 initially declined by~2-fold at 1 h of exposure to high-temperature conditions but recovered at 3 h and subsequently rose (2-to 4-fold) at 9 h and 24 h compared to the control. SlPLATZ11 was downregulated by 2.3-fold at 1 h after heat exposure but significantly upregulated (1.6-to 2-fold) at 3 h and 24 h. Similarly, the transcript accumulation levels of SlPLATZ12 were attenuated (2-fold vs. the control) at 1 h following exposure to a high temperature but rose 2.1-fold at 3 h. On the contrary, SlPLATZ19 was significantly repressed by 1.6-to 3.3-fold at most stages of heat treatment ( Figure 7D).
The expression levels of SlPLATZ1 were remarkably reduced (1.5-to 3.4-fold relative to the control) in the early phases of ABA treatment but significantly elevated (1.5-to 2-fold) at the later stages. SlPLATZ12 was repressed by >1.8-fold at 1 h following exposure to ABA stress, whereas SlPLATZ12 was minimally upregulated (1.3-fold) at 3 h. By contrast, the transcripts of SlPLATZ19 were not significantly altered upon exposure to exogenous ABA application ( Figure 7E).

Weighted Gene Co-Expression Network (WGCNA) Analysis and Functional Enrichment Analysis
WGCNA was conducted to construct the co-expression network of SlPLATZ genes using RNA-seq data (Figure 8). A total of 150 genes were detected in the co-expression networks of SlPLATZ genes. Specifically, 116, 10, and 24 genes were involved in the gene co-expression networks of SlPLATZ17, SlPLATZ18, and SlPLATZ19, respectively. The KEGG enrichment analysis of the co-expressed genes indicated that the co-expressed genes in the network of SlPLATZ genes were implicated in diverse biological pathways, including plant-pathogen interaction, RNA transport, mismatch repair, MAPK signaling pathway, biosynthesis of secondary metabolites, fatty acid biosynthesis, biosynthesis of amino acids, phenylpropanoid biosynthesis, nitrogen metabolism, protein export, and so forth (Figure 9, Tables S11 and S12). However, certain co-expressed genes have not been annotated in any biological process. Intriguingly, SlSHMT2 (Solyc02g091560.3) and a RIN-G finger gene (Solyc08g081370) responsive to multiple abiotic stresses, SnRK1 (Solyc02g067030) and SlZF-31 (Solyc08g063040) associated with abiotic stress tolerance of tomato, a NAC gene (Solyc04g005610) responsive to drought and salinity treatments, an abiotic stressinducible gene SlGS2 (Solyc01g080280), AHL5 (Solyc08g008030) implicated in the resistance to the oomycete pathogen Phytophthora capsica, and Glutaredoxin (Solyc08g062970) responsive to pathogen infection were co-expressed with the hub gene SlPLATZ19. SlGS2 (Solyc01g080280) inducible by various abiotic stresses, a stress-inducible gene coding for Chaperonin 21 (Solyc12g009250), and SnRK1 (Solyc02g067030) associated with stress tolerance in tomato were also present in the co-expression network of SlPLATZ18. Furthermore, SlPLATZ17 showed co-expression with several abiotic stress-associated genes, such as SlAPRR5 (Solyc03g081240) inducible by seawater stress, SlTLP6 (Solyc04g071750) responsive to osmotic stress, a gene encoding nucleolar protein 6 (Solyc02g085230) inducible by aluminum stress, and low-temperature inducible genes SlGRAS4 (Solyc01g100200) and SlPLDβ1(Solyc08g080130). We also identified several genes inducible by pathogen infec-tion in the co-expression network of SlPLATZ17, such as SlPP2C48 (Solyc06g076100), and genes that codify a GRAS transcription factor (Solyc07g063940), a major latex-like protein (Solyc08g023660), and a leucine-rich repeat receptor-like serine/threonine protein kinase (Solyc11g017270) (Figure 8).

Subcellular Location Analysis of SlPLATZ Proteins
The subcellular locations of SlPLATZ proteins predicted by the WoLF-PSORT server indicated that they were localized to the nucleus or cytoplasm (Table 1). For the further verification of their subcellular localization, full-length cDNAs of three SlPLATZ genes (SlPLATZ11, SlPLATZ12, and SlPLATZ19) were fused to the green fluorescent protein (GFP) and transiently expressed in rice protoplasts. Confocal imaging of protein fluorescence revealed that the fluorescence signals of SlPLATZ11 were observed not only in the nucleus but also in a large part of the cytoplasm. SlPLATZ12 was exclusively localized to the nucleus, while SlPLATZ19 was localized to the nucleus with its strong fluorescence signals detected in the cytosol (Figure 10).

Discussion
PLATZ genes, identified specifically in the plant kingdom, are novel zinc-dependent DNA-binding transcription factors that constitute a multi-gene family across a wide variety of plant species [14,32,[36][37][38][39][40]. In support of previous findings, we identified 20 PLATZ genes in the tomato genome, revealing their existence as a multiple-gene family in tomato. This finding suggested that PLATZ genes might play pivotal biological roles in several diverse plant species.
Our phylogenetic analysis divided the PLATZ proteins from various plant species into two major clades (Figure 2). Of these clades, only clade I contained the PLATZ homologs from the green alga Chlamydomonas reinhardtii, suggesting that the proteins from clade I, especially groups 2 and 3, may be more primitive compared with those from clade II, and they might have evolved before the divergence of chlorophytes and streptophytes over one billion years ago [43]. The phylogenetic groups 4, 5, and 7 of clade II contained solely clusters of PLATZ homologous proteins from angiosperms, suggesting that the corresponding proteins from these groups might have arisen before the monocot-dicot split occurred (~ 200 million years ago) [44]. Intriguingly, group 2 of clade I predominantly consisted of the corresponding proteins from solanaceous crops (tomato and potato) together with the homolog from Arabidopsis (AtPLATZ6) and green alga (CrePLATZ3), suggesting that these members from Arabidopsis and solanaceous crops in group 2 may

Discussion
PLATZ genes, identified specifically in the plant kingdom, are novel zinc-dependent DNA-binding transcription factors that constitute a multi-gene family across a wide variety of plant species [14,32,[36][37][38][39][40]. In support of previous findings, we identified 20 PLATZ genes in the tomato genome, revealing their existence as a multiple-gene family in tomato. This finding suggested that PLATZ genes might play pivotal biological roles in several diverse plant species.
Our phylogenetic analysis divided the PLATZ proteins from various plant species into two major clades (Figure 2). Of these clades, only clade I contained the PLATZ homologs from the green alga Chlamydomonas reinhardtii, suggesting that the proteins from clade I, especially groups 2 and 3, may be more primitive compared with those from clade II, and they might have evolved before the divergence of chlorophytes and streptophytes over one billion years ago [43]. The phylogenetic groups 4, 5, and 7 of clade II contained solely clusters of PLATZ homologous proteins from angiosperms, suggesting that the corresponding proteins from these groups might have arisen before the monocot-dicot split occurred (~200 million years ago) [44]. Intriguingly, group 2 of clade I predominantly consisted of the corresponding proteins from solanaceous crops (tomato and potato) together with the homolog from Arabidopsis (AtPLATZ6) and green alga (CrePLATZ3), suggesting that these members from Arabidopsis and solanaceous crops in group 2 may have emerged from the same progenitor before the divergence of algae and streptophytes. Notably, the tomato PLATZ proteins formed a cluster with their orthologs from potato in groups 1, 4, and 7 with 100% bootstrap support, underlying the evolutionary conservation of these PLATZ members in the solanaceous family.
Conserved domain analysis revealed that the majority of tomato PLATZ proteins were composed of only a single PLATZ domain, while an additional domain was identified in a few of them (Figure 3). This finding indicated that the integration of additional domains in SlPLATZ proteins may contribute to the diversification and enlargement of the SlPLATZ gene family. In addition, the similarity between the first conserved region of the PLATZ domain with certain double zinc finger domains, namely LIM (C 2 HC 5 ) and RING (C 3 HC 4 ) [11], and the second region of the PLATZ domain with the GATA finger (C 2 C 2 ) [9] suggested that the conserved regions of the plant-specific PLATZ domain might have evolved and diversified from its ancestral zinc finger motifs through nucleotide substitution and minor insertion/deletion mutations.
Analysis of the conserved motifs and exon-intron composition in SlPLATZ genes revealed the structural similarity among the majority of SlPLATZ genes clustered in the same phylogenetic groups, suggesting their structural conservation in the more closely related SlPLATZ family members (Figures S1 and S2). The similar motif composition of paralogous and orthologous gene pairs from tomato, Arabidopsis, and rice suggested the important role of these conserved motifs in the expansion of the PLATZ gene family and the evolutionary conservation of those motifs during species evolution.
Gene duplication is an important evolutionary mechanism that enhances the adaption of organisms to diverse environments by generating additional genes [45,46]. The prediction of three segmentally duplicated gene pairs and four tandemly duplicated gene pairs in the SlPLATZ gene family suggested that both segmental and tandem gene duplication triggered the expansion of the tomato PLATZ gene family during evolution. Our finding is in line with the previous reports that identified several duplicated gene pairs responsible for the enlargement of the PLATZ gene family in multiple plant species [14,32,[38][39][40]. Group-specific gene duplication events were observed in the tomato PLATZ gene family since the duplicated genes of each pair belonged to the same phylogenetic groups ( Table 2).
The microsynteny analysis revealed more orthologous gene pairs between tomato and Arabidopsis and very few orthologs between tomato and rice ( Figure 4). This result showed that tomato is evolutionarily more closely related to the dicotyledonous model plant Arabidopsis, in comparison with the monocotyledonous model plant rice. Two SlPLATZ genes were collinear with those of both Arabidopsis and rice, suggesting that these orthologous genes in tomato, Arabidopsis, and rice may have existed before the divergence of the ancestral genes during the evolutionary process.
The upstream regions of numerous stress-related genes have been identified to contain a variety of stress-associated cis-regulatory elements [47][48][49]. The prediction of a diverse number of stress-responsive and hormone-related cis-acting elements in the promoter regions of many SlPLATZ genes ( Figure S6 and Table S4) supported the abiotic stress responsiveness of SlPLATZ genes. For instance, the presence of the ABA-responsive element ABRE and the low-temperature responsive element LTR in the promoter region of SlPLATZ1, the LTR element and the drought-responsive element MBS in that of SlPLATZ12, and the LTR element in that of SlPLATZ19 suggested that these cis elements might act as molecular switches in their response to these stress conditions. Our result was corroborated by the presence of stress-responsive cis-elements in the promoters of PLATZ family members in diverse plant species [38][39][40].
To further explore the molecular structure and functions of SlPLATZ proteins, a 3D structure prediction was conducted ( Figure 5). The comparative modeling of SlPLATZ proteins predicted their binding affinities with diverse ligands, including nucleic acids and peptide molecules. Furthermore, gene ontology (GO) annotations for the obtained SlPLATZ models suggested their putative functions in a variety of biological processes, including cellular defense response, transcriptional regulation, and protein metabolic processes (Table S10). This result suggested their crucial biological roles in tomato, including functioning as transcription factors.
Differential expression of genes in various tissues provides valuable clues about their functional diversity and possible roles in developmental processes. The tissue-specific expression and differential expression patterns of tomato PLATZ gene family members in different organs suggested their functional diversity in the development of tomato ( Figure 6). This result is consistent with the tissue-specific expression and varied expression in various organs of PLATZ genes in other plant species [36,39].
The mRNA transcripts of SlPLATZ5, SlPLATZ14, and SlPLATZ20 were undetectable in any organ examined, suggesting that they could be pseudogenes in the tomato genome or might be expressed only at certain critical developmental stages. In addition to providing mechanical strength to the aerial parts of plants, the stem mediates the long-distance movement of water and minerals to facilitate plant growth under both favorable and unfavorable conditions. The predominant expression of SlPLATZ12 in stems hinted at its putative function in the development of stems and the abiotic stress adaptation of tomato. Flowering and floral development, regulated by numerous floral genes and various environmental factors, are prerequisites for fruit set and important for crop yield [60]. The flower-specific expression of SlPLATZ17 indicated its potential function in the development of flower organs in tomato. Besides SlPLATZ17, the transcript levels of SlPLATZ11 and SlPLATZ19 were the most abundant in flowers, revealing their likely function in the development of flowers ( Figure 6).
SlPLATZ10 is a fruit-specific gene that was expressed preferentially in 1 cm fruits but not in any other organs. This finding suggested that it may have a specific function in the initial developmental phase in tomato. SlPLATZ1 showed the highest expression in 1 cm fruit, but its expression gradually declined at later fruit developmental phases, highlighting its putative function in different developmental stages of tomato fruit, especially in the cell division stage of tomato ( Figure 6).
Plant responses and adaptation to a wide range of abiotic stimuli are closely linked to differential gene expression mediated by a complex network of multiple transcription factors and a variety of stress-related genes via an ABA-dependent or ABA-independent pathway [42,61,62]. The functions of PLATZ transcription factors in plant abiotic stress tolerance and their responses to various stress and hormone treatments have been well documented [14,32,40]. Here we showed that tomato SlPLATZ genes (i.e., SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19) were also responsive to different abiotic stress treatments ( Figure 7A-E), suggesting their functional role in abiotic stress responses of tomato. Drought stress triggered significant alterations in the expression levels of SlPLATZ1, SlPLATZ12, and SlPLATZ19. The transcripts of SlPLATZ1 and SlPLATZ19 were markedly induced at most or all stages of drought treatment, respectively, whereas the expression of SlPLATZ12 was considerably suppressed at all time points following exposure to drought ( Figure 7A). Our finding on the responses of SlPLATZ genes to drought stress is well correlated with a previous report in which the overexpression of PLATZ genes from bamboo and Arabidopsis was observed to promote drought tolerance in the transgenic Arabidopsis lines with respect to the wild-type plants [23][24][25], and another report in which the expression of PLATZ-TF7 (Zm00001d051511) was induced in drought-tolerant maize plants under drought treatment [31].
On exposure to salinity stress, SlPLATZ genes (SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19) displayed remarkable changes in their transcription levels ( Figure 7B). The expression of SlPLATZ12 and SlPLATZ19 was up-or down-regulated under salt treatment, whereas that of SlPLATZ11 and SlPLATZ19 decreased (2-to 3-fold over the control) in the early phase and late phase of the saline stress period, respectively. These responses of SlPLATZ genes under salt treatment revealed their probable involvement in the salt stress adaptation of tomato. This result is consistent with a previous report wherein GhPLATZ1overexpressing transgenic Arabidopsis plants showed robust resistance to osmotic and salt stress [28].
In response to low temperature, the expression of SlPLATZ genes was significantly up-or down-regulated. The transcription levels of SlPLATZ11 and SlPLATZ12 were sharply elevated throughout the cold stress period. Conversely, SlPLATZ1 and SlPLATZ19 were significantly downregulated following cold exposure ( Figure 7C). These findings highlighted the fact that tomato PLATZ genes may also play a role in the adaptation of tomato to low-temperature stress. Our finding is supported by recent work reporting that the expression of most PhePLATZ genes is induced in response to cold treatment [32].
High-temperature stress also remarkably altered the expression profiles of SlPLATZ genes ( Figure 7D). SlPLATZ1, SlPLATZ11, and SlPLATZ12 exhibited similar responses to heat stress, with significantly downregulated transcript levels at the early stages of heat treatment, followed by sharply upregulated expression levels at the later stages. SlPLATZ19 was markedly downregulated in multiple phases of the stress period, indicating the putative functions of SlPLATZ genes in the heat stress response of tomato. This result is corroborated by a previous study reporting the upregulation of the PLATZ gene (Ca01g07220) in the heat-tolerant pepper cultivar compared with the sensitive cultivar [30].
The regulation of various stress-associated genes by the stress hormone ABA has been extensively investigated in order to enhance plant adaptation to unfavorable environmental stresses, including drought, salt, cold, and heat stress [63][64][65]. SlPLATZ1 showed remarkable responses to ABA treatment, while the expression of SlPLATZ19 was not significantly altered and that of SlPLATZ11 and SlPLATZ12 was responsive only at one time point following exposure to ABA treatment ( Figure 7E). This finding implies that SlPLATZ1 may have a possible function in the abiotic stress adaptation of tomato through an ABA-dependent pathway, whereas the function of the other genes (SlPLATZ11, SlPLATZ12, and SlPLATZ19) is unlikely to be related to ABA.
To further explore the functions of SlPLATZ proteins, their subcellular localization was empirically observed under a confocal microscope after the transient expression of the SlPLATZ-GFP fusion proteins in the rice protoplast. The nuclear-specific localization of SlPLATZ12 suggested its specific role in transcriptional regulation as a transcription factor. The localization of SlPLATZ11 and SlPLATZ19 in both the nucleus and the cytoplasm spotlighted the possibility that they might shuttle between the cytoplasm and the nucleus to play a regulatory role in various cellular signaling pathways in addition to their role in gene regulation as transcription factors. Our finding is consistent with the previous studies that reported the localization of PLATZ proteins in the nucleus or the nucleus and cytoplasm [26,39].
To gain a better understanding of the putative functions of SlPLATZ genes, the gene co-expression network analysis of SlPLATZs was performed using the RNA-seq data. The co-expressed genes identified in the networks of SlPLATZ17, SlPLATZ18, and SlPLATZ19 were implicated in various biological pathways, including stress adaptation in tomato, underlying their biological importance in tomato (Figures 8 and 9, Tables S11 and S12).

Identification and Sequence Analysis of PLATZ Genes in the Tomato Genome
We downloaded Arabidopsis thaliana PLATZ protein sequences from TAIR https: //www.Arabidopsis.org/ (accessed on 31 March 2022) and the HMM profile of SlPLATZ (PF04640) from the Pfam database http://pfam.xfam.org/ (accessed on 31 March 2022). Then, we performed basic local alignment search tool (BLAST) searches with default parameters of the Sol genomics database using AtPLATZ protein sequences and the HMM profile of SlPLATZ (PF04640) as queries. The resulting 20 non-redundant SlPLATZ protein sequences were validated using the SMART web tool http://smart.emblheidelberg.de/ (accessed on 1 April 2022), the NCBI CDD search https://www.ncbi.nlm.nih.gov/Structure/ bwrpsb/bwrpsb.cgi (accessed on 1 April 2022), and the Pfam database to reconfirm the existence of the PLATZ domain. The protein length, molecular weight, GRAVY values (grand average of hydropathicity index), and isoelectric points of the identified SlPLATZ proteins were investigated using the Expasy server http://cn.expasy.org/tools/protparam.html (accessed on 4 April 2022). The Open Reading Frame Finder tool https://www.ncbi.nlm. nih.gov/orffinder/ (accessed on 4 April 2022) was employed to determine the open reading frames of the SlPLATZ genes. Clustal Omega and ESPript web tools were used for multipleprotein sequence alignment [84,85]. The Multiple EM for Motif Elicitation (MEME) web server http://meme-suite.org/ (accessed on 18 April 2022) was employed to analyze the conserved motifs in the full-length protein sequences from Arabidopsis, tomato, and rice. The analysis conditions were set as follows: a maximum number of 10 motifs and a motif length between 6 and 50 amino acids. The exon-intron structure of SlPLATZ genes was investigated using the Gene Structure Display Server (GSDS) http://gsds.cbi.pku.edu.cn/ (accessed on 7 April 2022). The web tool "WoLF-PSORT" https://wolfpsort.hgc.jp/ (accessed on 6 May 2022) was used to predict the subcellular localization of the identified SlPLATZ proteins.

Phylogenetic Analysis of Tomato PLATZ Proteins
The full-length PLATZ protein sequences from 10 plant species were aligned using ClustalW, which was followed by phylogenetic analysis using the neighbor-joining (NJ) method with 1000 bootstrap replications in MEGA 6.0 [86]. The deduced PLATZ protein sequences of Chalamdomonas renhadi, Physcomitrium patens, Amborella trichopoda, tomato, potato, Arabidopsis, and Sorghum bicolor were retrieved from the Phytozome database http://www.phytozome.net (accessed on 31 March 2022), while those of rice, maize, and Brassica rapa were obtained from the literature [38]. The gene names and accession numbers used in the phylogenetic tree are shown in Table S1.

Chromosomal Locations, Gene Duplication, and Microsynteny Analysis
The chromosomal locations of SlPLATZ genes were investigated using the Sol genomic database and then visualized with TBtools. Gene duplications among SlPLATZ genes were analyzed with the one-step MCScanX program of TBtools software [87] and examined by BLASTP with an E-value < 10 −10 . The synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates of duplicated SlPLATZ gene pairs were investigated using the method of Nei and Gojobori [88]. The mode of selection was determined by the analyzed Ka/Ks ratio [89]. The divergence time (T) of the predicted duplicated gene pairs was estimated using the formula T = Ks/2r MYA (millions of years ago). Ks is the synonymous substitution rate per site, and r is considered the constant for dicotyledonous plants of 1.5 × 10 −8 substitutions per site per year [90]. The microsynteny analysis of PLATZ genes across tomato, Arabidopsis, and rice was conducted using a reciprocal BLAST search approach against the entire genomes of these species. The duplicated gene pairs were visualized using TBtools software [87].

Prediction of Cis-Acting Elements and miRNA Target Sites
The psRNATarget online tool http://plantgrn.noble.org/psRNATarget/analysis (accessed on 16 July 2022) was employed to predict the putative miRNA sequences in tomato PLATZ genes. The putative cis-regulatory elements located within 1500 bp upstream from the transcription start point of SlPLATZ genes were analyzed using the Plant-Care database http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 11 April 2022).

Homology Modeling of SlPLATZ Proteins
The protein sequences of SlPLATZ1-20 were used as input to perform the 3D structure prediction through the I-TASSER server. The 3D models of tomato PLATZ proteins were generated from multiple threading alignments with LOMETS and iterative TASSER assembly simulations [91]. The template analogs were identified, and the best-modeled structures were selected based on the maximum scores. The generated 3D models were refined using ModRefiner [92]. The final modeled 3D structures and ligand-binding sites of SlPLATZ proteins were visualized with Discovery Studio v.21.1. The putative functions of modeled PLATZ proteins were predicted via the I-TASSER server according to global and local similarity to template proteins curated in the PDB with known structures and functions.

Plant Sample Collection and Stress Treatments
Tomato seeds (Solanum lycopersicum L. cv. Ailsa Craig) were germinated in commercial horticultural soil for seedlings and maintained in a growth chamber controlled at 25 • C day/20 • C night, a 16-h light/8-h dark photoperiod, a relative humidity of 55-70%, and a light intensity of 300 µmol m −2 s −1 until they developed into 28-day-old seedlings. Following the collection of fresh roots, stems, and leaves from 28-day-old plants for organspecific expression analysis, the remaining plants were moved to a greenhouse adjusted at 25/20 • C day/night temperatures until they reached the reproductive stage for the harvesting of flower and fruit samples. For the expression profiling at different fruit developmental stages, the following samples were collected: (i) young fruits approximately 7 days after pollination and 1 cm in diameter (1 cm fruits), (ii) immature fruits approximately 21 days after pollination (IM fruits), (iii) mature green fruits approximately 35 days after pollination (MG fruits), (iv) fruits at the breaker stage when the green color of mature fruits alters to light yellow-orange (B fruits), and (v) fruits 5 days after the breaker stage (B5 fruits) [49].
Leaf samples from 28-day-old tomato seedlings treated with various abiotic stress stimuli (salt (NaCl), drought, heat, cold, and abscisic acid (ABA)) were collected to analyze the expression patterns of SlPLATZ genes in response to different abiotic stress conditions. The transcript levels of SlPLATZ genes were monitored at 0, 1, 3, 9, and 24 h after treatment with heat, cold, salt (NaCl), and abscisic acid (ABA), and at 0, 24, 48, 60, and 72 h after drought treatment. Drought stress was imposed on the tomato plants by withholding water for 72 h [93]. ABA treatment was applied to the plants by spraying leaves with 100 µM ABA [94]. To impose heat and cold stress, the tomato seedlings were incubated in a growth chamber adjusted at 40 • C and 4 • C, respectively. Salt stress was induced by immersing roots in a nutrient solution with 200 mM NaCl, and the plants in the nutrient solution without salt were used as the 0 h control. Tomato seedlings grown in soil under normal conditions (25 • C) served as the 0 h controls for heat, cold, drought, and ABA treatment [95]. All plant samples were collected from three biological replicates, immediately frozen in liquid nitrogen, and stored at −80 • C for further use.

Expression Profiling of SlPLATZ Genes by qRT−PCR
Total RNA was extracted from all samples using an RNeasy Mini kit (Qiagen, Hilden, Germany) and purified with an RNase-free DNase I kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of extracted RNA were measured using a NanoDrop ® 1000 spectrophotometer (Wilmington, DE, USA). A Superscript ® III First-Strand cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA) was used to synthesize cDNA from 1µg of total RNA in accordance with the manufacturer's protocols. The gene-specific primers for SlPLATZ genes were designed using Primer3 software http://frodo.wi.mit.edu/primer3/input.htm (accessed on 25 April 2022) (Table S2), and a melting curve analysis was conducted to validate the specificity of each primer set [95]. The 18S rRNA (F: AAAAGGTCGACGCGGGCT, R: CGACAGAAGGGACGAGAC) from tomato was used as a control gene for normalization [96]. qRT-PCR analysis was conducted in a total volume 10 µL reaction mixture comprising 1 µL (50 ng) cDNA, 2 µL forward and reverse primers of 5 pmol concentration, 5 µL of iTaq SYBR Green (Qiagen, Hilden, Germany), and 2 µL double distilled water. A LightCycler ® 96SW 1.1 (Roche, Germany) was used for amplification and determination of the Cq value of the samples with the following PCR parameters: pre-denaturation at 94 • C for 5 min, followed by 40 cycles at 94 • C for 15 s, annealing at 60 • C for 20 s, and extension at 72 • C for 30 s. The 2 −∆∆Ct method was employed to analyze the relative expression of each gene against each treatment [97].

Co-Expression Network Analysis of SlPLATZ Genes
Raw RNA-Seq data used for the gene co-expression network analysis were downloaded from the NCBI SRA database (Table S3). The quality assessment on raw sequence reads was performed using the FastQC toolkit [98]. The raw reads were cleaned by removing low-quality reads and adapter contamination. HISAT2 2.1.0 was used to map the cleaned reads against the tomato reference genome ITAG4.0 [99]. The number of reads that mapped on exons was analyzed with FeatureCounts 1.6.2 [100]. The expression was analyzed as the fragment per kb per million reads (FPKM) by using DESeq2 software [101]; removeBatchEffect processing was carried out to remove the batch effect in R [102]. For the construction of the co-expression network of SlPLATZ genes, weighted gene co-expression network analysis (WGCNA) was conducted using those genes with an FPKM greater than 1 in R [103]. The co-expression networks were constructed using Cytoscape https://cytoscape.org/ (accessed on 29 June 2022). We used the top co-expressed genes from the network figures for GO and KEGG analysis, which was accomplished using the KOBAS web tool http://kobas.cbi.pku.edu.cn/ (accessed on 11 July 2022).

Subcellular Localization
The coding regions of SlPLATZ genes were amplified with gene-specific primers (Table S13) and then cloned in the pGA3452 vector driven by the maize Ubi1 promoter to generate SlPATZ -GFP fusion proteins [104]. A vector that expresses an NLS-mRFP fusion protein served as a nuclear marker. The SlPLATZ-GFP fusion construct and the NLS-mRFP construct were introduced into rice Oc cell protoplasts using the electroporation method [105]. After incubation overnight at 28 • C in the dark for 12 to 16 h, the resulting transformants were examined for the detection of fluorescent signals under a confocal fluorescence microscope (BX61; Olympus, Tokyo, Japan) using bright field illumination, the GFP channel, and the RFP channel.

Conclusions
We identified 20 non-redundant SlPLATZ genes, most of which harbored two conserved regions essential for zinc-dependent DNA binding, in the tomato genome. Gene duplication analysis spotlighted the fact that both segmental and tandem gene duplication events are accountable for the expansion of the PLATZ gene family in tomato over the course of evolution. The potential functions of SlPLATZ genes in the abiotic stress acclimation of tomato were highlighted by expression profiling that revealed the responses of SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19 to diverse abiotic stresses and gene co-expression network analysis that uncovered the co-expression of several stress-related genes in the networks of SlPLATZ11, SlPLATZ17, and SlPLATZ19. The tissue-specific expression of SlPLATZ10 and SlPLATZ17, and the varied expression levels of SlPLATZ1, SlPLATZ11, SlPLATZ12, and SlPLATZ19 in diverse organs suggested their different functional roles in the development of tomato. The predominant localization of selected SlPLATZ-GFP fusion proteins in the nucleus revealed their putative function as transcription factors. Our results provide valuable information for the further functional elucidation of potential PLATZ transcription factors implicated in the development and abiotic stress tolerance of tomato.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants11223112/s1, The following is attached as supplementary material related to this article: Figure S1. Schematic depiction of the exon-intron distribution of the PLATZ genes in tomato. Figure S2. Schematic representation of 10 conserved motifs identified in the PLATZ proteins from tomato, Arabidopsis, and rice. Figure S3. Overview of conserved motifs of PLATZ proteins from tomato, Arabidopsis, and rice investigated by the MEME web tool. Figure S4. Chromosomal distribution of the SlPLATZ gene family. Figure S5. Schematic diagram of the segmentally duplicated SlPLATZ gene pairs identified in tomato. Figure S6. Prediction of hormone and stress related cis-regulatory elements in the promoter regions of SlPLATZ genes. Figure S7. Homology in the cDNA sequences of SlPLATZ genes. Table S1. List of cis-regulatory elements predicted in the promoters of tomato PLATZ genes. Table S2. miRNA sequence characteristics of PLATZ genes in tomato. Table S3. Parameters used in 3D modeling of SlPLATZ proteins. Table S4. Model validation data for the predicted structures of SlPLATZ proteins. Table S5. Elements of secondary structure predicted in SlPLATZ proteins. Table S6. Predicted ligands and ligand-binding residues in SlPLATZ proteins. Table S7. Gene ontology (GO) annotation for SlPLATZ proteins. Table S8. Significantly enriched KEGG pathways for SlPLATZ co-expressed genes. Table S9. Annotated KEGG pathways and gene ontology of SlPLATZ co-expressed genes. Table S10. List of the PLATZ protein sequences used to construct the phylogenetic tree. Table S11. The primer sequences used in qRT−PCR analysis of the SlPLATZ genes. Table S12. Information of samples used to analyze RNA-seq data. Table S13. The primer sequences used in subcellular localization analysis.