Abstract
Introduction:
The existing evidence indicates that atherosclerosis (AS) plays a pivotal role in the progression and exacerbation of cardiovascular diseases and their associated complications. Current diagnostic and therapeutic strategies for atherosclerosis are limited in their ability to facilitate early detection and personalized treatment. This study employs a systems biology approach to investigate the role of lactylation-related genes (LRGs) in the pathogenesis of atherosclerosis, while considering the well-established correlation between inflammatory responses and atherosclerosis development.
Methods:
In this study, we utilized datasets obtained from the Gene Expression Omnibus (GEO) as well as data from previous studies on lactylation-related genes (LRGs). Following this, we identified 17 lactylation related genes associate with atherosclerosis (AS-LRGs) from the GSE100927 dataset. Subsequently, we employed the validation dataset (GSE43292) to assess these 17 AS-LRGs, resulting in the identification of 12 more reliable candidate genes. These genes were further analyzed for functional enrichment through Gene Ontology (GO) annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, and gene set enrichment analysis (GSEA). To elucidate the potential utility of AS-LRGs in diagnosing high-risk plaques, we assessed their expression in both early and late stages of atherosclerosis, as well as in high- and low-risk plaques. We then constructed interaction networks to elucidate the potential regulatory relationships among LRGs, miRNAs, transcription factors, and drugs. Finally, we utilized the single sample Gene Set Enrichment Analysis (ssGSEA) method to investigate immune infiltration in AS and evaluate the levels of immune cell infiltration.
Results:
We identified 12 lactylation-related genes that are more reliably associated with atherosclerosis: five upregulated genes (LSP1, IKZF1, MNDA, RCC2, and WAS) and seven downregulated genes (CSRP2, PPP1CB, CSRP1, HEXIM1, CALD1, PDLIM1, and RANBP2).
Discussion:
This study elucidates the pivotal role of lactylation in atherosclerosis (AS) and establishes a robust foundation for future research into targeted therapies and clinical applications of the identified biomarkers.
1 Introduction
Atherosclerosis is a progressive condition characterized by diverse structural modifications in the walls of large and medium-sized arteries, ultimately resulting in the development of atherosclerotic plaques (Camaré et al., 2017). It poses a significant threat to both individual patients and public health, as it is one of the primary causes of cardiovascular diseases. As demonstrated by prior studies, atherosclerotic cardiovascular disease constitutes a significant global contributor to mortality (Barquera et al., 2015). Furthermore, atherosclerotic plaques serve as the pathophysiological foundation for nearly all arterial vascular diseases (Pi et al., 2024). Current diagnostic and therapeutic approaches for atherosclerosis include pharmacological treatments, lifestyle modifications, and surgical interventions; however, significant limitations remain, such as inaccurate diagnoses and suboptimal therapeutic outcomes, which necessitate further attention (Libby et al., 2019; Shamaki et al., 2022).
Lactylation, which is regulated by lactate concentrations, has recently been recognized as a novel component of the epigenetic landscape. This finding not only facilitates extensive research into lactate metabolism but also offers essential reference points for subsequent functional and mechanistic studies (Liu et al., 2022). Prior research has demonstrated that lactylation modification is crucial in the pathogenesis of various diseases, including cancer (Li F. et al., 2024) and cardiovascular disorders (Wang N. et al., 2022). Integrating the distribution and metabolic characteristics of cells within atherosclerotic plaques, lactylation modification may present a novel strategy for targeted interventions in atherosclerosis (Xu et al., 2023). However, the role of lactylation modification in atherosclerosis remains a topic of active debate within the scientific community, with varying perspectives on its significance and mechanisms. Recent studies have demonstrated that lactic acid-mediated lactylation modification promotes the development of atherosclerosis (Dong et al., 2024); however, other research suggests that under specific conditions, lactylation modification may exert an inhibitory effect on atherosclerotic progression (Zhang et al., 2024).
In this study, we systematically compiled published data related to lactylation and conducted an in-depth analysis of gene expression levels across multiple atherosclerotic databases to identify potential diagnostic markers and therapeutic targets. Additionally, we elucidated a range of immune cells associated with atherosclerosis through bioinformatics analysis, highlighting their significant roles in disease progression and their interactions with immune responses. In conclusion, the intersection of lactylation and atherosclerosis represents a promising frontier for exploration in contemporary biomedical research. By focusing on the expression and functional implications of lactylation-related genes (LRGs) in atherosclerosis (AS), this study aims to address existing knowledge gaps and lay the foundation for future investigations into targeted therapeutic strategies.
2 Materials and methods
2.1 Identification of differentially expressed genes between atherosclerotic and healthy arteries
A total of four AS-related gene expression matrices (GSE100927, GSE43292, GSE28829, and GSE163154) were obtained from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). In the GSE100927 dataset (Steenman et al., 2018), atherosclerotic plaques and control arteries without such lesions (obtained from deceased organ donors) were collected from the carotid, femoral, and infra-popliteal arteries. Total RNA was extracted from the samples using standard protocols and subsequently hybridized to microarrays for further analysis. The “limma (Ritchie et al., 2015)” R package was employed to normalize the data and identify differentially expressed genes (DEGs). The significance criteria were defined as −0.58 < log2 fold change (FC) < 0.58, accompanied by a p-value of less than 0.05. In R version 4.2.1, the “ggplot2” and “ComplexHeatmap (Gu et al., 2016)” packages were utilized to generate volcano plots and heatmaps for the differentially expressed genes (DEGs). A list of 332 lactylation-related genes (LRGs) was obtained from previously published reports (Cheng et al., 2023). The collected LRG data were then integrated with the results of differential expression analysis of genes (DEGs) between the disease and control groups in the GSE100927 dataset, and these genes were classified as AS-related LRGs (AS-LRGs). The results were visualized utilizing the “ggplot2” and “VennDiagram” packages.
2.2 Analysis of LRGs expression in the validation dataset and receiver operating characteristic (ROC) curves
The GSE43292 dataset was utilized as the validation dataset. In GSE43292 dataset (Ayari and Bricca, 2013), a total of 34 patients who underwent carotid endarterectomy at the University Hospital of Lyon participated in this study. Carotid endarterectomy samples were obtained in the operating room and immediately divided into two fragments: the atheroma plaque and macroscopically intact tissue. We evaluated the expression levels of AS-LRGs in both affected and control groups, using the Wilcoxon rank-sum test for statistical analysis. The ROC analysis on the dataset was performed using the “pROC (Robin et al., 2011)” package. Statistical analysis of paired samples was conducted using a paired t-test. All visualizations were generated using the “ggplot2” package.
2.3 Functional and pathway enrichment analysis
Functional and pathway enrichment analyses were conducted using the Gene Ontology (2015) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2022) databases, employing the R package “clusterProfiler (Yu et al., 2012)” based on significant AS-LRGs (P < 0.05).
2.4 Gene set enrichment analysis (GSEA)
Gene sets redefined from the MSigDB database (https://www.gsea-msigdb.org/gsea/msigdb/collections.jsp) were utilized to analyze the distribution patterns of genes within a gene expression profile. This profile is ranked based on the correlation between gene expression and phenotypic characteristics, thereby evaluating the contributions and associations of these genes with the phenotype (Subramanian et al., 2005). We retrieved the gene sets “c2.cp.kegg.v2022.1.Hs.symbols” and “c5.all.v2022.1.Hs.symbols” from the MSigDB database and utilized the “clusterProfiler (Yu et al., 2012)” package in R to perform the analysis.
2.5 Construction of interaction networks involving AS-LRGs, associated microRNAs, transcription factors, and drugs
The TarBase v9.0 (Skoufos et al., 2024) and ENCODE (2012) databases were employed to identify the miRNAs and transcription factors associated with AS-LRGs, which were subsequently visualized using the NetworkAnalyst (2012) platform. We utilized the Drug Gene Interaction Database (DGIdb) (Cannon et al., 2024) (https://www.dgidb.org) to predict the relationships between genes and drugs, which were visualized using Cytoscape (Shannon et al., 2003).
2.6 Analysis of immune infiltration
The immune infiltration within the dataset was assessed using the single-sample GSEA (ssGSEA) algorithm from the R package GSVA (Hänzelmann et al., 2013), incorporating 24 immune cell markers detailed in an Immunity article (Bindea et al., 2013). We conducted the Wilcoxon rank-sum test to examine the differences between the normal group and the AS group. Spearman’s correlation analysis was conducted to assess the relationship between AS-LRGs and the level of immune infiltration in AS tissues. All visualizations were generated utilizing “ggplot2” package.
2.7 Analysis of AS-LRGs expression in the GSE28829 and GSE163154 datasets, with ROC curve evaluation
In the GSE28829 database, Seventeen samples of atherosclerotic tissue were collected from carotid artery segments, including nine samples representing early lesions and eight samples representing advanced lesions, sourced from the Maastricht Pathology Tissue Collection for inclusion in this study (Döring et al., 2012). The GSE28829 database was utilized to analyze the differential expression of AS-LRGs between early and advanced lesions. In GSE163154 database, they conducted an analysis of the genetic profiles of atherosclerotic lesion segments categorized as low-risk or high-risk in patients undergoing carotid endarterectomy, based on the presence or absence of intraplaque hemorrhage (Jin et al., 2021). The GSE163154 database was employed to analyze the differential expression of AS-LRGs in the absence and presence of intraplaque hemorrhage. The statistical and presentation methodologies utilized in this study are consistent with those previously described.
3 Results
3.1 Identification of lactylation-related genes in atherosclerosis
Figure 1 provides a comprehensive overview of the workflow. The differentially expressed genes (DEGs) in the GSE100927 dataset were visualized using volcano plots (Figure 2A). As observed, there were 1,201 upregulated DEGs and 765 downregulated DEGs. A classification heatmap was created to clearly present these results (Figure 2B). Additionally, by integrating the LRGs with the upregulated and downregulated DEGs from the GSE100927 dataset, we identified 9 upregulated genes (FABP5, LSP1, IKZF1, MNDA, CRABP2, ARID3A, RCC2, WAS, and HMGA1) and 8 downregulated genes (CSRP2, PPP1CB, CSRP1, HIST1H1C, HEXIM1, CALD1, PDLIM1, and RANBP2) that are common to both DEGs and LRGs and are associated with atherosclerosis (Figures 2C, D). Subsequently, we assessed the expression levels of these genes in the GSE43292 dataset. Among the identified genes, FABP5 and HIST1H1C were not detected as independent samples in the GSE43292 dataset. Subsequently, the expression of the 8 upregulated genes and 7 downregulated genes in the GSE43292 dataset were illustrated using box plots (Figures 3A, B). In the validation dataset (GSE43292), the 5 upregulated genes identified as having diagnostic value were: LSP1, IKZF1, MNDA, RCC2, and WAS (AUC >0.7) (Figures 3C, D); the 7 downregulated genes identified as having diagnostic value were: CSRP2, PPP1CB, CSRP1, HEXIM1, CALD1, PDLIM1 and RANBP2 (AUC >0.7) (Figures 3E, F). Subsequently, we conducted paired-sample t-tests on the 5 upregulated genes and 7 downregulated genes within the validation dataset, resulting in significant differences for all genes (Figures 4A–L).
FIGURE 1
FIGURE 2
FIGURE 3
FIGURE 4
3.2 Functional enrichment analysis of AS-LRGs and gene set enrichment analysis (GSEA)
The gene ontology (GO) analysis revealed significant enrichment in a multitude of biological processes, cellular components, and molecular functions. The KEGG analysis of the AS-LRGs revealed significant enrichment across a diverse array of pathways. The results of the GO and KEGG analyses for AS-LRGs are presented in Figure 5A. To evaluate the impact of gene expression levels on AS, Gene Set Enrichment Analysis (GSEA) was conducted using the GSE100927 dataset to elucidate the relationships between gene expression and related biological processes, cellular components, and molecular functions. The findings indicated that in the GSE100927 dataset, the genes were primarily associated with contractile fiber, I band, adaptive immune response, leukocyte-mediated immunity, immune response-regulating signaling pathways, and other biological functions (Figure 5B). The biological pathways predominantly regulated by genes in the GSE100927 dataset include those associated with lysosome, allograft rejection, antigen processing and presentation, leishmania infection, autoimmune thyroid disease, natural killer cell mediated cytotoxicity, cell adhesion molecules, chemokine signaling pathway, cytokine-cytokine receptor interaction, graft versus host disease, asthma, and type I diabetes mellitus (Figures 5C–H).
FIGURE 5
3.3 Network analysis of AS-LRGs and their associated miRNAs, transcription factors, and drugs
We constructed an interaction network between AS-LRGs and miRNAs, encompassing 12 genes and 561 miRNAs (Figure 6A). The top five predictions of the AS-LRG-miRNA interaction network included: RANBP2 (targeted by 294 miRNAs), PPP1CB (targeted by 276 miRNAs), CALD1 (targeted by 170 miRNAs), RCC2 (targeted by 164 miRNAs) and HEXIM1 (targeted by 146 miRNAs). The interaction network between AS-LRGs and transcription factors (TFs) consisted of 10 genes and 212 transcription factors (Figure 6B). The top five predictions of the AS-LRG-TF interaction network encompassed: HEXIM1 (targeted by 176 TFs), PPP1CP (targeted by 28 TFs), CSRP1 (targeted by 21 TFs), CALD1 (targeted by 19 TFs) and WAS (targeted by 16 TFs). The AS-LRG–drug interaction network consisted of 11 drugs and 4 genes, as illustrated in Figure 6C.
FIGURE 6
3.4 Analysis of the relationship between immune infiltration and AS
To achieve a more comprehensive understanding of the role of immune cells in AS and their association with critical genes, we employed the single-sample GSEA methodology to assess immune infiltration within the GSE100927 dataset. Among the 24 immune cell subpopulations, we observed significant differences between AS-affected tissues and normal tissues. The violin plot illustrated that AS patients exhibited varying levels of the majority of immune cell types (Figure 7A). Five upregulated genes and seven downregulated genes were identified as being associated with each immune cell type, as illustrated in Figures 7B, C. The chart indicates that the majority of upregulated genes exhibited a positive correlation with most immune cell types, with the exception of NK cells. Conversely, the chart shows that the majority of downregulated genes exhibited a negative correlation with most immune cell types, with the exception of NK cells.
FIGURE 7
3.5 Analysis of AS-LRGs in early and advanced lesions
To evaluate the differential expression of AS-LRGs in early and advanced lesions, we utilized the GSE28829 dataset, which consists of 9 early lesion samples and 8 advanced lesion samples. The expression levels of the five upregulated genes and seven downregulated genes in the GSE28829 dataset were illustrated using box plots (Figures 8A, B). In the GSE28829 dataset, the five upregulated genes identified as having diagnostic value were LSP1, IKZF1, MNDA, RCC2, and WAS (AUC >0.7) (Figure 8C). The five downregulated genes identified as having diagnostic value were CSRP2, PPP1CB, CSRP1, HEXIM1, and CALD1 (AUC >0.7) (Figure 8D).
FIGURE 8
3.6 Analysis of AS-LRGs in the absence and presence of intraplaque hemorrhage
The GSE163154 dataset, which includes samples indicating the presence or absence of intraplaque hemorrhage, was utilized to assess the relationship between AS-LRGs and intraplaque hemorrhage. The expression levels of the five upregulated genes and seven downregulated genes in the GSE163154 dataset were illustrated using box plots (Figures 9A, B). In the GSE163154 dataset, the five upregulated genes identified as having diagnostic value were LSP1, IKZF1, MNDA, RCC2, and WAS (AUC >0.7) (Figures 9C, D). The six downregulated genes identified as having diagnostic value were CSRP2, PPP1CB, CSRP1, HEXIM1, CALD1 and PDLIM1 (AUC >0.7) (Figures 9E, F).
FIGURE 9
4 Discussion
Atherosclerosis, characterized by the development of fibrofatty lesions within the arterial wall, significantly contributes to global morbidity and mortality from related diseases. It accounts for the majority of myocardial infarctions and numerous strokes and is a leading cause of disabling peripheral artery disease (Libby et al., 2019). Lactate, traditionally considered a mere byproduct of glycolysis, has garnered increasing recognition in recent years for its significant biological functions (Ouyang et al., 2023). Similarly, the significance of lactate in atherosclerosis has increasingly been acknowledged. Strategies aimed at targeting lactate production, regulating its transport, and modulating circulating levels may emerge as promising therapeutic approaches for atherosclerosis in the future (Li X. et al., 2024).
This study investigates the differential expression of LRGs in the context of AS. Previous research has highlighted the critical role of lactylation in cellular metabolism and signal transduction (Wang J. et al., 2022); however, the specific contributions of LRGs to the pathophysiology of AS remain poorly understood. We conducted a comprehensive analysis of the functional enrichment of these AS-LRGs, along with their associated transcription factors (TFs), microRNAs (miRNAs), and drugs, aiming to establish a solid foundation for subsequent in-depth investigations into their mechanisms. Moreover, the analysis of immune responses revealed that LRGs play a complex and multifaceted role in modulating immune cell activity within the atherosclerotic microenvironment. The findings demonstrated a significant positive correlation between several upregulated LRGs and immune cell infiltration, particularly in macrophages and T cells, highlighting their critical role in the inflammatory processes characteristic of atherosclerosis (Kuznetsova et al., 2020; Saigusa et al., 2020).
Atherosclerotic plaques typically develop over a period spanning several years to decades (Vergallo and Crea, 2020).However, thrombotic complications associated with atherosclerotic disease can arise abruptly and frequently without prior indication (Libby, 2013). The transition of atherosclerotic plaque from a stable to an unstable phenotype is associated with an increased incidence of acute cardiovascular events (Jin et al., 2021). Consequently, the assessment of plaque stability is of paramount importance. In this study, we investigated the expression differences of five upregulated genes and seven downregulated AS-LRGs in early versus advanced plaques, as well as high-risk versus low-risk plaques, along with their corresponding diagnostic ROC curves. We identified five upregulated genes and five downregulated AS-LRGs, both of which demonstrated significant differences and diagnostic efficacy between early and advanced plaques. Similarly, we identified five upregulated genes and six downregulated AS-LRGs, which also exhibited significant differences and diagnostic efficacy between high-risk and low-risk plaques. Understanding the fundamental mechanisms governing the transition from low-risk to high-risk rupture-prone plaques, as well as from early to advanced stages in human atherosclerosis, will facilitate early diagnosis and enable targeted interventions for cardiovascular diseases associated with atherosclerosis. Our forthcoming investigations will comprehensively assess the roles of these genes in the progression of atherosclerosis, along with their associated clinical significance.
The limitations of this study encompass the reliance on publicly available gene expression datasets lacking wet-lab validation, a relatively small sample size that may restrict the generalizability of the findings, and the potential for batch effects arising from the utilization of multiple datasets. These considerations underscore the imperative for future research to include larger cohorts and experimental validation in order to reinforce the clinical significance of the identified biomarkers.
In conclusion, this study conducted a systematic analysis of the expression patterns of lactylation-related genes in atherosclerosis, elucidating their potential roles in disease progression. The identification of these genes not only deepens our understanding of the molecular mechanisms underlying atherosclerosis but also establishes a foundation for future research aimed at developing innovative diagnostic and therapeutic strategies. Ultimately, the insights derived from this work may significantly enhance early diagnosis and personalized treatment options for patients with atherosclerosis.
Statements
Data availability statement
All data and analyses presented in this study are sourced from public databases, with all references clearly cited in the text.
Ethics statement
All data utilized in this study were sourced from publicly available databases. The ethics statement has been revised in accordance with the specified requirements.
Author contributions
JQ: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. QZ: Conceptualization, Funding acquisition, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by a grant from the National Natural Science Foundation of China (Grant No. 82170496).
Acknowledgments
Thanks a lot for all members of our team.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
- AS
Atherosclerosis
- LRGs
Lactylation-related genes
- DEGs
Differentially expressed genes
- AS-LRGs
Lactylation-related genes associate with atherosclerosis
- GO
Gene Ontology
- KEGG
Kyoto Encyclopedia of Genes and Genomes
- ROC
Receiver Operating Characteristic
- GSEA
Gene set enrichment analysis
- ssGSEA
Single-sample GSEA
- TFs
Transcription factors
References
1
AyariH.BriccaG. (2013). Identification of two genes potentially associated in iron-heme homeostasis in human carotid plaque using microarray analysis. J. Biosci.38 (2), 311–315. 10.1007/s12038-013-9310-2
2
BarqueraS.Pedroza-TobíasA.MedinaC.Hernández-BarreraL.Bibbins-DomingoK.LozanoR.et al (2015). Global overview of the epidemiology of atherosclerotic cardiovascular disease. Arch. Med. Res.46 (5), 328–338. 10.1016/j.arcmed.2015.06.006
3
BindeaG.MlecnikB.TosoliniM.KirilovskyA.WaldnerM.ObenaufA. C.et al (2013). Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity39 (4), 782–795. 10.1016/j.immuni.2013.10.003
4
CamaréC.PucelleM.Nègre-SalvayreA.SalvayreR. (2017). Angiogenesis in the atherosclerotic plaque. Redox Biol.12, 18–34. 10.1016/j.redox.2017.01.007
5
CannonM.StevensonJ.StahlK.BasuR.CoffmanA.KiwalaS.et al (2024). DGIdb 5.0: rebuilding the drug-gene interaction database for precision medicine and drug discovery platforms. Nucleic Acids Res.52 (D1), D1227–d1235. 10.1093/nar/gkad1040
6
ChengZ.HuangH.LiM.LiangX.TanY.ChenY. (2023). Lactylation-related gene signature effectively predicts prognosis and treatment responsiveness in hepatocellular carcinoma. Pharm. (Basel)16 (5), 644. 10.3390/ph16050644
7
DongM.ZhangY.ChenM.TanY.MinJ.HeX.et al (2024). ASF1A-dependent P300-mediated histone H3 lysine 18 lactylation promotes atherosclerosis by regulating EndMT. Acta Pharm. Sin. B14 (7), 3027–3048. 10.1016/j.apsb.2024.03.008
8
DöringY.MantheyH. D.DrechslerM.LievensD.MegensR. T.SoehnleinO.et al (2012). Auto-antigenic protein-DNA complexes stimulate plasmacytoid dendritic cells to promote atherosclerosis. Circulation125 (13), 1673–1683. 10.1161/circulationaha.111.046755
9
ENCODE (2012). An integrated encyclopedia of DNA elements in the human genome. Nature489 (7414), 57–74. 10.1038/nature11247
10
Gene Ontology (2015). Gene ontology consortium: going forward. Nucleic Acids Res.43 (Database issue), D1049–D1056. 10.1093/nar/gku1179
11
GuZ.EilsR.SchlesnerM. (2016). Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics32 (18), 2847–2849. 10.1093/bioinformatics/btw313
12
HänzelmannS.CasteloR.GuinneyJ. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinforma.14, 7. 10.1186/1471-2105-14-7
13
JinH.GoossensP.JuhaszP.EijgelaarW.MancaM.KarelJ. M. H.et al (2021). Integrative multiomics analysis of human atherosclerosis reveals a serum response factor-driven network associated with intraplaque hemorrhage. Clin. Transl. Med.11 (6), e458. 10.1002/ctm2.458
14
KanehisaM.SatoY.KawashimaM. (2022). KEGG mapping tools for uncovering hidden features in biological data. Protein Sci.31 (1), 47–53. 10.1002/pro.4172
15
KuznetsovaT.PrangeK. H. M.GlassC. K.de WintherM. P. J. (2020). Transcriptional and epigenetic regulation of macrophages in atherosclerosis. Nat. Rev. Cardiol.17 (4), 216–228. 10.1038/s41569-019-0265-3
16
LiF.SiW.XiaL.YinD.WeiT.TaoM.et al (2024a). Positive feedback regulation between glycolysis and histone lactylation drives oncogenesis in pancreatic ductal adenocarcinoma. Mol. Cancer23 (1), 90. 10.1186/s12943-024-02008-9
17
LiX.ChenM.ChenX.HeX.LiX.WeiH.et al (2024b). TRAP1 drives smooth muscle cell senescence and promotes atherosclerosis via HDAC3-primed histone H4 lysine 12 lactylation. Eur. Heart J.45 (39), 4219–4235. 10.1093/eurheartj/ehae379
18
LibbyP. (2013). Mechanisms of acute coronary syndromes and their implications for therapy. N. Engl. J. Med.368 (21), 2004–2013. 10.1056/NEJMra1216063
19
LibbyP.BuringJ. E.BadimonL.HanssonG. K.DeanfieldJ.BittencourtM. S.et al (2019). Atherosclerosis. Nat. Rev. Dis. Prim.5 (1), 56. 10.1038/s41572-019-0106-z
20
LiuX.ZhangY.LiW.ZhouX. (2022). Lactylation, an emerging hallmark of metabolic reprogramming: current progress and open challenges. Front. Cell Dev. Biol.10, 972020. 10.3389/fcell.2022.972020
21
OuyangJ.WangH.HuangJ. (2023). The role of lactate in cardiovascular diseases. Cell Commun. Signal21 (1), 317. 10.1186/s12964-023-01350-7
22
PiH.WangG.WangY.ZhangM.HeQ.ZhengX.et al (2024). Immunological perspectives on atherosclerotic plaque formation and progression. Front. Immunol.15, 1437821. 10.3389/fimmu.2024.1437821
23
RitchieM. E.PhipsonB.WuD.HuY.LawC. W.ShiW.et al (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res.43 (7), e47. 10.1093/nar/gkv007
24
RobinX.TurckN.HainardA.TibertiN.LisacekF.SanchezJ. C.et al (2011). pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinforma.12, 77. 10.1186/1471-2105-12-77
25
SaigusaR.WinkelsH.LeyK. (2020). T cell subsets and functions in atherosclerosis. Nat. Rev. Cardiol.17 (7), 387–401. 10.1038/s41569-020-0352-5
26
ShamakiG. R.MarksonF.Soji-AyoadeD.AgwuegboC. C.BamgboseM. O.TamunoinemiB. M. (2022). Peripheral artery disease: a comprehensive updated review. Curr. Probl. Cardiol.47 (11), 101082. 10.1016/j.cpcardiol.2021.101082
27
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res.13 (11), 2498–2504. 10.1101/gr.1239303
28
SkoufosG.KakoulidisP.TastsoglouS.ZacharopoulouE.KotsiraV.MiliotisM.et al (2024). TarBase-v9.0 extends experimentally supported miRNA-gene interactions to cell-types and virally encoded miRNAs. Nucleic Acids Res.52 (D1), D304–d310. 10.1093/nar/gkad1071
29
SteenmanM.EspitiaO.MaurelB.GuyomarchB.HeymannM. F.PistoriusM. A.et al (2018). Identification of genomic differences among peripheral arterial beds in atherosclerotic and healthy arteries. Sci. Rep.8 (1), 3940. 10.1038/s41598-018-22292-y
30
SubramanianA.TamayoP.MoothaV. K.MukherjeeS.EbertB. L.GilletteM. A.et al (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U. S. A.102 (43), 15545–15550. 10.1073/pnas.0506580102
31
VergalloR.CreaF. (2020). Atherosclerotic plaque healing. N. Engl. J. Med.383 (9), 846–857. 10.1056/NEJMra2000317
32
WangJ.YangP.YuT.GaoM.LiuD.ZhangJ.et al (2022a). Lactylation of PKM2 suppresses inflammatory metabolic adaptation in pro-inflammatory macrophages. Int. J. Biol. Sci.18 (16), 6210–6225. 10.7150/ijbs.75434
33
WangN.WangW.WangX.MangG.ChenJ.YanX.et al (2022b). Histone lactylation boosts reparative gene activation post-myocardial infarction. Circ. Res.131 (11), 893–908. 10.1161/circresaha.122.320488
34
XuR.YuanW.WangZ. (2023). Advances in glycolysis metabolism of atherosclerosis. J. Cardiovasc Transl. Res.16 (2), 476–490. 10.1007/s12265-022-10311-3
35
YuG.WangL. G.HanY.HeQ. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. Omics16 (5), 284–287. 10.1089/omi.2011.0118
36
ZhangY.JiangH.DongM.MinJ.HeX.TanY.et al (2024). Macrophage MCT4 inhibition activates reparative genes and protects from atherosclerosis by histone H3 lysine 18 lactylation. Cell Rep.43 (5), 114180. 10.1016/j.celrep.2024.114180
Summary
Keywords
atherosclerosis, lactylation, immune infiltration, risk prediction, biomarkers
Citation
Qian J and Zhou Q (2025) Role of lactylation and immune infiltration in atherosclerosis: novel insights from bioinformatics analyses. Front. Genet. 16:1520325. doi: 10.3389/fgene.2025.1520325
Received
31 October 2024
Accepted
13 February 2025
Published
03 April 2025
Volume
16 - 2025
Edited by
Juexin Wang, Indiana University, Purdue University Indianapolis, United States
Reviewed by
Archana Prabahar, Cleveland State University, United States
Xin Shi, Shanghai Jiao Tong University, China
Updates
Copyright
© 2025 Qian and Zhou.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qing Zhou, zhouqing_penn@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.