ORIGINAL RESEARCH article

Front. Genet., 30 May 2022

Sec. Genetics of Aging

Volume 13 - 2022 | https://doi.org/10.3389/fgene.2022.865827

The Inflamm-Aging Model Identifies Key Risk Factors in Atherosclerosis

  • 1. Department of Biomedical Engineering, School of Intelligent Sciences, China Medical University, Shenyang, China

  • 2. Tumor Etiology and Screening Department of Cancer Institute and General Surgery, The First Affiliated Hospital of China Medical University, Shenyang, China

Article metrics

View details

15

Citations

5,4k

Views

2,1k

Downloads

Abstract

Background: Atherosclerosis, one of the main threats to human life and health, is driven by abnormal inflammation (i.e., chronic inflammation or oxidative stress) during accelerated aging. Many studies have shown that inflamm-aging exerts a significant impact on the occurrence of atherosclerosis, particularly by inducing an immune homeostasis imbalance. However, the potential mechanism by which inflamm-aging induces atherosclerosis needs to be studied more thoroughly, and there is currently a lack of powerful prediction models.

Methods: First, an improved inflamm-aging prediction model was constructed by integrating aging, inflammation, and disease markers with the help of machine learning methods; then, inflamm-aging scores were calculated. In addition, the causal relationship between aging and disease was identified using Mendelian randomization. A series of risk factors were also identified by causal analysis, sensitivity analysis, and network analysis.

Results: Our results revealed an accelerated inflamm-aging pattern in atherosclerosis and suggested a causal relationship between inflamm-aging and atherosclerosis. Mechanisms involving inflammation, nutritional balance, vascular homeostasis, and oxidative stress were found to be driving factors of atherosclerosis in the context of inflamm-aging.

Conclusion: In summary, we developed a model integrating crucial risk factors in inflamm-aging and atherosclerosis. Our computation pipeline could be used to explore potential mechanisms of related diseases.

Introduction

The burden of cardiovascular disease (CVD) has increased, leading to it being considered one of the most expensive diseases (i.e., seriously costing the economic burden) (Lopez et al., 2021). CVD is becoming an increasingly serious problem, with a high incidence rate and accounting for nearly one-third of total deaths (Stewart et al., 2017). Atherosclerosis is considered to be the major cause of CVD (Pahwa and Jialal, 2021). A chronic disease, atherosclerosis, involves the accumulation of plaques, and its major traditional risk factors (i.e., blood lipid proteasemia, diabetes, smoking, hypertension, and genetic abnormalities) have been reported by many epidemiological studies (Bergheanu et al., 2017; Pahwa and Jialal, 2021). However, atherosclerosis is a complex disease, and there are multiple cellular and molecular interactions involved in the progression of atherosclerosis that need to be further investigated and integrated (Lusis, 2000). Moreover, treatment is often difficult, expensive, invasive, and risky, and there is still no effective means of preventing atherosclerosis (Fogoros and Ali, 2020).

An increasing number of reports show that aging is a driving factor of atherosclerosis (El Assar et al., 2012; Wang and Bennett, 2012). Aging is closely related to endothelial dysfunction and arterial stiffness, which are considered to be early events leading to CVD. The aging process involves promotion of a series of risk factors (i.e., oxidative stress, endothelial dysfunction, and pro-inflammatory cytokines), leading to endothelial dysfunction and vascular system damage. In addition, cellular aging induces the release of microbubbles, further contributing to the development and calcification of atherosclerotic plaques (Alique et al., 2018). Therefore, the incidence of atherosclerosis increases with chronological age, and accelerated aging is the main risk factor for the development of atherosclerotic plaques (Kahn, 2021). Furthermore, atherosclerotic plaques represent a key index of cellular aging, which is characterized by reduced cell proliferation, increased DNA damage, and telomere shortening. In summary, a growing body of evidence indicates that atherosclerosis is promoted by cellular aging (Wang and Bennett, 2012).

Atherosclerosis is closely related to inflamm-aging, along with oxidative stress, endothelial dysfunction, and inflammation (Alique et al., 2018). Inflamm-aging is defined as a chronic inflammatory process during aging and mainly characterized by chronic progressive strengthening of the pro-inflammatory response (Xia et al., 2016). In other words, inflamm-aging promotes the body’s pro-inflammatory status with advancing aging, which is closely related to many aging diseases (Xia et al., 2016). A series of studies have shown that aging can promote atherosclerosis by damaging the connection between mitochondrial function and the intravascular inflammatory pathway (Tyrrell and Goldstein, 2020). For example, chronic inflammation is the main cause of age-related atherosclerosis, possibly exerting its effect through the IL-6 signaling pathway (Tyrrell and Goldstein, 2020). Furthermore, inflammatory factors released by senescent cells result in a senescence-associated secretory phenotype and even atherosclerosis (Childs et al., 2017). Selective targeting and elimination of these senescent cells (called heterolysis) has been shown to slow the growth of atherosclerotic lesions by reducing the release of inflammatory and adhesion factors (Childs et al., 2017).

Inflamm-aging induces an immune homeostasis imbalance during atherosclerosis development. As a chronic inflammatory disease, atherosclerosis is characterized by endothelial dysfunction and abnormal immune responses; immune cells and lipids are also involved (Karunakaran et al., 2015). The normal function of the adaptive immune system appears to decline with chronological age, and an advanced inflammatory response is promoted, leading to the imbalance of immune homeostasis (Sadighi Akha, 2018). Both innate and adaptive immune responses have been identified in atherosclerosis, and components of cholesterol-carrying low-density lipoprotein can further trigger abnormal inflammation (Keesey and Powley, 2008). Unsaturated fatty acids can also regulate both the differentiation and activation of T cells and then induce inflammation and oxidative stress (Hubler and Kennedy, 2016), and accumulation of dysregulated proteins can damage the normal function of the endoplasmic reticulum, which has an important role in oxidative stress and the immune response (So, 2018). In short, the nutritional balance is critical to immune homeostasis (Chandra, 1997). Chronic oxidative stress is also vital to immune/inflammatory cells during the inflamm-aging process (Fuente and Miquel, 2009). However, further and more systematic investigation is required to fully elucidate how atherosclerosis is triggered by inflamm-aging in the context of immune homeostasis imbalance. In order to study the potential mechanisms involved in inflamm-aging, immune homeostasis imbalance, and atherosclerosis, we hypothesized that the progression of atherosclerosis would be promoted by accelerated chronic inflammation during the aging process when immune homeostasis was disrupted (Figure 1A).

FIGURE 1

FIGURE 1

(A) Diagram of the hypothetical mechanism of atherosclerosis. (B) Workflow of our study.

At present, there are many omics profiles to be further analyzed; thus, proper prediction models need to be developed (Munger et al., 2021). Machine learning (ML), which is regarded as an extension of classical statistical modeling, can digest large amounts of data to identify high-order correlations and generate predictions (Alber et al., 2019; Zhavoronkov et al., 2019). For example, the supervised ML method can be used to understand the aging process (Fabris et al., 2017). Recently, an ML-derived score has been reported to show satisfactory performance in the prediction of cardiovascular events (Sánchez-Cabo et al., 2020). Some studies have also shown that ML can predict CVD and identify key markers related to inflamm-aging and CVD (Alalawi and Alsuwat, 2021; Munger et al., 2021; Sayed et al., 2021). Network analysis can also be used to find clusters of highly related genes and to associate modules with each other and with external sample traits (Deshpande et al., 2016). Recently, causality analysis was applied to CVD and found to be informative, using Mendelian randomization (MR) to integrate genomic and phenotypic information (Yazdani et al., 2016). Although many studies have identified CVD-related risk markers (Shadrina et al., 2020), the inflamm-aging mechanism in atherosclerosis (thus the casual relationship between inflamm-aging and atherosclerosis as well as relative potential risk factors) needs to be explored more thoroughly. Therefore, there is an urgent need to model atherosclerosis progression based on the inflamm-aging process, integrating the potential mechanisms related to immune homeostasis imbalance at a system level.

In order to further explore the inflamm-aging mechanism in atherosclerosis, we present a computational pipeline using ML and systems biology methods (Figure 1B), comprising the following steps (Lopez et al., 2021): we developed an improved model of inflamm-aging-related predictors (the inflamm-aging model) of atherosclerosis and identified relative risk markers (Stewart et al., 2017). An inflamm-aging score was calculated by integrating both aging and inflammatory markers (Pahwa and Jialal, 2021). MR was carried out to explore the causal relationship between inflamm-aging and atherosclerosis (Bergheanu et al., 2017). An inflamm-aging differential network was constructed by summarizing the (partial) correlation between each pair of genes based on the inflamm-aging score (Lusis, 2000). The Markov chain Monte Carlo (MCMC) method was used to further explore a sensitivity index for inflamm-aging, immune homeostasis imbalance, and atherosclerosis (Fogoros and Ali, 2020). Enrichment analysis and network analysis were carried out to study potential risk factors in atherosclerosis. As a result, the casual relationship between inflamm-aging and atherosclerosis was identified, and potential mechanisms related to nutritional balance, oxidative stress, and vessel homeostasis were integrated based on our computational pipeline.

Results

Modeling Inflamm-Aging Predictors and Identifying Relative Risk Markers

Gene expression data were downloaded from the Gene Expression Omnibus (GEO) database, including 926 samples and 11,313 genes (Supplementary Tables S1, S2). The RelifF algorithm was used to sort relative genes, and the k-nearest neighbors (kNN; k = 5 with the cosine distance) algorithm was selected to mode aging predictors and disease predictors (determined by 10-fold cross-validation; results in Supplementary Tables S3–S7). Other classification results are shown in Supplementary Tables S8–S10 based on the support vector machine (SVM), naive Bayes (NB), and ensemble learning (integrating 100 decision tree models) algorithms (Supplementary Tables S8–S10). The accuracy of the models was then validated using test data (Figure 2; Supplementary Figure S1; Table 1) including 126 atherosclerosis samples and 170 controls, where the aging/inflammatory/disease as well as relative parameters (selected by the cross-validation based on the training set, in Supplementary Tables S5) were also used in the test set. Furthermore, top-ranked aging, inflammatory, and disease markers were used in the improved inflamm-aging predictor (Supplementary Tables S5). The aging predictor accuracy was 0.7017, and 441 aging markers were identified; the prediction accuracy of the improved inflamm-aging model was 0.7037, with 186 inflammatory markers and 316 disease markers identified (see Materials and Methods; relative parameters in Supplementary Tables S5). The area under the receiver operating characteristic (ROC) curve values were 0.85682 and 0.74206 for the aging predictor and the improved inflamm-aging predictor, respectively. These results showed that our predictors were reliable and able to identify healthy young and old samples as well as distinguish between normal and disease samples with sufficient accuracy (Figure 2; Table 1).

FIGURE 2

FIGURE 2

Machine learning results. (A,B) Aging predictor from our previous study, selecting the number of aging markers. (C,D) Improved inflamm-aging predictor, selecting the number of disease markers. (A,C) Learning curve for the training dataset. (B,D) ROC curve for the test dataset.

TABLE 1

Accuracy on the training datasetAccuracy on the test dataset
Traditional aging predictor0.78810.6278
Aging predictor from our previous study0.83100.7017
Traditional disease predictor0.80960.6699
Disease predictor from our previous study0.80180.6801
Improved inflamm-aging predictor0.79890.7037

Accuracy of the predictor on the training and test datasets.

Various biomarkers of aging/disease have been associated with biological functions (Supplementary Table S5). For example, MMP8, the top aging marker, is a collagenase with proteolytic activity against matrix proteins and is an important participant in the function of the vascular system and the development of atherosclerosis (Laxton et al., 2009; Toba et al., 2017; Yu et al., 2019). Its expression in the heart increases with chronological age, enabling it to regulate signal transduction during the aging process by regulating the expression and activity of cytokines, chemokines, growth factors, hormones, and angiogenic factors (Mota et al., 2018). The top disease marker is KIR3DL1. Killer immunoglobulin-like receptors (KIR) encode receptors expressed on some subsets of natural killer cells and T lymphocytes (Niepiekło-Miniewska et al., 2014), which can regulate inflammatory response (Ross, 2014; Pugh et al., 2019). Among the 186 inflammatory markers identified, HSPA6 had the lowest false discovery rate (FDR = 4.38e-9). HSPA6 participates in cell protection (Noonan et al., 2007) and cell proliferation and affects the receptor binding affinity of IgG (Takai, 2002; Wang et al., 2020a). HSPA6 is also a putative target of endothelial nitric oxide synthase and is involved in the highly clinically relevant inhibition of vascular smooth muscle cell (VSMC) proliferation, which can repair vessels in human atherosclerosis (McCullagh et al., 2016). Thus, it could be speculated that VSMCs in the vascular system are closely related to atherosclerosis development in the context of inflamm-aging, consistent with previous research results (Bennett et al., 2016).

Comparison of Inflamm-Aging Scores Between Normal and Disease Samples

In order to investigate the accelerated inflamm-aging pattern in atherosclerosis, an inflamm-aging score was constructed by integrating 186 inflammatory markers and 441 aging markers based on the training data (details are given in the “Materials and Methods” section); the results are shown in Table 2. Both the median and mean of the inflamm-aging scores of atherosclerosis samples were greater than those of control samples. These results indicated an accelerated inflamm-aging pattern in atherosclerosis, consistent with previous findings (Lavin Plaza et al., 2020). To further verify this conclusion, the Kruskal–Wallis test (Figure 3A) was used to validate the significance of inflamm-aging in atherosclerosis. The results confirmed that atherosclerosis samples showed a significantly accelerated inflamm-aging pattern compared with normal samples (p = 1.6856e-05). In addition, when inflamm-aging scores were adjusted by chronological age (see Materials and Methods), the accelerated inflamm-aging pattern remained significant (p = 0.0024; Figure 3B; Table 2). A principal component analysis (PCA) plot is also shown in Supplementary Figure S2.

TABLE 2

Mean (control)Mean (disease)Median (control)Median (disease)
Original score0.01290.02300.01310.0250
Score adjusted by age0.00960.01610.00920.0162

Inflamm-aging scores of the disease and control groups.

FIGURE 3

FIGURE 3

Accelerated inflamm-aging pattern using the Kruskal–Wallis test. (A) Original inflamm-aging score. (B) Adjusted inflamm-aging score.

Exploring the Causal Relationship Between Inflamm-Aging and Atherosclerosis by MR

In order to study the inflamm-aging mechanisms in atherosclerosis more thoroughly, the MR method was used to distinguish causal relationships among aging, inflammation, and atherosclerosis, and partial correlation coefficients were calculated (see Materials and Methods for details). Crucial causal relationships related to inflamm-aging were identified in this way (Table 3, Supplementary Table S1, S2). For example, RPRM was the top aging marker in MR, related to the most (Takai, 2002) inflamm-aging markers. RPRM is a highly glycosylated cytoplasmic protein encoded by a p53-dependent putative tumor suppressor gene (Ohki et al., 2000). Its overexpression decreases cell proliferation, migration, and invasion and increases apoptosis (Figueroa et al., 2017). PECAM1 was the top disease marker, related to the maximum number (Merl-Pham et al., 2019) of aging markers. The expression of PECAM1 is associated with the response of various immune cell types (i.e., CD4 + T cells, B cells, CD8 + T cells, neutrophils, macrophages, and dendritic cells) as well as endothelial cells (Harry et al., 2008; Ye et al., 2021). The activation of endothelial PECAM1 has been shown to lead to nuclear translocation of nuclear factor κB (NF-κB) and expression of inflammatory and adhesion mediators, which promoted atherosclerosis (Harry et al., 2008). Polymorphism of PECAM1 is associated with the incidence rate of coronary atherosclerosis (Harry et al., 2008), indicating that inflamm-aging could accelerate atherosclerosis progression.

TABLE 3

Aging markerDisease markerDifference
BIRC2PLOD10.0650
CRBNIRF50.0617
ICAM2IRF50.0581
BIRC2CPNE10.0524
RPRMCEBPA0.0518
ICAM2NGB0.0514
ICAM2CEBPA0.0499
ICAM2PARP160.0496
ICAM2PECAM10.0491
ATP6V0CPOM1210.0483

Top 10 pairs with the largest absolute difference in correlation and partial correlation between “inflammation+ atherosclerosis−” and “inflammation+ atherosclerosis+”.

BIRC2–PLOD1 was the top aging–atherosclerosis pair (with the maximum differential K-S value of 0.06498; Table 3). BIRC2 is an E3 ubiquitin-protein ligase that regulates NF-κB signaling to maintain endothelial cell survival and blood vessel homeostasis during vascular development (Santoro et al., 2007; Samanta et al., 2020). PLOD1 can promote cell growth and the aerobic glycolysis process (Zhang et al., 2021). It also plays vital roles in atherosclerosis, including enhancing collagen fibril deposition (Yuan et al., 2021), modifying collagen (Wågsäter et al., 2013), and regulating the expression of extracellular matrix (ECM) (Merl-Pham et al., 2019). Thus, we conjectured that the vascular system would have an effect on atherosclerosis. In summary, key causal relationships between inflamm-aging and atherosclerosis were identified using MR.

The Inflamm-Aging Differential Network Revealed Key Inflamm-Aging Mechanisms in Atherosclerosis

In order to explore potential inflamm-aging mechanisms in atherosclerosis at a system level, an inflamm-aging differential network was constructed by calculating the (partial) correlation coefficient for each pair of genes in the training data and test data, respectively (details shown in Materials and Methods). The Pearson correlation coefficient of degree and degree frequency after logarithmic transformation was -0.8827 (p = 2.7871e-30; Supplementary Figure S3), with a small ratio of genes having a large degree. That is, the network had the scale-free characteristic. Furthermore, when the network was verified on the test set, the p-value was close to 0 using Fisher’s exact test. All these results indicate the reliability of the network.

The network markers indicated that endothelium, oxidative stress, and inflammation have important roles in atherosclerosis; the top-ranked genes (by degree) and their corresponding functions are shown in Table 4 (Yang et al., 2000; Zhang, 2013; Schisler et al., 2015; Mezzadra et al., 2017; Zheng et al., 2017; Guan et al., 2018; Lan et al., 2020; Nie et al., 2020; Zheng et al., 2020; Chen et al., 2021a; Cao et al., 2021; Covarrubias et al., 2021; Liu et al., 2021; Wang et al., 2021; Zheng et al., 2021). For example, the gene with the largest degree was PXN, which encodes a cytoskeleton protein. PXN is related to the attachment of actin membrane to ECM, tissue remodeling, and cell proliferation and survival. PXN also has important roles in focal adhesion, endothelial dysfunction, inflammation, and oxidative stress (Chen et al., 2021a).

TABLE 4

Gene symbolDegreeFunctionReference
PXN1171) Related to the attachment of actin membrane to ECMChen et al. (2021a)
2) Participates in tissue remodeling, cell proliferation, and survival
3) Important role in focal adhesion, endothelial dysfunction, inflammation, and oxidative stress
ZNF22110Induces cell apoptosisSchisler et al. (2015), Nie et al. (2020)
RBM101001) Regulates cell apoptosis, cell proliferation, cell invasion and metastasis, and inflammatory responseCao et al. (2021)
2) Affects the pathogenesis of atherosclerosis
NADK931) Important role in regulating cell aging and aging-related diseasesZhang (2013), Covarrubias et al. (2021)
2) Responsible for the production of NADP in the cytoplasmic matrix and participates in the counteraction of oxidative damage
ANKRD1290Participates in bone marrow mesenchymal stem cell differentiation, including cardiac muscle cells, nerve cells, blood cells, and myogenic cellsWang et al. (2021)
PJA289Promotes M1 macrophage polarization, M2 to M1 macrophage transformation, and the inflammatory responseZheng et al. (2021)
ATP5I89Induces oxidative stress and DNA damageLiu et al. (2021)
MAP3K3881) Necessary for angiogenesisYang et al. (2000), Zheng et al. (2017)
2) Related to endothelial cell proliferation and apoptosis and interacts with heart and myocardium
3) Participates in cell proliferation, differentiation, and apoptosis
4) Related to cellular aging
CMTM6841) Important role in regulating immune response and inflammatory activationMezzadra et al. (2017), Guan et al. (2018), Zheng et al. (2020)
2) Regulates the expression of PD-L1. PD-L1 can inhibit the function of T-cell inhibition
3) Promotes cell migration and invasion
SNRPG821) Arrests the cell cycleLan et al. (2020)
2) Activates p53 signaling

Top 10 markers with the highest degrees.

Sensitivity Analysis Further Highlighted the Immune Homeostasis Imbalance in Atherosclerosis

In order to further investigate the mechanisms underlying inflamm-aging, immune homeostasis imbalances, and atherosclerosis, the MCMC method was used to develop a sensitivity index for the relationship between inflamm-aging and atherosclerosis (based on the results obtained by MR). A bipartite graph of aging–disease pairs was constructed (see Materials and Methods, and Supplementary Figure S4).

The top 10 most sensitive pairs are shown in Table 5 (based on absolute differential frequencies). The relationship pair with the greatest difference (0.08172) was CRBN–MRPL40 (Table 5). CRBN functions as a potential regulator of cell homeostasis by mediating the ubiquitination of target substrates involved in processes including ion transport, the AMP-activated protein kinase (AMPK) signaling pathway, and cellular metabolisms (Shi and Chen, 2017; Jeon et al., 2021). It also protects cells from fat metabolism disorders (induced by high-fat levels) and negatively regulates CD4+ T cell activation (Shi and Chen, 2017), and has been implicated in a variety of diseases including CVD, obesity, and fatty liver (Shi and Chen, 2017). MRPL40 is related to mitochondrial functions (Li et al., 2019) and is crucial to ensuring ribosome translational fidelity and subsequent assembly of the oxidative phosphorylation complex (Jia et al., 2009). In short, these results indicate that oxidative stress plays an important part in the development of atherosclerosis.

TABLE 5

Aging markerDisease markerDifference
CRBNMRPL400.0817
RDXINPP5E0.0774
ENOPH1IRF50.0744
TACC2CEBPA0.0720
CKAP4CEBPA0.0678
MMP11ARPC20.0678
MCL1RBBP50.0645
HERC6UBE2G20.0645
PURAUBE2G20.0635
ATP6V0CPOM1210.0614

Top 10 pairs with the largest absolute difference between the frequency of disease inflammation and the frequency of the healthy inflammation in the relationship pairs with the shortest paths.

The aging marker related to the most (Keesey and Powley, 2008) disease markers (Supplementary Tables S13) was HEY2, the main sensor of Notch signal transduction. HEY2 has a key role in cardiac development; for example, it protects the heart from age-induced hypertrophic changes (Röning et al., 2022). It also regulates the proliferation of VASCs (Shirvani et al., 2007). PECAM1 was the inflamm-aging marker related to the maximum number (Alber et al., 2019) of aging markers (Supplementary Tables S14) in both the sensitivity analysis and MR (Harry et al., 2008; Ye et al., 2021). These results highlighted the critical role of the immune homeostasis imbalance in atherosclerosis progression.

Underlying Inflamm-Aging Mechanisms Based on Enrichment Analysis

In order to further explore the potential mechanism linking inflamm-aging and atherosclerosis, the shortest path of each aging–disease pair (as identified by MR or MCMC) in the inflamm-aging differential network was determined using the Dijkstra algorithm; then, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene ontology (GO) enrichment analysis with respect to biological process (BP) was performed on these shortest paths.

The top 10 KEGG pathways (those enriched in the shortest paths) are shown in Table 6 (Turi et al., 2019; Jia et al., 2021; Chothani et al., 2019; Yin et al., 2021; Manhas and Rath, 2020; Li et al., 2018; Nakajima et al., 2019; Hosseini et al., 2021; Critchley et al., 2018; Lathe et al., 2014; Margiotta, 2021; Nickel et al., 2002; Chen et al., 2020; Luo et al., 2020; Miller et al., 2021; Stadtman et al., 2005; Bonifácio et al., 2021; Blachier et al., 2020; Kumar et al., 2017; Récher, 2021; Klemm and Ruland, 2006). The most frequently enriched KEGG pathway was “ribosome” (enriched in eight shortest paths; Figure 4A). The rate of ribosome biogenesis has an important role in cell cycle progression and proliferation, which are closely related to coronary restenosis and atherosclerosis (Turi et al., 2019; Jia et al., 2021). Further, RNA translation and RNA transcription in the ribosome participate in fibration, which is vital to the development of atherosclerosis (Chothani et al., 2019). Advanced inflamm-aging impairs ribosome biogenesis, leading to autoantibody production against ribosomal proteins (Yin et al., 2021). Increased oxidative stress may also damage ribosomal RNA (Manhas and Rath, 2020). In short, the ribosome has an important role in the development of atherosclerosis.

TABLE 6

KEGGEnriched shortest pathScoreFunctionReference
Ribosome87.49491) Important role in cell cycle progression and proliferation, closely related to coronary restenosis and atherosclerosisChothani et al. (2019), Turi et al. (2019), Manhas and Rath (2020), Jia et al. (2021), Yin et al. (2021)
2) Regulates fibrosis, which is an important factor in the occurrence of atherosclerosis
3) With increasing age, abnormality of the immune system will affect the biogenesis of ribosomes and lead to the production of ribosomal protein autoantibodies
4) Increases oxidative stress during aging, eventually causing damage to ribosomal RNA
Parkinson’s disease76.5642Cardiovascular disease and Parkinson’s disease share common risk factorsLi et al. (2018)
Oxidative phosphorylation54.72131) Regulates adiponectin secretion in epicardial adipose tissue and has anti-atherosclerotic and anti-inflammatory effects on blood vesselsNakajima et al. (2019), Hosseini et al. (2021)
2) Triggers atherosclerosis by affecting mitochondrial functions
Huntington’s disease43.7964Myocardial dysfunction and vasoconstrictor dysfunction are involved in Huntington’s disease progressionCritchley et al. (2018)
Alzheimer’s disease32.8729Both atherosclerosis and Alzheimer’s disease are involved in inflammation, macrophage infiltration, and vascular system obstructionLathe et al. (2014)
SNARE interactions in vesicular transport10.99141) SNARE protein is an essential component that allows membrane fusion and can regulate vesicle fusionNickel et al. (2002), Chen et al. (2020), Margiotta (2021)
2) SNARE has a key role in cell homeostasis. Vesicle transport is the main mechanism of protein and lipid exchange between membrane-bound organelles in eukaryotic cells
3) Vesicle transport in mitochondria affects mitochondrial function, which is one of the characteristics of aging
Cytosolic DNA-sensing pathway10.98431) Promotes the expression of immune genesLuo et al. (2020), Miller et al. (2021)
2) Cytoplasmic DNA and cytoplasmic DNA-sensing adapter STING play a key part in aortic degeneration by promoting smooth muscle cell damage, macrophage MMP production, and ECM degeneration
3) Endogenous cytoplasmic DNA is a contributing factor to inflammation in the absence of pathogenic infection, which is associated with many chronic age-related diseases, including cancer, cardiovascular disease, and neurodegenerative diseases
Cysteine and methionine metabolism10.97851) The accumulation of oxidized methionine residues of protein is related to agingStadtman et al. (2005), Kumar et al. (2017), Blachier et al. (2020), Bonifácio et al. (2021)
2) Hyperhomocysteinemia and increased circulating levels of homocysteine (Hcy) are generally considered to be independent risk factors for peripheral atherosclerosis
3) Hcy can cause oxidative stress, inflammation, and endothelial dysfunction
4) Cysteine is an important source of energy and biomass and has a central role in cell metabolism
5) The metabolites of cysteine and methionine are related to fat metabolism and cardiovascular disease
Acute myeloid leukemia10.9752Inflammation can affect the progression of myeloid leukemiaRécher (2021)
FC epsilon RI signaling pathway10.9692Can induce the production of pro-inflammatory factors and mast cell degranulationKlemm and Ruland (2006)

Top 10 KEGG with the most numerous enriched paths.

FIGURE 4

FIGURE 4

Shortest paths for enrichment analysis. (A)a Maximum number of enriched paths for the KEGG pathway (“ribosome”); (B)KEGG pathway with the minimum FDR (“SNARE interactions in vesicular transport”); (C) maximum number of enriched paths for BP terms (“aminoglycoside antibiotic metabolic process” (GO:0030647)); (D) BP term with the minimum FDR (“vesicle-mediated transport between endosomal compartments” (GO:009,892)); (E) KEGG pathway with the minimum FDR based on sensitivity analysis (“cytosolic DNA-sensing pathway”); (F) BP term with the minimum FDR based on sensitivity analysis (“positive regulation of reactive oxygen species biosynthetic process” (GO:1903428)). Orange nodes, aging biomarkers; blue nodes, genes connecting aging biomarkers and atherosclerosis biomarkers; green nodes, atherosclerosis biomarkers; genes in the red square frames, genes with enriched functions.

The top enriched KEGG pathway (with the minimum FDR) was “SNARE interactions in vesicular transport” (p = 4.6346 and FDR = 0.0086; Figure 4B). The SNARE protein is an essential component in the regulation of vesicle fusion and plays a key part in maintaining cell homeostasis (Margiotta, 2021). Vesicle transport is the main mechanism of protein and lipid exchange between membrane-bound organelles in eukaryotic cells (Nickel et al., 2002); vesicle transport in mitochondria also affects mitochondrial function, which is one of the key characteristics of aging (Chen et al., 2020). Extracellular vesicles are involved in various processes including inflammation, coagulation, vascular dysfunction, angiogenesis, and aging, which contribute to the occurrence and progression of atherosclerotic thrombotic diseases (Konkoth et al., 2021).

The top 10 BP terms (those enriched in the most paths) are shown in Table 7 (Krause et al., 2016; Sainlos et al., 2003; Gonzalez and Spencer, 1998; Swirski and Nahrendorf, 2014; Meng et al., 2022; De Stefano et al., 2021; Kurek et al., 2021; Kiechl et al., 1998; Tokarev et al., 2009; Hu et al., 2015; Van-Assche et al., 2011; Shimizu et al., 1997; Moore et al., 2013; Patel et al., 2000). “Aminoglycoside antibiotic metabolic process” (GO:0030647) was the most enriched BP term (enriched in two shortest paths; Figure 4C). Aminoglycosides provide a favorable scaffold for the synthesis of various cationic lipids (Krause et al., 2016; Sainlos et al., 2003), and are used to treat and prevent endocarditis (Gonzalez and Spencer, 1998). “Vesicle-mediated transport between endosomal compartments” (GO:0098927) was the BP term with minimum FDR (p = 9.0282e-08, FDR = 6.7540e-04; Figure 4D). The main communication process between a cell and its environment is vesicle transport, through bidirectional trafficking among the Golgi, endosomes, and lysosomes (Tokarev et al., 2009). The endosomal–lysosomal system and endomembrane system are also closely related to both lipid and protein metabolism (Hu et al., 2015).

TABLE 7

BPEnriched shortest pathScoreFunctionReference
Aminoglycoside antibiotic metabolic process 003064721.91601) Aminoglycoside antibiotic works by inhibiting protein synthesisGonzalez and Spencer (1998), Sainlos et al. (2003), Krause et al. (2016)
2) Aminoglycosides provide a favorable scaffold for the synthesis of various cationic lipids
3) Aminoglycosides are used to treat and prevent endocarditis
Myeloid leukocyte activation GO:000227421.87431) Leukocytes can engender protective immunity and protect the host from damageSwirski and Nahrendorf (2014)
2) Important role in inflammation, immunity, and atherosclerosis
Glycoside metabolic process GO:001613721.86821) Glycosides are associated with inflammation, oxidative stress, lipid metabolism, and atherosclerosisDe Stefano et al. (2021), Kurek et al. (2021), Meng et al. (2022)
2) Flavonoid glycosides can protect vascular endothelial cells by inhibiting inflammation to restrict atherosclerosis
3) Luteolin-7-o-glucoside can reduce the activity of oxidative stress and inflammatory mechanisms in different physiological systems
Tertiary alcohol metabolic process GO:190264421.8638Alcohol affects the progress of atherosclerosisKiechl et al. (1998)
Vesicle-mediated transport between endosomal compartments GO:009892710.99931) The main communication process between the cell and its environmentTokarev et al. (2009), Hu et al. (2015)
2) Bi-directional trafficking between the Golgi, endosomes, and lysosomes connects the two major intracellular trafficking pathways
3) Endosomes gradually mature into lysosomes, and the acidification of late endosomes is accompanied by vesicle transport
4) Endosomal–lysosomal system and endomembrane system are related to lipid homeostasis and protein homeostasis
Positive regulation of reactive oxygen species biosynthetic process GO:190342810.99841) Induces LDL oxidation and foam cell formation and activates many redox-sensitive transcription factors including NF-κB and AP1Van-Assche et al. (2011)
2) Regulates the expression of multiple promoters/anti-inflammatory genes involved in atherosclerosis
Positive regulation of nitric oxide metabolic process GO:190440710.9984Nitric oxide (NO) is associated with atherosclerosis as it inhibits the proliferation of VSMCs, monocyte/macrophage adhesion, platelet aggregation and adhesion, and LDL oxidationShimizu et al. (1997)
Regulation of nitric oxide metabolic process GO:008016410.9982NO is associated with atherosclerosis as it inhibits the proliferation VSMCs, monocyte/macrophage adhesion, platelet aggregation and adhesion, and LDL oxidationShimizu et al. (1997)
Regulation of macrophage activation GO:004303010.99701) The number and phenotype of macrophages can affect the inflammatory state of plaquesMoore et al. (2013)
2) Lipoprotein uptake by macrophages is considered to be one of the earliest pathogenic events in new plaques
3) The M2 polarization pathway of macrophages can prevent atherosclerosis
Reactive nitrogen species metabolic process GO:200105710.99691) Affects endothelial functionPatel et al. (2000)
2) Low-level proteins and lipids modified by reactive nitrogen can promote the development of atherosclerosis through mechanisms involving signal transduction

Top 10 BP with the most numerous enriched paths.

Enrichment analysis was also performed based on the results of the sensitivity analysis. “Cytosolic DNA-sensing pathway” was the KEGG pathway with the minimum FDR in the sensitivity analysis (p = 8.4266e-05, FDR = 0.0157; Figure 4E); this can promote the expression of immune markers (Luo et al., 2020). Cytoplasmic DNA and cytoplasmic DNA-sensing adapter STING play key parts in aortic degeneration by promoting smooth muscle cell damage, matrix metalloproteinase (MMP) production in macrophages, and ECM degeneration (Luo et al., 2020). Endogenous cytoplasmic DNA is also a contributing factor to inflammation in the absence of pathogenic infection and is associated with many chronic age-related diseases, including CVD and neurodegenerative diseases (Miller et al., 2021). “Positive regulation of reactive oxygen species biosynthetic process” was the BP term with minimum FDR (GO:1903428, p = 4.3964e-07, FDR = 0.0016; Figure 4F). Reactive oxygen species (ROS) play a critical part in the regulation of the expression of multiple promoters and anti-inflammatory genes involved in atherosclerosis, via inducing low-density lipoprotein (LDL) oxidation, foam cell formation, and activation of many redox-sensitive transcription factors including NF-κB and activating protein 1 (Van-Assche et al., 2011).

In summary, a series of risk factors in atherosclerosis were identified, including inflammation, nutritional metabolism, energy metabolism, and the vascular system.

Network Markers Revealed Crucial Mechanisms Between Aging and Atherosclerosis

Network markers were also identified by calculating the betweenness of each aging–disease shortest path; the top 10 markers are shown in Table 8. The network marker with the greatest betweenness (Munger et al., 2021) was RPL35 (Figure 5A; Table 8). RPL35 is an important component of the 60S ribosomal subunit, with a key role in mRNA translation and protein synthesis (Wu et al., 2021). Gene betweenness was also calculated according to the sensitivity analysis. The top markers and their betweenness values are shown in Table 9. According to this analysis, IRAK1 and VAMP8 were the top network markers (betweenness = 8; Figures 5B, C; Table 9). IRAK1 is a threonine/serine kinase that is related to the pathogenesis of atherosclerosis (Rana et al., 2016). It regulates lipid accumulation and prevents the formation of macrophage foam cells (Rana et al., 2016), as well as participating in multiple IL-1- and Toll-like receptors (TLR) driven signaling processes and in the regulation of levels of immune pro-inflammatory cytokines (Singer et al., 2018). VAMP8 is relevant to macrophage degranulation and secretion of TNF-α (Pushparaj et al., 2009), and regulates lysosome–autophagosome interaction and exocytosis (Jones et al., 2012; Chen et al., 2021b). Its phosphorylation regulates lipid metabolism in the liver (Huang et al., 2021). VAMP8 recruits NOX2 into phagosomes, causing lipid oxidation and membrane damage, and also mediates the transport of NOX2 to endosomes and phagosomes, thereby promoting the induction of the cytolytic T cell immune response (Dingjan et al., 2017). In short, the identification of network markers showed that the occurrence of atherosclerosis was closely related to inflammation and immunity homeostasis imbalance.

TABLE 8

Gene symbolBetweennessp-value
RPL35210
IRAK1200
VAMP8190
LY86170
MS4A3170
IMP4170
GDI1160
PJA2140
CLN3140
C6orf48140

Top 10 genes with the highest betweenness in the aging acceleration network.

FIGURE 5

FIGURE 5

Network marker with the highest betweenness. (A) Based on MR; (B,C) based on MCMC.

TABLE 9

Gene symbolBetweennessp-value
IRAK180
VAMP880
RPL3560
RPL1860
PJA260
PARP660
JTB50
SARS50
HDGF50
PPP1R12A50

Top 10 genes with the highest betweenness in the aging acceleration network based on sensitivity analysis.

Discussion

CVD has become a serious health problem, the incidence of which increases with age. It is well known that both aging and inflammation have vital roles in atherosclerosis, which is one of the key risk factors for CVD. However, the exact relationships between inflamm-aging and atherosclerosis and the role of immune homeostasis have to be studied more systematically. It is particularly important to elucidate the causal relationship between inflamm-aging and atherosclerosis, in order to determine the mechanisms underlying atherosclerosis and identify potential therapeutic targets. A causal analysis was used to identify the relationship between aging and atherosclerosis using MR. An ML approach was used owing to its great advantages in processing large amounts of data and generating prediction models with appropriate accuracies. Network analysis and sensitivity analysis were also performed to integrate multiple factors and identify relative risk markers, which will be conducive to more systematic study of atherosclerosis.

Our results showed that immune homeostasis is vital for the development of atherosclerosis. Disorders of the immune system can lead to atherosclerosis, and several characteristics related to immune homeostasis imbalance were identified in our study. For example, “myeloid leukocyte activation” (GO:0002274) and “regulation macrophage activation” (GO:0043030) involve the effects of inflammation on atherosclerosis. Several risk markers related to the immune response that have key roles in atherosclerosis were identified, including the top aging marker MMP8, disease markers KIR3DL1 and MAGEA6, network markers based on degree (RBM10, PJA2, and CMTM6; Table 4), and network markers based on betweenness (IRAK1 and VAMP8). Other potential markers including HEY2 (identified by sensitivity analysis) and PECAM1 (identified by MR and sensitivity analysis) and the top causal pairs BIRC2–PLOD1 and CRBN–MRPL40 also suggested a role of immunity in atherosclerosis. In brief, the abnormal immune response is a key characteristic of atherosclerosis (Karunakaran et al., 2015); our results further confirmed that immune homeostasis imbalance is crucial in atherosclerosis in the context of inflamm-aging.

Vascular homeostasis imbalance was also identified as a factor in atherosclerosis progression. It has been reported that dysfunction of endothelium and VSMCs, an early marker of atherosclerosis, can disrupt the balance between vasoconstriction and vasodilation, thereby promoting the inflammatory response and exacerbating atherosclerosis (Nkoenig et al., 2022; Alberts et al., 2002; Davignon and Ganz, 2004). The top inflammation marker (HSPA6) and the top network markers according to a degree (PXN and MAP3K3) in our results indicated an effect of endothelial cells and VSMCs. Dysfunction of endothelium contributes to the recruitment of inflammatory cells, and VASC proliferation is related to upregulation of pro-inflammatory cytokines (Alique et al., 2018). Consequently, there is a key link between inflammation and vascular homeostasis. The top enriched KEGG pathway was the “cytosolic DNA-sensing pathway” (Figure 4E), indicating the potential relevance of inflammation and vascular homeostasis. In addition, PECAM1 enables endothelial cells to respond to fluid shear stress and regulates inflammatory signaling pathways (Privratsky and Newman, 2014). In the causal pair BIRC2–PLOD1, BIRC2 promotes endothelial cell survival and blood vascular homeostasis via regulation of inflammatory/immune signaling pathways (Santoro et al., 2007; Samanta et al., 2020). PLOD1 encodes lysyl hydroxylase-1, which is essential for ECM maturation. ECM proteolytic fragments affect a variety of functions and properties of inflammatory and immune cells (Adairkirk and Senior, 2008); ECM was also relevant to the integrity and stability of VSMCs (Nkoenig et al., 2022). It has been reported that inflamm-aging can lead to changes in vascular structure and aging of VSMCs (Chi et al., 2019). In short, vascular homeostasis has a key role in atherosclerosis and inflamm-aging.

The role of oxidative stress was further investigated. Oxidative stress increases ROS generation and/or reduces the imbalance of the body’s innate antioxidant defense system and is thus a key influence on atherosclerosis (Kattoor et al., 2017). Our results provided further evidence of the association of oxidative stress with atherosclerosis, for instance, the top enriched BPs were “positive regulation of reactive oxygen species biosynthetic process” (GO:190,342) and “oxidative phosphorylation” (Table 4). Oxidative stress caused by excessive production of ROS has become the driven risk factor in atherosclerosis (Kattoor et al., 2017). The network marker (based on degree; Table 4) NADK is related to the production of NADPH. NADPH oxidase is the main source of ROS (Ganguly et al., 2021). ROS also has important roles in inflammatory response, apoptosis, cell growth, and changes in vascular tone as well as in the oxidation of LDL cholesterol (Kattoor et al., 2017). For example, ROS can cause physical damage to plaques via stimulating smooth muscle cell migration and collagen deposition, leading to endothelial dysfunction (Ganguly et al., 2021). In summary, the interaction between oxidative stress and vascular homeostasis is vital to atherosclerosis. The relationship between oxidative stress and aging was also explored. Oxidative stress is well known to have an important influence on aging. One of the aging markers identified by our MR approach was RPRM, a p53-inducible gene (Xu et al., 2012). TP53 regulates the expression of various genes involved in the maintenance of homeostasis, including cell cycle regulation and redox homeostasis (i.e., production of antioxidant enzymes) (Gambino et al., 2013; Beyfuss and Hood, 2018). Overall, these findings indicate that oxidative stress is a key factor in atherosclerosis in the context of inflamm-aging.

Nutritional imbalance (i.e., lipid and protein imbalance) has a vital role in atherosclerosis. Abnormal nutrient metabolism can affect the occurrence and development of atherosclerosis, and it has been suggested that atherosclerosis could be prevented and treated by adjusting the proportion of nutrients in the diet (Wei et al., 2021). For example, vesicles represent the main transport pathway for nutritional macromolecules (i.e., lipids and protein) and thus have a vital role in the establishment of the nutritional imbalance that is related to the development of atherosclerosis. The terms identified in our enrichment analyses, including “vesicle-mediated transport between endosomal compartments” (GO:0098927; Table 7) and “SNARE interactions in vesicular transport” (Table 6), confirmed this viewpoint. Furthermore, lipids have an important role in maintaining nutritional balance. Atherosclerosis is considered to be an inflammatory disease with chronic lipid accumulations and poor adaptability of arterial intima (So, 2018). Our sensitivity analysis identified CRBN, a negative regulator of AMPK in vivo, which is related to the metabolism of glucose and fat and is considered to be a regulator of body metabolism and energy homeostasis (Lee et al., 2013). We also identified IRAK1 and VAMP8 (using betweenness), which are related to the metabolism of lipids (Pushparaj et al., 2009; Singer et al., 2018). There is substantial evidence that lipids are related to atherosclerosis progression, as a high cholesterol level is a major risk factor in atherosclerosis and CVD, and fatty acid metabolism is closely related to endothelial function and protein homeostasis (Malekmohammad et al., 2021; Graham et al., 2021; Hasan and Fischer, 2021). Protein homeostasis is also relevant to macrophages in atherosclerosis (Koga et al., 2011). Our results demonstrated the role of protein as a nutrient in atherosclerosis. For instance, the top network marker (based on betweenness), RPL35, has a role in protein translation and endoplasmic reticulum docking (Wu et al., 2021). In addition, “aminoglycoside antibiotic metabolic process” (GO:0030647; Table 7) and “glycoside metabolic process” (GO:0016137; Table 7) further indicated a key role related to both lipids and proteins in atherosclerosis progression. As a result, the nutritional imbalance was identified as one of the driving factors behind atherosclerosis.

In the theory of inflamm-aging, the aging process leads to a chronic progressive inflammatory state, involving both the activation of innate immunity and accelerated pro-inflammatory response. As a result, the dynamic balance between pro- and anti-inflammatory cytokines is disrupted, and atherosclerosis was induced by the resulting immune homeostasis imbalance (Kahn, 2021). Our results further confirmed the key role of inflamm-aging in atherosclerosis and also identified dysfunction of endothelium and vessels, oxidative stress, and nutritional imbalance as risk factors. Studies have shown interactions among inflammation, nutrition, and endothelial and vascular homeostasis as well as oxidative stress. For instance, lipid accumulation is involved in the production of pro-inflammatory cytokines (Malekmohammad et al., 2021) and inflammation can affect lipid metabolism in macrophages (Osada et al., 2017). Moreover, oxidized LDL (ox-LDL) can remain in blood vessels, inducing dysfunction of the endothelium, which further responds to inflammatory signals, with a crucial role in atherosclerosis development (Malekmohammad et al., 2021). In addition, increased ROS production promotes endothelial dysfunction, leading to remodeling, platelet aggregation, loss of vasodilation, inflammation, and smooth muscle cell growth (Higashi et al., 2014). In our analysis, the “ribosome” pathway was identified as the most enriched KEGG pathway (Figure 4A), which highlighted the interactions among nutrition, vessels, and oxidative stress during the aging process (Stein et al., 2022). In summary, our work integrated potential mechanisms involving the inflammatory response, nutrition imbalance, and vascular homeostasis imbalance as well as oxidative stress, consistent with and contributing further to the inflamm-aging theory (Figure 6).

FIGURE 6

FIGURE 6

Mechanism of atherosclerosis induced by inflamm-aging. Gray genes, aging makers; blue genes, disease markers; green genes, inflammation markers; purple genes, markers with high degrees; orange genes, network nodes with high betweenness. Green arrow, inflammation; blue arrow, immune homeostasis; orange arrow, oxidative stress; gray arrow, vascular homeostasis; yellow arrow, nutritional balance.

Conclusion

In this study, an ML method was used to model predictors of aging and disease and identify relative risk markers. The accelerated inflammatory pattern in atherosclerosis was verified by comparing inflamm-aging scores. The causal relationship between inflamm-aging and atherosclerosis was explored through MR, and then through sensitivity analysis, enrichment analysis, and network analysis; the results indicated that disorders of lipid homeostasis and dysfunction of endothelial cells and vascular homeostasis could promote the development of atherosclerosis. The role of oxidative stress was also confirmed. In summary, our work revealed that in the context of inflamm-aging, nutrition imbalance and vascular homeostasis, as well as oxidative stress, are risk factors that promote the development of atherosclerosis in a coordinated manner.

Materials and Methods

Data and Pre-Processing

All gene expression data were downloaded from the GEO database (

https://www.ncbi.nlm.nih.gov/geo/

) in NCBI, including the GSE7638, GSE9820, GSE12288, GSE20681, GSE22253, GSE37171, and GSE123696 datasets. These datasets were from seven different platforms: GPL571, GPL6255, GPL96, GPL4133, GPL6244, GPL570, and GPL15207. The pre-processing steps for obtaining gene expression profiles were as follows:

  • 1) Age and disease information of each sample (disease was encoded as 1, healthy as 0) was obtained, and samples without age were deleted.

  • 2) The gene expression matrix for each dataset was integrated by summarizing the probe number within the gene symbol.

  • 3) The total data matrix was integrated, and the missing gene expression values were filled with values of 0.

  • 4) Genes with zero expression values ≥ 30% were deleted.

  • 5) Gene expression profiles were transformed by the logarithmic transformation method in containing outliers (>3 times of standard deviation or <−3 times of standard deviation).

  • 6) Based on the mean and the standard deviation of gene expression for healthy aged individuals, z-score normalization was performed for both normal aged samples and atherosclerosis samples.

  • 7) The singular value decomposition method was used to eliminate inter-sample variation based on the top three principal components in healthy aged individuals.

  • 8) The z-score was then used to normalize all individuals based on the mean and standard deviation of healthy aged individuals.

As a result, a total of 926 samples were obtained, including 313 samples from healthy aged individuals (age >50 years, 210 training datasets +103 test datasets), 207 samples from healthy young individuals (age ≤50, 140 + 67), 340 samples from aged atherosclerosis patients (age >50, 235 + 105), and 66 samples from young atherosclerosis patients (age ≤50, 45 + 21), with a total of 11,313 gene symbols.

Modeling the Improved Inflamm-Aging Predictor of Atherosclerosis

Gene expression profiles were divided into a training dataset and test dataset at a ratio close to 2:1, and the inflamm-aging prediction model was established by integrating aging, inflammation, and disease (atherosclerosis) markers. Aging and disease models were first constructed, respectively; then, an improved inflamm-aging model was developed by including inflamm-aging markers and adjusting the disease predictor (

Supplementary Figure S5

). The following steps were performed for the inflamm-aging model:

  • 1) To summarize the interactions among key markers according to aging/disease, the gene expression profiles were transformed (or replaced) based on the results of our previous study (Wang et al., 2020b). The Pearson correlation coefficient was used to evaluate relevance and redundancy:

where the phenotype was set as 0 (control) versus 1 (disease) in the disease model, or as the transformation of the age in the aging model:

The interaction is summarized as follows:

where only the first log2 (

n

+1) interactions for each gene were included.

  • 2) The 11,313 gene profiles (training dataset) were sorted by the ReliefF algorithm.

  • 3) The kNN (k = 5 with the cosine distance), SVM, NB, and ensemble learning (established 100 decision tree models) algorithms were used to select key markers (in the aging or disease model) based on the training dataset and to evaluate the performance on the test dataset, with the help of 10-fold cross-validation. The kNN (k = 5 with the cosine distance) algorithm was selected, and 441 aging markers and 390 candidate disease markers were identified in this step (Supplementary Tables S3, S4).

  • 4) The inflamm-aging model was constructed by integrating the aging, inflammation, and disease markers as follows:

    • ① 568 candidate inflammatory markers were obtained based on the GO BPs using the gene set enrichment analysis (GSEA) platform (http://software.broadinstitute.org/gsea/downloads.jsp, version 7.4, keyword “INFLAMMATORY”); then, aging markers or candidate disease markers overlapping with candidate inflammatory markers were removed.

    • ② Aging markers and candidate disease markers were further filtered using the Kruskal–Wallis test (p < 0.05 and FDR < 0.2).

    • ③ Inflamm-aging markers were identified by transforming selected inflammatory markers combined with the aging markers [also based on a previous study (Wang et al., 2020b) as in step (Lopez et al., 2021)], where an interaction was retained if p < 0.05 and FDR < 0.2 using the Kruskal–Wallis test.

    • ④ Disease markers were also transformed based on the inflammatory markers.

    • ⑤ The kNN method, ReliefF algorithm, and 10-fold cross-validation were used to reselect the optimized integrative disease markers, where some of the top candidate 500 disease markers were replaced by disease markers related to inflamm-aging in step ④. Therefore, the improved inflamm-aging predictor contained 116 aging markers, 186 inflammatory markers, and 316 disease markers (Supplementary Table S5).

  • 5) The traditional aging or disease predictor using original gene datasets was used as a comparator to evaluate the prediction ability of the improved inflamm-aging model (also using kNN, ReliefF, and 10-fold cross-validation with the same details; results in Supplementary Table S6, S7).

  • 6) Other classifiers including SVM, NB, and the ensemble algorithm were also used for comparison of their prediction abilities.

Calculating the Inflamm-Aging Score

The steps to calculate the inflamm-aging score for each sample were as follows:

  • 1) The distance between healthy young and old neighbors was calculated using the original inflamm-aging score based on the inflammatory markers in the inflamm-aging model; the cosine distance was also used:

  • 2) The original score was adjusted by age:

where b is the regression coefficient for an aging score, and the transformed age is calculated as follows:

  • 3) Both the original and adjusted inflamm-aging scores were used to test the accelerated pattern in atherosclerosis based on the Kruskal–Wallis test.

  • 4) The atherosclerosis (disease) score was calculated based on the identified disease markers (with the same details as for the inflamm-aging score) in the improved inflamm-aging model for further (sensitivity/network) analysis.

Causality Analysis Based on MR

MR is a statistical method used to evaluate the causal relationship between risk factors and results based on observed data (

Burgess et al., 2019

;

Burgess et al., 2020

). The causal relationships among the instrumental variable, the risk factor, and the outcome variable were evaluated as follows.

  • 1) There was a correlation between the instrumental variable and the risk factor.

  • 2) There was no correlation between the instrumental variable and the confounding factor.

  • 3) There was no correlation between the instrumental variable and the outcome variable after deleting the effect from the risk factor.

In this work, aging markers were used as instrumental variables, inflammatory markers were used as risk factors, and inflamm-aging markers were used as outcome variables to explore the causal relationship between inflamm-aging and atherosclerosis. The steps were as follows:

  • 1) The training dataset was divided into three groups: low inflammation in health (“inflammation− atherosclerosis-”), high inflammation in health (“inflammation+ atherosclerosis−“), and low inflammation in disease (“inflammation+ atherosclerosis+”).

  • 2) For each group, the Pearson correlation coefficient between each aging marker and the inflamm-aging score was calculated to evaluate the correlation between the instrument variable and the risk factor. The relationship was retained if the correlation for “inflammation− atherosclerosis− was between those of “inflammation+ atherosclerosis−” and “inflammation+ atherosclerosis+” (meaning that the atherosclerosis was triggered by the inflammation).

  • 3) To filter out the correlation between the instrumental variable and confounding factors, a permutation test was performed by generating the simulated inflamm-aging score from the same number of randomly selected markers; this process was repeated 1,000 times, then the p-value was calculated as the proportion of occurrence times (larger than the real difference) of the absolute difference between “inflammation+ atherosclerosis−” and “inflammation+ atherosclerosis+” in 1,000 permutations. The relationship between each aging marker and the inflamm-aging score was retained if the permutation p < 0.05.

  • 4) If the gene was within both the aging marker set and the disease marker set, then it was removed.

  • 5) For each group, the Pearson correlation coefficient between each aging marker and disease marker was calculated to evaluate the correlation between the instrument variable and the outcome variable; and the partial Pearson correlation coefficient between each aging marker and disease marker based on the inflamm-aging score was calculated to evaluate the correlation between the instrument variable and the outcome variable without the background of the risk factor. The causal relationship between each aging marker and disease marker was retained if the difference between the correlation and the partial correlation for “inflammation− atherosclerosis− was between those for “inflammation+ atherosclerosis−” and “inflammation+ atherosclerosis+” (meaning that the atherosclerosis was triggered by the inflammation).

  • 6) To filter out the potential correlation between confounding factors and the outcome variable, a permutation test was performed by generating the simulated aging marker and disease marker from the same number of randomly selected markers, repeated 1,000 times; then, the permutation P-value was calculated using the same details as in step (Pahwa and Jialal, 2021).

  • 7) To filter out potential weak instrumental variables, the aging markers were further analyzed by evaluating the correlation between the inflamm-aging score and each aging marker, using the result of the point multiplication between the inflamm-aging score and each aging marker; then, the residual of the inflamm-aging score was tested between the atherosclerosis and control subgroup, using the Kruskal–Wallis test (p < 0.05 and FDR < 0.2).

where

b

is the regression coefficient.

  • 8) To filter out the effect of horizontal pleiotropy, the aging–disease causal relationship was further examined by comparing the correlation between each aging and disease marker, through the inflamm-aging score or otherwise. Steps ①–③ were used to calculate the correlations between instrumental variables and outcome variables without the background of the risk factor, and step ④ was used to calculate the correlations between instrumental variables and outcome variables without the context of the risk factor.

① The residual of each disease marker (“residual A”) is calculated based on the inflamm-aging score:where b1 is the regression coefficient.

② The residual of each aging marker (“residual B”) is calculated based on the inflamm-aging score:where b2 is the regression coefficient.

③ The aforementioned two residuals were further compared, and the residual of the disease marker is calculated (as “residual C”):where b3 is the regression coefficient.

④ The residual of the disease marker (“residual D”) was calculated based on the aging marker.

⑤ The difference (between “residual C” and “residual D”) was tested between the atherosclerosis and control subgroups, using the Kruskal–Wallis test (P < 0.05 and FDR<0.2).

Finally, the causal relationship between the aging marker and disease marker (through the inflamm-aging score) was retained. Thus, 1,340 aging–disease pairs were identified as causal relationships, including 116 aging markers and 312 inflamm-aging markers.

Constructing the Inflamm-Aging Differential Network

In order to further reveal the relationship between inflamm-aging and atherosclerosis, an inflamm-aging differential network was constructed by the following steps:

  • 1) To compare the relationship of each gene pair in the context of inflamm-aging, both the Pearson correlation coefficient for each pair of genes and the partial correlation coefficient based on the inflamm-aging score were calculated based on the atherosclerosis and control groups, respectively.

  • 2) The Benjamini–Hochberg FDR method was used to adjust the P-values of the correlation coefficient and the partial correlation coefficient.

  • 3) The differences in the correlations and partial correlations between the atherosclerosis and control groups were calculated.

  • 4) The edge between each gene pair was retained if the signs of the difference in values for the correlation and partial correlation were different, and values of p < 0.05 and FDR<0.1 were obtained in step (Pahwa and Jialal, 2021).

  • 5) The scale-free characteristics of the inflamm-aging acceleration differential networks were verified by the power-law distribution.

  • 6) The shortest path between each pair of aging and inflamm-aging markers was selected based on each aging acceleration differential network using the Djikstra algorithm.

  • 7) The network constructed based on the training data was used for further analysis (identifying shortest paths, exploring biological functions between aging and atherosclerosis, etc.), and the network constructed based on the test data was used to validate the training network.

Global Sensitivity Analysis Using the MCMC Method

In order to further explore the immune homeostasis imbalance in atherosclerosis in the context of inflamm-aging, an overall sensitivity analysis based on the MCMC method was carried out, where the aging–disease pairs identified by MR were used as candidate relationships to be further evaluated. The MCMC method is used for sampling from certain posterior distributions following a given probabilistic background in high-dimensional space. The key step in MCMC is to construct a Markov chain whose equilibrium distribution equals the target probability distribution. The steps are as follows.

  • 1) Construct a transition kernel of an ergodic Markov chain. In this study, the prior distribution for each parameter was the normal distribution based on all the identified aging inflammation and inflamm-aging markers for each group (“inflammation− atherosclerosis−”, “inflammation+ atherosclerosis−, and “inflammation+ atherosclerosis+”).

  • 2) Simulate the chain until it reaches equilibrium. The Metropolis–Hastings sampling method is used to determine whether the new sample (θ*) is acceptable based on the α value:

where

p

(

θn|X

) and

p

(

θ*|X

) are the posterior probabilities of the

n

th accepted sample and the new sample,

q

(

θn→θ*

) is the transition probability from the

n

th accepted sample to the new sample, and

q

(

θ*→θn

) is the transition probability from the new sample to the

n

th accepted sample.

In this work, the inflamm-aging score and disease score were used to evaluate the simulated samples, with 1,000 random samples used as candidate samples for each group.

  • 3) Perform global sensitivity analysis. In this study, the distribution (ratio) of each aging–disease pair is calculated as follows:

In addition, if the value for “inflammation− atherosclerosis− was between those of “inflammation+ atherosclerosis− and “inflammation+ atherosclerosis+”, then the aging–disease pair was identified as a sensitive relationship for further network analysis. As a result, 442 aging–disease pairs were identified as causal relationships, including 93 aging markers and 81 disease markers.

Enrichment Analysis

Gene functions with significant correlations were further explored through enrichment analysis of the shortest paths. GO terms and KEGG pathways were identified using the GSEA platform (http://software.broadinstitute.org/gsea/downloads.jsp, version 7.4). The hypergeometric distribution was used to test the enrichment of GO BP terms and KEGG pathways using the following formula:where N is the total number of genes in the gene set, M is the number of known genes (KEGG pathways or BP terms), n is the number of genes identified in each shortest path, and k is the number of common genes between known genes and candidate genes identified in each aging–disease shortest path. The p-value for each path was controlled using the Benjamini–Hochberg method. Finally, the paths with p < 0.05 and FDR<01 were retained. The enrichment score is calculated as follows:

Identifying Network Markers

The subnetwork with the shortest pathways among the selected aging–disease pairs was constructed (after MR or MCMC), and genes in the subnetwork were sorted by their betweenness in descending order. To test whether the top betweenness genes were hubs in the background network, we ran a permutation to count the occurrences of the top genes in the shortest paths between randomly selected genes (containing the same numbers of aging–disease pairs, 1,340 in MR or 442 in MCMC) when they had greater betweennesses than those in our study. We repeated this process 1,000 times, and the p-value was calculated as the proportion of occurrences of the top betweenness genes in 1,000 permutations.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

YH, YC, LY, JW, and YW performed the algorithm and analyzed the data; YH, YC, LY, and YW wrote the manuscript; XS and YW designed and sponsored the study. All authors read and approved the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (32000478 to YW) and the National Key R&D Program of China (2016YFC0901704, 2017YFA0505500, 2017YFC0907505, and 2017YFC0908405). The funders had no role in study design, data collection, and analysis, the decision to publish, or the preparation of the manuscript.

Acknowledgments

We thank Charlesworth Author Services (http://www.charlesworthauthorservices.com/∼Frontiers) for its linguistic assistance during the preparation of this manuscript.

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.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.865827/full#supplementary-material

Supplementary Figure S1

ML results for control methods. (A, B) Traditional aging predictor. (C,D) Traditional atherosclerosis predictor. (C,D) Disease predictor from our previous study. (A,C,E) Learning curve for the training dataset. (B,D,F) ROC curve for the test dataset.

Supplementary Figure S2

PCA plots of the inflamm-aging model based on (A) the inflamm-aging score and disease risk score; and (B) the first two principal components. Blue circle, normal test set; blue ‘+’, normal training set; red circle, atherosclerosis test set; red ‘+’, atherosclerosis training set.

Supplementary Figure S3

Degree distribution of the inflamm-aging differential network.

Supplementary Figure S4

Heatmap of the bipartite graph between aging markers and disease markers based on MCMC. The x-axis represents disease markers; the Y-axis represents aging markers.

Supplementary Figure S5

Workflow of the improved inflamm-aging model.

Supplementary Table S1

Detailed datasets used in this work.

Supplementary Table S2

Gene symbols used in this work.

Supplementary Table S3

Selected aging markers from our previous study, with their ReliefF weights.

Supplementary Table S4

Selected disease markers from our previous study, with their ReliefF weights.

Supplementary Table S5

Selected aging, inflammatory, and disease markers in the improved inflamm-aging model, with their P-values, FDRs, and ReliefF weights.

Supplementary Table S6

Selected aging markers obtained by the traditional method, with their ReliefF weights.

Supplementary Table S7

Selected disease markers obtained by the traditional method, with their ReliefF weights.

Abbreviations

AUC, area under the ROC curve; BP, biological process; CVD, cardiovascular disease; ECM, extracellular matrix; FDR, false discovery rate; GEO, Gene Expression Omnibus; GO, Gene Ontology; GSEA, gene set enrichment analysis; KEGG, Kyoto Encyclopedia of Genes and Genomes; KIR, killer immunoglobulin-like receptors; kNN, k-nearest neighbors; MCMC, Markov chain Monte Carlo; ML, machine learning; MMP, matrix metalloproteinase; MR, Mendelian randomization; NB, naive Bayes; OXPHOS, oxidative phosphorylation; PCA, principal components analysis; ROC, receiver operating characteristic; ROS, reactive oxygen species; SVM, support vector machine; TLR, toll-like receptors; VSMC, vascular smooth muscle cells.

References

  • 1

    AdairkirkT.SeniorR. (2008). Fragments of Extracellular Matrix as Mediators of Inflammation. Int. J. Biochem. Cell Biol.40 (6-7), 11011110. 10.1016/j.biocel.2007.12.005

  • 2

    AlalawiH. H.AlsuwatM. S. (2021). Detection of Cardiovascular Disease Using Machine Learning Classification Models. IJERT07, 14.

  • 3

    AlberM.Buganza TepoleA.CannonW. R.DeS.Dura-BernalS.GarikipatiK.et al (2019). Integrating Machine Learning and Multiscale Modeling-Perspectives, Challenges, and Opportunities in the Biological, Biomedical, and Behavioral Sciences. npj Digit. Med.2, 115. 10.1038/s41746-019-0193-y

  • 4

    AlbertsB.JohnsonA.LewisJ.RaffM.RobertsK.WalterP. (2002). Molecular Biology of the Cell. 4th edition. New York: Garland Science.

  • 5

    AliqueM.Ramírez-CarracedoR.BodegaG.CarracedoJ.RamírezR. (2018). Senescent Microvesicles: A Novel Advance in Molecular Mechanisms of Atherosclerotic Calcification. Int. J. Mol. Sci.19 (7), 2003. 10.3390/ijms19072003

  • 6

    BennettM. R.SinhaS.OwensG. K. (2016). Vascular Smooth Muscle Cells in Atherosclerosis. Circ. Res.118 (4), 692702. 10.1161/CIRCRESAHA.115.306361

  • 7

    BergheanuS. C.BoddeM. C.JukemaJ. W. (2017). Pathophysiology and Treatment of Atherosclerosis. Neth Heart J.25 (4), 231242. 10.1007/s12471-017-0959-2

  • 8

    BeyfussK.HoodD. A. (2018). A Systematic Review of P53 Regulation of Oxidative Stress in Skeletal Muscle. Redox Rep.23 (1), 100117. 10.1080/13510002.2017.1416773

  • 9

    BlachierF.AndriamihajaM.BlaisA. (2020). Sulfur-Containing Amino Acids and Lipid Metabolism. J. Nutr.150 (Suppl. 1), 2524S2531S. 10.1093/jn/nxaa243

  • 10

    BonifácioV. D. B.PereiraS. A.SerpaJ.VicenteJ. B. (2021). Cysteine Metabolic Circuitries: Druggable Targets in Cancer. Br. J. Cancer124, 862879. 10.1038/s41416-020-01156-1

  • 11

    BurgessS.Davey SmithG.DaviesN. M.DudbridgeF.GillD.GlymourM. M.et al (2019). Guidelines for Performing Mendelian Randomization Investigations. Wellcome Open Res.4, 186. 10.12688/wellcomeopenres.15555.1

  • 12

    BurgessS.FoleyC. N.AllaraE.StaleyJ. R.HowsonJ. M. M. (2020). A Robust and Efficient Method for Mendelian Randomization with Hundreds of Genetic Variants. Nat. Commun.11, 376. 10.1038/s41467-019-14156-4

  • 13

    CaoY.DiX.ZhangQ.LiR.WangK. (2021). RBM10 Regulates Tumor Apoptosis, Proliferation, and Metastasis. Front. Oncol.11, 603932. 10.3389/fonc.2021.603932

  • 14

    ChandraR. K. (1997). Nutrition and the Immune System: an Introduction. Am. J. Clin. Nutr.66 (2), 460S463S. 10.1093/ajcn/66.2.460S

  • 15

    ChenP.-L.HuangK.-T.ChengC.-Y.LiJ.-C.ChanH.-Y.LinT.-Y.et al (2020). Vesicular Transport Mediates the Uptake of Cytoplasmic Proteins into Mitochondria in Drosophila melanogaster. Nat. Commun.11, 2592. 10.1038/s41467-020-16335-0

  • 16

    ChenQ.HaoM.WangL.LiL.ChenY.ShaoX.et al (2021). Prefused Lysosomes Cluster on Autophagosomes Regulated by VAMP8. Cell Death Dis.12, 939. 10.1038/s41419-021-04243-0

  • 17

    ChenY.ZhaoH.XiaoY.ShenP.TanL.ZhangS.et al (2021). Pan-cancer Analysis Reveals an Immunological Role and Prognostic Potential of PXN in Human Cancer. Aging13 (12), 1624816266. 10.18632/aging.203154

  • 18

    ChiC.LiD. J.JiangY. J.TongJ.FuH.WuY. H.et al (2019). Vascular Smooth Muscle Cell Senescence and Age-Related Diseases: State of the Art. Biochim. Biophys. Acta Mol. Basis Dis.1865, 18101821. 10.1016/j.bbadis.2018.08.015

  • 19

    ChildsB. G.GluscevicM.BakerD. J.LabergeR.-M.MarquessD.DananbergJ.et al (2017). Senescent Cells: an Emerging Target for Diseases of Ageing. Nat. Rev. Drug. Discov.16, 718735. 10.1038/nrd.2017.116

  • 20

    ChothaniS.SchäferS.AdamiE.ViswanathanS.WidjajaA. A.LangleyS. R.et al (2019). Widespread Translational Control of Fibrosis in the Human Heart by RNA-Binding Proteins. Circulation140 (11), 937951. 10.1161/CIRCULATIONAHA.119.039596

  • 21

    CovarrubiasA. J.PerroneR.GrozioA.VerdinE. (2021). NAD+ Metabolism and its Roles in Cellular Processes during Ageing. Nat. Rev. Mol. Cell Biol.22 (2), 119141. 10.1038/s41580-020-00313-x

  • 22

    CritchleyB. J.IsalanM.MielcarekM. (2018). Neuro-Cardio Mechanisms in Huntington's Disease and Other Neurodegenerative Disorders. Front. Physiol.9. 10.3389/fphys.2018.00559

  • 23

    DavignonJ.GanzP. (2004). Role of Endothelial Dysfunction in Atherosclerosis. Circulation109 (23 Suppl. 1), III2732. 10.1161/01.CIR.0000131515.03336.f8

  • 24

    De StefanoA.CaporaliS.Di DanieleN.RovellaV.CardilloC.SchinzariF.et al (2021). Anti-Inflammatory and Proliferative Properties of Luteolin-7-O-Glucoside. Int. J. Mol. Sci.22 (3), 1321. 10.3390/ijms22031321

  • 25

    DeshpandeV.SharmaA.MukhopadhyayR.ThotaL. N. R.GhatgeM.VangalaR. K.et al (2016). Understanding the Progression of Atherosclerosis through Gene Profiling and Co-expression Network Analysis in Apob tm2Sgy Ldlr tm1Her Double Knockout Mice. Genomics107 (6), 239247. 10.1016/j.ygeno.2016.04.007

  • 26

    DingjanI.PaardekooperL. M.VerboogenD. R. J.von MollardG. F.Ter BeestM.van den BogaartG. (2017). VAMP8-mediated NOX2 Recruitment to Endosomes Is Necessary for Antigen Release. Eur. J. Cell Biol.96 (7), 705714. 10.1016/j.ejcb.2017.06.007

  • 27

    El AssarM.AnguloJ.VallejoS.PeiróC.Sánchez-FerrerC. F.Rodríguez-MañasL. (2012). Mechanisms Involved in the Aging-Induced Vascular Dysfunction. Front. Physio.3, 132. 10.3389/fphys.2012.00132

  • 28

    FabrisF.MagalhãesJ. P. d.FreitasA. A. (2017). A Review of Supervised Machine Learning Applied to Ageing Research. Biogerontology18 (2), 171188. 10.1007/s10522-017-9683-y

  • 29

    FigueroaR. J.Carrasco-AvinoG.WichmannI. A.LangeM.OwenG. I.SiekmannA. F.et al (2017). Reprimo Tissue-specific Expression Pattern Is Conserved between Zebrafish and Human. PLoS One12 (5), e0178274. 10.1371/journal.pone.0178274

  • 30

    FogorosR. N. (2020). “An Overview of Atherosclerosis,” in Verywell Health. Editor AliY. S..

  • 31

    FuenteM.MiquelJ. (2009). An Update of the Oxidation-Inflammation Theory of Aging: the Involvement of the Immune System in Oxi-Inflamm-Aging. Cpd15 (26), 30033026. 10.2174/138161209789058110

  • 32

    GambinoV.De MicheleG.VeneziaO.MigliaccioP.Dall'OlioV.BernardL.et al (2013). Oxidative Stress Activates a Specific P53 Transcriptional Response that Regulates Cellular Senescence and Aging. Aging Cell12, 34353445. 10.1111/acel.12060

  • 33

    GangulyU.KaurU.ChakrabartiS. S.SharmaP.AgrawalB. K.SasoL.et al (2021). Oxidative Stress, Neuroinflammation, and NADPH Oxidase: Implications in the Pathogenesis and Treatment of Alzheimer's Disease. Oxidative Med. Cell. Longev.2021, 119. 10.1155/2021/7086512

  • 34

    GonzalezL. S.SpencerJ. P. (1998). Aminoglycosides: A Practical Review. Am. Fam. Physician58 (8), 18111820.

  • 35

    GrahamS. E.ClarkeS. L.WuK. H. H.KanoniS.ZajacG. J. M.RamdasS.et al (2021). The Power of Genetic Diversity in Genome-wide Association Studies of Lipids. Nature600, 675. 10.1038/s41586-021-04064-3

  • 36

    GuanX.ZhangC.ZhaoJ.SunG.SongQ.JiaW. (2018). CMTM6 Overexpression Is Associated with Molecular and Clinical Characteristics of Malignancy and Predicts Poor Prognosis in Gliomas. EbioMedicine35, 233243. 10.1016/j.ebiom.2018.08.012

  • 37

    HarryB. L.SandersJ. M.FeaverR. E.LanseyM.DeemT. L.ZarbockA.et al (2008). Endothelial Cell PECAM-1 Promotes Atherosclerotic Lesions in Areas of Disturbed Flow in ApoE-Deficient Mice. Atvb28, 20032008. 10.1161/atvbaha.108.164707

  • 38

    HasanS. S.FischerA. (2021). The Endothelium: An Active Regulator of Lipid and Glucose Homeostasis. Trends Cell Biol.31 (1), 3749. 10.1016/j.tcb.2020.10.003

  • 39

    HigashiY.MaruhashiT.NomaK.KiharaY. (2014). Oxidative Stress and Endothelial Dysfunction: Clinical Evidence and Therapeutic Implications. Trends Cardiovasc Med.24 (4), 165169. 10.1016/j.tcm.2013.12.001

  • 40

    HosseiniZ.MarinelloM.DeckerC.SansburyB. E.SadhuS.GerlachB. D.et al (2021). Resolvin D1 Enhances Necroptotic Cell Clearance through Promoting Macrophage Fatty Acid Oxidation and Oxidative Phosphorylation. Atvb41, 10621075. 10.1161/atvbaha.120.315758

  • 41

    HuY.-B.DammerE. B.RenR.-J.WangG. (2015). The Endosomal-Lysosomal System: from Acidification and Cargo Sorting to Neurodegeneration. Transl. Neurodegener.4, 18. 10.1186/s40035-015-0041-1

  • 42

    HuangH.OuyangQ.ZhuM.YuH.MeiK.LiuR. (2021). mTOR-Mediated Phosphorylation of VAMP8 and SCFD1 Regulates Autophagosome Maturation. Nat. Commun.12, 6622. 10.1038/s41467-021-26824-5

  • 43

    HublerM. J.KennedyA. J. (2016). Role of Lipids in the Metabolism and Activation of Immune Cells. J. Nutr. Biochem.34, 17. 10.1016/j.jnutbio.2015.11.002

  • 44

    JeonS.YoonY.-S.KimH. K.HanJ.LeeK. M.SeolJ. E.et al (2021). Ablation of CRBN Induces Loss of Type I Collagen and SCH in Mouse Skin by Fibroblast Senescence via the P38 MAPK Pathway. Aging13 (5), 64066419. 10.18632/aging.202744

  • 45

    JiaF.WuQ.WangZ.ZhangM.YuanS.CheY.et al (2021). BOP1 Knockdown Attenuates Neointimal Hyperplasia by Activating P53 and Inhibiting Nascent Protein Synthesis. Oxidative Med. Cell. Longev.2021, 120. 10.1155/2021/5986260

  • 46

    JiaL.KaurJ.StuartR. A. (2009). Mapping of the Saccharomyces cerevisiae Oxa1-Mitochondrial Ribosome Interface and Identification of MrpL40, a Ribosomal Protein in Close Proximity to Oxa1 and Critical for Oxidative Phosphorylation Complex Assembly. Eukaryot. Cell8 (11), 17921802. 10.1128/ec.00219-09

  • 47

    JonesL. C.MoussaL.FulcherM. L.ZhuY.HudsonE. J.O'NealW. K.et al (2012). VAMP8 Is a Vesicle SNARE that Regulates Mucin Secretion in Airway Goblet Cells. J. Physiol.590 (Pt 3), 545562. 10.1113/jphysiol.2011.222091

  • 48

    KahnJ. (2021). Atherosclerosis and Aging — Insights into the Role of the Endothelial Glycocalyx in Cardiovascular Health. Available at: https://www.todaysgeriatricmedicine.com/news/ex_052419.shtml.

  • 49

    KarunakaranD.ThrushA. B.NguyenM.-A.RichardsL.GeoffrionM.SingaraveluR.et al (2015). Macrophage Mitochondrial Energy Status Regulates Cholesterol Efflux and Is Enhanced by Anti-miR33 in Atherosclerosis. Circ. Res.117 (3), 266278. 10.1161/CIRCRESAHA.117.305624

  • 50

    KattoorA. J.PothineniN. V. K.PalagiriD.MehtaJ. L. (2017). Oxidative Stress in Atherosclerosis. Curr. Atheroscler. Rep.19, 42. 10.1007/s11883-017-0678-6

  • 51

    KeeseyR. E.PowleyT. L. (2008). Body Energy Homeostasis. Appetite51 (3), 442445. 10.1016/j.appet.2008.06.009

  • 52

    KiechlS.WilleitJ.RunggerG.EggerG.OberhollenzerF.BonoraE. (1998). Alcohol Consumption and Atherosclerosis: What Is the Relation?Stroke29 (5), 900907. 10.1161/01.str.29.5.900

  • 53

    KlemmS.RulandJ. (2006). Inflammatory Signal Transduction from the FcεRI to NF-Κb. Immunobiology211 (10), 815820. 10.1016/j.imbio.2006.07.001

  • 54

    KogaH.KaushikS.CuervoA. M. (2011). Protein Homeostasis and Aging: the Importance of Exquisite Quality Control. Ageing Res. Rev.10 (2), 205215. 10.1016/j.arr.2010.02.001

  • 55

    KonkothA.SaraswatR.DubrouC.SabatierF.LeroyerA. S.LacroixR.et al (2021). Multifaceted Role of Extracellular Vesicles in Atherosclerosis. Atherosclerosis319, 121131. 10.1016/j.atherosclerosis.2020.11.006

  • 56

    KrauseK. M.SerioA. W.KaneT. R.ConnollyL. E. (2016). Aminoglycosides: An Overview. Cold Spring Harb. Perspect. Med.6 (6), a027029. 10.1101/cshperspect.a027029

  • 57

    KumarA.PalfreyH. A.PathakR.KadowitzP. J.GettysT. W.MurthyS. N. (2017). The Metabolism and Significance of Homocysteine in Nutrition and Health. Nutr. Metab. (Lond)14, 78. 10.1186/s12986-017-0233-z

  • 58

    KurekJ. M.KrólE.KrejpcioZ. (2021). Steviol Glycosides Supplementation Affects Lipid Metabolism in High-Fat Fed STZ-Induced Diabetic Rats. Nutrients13 (1), 112. 10.3390/nu13010112

  • 59

    LanY.LouJ.HuJ.YuZ.LyuW.ZhangB. (2020). Downregulation of SNRPG Induces Cell Cycle Arrest and Sensitizes Human Glioblastoma Cells to Temozolomide by Targeting Myc through a P53-dependent Signaling Pathway. Cancer Biol. Med.17 (1), 112131. 10.20892/j.issn.2095-3941.2019.0164

  • 60

    LatheR.SapronovaA.KotelevtsevY. (2014). Atherosclerosis and Alzheimer - Diseases with a Common Cause? Inflammation, Oxysterols, Vasculature. BMC Geriatr.14, 36. 10.1186/1471-2318-14-36

  • 61

    Lavin PlazaB.PhinikaridouA.AndiaM. E.PotterM.LorrioS.RashidI.et al (2020). Sustained Focal Vascular Inflammation Accelerates Atherosclerosis in Remote Arteries. Atvb40 (9), 21592170. 10.1161/atvbaha.120.314387

  • 62

    LaxtonR. C.HuY.DucheneJ.ZhangF.ZhangZ.LeungK.-Y.et al (2009). A Role of Matrix Metalloproteinase-8 in Atherosclerosis. Circulation Res.105 (9), 921929. 10.1161/CIRCRESAHA.109.200279

  • 63

    LeeK. M.YangS.-J.KimY. D.ChoiY. D.NamJ. H.ChoiC. S.et al (2013). Disruption of the Cereblon Gene Enhances Hepatic AMPK Activity and Prevents High-Fat Diet-Induced Obesity and Insulin Resistance in Mice. Diabetes62 (6), 18551864. 10.2337/db12-1030

  • 64

    LiJ.RyanS. K.DeboerE.CookK.FitzgeraldS.LachmanH. M.et al (2019). Mitochondrial Deficits in Human iPSC-Derived Neurons from Patients with 22q11.2 Deletion Syndrome and Schizophrenia. Transl. Psychiatry9, 302. 10.1038/s41398-019-0643-y

  • 65

    LiQ.WangC.TangH.ChenS.MaJ. (2018). Stroke and Coronary Artery Disease Are Associated with Parkinson's Disease. Can. J. Neurol. Sci.45 (5), 559565. 10.1017/cjn.2018.56

  • 66

    LiuP.ChengY.XuH.HuangH.TangS.SongF.et al (2021). TRESK Regulates Gm11874 to Induce Apoptosis of Spinal Cord Neurons via ATP5i Mediated Oxidative Stress and DNA Damage. Neurochem. Res.46 (8), 19701980. 10.1007/s11064-021-03318-w

  • 67

    LopezE. O.BallardB. D.JanA. (2021). Cardiovascular Disease. StatPearls [Internet].

  • 68

    LuoW.WangY.ZhangL.RenP.ZhangC.LiY.et al (2020). Critical Role of Cytosolic DNA and its Sensing Adaptor STING in Aortic Degeneration, Dissection, and Rupture. Circulation141 (1), 4266. 10.1161/CIRCULATIONAHA.119.041460

  • 69

    LusisA. J. (2000). Atherosclerosis. Nature.407, 233241. 10.1038/35025203

  • 70

    MalekmohammadK.BezsonovE. E.Rafieian-KopaeiM. (2021). Role of Lipid Accumulation and Inflammation in Atherosclerosis: Focus on Molecular and Cellular Mechanisms. Front. Cardiovasc. Med.8, 707529. 10.3389/fcvm.2021.707529

  • 71

    ManhasR.RathP. C. (2020). Ribosome, Protein Synthesis, and Aging. Models, Mol. Mech. Biogerontology2020, 6787. 10.1007/978-981-32-9005-1_4

  • 72

    MargiottaA. (2021). Role of SNAREs in Neurodegenerative Diseases. Cells10 (5), 991. 10.3390/cells10050991

  • 73

    McCullaghK. J.CooneyR.O'BrienT. (2016). Endothelial Nitric Oxide Synthase Induces Heat Shock Protein HSPA6 (HSP70B') in Human Arterial Smooth Muscle Cells. Nitric Oxide52, 4148. 10.1016/j.niox.2015.11.002

  • 74

    MengN.ChenK.WangY.HouJ.ChuW.XieS.et al (2022). Dihydrohomoplantagin and Homoplantaginin, Major Flavonoid Glycosides from Salvia Plebeia R. Br. Inhibit oxLDL-Induced Endothelial Cell Injury and Restrict Atherosclerosis via Activating Nrf2 Anti-oxidation Signal Pathway. Molecules27 (6), 1990. 10.3390/molecules27061990

  • 75

    Merl-PhamJ.BasakT.KnüppelL.RamanujamD.AthanasonM.BehrJ.et al (2019). Quantitative Proteomic Profiling of Extracellular Matrix and Site-specific Collagen Post-translational Modifications in an In Vitro Model of Lung Fibrosis. Matrix Biol. Plus1, 100005. 10.1016/j.mbplus.2019.04.002

  • 76

    MezzadraR.SunC.JaeL. T.Gomez-EerlandR.de VriesE.WuW.et al (2017). Identification of CMTM6 and CMTM4 as PD-L1 Protein Regulators. Nature549 (7670), 106110. 10.1038/nature23669

  • 77

    MillerK. N.VictorelliS. G.SalmonowiczH.DasguptaN.LiuT.PassosJ. F.et al (2021). Cytoplasmic DNA: Sources, Sensing, and Role in Aging and Disease. Cell184 (22), 55065526. 10.1016/j.cell.2021.09.034

  • 78

    MooreK. J.SheedyF. J.FisherE. A. (2013). Macrophages in Atherosclerosis: a Dynamic Balance. Nat. Rev. Immunol.13 (10), 709721. 10.1038/nri3520

  • 79

    MotaR.ParryT. L.YatesC. C.QiangZ.EatonS. C.MwizaJ. M.et al (2018). Increasing Cardiomyocyte Atrogin-1 Reduces Aging-Associated Fibrosis and Regulates Remodeling In Vivo. Am. J. Pathology188 (7), 16761692. 10.1016/j.ajpath.2018.04.007

  • 80

    MungerE.HickeyJ. W.DeyA. K.JafriM. S.KinserJ. M.MehtaN. N. (2021). Application of Machine Learning in Understanding Atherosclerosis: Emerging Insights. Apl. Bioeng.5 (1), 011505. 10.1063/5.0028986

  • 81

    NakajimaT.YokotaT.ShinguY.YamadaA.IbaY.UjihiraK.et al (2019). Impaired Mitochondrial Oxidative Phosphorylation Capacity in Epicardial Adipose Tissue Is Associated with Decreased Concentration of Adiponectin and Severity of Coronary Atherosclerosis. Sci. Rep.9, 3535. 10.1038/s41598-019-40419-7

  • 82

    NickelW.BrüggerB.WielandF. T. (2002). Vesicular Transport: the Core Machinery of COPI Recruitment and Budding. J. Cell Sci.115 (16), 32353240. 10.1242/jcs.115.16.3235

  • 83

    NieY.CaiS.YuanR.DingS.ZhangX.ChenL.et al (2020). Zfp422 Promotes Skeletal Muscle Differentiation by Regulating EphA7 to Induce Appropriate Myoblast Apoptosis. Cell Death Differ.27 (5), 16441659. 10.1038/s41418-019-0448-9

  • 84

    Niepiekło-MiniewskaW.ŻukN.DubisJ.KurpiszM.SenitzerD.HavrylyukA.et al (2014). Two New Cases of KIR3DP1, KIR2DL4-Negative Genotypes, One of Which Is Also Lacking KIR3DL2. Arch. Immunol. Ther. Exp.62 (5), 423429. 10.1007/s00005-014-0299-5

  • 85

    NkoenigS.CavusO.WilliamsJ.BernierM.TonnigesJ.SucharskiH.et al (2022). New Mechanistic Insights to PLOD1-Mediated Human Vascular Disease. Transl. Res.239, 1. 10.1016/j.trsl.2021.08.002

  • 86

    NoonanE. J.PlaceR. F.GiardinaC.HightowerL. E. (2007). Hsp70B′ Regulation and Function. Cell Stress Chaper12, 393402. 10.1379/csc-278e.1

  • 87

    OhkiR.NemotoJ.MurasawaH.OdaE.InazawaJ.TanakaN.et al (2000). Reprimo, a New Candidate Mediator of the P53-Mediated Cell Cycle Arrest at the G2 Phase. J. Biol. Chem.275 (30), 2262722630. 10.1074/jbc.C000235200

  • 88

    OsadaK.MaedaY.YoshinoT.NojimaD.BowlerC.TanakaT. (2017). Enhanced NADPH Production in the Pentose Phosphate Pathway Accelerates Lipid Accumulation in the Oleaginous Diatom Fistulifera Solaris. Algal Res.23, 126134. 10.1016/j.algal.2017.01.015

  • 89

    PahwaR.JialalI. (2021). Atherosclerosis. StatPearls [Internet].

  • 90

    PatelR. P.MoelleringD.Murphy-UllrichJ.JoH.BeckmanJ. S.Darley-UsmarV. M. (2000). Cell Signaling by Reactive Nitrogen and Oxygen Species in Atherosclerosis. Free Radic. Biol. Med.28 (12), 17801794. 10.1016/s0891-5849(00)00235-5

  • 91

    PrivratskyJ. R.NewmanP. J. (2014). PECAM-1: Regulator of Endothelial Junctional Integrity. Cell Tissue Res.355 (3), 607619. 10.1007/s00441-013-1779-3

  • 92

    PughJ.Nemat-GorganiN.DjaoudZ.GuethleinL. A.NormanP. J.ParhamP. (2019). In Vitro education of Human Natural Killer Cells by KIR3DL1. Life Sci. Alliance2 (6), e201900434. 10.26508/lsa.201900434

  • 93

    PushparajP. N.TayH. K.WangC.-C.HongW.MelendezA. J. (2009). VAMP8 Is Essential in Anaphylatoxin-Induced Degranulation, TNF-α Secretion, Peritonitis, and Systemic Inflammation. J. Immunol.183 (2), 14131418. 10.4049/jimmunol.0804061

  • 94

    RanaM.KumarA.TiwariR. L.SinghV.ChandraT.DikshitM.et al (2016). IRAK Regulates Macrophage Foam Cell Formation by Modulating Genes Involved in Cholesterol Uptake and Efflux. Bioessays38 (7), 591604. 10.1002/bies.201600085

  • 95

    RécherC. (2021). Clinical Implications of Inflammation in Acute Myeloid Leukemia. Front. Oncol.11, 623952. 10.3389/fonc.2021.623952

  • 96

    RöningT.MaggaJ.AnnaL.KerkeläR.KoivunenP.SerpiR. (2022). Activation of the Hypoxia Response Pathway Protects against Age-Induced Cardiac Hypertrophy. J. Mol. Cell Cardiol.164, 148155. 10.1016/j.yjmcc.2021.12.003

  • 97

    RossK. A. (2014). Coherent Somatic Mutation in Autoimmune Disease. PLoS ONE9 (7), e101093. 10.1371/journal.pone.0101093

  • 98

    Sadighi AkhaA. A. (2018). Aging and the Immune System: An Overview. J. Immunol. Methods463, 2126. 10.1016/j.jim.2018.08.005

  • 99

    SainlosM.BelmontP.VigneronJ.-P.LehnP.LehnJ.-M. (2003). Aminoglycoside-Derived Cationic Lipids for Gene Transfection: Synthesis of Kanamycin A Derivatives. Eur. J. Org. Chem.2003 (15), 27642774. 10.1002/ejoc.200300164

  • 100

    SamantaD.HuangT. Y.-T.ShahR.YangY.PanF.SemenzaG. L. (2020). BIRC2 Expression Impairs Anti-cancer Immunity and Immunotherapy Efficacy. Cell Rep.32 (8), 108073. 10.1016/j.celrep.2020.108073

  • 101

    Sánchez-CaboF.RosselloX.FusterV.BenitoF.ManzanoJ. P.SillaJ. C.et al (2020). Machine Learning Improves Cardiovascular Risk Definition for Young, Asymptomatic Individuals. J. Am. Coll. Cardiol.76 (14), 16741685. 10.1016/j.jacc.2020.08.017

  • 102

    SantoroM. M.SamuelT.MitchellT.ReedJ. C.StainierD. Y. R. (2007). Birc2 (cIap1) Regulates Endothelial Cell Integrity and Blood Vessel Homeostasis. Nat. Genet.39, 13971402. 10.1038/ng.2007.8

  • 103

    SayedN.HuangY.NguyenK.Krejciova-RajaniemiZ.GraweA. P.GaoT.et al (2021). An Inflammatory Aging Clock (iAge) Based on Deep Learning Tracks Multimorbidity, Immunosenescence, Frailty and Cardiovascular Aging. Nat. Aging1, 598615. 10.1038/s43587-021-00082-y

  • 104

    SchislerJ. C.GrevengoedT. J.PascualF.CooperD. E.EllisJ. M.PaulD. S.et al (2015). Cardiac Energy Dependence on Glucose Increases Metabolites Related to Glutathione and Activates Metabolic Genes Controlled by Mechanistic Target of Rapamycin. Jaha4 (2), e001136. 10.1161/JAHA.114.001136

  • 105

    ShadrinaA. S.ShashkovaT. I.TorgashevaA. A.SharapovS. Z.KlarićL.PakhomovE. D.et al (2020). Prioritization of Causal Genes for Coronary Artery Disease Based on Cumulative Evidence from Experimental and In Silico Studies. Sci. Rep.10, 10486. 10.1038/s41598-020-67001-w

  • 106

    ShiQ.ChenL. (2017). Cereblon: A Protein Crucial to the Multiple Functions of Immunomodulatory Drugs as Well as Cell Metabolism and Disease Generation. J. Immunol. Res.2017, 18. 10.1155/2017/9130608

  • 107

    ShimizuH.TaniguchiT.IshikawaY.YokoyamaM. (1997). Effects of Nitric Oxide on Cholesterol Metabolism in Macrophages. Atherosclerosis129 (2), 193198. 10.1016/s0021-9150(96)06032-7

  • 108

    ShirvaniS. M.MookanamparambilL.RamoniM. F.ChinM. T. (2007). Transcription Factor CHF1/Hey2 Regulates the Global Transcriptional Response to Platelet-Derived Growth Factor in Vascular Smooth Muscle Cells. Physiol. Genomics30, 61. 10.1152/physiolgenomics.00277.2006

  • 109

    SingerJ. W.FleischmanA.Al-FayoumiS.MascarenhasJ. O.YuQ.AgarwalA. (2018). Inhibition of Interleukin-1 Receptor-Associated Kinase 1 (IRAK1) as a Therapeutic Strategy. Oncotarget9 (70), 3341633439. 10.18632/oncotarget.26058

  • 110

    SoJ. S. (2018). Roles of Endoplasmic Reticulum Stress in Immune Responses. Mol. Cells41 (8), 705716. 10.14348/molcells.2018.0241

  • 111

    StadtmanE. R.Van RemmenH.RichardsonA.WehrN. B.LevineR. L. (2005). Methionine Oxidation and Aging. Biochimica Biophysica Acta (BBA) - Proteins Proteomics1703 (2), 135140. 10.1016/j.bbapap.2004.08.010

  • 112

    SteinK. C.Morales-PolancoF.van der LiendenJ.RainboltT. K.FrydmanJ. (2022). Ageing Exacerbates Ribosome Pausing to Disrupt Cotranslational Proteostasis. Nature601 (7894), 637642. 10.1038/s41586-021-04295-4

  • 113

    StewartJ.ManmathanG.WilkinsonP. (2017). Primary Prevention of Cardiovascular Disease: A Review of Contemporary Guidance and Literature. JRSM Cardiovasc. Dis.6, 204800401668721. 10.1177/2048004016687211

  • 114

    SwirskiF. K.NahrendorfM. (2014). Leukocyte Behavior in Atherosclerosis, Myocardial Infarction, and Heart Failure. Science339, 161. 10.1126/science.1230719

  • 115

    TakaiT. (2002). Roles of Fc Receptors in Autoimmunity. Nat. Rev. Immunol.2, 580592. 10.1038/nri856

  • 116

    TobaH.CannonP. L.YabluchanskiyA.IyerR. P.D’ArmientoJ.LindseyM. L. (2017). Transgenic Overexpression of Macrophage Matrix Metalloproteinase-9 Exacerbates Age-Related Cardiac Hypertrophy, Vessel Rarefaction, Inflammation, and Fibrosis. Am. J. Physiology-Heart Circulatory Physiology312, H375H383. 10.1152/ajpheart.00633.2016

  • 117

    TokarevA. A.AlfonsoA.SegevN. (2009). “Overview of Intracellular Compartments and Trafficking Pathways,” in Trafficking inside Cells. Molecular Biology Intelligence Unit (New York, NY: Springer). 10.1007/978-0-387-93877-6_1

  • 118

    TuriZ.LaceyM.MistrikM.MoudryP. (2019). Impaired Ribosome Biogenesis: Mechanisms and Relevance to Cancer and Aging. Aging11, 25122540. 10.18632/aging.101922

  • 119

    TyrrellD. J.GoldsteinD. R. (2020). Ageing and Atherosclerosis: Vascular Intrinsic and Extrinsic Factors and Potential Role of IL-6. Nat. Rev. Cardiol.18, 5868. 10.1038/s41569-020-0431-7

  • 120

    Van-AsscheT.HuygelenV.CrabtreeM. J.AntoniadesC. (2011). Gene Therapy Targeting Inflammation in Atherosclerosis. Curr. Pharm. Des.17 (37), 42104223. 10.2174/138161211798764799

  • 121

    WågsäterD.PaloschiV.HanemaaijerR.HultenbyK.BankR. A.Franco‐CerecedaA.et al (2013). Impaired Collagen Biosynthesis and Cross‐linking in Aorta of Patients with Bicuspid Aortic Valve. Jaha2 (1), e000034. 10.1161/JAHA.112.000034

  • 122

    WangJ. C.BennettM. (2012). Aging and Atherosclerosis. Circ. Res.111 (2), 245259. 10.1161/CIRCRESAHA.111.261388

  • 123

    WangL.HouJ.WangJ.ZhuZ.ZhangW.ZhangX.et al (2020). Regulatory Roles of HSPA6 in Actinidia Chinensis Planch. Root Extract (acRoots)‐inhibited Lung Cancer Proliferation. Clin. Transl. Med.10, 2. 10.1002/ctm2.46

  • 124

    WangX.ChenT.DengZ.GaoW.LiangT.QiuX.et al (2021). Melatonin Promotes Bone Marrow Mesenchymal Stem Cell Osteogenic Differentiation and Prevents Osteoporosis Development through Modulating Circ_0003865 that Sponges miR-3653-3p. Stem Cell Res. Ther.12, 150. 10.1186/s13287-021-02224-w

  • 125

    WangY.HuangT.LiY.ShaX. (2020). The Self-Organization Model Reveals Systematic Characteristics of Aging. Theor. Biol. Med. Model.17, 4. 10.1186/s12976-020-00120-z

  • 126

    WeiT.LiuJ.ZhangD.WangX.LiG.MaR.et al (2021). The Relationship between Nutrition and Atherosclerosis. Front. Bioeng. Biotechnol.9, 635504. 10.3389/fbioe.2021.635504

  • 127

    WuW.YuN.LiF.GaoP.LinS.ZhuY. (2021). RPL35 Promotes Neuroblastoma Progression via the Enhanced Aerobic Glycolysis. Am. J. Cancer Res.11 (11), 57015714.

  • 128

    XiaS.ZhangX.ZhengS.KhanabdaliR.KalionisB.WuJ.et al (2016). An Update on Inflamm-Aging: Mechanisms, Prevention, and Treatment. J. Immunol. Res.2016, 112. 10.1155/2016/8426874

  • 129

    XuM.KnoxA. J.MichaelisK. A.Kiseljak-VassiliadesK.Kleinschmidt-DeB. K.LilleheiK. O.et al (2012). Reprimo (RPRM) Is a Novel Tumor Suppressor in Pituitary Tumors and Regulates Survival, Proliferation, and Tumorigenicity. Endocrinology153 (7), 29632973. 10.1210/en.2011-2021

  • 130

    YangJ.BoermM.McCartyM.BucanaC.FidlerI. J.Z J YangY.et al (2000). Mekk3 Is Essential for Early Embryonic Cardiovascular Development. Nat. Genet.24 (3), 309313. 10.1038/73550

  • 131

    YazdaniA.YazdaniA.SamieiA.BoerwinkleE. (2016). Generating a Robust Statistical Causal Structure over 13 Cardiovascular Disease Risk Factors Using Genomics Data. J. Biomed. Inf.60 (C), 114119. 10.1016/j.jbi.2016.01.012

  • 132

    YeG.-C.LiuY.-F.HuangL.ZhangC.-Y.ShengY.-L.WuB.et al (2021). Key microRNAs and Hub Genes Associated with Poor Prognosis in Lung Adenocarcinoma. Aging13 (3), 37423762. 10.18632/aging.202337

  • 133

    YinJ.IbrahimS.PetersenF.YuX. (2021). Autoimmunomic Signatures of Aging and Age-Related Neurodegenerative Diseases Are Associated with Brain Function and Ribosomal Proteins. Front. Aging Neurosci.13, 679688. 10.3389/fnagi.2021.679688

  • 134

    YuJ.MursuE.TyppöM.LaaksonenS.VoipioH.-M.PesonenP.et al (2019). MMP-3 and MMP-8 in Rat Mandibular Condylar Cartilage Associated with Dietary Loading, Estrogen Level, and Aging. Archives Oral Biol.97, 238244. 10.1016/j.archoralbio.2018.10.037

  • 135

    YuanB.XuY.ZhengS. (2021). PLOD1 Acts as a Tumor Promoter in Glioma via Activation of the HSF1 Signaling Pathway. Mol. Cell. Biochem.477, 549. 10.1007/s11010-021-04289-w

  • 136

    ZhangR. (2013). MNADK, a Novel Liver-Enriched Mitochondrion-Localized NAD Kinase. Biol. Open2 (4), 432438. 10.1242/bio.20134259

  • 137

    ZhangY.WuY.SuX. (2021). PLOD1 Promotes Cell Growth and Aerobic Glycolysis by Regulating the SOX9/PI3K/Akt/mTOR Signaling Pathway in Gastric Cancer. Front. Bioscience-Landmark26, 322. 10.52586/4946

  • 138

    ZhavoronkovA.MamoshinaP.VanhaelenQ.Scheibye-KnudsenM.MoskalevA.AliperA. (2019). Artificial Intelligence for Aging and Longevity Research: Recent Advances and Perspectives. Ageing Res. Rev.49, 4966. 10.1016/j.arr.2018.11.003

  • 139

    ZhengL.XuH.DiY.ChenL.LiuJ.KangL.et al (2021). ELK4 Promotes the Development of Gastric Cancer by Inducing M2 Polarization of Macrophages through Regulation of the KDM5A-PJA2-KSR1 axis. J. Transl. Med.19, 342. 10.1186/s12967-021-02915-1

  • 140

    ZhengY.WangC.SongA.JiangF.ZhouJ.LiG.et al (2020). CMTM6 Promotes Cell Proliferation and Invasion in Oral Squamous Cell Carcinoma by Interacting with NRP1. Am. J. Cancer Res.10 (6), 16911709.

  • 141

    ZhengY.LiuH.KongY. (2017). miR-188 Promotes Senescence of Lineage-Negative Bone Marrow Cells by Targeting MAP3K3 Expression. FEBS Lett.591 (15), 22902298. 10.1002/1873-3468.12720

Summary

Keywords

inflamm-aging, atherosclerosis, immune homeostasis, causal analysis, sensitive analysis

Citation

He Y, Chen Y, Yao L, Wang J, Sha X and Wang Y (2022) The Inflamm-Aging Model Identifies Key Risk Factors in Atherosclerosis. Front. Genet. 13:865827. doi: 10.3389/fgene.2022.865827

Received

30 January 2022

Accepted

26 April 2022

Published

30 May 2022

Volume

13 - 2022

Edited by

Xiao Dong, University of Minnesota Twin Cities, United States

Reviewed by

Reza Jabal, Albert Einstein College of Medicine, United States

Zhi Liu, Nanjing Medical University, China

Updates

Copyright

*Correspondence: Yin Wang,

This article was submitted to Genetics of Aging, a section of the journal Frontiers in Genetics

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics