Characterization of Distinct Monocyte Subtypes and Immune Features Associated with HIV, Tuberculosis, and Coronary Artery Disease in a Ugandan Cohort Using Mass Cytometry

Authors

José Cobeña-Reyes¹, Celestine N. Wanjalla^2,3,4, Manuel G. Feria⁵, Joshua Simmons ², Tecla Temu ⁶, Cindy Nochowicz ², Sheikh Yasir Arafat^1,7, Cissy Kityo ⁸, Geofrey Erem ⁹, Christopher T. Longenecker¹⁰, Sandra Andorf^1,11,12,*, Moises A. Huaman^5,*

Affiliated Institutions

¹ Division of Biomedical Informatics, Cincinnati; Children’s Hospital Medical Center,
Cincinnati, Ohio
² Division of Infectious Diseases, Vanderbilt University Medical Center, Nashville, Tennessee
³ The Center for AIDS Health Disparities Research, Nashville, Tennessee
⁴Veteran’s Health Administration, Tennessee Valley Healthcare System, Nashville, Tennessee
⁵Division of Infectious Diseases, University of Cincinnati College of Medicine, Cincinnati, Ohio
⁶Department of Pathology, Mass General Brigham, Harvard Medical School, Boston, Massachusetts
⁷ Department of Biostatistics, Health Informatics, and Data Sciences, College of Medicine, University of Cincinnati, Cincinnati, Ohio
⁸Joint Clinical Research Center, Kampala, Uganda
⁹Makerere University, Kampala, Uganda
¹⁰Division of Cardiology, Department of Global Health, University of Washington, Seattle, Washington
¹¹Division of Allergy and Immunology, Cincinnati Children’s Hospital Medical Center, Cincinnati, Ohio
¹²Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio
^•Contributed equally as senior authors

Corresponding Author

Moises A. Huaman
huamanms@ucmail.uc.edu

Suggested Citation

Cobeña-Reyes J, Wanjalla CN, Feria MG, Simmons J , Temu T , Nochowicz C , Arafat SA, Kityo C , Erem G , Longenecker CT, Andorf S, Huaman MA. Characterization of Distinct Monocyte Subtypes and Immune Features Associated with HIV, Tuberculosis, and Coronary Artery Disease in a Ugandan Cohort Using Mass Cytometry. Pathogens and Immunity. 2025;11(1):14-38. doi: 10.20411/pai.v11i1.945

DOI

10.20411/pai.v11i1.945

Abstract

Background: Coronary artery disease (CAD), tuberculosis (TB), and HIV are major global health concerns. Individuals affected by one or more of these conditions often exhibit chronic inflammation and immune dysregulation, with monocytes playing a central role. Monocyte subsets are known to expand in individuals with HIV, TB, or CAD, but the mechanisms by which these cells contribute to inflammation and immune responses remain poorly understood.

Methods: We employed high-dimensional mass cytometry to characterize monocyte heterogeneity in 61 Ugandan adults with varying combinations of HIV, TB, and subclinical or overt CAD. An integrative approach was used, combining manual gating, unsupervised clustering, and machine learning to identify distinct monocyte phenotypes associated with CAD and TB. Monocyte activation markers soluble CD14 (sCD14) and sCD163 were measured in plasma. CAD was diagnosed by coronary computed tomography angiography. TB was determined by a questionnaire and interferon-gamma release assay (IGRA) testing.

Results: Participants’ demographics and clinical characteristics were similar by CAD or HIV/TB status. Median age was 61 years; 37.7% were female. People living with HIV and latent TB or prior active TB had higher sCD14 plasma levels compared with HIV/TB-negative individuals. Individuals with CAD showed reduced surface expression of the scavenger receptor CD163 on non-classical monocytes. Unsupervised clustering further revealed 2 distinct non-classical monocyte subsets associated with disease states: A CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ subset significantly depleted in individuals with CAD, and a CD86+ CX3CR1++ CD45RA++ GPR56- CD38- CXCR3- subset enriched in individuals with latent TB.

Conclusions: These findings underscore the complexity of the monocyte landscape in CAD progression, particularly in regions where HIV and TB are co-endemic. Our study reveals distinct alterations within 2 non-classical monocyte subpopulations associated with CAD and with HIV/TB, offering mechanistic insights that may support the development of precision biomarkers and immune-targeted therapies across these disease contexts.

Keywords

Coronary Artery Disease; Tuberculosis; HIV; Monocytes; Mass Cytometry; Immune Dysregulation; CD163; CX3CR1; Biomarker; Unsupervised Clustering; Elastic Net

Introduction

HIV, [1] tuberculosis (TB), [2] and cardiovascular diseases [3, 4] are major public health concerns and leading causes of morbidity and mortality worldwide. Individuals residing in Sub-Saharan Africa are particularly vulnerable [5–7], often facing the added challenge of living with one or more of these comorbid conditions simultaneously. People living with HIV (PWH) experience complex immune dysregulation despite appropriate antiretroviral therapy, which may contribute to various comorbidities such as coronary artery disease (CAD) [8] or increased susceptibility to other infectious diseases such as TB [9]. Similarly, people with TB exhibit altered immune responses that can exacerbate immune dysfunction and further elevate CAD risk [10, 11].

Monocytes are cells of the innate immune system that play a critical role in the response to infection and in regulating the immune response. Monocytes are classified based on their CD14 and CD16 expression as classical monocytes (CD14++ CD16-), intermediate monocytes (CD14++ CD16+) and non-classical monocytes (CD14+ CD16++). In PWH, TB, or CAD, monocyte populations undergo significant alterations in both frequency and function [12–15]. In CAD, lower levels of classical monocytes circulating in the blood are usually encountered [16]. Intermediate monocytes tend to aggregate and form monocyte-platelet aggregates in people with ST-elevation myocardial infarction due to coronary atherothrombosis [17]. Further, monocytes infiltrate plaque and differentiate into macrophages (ie, monocyte-derived macrophages). However, these macrophages exhibit a decreased ability to migrate, promoting atherogenesis [18, 19]. In TB infection, CD16+ monocyte subsets are expanded, reflecting alterations in the composition of circulating monocyte populations [20]. The expansion is, however, reversed by anti-TB treatment [20]. Classical monocytes that express inflammatory markers typically increase as the disease progresses, whereas the percentage of non-classical monocytes expressing anti-pathogen infection markers decreases [21]. In PWH, alterations in the distribution of monocyte subsets are also observed, including an increase in intermediate and non-classical monocytes [22]. Interestingly, intermediate monocyte subsets are preferentially infected with HIV due to the facilitation of viral entry [23, 24]. Moreover, there are indications that CD16+ subsets might serve as a reservoir for HIV, complicating viral suppression [25]. Understanding changes in the distribution and phenotype of monocyte subsets is crucial for uncovering disease mechanisms and potential therapeutic targets.

In this work, to characterize monocyte populations in individuals with HIV, TB, and CAD, a cohort of people living with these conditions was studied. We hypothesized that within the spectrum of HIV, TB, and CAD, there would be identifiable monocyte subpopulations that would exhibit distinct phenotypic or functional marker expression patterns. Using mass cytometry data, both manually gated monocyte subtypes and those identified through unsupervised clustering were examined, and their distributions were compared across groups defined by TB, HIV, and CAD status. Additionally, elastic net regression was employed to identify monocyte features most strongly associated with CAD status. This statistical approach enabled the identification of previously undefined immune features that may serve as correlates for disease progression and potential therapeutic intervention.

Given the critical role of monocyte subsets in the regulation of HIV, TB, and CAD, this study aimed to further uncover critical insights into how monocyte populations are altered in HIV, TB and CAD by leveraging high-parameter mass cytometry data and advanced statistical analysis.

Methods

Study Design and Participants

Mass cytometry data from manually gated monocytes isolated from peripheral blood mononuclear cells (PBMCs) of 61 participants were analyzed. All the participants were enrolled at the Joint Clinical Research Center in Kampala, Uganda, in the Ugandan study of HIV effects on the myocardium and atherosclerosis (mUTIMA). For this analysis, the deidentified dataset contained information on individuals living with or without HIV, TB, and CAD. The protocol of the parent study was previously described [26]. The study was approved by the University Hospitals Cleveland Medical Center Institutional Review Board, the Joint Clinical Research Centre Research Ethics Committee, and the Uganda National Council for Science and Technology. All participants signed written informed consent.

TB infection status was defined as latent tuberculosis infection (TBI+) or negative (TB-) based on comprehensive TB symptom questionnaires and interferon-gamma release assay (IGRA) testing (QuantiFERON®-TB). Previous active tuberculosis (TBpr) was defined by self-reported history of prior active TB and completion of TB treatment, which was confirmed by review of medical records. HIV status of HIV-negative (HIV-) participants was confirmed with a rapid HIV test (HIV 1/2 STAT-PAK®; FDA-approved). The Atherosclerotic Cardiovascular Disease (ASCVD) pooled cohort equations (PCE) were used to calculate 10-year ASCVD risk scores [27]. All participants underwent coronary computed tomography angiography to quantify the presence and severity of CAD. In addition, 2 soluble markers of monocyte activation were measured in plasma from all participants [28]; sCD14 and sCD163, using ELISA (R&D Systems).

Two categories were employed based on the CAD status: A binary classification defined as the presence or absence of coronary plaque (CAD-/+) and a categorization based on burden of disease using the segment involvement score (SIS), an integer-based measure that assigns a value of 1 to each coronary artery segment with detectable atherosclerotic plaque, irrespective of plaque size or degree of luminal obstruction [29]: CAD-, CAD+ minimum (CAD min): ≤ 2 SIS; and CAD+ greater than minimum (CAD+ >min): > 2 SIS. Five study groups were also defined by combining the TB and HIV status as HIV-TB-, HIV-TBI+, HIV+TB-, HIV+TBI+, and HIV+TBpr.

Mass Cytometry

All the CyTOF experiments were conducted at Vanderbilt University Medical Center through the Center for AIDS Research (CFAR) core. Reagents were purchased from Fluidigm, now Standard BioTools (see Supplementary Table 1 for more details). All samples were run in a Helios CyTOF 3.0 instrument.

A solution of 90% FBS (Fetal Bovine Serum) with 10% DMSO (dimethyl sulfoxide) was used for cryopreservation of PBMCs, as it preserves the expression of markers included in this panel. This panel has been used to stain both cryopreserved and freshly isolated PBMCs, without significant changes in marker expression.

Cryopreserved PBMCs were rapidly thawed in a 37 °C water bath and immediately transferred to 10 mL ice‐cold RPMI-1640 (no serum) containing 20 μL Nuclease S7 (Roche) to degrade extracellular nucleic acids. Cells were washed twice by centrifugation at 420 × g for 10 minutes at room temperature (RT) in 10 mL phosphate‐buffered saline (PBS; without Ca²⁺/Mg²⁺), decanted, and gently resuspended in the residual supernatant. Viable cells were counted by trypan‐blue exclusion on a hemocytometer, and aliquots of 3 × 10⁶ cells were transferred to 5 mL polystyrene tubes (Falcon) for staining.

Cells were washed once in 2 mL PBS/1% bovine serum albumin (BSA; Sigma-Aldrich) at 420 × g × 10 minutes (RT) and then stained with 2 μL Live/Dead Rhodium-103 (Fluidigm; 1:500 final) for 15 minutes at 37 °C. After 2 additional washes in PBS/1% BSA (as above), pellets were tapped dry, then resuspended in 50 μL staining buffer (PBS/1% BSA). A master mix of metal‐conjugated surface antibodies (see Supplementary Table 1) and EQ™ Four Element Calibration Beads (Fluidigm) was prepared to a volume of 40 μL per sample, brought to ~100 μL total staining volume with staining buffer, and added to each tube. Samples were incubated for 30 minutes at RT, washed twice in 2 mL PBS (420 × g × 10 minutes, RT), and fixed by adding 20 μL 16% paraformaldehyde (PFA; Thermo Fisher) to 200 μL cells (final 1.6% PFA), incubating 15 minutes at RT.

Fixed cells were washed once in 2 mL PBS (800 × g × 7 minutes, RT), aspirated, and permeabilized by adding 1 mL of ice-cold methanol (MeOH), then stored at –20 °C overnight. Cells were then labeled with 2 μL of 25 μM iridium DNA intercalator (Fluidigm) in the presence of 1.6% PFA (20 μL PFA added to 200 μL cells) for 20 minutes at RT, followed by an overnight incubation at 4 °C.

Immediately before acquisition, cells were washed once in PBS and once in Millipore‐grade H₂O (800 × g × 7 minutes), resuspended at 5 × 10⁵ cells/mL in double‐distilled H₂O, and spiked with equilibration (EQ) beads at a 1:10 bead‐to‐cell ratio (vortex beads vigorously before use). Finally, cell suspensions were passed through a 35 μm nylon filter (CellTrics, Sysmex) and acquired on the Helios at <500 events.

Monocyte Phenotyping: Manual Gating

Conventional monocyte subsets were manually gated from PBMCs into total monocyte (TM), classical monocyte (CM), intermediate monocyte (IM), and non-classical monocyte (nCM) categories based on CD14 and CD16 expression, as previously described [6] (Supplementary Figure 1).

For each of the monocyte subsets, several populations were identified by manual gating. For the TM, frequencies of the following populations were determined: CD14+, CD14+ HLA-DR+, CD14+ CX3CR1+, CD14+ CD86+, CD14+ CD163+. For the CM, the following subpopulations were considered: CD14+ CD16-, CD14+ CD16- CX3CR1+, CD14+ CD16- CD86+, CD14+ CD16- CD163+. Subpopulations based on the same markers were studied for IM (CD14+ CD16-, CD14+ CD16- CX3CR1, CD14+ CD16- CD86+, CD14+ CD16- CD163+) and nCM (CD14- CD16+, CD14- CD16+ CX3CR1+, CD14- CD16+ CD86+, CD14- CD16+ CD163+).

Additionally, the Median Signal Intensity (MSI) of markers HLA-DR, CX3CR1, and CD86 was computed for TM, CM, IM, and nCM. MSI of CD163 was computed for TM, IM, and nCM.

Monocyte Phenotyping: Unsupervised Clustering

FCS files containing the total monocytes (TM) obtained from manual gating were read using the flowCore [30] package from the Bioconductor [31] project, and transformed using the arcsinh function with a cofactor of 5 [32]. Ten thousand cells were sampled randomly from each file. If a file contained fewer than 10,000 cells, all cells were used. Of the files, 46 contained more than 10,000 cells, while 15 contained fewer. Unsupervised clustering via FlowSOM [33] with default parameters and a predetermined number of 12 clusters was used to identify phenotypically distinct monocyte populations. All markers were used for unsupervised clustering. For each of the samples, the percentage of cells assigned to each of these clusters (ie, monocyte populations) was calculated.

To determine the marker expression within each monocyte population, scaled median expression values [34] were visualized using a heatmap. To visualize 2D map projections of the cell populations, the stratified CAD/SIS study groups (CAD-, CAD+ min and CAD+ > min), were used. For each group, 80,000 cells were randomly sampled and plotted using Uniform Manifold Approximation Projections (UMAPs) with colors for the populations based on the unsupervised clustering overlayed [35].

Elastic Net Regularization

The features from manual gating and the population percentages from FlowSOM were used to fit a generalized linear model, incorporating a regularization step through elastic net [36], to identify the features most strongly associated with CAD status. In this model, due to the small sample size, the CAD-/+ classification was used as the predicted variable. In the context of CAD status, negative coefficients from the elastic net model are primarily associated with CAD- (absence of disease), while positive coefficients are associated with CAD+ (presence of disease). Elastic net regularization combines the penalties of Ridge and LASSO regression, enabling both variable selection and shrinkage, which helps reduce overfitting and isolate the most informative features. Alpha values ranging from 0 (Ridge Regression) to 1 (LASSO regression), in increments of 0.1 were tested. For each value of alpha, one hundred values of lambda were tested using a 10-fold cross-validation to select the optimal lambda based on the smallest mean misclassification error (MME). Coefficients of nonzero features, ie, variables retained in the regularized model, were visualized in a bar plot.

Statistical Analysis

Percentages of populations and MSI of the markers from the manual gating, and the percentage of populations obtained from unsupervised clustering were compared across the study groups using a Kruskal-Wallis test [37]. The P-values were adjusted to control the false discovery rate (FDR) using the Benjamini-Hochberg method [38]. For clusters with a significant Kruskal-Wallis test (FDR-adjusted P < 0.05), Wilcoxon tests [39] between each pair of study groups were performed. The P-value for categorical variables was calculated using Fisher’s exact test.

All statistical analyses were performed in R (version 4.4.2). Boxplots show the medians, the 1st and 3rd quartiles. Whiskers represent the minimal and maximal values after excluding outliers, defined as points beyond 1.5 times the interquartile range (IQR) from the 25th or 75th percentile. Individual data points representing each participant are shown. Plotting unless otherwise stated was done using the ggplot2 R package (version 3.5.1). R packages flowCore (version 2.18.0), FlowSOM (version 2.14.0), glmnet (version 4.1-8), ggpubr (version 0.6.0), umap (version 0.2.10.0), corrplot (version 0.95), and ComplexHeatmap (version 2.22.0) were also employed.

Results

Study Population and Baseline Characteristics by CAD Satus and HIV/TB Status

A total of 61 participants were initially classified according to coronary artery disease (CAD) status as CAD+ or CAD-. Among CAD+ participants, disease burden was further stratified using the segment involvement score (SIS), which quantifies the number of coronary segments containing atherosclerotic plaque. This resulted in three final groups: CAD- (n=40), CAD+ with minimal involvement (CAD+ min; SIS ≤ 2, n=13), and CAD+ with greater involvement (CAD+ >min; SIS > 2, n=8). Sociodemographic characteristics, HIV history, and clinical variables were largely comparable across CAD and CAD/SIS groups. The only significant difference observed was the 10-year ASCVD risk score, which was highest among CAD+ participants overall (Supplementary Table 2) and particularly elevated in the CAD+ >min subgroup (Table 1).

In parallel, we characterized the study population by HIV and TB status, as these are 2 major infectious diseases that could drive monocyte alterations within our study setting. This yielded 5 HIV/TB categories: HIV-TB- (n=13), HIV-TBI+ (n=22), HIV+TB- (n=9), HIV+TBI+ (n=9), HIV+TBpr (n=8). Overall, no major differences were observed in baseline sociodemographic or clinical characteristics across HIV/TB groups (Table 2). Subgroup analyses within the TB categories (TB-, TBI+, and TBpr) and the HIV categories (HIV+ and HIV- similarly revealed no statistically significant baseline characteristics (Supplementary Tables 3 and 4).

Table 1. Demographics and Clinical Characteristics of the Cohort Using 3 CAD/SIS Groups: CAD-, CAD+ min, and CAD+ >min.

Continuous variables are presented as median (interquartile range). P-values for continuous variables were obtained using the Kruskal-Wallis test. The P-value for the only categorical variable was calculated using Fisher’s exact test. P-values were adjusted for false discovery rate (FDR) using the Benjamini-Hochberg (BH) procedure. None of the participants were current smokers.
*n=59

Table 2. Demographics and Clinical Characteristics of the Cohort Classified Using 5 HIV/TB Groups: HIV-TB-, HIV-TBI+, HIV+TB-, HIV+TBI+, HIV+TBpr.

Continuous variables are presented as median (interquartile range). P-values for continuous variables were obtained using the Kruskal-Wallis test, adjusted for false discovery rate (FDR) using the Benjamini-Hochberg (BH) procedure. The P-value for the only categorical variable was calculated using Fisher’s exact test. P-values were adjusted for false discovery rate (FDR) using the Benjamini-Hochberg (BH) procedure. None of the participants were current smokers.
*n=59

Soluble Markers of Monocyte Activation by CAD Status and HIV/TB Status

Soluble monocyte activation markers showed no significant differences across CAD groups (sCD14: P = 0.396 and sCD163: P = 0.951). Interestingly, significant differences emerged when comparing soluble monocyte activation markers, particularly sCD14, by HIV/TB status. Individuals living with HIV had higher sCD14 levels compared to HIV- participants (P = 0.000498, Supplementary Figure 2, top left panel). Within the TB groups, participants in the TBpr group exhibited elevated sCD14 levels relative to both TBI+ and TB- groups (P = 0.002; Supplementary Figure 2, top right panel). No statistically significant differences in sCD163 were observed across HIV and TB groups (Supplementary Figure 2, bottom panels). Additional differences were observed across the combined TB/HIV groups (Supplementary Figure 3), overall showing a gradual increase of sCD14 through the HIV/TB spectrum, with the lowest sCD14 plasma levels in the HIV-TB- group, and the highest sCD14 levels in the HIV+TBpr group. Across the CAD groups, sCD163 and sCD14 levels did not differ significantly (Supplementary Figure 4).

CAD+ >Min Participants Show Lower MSI of CD163 Than CAD- and CAD+ Min in Manually Gated nCM

We then moved to our planned characterization of monocyte populations by CAD status and HIV/TB status. We hypothesized that at least one monocyte subpopulation would differ significantly in abundance across the 5 HIV/TB groups as well as across the CAD groups. It was further hypothesized that such a population would exhibit distinct phenotypic or functional marker expression patterns compared with all other study groups.

Frequencies of different subpopulations and MSI of several markers in TM, CM, IM, and nCM, respectively, based on manual gating, were compared across the CAD-/+ and CAD/SIS groups (Supplementary Figures 5 – 12). Statistically significant differences were observed in the nCM (Supplementary Figure 8, FDR-adjusted P = 0.020). Individuals in the CAD+ >min group had a lower MSI of CD163 compared to both the CAD- and CAD+ min groups (P = 0.023 in both cases). No statistically significant differences in the manually gated populations were observed across the 5 TB/HIV nor in the CAD-/+ groups (Supplementary Figures 9 – 16).

Unsupervised Clustering Identifies Distinct Subpopulations

Unsupervised clustering via FlowSOM was used to separate the total monocytes into 12 clusters. Scaled median marker expressions were used to classify these into 2 CM, 5 IM, and 5 nCM populations (Supplementary Figure 17). The 2 CM populations (clusters 1 and 3) were merged into a single population, as relevant markers such as CD14, CD16, and HLA-DR were similarly expressed. Four IM subpopulations (clusters 4 and 5; and clusters 2 and 9) were also combined into 2 populations. Three nCM populations (clusters 6, 7, and 10) were merged due to similar characteristics. This resulted in a final list of 7 clusters, corresponding to 1 CM, 3 IM, and 3 nCM populations, for downstream analyses (Figure 1). The heatmap of scaled marker expression of each population, along with their percentages in all cells, is shown in Figure 1A. Visualization of the clustering results in a UMAP (Figure 1B) showed that most populations clustered closely together. However, 2 of the nCM populations were slightly (purple population in the plot) or markedly (orange population) separated from the rest.

The CM population comprises the highest percentage of monocytes (74.9%), which is also visually observed in most samples (Figure 1C). Among the 3 IM populations, the HLA-DR++ CD45RA+ CD38+ GPR56- is the largest (4.9%), followed by the CD86dim CX3CR1+ CD45RAdim GPR56- CXCR3- (2.5%) and by HLA-DR+ CD45RA++ CD38++ GPR56++ (0.3%). The nCM portion includes the CD86+ CX3CR1++ CD45RA++ GPR56- CD38- CXCR3+ population (14.1%), which is the largest after the CM. The other 2 nCM populations are CD86- CX3CR1+ CD45RA++ GPR56++ CXCR3dim (3.1%), which is located separately from the other populations in the corner of the UMAP (Figure 1B), possibly due to the lack of expression in CD14 and CD86; and the population CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ (0.1%).

Figure 1. Monocyte subpopulations identification via FlowSOM and UMAP. A) Scaled median marker expressions of the monocyte populations found using FlowSOM. Seven populations (1 classical, 3 intermediate, and 3 non-classical monocytes) were identified after merging similar populations (see Supplementary Figure 13). The bar plot on the right indicates the percentage of the population with respect to all cells, ie, all manually gated monocytes used during clustering, in all the samples. B) Uniform Manifold Approximation and Projections (UMAPs) of 240,000 monocytes randomly selected from the cells used in the unsupervised clustering. The populations obtained from the FlowSOM analysis were overlaid. C) Percentage of the monocyte populations obtained from FlowSOM in each of the samples in the cohort.
CM, classical monocytes; IM intermediate monocytes; nCM, non-classical monocytes

Percentages of 2 Different nCM Populations from the Unsupervised Clustering Differ Between the Different CAD/SIS and TB/HIV Study Groups

Next, the frequencies of the 7 populations derived from unsupervised clustering, calculated as percentages of the total monocyte population used for clustering, were compared across our
study groups.

Hypothesis testing was conducted to identify differences across the CAD/SIS groups (CAD-, CAD+ min, and CAD+ >min groups). The CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ nCM cluster (purple population in the UMAPs, Figure 2A) showed a statistically significant difference in the percentage of cells among the 3 study groups (FDR-adjusted P = 0.033). In post hoc analyses, participants from the CAD+ >min group had lower proportions of this population compared to the CAD- (P = 0.033) and CAD+ min groups (P = 0.022, Figure 2B). CAD- and CAD+ min participants showed comparable proportions of this population among total monocytes. There were no significant differences across the CAD-/+ groups (Supplementary Figure 18).

Then, the proportion of each of the 7 FlowSOM-derived populations was compared across the 5 HIV/TB study groups. In this analysis, a different nCM population showed statistically significant differences between the groups. The nCM cluster characterized as CD86+ CX3CR1++ CD45RA++ GPR56- CXCR3- (yellow population in the UMAPs, Figures 1 and Supplementary Figure 16) was present at significantly higher proportions in the HIV-TBI+ group compared to all the other 4 TB/HIV clinical groups (Figure 3). Interestingly, the same trend was observed in the manual gating results (Supplementary Figure 16), where the HIV-TBI+ group showed the highest median of both percentage of CX3CR1+ and MSI of CX3CR1 in nCM. The trends might be a sign of inflammation in people who are TBI+ , in agreement with previous research [20].

Non-Classical Monocyte CX3CR1 Expression and Subset Composition Correlate With
Cardiovascular Risk

To identify monocyte characteristics associated with CAD, elastic net regression analysis was performed using CAD-/+ status as the predicted variable. The regression model included 32 features (MSI and percent of populations) obtained by manually gating, and the percentages of the 7 populations identified in the FlowSOM clustering.

During model tuning, an alpha value of 0.6 was selected based on evaluation across 100 lambda values per alpha, as this combination yielded the lowest mean misclassification error (Supplementary Figure 20). As a result, the fitted model was equivalent to an elastic net penalization that combines L1 and L2 regularization, shrinking coefficients of less important features to zero and effectively selecting a sparse set of predictors. Figure 4 shows the coefficients of the nonzero features retained in the final model.

The top 3 variables that were associated with CAD+ were the MSI of CX3CR1 and CD86 in nCM, and the percent of IM CD14+ CD16+. The top 3 that were associated with CAD- were the percent of nCM CD86+ CX3CR1++ CD45RA++ GPR56- CD38- CXCR3- identified in the unsupervised clustering, the percent of TM CD14+ CD163+, and the MSI of HLA-DR in nCM. The TM CD14+ CD163 agrees with the results obtained from the unsupervised clustering shown in Figure 2.

Figure 2. Differences in monocyte subpopulations obtained with FlowSOM across CAD/SIS Groups. A) Uniform Manifold Approximation and Projections (UMAPs) stratified by CAD/SIS study groups. 80000 cells were randomly selected per CAD/SIS group, and clusters from the FlowSOM analysis were overlaid. Contour lines are of all cells sampled. B) Statistical differences of percentages of cells for each monocyte population obtained from unsupervised clustering by FlowSOM across the three CAD/SIS study groups. Kruskal-Wallis (KW) tests were performed as an omnibus test to assess overall differences across the groups and P-values were adjusted for False Discovery Rate (FDR) using the Benjamini-Hochberg method. The population nCM CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ was statistically significant (FDR-adjusted P = 0.033) and was further analyzed by post-hoc Wilcoxon tests across the 3 groups. The percentage in the CAD+ >min study group was statistically lower than the CAD+ min (P = 0.022) and CAD- (P = 0.033).

Figure 3. Differences in monocyte subpopulations obtained with FlowSOM across HIV/TB Groups. Statistical differences of percentages of cells for each monocyte population obtained from unsupervised clustering by FlowSOM across the 5 HIV/TB study groups. Kruskal-Wallis (KW) tests were performed as an omnibus test and P-values were adjusted for False Discovery Rate (FDR) using Benjamini-Hochberg method. Among the populations tested, only nCM CD86+ CX3CR1++ CD45RA++ GPR56- CD38- CXCR3- reached statistical significance (FDR-adjusted P < 0.05) and was further subjected to post-hoc Wilcoxon tests across the 5 groups. The small colored rectangles on top left corners follow the colors used in the UMAPs in Figures 1 and 2.

Figure 4. Monocyte features associated with CAD status via regularized regression. Nonzero features obtained from elastic net regression and their coefficients associated with the CAD+/- status. The coefficients were obtained after testing alpha values from 0 to 1 in 0.1 increments. The alpha value that led to the smallest MSE was 0.6 (elastic net regularization). Manual gating and unsupervised clustering variables were used as features. (UC) indicates that the variable was obtained from unsupervised clustering. A total of 16 variables were retained as important predictors of the CAD status after regularization. Seven and nine variables were positively and negatively associated, respectively.

Discussion

In this study, mass cytometry and advanced statistical analysis were employed to perform immunophenotyping of monocyte populations and study their associations with infectious and inflammatory diseases, specifically HIV, tuberculosis, and cardiovascular disease. Elastic net regression to identify key features that are most strongly correlated with CAD-/+ status was also applied.

The higher sCD14 levels in individuals with HIV likely reflect chronic monocyte activation driven by microbial translocation and systemic inflammation. The stepwise increase across TB categories suggests that M. tuberculosis exposure further amplifies this activation. The combined gradient (HIV+TBpr > HIV+TBI+ > HIV+ TB- > HIV-TBI+ > HIV-TB-) points to additive effects of HIV and TB, supporting sCD14 as an indicator of cumulative innate immune activation in co-endemic settings [40].

The analysis of the manually gated data showed a statistically significant difference across the CAD/SIS groups regarding the median signal intensity of CD163 in nCM (Figure S8). CD163 has previously been shown to be a potential predictor of mortality and heart failure in people with CAD [41]. Previous studies have indicated that CD163 is characteristic of anti-inflammatory monocytes and has been identified as a marker related to CAD [42]. The downregulation of CD163 MSI in nCMs among individuals with CAD may occur under specific conditions. For instance, inflammation could promote CD163 shedding, producing sCD163 while reducing detectable CD163 surface expression [43]. Furthermore, in vitro studies show CD163 downregulation in monocytes and macrophages has been associated with inflammation mediators such as TNF-a, IFN-g, and LPS [44].

When comparing the populations obtained from the unsupervised clustering across the CAD/SIS groups, the nCM population CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ was basically absent in the CAD+ >min group. CD86 is a co-stimulatory molecule and plays a crucial role in T cell activation by binding to CD28 on T cells, providing the necessary second signal for full T cell activation. CD86 expression suggests that the subset might be involved in interaction with T cells, likely promoting an immune response [45]. CX3CR1 is a chemokine receptor predominantly expressed on nCMs, which actively patrol the vascular endothelium under both homeostatic and inflammatory conditions [46]. CD45RA is a characteristic marker of naïve T cells and is less well-characterized in monocyte populations [47]. GPR56 is an adhesion protein-coupled receptor that primarily functions as an inhibitor in natural killer cells. Studies involving GPR56+ monocytes are limited. However, due to the similarities between GPR56 and EMR2, it is likely that the function of GPR56+ is similar to EMR2+ monocytes that are involved in the activation and migration of innate cells during systemic inflammation responses [48]. CXCR3 is another chemokine receptor present in monocytes that can infiltrate CAD-related lesions. CXCR3+ monocytes are associated with plaque instability and adverse cardiovascular events [49]. It is possible that depletion of this nCM population with low expression of CD86 and CX3CR1 could be an indication of increased relative abundance of highly expressing CD86 and CX3CR1 nCM in the circulation of people with CAD. Alternatively, the lower percentage of nCM with low expression of CD86 and CX3CR1 in people with CAD might indicate that most of these cells may have already migrated to the tissues and are no longer in the bloodstream, thus contributing to atherosclerotic lesion development.

Across the TB/HIV groups, the HIV-TBI+ group exhibited a greater percentage of nCM co-expressing CD86+ and CX3CR1++ compared to the other groups. The CX3CR1 upregulation suggests enhanced chemotactic potential and tissue patrolling capacity, possibly reflecting a more activated, vigilant monocyte phenotype even in latent infection [6]. Further, a recent study in mice infected with TB showed that the absence of chemokine receptors CX3CR1 was associated with an altered positioning of derived dendritic cells in mediastinal lymph nodes [50]. The expression of CD86 suggests a role of this population in antigen presentation. Previously it has been reported that circulating nCMs from HIV-TBI+ individuals exhibit a higher production of interleukin (IL)-6 and tumor necrosis factor (TNF)-α pro-inflammatory cytokines, compared to HIV-TBI- control individuals [6]. Future studies can define the functional capacity of this newly classified CD86+ CX3CR1++ CXCR3+ nCM subpopulation and their role in TB infection control and persistent inflammation.

It is noteworthy that in the HIV–TBI+ group, we found an enrichment of nCM subsets highly expressing CX3CR1 and CD86, whereas the CAD+>min group had a significant decrease of nCMs with low CD86 and CX3CR1 expression. These 2 nCM groups appear to reflect distinct biological trajectories. In the HIV-TBI+ group, the nCM subset shows elevated CX3CR1 and CD86 expression, consistent with a patrolling/monitoring phenotype. In contrast, the nCM CD86dim CX3CR1dim CD45RA+ GPR56+ CXCR3+ cluster is significantly reduced in individuals in the CAD+ >min group. This depletion suggests that CAD may be associated with lower circulating frequencies of CD86dim CX3CR1dim nCMs, potentially reflecting accelerated recruitment, activation, or apoptosis of this patrolling subset in the setting of chronic endothelial injury and systemic inflammation [51–53]. Given that CAD is characterized by active recruitment of CX3CR1+ CD86+ nCMs to inflamed vasculature, the CD86dim CX3CR1dim subset may also be mobilized out of circulation. Moreover, the CAD inflammatory milieu is known to shift monocytes toward more activated phenotypes (eg, CD86+) while selectively depleting less-activated populations [52]. This finding aligns with prior observations that chronic inflammatory diseases, including CAD, remodel hematopoiesis and skew circulating myeloid cell distributions. Collectively, these findings suggest that CX3CR1- and CD86-expressing nCM subsets participate in shared inflammatory pathways that may actively contribute to disease pathogenesis in both chronic infectious and cardiovascular settings. In CAD, depletion of circulating CD86^dim CX3CR1^dim nCMs may reflect enhanced recruitment of these cells to inflamed endothelium, where CX3CR1-mediated adhesion, antigen presentation, and local cytokine production could promote vascular inflammation and atherosclerotic lesion development. In contrast, enrichment of CX3CR1^+ CD86^+ non-classical monocytes in HIV-TBI+ individuals suggests a sustained patrolling and antigen-presenting phenotype that may contribute to persistent immune activation and immunopathology during latent infection. Further validation in longitudinal cohorts will be required to determine whether these subsets can inform disease monitoring, prognosis, or targeted immunomodulatory strategies.

The nonzero features obtained in the elastic net regression provided information about the relevant features associated with the CAD-/+ status. For instance, among the variables strongly associated with CAD+ status, the MSI of CX3XR1 in nCM had the coefficient with the greatest absolute value, which is an indication of a strong immune response and a sign of inflammation in HIV/TB coinfections [54]. The next positive coefficient was the MSI of CD86 in nCM. It has been shown that expression of CD86 increases as CAD progresses [55]. The CD86+ nCM are potentially geared toward antigen presentation and T-cell stimulation [56]. The percentage of IMs was also positively associated with CAD+. IM are pro-inflammatory and atherogenic [56]. Also, IMs are associated with increased carotid intima-media thickness, a marker of atherosclerosis, and predicted cardiovascular risk even in general population cohorts [57]. Among the features that were identified as important predictors of CAD-, the nCM population CD86+ CX3CR1++ CD45RA++ GPR56- CD38- CXCR3- was the most strongly inversely associated with CAD. Moderate CD86 expression in nCMs suggests a role in immune surveillance, rather than driving full activation and inflammation [58]. In general, CX3CR1 is strongly expressed on nCMs. CX3CR1+ nCM typically patrol the endothelium and contribute to tissue repair and vascular stability [46, 59]. The other 2 variables associated with CAD- are the percentage of TM CD14+ CD163+ and the MSI of HLA-DR in nCM. The first of these findings suggests that people in the CAD- group have higher percentages of CD163+monocytes. These cells are likely anti-inflammatory [60]. Angiotensin or context-specific studies link high expression of HLA-DR to monocytes with enhanced antigen presentation and immune surveillance traits, especially in nCMs. One study in HIV populations found that CCR9dim HLA-DR+ nCM were associated with absence of subclinical atherosclerosis, suggesting their HLA-DR+ phenotype is protective [61].

There are several limitations in this study. For instance, the total number of participants in the sample was small for some of the study groups. Therefore, there might not have been enough power to detect subtle differences in monocyte biomarkers between all groups. Moreover, our cohort includes individuals who live with multiple comorbid conditions, making the analysis even more challenging, a common limitation of human-based clinical and translational studies. The limited sample size does not allow for a detailed analysis or adjustment of all combinations of conditions simultaneously. Still, our identification of specific monocyte subpopulations associated with TB and CAD can be further investigated and validated in future studies. The unsupervised clustering analysis was performed on the manually gated monocyte population, which limited the number of cells that could be included per sample. We randomly selected 10,000 cells per sample or used all cells if the sample contained less than 10,000 cells. This was the case for 15 samples, and the minimum number of cells in a sample was 1,333. To address potential variability based on the limited sample size, we performed repeated runs of the unsupervised clustering.

The distinct monocyte subpopulations found in this study deserve more study to understand how these populations interact with the diseases as they evolve. This study provides a framework to investigate the functionality of unusual monocyte populations in diseases that affect a large percentage of the human population.

Acknowledgments

Thanks to the study participants.

Funding

This work was in part supported by the National Heart, Lung, and Blood Institute (grant number 1R01HL156779 to MAH, K23 HL123341 to CTL, K01 HL147723 to TMT, K23 HL156759 to CNW), the National Center for Advancing Translational Sciences (R03TR004097 to MAH), and the National Institute of Arthritis and Musculoskeletal and Skin Diseases (P30AR070549 to SA) at the National Institutes of Health, Doris Duke CSDA 2021193 (CNW), Burroughs Wellcome Fund 1021480 (CNW).

Disclaimer

The contents are solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health or the institutions with which the authors are affiliated. The funding source had no role in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to submit the report for publication.

Potential Conflicts of Interest

The authors report no relevant conflicts of interest related to this article.

Supplementary Data

Supplementary materials are available at the Pathogens and Immunity website. Supplementary data may be provided by the authors to benefit the reader. Supplementary data are not copyedited and are the sole responsibility of the authors. Questions or comments related to supplementary materials should be addressed to the corresponding author.

Supplementary Tables and Figures

References

1. Govender RD, Hashim MJ, Khan MA, Mustafa H, Khan G. Global Epidemiology of Hiv/Aids: A Resurgence in North America and Europe. J Epidemiol Glob Health. 2021;11(3):296-301. doi: 10.2991/jegh.k.210621.001. PubMed PMID: 34270183; PMCID: PMC8435868.

2. Berida T, Lindsley CW. Move over Covid, Tuberculosis Is Once Again the Leading Cause of Death from a Single Infectious Disease. J Med Chem. 2024;67(24):21633-40. doi: 10.1021/acs.jmedchem.4c02876. PubMed PMID: 39652566.

3. Gaidai O, Cao Y, Loginov S. Global Cardiovascular Diseases Death Rate Prediction. Curr Probl Cardiol. 2023;48(5):101622. doi: 10.1016/j.cpcardiol.2023.101622. PubMed PMID: 36724816.

4. Siedner MJ, Ghoshhajra B, Erem G, Nassanga R, Randhawa M, Ochieng A, Acan M, Lu MT, Thondapu V, Takigami A, Reynolds Z, Atwiine F, Tindimwebwa E, Gilbert RF, Passell E, Sagar S, Tong Y, Ntusi NAB, Tsai AC, Bibangambah P, Gaziano T, Hoeppner SS, Longenecker CT, Okello S, Asiimwe S. Epidemiology of Coronary Atherosclerosis among People Living with Hiv in Uganda : A Cross-Sectional Study. Ann Intern Med. 2025;178(4):468-78. doi: 10.7326/annals-24-02233. PubMed PMID: 40073231; PMCID: PMC12022967.

5. Hyle EP, Mayosi BM, Middelkoop K, Mosepele M, Martey EB, Walensky RP, Bekker LG, Triant VA. The Association between Hiv and Atherosclerotic Cardiovascular Disease in Sub-Saharan Africa: A Systematic Review. BMC Public Health. 2017;17(1):954. doi: 10.1186/s12889-017-4940-1. PubMed PMID: 29246206; PMCID: PMC5732372.

6. Feria MG, Chang C, Ticona E, Moussa A, Zhang B, Ballena I, Azañero R, Ticona C, De Cecco CN, Fichtenbaum CJ, O’Donnell RE, La Rosa A, Sanchez J, Andorf S, Atehortua L, Katz JD, Chougnet CA, Deepe GS, Jr., Huaman MA. Pro-Inflammatory Alterations of Circulating Monocytes in Latent Tuberculosis Infection. Open Forum Infect Dis. 2022;9(12):ofac629. doi: 10.1093/ofid/ofac629. PubMed PMID: 36570965; PMCID: PMC9772871.

7. Okello S, Amir A, Bloomfield GS, Kentoffio K, Lugobe HM, Reynolds Z, Magodoro IM, North CM, Okello E, Peck R, Siedner MJ. Prevention of Cardiovascular Disease among People Living with Hiv in Sub-Saharan Africa. Prog Cardiovasc Dis. 2020;63(2):149-59. doi: 10.1016/j.pcad.2020.02.004. PubMed PMID: 32035126; PMCID: PMC7237320.

8. Obare LM, Temu T, Mallal SA, Wanjalla CN. Inflammation in Hiv and Its Impact on Atherosclerotic Cardiovascular Disease. Circ Res. 2024;134(11):1515-45. doi: 10.1161/circresaha.124.323891. PubMed PMID: 38781301; PMCID: PMC11122788.

9. Rambaran S, Maseko TG, Lewis L, Hassan-Moosa R, Archary D, Ngcapu S, Garrett N, McKinnon LR, Padayatchi N, Naidoo K, Sivro A. Blood Monocyte and Dendritic Cell Profiles among People Living with Hiv with Mycobacterium Tuberculosis Co-Infection. BMC Immunol. 2023;24(1):21. doi: 10.1186/s12865-023-00558-z. PubMed PMID: 37480005; PMCID: PMC10362598.

10. Huaman MA, Henson D, Ticona E, Sterling TR, Garvy BA. Tuberculosis and Cardiovascular Disease: Linking the Epidemics. Trop Dis Travel Med Vaccines. 2015;1. doi: 10.1186/s40794-015-0014-5. PubMed PMID: 26835156; PMCID: PMC4729377.

11. Chidambaram V, Ruelas Castillo J, Kumar A, Wei J, Wang S, Majella MG, Gupte A, Wang JY, Karakousis PC. The Association of Atherosclerotic Cardiovascular Disease and Statin Use with Inflammation and Treatment Outcomes in Tuberculosis. Sci Rep. 2021;11(1):15283. doi: 10.1038/s41598-021-94590-x. PubMed PMID: 34315941; PMCID: PMC8316554.

12. Knudsen AD, Bouazzi R, Afzal S, Gelpi M, Benfield T, Høgh J, Thomsen MT, Trøseid M, Nordestgaard BG, Nielsen SD. Monocyte Count and Soluble Markers of Monocyte Activation in People Living with Hiv and Uninfected Controls. BMC Infect Dis. 2022;22(1):451. doi: 10.1186/s12879-022-07450-y. PubMed PMID: 35546661; PMCID: PMC9097376.

13. Wallis ZK, Williams KC. Monocytes in Hiv and Siv Infection and Aging: Implications for Inflamm-Aging and Accelerated Aging. Viruses. 2022;14(2). doi: 10.3390/v14020409. PubMed PMID: 35216002; PMCID: PMC8880456.

14. Lyu M, Xu G, Zhou J, Reboud J, Wang Y, Lai H, Chen Y, Zhou Y, Zhu G, Cooper JM, Ying B. Single-Cell Sequencing Reveals Functional Alterations in Tuberculosis. Adv Sci (Weinh). 2024;11(11):e2305592. doi: 10.1002/advs.202305592. PubMed PMID: 38192178; PMCID: PMC10953544.

15. Ji H, Li Y, Fan Z, Zuo B, Jian X, Li L, Liu T. Monocyte/Lymphocyte Ratio Predicts the Severity of Coronary Artery Disease: A Syntax Score Assessment. BMC Cardiovasc Disord. 2017;17(1):90. doi: 10.1186/s12872-017-0507-4. PubMed PMID: 28359298; PMCID: PMC5374608.

16. Moroni F, Ammirati E, Norata GD, Magnoni M, Camici PG. The Role of Monocytes and Macrophages in Human Atherosclerosis, Plaque Neoangiogenesis, and Atherothrombosis. Mediators Inflamm. 2019;2019:7434376. doi: 10.1155/2019/7434376. PubMed PMID: 31089324; PMCID: PMC6476044.

17. Tapp LD, Shantsila E, Wrigley BJ, Pamukcu B, Lip GY. The Cd14++Cd16+ Monocyte Subset and Monocyte-Platelet Interactions in Patients with St-Elevation Myocardial Infarction. J Thromb Haemost. 2012;10(7):1231-41. doi: 10.1111/j.1538-7836.2011.04603.x. PubMed PMID: 22212813.

18. Eligini S, Cosentino N, Fiorelli S, Fabbiocchi F, Niccoli G, Refaat H, Camera M, Calligaris G, De Martini S, Bonomi A, Veglia F, Fracassi F, Crea F, Marenzi G, Tremoli E. Biological Profile of Monocyte-Derived Macrophages in Coronary Heart Disease Patients: Implications for Plaque Morphology. Sci Rep. 2019;9(1):8680. doi: 10.1038/s41598-019-44847-3. PubMed PMID: 31213640; PMCID: PMC6581961.

19. Amengual J, Barrett TJ. Monocytes and Macrophages in Atherogenesis. Curr Opin Lipidol. 2019;30(5):401-8. doi: 10.1097/mol.0000000000000634. PubMed PMID: 31361625; PMCID: PMC7809604.

20. Sampath P, Moideen K, Ranganathan UD, Bethunaickan R. Monocyte Subsets: Phenotypes and Function in Tuberculosis Infection. Front Immunol. 2018;9:1726. doi: 10.3389/fimmu.2018.01726. PubMed PMID: 30105020; PMCID: PMC6077267.

21. Ma R, Yang W, Guo W, Zhang H, Wang Z, Ge Z. Single-Cell Transcriptome Analysis Reveals the Dysregulated Monocyte State Associated with Tuberculosis Progression. BMC Infect Dis. 2025;25(1):210. doi: 10.1186/s12879-025-10612-3. PubMed PMID: 39939918; PMCID: PMC11823163.

22. Bai R, Li Z, Lv S, Wang R, Hua W, Wu H, Dai L. Persistent Inflammation and Non-Aids Comorbidities during Art: Coming of the Age of Monocytes. Front Immunol. 2022;13:820480. doi: 10.3389/fimmu.2022.820480. PubMed PMID: 35479083; PMCID: PMC9035604.

23. Wong ME, Jaworowski A, Hearps AC. The Hiv Reservoir in Monocytes and Macrophages. Front Immunol. 2019;10:1435. doi: 10.3389/fimmu.2019.01435. PubMed PMID: 31297114; PMCID: PMC6607932.

24. Williams DW, Veenstra M, Gaskill PJ, Morgello S, Calderon TM, Berman JW. Monocytes Mediate Hiv Neuropathogenesis: Mechanisms That Contribute to Hiv Associated Neurocognitive Disorders. Curr HIV Res. 2014;12(2):85-96. doi: 10.2174/1570162x12666140526114526. PubMed PMID: 24862333; PMCID: PMC4351961.

25. Ellery PJ, Tippett E, Chiu YL, Paukovics G, Cameron PU, Solomon A, Lewin SR, Gorry PR, Jaworowski A, Greene WC, Sonza S, Crowe SM. The Cd16+ Monocyte Subset Is More Permissive to Infection and Preferentially Harbors Hiv-1 in Vivo. J Immunol. 2007;178(10):6581-9. doi: 10.4049/jimmunol.178.10.6581. PubMed PMID: 17475889.

26. Longenecker CT, Bogorodskaya M, Margevicius S, Nazzinda R, Bittencourt MS, Erem G, Nalukwago S, Huaman MA, Ghoshhajra BB, Siedner MJ, Juchnowski SM, Zidar DA, McComsey GA, Kityo C. Sex Modifies the Association between Hiv and Coronary Artery Disease among Older Adults in Uganda. J Int AIDS Soc. 2022;25(1):e25868. doi: 10.1002/jia2.25868. PubMed PMID: 34995413; PMCID: PMC8741262.

27. Goff DC, Jr., Lloyd-Jones DM, Bennett G, Coady S, D’Agostino RB, Gibbons R, Greenland P, Lackland DT, Levy D, O’Donnell CJ, Robinson JG, Schwartz JS, Shero ST, Smith SC, Jr., Sorlie P, Stone NJ, Wilson PW, Jordan HS, Nevo L, Wnek J, Anderson JL, Halperin JL, Albert NM, Bozkurt B, Brindis RG, Curtis LH, DeMets D, Hochman JS, Kovacs RJ, Ohman EM, Pressler SJ, Sellke FW, Shen WK, Smith SC, Jr., Tomaselli GF. 2013 Acc/Aha Guideline on the Assessment of Cardiovascular Risk: A Report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines. Circulation. 2014;129(25 Suppl 2):S49-73. doi: 10.1161/01.cir.0000437741.48606.98. PubMed PMID: 24222018.

28. Alencherry B, Erem G, Mirembe G, Ssinabulya I, Yun CH, Hung CL, Siedner MJ, Bittencourt M, Kityo C, McComsey GA, Longenecker CT. Coronary Artery Calcium, Hiv and Inflammation in Uganda Compared with the USA. Open Heart. 2019;6(1):e001046. doi: 10.1136/openhrt-2019-001046. PubMed PMID: 31218009; PMCID: PMC6546194.

29. Min JK, Shaw LJ, Devereux RB, Okin PM, Weinsaft JW, Russo DJ, Lippolis NJ, Berman DS, Callister TQ. Prognostic Value of Multidetector Coronary Computed Tomographic Angiography for Prediction of All-Cause Mortality. J Am Coll Cardiol. 2007;50(12):1161-70. doi: 10.1016/j.jacc.2007.03.067. PubMed PMID: 17868808.

30. Hahne F, LeMeur N, Brinkman RR, Ellis B, Haaland P, Sarkar D, Spidlen J, Strain E, Gentleman R. Flowcore: A Bioconductor Package for High Throughput Flow Cytometry. BMC Bioinformatics. 2009;10:106. doi: 10.1186/1471-2105-10-106. PubMed PMID: 19358741; PMCID: PMC2684747.

31. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J. Bioconductor: Open Software Development for Computational Biology and Bioinformatics. Genome Biol. 2004;5(10):R80. doi: 10.1186/gb-2004-5-10-r80. PubMed PMID: 15461798; PMCID: PMC545600.

32. Bendall SC, Simonds EF, Qiu P, Amir el AD, Krutzik PO, Finck R, Bruggner RV, Melamed R, Trejo A, Ornatsky OI, Balderas RS, Plevritis SK, Sachs K, Pe’er D, Tanner SD, Nolan GP. Single-Cell Mass Cytometry of Differential Immune and Drug Responses across a Human Hematopoietic Continuum. Science. 2011;332(6030):687-96. doi: 10.1126/science.1198704. PubMed PMID: 21551058; PMCID: PMC3273988.

33. Van Gassen S, Callebaut B, Van Helden MJ, Lambrecht BN, Demeester P, Dhaene T, Saeys Y. Flowsom: Using Self-Organizing Maps for Visualization and Interpretation of Cytometry Data. Cytometry A. 2015;87(7):636-45. doi: 10.1002/cyto.a.22625. PubMed PMID: 25573116.

34. Nowicka M, Krieg C, Crowell HL, Weber LM, Hartmann FJ, Guglietta S, Becher B, Levesque MP, Robinson MD. Cytof Workflow: Differential Discovery in High-Throughput High-Dimensional Cytometry Datasets. F1000Res. 2017;6:748. doi: 10.12688/f1000research.11622.3. PubMed PMID: 28663787; PMCID: PMC5473464.

35. Leland McInnes JH, James Melville. Umap: Uniform Manifold Approximation and Projection for Dimension Reduction [website]. arXiv2020. Available from: https://arxiv.org/abs/1802.03426.

36. Zou H, Hastie T. Regularization and Variable Selection Via the Elastic Net. Journal of the Royal Statistical Society Series B: Statistical Methodology. 2005;67(2):301-20. doi: 10.1111/j.1467-9868.2005.00503.x.

37. William H. Kruskal WAW. Use of Ranks in One-Criterion Variance Analysis [Website]. Journal of the American Statistical Association1952. Available from: https://webspace.ship.edu/pgmarr/Geo441/Readings/Kruskal%20and%20Wallis%201952%20-%20Use%20of%20Ranks%20in%20One-Criterion%20Variance%20Analysis.pdf.

38. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society: Series B (Methodological). 1995;57(1):289-300. doi: https://doi.org/10.1111/j.2517-6161.1995.tb02031.x.

39. Wilcoxon F. Individual Comparisons by Ranking Methods. In: Kotz S, Johnson NL, editors. Breakthroughs in Statistics: Methodology and Distribution. New York, NY: Springer New York; 1992. p. 196-202.

40. Lawn SD, Labeta MO, Arias M, Acheampong JW, Griffin GE. Elevated Serum Concentrations of Soluble Cd14 in Hiv- and Hiv+ Patients with Tuberculosis in Africa: Prolonged Elevation during Anti-Tuberculosis Treatment. Clin Exp Immunol. 2000;120(3):483-7. doi: 10.1046/j.1365-2249.2000.01246.x. PubMed PMID: 10844527; PMCID: PMC1905566.

41. Durda P, Raffield LM, Lange EM, Olson NC, Jenny NS, Cushman M, Deichgraeber P, Grarup N, Jonsson A, Hansen T, Mychaleckyj JC, Psaty BM, Reiner AP, Tracy RP, Lange LA. Circulating Soluble Cd163, Associations with Cardiovascular Outcomes and Mortality, and Identification of Genetic Variants in Older Individuals: The Cardiovascular Health Study. J Am Heart Assoc. 2022;11(21):e024374. doi: 10.1161/jaha.121.024374. PubMed PMID: 36314488; PMCID: PMC9673628.

42. Aristoteli LP, Møller HJ, Bailey B, Moestrup SK, Kritharides L. The Monocytic Lineage Specific Soluble Cd163 Is a Plasma Marker of Coronary Atherosclerosis. Atherosclerosis. 2006;184(2):342-7. doi: 10.1016/j.atherosclerosis.2005.05.004. PubMed PMID: 15979079.

43. Moreno JA, Ortega-Gómez A, Delbosc S, Beaufort N, Sorbets E, Louedec L, Esposito-Farèse M, Tubach F, Nicoletti A, Steg PG, Michel JB, Feldman L, Meilhac O. In Vitro and in Vivo Evidence for the Role of Elastase Shedding of Cd163 in Human Atherothrombosis. Eur Heart J. 2012;33(2):252-63. doi: 10.1093/eurheartj/ehr123. PubMed PMID: 21606088.

44. Etzerodt A, Moestrup SK. Cd163 and Inflammation: Biological, Diagnostic, and Therapeutic Aspects. Antioxid Redox Signal. 2013;18(17):2352-63. doi: 10.1089/ars.2012.4834. PubMed PMID: 22900885; PMCID: PMC3638564.

45. Subauste CS, de Waal Malefyt R, Fuh F. Role of Cd80 (B7.1) and Cd86 (B7.2) in the Immune Response to an Intracellular Pathogen1. The Journal of Immunology. 1998;160(4):1831-40. doi: 10.4049/jimmunol.160.4.1831.

46. Thomas G, Tacke R, Hedrick CC, Hanna RN. Nonclassical Patrolling Monocyte Function in the Vasculature. Arterioscler Thromb Vasc Biol. 2015;35(6):1306-16. doi: 10.1161/atvbaha.114.304650. PubMed PMID: 25838429; PMCID: PMC4441550.

47. Tahir S, Steffens S. Nonclassical Monocytes in Cardiovascular Physiology and Disease. Am J Physiol Cell Physiol. 2021;320(5):C761-c70. doi: 10.1152/ajpcell.00326.2020. PubMed PMID: 33596150.

48. Yona S, Lin HH, Dri P, Davies JQ, Hayhoe RP, Lewis SM, Heinsbroek SE, Brown KA, Perretti M, Hamann J, Treacher DF, Gordon S, Stacey M. Ligation of the Adhesion-Gpcr Emr2 Regulates Human Neutrophil Function. Faseb j. 2008;22(3):741-51. doi: 10.1096/fj.07-9435com. PubMed PMID: 17928360.

49. Altara R, Mallat Z, Booz GW, Zouein FA. The Cxcl10/Cxcr3 Axis and Cardiac Inflammation: Implications for Immunotherapy to Treat Infectious and Noninfectious Diseases of the Heart. J Immunol Res. 2016;2016:4396368. doi: 10.1155/2016/4396368. PubMed PMID: 27795961; PMCID: PMC5066021.

50. Mohapatra A, Howard Z, Ernst JD. Ccr2 Recruits Monocytes to the Lung, While Cx3cr1 Modulates Positioning of Monocyte-Derived Cd11c (Pos) Cells in the Lymph Node during Pulmonary Tuberculosis. bioRxiv. 2025. doi: 10.1101/2025.02.07.637199. PubMed PMID: 39974908; PMCID: PMC11839135.

51. Landsman L, Bar-On L, Zernecke A, Kim KW, Krauthgamer R, Shagdarsuren E, Lira SA, Weissman IL, Weber C, Jung S. Cx3cr1 Is Required for Monocyte Homeostasis and Atherogenesis by Promoting Cell Survival. Blood. 2009;113(4):963-72. doi: 10.1182/blood-2008-07-170787. PubMed PMID: 18971423.

52. Chavakis T, Mitroulis I, Hajishengallis G. Hematopoietic Progenitor Cells as Integrative Hubs for Adaptation to and Fine-Tuning of Inflammation. Nat Immunol. 2019;20(7):802-11. doi: 10.1038/s41590-019-0402-5. PubMed PMID: 31213716; PMCID: PMC6709414.

53. Mildner A, Kim KW, Yona S. Unravelling Monocyte Functions: From the Guardians of Health to the Regulators of Disease. Discov Immunol. 2024;3(1):kyae014. doi: 10.1093/discim/kyae014. PubMed PMID: 39430099; PMCID: PMC11486918.

54. Guo N, Chen Y, Su B, Yang X, Zhang Q, Song T, Wu H, Liu C, Liu L, Zhang T. Alterations of Ccr2 and Cx3cr1 on Three Monocyte Subsets during Hiv-1/Treponema Pallidum Coinfection. Front Med (Lausanne). 2020;7:272. doi: 10.3389/fmed.2020.00272. PubMed PMID: 32626718; PMCID: PMC7314900.

55. Sbrana S, Campolo J, Clemente A, Bastiani L, Cecchettini A, Ceccherini E, Caselli C, Neglia D, Parodi O, Chiappino D, Smit JM, Scholte AJ, Pelosi G, Rocchiccioli S. Blood Monocyte Phenotype Fingerprint of Stable Coronary Artery Disease: A Cross-Sectional Substudy of Smartool Clinical Trial. Biomed Res Int. 2020;2020:8748934. doi: 10.1155/2020/8748934. PubMed PMID: 32802883; PMCID: PMC7403909.

56. Patel VK, Williams H, Li SCH, Fletcher JP, Medbury HJ. Monocyte Inflammatory Profile Is Specific for Individuals and Associated with Altered Blood Lipid Levels. Atherosclerosis. 2017;263:15-23. doi: 10.1016/j.atherosclerosis.2017.05.026. PubMed PMID: 28570862.

57. SahBandar IN, Ndhlovu LC, Saiki K, Kohorn LB, Peterson MM, D’Antoni ML, Shiramizu B, Shikuma CM, Chow DC. Relationship between Circulating Inflammatory Monocytes and Cardiovascular Disease Measures of Carotid Intimal Thickness. J Atheroscler Thromb. 2020;27(5):441-8. doi: 10.5551/jat.49791. PubMed PMID: 31588100; PMCID: PMC7242227.

58. Padgett LE, Araujo DJ, Hedrick CC, Olingy CE. Functional Crosstalk between T Cells and Monocytes in Cancer and Atherosclerosis. J Leukoc Biol. 2020;108(1):297-308. doi: 10.1002/jlb.1mir0420-076r. PubMed PMID: 32531833; PMCID: PMC8006924.

59. Kim KW, Ivanov S, Williams JW. Monocyte Recruitment, Specification, and Function in Atherosclerosis. Cells. 2020;10(1). doi: 10.3390/cells10010015. PubMed PMID: 33374145; PMCID: PMC7823291.

60. Li B, Shaikh F, Zamzam A, Abdin R, Qadura M. Inflammatory Biomarkers to Predict Major Adverse Cardiovascular Events in Patients with Carotid Artery Stenosis. Medicina (Kaunas). 2024;60(6). doi: 10.3390/medicina60060997. PubMed PMID: 38929614; PMCID: PMC11205582.

61. Wiche Salinas TR, Zhang Y, Gosselin A, Rosario NF, El-Far M, Filali-Mouhim A, Routy JP, Chartrand-Lefebvre C, Landay AL, Durand M, Tremblay CL, Ancuta P. Alterations in Th17 Cells and Non-Classical Monocytes as a Signature of Subclinical Coronary Artery Atherosclerosis during Art-Treated Hiv-1 Infection. Cells. 2024;13(2). doi: 10.3390/cells13020157. PubMed PMID: 38247848; PMCID: PMC10813976.

Footnotes

Submitted December 22, 2025 | Accepted January 2, 2026 | Published January 29, 2026

Copyright

Copyright © 2026 The Authors. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International License.