Umeå University's logo

Multiclass data sets and large-scale studies are increasingly common in omics sciences, drug discovery, and clinical research due to advancements in analytical platforms. Efficiently handling these data sets and discerning subtle differences across multiple classes remains a significant challenge. In metabolomics, two-class orthogonal projection to latent structures discriminant analysis (OPLS-DA) models are widely used due to their strong discrimination capabilities and ability to provide interpretable information on class differences. However, these models face challenges in multiclass settings. A common solution is to transform the multiclass comparison into multiple two-class comparisons, which, while more effective than a global multiclass OPLS-DA model, unfortunately results in a manual, time-consuming model-building process with complicated interpretation. Here, we introduce an extension of OPLS-DA for data-driven multiclass classification: orthogonal partial least squares-hierarchical discriminant analysis (OPLS-HDA). OPLS-HDA integrates hierarchical cluster analysis (HCA) with the OPLS-DA framework to create a decision tree, addressing multiclass classification challenges and providing intuitive visualization of interclass relationships. To avoid overfitting and ensure reliable predictions, we use cross-validation during model building. Benchmark results show that OPLS-HDA performs competitively across diverse data sets compared to eight established methods. This method represents a significant advancement, offering a powerful tool to dissect complex multiclass data sets. With its versatility, interpretability, and ease of use, OPLS-HDA is an efficient approach to multiclass data analysis applicable across various fields.

Biopharmaceuticals are medical compounds derived from biological sources and are often manufactured by living cells, primarily Chinese hamster ovary (CHO) cells. CHO cells display variation among cell clones, leading to growth and productivity differences that influence the product's quantity and quality. The biological and environmental factors behind these differences are not fully understood. To identify metabolites with a consistent relationship to productivity or cell death over time, we analyzed the extracellular metabolome of 11 CHO clones with different growth and productivity characteristics over 14 days. However, in bioreactor processes, metabolic profiles and process variables are both strongly time-dependent, confounding the metabolite-process variable relationship. To address this, we customized an existing hierarchical approach for handling time dependency to highlight metabolites with a consistent correlation to a process variable over a selected timeframe. We benchmarked this new method against conventional orthogonal partial least squares (OPLS) models. Our hierarchical method highlighted several metabolites consistently related to productivity or cell death that the conventional method missed. These metabolites were biologically relevant; most were known already, but some that had not been reported in CHO literature before, such as 3-methoxytyrosine and succinyladenosine, had ties to cell death in studies with other cell types. The metabolites showed an inverse relationship with the response variables: those positively correlated with productivity were typically negatively correlated with the death rate, or vice versa. For both productivity and cell death, the citrate cycle and adjacent pathways (pyruvate, glyoxylate, pantothenate) were among the most important. In summary, we have proposed a new method to analyze time-dependent omics data in bioprocess production. This approach allowed us to identify metabolites tied to cell death and productivity that were not detected with traditional models.

Cell Painting is an established community-based microscopy-assay platform that provides high-throughput, high-content data for biological readouts. In November 2022, the JUMP-Cell Painting Consortium released the largest publicly available Cell Painting dataset with CellProfiler features, comprising more than 2 billion cell images. This dataset is designed for predicting the activity and toxicity of 115k drug compounds, with the aim to make cell images as computable as genomes and transcriptomes. In this context, our paper introduces a scalable and computationally efficient data analytics workflow created to meet the needs of researchers. This data-driven workflow facilitates the comparison of drug treatment effects through significant and biologically relevant insights. The workflow consists of two parts: first, the Equivalence score (Eq. score), a straightforward yet sophisticated metric highlighting relevant deviations from negative controls based on cell image morphology; second, the scalability of the workflow, by utilizing the Eq. scores on a large scale to predict and classify the subtle morphological changes in cell image profiles. By doing so, we show classification improvements compared to using the raw CellProfiler features on the CPJUMP1-pilot dataset on three types of perturbations. We hope that our workflow's contributions will enhance drug screening efficiency and streamline the drug development process. As this process is resource-intensive, every incremental improvement is valuable. Through our collective efforts in advancing the understanding of high-throughput image-based data, we aim to reduce both the time and cost of developing new, life-saving treatments.

Background: Pancreatic ductal adenocarcinoma (pancreatic cancer) is often detected at late stages resulting in poor overall survival. To improve survival, more patients need to be diagnosed early when curative surgery is feasible. We aimed to identify circulating metabolites that could be used as early pancreatic cancer biomarkers.

Methods: We performed metabolomics by liquid and gas chromatography-mass spectrometry in plasma samples from 82 future pancreatic cancer patients and 82 matched healthy controls within the Northern Sweden Health and Disease Study (NSHDS). Logistic regression was used to assess univariate associations between metabolites and pancreatic cancer risk. Least absolute shrinkage and selection operator (LASSO) logistic regression was used to design a metabolite-based risk score. We used receiver operating characteristic (ROC) analyses to assess the discriminative performance of the metabolite-based risk score.

Results: Among twelve risk-associated metabolites with a nominal P value <0.05, we defined a risk score of three metabolites [indoleacetate, 3-hydroxydecanoate (10:0-OH), and retention index (RI): 2,745.4] using LASSO. A logistic regression model containing these three metabolites, age, sex, body mass index (BMI), smoking status, sample date, fasting status, and carbohydrate antigen 19-9 (CA 19-9) yielded an internal area under curve (AUC) of 0.784 [95% confidence interval (CI): 0.714–0.854] compared to 0.681 (95% CI: 0.597–0.764) for a model without these metabolites (P value =0.007). Seventeen metabolites were significantly associated with pancreatic cancer survival [false discovery rate (FDR) <0.1].

Conclusions: Indoleacetate, 3-hydroxydecanoate (10:0-OH), and RI: 2,745.4 were identified as the top candidate biomarkers for early detection. However, continued efforts are warranted to determine the usefulness of these metabolites as early pancreatic cancer biomarkers.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease with poor survival. Novel biomarkers are urgently needed to improve the outcome through early detection. Here, we aimed to discover novel biomarkers for early PDAC detection using multi-omics profiling in pre-diagnostic plasma samples biobanked after routine health examinations.

A nested case-control study within the Northern Sweden Health and Disease Study was designed. Pre-diagnostic plasma samples from 37 future PDAC patients collected within 2.3 years before diagnosis and 37 matched healthy controls were included. We analyzed metabolites using liquid chromatography mass spectrometry and gas chromatography mass spectrometry, microRNAs by HTG edgeseq, proteins by multiplex proximity extension assays, as well as three clinical biomarkers using milliplex technology. Supervised and unsupervised multi-omics integration were performed as well as univariate analyses for the different omics types and clinical biomarkers. Multiple hypothesis testing was corrected using Benjamini-Hochberg's method and a false discovery rate (FDR) below 0.1 was considered statistically significant.

Carbohydrate antigen (CA) 19-9 was associated with PDAC risk (OR [95 % CI] = 3.09 [1.31–7.29], FDR = 0.03) and increased closer to PDAC diagnosis. Supervised multi-omics models resulted in poor discrimination between future PDAC cases and healthy controls with obtained accuracies between 0.429–0.500. No single metabolite, microRNA, or protein was differentially altered (FDR < 0.1) between future PDAC cases and healthy controls.

CA 19-9 levels increase up to two years prior to PDAC diagnosis but extensive multi-omics analysis including metabolomics, microRNAomics and proteomics in this cohort did not identify novel early biomarkers for PDAC.

Purpose: Peritoneal infection, due to anastomotic leakage, after resection for colorectal cancer have been shown to associate with increased cancer recurrence and mortality, as well as cardiovascsular morbidity. Alterations in circulating protein levels could help shed light on the underlying mechanisms, prompting this exploratory study of 64 patients operated for colorectal cancer with anastomosis. Methods: Thirty-two cases who suffered a postoperative peritoneal infection were matched with 32 controls who had a complication-free postoperative stay. Proteins in serum samples at their first postoperative visit and at one year after surgery were analysed using proximity extension assays and enzyme-linked immunosorbent assays. Multivariate projection methods, adjusted for multiple testing, were used to compare levels between groups, and enrichment and network analyses were performed. Results: Seventy-seven proteins, out of 270 tested, were differentially expressed at a median sampling time of 41 days postoperatively. These proteins were all normalised one year after surgery. Many of the differentially expressed top hub proteins have known involvement in cancer progression, survival, invasiveness and metastasis. Over-represented pathways were related to cardiomyopathy, cell-adhesion, extracellular matrix, phosphatidylinositol-3-kinase/Akt (PI3K-Akt) and transforming growth factor beta (TGF-β) signaling. Conclusion: These affected proteins and pathways could provide clues as to why patients with peritoneal infection might suffer increased cancer recurrence, mortality and cardiovascular morbidity.

Pancreatic ductal adenocarcinoma (PDAC) is a major cause of cancer death that typically presents at an advanced stage. No reliable markers for early detection presently exist. The prominent tumor stroma represents a source of circulating biomarkers for use together with cancer cell-derived biomarkers for earlier PDAC diagnosis. CA19-9 and CEA (cancer cell-derived biomarkers), together with endostatin and collagen IV (stroma-derived) were examined alone, or together, by multivariable modelling, using pre-diagnostic plasma samples (n = 259 samples) from the Northern Sweden Health and Disease Study biobank. Serial samples were available for a subgroup of future patients. Marker efficacy for future PDAC case prediction (n = 154 future cases) was examined by both cross-sectional (ROC analysis) and longitudinal analyses. CA19-9 performed well at, and within, six months to diagnosis and multivariable modelling was not superior to CA19-9 alone in cross-sectional analysis. Within six months to diagnosis, CA19-9 (AUC = 0.92) outperformed the multivariable model (AUC = 0.81) at a cross-sectional level. At diagnosis, CA19-9 (AUC = 0.995) and the model (AUC = 0.977) performed similarly. Longitudinal analysis revealed increases in CA19-9 up to two years to diagnosis which indicates a window of opportunity for early detection of PDAC.

BACKGROUND: Gliomas are complex tumors with several genetic aberrations and diverse metabolic programs contributing to their aggressive phenotypes and poor prognoses. This study defines key metabolic features that can be used to differentiate between glioma subtypes, with potential for improved diagnostics and subtype targeted therapy.

METHODS: Cross-platform global metabolomic profiling coupled with clinical, genetic, and pathological analysis of glioma tissue from 224 tumors - oligodendroglioma (n=31), astrocytoma (n=31) and glioblastoma (n=162) - were performed. Identified metabolic phenotypes were evaluated in accordance with the WHO classification, IDH-mutation, 1p/19q-codeletion, WHO-grading 2-4, and MGMT promoter methylation.

RESULTS: Distinct metabolic phenotypes separate all six analyzed glioma subtypes. IDH-mutated subtypes, expressing 2-hydroxyglutaric acid, were clearly distinguished from IDH-wildtype subtypes. Considerable metabolic heterogeneity outside of the mutated IDH pathway were also evident, with key metabolites being high expression of glycerophosphates, inositols, monosaccharides and sugar alcohols and low levels of sphingosine and lysoglycerophospholipids in IDH-mutants. Among the IDH-mutated subtypes, we observed high levels of amino acids, especially glycine and 2-aminoadipic acid, in grade 4 glioma, and N-acetyl aspartic acid in low-grade astrocytoma and oligodendroglioma. Both IDH-wildtype and mutated oligodendroglioma and glioblastoma were characterized by high levels of acylcarnitines, likely driven by rapid cell growth and hypoxic features. We found elevated levels of 5-HIAA in gliosarcoma and a subtype of oligodendroglioma not yet defined as a specific entity, indicating a previously not described role for the serotonin pathway linked to glioma with bimorphic tissue.

CONCLUSION: Key metabolic differences exist across adult glioma subtypes.

BACKGROUND: Colorectal anastomotic leakage can be considered a process of failed wound healing, for which related biomarkers might be a promising research area to decrease leak rates.

METHODS: Patients who had elective surgery with a primary anastomosis for non-metastatic colorectal cancer, at two university hospitals between 1 January 2010 and 31 December 2015 were included. Patients with an anastomotic leak were identified and matched (1:1) to complication-free controls on the basis of sex, age, tumour stage, tumour location, and operating hospital. Preoperative blood samples were analysed by use of protein panels associated with systemic or enteric inflammation by proteomics, and enzyme-linked immunosorbent assays. Multivariable projection methods were used in the statistical analyses and adjusted for multiple comparisons to reduce false positivity. Rectal cancer tissue samples were evaluated with immunohistochemistry to determine local expression of biomarkers that differed significantly between cases and controls.

RESULTS: Out of 726 patients undergoing resection, 41 patients with anastomotic leakage were matched to 41 controls. Patients with rectal cancer with leakage displayed significantly elevated serum levels of 15 proteins related to inflammation. After controlling for a false discovery rate, levels of C-X-C motif chemokine 6 (CXCL6) and C-C motif chemokine 11 (CCL11) remained significant. In patients with colonic cancer with leakage, levels of high-sensitivity C-reactive protein (hs-CRP) were increased before surgery. Local expression of CXCL6 and CCL11, and their receptors, were similar in rectal tissues between cases and controls.

CONCLUSION: Patients with anastomotic leakage could have an upregulated inflammatory response before surgery, as expressed by elevated serological levels of CXCL6 and CCL11 for rectal cancer and hs-CRP levels in patients with colonic cancer respectively.

Circulating type IV collagen (cCOL IV) is a potential biomarker for patients with colorectal liver metastases (CLM) who present with elevated levels of COL IV in both CLM tissue and circulation. This study aimed to establish the cellular origin of elevated levels of COL IV and analyze circulating COL IV in CLM patients. The cellular source was established through in situ hybridization, immunohistochemical staining, and morphological evaluation. Cellular expression in vitro was assessed by immunofluorescence. Tissue expression of COL IV-degrading matrix metalloproteinases (MMPs)-2, -7, -9, and -13 was studied with immunohistochemical staining. Plasma levels of COL IV in CLM patients and healthy controls were analyzed with ELISA. This study shows that cancer-associated fibroblasts (CAFs) express COL IV in the stroma of CLM and that COL IV is expressed in vitro by fibroblasts but not by tumor cells. MMP-2, -7, -9, and -13 are expressed in CLM tissue, mainly by hepatocytes and immune cells, and circulating COL IV is significantly elevated in CLM patients compared with healthy controls. Our study shows that stromal cells, not tumor cells, produce COL IV in CLM, and that circulating COL IV is elevated in patients with CLM.