Kinannote, a computer program to identify and classify members of the eukaryotic protein kinase superfamily

KiNext: a portable and scalable workflow for the identification and classification of protein kinases

BACKGROUND

Protein kinases are a diverse superfamily of proteins common to organisms across the tree of life that are typically involved in signal transduction, allowing organisms to sense and respond to biotic or abiotic environmental factors. They have important roles in organismal physiology, including development, reproduction, acclimation to environmental stress, while their dysregulation can lead to disease, including several forms of cancer. Identifying the complement of protein kinases (the kinome) of any organism is useful for understanding its physiological capabilities, limitations and adaptations to environmental stress. The increasing availability of genomes makes it now possible to examine and compare the kinomes across a broad diversity of organisms. Here we present a pipeline respecting the FAIR principles (findable, accessible, interoperable and reusable) that facilitates the search and identification of protein kinases from a predicted proteome, and classifies them according to group of serine/threonine/tyrosine protein kinases present in eukaryotes.

RESULTS

KiNext is a Nextflow pipeline that regroups a number of existing bioinformatic tools to search for and classify the protein kinases of an organism in a reproducible manner, starting from a set of amino acid sequences. Conventional eukaryotic protein kinases (ePKs) and atypical protein kinases (aPKs) are identified by using Hidden Markov Models (HMMs) generated from the catalytic domains of kinases. Furthermore, KiNext categorizes ePKs into the eight kinase groups by employing dedicated Hidden Markov Models (HMMs) tailored for each group. The performance of the KiNext pipeline was validated against previously identified kinomes obtained with other tools that were already published for two marine species, the Pacific oyster Crassostrea gigas and the unicellular green alga Ostreoccocus tauri. KiNext outperformed previous results by finding previously unidentified kinases and by attributing a large proportion of previously unclassified kinases to a group in both species. These results demonstrate improvements in kinase identification and classification, all while providing traceability and reproducibility of results in a FAIR pipeline. The default HMM models provided with KiNext are most suitable for eukaryotes, but the pipeline can be easily modified to include HMM models for other taxa of interest.

CONCLUSION

The KiNext pipeline enables efficient and reproducible identification of kinomes based on predicted amino acid sequences (i.e. proteomes). KiNext was designed to be easy to use, automated, portable and scalable.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Identifying potential drug targets in the kinomes of two monogenean species

Protein kinases are enzymes involved in essential biological processes such as signal transduction, transcription, metabolism, and the cell cycle. Human kinases are targets for several drugs approved by the US Food and Drug Administration. Therefore, the identification and classification of kinases in other organisms, including pathogenic parasites, is an interesting subject of study. Monogeneans are platyhelminths, mainly ectoparasites, capable of causing health problems in farmed fish. Although some genomes and transcriptomes are available for monogenean species, their full repertoire of kinases is unknown. The aim of this study was to identify and classify the putative kinases in the transcriptomes of two monogeneans, Rhabdosynochus viridisi and Scutogyrus longicornis, and then to predict potential monogenean drug targets (MDTs) and selective inhibitor drugs using computational approaches. Monogenean kinases having orthologs in the lethal phenotype of C. elegans but not in fish or humans were considered MDTs. A total of 160 and 193 kinases were identified in R. viridisi and S. longicornis, respectively. Of these, 22 kinases, belonging mainly to the major groups CAMK, AGC, and TK, were classified as MDTs, five of which were evaluated further. Molecular docking analysis indicated that dihydroergotamine, ergotamine, and lomitapide have the highest affinity for the kinases BRSK and MEKK1. These well-known drugs could be evaluated in future studies for potential repurposing as anti-monogenean agents. The present study contributes valuable data for the development of new antiparasitic candidates for finfish aquaculture.

Dark kinase annotation, mining, and visualization using the Protein Kinase Ontology

The Protein Kinase Ontology (ProKinO) is an integrated knowledge graph that conceptualizes the complex relationships among protein kinase sequence, structure, function, and disease in a human and machine-readable format. In this study, we have significantly expanded ProKinO by incorporating additional data on expression patterns and drug interactions. Furthermore, we have developed a completely new browser from the ground up to render the knowledge graph visible and interactive on the web. We have enriched ProKinO with new classes and relationships that capture information on kinase ligand binding sites, expression patterns, and functional features. These additions extend ProKinO's capabilities as a discovery tool, enabling it to uncover novel insights about understudied members of the protein kinase family. We next demonstrate the application of ProKinO. Specifically, through graph mining and aggregate SPARQL queries, we identify the p21-activated protein kinase 5 (PAK5) as one of the most frequently mutated dark kinases in human cancers with abnormal expression in multiple cancers, including a previously unappreciated role in acute myeloid leukemia. We have identified recurrent oncogenic mutations in the PAK5 activation loop predicted to alter substrate binding and phosphorylation. Additionally, we have identified common ligand/drug binding residues in PAK family kinases, underscoring ProKinO's potential application in drug discovery. The updated ontology browser and the addition of a web component, ProtVista, which enables interactive mining of kinase sequence annotations in 3D structures and Alphafold models, provide a valuable resource for the signaling community. The updated ProKinO database is accessible at https://prokino.uga.edu.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Integrated analysis of SR-like protein kinases Sky1 and Sky2 links signaling networks with transcriptional regulation in Candida albicans

Fungal infections are a major global health burden where Candida albicans is among the most common fungal pathogen in humans and is a common cause of invasive candidiasis. Fungal phenotypes, such as those related to morphology, proliferation and virulence are mainly driven by gene expression, which is primarily regulated by kinase signaling cascades. Serine-arginine (SR) protein kinases are highly conserved among eukaryotes and are involved in major transcriptional processes in human and S. cerevisiae. Candida albicans harbors two SR protein kinases, while Sky2 is important for metabolic adaptation, Sky1 has similar functions as in S. cerevisiae. To investigate the role of these SR kinases for the regulation of transcriptional responses in C. albicans, we performed RNA sequencing of sky1Δ and sky2Δ and integrated a comprehensive phosphoproteome dataset of these mutants. Using a Systems Biology approach, we study transcriptional regulation in the context of kinase signaling networks. Transcriptomic enrichment analysis indicates that pathways involved in the regulation of gene expression are downregulated and mitochondrial processes are upregulated in sky1Δ. In sky2Δ, primarily metabolic processes are affected, especially for arginine, and we observed that arginine-induced hyphae formation is impaired in sky2Δ. In addition, our analysis identifies several transcription factors as potential drivers of the transcriptional response. Among these, a core set is shared between both kinase knockouts, but it appears to regulate different subsets of target genes. To elucidate these diverse regulatory patterns, we created network modules by integrating the data of site-specific protein phosphorylation and gene expression with kinase-substrate predictions and protein-protein interactions. These integrated signaling modules reveal shared parts but also highlight specific patterns characteristic for each kinase. Interestingly, the modules contain many proteins involved in fungal morphogenesis and stress response. Accordingly, experimental phenotyping shows a higher resistance to Hygromycin B for sky1Δ. Thus, our study demonstrates that a combination of computational approaches with integration of experimental data can offer a new systems biological perspective on the complex network of signaling and transcription. With that, the investigation of the interface between signaling and transcriptional regulation in C. albicans provides a deeper insight into how cellular mechanisms can shape the phenotype.

Deciphering the kinome of Theileria annulata for identification of drug targets and anti-theilerial drug

Tropical theileriosis is one of the major parasitic diseases of ruminants. It is a tick-borne disease caused by an apicomplexan parasite, Theileria annulata. In the infected cells, these parasites induce phenotypes similar to cancerous cells. Among the most critical changes induced by the parasite are immortalization, hyperproliferation, and dissemination. The proliferative signal in the T. annulata transformed cells are provided by different kinases such as mitogen-activated protein kinases, SRC family kinases, casein kinase-2, and phosphatidylinositide 3-kinase. Deregulation of protein kinases in cancer is also well known. Targeting protein kinases in a cancerous cell is one of the most common methods in cancer therapy. Here, we revisited the kinome of T. annulata and studied its evolutionary relationship with other piroplasms. This analysis revealed that T. annulata kinome encodes 54 protein kinases. Based on our analysis, 12 of these 54 kinases were identified for the first time in the T. annulata proteome. Three protein kinases, TA16570, TA09820, and TA07000, had <40% identity with Bos taurus and >40% identity with the previously identified potential drug targets present in the Therapeutic Target Database (TTD). These 3 proteins were predicted to be essential for the survival of T. annulata and were selected as drug targets. Screening these drug targets in the Protein Kinase Inhibitor Database (PKID) led to shortlisting of 5 drugs. Only Dabrafenib, out of these 5 drugs, could bind to the ATP binding site (in silico) of the Calcium Dependent Protein Kinase 3 of both T. annulata and Theileria parva. Further, dabrafenib could inhibit the proliferation of T. annulata infected bovine leucocytes in 6 days proliferation assay with the IC50 value of 0.66 µM. Also, this drug did not have a cytotoxic effect on bovine peripheral blood mononuclear cells. In summary, the analysis of T. annulata kinome led to the identification of dabrafenib as a potential drug for treating theileriosis.

Genomic analysis reveals cryptic diversity in aphelids and sheds light on the emergence of Fungi

Over the past decade, molecular phylogenetics has reshaped our understanding of the fungal tree of life by unraveling a hitherto elusive diversity of the protistan relatives of Fungi. Aphelida constitutes one of these novel deep branches that precede the emergence of osmotrophic fungal lifestyle and hold particular significance as the pathogens of algae. Here, we obtain and analyze the genomes of aphelid species Amoeboaphelidium protococcarum and Amoeboaphelidium occidentale. Genomic data unmask the vast divergence between these species, hidden behind their morphological similarity, and reveal hybrid genomes with a complex evolutionary history in two strains of A. protococcarum. We confirm the proposed sister relationship between Aphelida and Fungi using phylogenomic analysis and chart the reduction of characteristic proteins involved in phagocytic activity in the evolution of Holomycota. Annotation of aphelid genomes demonstrates the retention of actin nucleation-promoting complexes associated with phagocytosis and amoeboid motility and also reveals a conspicuous expansion of receptor-like protein kinases, uncharacteristic of fungal lineages. We find that aphelids possess multiple carbohydrate-processing enzymes that are involved in fungal cell wall synthesis but do not display rich complements of algal cell-wall-processing enzymes, suggesting an independent origin of fungal plant-degrading capabilities. Aphelid genomes show that the emergence of Fungi from phagotrophic ancestors relied on a common cell wall synthetic machinery but required a different set of proteins for digestion and interaction with the environment.

Update and elucidation of Plasmodium kinomes: Prioritization of kinases as potential drug targets for malaria

Malaria is a tropical disease caused by Plasmodium spp. and transmitted by the bite of infected Anopheles mosquitoes. Protein kinases (PKs) play key roles in the life cycle of the etiological agent of malaria, turning these proteins attractive targets for antimalarial drug discovery campaigns. As part of an effort to understand parasite signaling functions, we report the results of a bioinformatics pipeline analysis of PKs of eight Plasmodium species. To date, no P. malariae and P. ovale kinome assemble has been conducted. We classified, curated and annotated predicted kinases to update P. falciparum, P. vivax, P. yoelii, P. berghei, P. chabaudi, and P. knowlesi kinomes published to date, as well as report for the first time the kinomes of P. malariae and P. ovale. Overall, from 76 to 97 PKs were identified among all Plasmodium spp. kinomes. Most of the kinases were assigned to seven of nine major kinase groups: AGC, CAMK, CMGC, CK1, STE, TKL, OTHER; and the Plasmodium-specific group FIKK. About 30% of kinases have been deeply classified into group, family and subfamily levels and only about 10% remained unclassified. Furthermore, updating and comparing the kinomes of P. vivax and P. falciparum allowed for the prioritization and selection of kinases as potential drug targets that could be explored for discovering new drugs against malaria. This integrated approach resulted in the selection of 37 protein kinases as potential targets and the identification of investigational compounds with moderate in vitro activity against asexual P. falciparum (3D7 and Dd2 strains) stages that could serve as starting points for the search of potent antimalarial leads in the future.

Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

Deciphering the Crosstalk Mechanisms of Wheat-Stem Rust Pathosystem: Genome-Scale Prediction Unravels Novel Host Targets

Triticum aestivum (wheat), a major staple food grain, is affected by various biotic stresses. Among these, fungal diseases cause about 15-20% of yield loss, worldwide. In this study, we performed a comparative analysis of protein-protein interactions between two Puccinia graminis races (Pgt 21-0 and Pgt Ug99) that cause stem (black) rust in wheat. The available molecular techniques to study the host-pathogen interaction mechanisms are expensive and labor-intensive. We implemented two computational approaches (interolog and domain-based) for the prediction of PPIs and performed various functional analysis to determine the significant differences between the two pathogen races. The analysis revealed that T. aestivum-Pgt 21-0 and T. aestivum-Pgt Ug99 interactomes consisted of ∼90M and ∼56M putative PPIs, respectively. In the predicted PPIs, we identified 115 Pgt 21-0 and 34 Pgt Ug99 potential effectors that were highly involved in pathogen virulence and development. Functional enrichment analysis of the host proteins revealed significant GO terms and KEGG pathways such as O-methyltransferase activity (GO:0008171), regulation of signal transduction (GO:0009966), lignin metabolic process (GO:0009808), plastid envelope (GO:0009526), plant-pathogen interaction pathway (ko04626), and MAPK pathway (ko04016) that are actively involved in plant defense and immune signaling against the biotic stresses. Subcellular localization analysis anticipated the host plastid as a primary target for pathogen attack. The highly connected host hubs in the protein interaction network belonged to protein kinase domain including Ser/Thr protein kinase, MAPK, and cyclin-dependent kinase. We also identified 5,577 transcription factors in the interactions, associated with plant defense during biotic stress conditions. Additionally, novel host targets that are resistant to stem rust disease were also identified. The present study elucidates the functional differences between Pgt 21-0 and Pgt Ug99, thus providing the researchers with strain-specific information for further experimental validation of the interactions, and the development of durable, disease-resistant crop lines.

Kinome analyses of Candida albicans, C. parapsilosis and C. tropicalis enable novel kinases as therapeutic drug targets in candidiasis

Candida spp. have attracted considerable attention as they cause serious human diseases in immunocompromised individuals. The genomes of the pathogenic Candida spp. have been sequenced, but systemic characterizations of their kinomes are yet to be reported. As in various eukaryotes, the protein kinases play crucial regulatory roles in pathogenicity of Candida. Increased frequency of antifungal resistance in Candida spp. requires significant attention to explore novel therapeutic molecules for their control. The present in-silico study involves novel bioinformatics strategies to identify the kinase proteins and their potential drug targets with the purpose to combat fungal infections. The study reports 103, 107 and 106 kinase proteins from 3 Candida spp., C. albicans, C. parapsilosis and C. tropicalis, respectively. Moreover, 79 common kinase proteins were identified, of which 54 proteins play essential roles in Candida spp. and 42 proteins were human non-homologues. Among the essential and human non-homologous protein kinases, 9 were found to be common essential human non-homologues, of which 6 are uniquely present in Candida. These 6 protein kinases namely, Hsl1, Npr1, Ptk2, Kin2, Ksp1 and orf19.3854 (CAALFM_CR06040WA) are involved in various molecular and cellular processes regulating virulence or pathogenicity. Further, these 6 kinases are prioritized as potential drug targets and explored for discovering new lead compounds against candidiasis. The drug repurposing approach for these 6 kinases show 13 approved drugs and investigational compounds that might play substantial inhibitory roles during combating candidiasis.

Sequence and Structure-Based Analysis of Specificity Determinants in Eukaryotic Protein Kinases

Protein kinases lie at the heart of cell-signaling processes and are often mutated in disease. Kinase target recognition at the active site is in part determined by a few amino acids around the phosphoacceptor residue. However, relatively little is known about how most preferences are encoded in the kinase sequence or how these preferences evolved. Here, we used alignment-based approaches to predict 30 specificity-determining residues (SDRs) for 16 preferences. These were studied with structural models and were validated by activity assays of mutant kinases. Cancer mutation data revealed that kinase SDRs are mutated more frequently than catalytic residues. We have observed that, throughout evolution, kinase specificity has been strongly conserved across orthologs but can diverge after gene duplication, as illustrated by the G protein-coupled receptor kinase family. The identified SDRs can be used to predict kinase specificity from sequence and aid in the interpretation of evolutionary or disease-related genomic variants.

Molecular evidence of hybridization between pig and human Ascaris indicates an interbred species complex infecting humans

Human ascariasis is a major neglected tropical disease caused by the nematode Ascaris lumbricoides. We report a 296 megabase (Mb) reference-quality genome comprised of 17,902 protein-coding genes derived from a single, representative Ascaris worm. An additional 68 worms were collected from 60 human hosts in Kenyan villages where pig husbandry is rare. Notably, the majority of these worms (63/68) possessed mitochondrial genomes that clustered closer to the pig parasite Ascaris suum than to A. lumbricoides. Comparative phylogenomic analyses identified over 11 million nuclear-encoded SNPs but just two distinct genetic types that had recombined across the genomes analyzed. The nuclear genomes had extensive heterozygosity, and all samples existed as genetic mosaics with either A. suum-like or A. lumbricoides-like inheritance patterns supporting a highly interbred Ascaris species genetic complex. As no barriers appear to exist for anthroponotic transmission of these 'hybrid' worms, a one-health approach to control the spread of human ascariasis will be necessary.

High-quality nuclear genome for Sarcoptes scabiei-A critical resource for a neglected parasite

The parasitic mite Sarcoptes scabiei is an economically highly significant parasite of the skin of humans and animals worldwide. In humans, this mite causes a neglected tropical disease (NTD), called scabies. This disease results in major morbidity, disability, stigma and poverty globally and is often associated with secondary bacterial infections. Currently, anti-scabies treatments are not sufficiently effective, resistance to them is emerging and no vaccine is available. Here, we report the first high-quality genome and transcriptomic data for S. scabiei. The genome is 56.6 Mb in size, has a a repeat content of 10.6% and codes for 9,174 proteins. We explored key molecules involved in development, reproduction, host-parasite interactions, immunity and disease. The enhanced 'omic data sets for S. scabiei represent comprehensive and critical resources for genetic, functional genomic, metabolomic, phylogenetic, ecological and/or epidemiological investigations, and will underpin the design and development of new treatments, vaccines and/or diagnostic tests.

Using phosphoproteomics data to understand cellular signaling: a comprehensive guide to bioinformatics resources

Mass spectrometry-based phosphoproteomics is becoming an essential methodology for the study of global cellular signaling. Numerous bioinformatics resources are available to facilitate the translation of phosphopeptide identification and quantification results into novel biological and clinical insights, a critical step in phosphoproteomics data analysis. These resources include knowledge bases of kinases and phosphatases, phosphorylation sites, kinase inhibitors, and sequence variants affecting kinase function, and bioinformatics tools that can predict phosphorylation sites in addition to the kinase that phosphorylates them, infer kinase activity, and predict the effect of mutations on kinase signaling. However, these resources exist in silos and it is challenging to select among multiple resources with similar functions. Therefore, we put together a comprehensive collection of resources related to phosphoproteomics data interpretation, compared the use of tools with similar functions, and assessed the usability from the standpoint of typical biologists or clinicians. Overall, tools could be improved by standardization of enzyme names, flexibility of data input and output format, consistent maintenance, and detailed manuals.

iEKPD 2.0: an update with rich annotations for eukaryotic protein kinases, protein phosphatases and proteins containing phosphoprotein-binding domains

Here, we described the updated database iEKPD 2.0 (http://iekpd.biocuckoo.org) for eukaryotic protein kinases (PKs), protein phosphatases (PPs) and proteins containing phosphoprotein-binding domains (PPBDs), which are key molecules responsible for phosphorylation-dependent signalling networks and participate in the regulation of almost all biological processes and pathways. In total, iEKPD 2.0 contained 197 348 phosphorylation regulators, including 109 912 PKs, 23 294 PPs and 68 748 PPBD-containing proteins in 164 eukaryotic species. In particular, we provided rich annotations for the regulators of eight model organisms, especially humans, by compiling and integrating the knowledge from 100 widely used public databases that cover 13 aspects, including cancer mutations, genetic variations, disease-associated information, mRNA expression, DNA & RNA elements, DNA methylation, molecular interactions, drug-target relations, protein 3D structures, post-translational modifications, protein expressions/proteomics, subcellular localizations and protein functional annotations. Compared with our previously developed EKPD 1.0 (∼0.5 GB), iEKPD 2.0 contains ∼99.8 GB of data with an ∼200-fold increase in data volume. We anticipate that iEKPD 2.0 represents a more useful resource for further study of phosphorylation regulators.

Evolution of protein kinase substrate recognition at the active site

Protein kinases catalyse the phosphorylation of target proteins, controlling most cellular processes. The specificity of serine/threonine kinases is partly determined by interactions with a few residues near the phospho-acceptor residue, forming the so-called kinase-substrate motif. Kinases have been extensively duplicated throughout evolution, but little is known about when in time new target motifs have arisen. Here, we show that sequence variation occurring early in the evolution of kinases is dominated by changes in specificity-determining residues. We then analysed kinase specificity models, based on known target sites, observing that specificity has remained mostly unchanged for recent kinase duplications. Finally, analysis of phosphorylation data from a taxonomically broad set of 48 eukaryotic species indicates that most phosphorylation motifs are broadly distributed in eukaryotes but are not present in prokaryotes. Overall, our results suggest that the set of eukaryotes kinase motifs present today was acquired around the time of the eukaryotic last common ancestor and that early expansions of the protein kinase fold rapidly explored the space of possible target motifs.

Unveiling the Kinomes of Leishmania infantum and L. braziliensis Empowers the Discovery of New Kinase Targets and Antileishmanial Compounds

Leishmaniasis is a neglected tropical disease caused by parasites of the genus Leishmania (NTD) endemic in 98 countries. Although some drugs are available, current treatments deal with issues such as toxicity, low efficacy, and emergence of resistance. Therefore, there is an urgent need to identify new targets for the development of new antileishmanial drugs. Protein kinases (PKs), which play an essential role in many biological processes, have become potential drug targets for many parasitic diseases. A refined bioinformatics pipeline was applied in order to define and compare the kinomes of L. infantum and L. braziliensis, species that cause cutaneous and visceral manifestations of leishmaniasis in the Americas, the latter being potentially fatal if untreated. Respectively, 224 and 221 PKs were identified in L. infantum and L. braziliensis overall. Almost all unclassified eukaryotic PKs were assigned to six of nine major kinase groups and, consequently, most have been classified into family and subfamily. Furthermore, revealing the kinomes for both Leishmania species allowed for the prioritization of potential drug targets that could be explored for discovering new drugs against leishmaniasis. Finally, we used a drug repurposing approach and prioritized seven approved drugs and investigational compounds to be experimentally tested against Leishmania. Trametinib and NMS-1286937 inhibited the growth of L. infantum and L. braziliensis promastigotes and amastigotes and therefore might be good candidates for the drug repurposing pipeline.

Genomic insights into multidrug-resistance, mating and virulence in Candida auris and related emerging species

Candida auris is an emergent multidrug-resistant fungal pathogen causing increasing reports of outbreaks. While distantly related to C. albicans and C. glabrata, C. auris is closely related to rarely observed and often multidrug-resistant species from the C. haemulonii clade. Here, we analyze near complete genome assemblies for the four C. auris clades and three related species, and map intra- and inter-species rearrangements across the seven chromosomes. Using RNA-Seq-guided gene predictions, we find that most mating and meiosis genes are conserved and that clades contain either the MTLa or MTLα mating loci. Comparing the genomes of these emerging species to those of other Candida species identifies genes linked to drug resistance and virulence, including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11. Gene expression analysis identifies transporters and metabolic regulators specific to C. auris and those conserved with related species which may contribute to differences in drug response in this emerging fungal clade.

Candida auris is an emergent fungal pathogen that is resistant to multiple antifungals. Here, Muñoz et al. analyse genomic sequences for isolates from each of the four major C. auris clades and for three related species, and identify genetic features associated with virulence, antifungal resistance and mating.

Computationally-guided drug repurposing enables the discovery of kinase targets and inhibitors as new schistosomicidal agents

The development of novel therapeutics is urgently required for diseases where existing treatments are failing due to the emergence of resistance. This is particularly pertinent for parasitic infections of the tropics and sub-tropics, referred to collectively as neglected tropical diseases, where the commercial incentives to develop new drugs are weak. One such disease is schistosomiasis, a highly prevalent acute and chronic condition caused by a parasitic helminth infection, with three species of the genus Schistosoma infecting humans. Currently, a single 40-year old drug, praziquantel, is available to treat all infective species, but its use in mass drug administration is leading to signs of drug-resistance emerging. To meet the challenge of developing new therapeutics against this disease, we developed an innovative computational drug repurposing pipeline supported by phenotypic screening. The approach highlighted several protein kinases as interesting new biological targets for schistosomiasis as they play an essential role in many parasite’s biological processes. Focusing on this target class, we also report the first elucidation of the kinome of Schistosoma japonicum, as well as updated kinomes of S. mansoni and S. haematobium. In comparison with the human kinome, we explored these kinomes to identify potential targets of existing inhibitors which are unique to Schistosoma species, allowing us to identify novel targets and suggest approved drugs that might inhibit them. These include previously suggested schistosomicidal agents such as bosutinib, dasatinib, and imatinib as well as new inhibitors such as vandetanib, saracatinib, tideglusib, alvocidib, dinaciclib, and 22 newly identified targets such as CHK1, CDC2, WEE, PAKA, MEK1. Additionally, the primary and secondary targets in Schistosoma of those approved drugs are also suggested, allowing for the development of novel therapeutics against this important yet neglected disease.

The rise of resistance through the intensive use of drugs targeted to treat specific infectious diseases means that new therapeutics are continually required. Diseases common in the tropics and sub-tropics, classified as neglected tropical diseases, suffer from a lack of new drug treatments due to the difficulty in developing new drugs and the lack of market incentive. One such disease is schistosomiasis, a major human helminth disease caused by worms from the genus Schistosoma. It is currently treated by a 40-year old drug, praziquantel, but its widespread use has led to signs of drug-resistance emerging, with no alternative effective treatments available. To meet this challenge, we have developed an innovative computational drug repurposing pipeline supported by experimental phenotypic screening. Protein kinases emerged from our pipeline as interesting new biological targets. Given that many human kinase inhibitors have been successfully applied specially in cancer therapy and kinases have conserved structures and functions, we also undertook a detailed analysis of the kinases present in all infective Schistosoma species and human host. This allowed identification of new Schistosoma-specific kinase targets and suggest approved drugs to be used for treating schistosomiasis as well as opening new avenues to treat this neglected disease.

Kinome Expansion in the Fusarium oxysporum Species Complex Driven by Accessory Chromosomes

Isolates of Fusarium oxysporum are adapted to survive a wide range of host and nonhost conditions. In addition, F. oxysporum was recently recognized as the top emerging opportunistic fungal pathogen infecting immunocompromised humans. The sensory and response networks of these fungi undoubtedly play a fundamental role in establishing the adaptability of this group. We have examined the kinomes of 12 F. oxysporum isolates and highlighted kinase families that distinguish F. oxysporum from other fungi, as well as different isolates from one another. The amplification of kinases involved in environmental signal relay and regulating downstream cellular responses clearly sets Fusarium apart from other Ascomycetes. Although the functions of many of these kinases are still unclear, their specific proliferation highlights them as a result of the evolutionary forces that have shaped this species complex and clearly marks them as targets for exploitation in order to combat disease.

The Fusarium oxysporum species complex (FOSC) is a group of soilborne pathogens causing severe disease in more than 100 plant hosts, while individual strains exhibit strong host specificity. Both chromosome transfer and comparative genomics experiments have demonstrated that lineage-specific (LS) chromosomes contribute to the host-specific pathogenicity. However, little is known about the functional importance of genes encoded in these LS chromosomes. Focusing on signaling transduction, this study compared the kinomes of 12 F. oxysporum isolates, including both plant and human pathogens and 1 nonpathogenic biocontrol strain, with 7 additional publicly available ascomycete genomes. Overall, F. oxysporum kinomes are the largest, facilitated in part by the acquisitions of the LS chromosomes. The comparative study identified 99 kinases that are present in almost all examined fungal genomes, forming the core signaling network of ascomycete fungi. Compared to the conserved ascomycete kinome, the expansion of the F. oxysporum kinome occurs in several kinase families such as histidine kinases that are involved in environmental signal sensing and target of rapamycin (TOR) kinase that mediates cellular responses. Comparative kinome analysis suggests a convergent evolution that shapes individual F. oxysporum isolates with an enhanced and unique capacity for environmental perception and associated downstream responses.

IMPORTANCE Isolates of Fusarium oxysporum are adapted to survive a wide range of host and nonhost conditions. In addition, F. oxysporum was recently recognized as the top emerging opportunistic fungal pathogen infecting immunocompromised humans. The sensory and response networks of these fungi undoubtedly play a fundamental role in establishing the adaptability of this group. We have examined the kinomes of 12 F. oxysporum isolates and highlighted kinase families that distinguish F. oxysporum from other fungi, as well as different isolates from one another. The amplification of kinases involved in environmental signal relay and regulating downstream cellular responses clearly sets Fusarium apart from other Ascomycetes. Although the functions of many of these kinases are still unclear, their specific proliferation highlights them as a result of the evolutionary forces that have shaped this species complex and clearly marks them as targets for exploitation in order to combat disease.

Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups

Despite the substantial amount of genomic and transcriptomic data available for a wide range of eukaryotic organisms, most genomes are still in a draft state and can have inaccurate gene predictions. To gain a sound understanding of the biology of an organism, it is crucial that inferred protein sequences are accurately identified and annotated. However, this can be challenging to achieve, particularly for organisms such as parasitic worms (helminths), as most gene prediction approaches do not account for substantial phylogenetic divergence from model organisms, such as Caenorhabditis elegans and Drosophila melanogaster, whose genomes are well-curated. In this paper, we describe a bioinformatic strategy for the curation of gene families and subsequent annotation of encoded proteins. This strategy relies on pairwise gene curation between at least two closely related species using genomic and transcriptomic data sets, and is built on recent work on kinase complements of parasitic worms. Here, we discuss salient technical aspects of this strategy and its implications for the curation of protein families more generally.

Advances in kinome research of parasitic worms - implications for fundamental research and applied biotechnological outcomes

Protein kinases are enzymes that play essential roles in the regulation of many cellular processes. Despite expansions in the fields of genomics, transcriptomics and bioinformatics, there is limited information on the kinase complements (kinomes) of most eukaryotic organisms, including parasitic worms that cause serious diseases of humans and animals. The biological uniqueness of these worms and the draft status of their genomes pose challenges for the identification and classification of protein kinases using established tools. In this article, we provide an account of kinase biology, the roles of kinases in diseases and their importance as drug targets, and drug discovery efforts in key socioeconomically important parasitic worms. In this context, we summarise methods and resources commonly used for the curation, identification, classification and functional annotation of protein kinase sequences from draft genomes; review recent advances made in the characterisation of the worm kinomes; and discuss the implications of these advances for investigating kinase signalling and developing small-molecule inhibitors as new anti-parasitic drugs.

Genome analysis reveals evolutionary mechanisms of adaptation in systemic dimorphic fungi

Dimorphic fungal pathogens cause a significant human disease burden and unlike most fungal pathogens affect immunocompetent hosts. To examine the origin of virulence of these fungal pathogens, we compared genomes of classic systemic, opportunistic, and non-pathogenic species, including Emmonsia and two basal branching, non-pathogenic species in the Ajellomycetaceae, Helicocarpus griseus and Polytolypa hystricis. We found that gene families related to plant degradation, secondary metabolites synthesis, and amino acid and lipid metabolism are retained in H. griseus and P. hystricis. While genes involved in the virulence of dimorphic pathogenic fungi are conserved in saprophytes, changes in the copy number of proteases, kinases and transcription factors in systemic dimorphic relative to non-dimorphic species may have aided the evolution of specialized gene regulatory programs to rapidly adapt to higher temperatures and new nutritional environments. Notably, both of the basal branching, non-pathogenic species appear homothallic, with both mating type locus idiomorphs fused at a single locus, whereas all related pathogenic species are heterothallic. These differences revealed that independent changes in nutrient acquisition capacity have occurred in the Onygenaceae and Ajellomycetaceae, and underlie how the dimorphic pathogens have adapted to the human host and decreased their capacity for growth in environmental niches.

Basidiomycete-specific PsCaMKL1 encoding a CaMK-like protein kinase is required for full virulence of Puccinia striiformis f. sp. tritici

Calcium/calmodulin-dependent kinases (CaMKs) are Ser/Thr protein kinases (PKs) that respond to changes in cytosolic free Ca²⁺ and play diverse roles in eukaryotes. In fungi, CAMKs are generally classified into four families CAMK1, CAMKL, RAD53 and CAMK-Unique. Among these, CAMKL constitutes the largest family. In some fungal plant pathogens, members of the CaMKL family have been shown to be responsible for pathogenesis. However, little is known about their role(s) in rust fungi. In this study, we functionally characterized a novel PK gene, PsCaMKL1, from Puccinia striiformis f. sp. tritici (Pst). PsCaMKL1 belongs to a group of PKs that is evolutionarily specific to basidiomyceteous fungi. PsCaMKL1 shows little intra-species polymorphism between Pst isolates. PsCaMKL1 transcripts are highly elevated at early infection stages, whereas gene expression is downregulated in barely germinated urediospores under KN93 treatment. Overexpression of PsCaMKL1 in fission yeast increased resistance to environmental stresses. Knock down of PsCaMKL1 using host-induced gene silencing (HIGS) reduced the virulence of Pst accompanied by reactive oxygen species (ROS) accumulation and a hypersensitive response. These results suggest that PsCaMKL1 is a novel pathogenicity factor that exerts it virulence function by regulating ROS production in wheat.

Whipworm kinomes reflect a unique biology and adaptation to the host animal

Roundworms belong to a diverse phylum (Nematoda) which is comprised of many parasitic species including whipworms (genus Trichuris). These worms have adapted to a biological niche within the host and exhibit unique morphological characteristics compared with other nematodes. Although these adaptations are known, the underlying molecular mechanisms remain elusive. The availability of genomes and transcriptomes of some whipworms now enables detailed studies of their molecular biology. Here, we defined and curated the full complement of an important class of enzymes, the protein kinases (kinomes) of two species of Trichuris, using an advanced and integrated bioinformatic pipeline. We investigated the transcription of Trichuris suis kinase genes across developmental stages, sexes and tissues, and reveal that selectively transcribed genes can be linked to central roles in developmental and reproductive processes. We also classified and functionally annotated the curated kinomes by integrating evidence from structural modelling and pathway analyses, and compared them with other curated kinomes of phylogenetically diverse nematode species. Our findings suggest unique adaptations in signalling processes governing worm morphology and biology, and provide an important resource that should facilitate experimental investigations of kinases and the biology of signalling pathways in nematodes.

Novel taxa of thermally dimorphic systemic pathogens in the Ajellomycetaceae (Onygenales)

Recent discoveries of novel systemic fungal pathogens with thermally dimorphic yeast-like phases have challenged the current taxonomy of the Ajellomycetaceae, a family currently comprising the genera Blastomyces, Emmonsia, Emmonsiellopsis, Helicocarpus, Histoplasma, Lacazia and Paracoccidioides. Our morphological, phylogenetic and phylogenomic analyses demonstrated species relationships and their specific phenotypes, clarified generic boundaries and provided the first annotated genome assemblies to support the description of two new species. A new genus, Emergomyces, accommodates Emmonsia pasteuriana as type species, and the new species Emergomyces africanus, the aetiological agent of case series of disseminated infections in South Africa. Both species produce small yeast cells that bud at a narrow base at 37°C and lack adiaspores, classically associated with the genus Emmonsia. Another novel dimorphic pathogen, producing broad-based budding cells at 37°C and occurring outside North America, proved to belong to the genus Blastomyces, and is described as Blastomyces percursus.

Analyses of Compact Trichinella Kinomes Reveal a MOS-Like Protein Kinase with a Unique N-Terminal Domain

Parasitic worms of the genus Trichinella (phylum Nematoda; class Enoplea) represent a complex of at least twelve taxa that infect a range of different host animals, including humans, around the world. They are foodborne, intracellular nematodes, and their life cycles differ substantially from those of other nematodes. The recent characterization of the genomes and transcriptomes of all twelve recognized taxa of Trichinella now allows, for the first time, detailed studies of their molecular biology. In the present study, we defined, curated, and compared the protein kinase complements (kinomes) of Trichinella spiralis and T. pseudospiralis using an integrated bioinformatic workflow employing transcriptomic and genomic data sets. We examined how variation in the kinome might link to unique aspects of Trichinella morphology, biology, and evolution. Furthermore, we utilized in silico structural modeling to discover and characterize a novel, MOS-like kinase with an unusual, previously undescribed N-terminal domain. Taken together, the present findings provide a basis for comparative investigations of nematode kinomes, and might facilitate the identification of Enoplea-specific intervention and diagnostic targets. Importantly, the in silico modeling approach assessed here provides an exciting prospect of being able to identify and classify currently unknown (orphan) kinases, as a foundation for their subsequent structural and functional investigation.

De Novo Analysis of the Transcriptome of Meloidogyne enterolobii to Uncover Potential Target Genes for Biological Control

Meloidogyne enterolobii is one of the obligate biotrophic root-knot nematodes that has the ability to reproduce on many economically-important crops. We carried out de novo sequencing of the transcriptome of M. enterolobii using Roche GS FLX and obtained 408,663 good quality reads that were assembled into 8193 contigs and 31,860 singletons. We compared the transcripts in different nematodes that were potential targets for biological control. These included the transcripts that putatively coded for CAZymes, kinases, neuropeptide genes and secretory proteins and those that were involved in the RNAi pathway and immune signaling. Typically, 75 non-membrane secretory proteins with signal peptides secreted from esophageal gland cells were identified as putative effectors, three of which were preliminarily examined using a PVX (pGR107)-based high-throughput transient plant expression system in Nicotiana benthamiana (N. benthamiana). Results showed that these candidate proteins suppressed the programmed cell death (PCD) triggered by the pro-apoptosis protein BAX, and one protein also caused necrosis, suggesting that they might suppress plant immune responses to promote pathogenicity. In conclusion, the current study provides comprehensive insight into the transcriptome of M. enterolobii for the first time and lays a foundation for further investigation and biological control strategies.

The Haemonchus contortus kinome--a resource for fundamental molecular investigations and drug discovery

BACKGROUND

Protein kinases regulate a plethora of essential signalling and other biological pathways in all eukaryotic organisms, but very little is known about them in most parasitic nematodes.

METHODS

Here, we defined, for the first time, the entire complement of protein kinases (kinome) encoded in the barber's pole worm (Haemonchus contortus) through an integrated analysis of transcriptomic and genomic datasets using an advanced bioinformatic workflow.

RESULTS

We identified, curated and classified 432 kinases representing ten groups, 103 distinct families and 98 subfamilies. A comparison of the kinomes of H. contortus and Caenorhabditis elegans (a related, free-living nematode) revealed considerable variation in the numbers of casein kinases, tyrosine kinases and Ca(2+)/calmodulin-dependent protein kinases, which likely relate to differences in biology, habitat and life cycle between these worms. Moreover, a suite of kinase genes was selectively transcribed in particular developmental stages of H. contortus, indicating central roles in developmental and reproductive processes. In addition, using a ranking system, drug targets (n = 13) and associated small-molecule effectors (n = 1517) were inferred.

CONCLUSIONS

The H. contortus kinome will provide a useful resource for fundamental investigations of kinases and signalling pathways in this nematode, and should assist future anthelmintic discovery efforts; this is particularly important, given current drug resistance problems in parasitic nematodes.

Defining the Schistosoma haematobium kinome enables the prediction of essential kinases as anti-schistosome drug targets

The blood fluke Schistosoma haematobium causes urogenital schistosomiasis, a neglected tropical disease (NTD) that affects more than 110 million people. Treating this disease by targeted or mass administration with a single chemical, praziquantel, carries the risk that drug resistance will develop in this pathogen. Therefore, there is an imperative to search for new drug targets in S. haematobium and other schistosomes. In this regard, protein kinases have potential, given their essential roles in biological processes and as targets for drugs already approved by the US Food and Drug Administration (FDA) for use in humans. In this context, we defined here the kinome of S. haematobium using a refined bioinformatic pipeline. We classified, curated and annotated predicted kinases, and assessed the developmental transcription profiles of kinase genes. Then, we prioritised a panel of kinases as potential drug targets and inferred chemicals that bind to them using an integrated bioinformatic pipeline. Most kinases of S. haematobium are very similar to those of its congener, S. mansoni, offering the prospect of designing chemicals that kill both species. Overall, this study provides a global insight into the kinome of S. haematobium and should assist the repurposing or discovery of drugs against schistosomiasis.

Genome sequence of Valsa canker pathogens uncovers a potential adaptation of colonization of woody bark

Canker caused by ascomycetous Valsa species are among the most destructive diseases of woody plants worldwide. These pathogens are distinct from other pathogens because they only effectively attack tree bark in the field. To unravel the potential adaptation mechanism of bark colonization, we examined the genomes of Valsa mali and Valsa pyri that preferentially infect apple and pear, respectively. We reported the 44.7 and 35.7 Mb genomes of V. mali and V. pyri, respectively. We also identified the potential genomic determinants of wood colonization by comparing them with related cereal pathogens. Both genomes encode a plethora of pathogenicity-related genes involved in plant cell wall degradation and secondary metabolite biosynthesis. In order to adapt to the nutrient limitation and low pH environment in bark, they seem to employ membrane transporters associated with nitrogen uptake and secrete proteases predominantly with acidic pH optima. Remarkably, both Valsa genomes are especially suited for pectin decomposition, but are limited in lignocellulose and cutin degradation. Besides many similarities, the two genomes show distinct variations in many secondary metabolism gene clusters. Our results show a potential adaptation of Valsa canker pathogens to colonize woody bark. Secondary metabolism gene clusters are probably responsible for this host specificity.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Dimorphic fungal pathogens including Blastomyces are the cause of major fungal diseases in North and South America. The genus Emmonsia includes species infecting small mammals as well as a newly emerging pathogenic species recently reported in HIV-positive patients in South Africa. Here, we synthesize both genome sequencing of four isolates of Blastomyces and two species of Emmonsia as well as deep sequencing of Blastomyces RNA to draw major new insights into the evolution of this group and the pathogen response to infection. We investigate the trajectory of genome evolution of this group, characterizing the phylogenetic relationships of these species, a remarkable genome expansion that formed large isochore-like regions of low GC content in Blastomyces, and variation of gene content, related to host interaction, among the dimorphic fungal pathogens. Using RNA-Seq, we profile the response of Blastomyces to macrophage and mouse pulmonary infection, identifying key pathways and novel virulence factors. The identification of key fungal genes involved in adaptation to the host suggests targets for further study and therapeutic intervention in Blastomyces and related dimorphic fungal pathogens.

The Arabidopsis Kinome: phylogeny and evolutionary insights into functional diversification

Background

Protein kinases constitute a particularly large protein family in Arabidopsis with important functions in cellular signal transduction networks. At the same time Arabidopsis is a model plant with high frequencies of gene duplications. Here, we have conducted a systematic analysis of the Arabidopsis kinase complement, the kinome, with particular focus on gene duplication events. We matched Arabidopsis proteins to a Hidden-Markov Model of eukaryotic kinases and computed a phylogeny of 942 Arabidopsis protein kinase domains and mapped their origin by gene duplication.

Results

The phylogeny showed two major clades of receptor kinases and soluble kinases, each of which was divided into functional subclades. Based on this phylogeny, association of yet uncharacterized kinases to families was possible which extended functional annotation of unknowns. Classification of gene duplications within these protein kinases revealed that representatives of cytosolic subfamilies showed a tendency to maintain segmentally duplicated genes, while some subfamilies of the receptor kinases were enriched for tandem duplicates. Although functional diversification is observed throughout most subfamilies, some instances of functional conservation among genes transposed from the same ancestor were observed. In general, a significant enrichment of essential genes was found among genes encoding for protein kinases.

Conclusions

The inferred phylogeny allowed classification and annotation of yet uncharacterized kinases. The prediction and analysis of syntenic blocks and duplication events within gene families of interest can be used to link functional biology to insights from an evolutionary viewpoint. The approach undertaken here can be applied to any gene family in any organism with an annotated genome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-548) contains supplementary material, which is available to authorized users.

EKPD: a hierarchical database of eukaryotic protein kinases and protein phosphatases

We present here EKPD (http://ekpd.biocuckoo.org), a hierarchical database of eukaryotic protein kinases (PKs) and protein phosphatases (PPs), the key molecules responsible for the reversible phosphorylation of proteins that are involved in almost all aspects of biological processes. As extensive experimental and computational efforts have been carried out to identify PKs and PPs, an integrative resource with detailed classification and annotation information would be of great value for both experimentalists and computational biologists. In this work, we first collected 1855 PKs and 347 PPs from the scientific literature and various public databases. Based on previously established rationales, we classified all of the known PKs and PPs into a hierarchical structure with three levels, i.e. group, family and individual PK/PP. There are 10 groups with 149 families for the PKs and 10 groups with 33 families for the PPs. We constructed 139 and 27 Hidden Markov Model profiles for PK and PP families, respectively. Then we systematically characterized ∼50,000 PKs and >10,000 PPs in eukaryotes. In addition, >500 PKs and >400 PPs were computationally identified by ortholog search. Finally, the online service of the EKPD database was implemented in PHP + MySQL + JavaScript.