In accordance with recent opinions, we distinguish two main layers of plant immunity, i.e. pathogenic pattern-triggered immunity (PTI) and effector-triggered immunity (ETI). These layers are frequently accompanied by activation of kinase cascades, Ca2+ influxes, generation of reactive oxygen species, transcriptional reprogramming, phytohormone signalling, proteasome degradation pathways and more specifically (i.e. in special cases of virus infections) also RNA silencing machinery  . PTI is initiated by interactions of pathogenic motifs (patterns) with cell surface signaling molecules exposed to the cell wall mostly belonging to the superfamily of receptor like kinases (RLK;    ). Among familiar groups of RLK, serine/threonine kinases receptor-like kinases (STRK) are frequent. In addition to their specific catalytic domains, these STRK can contain binding sites formed by e.g. lectin or leucine-rich repeat (LRR) domains    . ETI is mostly mediated by intracellular receptors called nucleotide binding sites/LRR proteins (NLR). Molecules of NLR are composed of i) nucleotide-binding LRR domain and ii) alternatively of either coiled-coil (CC) or Toll-like/interleukin 1 receptor (TIR) domain   .
The domains mentioned above, i.e. catalytic domains or STRK and binding TIR or LRR domains, also take part in mechanisms of metazoan immunity  . In accordance with this parallel occurrence, we can pose a question of whether at least some ancestor-like traces of typical animal superfamilies can be found in higher plants (Embryophyta), as consequences of horizontal transfer or co-evolution of superfamiliar ancestors. This concerns first of all traces of 800 million-years-old immunoglobulin (Ig) superfamily representing a very important group of immune proteins in most of metazoans including spongi and vertebrates    .
In our four preceding papers    , we hypothesized step-by-step a possible role of antibody-like phosphorylation sites (ALPS) or the corresponding nucleotide repeats in the evolution of antigen receptors and their ancestors. Consequently, ALPS sequences were used together with specifically restricted segments of conserved Ig domain sequences (Ig-cd) to compose four multiple protein sequence queries (MPSQ) inputted in our starting procedural step including BLAST searches. More precisely, the corresponding five-step selection procedures and the following data analysis were performed as described in Figure 1. Concerning superior (“non-refutable”) conserved domain similarities (cds) of Ig-cd (i.e. Ig-cds) present in the set NRI2 (cf. Figure 1), we found cds-recording files (cds-files) representing protein sequences attaining i) significant Ig-cds (p < 0.01), ii) quasi-significant Ig-cds (0.01 ≤ p < 0.1) and iii) less significant Ig-cds supported by FFAS-fold searches or literature contexts. Two types of Ig-cds were distinguished in NRI set. Dominant Ig-cds attained superior regional evaluation of cds. Recessive Ig-cds co-located with cds of non-Ig domains achieving higher bit score but sometimes an interesting context. Two significant dominant Ig-cds included bacterial Ig domains. In addition to Ig-cds evaluation, we indicated multiple associations of robust cds of catalytic-kinase-, lectin-, LRR- and TIR- related-domains, broadly occurring in the selected cds-files.
Figure 1. Five-step selection and the following analysis of Ig-similar and Ig-domain-related sequences. A concise methodological overview related to the content of our paper can be seen here. E—Expects, S—sample sizes.
2.1. Overall Description of Softwares and Procedures
2.2. Multiple Protein Sequence Queries (MPSQ)
Conserved variable Ig domains (IgV-cd) were used to form two MPSQ (MQI), i.e. MQI1 and MQI2. Segments forming our MQI composed the sequence block presented in Figure 2 of our previous paper . Besides local high sequence similarities in the selected sequence sub-blocks, additional phylogenic parameters decided about the final restriction of MQI segments. MQI1 included segments of sub-block (pm3) restricted based on fold relationships (block positions 15 - 42), the corresponding consensus segment and accompanying pattern WXXQXP. MQI2 comprised pattern DX(3)YXC and the segments containing common amino acids L,D,Y,C (block positions 81 - 97) restricting also the corresponding conserved sub-block of IgV-cd-related invertebrate sequences . Consensus was not used to form MQI2 due to its identity with the participating segment of cd00099. Two ALPS-related MPSQ (MQP), i.e. MQP1 and MQP2, were composed of the two different groups of ALPS achieving superior evaluation in database confirmation or prediction of ALPS. For the ALPS groups and protocol describing MPSQ formation from segment sequences and spacers see .
Figure 2. Quasi-ternary approaches (Q3A). In accordance with Methods, Q3A added a weak condition (quasi-third dimension) to the selection of common sequence IDs coming from certain pairs of specific BLAST mega-records obtained always with two different MPSQ. This weak condition was realized via four alternative ways, when requiring: i) additional participation of PHI-BLAST-associated sequence pattern in the selection process (this condition restricts the collection of F-sets), ii) occurrence of selected sequence ID in two mega-records of searches differing only in adjusted matrix (double matrix sets, i.e. DM-sets), iii) the presence of one or two words or abbreviations in molecular titles associated with the pre-selected sequence IDs (TI-sets) and iv) additional combinations of set fusion and conjugation in cases considered as spreading or analytical for our data processing (sets CUFC, NPoP and ter+). Ø—BLAST mega-records are formally assumed but in fact do not exist (this status followed from the absence of sequence patterns associated with MQP1 or MQP2 and impossibility to perform PHI-BLAST with the matrix PAM250); panI1, panI2—sets prepared by conjugation of five sets (cf. Section 2.6); q3dim—condition representing “quasi-third” dimension; oP—sets obtained with ALPS-related MPSQ (i.e. MQP1 or MQP2); -P, -NP, -NPoP—disjunctive subsets obtained only with, only without or simultaneously with and without ALPS-related MPSQ, respectively (similarly to the unique subset—NPoP, the set NPoP exists); r_partial—only incomplete list of r-values was used due to incompatibilities described in explanation of abbreviation Ø.
2.3. Two Different Restrictions of Sequence IDs Coming from Top Sequence Items of BLAST Records
Five parallel “sample records” differing in limiting Expect values or the download date but belonging to the same type of BLAST records were obtained in each of sixteen possible cases (cf. Figure 2). “Sample records” differed in the numbers of extracted items when forming top10 and t100 samples composed of upper 10 or 100 sequence items, respectively. These records were specifically non-redundantly fused forming S[i, r]-related sets top10[i, r] and t100[i, r] keeping sixteen lists corresponding to different types of multiplied BLAST records (cf. Figure 2). Terminal top10-NP subset (related to MQI1 or MQI2 inputs) and top10-P subset (related to MQP1 or MQP2 inputs) were prepared by selective non-redundant fusion of the corresponding top10[i, r]. Subsets of t100[i, r] were processed similarly to quasi-ternary DM-sets D[i, j, z] (described in Figure 2) using the same PHP program for the case of a less complex formula comprising two dimensional D[k, z] (cf. also WP1.1):
This means that non-redundant fusion of ALPS-unrelated (k = 1 - 10) and ALPS-related (k = 11 - 16) C[k] formed immediately disjunctive subsets t100dm-NP and t100dm-P, respectively.
2.4. BLAST-Derived Mega-Records Necessary for Quasi-Ternary and Multi-Cumulative Approaches
The aim of the searches was to obtain as many as possible long BLAST records of sequence items. Consequently, the productive Expect limits achieved the value range 2 × 105 - 2 × 107. This determined the 20000 sequence IDs for each BLASTP record, whereas PHI-BLAST searches yielded lower maximum numbers of sequence IDs, i.e. 6899-15.300. Each of the sixteen types of mega-records (cf. Figure 2) was downloaded in at least three versions differing in the date of download or limiting Expects.
2.5. Quasi-Ternary Approaches (Q3A) in Selection of Sequence IDs
Q3A represented compromising solution between selections of too extended sets determined by certain conjugations of two starting sets and too poor sets restricted by the conjugation of three sets. In accordance with Figure 2, Q3A comprised four manners how to perform weakened third dimension (quasi-third dimension) of SW architecture. For additional details see Figure 2.
2.6. Selection of Sequence IDs Based on Multi-Cumulative Approach
We selected sequence IDs simultaneously found by the same MPSQ in five required BLAST mega-records (cf. Figure 2), i.e. in three BLASTP and two PHI-BLAST searches. Due to the absence of patterns associated with MQP1 and MQP2, i.e. disability of these MPSQ to perform PHI-BLAST searches, only the lists related to MQI1 or MQI2 determined the sets of cds-files entitled panI1 and panI2, respectively.
2.7. Regularity of cds-Related Subsets
To avoid a false selection of cds-files containing only domain references but not searched cds, three keyword candidates (i.e. Pssm-ID, accession and name) were tested for numbers of selected cds-files. At least two of these numbers had to achieve minimum values to confirm regularity of the keywords associated with the minima (cf. WP2.1).
2.8. Two Set Families Were Derived When Removing Sequence Redundancy of cds-Files
Current conserved domain searches (see above) determining cds restricted by p < 0.01 and sample size 500 were performed with all items pre-selected by the preceding procedures (Figure 2). This resulted in a class of sets with starting sets of cds-files. The subsets of cds-files with positive domain occurrence of valid cds were reduced when removing sequence redundancy (see above and WP1.2) determining NS1a family of non-redundant sets of cds-files (i.e. NS1aF sets). Sequence redundancy was also checked in the cds-files without valid cds under modified conditions represented by limit E = 100 and sample size 2000. This restricted the NS2b family of cds-file subsets (NS2bF) and the corresponding non- redundant subsets of empty files of family NS1b (NS1bF). The following fusion of the corresponding sets of NS1aF and NS1bF then yielded sets of the reference family of NS1 sets (NS1F). Similarly to the formation of NS2bF and NS1F, we determined sets of NS2a family (NS2aF) including the same sequence IDs like NS1aF sets and then completed the assembly of sets of NS2 family (NS2F).
2.9. Three-Step Selection of Model Tyrosine Kinase Domains Forming Robust cds with Sequences of Plant Proteins and the Proteins of Early Metazoans Containing Variable Ig Domains in Addition
In the first step, we searched for NS1F-related files including receptor terms: “Ig) domain”, “Ig domain”, “Ig)-like”, “Ig-like”, “Immunoglobulin”, “B-cell receptor” and “T-cell receptor”. In the second step, the records of conserved domains extracted from the files were reduced to those including required cell type or cell-type-associated process, i.e. looking for the terms: “B-cell”, “T-cell”, “lymphocyte”, “lymphoma”, “amutoimmune”, “leukemia”, “macrophage”, “phagocyte”, “immune system”, “immunity”, “hematopoietic”. In the third step, we kept only such regular records of conserved domains whose names compose cds-files of three well defined signaling molecules containing both IgV-cd and tyrosine-kinase activity (RTK, SRTK, GCTK2) endogenously expressed in the immunologically important model living fossil Geodia cydonium  . Since cd05034 denotes the family of Src kinase-like protein tyrosine kinases including four kinase subfamilies important for specific immunity (Lck, Blk, Lyn and Fyn), we substituted the preceding selection procedures with the knowledge about the linkage of this domain to specific immunity and fossil IgV-cd mentioned above and complemented then the former set of nine selected cdigvtk-related accessions with cd05034 (for the list of ten selected cdigvtk see Results). Four types of strategies comprising searches with ten selected regular domain keywords (see above) of cdigvtk (dci) were employed, i.e. i) enumeration of the most frequently selected subsets of the same cds-files (mfs) when using dci, ii) identification of item number determined by any of dci (cdigvtk(max)) and iii) enumeration of cds-files containing all dci, i.e. total coincidence of dci (cdigvtk(tc)).
2.10. Statistical Evaluations
For multiple notes to statistical processing including enumeration of odds ratio values (OR and OR* for 2 × 2 tables including zero values) and t-test see our preceding paper , its important associated sources   and certain improvements mentioned in the webpage (WP2.2-4). Maximum Expect limit (E[max]) for significantly “non-refutable” Ig-domain related (NRI) cds was derived based on: i) significance limit related to inverse-phenomenon (for p < w max = 1 − w, i.e. all E achieving E ≥ E[max] can be seen as significantly refutable), ii) current BLAST formula determining relationships between Expect and p-values  and iii) the value w = 0.01 limiting significance in CDD searches:
E[max] enabled us to assemble of Ig-cds-related set selected from all sets of NSI2F. This set was called NRI2 set (cf. Figure 1). For evaluation of Expect-related specificity and overall hierarchy of levels classifying cds in our paper see WP2.5 and WP2.6, respectively.
2.11. Terms, Acronyms, Abbreviations, Texts Denoted by WP-Associated References and the Reasons for Color Grading in Map-Like Pictures or Absence of Phylogenic Trees
The corresponding information composes our web page supplement, i.e. a pdf-file accessible in the corresponding section of our web page http://www.papersatellitesjk.com or via e-mail correspondence.
3.1. Selected Records, Sets and Subsets
The starting set included sixteen types of BLAST search records containing 280.423 sequence IDs (for details see WP3.1-2; cf. Section 2.11). The selection procedures and the following removal of sequence redundancy resulted in 1323 cds-files present in 17 sets (composing the family of NS1 sets, i.e. NS1F). The same sequence IDs and sets then similarly composed the family of NS2 sets (NS2F) when using other limits for Expects (E ≤ 100) and sample sizes (s ≤ 2000; cf. Section 2.8). In accordance with the usage of ALPS-related MQP in selection of cds-files, we distinguished between ALPS-related, ALPS-unrelated, ALPS-diversified and undiversified intertwined sets (ter+, NPoP). In fact, six ALPS-diversified sets were derived here. These sets contained pairs of disjunctive ALPS-related and ALPS-unrelated subsets (denoted -P and -NP, respectively), whereas only unique additional undiversified subset was found in the case of resis set (subset resis-NPoP including unique cds-file). In addition to these primarily derived sets and subsets, 238 cds-files in different NS2F sets included conserved Ig-domain similarities (Ig-cds) in CDD records. 223 of these cds-files were non-redundant with respect to the recorded sequences. Forty-nine of the pre-selected 223 cds-files of NS2F belonged to those containing “non-refutable” Ig-domain-related cds (for the corresponding Ig-cds restriction see Section 2.10), i.e. they formed NRI2 set. The same sequence IDs like those present in cds-files of NRI2 enabled then to assemble NRI1 set from sets of NS1F origin. Both NRI sets then consisted of three subsets, i.e. NRI-NP, NRI-NP and NRI-nd, where the latest name denotes cds-files coming only from undiversified intertwined sets of NSF2. For the numbers of cds-files in all sets and subsets described here see Figure 3. For names of sets and subsets see Figure 2 and WP5.
3.2. Cds-Files, Whose Descriptions of Individual Cds Contain the Term Immunoglobulin
The robust cds with catalytic tyrosine kinases domains frequently included term “immunoglobulin” in their descriptions. This concerned also most of descriptions of cds with evolutionary important model IgV-domain-associated catalytic tyrosine kinases domains from Geodia cydonium (cdigvtk; preselected here according to Methods). The corresponding coincidence with the term immunoglobulin was demonstrated using cds of domain cd05034 (representing Src-kinase-like family) achieving the highest score among cdigvtk-derived cds (Figure 3 and Figure 4). In spite of these favorable results and the robust score values evaluating cds with cdigvtk (55 - 140 bits), the other kinases (serine/threonine kinases without immunoglobulin-related contexts) achieved considerably higher scores in their cds with the same sequence regions (Figure 4). In addition, the hierarchy of cds demonstrated in Figure 4 appeared to be consistent with the hierarchy of the cds-file- related subsets composing the receptor-kinase-domain-rich set RoK representing the set attaining the most statistically deviated occurrences of cdigvtk (cf. Table 1) and ALPS (104 cds-files; p < 0.05 in t-test):
where square brackets include numbers of cds-files (NCF) specifically double-
Figure 3. Introductory evaluation of the selected Ig-similar and Ig-domain-related sets and subsets of cds-files. The heights of columns shown here illustrate the numbers of sequence IDs and the corresponding cds-files (the files recording sequence-ID-derived conserved domain similarities, i.e. cds) in the selected sets and subsets. Set/subset-related clusters of columns are arranged descendently according to the heights of their second columns. Terms res and t100 substitute here set names resis and t100dm, respectively. For explanation of full set and subset names see Figure 2 or WP5. Sections (a) and (b): Sizes of sets and subsets. Colors: black—primarily derived sets and subsets; gray—sets and subsets obtained after a two-step procedure abolishing sequence redundancy (cf. Section 2.8, Figure 2 and WP1.2). Sections (c) and (d): Coincidences of term “immunoglobulin” and presence of cds with cd05034 in the selected sets and subsets. Only positive subset occurrences were displayed here. Extremely high and almost the same NCF (i.e. 113 for term “immunoglobulin” and 112 for other two cases found in RoK set and 99 resulting in all cases of RoK-P subset) were not comprised into the graph due to worsened resolution of complete presentation. Colors: black—searched term immunoglobulin; dark gray—cds with cd05034 (cf. Sections 2.9 and 3.2); light gray—coincidence of both indicated items in each enumerated cds-file.
restricted in cases of domain keywords (cf. Sections 2.7 and WP2.1); cdigvtk(max), cdigvtk(tc) denote subsets including cds with any of cdigvtk or all cdigvtk, respectively; cd05034, cd05102, cd14066, PLN00113 represent subsets including cds with catalytic domains of kinases, i.e. Src-kinase-like family, tyrosine kinase of VEGF receptor 3, STKc-IRAK (i.e. interleukin-1 receptor associated serine/threonine kinases; mostly dominant cds co-locating with cdigvtk) or LRR-associated kinases (achieving global maxima in Figure 4), respectively; I3R constitutes the search for two types of global keywords of plant immunity, i.e. conjugation of strategy immun3 (cf. Figure 5) and the term resistance (I3R is significantly associated with RoK as follows from specifically increased NCF and NCF-derived fractions of cds files, i.e. FCF; for both NCF and FCF holds p < 0.001; for details see WP3.4); Ig denotes a subset of cds-files containing the term immunoglobulin.
The only difference in a single cds-file can be seen in Formula (3), when comparing large RoK-derived subsets Ig with the cdigvtk(tc). The presented subset hierarchy in RoK set and data in Figure 4 moreover suggest phylogenic importance of two kinases, i.e. PLN00113 and cd14066 (see Section 4.3).
Figure 4. Statistically based comparison of co-locating cds indicating RTK- or STRK-related catalytic domains. The extents of reciprocally co-locating robust cds with catalytic kinase domains were evaluated in two preformed model subsets of cds-files limited by E ≤ 0.01 and content of the term immunoglobulin (i.e. RoK-rc and NRI-mc; for details see below). Bit score of cds with cdigvtk (i.e. catalytic tyrosine kinase domains accompanying IgV domains in molecules of Metazoa; cf. Section 2.9) were used as compared model bit scores (MBS(i)) to perform united t-test statistics concerning subset-related sample sizes 117 and 63 (see WP2.2; cf. also ). This enabled us to made two color maps of statistical relationships. As follows from these maps and cds-records, PLN00113 achieved extreme variability of bit score values in its cds (maximum SD and Rel-SD) accompanied by considerable differences in cds lengths (data not shown) and maximum bit score values in both mentioned subsets (cf. black elements). In contrast, minimum Rel-SD was found in rows of cd14066 mostly superior in individual cases of evaluated sequence IDs. Regular records Ras (cd00882), IRAK3 (pseudogene), Lyk3 and Sek (neither found among CDD sequence IDs) were not found in our cds-files. Abbreviations: b, u—bottom and upper bit score values within cdigvtk set in evaluated cds-file; BS—bit score; cdigvtk/out—domain cds excluded from the compared model set due to repeated Tie specificity; gray background—slight differences in mean-score-derived order were observed in RoK-rc and NRI-mc subsets; in interval of mean(MBS(i))—the evaluated value occurs within minimum interval containing all mean values of MBS(i), i.e. Mean BS values occurring in cdigvtk-related bottom parts of both maps; Max(max), Min(max)—global picture-related maximum or minimum of column-related maxima, respectively; Mean—mean score value in color representation; NRI-mc—multi-connective files (cf. the Section 3.3) representing independently also all cds-files containing the observed kinase domains present in NRI1; Rel-SD—relative standard deviation (Rel-SD = SD/Mean); RoK-rc—cds-files randomly selected from RoK subset of cds-files including the term immunoglobulin but not any of the five files composing NRI-mc subset (number of selected sequences corresponded to the formula n = ceil(N/10), where N = 121 denoted number of cds-files in RoK set); STK/STKc—catalytic domains specifically phosphorylating mainly serines or threonines (cf. Section 4.3); SD—standard deviation of selective mean value; slightly > max(mean(MBS(i))), slightly < min(mean(MBS(i)))—non-significantly different values slightly higher or lower than maximum or minimum mean(MBS(i)) values (see above), respectively; ↓, ↑—bottom and upper number limits following from the adequate two-tailed t-test statistics were used, respectively. For additional abbreviations see WP4 or WP5.
3.3. Monitoring of Conserved Domain Similarities and Terms in NS1F-Related Sets of cds-Files
In accordance with the results in Table 1, we introduced here five selected color
Table 1. Occurrences of cdigvtk-derived cds in focus of simple comparative approach.
Current t-test evaluation was used to indicate remote (statistically deviated) occurrences of cds-files with cds including cdigvtk, when evaluating numbers or fractions of such cds-files. Collection of PHI-BLAST derived F-sets (cf. Figure 2) represented here a model reference entity due to lucid and simple definition of its non-tendentiously selected sets. Abbreviations: all_num—number of cds-files (NCF) in an evaluated set; cdigvtk—tyrosine kinases associated with IgV domains broadly occurring in Metazoa (cf. Section 2.9); cdigvtk(max), cdigvtk(tc)—NCF containing at least one or all cdigvtk-derived cds, respectively (cf. WP2.1); mc—number of most frequently selected subsets (mfs) with the same file content, when using ten different cdigvtk (the expression 2x denotes the existence of two different maximum subsets); mc-, mc+—frequency of subsets which can be exclusively derived when either diminishing or extending the mfs, respectively; mc+/-—frequency of subsets which can be derived only when simultaneously diminishing and extending the mfs; bold in the column all_num—sufficient sample sizes. Colors in backround: gray—compared F-sets; black—significantly different values. For abbreviations see also Figure 2, WP4 or WP5.
maps, i.e. four maps of NS1F-sets achieving superior values in their upper sections and Ig-domain associated set NRI1. These maps demonstrated also certain relationships between the subsets within the observed sets (Figure 5; for color grading see Sections 2.11 and WP1.3).
Four alternative search strategies each looking for the occurrence of the same term written by two alternative entries (S2E) were described in Figure 5. The numbers of cds-files obtained with S2E antiviral2 in the RoK set (104) significantly differed from the corresponding numbers achieved in other sets of NS1F (p < 0.001 in t-test). 97 of these 104 cds-files contained cds of all cdigvtk. The records containing cdigvtk and terms of S2E antiviral2 were also present in all seven cds-files composing the RoK subset restricted with S2E antifungal2. In contrast,
Figure 5. Monitoring of cds- and term-related subsets of cds-files in four selected NS1F-related sets and NRI1. Three lengthwise sections are separated by rows with “**“ in each of five set-related pictures forming this figure. The first sections describe the occurrence (number) of Ig-related terms (terms and cds-associated domain identifiers) in cds-files of the evaluated sets. The second sections concern immunologically important terms, whereas the third sections comprises interesting or somewhere frequent terms. For our solution of problems with dissociation of conserved domain keywords/identifiers from cds records see Section 2.7 and WP2.1. Icons defining the extent of conjugations between subset pairs defined by rows and columns: ■—subset identities; #—at least 80% of cds-files are present in both subsets; Θ—one subset is fully included in the counterpart subset (cf. also all_num values). Abbreviations: all, al_num—color and numbers related to all items in the corresponding row or column, respectively; antifungal2, antibacterial2, antimicrobial2, antiviral2—all terms ahead of number 2 were scanned two times (with or without dash after term anti) to get a fusion subset; anti-all4—all four preceding sets were fused to a unique set; bac_ig—bacterial Ig-like domains whose accessions begin with big_ or BID_; cdigvtk(max), cdigvtk(tc)—the presence of any or all cds recording cdigvtk domains was required, respectively (cf. 2.9); cl27891—domain of the third LRR superfamily; death4—the subset of cds-files comprises all the following four terms: programmed cell death, apoptosis, necrosis, caspase; GST_N—glutathione S-transferase family; immun3—fusion of three subsets containing terms: immune, immunity, immunol; IG_FLMN, ig_v—conventional abbreviations denoting filamin-type Ig domains or multiple different variable Ig domains, respectively; LRR6x—only six leucine-rich repeat domains (LRR) formed cds-records in all our sets yielding this fusion subset, i.e. COG4886, pfam13855, cl27891, pfam12799, pfam14580 and pfam07725 (cf. WP1.4 and WP3.3); lrr* = PLN00113—the term “leucine-rich repeat” in fact restricts the same cds-files like cds with PLN00113 but much less frequently cds with LRR domains; pep U pro—a set with the terms peptidase or protease; pfam01453—domain of mannose specific lectins; TIR domain—toll-like/interleukin-1 receptor domain. For additional abbreviations see WP4 and WP5.
a unique set-derived subset of significantly increased size (NCF = 17; p < 0.01 in t-test) was selected with S2E antifungal2, i.e. a subset of F2oP1, where sixteen cds-files included robust cds with the domain of antifungal thaumatin-like proteins. This subset comprised a unique cds-file independently containing all cdigvtk-related accessions and terms of S2E antiviral2 but not the terms of S2E antimicrobial2. The cds-file of the same content like this file composed also the RoK set. The terms of S2E antimicrobial2 were not observed in RoK, but it still exhibited significantly increased occurrence in the set t100dm set (NCF = 18; p < 0.01 in t-test). The corresponding subset of exclusive term occurrence, i.e. t100dm-NP, did not include cds-files with cdigvtk-related accessions and the other compared terms of S2E antiviral2 or antifungal2. However, only six of the selected cds-files realized in fact the aim of this S2E indicating robust cds with ABC-type antimicrobial peptide transport system.
Twenty cds-files with TIR domains (item tir dom in Figure 5) were found in the resis set composed of cds-records generated by sequence IDs which contained the term resistance in the molecular title (cf. Figure 2). Other sets significantly differed from such extreme TIR domain occurrence in the number and density values (at least p < 0.001 in t-test; cf. Figure 5; the corresponding NCF achieved only values in interval 0 - 3 in other sets). Significantly increased values of densities concerned occurrences of LRR6x- and NB-ARC-related cds-files in TIR domain-derived subset of resis set (both 18 of 19 cds-files) and NB-ARC domains (35 of 38 cds-files) in the subset restricted with LRR6x (including all found LRR domains; Table 2). Simultaneous occurrence of NB-ARC, LRR6x and TIR domains in cds-files indicated the presence of 18 NLR receptors in the resis set. Additional distinct six sequences (cds-files) alternatively comprised NB-ARC, LRR6x and coiled-coil domains in accordance with alternative NLR composition described in Introduction.
Most of cds-files containing the term “-lectin” comprised the domain with accession pfam01453 (i.e. mannose specific lectin). Term lrr* (leucine-rich repeats) revealed mostly kinase domains associated with LRR but much less frequently LRR domains (for details concerning distribution of LRR domains see WP3.3 and Figure 5).
Seven multi-connective cds-files (five of RoK origin) were independently found in NRI1-set when: i) individually searching for the terms antiviral, killer, lymphocyte, macrophage, or the term programmed cell death; ii) combining the term resistance and strategy immun3 (cf. I3R mentioned above) or using conjugation strategy death4; and iii) looking for the occurrence of most frequent cdigvtk cd05034 or co-locating cds with cd14066 (Figure 5). The cds-records with TIR domains and aaa+ (superfamily of ATP-ases) were found neither in NRI1 nor in NRI2 set (Figure 5 and WP2.1). Apparent contradictions between NRI-related data present in Figure 5 and last table in fact followed only from the difference in sample sizes restricting content of cds-files composing NRI1 and NRI2 (see Figure 1).
Table 2. Statistical analysis of domain-(cds-) and term-related associations using 2 × 2 tables.
aSets or pairs of complementary subsets were compared with respect to the differences in occurrence of a certain query terms (see Q below). _—conjunction; _n—conjunction with a complementary subset; AQ+, BQ+—subsets of sets/subsets A or B positive for presence of Q or associated cds (if Q was a domain identifier), respectively; cd03053—Phi subfamily of glutathione S-transferases; cd04982—IgV domain of the gamma chain of gamma/delta TCR; NB-ARC—domain pfam00931 forming NLR proteins (cf. Introduction); Q—terms used as search queries, i.e. as: i) regular conserved-domain-derived keywords associated with cds or domain-related terms (cf. Sections 2.7, WP1.5 and WP2.1) or ii) functionally related terms; resub—residual complementary subset of compared set. Specifically shortened versions of abbreviations otherwise used in this paper: antifu2—antifungal/anti-fungal; antivir2—antiviral/anti-viral; C(m), C(t)—cdigvtk(max), cdigvtk(tc) described in Table 1, respectively; C*—case, if C(m) = C(t). For additional abbreviations see Figure 2 or Figure 5 and WP4 or WP5. bOdds ratio values were enumerated conventionally or as zero-associated Bayesian OR* (cf. Section 2.10). Statistical evaluation of the corresponding 2 × 2 tables was performed by Fisher’s exact test. Two tailed exact Fisher’s test can be always assumed as valid, whereas one-tailed test consists in the usage of additional orienting conditions (sometimes depends on interpretation or context). #—for the details see the Section 3.6; aE-n = a × 10−n; One-t, Two-t—values following from one-tailed and two-tailed Fisher’s test, respectively. cSemi-quantitative odds-ratio-related evaluation represents lucid insight into the presented data. Different comments to one and two tailed test are separated by “/”, if such difference exists. Classification of significance levels: poor—insufficiently low or none validity (p ≥ 0.1); q—quasi-significant (in this case 0.05 ≤ p < 0.1; cf. WP2.6); si—minimally significant (0.005 ≤ p < 0.05); si2—of improved significance (10−6 ≤ p < 0.005; cf. PSI BLAST limit); ssl—of superior significance (p < 10−6). Classification of linkages according to OR values: C—causal linkage (OR ≥ 10); inv—inverse linkage (OR < 1/1.4); QC—quasi-causal linkage (4.5 ≤ OR < 10; interval added here by authors to the current classification; cf. WP2.3); S—strong linkage (2 ≤ OR < 4.5); W = weak (1.4 ≤ OR < 2). dThe values ξ and η characterizing model mean set were enumerated here with the help of geometrical mean value (GM) of ratios derived from numbers of positivities and negativities typical for each topically compared/comparable set, respectively (cf. WP2.4; ξ = 65.58424789; η = 25.91575211).
3.4. Statistical Analysis of Linkages between Sets and Subsets of NS1F Based on Evaluation of 2 × 2 Tables
Three types of strong significant statistical associations following from robust cds can be seen in the different segments of Table 2: i) association of cdigvtk with RoK set or their related linkage to dominant STRK domain 14066 (segments i) and v)), ii) total disjunction of lectin (pfam01453) and LRR domains in the RoK set (segment iii)); for similar situation in the other sets see Figure 5), and iii) an extreme difference between the absence of cds with cd27891 and TIR domains in the RoK set (contrasting in case of LRR domain cl27891 with the presence of forty-three cds formed by distinct LRR in RoK) and their abundant co-occurrences in individual cds-files of resis set (segment v) and Figure 5). Though the linkages between cdigvtk and ALPS-related segments were mostly strong, they achieved only boundary line of significance (segment ii)).
3.5. Statistics of ALPS Occurrence
The occurrence of ALPS-selected cds-files (ASC; ASC represents molecules or sequence IDs selected with the assistance of sequences of ALPS forming two sequence queries MQP; cf. Section 2.2) was evaluated with respect to the number of selected cds-files (four sets with oP in the name and six subsets with -P in the name) and their relative occurrence (six pairs of subsets). Two values were significantly higher in t-test than the values in the compared collection of the residual sets (cf. WP2.2), i.e. i) maximum ratio (13) between ALPS-positive and such negative samples in CUFC set (p < 0.01) and ii) the number (104) of ASC in RoK set (p < 0.05). Significantly deviated maximum odds ratio (OR) and the best significance level of OR evaluation, selected CUFC and RoK sets as significantly ALPS-associated sets, respectively. For the data and methodology concerning model odds ratios see the sixth section of Table 2 and WP2.4, respectively.
3.6. Conserved Ig Domain Similarities in the Maximum Expect Extension of NSF2 and NRI2 Sets
Figure 6 represented not only a map of similarities in broad range of Expects but also starting point for screening of the corresponding Ig-cds (see also WP3.5). Among others the domain cd04982 (IgV_TCR_gamma) achieved significantly increased NCF (p < 0.05) in the t100dm-NP subset of NS2F (NCF = 18 holds also for the set t100dm). Twelve of the corresponding eighteen cd04982-related sequence IDs restricted cds-files specifically including cds with cd03053 (glutathione S-transferase GST_N_Phi frequent in Arabidopsis and Oryza genomes) in NS1F-related t100dm-NP set. This indicated a significant molecular association of short TCR-gamma-related peptides and cds with cd03053 (Table 2). These immunologically interesting peptides were found among others in certain species of vegetable and tobacco (for further comments see WP3.6).
3.7. More Detailed Description of the Selected “Non-Refutable” Ig-cds Forming NRI2 Set
Only seven sequences included “non-refutable” Ig-cds, which were called here as dominant Ig-cds, i.e. as Ig-cds that achieved the highest score and minimum Expect among cds co-locating with the same segments of individually evaluated sequences (cf. Table 3). Three and four of the selected sequences contained bacterial and metazoan dominant Ig-cds, respectively. Two significant and one quasi-significant dominant Ig-cds including bacterial Ig domains occurred in different sequence regions of nuclear pore complex protein GP21 (sequence ID XP_010248630.1; Table 3). All these three segments restricted by Ig-cds were approved by our Ig-fold evaluation. Another protein LOC105056499 (XP_010937019.2) was unique one selected by three different procedures. Its Ig-cds was composed
Figure 6. Occurrence of Ig-cds, Ig-cds-overlapping cds and Ig-domain-related terms in cds-files forming “expanded” NS2F and NRI2 sets. The cds limited by Expect values E ≤ 100, i.e. cds at least comparable with dense similarities between segments achieving lengths of secondary structures (cf. WP2.6), are described here. For distribution of re-selected “non-refutable” Ig-cds see Table 3. Abbreviations: big_/BID_—terms denoting major part of bacterial domains; cl26464—atrophin-1 family; igv_—cluster of IgV-related terms; NS2F-Ig, NRI2-Ig—content of Ig-cds and related terms or term clusters within non-redundant set NS2F-Ig (comprising 223 cds-files extracted from all sets of NS2F as subset of cds-files including Ig-cds records) and NRI2 described in Section 3.1, respectively; P1453—pfam01453 (domains of mannose specific lectins); P5938—pfam05938, i.e. domain of plant self-incompatibility protein S1; sma406—smart00406, i.e. IgV domain frequent in non-chordate Metazoa. For additional information and abbreviations see Figure 5 or Sections 2.11, WP1.5, and WP4.
of typical metazoan and dominant Ig domain Ig-cd Ig1_ IL1R_like. The segment restricted by Ig-cds achieved maximum fold score 6.81 in case of Ig light chain, but not prevailing number of required Ig folds (rule Q1 in Table 3). Nevertheless, the Ig-cds-related specificity of 96.3% was enumerated based on Expect values of all cds records overlapping the sequence segment in the cds-file of NRI2 origin (limited by E ≤ 100; for details see WP2.5). Other three dominant typical metazoan Ig-cds with sequences XP_021903093.1, XP_ 017226294.1, XP_022544643.1 achieved higher Expects than XP_010937019.2, but were fully approved by our folding analysis with FFAS03 (Table 3). Lower fraction (23.8%) of FFAS03-approved sequences was selected among forty-two “non-refutable” Ig-cds called here as recessive Ig-cds and representing Ig-cds co-locating with more valid non-Ig cds. Three types of interesting phenomena were observed in the case of recessive “non-refutable” Ig-cds: i) FFAS- and CDD-confirmed chimera of filamin domains and Ig-cds (see Table 3 and Section 4.1), ii) the existence of a plant self-incompatibility protein looking like a functional and structural analogue of plant T-cell receptor (E = 5.0 for co-locating weak cds with IgV-TCR_gamma) and iii) extremely frequent co-locations of weak Ig-cds (4.605 ≤ E ≤ 100) with “non-refutable” Ig-cds. The last phenomenon concerned mainly thirteen and fourteen weak Ig-cds co-locating with five and two also co-locating “non-refutable” Ig-cds found in cds-records of sequences XP_021662681.1 (ALE2 STRK) and XP_009381707.1 (potassium transporter) by CDD searches, respectively (cf. Table 3).
Table 3. Three groups of re-selected plant protein sequences forming “non-refutable” Ig-cds.
aSequences and their cds files origin. Search strategies are identified here by names of the selected sets (for procedures see Figure 2). Selected sequences are arranged ascendantly according to minimum or unique Expects concerning their Ig-cds. For additional comments to the selected molecules see Results and Discussion. Abbreviations in names of proteins: i.—isoform number; p.—protein; RK—receptor kinase; rest obs—restored obsolete file; STRK— serine/threonine receptor-like protein kinase; u.p.—uncharacterized protein. bThe displayed conserved domain similarities (cds) determined by CDD searches. Framed pairs or triplets of intervals—co-locating recessive Ig-cds and dominant non-Ig cds (cf. Section 3.7); big, BID—cds of Ig-cd found in bacteria; _CAML—_CAM-like; E-n in Expect values—the order related values 10−n; 2ig, 3ig—sequence segments with Ig domain-related cds achieving local chain maxima of the second or third site-related order, respectively; dii—non-Ig, non-lectin and non-max domain including less similar (recessive) Ig-cds; Ig2_LILR_KL—Ig2_LILR_KIR_like; IG_FLMN—filamin-type Ig-cd; iit, sit—records of associated domains with robust mostly inferior or superior cdigvtk similarities, i.e. cd05102 and cd05034, respectively; max, mig, mle—maximum cds, Ig-cds or cds of lectin domain found in cds-file of the referred sequence, respectively; SF—super family/superfamily. For more detailed information about conserved domains see special option in menu on NCBI web page. For additional abbreviations see WP4 or WP5. cParameters of existing co-locating “non-refutable” Ig-related cds. CA—numbers of all co-locating “non-refutable” or solely found Ig-cds; CA positions—sequence positions of the minimum segments where all co-locating “non-refutable” Ig-cds occur. dFFAS evaluation. Symmetrical sequence extension or fusion of substantially overlapping segments located in very close positions were specifically used in the cases of too short Ig-cds-derived sequence segments (length < 25 amino acids) for FFAS evaluation. Except for the rule Q3, other five selection rules (see below) included common requirement for Ig-related fold FFAS03 scores S_ig, i.e. occurrence of S_ig ≥ 5 in at least two separate database outputs. FM—FFAS derived maximum S_ig; FR—FFAS-associated selection rules, i.e. rules R necessary for sequence occurrences in the first section of the table and twilight-zone-related quasi-rules Q; Q1—7.5 > S_Ig ≥ 5.0 and S_ig achieves maximum value; Q2—10.0 > S_ig ≥ 5.0, S_ig achieves the second or the third superior values and prevailing number of S_ig is observed among all scores ≥ 5; Q3—like in R1, R2 or R3 but S_ig ≥ 5 occur only in one database record; R1—S_ig ≥ 10; R2—10 > S_ig ≥ 7.5 and S_ig achieves maximum value; R3—7.5 > S_ig ≥ 5, S_ig achieves maximum and prevailing number of S_ig is observed among all scores ≥ 5. eMaximal HMMER-derived similarities between the displayed segments restricted by Ig-cds and non-plant sequences. HMMER evaluation achieved lower Expect values than the preceding BLAST and CDD. This let us to observe the distribution of minimum Expect values (E) enumerated by HMMER (characterizing similarities of each “non-refutable” Ig-cds-related segment selected by CDD) by means of fifteen differently restricted histograms and to estimate thus intervals important for our evaluation (see WP3.9 and WP3.10). This effort yielded two intervals indicating grading positive evaluation, i.e. i) 10−21 > E ≥ 10−26 and ii) E < 10−26, gradually restricting suspicious and hot candidate segments for recent horizontal transfer, respectively. HME—HMMER-derived Expects; HMG—generum origin of the sequence achieving maximum HMMER-derived similarity; HMT—broader taxonomic group including HMG, bold numbers—hot candidate segments for recent horizontal transfer. Abbreviated expression in the column HMG: C. Magneobact.—Candidatus Magnetobacterium. fFFAS-mediated selection. Three selection rules (R1, R2, R3) described in footnote d were alternatively required to display the corresponding data in this table section. gIg-cds added by authors. These items forming also NRI sets appear to be interesting with respect to their literature or structural contexts (cf. Sections 4.1 and 4.3).
In addition to FFAS and CDD studies, the searches for maximum sequence similarities of non-plant sequences to the Ig-cds-derived plant segments dealt with here were performed using HMMER. Six non-plant sequences, including the two ones of fungal origin, were classified here as hot candidates for recent horizontal transfer exchanges (Table 3).
4.1. Selected Ig-Domain-Related Sequence Segments
In fact, at least three main groups of Ig-cds can be considered with respect to out data as ensues from the three following paragraphs.
The presence of two dominant cds recording significant bacterial Ig domains was indicated in Results and Table 3 in case of protein GP210 (XP_010248630.1;  ). Bacterial Ig domains perhaps evolved from eukaryotic Ig domains . In accordance with Table 3, bacterial Ig-cds selected here were not classified as positive with respect to recent horizontal transfer though such earlier transfer probably occurred via bacteria as vectors. These bacterial vectors are usually delivered to plants by means of insects and Nematoda worms  .
According to our data, the sequence XP_010937019.2 encoding protein found in oil palm (Elaeis guineensis;  ) formed Ig-cds with domain Ig1_IL1R_like. This Ig-cds represented maximum dominant Ig-cds with metazoan domains in Table 3. The sequence could be interesting from the point of view of vanished structures resulted via long-lasting or intensive mutation changes. Hence all cds of the sequence achieved only very low scores. Since Ig domains or Ig folds are usually resistant to mutations, the indication of superior cds with Ig-domain within the sequence appeared to be consistent with both Ig-domain relationship and vanished structure classification. Another interesting context of the discussed superior Ig-cds consists in its functional relationship to interleukin-1 receptor (IL1R). Hence two additional and even robust IL1R-related cds were also observed in the selected cds-files, i.e. cds with Toll/IL1R (TIR) and cd14066 (IL1R associated kinase) domains. Do these IL1R-related cds represent traces of very old ancestor IL1R-mediated regulation or only different phylogenic integrations of effective structures?
The occurrence of recessive Ig-cds (i.e. Ig-cds co-locating within more valid non-Ig domains) may have four reasons: common ancestor origin, transfer of transposons and transposon-like elements, other recombination events (cf. ), or a long-lasting convergent competition involving recognition of the same ligand/ligands or related set of ligands. For instance, cds of the multi-interactive filamin domain includes here significant recessive Ig-cds recoding Ig-FLMN co-locating with “non-refutable” bacterial Ig-cds with BID_1 restricting the segments of Ig-fold-related score of maximal value in Table 3 (see XP_019237668.1 in Table 3; for associated sequencing project see    ). As known, filamin domain is composed of Ig-like beta-sandwich fold . Consequently, the presented facts could suggest a common molecular ancestor origin of filamin and Ig domains. Filamin domains are among others present in protozoans . In plants, they even contribute to evolution of cell surface proteins in lineages from a common ancestor to glaucophytes, rhodophytes and viridiplants . Since mannose-specific lectin domains prevail among the selected lectin domains, the possible structural relationship between Ig-like beta-sandwich fold and beta-Prism I fold found in plant mannose- and galactose-specific lectins  could explain the occurrence of Ig-like structures within lectin domains. The consistent presence of Ig domains in channels (cf. Table 3) has as yet been described only in mammals  . Common molecular origin of Ig superfamily and self-incompatibility protein (SIP) recorded in Table 3 represents an important and interesting clue for authors working with these proteins because structural models were developed  . Similarly to compared TCR, SIP is involved in cell death events . The inclusions of Ig-cds within several catalytic or transporter domains described also in Table 3 appear to be interesting from the point of view of conformational flexibility of Ig domains (cf. ), though the relationship to Ig domains remains unclear. The segment restricted by cds with PLN00113 in LRR STRK on the other hand surprisingly forms cds with cdigvtk domains and Ig-cds located in its different subsegments (Table 3). For comments to important significant Ig-cds of ALE2 see next section.
In contrast to high occurrence of cds-files recording cds with cdigvtk and STRK-related catalytic domains in RoK set, only low file fraction of such cds forms the NRI sets including top Ig-cds (Table 3). Nevertheless relatively increased density but low number of cds with cdigvtk was observed in cds-files of ALPS-related fraction of NRI (Table 1; cf. also Sections 4.2 and 5). This indicates considerable independence between occurrences of top Ig-cds and cds of the discussed kinase domains in sequences of higher plant protein origin. For comments to incomplete Ig-like domain folds composing also “non-refutable” Ig-cds-derived segments desribed here see the paper of Berisio  (cf. also WP3.7).
4.2. Sequence Similarities of Antibody-Like Phosphorylation Sites
As is well known, the major part of phosphorylation sites with specificity corresponding to ALPS exists in plants    . Significantly increased number and relative occurrence value concerned ALPS-selected cds-files (ASC; cf. Results) forming in RLK-related RoK set. Another, functionally related resis set (including mostly NLR proteins) then contained 71.4% of ASC whereas NRI set only 38.8%. However, it is a question whether indeed functional plant phosphorylation sites form sequence similarities with mouse and human ALPS. Though eukaryotic phosphosites are generally more conserved than their non-phosphorylated counterparts with similar structural constrains , low sequence differences can be even observed between mouse and humans orthologues . A possible alternative mechanism maintaining at least some ALPS-related structures consists in independent participation of some structural segments of ALPS in recognition (antibody-like property of ALPS). This would be in accordance with: i) assumed related effects of similarly long LRR    and ii) prevalence of ALPS-derived sequence IDs to molecules involved in plant immunity, i.e. mainly in pathogen recognition. The latter possibility would follow from:: i) significantly increased occurrence of ALPS-related cds-files in RoK set (containing PTI-related STRK involved in pathogen recognition) accompanied moreover by significant association of these files with the term antiviral (Results and Table 2) and ii) considerable occurrence of ALPS-derived sequences in the resis set containing otherwise distinct ETI-related sequences of NLR-proteins involved also in plant resistance.
Robust cds between ALPS-selected plant sequence of ALE2 (described in ) and domain of Atrophin-1 superfamily (cl26464) includes significant recessive Ig-cds with domain cd05722 (i.e. Ig1_Neogenin; Table 3) supported by a maximum number of co-locating “non-refutable” and weak Ig-cds among the cds-files of NRI2 (cf. Section 3.7). In accordance with BLASTP comparison, the corresponding Ig-cds-related segment of ALE2 contains moreover four traces and overlaps one trace of ALPS-related segments (cf. WP3.8). Two segments even look like as two repeats occurring at amino acid positions 305 - 328 and 329 - 359 of ALE2 (cf. WP3.8). Possible role of ALPS in ageing was previously discussed . Similarly members of Atrophin-1 superfamily (domain co-located here) participate in a progressive neurodegenerative disorder in vertebrates . This possible agreement raises the question whether the described coincidence of the three types of sequence similarities mentioned here can be interesting for the phylogeny of ageing.
4.3. RoK Set as a Set of Phylogenic Interest
The ALPS-associated RoK set is composed of molecules selected by combined procedures and reselected according to the content of the terms receptor and kinase in molecular title (cf. Figure 2). In accordance with this selection strategy and Results, this set contained mainly serine/threonine receptor/receptor-like kinases (STRK). Many plant STRK exhibit additional tyrosine specificity forming a group of STY (i.e. Ser/Thr/Tyr) kinases and also the corresponding sequence chimerism . We can even think about deep evolution of STY kinases leading to the last universal common ancestor (LUCA) of Archaea, Bacteria and Eukarya  . In accordance with this opinion, STY kinases achieve structural relationship to animal non-receptor tyrosine kinases, Src, Abl, Lyn, Fes, Sek, Kin and Ras as well as receptor-like kinase Lyk3 suggesting thus common superfamily origin of the compared kinases  . In comparison with these data, our subsets of cds-files specifically compared in Figure 4 did not contain the catalytic domains of Ras and Lyk3 kinases.
Occurrence of global score maxima of PLN00113 (Figure 4), superiority of its subset in RoK set (Formula (4)), its lower cds densities and lower frequencies of individual superior cds than in case of cd14066 (Figure 4) appear to be in agreement with a broader and more distant, i.e. ancient ancestor-like structural context of PLN00113 (LRR-associated catalytic kinase domain) common with respect to the evaluated kinase domains. On the other hand cd14066 (interleukin-1 receptor associated catalytic serine/threonine kinase domain STRK; i.e. STKc-IRAK), achieves the densest cds and most frequently superior kinase-related scores in individual cds-files. These scores represent sometimes even score maxima among all cds in evaluated cds-file and determine extremely low relative standard deviations (Figure 4). This extreme data can indicate closer but later structural relationship of cd14066 to ancestor structure than in case of domain PLN00113. The described cds-derived domain chimerism appears to be interesting in two aspects, i.e. i) possible structural divergence between protozoan or early metazoan RTK (including also IgV-associated cdigvtk evaluated here) and plant STY kinases     and ii) inclusion of cds recording NRI-restricted domain Ig3_L1-CAM_like (28.94 bits) and cdigvtk (co-locating with cd14066) in different sites of longer and extremely robust cds (497.45 bits) of domain PLN00113 (Table 3).
4.4. Structures in the Light of Statistics
Removal of sequence redundancy diminished the number of sequence IDs from 2217 to 1323. In fact, this processing reflected the presence of sequences encoding duplicated isoforms and alternative splicing products. In agreement with the monitored associations between variously selected sets and subsets, many significant domain linkages and disjunctions also exist (cf. Table 2 and the corresponding section of Results) as possible consequences of diversified selective processes of domain shuffling (cf. ). In accordance with statistics in Table 2 and the relationships described in Results, we can among others pose question, whether the third superfamily of LRR domains represented by cl27891 is specifically involved in recognition by ETI-related NLR composed also of TIR domains (cf. ). Several significantly increased occurrences of terms in the sets were partially explained here. This comprised the occurrences of: i) term immunoglobulin, ii) terms linked in strategy I3R (related to plant immunity) in RoK set and iii) the terms antifungal2 in the set F2oP1 (see Section 3.2). Statistically predicted hot- candidate partners for recent horizontal transfer among genes of fungal origin displayed in Table 3 represent meanwhile a question for further analysis and critical comments exceeding extent of this paper.
The described procedures selected many non-Ig molecules important for plant immunity. To explain this relationship for cases of Ig-similar proteins and proteins with weak Ig-cds, we considered mainly similarities following from motif-motif, motif-secondary structure or ligand-secondary structure interactions often generated via convergent mutation changes, repeat effects and recombination changes on DNA level      .
Due to phylogenic distance of higher plants from animals including most of Ig-domains, we enlarged the limiting Expect values to those approximating i) significant refusing of domain similarity (Section 2.10) and among others also ii) dense similarities of short fold lengths (cf. WP2.6). Since uncertainty of such weakened restriction increased, complementary evaluation of cds specificity (uniquely proved here in case of quasi-significant Ig-cds of XP_010937019.2; cf. Sections 3.7 and WP2.5) or usage of alternative online accessible methods (cf. FFAS- fold-derived analysis in Table 3) became to be important. The function of the selected Ig-cds-related segments displayed in Table 3 remains unclear. We can only remind the discussed possible relationships to interleukin-1-receptor-related structures and structures involved in ageing or self-recognition interesting for further structural investigation.
Several extremely statistically strong and significant domain linkages or disjunctions in multi-domain molecules follow from our data. Nevertheless, it remains to determine whether all these linkages are independent of our Ig-related selection machinery. This could be decided in future revision of the described associations on larger sets keeping only domain relationships. Similarly, more detailed phylogenic evaluation of relationships between cds including cd14066, PLN00113 and cdigvtk segments could be important for the studies concerning deep evolution of SRTK, RTK and perhaps also Ig-domains (see Figure 4 and below).
In accordance with data presented in Table 2 and Table 3, an alternative role of ALPS-related segments in recognition was considered in Section 4.2. This raises the question, whether phylogenic (inter-genic or intragenic) interactions between at least certain genes encoding i) Ig domain ancestors or Ig domains and ii) STRK could include transfer of ALPS-related structures (or the corresponding repeats; cf. Section 4.2) from STRK improving sometimes recognition mediated pre-Ig or Ig domains.
 Brustolini, O.J., Silva, J.C., Sakamoto, T. and Fontes, E.P. (2017) Bioinformatics Analysis of the Receptor-Like Kinase (RLK) Superfamily. Methods in Molecular Biology, 1578, 123-132.
 Dezhsetan, S. (2017) Genome Scanning for Identification and Mapping of Receptor-Like Kinase (RLK) Gene Superfamily in Solanum tuberosum. Physiology and Molecular Biology of Plants: An International Journal of Functional Plant Biology, 23, 755-765.
 Shchennikova, A.V., Kochieva, E.Z., Beletsky, A.V., Filyushin, M.A., Shulga, O.A., Ravin, N.V. and Skryabin, K.G. (2017) Identification and Expression Analysis of Receptor-Like Kinase Gene ERECTA in Mycoheterotrophic Plant Monotropa hypopitys. Molecular Biology (Moscow), 51, 780-786.
 Bacete, L., Mélida, H., Miedes, E. and Molina, A. (2018) Plant Cell Wall-Mediated Immunity: Cell Wall Changes Trigger Disease Resistance Responses. The Plant Journal, 93, 614-636.
 Hervé, C., Dabos, P., Galaud, J.P., Rougé, P. and Lescure, B. (1996) Characterization of an Arabidopsis thaliana Gene That Defines a New Class of Putative Plant Receptor Kinases with an Extracellular Lectin-Like Domain. Journal of Molecular Biology, 258, 778-788.
 Afzal, A.J., Wood, A.J. and Lightfoot, D.A. (2008) Plant Receptor-Like Serine Threonine Kinases: Roles in Signaling and Plant Defense. Molecular Plant-Microbe Interactions, 21, 507-517.
 Liu, P.L., Du, L., Huang, Y., Gao, S.M. and Yu, M. (2017) Origin and Diversification of Leucine-Rich Repeat Receptor-Like Protein Kinase (LRR-RLK) Genes in Plants. BMC Evolutionary Biology, 17, 47.
 Ma, N., Liu, C., Li, H., Wang, J., Zhang, B., Lin, J. and Chang, Y. (2018) Genome-Wide Identification of Lectin Receptor Kinases in Pear: Functional Characterization of the L-Type LecRLK Gene PbLRK138. Gene, 661, 11-21.
 Di Gaspero, G. and Cipriani, G. (2003) Nucleotide Binding Site/Leucine-Rich Repeats, Pto-Like and Receptor-Like Kinases Related to Disease Resistance in Grapevine. Molecular Genetics and Genomics, 269, 612-623.
 Bonardi, V., Cherkis, K., Nishimura, M.T. and Dangl, J.L. (2012) A New Eye on NLR Proteins: Focused on Clarity or Diffused by Complexity? Current Opinion in Immunology, 24, 41-50.
 Burdett, H., Bentham, A.R., Williams, S.J., Dodds, P.N., Anderson, P.A., Banfield, M.J. and Kobe, B. (2019) The Plant “Resistosome”: Structural Insights into Immune Signaling. Cell Host and Microbe, 26, 193-201.
 Johnson, G.B., Brunn, G.J., Tang, A.H. and Platt, J.L. (2003) Evolutionary Clues to the Functions of the Toll-Like Family as Surveillance Receptors. Trends in Immunology, 24, 19-24.
 Schäcke, H., Rinkevich, B., Gamulin, V., Müller, I.M. and Müller, W.E. (1994) Immunoglobulin-Like Domain Is Present in the Extracellular Part of the Receptor Tyrosine Kinase from the Marine Sponge Geodia cydonium. Journal of Molecular Recognition, 7, 273-276.
 Blumbach, B., Diehl-Seifert, B., Seack, J., Steffen, R., Müller, I.M. and Müller, W.E. (1999) Cloning and Expression of New Receptors Belonging to the Immunoglobulin Superfamily from the Marine Sponge Geodia cydonium. Immunogenetics, 49, 751-763.
 Müller, W.E.G. (2001) Review: How Was Metazoan Threshold Crossed? The Hypothetical Urmetazoa. Comparative Biochemistry and Physiology Part A, 129, 433-460.
 Grassot, J., Gouy, M., Perrière, G. and Mouchiroud, G. (2006) Origin and Molecular Evolution of Receptor Tyrosine Kinases with Immunoglobulin-Like Domains. Molecular Biology and Evolution, 23, 1232-124.
 Kubrycht, J., Borecky, J. and Sigler, K. (2002) Sequence Similarities of Protein Kinase Peptide Substrates and Inhibitors: Comparison of Their Primary Structures with Immunoglobulin Repeats. Folia Microbiologica, 47, 319-358.
 Kubrycht, J., Sigler, K., Růžička, M., Souček, P., Borecky, J. and Ježek, P. (2006) Ancient Phylogenetic Beginnings of Immunoglobulin Hypermutation. Journal of Molecular Evolution, 63, 691-706.
 Kubrycht, J., Sigler, K., Souček, P. and Hudeček, J. (2016) Antibody-Like Phosphorylation Sites in Focus of Statistically Based Bilingual Approach. Computational Molecular Bioscience, 6, 1-22.
 Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W. and Lipman, D.J. (1997) Gapped BLAST and PSI-BLAST: A New Generation of Protein Database Search Programs. Nucleic Acids Research, 25, 3389-3402.
 Zhang, Z., Schäffer, A.A., Miller, W., Madden, T.L., Lipman, D.J., Koonin, E.V. and Altschul, S.F. (1998) Protein Sequence Similarity Searches Using Patterns as Seeds. Nucleic Acids Research, 26, 3986-3990.
 Marchler-Bauer, A., Panchenko, A.R., Shoemaker, B.A., Thiessen, P.A., Geer, L.Y. and Bryant, S.H. (2002) CDD: A Database of Conserved Domain Alignments with Links to Domain Three-Dimensional Structure. Nucleic Acids Research, 30, 281-283.
 Marchler-Bauer, A., Anderson, J.B., DeWeese-Scott, C., Fedorova, N.D., Geer, L.Y., He, S., Hurwitz, D.I., Jackson, J.D., Jacobs, A.R., Lanczycki, C.J., et al. (2003) CDD: A Curated Entrez Database of Conserved domain Alignments. Nucleic Acids Research, 31, 383-387.
 Marchler-Bauer, A., Lu, S., Anderson, J.B., Chitsaz, F., Derbyshire, M.K., DeWeese-Scott, C., Fong, J.H., Geer, L.Y., Geer, R.C., Gonzales, N.R., et al. (2011) CDD: A Conserved Domain Database for the Functional Annotation of Proteins. Nucleic Acids Research, 39, D225-D229.
 Marchler-Bauer, A., Derbyshire, M.K., Gonzales, N.R., Lu, S., Chitsaz, F., Geer, L.Y., Geer, R.C., He, J., Gwadz, M., Hurwitz, D.I., et al. (2015) CDD: NCBI’s Conserved Domain Database. Nucleic Acids Research, 43, D222-D226.
 Marchler-Bauer, A., Zheng, C., Chitsaz, F., Derbyshire, M.K., Geer, L.Y., Geer, R.C., Gonzales, N.R., Gwadz, M., Hurwitz, D.I., Lanczycki, C.J., et al. (2013) CDD: Conserved Domains and Protein Three-Dimensional Structure. Nucleic Acids Research, 41, D348-D352.
 Jaroszewski, L., Rychlewski, L., Li, Z., Li, W. and Godzik, A. (2005) FFAS03: A Server for Profile-Profile Sequence Alignments. Nucleic Acids Research, 33, W284-W288.
 Dembo, A., Karlin, S. and Zeitouni, O. (1994) Limit Distribution of Maximal Non-Aligned Two-Sequence Segmental Score. The Annals of Probability, 22, 2022-2039.
 Ming, R., VanBuren, R., Liu, Y., Yang, M., Han, Y., Li, L.T., Zhang, Q., Kim, M.J., Schatz, M.C., Campbell, M., et al. (2013) Genome of the Long-Living Sacred Lotus (Nelumbo nucifera Gaertn.). Genome Biology, 14, R41.
 Vimolmangkang, S., Deng, X., Owiti, A., Meelaph, T., Ogutu, C. and Han, Y. (2016) Evolutionary Origin of the NCSI Gene Subfamily Encoding Norcoclaurine Synthase Is Associated with the Biosynthesis of Benzylisoquinoline Alkaloids in Plants. Scientific Reports, 6, Article No. 26323.
 Frago, E., Dicke, M. and Godfray, H.C. (2012) Insect Symbionts as Hidden Players in Insect-Plant Interactions. Trends in Ecology and Evolution, 27, 705-711.
 Zhu, S. and Gao, B. (2014) Nematode-Derived Drosomycin-Type Antifungal Peptides Provide Evidence for Plant-to-Ecdysozoan Horizontal Transfer of a Disease Resistance Gene. Nature Communications, 5, 3154.
 Low, E.T., Alias, H., Boon, S.H., Shariff, E.M., Tan, C.Y., Ooi, L.C., Cheah, S.C., Raha, A.R., Wan, K.L. and Singh, R. (2008) Oil Palm (Elaeis guineensis Jacq.) Tissue Culture ESTs: Identifying Genes Associated with Callogenesis and Embryogenesis. BMC Plant Biology, 8, Article No. 62.
 Uthaipaisanwong, P., Chanprasert, J., Shearman, J.R., Sangsrakru, D., Yoocha, T., Jomchai, N., Jantasuriyarat, C., Tragoonrung, S. and Tangphatsornruang, S. (2012) Characterization of the Chloroplast Genome Sequence of Oil Palm (Elaeis guineensis Jacq.). Gene, 500, 172-180.
 Bohlmann, J., Stauber, E.J., Krock, B., Oldham, N.J., Gershenzon, J. and Baldwin, I.T. (2002) Gene Expression of 5-Epi-Aristolochene Synthase and Formation of Capsidiol in Roots of Nicotiana attenuata and N. sylvestris. Phytochemistry, 60, 109-116.
 Yukawa, M., Tsudzuki, T. and Sugiura, M. (2006) The Chloroplast Genome of Nicotiana sylvestris and Nicotiana tomentosiformis: Complete Sequencing Confirms that the Nicotiana sylvestris Progenitor is the Maternal Genome Donor of Nicotiana tabacum. Molecular Genetics and Genomics, 275, 367-373.
 Sierro, N., Battey, J.N., Ouadi, S., Bovet, L., Goepfert, S., Bakaher, N., Peitsch, M.C. and Ivanov, N.V. (2013) Reference Genomes and Transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis. Genome Biology, 14, R60.
 Light, S., Sagit, R., Ithychanda, S.S., Qin, J. and Elofsson, A. (2012) The Evolution of Filamin-a Protein Domain Repeat Perspective. Journal of Structural Biology, 179, 289-298.
 Becker, B., Doan, J.M., Wustman, B., Carpenter, E.J., Chen, L., Zhang, Y., Wong, G.K. and Melkonian, M. (2015) The Origin and Evolution of the Plant Cell Surface: Algal Integrin-Associated Proteins and a New Family of Integrin-Like Cytoskeleton-ECM Linker Proteins. Genome Biology and Evolution, 7, 1580-1589.
 Chandran, T., Sharma, A. and Vijayan, M. (2013) Generation of Ligand Specificity and Modes of Oligomerization in β-Prism I Fold Lectins. Advances in Protein Chemistry and Structural Biology. 92, 135-178.
 Fallen, K., Banerjee, S., Sheehan, J., Addison, D., Lewis, L.M., Meiler, J. and Denton, J.S. (2009) The Kir Channel Immunoglobulin Domain Is Essential for Kir1.1 (ROMK) Thermodynamic Stability, Trafficking and Gating. Channels (Austin), 3, 57-68.
 Yereddi, N.R., Cusdin, F.S., Namadurai, S., Packman, L.C., Monie, T.P., Slavny, P., Clare, J.J., Powell, A.J. and Jackson, A.P. (2013) The Immunoglobulin Domain of the Sodium Channel β3 Subunit Contains a Surface-Localized Disulfide Bond that Is Required for Homophilic Binding. FASEB Journal, 27, 568-580.
 Ashkani, J. and Rees, D.J. (2016) A Comprehensive Study of Molecular Evolution at the Self-Incompatibility Locus of Rosaceae. Journal of Molecular Evolution, 82, 128-145.
 Nasrallah, J.B. (2017) Plant Mating Systems: Self-Incompatibility and Evolutionary Transitions to Self-Fertility in the Mustard Family. Current Opinion in Genetics and Development, 47, 54-60.
 Wilkins, K.A., Poulter, N.S. and Franklin-Tong, V.E. (2014) Taking One for the Team: Self-Recognition and Cell Suicide in Pollen. Journal of Experimental Botany, 65, 1331-1342.
 Berisio, R., Ciccarelli, L., Squeglia, F., De Simone, A. and Vitagliano, L. (2012) Structural and Dynamic Properties of Incomplete Immunoglobulin-Like Fold Domains. Protein and Peptide Letters, 19, 1045-1053.
 Lee, T.Y., Bretaña, N.A. and Lu, C.T. (2011) PlantPhos: Using Maximal Dependence Decomposition to Identify Plant Phosphorylation Sites with Substrate Site Specificity. BMC Bioinformatics, 12, Article No. 261.
 Yoshiyama, K.O., Kobayashi, J., Ogita, N., Ueda, M., Kimura, S., Maki, H. and Umeda, M. (2013) ATM-Mediated Phosphorylation of SOG1 Is Essential for the DNA Damage Response in Arabidopsis. EMBO Reports, 14, 817-822.
 Zulawski, M., Braginets, R. and Schulze, W.X. (2013) PhosPhAt Goes Kinases Searchable Protein Kinase Target Information in the Plant Phosphorylation Site Database PhosPhAt. Nucleic Acids Research, 41, D1176-D1184.
 Takagi, M., Sakamoto, T., Suzuki, R., Nemoto, K., Obayashi, T., Hirakawa, T., Matsunaga, T.M., Kurihara, D., Nariai, Y., Urano, T., Sawasaki, T. and Matsunaga, S. (2016) Plant Aurora Kinases Interact with and Phosphorylate Transcription Factors. Journal of Plant Research, 129, 1165-1178.
 Gnad, F., Forner, F., Zielinska, D.F., Birney, E., Gunawardena, J. and Mann, M. (2010) Evolutionary Constraints of Phosphorylation in Eukaryotes, Prokaryotes, and Mitochondria. Molecular and Cellular Proteomics, 9, 2642-2653.
 Hirano, M., Guo, P., McCurley, N., Schorpp, M., Das, S., Boehm, T. and Cooper, M.D. (2013) Evolutionary Implications of a Third Lymphocyte Lineage in Lampreys. Nature, 501, 435-438.
 Li, D., Deng, Z., Qin, B., Liu, X. and Men, Z. (2012) De Novo Assembly and Characterization of Bark Transcriptome Using Illumina Sequencing and Development of EST-SSR Markers in Rubber Tree (Hevea brasiliensis Muell. Arg.). BMC Genomics, 13, Article No. 192.
 Vázquez, N., Rocha, S., López-Fernández, H., Torres, A., Camacho, R., Fdez-Riverola, F., Vieira, J., Vieira, C.P. and Reboiro-Jato, M. (2019) EvoPPI 1.0: A Web Platform for Within- and Between-Species Multiple Interactome Comparisons and Application to Nine PolyQ Proteins Determining Neurodegenerative Diseases. Interdisciplinary Sciences, Computational Life Sciences, 11, 45-56.
 Hirayama, T. and Oka, A. (1992) Novel Protein Kinase of Arabidopsis thaliana (APK1) That Phosphorylates Tyrosine, Serine and Threonine. Plant Molecular Biology, 20, 653-662.
 Hanks, S.K., Quinn, A.M. and Hunter, T. (1988) The Protein Kinase Family: Conserved Features and Deduced Phylogeny of the Catalytic Domains. Science, 241, 42-52.
 Stancik, I.A., Šestak, M.S., Ji, B., Axelson-Fisk, M., Franjevic, D., Jers, C., Domazet-Lošo, T. and Mijakovic, I. (2018) Serine/Threonine Protein Kinases from Bacteria, archaea and Eukarya Share a Common Evolutionary Origin Deeply Rooted in the Tree of Life. Journal of Molecular Biology, 430, 27-32.
 Rudrabhatla, P., Reddy, M.M. and Rajasekharan, R. (2006) Genome-Wide Analysis and Experimentation of Plant Serine/Threonine/Tyrosine-Specific Protein Kinases. Plant Molecular Biology, 60, 293-319.
 Klaus-Heisen, D., Nurisso, A., Pietraszewska-Bogiel, A., Mbengue, M., Camut, S., Timmers, T., Pichereaux, C., Rossignol, M., Gadella, T.W., Imberty, A., Lefebvre, B. and Cullimore, J.V. (2011) Structure-Function Similarities between a Plant Receptor-Like Kinase and the Human Interleukin-1 Receptor-Associated Kinase-4. Journal of Biological Chemistry, 286, 11202-11210.
 King, N. and Carroll, S.B. (2001) A Receptor Tyrosine Kinase from Choanoflagellates: Molecular Insights into Early Animal Evolution. Proceedings of National Academy of Sciences of the United States of America, 98, 15032-15037.
 Schultheiss, K.P., Craddock, B.P., Tong, M., Seeliger, M. and Miller, W.T. (2013) Metazoan-Like Signaling in a Unicellular Receptor Tyrosine Kinase. BMC Biochemistry, 14, Article No. 4.
 Mohanty, S., Oruganty, K., Kwon, A., Byrne, D.P., Ferries, S., Ruan, Z., Hanold, L.E., Katiyar, S., Kennedy, E.J., Eyers, P.A. and Kannan, N. (2016) Hydrophobic Core Variations Provide a Structural Framework for Tyrosine Kinase Evolution and Functional Specialization. PLoS Genetics, 12, e1005885.
 Hakim, I., Amariglio, N., Grossman, Z., Simoni-Brok, F., Ohno, S. and Rechavi, G. (1994) The Genome of the I Human Transposable Repetitive Elements Is Composed of a Basic Motif Homologous to an Ancestral Immunoglobulin Gene Sequence. Proceedings of National Academy of Sciences of the United States of America, 91, 7967-7969.
 Pang, E. and Lin, K. (2010) Yeast Protein-Protein Interaction Binding Sites: Prediction from the Motif-Motif, Motif-Domain and Domain-Domain Levels. Molecular Biosystems, 6, 2164-2173.
 Wang, F., Liu, M., Song, B., Li, D., Pei, H., Guo, Y., Huang, J. and Zhang, D. (2012) Prediction and Characterization of Protein-Protein Interaction Networks in Swine. Proteome Science, 10, Article No. 2.
 Kubrycht, J., Sigler, K. and Souček, P. (2012) Virtual Interactomics of Proteins from Biochemical Standpoint. Molecular Biology International, 2012, Article ID: 976385.