BACKGROUND
plant cells constantly react to multiple signals that come from the local environment, neighboring cells, or even from other organisms.
depending on the stimuli, plant cells may expand, divide, differentiate, synthesize compounds, prepare against pathogen infection, or induce necrosis
<cit> .
to perceive and receive these signals, plant cells possess complex systems of transmembrane receptor proteins that facilitate communication between the intracellular environment and the outside world.
one of the largest groups of these receptors is the receptor-like kinase  superfamily, which contains over  <dig> members in arabidopsis
<cit> .
rlks are structurally organized into an extracellular domain that can be highly divergent, followed by a transmembrane segment and a conserved intracellular serine/threonine kinase domain.
most rlks are localized in the plasma membrane, although there are also rlk members that are found in the cytoplasm.
in this case, rlks do not possess either an extracellular region or a transmembrane domain and are called receptor-like cytoplasmic kinases .
analyses of arabidopsis rlks by structural comparison of their extracellular region and phylogenetic analysis of their kinase domain revealed that they can be divided into over  <dig> subfamilies
<cit> .
several distinct rlks have been studied in the past decade, and a common theme that has emerged is that binding of a specific signal molecule to their extracellular domain is required to initiate a signal transduction cascade
<cit> .
generally, ligand-receptor interactions at the extracellular domain of rlks initiate the propagation of the signal through the membrane by inducing a conformational change in the receptor kinase domain, which allows interactions with other rlks resulting in homo- or heterodimers.
dimerized rlks are then transphosphorylated by their cytoplasmic kinase domain, leading to both activation of the kinase and establishment of docking sites for phosphorylation of downstream phosphorylation targets
<cit> .
this activation mechanism of plant rlks is similar to that of signal transduction mediated by receptor tyrosine kinases in animal cells, which share a common origin with plant serine/threonine kinases
<cit> .
functional analysis of rlks indicates that the majority of them are associated with plant development or defense response, but there are also rlks involved in cell wall attachment , plant-bacterial symbiotic interactions  and self-incompatibility .
among all rlks, those bearing a leucine-rich repeat  domains are overrepresented in the rlk superfamily, comprising over 38% of arabidopsis rlks, which are distributed into  <dig> subgroups .
the lrr domains in these receptors vary in number  and in the distribution pattern of the lrrs along the extracellular region.
examples of well-known lrr-rlks include clavata <dig>  which controls the size of stem cells in the apical meristem by forming a heterodimer with clavata <dig> and then interacting with clavata <dig> through the extracellular domain
<cit> , and brassinosteroid insensitive- <dig>
<cit> , which perceives brassinosteroids and interacts with it receptor partner, bak <dig>
<cit> .
other functions associated with lrr-rlks include morphogenesis
<cit> , embryogenesis
<cit> , pollen self-incompatibility
<cit>  and responses to environmental signals
<cit> .
in addition, some lrr-rlks are known to function as regulators of defense response to bacterial pathogen
<cit> , necrotrophic fungus
<cit>  and viral infection
<cit> .
most of the characterized rlks are from model plants such as arabidopsis and medicago truncatula, but significant efforts have been made to expand these studies to relevant field crops.
large-scale comparative analyses of arabidopsis rlks with rice
<cit>  and soybean
<cit>  rlks identified over  <dig> kinase proteins in rice and  <dig> in soybean belonging to the rlk superfamily; almost all members were grouped into previously determined arabidopsis rlk subfamilies.
the rlk subfamilies with developmental function have conserved size, whereas those involved in defense response have expanded their members, mainly by tandem duplication
<cit> .
although tomato  is one of the most consumed and cultivated field crops in the world, a large-scale phylogenetic analysis of tomato rlks has not yet been performed, and few members of the tomato rlk/pelle family  have been studied and characterized.
these members include pto
<cit> , pti <dig>
<cit> , and bti <dig>
<cit> , which interact with pseudomonas syringae elicitors; tark <dig>
<cit> , which interacts with the xanthomonas campestris elicitor; tpk1b
<cit> , whose expression is induced by mechanical wounding and oxidative stress; and sr <dig>
<cit> , which is the atbri <dig> ortholog and binds to systemin to respond to wounding or herbivore attack, although there is some debate about the function of this receptor
<cit> .
another well-studied rlk in tomato is nik , which interacts with nuclear shuttle protein  of geminivirus during infection
<cit> .
three homologs of nik in arabidopsis  have also been shown to interact with nsp through their kinase domain
<cit> .
this interaction causes inhibition of the kinase activity of niks and hence prevents the activation of the signal transduction cascade that evokes a plant defense response
<cit> .
these rlks are members of the lrrii subfamily that also contains the somatic embryogenesis receptor kinases
<cit> .
with the completion of the tomato genome sequencing along with the annotation of the encoded proteins
<cit> , it has become possible to study the rlk superfamily in this species using a large-scale phylogenetic approach.
according to genomic analyses, the tomato genome was predicted to have approximately  <dig> megabases of dna and encode  <dig>  proteins.
in this investigation, we identified and classified all putative tomato rlks by comparison with previously described arabidopsis rlks
<cit> .
we also showed that the tomato rlk members of lrrii subfamily, which comprises nik and serk genes, share similar biochemical activity , genomic structure and partial overlapping expression profiles with the arabidopsis orthologs.
our results provide a framework for understanding rlk function in tomato and reveal that some tomato and arabidopsis lrrii-rlk orthologs may play similar roles in antiviral defense.
RESULTS
the tomato rlk superfamily
the identification of the rlk superfamily members in tomato was initially performed by a batch blast search against a tomato protein database  using the kinase sequences of representative arabidopsis rlks as queries.
this analysis retrieved  <dig> tomato proteins that seemed to be rlks.
all of these retrieved tomato proteins were submitted for annotation of their domain structure using smart
<cit>   and pfam
<cit>   databases.
four proteins that did not bear a kinase domain were not considered for further analysis.
the remaining  <dig> proteins were used for phylogenetic analysis based on their kinase domain sequences.
for this analysis, we included all arabidopsis rlks to compare with tomato rlks and used representative proteins of other kinase families of arabidopsis and human as outgroups .
all arabidopsis rlks were placed in a major cluster together with  <dig> tomato proteins that were identified as members of the rlk superfamily .
the other  <dig> proteins were clustered with outgroups; consequently, they were not considered to be members of rlk superfamily .
the size of the tomato rlk superfamily  was similar to that of the arabidopsis rlk superfamily .
furthermore, almost all tomato rlks  were clustered with at least one arabidopsis rlk.
therefore, the tomato rlk superfamily was divided into the same  <dig> subfamilies as described previously for arabidopsis
<cit> .
as in arabidopsis, in which  <dig> out of all  <dig> rlks belong to leucine-rich repeat  subfamilies, tomato lrr subfamilies were the most abundant and contained  <dig> proteins.
another large rlk subfamily was rlck, which included  <dig> members in tomato, almost the same number as in arabidopsis .
among the  <dig> tomato rlks that were not clustered in the same branches as arabidopsis rlks, ten proteins were quite small and lacked a typical rlk structure, but the other six proteins had a clear rlk structure and as such were considered to be tomato-specific rlks.
among those six tomato-specific rlks, solyc03g <dig> contained a legume lectin domain similar to members of the lectin subfamily, and solyc02g <dig> harbored an amino oxidase domain , which is not found in any arabidopsis rlks.
the remaining four proteins did not have any predicted protein domains in their extracellular region.
rlk superfamily profiles in both species are summarized in figure
<dig>
although the rlk superfamilies of arabidopsis and tomato share common features, a close inspection reveals some interesting differences between them.
a comparison of membership size of each subfamily revealed some differences between the species.
to infer the number of rlks in their common ancestral and the occurrence of duplication/deletion events after the divergence of both species, we used a reconciliation method
<cit>  to compare the rlk superfamily tree  with a species tree generated at ncbi taxonomy browser .
then, the frequency of duplication or deletion events that occurred in each rlk subfamily was statistically analyzed.
we identified twelve rlk subfamilies from these plant species that have significantly expanded/reduced in their size .
in tomato, the lrk10l <dig>  sd <dig>  sd2b and lrrxii subfamilies displayed significant number of duplication, while in arabidopsis, expansion was observed in duf <dig>  l-lec, lrk10l <dig>  sd <dig>  wak, lrria and rlckxii/xiii subfamilies.
conversely, lysmii in arabidopsis and lrria and rlckxii in tomato reduced in size.
another distinct aspect of rlks in these plant species refers to the lack of a common domain structure in the extracellular region of rlk members of the lrk10l- <dig> subfamily.
whereas the extracellular region of arabidopsis lrk10l- <dig> members harbors diverse structures, such as thaumatin, glycerophosphoryl diester phosphodiesterase family  or malectin domains, the tomato rlks of this same subfamily do not contain any predicted domain structure in their n-terminal region.
further analyses were performed to predict function associated with expansion/reduction patterns.
as rlks are frequently associated to defense or developmental processes, we performed a search using the gene ontology  terms
<cit>  for functionally annotated rlks in those categories  and statistically compared the proportion of annotated genes in each rlk subfamily with the proportion of annotated genes in the whole rlk superfamily .
compelling evidence in the literature has demonstrated that defense-related genes have high duplication rate and are organized in tandem repeats
<cit> .
we also identified the rlks organized in tandem repeats and determined their frequency in each rlk subfamilies .
these analyses demonstrated that most of the rlk subfamilies, which expanded after arabidopsis and tomato species divergence, had their genes organized in tandem repeats.
functional annotation analysis of their genes revealed that the go terms of  <dig> out of  <dig> genes in tandem repeats were associated with defense response.
rlks annotated as developmental-related were overrepresented in crrlk1- <dig>  perk, lrrviia, lrrxi and lrrxiii subfamilies and all those subfamilies did not expand or present genes in tandem arrays.
the subfamilies that underwent reduction in their size were not associated with defense- or development-related functions, except for the lrria subfamily, which may be related to defense response.
consistent with the involvement of members of the lrrii subfamily in defense and development, the lrrii subfamily of tomato and arabidopsis had significantly high number of annotated members in either defense or developmental categories.
motif prediction, genomic structure and phylogenetic analysis of the lrrii subfamily
compelling evidence in the literature has revealed a fundamental role for members of the arabidopsis lrrii-rlk subfamily as co-receptors for transducing developmental and defense signals
<cit> .
the potential of the members of this subfamily as co-receptors involved in the activation of rlk-mediated signal transduction prompted us to perform a comprehensive analysis of the tomato lrrii-rlk subfamily to uncover related functions in tomatoes.
based on the phylogenetic tree of all members of rlk superfamily, the tomato lrrii-rlk subfamily encompassed  <dig> proteins.
the members of this group from both plant species have over  <dig> amino acids on average.
phylogenetic analysis of this group using full-length protein sequences resulted in a tree with three well-resolved clusters; the tomato and arabidopsis proteins were found in all clusters, although they had distinct sizes .
these clades were termed nik, serk and lrriic based on annotation of the arabidopsis members in each cluster.
the nik clade is formed by seven tomato members and six arabidopsis members, including the three atnik genes.
the serk clade clustered the five well-characterized serks in arabidopsis and three members of the tomato subfamily.
the lrriic clade consisted of three tomato proteins and three arabidopsis proteins whose functions are unknown.
motif prediction analysis on these proteins revealed that tomato and arabidopsis lrrii-rlk members display similar protein domains organized in the same fashion .
the consensus structural organization of the conserved domains between both species included an n-terminal signal peptide followed by a leucine zipper, five lrrs at the extracellular side and a transmembrane domain separating the n-terminal portion from the cytoplasmic c-terminal kinase domain.
among the serk genes of both species, there was also a proline-rich domain  localized between the last lrr and the transmembrane domain .
sequence alignment of the lrrii rlks showed several conserved amino acid positions among members of both plant species.
exon/intron boundaries were also well conserved.
variation at the sequence level was observed within the spp and signal peptide recognition domains.
among the lrrii clades, the proteins comprising the serk clade were more conserved with a larger number of conserved positions compared with the predicted proteins of the lrriic and nik clades .
genomic structure analysis revealed that in general lrrii genes are organized into  <dig> exons .
genes that varied in this number displayed fused exons, including atnik <dig>  at5g <dig>  and slnik <dig>  or had deleted exons, such as at5g <dig> , solyc02g <dig> .1 and solyc05g <dig> .
<dig>  intronic regions were larger in tomato members than in arabidopsis members and in serk genes compared with genes from other clades .
<cit> , respectively.
no expression data were obtained for atserk <dig> .
tissues with high mean-expression are summarized for each gene.
orthologous genes that had similar expression profiles are delimited by brackets.
expression analysis of lrrii subfamily genes in different tissues
we examined the expression profiles of lrrii-rlk genes in different tomato tissues, including leaf, stem, root, flower, cotyledon and hypocotyl, by real-time pcr.
the results are presented in figure
<dig> and summarized in figure
3b.
almost all analyzed genes exhibited expression in at least one organ, except for solyc05g <dig> .
<dig> and solyc02g <dig> .
<dig> in the nik clade, which had very low expression in all tissues tested.
the tomato genes in the serk clade were expressed more highly in leaf and cotyledon tissue than in stem or flower tissue.
the nik clade genes were highly expressed in diverse organs, such as leaves, flowers and roots.
the lrriic group encompassed genes with a higher level of expression in cotyledon, flower and leaf tissues and with lower expression in stem tissue.
the expression data for lrrii subfamily members from arabidopsis  were extracted from atgenexpress
<cit>  to examine whether there is some correspondence in the expression profiles between orthologous pairs of tomato and arabidopsis genes.
statistical analyses of correlation between tomato and arabidopsis expression data could not be performed because these data have been generated by different methods  and hence they have different units.
nevertheless, a subjective comparison of the expression analysis from both plants revealed that the majority of the orthologous genes displayed partially but not entirely overlapping expression profiles .
the orthologous groups that presented similar expression profile were atserk1/atserk <dig> and solyc04g <dig>  which had high expression levels in leaves, atnik <dig> and slnik <dig>  which were lowly expressed in stem and cotyledon tissues, and atnik <dig> and slnik <dig>  which were most highly expressed in the leaf.
interactions between representatives of the lrrii subfamily and nsp of toysv
we have previously shown that nsp from begomovirus interacts with members of the lrrii-rlk subfamily, such as atnik <dig>  atnik <dig> and atnik <dig>  to suppress host defense, and it interacts with a member of the perk-like rlk subfamily, nsp-associated kinase , to potentiate virus infection
<cit> .
either nik from tomato and nsak from arabidopsis were isolated by their capacity to interact with nsp through two-hybrid screening
<cit> .
the nsp interactions with the arabidopsis atnik <dig>  atnik <dig> and atnik <dig> and nsak were further confirmed by yeast-two hybrid and in vitro pull down assays
<cit> .
we have recently shown by bimolecular fluorescence complementation  assay that nsp also interacts with nik in vivo.
because begomovirus negatively impacts tomato cultivation worldwide, we selected representatives of tomato rlks from the lrrii subfamily and examined their capability to interact with nsp of tomato yellow spot virus  similar to the interaction observed with arabidopsis nik genes
<cit> .
yeast two-hybrid experiments were performed using the toysv-nsp  as prey and kinase domains of slnik <dig> , slnik <dig> , slnik <dig>  and slbak <dig>  as bait.
we also analyzed a perk representative  that is similar to the nsp-interactor perk-like gene of arabidopsis .
a tomato gene  from the lrrxii subfamily, homolog of the arabidopsis ef-tu receptor , was used as a negative control.
interactions between the viral nsp and host proteins were detected after co-transforming the yeast cells with both bait and prey plasmids and monitoring for histidine prototrophy.
nsp was found to interact with the kinase domains of slnik <dig>  slnik <dig>  slnik <dig> and slbak <dig> or with the kinase domain of the perk representative slnsak .
the nsp interactions were specific to the tomato lrrii-rlk orthologs and to perk-like slnsak because the his marker gene was not activated in yeast cells co-transformed with tynsp-p <dig>  and with either the empty vector or slefr-p <dig> .
furthermore, co-transformation of yeast with the nsp interactors fused to the gal4-binding domain and the empty vector expressing the gal4-activating domain alone also failed to activate the his marker gene .
these results expanded our previous observation that nsp-nik complex formation was neither virus-specific nor host-specific
<cit> .
they also suggest that slnsak is a nsp target during begomovirus infection in tomato.
certainly, the in vivo demonstration of these interactions will further support these interpretations.
discussion
the structure of the tomato rlk superfamily and the proposed evolution of rlk superfamily in plants
to date, the phylogenetic and structural characterization of the rlk superfamily has been limited to the following plant species: moss, rice, poplar, soybean and arabidopsis
<cit> .
the size of these families ranges from  <dig> to  <dig> proteins , and their extracellular regions bear a great variety of protein domain structures.
in the present investigation, we characterized and generated a complete list of the tomato rlk superfamily members , identifying  <dig> rlks, which falls in the size range of the arabidopsis  and soybean  superfamilies.
evolutionary analyses of the rlk superfamily has suggested that the rlk structure was established prior to the divergence of land plants from algae because proteins with rlk configurations were discovered in the unicellular algae chlamydomonas reinhardtii <cit> .
comparative analysis of rlks among moss, rice, poplar and arabidopsis revealed that the rlk superfamily underwent expansion in the beginning of the land plant lineage, after the divergence of angiosperm and bryophyte and independently during diversification of each angiosperm lineage.
the most dramatic expansion was observed in the rice and poplar lineages, which have almost twice as many as rlk members as arabidopsis
<cit> .
this evolutionary scenario has not been changed by inclusion of data regarding the soybean
<cit>  and tomato  rlk superfamily expansion .
the phylogenetic tree of the members of the rlk superfamilies in tomato and arabidopsis revealed that most of the rlk subfamilies have maintained approximately the same number of rlk members between these species.
exceptions were observed for the duf <dig>  l-lec, lrk10l <dig>  sd <dig>  sd2b, wak, lrria, lrrxii, lrrxiib, rlckxii/xiii subfamilies, in which specific and extensive expansion was observed in one of the two plant species, as well as for the lysmii, lrria and rlckxii, in which specific reduction was observed .
functional annotations of some arabidopsis rlks and the number of genes in tandem repeat that compose those subfamilies indicated a predominance of genes clustered in tandem array and defense-related rlks .
some of those subfamilies, such as duf <dig>  l-lec, sd <dig>  sd2b, wak and lrria, had been shown previously to be overrepresented in microarray analysis of arabidopsis under different stress conditions
<cit> .
taken together, these results are consistent with the previous assumption that specific expansions of rlk genes have occurred more frequently for those rlks associated with defense response.
specific expansion of tomato rlks compared to arabidopsis occurred in the lrk10l <dig>  lrrxii, sd <dig> and sd2b subfamilies.
interestingly, all of those rlk subfamilies, except for sd <dig>  were also overrepresented in rice when compared with arabidopsis
<cit> .
the lrrxii subfamily comprises the ef-tu receptor  and flagellin-sensitive  <dig>  in arabidopsis, and xa <dig> in rice, all of which are associated with defense responses.
the expansion of the lrrxii subfamily in cultivated plants such as tomato and rice has previously been suggested to be associated with the accumulation of resistance genes by intense breeding programs
<cit> .
likewise, the sd subfamily also has representatives involved in defense response, such as rk <dig> and rk3
<cit> , but this subfamily is also strongly associated with the self-incompatibility process
<cit> .
arabidopsis is known to be self-compatible, while rice, tomato and the close relative of a. thaliana, a. lyrata, are self-incompatible.
thus, the specific deletion of some representatives in a. thaliana could have contributed to the generation of the self-compatible mode in this species.
in contrast, some rlk subfamilies were expanded specifically in arabidopsis.
those include duf <dig>  l-lec, lrk10l <dig>  sd <dig>  wak, lrria and rlckxii/xiii.
except for sd <dig> and lrk10l <dig> subfamilies, no significant expansion were observed in these rlk subfamilies in tomato, although all of them contain rlks that are involved in defense response, such as fls2-induced receptor-kinase  <dig>
<cit> , lectin receptor kinase  <dig>
<cit> , pathogenesis-related 5-like receptor kinase
<cit> , resistant to fusarium oxysporium  <dig>
<cit>  and cysteine-rich rlk  <dig>
<cit> .
another relevant distinction between tomato and arabidopsis rlk subfamilies derives from the diversification of the extracellular domain patterns of the lrk10l <dig> subfamily representatives.
arabidopsis members of the lrk10l <dig> subfamily have unique domain structures, such as gdpd, thaumatin and malectin domains, while the tomato members do not possess any characterized domain structure in their extracellular region.
likewise, these domain structures have not been found in rice, poplar or moss rlks, indicating that within the rlk superfamily they are specific to arabidopsis.
we also identified a tomato-specific rlk that possesses an amino oxidase domain in its extracellular region.
these rlks may respond to molecular signals not perceived by other plants.
although gaining a novel protein structure could increase the repertoire of signals perceived by plants, the small number of lineage-specific rlks in tomato, as was also reported in rice and poplar, further substantiates the hypothesis that the expansion of existing rlk kinase subfamilies is the major mechanism of evolution of these proteins.
the members of the rlk superfamily are involved in diverse biological processes at all steps of plant development.
thus, the gain or loss of a rlk gene could have serious repercussions on plant phenotype.
the specific profiles of the rlk superfamily found in tomato and arabidopsis are certainly responsible for several differences between these plants, such as morphology, reproduction and, importantly, responsiveness to different stress conditions.
among tomato and arabidopsis, few rlk subfamilies have undergone specific expansion or reduction after their speciation.
this scarcity may indicate that variation in rlk superfamily profiles in both plants appeared recently.
the domestication process that tomatoes underwent could have been a significant factor contributing to this variation because some rlks have been directly linked to traits targeted by artificial selection, such as disease resistance and growth.
nevertheless, most of the specific expansions of rlk subfamilies observed in arabidopsis were also associated with defense response and had occurred at similar rates .
these data suggest that, independently of artificial selection, arabidopsis had also expanded and developed a specific machinery against abiotic or biotic stress response, which argues against the assumption that artificial selection leads to resistant genes accumulation.
to further examine the influence of artificial selection on the repertoire of plant rlks, new genetic resources for closely related wild plants are necessary.
functional expression analysis of the lrr ii subfamily members in tomato and arabidopsis
the lrrii subfamily contains rlks with dual functions in development and defense response
<cit> .
characterized members of this subfamily include  serk genes, which are associated with diverse processes, such as brassinosteroid signaling, flagellin, cell death, light and pathogen-associated molecular pattern  responses
<cit> , and  nik genes, which interact directly with geminivirus nsp during viral infection
<cit> .
phylogenetic and protein structure analyses on lrrii subfamily members of tomato and arabidopsis demonstrated that this group is highly conserved between these species.
in rice, in which the rlk superfamily has undergone a large expansion, the lrrii subfamily members are also conserved in number and sequence, indicating that biochemical pathways regulated by lrrii-rlks have essential and conserved roles in angiosperm species.
although lrrii members have well-conserved amino acid sequences among various species
<cit> , expression analysis of the members of the tomato and arabidopsis lrrii subfamilies demonstrated that only a few of the orthologous pairs resemble in their expression profiles.
by analogy with some evidence in the literature from other plant species, one may envision that these orthologous genes could have functionally diverged after the speciation event separating tomato and arabidopsis.
functional divergence in orthologous genes is not an uncommon event in both plants
<cit>  and animals
<cit> .
for example, the crabs claw transcription factor in arabidopsis is expressed in the carpel primordial abaxial region and in floral nectarines and regulates carpel morphology and nectar development, whereas its orthologous in rice, drooping leaf , is expressed in the whole carpel primordium and in central undifferentiated cells of leaves, where it regulates carpel identity and midrib development
<cit> .
the expression of orthologous genes has also been shown to vary differently in response to a stress condition.
in barley, which is tolerant to salinity, the expression of genes involved in root development, such as constans-like  <dig> , is suppressed by high salinity, whereas the expression of the rice orthologous is unchanged under the same stress condition
<cit> .
likewise, a large fraction of orthologous pairs of rice and arabidopsis genes with receptor activity do not display conserved co-expression
<cit> .
therefore, different patterns of expression between the orthologous genes in the lrrii subfamilies from tomato and arabidopsis may be a result of functional divergence that occurred between these genes.
functional divergence in receptor proteins with a developmental function may lead to a dramatic change in the plant phenotype because plant development is heavily guided by external signals.
for example, a tissue that displays high expression of certain rlks is likely to be more sensitive to perception of rlk-specific sensing signals, leading to a rapid and effective response.
in contrast, reduced expression of an orthologous gene from a different species in the same tissue would decrease the effectiveness and delay the signal perception and response.
this difference in the cell responsiveness to a specific signal could represent the differential timing of biochemical reactions that are regulated by this signal.
in developmental process, small differences in reaction time may be sufficient to generate a distinct phenotype in the plant.
in contrast, the orthologous pairs serk1/serk2/solyc04g <dig> .
<dig>  which display similar expression profiles , also contain the most conserved extracellular and intracellular domains .
the other two orthologous pairs, nik1/slnik <dig> and nik3/slnik <dig>  which also displayed similar expression profiles, also had highly conserved extracellular and intracellular regions.
in both orthologous pairs, sequence identity was approximately 65% in the extracellular regions and approximately 80% in the intracellular regions .
this finding may indicate a tight conservation of function between these members of the arabidopsis and tomato lrrii rlk subfamilies.
conservation of geminivirus interactions with members of the rlk family in tomato
although most of lrrii subfamily orthologous pairs exhibited functional divergence, we showed that the tomato orthologs of the lrrii-rlks members nik <dig>  nik <dig> and nik <dig> retain the capacity to interact with geminivirus nsp in yeast
<cit> .
at least for the nik <dig> and nik <dig> ortholog pairs, the functional conservation associated with specific protein-protein interactions may be linked to the high conservation of their nsp-interacting kinase domain  and similarity of expression profiles .
the current model of nik-mediated defense response posits that the immune receptor protects plant against geminiviruses by phosphorylating the ribosomal protein l <dig>
<cit> .
phosphorylation of rpl <dig> by nik redirects the ribosomal protein to the nucleus, where it may mount a defense mechanism to prevent viral proliferation.
during geminivirus infection, nsp interacts with the kinase domain of niks to inhibit their kinase activity, preventing activation of the defense response.
despite the high similarity between nik genes and serk genes, atbak1/serk <dig> and atserk <dig> do not functionally replace the atnik <dig> role in transducing an antiviral signaling response and do not interact with the viral nsp
<cit> .
in contrast, we found that the atbak <dig> ortholog from tomato interacts with nsp in yeast.
although the functional relevance of this interaction in planta remains to be determined, it is worth noting that the expression profiles of the bak <dig> orthologs are not similar, as would be expected for functionally divergent orthologs.
although both orthologs are ubiquitously expressed in the cognate plant species, they are expressed to different extents in distinct organs.
whereas atbak <dig> expression is quantitatively similar and relatively low in all organs analyzed, its ortholog from tomato displays a higher level of expression in the cotyledons, hypocotyls and leaves, where geminivirus infection largely takes place.
therefore, the expression profiles of the nsp interactors  seemed to be related with the onset of geminivirus infection.
due to the high expression of the atbak <dig> tomato orthologous in leaves, one may envision the existence of evolving selective pressures to diverge the corresponding nsp-interacting domains of the bak <dig> orthologs towards functional fitness with regard to geminivirus infection.
in contrast to nik receptors, which are inhibited by nsp interaction, atnsak, a member of perk subfamily, interacts with nsp and phosphorylates the viral protein in vitro <cit> .
loss of nsak function enhances tolerance to geminivirus infection, indicating that atnsak is a positive contributor to geminivirus infection in arabidopsis.
here, we showed that the nsak tomato orthologous retains its capacity to interact with viral nsp.
this demonstrates that specific members of the rlk family have conserved defense functions  or compatibility functions  in response to viral infection.
due to the emergence of new species of tomato-infecting begomoviruses that rapidly evolve through recombination or pseudo-recombination to produce divergent genome sequences that gives the virus an advantage over its host’s recognition system, a survey of the interactions between nsps from distinct tomato-infecting geminiviruses and slniks and slnsak may add insights into the co-evolution of the viral protein and host defense/compatibility functions.
CONCLUSIONS
the rlk superfamily is a large and diverse group of transmembrane receptors that enables plants to perceive a diverse array of signals at the cell surface, creating an efficient mechanism for cell-environment communication.
in this investigation, we generated a complete list of the members of the tomato rlk superfamily, which is made-up of  <dig> proteins.
the tomato rlk sequences exhibited a typical receptor-like kinase configuration and almost all of them were phylogenetically clustered with at least one member of the arabidopsis rlk superfamily.
therefore, the tomato rlk superfamily is similarly organized, with the same number and identity of subfamilies as previously defined for arabidopsis rlks.
among the  <dig> rlk subfamilies, twelve showed specific and extensive expansion or reduction in the number of their rlk members, which may be a reflection of lineage-specific responses to various biotic and abiotic stresses.
the intense breeding programs tomatoes have been subjected to may also have contributed to the establishment of the current rlk superfamily profile in this species.
this comprehensive analysis comparing the complete repertory of arabidopsis and tomato rlks may provide a framework to rationalize future functional studies of the members of this family.
phylogenetic and structural analyses of lrrii subfamily members from both tomato and arabidopsis reveal a well-conserved group both in terms of sequence and protein domain organization.
as a consequence, the tomato lrrii-rlk subfamily is organized into the same three with phylogenetically supported clades, serk, nik and lrriic clusters.
nevertheless, a comparison of the expression between orthologous genes of this subfamily demonstrated that the majority of the orthologous pairs did not share a similar expression profile, indicating that these orthologous lrrii-rlks may have undergone functional divergence.
this finding is supported by the observation that, in contrast to the arabidopsis atbak <dig>  slbak <dig> interacts with the geminivirus nsp and is highly expressed in leaves and the cotyledon.
this pattern of slbak <dig> expression is consistent with the pattern of infection by tomato-infecting begomoviruses, which infect leaf tissues and move through the phloem but do not invade roots.
additionally, as immune receptors, the orthologous pairs nik <dig> and nik <dig> displayed both the capacity to interact with the begomovirus virulence factor nsp and expression profiles that parallel the onset of begomovirus infection.
evidence for functional conservation between nik <dig> orthologs has been previously provided with the demonstration that nik <dig> from arabidopsis is capable of protecting tomato plants against tomato-infecting begomovirus
<cit> .
collectively, our results indicate that nik orthologs retain similar functions as defense receptors to protect plant cells against viral attack.
therefore, nik-mediated antiviral signaling likely also operates in tomato, suggesting that the tomato niks may be good candidate targets for engineering resistance against tomato-infecting begomoviruses.

 plant cells constantly react to multiple signals that come from the local environment, neighboring cells, or even from other organisms. depending on the stimuli, plant cells may expand, divide, differentiate, synthesize compounds, prepare against pathogen infection, or induce necrosis. to perceive and receive these signals, plant cells possess complex systems of transmembrane receptor proteins that facilitate communication between the intracellular environment and the outside world. one of the largest groups of these receptors is the receptor-like kinase superfamily, which contains over members in arabidopsis. rlks are structurally organized into an extracellular domain that can be highly divergent, followed by a transmembrane segment and a conserved intracellular serine/threonine kinase domain. most rlks are localized in the plasma membrane, although there are also rlk members that are found in the cytoplasm. in this case, rlks do not possess either an extracellular region or a transmembrane domain and are called receptor-like cytoplasmic kinases. analyses of arabidopsis rlks by structural comparison of their extracellular region and phylogenetic analysis of their kinase domain revealed that they can be divided into over subfamilies. several distinct rlks have been studied in the past decade, and a common theme that has emerged is that binding of a specific signal molecule to their extracellular domain is required to initiate a signal transduction cascade. generally, ligand-receptor interactions at the extracellular domain of rlks initiate the propagation of the signal through the membrane by inducing a conformational change in the receptor kinase domain, which allows interactions with other rlks resulting in homo- or heterodimers. dimerized rlks are then transphosphorylated by their cytoplasmic kinase domain, leading to both activation of the kinase and establishment of docking sites for phosphorylation of downstream phosphorylation targets. this activation mechanism of plant rlks is similar to that of signal transduction mediated by receptor tyrosine kinases in animal cells, which share a common origin with plant serine/threonine kinases. functional analysis of rlks indicates that the majority of them are associated with plant development or defense response, but there are also rlks involved in cell wall attachment, plant-bacterial symbiotic interactions and self-incompatibility. among all rlks, those bearing a leucine-rich repeat domains are overrepresented in the rlk superfamily, comprising over 38% of arabidopsis rlks, which are distributed into subgroups. the lrr domains in these receptors vary in number and in the distribution pattern of the lrrs along the extracellular region. examples of well-known lrr-rlks include clavata which controls the size of stem cells in the apical meristem by forming a heterodimer with clavata and then interacting with clavata through the extracellular domain, and brassinosteroid insensitive-, which perceives brassinosteroids and interacts with it receptor partner, bak. other functions associated with lrr-rlks include morphogenesis, embryogenesis, pollen self-incompatibility and responses to environmental signals. in addition, some lrr-rlks are known to function as regulators of defense response to bacterial pathogen, necrotrophic fungus and viral infection. most of the characterized rlks are from model plants such as arabidopsis and medicago truncatula, but significant efforts have been made to expand these studies to relevant field crops. large-scale comparative analyses of arabidopsis rlks with rice and soybean rlks identified over kinase proteins in rice and in soybean belonging to the rlk superfamily; almost all members were grouped into previously determined arabidopsis rlk subfamilies. the rlk subfamilies with developmental function have conserved size, whereas those involved in defense response have expanded their members, mainly by tandem duplication. although tomato is one of the most consumed and cultivated field crops in the world, a large-scale phylogenetic analysis of tomato rlks has not yet been performed, and few members of the tomato rlk/pelle family have been studied and characterized. these members include pto, pti, and bti, which interact with pseudomonas syringae elicitors; tark, which interacts with the xanthomonas campestris elicitor; tpk1b, whose expression is induced by mechanical wounding and oxidative stress; and sr, which is the atbri ortholog and binds to systemin to respond to wounding or herbivore attack, although there is some debate about the function of this receptor. another well-studied rlk in tomato is nik, which interacts with nuclear shuttle protein of geminivirus during infection. three homologs of nik in arabidopsis have also been shown to interact with nsp through their kinase domain. this interaction causes inhibition of the kinase activity of niks and hence prevents the activation of the signal transduction cascade that evokes a plant defense response. these rlks are members of the lrrii subfamily that also contains the somatic embryogenesis receptor kinases. with the completion of the tomato genome sequencing along with the annotation of the encoded proteins, it has become possible to study the rlk superfamily in this species using a large-scale phylogenetic approach. according to genomic analyses, the tomato genome was predicted to have approximately megabases of dna and encode proteins. in this investigation, we identified and classified all putative tomato rlks by comparison with previously described arabidopsis rlks. we also showed that the tomato rlk members of lrrii subfamily, which comprises nik and serk genes, share similar biochemical activity, genomic structure and partial overlapping expression profiles with the arabidopsis orthologs. our results provide a framework for understanding rlk function in tomato and reveal that some tomato and arabidopsis lrrii-rlk orthologs may play similar roles in antiviral defense.  the tomato rlk superfamily the identification of the rlk superfamily members in tomato was initially performed by a batch blast search against a tomato protein database using the kinase sequences of representative arabidopsis rlks as queries. this analysis retrieved tomato proteins that seemed to be rlks. all of these retrieved tomato proteins were submitted for annotation of their domain structure using smart and pfam databases. four proteins that did not bear a kinase domain were not considered for further analysis. the remaining proteins were used for phylogenetic analysis based on their kinase domain sequences. for this analysis, we included all arabidopsis rlks to compare with tomato rlks and used representative proteins of other kinase families of arabidopsis and human as outgroups. all arabidopsis rlks were placed in a major cluster together with tomato proteins that were identified as members of the rlk superfamily. the other proteins were clustered with outgroups; consequently, they were not considered to be members of rlk superfamily. the size of the tomato rlk superfamily was similar to that of the arabidopsis rlk superfamily. furthermore, almost all tomato rlks were clustered with at least one arabidopsis rlk. therefore, the tomato rlk superfamily was divided into the same subfamilies as described previously for arabidopsis. as in arabidopsis, in which out of all rlks belong to leucine-rich repeat subfamilies, tomato lrr subfamilies were the most abundant and contained proteins. another large rlk subfamily was rlck, which included members in tomato, almost the same number as in arabidopsis. among the tomato rlks that were not clustered in the same branches as arabidopsis rlks, ten proteins were quite small and lacked a typical rlk structure, but the other six proteins had a clear rlk structure and as such were considered to be tomato-specific rlks. among those six tomato-specific rlks, solyc03g contained a legume lectin domain similar to members of the lectin subfamily, and solyc02g harbored an amino oxidase domain, which is not found in any arabidopsis rlks. the remaining four proteins did not have any predicted protein domains in their extracellular region. rlk superfamily profiles in both species are summarized in figure  although the rlk superfamilies of arabidopsis and tomato share common features, a close inspection reveals some interesting differences between them. a comparison of membership size of each subfamily revealed some differences between the species. to infer the number of rlks in their common ancestral and the occurrence of duplication/deletion events after the divergence of both species, we used a reconciliation method to compare the rlk superfamily tree with a species tree generated at ncbi taxonomy browser. then, the frequency of duplication or deletion events that occurred in each rlk subfamily was statistically analyzed. we identified twelve rlk subfamilies from these plant species that have significantly expanded/reduced in their size. in tomato, the lrk10l sd sd2b and lrrxii subfamilies displayed significant number of duplication, while in arabidopsis, expansion was observed in duf l-lec, lrk10l sd wak, lrria and rlckxii/xiii subfamilies. conversely, lysmii in arabidopsis and lrria and rlckxii in tomato reduced in size. another distinct aspect of rlks in these plant species refers to the lack of a common domain structure in the extracellular region of rlk members of the lrk10l- subfamily. whereas the extracellular region of arabidopsis lrk10l- members harbors diverse structures, such as thaumatin, glycerophosphoryl diester phosphodiesterase family or malectin domains, the tomato rlks of this same subfamily do not contain any predicted domain structure in their n-terminal region. further analyses were performed to predict function associated with expansion/reduction patterns. as rlks are frequently associated to defense or developmental processes, we performed a search using the gene ontology terms for functionally annotated rlks in those categories and statistically compared the proportion of annotated genes in each rlk subfamily with the proportion of annotated genes in the whole rlk superfamily. compelling evidence in the literature has demonstrated that defense-related genes have high duplication rate and are organized in tandem repeats. we also identified the rlks organized in tandem repeats and determined their frequency in each rlk subfamilies. these analyses demonstrated that most of the rlk subfamilies, which expanded after arabidopsis and tomato species divergence, had their genes organized in tandem repeats. functional annotation analysis of their genes revealed that the go terms of out of genes in tandem repeats were associated with defense response. rlks annotated as developmental-related were overrepresented in crrlk1- perk, lrrviia, lrrxi and lrrxiii subfamilies and all those subfamilies did not expand or present genes in tandem arrays. the subfamilies that underwent reduction in their size were not associated with defense- or development-related functions, except for the lrria subfamily, which may be related to defense response. consistent with the involvement of members of the lrrii subfamily in defense and development, the lrrii subfamily of tomato and arabidopsis had significantly high number of annotated members in either defense or developmental categories. motif prediction, genomic structure and phylogenetic analysis of the lrrii subfamily compelling evidence in the literature has revealed a fundamental role for members of the arabidopsis lrrii-rlk subfamily as co-receptors for transducing developmental and defense signals. the potential of the members of this subfamily as co-receptors involved in the activation of rlk-mediated signal transduction prompted us to perform a comprehensive analysis of the tomato lrrii-rlk subfamily to uncover related functions in tomatoes. based on the phylogenetic tree of all members of rlk superfamily, the tomato lrrii-rlk subfamily encompassed proteins. the members of this group from both plant species have over amino acids on average. phylogenetic analysis of this group using full-length protein sequences resulted in a tree with three well-resolved clusters; the tomato and arabidopsis proteins were found in all clusters, although they had distinct sizes. these clades were termed nik, serk and lrriic based on annotation of the arabidopsis members in each cluster. the nik clade is formed by seven tomato members and six arabidopsis members, including the three atnik genes. the serk clade clustered the five well-characterized serks in arabidopsis and three members of the tomato subfamily. the lrriic clade consisted of three tomato proteins and three arabidopsis proteins whose functions are unknown. motif prediction analysis on these proteins revealed that tomato and arabidopsis lrrii-rlk members display similar protein domains organized in the same fashion. the consensus structural organization of the conserved domains between both species included an n-terminal signal peptide followed by a leucine zipper, five lrrs at the extracellular side and a transmembrane domain separating the n-terminal portion from the cytoplasmic c-terminal kinase domain. among the serk genes of both species, there was also a proline-rich domain localized between the last lrr and the transmembrane domain. sequence alignment of the lrrii rlks showed several conserved amino acid positions among members of both plant species. exon/intron boundaries were also well conserved. variation at the sequence level was observed within the spp and signal peptide recognition domains. among the lrrii clades, the proteins comprising the serk clade were more conserved with a larger number of conserved positions compared with the predicted proteins of the lrriic and nik clades. genomic structure analysis revealed that in general lrrii genes are organized into exons. genes that varied in this number displayed fused exons, including atnik at5g and slnik or had deleted exons, such as at5g, solyc02g.1 and solyc05g. intronic regions were larger in tomato members than in arabidopsis members and in serk genes compared with genes from other clades., respectively. no expression data were obtained for atserk. tissues with high mean-expression are summarized for each gene. orthologous genes that had similar expression profiles are delimited by brackets. expression analysis of lrrii subfamily genes in different tissues we examined the expression profiles of lrrii-rlk genes in different tomato tissues, including leaf, stem, root, flower, cotyledon and hypocotyl, by real-time pcr. the results are presented in figure  and summarized in figure  3b. almost all analyzed genes exhibited expression in at least one organ, except for solyc05g. and solyc02g. in the nik clade, which had very low expression in all tissues tested. the tomato genes in the serk clade were expressed more highly in leaf and cotyledon tissue than in stem or flower tissue. the nik clade genes were highly expressed in diverse organs, such as leaves, flowers and roots. the lrriic group encompassed genes with a higher level of expression in cotyledon, flower and leaf tissues and with lower expression in stem tissue. the expression data for lrrii subfamily members from arabidopsis were extracted from atgenexpress to examine whether there is some correspondence in the expression profiles between orthologous pairs of tomato and arabidopsis genes. statistical analyses of correlation between tomato and arabidopsis expression data could not be performed because these data have been generated by different methods and hence they have different units. nevertheless, a subjective comparison of the expression analysis from both plants revealed that the majority of the orthologous genes displayed partially but not entirely overlapping expression profiles. the orthologous groups that presented similar expression profile were atserk1/atserk and solyc04g which had high expression levels in leaves, atnik and slnik which were lowly expressed in stem and cotyledon tissues, and atnik and slnik which were most highly expressed in the leaf. interactions between representatives of the lrrii subfamily and nsp of toysv we have previously shown that nsp from begomovirus interacts with members of the lrrii-rlk subfamily, such as atnik atnik and atnik to suppress host defense, and it interacts with a member of the perk-like rlk subfamily, nsp-associated kinase, to potentiate virus infection. either nik from tomato and nsak from arabidopsis were isolated by their capacity to interact with nsp through two-hybrid screening. the nsp interactions with the arabidopsis atnik atnik and atnik and nsak were further confirmed by yeast-two hybrid and in vitro pull down assays. we have recently shown by bimolecular fluorescence complementation assay that nsp also interacts with nik in vivo. because begomovirus negatively impacts tomato cultivation worldwide, we selected representatives of tomato rlks from the lrrii subfamily and examined their capability to interact with nsp of tomato yellow spot virus similar to the interaction observed with arabidopsis nik genes. yeast two-hybrid experiments were performed using the toysv-nsp as prey and kinase domains of slnik, slnik, slnik and slbak as bait. we also analyzed a perk representative that is similar to the nsp-interactor perk-like gene of arabidopsis. a tomato gene from the lrrxii subfamily, homolog of the arabidopsis ef-tu receptor, was used as a negative control. interactions between the viral nsp and host proteins were detected after co-transforming the yeast cells with both bait and prey plasmids and monitoring for histidine prototrophy. nsp was found to interact with the kinase domains of slnik slnik slnik and slbak or with the kinase domain of the perk representative slnsak. the nsp interactions were specific to the tomato lrrii-rlk orthologs and to perk-like slnsak because the his marker gene was not activated in yeast cells co-transformed with tynsp-p and with either the empty vector or slefr-p. furthermore, co-transformation of yeast with the nsp interactors fused to the gal4-binding domain and the empty vector expressing the gal4-activating domain alone also failed to activate the his marker gene. these results expanded our previous observation that nsp-nik complex formation was neither virus-specific nor host-specific. they also suggest that slnsak is a nsp target during begomovirus infection in tomato. certainly, the in vivo demonstration of these interactions will further support these interpretations. discussion the structure of the tomato rlk superfamily and the proposed evolution of rlk superfamily in plants to date, the phylogenetic and structural characterization of the rlk superfamily has been limited to the following plant species: moss, rice, poplar, soybean and arabidopsis. the size of these families ranges from to proteins, and their extracellular regions bear a great variety of protein domain structures. in the present investigation, we characterized and generated a complete list of the tomato rlk superfamily members, identifying rlks, which falls in the size range of the arabidopsis and soybean superfamilies. evolutionary analyses of the rlk superfamily has suggested that the rlk structure was established prior to the divergence of land plants from algae because proteins with rlk configurations were discovered in the unicellular algae chlamydomonas reinhardtii. comparative analysis of rlks among moss, rice, poplar and arabidopsis revealed that the rlk superfamily underwent expansion in the beginning of the land plant lineage, after the divergence of angiosperm and bryophyte and independently during diversification of each angiosperm lineage. the most dramatic expansion was observed in the rice and poplar lineages, which have almost twice as many as rlk members as arabidopsis. this evolutionary scenario has not been changed by inclusion of data regarding the soybean and tomato rlk superfamily expansion. the phylogenetic tree of the members of the rlk superfamilies in tomato and arabidopsis revealed that most of the rlk subfamilies have maintained approximately the same number of rlk members between these species. exceptions were observed for the duf l-lec, lrk10l sd sd2b, wak, lrria, lrrxii, lrrxiib, rlckxii/xiii subfamilies, in which specific and extensive expansion was observed in one of the two plant species, as well as for the lysmii, lrria and rlckxii, in which specific reduction was observed. functional annotations of some arabidopsis rlks and the number of genes in tandem repeat that compose those subfamilies indicated a predominance of genes clustered in tandem array and defense-related rlks. some of those subfamilies, such as duf l-lec, sd sd2b, wak and lrria, had been shown previously to be overrepresented in microarray analysis of arabidopsis under different stress conditions. taken together, these results are consistent with the previous assumption that specific expansions of rlk genes have occurred more frequently for those rlks associated with defense response. specific expansion of tomato rlks compared to arabidopsis occurred in the lrk10l lrrxii, sd and sd2b subfamilies. interestingly, all of those rlk subfamilies, except for sd were also overrepresented in rice when compared with arabidopsis. the lrrxii subfamily comprises the ef-tu receptor and flagellin-sensitive in arabidopsis, and xa in rice, all of which are associated with defense responses. the expansion of the lrrxii subfamily in cultivated plants such as tomato and rice has previously been suggested to be associated with the accumulation of resistance genes by intense breeding programs. likewise, the sd subfamily also has representatives involved in defense response, such as rk and rk3, but this subfamily is also strongly associated with the self-incompatibility process. arabidopsis is known to be self-compatible, while rice, tomato and the close relative of a. thaliana, a. lyrata, are self-incompatible. thus, the specific deletion of some representatives in a. thaliana could have contributed to the generation of the self-compatible mode in this species. in contrast, some rlk subfamilies were expanded specifically in arabidopsis. those include duf l-lec, lrk10l sd wak, lrria and rlckxii/xiii. except for sd and lrk10l subfamilies, no significant expansion were observed in these rlk subfamilies in tomato, although all of them contain rlks that are involved in defense response, such as fls2-induced receptor-kinase, lectin receptor kinase, pathogenesis-related 5-like receptor kinase, resistant to fusarium oxysporium and cysteine-rich rlk. another relevant distinction between tomato and arabidopsis rlk subfamilies derives from the diversification of the extracellular domain patterns of the lrk10l subfamily representatives. arabidopsis members of the lrk10l subfamily have unique domain structures, such as gdpd, thaumatin and malectin domains, while the tomato members do not possess any characterized domain structure in their extracellular region. likewise, these domain structures have not been found in rice, poplar or moss rlks, indicating that within the rlk superfamily they are specific to arabidopsis. we also identified a tomato-specific rlk that possesses an amino oxidase domain in its extracellular region. these rlks may respond to molecular signals not perceived by other plants. although gaining a novel protein structure could increase the repertoire of signals perceived by plants, the small number of lineage-specific rlks in tomato, as was also reported in rice and poplar, further substantiates the hypothesis that the expansion of existing rlk kinase subfamilies is the major mechanism of evolution of these proteins. the members of the rlk superfamily are involved in diverse biological processes at all steps of plant development. thus, the gain or loss of a rlk gene could have serious repercussions on plant phenotype. the specific profiles of the rlk superfamily found in tomato and arabidopsis are certainly responsible for several differences between these plants, such as morphology, reproduction and, importantly, responsiveness to different stress conditions. among tomato and arabidopsis, few rlk subfamilies have undergone specific expansion or reduction after their speciation. this scarcity may indicate that variation in rlk superfamily profiles in both plants appeared recently. the domestication process that tomatoes underwent could have been a significant factor contributing to this variation because some rlks have been directly linked to traits targeted by artificial selection, such as disease resistance and growth. nevertheless, most of the specific expansions of rlk subfamilies observed in arabidopsis were also associated with defense response and had occurred at similar rates. these data suggest that, independently of artificial selection, arabidopsis had also expanded and developed a specific machinery against abiotic or biotic stress response, which argues against the assumption that artificial selection leads to resistant genes accumulation. to further examine the influence of artificial selection on the repertoire of plant rlks, new genetic resources for closely related wild plants are necessary. functional expression analysis of the lrr ii subfamily members in tomato and arabidopsis the lrrii subfamily contains rlks with dual functions in development and defense response. characterized members of this subfamily include serk genes, which are associated with diverse processes, such as brassinosteroid signaling, flagellin, cell death, light and pathogen-associated molecular pattern responses, and nik genes, which interact directly with geminivirus nsp during viral infection. phylogenetic and protein structure analyses on lrrii subfamily members of tomato and arabidopsis demonstrated that this group is highly conserved between these species. in rice, in which the rlk superfamily has undergone a large expansion, the lrrii subfamily members are also conserved in number and sequence, indicating that biochemical pathways regulated by lrrii-rlks have essential and conserved roles in angiosperm species. although lrrii members have well-conserved amino acid sequences among various species, expression analysis of the members of the tomato and arabidopsis lrrii subfamilies demonstrated that only a few of the orthologous pairs resemble in their expression profiles. by analogy with some evidence in the literature from other plant species, one may envision that these orthologous genes could have functionally diverged after the speciation event separating tomato and arabidopsis. functional divergence in orthologous genes is not an uncommon event in both plants and animals. for example, the crabs claw transcription factor in arabidopsis is expressed in the carpel primordial abaxial region and in floral nectarines and regulates carpel morphology and nectar development, whereas its orthologous in rice, drooping leaf, is expressed in the whole carpel primordium and in central undifferentiated cells of leaves, where it regulates carpel identity and midrib development. the expression of orthologous genes has also been shown to vary differently in response to a stress condition. in barley, which is tolerant to salinity, the expression of genes involved in root development, such as constans-like, is suppressed by high salinity, whereas the expression of the rice orthologous is unchanged under the same stress condition. likewise, a large fraction of orthologous pairs of rice and arabidopsis genes with receptor activity do not display conserved co-expression. therefore, different patterns of expression between the orthologous genes in the lrrii subfamilies from tomato and arabidopsis may be a result of functional divergence that occurred between these genes. functional divergence in receptor proteins with a developmental function may lead to a dramatic change in the plant phenotype because plant development is heavily guided by external signals. for example, a tissue that displays high expression of certain rlks is likely to be more sensitive to perception of rlk-specific sensing signals, leading to a rapid and effective response. in contrast, reduced expression of an orthologous gene from a different species in the same tissue would decrease the effectiveness and delay the signal perception and response. this difference in the cell responsiveness to a specific signal could represent the differential timing of biochemical reactions that are regulated by this signal. in developmental process, small differences in reaction time may be sufficient to generate a distinct phenotype in the plant. in contrast, the orthologous pairs serk1/serk2/solyc04g. which display similar expression profiles, also contain the most conserved extracellular and intracellular domains. the other two orthologous pairs, nik1/slnik and nik3/slnik which also displayed similar expression profiles, also had highly conserved extracellular and intracellular regions. in both orthologous pairs, sequence identity was approximately 65% in the extracellular regions and approximately 80% in the intracellular regions. this finding may indicate a tight conservation of function between these members of the arabidopsis and tomato lrrii rlk subfamilies. conservation of geminivirus interactions with members of the rlk family in tomato although most of lrrii subfamily orthologous pairs exhibited functional divergence, we showed that the tomato orthologs of the lrrii-rlks members nik nik and nik retain the capacity to interact with geminivirus nsp in yeast. at least for the nik and nik ortholog pairs, the functional conservation associated with specific protein-protein interactions may be linked to the high conservation of their nsp-interacting kinase domain and similarity of expression profiles. the current model of nik-mediated defense response posits that the immune receptor protects plant against geminiviruses by phosphorylating the ribosomal protein l. phosphorylation of rpl by nik redirects the ribosomal protein to the nucleus, where it may mount a defense mechanism to prevent viral proliferation. during geminivirus infection, nsp interacts with the kinase domain of niks to inhibit their kinase activity, preventing activation of the defense response. despite the high similarity between nik genes and serk genes, atbak1/serk and atserk do not functionally replace the atnik role in transducing an antiviral signaling response and do not interact with the viral nsp. in contrast, we found that the atbak ortholog from tomato interacts with nsp in yeast. although the functional relevance of this interaction in planta remains to be determined, it is worth noting that the expression profiles of the bak orthologs are not similar, as would be expected for functionally divergent orthologs. although both orthologs are ubiquitously expressed in the cognate plant species, they are expressed to different extents in distinct organs. whereas atbak expression is quantitatively similar and relatively low in all organs analyzed, its ortholog from tomato displays a higher level of expression in the cotyledons, hypocotyls and leaves, where geminivirus infection largely takes place. therefore, the expression profiles of the nsp interactors seemed to be related with the onset of geminivirus infection. due to the high expression of the atbak tomato orthologous in leaves, one may envision the existence of evolving selective pressures to diverge the corresponding nsp-interacting domains of the bak orthologs towards functional fitness with regard to geminivirus infection. in contrast to nik receptors, which are inhibited by nsp interaction, atnsak, a member of perk subfamily, interacts with nsp and phosphorylates the viral protein in vitro. loss of nsak function enhances tolerance to geminivirus infection, indicating that atnsak is a positive contributor to geminivirus infection in arabidopsis. here, we showed that the nsak tomato orthologous retains its capacity to interact with viral nsp. this demonstrates that specific members of the rlk family have conserved defense functions or compatibility functions in response to viral infection. due to the emergence of new species of tomato-infecting begomoviruses that rapidly evolve through recombination or pseudo-recombination to produce divergent genome sequences that gives the virus an advantage over its host’s recognition system, a survey of the interactions between nsps from distinct tomato-infecting geminiviruses and slniks and slnsak may add insights into the co-evolution of the viral protein and host defense/compatibility functions.  the rlk superfamily is a large and diverse group of transmembrane receptors that enables plants to perceive a diverse array of signals at the cell surface, creating an efficient mechanism for cell-environment communication. in this investigation, we generated a complete list of the members of the tomato rlk superfamily, which is made-up of proteins. the tomato rlk sequences exhibited a typical receptor-like kinase configuration and almost all of them were phylogenetically clustered with at least one member of the arabidopsis rlk superfamily. therefore, the tomato rlk superfamily is similarly organized, with the same number and identity of subfamilies as previously defined for arabidopsis rlks. among the rlk subfamilies, twelve showed specific and extensive expansion or reduction in the number of their rlk members, which may be a reflection of lineage-specific responses to various biotic and abiotic stresses. the intense breeding programs tomatoes have been subjected to may also have contributed to the establishment of the current rlk superfamily profile in this species. this comprehensive analysis comparing the complete repertory of arabidopsis and tomato rlks may provide a framework to rationalize future functional studies of the members of this family. phylogenetic and structural analyses of lrrii subfamily members from both tomato and arabidopsis reveal a well-conserved group both in terms of sequence and protein domain organization. as a consequence, the tomato lrrii-rlk subfamily is organized into the same three with phylogenetically supported clades, serk, nik and lrriic clusters. nevertheless, a comparison of the expression between orthologous genes of this subfamily demonstrated that the majority of the orthologous pairs did not share a similar expression profile, indicating that these orthologous lrrii-rlks may have undergone functional divergence. this finding is supported by the observation that, in contrast to the arabidopsis atbak slbak interacts with the geminivirus nsp and is highly expressed in leaves and the cotyledon. this pattern of slbak expression is consistent with the pattern of infection by tomato-infecting begomoviruses, which infect leaf tissues and move through the phloem but do not invade roots. additionally, as immune receptors, the orthologous pairs nik and nik displayed both the capacity to interact with the begomovirus virulence factor nsp and expression profiles that parallel the onset of begomovirus infection. evidence for functional conservation between nik orthologs has been previously provided with the demonstration that nik from arabidopsis is capable of protecting tomato plants against tomato-infecting begomovirus. collectively, our results indicate that nik orthologs retain similar functions as defense receptors to protect plant cells against viral attack. therefore, nik-mediated antiviral signaling likely also operates in tomato, suggesting that the tomato niks may be good candidate targets for engineering resistance against tomato-infecting begomoviruses. 
despite the relevance of the rlk family and the completion of the tomato genome sequencing, the tomato rlk family has not yet been characterized, and a framework for functional predictions of the members of the family is lacking. to generate a complete list of all the members of the tomato rlk family, we performed a phylogenetic analysis using the arabidopsis family as a template. the tomato rlk superfamily is made-up of proteins that form a monophyletic tree with the arabidopsis rlks and is divided into subfamilies. therefore, the tomato rlk family shares functional and structural conservation with arabidopsis. in addition to suggesting a functional analogy between these phylogenetically classified orthologs, these results expand our previous observation that nsp-nik interactions are neither virus-specific nor host-specific. nevertheless, the majority of the orthologous pairs did not display similar conservation in the gene expression profile, indicating that these orthologs may have diverged in function after speciation. a total of rlks were identified in the tomato genome, which were organized into the same subfamily clades as arabidopsis rlks. we also characterized the lrrii-rlk family by phylogeny, genomic analysis, expression profile and interaction with the virulence factor from begomoviruses, the nuclear shuttle protein. the tomato orthologs of nsp interactors, slniks and slnsak, interacted specifically with nsp in yeast and displayed an expression pattern consistent with the pattern of geminivirus infection. based on the fact that members of the arabidopsis lrrii subfamily interact with the begomovirus nuclear shuttle protein, we examined whether the tomato orthologs of nik, bak and nsak genes interact with nsp of tomato yellow spot virus. therefore, nik-mediated antiviral signaling is also likely to operate in tomato, suggesting that tomato niks may be good targets for engineering resistance against tomato-infecting begomoviruses. for the lrrii-rlk members slnik and slnik we observed functions analogous to those of their arabidopsis counterparts with respect to protein-protein interactions and similar expression profiles, which predominated in tissues that support high efficiency of begomovirus infection.
 receptor-like kinases play key roles during development and in responses to the environment. despite the relevance of the rlk family and the completion of the tomato genome sequencing, the tomato rlk family has not yet been characterized, and a framework for functional predictions of the members of the family is lacking.  to generate a complete list of all the members of the tomato rlk family, we performed a phylogenetic analysis using the arabidopsis family as a template. a total of rlks were identified in the tomato genome, which were organized into the same subfamily clades as arabidopsis rlks. only eight of rlk subfamilies exhibited specific expansion/reduction compared to their arabidopsis counterparts. we also characterized the lrrii-rlk family by phylogeny, genomic analysis, expression profile and interaction with the virulence factor from begomoviruses, the nuclear shuttle protein. the lrrii subfamily members from tomato and arabidopsis were highly conserved in both sequence and structure. nevertheless, the majority of the orthologous pairs did not display similar conservation in the gene expression profile, indicating that these orthologs may have diverged in function after speciation. based on the fact that members of the arabidopsis lrrii subfamily interact with the begomovirus nuclear shuttle protein, we examined whether the tomato orthologs of nik, bak and nsak genes interact with nsp of tomato yellow spot virus. the tomato orthologs of nsp interactors, slniks and slnsak, interacted specifically with nsp in yeast and displayed an expression pattern consistent with the pattern of geminivirus infection. in addition to suggesting a functional analogy between these phylogenetically classified orthologs, these results expand our previous observation that nsp-nik interactions are neither virus-specific nor host-specific.  the tomato rlk superfamily is made-up of proteins that form a monophyletic tree with the arabidopsis rlks and is divided into subfamilies. few subfamilies have undergone expansion/reduction, and only six proteins were lineage-specific. therefore, the tomato rlk family shares functional and structural conservation with arabidopsis. for the lrrii-rlk members slnik and slnik we observed functions analogous to those of their arabidopsis counterparts with respect to protein-protein interactions and similar expression profiles, which predominated in tissues that support high efficiency of begomovirus infection. therefore, nik-mediated antiviral signaling is also likely to operate in tomato, suggesting that tomato niks may be good targets for engineering resistance against tomato-infecting begomoviruses. receptor-like kinaseniknspbak1virus infectionprotein-protein interactionfunctional divergenceplant signalingsolanum lycopersicum 