Origins, Lineage-Specific Expansions, and Multiple Losses of Tyrosine Kinases in Eukaryotes

The draft genome of Cochliopodium minus reveals a complete meiosis toolkit and provides insight into the evolution of sexual mechanisms in Amoebozoa

To date, genomic analyses in amoebozoans have been mostly limited to model organisms or medically important lineages. Consequently, the vast diversity of Amoebozoa genomes remain unexplored. A draft genome of Cochliopodium minus, an amoeba characterized by extensive cellular and nuclear fusions, is presented. C. minus has been a subject of recent investigation for its unusual sexual behavior. Cochliopodium's sexual activity occurs during vegetative stage making it an ideal model for studying sexual development, which is sorely lacking in the group. Here we generate a C. minus draft genome assembly. From this genome, we detect a substantial number of lateral gene transfer (LGT) instances from bacteria (15%), archaea (0.9%) and viruses (0.7%) the majority of which are detected in our transcriptome data. We identify the complete meiosis toolkit genes in the C. minus genome, as well as the absence of several key genes involved in plasmogamy and karyogamy. Comparative genomics of amoebozoans reveals variation in sexual mechanism exist in the group. Similar to complex eukaryotes, C. minus (some amoebae) possesses Tyrosine kinases and duplicate copies of SPO11. We report a first example of alternative splicing in a key meiosis gene and draw important insights on molecular mechanism of sex in C. minus using genomic and transcriptomic data.

Molecular and Chemical Screening for Inherent Disease Resistance Factors of Norway Spruce (Picea abies) Clones Against Conifer Stem Rot Pathogen Heterobasidion parviporum

Root and stem rot of conifer trees caused by Heterobasidion annosum species complex leads to huge economic losses in Europe, yet not much is known about the molecular and chemical basis for host resistance. To identify inherent chemical or molecular markers in clones found to be either resistant or susceptible, we sampled needle tissues of all the clones before pathogen inoculation. We conducted a short-term resistance screening by using the pathogen H. parviporum to inoculate 70 Norway spruce clones. Based on lesion size, subsets of highly susceptible and resistant clones were further analyzed. Terpene detection and RNA sequencing were performed to explore inherent variations in genotypes differing in resistance to pathogenic challenge at chemical and transcriptional levels. A negative correlation emerged between resistance and growth. Terpene profiles of resistant clones showed higher content of monoterpenes and sesquiterpenes, with concomitant increased transcript abundance of genes involved in the terpenoid pathway. A set of upregulated genes relevant to flavonoid biosynthesis was observed in resistant genotypes, whereas higher transcripts of lignin biosynthetic genes were prevalent in susceptible clones. Genes involved in flavonoid and lignin biosynthesis as well as terpene content may have a role in facilitating resistance of Norway spruce against H. parviporum. Our results provide strong support on the feasibility of sampling needle tissues before pathogen inoculation, and the approach could be of value for large-scale screening of novel biomarkers for durable resistance. The additional insights could form a basis for further research on resistance screening in this pathosystem.

Evolution of functional diversity in the holozoan tyrosine kinome

The emergence of multicellularity is strongly correlated with the expansion of tyrosine kinases, a conserved family of signaling enzymes that regulates pathways essential for cell-to-cell communication. Although tyrosine kinases have been classified from several model organisms, a molecular-level understanding of tyrosine kinase evolution across all holozoans is currently lacking. Using a hierarchical sequence constraint-based classification of diverse holozoan tyrosine kinases, we construct a new phylogenetic tree that identifies two ancient clades of cytoplasmic and receptor tyrosine kinases separated by the presence of an extended insert segment in the kinase domain connecting the D and E-helices. Present in nearly all receptor tyrosine kinases, this fast-evolving insertion imparts diverse functionalities, such as post-translational modification sites and regulatory interactions. Eph and EGFR receptor tyrosine kinases are two exceptions which lack this insert, each forming an independent lineage characterized by unique functional features. We also identify common constraints shared across multiple tyrosine kinase families which warrant the designation of three new subgroups: Src module (SrcM), insulin receptor kinase-like (IRKL), and fibroblast, platelet-derived, vascular, and growth factor receptors (FPVR). Subgroup-specific constraints reflect shared autoinhibitory interactions involved in kinase conformational regulation. Conservation analyses describe how diverse tyrosine kinase signaling functions arose through the addition of family-specific motifs upon subgroup-specific features and coevolving protein domains. We propose the oldest tyrosine kinases, IRKL, SrcM, and Csk, originated from unicellular premetazoans and were coopted for complex multicellular functions. The increased frequency of oncogenic variants in more recent tyrosine kinases suggests that lineage-specific functionalities are selectively altered in human cancers.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Genomic Survey of Tyrosine Kinases Repertoire in Electrophorus electricus With an Emphasis on Evolutionary Conservation and Diversification

Tyrosine kinases (TKs) play key roles in the regulation of multicellularity in organisms and involved primarily in cell growth, differentiation, and cell-to-cell communication. Genome-wide characterization of TKs has been conducted in many metazoans; however, systematic information regarding this superfamily in Electrophorus electricus (electric eel) is still lacking. In this study, we identified 114 TK genes in the E electricus genome and investigated their evolution, molecular features, and domain architecture using phylogenetic profiling to gain a better understanding of their similarities and specificity. Our results suggested that the electric eel TK (EeTK) repertoire was shaped by whole-genome duplications (WGDs) and tandem duplication events. Compared with other vertebrate TKs, gene members in Jak, Src, and EGFR subfamily duplicated specifically, but with members lost in Eph, Axl, and Ack subfamily in electric eel. We also conducted an exhaustive survey of TK genes in genomic databases, identifying 1674 TK proteins in 31 representative species covering all the main metazoan lineages. Extensive evolutionary analysis indicated that TK repertoire in vertebrates tended to be remarkably conserved, but the gene members in each subfamily were very variable. Comparative expression profile analysis showed that electric organ tissues and muscle shared a similar pattern with specific highly expressed TKs (ie, epha7, musk, jak1, and pdgfra), suggesting that regulation of TKs might play an important role in specifying an electric organ identity from its muscle precursor. We further identified TK genes exhibiting tissue-specific expression patterns, indicating that members in TKs participated in subfunctionalization representing an evolutionary divergence required for the performance of different tissues. This work generates valuable information for further gene function analysis and identifying candidate TK genes reflecting their unique tissue-function specializations in electric eel.

Specificity and regulation of phosphotyrosine signaling through SH2 domains

Phosphotyrosine (pY) signaling is instrumental to numerous cellular processes. pY recognition occurs through specialized protein modules, among which the Src-homology 2 (SH2) domain is the most common. SH2 domains are small protein modules with an invariant fold, and are present in more than a hundred proteins with different function. Here we ask the question of how such a structurally conserved, small protein domain can recognize distinct phosphopeptides with the breath of binding affinity, specificity and kinetic parameters necessary for proper control of pY-dependent signaling and rapid cellular response. We review the current knowledge on structure, thermodynamics and kinetics of SH2-phosphopeptide complexes and conclude that selective phosphopeptide recognition is governed by both structure and dynamics of the SH2 domain, as well as by the kinetics of the binding events. Further studies on the thermodynamic and kinetic properties of SH2-phosphopeptide complexes, beyond their structure, are required to understand signaling regulation.

Evaluation of oxidative stress responses in human circulating blood cells after imatinib mesylate treatment - Implications to its mechanism of action

Imatinib mesylate (IM) is the first developed protein kinase inhibitor and recently it has topped consumption rates among targeted and total anticancer drugs. Although there are indications that IM possesses cyto/genotoxic activities against normal non-target cells as well, there is a lack of information regarding the underlying mechanism involved in those actions. Therefore, we aimed to evaluate the response of human circulating blood cells towards oxidative stress after IM treatment (0.0001–10 µg/mL) in vitro. Based on the results, IM had an influence on all of the oxidative stress parameters tested. Lower concentrations of IM induced an increase of glutathione level, following its decrease at higher IM concentrations indicating impairment in oxidative stress defences. Concomitant to a glutathione decrease, an increase of malondialdehyde and protein carbonyls level was observed indicating oxidative damage of lipids and proteins. The observed effects overlapped with the observed formation of oxidative base damage detected by formamidopyrimidine-DNA glycosylase modified-comet assay indicating that IM managed to induce oxidative DNA damage. Our results provide novelty in their mechanistic approach to IM-induced toxicity in non-target cells and suggest that IM can affect blood cells and induce oxidative stress.

Engineered SH2 domains with tailored specificities and enhanced affinities for phosphoproteome analysis

Protein phosphorylation is the most abundant post-translational modification in cells. Src homology 2 (SH2) domains specifically recognize phosphorylated tyrosine (pTyr) residues to mediate signaling cascades. A conserved pocket in the SH2 domain binds the pTyr side chain and the EF and BG loops determine binding specificity. By using large phage-displayed libraries, we engineered the EF and BG loops of the Fyn SH2 domain to alter specificity. Engineered SH2 variants exhibited distinct specificity profiles and were able to bind pTyr sites on the epidermal growth factor receptor, which were not recognized by the wild-type Fyn SH2 domain. Furthermore, mass spectrometry showed that SH2 variants with additional mutations in the pTyr-binding pocket that enhanced affinity were highly effective for enrichment of diverse pTyr peptides within the human proteome. These results showed that engineering of the EF and BG loops could be used to tailor SH2 domain specificity, and SH2 variants with diverse specificities and high affinities for pTyr residues enabled more comprehensive analysis of the human phosphoproteome. STATEMENT: Src Homology 2 (SH2) domains are modular domains that recognize phosphorylated tyrosine embedded in proteins, transducing these post-translational modifications into cellular responses. Here we used phage display to engineer hundreds of SH2 domain variants with altered binding specificities and enhanced affinities, which enabled efficient and differential enrichment of the human phosphoproteome for analysis by mass spectrometry. These engineered SH2 domain variants will be useful tools for elucidating the molecular determinants governing SH2 domains binding specificity and for enhancing analysis and understanding of the human phosphoproteome.

Genome-Wide Identification and Characterization of Tyrosine Kinases in the Silkworm, Bombyx mori

The tyrosine kinases (TKs) are important parts of metazoan signaling pathways and play significant roles in cell growth, development, apoptosis and disease. Genome-wide characterization of TKs has been conducted in many metazoans, however, systematic information about this family in Lepidoptera is still lacking. We retrieved 33 TK-encoding genes in silkworm and classified them into 25 subfamilies by sequence analysis, without members in AXL, FRK, PDGFR, STYK1 and TIE subfamilies. Although domain sequences in each subfamily are conserved, TKs in vertebrates tend to be remarkably conserved and stable. Our results of phylogenetic analysis supported the previous conclusion for the second major expansion of TK family. Gene-Ontology (GO) analysis revealed that a higher proportion of BmTKs played roles in binding, catalysis, signal transduction, metabolism, biological regulation and response to stimulus, compared to all silkworm genes annotated in GO. Moreover, the expression profile analysis of BmTKs among multiple tissues and developmental stages demonstrated that many genes exhibited stage-specific and/or sex-related expression during embryogenesis, molting and metamorphosis, and that 8 BmTKs presented tissue-specific high expression. Our study provides systematic description of silkworm tyrosine kinases, and may also provide further insights into metazoan TKs and assist future studies addressing their functions.

Auxin and Gibberellins Are Required for the Receptor-Like Kinase ERECTA Regulated Hypocotyl Elongation in Shade Avoidance in Arabidopsis

Plants use shade avoidance strategy to escape the canopy shade when grown under natural conditions. Previous studies showed that the Arabidopsis receptor-like kinase ERECTA (ER) is involved in shade avoidance syndrome. However, the mechanisms of ER in modulating SAR by promoting hypocotyl elongation are unknown yet. Here, we report that ER regulated hypocotyl elongation in shade avoidance requires auxin and gibberellins (GAs). The T-DNA insertional ER mutant er-3 shows a less hypocotyl length than that in Col-0 wild type. Promoter::GUS staining analysis shows that ER and its paralogous genes ERECTA-LIKE1 (ERL1) and ERECTA-LIKE2 (ERL2) are differentially expressed in the seedlings, of which only ER is most obviously upregulated in the hypocotyl by shade treatment. Exogenous feeding assay by using media-application with vertical-grown of Arabidopsis seedlings showed that the hypocotyl length of er-3 is partially promoted by indol-3-acetic acid (IAA), while it is relatively insensitive of er-3 to various concentrations of IAA than that of Col-0. Hypocotyl elongation of er-3 is promoted similar to that of Col-0 by high temperature in the white light condition, but the elongation was not significantly affected by the treatment of the auxin transport inhibitor 1-N-naphthylphthalamic acid (NPA). Exogenous GA3 increased the hypocotyl elongation of both er-3 and the wild type in the shade condition, and the GA3 biosynthesis inhibitor paclobutrazol (PAC) severely inhibits the hypocotyl elongation of Col-0 and er-3. Further analysis showed that auxin biosynthesis inhibitors yucasin and L-kynurenine remarkably inhibited the hypocotyl elongation of er-3 while yucasin shows a more severe inhibition to er-3 than Col-0. Relative expression of genes regulating auxin homeostasis and signaling, and GA homeostasis is less in er-3 than that in Col-0. Furthermore, genetic evidences show that ER regulated hypocotyl elongation is dependent of PHYTOCHROME B (PHYB). Overall, we propose that ER regulated shade avoidance by promoting hypocotyl elongation is PHYB-dependent and requires auxin and GAs.

Phosphotyrosine signalling and the origin of animal multicellularity

The evolution of multicellular animals (i.e. metazoans) from a unicellular ancestor is one of the most important yet least understood evolutionary transitions. Historically, given its indispensable functions in intercellular communication and exclusive presence in metazoans, phosphotyrosine (pTyr) signalling was considered a metazoan-specific evolutionary innovation that might have contributed to the origin of metazoan multicellularity. However, recent studies have led to a new understanding of pTyr signalling evolution and its role in the metazoan origin. Sequence analyses have unravelled a much earlier emergence of pTyr signalling in eukaryotic evolution. Even so, several distinct properties of holozoan pTyr signalling may have paved the way for a hypothesized functional transition of pTyr signalling at the multicellular origin, from environmental sensing to intercellular communication, and for it to evolve as a powerful intercellular signalling system for multicellularity. Biochemical analyses of premetazoan pTyr signalling components have further revealed the premetazoan origin of many key features of metazoan pTyr signalling, and the metazoan establishment of others, including the Csk-mediated negative regulation of the activity of Src, a conserved tyrosine kinase in the Holozoa. Finally, potential future directions are discussed, with a stress on the biological functions of premetazoan pTyr signalling via newly developed gene manipulation tools in non-animal holozoans.

Early emergence of negative regulation of the tyrosine kinase Src by the C-terminal Src kinase

Stringent regulation of tyrosine kinase activity is essential for normal cellular function. In humans, the tyrosine kinase Src is inhibited via phosphorylation of its C-terminal tail by another kinase, C-terminal Src kinase (Csk). Although Src and Csk orthologs are present across holozoan organisms, including animals and protists, the Csk-Src negative regulatory mechanism appears to have evolved gradually. For example, in choanoflagellates, Src and Csk are both active, but the negative regulatory mechanism is reportedly absent. In filastereans, a protist clade closely related to choanoflagellates, Src is active, but Csk is apparently inactive. In this study, we use a combination of bioinformatics, in vitro kinase assays, and yeast-based growth assays to characterize holozoan Src and Csk orthologs. We show that, despite appreciable differences in domain architecture, Csk from Corallochytrium limacisporum, a highly diverged holozoan marine protist, is active and can inhibit Src. However, in comparison with other Csk orthologs, Corallochytrium Csk displays broad substrate specificity and inhibits Src in an activity-independent manner. Furthermore, in contrast to previous studies, we show that Csk from the filasterean Capsaspora owczarzaki is active and that the Csk-Src negative regulatory mechanism is present in Csk and Src proteins from C. owczarzaki and the choanoflagellate Monosiga brevicollis Our results suggest that negative regulation of Src by Csk is more ancient than previously thought and that it might be conserved across all holozoan species.

Identification and characterization of tyrosine kinases in anole lizard indicate the conserved tyrosine kinase repertoire in vertebrates

The tyrosine kinases (TKs) play principal roles in regulation of multicellular aspects of the organism and are implicated in many cancer types and congenital disorders. The anole lizard has recently been introduced as a model organism for laboratory-based studies of organismal function and field studies of ecology and evolution. However, the TK family of anole lizard has not been systematically identified and characterized yet. In this study, we identified 82 TK-encoding genes in the anole lizard genome and classified them into 28 subfamilies through phylogenetic analysis, with no member from ROS and STYK1 subfamilies identified. Although TK domain sequences and domain organization in each subfamily were conserved, the total number of TKs in different species was much variable. In addition, extensive evolutionary analysis in metazoans indicated that TK repertoire in vertebrates tends to be remarkably stable. Phylogenetic analysis of Eph subfamily indicated that the divergence of EphA and EphB occurred prior to the whole genome duplication (WGD) but after the split of Urochordates and vertebrates. Moreover, the expression pattern analysis of lizard TK genes among 9 different tissues showed that 14 TK genes exhibited tissue-specific expression and 6 TK genes were widely expressed. Comparative analysis of TK expression suggested that the tissue specifically expressed genes showed different expression pattern but the widely expressed genes showed similar pattern between anole lizard and human. These results may provide insights into the evolutionary diversification of animal TK genes and would aid future studies on TK protein regulation of key growth and developmental processes.

Hundreds of dual-stage antimalarial molecules discovered by a functional gametocyte screen

Plasmodium falciparum stage V gametocytes are responsible for parasite transmission, and drugs targeting this stage are needed to support malaria elimination. We here screen the Tres Cantos Antimalarial Set (TCAMS) using the previously developed P. falciparum female gametocyte activation assay (Pf FGAA), which assesses stage V female gametocyte viability and functionality using Pfs25 expression. We identify over 400 compounds with activities <2 μM, chemically classified into 57 clusters and 33 singletons. Up to 68% of the hits are chemotypes described for the first time as late-stage gametocyte-targeting molecules. In addition, the biological profile of 90 compounds representing the chemical diversity is assessed. We confirm in vitro transmission-blocking activity of four of the six selected molecules belonging to three distinct scaffold clusters. Overall, this TCAMS gametocyte screen provides 276 promising antimalarial molecules with dual asexual/sexual activity, representing starting points for target identification and candidate selection.

There is a need for Plasmodium transmission blocking drugs for malaria elimination. Here, Miguel-Blanco et al. screen >10,000 compounds against stage V female gametocytes, identify active compounds belonging to 57 chemotypes and confirm transmission blocking activity of four selected compounds in vitro.

Classification and Lineage Tracing of SH2 Domains Throughout Eukaryotes

Today there exists a rapidly expanding number of sequenced genomes. Cataloging protein interaction domains such as the Src Homology 2 (SH2) domain across these various genomes can be accomplished with ease due to existing algorithms and predictions models. An evolutionary analysis of SH2 domains provides a step towards understanding how SH2 proteins integrated with existing signaling networks to position phosphotyrosine signaling as a crucial driver of robust cellular communication networks in metazoans. However organizing and tracing SH2 domain across organisms and understanding their evolutionary trajectory remains a challenge. This chapter describes several methodologies towards analyzing the evolutionary trajectory of SH2 domains including a global SH2 domain classification system, which facilitates annotation of new SH2 sequences essential for tracing the lineage of SH2 domains throughout eukaryote evolution. This classification utilizes a combination of sequence homology, protein domain architecture and the boundary positions between introns and exons within the SH2 domain or genes encoding these domains. Discrete SH2 families can then be traced across various genomes to provide insight into its origins. Furthermore, additional methods for examining potential mechanisms for divergence of SH2 domains from structural changes to alterations in the protein domain content and genome duplication will be discussed. Therefore a better understanding of SH2 domain evolution may enhance our insight into the emergence of phosphotyrosine signaling and the expansion of protein interaction domains.

Genomic and Post-Translational Modification Analysis of Leucine-Rich-Repeat Receptor-Like Kinases in Brassica rapa

Among several receptor-like kinases (RLKs), leucine-rich-repeat receptor-like kinases (LRR-RLKs) are a major group of genes that play crucial roles in growth, development and stress responses in plant systems. Given that they have several functional roles, it is important to investigate their roles in Brassica rapa. In the present study, 303 LRR-RLKs were identified in the genome of B. rapa and comparative phylogenetic analysis of 1213 combined LRR-RLKs of B. rapa, Arabidopsis thaliana, Oryza sativa and Populus trichocarpa helped us to categorize the gene family into 15 subfamilies based on their sequence and structural similarities. The chromosome localizations of 293 genes allowed the prediction of duplicates, and motif conservation and intron/exon patterns showed differences among the B. rapa LRR-RLK (BrLRR-RLK) genes. Additionally, computational function annotation and expression analysis was used to predict their possible functional roles in the plant system. Biochemical results for 11 selected genes showed variations in phosphorylation activity. Interestingly, BrBAK1 showed strong auto-phosphorylation and trans-phosphorylation on its tyrosine and threonine residues compared with AtBAK1 in previous studies. The AtBAK1 receptor kinase is involved in plant growth and development, plant innate immunity, and programmed cell death, and our results suggest that BrBAK1 might also be involved in the same functions. Another interesting result was that BrBAK1, BrBRI1, BrPEPR1 and BrPEPR2 showed activity with both anti-phosphotyrosine and anti-phosphothreonine antibodies, indicating that they might have dual-specificity kinase activity. This study provides comprehensive results for the BrLRR-RLKs, revealing expansion of the gene family through gene duplications, structural similarities and variations among the genes, and potential functional roles according to gene ontology, transcriptome profiling and biochemical analysis.

Genome of Rhodnius prolixus, an insect vector of Chagas disease, reveals unique adaptations to hematophagy and parasite infection

Rhodnius prolixus not only has served as a model organism for the study of insect physiology, but also is a major vector of Chagas disease, an illness that affects approximately seven million people worldwide. We sequenced the genome of R. prolixus, generated assembled sequences covering 95% of the genome (∼ 702 Mb), including 15,456 putative protein-coding genes, and completed comprehensive genomic analyses of this obligate blood-feeding insect. Although immune-deficiency (IMD)-mediated immune responses were observed, R. prolixus putatively lacks key components of the IMD pathway, suggesting a reorganization of the canonical immune signaling network. Although both Toll and IMD effectors controlled intestinal microbiota, neither affected Trypanosoma cruzi, the causal agent of Chagas disease, implying the existence of evasion or tolerance mechanisms. R. prolixus has experienced an extensive loss of selenoprotein genes, with its repertoire reduced to only two proteins, one of which is a selenocysteine-based glutathione peroxidase, the first found in insects. The genome contained actively transcribed, horizontally transferred genes from Wolbachia sp., which showed evidence of codon use evolution toward the insect use pattern. Comparative protein analyses revealed many lineage-specific expansions and putative gene absences in R. prolixus, including tandem expansions of genes related to chemoreception, feeding, and digestion that possibly contributed to the evolution of a blood-feeding lifestyle. The genome assembly and these associated analyses provide critical information on the physiology and evolution of this important vector species and should be instrumental for the development of innovative disease control methods.

Phylogenetic Profiling for Probing the Modular Architecture of the Human Genome

Information about functional connections between genes can be derived from patterns of coupled loss of their homologs across multiple species. This comparative approach, termed phylogenetic profiling, has been successfully used to infer genetic interactions in bacteria and eukaryotes. Rapid progress in sequencing eukaryotic species has enabled the recent phylogenetic profiling of the human genome, resulting in systematic functional predictions for uncharacterized human genes. Importantly, groups of co-evolving genes reveal widespread modularity in the underlying genetic network, facilitating experimental analyses in human cells as well as comparative studies of conserved functional modules across species. This strategy is particularly successful in identifying novel metabolic proteins and components of multi-protein complexes. The targeted sequencing of additional key eukaryotes and the incorporation of improved methods to generate and compare phylogenetic profiles will further boost the predictive power and utility of this evolutionary approach to the functional analysis of gene interaction networks.

Identification and expression analysis of the LRR-RLK gene family in tomato (Solanum lycopersicum) Heinz 1706

As the largest subfamily of receptor-like kinases (RLKs), leucine-rich repeat receptor-like kinases (LRR-RLKs) regulate the growth, development, and stress responses of plants. Through a reiterative process of sequence analysis and re-annotation, 234 LRR-RLK genes were identified in the genome of tomato (Solanum lycopersicum) ‘Heinz 1706’, which were further grouped into 10 major groups based on their sequence similarity. In comparison to the significant role of tandem duplication in the expansion process of this gene family in other species, only approximately 12% (29 out of 234) of SlLRR-RLK genes arose from tandem duplication. Using the multiple expectation maximization for motif elicitation (MEME) method, the motif composition and arrangement were found to be variably conserved within each SlLRR-RLK group, indicating their different extent of functional divergence. Expression profiling analyses by qRT-PCR data revealed that SlLRR-RLK genes were differentially expressed in various tomato organs and tissues, and some SlLRR-RLK genes exhibited preferential expression in fruits at distinct developmental stages, suggesting that SlLRR-RLK may take important roles in fruit development and ripening process. The results of this study provide an overview of the LRR-RLK gene family in tomato Heinz 1706, one important species of Solanaceae, and will be helpful for future functional analysis of this important protein family in fleshy fruit-bearing species.

Testing whether metazoan tyrosine loss was driven by selection against promiscuous phosphorylation

Protein tyrosine phosphorylation is a key regulatory modification in metazoans, and the corresponding kinase enzymes have diversified dramatically. This diversification is correlated with a genome-wide reduction in protein tyrosine content, and it was recently suggested that this reduction was driven by selection to avoid promiscuous phosphorylation that might be deleterious. We tested three predictions of this intriguing hypothesis. 1) Selection should be stronger on residues that are more likely to be phosphorylated due to local solvent accessibility or structural disorder. 2) Selection should be stronger on proteins that are more likely to be promiscuously phosphorylated because they are abundant. We tested these predictions by comparing distributions of tyrosine within and among human and yeast orthologous proteins. 3) Selection should be stronger against mutations that create tyrosine versus remove tyrosine. We tested this prediction using human population genomic variation data. We found that all three predicted effects are modest for tyrosine when compared with the other amino acids, suggesting that selection against deleterious phosphorylation was not dominant in driving metazoan tyrosine loss.

Evolution of the JAK-STAT pathway

The JAK-STAT pathway represents a finely tuned orchestra capable of rapidly facilitating an exquisite symphony of responses from a complex array of extracellular signals. This review explores the evolution of the JAK-STAT pathway: the origins of the individual domains from which it is constructed, the formation of individual components from these basic building blocks, the assembly of the components into a functional pathway, and the subsequent reiteration of this basic template to fulfill a variety of roles downstream of cytokine receptors.

Malaria protein kinase CK2 (PfCK2) shows novel mechanisms of regulation

Casein kinase 2 (protein kinase CK2) is a conserved eukaryotic serine/theronine kinase with multiple substrates and roles in the regulation of cellular processes such as cellular stress, cell proliferation and apoptosis. Here we report a detailed analysis of the Plasmodium falciparum CK2, PfCK2, demonstrating that this kinase, like the mammalian orthologue, is a dual specificity kinase able to phosphorylate at both serine and tyrosine. However, unlike the human orthologue that is auto-phosphorylated on tyrosine within the activation loop, PfCK2 shows no activation loop auto-phosphorylation but rather is auto-phosphorylated at threonine 63 within subdomain I. Phosphorylation at this site in PfCK2 is shown here to regulate the intrinsic kinase activity of PfCK2. Furthermore, we generate an homology model of PfCK2 in complex with the known selective protein kinase CK2 inhibitor, quinalizarin, and in so doing identify key co-ordinating residues in the ATP binding pocket that could aid in designing selective inhibitors to PfCK2.

Identification of a fungi-specific lineage of protein kinases closely related to tyrosine kinases

Tyrosine kinases (TKs) specifically catalyze the phosphorylation of tyrosine residues in proteins and play essential roles in many cellular processes. Although TKs mainly exist in animals, recent studies revealed that some organisms outside the Opisthokont clade also contain TKs. The fungi, as the sister group to animals, are thought to lack TKs. To better understand the origin and evolution of TKs, it is important to investigate if fungi have TK or TK-related genes. We therefore systematically identified possible TKs across the fungal kingdom by using the profile hidden Markov Models searches and phylogenetic analyses. Our results confirmed that fungi lack the orthologs of animal TKs. We identified a fungi-specific lineage of protein kinases (FslK) that appears to be a sister group closely related to TKs. Sequence analysis revealed that members of the FslK clade contain all the conserved protein kinase sub-domains and thus are likely enzymatically active. However, they lack key amino acid residues that determine TK-specific activities, indicating that they are not true TKs. Phylogenetic analysis indicated that the last common ancestor of fungi may have possessed numerous members of FslK. The ancestral FslK genes were lost in Ascomycota and Ustilaginomycotina and Pucciniomycotina of Basidiomycota during evolution. Most of these ancestral genes, however, were retained and expanded in Agaricomycetes. The discovery of the fungi-specific lineage of protein kinases closely related to TKs helps shed light on the origin and evolution of TKs and also has potential implications for the importance of these kinases in mushroom fungi.

Earliest Holozoan expansion of phosphotyrosine signaling

Phosphotyrosine (pTyr) signaling is involved in development and maintenance of metazoans’ multicellular body through cell-to-cell communication. Tyrosine kinases (TKs), tyrosine phosphatases, and other proteins relaying the signal compose the cascade. Domain architectures of the pTyr signaling proteins are diverse in metazoans, reflecting their complex intercellular communication. Previous studies had shown that the metazoan-type TKs, as well as other pTyr signaling proteins, were already diversified in the common ancestor of metazoans, choanoflagellates, and filastereans (which are together included in the clade Holozoa) whereas they are absent in fungi and other nonholozoan lineages. However, the earliest-branching holozoans Ichthyosporea and Corallochytrea, as well as the two fungi-related amoebae Fonticula and Nuclearia, have not been studied. Here, we analyze the complete genome sequences of two ichthyosporeans and Fonticula, and RNAseq data of three additional ichthyosporeans, one corallochytrean, and Nuclearia. Both the ichthyosporean and corallochytrean genomes encode a large variety of receptor TKs (RTKs) and cytoplasmic TKs (CTKs), as well as other pTyr signaling components showing highly complex domain architectures. However, Nuclearia and Fonticula have no TK, and show much less diversity in other pTyr signaling components. The CTK repertoires of both Ichthyosporea and Corallochytrea are similar to those of Metazoa, Choanoflagellida, and Filasterea, but the RTK sets are totally different from each other. The complex pTyr signaling equipped with positive/negative feedback mechanism likely emerged already at an early stage of holozoan evolution, yet keeping a high evolutionary plasticity in extracellular signal reception until the co-option of the system for cell-to-cell communication in metazoans.

A two-hybrid-receptor assay demonstrates heteromer formation as switch-on for plant immune receptors

Receptor kinases sense extracellular signals and trigger intracellular signaling and physiological responses. However, how does signal binding to the extracellular domain activate the cytoplasmic kinase domain? Activation of the plant immunoreceptor Flagellin sensing2 (FLS2) by its bacterial ligand flagellin or the peptide-epitope flg22 coincides with rapid complex formation with a second receptor kinase termed brassinosteroid receptor1 associated kinase1 (BAK1). Here, we show that the receptor pair of FLS2 and BAK1 is also functional when the roles of the complex partners are reversed by swapping their cytosolic domains. This reciprocal constellation prevents interference by redundant partners that can partially substitute for BAK1 and demonstrates that formation of the heteromeric complex is the molecular switch for transmembrane signaling. A similar approach with swaps between the Elongation factor-Tu receptor and BAK1 also resulted in a functional receptor/coreceptor pair, suggesting that a "two-hybrid-receptor assay" is of more general use for studying heteromeric receptor complexes.

Role of class III phosphatidylinositol 3-kinase during programmed nuclear death of Tetrahymena thermophila

Programmed nuclear death (PND) in the ciliate protozoan Tetrahymena thermophila is a novel type of autophagy that occurs during conjugation, in which only the parental somatic macronucleus is destined to die and is then eliminated from the progeny cytoplasm. Other coexisting nuclei, however, such as new micro- and macronuclei are unaffected. PND starts with condensation in the nucleus followed by apoptotic DNA fragmentation, lysosomal acidification, and final resorption. Because of the peculiarity in the process and the absence of some ATG genes in this organism, the mechanism of PND has remained unclear. In this study, we focus on the role of class III phosphatidylinositol 3-kinase (PtdIns3K, corresponding to yeast Vps34) in order to identify central regulators of PND. We identified the sole Tetrahymena thermophila ortholog (TtVPS34) to yeast Vps34 and human PIK3C3 (the catalytic subunit of PtdIns3K), through phylogenetic analysis, and generated the gene knockdown mutant for functional analysis. Loss of TtVPS34 activity prevents autophagosome formation on the parental macronucleus, and this nucleus escapes from the lysosomal pathway. In turn, DNA fragmentation and final resorption of the nucleus are drastically impaired. These phenotypes are similar to the situation in the ATG8Δ mutants of Tetrahymena, implying an inextricable link between TtVPS34 and TtATG8s in controlling PND as well as general macroautophagy. On the other hand, TtVPS34 does not appear responsible for the nuclear condensation and does not affect the progeny nuclear development. These results demonstrate that TtVPS34 is critically involved in the nuclear degradation events of PND in autophagosome formation rather than with an involvement in commitment to the death program.

Genome-wide identification, characterization and expression analysis of populus leucine-rich repeat receptor-like protein kinase genes

Background

Leucine-rich repeat receptor-like kinases (LRR-RLKs) comprise the largest group within the receptor-like kinase (RLK) superfamily in plants. This gene family plays critical and diverse roles in plant growth, development and stress response. Although the LRR-RLK families in Arabidopsis and rice have been previously analyzed, no comprehensive studies have been performed on this gene family in tree species.

Results

In this work, 379 LRR-RLK genes were retrieved from the Populus trichocarpa genome and further grouped into 14 subfamilies based on their structural and sequence similarities. Approximately 82% (312 out of 379) of the PtLRR-RLK genes are located in segmental duplication blocks indicating the role of duplication process in the expansion of this gene family. The conservation and variation in motif composition and intron/exon arrangement among PtLRR-RLK subfamilies were analyzed to provide additional support for their phylogenetic relationship and more importantly to indicate the potential divergence in their functions. Expression profiling of PtLRR-RLKs showed that they were differentially expressed in different organs and tissues and some PtLRR-RLKs were specifically expressed in meristem tissues, which indicated their potential involvement in tissue development and differentiation. For most AtLRR-RLKs with defined functions, Populus homologues exhibiting similar expression patterns could be identified, which might indicate the functional conservation during evolution. Among 12 types of environmental cues analyzed by the genome-wide microarray data, PtLRR-RLKs showed specific responses to shoot organogenesis, wounding, low ammonium feeding, hypoxia and seasonal dormancy, but not to drought, re-watering after drought, flooding, AlCl₃ treatment and bacteria or fungi treatments.

Conclusions

This study provides the first comprehensive genomic analysis of the Populus LRR-RLK gene family. Segmental duplication contributes significantly to the expansion of this gene family. Populus and Arabidopsis LRR-RLK homologues not only share similar genetic structures but also exhibit comparable expression patterns which point to the possible functional conservation of these LRR-RLKs in two model systems. Transcriptome profiling provides the first insight into the functional divergence among PtLRR-RLK gene subfamilies and suggests that they might take important roles in growth and adaptation of tree species.

Genomic survey of premetazoans shows deep conservation of cytoplasmic tyrosine kinases and multiple radiations of receptor tyrosine kinases

The evolution of multicellular metazoans from a unicellular ancestor is one of the most important advances in the history of life. Protein tyrosine kinases play important roles in cell-to-cell communication, cell adhesion, and differentiation in metazoans; thus, elucidating their origins and early evolution is crucial for understanding the origin of metazoans. Although tyrosine kinases exist in choanoflagellates, few data are available about their existence in other premetazoan lineages. To unravel the origin of tyrosine kinases, we performed a genomic and polymerase chain reaction (PCR)-based survey of the genes that encode tyrosine kinases in the two described filasterean species, Capsaspora owczarzaki and Ministeria vibrans, the closest relatives to the Metazoa and Choanoflagellata clades. We present 103 tyrosine kinase-encoding genes identified in the whole genome sequence of C. owczarzaki and 15 tyrosine kinase-encoding genes cloned by PCR from M. vibrans. Through detailed phylogenetic analysis, comparison of the organizations of the protein domains, and resequencing and revision of tyrosine kinase sequences previously found in some whole genome sequences, we demonstrate that the basic repertoire of metazoan cytoplasmic tyrosine kinases was established before the divergence of filastereans from the Metazoa and Choanoflagellata clades. In contrast, the receptor tyrosine kinases diversified extensively in each of the filasterean, choanoflagellate, and metazoan clades. This difference in the divergence patterns between cytoplasmic tyrosine kinases and receptor tyrosine kinases suggests that receptor tyrosine kinases that had been used for receiving environmental cues were subsequently recruited as a communication tool between cells at the onset of metazoan multicellularity.

TEC and MAPK Kinase Signalling Pathways in T helper (TH) cell Development, TH2 Differentiation and Allergic Asthma

Significant advances in our understanding of the signalling events during T cell development and differentiation have been made in the past few decades. It is clear that ligation of the T cell receptor (TCR) triggers a series of proximal signalling cascades regulated by an array of protein kinases. These orchestrated and highly regulated series of events, with differential requirements of particular kinases, highlight the disparities between αβ+CD4+ T cells. Throughout this review we summarise both new and old studies, highlighting the role of Tec and MAPK in T cell development and differentiation with particular focus on T helper 2 (TH2) cells. Finally, as the allergy epidemic continues, we feature the role played by TH2 cells in the development of allergy and provide a brief update on promising kinase inhibitors that have been tested in vitro, in pre-clinical disease models in vivo and into clinical studies.

The SH2 domain-containing proteins in 21 species establish the provenance and scope of phosphotyrosine signaling in eukaryotes

The Src homology 2 (SH2) domains are participants in metazoan signal transduction, acting as primary mediators for regulated protein-protein interactions with tyrosine-phosphorylated substrates. Here, we describe the origin and evolution of SH2 domain proteins by means of sequence analysis from 21 eukaryotic organisms from the basal unicellular eukaryotes, where SH2 domains first appeared, through the multicellular animals and increasingly complex metazoans. On the basis of our results, SH2 domains and phosphotyrosine signaling emerged in the early Unikonta, and the numbers of SH2 domains expanded in the choanoflagellate and metazoan lineages with the development of tyrosine kinases, leading to rapid elaboration of phosphotyrosine signaling in early multicellular animals. Our results also indicated that SH2 domains coevolved and the number of the domains expanded alongside protein tyrosine kinases and tyrosine phosphatases, thereby coupling phosphotyrosine signaling to downstream signaling networks. Gene duplication combined with domain gain or loss produced novel SH2-containing proteins that function within phosphotyrosine signaling, which likely have contributed to diversity and complexity in metazoans. We found that intra- and intermolecular interactions within and between SH2 domain proteins increased in prevalence along with organismal complexity and may function to generate more highly connected and robust phosphotyrosine signaling networks.

Global kinomic and phospho-proteomic analyses of the human malaria parasite Plasmodium falciparum

The role of protein phosphorylation in the life cycle of malaria parasites is slowly emerging. Here we combine global phospho-proteomic analysis with kinome-wide reverse genetics to assess the importance of protein phosphorylation in Plasmodium falciparum asexual proliferation. We identify 1177 phosphorylation sites on 650 parasite proteins that are involved in a wide range of general cellular activities such as DNA synthesis, transcription and metabolism as well as key parasite processes such as invasion and cyto-adherence. Several parasite protein kinases are themselves phosphorylated on putative regulatory residues, including tyrosines in the activation loop of PfGSK3 and PfCLK3; we show that phosphorylation of PfCLK3 Y526 is essential for full kinase activity. A kinome-wide reverse genetics strategy identified 36 parasite kinases as likely essential for erythrocytic schizogony. These studies not only reveal processes that are regulated by protein phosphorylation, but also define potential anti-malarial drug targets within the parasite kinome.

Integrative genomic approaches highlight a family of parasite-specific kinases that regulate host responses

Apicomplexan parasites release factors via specialized secretory organelles (rhoptries, micronemes) that are thought to control host cell responses. In order to explore parasite-mediated modulation of host cell signaling pathways, we exploited a phylogenomic approach to characterize the Toxoplasma gondii kinome, defining a 44 member family of coccidian-specific secreted kinases, some of which have been previously implicated in virulence. Comparative genomic analysis suggests that "ROPK" genes are under positive selection, and expression profiling demonstrates that most are differentially expressed between strains and/or during differentiation. Integrating diverse genomic-scale analyses points to ROP38 as likely to be particularly important in parasite biology. Upregulating expression of this previously uncharacterized gene in transgenic parasites dramatically suppresses transcriptional responses in the infected cell. Specifically, parasite ROP38 downregulates host genes associated with MAPK signaling and the control of apoptosis and proliferation. These results highlight the value of integrative genomic approaches in prioritizing candidates for functional validation.

The kinome of Phytophthora infestans reveals oomycete-specific innovations and links to other taxonomic groups

BACKGROUND

Oomycetes are a large group of economically and ecologically important species. Its most notorious member is Phytophthora infestans, the cause of the devastating potato late blight disease. The life cycle of P. infestans involves hyphae which differentiate into spores used for dispersal and host infection. Protein phosphorylation likely plays crucial roles in these stages, and to help understand this we present here a genome-wide analysis of the protein kinases of P. infestans and several relatives. The study also provides new insight into kinase evolution since oomycetes are taxonomically distant from organisms with well-characterized kinomes.

RESULTS

Bioinformatic searches of the genomes of P. infestans, P. ramorum, and P. sojae reveal they have similar kinomes, which for P. infestans contains 354 eukaryotic protein kinases (ePKs) and 18 atypical kinases (aPKs), equaling 2% of total genes. After refining gene models, most were classifiable into families seen in other eukaryotes. Some ePK families are nevertheless unusual, especially the tyrosine kinase-like (TKL) group which includes large oomycete-specific subfamilies. Also identified were two tyrosine kinases, which are rare in non-metazoans. Several ePKs bear accessory domains not identified previously on kinases, such as cyclin-dependent kinases with integral cyclin domains. Most ePKs lack accessory domains, implying that many are regulated transcriptionally. This was confirmed by mRNA expression-profiling studies that showed that two-thirds vary significantly between hyphae, sporangia, and zoospores. Comparisons to neighboring taxa (apicomplexans, ciliates, diatoms) revealed both clade-specific and conserved features, and multiple connections to plant kinases were observed. The kinome of Hyaloperonospora arabidopsidis, an oomycete with a simpler life cycle than P. infestans, was found to be one-third smaller. Some differences may be attributable to gene clustering, which facilitates subfamily expansion (or loss) through unequal crossing-over.

CONCLUSION

The large sizes of the Phytophthora kinomes imply that phosphorylation plays major roles in their life cycles. Their kinomes also include many novel ePKs, some specific to oomycetes or shared with neighboring groups. Little experimentation to date has addressed the biological functions of oomycete kinases, but this should be stimulated by the structural, evolutionary, and expression data presented here. This may lead to targets for disease control.

The evolutionary mechanics of domain organization in proteomes and the rise of modularity in the protein world

Protein domains are compact evolutionary units of structure and function that usually combine in proteins to produce complex domain arrangements. In order to study their evolution, we reconstructed genome-based phylogenetic trees of architectures from a census of domain structure and organization conducted at protein fold and fold-superfamily levels in hundreds of fully sequenced genomes. These trees defined timelines of architectural discovery and revealed remarkable evolutionary patterns, including the explosive appearance of domain combinations during the rise of organismal lineages, the dominance of domain fusion processes throughout evolution, and the late appearance of a new class of multifunctional modules in Eukarya by fission of domain combinations. Our study provides a detailed account of the history and diversification of a molecular interactome and shows how the interplay of domain fusions and fissions defines an evolutionary mechanics of domain organization that is fundamentally responsible for the complexity of the protein world.

Genome analysis of the unicellular green alga Chlamydomonas reinhardtii Indicates an ancient evolutionary origin for key pattern recognition and cell-signaling protein families

The evolution of specific cell signaling and adhesion domains may have played an important role in the transition to a multicellular existence in the metazoans. Genomic analysis indicates that several signaling domains predominately found in animals are also present in the unicellular green alga, Chlamydomonas reinhardtii. A large group of proteins is present, containing scavenger receptor cysteine-rich (SRCR) and C-type lectin domains, which function in ligand binding and play key roles in the innate immune system of animals. Chlamydomonas also contains a large family of putative tyrosine kinases, suggesting an important role for phosphotyrosine signaling in the green algae. These important signaling domains may therefore be widespread among eukaryotes and most probably evolved in ancestral eukaryotes before the divergence of the Opisthokonts (the animal and fungal lineage).

The evolutionarily conserved arrangement of domains in SRC family kinases is important for substrate recognition

The SH3-SH2-kinase domain arrangement in nonreceptor tyrosine kinases has been conserved throughout evolution. For Src family kinases, the relative positions of the domains are important for enzyme regulation; they permit the assembly of Src kinases into autoinhibited conformations. The SH3 and SH2 domains of Src family kinases have an additional role in determining the substrate specificity of the kinase. We addressed the question of whether the domain arrangement of Src family kinases has a role in substrate specificity by producing mutants with alternative arrangements. Our results suggest that changes in the positions of domains can lead to specific changes in the phosphorylation of Sam68 and Cas by Src. Phosphorylation of Cas by several mutants triggers downstream signaling leading to cell migration. The placement of the SH2 domain with respect to the catalytic domain of Src appears to be especially important for proper substrate recognition, while the placement of the SH3 domain is more flexible. The results suggest that the involvement of the SH3 and SH2 domains in substrate recognition is one reason for the strict conservation of the SH3-SH2-kinase architecture.

Gene expansion and retention leads to a diverse tyrosine kinase superfamily in amphioxus

Tyrosine kinase (TK) proteins play a central role in cellular behavior and development of animals. The expansion of this superfamily is regarded as a key event in the evolution of the complex signaling pathways and gene networks of metazoans and is a prominent example of how shuffling of protein modules may generate molecular novelties. Using the intron/exon structure within the TK domain (TK intron code) as a complementary tool for the assignment of orthology and paralogy, we identified and studied the 118 TK proteins of the amphioxus Branchiostoma floridae genome to elucidate TK gene family evolution in metazoans and chordates in particular. Unlike all characterized metazoans to date, amphioxus has members of all known widespread TK families, with not a single loss. Putting amphioxus TKs in an evolutionary context, including new data from the cnidarian Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus, and the ascidian Ciona intestinalis, we suggest new evolutionary histories for different TK families and draw a new global picture of gene loss/gain in the different phyla. Surprisingly, our survey also detected an unprecedented expansion of a group of closely related TK families, including TIE, FGFR, PDGFR, and RET, due most probably to massive gene duplication and exon shuffling. Based on their highly similar intron/exon structure at the TK domain, we suggest that this group of TK families constitute a superfamily of TK proteins, which we termed EXpanding TK, after their seemingly unique propensity to gene duplication and exon shuffling, not only in amphioxus but also across all metazoan groups. Due to this extreme tendency to both retention and expansion of TK genes, amphioxus harbors the richest and most diverse TK repertoire among all metazoans studied so far, retaining most of the gene complement of its ancestors, but having evolved its own repertoire of genetic novelties.

Plant receptor-like serine threonine kinases: roles in signaling and plant defense

Plants are hosts to a wide array of pathogens from all kingdoms of life. In the absence of an active immune system or combinatorial diversifications that lead to recombination-driven somatic gene flexibility, plants have evolved different strategies to combat both individual pathogen strains and changing pathogen populations. The receptor-like kinase (RLK) gene-family expansion in plants was hypothesized to have allowed accelerated evolution among domains implicated in signal reception, typically a leucine-rich repeat (LRR). Under that model, the gene-family expansion represents a plant-specific adaptation that leads to the production of numerous and variable cell surface and cytoplasmic receptors. More recently, it has emerged that the LRR domains of RLK interact with a diverse group of proteins leading to combinatorial variations in signal response specificity. Therefore, the RLK appear to play a central role in signaling during pathogen recognition, the subsequent activation of plant defense mechanisms, and developmental control. The future challenges will include determinations of RLK modes of action, the basis of recognition and specificity, which cellular responses each receptor mediates, and how both receptor and kinase domain interactions fit into the defense signaling cascades. These challenges will be complicated by the limited information that may be derived from the primary sequence of the LRR domain. The review focuses upon implications derived from recent studies of the secondary and tertiary structures of several plant RLK that change understanding of plant receptor function and signaling. In addition, the biological functions of plant and animal RLK-containing receptors were reviewed and commonalities among their signaling mechanisms identified. Further elucidated were the genomic and structural organizations of RLK gene families, with special emphasis on RLK implicated in resistance to disease and development.

Evolvable signaling networks of receptor tyrosine kinases: relevance of robustness to malignancy and to cancer therapy

Robust biological signaling networks evolved, through gene duplications, from simple, relatively fragile cascades. Architectural features such as layered configuration, branching and modularity, as well as functional characteristics (e.g., feedback control circuits), enable fail-safe performance in the face of internal and external perturbations. These universal features are exemplified here using the receptor tyrosine kinase (RTK) family. The RTK module is richly mutated and overexpressed in human malignancies, and pharmaceutical interception of its signaling effectively retards growth of specific tumors. Therapy-induced interception of RTK-signaling pathways and the common evolvement of drug resistance are respectively considered here as manifestations of fragility and plasticity of robust networks. The systems perspective we present views pathologies as hijackers of biological robustness and offers ways for identifying fragile hubs, as well as strategies to overcome drug resistance.

Large-scale identification and evolution indexing of tyrosine phosphorylation sites from murine brain

Metazoans employ reversible tyrosine phosphorylation to regulate innumerable biological processes. Thus, the large-scale identification of tyrosine phosphorylation sites from primary tissues is an essential step toward a molecular systems understanding of dynamic regulation in vivo. The relative paucity of phosphotyrosine has greatly limited its identification in large-scale phosphoproteomic experiments. However, using antiphosphotyrosine peptide immunoprecipitations, we report the largest study to date of tyrosine phosphorylation sites from primary tissue, identifying 414 unique tyrosine phosphorylation sites from murine brain. To measure the conservation of phosphorylated tyrosines and their surrounding residues, we constructed a computational pipeline and identified patterns of conservation within the signature of phosphotyrosine.

Classification and functional annotation of eukaryotic protein kinases

Reversible protein phosphorylation by protein kinases and phosphatases is a ubiquitous signaling mechanism in all eukaryotic cells. A multilevel hidden Markov model library is presented which is able to classify protein kinases into one of 12 families, with a misclassification rate of zero on the characterized kinomes of H. sapiens, M. musculus, D. melanogaster, C. elegans, S. cerevisiae, D. discoideum, and P. falciparum. The Library is shown to outperform BLASTP and a general Pfam hidden Markov model of the kinase catalytic domain in the retrieval and family-level classification of protein kinases. The application of the Library to the 38 unclassified kinases of yeast enriches the yeast kinome in protein kinases of the families AGC (5), CAMK (17), CMGC (4), and STE (1), thereby raising the family-level classification of yeast conventional protein kinases from 66.96 to 90.43%. The application of the Library to 21 eukaryotic genomes shows seven families (AGC, CAMK, CK1, CMGC, STE, PIKK, and RIO) to be present in all genomes analyzed, and so is likely to be essential to eukaryotes. Putative tyrosine kinases (TKs) are found in the plants A. thaliana (2), O. sativa ssp. Indica (6), and O. sativa ssp. Japonica (7), and in the amoeba E. histolytica (7). To our knowledge, TKs have not been predicted in plants before. This also suggests that a primitive set of TKs might have predated the radiation of eukaryotes. Putative tyrosine kinase-like kinases (TKLs) are found in the fungi C. neoformans (2), P. chrysosporium (4), in the Apicomplexans C. hominis (4), P. yoelii (4), and P. falciparum (6), the amoeba E. histolytica (109), and the alga T. pseudonana (6). TKLs are found to be abundant in plants (776 in A. thaliana, 1010 in O. sativa ssp. Indica, and 969 in O. sativa ssp. Japonica). TKLs might have predated the radiation of eukaryotes too and have been lost secondarily from some fungi. The application of the Library facilitates the annotation of kinomes and has provided novel insights on the early evolution and subsequent adaptations of the various protein kinase families in eukaryotes.

Domain architecture comparison for multidomain homology identification

Homology identification is the first step for many genomic studies. Current methods, based on sequence comparison, can result in a substantial number of mis-assignments due to the similarity of homologous domains in otherwise unrelated sequences. Here we propose methods to detect homologs through explicit comparison of protein domain content. We developed several schemes for scoring the homology of a pair of protein sequences based on methods used in the field of information retrieval. We evaluate the proposed methods and methods used in the literature using a benchmark of fifteen sequence families of known evolutionary history. The results of these studies demonstrate the effectiveness of comparing domain architectures using these similarity measures. We also demonstrate the importance of both weighting promiscuous domains and of compensating for the statistical effect of having a large number of domains in a protein. Using logistic regression, we demonstrate the benefit of combining similarity measures based on domain content with sequence similarity measures.

Circadian rhythms lit up in Chlamydomonas

Recent work on the circadian clock of the unicellular green alga Chlamydomonas reinhardtii strengthens its standing as a convenient model system for circadian study. It was shown to be amenable to molecular engineering using a luciferase-based real-time reporter for circadian rhythms. Together with the completed draft genomic sequence, the new system opens the door for genome-scale forward and reverse genetic analysis.

Comparative analysis of the kinomes of three pathogenic trypanosomatids: Leishmania major, Trypanosoma brucei and Trypanosoma cruzi

Background

The trypanosomatids Leishmania major, Trypanosoma brucei and Trypanosoma cruzi cause some of the most debilitating diseases of humankind: cutaneous leishmaniasis, African sleeping sickness, and Chagas disease. These protozoa possess complex life cycles that involve development in mammalian and insect hosts, and a tightly coordinated cell cycle ensures propagation of the highly polarized cells. However, the ways in which the parasites respond to their environment and coordinate intracellular processes are poorly understood. As a part of an effort to understand parasite signaling functions, we report the results of a genome-wide analysis of protein kinases (PKs) of these three trypanosomatids.

Results

Bioinformatic searches of the trypanosomatid genomes for eukaryotic PKs (ePKs) and atypical PKs (aPKs) revealed a total of 176 PKs in T. brucei, 190 in T. cruzi and 199 in L. major, most of which are orthologous across the three species. This is approximately 30% of the number in the human host and double that of the malaria parasite, Plasmodium falciparum. The representation of various groups of ePKs differs significantly as compared to humans: trypanosomatids lack receptor-linked tyrosine and tyrosine kinase-like kinases, although they do possess dual-specificity kinases. A relative expansion of the CMGC, STE and NEK groups has occurred. A large number of unique ePKs show no strong affinity to any known group. The trypanosomatids possess few ePKs with predicted transmembrane domains, suggesting that receptor ePKs are rare. Accessory Pfam domains, which are frequently present in human ePKs, are uncommon in trypanosomatid ePKs.

Conclusion

Trypanosomatids possess a large set of PKs, comprising approximately 2% of each genome, suggesting a key role for phosphorylation in parasite biology. Whilst it was possible to place most of the trypanosomatid ePKs into the seven established groups using bioinformatic analyses, it has not been possible to ascribe function based solely on sequence similarity. Hence the connection of stimuli to protein phosphorylation networks remains enigmatic. The presence of numerous PKs with significant sequence similarity to known drug targets, as well as a large number of unusual kinases that might represent novel targets, strongly argue for functional analysis of these molecules.

Identification of a gene in the process of being lost from the genus Agrostis

Lineage-specific gene loss is considered one of the processes contributing to speciation and genome diversity. Such gene loss has been inferred from interspecies comparisons of orthologous DNA segments. Examples of intraspecific gene loss are rare. Here we report identification of a gene, designated Crs-1 (creeping specific-1), that appears to be in the process of being lost from heterozygous populations of the species creeping bentgrass (Agrostis stolonifera). The Crs-1 gene encodes a protein with an N-terminal dirigent protein domain and a C-terminal lectin domain and is similar to the maize (Zea mays) beta-glucosidase aggregating factor. Most individual creeping bentgrass plants examined are lacking Crs-1. Some individuals are hemizygous for the Crs-1 locus, indicating major haplotype noncolinearity at that locus. Crs-1 was not detected in several other Agrostis species, indicating it is being lost from the genus. The Crs-1 locus in creeping bentgrass provides a rare example of the evolutionary process of gene loss occurring within a plant species.

The evolution of Ca2+ signalling in photosynthetic eukaryotes

It is likely that cytosolic Ca2+ elevations have played a part in eukaryotic signal transduction for about the last 2 Gyr, being mediated by a group of molecules which are collectively known as the [Ca2+]cyt signalling toolkit. Different eukaryotes often display strikingly similar [Ca2+]cyt signalling elevations, which may reflect conservation of toolkit components (homology) or similar constraints acting on different toolkits (homoplasy). Certain toolkit components, which are presumably ancestral, are shared by plants and animals, but some components are unique to photosynthetic organisms. We propose that the structure of modern plant [Ca2+]cyt signalling toolkits may be explained by their modular adaptation from earlier pathways.

Src proteins/src genes: from sponges to mammals

The genome of marine sponge Suberites domuncula, a member of the most ancient and most simple metazoan phylum Porifera, encodes at least five genes for Src-type proteins, more than, i.e., Caenorhabditis elegans or Drosophila melanogaster (two in each). Three proteins, SRC1SD, SRC2SD and SRC3SD, were fully characterized. The overall homology (identity+similarity) among the three S. domuncula Srcs (68-71%) is much lower than the sequence conservation between orthologous Src proteins from freshwater sponges (82-85%). It is therefore very likely that several src genes/proteins were already present in the genome of Urmetazoa, the hypothetical metazoan ancestor. We have identified in the S. domuncula expressed sequence tags (ESTs) database further Src homology 2 (SH2) and 3 (SH3) domains that are unrelated to protein tyrosine kinases (PTKs). Src-related SH2 and SH3 domains from different species are much more conserved than SH2 and SH3 domains from different proteins in the same organism (S. domuncula), supporting the view that the common, ancestral src gene was already a multidomain protein composed of SH3, SH2 and tyrosine kinase (TK) domains. Two S. domuncula src genes were fully sequenced: src1SD gene has six and src2SD gene only one intron in front of SH2 domain, located at the same position in both genes. All vertebrate src genes, from fish to human, originated from the same ancestral gene, because they all have 10 introns at conserved positions. However, src genes in invertebrates have fewer introns that are located at different positions. Only the intron in front of the SH2 domain is present at the absolutely conserved position (and phase) in all known src genes, indicating that at least this intron was already present in the ancestral gene, common to all Metazoa. Our results also suggest that TK domain in this ancestral src was encoded on a single exon.