Domain Architecture Analysis Research Articles

SMART: recent updates, new developments and status in 2020.

SMART (Simple Modular Architecture Research Tool) is a web resource (https://smart.embl.de) for the identification and annotation of protein domains and the analysis of protein domain architectures. SMART version 9 contains manually curated models for more than 1300 protein domains, with a topical set of 68 new models added since our last update article (1). All the new models are for diverse recombinase families and subfamilies and as a set they provide a comprehensive overview of mobile element recombinases namely transposase, integrase, relaxase, resolvase, cas1 casposase and Xer like cellular recombinase. Further updates include the synchronization of the underlying protein databases with UniProt (2), Ensembl (3) and STRING (4), greatly increasing the total number of annotated domains and other protein features available in architecture analysis mode. Furthermore, SMART’s vector-based protein display engine has been extended and updated to use the latest web technologies and the domain architecture analysis components have been optimized to handle the increased number of protein features available.

TREND: a platform for exploring protein function in prokaryotes based on phylogenetic, domain architecture and gene neighborhood analyses.

Key steps in a computational study of protein function involve analysis of (i) relationships between homologous proteins, (ii) protein domain architecture and (iii) gene neighborhoods the corresponding proteins are encoded in. Each of these steps requires a separate computational task and sets of tools. Currently in order to relate protein features and gene neighborhoods information to phylogeny, researchers need to prepare all the necessary data and combine them by hand, which is time-consuming and error-prone. Here, we present a new platform, TREND (tree-based exploration of neighborhoods and domains), which can perform all the necessary steps in automated fashion and put the derived information into phylogenomic context, thus making evolutionary based protein function analysis more efficient. A rich set of adjustable components allows a user to run the computational steps specific to his task. TREND is freely available at http://trend.zhulinlab.org.

Molecular characterization and induced expression analysis of the terminal complement component C9 in rohu,Labeo rohita

The complement component 9 (C9) plays a significant role in the formation of membrane attack complex (MAC) on the targeted cell surface. The current study is dealt with molecular characterization of C9 gene from rohu, Labeo rohita, an important cultured carp species in India. An open reading frame (ORF) of 1998 bp was amplified by polymerase chain reaction (PCR) that encodes a polypeptide of 666 amino acids having a signal peptide of 19 amino acids and a mature peptide of 647 amino acids. The SMART domain architecture analysis revealed two thrombospondin type-1 domains (TSP1), a low-density lipoprotein receptor domain class A (LDLa), a membrane attack complex and perforin (MACPF) domain, and an epidermal growth factor (EGF)-like domain. Multiple sequence alignment and evolutionary analysis revealed a primitive C9 sequence of rohu with maximum similarity and clustering with common carp, grass carp and zebrafish. C9 was highly expressed in liver and constitutively expressed in wide array of tissues except in eye of rohu juveniles. It was expressed during early developmental days of rohu including in milt. A variable level of up-regulation in C9 expression was noticed upon poly I:C induction, Aeromonas hydrophila and Argulus siamensis infections in liver, spleen and gill tissues of rohu at different time points. A constitutive expression of C9 in different stages of rohu during the ontogeny and in response to pathogen exposures along with high degree of sequence homology with other fish species proved it as an important primitive immune molecule of the complement system lytic pathway.

Small heat shock proteins in the amphibian Pelophylax bergeri: Cloning and characterization of Hsp27 and Hsp30 cDNAs and their expression analysis in ex vivo skin exposed to abiotic stresses

Small Heat Shock Proteins (sHSP) are molecular chaperones that play an essential role in maintaining protein homeostasis and promoting cell survival. In this work, for the first time, multiple cDNAs encoding for small Hsp27 and Hsp30, designed, respectively, as PbHsp27-(1-2) and PbHsp30-(1-5), were cloned and characterized in the amphibian Pelophylax bergeri, which is a suitable model for studying biological responses to environmental perturbations. Domain architecture analysis showed that PbHsp27 and PbHsp30 cDNAs displayed the typical signature motifs of the sHSP family such as the conserved α-crystallin domain flanked by variable N-terminal and C-terminal regions. Phylogenetic analysis revealed that PbHsp27 and PbHsp30 clustered, respectively, with Hsp27 and Hsp30 members of other vertebrates, but more closely with amphibians. Overall PbHsp27 and PbHsp30 transcriptional activity, analyzed by qRT-PCR, evidenced that, in ex vivo skin exposed to thermal shock and cadmium treatment, PbHsp27 and PbHsp30 mRNAs were inducible and regulated differently. This study provides the basis for future research on the potential use of PbHsp27 and PbHsp30 as biomarkers of proteotoxic stress in amphibians.

Changes in the Arabidopsis RNA-binding proteome reveal novel stress response mechanisms

BackgroundRNA-binding proteins (RBPs) are increasingly recognized as regulatory component of post-transcriptional gene expression. RBPs interact with mRNAs via RNA-binding domains and these interactions affect RNA availability for translation, RNA stability and turn-over thus affecting both RNA and protein expression essential for developmental and stimulus specific responses. Here we investigate the effect of severe drought stress on the RNA-binding proteome to gain insights into the mechanisms that govern drought stress responses at the systems level.ResultsLabel-free mass spectrometry enabled the identification 567 proteins of which 150 significantly responded to the drought-induced treatment. A gene ontology analysis revealed enrichment in the “RNA binding” and “RNA processing” categories as well as biological processes such as “response to abscisic acid” and “response to water deprivation”. Importantly, a large number of the stress responsive proteins have not previously been identified as RBPs and include proteins in carbohydrate metabolism and in the glycolytic and citric acid pathways in particular. This suggests that RBPs have hitherto unknown roles in processes that govern metabolic changes during stress responses. Furthermore, a comparative analysis of RBP domain architectures shows both, plant specific and common domain architectures between plants and animals. The latter could be an indication that RBPs are part of an ancient stress response.ConclusionThis study establishes mRNA interactome capture technique as an approach to study stress signal responses implicated in environmental changes. Our findings denote RBP changes in the proteome as critical components in plant adaptation to changing environments and in particular drought stress protein-dependent changes in RNA metabolism.

Two cyp17 genes perform different functions in the sex hormone biosynthesis and gonadal differentiation in Japanese flounder (Paralichthys olivaceus)

P450c17, a key enzyme in the steroid generation pathway, plays an important role in the production of sex steroid and cortisol. In this study, two cyp17 gene isoforms, Pocyp17-I and Pocyp17-II were isolated from Paralichthys olivaceus gonads. Domain architecture analysis of Pocyp17-I and Pocyp17-II revealed that they had three regions important to enzymatic function. Structural analysis showed that Pocyp17-I and Pocyp17-II had 8 and 9 exons respectively, and the difference was caused by the insertion of an extra intron (intron1) in the latter. Quantitative real-time polymerase chain reaction results indicated that the expression of these two genes showed sexually dimorphism that Pocyp17-I and Pocyp17-II were highest expressed in testis and ovary, respectively. The in situ hybridization analysis of gonads indicated that Pocyp17-I and Pocyp17-II mRNA were both detected in oocytes, spermatocytes and Sertoli cells. After injection of androgen and estrogen (17α-methyltestosterone, 17β-estradiol) of different concentrations, the expression level of Pocyp17-I decreased significantly (P < 0.01), whereas estrogen had no influence on Pocyp17-II, but androgen upregulated the expression of Pocyp17-II (P < 0.05). Moreover, Pocyp17-I expression level was down-regulated significantly by NR0b1 but up-regulated by NR5a2 (P < 0.05), whereas Pocyp17-II expression level was down-regulated significantly by NR0b1 and NR5a2 (P < 0.05). All these results demonstrated that there were differences in expression patterns, feedback actions of sex hormones and transcriptional regulations between cyp17-I and cyp17-II, which revealed that cyp17-I and cyp17-II might perform different functions in sex hormones biosynthesis and gonadal differentiation in Japanese flounder.

Messages From the Past: New Insights in Plant Lectin Evolution.

Lectins are a large and diverse class of proteins, found in all kingdoms of life. Plants are known to express different types of carbohydrate-binding proteins, each containing at least one particular lectin domain which enables them to specifically recognize and bind carbohydrate structures. The group of plant lectins is heterogeneous in terms of structure, biological activity and function. Lectins control various aspects of plant development and defense. Some lectins facilitate recognition of exogenous danger signals or play a role in endogenous signaling pathways, while others are considered as storage proteins or involved in symbiotic relationships. In this study, we revisit the origin of the different plant lectin families in view of the recently reshaped tree of life. Due to new genomic sampling of previously unknown microbial lineages, the tree of life has expanded and was reshaped multiple times. In addition, more plant genomes especially from basal Phragmoplastophyta, bryophytes, and Salviniales (e.g., Chara braunii, Marchantia polymorpha, Physcomitrella patens, Azolla filiculoides, and Salvinia cucullata) have been analyzed, and annotated genome sequences have become accessible. We searched 38 plant genome sequences including core eudicots, monocots, gymnosperms, fern, lycophytes, bryophytes, charophytes, chlorophytes, glaucophytes, and rhodophytes for lectin motifs, performed an extensive comparative analysis of lectin domain architectures, and determined the phylogenetic and evolutionary history of lectins in the plant lineage. In conclusion, we describe the conservation of particular domains in plant lectin sequences obtained from algae to higher plants. The strong conservation of several lectin motifs highlights their significance for plants.

Classification of human Herpesviridae proteins using Domain-architecture Aware Inference of Orthologs (DAIO)

We developed a computational approach called Domain-architecture Aware Inference of Orthologs (DAIO) for the analysis of protein orthology by combining phylogenetic and protein domain-architecture information. Using DAIO, we performed a systematic study of the proteomes of all human Herpesviridae species to define Strict Ortholog Groups (SOGs). In addition to assessing the taxonomic distribution for each protein based on sequence similarity, we performed a protein domain-architecture analysis for every protein family and computationally inferred gene duplication events. While many herpesvirus proteins have evolved without any detectable gene duplications or domain rearrangements, numerous herpesvirus protein families do exhibit complex evolutionary histories. Some proteins acquired additional domains (e.g., DNA polymerase), whereas others show a combination of domain acquisition and gene duplication (e.g., betaherpesvirus US22 family), with possible functional implications. This novel classification system of SOGs for human Herpesviridae proteins is available through the Virus Pathogen Resource (ViPR, www.viprbrc.org).

Genome-wide characterization of respiratory burst oxidase homolog (RBOH) in soybean [(glycine max l. Merr.)

Plant respiratory burst oxidase homologs (RBOHs), also known as NADPH oxidase (NOX) have involved as a key producers of reactive oxygen species (ROS) in a wide range of biological processes, including plant growth, development, and abiotic stress responses. In this work, a total 17 soybean RBOH encoding genes (GmRBOHs) were identified. Furthermore, their domain architecture, chromosomal distribution, gene structure, gene duplication, and phylogenetic analysis were carried out. Phylogenetic tree divided the RBOHs from soybean, another legume plants (Phaseolus vulgaris), and Arabidopsis into six subfamilies, with soybean RBOH genes distributed in subfamily I to VI. All RBOH proteins of soybean contained four major domains, which focused on EF-hands domains with Ca2+-binding residues provided more insights into function of RBOH proteins in soybean. The differential expression pattern using RNA-seq of 17 RBOH genes in soybean indicated that most of the genes express in root and nodule; among them, GmRboh9 and GmRboh12 have the highest expression. The analysis presented provides a scientific basis for further research on function of RBOH genes in soybean.

In silico analysis of onion chitinases using transcriptome data.

Chitinases are glycoside hydrolase (GH) family of proteins having multifaceted roles in plants. It is of interest to identify and characterize chitinase-encoding genes from the popular bulbous plant onion (Allium cepa L.). We have used the EST sequences for onion chitinases to elucidate its functional features using sequence, structure and functional analysis. These contigs belong to the GH19 chitinases family according to domain architecture analysis. They have highly conserved chitinase motifs including motifs exclusive to plant chitinases as implied using the MEME based structural characterization. Estimation of biochemical properties suggested that these proteins have features to form stable and hydrophilic proteins capable of localizing extracellular and in vacuoles. Further, they have multiple cellular processes including defense role as inferred by DeepGO function prediction. Phylogenetic analysis grouped them as class I and class VII plant chitinase, with possible abundance of class I chitinase in onion. These observations help in the isolation and functional validation of onion chitinases.

Genome-wide identification of chitinase and chitin deacetylase gene families in the oriental fruit fly, Bactrocera dorsalis (Hendel)

Chitinases (Chts) and chitin deacetylases (CDAs) are important enzymes required for chitin metabolism in insects. In this study, 12 Cht-related genes (including seven Cht genes and five imaginal disc growth factor genes) and 6 CDA genes (encoding seven proteins) were identified in Bactrocera dorsalis using genome-wide searching and transcript profiling. Based on the conserved sequences and phylogenetic relationships, 12 Cht-related proteins were clustered into eight groups (group I–V and VII–IX). Further domain architecture analysis showed that all contained at least one chitinase catalytic domain, however, only four (BdCht5, BdCht7, BdCht8 and BdCht10) possessed chitin-binding domains. The subsequent phylogenetic analysis revealed that seven CDAs were clustered into five groups (group I–V), and all had one chitin deacetylase catalytic domain. However, only six exhibited chitin-binding domains. Finally, the development- and tissue-specific expression profiling showed that transcript levels of the 12 Cht-related genes and 6 CDA genes varied considerably among eggs, larvae, pupae and adults, as well as among different tissues of larvae and adults. Our findings illustrate the structural differences and expression patterns of Cht and CDA genes in B. dorsalis, and provide important information for the development of new pest control strategies based on these vital enzymes.

Genome-wide survey of remote homologues for protein domain superfamilies of known structure reveals unequal distribution across structural classes.

Domains are the basic building blocks of proteins which can combine to give rise to different domain architectures. Annotation of domains in a sequence is the first step towards understanding the biological function. Since there are a limited number of folds and evolutionarily related proteins have a similar structure, function can be inferred through remote homology. Computational sequence searches were performed for remote homologues on genomes of around ∼160 000 different organisms, starting from nearly 11 000 superfamily queries of known structure. Case studies revealed that most of the associated domains are involved in the same biological process. Using all the proteins predicted to have at least one structural domain, a coverage of 61% of Pfam families was achieved which is higher than the existing methods (43.36% by SIFTS). Taxonomic analysis of the proteins revealed 493 superfamilies in all the major kingdoms of life and a few lateral gene transfers between viruses and cellular organisms. The distribution of remote homologues across different classes, folds and superfamilies was studied and reveals that sequences are unequally distributed across structural classes. Finally, domain architectures were computed for the homologues and these data were compiled for each superfamily and organism.

Phylogenetic analysis of the core histone doublet and DNA topo II genes of Marseilleviridae: evidence of proto-eukaryotic provenance

BackgroundWhile the genomes of eukaryotes and Archaea both encode the histone-fold domain, only eukaryotes encode the core histone paralogs H2A, H2B, H3, and H4. With DNA, these core histones assemble into the nucleosomal octamer underlying eukaryotic chromatin. Importantly, core histones for H2A and H3 are maintained as neofunctionalized paralogs adapted for general bulk chromatin (canonical H2 and H3) or specialized chromatin (H2A.Z enriched at gene promoters and cenH3s enriched at centromeres). In this context, the identification of core histone-like “doublets” in the cytoplasmic replication factories of the Marseilleviridae (MV) is a novel finding with possible relevance to understanding the origin of eukaryotic chromatin. Here, we analyze and compare the core histone doublet genes from all known MV genomes as well as other MV genes relevant to the origin of the eukaryotic replisome.ResultsUsing different phylogenetic approaches, we show that MV histone domains encode obligate H2B-H2A and H4-H3 dimers of possible proto-eukaryotic origin. MV core histone moieties form sister clades to each of the four eukaryotic clades of canonical and variant core histones. This suggests that MV core histone moieties diverged prior to eukaryotic neofunctionalizations associated with paired linear chromosomes and variant histone octamer assembly. We also show that MV genomes encode a proto-eukaryotic DNA topoisomerase II enzyme that forms a sister clade to eukaryotes. This is a relevant finding given that DNA topo II influences histone deposition and chromatin compaction and is the second most abundant nuclear protein after histones.ConclusionsThe combined domain architecture and phylogenomic analyses presented here suggest that a primitive origin for MV histone genes is a more parsimonious explanation than horizontal gene transfers + gene fusions + sufficient divergence to eliminate relatedness to eukaryotic neofunctionalizations within the H2A and H3 clades without loss of relatedness to each of the four core histone clades. We thus suggest MV histone doublet genes and their DNA topo II gene possibly were acquired from an organism with a chromatinized replisome that diverged prior to the origin of eukaryotic core histone variants for H2/H2A.Z and H3/cenH3. These results also imply that core histones were utilized ancestrally in viral DNA compaction and/or protection from host endonucleases.

Protein kinase KIN10 from Arabidopsis thaliana as a potential regulator of primary microtubule nucleation centers in plants

The nearest plant homologues of animal protein kinase BRSK were identified using the methods of classical and structural bioinformatics. The selection was performed based on the sequence comparison, results of phylogenetic clustering, and analysis of domain architecture. Spatial structures of human BRSK1 and KIN10 from A. thaliana were compared. The relationship between KIN10 and the regulation of primary microtubule nucleation centers in A. thaliana was revealed. Obvious homology of plant KIN10 and mammalian BRSK1 evidence to suggest that this plant protein kinase is associated with the regulation of the structure and function of primary microtubule nucleation centers and is able to phosphorylate γ-tubulin from Arabidopsis (TUBG1 and TUBG2) at Ser131, affecting the γTuSC monomer structure as well as the γTuRC complex assembly. The effect of the modification on the TUBG1-GACP3 interaction was suggested.

20 years of the SMART protein domain annotation resource.

SMART (Simple Modular Architecture Research Tool) is a web resource (http://smart.embl.de) for the identification and annotation of protein domains and the analysis of protein domain architectures. SMART version 8 contains manually curated models for more than 1300 protein domains, with approximately 100 new models added since our last update article (1). The underlying protein databases were synchronized with UniProt (2), Ensembl (3) and STRING (4), doubling the total number of annotated domains and other protein features to more than 200 million. In its 20th year, the SMART analysis results pages have been streamlined again and its information sources have been updated. SMART’s vector based display engine has been extended to all protein schematics in SMART and rewritten to use the latest web technologies. The internal full text search engine has been redesigned and updated, resulting in greatly increased search speed.

Systematic classification of the His-Me finger superfamily

The His-Me finger endonucleases, also known as HNH or ββα-metal endonucleases, form a large and diverse protein superfamily. The His-Me finger domain can be found in proteins that play an essential role in cells, including genome maintenance, intron homing, host defense and target offense. Its overall structural compactness and non-specificity make it a perfectly-tailored pathogenic module that participates on both sides of inter- and intra-organismal competition. An extremely low sequence similarity across the superfamily makes it difficult to identify and classify new His-Me fingers. Using state-of-the-art distant homology detection methods, we provide an updated and systematic classification of His-Me finger proteins. In this work, we identified over 100 000 proteins and clustered them into 38 groups, of which three groups are new and cannot be found in any existing public domain database of protein families. Based on an analysis of sequences, structures, domain architectures, and genomic contexts, we provide a careful functional annotation of the poorly characterized members of this superfamily. Our results may inspire further experimental investigations that should address the predicted activity and clarify the potential substrates, to provide more detailed insights into the fundamental biological roles of these proteins.

Genome-wide identification and expression analyses of the homeobox transcription factor family during ovule development in seedless and seeded grapes

Seedless grapes are of considerable importance for the raisin and table grape industries. Previous transcriptome analyses of seed development in grape revealed that genes encoding homeobox transcription factors were differentially regulated in seedless compared with seeded grape during seed development. In the present study, we identified a total of 73 homeobox-like genes in the grapevine genome and analyzed the genomic content and expression profiles of these genes. Based on domain architecture and phylogenetic analyses grape homeobox genes can be classified into eleven subfamilies. An analysis of the exon-intron structures and conserved motifs provided further insight into the evolutionary relationships between these genes. Evaluation of synteny indicated that segmental and tandem duplications have contributed greatly to the expansion of the grape homeobox gene superfamily. Synteny analysis between the grape and Arabidopsis genomes provided a potential functional relevance for these genes. The tissue-specific expression patterns of homeobox genes suggested roles in both vegetative and reproductive tissues. Expression profiling of these genes during the course of ovule development in seeded and seedless cultivars suggested a potential role in ovule abortion associated with seedlessness. This study will facilitate the functional analysis of these genes and provide new resources for molecular breeding of seedless grapes.

Identification of a homogenous structural basis for oligomerization by retroviral Rev-like proteins

BackgroundRev-like proteins are post-transcriptional regulatory proteins found in several retrovirus genera, including lentiviruses, betaretroviruses, and deltaretroviruses. These essential proteins mediate the nuclear export of incompletely spliced viral RNA, and act by tethering viral pre-mRNA to the host CRM1 nuclear export machinery. Although all Rev-like proteins are functionally homologous, they share less than 30% sequence identity. In the present study, we computationally assessed the extent of structural homology among retroviral Rev-like proteins within a phylogenetic framework.ResultsWe undertook a comprehensive analysis of overall protein domain architecture and predicted secondary structural features for representative members of the Rev-like family of proteins. Similar patterns of α-helical domains were identified for Rev-like proteins within each genus, with the exception of deltaretroviruses, which were devoid of α-helices. Coiled-coil oligomerization motifs were also identified for most Rev-like proteins, with the notable exceptions of HIV-1, the deltaretroviruses, and some small ruminant lentiviruses. In Rev proteins of primate lentiviruses, the presence of predicted coiled-coil motifs segregated within specific primate lineages: HIV-1 descended from SIVs that lacked predicted coiled-coils in Rev whereas HIV-2 descended from SIVs that contained predicted coiled-coils in Rev. Phylogenetic ancestral reconstruction of coiled-coils for all Rev-like proteins predicted a single origin for the coiled-coil motif, followed by three losses of the predicted signal. The absence of a coiled-coil signal in HIV-1 was associated with replacement of canonical polar residues with non-canonical hydrophobic residues. However, hydrophobic residues were retained in the key ‘a’ and ‘d’ positions, and the α-helical region of HIV-1 Rev oligomerization domain could be modeled as a helical wheel with two predicted interaction interfaces. Moreover, the predicted interfaces mapped to the dimerization and oligomerization interfaces in HIV-1 Rev crystal structures. Helical wheel projections of other retroviral Rev-like proteins, including endogenous sequences, revealed similar interaction interfaces that could mediate oligomerization.ConclusionsSequence-based computational analyses of Rev-like proteins, together with helical wheel projections of oligomerization domains, reveal a conserved homogeneous structural basis for oligomerization by retroviral Rev-like proteins.

Genome-Wide Screening for Lectin Motifs in Arabidopsis thaliana.

For more than three decades, served as a model for plant biology research. At present only a few protein families have been studied in detail in . This study focused on all sequences with lectin motifs in the genome of . Based on amino acid sequence similarity (BLASTp searches), 217 putative lectin genes were retrieved belonging to 9 out of 12 different lectin families. The domain organization and genomic distribution for each lectin family were analyzed. Domain architecture analysis revealed that most of these lectin gene sequences are linked to other domains, often belonging to protein families with catalytic activity. Many protein domains identified are known to play a role in stress signaling and defense, suggesting a major contribution of the putative lectins in development and plant defense. This genome-wide screen for different lectin motifs will help to unravel the functional characteristics of lectins. In addition, phylogenetic trees and WebLogos were created and showed that most lectin sequences that share the same domain architecture evolved together. Furthermore, the amino acids responsible for carbohydrate binding are largely conserved. Our results provide information about the evolutionary relationships and functional divergence of the lectin motifs in .

Genome-wide analysis of lectin receptor-like kinases in Populus.

BackgroundReceptor-like kinases (RLKs) belong to a large protein family with over 600 members in Arabidopsis and over 1000 in rice. Among RLKs, the lectin receptor-like kinases (LecRLKs) possess a characteristic extracellular carbohydrate-binding lectin domain and play important roles in plant development and innate immunity. There are 75 and 173 LecRLKs in Arabidopsis and rice, respectively. However, little is known about LecRLKs in perennial woody plants.ResultsHere we report the genome-wide analysis of classification, domain architecture and expression of LecRLKs in the perennial woody model plant Populus. We found that the LecRLK family has expanded in Populus to a total of 231, including 180 G-type, 50 L-type and 1 C-type LecRLKs. Expansion of the Populus LecRLKs (PtLecRLKs) occurred partially through tandem duplication. Based on domain architecture and orientation features, we classified PtLecRLKs into eight different classes. RNA-seq-based transcriptomics analysis revealed diverse expression patterns of PtLecRLK genes among leaves, stems, roots, buds and reproductive tissues and organs.ConclusionsThis study offers a comprehensive view of LecRLKs in the perennial woody model plant Populus and provides a foundation for functional characterization of this important family of receptor-like kinases.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3026-2) contains supplementary material, which is available to authorized users.