Correction: Functional studies of plant transcription factors and their relevance in the plant root-knot nematode interaction

Jose Domínguez FigueroaAlmudena Gómez-RojasCarolina Escobar^*

• ENVIRONMENTAL SCIENCES, Faculty of Environmental Sciences and Biochemistry, University of Castilla-la Mancha, Toledo, Spain

Root-knot nematodes are polyphagous parasitic nematodes that cause severe losses in the agriculture worldwide. They enter the root in the elongation zone and subtly migrate to the root meristem where they reach the vascular cylinder and establish a feeding site called gall. Inside the galls they induce a group of transfer cells that serve to nurture them along their parasitic stage, the giant cells. Galls and giant cells develop through a process of post-embryogenic organogenesis that involves manipulating different genetic regulatory networks within the cells, some of them through hijacking some molecular transducers of established plant developmental processes, such as lateral root formation or root regeneration. Galls/giant cells formation involves different mechanisms orchestrated by the nematode´s effectors that generate diverse plant responses in different plant tissues, some of them include sophisticated mechanisms to overcome plant defenses. Yet, the plant-nematode interaction is normally accompanied to dramatic transcriptomic changes within the galls and giant cells. It is therefore expected a key regulatory role of plant-transcription factors, coordinating both, the new organogenesis process induced by the RKNs and the plant response against the nematode. Knowing the role of plant-transcription factors participating in this process becomes essential for a clear understanding of the plant-RKNs interaction and provides an opportunity for the future development and design of directed control strategies. In this review, we present the existing knowledge of the TFs with a functional role in the plant-RKN interaction through a comprehensive analysis of current scientific literature and available transcriptomic data. Key words: plant-RKNs interaction, galls, giant cells, transcription factors, new organogenesis, plant defense, plant-development   INTRODUCTION Plant-parasitic nematodes have the ability to infect a wide range of host plants from which they feed depleting their resources, resulting in significant economic losses in agricultural production worldwide (Singh et al., 2015; Kikuchi et al., 2017). Among these destructive pathogens, the endoparasitic Root-Knot Nematodes (RKNs; Meloidogyne spp.) are one of the most economically impactful (Elling, 2013). RKNs, use their stylet and a diverse range of effectors to invade the plant roots and initiate the formation of specialized feeding cells known as giant cells (GCs). These GCs are contained within a novel pseudo-organ called gall that constitutes their feeding site (Escobar et al., 2015). While the significant role of the pericycle in gall formation is well-established from experiments with transgenic lines that induce chemical ablation, the precise origin of the GCs precursor cells remains not fully understood. However, some evidence point to their origin from precursor cells of the pericycle, xylem and/or vascular cambium (Cabrera et al., 2014b, Olmo et al., 2017; Olmo et al., 2020). GCs undergo mitosis accompanied by incomplete cytokinesis and DNA endoreduplication forming a multinucleated cell with greatly increased volume and a dense cytosol. Moreover, GCs also show fragmented vacuoles, undergo cell wall modifications, and ultimately develop membrane invaginations, becoming transfer cells to nourish the nematode (de Almeida Engler & Favery, 2011; Cabrera et al., 2014b; Escobar et al., 2015). RKNs employ sophisticated mechanisms to overcome plant defenses and modulate the host biochemistry and physiology (Kikuchi et al., 2017). They manipulate different genetic regulatory programs of the plant cells, including the cell cycle, various developmental programs, and stress responses, in order to undergo new post-embryogenic organogenesis leading to gall formation. For instance, Meloidogyne javanica infection alters pathways involved in de novo organogenesis leading to feeding site formation by interfering with auxins signaling cascades (Cabrera et al., 2014b). Consequently, a substantial transcriptional response is triggered (e.g., in Arabidopsis thaliana; Jammes et al., 2005; Fuller et al., 2007; Barcala et al., 2010; Silva et al., 2022). Due to the mentioned dramatic transcriptomic changes described in galls, a pivotal regulatory role of transcription factors (TFs) is therefore expected, coordinating both, the new organogenesis process induced by the RKNs and the plant response against the nematode. Therefore, understanding the role of TFs during RKN infection is essential for a deeper understanding of the plant-nematode interaction and for the development of future control strategies. In this review, we have explored the existing knowledge on TFs involved in the plant-RKN interaction through a comprehensive analysis of current scientific literature and available transcriptomic data, the latter focused on Arabidopsis as considerable transcriptomic data is available and it was shown to be a good model of the plant-RKN interaction (Gheysen and Fenoll, 2011). Our aim is to compile the existing knowledge regarding the crucial role of TFs in the orchestration of the transcriptional response activated within the plant after RKN infection. Transcriptional profiling of Transcription Factors families in Arabidopsis Several transcriptomic analyses have been performed to investigate mRNA population changes during RKNs establishment and gall formation in Arabidopsis plants (Jammes et al., 2005; Fuller et al., 2007; Barcala et al., 2010; Silva et al., 2022). These studies have provided valuable insights into the genetic and transcriptional dynamics associated with gall development. Functional classification of differentially expressed genes (DEGs) revealed RNA-related pathways as one of the groups with a high number of DEGs (Barcala et al., 2010), which are mostly involved in biological processes such as transcriptional regulation. Therefore, we analyzed the data contained in the NEMATIC database (NEMatode-Arabidopsis-Transcriptomic-Interaction-Tool; Cabrera et al., 2014a) that includes the most representative transcriptomic experiments of the RKN interaction in Arabidopsis, and the data from a recent RNAseq of galls 3 days post-infection (dpi) in Arabidopsis (Silva et al., 2022). These data show that of the 1717 annotated TF loci in the Arabidopsis genome based on the criteria of the Plant Transcription Factor Data Base (PlantTFDB; Jin et al., 2014), 834 TFs are differentially expressed (DE) in one or more experiments, that correspond to approximately 49% of the total known TFs in Arabidopsis. Among the 58 families classified according to PlantTFDB (Jin et al., 2014), 53 of them have DE members at some stage of gall formation (91%; Figure 1A), 52 TFs families at early stage (3 dpi) and 37 at medium-late stages (7, 14 and 21 dpi), (89% and 64%, respectively; Figure 1A). In GCs at 3 dpi, 29 TF family members were DE (50%; Figure 1B). Only 5 TFs families did not show DE members in either GCs or any of the gall stages. This indicates that most of the TF families are DE at same stage of gall and/or GCs formation, which presumably should have a great impact in the dramatic transcriptional changes described in galls (see introduction). Figure 1C-E also shows the percentage of DE TFs within the top 36 TF families with the highest number of DE members in early and mid-late stage galls and GCs. Six of these belong to TF superfamilies in which all TF members were included. The predominant families in all three transcriptomes were MYB, bHLH, ERF, NAC and WRKY. The role of several members of these TF families during the plant-nematode interaction was analysed and is our focus throughout the manuscript. While all these data indicate a substantial involvement of TFs in regulating the transcriptional responses of plants to RKN infection, it is worth noting that the functional roles of only over 30 Arabidopsis TFs and about a dozen in tomato (Solanum lycopersicum) have been investigated (Table 1). Therefore, our understanding of the regulatory networks orchestrated by TFs in plant-RKN interaction remains rather limited. Plant Transcription Factors with a role in plant-defense during the RKNs interaction. Plants have developed a multitude of defense mechanisms to counter potential pathogen attacks. The two primary plant immune responses are PAMP-triggered immunity (PTI), which is initiated by the recognition of receptors that recognize pathogen-associated molecular patterns (PAMPs), and effector-triggered immunity (ETI), in which pathogen effectors are recognized by plant resistance proteins known as R proteins (Peng et al., 2018). Both signaling pathways trigger similar molecular processes such as MAPK cascades, the production of reactive oxygen species, secondary metabolites and an increase in the biosynthesis of hormones such as salicylic acid (SA), jasmonic acid (JA), abscisic acid (ABA) and ethylene (ETH; Isah, 2019; Peng et al., 2018). However, PTI induces rapid and transient activation of MAPKs to enhance local immune responses without triggering plant cell death, whereas ETI results in prolonged and sustained MAPK activity, usually leading to hypersensitive response and programmed cell death (Tsuda et al., 2013). However, both PTI and ETI pathways seem to be interconnected as it has recently been described that the activation of either PTI or ETI alone is not sufficient for effective resistance to the bacterial plant-pathogen Pseudomonas syringae. Thus, both immune responses mutually potentiate to activate strong defences against pathogens (Ngou et al., 2021). In any case, plant defense involves a complex interconnected signaling network that ensures a precise transcriptional response. Consequently, several TFs have been identified as crucial for fine-tuning the plant's transcriptional immune response (Birkenbihl et al., 2017). In this respect, TFs from different families, such as WRKY, MYB, AP2 and bZIP, typically induced in response to various biotic and abiotic stresses (Ambawat et al., 2013; Jiang et al., 2017) were also differentially expressed during the RKNs interaction (e.g., in Arabidopsis thaliana; Jammes et al., 2005; Fuller et al., 2007; Silva et al., 2022) including the GCs (Barcala et al., 2010). However, plant-parasitic nematodes, like other pathogens such as bacteria, fungi, oomycetes and some insects, secrete proteins and small molecules, called effectors, to suppress or evade host defense responses and alter host cell structure and function to their advantage, thereby facilitating nematode establishment (Rutter et al., 2022). This section of the review explores the existing literature on defense-related transcription factors (TFs) and their role in the context of the plant-RKNs interaction. The WRKY family is one of the largest TF families found exclusively in plants. Its members play critical roles in various plant processes, encompassing growth, development, abiotic and biotic stress responses, and plant innate immunity (Wani et al., 2021). Some of them are key components in pathways responsible for PTI and ETI activation (Ribeiro et al., 2022). In an RNAseq data of tomato roots after M. javanica infection (15 dpi; Chinnapandi et al., 2017), several WRKYs were identified as negative regulators of the defence response. Among the up-regulated genes, SlWRKY45 was further studied using a promoter::GUS reporter line. The line showed considerable activation at 5 dpi, which continued through feeding-site development and gall maturation (15 and 28 dpi; Table 1), in line with the RNAseq data. Two independent overexpressing lines of 35S:SlWRKY45 showed an increase in nematode infection and GCs area compared to the control lines, although the number of GCs within the galls was not affected (Chinnapandi et al., 2017, Table 1). Consistently, qRT-PCR revealed a down-regulation of defence-related genes, which are typical markers of JA and SA-mediated pathways, such as those encoding pathogenesis-related (PR-1) and proteinase inhibitor II (Pin2) proteins respectively, which could explain the increased nematode infection in the overexpressing line. In this respect, it has been recently described that SlWRKY45 interact with JA-ZIM domain family proteins that are key repressors of the JA signalling, and it is also able to bind and inhibit the activity of the promoter of the JA biosynthesis gene ALLENE OXIDE CYCLASE (AOC) (Huang et al., 2022). All of this is consistent with the attenuated resistance to Meloidogyne incognita of SlWRKY45 overexpression and confirm its role as a negative regulator fo the defense response (PTI) against Meloidogyne spp. Additionally, SlWRKY45 is upregulated by cytokinins, and its overexpression caused the repression of the cytokinin response factor 1 (CRF1) and CRF6 (Chinnapandi et al., 2017). RKNs and cyst nematodes (CNs) have the ability to synthesize and secrete cytokinins, as noted by De Meutter et al. (2003). Furthermore, it has been demonstrated that the CN Heterodera shachtii has a functional cytokinin-synthesizing isopentenyltransferase gene which is essential for virulence and feeding site expansion (Siddique, et al., 2015).The secretion of nematode cytokinins could potentially disrupt the balance of plant hormones and cytokinin signalling. Therefore, SlWRK45 may play a crucial role in coordinating hormone signals that promote nematode development within the root tissue. Similarly, recent studies have identified SlWRKY16 and SlWRKY31 as negative regulators of plant immunity and defence (PTI) in tomato. Both genes were induced during nematode infection until late infection stages (28dpi). Overexpression of these genes in tomato lines using R. rhizogenes-mediated transformation under the control of the CaMV35S promoter resulted in increased susceptibility to M. javanica, as evidenced by enhanced galling and reproduction parameters (Kumar et al., 2023; Table 1). In contrast to SlWRKY45, 16, and 31, which act as negative regulators of the plant defence response to RKNs, other members of the WRKY family have been described as positive regulators of defence against RKNs. For instance, WRKY11, whose expression is induced 24 hours after infection with M. incognita, and WRKY17, which can function in partial redundancy with WRKY11, are associated with the activation of basal defence mechanisms (PTI). The Arabidopsis lines wrky11, wrky17, and wrky11/wrky17 exhibited increased susceptibility to M. incognita. This was evident from the significantly higher number of galls observed 4 weeks after inoculation in both the single and double mutant lines compared to the wild-type plants. However, there were no significant differences between the single mutants and the double mutant, indicating that these two TFs do not function redundantly in this pathogenic interaction (Teixeira et al., 2016; Table 1). Similarly, mutant lines wrky11 and wrky17 showed more susceptibility in Arabidopsis infected with the CN Heterodera schachtii (Ali et al., 2014), indicating commonalities in the basal resistance mechanisms between both plant-(RKNs and CNs) interactions. Additionally, WRKY11pro::GUS lines showed that the WRKY11 promoter was activated in the root elongation zone and root tip 24 hai with M. incognita (Table 1). Furthermore, in assays based on treatments with crude extracts of J2 larvae, GUS activity was also detected in roots, mainly in the elongation zone where RKNs invade roots. The promoter activity was detected in the absence of any mechanical damage produced by RKNs penetration and intercellular migration (Teixeira et al., 2016). This suggests that the gene is an early responder to the presence of the nematode. Similarly, GUS staining was observed restricted to the root elongation zone of MYB51pro::GUS plants early after infection with M. incognita (Table 1; Figure 2). MYB51 is a member of the MYB Transcription Factor family and, together with MYB34, regulates glucosinolate biosynthesis. Accordingly, the Arabidopsis myb51 myb34 double mutant, which is completely impaired in glucosinolate production (Frerigmann & Gigolashvili, 2014), displayed a higher number of galls compared to wild-type plants in M. incognita infection assays. This increased susceptibility of the double mutant point to a positive role of glucosinolates in RKNs defense (Teixeira et al., 2016; Table 1; Figure 2). Importantly, WRKY11 and MYB51 are involved in BRASSINOSTEROID INSENSITIVE-ASSOCIATED KINASE 1 (BAK1)-dependent PTI responses activated by RKN infection, however, another BAK-independent immune signaling pathway was detected probably involved in the camalexin biosynthetic pathway (Teixeira et al., 2016). It is well-established that nematodes can suppress defense-related genes in feeding cells. In this respect, transcriptomic data of Arabidopsis and tomato micro-dissected GCs induced by M. javanica at early infection stages (3dpi) showed a down-regulation of MYB34 in Arabidopsis and its ortholog in tomato, as reported by Barcala et al. (2010) and Portillo et al. (2013), respectively. In contrast, the RNAseq analysis conducted by Silva et al. (2022) at the same infection stage, showed that MYB34 is induced in whole galls compared to control non-infected root segments. These findings suggest that RKNs may inhibit glucosinolate production in feeding cells, while a defense response is likely maintained in the remaining gall tissues. Other WRKY members related to the basal defense against RKNs as positive regulators are AtWRKY72 and its orthologs in Solanum lycopersicum, SlWRKY72a and SlWRKY72b (Bhattarai et al., 2010; Table 1; Figure 2). Arabidopsis WRKY72 mutant lines showed a significant increase in egg masses compared to Col-0. Tomato roots of the susceptible cultivar (cv) Moneymaker (mi-1/mi-1), with SlWRKY72a and/or SlWRKY72b transient silenced or co-silenced showed similar results, confirming the involvement of SlWRKY72 in basal defense (PTI) in both plant species (Bhattarai et al., 2010; Table 1). In addition, data obtained with the tomato resistant cv. Motelle (Mi-1/Mi-1) indicated that SlWRKY72 is also involved in gene-for-gene resistance (ETI) during RKN infection of tomato roots. Thus, a significant increase of SlWRKY72a expression level was observed in response to M. incognita 12 and 36 hours after inoculation (hai), as well as of SlWRKY72b expression at 12 hai in tomato roots of the resistant cv. Motelle (Mi-1/Mi-1), whereas this was not observed in the susceptible cv. Moneymaker (mi-1/mi-1Bhattarai et al., 2010). Functional confirmation of the participation of both TFs in tomato gene-for-gene resistance was obtained by transient silencing or co-silencing SlWRKY72a and SlWRKY72b in the resistant Motelle roots. This resulted in an increased susceptibility to RKNs, while no infection was observed in the control Motelle roots. The expression changes of another WRKY family member in tomato, SlWRKY70, also suggest a putative role in Mi-1-mediated resistance (ETI) , as its mRNA was up-regulated in both susceptible and non-susceptible lines and RKNs were able to infect and reproduce in WRKY70 transiently silenced tomato roots of the resistant cv. Motelle, whereas no infection and reproduction was observed in non-silenced Motelle tomato roots (Atamian et al., 2012; Table 1; Figure 2). In addition, WRKY70 has been described in Arabidopsis to mediate basal and R gene defence in response to aphids (Knoth et al., 2007). Therefore, both WRKY70 and WRKY72 showed a conserved role as activators of defense in different plant species in response to different pests and pathogens (Knoth et al., 2007; Bhattarai et al., 2010; Atamian et al., 2012). However, they appear to mediate defense through different signaling pathways; in Arabidopsis, WRKY70 acts as a mediator between the SA and JA defense pathways during biotic stress (Li et al., 2006); a similar scenario might be happening in tomato as SlWRKY70 is up-regulated during the first hours after SA treatment, but repressed after Methyl jasmonate (MeJA) application (Atamian et al., 2012; Table 1). On the other hand, WRKY72 seems to be involved in plant defence against pathogens in a SA independent manner as genes deregulated in Atwrky72 did not respond to SA analogues (Bhattarai et al 2010). During the interaction with RKN, other members of the WRKY transcription factor family in tomato, such as WRKY3 and 35, were induced and functional, i.e., SlWRKY3pro::GUS and SlWRKY35pro::GUS lines showed GUS signal at early infection stages with M. javanica (2 and 5 dpi). Furthermore, the overexpression of SlWRKY3 in transgenic hairy root lines led to a decrease in the reproduction of M. javanica. This was accompanied by an increase in the accumulation of defence molecules from the shikimate and oxylipin pathways. On the other hand, SlWRKY3 RNAi lines promoted reproduction compared to wild-type plants, confirming its role in basal defence (PTI) against RKNs (Chinappandi et al., 2019; Table 1; Figure 2). Moreover, SlWRKY80 was also recently identified as a positive regulator that affects Mi-1-mediated resistance (ETI). In Motelle lines carrying the Mi-1 resistance gene, SlWRKY80 transcript levels increased by more than 3.21 and 4.56-fold at 3 and 6 dpi, respectively, compared to the control reference M82 variety. The expression of SlWRKY80 was also significantly induced by the defence hormones JA and SA in tomato Motelle. Furthermore, virus-induced gene silencing (VIGS) assays demonstrated that silencing SlWRKY80 in the resistant Motelle tomato resulted in a significant increase in the number of egg masses in the roots of individual plants and a significant decrease in its resistance level. This confirms the role of SlWRKY80 as a positive regulator of Mi-1-mediated tomato resistance (Nie et al., 2023; Table 1; Figure 2). Finally, it is possible that WRKY family members are involved in defence responses induced during the infection by RKNs in monocotyledonous species. This is supported by a transcriptomic study of galls at 3 and 7 dpi of Oryza sativa infected with M. graminicola that revealed upregulation of OsWRKY34, OsWRKY36 and OsWRKY62 (Kyndt et al., 2012; Table 1; Figure 2). Yet, further analysis is required to determine the functional implications of these TFs in rice responses against RKNs. Plant transcription factors with a dual link to stress caused by RKNs and plant developmental programs Although, the structure of this review makes a sharp separation between those TFs involved in plant defense and those involved in plant development, the boundaries are not always clear. In recent years, a connection between stress signaling and developmental programs, such as root regeneration, has been extensively described (Ikeuchi et al., 2019). One example of the RKNs interaction involves the ERF109 and ERF115 transcription factors, which belong to the Ethylene Responsive Factor (ERF) subfamily within the APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) superfamily (Sakuma et al., 2002; Wu et al., 2022). These transcription factors are involved in a core molecular network triggered by wound-induced JA, which induces stem cell activation and regeneration of Arabidopsis thaliana roots (Zhou et al., 2019). ERF109 and ERF115 play a crucial role in maintaining the quiescence of the root stem organizer cells, also known as the quiescent center (QC). ERF109 is activated transcriptionally within minutes in response to JA and wounding and operates upstream of ERF115. Conversely, ERF115 operates upstream of the protein complex RBR-SCR, which regulates the asymmetric cell divisions of root stem cells, QC quiescence, and the activation of the QC regulatory protein WOX5. In addition, it is worth noting that ERF115 is activated not only by JA but also by auxin signaling, which is crucial in galls (Cabrera et al., 2014b; Olmo et al., 2020). It is interesting to observe that GUS staining assays revealed the induction of both ERF109 and ERF115 promoters after the infection of the M. incognita as early as 1 dpi. Furthermore, time course experiments have shown that ERF109 is induced during nematode penetration, while ERF115 is strongly induced in vascular and/or endodermal cells at all stages from penetration and feeding site initiation until gall formation. Interestingly the repressor activity of ERF115 in the dominant-negative ERF115-SRDX line resulted in a loss of root growth recovery capacity after infection with M. incognita compared to Col-0. Furthermore, the erf109 mutant and ERF115-SRDX line exhibited reduced susceptibility to infection and developed fewer egg masses compared to Col-0 seedlings 7 weeks after inoculation. Consistently, gall formation in ERF115-SRDX roots progressed slower than in Col-0, and DNA synthesis at feeding sites was less active compared to Col-0 (Zhou et al., 2019). Therefore, the JA-induced regeneration pathway, with ERF109 and ERF115 acting as regulators, stimulates root growth following nematode infection and contributes to the reproductive success of M. incognita (Zhou et al., 2019; Table 1, Figure 2). In this respect, it is relevant to mention that the damage caused by H. schachtii during invasion also activates a jasmonate-dependent ERF109 pathway, promoting secondary root formation (Guarneri et al., 2022). Therefore, new root-growth and/or regeneration in both RKNs and CNs interaction seems to be mediated by common JA responsive transducers as ERF109. Furthermore, it has recently been reported that ERF115 interacts with PAT1, a transcription factor belonging to the PHYTOCHROME A TRANSDUCTION 1 (PAT1) GRAS subfamily, forming heterodimers. Additionally, erf115, pat1-2, and erf115/pat1-2 double mutants showed GCs often with less cytoplasm, fewer nuclei, and less ploidy than the wild-type lines, which probably caused a delay in nematode development. Furthermore, overexpression of ERF115 (ERF115OE), ERF114, and PAT1 resulted in accelerated gall induction and downstream activation of key players in the regenerative pathway, possibly related to the high cell division rates observed, particularly in the ERF115OE line, following mechanical stress induced by RKNs. In conclusion, the ERF115-PAT1 complex contributes to the regenerative potential of nematode-induced galls by facilitating tissue healing, thereby maintaining the gall's functionality until maturation and nematode reproduction (Ribeiro et al., 2023; Table 1; Figure 2). ERF6 is another stress-related gene that does not appear to regulate plant defenses. It is a member of the Ethylene Responsive Factors (ERF) family and has a unique AP2/ERF domain (Nakano et al., 2006). It was identified from a quantitative trait loci (QTL) study that investigated the relationship between allelic variation in specific loci (QTLs) and the susceptibility of Arabidopsis to M. incognita (Warmerdam et al., 2019). qPCR analysis showed a significant down-regulation of ERF6 in Arabidopsis wild type seedling roots at 7 dpi in association with nematode infections. A reproduction test was also conducted on the erf6-1 mutant line infected with M. incognita at 7 dpi, which resulted in a significant increase of 28% in egg masses per plant compared to the wild type line. In addition, microarray analysis revealed that there were 489 differentially expressed genes in the roots of erf6-1 nematode-infected plants compared to infected wild-type plants (Warmerdam et al., 2019). Previous studies have shown that ERF6 is phosphorylated by MPK3/MPK6 during biotic stress, which activates the expression of Jas/ETH defense genes such as PDF1.2a and PDF1.2b, thereby enhancing Arabidopsis' defense against fungal infections (Meng et al., 2013; Wang et al., 2022). However, in erf6 RKNs-infected roots, the expression of PDF1.1 and PDF1.2 was not significantly altered compared to wild-type infected plants. This is similar to the expression of other pathogenesis-related genes, such as PR1, PR2, PR3, and PR4, which are known to be regulated by ERFs. These findings suggest that ERF6 does not suppress plant defenses during the RKNs interaction. Many of the genes that are differentially regulated in the roots of nematode-infected erf6-1 plants at 7 dpi are involved in responses to abiotic stresses, such as osmotic stress. This suggests that ERF6 regulates abiotic stress responses in the plant-nematode interaction (Warmerdam et al., 2019; Table 1; Figure 2). Plant transcription factors induced during RKNs infection with impact in plant-development The root-knot nematodes induce feeding cells, GCs, with enlarged nuclei within their heterogeneous feeding sites or feeding organs called galls, which indicates increased DNA replication cycles (De Almeida Engler et al., 1999). This process is called endoreduplication and occurs when successive phases of DNA synthesis (S) follow each other without intervening mitosis or cytokinesis. Endopolyploidy is observed in differentiated and enlarged plant cells, such as Arabidopsis trichomes, endosperm, and fruit (Sabelli et al., 2007). Somatic polyploidy is particularly prevalent in higher plants. A high-resolution DNA endoploidy map of the developing Arabidopsis root has revealed the importance of endoreduplication for the expression of genes encoding cell-wall modifying enzymes. These enzymes are crucial during GC development, suggesting that these responses may serve as a buffering system for stress conditions (Wieczoreck, 2015; Bhosale et al., 2018). In this respect, the TF E2Fe/DEL1 is an inhibitor of endoreduplication that maintains the mitotic state of proliferating cells by suppressing transcription of genes necessary to enter the endocycle (Dimova and Dyson, 2005; Inze et al., 2006). The timing of cell cycle exit and onset of endoreduplication is determined by the levels of E2Fe/DEL1, which control anaphase-promoting activator genes such as CCCS52A2 (Lammens et al., 2008). Arabidopsis has three DEL genes (DEL1, DEL2 and DEL3). Loss of DEL1 function results in increased ploidy, while ectopic expression of DEL1 results in decreased endoreduplication levels and cells are prone to rapid expansion (Vlieghe et al., 2005; Ramirez-Parra et al., 2004; Lammens et al., 2008). The role of E2Fe/DEL in GCs formation induced by RKNs was analyzed. Ectopic expression of DEL1 resulted in morphological changes of the GCs within the galls, as the GCs were smaller with profuse cell wall invaginations and smaller nuclei than in the wild type line at 21 dpi. Furthermore, the DEL1 overexpressing line exhibited a significant reduction in the number of M. incognita egg masses due to an induction of the mitotic state, resulting in severe impairment of reproduction. Conversely, the del1-1 loss of function line displayed small and malformed GCs with reduced mitotic activity and little cytoplasm, possibly due to a premature initiation of the endocycle (De Almeida-Engler et al., 2012; Table 1; Figure 2). The results indicate that multiple nuclei resulting from acytokinetic mitotic events are not sufficient to drive GC expansion. Therefore, during the plant-RKNs interaction, there is a reprogramming of the plant cell cycle machinery, inducing mitotic cycles in GCs followed by repeated endoreduplication cycles, both of which are necessary for correct GC development (De Almeida-Engler et al., 2012; Table 1; Figure 2). Although only data from the variation of expression levels during RKNs infection and in situ localization of its transcripts are available, the orthologue DEL1 gene from Arabidopsis in woody plants, such as Prunus sogdiana, PsoDEL1, also appears to be negatively correlated with endoreduplication and growth of GCs. PsoDEL1 exhibited weak expression in the feeding sites during the early stages of infection (3 dpi). As the infection progressed, the hybridization signal was barely detected at the feeding site and within the GCs at 7, 14 and 21 dpi (Xiao et al., 2020; Table 1; Figure 2). Therefore, it is highly probable that the role of DEL1 is conserved in distant plant species during feeding site and GCs formation. Interestingly, E2Fe/DEL1 also plays a role in SA biosynthesis as a transcriptional repressor of a member of the isochorismate pathway, ENHANCED DISEASE SUSCEPTIBILITY 5 (EDS5), that encodes a SA transporter required for elevated SA immunity in Arabidopsis (Chandran et al., 2014). Repression of genes involved in plant defenses is a characteristic of the compatible interaction between RKNs and dicotyledonous or monocotyledonous species, such as Arabidopsis and rice, particularly in GCs (Barcala et al., 2010; Ji et al., 2013). Furthermore, the identification of M. incognita effectors, such as Mi-ISC-1, confirms that the nematode actively deploys a functional isochorismatase to suppress SA-mediated plant defences by altering the isochorismate synthase pathway for SA biosynthesis, favoring parasitism (Qin et al., 2022). The del1-1 mutant showed SA accumulation in 7 dpi Arabidopsis galls, while in uninfected roots, the SA levels did not change significantly compared to the control background Col-0. Therefore, DEL1 seems to repress SA biosynthesis in RKNs-induced galls in Arabidopsis (Nakagami et al., 2020; Table 1; Figure 2). As a result, del1-1 mutant galls at 5 dpi exhibited more intense lignin staining than wild type galls, and the expression patterns of genes encoding enzymes related to lignin biosynthesis, such as 4-coumarate: CoA ligase 1 (4CL1), 4CL2, alcohol dehydrogenase 5 (CAD5), phenylalanine ammonia-lyase 1 (PAL1), PAL2, and cinnamate 4-hydroxylase (C4H) were significantly up-regulated as compared to wild type galls (Nakagami et al., 2020). In conclusion, the loss of function of DEL1/E2Fe in Arabidopsis galls leads to enhanced resistance. However, it also causes growth inhibition, likely due to excessive lignification and/or SA accumulation in RKNs-induced galls. Therefore, DEL1 may mediate a balance between growth and defense by limiting the accumulation of SA in the infection sites (Nakagami et al., 2020), similar to what was reported during fungal infection in leaves (Chandran et al., 2014). Plant MYB3R TFs, which are homologous to Myb oncoproteins, are also involved in controlling mitosis and cytokinesis progression by recognizing Mitotic-specific activator (MSA) elements present in genes expressed during the G2 and M-phase, such as B-cyclins (Ito et al., 1998, 2001; Menges et al., 2005). The MYB transcription factor family is a large family found in all eukaryotes and is involved in various processes that control plant development, metabolism, and responses to biotic and abiotic stress (Dubos et al., 2010). In Arabidopsis, five MYB3R genes have been identified. Among them, MYB3R1 and MYB3R4 function redundantly as activators, with MYB3R4 contributing more to the activation of G2/M phase-specific genes (Haga et al., 2007, 2011). However, MYB3R1 has a dual function as it can act as a repressor along with MYB3R3 and MYB3R5. Meanwhile, MYB3R4 is only expressed during the G2/M phase, whereas the repressor-type MYB3Rs are active in post-mitotic cells and proliferative cells outside the G2/M phase. In Arabidopsis, MYB3R4::GUS lines showed a GUS signal in the vasculature, where a weak expression of MYB3R5::GUS was detected. Following M. incognita infection, the MYB3R4::GUS line exhibited GUS signal in the centre of 3, 5 and 7 dpi galls, while the MYB3R5::GUS line showed two strands of signal surrounding the centre of 3 and 5 dpi galls that decreased and became patchy at 7 dpi. In contrast, no signal was detected in the roots of the MYB3R3::GUS line, whether infected or uninfected. The loss of function lines myb3r1, myb3r5 and myb3r1/3/5 showed no effect on nematode infection. In contrast, myb3r3, myb3r4, myb3r1/4 and myb3r3/5 showed a significant reduction in the number of galls compared to wild type plants (Suzuki et al., 2021a; Table 1; Figure 2). It is known that the myb3r1/4 mutant presents aberrant cytokinesis and down-regulation of cell cycle genes (Haga et al., 2007, 2011), and MYB3R4 is involved in endoreduplication, acting as an activator or repressor depending on its phosphorylation state (Chandran et al., 2010). In addition, MYB3Rs proteins can form complexes during cell cycle progression. For example, MYB3R4 interacts with RBR1 (Retinoblastoma-related) and E2FB, while MYB3R3 interacts with RBR1 and E2FC, which are necessary for endoreduplication (Del Pozo et al., 2002, 2006). In this respect, as mentioned before, the overexpression of a transcription factor of the DP-E2F-like family (E2Fe/DEL1) that maintains the mitotic state of proliferating cells, caused increased mitotic activity and consequent endocycle inhibition in the galls formed by RKNs. Thus, the feeding cells within the galls showed multiple nuclei and inhibited cell expansion affecting nematode development (De Almeida Engler et al., 2012; 2013).Therefore, although a direct interaction of MYB3Rs proteins with E2F members has not yet been described in the RKN interaction, the role of MYB3Rs activators in the RKNs-interaction may be related to the regulation of key cellular processes during cell cycle progression in galls/GCs development. One of the characteristics observed in the transcriptomes of Arabidopsis galls and micro-dissected GCs are the high number of genes included in categories such as development or hormone metabolism, both directly related (Barcala et al., 2010; Silva et al., 2022). Experimental data has confirmed that gall and GCs formation share TFs that are molecular components of transduction pathways involved in lateral root and callus formation, as well as other plant developmental processes such as tuberization, nodulation, fruit development, and flowering time (Cabrera et al., 2014b; Diaz-Manzano et al., 2018; Medina et al., 2017; Olmo et al., 2019, 2020). One of the initial examples discussed is LBD16, a member of the LATERAL ORGAN BOUNDARIES-DOMAIN TF family. LBD16 is a crucial component of the auxin pathway that leads to cell divisions in the xylem pole pericycle, which are necessary for lateral root formation (Goh et al., 2019) and is also involved in pluripotency acquisition in callus cells (Liu et al., 2018). LBD16 is activated early during nematode establishment in xylem pole pericycle cells near the nematode head and by nematode secretions in protoplast. Within the forming galls its expression was maintained till medium infection stages (11 dpi) as indicated by a promoter-GUS fusion. It also showed a crucial function during RKNs establishment, but not during the establishment of CNs as loss of function lines of LBD16 (lbd16-1, 35S::LBD16-SRDX; pLBD16::LBD16-SRDX) showed significant less infections by M. javanica than the control wild type line. LBD16 is also important for the GCs and galls development, as smaller galls and less expanded GCs were observed than in Col-0 in some of those mentioned loss of function lines. Interestingly, LBD16 is regulated by auxins in galls as also described during lateral root formation (Cabrera et al., 2014b; Table 1; Figure 2). Unexpectedly, LBD16 was locally induced in the vascular tissue of leaves after RKNs infection, as it was proven that M. javanica is able to establish, induce GCs, and reproduce in in vitro cultured Arabidopsis leaves. LBD16 is also essential for feeding site formation in leaves, as evidenced by the inability of RKNs to establish in the 35S::LBD16-SRDX line, which contains the LBD16 coding sequence fused to the transcriptional repressor domain SRDX driven by the 35S promoter (Olmo et al., 2017; see Table 1). Thus, LBD16 appears to be a conserved molecular hub connecting developmental signals with those necessary for RKNs feeding site formation in Arabidopsis. A role for LBD16 in feeding site formation induced by CNs has not been described. However, LBD16 is induced in a WOX11-dependent manner in the primordia of adventitious lateral roots that are promoted after H. schchatii infection (Willing et al., 2024). This finding connects the plant-responses to both, RKNs and CNs, to the activation of critical transducers of root developmental programs. ABERRANT LATERAL ROOT FORMATION 4 (ALF4) is another gene relevant to lateral root formation and gall development. It encodes a nuclear-localized protein that is not a transcription factor. However, due to its localization and participation upstream of the auxin signaling pathways leading to lateral root formation, we are mentioning it (DiDonato et al., 2004; Bagchi et al., 2017). It is also involved in developmental processes such as vascular vessels reconnection in grafting, hormone-induced callus formation or de novo root organogenesis from leaf explants (Melnyk et al., 2015; Sujimoto et al., 2010). ALF4 was induced at very early infection stages of infection by M. javanica (1dpi), as indicated by the activation of a pALF4:.GUS construct. The GUS signal increased at 4 dpi and it was maintained at medium stages of gall development (7–10 dpi), but disappeared at 14 dpi. ALF4 is necessary for the proper development of galls and GCs formed by Meloidogyne spp in Arabidopsis, as the mutant alf4-1 presents aberrant galls and GCs with severe structural abnormalities that cause a dramatic reduction in the nematode’s reproduction (Olmo et al., 2017). PUCHI, a member of the ERF Transcription factor family, is also involved in the formation of new organs such as lateral roots or floral formation (Hirota et al., 2007; Karim et al., 2009). It is activated by auxin through LBD16, controlling lateral root primordium patterning (Hirota et al., 2007; Goh et al., 2019). Interestingly, it is up-regulated in galls at early-mid stages (3, 5 and 7 dpi; Yamaguchi et al., 2017; Table 1; Figure 2) in line with its promoter activity (Suzuki et al., 2021b; Table 1; Figure 2). The expression peaked at 5 dpi, but there was no significant difference in the number of galls at 2 dpi or the number of egg masses in the mutant line puchi-1 compared to the control wild type line. These results suggest that PUCHI does not play a significant role in nematode invasion, gall formation, or nematode reproduction (Suzuki et al., 2021b; Table 1; Figure 2). However, the function of PUCHI during nematode infection may be related to cell wall morphology. This is suggested by the observation that the puchi-1 mutant line displayed aberrant giant cells (GCs) with dramatic protrusions and invaginations containing thick cell walls that were not present in galls from the wild type line (Suzuki et al., 2021b). Trinh et al. (2019) reported that PUCHI controls the biosynthesis pathway of very long-chain fatty acid (VLCFA) during lateral root formation. Additionally, Uemura et al. (2023) found that VLCFAs can modify the cell wall through the activation of MYB93, which regulates cell wall genes. RNAseq and promoter::GUS activation assays revealed up-regulation of genes encoding 3-KETOACYL-COA SYNTHASE 1 (KCS1) and KCS20, enzymes implicated in very-long-chain fatty acid (VLCFA) synthesis in galls. Their expression peaked at around 5 dpi, similar to that of PUCHI. Additionally, GCs from the mutants kcs1-5 and puchi-1 exhibited a similar phenotype with thicker walls and protuberances compared to wild-type galls. Therefore, the observed phenotype in the puchi-1 mutant may be attributed to modifications in the VLCFA composition of the cell wall and cell membrane of the GCs (Suzuki et al., 2021b). However, PUCHI does not significantly affect nematode infectivity or reproductive parameters. The formation of galls by RKNs is a process of post-embryonic new organogenesis as new structures specialized for nematode nourishment are induced by the nematode into the vascular cylinder of the host plant roots. The study of two TFs involved in common signaling pathways for lateral root formation, AUXIN-RESPONSIVE-FACTOR-5/ARF5, a key factor for root stem-cell niche regeneration and lateral root initiation, and GATA-TRANSCRIPTION FACTOR-23/GATA23 that specifies pluripotent founder cells during lateral root formation (De Rybel et al., 2010; De Smet et al., 2010; Efroni et al., 2016) shed light on the plant transduction pathways used or hijacked by the nematode to achieve those dramatic reprogramming events. The impact on nematode infection, galls, and GCs development was significant in the arf5-2 mutant, as well as in inducible knockout lines for ARF5, and in a knockdown line of GATA23 as compared to wild type lines. pGATA23::GUS, was induced at early infection stages, 3 dpi-7 dpi, but at 14 dpi no signal was detected, whereas ARF5::GUS was active in a shorter window, i.e., at 3 dpi a clear signal was detected that faded at 7dpi, this confirmed their induction at early-mid infection according to their putative roles during galls/GCs formation(Olmo et al., 2020; Table 1; Figure 2). Therefore, the results suggest that transient pluripotency reprogramming, which leads to lateral root founder cell-like specification and root regeneration, is also necessary for gall/GCs organogenesis. In contrast, other TFs that are the main upstream transducers during lateral root development, such as ARF7 and ARF19 (Okushima et al., 2007), did not exhibit a significant role or specific expression pattern during gall/GCs formation (Olmo et al., 2020; Table 1; Figure 2). However, the regulation of another auxin TF in galls, ARF3, was also similar to that of lateral root growth (Marin et al., 2010; Cabrera et al., 2016). ARF3 is a TF that participates in a regulatory module. In this module, miR390a controls the biogenesis of TAS3-derived tasiRNAs that regulate the auxin responsive factors ARF2, ARF3 and ARF4 by degrading their transcripts and controlling lateral root growth (Marin et al., 2010). Two sensor lines (pARF3:ARF3-GUS and a tasiRNA-resistant ARF3 line, pARF3:ARF3m-GUS) indicated the binding of TAS3-derived tasiRNAs to the ARF3 sequence in galls. The results strongly suggest that the promoters of miR390 and TAS3 are active, and their products are functional in galls, repressing ARF3 (Cabrera et al., 2016; Table 1). Therefore, silencing of ARF3 seems to be important during gall development and establishment of RKNs. Recently, other ARFs, such as ARF8A and ARF8B, have been studied in tomato. These were induced during the early to mid-infection stages (7-14 dpi) in galls and GCs of tomato transgenic lines pARF8A:GUS and pARF8B:GUS. The up-regulation of ARF8A/B transcripts in galls compared to uninfected roots in transcriptomic analysis (RNAseq) is due to the high activity of their promoter combined with reduced silencing by miR167. Furthermore, the mutant lines slarf8a, slarf8ab, and slarf8ab showed severely compromised infection and reproduction of M. incognita. In addition, expression of ARF8A and ARF8B is required for correct giant development as the former mutant lines showed giant cells significantly smaller than in the wild type line (Noureddine et al., 2023; Table 1; Figure 2). All these data, support a key role for ARF8s in feeding site formation. Following the robust hypothesis that gall/GCs formation is a new organogenesis process, and the described similarities with callus formation, it is important to note that callus formation involves the differentiation of pericycle or pericycle-like cells in a process that resembles root tip organization. Thus, crucial root meristem (RAM) TFs marker genes, namely SHORTROOT/SHR, SCARECROW/SCR, SCHIZORIZA/SCZ, and WUSCHEL-RELATED-HOMEOBOX-5/WOX5 are expressed (Sugimoto et al., 2010). It also requires the ectopic activation of a lateral root developmental program and consequently the expression of LBD genes (Sugimoto et al., 2010). It is noteworthy that these genes were induced very early during gall formation (2-5 dpi), but no signal was detected 7-8 dpi and most of them also played important roles in the establishment of RKNs. Thus, the activation of plant developmental programs that promote transient pluripotency/stemness leads to the generation of quiescent center and meristematic-like cell identities within the vascular cylinder of galls (Olmo et al., 2020; Table 1; Figure 2). Moreover, a process of new organogenesis also involves revascularisation, which is crucial for maintaining GCs growth as they are symplastically isolated specialized transfer cells (Rodiuc et al., 2014; Hoth et al., 2008). Phloem formation is induced during gall development (Bartlem et al., 2014). APL (Altered phloem development), a MYB coiled-coil-type TF involved in phloem identity acquisition, is expressed in protophloem, metaphloem and companion cells (Bonke et al., 2003). It is induced early after infection with M. incognita and M. javanica in Arabidopsis, as shown by an APL::GUS line with a strand signal in 3 dpi galls that increases at 5 dpi (Suzuki et al., 2021; Table 1; Figure 2). However, functional tests are still needed to confirm its role during gall formation. It is known that the balance between cell proliferation and cell differentiation in the procambium is regulated downstream of the receptor and kinase cascade by the WOX4 TF. Cyst nematodes have been shown to modulate the procambial cell proliferation of feeding cell formation probably by mimicking the plant B-type CLE TDIF (tracheary element differentiation inhibitory factor) peptide that is encoded by two genes CLE41 and CLE44 in Arabidopsis, and by taking control of the (TDIF RECEPTOR/PHLOEM INTERCALATED WITH XYLEM (TDR/PXY)-WOX4 signaling pathway (Guo et al., 2017). In this regard, WOX4 and ATHB8, typical procambium marker genes, were induced in M. incognita galls at 3, 5 and 7 dpi. However, the analysis of athb8-11 and wox4-1 loss-of-function mutants did not cause any visible effect on the infection parameters or gall diameter. These genes usually function redundantly; therefore, single mutations were probably not sufficient to prevent procambial cell formation (Yamaguchi et al., 2017; Table 1; Figure 2). Nevertheless, the connection between gall formation and different developmental pathways is evident. The riboregulator miRNA172 post-transcriptionally targets a small group of regulatory repressor genes, including APETALA2 (AP2) and AP2-like genes, such as TARGET OF EARLY ACTIVATION TAGGED 1 (TOE1). These miRNA172-targeted AP2-like TFs are involved in controlling several developmental processes, such as plant aging, flowering time, tuber formation, fruit growth, and nodulation (Martin et al., 2009; Zhu & Helliwell, 2010; Yan et al., 2013; Wang et al., 2014; Ripoll et al., 2015). Functional analysis of RKNs infective and reproductive parameters was conducted on Arabidopsis lines with altered activities based on 35S::MIMICRY172 (MIM172), 35S::TARGET OF EARLY ACTIVATION TAGGED 1 (TOE1)-miR172-resistant (35S::TOE1R) and mutant (flowering locus T-10 (ft-10)) during gall and GCs development. The results indicated that the regulatory module miRNA172/TOE1/FT plays a crucial role during GCs and gall development (Díaz-Manzano et al., 2018; Table 1; Figure 2). Therefore, the repression of TOE1 by miRNA172 is relevant for the normal establishment of RKNs and the formation of galls/GCs. Interestingly, the SPL7/MIR408-UCC2/MIR398-CSD1 copper module (Griffiths-Jones et al., 2007) is also functional and active within galls (Noureddine et al., 2022; Table 1; Figure 2).). Loss of function lines of miR398b/c and miR408 in Arabidopsis, resulted in fewer galls and smaller infection sites as compared with the control lines. These findings together with the expression data of two microRNA families, miR398 and miR408, upregulated in galls, similarly to that of the TF SLP7, strongly suggest that the expression of MIR408 and MIR398B and -C is activated by SPL7 in response to a decrease in copper concentration in galls (Noureddine et al., 2022). The role of this module might be related to its involvement in lignin metabolism (Reyt et al., 2020) as the cell wall suffer dramatic changes during gall and GCs development and numerous cell wall modifying enzymes are activated (Wieczoreck, 2015). However, further research will be needed to elucidate it. CONCLUSION Despite the abundance of DEGs encoding plant TFs during the RKNs interaction in several plant species, that for example in the Arabidopsis transcriptomes cover most of the TFs families identified within the genome (Figure 1), their biological function during RKNs infection is still poorly understood. The functional role of only around 40 TFs have been assessed (Table 1). Most of the data was obtained from plant lines with altered expression for each TF, mainly loss of function lines, that were infected, and significant differences either in the infection rate, gall formation, gall or GCs development compared to their control wild type lines were encountered. As a result, a clear phenotype during RKNs infection was identified. Two main groups of functionalities can be identified: TFs related to plant defenses, whose downstream targets are defense-related genes (see Table 1 and Figure 2), and TFs involved in plant developmental programs, such as lateral root and/or callus formation, root apical meristem, or root regeneration. These are presumably hijacked by nematodes for their own benefit, including some TFs with a role in basic cellular functions, such as cell cycle control (Table 1 and Figure 2). Interestingly, there is increasing evidence of the involvement of TFs with a dual role in plant development and defense and/or as integrator hubs between stress signals and developmental signals, such as DEL1, ERF109, ERF115, ERF114 and ERF6 (Table1; Figure 2). However, the signal transduction pathways regulated by those TFs during RKNs infection are mostly unknown or only partially understood. Nevertheless, few regulatory modules involving TFs have been fully or partially proven in the interaction between RKNs and plants. These include the miRNA172/TOE1/FT module (Díaz-Manzano et al., 2018) and the SPL7/MIR408-UCC2/MIR398-CSD1 copper module (Noureddine et al., 2022). Clearly, further studies are needed to increase our knowledge in the regulatory networks driven by plant TFs modified by the nematode and as a plant response during the RKN interactions. AUTHOR CONTRIBUTIONS JD-F and AG-R performed specialized literature search and participated in the main writing of the manuscript, as well as in the figures and table design. CE coordinated and participated in the writing of the manuscript and in the correction of the manuscript. All the authors contributed to the critical reading of the manuscript. FUNDING The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Spanish Government (PID2019-105924RBI00 and PID2022-138989OB-I00, Ministerio de Ciencia e Innovacion, MCIN/AEI/10.13039/501100011033 FEDER, UE), CPP2021-008347 2022-2024, MICIN/AEI /10.13039/ 501100011033 and by the European Union Next Generation EU/PRTR, by the Castilla-La Mancha Government and EU FEDER (SBPLY/17/180501/000287; SBPLY/21/180501/000033), to CE. AG-R was a recipient of a grant from Plan Propio I+D+I of UCLM co-financed by FSE [2019/5964] and JD-F received funding from Ministerio de Universidades-European Union in the frame of NextGeneration EU RD 289/2021 (Universidad Politécnica de Madrid).