# Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers.

Adrien Leite Pereira, Nicolas Tchitchek, Olivier Lambotte, Roger Le Grand, Antonio Cosma

<figure><img src="image_1.png" /></figure>

## To cite this version:

Adrien Leite Pereira, Nicolas Tchitchek, Olivier Lambotte, Roger Le Grand, Antonio Cosma. Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers. Frontiers in Immunology, 2019, 10, pp.1777. 10.3389/fimmu.2019.01777. hal-03961014

HAL Id: hal-03961014

https://hal.sorbonne-universite.fr/hal-03961014

Submitted on 5 Mar 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers

Adrien Leite Pereira¹†, Nicolas Tchitchek ¹*†, Olivier Lambotte ¹,²,³, Roger Le Grand¹ and Antonio Cosma¹

¹CEA-Université Paris Sud 11-INSERM U1184, Immunology of Viral Infections and Autoimmune Diseases, IDMIT Infrastructure, Fontenay-aux-Roses, France, ²APHP, Service de Médecine Interne-Immunologie Clinique, Hôpitaux Universitaires Paris Sud, Le Kremlin-Bicêtre, France, ³Université Paris Sud, Le Kremlin-Bicêtre, France

OPEN ACCESS

**Edited by:**
Helen Marie McGuire,
The University of Sydney, Australia

**Reviewed by:**
Juan C. Hernandez,
Cooperative University of
Colombia, Colombia
Nabila Seddiki,
Vaccine Research Institute
(VRI), France

*Correspondence:
Nicolas Tchitchek
nicolas.tchitchek@gmail.com

† These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Viral Immunology,
a section of the journal
Frontiers in Immunology

**Received:** 15 April 2019
**Accepted:** 15 July 2019
**Published:** 06 August 2019

**Citation:**
Leite Pereira A, Tchitchek N, Lambotte O, Le Grand R and Cosma A (2019) Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers. *Front. Immunol.* 10:1777. doi: 10.3389/fimmu.2019.01777

**Motivation:** Mass cytometry is a technique used to measure the intensity levels of proteins expressed by cells, at a single cell resolution. This technique is essential to characterize the phenotypes and functions of immune cell populations, but is currently limited to the measurement of 40 cell markers that restricts the characterization of complex diseases. However, algorithms and multi-tube cytometry techniques have been designed for combining phenotypic information obtained from different cytometric panels. The characterization of chronic HIV infection represents a good study case for multi-tube mass cytometry as this disease triggers a complex interactions network of more than 70 cell markers.

**Method:** We collected whole blood from non-viremic HIV-infected patients on combined antiretroviral therapies and healthy donors. Leukocytes from each individual were stained using three different mass cytometry panels, which consisted of 35, 32, and 33 cell markers. For each patient and using the CytoBackBone algorithm, we combined phenotypic information from three different antibody panels into a single cytometric profile, reaching a phenotypic resolution of 72 markers. These high-resolution cytometric profiles were analyzed using SPADE and viSNE algorithms to decipher the immune response to HIV.

**Results:** We detected an upregulation of several proteins in HIV-infected patients relative to healthy donors using our profiling of 72 cell markers. Among them, CD11a and CD11b were upregulated in PMNs, monocytes, mDCs, NK cells, and T cells. CD11b was also upregulated on pDCs. Other upregulated proteins included: CD38 on PMNs, monocytes, NK cells, basophils, B cells, and T cells; CD83 on monocytes, mDCs, B cells, and T cells; and TLR2, CD32, and CD64 on PMNs and monocytes. These results were validated using a mass cytometry panel of 25 cells markers.

**Impacts:** We demonstrate here that multi-tube cytometry can be applied to mass cytometry for exploring, at an unprecedented level of details, cell populations impacted
# Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers

Adrien Leite Pereira¹†, Nicolas Tchitchek ¹*†, Olivier Lambotte ¹,²,³, Roger Le Grand¹ and Antonio Cosma¹

¹CEA-Université Paris Sud 11-INSERM U1184, Immunology of Viral Infections and Autoimmune Diseases, IDMIT Infrastructure, Fontenay-aux-Roses, France, ²APHP, Service de Médecine Interne-Immunologie Clinique, Hôpitaux Universitaires Paris Sud, Le Kremlin-Bicêtre, France, ³Université Paris Sud, Le Kremlin-Bicêtre, France

OPEN ACCESS

**Edited by:**
Helen Marie McGuire,
The University of Sydney, Australia

**Reviewed by:**
Juan C. Hernandez,
Cooperative University of
Colombia, Colombia
Nabila Seddiki,
Vaccine Research Institute
(VRI), France

*Correspondence:
Nicolas Tchitchek
nicolas.tchitchek@gmail.com

† These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Viral Immunology,
a section of the journal
Frontiers in Immunology

**Received:** 15 April 2019
**Accepted:** 15 July 2019
**Published:** 06 August 2019

**Citation:**
Leite Pereira A, Tchitchek N, Lambotte O, Le Grand R and Cosma A (2019) Characterization of Leukocytes From HIV-ART Patients Using Combined Cytometric Profiles of 72 Cell Markers. *Front. Immunol.* 10:1777. doi: 10.3389/fimmu.2019.01777

**Motivation:** Mass cytometry is a technique used to measure the intensity levels of proteins expressed by cells, at a single cell resolution. This technique is essential to characterize the phenotypes and functions of immune cell populations, but is currently limited to the measurement of 40 cell markers that restricts the characterization of complex diseases. However, algorithms and multi-tube cytometry techniques have been designed for combining phenotypic information obtained from different cytometric panels. The characterization of chronic HIV infection represents a good study case for multi-tube mass cytometry as this disease triggers a complex interactions network of more than 70 cell markers.

**Method:** We collected whole blood from non-viremic HIV-infected patients on combined antiretroviral therapies and healthy donors. Leukocytes from each individual were stained using three different mass cytometry panels, which consisted of 35, 32, and 33 cell markers. For each patient and using the CytoBackBone algorithm, we combined phenotypic information from three different antibody panels into a single cytometric profile, reaching a phenotypic resolution of 72 markers. These high-resolution cytometric profiles were analyzed using SPADE and viSNE algorithms to decipher the immune response to HIV.

**Results:** We detected an upregulation of several proteins in HIV-infected patients relative to healthy donors using our profiling of 72 cell markers. Among them, CD11a and CD11b were upregulated in PMNs, monocytes, mDCs, NK cells, and T cells. CD11b was also upregulated on pDCs. Other upregulated proteins included: CD38 on PMNs, monocytes, NK cells, basophils, B cells, and T cells; CD83 on monocytes, mDCs, B cells, and T cells; and TLR2, CD32, and CD64 on PMNs and monocytes. These results were validated using a mass cytometry panel of 25 cells markers.

**Impacts:** We demonstrate here that multi-tube cytometry can be applied to mass cytometry for exploring, at an unprecedented level of details, cell populations impacted
by complex diseases. We showed that the monocyte and PMN populations were strongly affected by the HIV infection, as CD11a, CD11b, CD32, CD38, CD64, CD83, CD86, and TLR2 were upregulated in these populations. Overall, these results demonstrate that HIV induced a specific environment that similarly affected multiple immune cells.

**Keywords:** high-dimensional cytometry, cytometric profiles merging, HIV infection, chronic inflammation, biomarkers

## INTRODUCTION

Cytometry is an experimental technique used to measure the intensity levels of proteins expressed by cells, at a single cell resolution. However, none of these technologies can currently handle the simultaneous study of more than 50 markers (1). This restriction is particularly limiting to deeply characterize immune cell populations in complex inflammatory diseases, such as HIV infection.

Bioinformatics approaches have been successfully developed and applied to extend the number of observable cell markers per cytometric profiles through the merging of phenotypic information obtained from different cytometry profiles (generated from different cytometry panels). In such approaches, phenotypic information from the different profiles is combined into a single profile based on the expression of a common set of shared markers (2–4). Especially, the CytoBackBone algorithm, which is based on a Nearest Neighbor technique, can be used to extend our vision of complex immune diseases (5). However, no application of these bioinformatics approaches was performed in the context of infectious diseases using mass cytometry data.

HIV infection is still a major sexually transmissible disease. Although detected by different immune receptors, HIV is not cleared by the immune system. A global overview of the cell interaction network triggered in HIV infections is required to understand the immune system dysfunctions and to set up new therapeutic treatments.

HIV preferentially infects CD4⁺ T cells. In the long term, HIV induces a decrease of T cells number both in the blood and in the lymphoid organs. This situation then triggers a severe immune deficiency, that leads to Acquired Immunodeficiency Syndrome (AIDS) (6, 7). The monocyte compartment is also heavily impacted by HIV-infection. Indeed, this disease leads to an increase in the number of CD14low CD16high circulating monocytes (called non-classical monocytes), which is to the detriment of CD14high CD16low monocytes (8). Non-classical monocytes have been characterized as pro-inflammatory cells, as they are able to produce TNF-α, IL-1α, IL1-β, and IP-10 (9–11). In addition, they also overexpress costimulatory molecules such as CD80 and CD86, which again characterizes their activation state (12). Furthermore, the persistence of the virus in the organism results in chronic inflammation, which is also observed in treated HIV-infected patients (13, 14). Interestingly, the persistence of chronic inflammation could lead to an exhaustion of leukocytes (15). This exhaustion is characterized by a delay in the TLR4- or TLR7/8-dependent cytokine productions in monocytes and dendritic cells, as well as by the loss of functional and proliferative capacities of effector T cells. Moreover, several

negative regulators of immune activation, such as PD-1, LAG-3, Tim-3, and CTLA-4, are preferentially upregulated on T cells during HIV infection (16–18). Overall, a large set of immune cell populations are involved in this disease.

Many proteins involved in chronic inflammation in HIV-infected subjects have been described (19, 20). Indeed, the literature provides more than 70 cell markers involved in this disease (19, 20). Nonetheless, it is likely that a large number of cell markers are yet to be discovered. Also, a single-cell technology able to simultaneously monitor the expression of a complex network of proteins, involving more than 70 cell markers, is needed to better understand the inflammatory processes generated by HIV infection.

Here, we used the CytoBackBone algorithm (5) to characterize innate and adaptive immune cells obtained from treated HIV-infected patients and healthy donors using combined mass cytometric profiles of 72 cell markers. CytoBackBone is multi-tube cytometry algorithm based on a nearest neighbor approach. Especially, this algorithm uses the notion of acceptable and non-ambiguous neighbors to optimize the quality of the merging of cytometric profiles. Thanks to this notion, the merging performed by CytoBackBone are more stringent and noise-free compared to other approaches. Using this innovative multi-tube cytometry approach, we observed important phenotypic differences between treated HIV-infected patients and healthy donors in the whole leukocyte populations. Especially, we show that monocytes and polymorphonuclear cells (PMNs) are deeply impacted in treated HIV-infected patients. We confirmed these findings using cytometric profiles obtained from a single mass cytometry panel.

## MATERIALS AND METHODS

### Ethics

This experiment was approved by the Comité de Protection des Personnes (CPP) Ile de France VII, under protocol number PP 14-003. All subjects gave written informed consent to participate in this study.

### Whole Blood Collection

Whole blood from healthy ($n = 3$) and HIV-1 ART-treated non-viremic donors (undetectable plasma RNA, $n = 3$) was collected in lithium heparin tubes by the Etablissement Français du Sang (EFS, Hôpital Saint Louis, Paris, France) and Hôpital du Kremlin Bicêtre, respectively. Information concerning the gender, current age, contamination pathway, viral load, year of detection of the HIV infection, starting year of ARV treatment, and the type and the duration of treatment is provided for each HIV-infected
by complex diseases. We showed that the monocyte and PMN populations were strongly affected by the HIV infection, as CD11a, CD11b, CD32, CD38, CD64, CD83, CD86, and TLR2 were upregulated in these populations. Overall, these results demonstrate that HIV induced a specific environment that similarly affected multiple immune cells.

**Keywords:** high-dimensional cytometry, cytometric profiles merging, HIV infection, chronic inflammation, biomarkers

# INTRODUCTION

Cytometry is an experimental technique used to measure the intensity levels of proteins expressed by cells, at a single cell resolution. However, none of these technologies can currently handle the simultaneous study of more than 50 markers (1). This restriction is particularly limiting to deeply characterize immune cell populations in complex inflammatory diseases, such as HIV infection.

Bioinformatics approaches have been successfully developed and applied to extend the number of observable cell markers per cytometric profiles through the merging of phenotypic information obtained from different cytometry profiles (generated from different cytometry panels). In such approaches, phenotypic information from the different profiles is combined into a single profile based on the expression of a common set of shared markers (2–4). Especially, the CytoBackBone algorithm, which is based on a Nearest Neighbor technique, can be used to extend our vision of complex immune diseases (5). However, no application of these bioinformatics approaches was performed in the context of infectious diseases using mass cytometry data.

HIV infection is still a major sexually transmissible disease. Although detected by different immune receptors, HIV is not cleared by the immune system. A global overview of the cell interaction network triggered in HIV infections is required to understand the immune system dysfunctions and to set up new therapeutic treatments.

HIV preferentially infects CD4⁺ T cells. In the long term, HIV induces a decrease of T cells number both in the blood and in the lymphoid organs. This situation then triggers a severe immune deficiency, that leads to Acquired Immunodeficiency Syndrome (AIDS) (6, 7). The monocyte compartment is also heavily impacted by HIV-infection. Indeed, this disease leads to an increase in the number of CD14low CD16high circulating monocytes (called non-classical monocytes), which is to the detriment of CD14high CD16low monocytes (8). Non-classical monocytes have been characterized as pro-inflammatory cells, as they are able to produce TNF-α, IL-1α, IL1-β, and IP-10 (9–11). In addition, they also overexpress costimulatory molecules such as CD80 and CD86, which again characterizes their activation state (12). Furthermore, the persistence of the virus in the organism results in chronic inflammation, which is also observed in treated HIV-infected patients (13, 14). Interestingly, the persistence of chronic inflammation could lead to an exhaustion of leukocytes (15). This exhaustion is characterized by a delay in the TLR4- or TLR7/8-dependent cytokine productions in monocytes and dendritic cells, as well as by the loss of functional and proliferative capacities of effector T cells. Moreover, several

negative regulators of immune activation, such as PD-1, LAG-3, Tim-3, and CTLA-4, are preferentially upregulated on T cells during HIV infection (16–18). Overall, a large set of immune cell populations are involved in this disease.

Many proteins involved in chronic inflammation in HIV-infected subjects have been described (19, 20). Indeed, the literature provides more than 70 cell markers involved in this disease (19, 20). Nonetheless, it is likely that a large number of cell markers are yet to be discovered. Also, a single-cell technology able to simultaneously monitor the expression of a complex network of proteins, involving more than 70 cell markers, is needed to better understand the inflammatory processes generated by HIV infection.

Here, we used the CytoBackBone algorithm (5) to characterize innate and adaptive immune cells obtained from treated HIV-infected patients and healthy donors using combined mass cytometric profiles of 72 cell markers. CytoBackBone is multi-tube cytometry algorithm based on a nearest neighbor approach. Especially, this algorithm uses the notion of acceptable and non-ambiguous neighbors to optimize the quality of the merging of cytometric profiles. Thanks to this notion, the merging performed by CytoBackBone are more stringent and noise-free compared to other approaches. Using this innovative multi-tube cytometry approach, we observed important phenotypic differences between treated HIV-infected patients and healthy donors in the whole leukocyte populations. Especially, we show that monocytes and polymorphonuclear cells (PMNs) are deeply impacted in treated HIV-infected patients. We confirmed these findings using cytometric profiles obtained from a single mass cytometry panel.

# MATERIALS AND METHODS

## Ethics

This experiment was approved by the Comité de Protection des Personnes (CPP) Ile de France VII, under protocol number PP 14-003. All subjects gave written informed consent to participate in this study.

## Whole Blood Collection

Whole blood from healthy (n = 3) and HIV-1 ART-treated non-viremic donors (undetectable plasma RNA, n = 3) was collected in lithium heparin tubes by the Etablissement Français du Sang (EFS, Hôpital Saint Louis, Paris, France) and Hôpital du Kremlin Bicêtre, respectively. Information concerning the gender, current age, contamination pathway, viral load, year of detection of the HIV infection, starting year of ARV treatment, and the type and the duration of treatment is provided for each HIV-infected
<table><caption>TABLE 1 | Characteristics of HIV-infected patients and healthy donors.</caption><thead><tr><th>Individuals</th><th>Gender</th><th>Current age</th><th>Contamination</th><th>Number of CD4<sup>+</sup> T cells (cells/mm<sup>3</sup>)</th><th>CD4<sup>+</sup> T cells nadir (cells/mm<sup>3</sup>)</th><th>Viral load (copie/ml)</th><th>Detection (age of patient)</th><th>Treatment started - Updated treatment</th><th>Current treatment</th><th>HCV co-infection</th></tr></thead><tbody><tr><td>PAT-1</td><td>Male</td><td>54</td><td>Sexual (Homo)</td><td>543</td><td>204</td><td>&lt;40</td><td>1994 (32 years)</td><td>1997 - 2015</td><td>Tivicay Truvada</td><td>No</td></tr><tr><td>PAT-2</td><td>Female</td><td>54</td><td>Sexual (Hetero)</td><td>1376</td><td>167</td><td>&lt;40</td><td>1991 (26 years)</td><td>1994 - 2014</td><td>Isentress Kivexa</td><td>No</td></tr><tr><td>PAT-3</td><td>Female</td><td>70</td><td>Sexual (Hetero)</td><td>404</td><td>170</td><td>&lt;40</td><td>2001 (55 years)</td><td>2008 - 2010</td><td>Isentress Truvada</td><td>No</td></tr><tr><td>PAT-4</td><td>Male</td><td>51</td><td>Sexual (Homo)</td><td>1,451</td><td>130</td><td>&lt;40</td><td>1990 (24 years)</td><td>1997 - 2013</td><td>Kivexa Norvir Reyataz</td><td>No</td></tr><tr><td>PAT-5</td><td>Male</td><td>58</td><td>Sexual (Hetero)</td><td>324</td><td>13</td><td>&lt;40</td><td>2002 (43 years)</td><td>2002 - 2016</td><td>Kivexa Norvir Reyataz</td><td>No</td></tr><tr><td>PAT-6</td><td>Male</td><td>48</td><td>Sexual (Homo)</td><td>744</td><td>243</td><td>&lt;40</td><td>1991 (22 years)</td><td>1998 - 2015</td><td>Eviplera</td><td>No</td></tr><tr><td>HEA-1</td><td>Male</td><td>48</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-2</td><td>Male</td><td>63</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-3</td><td>Female</td><td>28</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-4</td><td>Male</td><td>59</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-5</td><td>Male</td><td>56</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-6</td><td>Male</td><td>56</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-7</td><td>Female</td><td>59</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-8</td><td>Female</td><td>25</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-9</td><td>Male</td><td>27</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Information about the gender, current age, contamination pathway, viral load, year of detection of the HIV infection, starting year of the ARV treatment, updated year of the ARV treatment, type treatment, the duration of treatments, and the HCV co-infection status were provided for each HIV-infected patient. Information about the gender and the age were provided for each healthy donor. Information presented in this table corresponds to information obtained at the time of the blood samplings.

patient in Table 1. The gender and current age of each healthy donor are also provided.

To validate the results with a single cytometric panel, blood from the three previous HIV-infected subjects, three new HIV-infected patients (n = 6), and six new healthy subjects (n = 6) was collected. Information concerning the gender and the current age is provided for each HIV-infected patient and each healthy donor in Table 1. In addition, information concerning contamination pathway, viral load, year of detection of the HIV infection, starting year of ARV treatment, and type and duration of treatment is also provided for each HIV-infected patient.

## Sample Processing for Mass Cytometry Data

Blood samples were processed according to a previously described protocol (21). The cells (from 1 ml blood) were mixed with 10 ml fixation mixture (FM) in 50-ml plastic tubes and incubated for 10 min at 4°C. After centrifugation at 800 x g for 5 min at room temperature (RT), red cells were lysed by adding 10 ml Milli-Q water at RT for 20 min, without agitation. After two

washes with 1X DPBS, cells were counted and stored at -80°C in FM at a final concentration of 15 × 10⁶ cells/ml and distributed into aliquots containing 3 × 10⁶ cells. FM used to fix and store the cells was prepared the day before the experiments and conserved at 4°C. The 5% formaldehyde FM solution was prepared from 36% paraformaldehyde (VWR BDH Prolabo, Fontenay-sous-Bois) and contained 18.5% glycerol (Sigma-Aldrich, Lyon, France) in 1X-Dulbecco's phosphate buffered saline (DPBS), without CaCl₂ or MgCl₂, pH 7.4 (Gibco by life Technologies, Villebon-Sur-Yvette, France). This solution allowed freezing and recovery of all blood leukocytes, especially polymorphonuclear cells, which are highly labile and cryopreservation-sensitive. Healthy and HIV-infected samples used for the multi-tube 72-marker experiment were cryopreserved for a maximum of 12 days.

## Staining Protocols for Mass Cytometry Data

For each sample, 3 × 10⁶ cryopreserved fixed cells were washed twice with staining buffer [PBS/0.5% BSA, prepared by mixing
<table><caption>**TABLE 1** | Characteristics of HIV-infected patients and healthy donors.</caption><thead><tr><th>Individuals</th><th>Gender</th><th>Current age</th><th>Contamination</th><th>Number of CD4<sup>+</sup> T cells (cells/mm<sup>3</sup>)</th><th>CD4<sup>+</sup> T cells nadir (cells/mm<sup>3</sup>)</th><th>Viral load (copie/ml)</th><th>Detection (age of patient)</th><th>Treatment started - Updated treatment</th><th>Current treatment</th><th>HCV co-infection</th></tr></thead><tbody><tr><td>PAT-1</td><td>Male</td><td>54</td><td>Sexual (Homo)</td><td>543</td><td>204</td><td>&lt;40</td><td>1994 (32 years)</td><td>1997 - 2015</td><td>Tivicay Truvada</td><td>No</td></tr><tr><td>PAT-2</td><td>Female</td><td>54</td><td>Sexual (Hetero)</td><td>1376</td><td>167</td><td>&lt;40</td><td>1991 (26 years)</td><td>1994 - 2014</td><td>Isentress Kivexa</td><td>No</td></tr><tr><td>PAT-3</td><td>Female</td><td>70</td><td>Sexual (Hetero)</td><td>404</td><td>170</td><td>&lt;40</td><td>2001 (55 years)</td><td>2008 - 2010</td><td>Isentress Truvada</td><td>No</td></tr><tr><td>PAT-4</td><td>Male</td><td>51</td><td>Sexual (Homo)</td><td>1,451</td><td>130</td><td>&lt;40</td><td>1990 (24 years)</td><td>1997 - 2013</td><td>Kivexa Norvir Reyataz</td><td>No</td></tr><tr><td>PAT-5</td><td>Male</td><td>58</td><td>Sexual (Hetero)</td><td>324</td><td>13</td><td>&lt;40</td><td>2002 (43 years)</td><td>2002 - 2016</td><td>Kivexa Norvir Reyataz</td><td>No</td></tr><tr><td>PAT-6</td><td>Male</td><td>48</td><td>Sexual (Homo)</td><td>744</td><td>243</td><td>&lt;40</td><td>1991 (22 years)</td><td>1998 - 2015</td><td>Eviplera</td><td>No</td></tr><tr><td>HEA-1</td><td>Male</td><td>48</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-2</td><td>Male</td><td>63</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-3</td><td>Female</td><td>28</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-4</td><td>Male</td><td>59</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-5</td><td>Male</td><td>56</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-6</td><td>Male</td><td>56</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-7</td><td>Female</td><td>59</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-8</td><td>Female</td><td>25</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>HEA-9</td><td>Male</td><td>27</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Information about the gender, current age, contamination pathway, viral load, year of detection of the HIV infection, starting year of the ARV treatment, updated year of the ARV treatment, type treatment, the duration of treatments, and the HCV co-infection status were provided for each HIV-infected patient. Information about the gender and the age were provided for each healthy donor. Information presented in this table corresponds to information obtained at the time of the blood samplings.

patient in **Table 1**. The gender and current age of each healthy donor are also provided.

To validate the results with a single cytometric panel, blood from the three previous HIV-infected subjects, three new HIV-infected patients (n = 6), and six new healthy subjects (n = 6) was collected. Information concerning the gender and the current age is provided for each HIV-infected patient and each healthy donor in **Table 1**. In addition, information concerning contamination pathway, viral load, year of detection of the HIV infection, starting year of ARV treatment, and type and duration of treatment is also provided for each HIV-infected patient.

## Sample Processing for Mass Cytometry Data

Blood samples were processed according to a previously described protocol (21). The cells (from 1 ml blood) were mixed with 10 ml fixation mixture (FM) in 50-ml plastic tubes and incubated for 10 min at 4°C. After centrifugation at 800 x g for 5 min at room temperature (RT), red cells were lysed by adding 10 ml Milli-Q water at RT for 20 min, without agitation. After two

washes with 1X DPBS, cells were counted and stored at -80°C in FM at a final concentration of 15 × 10⁶ cells/ml and distributed into aliquots containing 3 × 10⁶ cells. FM used to fix and store the cells was prepared the day before the experiments and conserved at 4°C. The 5% formaldehyde FM solution was prepared from 36% paraformaldehyde (VWR BDH Prolabo, Fontenay-sous-Bois) and contained 18.5% glycerol (Sigma-Aldrich, Lyon, France) in 1X-Dulbecco's phosphate buffered saline (DPBS), without CaCl₂ or MgCl₂, pH 7.4 (Gibco by life Technologies, Villebon-Sur-Yvette, France). This solution allowed freezing and recovery of all blood leukocytes, especially polymorphonuclear cells, which are highly labile and cryopreservation-sensitive. Healthy and HIV-infected samples used for the multi-tube 72-marker experiment were cryopreserved for a maximum of 12 days.

## Staining Protocols for Mass Cytometry Data

For each sample, 3 × 10⁶ cryopreserved fixed cells were washed twice with staining buffer [PBS/0.5% BSA, prepared by mixing
1X DPBS modified (Gibco by Life Technologies) with 0.5% BSA (Sigma-Aldrich, Lyon, France)] and labeled with conjugated antibodies according to the following procedure. Cells were incubated at 4°C for 30 min with a mix of the metal-labeled surface antibodies (Abs) in staining buffer. After two washes with 1X DPBS, cells were incubated in fixation solution (PBS/1.6% PFA, prepared by diluting 16% paraformaldehyde (PFA; Electron Microscopy Sciences Hartfield, USA) in DPBS 10X and Milli-Q water) at RT for 20 min and permeabilized with 1X Perm/wash (BD) buffer at RT for 10 min. Staining with metal-labeled intracellular Abs and an iridium nucleic acid intercalator in 1X Perm/Wash was carried out as for extracellular staining, but without a new permeabilization step. Cells were stored overnight with 0.1 µM iridium nucleic acid intercalator in fixation mixture. The following day, cells were washed with Milli-Q water, resuspended in 1 ml Milli-Q water and filtered using a 35-µm nylon mesh cell strainer (BD Biosciences), before the addition of EQ Four-Element Calibration Beads (Fluidigm, San Francisco, USA), according to the manufacturer's instructions. Acquisition of each sample was performed two times in succession on a CyTOF or Helios instrument (Fluidigm). No autosampler was used for sample acquisition. Moreover, the dual sample loop system was used for sample acquisition performed on the CyTOF. CyTOF settings were parametrized following the quality control of the instrument.

## Processing of Mass Cytometry Data

Data were acquired using EQTM Four-Element Calibration Beads, normalized using Rachel Finck's MATLAB normalizer, and concatenated using the FCS file concatenation tool (Cytobank, Mountain View, CA). Cytometric profiles were manually gated to exclude the EQTM Four-Element Calibration Beads, select singlets, and gate out non-specific background generated by metal-conjugated Ab binding to eosinophils (22), as represented in **Supplementary Figure 1**. The arcsinh (cofactor = 5) transformation was applied on marker expressions.

## Merging of Cytometry Profiles Using CytoBackBone

The CytoBackBone algorithm (5) was used to merge cytometric profiles from the different panels. In details, CytoBackBone is an nearest neighbor-based algorithm that combines phenotypic information of different cytometric profiles obtained from different cytometry panels. In our approach, CytoBackBone combines marker expression information of two cells from different cytometric profiles if, and only if, these two cells are acceptable and non-ambiguous nearest neighbors. In details, each time two cells from different cytometric profiles are identified as acceptable and non-ambiguous nearest neighbors, CytoBackBone merges them and a new cell with a combined phenotype is created. This new cell is added to the resulting cytometric profile and its phenotype corresponds to the average marker expression for the set of common markers. For the non-common markers (i.e., markers that do not belong to the backbone composition), the phenotype of this cell corresponds to the specific marker expression. CytoBackBone is implemented in R and the source code can be found at: https://

github.com/tchitchek-lab/CytoBackBone. Each merge realized by CytoBackBone was calculated based on the set of common markers between the cytometric profiles to be merged (CD1c, CD3, CD11c, CD14, CD16, CD19, CD32, CD64, CD66, CD86, CD123, CD141, Granzyme B, and HLADR). A distance threshold of 3 was used for all merging.

## Analysis of Merged Mass Cytometric Profiles

SPADE analyses were generated using the public version of the SPADE R package. viSNE representations were generated based on the Barnes-Hut implementation of the t-SNE algorithm. Dot-plot visualizations were performed using Cytobank software (Mountain View, CA). SPADE analyses were set up and carried out with a down-sampling parameter of 5%. A random pre-downsampling was used to select 67,000 cells from each sample (67,000 corresponded to the number of cells contained in the smallest sample) before the SPADE analysis. Full upsampling was eventually performed after the cell cluster identification by SPADE. Categorical marker heatmaps were generated using the SPADEVizR package (23). Marker expression categories represented in the heatmap were generated based on the intensity expression range of each marker. Marker expression ranges were calculated based on the 5 and 95th percentiles of marker expression. Then, each range was divided into five uniform categories. The categorization of marker expression was computed based on the means of the individual SPADE-expression medians of each marker. These categories represented negative, low, medium, high, and very high marker expression using a color scale from white to dark red, which was used to produce the heatmap. Cell clusters with <50 cells (for the whole dataset) were annotated as unassigned and were colored in gray in the categorical heatmap representation. Indeed, the phenotypes of these clusters cannot be evaluated based on the mean expression of marker expressions. Hierarchical clustering of markers and clusters was performed using the Euclidian distance and was based on the complete linkage method. Statistical analyses of cell cluster abundances were obtained using a non-parametric permutation test (available in the Deducer R package).

# RESULTS

## Phenotypic Characterization of Leukocytes From HIV-Infected Patients and Healthy Donor Using Merged Cytometric Profiles of 72 Cell Markers

To characterize the phenotypes of leukocytes obtained from HIV-infected patients and healthy donors, we used three different mass cytometry panels. Phenotypic information from the different cytometry profiles was combined using the CytoBackBone algorithm. Such multidimensional phenotype characterization was performed based on the simultaneous expression of 72 markers.

We collected whole blood from three non-viremic HIV-infected patients on combined antiretroviral therapies (named
1X DPBS modified (Gibco by Life Technologies) with 0.5% BSA (Sigma-Aldrich, Lyon, France)] and labeled with conjugated antibodies according to the following procedure. Cells were incubated at 4°C for 30 min with a mix of the metal-labeled surface antibodies (Abs) in staining buffer. After two washes with 1X DPBS, cells were incubated in fixation solution (PBS/1.6% PFA, prepared by diluting 16% paraformaldehyde (PFA; Electron Microscopy Sciences Hartfield, USA) in DPBS 10X and Milli-Q water) at RT for 20 min and permeabilized with 1X Perm/wash (BD) buffer at RT for 10 min. Staining with metal-labeled intracellular Abs and an iridium nucleic acid intercalator in 1X Perm/Wash was carried out as for extracellular staining, but without a new permeabilization step. Cells were stored overnight with 0.1 µM iridium nucleic acid intercalator in fixation mixture. The following day, cells were washed with Milli-Q water, resuspended in 1 ml Milli-Q water and filtered using a 35-µm nylon mesh cell strainer (BD Biosciences), before the addition of EQ Four-Element Calibration Beads (Fluidigm, San Francisco, USA), according to the manufacturer's instructions. Acquisition of each sample was performed two times in succession on a CyTOF or Helios instrument (Fluidigm). No autosampler was used for sample acquisition. Moreover, the dual sample loop system was used for sample acquisition performed on the CyTOF. CyTOF settings were parametrized following the quality control of the instrument.

## Processing of Mass Cytometry Data

Data were acquired using EQTM Four-Element Calibration Beads, normalized using Rachel Finck's MATLAB normalizer, and concatenated using the FCS file concatenation tool (Cytobank, Mountain View, CA). Cytometric profiles were manually gated to exclude the EQTM Four-Element Calibration Beads, select singlets, and gate out non-specific background generated by metal-conjugated Ab binding to eosinophils (22), as represented in **Supplementary Figure 1**. The arcsinh (cofactor = 5) transformation was applied on marker expressions.

## Merging of Cytometry Profiles Using CytoBackBone

The CytoBackBone algorithm (5) was used to merge cytometric profiles from the different panels. In details, CytoBackBone is an nearest neighbor-based algorithm that combines phenotypic information of different cytometric profiles obtained from different cytometry panels. In our approach, CytoBackBone combines marker expression information of two cells from different cytometric profiles if, and only if, these two cells are acceptable and non-ambiguous nearest neighbors. In details, each time two cells from different cytometric profiles are identified as acceptable and non-ambiguous nearest neighbors, CytoBackBone merges them and a new cell with a combined phenotype is created. This new cell is added to the resulting cytometric profile and its phenotype corresponds to the average marker expression for the set of common markers. For the non-common markers (i.e., markers that do not belong to the backbone composition), the phenotype of this cell corresponds to the specific marker expression. CytoBackBone is implemented in R and the source code can be found at: https://

github.com/tchitchek-lab/CytoBackBone. Each merge realized by CytoBackBone was calculated based on the set of common markers between the cytometric profiles to be merged (CD1c, CD3, CD11c, CD14, CD16, CD19, CD32, CD64, CD66, CD86, CD123, CD141, Granzyme B, and HLADR). A distance threshold of 3 was used for all merging.

## Analysis of Merged Mass Cytometric Profiles

SPADE analyses were generated using the public version of the SPADE R package. viSNE representations were generated based on the Barnes-Hut implementation of the t-SNE algorithm. Dot-plot visualizations were performed using Cytobank software (Mountain View, CA). SPADE analyses were set up and carried out with a down-sampling parameter of 5%. A random pre-downsampling was used to select 67,000 cells from each sample (67,000 corresponded to the number of cells contained in the smallest sample) before the SPADE analysis. Full upsampling was eventually performed after the cell cluster identification by SPADE. Categorical marker heatmaps were generated using the SPADEVizR package (23). Marker expression categories represented in the heatmap were generated based on the intensity expression range of each marker. Marker expression ranges were calculated based on the 5 and 95th percentiles of marker expression. Then, each range was divided into five uniform categories. The categorization of marker expression was computed based on the means of the individual SPADE-expression medians of each marker. These categories represented negative, low, medium, high, and very high marker expression using a color scale from white to dark red, which was used to produce the heatmap. Cell clusters with <50 cells (for the whole dataset) were annotated as unassigned and were colored in gray in the categorical heatmap representation. Indeed, the phenotypes of these clusters cannot be evaluated based on the mean expression of marker expressions. Hierarchical clustering of markers and clusters was performed using the Euclidian distance and was based on the complete linkage method. Statistical analyses of cell cluster abundances were obtained using a non-parametric permutation test (available in the Deducer R package).

# RESULTS

## Phenotypic Characterization of Leukocytes From HIV-Infected Patients and Healthy Donor Using Merged Cytometric Profiles of 72 Cell Markers

To characterize the phenotypes of leukocytes obtained from HIV-infected patients and healthy donors, we used three different mass cytometry panels. Phenotypic information from the different cytometry profiles was combined using the CytoBackBone algorithm. Such multidimensional phenotype characterization was performed based on the simultaneous expression of 72 markers.

We collected whole blood from three non-viremic HIV-infected patients on combined antiretroviral therapies (named
<table><caption>**TABLE 2** | Antibodies and cell markers used to create combined mass cytometry profiles including 72 cell markers.</caption><thead><tr><th rowspan="2">Metals</th><th colspan="2">Panel #A</th><th colspan="2">Panel #B</th><th colspan="2">Panel #C</th></tr><tr><th>Antibodies</th><th>Clones</th><th>Antibodies</th><th>Clones</th><th>Antibodies</th><th>Clones</th></tr></thead><tbody><tr><td>Pr141</td><td>CD66</td><td>TET2</td><td>CD66</td><td>TET2</td><td>CD66</td><td>TET2</td></tr><tr><td>Nd142</td><td>HLADR</td><td>L243 (G46-6)</td><td>HLADR</td><td>L243 (G46-6)</td><td>HLADR</td><td>L243 (G46-6)</td></tr><tr><td>Nd143</td><td>CD3</td><td>UCHT1</td><td>CD3</td><td>UCHT1</td><td>CD3</td><td>UCHT1</td></tr><tr><td>Nd144</td><td>CD64</td><td>10.1.1</td><td>CD64</td><td>10.1.1</td><td>CD64</td><td>10.1.1</td></tr><tr><td>Nd145</td><td>CD8a</td><td>37006</td><td>CCR7</td><td>150503</td><td>CD209</td><td>DCN47.5</td></tr><tr><td>Nd146</td><td>IL-6</td><td>MQ2-13A5</td><td>NF-ATC</td><td>7A6</td><td>CD62L</td><td>DREG-56</td></tr><tr><td>Sm147</td><td>Granzyme A</td><td>CTLA-3</td><td>CD335</td><td>9 E 2</td><td>CD45</td><td>D058-1283</td></tr><tr><td>Nd148</td><td>IL-1β</td><td>H1b-98</td><td>IL-4</td><td>7A3-3</td><td>CD137</td><td>4B4-1</td></tr><tr><td>Sm149</td><td>CD14</td><td>M5E2</td><td>CD14</td><td>M5E2</td><td>CD14</td><td>M5E2</td></tr><tr><td>Nd150</td><td>CD123</td><td>7G3</td><td>CD123</td><td>7G3</td><td>CD123</td><td>7G3</td></tr><tr><td>Eu151</td><td>IL-8</td><td>NAPII</td><td>CD107a</td><td>H4A3</td><td>-</td><td>-</td></tr><tr><td>Sm152</td><td>CD16</td><td>B73.1</td><td>CD16</td><td>B73.1</td><td>CD16</td><td>B73.1</td></tr><tr><td>Eu153</td><td>CD23</td><td>M-L233</td><td>CD154</td><td>TRAP1</td><td>CXCR3 (CD183)</td><td>1C6</td></tr><tr><td>Sm154</td><td>CD86</td><td>2,331 (FUN-1)</td><td>CD86</td><td>2,331 (FUN-1)</td><td>CD86</td><td>2331 (FUN-1)</td></tr><tr><td>Gd155</td><td>CD32</td><td>2 E1</td><td>CD32</td><td>2 E1</td><td>CD32</td><td>2 E1</td></tr><tr><td>Gd156</td><td>MIP-1β</td><td>D21-1351</td><td>CD54</td><td>LB-2</td><td>CD56</td><td>AF12-7H3</td></tr><tr><td>Gd158</td><td>IP-10</td><td>6D4</td><td>IL-2</td><td>N7.48A</td><td>ITGβ7</td><td>FIB504</td></tr><tr><td>Tb159</td><td>TNF-α</td><td>MAb11</td><td>-</td><td>-</td><td>Bcl-6</td><td>K112-91</td></tr><tr><td>Gd160</td><td>IL-1α</td><td>364/3B3-14</td><td>CD69</td><td>FN50</td><td>CD83</td><td>HB15E</td></tr><tr><td>Dy161</td><td>CD141</td><td>MAB3947</td><td>CD141</td><td>MAB3947</td><td>CD141</td><td>MAB3947</td></tr><tr><td>Dy162</td><td>IL-12</td><td>C8,6</td><td>Ki67</td><td>B56</td><td>CD279</td><td>MIH4</td></tr><tr><td>Dy163</td><td>CD1c</td><td>AF5916</td><td>CD1c</td><td>AF5916</td><td>CD1c</td><td>AF5916</td></tr><tr><td>Dy164</td><td>CXCR4</td><td>12G5</td><td>CD25</td><td>BC96</td><td>CD127</td><td>HIL-7R-M21</td></tr><tr><td>Ho165</td><td>TLR2</td><td>REA109</td><td>CD11a</td><td>HI111</td><td>IL-10R</td><td>3F9</td></tr><tr><td>Er166</td><td>CCR5</td><td>3A9</td><td>CD11b</td><td>ICRF44</td><td>CD27</td><td>M-T271</td></tr><tr><td>Er167</td><td>CD28</td><td>CD28.2</td><td>CD38</td><td>AT1</td><td>IgG2c</td><td>A23-1</td></tr><tr><td>Er168</td><td>CD11c</td><td>B-ly6</td><td>CD11c</td><td>B-ly6</td><td>CD11c</td><td>B-ly6</td></tr><tr><td>Tm169</td><td>IFN-α</td><td>LT27:295</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Er170</td><td>CD45RA</td><td>T6D11</td><td>IgG2b</td><td>27-35</td><td>CCL5</td><td>2D5</td></tr><tr><td>Yb171</td><td>IFN-γ</td><td>25723</td><td>IL-10</td><td>JES3-9D7</td><td>IgG1 K</td><td>MOPC-21</td></tr><tr><td>Yb172</td><td>CD4</td><td>L200</td><td>MyD88</td><td>RB2101</td><td>IL-2RA</td><td>24212</td></tr><tr><td>Yb173</td><td>Granzyme B</td><td>GB11</td><td>Granzyme B</td><td>GB11</td><td>Granzyme B</td><td>GB11</td></tr><tr><td>Yb174</td><td>CD19</td><td>HIB19</td><td>CD19</td><td>HIB19</td><td>CD19</td><td>HIB19</td></tr><tr><td>Yb175</td><td>IL-1RA</td><td>AS17</td><td>-</td><td>-</td><td>IgG2a</td><td>R35-95</td></tr><tr><td>Yb176</td><td>MCP-1</td><td>5D3-F7</td><td>Perforin</td><td>dG9-DTAG9</td><td>nfKB</td><td>K10-892.12.50</td></tr><tr><td>Ir191</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Ir193</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

The metal isotopes, markers, and antibody clones are indicated for each panel. Backbone markers are indicated in bold.

PAT-1 to PAT-3) and three healthy donors (named HEA-1 to HEA-3), as indicated in **Table 1**. Following red cell lysis, cells from each sample were stained using three mass cytometry panels, named #A, #B, #C, which consisted of 35, 32, and 33 markers, respectively (**Table 2**). These three panels shared 14 common markers.

For each patient, phenotypic information from various mass cytometry panels was merged using the CytoBackBone algorithm. Thus, for each of the six subjects, we generated a merged cytometric profile of 72 cell markers.

For each cytometric profile, the number of cells obtained before and after merging is shown in **Table 3**. The average of cells present in merged files was 142,265 cells. Furthermore, cytometric profiles related to healthy subjects were composed of 157,083 merged cells (ranging from 103,088 to 229,058 cells), whereas cytometric profiles related to HIV patients were composed of 127,446 merged cells (ranging from 67,280 to 211,223 cells). Among the three cytometric profiles to combine for each individual, more of 97% of cells included in cytometric profiles having the lowest number of cells
<table><caption>**TABLE 2** | Antibodies and cell markers used to create combined mass cytometry profiles including 72 cell markers.</caption><thead><tr><th rowspan="2">Metals</th><th colspan="2">Panel #A</th><th colspan="2">Panel #B</th><th colspan="2">Panel #C</th></tr><tr><th>Antibodies</th><th>Clones</th><th>Antibodies</th><th>Clones</th><th>Antibodies</th><th>Clones</th></tr></thead><tbody><tr><td>Pr141</td><td>CD66</td><td>TET2</td><td>CD66</td><td>TET2</td><td>CD66</td><td>TET2</td></tr><tr><td>Nd142</td><td>HLADR</td><td>L243 (G46-6)</td><td>HLADR</td><td>L243 (G46-6)</td><td>HLADR</td><td>L243 (G46-6)</td></tr><tr><td>Nd143</td><td>CD3</td><td>UCHT1</td><td>CD3</td><td>UCHT1</td><td>CD3</td><td>UCHT1</td></tr><tr><td>Nd144</td><td>CD64</td><td>10.1.1</td><td>CD64</td><td>10.1.1</td><td>CD64</td><td>10.1.1</td></tr><tr><td>Nd145</td><td>CD8a</td><td>37006</td><td>CCR7</td><td>150503</td><td>CD209</td><td>DCN47.5</td></tr><tr><td>Nd146</td><td>IL-6</td><td>MQ2-13A5</td><td>NF-ATC</td><td>7A6</td><td>CD62L</td><td>DREG-56</td></tr><tr><td>Sm147</td><td>Granzyme A</td><td>CTLA-3</td><td>CD335</td><td>9 E 2</td><td>CD45</td><td>D058-1283</td></tr><tr><td>Nd148</td><td>IL-1β</td><td>H1b-98</td><td>IL-4</td><td>7A3-3</td><td>CD137</td><td>4B4-1</td></tr><tr><td>Sm149</td><td>CD14</td><td>M5E2</td><td>CD14</td><td>M5E2</td><td>CD14</td><td>M5E2</td></tr><tr><td>Nd150</td><td>CD123</td><td>7G3</td><td>CD123</td><td>7G3</td><td>CD123</td><td>7G3</td></tr><tr><td>Eu151</td><td>IL-8</td><td>NAPII</td><td>CD107a</td><td>H4A3</td><td>-</td><td>-</td></tr><tr><td>Sm152</td><td>CD16</td><td>B73.1</td><td>CD16</td><td>B73.1</td><td>CD16</td><td>B73.1</td></tr><tr><td>Eu153</td><td>CD23</td><td>M-L233</td><td>CD154</td><td>TRAP1</td><td>CXCR3 (CD183)</td><td>1C6</td></tr><tr><td>Sm154</td><td>CD86</td><td>2,331 (FUN-1)</td><td>CD86</td><td>2,331 (FUN-1)</td><td>CD86</td><td>2331 (FUN-1)</td></tr><tr><td>Gd155</td><td>CD32</td><td>2 E1</td><td>CD32</td><td>2 E1</td><td>CD32</td><td>2 E1</td></tr><tr><td>Gd156</td><td>MIP-1β</td><td>D21-1351</td><td>CD54</td><td>LB-2</td><td>CD56</td><td>AF12-7H3</td></tr><tr><td>Gd158</td><td>IP-10</td><td>6D4</td><td>IL-2</td><td>N7.48A</td><td>ITGβ7</td><td>FIB504</td></tr><tr><td>Tb159</td><td>TNF-α</td><td>MAb11</td><td>-</td><td>-</td><td>Bcl-6</td><td>K112-91</td></tr><tr><td>Gd160</td><td>IL-1α</td><td>364/3B3-14</td><td>CD69</td><td>FN50</td><td>CD83</td><td>HB15E</td></tr><tr><td>Dy161</td><td>CD141</td><td>MAB3947</td><td>CD141</td><td>MAB3947</td><td>CD141</td><td>MAB3947</td></tr><tr><td>Dy162</td><td>IL-12</td><td>C8,6</td><td>Ki67</td><td>B56</td><td>CD279</td><td>MIH4</td></tr><tr><td>Dy163</td><td>CD1c</td><td>AF5916</td><td>CD1c</td><td>AF5916</td><td>CD1c</td><td>AF5916</td></tr><tr><td>Dy164</td><td>CXCR4</td><td>12G5</td><td>CD25</td><td>BC96</td><td>CD127</td><td>HIL-7R-M21</td></tr><tr><td>Ho165</td><td>TLR2</td><td>REA109</td><td>CD11a</td><td>HI111</td><td>IL-10R</td><td>3F9</td></tr><tr><td>Er166</td><td>CCR5</td><td>3A9</td><td>CD11b</td><td>ICRF44</td><td>CD27</td><td>M-T271</td></tr><tr><td>Er167</td><td>CD28</td><td>CD28.2</td><td>CD38</td><td>AT1</td><td>IgG2c</td><td>A23-1</td></tr><tr><td>Er168</td><td>CD11c</td><td>B-ly6</td><td>CD11c</td><td>B-ly6</td><td>CD11c</td><td>B-ly6</td></tr><tr><td>Tm169</td><td>IFN-α</td><td>LT27:295</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Er170</td><td>CD45RA</td><td>T6D11</td><td>IgG2b</td><td>27-35</td><td>CCL5</td><td>2D5</td></tr><tr><td>Yb171</td><td>IFN-γ</td><td>25723</td><td>IL-10</td><td>JES3-9D7</td><td>IgG1 K</td><td>MOPC-21</td></tr><tr><td>Yb172</td><td>CD4</td><td>L200</td><td>MyD88</td><td>RB2101</td><td>IL-2RA</td><td>24212</td></tr><tr><td>Yb173</td><td>Granzyme B</td><td>GB11</td><td>Granzyme B</td><td>GB11</td><td>Granzyme B</td><td>GB11</td></tr><tr><td>Yb174</td><td>CD19</td><td>HIB19</td><td>CD19</td><td>HIB19</td><td>CD19</td><td>HIB19</td></tr><tr><td>Yb175</td><td>IL-1RA</td><td>AS17</td><td>-</td><td>-</td><td>IgG2a</td><td>R35-95</td></tr><tr><td>Yb176</td><td>MCP-1</td><td>5D3-F7</td><td>Perforin</td><td>dG9-DTAG9</td><td>nfKB</td><td>K10-892.12.50</td></tr><tr><td>Ir191</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Ir193</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

The metal isotopes, markers, and antibody clones are indicated for each panel. Backbone markers are indicated in bold.

PAT-1 to PAT-3) and three healthy donors (named HEA-1 to HEA-3), as indicated in **Table 1**. Following red cell lysis, cells from each sample were stained using three mass cytometry panels, named #A, #B, #C, which consisted of 35, 32, and 33 markers, respectively (**Table 2**). These three panels shared 14 common markers.

For each patient, phenotypic information from various mass cytometry panels was merged using the CytoBackBone algorithm. Thus, for each of the six subjects, we generated a merged cytometric profile of 72 cell markers.

For each cytometric profile, the number of cells obtained before and after merging is shown in **Table 3**. The average of cells present in merged files was 142,265 cells. Furthermore, cytometric profiles related to healthy subjects were composed of 157,083 merged cells (ranging from 103,088 to 229,058 cells), whereas cytometric profiles related to HIV patients were composed of 127,446 merged cells (ranging from 67,280 to 211,223 cells). Among the three cytometric profiles to combine for each individual, more of 97% of cells included in cytometric profiles having the lowest number of cells
<table><caption>**TABLE 3** | Number of events and cells obtained, before and after cell merging, for samples used to characterize specific cell populations from healthy subjects and HIV-infected patients.</caption><thead><tr><th>Sample ID</th><th>Stimulation</th><th>Profile</th><th>Number of events</th><th>Number of cells</th><th>Number of excluded cells</th></tr></thead><tbody><tr><td>HEA-1</td><td>CTRL</td><td>#A</td><td>443,527</td><td>237,039</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#B</td><td>321,545</td><td>192,201</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#C</td><td>245,550</td><td>139,218</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>139,104</td><td>151,146</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#A</td><td>464,263</td><td>249,773</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#B</td><td>463,584</td><td>249,003</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#C</td><td>413,188</td><td>229,388</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>229,058</td><td>40,990</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#A</td><td>304,750</td><td>128,604</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#B</td><td>243,011</td><td>112,246</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#C</td><td>219,468</td><td>103,360</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>103,088</td><td>34,946</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#A</td><td>395,659</td><td>221,756</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#B</td><td>275,968</td><td>162,754</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#C</td><td>177,213</td><td>104,389</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>103,834</td><td>177,397</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#A</td><td>424,726</td><td>254,188</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#B</td><td>380,247</td><td>229,615</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#C</td><td>333,418</td><td>211,915</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>211,223</td><td>62,049</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#A</td><td>132,587</td><td>77,907</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#B</td><td>108,431</td><td>68,757</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#C</td><td>110,803</td><td>68,749</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>67,280</td><td>13,573</td></tr></tbody></table>

The number of cells obtained for each profile before and after cells merging was indicated. The profile #A+B+C represents the combined cytometric profiles obtained by the merging of profiles #A, #B, and #C. For each combined profile, the number of cells excluded during the merging was indicated. Samples with names starting with HEA correspond to healthy subjects, and samples with names starting with PAT correspond to HIV-infected patients.

have were merged. Thus, the processes of merging were thus efficient.

## HIV-Infection Induces Deep Phenotypic Modifications in Leukocytes

The SPADE algorithm was then used to identify 500 cell populations in the whole dataset of the six merged cytometric profiles. The 72 markers were used as clustering markers to capture the maximal cell diversity (Supplementary Figure 2A). For each cluster, the expression of CD3, CD11c, CD14, CD16, CD19, CD64, CD66, CD123, Granzyme B, and HLADR was represented (Supplementary Figure 2B).

We then generated a categorical heatmap representing the relative marker expression for each cluster to interpret the cell-cluster phenotypes (Figure 1). T cell, B cell, PMN, NK cell, monocyte, and dendritic cell clusters were directly annotated on the heatmap according to the expression of CD3, CD11c, CD14, CD16, CD19, CD66, CD123, Granzyme B, and HLADR (Figure 1A). Clusters annotated as unassigned in the heatmap correspond to cell clusters having <50 cells from the whole dataset. In these situations, the number of associated cells is not enough to properly define an MSI for each marker and to

associate a population to these cell clusters. The PMN population had the highest diversity and represented more than 50% of the cell clusters. In contrast, the pDC population represented only 0.4% of the cell clusters.

The 500 identified clusters were classified into two categories according to their cell abundance. The first category, named Enrichment Trend Clusters (ETCs), corresponded to clusters for which the cell abundance was specific to either HIV patients or healthy donors (cell abundance fold-change >4). The second category corresponded to other cell clusters (clusters without a specific enrichment trend). The trends of cell-cluster abundance are displayed at the bottom of the heatmap (Figure 1B). The sum of the associated cells from HIV-infected patients or healthy donors is also indicated for each cluster (Figure 1C).

The localization of the ETCs in the heatmap revealed groups of T cell, B cell, PMN, NK cell, and monocyte clusters that were homogeneously enriched either in HIV-infected patients or healthy donors. ETCs specific to the HIV condition were almost exclusively composed of cells from HIV-infected patients, whereas the number of cells from healthy subjects was negligible. Inversely, the enriched clusters for the healthy (HEA) condition were almost exclusively composed of cells from healthy donors,
<table><caption>**TABLE 3** | Number of events and cells obtained, before and after cell merging, for samples used to characterize specific cell populations from healthy subjects and HIV-infected patients.</caption><thead><tr><th>Sample ID</th><th>Stimulation</th><th>Profile</th><th>Number of events</th><th>Number of cells</th><th>Number of excluded cells</th></tr></thead><tbody><tr><td>HEA-1</td><td>CTRL</td><td>#A</td><td>443,527</td><td>237,039</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#B</td><td>321,545</td><td>192,201</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#C</td><td>245,550</td><td>139,218</td><td>-</td></tr><tr><td>HEA-1</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>139,104</td><td>151,146</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#A</td><td>464,263</td><td>249,773</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#B</td><td>463,584</td><td>249,003</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#C</td><td>413,188</td><td>229,388</td><td>-</td></tr><tr><td>HEA-2</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>229,058</td><td>40,990</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#A</td><td>304,750</td><td>128,604</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#B</td><td>243,011</td><td>112,246</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#C</td><td>219,468</td><td>103,360</td><td>-</td></tr><tr><td>HEA-3</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>103,088</td><td>34,946</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#A</td><td>395,659</td><td>221,756</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#B</td><td>275,968</td><td>162,754</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#C</td><td>177,213</td><td>104,389</td><td>-</td></tr><tr><td>PAT-1</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>103,834</td><td>177,397</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#A</td><td>424,726</td><td>254,188</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#B</td><td>380,247</td><td>229,615</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#C</td><td>333,418</td><td>211,915</td><td>-</td></tr><tr><td>PAT-2</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>211,223</td><td>62,049</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#A</td><td>132,587</td><td>77,907</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#B</td><td>108,431</td><td>68,757</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#C</td><td>110,803</td><td>68,749</td><td>-</td></tr><tr><td>PAT-3</td><td>CTRL</td><td>#A+B+C</td><td>-</td><td>67,280</td><td>13,573</td></tr></tbody></table>

The number of cells obtained for each profile before and after cells merging was indicated. The profile #A+B+C represents the combined cytometric profiles obtained by the merging of profiles #A, #B, and #C. For each combined profile, the number of cells excluded during the merging was indicated. Samples with names starting with HEA correspond to healthy subjects, and samples with names starting with PAT correspond to HIV-infected patients.

have were merged. Thus, the processes of merging were thus efficient.

## HIV-Infection Induces Deep Phenotypic Modifications in Leukocytes

The SPADE algorithm was then used to identify 500 cell populations in the whole dataset of the six merged cytometric profiles. The 72 markers were used as clustering markers to capture the maximal cell diversity (Supplementary Figure 2A). For each cluster, the expression of CD3, CD11c, CD14, CD16, CD19, CD64, CD66, CD123, Granzyme B, and HLADR was represented (Supplementary Figure 2B).

We then generated a categorical heatmap representing the relative marker expression for each cluster to interpret the cell-cluster phenotypes (Figure 1). T cell, B cell, PMN, NK cell, monocyte, and dendritic cell clusters were directly annotated on the heatmap according to the expression of CD3, CD11c, CD14, CD16, CD19, CD66, CD123, Granzyme B, and HLADR (Figure 1A). Clusters annotated as unassigned in the heatmap correspond to cell clusters having <50 cells from the whole dataset. In these situations, the number of associated cells is not enough to properly define an MSI for each marker and to

associate a population to these cell clusters. The PMN population had the highest diversity and represented more than 50% of the cell clusters. In contrast, the pDC population represented only 0.4% of the cell clusters.

The 500 identified clusters were classified into two categories according to their cell abundance. The first category, named Enrichment Trend Clusters (ETCs), corresponded to clusters for which the cell abundance was specific to either HIV patients or healthy donors (cell abundance fold-change >4). The second category corresponded to other cell clusters (clusters without a specific enrichment trend). The trends of cell-cluster abundance are displayed at the bottom of the heatmap (Figure 1B). The sum of the associated cells from HIV-infected patients or healthy donors is also indicated for each cluster (Figure 1C).

The localization of the ETCs in the heatmap revealed groups of T cell, B cell, PMN, NK cell, and monocyte clusters that were homogeneously enriched either in HIV-infected patients or healthy donors. ETCs specific to the HIV condition were almost exclusively composed of cells from HIV-infected patients, whereas the number of cells from healthy subjects was negligible. Inversely, the enriched clusters for the healthy (HEA) condition were almost exclusively composed of cells from healthy donors,
<figure><img src="image_1.png" /></figure>

whereas the number of cells from HIV-infected patients was
negligible. We identified two pDC clusters: one enriched in
the HEA condition, and another in the HIV condition. We

compared the phenotypes of CD4T cells, CD8T cells, B
cells, PMNs, NK cells, monocytes, and pDCs of HIV-infected
patients and healthy donors using this heatmap representing
<figure><img src="image_3.png" /></figure>

whereas the number of cells from HIV-infected patients was
negligible. We identified two pDC clusters: one enriched in
the HEA condition, and another in the HIV condition. We

compared the phenotypes of CD4T cells, CD8T cells, B
cells, PMNs, NK cells, monocytes, and pDCs of HIV-infected
patients and healthy donors using this heatmap representing
<figure><img src="image_2.png" /></figure>

**FIGURE 2** | Dot plots representing the phenotypes of mDCs and basophils from HIV-infected patients and healthy donors. Cell clusters corresponding to mDCs were isolated based on the same SPADE analysis presented in **Figure 1**. The same approach was used to isolate basophils. (A) Biplot representations showing the co-expression of HLADR and either CD11a, CD11b, CD83, CD86, or TLR2 by mDC populations. (B) Biplot representations showing the co-expression of CD123 and either CD11b, CD38, or CD62L by basophil populations. Blue lines were added to better visualize the difference in marker expression between healthy donors and HIV-infected patients. MSI represents the mean (MSI) of the markers in the X axes. "BB" in subscript indicates backbone markers.

the level of 72 cell markers. No specific groups of basophils or mDCs were associated with a specific individual group. The cell clusters of these populations were thus annotated on the SPADE tree and the cells directly analyzed using biplot representations (Figures 2A,B). The CytoBackBone algorithm enabled to observe with a high level of detail differences in the levels of several cell markers. The results and statistical analyses performed here are summarized in **Supplementary Table 1**.

Most CD8T cells from healthy subjects were CXCR4<sup>mid</sup>, CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD27<sup>low</sup>, CD38<sup>neg</sup>, CD83<sup>neg</sup>, Granzyme A<sup>neg</sup>, and Granzyme B<sup>neg/low</sup>. In contrast, CD8T cells from HIV-infected patients were CXCR4<sup>high</sup>, CD11a<sup>high</sup>, CD11b<sup>mid</sup>, CD27<sup>mid</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, Granzyme A<sup>high</sup>, and Granzyme B<sup>high</sup> or CXCR4<sup>high</sup>, CD11a<sup>mid</sup>, CD11b<sup>mid</sup>, CD27<sup>high</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, Granzyme A<sup>mid</sup>, and Granzyme B<sup>mid</sup>. Although low in abundance, we observed a CD32<sup>+</sup> subpopulation in HIV-infected patients. Additionally, CD8T cells from HIV-infected and healthy individuals were mainly CCR7<sup>neg</sup>, although a set was CCR7<sup>high</sup> in healthy subjects.

CD4T cells from healthy subjects were CD11a<sup>low</sup>, CD27<sup>low</sup>, CD38<sup>neg</sup>, CD83<sup>neg</sup>, and Granzyme B<sup>neg/low</sup>. In contrast, CD4T cells from HIV-infected patients were CD11a<sup>mid</sup>, CD27<sup>mid</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, and Granzyme B<sup>mid</sup>. Although low in abundance, we observed a CD32<sup>+</sup> subpopulation in HIV-infected patients.

B cells from HIV-infected patients had higher expression of CXCR4, HLADR, CD38, and CD83 than healthy subjects. Indeed, B cells from HIV-infected patients were CXCR4<sup>high</sup>, HLADR<sup>bright</sup>, CD38<sup>mid</sup>, and CD83<sup>high</sup>, whereas B cells from healthy individuals were CXCR4<sup>mid</sup>, HLADR<sup>mid/high</sup>, CD38<sup>neg</sup>, and CD83<sup>neg</sup>.

NK cells from healthy donors were CXCR4<sup>mid</sup>, CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD38<sup>neg</sup>, Granzyme A<sup>neg</sup>, Granzyme B<sup>neg</sup>, and Perforin<sup>neg</sup>, whereas they were CXCR4<sup>high</sup>, CD11a<sup>high</sup>, CD11b<sup>mid</sup>, CD38<sup>high</sup>, Granzyme A<sup>high</sup>, Granzyme B<sup>high</sup>, and Perforin<sup>high</sup> in HIV-infected patients.

PMNs from HIV-infected patients were CD11a<sup>high</sup>, CD11b<sup>high</sup>, CD32<sup>high</sup>, CD38<sup>high</sup>, CD64<sup>high</sup>, and TLR2<sup>mid</sup>.
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Dot plots representing the phenotypes of mDCs and basophils from HIV-infected patients and healthy donors. Cell clusters corresponding to mDCs were isolated based on the same SPADE analysis presented in **Figure 1**. The same approach was used to isolate basophils. (A) Biplot representations showing the co-expression of HLADR and either CD11a, CD11b, CD83, CD86, or TLR2 by mDC populations. (B) Biplot representations showing the co-expression of CD123 and either CD11b, CD38, or CD62L by basophil populations. Blue lines were added to better visualize the difference in marker expression between healthy donors and HIV-infected patients. MSI represents the mean (MSI) of the markers in the X axes. "BB" in subscript indicates backbone markers.

the level of 72 cell markers. No specific groups of basophils or mDCs were associated with a specific individual group. The cell clusters of these populations were thus annotated on the SPADE tree and the cells directly analyzed using biplot representations (Figures 2A,B). The CytoBackBone algorithm enabled to observe with a high level of detail differences in the levels of several cell markers. The results and statistical analyses performed here are summarized in **Supplementary Table 1**.

Most CD8T cells from healthy subjects were CXCR4<sup>mid</sup>, CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD27<sup>low</sup>, CD38<sup>neg</sup>, CD83<sup>neg</sup>, Granzyme A<sup>neg</sup>, and Granzyme B<sup>neg/low</sup>. In contrast, CD8T cells from HIV-infected patients were CXCR4<sup>high</sup>, CD11a<sup>high</sup>, CD11b<sup>mid</sup>, CD27<sup>mid</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, Granzyme A<sup>high</sup>, and Granzyme B<sup>high</sup> or CXCR4<sup>high</sup>, CD11a<sup>mid</sup>, CD11b<sup>mid</sup>, CD27<sup>high</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, Granzyme A<sup>mid</sup>, and Granzyme B<sup>mid</sup>. Although low in abundance, we observed a CD32<sup>+</sup> subpopulation in HIV-infected patients. Additionally, CD8T cells from HIV-infected and healthy individuals were mainly CCR7<sup>neg</sup>, although a set was CCR7<sup>high</sup> in healthy subjects.

CD4T cells from healthy subjects were CD11a<sup>low</sup>, CD27<sup>low</sup>, CD38<sup>neg</sup>, CD83<sup>neg</sup>, and Granzyme B<sup>neg/low</sup>. In contrast, CD4T cells from HIV-infected patients were CD11a<sup>mid</sup>, CD27<sup>mid</sup>, CD38<sup>mid</sup>, CD83<sup>mid</sup>, and Granzyme B<sup>mid</sup>. Although low in abundance, we observed a CD32<sup>+</sup> subpopulation in HIV-infected patients.

B cells from HIV-infected patients had higher expression of CXCR4, HLADR, CD38, and CD83 than healthy subjects. Indeed, B cells from HIV-infected patients were CXCR4<sup>high</sup>, HLADR<sup>bright</sup>, CD38<sup>mid</sup>, and CD83<sup>high</sup>, whereas B cells from healthy individuals were CXCR4<sup>mid</sup>, HLADR<sup>mid/high</sup>, CD38<sup>neg</sup>, and CD83<sup>neg</sup>.

NK cells from healthy donors were CXCR4<sup>mid</sup>, CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD38<sup>neg</sup>, Granzyme A<sup>neg</sup>, Granzyme B<sup>neg</sup>, and Perforin<sup>neg</sup>, whereas they were CXCR4<sup>high</sup>, CD11a<sup>high</sup>, CD11b<sup>mid</sup>, CD38<sup>high</sup>, Granzyme A<sup>high</sup>, Granzyme B<sup>high</sup>, and Perforin<sup>high</sup> in HIV-infected patients.

PMNs from HIV-infected patients were CD11a<sup>high</sup>, CD11b<sup>high</sup>, CD32<sup>high</sup>, CD38<sup>high</sup>, CD64<sup>high</sup>, and TLR2<sup>mid</sup>.
However, PMNs from healthy donors were CD11a<sup>low</sup>, CD11b<sup>low</sup>, CD32<sup>neg</sup>, CD38<sup>neg</sup>, CD64<sup>neg</sup>, and TLR2<sup>low</sup>. PMNs from healthy donors also exhibited a set of CCR7<sup>high</sup> clusters. Basophils from HIV patients were CD11b<sup>high</sup>, CD38<sup>mid</sup>, CD62L<sup>high</sup>, whereas they were CD11b<sup>neg</sup>, CD38<sup>neg</sup>, and CD62L<sup>neg</sup> in healthy donors.

Monocytes from healthy donors were CD11a<sup>mid</sup>, CD11b<sup>neg</sup>, CD11c<sup>mid</sup>, CD32<sup>neg</sup>, CD38<sup>neg</sup>, TLR2<sup>mid</sup>, and MCP1<sup>neg</sup>, whereas they were CD11a<sup>high</sup>, CD11b<sup>high</sup>, CD11c<sup>high</sup>, CD32<sup>high</sup>, CD38<sup>high</sup>, TLR2<sup>high</sup>, and MCP1<sup>mid/high</sup> in HIV-infected patients. Although CD64, CD83, CD86, and HLADR were highly expressed by both HIV-infected patients and healthy donors, their levels were higher in HIV-infected patients.

The mDC population can be split into two parts based on the expression of HLADR. However, the modifications of protein expression levels in HIV-infected patients were similar for the two mDC populations. mDCs from HIV-infected patients were CD11a<sup>mid</sup>, CD11b<sup>mid</sup>, CD83<sup>bright</sup>, CD86<sup>bright</sup>, TLR2<sup>high</sup>, and HLADR<sup>bright</sup>, whereas those from healthy donors were CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD83<sup>mid</sup>, CD86<sup>mid</sup>, TLR2<sup>mid</sup>, and HLADR<sup>high</sup>.

pDCs from HIV-infected patients were CD11b<sup>mid</sup> and Granzyme B<sup>mid</sup>, unlike pDCs from healthy donors which were CD11b<sup>neg</sup> and Granzyme B<sup>neg</sup>. Although HLADR was highly expressed in both HIV-infected patients and healthy donors, its level was higher in HIV-infected individuals.

Here, we have described HIV infection-dependent phenotypic changes in T cells, B cells, NK cells, PMNs, monocytes, dendritic cells, and basophils. Overall, the monocyte and PMN populations appeared to be the cell populations most highly affected by HIV infection.

## Confirmation of HIV Infection-Dependent Phenotypic Changes Observed in Monocytes

We next intended to confirm several phenotypic changes obtained with the combined cytometry profiles of 72 markers. To do so, we designed a single mass cytometry panel of 25 cell markers targeting the monocyte population, which was used to stain six new healthy donors (named HEA-4 to HEA-9), three previously used HIV-ART patients (named PAT-1 to PAT-3) and three new HIV-ART patients (named PAT-4 to PAT-6).

The monocyte population was used for validation based on two main reasons. First, monocytes were the most impacted population by HIV and were thus of strong relevance in term of disease characterization. Second, we previously predicted the upregulation of CD11b, CD32, CD64, CD83, CD86, and TLR2 in monocytes from HIV-infected patients. However, these markers were located on different panels. CD32, CD64, and CD86 were present in the backbone, whereas TLR2 was present in panel #A, CD11b in panel #B, and CD83 in panel #C (Table 2). Therefore, the monocyte population represented a good opportunity to validate our results, as these cells expressed these six markers. Thus, we created a single mass cytometry panel of 25 markers that contained these six specific markers (Supplementary Table 2).

Leukocytes from the original three HIV-infected patients (named PAT-1 to PAT-3), three new HIV-infected patients

(named PAT-4 to PAT-6), and six new healthy donors (named HEA-4 to HEA-9) (Table 1) were stained using this panel (Supplementary Table 2). After the acquisition of the whole dataset, SPADE analysis was performed. The SPADE analysis was parameterized to obtain 100 clusters using 5% down-sampling, and the clustering based on the levels of CD1c, CD3, CD4, CD8, CD11c, CD14, CD16, CD19, CD32, CD64, CD66, CD86, CD123, CD141, HLADR, and Granzyme B. Then, we annotated the SPADE cell clusters based on the levels of HLADR, CD3, CD11c, CD14, CD16, CD19, CD64, CD66, and CD123. Thereafter, monocytes were computationally isolated. Finally, we visualized the upregulation and co-expression of CD11b, CD32, CD64, CD83, CD86, and TLR2 on monocytes by performing a viSNE analysis (24).

The viSNE visualization, based on the expression of CD11b, CD32, CD64, CD83, CD86, and TLR2, allowed the perfect separation of monocytes from healthy and HIV-infected individuals (Figure 3A). This result demonstrates that the phenotype of monocytes from healthy donors was different from that from HIV-infected individuals. The use of statistical comparisons based on permutation tests allowed us to show that monocytes from HIV-infected patients significantly upregulated (MSI) CD11b ($p = 0.0011$), CD32 ($p \le 0.0001$), CD64 ($p = 0.0028$), CD83 ($p = 0.0016$), CD86 ($p = 0.0410$), and TLR2 ($p = 0.0008$) (Figure 3B). In addition, the viSNE visualization allowed us to show that most monocytes from HIV-infected patients were CD11b<sup>high</sup>, CD32<sup>high</sup>, CD64<sup>bright</sup>, CD83<sup>bright</sup>, CD86<sup>bright</sup>, and TLR2<sup>high</sup>. However, there was a subpopulation of monocytes with a CD64<sup>bright</sup> and TLR2<sup>neg</sup> phenotype.

Results from the validation experiment were compared to those obtained from the combined mass cytometry profiles using classical dot plots (Figures 4A,B). We observed the same trends of co-expression but no statistical analyses were performed, as these two sets of data were obtained using different panels and were unbalanced in terms of the number of samples.

## DISCUSSION

So far, multi-tube cytometry approaches have not been used in the context of chronic inflammation characterization. In addition to the technical difficulties to perform these experiments using mass cytometry panels, the analysis of such complex data requires further bioinformatics approaches. Here, we addressed these analysis challenges using the SPADE and SPADEVizR pipelines. Even if this study is based on a restricted number of samples from HIV-infected patients, the methodology that we describe here demonstrates that multi-tube cytometry can be applied for characterizing complex diseases.

Here, we characterized the HIV inflammation at an unprocessed phenotypical resolution of 72 cell marker. We detected the upregulation of several proteins in HIV-infected patients. Among them, CD11a and CD11b were upregulated in PMNs, monocytes, mDCs, NK cells, and T cells. CD11b was also upregulated on pDCs. Other upregulated proteins included: CD38 on PMNs, monocytes, NK cells, basophils, B cells, and T cells; CD83 on monocytes, mDCs, B cells, and T cells; and TLR2,
However, PMNs from healthy donors were CD11a<sup>low</sup>, CD11b<sup>low</sup>, CD32<sup>neg</sup>, CD38<sup>neg</sup>, CD64<sup>neg</sup>, and TLR2<sup>low</sup>. PMNs from healthy donors also exhibited a set of CCR7<sup>high</sup> clusters. Basophils from HIV patients were CD11b<sup>high</sup>, CD38<sup>mid</sup>, CD62L<sup>high</sup>, whereas they were CD11b<sup>neg</sup>, CD38<sup>neg</sup>, and CD62L<sup>neg</sup> in healthy donors.

Monocytes from healthy donors were CD11a<sup>mid</sup>, CD11b<sup>neg</sup>, CD11c<sup>mid</sup>, CD32<sup>neg</sup>, CD38<sup>neg</sup>, TLR2<sup>mid</sup>, and MCP1<sup>neg</sup>, whereas they were CD11a<sup>high</sup>, CD11b<sup>high</sup>, CD11c<sup>high</sup>, CD32<sup>high</sup>, CD38<sup>high</sup>, TLR2<sup>high</sup>, and MCP1<sup>mid/high</sup> in HIV-infected patients. Although CD64, CD83, CD86, and HLADR were highly expressed by both HIV-infected patients and healthy donors, their levels were higher in HIV-infected patients.

The mDC population can be split into two parts based on the expression of HLADR. However, the modifications of protein expression levels in HIV-infected patients were similar for the two mDC populations. mDCs from HIV-infected patients were CD11a<sup>mid</sup>, CD11b<sup>mid</sup>, CD83<sup>bright</sup>, CD86<sup>bright</sup>, TLR2<sup>high</sup>, and HLADR<sup>bright</sup>, whereas those from healthy donors were CD11a<sup>low</sup>, CD11b<sup>neg</sup>, CD83<sup>mid</sup>, CD86<sup>mid</sup>, TLR2<sup>mid</sup>, and HLADR<sup>high</sup>.

pDCs from HIV-infected patients were CD11b<sup>mid</sup> and Granzyme B<sup>mid</sup>, unlike pDCs from healthy donors which were CD11b<sup>neg</sup> and Granzyme B<sup>neg</sup>. Although HLADR was highly expressed in both HIV-infected patients and healthy donors, its level was higher in HIV-infected individuals.

Here, we have described HIV infection-dependent phenotypic changes in T cells, B cells, NK cells, PMNs, monocytes, dendritic cells, and basophils. Overall, the monocyte and PMN populations appeared to be the cell populations most highly affected by HIV infection.

## Confirmation of HIV Infection-Dependent Phenotypic Changes Observed in Monocytes

We next intended to confirm several phenotypic changes obtained with the combined cytometry profiles of 72 markers. To do so, we designed a single mass cytometry panel of 25 cell markers targeting the monocyte population, which was used to stain six new healthy donors (named HEA-4 to HEA-9), three previously used HIV-ART patients (named PAT-1 to PAT-3) and three new HIV-ART patients (named PAT-4 to PAT-6).

The monocyte population was used for validation based on two main reasons. First, monocytes were the most impacted population by HIV and were thus of strong relevance in term of disease characterization. Second, we previously predicted the upregulation of CD11b, CD32, CD64, CD83, CD86, and TLR2 in monocytes from HIV-infected patients. However, these markers were located on different panels. CD32, CD64, and CD86 were present in the backbone, whereas TLR2 was present in panel #A, CD11b in panel #B, and CD83 in panel #C (Table 2). Therefore, the monocyte population represented a good opportunity to validate our results, as these cells expressed these six markers. Thus, we created a single mass cytometry panel of 25 markers that contained these six specific markers (Supplementary Table 2).

Leukocytes from the original three HIV-infected patients (named PAT-1 to PAT-3), three new HIV-infected patients

(named PAT-4 to PAT-6), and six new healthy donors (named HEA-4 to HEA-9) (Table 1) were stained using this panel (Supplementary Table 2). After the acquisition of the whole dataset, SPADE analysis was performed. The SPADE analysis was parameterized to obtain 100 clusters using 5% down-sampling, and the clustering based on the levels of CD1c, CD3, CD4, CD8, CD11c, CD14, CD16, CD19, CD32, CD64, CD66, CD86, CD123, CD141, HLADR, and Granzyme B. Then, we annotated the SPADE cell clusters based on the levels of HLADR, CD3, CD11c, CD14, CD16, CD19, CD64, CD66, and CD123. Thereafter, monocytes were computationally isolated. Finally, we visualized the upregulation and co-expression of CD11b, CD32, CD64, CD83, CD86, and TLR2 on monocytes by performing a viSNE analysis (24).

The viSNE visualization, based on the expression of CD11b, CD32, CD64, CD83, CD86, and TLR2, allowed the perfect separation of monocytes from healthy and HIV-infected individuals (Figure 3A). This result demonstrates that the phenotype of monocytes from healthy donors was different from that from HIV-infected individuals. The use of statistical comparisons based on permutation tests allowed us to show that monocytes from HIV-infected patients significantly upregulated (MSI) CD11b ($p = 0.0011$), CD32 ($p \le 0.0001$), CD64 ($p = 0.0028$), CD83 ($p = 0.0016$), CD86 ($p = 0.0410$), and TLR2 ($p = 0.0008$) (Figure 3B). In addition, the viSNE visualization allowed us to show that most monocytes from HIV-infected patients were CD11b<sup>high</sup>, CD32<sup>high</sup>, CD64<sup>bright</sup>, CD83<sup>bright</sup>, CD86<sup>bright</sup>, and TLR2<sup>high</sup>. However, there was a subpopulation of monocytes with a CD64<sup>bright</sup> and TLR2<sup>neg</sup> phenotype.

Results from the validation experiment were compared to those obtained from the combined mass cytometry profiles using classical dot plots (Figures 4A,B). We observed the same trends of co-expression but no statistical analyses were performed, as these two sets of data were obtained using different panels and were unbalanced in terms of the number of samples.

## DISCUSSION

So far, multi-tube cytometry approaches have not been used in the context of chronic inflammation characterization. In addition to the technical difficulties to perform these experiments using mass cytometry panels, the analysis of such complex data requires further bioinformatics approaches. Here, we addressed these analysis challenges using the SPADE and SPADEVizR pipelines. Even if this study is based on a restricted number of samples from HIV-infected patients, the methodology that we describe here demonstrates that multi-tube cytometry can be applied for characterizing complex diseases.

Here, we characterized the HIV inflammation at an unprocessed phenotypical resolution of 72 cell marker. We detected the upregulation of several proteins in HIV-infected patients. Among them, CD11a and CD11b were upregulated in PMNs, monocytes, mDCs, NK cells, and T cells. CD11b was also upregulated on pDCs. Other upregulated proteins included: CD38 on PMNs, monocytes, NK cells, basophils, B cells, and T cells; CD83 on monocytes, mDCs, B cells, and T cells; and TLR2,
<figure><img src="image_3.png" /><figcaption>**FIGURE 3** | viSNE map representing the phenotype of leukocytes from HIV-infected patients and healthy donors using a single antibody panel. Leukocytes from the three original HIV-infected subjects, three new HIV-infected patients, and six new healthy subjects were collected. SPADE analysis was performed to isolate monocytes. Monocyte cells were computationally distributed into a two-dimensional space using the viSNE algorithm. Such viSNE representations were generated based on the expressions of CD11b, CD32, CD64, CD83, CD86, and TLR2. (A) The viSNE map is colored to represent cells from HIV-infected patients in orange and those from healthy individuals in purple. (B) The viSNE map is colored to represent the mean marker expression (MSI) of CD11b, CD32, CD64, CD83, CD86, and TLR2 for each cell. Comparisons of the MSIs from the individual groups were performed using non-parametric permutation tests.</figcaption></figure>

CD32, and CD64 on PMNs and monocytes. The monocyte and PMN population appeared to be strongly affected by the HIV infection, as CD11a, CD11b, CD32, CD38, CD64, CD83, CD86, and TLR2 were upregulated in these populations. Overall, these results show that HIV infection induced a specific environment that similarly affected various immune cells in these patients.

Following the detection of Pathogen Associated Molecular Patterns (PAMP) in tissues by cell populations expressing pattern recognition receptors (PRRs), different cytokines and chemokines are produced. The release of these inflammatory proteins leads to the recruitment of innate immune cells, such as neutrophils, monocytes, dendritic cells, and NK cells

(25, 26). Whereas, IL-8 and MIP-1β upregulate expression level of adhesion molecules on innate immune cells, cytokines such as IL-1α, IL-6, and TNF-α upregulate the expression levels of selectin molecules on endothelial cells. Interactions between adhesion molecules and selectin molecules trigger the diapedesis mechanism (27–29). In our study, we found that the HIV infection induced an overexpression of adhesion molecules CD11a or CD11b on the surface of neutrophils, monocytes, and dendritic cells. An overexpression of CD11a in NK cells was also detected. Thus, these results suggest that immune cells could be more easily recruited on inflammatory sites.
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | viSNE map representing the phenotype of leukocytes from HIV-infected patients and healthy donors using a single antibody panel. Leukocytes from the three original HIV-infected subjects, three new HIV-infected patients, and six new healthy subjects were collected. SPADE analysis was performed to isolate monocytes. Monocyte cells were computationally distributed into a two-dimensional space using the viSNE algorithm. Such viSNE representations were generated based on the expressions of CD11b, CD32, CD64, CD83, CD86, and TLR2. (A) The viSNE map is colored to represent cells from HIV-infected patients in orange and those from healthy individuals in purple. (B) The viSNE map is colored to represent the mean marker expression (MSI) of CD11b, CD32, CD64, CD83, CD86, and TLR2 for each cell. Comparisons of the MSIs from the individual groups were performed using non-parametric permutation tests.</figcaption></figure>

CD32, and CD64 on PMNs and monocytes. The monocyte and PMN population appeared to be strongly affected by the HIV infection, as CD11a, CD11b, CD32, CD38, CD64, CD83, CD86, and TLR2 were upregulated in these populations. Overall, these results show that HIV infection induced a specific environment that similarly affected various immune cells in these patients.

Following the detection of Pathogen Associated Molecular Patterns (PAMP) in tissues by cell populations expressing pattern recognition receptors (PRRs), different cytokines and chemokines are produced. The release of these inflammatory proteins leads to the recruitment of innate immune cells, such as neutrophils, monocytes, dendritic cells, and NK cells

(25, 26). Whereas, IL-8 and MIP-1β upregulate expression level of adhesion molecules on innate immune cells, cytokines such as IL-1α, IL-6, and TNF-α upregulate the expression levels of selectin molecules on endothelial cells. Interactions between adhesion molecules and selectin molecules trigger the diapedesis mechanism (27–29). In our study, we found that the HIV infection induced an overexpression of adhesion molecules CD11a or CD11b on the surface of neutrophils, monocytes, and dendritic cells. An overexpression of CD11a in NK cells was also detected. Thus, these results suggest that immune cells could be more easily recruited on inflammatory sites.
<figure><img src="image_4.png" /><figcaption>**FIGURE 4** | Confirmation of results generated by CytoBackBone. (A) Leukocytes were collected from HIV-infected patients and healthy subjects. Cells were stained with a single antibody panel and SPADE analysis was performed to isolate monocytes. The biplot representations show the co-expression of CD64 and TLR2, CD11b and CD32, and CD86 and CD83 in monocytes from patient PAT-1. (B) Leukocytes were collected from HIV-infected patients and healthy subjects. Cells were stained with three different antibody panels. Then, CytoBacBone was used to merge the different cytometry profiles. SPADE analysis was then performed to isolate monocytes. The biplot representations show the co-expression of CD64 and TLR2, CD11b and CD32, and CD86 and CD83 in monocytes from patient PAT-1.</figcaption></figure>

The HIV infection induces an exhaustion of different innate immune cells, such as neutrophils, monocytes, and dendritic cells. Indeed, we demonstrated that the production of IL-1α, IL-6, and TNF-α were delayed following different TLR stimulations in HIV-infected patients (15). These proteins are essential for the induction of selectin molecules on the surface of endothelial cells. Without these cytokines, leukocytes expressing high levels

of CD11a and CD11b cannot adhere to endothelial cells (27-29), which greatly limits the recruitment of immune cells to inflammatory sites.

Innate immunity is strongly activated in primary infection as depicted by the cytokine storm present in the first weeks of HIV infection (30). Numerous types of immune cells are involved in the acute inflammation seen in primary infection.
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Confirmation of results generated by CytoBackBone. (A) Leukocytes were collected from HIV-infected patients and healthy subjects. Cells were stained with a single antibody panel and SPADE analysis was performed to isolate monocytes. The biplot representations show the co-expression of CD64 and TLR2, CD11b and CD32, and CD86 and CD83 in monocytes from patient PAT-1. (B) Leukocytes were collected from HIV-infected patients and healthy subjects. Cells were stained with three different antibody panels. Then, CytoBacBone was used to merge the different cytometry profiles. SPADE analysis was then performed to isolate monocytes. The biplot representations show the co-expression of CD64 and TLR2, CD11b and CD32, and CD86 and CD83 in monocytes from patient PAT-1.</figcaption></figure>

The HIV infection induces an exhaustion of different innate immune cells, such as neutrophils, monocytes, and dendritic cells. Indeed, we demonstrated that the production of IL-1α, IL-6, and TNF-α were delayed following different TLR stimulations in HIV-infected patients (15). These proteins are essential for the induction of selectin molecules on the surface of endothelial cells. Without these cytokines, leukocytes expressing high levels

of CD11a and CD11b cannot adhere to endothelial cells (27-29), which greatly limits the recruitment of immune cells to inflammatory sites.

Innate immunity is strongly activated in primary infection as depicted by the cytokine storm present in the first weeks of HIV infection (30). Numerous types of immune cells are involved in the acute inflammation seen in primary infection.
Additionally, these cells are involved in the persistence of the chronic inflammation observed in treated HIV-infected patients. Moreover, a low neutrophil count is associated with the disease progression, a CD4⁺ T cells loss, and a high HIV viral load (31).

Neutrophils contribute to an anti-HIV response through the production of α-defensins, myeloperoxidase, and reactive oxygen species. Neutrophils also participate in the elimination of HIV-infected cells using antibody-dependent cell-mediated cytotoxic processes (ADCC) but are less effective than monocytes and NK cells. Finally, neutrophils and monocytes can perform antibody-dependent cellular phagocytosis of infected cells and immune complexes (32). Thus, phenotype impairment related to PMN and monocytes could be severely harmful. In our study, we observed strong upregulations of CD32 and CD64 on PMN and monocyte populations. CD32 and CD64 are two Fc receptors playing a major role in antibody-mediated processes such as phagocytosis. As phagocytosis is known to be affected by the HIV infection (33, 34), upregulation of these proteins could compensate for this impairment.

CD32 and CD64 also mediate ADCC, which is important in the response against HIV. ADCC is mostly performed by NK cells, but we did not observe an upregulation of Fc receptors on this population. However, monocytes and PMNs also carry out ADCC (35–37). The upregulation of CD32 and CD64 could suggest an increase of ADCC by these cells. Additionally, the increase of expression level of Fc receptors could facilitate the development and maintenance of chronic inflammation. Indeed, the engagement of these receptors leads to the production of cytokines, chemokines, and inflammatory mediators. Consequently, due to the increased expression of these Fc-receptors, associated signaling pathways are probably more expressed. Therefore, this increased expression would result in a greater release of inflammatory molecules, which would promote the maintenance of chronic inflammation.

We also detected an overexpression of MCP-1 on monocytes. MCP-1 is a chemokine that attracts monocytes, lymphocytes, and polymorphonuclear cells. The HIV-NEF protein is known to upregulate the expression of MCP-1 in astrocytes, the most abundant cell type in the brain (38). The detection of MCP-1 on monocytes can be explained by the fact that MCP-1 binds to the CCR2 receptor expressed on monocytes. This could lead to the recruitment of monocytes to the brain through the blood-brain barrier, which could subsequently enhance both the generation of the viral reservoir and neurological complications (39).

Alterations in cytokine expression profiles have reported in treated and non-treated HIV-infected individuals. Here, several cytokines were included in our antibody panels, such as IL-1α, IL-1β, IL-1RA, IL2, IL4, IL6, IL8, IL10, IL12, IP-10, IFN-α, IFN-γ, MIP-1β, TNF-α, and MCP-1. Except for MCP-1, no signal of these cytokines was detected in the different leukocyte populations. This no detection of cytokines is explained by the fact that leukocytes were not stimulated by ligands such as LPS or R848. Without such stimulation, the expression level of cytokines inside leukocytes is too low to be detected.

In our study, no non-classical monocytes (CD14<sup>low</sup> and CD16<sup>high</sup>) were found in treated HIV-infected patients, neither

before nor after the cell merging process. Indeed, all monocytes identified in this study were CD14<sup>high</sup> and CD16<sup>neg/low</sup>. However, we found that all monocytes up-regulated the expression of CD11a, CD11b, CD11c, CD32, CD38, TLR2, and MCP-1. Thus, these results suggested that both classical and non-classical monocytes were impacted.

Classical myeloid dendritic cells from HIV-infected patients exhibit altered functions including a weak capacity of maturation and inefficient antigen presentation to CD4⁺ T cells (40). Here, we observed a strong up-regulation of CD86 on monocyte and mDC populations, which is a protein essential for antigen presentation. The upregulation of this molecule on mDCs could compensate the dysregulations related to antigen presentation.

Furthermore, we found that the expression of CD83 was higher on monocytes and mDCs from HIV-infected patients compared to healthy subjects. CD83 is a costimulatory molecule and is upregulated after cell activation (41) and maturation (42). CD83-stimulated monocytes suppress T cell immune responses through the production of prostaglandin E2 (43). Upregulation of CD83 also highlights chronic activation of the immune system in ART-treated HIV-infected patients (42).

As pDCs are among the main cells producing IFN-I, they play a main role in the antiviral responses. In HIV-infected patients, the number of pDCs is known to be decreased compared to healthy patients (44). Moreover, the production of cytokines is also impacted by the infection, even in treated HIV-infected subjects (15). In our study, pDCs displayed a hyperactivation phenotype through an up-regulation of CD11b. Thus, we can suppose that pDCs chronically produce IFN-α, which can participate in the increase of myeloid-derived suppressor cells (MDSCs) counts.

Higher transcription of genes encoding TLR2, TLR3, TLR4, TLR6, TLR7, and TLR8 has been described in PBMCs from HIV infected individuals (45, 46). Here, we observed an overexpression of TLR2 in monocytes and PMNs from treated HIV-infected patients. As HIV is able to inhibit signaling pathways induced following the TLRs engagement, the upregulation of these receptors could compensate these inhibitions. In addition, the overexpression of TLRs could have a deleterious effect on the spread of HIV. Indeed, HIV is able to use TLR-dependent signaling pathways to facilitate its replication and propagation (47, 48). Importantly, as the microbial translocation is favored in HIV-infected patients, the multiplication of TLR2 activation could participate in the disease progression. Indeed, commensal bacteria (such as *Lactobacillus acidophilus*, *Prevotella melaninogenica*, *Prevotella bivia*, and *Mycobacterium smegmatis*) enhance the spread of HIV by activating TLR2-dependent signaling (49). Finally, the TLR2 activation is also associated with reactivation of latent viruses (50).

NK cells are critical antiviral effectors of the innate immune system, and deficiencies of this cell population are associated with an increased probability of HIV infection (51, 52). NK cells contribute to the elimination of HIV-infected cells during the acute phase. Indeed, NK cells can release cytotoxic granules, cytokines, and chemokines through the activation of the KIR receptors, the natural cytotoxicity receptors, the C-type lectin
Additionally, these cells are involved in the persistence of the chronic inflammation observed in treated HIV-infected patients. Moreover, a low neutrophil count is associated with the disease progression, a CD4+ T cells loss, and a high HIV viral load (31).

Neutrophils contribute to an anti-HIV response through the production of α-defensins, myeloperoxidase, and reactive oxygen species. Neutrophils also participate in the elimination of HIV-infected cells using antibody-dependent cell-mediated cytotoxic processes (ADCC) but are less effective than monocytes and NK cells. Finally, neutrophils and monocytes can perform antibody-dependent cellular phagocytosis of infected cells and immune complexes (32). Thus, phenotype impairment related to PMN and monocytes could be severely harmful. In our study, we observed strong upregulations of CD32 and CD64 on PMN and monocyte populations. CD32 and CD64 are two Fc receptors playing a major role in antibody-mediated processes such as phagocytosis. As phagocytosis is known to be affected by the HIV infection (33, 34), upregulation of these proteins could compensate for this impairment.

CD32 and CD64 also mediate ADCC, which is important in the response against HIV. ADCC is mostly performed by NK cells, but we did not observe an upregulation of Fc receptors on this population. However, monocytes and PMNs also carry out ADCC (35-37). The upregulation of CD32 and CD64 could suggest an increase of ADCC by these cells. Additionally, the increase of expression level of Fc receptors could facilitate the development and maintenance of chronic inflammation. Indeed, the engagement of these receptors leads to the production of cytokines, chemokines, and inflammatory mediators. Consequently, due to the increased expression of these Fc-receptors, associated signaling pathways are probably more expressed. Therefore, this increased expression would result in a greater release of inflammatory molecules, which would promote the maintenance of chronic inflammation.

We also detected an overexpression of MCP-1 on monocytes. MCP-1 is a chemokine that attracts monocytes, lymphocytes, and polymorphonuclear cells. The HIV-NEF protein is known to upregulate the expression of MCP-1 in astrocytes, the most abundant cell type in the brain (38). The detection of MCP-1 on monocytes can be explained by the fact that MCP-1 binds to the CCR2 receptor expressed on monocytes. This could lead to the recruitment of monocytes to the brain through the blood-brain barrier, which could subsequently enhance both the generation of the viral reservoir and neurological complications (39).

Alterations in cytokine expression profiles have reported in treated and non-treated HIV-infected individuals. Here, several cytokines were included in our antibody panels, such as IL-1α, IL-1β, IL-1RA, IL2, IL4, IL6, IL8, IL10, IL12, IP-10, IFN-α, IFN-γ, MIP-1β, TNF-α, and MCP-1. Except for MCP-1, no signal of these cytokines was detected in the different leukocyte populations. This no detection of cytokines is explained by the fact that leukocytes were not stimulated by ligands such as LPS or R848. Without such stimulation, the expression level of cytokines inside leukocytes is too low to be detected.

In our study, no non-classical monocytes (CD14^low^ and CD16^high^) were found in treated HIV-infected patients, neither

before nor after the cell merging process. Indeed, all monocytes identified in this study were CD14^high^ and CD16^neg/low^. However, we found that all monocytes up-regulated the expression of CD11a, CD11b, CD11c, CD32, CD38, TLR2, and MCP-1. Thus, these results suggested that both classical and non-classical monocytes were impacted.

Classical myeloid dendritic cells from HIV-infected patients exhibit altered functions including a weak capacity of maturation and inefficient antigen presentation to CD4+ T cells (40). Here, we observed a strong up-regulation of CD86 on monocyte and mDC populations, which is a protein essential for antigen presentation. The upregulation of this molecule on mDCs could compensate the dysregulations related to antigen presentation.

Furthermore, we found that the expression of CD83 was higher on monocytes and mDCs from HIV-infected patients compared to healthy subjects. CD83 is a costimulatory molecule and is upregulated after cell activation (41) and maturation (42). CD83-stimulated monocytes suppress T cell immune responses through the production of prostaglandin E2 (43). Upregulation of CD83 also highlights chronic activation of the immune system in ART-treated HIV-infected patients (42).

As pDCs are among the main cells producing IFN-I, they play a main role in the antiviral responses. In HIV-infected patients, the number of pDCs is known to be decreased compared to healthy patients (44). Moreover, the production of cytokines is also impacted by the infection, even in treated HIV-infected subjects (15). In our study, pDCs displayed a hyperactivation phenotype through an up-regulation of CD11b. Thus, we can suppose that pDCs chronically produce IFN-α, which can participate in the increase of myeloid-derived suppressor cells (MDSCs) counts.

Higher transcription of genes encoding TLR2, TLR3, TLR4, TLR6, TLR7, and TLR8 has been described in PBMCs from HIV infected individuals (45, 46). Here, we observed an overexpression of TLR2 in monocytes and PMNs from treated HIV-infected patients. As HIV is able to inhibit signaling pathways induced following the TLRs engagement, the upregulation of these receptors could compensate these inhibitions. In addition, the overexpression of TLRs could have a deleterious effect on the spread of HIV. Indeed, HIV is able to use TLR-dependent signaling pathways to facilitate its replication and propagation (47, 48). Importantly, as the microbial translocation is favored in HIV-infected patients, the multiplication of TLR2 activation could participate in the disease progression. Indeed, commensal bacteria (such as *Lactobacillus acidophilus*, *Prevotella melaninogenica*, *Prevotella bivia*, and *Mycobacterium smegmatis*) enhance the spread of HIV by activating TLR2-dependent signaling (49). Finally, the TLR2 activation is also associated with reactivation of latent viruses (50).

NK cells are critical antiviral effectors of the innate immune system, and deficiencies of this cell population are associated with an increased probability of HIV infection (51, 52). NK cells contribute to the elimination of HIV-infected cells during the acute phase. Indeed, NK cells can release cytotoxic granules, cytokines, and chemokines through the activation of the KIR receptors, the natural cytotoxicity receptors, the C-type lectin
receptors, signaling lymphocyte activation family receptors, and Fc receptors.

Although the ADCC is mostly performed by NK cells (51, 52), monocytes and PMNs are also able to do it. Here, we did not observe upregulations of CD32 and CD64 in NK cells from HIV-infected individuals, unlike in monocytes and PMNs. Thus, the capacity of monocytes and PMNs to perform ADCC could be increased in HIV-infected patients.

Here, alterations observed in NK cells from HIV-infected patients characterize a hyperactivation state, which leads to the persistence of chronic inflammation. This inflammation leads to alterations of the NK cells distribution and their functional capacities. Moreover, the chronic inflammation enhances the generation of anergic NK cells, reduces the ADCC activity, and contributes to a poor immunologic reconstitution of CD4 T cells in HIV-infected individuals (53). Finally, we showed that NK cells overexpressed CD38, which is related to progression to AIDS (54).

CD27 is a protein constitutively expressed by naive T cells. This protein plays a main role in proliferation, survival, and differentiation of T cells. Especially, CD27 promotes immune activation and enhances primary, secondary, memory, and recall responses toward viral infections (55). Following the activation of CD4 or CD8 T cells via the TCR/CD3 pathway, the expression of CD27 is increased in T cells. However, after prolonged activation, CD27 becomes gradually switched off (56). Therefore, a high level of CD27 on T cells is often considered as a marker of early activation. Herein, we demonstrated that the expressions of CD27 on CD4 and CD8 T cells were higher in HIV-1 donors compared to healthy donors. These results again describe a persistence of inflammation in HIV-infected patients, which could be associated either to the residual replication HIV or to HIV-dependent inflammatory mechanisms (such as higher microbial translocation). CD27 also plays a key role in the generation and long-term maintenance of T cell responses. The long-term expression of CD27 is associated with impaired maturation of T cells (55). Finally, as the constitutive engagement of CD27/CD70 promotes T cell exhaustion, high levels of CD27 observed in HIV patients could contribute to the loss of T cells effector functions (57).

We found that the expression of CD83 on T cells and B cells from HIV-infected patients was also upregulated compared to healthy donors. The expression of CD83 on CD4⁺/CD25⁺ T cells confers immunosuppressive functions (58), whereas the upregulation of CD83 on murine B cells has a regulatory role in humoral responses (59). Finally, an up-regulation of CD38 was observed on T cells and B cells from HIV-ART patients relative to control individuals, which is associated with HIV disease progression (54). Together, these results support that chronic inflammation could induce important dysfunctions in adaptive immune responses.

NK cells and CD8 T cells naturally express high levels of Granzymes or Perforin. High expressions of Granzyme B and Perforin in CD8 T cells and NK cells are a signature of immune activation, which could be associated here to the persistence of HIV. Indeed, HIV replication is known to induce the production of these cytotoxic proteins in these two cell

populations. Importantly, an increase in the number of T cells secreting Granzyme B is associated with reduced viral reservoirs in HIV infection. Thus, HIV infection could limit the secretion of Granzymes and Perforin, and could lead to an accumulation of granules in NK or CD8 T cells. Therefore, the accumulation of granules could lead to the loss of cytolytic activity for NK and CD8 T cells. Finally, the high level of granzymes or perforin observed in HIV-infected patients could be associated with Cytomegalovirus (CMV). Indeed, nearly 90% of HIV-infected patients are CMV positive contrasting with about 60% in healthy donors.

Because HIV-1 infected patients were ART-treated, we suppose that the replication of HIV is not involved in the reported alterations. However, the virus could still play a role, as increasing pieces of evidence point to the persistence of HIV in tissues despite undetectable viral load in plasma (60, 61). Here, we suppose that the phenotypical modifications observed in treated HIV-infected patients are linked to the persistence of inflammation. Indeed, chronic inflammation is associated with a hyperactivation and exhaustion of innate immune cells. This inflammation could be induced by a residual replication of HIV, the presence of HIV in reservoir cells, and an increasing microbial translocation.

Although our combined cytometry panel is composed of 72 cells markers, several important markers for phenotypical annotation of cell subsets are unfortunately absent. Indeed, markers such as CD10, CXCL13, CD278 (follicular T cells), NKG2A, NKp44 (NK cells), T-bet, GATA, RORγt, IL-17 (ILC), or FoxP3 (T regs) were not included in our antibodies panels. Thus, further research must be conducted to extend this approach to characterize all cutting-edge subpopulations involved in HIV.

In the viSNE analysis, monocytes were partly CD64<sup>bright</sup> and TLR2<sup>neg</sup> in the single-tube experiment. However, we did not observe this result in the heatmap from the multi-tube 72-marker experiment. The multi-tube 72-marker experiment was performed with three HIV-ART patients (named PAT-1 to PAT-3), whereas the single-tube experiment was performed with six HIV-ART patients (including PAT-1 to PAT-3, but also PAT-4 to PAT-6). After deeper analysis, we found this CD64<sup>bright</sup> and TLR2<sup>neg</sup> monocyte population to be detectable in patients PAT-2, PAT-4, PAT-5, and PAT-6. Thus, it was normal to not observe this population on the 72 cell-marker heatmap, as monocytes from patients PAT-1 to PAT-3 were mostly CD64<sup>bright</sup> and TLR2<sup>high</sup>.

Overall, these analyses lead to the same conclusions: the effect of HIV infection on most innate and adaptive immune cell types are still largely present in the six HIV-infected patients, with a prolonged undetectable viral load and, for most, CD4 T cell counts >500/mm³. These results are in agreement with several publications showing persistent chronic immune activation in patients on prolonged ART for whom increased levels of soluble biomarkers (such as sCD14 and IP10) were detected (62, 63).

## DATA AVAILABILITY

Raw and merged cytometric profiles used to characterize the cell phenotypes from healthy subjects and HIV-infected patients as
receptors, signaling lymphocyte activation family receptors, and Fc receptors.

Although the ADCC is mostly performed by NK cells (51, 52), monocytes and PMNs are also able to do it. Here, we did not observe upregulations of CD32 and CD64 in NK cells from HIV-infected individuals, unlike in monocytes and PMNs. Thus, the capacity of monocytes and PMNs to perform ADCC could be increased in HIV-infected patients.

Here, alterations observed in NK cells from HIV-infected patients characterize a hyperactivation state, which leads to the persistence of chronic inflammation. This inflammation leads to alterations of the NK cells distribution and their functional capacities. Moreover, the chronic inflammation enhances the generation of anergic NK cells, reduces the ADCC activity, and contributes to a poor immunologic reconstitution of CD4 T cells in HIV-infected individuals (53). Finally, we showed that NK cells overexpressed CD38, which is related to progression to AIDS (54).

CD27 is a protein constitutively expressed by naive T cells. This protein plays a main role in proliferation, survival, and differentiation of T cells. Especially, CD27 promotes immune activation and enhances primary, secondary, memory, and recall responses toward viral infections (55). Following the activation of CD4 or CD8 T cells via the TCR/CD3 pathway, the expression of CD27 is increased in T cells. However, after prolonged activation, CD27 becomes gradually switched off (56). Therefore, a high level of CD27 on T cells is often considered as a marker of early activation. Herein, we demonstrated that the expressions of CD27 on CD4 and CD8 T cells were higher in HIV-1 donors compared to healthy donors. These results again describe a persistence of inflammation in HIV-infected patients, which could be associated either to the residual replication HIV or to HIV-dependent inflammatory mechanisms (such as higher microbial translocation). CD27 also plays a key role in the generation and long-term maintenance of T cell responses. The long-term expression of CD27 is associated with impaired maturation of T cells (55). Finally, as the constitutive engagement of CD27/CD70 promotes T cell exhaustion, high levels of CD27 observed in HIV patients could contribute to the loss of T cells effector functions (57).

We found that the expression of CD83 on T cells and B cells from HIV-infected patients was also upregulated compared to healthy donors. The expression of CD83 on CD4⁺/CD25⁺ T cells confers immunosuppressive functions (58), whereas the upregulation of CD83 on murine B cells has a regulatory role in humoral responses (59). Finally, an up-regulation of CD38 was observed on T cells and B cells from HIV-ART patients relative to control individuals, which is associated with HIV disease progression (54). Together, these results support that chronic inflammation could induce important dysfunctions in adaptive immune responses.

NK cells and CD8 T cells naturally express high levels of Granzymes or Perforin. High expressions of Granzyme B and Perforin in CD8 T cells and NK cells are a signature of immune activation, which could be associated here to the persistence of HIV. Indeed, HIV replication is known to induce the production of these cytotoxic proteins in these two cell

populations. Importantly, an increase in the number of T cells secreting Granzyme B is associated with reduced viral reservoirs in HIV infection. Thus, HIV infection could limit the secretion of Granzymes and Perforin, and could lead to an accumulation of granules in NK or CD8 T cells. Therefore, the accumulation of granules could lead to the loss of cytolytic activity for NK and CD8 T cells. Finally, the high level of granzymes or perforin observed in HIV-infected patients could be associated with Cytomegalovirus (CMV). Indeed, nearly 90% of HIV-infected patients are CMV positive contrasting with about 60% in healthy donors.

Because HIV-1 infected patients were ART-treated, we suppose that the replication of HIV is not involved in the reported alterations. However, the virus could still play a role, as increasing pieces of evidence point to the persistence of HIV in tissues despite undetectable viral load in plasma (60, 61). Here, we suppose that the phenotypical modifications observed in treated HIV-infected patients are linked to the persistence of inflammation. Indeed, chronic inflammation is associated with a hyperactivation and exhaustion of innate immune cells. This inflammation could be induced by a residual replication of HIV, the presence of HIV in reservoir cells, and an increasing microbial translocation.

Although our combined cytometry panel is composed of 72 cells markers, several important markers for phenotypical annotation of cell subsets are unfortunately absent. Indeed, markers such as CD10, CXCL13, CD278 (follicular T cells), NKG2A, NKp44 (NK cells), T-bet, GATA, RORγt, IL-17 (ILC), or FoxP3 (T regs) were not included in our antibodies panels. Thus, further research must be conducted to extend this approach to characterize all cutting-edge subpopulations involved in HIV.

In the viSNE analysis, monocytes were partly CD64<sup>bright</sup> and TLR2<sup>neg</sup> in the single-tube experiment. However, we did not observe this result in the heatmap from the multi-tube 72-marker experiment. The multi-tube 72-marker experiment was performed with three HIV-ART patients (named PAT-1 to PAT-3), whereas the single-tube experiment was performed with six HIV-ART patients (including PAT-1 to PAT-3, but also PAT-4 to PAT-6). After deeper analysis, we found this CD64<sup>bright</sup> and TLR2<sup>neg</sup> monocyte population to be detectable in patients PAT-2, PAT-4, PAT-5, and PAT-6. Thus, it was normal to not observe this population on the 72 cell-marker heatmap, as monocytes from patients PAT-1 to PAT-3 were mostly CD64<sup>bright</sup> and TLR2<sup>high</sup>.

Overall, these analyses lead to the same conclusions: the effect of HIV infection on most innate and adaptive immune cell types are still largely present in the six HIV-infected patients, with a prolonged undetectable viral load and, for most, CD4 T cell counts >500/mm³. These results are in agreement with several publications showing persistent chronic immune activation in patients on prolonged ART for whom increased levels of soluble biomarkers (such as sCD14 and IP10) were detected (62, 63).

DATA AVAILABILITY

Raw and merged cytometric profiles used to characterize the cell phenotypes from healthy subjects and HIV-infected patients as
well as the cytometric profiles obtained from the confirmatory
single panel are available in the FlowRepository database through
the accession IDs: FR-FCM-Z26Z and FR-FCM-Z26Y.

ETHICS STATEMENT

This experiment was approved by the Comité de Protection des
Personnes (CPP) Ile de France VII, under protocol number PP
14-003. All subjects gave written informed consent to participate
in this study.

AUTHOR CONTRIBUTIONS

AL, NT, and AC: conceptualization and methodology, investigation. AL, NT, AC, RL, and OL: validation and writing—review and editing. AL and NT: formal analysis and writing—original draft. AL and OL: resources. AC and RL: funding acquisition. AC, RL, and OL: supervision.

REFERENCES

1. Bendall SC, Nolan GP, Roederer M, Chattopadhyay PK. A deep profiler's guide to cytometry. *Trends Immunol.* (2012) 33:323-32. doi: 10.1016/j.it.2012.02.010

2. O'Neill K, Aghaepour N, Parker J, Hogge D, Karsan A, Dalal B, Brinkman RR. Deep profiling of multitube flow cytometry data. *Bioinformatics*. (2015) 31:1623-31. doi: 10.1093/bioinformatics/btv008

3. Pedreira CE, Costa ES, Barrena S, Lecrevisse Q, Almeida J, Van Dongen JJM, et al. Generation of flow cytometry data files with a potentially infinite number of dimensions. *Cytometry Part A*. (2008) 73:834-46. doi: 10.1002/cyto.a.20608

4. Lee G, Finn W, Scott C. Statistical file matching of flow cytometry data. *J Biomed Inform.* (2011) 44:663-76. doi: 10.1016/j.jbi.2011.03.004

5. Leite Pereira A, Lambotte O, Le Grand R, Cosma A, Tchitchek N. CytoBackBone: an algorithm for merging of phenotypic information from different cytometric profiles. *Bioinformatics*. (2019). doi: 10.1093/bioinformatics/btz212. [Epub ahead of print].

6. Berger EA, Murphy PM, Farber JM. Chemokine receptors as HIV-1 coreceptors: roles in viral entry, tropism, and disease. *Ann Rev Immunol.* (2002) 17:657-700. doi: 10.1146/annurev.immunol.17.1.657

7. Simon V, Ho D, Karim Q. HIV/AIDS epidemiology, pathogenesis, prevention and treatment. *Lancet*. (2010) 368:489-504. doi: 10.1016/S0140-6736(06)69157-5

8. Campbell JH, Hearps AC, Martin GE, Williams KC, Crowe SM. The importance of monocytes and macrophages in HIV pathogenesis,treatment,and cure. *AIDS*. (2014) 28:2175-87. doi: 10.1097/QAD.0000000000000408

9. Simmons RP, Scully EP, Groden EE, Arnold KB, Chang JJ, Lane K, et al. HIV-1 infection induces strong production of IP-10 through TLR7/9-dependent pathways. *AIDS*. (2013) 27:2505-17. doi: 10.1097/01.aids.0000432455.06476.bc

10. Dutertre CA, Amraoui S, DeRosa A, Jourdain JP, Vimeux L, Goguet M, et al. Pivotal role of M-DC8(+) monocytes from viremic HIV-infected patients in TNFalpha overproduction in response to microbial products. *Blood*. (2012) 120:2259-68. doi: 10.1182/blood-2012-03-418681

11. Guo H, Gao J, Taxman DJ, Ting JPY, Su L. HIV-1 infection induces interleukin-1β production via TLR8 protein-dependent and NLRP3 inflammasome mechanisms in human monocytes. *J Biol Chem.* (2014) 289:21716-26. doi: 10.1074/jbc.M114.566620

12. Boasso A, Hardy AW, Landay AL, Martinson JL, Anderson SA, Dolan MJ, et al. PDL-1 upregulation on monocytes and T cells by HIV via type I interferon: Restricted expression of type I interferon receptor

FUNDING

The IDMIT infrastructure is supported by the French government Programme d'Investissements d'Avenir (PIA) under grant ANR-11-INBS-0008 and grant ANR-10-EQPX-02-01 (FlowCyTech facility). NT was supported by fellowships from the ANRS.

ACKNOWLEDGMENTS

We warmly thank the HIV-infected patients and
healthy donors. We are grateful to the physicians at
the hospitals.

SUPPLEMENTARY MATERIAL

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

by CCR5-expressing leukocytes. *Clin Immunol.* (2008) 129:132-44.
doi: 10.1016/j.clim.2008.05.009

13. Morlat P, Roussillon C, Henard S, Salmon D, Bonnet F, Cacoub P, et al. Causes of death among HIV-infected patients in France in 2010 (national survey). *AIDS*. (2014) 28:1181-91. doi: 10.1097/QAD.0000000000000222

14. Brenchley JM, Price DA, Schacker TW, Asher TE, Silvestri G, Rao S, et al. Microbial translocation is a cause of systemic immune activation in chronic HIV infection. *Nat Med.* (2007) 12:1365-71. doi: 10.1038/nm1511

15. Leite Pereira A, Tchitchek N, Marcos Lopez E, Lambotte O, Le Grand R, Cosma A. A high-resolution mass cytometry analysis reveals a delay of cytokines production after TLR4 or TLR7/8 engagements in HIV-1 infected humans. *Cytokine*. (2018) 111:97-105. doi: 10.1016/j.cytok.2018.08.018

16. Velu V, Shetty RD, Larsson M, Shankar EM. Role of PD-1 co-inhibitory pathway in HIV infection and potential therapeutic options. *Retrovirology*. (2015) 12:1-17. doi: 10.1186/s12977-015-0144-x

17. Porichis F, Hart MG, Zupkosky J, Barblu L, Kwon DS, McMullen A, et al. Differential Impact of PD-1 and/or interleukin-10 blockade on HIV-1-specific CD4 T cell and antigen-presenting cell functions. *J Virol.* (2014) 88:2508-18. doi: 10.1128/JVI.02034-13

18. Trono D, Van Lint C, Rouzioux C, Verdin E, Barre-Sinoussi F, Chun T-W, et al. HIV Persistence and the prospect of long-term drug-free remissions for HIV-infected individuals. *Science*. (2010) 329:174-80. doi: 10.1126/science.1191047

19. Kuller LH, Tracy R, Belloso W, De Wit S, Drummond F, Lane HC, et al. Inflammatory and coagulation biomarkers and mortality in patients with HIV infection. *PLoS Med*. (2008) 5:1496-508. doi: 10.1371/journal.pmed.0050203

20. Rodger AJ, Fox Z, Lundgren JD, Kuller L, Boesecke C, Skoutelis A, et al. Activation and coagulation biomarkers are independent predictors of the development of opportunistic disease in patients with HIV infection. *J Infect Dis.* (2009) 200:973-83. doi: 10.1086/605447

21. Elhmouzi-Younes J, Palgen JL, Tchitchek N, Delandre S, Namet I, Bodinham CL, et al. In depth comparative phenotyping of blood innate myeloid leukocytes from healthy humans and macaques using mass cytometry. *Cytometry Part A*. (2017) 91:969-82. doi: 10.1002/cyto.a.23107

22. Rahman AH, Tordesillas L, Berin MC. Heparin reduces nonspecific eosinophil staining artifacts in mass cytometry experiments. *Cytometry Part A*. (2016) 89:601-7. doi: 10.1002/cyto.a.22826

23. Gautreau G, Pejoski D, Cosma A, Le Grand R, Beignon AS, Tchitchek N. SPADEVizR: an R package for visualization, analysis
well as the cytometric profiles obtained from the confirmatory
single panel are available in the FlowRepository database through
the accession IDs: FR-FCM-Z26Z and FR-FCM-Z26Y.

## ETHICS STATEMENT

This experiment was approved by the Comité de Protection des
Personnes (CPP) Ile de France VII, under protocol number PP
14-003. All subjects gave written informed consent to participate
in this study.

## AUTHOR CONTRIBUTIONS

AL, NT, and AC: conceptualization and methodology, investigation. AL, NT, AC, RL, and OL: validation and writing—review and editing. AL and NT: formal analysis and writing—original draft. AL and OL: resources. AC and RL: funding acquisition. AC, RL, and OL: supervision.

## REFERENCES

1. Bendall SC, Nolan GP, Roederer M, Chattopadhyay PK. A deep profiler's guide to cytometry. *Trends Immunol.* (2012) 33:323-32. doi: 10.1016/j.it.2012.02.010

2. O'Neill K, Aghaepour N, Parker J, Hogge D, Karsan A, Dalal B, Brinkman RR. Deep profiling of multitube flow cytometry data. *Bioinformatics*. (2015) 31:1623-31. doi: 10.1093/bioinformatics/btv008

3. Pedreira CE, Costa ES, Barrena S, Lecrevisse Q, Almeida J, Van Dongen JJM, et al. Generation of flow cytometry data files with a potentially infinite number of dimensions. *Cytometry Part A*. (2008) 73:834-46. doi: 10.1002/cyto.a.20608

4. Lee G, Finn W, Scott C. Statistical file matching of flow cytometry data. *J Biomed Inform.* (2011) 44:663-76. doi: 10.1016/j.jbi.2011.03.004

5. Leite Pereira A, Lambotte O, Le Grand R, Cosma A, Tchitchek N. CytoBackBone: an algorithm for merging of phenotypic information from different cytometric profiles. *Bioinformatics*. (2019). doi: 10.1093/bioinformatics/btz212. [Epub ahead of print].

6. Berger EA, Murphy PM, Farber JM. Chemokine receptors as HIV-1 coreceptors: roles in viral entry, tropism, and disease. *Ann Rev Immunol.* (2002) 17:657-700. doi: 10.1146/annurev.immunol.17.1.657

7. Simon V, Ho D, Karim Q. HIV/AIDS epidemiology, pathogenesis, prevention and treatment. *Lancet*. (2010) 368:489-504. doi: 10.1016/S0140-6736(06)69157-5

8. Campbell JH, Hearps AC, Martin GE, Williams KC, Crowe SM. The importance of monocytes and macrophages in HIV pathogenesis,treatment,and cure. *AIDS*. (2014) 28:2175-87. doi: 10.1097/QAD.0000000000000408

9. Simmons RP, Scully EP, Groden EE, Arnold KB, Chang JJ, Lane K, et al. HIV-1 infection induces strong production of IP-10 through TLR7/9-dependent pathways. *AIDS*. (2013) 27:2505-17. doi: 10.1097/01. aids.0000432455.06476.bc

10. Dutertre CA, Amraoui S, DeRosa A, Jourdain JP, Vimeux L, Goguet M, et al. Pivotal role of M-DC8(+) monocytes from viremic HIV-infected patients in TNFalpha overproduction in response to microbial products. *Blood*. (2012) 120:2259-68. doi: 10.1182/blood-2012-03-418681

11. Guo H, Gao J, Taxman DJ, Ting JPY, Su L. HIV-1 infection induces interleukin-1β production via TLR8 protein-dependent and NLRP3 inflammasome mechanisms in human monocytes. *J Biol Chem.* (2014) 289:21716-26. doi: 10.1074/jbc.M114.566620

12. Boasso A, Hardy AW, Landay AL, Martinson JL, Anderson SA, Dolan MJ, et al. PDL-1 upregulation on monocytes and T cells by HIV via type I interferon: Restricted expression of type I interferon receptor

## FUNDING

The IDMIT infrastructure is supported by the French government Programme d'Investissements d'Avenir (PIA) under grant ANR-11-INBS-0008 and grant ANR-10-EQPX-02-01 (FlowCyTech facility). NT was supported by fellowships from the ANRS.

## ACKNOWLEDGMENTS

We warmly thank the HIV-infected patients and healthy donors. We are grateful to the physicians at the hospitals.

## SUPPLEMENTARY MATERIAL

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

by CCR5-expressing leukocytes. *Clin Immunol.* (2008) 129:132-44. doi: 10.1016/j.clim.2008.05.009

13. Morlat P, Roussillon C, Henard S, Salmon D, Bonnet F, Cacoub P, et al. Causes of death among HIV-infected patients in France in 2010 (national survey). *AIDS*. (2014) 28:1181-91. doi: 10.1097/QAD.0000000000000222

14. Brenchley JM, Price DA, Schacker TW, Asher TE, Silvestri G, Rao S, et al. Microbial translocation is a cause of systemic immune activation in chronic HIV infection. *Nat Med.* (2007) 12:1365-71. doi: 10.1038/nm1511

15. Leite Pereira A, Tchitchek N, Marcos Lopez E, Lambotte O, Le Grand R, Cosma A. A high-resolution mass cytometry analysis reveals a delay of cytokines production after TLR4 or TLR7/8 engagements in HIV-1 infected humans. *Cytokine*. (2018) 111:97-105. doi: 10.1016/j.cyto.2018. 08.018

16. Velu V, Shetty RD, Larsson M, Shankar EM. Role of PD-1 co-inhibitory pathway in HIV infection and potential therapeutic options. *Retrovirology*. (2015) 12:1-17. doi: 10.1186/s12977-015-0144-x

17. Porichis F, Hart MG, Zupkosky J, Barblu L, Kwon DS, McMullen A, et al. Differential Impact of PD-1 and/or interleukin-10 blockade on HIV-1-specific CD4 T cell and antigen-presenting cell functions. *J Virol*. (2014) 88:2508-18. doi: 10.1128/JVI.02034-13

18. Trono D, Van Lint C, Rouzioux C, Verdin E, Barre-Sinoussi F, Chun T-W, et al. HIV Persistence and the prospect of long-term drug-free remissions for HIV-infected individuals. *Science*. (2010) 329:174-80. doi: 10.1126/science.1191047

19. Kuller LH, Tracy R, Belloso W, De Wit S, Drummond F, Lane HC, et al. Inflammatory and coagulation biomarkers and mortality in patients with HIV infection. *PLoS Med*. (2008) 5:1496-508. doi: 10.1371/journal.pmed.0050203

20. Rodger AJ, Fox Z, Lundgren JD, Kuller L, Boesecke C, Skoutelis A, et al. Activation and coagulation biomarkers are independent predictors of the development of opportunistic disease in patients with HIV infection. *J Infect Dis.* (2009) 200:973-83. doi: 10.1086/ 605447

21. Elhmouzi-Younes J, Palgen JL, Tchitchek N, Delandre S, Namet I, Bodinham CL, et al. In depth comparative phenotyping of blood innate myeloid leukocytes from healthy humans and macaques using mass cytometry. *Cytometry Part A*. (2017) 91:969-82. doi: 10.1002/cyto.a.23107

22. Rahman AH, Tordesillas L, Berin MC. Heparin reduces nonspecific eosinophil staining artifacts in mass cytometry experiments. *Cytometry Part A*. (2016) 89:601-7. doi: 10.1002/cyto.a. 22826

23. Gautreau G, Pejoski D, Cosma A, Le Grand R, Beignon AS, Tchitchek N. SPADEVizR: an R package for visualization, analysis
and integration of SPADE results. *Bioinformatics*. (2016) 33:779-81.
doi: 10.1093/bioinformatics/btw708

24. Amir ED, Davis KL, Tadmor MD, Simonds EF, Jacob H. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. *Nat Biotechnol.* (2014) 31:545–52. doi: 10.1038/nbt.2594

25. Kawai T, Akira S. The roles of TLRs, RLRs and NLRs in pathogen recognition. *Int Immunol.* (2009) 21:317–37. doi: 10.1093/intimm/dxp017

26. Uematsu S, Akira S. Toll-like receptors and innate immunity. *Adv Immunol.* (2006) 78:1–56. doi: 10.2222/jsv.56.1

27. Aziz KE, Wakefield D. Modulation of endothelial cell expression of ICAM-1, E-selectin, and VCAM-1 by beta-estradiol, progesterone, and dexamethasone. *Cell Immunol.* (1996) 167:79–85. doi: 10.1006/cimm.1996.0010

28. Zhang H, Park Y, Wu J, Chen X, Lee S, Yang J, Dellsperger KC, Zhang C. Role of TNF-alpha in vascular dysfunction. *Clin Sci.* (2009) 116:219–30. doi: 10.1042/CS20080196

29. Donati KDG, Rabagliati R, Iacoviello L, Cauda R. Review HIV infection, HAART, and endothelial adhesion molecules: current perspectives. *Lancet Infect Dis.* (2004) 4:213–22. doi: 10.1016/S1473-3099(04) 00971-5

30. Stacey AR, Norris PJ, Qin L, Haygreen EA, Taylor E, Heitman J, et al. Induction of a striking systemic cytokine cascade prior to peak viremia in acute human immunodeficiency virus type 1 infection, in contrast to more modest and delayed responses in acute hepatitis B and C virus infections. *J Virol.* (2009) 83:3719–33. doi: 10.1128/JVI. 01844-08

31. Shi X, Sims MD, Hanna MM, Xie M, Gulick PG, Zheng Y, Basson MD, et al. Neutropenia during HIV infection: adverse consequences and remedies. *Int Rev Immunol.* (2014) 33:511–36. doi: 10.3109/08830185.2014. 893301

32. Hensley-McBain T, Klatt NR. The dual role of neutrophils in HIV infection. *Curr HIV/AIDS Rep.* (2018) 15:206–21. doi: 10.1007/s11904-018- 0370-7

33. Mazzolini J, Herit F, Niedergang F. Inhibition of phagocytosis in HIV-1 - infected macrophages relies on Nef-dependent alteration of focal delivery of recycling compartments. *Blood.* (2017) 115:4226–37. doi: 10.1182/blood-2009-12-259473

34. Debaisieux S, Lachambre S, Gross A, Mettling C, Besteiro S, Yezid H, et al. HIV-1 Tat inhibits phagocytosis by preventing the recruitment of Cdc42 to the phagocytic cup. *Nat Commun.* (2015) 6:6211. doi: 10.1038/ncomms7211

35. Garagiola DM, Huard TK, Lobuglio AF. Comparison of monocyte and alveolar macrophage dependent cellular cytotoxicity and fc-receptor activity. *Cell Immunol.* (1981) 370:359–70. doi: 10.1016/0008-8749(81) 90487-1

36. Yeap WH, Wong KL, Shimasaki N, Chi E, Teo Y. CD16 is indispensable for antibody- dependent cellular cytotoxicity by human monocytes. *Nat Publish Group.* (2016) 6:34310. doi: 10.1038/srep34310

37. Young NJ De, Gill PG. Monocyte antibody-dependent cellular cytotoxicity in cancer patients. *Cancer Immunol Immunother.* (1984) 18:54–8. doi: 10.1007/BF00205400

38. Lehmann MH, Masanetz S, Kramer S, Erfle V. HIV-1 Nef upregulates CCL2/MCP-1 expression in astrocytes in a myristoylation- and calmodulin- dependent manner. *J Cell Sci.* (2006) 119:4520–30. doi: 10.1242/jcs.03231

39. Marban C, Forouzanfar F, Ait-Ammar A, Fahmi F, El Mekdad H, Daouad F, et al. Targeting the brain reservoirs: Toward an HIV cure. *Front Immunol.* (2016) 7:00397. doi: 10.3389/fimmu.2016.00397

40. Miller E, Bhardwaj N, Sinai M. Dendritic cell dysregulation during HIV-1 infection. *Immunol Rev.* (2017) 254:170–89. doi: 10.1111/imr. 12082

41. Aerts-Toegaert C, Heirman C, Tuyaerts S, Corthals J, Aerts JL, Bonehill A, et al. CD83 expression on dendritic cells and T cells: correlation with effective immune responses. *Eur J Immunol.* (2007) 37:686–95. doi: 10.1002/eji.200636535

42. Harman AN, Wilkinson J, Bye CR, Stern JL, Nicholle M, Lai J, et al. HIV Induces maturation of monocyte-derived dendritic cells and langerhans cells. *J Immunol.* (2006) 177:7103–113. doi: 10.4049/jimmunol.177. 10.7103

43. Chen L, Zhu Y, Zhang G, Gao C, Zhong W, Zhang X. CD83-stimulated monocytes suppress T-cell immune responses through production of prostaglandin E2. *Proc Natl Acad Sci USA.* (2011) 108:18778–83. doi: 10.1073/pnas.1018994108

44. Barblu L, MacHmach K, Gras C, Delfraissy JF, Boufassa F, Leal M, et al. Plasmacytoid dendritic cells (pDCs) from HIV controllers produce interferon-and differentiate into functional killer pDCs under HIV activation. *J Infect Dis.* (2012) 206:790–801. doi: 10.1093/infdis/ jis384

45. Hernández JC, Stevenson M, Latz E, Urcuqui-Inchima S. HIV type 1 infection up-regulates TLR2 and TLR4 expression and function *in vivo* and *in vitro*. *AIDS Res Hum Retroviruses*. (2012) 28:1313–28. doi: 10.1089/aid.2011. 0297

46. Lester RT, Yao XD, Ball TB, McKinnon LR, Kaul R, Wachihi C, et al. Toll-like receptor expression and responsiveness are increased in viraemic HIV-1 infection. *AIDS*. (2008) 22:685–94. doi: 10.1097/QAD.0b013e3282f4de35

47. Martinelli E, Cicala C, Van Ryk D, Goode DJ, Macleod K, Arthos J, Fauci AS. HIV-1 gp120 inhibits TLR9-mediated activation and IFN-[alpha] secretion in plasmacytoid dendritic cells. *Proc Natl Acad Sci USA.* (2007) 104:3396–401. doi: 10.1073/pnas.0611353104

48. Reuven EM, Ali M, Rotem E, Schwarzter R, Gramatica A, Futerman AH, et al. The HIV-1 envelope transmembrane domain binds TLR2 through a distinct dimerization motif and inhibits TLR2-mediated responses. *PLoS Pathog.* (2014) 10:e1004248. doi: 10.1371/journal.ppat. 1004248

49. Ahmed N, Hayashi T, Hasegawa A, Furukawa H, Okamura N, Chida T, et al. Suppression of human immunodeficiency virus type 1 replication in macrophages by commensal bacteria preferentially stimulating Toll-like receptor 4. *J General Virol.* (2010) 91:2804–13. doi: 10.1099/vir.0. 022442-0

50. Novis CL, Archin NM, Buzon MJ, Verdin E, Round JL, Lichterfeld M, et al. Reactivation of latent HIV-1 in central memory CD4+ T cells through TLR-1/2 stimulation. *Retrovirology*. (2013) 10:119. doi: 10.1186/1742-4690- 10-119

51. Mikulak J, Oriolo F, Zaghi E, Di Vito C, Mavilio D. Natural killer cells in HIV-1 infection and therapy. *AIDS*. (2017) 31:2317–30. doi: 10.1097/QAD.0000000000001645

52. Scully E, Alter G. NK Cells in HIV disease. *Curr HIV/AIDS Rep.* (2016) 13:85–94. doi: 10.1007/s11904-016-0310-3

53. Schafer JL, Li H, Evans TI, Estes JD, Reeves RK. Accumulation of Cytotoxic CD16 + NK cells in simian immunodeficiency virus-infected lymph nodes associated with *in situ* differentiation and functional anergy. *J Virol.* (2015) 89:6887–94. doi: 10.1128/JVI.00660-15

54. Langford SE, Ananworanich J, Cooper DA. Predictors of disease progression in HIV infection: a review. *AIDS Res Ther.* (2007) 4:11. doi: 10.1186/1742-6405-4-11

55. Grant EJ, Nüssing S, Sant S, Clemens EB, Kedzierska K. The role of CD27 in anti-viral T-cell immunity. *Curr Opin Virol.* (2017) 22:77–88. doi: 10.1016/j.coviro.2016.12.001

56. Hintzen RQ, de Jong R, Lens SM, Brouwer M, Baars P, van Lier RA. Regulation of CD27 expression on subsets of mature T-lymphocytes. *J Immunol.* (1993) 151:2426–35.

57. Borst J. Targeting the T-cell co-stimulatory CD27 / CD70 pathway in cancer immunotherapy : rationale and potential. *Immunotherapy.* (2015) 7:655–67. doi: 10.2217/imt.15.32

58. Reinwald S, Wiethe C, Westendorf AM, Breloer M, Probst-Kepper M, Fleischer B, et al. CD83 Expression in CD4+ T cells modulates inflammation and autoimmunity. *J Immunol.* (2008) 180:5890–7. doi: 10.4049/jimmunol.180. 9.5890

59. Breloer M, Kretschmer B, Lüthje K, Ehrlich S, Ritter U, Bickert T, et al. CD83 is a regulator of murine B cell function *in vivo*. *Eur J Immunol.* (2007) 37:634–48. doi: 10.1002/eji.200636852

60. Damouche A, Lazonne T, Avettand-f V, Huot N, Dejucq-rainsford N, Satie A, et al. Adipose tissue is a neglected viral reservoir and an inflammatory site during chronic HIV and SIV infection. *PLoS Pathog.* (2015) 11:1–28. doi: 10.1371/journal.ppat.1005153
and integration of SPADE results. *Bioinformatics*. (2016) 33:779-81.
doi: 10.1093/bioinformatics/btw708

24. Amir ED, Davis KL, Tadmor MD, Simonds EF, Jacob H. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. *Nat Biotechnol.* (2014) 31:545–52. doi: 10.1038/nbt.2594

25. Kawai T, Akira S. The roles of TLRs, RLRs and NLRs in pathogen recognition. *Int Immunol.* (2009) 21:317–37. doi: 10.1093/intimm/dpx017

26. Uematsu S, Akira S. Toll-like receptors and innate immunity. *Adv Immunol.* (2006) 78:1–56. doi: 10.2222/jsv.56.1

27. Aziz KE, Wakefield D. Modulation of endothelial cell expression of ICAM-1, E-selectin, and VCAM-1 by beta-estradiol, progesterone, and dexamethasone. *Cell Immunol.* (1996) 167:79–85. doi: 10.1006/cimm.1996.0010

28. Zhang H, Park Y, Wu J, Chen X, Lee S, Yang J, Dellsperger KC, Zhang C. Role of TNF-alpha in vascular dysfunction. *Clin Sci.* (2009) 116:219–30. doi: 10.1042/CS20080196

29. Donati KDG, Rabagliati R, Iacoviello L, Cauda R. Review HIV infection, HAART, and endothelial adhesion molecules: current perspectives. *Lancet Infect Dis.* (2004) 4:213–22. doi: 10.1016/S1473-3099(04) 00971-5

30. Stacey AR, Norris PJ, Qin L, Haygreen EA, Taylor E, Heitman J, et al. Induction of a striking systemic cytokine cascade prior to peak viremia in acute human immunodeficiency virus type 1 infection, in contrast to more modest and delayed responses in acute hepatitis B and C virus infections. *J Virol.* (2009) 83:3719–33. doi: 10.1128/JVI. 01844-08

31. Shi X, Sims MD, Hanna MM, Xie M, Gulick PG, Zheng Y, Basson MD, et al. Neutropenia during HIV infection: adverse consequences and remedies. *Int Rev Immunol.* (2014) 33:511–36. doi: 10.3109/08830185.2014. 893301

32. Hensley-McBain T, Klatt NR. The dual role of neutrophils in HIV infection. *Curr HIV/AIDS Rep.* (2018) 15:206–21. doi: 10.1007/s11904-018- 0370-7

33. Mazzolini J, Herit F, Niedergang F. Inhibition of phagocytosis in HIV-1 - infected macrophages relies on Nef-dependent alteration of focal delivery of recycling compartments. *Blood.* (2017) 115:4226–37. doi: 10.1182/blood-2009-12-259473

34. Debaisieux S, Lachambre S, Gross A, Mettling C, Besteiro S, Yezid H, et al. HIV-1 Tat inhibits phagocytosis by preventing the recruitment of Cdc42 to the phagocytic cup. *Nat Commun.* (2015) 6:6211. doi: 10.1038/ncomms7211

35. Garagiola DM, Huard TK, Lobuglio AF. Comparison of monocyte and alveolar macrophage dependent cellular cytotoxicity and fc-receptor activity. *Cell Immunol.* (1981) 370:359–70. doi: 10.1016/0008-8749(81) 90487-1

36. Yeap WH, Wong KL, Shimasaki N, Chi E, Teo Y. CD16 is indispensable for antibody- dependent cellular cytotoxicity by human monocytes. *Nat Publish Group.* (2016) 6:34310. doi: 10.1038/srep34310

37. Young NJ De, Gill PG. Monocyte antibody-dependent cellular cytotoxicity in cancer patients. *Cancer Immunol Immunother.* (1984) 18:54–8. doi: 10.1007/BF00205400

38. Lehmann MH, Masanetz S, Kramer S, Erfle V. HIV-1 Nef upregulates CCL2/MCP-1 expression in astrocytes in a myristoylation- and calmodulin-dependent manner. *J Cell Sci.* (2006) 119:4520–30. doi: 10.1242/jcs.03231

39. Marban C, Forouzanfar F, Ait-Ammar A, Fahmi F, El Mekdad H, Daouad F, et al. Targeting the brain reservoirs: Toward an HIV cure. *Front Immunol.* (2016) 7:00397. doi: 10.3389/fimmu.2016.00397

40. Miller E, Bhardwaj N, Sinai M. Dendritic cell dysregulation during HIV-1 infection. *Immunol Rev.* (2017) 254:170–89. doi: 10.1111/imr. 12082

41. Aerts-Toegaert C, Heirman C, Tuyaerts S, Corthals J, Aerts JL, Bonehill A, et al. CD83 expression on dendritic cells and T cells: correlation with effective immune responses. *Eur J Immunol.* (2007) 37:686–95. doi: 10.1002/eji.200636535

42. Harman AN, Wilkinson J, Bye CR, Stern JL, Nicholle M, Lai J, et al. HIV Induces maturation of monocyte-derived dendritic cells and langerhans cells. *J Immunol.* (2006) 177:7103–113. doi: 10.4049/jimmunol.177. 10.7103

43. Chen L, Zhu Y, Zhang G, Gao C, Zhong W, Zhang X. CD83-stimulated monocytes suppress T-cell immune responses through production of prostaglandin E2. *Proc Natl Acad Sci USA.* (2011) 108:18778–83. doi: 10.1073/pnas.1018994108

44. Barblu L, MacHmach K, Gras C, Delfraissy JF, Boufassa F, Leal M, et al. Plasmacytoid dendritic cells (pDCs) from HIV controllers produce interferon-and differentiate into functional killer pDCs under HIV activation. *J Infect Dis.* (2012) 206:790–801. doi: 10.1093/infdis/ jis384

45. Hernández JC, Stevenson M, Latz E, Urcuqui-Inchima S. HIV type 1 infection up-regulates TLR2 and TLR4 expression and function *in vivo* and *in vitro*. *AIDS Res Hum Retroviruses*. (2012) 28:1313–28. doi: 10.1089/aid.2011. 0297

46. Lester RT, Yao XD, Ball TB, McKinnon LR, Kaul R, Wachihi C, et al. Toll-like receptor expression and responsiveness are increased in viraemic HIV-1 infection. *AIDS*. (2008) 22:685–94. doi: 10.1097/QAD.0b013e3282f4de35

47. Martinelli E, Cicala C, Van Ryk D, Goode DJ, Macleod K, Arthos J, Fauci AS. HIV-1 gp120 inhibits TLR9-mediated activation and IFN-[alpha] secretion in plasmacytoid dendritic cells. *Proc Natl Acad Sci USA.* (2007) 104:3396–401. doi: 10.1073/pnas.0611353104

48. Reuven EM, Ali M, Rotem E, Schwarzter R, Gramatica A, Futerman AH, et al. The HIV-1 envelope transmembrane domain binds TLR2 through a distinct dimerization motif and inhibits TLR2-mediated responses. *PLoS Pathog.* (2014) 10:e1004248. doi: 10.1371/journal.ppat. 1004248

49. Ahmed N, Hayashi T, Hasegawa A, Furukawa H, Okamura N, Chida T, et al. Suppression of human immunodeficiency virus type 1 replication in macrophages by commensal bacteria preferentially stimulating Toll-like receptor 4. *J General Virol.* (2010) 91:2804–13. doi: 10.1099/vir.0. 022442-0

50. Novis CL, Archin NM, Buzon MJ, Verdin E, Round JL, Lichterfeld M, et al. Reactivation of latent HIV-1 in central memory CD4+ T cells through TLR-1/2 stimulation. *Retrovirology*. (2013) 10:119. doi: 10.1186/1742-4690- 10-119

51. Mikulak J, Oriolo F, Zaghi E, Di Vito C, Mavilio D. Natural killer cells in HIV-1 infection and therapy. *AIDS*. (2017) 31:2317–30. doi: 10.1097/QAD.0000000000001645

52. Scully E, Alter G. NK Cells in HIV disease. *Curr HIV/AIDS Rep.* (2016) 13:85–94. doi: 10.1007/s11904-016-0310-3

53. Schafer JL, Li H, Evans TI, Estes JD, Reeves RK. Accumulation of Cytotoxic CD16 + NK cells in simian immunodeficiency virus-infected lymph nodes associated with *in situ* differentiation and functional anergy. *J Virol.* (2015) 89:6887–94. doi: 10.1128/JVI.00660-15

54. Langford SE, Ananworanich J, Cooper DA. Predictors of disease progression in HIV infection: a review. *AIDS Res Ther.* (2007) 4:11. doi: 10.1186/1742-6405-4-11

55. Grant EJ, Nüssing S, Sant S, Clemens EB, Kedzierska K. The role of CD27 in anti-viral T-cell immunity. *Curr Opin Virol.* (2017) 22:77–88. doi: 10.1016/j.coviro.2016.12.001

56. Hintzen RQ, de Jong R, Lens SM, Brouwer M, Baars P, van Lier RA. Regulation of CD27 expression on subsets of mature T-lymphocytes. *J Immunol.* (1993) 151:2426–35.

57. Borst J. Targeting the T-cell co-stimulatory CD27 / CD70 pathway in cancer immunotherapy : rationale and potential. *Immunotherapy.* (2015) 7:655–67. doi: 10.2217/imt.15.32

58. Reinwald S, Wiethe C, Westendorf AM, Breloer M, Probst-Kepper M, Fleischer B, et al. CD83 Expression in CD4+ T cells modulates inflammation and autoimmunity. *J Immunol.* (2008) 180:5890–7. doi: 10.4049/jimmunol.180. 9.5890

59. Breloer M, Kretschmer B, Lüthje K, Ehrlich S, Ritter U, Bickert T, et al. CD83 is a regulator of murine B cell function *in vivo*. *Eur J Immunol.* (2007) 37:634–48. doi: 10.1002/eji.200636852

60. Damouche A, Lazonne T, Avettand-f V, Huot N, Dejucq-rainsford N, Satie A, et al. Adipose tissue is a neglected viral reservoir and an inflammatory site during chronic HIV and SIV infection. *PLoS Pathog.* (2015) 11:1–28. doi: 10.1371/journal.ppat.1005153
61. Estes JD, Kityo C, Ssali F, Swainson L, Makamdop KN, Del Prete GQ, et al. Defining total body AIDS virus burden: implications for curative strategies. *Nat Med*. (2018) 23:1271–6. doi: 10.1038/nm.4411

62. Hattab S, Guiguet M, Carcelain G, Fourati S, Guihot A, Autran B, et al. Soluble biomarkers of immune activation and inflammation in HIV infection: impact of 2 years of effective first-line combination antiretroviral therapy. *HIV Med*. (2015) 16:553–62. doi: 10.1111/hiv.12257

63. Noel N, Boufassa F, Lécuroux C, Saez-Cirion A, Bourgeois C, Dunyach-Remy C, et al. Elevated IP10 levels are associated with immune activation and low CD4+ T-cell counts in HIV controller patients. *AIDS*. (2014) 28:467–76. doi: 10.1097/QAD.0000000000000174

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Leite Pereira, Tchitchek, Lambotte, Le Grand and Cosma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
61. Estes JD, Kityo C, Ssali F, Swainson L, Makamdop KN, Del Prete GQ, et al. Defining total body AIDS virus burden: implications for curative strategies. *Nat Med*. (2018) 23:1271–6. doi: 10.1038/nm.4411

62. Hattab S, Guiguet M, Carcelain G, Fourati S, Guihot A, Autran B, et al. Soluble biomarkers of immune activation and inflammation in HIV infection: impact of 2 years of effective first-line combination antiretroviral therapy. *HIV Med*. (2015) 16:553–62. doi: 10.1111/hiv.12257

63. Noel N, Boufassa F, Lécuroux C, Saez-Cirion A, Bourgeois C, Dunyach-Remy C, et al. Elevated IP10 levels are associated with immune activation and low CD4+ T-cell counts in HIV controller patients. *AIDS*. (2014) 28:467–76. doi: 10.1097/QAD.0000000000000174

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Leite Pereira, Tchitchek, Lambotte, Le Grand and Cosma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.