Predicting immune checkpoint therapy response in three independent metastatic melanoma cohorts.

While Immune checkpoint inhibition (ICI) therapy shows significant efficacy in metastatic melanoma, only about 50% respond, lacking reliable predictive methods. We introduce a panel of six proteins aimed at predicting response to ICI therapy. Evaluating previously reported proteins in two untreated melanoma cohorts, we used a published predictive model (EaSIeR score) to identify potential proteins distinguishing responders and non-responders. Six proteins initially identified in the ICI cohort correlated with predicted response in the untreated cohort. Additionally, three proteins correlated with patient survival, both at the protein, and at the transcript levels, in an independent immunotherapy treated cohort. Our study identifies predictive biomarkers across three melanoma cohorts, suggesting their use in therapeutic decision-making.

The MTIST platform: a microbiome time series inference standardized test.

The human gut microbiome is a promising therapeutic target, but interventions are hampered by our limited understanding of microbial ecosystems. Here, we present a platform to develop, evaluate, and score approaches to learn ecological interactions from microbiome time series data. The microbiome time series inference standardized test (MTIST) comprises: a simulation framework for the in silico generation of microbiome study data akin to what is obtained with quantitative next-generation sequencing approaches, a compilation of a large curated data set generated by the simulation framework representing 648 simulated microbiome studies containing 18,360 time series, with a total of 2,182,800 species abundance measurements, and a scoring method to rank ecological inference algorithms. We use the MTIST platform to rank five implementations of microbiome inference approaches, revealing that while all algorithms performed well on ecosystems with few species (3 and 10), all algorithms failed to infer most interaction in a large ecosystem with 100 member species. However, we do find that the strongest interactions within a large ecosystem are inferred with higher success by all algorithms. Finally, we use the MTIST platform to compare different microbiome study designs, characterizing tradeoffs between samples per subject and number of subjects. Interestingly, we find that when only few samples can be collected per subject, ecological inference is most successful when these samples are collected with highest feasible temporal frequency. Taken together, we provide a computational tool to aid the development of better microbiome ecosystem inference approaches, which will be crucial towards the development of reliable and predictable therapeutic approaches that target the microbiome ecosystem.

Modulation of GPR133 (ADGRD1) signaling by its intracellular interaction partner extended synaptotagmin 1.

GPR133 (ADGRD1) is an adhesion G-protein-coupled receptor that signals through Gαs/cyclic AMP (cAMP) and is required for the growth of glioblastoma (GBM), an aggressive brain malignancy. The regulation of GPR133 signaling is incompletely understood. Here, we use proximity biotinylation proteomics to identify ESYT1, a Ca2+-dependent mediator of endoplasmic reticulum-plasma membrane bridge formation, as an intracellular interactor of GPR133. ESYT1 knockdown or knockout increases GPR133 signaling, while its overexpression has the opposite effect, without altering GPR133 levels in the plasma membrane. The GPR133-ESYT1 interaction requires the Ca2+-sensing C2C domain of ESYT1. Thapsigargin-mediated increases in cytosolic Ca2+ relieve signaling-suppressive effects of ESYT1 by promoting ESYT1-GPR133 dissociation. ESYT1 knockdown or knockout in GBM slows tumor growth, suggesting tumorigenic functions of ESYT1. Our findings demonstrate a mechanism for the modulation of GPR133 signaling by increased cytosolic Ca2+, which reduces the signaling-suppressive interaction between GPR133 and ESYT1 to raise cAMP levels.

Few-shot genes selection: subset of PAM50 genes for breast cancer subtypes classification.

BACKGROUND: In recent years, researchers have made significant strides in understanding the heterogeneity of breast cancer and its various subtypes. However, the wealth of genomic and proteomic data available today necessitates efficient frameworks, instruments, and computational tools for meaningful analysis. Despite its success as a prognostic tool, the PAM50 gene signature's reliance on many genes presents challenges in terms of cost and complexity. Consequently, there is a need for more efficient methods to classify breast cancer subtypes using a reduced gene set accurately. RESULTS: This study explores the potential of achieving precise breast cancer subtype categorization using a reduced gene set derived from the PAM50 gene signature. By employing a "Few-Shot Genes Selection" method, we randomly select smaller subsets from PAM50 and evaluate their performance using metrics and a linear model, specifically the Support Vector Machine (SVM) classifier. In addition, we aim to assess whether a more compact gene set can maintain performance while simplifying the classification process. Our findings demonstrate that certain reduced gene subsets can perform comparable or superior to the full PAM50 gene signature. CONCLUSIONS: The identified gene subsets, with 36 genes, have the potential to contribute to the development of more cost-effective and streamlined diagnostic tools in breast cancer research and clinical settings.

CDK-independent role of D-type cyclins in regulating DNA mismatch repair.

Although mismatch repair (MMR) is essential for correcting DNA replication errors, it can also recognize other lesions, such as oxidized bases. In G0 and G1, MMR is kept in check through unknown mechanisms as it is error-prone during these cell cycle phases. We show that in mammalian cells, D-type cyclins are recruited to sites of oxidative DNA damage in a PCNA- and p21-dependent manner. D-type cyclins inhibit the proteasomal degradation of p21, which competes with MMR proteins for binding to PCNA, thereby inhibiting MMR. The ability of D-type cyclins to limit MMR is CDK4- and CDK6-independent and is conserved in G0 and G1. At the G1/S transition, the timely, cullin-RING ubiquitin ligase (CRL)-dependent degradation of D-type cyclins and p21 enables MMR activity to efficiently repair DNA replication errors. Persistent expression of D-type cyclins during S-phase inhibits the binding of MMR proteins to PCNA, increases the mutational burden, and promotes microsatellite instability.

Identifying Ferroptosis Inducers, HDAC, and RTK Inhibitor Sensitivity in Melanoma Subtypes through Unbiased Drug Target Prediction.

The utilization of PD1 and CTLA4 inhibitors has revolutionized the treatment of malignant melanoma (MM). However, resistance to targeted and immune-checkpoint-based therapies still poses a significant problem. Here we mine large scale MM proteogenomic data integrating it with MM cell line dependency screen, and drug sensitivity data to identify druggable targets and forecast treatment efficacy and resistance. Leveraging protein profiles from established MM subtypes and molecular structures of 82 cancer treatment drugs, we identified nine candidate hub proteins, mTOR, FYN, PIK3CB, EGFR, MAPK3, MAP4K1, MAP2K1, SRC and AKT1, across five distinct MM subtypes. These proteins serve as potential drug targets applicable to one or multiple MM subtypes. By analyzing transcriptomic data from 48 publicly accessible melanoma cell lines sourced from Achilles and CRISPR dependency screens, we forecasted 162 potentially targetable genes. We also identified genetic resistance in 260 genes across at least one melanoma subtype. In addition, we employed publicly available compound sensitivity data (Cancer Therapeutics Response Portal, CTRPv2) on the cell lines to assess the correlation of compound effectiveness within each subtype. We have identified 20 compounds exhibiting potential drug impact in at least one melanoma subtype. Remarkably, employing this unbiased approach, we have uncovered compounds targeting ferroptosis, that demonstrate a striking 30x fold difference in sensitivity among different subtypes. This implies that the proteogenomic classification of melanoma has the potential to predict sensitivity to ferroptosis compounds. Our results suggest innovative and novel therapeutic strategies by stratifying melanoma samples through proteomic profiling, offering a spectrum of novel therapeutic interventions and prospects for combination therapy. Highlights: (1) Proteogenomic subtype classification can define the landscape of genetic dependencies in melanoma (2) Nine proteins from molecular subtypes were identified as potential drug targets for specified MM patients (3) 20 compounds identified that show potential effectiveness in at least one melanoma subtype (4) Proteogenomics can predict specific ferroptosis inducers, HDAC, and RTK Inhibitor sensitivity in melanoma subtypes.

D-type cyclins regulate DNA mismatch repair in the G1 and S phases of the cell cycle, maintaining genome stability.

The large majority of oxidative DNA lesions occurring in the G1 phase of the cell cycle are repaired by base excision repair (BER) rather than mismatch repair (MMR) to avoid long resections that can lead to genomic instability and cell death. However, the molecular mechanisms dictating pathway choice between MMR and BER have remained unknown. Here, we show that, during G1, D-type cyclins are recruited to sites of oxidative DNA damage in a PCNA- and p21-dependent manner. D-type cyclins shield p21 from its two ubiquitin ligases CRL1SKP2 and CRL4CDT2 in a CDK4/6-independent manner. In turn, p21 competes through its PCNA-interacting protein degron with MMR components for their binding to PCNA. This inhibits MMR while not affecting BER. At the G1/S transition, the CRL4AMBRA1-dependent degradation of D-type cyclins renders p21 susceptible to proteolysis. These timely degradation events allow the proper binding of MMR proteins to PCNA, enabling the repair of DNA replication errors. Persistent expression of cyclin D1 during S-phase increases the mutational burden and promotes microsatellite instability. Thus, the expression of D-type cyclins inhibits MMR in G1, whereas their degradation is necessary for proper MMR function in S.

The expression profile and tumorigenic mechanisms of CD97 (ADGRE5) in glioblastoma render it a targetable vulnerability.

Glioblastoma (GBM) is the most common and aggressive primary brain malignancy. Adhesion G protein-coupled receptors (aGPCRs) have attracted interest for their potential as treatment targets. Here, we show that CD97 (ADGRE5) is the most promising aGPCR target in GBM, by virtue of its de novo expression compared to healthy brain tissue. CD97 knockdown or knockout significantly reduces the tumor initiation capacity of patient-derived GBM cultures (PDGCs) in vitro and in vivo. We find that CD97 promotes glycolytic metabolism via the mitogen-activated protein kinase (MAPK) pathway, which depends on phosphorylation of its C terminus and recruitment of ß-arrestin. We also demonstrate that THY1/CD90 is a likely CD97 ligand in GBM. Lastly, we show that an anti-CD97 antibody-drug conjugate selectively kills tumor cells in vitro. Our studies identify CD97 as a regulator of tumor metabolism, elucidate mechanisms of receptor activation and signaling, and provide strong scientific rationale for developing biologics to target it therapeutically in GBM.

Proteogenomics Reveal the Overexpression of HLA-I in Cancer.

An accurate quantification of HLA class I gene expression is important in understanding the interplay with the tumor microenvironment of antitumor cytotoxic T cell activities. Because HLA-I sequences are highly variable, standard RNAseq and mass spectrometry-based quantification workflows using common genome and protein sequence references do not provide HLA-I allele specific quantifications. Here, we used personalized HLA-I nucleotide and protein reference sequences based on the subjects' HLA-I genotypes and surveyed tumor and adjacent normal samples from patients across nine cancer types. Mass spectrometry using data dependent acquisition data was validated to be sufficient to estimate HLA-A protein expression at the allele level. We found that HLA-I proteins were present in significantly higher levels in tumors compared to adjacent normal tissues from 41 to 63% of head and neck squamous cell carcinoma, uterine corpus endometrial carcinoma, and clear cell renal cell carcinoma patients, and this was driven by increased levels of HLA-I gene transcripts. Most immune cell types are universally enriched in HLA-I high tumors, while endothelial and neuronal cells showed divergent relationships with HLA-I. Pathway analysis revealed that tumor senescence and autophagy activity influence the level of HLA-I proteins in glioblastoma. Genes correlated to HLA-I protein expression are mostly the ones directly involved in HLA-I function in immune response and cell death, while glycosylation genes are exclusively co-expressed with HLA-I at the protein level.

Engineered transcription-associated Cas9 targeting in eukaryotic cells.

DNA targeting Class 2 CRISPR-Cas effector nucleases, including the well-studied Cas9 proteins, evolved protospacer-adjacent motif (PAM) and guide RNA interactions that sequentially license their binding and cleavage activities at protospacer target sites. Both interactions are nucleic acid sequence specific but function constitutively; thus, they provide intrinsic spatial control over DNA targeting activities but naturally lack temporal control. Here we show that engineered Cas9 fusion proteins which bind to nascent RNAs near a protospacer can facilitate spatiotemporal coupling between transcription and DNA targeting at that protospacer: Transcription-associated Cas9 Targeting (TraCT). Engineered TraCT is enabled when suboptimal PAM interactions limit basal activity in vivo and when one or more nascent RNA substrates are still tethered to the actively transcribing target DNA in cis. We further show that this phenomenon can be exploited for selective editing at one of two identical targets in distinct gene loci, or, in diploid allelic loci that are differentially transcribed. Our work demonstrates that temporal control over Cas9's targeting activity at specific DNA sites may be engineered without modifying Cas9's core domains and guide RNA components or their expression levels. More broadly, it establishes RNA binding in cis as a mechanism that can conditionally stimulate CRISPR-Cas DNA targeting in eukaryotes.

MED19 encodes two unique protein isoforms that confer prostate cancer growth under low androgen through distinct gene expression programs.

MED19, a component of the mediator complex and a co-regulator of the androgen receptor (AR), is pivotal in prostate cancer cell proliferation. MED19 has two isoforms: a full-length "canonical" and a shorter "alternative" variant. Specific antibodies were developed to investigate these isoforms. Both exhibit similar expression in normal prostate development and adult prostate tissue, but the canonical isoform is elevated in prostate adenocarcinomas. Overexpression of canonical MED19 in LNCaP cells promotes growth under conditions of androgen deprivation in vitro and in vivo, mirroring earlier findings with alternative MED19-overexpressing LNCaP cells. Interestingly, alternative MED19 cells displayed strong colony formation in clonogenic assays under conditions of androgen deprivation, while canonical MED19 cells did not, suggesting distinct functional roles. These isoforms also modulated gene expression differently. Canonical MED19 triggered genes related to extracellular matrix remodeling while suppressing those involved in androgen-inactivating glucuronidation. In contrast, alternative MED19 elevated genes tied to cell movement and reduced those associated with cell adhesion and differentiation. The ratio of MED19 isoform expression in prostate cancers shifts with the disease stage. Early-stage cancers exhibit higher canonical MED19 expression than alternative MED19, consistent with canonical MED19's ability to promote cell proliferation under androgen deprivation. Conversely, alternative MED19 levels were higher in later-stage metastatic prostate cancer than in canonical MED19, reflecting alternative MED19's capability to enhance cell migration and autonomous cell growth. Our findings suggest that MED19 isoforms play unique roles in prostate cancer progression and highlights MED19 as a potential therapeutic target for both early and late-stage prostate cancer.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

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.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Integrative multi-omic cancer profiling reveals DNA methylation patterns associated with therapeutic vulnerability and cell-of-origin.

DNA methylation plays a critical role in establishing and maintaining cellular identity. However, it is frequently dysregulated during tumor development and is closely intertwined with other genetic alterations. Here, we leveraged multi-omic profiling of 687 tumors and matched non-involved adjacent tissues from the kidney, brain, pancreas, lung, head and neck, and endometrium to identify aberrant methylation associated with RNA and protein abundance changes and build a Pan-Cancer catalog. We uncovered lineage-specific epigenetic drivers including hypomethylated FGFR2 in endometrial cancer. We showed that hypermethylated STAT5A is associated with pervasive regulon downregulation and immune cell depletion, suggesting that epigenetic regulation of STAT5A expression constitutes a molecular switch for immunosuppression in squamous tumors. We further demonstrated that methylation subtype-enrichment information can explain cell-of-origin, intra-tumor heterogeneity, and tumor phenotypes. Overall, we identified cis-acting DNA methylation events that drive transcriptional and translational changes, shedding light on the tumor's epigenetic landscape and the role of its cell-of-origin.

Deep learning integrates histopathology and proteogenomics at a pan-cancer level.

We introduce a pioneering approach that integrates pathology imaging with transcriptomics and proteomics to identify predictive histology features associated with critical clinical outcomes in cancer. We utilize 2,755 H&E-stained histopathological slides from 657 patients across 6 cancer types from CPTAC. Our models effectively recapitulate distinctions readily made by human pathologists: tumor vs. normal (AUROC = 0.995) and tissue-of-origin (AUROC = 0.979). We further investigate predictive power on tasks not normally performed from H&E alone, including TP53 prediction and pathologic stage. Importantly, we describe predictive morphologies not previously utilized in a clinical setting. The incorporation of transcriptomics and proteomics identifies pathway-level signatures and cellular processes driving predictive histology features. Model generalizability and interpretability is confirmed using TCGA. We propose a classification system for these tasks, and suggest potential clinical applications for this integrated human and machine learning approach. A publicly available web-based platform implements these models.

PhosphoDisco: A Toolkit for Co-regulated Phosphorylation Module Discovery in Phosphoproteomic Data.

Kinases are key players in cancer-relevant pathways and are the targets of many successful precision cancer therapies. Phosphoproteomics is a powerful approach to study kinase activity and has been used increasingly for the characterization of tumor samples leading to the identification of novel chemotherapeutic targets and biomarkers. Finding co-regulated phosphorylation sites which represent potential kinase-substrate sets or members of the same signaling pathway allows us to harness these data to identify clinically relevant and targetable alterations in signaling cascades. Unfortunately, studies have found that databases of co-regulated phosphorylation sites are only experimentally supported in a small number of substrate sets. To address the inherent challenge of defining co-regulated phosphorylation modules relevant to a given dataset, we developed PhosphoDisco, a toolkit for determining co-regulated phosphorylation modules. We applied this approach to tandem mass spectrometry based phosphoproteomic data for breast and non-small cell lung cancer and identified canonical as well as putative new phosphorylation site modules. Our analysis identified several interesting modules in each cohort. Among these was a new cell cycle checkpoint module enriched in basal breast cancer samples and a module of PRKC isozymes putatively co-regulated by CDK12 in lung cancer. We demonstrate that modules defined by PhosphoDisco can be used to further personalized cancer treatment strategies by establishing active signaling pathways in a given patient tumor or set of tumors, and in providing new ways to classify tumors based on signaling activity.