NOREVA: normalization and evaluation of MS-based metabolomics data

Explainable machine learning model for identifying key gut microbes and metabolites biomarkers associated with myasthenia gravis

Diagnostic markers for myasthenia gravis (MG) are limited; thus, innovative approaches are required for supportive diagnosis and personalized care. Gut microbes are associated with MG pathogenesis; however, few studies have adopted machine learning (ML) to identify the associations among MG, gut microbiota, and metabolites. In this study, we developed an explainable ML model to predict biomarkers for MG diagnosis. We enrolled 19 MG patients and 10 non-MG individuals. Stool samples were collected and microbiome assessment was performed using 16S rRNA sequencing. Untargeted metabolic profiling was conducted to identify fecal amplicon significant variants (ASVs) and metabolites. We developed an explainable ML model in which the top ASVs and metabolites are combined to identify the best predictive performance. This model uses the SHapley Additive exPlanations method to generate both global and personalized explanations. Fecal microbe-metabolite composition differed significantly between groups. The key bacterial families were Lachnospiraceae and Ruminococcaceae, and the top three features were Lachnospiraceae, inosine, and methylhistidine. An ML model trained with the top 1 % ASVs and top 15 % metabolites combined outperformed all other models. Personalized explanations revealed different patterns of microbe-metabolite contributions in patients with MG. The integration of the microbiota-metabolite features and the development of an explainable ML framework can accurately identify MG and provide personalized explanations, revealing the associations between gut microbiota, metabolites, and MG. An online calculator employing this algorithm was developed that provides a streamlined interface for MG diagnosis screening and conducting personalized evaluations.

Multi-scale DNA language model improves 6 mA binding sites prediction

DNA methylation at the N6 position of adenine (N6-methyladenine, 6 mA), which refers to the attachment of a methyl group to the N6 site of the adenine (A) of DNA, is an important epigenetic modification in prokaryotic and eukaryotic genomes. Accurately predicting the 6 mA binding sites can provide crucial insights into gene regulation, DNA repair, disease development and so on. Wet experiments are commonly used for analyzing 6 mA binding sites. However, they suffer from high cost and expensive time. Therefore, various deep learning methods have been widely used to predict 6 mA binding sites recently. In this study, we develop a framework based on multi-scale DNA language model named "iDNA6mA-MDL". "iDNA6mA-MDL" integrates multiple kmers and the nucleotide property and frequency method for feature embedding, which can capture a full range of DNA sequence context information. At the prediction stage, it also leverages DNABERT to compensate for the incomplete capture of global DNA information. Experiments show that our framework obtains average AUC of 0.981 on a classic 6 mA rice gene dataset, going beyond all existing advanced models under fivefold cross-validations. Moreover, "iDNA6mA-MDL" outperforms most of the popular state-of-the-art methods on another 11 6 mA datasets, demonstrating its effectiveness in 6 mA binding sites prediction.

DrugMAP 2.0: molecular atlas and pharma-information of all drugs

The escalating costs and high failure rates have decelerated the pace of drug development, which amplifies the research interests in developing combinatorial/repurposed drugs and understanding off-target adverse drug reaction (ADR). In other words, it is demanded to delineate the molecular atlas and pharma-information for the combinatorial/repurposed drugs and off-target interactions. However, such invaluable data were inadequately covered by existing databases. In this study, a major update was thus conducted to the DrugMAP, which accumulated (a) 20831 combinatorial drugs and their interacting atlas involving 1583 pharmacologically important molecules; (b) 842 repurposed drugs and their interacting atlas with 795 molecules; (c) 3260 off-targets relevant to the ADRs of 2731 drugs and (d) various types of pharmaceutical information, including diverse ADMET properties, versatile diseases, and various ADRs/off-targets. With the growing demands for discovering combinatorial/repurposed therapies and the rapidly emerging interest in AI-based drug discovery, DrugMAP was highly expected to act as an indispensable supplement to existing databases facilitating drug discovery, which was accessible at: https://idrblab.org/drugmap/.

Application of artificial intelligence in drug design: A review

Artificial intelligence (AI) is a field of computer science that involves acquiring information, developing rule bases, and mimicking human behaviour. The fundamental concept behind AI is to create intelligent computer systems that can operate with minimal human intervention or without any intervention at all. These rule-based systems are developed using various machine learning and deep learning models, enabling them to solve complex problems. AI is integrated with these models to learn, understand, and analyse provided data. The rapid advancement of Artificial Intelligence (AI) is reshaping numerous industries, with the pharmaceutical sector experiencing a notable transformation. AI is increasingly being employed to automate, optimize, and personalize various facets of the pharmaceutical industry, particularly in pharmacological research. Traditional drug development methods areknown for being time-consuming, expensive, and less efficient, often taking around a decade and costing billions of dollars. The integration of artificial intelligence (AI) techniques addresses these challenges by enabling the examination of compounds with desired properties from a vast pool of input drugs. Furthermore, it plays a crucial role in drug screening by predicting toxicity, bioactivity, ADME properties (absorption, distribution, metabolism, and excretion), physicochemical properties, and more. AI enhances the drug design process by improving the efficiency and accuracy of predicting drug behaviour, interactions, and properties. These approaches further significantly improve the precision of drug discovery processes and decrease clinical trial costs leading to the development of more effective drugs.

Multi-scale variational autoencoder for imputation of missing values in untargeted metabolomics using whole-genome sequencing data

BACKGROUND

Missing data is a common challenge in mass spectrometry-based metabolomics, which can lead to biased and incomplete analyses. The integration of whole-genome sequencing (WGS) data with metabolomics data has emerged as a promising approach to enhance the accuracy of data imputation in metabolomics studies.

METHOD

In this study, we propose a novel method that leverages the information from WGS data and reference metabolites to impute unknown metabolites. Our approach utilizes a multi-scale variational autoencoder to jointly model the burden score, polygenetic risk score (PGS), and linkage disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature extraction and missing metabolomics data imputation. By learning the latent representations of both omics data, our method can effectively impute missing metabolomics values based on genomic information.

RESULTS

We evaluate the performance of our method on empirical metabolomics datasets with missing values and demonstrate its superiority compared to conventional imputation techniques. Using 35 template metabolites derived burden scores, PGS and LD-pruned SNPs, the proposed methods achieved R2-scores > 0.01 for 71.55 % of metabolites.

CONCLUSION

The integration of WGS data in metabolomics imputation not only improves data completeness but also enhances downstream analyses, paving the way for more comprehensive and accurate investigations of metabolic pathways and disease associations. Our findings offer valuable insights into the potential benefits of utilizing WGS data for metabolomics data imputation and underscore the importance of leveraging multi-modal data integration in precision medicine research.

Emotion analysis and recognition in 3D space using classifier-dependent feature selection in response to tactile enhanced audio-visual content using EEG

Traditional media such as text, images, audio, and video primarily target specific senses like vision and hearing. In contrast, multiple sensorial media aims to create immersive experiences by integrating additional sensory modalities such as touch, smell, and taste where applicable. Tactile enhanced audio-visual content leverages the sense of touch in addition to visual and auditory stimuli, aiming to create a more immersive and engaging interaction for users. Previously, tactile enhanced content has been explored in 2D emotional space (valence and arousal). In this paper, EEG data against tactile enhanced audio-visual content is labeled based on a self-assessment manikin scale in 3 dimensions i.e., valence, arousal, and dominance. Statistical significance (with a 95% confidence interval) is also established based on gathered scores, highlighting a significant difference in the arousal and dominance dimension of traditional media and tactile enhanced media. A new methodology is proposed using classifier-dependent feature selection approach to classify valence, arousal, and dominance states using three different classifiers. A highest accuracy of 75%, 73.8%, and 75% is achieved for classifying valence, arousal, and dominance states, respectively. The proposed scheme outperforms previous emotion recognition based studies in response to enhanced multimedia content in terms of accuracy, F-score, and other error parameters.

Antiprotozoal peptide prediction using machine learning with effective feature selection techniques

Background

Protozoal pathogens pose a considerable threat, leading to notable mortality rates and the ongoing challenge of developing resistance to drugs. This situation underscores the urgent need for alternative therapeutic approaches. Antimicrobial peptides stand out as promising candidates for drug development. However, there is a lack of published research focusing on predicting antimicrobial peptides specifically targeting protozoal pathogens. In this study, we introduce a successful machine learning-based framework designed to predict potential antiprotozoal peptides effective against protozoal pathogens.

Objective

The primary objective of this study is to classify and predict antiprotozoal peptides using diverse negative datasets.

Methods

A comprehensive literature review was conducted to gather experimentally validated antiprotozoal peptides, forming the positive dataset for our study. To construct a robust machine learning classifier, multiple negative datasets were incorporated, including (i) non-antimicrobial, (ii) antiviral, (iii) antibacterial, (iv) antifungal, and (v) antimicrobial peptides excluding those targeting protozoal pathogens. Various compositional features of the peptides were extracted using the pfeature algorithm. Two feature selection methods, SVC-L1 and mRMR, were employed to identify highly relevant features crucial for distinguishing between the positive and negative datasets. Additionally, five popular classifiers i.e. Decision Tree, Random Forest, Support Vector Machine, Logistic Regression, and XGBoost were used to build efficient decision models.

Results

XGBoost was the most effective in classifying antiprotozoal peptides from each negative dataset based on the features selected by the mRMR feature selection method. The proposed machine learning framework efficiently differentiate the antiprotozoal peptides from (i) non-antimicrobial (ii) antiviral (iii) antibacterial (iv) antifungal and (v) antimicrobial with accuracy of 97.27 %, 93.64 %, 86.36 %, 90.91 %, and 89.09 % respectively on the validation dataset.

Conclusion

The models are incorporated in a user-friendly web server (www.soodlab.com/appred) to predict the antiprotozoal activity of given peptides.

Differential immunometabolic responses to Delta and Omicron SARS-CoV-2 variants in golden syrian hamsters

Delta (B.1.617.2) and Omicron (B.1.1.529) variants of SARS-CoV-2 represents unique clinical characteristics. However, their role in altering immunometabolic regulations during acute infection remains convoluted. Here, we evaluated the differential immunopathogenesis of Delta vs. Omicron variants in Golden Syrian hamsters (GSH). The Delta variant resulted in higher virus titers in throat swabs and the lungs and exhibited higher lung damage with immune cell infiltration than the Omicron variant. The gene expression levels of immune mediators and metabolic enzymes, Arg-1 and IDO1 in the Delta-infected lungs were significantly higher compared to Omicron. Further, Delta/Omicron infection perturbed carbohydrates, amino acids, nucleotides, and TCA cycle metabolites and was differentially regulated compared to uninfected lungs. Collectively, our data provide a novel insight into immunometabolic/pathogenic outcomes for Delta vs. Omicron infection in the GSH displaying concordance with COVID-19 patients associated with inflammation and tissue injury during acute infection that offered possible new targets to develop potential therapeutics.

Targeted saliva metabolomics in Sjögren's syndrome

OBJECTIVE

Sjögren's Syndrome (SS) is a chronic inflammatory autoimmune exocrinopathy, and although, the role of metabolism in the autoimmune responses has been discussed in diseases such as lupus erythematosus, rheumatoid arthritis, psoriasis and scleroderma. There is a lack of information regarding the metabolic implications of SS. Considering that the disease affects primarily salivary glands; the aim of this study is to evaluate the metabolic changes in the salivary glands' microenvironment using a targeted metabolomics approach.

METHODS

The saliva from 10 patients diagnosed with SS by the American-European consensus and 10 healthy volunteers was analyzed in an Ultra-high Performance Liquid Chromatograph Coupled Mass Spectrometry (UPLC-MS).

RESULTS

The results showed an increased concentration in SS of metabolites involved in oxidative stress such as lactate, alanine and malate, and amino acids involved in the growth and proliferation of T-cells, such as arginine, leucine valine and isoleucine.

CONCLUSIONS

These results revealed that is possible to differentiate the metabolic profile of SS and healthy individuals using a small amount of saliva, which in its turn may reflect the cellular changes observed in the microenvironments of damaged salivary glands from these patients.

Ablation of Sam50 is associated with fragmentation and alterations in metabolism in murine and human myotubes

The sorting and assembly machinery (SAM) Complex is responsible for assembling β-barrel proteins in the mitochondrial membrane. Comprising three subunits, Sam35, Sam37, and Sam50, the SAM complex connects the inner and outer mitochondrial membranes by interacting with the mitochondrial contact site and cristae organizing system complex. Sam50, in particular, stabilizes the mitochondrial intermembrane space bridging (MIB) complex, which is crucial for protein transport, respiratory chain complex assembly, and regulation of cristae integrity. While the role of Sam50 in mitochondrial structure and metabolism in skeletal muscle remains unclear, this study aims to investigate its impact. Serial block-face-scanning electron microscopy and computer-assisted 3D renderings were employed to compare mitochondrial structure and networking in Sam50-deficient myotubes from mice and humans with wild-type (WT) myotubes. Furthermore, autophagosome 3D structure was assessed in human myotubes. Mitochondrial metabolic phenotypes were assessed using Gas Chromatography-Mass Spectrometry-based metabolomics to explore differential changes in WT and Sam50-deficient myotubes. The results revealed increased mitochondrial fragmentation and autophagosome formation in Sam50-deficient myotubes compared to controls. Metabolomic analysis indicated elevated metabolism of propanoate and several amino acids, including ß-Alanine, phenylalanine, and tyrosine, along with increased amino acid and fatty acid metabolism in Sam50-deficient myotubes. Furthermore, impairment of oxidative capacity was observed upon Sam50 ablation in both murine and human myotubes, as measured with the XF24 Seahorse Analyzer. Collectively, these findings support the critical role of Sam50 in establishing and maintaining mitochondrial integrity, cristae structure, and mitochondrial metabolism. By elucidating the impact of Sam50-deficiency, this study enhances our understanding of mitochondrial function in skeletal muscle.

An integrated in silico approach for the identification of novel potential drug target and chimeric vaccine against Neisseria meningitides strain 331401 serogroup X by subtractive genomics and reverse vaccinology

Neisseria meningitidis, commonly known as the meningococcus, leads to substantial illness and death among children and young adults globally, revealing as either epidemic or sporadic meningitis and/or septicemia. In this study, we have designed a novel peptide-based chimeric vaccine candidate against the N. meningitidis strain 331,401 serogroup X. Through rigorous analysis of subtractive genomics, two essential cytoplasmic proteins, namely UPI000012E8E0(UDP-3-O-acyl-GlcNAc deacetylase) and UPI0000ECF4A9(UDP-N-acetylglucosamine acyltransferase) emerged as potential drug targets. Additionally, using reverse vaccinology, the outer membrane protein UPI0001F4D537 (Membrane fusion protein MtrC) identified by subcellular localization and recognized for its known indispensable role in bacterial survival was identified as a novel chimeric vaccine target. Following a careful comparison of MHC-I, MHC-II, T-cell, and B-cell epitopes, three epitopes derived from UPI0001F4D537 were linked with three types of linkers-GGGS, EAAAK, and the essential PADRE-for vaccine construction. This resulted in eight distinct vaccine models (V1-V8). Among them V1 model was selected as the final vaccine construct. It exhibits exceptional immunogenicity, safety, and enhanced antigenicity, with 97.7 % of its residues in the Ramachandran plot's most favored region. Subsequently, the vaccine structure was docked with the TLR4/MD2 complex and six different HLA allele receptors using the HADDOCK server. The docking resulted in the lowest HADDOCK score of 39.3 ± 9.0 for TLR/MD2. Immune stimulation showed a strong immune response, including antibodies creation and the activation of B-cells, T Cytotoxic cells, T Helper cells, Natural Killer cells, and interleukins. Furthermore, the vaccine construct was successfully expressed in the Escherichia coli system by reverse transcription, optimization, and ligation in the pET-28a (+) vector for the expression study. The current study proposes V1 construct has the potential to elicit both cellular and humoral responses, crucial for the developing an epitope-based vaccine against N. meningitidis strain 331,401 serogroup X.

Identification of a G-protein coupled receptor-related gene signature through bioinformatics analysis to construct a risk model for ovarian cancer prognosis

BACKGROUND

Ovarian cancer (OV) is a common malignant tumor of the female reproductive system with a 5-year survival rate of ∼30 %. Inefficient early diagnosis and prognosis leads to poor survival in most patients. G protein-coupled receptors (GPCRs, the largest family of human cell surface receptors) are associated with OV. We aimed to identify GPCR-related gene (GPCRRG) signatures and develop a novel model to predict OV prognosis.

METHOD

We downloaded data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Prognostic GPCRRGs were screened using least absolute shrinkage and selection operator (LASSO) Cox regression analysis, and a prognostic model was constructed. The predictive ability of the model was evaluated by Kaplan-Meier (K-M) survival analysis. The levels of GPCRRGs were examined in normal and OV cell lines using quantitative reverse-Etranscription polymerase chain reaction. The immunological characteristics of the high- and low-risk groups were analyzed using single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT.

RESULTS

Based on the risks scores, 17 GPCRRGs were associated with OV prognosis. CXCR4, GPR34, LGR6, LPAR3, and RGS2 were significantly expressed in three OV datasets and enabled accurate OV diagnosis. K-M analysis of the prognostic model showed that it could differentiate high- and low-risk patients, which correspond to poorer and better prognoses, respectively. GPCRRG expression was correlated with immune infiltration rates.

CONCLUSIONS

Our prognostic model elaborates on the roles of GPCRRGs in OV and provides a new tool for prognosis and immune response prediction in patients with OV.

scHD4E: Novel ensemble learning-based differential expression analysis method for single-cell RNA-sequencing data

Differential expression (DE) analysis between cell types for scRNA-seq data by capturing its complicated features is crucial. Recently, different methods have been developed for targeting the scRNA-seq data analysis based on different modeling frameworks, assumptions, strategies and test statistic in considering various data features. The scDEA is an ensemble learning-based DE analysis method developed recently, yielding p-values using Lancaster's combination, generated by 12 individual DE analysis methods, and producing more accurate and stable results than individual methods. The objective of our study is to propose a new ensemble learning-based DE analysis method, scHD4E, using top performers in only 4 separate methods. The top performer 4 methods have been selected through an evaluation process using six real scRNA-seq data sets. We conducted comprehensive experiments for five experimental data sets to evaluate our proposed method based on the sample size effects, batch effects, type I error control, gene ontology enrichment analysis, runtime, identified matched DE genes, and semantic similarity measurement between methods. We also perform similar analyses (except the last 3 terms) and compute performance measures like accuracy, F1 score, Mathew's correlation coefficient etc. for a simulated data set. The results show that scHD4E is performs better than all the individual and scDEA methods in all the above perspectives. We expect that scHD4E will serve the modern data scientists for detecting the DEGs in scRNA-seq data analysis. To implement our proposed method, a Github R package scHD4E and its shiny application has been developed, and available in the following links: https://github.com/bbiswas1989/scHD4E and https://github.com/bbiswas1989/scHD4E-Shiny.

A diverse proteome is present and enzymatically active in metabolite extracts

Metabolite extraction is the critical first-step in metabolomics experiments, where it is generally regarded to inactivate and remove proteins. Here, arising from efforts to improve extraction conditions for polar metabolomics, we discover a proteomic landscape of over 1000 proteins within metabolite extracts. This is a ubiquitous feature across several common extraction and sample types. By combining post-resuspension stable isotope addition and enzyme inhibitors, we demonstrate in-extract metabolite interconversions due to residual transaminase activity. We extend these findings with untargeted metabolomics where we observe extensive protein-mediated metabolite changes, including in-extract formation of glutamate dipeptide and depletion of total glutathione. Finally, we present a simple extraction workflow that integrates 3 kDa filtration for protein removal as a superior method for polar metabolomics. In this work, we uncover a previously unrecognized, protein-mediated source of observer effects in metabolomics experiments with broad-reaching implications across all research fields using metabolomics and molecular metabolism.

Enhancing biomass conversion to bioenergy with machine learning: Gains and problems

The growing concerns about environmental sustainability and energy security, such as exhaustion of traditional fossil fuels and global carbon footprint growth have led to an increasing interest in alternative energy sources, especially bioenergy. Recently, numerous scenarios have been proposed regarding the use of bioenergy from different sources in the future energy systems. In this regard, one of the biggest challenges for scientists is managing, modeling, decision-making, and future forecasting of bioenergy systems. The development of machine learning (ML) techniques can provide new opportunities for modeling, optimizing and managing the production, consumption and environmental effects of bioenergy. However, researchers in bioenergy fields have not widely utilized the ML concepts and practices. Therefore, a comparative review of the current ML techniques used for bioenergy productions is presented in this paper. This review summarizes the common issues and difficulties existing in integrating ML with bioenergy studies, and discusses and proposes the possible solutions. Additionally, a detailed discussion of the appropriate ML application scenarios is also conducted in every sector of the entire bioenergy chain. This indicates the modernized conversion processes supported by ML techniques are imperative to accurately capture process-level subtleties, and thus improving techno-economic resilience and socio-ecological integrity of bioenergy production. All the efforts are believed to help in sustainable bioenergy production with ML technologies for the future.

Precision healthcare: A deep dive into machine learning algorithms and feature selection strategies for accurate heart disease prediction

This paper presents a comprehensive exploration of machine learning algorithms (MLAs) and feature selection techniques for accurate heart disease prediction (HDP) in modern healthcare. By focusing on diverse datasets encompassing various challenges, the research sheds light on optimal strategies for early detection. MLAs such as Decision Trees (DT), Random Forests (RF), Support Vector Machines (SVM), Gaussian Naive Bayes (NB), and others were studied, with precision and recall metrics emphasized for robust predictions. Our study addresses challenges in real-world data through data cleaning and one-hot encoding, enhancing the integrity of our predictive models. Feature extraction techniques-Recursive Feature Extraction (RFE), Principal Component Analysis (PCA), and univariate feature selection-play a crucial role in identifying relevant features and reducing data dimensionality. Our findings showcase the impact of these techniques on improving prediction accuracy. Optimized models for each dataset have been achieved through grid search hyperparameter tuning, with configurations meticulously outlined. Notably, a remarkable 99.12 % accuracy was achieved on the first Kaggle dataset, showcasing the potential for accurate HDP. Model robustness across diverse datasets was highlighted, with caution against overfitting. The study emphasizes the need for validation of unseen data and encourages ongoing research for generalizability. Serving as a practical guide, this research aids researchers and practitioners in HDP model development, influencing clinical decisions and healthcare resource allocation. By providing insights into effective algorithms and techniques, the paper contributes to reducing heart disease-related morbidity and mortality, supporting the healthcare community's ongoing efforts.

Biological sex impacts oxidative stress in skeletal muscle in a porcine heat stress model

Oxidative stress contributes to heat stress (HS)-mediated alterations in skeletal muscle; however, the extent to which biological sex mediates oxidative stress during HS remains unknown. We hypothesized muscle from males would be more resistant to oxidative stress caused by HS than muscle from females. To address this, male and female pigs were housed in thermoneutral conditions (TN; 20.8 ± 1.6°C; 62.0 ± 4.7% relative humidity; n = 8/sex) or subjected to HS (39.4 ± 0.6°C; 33.7 ± 6.3% relative humidity) for 1 (HS1; n = 8/sex) or 7 days (HS7; n = 8/sex) followed by collection of the oxidative portion of the semitendinosus. Although HS increased muscle temperature, by 7 days, muscle from heat-stressed females was cooler than muscle from heat-stressed males (0.3°C; P < 0.05). Relative protein abundance of 4-hydroxynonenal (4-HNE)-modified proteins increased in HS1 females compared with TN (P = 0.05). Furthermore, malondialdehyde (MDA)-modified proteins and 8-hydroxy-2'-deoxyguanosine (8-OHdG) concentration, a DNA damage marker, was increased in HS7 females compared with TN females (P = 0.05). Enzymatic activities of catalase and superoxide dismutase (SOD) remained similar between groups; however, glutathione peroxidase (GPX) activity decreased in HS7 females compared with TN and HS1 females (P ≤ 0.03) and HS7 males (P = 0.02). Notably, HS increased skeletal muscle Ca2+ deposition (P = 0.05) and was greater in HS1 females compared with TN females (P < 0.05). Heat stress increased sarco(endo)plasmic reticulum Ca2+ ATPase (SERCA)2a protein abundance (P < 0.01); however, Ca2+ ATPase activity remained similar between groups. Overall, despite having lower muscle temperature, muscle from heat-stressed females had increased markers of oxidative stress and calcium deposition than muscle from males following identical environmental exposure.NEW & NOTEWORTHY Heat stress is a global threat to human health and agricultural production. We demonstrated that following 7 days of heat stress, skeletal muscle from females was more susceptible to oxidative stress than muscle from males in a porcine model, despite cooler muscle temperatures. The vulnerability to heat stress-induced oxidative stress in females may be driven, at least in part, by decreased antioxidant capacity and calcium dysregulation.

Discovery of the potential biomarkers for early diagnosis of endometrial cancer via integrating metabolomics and transcriptomics

Endometrial cancer (EC) is one of the most common malignant tumors in women, and the increasing incidence and mortality pose a serious threat to the public health. Early diagnosis of EC could prolong the survival period and optimize the survivorship, greatly alleviating patients' suffering and social medical pressure. In this study, we collected urine and serum samples from the recruited patients, analyzed the samples using LC-MS approach, and identified the differential metabolites through metabolomic analysis. Then, the differentially expressed genes were identified through the systematic transcriptomic analysis of EC-related dataset from Gene Expression Omnibus (GEO), followed by network profiling of metabolic-reaction-enzyme-gene. In this experiment, a total of 83 differential metabolites and 19 hub genes were discovered, of which 10 different metabolites and 3 hub genes were further evaluated as more potential biomarkers based on network analysis. According to the KEGG enrichment analysis, the potential biomarkers and gene-encoded proteins were found to be involved in the arginine and proline metabolism, histidine metabolism, and pyrimidine metabolism, which was of significance for the early diagnosis of EC. In particular, the combination of metabolites (histamine, 1-methylhistamine, and methylimidazole acetaldehyde) as well as the combination of RRM2, TYMS and TK1 exerted more accurate discrimination abilities between EC and healthy groups, providing more criteria for the early diagnosis of EC.

Characteristics of microbiome-derived metabolomics according to the progression of alcoholic liver disease

BACKGROUND AND AIM

The prevalence and severity of alcoholic liver disease (ALD) are increasing. The incidence of alcohol-related cirrhosis has risen up to 2.5%. This study aimed to identify novel metabolite mechanisms involved in the development of ALD in patients. The use of gut microbiome-derived metabolites is increasing in targeted therapies. Identifying metabolic compounds is challenging due to the complex patterns that have long-term effects on ALD. We investigated the specific metabolite signatures in ALD patients.

METHODS

This study included 247 patients (heathy control, HC: n = 62, alcoholic fatty liver, AFL; n = 25, alcoholic hepatitis, AH; n = 80, and alcoholic cirrhosis, AC, n = 80) identified, and stool samples were collected. 16S rRNA sequencing and metabolomics were performed with MiSeq sequencer and liquid chromatography coupled to time-of-flight-mass spectrometry (LC-TOF-MS), respectively. The untargeted metabolites in AFL, AH, and AC samples were evaluated by multivariate statistical analysis and metabolic pathotypic expression. Metabolic network classifiers were used to predict the pathway expression of the AFL, AH, and AC stages.

RESULTS

The relative abundance of Proteobacteria was increased and the abundance of Bacteroides was decreased in ALD samples (p = 0.001) compared with that in HC samples. Fusobacteria levels were higher in AH samples (p = 0.0001) than in HC samples. Untargeted metabolomics was applied to quantitatively screen 103 metabolites from each stool sample. Indole-3-propionic acid levels are significantly lower in AH and AC (vs. HC, p = 0.001). Indole-3-lactic acid (ILA: p = 0.04) levels were increased in AC samples. AC group showed an increase in indole-3-lactic acid (vs. HC, p = 0.040) level. Compared with that in HC samples, the levels of short-chain fatty acids (SCFAs: acetic acid, butyric acid, propionic acid, iso-butyric acid, and iso-valeric acid) and bile acids (lithocholic acids) were significantly decreased in AC. The pathways of linoleic acid metabolism, indole compounds, histidine metabolism, fatty acid degradation, and glutamate metabolism were closely associated with ALD metabolism.

CONCLUSIONS

This study identified that microbial metabolic dysbiosis is associated with ALD-related metabolic dysfunction. The SCFAs, bile acids, and indole compounds were depleted during ALD progression.

CLINICAL TRIAL

Clinicaltrials.gov, number NCT04339725.

FERREG: ferroptosis-based regulation of disease occurrence, progression and therapeutic response

Ferroptosis is a non-apoptotic, iron-dependent regulatory form of cell death characterized by the accumulation of intracellular reactive oxygen species. In recent years, a large and growing body of literature has investigated ferroptosis. Since ferroptosis is associated with various physiological activities and regulated by a variety of cellular metabolism and mitochondrial activity, ferroptosis has been closely related to the occurrence and development of many diseases, including cancer, aging, neurodegenerative diseases, ischemia-reperfusion injury and other pathological cell death. The regulation of ferroptosis mainly focuses on three pathways: system Xc-/GPX4 axis, lipid peroxidation and iron metabolism. The genes involved in these processes were divided into driver, suppressor and marker. Importantly, small molecules or drugs that mediate the expression of these genes are often good treatments in the clinic. Herein, a newly developed database, named 'FERREG', is documented to (i) providing the data of ferroptosis-related regulation of diseases occurrence, progression and drug response; (ii) explicitly describing the molecular mechanisms underlying each regulation; and (iii) fully referencing the collected data by cross-linking them to available databases. Collectively, FERREG contains 51 targets, 718 regulators, 445 ferroptosis-related drugs and 158 ferroptosis-related disease responses. FERREG can be accessed at https://idrblab.org/ferreg/.

User-friendly relative quantification procedure for gas chromatography/mass spectrometry-based plasma metabolome analysis

RATIONALE

Recently, metabolome analysis has been applied to a variety of research fields, but differences between batches or facilities can cause discrepancies in the results of such analyses. To resolve these issues using comprehensive metabolome analysis, in which it is difficult to perform quantitative analyses of all detected metabolites, internal standard compounds are used to obtain relative metabolite levels. This study investigated gas chromatography/mass spectrometry-based plasma metabolome analysis methods that are superior to relative quantification using internal standard compounds.

METHODS

In experiment I, four analyses were performed under different analytical conditions at one facility, and then the data from the four analyses were compared. In experiment II, the same samples were analyzed at three facilities, and then the data from the three facilities were compared.

RESULTS

Regarding the relative values obtained through comparisons with the internal standard compound, differences in the analytical results were observed among the four analytical conditions in experiment I and among the three facilities in experiment II, and the differences observed among the three facilities (experiment II) were larger. When correction was performed using plasma as a quality control, which is the procedure suggested in this study, these differences were markedly ameliorated.

CONCLUSION

The suggested procedure involves the analysis of a plasma standard as a quality control for each batch and the calculation of relative target plasma to quality-control plasma values for each metabolite. This is an easy and low-cost method and could be readily employed by researchers during comprehensive plasma metabolome analysis.

Rat Model of Type 2 Diabetes Mellitus Recapitulates Human Disease in the Anterior Segment of the Eye

Changes in the anterior segment of the eye due to type 2 diabetes mellitus (T2DM) are not well-characterized, in part due to the lack of a reliable animal model. This study evaluates changes in the anterior segment, including crystalline lens health, corneal endothelial cell density, aqueous humor metabolites, and ciliary body vasculature, in a rat model of T2DM compared with human eyes. Male Sprague-Dawley rats were fed a high-fat diet (45% fat) or normal diet, and rats fed the high-fat diet were injected with streptozotocin i.p. to generate a model of T2DM. Cataract formation and corneal endothelial cell density were assessed using microscopic analysis. Diabetes-related rat aqueous humor alterations were assessed using metabolomics screening. Transmission electron microscopy was used to assess qualitative ultrastructural changes ciliary process microvessels at the site of aqueous formation in the eyes of diabetic rats and humans. Eyes from the diabetic rats demonstrated cataracts, lower corneal endothelial cell densities, altered aqueous metabolites, and ciliary body ultrastructural changes, including vascular endothelial cell activation, pericyte degeneration, perivascular edema, and basement membrane reduplication. These findings recapitulated diabetic changes in human eyes. These results support the use of this model for studying ocular manifestations of T2DM and support a hypothesis postulating blood-aqueous barrier breakdown and vascular leakage at the ciliary body as a mechanism for diabetic anterior segment pathology.

Is fragment-based graph a better graph-based molecular representation for drug design? A comparison study of graph-based models

Graph Neural Networks (GNNs) have gained significant traction in various sectors of AI-driven drug design. Over recent years, the integration of fragmentation concepts into GNNs has emerged as a potent strategy to augment the efficacy of molecular generative models. Nonetheless, challenges such as symmetry breaking and potential misrepresentation of intricate cycles and undefined functional groups raise questions about the superiority of fragment-based graph representation over traditional methods. In our research, we undertook a rigorous evaluation, contrasting the predictive prowess of eight models-developed using deep learning algorithms-across 12 benchmark datasets that span a range of properties. These models encompass established methods like GCN, AttentiveFP, and D-MPNN, as well as innovative fragment-based representation techniques. Our results indicate that fragment-based methodologies, notably PharmHGT, significantly improve model performance and interpretability, particularly in scenarios characterized by limited data availability. However, in situations with extensive training, fragment-based molecular graph representations may not necessarily eclipse traditional methods. In summation, we posit that the integration of fragmentation, as an avant-garde technique in drug design, harbors considerable promise for the future of AI-enhanced drug design.

MultiClassMetabo: A Superior Classification Model Constructed Using Metabolic Markers in Multiclass Metabolomics

Multiclass metabolomics has become a popular technique for revealing the mechanisms underlying certain physiological processes, different tumor types, or different therapeutic responses. In multiclass metabolomics, it is highly important to uncover the underlying biological information on biosamples by identifying the metabolic markers with the most associations and classifying the different sample classes. The classification problem of multiclass metabolomics is more difficult than that of the binary problem. To date, various methods exist for constructing classification models and identifying metabolic markers consisting of well-established techniques and newly emerging machine learning algorithms. However, how to construct a superior classification model using these methods remains unclear for a given multiclass metabolomic data set. Herein, MultiClassMetabo has been developed for constructing a superior classification model using metabolic markers identified in multiclass metabolomics. MultiClassMetabo can enable online services, including (a) identifying metabolic markers by marker identification methods, (b) constructing classification models by classification methods, and (c) performing a comprehensive assessment from multiple perspectives to construct a superior classification model for multiclass metabolomics. In summary, MultiClassMetabo is distinguished for its capability to construct a superior classification model using the most appropriate method through a comprehensive assessment, which makes it an important complement to other available tools in multiclass metabolomics. MultiClassMetabo can be accessed at http://idrblab.cn/multiclassmetabo/.

Buffering Mitigates Chondrocyte Oxidative Stress, Metabolic Dysfunction, and Death Induced by Normal Saline: Formulation of a Novel Arthroscopic Irrigant

For decades, surgeons have utilized 0.9% normal saline (NS) for joint irrigation to improve visualization during arthroscopic procedures. This continues despite mounting evidence that NS exposure impairs chondrocyte metabolism and compromises articular cartilage function. We hypothesized that chondrocyte oxidative stress induced by low pH is the dominant factor driving NS toxicity, and that buffering NS to increase its pH would mitigate these effects. Effects on chondrocyte viability, reactive oxygen species (ROS) production, and overall metabolic function were assessed. Even brief exposure to NS caused cell death, ROS overproduction, and disruption of glycolysis, pentose phosphate, and tricarboxylic acid (TCA) cycle pathways. NS also stimulated ROS overproduction in synovial cells that could adversely alter the synovial function and subsequently the entire joint health. Buffering NS with 25 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) significantly increased chondrocyte viability, reduced ROS production, and returned metabolite levels to near control levels while also reducing ROS production in synovial cells. These results confirm that chondrocytes and synoviocytes are vulnerable to insult from the acidic pH of NS and demonstrate that adding a buffering agent to NS averts many of its most harmful effects.

VARIDT 3.0: the phenotypic and regulatory variability of drug transporter

The phenotypic and regulatory variability of drug transporter (DT) are vital for the understanding of drug responses, drug-drug interactions, multidrug resistances, and so on. The ADME property of a drug is collectively determined by multiple types of variability, such as: microbiota influence (MBI), transcriptional regulation (TSR), epigenetics regulation (EGR), exogenous modulation (EGM) and post-translational modification (PTM). However, no database has yet been available to comprehensively describe these valuable variabilities of DTs. In this study, a major update of VARIDT was therefore conducted, which gave 2072 MBIs, 10 610 TSRs, 46 748 EGRs, 12 209 EGMs and 10 255 PTMs. These variability data were closely related to the transportation of 585 approved and 301 clinical trial drugs for treating 572 diseases. Moreover, the majority of the DTs in this database were found with multiple variabilities, which allowed a collective consideration in determining the ADME properties of a drug. All in all, VARIDT 3.0 is expected to be a popular data repository that could become an essential complement to existing pharmaceutical databases, and is freely accessible without any login requirement at: https://idrblab.org/varidt/.

TheMarker: a comprehensive database of therapeutic biomarkers

Distinct from the traditional diagnostic/prognostic biomarker (adopted as the indicator of disease state/process), the therapeutic biomarker (ThMAR) has emerged to be very crucial in the clinical development and clinical practice of all therapies. There are five types of ThMAR that have been found to play indispensable roles in various stages of drug discovery, such as: Pharmacodynamic Biomarker essential for guaranteeing the pharmacological effects of a therapy, Safety Biomarker critical for assessing the extent or likelihood of therapy-induced toxicity, Monitoring Biomarker indispensable for guiding clinical management by serially measuring patients' status, Predictive Biomarker crucial for maximizing the clinical outcome of a therapy for specific individuals, and Surrogate Endpoint fundamental for accelerating the approval of a therapy. However, these data of ThMARs has not been comprehensively described by any of the existing databases. Herein, a database, named 'TheMarker', was therefore constructed to (a) systematically offer all five types of ThMAR used at different stages of drug development, (b) comprehensively describe ThMAR information for the largest number of drugs among available databases, (c) extensively cover the widest disease classes by not just focusing on anticancer therapies. These data in TheMarker are expected to have great implication and significant impact on drug discovery and clinical practice, and it is freely accessible without any login requirement at: https://idrblab.org/themarker.

Data processing solutions to render metabolomics more quantitative: case studies in food and clinical metabolomics using Metabox 2.0

In classic semiquantitative metabolomics, metabolite intensities are affected by biological factors and other unwanted variations. A systematic evaluation of the data processing methods is crucial to identify adequate processing procedures for a given experimental setup. Current comparative studies are mostly focused on peak area data but not on absolute concentrations. In this study, we evaluated data processing methods to produce outputs that were most similar to the corresponding absolute quantified data. We examined the data distribution characteristics, fold difference patterns between 2 metabolites, and sample variance. We used 2 metabolomic datasets from a retail milk study and a lupus nephritis cohort as test cases. When studying the impact of data normalization, transformation, scaling, and combinations of these methods, we found that the cross-contribution compensating multiple standard normalization (ccmn) method, followed by square root data transformation, was most appropriate for a well-controlled study such as the milk study dataset. Regarding the lupus nephritis cohort study, only ccmn normalization could slightly improve the data quality of the noisy cohort. Since the assessment accounted for the resemblance between processed data and the corresponding absolute quantified data, our results denote a helpful guideline for processing metabolomic datasets within a similar context (food and clinical metabolomics). Finally, we introduce Metabox 2.0, which enables thorough analysis of metabolomic data, including data processing, biomarker analysis, integrative analysis, and data interpretation. It was successfully used to process and analyze the data in this study. An online web version is available at http://metsysbio.com/metabox.

Biological mass spectrometry enables spatiotemporal 'omics: From tissues to cells to organelles

Biological processes unfold across broad spatial and temporal dimensions, and measurement of the underlying molecular world is essential to their understanding. Interdisciplinary efforts advanced mass spectrometry (MS) into a tour de force for assessing virtually all levels of the molecular architecture, some in exquisite detection sensitivity and scalability in space-time. In this review, we offer vignettes of milestones in technology innovations that ushered sample collection and processing, chemical separation, ionization, and 'omics analyses to progressively finer resolutions in the realms of tissue biopsies and limited cell populations, single cells, and subcellular organelles. Also highlighted are methodologies that empowered the acquisition and analysis of multidimensional MS data sets to reveal proteomes, peptidomes, and metabolomes in ever-deepening coverage in these limited and dynamic specimens. In pursuit of richer knowledge of biological processes, we discuss efforts pioneering the integration of orthogonal approaches from molecular and functional studies, both within and beyond MS. With established and emerging community-wide efforts ensuring scientific rigor and reproducibility, spatiotemporal MS emerged as an exciting and powerful resource to study biological systems in space-time.

Evaluating Machine Learning Methods of Analyzing Multiclass Metabolomics

Multiclass metabolomic studies have become popular for revealing the differences in multiple stages of complex diseases, various lifestyles, or the effects of specific treatments. In multiclass metabolomics, there are multiple data manipulation steps for analyzing raw data, which consist of data filtering, the imputation of missing values, data normalization, marker identification, sample separation, classification, and so on. In each step, several to dozens of machine learning methods can be chosen for the given data set, with potentially hundreds or thousands of method combinations in the whole data processing chain. Therefore, a clear understanding of these machine learning methods is helpful for selecting an appropriate method combination for obtaining stable and reliable analytical results of specific data. However, there has rarely been an overall introduction or evaluation of these methods based on multiclass metabolomic data. Herein, detailed descriptions of these machine learning methods in multiple data manipulation steps are reviewed. Moreover, an assessment of these methods was performed using a benchmark data set for multiclass metabolomics. First, 12 imputation methods for imputing missing values were evaluated based on the PSS (Procrustes statistical shape analysis) and NRMSE (normalized root-mean-square error) values. Second, 17 normalization methods for processing multiclass metabolomic data were evaluated by applying the PMAD (pooled median absolute deviation) value. Third, different methods of identifying markers of multiclass metabolomics were evaluated based on the CWrel (relative weighted consistency) value. Fourth, nine classification methods for constructing multiclass models were assessed using the AUC (area under the curve) value. Performance evaluations of machine learning methods are highly recommended to select the most appropriate method combination before performing the final analysis of the given data. Overall, detailed descriptions and evaluation of various machine learning methods are expected to improve analyses of multiclass metabolomic data.

Workflow for Evaluating Normalization Tools for Omics Data Using Supervised and Unsupervised Machine Learning

To achieve high quality omics results, systematic variability in mass spectrometry (MS) data must be adequately addressed. Effective data normalization is essential for minimizing this variability. The abundance of approaches and the data-dependent nature of normalization have led some researchers to develop open-source academic software for choosing the best approach. While these tools are certainly beneficial to the community, none of them meet all of the needs of all users, particularly users who want to test new strategies that are not available in these products. Herein, we present a simple and straightforward workflow that facilitates the identification of optimal normalization strategies using straightforward evaluation metrics, employing both supervised and unsupervised machine learning. The workflow offers a "DIY" aspect, where the performance of any normalization strategy can be evaluated for any type of MS data. As a demonstration of its utility, we apply this workflow on two distinct datasets, an ESI-MS dataset of extracted lipids from latent fingerprints and a cancer spheroid dataset of metabolites ionized by MALDI-MSI, for which we identified the best-performing normalization strategies.

3D reconstruction of murine mitochondria reveals changes in structure during aging linked to the MICOS complex

During aging, muscle gradually undergoes sarcopenia, the loss of function associated with loss of mass, strength, endurance, and oxidative capacity. However, the 3D structural alterations of mitochondria associated with aging in skeletal muscle and cardiac tissues are not well described. Although mitochondrial aging is associated with decreased mitochondrial capacity, the genes responsible for the morphological changes in mitochondria during aging are poorly characterized. We measured changes in mitochondrial morphology in aged murine gastrocnemius, soleus, and cardiac tissues using serial block-face scanning electron microscopy and 3D reconstructions. We also used reverse transcriptase-quantitative PCR, transmission electron microscopy quantification, Seahorse analysis, and metabolomics and lipidomics to measure changes in mitochondrial morphology and function after loss of mitochondria contact site and cristae organizing system (MICOS) complex genes, Chchd3, Chchd6, and Mitofilin. We identified significant changes in mitochondrial size in aged murine gastrocnemius, soleus, and cardiac tissues. We found that both age-related loss of the MICOS complex and knockouts of MICOS genes in mice altered mitochondrial morphology. Given the critical role of mitochondria in maintaining cellular metabolism, we characterized the metabolomes and lipidomes of young and aged mouse tissues, which showed profound alterations consistent with changes in membrane integrity, supporting our observations of age-related changes in muscle tissues. We found a relationship between changes in the MICOS complex and aging. Thus, it is important to understand the mechanisms that underlie the tissue-dependent 3D mitochondrial phenotypic changes that occur in aging and the evolutionary conservation of these mechanisms between Drosophila and mammals.

A benchmark study of data normalisation methods for PTR-TOF-MS exhaled breath metabolomics

Volatilomics is the branch of metabolomics dedicated to the analysis of volatile organic compounds in exhaled breath for medical diagnostic or therapeutic monitoring purposes. Real-time mass spectrometry (MS) technologies such as proton transfer reaction (PTR) MS are commonly used, and data normalisation is an important step to discard unwanted variation from non-biological sources, as batch effects and loss of sensitivity over time may be observed. As normalisation methods for real-time breath analysis have been poorly investigated, we aimed to benchmark known metabolomic data normalisation methods and apply them to PTR-MS data analysis. We compared seven normalisation methods, five statistically based and two using multiple standard metabolites, on two datasets from clinical trials for COVID-19 diagnosis in patients from the emergency department or intensive care unit. We evaluated different means of feature selection to select the standard metabolites, as well as the use of multiple repeat measurements of ambient air to train the normalisation methods. We show that the normalisation tools can correct for time-dependent drift. The methods that provided the best corrections for both cohorts were probabilistic quotient normalisation and normalisation using optimal selection of multiple internal standards. Normalisation also improved the diagnostic performance of the machine learning models, significantly increasing sensitivity, specificity and area under the receiver operating characteristic (ROC) curve for the diagnosis of COVID-19. Our results highlight the importance of adding an appropriate normalisation step during the processing of PTR-MS data, which allows significant improvements in the predictive performance of statistical models.Clinical trials: VOC-COVID-Diag (EudraCT 2020-A02682-37); RECORDS trial (EudraCT 2020-000296-21).

Ablation of Sam50 is associated with fragmentation and alterations in metabolism in murine and human myotubes

The Sorting and Assembly Machinery (SAM) Complex is responsible for assembling β-barrel proteins in the mitochondrial membrane. Comprising three subunits, Sam35, Sam37, and Sam50, the SAM complex connects the inner and outer mitochondrial membranes by interacting with the mitochondrial contact site and cristae organizing system (MICOS) complex. Sam50, in particular, stabilizes the mitochondrial intermembrane space bridging (MIB) complex, which is crucial for protein transport, respiratory chain complex assembly, and regulation of cristae integrity. While the role of Sam50 in mitochondrial structure and metabolism in skeletal muscle remains unclear, this study aims to investigate its impact. Serial block-face-scanning electron microscopy (SBF-SEM) and computer-assisted 3D renderings were employed to compare mitochondrial structure and networking in Sam50-deficient myotubes from mice and humans with wild-type (WT) myotubes. Furthermore, autophagosome 3D structure was assessed in human myotubes. Mitochondrial metabolic phenotypes were assessed using Gas Chromatography-Mass Spectrometry-based metabolomics to explore differential changes in WT and Sam50-deficient myotubes. The results revealed increased mitochondrial fragmentation and autophagosome formation in Sam50-deficient myotubes compared to controls. Metabolomic analysis indicated elevated metabolism of propanoate and several amino acids, including ß-Alanine, phenylalanine, and tyrosine, along with increased amino acid and fatty acid metabolism in Sam50-deficient myotubes. Furthermore, impairment of oxidative capacity was observed upon Sam50 ablation in both murine and human myotubes, as measured with the XF24 Seahorse Analyzer. Collectively, these findings support the critical role of Sam50 in establishing and maintaining mitochondrial integrity, cristae structure, and mitochondrial metabolism. By elucidating the impact of Sam50-deficiency, this study enhances our understanding of mitochondrial function in skeletal muscle.

A critical review of machine-learning for "multi-omics" marine metabolite datasets

During the last decade, genomic, transcriptomic, proteomic, metabolomic, and other omics datasets have been generated for a wide range of marine organisms, and even more are still on the way. Marine organisms possess unique and diverse biosynthetic pathways contributing to the synthesis of novel secondary metabolites with significant bioactivities. As marine organisms have a greater tendency to adapt to stressed environmental conditions, the chance to identify novel bioactive metabolites with potential biotechnological application is very high. This review presents a comprehensive overview of the available "-omics" and "multi-omics" approaches employed for characterizing marine metabolites along with novel data integration tools. The need for the development of machine-learning algorithms for "multi-omics" approaches is briefly discussed. In addition, the challenges involved in the analysis of "multi-omics" data and recommendations for conducting "multi-omics" study were discussed.

Exploration and optimization of extraction, analysis and data normalization strategies for mass spectrometry-based DNA adductome mapping and modeling

Although interest in characterizing DNA damage by means of DNA adductomics has substantially grown, the field of DNA adductomics is still in its infancy, with room for optimization of methods for sample analysis, data processing and DNA adduct identification. In this context, the first objective of this study was to evaluate the use of hydrophilic interaction (HILIC) vs. reversed phase liquid chromatography (RPLC) coupled to high resolution mass spectrometry (HRMS) and thermal acidic vs. enzymatic hydrolysis of DNA followed by DNA adduct purification and enrichment using solid-phase extraction (SPE) or fraction collection for DNA adductome mapping. The second objective was to assess the use of total ion count (TIC) and median intensity (MedI) normalization compared to QC (quality control), iQC (internal QC) and quality control-based robust locally estimated scatterplot smoothing (LOESS) signal correction (QC-RLSC) normalization for processing of the acquired data. The results demonstrate that HILIC compared to RPLC allowed better modeling of the tentative DNA adductome, particularly in combination with thermal acidic hydrolysis and SPE (more valid models, with an average Q²(Y) and R²(Y) of 0.930 and 0.998, respectively). Regarding the need for data normalization and the management of (limited) system instability and signal drift, QC normalization outperformed TIC, MedI, iQC and LOESS normalization. As such, QC normalization can be put forward as the default data normalization strategy. In case of momentous signal drift and/or batch effects however, comparison to other normalization strategies (like e.g. LOESS) is recommended. In future work, further optimization of DNA adductomics may be achieved by merging of HILIC and RPLC datasets and/or application of 2D-LC, as well as the inclusion of Schiff base stabilization and/or fraction collection in the thermal acidic hydrolysis-SPE sample preparation workflow.

Deficiency of endothelial sirtuin1 in mice stimulates skeletal muscle insulin sensitivity by modifying the secretome

Downregulation of endothelial Sirtuin1 (Sirt1) in insulin resistant states contributes to vascular dysfunction. Furthermore, Sirt1 deficiency in skeletal myocytes promotes insulin resistance. Here, we show that deletion of endothelial Sirt1, while impairing endothelial function, paradoxically improves skeletal muscle insulin sensitivity. Compared to wild-type mice, male mice lacking endothelial Sirt1 (E-Sirt1-KO) preferentially utilize glucose over fat, and have higher insulin sensitivity, glucose uptake, and Akt signaling in fast-twitch skeletal muscle. Enhanced insulin sensitivity of E-Sirt1-KO mice is transferrable to wild-type mice via the systemic circulation. Endothelial Sirt1 deficiency, by inhibiting autophagy and activating nuclear factor-kappa B signaling, augments expression and secretion of thymosin beta-4 (Tβ4) that promotes insulin signaling in skeletal myotubes. Thus, unlike in skeletal myocytes, Sirt1 deficiency in the endothelium promotes glucose homeostasis by stimulating skeletal muscle insulin sensitivity through a blood-borne mechanism, and augmented secretion of Tβ4 by Sirt1-deficient endothelial cells boosts insulin signaling in skeletal muscle cells.

Rm-LR: A long-range-based deep learning model for predicting multiple types of RNA modifications

Recent research has highlighted the pivotal role of RNA post-transcriptional modifications in the regulation of RNA expression and function. Accurate identification of RNA modification sites is important for understanding RNA function. In this study, we propose a novel RNA modification prediction method, namely Rm-LR, which leverages a long-range-based deep learning approach to accurately predict multiple types of RNA modifications using RNA sequences only. Rm-LR incorporates two large-scale RNA language pre-trained models to capture discriminative sequential information and learn local important features, which are subsequently integrated through a bilinear attention network. Rm-LR supports a total of ten RNA modification types (m6A, m1A, m5C, m5U, m6Am, Ψ, Am, Cm, Gm, and Um) and significantly outperforms the state-of-the-art methods in terms of predictive capability on benchmark datasets. Experimental results show the effectiveness and superiority of Rm-LR in prediction of various RNA modifications, demonstrating the strong adaptability and robustness of our proposed model. We demonstrate that RNA language pretrained models enable to learn dense biological sequential representations from large-scale long-range RNA corpus, and meanwhile enhance the interpretability of the models. This work contributes to the development of accurate and reliable computational models for RNA modification prediction, providing insights into the complex landscape of RNA modifications.

malbacR: A Package for Standardized Implementation of Batch Correction Methods for Omics Data

Mass spectrometry is a powerful tool for identifying and analyzing biomolecules such as metabolites and lipids in complex biological samples. Liquid chromatography and gas chromatography mass spectrometry studies quite commonly involve large numbers of samples, which can require significant time for sample preparation and analyses. To accommodate such studies, the samples are commonly split into batches. Inevitably, variations in sample handling, temperature fluctuation, imprecise timing, column degradation, and other factors result in systematic errors or biases of the measured abundances between the batches. Numerous methods are available via R packages to assist with batch correction for omics data; however, since these methods were developed by different research teams, the algorithms are available in separate R packages, each with different data input and output formats. We introduce the malbacR package, which consolidates 11 common batch effect correction methods for omics data into one place so users can easily implement and compare the following: pareto scaling, power scaling, range scaling, ComBat, EigenMS, NOMIS, RUV-random, QC-RLSC, WaveICA2.0, TIGER, and SERRF. The malbacR package standardizes data input and output formats across these batch correction methods. The package works in conjunction with the pmartR package, allowing users to seamlessly include the batch effect correction in a pmartR workflow without needing any additional data manipulation.

Denoising Autoencoder Normalization for Large-Scale Untargeted Metabolomics by Gas Chromatography-Mass Spectrometry

Large-scale metabolomics assays are widely used in epidemiology for biomarker discovery and risk assessments. However, systematic errors introduced by instrumental signal drifting pose a big challenge in large-scale assays, especially for derivatization-based gas chromatography-mass spectrometry (GC-MS). Here, we compare the results of different normalization methods for a study with more than 4000 human plasma samples involved in a type 2 diabetes cohort study, in addition to 413 pooled quality control (QC) samples, 413 commercial pooled plasma samples, and a set of 25 stable isotope-labeled internal standards used for every sample. Data acquisition was conducted across 1.2 years, including seven column changes. In total, 413 pooled QC (training) and 413 BioIVT samples (validation) were used for normalization comparisons. Surprisingly, neither internal standards nor sum-based normalizations yielded median precision of less than 30% across all 563 metabolite annotations. While the machine-learning-based SERRF algorithm gave 19% median precision based on the pooled quality control samples, external cross-validation with BioIVT plasma pools yielded a median 34% relative standard deviation (RSD). We developed a new method: systematic error reduction by denoising autoencoder (SERDA). SERDA lowered the median standard deviations of the training QC samples down to 16% RSD, yielding an overall error of 19% RSD when applied to the independent BioIVT validation QC samples. This is the largest study on GC-MS metabolomics ever reported, demonstrating that technical errors can be normalized and handled effectively for this assay. SERDA was further validated on two additional large-scale GC-MS-based human plasma metabolomics studies, confirming the superior performance of SERDA over SERRF or sum normalizations.

The mitochondrial pyruvate carrier complex potentiates the efficacy of proteasome inhibitors in multiple myeloma

Multiple myeloma (MM) is a hematological malignancy that emerges from antibody-producing plasma B cells. Proteasome inhibitors, including the US Food and Drug Administration-approved bortezomib (BTZ) and carfilzomib (CFZ), are frequently used for the treatment of patients with MM. Nevertheless, a significant proportion of patients with MM are refractory or develop resistance to this class of inhibitors, which represents a significant challenge in the clinic. Thus, identifying factors that determine the potency of proteasome inhibitors in MM is of paramount importance to bolster their efficacy in the clinic. Using genome-wide CRISPR-based screening, we identified a subunit of the mitochondrial pyruvate carrier (MPC) complex, MPC1, as a common modulator of BTZ response in 2 distinct human MM cell lines in vitro. We noticed that CRISPR-mediated deletion or pharmacological inhibition of the MPC complex enhanced BTZ/CFZ-induced MM cell death with minimal impact on cell cycle progression. In fact, targeting the MPC complex compromised the bioenergetic capacity of MM cells, which is accompanied by reduced proteasomal activity, thereby exacerbating BTZ-induced cytotoxicity in vitro. Importantly, we observed that the RNA expression levels of several regulators of pyruvate metabolism were altered in advanced stages of MM for which they correlated with poor patient prognosis. Collectively, this study highlights the importance of the MPC complex for the survival of MM cells and their responses to proteasome inhibitors. These findings establish mitochondrial pyruvate metabolism as a potential target for the treatment of MM and an unappreciated strategy to increase the efficacy of proteasome inhibitors in the clinic.

3-hydroxykynurenine is a ROS-inducing cytotoxic tryptophan metabolite that disrupts the TCA cycle

Tryptophan is an essential amino acid that is extensively characterized as a regulator of cellular function through its metabolism by indoleamine 2,3-deoxygenase (IDO) into the kynurenine pathway. However, despite decades of research on tryptophan metabolism, the metabolic regulatory roles of it and its metabolites are not well understood. To address this, we performed an activity metabolomics screen of tryptophan and most of its known metabolites in cell culture. We discovered that treatment of human colon cancer cells (HCT116) with 3-hydroxykynurenine (3-HK), a metabolite of kynurenine, potently disrupted TCA cycle function. Citrate and aconitate levels were increased, while isocitrate and all downstream TCA metabolites were decreased, suggesting decreased aconitase function. We hypothesized that 3HK or one of its metabolites increased reactive oxygen species (ROS) and inhibited aconitase activity. Accordingly, we observed almost complete depletion of reduced glutathione and a decrease in total glutathione levels. We observed a dose-dependent decrease in cell viability after 48 hours of 3HK treatment. These data suggest that raising the intracellular levels of 3HK could be sufficient to induce ROS-mediated apoptosis. We modulated the intracellular levels of 3HK by combined induction of IDO and knockdown of kynureninase (KYNU) in HCT116 cells. Cell viability decreased significantly after 48 hours of KYNU knockdown compared to controls, which was accompanied by increased ROS production and Annexin V staining revealing apoptosis. Finally, we identify xanthommatin production from 3-HK as a candidate radical-producing, cytotoxic mechanism. Our work indicates that KYNU may be a target for disrupting tryptophan metabolism. Interestingly, many cancers exhibit overexpression of IDO, providing a cancer-specific metabolic vulnerability that could be exploited by KYNU inhibition.

MetaPro: a web-based metabolomics application for LC-MS data batch inspection and library curation

INTRODUCTION

Metabolomics analysis based on liquid chromatography-mass spectrometry (LC-MS) has been a prevalent method in the metabolic field. However, accurately quantifying all the metabolites in large metabolomics sample cohorts is challenging. The analysis efficiency is restricted by the abilities of software in many labs, and the lack of spectra for some metabolites also hinders metabolite identification.

OBJECTIVES

Develop software that performs semi-targeted metabolomics analysis with an optimized workflow to improve quantification accuracy. The software also supports web-based technologies and increases laboratory analysis efficiency. A spectral curation function is provided to promote the prosperity of homemade MS/MS spectral libraries in the metabolomics community.

METHODS

MetaPro is developed based on an industrial-grade web framework and a computation-oriented MS data format to improve analysis efficiency. Algorithms from mainstream metabolomics software are integrated and optimized for more accurate quantification results. A semi-targeted analysis workflow is designed based on the concept of combining artificial judgment and algorithm inference.

RESULTS

MetaPro supports semi-targeted analysis workflow and functions for fast QC inspection and self-made spectral library curation with easy-to-use interfaces. With curated authentic or high-quality spectra, it can improve identification accuracy using different peak identification strategies. It demonstrates practical value in analyzing large amounts of metabolomics samples.

CONCLUSION

We offer MetaPro as a web-based application characterized by fast batch QC inspection and credible spectral curation towards high-throughput metabolomics data. It aims to resolve the analysis difficulty in semi-targeted metabolomics.

Instrumental Drift in Untargeted Metabolomics: Optimizing Data Quality with Intrastudy QC Samples

Untargeted metabolomics is an important tool in studying health and disease and is employed in fields such as biomarker discovery and drug development, as well as precision medicine. Although significant technical advances were made in the field of mass-spectrometry driven metabolomics, instrumental drifts, such as fluctuations in retention time and signal intensity, remain a challenge, particularly in large untargeted metabolomics studies. Therefore, it is crucial to consider these variations during data processing to ensure high-quality data. Here, we will provide recommendations for an optimal data processing workflow using intrastudy quality control (QC) samples that identifies errors resulting from instrumental drifts, such as shifts in retention time and metabolite intensities. Furthermore, we provide an in-depth comparison of the performance of three popular batch-effect correction methods of different complexity. By using different evaluation metrics based on QC samples and a machine learning approach based on biological samples, the performance of the batch-effect correction methods were evaluated. Here, the method TIGER demonstrated the overall best performance by reducing the relative standard deviation of the QCs and dispersion-ratio the most, as well as demonstrating the highest area under the receiver operating characteristic with three different probabilistic classifiers (Logistic regression, Random Forest, and Support Vector Machine). In summary, our recommendations will help to generate high-quality data that are suitable for further downstream processing, leading to more accurate and meaningful insights into the underlying biological processes.

Integrative computational modeling to unravel novel potential biomarkers in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a major health problem around the world. The management of this disease is complicated by the lack of noninvasive diagnostic tools and the few treatment options available. Better clinical outcomes can be achieved if HCC is detected early, but unfortunately, clinical signs appear when the disease is in its late stages. We aim to identify novel genes that can be targeted for the diagnosis and therapy of HCC. We performed a meta-analysis of transcriptomics data to identify differentially expressed genes and applied network analysis to identify hub genes. Fatty acid metabolism, complement and coagulation cascade, chemical carcinogenesis and retinol metabolism were identified as key pathways in HCC. Furthermore, we integrated transcriptomics data into a reference human genome-scale metabolic model to identify key reactions and subsystems relevant in HCC. We conclude that fatty acid activation, purine metabolism, vitamin D, and E metabolism are key processes in the development of HCC and therefore need to be further explored for the development of new therapies. We provide the first evidence that GABRP, HBG1 and DAK (TKFC) genes are important in HCC in humans and warrant further studies.

Clustering ensemble in scRNA-seq data analysis: Methods, applications and challenges

With the rapid development of single-cell RNA-sequencing techniques, various computational methods and tools were proposed to analyze these high-throughput data, which led to an accelerated reveal of potential biological information. As one of the core steps of single-cell transcriptome data analysis, clustering plays a crucial role in identifying cell types and interpreting cellular heterogeneity. However, the results generated by different clustering methods showed distinguishing, and those unstable partitions can affect the accuracy of the analysis to a certain extent. To overcome this challenge and obtain more accurate results, currently clustering ensemble is frequently applied to cluster analysis of single-cell transcriptome datasets, and the results generated by all clustering ensembles are nearly more reliable than those from most of the single clustering partitions. In this review, we summarize applications and challenges of the clustering ensemble method in single-cell transcriptome data analysis, and provide constructive thoughts and references for researchers in this field.

Methodology to identify a gene expression signature by merging microarray datasets

A vast number of microarray datasets have been produced as a way to identify differentially expressed genes and gene expression signatures. A better understanding of these biological processes can help in the diagnosis and prognosis of diseases, as well as in the therapeutic response to drugs. However, most of the available datasets are composed of a reduced number of samples, leading to low statistical, predictive and generalization power. One way to overcome this problem is by merging several microarray datasets into a single dataset, which is typically a challenging task. Statistical methods or supervised machine learning algorithms are usually used to determine gene expression signatures. Nevertheless, statistical methods require an arbitrary threshold to be defined, and supervised machine learning methods can be ineffective when applied to high-dimensional datasets like microarrays. We propose a methodology to identify gene expression signatures by merging microarray datasets. This methodology uses statistical methods to obtain several sets of differentially expressed genes and uses supervised machine learning algorithms to select the gene expression signature. This methodology was validated using two distinct research applications: one using heart failure and the other using autism spectrum disorder microarray datasets. For the first, we obtained a gene expression signature composed of 117 genes, with a classification accuracy of approximately 98%. For the second use case, we obtained a gene expression signature composed of 79 genes, with a classification accuracy of approximately 82%. This methodology was implemented in R language and is available, under the MIT licence, at https://github.com/bioinformatics-ua/MicroGES.

Normalization of organ-on-a-Chip samples for mass spectrometry based proteomics and metabolomics via Dansylation-based assay

Mass spectrometry based 'omics pairs well with organ-on-a-chip-based investigations, which often have limited cellular material for sampling. However, a common issue with these chip-based platforms is well-to-well or chip-to-chip variability in the proteome and metabolome due to factors such as plate edge effects, cellular asynchronization, effluent flow, and limited cell count. This causes high variability in the quantitative multi-omics analysis of samples, potentially masking true biological changes within the system. Solutions to this have been approached via data processing tools and post-acquisition normalization strategies such as constant median, constant sum, and overall signal normalization. Unfortunately, these methods do not adequately correct for the large variations, resulting in a need for increased biological replicates. The methods in this work utilize a dansylation based assay with a subset of labeled metabolites that allow for pre-acquisition normalization to better correlate the biological perturbations that truly occur in chip-based platforms. BCA protein assays were performed in tandem with a proteomics pipeline to achieve pre-acquisition normalization. The CN Bio PhysioMimix was seeded with primary hepatocytes and challenged with VX after six days of culture, and the metabolome and proteome were analyzed using the described normalization methods. A decreased coefficient of variation percentage is achieved, significant changes are observed through the proteome and metabolome, and better classification of biological replicates acquired because of these strategies.

Deep learning automates bidimensional and volumetric tumor burden measurement from MRI in pre- and post-operative glioblastoma patients

Tumor burden assessment by magnetic resonance imaging (MRI) is central to the evaluation of treatment response for glioblastoma. This assessment is, however, complex to perform and associated with high variability due to the high heterogeneity and complexity of the disease. In this work, we tackle this issue and propose a deep learning pipeline for the fully automated end-to-end analysis of glioblastoma patients. Our approach simultaneously identifies tumor sub-regions, including the enhancing tumor, peritumoral edema and surgical cavity in the first step, and then calculates the volumetric and bidimensional measurements that follow the current Response Assessment in Neuro-Oncology (RANO) criteria. Also, we introduce a rigorous manual annotation process which was followed to delineate the tumor sub-regions by the human experts, and to capture their segmentation confidences that are later used while training deep learning models. The results of our extensive experimental study performed over 760 pre-operative and 504 post-operative adult patients with glioma obtained from the public database (acquired at 19 sites in years 2021-2020) and from a clinical treatment trial (47 and 69 sites for pre-/post-operative patients, 2009-2011) and backed up with thorough quantitative, qualitative and statistical analysis revealed that our pipeline performs accurate segmentation of pre- and post-operative MRIs in a fraction of the manual delineation time (up to 20 times faster than humans). Volumetric measurements were in strong agreement with experts with the Intraclass Correlation Coefficient (ICC): 0.959, 0.703, 0.960 for ET, ED, and cavity. Similarly, automated RANO compared favorably with experienced readers (ICC: 0.681 and 0.866) producing consistent and accurate results. Additionally, we showed that RANO measurements are not always sufficient to quantify tumor burden. The high performance of the automated tumor burden measurement highlights the potential of the tool for considerably improving and simplifying radiological evaluation of glioblastoma in clinical trials and clinical practice.

scMAGS: Marker gene selection from scRNA-seq data for spatial transcriptomics studies

Single-Cell RNA sequencing (scRNA-seq) has provided unprecedented opportunities for exploring gene expression and thus uncovering regulatory relationships between genes at the single-cell level. However, scRNA-seq relies on isolating cells from tissues. Therefore, the spatial context of the regulatory processes is lost. A recent technological innovation, spatial transcriptomics, allows for the measurement of gene expression while preserving spatial information. An initial step in the spatial transcriptomic analysis is to identify the cell type, which requires a careful selection of cell-specific marker genes. For this purpose, currently, scRNA-seq data is used to select a limited number of marker genes from among all genes that distinguish cell types from each other. This study proposes scMAGS (single-cell MArker Gene Selection), a novel method for marker gene selection from scRNA-seq data for spatial transcriptomics studies. scMAGS uses a filtering step in which the candidate genes are identified before the marker gene selection step. For the selection of marker genes, cluster validity indices, the Silhouette index, or the Calinski-Harabasz index (for large datasets) are utilized. Experimental results showed that, in comparison to the existing methods, scMAGS is scalable, fast, and accurate. Even for large datasets with millions of cells, scMAGS could find the required number of marker genes in a reasonable amount of time with fewer memory requirements. scMAGS is made freely available at https://github.com/doganlab/scmags and can be downloaded from the Python Package Directory (PyPI) software repository with the command pip install scmags.