Bioinformatics tools and challenges in structural analysis of lipidomics MS/MS data

LipidSIM: Inferring mechanistic lipid biosynthesis perturbations from lipidomics with a flexible, low-parameter, Markov modeling framework

Lipid metabolism is a complex and dynamic system involving numerous enzymes at the junction of multiple metabolic pathways. Disruption of these pathways leads to systematic dyslipidemia, a hallmark of many pathological developments, such as nonalcoholic steatohepatitis and diabetes. Recent advances in computational tools can provide insights into the dysregulation of lipid biosynthesis, but limitations remain due to the complexity of lipidomic data, limited knowledge of interactions among involved enzymes, and technical challenges in standardizing across different lipid types. Here, we present a low-parameter, biologically interpretable framework named Lipid Synthesis Investigative Markov model (LipidSIM), which models and predicts the source of perturbations in lipid biosynthesis from lipidomic data. LipidSIM achieves this by accounting for the interdependency between the lipid species via the lipid biosynthesis network and generates testable hypotheses regarding changes in lipid biosynthetic reactions. This feature allows the integration of lipidomics with other omics types, such as transcriptomics, to elucidate the direct driving mechanisms of altered lipidomes due to treatments or disease progression. To demonstrate the value of LipidSIM, we first applied it to hepatic lipidomics following Keap1 knockdown and found that changes in mRNA expression of the lipid pathways were consistent with the LipidSIM-predicted fluxes. Second, we used it to study lipidomic changes following intraperitoneal injection of CCl4 to induce fast NAFLD/NASH development and the progression of fibrosis and hepatic cancer. Finally, to show the power of LipidSIM for classifying samples with dyslipidemia, we used a Dgat2-knockdown study dataset. Thus, we show that as it demands no a priori knowledge of enzyme kinetics, LipidSIM is a valuable and intuitive framework for extracting biological insights from complex lipidomic data.

A Critical Overview of HPLC-MS-Based Lipidomics in Determining Triacylglycerol and Phospholipid in Foods

With the current advancement in mass spectrometry (MS)-based lipidomics, the knowledge of lipidomes and their diverse roles has greatly increased, enabling a deeper understanding of the action of bioactive lipid molecules in plant- and animal-based foods. This review provides in-depth information on the practical use of MS techniques in lipidomics, including lipid extraction, adduct formation, MS analysis, data processing, statistical analysis, and bioinformatics. Moreover, this contribution demonstrates the effectiveness of MS-based lipidomics for identifying and quantifying diverse lipid species, especially triacylglycerols and phospholipids, in foods. Further, it summarizes the wide applications of MS-based lipidomics in food science, such as for assessing food processing methods, detecting food adulteration, and measuring lipid oxidation in foods. Thus, MS-based lipidomics may be a useful method for identifying the action of individual lipid species in foods.

Adipose triglyceride lipase-mediated lipid catabolism is essential for bronchiolar regeneration

The lung airways are constantly exposed to inhaled toxic substances, resulting in cellular damage that is repaired by local expansion of resident bronchiolar epithelial club cells. Disturbed bronchiolar epithelial damage repair lies at the core of many prevalent lung diseases, including chronic obstructive pulmonary disease, asthma, pulmonary fibrosis, and lung cancer. However, it is still not known how bronchiolar club cell energy metabolism contributes to this process. Here, we show that adipose triglyceride lipase (ATGL), the rate-limiting enzyme for intracellular lipolysis, is critical for normal club cell function in mice. Deletion of the gene encoding ATGL, Pnpla2 (also known as Atgl), induced substantial triglyceride accumulation, decreased mitochondrial numbers, and decreased mitochondrial respiration in club cells. This defect manifested as bronchiolar epithelial thickening and increased airway resistance under baseline conditions. After naphthalene‑induced epithelial denudation, a regenerative defect was apparent. Mechanistically, dysfunctional PPARα lipid-signaling underlies this phenotype because (a) ATGL was needed for PPARα lipid-signaling in regenerating bronchioles and (b) administration of the specific PPARα agonist WY14643 restored normal bronchiolar club cell ultrastructure and regenerative potential. Our data emphasize the importance of the cellular energy metabolism for lung epithelial regeneration and highlight the significance of ATGL-mediated lipid catabolism for lung health.

Identifying active xenobiotics in humans by use of a suspect screening technique coupled with lipidomic analysis

Lipidomic analysis has been proven to be a powerful technique to explore the underlying associations between xenobiotics and health status of organisms. Here, we established a strategy that combined the lipidomic analysis with high-throughput suspect contaminant screening technique with an aim to efficiently identify active xenobiotics in humans. Firstly, in the light of single liquid phase equilibrium of chloroform-methanol-water (15:14:2, v/v/v), we developed an efficient method that was able to simultaneously extract both polar and nonpolar lipids in serum samples. By use of this method, targeted and non-targeted lipid analyses were conducted for n = 120 serum samples collected from Wuxi city, China. Secondly, we established a suspect database containing 1450 contaminants that have been previously reported in human samples, and contaminants in this database were screened in the same batch of serum samples by use of high-resolution mass spectrometry (HR-MS). Thirdly, the underlying associations between suspect contaminants and lipids were explored and discussed, and we observed that levels of some lipids were statistically correlated with concentrations of numerous contaminants. Among these active contaminants, 23 ones were identified on the basis of HR MS¹ and MS² characteristics, and these contaminants belonged to the classes of phthalates, phenols, parabens, or perfluorinated compounds (PFCs). Three active xenobiotics were fully validated by comparison with authentic standards, and they were perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS), and diethyl phthalate (DEP). There were statistically significant changes in levels of triglyceride (TG), lysophosphocholine (LPC), and sphingomyelin (SM) as peak areas of xenobiotics increase. We also observed that, among target lipid molecules, 18:0 lysophosphatidylethanolamine (LPE(18:0)) was very sensitive, and this lipid responded to exposure of various contaminants. Our present study provides novel knowledge on potential alteration of lipid metabolism in humans following exposure to xenobiotics, and provides an efficient strategy for efficiently identifying active xenobiotics in humans.

Computational solutions in redox lipidomics - Current strategies and future perspectives

The high chemical diversity of lipids allows them to perform multiple biological functions ranging from serving as structural building blocks of biological membranes to regulation of metabolism and signal transduction. In addition to the native lipidome, lipid species derived from enzymatic and non-enzymatic modifications (the epilipidome) make the overall picture even more complex, as their functions are still largely unknown. Oxidized lipids represent the fraction of epilipidome which has attracted high scientific attention due to their apparent involvement in the onset and development of numerous human disorders. Development of high-throughput analytical methods such as liquid chromatography coupled on-line to mass spectrometry provides the possibility to address epilipidome diversity in complex biological samples. However, the main bottleneck of redox lipidomics, the branch of lipidomics dealing with the characterization of oxidized lipids, remains the lack of optimal computational tools for robust, accurate and specific identification of already discovered and yet unknown modified lipids. Here we discuss the main principles of high-throughput identification of lipids and their modified forms and review the main software tools currently available in redox lipidomics. Different levels of confidence for software assisted identification of redox lipidome are defined and necessary steps toward optimal computational solutions are proposed.

Qualitative analysis of phospholipids and their oxidised derivatives - used techniques and examples of their applications related to lipidomic research and food analysis

Phospholipids (PLs) are important biomolecules that not only constitute structural building blocks and scaffolds of cell and organelle membranes but also play a vital role in cell biochemistry and physiology. Moreover, dietary exogenous PLs are characterised by high nutritional value and other beneficial health effects, which are confirmed by numerous epidemiological studies. For this reason, PLs are of high interest in lipidomics that targets both the analysis of membrane lipid distribution as well as correlates composition of lipids with their effects on functioning of cells, tissues and organs. Lipidomic assessments follow-up the changes occurring in living organisms, such as free radical attack and oxidative modifications of the polyunsaturated fatty acids (PUFAs) build in PL structures. Oxidised PLs (oxPLs) can be generated exogenously and supplied to organisms with processed food or formed endogenously as a result of oxidative stress. Cellular and tissue oxPLs can be a biomarker predictive of the development of numerous diseases such as atherosclerosis or neuroinflammation. Therefore, suitable high-throughput analytical techniques, which enable comprehensive analysis of PL molecules in terms of the structure of hydrophilic group, fatty acid (FA) composition and oxidative modifications of FAs, have been currently developed. This review addresses all aspects of PL analysis, including lipid isolation, chromatographic separation of PL classes and species, as well as their detection. The bioinformatic tools that enable handling of a large amount of data generated during lipidomic analysis are also discussed. In addition, imaging techniques such as confocal microscopy and mass spectrometry imaging for analysis of cellular lipid maps, including membrane PLs, are presented.

Spatially resolved endogenous improved metabolite detection in human osteoarthritis cartilage by matrix assisted laser desorption ionization mass spectrometry imaging

Osteoarthritis (OA) is one of the most common musculoskeletal diseases, characterized by the progressive deterioration of articular cartilage. Although the disease has been well studied in the past few years, the endogenous metabolic composition and more importantly the spatial information of these molecules in cartilage is still poorly understood. Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI) has been previously used for the investigation of the bimolecular distribution of proteins and lipids through the in situ analysis of cartilage tissue sections. MALDI-MSI as a tool to detect metabolites remains challenging, as these species have low abundance and degrade rapidly. In this work, we present a complete methodology, from sample preparation to data analysis for the detection of endogenous metabolites on cartilage by MSI. Our results demonstrate for the first time the ability to detect small molecules in fragile, challenging tissues through an optimized protocol, and render MSI as a tool towards a better understanding of OA.

Extraction of Plant Lipids for LC-MS-Based Untargeted Plant Lipidomics

Lipids are defined as hydrophobic or amphipathic small molecules which consist of a number of structurally and functionally distinct molecules that span from nonpolar to neutral to polar compounds. Lipidomics is the comprehensive analysis of all lipids in a biological system. Changes in lipid metabolism and composition, as well as of distinct lipid species have been linked with altered plant growth, development, and responses to environmental stresses including salinity. Recently, improved liquid chromatography mass spectrometry (LC-MS)-based techniques have provided the rapid expansion of lipidomics research. Sample preparation and lipid extraction are important steps in lipidomics, and this chapter describes important considerations in lipid monophasic and biphasic extractions from plant tissues prior to untargeted plant lipidomics approaches with LC-MS.

LipidHunter Identifies Phospholipids by High-Throughput Processing of LC-MS and Shotgun Lipidomics Datasets

Lipids are dynamic constituents of biological systems, rapidly responding to any changes in physiological conditions. Thus, there is a large interest in lipid-derived markers for diagnostic and prognostic applications, especially in translational and systems medicine research. As lipid identification remains a bottleneck of modern untargeted lipidomics, we developed LipidHunter, a new open source software for the high-throughput identification of phospholipids in data acquired by LC-MS and shotgun experiments. LipidHunter resembles a workflow of manual spectra annotation. Lipid identification is based on MS/MS data analysis in accordance with defined fragmentation rules for each phospholipid (PL) class. The software tool matches product and neutral loss signals obtained by collision-induced dissociation to a user-defined white list of fatty acid residues and PL class-specific fragments. The identified signals are tested against elemental composition and bulk identification provided via LIPID MAPS search. Furthermore, LipidHunter provides information-rich tabular and graphical reports allowing to trace back key identification steps and perform data quality control. Thereby, 202 discrete lipid species were identified in lipid extracts from rat primary cardiomyocytes treated with a peroxynitrite donor. Their relative quantification allowed the monitoring of dynamic reconfiguration of the cellular lipidome in response to mild nitroxidative stress. LipidHunter is available free for download at https://bitbucket.org/SysMedOs/lipidhunter .

A lipidomics investigation into the intervention of celastrol in experimental colitis

Celastrol is well known for its anti-inflammatory and anti-cancer effects. In this study, the efficacy of celastrol against dextran sulfate sodium (DSS)-induced inflammatory bowel disease (IBD) in mice was established and the mechanism was investigated using lipidomics. Celastrol treatment significantly alleviated DSS-induced colitis in mice, as revealed by the body weight, colon length, scores of rectal bleeding and diarrhea, serum TNF-α level, and histological analysis results. Lipidomics analysis based on UPLC/MS revealed characteristic changes in the metabolic profiles of the colitis mice, with altered levels of lipid markers associated with IBD, including LPC18 : 0, LPC18 : 1, LPC18 : 2, sphingomyelin (SM), and increased LPC18 : 0/LPC18 : 1 and LPC18 : 0/LPC18 : 2 ratios. For the celastrol-treated colitis mice, however, levels of the above lipid markers were restored, together with recovered saturated LPC/unsaturated LPC ratios. Accordingly, using GC-MS analysis, increased stearic acid (C18 : 0)/oleic acid (C18 : 1) and stearic acid (C18 : 0)/linoleic acid (C18 : 2) ratios were observed in colitis mice, which were later recovered after celastrol treatment. Quantitative real-time PCR analysis revealed that the liver expression of stearoyl-coenzyme A desaturase 1 (SCD1), the key enzyme controlling the desaturation of saturated fatty acid, was dramatically inhibited in IBD mice, and was obviously recovered after celastrol treatment. These results suggest that the increased saturated LPC/unsaturated LPC (and saturated fatty acid/unsaturated fatty acid) ratios associated with SCD1 down-regulation could be regarded as biomarkers of colitis, and celastrol alleviates DSS-induced colitis partially via up-regulation of SCD1, restoring the altered balance between stearic acid- and oleic acid-derived lipid species, which plays an important role in alleviating colitis. In all, this study provided the scientific basis for further development of celastrol in treating IBD.

Plant lipidomics at the crossroads: From technology to biology driven science

The identification and quantification of lipids from plant tissues have become commonplace and many researchers now incorporate lipidomics approaches into their experimental studies. Plant lipidomics research continues to involve technological developments such as those in mass spectrometry imaging, but in large part, lipidomics approaches have matured to the point of being accessible to the novice. Here we review some important considerations for those planning to apply plant lipidomics to their biological questions, and offer suggestions for appropriate tools and practices. This article is part of a Special Issue entitled: BBALIP_Lipidomics Opinion Articles edited by Sepp Kohlwein.

On Mass Ambiguities in High-Resolution Shotgun Lipidomics

Mass-spectrometry-based lipidomics aims to identify as many lipid species as possible from complex biological samples. Due to the large combinatorial search space, unambiguous identification of lipid species is far from trivial. Mass ambiguities are common in direct-injection shotgun experiments, where an orthogonal separation (e.g., liquid chromatography) is missing. Using the rich information within available lipid databases, we generated a comprehensive rule set describing mass ambiguities, while taking into consideration the resolving power (and its decay) of different mass analyzers. Importantly, common adduct species and isotopic peaks are accounted for and are shown to play a major role, both for perfect mass overlaps due to identical sum formulas and resolvable mass overlaps. We identified known and hitherto unknown mass ambiguities in high- and ultrahigh resolution data, while also ranking lipid classes by their propensity to cause ambiguities. On the basis of this new set of ambiguity rules, guidelines and recommendations for experimentalists and software developers of what constitutes a solid lipid identification in both MS and MS/MS were suggested. For researchers new to the field, our results are a compact source of ambiguities which should be accounted for. These new findings also have implications for the selection of internal standards, peaks used for internal mass calibration, optimal choice of instrument resolution, and sample preparation, for example, in regard to adduct ion formation.

α/β Hydrolase Domain-containing 6 (ABHD6) Degrades the Late Endosomal/Lysosomal Lipid Bis(monoacylglycero)phosphate

α/β Hydrolase domain-containing 6 (ABHD6) can act as monoacylglycerol hydrolase and is believed to play a role in endocannabinoid signaling as well as in the pathogenesis of obesity and liver steatosis. However, the mechanistic link between gene function and disease is incompletely understood. Here we aimed to further characterize the role of ABHD6 in lipid metabolism. We show that mouse and human ABHD6 degrade bis(monoacylglycero)phosphate (BMP) with high specific activity. BMP, also known as lysobisphosphatidic acid, is enriched in late endosomes/lysosomes, where it plays a key role in the formation of intraluminal vesicles and in lipid sorting. Up to now, little has been known about the catabolism of this lipid. Our data demonstrate that ABHD6 is responsible for ∼90% of the BMP hydrolase activity detected in the liver and that knockdown of ABHD6 increases hepatic BMP levels. Tissue fractionation and live-cell imaging experiments revealed that ABHD6 co-localizes with late endosomes/lysosomes. The enzyme is active at cytosolic pH and lacks acid hydrolase activity, implying that it degrades BMP exported from acidic organelles or de novo-formed BMP. In conclusion, our data suggest that ABHD6 controls BMP catabolism and is therefore part of the late endosomal/lysosomal lipid-sorting machinery.

Oxidative lipidomics coming of age: advances in analysis of oxidized phospholipids in physiology and pathology

SIGNIFICANCE

Oxidized phospholipids are now well recognized as markers of biological oxidative stress and bioactive molecules with both pro-inflammatory and anti-inflammatory effects. While analytical methods continue to be developed for studies of generic lipid oxidation, mass spectrometry (MS) has underpinned the advances in knowledge of specific oxidized phospholipids by allowing their identification and characterization, and it is responsible for the expansion of oxidative lipidomics.

RECENT ADVANCES

Studies of oxidized phospholipids in biological samples, from both animal models and clinical samples, have been facilitated by the recent improvements in MS, especially targeted routines that depend on the fragmentation pattern of the parent molecular ion and improved resolution and mass accuracy. MS can be used to identify selectively individual compounds or groups of compounds with common features, which greatly improves the sensitivity and specificity of detection. Application of these methods has enabled important advances in understanding the mechanisms of inflammatory diseases such as atherosclerosis, steatohepatitis, leprosy, and cystic fibrosis, and it offers potential for developing biomarkers of molecular aspects of the diseases.

CRITICAL ISSUES AND FUTURE DIRECTIONS

The future in this field will depend on development of improved MS technologies, such as ion mobility, novel enrichment methods and databases, and software for data analysis, owing to the very large amount of data generated in these experiments. Imaging of oxidized phospholipids in tissue MS is an additional exciting direction emerging that can be expected to advance understanding of physiology and disease.

Morphogenesis checkpoint kinase Swe1 is the executor of lipolysis-dependent cell-cycle progression

Cell growth and division requires the precise duplication of cellular DNA content but also of membranes and organelles. Knowledge about the cell-cycle-dependent regulation of membrane and storage lipid homeostasis is only rudimentary. Previous work from our laboratory has shown that the breakdown of triacylglycerols (TGs) is regulated in a cell-cycle-dependent manner, by activation of the Tgl4 lipase by the major cyclin-dependent kinase Cdc28. The lipases Tgl3 and Tgl4 are required for efficient cell-cycle progression during the G1/S (Gap1/replication phase) transition, at the onset of bud formation, and their absence leads to a cell-cycle delay. We now show that defective lipolysis activates the Swe1 morphogenesis checkpoint kinase that halts cell-cycle progression by phosphorylation of Cdc28 at tyrosine residue 19. Saturated long-chain fatty acids and phytosphingosine supplementation rescue the cell-cycle delay in the Tgl3/Tgl4 lipase-deficient strain, suggesting that Swe1 activity responds to imbalanced sphingolipid metabolism, in the absence of TG degradation. We propose a model by which TG-derived sphingolipids are required to activate the protein phosphatase 2A (PP2A(Cdc55)) to attenuate Swe1 phosphorylation and its inhibitory effect on Cdc28 at the G1/S transition of the cell cycle.

PPARα-dependent exacerbation of experimental colitis by the hypolipidemic drug fenofibrate

Fibrates, such as fenofibrate, are peroxisome proliferator-activated receptor-α (PPARα) agonists and have been used for several decades as hypolipidemic agents in the clinic. However, contradictory observations exist on the role of fibrates in host response to acute inflammation, with unclear mechanisms. The role of PPARα in colitis was assessed using fenofibrate and Ppara-null mice. Wild-type or Ppara-null mice were subjected to acute colitis under three distinct protocols, dextran sulfate sodium, trinitrobenzenesulfonic acid, and Salmonella Typhi. Serum and colon lipidomics were analyzed to characterize the metabolic profiles by ultra-performance liquid chromatography-coupled with electrospray ionization quadrupole time-of-flight mass spectrometry. Messenger RNAs of PPARα target genes and genes involved in inflammation were determined by qunatitative PCR analysis. Fenofibrate treatment exacerbated inflammation and tissue injury in acute colitis, and this was dependent on PPARα activation. Lipidomics analysis revealed that bioactive sphingolipids, including sphingomyelins (SM) and ceramides, were significantly increased in the colitis group compared with the control group; this was further potentiated following fenofibrate treatment. In the colon, fenofibrate did not reduce the markedly increased expression of mRNA encoding TNFα found in the acute colitis model, while it decreased hydrolysis and increased synthesis of SM, upregulated RIPK3-dependent necrosis, and elevated mitochondrial fatty acid β-oxidation, which were possibly related to the exacerbated colitis.

Bile acid signaling in lipid metabolism: metabolomic and lipidomic analysis of lipid and bile acid markers linked to anti-obesity and anti-diabetes in mice

Bile acid synthesis is the major pathway for catabolism of cholesterol. Cholesterol 7α-hydroxylase (CYP7A1) is the rate-limiting enzyme in the bile acid biosynthetic pathway in the liver and plays an important role in regulating lipid, glucose and energy metabolism. Transgenic mice overexpressing CYP7A1 (CYP7A1-tg mice) were resistant to high-fat diet (HFD)-induced obesity, fatty liver, and diabetes. However the mechanism of resistance to HFD-induced obesity of CYP7A1-tg mice has not been determined. In this study, metabolomic and lipidomic profiles of CYP7A1-tg mice were analyzed to explore the metabolic alterations in CYP7A1-tg mice that govern the protection against obesity and insulin resistance by using ultra-performance liquid chromatography-coupled with electrospray ionization quadrupole time-of-flight mass spectrometry combined with multivariate analyses. Lipidomics analysis identified seven lipid markers including lysophosphatidylcholines, phosphatidylcholines, sphingomyelins and ceramides that were significantly decreased in serum of HFD-fed CYP7A1-tg mice. Metabolomics analysis identified 13 metabolites in bile acid synthesis including taurochenodeoxycholic acid, taurodeoxycholic acid, tauroursodeoxycholic acid, taurocholic acid, and tauro-β-muricholic acid (T-β-MCA) that differed between CYP7A1-tg and wild-type mice. Notably, T-β-MCA, an antagonist of the farnesoid X receptor (FXR) was significantly increased in intestine of CYP7A1-tg mice. This study suggests that reducing 12α-hydroxylated bile acids and increasing intestinal T-β-MCA may reduce high fat diet-induced increase of phospholipids, sphingomyelins and ceramides, and ameliorate diabetes and obesity. This article is part of a Special Issue entitled Linking transcription to physiology in lipodomics.

Lipidomics: analysis of the lipid composition of cells and subcellular organelles by electrospray ionization mass spectrometry

Lipidomics aims to quantitatively define lipid classes, including their molecular species, in biological systems. Lipidomics has experienced rapid progress, mainly because of continuous technical advances in instrumentation that are now enabling quantitative lipid analyses with an unprecedented level of sensitivity and precision. The still-growing category of lipids includes a broad diversity of chemical structures with a wide range of physicochemical properties. Reflecting this diversity, different methods and strategies are being applied to the quantification of lipids. Here, I review state-of-the-art electrospray ionization tandem mass spectrometric approaches and direct infusion to quantitatively assess lipid compositions of cells and subcellular fractions. Finally, I discuss a few examples of the power of mass spectrometry-based lipidomics in addressing cell biological questions.

Emerging targets in lipid-based therapy

The use of prostaglandins and NSAIDS in the clinic has proven that lipid mediators and their associated pathways make attractive therapeutic targets. When contemplating therapies involving lipid pathways, several basic agents come to mind. There are the enzymes and accessory proteins that lead to the metabolism of lipid substrates, provided through diet or through actions of lipases, the subsequent lipid products, and finally the lipid sensors or receptors. There is abundant evidence that molecules along this lipid continuum can serve as prognostic and diagnostic indicators and are in fact viable therapeutic targets. Furthermore, lipids themselves can be used as therapeutics. Despite this, the vernacular dialog pertaining to "biomarkers" does not routinely include mention of lipids, though this is rapidly changing. Collectively these agents are becoming more appreciated for their respective roles in diverse disease processes from cancer to preterm labor and are receiving their due appreciation after decades of ground work in the lipid field. By relating examples of disease processes that result from dysfunction along the lipid continuum, as well as examples of lipid therapies and emerging technologies, this review is meant to inspire further reading and discovery.

Omics strategies for revealing Yersinia pestis virulence

Omics has remarkably changed the way we investigate and understand life. Omics differs from traditional hypothesis-driven research because it is a discovery-driven approach. Mass datasets produced from omics-based studies require experts from different fields to reveal the salient features behind these data. In this review, we summarize omics-driven studies to reveal the virulence features of Yersinia pestis through genomics, trascriptomics, proteomics, interactomics, etc. These studies serve as foundations for further hypothesis-driven research and help us gain insight into Y. pestis pathogenesis.