A New Strategy for Analyzing Time-Series Data Using Dynamic Networks: Identifying Prospective Biomarkers of Hepatocellular Carcinoma

Harnessing the Power of AI for Enhanced Diagnosis and Treatment of Hepatocellular Carcinoma

Since its advent, artificial intelligence (AI) has been continuously researched, and substantial progress has been made in many fields, such as the diagnosis and therapies for cancer. Due to the advantages of high efficiency, rapidity, and precision, AI has been increasingly adopted in the medical field to improve patient prognosis and the efficiency of medical procedures. Thus, AI technology has become a powerful driving force and support mechanism in the medical and health industry. Hepatocellular carcinoma (HCC) is the sixth most common cancer and the fourth leading cause of cancer mortality worldwide, with a 5-year relative survival rate of 18%. The early diagnosis and postsurgical treatment of HCC are key factors related to its prognosis, which provides an opportunity for the application of AI in HCC diagnosis and treatment. This review will summarize the application of AI in the diagnosis and treatment of HCC to provide prospects for deeper and wider applications.

Analyzing omics data based on sample network

Identifying valuable features from complex omics data is of great significance for disease diagnosis study. This paper proposes a new feature selection algorithm based on sample network (FS-SN) to mine important information from omics data. The sample network is constructed according to the sample neighbor relationship at the molecular (feature) expression level, and the distinguishing ability of the feature is evaluated based on the topology of the sample network. The sample network established on a feature with a strong discriminating ability tends to have many edges between the same group samples and few edges between the different group samples. At the same time, FS-SN removes redundant features according to the gravitational interaction between features. To show the validation of FS-SN, it was compared on ten public datasets with ERGS, mRMR, ReliefF, ATSD-DN, and INDEED which are efficient in omics data analysis. Experimental results show that FS-SN performed better than the compared methods in accuracy, sensitivity and specificity in most cases. Hence, FS-SN making use of the topology of the sample network is effective for analyzing omics data, it can identify key features that reflect the occurrence and development of diseases, and reveal the underlying biological mechanism.

Multiomics characterization of fatty acid metabolism for the clinical management of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a prevalent malignancy and there is a lack of effective biomarkers for HCC diagnosis. Living organisms are complex, and different omics molecules interact with each other to implement various biological functions. Genomics and metabolomics, which are the top and bottom of systems biology, play an important role in HCC clinical management. Fatty acid metabolism is associated with malignancy, prognosis, and immune phenotype in cancer, which is a potential hallmark in malignant tumors. In this study, the genes and metabolites related to fatty acid metabolism were thoroughly investigated by a dynamic network construction algorithm named EWS-DDA for the early diagnosis and prognosis of HCC. Three gene ratios and eight metabolite ratios were identified by EWS-DDA as potential biomarkers for HCC clinical management. Further analysis using biological analysis, statistical analysis and document validation in the discovery and validation sets suggested that the selected potential biomarkers had great clinical prognostic value and helped to achieve effective early diagnosis of HCC. Experimental results suggested that in-depth evaluation of fatty acid metabolism from different omics viewpoints can facilitate the further understanding of pathological alterations associated with HCC characteristics, improving the performance of early diagnosis and clinical prognosis.

Gut microbiota and holobiont metabolome composition of the medaka fish (Oryzias latipes) are affected by a short exposure to the cyanobacterium Microcystis aeruginosa

Blooms of toxic cyanobacteria are a common stress encountered by aquatic fauna. Evidence indicates that long-lasting blooms affect fauna-associated microbiota. Because of their multiple roles, host-associated microbes are nowadays considered relevant to ecotoxicology, yet the respective timing of microbiota versus functional changes in holobionts response needs to be clarified. The response of gut microbiota and holobiont's metabolome to exposure to a dense culture of Microcystis aeruginosa was investigated as a microcosm-simulated bloom in the model fish species Oryzias latipes (medaka). Both gut microbiota and gut metabolome displayed significant composition changes after only 2 days of exposure. A dominant symbiont, member of the Firmicutes, plummeted whereas various genera of Proteobacteria and Actinobacteriota increased in relative abundance. Changes in microbiota composition occurred earlier and faster compared to metabolome composition. Liver and muscle metabolome were much less affected than guts, supporting that the gut and associated microbiota are in the front row upon exposure. This study highlights that even short cyanobacterial blooms, that are increasingly frequent, trigger changes in microbiota composition and holobiont metabolome. It emphasizes the relevance of multi-omics approaches to explore organism's response to an ecotoxicological stress.

A Time-Series Metabolomic Analysis of SARS-CoV-2 Infection in a Ferret Model

The global threat of COVID-19 has led to an increased use of metabolomics to study SARS-CoV-2 infections in animals and humans. In spite of these efforts, however, understanding the metabolome of SARS-CoV-2 during an infection remains difficult and incomplete. In this study, metabolic responses to a SAS-CoV-2 challenge experiment were studied in nasal washes collected from an asymptomatic ferret model (n = 20) at different time points before and after infection using an LC-MS-based metabolomics approach. A multivariate analysis of the nasal wash metabolome data revealed several statistically significant features. Despite no effects of sex or interaction between sex and time on the time course of SARS-CoV-2 infection, 16 metabolites were significantly different at all time points post-infection. Among these altered metabolites, the relative abundance of taurine was elevated post-infection, which could be an indication of hepatotoxicity, while the accumulation of sialic acids could indicate SARS-CoV-2 invasion. Enrichment analysis identified several pathways influenced by SARS-CoV-2 infection. Of these, sugar, glycan, and amino acid metabolisms were the key altered pathways in the upper respiratory channel during infection. These findings provide some new insights into the progression of SARS-CoV-2 infection in ferrets at the metabolic level, which could be useful for the development of early clinical diagnosis tools and new or repurposed drug therapies.

Identification of Endoplasmic Reticulum Stress-Related Biomarkers of Periodontitis Based on Machine Learning: A Bioinformatics Analysis

Objective

To screen for potential endoplasmic reticulum stress- (ERS-) related biomarkers of periodontitis using machine learning methods and explore their relationship with immune cells.

Methods

Three datasets of periodontitis (GSE10334, GES16134, and GES23586) were obtained from the Gene Expression Omnibus (GEO), and the samples were randomly assigned to the training set or the validation set. ERS-related differentially expressed genes (DEGs) between periodontitis and healthy periodontal tissues were screened and analyzed for GO, KEGG, and DO enrichment. Key DEGs were screened by two machine learning algorithms, LASSO regression and support vector machine-recursive feature elimination (SVM-RFE); then, the potential biomarkers were identified through validation. The infiltration of immune cells of periodontitis was calculated using the CIBERSORT algorithm, and the correlation between immune cells and potential biomarkers was specifically analyzed through the Spearman method.

Results

We obtained 36 ERS-related DEGs of periodontitis from the training set, from which 11 key DEGs were screened by further machine learning. SERPINA1, ERLEC1, and VWF showed high diagnostic values (AUC > 0.85), so they were considered as potential biomarkers for periodontitis. According to the results of the immune cell infiltration analysis, these three potential biomarkers showed marked correlations with plasma cells, neutrophils, resting dendritic cells, resting mast cells, and follicular helper T cells.

Conclusions

Three ERS-related genes, SERPINA1, ERLEC1, and VWF, showed valuable biomarker potential for periodontitis, which provide a target base for future studies on early diagnosis and treatment of periodontitis.

A Novel Method for Constructing Classification Models by Combining Different Biomarker Patterns

Different biomarker patterns, such as those of molecular biomarkers and ratio biomarkers, have their own merits in clinical applications. In this study, a novel machine learning method used in biomedical data analysis for constructing classification models by combining different biomarker patterns (CDBP)is proposed. CDBP uses relative expression reversals to measure the discriminative ability of different biomarker patterns, and selects the pattern with the higher score for classifier construction. The decision boundary of CDBP can be characterized in simple and biologically meaningful manners. The CDBP method was compared with eight state-of-the-art methods on eight gene expression datasets to test its performance. CDBP, with fewer features or ratio features, had the highest classification performance. Subsequently, CDBP was employed to extract crucial diagnostic information from a rat hepatocarcinogenesis metabolomics dataset. The potential biomarkers selected by CDBP provided better classification of hepatocellular carcinoma (HCC)and non-HCC stages than previous works in the animal model. The statistical analyses of these potential biomarkers in an independent human dataset confirmed their discriminative abilities of different liver diseases. These experimental results highlight the potential of CDBP for biomarker identification from high-dimensional biomedical datasets and demonstrate that it can be a useful tool for disease classification.

Data analysis methods for defining biomarkers from omics data

Omics mainly includes genomics, epigenomics, transcriptomics, proteomics and metabolomics. The rapid development of omics technology has opened up new ways to study disease diagnosis and prognosis and to define prospective information of complex diseases. Since omics data are usually large and complex, the method used to analyze the data and to define important information is crucial in omics study. In this review, we focus on advances in biomarker discovery methods based on omics data in the last decade, and categorize them as individual feature analysis, combinatorial feature analysis and network analysis. We also discuss the challenges and perspectives in this field.

A computational strategy for metabolic network construction based on the overlapping ratio: Study of patients' metabolic responses to different dialysis patterns

BACKGROUND

Uremia is a worldwide epidemic disease and poses a serious threat to human health. Both maintenance hemodialysis (HD) and maintenance high flux hemodialysis (HFD) are common treatments for uremia and are generally used in clinical applications. In-depth exploration of patients' metabolic responses to different dialysis patterns can facilitate the understanding of pathological alterations associated with uremia and the effects of different dialysis methods on uremia, which may be used for future personalized therapy. However, due to variations of multiple factors (i.e., genetic, epigenetic and environment) in the process of disease treatments, identification of the similarities and differences in plasma metabolite changes in uremic patients in response to HD and HFD remains challenging.

METHODS

In this study, a computational strategy for metabolic network construction based on the overlapping ratio (MNC-OR) was proposed for disease treatment effect research. In MNC-OR, the overlapping ratio was introduced to measure metabolic reactions and to construct metabolic networks for analysis of different treatment options. Then, MNC-OR was employed to analyze HD-pattern-dependent changes in plasma metabolites to explore the pathological alterations associated with uremia and the effectiveness of different dialysis patterns (i.e., HD and HFD) on uremia. Based on the networks constructed by MNC-OR, two network analysis techniques, namely, similarity analysis and difference analysis of network topology, were used to find the similarity and differences in metabolic signals in patients under treatment with either HD or HFD, which can facilitate the understanding of pathological alterations associated with uremia and provide the guidance for personalized dialysis therapy.

RESULTS

Similarity analysis of network topology suggested that abnormal energy metabolism, gut metabolism and pyrimidine metabolism might occur in uremic patients, and maintenance of both HFD and HD therapies have beneficial effects on uremia. Then, difference analysis of network topology was employed to extract the crucial information related to HD-pattern-dependent changes in plasma metabolites. Experimental results indicated that the amino acid metabolism was closer to the normal status in HFD-treated patients; however, in HD-treated patients, the ability of antioxidation showed greater reduction, and the protein O-GlcNAcylation level was higher. Our findings demonstrate the potential of MNC-OR for explaining the metabolic similarities and differences of patients in response to different dialysis methods, thereby contributing to the guidance of personalized dialysis therapy.

Bayesian approach for predicting responses to therapy from high-dimensional time-course gene expression profiles

Background

Historical and updated information provided by time-course data collected during an entire treatment period proves to be more useful than information provided by single-point data. Accurate predictions made using time-course data on multiple biomarkers that indicate a patient’s response to therapy contribute positively to the decision-making process associated with designing effective treatment programs for various diseases. Therefore, the development of prediction methods incorporating time-course data on multiple markers is necessary.

Results

We proposed new methods that may be used for prediction and gene selection via time-course gene expression profiles. Our prediction method consolidated multiple probabilities calculated using gene expression profiles collected over a series of time points to predict therapy response. Using two data sets collected from patients with hepatitis C virus (HCV) infection and multiple sclerosis (MS), we performed numerical experiments that predicted response to therapy and evaluated their accuracies. Our methods were more accurate than conventional methods and successfully selected genes, the functions of which were associated with the pathology of HCV infection and MS.

Conclusions

The proposed method accurately predicted response to therapy using data at multiple time points. It showed higher accuracies at early time points compared to those of conventional methods. Furthermore, this method successfully selected genes that were directly associated with diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04052-4.

Emerging computational tools and models for studying gut microbiota composition and function

The gut microbiota and its metabolites play critical roles in human health and disease. Advances in high-throughput sequencing, mass spectrometry, and other omics assay platforms have improved our ability to generate large volumes of data exploring the temporal variations in the compositions and functions of microbial communities. To elucidate mechanisms, methods and tools are needed that can rigorously model the dependencies within time-series data. Longitudinal data are often sparse and unevenly sampled, and nontrivial challenges remain in determining statistical significance, normalization across different data types, and model validation. In this review, we highlight recent developments in models and software tools for the analysis of time series microbiome and metabolome data, as well as integration of these data.

Breath Metabolomics Provides an Accurate and Noninvasive Approach for Screening Cirrhosis, Primary, and Secondary Liver Tumors

Hepatocellular carcinoma (HCC) and secondary liver tumors, such as colorectal cancer liver metastases are significant contributors to the overall burden of cancer-related morality. Current biomarkers, such as alpha-fetoprotein (AFP) for HCC, result in too many false negatives, necessitating noninvasive approaches with improved sensitivity. Volatile organic compounds (VOCs) detected in the breath of patients can provide valuable insight into disease processes and can differentiate patients by disease status. Here, we investigate whether 22 VOCs from the breath of 296 patients can distinguish those with no liver disease (n = 54), cirrhosis (n = 30), HCC (n = 112), pulmonary hypertension (n = 49), or colorectal cancer liver metastases (n = 51). This work extends previous studies by evaluating the ability for VOC signatures to differentiate multiple diseases in a large cohort of patients. Pairwise disease comparisons demonstrated that most of the VOCs tested are present in significantly different relative abundances (false discovery rate P < 0.1), indicating broad impacts on the breath metabolome across diseases. A predictive model developed using random forest machine learning and cross validation classified patients with 85% classification accuracy and 75% balanced accuracy. Importantly, the model detected HCC with 73% sensitivity compared with 53% for AFP in the same cohort. An added value of this approach is that influential VOCs in the predictive model may provide insight into disease etiology. Acetaldehyde and acetone, both of which have roles in tumor promotion, were considered important VOCs for differentiating disease groups in the predictive model and were increased in patients with cirrhosis and HCC compared to patients with no liver disease (false discovery rate P < 0.1). Conclusion: The use of machine learning and breath VOCs shows promise as an approach to develop improved, noninvasive screening tools for chronic liver disease and primary and secondary liver tumors.

Identifying Early Warning Signals for the Sudden Transition from Mild to Severe Tobacco Etch Disease by Dynamical Network Biomarkers

Complex systems exhibit critical thresholds at which they transition among alternative phases. Complex systems theory has been applied to analyze disease progression, distinguishing three stages along progression: (i) a normal noninfected state; (ii) a predisease state, in which the host is infected and responds and therapeutic interventions could still be effective; and (iii) an irreversible state, where the system is seriously threatened. The dynamical network biomarker (DNB) theory sought for early warnings of the transition from health to disease. Such DNBs might range from individual genes to complex structures in transcriptional regulatory or protein-protein interaction networks. Here, we revisit transcriptomic data obtained during infection of tobacco plants with tobacco etch potyvirus to identify DNBs signaling the transition from mild/reversible to severe/irreversible disease. We identified genes showing a sudden transition in expression along disease categories. Some of these genes cluster in modules that show the properties of DNBs. These modules contain both genes known to be involved in response to pathogens (e.g., ADH2, CYP19, ERF1, KAB1, LAP1, MBF1C, MYB58, PR1, or TPS5) and other genes not previously related to biotic stress responses (e.g., ABCI6, BBX21, NAP1, OSM34, or ZPN1).

Juice Index: an integrated Sauvignon blanc grape and wine metabolomics database shows mainly seasonal differences

INTRODUCTION

Although Sauvignon Blanc (SB) grapes are cultivated widely throughout New Zealand, wines from the Marlborough region are most famous for their typical varietal combination of tropical and vegetal aromas. These wines differ in composition from season to season as well as among locations within the region, which makes the continual production of good quality wines challenging. Here, we developed a unique database of New Zealand SB grape juices and wines to develop tools to help winemakers to make blending decisions and assist in the development of new wine styles.

METHODS

About 400 juices were collected from different regions in New Zealand over three harvest seasons (2011-2013), which were then fermented under controlled conditions using a commercial yeast strain Saccharomyces cerevisiae EC1118. Comprehensive metabolite profiling of these juices and wines by gas chromatography-mass spectrometry (GC-MS) was combined with their detailed oenological parameters and associated meteorological data.

RESULTS

These combined metabolomics data clearly demonstrate that seasonal variation is more prominent than regional difference in both SB grape juices and wines, despite almost universal use of vineyard irrigation to mitigate seasonal rainfall and evapotranspiration differences, Additionally, we identified a group of juice metabolites that play central roles behind these variations, which may represent chemical signatures for juice and wine quality assessment.

CONCLUSION

This database is the first of its kind in the world to be available for the wider scientific community and offers potential as a predictive tool for wine quality and innovation when combined with mathematical modelling.

Device-measured physical activity data for classification of patients with ventricular arrhythmia events: A pilot investigation

Low levels of physical activity are associated with increased mortality risk, especially in cardiac patients, but most studies are based on self-report. Cardiac implantable electronic devices (CIEDs) offer an opportunity to collect data for longer periods of time. However, there is limited agreement on the best approaches for quantification of activity measures due to the time series nature of the data. We examined physical activity time series data from 235 subjects with CIEDs and at least 365 days of uninterrupted measures. Summary statistics for raw daily physical activity (minutes/day), including statistical moments (e.g., mean, standard deviation, skewness, kurtosis), time series regression coefficients, frequency domain components, and forecasted predicted values, were calculated for each individual, and used to predict occurrence of ventricular tachycardia (VT) events as recorded by the device. In unsupervised analyses using principal component analysis, we found that while certain features tended to cluster near each other, most provided a reasonable spread across activity space without a large degree of redundancy. In supervised analyses, we found several features that were associated with the outcome (P < 0.05) in univariable and multivariable approaches, but few were consistent across models. Using a machine-learning approach in which the data was split into training and testing sets, and models ranging in complexity from simple univariable logistic regression to ensemble decision trees were fit, there was no improvement in classification of risk over naïve methods for any approach. Although standard approaches identified summary features of physical activity data that were correlated with risk of VT, machine-learning approaches found that none of these features provided an improvement in classification. Future studies are needed to explore and validate methods for feature extraction and machine learning in classification of VT risk based on device-measured activity.

Lipid profiling of pre-treatment plasma reveals biomarker candidates associated with response rates and hand-foot skin reactions in sorafenib-treated patients

Sorafenib is a multi-kinase inhibitor for treatment of advanced hepatocellular carcinoma (HCC). Beyond its clinical benefit against advanced HCC, the efficacy and safety of sorafenib chemotherapy are critical concerns. In this study, we addressed the lipid profiles associated with the efficacy and safety of sorafenib chemotherapy. Plasma samples from HCC patients before sorafenib chemotherapy (N = 44) were collected and subjected to lipidomic analysis. We measured the levels of 176 lipids belonging to 8 classes of phosphoglycerolipids, 2 classes of sphingolipids, 3 classes of neutral lipids, and 4 other classes of lipids. To characterize lipids associated with efficacy, we compared the responder group (N = 21; partial response and stable disease) with non-responder group (N = 22; progressive disease). To characterize lipids associated with hand-foot skin reaction (HFSR), we compared the susceptible group (N = 12; grade 2 and 3) with non-susceptible group (N = 32; grade 0 and 1). The levels of 8 lipids, including phosphatidylcholine (PC)[34:2], PC[34:3]a, PC[35:2], PC[36:4]a, PC[34:3e], acylcarnitine (Car)[18:0], cholesterol ester[20:2], and diacylglycerol (DG)[34:2], were significantly lower in the responder group, and 6 out of 8 these lipids contained FA(18:2). In addition, the levels of 7 lipids (Car[12:0], Car[18:0], Car[18:1], Car[20:1] and fatty acid amides (FAA[16:0], FAA[18:0], and FAA[18:1]b)) were significantly lower in the group susceptible to HFSR. Our comprehensive lipidomics study using samples from sorafenib-treated patients with HCC revealed that significant differences in the lipid profiles of pre-treatment plasma were associated with sorafenib efficacy and sorafenib-induced HFSR. Validation using another set of patient plasma samples and elucidating the molecular basis of these changes will lead to better treatment with sorafenib chemotherapy.

Deciphering hepatocellular carcinoma through metabolomics: from biomarker discovery to therapy evaluation

Hepatocellular carcinoma (HCC) is the third most common cause of death from cancer, with increasing prevalence worldwide. The mortality rate of HCC is similar to its incidence rate, which reflects its poor prognosis. At present, the diagnosis of HCC is still mostly dependent on invasive biopsy, imaging methods, and serum α-fetoprotein (AFP) testing. Because of the asymptomatic nature of early HCC, biopsy and imaging methods usually detect HCC at the middle-late stages. AFP has limited sensitivity and specificity, as many other nonmalignant liver diseases can also result in a very high serum level of AFP. Therefore, better biomarkers with higher sensitivity and specificity at earlier stages are greatly needed. Since metabolic reprogramming is an essential hallmark of cancer and the liver is the metabolic hub of living systems, it is useful to investigate HCC from a metabolic perspective. As a noninvasive and nondestructive approach, metabolomics provides holistic information on dynamically metabolic responses of living systems to both endogenous and exogenous factors. Therefore, it would be conducive to apply metabolomics in investigating HCC. In this review, we summarize recent metabolomic studies on HCC cellular, animal, and clinicopathologic models with attention to metabolomics as a biomarker in cancer diagnosis. Recent applications of metabolomics with respect to therapeutic and prognostic evaluation of HCC are also covered, with emphasis on the potential of treatment by drugs from natural products. In the last section, the current challenges and trends of future development of metabolomics on HCC are discussed. Overall, metabolomics provides us with novel insight into the diagnosis, prognosis, and therapeutic evaluation of HCC.

Elastic net regularized regression for time-series analysis of plasma metabolome stability under sub-optimal freezing condition

In this paper, the stability of the plasma metabolome at −20 °C for up to 30 days was evaluated using liquid chromatography-high resolution mass spectrometric metabolomics analysis. To follow the time-series deterioration of the plasma metabolome, the use of an elastic net regularized regression model for the prediction of storage time at −20 °C based on the plasma metabolomic profile, and the selection and ranking of metabolites with high temporal changes was demonstrated using the glmnet package in R. Out of 1229 (positive mode) and 1483 (negative mode) metabolite features, the elastic net model extracted 32 metabolites of interest in both positive and negative modes. L-gamma-glutamyl-L-(iso)leucine (tentative identification) was found to have the highest time-dependent change and significantly increased proportionally to the storage time of plasma at −20 °C (R² = 0.6378 [positive mode], R² = 0.7893 [negative mode], p-value < 0.00001). Based on the temporal profiles of the extracted metabolites by the model, results show only minimal deterioration of the plasma metabolome at −20 °C up to 1 month. However, majority of the changes appeared at around 12–15 days of storage. This allows scientists to better plan logistics and storage strategies for samples obtained from low-resource settings, where −80 °C storage is not guaranteed.

A Computational Method of Defining Potential Biomarkers based on Differential Sub-Networks

Analyzing omics data from a network-based perspective can facilitate biomarker discovery. To improve disease diagnosis and identify prospective information indicating the onset of complex disease, a computational method for identifying potential biomarkers based on differential sub-networks (PB-DSN) is developed. In PB-DSN, Pearson correlation coefficient (PCC) is used to measure the relationship between feature ratios and to infer potential networks. A differential sub-network is extracted to identify crucial information for discriminating different groups and indicating the emergence of complex diseases. Subsequently, PB-DSN defines potential biomarkers based on the topological analysis of these differential sub-networks. In this study, PB-DSN is applied to handle a static genomics dataset of small, round blue cell tumors and a time-series metabolomics dataset of hepatocellular carcinoma. PB-DSN is compared with support vector machine-recursive feature elimination, multivariate empirical Bayes statistics, analyzing time-series data based on dynamic networks, molecular networks based on PCC, PinnacleZ, graph-based iterative group analysis, KeyPathwayMiner and BioNet. The better performance of PB-DSN not only demonstrates its effectiveness for the identification of discriminative features that facilitate disease classification, but also shows its potential for the identification of warning signals.

Extracellular Microbial Metabolomics: The State of the Art

Microorganisms produce and secrete many primary and secondary metabolites to the surrounding environment during their growth. Therefore, extracellular metabolites provide important information about the changes in microbial metabolism due to different environmental cues. The determination of these metabolites is also comparatively easier than the extraction and analysis of intracellular metabolites as there is no need for cell rupture. Many analytical methods are already available and have been used for the analysis of extracellular metabolites from microorganisms over the last two decades. Here, we review the applications and benefits of extracellular metabolite analysis. We also discuss different sample preparation protocols available in the literature for both types (e.g., metabolites in solution and in gas) of extracellular microbial metabolites. Lastly, we evaluate the authenticity of using extracellular metabolomics data in the metabolic modelling of different industrially important microorganisms.