The Identification of Candidate Biomarkers and Pathways in Atherosclerosis by Integrated Bioinformatics Analysis

Unveiling Atherosclerotic Plaque Heterogeneity and SPP1+/VCAN+ Macrophage Subtype Prognostic Significance Through Integrative Single-Cell and Bulk-Seq Analysis

Background

Dysregulated macrophages are important causes of Atherosclerosis (AS) formation and increased plaque instability, but the heterogeneity of these plaques and the role of macrophage subtypes in plaque instability have yet to be clarified.

Methods

This study integrates single-cell and bulk-seq data to analyze atherosclerotic plaques. Unsupervised clustering was used to reveal distinct plaque subtypes, while survival analysis and gene set variation analysis (GSVA) methods helped in understanding their clinical outcomes. Enrichment of differential expression of macrophage genes (DEMGs) score and pseudo-trajectory analysis were utilized to explore the biological functions and differentiation stages of macrophage subtypes in AS progression. Additionally, CellChat and the BayesPrism deconvolution method were used to elucidate macrophage subtype interaction and their prognostic significance at single-cell resolution. Finally, the expression of biomarkers was validated in mouse experiments.

Results

Three distinct AS plaque subtypes were identified, with cluster 3 plaque subtype being particularly associated with higher immune infiltration and poorer prognosis. The DEMGs score exhibited a significant elevation in three macrophage subtypes (SPP1+/VCAN+ macrophages, IL1B+ macrophages, and FLT3LG+ macrophages), associated with cluster 3 plaque subtype and highlighted the prognostic significance of these subtypes. Activation trajectory of the macrophage subtypes is divided into three states (Pre-branch, Cell fate 1, and Cell fate 2), and Cell fate 2 (SPP1+/VCAN+ macrophages, IL1B+ macrophages, and FLT3LG+ macrophages dominant) exhibiting the highest DEMGs score, distinct interactions with other cell components, and relating to poorer prognosis of ischemic events. This study also uncovered a unique SPP1+/VCAN+ macrophage subtype, rare in quantity but significant in influencing AS progression. Machine learning algorithms identified 10 biomarkers crucial for AS diagnosis. The validation of these biomarkers was performed using Mendelian Randomization analysis and in vitro methods, supporting their relevance in AS pathology.

Conclusion

Our study provides a comprehensive view of AS plaque heterogeneity and the prognostic significance of macrophage subtypes in plaque instability.

Associations of genetic markers of diabetes mellitus with carotid atherosclerosis: a community-based case-control study

BACKGROUND

Diabetes mellitus (DM) is a well-established determinant of atherosclerosis and cardiovascular diseases (CVD). Recently, genome-wide association studies (GWAS) identified several single nucleotide polymorphism (SNP) significantly correlated with DM. The study aimed to explore the relationships of the top significant DM SNPs with carotid atherosclerosis (CA).

METHODS

We used a case-control design and randomly selected 309 cases and 439 controls with and without, respectively, carotid plaque (CP) from a community-based cohort. Eight recent GWAS on DM in East Asians reported hundreds of SNPs with genome-wide significance. The study used the top significant DM SNPs, with a p-value < 10^-16, as the candidate genetic markers of CA. The independent effects of these DM SNPs on CA were assessed by multivariable logistic regression analyses to control the effects of conventional cardio-metabolic risk factors.

RESULTS

Multivariable analyses showed that, 9 SNPs, including rs4712524, rs1150777, rs10842993, rs2858980, rs9583907, rs1077476, rs7180016, rs4383154, and rs9937354, showed promising associations with the presence of carotid plaque (CP). Among them, rs9937354, rs10842993, rs7180016, and rs4383154 showed significantly independent effects. The means (SD) of the 9-locus genetic risk score (9-GRS) of CP-positive and -negative subjects were 9.19 (1.53) and 8.62 (1.63), respectively (p < 0.001). The corresponding values of 4-locus GRS (4-GRS) were 4.02 (0.81) and. 3.78 (0.92), respectively (p < 0.001). The multivariable-adjusted odds ratio of having CP for per 1.0 increase in 9-GRS and 4-GRS were 1.30 (95% CI 1.18-1.44; p = 4.7 × 10^-7) and 1.47 (95% CI 1.74-9.40; p = 6.1 × 10^-5), respectively. The means of multi-locus GRSs of DM patients were similar to those of CP-positive subjects and higher than those of CP-negative or DM-negative subjects.

CONCLUSIONS

We identified 9 DM SNPs showing promising associations with CP. The multi-locus GRSs may be used as biomarkers for the identification and prediction of high-risks subjects for atherosclerosis and atherosclerotic diseases. Future studies on these specific SNPs and their associated genes may provide valuable information for the preventions of DM and atherosclerosis.

Identification of immune-related genes in diagnosing atherosclerosis with rheumatoid arthritis through bioinformatics analysis and machine learning

Background

Increasing evidence has proven that rheumatoid arthritis (RA) can aggravate atherosclerosis (AS), and we aimed to explore potential diagnostic genes for patients with AS and RA.

Methods

We obtained the data from public databases, including Gene Expression Omnibus (GEO) and STRING, and obtained the differentially expressed genes (DEGs) and module genes with Limma and weighted gene co-expression network analysis (WGCNA). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis, the protein-protein interaction (PPI) network, and machine learning algorithms [least absolute shrinkage and selection operator (LASSO) regression and random forest] were performed to explore the immune-related hub genes. We used a nomogram and receiver operating characteristic (ROC) curve to assess the diagnostic efficacy, which has been validated with GSE55235 and GSE73754. Finally, immune infiltration was developed in AS.

Results

The AS dataset included 5,322 DEGs, while there were 1,439 DEGs and 206 module genes in RA. The intersection of DEGs for AS and crucial genes for RA was 53, which were involved in immunity. After the PPI network and machine learning construction, six hub genes were used for the construction of a nomogram and for diagnostic efficacy assessment, which showed great diagnostic value (area under the curve from 0.723 to 1). Immune infiltration also revealed the disorder of immunocytes.

Conclusion

Six immune-related hub genes (NFIL3, EED, GRK2, MAP3K11, RMI1, and TPST1) were recognized, and the nomogram was developed for AS with RA diagnosis.

Genome-wide associated variants of subclinical atherosclerosis among young people with HIV and gene-environment interactions

BACKGROUND

Genome-wide association studies (GWAS) have identified some variants associated with subclinical atherosclerosis (SCA) in general population but lacking sufficient validation. Besides traditional risk factors, whether and how would genetic variants associate with SCA among people with HIV (PWH) remains to be elucidated.

METHOD

A large original GWAS and gene-environment interaction analysis of SCA were conducted among Chinese PWH (n = 2850) and age/sex-matched HIV-negative controls (n = 5410). Subgroup analyses by age and functional annotations of variants were also performed.

RESULTS

Different from HIV-negative counterparts, host genome had a greater impact on young PWH rather than the elders: one genome-wide significant variant (rs77741796, P = 2.20 × 10^-9) and eight suggestively significant variants (P < 1 × 10^-6) were identified to be specifically associated with SCA among PWH younger than 45 years. Seven genomic loci and 15 genes were mapped to play a potential role on SCA among young PWH, which were enriched in the biological processes of atrial cardiac muscle cell membrane repolarization and molecular function of protein kinase A subunit binding. Furthermore, genome-wide interaction analyses revealed significant HIV-gene interactions overall as well as gene-environment interactions with alcohol consumption, tobacco use and obesity among PWH. The identified gene-environment interaction on SCA among PWH might be useful for discovering high-risk individuals for the prevention of SCA, particularly among those with tobacco use and alcohol consumption.

CONCLUSION

The present study provides new clues for the genetic contribution of SCA among young PWH and is the starting point of precision intervention targeting HIV-related atherosclerosis.

Identifying potential signatures for atherosclerosis in the context of predictive, preventive, and personalized medicine using integrative bioinformatics approaches and machine-learning strategies

Background

Atherosclerosis is a major contributor to morbidity and mortality worldwide. Although several molecular markers associated with atherosclerosis have been developed in recent years, the lack of robust evidence hinders their clinical applications. For these reasons, identification of novel and robust biomarkers will directly contribute to atherosclerosis management in the context of predictive, preventive, and personalized medicine (PPPM). This integrative analysis aimed to identify critical genetic markers of atherosclerosis and further explore the underlying molecular immune mechanism attributing to the altered biomarkers.

Methods

Gene Expression Omnibus (GEO) series datasets were downloaded from GEO. Firstly, differential expression analysis and functional analysis were conducted. Multiple machine-learning strategies were then employed to screen and determine key genetic markers, and receiver operating characteristic (ROC) analysis was used to assess diagnostic value. Subsequently, cell-type identification by estimating relative subsets of RNA transcript (CIBERSORT) and a single-cell RNA sequencing (scRNA-seq) data were performed to explore relationships between signatures and immune cells. Lastly, we validated the biomarkers' expression in human and mice experiments.

Results

A total of 611 overlapping differentially expressed genes (DEGs) included 361 upregulated and 250 downregulated genes. Based on the enrichment analysis, DEGs were mapped in terms related to immune cell involvements, immune activating process, and inflaming signals. After using multiple machine-learning strategies, dehydrogenase/reductase 9 (DHRS9) and protein tyrosine phosphatase receptor type J (PTPRJ) were identified as critical biomarkers and presented their high diagnostic accuracy for atherosclerosis. From CIBERSORT analysis, both DHRS9 and PTPRJ were significantly related to diverse immune cells, such as macrophages and mast cells. Further scRNA-seq analysis indicated DHRS9 was specifically upregulated in macrophages of atherosclerotic lesions, which was confirmed in atherosclerotic patients and mice.

Conclusions

Our findings are the first to report the involvement of DHRS9 in the atherogenesis, and the proatherogenic effect of DHRS9 is mediated by immune mechanism. In addition, we confirm that DHRS9 is localized in macrophages within atherosclerotic plaques. Therefore, upregulated DHRS9 could be a novel potential target for the future predictive diagnostics, targeted prevention, patient stratification, and personalization of medical services in atherosclerosis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-022-00289-y.

Proteomics-based evaluation of the mechanism underlying vascular injury via DNA interstrand crosslinks, glutathione perturbation, mitogen-activated protein kinase, and Wnt and ErbB signaling pathways induced by crotonaldehyde

Crotonaldehyde (CRA)—one of the major environmental pollutants from tobacco smoke and industrial pollution—is associated with vascular injury (VI). We used proteomics to systematically characterize the presently unclear molecular mechanism of VI and to identify new related targets or signaling pathways after exposure to CRA. Cell survival assays were used to assess DNA damage, whereas oxidative stress was determined using colorimetric assays and by quantitative fluorescence study; additionally, cyclooxygenase-2, mitogen-activated protein kinase pathways, Wnt3a, β-catenin, phospho-ErbB2, and phospho-ErbB4 were assessed using ELISA. Proteins were quantitated via tandem mass tag-based liquid chromatography-mass spectrometry and bioinformatics analyses, and 34 differentially expressed proteins were confirmed using parallel reaction monitoring, which were defined as new indicators related to the mechanism underlying DNA damage; glutathione perturbation; mitogen-activated protein kinase; and the Wnt and ErbB signaling pathways in VI based on Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and protein–protein interaction network analyses. Parallel reaction monitoring confirmed significant (p < 0.05) upregulation (> 1.5-fold change) of 23 proteins and downregulation (< 0.667-fold change) of 11. The mechanisms of DNA interstrand crosslinks; glutathione perturbation; mitogen-activated protein kinase; cyclooxygenase-2; and the Wnt and ErbB signaling pathways may contribute to VI through their roles in DNA damage, oxidative stress, inflammation, vascular dysfunction, endothelial dysfunction, vascular remodeling, coagulation cascade, and the newly determined signaling pathways. Moreover, the Wnt and ErbB signaling pathways were identified as new disease pathways involved in VI. Taken together, the elucidated underlying mechanisms may help broaden existing understanding of the molecular mechanisms of VI induced by CRA.

Turning Hot into Cold: Immune Microenvironment Reshaping for Atherosclerosis Attenuation Based on pH-Responsive shSiglec-1 Delivery System

Current atherosclerosis treatment is based on a combination of cholesterol-lowering medication and low-fat diets; however, the clinical effect is unsatisfactory. It has been shown that the level of immune cell infiltration and pro-inflammatory factors in the atherosclerotic immune microenvironment (AIM) play important roles in the development and progression of atherosclerosis. Therefore, we hypothesized that reshaping "hot AIM" into "cold AIM" could attenuate atherosclerosis. For this purpose, we designed a pH-responsive and charge-reversible nanosystem, referred to as Au-PEI/shSiglec-1/PEI-acetylsalicylic acid (ASPA NPs) to effectively deliver shSiglec-1, which blocked the interactions between macrophages with CD8⁺ T/NKT cells, thus inhibiting immune cell infiltration. Further, we demonstrated that acetylsalicylic acid (ASA), detached from the pH-responsive PEI-ASA polymer, and inhibited lipid accumulation in macrophage, thereby decreasing the lipid antigen presentation. Additionally, reduced macrophage-produced inflammatory factors by ASA and low CD8⁺ T/NKT cell infiltration levels synergistically inhibit Th17 cell differentiation, thus further dramatically attenuating inflammation in AIM by decreasing the IL-17A production. Eventually, ASPA NPs efficiently reshaped AIM by inhibiting immune cell infiltration, lipid antigen presentation, and pro-inflammation, which provided a feasible therapeutic strategy for atherosclerosis immunotherapy.

Multifaceted Roles of the E3 Ubiquitin Ligase RING Finger Protein 115 in Immunity and Diseases

Ubiquitination is a post-translational modification that plays essential roles in various physiological and pathological processes. Protein ubiquitination depends on E3 ubiquitin ligases that catalyze the conjugation of ubiquitin molecules on lysine residues of targeted substrates. RING finger protein 115 (RNF115), also known as breast cancer associated gene 2 (BCA2) and Rab7-interacting RING finger protein (Rabring7), has been identified as a highly expressed protein in breast cancer cells and tissues. Later, it has been demonstrated that RNF115 catalyzes ubiquitination of a series of proteins to modulate a number of signaling pathways, and thereby regulates viral infections, autoimmunity, cell proliferation and death and tumorigenesis. In this review, we introduce the identification, expression and activity regulation of RNF115, summarize the substrates and functions of RNF115 in different pathways, and discuss the roles of RNF115 as a biomarker or therapeutic target in diseases.

Integrated Bioinformatic Analysis of the Shared Molecular Mechanisms Between Osteoporosis and Atherosclerosis

BACKGROUND

Osteoporosis and atherosclerosis are common in the elderly population, conferring a heavy worldwide burden. Evidence links osteoporosis and atherosclerosis but the exact underlying common mechanism of its occurrence is unclear. The purpose of this study is to further explore the molecular mechanism between osteoporosis and atherosclerosis through integrated bioinformatic analysis.

METHODS

The microarray data of osteoporosis and atherosclerosis in the Gene Expression Omnibus (GEO) database were downloaded. The Weighted Gene Co-Expression Network Analysis (WGCNA) and differentially expressed genes (DEGs) analysis were used to identify the co-expression genes related to osteoporosis and atherosclerosis. In addition, the common gene targets of osteoporosis and atherosclerosis were analyzed and screened through three public databases (CTD, DISEASES, and GeneCards). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by Metascape. Then, the common microRNAs (miRNAs) in osteoporosis and atherosclerosis were screened out from the Human microRNA Disease Database (HMDD) and the target genes of whom were predicted through the miRTarbase. Finally, the common miRNAs-genes network was constructed by Cytoscape software.

RESULTS

The results of common genes analysis showed that immune and inflammatory response may be a common feature in the pathophysiology of osteoporosis and atherosclerosis. Six hub genes (namely, COL1A1, IBSP, CTSD, RAC2, MAF, and THBS1) were obtained via taking interaction of different analysis results. The miRNAs-genes network showed that has-let-7g might play an important role in the common mechanisms between osteoporosis and atherosclerosis.

CONCLUSION

This study provides new sights into shared molecular mechanisms between osteoporosis and atherosclerosis. These common pathways and hub genes may offer promising clues for further experimental studies.