Structural and evolutionary adaptation of rhoptry kinases and pseudokinases, a family of coccidian virulence factors

BACKGROUND: The widespread protozoan parasite Toxoplasma gondii interferes with host cell functions by exporting the contents of a unique apical organelle, the rhoptry. Among the mix of secreted proteins are an expanded, lineage-specific family of protein kinases termed rhoptry kinases (ROPKs), several of which have been shown to be key virulence factors, including the pseudokinase ROP5. The extent and details of the diversification of this protein family are poorly understood. RESULTS: In this study, we comprehensively catalogued the ROPK family in the genomes of Toxoplasma gondii, Neospora caninum and Eimeria tenella, as well as portions of the unfinished genome of Sarcocystis neurona, and classified the identified genes into 42 distinct subfamilies. We systematically compared the rhoptry kinase protein sequences and structures to each other and to the broader superfamily of eukaryotic protein kinases to study the patterns of diversification and neofunctionalization in the ROPK family and its subfamilies. We identified three ROPK sub-clades of particular interest: those bearing a structurally conserved N-terminal extension to the kinase domain (NTE), an E. tenella-specific expansion, and a basal cluster including ROP35 and BPK1 that we term ROPKL. Structural analysis in light of the solved structures ROP2, ROP5, ROP8 and in comparison to typical eukaryotic protein kinases revealed ROPK-specific conservation patterns in two key regions of the kinase domain, surrounding a ROPK-conserved insert in the kinase hinge region and a disulfide bridge in the kinase substrate-binding lobe. We also examined conservation patterns specific to the NTE-bearing clade. We discuss the possible functional consequences of each. CONCLUSIONS: Our work sheds light on several important but previously unrecognized features shared among rhoptry kinases, as well as the essential differences between active and degenerate protein kinases. We identify the most distinctive ROPK-specific features conserved across both active kinases and pseudokinases, and discuss these in terms of sequence motifs, evolutionary context, structural impact and potential functional relevance. By characterizing the proteins that enable these parasites to invade the host cell and co-opt its signaling mechanisms, we provide guidance on potential therapeutic targets for the diseases caused by coccidian parasites

Bio.Phylo: A Unified Toolkit for Processing, Analyzing and Visualizing Phylogenetic Trees in Biopython

Background: Ongoing innovation in phylogenetics and evolutionary biology has been accompanied by a proliferation of software tools, data formats, analytical techniques and web servers. This brings with it the challenge of integrating phylogenetic and other related biological data found in a wide variety of formats, and underlines the need for reusable software that can read, manipulate and transform this information into the various forms required to build computational pipelines. Results: We built a Python software library for working with phylogenetic data that is tightly integrated with Biopython, a broad-ranging toolkit for computational biology. Our library, Bio.Phylo, is highly interoperable with existing libraries, tools and standards, and is capable of parsing common file formats for phylogenetic trees, performing basic transformations and manipulations, attaching rich annotations, and visualizing trees. We unified the modules for working with the standard file formats Newick, NEXUS and phyloXML behind a consistent and simple API, providing a common set of functionality independent of the data source. Conclusions: Bio.Phylo meets a growing need in bioinformatics for working with heterogeneous types of phylogenetic data. By supporting interoperability with multiple file formats and leveraging existing Biopython features, this library simplifies the construction of phylogenetic workflows. We also provide examples of the benefits of building a community around a shared open-source project. Bio.Phylo is included with Biopython, available through the Biopython website, http://biopython.org

An evolutionary perspective on the kinome of malaria parasites

Malaria parasites belong to an ancient lineage that diverged very early from the main branch of eukaryotes. The approximately 90-member plasmodial kinome includes a majority of eukaryotic protein kinases that clearly cluster within the AGC, CMGC, TKL, CaMK and CK1 groups found in yeast, plants and mammals, testifying to the ancient ancestry of these families. However, several hundred millions years of independent evolution, and the specific pressures brought about by first a photosynthetic and then a parasitic lifestyle, led to the emergence of unique features in the plasmodial kinome. These include taxon-restricted kinase families, and unique peculiarities of individual enzymes even when they have homologues in other eukaryotes. Here, we merge essential aspects of all three malaria-related communications that were presented at the Evolution of Protein Phosphorylation meeting, and propose an integrated discussion of the specific features of the parasite's kinome and phosphoproteome

Genomic profiling of malignant peritoneal mesothelioma reveals recurrent alterations in epigenetic regulatory genes BAP1, SETD2, and DDX3X.

Malignant mesothelioma is a rare cancer that arises from the mesothelial cells that line the pleural cavity and less commonly from the peritoneal lining of the abdomen and pelvis. Most pleural mesotheliomas arise in patients with a history of asbestos exposure, whereas the association of peritoneal mesotheliomas with exposure to asbestos and other potential carcinogens is less clear, suggesting that the genetic alterations that drive malignant peritoneal mesothelioma may be unique from those in pleural mesothelioma. Treatment options for all malignant mesotheliomas are currently limited, with no known targeted therapies available. To better understand the molecular pathogenesis of malignant peritoneal mesothelioma, we sequenced 510 cancer-related genes in 13 patients with malignant mesothelioma arising in the peritoneal cavity. The most frequent genetic alteration was biallelic inactivation of the BAP1 gene, which occurred in 9/13 cases, with an additional two cases demonstrating monoallelic loss of BAP1. All 11 of these cases demonstrated loss of BAP1 nuclear staining by immunohistochemistry, whereas two tumors without BAP1 alteration and all 42 cases of histologic mimics in peritoneum (8 multilocular peritoneal inclusion cyst, 6 well-differentiated papillary mesothelioma of the peritoneum, 16 adenomatoid tumor, and 12 low-grade serous carcinoma of the ovary) demonstrated intact BAP1 nuclear staining. Additional recurrently mutated genes in this cohort of malignant peritoneal mesotheliomas included NF2 (3/13), SETD2 (2/13), and DDX3X (2/13). While these genes are known to be recurrently mutated in pleural mesotheliomas, the frequencies are distinct in peritoneal mesotheliomas, with nearly 85% of peritoneal tumors harboring BAP1 alterations versus only 20-30% of pleural tumors. Together, these findings demonstrate the importance of epigenetic modifiers including BAP1, SETD2, and DDX3X in mesothelial tumorigenesis and suggest opportunities for targeted therapies

Adenomatoid tumors of the male and female genital tract are defined by TRAF7 mutations that drive aberrant NF-kB pathway activation.

Adenomatoid tumors are the most common neoplasm of the epididymis, and histologically similar adenomatoid tumors also commonly arise in the uterus and fallopian tube. To investigate the molecular pathogenesis of these tumors, we performed genomic profiling on a cohort of 31 adenomatoid tumors of the male and female genital tracts. We identified that all tumors harbored somatic missense mutations in the TRAF7 gene, which encodes an E3 ubiquitin ligase belonging to the family of tumor necrosis factor receptor-associated factors (TRAFs). These mutations all clustered into one of five recurrent hotspots within the WD40 repeat domains at the C-terminus of the protein. Functional studies in vitro revealed that expression of mutant but not wild-type TRAF7 led to increased phosphorylation of nuclear factor-kappa B (NF-kB) and increased expression of L1 cell adhesion molecule (L1CAM), a marker of NF-kB pathway activation. Immunohistochemistry demonstrated robust L1CAM expression in adenomatoid tumors that was absent in normal mesothelial cells, malignant peritoneal mesotheliomas and multilocular peritoneal inclusion cysts. Together, these studies demonstrate that adenomatoid tumors of the male and female genital tract are genetically defined by TRAF7 mutation that drives aberrant NF-kB pathway activation

Co-Conserved Features Associated with cis Regulation of ErbB Tyrosine Kinases

BACKGROUND: The epidermal growth factor receptor kinases, or ErbB kinases, belong to a large sub-group of receptor tyrosine kinases (RTKs), which share a conserved catalytic core. The catalytic core of ErbB kinases have functionally diverged from other RTKs in that they are activated by a unique allosteric mechanism that involves specific interactions between the kinase core and the flanking Juxtamembrane (JM) and COOH-terminal tail (C-terminal tail). Although extensive studies on ErbB and related tyrosine kinases have provided important insights into the structural basis for ErbB kinase functional divergence, the sequence features that contribute to the unique regulation of ErbB kinases have not been systematically explored. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we use a Bayesian approach to identify the selective sequence constraints that most distinguish ErbB kinases from other receptor tyrosine kinases. We find that strong ErbB kinase-specific constraints are imposed on residues that tether the JM and C-terminal tail to key functional regions of the kinase core. A conserved RIxKExE motif in the JM-kinase linker region and a glutamine in the inter-lobe linker are identified as two of the most distinguishing features of the ErbB family. While the RIxKExE motif tethers the C-terminal tail to the N-lobe of the kinase domain, the glutamine tethers the C-terminal tail to hinge regions critical for inter-lobe movement. Comparison of the active and inactive crystal structures of ErbB kinases indicates that the identified residues are conformationally malleable and can potentially contribute to the cis regulation of the kinase core by the JM and C-terminal tail. ErbB3, and EGFR orthologs in sponges and parasitic worms, diverge from some of the canonical ErbB features, providing insights into sub-family and lineage-specific functional specialization. CONCLUSION/SIGNIFICANCE: Our analysis pinpoints key residues for mutational analysis, and provides new clues to cancer mutations that alter the canonical modes of ErbB kinase regulation

Structural and evolutionary divergence of eukaryotic protein kinases in Apicomplexa

Abstract Background The Apicomplexa constitute an evolutionarily divergent phylum of protozoan pathogens responsible for widespread parasitic diseases such as malaria and toxoplasmosis. Many cellular functions in these medically important organisms are controlled by protein kinases, which have emerged as promising drug targets for parasitic diseases. However, an incomplete understanding of how apicomplexan kinases structurally and mechanistically differ from their host counterparts has hindered drug development efforts to target parasite kinases. Results We used the wealth of sequence data recently made available for 15 apicomplexan species to identify the kinome of each species and quantify the evolutionary constraints imposed on each family of apicomplexan kinases. Our analysis revealed lineage-specific adaptations in selected families, namely cyclin-dependent kinase (CDK), calcium-dependent protein kinase (CDPK) and CLK/LAMMER, which have been identified as important in the pathogenesis of these organisms. Bayesian analysis of selective constraints imposed on these families identified the sequence and structural features that most distinguish apicomplexan protein kinases from their homologs in model organisms and other eukaryotes. In particular, in a subfamily of CDKs orthologous to Plasmodium falciparum crk-5, the activation loop contains a novel PTxC motif which is absent from all CDKs outside Apicomplexa. Our analysis also suggests a convergent mode of regulation in a subset of apicomplexan CDPKs and mammalian MAPKs involving a commonly conserved arginine in the αC helix. In all recognized apicomplexan CLKs, we find a set of co-conserved residues involved in substrate recognition and docking that are distinct from metazoan CLKs. Conclusions We pinpoint key conserved residues that can be predicted to mediate functional differences from eukaryotic homologs in three identified kinase families. We discuss the structural, functional and evolutionary implications of these lineage-specific variations and propose specific hypotheses for experimental investigation. The apicomplexan-specific kinase features reported in this study can be used in the design of selective kinase inhibitors.</p

CNVkit: Genome-Wide Copy Number Detection and Visualization from Targeted DNA Sequencing.

Germline copy number variants (CNVs) and somatic copy number alterations (SCNAs) are of significant importance in syndromic conditions and cancer. Massively parallel sequencing is increasingly used to infer copy number information from variations in the read depth in sequencing data. However, this approach has limitations in the case of targeted re-sequencing, which leaves gaps in coverage between the regions chosen for enrichment and introduces biases related to the efficiency of target capture and library preparation. We present a method for copy number detection, implemented in the software package CNVkit, that uses both the targeted reads and the nonspecifically captured off-target reads to infer copy number evenly across the genome. This combination achieves both exon-level resolution in targeted regions and sufficient resolution in the larger intronic and intergenic regions to identify copy number changes. In particular, we successfully inferred copy number at equivalent to 100-kilobase resolution genome-wide from a platform targeting as few as 293 genes. After normalizing read counts to a pooled reference, we evaluated and corrected for three sources of bias that explain most of the extraneous variability in the sequencing read depth: GC content, target footprint size and spacing, and repetitive sequences. We compared the performance of CNVkit to copy number changes identified by array comparative genomic hybridization. We packaged the components of CNVkit so that it is straightforward to use and provides visualizations, detailed reporting of significant features, and export options for integration into existing analysis pipelines. CNVkit is freely available from https://github.com/etal/cnvkit