Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Integrated Proteogenomic Characterization across Major Histological Types of Pediatric Brain Cancer.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

A new molecular signature method for prediction of driver cancer pathways from transcriptional data.

Assigning cancer patients to the most effective treatments requires an understanding of the molecular basis of their disease. While DNA-based molecular profiling approaches have flourished over the past several years to transform our understanding of driver pathways across a broad range of tumors, a systematic characterization of key driver pathways based on RNA data has not been undertaken. Here we introduce a new approach for predicting the status of driver cancer pathways based on signature functions derived from RNA sequencing data. To identify the driver cancer pathways of interest, we mined DNA variant data from TCGA and nominated driver alterations in seven major cancer pathways in breast, ovarian and colon cancer tumors. The activation status of these driver pathways were then characterized using RNA sequencing data by constructing classification signature functions in training datasets and then testing the accuracy of the signatures in test datasets. The signature functions differentiate well tumors with nominated pathway activation from tumors with no signs of activation: average AUC equals to 0.83. Our results confirm that driver genomic alterations are distinctively displayed at the transcriptional level and that the transcriptional signatures can generally provide an alternative to DNA sequencing methods in detecting specific driver pathways.

Genomic Analysis of Uterine Lavage Fluid Detects Early Endometrial Cancers and Reveals a Prevalent Landscape of Driver Mutations in Women without Histopathologic Evidence of Cancer: A Prospective Cross-Sectional Study.

BACKGROUND: Endometrial cancer is the most common gynecologic malignancy, and its incidence and associated mortality are increasing. Despite the immediate need to detect these cancers at an earlier stage, there is no effective screening methodology or protocol for endometrial cancer. The comprehensive, genomics-based analysis of endometrial cancer by The Cancer Genome Atlas (TCGA) revealed many of the molecular defects that define this cancer. Based on these cancer genome results, and in a prospective study, we hypothesized that the use of ultra-deep, targeted gene sequencing could detect somatic mutations in uterine lavage fluid obtained from women undergoing hysteroscopy as a means of molecular screening and diagnosis. METHODS AND FINDINGS: Uterine lavage and paired blood samples were collected and analyzed from 107 consecutive patients who were undergoing hysteroscopy and curettage for diagnostic evaluation from this single-institution study. The lavage fluid was separated into cellular and acellular fractions by centrifugation. Cellular and cell-free DNA (cfDNA) were isolated from each lavage. Two targeted next-generation sequencing (NGS) gene panels, one composed of 56 genes and the other of 12 genes, were used for ultra-deep sequencing. To rule out potential NGS-based errors, orthogonal mutation validation was performed using digital PCR and Sanger sequencing. Seven patients were diagnosed with endometrial cancer based on classic histopathologic analysis. Six of these patients had stage IA cancer, and one of these cancers was only detectable as a microscopic focus within a polyp. All seven patients were found to have significant cancer-associated gene mutations in both cell pellet and cfDNA fractions. In the four patients in whom adequate tumor sample was available, all tumor mutations above a specific allele fraction were present in the uterine lavage DNA samples. Mutations originally only detected in lavage fluid fractions were later confirmed to be present in tumor but at allele fractions significantly less than 1%. Of the remaining 95 patients diagnosed with benign or non-cancer pathology, 44 had no significant cancer mutations detected. Intriguingly, 51 patients without histopathologic evidence of cancer had relatively high allele fraction (1.0%-30.4%), cancer-associated mutations. Participants with detected driver and potential driver mutations were significantly older (mean age mutated = 57.96, 95% confidence interval [CI]: 3.30-∞, mean age no mutations = 50.35; p-value = 0.002; Benjamini-Hochberg [BH] adjusted p-value = 0.015) and more likely to be post-menopausal (p-value = 0.004; BH-adjusted p-value = 0.015) than those without these mutations. No associations were detected between mutation status and race/ethnicity, body mass index, diabetes, parity, and smoking status. Long-term follow-up was not presently available in this prospective study for those women without histopathologic evidence of cancer. CONCLUSIONS: Using ultra-deep NGS, we identified somatic mutations in DNA extracted both from cell pellets and a never previously reported cfDNA fraction from the uterine lavage. Using our targeted sequencing approach, endometrial driver mutations were identified in all seven women who received a cancer diagnosis based on classic histopathology of tissue curettage obtained at the time of hysteroscopy. In addition, relatively high allele fraction driver mutations were identified in the lavage fluid of approximately half of the women without a cancer diagnosis. Increasing age and post-menopausal status were associated with the presence of these cancer-associated mutations, suggesting the prevalent existence of a premalignant landscape in women without clinical evidence of cancer. Given that a uterine lavage can be easily and quickly performed even outside of the operating room and in a physician's office-based setting, our findings suggest the future possibility of this approach for screening women for the earliest stages of endometrial cancer. However, our findings suggest that further insight into development of cancer or its interruption are needed before translation to the clinic.

Multi-scale signaling and tumor evolution in high-grade gliomas.

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Association between Incidental Pelvic Inflammation and Aggressive Prostate Cancer.

The impact of pelvic inflammation on prostate cancer (PCa) biology and aggressive phenotype has never been studied. Our study objective was to evaluate the role of pelvic inflammation on PCa aggressiveness and its association with clinical outcomes in patients following radical prostatectomy (RP). This study has been conducted on a retrospective single-institutional consecutive cohort of 2278 patients who underwent robot-assisted laparoscopic prostatectomy (RALP) between 01/2013 and 10/2019. Data from 2085 patients were analyzed to study the association between pelvic inflammation and adverse pathology (AP), defined as Gleason Grade Group (GGG) > 2 and ≥ pT3 stage, at resection. In a subset of 1997 patients, the association between pelvic inflammation and biochemical recurrence (BCR) was studied. Alteration in tumor transcriptome and inflammatory markers in patients with and without pelvic inflammation were studied using microarray analysis, immunohistochemistry, and culture supernatants derived from inflamed sites used in functional assays. Changes in blood inflammatory markers in the study cohort were analyzed by O-link. In univariate analyses, pelvic inflammation emerged as a significant predictor of AP. Multivariate cox proportional-hazards regression analyses showed that high pelvic inflammation with pT3 stage and positive surgical margins significantly affected the time to BCR (p ≤ 0.05). PCa patients with high inflammation had elevated levels of pro-inflammatory cytokines in their tissues and in blood. Genes involved in epithelial-to-mesenchymal transition (EMT) and DNA damage response were upregulated in patients with pelvic inflammation. Attenuation of STAT and IL-6 signaling decreased tumor driving properties of conditioned medium from inflamed sites. Pelvic inflammation exacerbates the progression of prostate cancer and drives an aggressive phenotype.

Proteogenomic and metabolomic characterization of human glioblastoma.

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.

New statistical potential for quality assessment of protein models and a survey of energy functions.

BACKGROUND: Scoring functions, such as molecular mechanic forcefields and statistical potentials are fundamentally important tools in protein structure modeling and quality assessment. RESULTS: The performances of a number of publicly available scoring functions are compared with a statistical rigor, with an emphasis on knowledge-based potentials. We explored the effect on accuracy of alternative choices for representing interaction center types and other features of scoring functions, such as using information on solvent accessibility, on torsion angles, accounting for secondary structure preferences and side chain orientation. Partially based on the observations made, we present a novel residue based statistical potential, which employs a shuffled reference state definition and takes into account the mutual orientation of residue side chains. Atom- and residue-level statistical potentials and Linux executables to calculate the energy of a given protein proposed in this work can be downloaded from http://www.fiserlab.org/potentials. CONCLUSIONS: Among the most influential terms we observed a critical role of a proper reference state definition and the benefits of including information about the microenvironment of interaction centers. Molecular mechanical potentials were also tested and found to be over-sensitive to small local imperfections in a structure, requiring unfeasible long energy relaxation before energy scores started to correlate with model quality.

Improved scoring function for comparative modeling using the M4T method.

Improvements in comparative protein structure modeling for the remote target-template sequence similarity cases are possible through the optimal combination of multiple template structures and by improving the quality of target-template alignment. Recently developed MMM and M4T methods were designed to address these problems. Here we describe new developments in both the alignment generation and the template selection parts of the modeling algorithms. We set up a new scoring function in MMM to deliver more accurate target-template alignments. This was achieved by developing and incorporating into the composite scoring function a novel statistical pairwise potential that combines local and non-local terms. The non-local term of the statistical potential utilizes a shuffled reference state definition that helped to eliminate most of the false positive signal from the background distribution of pairwise contacts. The accuracy of the scoring function was further increased by using BLOSUM mutation table scores.

Nonprotein based enrichment method to analyze peptide cross-linking in protein complexes.

Cross-linking analysis of protein complexes and structures by tandem mass spectrometry (MS/MS) has advantages in speed, sensitivity, specificity, and the capability of handling complicated protein assemblies. However, detection and accurate assignment of the cross-linked peptides are often challenging due to their low abundance and complicated fragmentation behavior in collision-induced dissociation (CID). To simplify the MS analysis and improve the signal-to-noise ratio of the cross-linked peptides, we developed a novel peptide enrichment strategy that utilizes a cross-linker with a cryptic thiol group and using beads modified with a photocleavable cross-linker. The functional cross-linkers were designed to react with the primary amino groups in proteins. Human serum albumin was used as a model protein to detect intra- and intermolecular cross-linkages. Use of this protein-free selective retrieval method eliminates the contamination that can result from avidin-biotin based retrieval systems and simplifies data analysis. These features may make the method suitable to investigate protein-protein interactions in biological samples.

In vitro evaluation, biodistribution and scintigraphic imaging in mice of radiolabeled anthrax toxins.

INTRODUCTION: There is a lot of interest towards creating therapies and vaccines for Bacillus anthracis, a bacterium which causes anthrax in humans and which spores can be made into potent biological weapons. Systemic injection of lethal factor (LF), edema factor (EF) and protective antigen (PA) in mice produces toxicity, and this protocol is commonly used to investigate the efficacy of specific antibodies in passive protection and vaccine studies. Availability of toxins labeled with imageable radioisotopes would allow to demonstrate their tissue distribution after intravenous injection at toxin concentration that are below pharmacologically significant to avoid masking by toxic effects. METHODS: LF, EF and PA were radiolabeled with (188)Re and (99m)Tc, and their performance in vitro was evaluated by macrophages and Chinese hamster ovary cells toxicity assays and by binding to macrophages. Scintigraphic imaging and biodistribution of intravenously (IV) injected (99m)Tc-and (123)I-labeled toxins was performed in BALB/c mice. RESULTS: Radiolabeled toxins preserved their biological activity. Scatchard-type analysis of the binding of radiolabeled PA to the J774.16 macrophage-like cells revealed 6.6 x 10(4) binding sites per cell with a dissociation constant of 6.7 nM. Comparative scintigraphic imaging of mice injected intravenously with either (99m)Tc-or (123)I-labeled PA, EF and LF toxins demonstrated similar biodistribution patterns with early localization of radioactivity in the liver, spleen, intestines and excretion through kidneys. The finding of renal excretion shortly after IV injection strongly suggests that toxins are rapidly degraded which could contribute to the variability of mouse toxigenic assays. Biodistribution studies confirmed that all three toxins concentrated in the liver and the presence of high levels of radioactivity again implied rapid degradation in vivo. CONCLUSIONS: The availability of (188)Re and (99m)Tc-labeled PA, LF and EF toxins allowed us to confirm the number of PA binding sites per cell, to provide an estimate of the dissociation constant of PA for its receptor and to demonstrate tissue distribution of toxins in mice after intravenous injection.

Effects of amino acid composition, finite size of proteins, and sparse statistics on distance-dependent statistical pair potentials.

Statistical distance dependent pair potentials are frequently used in a variety of folding, threading, and modeling studies of proteins. The applicability of these types of potentials is tightly connected to the reliability of statistical observations. We explored the possible origin and extent of false positive signals in statistical potentials by analyzing their distance dependence in a variety of randomized protein-like models. While on average potentials derived from such models are expected to equal zero at any distance, we demonstrate that systematic and significant distortions exist. These distortions originate from the limited statistical counts in local environments of proteins and from the limited size of protein structures at large distances. We suggest that these systematic errors in statistical potentials are connected to the dependence of amino acid composition on protein size and to variation in protein sizes. Additionally, atom-based potentials are dominated by a false positive signal that is due to correlation among distances measured from atoms of one residue to atoms of another residue. The significance of residue-based pairwise potentials at various spatial pair separations was assessed in this study and it was found that as few as approximately 50% of potential values were statistically significant at distances below 4 A, and only at most approximately 80% of them were significant at larger pair separations. A new definition for reference state, free of the observed systematic errors, is suggested. It has been demonstrated to generate statistical potentials that compare favorably to other publicly available ones.

Acetylsalicylic acid (aspirin) reduces damage to reconstituted human tissues infected with Candida species by inhibiting extracellular fungal lipases.

A reconstituted human tissue model was used to mimic Candida albicans and Candida parapsilosis infection in order to investigate the protective effects of acetylsalicylic acid (aspirin, ASA). We found that therapeutic concentrations of ASA reduced tissue damage in the in vitro infection model. We further evaluated the lipase inhibitory effects of ASA by investigating the growth of C. albicans, C. parapsilosis and C. parapsilosis lipase negative (Deltacplip1-2/Deltacplip1-2) mutants in a lipid rich minimal medium supplemented with olive oil and found that a therapeutic concentration of ASA inhibited the growth of wild type fungi. The lipase inhibitors quinine and ebelactone B were also shown to reduce growth and protect against tissue damage from Candida species, respectively. A lipolytic activity assay also showed that therapeutic concentrations of ASA inhibited C. antarctica and C. cylindracea purified lipases obtained through a commercial kit. The relationship between ASA and lipase was characterized through a computed structural model of the Lipase-2 protein from C. parapsilosis in complex with ASA. The results suggest that development of inhibitors of fungal lipases could result in broad-spectrum therapeutics, especially since fungal lipases are not homologous to their human analogues.

Computational analysis and experimental validation of gene predictions in Toxoplasma gondii.

BACKGROUND: Toxoplasma gondii is an obligate intracellular protozoan that infects 20 to 90% of the population. It can cause both acute and chronic infections, many of which are asymptomatic, and, in immunocompromised hosts, can cause fatal infection due to reactivation from an asymptomatic chronic infection. An essential step towards understanding molecular mechanisms controlling transitions between the various life stages and identifying candidate drug targets is to accurately characterize the T. gondii proteome. METHODOLOGY/PRINCIPAL FINDINGS: We have explored the proteome of T. gondii tachyzoites with high throughput proteomics experiments and by comparison to publicly available cDNA sequence data. Mass spectrometry analysis validated 2,477 gene coding regions with 6,438 possible alternative gene predictions; approximately one third of the T. gondii proteome. The proteomics survey identified 609 proteins that are unique to Toxoplasma as compared to any known species including other Apicomplexan. Computational analysis identified 787 cases of possible gene duplication events and located at least 6,089 gene coding regions. Commonly used gene prediction algorithms produce very disparate sets of protein sequences, with pairwise overlaps ranging from 1.4% to 12%. Through this experimental and computational exercise we benchmarked gene prediction methods and observed false negative rates of 31 to 43%. CONCLUSIONS/SIGNIFICANCE: This study not only provides the largest proteomics exploration of the T. gondii proteome, but illustrates how high throughput proteomics experiments can elucidate correct gene structures in genomes.