Faculty Bibliography

The goal of this talk is to introduce an integrated quantitative proteogenomic approach to comprehensively map proteomic information back to their encoding genes. We seek evidence from mass spectrometry-based large-scale proteomic data of patient populations in conjunction with patient-centric next-generation sequencing data and unbiased sequencing strategies to study breast cancer (BC) subtypes from a genomic context. We have obtained global and phosphoproteomic data with matching next generation sequencing data for 18 patient-derived xenografts (PDXs) representing the major clinical subtypes of BC. Our workflow starts with the creation of several protein sequence databases that serve as a template for mass spectrometry database identifications. These databases include 1) completely annotated reference protein sequences, 2) patient-specific databases that were created using next generation sequencing data, 3) isoform databases that contain all possible splicing combinations, and 4) amino acid sequence database resulting from a six-frame translation of the entire human reference and customized genomes. All mass spectrometry raw data are searched against the databases for obtaining identifications at the peptide level, and assembly of peptides for quantification using taxonomy-based label-free quantitation (LFQ) that can specifically quantify unique human peptide sequences found in PDXs. The peptides are then mapped to the human genome and visualized using a genome browser. Quantitative changes across PDXs are presented at the protein level or at the isoform level via peptide role-up to specific exons and visualized as a quantitative data track. By combining search results from these databases we obtain a comprehensive view of our PDXs. The complementary nature of the databases enable greater proteomic depth, i.e. databases with complete splicing combinations capture proteomic evidence when patient-specific databases fail due to possible erroneous RNA-seq reads. Similarly, 6-frame translated amino acid databases can capture potentially novel coding regions but are unable to detect splicing. Peptide maps are obtained for individual genes or specific protein isoforms covering both knowledge-driven and novel genomic annotation types. We compile peptides carrying variants, splice junctions, fusions, and new coding regions specific to each PDX or in common with a specific BC subtype. Majority of data is mapped back to the genome loci using unmodified peptides via global proteomics while phosphopeptides that contain variants and splicing are also mapped in a similar manner using phosphoproteomic data. We have currently annotated 455 novel proteogenomic hits covering many examples outlined above for genes related to breast cancer and show how these can be specifically identified and in some instances differentially quantified in the PDX models

A significant fraction of cancers have a heritable component, and require an interplay between somatic and germline variants. Common and rare germline variants have been investigated by previous GWAS and family based studies. However, a comprehensive analysis of both somatic and germline variants in cancer using high throughput sequencing data to discover genetic variants of functional relevance is lacking. Herein, we investigated the potential role of somatic and germline variants from over 20 major cancer types from large-scale studies such as TCGA and ICGC and discovered thousands of somatic and germline variants in cancer genomes. To link these genetic variants to cancer phenotypes, we analyzed the proteomics data in breast, ovarian and colorectal cancers generated by the Clinical Proteomic Tumour Analysis Consortium (CPTAC) using a PepScan pipeline that can detect whether a genomic variant is observed at the peptide level. The pipeline incorporated QUILTS to construct patient-specific protein database and MS-GF+ to identify peptide sequences in the database from MS spectra. Our analysis validated roughly 2% of non-synonymous genetic variants in peptides with matched spectra. We correlated the effect of genetic mutations on proteomic subtypes based on global protein expression levels. Additionally, we assessed the role of genetic variations in kinase genes using phosphoproteome profiles, and identified downstream markers that may be candidate targets for diagnosis or treatment. This analysis also helped us prioritize kinase variants that are likely functional candidates for experimental validation. In conclusion, the comprehensive study of genetic variants utilizing an integrated proteogenomic approach revealed genetic variants with potential functional impacts in cancer

Nanobodies are single-domain antibodies derived from the variable regions of Camelidae atypical immunoglobulins. They show promise as high-affinity reagents for research, diagnostics and therapeutics owing to their high specificity, small size ( approximately 15 kDa) and straightforward bacterial expression. However, identification of repertoires with sufficiently high affinity has proven time consuming and difficult, hampering nanobody implementation. Our approach generates large repertoires of readily expressible recombinant nanobodies with high affinities and specificities against a given antigen. We demonstrate the efficacy of this approach through the production of large repertoires of nanobodies against two antigens, GFP and mCherry, with Kd values into the subnanomolar range. After mapping diverse epitopes on GFP, we were also able to design ultrahigh-affinity dimeric nanobodies with Kd values as low as approximately 30 pM. The approach presented here is well suited for the routine production of high-affinity capture reagents for various biomedical applications.

AIMS: It is well-known that connexin43 (Cx43) forms gap junctions. We recently showed that Cx43 is also part of a protein interacting network that regulates excitability. Cardiac-specific truncation of Cx43 C-terminus (mutant "Cx43D378stop") led to lethal arrhythmias. Cx43D378stop localized to the intercalated disc (ID); cell-cell coupling was normal, but there was significant sodium current (INa) loss. We proposed that the microtubule plus-end is at the crux of the Cx43-INa relation. Yet, specific localization of relevant molecular players was prevented due to the resolution limit of fluorescence microscopy. Here, we use nanoscale imaging to establish: a) the morphology of clusters formed by the microtubule plus-end tracking protein "end binding 1" (EB1), b) their position, and that of sodium channel alpha-subunit NaV1.5, relative to N-cadherin rich sites, c) the role of Cx43 C-terminus on the above-mentioned parameters and on the location-specific function of INa. METHODS AND RESULTS: Super-resolution fluorescence localization microscopy in murine adult cardiomyocytes revealed EB1 and NaV1.5 as distinct clusters preferentially localized to N-cadherin-rich sites. Extent of co-localization decreased in Cx43D378stop cells. Macropatch and scanning patch clamp showed reduced INa exclusively at cell end, without changes in unitary conductance. Experiments in Cx43-modified HL1 cells confirmed the relation between Cx43, INa and microtubules. CONCLUSIONS: NaV1.5 and EB1 localization at cell end is Cx43-dependent. Cx43 is part of a molecular complex that determines capture of the microtubule plus-end at the ID, facilitating cargo delivery. These observations link excitability and electrical coupling through a common molecular mechanism.

Extensive genomic characterization of human cancers presents the problem of inference from genomic abnormalities to cancer phenotypes. To address this problem, we analysed proteomes of colon and rectal tumours characterized previously by The Cancer Genome Atlas (TCGA) and perform integrated proteogenomic analyses. Somatic variants displayed reduced protein abundance compared to germline variants. Messenger RNA transcript abundance did not reliably predict protein abundance differences between tumours. Proteomics identified five proteomic subtypes in the TCGA cohort, two of which overlapped with the TCGA 'microsatellite instability/CpG island methylation phenotype' transcriptomic subtype, but had distinct mutation, methylation and protein expression patterns associated with different clinical outcomes. Although copy number alterations showed strong cis- and trans-effects on mRNA abundance, relatively few of these extend to the protein level. Thus, proteomics data enabled prioritization of candidate driver genes. The chromosome 20q amplicon was associated with the largest global changes at both mRNA and protein levels; proteomics data highlighted potential 20q candidates, including HNF4A (hepatocyte nuclear factor 4, alpha), TOMM34 (translocase of outer mitochondrial membrane 34) and SRC (SRC proto-oncogene, non-receptor tyrosine kinase). Integrated proteogenomic analysis provides functional context to interpret genomic abnormalities and affords a new paradigm for understanding cancer biology.

Protein abundance and phosphorylation convey important information about pathway activity and molecular pathophysiology in diseases including cancer, providing biological insight, informing drug and diagnostic development, and guiding therapeutic intervention. Analyzed tissues are usually collected without tight regulation or documentation of ischemic time. To evaluate the impact of ischemia, we collected human ovarian tumor and breast cancer xenograft tissue without vascular interruption and performed quantitative proteomics and phosphoproteomics after defined ischemic intervals. While the global expressed proteome and most of the >25,000 quantified phosphosites were unchanged after 60 minutes, rapid phosphorylation changes were observed in up to 24% of the phosphoproteome, representing activation of critical cancer pathways related to stress response, transcriptional regulation and cell death. Both pan-tumor and tissue-specific changes were observed. The demonstrated impact of pre-analytical tissue ischemia on tumor biology mandates caution in interpreting stress-pathway activation in such samples, and motivates reexamination of collection protocols for phosphoprotein analysis.

One subproject within the global Chromosome 19 Consortium is to define chromosome 19 gene and protein expression in glioma-derived cancer stem cells (GSCs). Chromosome 19 is notoriously linked to glioma by 1p/19q codeletions, and clinical tests are established to detect that specific aberration. GSCs are tumor-initiating cells and are hypothesized to provide a repository of cells in tumors that can self-replicate and be refractory to radiation and chemotherapeutic agents developed for the treatment of tumors. In this pilot study, we performed RNA-Seq, label-free quantitative protein measurements in six GSC lines, and targeted transcriptomic analysis using a chromosome 19-specific microarray in an additional six GSC lines. The data have been deposited to the ProteomeXchange with identifier PXD000563. Here we present insights into differences in GSC gene and protein expression, including the identification of proteins listed as having no or low evidence at the protein level in the Human Protein Atlas, as correlated to chromosome 19 and GSC subtype. Furthermore, the upregulation of proteins downstream of adenovirus-associated viral integration site 1 (AAVS1) in GSC11 in response to oncolytic adenovirus treatment was demonstrated. Taken together, our results may indicate new roles for chromosome 19, beyond the 1p/19q codeletion, in the future of personalized medicine for glioma patients.