Nathan Edwards

Proteogenomic characterization of human colon and rectal cancer

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.

Integrated proteogenomic characterization of human high-grade serous ovarian cancer

To provide a detailed analysis of the molecular components and underlying mechanisms associated with ovarian cancer, we performed a comprehensive mass-spectrometry-based proteomic characterization of 174 ovarian tumors previously analyzed by The Cancer Genome Atlas (TCGA), of which 169 were high-grade serous carcinomas (HGSCs). Integrating our proteomic measurements with the genomic data yielded a number of insights into disease, such as how different copy-number alternations influence the proteome, the proteins associated with chromosomal instability, the sets of signaling pathways that diverse genome rearrangements converge on, and the ones most associated with short overall survival. Specific protein acetylations associated with homologous recombination deficiency suggest a potential means for stratifying patients for therapy. In addition to providing a valuable resource, these findings provide a view of how the somatic genome drives the cancer proteome and associations between protein and post-translational modification levels and clinical outcomes in HGSC. VIDEO ABSTRACT.

Whole-genome shotgun assembly and comparison of human genome assemblies

We report a whole-genome shotgun assembly (called WGSA) of the human genome generated at Celera in 2001. The Celera-generated shotgun data set consisted of 27 million sequencing reads organized in pairs by virtue of end-sequencing 2-kbp, 10-kbp, and 50-kbp inserts from shotgun clone libraries. The quality-trimmed reads covered the genome 5.3 times, and the inserts from which pairs of reads were obtained covered the genome 39 times. With the nearly complete human DNA sequence [National Center for Biotechnology Information (NCBI) Build 34] now available, it is possible to directly assess the quality, accuracy, and completeness of WGSA and of the first reconstructions of the human genome reported in two landmark papers in February 2001 [Venter, J. C., Adams, M. D., Myers, E. W., Li, P. W., Mural, R. J., Sutton, G. G., Smith, H. O., Yandell, M., Evans, C. A., Holt, R. A., et al. (2001) Science 291, 1304–1351; International Human Genome Sequencing Consortium (2001) Nature 409, 860–921]. The analysis of WGSA shows 97% order and orientation agreement with NCBI Build 34, where most of the 3% of sequence out of order is due to scaffold placement problems as opposed to assembly errors within the scaffolds themselves. In addition, WGSA fills some of the remaining gaps in NCBI Build 34. The early genome sequences all covered about the same amount of the genome, but they did so in different ways. The Celera results provide more order and orientation, and the consortium sequence provides better coverage of exact and nearly exact repeats.

Deconvolution and Database Search of Complex Tandem Mass Spectra of Intact Proteins A COMBINATORIAL APPROACH

Top-down proteomics studies intact proteins, enabling new opportunities for analyzing post-translational modifications. Because tandem mass spectra of intact proteins are very complex, spectral deconvolution (grouping peaks into isotopomer envelopes) is a key initial stage for their interpretation. In such spectra, isotopomer envelopes of different protein fragments span overlapping regions on the m/z axis and even share spectral peaks. This raises both pattern recognition and combinatorial challenges for spectral deconvolution. We present MS-Deconv, a combinatorial algorithm for spectral deconvolution. The algorithm first generates a large set of candidate isotopomer envelopes for a spectrum, then represents the spectrum as a graph, and finally selects its highest scoring subset of envelopes as a heaviest path in the graph. In contrast with other approaches, the algorithm scores sets of envelopes rather than individual envelopes. We demonstrate that MS-Deconv improves on Thrash and Xtract in the number of correctly recovered monoisotopic masses and speed. We applied MS-Deconv to a large set of top-down spectra from Yersinia rohdei (with a still unsequenced genome) and further matched them against the protein database of related and sequenced bacterium Yersinia enterocolitica. MS-Deconv is available at http://proteomics.ucsd.edu/Software.html.

A Survey of Computational Methods for Determining Haplotypes

It is widely anticipated that the study of variation in the human genome will provide a means of predicting risk of a variety of complex diseases. Single nucleotide polymorphisms (SNPs) are the most common form of genomic variation. Haplotypes have been suggested as one means for reducing the complexity of studying SNPs. In this paper we review some of the computational approaches that have been taking for determining haplotypes and suggest new approaches.

A unique Mycobacterium ESX-1 protein co-secretes with CFP-10/ESAT-6 and is necessary for inhibiting phagosome maturation.

The ESX‐1 secretion system plays a critical role in the virulence of Mycobacterium tuberculosis and M. marinum. To date, three proteins are known to be secreted by ESX‐1 and necessary for virulence, two of which are CFP‐10 and ESAT‐6. The ESX‐1 secretion and the virulence mechanisms are not well understood. In this study, we have examined the M. marinum secretomes and identified four proteins specific to ESX‐1. Two of those are CFP‐10 and ESAT‐6, and the other two are novel: MM1553 (homologous to Rv3483c) and Mh3881c (homologous to Rv3881c). We have shown that Mh3881c, CFP‐10 and ESAT‐6 are co‐dependent for secretion. Mh3881c is being cleaved at close to the C‐terminus during secretion, and the C‐terminal portion is critical to the co‐dependent secretion, the ESAT‐6 cellular levels, and interaction with ESAT‐6. The co‐dependent secretion is required for M. marinum intracellular growth in macrophages, where the Mh3881c C‐terminal portion plays a critical role. The role of the co‐dependent secretion in intracellular growth correlates with its role in inhibiting phagosome maturation. Both the secretion and the virulence defects of the Mh3881c mutant are complemented by Mh3881c or its M. tuberculosis homologue Rv3881c, suggesting that in M. tuberculosis, Rv3881c has similar functions.

Semi-Automated Identification of N-Glycopeptides by Hydrophilic Interaction Chromatography, nano-Reverse-Phase LC-MS/MS, and Glycan Database Search

Glycoproteins fulfill many indispensable biological functions, and changes in protein glycosylation have been observed in various diseases. Improved analytical methods are needed to allow a complete characterization of this complex and common post-translational modification. In this study, we present a workflow for the analysis of the microheterogeneity of N-glycoproteins that couples hydrophilic interaction and nanoreverse-phase C18 chromatography to tandem QTOF mass spectrometric analysis. A glycan database search program, GlycoPeptideSearch, was developed to match N-glycopeptide MS/MS spectra with the glycopeptides comprised of a glycan drawn from the GlycomeDB glycan structure database and a peptide from a user-specified set of potentially glycosylated peptides. Application of the workflow to human haptoglobin and hemopexin, two microheterogeneous N-glycoproteins, identified a total of 57 distinct site-specific glycoforms in the case of haptoglobin and 14 site-specific glycoforms of hemopexin. Using glycan oxonium ions and peptide-characteristic glycopeptide fragment ions and by collapsing topologically redundant glycans, the search software was able to make unique N-glycopeptide assignments for 51% of assigned spectra, with the remaining assignments primarily representing isobaric topological rearrangements. The optimized workflow, coupled with GlycoPeptideSearch, is expected to make high-throughput semiautomated glycopeptide identification feasible for a wide range of users.

High-throughput middle-down analysis using an orbitrap.

This report demonstrates the application of a capillary LC-LTQ-orbitrap system to provide automated middle-down analysis of proteolytic peptides in the mass range 3000 to 10,000 Da. The novel workflow combines an underutilized method in the orbitrap-high resolution, mass-accurate product ion measurements-with software tailored to search such data (ProSightPC 2.0) and an Asp-selective chemical cleavage approach that generates peptides across an extended mass range. The strategy using high resolution mass measurements on both precursor and product ions is analogous to that widely used on FT-ICR analyzers. The approach is demonstrated in an analysis of the highly basic ribosomal proteome isolated from human MCF7 cancer cells.

Top-Down Identification of Protein Biomarkers in Bacteria with Unsequenced Genomes

MALDI mass spectrometry-based systems for rapid characterization of microorganisms in biodefense or medical diagnostics usually detect intact proteins in the 5000-20,000 Da range. To evaluate the reliability of species discrimination, and also for forensic applications, it is important that these biomarker proteins be identified. In the present study we apply high resolution tandem mass analysis on an Orbitrap and top-down bioinformatics to identify major biomarker proteins observed in MALDI spectra of intact bacteria for which little genomic or protein sequence information is available. The strategy depends on recognition of proteins with very high homology in related (sequenced) species, making it possible to place unsequenced organisms in their correct phylogenetic context. We show that this rapid proteomics based approach to phylogenetic characterization produces similar results to the traditional techniques, and may even be applied to target organisms of undetermined taxonomy. We further discuss important issues in combining genomics/proteomics databases and MALDI MS for the rapid characterization of microorganisms.

Site-specific Glycoforms of Haptoglobin in Liver Cirrhosis and Hepatocellular Carcinoma

Haptoglobin is a liver-secreted glycoprotein with four N-glycosylation sites. Its glycosylation was reported to change in several cancer diseases, which prompted us to examine site-specific glycoforms of haptoglobin in liver cirrhosis and hepatocellular carcinoma. To this end, we have used two-dimensional separation composed of hydrophilic interaction and nano-reverse phase chromatography coupled to QTOF mass spectrometry of the enriched glycopeptides. Our results show increased fucosylation of haptoglobin in liver disease with up to six fucoses associated with specific glycoforms of one glycopeptide. Structural analysis using exoglycosidase treatment and MALDI-MS/MS of detached permethylated glycans led to the identification of Lewis Y-type structures observed particularly in the pooled hepatocellular carcinoma sample. To confirm the increase of the Lewis Y structures observed by LC-MS, we have used immunoaffinity detection with Lewis Y-specific antibodies. The presence of multiply fucosylated Lewis Y glycoforms of haptoglobin in the disease context could have important functional implications.