Genome-Wide Analysis of Calmodulin Binding Transcription Activator (CAMTA) Gene Family in Peach (Prunus persica L. Batsch) and Ectopic Expression of PpCAMTA1 in Arabidopsis camta2,3 Mutant Restore Plant Development

Calmodulin-binding transcription activator (CAMTA) is a transcription factor family containing calmodulin (CaM) binding sites and is involved in plant development. Although CAMTAs in Arabidopsis have been extensively investigated, the functions of CAMTAs remain largely unclear in peaches. In this study, we identified five peach CAMTAs which contained conserved CG-1 box, ANK repeats, CaM binding domain (CaMBD) and IQ motifs. Overexpression in tobacco showed that PpCAMTA1/2/3 were located in the nucleus, while PpCAMTA4 and PpCAMTA5 were located in the plasma membrane. Increased expression levels were observed for PpCAMTA1 and PpCAMTA3 during peach fruit ripening. Expression of PpCAMTA1 was induced by cold treatment and was inhibited by ultraviolet B irradiation (UV-B). Driven by AtCAMTA3 promoter, PpCAMTA1/2/3 were overexpressed in Arabidopsis mutant. Here, we characterized peach PpCAMTA1, representing an ortholog of AtCAMTA3. PpCAMTA1 expression in Arabidopsis complements the developmental deficiencies of the camta2,3 mutant, and restored the plant size to the wild type level. Moreover, overexpressing PpCAMTA1 in camta2,3 mutant inhibited salicylic acid (SA) biosynthesis and expression of SA-related genes, resulting in a susceptibility phenotype to Pst DC3000. Taken together, our results provide new insights for CAMTAs in peach fruit and indicate that PpCAMTA1 is associated with response to stresses during development.


Introduction
Calcium ion (Ca 2+ ), as a best characterized second messenger, is involved in diverse physiological processes in plants [1]. For fruit, the complex relationship between calcium, water, cell walls and signaling pathways make calcium an essential player in fruit physiology and development [2]. Calmodulin (CaM) is the important Ca 2+ receptor [3], the expression of CaMs in papaya fruit showed responses to temperature stresses and during fruit development [4]. CaM-binding transcription activator (CAMTA), also named as signal responsive (SR), receives the unambiguous alerts of calcium influxes by CaMs, which is considered to be the largest and best characterized CaM-binding TF family for Ca 2+ signaling in plants [5]. Tobacco ethylene response gene (NtER1) is the earliest CAMTA gene identified in plants [6]. Later, at least 465 CAMTA genes were identified from 112 plants [7]. Plant CAMTA proteins have several conserved functional domains, including a DNA-binding domain (CG-1), a transcription-associated immunoglobulin domain peach CAMTAs under the control of the native AtCAMTA3 promoter in camta2,3 mutant. Our result showed that overexpression of PpCAMTA1 complements the developmental defects of camta2,3 mutant.

Genome-Wide Identification and Protein Properties Analysis of CAMTAs in Peach
Five CAMTA members were identified in the peach genome database using BLASTP program with six Arabidopsis CAMTA protein sequences as queries. Same number of CAMTAs was also identified in various genome assemblies of different peach varieties (Table S1). In this study, the CDS sequences of these five genes were cloned and sequenced, then used for subsequent analysis. The information for the peach CAMTAs gene ID, protein length, molecular weight (MW), isoelectric point (pI) and grand average of hydropathicity (GRAVY values) were listed in Table 1. The peach CAMTA proteins varied in length from 914 to 1131 amino acids (aa), which resulted in variations in MW. The pI values of five CAMTAs range from 5.58 to 6.69, and the instability index (II) varies between 39.54 and 44.84. These proteins were considered as hydrophilic with the negative GRAVY values (Table 1). Next, we extend the information of CAMTA genes to other Rosaceae species. The results showed that fruit crops from Prunus species (sweet cherry, almond and apricot) have similar CAMTAs number with peaches, which were 5, 5 and 4, respectively (Tables 2 and S2). The number of CAMTA genes of strawberry from Rosaceae family is 16, which is more than three times than that of peach ( Table 2 and Table S2). Here, we identified a total of 53 CAMTAs from Rosaceae family fruit crops, including strawberry, pear, apple, sweet berry, almond, apricot and peach.

Phylogenetic, Conserved Domain and Motifs Analysis of CAMTA Genes in Peach
The phylogenetic tree was constructed to investigate the evolutionary relationship among Arabidopsis CAMTAs and peach CAMTAs ( Figure 1A). Based on the phylogenetic tree, CAMTAs were divided into three groups. PpCAMTA1, PpCAMTA2 and PpCAMTA3 were distributed in a same clade with AtCAMTA1 to AtCAMTA3. PpCAMTA4 showed high sequence similarity with AtCAMTA4, while PpCAMTA5 was grouped with At-CAMTA5 and AtCAMTA6. Furthermore, the conserved domains of the peach CAMTA proteins were visualized ( Figure 1B). Analysis of protein structures showed that they contained the essential conserved domains, including CG-1 box, ANK repeats domain, CaM binding domain (CaMBD) and IQ motifs. PpCAMTA2 and PpCAMTA5 were found to contain no TIG domain and were grouped into non-TIG type. The CaMBD amino acids of peach CAMTAs have high sequence similarity with their homologs in Arabidopsis ( Figure 1C). The greatest similarity of the sequences was found between PpCAMTA4 and AtCAMTA4 with only two amino acids difference. The sequence logo for CaMBD in these CAMTAs was shown below the alignment ( Figure 1C).
between PpCAMTA4 and AtCAMTA4 with only two amino acids difference. quence logo for CaMBD in these CAMTAs was shown below the alignment (Figu To further examine the structural features of peach CAMTAs, the conserved were analyzed using the online server, MEME Suite. A total of 13 motifs were fo peach CAMTAs ( Figure S1). Among them, 11 motifs existed in all peach CAMT some motifs corresponded to protein conserved domains. For example, motifs 1, 5 were linked to the conserved CG-1 domain, motif 2 and 13 partly overlapped w domain and CaMBD domain. Conserved domains and motifs indicated the conse of these identified peach CAMTA family members.  To further examine the structural features of peach CAMTAs, the conserved motifs were analyzed using the online server, MEME Suite. A total of 13 motifs were found in peach CAMTAs ( Figure S1). Among them, 11 motifs existed in all peach CAMTAs, and some motifs corresponded to protein conserved domains. For example, motifs 1, 5 and 10 were linked to the conserved CG-1 domain, motif 2 and 13 partly overlapped with IQ domain and CaMBD domain. Conserved domains and motifs indicated the conservatism of these identified peach CAMTA family members.

Exon-Intron Structure Analysis and Chromosomal Locations of Peach CAMTAs
To get gene structure information of peach CAMTAs, the intron-exon structures were analyzed. As shown in Figure 2A, PpCAMTA2, PpCAMTA3 and PpCAMTA5 had 12 introns, while PpCAMTA1 and PpCAMTA4 contained 11 introns. Different numbers of introns were consistent with that of most plant species (10 to 13 introns) [7]. For the intron phases of five peach CAMTA genes, an intron phase pattern of '1102-0110-0200 was observed in PpCAMTA2 and PpCAMTA3, fitting with one of the most common distribution modes across plant CAMTA genes. However, PpCAMTA1, PpCAMTA4 and PpCAMTA5 showed different patterns with '1102-010-0200', '1102-001-0200' and '1102-0010-0200', respectively. Similar to previous studies [7], our results showed that the lengths of the exons corresponding to intron phase pattern '0200' were significantly longer than those corresponding to intron phase pattern '1102'. Chromosomal locations of peach CAMTA genes revealed that the genes were unevenly distributed on three chromosomes ( Figure 2B). Most CAMTA members (PpCAMTA1, 4 and 5) were located on chromosome 1, PpCAMTA3 on chromosome 6, PpCAMTA2 on chromosome 8.
To get gene structure information of peach CAMTAs, the intron-exon structures were analyzed. As shown in Figure 2A, PpCAMTA2, PpCAMTA3 and PpCAMTA5 had 12 introns, while PpCAMTA1 and PpCAMTA4 contained 11 introns. Different numbers of introns were consistent with that of most plant species (10 to 13 introns) [7]. For the intron phases of five peach CAMTA genes, an intron phase pattern of '1102-0110-0200′ was observed in PpCAMTA2 and PpCAMTA3, fitting with one of the most common distribution modes across plant CAMTA genes. However, PpCAMTA1, PpCAMTA4 and PpCAMTA5 showed different patterns with '1102-010-0200', '1102-001-0200' and '1102-0010-0200', respectively. Similar to previous studies [7], our results showed that the lengths of the exons corresponding to intron phase pattern '0200' were significantly longer than those corresponding to intron phase pattern '1102'. Chromosomal locations of peach CAMTA genes revealed that the genes were unevenly distributed on three chromosomes ( Figure 2B). Most CAMTA members (PpCAMTA1, 4 and 5) were located on chromosome 1, PpCAMTA3 on chromosome 6, PpCAMTA2 on chromosome 8.

Promoter Cis-Acting Regulatory Elements of Peach CAMTA Genes
Cis-acting regulatory elements are specific motifs existing in the gene promoter regions, serve as binding sites for TFs to regulate gene transcription [34]. To gain more information about the response of peach CAMTA genes to stress and hormonal signals, cisacting regulatory elements within 2000 bp upstream promoter region of peach CAMTA

Promoter Cis-Acting Regulatory Elements of Peach CAMTA Genes
Cis-acting regulatory elements are specific motifs existing in the gene promoter regions, serve as binding sites for TFs to regulate gene transcription [34]. To gain more information about the response of peach CAMTA genes to stress and hormonal signals, cis-acting regulatory elements within 2000 bp upstream promoter region of peach CAMTA genes were predicted using the PlantCARE database. According to the annotations, we divided the predicted regulatory motifs into three categories: "growth and development", "phytohormone responsive" and "abiotic and biotic responsive". The promoters of peach CAMTA genes contained abundant stress-related and hormone-related regulatory motifs ( Figure 3A). For instance, PpCAMTA1 had TC-rich repeat and WUN-motif which respond to defense and stress. genes were predicted using the PlantCARE database. According to the annotations, we divided the predicted regulatory motifs into three categories: "growth and development", "phytohormone responsive" and "abiotic and biotic responsive". The promoters of peach CAMTA genes contained abundant stress-related and hormone-related regulatory motifs ( Figure 3A). For instance, PpCAMTA1 had TC-rich repeat and WUN-motif which respond to defense and stress.  The promoter regions of all peach CAMTA genes contained a large number of 'light responsive' motifs with 59 totally, including G-box, TCT motif, Box4, etc. Except for AtCAMTA3, the other four gene promoters had more than 50% regulatory motifs responding to light. Besides, the motifs involved in phytohormone responses were found in these promoters, including ABRE (abscisic acid-responsive elements), CGTCA-motif and TGACG-motif (methyl jasmonate-responsive elements), TCA-motif (salicylic acidresponsive elements) and P-box (gibberellin-responsive element) ( Figure 3A). Especially for the PpCAMTA3 promoter, the maximum number of 'methyl jasmonate (MeJA) responsive' motifs was found with 18. In addition, regulatory motifs that were critical for plant growth and development response were discovered in the promoters of peach CAMTA genes. Among these five peach CAMTAs, the promoter of PpCAMTA3 has the largest number of regulatory motifs with 43, followed by PpCAMTA1 and PpCAMTA2 with 32 and 24, respectively ( Figure 3B).

Expression of CAMTAs in Response to Cold, UV-B and MeJA Treatment in Peach Fruit
To explore the functional relevance of peach CAMTA genes, we next examined the expression patterns of CAMTAs under various treatments based on RNA-seq data. Pp-CAMTA1 was significantly induced by low temperature ( Figure 4A), the transcript level started increasing after 1-day cold treatment and kept rising until 28 days, about 2-fold higher than that of the initial point. On the contrary, the transcript levels of PpCAMTA2 and PpCAMTA5 were inhibited by low temperature (inhibited by 29 % and 67 % after 28 days, respectively). The expression of PpCAMTA3 and PpCAMTA4 were induced by low temperature and peaked at 1 d followed by a decline during cold storage. moter region of peach CAMTA genes.
The promoter regions of all peach CAMTA genes contained a large number of 'light responsive' motifs with 59 totally, including G-box, TCT motif, Box4, etc. Except for AtCAMTA3, the other four gene promoters had more than 50% regulatory motifs responding to light. Besides, the motifs involved in phytohormone responses were found in these promoters, including ABRE (abscisic acid-responsive elements), CGTCA-motif and TGACG-motif (methyl jasmonate-responsive elements), TCA-motif (salicylic acid-responsive elements) and P-box (gibberellin-responsive element) ( Figure 3A). Especially for the PpCAMTA3 promoter, the maximum number of 'methyl jasmonate (MeJA) responsive' motifs was found with 18. In addition, regulatory motifs that were critical for plant growth and development response were discovered in the promoters of peach CAMTA genes. Among these five peach CAMTAs, the promoter of PpCAMTA3 has the largest number of regulatory motifs with 43, followed by PpCAMTA1 and PpCAMTA2 with 32 and 24, respectively ( Figure 3B).

Expression of CAMTAs in Response to Cold, UV-B and MeJA Treatment in Peach Fruit
To explore the functional relevance of peach CAMTA genes, we next examined the expression patterns of CAMTAs under various treatments based on RNA-seq data. PpCAMTA1 was significantly induced by low temperature (Figure 4A), the transcript level started increasing after 1-day cold treatment and kept rising until 28 days, about 2-fold higher than that of the initial point. On the contrary, the transcript levels of PpCAMTA2 and PpCAMTA5 were inhibited by low temperature (inhibited by 29 % and 67 % after 28 days, respectively). The expression of PpCAMTA3 and PpCAMTA4 were induced by low temperature and peaked at 1 d followed by a decline during cold storage. Transcript levels of PpCAMTA1 in UV-B treated group were reduced by approximately 50% relative to those of the control ( Figure 4B). PpCAMTA2 and PpCAMTA5 were suppressed slightly by UV-B. In contrast, transcript level of PpCAMTA4 in UV-B treated peach fruit was approximately 2-fold higher than that of control at 6 h, and decreased to the control level at 48 h. With regard to MeJA treatment, PpCAMTA3 expression was inhibited by 20% at 1 d ( Figure 4C). This was consistent with the rich MeJA-responsive motifs in PpCAMTA3 promoter (Figure 3). PpCAMTA4 and PpCAMTA5 were slightly induced after being treated by MeJA for 3 days. Transcript levels of PpCAMTA1 and PpCAMTA2 remained unchanged after MeJA treatment ( Figure 4C).
To confirm expression patterns of peach CAMTAs to various treatments, RT-qPCR analysis was performed. In agreement with RNA-seq results ( Figure 4B), induced transcript levels of PpCAMTA4 by UV-B were observed at 6 h as well ( Figure S2A). UV-B treatment significantly reduced transcript levels of PpCAMTA1 at both 6 h and 48 h ( Figures 4B and S2A). Regarding to MeJA treatment, transcript level of PpCAMTA3 was reduced significantly at 1 d ( Figure S2B), in agreement with RNA-seq results ( Figure 4C). Significant induction of PpCAMTA4 expression was detected at 3 d after MeJA treatment for both RNA-seq and RT-qPCR ( Figures 4C and S2B). Overall, similar gene expression patterns were observed between RT-qPCR and RNA-seq. Our results indicated that peach fruit CAMTAs responded to various stimuli in a different manner.

Expression of CAMTAs in Different Organs and During Fruit Development
To investigate the spatial and temporal expression patterns of peach CAMTA genes, their transcript levels in leaf, flower and fruit at different developmental stages were analyzed using RNA-seq. All peach CAMTA genes showed the highest transcript levels in ripening peach fruit comparing to leaf and flower ( Figure 5A). Notably, the transcript level of PpCAMTA4 in ripe fruit was 3.3 and 3.5-fold of that in leaf and flower, respectively. The expression of PpCAMTA5 in flowers remained at the similar level as in ripe fruit. At the fruit development stages, PpCAMTA1 was constantly accumulated with fruit ripening, peaked at 6 DAH when transcript level was about 2-fold higher than that at 34 DAB ( Figure 5B). A similar pattern was observed for PpCAMTA3. The transcript level of PpCAMTA2 increased slightly and reached the maximum at 71 DAB. The sustained decline was found in the transcript level of PpCAMTA5, which was inhibited by approximately 60% at 6 DAH compared with that at 34 DAB. Transcript levels of PpCAMTA4 did not appear to be affected by fruit development. Taken together, peach CAMTAs exhibited fruit-specific and ripening-associated expression patterns, PpCAMTA1 and PpCAMTA3 were enriched in ripe fruit, suggesting that they may play important roles in peach fruit ripening.

Subcellular Localization of Peach CAMTA Proteins
As CAMTA family members are transcription factors that contain the nuclear localization signal (NLS) domain, they were considered to be located in cell nucleus [23]. The NLS domains were predicted in peach CAMTA proteins ( Figure 1B). In order to further understand the role of peach CAMTAs, subcellular localizations of these proteins were carried out by transiently expressing the PpCAMTAs-GFP constructs in tobacco leaves. For PpCAMTA1, PpCAMTA2 and PpCAMTA3, fluorescence microscopy showed distinct bright spots detected on the cell nucleus ( Figure 6A), suggesting their nuclear localization. However, the green fluorescence from PpCAMTA4-GFP and PpCAMTA5-GFP proteins co-localized with the red fluorescence from mCherry-H + -ATPase in the plasma membrane of tobacco leaves ( Figure 6B). This extranuclear localization of CAMTA proteins has also been reported in Arabidopsis [19]. Taken together, our results showed that PpCAMTA1, PpCAMTA2 and PpCAMTA3 were in the nucleus, PpCAMTA4 and PpCAMTA5 were located in the plasma membrane of the cell.

Subcellular Localization of Peach CAMTA Proteins
As CAMTA family members are transcription factors that contain the n zation signal (NLS) domain, they were considered to be located in cell nucl NLS domains were predicted in peach CAMTA proteins ( Figure 1B). In ord understand the role of peach CAMTAs, subcellular localizations of these p carried out by transiently expressing the PpCAMTAs-GFP constructs in to For PpCAMTA1, PpCAMTA2 and PpCAMTA3, fluorescence microscopy sho bright spots detected on the cell nucleus ( Figure 6A), suggesting their nuclear However, the green fluorescence from PpCAMTA4-GFP and PpCAMTA5co-localized with the red fluorescence from mCherry-H + -ATPase in the plasm of tobacco leaves ( Figure 6B). This extranuclear localization of CAMTA pro been reported in Arabidopsis [19]. Taken together, our results showed that PpCAMTA2 and PpCAMTA3 were in the nucleus, PpCAMTA4 and PpCA located in the plasma membrane of the cell.

Overexpressing PpCAMTA1 Restored the Plant Size of Camta2,3 Mutant
The evolutionary tree has demonstrated the close phylogenetic relationships between peach CAMTA1/2/3 and Arabidopsis CAMTA1/2/3 ( Figure 1A). In order to answer the question of whether PpCAMTA1/2/3 represents the functional homologs of AtCAMTA1/2/3, we expressed PpCAMTA1, PpCAMTA2 and PpCAMTA3 in the camta2,3 mutant, respectively ( Figure S3). We chose camta2,3 as the genetic background for the same reasons as the previous study [15], camta2,3 double mutant shows greater changes in SA biosynthesis and expression of SA pathway genes than those in the camta3 single mutant, while camta1,2,3 triple mutant is tiny in size and difficult to work with.
To avoid the excessively strong expression caused by the CaMV35S promoter, we replaced the 35S promoter of the pBI121 vector with the endogenous promoter of AtCAMTA3 to ensure the natural expression level ( Figure 7A). The results showed that overexpression of these three peach CAMTA genes complemented the plant size of camta2,3 to different extent. As shown in Figure 7B, the plant size of OE-PpCAMTA1 transgenic lines were larger than camta2,3 and closed to that of wild type (WT). Although OE-PpCAMTA2 and OE-PpCAMTA3 transgenic lines did not complement the plant size to the level of wild type, expression of these two genes did result in larger plant size than that of camta2,3 mutant ( Figure 7B). In brief, peach PpCAMTA1 could largely rescue the plant size defect phenotype caused by camta2,3 mutation, PpCAMTA2 and PpCAMTA3 partially restored plant development after overexpressing in camta2,3 mutant.

Overexpressing PpCAMTA1 Restored the Plant Size of Camta2,3 Mutant
The evolutionary tree has demonstrated the close phylogenetic relationships b peach CAMTA1/2/3 and Arabidopsis CAMTA1/2/3 ( Figure 1A). In order to ans question of whether PpCAMTA1/2/3 represents the functional homol AtCAMTA1/2/3, we expressed PpCAMTA1, PpCAMTA2 and PpCAMTA3 in the c mutant, respectively ( Figure S3). We chose camta2,3 as the genetic background same reasons as the previous study [15], camta2,3 double mutant shows greater in SA biosynthesis and expression of SA pathway genes than those in the camta

Overexpressing PpCAMTA1 Reduces SA Biosynthesis and Weakens Plant Resistance to Pathogen of Camta2,3 Mutant
Our results showed that overexpressing PpCAMTA1 recovered camta2,3 mutant size close to that of WT ( Figure 7B), indicating that PpCAMTA1 is a candidate gene with a role similar to Arabidopsis CAMTA3 during plant development. Arabidopsis CAMTA1, CAMTA2 and CAMTA3 repress SA biosynthesis and ICS1 expression by regulating tran-scription factors at ambient temperature [14,17]. In addition, CAMTA3 regulates SARD1 and CBP60g expression [18]. Therefore, we next explored the role of PpCAMTA1 in plant immune response. genic lines were larger than camta2,3 and closed to that of wild type (WT). Althou PpCAMTA2 and OE-PpCAMTA3 transgenic lines did not complement the plant siz level of wild type, expression of these two genes did result in larger plant size th of camta2,3 mutant ( Figure 7B). In brief, peach PpCAMTA1 could largely rescue th size defect phenotype caused by camta2,3 mutation, PpCAMTA2 and PpCAMTA3 p restored plant development after overexpressing in camta2,3 mutant.

Overexpressing PpCAMTA1 Reduces SA Biosynthesis and Weakens Plant Resistance Pathogen of Camta2,3 Mutant
Our results showed that overexpressing PpCAMTA1 recovered camta2,3 muta close to that of WT ( Figure 7B), indicating that PpCAMTA1 is a candidate gene with similar to Arabidopsis CAMTA3 during plant development. Arabidopsis CA Significantly higher SA content and related genes expression levels were observed in camta2,3 mutant than those in WT (Figure 8). After overexpressing PpCAMTA1, transgenic plants resulted in at least 52 % reduction in SA content in relative to that of camta2,3 mutant ( Figure 8A). Correspondingly, transcript levels of AtICS1 (a gene responsible for SA synthesis) were reduced by approximately 90% in OE-PpCAMTA1 lines when compared to that in camta2,3 mutant ( Figure 8B). Besides SA content and SA synthesis related gene expression, expression levels for SA related transcription factors were analyzed as well.
Transcription levels of AtCBP60g in transgenic lines were about 5.3%, 7.3% and 14.2% of that in camta2,3, respectively ( Figure 8C). Moreover, in transgenic lines after overexpressing PpCAMTA1, expression levels of AtSARD1, AtEDS1 and AtPAD4 were significantly lower than those in camta2,3 mutant as well. It is noteworthy that both SA content and expression levels for genes associated with SA synthesis and signaling were similar between WT and transgenic plants after overexpressing PpCAMTA1. tant ( Figure 8A). Correspondingly, transcript levels of AtICS1 (a gene responsible for SA synthesis) were reduced by approximately 90% in OE-PpCAMTA1 lines when compared to that in camta2,3 mutant ( Figure 8B). Besides SA content and SA synthesis related gene expression, expression levels for SA related transcription factors were analyzed as well. Transcription levels of AtCBP60g in transgenic lines were about 5.3%, 7.3% and 14.2% of that in camta2,3, respectively ( Figure 8C). Moreover, in transgenic lines after overexpressing PpCAMTA1, expression levels of AtSARD1, AtEDS1 and AtPAD4 were significantly lower than those in camta2,3 mutant as well. It is noteworthy that both SA content and expression levels for genes associated with SA synthesis and signaling were similar between WT and transgenic plants after overexpressing PpCAMTA1.  Given the above phenotype, Pst DC3000 (Pseudomonas syringae pv. tomato DC3000) was used to measure the effect of PpCAMTA1 on pathogen resistance. Three days after pathogen infection, camta2,3 mutant plants exhibited scarce chlorotic spots on their leaves, while wild type and OE-PpCAMTA1 lines exhibited severity chlorotic spots on leaves, showed more susceptibility to the Pst strain ( Figure 9A). Gene expression was also shown to induce stronger AtPR1 in camta2,3 mutants after 1 d of inoculation when compared with OE-PpCAMTA1 lines ( Figure 9B). The expression of SA-related genes mentioned above showed similar trends among these plants (Table S6). Taken together, these results revealed that PpCAMTA1 negatively regulated SA biosynthesis and genes related to SA pathway, which may act as a negative regulator of plant immunity in peach.
showed more susceptibility to the Pst strain ( Figure 9A). Gene expression was also s to induce stronger AtPR1 in camta2,3 mutants after 1 d of inoculation when compare OE-PpCAMTA1 lines ( Figure 9B). The expression of SA-related genes mentioned showed similar trends among these plants (Table S6). Taken together, these resu vealed that PpCAMTA1 negatively regulated SA biosynthesis and genes related pathway, which may act as a negative regulator of plant immunity in peach.

Discussion
As important transcription factors containing CaM binding sites, the functi CAMTAs have been reported widely in recent years [7, 23,35]. Although compreh genome-wide analysis of CAMTAs have been reported in horticultural crops [26][27][28][29] tions of CAMTAs in fruit remain largely unknown. Peach is considered as a mode of Rosaceae family [31]. In this study, we identified five peach CAMTA genes. Similar ber of CAMTA genes were also observed for sweet cherry, apricot and almond. R family fruit crops such as apple and strawberry have more CAMTA gene member number of CAMTA genes varing in different plant species is suggested to be asso with evolution although the process is complex [7].
Pathogen infection and disease development are major factors that cause posth decay, leading to a limitation in the storage period and marketing life of peach fru Improving immunity resistance is an efficient and preferred strategy for reducing d development in fruit [37]. Several studies have shown that application of calcium and after harvesting delayed the development of peach fruit disease, and thus prol storage period of peach fruit [38][39][40]. Therefore, it is necessary to identify fruit im related transcription factors on calcium signaling pathways. In Arabid

Discussion
As important transcription factors containing CaM binding sites, the functions of CAMTAs have been reported widely in recent years [7, 23,35]. Although comprehensive genome-wide analysis of CAMTAs have been reported in horticultural crops [26][27][28][29], functions of CAMTAs in fruit remain largely unknown. Peach is considered as a model plant of Rosaceae family [31]. In this study, we identified five peach CAMTA genes. Similar number of CAMTA genes were also observed for sweet cherry, apricot and almond. Rosaceae family fruit crops such as apple and strawberry have more CAMTA gene members. The number of CAMTA genes varing in different plant species is suggested to be associated with evolution although the process is complex [7].
Pathogen infection and disease development are major factors that cause postharvest decay, leading to a limitation in the storage period and marketing life of peach fruit [36]. Improving immunity resistance is an efficient and preferred strategy for reducing disease development in fruit [37]. Several studies have shown that application of calcium before and after harvesting delayed the development of peach fruit disease, and thus prolonged storage period of peach fruit [38][39][40]. Therefore, it is necessary to identify fruit immune-related transcription factors on calcium signaling pathways. In Arabidopsis, At-CAMTA1/2/3 are components derived from calcium signaling pathway and play negative roles in plant immune responses. For instance, mutant camta2,3 has increased SA content and higher expression levels of SA-regulated immunity genes than those in wild type Arabidopsis [14]. Overexpressing PpCAMTA1 in the camta2,3 mutants resulted in lower SA content, contributing to more susceptible to the virulent Pst DC3000 strain. Notably, overexpression of PpCAMTA1 restored plant size, close to that of WT plant. These results indicated that PpCAMTA1 could functionally complement Arabidopsis CAMTA2/3 by overexpressing peach gene in cmata2,3 mutant. During peach fruit storage at ambient temperature, transcript levels of PpCAMTA1 increased, accompanied with fruit postharvest ripening and senescence. Similarly, increased expression levels were also observed for tomato SlSR1/1L [10] and durian DzCAMTA3/8 during postharvest fruit ripening [30]. Induced expression of tobacco CAMTA member NtER1 was suggested to be associated with ethylene induced senescence for both leaf and flower [6]. According to expression pattern of PpCAMTA1 and transgenic Arabidopsis phenotype in camta2,3 genetic background, our results suggested that increased PpCAMTA1 expression was associated with reduced peach fruit immune response during postharvest storage, thus contributing to fruit senescence and seed dispersal.
Apart from biological stress, expression of plant CAMTAs could also be affected by abiotic stresses [23]. However, effect of abiotic stresses on expression of CAMTAs has not been well studied in fruit species. In our study, we investigated peach CAMTA genes expression profiles during peach fruit ripening and in response to stresses. The present study revealed that PpCAMTA1 expression was induced by cold treatment in peach fruit, while the expression of the left four members was inhibited during postharvest cold storage for 28 d. In tomato fruit, transcript levels of SlSR2 and SlSR3L were induced by cold treatment [10]. Transcript levels of Arabidopsis AtCAMTA3 increased and stayed at high level during exposure to low temperature for 5 weeks [14]. These results showed that induced expression of CAMTAs was conserved both in plants and in fruit species, indicating an important role of peach PpCAMTA1 during fruit cold storage. UV-B irradiation for Arabidopsis seedlings induced expression of CAMTAs [41]. Here, we showed that expression levels of PpCAMTA1, PpCAMTA2 and PpCAMTA5 were inhibited by UV-B. Based on RT-qPCR results, expression level of PpCAMTA3 was also repressed by UV-B ( Figure S2A). These results revealed different response of CAMTA expression patterns between Arabidopsis leaf and peach fruit upon to UV-B treatment. Our previous study showed that UV-B had great effects on contents of both anthocyanin and volatile linalool in peach fruit [42,43], suggesting that peach CAMTAs is associated with the formation of fruit secondary metabolites.
Most previous reports showed that CAMTAs belong to nuclear protein [9,44]. Arabidopsis CAMTA2, CAMTA3 and CAMTA5 were located in nuclei [19]. The present study showed that peach CAMTAs were located both in nucleus and plasma membrane. For instance, PpCAMTA1, PpCAMTA2 and PpCAMTA3 were located in nucleus, and Pp-CAMTA4 and PpCAMTA5 were located in the plasma membrane. Phylogenetic analysis showed that PpCAMTA4 and PpCAMTA5 were on different clades from PpCAMTA1-3 due to difference in sequence. A previous study analyzed the subcellular localization of 465 CAMTA proteins, of which 390 proteins were predicted to located in the nucleus, and the remaining CAMTAs were localized in the cytosol, chloroplast and plasma membrane [7]. Some CAMTAs contain nuclear localization signal sequence which allows these proteins to enter the nucleus to perform their functions. Moreover, these extranuclear localizations of CAMTAs have also been reported in Arabidopsis. For example, destabilization and nuclear export were observed for AtCAMTA3 after protein phosphorylation modification caused by mitogen-activated protein kinases (MAPKs) MPK3 and MPK6 [45]. These two MAPKs were induced by bacterial flg22 elicitor in Arabidopsis. Interestingly, Arabidopsis AtCAMTA1, AtCAMTA4, and AtCAMTA6 were suggested to be unstable at room temperature or 4 • C [19]. These results indicated that it is not surprising that PpCAMTA4/5 were located outside the nucleus. The mechanisms need to be further investigated for the different location of PpCAMTAs in future.
Among five peach CAMTA gene members, PpCAMTA1, PpCAMTA2 and PpCAMPT3 were grouped with Arabidopsis AtCAMTA3 which is involved in plant development. To test if peach CAMTA genes could functionally complement the defect phenotype in Arabidopsis mutant, PpCAMTA1, PpCAMTA2 and PpCAMPT3 were overexpressed driven by the native AtCAMTA3 promoters. The second reason to choose the native At-CAMTA3 promoter is to avoid the excessively strong expression caused by the CaMV35S promoter. We found that overexpression of these three peach CAMTA genes complemented the plant size of camta2,3 to different extent. Peach PpCAMTA1 rescued the developmental defect phenotype caused by camta2,3 mutation, restored the plant size to the WT level. Notably, only partially restored plant development phenotype was observed for Pp-CAMTA2 and PpCAMTA3. However, we cannot exclude roles of other peach CAMTA genes during plant development. The identification of functions of each peach CAMTA gene will require further research. Our dataset in hand suggest that PpCAMTA1 is associated with plant development, including fruit ripening and postharvest responses to biotic and abiotic stresses.

Plant Material and Treatment
Peach (Prunus persica L. Batsch cv. Hujingmilu) fruit at different developmental stages, fully expanded leaves and flowers at bloom stage were used for gene expression analysis according to our previous study [46]. After harvest at mature stage (108 days after bloom, DAB), peach fruits were held at 20 • C for 3 d (3 days after harvest, DAH) and 6 d (6 DAH). Treatments with UV-B irradiation [42] and methyl jasmonate (MeJA) [47] were used in this study. For cold treatment, peach fruits were stored at 0 • C up to 28 days. Three biological replicates with five fruits each were sampled, frozen in liquid nitrogen immediately and stored at −80 • C for further analysis.

Identification of Peach CAMTA Genes
The peach CAMTA genes were identified by the following approaches. Firstly, searching in the peach genome database at Phytozome v12.1 website (https://phytozome.jgi.doe. gov, accessed on 24 February 2022) based on gene annotations. Secondly, sequences of Arabidopsis CAMTA members were used as queries to search against the peach genome database with BLASTN and BLASTP programs. Thirdly, sequences of peach CAMTAs were checked for CAMTA conserved protein domains by using HMMER (http://www.ebi.ac.uk, accessed on 24 February 2022), Pfam database (http://pfam.janelia.org/, accessed on 24 February 2022) and SMART (http://smart.embl-heidelberg.de/, accessed on 24 February 2022). Peach CAMTA genes were cloned and sequenced for the subsequent analysis. Finally, sequence information of CAMTA genes was also extended to various genome assemblies in Prunus persica and Rosaceae genome whose sequence were gained from GDR website (https://www.rosaceae.org/, accessed on 22 August 2022).

Physicochemical Property Analysis and Phylogenetic Analysis
Amino acid properties and physicochemical traits of the given peach CAMTA proteins were calculated using the ProtParam tool in the ExPASy web server (https://web.expasy. org/protparam/, accessed on 28 February 2022). The full-length amino acid sequences of CAMTAs in Arabidopsis and peach were used to perform the multiple sequence alignments using Clustal W with default parameters. The phylogenetic analysis was constructed using the neighbor-joining (NJ) method with 1000 bootstrap replications by MEGA 11. Peach CAMTA genes were named based on the genetic relationship with Arabidopsis CAMTA genes.

Exon-Intron Organization and Chromosomal Map Construction
The distribution patterns of exon-intron of peach CAMTA genes were obtained from the peach genome database and displayed using the online website Gene Structure Display Server v2.0 (GSDS, http://gsds.cbi.pku.edu.cn/, accessed on 11 March 2022). The locations of CAMTA genes on the chromosomes were visualized by MapChart (v2.3) (Voorrips, R.E., Wageningen, The Netherlands).

Conserved Protein Domains Distribution and Cis-Regulatory Element Analysis
Bioinformatic analysis of the conserved protein domains was conducted through the online SMART tool (http://smart.embl-heidelberg.de/, accessed on 28 February 2022). CaM binding domain (CaMBD) was predicted by the Calmodulin Target Database (http://calcium.uhnres.utoronto.ca/ctdb/ctdb/, accessed on 28 February 2022). A schematic representation of conserved protein domain structures was constructed using DOG 2.0 software (Jian Ren, Hefei, China) [48]. Online server MEME (Multiple Em for Motif Elicitation) was used to carry out motif analysis. For promoter cis-regulatory elements analysis, putative promoter sequences (2 kb from TSS) of the peach CAMTA genes were downloaded from the peach genome database, and then identification of possible cis-acting regulatory elements was performed by PlantCARE (http://bioinformatics.psb.ugent.be/webtools/ plantcare/html/, accessed on 11 March 2022).

Subcellular Localization Analysis
Subcellular localization analysis was conducted according to our previous study [49]. The recombined 35S-PpCAMTA-eGFP vectors were constructed using primers listed in Table S3. To identify the nuclear location, the vector was infiltrated into transgenic tobacco plants expressing a red fluorescent nuclear marker (Nucleus-RFP). To verify the plasma membrane localization, mCherry-H + -ATPase with red fluorescence was used as a marker of plasma membrane [50]. Leaves were detached for analysis using a confocal laser scanning microscope (LSM 780; Carl Zeiss, Oberkochen, Germany).

RNA Extraction and Gene Expression Analysis
Total RNA of peach fruit was extracted according to our previous study [51], and libraries for high-throughput Illumina strand-specific RNA-sequencing were prepared as described previously [52]. Reads per kilobase of exon model per million mapped reads (RPKM) based on the length of the gene and the number of reads mapped to this gene was used to express transcript abundance [42,46]. Three biological replicates for each sampling time point were performed for RNA-Seq.
Total RNA of Arabidopsis plant leaves was extracted by using the reagent RNAiso Plus (Takara) according to its instructions. For quantitative reverse-transcription PCR (RT-qPCR) analysis, HiScript II Q RT SuperMix for qPCR (+gDNA wiper) and ChamQ Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China) were used. Each RT-qPCR analysis contained three biological replicates. Oligonucleotide primers for RT-qPCR analysis were listed in Table S4.

Stable Overexpression of Arabidopsis
The transformation was conducted in four-week-old seedlings of Arabidopsis plants according to the floral dip method. The camta2,3 double mutant T-DNA insertion mutants were in Col-0 background and got from Kim [14]. The modified pBI121-PpCAMTA1/2/3 constructs were used to obtain overexpression (OE) lines in the camta2,3 genetic background. Notably, the overexpression of peach PpCAMTA1/2/3 was driven by the AtCAMTA3 promoter rather than the CaMV35S promoter. Primers were listed in Table S5. Transgenic seeds were selected on 1/2 Murashige and Skoog (MS) medium with 50 mg/L kanamycin. Wild type (WT) and mutants were planted in a growth chamber with 22 • C under 16 h light (100 µmol m −2 s −1 constant light) and 8 h dark. After 1-month growth, rosette-like basal leaves of the above-ground part were sampled for subsequent analysis.

Pathogen Infection Assay
For infection with Pst DC3000 (Pseudomonas syringae pv. tomato DC3000), some modifications were made by referring to the methods of previous studies [18]. Briefly, leaves of 4-week-old plants were sprayed with the bacteria at a dose of OD 600 = 0.002 in 10 mM MgCl 2 . Infected leaves were covered with a clean dome to maintain high humidity, grown at 22 • C under 16/8-h light/dark cycles in a growth chamber. On 1 day after inoculation, the leaf samples were collected for RNA extraction and gene expression analysis, and plants scored 3 days for phenotypic observation.

Quantification of Salicylic Acid (SA) Content
Free SA content was extracted and measured using the previous methods with some modifications [53]. A total of 0.2 g fresh weight of leaf tissue powder was homogenized and extracted for 4 h at 65 • C using 1.2 mL of 90% (v/v) methanol. After centrifugation, the supernatants were dried in a speed vacuum at 60 • C. The residue was resuspended using 600 µL 5% (v/v) trichloroacetic acid and sonicated for 10 min. Then, 2 vol (1.2 mL) of ethylacetate-cyclopentane-isopropanol (50:50:1) were added to extract SA. The superorganic phase containing the SA was dried in a speed vacuum with 60 • C. The residue was dissolved in 500 µL methanol. SA was detected using a fluorescence detector in HPLC analysis, the excitation wavelength is 295-nm, and the emission wavelength is 405-nm.

Statistical Analysis
A completely randomized design was used in the experiments. Mean value and standard errors (SE) were calculated by Microsoft Excel, figures were produced by OriginPro 9 (Microcal Software Inc., Northampton, MA, USA). The two-sample significance test was calculated using unpaired Student's t-test (*, p < 0.05 and **, p < 0.01). For multiple samples significance test, ANOVA analysis followed by a post-hoc Tukey HSD test was used to calculate p value, and letters were designated for significantly different values (SPSS 19.0; SPSS Inc., Chicago, IL, USA).

Conclusions
In this study, a total of five CAMTA genes were identifed in peach (Prunus persica L. Batsch) at the level of gene structure, sequence characteristics, promoter cis-acting elements and expression patterns. Increased transcription levels were detected for PpCAMTA1 during fruit development and ripening. Our hypothesis that PpCAMTA1 is an ortholog of AtCAMTA3 was experimentally confirmed by overexpressing PpCAMTA1 under the control of the native AtCAMTA3 promoter in Arabidopsis mutant system. Overexpressing PpCAMTA1 can complement the developmental deficiencies in Arabidopsis camta2,3, indicating an association of PpCAMTA1 during developmental process such as fruit ripening and postharvest responses to biotic and abiotic stresses.

Data Availability Statement:
The relevant datasets can be found in the NCBI with accession number PRJNA576753, PRJNA574004, PRJNA852345 and SRP103523.