Epidemiologic Trends In Children With Toxicologic Exposures Requiring Intensive Care Before and During the COVID-19 Pandemic.

OBJECTIVE: Pediatric poison exposures are a common reason for pediatric intensive care unit (PICU) -admission. The purpose of this study was to examine the exposure trends and patient outcomes in 2018-2019 compared with 2020-2021 amidst the COVID-19 pandemic. METHODS: This was a retrospective cohort study of patients 18 years of age or younger with a suspected toxicologic exposure from January 2018 to March 2021. The primary endpoint was rate of PICU admissions between the 2 cohorts. Secondary endpoints included medical outcome stratified by severity, PICU length of stay, and need for mechanical ventilation. RESULTS: Our study included a total of 340 patients with median age 14.5 (IQR, 11.9-16.1) years. There was no significant difference in age, sex, or race between the 2 cohorts. The percentage of patients admitted to the PICU for poison exposures was significantly higher in the COVID-19 cohort compared with the pre-COVID-19 cohort (8.4% vs 3.7%, p < 0.01). Severity of medical outcomes differed between the groups; the COVID-19 cohort had more extreme clinical presentations of no effect or death (p < 0.01). No significant difference was found among the remaining secondary outcomes. Classes of substances ingested were comparable with baseline poison center data. CONCLUSIONS: Poisoning-related PICU admissions occurred at more than twice the pre-pandemic rate. This may emphasize the effect of the COVID-19 pandemic on pediatric access and exposure to poisons.

Proteomic analysis of the urothelial cancer landscape.

Urothelial bladder cancer (UC) has a wide tumor biological spectrum with challenging prognostic stratification and relevant therapy-associated morbidity. Most molecular classifications relate only indirectly to the therapeutically relevant protein level. We improve the pre-analytics of clinical samples for proteome analyses and characterize a cohort of 434 samples with 242 tumors and 192 paired normal mucosae covering the full range of UC. We evaluate sample-wise tumor specificity and rank biomarkers by target relevance. We identify robust proteomic subtypes with prognostic information independent from histopathological groups. In silico drug prediction suggests efficacy of several compounds hitherto not in clinical use. Both in silico and in vitro data indicate predictive value of the proteomic clusters for these drugs. We underline that proteomics is relevant for personalized oncology and provide abundance and tumor specificity data for a large part of the UC proteome ( www.cancerproteins.org ).

Cystatin C loaded in brain-derived extracellular vesicles rescues synapses after ischemic insult in vitro and in vivo.

Synaptic loss is an early event in the penumbra area after an ischemic stroke. Promoting synaptic preservation in this area would likely improve functional neurological recovery. We aimed to detect proteins involved in endogenous protection mechanisms of synapses in the penumbra after stroke and to analyse potential beneficial effects of these candidates for a prospective stroke treatment. For this, we performed Liquid Chromatography coupled to Mass Spectrometry (LC-MS)-based proteomics of synaptosomes isolated from the ipsilateral hemispheres of mice subjected to experimental stroke at different time points (24 h, 4 and 7 days) and compared them to sham-operated mice. Proteomic analyses indicated that, among the differentially expressed proteins between the two groups, cystatin C (CysC) was significantly increased at 24 h and 4 days following stroke, before returning to steady-state levels at 7 days, thus indicating a potential transient and intrinsic rescue mechanism attempt of neurons. When CysC was applied to primary neuronal cultures subjected to an in vitro model of ischemic damage, this treatment significantly improved the preservation of synaptic structures. Notably, similar effects were observed when CysC was loaded into brain-derived extracellular vesicles (BDEVs). Finally, when CysC contained in BDEVs was administered intracerebroventricularly to stroked mice, it significantly increased the expression of synaptic markers such as SNAP25, Homer-1, and NCAM in the penumbra area compared to the group supplied with empty BDEVs. Thus, we show that CysC-loaded BDEVs promote synaptic protection after ischemic damage in vitro and in vivo, opening the possibility of a therapeutic use in stroke patients.

HERC5 downregulation in non-small cell lung cancer is associated with altered energy metabolism and metastasis.

BACKGROUND: Metastasis is the leading cause of cancer-related death in non-small cell lung cancer (NSCLC) patients. We previously showed that low HERC5 expression predicts early tumor dissemination and a dismal prognosis in NSCLC patients. Here, we performed functional studies to unravel the mechanism underlying the "metastasis-suppressor" effect of HERC5, with a focus on mitochondrial metabolism pathways. METHODS: We assessed cell proliferation, colony formation potential, anchorage-independent growth, migration, and wound healing in NSCLC cell line models with HERC5 overexpression (OE) or knockout (KO). To study early tumor cell dissemination, we used these cell line models in zebrafish experiments and performed intracardial injections in nude mice. Mass spectrometry (MS) was used to analyze protein changes in whole-cell extracts. Furthermore, electron microscopy (EM) imaging, cellular respiration, glycolytic activity, and lactate production were used to investigate the relationships with mitochondrial energy metabolism pathways. RESULTS: Using different in vitro NSCLC cell line models, we showed that NSCLC cells with low HERC5 expression had increased malignant and invasive properties. Furthermore, two different in vivo models in zebrafish and a xenograft mouse model showed increased dissemination and metastasis formation (in particular in the brain). Functional enrichment clustering of MS data revealed an increase in mitochondrial proteins in vitro when HERC5 levels were high. Loss of HERC5 leads to an increased Warburg effect, leading to improved adaptation and survival under prolonged inhibition of oxidative phosphorylation. CONCLUSIONS: Taken together, these results indicate that low HERC5 expression increases the metastatic potential of NSCLC in vitro and in vivo. Furthermore, HERC5-induced proteomic changes influence mitochondrial pathways, ultimately leading to alterations in energy metabolism and demonstrating its role as a new potential metastasis suppressor gene.

Temporal change of DNA methylation subclasses between matched newly diagnosed and recurrent glioblastoma.

The longitudinal transition of phenotypes is pivotal in glioblastoma treatment resistance and DNA methylation emerged as an important tool for classifying glioblastoma phenotypes. We aimed to characterize DNA methylation subclass heterogeneity during progression and assess its clinical impact. Matched tissues from 47 glioblastoma patients were subjected to DNA methylation profiling, including CpG-site alterations, tissue and serum deconvolution, mass spectrometry, and immunoassay. Effects of clinical characteristics on temporal changes and outcomes were studied. Among 47 patients, 8 (17.0%) had non-matching classifications at recurrence. In the remaining 39 cases, 28.2% showed dominant DNA methylation subclass transitions, with 72.7% being a mesenchymal subclass. In general, glioblastomas with a subclass transition showed upregulated metabolic processes. Newly diagnosed glioblastomas with mesenchymal transition displayed increased stem cell-like states and decreased immune components at diagnosis and exhibited elevated immune signatures and cytokine levels in serum. In contrast, tissue of recurrent glioblastomas with mesenchymal transition showed increased immune components but decreased stem cell-like states. Survival analyses revealed comparable outcomes for patients with and without subclass transitions. This study demonstrates a temporal heterogeneity of DNA methylation subclasses in 28.2% of glioblastomas, not impacting patient survival. Changes in cell state composition associated with subclass transition may be crucial for recurrent glioblastoma targeted therapies.

Fatty acid synthesis suppresses dietary polyunsaturated fatty acid use.

Dietary polyunsaturated fatty acids (PUFA) are increasingly recognized for their health benefits, whereas a high production of endogenous fatty acids - a process called de novo lipogenesis (DNL) - is closely linked to metabolic diseases. Determinants of PUFA incorporation into complex lipids are insufficiently understood and may influence the onset and progression of metabolic diseases. Here we show that fatty acid synthase (FASN), the key enzyme of DNL, critically determines the use of dietary PUFA in mice and humans. Moreover, the combination of FASN inhibition and PUFA-supplementation decreases liver triacylglycerols (TAG) in mice fed with high-fat diet. Mechanistically, FASN inhibition causes higher PUFA uptake via the lysophosphatidylcholine transporter MFSD2A, and a diacylglycerol O-acyltransferase 2 (DGAT2)-dependent incorporation of PUFA into TAG. Overall, the outcome of PUFA supplementation may depend on the degree of endogenous DNL and combining PUFA supplementation and FASN inhibition might be a promising approach to target metabolic disease.

Application of sample displacement batch chromatography for fractionation of proteoforms.

Fractionation of proteoforms is currently the most challenging topic in the field of proteoform analysis. The need for considering the existence of proteoforms in experimental approaches is not only important in Life Science research in general but especially in the manufacturing of therapeutic proteins (TPs) like recombinant therapeutic antibodies (mAbs). Some of the proteoforms of TPs have significantly decreased actions or even cause side effects. The identification and removal of proteoforms differing from the main species, having the desired action, is challenging because the difference in the composition of atoms is often very small and their concentration in comparison to the main proteoform can be low. In this study, we demonstrate that sample displacement batch chromatography (SDBC) is an easy-to-handle, economical, and efficient method for fractionating proteoforms. As a model sample a commercial ovalbumin fraction was used, containing many ovalbumin proteoforms. The most promising parameters for the SDBC were determined by a screening approach and applied for a 10-segment fractionation of ovalbumin with cation exchange chromatography resins. Mass spectrometry of intact proteoforms was used for characterizing the SDBC fractionation process. By SDBC, a significant separation of different proteoforms was obtained.

Protocol for differential multi-omic analyses of distinct cell types in the mouse cerebral cortex.

Here, we present a protocol for differential multi-omic analyses of distinct cell types in the developing mouse cerebral cortex. We describe steps for in utero electroporation, subsequent flow-cytometry-based isolation of developing mouse cortical cells, bulk RNA sequencing or quantitative liquid chromatography-tandem mass spectrometry, and bioinformatic analyses. This protocol can be applied to compare the proteomes and transcriptomes of developing mouse cortical cell populations after various manipulations (e.g., epigenetic). For complete details on the use and execution of this protocol, please refer to Meka et al. (2022).1.

Integrated proteomics spotlight the proteasome as a therapeutic vulnerability in embryonal tumors with multilayered rosettes.

BACKGROUND: Embryonal tumors with multilayered rosettes (ETMR) are rare malignant embryonal brain tumors. The prognosis of ETMR is poor and novel therapeutic approaches are desperately needed. Comprehension of ETMR tumor biology is currently based on only few previous molecular studies, which mainly focused on the analyses of nucleic acids. In this study, we explored integrated ETMR proteomics. METHODS: Using mass spectrometry, proteome data were acquired from 16 ETMR and the ETMR cell line BT183. Proteome data were integrated with case-matched global DNA methylation data, publicly available transcriptome data, and proteome data of further embryonal and pediatric brain tumors. RESULTS: Proteome-based cluster analyses grouped ETMR samples according to histomorphology, separating neuropil-rich tumors with neuronal signatures from primitive tumors with signatures relating to stemness and chromosome organization. Integrated proteomics showcased that ETMR and BT183 cells harbor proteasome regulatory proteins in abundance, implicating their strong dependency on the proteasome machinery to safeguard proteostasis. Indeed, in vitro assays using BT183 highlighted that ETMR tumor cells are highly vulnerable toward treatment with the CNS penetrant proteasome inhibitor Marizomib. CONCLUSIONS: In summary, histomorphology stipulates the proteome signatures of ETMR, and proteasome regulatory proteins are pervasively abundant in these tumors. As validated in vitro, proteasome inhibition poses a promising therapeutic option in ETMR.

Direct 3D Sampling of the Embryonic Mouse Head: Layer-wise Nanosecond Infrared Laser (NIRL) Ablation from Scalp to Cortex for Spatially Resolved Proteomics.

Common workflows in bottom-up proteomics require homogenization of tissue samples to gain access to the biomolecules within the cells. The homogenized tissue samples often contain many different cell types, thereby representing an average of the natural proteome composition, and rare cell types are not sufficiently represented. To overcome this problem, small-volume sampling and spatial resolution are needed to maintain a better representation of the sample composition and their proteome signatures. Using nanosecond infrared laser ablation, the region of interest can be targeted in a three-dimensional (3D) fashion, whereby the spatial information is maintained during the simultaneous process of sampling and homogenization. In this study, we ablated 40 µm thick consecutive layers directly from the scalp through the cortex of embryonic mouse heads and analyzed them by subsequent bottom-up proteomics. Extra- and intracranial ablated layers showed distinct proteome profiles comprising expected cell-specific proteins. Additionally, known cortex markers like SOX2, KI67, NESTIN, and MAP2 showed a layer-specific spatial protein abundance distribution. We propose potential new marker proteins for cortex layers, such as MTA1 and NMRAL1. The obtained data confirm that the new 3D tissue sampling and homogenization method is well suited for investigating the spatial proteome signature of tissue samples in a layerwise manner. Characterization of the proteome composition of embryonic skin and bone structures, meninges, and cortex lamination in situ enables a better understanding of molecular mechanisms of development during embryogenesis and disease pathogenesis.

PKA mediates modality-specific modulation of the mechanically gated ion channel PIEZO2.

PKA is a downstream effector of many inflammatory mediators that induce pain hypersensitivity by increasing the mechanosensitivity of nociceptive sensory afferent. Here, we examine the molecular mechanism underlying PKA-dependent modulation of the mechanically activated ion channel PIEZO2, which confers mechanosensitivity to many nociceptors. Using phosphorylation site prediction algorithms, we identified multiple putative and highly conserved PKA phosphorylation sites located on intracellular intrinsically disordered regions of PIEZO2. Site-directed mutagenesis and patch-clamp recordings showed that substitution of one or multiple putative PKA sites within a single intracellular domain does not alter PKA-induced PIEZO2 sensitization, whereas mutation of a combination of nine putative sites located on four different intracellular regions completely abolishes PKA-dependent PIEZO2 modulation, though it remains unclear whether all or just some of these nine sites are required. By demonstrating that PIEZO1 is not modulated by PKA, our data also reveal a previously unrecognized functional difference between PIEZO1 and PIEZO2. Moreover, by demonstrating that PKA only modulates PIEZO2 currents evoked by focal mechanical indentation of the cell, but not currents evoked by pressure-induced membrane stretch, we provide evidence suggesting that PIEZO2 is a polymodal mechanosensor that engages different protein domains for detecting different types of mechanical stimuli.

Cutaneous manifestations in Costello syndrome: HRAS p.Gly12Ser affects RIN1-mediated integrin trafficking in immortalized epidermal keratinocytes.

Heterozygous germline missense variants in the HRAS gene underlie Costello syndrome (CS). The molecular basis for cutaneous manifestations in CS is largely unknown. We used an immortalized human cell line, HaCaT keratinocytes, stably expressing wild-type or CS-associated (p.Gly12Ser) HRAS and defined RIN1 as quantitatively most prominent, high-affinity effector of active HRAS in these cells. As an exchange factor for RAB5 GTPases, RIN1 is involved in endosomal sorting of cell-adhesion integrins. RIN1-dependent RAB5A activation was strongly increased by HRASGly12Ser, and HRAS-RIN1-ABL1/2 signaling was induced in HRASWT- and HRASGly12Ser-expressing cells. Along with that, HRASGly12Ser expression decreased total integrin levels and enriched ß1 integrin in RAB5- and EEA1-positive early endosomes. The intracellular level of active ß1 integrin was increased in HRASGly12Ser HaCaT keratinocytes due to impaired recycling, whereas RIN1 disruption raised ß1 integrin cell surface distribution. HRASGly12Ser induced co-localization of ß1 integrin with SNX17 and RAB7 in early/sorting and late endosomes, respectively. Thus, by retaining ß1 integrin in intracellular endosomal compartments, HRAS-RIN1 signaling affects the subcellular availability of ß1 integrin. This may interfere with integrin-dependent processes as we detected for HRASGly12Ser cells spreading on fibronectin. We conclude that dysregulation of receptor trafficking and integrin-dependent processes such as cell adhesion are relevant in the pathobiology of CS.

Tissue Sampling and Homogenization with NIRL Enables Spatially Resolved Cell Layer Specific Proteomic Analysis of the Murine Intestine.

For investigating the molecular physiology and pathophysiology in organs, the most exact data should be obtained; if not, organ-specific cell lines are analyzed, or the whole organ is homogenized, followed by the analysis of its biomolecules. However, if the morphological organization of the organ can be addressed, then, in the best case, the composition of molecules in single cells of the target organ can be analyzed. Laser capture microdissection (LCM) is a technique which enables the selection of specific cells of a tissue for further analysis of their molecules. However, LCM is a time-consuming two-dimensional technique, and optimal results are only obtained if the tissue is fixed, e.g., by formalin. Especially for proteome analysis, formalin fixation reduced the number of identifiable proteins, and this is an additional drawback. Recently, it was demonstrated that sampling of fresh-frozen (non-fixed) tissue with an infrared-laser is giving higher yields with respect to the absolute protein amount and number of identifiable proteins than conventional mechanical homogenization of tissues. In this study, the applicability of the infrared laser tissue sampling for the proteome analysis of different cell layers of murine intestine was investigated, using LC-MS/MS-based differential quantitative bottom-up proteomics. By laser ablation, eight consecutive layers of colon tissue were obtained and analyzed. However, a clear distinguishability of protein profiles between ascending, descending, and transversal colon was made, and we identified the different intestinal-cell-layer proteins, which are cell-specific, as confirmed by data from the Human Protein Atlas. Thus, for the first time, sampling directly from intact fresh-frozen tissue with three-dimensional resolution is giving access to the different proteomes of different cell layers of colon tissue.

HarmonizR enables data harmonization across independent proteomic datasets with appropriate handling of missing values.

Dataset integration is common practice to overcome limitations in statistically underpowered omics datasets. Proteome datasets display high technical variability and frequent missing values. Sophisticated strategies for batch effect reduction are lacking or rely on error-prone data imputation. Here we introduce HarmonizR, a data harmonization tool with appropriate missing value handling. The method exploits the structure of available data and matrix dissection for minimal data loss, without data imputation. This strategy implements two common batch effect reduction methods-ComBat and limma (removeBatchEffect()). The HarmonizR strategy, evaluated on four exemplarily analyzed datasets with up to 23 batches, demonstrated successful data harmonization for different tissue preservation techniques, LC-MS/MS instrumentation setups, and quantification approaches. Compared to data imputation methods, HarmonizR was more efficient and performed superior regarding the detection of significant proteins. HarmonizR is an efficient tool for missing data tolerant experimental variance reduction and is easily adjustable for individual dataset properties and user preferences.

Circulating Cellular Communication Network Factor 1 Protein as a Sensitive Liquid Biopsy Marker for Early Detection of Breast Cancer.

BACKGROUND: Despite recent progress in liquid biopsy technologies, early blood-based detection of breast cancer is still a challenge. METHODS: We analyzed secretion of the protein cellular communication network factor 1 (CCN1, formerly cysteine-rich angiogenic inducer 61) in breast cancer cell lines by an enzyme-linked immunosorbent assay (ELISA). Soluble CCN1 in the plasma (2.5 µL) of 544 patients with breast cancer and 427 healthy controls was analyzed by ELISA. The breast cancer samples were acquired at the time of primary diagnosis prior to neoadjuvant therapy or surgery. A classifier was established on a training cohort of patients with breast cancer and age-adapted healthy controls and further validated on an independent cohort comprising breast cancer patients and healthy controls. Samples from patients with benign breast diseases were investigated as additional controls. Samples from patients with acute heart diseases (n = 127) were investigated as noncancer controls. The diagnostic accuracy was determined by receiver operating characteristic using the parameters area under the curve, sensitivity, and specificity. RESULTS: CCN1 was frequently secreted by breast cancer cell lines into the extracellular space. Subsequent analysis of clinical blood samples from patients with breast cancer and age-adjusted healthy controls revealed an overall specificity of 99.0% and sensitivity of 80.0% for cancer detection. Remarkably, 81.5% of small T1 cancers were already CCN1-positive, while CCN1 concentrations in patients with benign breast lesions were below the threshold for breast cancer detection. CONCLUSIONS: Circulating CCN1 is a potentially novel blood biomarker for the detection of breast cancer at the earliest invasive stage.

Investigation of the Proteomes of the Truffles Tuber albidum pico, T. aestivum, T. indicum, T. magnatum, and T. melanosporum.

Truffles of the Tuber species are known as expensive foods, mainly for their distinct aroma and taste. This high price makes them a profitable target of food fraud, e.g., the misdeclaration of cheaper truffle species as expensive ones. While many studies investigated truffles on the metabolomic level or the volatile organic compounds extruded by them, research at the proteome level as a phenotype determining basis is limited. In this study, a bottom-up proteomic approach based on LC-MS/MS measurements in data-independent acquisition mode was performed to analyze the truffle species Tuber aestivum, Tuber albidum pico, Tuber indicum, Tuber magnatum, and Tuber melanosporum, and a protein atlas of the investigated species was obtained. The yielded proteomic fingerprints are unique for each of the of the five truffle species and can now be used in case of suspected food fraud. First, a comprehensive spectral library containing 9000 proteins and 50,000 peptides was generated by two-dimensional liquid chromatography coupled to tandem mass spectrometry (2D-LC-MS/MS). Then, samples of the truffle species were analyzed in data-independent acquisition (DIA) proteomics mode yielding 2715 quantified proteins present in all truffle samples. Individual species were clearly distinguishable by principal component analysis (PCA). Quantitative proteome fingerprints were generated from 2066 ANOVA significant proteins, and side-by-side comparisons of truffles were done by T-tests. A further aim of this study was the annotation of functions for the identified proteins. For Tuber magnatum and Tuber melanosporum conclusive links to their superior aroma were found by enrichment of proteins responsible for sulfur-metabolic processes in comparison with other truffles. The obtained data in this study may serve as a reference library for food analysis laboratories in the future to tackle food fraud by misdeclaration of truffles. Further identified proteins with their corresponding abundance values in the different truffle species may serve as potential protein markers in the establishment of targeted analysis methods. Lastly, the obtained data may serve in the future as a basis for deciphering the biochemistry of truffles more deeply as well, when protein databases of the different truffle species will be more complete.

Co-activation of Sonic hedgehog and Wnt signaling in murine retinal precursor cells drives ocular lesions with features of intraocular medulloepithelioma.

Intraocular medulloepithelioma (IO-MEPL) is a rare embryonal ocular neoplasm, prevalently occurring in children. IO-MEPLs share histomorphological features with CNS embryonal tumors with multilayered rosettes (ETMRs), referred to as intracranial medulloepitheliomas. While Sonic hedgehog (SHH) and WNT signaling pathways are crucial for ETMR pathogenesis, the impact of these pathways on human IO-MEPL development is unclear. Gene expression analyses of human embryonal tumor samples revealed similar gene expression patterns and significant overrepresentation of SHH and WNT target genes in both IO-MEPL and ETMR. In order to unravel the function of Shh and Wnt signaling for IO-MEPL pathogenesis in vivo, both pathways were activated in retinal precursor cells in a time point specific manner. Shh and Wnt co-activation in early Sox2- or Rax-expressing precursor cells resulted in infiltrative ocular lesions that displayed extraretinal expansion. Histomorphological, immunohistochemical, and molecular features showed a strong concordance with human IO-MEPL. We demonstrate a relevant role of WNT and SHH signaling in IO-MEPL and report the first mouse model to generate tumor-like lesions with features of IO-MEPL. The presented data may be fundamental for comprehending IO-MEPL initiation and developing targeted therapeutic approaches.

Tissue Sampling and Homogenization in the Sub-Microliter Scale with a Nanosecond Infrared Laser (NIRL) for Mass Spectrometric Proteomics.

It was recently shown that ultrashort pulse infrared (IR) lasers, operating at the wavelength of the OH vibration stretching band of water, are highly efficient for sampling and homogenizing biological tissue. In this study we utilized a tunable nanosecond infrared laser (NIRL) for tissue sampling and homogenization with subsequent liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis for mass spectrometric proteomics. For the first time, laser sampling was performed with murine spleen and colon tissue. An ablation volume of 1.1 × 1.1 × 0.4 mm³ (approximately 0.5 µL) was determined with optical coherence tomography (OCT). The results of bottom-up proteomics revealed proteins with significant abundance differences for both tissue types, which are in accordance with the corresponding data of the Human Protein Atlas. The results demonstrate that tissue sampling and homogenization of small tissue volumes less than 1 µL for subsequent mass spectrometric proteomics is feasible with a NIRL.

An Integrated Strategy Reveals Complex Glycosylation of Erythropoietin Using Mass Spectrometry.

The characterization of therapeutic glycoproteins is challenging due to the structural heterogeneity of the therapeutic protein glycosylation. This study presents an in-depth analytical strategy for glycosylation of first-generation erythropoietin (epoetin beta), including a developed mass spectrometric workflow for N-glycan analysis, bottom-up mass spectrometric methods for site-specific N-glycosylation, and a LC-MS approach for O-glycan identification. Permethylated N-glycans, peptides, and enriched glycopeptides of erythropoietin were analyzed by nanoLC-MS/MS, and de-N-glycosylated erythropoietin was measured by LC-MS, enabling the qualitative and quantitative analysis of glycosylation and different glycan modifications (e.g., phosphorylation and O-acetylation). The newly developed Python scripts enabled the identification of 140 N-glycan compositions (237 N-glycan structures) from erythropoietin, especially including 8 phosphorylated N-glycan species. The site-specificity of N-glycans was revealed at the glycopeptide level by pGlyco software using different proteases. In total, 114 N-glycan compositions were identified from glycopeptide analysis. Moreover, LC-MS analysis of de-N-glycosylated erythropoietin species identified two O-glycan compositions based on the mass shifts between non-O-glycosylated and O-glycosylated species. Finally, this integrated strategy was proved to realize the in-depth glycosylation analysis of a therapeutic glycoprotein to understand its pharmacological properties and improving the manufacturing processes.