BACKGROUND
ralstonia solanacearum is a widely distributed soil-borne phytopathogen belonging to the beta subdivision of proteobacteria  <cit> .
it causes lethal bacterial wilt of more than  <dig> plant species, including economically important crops  <cit> .
among the pathogenicity determinants of this bacterium, the type iii secretion system  plays a crucial role because mutants unable to produce this specialized secretion machinery are unable to cause disease on plants  <cit> .
this t3ss ensures the direct translocation of type iii effector  proteins from the bacterium to the plant cell cytosol  <cit> .
these t3e are presumed to perturb host cell processes and modulate plant innate immunity to allow bacterial infection  <cit> .
phylogenetic analyses of ralstonia strains causing wilt diseases revealed an extensive diversity  <cit>  and this group of organisms is now commonly called the r. solanacearum species complex   <cit> .
this species complex includes strains with broad and narrow host ranges with different geographic origins.
based on phylogenetic analyses and on comparative genomic hybridization, the rssc has been classified in four phylogenetic groups called phylotypes, which reflect their origins as follows: asia , the americas , africa  or indonesia   <cit> .
to date,  <dig> strains belonging to the rssc have been completely sequenced.
pioneering studies have established that t3e repertoires are highly variable among strains and shape the host range of bacterial pathogens  <cit> .
first exhaustive inventories of rssc t3e using different in silico or experimental approaches were made in phylotype  <dig> strains gmi <dig>  <cit>  and rs <dig>  <cit> .
gmi <dig> and rs <dig> have almost identical repertoires that comprise  <dig> and  <dig> t3e for which t3ss-dependent plant cell targeting have been experimentally validated in rs <dig>  <cit> .
a feature of these repertoires is the existence of multigenic t3e families  <cit> .
functional studies have been carried out on members of the gala family, which are proteins with f-box and leucine rich repeat domains collectively required for full virulence  <cit> , and members of the popp family, which includes the avirulence proteins popp <dig>  <cit>  and popp <dig>  the latter possessing acetyltransferase activity  <cit> .
recently a functional analysis of the awr family demonstrated that some awr t3e induce cell death necrotic reactions on plants and are required for full virulence  <cit> .
the genome sequence data from strains representative of the biodiversity of the rssc opens the way towards understanding the evolutionary processes that structured their t3e gene repertoire.
this will also provide clues towards defining what makes a given strain more aggressive than others on a specific host.
however such comparative genomic approaches are actually hampered by the fact that t3e inventories in multiple strains have not been accurately established: several t3e genes have been overlooked by automatic annotation programs and/or have been incorrectly predicted.
moreover, the lack of a unified nomenclature for rssc t3e is confusing for a non-expert since many t3e genes from rssc strains have different names in the published literature .
this doesn’t help the already difficult task of identifying orthologous and paralogous genes in strains harboring between  <dig> to  <dig> t3e genes.
this work presents an integrative and comprehensive database for the t3e of the rssc.
this database is a compendium of manually re-annotated genes across  <dig> sequenced strains and ordered with a novel and unifying nomenclature.
this database is publicly available for browsing and retrieving data and information.
our analyses on this particular gene set at the forefront of the interaction between the bacteria and its host, provides new insight into their evolutionary history and their potential contribution to host specificity
RESULTS
ralstonia solanacearum t3e database
inventory and re-annotation of t3e genes in the rssc
our goal is to provide a comprehensive and an as exhaustive as possible inventory of t3e in the rssc as a public database from which curated information can be retrieved.
to this end, we manually curated and compiled the t3e genes from eleven sequenced strains representative of the genetic diversity of the rssc .
the workflow of the retrieval and annotation of the t3e genes from the rssc genomes as well as the main outputs of this analysis are shown in figure  <dig>  the inventory of t3e in the published rssc genome sequences was primarily based on homology searches with the established repertoires of strains gmi <dig>  <cit>  and rs <dig>  <cit> .
identification of additional t3e was conducted using criteria defined previously  <cit>  to mine the gmi <dig> genome:  homology to known t3e in other bacterial species  presence of a hrpii box in the promoter region since 52/ <dig> t3e gene promoters harbor this cis-regulatory element in gmi <dig>  <cit> ,  existence of specific amino acid distribution biases in the 50 n-terminal domain  <cit> .
these two latter criteria were hampered by the fact that many t3e genes have wrongly annotated start codons.
hence all the genes possessing a putative hrpii box were inspected for potential start codon errors before being included in the t3e annotation workflow .
this process led to the “discovery” of twenty new t3e genes , and the re-annotation of 34% of the existing rssc t3e genes.
altogether these changes affect 39% of the rssc t3e dataset  submitted or already present in genbank to date.
identification of t3e candidate genes in rssc strains
a mining of the genome of nine rssc strains from phylotypes  <dig>   <dig> and  <dig> for previously undescribed t3e gene families was performed based on the criteria listed above  <cit> .
in this process, we only kept the t3e candidates strictly fitting with both criteria  and  described above.
this search yielded  <dig> rssc t3e candidates, for which t3ss-dependent translocation is not yet demonstrated.
these  <dig> hypothetical t3e gene families are listed in the additional file  <dig> as well as in the rssc-t3e database.
most of the corresponding genes did not display homology to any other known proteins, except for families rssc-t3e-hyp <dig>  hyp <dig> and hyp <dig> having homologues only in acidovorax spp or xanthomonas spp, which are both plant pathogenic bacteria.
pseudogenes
in many cases, t3e genes appeared to have frameshift mutations or to be split into several independent open reading frames on the assembled genomes.
this could be due to mutations leading to gene inactivation or, more probably, to sequence and assembly errors in the available genome sequences.
it should be noted that there are important differences in terms of quality in the available assembled genomes .
in some other cases, genome sequence gaps resulted in incomplete t3e gene prediction.
many genes encoding t3e with internal repeats are often predicted as truncated or incomplete, probably due to the difficulty to assemble repeat-containing short sequence reads .
frameshift-mutated and incomplete t3e genes were included in the rssc-t3e database and are distinguished by the prefix fs  before the gene name.
future re-sequencing should verify the current pseudogene status of these genes.
probable non-functional pseudogenes are also listed in the rssc-t3e database .
these pseudogenes correspond to genes or gene fragments which are either gene remnants, open reading frames disrupted by a transposable element insertion or frameshift mutated genes confirmed after re-sequencing.
the number of predicted pseudogenes varies from one to eight among the eleven strains analyzed .
however, the formal distinction between a pseudogene and a functional gene is difficult to establish without experimental validation  <cit> .
in some cases, the absence of specific domains  raises the question of the functionality of the corresponding protein.
the rssc-t3e database interface
the dataset corresponding to the lists and expert annotation of validated and candidate t3e in the  <dig> sequenced strains representative of the  <dig> rssc phylotypes were compiled in a web interface named “ralstonia t3e”  designed to provide the user with a convenient and straightforward access to all the underlying data.
the home page provides a synthetic table displaying the distribution of the  <dig> t3e gene families in the rssc strains under the proposed nomenclature .
this table summarises for each strain whether a gene member is present , absent, or is predicted as being not functional .
a specific colour code also indicates genes with putative frameshift mutations.
this information is also available as a table in the additional file  <dig>  the clickable t3e genes provide a link to multifasta files of the curated nucleotide and protein sequences as well as view of the corresponding dna and protein alignments  <cit> .
tab-style navigation provides a link to the  <dig> t3e candidate genes as well as a link to different services like “scanyourgenome” , pat scan, hmscan and blast.
proposed guidelines for the nomenclature of t3es in rssc strains
the recent availability of complete genome sequences for a number of rscc strains has led to a significant increase in the rate of t3e discovery.
however, the absence of a systematic nomenclature has resulted in multiple names being assigned to the same t3e gene.
some genes were named as brg   <cit>  or hpx   <cit>  genes based on regulation studies/screens or as rip  genes  <cit> .
we propose the usage of the generic term of rip for renaming all the t3e genes in the rssc, a term previously used after demonstration of the translocation of these effectors into plant cells  <cit> .
this new nomenclature should follow the rules defined previously for naming the p. syringae t3e  <cit> ; such as: ripxy#strain, wherein the gene is indicated by alphabetic characters, paralogous genes in numerically characters, and the strain in subscript.
the proposed attribution of this novel nomenclature to known translocated rssc t3e is presented in table  <dig> .
r. solanacearum
species complex and proposal for a unified nomenclature
a representative gene member for each family is provided  with other names published in the literature.
homologues t3e from pseudomonas syringae sp.
or xanthomonas sp.
are indicated.
the last column lists t3e for which type  <dig> secretion system-dependent secretion or translocation was experimentally demonstrated.
after identifying groups of homologous genes by reciprocal best hit in the curated list of rssc likely t3e genes, we concentrated our effort in grouping the different genes in orthologous groups and naming then accordingly.
three situation can occur:  a single hit  in each strain, with conservation of synteny on the genome;  a single hit  in each strain, but with a breach of synteny for at least one of the homologous genes;  multiple hits  in different strains.
in the first case a single orthologous group is defined irrespective of the pairwise identity between the orthologous genes.
this can be exemplified by ripb a single gene present in all strains with pairwise amino acid identity ranging from  <dig> to 100%.
another case is ripu also a single gene present in all strains with a strict conservation of synteny, but with surprising divergent members .
even though it is likely that ripu has evolved different functions in the different strains, based on the likely common ancestral origin suggested by the conservation of synteny  <cit> , we advocate for keeping a single orthology group.
in the second situation, an apparent single orthologous group exists but differences in synteny support a scenario of gene duplication followed by gene loss or lateral gene transfer between strains.
here we favour synteny as a ruler for ortholog definition  <cit> .
this is exemplified by ripo <dig> and ripo <dig>  the latter being present only in the strain r <dig>  devoid of ripo <dig>
finally when there are strains with two or more paralogous genes, again we favour the synteny rule to identify groups of orthology  <cit> .
a careful phylogenetic reconstruction for these homologous genes across the whole species complex  illustrates the accuracy of the orthology attributions  <cit> .
these phylogenetic trees also highlighted the existence of two paralogs in several strains that clearly belong to a clade defined as an orthologous group .
we believe that these paralogs result from strain specific  recent gene duplication.
we thus choose to name these genes in a way that indicates their recent evolution: e.g.
ripa5_1molk <dig> and ripa5_2molk2; ripf1_1cmr <dig> and ripf1_2cmr <dig> etc.…the rule of synteny is conserved since we verified that all these genes have indeed a conserved synteny .
suggested name reassignment of previously characterized r. solanacearum t3e
whenever possible the proposed new nomenclature conserves the original letter designations used in previous annotation e.g.
ripp <dig> is popp <dig>  <cit> ; ripp <dig> is popp <dig>  <cit> ; ripaa is avra  <cit> .
in the case of paralogous genes, the names are, for instance: ripg <dig>  ripg <dig>  …to ripg <dig> for the gala gene family  <cit> ; ripa <dig>  ripa <dig>  …to ripa <dig> for the awr family  <cit> .
in a few cases, there is evidence for recent t3e gene duplications resulting in two or more gene copies in a single given strain, e.g.
strain psi <dig> harbors  <dig> copies of ripg <dig>  <cit>  and  <dig> copies of riph2: these were renamed ripg1_ <dig>  ripg1_ <dig>  ripg1_ <dig> and riph2_ <dig>  riph2_ <dig>  respectively, to differentiate them from the other riph and ripg genes in this strain .
in addition, a rip name is proposed for the  <dig> t3e previously identified as pop  <cit>  or avr  <cit> .
the pop designation is historical and was formerly coined when r. solanacearum was known as pseudomonas solanacearum <cit> , the “avr“ term was solely used for the avra avirulence protein identified in  <dig>  <cit> .
these designations can be confusing because the pop term has also been used to name some pseudomonas aeruginosa t3e  <cit>  and avra also refer to an unrelated t3e from salmonella species  <cit> .
“scanyourgenome” a bioinformatic tool for detecting t3e orthologs
in order to swiftly analyse the t3e content of newly produced genome sequences, we developed a protocol for the identification of putative effector candidates.
this pipeline is based on a de novo effectome prediction using t3e models.
then each candidate is tested using different methods with decreasing stringency to assign them to the most probable known effector gene .
this protocol was first tested on reference genomes used above for manual annotation of the t3e genes in order to calibrate the detection parameters  before using it for predicting t3e in the recently published draft genomes of strains k <dig>  <cit> , fqy_ <dig>  <cit>  and y <dig>  <cit> .
this analysis yielded a prediction of  <dig>   <dig> and  <dig> potential t3e encoding genes encoded respectively by the k <dig>  fqy_ <dig> and y <dig> genomes, .
the gene model prediction takes into account possible frameshifts, also when the gene is shorter than 80% of the average length of the other alleles of this rip gene, the predicted gene is tagged as potential pseudogene.
both frameshift and pseudogene annotations appear in the prediction.
this orthology search engine and the consequent rip assignment are available to the community for queries of draft or complete genome sequences.
for shorter gene sequences a more straightforward blast is advised.
the advantage of a sliding scale of orthology detection is the possibility to unequivocally assign each potential t3e gene to a specific orthologous group.
whenever a new candidate t3e gene, experimentally validated as being secreted or translocated into plant cells, will not retrieve an already labelled orthologous rip family, this gene will be assigned the next available rip code.
evolutionary dynamics of rip genes
classification of paralogous rip genes
a specific feature of r. solanacearum t3es is the abundance of paralogous rip genes in all the strains sequenced to date.
some of these paralogous genes are well represented in strains from the four phylotypes, hence they probably originated from ancient duplications in the common ancestor of these diverse strains.
this was well documented for the ripg1-g <dig>  <cit>  and the ripa1-a <dig>  <cit>  paralogous gene families and is probably also true for riph1-h <dig> and rips1-s <dig>  although all strains contain members of these paralogous family, the likely ancient duplications doesn’t exclude some phylotype specificities explained by loss or more simply by recent duplications e.g.
ripa <dig> and rips <dig> seem to be specific to phylotype  <dig>  ripg <dig> is only found in cmr <dig>  the sole representative of phylotype 3; and riph <dig> seems to be specific of the phylotype  <dig> strains .
a second group of paralogous rip genes is characterised by a smaller number  of paralogous sequences in a given strain.
phylogenetic analyses were used to estimate the evolutionary relationships between paralogues using sequence data from the  <dig> rssc representative strains.
we defined eight additional rip genes  .
several of these paralogous genes, such as ripc <dig> or ripo <dig>  seem to differ significantly from ripc <dig> and ripo <dig> respectively and could have originated through lateral gene transfer  since homologous genes exist in other bacterial species.
for the gene families present in most of the rssc strains , the genes are located in each strain in a similar genomic context, an observation which also supports a common evolutionary origin.
but distribution of some paralogs can be variable among strains: .i.e.
ripe <dig> seems to be ubiquitously present whereas ripe <dig> is absent in phylotype  <dig> strains.
protein sequence analyses indicated that ripar, ripaw, ripv <dig>  ripv <dig> and ripbg contain putative ubiquitin-ligase domains , likewise, ripj, ripk, ripae, ripbc, ripp <dig> and ripp <dig> could all potentially display acetyltransferase activity .
notwithstanding this apparent functional conservation, the sequences of these t3e genes have diverged significantly and can’t be assigned in orthologous goups.
it has to be noted that the numerical identification of the two ripp <dig> and ripp <dig>  and the pseudogene ripp3gmi <dig> is used in reference to their previous names popp <dig>  <cit>  , popp <dig>  <cit>   and popp <dig>  <cit> .
this is an exception to the previous rule as we don’t consider these to be paralogs.
horizontally acquired rip genes
the detection of horizontal gene transfer  events in a given bacterial genome can be performed retrospectively through bioinformatics-based comparative analyses  <cit> .
a frequent hallmark of genes with an extrinsic origin is the difference in gc content of these genes compared of the mean content of the host genome  <cit> .
thirteen rip genes exhibit a mean gc% below 60%  .
in several cases, the t3e gene is physically associated with insertion sequence elements , integrases  or are part of prophage sequences integrated in the genome .
from these observations, we can assume that bacteriophage-mediated transfer appears to be an efficient mean for lateral transfer of these t3e in the rssc.
phylogenetic analyses also provided interesting insights into possible hgt with other bacterial plant pathogens.
for example, ripc2cfbp <dig>  outgroup of the ripc <dig> gene family, could derive from the xopc t3e from xanthomonas spp .
furthermore, the low gc content of ripc2
cfbp2957
supports the hypothesis of an hgt, with the possibility of a shared common ancestor between ripc2
cfbp2957
and xopc.
similar observations can be made with ripo2
r.syzygii r24
, ripaf2
r.syzygii r24
, ripe <dig> , ripp <dig> , ripax <dig>  and riph <dig> , see additional file  <dig>  together with riptal , already suspected of inter-species transfer  <cit> , this analysis thus provided a total of seven t3e genes that could have been acquired through hgt.
evidence of phylogenetic incongruences
examination of the intra-family phylogenetic relationships of t3e genes distributed in the nine rssc sequenced strains revealed in some cases incongruences with the species phylogenetic tree.
this can be illustrated by individual rip contradicting the species phylogeny like ripg7cmr <dig> <cit> , ripd cmr <dig>  riph2_1po <dig> and ripax1po <dig>  which could be indicative of rapidly evolving or horizontally acquired genes .
some other conflicting phylogenies can’t be directly associated with a single divergent gene.
this is the case for ripi, ripu and ripac which are present in most of the rssc sequenced strains  but with great sequences divergence .
the only strong evidence for them being orthologs is the fact that ripac and ripu genes are located in two highly syntenic regions with their respective flanking genes strictly following the species phylogeny .
this suggests that ripu and ripac evolved faster in some strains  resulting in this particular high sequence polymorphism  <cit> .
another example of discrepancy between species and gene phylogeny is for ripaa.
here the increased polymorphism is correlated with the presence of a hypervariable domain consisting of variable number of tandem repeats  <cit> .
several rip genes underwent selection and recombination
after excluding from the datasets the likely pseudogenes, all rip genes with more than  <dig> orthologs  were analysed for traces of recurrent diversifying positive selection.
the analysis performed here was carried out like described previously  <cit> , except that gene phylogenies were inferred using one-ratio codon model m <dig>  <cit> .
the full results are displayed in the additional file  <dig>  considering that some of the datasets were rather small we concentrated on identifying rip genes with strong indications of positive selection.
this was the case for the nine following rip genes: ripaa, ripaj, ripat, ripaw, ripbd, ripd, ripg <dig>  riph <dig> and rips <dig> with three out of three likelihood ration tests  for positive selection being significant .
six out of these  <dig> rip genes have an estimated proportion of sites under positive selection higher than 5%, with the highest level reached for ripaj and ripg <dig>  in agreement with a previous analysis  <cit> .
t3e gene
n
e
r
a
b
c
avalues supporting evidence for recombination are shown in bold.
bfor the presented  <dig> genes all three lrts for positive selection were significant, as well as the lrt comparing m <dig> vs m <dig> supporting strong variability of selection pressure among sites.
codon models are as described in  <cit> .
cestimates of selection regimes are according to model m <dig> if lrt comparing m8a and m <dig> was significant.
otherwise, selection regimes are reported according to model m8a.
for strict and relaxed negative selection, the average omega value over respective selection classes is shown.
note that percentages for the four categories do not always add up to 100% due to rounding.
importantly, the presence of a high degree of recombination can hamper lrts for positive diversifying selection, leading to false positives  <cit> .
however inference of recombination can also be affected by selection forces  <cit> .
this is why we systematically analysed all data for evidence of recombination .
table  <dig> also displays the results of tests for recombination for the nine previously identified rip genes.
among these, only two  could also be affected by recombination, while for ripaa the evidence of recombination is not clear-cut.
the interplay between selection and recombination was already disentangled previously for ripg <dig>  <cit> , with the conclusion that there is indeed a strong likelihood of positive selection acting on this gene.
here we won’t address the question further for ripaa and ripaw but a future analysis with more allelic variants should be informative.
it is interesting to note that in the multigene paralogous families there seems to be one member under positive selection: riph <dig>  rips <dig>  ripg <dig>  when we consider only  <dig> out of  <dig> lrts for positive selection , we can define  <dig> more rip coding sequences with evidence for positive selection, out of which  <dig> belong to the above-mentioned paralogous families .
it is tempting to speculate that after duplications some of the paralogous genes could have undergone sub- or neo-functionalisation allowing the cognate rip proteins to adapt to evolving plant targets or evade from host immunity.
comparative genomics and functional implications
the rssc t3e core set: a large group of conserved effectors
the establishment of a near-complete t3e repertoire in strains representative of the large phylogenetic diversity of the rscc allows a more specific and accurate comparison than those based on comparative genomic hybridizations  <cit> .
we performed t3e repertoire comparisons using the following criteria:  rip genes listed as pseudogenes in the database were considered non-functional but those listed as containing frameshifts were considered as functional genes.
the assumption that all the frameshifts are due to sequencing errors is probably an overestimation.
since we can’t validate this experimentally, and considering that the number of frameshifts identified is inversely correlated with the genomic sequence quality, we will keep this assumption.
this is exemplified with gm <dig> and cfbp <dig> high quality genomes, not containing a single frameshift mutation in their t3e genes.
the  <dig> hypothetical t3e newly identified in the different strains were also included in the repertoire for comparisons.
the rssc is divided in three main phylogenetic clades corresponding to phylotypes  <dig>   <dig> and   <cit> .
a first comparison showed that  <dig> rip gene families are present in the  <dig> strains studied.
when the presence requirement is lowered at  <dig> out of the  <dig> strains, the number of gene families reaches  <dig> .
considering that the event of loss of specific t3e genes in some strain lineages is possible , we believe that these  <dig> t3e are a good estimation of the subset of t3e probably present in the ancestral r. solanacearum strain.
interestingly,  <dig> out of  <dig> t3e genes families showing a strong signal of diversifying selection also belong to the core effector group .
it is also interesting to notice that distinct members of paralogous rip families  are also conserved among the  <dig> analyzed strains, indicating that duplications followed by differential evolution of these genes took place early before phylotype divergence  <cit> .
the estimate of  <dig> core t3e certainly reflects the abundance of t3e in the r. solanacearum and, considering its genetic diversity as a species complex, appears significantly higher than the core list identified in p. syringae which is only  <dig> among  <dig> strains  <cit> .
r. solanacearum ancestor presumably possessed more than  <dig> t3e, which were possibly acquired from the bacterial and phage communities in the soil or aquatic reservoirs.
rip
genes.
circles group  genes under strong positive selection,  and  genes belonging to the core group of t3e conserved in  <dig> out of the  <dig> rssc genome sequences.
t3e repertoire comparisons provide no clues on host specificity determinants
phylotypes 1– <dig>   <dig> and  <dig> are the main genetic groups structuring the rssc  <cit> , a comparison of the t3e repertoires  from gmi <dig> , cfbp <dig>  and psi <dig> , representing the three species clades and all isolated from tomato, reveals a diversity of  <dig> t3e genes, almost half of which  are conserved among the three strains whereas one third  appears to be strain-specific .
this confirms that a majority of t3e are widely conserved in this species complex but also shows that the strain repertoires are also diversified, as observed in p. syringae <cit>  or xanthomonas sp.
<cit> .
r. solanacearum strains exhibit great variations in host range  <cit>  and it is tempting to speculate that t3e repertoires shape these host range capabilities.
in order to tentatively identify candidate genes involved in host specificity, we performed t3e repertoire comparisons within specific phylogenetic groups such as phylotype  <dig> or  <dig> using strains with marked host range differences .
these comparisons identified strain-specific genes but did not pinpoint strong host-specificity candidates.
indeed, none of the molk <dig> specific t3e is common with those of the bdbr <dig> strain which is also pathogenic on banana; the same is true for potato-associated t3e genes from the po <dig> and uw551/ipo <dig> strains.
although more genomic sequences of rssc strains are needed to perform robust associations between host range and t3e repertoires, these observations already suggest that host-range maybe controlled by multiple or differential combinations of t3e determinants, or determinants others than t3e, or that differences in t3e protein sequence or gene expression might also be involved  <cit> .
similar observations were reported for comparison of p. syringae pathovars t3e repertoires  <cit> , thus reinforcing the idea that a complex genetic basis underlies host range evolution in plant pathogens.
finally, intra-phylotype comparisons suggest that the proportion of conserved t3e is higher in phylotype  <dig> than in phylotype  <dig> strains .
although phylotype  <dig> strains bdbr <dig> and r <dig> have undergone gene reduction potentially affecting this comparison, we still believe that this difference reflects the highest genetic diversity within phylotype  <dig>  <cit>  and could also be associated with the diverse lifestyle among phylotype  <dig> strains  <cit> .
identification of novel t3e gene harboring putative ubiquitin-ligase domains
molecular functions of most r. solanacearum t3e remain unknown, and more than half of the repertoire corresponds to proteins with no structural motif or domain suggestive of function.
the search for functional motifs identified two t3e proteins, ripar and ripaw, carrying a c-terminal domain structurally related to the shigella flexneri ipah ubiquitin ligase domain  <cit> .
although the overall similarity between ipah and ripar/ripaw is low, these r. solanacearum t3e have a c-terminal domain with a predicted structure consisting of  <dig> alpha-helices as determined for ipah family proteins  <cit> .
most of the highly conserved residues in the ipah family, including a highly conserved cysteine residue essential for activity  <cit> , are conserved in ripar and ripaw see sequence alignment in additional file  <dig>  considering the previously identified t3e ripv, a salmonella ssph <dig> homologue  <cit> , and the ripg family members  <cit> , r. solanacearum potentially harbors a total of  <dig> t3e endowed with potential ubiquitin-ligase activity.
this highlights the probable central mechanism consisting in subversion of the host’s ubiquitination system by t3e during plant pathogenesis  <cit> .
the specific case of the ripf translocon proteins
the rssc t3e list include ripf proteins  as substrates of the t3ss since they were identified as translocated into plant cells using the adenylate cyclase reporter assay  <cit> .
ripf proteins are required for the translocation of other t3e and are t3ss translocator proteins presumably acting at the tip of the hrp pilus and inserting into host cell membranes to permit t3e translocation  <cit> .
contrary to the structural components of the t3ss  which are strongly conserved among all the strains from the rssc analyzed to date, a comparative analysis of ripf revealed major differences among the currently sequenced rssc strains.
strains belonging to phylotypes  <dig>   <dig> and  <dig> possess two ripf whereas strains from phylotype  <dig> have only one  as xanthomonas spp.
in phylotypes  <dig> and  <dig> the second gene, formerly named popf <dig>  <cit> , is phylogenetically close to the first one named popf <dig>  however in phylotype  <dig>  the second gene product belongs to a distinct phylogenetic branch, suggesting an ancient divergence from the other ripf1/popf <dig> lineage.
these observations incited us to rename gmi <dig> popf <dig> as ripf1_ <dig>  whereas ripf <dig> is proposed to designate the gene from the phylotype  <dig> .
this peculiar evolutionary history of the ripf family makes this one of the most stringent discriminating probe among all rip genes for distinguishing the three main phylotype groups of the rssc.
the biological implications of this gene duplication of the ripf translocator in some rssc lineages and the structural divergence between the ripf1/ripf <dig> family members are unknown.
in gmi <dig>  ripf1_ <dig> has a major role in t3e translocation in tomato and tobacco whereas ripf1_ <dig> plays a minor role in this process on these hosts  <cit> .
the specific involvement of ripf <dig> and ripf <dig> in pathogenicity of phylotype  <dig> strains will need to be addressed in future studies.
CONCLUSIONS
t3e are essential to r. solanacearum pathogenesis but progress in understanding of their relative contribution to disease through reverse genetic approaches has been hampered by the evidence of functional redundancies, due to the existence of large t3e repertoires.
in this study, we have undertaken groundwork for a global inventory of r. solanacearum t3e at the species level in order to provide to the community a curated dataset, tools and a rationalized nomenclature that should pave the way for future work on rssc effectomics.
we conducted a large scale approach aimed at the identification, expert annotation and phylogenetic analyses of t3e from the rssc, a species complex showing considerable genomic diversity  <cit>  and responsible for one of the most devastating bacterial disease of plants worldwide  <cit> .
our search yielded a total of  <dig> t3e rip genes and  <dig> additional candidate t3e genes distributed among the  <dig> genomes analyzed in this study.
this total of more than  <dig> predicted t3es is significantly higher than the t3e inventories from other bacterial plant pathogens.
indeed, in p. syringae, genome analysis of  <dig> phylogenetically diverse isolates revealed the existence of  <dig> t3e genes  <cit>   whereas this number is estimated to  <dig> in xanthomonas spp  <cit> .
these comparisons highlight the great diversity of t3e genes present in the rssc and the apparent complexity of t3ss-dependent pathogenesis in this species complex.
the rssc t3e also appears to be highly dynamic, as evidenced by the number of t3e under positive selection indicative of possible neo-functionalization or the number of t3e pseudogenes identified in this study.
in particular, there is an obvious tendency to t3e gene decay in r. syzygii which is correlated with the genome reduction in this strain  <cit> .
r. syzygii is an exception among the rssc since it is strictly limited to clover tree, the t3e repertoire reduction in this strain may be a consequence of this host specialization.
on the other hand, the cornucopia of t3e identified in r. solanacearum and other related pathogenic beta-proteobacteria is probably a factor explaining the exceptional adaptation of these pathogens to such a wide diversity of hosts.
importantly, phylogenetic analyses allowed the definition of novel t3e genes, resulting in the definition of new rip genes orthologous group or paralogs of already identified rip genes.
it is conceivable that these newly defined groups correspond to t3e genes with novel functional specificities.
our analysis should also be helpful for refined functional studies:  the ripf1-ripf <dig> translocon proteins appear as major discriminating determinants among the main lineages of the rssc and this probably reflects a fundamental evolutionary divergence  global comparisons of repertoires among genetically diverse strains identified a set of 20– <dig> core t3e widely distributed in the species which could presumably be considered as ancestral t3e important in the interaction of the pathogen with its hosts, and  the identification of t3e displaying a positive selection pattern may provide hints on the determinants evolving under plant selection pressure,  our bioinformatics pipeline is dedicated to rapidly predict and assign rip identifiers to all homologous t3e genes in newly sequenced strains of the rssc.

 ralstonia solanacearum is a widely distributed soil-borne phytopathogen belonging to the beta subdivision of proteobacteria. it causes lethal bacterial wilt of more than plant species, including economically important crops. among the pathogenicity determinants of this bacterium, the type iii secretion system plays a crucial role because mutants unable to produce this specialized secretion machinery are unable to cause disease on plants. this t3ss ensures the direct translocation of type iii effector proteins from the bacterium to the plant cell cytosol. these t3e are presumed to perturb host cell processes and modulate plant innate immunity to allow bacterial infection. phylogenetic analyses of ralstonia strains causing wilt diseases revealed an extensive diversity and this group of organisms is now commonly called the r. solanacearum species complex. this species complex includes strains with broad and narrow host ranges with different geographic origins. based on phylogenetic analyses and on comparative genomic hybridization, the rssc has been classified in four phylogenetic groups called phylotypes, which reflect their origins as follows: asia, the americas, africa or indonesia. to date, strains belonging to the rssc have been completely sequenced. pioneering studies have established that t3e repertoires are highly variable among strains and shape the host range of bacterial pathogens. first exhaustive inventories of rssc t3e using different in silico or experimental approaches were made in phylotype strains gmi and rs. gmi and rs have almost identical repertoires that comprise and t3e for which t3ss-dependent plant cell targeting have been experimentally validated in rs. a feature of these repertoires is the existence of multigenic t3e families. functional studies have been carried out on members of the gala family, which are proteins with f-box and leucine rich repeat domains collectively required for full virulence, and members of the popp family, which includes the avirulence proteins popp and popp the latter possessing acetyltransferase activity. recently a functional analysis of the awr family demonstrated that some awr t3e induce cell death necrotic reactions on plants and are required for full virulence. the genome sequence data from strains representative of the biodiversity of the rssc opens the way towards understanding the evolutionary processes that structured their t3e gene repertoire. this will also provide clues towards defining what makes a given strain more aggressive than others on a specific host. however such comparative genomic approaches are actually hampered by the fact that t3e inventories in multiple strains have not been accurately established: several t3e genes have been overlooked by automatic annotation programs and/or have been incorrectly predicted. moreover, the lack of a unified nomenclature for rssc t3e is confusing for a non-expert since many t3e genes from rssc strains have different names in the published literature. this doesn’t help the already difficult task of identifying orthologous and paralogous genes in strains harboring between to t3e genes. this work presents an integrative and comprehensive database for the t3e of the rssc. this database is a compendium of manually re-annotated genes across sequenced strains and ordered with a novel and unifying nomenclature. this database is publicly available for browsing and retrieving data and information. our analyses on this particular gene set at the forefront of the interaction between the bacteria and its host, provides new insight into their evolutionary history and their potential contribution to host specificity  ralstonia solanacearum t3e database inventory and re-annotation of t3e genes in the rssc our goal is to provide a comprehensive and an as exhaustive as possible inventory of t3e in the rssc as a public database from which curated information can be retrieved. to this end, we manually curated and compiled the t3e genes from eleven sequenced strains representative of the genetic diversity of the rssc. the workflow of the retrieval and annotation of the t3e genes from the rssc genomes as well as the main outputs of this analysis are shown in figure  the inventory of t3e in the published rssc genome sequences was primarily based on homology searches with the established repertoires of strains gmi and rs. identification of additional t3e was conducted using criteria defined previously to mine the gmi genome: homology to known t3e in other bacterial species presence of a hrpii box in the promoter region since 52/ t3e gene promoters harbor this cis-regulatory element in gmi, existence of specific amino acid distribution biases in the 50 n-terminal domain. these two latter criteria were hampered by the fact that many t3e genes have wrongly annotated start codons. hence all the genes possessing a putative hrpii box were inspected for potential start codon errors before being included in the t3e annotation workflow. this process led to the “discovery” of twenty new t3e genes, and the re-annotation of 34% of the existing rssc t3e genes. altogether these changes affect 39% of the rssc t3e dataset submitted or already present in genbank to date. identification of t3e candidate genes in rssc strains a mining of the genome of nine rssc strains from phylotypes and for previously undescribed t3e gene families was performed based on the criteria listed above. in this process, we only kept the t3e candidates strictly fitting with both criteria and described above. this search yielded rssc t3e candidates, for which t3ss-dependent translocation is not yet demonstrated. these hypothetical t3e gene families are listed in the additional file as well as in the rssc-t3e database. most of the corresponding genes did not display homology to any other known proteins, except for families rssc-t3e-hyp hyp and hyp having homologues only in acidovorax spp or xanthomonas spp, which are both plant pathogenic bacteria. pseudogenes in many cases, t3e genes appeared to have frameshift mutations or to be split into several independent open reading frames on the assembled genomes. this could be due to mutations leading to gene inactivation or, more probably, to sequence and assembly errors in the available genome sequences. it should be noted that there are important differences in terms of quality in the available assembled genomes. in some other cases, genome sequence gaps resulted in incomplete t3e gene prediction. many genes encoding t3e with internal repeats are often predicted as truncated or incomplete, probably due to the difficulty to assemble repeat-containing short sequence reads. frameshift-mutated and incomplete t3e genes were included in the rssc-t3e database and are distinguished by the prefix fs before the gene name. future re-sequencing should verify the current pseudogene status of these genes. probable non-functional pseudogenes are also listed in the rssc-t3e database. these pseudogenes correspond to genes or gene fragments which are either gene remnants, open reading frames disrupted by a transposable element insertion or frameshift mutated genes confirmed after re-sequencing. the number of predicted pseudogenes varies from one to eight among the eleven strains analyzed. however, the formal distinction between a pseudogene and a functional gene is difficult to establish without experimental validation. in some cases, the absence of specific domains raises the question of the functionality of the corresponding protein. the rssc-t3e database interface the dataset corresponding to the lists and expert annotation of validated and candidate t3e in the sequenced strains representative of the rssc phylotypes were compiled in a web interface named “ralstonia t3e” designed to provide the user with a convenient and straightforward access to all the underlying data. the home page provides a synthetic table displaying the distribution of the t3e gene families in the rssc strains under the proposed nomenclature. this table summarises for each strain whether a gene member is present, absent, or is predicted as being not functional. a specific colour code also indicates genes with putative frameshift mutations. this information is also available as a table in the additional file the clickable t3e genes provide a link to multifasta files of the curated nucleotide and protein sequences as well as view of the corresponding dna and protein alignments. tab-style navigation provides a link to the t3e candidate genes as well as a link to different services like “scanyourgenome”, pat scan, hmscan and blast. proposed guidelines for the nomenclature of t3es in rssc strains the recent availability of complete genome sequences for a number of rscc strains has led to a significant increase in the rate of t3e discovery. however, the absence of a systematic nomenclature has resulted in multiple names being assigned to the same t3e gene. some genes were named as brg or hpx genes based on regulation studies/screens or as rip genes. we propose the usage of the generic term of rip for renaming all the t3e genes in the rssc, a term previously used after demonstration of the translocation of these effectors into plant cells. this new nomenclature should follow the rules defined previously for naming the p. syringae t3e; such as: ripxy#strain, wherein the gene is indicated by alphabetic characters, paralogous genes in numerically characters, and the strain in subscript. the proposed attribution of this novel nomenclature to known translocated rssc t3e is presented in table . r. solanacearum species complex and proposal for a unified nomenclature a representative gene member for each family is provided with other names published in the literature. homologues t3e from pseudomonas syringae sp. or xanthomonas sp. are indicated. the last column lists t3e for which type secretion system-dependent secretion or translocation was experimentally demonstrated. after identifying groups of homologous genes by reciprocal best hit in the curated list of rssc likely t3e genes, we concentrated our effort in grouping the different genes in orthologous groups and naming then accordingly. three situation can occur: a single hit in each strain, with conservation of synteny on the genome; a single hit in each strain, but with a breach of synteny for at least one of the homologous genes; multiple hits in different strains. in the first case a single orthologous group is defined irrespective of the pairwise identity between the orthologous genes. this can be exemplified by ripb a single gene present in all strains with pairwise amino acid identity ranging from to 100%. another case is ripu also a single gene present in all strains with a strict conservation of synteny, but with surprising divergent members. even though it is likely that ripu has evolved different functions in the different strains, based on the likely common ancestral origin suggested by the conservation of synteny, we advocate for keeping a single orthology group. in the second situation, an apparent single orthologous group exists but differences in synteny support a scenario of gene duplication followed by gene loss or lateral gene transfer between strains. here we favour synteny as a ruler for ortholog definition. this is exemplified by ripo and ripo the latter being present only in the strain r devoid of ripo finally when there are strains with two or more paralogous genes, again we favour the synteny rule to identify groups of orthology. a careful phylogenetic reconstruction for these homologous genes across the whole species complex illustrates the accuracy of the orthology attributions. these phylogenetic trees also highlighted the existence of two paralogs in several strains that clearly belong to a clade defined as an orthologous group. we believe that these paralogs result from strain specific recent gene duplication. we thus choose to name these genes in a way that indicates their recent evolution: e.g. ripa5_1molk and ripa5_2molk2; ripf1_1cmr and ripf1_2cmr etc.…the rule of synteny is conserved since we verified that all these genes have indeed a conserved synteny. suggested name reassignment of previously characterized r. solanacearum t3e whenever possible the proposed new nomenclature conserves the original letter designations used in previous annotation e.g. ripp is popp; ripp is popp; ripaa is avra. in the case of paralogous genes, the names are, for instance: ripg ripg …to ripg for the gala gene family; ripa ripa …to ripa for the awr family. in a few cases, there is evidence for recent t3e gene duplications resulting in two or more gene copies in a single given strain, e.g. strain psi harbors copies of ripg and copies of riph2: these were renamed ripg1_ ripg1_ ripg1_ and riph2_ riph2_ respectively, to differentiate them from the other riph and ripg genes in this strain. in addition, a rip name is proposed for the t3e previously identified as pop or avr. the pop designation is historical and was formerly coined when r. solanacearum was known as pseudomonas solanacearum, the “avr“ term was solely used for the avra avirulence protein identified in. these designations can be confusing because the pop term has also been used to name some pseudomonas aeruginosa t3e and avra also refer to an unrelated t3e from salmonella species. “scanyourgenome” a bioinformatic tool for detecting t3e orthologs in order to swiftly analyse the t3e content of newly produced genome sequences, we developed a protocol for the identification of putative effector candidates. this pipeline is based on a de novo effectome prediction using t3e models. then each candidate is tested using different methods with decreasing stringency to assign them to the most probable known effector gene. this protocol was first tested on reference genomes used above for manual annotation of the t3e genes in order to calibrate the detection parameters before using it for predicting t3e in the recently published draft genomes of strains k, fqy_ and y. this analysis yielded a prediction of and potential t3e encoding genes encoded respectively by the k fqy_ and y genomes,. the gene model prediction takes into account possible frameshifts, also when the gene is shorter than 80% of the average length of the other alleles of this rip gene, the predicted gene is tagged as potential pseudogene. both frameshift and pseudogene annotations appear in the prediction. this orthology search engine and the consequent rip assignment are available to the community for queries of draft or complete genome sequences. for shorter gene sequences a more straightforward blast is advised. the advantage of a sliding scale of orthology detection is the possibility to unequivocally assign each potential t3e gene to a specific orthologous group. whenever a new candidate t3e gene, experimentally validated as being secreted or translocated into plant cells, will not retrieve an already labelled orthologous rip family, this gene will be assigned the next available rip code. evolutionary dynamics of rip genes classification of paralogous rip genes a specific feature of r. solanacearum t3es is the abundance of paralogous rip genes in all the strains sequenced to date. some of these paralogous genes are well represented in strains from the four phylotypes, hence they probably originated from ancient duplications in the common ancestor of these diverse strains. this was well documented for the ripg1-g and the ripa1-a paralogous gene families and is probably also true for riph1-h and rips1-s although all strains contain members of these paralogous family, the likely ancient duplications doesn’t exclude some phylotype specificities explained by loss or more simply by recent duplications e.g. ripa and rips seem to be specific to phylotype ripg is only found in cmr the sole representative of phylotype 3; and riph seems to be specific of the phylotype strains. a second group of paralogous rip genes is characterised by a smaller number of paralogous sequences in a given strain. phylogenetic analyses were used to estimate the evolutionary relationships between paralogues using sequence data from the rssc representative strains. we defined eight additional rip genes. several of these paralogous genes, such as ripc or ripo seem to differ significantly from ripc and ripo respectively and could have originated through lateral gene transfer since homologous genes exist in other bacterial species. for the gene families present in most of the rssc strains, the genes are located in each strain in a similar genomic context, an observation which also supports a common evolutionary origin. but distribution of some paralogs can be variable among strains:.i.e. ripe seems to be ubiquitously present whereas ripe is absent in phylotype strains. protein sequence analyses indicated that ripar, ripaw, ripv ripv and ripbg contain putative ubiquitin-ligase domains, likewise, ripj, ripk, ripae, ripbc, ripp and ripp could all potentially display acetyltransferase activity. notwithstanding this apparent functional conservation, the sequences of these t3e genes have diverged significantly and can’t be assigned in orthologous goups. it has to be noted that the numerical identification of the two ripp and ripp and the pseudogene ripp3gmi is used in reference to their previous names popp, popp and popp. this is an exception to the previous rule as we don’t consider these to be paralogs. horizontally acquired rip genes the detection of horizontal gene transfer events in a given bacterial genome can be performed retrospectively through bioinformatics-based comparative analyses. a frequent hallmark of genes with an extrinsic origin is the difference in gc content of these genes compared of the mean content of the host genome. thirteen rip genes exhibit a mean gc% below 60%. in several cases, the t3e gene is physically associated with insertion sequence elements, integrases or are part of prophage sequences integrated in the genome. from these observations, we can assume that bacteriophage-mediated transfer appears to be an efficient mean for lateral transfer of these t3e in the rssc. phylogenetic analyses also provided interesting insights into possible hgt with other bacterial plant pathogens. for example, ripc2cfbp outgroup of the ripc gene family, could derive from the xopc t3e from xanthomonas spp. furthermore, the low gc content of ripc2 cfbp2957 supports the hypothesis of an hgt, with the possibility of a shared common ancestor between ripc2 cfbp2957 and xopc. similar observations can be made with ripo2 r.syzygii r24, ripaf2 r.syzygii r24, ripe, ripp, ripax and riph, see additional file together with riptal, already suspected of inter-species transfer, this analysis thus provided a total of seven t3e genes that could have been acquired through hgt. evidence of phylogenetic incongruences examination of the intra-family phylogenetic relationships of t3e genes distributed in the nine rssc sequenced strains revealed in some cases incongruences with the species phylogenetic tree. this can be illustrated by individual rip contradicting the species phylogeny like ripg7cmr, ripd cmr riph2_1po and ripax1po which could be indicative of rapidly evolving or horizontally acquired genes. some other conflicting phylogenies can’t be directly associated with a single divergent gene. this is the case for ripi, ripu and ripac which are present in most of the rssc sequenced strains but with great sequences divergence. the only strong evidence for them being orthologs is the fact that ripac and ripu genes are located in two highly syntenic regions with their respective flanking genes strictly following the species phylogeny. this suggests that ripu and ripac evolved faster in some strains resulting in this particular high sequence polymorphism. another example of discrepancy between species and gene phylogeny is for ripaa. here the increased polymorphism is correlated with the presence of a hypervariable domain consisting of variable number of tandem repeats. several rip genes underwent selection and recombination after excluding from the datasets the likely pseudogenes, all rip genes with more than orthologs were analysed for traces of recurrent diversifying positive selection. the analysis performed here was carried out like described previously, except that gene phylogenies were inferred using one-ratio codon model m. the full results are displayed in the additional file considering that some of the datasets were rather small we concentrated on identifying rip genes with strong indications of positive selection. this was the case for the nine following rip genes: ripaa, ripaj, ripat, ripaw, ripbd, ripd, ripg riph and rips with three out of three likelihood ration tests for positive selection being significant. six out of these rip genes have an estimated proportion of sites under positive selection higher than 5%, with the highest level reached for ripaj and ripg in agreement with a previous analysis. t3e gene n e r a b c avalues supporting evidence for recombination are shown in bold. bfor the presented genes all three lrts for positive selection were significant, as well as the lrt comparing m vs m supporting strong variability of selection pressure among sites. codon models are as described in. cestimates of selection regimes are according to model m if lrt comparing m8a and m was significant. otherwise, selection regimes are reported according to model m8a. for strict and relaxed negative selection, the average omega value over respective selection classes is shown. note that percentages for the four categories do not always add up to 100% due to rounding. importantly, the presence of a high degree of recombination can hamper lrts for positive diversifying selection, leading to false positives. however inference of recombination can also be affected by selection forces. this is why we systematically analysed all data for evidence of recombination. table  also displays the results of tests for recombination for the nine previously identified rip genes. among these, only two could also be affected by recombination, while for ripaa the evidence of recombination is not clear-cut. the interplay between selection and recombination was already disentangled previously for ripg, with the conclusion that there is indeed a strong likelihood of positive selection acting on this gene. here we won’t address the question further for ripaa and ripaw but a future analysis with more allelic variants should be informative. it is interesting to note that in the multigene paralogous families there seems to be one member under positive selection: riph rips ripg when we consider only out of lrts for positive selection, we can define more rip coding sequences with evidence for positive selection, out of which belong to the above-mentioned paralogous families. it is tempting to speculate that after duplications some of the paralogous genes could have undergone sub- or neo-functionalisation allowing the cognate rip proteins to adapt to evolving plant targets or evade from host immunity. comparative genomics and functional implications the rssc t3e core set: a large group of conserved effectors the establishment of a near-complete t3e repertoire in strains representative of the large phylogenetic diversity of the rscc allows a more specific and accurate comparison than those based on comparative genomic hybridizations. we performed t3e repertoire comparisons using the following criteria: rip genes listed as pseudogenes in the database were considered non-functional but those listed as containing frameshifts were considered as functional genes. the assumption that all the frameshifts are due to sequencing errors is probably an overestimation. since we can’t validate this experimentally, and considering that the number of frameshifts identified is inversely correlated with the genomic sequence quality, we will keep this assumption. this is exemplified with gm and cfbp high quality genomes, not containing a single frameshift mutation in their t3e genes. the hypothetical t3e newly identified in the different strains were also included in the repertoire for comparisons. the rssc is divided in three main phylogenetic clades corresponding to phylotypes and. a first comparison showed that rip gene families are present in the strains studied. when the presence requirement is lowered at out of the strains, the number of gene families reaches. considering that the event of loss of specific t3e genes in some strain lineages is possible, we believe that these t3e are a good estimation of the subset of t3e probably present in the ancestral r. solanacearum strain. interestingly, out of t3e genes families showing a strong signal of diversifying selection also belong to the core effector group. it is also interesting to notice that distinct members of paralogous rip families are also conserved among the analyzed strains, indicating that duplications followed by differential evolution of these genes took place early before phylotype divergence. the estimate of core t3e certainly reflects the abundance of t3e in the r. solanacearum and, considering its genetic diversity as a species complex, appears significantly higher than the core list identified in p. syringae which is only among strains. r. solanacearum ancestor presumably possessed more than t3e, which were possibly acquired from the bacterial and phage communities in the soil or aquatic reservoirs. rip genes. circles group genes under strong positive selection, and genes belonging to the core group of t3e conserved in out of the rssc genome sequences. t3e repertoire comparisons provide no clues on host specificity determinants phylotypes 1– and are the main genetic groups structuring the rssc, a comparison of the t3e repertoires from gmi, cfbp and psi, representing the three species clades and all isolated from tomato, reveals a diversity of t3e genes, almost half of which are conserved among the three strains whereas one third appears to be strain-specific. this confirms that a majority of t3e are widely conserved in this species complex but also shows that the strain repertoires are also diversified, as observed in p. syringae or xanthomonas sp.. r. solanacearum strains exhibit great variations in host range and it is tempting to speculate that t3e repertoires shape these host range capabilities. in order to tentatively identify candidate genes involved in host specificity, we performed t3e repertoire comparisons within specific phylogenetic groups such as phylotype or using strains with marked host range differences. these comparisons identified strain-specific genes but did not pinpoint strong host-specificity candidates. indeed, none of the molk specific t3e is common with those of the bdbr strain which is also pathogenic on banana; the same is true for potato-associated t3e genes from the po and uw551/ipo strains. although more genomic sequences of rssc strains are needed to perform robust associations between host range and t3e repertoires, these observations already suggest that host-range maybe controlled by multiple or differential combinations of t3e determinants, or determinants others than t3e, or that differences in t3e protein sequence or gene expression might also be involved. similar observations were reported for comparison of p. syringae pathovars t3e repertoires, thus reinforcing the idea that a complex genetic basis underlies host range evolution in plant pathogens. finally, intra-phylotype comparisons suggest that the proportion of conserved t3e is higher in phylotype than in phylotype strains. although phylotype strains bdbr and r have undergone gene reduction potentially affecting this comparison, we still believe that this difference reflects the highest genetic diversity within phylotype and could also be associated with the diverse lifestyle among phylotype strains. identification of novel t3e gene harboring putative ubiquitin-ligase domains molecular functions of most r. solanacearum t3e remain unknown, and more than half of the repertoire corresponds to proteins with no structural motif or domain suggestive of function. the search for functional motifs identified two t3e proteins, ripar and ripaw, carrying a c-terminal domain structurally related to the shigella flexneri ipah ubiquitin ligase domain. although the overall similarity between ipah and ripar/ripaw is low, these r. solanacearum t3e have a c-terminal domain with a predicted structure consisting of alpha-helices as determined for ipah family proteins. most of the highly conserved residues in the ipah family, including a highly conserved cysteine residue essential for activity, are conserved in ripar and ripaw see sequence alignment in additional file considering the previously identified t3e ripv, a salmonella ssph homologue, and the ripg family members, r. solanacearum potentially harbors a total of t3e endowed with potential ubiquitin-ligase activity. this highlights the probable central mechanism consisting in subversion of the host’s ubiquitination system by t3e during plant pathogenesis. the specific case of the ripf translocon proteins the rssc t3e list include ripf proteins as substrates of the t3ss since they were identified as translocated into plant cells using the adenylate cyclase reporter assay. ripf proteins are required for the translocation of other t3e and are t3ss translocator proteins presumably acting at the tip of the hrp pilus and inserting into host cell membranes to permit t3e translocation. contrary to the structural components of the t3ss which are strongly conserved among all the strains from the rssc analyzed to date, a comparative analysis of ripf revealed major differences among the currently sequenced rssc strains. strains belonging to phylotypes and possess two ripf whereas strains from phylotype have only one as xanthomonas spp. in phylotypes and the second gene, formerly named popf, is phylogenetically close to the first one named popf however in phylotype the second gene product belongs to a distinct phylogenetic branch, suggesting an ancient divergence from the other ripf1/popf lineage. these observations incited us to rename gmi popf as ripf1_ whereas ripf is proposed to designate the gene from the phylotype. this peculiar evolutionary history of the ripf family makes this one of the most stringent discriminating probe among all rip genes for distinguishing the three main phylotype groups of the rssc. the biological implications of this gene duplication of the ripf translocator in some rssc lineages and the structural divergence between the ripf1/ripf family members are unknown. in gmi ripf1_ has a major role in t3e translocation in tomato and tobacco whereas ripf1_ plays a minor role in this process on these hosts. the specific involvement of ripf and ripf in pathogenicity of phylotype strains will need to be addressed in future studies.  t3e are essential to r. solanacearum pathogenesis but progress in understanding of their relative contribution to disease through reverse genetic approaches has been hampered by the evidence of functional redundancies, due to the existence of large t3e repertoires. in this study, we have undertaken groundwork for a global inventory of r. solanacearum t3e at the species level in order to provide to the community a curated dataset, tools and a rationalized nomenclature that should pave the way for future work on rssc effectomics. we conducted a large scale approach aimed at the identification, expert annotation and phylogenetic analyses of t3e from the rssc, a species complex showing considerable genomic diversity and responsible for one of the most devastating bacterial disease of plants worldwide. our search yielded a total of t3e rip genes and additional candidate t3e genes distributed among the genomes analyzed in this study. this total of more than predicted t3es is significantly higher than the t3e inventories from other bacterial plant pathogens. indeed, in p. syringae, genome analysis of phylogenetically diverse isolates revealed the existence of t3e genes whereas this number is estimated to in xanthomonas spp. these comparisons highlight the great diversity of t3e genes present in the rssc and the apparent complexity of t3ss-dependent pathogenesis in this species complex. the rssc t3e also appears to be highly dynamic, as evidenced by the number of t3e under positive selection indicative of possible neo-functionalization or the number of t3e pseudogenes identified in this study. in particular, there is an obvious tendency to t3e gene decay in r. syzygii which is correlated with the genome reduction in this strain. r. syzygii is an exception among the rssc since it is strictly limited to clover tree, the t3e repertoire reduction in this strain may be a consequence of this host specialization. on the other hand, the cornucopia of t3e identified in r. solanacearum and other related pathogenic beta-proteobacteria is probably a factor explaining the exceptional adaptation of these pathogens to such a wide diversity of hosts. importantly, phylogenetic analyses allowed the definition of novel t3e genes, resulting in the definition of new rip genes orthologous group or paralogs of already identified rip genes. it is conceivable that these newly defined groups correspond to t3e genes with novel functional specificities. our analysis should also be helpful for refined functional studies: the ripf1-ripf translocon proteins appear as major discriminating determinants among the main lineages of the rssc and this probably reflects a fundamental evolutionary divergence global comparisons of repertoires among genetically diverse strains identified a set of 20– core t3e widely distributed in the species which could presumably be considered as ancestral t3e important in the interaction of the pathogen with its hosts, and the identification of t3e displaying a positive selection pattern may provide hints on the determinants evolving under plant selection pressure, our bioinformatics pipeline is dedicated to rapidly predict and assign rip identifiers to all homologous t3e genes in newly sequenced strains of the rssc. 
pathogenicity of r. solanacearum relies on the type iii secretion system that injects type iii effector proteins into plant cells. we provide the catalogue of t3e in the r. solanacearum species complex, as well as candidates in newly sequenced strains. establishment of t3e repertoires in strains representatives of the r. solanacearum biodiversity allowed determining a set of t3e present in all the strains but provided no clues on host specificity determinants. systematical analyses were performed to detect lateral t3e gene transfer events and identify t3e genes under positive selection. ralstonia solanacearum is a soil-borne beta-proteobacterium that causes bacterial wilt disease in many food crops and is a major problem for agriculture in intertropical regions. type iii effectorralstonia solanacearumselectionhorizontal gene transferhost specificity our analyses also pinpoint the ripf translocon proteins as major discriminating determinants among the phylogenetic lineages. t3e collectively perturb host cell processes and modulate plant immunity to enable bacterial infection. this curated t3e catalog is available on a public website and a bioinformatic pipeline has been designed to rapidly predict t3e genes in newly sequenced strains. r. solanacearum is a heterogeneous species, both phenotypically and genetically, and is considered as a species complex.
 ralstonia solanacearum is a soil-borne beta-proteobacterium that causes bacterial wilt disease in many food crops and is a major problem for agriculture in intertropical regions. r. solanacearum is a heterogeneous species, both phenotypically and genetically, and is considered as a species complex. pathogenicity of r. solanacearum relies on the type iii secretion system that injects type iii effector proteins into plant cells. t3e collectively perturb host cell processes and modulate plant immunity to enable bacterial infection.  we provide the catalogue of t3e in the r. solanacearum species complex, as well as candidates in newly sequenced strains. t3e orthologous groups were defined on phylogenetic bases and ordered using a uniform nomenclature. this curated t3e catalog is available on a public website and a bioinformatic pipeline has been designed to rapidly predict t3e genes in newly sequenced strains. systematical analyses were performed to detect lateral t3e gene transfer events and identify t3e genes under positive selection. our analyses also pinpoint the ripf translocon proteins as major discriminating determinants among the phylogenetic lineages.  establishment of t3e repertoires in strains representatives of the r. solanacearum biodiversity allowed determining a set of t3e present in all the strains but provided no clues on host specificity determinants. the definition of a standardized nomenclature and the optimization of predictive tools will pave the way to understanding how variation of these repertoires is correlated to the diversification of this species complex and how they contribute to the different strain pathotypes. type iii effectorralstonia solanacearumselectionhorizontal gene transferhost specificity 