- 1Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
- 2Center for Plant Systems Biology, VIB, Ghent, Belgium
- 3iRIP Unit, Laboratory of Microbiology, Department of Biochemistry and Microbiology, Ghent University, Ghent, Belgium
The type III secretion system (T3SS) is a well-studied pathogenicity determinant of many bacteria through which effectors (T3Es) are translocated into the host cell, where they exercise a wide range of functions to deceive the host cell’s immunity and to establish a niche. Here we look at the different approaches that are used to functionally characterize a T3E. Such approaches include host localization studies, virulence screenings, biochemical activity assays, and large-scale omics, such as transcriptomics, interactomics, and metabolomics, among others. By means of the phytopathogenic Ralstonia solanacearum species complex (RSSC) as a case study, the current advances of these methods will be explored, alongside the progress made in understanding effector biology. Data obtained by such complementary methods provide crucial information to comprehend the entire function of the effectome and will eventually lead to a better understanding of the phytopathogen, opening opportunities to tackle it.
1. Introduction
1.1. The perpetrator: The Ralstonia disease
Ralstonia solanacearum is a bacterial phytopathogen that poses serious threats to agriculture due to its broad host range and its worldwide distribution. The disease, commonly known as bacterial wilt or brown rot, has been reported in more than 65 countries, affecting more than 200 plant species from over 50 botanical families, among which various solanaceous crops, such as tomato (Solanum lycopersicum) and potato (Solanum tuberosum) and a wide range of ornamentals (Mansfield et al., 2012; European and Mediterranean Plant Protection Organization, 2023). For these reasons, it was ranked as the second most important bacterial phytopathogen in molecular plant pathology worldwide (Mansfield et al., 2012).
Classification of the genetically diverse R. solanacearum species complex (RSSC) has been very dynamic over the years. In recent literature, the main classification used has been the phylotype distinction, because it provides evolutionary relationships and is based on the similarities in the DNA sequences of the 16S-23S internally transcribed spacer region, besides the hypersensitive response and pathogenesis B (hrpB), endoglucanase (egl), and mutator S (mutS) genes (Guidot et al., 2007; Remenant et al., 2010). The RSSC has been classified into four different phylotypes, reflecting its geographic origin, namely phylotype I, II, III, and IV strains, originating from Asia, America, Africa, and Indonesia, Australia, and Japan, respectively. Phylotype II strains are further subclassified into IIa and IIb. A detailed comparison of the genomes by average nucleotide identity (ANI) from different phylotypes revealed that genetic distances between strains were large enough to consider a reclassification into three distinct species, namely one containing phylotype I and III strains, termed R. pseudosolanacearum, another species consisting of phylotype II strains, termed R. solanacearum, and a third species encompassing strains from phylotype IV, termed R. syzygii (Remenant et al., 2010).
Instead of the categorization into phylotypes, the RSSC classification into five different races reflects its host range and pathogenicity (Buddenhagen et al., 1962; Denny, 2006). Race 1 strains can infect tobacco (Nicotiana tabacum), tomato, potato, eggplant (Solanum melongena), and diploid banana (Musa sp.), preferring temperatures of 35–37°C, just like race 2 strains that can infect triploid banana plants and Heliconia spp. and are responsible for the so-called “Moko disease” of banana. As such, race 1 and 2 strains are more widespread in tropical, subtropical, and warm temperate areas, whereas race 3 strains with a more limited host range can infect potatoes, tomatoes and weeds at moderate temperatures (27°C). Race 4 and race 5 strains can infect ginger (Zingiber officinale) and mulberry trees (Morus sp.) in China, respectively.
Ralstonia solanacearum owes its success to its ability to survive for an extended period in water, soil, and plant debris. To locate potential host roots, the pathogen uses chemotaxis, by which external signals, such as root exudates, influence bacterial motility and direction. Upon recognition, the bacteria stick to the root surface with polysaccharides, adhesin proteins, and type IV pili. From lateral root emergence sites or natural wounds, bacteria enter the roots and move inward toward the xylem vessels where they start growing, and are then transported systemically, either along with the sap flow or by twitching motility, eventually obstructing the flow of xylem sap (Planas-Marquès et al., 2020). The result is wilting and ultimately death of the host, subsequently releasing the pathogen into the environment, ready for the next infection.
Just like many other Gram-negative bacteria, the RSSC employs a type III secretion system (T3SS) that functions as a molecular syringe, pumping an array of virulence factors, designated T3SS effectors (T3Es) into the host cells to hijack the cellular signaling and manipulate the host’s immunity, with the pathogen’s unique infection strategy as a consequence. In this manner, the plant response is ineffective in preventing the disease, while the bacteria use its host as replicative niche and nutrient source.
1.2. The motive: The role of type III effectors in immunity and host susceptibility
Plants possess a multilayered immune system to deal with neighboring organisms that aim at using the plant as a host for nutrients and/or niche. As a first line of defense, conserved microbially associated molecular patterns (MAMPs), such as flagellin or lipopolysaccharides, are detected by specific pattern recognition receptors (PRRs) expressed on the surface of host membranes, resulting in an innate immunity response. These receptors are either plasma membrane-bound receptor-like kinases (RLKs) or receptor-like proteins (RLPs) that control a signaling cascade once the respective MAMPs are recognized. This innate immune response is referred to as MAMP-triggered immunity (MTI) or, in the case of pathogens, pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI). Recognition of microbial patterns results in different defense responses, such as callose deposition, oxidative burst, cell wall strengthening, and expression of pathogenesis-related proteins.
At the molecular level, various modes of actions have been attributed to T3Es, among which MTI pathway inhibition (Büttner, 2016), and nutrient acquisition (Xian et al., 2020). Many T3Es target PRRs to shut down host defense pathways. Multiple T3Es from the same pathogen can even target the same PRR. The phytopathogen Pseudomonas syringae pv. tomato (Pto) DC3000, for example, encodes four effectors (AvrPto, AvrPtoB, HopF2, and HopQ1) that all (in)directly target the same leucine-rich repeat (LRR)-RLK FLAGELLIN-SENSITIVE 2 (FLS2), responsible for detection of flagellin. Together, these T3Es block downstream the FLS2 signaling by (i) hindering interaction with its coreceptor BRASSINOSTEROID RECEPTOR-ASSOCIATED KINASE 1 (BAK1; AvrPto, AvrPtoB, and HopF2; Shan et al., 2008; Zhou et al., 2014), (ii) preventing phosphorylation of its downstream signaling component BOTRYTIS-INDUCED KINASE 1 (BIK1; AvrPto; Xiang et al., 2010), (iii) targeting FLS2 for degradation (AvrPto; Göhre et al., 2008), or (iv) suppressing FLS2 accumulation in a cytokinin-dependent fashion (HopQ1; Hann et al., 2014). AvrPto also associates with ELONGATION FACTOR TU RECEPTOR (EFR), the PRR recognizing the bacterial elongation factor Tu (EF-Tu), to block MTI (Xiang et al., 2008). Interestingly, expression of the Arabidopsis thaliana EFR gene in tomato and tobacco (Nicotiana benthamiana) enhances resistance to various pathogens, such as Pseudomonas, Agrobacterium, Xanthomonas, and Ralstonia (Lacombe et al., 2010). Thus, signaling pathways downstream of the EF-Tu recognition are probably conserved between plant species and are interesting targets for disease resistance engineering. Notably, flagellin does not represent a major defense elicitor in R. solanacearum cells (Pfund et al., 2004), since the presence of polymorphisms in R. solanacearum flg22 avoids detection by FLS receptors from multiple plants such as Arabidopsis and members of the Solanaceae family (Wei et al., 2020). This finding, further exemplifies the evolutionary arms race occurring at the level of MTI. Additionally, T3Es may alter plant metabolism to obtain nutrient sources. The R. solanacearum T3E RipI, for example, interacts with and promotes activity of host glutamate decarboxylases that catalyze the biosynthesis of gamma-aminobutyric acid, a plant non-proteinogenic amino acid that can be used as nutritional source by the pathogen (Xian et al., 2020).
Plants can also distinguish T3Es either directly or indirectly by means of resistance proteins, enhancing the immune response, designated the effector-triggered immunity (ETI). Many resistance proteins contain a nucleotide-binding site/leucine-rich repeat (NB-LRR) structure. In ETI, pathogen spread is prevented by a hypersensitivity response (HR), i.e., a rapid localized cell death induced by reactive oxygen species (ROS) at the infection site. In A. thaliana Wassilewskija-2 (Ws-2), the RESISTANCE TO RALSTONIA SOLANACEARUM 1 (RRS-1) and RESISTANCE TO PSEUDOMONAS SYRINGAE 4 (RPS-4) proteins work as a dual resistance protein complex that recognizes the T3E Ralstonia-injected protein P2 (RipP2) of R. solanacearum (Narusaka et al., 2009). Moreover, the paired immune receptors RRS-1/RPS4 expressed in tomato triggers resistance to R. solanacearum and Pto DC3000 (Narusaka et al., 2013). Subsequently, coevolution has allowed many pathogens to develop mechanisms to overcome this immunity, for instance through the production of new ETI-suppressing effectors (Block et al., 2014).
As single T3E mutants rarely show an altered phenotype, T3Es have been hypothesized to form a robust interconnected network, in which the effectors can be seen as nodes and the interaction between the effectors or their cellular partners as edges, providing the bacterium with a certain flexibility (Sanchez-Garrido et al., 2021). However, this network robustness is limited, because individual effectors can become essential when part of the network is perturbed (Ruano-Gallego et al., 2021). This “context-dependent effector essentiality” illustrates the codependency of effectors. Effector-effector interactions may be involved in shaping this interconnected network. These interactions can be either antagonistic, where effector function or recognition by the host is reduced or prohibited, or synergistic, where the interaction promotes effector activity (Urbanus et al., 2016). This extra layer of regulation, also referred to as metaeffector activity, provides the pathogen with additional strategies to establish host susceptibility (Martel et al., 2022). Therefore, investigation of T3Es, their interactions, and functions in planta in a context-dependent manner is an important research area to understand the coevolution of effector-host interactions and to discover potential plant targets for resistance breeding. Although extensively studied, adaptation of breeding strategies to prevent phytopathogen attacks are still lacking for the RSSC, mainly because of its worldwide geographic distribution and acclimatization, broad host range, ability to survive in soil and water for extended time periods, and killer potential.
Whereas some T3E functions have already been discovered for different phytopathogens, many are still unexplored. Although the RSSC encodes the highest number of T3Es of all bacterial phytopathogens reported to date, the functions of only a few Ralstonia effectors have been described, in contrast to the extensively characterized T3Es of other well-studied bacterial phytopathogens, such as P. syringae and Xanthomonas spp. (reviewed in Schwartz et al., 2015). The reasons for the lack of knowledge on T3Es from the RSSC are most probably the size of the type III effectome (the repertoire of T3Es from a microbe) and its extraordinary diverse host range (Hayward, 1991).
1.3. Beyond a reasonable doubt: Challenges to overcome by the characterization of Ralstonia T3Es
As aforementioned, the description of the type III effectome of phytopathogens is complicated by the high number of T3Es typically encoded by phytopathogens. Effectors are usually grouped into families, referring to vertically inherited effectors indicative of an ancestral horizontal gene acquisition. To date, 112 curated T3E families have been reported for 155 strains of the RSSC with an average of 46–71 T3Es per strain.1 This number is higher than that of other well-known bacterial phytopathogens, such as Xanthomonas and P. syringae that encode, on average, 30 T3Es per strain (Schwartz et al., 2015; http://www.pseudomonas-syringae.org/). Basically, the functional characterization of the effectome has been lagging because of the large number of strains in the RSSC, each with its unique set of T3Es. To simplify matters, an interesting route could be to investigate the T3Es shared between strains, or the “core” effectome, that possibly represent the T3Es required for the successful pathogen infection strategy. For example, eight T3Es have been shown to be shared across 84 strains from four different phylotypes (Sabbagh et al., 2019). Nonetheless, “unique” effectors identified in phytopathogenic bacteria might provide more insight into specific host adaptation mechanisms. For instance, 11 and 9 T3Es are unique to phylotype I and III and to phylotype II strains, respectively (Sabbagh et al., 2019).
Alternatively, instead of focusing on the similarities or differences in T3E composition between strains, T3Es shared between strains and able to infect the same host could be investigated. However, because the RSSC has a broad host range that is difficult to define, such an approach is often complex. Nevertheless, the genomes of 25 RSSC strains grouped in different pathovars based on their host specificity were compared (Cho et al., 2019) and analysis of the pan-genome orthologous group and T3E repertoires revealed that, in contrast to their nonpathogenic variants, the pathogenic strains of tomato, eggplant, and pepper (Capsicum annuum) shared 8, 7, and 34 T3Es, respectively, but also that the T3Es RipS3 and RipH3 were found only in tomato-pathogenic strains and RipAC exclusively found in eggplant-pathogenic strains (Cho et al., 2019). From a dataset of 102 T3E sequences from the RSSC, the genomes of 140 distinct R. solanacearum strains were scanned for the presence of T3E sequences and 29 T3Es were common and 15 T3Es unique between Ralstonia strains isolated from tomato and eggplant, indicating that, although many T3Es are shared, different effectors might be required for different hosts (Sabbagh et al., 2019), a concept generally proven for other model species as well (Lindeberg et al., 2009; Schwartz et al., 2015). As the host range of the RSSC is very diverse, additional studies on plant host-specific T3E repertoires are required to further unravel T3E-host interdependencies. Recent research focusing on T3Es of mammalian pathogens even show that distinct effector networks are required depending on the occupied environmental host niche and context (Chen et al., 2021; Ruano-Gallego et al., 2021). Furthermore, because R. solanacearum occupies different parts of the plant during colonization and disease progression (Planas-Marquès et al., 2020), an interesting approach could be to investigate effector expression in a temporal fashion.
These examples indicate that the effectome size, the large strain diversity, and the broad host range result in a vast number of effector-strain-host combinations, especially complicating T3E research for the RSSC. Focus on unique or core T3Es or on T3Es associated with certain hosts can help the researcher in selecting the T3E to investigate.
In this review, we will examine the different manners to typify a T3E and, additionally, the effectome, based on previously identified T3Es and in view of reported methodologies with the RSSC as a case study. Moreover, we will highlight the current challenges toward the full and functional characterization of the RSSC effectome.
2. Different forensic disciplines for the functional identification of a T3E
2.1. Laying the foundation of the investigation: T3E prediction
Different approaches to predict T3Es can be divided into genome-dependent and genome-independent methods: the former include various computational prediction programs by which new T3Es are expected based on available information on T3E sequences from well-known plant and animal microbes, whereas the latter rely not solely on available genetic information, but also enable the identification of effectors based on functional analysis.
The computational prediction methods use in silico information on the intrinsic genetic elements from previously identified T3Es, not necessarily of the same species, and protein features, such as (positional) amino acid representation or the presence of translocation signal sequences, etc. For instance, T3Es frequently display a high G/C content and are often acquired by horizontal gene transfer. Additionally, T3E expression is typically under the control of a master regulator, denoted by specific regions in the promoter. Specifically for the RSSC strains, a hrpII box motif (TTCGn16TTCG) in the promoter region (Cunnac et al., 2004), necessary for HrpB-dependent activation, is commonly used as a criterion for RSSC T3E prediction, alongside homology to known T3Es and the presence of N-terminal specific export patterns (Peeters et al., 2013). Approximately 75% of the identified T3Es from R. solanacearum strain GMI1000 contain such a hrpII box motif in the corresponding promoter region (Cunnac et al., 2004). Likewise, 102 T3E-encoding genes and 16 hypothetical T3E genes in 12 evolutionarily distant strains from the RSSC have been detected (Sabbagh et al., 2019). By means of information on intrinsic genetic properties, different online tools were built allowing users to “blast” a query to the T3E database or scan a genome from the RSSC for the presence of putative T3Es (https://iant.toulouse.inra.fr/bacteria/annotation/site/prj/T3Ev3/# -Tools). As a consequence, T3E sequences in a genome could be screened when a new Ralstonia genome had been sequenced. The tool was, for example, used to identify T3Es from four different strains of R. solanacearum isolated either from blueberries (Vaccinium), hybrid tea roses, rocktrumpet (Mandevilla), or African daisies (Osteospermum; Bocsanczy et al., 2022). Many other tools for T3E prediction as well as T3E databases are also available, but are not restricted to the RSSC (for instance, The EuroXanth DocuWiki: https://internet.myds.me/dokuwiki/doku.php?id=bacteria:t3e:software). Some examples of the T3E prediction tools listed are, among others, Effectidor, EffectiveT3, and T3SPs (Arnold et al., 2009; Yang et al., 2013; Wagner et al., 2022).
Although useful, one should be aware that genome-dependent approaches can still produce false positive or false negative predictions. Indeed, not all characterized T3Es carry a conserved N-terminal secretion signal and many display a high sequence diversity, complicating homology-based database searches. Noteworthy, predictions obtained from different annotation tools concur moderately, even more so, when the accuracy of translation initiation site (TIS) predictions is considered. These findings are consistent with recent proteogenomic endeavors, revealing multiple translation initiation evidence that occurs at nonannotated start sites, besides showing extensive translation outside the annotated protein-coding regions. Start codon plurality within the same coding sequence, which - in the case of in-frame start codons - may result in translation of N-terminal proteoforms (i.e., molecular protein forms of a single gene originating from alternative TIS selection; Willems et al., 2022). As a representative example of T3E proteoform expression, an N-terminally extended proteoform of the Salmonella enterica subspecies enterica serovar Typhimurium (S. Typhimurium) deubiquitinase T3E SseL has been reported (Figure 1; Ndah et al., 2017; Fijalkowski et al., 2022), meaning that the correct annotation of the expressed open reading frame (ORF) is of critical importance, because it serves as the starting point for further follow-up studies.
Figure 1. Retapamulin-assisted ribosome profiling (Ribo-RET) profile of S. typhimurium T3E SseL pointing to translation initiation at a newly identified upstream TIS. Bacterial ribosomes are specifically arrested at translation start sites by the antibiotic retapamulin. Here, translation initiation at a newly identified upstream TIS resulted in translation of an N-terminal extended SseL proteoform (Chromosome: 2392442–2393461). The annotated and newly discovered TISs are marked with a blue and orange arrowhead, respectively. Data from Fijalkowski et al. (2022).
Besides in silico prediction, tools based on functional characterization might be helpful in the search for new effectors. As the expression of many T3Es from the RSSC is regulated by the HrpB regulator, microarrays have been used to identify putative T3Es by comparing the expression of a wild type with a hrpB− mutant strain. In this manner, 143 HrpB-upregulated genes were detected in R. solanacearum GMI1000 by means of a microarray consisting of 5,074 of the predicted coding sequences, thereby extending the repertoire with 26 new candidate T3Es and exemplifying that gene expression can also be used for the identification of T3E genes (Occhialini et al., 2005).
Another genome-independent approach utilizes a truncated version of the adenylate cyclase domain from a calmodulin-dependent adenylate cyclase (CyaA’) of the cyclolysin toxin of Bordetella pertussis. This method requires, in contrast to genome-based T3E prediction, no prior knowledge of the genome or its encoded T3Es. Many T3Es from the RSSC have been discovered in this manner (Mukaihara and Tamura, 2009; Mukaihara et al., 2010) and a detailed protocol has been reported (Lonjon et al., 2018). In short, the CyaA’ reporter assay relies on the random bacterial genome insertion of an N-terminal-truncated cyaA’ gene flanked by transposon sequences (Figure 2). If by chance the cyaA’ gene is inserted in-frame into an ORF encoding (part of) a T3E and if the insertion still permits translocation into the plant cell, the translocated effector-cyaA’ fusion leads to the accumulation of cyclic adenosyl monophosphate (cAMP) in the plant. Hence, because intrinsic cAMP levels in plant cells are low, an increase in cAMP would hint at the functional translocation of an effector-cyaA’-fusion protein. The upstream region of the transposon can then be queried from these mutants and checked for the presence of a hrpII box, supporting the discovery of an effector. Other methods to monitor effector translocation have been reviewed by Braet et al. (2022).
Figure 2. CyaA’ translocation assay to discover T3Es. By means of transformation, the N-terminal part of the calmodulin-dependent adenylate cyclase gene (cyaA’) is randomly inserted into the genome of the bacterium. If genomic insertion disturbs effector translocation to the plant cell via the T3SS, no change in cyclic AMP (cAMP) occurs (left). However, if random insertion allows translocation to the plant cell, CyaA’ converts ATP to cAMP (right). The genome of bacteria that induced plant cAMP levels can then be queried for the presence of an effector by use of primers binding to the transposon elements flanking the cyaA’ gene (black).
2.2. Unmasking the suspects: T3E validation by in vitro secretion or in vivo translocation
Once effectors have been predicted, functional validation is required. Secretion or translocation of putative T3Es should always be verified before further in-depth analysis. In vitro secretion of T3E candidates is often investigated by the analysis of a translational fusion of the T3E with an epitope tag. Under T3E secretion-inducing conditions in bacterial cultures, secreted T3Es can be detected in the bacterial supernatant by immunoblotting (Lohou et al., 2014). In view of the time needed for cloning and/or endogenous tagging, this analysis is time consuming when several effectors would be analyzed simultaneously. For a more global assessment of T3E secretomes, proteomics approaches can be used. In this manner, 35 T3Es were identified, among which the previously unidentified T3E RipBJ in the supernatant of R. solanacearum strain GMI1000, when compared to secretome profiles of T3SS-defective mutants (Lonjon et al., 2016).
To make T3E prediction and validation more comprehensive, a combined transcriptomic and proteomic approach can be selected, as done for the symbiotic bacterium Bradyrhizobium vignae strain ORS3257. Here, the transcriptome of a wild-type strain was compared to a T3E-transcriptional regulator mutant sampled under either T3E-inducing or noninducing conditions whereafter the supernatant culture of these strains was compared by shotgun proteomics. Combined with an in silico approach, 36 putative T3Es were identified in at least one of the approaches (Busset et al., 2021). This method could also work for the RSSC by comparison of the transcriptome and proteome of the wild type and a transcriptional regulator (hrpB) mutant strain.
Instead of the secretion assessment in culture media, in vivo translocation to the host plant cell of putative T3Es can be analyzed, although such an analysis might be complicated by the frequent nonsynchronous nature of infection and the time-dependent differences of T3E translocation, hence limiting the sensitivity of such in vivo methods. Nevertheless, CyaA’ reporters can, for example, also be used for the validation of T3E translocation into plant cells (Sory and Cornelis, 1994). So, for instance, the translocation into tobacco (N. tabacum) was validated for RipBJ by measuring cAMP levels of GMI1000 compared with a T3SS-impaired (hrcV−) mutant strain carrying RipBJ-CyaA’ (Lonjon et al., 2016).
2.3. Gathering leads: Online tools to aid in T3E characterization
Once a T3E has been identified, metadata can be gathered with several freely available online tools to help decide on further experimental steps. For example, it can be helpful to know beforehand at which subcellular host localization the translocated T3E is expected or whether the T3E has any relatedness with other better-defined T3Es. The presence of certain protein domains (e.g., enzymatic or protein interaction motifs/domains), potentially identified through structural similarity with other proteins, may assist in designing experiments for the characterization of the biochemical function and activity of the T3E under study. Here, we list some online tools that can be used for these purposes (Table 1). Although many other online prediction tools exist, they mostly provide comparable outputs.
Some databases, such as InterPro or Uniprot, incorporate several different of the aforementioned prediction tools to provide an overview of the current knowledge on a certain entry, but they do not allow do novo sequence interpretations, albeit other tools (indicated by superscript 2 in Table 1) do look at intrinsic properties of the submitted sequence, whether or not supported by homology searches. Noteworthy, regardless of the program used, the output merely concerns predictions and requires experimental validation. For instance, application of the Phyre2 software to analyze the protein sequence of the R. solanacearum T3E RipAK provided information on putative functional domains and allowed the set-up of experiments by means of domain deletion mutants (Sun et al., 2017), resulting in the finding that the full-length RipAK is required for HR suppression in tobacco. Phyre2 was also used to investigate the T3E RipAY and revealed a high structural similarity with γ-glutamyl cyclotransferase (GGCT), as well as conserved catalytic key GGCT residues, leading to the discovery of the RipAY GGCT activity through mutational analyses of corresponding putative catalytic residues (Fujiwara et al., 2016).
An upcoming tool is AlphaFold and the more recent protein complex prediction algorithm AlphaFold-Multimer. AlphaFold uses a neural network to predict the 3D structure of individual protein chains based on information on the amino acid sequence, multiple sequence alignments, and homology to other proteins (Senior et al., 2020). Although currently it is still not possible to submit de novo protein sequences for 3D structure prediction in AlphaFold, this can be overcome by use of the source code available via Collabfold online. AlphaFold-Multimer was developed, extending AlphaFold to multiple chains as well (Evans et al., 2021), possibly allowing the mapping of effector-host protein interactions in the future. By means of DALI, a protein structure can be compared with protein structures in the Protein Data Bank (PDB) or with a species subset (such as Arabidopsis) in the AlphaFold database provided a protein structure has been predicted or defined and a PDB file of the protein is available (Holm, 2020). In this manner, it can become clear which other proteins share a similar 3D structure to use as a basis for further experiments.
2.4. Evidence collection: Experimental procedures to characterize a T3E
The method of choice for a functional characterization often depends on the field of interest, the background of the researcher, and the available resources. Nevertheless, different methods may offer complementary information to obtain a more complete view on T3E functioning (Figure 3). Whereas some methods, such as virulence assays, localization and/or protein structural analyses give no direct evidence of the effector’s function, they might provide a starting point for further characterization by other approaches, such as proteomics, biochemical activity assays, and metabolomics, among others (Figure 3).
Figure 3. Different experimental methods to characterize a T3E. Various methods can be applied in a complementary fashion to elucidate T3E functioning. Commonly used experimental methods are listed, but this list is not exhaustive.
To our knowledge, at present, only 14 T3Es from the RSSC have a known in planta function described and investigated with the aforementioned methods. In Table 2, the function(s) of these T3Es are summarized and as a representative example, in Table 3, the methods used to functionally characterize the T3E RipAB are listed. Other (extensively) characterized T3Es from the RSSC are listed together with the study methods in Supplementary Table 1. In this review, we will further discuss some of the commonly applied and upcoming state-of-the-art methods to functionally characterize T3Es.
2.4.1. Securing the crime scene: Subcellular localization
A reporter-tagged T3E is often expressed in planta to check full-length expression of the fusion construct and, additionally, to provide information on its subcellular localization. For phytopathogens, the go-to transient expression system is tobacco infiltration with an Agrobacterium tumefaciens strain carrying a binary vector for the selection and expression of a tagged effector gene. Other in planta expression systems might also be used, such as expression in transgenic roots obtained through A. rhizogenes transformation or transgenic A. tumefaciens-transformed A. thaliana (Deslandes et al., 2003; Tasset et al., 2010; Ron et al., 2014; Sun et al., 2017). Lastly, infection with a bacterial mutant strain encoding an endogenously fluorescently tagged T3E (by means of allelic replacement or recombination) or transformed with an alike T3E expression construct can provide insight into host effector localization in a physiological setting. In the case of fluorescent or luminescent reporter expression constructs, confocal microscopy or bioluminescence microscopy, potentially supported by the expression of subcellular markers, has been utilized to localize T3Es in different subcellular compartments (Denne et al., 2021). Whereas no direct information on T3E functioning is obtained this way, it can provide contextual information for the interpretation of the results acquired from other methods applied. RipE1, for example, localizes to the cytoplasm and nucleus of N. benthamiana cells (Nakano and Mukaihara, 2019). A yeast 2-hybrid (Y2H) screen with a cDNA library revealed an interaction of RipE1 with JASMONATE-ZIM DOMAIN (JAZ) proteins, known to be predominantly localized in the cell nucleus, thus matching the subcellular RipE1 localization (Nakano and Mukaihara, 2019).
Whereas native delivery of the T3E reporter by the phytopathogen is desirable over its (over)expression in the host, native translocation is limited, for instance, by the size or mechanical stability (i.e., the ease of protein unfolding under force) of the reporter tag used. This restriction was shown in the context of the green fluorescent protein (GFP) that displayed a comparable thermodynamic stability, but a higher mechanical stability compared to two Salmonella effectors (SptP and SopE2), implying that mechanical instability enhances secretion through the T3SS (LeBlanc et al., 2021; Braet et al., 2022). Some efforts have been made to overcome the passage drawbacks through the T3SS. For instance, the T3E RipP2 from R. solanacearum had been cloned to the C-terminal β strand of GFP (GFP11), instead of the commonly used full-length 11 β barrel strand-consisting GFP (Henry et al., 2017). The complementary part of GFP (GFP1-10) was expressed under a constitutive promoter in transformed Arabidopsis plants. After native delivery of the RipP2-GFP11 protein upon infection and complementation with GFP1-10 expressed in a ripP2-deletion mutant, a GFP signal was observed in the nucleus of Arabidopsis cells (Henry et al., 2017), a methodology commonly referred to as split-GFP. Other methods to study T3E subcellular localization, whether or not by native delivery, have recently been extensively reviewed (O'Boyle et al., 2018; Braet et al., 2022).
For the RSSC, the subcellular host localization of many T3Es has been reported. Recently, subcellular localization of a subset of 19 T3Es was assessed in tomato leaves and hairy roots, underlining the diverse destinations of different T3Es, including cell periphery, nucleus, tonoplast, and peroxisomes (Denne et al., 2021) and highlighting the importance of in planta localization studies, rather than relying solely on in silico localization predictions. Interestingly, subcellular localization prediction tools were accurate when the T3E was predicted to localize in the nucleus or plasma membrane, but less accurate for observed organellar or cytoskeleton localizations.
2.4.2. Conducting a primary survey: Virulence and effect on plant immunity
Many T3Es have been shown to be involved in avoidance or suppression of the innate immune response. This function can be analyzed through several assays, such as an infection assay, in which, in the case of R. solanacearum, wilting symptoms are compared between a wild type and a (multiple) T3E deletion mutant. For example, for R. solanacearum strain OE1-1, a multiple effector mutant lacking 42 T3Es was generated and the disease progression of infected tobacco plants measured in a virulence assay, revealing that some T3E families, such as the RipA T3E family members, had a higher impact on virulence than others, e.g., RipG or RipH T3E families (Lei et al., 2020). Instead of examining the effector virulence through effector deletion (loss-of-function), the involvement of single effectors can also be tested in a gain-of-function experiment in an effectorless polymutant of P. syringae DC3000 strain (Wei et al., 2018), as was done to study the effect of five core Ralstonia T3Es (RipAE, RipAQ, RipC1, RipU and RipW) on plant immune responses (Cong et al., 2022). Here, RipU was shown to reduce the ROS burst, and to upregulate MAPK cascades in tobacco leaves. In addition to disease progression measurement, bacterial enumeration and, thus, replication inside the host are also used as a proxy to estimate the influence of one or more T3Es on virulence (Deslandes et al., 2003; Sun et al., 2017; Yu et al., 2020).
The involvement of a T3E in virulence can also be investigated by overexpression of the T3E within plant cells and by determination of the effect on mechanisms known to be implicated in immunity. Most often the development of ensuing immunity-related effects is evaluated, for instance, by performing cell death assays in tobacco leaves with an overexpression construct after agro-infiltration, HR induction assessment, ion leakage assays, electrolyte leakage monitoring from leaves as a measure for the cell death severity, or ROS assays (Table 3; Sang et al., 2018, 2020; Nakano and Mukaihara, 2019; Yu et al., 2020). Suppression of PTI by a T3E can also be examined in a PTI inhibition assay, based on the principle that previous PTI induction by a nonpathogenic bacterium can dampen the localized ETI necrotic response elicited by a pathogenic bacterial strain. However, when the nonpathogenic bacterium delivers a T3E that suppresses PTI, an ETI necrotic response occurs where the pathogenic strain was infiltrated (Badel et al., 2013; Le Roux et al., 2015). Additionally, expression of PTI/ETI marker genes, such as FLG22-INDUCED RECEPTOR-LIKE KINASE 1 (FLK1) or PATHOGENESIS-RELATED GENE 1 (PR1), is also often monitored (for instance, by reverse transcription-quantitative polymerase chain reaction (qRT-PCR)). PTI marker expression levels are regularly compared between leaves treated with the flagellin-22 (flg22) peptide (a bacterial PAMP, positive control) and leaves expressing the T3E or a mutant version (Le Roux et al., 2015).
2.4.3. Checking for fingerprints: T3E expression profiling and the effect on the host transcriptome
The impact of T3Es on the host transcriptome is nowadays frequently investigated by means of RNA-sequencing (RNA-seq), although previously, microarrays had been used as well. In addition, the spatiotemporal expression of T3Es can also provide necessary information on the infection process itself and on the putative function of the T3E.
For instance, the spatiotemporal expression of T3Es examined by an RNA-seq approach in potato (de Pedro-Jové et al., 2021), retrieved RNA from potato plants infected with the cold-adapted R. solanacearum strain UY031 at three different infection stages, namely at an early, middle, and late stage, from the apoplast, xylem of asymptotic plants, and xylem of wilted plants, respectively. After normalization for bacterial counts, the T3E expression generally became more prominent at the middle and late stages, challenging the view that T3Es would mainly be important at early infection stages, although sensitivity at the earliest stages of infection might still be (too) limited. Nevertheless, some T3Es were found to have exceptional expression patterns, such as RipD and RipAD, both expressed at all infection stages (de Pedro-Jové et al., 2021).
The host can be monitored alongside the pathogen’s transcriptome with dual RNA-seq. This approach was used to observe both the T3E of R. solanacearum strain Rs-SY1 and the pepper gene expression during infection (Du et al., 2021). Compared to noninfected pepper plants, many pepper genes were differentially expressed 1, 3, 5, and 7 days after inoculation. Among others, at all stages, the genome ontology (GO) terms were enriched in the upregulated differentially expressed genes “defense response to bacterium,” “ethylene biosynthetic process” and “isoprenoid biosynthetic process” and in the downregulated differentially expressed genes “photosynthesis, light harvesting” and “response to abiotic stimulus.” Only genes involved in “sesquiterpenoid and triterpenoid biosynthesis” were gradually upregulated with disease progression. After mining the genome of R. solanacearum Rs-SY1, 84 T3E coding sequences were predicted, of which only seven had a higher expression in early infected pepper hypocotyls than that in the control, i.e., nutrient agar-grown R. solanacearum, hinting at the potential involvement of these T3Es at the onset of the disease (Du et al., 2021). Additionally, because only 7 out of 84 predicted T3Es were differentially expressed, the question was raised whether the other T3Es might have a function at later infection stages.
A transcriptomic approach can also investigate one or more T3Es under controlled expression instead of the broad T3E expression profiles during infection. The effect of T3E expression on transcriptomic profiles was studied in yeast cells ectopically expressing the T3E RipA5 (Popa et al., 2016). After measuring the genome-wide transcriptomic changes at several timepoints (2, 4, and 6 h) after induction of the RipA5 expression by means of DNA microarrays, 319 and 447 genes were discovered that were at least 2-fold induced and repressed, among which many nitrogen catabolite repression genes and ribosomal protein-encoding genes and genes involved in ribosome biogenesis, respectively. Intriguingly, this expression profile was fairly similar to that of rapamycin-treated yeast, a known inhibitor of the TOR COMPLEX 1 (TORC1) pathway. After some follow-up experiments, RipA5 could be shown to impact the TOR pathway in planta by acting as a TOR inhibitor in yeast and plant cells (Popa et al., 2016). A study on in planta function of the R. solanacearum T3E RipAB by RNA-seq (Zheng et al., 2019) revealed 417 differentially upregulated genes by comparison of three effector-expressing potato lines versus a control cultivar. A GO term enrichment analysis pointed to the involvement of RipAB in Ca2+ signaling (Zheng et al., 2019). In another experiment, RNA-seq carried out at different timepoints post inoculation (2 or 4 h) of N. tabacum cv. Xanthi leaves with wild-type R. solanacearum or with a ripAK deletion mutant identified regulated genes, including upregulation of photorespiration-related genes and downregulation of plant immunity-related genes, hinting at the involvement of RipAK in the modulation of early plant responses upon infection (Sun et al., 2017).
Transcription activator-like effectors (TALEs) also represent interesting targets for transcriptome-based analyses. TALEs represent a family of T3Es in the Xanthomonas genus, but also have one representative RSSC member, RipTAL1 (aka Brg11), that serve as transcription factors and activate the expression of certain virulence- or susceptibility-related genes by binding specific sites in target promoter regions. By differential RNA-seq of tomato plants infected for 24 h with X. euvesicatoria encoding either the R. solanacearum T3E RipTAL1 or RipTAL1 without its DNA-binding domain (DBD) under control of a constitutive promoter, the tomato ARGININE DECARBOXYLASE (ADC) genes SlADC1 and SlADC2 were identified as upregulated targets of wild-type RipTAL1 (Wu et al., 2019). Subsequently RipTAL1 was found to increase ADC activity and to boost plant polyamine levels, thereby potentially attenuating growth of R. solanacearum competitors.
2.4.4. Interrogating witnesses: Proteomics and interactomics approaches to study effector functions
For the inspection of differing protein expression profiles, a shotgun proteomics approach can be taken. Such a method was used to monitor proteome changes in early and later stages of P. syringae infection (Fan et al., 2019). Similarly, it can be applied to investigate the effect of a microbe on the posttranslational modifications of the host (Walley et al., 2018).
However, these methods do not provide information on the direct function of a T3E. As many T3Es exert their function by binding with host proteins or metabolites, knowledge on their interactions is thus of fundamental importance. In this section, we will focus on different methods that can be used to elucidate effector-host protein–protein interactions (EH-PPI), each with their advantages and disadvantages as recently extensively reviewed (Struk et al., 2019; De Meyer et al., 2020).
Although certainly valuable, functional protein microarrays have rarely been used to examine PPIs, mainly due to the complexity and the cost of the microarray design. In a protein microarray, individual proteins are immobilized on a microarray and a purified protein of interest is applied, resulting in a readable signal. In this manner, a proteome array probed with 15,000 proteins (56% proteome coverage) has been designed for Arabidopsis (Popescu et al., 2007; Manohar et al., 2014). Even though microarrays allow the detection of weak PPIs, the native context is missing, resulting in high false positive and negative rates because of the absence of chaperones (and thus also correct folding), posttranslational modifications, etc. Nevertheless, the interaction between a large set of human proteins and several type IV effectors from Legionella pneumophila (Yu et al., 2015, 2018) was successfully studied using protein microarrays, but, to our knowledge, proteome microarrays to investigate the interaction between T3Es and plant proteins are yet to be reported.
More frequently, Y2H-sequencing (Y2H-seq) and mass spectrometry (MS)-based approaches, such as affinity purification (AP) or biotin identification (BioID) are used to identify EH-PPIs. In Y2H, a bait protein (such as a T3E) is fused to a DBD and the prey (such as plant protein) to an activation domain (AD) or vice versa. Upon direct interaction, the DBD and AD are brought together, resulting in the reconstitution of a transcription factor that drives the reporter gene expression, often a gene responsible for the production of an essential amino acid. Conventional Y2H is frequently used for PPI validation, but also as screen by means of cDNA libraries, as in the case of Y2H-seq. In Y2H-seq, the bait (for instance, a T3E-expressing) yeast strain is supertransformed with a cDNA prey library, allowing the discovery of new PPIs. The cDNA library can be obtained by cDNA synthesis from the RNA of, for example, Arabidopsis plants infected with the phytopathogen of interest to reflect the expression profile more closely upon infection. cDNA libraries should preferably represent the expression profile upon infection, but might also be made after treatment of the host plant with biotic stress-related phytohormones, such as jasmonate or salicylic acid. The ensuing yeast colonies after supertransformation and selection are then pooled, the extracted plasmid DNA sequenced, and the resulting sequencing reads mapped to the host genome (in this case Arabidopsis; Erffelinck et al., 2018). The results of a large-scale Y2H-seq screening experiment for T3Es from two distinct vascular phytopathogens (R. solanacearum and X. campestris) against a cDNA library of Arabidopsis (González-Fuente et al., 2020) have been made public via an online interactive database2 and show that T3Es from distinct pathogens share some interactors, commonly referred to as plant effector hubs, whereas others are specific to the pathogen. Although Y2H-seq is a high-throughput and relatively cheap method, it is also linked to high rates of false positives and negatives due to the yeast background, protein folding artefacts, underrepresentation of full-length prey proteins encoded in cDNA libraries -especially large proteins-, and membrane protein targets.
For these reasons, in planta methods, such as (tandem) AP/IP or BioID coupled to MS-based methods are often preferred. The increasing sensitivity of mass spectrometers allows label-free quantification (LFQ), thereby facilitating experimental setups. Postmetabolic labeling approaches, such as tandem mass tag (TMT) labeling, are commonly used for multiplexed quantification. A frequently used IP-MS method is GFP trapping, in which GFP is fused to the protein of interest and antibodies against GFP precipitate the interaction complex. In this manner, interaction partners of several T3Es have been identified, such as the h-type thioredoxin targets (NbTRX-h9, NbTRX-h10, NbTRX-h11, and NbTRX-h14) of the R. solanacearum T3E RipAY in tobacco plants (Sang et al., 2018). By contrast, BioID relies on the use of a promiscuous biotin ligase that biotinylates primary amines of vicinal proteins after biotin application. When fused to the bait of interest, the proxeome (proteins in the bait proximity) can be identified by streptavidin-based purification. Over the years, different versions of the biotin ligase have been engineered: Escherichia coli BirA* (R118G; 35.3 kDa), Aquifex aeolicus BioID2 (27 kDa), and Bacillus subtilis BASU (28 kDa; Choi-Rhee et al., 2004; Kim et al., 2016; Ramanathan et al., 2018). However, these biotin-dependent proximity labeling tags show slow kinetics, require long biotin labeling times, and thrive at an optimal temperature of 37°C, making them rather unsuited for usage in plants. Nevertheless, despite poor labeling kinetics, the BirA* tag was successfully applied in rice (Oryza sativa) protoplasts to study the rice transcription factor OsFD2 (Lin et al., 2017) and later in transgenic Arabidopsis to investigate the proxeome of HopF2, a T3E from P. syringae (Khan et al., 2018). An AP-MS to complement the BioID results revealed an overlap of 58% (11/19) of the putative targets identified, although some previously reported interactors of HopF2 were missed with the BioID approach, probably because of the inherent BioID limitations, such as unavailability of free lysine residues in interacting preys, etc. (Khan et al., 2018). To overcome some of the challenges of the BirA* tag, such as a long labeling time, high biotin concentrations, and optimal temperature of 37°C, two new versions were created that enabled more efficient in planta proximity labeling. By a directed evolution-based approach, TurboID and its shorter version miniTurboID were established, allowing faster labeling and working temperatures around 30°C or less (Branon et al., 2018). Compared to the other proximity-labeling tags, TurboID was the most effective in different plant models, working at temperatures ranging from 22 to 28°C, and necessitating the addition of biotin for shorter periods of time (~2 h) for the efficient capture of plasma membrane interactomes (Arora et al., 2020). TurboID was, for example, used in tomato hairy root cultures to investigate the R. rhizogenes effector RolB, leading to the identification of TOPLESS and Novel Interactor of JAZ (NINJA) as direct interactors of RolB (Gryffroy et al., 2023). UltraID (19.7 kDa), a directed evolution variant of BioID2, has recently been developed, requiring even shorter labeling times, and showing less background biotinylation in the absence of exogenous biotin (Kubitz et al., 2022) in the case of mammalian cell culture, yeast, and bacteria. The use of ultraID in planta, however, has not been reported yet.
When these EH-PPI methods are used, a list of putative interactors is acquired. The read count or the peptide count/protein intensity (LFQ) and the expression ratios obtained from the Y2H-seq or MS-based methods can be indicative for the reliability (and proximity) of a certain hit. Although Y2H-seq, AP-MS, and BioID are complementary, the overlap of identified proteins is rather limited, mainly, as mentioned, because of the method inherent limitations and the fact that Y2H returns binary interactions, whereas AP-MS and BioID also can provide proximal interactions. Therefore, further validation of direct and indirect interactions remains necessary. For direct interactions, the highest scoring hits can be verified with conventional Y2H, co-immunoprecipitation (Co-IP), Förster Resonance Energy Transfer-Fluorescence-lifetime imaging microscopy (FRET-FLIM), bimolecular fluorescence complementation (BiFC), or split-luciferase, each with their advantages and disadvantages (Struk et al., 2019). Validation of indirect interactions requires alternative follow-up studies. Recently, a new validation technique for application in plants, designated knocksideways (KSP) is based on the heterodimerization of two domains, the human FK506-binding protein (FKBP) and the FKBP12-rapamycin-binding domain of mTOR (FRB) upon the addition of rapamycin (Winkler et al., 2021). FKBP is typically fused to the bait protein and FRB to an organellar marker of choice that differs from the subcellular localization of the bait protein. Upon rapamycin treatment, both domains reconstitute at the chosen organelle, unless the delocalization of the bait protein is hindered by interaction with a prey protein of interest. All three studied proteins are labeled with different fluorophores, so the (change in) localization can be monitored. An image analysis script was developed to quantify the interactions. Thus, KSP allows in planta localization of bait and prey fluorophores as well as validation of multiprotein interactions and quantification of the interaction.
2.4.5. Criminal profiling: Metabolomics
Besides their interaction with host proteins and/or alteration of host RNA levels, T3Es can also modify the host metabolism. Metabolomic approaches are frequently useful to describe microbe-host interactions. For example, an untargeted gas chromatography/MS metabolomics experiment revealed that 22 metabolites were enriched in the xylem upon infection of tomato plants with R. solanacearum, of which eight could be used as sole carbon or nitrogen source (Lowe-Power et al., 2018). Although in microbe-host metabolomic experiments, it is often challenging to determine which organism produces the identified metabolite, R. solanacearum was shown to be responsible for the production of putrescine, one of the 22 enriched metabolites found, because the pathogen harbors the corresponding biosynthesis and export genes. Putrescine was required for proliferation inside the xylem and an increased disease progression (Lowe-Power et al., 2018). Several attempts have been made to distinguish more easily the host from the microbe metabolites (Allwood et al., 2010; Pang et al., 2018). By labeling the metabolites from P. syringae DC3000 with stable heavy isotopes, Arabidopsis guard cell metabolites could be separated from microbial metabolites (Pang et al., 2018). This approach also allowed the separation of metabolic profiles at different timepoints after infection. Nevertheless, the role of T3Es in shaping the host metabolite landscape is still poorly understood. The use of single (or multiple) T3E deletion mutants could be insightful to determine their individual (or combined) roles in the generation of metabolomic changes during infection, but very few untargeted metabolomics have been done for single phytobacterial T3Es thus far. One example, however, is coumaroyl tyramine, a metabolic compound produced in the phenylpropanoid pathway that accumulates by the T3E WtsE from the maize (Zea mays) pathogen Pantoea stewartii spp. stewartia (Asselin et al., 2015). This compound was identified through LC–MS/MS by comparing the metabolite profiles of maize seedling leaves infected either with the wild-type or with the wtsE-defective strain. Instead of a metabolomic approach, many more targeted approaches have been utilized to determine the effect of a T3E on specific host metabolites, as, for instance, accumulation of salicylic acid, jasmonic acid and glutathione by the R. solanacearum T3E RipE1 was previously reported (Sang et al., 2020). After the identification of accumulated/depleted metabolites caused by T3E expression, a biochemical activity assay could be set up to pinpoint whether this T3E is directly (enzymatically) responsible for the accumulation/depletion of the compound.
2.4.6. Determining the modus operandi: Biochemical activity assays
Driven by predictions and gathered data, experiments can be set up to pinpoint the exact function of the T3E of interest in planta. In the case of predicted protein domains with a known enzymatic function, targeted biochemical assays can be carried out. For instance, a ChaC domain that displays a GGCT activity in yeast and mammalian cells, was also found in the RipAY protein sequence and prompted the analysis of its GGCT activity and its ability to degrade glutathione (Fujiwara et al., 2016), revealing that RipAY had a robust GGCT activity in the presence of yeast or plant thioredoxins. Later on, RipAY was shown to interact with h-type thioredoxins in planta (Sang et al., 2018). Similarly, the discovery of the RipP2 acylation (i.e., (auto)acetylation) activity was steered by its homology to the YopJ-like family of effectors with a wide range of activities, such as (de-)sumoylation (AvrRxv and AvrXv4 from X. campestris pv. vesicatoria), de-ubiquitination and acetylation (YopJ from Yersinia spp.; Mukherjee et al., 2006, 2007; Roden et al., 2007), leading to the determination of the exact acetylation targets of RipP2 (Tasset et al., 2010; Le Roux et al., 2015). A third example is the R. solanacearum T3E RipTPS that presents homology with proteins harboring trehalose-6-phosphate (T6P) synthase (TPS) activity and really catalyzes T6P synthesis in yeast cells (Poueymiro et al., 2014). Besides relying on homology of protein domains to design follow-up experiments, known or newly identified EH-PPIs can also inform the selection of biochemical activity assays. For example, a catalase activity assay revealed that RipAK interacts with host catalases and pyruvate decarboxylases, thereby inhibiting host catalase activity and suppressing the plant’s HR (Sun et al., 2017; Wang et al., 2021).
2.4.7. Generating a facial composite: T3E and T3E-host complexed protein structures
The examples mentioned above indicated that information on the protein structure can aid successful functional characterization of T3Es. Determination of the structure of a protein (complex) is often a challenging task (Fujiwara et al., 2020). Nevertheless, the protein structures of some T3Es have been elucidated. So, the structure of Pto DC3000 T3E AvrPto has been solved in solution by nucleic magnetic resonance (NMR) spectroscopy (Lee et al., 2004) and the crystal structure of AvrB was determined by single wavelength anomalous dispersion (SAD; Wulf et al., 2004). For the RSSC, only the crystal structures of RipP2 (apo or in complex) have been described (Xia et al., 2021). Complexed protein structures are also interesting for the characterization of the protein function. In one case, the crystal structure was determined for RipP2, a known acetyltransferase, in complex with its host cofactor inositol hexaphosphate (IP6), acetyl-coenzyme A (AcCoA), and the host resistance protein harboring a WRKY motif (RRS1-RWRKY; Zhang et al., 2017). The WRKYGQK motif of the RRS1-RWRKY substrate drives RipP2 to a lysine residue in the acetylation active site pocket. Additionally, IP6 enhanced AcCoA and WRKY binding with RipP2, providing novel insights into the regulation of the RipP2 activity. To our knowledge, other RSSC T3E protein structures are still to be reported.
A potential solution for the challenging protein structure determination has recently been developed. The program, designated AlphaFold, relies on artificial intelligence to predict protein structures and has a high prediction accuracy (Senior et al., 2020). The AlphaFold database is currently vastly expanding, including additional proteomes and catalogued protein structures for even more accurate predictions. At the publication time, the AlphaFold source code is accessible online (Collabfold) and can be used to predict the protein structure of any protein of interest, as was done for the Ralstonia T3E RipE1 that harbors a cysteine protease domain (Tsakiri et al., 2022). More recently, the algorithm AlphaFold-Multimer has been created that could be used for the modelling of (effector-host protein) interactions (Evans et al., 2021).
3. Conclusion
In this review, the RSSC was used as a case study to exemplify the current advances made using both conventional and state-of-the-art experimental techniques that have helped to functionally characterize T3Es. An overview of the methodologies used to elucidate the function of RipP2, one of the best described T3Es from the RSSC, is illustrated (Figure 4).
Figure 4. The route for the functional characterization of RipP2. Different laboratories have contributed to the elucidation of the putative function of the T3E RipP2, formally known as PopP2, by using the different methods discussed in this review. Starting from its homology with other known acetyltransferases, the biochemical activity of the T3E (acetylation of WRKY transcription factors, including RRS1) was found as well as the catalytic residues responsible for this function. Eventually, the apo crystal structure of RipP2 or RipP2 in complex with IP6, acetyl-CoA, and the WRKY domain from RRS1-R was determined. This figure provides an overview of the results found by Deslandes et al. (2003), Tasset et al. (2010), Le Roux et al. (2015), Xiou et al. (2015), Zhang et al. (2017), Xia et al. (2021), and Huh (2022).
Clearly, complementary discoveries add onto one another and accelerate the functional characterization of a protein. In the case of RipP2, the homology with acetyltransferases and the link with the RRS1-R protein in Arabidopsis has led to the eventual uncovering of the trans- and autoacetylation function of RipP2 (Figure 4; Deslandes et al., 2003; Tasset et al., 2010). This finding was mined for further identification of the RipP2 catalytic triad, for mapping the acetyl-receiving residues of host WRKY transcription factors, and elucidation of the apo crystal structure of RipP2 and RipP2 in complex with IP6, acetyl-CoA, and the WRKY domain from RRS1-R (Le Roux et al., 2015; Xiou et al., 2015; Zhang et al., 2017; Xia et al., 2021). In turn, available structural data might be helpful for homology modelling and functional characterization of a different effector protein.
Regarding the RSSC, still much research is required to fully functionally characterize its effectome. As mentioned, this process is hindered by the effectome size, the strain diversity, and the broad host range of the pathogen. Ample efforts were undertaken to optimize T3E prediction and to understand its effect on the HR response, for instance by tobacco infiltration of an effectorless mutant strain. Further, increasing knowledge on RSSC effector functioning could be obtained by performing complementary transcriptome, proteome, or metabolome studies in different host plants. The exact functions and host targets of various T3Es remain unsolved, although multiple T3Es from the RSSC have a (predicted) enzymatic function: RipE1 has a protease activity, targeting JAZ proteins; RipG T3E family members mimic eukaryotic E3 ligases to target unknown host proteins for degradation; RipAR, RipAW, and RipV2 have evolved novel E3 ubiquitin ligase domains and RipN is a nudix hydrolase which can hydrolyze a wide range of organic pyrophosphates (for a review, see Schreiber et al., 2021).
To elucidate their (enzymatic) targets, the need for specialized methods is high. For the identification of protease substrates, for example, a “degradomics” approach could be opted. This subfield of proteomics compares proteomes with and without proteolytic activity, but this approach, remains to be explored in the context of effector-plant interactions (Mooney et al., 2021). It would also be interesting to study whether T3Es behave differently, not only in different hosts, but also at various stages, or even in different environmental niches during infection. To this end, a time-course (multi)-omics experiment for different hosts could be envisaged, but undoubtedly in the most natural environment as possible for the different host infections. By native delivery of a tagged T3E upon infection, for example, the EH-PPI network exploited by the T3E could be investigated at different timepoints and/or in different hosts. As effectors seem to function in tightly connected networks, information on their potential cooperativity in the plant cell could also improve our understanding of bacterial pathogen-mediated disease development in plants (Sanchez-Garrido et al., 2021). Finally, collaborative efforts of the scientific community will allow the full functional description of the effectome and help understand the disease in general.
It is important to remain vigilant when characterizing a T3E, because new viewpoints and techniques could modify the perception of previous work. Techniques previously mainly used in mammalian cells, such as BioID, are applied to the plant field and have more recently successfully been used for the functional characterization of T3Es (Khan et al., 2018; González-Fuente et al., 2020). Technology has also progressed, such as the development of new BioID modules (e.g., TurboID and ultraID) that allow the use of these methods in broader and more natural (in vivo) settings (Branon et al., 2018; Kubitz et al., 2022). Still, scientists should always be on the lookout for generation of new experimental techniques that can overcome (part of) the shortcomings of others. The use of single-cell omics profiling, for example, is already being used for plant systems and a single-cell omics perspective of T3E research would certainly result in a more detailed comprehension of the disease (Clark et al., 2022), because different host cell types might respond differently to the pathogen or to its deployed effectome. As R. solanacearum moves through distinct cell types of the plant (Planas-Marquès et al., 2020), single-cell omics profiling could also be used in a similar way to look at T3E expression and translocation, thereby more accurately distinguishing early from late T3E functioning. Altogether, the insights gained into T3E functioning by complementary approaches (e.g., AP-MS and BioID) should always be explored and combined to obtain a more detailed understanding of T3E functioning.
Author contributions
JDR: conceptualization, investigation, figure preparations, and writing—review and editing. PVD: supervision, conceptualization, investigation, figure preparations, and writing—review and editing. SG: supervision, conceptualization, investigation, resources, and writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (PROPHECY grant agreement No 803972) from the Research Foundation—Flanders (project number G.0511.20N to PVD) and by the Special Research Fund – Concerted Research Actions (project number BOF18-GOA-013 to SG). JDR is a predoctoral fellow of the Research Foundation-Flanders (Strategic Basic Research fellowship no. 1S83919N).
Acknowledgments
We thank Martine De Cock for her help in preparing the manuscript. Figures 2–4 were created with BioRender.com.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1113442/full#supplementary-material
Footnotes
References
Allwood, J. W., Clarke, A., Goodacre, R., and Mur, L. A. J. (2010). Dual metabolomics: a novel approach to understanding plant–pathogen interactions. Phytochemistry 71, 590–597. doi: 10.1016/j.phytochem.2010.01.006
Altenhoff, A. M., Train, C.-M., Gilbert, K. J., Mediratta, I., Mendes de Farias, T., Moi, D., et al. (2021). OMA orthology in 2021: website overhaul, conserved isoforms, ancestral gene order and more. Nucleic Acids Res. 49, D373–D379. doi: 10.1093/nar/gkaa1007
Arnold, R., Brandmaier, S., Kleine, F., Tischler, P., Heinz, E., Behrens, S., et al. (2009). Sequence-based prediction of type III secreted proteins. PLoS Pathog. 5:e1000376. doi: 10.1371/journal.ppat.1000376
Arora, D., Abel, N. B., Liu, C., Van Damme, P., Yperman, K., Eeckhout, D., et al. (2020). Establishment of proximity-dependent biotinylation approaches in different plant model systems. Plant Cell 32, 3388–3407. doi: 10.1105/tpc.20.00235
Asselin, J. A. E., Lin, J., Perez-Quintero, A. L., Gentzel, I., Majerczak, D., Opiyo, S. O., et al. (2015). Perturbation of maize phenylpropanoid metabolism by an AvrE family type III effector from Pantoea stewartii. Plant Physiol. 167, 1117–1135. doi: 10.1104/pp.114.253120
Badel, J. L., Piquerez, S. J. M., Greenshields, D., Rallapalli, G., Fabro, G., Ishaque, N., et al. (2013). In planta effector competition assays detect Hyaloperonospora arabidopsidis effectors that contribute to virulence and localize to different plant subcellular compartments. Mol. Plant-Microbe Interact. 26, 745–757. doi: 10.1094/MPMI-06-12-0154-R
Bernhofer, M., Dallago, C., Karl, T., Satagopam, V., Heinzinger, M., Littmann, M., et al. (2021). PredictProtein - predicting protein structure and function for 29 years. Nucleic Acids Res. 49, W535–W540. doi: 10.1093/nar/gkab354
Block, A., Toruño, T. Y., Elowsky, C. G., Zhang, C., Steinbrenner, J., Beynon, J., et al. (2014). The Pseudomonas syringae type III effector HopD1 suppresses effector-triggered immunity, localizes to the endoplasmic reticulum, and targets the Arabidopsis transcription factor NTL9. New Phytol. 201, 1358–1370. doi: 10.1111/nph.12626
Bocsanczy, A. M., Bonants, P., van der Wolf, J., Bergsma-Vlami, M., and Norman, D. J. (2022). Identification of candidate type 3 effectors that determine host specificity associated with emerging Ralstonia pseudosolanacearum strains. Eur. J. Plant Pathol. 163, 35–50. doi: 10.1007/s10658-021-02455-w
Braet, J., Catteeuw, D., and Van Damme, P. (2022). Recent advancements in tracking bacterial effector protein translocation. Microorganisms 10:260. doi: 10.3390/microorganisms10020260
Branon, T. C., Bosch, J. A., Sanchez, A. D., Udeshi, N. D., Svinkina, T., Carr, S. A., et al. (2018). Efficient proximity labeling in living cells and organisms with TurboID. Nat. Biotechnol. 36, 880–887. doi: 10.1038/nbt.4201
Buddenhagen, I., Sequeira, L., and Kelman, A. (1962). Designation of races in Pseudomonas solanacearum. Phytopathology 52:726. (Abstract)
Busset, N., Gully, D., Teulet, A., Fardoux, J., Camuel, A., Cornu, D., et al. (2021). The type III effectome of the symbiotic Bradyrhizobium vignae strain ORS3257. Biomol. Ther. 11:1592. doi: 10.3390/biom11111592
Büttner, D. (2016). Behind the lines – actions of bacterial type III effector proteins in plant cells. FEMS Microbiol. Rev. 40, 894–937. doi: 10.1093/femsre/fuw026
Cao, P., Chen, J., Wang, R., Zhao, M., Zhang, S., An, Y., et al. (2022). A conserved type III effector RipB is recognized in tobacco and contributes to Ralstonia solanacearum virulence in susceptible host plants. Biochem. Biophys. Res. Commun. 631, 18–24. doi: 10.1016/J.BBRC.2022.09.062
Chen, D., Burford, W. B., Pham, G., Zhang, L., Alto, L. T., Ertelt, J. M., et al. (2021). Systematic reconstruction of an effector-gene network reveals determinants of Salmonella cellular and tissue tropism. Cell Host Microbe 29, 1531–1544.e9. doi: 10.1016/j.chom.2021.08.012
Cho, H., Song, E.-S., Heu, S., Baek, J., Lee, Y. K., Lee, S., et al. (2019). Prediction of host-specific genes by pan-genome analyses of the Korean Ralstonia solanacearum species complex. Front. Microbiol. 10:506. doi: 10.3389/fmicb.2019.00506
Choi, Y., and Chan, A. P. (2015). PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 31, 2745–2747. doi: 10.1093/bioinformatics/btv195
Choi-Rhee, E., Schulman, H., and Cronan, J. E. (2004). Promiscuous protein biotinylation by Escherichia coli biotin protein ligase. Protein Sci. 13, 3043–3050. doi: 10.1110/ps.04911804
Clark, N. M., Elmore, J. M., and Walley, J. W. (2022). To the proteome and beyond: advances in single-cell omics profiling for plant systems. Plant Physiol. 188, 726–737. doi: 10.1093/plphys/kiab429
Cong, S., Li, J. Z., Xiong, Z. Z., and Wei, H. L. (2022). Diverse interactions of five core type III effectors from Ralstonia solanacearum with plants. J. Genet. Genomics. doi: 10.1016/J.JGG.2022.04.018
Cunnac, S., Boucher, C., and Genin, S. (2004). Characterization of the cis-acting regulatory element controlling HrpB-mediated activation of the type III secretion system and effector genes in Ralstonia solanacearum. J. Bacteriol. 186, 2309–2318. doi: 10.1128/JB.186.8.2309-2318.2004
Dawson, N. L., Sillitoe, I., Lees, J. G., Lam, S. D., and Orengo, C. A. (2017). CATH-Gene3D: generation of the resource and its use in obtaining structural and functional annotations for protein sequences. Methods Mol. Biol. 1558, 79–110. doi: 10.1007/978-1-4939-6783-4_4
De Meyer, M., De Ryck, J., Goormachtig, S., and Van Damme, P. (2020). Keeping in touch with type-III secretion system effectors: mass spectrometry-based proteomics to study effector–host protein–protein interactions. Int. J. Mol. Sci. 21:6891. doi: 10.3390/ijms21186891
de Pedro-Jové, R., Puigvert, M., Sebastià, P., Macho, A. P., Monteiro, J. S., Coll, N. S., et al. (2021). Dynamic expression of Ralstonia solanacearum virulence factors and metabolism-controlling genes during plant infection. BMC Genomics 22:170. doi: 10.1186/s12864-021-07457-w
del Toro, N., Shrivastava, A., Ragueneau, E., Meldal, B., Combe, C., Barrera, E., et al. (2022). The IntAct database: efficient access to fine-grained molecular interaction data. Nucleic Acids Res. 50, D648–D653. doi: 10.1093/nar/gkab1006
Denne, N. L., Hiles, R. R., Kyrysyuk, O., Iyer-Pascuzzi, A. S., and Mitra, R. M. (2021). Ralstonia solanacearum effectors localize to diverse organelles in solanum hosts. Phytopathology 111, 2213–2226. doi: 10.1094/PHYTO-10-20-0483-R
Denny, T. (2006). “Plant pathogenic Ralstonia species” in Plant-associated bacteria. ed. S. S. Gnanamanickam (Dordrecht: Springer Netherlands), 573–644.
Deslandes, L., Olivier, J., Peeters, N., Feng, D. X., Khounlotham, M., Boucher, C., et al. (2003). Physical interaction between RRS1-R, a protein conferring resistance to bacterial wilt, and PopP2, a type III effector targeted to the plant nucleus. Proc. Natl. Acad. Sci. U. S. A. 100, 8024–8029. doi: 10.1073/pnas.1230660100
Dong, X., Lu, X., and Zhang, Z. (2015). BEAN 2.0: an integrated web resource for the identification and functional analysis of type III secreted effectors. Database 2015:bav64. doi: 10.1093/database/bav064
Du, H., Yang, J., Chen, B., Zhang, X., Xu, X., Wen, C., et al. (2021). Dual RNA-seq reveals the global transcriptome dynamics of Ralstonia solanacearum and pepper (Capsicum annuum) hypocotyls during bacterial wilt pathogenesis. Phytopathology 112, 630–642. doi: 10.1094/PHYTO-01-21-0032-R
Erffelinck, M.-L., Ribeiro, B., Perassolo, M., Pauwels, L., Pollier, J., Storme, V., et al. (2018). A user-friendly platform for yeast two-hybrid library screening using next generation sequencing. PLoS One 13:e0201270. doi: 10.1371/journal.pone.0201270
European and Mediterranean Plant Protection Organization (2023). EPPO Global Database [Online]. Paris: European Union. Available at: https://gd.eppo.int (Accessed January 25, 2023).
Evans, R., O’Neill, M., Pritzel, A., Antropova, N., Senior, A., Green, T., et al. (2021). Protein complex prediction with AlphaFold-Multimer. bioRxiv 2021:10.04.463034. doi: 10.1101/2021.10.04.463034
Fan, K.-T., Wang, K.-H., Chang, W.-H., Yang, J.-C., Yeh, C.-F., Cheng, K.-T., et al. (2019). Application of data-independent acquisition approach to study the proteome change from early to later phases of tomato pathogenesis responses. Int. J. Mol. Sci. 20:863. doi: 10.3390/ijms20040863
Fijalkowski, I., Willems, P., Jonckheere, V., Simoens, L., and Van Damme, P. (2022). Hidden in plain sight: challenges in proteomics detection of small ORF-encoded polypeptides. microLife 3:uqac005. doi: 10.1093/femsml/uqac005
Fujiwara, S., Ikejiri, A., Tanaka, N., and Tabuchi, M. (2020). Characterization of the mechanism of thioredoxin-dependent activation of γ-glutamylcyclotransferase, RipAY, from Ralstonia solanacearum. Biochem. Biophys. Res. Commun. 523, 759–765. doi: 10.1016/j.bbrc.2019.12.092
Fujiwara, S., Kawazoe, T., Ohnishi, K., Kitagawa, T., Popa, C., Valls, M., et al. (2016). RipAY, a plant pathogen effector protein, exhibits robust γ-glutamyl cyclotransferase activity when stimulated by eukaryotic thioredoxins. J. Biol. Chem. 291, 6813–6830. doi: 10.1074/jbc.M115.678953
Göhre, V., Spallek, T., Häweker, H., Mersmann, S., Mentzel, T., Boller, T., et al. (2008). Plant pattern-recognition receptor FLS2 is directed for degradation by the bacterial ubiquitin ligase AvrPtoB. Curr. Biol. 18, 1824–1832. doi: 10.1016/j.cub.2008.10.063
González-Fuente, M., Carrère, S., Monachello, D., Marsella, B. G., Cazalé, A.-C., Zischek, C., et al. (2020). EffectorK, a comprehensive resource to mine for Ralstonia, Xanthomonas, and other published effector interactors in the Arabidopsis proteome. Mol. Plant Pathol. 21, 1257–1270. doi: 10.1111/mpp.12965
Gryffroy, L., Ceulemans, E., Pérez, N. M., Venegas-Molina, J., Jaramillo-Madrid, A. C., Rodrigues, S. D., et al. (2023). Rhizogenic Agrobacterium protein RolB interacts with the TOPLESS repressor proteins to reprogram plant immunity and development. Proc. Natl. Acad. Sci. 120:e2210300120. doi: 10.1073/PNAS.2210300120
Guidot, A., Prior, P., Schoenfeld, J., Carrère, S., Genin, S., and Boucher, C. (2007). Genomic structure and phylogeny of the plant pathogen Ralstonia solanacearum inferred from gene distribution analysis. J. Bacteriol. 189, 377–387. doi: 10.1111/mpp.12965
Hann, D. R., Domínguez-Ferreras, A., Motyka, V., Dobrev, P. I., Schornack, S., Jehle, A., et al. (2014). The pseudomonas type III effector HopQ1 activates cytokinin signaling and interferes with plant innate immunity. New Phytol. 201, 585–598. doi: 10.1111/nph.12544
Hayward, A. C. (1991). Biology and epidemiology of bacterial wilt caused by Pseudomonas solanacearum. Annu. Rev. Phytopathol. 29, 65–87. doi: 10.1146/annurev.py.29.090191.000433
Henry, E., Toruño, T. Y., Jauneau, A., Deslandes, L., and Coaker, G. (2017). Direct and indirect visualization of bacterial effector delivery into diverse plant cell types during infection. Plant Cell 29, 1555–1570. doi: 10.1105/tpc.17.00027
Holm, L. (2020). DALI and the persistence of protein shape. Protein Sci. 29, 128–140. doi: 10.1002/pro.3749
Huerta-Cepas, J., Szklarczyk, D., Forslund, K., Cook, H., Heller, D., Walter, M. C., et al. (2016). eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, D286–D293. doi: 10.1093/nar/gkv1248
Huh, S. U. (2022). PopP2 interacts with PAD4 in an acetyltransferase activity-dependent manner and affects plant immunity 16. doi: 10.1080/15592324.2021.2017631,
Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. E. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053
Khan, M., Youn, J.-Y., Gingras, A.-C., Subramaniam, R., and Desveaux, D. (2018). In planta proximity dependent biotin identification (BioID). Sci. Rep. 8:9212. doi: 10.1038/s41598-018-27500-3
Kim, D. I., Jensen, S. C., Noble, K. A., Kc, B., Roux, K. H., Motamedchaboki, K., et al. (2016). An improved smaller biotin ligase for BioID proximity labeling. Mol. Biol. Cell 27, 1188–1196. doi: 10.1091/mbc.E15-12-0844
Kubitz, L., Bitsch, S., Zhao, X., Schmitt, K., Deweid, L., Roehrig, A., et al. (2022). Engineering of ultraID, a compact and hyperactive enzyme for proximity-dependent biotinylation in living cells. Commun. Biol. 5:657. doi: 10.1038/s42003-022-03604-5
Lacombe, S., Rougon-Cardoso, A., Sherwood, E., Peeters, N., Dahlbeck, D., van Esse, H. P., et al. (2010). Interfamily transfer of a plant pattern-recognition receptor confers broad-spectrum bacterial resistance. Nat. Biotechnol. 28, 365–369. doi: 10.1038/nbt.1613
Le Roux, C., Huet, G., Jauneau, A., Camborde, L., Trémousaygue, D., Kraut, A., et al. (2015). A receptor pair with an integrated decoy converts pathogen disabling of transcription factors to immunity. Cells 161, 1074–1088. doi: 10.1016/j.cell.2015.04.025
LeBlanc, M.-A., Fink, M. R., Perkins, T. T., and Sousa, M. C. (2021). Type III secretion system effector proteins are mechanically labile. Proc. Natl. Acad. Sci. U. S. A. 118:e2019566118. doi: 10.1073/pnas.2019566118
Lee, C. C., Wood, M. D., Ng, K., Andersen, C. B., Liu, Y., Luginbühl, P., et al. (2004). Crystal structure of the type III effector AvrB from Pseudomonas syringae. Structure 12, 487–494. doi: 10.1016/j.str.2004.02.013
Lei, N., Chen, L., Kiba, A., Hikichi, Y., Zhang, Y., and Ohnishi, K. (2020). Super-multiple deletion analysis of type III effectors in Ralstonia solanacearum OE1-1 for full virulence toward host plants. Front. Microbiol. 11:1683. doi: 10.3389/fmicb.2020.01683
Lin, Q., Zhou, Z., Luo, W., Fang, M., Li, M., and Li, H. (2017). Screening of proximal and interacting proteins in rice protoplasts by proximity-dependent biotinylation. Front. Plant Sci. 8:749. doi: 10.3389/fpls.2017.00749
Lindeberg, M., Cunnac, S., and Collmer, A. (2009). The evolution of Pseudomonas syringae host specificity and type III effector repertoires. Mol. Plant Pathol. 10, 767–775. doi: 10.1111/j.1364-3703.2009.00587.x
Lohou, D., Turner, M., Lonjon, F., Cazalé, A.-C., Peeters, N., Genin, S., et al. (2014). HpaP modulates type III effector secretion in Ralstonia solanacearum and harbours a substrate specificity switch domain essential for virulence. Mol. Plant Pathol. 15, 601–614. doi: 10.1111/mpp.12119
Lonjon, F., Peeters, N., Genin, S., and Vailleau, F. (2018). In vitro and in vivo secretion/translocation assays to identify novel Ralstonia solanacearum type 3 effectors. Methods Mol. Biol. 1734, 209–222. doi: 10.1007/978-1-4939-7604-1_17
Lonjon, F., Turner, M., Henry, C., Rengel, D., Lohou, D., van de Kerkhove, Q., et al. (2016). Comparative secretome analysis of Ralstonia solanacearum type 3 secretion-associated mutants reveals a fine control of effector delivery, essential for bacterial pathogenicity. Mol. Cell. Proteomics 15, 598–613. doi: 10.1074/mcp.M115.051078
Lowe-Power, T. M., Hendrich, C. G., von Roepenack-Lahaye, E., Li, B., Wu, D., Mitra, R., et al. (2018). Metabolomics of tomato xylem sap during bacterial wilt reveals Ralstonia solanacearum produces abundant putrescine, a metabolite that accelerates wilt disease. Environ. Microbiol. 20, 1330–1349. doi: 10.1111/1462-2920.14020
Manohar, M., Tian, M., Moreau, M., Park, S.-W., Choi, H. W., Fei, Z., et al. (2014). Identification of multiple salicylic acid-binding proteins using two high throughput screens. Front. Plant Sci. 5:777. doi: 10.3389/fpls.2014.00777
Mansfield, J., Genin, S., Magori, S., Citovsky, V., Sriariyanum, M., Ronald, P., et al. (2012). Top 10 plant pathogenic bacteria in molecular plant pathology. Mol. Plant Pathol. 13, 614–629. doi: 10.1111/j.1364-3703.2012.00804.x
Marchler-Bauer, A., Derbyshire, M. K., Gonzales, N. R., Lu, S., Chitsaz, F., Geer, L. Y., et al. (2015). CDD: NCBI's conserved domain database. Nucleic Acids Res. 43, D222–D226. doi: 10.1093/nar/gku1221
Martel, A., Laflamme, B., Breit-McNally, C., Wang, P., Lonjon, F., Desveaux, D., et al. (2022). Metaeffector interactions modulate the type III effector-triggered immunity load of Pseudomonas syringae. PLOS Pathogens 18:e1010541. doi: 10.1371/JOURNAL.PPAT.1010541
Mitchell, A. L., Attwood, T. K., Babbitt, P. C., Blum, M., Bork, P., Bridge, A., et al. (2019). InterPro in 2019: improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 47, D351–D360. doi: 10.1093/nar/gky1100
Mooney, B. C., Mantz, M., Graciet, E., and Huesgen, P. F. (2021). Cutting the line: manipulation of plant immunity by bacterial type III effector proteases10.1093/jxb/erab095. J. Exp. Bot. 72, 3395–3409. doi: 10.1093/jxb/erab095
Mukaihara, T., and Tamura, N. (2009). Identification of novel Ralstonia solanacearum type III effector proteins through translocation analysis of hrpB-regulated gene products. Microbiology 155, 2235–2244. doi: 10.1099/mic.0.027763-0
Mukaihara, T., Tamura, N., and Iwabuchi, M. (2010). Genome-wide identification of a large repertoire of Ralstonia solanacearum type III effector proteins by a new functional screen. Mol. Plant-Microbe Interact. 23, 251–262. doi: 10.1094/MPMI-23-3-0251
Mukherjee, S., Hao, Y.-H., and Orth, K. (2007). A newly discovered post-translational modification – the acetylation of serine and threonine residues. Trends Biochem. Sci. 32, 210–216. doi: 10.1016/j.tibs.2007.03.007
Mukherjee, S., Keitany, G., Li, Y., Wang, Y., Ball, H. L., Goldsmith, E. J., et al. (2006). Yersinia YopJ acetylates and inhibits kinase activation by blocking phosphorylation. Science 312, 1211–1214. doi: 10.1126/science.1126867
Nakano, M., and Mukaihara, T. (2019). Comprehensive identification of PTI suppressors in type III effector repertoire reveals that Ralstonia solanacearum activates jasmonate signaling at two different steps. Int. J. Mol. Sci. 20:5992. doi: 10.3390/ijms20235992
Nakano, M., Ichinose, Y., and Mukaihara, T. (2021). Ralstonia solanacearum Type III Effector RipAC Targets SGT1 to Suppress Effector-Triggered Immunity. Plant Cell Physiol. 61, 2067–2076. doi: 10.1093/PCP/PCAA122
Narusaka, M., Kubo, Y., Hatakeyama, K., Imamura, J., Ezura, H., Nanasato, Y., et al. (2013). Interfamily transfer of dual NB-LRR genes confers resistance to multiple pathogens. PLoS One 8:e55954. doi: 10.1371/journal.pone.0055954
Narusaka, M., Shirasu, K., Noutoshi, Y., Kubo, Y., Shiraishi, T., Iwabuchi, M., et al. (2009). RRS1 and RPS4 provide a dual resistance-gene system against fungal and bacterial pathogens. Plant J. 60, 218–226. doi: 10.1111/j.1365-313X.2009.03949.x
Ndah, E., Jonckheere, V., Giess, A., Valen, E., Menschaert, G., and Van Damme, P. (2017). REPARATION: ribosome profiling assisted (re-)annotation of bacterial genomes. Nucleic Acids Res. 45:e168. doi: 10.1093/nar/gkx758
O'Boyle, N., Connolly, J. P. R., and Roe, A. J. (2018). Tracking elusive cargo: illuminating spatio-temporal type 3 effector protein dynamics using reporters. Cell. Microbiol. 20:e12797. doi: 10.1111/cmi.12797
Occhialini, A., Cunnac, S., Reymond, N., Genin, S., and Boucher, C. (2005). Genome-wide analysis of gene expression in Ralstonia solanacearum reveals that the hrpB gene acts as a regulatory switch controlling multiple virulence pathways. Mol. Plant-Microbe Interact. 18, 938–949. doi: 10.1094/MPMI-18-0938
Pagni, M., Ioannidis, V., Cerutti, L., Zahn-Zabal, M., Jongeneel, C. V., Hau, J., et al. (2007). MyHits: improvements to an interactive resource for analyzing protein sequences. Nucleic Acids Res. 35, W433–W437. doi: 10.1093/nar/gkm352
Pang, Q., Zhang, T., Wang, Y., Kong, W., Guan, Q., Yan, X., et al. (2018). Metabolomics of early stage plant cell–microbe interaction using stable isotope labeling. Front. Plant Sci. 9:760. doi: 10.3389/fpls.2018.00760
Peeters, N., Carrère, S., Anisimova, M., Plener, L., Cazalé, A.-C., and Genin, S. (2013). Repertoire, unified nomenclature and evolution of the type III effector gene set in the Ralstonia solanacearum species complex. BMC Genomics 14:859. doi: 10.1186/1471-2164-14-859
Penel, S., Arigon, A.-M., Dufayard, J.-F., Sertier, A.-S., Daubin, V., Duret, L., et al. (2009). Databases of homologous gene families for comparative genomics. BMC Bioinform. 10:S3. doi: 10.1186/1471-2105-10-S6-S3
Pfund, C., Tans-Kersten, J., Dunning, F. M., Alonso, J. M., Ecker, J. R., Allen, C., et al. (2004). Flagellin is not a major defense elicitor in Ralstonia solanacearum cells or extracts applied to Arabidopsis thaliana. Mol. Plant-Microbe Interact. 17, 696–706. doi: 10.1094/MPMI.2004.17.6.696
Pieper, U., Webb, B. M., Barkan, D. T., Schneidman-Duhovny, D., Schlessinger, A., Braberg, H., et al. (2011). ModBase, a database of annotated comparative protein structure models, and associated resources. Nucleic Acids Res. 39, D465–D474. doi: 10.1093/nar/gkq1091
Planas-Marquès, M., Kressin, J. P., Kashyap, A., Panthee, D. R., Louws, F. J., Coll, N. S., et al. (2020). Four bottlenecks restrict colonization and invasion by the pathogen Ralstonia solanacearum in resistant tomato. J. Exp. Bot. 71, 2157–2171. doi: 10.1093/jxb/erz562
Popa, C., Li, L., Gil, S., Tatjer, L., Hashii, K., Tabuchi, M., et al. (2016). The effector AWR5 from the plant pathogen Ralstonia solanacearum is an inhibitor of the TOR signalling pathway. Sci. Rep. 6:27058. doi: 10.1038/srep27058
Popescu, S. C., Popescu, G. V., Bachan, S., Zhang, Z., Seay, M., Gerstein, M., et al. (2007). Differential binding of calmodulin-related proteins to their targets revealed through high-density Arabidopsis protein microarrays. Proc. Natl. Acad. Sci. U. S. A. 104, 4730–4735. doi: 10.1073/pnas.061161510
Poueymiro, M., Cazalé, A. C., François, J. M., Parrou, J. L., Peeters, N., and Genin, S. (2014). A Ralstonia solanacearum type III effector directs the production of the plant signal metabolite trehalose-6-phosphate. MBio 5:e02065. doi: 10.1128/mBio.02065-14
Qi, P., Huang, M., Hu, X., Zhang, Y., Wang, Y., Li, P., et al. (2022). A Ralstonia solanacearum effector targets TGA transcription factors to subvert salicylic acid signaling. Plant Cell 34, 1666–1683. doi: 10.1093/plcell/koac015
Ramanathan, M., Majzoub, K., Rao, D. S., Neela, P. H., Zarnegar, B. J., Mondal, S., et al. (2018). RNA–protein interaction detection in living cells. Nat. Methods 15, 207–212. doi: 10.1038/nmeth.4601
Remenant, B., Coupat-Goutaland, B., Guidot, A., Cellier, G., Wicker, E., Allen, C., et al. (2010). Genomes of three tomato pathogens within the Ralstonia solanacearum species complex reveal significant evolutionary divergence. BMC Genomics 11:379. doi: 10.1186/1471-2164-11-379
Roden, J., Eardley, L., Hotson, A., Cao, Y., and Mudgett, M. B. (2007). Characterization of the Xanthomonas AvrXv4 effector, a SUMO protease translocated into plant cells. Mol. Plant-Microbe Interact. 17, 633–643. doi: 10.1094/MPMI.2004.17.6.633
Ron, M., Kajala, K., Pauluzzi, G., Wang, D., Reynoso, M. A., Zumstein, K., et al. (2014). Hairy root transformation using Agrobacterium rhizogenes as a tool for exploring cell type-specific gene expression and function using tomato as a model. Plant Physiol. 166, 455–469. doi: 10.1104/pp.114.239392
Roy, A., Kucukural, A., and Zhang, Y. (2010). I-TASSER: a unified platform for automated protein structure and function prediction. Nat. Protoc. 5, 725–738. doi: 10.1038/nprot.2010.5
Ruano-Gallego, D., Sanchez-Garrido, J., Kozik, Z., Núñez-Berrueco, E., Cepeda-Molero, M., Mullineaux-Sanders, C., et al. (2021). Type III secretion system effectors form robust and flexible intracellular virulence networks. Science 371:eabc9531. doi: 10.1126/science.abc9531
Sabbagh, C. R. R., Carrere, S., Lonjon, F., Vailleau, F., Macho, A. P., Genin, S., et al. (2019). Pangenomic type III effector database of the plant pathogenic Ralstonia spp. PeerJ 7:e7346. doi: 10.7717/peerj.7346
Sanchez-Garrido, J., Ruano-Gallego, D., Choudhary, J. S., and Frankel, G. (2021). The type III secretion system effector network hypothesis. Trends Microbiol. 30, 524–533. doi: 10.1016/j.tim.2021.10.007
Sang, Y., Wang, Y., Ni, H., Cazalé, A.-C., She, Y.-M., Peeters, N., et al. (2018). The Ralstonia solanacearum type III effector RipAY targets plant redox regulators to suppress immune responses. Mol. Plant Pathol. 19, 129–142. doi: 10.1111/mpp.12504
Sang, Y., Yu, W., Zhuang, H., Wei, Y., Derevnina, L., Yu, G., et al. (2020). Intra-strain elicitation and suppression of plant immunity by Ralstonia solanacearum type-III effectors in Nicotiana benthamiana. Plant Commun. 1:100025. doi: 10.1016/j.xplc.2020.100025
Schreiber, K. J., Chau-Ly, I. J., and Lewis, J. D. (2021). What the wild things do: mechanisms of plant host manipulation by bacterial type III-secreted effector proteins. Microorganisms 9:1029. doi: 10.3390/microorganisms9051029
Schwartz, A. R., Potnis, N., Timilsina, S., Wilson, M., Patané, J., Martins, J., et al. (2015). Phylogenomics of Xanthomonas field strains infecting pepper and tomato reveals diversity in effector repertoires and identifies determinants of host specificity. Front. Microbiol. 6:535. doi: 10.3389/fmicb.2015.00535
Senior, A. W., Evans, R., Jumper, J., Kirkpatrick, J., Sifre, L., Green, T., et al. (2020). Improved protein structure prediction using potentials from deep learning. Nature 577, 706–710. doi: 10.1038/s41586-019-1923-7
Shan, L., He, P., Li, J., Heese, A., Peck, S. C., Nürnberger, T., et al. (2008). Bacterial effectors target the common signaling partner BAK1 to disrupt multiple MAMP receptor-signaling complexes and impede plant immunity. Cell Host Microbe 4, 17–27. doi: 10.1016/j.chom.2008.05.017
Sory, M.-P., and Cornelis, G. R. (1994). Translocation of a hybrid YopE-adenylate cyclase from Yersinia enterocolitica into HeLa cells. Mol. Microbiol. 14, 583–594. doi: 10.1111/j.1365-2958.1994.tb02191.x
Sperschneider, J., Catanzariti, A.-M., DeBoer, K., Petre, B., Gardiner, D. M., Singh, K. B., et al. (2017). LOCALIZER: subcellular localization prediction of both plant and effector proteins in the plant cell. Sci. Rep. 7:44598. doi: 10.1038/srep44598
Struk, S., Jacobs, A., Sánchez Martín-Fontecha, E., Gevaert, K., Cubas, P., and Goormachtig, S. (2019). Exploring the protein–protein interaction landscape in plants. Plant Cell Environ. 42, 387–409. doi: 10.1111/pce.13433
Sun, Y., Li, P., Deng, M., Shen, D., Dai, G., Yao, N., et al. (2017). The Ralstonia solanacearum effector RipAK suppresses plant hypersensitive response by inhibiting the activity of host catalases. Cell. Microbiol. 19:e12736. doi: 10.1111/cmi.12736
Sun, Y., Li, P., Shen, D., Wei, Q., He, J., and Lu, Y. (2019). The Ralstonia solanacearum effector RipN suppresses plant PAMP-triggered immunity, localizes to the endoplasmic reticulum and nucleus, and alters the NADH/NAD+ ratio in Arabidopsis. Mol. Plant Pathol. 20, 533–546. doi: 10.1111/mpp.12773
Sun, T., Wu, W., Wu, H., Rou, W., Zhou, Y., Zhuo, T., et al. (2020). Ralstonia solanacearum elicitor RipX Induces Defense Reaction by Suppressing the Mitochondrial atpA Gene in Host Plant. Int. J. Mol. Sci. 21:2000. doi: 10.3390/IJMS21062000
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Tasset, C., Bernoux, M., Jauneau, A., Pouzet, C., Brière, C., Kieffer-Jacquinod, S., et al. (2010). Autoacetylation of the Ralstonia solanacearum effector PopP2 targets a lysine residue essential for RRS1-R-mediated immunity in Arabidopsis. PLoS Pathog. 6:e1001202. doi: 10.1371/journal.ppat.1010368
Tsakiri, D., Kotsaridis, K., Michalopoulou, V. A., Kokkinidis, M., and Sarris, P. F. (2022). Ralstonia solanacearum core effector RipE1 interacts and cleaves the Arabidopsis exocyst component Exo70B1. BioRxiv 2022:31.506019. doi: 10.1101/2022.08.31.506019
Urbanus, M. L., Quaile, A. T., Stogios, P. J., Morar, M., Rao, C., Leo, R.Di, et al. (2016). Diverse mechanisms of metaeffector activity in an intracellular bacterial pathogen, Legionella pneumophila. Molecular Systems Biology, 12:893. doi: 10.15252/MSB.20167381
Wagner, N., Avram, O., Gold-Binshtok, D., Zerah, B., Teper, D., and Pupko, T. (2022). Effectidor: an automated machine-learning-based web server for the prediction of type-III secretion system effectors. Bioinformatics 38, 2341–2343. doi: 10.1093/bioinformatics/btac087
Walley, J. W., Shen, Z., McReynolds, M. R., Schmelz, E. A., and Briggs, S. P. (2018). Fungal-induced protein hyperacetylation in maize identified by acetylome profiling. Proc. Natl. Acad. Sci. U. S. A. 115, 210–215. doi: 10.1073/pnas.1717519115
Wang, Y., Zhao, A., Morcillo, R. J. L., Yu, G., Xue, H., Rufian, J. S., et al. (2021). A bacterial effector protein uncovers a plant metabolic pathway involved in tolerance to bacterial wilt disease. Molecular Plant 14, 1281–1296. doi: 10.1016/J.MOLP.2021.04.014
Waterhouse, A., Bertoni, M., Bienert, S., Studer, G., Tauriello, G., Gumienny, R., et al. (2018). SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 46, W296–W303. doi: 10.1093/nar/gky427
Wei, H.-L., Zhang, W., and Collmer, A. (2018). Modular study of the type III effector repertoire in Pseudomonas syringae pv. tomato DC3000 reveals a matrix of effector interplay in pathogenesis. Cell Rep. 23, 1630–1638. doi: 10.1016/j.celrep.2018.04.037
Wei, Y., Balaceanu, A., Rufian, J. S., Segonzac, C., Zhao, A., Morcillo, R. J. L., et al. (2020). An immune receptor complex evolved in soybean to perceive a polymorphic bacterial flagellin. Nature Communications 11, 1–11. doi: 10.1038/s41467-020-17573-y
Willems, P., Ndah, E., Jonckheere, V., Van Breusegem, F., and Van Damme, P. (2022). To New Beginnings: Riboproteogenomics Discovery of N-Terminal Proteoforms in Arabidopsis Thaliana. Front. Plant Sci. 12. doi: 10.3389/FPLS.2021.778804
Winkler, J., Mylle, E., De Meyer, A., Pavie, B., Merchie, J., Grones, P., et al. (2021). Visualizing protein–protein interactions in plants by rapamycin-dependent delocalization. Plant Cell 33, 1101–1117. doi: 10.1093/plcell/koab004
Wu, W., Luo, X., Chen, X., Wang, L., Wang, K., Tian, S., et al. (2022). Ralstonia solanacearum effector RipAF1 ADP-ribosylates host FBN1 to induce resistance against bacterial wilt. BioRxiv 2022:04.487053. doi: 10.1101/2022.04.04.487053
Wu, D., von Roepenack-Lahaye, E., Buntru, M., de Lange, O., Schandry, N., Pérez-Quintero, A. L., et al. (2019). A plant pathogen type III effector protein subverts translational regulation to boost host polyamine levels. Cell Host Microbe 26, 638–649.e5. doi: 10.1016/j.chom.2019.09.014
Wulf, J., Pascuzzi, P. E., Fahmy, A., Martin, G. B., and Nicholson, L. K. (2004). The solution structure of type III effector protein AvrPto reveals conformational and dynamic features important for plant pathogenesis. Structure 12, 1257–1268. doi: 10.1016/j.str.2004.04.017
Xia, Y., Zou, R., Escouboué, M., Zhong, L., Zhu, C., Pouzet, C., et al. (2021). Secondary-structure switch regulates the substrate binding of a YopJ family acetyltransferase. Nat. Commun. 12:5969. doi: 10.1016/j.chom.2019.09.014
Xian, L., Yu, G., Wei, Y., Xue, H., Morcillo, R. J. L., Correspondence, A. P. M., et al. (2020). A Bacterial Effector Protein Hijacks Plant Metabolism to Support Pathogen Nutrition. Cell Host and Microbe 28, 548–557.e7. doi: 10.1016/j.chom.2020.07.003
Xiang, T., Zong, N., Zhang, J., Chen, J., Chen, M., and Zhou, J.-M. (2010). BAK1 is not a target of the Pseudomonas syringae effector AvrPto. Mol. Plant-Microbe Interact. 24, 100–107. doi: 10.1094/MPMI-04-10-0096
Xiang, T., Zong, N., Zou, Y., Wu, Y., Zhang, J., Xing, W., et al. (2008). Pseudomonas syringae effector AvrPto blocks innate immunity by targeting receptor kinases. Curr. Biol. 18, 74–80. doi: 10.1016/j.cub.2007.12.020
Xiou, X., Bihao, C., Guannan, L., Jianjun, L., Qinghua, C., Jin, J., et al. (2015). Functional characterization of a putative bacterial wilt resistance gene (RE-bw) in eggplant. Plant Mol. Biol. Report. 33, 1058–1073. doi: 10.1007/s11105-014-0814-1
Xue, Y., Ren, J., Gao, X., Jin, C., Wen, L., and Yao, X. (2008). GPS 2.0, a tool to predict kinase-specific phosphorylation sites in hierarchy. Mol. Cell. Proteomics 7, 1598–1608. doi: 10.1074/MCP.M700574-MCP200
Yang, X., Guo, Y., Luo, J., Pu, X., and Li, M. (2013). Effective identification of Gram-negative bacterial type III secreted effectors using position-specific residue conservation profiles. PLoS One 8:e84439. doi: 10.1371/journal.pone.0084439
Yu, X., Decker, K. B., Barker, K., Neunuebel, M. R., Saul, J., Graves, M., et al. (2015). Host-pathogen interaction profiling using self-assembling human protein arrays. J. Proteome Res. 14, 1920–1936. doi: 10.1021/pr5013015
Yu, G., Derkacheva, M., Rufian, J. S., Brillada, C., Kowarschik, K., Jiang, S., et al. (2022). The Arabidopsis E3 ubiquitin ligase PUB4 regulates BIK1 and is targeted by a bacterial type-III effector. The EMBO J. 41:e107257. doi: 10.15252/EMBJ.2020107257
Yu, X., Noll, R. R., Romero Dueñas, B. P., Allgood, S. C., Barker, K., Caplan, J. L., et al. (2018). Legionella effector AnkX interacts with host nuclear protein PLEKHN1. BMC Microbiol. 18:5. doi: 10.1186/s12866-017-1147-7
Yu, G., Xian, L., Xue, H., Yu, W., Rufian, J. S., Sang, Y., et al. (2020). A bacterial effector protein prevents MAPK-mediated phosphorylation of SGT1 to suppress plant immunity. PLoS Pathog. 16:e1008933. doi: 10.1371/journal.ppat.1008933
Zhang, Z.-M., Ma, K.-W., Gao, L., Hu, Z., Schwizer, S., Ma, W., et al. (2017). Mechanism of host substrate acetylation by a YopJ family effector. Nat. Plants 3:17115. doi: 10.1038/nplants.2017.115
Zheng, X., Li, X., Wang, B., Cheng, D., Li, Y., Li, W., et al. (2019). A systematic screen of conserved Ralstonia solanacearum effectors reveals the role of RipAB, a nuclear-localized effector that suppresses immune responses in potato. Mol. Plant Pathol. 20, 547–561. doi: 10.1111/mpp.12774
Zhuo, T., Wang, X., Chen, Z., Cui, H., Zeng, Y., Chen, Y., et al. (2020). The Ralstonia solanacearum effector RipI induces a defence reaction by interacting with the bHLH93 transcription factor in Nicotiana benthamiana. Mol. Plant Pathol. 21, 999–1004. doi: 10.1111/MPP.12937
Keywords: type III effector, Ralstonia solanacearum, prediction, functional characterization, immunity, pathogenicity
Citation: De Ryck J, Van Damme P and Goormachtig S (2023) From prediction to function: Current practices and challenges towards the functional characterization of type III effectors. Front. Microbiol. 14:1113442. doi: 10.3389/fmicb.2023.1113442
Edited by:
Johannes Stuttmann, Martin Luther University of Halle-Wittenberg, GermanyReviewed by:
Gongyou Chen, Shanghai Jiao Tong University, ChinaJose Sebastian Rufian, Universidad de Málaga, Spain
Copyright © 2023 De Ryck, Van Damme and Goormachtig. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sofie Goormachtig, sogoo@psb.ugent.be