Proteomic Characterization Identifies Clinically Relevant Subgroups of Gastrointestinal Stromal Tumors.

BACKGROUND & AIMS: Gastrointestinal stromal tumor (GIST) is the most common mesenchymal tumor of the gastrointestinal tract, and it has high metastatic and recurrence rates. We aimed to characterize the proteomic features of GIST to understand biological processes and treatment vulnerabilities. METHODS: Quantitative proteomics and phosphoproteomics analyses were performed on 193 patients with GIST to reveal the biological characteristics of GIST. Data-driven hypotheses were tested by performing functional experiments using both GIST cell lines and xenograft mouse models. RESULTS: Proteomic analysis revealed differences in the molecular features of GISTs from different locations or with different histological grades. MAPK7 was identified and functionally proved to be associated with tumor cell proliferation in GIST. Integrative analysis revealed that increased SQSTM1 expression inhibited the patient response to imatinib mesylate. Proteomics subtyping identified 4 clusters of tumors with different clinical and molecular attributes. Functional experiments confirmed the role of SRSF3 in promoting tumor cell proliferation and leading to poor prognosis. CONCLUSIONS: Our study provides a valuable data resource and highlights potential therapeutic approaches for GIST.

Blood molecular markers associated with COVID-19 immunopathology and multi-organ damage.

COVID-19 is characterized by dysregulated immune responses, metabolic dysfunction and adverse effects on the function of multiple organs. To understand host responses to COVID-19 pathophysiology, we combined transcriptomics, proteomics, and metabolomics to identify molecular markers in peripheral blood and plasma samples of 66 COVID-19-infected patients experiencing a range of disease severities and 17 healthy controls. A large number of expressed genes, proteins, metabolites, and extracellular RNAs (exRNAs) exhibit strong associations with various clinical parameters. Multiple sets of tissue-specific proteins and exRNAs varied significantly in both mild and severe patients suggesting a potential impact on tissue function. Chronic activation of neutrophils, IFN-I signaling, and a high level of inflammatory cytokines were observed in patients with severe disease progression. In contrast, COVID-19-infected patients experiencing milder disease symptoms showed robust T-cell responses. Finally, we identified genes, proteins, and exRNAs as potential biomarkers that might assist in predicting the prognosis of SARS-CoV-2 infection. These data refine our understanding of the pathophysiology and clinical progress of COVID-19.

Proteogenomic characterization of cholangiocarcinoma.

BACKGROUND AND AIMS: Cholangiocarcinoma (CCA) is a highly heterogeneous cancer with limited understanding and few effective therapeutic approaches. We aimed at providing a proteogenomic CCA characterization to inform biological processes and treatment vulnerabilities. APPROACH AND RESULTS: Integrative genomic analysis with functional validation uncovered biological perturbations downstream of driver events including DPCR1 , RBM47 mutations, SH3BGRL2 copy number alterations, and FGFR2 fusions in CCA. Proteomic clustering identified three subtypes with distinct clinical outcomes, molecular features, and potential therapeutics. Phosphoproteomics characterized targetable kinases in CCA, suggesting strategies for effective treatment with CDK and MAPK inhibitors. Patients with CCA with HBV infection showed increased antigen processing and presentation (APC) and T cell infiltration, conferring a favorable prognosis compared with those without HBV infection. The characterization of extrahepatic CCA recommended the feasible application of vascular endothelial-derived growth factor inhibitors. Multiomics profiling presented distinctive molecular characteristics of the large bile duct and the small bile duct of intrahepatic CCA. The immune landscape further revealed diverse tumor immune microenvironments, suggesting immune subtypes C1 and C5 might benefit from immune checkpoint therapy. TCN1 was identified as a potential CCA prognostic biomarker, promoting cell growth by enhancing vitamin B12 metabolism. CONCLUSIONS: We characterized the proteogenomic landscape of 217 CCAs with 197 paired normal adjacent tissues and identified their subtypes and potential therapeutic targets. The multiomics analyses with other databases and some functional validations have indicated strategies regarding the clinical, biological, and therapeutic approaches to the management of CCA.

High homocysteine is associated with idiopathic normal pressure hydrocephalus in deep perforating arteriopathy: a cross-sectional study.

BACKGROUND AND OBJECTIVE: The pathogenesis and pathophysiology of idiopathic normal pressure hydrocephalus (iNPH) remain unclear. Homocysteine may reduce the compliance of intracranial arteries and damage the endothelial function of the blood-brain barrier (BBB), which may be the underlying mechanism of iNPH. The overlap cases between deep perforating arteriopathy (DPA) and iNPH were not rare for the shared risk factors. We aimed to investigate the relationship between serum homocysteine and iNPH in DPA. METHODS: A total of 41 DPA patients with iNPH and 49 DPA patients without iNPH were included. Demographic characteristics, vascular risk factors, laboratory results, and neuroimaging data were collected. Multivariable logistic regression analysis was performed to investigate the relationship between serum homocysteine and iNPH in DPA patients. RESULTS: Patients with iNPH had significantly higher homocysteine levels than those without iNPH (median, 16.34 mmol/L versus 14.28 mmol/L; P = 0.002). There was no significant difference in CSVD burden scores between patients with iNPH and patients without iNPH. Univariate logistic regression analysis demonstrated that patients with homocysteine levels in the Tertile3 were more likely to have iNPH than those in the Tertile1 (OR, 4.929; 95% CI, 1.612-15.071; P = 0.005). The association remained significant after multivariable adjustment for potential confounders, including age, male, hypertension, diabetes mellitus, atherosclerotic cardiovascular disease (ASCVD) or hypercholesterolemia, and eGFR level. CONCLUSION: Our study indicated that high serum homocysteine levels were independently associated with iNPH in DPA. However, further research is needed to determine the predictive value of homocysteine and to confirm the underlying mechanism between homocysteine and iNPH.

Chinese newborn screening for the incidence of G6PD deficiency and variant of G6PD gene from 2013 to 2017.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is one of the most common X-linked enzymopathies caused by G6PD gene variant. We aimed to provide the characteristics of G6PD deficiency and G6PD gene variant distribution in a large Chinese newborn screening population. We investigated the prevalence of G6PD in China from 2013 to 2017. Then, we examined G6PD activity and G6PD gene in representative Chinese birth cohort to explore the distribution of G6PD gene variant in 2016. We then performed multicolor melting curve analysis to classify G6PD gene variants in 10,357 neonates with activity-confirmed G6PD deficiency, and DNA Sanger sequencing for G6PD coding exons if hot site variants were not found. The screened population, organizations, and provinces of G6PD deficiency were increased from 2013 to 2017 in China. The top five frequency of G6PD gene variants were c.1376G>T, c.1388G>A, c.95A>G, c.1024C>T, and c.871G>A and varied in different provinces, with regional and ethnic features, and four pathogenic variant sites (c.152C>T, c.290A>T, c.697G>C, and c.1285A>G) were first reported. G6PD deficiency mainly occurs in South China, and the frequency of G6PD gene variant varies in different regions and ethnicities.

rPTMDetermine: A Fully Automated Methodology for Endogenous Tyrosine Nitration Validation, Site-Localization, and Beyond.

We present herein rPTMDetermine, an adaptive and fully automated methodology for validation of the identification of rarely occurring post-translational modifications (PTMs), using a semisupervised approach with a linear discriminant analysis (LDA) algorithm. With this strategy, verification is enhanced through similarity scoring of tandem mass spectrometry (MS/MS) comparisons between modified peptides and their unmodified analogues. We applied rPTMDetermine to (1) perform fully automated validation steps for modified peptides identified from an in silico database and (2) retrieve potential yet-to-be-identified modified peptides from raw data (that had been missed through conventional database searches). In part (1), 99 of 125 3-nitrotyrosyl-containing (nitrated) peptides obtained from a ProteinPilot search were validated and localized. Twenty nitrated peptides were falsely assigned because of incorrect monoisotopic peak assignments, leading to erroneous identification of deamidation and nitration. Five additional nitrated peptides were, however, validated after performing nonmonoisotopic peak correction. In part (2), an additional 236 unique nitrated peptides were retrieved and localized, containing 113 previously unreported nitration sites; 25 endogenous nitrated peptides with novel sites were selected and verified by comparison with synthetic analogues. In summary, we identified and confidently validated 296 unique nitrated peptides-collectively representing the largest number of endogenously identified 3-nitrotyrosyl-containing peptides from the cerebral cortex proteome of a Macaca fascicularis model of stroke. Furthermore, we harnessed the rPTMDetermine strategy to complement conventional database searching and enhance the confidence of assigning rarely occurring PTMs, while recovering many missed peptides. In a final demonstration, we successfully extended the application of rPTMDetermine to peptides featuring tryptophan oxidation.

A reference peptide database for proteome quantification based on experimental mass spectrum response curves.

Motivation: Mass spectrometry (MS) based quantification of proteins/peptides has become a powerful tool in biological research with high sensitivity and throughput. The accuracy of quantification, however, has been problematic as not all peptides are suitable for quantification. Several methods and tools have been developed to identify peptides that response well in mass spectrometry and they are mainly based on predictive models, and rarely consider the linearity of the response curve, limiting the accuracy and applicability of the methods. An alternative solution is to select empirically superior peptides that offer satisfactory MS response intensity and linearity in a wide dynamic range of peptide concentration. Results: We constructed a reference database for proteome quantification based on experimental mass spectrum response curves. The intensity and dynamic range of over 2 647 773 transitions from 121 318 peptides were obtained from a set of dilution experiments, covering 11 040 gene products. These transitions and peptides were evaluated and presented in a database named SCRIPT-MAP. We showed that the best-responder (BR) peptide approach for quantification based on SCRIPT-MAP database is robust, repeatable and accurate in proteome-scale protein quantification. This study provides a reference database as well as a peptides/transitions selection method for quantitative proteomics. Availability and implementation: SCRIPT-MAP database is available at http://www.firmiana.org/responders/. Supplementary information: Supplementary data are available at Bioinformatics online.

Longitudinal plasma proteome profiling reveals the diversity of biomarkers for diagnosis and cetuximab therapy response of colorectal cancer.

Cetuximab therapy is the major treatment for colorectal cancer (CRC), but drug resistance limits its effectiveness. Here, we perform longitudinal and deep proteomic profiling of 641 plasma samples originated from 147 CRC patients (CRCs) undergoing cetuximab therapy with multi-course treatment, and 90 healthy controls (HCs). COL12A1, THBS2, S100A8, and S100A9 are screened as potential proteins to distinguish CRCs from HCs both in plasma and tissue validation cohorts. We identify the potential biomarkers (RRAS2, MMP8, FBLN1, RPTOR, and IMPDH2) for the initial response prediction. In a longitudinal setting, we identify two clusters with distinct fluctuations and construct the model with high accuracy to predict the longitudinal response, further validated in the independent cohort. This study reveals the heterogeneity of different biomarkers for tumor diagnosis, the initial and longitudinal response prediction respectively in the first course and multi-course cetuximab treatment, may ultimately be useful in monitoring and intervention strategies for CRC.

Proteomic characterization of esophageal squamous cell carcinoma response to immunotherapy reveals potential therapeutic strategy and predictive biomarkers.

Immunotherapy is the first-line therapy for esophageal squamous cell carcinoma (ESCC), yet many patients do not respond due to drug resistance and the lack of reliable predictive markers. We collected 73 ESCC patients (including discovery cohort and validation cohort) without immune thrombocytopenia and undergoing anti-PD1 immunotherapy. Proteomic and phosphoproteomic analysis of 73 ESCC treatment-naive samples by mass spectrometry-based label-free quantification were applied to explore the potential resistant and sensitive mechanisms, and identify predictive markers of ESCC immunotherapy. Comparative analysis found the pathways related to immune and mitochondrial functions were associated with ESCC immunotherapy sensitivity; while platelet activation bioprocess showed negative correlation with CD8+ T cells and related to ESCC immunotherapy non-sensitivity. Finally, we identified 10 ESCC immunotherapy predictive biomarkers with high accuracy (≥ 0.90) to predict the immunotherapeutic response, which was validated in the independent cohort.

Comprehensive Proteogenomic Profiling Reveals the Molecular Characteristics of Colorectal Cancer at Distinct Stages of Progression.

Colorectal cancer (CRC) is the second most common malignant tumor world-wide. Analysis of the changes that occur during CRC progression could provide insights into the molecular mechanisms driving CRC development and identify improved treatment strategies. Here, we performed an integrated multi-omics analysis of 435 trace-tumor-samples from 148 colorectal cancer (CRC) patients, covering non-tumor (NT), intraepithelial neoplasia (IEN), infiltration (IFT), and advanced-stage CRC (A-CRC) phases. Proteogenomics analyses demonstrated that KRAS and BRAF mutations were mutually exclusive and elevated oxidation phosphorylation in the IEN phase. Chr17q loss and chr20q gain were also mutually exclusive, occurred predominantly in the IEN and IFT phases, respectively, and impacted the cell cycle. Mutation of TP53 was frequent in the A-CRC phase and associated with tumor microenvironment, including increased extracellular matrix rigidity and stromal infiltration. Analysis of the profiles of CRC based on CMS and CRIS classifications revealed the progression paths of each subtype and indicated that microsatellite instability was associated with specific subtype classifications. Additional comparison of molecular characteristics of CRC based on location showed that ANKRD22 amplification by chr10q23.31 gain enhanced glycolysis in the right-sided CRC. The AOM/DSS-induced CRC carcinogenesis mouse model in mice indicated that DDX5 deletion due to chr17q loss promoted CRC development, consistent with the findings from the patient samples. Collectively, this study provides an informative resource for understanding the driving events of different stages of CRC and identifying the potential therapeutic targets.

Proteomic characterization of the colorectal cancer response to chemoradiation and targeted therapies reveals potential therapeutic strategies.

Chemoradiation and targeted therapies are the major treatments for colorectal cancer (CRC); however, molecular properties associated with therapy resistance are incompletely characterized. Here, we profile the proteome of 254 tumor tissues from patients with CRC undergoing chemotherapy, chemoradiation, or chemotherapy combined with targeted therapy. Proteome-based classification reveals four subtypes featured with distinct biological and therapeutic characteristics. The integrative analysis of CRC cell lines and clinical samples indicates that immune regulation is significantly associated with drug sensitivity. HSF1 can increase DNA damage repair and cell cycle, thus inducing resistance to radiation, while high expression of HDAC6 is negatively associated with response of cetuximab. Furthermore, we develop prognostic models with high accuracy to predict the therapeutic response, further validated by parallel reaction monitoring (PRM) assay in an independent validation cohort. This study provides a rich resource for investigating the mechanisms and indicators of chemoradiation and targeted therapy in CRC.

Proteogenomics of different urothelial bladder cancer stages reveals distinct molecular features for papillary cancer and carcinoma in situ.

The progression of urothelial bladder cancer (UC) is a complicated multi-step process. We perform a comprehensive multi-omics analysis of 448 samples from 190 UC patients, covering the whole spectrum of disease stages and grades. Proteogenomic integration analysis indicates the mutations of HRAS regulated mTOR signaling to form urothelial papilloma rather than papillary urothelial cancer (PUC). DNA damage is a key signaling pathway in the progression of carcinoma in situ (CIS) and related to APOBEC signature. Glucolipid metabolism increase and lower immune cell infiltration are associated with PUC compared to CIS. Proteomic analysis distinguishes the origins of invasive tumors (PUC-derived and CIS-derived), related to distinct clinical prognosis and molecular features. Additionally, loss of RBPMS, associated with CIS-derived tumors, is validated to increase the activity of AP-1 and promote metastasis. This study reveals the characteristics of two distinct branches (PUC and CIS) of UC progression and may eventually benefit clinical practice.

Comprehensive proteogenomic characterization of early duodenal cancer reveals the carcinogenesis tracks of different subtypes.

The subtypes of duodenal cancer (DC) are complicated and the carcinogenesis process is not well characterized. We present comprehensive characterization of 438 samples from 156 DC patients, covering 2 major and 5 rare subtypes. Proteogenomics reveals LYN amplification at the chromosome 8q gain functioned in the transmit from intraepithelial neoplasia phase to infiltration tumor phase via MAPK signaling, and illustrates the DST mutation improves mTOR signaling in the duodenal adenocarcinoma stage. Proteome-based analysis elucidates stage-specific molecular characterizations and carcinogenesis tracks, and defines the cancer-driving waves of the adenocarcinoma and Brunner's gland subtypes. The drug-targetable alanyl-tRNA synthetase (AARS1) in the high tumor mutation burden/immune infiltration is significantly enhanced in DC progression, and catalyzes the lysine-alanylation of poly-ADP-ribose polymerases (PARP1), which decreases the apoptosis of cancer cells, eventually promoting cell proliferation and tumorigenesis. We assess the proteogenomic landscape of early DC, and provide insights into the molecular features corresponding therapeutic targets.

Integrative proteogenomic characterization of early esophageal cancer.

Esophageal squamous cell carcinoma (ESCC) is malignant while the carcinogenesis is still unclear. Here, we perform a comprehensive multi-omics analysis of 786 trace-tumor-samples from 154 ESCC patients, covering 9 histopathological stages and 3 phases. Proteogenomics elucidates cancer-driving waves in ESCC progression, and reveals the molecular characterization of alcohol drinking habit associated signatures. We discover chromosome 3q gain functions in the transmit from nontumor to intraepithelial neoplasia phases, and find TP53 mutation enhances DNA replication in intraepithelial neoplasia phase. The mutations of AKAP9 and MCAF1 upregulate glycolysis and Wnt signaling, respectively, in advanced-stage ESCC phase. Six major tracks related to different clinical features during ESCC progression are identified, which is validated by an independent cohort with another 256 samples. Hyperphosphorylated phosphoglycerate kinase 1 (PGK1, S203) is considered as a drug target in ESCC progression. This study provides insight into the understanding of ESCC molecular mechanism and the development of therapeutic targets.

Proteogenomics of diffuse gliomas reveal molecular subtypes associated with specific therapeutic targets and immune-evasion mechanisms.

Diffuse gliomas are devastating brain tumors. Here, we perform a proteogenomic profiling of 213 retrospectively collected glioma tumors. Proteogenomic analysis reveals the downstream biological events leading by EGFR-, IDH1-, TP53-mutations. The comparative analysis illustrates the distinctive features of GBMs and LGGs, indicating CDK2 inhibitor might serve as a promising drug target for GBMs. Further proteogenomic integrative analysis combined with functional experiments highlight the cis-effect of EGFR alterations might lead to glioma tumor cell proliferation through ERK5 medicates nucleotide synthesis process. Proteome-based stratification of gliomas defines 3 proteomic subgroups (S-Ne, S-Pf, S-Im), which could serve as a complement to WHO subtypes, and would provide the essential framework for the utilization of specific targeted therapies for particular glioma subtypes. Immune clustering identifies three immune subtypes with distinctive immune cell types. Further analysis reveals higher EGFR alteration frequencies accounts for elevation of immune check point protein: PD-L1 and CD70 in T-cell infiltrated tumors.

Plasma proteomic profiling discovers molecular features associated with upper tract urothelial carcinoma.

Upper tract urothelial carcinoma (UTUC) is often diagnosed late and exhibits poor prognosis. Limited data are available on potential non-invasive biomarkers for disease monitoring. Here, we investigate the proteomic profile of plasma in 362 UTUC patients and 239 healthy controls. We present an integrated tissue-plasma proteomic approach to infer the signature proteins for identifying patients with muscle-invasive UTUC. We discover a protein panel that reflects lymph node metastasis, which is of interest in identifying UTUC patients with high risk and poor prognosis. We also identify a ten-protein classifier and establish a progression clock predicting progression-free survival of UTUC patients. Finally, we further validate the signature proteins by parallel reaction monitoring assay in an independent cohort. Collectively, this study portrays the plasma proteomic landscape of a UTUC cohort and provides a valuable resource for further biological and diagnostic research in UTUC.

Proteogenomics of clear cell renal cell carcinoma response to tyrosine kinase inhibitor.

The tyrosine kinase inhibitor (TKI) Sunitinib is one the therapies approved for advanced renal cell carcinoma. Here, we undertake proteogenomic profiling of 115 tumors from patients with clear cell renal cell carcinoma (ccRCC) undergoing Sunitinib treatment and reveal the molecular basis of differential clinical outcomes with TKI therapy. We find that chromosome 7q gain-induced mTOR signaling activation is associated with poor therapeutic outcomes with Sunitinib treatment, whereas the aristolochic acid signature and VHL mutation synergistically caused enhanced glycolysis is correlated with better prognosis. The proteomic and phosphoproteomic analysis further highlights the responsibility of mTOR signaling for non-response to Sunitinib. Immune landscape characterization reveals diverse tumor microenvironment subsets in ccRCC. Finally, we construct a multi-omics classifier that can detect responder and non-responder patients (receiver operating characteristic-area under the curve, 0.98). Our study highlights associations between ccRCC molecular characteristics and the response to TKI, which can facilitate future improvement of therapeutic responses.

Differential Paralog-Specific Expression of Multiple Small Subunit Proteins Cause Variations in Rpl42/eL42 Incorporation in Ribosome in Fission Yeast.

Ribosomes within a cell are commonly viewed as biochemically homogenous RNA-protein super-complexes performing identical functions of protein synthesis. However, recent evidence suggests that ribosomes may be a more dynamic macromolecular complex with specialized roles. Here, we present extensive genetic and molecular evidence in the fission yeast S. pombe that the paralogous genes for many ribosomal proteins (RPs) are functionally different, despite that they encode the same ribosomal component, often with only subtle differences in the sequences. Focusing on the rps8 paralog gene deletions rps801d and rps802d, we showed that the mutant cells differ in the level of Rpl42p in actively translating ribosomes and that their phenotypic differences reside in the Rpl42p level variation instead of the subtle protein sequence difference between Rps801p and Rps802p. Additional 40S ribosomal protein paralog pairs also exhibit similar phenotypic differences via differential Rpl42p levels in actively translating ribosomes. Together, our work identifies variations in the Rpl42p level as a potential form of ribosome heterogeneity in biochemical compositions and suggests a possible connection between large and small subunits during ribosome biogenesis that may cause such heterogeneity. Additionally, it illustrates the complexity of the underlying mechanisms for the genetic specificity of ribosome paralogs.

Integrated proteogenomic characterization of urothelial carcinoma of the bladder.

BACKGROUND: Urothelial carcinoma (UC) is the most common pathological type of bladder cancer, a malignant tumor. However, an integrated multi-omics analysis of the Chinese UC patient cohort is lacking. METHODS: We performed an integrated multi-omics analysis, including whole-exome sequencing, RNA-seq, proteomic, and phosphoproteomic analysis of 116 Chinese UC patients, comprising 45 non-muscle-invasive bladder cancer patients (NMIBCs) and 71 muscle-invasive bladder cancer patients (MIBCs). RESULT: Proteogenomic integration analysis indicated that SND1 and CDK5 amplifications on chromosome 7q were associated with the activation of STAT3, which was relevant to tumor proliferation. Chromosome 5p gain in NMIBC patients was a high-risk factor, through modulating actin cytoskeleton implicating in tumor cells invasion. Phosphoproteomic analysis of tumors and morphologically normal human urothelium produced UC-associated activated kinases, including CDK1 and PRKDC. Proteomic analysis identified three groups, U-I, U-II, and U-III, reflecting distinct clinical prognosis and molecular signatures. Immune subtypes of UC tumors revealed a complex immune landscape and suggested the amplification of TRAF2 related to the increased expression of PD-L1. Additionally, increased GARS, related to subtype U-II, was validated to promote pentose phosphate pathway by inhibiting activities of PGK1 and PKM2. CONCLUSIONS: This study provides a valuable resource for researchers and clinicians to further identify molecular pathogenesis and therapeutic opportunities in urothelial carcinoma of the bladder.