Eric Talevich

The Genetic Evolution of Melanoma from Precursor Lesions

BACKGROUND The pathogenic mutations in melanoma have been largely catalogued; however, the order of their occurrence is not known. METHODS 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. RESULTS 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. CONCLUSIONS 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.).

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.

Exome sequencing of desmoplastic melanoma identifies recurrent NFKBIE promoter mutations and diverse activating mutations in the MAPK pathway

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.

Global Analysis of Protein Expression and Phosphorylation of Three Stages of Plasmodium falciparum Intraerythrocytic Development

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.

Bio.Phylo: A unified toolkit for processing, analyzing and visualizing phylogenetic trees in Biopython

BackgroundOngoing 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.ResultsWe 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.ConclusionsBio.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.

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

BackgroundThe 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.ResultsIn 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.ConclusionsOur 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.

Activating MET kinase rearrangements in melanoma and Spitz tumours

Oncogenic gene fusions have been identified in many cancers and many serve as biomarkers or targets for therapy. Here we identify six different melanocytic tumors with genomic rearrangements of MET fusing the kinase domain of MET in-frame to six different N-terminal partners. These tumors 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.

Targeted next-generation sequencing of pediatric neuro-oncology patients improves diagnosis, identifies pathogenic germline mutations, and directs targeted therapy

Background Molecular profiling is revolutionizing cancer diagnostics and leading to personalized therapeutic approaches. Herein we describe our clinical experience performing targeted sequencing for 31 pediatric neuro-oncology 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 profiling 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 profiling 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 identified in 19 patients (61%). Additionally, novel likely pathogenic alterations were identified 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 significant impact of molecular profiling on diagnosis and treatment of pediatric brain tumors and confirms its feasibility for use at the time of diagnosis or recurrence.

Structural and evolutionary divergence of eukaryotic protein kinases in Apicomplexa

BackgroundThe 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.ResultsWe 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.ConclusionsWe 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.

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 parasites kinome and phosphoproteome.