Skip to main content

Genome-wide characterization and expression profiling of E2F/DP gene family members in response to abiotic stress in tomato (Solanum lycopersicum L.)

Abstract

Background

E2F/DP (Eukaryotic 2 transcription factor/dimerization partner) family proteins play an essential function in the cell cycle development of higher organisms. E2F/DP family genes have been reported only in a few plant species. However, comprehensive genome-wide characterization analysis of the E2F/DP gene family of Solanum lycopersicum has not been reported so far.

Results

This study identified eight nonredundant SlE2F/DP genes that were classified into seven groups in the phylogenetic analysis. All eight genes had a single E2F-TDP domain and few genes had additional domains. Two segmental duplication gene pairs were observed within tomato, in addition to cis-regulatory elements, miRNA target sites and phosphorylation sites which play an important role in plant development and stress response in tomato. To explore the three-dimensional (3D) models and gene ontology (GO) annotations of SlE2F/DP proteins, we pointed to their putative transporter activity and their interaction with several putative ligands. The localization of SlE2F/DP-GFP fused proteins in the nucleus and endoplasmic reticulum suggested that they may act in other biological functions. Expression studies revealed the differential expression pattern of most of the SlE2F/DP genes in various organs. Moreover, the expression of E2F/DP genes against abiotic stress, particularly SlE2F/DP2 and/or SlE2F/DP7, was upregulated in response to heat, salt, cold and ABA treatment. Furthermore, the co-expression analysis of SlE2F/DP genes with multiple metabolic pathways was co-expressed with defence genes, transcription factors and so on, suggested their crucial role in various biological processes.

Conclusions

Overall, our findings provide a way to understand the structure and function of SlE2F/DP genes; it might be helpful to improve fruit development and tolerance against abiotic stress through marker-assisted selection or transgenic approaches.

Peer Review reports

Background

Due to various biotic and abiotic stresses, plants have evolved complex signalling pathways for their defence mechanism [1]. The innate mechanism in plants induces the release of secondary metabolism and phytohormones, signalling receptors and signal transducers to active defence pathways [2]. Some signalling pathways are reported to be common in plants against both abiotic and biotic stress [3]. Different abiotic stresses, such as heat, cold, drought, and salt, trigger plant sensors by a single or simultaneous attack to activate the resistance mechanism [4, 5]. Many of the transcription factor (TFs) regulate gene expression by binding to cis-acting elements in the promoter region against abiotic stress [6, 7]. Transcription factors (TF) induce the genes involved in cell signalling, plant hormone metabolism, pathogenesis-related protein and in primary metabolism and secondary metabolism of defence (SAR). Till now, around 320k TFs have been reported from 165 different plant species. The most-known TF families are WRKY, MBS, MADS, ARE, bZIP, GRAS etc [7,8,9]. In addition, some sub-category TF families, SNF7, CPP, and E2F/DP are also reported to be involved in defence against heat stress in Arabidopsis and rice [10].

In plants and animals, the Eukaryotic 2 transcription factor/dimerization partner transcription factor family (henceforth E2F/DP) is reportedly involved in the cell cycle and its development. The E2F/DP transcription factors in plants are categorized into E2F, DP, and DEL (DP-E2F-like) groups based on their conserved domains. The E2F gene group contains four functional domains: an RBR-binding domain, a DNA-binding domain, a ‘marked box’ domain and a leucine zipper dimerization domain [11]. However, the members of the DP group contain a DNA-binding domain and/or leucine zipper dimerization domain when compared with the E2F genes. The DEL group of genes has only a DNA-binding domain [11].

The E2F/DP TFs are reportedly involved in cell growth and proliferation in mammalian E2F signalling pathways [12]. Likewise, in Arabidopsis plants the E2F signalling pathway functions as a core regulator with cyclin-dependent kinase inhibitors (CKIs), retinoblastoma (RBs), cyclins, and cyclin-dependent kinases (CDKs). In plants and mammals, the E2F transcription factor binds to the TTTCCCGCC motif with or without slight modification in its promoter region of the target genes [13]. Another study of E2F family gene structural analysis revealed that most of the genes shared similar intron numbers/lengths within the same groups.

In Arabidopsis, the E2Fa and E2Fb genes bind with DP to regulate their gene expression [14]. Similarly, in Moso bamboo, the PheE2F/DP genes were reportedly involved in leaf and root growth, genome integrity and viability, and abiotic stress response [14,15,16,17,18]. Most of the PheE2F/DP genes were upregulated against drought and salt stress at different time points. Similarly, four genes (E2Fa, E2Fb, E2Fc and E2FDPa) showed more than 2-fold against salt stress in Medicago truncatula [18, 19]. In wheat (Triticum aestivum), TaE2F1I-19 and TaDP3III-15 were upregulated against drought stress; TaDP2III-3 and TaDEL2II-27 genes were expressed against salt stress and heat stress; TaE2F1I-19 and TaDP3III-15 were significantly upregulated against cold stress [20].

Tomato (Solanum lycopersicum) is a model vegetable crop cultivated worldwide, and it stands in the fourth position among the leading vegetables in terms of production [21]. Tomato is grown on an area of 5.06 million hectares worldwide with an annual production of 182.05 million tons [22]. Biotic and abiotic stress affects a significant amount of tomato production, which directly affects fruit development and ripening [23]. Some of the major biotic stresses reported in tomatoes are bacterial wilt and spot, cucumber mosaic virus, cotton bollworm, nematodes and so on [21]. Similarly, some of the reported abiotic stresses in tomatoes are temperature, light, salt, cold, and drought [24]. It has been estimated that around 70% of yield losses yearly are due to abiotic and biotic stress in tomato [24]. It is a rich source of nutrients such as lycopene, lutein, zeaxanthin, potassium, ascorbic acid, and β-carotene [25]. Advanced biotechnological tools have been introduced to new dimensional omics disciplines such as transcriptomics, proteomics and metabolic studies. Several genome-wide characterizations of genes were reported after the release of the tomato genome sequence in 2012 [26].

The genome-wide identification and characterization of gene families will provide the following parameters: chromosomal distribution, number of introns and exons, phylogenetic relationship, gene duplication, expression profiles and so forth that can be used for tomato crop improvement [27]. A genome-wide characterization of the E2F/DP gene family has been performed in a few plant species. Nevertheless, no such efficient investigation of the evolutionary relationships or characteristics of the E2F/DP family genes has been performed in solanaceous crops. Therefore, in the present study, we studied the gene expression of SlE2F/DP genes in different organs and against five abiotic stresses at different time points in tomato. In addition to the gene structure, cis-acting elements, miRNA prediction in SlE2F/DP proteins, subcellular localization assays and gene co-expression network analysis using RNA sequencing data were also carried out to understand their probable functional role in tomato. Our findings aimed to provide a basis for further functional verification of potential E2F/DP transcription factors in tomatoes for crop improvement programs.

Results

In silico identification ofE2F/DP in tomato

This study identified eight E2F/DP tomato genes and named them SlE2F/DP1–SlE2F/DP8 based on their chromosomal locations. The ORFs of SlE2F/DP showed significant variation in length, ranging from 588 bp (SlE2F/DP6) to 2554 bp (SlE2F/DP5), with a mean of 1309 bp. Similarly, the lengths of the E2F/DP proteins varied from 195 to 847 amino acids (aa) for SlE2F/DP6 and SlE2F/DP5, respectively, with a mean of 435 aa. In addition, the predicted molecular weights (MW) vary from 22.2 to 93.2 kDa. The isoelectric point (pI) values of SlE2F/DP proteins varied between 4.68 and 8.92, indicating that E2F/DPs may be acidic or basic. The GRAVY values ranged from − 0.505 to − 0.734, for SlE2F/DP proteins indicating that all these proteins are hydrophilic in nature. Additional information such as subcellular location, chromosome location, and ORF, of these SlE2F/DP family members is given in Table 1.

Table 1 List of identified tomato E2F/DP gene family members with their corresponding encoded protein information

Phylogenetic analysis of tomato E2F/DP family proteins

The protein sequence homology ranged from 9 to 71% among SlE2F/DP proteins (Table S1), and those proteins that were clustered in the same phylogenetic groups had higher sequence similarity among the nine different plant species (Table S13). The phylogenetic analysis classified diverse plant E2F/DP family proteins into seven different groups based on their phylogenetic association (Fig. 1). Mainly, all the proteins of monocots (stiff brome, maize, rice, and wheat) and dicots (black cottonwood, tomato, soybean, Arabidopsis, and potato) split into two different groups in the phylogenetic analysis. The eight tomato E2F/DPs were distributed in four different groups, with a similar number of proteins (SlE2F/DP8 and SlE2F/DP4) in group I and another two proteins (SlE2F/DP1 and SlE2F/DP5) in group III. The other two groups had two proteins each: group IV with SlE2F/DP6 and SlE2F/DP7, and group VI confined with SlE2F/DP2 and SlE2F/DP3 proteins. All SlE2F/DP protein members preferentially clustered with their homologs from potato, which is an evolutionarily closely related species. All the dicotyledonous species were clustered into groups 1, 2 and 6, and monocot homologs were clustered into groups 3, 5 and 7. In contrast, group four clustered both monocot and dicot E2F/DP homologs in the evolutionary relationship.

Fig. 1
figure 1

Phylogenetic analysis of E2F/DP proteins from tomato and different plant species. The neighbor-joining tree was constructed with the full length E2F/DP proteins using ClustalW and MEGA11 with 1000 bootstrap replicates. A species abbreviation was provided prior to each E2F/DP protein name: Sl, Solanum lycopersicum; St, Solanum tuberosum; At, Arabidopsis thaliana; Gm, Glycine max; Os, Oryza sativa; Bd, Brachypodium distachyon; Zm, Zea mays; Pt, Populus trichocarpa and Ta, Triticum aestivum

Gene structure, conserved motif and domain analysis of tomato E2F/DP genes

Exon-intron analysis revealed structural divergence in the tomato E2F/DP gene family (Fig. S1). The number of exons present in the SlE2F/DP genes varied from seven to twenty-eight, with a mean of thirteen. However, the exon-intron composition of the E2F/DP genes within the same phylogenetic group was similar. Likewise, most of the SlE2F/DP homologs clustered in the same phylogenetic groups contained similar conserved motif distributions, (Fig. S2, Fig. S3). This study also identified several motifs unique to some E2F/DP members, such as motif 5, motif 8 and motif 9 for SlE2F/DP2 and SlE2F/DP3, motif 4 for SlE2F/DP6 and SlE2F/DP7, and motif 2 and motif 10 for SlE2F/DP1, SlE2F/DP4, SlE2F/DP5 and SlE2F/DP8 genes. The conserved motifs predicted among the homologous proteins from three different species (Arabidopsis, tomato, and rice) were clustered in the same phylogenetic groups. Motifs 1 and 6 were present in all E2F/DP homologous proteins of Arabidopsis, tomato, and rice, indicating that these motifs were conserved in both dicot and monocot plants during evolution. The domain architecture of SlE2F/DP genes were assorted, changing from the presence of the single E2F_TDP domain to the composition of other additional domains, such as the E2F_CC-MB domain, E2F_DD domains, DP domains and DP_DD superfamily domain (Fig. 2). The E2F_CC-MB domains were identified only in the SlE2F/DP4 and SlE2F/DP5 genes. Similarly, another domain, E2F_DD, was identified only in the SlE2F/DP1 and SlE2F/DP8 genes. Other types of domains, DP and DP_DD, were only identified in the SlE2F/DP6 and SlE2F/DP7 genes. Sequence alignments of SlE2F/DP family members with the human homolog E2F8 protein revealed that the E2F/DP family proteins from a human and a dicotyledonous tomato plant shared five α-helices, β-turns labeled as TT and two β-sheets (Fig. 3).

Fig. 2
figure 2

Schematic depiction of the domain organization of SlE2F/DP proteins. E2F_TDP domain, E2F_CC-MB domain, E2F_DD domains, DP domains and DP_DD superfamily domain identified in SlE2F/DP proteins are shown

Fig. 3
figure 3

Alignment of E2F/DP domains of SlE2F/DP proteins with that of the typical human E2F8. The secondary structural elements determined by the Espript 3.0 web tool are indicated above the alignment

Chromosomal position, gene duplication, and microsynteny analysis of SlE2F/DP genes

SlE2F/DP genes were evenly distributed on 8 of the 12 chromosomes (Chr) of tomato, and they were placed near the distal portions of the chromosomes (Fig. S4). Two segmentally duplicated gene pairs (SlE2F/DP1 and SlE2F/DP5 and SlE2F/DP2 and SlE2F/DP3) were predicted within the tomato genome, one duplicate pair clustered into phylogenetic group VI, whereas other pairs distributed in group II and group I (Fig. 1, Fig. S5). However, all the duplicated genes of each pair were located on different chromosomes with a single SlE2F/DP gene. This study did not detect tandem duplication events because all the genes were mapped on different chromosomes. The Ka/Ks ratio was less than one for the two duplicated gene pairs in the tomato E2F/DP gene family suggested that these genes duplicated due to intense purifying selection over the course of evolution. Similarly, the Ka/Ks ratio was less than one for the duplicated gene pairs between tomato and Arabidopsis (Table 2). It has been estimated that around 13.15–14.76 million years ago these duplicated gene pairs separated. The time estimation was ranged between 49.51 and 142.41 million years ago for the duplicated gene pairs of Arabidopsis and tomato. The comparative microsynteny analysis was also conducted among Arabidopsis, tomato, and rice E2F/DP orthologous genes to distinguish the evolutionary correlation across the genomes (Fig. 4). This analysis identified six orthologous gene pairs between Arabidopsis and tomato, whereas three pairs between tomato and rice and three pairs between Arabidopsis and rice.

Table 2 Predicted Ka/Ks ratio of the duplicated SlE2F/DP and AtE2F/DP gene pair along with its divergence time
Fig. 4
figure 4

Microsyntenic relationship of E2F/DP genes across Arabidopsis, tomato, and rice. The chromosomes of the three species are depicted by different colors: Arabidopsis, pink; tomato, yellow; and rice, blue. All chromosomes are showed with the scale in megabase pairs (Mbp). The duplicated SlE2F/DP genes in tomato genome are indicated by blue lines,

Prediction of Cis-regulatory elements, microrna (mirna) target sites and phosphorylation sites

Hormonal regulation is essential for activating gene expression to control plant stress responses [28]. This study identified many cis-regulatory elements related to phytohormones and abiotic stress responses in the promoter regions of tomato E2F/DP genes. They are drought‐responsive MYB‐binding site (MBS) and MYB elements, stress/defense‐related elements (TC‐rich repeats), jasmonic acid‐responsive elements (CGTCA‐motif), related to the SA response (TCA‐elements), auxin responsive elements (TGA‐element and TGA‐box), ABA‐responsive elements (ABRE), low‐temperature‐responsive elements (LTR), and GA responsive elements (P‐box, GARE‐ and TATC) (Fig. S6, Table S2). In addition to cis-acting elements, miRNAs also regulate stress-related genes against biotic and abiotic stress in plants [29, 30]. This study identified six miRNA target sites related to abiotic stress tolerance in tomato in SlE2F/DP1, SlE2F/DP2, SlE2F/DP4 and SlE2F/DP8 genes. The remaining four genes were not predicted to have miRNA target sites (Table S3). Post‐translational regulation of stress‐related proteins via phosphorylation is important in plant stress responses. This study identified several phosphorylation sites, including CKI and CKII, in tomato SlE2F/DP proteins (Table S4).

Comparative modeling of tomato E2F/DP proteins

All these SlE2F/DP proteins were modeled for the 3D protein structures using the online server I-TASSER along with the C-score, sequence similarity percent, coverage and binding sites (Tables S5 and S6). Using the binding sites, rotating 3D models were created using the software Discovery Studio v.21.1. The number of secondary structural components of tomato E2F/DP proteins was 5–17 for α-helixes, 6–14 for β‐strands, and 13–29 for coils (Table S5). However, the number of secondary structures differed, one protein had few α‐helixes, β‐strands and coils (5, 7 and 13 for SlE2F/DP6), whereas, other proteins had high α‐helixes and coils (17 and 29 for SlE2F/DP2). Along with the secondary structure, the C‐scores and TM‐scores were predicted for all the protein models in Tables S6 and Table S7 [31]. This study identified the putative ligand‐binding sites that can interact with diverse molecules in all predicted models (Fig. 5). To understand the molecular functions of SlE2F/DP family members, gene ontology (GO) terms were used to determine the binding ability to a variety of ligands, transporters and transferase activity using the online server I‐TASSER (Table S8).

Fig. 5
figure 5

Predicted three-dimensional homology structure of tomato E2F/DP proteins. The final 3D structures of SlE2F/DP proteins were built by Discovery Studio v.21.1. The secondary structural components: α‐helices (red), β‐sheets (cyan), coils (green), and loops (gray) as well as the top four putative binding sites: site 1 (yellow sphere), site 2 (green sphere), site 3 (red sphere), and site 4 (blue sphere) are indicated in the predicted 3D models of (A) SlE2F/DP1; (B) SlE2F/DP2; (C) SlE2F/DP3; (D) SlE2F/DP4; (E) SlE2F/DP5; (F) SlE2F/DP6; (G) SlE2F/DP7 and (H) SlE2F/DP8

Co-expression analysis

Co-expression profiling was performed for the SlE2F/DP genes using RNA sequencing data. The results showed that a total of 166 genes were co‐expressed with SlE2F/DP genes (Fig. 6). The hub genes, SlE2F/DP3, SlE2F/DP5, SlE2F/DP7, and SlE2F/DP8 were co‐expressed with 102, 7, 24, and 33 genes, respectively. The KEGG enrichment analysis of these co-expressed genes was concerned with biological pathways, including nucleotide excision repair, pyrimidine metabolism, base excision repair, amino sugar and nucleotide sugar metabolism, basal transcription factors, biosynthesis of secondary metabolites, one carbon pool by folate, phenylpropanoid biosynthesis and Folate biosynthesis (Fig. S7, Fig. S8, Table S9). However, some of the co-expressed genes were not annotated in any of these pathways. In this study, a few of the defence response genes namely HSP (Solyc01g100230), NBS-LRR (Solyc08g005510), F-box-LRR (Solyc10g085600), NAC domain (Solyc02g036430), WD40 domain (Solyc03g114690), and B3 domain (Solyc06g007530) were co-expressed with the SlE2F/DP3 gene. Furthermore, the SlE2F/DP5 gene was co-expressed with a set of genes such as thymidylate synthase (Solyc01g109830), which is associated with tomato fruit ripening and AEP (Solyc12g098160). Similarly, another HSP20 gene (Solyc01g102960) and GDSL (Solyc11g006250) gene involved in plant growth and development were co-expressed with the SlE2F/DP7 gene. Likewise, a few transcription factors namely zinc finger C3HC4 and CCCH domain (Solyc05g051800 and Solyc06g054600), serine/threonine gene (Solyc09g011320) reported against biotic and abiotic stress, cyclin related gene (Solyc04g072880) involved in the signalling control of stress tolerance were co-expressed with the hub gene SlE2F/DP8.

Fig. 6
figure 6

Weighted gene co-expression network analysis (WGCNA) of SlE2F/DP genes. A–D: the co‐expressed genes in the network of (A) SlE2F/DP3, (B) SlE2F/DP5, (C) SlE2F/DP7 and (D) SlE2F/DP8. The SlE2F/DP genes are marked in red

Subcellular localization of SlE2F/DP proteins

The SlE2F/DP family genes were predicted for localization analysis using the wolfsport server. The E2F/DP proteins were predicted in various parts of the cell, such as the nucleus, ER, and cytoplasm (Table S10). Furthermore, to validate the insilico predicted results, subcellular localization was carried out by expressing these SlE2F/DP proteins fused with sGFP in rice protoplasts. The results showed that SlE2F/DP4, SlE2F/DP6 and SlE2F/DP7 had florescence signal in the ER whereas, SlE2F/DP3 was localized in the nucleus (Fig. 7).

Fig. 7
figure 7

Subcellular localization of SlE2F/DP proteins. SlE2F/DPs-SGFP fusion constructs were used to analyze the localization of SlE2F/DP3, SlE2F/DP4, SlE2F/DP6 and SlE2F/DP7 and the fluorescence signals were visualized with the confocal microscope. The NLS-mRFP construct and ER-mRFP was utilized as an Nucleus and ER localization marker. Scale bars = 10 μm

Expression profiling of tomatoE2F/DPgenes in different organs

Further to understand the functional role of these SlE2F/DP genes in different developing tomato organs, expression analysis was conducted (roots, stems, leaves, flowers, fruits: 1 cm, IM, MG, B, and B5) of the Alisa Craig genotype through RT-qPCR. (Fig. 8). The expression levels of eight SlE2F/DP genes in tomato organs revealed differential expression patterns (Fig. 8). The duplicated gene pair, SlE2F/DP1 and SlE2F/DP5, showed high expression levels in roots (> 4-fold) compared to the control (leaf) (Fig. 8). Similarly, another duplicated gene pair SlE2F/DP2 and SlE2F/DP3, also showed higher expression in roots (> 8‐fold) than the control. The SlE2F/DP1 gene, showed higher relative expression in B5 (12 fold change) followed by B fruit (10 fold), roots (7 fold change), 1 cm fruits (5 fold change), MG fruit (3 fold change), IM fruit (3 fold change), stem (2 fold change), but the expression in flowers was lower than that in the control (Fig. 8). Three genes SlE2F/DP2, SlE2F/DP4 and SlE2F/DP5 were significantly upregulated in all the samples (Fig. 8). SlE2F/DP2 had higher expression in the root (9 fold) and B5 (8 fold), followed by 1 cm fruit (5 fold) and B (4 fold), stem and MG (3 fold) and flower and IM (2 fold) compared to its respective control (Fig. 8). SlE2F/DP4 gene expressed highly in 1 cm fruit, root, B5, stem, flower, B, MG, IM by 18, 17,14,10,9, 6.7, 6.2 and 6.2 fold values (Fig. 8). Similarly, SlE2F/DP5 gene had high expression in root (7.7 fold) followed by 1 cm fruit, B5, B, flower, stem, IM and MG (6.6, 6.4, 6.0,4.8, 4,3.1 and 2.3 fold) (Fig. 8). The gene SlE2F/DP3 showed high expression in root sample (11 fold) followed by 1 cm fruit (8 fold), stem (4 fold) compared with their control. However, its expression level was less than the control in other samples. Similarly, SlE2F/DP6 expressed more than 2-fold in the root (4.7) and B fruit (2.0), but the expression was less than the control in other samples. SlE2F/DP7 showed high relative expression in B5, root, B, stem and 1 cm fruit (12, 6.9, 5.8, 3.1 and 2.5) and the remaining samples were less than 1 fold. SlE2F/DP8 gene showed 4.9 fold in the root, and in five samples (stem, 1 cm, MG, B, B5) more than 2 fold change than the control sample.

Fig. 8
figure 8

Expression profiles of SlE2F/DP genes in various organs: leaves, stems, roots, flowers, 1 cm fruits, immature fruits (IM), mature green fruits (MG), breaker fruits (B), and fruits 5 days after the breaker stage (B5). The standard deviations of the means of three independent biological replicates are represented by the error bars. The different asterisks above the bars indicate the significant variations between the control samples (leaves) and the other organs by Student’s t-test with p‐values less than 0.05 for *, 0.01 for **, and 0.001 for ***, respectively

Expression analysis of SlE2F/DP genes against abiotic stresses and phytohormone treatment

The transcript profiles of tomato E2F/DP genes were studied in leaf samples exposed to different abiotic stresses (heat, drought, cold, and salt) and phytohormone (ABA) treatments through qRT-PCR (Fig. 9A–E). Six genes (SlE2F/DP1, SlE2F/DP3, SlE2F/DP4, SlE2F/DP5, SlE2F/DP6 and SlE2F/DP8) were significantly downregulated under heat stress at 40 oC (0.30‐ to 1.44‐fold). Interestingly, two genes SlE2F/DP2 and SlE2F/DP7 were upregulated by 2.30-fold change at 1 h and 3 h under heat stress (Fig. 9A).

Fig. 9
figure 9

Expression profiling of SlE2F/DP genes under different abiotic stresses: (a) heat, (B) salt, (C) cold, (D) drought, and (E) ABA treatment. Error bars represent standard deviations of the means of three independent biological replicates. The various asterisk marks (* for p-value < 0.05, ** for p‐value < 0.01, and *** for p‐value < 0.001) indicate statistically significant compared to its respective control using the student t‐test

Salt stress triggered the expression level of the SlE2F/DP genes at various time-points (Fig. 9B). Compared to the control, SlE2F/DP2 was upregulated (2.12-fold) at 24 h after salt treatment. Beside SlE2F/DP7 gene expression was repressed during the early hours of salt stress but gradually increased after 3 h (3.8 fold) and maintained its expression level until 24 h (Fig. 9B). In contrast, SlE2F/DP8 gene was upregulated by 2.7 fold change at 3 h and gradually decreased at later stages after salt stress.

Seven SlE2F/DP genes were significantly upregulated at several time points against the cold stress (4 oC). SlE2F/DP2 gene showed significant upregulation of 2-fold change at 3 h of exposure to cold stress. Similarly, SlE2F/DP3 and SlE2F/DP5 genes showed upregulation (> 2 fold) at 1 h and 3 h against the cold stress. In contrast, SlE2F/DP4 gene showed > 2 fold upregulation at the early stages and gradually decreased at 24 h under cold stress. Remaining four genes SlE2F/DP1, SlE2F/DP6, SlE2F/DP7 and SlE2F/DP8 were downregulated (0.61 to 1.79 fold) compared to the control (Fig. 9C).

Drought stress-regulated the expression levels of the SLE2F/DP genes. Four genes SlE2F/DP1, SlE2F/DP3, SlE2F/DP4 and SlE2F/DP6 were significantly downregulated under drought stress. Another four genes SlE2F/DP2, SlE2F/DP5, SlE2F/DP7 and SlE2F/DP8 were upregulated by 1.1 to 1.4 fold under drought treatment compared with the respective control (Fig. 9D).

The SlE2F/DP genes expression level was similar to control after the ABA hormone treatment. Only SlE2F/DP7 gene showed 2-fold upregulation at 3 h compared to the control in response to ABA treatment. Six genes (SlE2F/DP1, SlE2F/DP2, SlE2F/DP4, SlE2F/DP5, SlE2F/DP6 and SlE2F/DP8) were expressed from 1.1 to 1.8 fold at various time points under ABA treatment compared with the respective control. In contrast, SlE2F/DP3 gene was downregulated at all the five time points tested (Fig. 9E).

Discussion

E2F/DP transcription factors (TFs) have been identified in both plants and mammals. The classification of plant E2F/DP genes is generally based on their phylogenetic analysis, domains, motifs, and subcellular localization, and their role against abiotic stress has been reported in few plant species [19, 20, 32], but not in any solanaceous crop. The present study identified and characterized eight E2F/DP genes of tomato. The number of E2F/DP genes was observed to vary from one species to another within the plant kingdom. Hence, these genes may play an important biological role in this vegetable crop.

Based on the phylogenetic analysis the E2F/DP homologs mostly from the dicot species, and the monocot species, clustered separately suggested that these proteins may be separated during the divergence of monocots and dicots (Fig. 1). However, in group II, only the tomato and potato proteins were clustered, revealing that these proteins were particular to solanum families among the dicots. In contrast, group IV clustered both the monocot and dicot plant species, suggested that they originated from the common ancestor angiosperm before the divergence of monocot and dicot (200 million years ago) [33]. Notably, all the SlE2F/DP proteins formed a close association with potatoes E2F/DP proteins, underlying the evolutionary conservation of these E2F/DP members in the solanaceous family.

The conserved domain analysis results suggested that all the SlE2F/DP proteins had E2F-TDP domain along with or without additional domains (Fig. 2). These findings suggested that the combination of additional domains may have evolved and differed from the ancestral E2F/DP domains. All these SlE2F/DP genes were analysed for the conserved motifs and gene structure prediction suggested that structural differences may be particular to each gene family (Fig. S1, Fig. S2). The similar motifs of paralogous and orthologous gene pairs from tomato, rice and Arabidopsis suggested that these conserved motifs remain similar even after species evolution between monocots and dicots of SlE2F/DP family genes.

During the evolutionary process, gene duplication events are the primary sources of genetic novelty [34]. This segmental duplication increases the number of genes within a particular gene family, which enhances plants adaptation to various environmental stresses [35, 36]. The two duplication pairs in the SlE2F/DP gene family suggested that segmental gene duplication triggered the expansion of the tomato E2F/DP gene family during evolution. Multiple segmental duplication events were reported in the E2F/DP gene family in wheat [20]. Another study reported that 12 duplicated gene pairs were observed in the E2F/DP gene family in Moso bamboo [32]. In contrast, no gene duplication was observed in the E2F/DP gene family in Medicago truncatula [19].

Microsynteny analysis revealed six segmental gene pairs between tomato and Arabidopsis and divergence time was calculated for dicot (Table 2) indicate that tomato has the closest relationship with the dicotyledonous model plant Arabidopsis, than the monocotyledonous model plant rice (Fig. 4). Interestingly, the three tomato genes were commonly duplicated between tomato and Arabidopsis as well as between tomato and rice. The duplication of gene pairs observed between Arabidopsis and rice suggested that these genes could have been present in Arabidopsis and rice before the separation of ancestral genes during the evolutionary process. Moreover, three tomato genes were commonly duplicated between the dicotyledonous plant Arabidopsis and the monocot rice plant suggested that these genes play an important role in evolutionary changes (Fig. 4).

The cis-regulatory elements in the promoter region has been reported to be associated with stress-related genes against biotic and abiotic stress [37, 38]. In the present study, we also predicted numerous hormone-related and abiotic stress-responsive cis-acting elements in the promoter regions of 8 SlE2F/DP genes (Fig. S6 and Table S2). Low-temperature responsive element (LTR) motifs were present in the promoter regions of SlE2F/DP1 and SlE2F/DP4. All the 8 SlE2F/DP promoter regions had more than two light response (TCT, GT1, GATA, GA, G-box, I-box and ACE) motifs. Similarly, all 23 PheE2F/DP genes had light response motifs in their promoter region [32]. The drought-responsive element MBS motif was detected in the promoters of SlE2F/DP6, SlE2F/DP7 and SlE2F/DP8 suggested their role in abiotic stress. These MBS elements were reported in HVA22 gene and PLATZ transcription factor against abiotic stress [39, 40]. The ABRE cis-elements and the Meja responsive motif is abundant in SlE2F/DP genes. In addition, auxin, salicylic acid and gibberellin-responsive elements were also observed in some of the SlE2F/DP genes, suggested that they may be stimulated by environmental stress. This type of hormone cross-talk regulatory element was also reported in E2F/DP genes of Moso Bamboo [32].

miRNAs can regulate gene expression either by inhibiting transcription and translation or by suppressing transcription and translation [41]. In the current study, six tomato miRNAs (sly-miR172a, sly-miR172b, sly-miR9476-5p, sly-miR9472-3p, sly-miR9471a-3p and sly-miR9471b-3p) were targeted to cleave four SlE2F/DP genes (Table S3). Sly-miR172a and Sly-miR172b were reported to be involved in tomato fruit development, ripening and against salt and drought stress [42, 43]. Similarly, other miRNAs, namely, Sly-miR9476 and Sly-miR9472 were reported to be involved in the regulation of defence genes DRP, HD-ZIP, MYB, and NAC against drought stress [44]. Another sly-miR9472 was reported to be involved in drought-stress through RNA sequencing and qRT-PCR [45]. Interestingly, SlE2F/DP1 and SlE2F/DP8 were predicted to have targets for multiple miRNAs. These results suggested that multiple miRNAs could regulate the gene expression of SlE2F/DP genes in tomatoes.

Protein phosphorylation plays an important role in post-translational modification by altering the protein structure or its function. To activate or deactivate the proteins, the amino acids are phosphorylated by kinases and phosphatases that are involved in signalling pathways and stress responses in plants [46, 47]. This study predicted phosphorylation sites in all tomato E2F/DP proteins (Table S4). Protein kinases are reported against cold and heat stress in tomato, Arabidopsis, wheat, and soybean [48, 49]. The phosphorylation sites S686, T688, S689, and T705 act as targets for mutagenesis and are required to maintain Xa21 protein stability of rice bacterial blight disease [50]. Similarly, phosphorylation sites T867, S878, T1040, and T1072 is also essential for the function of the FLS2 gene against bacterial Pseudomonas syringae in Arabidopsis [51]. In the recent years, around 1098 phosphorylation sites were reported in response to drought or heat stress in the grass species [52]. However, the exact role of these phosphorylation sites can be elucidated after conducting additional experimental studies.

To understand the possible interaction with other molecules and ligand-binding sites, a three-dimensional structure of SlE2F/DP proteins was predicted (Fig. 5). Most of the E2F/DP proteins had a higher number of coils than α-helixes and β‐strands (Table S5) suggested its important role in the molecular structure during evolution. Although it has yet to be experimentally validated, the SlE2F/DP proteins were predicted to have a 3D modelling structure with their binding molecules such as ligand binding sites, ions, including peptide molecules and DNA that can change cell function.

The transcription factors are reported to be localized in the nucleus. The subcellular localization results showed that only SlE2F/DP3 localized in the nucleus and three proteins SlE2F/DP4, SlE2F/DP6 and SlE2F/DP7, were predominantly localized in the ER (Fig. 7). Even though, TFs are reported to be localized in the nucleus, some TFs (bZIP and NAC) families are not localized in the nucleus during their initial synthesis [53]. In the genome wide characterization of MYB TFs in Physcomitrella paten, one gene was localized in the ER [52]. Similarly, the CREB3 (bZIP) family of TFs are ER-localized proteins [54]. In contrast, to our results E2F/DP proteins were expressed in the nucleus and cytoplasm in wheat [20]. These results suggested that SlE2F/DP3 may act as a transcription repressor and that other E2F/DP proteins may act as co-activators which is required for E2F activity or in regulating signalling pathways [55].

Gene expression studies of SlE2F/DP family genes were tested in various tissues to understand their functional role and diversity in tomatoes. Interestingly, many of the E2F/DP family genes namely SlE2F/DP3, SlE2F/DP6 and SlE2F/DP8 were predominately expressed in the vegetative tissues (roots), while other genes, SlE2F/DP1 and SlE2F/DP7 were highly expressed in the reproductive tissues (B5), suggested its importance in the tissue-specific functions and development in tomato. The remaining three genes SlE2F/DP2, SlE2F/DP4 and SlE2F/DP5 showed upregulation in both vegetative and reproductive tissues. All the eight tomato E2F/DP genes that were highly expressed in roots and stems, suggested their involvement in root growth is to uptake and supply nutrition and water to other parts of plants. A similar trend was observed in most of the E2F/DP gene expression studies in Medicago truncatula [19]. In contrast, the most of the E2F/DP genes are highly expressed in leaf tissues in wheat [20].

The role of E2F/DP genes in fruit development has not been studied in any vegetable crop, therefore our study is the first analysis of E2F/DP genes in the reproductive stages of tomato. Out of the eight E2F/DP genes, four genes namely SlE2F/DP1, SlE2F/DP2, SlE2F/DP4 and SlE2F/DP5 were expressed in both 1-cm and breaker stages suggested that these genes may be involved in both cell division and ripening processes. Another three genes namely SlE2F/DP1, SlE2F/DP2 and SlE2F/DP7 were highly expressed in the B5 followed by the B revealing that these genes may play an important role in the fruit ripening processes. The higher expression levels (≥ 5 fold change) of SlE2F/DP1, SlE2F/DP3, SlE2F/DP4 and SlE2F/DP5 in 1-cm fruit suggested that they may be involved in the cell division phase of tomato fruit development. In another study, SlE2F/DP3 (Solyc03g113760) was reported to be involved in tomato fruit development [56]. In contrast, SlE2F/DP6 gene showed lower expression levels in all developmental stages of fruit. These results suggested that E2F/DP genes play a regulatory role in fruit development, enlargement and ripening, which may lay the foundation for functional validation of this gene family in tomatoes.

Plants have been exposed to many wide ranges of environmental stress and biotic stress. Plant defence against abiotic stress was regulated by a complex network of genes including numerous transcription factors, defence genes along with its downstream genes reported to be involved in the ABA-mediated pathway [57]. A few reports of E2F/DP transcription factors in plant abiotic stress have been reported in different plant species [19, 20, 32].

In the current study, tomato SlE2F/DP genes displayed different expression levels against different abiotic stress (Fig. 9A-E). Out of the eight genes tested, two genes showed significant upregulation against heat stress (Fig. 9A), which agrees with the earlier finding of E2F/DP genes in wheat in response to heat stress [20]. Salt stress triggered the expression level in three genes SlE2F/DP2, SlE2F/DP7 and SlE2F/DP8 (Fig. 9B). These results suggested their possible role in salt tolerance in tomato. Similarly, most of the E2F/DP genes were expressed significantly in Moso bamboo and barrel clover upon exposure to salt stress [19, 32]. In contrast, the E2F/DP genes were expressed less than two-fold at different time points in wheat [20]. The expression level of four SlE2F/DP genes (SlE2F/DP2, SlE2F/DP3, SlE2F/DP4 and SlE2F/DP5) were induced following cold stress (Fig. 9C). These results were consistent with previous reports of three wheat E2F/DP genes (TaE2F11-9, TaDP3111-15 and TaDEL211-27) in response to cold stress [20]. On exposure to drought stress, the expression level of SlE2F/DP genes was induced ≤ 1.2 fold at any time point. Our results were consistent with a previous report of E2F/DP genes in wheat [20]. The phytohormone ABA is well-studied and known for its role in regulating many stress related genes against various environmental stresses [58, 59]. Our results showed that SlE2F/DP7 was upregulated only at one-time point in tomato plants following ABA treatment (Fig. 9E) suggested its possible role in drought tolerance. To understand the defence mechanism against rice blast infection, 14 different TFs were used for that study. Interestingly, the E2F/DP TF gene of rice (LOC_Os12g06200.1) was upregulated at the early stage against the blast infection [60]. Till now, no reports of E2F/DP genes treated with ABA phytohormones in other plant species. Besides, key cis-regulatory elements related to hormones-responsive, including ABA, auxin, SA, and jasmonic acid were observed in upstream site of SlE2F/DP genes. It was stated that phytohormones such as ABA, SA, and jasmonic acid are associated with plant cell signalling pathways linked to response of biotic and abiotic stresses [61]. Key stress signalling pathways and ion transporter genes are controlled by ABA signalling in response to stresses [62]. Our results disclose that SlE2F/DP genes are probably being in cellular cascades related to ABA.

To gain a better understanding of the putative functions of SlE2F/DP genes, gene co-expression network analysis of SlE2F/DP was performed using the RNA-seq data (Fig. 6). Several genes involved in the co-expression network interaction were related to abiotic stress tolerance. SlE2F/DP3 co-expressed genes such as SlHSP SlLRR, SlF-box-LRR, NAC domain gene and B3 DNA binding domain were reported against biotic and abiotic stress [19, 34, 63, 64]. Similarly, another gene thymidylate synthase reported to control the expression of folate biosynthesis genes in tomato fruit was co-expressed with the SLE2F/DP5 gene [65]. In addition, another SlHSP20 involved in heat tolerance in tomato and SlGDSL gene regulates fruit development in tomato [66, 67] were detected in the co-expression analysis of the SlE2F/DP7. Likely, SlE2F/DP8 gene established a co-expression network with several genes such as Zinc finger C3HC4, CCCH domain, serine/threonine and cyclin-related gene are reported to be involved in both biotic and abiotic stress and cell division in tomato [68,69,70]. These results suggested that these transcription factor E2F/DP genes may play an important role in abiotic stress tolerance in tomato.

Conclusion

In the present study, we identified eight SlE2F/DP genes in tomato. These genes were clustered into seven phylogenetic groups based on their domain, gene structure and motif analyses. Gene duplication analyses spotted the two segmental gene duplication events were responsible for the expansion of the E2F/DP gene family throughout the evolution in tomato. We hypothesize that tomato E2F/DP genes play an important role in development process. In addition, E2F/DP genes were associated with cis-regulatory elements, miRNA target sites and phosphorylation sites in their sequences in regulating gene expression. Expression profiling showed that most of the SlE2F/DP genes are highly expressed in both vegetative and reproductive organs. Many of the SlE2F/DP genes were significantly induced under different abiotic stress concluding that their possible function in tolerance against this abiotic stress. Co-expression network analysis highlighted the biological importance of SlE2F/DP genes in tomato and predicted the hub genes co-expressed with SlE2F/DP genes revealing their possible role in development and abiotic stress tolerance. Our findings provide the valuable information for the further functional exploration of the SlE2F/DP gene family and identify the potential candidate genes for the genetic improvement of tomato.

Methods

Identification and sequence analysis of E2FDP genes in the tomato genome

We identified the E2F/DP genes in the Plant Transcription Factor Database (http://www.Planttfdb.gao-lab.org) using the keyword search “E2F/DP’’ [71]. Simultaneously, E2F/DP protein sequences were downloaded from the TAIR database (https://www.Arabidopsis.org) of Arabidopsis thaliana. We performed a blast search using AtE2F/DP protein sequences in the Sol genomics database (http://www.solgenomics.net/tools/blast) with its default parameters [72]. The HMM profile of SlE2F/DP1 (PF16421) was taken from the Pfam (http://pfam.xfam.org/). The eight SlE2F/DP protein sequences were used for checking the presence of domains using the two online webserver SMART tool (http://smart.emblheidelberg.de/) and NCBI CDD webtool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi). It is essential to observe the molecular weight, protein length, grand average of hydropathicity index (GRAVY), and isoelectric points of the SlE2FDP proteins using the online webserver Expasy (http://cn.expasy.org/tools/protparam.html) [73]. To visualize the exon-intron structure of SlE2FDP genes the genomics and coding sequence data were loaded to the Gene Structure Display Server-2.0 (GSDS-2.0) web server (http://gsds.cbi.pku.edu.cn/) [74]. The open reading frames of the SlE2FDP genes were analysed using the Open Reading Frame Finder tool (https://www.ncbi.nlm.nih.gov/orffinder/). To identify the conserved domains, full length protein sequences of E2F/DP from Arabidopsis, tomato, and rice were analysed using the Multiple EM for Motif Elicitation (MEME) web server (http://meme-suite.org/) with the parameter of 10 motifs with a length of 6 to 50 [75]. Multiple alignment was performed using the protein sequences of SlE2FDP genes by the online webtools Clustal Omega and ESPript [76, 77]. The subcellular localization was analysed for SlE2FDP proteins using the webtool “WoLF-PSORT” (https://wolfpsort.hgc.jp/) [78]. The sequence homology of SlE2FDP proteins was analysed using the webserver “Immunomedicine Group” (http://imed.med.ucm.es/Tools/sias.html).

Phylogenetic analysis of tomato E2F/DP proteins

The E2F/DP protein sequences were downloaded from the nine different plant species and aligned using the online server ClustalW and then the aligned file was used for phylogenetic analysis using the method the neighbour-joining (NJ) with the 1000 bootstrap replications in MEGA 6.0 software [79]. The nine plant species sequences of Arabidopsis thaliana, Brachypodium distachyon, Glycine max, Oryza sativa, Populus trichocarpa, Solanum lycopersicum, Solanum tuberosum, Triticum aestivum, and Zea mays sequences were taken from the Phytozome database http://www.phytozome.net. All the sequences, gene names and accession number details used in this study are given in Table S13.

Chromosomal locations, gene duplication, and microsynteny analysis

For all the eight SlE2F/DP genes the chromosomal locations were collected from the Sol genomic database (http://www.solgenomics.net). With the gene physical position and chromosome details these genes were mapped on their respective chromosomes using the online server MapGene2Chrom web v2 (http://mg2c.iask.in/mg2c_v2.0/). The duplicated gene pairs among the eight SlE2F/DP genes were predicted using the one-step MCScanX program of TBtools software [80] and followed by BLASTP program with an E-value < 10 − 10. The duplicated gene pairs, synonymous (Ks) and non-synonymous (Ka) ratio were calculated with the Nei and Gojobori (NG) program of TB tool software [81]. Furthermore, the mode of selection was calculated based on the Ka/Ks ratio [80]. This duplicated gene pair’s divergence time (T) was calculated with the standard formula T = Ks/2r MYA (millions of years ago). Ks-stands for synonymous substitution rate per site and r stands for 1.5 × 10 − 8 substitutions per site/year for all the dicot species [81]. The E2F/DP genes from the tomato, Arabidopsis, and rice was performed for the synteny analysis using the reciprocal BLAST search program of TB tool software. The duplicated gene pairs across the three genomes were visualized using the circos program of the TB tool software [80].

Prediction of cis-acting elements, mirna target sites, phosphorylation sites

The putative miRNA present in the coding region of E2F/DP genes was predicted using the online tool psRNATarget (http://plantgrn.noble.org/psRNATarget/analysis). Likewise, the cis-acting elements existing in the promoter (1500 bp) of SlE2F/DP genes were analysed using the online webserver Plant-Care database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) [82]. The phosphorylation sites (Ser/Thr/Tyr) was predicted in the SlE2F/DP genes using the online web-based tool NetPhos 3.1 [83].

Homology modelling of SlE2F/DP proteins

The 3D structure prediction and function of the eight SlE2F/DP proteins was analysed using the web server I-TASSER. These eight SlE2F/DP 3D models were generated using multiple threading alignments with the iterative TASSER and LOMETS assembly programs of the I-TASSER web server [29]. Based on the template, the protein analogs were identified and based on the maximum score and homology in PDB with known structures, the best 3D structures were selected and refined using the ModRefiner program [ 84]. The 3D structure of SlE2F/DP proteins along with their ligand-binding sites were visualized using Discovery Studio v.21.1 software.

Plant sample collection and stress treatments

Tomato seeds of the Ailsa Craig genotype (Solanum lycopersicum L.) were collected from the Giovannoni laboratory at the Boyce Thompson Institute. The seeds were germinated in the nursery tray and maintained in a growth chamber under the optimum conditions at 25°C day/20°C night, with 16-h light/8-h dark photoperiod, 55–70% humidity, and 300 µmol m− 2S− 1 of light intensity. Twenty-eight-day old seedlings were used for collecting fresh root, stem, and leaf samples. The remaining plants were moved to the greenhouse under the controlled conditions at 25/20°C day/night temperatures. When the plants reached the reproductive stage the flower and fruit samples were collected for organ specific expression studies. The fruit sample was collected at different development stages (i) after 7 days of pollination with 1 cm diameter (1 cm fruits), (ii) after 21 days of pollination immature fruits (IM fruits), (iii) after 35 days of pollination mature green (MG fruits), (iv) when the fruits reach the breaker stage (B fruits), and (v) after 5 days of breaker stage (B5 fruits) [85].

When the tomato seedlings were 28 days old, the plants were treated with different abiotic stress treatments such as abscisic acid (ABA), cold, drought, heat and salt (NaCl). The leaf samples were collected at five different time points 0, 1, 3, 9, and 24 h after heat at 40˚C, cold at 4˚C, salt (200 mM NaCl) and 100 µM ABA spraying on the leaf samples [37]. In the case of drought treatment, the tomato plants were allowed to withhold water and leaf samples at different time points 0, 24, 48, 60, and 72 h were collected [ 84]. For the control samples (0 h) the tomato seedlings were grown in soil under normal conditions at (25°C) against heat, cold, drought, and ABA stress [ 85]. All the samples were collected for three biological replicates, crushed in liquid nitrogen, and stored at − 80°C.

Expression profiling of SlE2FDP genes by qRT − PCR

Total RNA was isolated from treated and control samples using an RNeasy Mini kit (Qiagen, Germany), as per the manufacturer’s guidelines. The concentration and the purity of RNA samples were checked using NanoDrop® 1000 spectrophotometer (Wilmington, USA). cDNA was synthesized for all the samples with Superscript® III First-Strand kit (Invitrogen, USA) as per the manufacturer’s guidelines. The primers for SlE2F/DP genes were designed using Primer3 online software http://frodo.wi.mit.edu/primer3/input.html (Table S12), and each primer pair was validated for specificity using a melting curve analysis [82]. The 18 S rRNA (F: AAAAGGTCGACGCGGGCT, R: CGACAGAAGGGACGAGAC) gene from tomato was used as a normalization gene [86]. qRT-PCR reaction volume of 10 µL contained 1 µL (50 ng) cDNA, 2 µL of 5 pmol concentration primers, 5 µL of iTaq SYBR Green (Qiagen, Hilden, Germany), and 2 µL of double distilled water. The PCR was performed using Quantstudio 5 ® 96 (Applied Biosystems, USA) with following conditions: 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 relative expression values were calculated using the 2−∆∆Ct method [87].

Co-expression network analysis of SlE2FDP genes

The RNA-seq data were downloaded from NCBI database (https://www.ncbi.nlm.nih.gov/sra) for the co-expression analysis. The following accessions SRR7652567, SRR7652566, SRR7652565, SRR7652564, SRR7652571, SRR7652570, SRR7652569, SRR7652568, SRR7652563, SRR7652562, SRR7652573, SRR7652572, SRR12026415, SRR12026416, SRR12026417, SRR12026418, SRR12026419, SRR12026420, SRR12026421, SRR12026422, SRR12026423, SRR12026420, SRR12026425, SRR12026426, SRR21847247, SRR21847273, SRR21847264, SRR21847257, SRR21847271, SRR21847272, SRR21847267, SRR21847266 and SRR21847265 were downloaded. The raw sequence quality was checked using the FastQC tool and the low-quality and adapter sequences were removed if any [88]. The good quality reads were mapped on the tomato genome ITAG4.0 using HISAT2 2.1.0 Galaxy platform web server [89, 90]. The reads mapped on the exon were counted using FeatureCounts [91]. The expression of the reads was analysed as fragments per kb per million reads (FPKM) by using DESeq2 software [92]. The co-expression network of SlE2FDP genes and WGCNA using the IDEP web tool with genes greater than 1 FPKM [93]. The visualization was constructed using Cytoscape (https://cytoscape.org/). The top ten co-expressed genes were used for GO and KEGG analysis [94], using the KOBAS web tool http://kobas.cbi.pku.edu.cn/.

Subcellular localization

The coding regions of SlE2F/DP genes were amplified with gene-specific primers (Table S11) and then cloned into the sGFP-tagged vector pGA3452 driven by the maize ubiquitin1 promoter to generate transiently expressed SlE2F/DP–sGFP fusion proteins [95]. A vector expressing NLS-mRFP and OsAsp1-mRFP fusion proteins served as nuclear and ER markers, respectively [96]. The SlE2F/DP–sGFP fusion construct and the NLS-mRFP and OsAsp1-mRFP constructs were introduced into rice cell protoplasts using PEG-mediated transformation [97]. The samples were incubated overnight at 28°C for a period of 12 to 16 h in the dark condition. After incubation the samples were observed under a confocal microscope (BX61; Olympus, Tokyo, Japan) for GFP and RFP channel fluorescent signals.

Statistical analysis

The data were analysed using Quantstudio 5 ® 96 inbuild software (Applied Biosystems, USA) with its default baseline and threshold values. Three replicate mean values were tested with Student’s t-tests on MS-Excel 2010 software. Asterisks (*, **, and ***) indicate significance at p-values < 0.05, < 0.01, and < 0.001, respectively.

Data availability

We declare that the dataset(s) required to reproduce the results of this article are included in the article and additional file(s) available in the journal webpage.

Abbreviations

E2F/DP :

Eukaryotic 2 Transcription Factor/Dimerization Partner

3D:

Three Dimensional

GO:

Go Ontology

TF:

Template Modeling

RMSD:

Root-Mean-Square Deviation

C-score:

Confident Score

ABA:

Abscisic Acid

WGCNA:

Weighted Gene Co-expression Network Analysis

References

  1. Xiao J, Cheng H, Li X, Xiao J, Xu C, Wang S. Rice WRKY13 regulates cross talk between abiotic and biotic stress signaling pathways by selective binding to different cis-elements. Plant Physiol. 2013;163:1868–82. https://doi.org/10.1104/pp.113.226019.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Yokotani N, Sato Y, Tanabe S, Chujo T, Shimizu T, Okada K, Yamane H, Shimono M, Sugano S, Takatsuji H, Kaku H. WRKY76 is a rice transcriptional repressor playing opposite roles in blast disease resistance and cold stress tolerance. J Exp Bot. 2013;64(16):5085–97. https://doi.org/10.1093/jxb/ert298.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Walley JW, Dehesh K. Molecular mechanisms regulating rapid stress signaling networks in Arabidopsis. J Integr Plant Biol. 2010;52:354–9. https://doi.org/10.1111/j.1744-7909.2010.00940.x.

    Article  CAS  PubMed  Google Scholar 

  4. Vinocur B, Altman A. Recent advances in engineering plant tolerance to abiotic stress: achievements and limitations. Curr Opin Biotechnol. 2005;16:123–32. https://doi.org/10.1016/j.copbio.2005.02.001.

    Article  CAS  PubMed  Google Scholar 

  5. Atkinson NJ, Urwin E. The integration of plant biotic and abiotic stresses: from genes to the field. J Exp Bot. 2012;63:3523–44. https://doi.org/10.1093/jxb/ers100.

    Article  CAS  PubMed  Google Scholar 

  6. Liu JH, Peng T, Dai W. Critical cis-acting elements and interacting transcription factors: key players associated with abiotic stress responses in plants. Plant Mol Biol Rep. 2014;32:303–17. https://doi.org/10.1007/s11105-013-0667-z.

    Article  CAS  Google Scholar 

  7. Zhang J, Liu JF, Zhao TT, Ren J, Xu XY. Research progress of plant C2H2 zinc finger protein. Mol Plant Breed. 2018;29:1133–41. https://doi.org/10.3969/gab.029.001133.

    Article  Google Scholar 

  8. Shen Z, Lin Y, Zou Q. Transcription factors–DNA interactions in rice: identification and verification. Brief Bioinf. 2020;21(3):946–56. https://doi.org/10.1093/bib/bbz045.

    Article  CAS  Google Scholar 

  9. Sidhu NS, Pruthi G, Singh S, Bishnoi R, Singla D. Genome-wide identification and analysis of GRAS transcription factors in the bottle gourd genome. Sci Rep. 2020;10(1):1–9. https://doi.org/10.1038/s41598-020-71240-2.

    Article  CAS  Google Scholar 

  10. Ueda K, Matsuura H, Yamaguchi M, Demura T, Kato K. Genome-wide analyses of changes in translation state caused by elevated temperature in Oryza sativa. Plant Cell Physiol. 2012;53:1481–91. https://doi.org/10.1093/pcp/pcs092.

    Article  CAS  PubMed  Google Scholar 

  11. Ramirez-Parra E, Lopez-Matas MA, Frundt C, Gutierrez C. Role of an atypical E2F transcription factor in the control of Arabidopsis cell growth and differentiation. Plant Cell. 2004;16(9):2350–63. https://doi.org/10.1105/tpc.104.023978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Polager S, Ginsberg D. E2F - at the crossroads of life and death. Trends Cell Biol. 2008;18(11):528–35. https://doi.org/10.1016/j.tcb.2008.08.003.

    Article  CAS  PubMed  Google Scholar 

  13. Ramirez-Parra E, Friindt C, Gutierrez C. A genome-wide identification of E2F-regulated genes in Arabidopsis. Plant J. 2003;33(4):801–11. https://doi.org/10.1046/j.1365-313X.2003.01662.x.

    Article  CAS  PubMed  Google Scholar 

  14. Sozzani R, Maggio C, Varotto S, Canova S, Bergonioux C, Albani D, Cella R, Interplay. between Arabidopsis activating factors E2Fb and E2Fa in cell cycle progression and development. Plant Physiol. 2006;140(4):1355–66. https://doi.org/10.1104/pp.106.077990.

  15. Bénédicte D, Elena RP, Qi X, Chua NH, Crisanto G. Cell type-specific role of the retinoblastoma/E2F pathway during Arabidopsis leaf development. Plant Physiol. 2005;140(1):67–80. https://doi.org/10.1104/pp.105.071027.

    Article  CAS  Google Scholar 

  16. Erika Ő, Csaba P, Binish M, Aladár P, Tünde L. E2FB interacts with RETINOBLASTOMA RELATED and regulates cell proliferation during leaf development. Plant Physiol. 2020;182(1):518–33. https://doi.org/10.1104/pp.19.00212.

    Article  CAS  Google Scholar 

  17. Horvath BM, Kourova H, Nagy S, Nemeth E, Magyar Z, Papdi C, Ahmad Z, Sanchez-Perez GF, Perilli S, Blilou I. Arabidopsis RETINOBLASTOMA RELATED directly regulates DNA damage responses through functions beyond cell cycle control. EMBO J. 2017;36(9):1261–78. https://doi.org/10.15252/embj.201694561.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Li D, Su Z, Dong J, Wang T. An expression database for roots of the model legume Medicago Truncatula under salt stress. BMC Genomics. 2009;10(1):517. https://doi.org/10.1186/1471-2164-10-517.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ma TY, Li ZW, Zhang SY, Liang GT, Guo J. Identification and expression analysis of the E2F/DP genes under salt stress in Medicago truncatula. Genes Genom. 2014;36:819–28. https://doi.org/10.1007/s13258-014-0218-5.

    Article  CAS  Google Scholar 

  20. Yu Y, Wang XL, Zhang YR, Song TQ, Zhang SX, Zhang XK, Chen DS. Genome-wide identification and characterization of the E2F/DP Transcription Factor Family in Triticum aestivum L. Russ J Plant Physiol. 2022;69(2):21. https://doi.org/10.1134/S1021443722020224.

    Article  CAS  Google Scholar 

  21. Liu W, Liu K, Chen D, Zhang Z, Li B, El-Mogy MM, Tian S, Chen T. Solanum lycopersicum, a Model Plant for the studies in Developmental Biology. Stress Biol Food Sci Foods. 2022;11(16):2402. https://doi.org/10.3390/foods11162402.

    Article  CAS  Google Scholar 

  22. Food and Agriculture Organization of the United Nations. 2021. Available online: http://www.fao.org/home/en/.

  23. Bai Y, Kissoudis C, Yan Z, Visser RGF, van der Linden G. Plant behaviour under combined stress: tomato responses to combined salinity and pathogen stress. Plant J. 2018;93(4):781–93. https://doi.org/10.1111/tpj.13800.

    Article  CAS  PubMed  Google Scholar 

  24. Krishna R, Karkute SG, Ansari WA, Jaiswal DK, Verma JP, Singh M. Transgenic tomatoes for abiotic stress tolerance: status and way ahead. 3 Biotech. 2019;9(4):143. https://doi.org/10.1007/s13205-019-1665-0.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Costa JM, Heuvelink EP. The global tomato industry. In Tomatoes 2018;1–26. Wallingford UK: CABI.

  26. Sato S, Tabata S, Hirakawa H, Asamizu E, Shirasawa K, Isobe S, Kaneko T, Nakamura Y, Shibata D, Aoki K, Egholm M. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485(7400):635–41. https://doi.org/10.1038/nature11119.

    Article  CAS  Google Scholar 

  27. Safder I, Shao G, Sheng Z, Hu P, Tang S. Genome-wide identification studies–A primer to explore new genes in plant species. Plant Biol. 2022;24:9–22. https://doi.org/10.1111/plb.13340.

    Article  CAS  PubMed  Google Scholar 

  28. Verma V, Ravindran P, Kumar PP. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016;16:86. https://doi.org/10.1186/s12870-016-0771-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Filipowicz W, Bhattacharyya SN, Sonenberg N. Mechanisms of post-transcriptional regulation by microRNAs: are the answers insight? Nat Rev Genet. 2008;9:102–14. https://doi.org/10.1038/nrg2290.

    Article  CAS  PubMed  Google Scholar 

  30. López-Galiano MJ, Sentandreu V, Martínez‐Ramírez AC, Rausell C, Real MD, Camañes G, Ruiz‐Rivero O, Crespo‐Salvador O, García‐Robles I. Identification of stress associated microRNAs in Solanum lycopersicum by high‐throughput sequencing. Genes. 2019;10(6):475. https://doi.org/10.3390/genes10060475.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER suite: protein structure and function prediction. Nat Methods. 2015;12:7–8. https://doi.org/10.1038/nmeth.3213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Li L, Shi Q, Li Z, Gao J. Genome-wide identification and functional characterization of the PheE2F/DP gene family in Moso bamboo. BMC Plant Biol. 2021;21(1):1–5. https://doi.org/10.1186/s12870-021-02924-8.

    Article  CAS  Google Scholar 

  33. Wolfe KH, Gouy M, Yang YW, Sharp PM, Li WH. Date of the monocot-dicot divergence estimated from chloroplast DNA sequence data. Proc Natl Acad Sci USA. 1989;86:6201–5. https://doi.org/10.1073/pnas.86.16.6201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Zhu M, Chen G, Zhang J, Zhang Y, Xie Q, Zhao Z, Pan Y, Hu Z. The abiotic stress-responsive NAC-type transcription factor SlNAC4 regulates salt and drought tolerance and stress-related genes in tomato (Solanum lycopersicum). Plant Cell Rep. 2014;33(11):1851–63. https://doi.org/10.1007/s00299-014-1662-z.

    Article  CAS  PubMed  Google Scholar 

  35. Cannon SB, Mitra A, Baumgarten A, Young ND, May G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004;4:10. https://doi.org/10.1186/1471-2229-4-10.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Li Z, Jiang H, Zhou L, Deng L, Lin Y, Peng X, Yan H, Cheng B. Molecular evolution of the HD-ZIP I gene family in legume genomes. Gene. 2014;533:218–28. https://doi.org/10.1016/j.gene.2013.09.084.

    Article  CAS  PubMed  Google Scholar 

  37. Priest HD, Filichkin SA, Mockler TC. Cis-regulatory elements in plant cell signaling. Curr Opin Plant Biol. 2009;12:643–9. https://doi.org/10.1016/j.pbi.2009.07.016.

    Article  CAS  PubMed  Google Scholar 

  38. Biłas R, Szafran K, Hnatuszko-Konka K, Kononowicz AK. Cis-regulatory elements used to control gene expression in plants. Plant Cell Tissue Organ Cult. 2016;127:269–87. https://doi.org/10.1007/s11240-016-1057-7.

    Article  CAS  Google Scholar 

  39. Wai AH, Waseem M, Cho LH, Kim ST, Lee DJ, Kim CK, Chung MY. Comprehensive genome-wide analysis and expression pattern profiling of the SlHVA22 gene family unravels their likely involvement in the abiotic stress adaptation of tomato. Int J Mol Sci. 2022;23(20):12222. https://doi.org/10.3390/ijms232012222.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Wai AH, Rahman MM, Waseem M, Cho LH, Naing AH, Jeon JS, Lee DJ, Kim CK, Chung MY. Comprehensive genome-wide analysis and expression pattern profiling of PLATZ Gene Family members in Solanum Lycopersicum L. under multiple abiotic stresses. Plants. 2022;11(22):3112. https://doi.org/10.3390/plants11223112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol. 2006;57:19–53. https://doi.org/10.1146/annurev.arplant.57.032905.105218.

    Article  CAS  PubMed  Google Scholar 

  42. Chung MY, Nath UK, Vrebalov J, Gapper N, Lee JM, Lee DJ, Kim CK, Giovannoni J. Ectopic expression of miRNA172 in tomato (Solanum lycopersicum) reveals novel function in fruit development through regulation of an AP2 transcription factor. BMC Plant Biol. 2020;20(1):283. https://doi.org/10.1186/s12870-020-02489-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Li W, Wang T, Zhang Y, Li Y. Overexpression of soybean miR172c confers tolerance to water deficit and salt stress, but increases ABA sensitivity in transgenic Arabidopsis thaliana. J Exp Bot. 2016;67(1):175–94. https://doi.org/10.1093/jxb/erv450.

    Article  CAS  PubMed  Google Scholar 

  44. Candar-Cakir B, Arican E, Zhang B. Small RNA and degradome deep sequencing reveals drought-and tissue-specific micrornas and their important roles in drought-sensitive and drought-tolerant tomato genotypes. Plant Biotechnol J. 2016;14(8):1727–46. https://doi.org/10.1111/pbi.12533.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Liu M, Yu H, Zhao G, Huang Q, Lu Y, Ouyang B. Identification of drought-responsive microRNAs in tomato using high-throughput sequencing. Funct Integr Genomics. 2018;18(1):67–78. https://doi.org/10.1007/s10142-017-0575-7.

    Article  CAS  PubMed  Google Scholar 

  46. Ardito F, Giuliani M, Perrone D, Troiano G, Lo Muzio L. The crucial role of protein phosphorylation in cell signaling and its use as targeted therapy. Int J Mol Med. 2017;40:271–80. https://doi.org/10.3892/ijmm.2017.3036.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Hashiguchi A, Komatsu S. Impact of post-translational modifications of crop proteins under abiotic stress. Proteomes. 2016;4:42. https://doi.org/10.3390/proteomes4040042.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Hsu CC, Zhu Y, Arrington JV, Paez JS, Wang P, Zhu P, Chen IH, Zhu JK, Tao WA. Universal plant phosphoproteomics workflow and its application to tomato signaling in response to cold stress. Mol Cell Proteom. 2018;17:2068–80. https://doi.org/10.1074/mcp.TIR118.000702.

    Article  CAS  Google Scholar 

  49. Vu LD, Xu X, Zhu T, Pan L, van Zanten M, de Jong D, Wang Y, Vanremoortele T, Locke AM, Van De Cotte B, De Winne N. The membrane-localized protein kinase MAP4K4/TOT3 regulates thermomorphogenesis. Nat Commun. 2021;12(1):2842. https://doi.org/10.1038/s41467-021-23112-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Xu WH, Wang YS, Liu GZ, Chen X, Tinjuangjun P, Pi LY, Song WY. The autophosphorylated Ser686, Thr688, and Ser689 residues in the intracellular juxtamembrane domain of XA21 are implicated in stability control of rice receptor-like kinase. Plant J. 2006;45:740–51. https://doi.org/10.1111/j.1365-313X.2005.02638.x.

    Article  CAS  PubMed  Google Scholar 

  51. Robatzek S, Chinchilla D, Boller T. Ligand-induced endocytosis of the pattern recognition receptor FLS2 in Arabidopsis. Genes Dev. 2006;20:537–42. https://doi.org/10.1101/gad.366506. http://www.genesdev.org/cgi/doi/.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Zhang X, Zhuang L, Liu Y, Yang Z, Huang B. Protein phosphorylation associated with drought priming-enhanced heat tolerance in a temperate grass species. Hortic Res. 2020;7:207. https://doi.org/10.1038/s41438-020-00440-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Pu X, Yang L, Liu L, Dong X, Chen S, Chen Z, Liu G, Jia Y, Yuan W, Liu L. Genome-wide analysis of the MYB transcription factor superfamily in Physcomitrella patens. Int J Mol Sci. 2020;21(3):p975. https://doi.org/10.3390/ijms21030975.

    Article  CAS  Google Scholar 

  54. Sampieri L, Di Giusto P, Alvarez C. CREB3 transcription factors: ER-golgi stress transducers as hubs for cellular homeostasis. Front Cell Dev Biol. 2019;7:p123. https://doi.org/10.3389/fcell.2019.00123.

    Article  Google Scholar 

  55. Sánchez-Camargo VA, Suárez-Espinoza C, Romero-Rodríguez S, Garza-Aguilar SM, Stam M, García-Ramírez E, Lara-Núñez A, Vázquez-Ramos JM. Maize E2F transcription factors. Expression, association to promoters of S-phase genes and interaction with the RBR1 protein in chromatin during seed germination. Plant Sci. 2020;296:110491. https://doi.org/10.1016/j.plantsci.2020.110491.

    Article  CAS  PubMed  Google Scholar 

  56. Liu S, Zhang Y, Feng Q, Qin L, Pan C, Lamin-Samu AT, Lu G. Tomato AUXIN RESPONSE FACTOR 5 regulates fruit set and development via the mediation of auxin and gibberellin signaling. Sci Rep. 2018;8(1):2971. https://doi.org/10.1038/s41598-018-21315-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Qin F, Shinozaki K, Yamaguchi-Shinozaki K. Achievements and challenges in understanding plant abiotic stress responses and tolerance. Plant Cell Physiol. 2011;52:1569–82. https://doi.org/10.1093/pcp/pcr106.

    Article  CAS  PubMed  Google Scholar 

  58. Huang YC, Niu CY, Yang CR, Jinn TL. The heat stress factor HSFA6b connects ABA signaling and ABA-mediated heat responses. Plant Physiol. 2016;172:1182–99. https://doi.org/10.1104/pp.16.00860.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Suzuki N, Bassil E, Hamilton JS, Inupakutika MA, Zandalinas SI, Tripathy D, Luo Y, Dion E, Fukui G, Kumazaki A, Nakano R. ABA is required for plant acclimation to a combination of salt and heat stress. PLoS ONE. 2016;11:e147625. https://doi.org/10.1371/journal.pone.0147625.

    Article  CAS  Google Scholar 

  60. Li W, Liu Y, Wang J, He M, Zhou X, Yang C, Yuan C, Wang J, Chern M, Yin J, Chen W, Ma B, Wang Y, Qin P, Li S, Ronald P, Chen X. The durably resistant rice cultivar Digu activates defence gene expression before the full maturation of Magnaporthe Oryzae appressorium. Mol Plant Pathol. 2016;17(3):354–68. https://doi.org/10.1111/mpp.12286.

    Article  CAS  PubMed  Google Scholar 

  61. Heidari P, Reza Amerian M, Barcaccia G. Hormone profiles and antioxidant activity of cultivated and wild tomato seedlings under low-temperature stress. Agron. 2021;11(6):1146. https://doi.org/10.3390/agronomy11061146.

    Article  CAS  Google Scholar 

  62. Hashemipetroudi SH, Arab M, Heidari P, Kuhlmann M. Genome-wide analysis of the laccase (LAC) gene family in Aeluropus littoralis: a focus on identification, evolution and expression patterns in response to abiotic stresses and ABA treatment. Front Plant Sci. 2023;14:1112354. https://doi.org/10.3389/fpls.2023.1112354/full. https://www.frontiersin.org/articles/.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Ahmad MZ, Shah Z, Ullah A, Ahmed S, Ahmad B, Khan A. Genome wide and evolutionary analysis of heat shock protein 70 proteins in tomato and their role in response to heat and drought stress. Mol Biol Rep. 2022;49(12):11229–41. https://doi.org/10.1007/s11033-022-07734-1.

    Article  CAS  PubMed  Google Scholar 

  64. Kandeel M, Morsy MA, Abd El-Lateef HM, Marzok M, El-Beltagi HS, Al Khodair KM, Albokhadaim I, Venugopala KN. Genome-wide identification of B3 DNA-Binding Superfamily Members (ABI, HIS, ARF, RVL, REM) and their involvement in stress responses and development in Camelina sativa. Agronomy. 2023;13(3):648. https://doi.org/10.3390/agronomy13030648.

    Article  CAS  Google Scholar 

  65. Waller JC, Akhtar TA, Lara-Núñez A, Gregory JF 3rd, McQuinn RP, Giovannoni JJ, Hanson AD. Developmental and feedforward control of the expression of folate biosynthesis genes in tomato fruit. Mol Plant. 2010;3(1):66–77. https://doi.org/10.1093/mp/ssp057.

    Article  CAS  PubMed  Google Scholar 

  66. Ding H, Mo S, Qian Y, Yuan G, Wu X, Ge C. Integrated proteome and transcriptome analyses revealed key factors involved in tomato (Solanum lycopersicum) under high temperature stress. Food Energy Secur. 2020;9(4):pe239. https://doi.org/10.1002/fes3.239.

    Article  Google Scholar 

  67. Lashbrooke J, Adato A, Lotan O, Alkan N, Tsimbalist T, Rechav K, Fernandez-Moreno JP, Widemann E, Grausem B, Pinot F, Granell A, Costa F, Aharoni A. The Tomato MIXTA-Like transcription factor coordinates Fruit Epidermis Conical Cell Development and Cuticular lipid biosynthesis and assembly. Plant Physiol. 2015;169(4):2553–71. https://doi.org/10.1104/pp.15.01145.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Zhao T, Wu T, Zhang J, Wang Z, Pei T, Yang H, Li J, Xu X. Genome-wide analyses of the genetic screening of C2H2-Type zinc finger transcription factors and Abiotic and biotic stress responses in Tomato (Solanum lycopersicum) based on RNA-Seq Data. Front Genet. 2020;11:540. https://doi.org/10.3389/fgene.2020.00540.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Olivieri F, Graci S, Francesca S, Rigano MM, Barone A. Accelerating the development of heat tolerant tomato hybrids through a multi-traits evaluation of parental lines combining phenotypic and genotypic analysis. Plants. 2021;10(10):2168. https://doi.org/10.3390/plants10102168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Zhang T, Wang X, Lu Y, Cai X, Ye Z, Zhang J. Genome-wide analysis of the cyclin gene family in tomato. Inter J Mol Sci. 2013;15(1):120–40. https://doi.org/10.3390/ijms15010120.

    Article  CAS  Google Scholar 

  71. Tian F, Yang DC, Meng YQ, Jin JP, Gao G. PlantRegMap: charting functional regulatory maps in plants. Nucleic Acids Res. 2019;48(1):1104–13. https://doi.org/10.1093/nar/gkz1020.

    Article  CAS  Google Scholar 

  72. Mueller LA, Solow TH, Taylor N, Skwarecki B, Buels R, Binns J, Lin C, Wright MH, Ahrens R, Wang Y, Herbst EV. The, SOL Genomics, Network. A comparative resource for Solanaceae biology and beyond. Plant Physiol. 2005;138:1310–7. https://doi.org/10.1104/pp.105.060707.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. GasteigerE,Hoogland CGA, Wilkins MR, Appel RD, Bairoch A. Protein identification and analysis tools on the ExPASy server. The Proteomics protocols Handbook. Totova, NJ, USA: Humana; 2005. pp. 571–607. https://doi.org/10.1385/1-59259-890-0:571.

    Chapter  Google Scholar 

  74. Guo AY, Zhu QH, Chen X, Luo JC. GSDS: a gene structure display server. Yi Chuan. 2007;29(8):1023–6. https://doi.org/10.1360/yc-007-1023.

    Article  CAS  PubMed  Google Scholar 

  75. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:202–8. https://doi.org/10.1093/nar/gkp335.

    Article  CAS  Google Scholar 

  76. Robert X, Gouet P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014;42:320–4. https://doi.org/10.1093/nar/gku316.

    Article  CAS  Google Scholar 

  77. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thomposon JD. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal, Omega. Mol Syst Biol. 2011;7(1):539. https://doi.org/10.1038/msb.2011.75.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K. WoLF PSORT: protein localization predictor. Nucleic Acids Res. 2007;35:585–7. https://doi.org/10.1093/nar/gkm259.

    Article  Google Scholar 

  79. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38:3022–7. https://doi.org/10.1093/molbev/msab120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, Xia R. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.

    Article  CAS  PubMed  Google Scholar 

  81. Koch MA, Haubold B, Mitchell-Olds T. Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol Biol Evol. 2000;17(10):1483–98. https://doi.org/10.1093/oxfordjournals.molbev.a026248.

    Article  CAS  PubMed  Google Scholar 

  82. Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouzé P, Rombauts S. PlantCARE., a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30(1):325–7. https://doi.org/10.1093/nar/30.1.325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Blom N, Gammeltoft S, Brunak S. Sequence- and structure‐based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999;294(5):1351–62. https://doi.org/10.1006/jmbi.1999.3310.

    Article  CAS  PubMed  Google Scholar 

  84. Xu D, Zhang Y. Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization. Biophys J. 2011;101:2525–34. https://doi.org/10.1016/j.bpj.2011.10.024.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Wai AH, Cho LH, Peng X, Waseem M, Lee DJ, Lee JM, Kim CK, Chung MY. Genome-wide identification and expression profiling of Alba gene family members in response to abiotic stress in tomato (Solanum lycopersicum L). BMC Plant Biol. 2021;21:530. https://doi.org/10.1186/s12870-021-03310-0.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Balestrini R, Gómez-Ariza J, Lanfranco L, Bonfante P. Laser microdissection reveals that transcripts for five plant and one fungal phosphate transporter genes are contemporaneously present in arbusculated cells. Mol Plant-Microbe Interact. 2007;20:1055–62. https://doi.org/10.1094/MPMI-20-9-1055.

    Article  CAS  PubMed  Google Scholar 

  87. Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8. https://doi.org/10.1038/nprot.2008.73.

    Article  CAS  PubMed  Google Scholar 

  88. Andrews S. FastQC A quality control tool for high throughput sequence data. In Babraham Bioinformatics; Babraham Institute: Cambridge, UK, 2010.

  89. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019;37:907–15. https://doi.org/10.1038/s41587-019-0201-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Vandel J, Gheeraert C, Staels B, Eeckhoute J, Lefebvre P, Dubois-Chevalier J. GIANT: galaxy-based tool for interactive analysis of transcriptomic data. Sci Rep. 2020;10:19835. https://doi.org/10.1038/s41598-020-76769-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Liao Y, Smyth GK, Shi W. FeatureCounts: an efficient general-purpose read summarization program. Bioinformatics. 2014;30:923–30.

    Article  CAS  PubMed  Google Scholar 

  92. Love MI, Huber W, Anders S. Moderated estimation of Fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:1–21. https://doi.org/10.1186/s13059-014-0550-8.

    Article  CAS  Google Scholar 

  93. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:559. https://doi.org/10.1186/1471-2105-9-559.

    Article  CAS  Google Scholar 

  94. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30. https://doi.org/10.1093/nar/28.1.27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  95. Kim SR, Lee DY, Yang JI, Moon S, An G. Cloning vectors for rice. J Plant Biol. 2009;52:73–8. https://doi.org/10.1007/s12374-008-9008-4.

    Article  CAS  Google Scholar 

  96. Chang S, Chen Y, Jia S, Li Y, Liu K, Lin Z, Wang H, Chu Z, Liu J, Xi C, Zhao H. Auxin apical dominance governed by the OsAsp1-OsTIF1 complex determines distinctive rice caryopses development on different branches. PLoS Genet. 2020;16(10):e1009157. https://doi.org/10.1371/journal.pgen.1009157.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  97. Page MT, Parry MA, Carmo-Silva E. A high‐throughput transient expression system for rice. 2019; 42(7): 2057:2064. https://doi.org/10.1111/pce.13542.

Download references

Acknowledgements

Not applicable.

Funding

Not Applicable.

Author information

Authors and Affiliations

Authors

Contributions

M.-Y.C., D.-j.L. and C-K.K. supervised and designed the experiments. D.D., L.-H.C and D.K. carried out the experiments. D.D. performed the Insilico analysis and data analysis. D.D. and M.-Y.C. prepared the manuscript. A.H.K.R. and C.-K.K. critically edited the manuscript. All authors made contribution to the manuscript. All the authors read and approved the final manuscript.M.-Y.C., and C-K.K. are equally corresponding authors.

Corresponding authors

Correspondence to Chang-Kil Kim or Mi-Young Chung.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Divya, D., Robin, A.H.K., Cho, LH. et al. Genome-wide characterization and expression profiling of E2F/DP gene family members in response to abiotic stress in tomato (Solanum lycopersicum L.). BMC Plant Biol 24, 436 (2024). https://doi.org/10.1186/s12870-024-05107-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05107-3

Keywords