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.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

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.

Proteogenomic characterization of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Integrated Proteogenomic Characterization of HBV-Related Hepatocellular Carcinoma.

We performed the first proteogenomic characterization of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) using paired tumor and adjacent liver tissues from 159 patients. Integrated proteogenomic analyses revealed consistency and discordance among multi-omics, activation status of key signaling pathways, and liver-specific metabolic reprogramming in HBV-related HCC. Proteomic profiling identified three subgroups associated with clinical and molecular attributes including patient survival, tumor thrombus, genetic profile, and the liver-specific proteome. These proteomic subgroups have distinct features in metabolic reprogramming, microenvironment dysregulation, cell proliferation, and potential therapeutics. Two prognostic biomarkers, PYCR2 and ADH1A, related to proteomic subgrouping and involved in HCC metabolic reprogramming, were identified. CTNNB1 and TP53 mutation-associated signaling and metabolic profiles were revealed, among which mutated CTNNB1-associated ALDOA phosphorylation was validated to promote glycolysis and cell proliferation. Our study provides a valuable resource that significantly expands the knowledge of HBV-related HCC and may eventually benefit clinical practice.

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.

Epigenetic regulation during cancer transitions across 11 tumour types.

Chromatin accessibility is essential in regulating gene expression and cellular identity, and alterations in accessibility have been implicated in driving cancer initiation, progression and metastasis1-4. Although the genetic contributions to oncogenic transitions have been investigated, epigenetic drivers remain less understood. Here we constructed a pan-cancer epigenetic and transcriptomic atlas using single-nucleus chromatin accessibility data (using single-nucleus assay for transposase-accessible chromatin) from 225 samples and matched single-cell or single-nucleus RNA-sequencing expression data from 206 samples. With over 1 million cells from each platform analysed through the enrichment of accessible chromatin regions, transcription factor motifs and regulons, we identified epigenetic drivers associated with cancer transitions. Some epigenetic drivers appeared in multiple cancers (for example, regulatory regions of ABCC1 and VEGFA; GATA6 and FOX-family motifs), whereas others were cancer specific (for example, regulatory regions of FGF19, ASAP2 and EN1, and the PBX3 motif). Among epigenetically altered pathways, TP53, hypoxia and TNF signalling were linked to cancer initiation, whereas oestrogen response, epithelial-mesenchymal transition and apical junction were tied to metastatic transition. Furthermore, we revealed a marked correlation between enhancer accessibility and gene expression and uncovered cooperation between epigenetic and genetic drivers. This atlas provides a foundation for further investigation of epigenetic dynamics in cancer transitions.

SARS-CoV-2 disrupts host epigenetic regulation via histone mimicry.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and caused the devastating global pandemic of coronavirus disease 2019 (COVID-19), in part because of its ability to effectively suppress host cell responses1-3. In rare cases, viral proteins dampen antiviral responses by mimicking critical regions of human histone proteins4-8, particularly those containing post-translational modifications required for transcriptional regulation9-11. Recent work has demonstrated that SARS-CoV-2 markedly disrupts host cell epigenetic regulation12-14. However, how SARS-CoV-2 controls the host cell epigenome and whether it uses histone mimicry to do so remain unclear. Here we show that the SARS-CoV-2 protein encoded by ORF8 (ORF8) functions as a histone mimic of the ARKS motifs in histone H3 to disrupt host cell epigenetic regulation. ORF8 is associated with chromatin, disrupts regulation of critical histone post-translational modifications and promotes chromatin compaction. Deletion of either the ORF8 gene or the histone mimic site attenuates the ability of SARS-CoV-2 to disrupt host cell chromatin, affects the transcriptional response to infection and attenuates viral genome copy number. These findings demonstrate a new function of ORF8 and a mechanism through which SARS-CoV-2 disrupts host cell epigenetic regulation. Further, this work provides a molecular basis for the finding that SARS-CoV-2 lacking ORF8 is associated with decreased severity of COVID-19.

Novel mutations in FLVCR1 cause tremors, sensory neuropathy with retinitis pigmentosa.

The mutations of the feline leukemia virus subgroup C receptor-related protein 1 (FLVCR1) cause ataxia with retinitis pigmentosa. Recent studies indicated a large variation in the phenotype of FLVCR1-associated diseases. In this report, we describe an adult male who manifested first with tremors in his third decade, followed by retinitis pigmentosa, sensory ataxia, and sensory neuropathy in his fourth decade. While retinitis pigmentosa and sensory ataxia are well-recognized features of FLVCR1-associated disease, tremor is rarely described. Whole-exome sequencing revealed novel compound heterozygous pathogenic FLVCR1 variants: c.498 G > A; p.(Trp166*) and c.369 T > G; p.(Phe123Leu). In addition, we have highlighted the ultrastructural abnormalities of the sural biopsy in this patient.

Adenovirus prevents dsRNA formation by promoting efficient splicing of viral RNA.

Eukaryotic cells recognize intracellular pathogens through pattern recognition receptors, including sensors of aberrant nucleic acid structures. Sensors of double-stranded RNA (dsRNA) are known to detect replication intermediates of RNA viruses. It has long been suggested that annealing of mRNA from symmetrical transcription of both top and bottom strands of DNA virus genomes can produce dsRNA during infection. Supporting this hypothesis, nearly all DNA viruses encode inhibitors of dsRNA-recognition pathways. However, direct evidence that DNA viruses produce dsRNA is lacking. Contrary to dogma, we show that the nuclear-replicating DNA virus adenovirus (AdV) does not produce detectable levels of dsRNA during infection. In contrast, abundant dsRNA is detected within the nucleus of cells infected with AdV mutants defective for viral RNA processing. In the presence of nuclear dsRNA, the cytoplasmic dsRNA sensor PKR is relocalized and activated within the nucleus. Accumulation of viral dsRNA occurs in the late phase of infection, when unspliced viral transcripts form intron/exon base pairs between top and bottom strand transcripts. We propose that DNA viruses actively limit dsRNA formation by promoting efficient splicing and mRNA processing, thus avoiding detection and restriction by host innate immune sensors of pathogenic nucleic acids.

Zika virus employs the host antiviral RNase L protein to support replication factory assembly.

Infection with the flavivirus Zika virus (ZIKV) can result in tissue tropism, disease outcome, and route of transmission distinct from those of other flaviviruses; therefore, we aimed to identify host machinery that exclusively promotes the ZIKV replication cycle, which can inform on differences at the organismal level. We previously reported that deletion of the host antiviral ribonuclease L (RNase L) protein decreases ZIKV production. Canonical RNase L catalytic activity typically restricts viral infection, including that of the flavivirus dengue virus (DENV), suggesting an unconventional, proviral RNase L function during ZIKV infection. In this study, we reveal that an inactive form of RNase L supports assembly of ZIKV replication factories (RFs) to enhance infectious virus production. Compared with the densely concentrated ZIKV RFs generated with RNase L present, deletion of RNase L induced broader subcellular distribution of ZIKV replication intermediate double-stranded RNA (dsRNA) and NS3 protease, two constituents of ZIKV RFs. An inactive form of RNase L was sufficient to contain ZIKV genome and dsRNA within a smaller RF area, which subsequently increased infectious ZIKV release from the cell. Inactive RNase L can interact with cytoskeleton, and flaviviruses remodel cytoskeleton to construct RFs. Thus, we used the microtubule-stabilization drug paclitaxel to demonstrate that ZIKV repurposes RNase L to facilitate the cytoskeleton rearrangements required for proper generation of RFs. During infection with flaviviruses DENV or West Nile Kunjin virus, inactive RNase L did not improve virus production, suggesting that a proviral RNase L role is not a general feature of all flavivirus infections.

SARS-CoV-2 induces double-stranded RNA-mediated innate immune responses in respiratory epithelial-derived cells and cardiomyocytes.

Coronaviruses are adept at evading host antiviral pathways induced by viral double-stranded RNA, including interferon (IFN) signaling, oligoadenylate synthetase-ribonuclease L (OAS-RNase L), and protein kinase R (PKR). While dysregulated or inadequate IFN responses have been associated with severe coronavirus infection, the extent to which the recently emerged SARS-CoV-2 activates or antagonizes these pathways is relatively unknown. We found that SARS-CoV-2 infects patient-derived nasal epithelial cells, present at the initial site of infection; induced pluripotent stem cell-derived alveolar type 2 cells (iAT2), the major cell type infected in the lung; and cardiomyocytes (iCM), consistent with cardiovascular consequences of COVID-19 disease. Robust activation of IFN or OAS-RNase L is not observed in these cell types, whereas PKR activation is evident in iAT2 and iCM. In SARS-CoV-2-infected Calu-3 and A549ACE2 lung-derived cell lines, IFN induction remains relatively weak; however, activation of OAS-RNase L and PKR is observed. This is in contrast to Middle East respiratory syndrome (MERS)-CoV, which effectively inhibits IFN signaling and OAS-RNase L and PKR pathways, but is similar to mutant MERS-CoV lacking innate immune antagonists. Remarkably, OAS-RNase L and PKR are activated in MAVS knockout A549ACE2 cells, demonstrating that SARS-CoV-2 can induce these host antiviral pathways despite minimal IFN production. Moreover, increased replication and cytopathic effect in RNASEL knockout A549ACE2 cells implicates OAS-RNase L in restricting SARS-CoV-2. Finally, while SARS-CoV-2 fails to antagonize these host defense pathways, which contrasts with other coronaviruses, the IFN signaling response is generally weak. These host-virus interactions may contribute to the unique pathogenesis of SARS-CoV-2.

MicroRNA-7a-5p ameliorates diabetic peripheral neuropathy by regulating VDAC1/JNK/c-JUN pathway.

AIMS: The pathogenesis of diabetic peripheral neuropathy (DPN) is complex, and its treatment is extremely challenging. MicroRNA-7a-5p (miR-7a-5p) has been widely reported to alleviate apoptosis and oxidative stress in various diseases. This study aimed to investigate the mechanism of miR-7a-5p in DPN. METHODS: DPN cell model was constructed with high-glucose-induced RSC96 cells. Cell apoptosis and viability were detected by flow cytometry analysis and cell counting kit-8 (CCK-8) assay respectively. The apoptosis and Jun N-terminal kinase (JNK)/c-JUN signalling pathway-related proteins expression were detected by Western blotting. The intracellular calcium content and oxidative stress levels were detected by flow cytometry and reagent kits. Mitochondrial membrane potential was evaluated by tetrechloro-tetraethylbenzimidazol carbocyanine iodide (JC-1) staining. The targeting relationship between miR-7a-5p and voltage-dependent anion-selective channel protein 1 (VDAC1) was determined by RNA pull-down assay and dual-luciferase reporter gene assay. The streptozotocin (STZ) rat model was constructed to simulate DPN in vivo. The paw withdrawal mechanical threshold (PTW) was measured by Frey capillary line, and the motor nerve conduction velocity (MNCV) was measured by electromyography. RESULTS: MiR-7a-5p expression was decreased, while VDAC1 expression was increased in HG-induced RSC96 cells and STZ rats. In HG-induced RSC96 cells, miR-7a-5p overexpression promoted cell proliferation, inhibited apoptosis, down-regulated calcium release, improved mitochondrial membrane potential and repressed oxidative stress response. MiR-7a-5p negatively regulated VDAC1 expression. VDAC1 knockdown improved cell proliferation activity, suppressed cell apoptosis and mitochondrial dysfunction by inhibiting JNK/c-JUN pathway activation. MiR-7a-5p overexpression raised PTW, restored MNCV and reduced oxidative stress levels and nerve cell apoptosis in STZ rats. CONCLUSION: MiR-7a-5p overexpression ameliorated mitochondrial dysfunction and inhibited apoptosis in DPN by regulating VDAC1/JNK/c-JUN pathway.

Lysine demethylase KDM5B down-regulates SIRT3-mediated mitochondrial glucose and lipid metabolism in diabetic neuropathy.

BACKGROUND: Diabetic peripheral neuropathy (DPN) is a common neurological complication of diabetes mellitus without efficient interventions. Both lysine demethylase 5B (KDM5B) and sirtuin-3 (SIRT3) have been found to regulate islet function and glucose homeostasis. KDM5B was predicted to bind to the SIRT3 promoter by bioinformatics. Here, we investigated whether KDM5B affected DPN development via modulating SIRT3. METHODS: The db/db mice and high glucose-stimulated Schwann cells (RSC96) were used as in vivo and in vitro models of DPN, respectively. Glucose level, glucose and insulin tolerance of mice were measured. Neurological function was evaluated by motor nerve conduction velocity (MNCV), tactile allodynia assay and thermal sensitivity assay. Adenosine triphosphate level, oxygen consumption rate, extracellular acidification rate, ß-oxidation rate, acetyl-CoA level, acetylation levels and activities of long-chain acyl CoA dehydrogenase (LCAD) and pyruvate dehydrogenase (PDH) were detected. Methyl thiazolyl tetrazolium assay was adopted to determine cell viability. Reactive oxygen species (ROS) production was detected by MitoSox staining. Western blotting for measuring target protein levels. Molecular mechanisms were investigated by co-immunoprecipitine (Co-IP), chromatin immunoprecipitation (ChIP) and luciferase reporter assay. RESULTS: KDM5B was up-regulated, while SIRT3 was down-regulated in DPN models. SIRT3 overexpression or AMPK activation ameliorated mitochondrial metabolism dysfunction and ROS overproduction during DPN. KDM5B overexpression triggered mitochondrial metabolism disorder and oxidative stress via directly transcriptional inhibiting SIRT3 expression by demethylating H3K4me3 or indirectly repressing AMPK pathway-regulated SIRT3 expression. CONCLUSION: KDM5B contributes to DPN via regulating SIRT3-mediated mitochondrial glucose and lipid metabolism. KDM5B inhibition may be an effective intervention for DPN.