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.     

Signal transduction protein array analysis links LRRK2 to Ste20 kinases and PKC zeta that modulate neuronal plasticity

Background

Dominant mutations in leucine-rich repeat kinase 2 (LRRK2) are the most common genetic cause of Parkinson''s disease, however, the underlying pathogenic mechanisms are poorly understood. Several in vitro studies have shown that the most frequent mutation, LRRK2(G2019S), increases kinase activity and impairs neuronal survival. LRRK2 has been linked to the mitogen-activated protein kinase kinase kinase family and the receptor-interacting protein kinases based on sequence similarity within the kinase domain and in vitro substrate phosphorylation.

Methodology/Principal Findings

We used an unbiased proteomic approach to identify the kinase signaling pathways wherein LRRK2 may be active. By incubation of protein microarrays containing 260 signal transduction proteins we detected four arrayed Ste20 serine/threonine kinase family members (TAOK3, STK3, STK24, STK25) as novel LRRK2 substrates and LRRK2 interacting proteins, respectively. Moreover, we found that protein kinase C (PKC) zeta binds and phosphorylates LRRK2 both in vitro and in vivo.

Conclusions/Significance

Ste20 kinases and PKC zeta contribute to neuronal Tau phosphorylation, neurite outgrowth and synaptic plasticity under physiological conditions. Our data suggest that these kinases may also be involved in synaptic dysfunction and neurite fragmentation in transgenic mice and in human PD patients carrying toxic gain-of-function LRRK2 mutations.     

Kinome analysis of receptor-induced phosphorylation in human natural killer cells

Background

Natural killer (NK) cells contribute to the defense against infected and transformed cells through the engagement of multiple germline-encoded activation receptors. Stimulation of the Fc receptor CD16 alone is sufficient for NK cell activation, whereas other receptors, such as 2B4 (CD244) and DNAM-1 (CD226), act synergistically. After receptor engagement, protein kinases play a major role in signaling networks controlling NK cell effector functions. However, it has not been characterized systematically which of all kinases encoded by the human genome (kinome) are involved in NK cell activation.

Results

A kinase-selective phosphoproteome approach enabled the determination of 188 kinases expressed in human NK cells. Crosslinking of CD16 as well as 2B4 and DNAM-1 revealed a total of 313 distinct kinase phosphorylation sites on 109 different kinases. Phosphorylation sites on 21 kinases were similarly regulated after engagement of either CD16 or co-engagement of 2B4 and DNAM-1. Among those, increased phosphorylation of FYN, KCC2G (CAMK2), FES, and AAK1, as well as the reduced phosphorylation of MARK2, were reproducibly observed both after engagement of CD16 and co-engagement of 2B4 and DNAM-1. Notably, only one phosphorylation on PAK4 was differentally regulated.

Conclusions

The present study has identified a significant portion of the NK cell kinome and defined novel phosphorylation sites in primary lymphocytes. Regulated phosphorylations observed in the early phase of NK cell activation imply these kinases are involved in NK cell signaling. Taken together, this study suggests a largely shared signaling pathway downstream of distinct activation receptors and constitutes a valuable resource for further elucidating the regulation of NK cell effector responses.     

ProKinO: an ontology for integrative analysis of protein kinases in cancer

Background

Protein kinases are a large and diverse family of enzymes that are genomically altered in many human cancers. Targeted cancer genome sequencing efforts have unveiled the mutational profiles of protein kinase genes from many different cancer types. While mutational data on protein kinases is currently catalogued in various databases, integration of mutation data with other forms of data on protein kinases such as sequence, structure, function and pathway is necessary to identify and characterize key cancer causing mutations. Integrative analysis of protein kinase data, however, is a challenge because of the disparate nature of protein kinase data sources and data formats.

Results

Here, we describe ProKinO, a protein kinase-specific ontology, which provides a controlled vocabulary of terms, their hierarchy, and relationships unifying sequence, structure, function, mutation and pathway information on protein kinases. The conceptual representation of such diverse forms of information in one place not only allows rapid discovery of significant information related to a specific protein kinase, but also enables large-scale integrative analysis of protein kinase data in ways not possible through other kinase-specific resources. We have performed several integrative analyses of ProKinO data and, as an example, found that a large number of somatic mutations (∼288 distinct mutations) associated with the haematopoietic neoplasm cancer type map to only 8 kinases in the human kinome. This is in contrast to glioma, where the mutations are spread over 82 distinct kinases. We also provide examples of how ontology-based data analysis can be used to generate testable hypotheses regarding cancer mutations.

Conclusion

We present an integrated framework for large-scale integrative analysis of protein kinase data. Navigation and analysis of ontology data can be performed using the ontology browser available at: http://vulcan.cs.uga.edu/prokino.     

Mutational Profiling of Kinases in Human Tumours of Pancreatic Origin Identifies Candidate Cancer Genes in Ductal and Ampulla of Vater Carcinomas

Background

Protein kinases are key regulators of cellular processes (such as proliferation, apoptosis and invasion) that are often deregulated in human cancers. Accordingly, kinase genes have been the first to be systematically analyzed in human tumors leading to the discovery that many oncogenes correspond to mutated kinases. In most cases the genetic alterations translate in constitutively active kinase proteins, which are amenable of therapeutic targeting. Tumours of the pancreas are aggressive neoplasms for which no effective therapeutic strategy is currently available.

Methodology/Principal Findings

We conducted a DNA-sequence analysis of a selected set of 35 kinase genes in a panel of 52 pancreatic exocrine neoplasms, including 36 pancreatic ductal adenocarcinoma, and 16 ampulla of Vater cancer. Among other changes we found somatic mutations in ATM, EGFR, EPHA3, EPHB2, and KIT, none of which was previously described in cancers.

Conclusions/Significance

Although the alterations identified require further experimental evaluation, the localization within defined protein domains indicates functional relevance for most of them. Some of the mutated genes, including the tyrosine kinases EPHA3 and EPHB2, are clearly amenable to pharmacological intervention and could represent novel therapeutic targets for these incurable cancers.     

Genome-Wide Analysis of the Phosphoinositide Kinome from Two Ciliates Reveals Novel Evolutionary Links for Phosphoinositide Kinases in Eukaryotic Cells

Background

The complexity of phosphoinositide signaling in higher eukaryotes is partly due to expansion of specific families and types of phosphoinositide kinases (PIKs) that can generate all phosphoinositides via multiple routes. This is particularly evident in the PI3Ks and PIPKs, and it is considered an evolutionary trait associated with metazoan diversification. Yet, there are limited comprehensive studies on the PIK repertoire of free living unicellular organisms.

Methodology/Principal Findings

We undertook a genome-wide analysis of putative PIK genes in two free living ciliated cells, Tetrahymena and Paramecium. The Tetrahymena thermophila and Paramecium tetraurelia genomes were probed with representative kinases from all families and types. Putative homologs were verified by EST, microarray and deep RNA sequencing database searches and further characterized for domain structure, catalytic efficiency, expression patterns and phylogenetic relationships. In total, we identified and characterized 22 genes in the Tetrahymena thermophila genome and 62 highly homologues genes in Paramecium tetraurelia suggesting a tight evolutionary conservation in the ciliate lineage. Comparison to the kinome of fungi reveals a significant expansion of PIK genes in ciliates.

Conclusions/Significance

Our study highlights four important aspects concerning ciliate and other unicellular PIKs. First, ciliate-specific expansion of PI4KIII-like genes. Second, presence of class I PI3Ks which, at least in Tetrahymena, are associated with a metazoan-type machinery for PIP₃ signaling. Third, expansion of divergent PIPK enzymes such as the recently described type IV transmembrane PIPKs. Fourth, presence of possible type II PIPKs and presumably inactive PIKs (hence, pseudo-PIKs) not previously described. Taken together, our results provide a solid framework for future investigation of the roles of PIKs in ciliates and indicate that novel functions and novel regulatory pathways of phosphoinositides may be more widespread than previously thought in unicellular organisms.     

Novel and unique domains in aminoacyl-tRNA synthetases from human fungal pathogens Aspergillus niger,Candida albicans and Cryptococcus neoformans

Background

Some species of fungi can cause serious human diseases, particularly to immuno-compromised individuals. Opportunistic fungal infections are a leading cause of mortality, and present an emerging challenge that requires development of new and effective therapeutics. Aminoacyl-tRNA synthetases (aaRSs) are indispensable components of cellular protein translation machinery and can be targeted for discovery of novel anti-fungal agents.

Results

Validation of aaRSs as potential drug targets in pathogenic microbes prompted us to investigate the genomic distribution of aaRSs within three fungi that infect humans – A. niger, C. albicans and C. neoformans. Hidden Markov Models were built for aaRSs and related proteins to search for homologues in these fungal genomes. Here, we provide a detailed and comprehensive annotation for 3 fungal genome aaRSs and their associated proteins. We delineate predicted localizations, subdomain architectures and prevalence of unusual motifs within these aaRSs. Several fungal aaRSs have unique domain appendages of unknown function e.g. A. niger AsxRS and C. neoformans TyrRS have additional domains that are absent from human homologs.

Conclusions

Detailed comparisons of fungal aaRSs with human homologs suggest key differences that could be exploited for specific drug targeting. Our cataloging and structural analyses provide a comprehensive foundation for experimentally dissecting fungal aaRSs that may enable development of new anti-fungal agents.

Electronic supplementary material

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

Findings from the Peutz-Jeghers Syndrome Registry of Uruguay

Background

Peutz-Jeghers syndrome (PJS) is characterized by intestinal polyposis, mucocutaneous pigmentation and an increased cancer risk, usually caused by mutations of the STK11 gene. This study collected epidemiological, clinical and genetic data from all Uruguayan PJS patients.

Methods

Clinical data were obtained from public and private medical centers and updated annually. Sequencing of the STK11 gene in one member of each family was performed.

Results and discussion

25 cases in 11 unrelated families were registered (15 males, 10 females). The average age of diagnosis and death was 18 and 41 years respectively. All patients had characteristic PJS pigmentation and gastrointestinal polyps. 72% required urgent surgery due to intestinal obstruction. 3 families had multiple cases of seizure disorder, representing 20% of cases. 28% developed cancer and two patients had more than one cancer. An STK11 mutation was found in 8 of the 9 families analyzed. A unique M136K missense mutation was noted in one family. Comparing annual live births and PJS birth records from 1970 to 2009 yielded an incidence of 1 in 155,000.

Conclusion

The Uruguayan Registry for Peutz-Jeghers patients showed a high chance of emergent surgery, epilepsy, cancer and shortened life expectancy. The M136K missense mutation is a newly reported STK 11 mutation.     

Classification of Protein Kinases on the Basis of Both Kinase and Non-Kinase Regions

Background

Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multi-domain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies.

Methodology

Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica.

Conclusions

The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multi-domain architecture.     

Characterization of the MLO gene family in Rosaceae and gene expression analysis in Malus domestica

Background

Powdery mildew (PM) is a major fungal disease of thousands of plant species, including many cultivated Rosaceae. PM pathogenesis is associated with up-regulation of MLO genes during early stages of infection, causing down-regulation of plant defense pathways. Specific members of the MLO gene family act as PM-susceptibility genes, as their loss-of-function mutations grant durable and broad-spectrum resistance.

Results

We carried out a genome-wide characterization of the MLO gene family in apple, peach and strawberry, and we isolated apricot MLO homologs through a PCR-approach. Evolutionary relationships between MLO homologs were studied and syntenic blocks constructed. Homologs that are candidates for being PM susceptibility genes were inferred by phylogenetic relationships with functionally characterized MLO genes and, in apple, by monitoring their expression following inoculation with the PM causal pathogen Podosphaera leucotricha.

Conclusions

Genomic tools available for Rosaceae were exploited in order to characterize the MLO gene family. Candidate MLO susceptibility genes were identified. In follow-up studies it can be investigated whether silencing or a loss-of-function mutations in one or more of these candidate genes leads to PM resistance.

Electronic supplementary material

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

Rana grylio Virus (RGV) 50L Is Associated with Viral Matrix and Exhibited Two Distribution Patterns

Background

The complete genome of Rana grylio virus (RGV) was sequenced and analyzed recently, which revealed that RGV 50L had homologues in many iridoviruses with different identities; however, the characteristics and functions of 50L have not been studied yet.

Methodology/Principal Findings

We cloned and characterized RGV50L, and revealed 50L functions in virus assembly and gene regulation. 50L encoded a 499-amino acid structural protein of about 85 kDa in molecular weight and contained a nuclear localization signal (NLS) and a helix- extension-helix motif. Drug inhibition assay demonstrated that 50L was an immediate-early (IE) gene. Immuno-fluorescence assay revealed that 50L appeared early and persisted in RGV-infected cells following two distribution patterns. One pattern was that 50L exhibited a cytoplasm-nucleus- viromatrix distribution pattern, and mutagenesis of the NLS motif revealed that localization of 50L in the nucleus was NLS-dependent; the other was that 50L co-localized with viral matrix which plays important roles in virus assembly and the life circle of viruses.

Conclusions/Significance

RGV 50L is a novel iridovirus IE gene encoded structural protein which plays important roles in virus assembly.     

Rust secreted protein Ps87 is conserved in diverse fungal pathogens and contains a RXLR-like motif sufficient for translocation into plant cells

Background

Effector proteins of biotrophic plant pathogenic fungi and oomycetes are delivered into host cells and play important roles in both disease development and disease resistance response. How obligate fungal pathogen effectors enter host cells is poorly understood. The Ps87 gene of Puccinia striiformis encodes a protein that is conserved in diverse fungal pathogens. Ps87 homologs from a clade containing rust fungi are predicted to be secreted. The aim of this study is to test whether Ps87 may act as an effector during Puccinia striiformis infection.

Methodology/Principal Findings

Yeast signal sequence trap assay showed that the rust protein Ps87 could be secreted from yeast cells, but a homolog from Magnaporthe oryzae that was not predicted to be secreted, could not. Cell re-entry and protein uptake assays showed that a region of Ps87 containing a conserved RXLR-like motif [K/R]RLTG was confirmed to be capable of delivering oomycete effector Avr1b into soybean leaf cells and carrying GFP into soybean root cells. Mutations in the Ps87 motif (KRLTG) abolished the protein translocation ability.

Conclusions/Significance

The results suggest that Ps87 and its secreted homologs could utilize similar protein translocation machinery as those of oomycete and other fungal pathogens. Ps87 did not show direct suppression activity on plant defense responses. These results suggest Ps87 may represent an “emerging effector” that has recently acquired the ability to enter plant cells but has not yet acquired the ability to alter host physiology.     

The Repertoire of Heterotrimeric G Proteins and RGS Proteins in Ciona intestinalis

Background

Heterotrimeric G proteins and regulators of G protein signaling (RGS) proteins are key downstream interacting partners in the G protein coupled receptor (GPCR) signaling pathway. The highly versatile GPCR transmembrane signaling system is a consequence of the coupling of a diverse set of receptors to downstream partners that include multiple subforms of G proteins and regulatory proteins including RGS proteins, among others. While the GPCR repertoire of Ciona intestinalis, representing the basal chordate is known, the repertoire of the heterotrimeric G proteins and RGS proteins is unknown.

Methodology/Principal Findings

In the present study, we performed an in-silico genome-wide search of C. intestinalis for its complement of G proteins and RGS proteins. The identification of several one-to-one orthologs of human G proteins at the levels of families, subfamilies and types and of homologs of the human RGS proteins suggests an evolutionarily conserved structure function relationship of the GPCR signaling mechanism in the chordates.

Conclusions

The C. intestinalis genome encodes a highly conserved, albeit, limited repertoire of the heterotrimeric G protein complexes with the size of subunit types comparable with that in lower eukaryotes.     

Genome-wide exploration of the molecular evolution and regulatory network of mitogen-activated protein kinase cascades upon multiple stresses in Brachypodium distachyon

Background

Brachypodium distachyon is emerging as a widely recognized model plant that has very close relations with several economically important Poaceae species. MAPK cascade is known to be an evolutionarily conserved signaling module involved in multiple stresses. Although the gene sequences of MAPK and MAPKK family have been fully identified in B. distachyon, the information related to the upstream MAPKKK gene family especially the regulatory network among MAPKs, MAPKKs and MAPKKKs upon multiple stresses remains to be understood.

Results

In this study, we have identified MAPKKKs which belong to the biggest gene family of MAPK cascade kinases. We have systematically investigated the evolution of whole MAPK cascade kinase gene family in terms of gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication analysis. Our results showed that most BdMAPK cascade kinases were located at the low-CpG-density region, and the clustered members in each group shared similar structures of the genes and proteins. Synteny analysis showed that 62 or 21 pairs of duplicated orthologs were present between B. distachyon and Oryza sativa, or between B. distachyon and Arabidopsis thaliana respectively. Gene expression data revealed that BdMAPK cascade kinases were rapidly regulated by stresses and phytohormones. Importantly, we have constructed a regulation network based on co-expression patterns of the expression profiles upon multiple stresses performed in this study.

Conclusions

BdMAPK cascade kinases were involved in the signaling pathways of multiple stresses in B. distachyon. The network of co-expression regulation showed the most of duplicated BdMAPK cascade kinase gene orthologs demonstrated their convergent function, whereas few of them developed divergent function in the evolutionary process. The molecular evolution analysis of identified MAPK family genes and the constructed MAPK cascade regulation network under multiple stresses provide valuable information for further investigation of the functions of BdMAPK cascade kinase genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1452-1) contains supplementary material, which is available to authorized users.     

UDP-glucose dehydrogenase modulates proteoglycan synthesis in articular chondrocytes: its possible involvement and regulation in osteoarthritis

Introduction

The objective of this study was to investigate the possible role of UDP-glucose dehydrogenase (UGDH) in osteoarthritis (OA) and uncover whether, furthermore how interleukin-1beta (IL-1β) affects UGDH gene expression.

Methods

UGDH specific siRNAs were applied to determine the role of UGDH in proteoglycan (PG) synthesis in human articular chondrocytes. Protein levels of UGDH and Sp1 in human and rat OA cartilage were detected. Then, human primary chondrocytes were treated with IL-1β to find out whether and how IL-1β could regulate the gene expression of UGDH and its trans-regulators, that is Sp1, Sp3 and c-Krox. Finally, p38 mitogen-activated protein kinase (MAPK) inhibitor SB203580 and stress-activated protein kinase/c-Jun N-terminal kinase (SAP/JNK) inhibitor SP600125 were used to pick out the pathway that mediated the IL-1β-modulated PGs synthesis and gene expression of UGDH, Sp1, Sp3 and c-Krox.

Results

UGDH specific siRNAs markedly inhibited UGDH mRNA and protein expression, and thus led to an obvious suppression of PGs synthesis in human articular chondrocytes. UGDH protein level in human and rat OA cartilage were much lower than the corresponding controls and negatively correlated to the degree of OA. Decrease in Sp1 protein level was also observed in human and rat OA cartilage respectively. Meanwhile, IL-1β suppressed UGDH gene expression in human articular chondrocytes in the late phase, which also modulated gene expression of Sp1, Sp3 and c-Krox and increased both Sp3/Sp1 and c-Krox/Sp1 ratio. Moreover, the inhibition of SAP/JNK and p38 MAPK pathways both resulted in an obvious attenuation of the IL-1β-induced suppression on the UGDH gene expression.

Conclusions

UGDH is essential in the PGs synthesis of articular chondrocytes, while the suppressed expression of UGDH might probably be involved in advanced OA, partly due to the modulation of p38 MAPK and SAP/JNK pathways and its trans-regulators by IL-1β.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-014-0484-2) contains supplementary material, which is available to authorized users.