Proteogenomic Analysis of Human Colon Cancer Reveals New Therapeutic Opportunities.

We performed the first proteogenomic study on a prospectively collected colon cancer cohort. Comparative proteomic and phosphoproteomic analysis of paired tumor and normal adjacent tissues produced a catalog of colon cancer-associated proteins and phosphosites, including known and putative new biomarkers, drug targets, and cancer/testis antigens. Proteogenomic integration not only prioritized genomically inferred targets, such as copy-number drivers and mutation-derived neoantigens, but also yielded novel findings. Phosphoproteomics data associated Rb phosphorylation with increased proliferation and decreased apoptosis in colon cancer, which explains why this classical tumor suppressor is amplified in colon tumors and suggests a rationale for targeting Rb phosphorylation in colon cancer. Proteomics identified an association between decreased CD8 T cell infiltration and increased glycolysis in microsatellite instability-high (MSI-H) tumors, suggesting glycolysis as a potential target to overcome the resistance of MSI-H tumors to immune checkpoint blockade. Proteogenomics presents new avenues for biological discoveries and therapeutic development.

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.

Mass Spectrometry-Based Proteogenomics: New Therapeutic Opportunities for Precision Medicine.

Proteogenomics refers to the integration of comprehensive genomic, transcriptomic, and proteomic measurements from the same samples with the goal of fully understanding the regulatory processes converting genotypes to phenotypes, often with an emphasis on gaining a deeper understanding of disease processes. Although specific genetic mutations have long been known to drive the development of multiple cancers, gene mutations alone do not always predict prognosis or response to targeted therapy. The benefit of proteogenomics research is that information obtained from proteins and their corresponding pathways provides insight into therapeutic targets that can complement genomic information by providing an additional dimension regarding the underlying mechanisms and pathophysiology of tumors. This review describes the novel insights into tumor biology and drug resistance derived from proteogenomic analysis while highlighting the clinical potential of proteogenomic observations and advances in technique and analysis tools.

Molecular-Level Dysregulation of Insulin Pathways and Inflammatory Processes in Peripheral Blood Mononuclear Cells by Circadian Misalignment.

Circadian misalignment due to night work has been associated with an elevated risk for chronic diseases. We investigated the effects of circadian misalignment using shotgun protein profiling of peripheral blood mononuclear cells taken from healthy humans during a constant routine protocol, which was conducted immediately after participants had been subjected to a 3-day simulated night shift schedule or a 3-day simulated day shift schedule. By comparing proteomic profiles between the simulated shift conditions, we identified proteins and pathways that are associated with the effects of circadian misalignment and observed that insulin regulation pathways and inflammation-related proteins displayed markedly different temporal patterns after simulated night shift. Further, by integrating the proteomic profiles with previously assessed metabolomic profiles in a network-based approach, we found key associations between circadian dysregulation of protein-level pathways and metabolites of interest in the context of chronic metabolic diseases. Endogenous circadian rhythms in circulating glucose and insulin differed between the simulated shift conditions. Overall, our results suggest that circadian misalignment is associated with a tug of war between central clock mechanisms controlling insulin secretion and peripheral clock mechanisms regulating insulin sensitivity, which may lead to adverse long-term outcomes such as diabetes and obesity. Our study provides a molecular-level mechanism linking circadian misalignment and adverse long-term health consequences of night work.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection.

BACKGROUND: Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors. METHODS: We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers. RESULTS: We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors. CONCLUSIONS: This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Proteomic Determinants of Variation in Cholesterol Efflux: Observations from the Dallas Heart Study.

High-density lipoproteins (HDLs) are promising targets for predicting and treating atherosclerotic cardiovascular disease (ASCVD), as they mediate removal of excess cholesterol from lipid-laden macrophages that accumulate in the vasculature. This functional property of HDLs, termed cholesterol efflux capacity (CEC), is inversely associated with ASCVD. HDLs are compositionally diverse, associating with >250 different proteins, but their relative contribution to CEC remains poorly understood. Our goal was to identify and define key HDL-associated proteins that modulate CEC in humans. The proteomic signature of plasma HDL was quantified in 36 individuals in the multi-ethnic population-based Dallas Heart Study (DHS) cohort that exhibited persistent extremely high (>=90th%) or extremely low CEC (<=10th%) over 15 years. Levels of apolipoprotein (Apo)A-I associated ApoC-II, ApoC-III, and ApoA-IV were differentially correlated with CEC in high (r = 0.49, 0.41, and -0.21 respectively) and low (r = -0.46, -0.41, and 0.66 respectively) CEC groups (p for heterogeneity (pHet) = 0.03, 0.04, and 0.003 respectively). Further, we observed that levels of ApoA-I with ApoC-III, complement C3 (CO3), ApoE, and plasminogen (PLMG) were inversely associated with CEC in individuals within the low CEC group (r = -0.11 to -0.25 for subspecies with these proteins vs. r = 0.58 to 0.65 for subspecies lacking these proteins; p < 0.05 for heterogeneity). These findings suggest that enrichment of specific proteins on HDLs and, thus, different subspecies of HDLs, differentially modulate the removal of cholesterol from the vasculature.

Proteomic and phosphoproteomic measurements enhance ability to predict ex vivo drug response in AML.

Acute Myeloid Leukemia (AML) affects 20,000 patients in the US annually with a five-year survival rate of approximately 25%. One reason for the low survival rate is the high prevalence of clonal evolution that gives rise to heterogeneous sub-populations of leukemic cells with diverse mutation spectra, which eventually leads to disease relapse. This genetic heterogeneity drives the activation of complex signaling pathways that is reflected at the protein level. This diversity makes it difficult to treat AML with targeted therapy, requiring custom patient treatment protocols tailored to each individual's leukemia. Toward this end, the Beat AML research program prospectively collected genomic and transcriptomic data from over 1000 AML patients and carried out ex vivo drug sensitivity assays to identify genomic signatures that could predict patient-specific drug responses. However, there are inherent weaknesses in using only genetic and transcriptomic measurements as surrogates of drug response, particularly the absence of direct information about phosphorylation-mediated signal transduction. As a member of the Clinical Proteomic Tumor Analysis Consortium, we have extended the molecular characterization of this cohort by collecting proteomic and phosphoproteomic measurements from a subset of these patient samples (38 in total) to evaluate the hypothesis that proteomic signatures can improve the ability to predict response to 26 drugs in AML ex vivo samples. In this work we describe our systematic, multi-omic approach to evaluate proteomic signatures of drug response and compare protein levels to other markers of drug response such as mutational patterns. We explore the nuances of this approach using two drugs that target key pathways activated in AML: quizartinib (FLT3) and trametinib (Ras/MEK), and show how patient-derived signatures can be interpreted biologically and validated in cell lines. In conclusion, this pilot study demonstrates strong promise for proteomics-based patient stratification to assess drug sensitivity in AML.

Hypergraph models of biological networks to identify genes critical to pathogenic viral response.

BACKGROUND: Representing biological networks as graphs is a powerful approach to reveal underlying patterns, signatures, and critical components from high-throughput biomolecular data. However, graphs do not natively capture the multi-way relationships present among genes and proteins in biological systems. Hypergraphs are generalizations of graphs that naturally model multi-way relationships and have shown promise in modeling systems such as protein complexes and metabolic reactions. In this paper we seek to understand how hypergraphs can more faithfully identify, and potentially predict, important genes based on complex relationships inferred from genomic expression data sets. RESULTS: We compiled a novel data set of transcriptional host response to pathogenic viral infections and formulated relationships between genes as a hypergraph where hyperedges represent significantly perturbed genes, and vertices represent individual biological samples with specific experimental conditions. We find that hypergraph betweenness centrality is a superior method for identification of genes important to viral response when compared with graph centrality. CONCLUSIONS: Our results demonstrate the utility of using hypergraphs to represent complex biological systems and highlight central important responses in common to a variety of highly pathogenic viruses.

leapR: An R Package for Multiomic Pathway Analysis.

A generalized goal of many high-throughput data studies is to identify functional mechanisms that underlie observed biological phenomena, whether they be disease outcomes or metabolic output. Increasingly, studies that rely on multiple sources of high-throughput data (genomic, transcriptomic, proteomic, metabolomic) are faced with a challenge of summarizing the data to generate testable hypotheses. However, this requires a time-consuming process to evaluate numerous statistical methods across numerous data sources. Here, we introduce the leapR package, a framework to rapidly assess biological pathway activity using diverse statistical tests and data sources, allowing facile integration of multisource data. The leapR package with a user manual and example workflow is available for download from GitHub (https://github.com/biodataganache/leapR).

Night shift schedule causes circadian dysregulation of DNA repair genes and elevated DNA damage in humans.

Circadian disruption has been identified as a risk factor for health disorders such as obesity, cardiovascular disease, and cancer. Although epidemiological studies suggest an increased risk of various cancers associated with circadian misalignment due to night shift work, the underlying mechanisms have yet to be elucidated. We sought to investigate the potential mechanistic role that circadian disruption of cancer hallmark pathway genes may play in the increased cancer risk in shift workers. In a controlled laboratory study, we investigated the circadian transcriptome of cancer hallmark pathway genes and associated biological pathways in circulating leukocytes obtained from healthy young adults during a 24-hour constant routine protocol following 3 days of simulated day shift or night shift. The simulated night shift schedule significantly altered the normal circadian rhythmicity of genes involved in cancer hallmark pathways. A DNA repair pathway showed significant enrichment of rhythmic genes following the simulated day shift schedule, but not following the simulated night shift schedule. In functional assessments, we demonstrated that there was an increased sensitivity to both endogenous and exogenous sources of DNA damage after exposure to simulated night shift. Our results suggest that circadian dysregulation of DNA repair may increase DNA damage and potentiate elevated cancer risk in night shift workers.

An Integrative Analysis of Tumor Proteomic and Phosphoproteomic Profiles to Examine the Relationships Between Kinase Activity and Phosphorylation.

Phosphorylation of proteins is a key way cells regulate function, both at the individual protein level and at the level of signaling pathways. Kinases are responsible for phosphorylation of substrates, generally on serine, threonine, or tyrosine residues. Though particular sequence patterns can be identified that dictate whether a residue will be phosphorylated by a specific kinase, these patterns are not highly predictive of phosphorylation. The availability of large scale proteomic and phosphoproteomic data sets generated using mass-spectrometry-based approaches provides an opportunity to study the important relationship between kinase activity, substrate specificity, and phosphorylation. In this study, we analyze relationships between protein abundance and phosphopeptide abundance across more than 150 tumor samples and show that phosphorylation at specific phosphosites is not well correlated with overall kinase abundance. However, individual kinases show a clear and statistically significant difference in correlation among known phosphosite targets for that kinase and randomly selected phosphosites. We further investigate relationships between phosphorylation of known activating or inhibitory sites on kinases and phosphorylation of their target phosphosites. Combined with motif-based analysis, this approach can predict novel kinase targets and show which subsets of a kinase's target repertoire are specifically active in one condition versus another.

Phosphoproteome Analysis Reveals Estrogen-ER Pathway as a Modulator of mTOR Activity Via DEPTOR.

ER-positive breast tumors represent ∼70% of all breast cancer cases. Although their treatment with endocrine therapies is effective in the adjuvant or recurrent settings, the development of resistance compromises their effectiveness. The binding of estrogen to ERα, a transcription factor, triggers the regulation of the target genes (genomic pathway). Additionally, a cytoplasmic fraction of estrogen-bound ERα activates oncogenic signaling pathways such as PI3K/AKT/mTOR (nongenomic pathway). The upregulation of the estrogenic and the PI3K/AKT/mTOR signaling pathways are frequently associated with a poor outcome. To better characterize the connection between these two pathways, we performed a phosphoproteome analysis of ER-positive MCF7 breast cancer cells treated with estrogen or estrogen and the mTORC1 inhibitor rapamycin. Many proteins were identified as estrogen-regulated mTORC1 targets and among them, DEPTOR was selected for further characterization. DEPTOR binds to mTOR and inhibits the kinase activity of both mTOR complexes mTORC1 and mTORC2, but mitogen-activated mTOR promotes phosphorylation-mediated DEPTOR degradation. Although estrogen enhances the phosphorylation of DEPTOR by mTORC1, DEPTOR levels increase in estrogen-stimulated cells. We demonstrated that DEPTOR accumulation is the result of estrogen-ERα-mediated transcriptional upregulation of DEPTOR expression. Consequently, the elevated levels of DEPTOR partially counterbalance the estrogen-induced activation of mTORC1 and mTORC2. These results underscore the critical role of estrogen-ERα as a modulator of the PI3K/AKT/mTOR signaling pathway in ER-positive breast cancer cells. Additionally, these studies provide evidence supporting the use of dual PI3K/mTOR or dual mTORC1/2 inhibitors in combination with endocrine therapies as a first-line treatment option for the patients with ER-positive advanced breast cancer.

Unified feature association networks through integration of transcriptomic and proteomic data.

High-throughput multi-omics studies and corresponding network analyses of multi-omic data have rapidly expanded their impact over the last 10 years. As biological features of different types (e.g. transcripts, proteins, metabolites) interact within cellular systems, the greatest amount of knowledge can be gained from networks that incorporate multiple types of -omic data. However, biological and technical sources of variation diminish the ability to detect cross-type associations, yielding networks dominated by communities comprised of nodes of the same type. We describe here network building methods that can maximize edges between nodes of different data types leading to integrated networks, networks that have a large number of edges that link nodes of different-omic types (transcripts, proteins, lipids etc). We systematically rank several network inference methods and demonstrate that, in many cases, using a random forest method, GENIE3, produces the most integrated networks. This increase in integration does not come at the cost of accuracy as GENIE3 produces networks of approximately the same quality as the other network inference methods tested here. Using GENIE3, we also infer networks representing antibody-mediated Dengue virus cell invasion and receptor-mediated Dengue virus invasion. A number of functional pathways showed centrality differences between the two networks including genes responding to both GM-CSF and IL-4, which had a higher centrality value in an antibody-mediated vs. receptor-mediated Dengue network. Because a biological system involves the interplay of many different types of molecules, incorporating multiple data types into networks will improve their use as models of biological systems. The methods explored here are some of the first to specifically highlight and address the challenges associated with how such multi-omic networks can be assembled and how the greatest number of interactions can be inferred from different data types. The resulting networks can lead to the discovery of new host response patterns and interactions during viral infection, generate new hypotheses of pathogenic mechanisms and confirm mechanisms of disease.

Proteome Profiling Outperforms Transcriptome Profiling for Coexpression Based Gene Function Prediction.

Coexpression of mRNAs under multiple conditions is commonly used to infer cofunctionality of their gene products despite well-known limitations of this "guilt-by-association" (GBA) approach. Recent advancements in mass spectrometry-based proteomic technologies have enabled global expression profiling at the protein level; however, whether proteome profiling data can outperform transcriptome profiling data for coexpression based gene function prediction has not been systematically investigated. Here, we address this question by constructing and analyzing mRNA and protein coexpression networks for three cancer types with matched mRNA and protein profiling data from The Cancer Genome Atlas (TCGA) and the Clinical Proteomic Tumor Analysis Consortium (CPTAC). Our analyses revealed a marked difference in wiring between the mRNA and protein coexpression networks. Whereas protein coexpression was driven primarily by functional similarity between coexpressed genes, mRNA coexpression was driven by both cofunction and chromosomal colocalization of the genes. Functionally coherent mRNA modules were more likely to have their edges preserved in corresponding protein networks than functionally incoherent mRNA modules. Proteomic data strengthened the link between gene expression and function for at least 75% of Gene Ontology (GO) biological processes and 90% of KEGG pathways. A web application Gene2Net (http://cptac.gene2net.org) developed based on the three protein coexpression networks revealed novel gene-function relationships, such as linking ERBB2 (HER2) to lipid biosynthetic process in breast cancer, identifying PLG as a new gene involved in complement activation, and identifying AEBP1 as a new epithelial-mesenchymal transition (EMT) marker. Our results demonstrate that proteome profiling outperforms transcriptome profiling for coexpression based gene function prediction. Proteomics should be integrated if not preferred in gene function and human disease studies.

Improving network inference algorithms using resampling methods.

BACKGROUND: Relatively small changes to gene expression data dramatically affect co-expression networks inferred from that data which, in turn, can significantly alter the subsequent biological interpretation. This error propagation is an underappreciated problem that, while hinted at in the literature, has not yet been thoroughly explored. Resampling methods (e.g. bootstrap aggregation, random subspace method) are hypothesized to alleviate variability in network inference methods by minimizing outlier effects and distilling persistent associations in the data. But the efficacy of the approach assumes the generalization from statistical theory holds true in biological network inference applications. RESULTS: We evaluated the effect of bootstrap aggregation on inferred networks using commonly applied network inference methods in terms of stability, or resilience to perturbations in the underlying expression data, a metric for accuracy, and functional enrichment of edge interactions. CONCLUSION: Bootstrap aggregation results in improved stability and, depending on the size of the input dataset, a marginal improvement to accuracy assessed by each method's ability to link genes in the same functional pathway.

Time-resolved proteome profiling of normal lung development.

Biochemical networks mediating normal lung morphogenesis and function have important implications for ameliorating morbidity and mortality in premature infants. Although several transcript-level studies have examined normal lung development, corresponding protein-level analyses are lacking. Here we performed proteomics analysis of murine lungs from embryonic to early adult ages to identify the molecular networks mediating normal lung development. We identified 8,932 proteins, providing a deep and comprehensive view of the lung proteome. Analysis of the proteomics data revealed discrete modules and the underlying regulatory and signaling network modulating their expression during development. Our data support the cell proliferation that characterizes early lung development and highlight responses of the lung to exposure to a nonsterile oxygen-rich ambient environment and the important role of lipid (surfactant) metabolism in lung development. Comparison of dynamic regulation of proteomic and recent transcriptomic analyses identified biological processes under posttranscriptional control. Our study provides a unique proteomic resource for understanding normal lung formation and function and can be freely accessed at Lungmap.net.

Network analysis of transcriptomics expands regulatory landscapes in Synechococcus sp. PCC 7002.

Cyanobacterial regulation of gene expression must contend with a genome organization that lacks apparent functional context, as the majority of cellular processes and metabolic pathways are encoded by genes found at disparate locations across the genome and relatively few transcription factors exist. In this study, global transcript abundance data from the model cyanobacterium Synechococcus sp. PCC 7002 grown under 42 different conditions was analyzed using Context-Likelihood of Relatedness (CLR). The resulting network, organized into 11 modules, provided insight into transcriptional network topology as well as grouping genes by function and linking their response to specific environmental variables. When used in conjunction with genome sequences, the network allowed identification and expansion of novel potential targets of both DNA binding proteins and sRNA regulators. These results offer a new perspective into the multi-level regulation that governs cellular adaptations of the fast-growing physiologically robust cyanobacterium Synechococcus sp. PCC 7002 to changing environmental variables. It also provides a methodological high-throughput approach to studying multi-scale regulatory mechanisms that operate in cyanobacteria. Finally, it provides valuable context for integrating systems-level data to enhance gene grouping based on annotated function, especially in organisms where traditional context analyses cannot be implemented due to lack of operon-based functional organization.

Reproducibility of Differential Proteomic Technologies in CPTAC Fractionated Xenografts.

The NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) employed a pair of reference xenograft proteomes for initial platform validation and ongoing quality control of its data collection for The Cancer Genome Atlas (TCGA) tumors. These two xenografts, representing basal and luminal-B human breast cancer, were fractionated and analyzed on six mass spectrometers in a total of 46 replicates divided between iTRAQ and label-free technologies, spanning a total of 1095 LC-MS/MS experiments. These data represent a unique opportunity to evaluate the stability of proteomic differentiation by mass spectrometry over many months of time for individual instruments or across instruments running dissimilar workflows. We evaluated iTRAQ reporter ions, label-free spectral counts, and label-free extracted ion chromatograms as strategies for data interpretation (source code is available from http://homepages.uc.edu/~wang2x7/Research.htm ). From these assessments, we found that differential genes from a single replicate were confirmed by other replicates on the same instrument from 61 to 93% of the time. When comparing across different instruments and quantitative technologies, using multiple replicates, differential genes were reproduced by other data sets from 67 to 99% of the time. Projecting gene differences to biological pathways and networks increased the degree of similarity. These overlaps send an encouraging message about the maturity of technologies for proteomic differentiation.

The landscape of viral proteomics and its potential to impact human health.

INTRODUCTION: Advances in mass spectrometry-based proteomic technologies are enhancing studies of viral pathogenesis. Identification and quantification of host and viral proteins and modifications in cells and extracellular fluids during infection provides useful information about pathogenesis, and will be critical for directing clinical interventions and diagnostics. AREAS COVERED: Herein we review and discuss a broad range of global proteomic studies conducted during viral infection, including those of cellular responses, protein modifications, virion packaging, and serum proteomics. We focus on viruses that impact human health and focus on experimental designs that reveal disease processes and surrogate markers. Expert commentary: Global proteomics is an important component of systems-level studies that aim to define how the interaction of humans and viruses leads to disease. Viral-community resource centers and strategies from other fields (e.g., cancer) will facilitate data sharing and platform-integration for systems-level analyses, and should provide recommended standards and assays for experimental designs and validation.