Search

Plotting a) the median spread of coverage per-gene, from the pooled reference data and b) the per-gene standard deviations for duplicated samples shows that some genes are less reliable than others.

Genomic regions used for the targeted panel sequencing in this study.

Purpose: Precise detection of copy number aberrations (CNA) from tumor biopsies is critically important to the treatment of metastatic prostate cancer. The use of targeted panel next-generation sequencing (NGS) is inexpensive, high throughput, and easily feasible, allowing single-nucleotide variant calls, but CNA estimation from this remains challenging.Experimental Design: We evaluated CNVkit for CNA identification from amplicon-based targeted NGS in a cohort of 110 fresh castration-resistant prostate cancer biopsies and used capture-based whole-exome sequencing (WES), array comparative genomic hybridization (aCGH), and FISH to explore the viability of this approach.Results: We showed that this method produced highly reproducible CNA results (r = 0.92), with the use of pooled germline DNA as a coverage reference supporting precise CNA estimation. CNA estimates from targeted NGS were comparable with WES (r = 0.86) and aCGH (r = 0.7); for key selected genes (BRCA2, MYC, PIK3CA, PTEN, and RB1), CNA estimation correlated well with WES (r = 0.91) and aCGH (r = 0.84) results. The frequency of CNAs in our population was comparable with that previously described (i.e., deep deletions: BRCA2 4.5%; RB1 8.2%; PTEN 15.5%; amplification: AR 45.5%; gain: MYC 31.8%). We also showed, utilizing FISH, that CNA estimation can be impacted by intratumor heterogeneity and demonstrated that tumor microdissection allows NGS to provide more precise CNA estimates.Conclusions: Targeted NGS and CNVkit-based analyses provide a robust, precise, high-throughput, and cost-effective method for CNA estimation for the delivery of more precise patient care. Clin Cancer Res; 23(20); 6070–7. ©2017 AACR.

Supplementary legend

Importance Cerebrospinal fluid (CSF) cytologic testing and flow cytometry are insensitive for diagnosing neoplasms of the central nervous system (CNS). Such clinical phenotypes can mimic infectious and autoimmune causes of meningoencephalitis. Objective To ascertain whether CSF metagenomic next-generation sequencing (mNGS) can identify aneuploidy, a hallmark of malignant neoplasms, in difficult-to-diagnose cases of CNS malignant neoplasm. Design, setting, and participants Two case-control studies were performed at the University of California, San Francisco (UCSF). The first study used CSF specimens collected at the UCSF Clinical Laboratories between July 1, 2017, and December 31, 2019, and evaluated test performance in specimens from patients with a CNS malignant neoplasm (positive controls) or without (negative controls). The results were compared with those from CSF cytologic testing and/or flow cytometry. The second study evaluated patients who were enrolled in an ongoing prospective study between April 1, 2014, and July 31, 2019, with presentations that were suggestive of neuroinflammatory disease but who were ultimately diagnosed with a CNS malignant neoplasm. Cases of individuals whose tumors could have been detected earlier without additional invasive testing are discussed. Main outcomes and measures The primary outcome measures were the sensitivity and specificity of aneuploidy detection by CSF mNGS. Secondary subset analyses included a comparison of CSF and tumor tissue chromosomal abnormalities and the identification of neuroimaging characteristics that were associated with test performance. Results Across both studies, 130 participants were included (median [interquartile range] age, 57.5 [43.3-68.0] years; 72 men [55.4%]). The test performance study used 125 residual laboratory CSF specimens from 47 patients with a CNS malignant neoplasm and 56 patients with other neurological diseases. The neuroinflammatory disease study enrolled 12 patients and 17 matched control participants. The sensitivity of the CSF mNGS assay was 75% (95% CI, 63%-85%), and the specificity was 100% (95% CI, 96%-100%). Aneuploidy was detected in 64% (95% CI, 41%-83%) of the patients in the test performance study with nondiagnostic cytologic testing and/or flow cytometry, and in 55% (95% CI, 23%-83%) of patients in the neuroinflammatory disease study who were ultimately diagnosed with a CNS malignant neoplasm. Of the patients in whom aneuploidy was detected, 38 (90.5%) had multiple copy number variations with tumor fractions ranging from 31% to 49%. Conclusions and relevance This case-control study showed that CSF mNGS, which has low specimen volume requirements, does not require the preservation of cell integrity, and was orginally developed to diagnose neurologic infections, can also detect genetic evidence of a CNS malignant neoplasm in patients in whom CSF cytologic testing and/or flow cytometry yielded negative results with a low risk of false-positive results.

Summary: BRAF fusions are detected in numerous neoplasms, but their clinical management remains unresolved. We identified six melanoma lines harboring BRAF fusions representative of the clinical cases reported in the literature. Their unexpected heterogeneous responses to RAF and MEK inhibitors could be categorized upon specific features of the fusion kinases. Higher expression level correlated with resistance, and fusion partners containing a dimerization domain promoted paradoxical activation of the mitogen-activated protein kinase (MAPK) pathway and hyperproliferation in response to first- and second-generation RAF inhibitors. By contrast, next-generation αC-IN/DFG-OUT RAF inhibitors blunted paradoxical activation across all lines and had their therapeutic efficacy further increased in vitro and in vivo by combination with MEK inhibitors, opening perspectives in the clinical management of tumors harboring BRAF fusions. : Botton et al. shed light on the heterogeneity of BRAF fusions encountered in melanocytic tumors and characterize features influencing their signaling and drug response. These findings unveil the singularities of BRAF fusions and establish general principles to guide their clinical management in melanoma and other malignancies. Keywords: BRAF fusion, melanoma, paradoxical activation, RAF inhibitor, MEK inhibitor, sequencing, rearrangement, translocation, kinase, pre-clinical

Combined hepatocellular-cholangiocarcinomas (CHC) are mixed tumours with both hepatocellular carcinoma (HCC) and cholangiocarcinoma (CC) components. CHC prognosis is similar to intrahepatic CC (ICC) and worse than HCC; staging and treatment generally follow ICC algorithms. However, the molecular biology of CHC remains poorly characterised. We performed capture-based next-generation sequencing of 20 CHC and, for comparison, 10 ICC arising in cirrhosis. Intratumour heterogeneity was assessed by separately sequencing the HCC and CC components of nine CHC. CHC demonstrated molecular profiles similar to HCC, even in the CC component. CHC harboured recurrent alterations in TERT (80%), TP53 (80%), cell cycle genes (40%; CCND1, CCNE1, CDKN2A), receptor tyrosine kinase/Ras/PI3-kinase pathway genes (55%; MET, ERBB2, KRAS, PTEN), chromatin regulators (20%; ARID1A, ARID2) and Wnt pathway genes (20%; CTNNB1, AXIN, APC). No CHC had alterations in IDH1, IDH2, FGFR2 or BAP1, genes typically mutated in ICC. TERT promoter mutations were consistently identified in both HCC and CC components, supporting TERT alteration as an early event in CHC evolution. TP53 mutations were present in both components in slightly over half the TP53-altered cases. By contrast, focal amplifications of CCND1, MET and ERRB2, as well as Wnt pathway alterations, were most often exclusive to one component, suggesting that these are late events in CHC evolution. ICC in cirrhosis demonstrated alterations similar to ICC in non-cirrhotic liver, including in IDH1 or IDH2 (30%), CDKN2A (40%), FGFR2 (20%), PBRM1 (20%), ARID1A (10%) and BAP1 (10%). TERT promoter and TP53 mutation were present in only one ICC each. Our data demonstrate that CHC genetics are distinct from ICC (even in cirrhosis) and similar to HCC, which has diagnostic utility and implications for treatment. Copyright © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

ImportanceMalignant neoplasms of the central nervous system (CNS) are frequently not detected by cerebrospinal fluid (CSF) flow cytometry or cytology, and clinical phenotypes can overlap with inflammatory meningoencephalitis.ObjectiveTo determine whether an existing CSF metagenomic next-generation sequencing (mNGS) assay can identify a hallmark of malignant neoplasms — aneuploidy — in difficult-to-diagnose cases of CNS malignancy.DesignTwo retrospective, case-control studies included a total of 155 samples from patients with an eventual diagnosis of a CNS malignancy (n=59 patients) and controls with other CNS diseases (n=73 patients). The first study was used to evaluate test performance in positive and negative controls. The second study was used to assess the potential utility of aneuploidy detection in patients whose CSF was sent for mNGS because of suspected neuroinflammatory disease who were ultimately found to have a CNS malignancy.SettingThis is a single site study at a large tertiary care center, University of California San Francisco, that enrolled from 2014 to 2019.ParticipantsThe test performance case-control study enrolled positive control patients with a CNS malignancy (n=47 patients) and negative controls with other neurologic diseases (n=56 patients) who had had CSF flow cytometry and/or cytology performed. The second case-control study enrolled patients with suspected neuroinflammatory disease who were ultimately diagnosed with a CNS malignancy (n=12) and other neurologic disease controls (n=17).Main Outcome(s) and Measure(s)The primary outcome measures were the performance characteristics of detecting aneuploidy in CSF by a cell-free DNA mNGS assay compared to cytology and/or flow cytometry and the tumor fraction in CSF from patients with CNS malignancies.ResultsAcross the two case-control studies, the overall sensitivity of the CSF mNGS assay for detecting aneuploidy in patients ultimately diagnosed with a CNS malignancy was 75% (63-96%, 95% CI), and specificity was 100% (96-100%, 95% CI). Notably, CSF mNGS detected aneuploidy in 64% of the non-diagnostic cytology and flow cytometry cases in the test performance study and in 55% of the cases with suspected neuroinflammatory disease who were ultimately diagnosed with a CNS malignancy. Of the cases in whom aneuploidy was detected, 90% had multiple chromosomal copy number variants with tumor fractions ranging from 31% to 49%.Conclusions and RelevanceMetagenomic NGS of CSF, originally designed to diagnose neurologic infections, detects evidence of CNS malignancies (i.e., aneuploidy) in cases where CSF flow cytometry and/or cytology were negative with a low risk of false positive results.3 Key PointsQuestionCan CSF metagenomic NGS, a test designed to diagnose infections, also detect genetic signatures of cancer in patients with suspected neuroinflammatory disease?FindingsAcross two case-control studies of patients with negative CSF cytology and/or flow cytometry, CSF mNGS detected genetic evidence for a malignancy with a sensitivity of 75% (63-85%, 95% CI) and specificity of 100% (96-100%, 95% CI).MeaningCSF mNGS, an assay with low sample volume requirements that does not require the preservation of cell integrity, has the potential to complement cancer detection by CSF flow cytometry and cytology.

Background Metagenomic next-generation sequencing (mNGS) of body fluids is an emerging approach to identify occult pathogens in undiagnosed patients. We hypothesized that metagenomic testing can be simultaneously used to detect malignant neoplasms in addition to infectious pathogens. Methods From two independent studies (n = 205), we used human data generated from a metagenomic sequencing pipeline to simultaneously screen for malignancies by copy number variation (CNV) detection. In the first case-control study, we analyzed body fluid samples (n = 124) from patients with a clinical diagnosis of either malignancy (positive cases, n = 65) or infection (negative controls, n = 59). In a second verification cohort, we analyzed a series of consecutive cases (n = 81) sent to cytology for malignancy workup that included malignant positives (n = 32), negatives (n = 18), or cases with an unclear gold standard (n = 31). Results The overall CNV test sensitivity across all studies was 87% (55 of 63) in patients with malignancies confirmed by conventional cytology and/or flow cytometry testing and 68% (23 of 34) in patients who were ultimately diagnosed with cancer but negative by conventional testing. Specificity was 100% (95% CI 95–100%) with no false positives detected in 77 negative controls. In one example, a patient hospitalized with an unknown pulmonary illness had non-diagnostic lung biopsies, while CNVs implicating a malignancy were detectable from bronchoalveolar fluid. Conclusions Metagenomic sequencing of body fluids can be used to identify undetected malignant neoplasms through copy number variation detection. This study illustrates the potential clinical utility of a single metagenomic test to uncover the cause of undiagnosed acute illnesses due to cancer or infection using the same specimen.

The deubiquitinating enzyme BAP1 is mutated in a hereditary cancer syndrome with a high risk of mesothelioma and melanocytic tumors. Here, we show that Bap1 deletion in melanocytes cooperates with the constitutively active, oncogenic form of BRAF (BRAFV600E ) and UV to cause melanoma in mice, albeit at very low frequency. In addition, Bap1-null melanoma cells derived from mouse tumors are more aggressive and colonize and grow at distant sites more than their wild-type counterparts. Molecularly, Bap1-null melanoma cell lines have increased DNA damage measured by γH2aX and hyperubiquitination of histone H2a. Therapeutically, these Bap1-null tumors are completely responsive to BRAF- and MEK-targeted therapies. Therefore, BAP1 functions as a tumor suppressor and limits tumor progression in melanoma.

Chordoid glioma is a rare brain tumor thought to arise from specialized glial cells of the lamina terminalis along the anterior wall of the third ventricle. Despite being histologically low-grade, chordoid gliomas are often associated with poor outcome, as their stereotypic location in the third ventricle makes resection challenging and efficacious adjuvant therapies have not been developed. Here we performed genomic profiling on 13 chordoid gliomas and identified a recurrent D463H missense mutation in PRKCA in all tumors, which localizes in the kinase domain of the encoded protein kinase C alpha (PKCα). Expression of mutant PRKCA in immortalized human astrocytes led to increased phospho-ERK and anchorage-independent growth that could be blocked by MEK inhibition. These studies define PRKCA as a recurrently mutated oncogene in human cancer and identify a potential therapeutic vulnerability in this uncommon brain tumor., Chordoid glioma is a rare low-grade brain tumor that originates from the anterior wall of the third ventricle where surgical resection is challenging; the clinical outcome of patients after subtotal resection or disease recurrence is poor. Here the authors identify a recurrent missense mutation in PRKCA that may serve as a potential therapeutic target in this uncommon brain cancer.

Malignant phyllodes tumors of the breast are poorly understood rare neoplasms with potential for aggressive behavior. Few efficacious treatment options exist for progressed or metastatic disease. The molecular features of malignant phyllodes tumors are poorly defined, and a deeper understanding of the genetics of these tumors may shed light on pathogenesis and progression and potentially identify novel treatment approaches. We sequenced 510 cancer-related genes in 10 malignant phyllodes tumors, including 5 tumors with liposarcomatous differentiation and 1 with myxoid chondrosarcoma-like differentiation. Intratumoral heterogeneity was assessed by sequencing two separate areas in 7 tumors, including non-heterologous and heterologous components of tumors with heterologous differentiation. Activating hotspot mutations in FGFR1 were identified in 2 tumors. Additional recurrently mutated genes included TERT promoter (6/10), TP53 (4/10), PIK3CA (3/10), MED12 (3/10), SETD2 (2/10) and KMT2D (2/10). Together, genomic aberrations in FGFR/EGFR PI-3 kinase and RAS pathways were identified in 8 (80%) tumors and included mutually exclusive and potentially actionable activating FGFR1, PIK3CA and BRAF V600E mutations, inactivating TSC2 mutation, EGFR amplification and PTEN loss. Seven (70%) malignant phyllodes tumors harbored TERT aberrations (six promoter mutations, one amplification). For comparison, TERT promoter mutations were identified by Sanger sequencing in 33% borderline (n=12) and no (0%, n=8) benign phyllodes tumors (P=0.391 and P=0.013 vs malignant tumors, respectively). Genetic features specific to liposarcoma, including CDK4/MDM2 amplification, were not identified. Copy number analysis revealed intratumoral heterogeneity and evidence for divergent tumor evolution in malignant phyllodes tumors with and without heterologous differentiation. Tumors with liposarcomatous differentiation revealed more chromosomal aberrations in non-heterologous components compared with liposarcomatous components. EGFR amplification was heterogeneous and present only in the non-heterologous component of one tumor with liposarcomatous differentiation. The results identify novel pathways involved in the pathogenesis of malignant phyllodes tumors, which significantly increase our understanding of tumor biology and have potential clinical impact.

RNA-sequencing is most commonly used to measure gene expression, but it is possible to extract genotypic information from RNA-sequencing data, too. Point mutations and translocations can be detected when they occur in expressed genes, however, there are few software solutions to infer copy number information from RNA-sequencing data. This is because a gene’s expression is dictated by a number of variables, including, but not limited to, copy number variation. Here, we report new functionalities within the software package CNVkit that enable copy number inference from RNA-sequencing data. First, CNVkit removes technical variation in gene expression associated with GC-content and transcript length. Next, CNVkit assigns a weight, dictated by several variables, to each transcript with the net effect of preferentially inferring copy number from highly and stably expressed genes. We benchmarked our approach on 105 melanomas from The Cancer Genome Atlas project and observed a high degree of concordance (R = 0.739) between our estimates and those from array comparative genomic hybridization (aCGH) on the same samples. After initial configuration, the software requires few inputs, is able to process a batch of up to 100 samples in less than ten minutes, and can be used in conjunction with pre-existing features of CNVkit, including visualization tools. Overall, we present a rapid, user-friendly software solution to infer copy number information from gene expression data.

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.

Desmoplastic melanoma is an uncommon variant of melanoma with sarcomatous histology, distinct clinical behavior and unknown pathogenesis. We performed low-coverage genome and high-coverage exome sequencing of 20 desmoplastic melanomas, followed by targeted sequencing of 293 genes in a validation cohort of 42 cases. A high mutation burden (median of 62 mutations/Mb) ranked desmoplastic melanoma among the most highly mutated cancers. Mutation patterns strongly implicate ultraviolet radiation as the dominant mutagen, indicating a superficially located cell of origin. Newly identified alterations included recurrent promoter mutations of NFKBIE, encoding NF-κB inhibitor ɛ (IκBɛ), in 14.5% of samples. Common oncogenic mutations in melanomas, in particular in BRAF (encoding p.Val600Glu) and NRAS (encoding p.Gln61Lys or p.Gln61Arg), were absent. Instead, other genetic alterations known to activate the MAPK and PI3K signaling cascades were identified in 73% of samples, affecting NF1, CBL, ERBB2, MAP2K1, MAP3K1, BRAF, EGFR, PTPN11, MET, RAC1, SOS2, NRAS and PIK3CA, some of which are candidates for targeted therapies.

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 which 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 2 cases demonstrating monoallelic loss of BAP1. All 11 of these cases demonstrated loss of BAP1 nuclear staining by immunohistochemistry, whereas the 2 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.

Secretory carcinomas of the breast are rare tumors with distinct histologic features, recurrent t(12;15)(p13;q25) translocation resulting in ETV6-NTRK3 gene fusion and indolent clinical behavior. Mammary analog secretory carcinomas arising in other sites are histopathologically similar to the breast tumors and also harbor ETV6-NTRK3 fusions. Breast secretory carcinomas are often triple (estrogen and progesterone receptor, HER2) negative with a basal-like immunophenotype. However, genomic studies are lacking, and whether these tumors share genetic features with other basal and/or triple negative breast cancers is unknown. Aside from shared ETV6-NTRK3 fusions, the genetic relatedness of secretory carcinomas arising in different sites is also uncertain. We immunoprofiled and sequenced 510 cancer-related genes in nine breast secretory carcinomas and six salivary gland mammary analog secretory carcinomas. Immunoprofiles of breast and salivary gland secretory carcinomas were similar. All the tumors showed strong diffuse MUC4 expression (n=15), and SOX10 was positive in all nine breast and in five out of six salivary gland tumors. All breast secretory carcinomas were triple negative or weakly ER-positive, and all tumors at both the sites expressed CK5/6 and/or EGFR, consistent with a basal-like phenotype. Sequencing revealed classic ETV6-NTRK3 fusion genes in all cases, including in carcinoma in situ of one breast tumor. Translocations were reciprocal and balanced in six out of nine breast and three out of six salivary gland tumors and were complex in three others. In contrast to most breast basal carcinomas, the mutational burden of secretory carcinomas was very low, and no additional pathogenic aberrations were identified in genes typically mutated in breast cancer. Five (56%) breast and two (33%) salivary gland tumors had simple genomes without copy number changes; the remainder had very few changes, averaging 1.3 per tumor. The ETV6-NTRK3 derivative chromosome was duplicated in one breast and one salivary gland tumor, and was the only copy number change in the latter. The findings highlight breast secretory carcinoma as a subtype more closely related to mammary analog secretory carcinoma than to basal/triple negative breast cancers of no special type. Lack of pathogenic mutations in common cancer-related genes suggests that ETV6-NTRK3 alone may suffice to drive these tumors and likely helps explain their indolent behavior.

Protein kinases represent a large and diverse family of evolutionarily related proteins that are abnormally regulated in human cancers. Although genome sequencing studies have revealed thousands of variants in protein kinases, translating “big” genomic data into biological knowledge remains a challenge. Here, we describe an ontological framework for integrating and conceptualizing diverse forms of information related to kinase activation and regulatory mechanisms in a machine readable, human understandable form. We demonstrate the utility of this framework in analyzing the cancer kinome, and in generating testable hypotheses for experimental studies. Through the iterative process of aggregate ontology querying, hypothesis generation and experimental validation, we identify a novel mutational hotspot in the αC-β4 loop of the kinase domain and demonstrate the functional impact of the identified variants in epidermal growth factor receptor (EGFR) constitutive activity and inhibitor sensitivity. We provide a unified resource for the kinase and cancer community, ProKinO, housed at http://vulcan.cs.uga.edu/prokino.

© 2017 The Author. Background. Molecular profling is revolutionizing cancer diagnostics and leading to personalized therapeutic approaches. Herein we describe our clinical experience performing targeted sequencing for 31 pediatric neurooncology patients. Methods. We sequenced 510 cancer-associated genes from tumor and peripheral blood to identify germline and somatic mutations, structural variants, and copy number changes. Results. Genomic profling was performed on 31 patients with tumors including 11 high-grade gliomas, 8 medulloblastomas, 6 low-grade gliomas, 1 embryonal tumor with multilayered rosettes, 1 pineoblastoma, 1 uveal ganglioneuroma, 1 choroid plexus carcinoma, 1 chordoma, and 1 high-grade neuroepithelial tumor. In 25 cases (81%), results impacted patient management by: (i) clarifying diagnosis, (ii) identifying pathogenic germline mutations, or (iii) detecting potentially targetable alterations. The pathologic diagnosis was amended after genomic profling for 6 patients (19%), including a high-grade glioma to pilocytic astrocytoma, medulloblastoma to pineoblastoma, ependymoma to high-grade glioma, and medulloblastoma to CNS high-grade neuroepithelial tumor with BCOR alteration. Multiple patients had pathogenic germline mutations, many of which were previously unsuspected. Potentially targetable alterations were identifed in 19 patients (61%). Additionally, novel likely pathogenic alterations were identifed in 3 cases: an in-frame RAF1 fusion in a BRAF wild-type pleomorphic xanthoastrocytoma, an inactivating ASXL1 mutation in a histone H3 wild-type diffuse pontine glioma, and an in-frame deletion within exon 2 of MAP2K1 in a low-grade astrocytic neoplasm. Conclusions. Our experience demonstrates the signifcant impact of molecular profling on diagnosis and treatment of pediatric brain tumors and confrms its feasibility for use at the time of diagnosis or recurrence.

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.

Background Many prokaryotic kinases that phosphorylate small molecule substrates, such as antibiotics, lipids and sugars, are evolutionarily related to Eukaryotic Protein Kinases (EPKs). These Eukaryotic-Like Kinases (ELKs) share the same overall structural fold as EPKs, but differ in their modes of regulation, substrate recognition and specificity—the sequence and structural determinants of which are poorly understood. Results To better understand the basis for ELK specificity, we applied a Bayesian classification procedure designed to identify sequence determinants responsible for functional divergence. This reveals that a large and diverse family of aminoglycoside kinases, characterized members of which are involved in antibiotic resistance, fall into major sub-groups based on differences in putative substrate recognition motifs. Aminoglycoside kinase substrate specificity follows simple rules of alternating hydroxyl and amino groups that is strongly correlated with variations at the DFG + 1 position. Conclusions Substrate specificity determining features in small molecule kinases are mostly confined to the catalytic core and can be identified based on quantitative sequence and crystal structure comparisons. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0576-x) contains supplementary material, which is available to authorized users.

The pathogenic mutations in melanoma have been largely catalogued; however, the order of their occurrence is not known.We sequenced 293 cancer-relevant genes in 150 areas of 37 primary melanomas and their adjacent precursor lesions. The histopathological spectrum of these areas included unequivocally benign lesions, intermediate lesions, and intraepidermal or invasive melanomas.Precursor lesions were initiated by mutations of genes that are known to activate the mitogen-activated protein kinase pathway. Unequivocally benign lesions harbored BRAF V600E mutations exclusively, whereas those categorized as intermediate were enriched for NRAS mutations and additional driver mutations. A total of 77% of areas of intermediate lesions and melanomas in situ harbored TERT promoter mutations, a finding that indicates that these mutations are selected at an unexpectedly early stage of the neoplastic progression. Biallelic inactivation of CDKN2A emerged exclusively in invasive melanomas. PTEN and TP53 mutations were found only in advanced primary melanomas. The point-mutation burden increased from benign through intermediate lesions to melanoma, with a strong signature of the effects of ultraviolet radiation detectable at all evolutionary stages. Copy-number alterations became prevalent only in invasive melanomas. Tumor heterogeneity became apparent in the form of genetically distinct subpopulations as melanomas progressed.Our study defined the succession of genetic alterations during melanoma progression, showing distinct evolutionary trajectories for different melanoma subtypes. It identified an intermediate category of melanocytic neoplasia, characterized by the presence of more than one pathogenic genetic alteration and distinctive histopathological features. Finally, our study implicated ultraviolet radiation as a major factor in both the initiation and progression of melanoma. (Funded by the National Institutes of Health and others.).

Prostate cancer is a highly heterogeneous disease with distinct genomic underpinnings and generally presents a modest number of genome wide somatic single point mutations and indels at approximately one or two per megabase, but frequently has large-scale somatic copy number and structural alterations. We recently showed that 25-30% of prostate tumors bear defects in double-strand DNA repair genes that render them sensitive to treatment with PARP inhibition, and that these are commonly mediated through copy number events. Accurately estimate gene copy number variation (CNV) is crucial to understanding the genomic background of advanced prostate tumours. However, it has been challenging in targeted amplicon-based next generation sequencing (NGS). Here, we describe the analytical validation and clinical qualification of a CNV assessment, based on targeted amplicon sequencing of a focused biomarker gene panel dedicated to identifying DNA repair defects, with an aim to identify patients that may benefit from PARP inhibitor treatment. Our targeted NGS protocol is used to profile a cohort of 110 castration-resistant prostate cancer patients. The key findings include AR amplification (Log2 > 2) in 45.5% (50/110) of samples and commonly shared loss of BRCA2 and RB1 (88% of samples reporting any BRCA2 loss also had some loss of RB1, while 33.8% of samples with RB1 loss had loss of BRCA2). Moreover, we validated our results using other CNV estimation platforms to establish that: 1) the protocol produces highly reproducible results (r = 0.98); 2) in the covered regions, CNV estimated from targeted amplicon NGS is comparable with capture based whole exome sequencing (r = 0.86) as well as array CGH (r = 0.81); and 3) CNV estimation can be validated by orthogonal methods such as ddPCR and FISH. The approach has the benefit of estimating CNV from tumour-only sample by using pooled (n=34) germline sequencing data as a reference. In this study, we have developed an accurate, affordable and rapid approach to estimate copy number aberrations using targeted gene panel sequencing. We demonstrated high reproducibility, validating our results with different orthogonal copy number estimation platforms. And this protocol has potential to be integrated with other protocols for patient assessment and precision medicine. Citation Format: George Seed, Wei Yuan, Joaquin Mateo, Suzanne Carreira, Maryou Lambros, Gunther Boysen, Roberta Ferraldeschi, Susana Miranda, Ines Figueiredo, Ruth Riisnaes, Mateus Crespo, Daniel Nava Rodrigues, Eric Talevich, Dan Robinson, Priya Kunju, Yi-mi Wu, Robert Lonigro, Shahneen Sandhu, Arul Chinnayan, Johann De Bono. Copy number estimation from targeted amplicon-based next-generation sequencing of castration-resistant prostate cancer biopsies: analytic validation and clinical qualification for a iPARP clinical trial [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-044. doi:10.1158/1538-7445.AM2017-LB-044

11575 Background: Targeted next-generation DNA sequencing of paired tumor and normal DNA samples allows for detection of biologically relevant variants in the tumor with significantly greater accuracy than tumor-only sequencing. In addition, this approach presents with the opportunity of unveiling previously unknown cancer predisposition traits in a patient’s germline DNA. Methods: We sought to determine the rate of pathogenic germline alterations in 546 consecutive pediatric and adult patients who underwent molecular profiling using the UCSF 500 assay, a hybrid capture-based DNA sequencing assay targeting the coding regions of ~500 cancer-related genes, TERT promoter, select introns from 40 genes (for detection of gene fusions and other structural variants), and intergenic regions at regular intervals along each chromosome (for chromosomal copy number and LOH assessment). Results: Pathogenic germline alterations were found in 89/546 patients (16.3%), including 25 pediatric and 64 adult cases. Germline variants were identified in 37 genes with MUTYH (n = 15, 17%), CHEK2 (n = 10, 11%), BRCA2 (n = 9, 10%), BRCA1 (n = 5, 6%), TP53 (n = 5, 5%), and APC (n = 4, 5%) being altered most frequently. Loss of heterozygosity of genes affected in the germline was seen in tumor DNA 37 (42%) cases, highlighting their likely role as drivers of tumorigenesis. Clinically relevant germline findings not associated with increased cancer risk were identified in 6 (7%) of the cases, for example a COL1A1 mutation associated with Ehlers-Danlos Syndrome. In In 73 (82%) of the cases, pathogenic germline alterations were new findings and genetic counseling was recommended. A possible, previously unknown, role of germline mutations was found in some instances, including a TSC2 germline mutation in a patient with hybrid oncocytoma/chromophobe tumor (HOCT) with loss of the normal TSC2 allele in the tumor. Conclusions: Our data suggest that paired tumor/normal DNA analysis uncovers actionable heritable traits in a substantial fraction of patients and represents the preferred approach to analyzing malignancies in children and adults.

Germline copy number variants (CNVs) and somatic copy number alterations (SCNAs) are of significant importance in syndromic conditions and cancer. Massive 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 compatibility with other software. Availability: http://github.com/etal/cnvkit

Cancer is a genetic disease that develops through a series of somatic mutations, a subset of which drive cancer progression. Although cancer genome sequencing studies are beginning to reveal the mutational patterns of genes in various cancers, identifying the small subset of “causative” mutations from the large subset of “non-causative” mutations, which accumulate as a consequence of the disease, is a challenge. In this article, we present an effective machine learning approach for identifying cancer-associated mutations in human protein kinases, a class of signaling proteins known to be frequently mutated in human cancers. We evaluate the performance of 11 well known supervised learners and show that a multiple-classifier approach, which combines the performances of individual learners, significantly improves the classification of known cancer-associated mutations. We introduce several novel features related specifically to structural and functional characteristics of protein kinases and find that the level of conservation of the mutated residue at specific evolutionary depths is an important predictor of oncogenic effect. We consolidate the novel features and the multiple-classifier approach to prioritize and experimentally test a set of rare unconfirmed mutations in the epidermal growth factor receptor tyrosine kinase (EGFR). Our studies identify T725M and L861R as rare cancer-associated mutations inasmuch as these mutations increase EGFR activity in the absence of the activating EGF ligand in cell-based assays., Author Summary Cancer progresses by accumulation of mutations in a subset of genes that confer growth advantage. The 518 protein kinase genes encoded in the human genome, collectively called the kinome, represent one of the largest families of oncogenes. Targeted sequencing studies of many different cancers have shown that the mutational landscape comprises both cancer-causing “driver” mutations and harmless “passenger” mutations. While the frequent recurrence of some driver mutations in human cancers helps distinguish them from the large number of passenger mutations, a significant challenge is to identify the rare “driver” mutations that are less frequently observed in patient samples and yet are causative. Here we combine computational and experimental approaches to identify rare cancer-associated mutations in Epidermal Growth Factor receptor kinase (EGFR), a signaling protein frequently mutated in cancers. Specifically, we evaluate a novel multiple-classifier approach and features specific to the protein kinase super-family in distinguishing known cancer-associated mutations from benign mutations. We then apply the multiple classifier to identify and test the functional impact of rare cancer-associated mutations in EGFR. We report, for the first time, that the EGFR mutations T725M and L861R, which are infrequently observed in cancers, constitutively activate EGFR in a manner analogous to the frequently observed driver mutations.

During asexual intraerythrocytic development, Plasmodium falciparum diverges from the paradigm of the eukaryotic cell cycles by undergoing multiple rounds of DNA replication and nuclear division without cytokinesis. A better understanding of the molecular switches that coordinate a myriad of events for the progression of the parasite through the intraerythrocytic developmental stages will be of fundamental importance for rational design of intervention strategies. To achieve this goal, we performed isobaric tag-based quantitative proteomics and phosphoproteomics analyses of three developmental stages in the Plasmodium asexual cycle and identified 2767 proteins, 1337 phosphoproteins, and 6293 phosphorylation sites. Approximately 34% of identified proteins and 75% of phosphorylation sites exhibit changes in abundance as the intraerythrocytic cycle progresses. Our study identified 43 distinct phosphorylation motifs and a range of potential MAPK/CDK substrates. Further analysis of phosphorylated kinases identified 30 protein kinases with 126 phosphorylation sites within the kinase domain or in N- or C-terminal tails. Many of these phosphorylations are likely CK2-mediated. We define the constitutive and regulated expression of the Plasmodium proteome during the intraerythrocytic developmental cycle, offering an insight into the dynamics of phosphorylation during asexual cycle progression. Our system-wide comprehensive analysis is a major step toward defining kinase-substrate pairs operative in various signaling networks in the parasite.

The pathogenic mutations in melanoma have largely been catalogued, but the order of their occurrence is not known. We identified 82 cases of melanocytic neoplasia with two or more histopathologically distinct precursor and descendent portions. These portions spanned several histopathologic stages including benign nevi, intermediate neoplasms such as dysplastic nevi, melanoma in situ, invasive melanoma, and metastatic melanoma. In total, 277 histopathogically distinct areas were microdissected and sequenced using a panel of several hundred cancer-genes at an average of 275-fold coverage. In all cases precursor lesions harbored mutations known to activate the MAPK pathway, most commonly affecting BRAF or NRAS, and thus these were considered initiating mutations. TERT promoter mutations were the earliest secondary mutations to occur, emerging in intermediate neoplasms and in situ melanomas. Bi-allelic CDKN2A mutations and/or deletions arose predominantly at the transition to invasive melanoma. TP53 and PTEN mutations occurred comparatively later in thicker primary melanomas. In the 12 cases in which matched primaries and metastases were sequenced, no pathogenic mutations were specifically associated with the transition to metastatic melanoma, suggesting that metastatic capability was already present in the primary tumors. Unequivocally benign neoplasms, which mostly consisted of conventional nevi, disproportionately harbored BRAFV600E mutations as the only apparent pathogenic alteration. Lesions that were histopathologically intermediate, such as dysplastic nevi, had two or more pathogenic mutations and were enriched for NRAS mutations. These findings provide genetic support for the existence of an intermediate stage of melanocytic neoplasia, resolving a long-standing controversy. We also found that melanomas with different initiating mutations evolve through distinct evolutionary trajectories linked to specific histopathologic precursors. A mutational signature of UV-radiation-induced DNA damage predominated in lesions from all histopathologic stages, implicating UV-radiation in both the initiation and progression of melanoma. The point mutation burden increased steadily from one progression stage to the next, stabilizing once melanomas became invasive. By contrast, copy number alterations only became prevalent at the transition to invasive melanoma. Melanomas evolve from benign lesions linearly, whereas, they show patterns of branched evolution when they evolve from intermediate lesions; this indicates that evolution may accelerate once a melanocytic neoplasm reaches an intermediate stage. In summary, our study outlines a framework of the genetic evolution of melanocytic neoplasms from precursor lesions and unveils the rate-limiting homeostatic factors that become disrupted as neoplasms progress from one stage to the next, as well as the pathogenic factors driving this evolution. Citation Format: Alan H. Shain, Richard Yu, Iwei Yeh, Jamal Benhamida, Ivanka Kovalyshyn, Aravindhan Sriharan, Eric Talevich, Reinhard Dummer, Jeffrey North, Laura Pincus, Beth Ruben, William Rickaby, Corrado D’Arrigo, Alistair Robson, Robert Judson, Nancy Joseph, Boris Bastian. The genetic evolution of melanoma. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2372.

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.

[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., We acknowledge the role of Google Summer of Code 2009 and the National Evolutionary Synthesis Center (NESCent) Phyloinformatics program in sponsoring and initiating this project. BI is supported by a FI-DGR from AGAUR (Catalonia, Spain). Harvard Bioinformatics Core provided for publication costs and moral support.

Desmoplastic melanomas (DMs) comprise 4% of the overall melanoma burden and have a 5-year survival rate of 85%. DMs are dermal tumors characterized by spindled melanocytes situated within abundant desomplastic stroma. These unusual histological features commonly lead to misdiagnosis. Currently, there are no known genetic drivers. A better understanding of the underlying biology of desmoplastic melanoma would provide biomarkers and therapeutic opportunities. Towards this goal, we performed low-coverage genome and high-coverage exome sequencing of 20 DMs in a discovery cohort, followed by targeted sequencing of 293 candidate genes on a validation cohort of 42 cases. Additionally, high-resolution aCGH was performed on samples from both cohorts. A high mutation burden (median 62 mutations/Mb) ranked desmoplastic melanoma among the most highly mutated cancers sequenced to date. Mutation patterns strongly indicate that UV-radiation is the dominant mutagen and implicate a superficially located cell of origin despite their predominantly intradermal presentation. Novel alterations included recurrent promoter mutations and amplification of NF-kappa B inhibitor epsilon, NFKBIE (IkBϵ) in 14.5% of samples. The promoter mutations typically affect both alleles and occur over a highly conserved DNA region. The mutations are predicted to disrupt a canonical Ets Like Factor 1 (ELF1) binding site. In total, these data imply aberrant NF-kappa B signaling as a pathogenic feature of desmoplastic melanoma. Commonly mutated oncogenes in melanomas, in particular BRAF V600E and NRAS Q61K/R, were absent. Instead, other genetic alterations known to activate the MAPK and PI3K signaling cascades were identified in 73% of samples, affecting NF1, CBL, ERBB2, MAP2K1, MAP3K1, BRAF, EGFR, PTPN11, MET, RAC1, SOS2, NRAS, and PIK3CA. Rb and p53 pathway alterations occurred respectively in 71% and 66% of tumors, affecting RB1, FBXW7, CDK4, PPP6C, CCND1, CDKN2A, TP53, and MDM2. Finally, TERT promoter mutations or amplifications occurred in 90% of tumors. The consequences of the mutations on protein expression levels was confirmed by immunostaining for NF1, EGFR, Rb, CDK4, CCND1, p16, p53, and Mdm2. Collectively, many of these oncogenic mutations are potentially druggable. In conclusion, desmoplastic melanomas harbor distinct genetic alterations that explain their unique biology, and this study illuminates genetic biomarkers and nominates targets for therapeutic intervention. Citation Format: Alan H. Shain, Maria Garrido, Thomas Botton, Eric Talevich, Iweh Yeh, Zack Sanborn, Jongsuk Chung, Nicholas Wang, Hojabr Kakavand, Graham Mann, John Thompson, Thomas Wiesner, Ritu Roy, Adam Olshen, Alexander Gagnon, Joe Gray, Nam Huh, Joe Hur, Klaus Busam, Richard Scolyer, Raymond Cho, Rajmohan Murali, Boris Bastian. Exome sequencing of desmoplastic melanoma reveals recurrent NFKBIE promoter mutations and diverse MAPK/PI3K pathway activating mutations. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2968. doi:10.1158/1538-7445.AM2015-2968

Oncogenic gene fusions have been identified in many cancers and many serve as biomarkers or targets for therapy. Here we identify six different melanocytic tumours with genomic rearrangements of MET fusing the kinase domain of MET in-frame to six different N-terminal partners. These tumours lack activating mutations in other established melanoma oncogenes. We functionally characterize two of the identified fusion proteins (TRIM4-MET and ZKSCAN1-MET) and find that they constitutively activate the mitogen-activated protein kinase (MAPK), phosphoinositol-3 kinase (PI3K) and phospholipase C gamma 1 (PLCγ1) pathways. The MET inhibitors cabozantinib (FDA-approved for progressive medullary thyroid cancer) and PF-04217903 block their activity at nanomolar concentrations. MET fusion kinases thus provide a potential therapeutic target for a rare subset of melanoma for which currently no targeted therapeutic options currently exist.

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.

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.