Materials and Methods
Cell lines – Chronically infected HIV-1 promonocytic (U1, ARP-165) cell lines were originally derived by limiting dilution cloning of U937 surviving an acute infection with HIV-1 by Folks et al.51. The cells were obtained through the NIH HIV Reagent Program, Division of AIDS, NIAID, NIH: Anti-Human CD34 Hybridoma (PR18), ARP-165, contributed by Dr. Thomas Folks. The U1 cells were cultured in RPMI 1640 containing 2.0mM of L-glutamine and 10% heat inactivated fetal bovine serum (Lonza Bioscience, Morrisville, NC) and maintained in log phase with >98% viability prior to plating at 2x106/ml in 24 well plates and stimulation.
Treatments in vitro – All concentrations were optimized prior to experimentation. Dopamine hydrochloride (H8502, Sigma Aldrich, St. Louis, MO) was added at 1 and 10μM. Latency reversal agents were Ingenol-3-angelate (PEP005) at 0.1 and 10uM (Sigma-Aldrich, St.Louis, MO), Bromodomain and extraterminal domain inhibitor 151 (BETi) at 1 and 10uM (Tocris, Minneapolis, MN) and Phorbol-myristate acetate (PMA) at 10 and 100ng/ml (Sigma-Aldrich). Selective DA receptor agonists SKF38393 (DRD1/5), Ropinirole (DRD2), Pramipexole (DRD3), and PD168077 (DRD4) (all from Tocris) were added at 1 and 10uM. Selective DRD antagonists SCH23390 (DRD1/5), Haloperidol hydrochloride (DRD2), PNU177864 (DRD3) and L745870 hydrochloride (DRD4) (all from Tocris), were added at 1uM to cultures, alone or combined with 1 and 10uM DA. Recombinant human heterodimer MRP8/14 proteins (S100A8/S100A9, R&D Systems, Minneapolis, MN) was used at 0.5, 1 and 5 uM. RAGE antagonist FPS-ZM1 (Tocris) was used alone or simultaneously with other stimuli, at 1 and 5 uM. All stimulations, with single or combined compounds, were tested at 24hrs.
Quantification of HIV replication – HIV p24 was measured in the supernatants 24hrs after stimulation using ELISA kit (BioXpress International, Frederick, MD), and viral RNA was measured by digital droplet PCR using normalized msTatRev and skGag primers, performed at the University of California San Diego Translational Virology Core.
Detection of MRP8/14 in culture supernatants – Culture supernatants were collected after stimulation, and levels of MRP8/14 were measured using LEGENDplex 8.0 microbeads immunoassay (Biolegend, San Diego, CA) by flow cytometry52-54, using Human MRP8/14 Capture beads A5 13X, following manufacturer’s instructions, and acquired in a CytoFlex S (Beckman Coulter, Brea, CA), and analyzed in cloud-based LEGENDplex software (Biolegend).
Detection of MRP8/14 by ELISA in protein preparations – Following stimulation, the cells were collected from tissue culture plates, washed once with ice cold PBS and then homogenized in ice-cold RIPA buffer containing 1 mM protease inhibitor PMSF and 10% Complete phosphatase inhibitor (Roche, Indianapolis, IN, USA). Protein concentration in the lysate was determined with the BCA assay (Pierce Biotechnology, Waltham, MA). Cytosolic and membrane fractions were separated with MEM-PER Plus kit (ThermoFisher, Waltham, MA). MRP8/14 levels were tested using Legend Max ELISA kit (BioLegend) following the protocol from the manufacturer.
Transcriptional Profiles – The total cellular RNA was extracted from the U1 cell cultures at a density of 106/ml, 12 hours following exposure to DA or selective agonists. Extraction was performed using the RNAasey Mini Kit (Qiagen, Germantown, MD) according to the manufacturer’s protocol. Total RNA was quantified by measuring the total absorbance at OD260/OD280 nm in Nanodrop (Thermo Fisher Scientific, Watham, MA) and the quality of total RNA was monitored using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA).
Transcriptional changes were examined in biological triplicates of DA, DRD1/5 and DRD4 agonists, using the Agilent Whole Human genome platform, performed by Phalanx Biotech (San Diego, CA). Normalized values were compared using log Fold Change between conditions, with alpha=0.05.
Systems Biology and Visualization – Transcriptional changes caused by DA and DRD agonists, compared to vehicle controls, were visualized using Genemania55,56 in Cytoscape 3.3 interface57,58, with links through Pathway Commons algorithm (Blue lines), and public databases of shared protein domains (Red lines). Blue shapes represent down-regulation, and yellow shapes represent upregulation. Circles represent p<0.05 compared to vehicle.
Human cohorts - Participants included 102 adults enrolled by the University of California San Diego’s HIV Neurobehavioral Research Program (HNRP) and Translational Methamphetamine Research Center (TMARC) under informed consent and approved protocols. Cognitive scores were collected using previously described protocols
59. Peripheral blood leukocytes were collected during visits and archived. The subjects selected for this study were males, between 35 – 49 years old, due to characteristics of the cohort and for decreasing confounders. The participants were divided based on HIV serostatus (HIV+/-) and on Meth use (METH+/-). METH+ was defined as met lifetime DSM-IV criteria for methamphetamine abuse or dependence. A cross sectional design assembled the following groups: HIV-METH- (n=27), HIV+METH- (n=25), HIV-METH+ (n=25) and HIV+METH+ (n=25). Exclusion criteria was history of non-HIV-related neurological, medical or psychiatric disorders that affect brain function (e.g., schizophrenia, traumatic brain injury, epilepsy), learning disabilities, or dementia. Major depressive disorder (MDD) and polysubstance use were minimized but not excluded due to the high prevalence of depression and polysubstance use in Meth users.
Table 1 shows demographics, health and psychiatric characteristics. CSF and plasma viral loads were measured by RT-PCR in a CLIA-certified laboratory. CD4 and CD8 T-lymphocyte counts were measured by flow cytometry for PWH. Nadir CD4 levels were taken from medical records, study-obtained values, or self-report.
Table 1.
– Cohort demographic, health and psychiatric characteristics.
Table 1.
– Cohort demographic, health and psychiatric characteristics.
|
HIV-/METH- |
HIV+/METH- |
HIV-/METH+ |
HIV+/METH+ |
ANOVA p values |
N |
27 |
25 |
25 |
25 |
|
|
Mean |
STD |
Mean |
STD |
Mean |
STD |
Mean |
STD |
|
Age |
37.68 |
9.17 |
38.8 |
6.94 |
36.16 |
9.26 |
36.32 |
6.44 |
0.4980 |
Education |
13.68 |
2.37 |
13.24 |
2.61 |
12.7 |
2.45 |
12.08 |
2.81 |
0.0882 |
Global T score |
49.70 |
6.64 |
46.19 |
5.66 |
45.94 |
5.7 |
46.17 |
6.77 |
0.0524 |
CD4 Nadir |
815 |
205 |
287 |
219 |
918 |
184 |
311 |
.23 |
<0.0001 |
CD4/CD8 Ratio |
1.84 |
1.29 |
0.53 |
0.37 |
2.48 |
1.28 |
0.56 |
0.26 |
<0.0001 |
Duration of infection (yrs) |
NA |
NA |
7.26 |
5.73 |
NA |
NA |
6.18 |
6.98 |
0.4431 |
% Black |
7 |
3 |
3 |
0 |
|
% Hispanic |
5 |
6 |
2 |
5 |
|
% Asian |
1 |
1 |
1 |
2 |
|
% White |
14 |
15 |
19 |
18 |
|
Detectable Plasma VL (% of total) |
NA |
60 |
NA |
56 |
0.0774* |
Detectable CSF VL (% of total) |
NA |
48 |
NA |
58 |
0.0293* |
LT Alcohol dep (% of total) |
60 |
45 |
68 |
52 |
0.0407 |
LT Cocaine dep (% of total) |
20 |
32 |
36 |
24 |
|
LT Opioid dep (% of total) |
4 |
4.5 |
24 |
12 |
<0.0001 |
LT MDD (% of total) |
16 |
72 |
48 |
52 |
|
Peripheral Leukocytes specimens – Archived specimens were thawed in fetal bovine serum, washed and viability was counted using disposable hemocytometers (Bulldog Bio, Portsmouth, NH). Cell numbers were adjusted for flow cytometry.
Flow Cytometry - Washed cells had adjusted concentrations of 105/100ul and were resuspended in HBSS without phenol red, containing 2% fetal bovine serum and 0.2% sodium azide, and stained with pre-defined concentrations of subset, function and activation-specific antibodies. The surface expression of CCR5, CD38, PDL1, MRP8/14 and RAGE was performed on cells that were incubated from 2 to 24hrs with DA or DRD agonists, or on straight cryopreserved human leukocytes from human subjects. The cells were stained with a PerCP-labeled anti-CD11b (clone M1/70, Biolegend, San Diego, CA), APC-labeled anti-human CCR5 (CD195, clone 3A9, BD Pharmingen, San Diego, CA), PE-Cy5-labelled anti-human CD38 (HIT2, Thermo Fisher Scientific), PE-Cy7-labelled anti human PDL1 (Thermo Fisher Scientific), FITC-labelled anti-human MRP8/14 (Origene, Rockville, MD) and PE-labelled anti-RAGE EPR21171 (Abcam, Wantham, MA). The detection of intracellular p24 was performed with All the cells were fixed in 4% paraformaldehyde and kept at 4oC and protected from light until acquisition. Cell acquisition was performed in a CytoFlex S benchtop platform (Beckman Coulter, Indianapolis, IN), and analyzed using FlowJo software (FlowJo LLC, Ashland, OR).
Brain specimens– Formalin-fixed Paraffin-embedded tissue sections of human post-mortem basal ganglia specimens were provided by the National NeuroAIDS Tissue Consortium (NNTC) upon request number R593, with Institutional IRB approval (IRB-18-001-MCM). The specimens were selected among all male HIV+ cases (receiving ART and divided in groups following the same criteria described above, and a post-mortem interval below 48hrs (n=5/group).
Immunohistochemistry - Seven µm sections were stained using antibody anti MRP8/14 (ab17050, mouse anti-human monoclonal 27E10, Abcam, Cambridge, UK), performed as previously described 60. Colorimetric development was performed with NovaRed chromogen (Vector Laboratories, Burlingame, CA and AbD Serotec, Raleigh, NC respectively) followed by Gill’s hematoxylin counterstain (Sigma-Aldrich).
RNAscope In Situ Hybridization for HIV (vRNA) Detection in Brain specimens - RNAScope 2.5 HD assay (Advanced Cell Diagnostics, ACD) was performed. Briefly, pre-treatment was performed with RNAscope 2.5 HD Detection Kit (RED) (Cat# 322360), RNAscope 2.5 Pretreat Reagents-H202 and Protease Plus (Cat#322330), RNAscope Target Retrieval (Cat#322000) and RNAscope Wash Buffer (Cat#310091), following manual instructions. Probe set V-HIV1-clade B-C3 (Cat#425531-C3) targeting different segments within the gag-pol region.
Image processing and quantification - Images were captured using the light feature of a Zeiss AXIO Observer.Z1 (Carl Zeiss AG, Oberkochen, Germany). Digitalized images in tiff format were opened in Fiji/ImageJ (National Institute of Health, USA) and converted into 8-bit binary masks for quantification and normalization to the total area.
Statistical Analysis – For gene expression analysis, scores were calculated as the log2 normalized expression of each gene. Variance due to noise for each score was estimated in a linear mixed model. Pairwise comparisons and false discovery rate were calculated for each group. Fitted value models were used to determine variable and intercept effects. All biomarker analysis was performed using ANOVA followed by multiple comparisons performed using Tukey HD. Predictor screening was used to determine specific effects of individual variables, named HIV and METH, and their interactions. In human specimens, mixed models were used to identify interactive effects of demographics, substance use, detectable cerebrospinal fluid (CSF) viral load, and cognitive Global Deficit scores (GDS) on cell surface markers. All statistical analysis was performed in JMP Pro 15.2.0 software (SAS Institute Inc., Cary, NC, USA).
Results
- 1.
HIV latency in myeloid cells is partially reverted by DA stimulation via DRD1.
U1 cells, an established model of HIV latency in the myeloid compartment
51, were used to examine how the neurotransmitter DA affects innate immune phenotypes associated with latency reversal and increased viral replication, mimicking neuroimmune effects of Meth in the brain. Compared to latency reversal agents such as PEP005, BETi and PMA, which act on 80 to 100% of all the cells in the culture, DA stimulation at concentrations that mimic levels found in the brain environment during substance abuse
36 consistently promoted replication in 10-15% of the U1 cells, as determined by levels of p24 in the supernatant after 24hrs (
Figure 1A). The same effect was obtained by the selective D1-like receptor (DRD1/5) agonist SKF38393 (SKF), even combined with the DRD2 agonist Ropinirole (Rop), which had no effect on p24 levels, or with DRD3 agonist Pramipexole (Pra). On the other hand, selective DRD4 stimulation with PD168077 (PD) did not revert latency and prevented reversal by SKF (
Figure 2B). Conversely, the stimulation with DA in the presence of a DRD4 antagonist L745 significantly increased p24, while other antagonists had limited or no effect on DA-induced latency reversal (
Figure 1C). Together, the results in
Figure 1B suggest that D1-like receptors (DRD1/5) and DRD4 have opposing effects on latency. These effects were further characterized.
Figure 1.
– Effects of DA and DRDs on p24 levels- HIV p24 levels were measured by ELISA in U1 cell culture supernatants 24hrs following the stimulation with A) DA and classic latency reversal agents PEP005 (0.1 and 10uM), BETi (1 and 10uM) and PMA (10 and 100ng/ml), B) Selective DRD agonists SKF38393 (SKF, DRD1/5, 1 and 10uM), ropinirole (Rop, DRD2, 1 and 10uM), pramipexole (Pra, DRD3, 1 and 10uM), and PD168077 (PD, DRD4, 1 and 10uM), and C) Selective DRD antagonists SCH23390 (SCH, DRD1/5), haloperidol hydrochloride (Hal, DRD2), PNU177864 (PNU, DRD3) and L745870 hydrochloride (L745, DRD4), were added at 1uM to cultures, alone or combined with 1 and 10uM DA. Representative experiments out of 3 independent assays performed in quadruplicate show the mean ± SD. NT=Non-treated vehicle control. *p<0.05 in ANOVA followed by Tukey HSD.
Figure 1.
– Effects of DA and DRDs on p24 levels- HIV p24 levels were measured by ELISA in U1 cell culture supernatants 24hrs following the stimulation with A) DA and classic latency reversal agents PEP005 (0.1 and 10uM), BETi (1 and 10uM) and PMA (10 and 100ng/ml), B) Selective DRD agonists SKF38393 (SKF, DRD1/5, 1 and 10uM), ropinirole (Rop, DRD2, 1 and 10uM), pramipexole (Pra, DRD3, 1 and 10uM), and PD168077 (PD, DRD4, 1 and 10uM), and C) Selective DRD antagonists SCH23390 (SCH, DRD1/5), haloperidol hydrochloride (Hal, DRD2), PNU177864 (PNU, DRD3) and L745870 hydrochloride (L745, DRD4), were added at 1uM to cultures, alone or combined with 1 and 10uM DA. Representative experiments out of 3 independent assays performed in quadruplicate show the mean ± SD. NT=Non-treated vehicle control. *p<0.05 in ANOVA followed by Tukey HSD.
Figure 2.
– Surface markers of population and activation status in latent U1 cells treated with DA and DRD agonists. The expression of population and activation markers was measured on the U1 cell surface by flow cytometry upon exposure to DA and selective receptor agonists. U1 cells were either (A, G, M) in control conditions (vehicle/non-treated), (B, H, N) treated with 10uM of DA, or treated 1uM of DRD selective agonists (C, I, O) SKF38393 (DRD1/5), (D, J, P) Ropinirole (DRD2), (E, K, Q) Pramipexole (DRD3), or (F, L, R) PD168077 (DRD4). The plots are representative of 3 independent experiments. (A-L) The CD14+ CD16low subset in shown in dark blue and the CD14+CD16high activated subset in shown light blue, displaying a three-dimensional distribution of CCR5 (x-axis, PDL1(y-axis), and CD38 (z-axis). (M-R) Flat scatter plots of surface CCR5 (x-axis) and CD38 (y-axis) in ungated live cells. A gate on the CCR5+CD38+ was used to estimate the percentage of the total cells expressing these markers.
Figure 2.
– Surface markers of population and activation status in latent U1 cells treated with DA and DRD agonists. The expression of population and activation markers was measured on the U1 cell surface by flow cytometry upon exposure to DA and selective receptor agonists. U1 cells were either (A, G, M) in control conditions (vehicle/non-treated), (B, H, N) treated with 10uM of DA, or treated 1uM of DRD selective agonists (C, I, O) SKF38393 (DRD1/5), (D, J, P) Ropinirole (DRD2), (E, K, Q) Pramipexole (DRD3), or (F, L, R) PD168077 (DRD4). The plots are representative of 3 independent experiments. (A-L) The CD14+ CD16low subset in shown in dark blue and the CD14+CD16high activated subset in shown light blue, displaying a three-dimensional distribution of CCR5 (x-axis, PDL1(y-axis), and CD38 (z-axis). (M-R) Flat scatter plots of surface CCR5 (x-axis) and CD38 (y-axis) in ungated live cells. A gate on the CCR5+CD38+ was used to estimate the percentage of the total cells expressing these markers.
- 2.
DA stimulation via DRD1 triggers a novel CD16hi subset and novel transcriptional signatures.
Key surface markers that have been previously described as reporters of activation and inflammation, or in association with changes in replication and myeloid response to DA
44, were measured by flow cytometry, with the goal of identifying phenotypic changes that parallel DA-induced increases in p24 levels, which were also present in DRD1/5 selective stimulation. These markers were CD11b, CD14 and CD16 (population and activation markers
61), CCR5
44, CD38
62-64 and PDL1
65 (inflammation, viral replication, activation and immune senescence). We found that both DA and the D1-like receptor (DRD1/5) agonist increased the levels of CD16 in a subset of cells, which also expressed CCR5. This subset of cells activated by DA comprised 3-5% of the total cells and also co-expressed higher levels of CD38, compared to controls (
Figure 2B, 2H and 2N). This observation has suggested that exposure of latent cells to DA promotes the development of a subset of latent cells with pro-inflammatory markers (
Figure 2). Importantly, cells with a similar phenotype were also identified in U1 cells stimulated by the D1-like receptor agonist SKF38393 (
Figure 2N and 2O), but not by other DA receptor agonists (
Figure 2P, 2Q and 2R).
One distinction between DA and D1-like receptors (DRD1/5) selective stimulation was related to PDL1 levels, which were predominantly low in CD16high cells following DA, but high in the same subset following SKF38393 stimulation (
Figure 2B and 2C). However, when U1 cells were gated in the CCR5+ population rather than CD16 (
Figure 3), we observed that both DA and SKF38393, as well as the stimulation with the DRD2 selective agonist Ropinirole, showed increased CD38 (
Figure 3B, 3C, and 3D), as well as PDL1 expression (
Figure 3H, 3I and 3J), in association with CD16hi activation. Cells stimulated with the DRD4 agonist showed a minimal increase in CCR5 or CD38 (
Figure 2F, 2L and 2R), and also had low CD38 and PDL1 expression within the CCR5-gated CD16hi subset (
Figure 3F and 3L).
Figure 3.
– Surface markers of inflammation in gated CCR5+ U1 cells treated with DA and DRD agonists. The expression of population and activation markers was measured by flow cytometry on U1 cells upon exposure to DA and selective receptor agonists. The markers examined in combination were (A, B, C, D, E, F) CD16 vs. CD38 or (G, H, I, J, K, L) CD16 vs. PDL1. U1 cells were either (A, G) in control conditions (vehicle/non-treated), (B, H) treated with 10uM of DA, or treated 1uM of DRD selective agonists (C, I) SKF38393 (DRD1/5), (D, J) Ropinirole (DRD2), (E, K) Pramipexole (DRD3), or (F, L) PD168077 (DRD4). The plots are representative of 3 independent experiments.
Figure 3.
– Surface markers of inflammation in gated CCR5+ U1 cells treated with DA and DRD agonists. The expression of population and activation markers was measured by flow cytometry on U1 cells upon exposure to DA and selective receptor agonists. The markers examined in combination were (A, B, C, D, E, F) CD16 vs. CD38 or (G, H, I, J, K, L) CD16 vs. PDL1. U1 cells were either (A, G) in control conditions (vehicle/non-treated), (B, H) treated with 10uM of DA, or treated 1uM of DRD selective agonists (C, I) SKF38393 (DRD1/5), (D, J) Ropinirole (DRD2), (E, K) Pramipexole (DRD3), or (F, L) PD168077 (DRD4). The plots are representative of 3 independent experiments.
To further characterize changes associated with the onset of a CD16hi novel population, we screened transcriptional signatures triggered by DA via DRD1/5 and in correlation with latency reversal, as well as via DRD4 using gene array. We focused on changes that were common between DA and the DRD1/5 agonist SKF38393, which differed from changes caused by the DRD4 agonist PD168077, which presented significant and opposing effects on p24 levels. Genes that were significantly changed were further tested for their involvement in latency reversal by DA, and as signatures of HIV uncontrolled replication in the context of Meth.
To further characterize changes in U1 cells stimulated with DA, or DRD1/5 and DRD4 selective agonists in parallel with changes in p24, we examined global changes in the transcriptome, 12hrs after stimulation. We found that DA stimulation significantly upregulated eight genes in U1 cells. Of these, the calcium binding and macrophage inhibitory factor related proteins S100A8 and S100A9, known respectively as myeloid-related protein (MRP) 8 and 14, were significantly increased by both DA and DRD1/5 stimulation (
Figure 4A). DA increased MRP8 by 1.67-fold (p=0.007) and DRD1/5 stimulation increased it by 1.89-fold (p=0.005) (
Figure 4A). MRP14 was increased by 1.58-fold by both DA and DRD1/5 stimulation (p=0.021 and p=0.029, respectively. Importantly, MRP8 and MRP14 form a heterodimer that was previously described in plasma as a marker of aggravated HIV infection
66, in neurological disorders
67,68, and in pro-inflammatory responses by endothelial cells
69. Other genes increased by DA included PDE4B (Phosphodiesterase 4B, 1.67-fold), CCL2 (1.52-fold), FTH1 (Ferritin heavy polypeptide 1, 1.66-fold), KYNU (kynureninase, 1.68-fold) and HMOX (Heme oxygenase, 2.47-fold) (
Figure 4A). DRD1/5 increased RAS dexamethasone-induced protein 1 (RASD1) by 2.8-fold (p=0.00002) as an exclusive signature. DRD4 upregulated 16 exclusive genes: Retinoic acid receptor gamma (RARG), NADH dehydrogenase (NDUFAF5, 1.8-fold), midline 2 (MID2, 1.5-fold), ribosomal protein L10 (RPL10, 1.65-fold), RAR-related orphan receptor A (RORA, 1.64-fold), mucin 22 (MUC22, 1.53-fold), MT-RNR2-like 4 (MTRNR2L4, 1.7-fold), translin-associated factor X interacting protein 1 (TSNAXIP1, 2.03-fold), phospholipase A2, group IIE (PLA2G2E, 1.52-fold), ankyrin repeat and SOCS box containing 16 (ASB16, 1.52-fold), disrupted in schizophrenia 1 (DISC1, 1.64-fold), helicase with zinc finger 2, transcriptional coactivator (HELZ2, 1.58-fold), chemokine C-C motif 6 (CCR6, 1.55-fold), DPY19L2, Kruppel-like factor 11 (KLF11, 1.56-fold) and chromosome 14 open reading frame 37 (C14orf37, 1.51-fold). Growth factor receptor-bound protein 10 (GRB10) was increased by both DRD1/5 and DRD4 stimulation (5.02-fold, p=0.03 and 10.03-fold, p=0.000006), but not by DA. The enhancer of mRNA decapping 3 (EDC3) was significantly increased by DA, DRD1/5 as well as DRD4 (2.24, 5.39 and 5.26-fold, respectively). Apart from co-expression, pathway-based or functional interactions between these genes that were increased exclusively by DRD4 could not be identified by GeneMania. However, a transcription factor usage prediction using iRegulon suggested that these changes may have occurred with significant contribution of the transcription factors FOXA2 (7.6%), AP4 (7.4%) and GATA1 (7%). Using both DAVID Bioinformatics Database
44,70 and Gene Ontology Resource
71, genes upregulated by DA were annotated to leukocyte migration (p=0.00000002), positive regulation of inflammation (p=0.0003), inflammatory response (p=0.0008), regulation of NFkB transcription factor activity (p=0.00094). Genes upregulated by DRD4 selective stimulation were functionally annotated to transcription initiation from RNA polymerase (p=3.3E-25).
Figure 4.
– Venn diagrams showing gene signatures identified in response to DA and DA receptor agonists and in correlation with replication (DA and DRD1/5 agonist SKF38393) or latency (DRD4 agonist PD168077). Following group comparisons for each individual gene, log10 fold change between DA and Vehicle (DA), SKF38393 and Vehicle (DRD1/5), or PD168077 and Vehicle (DRD4). Significantly changed genes (FDR p<0.05) were selected and a cut off of 1.5-fold was used to identify genes that were A) Upregulated or B) Downregulated by stimuli, compared to controls. Venn diagrams indicate the genes fitting in these criteria in each group, and genes that overlapped between stimuli. When the number of genes in a given comparison was too large to be listed, the number of genes was scripted and the 10 most downregulated were listed. Systems biology and visualization strategies were used to determine the relationship and interactions between the genes downregulated by DRD1/5 and DRD4 as well as by both. The legend shows connectors’ colors indicating the different levels of interactions between genes. C) Interactions between the genes downregulated exclusively by DRD1/5 selective stimulation. D) Interactions between genes exclusively downregulated by DRD4. E) Interactions between genes commonly downregulated by DRD1/5 and DRD4 selective stimulation. Genes that are changed but do not have reported interactions do not show connectors. A complete list of differences can be found in supplementary
Table S1.
Figure 4.
– Venn diagrams showing gene signatures identified in response to DA and DA receptor agonists and in correlation with replication (DA and DRD1/5 agonist SKF38393) or latency (DRD4 agonist PD168077). Following group comparisons for each individual gene, log10 fold change between DA and Vehicle (DA), SKF38393 and Vehicle (DRD1/5), or PD168077 and Vehicle (DRD4). Significantly changed genes (FDR p<0.05) were selected and a cut off of 1.5-fold was used to identify genes that were A) Upregulated or B) Downregulated by stimuli, compared to controls. Venn diagrams indicate the genes fitting in these criteria in each group, and genes that overlapped between stimuli. When the number of genes in a given comparison was too large to be listed, the number of genes was scripted and the 10 most downregulated were listed. Systems biology and visualization strategies were used to determine the relationship and interactions between the genes downregulated by DRD1/5 and DRD4 as well as by both. The legend shows connectors’ colors indicating the different levels of interactions between genes. C) Interactions between the genes downregulated exclusively by DRD1/5 selective stimulation. D) Interactions between genes exclusively downregulated by DRD4. E) Interactions between genes commonly downregulated by DRD1/5 and DRD4 selective stimulation. Genes that are changed but do not have reported interactions do not show connectors. A complete list of differences can be found in supplementary
Table S1.
We also observed downregulation of genes by DA and in particular by DRD1/5 or DRD4 selective stimulation (
Figure 4B). DA significantly downregulated a total of 9 genes (p<0.05), of which 2 were exclusive to DA [sulfite oxidase (SUOX, -1.5-fold) and cytochrome P450, family 17, subfamily A, polypeptide 1 (CYP17A1, -1.5-fold)]; 2 were common between DA and DRD1/5 stimulation [olfactory receptor 51I1 (OR51l1, -2.12- and -1.97-fold, respectively) and CUB and sushi multiple domains 2 (CSMD2, -1.5-, and -1.58-fold, respectively)]; 3 were common between DA and DRD4 [Leucine Rich Repeat Containing 41 (LRRC41, -3.25- and -1.79-fold, respectively), Solute Carrier Organic Anion Transporter Family Member 1A2 (SLCO1A2, -1.65- and -1.91-fold, respectively) and Solute Carrier Family 4 Member 1 (SLC4A1, -1.65- and -2.39-fold, respectively)]; and 2 were upregulated by all 2 stimuli [olfactory receptor 52E4 (OR52E4, -1.9-fold, -1.52-fold and -1.94-fold, respectively) and Ras protein-specific guanine nucleotide-releasing factor 1 (RASGRF1, -1.84-, -1.93- and -2.39-fold, respectively)](
Figure 4B). Interestingly, DRD1/5 stimulation significantly downregulated 92 exclusive genes, which were connected by co-expression, physical and genetic interactions, co-localization and pathway, color-coded as in legend (
Figure 4C) and assigned to the Rap1-signaling pathway (p=0.0032). DRD4 downregulated 64 exclusive genes, which also showed interactions (
Figure 4D) and were functionally assigned to regulation of transcription from the RNA Pol II promoter (p=0.002). Moreover, the downregulation of 38 genes was common to both DRD1/5 and DRD4, with identifiable interactions (
Figure 4E) and annotated to long term depression (p=0.0072). Among the genes downregulated by DRD1/5, DRD4 or both,
Figure 4B shows a short list of genes with a role in inflammation. See supplementary
Table S1 for a complete list of all significantly affected genes.
Data integration, analysis and visualization were performed in Cytoscape 3.7 software as previously described
44,70. Pathway-based interactions were sorted using GeneMania
55,56 to identify significantly affected gene networks annotated to biological processes and pathways, confirmed using DAVID Bioinformatics Database and Gene Ontology Resources, as described
44,70.
Figure 5 shows the two most significant networks with overrepresented gene changes as a result of DA, DRD1/5 and DRD4 agonist stimulations combined, compared to vehicle, although with different behaviors. These networks had upregulated (yellow shades) and downregulated genes (blue shades). Genes involved in response to chemical stimulus and control of transcription through RNA polymerase (p<0.0004) (
Figure 5A, 5B and 5C) and inflammation and chemokine receptors (p<0.00001) (
Figure 5D, 5E and 5F) exhibited similarities between DA and DRD1/5 stimulation, but also differences. Other identified pathways included cell motility (p<0.0004) and neurotransmission and cardiac function (p<0.0005), with overlapping representation with the networks in
Figure 5. The pathway annotated to response to stimulus and regulation of transcription (
Figure 5A, 5B and 5C) showed differences in the expression of specific genes that paralleled reversed latency by DA and DRD1/5 stimulation and replicated the visualization in
Figure 4, with increased levels of S100A8 and S100A9, or MRP8 and MRP14, by DA and DRD1/5 (
Figure 5A and 5B). The visualization of gene interactions in
Figure 5A, 5B and 5C suggested that the transcriptional levels MRP8 and MRP14 were inversely correlated with the transcription of retinoic acid receptor gamma (RARG) levels. Regarding the pathway annotated to inflammation and chemokine receptor (
Figure 5D, 5E and 5F) expression patterns were mixed.
Figure 5.
– Pathways significantly affected by DA and by selective DRD1/5 and DRD4 agonists. U1 cells stimulated with (A, D) 10uM DA, or (B, E) D1-like DRD (DRD1/5) agonist SKF389393, or (C, F) with the DRD4 agonist PD168077 were compared to vehicle treated, non-stimulated cells. Transcriptional changes were examined by Agilent Whole Human genome platform. Normalized values were compared using log Fold Change between conditions, and then visualized using Genemania in Cytoscape 3.3 interface, with pathway-based (Blue lines) and shared protein domains-based interactions (Red lines). Blue shapes represent down-regulation, and yellow shapes represent upregulation. Circles are p<0.05 compared to vehicle. Two main gene networks were identified and annotated to (A, B, C) Response to stimulus and control of transcription, and (D, E, F) Inflammation.
Figure 5.
– Pathways significantly affected by DA and by selective DRD1/5 and DRD4 agonists. U1 cells stimulated with (A, D) 10uM DA, or (B, E) D1-like DRD (DRD1/5) agonist SKF389393, or (C, F) with the DRD4 agonist PD168077 were compared to vehicle treated, non-stimulated cells. Transcriptional changes were examined by Agilent Whole Human genome platform. Normalized values were compared using log Fold Change between conditions, and then visualized using Genemania in Cytoscape 3.3 interface, with pathway-based (Blue lines) and shared protein domains-based interactions (Red lines). Blue shapes represent down-regulation, and yellow shapes represent upregulation. Circles are p<0.05 compared to vehicle. Two main gene networks were identified and annotated to (A, B, C) Response to stimulus and control of transcription, and (D, E, F) Inflammation.
- 3.
MRP8 and MRP14 are DA signatures that reverses HIV latency via RAGE
Given that both DA and the DRD1/5 agonist increased the transcription of MRP8 and MRP14, we examined whether the heterodimer is induced at the protein level in U1 cells following DA exposure (
Figure 6)
72,73. We found that DA both at 1uM and at 10uM, was able to significantly increase the secretion of MRP8/14 in the U1 supernatant, detectable at 12 hrs after stimulation. Interestingly, the DRD1/5 agonist at the lowest concentration of 0.1uM, and the DRD4 agonist at the highest 1uM concentration, were able to trigger the release of MRP8/14 (
Figure 6A). Other selective agonists did not have an effect (not shown). Similarly, the concentration of MRP8/14 in the cytoplasm was increased by DA, via both receptors (
Figure 6B). The increase in MRP8/14 complex was also observed on the cell surface following DA stimulation (both 1 and 10uM), as well as with both DRD1/5 and DRD4 selective agonists (
Figure 6B, and 6C), but not with other DRD agonists (not shown). Interestingly, the stimulation via DRD4 significantly decreased the surface expression of the MRP8/14 receptor RAGE
74 (
Figure 6F). Neither DA or the DRD1/5 agonist SKF389393 affect the expression of the RAGE on the cell surface (
Figure 6F). To test the significance of these findings for latency, we examined whether MRP8/14 was able to increase expression p24 in U1 cells, mimicking the effects of DA and SKF38393 on p24 (
Figure 6G). For that, U1 cells were exposed to different doses of the recombinant MRP8/14. We also tested the effect of MRP8/14 via its two described receptors, the receptor for advanced glycation end products or RAGE
74 and toll-like receptor 4 or TLR4
75,76. Indeed, MRP8/14 was able to increase the levels of p24 detectable at 24hrs, in a dose-dependent manner (
Figure 6G). The inhibition of RAGE by a selective antagonist FPS-ZM1
77 prevented the effect of MRP8/14 on p24 also in a dose-dependent manner (
Figure 6G), without affecting cell viability (not shown). The selective TLR4 inhibitor TAK242, on the other hand, did not inhibit MRP8/14-induced increase in p24, unless at a higher dose of 10uM, when it also decreased cell viability from 98% (±7) to 56% (±10) (p=0.003). This suggested that MRP8/14 are the transcriptional signatures induced by DA that influence latency, reversing it via binding to RAGE. Interestingly, although both DRD1/5 and DRD4 selective agonists increased secreted, cytosol and surface MRP8/14, DRD4 selective stimulation caused a drastic decrease in RAGE expression (
Figure 6F), indicating its ability to prevent activation and latency reversal via RAGE. This observation could provide a mechanistic link to explain relative differences in the latency response to DA and to why DRD4 selective stimulation may prevent latency reversal, as seen in
Figure 1.
Figure 6.
–MRP8/14 is upregulated by DA and reverses latency via RAGE. A) Levels of MRP8/14 in culture supernatants measured by bead-based immunoassay, 24hrs after stimulation with DA, the DRD1/5 agonist SKF38393 (SKF), or the DRD4 agonist PD168077 (PD) at the indicated concentrations, normalized by volume and calculated using standard curve. B) Concentration of MRP8/14 in the cytosolic fraction extracted 24 hrs after stimulation, normalized to total protein. (C-F) By flow cytometry, the surface levels of DA-induced biomarkers were detectable. Representative histograms indicate the geometric mean fluorescence (GeoMean) of (C, D) MRP8/14 and (E, F) RAGE in cells stimulated with (C, E) 10uM DA (D, F) with the DRD1/5 agonist SKF389393, or with the DRD4 agonist PD168077, compared to controls, as indicated in the legend. G) HIV p24 levels in the culture supernatant measured by ELISA 24 hrs after stimulation with recombinant MRP8/14 at indicated doses, in the presence or absence of FPS-ZM1 (RAGE antagonist) or TAK242 (TLR4 antagonist). Tukey HD multiple comparisons, * p<0.05 compared to non-treated (NT) control, # p<0.05 compared to 1uM MRP8/14.
Figure 6.
–MRP8/14 is upregulated by DA and reverses latency via RAGE. A) Levels of MRP8/14 in culture supernatants measured by bead-based immunoassay, 24hrs after stimulation with DA, the DRD1/5 agonist SKF38393 (SKF), or the DRD4 agonist PD168077 (PD) at the indicated concentrations, normalized by volume and calculated using standard curve. B) Concentration of MRP8/14 in the cytosolic fraction extracted 24 hrs after stimulation, normalized to total protein. (C-F) By flow cytometry, the surface levels of DA-induced biomarkers were detectable. Representative histograms indicate the geometric mean fluorescence (GeoMean) of (C, D) MRP8/14 and (E, F) RAGE in cells stimulated with (C, E) 10uM DA (D, F) with the DRD1/5 agonist SKF389393, or with the DRD4 agonist PD168077, compared to controls, as indicated in the legend. G) HIV p24 levels in the culture supernatant measured by ELISA 24 hrs after stimulation with recombinant MRP8/14 at indicated doses, in the presence or absence of FPS-ZM1 (RAGE antagonist) or TAK242 (TLR4 antagonist). Tukey HD multiple comparisons, * p<0.05 compared to non-treated (NT) control, # p<0.05 compared to 1uM MRP8/14.
- 4.
DA signatures expressed on peripheral innate immune cell surface are biomarkers of HIV infection and detectable CSF viral load in the context of Meth use
We examined the translational value of our in vitro findings by examining whether the observed signatures triggered by DA could be identified in human subjects and whether they can distinguish HIV+ individuals that are Meth users compared to non-Meth users (
Figure 7). Moreover, given the involvement of these biomarkers in latency reversal, we compared individuals with detectable or undetectable viral load in the cerebrospinal fluid (CSF) at the time-point of specimen collection. For that, specimens from 102 individuals from the Translational Methamphetamine Research Center (TMARC) cohorts (UCSD) were distributed in 4 groups as shown in
Table 1, with HIV-METH- (n=27), HIV+METH- (n=25), HIV-METH+ (n=25) and HIV+METH+ (n=25) subjects, homogeneous for sex (all males due to characteristics of the cohort and sample size), age and education to minimize confounders.
Table 1 shows demographics and other characteristics of the cohort, including the percentage of individuals with detectable viral load in plasma and CSF. We tested the expression of CD38, CCR5, PDL1, RAGE and MRP8/14 by flow cytometry on the surface of peripheral blood leukocytes.
Figure 7.
– Characterization of peripheral monocytes from HIV+ Meth users using surface expression of DA signature proteins. Using flow cytometry, we measured the expression of DA signatures and inflammatory markers in peripheral blood monocytes from HIV+ and HIV- subjects that were Meth users (METH+) or not (METH-). The monocytes were gated using CD11b, CD14 and CD16 antibodies, and the geometric mean fluorescence of (A) CD38, (B) CCR5, (C) PDL1, (D) MRP8/14, and (E) RAGE were calculated in FlowJo software. Multiple comparisons were performed using Tukey HD, with indicated p values. Dots indicate outliers. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
Figure 7.
– Characterization of peripheral monocytes from HIV+ Meth users using surface expression of DA signature proteins. Using flow cytometry, we measured the expression of DA signatures and inflammatory markers in peripheral blood monocytes from HIV+ and HIV- subjects that were Meth users (METH+) or not (METH-). The monocytes were gated using CD11b, CD14 and CD16 antibodies, and the geometric mean fluorescence of (A) CD38, (B) CCR5, (C) PDL1, (D) MRP8/14, and (E) RAGE were calculated in FlowJo software. Multiple comparisons were performed using Tukey HD, with indicated p values. Dots indicate outliers. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
The value of these markers varied to indicate effects of HIV, of Meth, or their interaction (
Figure 7), as well as to distinguish the group with detectable CSF viral load, particularly within the HIV+ METH+ group (
Figure 8). For instance, CD38 expression was overall higher in subjects of the HIV+METH+ group with detectable CSF viral load, compared to undetectable viral load within the same group (p=0.043) (
Figure 8A), although this marker was not overall affected by HIV or Meth (
Figure 7A). While CCR5 did not distinguish any variables (
Figure 7B and 8B), PDL1 expression was significantly decreased by HIV infection, regardless of Meth status (
Figure 7C). MRP8/14 was lower in CD16hi cells from HIV- Meth users (
Figure 7D). In HIV+ subjects, an interactive effect of HIV, Meth and detectable CSF viral load was indicated, with a trend to higher expression of MRP8/14 in Meth users with detectable viral load. A similar trend was observed with the expression of RAGE (
Figure 8E).
Figure 8.
– Surface expression of DA signatures on monocytes from HIV+ Meth and CSF viral load. Using flow cytometry, we measured the expression of DA signatures and inflammatory markers in peripheral blood monocytes from HIV+ subjects that were Meth users (METH+) or not (METH-), further divided into groups with detectable or undetectable CSF viral load. The monocytes were gated using CD11b, CD14 and CD16 antibodies, and the geometric mean fluorescence of (A) CD38, (B) CCR5, (C) PDL1, (D) MRP8/14, and (E) RAGE were calculated in FlowJo software. Multiple comparisons were performed using Tukey HD, with indicated p values. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
Figure 8.
– Surface expression of DA signatures on monocytes from HIV+ Meth and CSF viral load. Using flow cytometry, we measured the expression of DA signatures and inflammatory markers in peripheral blood monocytes from HIV+ subjects that were Meth users (METH+) or not (METH-), further divided into groups with detectable or undetectable CSF viral load. The monocytes were gated using CD11b, CD14 and CD16 antibodies, and the geometric mean fluorescence of (A) CD38, (B) CCR5, (C) PDL1, (D) MRP8/14, and (E) RAGE were calculated in FlowJo software. Multiple comparisons were performed using Tukey HD, with indicated p values. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
We also tested whether the expression of combined markers could detect subject groups more efficiently than single markers (
Figure 9). We found that HIV+METH+ subjects with undetectable CSF viral loads had significantly lower CCR5 expression in the CD11b-gated CD16hi CD38hi subset compared to HIV+METH+ subjects with detectable viral load (p=0.0021) and to HIV+METH- subjects with undetectable viral load (
Figure 9A). We also found that the HIV+METH+ group with detectable CSF viral load had significantly lower PDL1 expression in the same cell subset compared to subjects with undetectable viral load in the CSF (p=0.003) or with METH- groups (p=0.03 and p<0.0001, detectable and undetectable CSF viral loads respectively), (
Figure 9B).
Figure 9.
– Combined expression of surface markers on peripheral monocytes. Using flow cytometry, CD11b+ CD16hi activated cells were gated and further analyzed for their expression levels of CD38 combined with A) CCR5 and B) PDL1. Bar colors represent expression level ranges as indicated in legends. Multiple comparisons were performed using Tukey HD, with indicated p values. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
Figure 9.
– Combined expression of surface markers on peripheral monocytes. Using flow cytometry, CD11b+ CD16hi activated cells were gated and further analyzed for their expression levels of CD38 combined with A) CCR5 and B) PDL1. Bar colors represent expression level ranges as indicated in legends. Multiple comparisons were performed using Tukey HD, with indicated p values. Mixed models indicated effects of HIV, METH and CSF viral load interactions.
- 5.
MRP8/14 and HIV RNA are increased in HIV-infected post-mortem brain specimens in association with Meth use
We also examined whether MRP8/14+ cells could be identified in the brains of HIV+ Meth users, using post-mortem brain specimens of basal ganglia obtained from the National NeuroAIDS Tissue Consortium (NNTC). For that, HIV-METH- controls (n=5), HIV+METH- (n=5), HIV-METH+ (n=5) and HIV+METH+ (n=5) specimens (5 sections per subject), selected to display minimal confounders and no reported polysubstance use, were stained using immunohistochemistry to detect MRP8/14+ cells. Representative sections from each group (
Figure 10A) and the assessment of MRP8/14+ cell numbers per area (
Figure 10B) indicated that both HIV and Meth alone significantly increased 2 and 3-fold MRP8/14 cell numbers, respectively, compared to controls, while HIV+METH+ had a significant >4-fold increase in this subset. This suggests that Meth use is a factor in further increasing MRP8/14+ cells in HIV+ brains. Moreover, by RNAscope we observed an increase in detectable viral RNA in correlation with Meth use, in basal ganglia (
Figure 10C). This indicates that MRP8/14 may be a brain biomarker linked to viral replication in the context of Meth use in humans.
Figure 10.
– Detection of MRP8/14 in human brain specimens of HIV+ Meth users. Post-mortem formalin-fixed paraffin-embedded human brain specimen sections were obtained from NNTC. Immunohistochemistry was used to detect MRP8/14-positive cells in the basal ganglia. A) Representative sections from HIV-METH-, HIV+METH-, HIV-METH+ and HIV+METH+ brains. B) The average number of MRP8/14-positive events was counted in 25 fields per section, in 5 sections per subject, 5 subjects per group, using 40x magnification. C) RNAscope in situ hybridization was used to detect HIV RNA in the same sections. The pictures show a representative section from HIV+METH- and HIV+METH+ groups.
Figure 10.
– Detection of MRP8/14 in human brain specimens of HIV+ Meth users. Post-mortem formalin-fixed paraffin-embedded human brain specimen sections were obtained from NNTC. Immunohistochemistry was used to detect MRP8/14-positive cells in the basal ganglia. A) Representative sections from HIV-METH-, HIV+METH-, HIV-METH+ and HIV+METH+ brains. B) The average number of MRP8/14-positive events was counted in 25 fields per section, in 5 sections per subject, 5 subjects per group, using 40x magnification. C) RNAscope in situ hybridization was used to detect HIV RNA in the same sections. The pictures show a representative section from HIV+METH- and HIV+METH+ groups.
Discussion
This is study is the first report of a link between DA-induced signatures in innate immune cells and phenotypes of HIV latency reversal as biomarkers of CNS infection in the context of Meth use disorder.
The in vitro latent promonocyte cell line model has shown that DA has a moderate but significant effect on reversing HIV proviral latency while promoting the development of a subset of cells expressing pro-inflammatory markers. The cell subset that appears in correlation with higher p24 is CCR5high, CD38high and PDL1low, within the CD16high activated subset. We have previously identified cells with similar characteristics in the brains of SIV-infected rhesus macaques that were chronically treated with Meth, in correlation with a higher brain viral load4. The percentage of cells that upregulate these markers matches the calculated percentage of cells replicating virus following DA or DRD1/5-selective stimulation. However, it is not clear from these experiments, whether they correspond to the same cells. Yet, these are molecules with described consequences to inflammation and in HIV infection, as well as in neurological disorders. For instance, the increase in CCR5, which increased expression in peripheral cells has been described as a biomarker of Meth4,44,78,79.
CD38 is a multifunctional enzyme expressed on the cell surface that plays a major role in the hydrolysis of NAD, and consequently regulates Sirt1 enzymatic activity80,81, 21. We have shown evidence that latently infected myeloid cells exposed to DA experience the onset of a novel subset expressing high levels of CD38 and CCR5 in correlation with latency reversal. PDL1 on the other hand, indicates cell senescence82. A subset with low levels of PDL1 may have a long-living phenotype, indicating reservoir retention. We also show that a subset of CCR5high, CD38high and PDL1low, within the CD16high subset is increased in peripheral blood from PWH that are Meth users.
Transcriptionally, the stimulations with DA, DRD1/5 or DRD4 agonists, predominantly caused downregulation of genes rather than upregulation within the latently infected myeloid cells. Interestingly, we have previously observed that Meth administration to mice triggers strong suppression of gene expression in the brain, although it seems to be counteracted by the expression of Tat in the astrocytic compartment 83. Here, the in vitro data suggests that in the isolated latently infected compartment, DA as a major neurotransmitter in the brains of Meth users may contribute to retaining low gene transcription profiles while stimulating pathways that favor viral replication secondarily or as a result of inflammation. DA increased few markers involved in stimulus response and transcriptional regulation, which also have described effects on inflammatory outcomes. Two of these transcriptional signatures were associated with latency reversal via DRD1/5 in our model, MRP8 and MRP14, which form a heterodimer at the protein level73, also known as Calprotectin. Interestingly, the expression of MRP14 and the heterodimer have been previously shown to serve as a signature of a macrophage subset that is not infected but is enriched in the perivascular domain in HIV and SIV encephalitis84. We have validated the role of DA in the upregulation of the MRP8/14 heterodimer both secreted and on the cell surface via specific receptors, indicating that this may be a marker of CNS disorder in drug users, as well as a factor triggering replication in latent cells via binding to its receptor, RAGE on infected cells.
Overall, few reports link MRP8 and/or MRP14, alone or as heterodimer, with aggravated HIV pathogenesis 66,85, associated to CD4 T cells86, or associated with early monocyte migration84. A recent report has linked MRP8 (S100A8) expression to active HIV transcription and metabolism in cultured urethral macrophages, along to other inflammatory markers such as IL1R and MMP787. One report has found direct effects of MRP8 in HIV replication88. One report links the expression of the heterodimer in mucosa with viral shedding independent from plasma viral load89, suggesting localized effects. On the other hand, reports of no effects of complex on HIV infection, or suppression of p24 levels by subunits, have been performed in models where infected cells are differentiated and activated at baseline 90. These approaches differ significantly compared to us, as here we attempted to address effects on infected cells exhibiting a latent phenotype, relevant in the post-ART era. It also suggests that MRP8 and MRP14 may serve as latency reversal cytokines but may not further enhance ongoing replication, or may even be beneficial in the context of ART. These markers have been identified in other infections, such as Sars-Cov2, in association with high inflammatory levels90. In conditions of co-infection or co-morbidities, these cytokines could act differently. Moreover, it has been reported that MRP8/14 complex is increased in the cerebrospinal fluid of HIV+ subjects exhibiting co-infections91. It has never been examined in the context of SUD. It remains to be examined how DA affects actively replicating cells. In addition, the effects of MRP8/14 may differ by activation state in part due to binding to different receptors.
The receptor for MRP8/14-mediated effects on HIV replication in latent cells was identified as RAGE 92,93, a pattern recognition receptor that decreases in expression and function during aging, resulting in higher oxidative stress and with implications for vascular calcification, diabetes, dementia and other aging consequences94-96. Interestingly, RAGE has been described as a target for Alzheimer’s disease improvement and prevention of amyloid deposition in the CNS97-102, while increasing CCR5 expression and promoting leukocyte transmigration across the blood brain barrier103. Studies in TLR4 and RAGE knock-out models have suggested a major role of TLR4 as a receptor to MRP8/14 on myeloid cells 76. However, we did not find compelling evidence for TLR4 in latency reversal by MRP8/14 in our model. In HIV infection, low levels of soluble RAGE have been found in infected subjects with subclinical carotid atherosclerosis under ART 104, but also complex protective effects against cys-infection from myeloid dendritic cells due to CCR5 repression, in spite of the induction of pro-inflammatory cytokines105. In our model, where the ligands to RAGE are upregulated by DA and RAGE is involved in latency reversal, higher CCR5 expression is one of the detectable signatures in a subset of the cells, and during chronic Meth in vivo4,44,78. Perturbations to RAGE ligands such as monoamines, have been involved in CNS conditions, including aging and neurodegeneration, and impairment of antioxidant mechanisms106-110. An interesting study in flies has shown that feeding the animals with cocaine or Meth, as well as exposing them to DA, increases the levels of advanced glycation end products, RAGE ligands. Moreover, mutant flies lacking the dopaminergic transporter or DRD1 have shown accumulation of these ligands111, further suggesting the role of RAGE in DA-induced perturbations. Here we show that DA exposure promotes an increase of a myeloid-related complex MRP8/14 that serves as RAGE ligand in latently infected human myeloid cells. The results suggest that these molecules have a biomarker value and play a role in HIV replication linked to inflammation. We have also shown that cells expressing these molecules are enriched in the brains of HIV+ individuals and even further in Meth users, suggesting a mechanism for DA-mediated increase of viral load in the brains of HIV+ individuals with substance use disorders. Relevant to the brain, MRP8 has been also previously implicated in neurological outcomes such as depression68, inflammation in the CNS86,112-115 and elsewhere72,116-119, as mediator and biomarker.
The differential effects of DRDs on the onset of subset populations in correlation with levels of p24 were intriguing, particularly the effects of DRD4 versus DRD1/5. We show evidence that the RAGE ligand MRP8/14 heterodimer is increased in vitro by DA and in vivo by Meth, and that RAGE surface expression on latent myeloid cells is decreased by DRD4, by means that may actively maintain latency. RAGE is also expressed at lower levels on monocytes from HIV+ Meth users with undetectable viral load. Interestingly, DRD4 is a gene with an unusually high degree of polymorphism, with some variants increasing the risk of substance use disorders120-123. The implications of this polymorphism for immune responses in the context of substance use or neurological disorders have not been examined.
Although it is known that immune functions interfere with neuronal functions and vice-versa, specific DRD subtypes, immune activation and neurological disease have been rarely linked. It has been shown that DA signals via DRD4 to promote Th2 responses, with implications for allergic reactions124. In microglia, a decrease in the inflammatory activation has been described in hyperdopaminergic conditions, attributed to a decrease in nitric oxide125 and of phagocytic capacity126, potentially mediated by DRD1 and DRD2127. It has been also suggested that DRD1 signaling in microglia induces degradation of the NLRP3 inflammasome in a cAMP-dependent fashion in vivo128. Here, we show that downregulation of genes in latent innate immune cells is more significant than upregulation, both via DRD1/5 and via DRD4, with distinct and common signatures. Yet, we also show that among few upregulated genes, pro-inflammatory ones are expressed. Similarly, the administration of Meth in mouse models of neuroHIV has led to the observation of drastic suppression of gene expression, counteracted by expression of the HIV accessory protein Tat and with focal increase of specific inflammatory pathways 83.
The combined expression of the markers increased by DA via DRD1/5 in a subset of cells, CCR5, CD38 and PDL1, along with MRP8/14 and RAGE, have indicated strong potential to distinguish PWH that are Meth users with detectable and undetectable CSF viral load. We did not find a correlation with plasma viral load, further suggesting the CNS-localized effect of the proposed mechanism of latency reversal involving upregulation of MRP8/14 and its binding to RAGE, triggered by high levels of DA in the context of Meth use. The analysis of males, which was due to cohort characteristics and to increase the power with an n=25, is a limitation. Yet, the study addresses the value of biomarkers considering demographic trends and a gender disparity in the population of individuals with SUD and HIV129.
Interactions between HIV and neurotransmitters induced in the brain by addictive stimulant substances must be better understood for addressing consequences, including the aggravated damage to areas bearing dopaminergic projections, such as the basal ganglia 130,131, and the neuro-pathologic effects of glial cells and inflammation 132. Specific signatures of virus-infected cells in the context of a hyperdopaminergic micro-environment may explain variabilities of inflammatory outcomes in the context of substance use and provide mechanisms to explain changes in replication as well as aggravation that are found in Meth users.