Combined protein and transcript single-cell RNA sequencing in human peripheral blood mononuclear cells

Abstract

Background

Cryopreserved peripheral blood mononuclear cells (PBMCs) are frequently collected and provide disease- and treatment-relevant data in clinical studies. Here, we developed combined protein (40 antibodies) and transcript single-cell (sc)RNA sequencing (scRNA-seq) in PBMCs.

Results

Among 31 participants in the Women’s Interagency HIV Study (WIHS), we sequenced 41,611 cells. Using Boolean gating followed by Seurat UMAPs (tool for visualizing high-dimensional data) and Louvain clustering, we identified 50 subsets among CD4+ T, CD8+ T, B, NK cells, and monocytes. This resolution was superior to flow cytometry, mass cytometry, or scRNA-seq without antibodies. Combined protein and transcript scRNA-seq allowed for the assessment of disease-related changes in transcriptomes and cell type proportions. As a proof-of-concept, we showed such differences between healthy and matched individuals living with HIV with and without cardiovascular disease.

Conclusions

In conclusion, combined protein and transcript scRNA sequencing is a suitable and powerful method for clinical investigations using PBMCs.

Background

PBMCs are a rich source of disease- and treatment-relevant information [6, 19, 36, 48, 57, 76, 84,85,86]. PBMCs can be analyzed without mechanical or enzymatic dissociation, which are known to alter cell surface markers and transcriptomes [79]. PBMCs can be cryopreserved without loss of viability. At the most basic level, lymphocytes and monocyte can be distinguished by morphology using automated cell counters (blood cell counter, CBC) [7]. The current practice is to use flow cytometry of between 8 and 30 markers simultaneously [41, 46, 51, 60, 64] or even 43 markers with spectral cytometry [66]. Mass cytometry [2, 23, 71, 80] allows for analysis of up to 50 markers, but no transcriptomes [23, 61]. Single-cell RNA sequencing (scRNA-seq) allows the interrogation of expressed genes [17, 45, 53, 72, 87] and surface markers (cellular indexing of transcriptional epitope sequencing, CITE-Seq) [67, 72, 87].

In immune cells, the correlation between mRNA and surface expression of any given surface marker is, in the majority of the cases, weak [42, 67]. The correlation between mRNA and protein levels in mammalian immune cells is around 40% [69]. This is because cell surface expression is not only determined by gene expression, but also by post-translational protein modifications [43], trafficking to the cell surface, protein stability, and proteolytic modifications. The capture efficiency of mRNA is not perfect, and mRNA drop-outs further weaken the correlation between gene and surface protein [87]. Cell types in PBMCs have been defined by flow cytometry, and the surface markers of the major cell types are very well known. Yet, it is difficult to call even major cell types by scRNA-seq. For example, CD4+ T cells were not resolved from CD8+ T cells and natural killer (NK) cells [83]. To take advantage of the extensive knowledge and vast literature in flow and mass cytometry, it is necessary to assess cell surface phenotype along with transcriptomes.

Currently, only two publications report combined single-cell transcriptomes and protein (antibodies) from patients with atherosclerosis. The first one [19] includes human carotid endarterectomy specimens and matched PBMCs where 1652 PBMCs from one patient (without PBMCs from a healthy control) were analyzed by 10x Genomics 3′ and CITE-seq [53, 72], using a panel of 21 oligonucleotide-tagged antibodies. The second paper analyzed human carotid arteries from explanted hearts with 10x Genomics and a panel of 11 oligonucleotide-tagged antibodies. It is focused on smooth muscle cells and fibromyocytes rather than immune cells [81]. A recent study reported the effect of HIV infection on PBMC transcriptomes [33], focusing on acute HIV infection (before antiretroviral therapy started) and reporting PBMC transcriptomes in four patients at 8 defined time points (average of 1976 PBMC transcriptomes per participant and condition). No scRNA-seq or CITE-seq studies of PBMCs of people living with chronic HIV infection have been reported and no single-cell studies of the interaction between HIV and cardiovascular disease (CVD) are available.

Here, we report transcriptomes and cell surface phenotypes of almost 42,000 PBMCs using the targeted scRNA-seq BD Rhapsody platform [17, 45] that simultaneously provides single-cell surface phenotype (40 monoclonal antibodies, mAbs) and transcriptomes (485 immune and inflammatory transcripts) in the same cells. As a proof-of-concept, we show significant differences in cell proportions and cell transcriptomes between healthy subjects and matched subjects living with HIV or cardiovascular disease from the WIHS cardiovascular sub-study. WIHS is an ongoing multi-center, prospective, observational cohort study of women with or at risk of HIV infection. PBMCs were cryopreserved on liquid N2, following strict standard operating procedures that ensured preservation of cell surface phenotype, viability, and transcriptomes.

Methods

Study characteristics and sample selection

The Women’s Interagency HIV Study (WIHS) was initiated in 1994 at six (now expanded to ten) US locations [25, 28]. It is an ongoing prospective study of over 4000 women with or at risk of HIV infection. Recruitment in the WIHS occurred in four phases (1994–1995, 2001–2002, 2010–2012, and 2013–2015) from HIV primary care clinics, hospital-based programs, and community outreach and support groups. Briefly, the WIHS involves semi-annual follow-up visits, during which participants undergo similar detailed examinations, specimen collection, and structured interviews assessing health behaviors, medical history, and medication use. All participants provided informed consent, and each site’s Institutional Review Board approved the studies.

All participants in the current analysis were part of a vascular sub-study nested within the WIHS [25, 28, 32]. The baseline visit for the vascular sub-study occurred between 2004 and 2006, and a follow-up visit occurred on average 7 years later. Participants underwent high-resolution B-mode carotid artery ultrasound to image six locations in the right carotid artery: the near and far walls of the common carotid artery, carotid bifurcation, and internal carotid artery. A standardized protocol was used at all sites [32], and measurements of the carotid artery focal plaque, a marker of subclinical atherosclerosis, were obtained at a centralized reading center (University of Southern California). Subclinical CVD (sCVD) was defined based on the presence of one or more carotid artery lesions [32].

From the initial 1865 participants in the WIHS vascular sub-study, 32 participants were selected for scRNA-seq analysis. sCVD was defined as the presence of carotid artery focal plaque at either vascular sub-study visit to define four groups of eight participants each: HIV−, HIV+CVD−, HIV+CVD+, and HIV+CVD+ on cholesterol reduce treatment (CRT). Because we were interested in the joint relationships of HIV infection and sCVD with surface marker and RNA expression by different cell subtypes, we selected matched samples based on HIV, sCVD, and CRT (mostly statins) (Additional file 1: Fig. S1). The latter was done because we found that CRT had a major impact on monocyte transcriptomes [16]. HIV infection status was ascertained by enzyme-linked immunosorbent assay (ELISA) and confirmed by Western blot. Non-sCVD participants with self-reported coronary heart disease or current lipid-lowering therapy use were excluded. Participants were formed in quartets matched by race/ethnicity (except one quartet), age (± 5 years) at the baseline vascular sub-study (except one quartet where the age difference was more but all the women were post-menopausal), visit number, smoking history, and date of specimen collection (within 1 year).

Demographic, clinical, and laboratory variables were assessed from the same study visit using standardized protocols. The median age at the baseline study visit was 55 years, and 96% of participants were either of Black race or Hispanic ethnicity. Most (86%) reported a history of smoking. Substance use was highly prevalent, with 43% of HIV+ and 50% of HIV− participants reporting either a history of injection drug use; current use of crack, cocaine, or heroin; or alcohol use (≥14 drinks per week). Among HIV+ participants, over 80% reported use of highly active antiretroviral therapy (HAART) at the time PBMCs were obtained, and 59% reported an undetectable HIV-1 RNA level. The median CD4+ T-cell count was 585 cells/μL (IQR 382–816) in HIV+ women without sCVD and 535 cells/μL (IQR 265–792) in HIV+ women with sCVD. From now on along the text, sCVD will be referred as CVD.

Preparation of PBMC samples for CITE-seq

To avoid batch effects, sixteen samples each were processed on the same day. PBMC tubes were thawed in a 37°C water bath and tubes filled with 8 mL of complete RPMI-1640 solution (cRPMI-1640) which contains human serum albumin, HEPES, sodium pyruvate, MEM-NEAA, penicillin-streptomycin, GlutaMax, and mercaptoethanol. Main reagents, manufacturers, and catalogue numbers are listed in Additional file 1: Table S1. The tubes were centrifuged at 400 ×g for 5 min and pellets resuspended in cold staining buffer (2 % fetal bovine serum (FBS) in phosphate-buffered saline (PBS)). Manual cell counting (trypan blue solution, 0.4%) was performed by diluting cell concentration to achieve 100–400 cells per hemocytometer count. Cells were aliquoted to a count of 1 million cells each and incubated on ice with Fc Block reagent (BD Biosciences, Additional file 1: Table S1) at a 1:20 dilution, centrifuged at 400 ×g for 5 min, resuspended in 180 μL of staining buffer, and transferred to their respective sample multiplexing kit tubes (BD Biosciences). The cells were incubated for 20 min at room temperature, transferred to 5-mL polystyrene tubes, washed 3 times, and centrifuged at 400 ×g for 5 min. The cells were resuspended in 400 μL of staining buffer and 2 μL of 0.3mM DRAQ7 and 2 μL of 2mM Calcein AM were added to each tube. The viability and cell count of each tube were determined using the BD Rhapsody scanner (Additional file 1: Table S2). Tube contents were pooled in equal proportions with total cell counts not to exceed 1 million cells. The tubes were then centrifuged at 400 ×g for 5 min and resuspended in a cocktail of 40 oligonucleotide-tagged antibodies (listed in Additional file 1: Table S3) (2 μL each antibody and 20 μL of staining buffer) on ice for 30–60 min per manufacturer’s recommendations. The tubes were then washed with 2 mL of staining buffer followed by centrifugation at 400 ×g for 5 min. This was repeated two more times for a total of 3 washes. The cells were then counted again using the scanner and loaded into Rhapsody nanowell plates (4 samples per plate).

Library preparation

Cells were loaded at 800–1000 cells/μL into the primed plate per the BD user guide. The beads were isolated with a magnet and the supernatant removed. Reverse transcription was performed at 37 °C on a thermomixer at 1200 rpm for 20 min. Exonuclease I was incubated at 37 °C on a thermomixer at 1200 rpm for 30 min and then immediately placed on a heat block at 80 °C for 20 min. The tube was placed on ice followed by supernatant removal while beads were on a magnet. The beads were resuspended in BD bead resuspension solution. The tubes were stored at 4 °C until further processing.

Per BD’s protocol, the reagents for PCR1 including the BD Human Immune Response Panel and a custom panel of ~100 genes (Additional file 1: Table S4) were added to the beads. Samples were aliquoted into strip PCR tubes and incubated for 10 cycles according to BD’s protocol for PCR1. A double size selection was performed to remove high genomic DNA fragments by adding 0.7× volume AMPure XP SPRI beads to the PCR products. After incubation, the supernatant is recovered and transferred to a new tube followed by purifying the supernatant with an additional 100 μL of AMPure XP beads. The content was eluted off the beads using 30 μL of BD elution buffer and then transferred to a 1.5-mL tube.

Pre-sequencing quality control (QC)

Each tube had 12 cycles of PCR performed according to BD’s user guide. The tubes were cleaned with AMPure XP beads at 0.8X for mRNA and 1.2X for sample tags. Two 200-μL washes per sample were performed during the clean-up using 80% ethanol. The cDNA was eluted off using BD elution buffer. QC and quantification was performed using Agilent TapeStation high sensitivity D1000 screen tape and Qubit double-stranded high sensitivity DNA test kit. The mRNA was then diluted, if necessary, to a concentration of 1.2–2.7 ng/μL and the antibody and sample tag libraries from PCR2 were diluted, if needed, to a concentration of 0.5–1.1 ng/μL. From each sample, 3 μL was added to a volume of 47 μL of reagents for PCR3 as described by BD’s user recommendations following the protocol and number of cycles listed, except for AbSeq, which had 9 cycles of PCR performed as determined by previous optimization. The three libraries were then cleaned with AMPure XP beads at 0.7X for AbSeq and 0.8X for sample tags. Samples were washed twice with 200 μL of 80% ethanol. The cDNA was eluted off the beads using BD’s elution buffer. Final QC and quantification was performed using TapeStation and Qubit kits and reagents.

Sequencing

The samples were pooled and sequenced to the following nominal depth recommended by BD: AbSeq: n × 1000 reads per cell, where n is the plexity of AbSeq used; mRNA: 20,000 reads per cell; sample tags: 600 reads per cell. Thus, a total of 60,600 reads per cell were desired for sequencing on the NovaSeq. The samples and specifications for pooling and sequencing depth, along with number of cells loaded onto each plate, were optimized for S1 and S2 100 cycle kits (Illumina) with the configuration of 67 × 8 × 50 bp. Once sequencing was complete, a FASTA file was generated by BD as a reference for our AbSeq and genes we targeted with these assays. The FASTA file and FASTQ files generated by the NovaSeq were uploaded to Seven Bridges Genomics pipeline, where the data was filtered and matrices and csv files were generated. This analysis generated draft transcriptomes and surface phenotypes of 54,078 cells (496 genes, 40 antibodies). Eleven genes were not expressed, leaving 485 genes for analysis.

Doublet removal

Based on the 4 sample tags used per plate, 8359 doublets were removed. The remaining 45,719 cells were analyzed using the Doublet Finder package on R (https://github.com/chris-mcginnis-ucsf/DoubletFinder) with the default doublet formation rate (7.5%). This removed another 3322 doublets, leaving 42,397 Cells. Finally, we removed all cells that had less than 128 (27) antibody molecules sequenced. This removed 786 noisy cells, resulting in 41,611 cell transcriptomes. All antibody data were CLR (centered log-ratio) normalized and converted to log2 scale. All transcripts were normalized by total UMIs in each cell and scaled up to 1000.

Thresholding

Preliminary experiments showed that each antibody had both specific and non-specific binding, as expected. To remove the non-specific signal, a threshold value separating noise from expression for each surface marker was obtained as follows (Additional file 1: Table S5). Density plots for expression of each surface marker in the main cell types (Ridgeline plots) were used to define the signal in a known negative cell population or by deconvolution of overlapping normal distributions (we used the function “normalmixEM” to deconvolute the overlapping distributions in the R package “mixtools”). The intersect of the two Gaussian distributions was then set as minimum expression threshold for the antibody, setting surface expression values below the threshold to zero. In combined protein and transcript panel single-cell sequencing, non-specific background staining is caused by incomplete Fc block and oligonucleotide-tagged antibody being trapped in the nanowell [72] as well as by incorrect titrations. Some antibodies like most antibodies to chemokine receptors have inherent background. The adjustment of antibody concentration improves signal and lowers the background [8]. Ridgeline plots of the thresholded protein expressions for each main cell type are shown in Additional file 1: Fig. S2A, which indicates how the thresholding worked on each protein expression. Based on Fig. S2A, CCR7 (CD197) antibody data were not included in the analysis.

Clustering

Prior to clustering the data based on antibodies, we ensured that the data were batch-corrected using the Harmony package. To prepare the data for clustering, we first reduced the dimensionality of the data using UMAP (Uniform Manifold Approximation and Projection) [1] to visualize the clusters. UMAP is a manifold learning technique that helps find the latent space in which the data lies within the higher dimension space by reducing the dimensions of data. It is a dimension reduction technique used for visualization. We use Louvain clustering [3] in order to cluster the data. The parameter “resolution,” which determines the quality of clustering, is set 0.8 for B cells, 1.0 for CD4+ T cells, 1.3 for CD8+ T cells, 0.5 for classical monocytes, 0.4 for intermediate monocytes, 0.4 for nonclassical monocytes, and 0.3 for NK cells. Subclustering of each major cell type was based on all non-negative antibodies (Additional file 1: Table S6). Gates were overlaid and used in all subsequent UMAP figures (cell numbers in each cluster in Additional file 1: Table S7).

Cluster assignment

Clustering was done just using antibodies. In CD4+ T cells, 4 of the initial clusters were further divided based on the expression of CD11c, CD56, CD25, CD127, CXCR3, and CCR2. CD8+ T cells had two clusters that were divided based on CD11c, CD16, and CXCR3 protein expression. One cluster from classical monocytes and one cluster from intermediate monocytes were further divided based on CD163 and CD152 expression, respectively. In nonclassical monocytes, one cluster showed differential expression of CD36 and CD152 expression and was divided in two. In B cells, one cluster was split because it showed differential expression of CD25 and CXCR3 within the cluster. Finally, two clusters from NK cells were split due to CD16, CD56, and CD11c expression.

Comparing gene expression among participant types

To determine differential expression (DE) among the four types of participants, we ran Wilcoxon rank sum test in the Seurat package [73] in R with no thresholds over avg_logFC, minimum fraction of cells required in the two populations being compared, minimum number of cells and minimum number of cells expressing a feature in either group. We filtered for adjusted p < 0.05 and compared HIV−, HIV+CVD−, HIV+CVD+, and HIV+CVD+CRT+. From this data, dotplots were generated using ggpubr package in R. Significant genes were selected based on an adjusted p-value (Bonferroni corrected) threshold of < 0.05 and pct.1 (percentage of cells where the feature is detected in the first group) value > 0.2. Significantly differentially expressed genes (exact p-values) for each major cell types are shown in Additional file 1: Table S8.

Comparing cell proportions

To find changes in proportions, we identified the cell numbers for each participant in each cluster (Additional file 1: Table S9). Statistical differences in cell proportions were calculated by log-odds ratio defined as p/(1−p) where p is the proportion of cells, followed by ANOVA and Tukey’s multiple comparison test between the four groups. For clarity, the data are presented as percentage of cells.

Correlation analysis

We correlated each antibody to its corresponding gene(s) using Spearman rank correlation and significance (R package). For each combination of gene-antibody, we discarded cells that had values below the corresponding threshold for that antibody as well as cells with zero counts for that gene. After this filter, any gene-antibody combination that had 10 cells or less was deemed insignificant. Finally, all non-significant (p-value > 0.05) were designated a nominal value of zero as the Spearman rank correlation coefficient and we selected only those genes or antibodies that had at least one correlation whose coefficient ≥ 0.25 or whose coefficient ≤ −0.25. All significant non-negative correlations are reported in Additional file 1: Table S10.

Random forest analysis

A machine learning (ML) approach was implemented to identify the genes that distinguish between disease groups. To accomplish this goal, the Random Forest (RF) model [37, 40] was trained with the normalized gene expression from 1000 randomly selected cells from each condition and variable importance scores of the genes were calculated. This procedure was repeated for 15 iterations and importance scores in each iteration were scaled to 0–100. A higher score indicated a higher power for classifying the disease groups.

Results

Identification of main cell types based on surface marker expression

To identify the major known cell types (Fig. 1A), we used ridgeline plots on CD3, CD4, CD8, CD14, CD16, CD19, CD56, CD123, and CD206. This approach defines (Fig. 1B):

*

B cells: CD19+ CD3−

*

T cells: CD19− CD3+

*

CD4+ T cells: CD4+ CD8− T cells

*

CD8+ T cells: CD8+ CD4− T cells

*

Monocytes (M): CD19−CD3−CD56−

*

Classical (CM): CD14+CD16−

*

Intermediate (INT): CD14+CD16+

*

Nonclassical (NCM): CD14−CD16+CD56−

*

NK cells (NK): CD4− CD56+ CD14− CD20− CD123− CD206−

[IMAGE OMITTED: SEE PDF]

This resulted in 2835 B cells, 11,019 CD4+ T cells, 10,865 CD8+ T cells, 5145 CM, 995 INT, 475 NCM, and 1801 NK cells. Each of these major cell types was then re-clustered separately, using Seurat [73] to construct UMAPs, an effective tool for visualizing high-dimensional data, with Louvain clustering (Fig. 1C). Like in flow or mass cytometry, we clustered on antibody staining only. This “preserves” the transcriptomes for investigations into disease- and treatment-related changes. Using this approach, we identified 13 CD4+ T cell subsets, 14 CD8+ T cell subsets, 7 CM subsets, 3 NCM subsets, 3 INT subsets, 7 B cell subsets, and 3 NK cell subsets (Fig. 1C). The age of each of the participants was projected onto the main cell type UMAPs to study the possible effect of this variable in the results. Similar age distribution was observed in CD4+ T cells, CD8+ cells, CM, INT, and NCM (Additional file 1: Fig. S2B). In the case of B cells and NK UMAPs, some of the clusters presented unequal age distribution, but those clusters did not show any remarkable finding in this study (Fig. 1D).

Cell subset calling using 40 surface markers

Next, we constructed heatmaps for those antibodies that were significantly differentially expressed in at least one subset and relevant for cell calling (Fig. 2). Based on Fig. S2A, CCR7 (CD197) antibody data were not included in the analysis. This information allowed us to call all CD4+ T and CD8+ T cell subsets in accordance with the published work. Among CD4+ T cells, CD2 was expressed in almost all cells, as expected. The high affinity IL2 receptor IL2RA (CD25) was expressed in about a third of the CD4+ T cells and was strikingly high in cluster 10, which was also low for IL7 receptor (CD127), defining cluster 10 as regulatory T cells (Tregs). CD45RA and RO were mutually exclusive, separating naive and antigen-experienced CD4+ T cells. CXCR3 (CD183) identifies T-helper-1 (Th1) cells and was highly expressed in cluster 1. Cluster 9 expressed CXCR5 (CD185) as the only chemokine receptor, suggesting it may contain follicular helper (Tfh) T cells. Cluster 9 was the only cluster with high expression of CD56, suggesting a CD56+ CD4+ T cell. Based on protein information, all CD4 T cell clusters were called (Fig. 2A). Most CD8+ T cells expressed CD2. Chemokine receptors (CD183, 184, 192, 194, 195, 197) were expressed on cells in clusters 5, 7, and 8. Cluster 6 was identified as NK-like (CD56+) T cells with a CD45RA+ terminally differentiated memory (Emra) phenotype (Fig. 2B).

[IMAGE OMITTED: SEE PDF]

Among monocytes, we were able to call all the classical monocyte subsets based on published data [23, 75], related to subsets described by mass cytometry [23]. All CM were CD11b+ (Fig. 2C). There were gradients of CD9, CD69, and CD184 expression. The scavenger receptor CD36 and the chemokine receptor CCR2 (CD192) were expressed in all classical monocytes. Some markers were striking in each of the CM subsets. For example, CM cluster 7 expressed high levels of CD142 (tissue factor), which has previously been implicated in people living with HIV [68], whereas cluster 5 expressed high CD137 and cluster 3 high expression of CD163 (hemoglobin-haptoglobin receptor). INT CD14+CD16+ monocytes have been considered pro-inflammatory and are known to be increased in people with HIV [24] and with CVD [62, 74]. All INT monocytes highly expressed the activation molecules CD69, CD9, and CD36 (Fig. 2C). Since INT subsets have not been described before, we propose a provisional naming suggestion (Fig. 2C) based on their most remarkable expressed markers. NCM formed 3 clusters. Strikingly, the expression of CD9 and CD36 was limited to cluster 2, suggesting that this cluster corresponds to the previously described CD9+CD36+ NCM [23]. CD11c, CD86, CD141, and CD152 were expressed in all NCMs (Fig. 2C).

Using 18 surface antibodies to subtype B cells obtained from women with coronary artery diseases, B cells were subtyped into 7 clusters (Fig. 2D) with 3 clusters of CD27+ B cells and 4 clusters of CD27− B cells. CD27 is known as a conventional memory B cell marker [31]. Therefore, clusters 4, 6, and 7 were classified as memory B cell populations. Clusters 6 and 7 were notable for high CD11c expression. CD11c+ memory B cells have been shown to be precursors of antibody secreting cells [22]. Cluster 4 was remarkable for CD25+, and CD25+ B cells are known as antigen presenting B cells [5]. Clusters 1, 2, 3, and 5 were CD27− B cells, likely a combination of naïve B cells and transitional B cells.

Among NK cells, we found 3 clusters. Cluster 1 contained HLADR+ mature NK cells with CD56-CD16high expression, an NK cell subset known to be elevated in chronic HIV infection [29]. Cluster 2 of NK cells was mature (CD56dim/CD16+). Cluster 3 contained immature (CD56brightCD16−) NK cells (Fig. 2E).

Changes in PBMC subset abundance by disease or treatment

Based on this data, it is possible to address shifts in cell proportion based on disease or treatment: HIV effect, comparing HIV− vs HIV+; the cardiovascular disease effect in women living with HIV, comparing HIV+ vs HIV+CVD+; and the effect of cholesterol reduction treatment, comparing HIV+ CVD+CRT− vs HIV+ CVD+ CRT+. We found significant differences in cell proportions (p) (calculated by log odds ratio, p/(1−p) followed by ANOVA and Tukey’s multiple comparison test) in one intermediate monocyte subsets, one nonclassical monocyte subset, one CD8+ T cell, and one B cell subset (Fig. 3). Among B cells, antigen presenting memory B cells (Fig. 3A) were severely lower in all WIHS participants with HIV with or without subclinical CVD. Within CD8+ T cells, cluster 7 (CXCR3+ CD8+ T cells) was reduced in patients with HIV and CVD with CRT treatment compared to those without cholesterol-reducing treatment (Fig. 3B). Strikingly, two subsets of monocytes showed significantly different abundances. Cluster 2 from intermediate monocytes (CD11bhi CD9lo) (Fig. 3C) was elevated in patients with HIV and cluster 3 from nonclassical monocytes (CD9hi CXCR4+) (Fig. 3D) was significantly elevated in WIHS participants living with HIV with subclinical CVD.

[IMAGE OMITTED: SEE PDF]

Differential gene expression in each of the clusters

Since the transcriptomic information was not used for UMAPs and clustering shown in Fig. 1C, we were able to run unbiased gene expression patterns in each cell subset against all other clusters within the same cell type. We filtered for genes that were significantly differentially expressed in at least one of the subsets (Fig. 4, Additional file 1: Table S8). This analysis revealed gene signatures for most subsets and confirmed the identity of the cell clusters identified by CITE-seq (Fig. 2) and expanded phenotype information. Cluster 4 CD36+ effector memory CD4+ T cells overexpressed GNLY, GZMA, GZMH, NKG7, and FGFBP2, which are known as effector memory (Em) signature genes. Cluster 10 expressed regulatory T cell (Treg) signature genes including FoxP3, TIGIT, CTLA4, and LGALS3. Cluster 9 Tfh Em CD4+ T cells expressed CXCR5, the classical chemokine receptor characteristic for Tfh cells [13]. CD56+ CD4+ T cells (cluster 11) expressed NK signatures such as HOPX, CTSW, KLRC4, and KLRK (Fig. 4A). Within CD8+ T cells, cluster 2 (naïve) overexpressed CCR7, SELL, or LEF1, representative genes of naïve cells. NK-like phenotype genes such as GNLY, KLRF1, GZMB, and FCGR3A were expressed in cluster 6 (CD56+ EMRA) (Fig. 4B). CTLA4+ CM (cluster 7) expressed CCL20, SOD2, or AQP9, classical monocyte markers (Fig. 4C). No significant differences were found for B cells and NK cells.

[IMAGE OMITTED: SEE PDF]

Transcriptomes shift with HIV, CVD, and cholesterol control

We hypothesized blood immune cell transcriptomes change with disease state. To test this possibility, we constructed dotplots for each of the disease status and compared HIV− with HIV+, HIV+ with HIV+ with CVD, and HIV+CVD+ with HIV+ CVD+ with CRT treatment (Fig. 5, Additional file 1: Fig. S3, Additional file 2: Data S1). Many genes in CD4+ and CD8+ T cell subsets showed significant differences. Some genes in monocyte and NK cell subsets showed significant differences. In CD4+ T, CD8+ T, and NK cells (Fig. 5A, B, F), IL-32 was highly significantly increased by CVD, but not in CVD+ women on CRT. IL-32 is an inflammatory cytokine that is known to be important in CVD [14, 34]. The transcription factor JUNB and the lymphocyte-specific protein tyrosine kinase LCK were upregulated in both, CD4+ T and CD8+ T cells from women with CVD, and reduced in the presence of CRT. JUMB promotes the development of inflammatory Th17 cells and restricts flexibility towards alternative effector and regulatory programs [9]. LCK is a key molecule in the activation of the TCR signaling and T cells. Several killer cell lectin receptors (KLRC4 and KLRK1) were also significantly upregulated in CVD, but downregulated in the presence of CRT in CD8+ T cells. RUNX is reduced in CVD but increased in women on CRT [35].

[IMAGE OMITTED: SEE PDF]

DUSP1 was highly overexpressed in classical monocytes of women with CVD (Fig. 5C). DUSP1 oxidation prolongs MAPK activation, ultimately resulting in enhanced inflammatory responses [30]. In addition to CCL3, CCL4, and DUSP2, IL1B, known to be highly relevant in atherosclerosis, was highly upregulated in CM of HIV+CVD+ women. TNFSF10 (TRAIL), TNFSF13 (APRIL), and TNFSF13B (BAFF), important B cell regulators, were upregulated in CM from women with CVD. The Toll-like receptor TLR2, which is known to be involved in atherosclerosis [15, 49, 55], was upregulated by CVD. In intermediate monocytes (Fig. 5D), CCL4, TNF, IL1B, FCGR3A, and PIK3AP1 were associated with CVD in the women that did not receive CRT whereas FCN1 was increased in NCM from women with CRT (Fig. 5E). In nonclassical monocytes, CD52 was downregulated in women with HIV, consistent with a previous report [77]. CD74 in B cells (Fig. 5G) was slightly downregulated in CVD women on CRT compared to CVD. DUSP1 in B cells was downregulated in HIV+CVD+ women compared to HIV+ CVD−. The clusters with genes that passed the filtering process (p-value <0.05, avg.Log2FC>0 or <0, and pct.1 > 0.2) from each of the main cell types are presented in Fig. S3.

Random Forest analysis to identify genes that distinguish between disease groups

We used the Random Forest Machine learning approach to identify the genes with the highest capability to distinguish between disease groups (Additional file 2: Data S2). We showed the top 50 ranked genes to separate HIV vs healthy (Fig. 6A), HIV vs HIV+CVD+ (Fig. 6B), and HIV+CVD+ vs HIV+CVD+CRT+ (Fig. 6C), with the importance ranked between 0 and 100. The histograms in Fig. 6A–C show overlaid ridgeline plots for key genes. KLF2 was the most important gene to separate HIV from non-HIV participants. Both in CVD and in CVD+CRT+ in HIV, IL32, and CD52 were highly ranked (Fig. 6B, C). Interestingly, these genes were upregulated in CVD and conversely were downregulated in CVD+CRT+. Many of the genes that were important in distinguishing between women with and without CVD, including DUSP1, DUSP2, CCL5, and LGALS1, were regulated in the opposite direction by CRT.

[IMAGE OMITTED: SEE PDF]

Discussion

In immunology, surface markers are widely used to define and distinguish cell types [11, 71, 80]. Flow cytometry is the discipline-defining method of immunology [60]. Similar to flow cytometry, in CyTOF, single-cell suspensions are stained with antibody panels to detect cellular antigens. Unlike CyTOF, scRNA-seq allows the detection of single-cell transcriptomes. Since the correlation between cell surface protein and mRNA expression is weak in immune cells [42], the transcriptome provides a valuable additional dimension. scRNA-seq without surface phenotype information has led to much frustration in the field, because the expression of many genes encoding well-known surface markers remains undetected in scRNA-seq [10, 44, 80]. It is still difficult to call cell types based on gene expression data alone, which emphasizes the need for cell surface phenotypes in addition to transcriptomes. Here, we correlated gene expression with cell surface expression for 41 pairs of genes and proteins. CD74 protein expression was well correlated with the expression of both the CD74 and the HLA-DRA genes. CD4 and CD16 surface and gene expression were reasonably well correlated across all cell types. A few other genes including CD14, CD16, IL-3 receptor (CD123), and CD27 were somewhat correlated with the surface expression of their proteins in some cell types. For most markers, we confirm weak correlations (Additional file 1: Table S10) [42, 67, 69], which illustrates the value of monitoring cell surface phenotype in scRNA-seq.

PBMCs can be analyzed without mechanical or enzymatic dissociation, which are known to alter cell surface markers and transcriptomes [79. ]. PBMCs are attractive for single-cell RNA sequencing (scRNA-seq) studies, because they are available in many clinical studies of specific populations with defined diseases and outcomes. The participants sampled for the present study were part of a sub-study nested within the WIHS [25, 28, 32], which provided detailed information on subclinical atherosclerosis. Participants underwent high-resolution B-mode carotid artery ultrasound to image six locations in the right carotid artery [28]. Although our study is not definitive, it is suggestive of significant changes in cell proportions and transcriptomes in subjects with cardiovascular disease.

scRNA-seq has been applied to human PBMCs in diseases including cancers [6, 84,85,86], inflammatory bowel disease [48, 76], and autoimmune disease [36, 57], as well as atherosclerosis [18, 19, 81]. The foundational paper for the 10x Genomics drop-Seq method [87] demonstrated the feasibility of using scRNA-seq on PBMCs. Other studies reported scPBMC transcriptomes in colorectal cancer [85], γδ T cells [54], liver cancer [86], in vitro Salmonella infection [4], and memory T cells [45]. Two publications have reported combined single-cell transcriptomes and proteomics from patients with atherosclerosis. Fernandez et al. [19] ran antibody sequencing in only one human plaque, but its power was revealed by identifying five distinct macrophage clusters. Wirka et al. [81] revealed that the process of smooth muscle cell phenotypic modulation in vivo can be altered by the expression of Tcf21, a gene causally associated with reduced risk of coronary artery disease. The loss of Tcf21 results in fewer fibromyocytes in the lesions and the protective fibrous cap [81].

In the current study, four clusters showed significantly different abundance of cells in the four groups of participants. One of them is an intermediate monocyte subset, which underscores the extraordinary importance of this cell type in chronic HIV infection [26, 47]. Intermediate monocyte numbers have previously been found increased in non-HIV subjects with peripheral artery occlusive disease [78] and significantly predicted cardiovascular events [27, 62, 63]. A nonclassical monocyte cluster was increased in CVD, which was reported previously [20, 39].

CD8+ T cells are abundant in atherosclerotic plaques in humans [21] and they are found in higher numbers than CD4+ cells [19, 21, 65]. In advanced human lesions, CD8+ T cells are mostly found in fibrous cap areas [52]. Cluster 7 of CD8 + T cells (CXCR3+) was reduced in women on CRT. CD8+ T cells express higher levels of CXCR3 in patients with symptomatic atherosclerosis [19]. We found IL-32 highly expressed in most T and NK cell clusters in CVD. IL-32 is a 27-kDa cytokine expressed in T cells and monocytes that is secreted after apoptosis [50]. It is an inflammatory cytokine that drives IL-1β, clinically important in CVD [59], TNF, IL-6, and IL-8 expression [14, 34, 50]. IL-32 activates the leukocyte surface protease PR3, which in turn triggers the G-protein-coupled receptor PAR2 [50] and is known to be important in viral infections [38, 50, 56, 70]. IL32 was upregulated in CVD and was downregulated in women with CVD and cholesterol reduce treatment. Since IL-32 appears to be CVD-specific, we advocate for future prospective studies in larger cohorts to determine whether IL32 mRNA is a useful biomarker.

KLF2 was the most important gene to separate HIV from non-HIV participants. Some previous studies showed that KLF2 is related to HIV infection [12, 58]. This gene and some others including CD74 and CD52 were also important to separate CVD from non-CVD in HIV, and CVD+CRT+ from CVD+, suggesting that these genes may be related to HIV and CVD. In fact, it has been reported that anti-CD52 antibodies might be effective in HIV individuals on antiretroviral therapy [77] as well as a potential diagnostic value evaluating antiretroviral efficiency [82]. The interaction between HIV and CVD observed here is novel.

Our discovery study will encourage prospective epidemiological studies to address which PBMC subsets and transcriptomes are best suited as clinical biomarkers to stratify risk and guide treatment in subjects with coronary or peripheral artery disease. The current findings have some limitations. They need to be extended to men (the current data is based on women) and other races and ethnicities (the current data is based on mostly African American and Hispanic women). Studies of CVD in non-smokers are also needed (the current data is based on smokers), and the age range needs to be broadened.

Conclusions

In conclusion, we demonstrate the utility of scRNA-seq with cell surface phenotype assessment in the same cells. The identification of 50 distinct clusters of CD4+ T and CD8+ T cells, B cells, NK cells, and monocytes helps to gain a deeper understanding of PBMCs, a rich and readily accessible source of biological and clinical information. The discovery of subsets of intermediate monocytes calls for identifying such subsets in model organisms to test their function in vivo.

Availability of data and materials

All data are available in the main text or the supplementary materials. The datasets presented in this study can be found in online repositories. The name of the repository and accession number for the data reported in this paper is Gene Expression Omnibus (GEO), GSE205320 [88].

Abbreviations

CBC:

Blood cell counter

CITE-seq:

Cellular indexing of transcriptional epitope sequencing

CLR:

Centered log-ratio

CM:

Classical monocytes

CRT:

Cholesterol reduce treatment

CVD:

Cardiovascular disease

DE:

Differential expression

ELISA:

Enzyme-linked immunosorbent assay

Em:

Effector memory

Emra:

Terminally differentiated memory

FBS:

Fetal bovine serum

HAART:

Highly active antiretroviral therapy

INT:

Intermediate

M:

Monocytes

mAbs:

Monoclonal antibodies

ML:

Machine learning

NCM:

Nonclassical monocyte

NK:

Natural killer

PBMCs:

Peripheral blood mononuclear cells

PBS:

Phosphate-buffered saline

QC:

Quality control

RF:

Random forest

scRNA-seq:

Single-cell RNA sequencing

sCVD:

Subclinical CVD

Tfh:

Follicular helper T cells

Tregs:

Regulatory T cells

UMAP:

Uniform manifold approximation and projection

WIHS:

Women’s Interagency HIV Study

References

Amezquita RA, Lun ATL, Becht E, Carey VJ, Carpp LN, Geistlinger L, et al. Orchestrating single-cell analysis with bioconductor. Nat Methods. 2020;17:137–45. https://doi.org/10.1038/s41592-019-0654-x.

Bekele Y, Lakshmikanth T, Chen Y, Mikes J, Nasi A, Petkov S, et al. Mass cytometry identifies distinct CD4+ T cell clusters distinguishing HIV-1-infected patients according to antiretroviral therapy initiation. JCI Insight. 2019;4:e125442.

Blondel VD, Guillaume J-L, Lambiotte R, Lefebvre E. Fast unfolding of communities in large networks. J Stat Mech Theory Exp. 2008;P10008. https://doi.org/10.1088/1742-5468/2008/10/p10008.

Bossel Ben-Moshe N, Hen-Avivi S, Levitin N, Yehezkel D, Oosting M, Joosten LAB, et al. Predicting bacterial infection outcomes using single cell RNA-sequencing analysis of human immune cells. Nat Commun. 2019;10:3266. https://doi.org/10.1038/s41467-019-11257-y Gene expression omnibus (GEO) https://identifiers.org/geo:GSE122084 (2019).

Brisslert M, Bokarewa M, Larsson P, Wing K, Collins LV, Tarkowski A. Phenotypic and functional characterization of human CD25+ B cells. Immunology. 2006;117:548–57. https://doi.org/10.1111/j.1365-2567.2006.02331.x.

Brown CC, Gudjonson H, Pritykin Y, Deep D, Lavallée V-P, Mendoza A, et al. Transcriptional basis of mouse and human dendritic cell heterogeneity. Cell. 2019;179:846–863.e24. https://doi.org/10.1016/j.cell.2019.09.035 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE137710 (2019) and https://identifiers.org/geo:GSE130201 (2019).

Buttarello M, Plebani M. Automated blood cell counts: state of the art. Am J Clin Pathol. 2008;130:104–16. https://doi.org/10.1309/EK3C7CTDKNVPXVTN.

Buus TB, Herrera A, Ivanova E, Mimitou E, Cheng A, Herati RS, et al. Improving oligo-conjugated antibody signal in multimodal single-cell analysis. Elife. 2021;10:e61973. https://doi.org/10.7554/eLife.61973 Figshare 10.6084/m9.figshare.c.5018987.v1, 2020.

Carr TM, Wheaton JD, Houtz GM, Ciofani M. JunB promotes Th17 cell identity and restrains alternative CD4(+) T-cell programs during inflammation. Nat Commun. 2017;8:301. https://doi.org/10.1038/s41467-017-00380-3 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE98414 (2017).

Cochain C, Vafadarnejad E, Arampatzi P, Pelisek J, Winkels H, Ley K, et al. Single-cell RNA-Seq reveals the transcriptional landscape and heterogeneity of aortic macrophages in murine atherosclerosis. Circ Res. 2018;122:1661–74. https://doi.org/10.1161/CIRCRESAHA.117.312509 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE97310, 2018.

Cole JE, Park I, Ahern DJ, Kassiteridi C, Danso Abeam D, Goddard ME, et al. Immune cell census in murine atherosclerosis: cytometry by time of flight illuminates vascular myeloid cell diversity. Cardiovasc Res. 2018;114:1360–71. https://doi.org/10.1093/cvr/cvy109.

Collins DR, Urbach JM, Racenet ZJ, Arshad U, Power KA, Newman RM, et al. Functional impairment of HIV-specific CD8(+) T cells precedes aborted spontaneous control of viremia. Immunity. 2021;54:2372–2384.e7. https://doi.org/10.1016/j.immuni.2021.08.007 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE168296 (2021).

Crotty S. T follicular helper cell biology: a decade of discovery and diseases. Immunity. 2019;50:1132–48. https://doi.org/10.1016/j.immuni.2019.04.011.

Damen MSMA, Popa CD, Netea MG, Dinarello CA, Joosten LAB. Interleukin-32 in chronic inflammatory conditions is associated with a higher risk of cardiovascular diseases. Atherosclerosis. 2017;264:83–91. https://doi.org/10.1016/j.atherosclerosis.2017.07.005.

Edfeldt K, Swedenborg J, Hansson GK, Yan Z. Expression of toll-like receptors in human atherosclerotic lesions: a possible pathway for plaque activation. Circulation. 2002;105:1158–61. https://doi.org/10.1161/01.CIR.0000012489.17433.31.

Ehinger E, Ghosheh Y, Pramod AB, Lin J, Hanna DB, Mueller K, et al. Classical monocyte transcriptomes reveal significant anti-inflammatory statin effect in women with chronic HIV. Cardiovasc Res. 2021;117:1166–77. https://doi.org/10.1093/cvr/cvaa188.

Fan HC, Fu GK, Fodor SPA. Expression profiling. Combinatorial labeling of single cells for gene expression cytometry. Science. 2015;347. https://doi.org/10.1126/science.1258367.

Fernandez DM, Giannarelli C. Immune cell profiling in atherosclerosis: role in research and precision medicine. Nat Rev Cardiol. 2022;19:43–58. https://doi.org/10.1038/s41569-021-00589-2.

Fernandez DM, Rahman AH, Fernandez NF, Chudnovskiy A, Amir E-AD, Amadori L, et al. Single-cell immune landscape of human atherosclerotic plaques. Nat Med. 2019;25:1576–88. https://doi.org/10.1038/s41591-019-0590-4.

Funderburg NT, Zidar DA, Shive C, Lioi A, Mudd J, Musselwhite LW, et al. Shared monocyte subset phenotypes in HIV-1 infection and in uninfected subjects with acute coronary syndrome. Blood. 2012;120:4599–608. https://doi.org/10.1182/blood-2012-05-433946.

Gewaltig J, Kummer M, Koella C, Cathomas G, Biedermann BC. Requirements for CD8 T-cell migration into the human arterial wall. Hum Pathol. 2008;39:1756–62. https://doi.org/10.1016/j.humpath.2008.04.018.

Golinski M-L, Demeules M, Derambure C, Riou G, Maho-Vaillant M, Boyer O, et al. CD11c(+) B cells are mainly memory cells, precursors of antibody secreting cells in healthy donors. Front Immunol. 2020;11:32. https://doi.org/10.3389/fimmu.2020.00032 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE112515 (2020).

Hamers AAJ, Dinh HQ, Thomas GD, Marcovecchio P, Blatchley A, Nakao CS, et al. Human monocyte heterogeneity as revealed by high-dimensional mass cytometry. Arterioscler Thromb Vasc Biol. 2019;39:25–36. https://doi.org/10.1161/ATVBAHA.118.311022.

Han J, Wang B, Han N, Zhao Y, Song C, Feng X, et al. CD14(high)CD16(+) rather than CD14(low)CD16(+) monocytes correlate with disease progression in chronic HIV-infected patients. J Acquir Immune Defic Syndr. 2009;52:553–9. https://doi.org/10.1097/qai.0b013e3181c1d4fe.

Hanna DB, Post WS, Deal JA, Hodis HN, Jacobson LP, Mack WJ, et al. HIV infection is associated with progression of subclinical carotid atherosclerosis. Clin Infect Dis. 2015;61:640–50. https://doi.org/10.1093/cid/civ325.

Hearps AC, Maisa A, Cheng W-J, Angelovich TA, Lichtfuss GF, Palmer CS, et al. HIV infection induces age-related changes to monocytes and innate immune activation in young men that persist despite combination antiretroviral therapy. AIDS. 2012;26:843–53. https://doi.org/10.1097/QAD.0b013e328351f756.

Heine GH, Ulrich C, Seibert E, Seiler S, Marell J, Reichart B, et al. CD14(++)CD16+ monocytes but not total monocyte numbers predict cardiovascular events in dialysis patients. Kidney Int. 2008;73:622–9. https://doi.org/10.1038/sj.ki.5002744.

Hodis HN, Mack WJ, Lobo RA, Shoupe D, Sevanian A, Mahrer PR, et al. Estrogen in the prevention of atherosclerosis. A randomized, double-blind, placebo-controlled trial. Ann Intern Med. 2001;135:939–53. https://doi.org/10.7326/0003-4819-135-11-200112040-00005.

Hong HS, Eberhard JM, Keudel P, Bollmann BA, Ahmad F, Ballmaier M, et al. Phenotypically and functionally distinct subsets contribute to the expansion of CD56-/CD16+ natural killer cells in HIV infection. AIDS. 2010;24:1823–34. https://doi.org/10.1097/QAD.0b013e32833b556f.

Hoppstädter J, Ammit AJ. Role of dual-specificity phosphatase 1 in glucocorticoid-driven anti-inflammatory responses. Front Immunol. 2019;10:1446. https://doi.org/10.3389/fimmu.2019.01446.

Kaminski DA, Wei C, Qian Y, Rosenberg AF, Sanz I. Advances in human B cell phenotypic profiling. Front Immunol. 2012;3:302. https://doi.org/10.3389/fimmu.2012.00302.

Kaplan RC, Kingsley LA, Gange SJ, Benning L, Jacobson LP, Lazar J, et al. Low CD4+ T-cell count as a major atherosclerosis risk factor in HIV-infected women and men. AIDS. 2008;22:1615–24. https://doi.org/10.1097/QAD.0b013e328300581d.

Kazer SW, Aicher TP, Muema DM, Carroll SL, Ordovas-Montanes J, Miao VN, et al. Integrated single-cell analysis of multicellular immune dynamics during hyperacute HIV-1 infection. Nat Med. 2020;26:511–8. https://doi.org/10.1038/s41591-020-0799-2.

Kim S-H, Han S-Y, Azam T, Yoon D-Y, Dinarello CA. Interleukin-32: a cytokine and inducer of TNFalpha. Immunity. 2005;22:131–42. https://doi.org/10.1016/j.immuni.2004.12.003.

Korinfskaya S, Parameswaran S, Weirauch MT, Barski A. Runx transcription factors in T cells-what is beyond thymic development? Front Immunol. 2021;12:701924. https://doi.org/10.3389/fimmu.2021.701924.

Kotliarov Y, Sparks R, Martins AJ, Mule MP, Lu Y, Goswami M, et al. Broad immune activation underlies shared set point signatures for vaccine responsiveness in healthy individuals and disease activity in patients with lupus. Nat Med. 2020:618–29. https://doi.org/10.1038/s41591-020-0769-8.

Kuhn M. Building predictive models in R using the caret package. J Stat Softw. 2008;28:1–26. https://doi.org/10.18637/jss.v028.i05.

Li W, Sun W, Liu L, Yang F, Li Y, Chen Y, et al. IL-32: a host proinflammatory factor against influenza viral replication is upregulated by aberrant epigenetic modifications during influenza a virus infection. J Immunol. 2010;185:5056–65. https://doi.org/10.4049/jimmunol.0902667.

Liang H, Xie Z, Shen T. Monocyte activation and cardiovascular disease in HIV infection. Cell Mol Immunol. 2017;14:960–2. https://doi.org/10.1038/cmi.2017.109.

Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2:18–22.

Liechti T, Roederer M. OMIP-060: 30-parameter flow cytometry panel to assess T cell effector functions and regulatory T cells. Cytom A. 2019;95:1129–34. https://doi.org/10.1002/cyto.a.23853.

Liu Y, Beyer A, Aebersold R. On the dependency of cellular protein levels on mRNA abundance. Cell. 2016;165:535–50. https://doi.org/10.1016/j.cell.2016.03.014.

Londino JD, Gulick DL, Lear TB, Suber TL, Weathington NM, Masa LS, et al. Post-translational modification of the interferon-gamma receptor alters its stability and signaling. Biochem J. 2017;474:3543–57. https://doi.org/10.1042/BCJ20170548.

Lundberg E, Fagerberg L, Klevebring D, Matic I, Geiger T, Cox J, et al. Defining the transcriptome and proteome in three functionally different human cell lines. Mol Syst Biol. 2010;6:450. https://doi.org/10.1038/msb.2010.106.

Mair F, Erickson JR, Voillet V, Simoni Y, Bi T, Tyznik AJ, et al. A targeted multi-omic analysis approach measures protein expression and low-abundance transcripts on the single-cell level. Cell Rep. 2020;31:107499. https://doi.org/10.1016/j.celrep.2020.03.063.

Mair F, Prlic M. OMIP-044: 28-color immunophenotyping of the human dendritic cell compartment. Cytom A. 2018;93:402–5. https://doi.org/10.1002/cyto.a.23331 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE135325 (2020).

Martin GE, Gouillou M, Hearps AC, Angelovich TA, Cheng AC, Lynch F, et al. Age-associated changes in monocyte and innate immune activation markers occur more rapidly in HIV infected women. PLoS One. 2013;8:e55279. https://doi.org/10.1371/journal.pone.0055279.

Martin JC, Chang C, Boschetti G, Ungaro R, Giri M, Grout JA, et al. Single-cell analysis of Crohn’s disease lesions identifies a pathogenic cellular module associated with resistance to anti-TNF therapy. Cell. 2019;178:1493–1508.e20. https://doi.org/10.1016/j.cell.2019.08.008 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE134809 (2019).

Mullick AE, Soldau K, Kiosses WB, Bell TA 3rd, Tobias PS, Curtiss LK. Increased endothelial expression of toll-like receptor 2 at sites of disturbed blood flow exacerbates early atherogenic events. J Exp Med. 2008;205:373–83. https://doi.org/10.1084/jem.20071096.

Nakayama M, Niki Y, Kawasaki T, Takeda Y, Ikegami H, Toyama Y, et al. IL-32-PAR2 axis is an innate immunity sensor providing alternative signaling for LPS-TRIF axis. Sci Rep. 2013;3:2960. https://doi.org/10.1038/srep02960.

Nettey L, Giles AJ, Chattopadhyay PK. OMIP-050: a 28-color/30-parameter fluorescence flow cytometry panel to enumerate and characterize cells expressing a wide array of immune checkpoint molecules. Cytom A. 2018;93:1094–6. https://doi.org/10.1002/cyto.a.23608.

Paul VSV, Paul CMP, Kuruvilla S. Quantification of various inflammatory cells in advanced atherosclerotic plaques. J Clin Diagn Res. 2016;10:EC35–8. https://doi.org/10.7860/JCDR/2016/19354.7879.

Peterson VM, Zhang KX, Kumar N, Wong J, Li L, Wilson DC, et al. Multiplexed quantification of proteins and transcripts in single cells. Nat Biotechnol. 2017;35:936–9. https://doi.org/10.1038/nbt.3973 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE100501 (2017).

Pizzolato G, Kaminski H, Tosolini M, Franchini D-M, Pont F, Martins F, et al. Single-cell RNA sequencing unveils the shared and the distinct cytotoxic hallmarks of human TCRVdelta1 and TCRVdelta2 gammadelta T lymphocytes. Proc Natl Acad Sci U S A. 2019;116:11906–15. https://doi.org/10.1073/pnas.1818488116 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE128223 (2019).

Quillard T, Araújo HA, Franck G, Shvartz E, Sukhova G, Libby P. TLR2 and neutrophils potentiate endothelial stress, apoptosis and detachment: implications for superficial erosion. Eur Heart J. 2015;36:1394–404. https://doi.org/10.1093/eurheartj/ehv044.

Rasool ST, Tang H, Wu J, Li W, Mukhtar MM, Zhang J, et al. Increased level of IL-32 during human immunodeficiency virus infection suppresses HIV replication. Immunol Lett. 2008;117:161–7. https://doi.org/10.1016/j.imlet.2008.01.007.

Reyes M, Vickers D, Billman K, Eisenhaure T, Hoover P, Browne EP, et al. Multiplexed enrichment and genomic profiling of peripheral blood cells reveal subset-specific immune signatures. Sci Adv. 2019;5:eaau9223. https://doi.org/10.1126/sciadv.aau9223 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE120442 (2019).

Richardson MW, Jadlowsky J, Didigu CA, Doms RW, Riley JL. Kruppel-like factor 2 modulates CCR5 expression and susceptibility to HIV-1 infection. J Immunol. 2012;189:3815–21. https://doi.org/10.4049/jimmunol.1201431.

Ridker PM, Everett BM, Thuren T, MacFadyen JG, Chang WH, Ballantyne C, et al. Antiinflammatory therapy with canakinumab for atherosclerotic disease. N Engl J Med. 2017;377:1119–31. https://doi.org/10.1056/NEJMoa1707914.

Robinson JP, Roederer M. History of science. Flow cytometry strikes gold. Science. 2015;350:739–40. https://doi.org/10.1126/science.aad6770.

Rodriguez L, Pekkarinen PT, Lakshmikanth T, Tan Z, Consiglio CR, Pou C, et al. Systems-level immunomonitoring from acute to recovery phase of severe COVID-19. Cell Rep Med. 2020;1:100078. https://doi.org/10.1016/j.xcrm.2020.100078.

Rogacev KS, Cremers B, Zawada AM, Seiler S, Binder N, Ege P, et al. CD14++CD16+ monocytes independently predict cardiovascular events: a cohort study of 951 patients referred for elective coronary angiography. J Am Coll Cardiol. 2012;60:1512–20. https://doi.org/10.1016/j.jacc.2012.07.019.

Rogacev KS, Seiler S, Zawada AM, Reichart B, Herath E, Roth D, et al. CD14++CD16+ monocytes and cardiovascular outcome in patients with chronic kidney disease. Eur Heart J. 2011;32:84–92. https://doi.org/10.1093/eurheartj/ehq371.

De Rosa SC, Brenchley JM, Roederer M. Beyond six colors: a new era in flow cytometry. Nat Med. 2003;9:112–7. https://doi.org/10.1038/nm0103-112.

Rossmann A, Henderson B, Heidecker B, Seiler R, Fraedrich G, Singh M, et al. T-cells from advanced atherosclerotic lesions recognize hHSP60 and have a restricted T-cell receptor repertoire. Exp Gerontol. 2008;43:229–37. https://doi.org/10.1016/j.exger.2007.11.009.

Sahir F, Mateo JM, Steinhoff M, Siveen KS. Development of a 43 color panel for the characterization of conventional and unconventional T-cell subsets, B cells, NK cells, monocytes, dendritic cells, and innate lymphoid cells using spectral flow cytometry. Cytom A. 2020. https://doi.org/10.1002/cyto.a.24288.

Saigusa R, Durant CP, Suryawanshi V, Ley K. Single-cell antibody sequencing in atherosclerosis research. Methods Mol Biol. 2022;2419:765–78. https://doi.org/10.1007/978-1-0716-1924-7_46.

Schechter ME, Andrade BB, He T, Richter GH, Tosh KW, Policicchio BB, et al. Inflammatory monocytes expressing tissue factor drive SIV and HIV coagulopathy. Sci Transl Med. 2017;9. https://doi.org/10.1126/scitranslmed.aam5441.

Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–42. https://doi.org/10.1038/nature10098.

Smith AJ, Toledo CM, Wietgrefe SW, Duan L, Schacker TW, Reilly CS, et al. The immunosuppressive role of IL-32 in lymphatic tissue during HIV-1 infection. J Immunol. 2011;186:6576–84. https://doi.org/10.4049/jimmunol.1100277 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE16363, 2009.

Spitzer MH, Nolan GP. Mass cytometry: single cells, many features. Cell. 2016;165:780–91. https://doi.org/10.1016/j.cell.2016.04.019.

Stoeckius M, Hafemeister C, Stephenson W, Houck-Loomis B, Chattopadhyay PK, Swerdlow H, et al. Simultaneous epitope and transcriptome measurement in single cells. Nat Methods. 2017;14:865–8. https://doi.org/10.1038/nmeth.4380 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE100866, 2017.

Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM 3rd, et al. Comprehensive integration of single-cell data. Cell. 2019;177:1888–1902.e21. https://doi.org/10.1016/j.cell.2019.05.031 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE128639 (2019).

Tapp LD, Shantsila E, Wrigley BJ, Pamukcu B, Lip GYH. The CD14++CD16+ monocyte subset and monocyte-platelet interactions in patients with ST-elevation myocardial infarction. J Thromb Haemost. 2012;10:1231–41. https://doi.org/10.1111/j.1538-7836.2011.04603.x.

Thomas GD, Hamers AAJ, Nakao C, Marcovecchio P, Taylor AM, McSkimming C, et al. Human blood monocyte subsets: a new gating strategy defined using cell surface markers identified by mass cytometry. Arterioscler Thromb Vasc Biol. 2017;37:1548–58. https://doi.org/10.1161/ATVBAHA.117.309145.

Uniken Venema WT, Voskuil MD, Vila AV, van der Vries G, Jansen BH, Jabri B, et al. Single-cell RNA sequencing of blood and ileal T cells from patients with Crohn’s disease reveals tissue-specific characteristics and drug targets. Gastroenterology. 2019;156:812–815.e22. https://doi.org/10.1053/j.gastro.2018.10.046.

Vojdeman FJ, Gaardbo JC, Hartling HJ, Gelpi M, Hove-Skovsgaard M, Pedersen AE, et al. Brief report: CD52 expression on CD4+ T cells in HIV-positive individuals on cART. J Acquir Immune Defic Syndr. 2018;77:217–20. https://doi.org/10.1097/QAI.0000000000001568.

Wildgruber M, Aschenbrenner T, Wendorff H, Czubba M, Glinzer A, Haller B, et al. The “intermediate” CD14(++)CD16(+) monocyte subset increases in severe peripheral artery disease in humans. Sci Rep. 2016;6:39483. https://doi.org/10.1038/srep39483.

Williams JW, Winkels H, Durant CP, Zaitsev K, Ghosheh Y, Ley K. Single cell RNA sequencing in atherosclerosis research. Circ Res. 2020;126:1112–26. https://doi.org/10.1161/CIRCRESAHA.119.315940.

Winkels H, Ehinger E, Vassallo M, Buscher K, Dinh HQ, Kobiyama K, et al. Atlas of the immune cell repertoire in mouse atherosclerosis defined by single-cell RNA-sequencing and mass cytometry. Circ Res. 2018;122:1675–88. https://doi.org/10.1161/CIRCRESAHA.117.312513.

Wirka RC, Wagh D, Paik DT, Pjanic M, Nguyen T, Miller CL, et al. Atheroprotective roles of smooth muscle cell phenotypic modulation and the TCF21 disease gene as revealed by single-cell analysis. Nat Med. 2019;25:1280–9. https://doi.org/10.1038/s41591-019-0512-5 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE131780 (2019).

Wu JQ, Dyer WB, Chrisp J, Belov L, Wang B, Saksena NK. Longitudinal microarray analysis of cell surface antigens on peripheral blood mononuclear cells from HIV+ individuals on highly active antiretroviral therapy. Retrovirology. 2008;5:24. https://doi.org/10.1186/1742-4690-5-24.

Zernecke A, Winkels H, Cochain C, Williams JW, Wolf D, Soehnlein O, et al. Meta-analysis of leukocyte diversity in atherosclerotic mouse aortas. Circ Res. 2020;127:402–26. https://doi.org/10.1161/CIRCRESAHA.120.316903.

Zhang L, Li Z, Skrzypczynska KM, Fang Q, Zhang W, O’Brien SA, et al. Single-cell analyses inform mechanisms of myeloid-targeted therapies in colon cancer. Cell. 2020;181:442–459.e29. https://doi.org/10.1016/j.cell.2020.03.048 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE146771 (2020).

Zhang Y, Zheng L, Zhang L, Hu X, Ren X, Zhang Z. Deep single-cell RNA sequencing data of individual T cells from treatment-naive colorectal cancer patients. Sci Data. 2019;6:131. https://doi.org/10.1038/s41597-019-0131-5 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE108989 (2018).

Zheng C, Zheng L, Yoo J-K, Guo H, Zhang Y, Guo X, et al. Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell. 2017;169:1342–1356.e16. https://doi.org/10.1016/j.cell.2017.05.035 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE98638 (2017).

Zheng GXY, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:1–12. https://doi.org/10.1038/ncomms14049 Gene expression omnibus (GEO) https://identifiers.org/geo:GSE93421 (2017).

Ley K, Kaplan RC, Hedrick C, Lanier LL, McNamara C, La Jolla Institute of Immunology. Combined protein and transcript single cell RNA sequencing in human peripheral blood mononuclear cells (human), NCBI GEO, 02.8, GSE205320. (2022).