Next Article in Journal
Molecular Characteristics of Cell Pyroptosis and Its Inhibitors: A Review of Activation, Regulation, and Inhibitors
Next Article in Special Issue
Alternative Splicing of TaHsfA2-7 Is Involved in the Improvement of Thermotolerance in Wheat
Previous Article in Journal
Differential Upregulation and Functional Activity of S1PR1 in Human Peripheral Blood Basophils of Atopic Patients
Previous Article in Special Issue
Comparative Ubiquitination Proteomics Revealed the Salt Tolerance Mechanism in Sugar Beet Monomeric Additional Line M14
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Transcriptional Networks Involved in De Novo Organ Formation in Tomato Hypocotyl Explants

by
Eduardo Larriba
,
Míriam Nicolás-Albujer
,
Ana Belén Sánchez-García
and
José Manuel Pérez-Pérez
*
Instituto de Bioingeniería, Universidad Miguel Hernández, 03202 Elche, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(24), 16112; https://doi.org/10.3390/ijms232416112
Submission received: 11 November 2022 / Revised: 9 December 2022 / Accepted: 12 December 2022 / Published: 17 December 2022
(This article belongs to the Special Issue Omics Study to Uncover Signalling and Gene Regulation in Plants)

Abstract

:
Some of the hormone crosstalk and transcription factors (TFs) involved in wound-induced organ regeneration have been extensively studied in the model plant Arabidopsis thaliana. In previous work, we established Solanum lycopersicum “Micro-Tom” explants without the addition of exogenous hormones as a model to investigate wound-induced de novo organ formation. The current working model indicates that cell reprogramming and founder cell activation requires spatial and temporal regulation of auxin-to-cytokinin (CK) gradients in the apical and basal regions of the hypocotyl combined with extensive metabolic reprogramming of some cells in the apical region. In this work, we extended our transcriptomic analysis to identify some of the gene regulatory networks involved in wound-induced organ regeneration in tomato. Our results highlight a functional conservation of key TF modules whose function is conserved during de novo organ formation in plants, which will serve as a valuable resource for future studies.

Graphical Abstract

1. Introduction

Although plants have an outstanding ability to repair damaged tissues of even regenerate lost organs, our current understanding of the molecular networks involved in their regulation [1,2,3] is mainly limited to studies using the Arabidopsis thaliana model system [4]. The first gene regulatory network (GRN) associated with cell reprogramming during plant regeneration was obtained by using an enhanced yeast one-hybrid screen with 252 TFs and 48 gene promoters likely to be relevant for regeneration in A. thaliana [5]. Their results revealed several regulatory nodes, including TFs from the APETALA 2/ETHYLENE RESPONSE FACTOR (AP2/ERF) and LATERAL ORGAN BOUNDARY/ASYMMETRIC LEAVES 2 DOMAIN (LBD) families, which might have overlapping targets associated with responses to stress, cytokinin, auxin and cell proliferation [5]. More recently, a web-based application and analytical platform, REGENOMICS, was implemented and contains all currently available regeneration-specific RNA-sequencing (RNA-Seq) datasets. REGENOMICS allows users to perform integrative analyses to generate co-expression networks and provide functional information about plant regeneration processes [6]. However, the current version of REGENOMICS is composed primarily of A. thaliana RNA-Seq datasets, making it difficult to find conserved factors involved in plant regeneration through multi-species analyses.
Several studies have identified various wound-induced TFs that have critical roles in regeneration [1,2,3]. One of these TFs belonging to the AP2/ERF family is WOUND INDUCED DEDIFFERENTIATION 1 (WIND1) and its closest homologs WIND2–WIND4. WIND TFs promote callus formation by activating CK response [7], and their ectopic expression can bypass both wounding and auxin pre-treatment, as well as increase de novo shoot regeneration from root explants of different plant species, including tomato [8,9]. WIND1 promotes shoot regeneration through the direct activation of the AP2/ERF TF ENHANCER OF SHOOT REGENERATION 1 (ESR1) [10]. Recent results from time-course transcriptome analysis upon WIND1 induction have revealed that WIND TFs function as key regulators of wound-induced responses by promoting dynamic transcriptional alterations of genes involved in cell reprogramming, vascular regeneration and defense responses [11].
In this study, we carried out a bioinformatics approach to characterize the expression of all TF genes during de novo organ formation in Solanum lycopersicum “Micro-Tom” using the well-established model of young hypocotyl explants developed in our laboratory [12,13,14]. We have annotated 1399 expressed genes putatively encoding TFs and studied their expression profiles using k-means [15] and STEM [16] analysis in the apical and basal regions of the hypocotyl and over time. Because known TF modules required for organ regeneration are conserved between distant plant species, we have established homologous GRNs for some of the TF genes that showed transient upregulation during cell reprogramming and/or tissue specification. New experimental data on cell cycle regulation and reactive oxygen species (ROS) distribution are also provided. Our results have revealed candidate TFs to regulate tissue regeneration that could serve to bridge the gap in our knowledge about tissue-specific wound-induced organ formation in tomato.

2. Results

2.1. Expression of Genes Encoding TFs during De Novo Organ Formation

From a previous RNA-Seq experiment [13], we identified 1399 expressed genes putatively encoding TFs that were assigned to 48 families (Table S1). We found that 125 of these genes (8.93%) were preferentially expressed in the apical region of hypocotyl explants between 0 and 192 h after excision (hae; Figure 1a), and 33 of them were specifically expressed in the apical region over time. We constructed a TF interaction network based on the results of STRING analyses using data from its A. thaliana orthologs and found significantly enriched gene ontologies (GOs) related to shoot meristem formation and leaf development (Figure 1b). Tomato orthologs of WUSCHEL (WUS; Solyc02g083950), AGAMOUS (AG; Solyc02g071730) and CUP-SHAPED COTYLEDON 2 (CUC2; Solyc07g062840), among others, were upregulated at specific time points in the apical region (Figure 1c and Figure S1a). On the other hand, 271 TF genes (19.37%) were preferentially expressed in the basal region of hypocotyl explants over time, and 101 of these genes were exclusively expressed in the basal region (Figure 1d and Figure S1b). Through STRING analysis, we identified genes belonging to the TF families of LBD and PLETHORA/AINTEGUMENTA-LIKE (PLT/AIL); some of its members are well-known regulators of the specification and growth of the root meristem (Figure 1e). Tomato orthologs of WUSCHEL-RELATED HOMEOBOX 5 (WOX5; Solyc03g096300), PLT1/2 (Solyc11g061750), PUCHI (Solyc01g067540), FEZ (Solyc05g007550) and SOMBRERO (SMB; Solyc12g017400) were upregulated during adventitious root (AR) formation in the basal region of hypocotyl explants. In addition, the expression of the superlocus SHOOTBORNE ROOTLESS (SBRL; Solyc09g066260 and Solyc09g066270), which controls the initiation of ARs and lateral roots (LRs) in tomato [17], showed higher expression in the basal region of hypocotyl explants after wounding (Figure 1f). These results confirm the functional relevance of our transcriptome dataset for identifying the transcriptional networks involved in de novo organ formation in tomato.

2.2. Analysis of TF Expression Profiles during De Novo Organ Formation

The 1296 TF genes that were expressed in the apical region of the hypocotyl were grouped into eight expression clusters. Clusters C1 to C4 included TF genes whose expression decreased over time, whereas TF gene expression in clusters C5 to C8 increased relative to excision time (Figure 2a). To investigate the expression profiles of TF genes, we performed a time-series analysis using STEM [16]. Ten different gene expression profiles were chosen for further study. Most of the constitutively downregulated TF genes were found in P0, P2, and P7 expression profiles, while P3 included specifically downregulated TF genes at 24–96 hae (Figure 2b,c). On the other hand, constitutively upregulated TFs were distributed across P13, P14 and P15 expression profiles (Figure 2b,d). Most of the upregulated TF genes at 24 hae were included in the P11 profile, whereas the P5 and P6 profiles included upregulated TF genes at the final time point of 192 hae (Figure 2b,e). Next, we selected 111 TF genes, whose orthologs in A. thaliana have previously been linked to regeneration (Table S2), to use them as markers for the characterization of the studied expression profiles (Figure 2c–e). In A. thaliana, WIND1 and ESR1 establish a central regulatory hub necessary for shoot regeneration [10]. Consistent with their role in shoot regeneration, tomato orthologs of WIND1 (Solyc04g072900) and ESR1/2 (Solyc05g013540) were upregulated at 24 hae in the apical region of the hypocotyl and were recovered from P7 and P14 profiles, respectively (Figure 2c,d). Tomato orthologs of know genes required for root development, such as SCARECROW (SCR; Solyc10g074680) or PLT5 (Solyc07g018290), were found to be downregulated in the apical region of hypocotyl explants during de novo shoot formation and were included in the P0 profile (Figure 2c,f).
Similar to what was found in the apical region, k-means analysis of the 1364 TF genes expressed in the basal region of the hypocotyl also yielded eight expression clusters (Figure 3a). Most of the downregulated TFs in clusters C1 to C3 were assigned by STEM analysis to gene expression profiles P0, P2, P3, and P7 (Figure 3b,c). In contrast, the upregulated TF genes in the basal region of the hypocotyl during AR formation (clusters C6 to C8) were assigned to P13, P14, and P15 expression profiles (Figure 3b,d). TF genes included in the P5 and P6 profiles also displayed dynamic expression patterns, which were characterized by downregulation at earlier time points and upregulation at later time points (Figure 3b,e). Most of these TF genes were found in the C4, C7 and C8 k-means clusters (Figure 3b). The C5 cluster contained upregulated TF genes at 24 hae, and almost all were included in the P11 profile (Figure 3b,e). Some previously known TF genes involved in cell reprogramming and plant regeneration (see above) were assigned to the studied STEM profiles (Figure 3c–e). Many know regulators of shoot regeneration, such as WIND1 and ESR1, were downregulated in the basal region of the hypocotyl (Figure 3c), in opposite contrast to that found for these genes in the apical region (see above). On the other hand, the expression of putative tomato orthologs of some key regulators of root specification and growth, such as PLT/AIL genes (PLT1/2, PLT3 (Solyc11g010710), and BABY BOOM (BBM; Solyc11g008560)), MONOPTEROS (MP; Solyc04g081240) and WOX8/9 (Solyc02g077390), among others, were upregulated over time specifically in the basal region of the hypocotyl (Figure 3f). PLT5 and SHORT-ROOT (SHR; Solyc02g092370) expression decreased from 24 to 96 h, and their steady-state levels recovered at 192 h, coinciding with AR emergence (Figure 3c,f). Among the genes deserving of further study is the putative tomato ortholog of WOX13, which was found in the P11 expression profile (Figure 3e,f). Taken together, our results have enabled us to identify putative TFs that might be involved in de novo AR and shoot formation in tomato hypocotyl explants after wounding.

2.3. GRNs Required for De Novo Organ Formation

In our previous work [13,14], we identified some GO categories that were specifically enriched in the apical or basal regions of the hypocotyl and at different times during de novo organ formation. To investigate the regulatory networks of TFs and their targets that contribute to the observed differences between de novo AR and shoot formation, we focused on specific STEM expression profiles. We searched for putative TF gene targets in each profile and performed GO ontology enrichment analyses between them (see Section 4). We reasoned that TF genes upregulated at 24 hae, most of which are included in the P11 profile in both regions (Figure 2e and Figure 3e), could be important for cell reprogramming and tissue specification during de novo organ formation.
In the apical region, we generated GRNs that included TFs of the P11 and P13 profiles and their putative target genes that were differentially expressed over time (see Section 4). Among the upregulated targets of TFs expressed in the P11 profile, we found GO enrichment for terms related to cell cycle regulation and glycolate/glyoxylate metabolism (i.e., photorespiration), while the downregulated targets associated with TF genes in this profile were enriched for GO terms related to the light-harvesting complex of photosynthesis and auxin responses (Figure S2a,c and Table S3). We found that putative tomato orthologs of ERF3 (Solyc07g064890), BASIC LEUCINE-ZIPPER 44 (bZIP44; Solyc01g109880), and ELONGATED HYPOCOTYL 5 (HY5; Solyc08g061130) connected different sets of core cell cycle regulators (Figure 4a,c). Most of the cell cycle targets of bZIP44 and HY5 overlapped (Figure 4a and Table S4), suggesting a cooperative function of these two TFs. In contrast, the Solyc03g116100 node, which encodes the putative tomato ortholog of MYB DOMAIN PROTEIN 94 (MYB94) and MYB96, connected some genes in the auxin pathway that were downregulated over time in the apical region (Figure 4a,d).
The P13 profile of constitutively upregulated TF genes in the apical region included TFs involved in shoot apical meristem specification, such as WUS (Figure 2d). We found that GO terms related to cell cycle regulation and oxidoreductase activity were specifically enriched among upregulated targets of TF genes assigned to this profile (Figure S2b,d and Table S3). One of the most prevalent cell cycle-related TF nodes in this profile was that of Solyc03g113760 (Figure 4b), which putatively encodes the DP-E2F-LIKE 1 (DEL1) protein known to inhibit the endocycle in A. thaliana [18]. Two other TF nodes, Solyc01g095460 and Solyc06g065820 (Figure 4e), encoding the putative orthologs of G-BOX BINDING FACTOR 3 and ETHYLENE AND SALT INDUCIBLE 3, included genes that might be involved in ROS regulation (Figure 4b,f). We found 25 putative targets involved in cell cycle regulation that were included in the GRNs of P11 and P13 profiles, most of which are involved in DNA replication and showed upregulation at 24 and 192 hae (Table S4).
In the basal region, most of the upregulated targets of TF genes within P11 expression profile (Figure 3e) were related to ribosome function, whereas the downregulated target genes were enriched for photosynthesis-related and oxidoreductase activity GO terms (Table S4). Of the many GRNs entangled in this profile (Figure S3a,c), we focused on those of Solyc01g067540 and Solyc05g054400, putatively encoding orthologs of PUCHI [19] and GATA15/17 [20], respectively, in A. thaliana. The node of GATA15/17 included genes encoding glutaredoxins that were deregulated over time, whereas the node of PUCHI included several genes encoding constituent proteins of the cytoplasmic ribosome (Figure 5a and Table S4). Some of the key regulators of root founder cell specification, such as SBRL or LBD16 [21] were assigned to the P15 profile (Figure 3d). From this profile, a complex GRN involving TF genes and their putative targets was constructed (Figure S3b,d and Table S4). The TF gene with the highest number of interactions with its putative deregulated targets was Solyc11g008560 (Figure 5b), which encodes the BBM TF that is required in cooperation with other PLT genes to maintain stem cell activity in A. thaliana roots [22]. Among the putative BBM targets identified, several cell cycle-related genes encoding cyclins (e.g., Solyc11g005090 (CycA1;1), Solyc02g082820 (CycB2;1/CycB2;2) and Solyc02g092980 (CycD3;3/CycD3;2)), were found to be upregulated after 24 h in the basal region of the hypocotyl (Figure 5c). Some BBM targets involved in cell wall biogenesis (e.g., cellulose synthases and pectin methylesterases) were also similarly regulated (Figure 5c). On the other hand, putative BBM targets involved in the transduction of various hormones, such as auxins (SlARF1, SlIAA14, SlIAA17 and SlIAA7) or abscisic acid, were strongly downregulated starting at 24 h (Figure 5d). We found several BBM targets encoding ribosomal proteins, which do not overlap with PUCHI targets encoding other ribosomal proteins (Table S4).

2.4. The Cell Cycle Is Differentially Regulated in the Apical and Basal Regions during de Novo Organ Formation

Since the results found above suggested a differential regulation of the cell cycle between the apical and basal regions of the hypocotyl, we measured the DNA content by flow cytometry (Figure 6a,b) and found significant differences in the ploidy profiles of the apical and the basal region of the hypocotyl over time (Figure 6c and Figure S4a,b). We estimated the mitotic index (MI) and endoreduplication index (EI) as previously described [23]. In the apical region, MI increased significantly at 24 and 96 hae relative to 0 hae, and MI values were restored to steady-state levels at 192 hae (Figure 6d). In the basal region at 0 hae, MI values were higher than those in the apical region and did not change significantly (Figure 6d). The apical region at 0 hae showed the lowest EI values, which were similar in both regions and over time (Figure 6e).
We found 172 expressed genes that are associated with core cell cycle functions, 97 of which were deregulated over time (Table S5). We classified these genes into 16 functional groups and performed an enrichment analysis with respect to their expression at 0 hae (Table S5 and Figure S4c). In the apical region, we found an enrichment of upregulated DEGs related to DNA replication, DNA repair, and genes involved in G1/S transition or mitotic regulation, as well as downregulated DEGs involved in inhibition of endoreduplication (Figure 6f). On the other hand, upregulated DEGs related to cell proliferation in the basal region were found to be enriched at 96 hae (Figure 6g). Some genes encoding central regulators of mitosis, such as B-type CDKs (e.g., Solyc04g082840 and Solyc10g074720), genes encoding regulatory subunits of the anaphase promoting complex or cyclosome (APC/C) (e.g., Solyc03g096870 and Solyc06g043150), or genes associated with the mitotic spindle checkpoint (e.g., Solyc03g007060 and Solyc08g066050), were expressed at different time points in the apical and basal regions (Table S5 and Figure 6g). Most genes encoding A-type and B-type mitotic cyclins were upregulated in the apical and basal regions after wounding, but at different times between regions. The Solyc04g078310 and Solyc12g088530 genes, encoding the CycA3 cyclins whose A. thaliana orthologs are required for cell cycle transitions [24], were specifically upregulated at 24 hae in the apical region (Figure 6f). The CycA and CycB genes were strongly downregulated in the basal region at 24 hae and increased their expression at later times (Figure 6g). In turn, the expression of Solyc07g052610, which encodes an atypical cyclin [25], was differentially upregulated in the apical and basal regions of the hypocotyl explants after wounding (Table S5).

2.5. Characterization of ROS Accumulation during De Novo Organ Formation

Previous results suggested the involvement of ROS homeostasis in wound-induced regeneration in tomato hypocotyls [14]. To assess ROS production, we visualized H2O2 accumulation during wound-induced organ formation by staining with 3,3’-diaminobenzidine (DAB). At 24 hae, an increase in DAB staining was observed in the apical region of the explants (Figure 7a). The DAB precipitate was confined to tissues near the wounding site at 96 hae, in both the apical and basal regions of the explants (Figure 7a, 96 hae). At 192 hae, DAB staining decreased substantially and was restricted to more distal cells in the basal region and to specific cells within the apical callus (Figure 7a, T8). We confirmed the DAB-specific staining of H2O2 by incubation with catalase (Figure 7a, 24 hae cat).
Due to these specific patterns of ROS accumulation, we identified 110 DEGs related to ROS homeostasis (Figure 7b and Table S6) that were assigned to different families based on the enzymatic activities of their products [26] (Figure S5). In the apical region, the expression of most genes encoding ROS-scavenging proteins was upregulated (Figure S5a,b). Furthermore, the expression of five genes encoding membrane-dependent NADPH oxidases of the RBOH family (Figure S5a,b), which are involved in stress-induced singlet oxygen production [27,28], was upregulated in both the apical and basal regions (Figure S5b). Since an excess of singlet oxygen could cause increased ROS levels eventually leading to cell death [29], we investigated cell death by trypan blue staining (Figure 7c). At 24 hae, we observed faint trypan blue staining in the distal end of the apical region and no trypan blue staining in the basal region (Figure 7c, 24 hae). However, at 96 hae, trypan blue stained the outer cell layer of the basal region and the more distal cell layer of the apical region of the explants (Figure 7c, 96 hae). Similar trypan blue staining was observed in distal cells within the apical and basal regions of the hypocotyl explants at 192 hae (Figure 7c, 192 hae), indicating that the degree of ROS production was not associated directly with cell death. Genes encoding peroxiredoxins (PRX), thioredoxins (TRX), ascorbate peroxidases (APX), and catalases (CAT) were differentially regulated between the apical and basal regions (Figure S5a,b), which could explain the differences observed in the accumulation of H2O2 after wounding (Figure 7a).
We identified 85 expressed genes encoding class III peroxidases (POX), 81 of which (95.3%) were deregulated to a greater degree than genes encoding other ROS-scavenging proteins (Table S6). We grouped these genes according to six different expression profiles (Figure 7d–g). Forty-one DEGs were constitutively upregulated in both the apical region and basal regions over time (e.g., Prx54, Prx68), and some of these genes (e.g., Prx92, Prx97) were downregulated in the basal region at 24 hae (clusters A and B, respectively; Figure 7d). Twenty-nine DEGs showed contrasting deregulation in the apical and basal regions, and 21 of these DEGs were upregulated primarily in the basal region (e.g., Prx100, Prx105, clusters C and D; Figure 7e), and eight DEGs were upregulated only in the apical region (e.g., Prx44, Prx95, cluster E; Figure 7f). Eleven DEGs were mainly downregulated in both tissues (e.g., Prx20, Prx87, cluster F; Figure 7g). Our results indicate a differential regulation of POX-encoding genes in the apical and basal regions of hypocotyl explants over time, which could explain the dynamic ROS homeostasis observed during wound-induced organ formation in tomato hypocotyl explants.

2.6. CDF3 Is Required for De Novo Shoot Formation in Tomato Hypocotyl Explants

Our results indicate a differential regulation of the GRNs involved in de novo organ formation in the apical and basal region of the hypocotyl. To search for TFs putatively involved in tissue-specific cell reprogramming, we selected TFs with contrasting expression levels whose orthologs are known to regulate environmentally regulated developmental traits. One of these TFs was CYCLING DOF FACTOR 3 (CDF3), encoded by Solyc03g115940, which was identified in two contrasting expression profiles over time: the P6 profile in the apical region (Figure 2e) and the P14 profile in the basal region (Figure 3d). Using PlantRegMap and our RNA-Seq dataset (see Section 4), we identified 1063 deregulated genes (773 and 757 DEGs in the apical and basal regions of the hypocotyl, respectively; Table S7) among previously described putative CDF3 targets [14]. Our GO enrichment analysis suggests that CDF3 regulates auxin-related genes and other TF genes (Table S8). Some of the TF genes that might be downregulated by CDF3 in the apical region of hypocotyl explants are Solyc12g010800, the tomato ortholog of AGL15, and Solyc01g096070, which putatively encode the ortholog of ARF11 and ARF18 in A. thaliana (Table S7). Regarding the basal region of the hypocotyl, we found enrichment in GO terms related to inorganic anion transport (Table S8).
To characterize the relevance of this TF node during de novo organ formation, we studied two overexpression (OE) lines (#5.2 and #11.2) and one RNAi line (#1.1) of CDF3, which were compared with the Moneymaker (MM) background (see Section 4). Although the RNAi line showed a slight delay in germination compared to the other lines, we found no significant difference in the number of LRs after root tip excision with respect to their genotype (p value = 0.486; Figure S6a,b). After whole root excision, AR emergence from the hypocotyl occurred slightly earlier in the OE lines than in the RNAi line and MM (p value < 0.05; Figure 8a). However, the RNAi line showed a significant increase (p value < 0.01) in the number of ARs at 5 days after excision (dae) compared to the CDF3 OE lines and MM (Figure 8b). At 10 dae, we excised the shoot apex and studied de novo shoot regeneration as previously described [14]. The RNAi line showed a slight but significant delay (p value < 0.01) in de novo shoot regeneration compared to the OE lines (Figure 8c), which resulted in a smaller root system of the RNAi line at the end of the experiment (Figure 8d). These results suggested that CDF3 levels should be tightly regulated during de novo organ formation in a tissue-specific context, acting as a negative regulator of AR initiation but as a positive regulator of de novo shoot formation. Due to the known redundancy among the CDF genes [30], further studies will be needed to determine the relevance of this DNA binding with one finger (DOF) subfamily of TFs in wound-induced cell reprogramming in tomato.

3. Discussion

In previous work from our laboratory, we aimed to study de novo organ formation in tomato using “Micro-Tom” hypocotyl explants after wounding [12]. The current working model involves spatial and temporal regulation of auxin-to-CK gradients in the apical and basal regions of the hypocotyl soon after wounding to regulate cell reprogramming and founder cell activation [14]. We also found that wounding induced extensive metabolic reprogramming of some cells in the apical region of the hypocotyls, highlighting a key role of photorespiration in supplying the energy and structural elements necessary for initial callus growth preceding de novo shoot formation [13]. In this work, we extended our RNA-Seq analysis to identify some of the GRNs involved in wound-induced de novo organ formation in this species.
We found that the expression of the tomato WUS homologue (Solyc02g083950) in the apical region of the hypocotyl after wounding defined a GRN that could be involved in the specification of the new shoots. Previous studies in tomato have shown that after laser ablation of the central zone of the shoot apical meristem, WUS was ectopically induced in neighboring cells some distance from the lesion, which anticipated the establishment of new regenerated shoots [31]. Furthermore, Solyc02g071730 and Solyc05g009380, putatively encoding AG and JAGGED (JAG) TFs, showed a WUS-like expression profile, and we hypothesize that they could contribute to de novo shoot formation. A weak gain-of-function WUS allele in tomato, termed locule number, affected a CArG element downstream of WUS coding region that is capable of binding to AG and downregulating WUS expression that, in turn, terminates meristem activity during flower development [32]. JAG loss-of-function alleles in tomato (i.e., lyrate mutants) caused a broad reduction in overall leaf growth, but increased leaflet dissection by interacting with both the KNOTTED-LIKE HOMEOBOX (KNOX) genes and the auxin transcriptional network [33]. We hypothesize that in tomato, similar to what was found in A. thaliana [34], JAG could antagonize meristem identity genes, such as WUS, in cells that are switching from meristem to organ primordium identity, and repressed the KIP RELATED PROTEIN (KRP) genes in these cells, thus releasing the constraint for DNA replication, a hypothesis than can now be experimentally validated.
K-means and STEM analysis allowed us to identify TF genes with unique expression patterns during de novo organ formation in tomato. We found putative orthologs of two conserved TF modules in A. thaliana that have been described as having key roles in wound-induced cell reprogramming and cell cycle reentry during shoot regeneration. On the one hand, Solyc04g072900 and Solyc05g013540, which encode the tomato orthologs of WIND1 and ESR1/2, respectively, were upregulated at 24 hae in the apical region. WIND1 promotes callus formation and shoot regeneration by activating ESR1 expression [10]. Our results suggest that the function of the WIND1-ESR1/2 module during shoot regeneration might be conserved beyond the Brassicaceae, although other experiments with known loss-of-function tomato ESR1/2 mutants (i.e., leafless) [35] wait for confirmation. Other conserved module included Solyc07g042260 (AUXIN RESPONSE FACTOR 7 (ARF7)/ARF19) and Solyc11g071300 (MYB DOMAIN PROTEIN 3R-4 (MYB3R4)), which in A. thaliana has been recently shown to regulate cell cycle activation of leaf mesophyll protoplasts in response to exogenous auxin [36]. Implicit GRNs were constructed with the target genes from the upregulated TF genes in selected STEM profiles in the apical region (P11 and P13). One such upregulated node was that of Solyc03g116100, which encodes the putative tomato ortholog of TFs MYB94 and MYB96 connecting some downregulated genes in the auxin pathway (PIN-FORMED 1, ARF4, etc.) at earlier time points during de novo shoot formation. In A. thaliana, MYB94 and MYB96 repress LBD29 expression, thus playing negative regulatory roles in auxin-induced and callus formation mediated by ARF7 and ARF19 TFs [37]. Therefore, our results suggest a functional conservation of key TF modules involved in the regulation of wound-induced de novo shoot formation, which need to be confirmed experimentally using available genetic tools (known mutant alleles and CRISPR/Cas9-mediated gene editing). Functional analysis of CDF3 in A. thaliana suggests a direct role for this DOF TF in drought and osmotic stress responses, likely through regulation of cell integrity, carbon and nitrogen metabolism, and ROS homeostasis [38]. The ectopic expression of CDF3 in tomato increases the photosynthetic rate and biomass production under salinity stress by facilitating nitrogen assimilation and sucrose availability [30,39]. We found that CDF3 RNAi lines showed a slight delay in de novo shoot regeneration from tomato hypocotyl explants, suggesting that it might also play a role during cell reprogramming. Our previous results [13] indicate that apical tissues in hypocotyl explants need to be reprogrammed into sink tissues; thus, CDF3 may be required for this process. Alternatively, CDF3 could indirectly regulate the auxin-to-CK gradient in apical tissues during de novo shoot regeneration [14] through its auxin-related TF targets.
Several TF genes whose putative orthologs are required for root development in A. thaliana (PLT1/2, BBM, FEZ and SMB) were specifically expressed in the basal region of the hypocotyl during AR formation, suggesting their functional conservation in wound-induced AR formation in tomato. Instead, two PLT genes, PLT3 and PLT7, were expressed in the apical and basal regions of the hypocotyl after 24 hae, which is consistent with their positive role in establishing pluripotency during de novo organ formation, as was previously shown in A. thaliana [40]. We also recovered the expression of the two LBD genes from the SBRL superlocus, Solyc09g066260 and Solyc09g066270, in the basal region of the hypocotyl [17], confirming that SBRL expression can be used as an early marker for AR specification. Two WOX genes, Solyc03g096300 and Solyc02g077390, encoding the putative tomato orthologs of WOX5 and WOX8/9, respectively, were also recovered as early markers for wound-induced AR formation in our study. Recent work identified a CRISPR/Cas9 cis-regulatory allelic series of the tomato WOX8/9 gene that exhibited alterations in inflorescence development, whereas wox8/9 null mutants are embryonically lethal, as in A. thaliana [40]. To confirm the functional role of WOX8/9 expression in wound-induced AR formation in tomato, a phenotypic analysis of the available wox8/9 promoter mutants is needed.
One of the few genes specifically upregulated in the basal region of the hypocotyl at 24 hae was PUCHI (Solyc01g067540), whose expression in A. thaliana is directly regulated by LBD16 in LR founder cells, downstream of the SLR/IAA14-ARF7-ARF19 module [19]. PUCHI is required for the correct timing of meristem establishment by repressing the premature expression of key meristematic genes in the early stages of LR development, such as PLT1, BBM, and WOX5 [41]. PLT1/2 (Solyc11g061750) and BBM (Solyc11g008560) were specifically expressed in the basal region of the hypocotyl explants after 24 hae. The upregulation of PLT1/2 and BBM was correlated with the downregulation of some of their Aux/IAA target genes, such as Solyc06g008590 (SlIAA17), Solyc06g053830 (SlIAA7) and Solyc09g083290 (SlIAA14), which are the putative orthologs of SLR/IAA14 required for specification of LR founder cells in A. thaliana [21]. These results suggest that the specification of AR founder cells in the basal region of tomato hypocotyls after wounding is based on the auxin-mediated SBRL-PUCHI module, while subsequent AR initiation occurs through local activation of the PLT1/2-BBM-WOX5 module. Our results also suggest a negative feedback loop involving PLT1/2-BBM that restricts the specification of AR founder cells in the basal region of the hypocotyl at later time points. Two tomato genes encoding ERF105 (Solyc03g093560) and ERF109 (Solyc10g050970) TFs were specifically expressed in the basal region of the hypocotyl and were strongly downregulated after wounding. While ERF105 has been linked to regulate cold acclimation responses [42], transcriptomic analysis of excised leaves in A. thaliana identified ERF109 as a direct target of jasmonic acid (JA) signaling to establish a transient pulse of auxin production necessary for cell cycle activation and de novo root formation [43,44]. These results are consistent with our findings of a steep methyl-JA gradient in the basal region of the hypocotyl soon after wounding (0-10 hae) that anticipated AR formation [14] and that could explain the observed downregulation of these two ERF TFs.
Our expression profiling analysis allowed us to identify commonly deregulated TF genes in the apical and basal regions of the hypocotyl, which might be involved in early events (i.e., cell dedifferentiation) during de novo organ formation. One of these regulators is Solyc02g082670 (WOX13), since its expression increased significantly at 24 hae in both tissues. A recent report identified WOX13 as a key regulator of callus formation in A. thaliana, downstream of auxin signaling and WIND1 [45,46]. A previous study showed that WOX13 orthologues in the moss Physcomitrium patens are involved in cellular reprogramming at wounding sites [47]. Parallels in WOX13 function between plant species suggest the existence of a conserved GRN involving WOX13 that regulates cell reprogramming after tissue damage. On the other hand, Solyc08g077110, which encodes the NAC DOMAIN CONTAINING PROTEIN 71 (NAC071) and NAC096 homologs in A. thaliana, was also upregulated after wounding in the apical and basal regions of the hypocotyl. In A. thaliana, NAC071 and NAC096 are induced by endogenous auxin in different regions near the wounding in cut stems and act redundantly to promote cell proliferation of cambial cells during tissue reunion [48,49]. Another key regulator of cambial cell division during de novo organ formation could be Solyc11g068960, which putatively encodes the LONESOME HIGHWAY (LHW) TF, which in A. thaliana, regulates formative divisions in vascular cells in roots, leaves, and stems, downstream of auxin [50,51]. Therefore, our transcriptomic analysis and implicit GRN reconstruction could provide candidate TF genes whose function is conserved during de novo organ formation in plants.
We found significant differences in the ploidy profile of cells in the basal and apical regions of the hypocotyl at the time of wounding, suggesting region-specific cell cycle regulation. Our results suggested that basal tissues contained a higher proportion of mitotic and endoreduplicated cells than apical ones at 0 hae, and that a small increase in these two cell populations occurred in the basal region of the hypocotyl during de novo root formation. On the other hand, the apical region was characterized by a strong regulation of the cell cycle after wounding, which is necessary for callus formation during de novo shoot initiation. MI and EI were significantly increased at 24 hae in the apical region of the hypocotyl, which coincided with the upregulation of key genes involved in DNA replication and mitotic regulation, such as those of CycA3 (Solyc04g078310 and Solyc12g088530) and CycP2;1 (Solyc07g052610) cyclins, whose orthologs in A. thaliana are respectively required to control formative cell divisions [52] and the G2/M transition of meristematic cells in response to sugar [25]. Two TF genes, Solyc04g081350 (E2FC) and Solyc03g113760 (E2FE, also known as DEL1), were specifically deregulated in the apical region of the hypocotyl over time (profiles P11 and P13). It has been proposed that E2FC acts as a transcriptional repressor necessary for the switch between cell division and endoreduplication, which could be important in establishing the developmental gradient of cell divisions in growing organs [53]. On the other hand, DEL1 preserves the mitotic state of proliferating cells by suppressing transcription of genes necessary to enter the endocycle [18]. Our GRN analysis identified a DEL1 regulatory node that includes upregulated target genes encoding key components of the replication licensing factor complex, such as ORC6 (Solyc05g007450), MCM4 (Solyc01g110130), MCM6 (Solyc08g082200) and MCM10 (Solyc01g111550), and a DNA primase (Solyc08g082200), among others. These results revealed an unprecedented role of tomato DEL1 TF in the regulation of G1/S transition during de novo organ formation, which deserves further attention.
In our previous study [13], we observed differential regulation of photosynthesis-related genes in the apical and basal regions of hypocotyls after wounding that resembled adaptation to high-light stress. We observed a local increase in ROS production in the apical region at 24 hae, which decreased at the later time point, probably due to enhanced expression of genes encoding PRX and TRX, which are required for chloroplast performance under stress conditions [54]. In this study, we found striking patterns of ROS accumulation in the apical and basal regions of hypocotyl explants over the studied time course, which were not directly correlated with the degree of cell death and therefore might have a signaling role during de novo organ formation. The accumulation of ROS in the basal region of the hypocotyl explants at 96 and 192 hae could be related to AR emergence [55]. During primary root and LR development in A. thaliana, ROS distribution (which is regulated by specific POX enzymes) is required for the transition between cell proliferation and differentiation within the meristem [56,57]. POX performs different functions in addition to its primary role in wound response [58,59]. Some of the POX-encoding genes were strongly deregulated in our RNA-Seq dataset. In A. thaliana, a subset of class III POX genes was targeted by LBD29 [60], a master regulator of auxin-induced callus formation [61]. Furthermore, the biological responses regulated by LBD29 closely resemble those regulated by hypoxia treatment, which has long been considered a crucial factor in determining cell fate in mammals [60]. We found several POX genes (Prx05, Prx11, Prx25, Prx26, Prx38 and Prx68) whose expression was upregulated in the apical region of hypocotyl explants and that were structurally related to the LBD targets described above, suggesting that a GRN involving LBD29-like TF might be conserved in tomato. In addition, the natural variation in hormone-induced shoot regeneration among A. thaliana accessions depended on the expression of a gene encoding a TRX that directly modulates ROS homeostasis [62]. Despite differences in wound responses between animals and plants, ROS is a common signal in both systems [63]. In cancer cells, an increase in ROS stabilizes hypoxia-inducible factor-1α (HIF1α), which enhances glucose catabolism and allows for tumor progression [64]. Furthermore, iPSC reprogramming is enhanced under hypoxic conditions through ROS-mediated activation of HIF1α [65]. A similar hypoxia-mediated mechanism may be present in plants, but negative regulators of ERF-VII TFs may be targeted instead [66]. However, a direct link among ROS production, metabolic reprogramming, and wound-induced organ regeneration in tomato hypocotyl explants awaits experimental confirmation.

4. Materials and Methods

4.1. Plant Material and Growth Conditions

Seedlings of the tomato cultivar “Micro-Tom” were grown in vitro as described elsewhere [12]. Moneymaker (MM), CDF3 RNAi #1.1, 35S::CDF3 #5.2 and 35S::CDF #11.2 lines have been described elsewhere [67]. For LR emergence, 3–4 mm of the root tips in 5-day-old seedlings were excised with a sharp scalpel, and the seedlings were grown for another 2 days; newly emerged LRs were then counted. Hypocotyl explants were obtained after removing the whole root system (2–3 mm above the hypocotyl–root junction) and the shoot apex (just below the cotyledons’ petiole) with a sharp scalpel, at the 100–101 growth stages [68] (0 dae). Hypocotyl explants were transferred to 65 × 150 mm (diameter × height) glass jars with 50 mL of standard growing medium [12]. ARs were periodically counted between 2 and 12 dae. Shoot regeneration stages were quantified as described previously [14]. Data were collected 31 days later (50 das). Three to four millimeters of the apical and the basal regions of the hypocotyl were collected at 0, 24, 96 and 192 hae. Three biological replicates, each consisting of 12 apical or basal fragments were obtained.

4.2. Tomato TF Identification and Annotation

For TF identification, ITAG4 protein sequences were retrieved from the SolGenomics database [69] and were used for prediction of TF genes at iTAK [70], PlantTFDB [71] and PlantTFCat [72] databases. The resulting TF gene list was manually curated through PubMed searches (Table S1). Protein sequences smaller than 90 amino acids or that did not contain a DNA binding domain identified by hmmscan in HMMER [73] were discarded. The A. thaliana orthologs were identified using a reciprocal best BLAST hit strategy (RBBH) with ProteinOrtho [74] using the Araport 11 [75] and tomato ITAG4 proteomes. For auxin-related genes, we used the annotation from [14].

4.3. RNA-Seq Analysis and GRN Generation

RNA-Seq data were obtained from the NCBI SRA database (SRR14598206–SRR14598225) and processed as described in [13]. Briefly, genes with levels below 1 cpm in the 3 libraries were discarded for further analyses, resulting in 22492 expressed genes. DEGs were filtered using a false discovery rate (FDR) < 0.01 and log2 fold change > |1|. Expression analyses were carried out using the iDEP 0.93 web application [15]. Expression profiles were generated using STEM [16]; apical and basal TF genes with ≥0.5 cpm were excluded. For hierarchical clustering and heatmap drawing, the Morpheus web-tool was used [76]. TF genes from selected STEM profiles were used to obtain their putative target genes from the PlantRegMap database [77]. GO enrichment analyses were performed by using the ShinyGO 0.76 web tool [78] with differentially expressed putative targets. GRNs of TFs and their putative targets were constructed using Cytoscape 3.2 [79], and interaction analyses were performed with the NetworkAnalyzer plugin. The annotation of peroxidases in tomato, as well as their assignment in families, was obtained from RedoxiBase [26].

4.4. Cell Cycle Analysis

For flow cytometry, 3–4 mm of the apical or the basal regions of the hypocotyl were collected at 0, 24, 96 and 192 hae. Nuclei were isolated using the chopping method described elsewhere [80]. The suspension was filtered over 40 µm nylon mesh, treated with RNase A (0.8 µg/mL) and stained with propidium iodide (PI; Sigma-Aldrich, St. Lois, MO, USA) (2 µg/mL) for 15 min. Nuclei were analyzed using a SH800C Cell Sorter (Sony Biotechnology, San Jose, CA, USA). Dot-Plot side scatter area (SSC-A) versus forward scatter area (FSC-A), and PI intensity area (PI-A) versus SSC-A were used to locate nuclear populations by particle side. A univariate histogram of PI (564–606 nm emission) was used to represent DNA content 2C, 4C, 8C and large ploidy nuclei on the PE-A axis. A minimum of four biological replicates including 10,000 nuclei were analyzed in each sample.

4.5. Cell Death and ROS Assays

For the microscopic examination of dying cells, hypocotyl explants were harvested at 0, 24, 96 and 192 hae and stained with hot acetic acid/trypan blue solution (0.5% weight/volume (w/v) in 45% acetic acid) for 3 min. After incubation for 1 h at room temperature, the hypocotyls were cleared with chloral hydrate solution (80 g chloral hydrate in 30 mL distilled water). For DAB staining, whole hypocotyl explants were harvested at 0, 24, 96 and 192 hae and immersed in a 1% w/v solution of DAB in Tris-HCl buffer, pH 4.8. After 1 h of vacuum infiltration, the samples were incubated in DAB solution overnight at room temperature. Then, the explants were cleared with chloral hydrate solution. To confirm the specific accumulation of H2O2, we incubated hypocotyl explants with 100 units/µL−1 catalase (Sigma-Aldrich) for 30 min before DAB staining. Thereafter, the samples were processed as indicated above.

5. Conclusions

Through k-means and STEM analysis of bulk RNA-Seq data, we have identified several TF genes with specific spatial and temporal expression patterns, such as WUS, CUC2, SBRL, PUCHI, and WOX5, to be used as reliable markers of the different stages during de novo organ formation in tomato mutants that show defects in tissue regeneration after wounding. We have identified several conserved TF modules, such as those of WIND1-ESR1/2-WOX13 and SLR/IAA14-ARF7/ARF19-SBRL, that are involved in the specification of callus and root identity during de novo shoot and root formation, respectively. In addition, candidate TFs that regulate cambial cell division in response to wounding have been identified. Our results also highlight the fact that a dynamic regulation of cell cycle and ROS could directly contribute to region-specific organ regeneration. Finally, our TF candidates need to be confirmed experimentally using all available genetic tools (i.e., known mutant alleles and CRISPR/Cas9-mediated gene editing).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232416112/s1. References [5,49,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123] are cited in Supplementary Materials.

Author Contributions

Conceptualization, J.M.P.-P. and E.L.; methodology, E.L., M.N.-A. and A.B.S.-G.; formal analysis, E.L. and J.M.P.-P.; investigation, E.L., M.N.-A. and A.B.S.-G.; data curation, E.L. and J.M.P.-P.; writing—original draft preparation, J.M.P.-P.; writing—review and editing, E.L. and M.N.-A.; visualization, E.L. and J.M.P.-P.; supervision, J.M.P.-P.; project administration, J.M.P.-P.; funding acquisition, J.M.P.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministerio de Ciencia e Innovación of Spain, grant numbers BIO2015-64255-R and RTI2018-096505-B-I00, the Conselleria d’Educació, Cultura i Sport of the Generalitat Valenciana, grant numbers IDIFEDER/2018/016 and PROMETEO/2019/117, and the European Regional Development Fund (ERDF) of the European Commission. MNA holds a research contract of the Conselleria d’Educació, Cultura i Sport of the Generalitat Valenciana, grant number EDGJID/2021/025.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw sequence files and read count files are publicly available in the NCBI’s BioProject PRJNA731333. Gene functional annotation is available in the Supplementary Material of this article. All other data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank María Gomáriz Clemente from the SADTI-UMH for her expert technical assistance. We are specially indebted to José Luis Micol (Instituto de Bioingeniería, Universidad Miguel Hernández) for the use of the flow cytometer and Joaquín Medina (Centro de Biotecnología y Genómica de Plantas (CBGP) UPM–INIA) for sharing seeds of CDF3 overexpression and RNAi tomato lines.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ikeuchi, M.; Ogawa, Y.; Iwase, A.; Sugimoto, K. Plant Regeneration: Cellular Origins and Molecular Mechanisms. Development 2016, 143, 1442–1451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ikeuchi, M.; Favero, D.S.; Sakamoto, Y.; Iwase, A.; Coleman, D.; Rymen, B.; Sugimoto, K. Molecular Mechanisms of Plant Regeneration. Annu. Rev. Plant Biol. 2019, 70, 377–406. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, L.; Huang, H. Genetic and Epigenetic Controls of Plant Regeneration. Curr. Top. Dev. Biol. 2014, 108, 1–33. [Google Scholar] [CrossRef]
  4. Mathew, M.M.; Prasad, K. Model Systems for Regeneration: Arabidopsis. Development 2021, 148, dev195347. [Google Scholar] [CrossRef] [PubMed]
  5. Ikeuchi, M.; Shibata, M.; Rymen, B.; Iwase, A.; Bgman, A.M.; Watt, L.; Coleman, D.; Favero, D.S.; Takahashi, T.; Ahnert, S.E.; et al. A Gene Regulatory Network for Cellular Reprogramming in Plant Regeneration. Plant Cell Physiol. 2018, 59, 770–782. [Google Scholar] [CrossRef] [PubMed]
  6. Bae, S.H.; Noh, Y.S.; Seo, P.J. REGENOMICS: A Web-Based Application for Plant REGENeration-Associated TranscriptOMICS Analyses. Comput. Struct. Biotechnol. J. 2022, 20, 3234. [Google Scholar] [CrossRef] [PubMed]
  7. Iwase, A.; Mitsuda, N.; Koyama, T.; Hiratsu, K.; Kojima, M.; Arai, T.; Inoue, Y.; Seki, M.; Sakakibara, H.; Sugimoto, K.; et al. The AP2/ERF Transcription Factor WIND1 Controls Cell Dedifferentiation in Arabidopsis. Curr. Biol. 2011, 21, 508–514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Iwase, A.; Mitsuda, N.; Ikeuchi, M.; Ohnuma, M.; Koizuka, C.; Kawamoto, K.; Imamura, J.; Ezura, H.; Sugimoto, K. Arabidopsis WIND1 Induces Callus Formation in Rapeseed, Tomato, and Tobacco. Plant Signal. Behav. 2013, 8, e27432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Iwase, A.; Mita, K.; Nonaka, S.; Ikeuchi, M.; Koizuka, C.; Ohnuma, M.; Ezura, H.; Imamura, J.; Sugimoto, K. WIND1-Based Acquisition of Regeneration Competency in Arabidopsis and Rapeseed. J. Plant Res. 2015, 128, 389–397. [Google Scholar] [CrossRef] [Green Version]
  10. Iwase, A.; Harashima, H.; Ikeuchi, M.; Rymen, B.; Ohnuma, M.; Komaki, S.; Morohashi, K.; Kurata, T.; Nakata, M.; Ohme-Takagi, M.; et al. WIND1 Promotes Shoot Regeneration through Transcriptional Activation of ENHANCER OF SHOOT REGENERATION1 in Arabidopsis. Plant Cell 2017, 29, 54–69. [Google Scholar] [CrossRef]
  11. Iwase, A.; Kondo, Y.; Laohavisit, A.; Takebayashi, A.; Ikeuchi, M.; Matsuoka, K.; Asahina, M.; Mitsuda, N.; Shirasu, K.; Fukuda, H.; et al. WIND Transcription Factors Orchestrate Wound-Induced Callus Formation, Vascular Reconnection and Defense Response in Arabidopsis. New Phytol. 2021, 232, 734–752. [Google Scholar] [CrossRef] [PubMed]
  12. Alaguero-Cordovilla, A.; Sánchez-García, A.B.; Ibáñez, S.; Albacete, A.; Cano, A.; Acosta, M.; Pérez-Pérez, J.M. An Auxin-Mediated Regulatory Framework for Wound-Induced Adventitious Root Formation in Tomato Shoot Explants. Plant Cell Environ. 2021, 44, 1642–1662. [Google Scholar] [CrossRef]
  13. Larriba, E.; Sánchez-García, A.B.; Martínez-Andújar, C.; Albacete, A.; Pérez-Pérez, J.M. Tissue-Specific Metabolic Reprogramming during Wound-Induced Organ Formation in Tomato Hypocotyl Explants. Int. J. Mol. Sci. 2021, 22, 10112. [Google Scholar] [CrossRef] [PubMed]
  14. Larriba, E.; Belén Sánchez-García, A.; Salud Justamante, M.; Martínez-Andújar, C.; Albacete, A.; Manuel Pérez-Pérez, J.; Tran, P.; Golam Mostofa, M. Dynamic Hormone Gradients Regulate Wound-Induced de Novo Organ Formation in Tomato Hypocotyl Explants. Int. J. Mol. Sci. 2021, 22, 11843. [Google Scholar] [CrossRef]
  15. Ge, S.X.; Son, E.W.; Yao, R. IDEP: An Integrated Web Application for Differential Expression and Pathway Analysis of RNA-Seq Data. BMC Bioinform. 2018, 19, 1–24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Ernst, J.; Bar-Joseph, Z. STEM: A Tool for the Analysis of Short Time Series Gene Expression Data. BMC Bioinform. 2006, 7, 1–11. [Google Scholar] [CrossRef] [Green Version]
  17. Omary, M.; Gil-Yarom, N.; Yahav, C.; Steiner, E.; Hendelman, A.; Efroni, I. A Conserved Superlocus Regulates Above- and Belowground Root Initiation. Science 2022, 375, eabf4368. [Google Scholar] [CrossRef]
  18. Vlieghe, K.; Boudolf, V.; Beemster, G.T.S.; Maes, S.; Magyar, Z.; Atanassova, A.; De Almeida Engler, J.; De Groodt, R.; Inzé, D.; De Veylder, L. The DP-E2F-like Gene DEL1 Controls the Endocycle in Arabidopsis Thaliana. Curr. Biol. 2005, 15, 59–63. [Google Scholar] [CrossRef]
  19. Goh, T.; Toyokura, K.; Yamaguchi, N.; Okamoto, Y.; Uehara, T.; Kaneko, S.; Takebayashi, Y.; Kasahara, H.; Ikeyama, Y.; Okushima, Y.; et al. Lateral Root Initiation Requires the Sequential Induction of Transcription Factors LBD16 and PUCHI in Arabidopsis Thaliana. New Phytol. 2019, 224, 749–760. [Google Scholar] [CrossRef]
  20. Ranftl, Q.L.; Bastakis, E.; Klermund, C.; Schwechheimer, C. LLM-Domain Containing B-GATA Factors Control Different Aspects of Cytokinin-Regulated Development in Arabidopsis Thaliana. Plant Physiol. 2016, 170, 2295–2311. [Google Scholar] [CrossRef]
  21. Goh, T.; Joi, S.; Mimura, T.; Fukaki, H. The Establishment of Asymmetry in Arabidopsis Lateral Root Founder Cells Is Regulated by LBD16/ASL18 and Related LBD/ASL Proteins. Development 2012, 139, 883–893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Aida, M.; Beis, D.; Heidstra, R.; Willemsen, V.; Blilou, I.; Galinha, C.; Nussaume, L.; Noh, Y.S.; Amasino, R.; Scheres, B. The PLETHORA Genes Mediate Patterning of the Arabidopsis Root Stem Cell Niche. Cell 2004, 119, 109–120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ochatt, S.J. Medicago Truncatula Handbook Flow Cytometry (Ploidy Determination, Cell Cycle Analysis, DNA Content per Nucleus). 2006. Available online: https://www.researchgate.net/publication/342002480_Medicago_truncatula_protocols_flow_cytometry_ploidy_determination_cell_cycle_analysis_DNA_content_per_nucleus (accessed on 11 November 2022).
  24. Takahashi, I.; Kojima, S.; Sakaguchi, N.; Umeda-Hara, C.; Umeda, M. Two Arabidopsis Cyclin A3s Possess G1 Cyclin-like Features. Plant Cell Rep. 2010, 29, 307–315. [Google Scholar] [CrossRef] [PubMed]
  25. Peng, L.; Skylar, A.; Chang, P.L.; Bisova, K.; Wu, X. CYCP2;1 Integrates Genetic and Nutritional Information to Promote Meristem Cell Division in Arabidopsis. Dev. Biol. 2014, 393, 160–170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Savelli, B.; Li, Q.; Webber, M.; Jemmat, A.M.; Robitaille, A.; Zamocky, M.; Mathé, C.; Dunand, C. RedoxiBase: A Database for ROS Homeostasis Regulated Proteins. Redox Biol. 2019, 26, 101247. [Google Scholar] [CrossRef]
  27. Devireddy, A.R.; Zandalinas, S.I.; Fichman, Y.; Mittler, R. Integration of Reactive Oxygen Species and Hormone Signaling during Abiotic Stress. Plant J. 2021, 105, 459–476. [Google Scholar] [CrossRef]
  28. Sagi, M.; Davydov, O.; Orazova, S.; Yesbergenova, Z.; Ophir, R.; Stratmann, J.W.; Fluhr, R. Plant Respiratory Burst Oxidase Homologs Impinge on Wound Responsiveness and Development in Lycopersicon Esculentum. Plant Cell 2004, 16, 616–628. [Google Scholar] [CrossRef] [Green Version]
  29. Gill, S.S.; Tuteja, N. Reactive Oxygen Species and Antioxidant Machinery in Abiotic Stress Tolerance in Crop Plants. Plant Physiol. Biochem. 2010, 48, 909–930. [Google Scholar] [CrossRef]
  30. Renau-Morata, B.; Carrillo, L.; Dominguez-Figueroa, J.; Vicente-Carbajosa, J.; Molina, R.V.; Nebauer, S.G.; Medina, J.; Doerner, P. CDF Transcription Factors: Plant Regulators to Deal with Extreme Environmental Conditions. J. Exp. Bot. 2020, 71, 3803–3815. [Google Scholar] [CrossRef]
  31. Reinhardt, D.; Frenz, M.; Mandel, T.; Kuhlemeier, C. Microsurgical and Laser Ablation Analysis of Interactions between the Zones and Layers of the Tomato Shoot Apical Meristem. Development 2003, 130, 4073–4083. [Google Scholar] [CrossRef]
  32. Rodríguez-Leal, D.; Lemmon, Z.H.; Man, J.; Bartlett, M.E.; Lippman, Z.B. Engineering Quantitative Trait Variation for Crop Improvement by Genome Editing. Cell 2017, 171, 470–480.e8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. David-Schwartz, R.; Koenig, D.; Sinha, N.R. LYRATE Is a Key Regulator of Leaflet Initiation and Lamina Outgrowth in Tomato. Plant Cell 2009, 21, 3093–3104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Schiessl, K.; Muiño, J.M.; Sablowski, R. Arabidopsis JAGGED Links Floral Organ Patterning to Tissue Growth by Repressing Kip-Related Cell Cycle Inhibitors. Proc. Natl. Acad. Sci. USA 2014, 111, 2830–2835. [Google Scholar] [CrossRef] [Green Version]
  35. Capua, Y.; Eshed, Y. Coordination of Auxin-Triggered Leaf Initiation by Tomato LEAFLESS. Proc. Natl. Acad. Sci. USA 2017, 114, 3246–3251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Sakamoto, Y.; Kawamura, A.; Suzuki, T.; Segami, S.; Maeshima, M.; Polyn, S.; De Veylder, L.; Sugimoto, K. Transcriptional Activation of Auxin Biosynthesis Drives Developmental Reprogramming of Differentiated Cells. Plant Cell 2022, 34, 4348–4365. [Google Scholar] [CrossRef] [PubMed]
  37. Dai, X.; Liu, N.; Wang, L.; Li, J.; Zheng, X.; Xiang, F.; Liu, Z. MYB94 and MYB96 Additively Inhibit Callus Formation via Directly Repressing LBD29 Expression in Arabidopsis Thaliana. Plant Sci. 2020, 293, 110323. [Google Scholar] [CrossRef]
  38. Domínguez-Figueroa, J.; Carrillo, L.; Renau-Morata, B.; Yang, L.; Molina, R.V.; Marino, D.; Canales, J.; Weih, M.; Vicente-Carbajosa, J.; Nebauer, S.G.; et al. The Arabidopsis Transcription Factor CDF3 Is Involved in Nitrogen Responses and Improves Nitrogen Use Efficiency in Tomato. Front. Plant Sci. 2020, 11, 1825. [Google Scholar] [CrossRef]
  39. Kareem, A.; Durgaprasad, K.; Sugimoto, K.; Du, Y.; Pulianmackal, A.J.; Trivedi, Z.B.; Abhayadev, P.V.; Pinon, V.; Meyerowitz, E.M.; Scheres, B.; et al. PLETHORA Genes Control Regeneration by a Two-Step Mechanism. Curr. Biol. 2015, 25, 1017–1030. [Google Scholar] [CrossRef] [Green Version]
  40. Hendelman, A.; Zebell, S.; Rodriguez-Leal, D.; Dukler, N.; Robitaille, G.; Wu, X.; Kostyun, J.; Tal, L.; Wang, P.; Bartlett, M.E.; et al. Conserved Pleiotropy of an Ancient Plant Homeobox Gene Uncovered by Cis-Regulatory Dissection. Cell 2021, 184, 1724–1739.e16. [Google Scholar] [CrossRef]
  41. Bellande, K.; Trinh, D.C.; Gonzalez, A.A.; Dubois, E.; Petitot, A.S.; Lucas, M.; Champion, A.; Gantet, P.; Laplaze, L.; Guyomarc’H, S. PUCHI Represses Early Meristem Formation in Developing Lateral Roots of Arabidopsis Thaliana. J. Exp. Bot. 2022, 73, 3496–3510. [Google Scholar] [CrossRef]
  42. Bolt, S.; Zuther, E.; Zintl, S.; Hincha, D.K.; Schmülling, T. ERF105 Is a Transcription Factor Gene of Arabidopsis Thaliana Required for Freezing Tolerance and Cold Acclimation. Plant Cell Environ. 2017, 40, 108–120. [Google Scholar] [CrossRef] [PubMed]
  43. Zhang, G.; Zhao, F.; Chen, L.; Pan, Y.; Sun, L.; Bao, N.; Zhang, T.; Cui, C.X.; Qiu, Z.; Zhang, Y.; et al. Jasmonate-Mediated Wound Signalling Promotes Plant Regeneration. Nat. Plants 2019, 5, 491–497. [Google Scholar] [CrossRef] [PubMed]
  44. Zhou, W.; Lozano-Torres, J.L.; Blilou, I.; Zhang, X.; Zhai, Q.; Smant, G.; Li, C.; Scheres, B. A Jasmonate Signaling Network Activates Root Stem Cells and Promotes Regeneration. Cell 2019, 177, 942–956.e14. [Google Scholar] [CrossRef] [PubMed]
  45. Tanaka, H.; Hashimoto, N.; Kawai, S.; Yumoto, E.; Shibata, K.; Tameshige, T.; Yamamoto, Y.; Sugimoto, K.; Asahina, M.; Ikeuchi, M. Auxin-Induced WUSCHEL-RELATED HOMEOBOX13 Mediates Asymmetric Activity of Callus Formation upon Cutting. Plant Cell Physiol. 2022. [Google Scholar] [CrossRef]
  46. Ikeuchi, M.; Iwase, A.; Ito, T.; Tanaka, H.; Favero, D.S.; Kawamura, A.; Sakamoto, S.; Wakazaki, M.; Tameshige, T.; Fujii, H.; et al. Wound-Inducible WUSCHEL-RELATED HOMEOBOX 13 Is Required for Callus Growth and Organ Reconnection. Plant Physiol. 2022, 188, 425–441. [Google Scholar] [CrossRef]
  47. Sakakibara, K.; Reisewitz, P.; Aoyama, T.; Friedrich, T.; Ando, S.; Sato, Y.; Tamada, Y.; Nishiyama, T.; Hiwatashi, Y.; Kurata, T.; et al. WOX13-like Genes Are Required for Reprogramming of Leaf and Protoplast Cells into Stem Cells in the Moss Physcomitrella Patens. Development 2014, 141, 1660–1670. [Google Scholar] [CrossRef] [Green Version]
  48. Asahina, M.; Azuma, K.; Pitaksaringkarn, W.; Yamazaki, T.; Mitsuda, N.; Ohme-Takagi, M.; Yamaguchi, S.; Kamiya, Y.; Okada, K.; Nishimura, T.; et al. Spatially Selective Hormonal Control of RAP2.6L and ANAC071 Transcription Factors Involved in Tissue Reunion in Arabidopsis. Proc. Natl. Acad. Sci. USA 2011, 108, 16128–16132. [Google Scholar] [CrossRef] [Green Version]
  49. Matsuoka, K.; Sato, R.; Matsukura, Y.; Kawajiri, Y.; Iino, H.; Nozawa, N.; Shibata, K.; Kondo, Y.; Satoh, S.; Asahina, M. Wound-Inducible ANAC071 and ANAC096 Transcription Factors Promote Cambial Cell Formation in Incised Arabidopsis Flowering Stems. Commun. Biol. 2021, 4, 1–12. [Google Scholar] [CrossRef]
  50. Vera-Sirera, F.; De Rybel, B.; Úrbez, C.; Kouklas, E.; Pesquera, M.; Álvarez-Mahecha, J.C.; Minguet, E.G.; Tuominen, H.; Carbonell, J.; Borst, J.W.; et al. A BHLH-Based Feedback Loop Restricts Vascular Cell Proliferation in Plants. Dev. Cell 2015, 35, 432–443. [Google Scholar] [CrossRef] [Green Version]
  51. De Rybel, B.; Möller, B.; Yoshida, S.; Grabowicz, I.; Barbier de Reuille, P.; Boeren, S.; Smith, R.S.; Borst, J.W.; Weijers, D. A BHLH Complex Controls Embryonic Vascular Tissue Establishment and Indeterminate Growth in Arabidopsis. Dev. Cell 2013, 24, 426–437. [Google Scholar] [CrossRef]
  52. Willems, A.; Heyman, J.; Eekhout, T.; Achon, I.; Pedroza-Garcia, J.A.; Zhu, T.; Li, L.; Vercauteren, I.; Van den Daele, H.; van de Cotte, B.; et al. The Cyclin CYCA3;4 Is a Postprophase Target of the APC/CCCS52A2 E3-Ligase Controlling Formative Cell Divisions in Arabidopsis. Plant Cell 2020, 32, 2979–2996. [Google Scholar] [CrossRef] [PubMed]
  53. Kobayashi, K.; Suzuki, T.; Iwata, E.; Nakamichi, N.; Suzuki, T.; Chen, P.; Ohtani, M.; Ishida, T.; Hosoya, H.; Müller, S.; et al. Transcriptional Repression by MYB3R Proteins Regulates Plant Organ Growth. EMBO J. 2015, 34, 1992–2007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Cejudo, F.J.; González, M.-C.; Pérez-Ruiz, J.M. Redox Regulation of Chloroplast Metabolism. Plant Physiol. 2020, 186, 9–21. [Google Scholar] [CrossRef] [PubMed]
  55. Mhimdi, M.; Pérez-Pérez, J.M. Understanding of Adventitious Root Formation: What Can We Learn From Comparative Genetics? Front. Plant Sci. 2020, 11, 1541. [Google Scholar] [CrossRef] [PubMed]
  56. Tsukagoshi, H.; Busch, W.; Benfey, P.N. Transcriptional Regulation of ROS Controls Transition from Proliferation to Differentiation in the Root. Cell 2010, 143, 606–616. [Google Scholar] [CrossRef] [Green Version]
  57. Fernández-Marcos, M.; Desvoyes, B.; Manzano, C.; Liberman, L.M.; Benfey, P.N.; del Pozo, J.C.; Gutierrez, C. Control of Arabidopsis Lateral Root Primordium Boundaries by MYB36. New Phytol. 2017, 213, 105–112. [Google Scholar] [CrossRef] [Green Version]
  58. Minibayeva, F.; Beckett, R.P.; Kranner, I. Roles of Apoplastic Peroxidases in Plant Response to Wounding. Phytochemistry 2015, 112, 122–129. [Google Scholar] [CrossRef]
  59. Shigeto, J.; Tsutsumi, Y. Diverse Functions and Reactions of Class III Peroxidases. New Phytol. 2016, 209, 1395–1402. [Google Scholar] [CrossRef] [Green Version]
  60. Xu, C.; Cao, H.; Xu, E.; Zhang, S.; Hu, Y. Genome-Wide Identification of Arabidopsis LBD29 Target Genes Reveals the Molecular Events behind Auxin-Induced Cell Reprogramming during Callus Formation. Plant Cell Physiol. 2018, 59, 749–760. [Google Scholar] [CrossRef]
  61. Fan, M.; Xu, C.; Xu, K.; Hu, Y. LATERAL ORGAN BOUNDARIES DOMAIN Transcription Factors Direct Callus Formation in Arabidopsis Regeneration. Cell Res. 2012, 22, 1169–1180. [Google Scholar] [CrossRef]
  62. Zhang, H.; Zhang, T.T.; Liu, H.; Shi, D.Y.; Wang, M.; Bie, X.M.; Li, X.G.; Zhang, X.S. Thioredoxin-Mediated ROS Homeostasis Explains Natural Variation in Plant Regeneration. Plant Physiol. 2018, 176, 2231–2250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Suzuki, N.; Mittler, R. Reactive Oxygen Species-Dependent Wound Responses in Animals and Plants. Free Radic. Biol. Med. 2012, 53, 2269–2276. [Google Scholar] [CrossRef] [PubMed]
  64. Nagao, A.; Kobayashi, M.; Koyasu, S.; Chow, C.C.T.; Harada, H. HIF-1-Dependent Reprogramming of Glucose Metabolic Pathway of Cancer Cells and Its Therapeutic Significance. Int. J. Mol. Sci. 2019, 20, 238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Hawkins, K.E.; Joy, S.; Delhove, J.M.K.M.; Kotiadis, V.N.; Fernandez, E.; Fitzpatrick, L.M.; Whiteford, J.R.; King, P.J.; Bolanos, J.P.; Duchen, M.R.; et al. NRF2 Orchestrates the Metabolic Shift during Induced Pluripotent Stem Cell Reprogramming. Cell Rep. 2016, 14, 1883–1891. [Google Scholar] [CrossRef] [Green Version]
  66. Pucciariello, C.; Perata, P. The Oxidative Paradox in Low Oxygen Stress in Plants. Antioxidants 2021, 10, 332. [Google Scholar] [CrossRef]
  67. Renau-Morata, B.; Molina, R.V.; Carrillo, L.; Cebolla-Cornejo, J.; Sánchez-Perales, M.; Pollmann, S.; Domínguez-Figueroa, J.; Corrales, A.R.; Flexas, J.; Vicente-Carbajosa, J.; et al. Ectopic Expression of CDF3 Genes in Tomato Enhances Biomass Production and Yield under Salinity Stress Conditions. Front. Plant Sci. 2017, 8, 660. [Google Scholar] [CrossRef] [Green Version]
  68. Feller, C.; Bleiholder, H.; Buhr, L.; Hack, H.; Hess, M.; Klose, R.; Meier, U.; Stauss, R.; Van den Boom, T.; Weber, E. Phänologische Entwicklungsstadien von Gemüsepflanzen: II. Fruchtgemüse Und Hülsenfrüchte. Nachrichtenbl. Deut. Pflanzenschutzd. 1995, 47, 217–232. [Google Scholar]
  69. Sol Genomics Network. Available online: https://solgenomics.net/organism/Solanum_lycopersicum/genome/ (accessed on 11 November 2022).
  70. Zheng, Y.; Jiao, C.; Sun, H.; Rosli, H.G.; Pombo, M.A.; Zhang, P.; Banf, M.; Dai, X.; Martin, G.B.; Giovannoni, J.J.; et al. ITAK: A Program for Genome-Wide Prediction and Classification of Plant Transcription Factors, Transcriptional Regulators, and Protein Kinases. Mol. Plant 2016, 9, 1667–1670. [Google Scholar] [CrossRef] [Green Version]
  71. Jin, J.; Tian, F.; Yang, D.C.; Meng, Y.Q.; Kong, L.; Luo, J.; Gao, G. PlantTFDB 4.0: Toward a Central Hub for Transcription Factors and Regulatory Interactions in Plants. Nucleic Acids Res. 2017, 45, D1040–D1045. [Google Scholar] [CrossRef] [Green Version]
  72. Dai, X.; Sinharoy, S.; Udvardi, M.; Zhao, P.X. PlantTFcat: An Online Plant Transcription Factor and Transcriptional Regulator Categorization and Analysis Tool. BMC Bioinform. 2013, 14, 321. [Google Scholar] [CrossRef] [Green Version]
  73. Hmmscan. Available online: https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan (accessed on 11 November 2022).
  74. Lechner, M.; Findeiß, S.; Steiner, L.; Marz, M.; Stadler, P.F.; Prohaska, S.J. Proteinortho: Detection of (Co-)Orthologs in Large-Scale Analysis. BMC Bioinform. 2011, 12, 124. [Google Scholar] [CrossRef] [PubMed]
  75. Araport 11 Genome Release. Available online: https://www.arabidopsis.org/download/ (accessed on 11 November 2022).
  76. Morpheus. Available online: https://software.broadinstitute.org/morpheus/ (accessed on 11 November 2022).
  77. Tian, F.; Yang, D.C.; Meng, Y.Q.; Jin, J.; Gao, G. PlantRegMap: Charting Functional Regulatory Maps in Plants. Nucleic Acids Res. 2020, 48, D1104–D1113. [Google Scholar] [CrossRef] [PubMed]
  78. Ge, S.X.; Jung, D.; Jung, D.; Yao, R. ShinyGO: A Graphical Gene-Set Enrichment Tool for Animals and Plants. Bioinformatics 2020, 36, 2628–2629. [Google Scholar] [CrossRef] [PubMed]
  79. Cytoscape. Available online: https://cytoscape.org/ (accessed on 11 November 2022).
  80. Galbraith, D.W.; Harkins, K.R.; Maddox, J.M.; Ayres, N.M.; Sharma, D.P.; Firoozabady, E. Rapid Flow Cytometric Analysis of the Cell Cycle in Intact Plant Tissues. Science 1983, 220, 1049–1051. [Google Scholar] [CrossRef] [PubMed]
  81. Ai, G.; Zhang, D.; Huang, R.; Zhang, S.; Li, W.; Ahiakpa, J.K.; Zhang, J. Genome-Wide Identification and Molecular Characterization of the Growth-Regulating Factors-Interacting Factor Gene Family in Tomato. Genes 2020, 11, 1435. [Google Scholar] [CrossRef] [PubMed]
  82. Chini, A.; Ben-Romdhane, W.; Hassairi, A.; Aboul-Soud, M.A.M. Identification of TIFY/JAZ family genes in Solanum lycopersicum and their regulation in response to abiotic stresses. PLoS ONE 2017, 12, e0177381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Gao, Y.; Gao, S.; Xiong, C.; Yu, G.; Chang, J.; Ye, Z.; Yang, C. Comprehensive analysis and expression profile of the homeodomain leucine zipper IV transcription factor family in tomato. Plant Physiol. Biochem. 2015, 96, 141–153. [Google Scholar] [CrossRef]
  84. Huang, S.; Gao, Y.; Liu, J.; Peng, X.; Niu, X.; Fei, Z.; Cao, S.; Liu, Y. Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum. Mol. Genet. Genom. 2012, 287, 495–513. [Google Scholar] [CrossRef]
  85. Huang, Z.; Van Houten, J.; Gonzalez, G.; Xiao, H.; van der Knaap, E. Genome-wide identification, phylogeny and expression analysis of SUN, OFP and YABBY gene family in tomato. Mol. Genet. Genom. 2013, 288, 111–129. [Google Scholar] [CrossRef]
  86. Hu, X.; Zhu, L.; Zhang, Y.; Xu, L.; Li, N.; Zhang, X.; Pan, Y. Genome-wide identification of C2H2 zinc-finger genes and their expression patterns under heat stress in tomato (Solanum lycopersicum L.). PeerJ 2019, 7, e7929. [Google Scholar] [CrossRef] [Green Version]
  87. Jin, J.F.; Wang, Z.Q.; He, Q.Y.; Wang, J.Y.; Li, P.F.; Xu, J.M.; Zheng, S.J.; Fan, W.; Yang, J.L. Genome-wide identification and expression analysis of the NAC transcription factor family in tomato (Solanum lycopersicum) during aluminum stress. BMC Genom. 2020, 21, 288. [Google Scholar] [CrossRef] [PubMed]
  88. Khatun, K.; Robin, A.H.K.; Park, J.-I.; Nath, U.K.; Kim, C.K.; Lim, K.-B.; Nou, I.S.; Chung, M.-Y. Molecular Characterization and Expression Profiling of Tomato GRF Transcription Factor Family Genes in Response to Abiotic Stresses and Phytohormones. Int. J. Mol. Sci. 2017, 18, 1056. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Li, D.; Fu, F.; Zhang, H.; Song, F. Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato (Solanum lycopersicum L.). BMC Genom. 2015, 16, 771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Li, S.; Li, K.; Ju, Z.; Cao, D.; Fu, D.; Zhu, H.; Zhu, B.; Luo, Y. Genome-wide analysis of tomato NF-Y factors and their role in fruit ripening. BMC Genom. 2016, 17, 36. [Google Scholar] [CrossRef] [Green Version]
  91. Li, Z.; Peng, R.; Tian, Y.; Han, H.; Xu, J.; Yao, Q. Genome-Wide Identification and Analysis of the MYB Transcription Factor Superfamily in Solanum lycopersicum. Plant Cell Physiol. 2016, 57, 1657–1677. [Google Scholar] [CrossRef] [Green Version]
  92. Li, X.; Hamyat, M.; Liu, C.; Ahmad, S.; Gao, X.; Guo, C.; Wang, Y.; Guo, Y. Identification and Characterization of the WOX Family Genes in Five Solanaceae Species Reveal Their Conserved Roles in Peptide Signaling. Genes 2018, 9, 260, Erratum in Genes 2018, 9, 457. [Google Scholar] [CrossRef] [Green Version]
  93. Niu, Y.; Zhao, T.; Xu, X.; Li, J. Genome-wide identification and characterization of GRAS transcription factors in tomato (Solanum lycopersicum). PeerJ 2017, 5, e3955. [Google Scholar] [CrossRef] [Green Version]
  94. Wang, J.; Hu, Z.; Zhao, T.; Yang, Y.; Chen, T.; Yang, M.; Yu, W.; Zhang, B. Genome-wide analysis of bHLH transcription factor and involvement in the infection by yellow leaf curl virus in tomato (Solanum lycopersicum). BMC Genom. 2015, 16, 39. [Google Scholar] [CrossRef] [Green Version]
  95. Wang, Y.; Zhang, J.; Hu, Z.; Guo, X.; Tian, S.; Chen, G. Genome-Wide Analysis of the MADS-Box Transcription Factor Family in Solanum lycopersicum. Int. J. Mol. Sci. 2019, 20, 2961. [Google Scholar] [CrossRef] [Green Version]
  96. Xu, R. Genome-wide analysis and identification of stress-responsive genes of the CCCH zinc finger family in Solanum lycopersicum. Mol. Genet. Genom. 2014, 289, 965–979. [Google Scholar] [CrossRef]
  97. Xing, H.; Jiang, Y.; Zou, Y.; Long, X.; Wu, X.; Ren, Y.; Li, Y.; Li, H.L. Genome-wide investigation of the AP2/ERF gene family in ginger: Evolution and expression profiling during development and abiotic stresses. BMC Plant Biol. 2021, 21, 561. [Google Scholar] [CrossRef] [PubMed]
  98. Yuan, Q.; Zhang, C.; Zhao, T.; Yao, M.; Xu, X. A genome-wide analysis of GATA transcription factor family in tomato and analysis of expression patterns. Int. J. Agric. Biol. 2018, 20, 1274–1282. [Google Scholar]
  99. Yu, C.; Cai, X.; Ye, Z.; Li, H. Genome-wide identification and expression profiling analysis of trihelix gene family in tomato. Biochem. Biophys. Res. Commun. 2015, 468, 653–659. [Google Scholar] [CrossRef] [PubMed]
  100. Ohashi-Ito, K.; Saegusa, M.; Iwamoto, K.; Oda, Y.; Katayama, H.; Kojima, M.; Sakakibara, H.; Fukuda, H. A bHLH complex activates vascular cell division via cytokinin action in root apical meristem. Curr. Biol. 2014, 24, 2053–2058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Aliaga Fandino, A.C.; Kim, H.; Rademaker, J.D.; Lee, J.Y. Reprogramming of the cambium regulators during adventitious root development upon wounding of storage tap roots in radish (Raphanus sativus L.). Biol. Open 2019, 8, bio039677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  102. Melnyk, C.W.; Gabel, A.; Hardcastle, T.J.; Robinson, S.; Miyashima, S.; Grosse, I.; Meyerowitz, E.M. Transcriptome dynamics at Arabidopsis graft junctions reveal an intertissue recognition mechanism that activates vascular regeneration. Proc. Natl. Acad. Sci. USA 2018, 115, E2447–E2456. [Google Scholar] [CrossRef] [Green Version]
  103. Han, G.; Wei, X.; Dong, X.; Wang, C.; Sui, N.; Guo, J.; Yuan, F.; Gong, Z.; Li, X.; Zhang, Y.; et al. Arabidopsis ZINC FINGER PROTEIN1 Acts Downstream of GL2 to Repress Root Hair Initiation and Elongation by Directly Suppressing bHLH Genes. Plant Cell 2020, 32, 206–225. [Google Scholar] [CrossRef] [Green Version]
  104. Chandler, J.W. Class VIIIb APETALA2 Ethylene Response Factors in Plant Development. Trends Plant Sci. 2018, 23, 151–162. [Google Scholar] [CrossRef]
  105. Sengupta, S.; Nag Chaudhuri, R. ABI3 plays a role in de-novo root regeneration from Arabidopsis thaliana callus cells. Plant Signal Behav. 2020, 15, 1794147. [Google Scholar] [CrossRef]
  106. Durgaprasad, K.; Roy, M.V.; Venugopal, M.A.; Kareem, A.; Raj, K.; Willemsen, V.; Mähönen, A.P.; Scheres, B.; Prasad, K. Gradient Expression of Transcription Factor Imposes a Boundary on Organ Regeneration Potential in Plants. Cell Rep. 2019, 29, 453–463.e3. [Google Scholar] [CrossRef] [Green Version]
  107. Tang, L.P.; Zhou, C.; Wang, S.S.; Yuan, J.; Zhang, X.S.; Su, Y.H. FUSCA3 interacting with LEAFY COTYLEDON2 controls lateral root formation through regulating YUCCA4 gene expression in Arabidopsis thaliana. New Phytol. 2017, 213, 1740–1754. [Google Scholar] [CrossRef] [PubMed]
  108. Ohashi-Ito, K.; Bergmann, D.C. Regulation of the Arabidopsis root vascular initial population by LONESOME HIGHWAY. Development 2007, 134, 2959–2968. [Google Scholar] [CrossRef] [PubMed]
  109. Ramirez-Parra, E.; Perianez-Rodriguez, J.; Navarro-Neila, S.; Gude, I.; Moreno-Risueno, M.A.; Del Pozo, J.C. The transcription factor OBP4 controls root growth and promotes callus formation. New Phytol. 2017, 213, 1787–1801. [Google Scholar] [CrossRef] [PubMed]
  110. Zhang, M.M.; Zhang, H.K.; Zhai, J.F.; Zhang, X.S.; Sang, Y.L.; Cheng, Z.J. ARF4 regulates shoot regeneration through coordination with ARF5 and IAA12. Plant Cell Rep. 2021, 40, 315–325. [Google Scholar] [CrossRef] [PubMed]
  111. Ye, B.B.; Shang, G.D.; Pan, Y.; Xu, Z.G.; Zhou, C.M.; Mao, Y.B.; Bao, N.; Sun, L.; Xu, T.; Wang, J.W. AP2/ERF Transcription Factors Integrate Age and Wound Signals for Root Regeneration. Plant Cell. 2020, 32, 226–241. [Google Scholar] [CrossRef] [Green Version]
  112. Canher, B.; Lanssens, F.; Zhang, A.; Bisht, A.; Mazumdar, S.; Heyman, J.; Wolf, S.; Melnyk, C.W.; De Veylder, L. The regeneration factors ERF114 and ERF115 regulate auxin-mediated lateral root development in response to mechanical cues. Mol. Plant. 2022, 15, 1543–1557. [Google Scholar] [CrossRef]
  113. De Lucas, M.; Pu, L.; Turco, G.; Gaudinier, A.; Morao, A.K.; Harashima, H.; Kim, D.; Ron, M.; Sugimoto, K.; Roudier, F.; et al. Transcriptional Regulation of Arabidopsis Polycomb Repressive Complex 2 Coordinates Cell-Type Proliferation and Differentiation. Plant Cell. 2016, 28, 2616–2631. [Google Scholar] [CrossRef] [Green Version]
  114. Efroni, I.; Mello, A.; Nawy, T.; Ip, P.L.; Rahni, R.; DelRose, N.; Powers, A.; Satija, R.; Birnbaum, K.D. Root Regeneration Triggers an Embryo-like Sequence Guided by Hormonal Interactions. Cell 2016, 165, 1721–1733. [Google Scholar] [CrossRef] [Green Version]
  115. Heyman, J.; Cools, T.; Canher, B.; Shavialenka, S.; Traas, J.; Vercauteren, I.; Van den Daele, H.; Persiau, G.; De Jaeger, G.; Sugimoto, K.; et al. The heterodimeric transcription factor complex ERF115-PAT1 grants regeneration competence. Nat. Plants 2016, 2, 16165. [Google Scholar] [CrossRef]
  116. Sabatini, S.; Heidstra, R.; Wildwater, M.; Scheres, B. SCARECROW is involved in positioning the stem cell niche in the Arabidopsis root meristem. Genes Dev. 2003, 17, 354–358. [Google Scholar] [CrossRef] [Green Version]
  117. Debernardi, J.M.; Tricoli, D.M.; Ercoli, M.F.; Hayta, S.; Ronald, P.; Palatnik, J.F.; Dubcovsky, J. A GRF-GIF chimeric protein improves the regeneration efficiency of transgenic plants. Nat. Biotechnol. 2020, 38, 1274–1279. [Google Scholar] [CrossRef] [PubMed]
  118. Joshi, S.; Paul, P.; Hartman, J.M.; Perry, S.E. AGL15 Promotion of Somatic Embryogenesis: Role and Molecular Mechanism. Front. Plant Sci. 2022, 13, 861556. [Google Scholar] [CrossRef] [PubMed]
  119. Motte, H.; Verstraeten, I.; Werbrouck, S.; Geelen, D. CUC2 as an early marker for regeneration competence in Arabidopsis root explants. J. Plant Physiol. 2011, 168, 1598–1601. [Google Scholar] [CrossRef] [PubMed]
  120. Gordon-Kamm, B.; Sardesai, N.; Arling, M.; Lowe, K.; Hoerster, G.; Betts, S.; Jones, A.T. Using Morphogenic Genes to Improve Recovery and Regeneration of Transgenic Plants. Plants 2019, 8, 38. [Google Scholar] [CrossRef] [Green Version]
  121. Yang, W.; Choi, M.H.; Noh, B.; Noh, Y.S. De Novo Shoot Regeneration Controlled by HEN1 and TCP3/4 in Arabidopsis. Plant Cell Physiol. 2020, 61, 1600–1613. [Google Scholar] [CrossRef]
  122. Zhai, N.; Xu, L. Pluripotency acquisition in the middle cell layer of callus is required for organ regeneration. Nat. Plants 2021, 7, 1453–1460. [Google Scholar] [CrossRef]
  123. Hassani, S.B.; Trontin, J.F.; Raschke, J.; Zoglauer, K.; Rupps, A. Constitutive Overexpression of a Conifer WOX2 Homolog Affects Somatic Embryo Development in Pinus pinaster and Promotes Somatic Embryogenesis and Organogenesis in Arabidopsis Seedlings. Front. Plant Sci. 2022, 13, 838421. [Google Scholar] [CrossRef]
Figure 1. Region-specific TFs identified during de novo organ formation in tomato hypocotyl explants after wounding. (a) Percentage of TFs specifically expressed in the apical region of the hypocotyl over time. (b) STRING interaction network of A. thaliana orthologs of tomato TFs expressed specifically in the apical region of the hypocotyl. (c) Heatmap of selected TFs expressed in the apical region of the hypocotyl explants. (d) Percentage of TFs specifically expressed in the basal region of the hypocotyl over time. (e) STRING interaction network of A. thaliana orthologs of tomato TFs expressed specifically in the basal region of the hypocotyl. (f) Heatmap of selected TF expressed in the basal region of the hypocotyl explants. Blue/white/red in c and f indicate normalized counts per million reads mapped (cpm) values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Figure 1. Region-specific TFs identified during de novo organ formation in tomato hypocotyl explants after wounding. (a) Percentage of TFs specifically expressed in the apical region of the hypocotyl over time. (b) STRING interaction network of A. thaliana orthologs of tomato TFs expressed specifically in the apical region of the hypocotyl. (c) Heatmap of selected TFs expressed in the apical region of the hypocotyl explants. (d) Percentage of TFs specifically expressed in the basal region of the hypocotyl over time. (e) STRING interaction network of A. thaliana orthologs of tomato TFs expressed specifically in the basal region of the hypocotyl. (f) Heatmap of selected TF expressed in the basal region of the hypocotyl explants. Blue/white/red in c and f indicate normalized counts per million reads mapped (cpm) values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Ijms 23 16112 g001
Figure 2. Expression profiles of TF genes during de novo shoot formation. (a) K-means analysis of the different TFs expressed in the apical region of the hypocotyl over time. (b) Percentage of k-means clusters associated with each expression profile. (ce) Selected expression profiles of downregulated TF genes identified in STEM. Some TFs whose A. thaliana orthologs are known to be involved in regeneration are indicated. The numbers in the left corner indicate the number of TF genes in each expression profile, and the numbers in the right corner indicate the p value; ns: not significant p value. (f) Heatmap of selected TF genes with differential expression in the apical region of the hypocotyl. Blue/white/red in (a,f) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Figure 2. Expression profiles of TF genes during de novo shoot formation. (a) K-means analysis of the different TFs expressed in the apical region of the hypocotyl over time. (b) Percentage of k-means clusters associated with each expression profile. (ce) Selected expression profiles of downregulated TF genes identified in STEM. Some TFs whose A. thaliana orthologs are known to be involved in regeneration are indicated. The numbers in the left corner indicate the number of TF genes in each expression profile, and the numbers in the right corner indicate the p value; ns: not significant p value. (f) Heatmap of selected TF genes with differential expression in the apical region of the hypocotyl. Blue/white/red in (a,f) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Ijms 23 16112 g002
Figure 3. Expression profiles of TF genes during de novo AR formation. (a) K-means analysis of the different TFs expressed in the basal region of the hypocotyl over time. (b) Percentage of k-means clusters associated with each expression profile. (ce) Selected expression profiles of downregulated TF genes identified in STEM. Some TFs whose A. thaliana orthologs are known to be involved in regeneration are indicated. The numbers in the left corner indicate the number of TF genes in each expression profile, and the numbers in the right corner indicate the p value; ns: not significant p value. (f) Heatmap of selected TF genes with differential expression in the basal region of the hypocotyl. Blue/white/red in (a,f) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Figure 3. Expression profiles of TF genes during de novo AR formation. (a) K-means analysis of the different TFs expressed in the basal region of the hypocotyl over time. (b) Percentage of k-means clusters associated with each expression profile. (ce) Selected expression profiles of downregulated TF genes identified in STEM. Some TFs whose A. thaliana orthologs are known to be involved in regeneration are indicated. The numbers in the left corner indicate the number of TF genes in each expression profile, and the numbers in the right corner indicate the p value; ns: not significant p value. (f) Heatmap of selected TF genes with differential expression in the basal region of the hypocotyl. Blue/white/red in (a,f) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S1.
Ijms 23 16112 g003
Figure 4. GRNs for putative TF genes and their deregulated targets associated with P11 and P13 profiles in the apical region of the hypocotyl. (a,b) Selected GRNs from P11 (a) and P13 (b) profiles. (cf) Hierarchical clustering over time of selected TF genes (c,e) and some of their putative targets (d,f) from P11 (c,d) and P13 (e,f) profiles. Blue/white/red in (cf) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Tables S1 and S4.
Figure 4. GRNs for putative TF genes and their deregulated targets associated with P11 and P13 profiles in the apical region of the hypocotyl. (a,b) Selected GRNs from P11 (a) and P13 (b) profiles. (cf) Hierarchical clustering over time of selected TF genes (c,e) and some of their putative targets (d,f) from P11 (c,d) and P13 (e,f) profiles. Blue/white/red in (cf) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Tables S1 and S4.
Ijms 23 16112 g004
Figure 5. GRNs for putative TF genes and their deregulated targets associated with P11 and P15 profiles in the basal region of the hypocotyl. (a,b) Selected GRNs from P11 (a) and P15 (b) profiles. (ce) Hierarchical clustering over time of selected putative targets from P11 (c,d) and P15 (e) profiles. Blue/white/red indicate in (ce) normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Tables S1 and S4.
Figure 5. GRNs for putative TF genes and their deregulated targets associated with P11 and P15 profiles in the basal region of the hypocotyl. (a,b) Selected GRNs from P11 (a) and P15 (b) profiles. (ce) Hierarchical clustering over time of selected putative targets from P11 (c,d) and P15 (e) profiles. Blue/white/red indicate in (ce) normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Tables S1 and S4.
Ijms 23 16112 g005
Figure 6. Regulation of the cell cycle during de novo organ formation. (ae) Ploidy analysis by flow cytometry. (a) Bivariate scatter plot of forward scatter height (FSC-H) versus propidium iodide (PI) fluorescence. (b) Histogram of the DNA content of a representative sample from the apical region at 0 hae. (c) Proportion of nuclei with different DNA content in the apical (white bars) or the basal (black bars) region at 0 hae; asterisks indicate significant differences (p value < 0.001) between the regions studied. (de) Mitotic index (MI; (d)) and endoreduplication index (EI; (e)) over time; letters in (d,e) indicate significant differences (Fisher’s least significant difference (LSD); p value < 0.001) between samples. (f,g) Expression profile of DEG involved in cell cycle regulation. A set of 172 putative core cell cycle genes involved in the different phase transitions were studied (Table S5). Blue/white/red in (f,g) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S5.
Figure 6. Regulation of the cell cycle during de novo organ formation. (ae) Ploidy analysis by flow cytometry. (a) Bivariate scatter plot of forward scatter height (FSC-H) versus propidium iodide (PI) fluorescence. (b) Histogram of the DNA content of a representative sample from the apical region at 0 hae. (c) Proportion of nuclei with different DNA content in the apical (white bars) or the basal (black bars) region at 0 hae; asterisks indicate significant differences (p value < 0.001) between the regions studied. (de) Mitotic index (MI; (d)) and endoreduplication index (EI; (e)) over time; letters in (d,e) indicate significant differences (Fisher’s least significant difference (LSD); p value < 0.001) between samples. (f,g) Expression profile of DEG involved in cell cycle regulation. A set of 172 putative core cell cycle genes involved in the different phase transitions were studied (Table S5). Blue/white/red in (f,g) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S5.
Ijms 23 16112 g006
Figure 7. Regulation of ROS production during de novo organ formation. (a) H2O2 accumulation visualized by DAB staining in hypocotyl explants over time. (b) Proposed pathways of ROS production and detoxification with the indicated key enzyme activities: RBOH, NADPH oxidase/respiratory burst oxidase homologue; SOD, superoxide dismutase; POX, peroxidase; GPX, glutathione peroxidase; PRX/TRX, peroxiredoxins/thioredoxins; APX, ascorbate peroxidases; CAT, catalases. (c) Localized cell death in hypocotyl explants was visualized by trypan blue staining. (df) Hierarchical clustering of DEGs encoding class III peroxidases. Blue/white/red in (dg) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S6. Scale bars (a,c): 3 mm.
Figure 7. Regulation of ROS production during de novo organ formation. (a) H2O2 accumulation visualized by DAB staining in hypocotyl explants over time. (b) Proposed pathways of ROS production and detoxification with the indicated key enzyme activities: RBOH, NADPH oxidase/respiratory burst oxidase homologue; SOD, superoxide dismutase; POX, peroxidase; GPX, glutathione peroxidase; PRX/TRX, peroxiredoxins/thioredoxins; APX, ascorbate peroxidases; CAT, catalases. (c) Localized cell death in hypocotyl explants was visualized by trypan blue staining. (df) Hierarchical clustering of DEGs encoding class III peroxidases. Blue/white/red in (dg) indicate normalized cpm values in each row, where red indicates the highest abundance of transcripts. Gene annotations are found in Table S6. Scale bars (a,c): 3 mm.
Ijms 23 16112 g007
Figure 8. Wound-induced AR formation and shoot regeneration after apical cutting in CDF3 tomato mutants. (a) AR emergence of CDF3 overexpression lines (blue and purple), CDF3 RNAi lines (orange) and Moneymaker (MM; green). (b) Rooting capacity of shoot explants at 5 dae. (c) Weight of the apical region of the hypocotyl after 31 dae. (d) Representative images of regenerated explants at 31 dae. Letters in (b,c) indicate significant differences (p value < 0.01) between genotypes. Asterisks in (a) indicate significant differences (Fisher’s LSD; p value < 0.05) between genotypes. Scale bars: 25 mm.
Figure 8. Wound-induced AR formation and shoot regeneration after apical cutting in CDF3 tomato mutants. (a) AR emergence of CDF3 overexpression lines (blue and purple), CDF3 RNAi lines (orange) and Moneymaker (MM; green). (b) Rooting capacity of shoot explants at 5 dae. (c) Weight of the apical region of the hypocotyl after 31 dae. (d) Representative images of regenerated explants at 31 dae. Letters in (b,c) indicate significant differences (p value < 0.01) between genotypes. Asterisks in (a) indicate significant differences (Fisher’s LSD; p value < 0.05) between genotypes. Scale bars: 25 mm.
Ijms 23 16112 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Larriba, E.; Nicolás-Albujer, M.; Sánchez-García, A.B.; Pérez-Pérez, J.M. Identification of Transcriptional Networks Involved in De Novo Organ Formation in Tomato Hypocotyl Explants. Int. J. Mol. Sci. 2022, 23, 16112. https://doi.org/10.3390/ijms232416112

AMA Style

Larriba E, Nicolás-Albujer M, Sánchez-García AB, Pérez-Pérez JM. Identification of Transcriptional Networks Involved in De Novo Organ Formation in Tomato Hypocotyl Explants. International Journal of Molecular Sciences. 2022; 23(24):16112. https://doi.org/10.3390/ijms232416112

Chicago/Turabian Style

Larriba, Eduardo, Míriam Nicolás-Albujer, Ana Belén Sánchez-García, and José Manuel Pérez-Pérez. 2022. "Identification of Transcriptional Networks Involved in De Novo Organ Formation in Tomato Hypocotyl Explants" International Journal of Molecular Sciences 23, no. 24: 16112. https://doi.org/10.3390/ijms232416112

APA Style

Larriba, E., Nicolás-Albujer, M., Sánchez-García, A. B., & Pérez-Pérez, J. M. (2022). Identification of Transcriptional Networks Involved in De Novo Organ Formation in Tomato Hypocotyl Explants. International Journal of Molecular Sciences, 23(24), 16112. https://doi.org/10.3390/ijms232416112

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop