- 1Departments of Neurological Sciences, University of Nebraska Medical Center, Omaha, NE, United States
- 2Genetics, Cell Biology and Anatomy, University of Nebraska Medical Center, Omaha, NE, United States
- 3Pharmacology and Experimental Neuroscience, University of Nebraska Medical Center, Omaha, NE, United States
- 4The Genomics Core Facility, University of Nebraska Medical Center, Omaha, NE, United States
- 5The Antiviral Pharmacology Laboratory, University of Nebraska Medical Center, Omaha, NE, United States
The twin pandemics of opioid abuse and HIV infection can have devastating effects on physiological systems, including on the brain. Our previous work found that morphine increased the viral reservoir in the brains of treated SIV-infected macaques. In this study, we investigated the interaction of morphine and SIV to identify novel host-specific targets using a multimodal approach. We probed systemic parameters and performed single-cell examination of the targets for infection in the brain, microglia and macrophages. Morphine treatment created an immunosuppressive environment, blunting initial responses to infection, which persisted during antiretroviral treatment. Antiretroviral drug concentrations and penetration into the cerebrospinal fluid and brain were unchanged by morphine treatment. Interestingly, the transcriptional signature of both microglia and brain macrophages was transformed to one of a neurodegenerative phenotype. Notably, the expression of osteopontin, a pleiotropic cytokine, was significantly elevated in microglia. This was especially notable in the white matter, which is also dually affected by HIV and opioids. Increased osteopontin expression was linked to numerous HIV neuropathogenic mechanisms, including those that can maintain a viral reservoir. The opioid morphine is detrimental to SIV/HIV infection, especially in the brain.
Introduction
The clinical outlook for people with HIV (PWH) has improved dramatically as increasingly effective combination antiretroviral therapy (cART) has become available (1–4). Since untreated HIV infection leads to the acquired immunodeficiency syndrome (AIDS), a focus on the immune abnormalities has been prominent. However, the central nervous system could also be adversely affected. Before the onset of successful treatment options, 8-15% of PWH in the United States (US) developed a severe cognitive disorder, HIV-associated dementia (HAD) (5, 6). In the post-cART era HAD has been reduced to 1-2% in the US and 5% globally (7, 8). Paradoxically, increased survival rates resulting from therapy have also resulted in an undesirable increase in the prevalence of minor neurocognitive impairments, collectively referred to as HIV-associated neurocognitive disorders (HAND), also known as neuroHIV (7, 8). Unfortunately, cART does not completely shield the brain from the effects of HIV-1 infection.
HIV, as well as the analogous nonhuman primate (NHP) virus simian immunodeficiency virus (SIV), enters the brain early after infection (9–13), and can set up a lifelong presence in the central nervous system (CNS) (14, 15). Many widely used antiretroviral drugs poorly penetrate the CNS, perhaps facilitating pathogenic effects of infection and the maintenance of a viral reservoir (16). HIV infection remains a significant problem: at the end of 2018, there were over one million people diagnosed with HIV infection in the US and an estimated over 160,000 people whose infections had not been diagnosed (17).
Over two million people in the US have opioid use disorder (OUD), a significant risk for HIV infection and transmission; in addition, PWH are frequently prescribed opioids for pain management (18). Opioid abuse and HIV have been described as dual linked health crises (19), and drug abuse could likely be one of the leading causes of HIV transmission in the United States. Injection drug use was found to be responsible, solely or in conjunction with other causes, for 11% of HIV infections in 2020 (20), and other cases can arise through risky sexual behavior due to drug use (21). While HIV infection has transformed from a death sentence to a manageable chronic condition, the comorbid condition of drug use, such as opioid abuse, in infected individuals is on the rise, and may lead in turn to increased neurologic and cognitive deficits in PWH (22). The CNS may be uniquely susceptible to the combined effects of opioid abuse and HIV. Among the various abused drugs, the use of opioids is highly prevalent in PWH, both from illicit drug use and the management of HIV-associated neuropathic pain. It is estimated that one-quarter to one-half of PWH are prescribed opioids, and PWH are prescribed higher doses than those without HIV (23). Intriguingly, brain regions that are targets for morphine (expressing high numbers of μ-opioid receptors), are also the regions with increased viral loads and predilection for HIV infection (24–26). Thus, the interplay of HIV and opioids raises concerns regarding the effect of opioids on HIV’s impact on the brain. Identifying CNS molecular aberrations underlying HIV and opioid abuse thus remains an unmet need in the field.
Regarding neuroHIV, the effects of substance abuse can vary (27, 28). Studies from both before and after the introduction of treatment have reported increased neuropathogenesis associated with opioid use (29, 30). Still, studies on opioid users as well as results on the effects of opioids on HIV studied through in vivo and in vitro systems have yielded varied results on whether opioids increase neuroHIV or factors attributed to HIV neuropathogenesis; these have been thoroughly reviewed recently (18, 22, 31). The different findings could be attributed to the pleiotropic effects of opioids, such as in the balance of opioid-enhanced infection by the virus versus opioid-induced dampening of inflammatory lesions. In addition, the conflicting findings in humans could be attributed to high variability in drug usage patterns, use of different opioids, multidrug use, differences in nutritional and socio-economic status, variable access to medical care, and potential inaccuracy in drug use data collection due to overdependence on self-reporting. In the current era issues such as access for treatment in opioid users due to economic and social issues, difficulties in treatment compliance due to drug use and its consequences such as incarceration, and comorbid conditions due to drug use are additional complicating issues clouding the ability to determine the effects of opioids in PWH.
NHP infection by SIV remains the best model for the study of HIV pathogenesis and treatment, including neuropathogenesis (32–37). NHPs are also an excellent animal model for studying the effects of drugs of abuse and their interactions with SIV infection, especially for effects on the brain (38–42). Such animal experiments allow better controlled conditions than observational studies in humans. However, studies on the effects of opioids during experimental infection of NHP by SIV in the absence of treatment have yielded inconsistent findings (39, 41, 43–46).
In SIV-infected NHP, as in people, the viral targets in the brain are primarily brain macrophages and microglia (47, 48). These key cells support productive HIV and SIV replication in the CNS and can serve as virus reservoirs in the brain (48–52). Infection of macrophages and microglia and the subsequent inflammation associated with infection can alter their function and phenotype, affecting other glia as well as neurons. In parallel, opioids have also been shown to impact macrophage and microglial functioning and phenotype (53–56).
We recently reported that in monkeys chronically administered the prototypic opioid morphine, followed by SIV infection and subsequent suppression of viremia below the limit of detection by cART, the size of the viral reservoir in the brain was increased relative to that in animals given saline instead of morphine (57). In the current study, we investigated longitudinal changes in response to SIV infection of rhesus monkeys in the presence of morphine compared to animals administered saline before and after cART treatment. Using single-cell and single-nucleus RNA sequencing (scRNA-seq and snRNA-seq), we examined the effect of morphine on brain-resident macrophages and microglia in these SIV-infected, cART-suppressed animals to better understand the factors responsible for the morphine-induced increased CNS reservoir and the role of opioids in neuropathogenesis.
Methods
Animals
The animals used in this study have been previously described (57). Briefly, the rhesus macaques tested negative for a panel of viral and bacterial pathogens before their use in the study. Animals in the morphine group received intramuscular injections of morphine, which were ramped-up over a 2-week period to a final dose of 6 mg/kg morphine administered twice daily on weekdays, and once daily on weekends. This dose was maintained for 7 weeks. Animals in the saline group received saline injections on the same schedule. Animals were then inoculated with SIVmac251. Five-weeks post-inoculation, all animals received a once-daily subcutaneous injection of 1 mL/kg of body weight containing 40 mg/mL emtricitabine (FTC), 20 mg/mL tenfovr (TFV) and 2.5 mg/mL dolutegravir (DTG), continued until the time of necropsy. At necropsy, blood-borne cells were cleared by perfusion with PBS containing 1U/mL heparin. Microglia/macrophage-enriched cellular isolation from the brain was performed using our previously described procedure of physical and enzymatic disassociation followed by Ficoll gradient purification (10, 58), and then cryopreserved in the presence of 10% DMSO/90% fetal bovine serum and stored over liquid nitrogen. Brain regions were dissected by a board-certified pathologist, frozen, and stored at -80° C.
Blood cell subsets, albumin, and flow cytometry
Complete blood counts (CBC) and chemistry/metabolic panels were performed in the CLIA-certified Pathology labs of Nebraska Medicine (Omaha, NE). Overall cell numbers and proportions of the general classes of cells were reported on the CBC; the albumin concentration was determined on the chemistry/metabolic panel.
Flow cytometry evaluation was completed on peripheral whole blood and assessed prior to viral inoculation and longitudinally throughout the study. Cells and staining reagents were kept on ice, centrifugation was conducted at 4° C and 3500 × g. Cells were first stained with UV blue Live dead assay (Invitrogen, Carlsbad, CA), washed, blocked using flow cytometry buffer (eBioscience, Waltham, MA), pelleted, and surface antigens stained with a cocktail containing brilliant stain buffer (BD Biosciences, San Jose, CA) and fluorochrome-conjugated monoclonal antibodies against CD11b-BV605 (Biolegend, San Diego, CA), CD14-ECD (Beckman Coulter, Brea CA), CD3-BV421, CD4-BV786, CD8-PeCy7, CD16-APCH7, and CD20-BV711 (all from BD Bioscience). Stained whole blood was washed with flow cytometry buffer, pelleted and red blood cells lysed with red blood cell lysis buffer (Sigma Aldrich, St. Louis, MO). Lysed samples were washed in flow cytometry buffer, pelleted, and fixed using 1% paraformaldehyde in flow cytometry buffer. Cells were stored overnight at 4° C before acquisition. All events were acquired using the BD Fortessa X450, and data were analyzed using BD FACSDiva v9.0 or FlowJo for Mac version 10.6 (BD Biosciences).
The lymphocyte percentage was determined as proportion of white blood cells measured on the day of the CBC draw. Lymphocyte percentage was multiplied by the white blood cell count to obtain absolute lymphocyte count per microliter, and further multiplied by the flow cytometric determination of the CD3+CD4+ and CD3+CD8+ proportions within the lymphocyte FSC/SSC gate on live cells to obtain absolute number of CD4+ and CD8+ T-cells, respectively, per microliter of blood.
Plasma biomarkers
Plasma, collected in EDTA-anticoagulated blood, was separated within an hour of collection and stored frozen at -80° C. Samples were sent to Myriad RBM (Austin, TX), a CLIA-certified lab, for multiplexed biomarker assay using the Human v2.0 Multianalyte profile (MAP) immunoassay. We previously reported the normal range for the 46 analytes that can be detected in rhesus macaques (59).
Antiretroviral quantitation
Concentrations of TFV, FTC, and DTG in plasma and CSF, and DTG and the intracellular triphosphates (TP), tenofovir-diphosphate (TFV-DP) and emtricitabine-triphosphate (FTC-TP) in tissue homogenates, were quantified by liquid chromatography with tandem mass spectrometry with previously described validated methods (60–62).
ScRNA-seq/snRNA-seq
From the 11 SIV-infected, cART-treated animals (6 in the morphine group, 5 in the saline group), 4 from each group were used for the single cell/nucleus analyses. For scRNA-seq of microglia, the saline group consisted of samples from animals 74T (referred to as CI47 in (57)), 77T (CK40), 82T (CI04), and 83T (CK03), and the morphine group samples from animals 73T (CK22), 75T (CK48), 78T (CI90), and 80T (CH89). The scRNA-seq data for the three of the members of the saline group, 74T, 82T, and 83T, were reported previously (63). For the snRNA-seq several animal samples were substituted within the groups due to the availability or suitability of samples for processing and analysis. For the frontal lobe white matter, 81T (CK07) was substituted for 82T for the saline group and 76T (CK31) for 80T in the morphine group. For the caudate those same substitutions were made, in addition morphine 79T (CI83) was substituted for 75T in the morphine group.
ScRNA-seq procedures were completed using our previously described methods (58, 64). Briefly, cryopreserved microglia were thawed and treated with 1% DNase (Sigma Aldrich), washed, stained with UV-blue live/dead assay (Invitrogen, Waltham, MA). Post stain, microglia were washed, resuspended in MACS separation buffer with 0.1% bovine serum albumin, isolates were then enriched with nonhuman primate anti-CD11b microbeads (Miltenyi, Gladbach, Germany), followed by washing and positive selection on MACS Separator LS columns. Recovered cells were stained anti-mouse/human CD11b antibody clone M1/70 (BV605; Biolegend, San Diego, CA), washed, and sorted based on size, singlets, live, CD11b+ events using an Aria2 flow cytometer (BD Biosciences, San Diego, CA).
For snRNA-seq, caudate tissue was chopped fine, transferred to a Dounce homogenizer with 300 µL of NP40 lysis buffer (consisting of 10 mM Tris-HCL (pH 7.4), 10 mM NaCl, 3 mM MgCl2, 1% nonidet P40 substitute, 1 mM DTT, 1 U/µL RNase inhibitor; all from Sigma Aldrich, St. Louis, MO). Samples were homogenized, rested on wet ice, washed in 1% BSA in dPBS, rested on wet ice, and filtered through 70 µm filter. Samples were volume adjusted to 1 mL and stained with 7-AAD staining solution (Sigma Aldrich). Nuclear events were captured based on size, singlets, and nuclear based on 7-AAD stain, using an Aria2 flow cytometer.
A similar method was used to isolate nuclei from white matter samples, except after filtration samples were adjusted to a total volume of 2 mL of 1% BSA in dPBS. Samples were divided into two 2 mL tubes. In each tube, 1 mL of filtered sample was mixed with 300 µL of debris removal solution (Miltenyi) overlaid with 1 mL 1% BSA in dPBS. Samples were centrifuged at 3000 x G for 10 minutes. Subsequent upper layers were discarded and lower layer, containing the nuclei, were washed in ice cold 1% BSA in dPBS. Nuclei isolates were washed again, filtered through 70 µm filter and final volume was adjusted to 1 mL. Samples were then stained with 7-AAD and sorted on an Aria 2 flow cytometer.
Post-processing, isolates were concentrated to approximately 1,000 cells or nuclei/μL and assessed by trypan blue for viability and concentration. Based on 10x Genomics parameters targeting 8000 cells, the calculated volume of cells or nuclei was loaded onto the 10x Genomics (Pleasanton, CA) Chromium GEM Chip and placed into Chromium Controller for cell capturing and library preparation. We used the 10x Genomics Single Cell 3’ GEM, Library, and Gel bead kit v3.1. The prepared libraries were sequenced using a NovaSeq6000 sequencer (Illumina, San Diego, CA) using Kit S1-100 Version 1.5 for scRNA-seq and Kit S2-100 Version 1.5 for snRNA-seq.
Bioinformatics
We generated FASTQ files and count matrices using 10x Genomics cellranger pipeline version 3.1.0. Specifically, we used cellranger mkfastq to convert the raw base calls from the Illumina NovaSeq6000 into FASTQ files. Then we used cellranger count to map the reads to the customized combination reference genome of Macaca mulatta (Mmul10) and the SIV genome divided into 5 regions (based on NCBI reference sequences M22262.1), along with filtering, barcode counting, and unique molecular identifier (UMI) counting (64). As a result, a gene-cell matrix is generated for each sample. DoubletFinder was used to remove doublets on each of the matrices (65).
For in-depth analysis, we used Partek (St. Louis. MO) Flow version 10.0.22. For scRNA-seq of microglia, we used single-cell QA/QC component to filter out cells with the following attributes: total UMI <400 or >20,000, gene count per cell <300 or >5000, and mitochondrial UMI proportion > 15%. Lymphocytes were removed (based on canonical T-cell and B-cell gene markers), and an average of 5723 cells per animal was identified for analysis. The remaining genes were then filtered by a rhesus monkey to human mapping list (64). Non-protein coding RNAs, mitochondria-encoded genes, ribosomal protein encoding genes, and any gene that had a maximum expression value ≤1 was removed, resulting in 12,308 genes for analysis.
For snRNA-seq we used single-cell QA/QC component to filter out cells with the following attributes: total UMI <1000 or >100,000, gene count per cell <700 or >10,000, and mitochondrial UMI proportion > 1%. An average of 5616 nuclei from frontal lobe white matter and 5906 nuclei from the caudate were identified for analysis from each animal. The remaining genes were then filtered by a rhesus monkey to human mapping list (64), and any gene with a maximum expression value ≤1 was removed, resulting in 14,143 and 14,379 genes for analysis for the frontal lobe white matter and caudate nucleus, respectively.
Normalization was performed using scTransform (66), followed by principal component analysis (PCA) and then Harmony to eliminate batch effects (67). Uniform manifold approximation and projection (UMAP) (68) and graph-based clustering were then performed, the latter using Louvain clustering (69) with a resolution of 0.3. ANOVA was used for differential expression analysis based on the expression data normalized by log2(counts per million plus 1). Differentially-expressed genes (DEGs) lists were defined as genes with a false discovery rate of <0.01, a fold change >|1.5|and a minimal expression (least square mean) of at least 50 in one of the conditions or clusters.
Gene set enrichment analysis (GSEA) (version 4.1.0) was performed on gene sets for disease-associated microglia (DAM), activation-response microglia (ARM), (both sets taken from Supplemental Table S1 in (70)), and neurodegeneration-related brain myeloid cells (from Supplemental Data S4 in (71)), converting the genes to rhesus gene symbols and eliminating any ribosomal protein genes. To perform GSEA we then used the pre-ranked gene lists from each of the clusters for those genes with an FDR <0.01 on differences between morphine and saline for each cluster, with 1000 permutations to calculate the p-value (72). Ingenuity pathway analysis (IPA) (73) was performed to generate canonical pathways based on the DEGs between morphine and saline for each cluster.
Data availability
The snRNA-seq data have been deposited in NCBI GEO, accession # GSE209606. The scRNA-seq data for animals 73T, 75T, 77T, 78T, and 80T have also been deposited in accession # GSE209606. The scRNA-seq data for 74T, 82T, and 83T were previously deposited, accession # GSE195574.
Statistics
Statistical analyses for studies other than those described in the bioinformatic section above were performed using Prism (version 9.4.0, GraphPad Software, San Diego, CA). Specific tests are indicated in each figure legend.
Ethical approval
No new animals were used in these studies, experiments were carried out on archived material from a prior study (57). In that study, macaques were housed in compliance with the Animal Welfare Act and the Guide for the Care and Use of Laboratory Animals in the nonhuman primate facilities of the Department of Comparative Medicine, UNMC, which has been accredited by the AAALAC international. The UNMC Institutional Animal Care and Use Committee reviewed and approved that study under protocol 16-073-FC.
Results
Acute cellular response to SIV is suppressed by morphine
In our recent study, we infected morphine-dependent (n=6) and saline-administered control (n=5) rhesus monkeys with SIV, followed by treatment with cART to suppress viremia, modeling the acquisition of HIV infection in those who use opioids. After 30 weeks (or more) of treatment with antiretrovirals, viremia was durably suppressed in the plasma and cerebrospinal fluid (CSF), and analysis of the brain viral reservoir revealed that it was higher in CD11b+ microglia/macrophages in morphine-dependent animals than in those given saline as a control (57). Numerous factors could drive this effect. To address these in the context of morphine, we first examined two general health measures and several aspects of blood cell composition in the two groups throughout infection (Figure 1). The viral loads in plasma and CSF are also shown for reference. Morphine-treated animals had a significantly lower weight gain than the saline-group, consistent with the known effects of morphine on appetite. Albumin, a general marker of nutrition, remained unchanged.
Figure 1 Longitudinal Assessments. Viral loads, general health, and blood cell composition in the saline (blue) and morphine (red) groups. The points represent the mean and the error bars the standard error of the mean. Viral inoculation is denoted by black arrow at week 0; cART was begun at week 5 and continued throughout the rest of the experiment. Biofluid viral loads were previously reported (57). The weight and albumin data begin at week -2, when the morphine group animals were on stable maintenance doses of morphine. The blood cell composition includes a pre-morphine value (and the equivalent time in the saline group). Statistical testing was done by 2-way ANOVA and the p-value indicated above for significant differences between the groups. When the ANOVA was significant, testing on the individual time points was performed using Sidak’s multiple comparison test. * = <0.05, ** = < 0.01, *** = < 0.001.
Striking findings were found in blood cells. The total number of white blood cells and the lymphocyte subset were significantly lower in the morphine group. All lymphocyte subsets examined: CD4+ T cells, CD8+ T cells, NK cells and B cells were lower in the morphine group. Strikingly the increase in CD8+ T cells and NK cells at the acute infection time point (2 weeks p.i.) in the saline group was significantly reduced in the morphine group. Monocytes showed the same pattern, lower in the morphine group, and the increase found in the saline group at 2 weeks p.i. largely ablated. Platelets, while dropping in the acute and post-acute infection stage, with recovery after cART administration, did not differ between the two groups. Thus the acute cellular immune reaction to infection, in terms of increases in lymphoid and myeloid cells, was significantly suppressed in response to morphine.
Chronic morphine creates an immunosuppressive environment
Next, we addressed the expression of several blood biomarkers of various physiological processes, including inflammation. Plasma from the six morphine-administered animals was assessed pre- and post-morphine ramp-up and maintenance administration (2-week ramp-up period and 7-week maintenance), before SIV infection. Interestingly, several molecules were altered by morphine (Figure 2, left graft for each analyte). Adiponectin, a hormone made by fat cells, functions to sensitize cells to insulin as well as anti-inflammatory effects (74), was significantly increased by morphine. In keeping with this finding on an anti-inflammatory molecule, a number of inflammatory markers were decreased following morphine treatment: interleukin 16 (IL-16), a proinflammatory cytokine that is a chemoattractant for a number of different classes of immune cells (75); matrix metalloprotein 2 (MMP-2), one of a family of MMPs that serve to degrade the extracellular matrix to enable cellular inflammation (76); tumor necrosis factor receptor 2 (TNFR2), one of the two TNF receptors, and its soluble form used as a biomarker for TNF activity in inflammation (77); and vascular cell adhesion molecule 1 (VCAM-1), an endothelial adhesion molecule for immune cells, and its soluble form used a biomarker for inflammation and endothelial activation (78).
Figure 2 Plasma biomarkers. A panel of plasma proteins was assessed before and after morphine administration in the morphine group, and then throughout infection in both the morphine (red) and saline (blue) groups. The left panel for each analyte contains the pre- and post-morphine (before SIV inoculation); statistical testing was performed by paired t-test. The right panels show values starting on the day of inoculation (week 0), and then week 2 (acute), week 5 (immediately before cART administration began), and then with viremia suppressed by cART, week 22 and week 35. The points represent the mean and the error bars the standard error of the mean. Statistical testing was done by 2-way ANOVA and reported for significant differences between the groups. When the ANOVA was significant, testing on the individual time points were performed using Sidak’s multiple comparison test. * = <0.05, ** = < 0.01, *** = < 0.001. LLOQ = lower limit of quantitation, shown as a green line. One-half the LLOQ value was used for nonquantifiable values <LLOQ.
Following this baseline period of morphine administration, both the morphine and saline groups of animals were inoculated with SIV. In addition to the pre-inoculation time point (week 0), biomarkers were examined during the acute phase of infection at the peak of viremia (2 weeks post-infection, p.i.), when the acute viremia was lowered by the host immune response (5 weeks p.i.), at two points following virus suppression by cART (22 and 35 weeks p.i.) (Figure 2, right graph for each analyte). For the longitudinal measurements, each of the above analytes was significantly different, with adiponectin increased, and MMP-2, TNFR2, and VCAM-1 decreased. Strikingly, several additional molecules were also altered over the time course and showed specific significant differences during the acute 2-week p.i. time point. C-reactive protein (CRP), a general marker for inflammation (79), increased dramatically in saline-administered animals, whereas the levels did not change over baseline in the morphine-administered animals. Ferritin, an iron-sequestering protein, is known to be elevated in inflammation due to cellular damage (80), and increased significantly in saline compared with morphine-administered animals. The interleukin 1 receptor antagonist (IL-1ra), produced in response to inflammation to help control proinflammatory IL-1 effects (81), was also elevated to a greater extent in the saline group. Interleukin 18 (IL-18), a pleiotropic proinflammatory cytokine in the IL-1 family (81), was also elevated acutely in the morphine group.
The expression of several additional molecules was found to be altered in the morphine-treated group compared with the saline-treated group over the course of infection. Complement factor 3 (C3), which is elevated in inflammation (82), showed a slight decrease in the morphine group, reaching significance at week 5 p.i. Macrophage inflammatory protein-1 beta (MIP-1β), a pro-inflammatory chemokine for a number of lymphocyte subsets and inhibitor of HIV and SIV infection through CCR5 (83), also showed slightly decreased levels in the morphine-treated group. In contrast, C-C motif chemokine ligand 18 (CCL18) was increased over the time course of infection by morphine, particularly in the chronic virus-suppressed time points. CCL18 has been found to have both immune activating and immune suppressing functions (84). A list of the analytes measured, and their level of expression is given in Supplemental Table 1.
In summary, morphine generally suppressed inflammatory markers before and after infection. These effects were magnified at the acute 2-week p.i. time point, a time point when the immune system is highly activated in response to the infection and high level of viremia. This decreased inflammatory reaction is consistent with the decreased cellular response above in the morphine group.
CNS and peripheral antiretroviral distribution and concentration are unchanged by morphine
While the morphine immunosuppression could help explain the increased SIV myeloid reservoir present in the brain and untoward effects of opioids on HIV pathogenesis and neuroHIV, it is likely that modulation of antiretroviral (ARV) levels in biofluids and tissues by morphine can affect the reservoir size. We thus sought to measure the concentrations of ARVs in the plasma and CSF as well as tissues. At necropsy, animals were given their daily subcutaneous injection of cART, followed 90 minutes later by bleeding and CSF sampling, with subsequent tissue harvest 20-30 minutes later. No difference was found between the ARV drug concentrations in the saline or morphine groups in any of the measured parameters (Figure 3A). Some differences, however, were found among the different brain regions (Figure 3B). For example, dolutegravir (DTG) was lower in the white matter than the putamen or pons. The active triphosphate of tenofovir (TFV-diphosphate, TFV-DP) was lower in the frontal cortex than the hippocampus, putamen, pons and white matter, and the putamen lower than pons and white matter. The active triphosphate of emtricitabine (FTC-TP) was lower in the frontal cortex and putamen than in the pons. When significant differences were found, the pons had higher levels for all three ARVs within the brain. Still, the levels in the brain were consistently lower than those in the rest of the body. Interestingly, examining CSF concentrations, FTC had the highest penetration from plasma to CSF, whereas for the active drug, when comparing spleen to brain concentrations, FTC-TP was the lowest of the three drugs (Figure 3C), consistent with our findings in rodents (85, 86). Although in the biofluids the CSF drug concentrations are much lower than those in plasma, the CSF drug concentrations exceeded the protein-free (as CSF has greatly reduced protein compared to plasma) 50% inhibitory concentration (IC50, as measured for HIV) for TFV (1.98-fold greater), FTC (8.59-fold greater) and the IC90 for DTG (13.3-fold greater). However the intracellular tissue concentrations of the active drugs were below the IC50/IC90 in the brain for the drugs, with higher values, many exceeding the IC50/IC90, in the lymphoid organs and gut (Supplemental Figure 1). Overall, there is no evidence that morphine-related differences in ARV brain concentrations affect the CNS viral reservoir or aspects of CNS function.
Figure 3 Antiretroviral drug measurements in biofluids and tissues at necropsy. (A) Concentrations of dolutegravir (DTG), tenofovir (TFV) and emtricitabine (FTC) in the plasma and CSF, 90 minutes after subcutaneous administration of the drugs. (B) Tissue concentrations (with brain regions indicated) of DTG and the active diphosphate (DP) metabolite of TFV (TFV-DP) and triphosphate (TP) metabolite of FTC (FTC-TP). Amounts of these entities were normalized to the amounts found in the spleen of each animal. Statistical testing was done by 2-way ANOVA for the brain region only. When the ANOVA was significant, testing on the individual regions was performed using Sidak’s multiple comparison test. * = <0.05, ** = < 0.01, *** = < 0.001, **** = <0.0001. (C) Ratio of the concentrations found in CSF relative to spleen (left) and of the average of the 5 brain regions relative to spleen (right).
Chronic morphine does not change the pattern of CNS myeloid cell subsets
Since SIV and HIV target the myeloid cells in the brain, and such cells can serve as a viral reservoir during treatment as well as drive neuropathogenesis, we next characterized the myeloid cells (microglia and macrophages) purified from the brain of four animals in each group by single-cell RNA sequencing (scRNA-seq). On average we examined ~5,700 such cells from each animal. No expression of SIV transcripts was identified in the cells. Graph-based clustering revealed the cells could be divided into six clusters (Figure 4A). Morphine did not alter the proportion of cells in the clusters, as the distribution of cells from each condition was similar between the clusters (Figure 4B). Using a panel of myeloid gene markers, based on gene expression profiles, clusters 1-4 and 6 most closely resemble microglia, with higher C-X3-C motif chemokine receptor 1 (CX3CR1), G protein-coupled receptor 34 (GPR34), and purinergic receptor P2Y12 (P2RY12). On the other hand, cluster 5 consists of macrophages, with higher levels of CD74 (the major histocompatibility class II invariant chain), CD163, major histocompatibility complex, class II, DR alpha (MAMU-DRA), as well as S100 calcium binding protein A4 (S100A4) expression (Figure 4C).
Figure 4 ScRNA-seq analysis of brain macrophages and microglia from the saline and morphine groups of animals. (A) Uniform Manifold Approximation and Projection (UMAP) representation of the cells, grouped in 6 clusters and differentially colored as determined by graph-based clustering. (B) Pie charts of the proportion of cells in each cluster for the saline and morphine groups. (C) Bubble plot of selected genes differentiating the clusters. The color indicates the least-square mean expression (log2), and the size of the circle indicates the percentage of cells in which gene expression was detected.
Expression values for the cells in the clusters were examined for differentially expressed genes (DEGs) versus the other clusters, some of the genes are illustrated in Figure 4C. Clusters 1 and 2 did not differ much in DEGs compared to the other clusters, and showed some differences with each other, with complement C1q B chain (C1QB) and beta-2-microglobulin (B2M) higher in cluster 1, and activating transcription factor 3 (ATF3) and FosB proto-oncogene, AP-1 transcription factor subunit (FOSB) higher in cluster 2. Cluster 3 was notable for higher levels of numerous genes, including increases in apolipoprotein C1 (APOC1), major histocompatibility complex, class II, DR alpha (MAMU-DRA), and secreted phosphoprotein 1 (SPP1). Cluster 4 had increased expression of C-C motif chemokine 3-like and 4 (CCL3L and CCL4), whereas in cluster 6 Rho family GTPase 3 (RND3) was greatly increased. In general clusters 1 and 2 expressed a homeostatic microglia pattern (87–91). The differences between cluster 1 and 2 may be stress-related, as ATF3 and FOSB are members of the AP-1 transcription factor family, induced in acute stress (92). ATF3 has been proposed to be induced as a negative-feedback loop to regulate macrophage inflammatory cytokine production (93). APOC1 and SPP1, prominent in cluster 3, are two of the genes characterizing the disease-associated microglia (DAM) subtype of microglia found in animal models of Alzheimer’s disease and other conditions (71, 88). Regarding cluster 4, CCL4 (as well as CCL3) expression characterizes a small population of microglia that increase with aging in mice (87). The high expression of RND3 is intriguing as it has been proposed to induce a neuroprotective phenotype in microglia (94). A complete list of DEGs between the clusters is presented in Supplemental Table 2.
Chronic morphine alters CNS myeloid cell gene expression leading to a neurodegeneration-related phenotype
We next examined the clusters for DEGs in the context of morphine administration. The microglial clusters showed appreciable similarity in the DEGs. In the five microglia clusters (clusters 1-4 and 6), 17 genes were upregulated in common. Microglia cluster 6 failed to exhibit morphine downregulated genes. Excluding cluster 6, there were 33 upregulated, and 22 downregulated, genes that were common in microglia clusters 1-4 (Supplemental Figure 2).
We used gene set enrichment analysis (GSEA) to assess whether morphine altered the phenotype of the clusters in a characteristic manner. We used curated gene lists (Supplemental Table 4) for DAM (88) and activation-response microglia (ARM) (89), as well as neurodegeneration-related brain myeloid cells (71). As shown in Table 1, morphine induced a neurodegeneration-related myeloid cell phenotype in all the microglia clusters except the smallest cluster, cluster 6, and a DAM and ARM phenotype in the largest cluster, cluster 1. Macrophage cluster 5 also took on a neurodegeneration-related phenotype in response to morphine.
Table 1 Gene set enrichment analysis. The enrichment score (ES) is the degree to which a set is overrepresented comparing the morphine with the saline group.
The DEGs from the morphine versus saline condition in each cluster were then used for Ingenuity Pathway Analysis (IPA). We focused on the canonical pathways, and the results are presented in Supplemental Table 5; selected pathways are shown in Table 2. Interestingly the phagosome formation pathway was increased in 3 of the microglia clusters, and several signaling pathways were increased in the microglia clusters, including those of neuroinflammation in cluster 3, TREM1 in cluster 4, and G-protein coupled receptor in cluster 6. Signaling pathways were also enriched in macrophage cluster 5, including those involved in leukocyte extravasation and neuroinflammation.
Several morphine-upregulated DEGs are illustrated in Figure 5. In the microglia clusters 1-4 and 6, transcripts for immediate early response 3 (IER3), a feedback inhibitor of NF-κB in the inflammatory response (95), were increased, as was the mRNA for receptor activity-modifying protein 1 (RAMP1), a component of the calcitonin gene-related peptide receptor complex, which has microglia-activating activity (96). Also increased was the sodium bicarbonate cotransporter 3 (SLC4A7), critical for phagosome acidification (97), and SPP1 (especially prominent in cluster 3), a pleiotropic cytokine with numerous links to CNS diseases (98). There was no such change in the expression of these genes in macrophage cluster 5. Instead, several genes were upregulated in the morphine group, including C-type lectin domain family 7 member A (CLEC7A), a myeloid-predominant pattern-recognition receptor implicated in protection from Alzheimer’s disease (99, 100), type II major histocompatibility genes such as MAMU-DQA1, which are increased in macrophage activation (101), and the intermediate filament protein gene vimentin (VIM), which is increased in macrophages differentiation and activation (102). None of these genes had appreciable expression in the microglial clusters.
Figure 5 Morphine-related DEGs expression in the different clusters. Bubble plot of selected DEGs between the saline and morphine cells within the clusters. The color indicates the least-square mean expression (log2), and the size of the circle indicates the percentage of cells in which gene expression was detected.
Osteopontin expression is increased in microglia, particularly in white matter microglia, by morphine in SIV-infected monkeys
In examining the DEGs, we were particularly intrigued by the increase in SPP1 (also known as osteopontin), a pleiotropic molecule produced by, and with numerous effects on, macrophage and microglia, as well as with links to neurodegenerative disorders and many other diseases (98, 103–107). Its expression was significantly increased in all microglia clusters, but not in the macrophage cluster (Table 3).
SPP1/osteopontin is known to be expressed by neurons, oligodendrocytes, and microglia, as well as CNS-infiltrating immune cells, with different patterns in different brain regions, and can change in reaction to stimuli (108–113). To address regional and cellular expression, we performed single-nucleus RNA sequencing (snRNA-seq) on two brain regions – frontal lobe white matter and the caudate nucleus. Again, four animals from each group were examined. An average of ~5,600 nuclei were examined from each animal from the white matter and ~5,900 nuclei from the caudate. Clustering and marker gene analysis of the white matter nuclei revealed 1 cluster of oligodendrocyte precursor cells (OPC), 3 clusters of oligodendrocytes, 1 cluster of mixed astrocytes and neurons, and 1 cluster of microglia (containing 4,761 nuclei, 10% of all nuclei) (Figure 6A). In the caudate, 1 cluster of OPC, 3 clusters of oligodendrocytes, 1 cluster of astrocytes, 3 clusters of neurons, and 1 cluster of microglia (containing 3,806 nuclei, 8.1% of all nuclei) were identified (Figure 6B). In both, cluster 4 contained microglia (P2RY12 and CSF1R positive).
Figure 6 SnRNA-seq of brain regions. (A) White Matter and (B) Caudate. Left in each are UMAP plots, right are violin plots of cellular identifying markers. OPC express carbonic anhydrase 10 (CA10), oligodendrocytes myelin oligodendrocyte glycoprotein (MOG), astrocytes aldehyde dehydrogenase 1 family member A1 (ALDH1A1), neurons RNA binding Fox-1 homolog 3 (RBFOX3, also known as neuronal nuclei, NeuN), and microglia purinergic receptor P2Y12 (P2RY12) and colony-stimulating factor-1 receptor (CSF1R). Values on the Y-axis are the least-square mean expression (log2).
The white matter and caudate microglia were then extracted and assessed by graph-based clustering (Figure 7A). Four clusters were found, there was a striking difference in the proportions of the microglia clusters between the white matter and caudate, whereas there were minimal differences in response to morphine treatment (Figure 7B). In each cluster, SPP1 expression was higher in the white matter microglia than the caudate microglia, and morphine only significantly upregulated SPP1 expression in clusters 1 (which comprised the majority of white matter microglial cells) and 3 in white matter microglia (Table 4; Figure 7C).
Figure 7 Microglia in White Matter and Caudate. (A) UMAP with graph-based clustering of microglia identified from the white matter and caudate microglia clusters. (B) Pie charts with areas proportional to the percentage of cells in each treatment condition from the different graph-based clusters of microglia in the white matter and caudate. (C) Violin plots of SPP1/osteopontin expression in microglia in the different clusters, separated by location and treatment. Values on the Y-axis are the least-square mean expression (log2).
Discussion
Substance use disorders, and more specifically abuse of opioids, are important comorbidities among PWH. The current study investigated interactive effects of opioids in the brains of SIV-infected rhesus macaques with suppressed viremia. The primary objective was to examine factors that can contribute to neuroHIV, including the increased viral reservoir found in morphine-treated animals. While HIV encephalitis dominated neuropathogenesis in the pre-treatment era, in the current era of efficacious treatment neuropathological studies have not uncovered the substrate of central nervous system (CNS) dysfunction (114, 115). While ample studies have focused on the characterization of latent virus reservoirs in the peripheral circulating CD4+ T cells and lymphoid tissues, there is a paucity of information on the mechanisms involved in the persistence of HIV reservoirs in the myeloid cells of the CNS.
To assess the impact of infection and cART in the context of morphine dependence, we evaluated the expression of biomarkers as well as lymphocyte subsets before infection, during the acute phase of infection at peak viremia (2 weeks p.i.), post-peak viremia (5 weeks p.i.) as well as at two additional time points post-virus suppression following cART (22 and 35 weeks p.i., which correspond to 17- and 30-weeks post-cART initiation, respectively). Before infection, morphine exposure itself led to an anti-inflammatory phenotype. For example, the expression of adiponectin, a hormone known to attenuate the inflammatory response, was significantly upregulated in the presence of morphine. Further, expression of several other inflammatory mediators and markers such as IL-16, MMP-2, TNFR2, and VCAM-1 were all downregulated in the plasma of uninfected morphine-dependent macaques relative to uninfected saline-treated animals.
Striking results were found during the acute infection phase, which is crucial in determining disease course in untreated infection. The increase in CD8+ T cells and NK cells, characteristic of response to infection in the saline group, was largely attenuated in the morphine group. Expression of CRP, a liver-produced marker of inflammation, increased dramatically in saline-treated SIV-infected animals, while its levels remained unchanged over baseline in the presence of morphine. Similarly, ferritin, an iron-sequestering protein that is elevated during inflammation, increased significantly in the saline versus morphine-dependent macaques. IL-1RA, IL-16, IL-18, MMP2, TNFR2, and VCAM-1 were all significantly lower at this acute time point in the morphine group compared with the saline group. These findings are in keeping with published reports implicating attenuation of the inflammatory immune response by morphine in the periphery of SIV-infected macaques (41). Along these lines, heroin abusers have been shown to exhibit reduced immune responses and increased incidence of infections (116, 117). Furthermore, chronic administration of relatively low doses of morphine in a cohort of rhesus monkeys over a 2-year period was shown to affect immunocompetence, suppressing NK-cell activity and decreasing the percentage of CD4+ circulating lymphocytes (118, 119). In addition a study on morphine-treated SIV-infected monkeys found that they also showed depressed weight gain as well as were more vulnerable to opportunistic infections with quicker progression to simian AIDS (120). Studies in a variety of systems confirm the immune suppressive properties of opioids (121).
Most of the effects of morphine found in the uninfected animals persisted even following cART in the infected, morphine-dependent animals compared with the infected saline group. It is likely that a muted inflammatory response during the initial infection period could likely contribute to increased seeding of virus reservoirs with decreased clearance of cells resistant to the cytotoxic effects of the virus, such as macrophages and microglia. We did not see any change in antiretroviral concentrations in either brain, or biofluids in the context of morphine, likely indicating that drug pharmacokinetics played a negligible role in the maintenance of CNS virus reservoirs in our model system.
Brain-resident myeloid cells – macrophages and microglia – are the primary HIV/SIV infected cells in the brain and comprise the brain reservoir during treatment (48, 49). These cells play multifaceted roles and can control CNS damage owing to their ability to elicit key immune defenses in the CNS and their inherent ability to release proinflammatory cytokines and chemokines in response to cellular insults (122, 123). In the CNS, the deleterious effects of infection include the generation of proinflammatory mediators within these infected cells, which then in turn cause bystander dysfunction and toxicity in uninfected neurons and glia. Drugs of abuse, including opioids such as morphine, also have distinct effects on brain microglia and macrophages.
To assess the effects of morphine on microglia and macrophages, we performed scRNA-seq on brain cells from the two groups of monkeys. The proportion of different microglia subtypes, identified by clustering and proportion of macrophages, did not differ between the groups. Both the microglia and the macrophages showed changes in gene expression between the saline-administered and morphine-dependent animals. While some differences were present between the microglia subtypes, the upregulated and downregulated genes showed commonalities. In the morphine-treated group, expression of neurodegeneration-related genes was increased in all microglia subtypes and macrophages. Enrichment of genes related to the DAM phenotype was present in one of the microglia clusters due to morphine. However, no specific DAM cluster of microglia was identified in the morphine group. While the DAM phenotype is prominent in mouse models of neurodegeneration, its translation to humans is inconsistent (71), and thus calls into question whether a specific DAM microglia cluster could be found in NHP. Pathway analysis was notable for increases in the phagosome formation pathway in microglia as well as the neuroinflammation signaling pathway in microglia and macrophages, as well as other signaling pathways in the different cellular clusters.
The DEGs in context of morphine were different in the macrophage cluster than those in the microglia clusters. In microglia clusters, morphine increased the expression of SPP1/osteopontin, which has been linked to HIV-associated neurocognitive disorders in several studies. While our microglia preparations for the scRNA-seq were derived from the whole brain, to obtain an independent assessment of morphine-related increase in SPP1 transcripts in microglia, we performed snRNA-seq on nuclei isolated from two brain regions, frontal lobe white matter and the caudate. In both regions we could identify a cluster of cells that represented microglia.
Somewhat surprisingly there was an anatomic location difference in SPP1 expression, both for the value for SPP1 expression in the saline group as well as the magnitude and direction of change in the morphine group. Much higher levels of SPP1 were found in white matter microglia than in caudate microglia. In the morphine group the SPP1 expression level was significantly increased in the white matter. Changes in the proportions of different clusters of microglia found in the snRNA-seq analysis, were also found between white matter and caudate, contributing to the differences in SPP1 expression.
Interestingly a study on human brains found higher SPP1 expression in white matter versus grey matter microglia, and the difference increased with aging (124). This study, as well as another in mice (125) also found white matter microglia had higher levels of type II major histocompatibility genes. Another study in mice found distinct clusters of microglia only within the white matter, denoted as white matter-associated microglia (WAMs), and such clusters were enriched SPP1 expression (126). This study found that WAMs increased with aging and expressed a subset of the DAM gene signature identified in Alzheimer’s disease models. In addition to WAMs, other differences reported between white and grey matter microglia include a higher density of microglia in the white matter and increased markers of antigen presentation and phagocytosis (127). Numerous abnormalities in white matter have been characterized in PWH and some are found to be associated with neurocognitive disorders (128, 129). Interestingly, opioids have also been shown to induce white matter damage (130).
The role of SPP1/osteopontin in HIV infection has been assessed in several studies. We first reported increased osteopontin mRNA expression in the brains of SIV-infected rhesus monkeys with encephalitis (131), and that osteopontin protein was increased in the plasma and CSF of animals with SIV encephalitis (132). We subsequently found increased osteopontin levels in the CSF of PWH and in the plasma of PWH with dementia (133). Others have also found correlations with osteopontin. In PWH, regardless of cART treatment, osteopontin was found to be elevated in the plasma (134). In the basal ganglia of SIV-infected monkeys, osteopontin mRNA was highly expressed compared to that in the uninfected animals or those with viral suppression by cART (135). In another study in PWH, osteopontin was found to be increased in the brain and CSF, and importantly, shown to enhance HIV replication in primary human macrophages (136, 137). Further, osteopontin can affect monocyte entry into, and egress from, the brain. Osteopontin is a chemotactic agent for macrophages and increases their survival (138–140). In addition, we found that osteopontin prevents monocyte recirculation after traversing a barrier, and reduces monocyte apoptosis, thus retaining monocytes that enter the brain (132). SPP1 is frequently found in gene lists of activated or disease-linked microglia subsets identified by scRNA-seq in mice and humans (141). Interestingly in mice, a microglia subset expressing SPP1/osteopontin has been found to execute proinflammatory responses (142).
Osteopontin actions are also regulated through several mechanisms. Alternate translational start sites lead to an intracellular and secreted forms (143). The intracellular form predominates in myeloid cells, and has distinct properties compared to the secreted form, which is abundantly produced in T cells (144). Secreted osteopontin can be cleaved by thrombin, leading to enhanced chemotactic activity on monocytes (145). Furthermore, genetic polymorphisms exist in the osteopontin gene, and alternate splicing can produce different isoforms of osteopontin (144). Post-translational phosphorylation also leads to differential activities on macrophage migration and macrophage activation (146). The role of these different forms of osteopontin, and their effect on brain microglia, macrophages, HIV, and the effects of drugs of abuse such as opioids remains a significant gap in our knowledge and awaits future study.
Using scRNA-seq of microglia from SIV-infected animals, we recently found that in untreated animals, microglia exhibited a proinflammatory state, and that viral suppression with cART restored the microglia to a homeostatic-like state (63). Here we find that in the setting of viral suppression with cART, morphine had numerous other effects on the microglia as well as brain macrophages. Strikingly, the 4 largest microglia clusters and the macrophage cluster all showed enriched expression of genes related to neurodegeneration. These genes were curated as a core set of genes exhibiting increased expression in CNS myeloid cells in multiple models of neurodegenerative disease. Interestingly, it was suggested that in neurodegeneration this increase in expression changes how these cells interact with their environment, since the majority of the genes in this classification were annotated with a Gene Ontology (GO) term associated with either the plasma membrane or extracellular space (71). While we did not find separate clusters of CNS myeloid cells in the morphine-treated group samples, all clusters of these cells were enriched for increased expression of neurodegeneration-related genes.
We note that this enrichment of genes related to a neurodegeneration phenotype is in keeping with the analysis identifying enrichment of genes related to a neuroinflammatory pathway. Many studies in human disease and animal models have found that activation of microglia leads to inflammatory molecules in the brain, and such neuroinflammation leads to damage and neurodegeneration (147). While our study focused on microglia, certainly other brain cells, such as astrocytes, also can contribute to the neuroinflammatory environment and result in neurodegenerative processes (148). This can be a two-way (or larger) interaction, as damaged neurons can also signal glia to activate and expressed inflammatory molecules. The ultimate reason for the differential effects of morphine may be intrinsic to microglia or may result from such interactions. While microglia and macrophages have been reported to express the mu opioid receptor (54, 149), the key receptor for effects of morphine, we only found rare microglia or macrophages with expression of the mu or other opioid receptors in our scRNA-seq and scRNA-seq data, although they could be detected in neurons in the caudate, as expected. Another scRNA-seq data set also showed only relatively rare expression of opioid receptors in microglia (87). This is likely due to sensitivity issues in single cell transcriptomic experiments, and the resolution of the target(s) for opioids in the differences found in microglia and brain macrophages will require additional research.
The enrichment for the DAM phenotype in one of the clusters of microglia and macrophages reveals a similar response to damage found in a variety of experimental animal models of CNS disease, and, to an extent, in human disease. In addition, we found that administration of morphine increased the expression levels of genes regulating the neuroinflammation signaling pathway in macrophages. Other myeloid signaling pathways, such as the integrin and triggering receptor expressed on myeloid cells 1 (TREM1) pathway, were also increased in different myeloid clusters and the macrophage cluster. These, along with an increase in the phagosome formation pathway in 3 clusters, connote responsiveness of the microglia and macrophages to morphine within the brains of SIV-infected NHP. While there is evidence of role of TREM1 in resistance to HIV-mediated cytopathogenesis in macrophages (150, 151), this effect, however, has not been examined thus far in the context of opioids.
Very few studies using scRNA-seq for assessing the opioid effects have been performed. In one study, cells from the nucleus accumbens (NAc) were studied following acute morphine treatment in mice. Four hours after morphine administration, a cluster of dopamine receptor 1a expressing neurons, astrocytes, and multiple clusters of oligodendrocytes, including OPC, were found to respond transcriptionally to the drug, with the most robust response found in the oligodendrocytes (152). While we did not examine the NAc, we note that we did not see similar responses in the oligodendrocytes from the two regions we examined in the brain, including the caudate nucleus, also part of the striatum. Additionally, our model was of chronic, not acute, morphine treatment. Another study examined peripheral blood mononuclear cells from opioid-dependent and control subjects (153). Here, in circulating naïve monocytes, suppression of antiviral genes was found, and the response to stimulation in monocytes and other blood cell types was suppressed. In vitro morphine treatment of blood immune cells yielded similar results. We indeed found evidence for such morphine-linked immunosuppression throughout the longitudinal course of infection in the morphine-treated group.
ScRNA-seq CNS studies in the context of HIV infection have also been limited, in part, understandably, to the inaccessibility of the brain during life. However the CSF has been studied and interesting results have been found. The first study examined the CSF cells from PWH virally suppressed on cART compared to uninfected controls and identified a small population of myeloid cells in the PWH subjects with increased expression of genes found in microglia and specifically in neurodegeneration-associated microglia clusters (154). This group then expanded these studies, and again identifying this rare cluster of microglia-like cells in the CSF, as well as finding CD4+ T cells expressing HIV in the CSF from these virally-suppressed subjects (155). Parallel examination of the CSF cells and the brain in SIV-infected rhesus macaques could help clarify the role of this unique subset of myeloid cells in the CSF in HIV/SIV infection, as well as any effect of opioids.
Limitations to our study include the lack of analysis of uninfected animals, in both control and morphine groups, and a relatively small sample size. Only male animals were used in this study, making it difficult to make any inferences on sex-specific differences. Furthermore, in our study, the length of morphine-dependence and duration of viral infection (which was treated early to suppress viremia) was much shorter than that observed in PWH. Numerous other molecules which could potentially affect the HIV reservoirs were differentially regulated by morphine in the brain macrophages and microglia and were not investigated further other than through GSEA and pathway analysis. We could not discern the effects of different forms of osteopontin analysis on the brain reservoir in this model. Other than brain macrophages and microglia, the actions of morphine on other cell types in the brain were not examined.
In summary, morphine, a prototypic opioid, exhibited marked effects on NHP before and after SIV infection. A general immune suppressive environment was found based on immune cell and soluble mediator analysis in the blood. Following sustained viral suppression, we analyzed brain macrophages and microglia, targets for SIV and HIV and mediators of both protective and destructive actions as well as the CNS reservoir for these viruses, by scRNA-seq and snRNA-seq. These findings revealed gene expression changes that could indicate transition to a neurodegenerative phenotype due to morphine administration, and changes that likely contribute to the basis of our prior findings of increased viral reservoirs in the brains of cART-suppressed, morphine-dependent animals. In summary, the overall immune suppressive environment, altered gene expression by CNS myeloid cells, and the actions of one of the morphine-induced genes, SPP1/osteopontin, yield novel insights into the interactions of the dual-linked health crises of HIV and opioids.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, GSE209606. https://www.ncbi.nlm.nih.gov/, GSE195574.
Ethics statement
The animal study was reviewed and approved by the UNMC Institutional Animal Care and Use Committee under protocol 16-073-FC.
Author contributions
Conceptualization: HF, SNB, SB. Investigation: HF, BM, SL, KE, PP, SC, AA, GK, SD. Data – curation: HF, BM, SD. Data – formal analysis: HF, MN, BD, CF. Project administration: HF, BM, BL, JE, CF. Funding acquisition: HF, SNB, SB, CF. Writing – original draft: HF. Writing – reviewing, editing and revision: HF, MN, BM, ML, KE, PP, SC, AA, GK, JE, CG, SD, CF, SNB, SB. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grant #s R01DA043164 (to HF, SB, and SNB) and U01DA053624 (to HF and SB) from the National Institute of Drug Abuse, P30MH062261 (to HF and SB) from the National Institute of Mental Health, and R01AI124965 (to CF) from the National Institute of Allergy and Infectious Diseases, all from the National Institutes of Health. The UNMC DNA Sequencing Core receives support from NIH grants P20GM103427 and P30CA036727. The UNMC Flow Cytometry Research Facility receives state funds from the Nebraska Research Initiative and the NIH grant P30CA036727. The UNMC Bioinformatics and Systems Biology Core is partly supported by NIH awards P20GM103427, P30CA036727, P30MH062261, and U54GM115458. Major instrumentation has been provided by the UNMC Office of the Vice Chancellor for Research, The University of Nebraska Foundation, the Nebraska Banker’s Fund, and by the NIH Shared Instrument Program.
Acknowledgments
We thank Kabita Pandey, Lindsey A. Knight, Alexandra Sheldon, Chase Ochs and Natasha Ferguson for experimental assistance. This paper is dedicated to the memory of Brenda Morsey.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1012884/full#supplementary-material
References
1. Buchacz K, Armon C, Palella FJ Jr., Novak RM, Fuhrer J, Tedaldi E, et al. The HIV outpatient study-25 years of HIV patient care and epidemiologic research. Open Forum Infect Dis (2020) 7(5):ofaa123. doi: 10.1093/ofid/ofaa123
2. Michaels SH, Clark R, Kissinger P. Declining morbidity and mortality among patients with advanced human immunodeficiency virus infection. N Engl J Med (1998) 339(6):405–6. doi: 10.1056/NEJM199808063390612
3. Moore RD, Chaisson RE. Natural history of HIV infection in the era of combination antiretroviral therapy. AIDS. (1999) 13(14):1933–42. doi: 10.1097/00002030-199910010-00017
4. Palella FJ Jr., Baker RK, Moorman AC, Chmiel JS, Wood KC, Brooks JT, et al. Mortality in the highly active antiretroviral therapy era: changing causes of death and disease in the HIV outpatient study. J Acquir Immune Defic Syndr (2006) 43(1):27–34. doi: 10.1097/01.qai.0000233310.90484.16
5. Janssen RS, Nwanyanwu OC, Selik RM, Stehr-Green JK. Epidemiology of human immunodeficiency virus encephalopathy in the united states. Neurology. (1992) 42(8):1472–6. doi: 10.1212/WNL.42.8.1472
6. McArthur JC, Hoover DR, Bacellar H, Miller EN, Cohen BA, Becker JT, et al. Dementia in AIDS patients: incidence and risk factors. Multicenter AIDS Cohort Study. Neurology. (1993) 43(11):2245–52. doi: 10.1212/WNL.43.11.2245
7. Heaton RK, Clifford DB, Franklin DR Jr., Woods SP, Ake C, Vaida F, et al. HIV-Associated neurocognitive disorders persist in the era of potent antiretroviral therapy: CHARTER study. Neurology. (2010) 75(23):2087–96. doi: 10.1212/WNL.0b013e318200d727
8. Wang Y, Liu M, Lu Q, Farrell M, Lappin JM, Shi J, et al. Global prevalence and burden of HIV-associated neurocognitive disorder: A meta-analysis. Neurology. (2020) 95(19):e2610–e21. doi: 10.1212/WNL.0000000000010752
9. Chakrabarti L, Hurtrel M, Maire MA, Vazeux R, Dormont D, Montagnier L, et al. Early viral replication in the brain of SIV-infected rhesus monkeys. Am J Pathol (1991) 139(6):1273–80.
10. Marcondes MC, Burudi EM, Huitron-Resendiz S, Sanchez-Alavez M, Watry D, Zandonatti M, et al. Highly activated CD8(+) T cells in the brain correlate with early central nervous system dysfunction in simian immunodeficiency virus infection. J Immunol (2001) 167(9):5429–38. doi: 10.4049/jimmunol.167.9.5429
11. Williams KC, Corey S, Westmoreland SV, Pauley D, Knight H, deBakker C, et al. Perivascular macrophages are the primary cell type productively infected by simian immunodeficiency virus in the brains of macaques: implications for the neuropathogenesis of AIDS. J Exp Med (2001) 193(8):905–15. doi: 10.1084/jem.193.8.905
12. Lackner AA, Vogel P, Ramos RA, Kluge JD, Marthas M. Early events in tissues during infection with pathogenic (SIVmac239) and nonpathogenic (SIVmac1A11) molecular clones of simian immunodeficiency virus. Am J Pathol (1994) 145(2):428–39.
13. Davis LE, Hjelle BL, Miller VE, Palmer DL, Llewellyn AL, Merlin TL, et al. Early viral brain invasion in iatrogenic human immunodeficiency virus infection. Neurology. (1992) 42(9):1736–9. doi: 10.1212/WNL.42.9.1736
14. Cochrane CR, Angelovich TA, Byrnes SJ, Waring E, Guanizo AC, Trollope GS, et al. Intact HIV proviruses persist in the brain despite viral suppression with ART. Ann Neurol (2022) 92(4):532–44. doi: 10.1002/ana.26456
15. Avalos CR, Abreu CM, Queen SE, Li M, Price S, Shirk EN, et al. Brain macrophages in simian immunodeficiency virus-infected, antiretroviral-suppressed macaques: a functional latent reservoir. MBio (2017) 8(4):e01186–17. doi: 10.1128/mBio.01186-17
16. Letendre S, Marquie-Beck J, Capparelli E, Best B, Clifford D, Collier AC, et al. Validation of the CNS penetration-effectiveness rank for quantifying antiretroviral penetration into the central nervous system. Arch Neurol (2008) 65(1):65–70. doi: 10.1001/archneurol.2007.31
17. CDC. Estimated HIV incidence and prevalence in the united states, 2014–2018. In: HIV Surveillance supplemental report 2020. USA: Center for Disease Control (2020). Available at: https://www.cdc.gov/hiv/pdf/library/reports/surveillance/cdc-hiv-surveillance-supplemental-report-vol-25-1.pdf.
18. Murphy A, Barbaro J, Martinez-Aguado P, Chilunda V, Jaureguiberry-Bravo M, Berman JW. The effects of opioids on HIV neuropathogenesis. Front Immunol (2019) 10:2445. doi: 10.3389/fimmu.2019.02445
19. Hodder SL, Feinberg J, Strathdee SA, Shoptaw S, Altice FL, Ortenzio L, et al. The opioid crisis and HIV in the USA: deadly synergies. Lancet. (2021) 397(10279):1139–50. doi: 10.1016/S0140-6736(21)00391-3
20. CDC. HIV Surveillance report. USA: Center for Disease Control (2022). Available at: https://www.cdc.gov/hiv/library/reports/hiv-surveillance.html.
21. CDC. HIV And substance use. USA: Center for Disease Control (2021). Available at: https://www.cdc.gov/hiv/basics/hiv-transmission/substance-use.html.
22. Fitting S, McRae M, Hauser KF. Opioid and neuroHIV comorbidity - current and future perspectives. J Neuroimmune Pharmacol (2020) 15(4):584–627. doi: 10.1007/s11481-020-09941-8
23. Cunningham CO. Opioids and HIV infection: From pain management to addiction treatment. Top Antivir Med (2018) 25(4):143–6.
24. Nath A, Anderson C, Jones M, Maragos W, Booze R, Mactutus C, et al. Neurotoxicity and dysfunction of dopaminergic systems associated with AIDS dementia. J Psychopharmacol (Oxford England). (2000) 14(3):222–7. doi: 10.1177/026988110001400305
25. Nath A, Hauser KF, Wojna V, Booze RM, Maragos W, Prendergast M, et al. Molecular basis for interactions of HIV and drugs of abuse. J Acquir Immune Defic Syndr (2002) 31 Suppl 2:S62–9. doi: 10.1097/00126334-200210012-00006
26. Hauser KF, El-Hage N, Buch S, Berger JR, Tyor WR, Nath A, et al. Molecular targets of opiate drug abuse in neuroAIDS. Neurotox Res (2005) 8(1-2):63–80. doi: 10.1007/BF03033820
27. Burdo TH, Katner SN, Taffe MA, Fox HS. Neuroimmunity, drugs of abuse, and neuroAIDS. J Neuroimmune Pharmacol (2006) 1(1):41–9. doi: 10.1007/s11481-005-9001-3
28. Chilunda V, Calderon TM, Martinez-Aguado P, Berman JW. The impact of substance abuse on HIV-mediated neuropathogenesis in the current ART era. Brain Res (2019) 1724:146426. doi: 10.1016/j.brainres.2019.146426
29. Bell JE, Brettle RP, Chiswick A, Simmonds P. HIV Encephalitis, proviral load and dementia in drug users and homosexuals with AIDS. effect of neocortical involvement. Brain. (1998) 121(Pt 11):2043–52.
30. Smith DB, Simmonds P, Bell JE. Brain viral burden, neuroinflammation and neurodegeneration in HAART-treated HIV positive injecting drug users. J Neurovirol. (2014) 20(1):28–38. doi: 10.1007/s13365-013-0225-3
31. Buch S, Periyasamy P, Thangaraj A, Sil S, Chivero ET, Tripathi A. Opioid-mediated HIV-1 immunopathogenesis. J Neuroimmune Pharmacol (2020) 15(4):628–42. doi: 10.1007/s11481-020-09960-5
32. Burudi EM, Fox HS. Simian immunodeficiency virus model of HIV-induced central nervous system dysfunction. Adv Virus Res (2001) 56:435–68. doi: 10.1016/S0065-3527(01)56035-2
33. Rahman MA, Robert-Guroff M. Accelerating HIV vaccine development using non-human primate models. Expert Rev Vaccines (2019) 18(1):61–73. doi: 10.1080/14760584.2019.1557521
34. Martins MA, Watkins DI. What is the predictive value of animal models for vaccine efficacy in humans? rigorous simian immunodeficiency virus vaccine trials can be instructive. Cold Spring Harb Perspect Biol (2018) 10(4):a029504. doi: 10.1101/cshperspect.a029504
35. Mallard J, Williams KC. Animal models of HIV-associated disease of the central nervous system. Handb Clin Neurol (2018) 152:41–53. doi: 10.1016/B978-0-444-63849-6.00004-9
36. Goulder P, Deeks SG. HIV Control: Is getting there the same as staying there? PloS Pathog (2018) 14(11):e1007222. doi: 10.1371/journal.ppat.1007222
37. Estes JD, LeGrand R, Petrovas C. Visualizing the immune system: Providing key insights into HIV/SIV infections. Front Immunol (2018) 9:423. doi: 10.3389/fimmu.2018.00423
38. Bokhari SM, Hegde R, Callen S, Yao H, Adany I, Li Q, et al. Morphine potentiates neuropathogenesis of SIV infection in rhesus macaques. J Neuroimmune Pharmacol (2011) 6(4):626–39. doi: 10.1007/s11481-011-9272-9
39. Kumar R, Orsoni S, Norman L, Verma AS, Tirado G, Giavedoni LD, et al. Chronic morphine exposure causes pronounced virus replication in cerebral compartment and accelerated onset of AIDS in SIV/SHIV-infected Indian rhesus macaques. Virology. (2006) 354(1):192–206. doi: 10.1016/j.virol.2006.06.020
40. Marcario JK, Pendyala G, Riazi M, Fleming K, Marquis J, Callen S, et al. Effects of morphine on behavioral task performance in SIV-infected rhesus macaques. J Neuroimmune Pharmacol (2016) 11(2):348–57. doi: 10.1007/s11481-016-9667-8
41. Marcario JK, Riazi M, Adany I, Kenjale H, Fleming K, Marquis J, et al. Effect of morphine on the neuropathogenesis of SIVmac infection in Indian rhesus macaques. J Neuroimmune Pharmacol (2008) 3(1):12–25. doi: 10.1007/s11481-007-9085-z
42. Marcondes MC, Flynn C, Watry DD, Zandonatti M, Fox HS. Methamphetamine increases brain viral load and activates natural killer cells in simian immunodeficiency virus-infected monkeys. Am J Pathol (2010) 177(1):355–61. doi: 10.2353/ajpath.2010.090953
43. Donahoe RM. Multiple ways that drug abuse might influence AIDS progression: clues from a monkey model. J Neuroimmunol. (2004) 147(1-2):28–32. doi: 10.1016/j.jneuroim.2003.10.011
44. Kumar R, Torres C, Yamamura Y, Rodriguez I, Martinez M, Staprans S, et al. Modulation by morphine of viral set point in rhesus macaques infected with simian immunodeficiency virus and simian-human immunodeficiency virus. J Virol (2004) 78(20):11425–8. doi: 10.1128/JVI.78.20.11425-11428.2004
45. Cornwell WD, Wagner W, Lewis MG, Fan X, Rappaport J, Rogers TJ. Effect of chronic morphine administration on circulating dendritic cells in SIV-infected rhesus macaques. J Neuroimmunol. (2016) 295-296:30–40. doi: 10.1016/j.jneuroim.2016.04.007
46. Chuang LF, Killam KF Jr., Chuang RY. Increased replication of simian immunodeficiency virus in CEM x174 cells by morphine sulfate. Biochem Biophys Res Commun (1993) 195(3):1165–73. doi: 10.1006/bbrc.1993.2167
47. Joseph SB, Arrildt KT, Sturdevant CB, Swanstrom R. HIV-1 target cells in the CNS. J Neurovirol. (2015) 21(3):276–89. doi: 10.1007/s13365-014-0287-x
48. Wallet C, De Rovere M, Van Assche J, Daouad F, De Wit S, Gautier V, et al. Microglial cells: The main HIV-1 reservoir in the brain. Front Cell Infect Microbiol (2019) 9:362. doi: 10.3389/fcimb.2019.00362
49. Abreu C, Shirk EN, Queen SE, Beck SE, Mangus LM, Pate KAM, et al. Brain macrophages harbor latent, infectious simian immunodeficiency virus. AIDS. (2019) 33 Suppl 2:S181–S8. doi: 10.1097/QAD.0000000000002269
50. Ko A, Kang G, Hattler JB, Galadima HI, Zhang J, Li Q, et al. Macrophages but not astrocytes harbor HIV DNA in the brains of HIV-1-Infected aviremic individuals on suppressive antiretroviral therapy. J Neuroimmune Pharmacol (2019) 14(1):110–9. doi: 10.1007/s11481-018-9809-2
51. Mohammadzadeh N, Roda W, Branton WG, Clain J, Rabezanahary H, Zghidi-Abouzid O, et al. Lentiviral infections persist in brain despite effective antiretroviral therapy and neuroimmune activation. mBio. (2021) 12(6):e0278421. doi: 10.1128/mBio.02784-21
52. Borrajo A, Svicher V, Salpini R, Pellegrino M, Aquaro S. Crucial role of central nervous system as a viral anatomical compartment for HIV-1 infection. Microorganisms. (2021) 9(12):2537. doi: 10.3390/microorganisms9122537
53. Hu S, Chao CC, Hegg CC, Thayer S, Peterson PK. Morphine inhibits human microglial cell production of, and migration towards, RANTES. J Psychopharmacol (Oxford England). (2000) 14(3):238–43. doi: 10.1177/026988110001400307
54. Wen S, Jiang Y, Liang S, Cheng Z, Zhu X, Guo Q. Opioids regulate the immune system: Focusing on macrophages and their organelles. Front Pharmacol (2021) 12:814241. doi: 10.3389/fphar.2021.814241
55. Yang Y, Sun Y, Hu R, Yan J, Wang Z, Li W, et al. Morphine promotes microglial activation by upregulating the EGFR/ERK signaling pathway. PloS One (2021) 16(9):e0256870. doi: 10.1371/journal.pone.0256870
56. Mali AS, Novotny J. Opioid receptor activation suppresses the neuroinflammatory response by promoting microglial M2 polarization. Mol Cell Neurosci (2022) 121:103744. doi: 10.1016/j.mcn.2022.103744
57. Acharya A, Olwenyi OA, Thurman M, Pandey K, Morsey BM, Lamberty B, et al. Chronic morphine administration differentially modulates viral reservoirs in SIVmac251 infected rhesus macaque model. J Virol (2020) 95(5):e01657–20. doi: 10.1101/2020.09.09.290700
58. Morsey B, Niu M, Dyavar SR, Fletcher CV, Lamberty BG, Emanuel K, et al. Cryopreservation of microglia enables single-cell RNA sequencing with minimal effects on disease-related gene expression patterns. iScience. (2021) 24(4):102357. doi: 10.1016/j.isci.2021.102357
59. Emanuel KM, Runner K, Brodnik ZD, Morsey BM, Lamberty BG, Johnson HS, et al. Deprenyl reduces inflammation during acute SIV infection. iScience. (2022) 25(5):104207. doi: 10.1016/j.isci.2022.104207
60. Fletcher CV, Staskus K, Wietgrefe SW, Rothenberger M, Reilly C, Chipman JG, et al. Persistent HIV-1 replication is associated with lower antiretroviral drug concentrations in lymphatic tissues. Proc Natl Acad Sci U S A. (2014) 111(6):2307–12. doi: 10.1073/pnas.1318249111
61. Fletcher CV, Podany AT, Thorkelson A, Winchester LC, Mykris T, Anderson J, et al. The lymphoid tissue pharmacokinetics of tenofovir disoproxil fumarate and tenofovir alafenamide in HIV-infected persons. Clin Pharmacol Ther (2020) 108(5):971–5. doi: 10.1002/cpt.1883
62. Fletcher CV, Kroon E, Schacker T, Pinyakorn S, Chomont N, Chottanapund S, et al. Persistent HIV transcription and variable antiretroviral drug penetration in lymph nodes during plasma viral suppression. AIDS. (2022) 36(7):985–90. doi: 10.1097/QAD.0000000000003201
63. Trease AJ, Niu M, Morsey B, Guda C, Byrareddy SN, Buch S, et al. Antiretroviral therapy restores the homeostatic state of microglia in SIV-infected rhesus macaques. J Leukoc Biol (2022) 112(5):969–81. doi: 10.1002/JLB.3HI0422-635R
64. Niu M, Morsey B, Lamberty BG, Emanuel K, Yu F, Leon-Rivera R, et al. Methamphetamine increases the proportion of SIV-infected Microglia/Macrophages, alters metabolic pathways, and elevates cell death pathways: A single-cell analysis. Viruses (2020) 12(11):1297. doi: 10.3390/v12111297
65. McGinnis CS, Murrow LM, Gartner ZJ. DoubletFinder: Doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst (2019) 8(4):329–37.e4. doi: 10.1016/j.cels.2019.03.003
66. Hafemeister C, Satija R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biol (2019) 20(1):296. doi: 10.1186/s13059-019-1874-1
67. Korsunsky I, Millard N, Fan J, Slowikowski K, Zhang F, Wei K, et al. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods (2019) 16(12):1289–96. doi: 10.1038/s41592-019-0619-0
68. Becht E, McInnes L, Healy J, Dutertre CA, Kwok IWH, Ng LG, et al. Dimensionality reduction for visualizing single-cell data using UMAP. Nat Biotechnol (2018) 37:38–44. doi: 10.1038/nbt.4314
69. Blondel VD, Guillaume J-L, Lambiotte R, Lefebre E. Fast unfolding of communities in large networks. J Stat Mech (2008) 2008(10):P10008. doi: 10.1088/1742-5468/2008/10/P10008
70. Thrupp N, Sala Frigerio C, Wolfs L, Skene NG, Fattorelli N, Poovathingal S, et al. Single-nucleus RNA-seq is not suitable for detection of microglial activation genes in humans. Cell Rep (2020) 32(13):108189. doi: 10.1016/j.celrep.2020.108189
71. Friedman BA, Srinivasan K, Ayalon G, Meilandt WJ, Lin H, Huntley MA, et al. Diverse brain myeloid expression profiles reveal distinct microglial activation states and aspects of alzheimer's disease not evident in mouse models. Cell Rep (2018) 22(3):832–47. doi: 10.1016/j.celrep.2017.12.066
72. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
73. QIAGEN. Ingenuity pathway analysis . Available at: https://digitalinsights.qiagen.com/IPA.
74. Ouchi N, Walsh K. Adiponectin as an anti-inflammatory factor. Clin Chim Acta (2007) 380(1-2):24–30. doi: 10.1016/j.cca.2007.01.026
75. Cruikshank WW, Kornfeld H, Center DM. Interleukin-16. J Leukoc Biol (2000) 67(6):757–66. doi: 10.1002/jlb.67.6.757
76. Hannocks MJ, Zhang X, Gerwien H, Chashchina A, Burmeister M, Korpos E, et al. The gelatinases, MMP-2 and MMP-9, as fine tuners of neuroinflammatory processes. Matrix Biol (2019) 75-76:102–13. doi: 10.1016/j.matbio.2017.11.007
77. Dibbs Z, Thornby J, White BG, Mann DL. Natural variability of circulating levels of cytokines and cytokine receptors in patients with heart failure: implications for clinical trials. J Am Coll Cardiol (1999) 33(7):1935–42. doi: 10.1016/S0735-1097(99)00130-8
78. Wellicome SM, Kapahi P, Mason JC, Lebranchu Y, Yarwood H, Haskard DO. Detection of a circulating form of vascular cell adhesion molecule-1: raised levels in rheumatoid arthritis and systemic lupus erythematosus. Clin Exp Immunol (1993) 92(3):412–8. doi: 10.1111/j.1365-2249.1993.tb03413.x
79. Luan YY, Yao YM. The clinical significance and potential role of c-reactive protein in chronic inflammatory and neurodegenerative diseases. Front Immunol (2018) 9:1302. doi: 10.3389/fimmu.2018.01302
80. Kell DB, Pretorius E. Serum ferritin is an important inflammatory disease marker, as it is mainly a leakage product from damaged cells. Metallomics. (2014) 6(4):748–73. doi: 10.1039/C3MT00347G
81. Dinarello CA. Overview of the IL-1 family in innate inflammation and acquired immunity. Immunol Rev (2018) 281(1):8–27. doi: 10.1111/imr.12621
82. Ricklin D, Reis ES, Mastellos DC, Gros P, Lambris JD. Complement component C3 - the "Swiss army knife" of innate immunity and host defense. Immunol Rev (2016) 274(1):33–58. doi: 10.1111/imr.12500
83. Menten P, Wuyts A, Van Damme J. Macrophage inflammatory protein-1. Cytokine Growth Factor Rev (2002) 13(6):455–81. doi: 10.1016/S1359-6101(02)00045-X
84. Schutyser E, Richmond A, Van Damme J. Involvement of CC chemokine ligand 18 (CCL18) in normal and pathological processes. J Leukoc Biol (2005) 78(1):14–26. doi: 10.1189/jlb.1204712
85. Dyavar SR, Kumar S, Gautam N, Podany AT, Winchester LC, Weinhold JA, et al. Intramuscular and subcutaneous administration of antiretroviral drugs, compared with oral, enhances delivery to lymphoid tissues in BALB/c mice. J Antimicrob Chemother (2021) 76(10):2651–8. doi: 10.1093/jac/dkab228
86. Dyavar SR, Gautam N, Podany AT, Winchester LC, Weinhold JA, Mykris TM, et al. Assessing the lymphoid tissue bioavailability of antiretrovirals in human primary lymphoid endothelial cells and in mice. J Antimicrob Chemother (2019) 74(10):2974–8. doi: 10.1093/jac/dkz273
87. Hammond TR, Dufort C, Dissing-Olesen L, Giera S, Young A, Wysoker A, et al. Single-cell RNA sequencing of microglia throughout the mouse lifespan and in the injured brain reveals complex cell-state changes. Immunity. (2019) 50(1):253–71.e6. doi: 10.1016/j.immuni.2018.11.004
88. Keren-Shaul H, Spinrad A, Weiner A, Matcovitch-Natan O, Dvir-Szternfeld R, Ulland TK, et al. A unique microglia type associated with restricting development of alzheimer's disease. Cell. (2017) 169(7):1276–90.e17. doi: 10.1016/j.cell.2017.05.018
89. Sala Frigerio C, Wolfs L, Fattorelli N, Thrupp N, Voytyuk I, Schmidt I, et al. The major risk factors for alzheimer's disease: Age, sex, and genes modulate the microglia response to abeta plaques. Cell Rep (2019) 27(4):1293–306.e6. doi: 10.1016/j.celrep.2019.03.099
90. Butovsky O, Weiner HL. Microglial signatures and their role in health and disease. Nat Rev Neurosci (2018) 19(10):622–35. doi: 10.1038/s41583-018-0057-5
91. Gosselin D, Skola D, Coufal NG, Holtman IR, Schlachetzki JCM, Sajti E, et al. An environment-dependent transcriptional network specifies human microglia identity. Science (2017) 356(6344):eaal3222. doi: 10.1126/science.aal3222
92. Chaum E, Yin J, Yang H, Thomas F, Lang JC. Quantitative AP-1 gene regulation by oxidative stress in the human retinal pigment epithelium. J Cell Biochem (2009) 108(6):1280–91. doi: 10.1002/jcb.22358
93. Gilchrist M, Thorsson V, Li B, Rust AG, Korb M, Roach JC, et al. Systems biology approaches identify ATF3 as a negative regulator of toll-like receptor 4. Nature. (2006) 441(7090):173–8. doi: 10.1038/nature04768
94. Neubrand VE, Forte-Lago I, Caro M, Delgado M. The atypical RhoGTPase RhoE/Rnd3 is a key molecule to acquire a neuroprotective phenotype in microglia. J Neuroinflammation. (2018) 15(1):343. doi: 10.1186/s12974-018-1386-z
95. Arlt A, Schafer H. Role of the immediate early response 3 (IER3) gene in cellular stress response, inflammation and tumorigenesis. Eur J Cell Biol (2011) 90(6-7):545–52. doi: 10.1016/j.ejcb.2010.10.002
96. An Q, Sun C, Li R, Chen S, Gu X, An S, et al. Calcitonin gene-related peptide regulates spinal microglial activation through the histone H3 lysine 27 trimethylation via enhancer of zeste homolog-2 in rats with neuropathic pain. J Neuroinflammation. (2021) 18(1):117. doi: 10.1186/s12974-021-02168-1
97. Sedlyarov V, Eichner R, Girardi E, Essletzbichler P, Goldmann U, Nunes-Hasler P, et al. The bicarbonate transporter SLC4A7 plays a key role in macrophage phagosome acidification. Cell Host Microbe (2018) 23(6):766–74.e5. doi: 10.1016/j.chom.2018.04.013
98. Rosmus DD, Lange C, Ludwig F, Ajami B, Wieghofer P. The role of osteopontin in microglia biology: Current concepts and future perspectives. Biomedicines. (2022) 10(4):840. doi: 10.3390/biomedicines10040840
99. Huysamen C, Brown GD. The fungal pattern recognition receptor, dectin-1, and the associated cluster of c-type lectin-like receptors. FEMS Microbiol Lett (2009) 290(2):121–8. doi: 10.1111/j.1574-6968.2008.01418.x
100. Ennerfelt H, Frost EL, Shapiro DA, Holliday C, Zengeler KE, Voithofer G, et al. SYK coordinates neuroprotective microglial responses in neurodegenerative disease. Cell (2022) 185(22):4135–4152.e22. doi: 10.1016/j.cell.2022.09.030
101. Celada A, Klemsz MJ, Maki RA. Interferon-gamma activates multiple pathways to regulate the expression of the genes for major histocompatibility class II I-a beta, tumor necrosis factor and complement component C3 in mouse macrophages. Eur J Immunol (1989) 19(6):1103–9. doi: 10.1002/eji.1830190621
102. Benes P, Maceckova V, Zdrahal Z, Konecna H, Zahradnickova E, Muzik J, et al. Role of vimentin in regulation of monocyte/macrophage differentiation. Differentiation. (2006) 74(6):265–76. doi: 10.1111/j.1432-0436.2006.00077.x
103. Cappellano G, Vecchio D, Magistrelli L, Clemente N, Raineri D, Barbero Mazzucca C, et al. The yin-yang of osteopontin in nervous system diseases: damage versus repair. Neural Regener Res (2021) 16(6):1131–7. doi: 10.4103/1673-5374.300328
104. Carecchio M, Comi C. The role of osteopontin in neurodegenerative diseases. J Alzheimers Dis (2011) 25(2):179–85. doi: 10.3233/JAD-2011-102151
105. Icer MA, Gezmen-Karadag M. The multiple functions and mechanisms of osteopontin. Clin Biochem (2018) 59:17–24. doi: 10.1016/j.clinbiochem.2018.07.003
106. Rittling SR. Osteopontin in macrophage function. Expert Rev Mol Med (2011) 13:e15. doi: 10.1017/S1462399411001839
107. Yim A, Smith C, Brown AM. Osteopontin/secreted phosphoprotein-1 harnesses glial-, immune-, and neuronal cell ligand-receptor interactions to sense and regulate acute and chronic neuroinflammation. Immunol Rev (2022) 311(1):224–33. doi: 10.1111/imr.13081
108. Kim SY, Choi YS, Choi JS, Cha JH, Kim ON, Lee SB, et al. Osteopontin in kainic acid-induced microglial reactions in the rat brain. Mol Cells (2002) 13(3):429–35.
109. Lee MY, Shin SL, Choi YS, Kim EJ, Cha JH, Chun MH, et al. Transient upregulation of osteopontin mRNA in hippocampus and striatum following global forebrain ischemia in rats. Neurosci Lett (1999) 271(2):81–4. doi: 10.1016/S0304-3940(99)00526-1
110. Shi L, Sun Z, Su W, Xu F, Xie D, Zhang Q, et al. Treg cell-derived osteopontin promotes microglia-mediated white matter repair after ischemic stroke. Immunity. (2021) 54(7):1527–42.e8. doi: 10.1016/j.immuni.2021.04.022
111. Silva K, Hope-Lucas C, White T, Hairston TK, Rameau T, Brown A. Cortical neurons are a prominent source of the proinflammatory cytokine osteopontin in HIV-associated neurocognitive disorders. J Neurovirol. (2015) 21(2):174–85. doi: 10.1007/s13365-015-0317-3
112. Sugiyama Y, Oishi T, Yamashita A, Murata Y, Yamamoto T, Takashima I, et al. Neuronal and microglial localization of secreted phosphoprotein 1 (osteopontin) in intact and damaged motor cortex of macaques. Brain Res (2019) 1714:52–64. doi: 10.1016/j.brainres.2019.02.021
113. Wang X, Louden C, Yue TL, Ellison JA, Barone FC, Solleveld HA, et al. Delayed expression of osteopontin after focal stroke in the rat. J Neurosci (1998) 18(6):2075–83. doi: 10.1523/JNEUROSCI.18-06-02075.1998
114. Gelman BB. Neuropathology of HAND with suppressive antiretroviral therapy: Encephalitis and neurodegeneration reconsidered. Curr HIV/AIDS Rep (2015) 12(2):272–9. doi: 10.1007/s11904-015-0266-8
115. Rosenthal J, Tyor W. Aging, comorbidities, and the importance of finding biomarkers for HIV-associated neurocognitive disorders. J Neurovirol. (2019) 25(5):673–85. doi: 10.1007/s13365-019-00735-0
116. Ocasio FM, Jiang Y, House SD, Chang SL. Chronic morphine accelerates the progression of lipopolysaccharide-induced sepsis to septic shock. J Neuroimmunol. (2004) 149(1-2):90–100. doi: 10.1016/j.jneuroim.2003.12.016
117. Vallejo R, de Leon-Casasola O, Benyamin R. Opioid therapy and immunosuppression: a review. Am J Ther (2004) 11(5):354–65. doi: 10.1097/01.mjt.0000132250.95650.85
118. Carr DJ, France CP. Immune alterations in morphine-treated rhesus monkeys. J Pharmacol Exp Ther (1993) 267(1):9–15. doi: 10.1007/978-1-4615-2980-4_6
119. Carr DJ, France CP. Immune alterations in chronic morphine-treated rhesus monkeys. Adv Exp Med Biol (1993) 335:35–9. doi: 10.1007/978-1-4615-2980-4_6
120. Killam KF, Chuang LF, Chuang RY. Opioid dependency and the progression of simian AIDS: opioid dependency and behavioral observations. Adv Exp Med Biol (1996) 402:43–51. doi: 10.1007/978-1-4613-0407-4_7
121. Eisenstein TK. The role of opioid receptors in immune system function. Front Immunol (2019) 10:2904. doi: 10.3389/fimmu.2019.02904
122. Ransohoff RM, El Khoury J. Microglia in health and disease. Cold Spring Harb Perspect Biol (2015) 8(1):a020560. doi: 10.1101/cshperspect.a020560
123. Prinz M, Jung S, Priller J. Microglia biology: One century of evolving concepts. Cell. (2019) 179(2):292–311. doi: 10.1016/j.cell.2019.08.053
124. Sankowski R, Bottcher C, Masuda T, Geirsdottir L, Sagar, Sindram E, et al. Mapping microglia states in the human brain through the integration of high-dimensional techniques. Nat Neurosci (2019) 22(12):2098–110. doi: 10.1038/s41593-019-0532-y
125. Hagemeyer N, Hanft KM, Akriditou MA, Unger N, Park ES, Stanley ER, et al. Microglia contribute to normal myelinogenesis and to oligodendrocyte progenitor maintenance during adulthood. Acta Neuropathol. (2017) 134(3):441–58. doi: 10.1007/s00401-017-1747-1
126. Safaiyan S, Besson-Girard S, Kaya T, Cantuti-Castelvetri L, Liu L, Ji H, et al. White matter aging drives microglial diversity. Neuron. (2021) 109(7):1100–17.e10. doi: 10.1016/j.neuron.2021.01.027
127. Amor S, McNamara NB, Gerrits E, Marzin MC, Kooistra SM, Miron VE, et al. White matter microglia heterogeneity in the CNS. Acta Neuropathol. (2022) 143(2):125–41. doi: 10.1007/s00401-021-02389-x
128. Stubbe-Drger B, Deppe M, Mohammadi S, Keller SS, Kugel H, Gregor N, et al. Early microstructural white matter changes in patients with HIV: a diffusion tensor imaging study. BMC Neurol (2012) 12:23. doi: 10.1186/1471-2377-12-23
129. Alakkas A, Ellis RJ, Watson CW, Umlauf A, Heaton RK, Letendre S, et al. White matter damage, neuroinflammation, and neuronal integrity in HAND. J Neurovirol. (2019) 25(1):32–41. doi: 10.1007/s13365-018-0682-9
130. Hampton WH, Hanik IM, Olson IR. Substance abuse and white matter: Findings, limitations, and future of diffusion tensor imaging research. Drug Alcohol Depend. (2019) 197:288–98. doi: 10.1016/j.drugalcdep.2019.02.005
131. Roberts ES, Zandonatti MA, Watry DD, Madden LJ, Henriksen SJ, Taffe MA, et al. Induction of pathogenic sets of genes in macrophages and neurons in NeuroAIDS. Am J Pathol (2003) 162(6):2041–57. doi: 10.1016/S0002-9440(10)64336-2
132. Burdo TH, Wood MR, Fox HS. Osteopontin prevents monocyte recirculation and apoptosis. J Leukoc Biol (2007) 81(6):1504–11. doi: 10.1189/jlb.1106711
133. Burdo TH, Ellis RJ, Fox HS. Osteopontin is increased in HIV-associated dementia. J Infect Dis (2008) 198(5):715–22. doi: 10.1086/590504
134. Chagan-Yasutan H, Saitoh H, Ashino Y, Arikawa T, Hirashima M, Li S, et al. Persistent elevation of plasma osteopontin levels in HIV patients despite highly active antiretroviral therapy. Tohoku J Exp Med (2009) 218(4):285–92. doi: 10.1620/tjem.218.285
135. Solis-Leal A, Siddiqui S, Wu F, Mohan M, Hu W, Doyle-Meyers LA, et al. Neuroinflammatory profiling in SIV-infected Chinese-origin rhesus macaques on antiretroviral therapy. Viruses (2022) 14(1):139. doi: 10.3390/v14010139
136. Brown A, Islam T, Adams R, Nerle S, Kamara M, Eger C, et al. Osteopontin enhances HIV replication and is increased in the brain and cerebrospinal fluid of HIV-infected individuals. J Neurovirol. (2011) 17(4):382–92. doi: 10.1007/s13365-011-0035-4
137. Eger C, Cirelli K, Budiaman J, Brown A. Noncontiguous protein interaction domains in osteopontin contribute to enhance HIV-1 replication. J Hum Virol Retrovirol. (2014) 1(1):7–14. doi: 10.15406/jhvrv.2014.01.00003
138. Marcondes MC, Poling M, Watry DD, Hall D, Fox HS. In vivo osteopontin-induced macrophage accumulation is dependent on CD44 expression. Cell Immunol (2008) 254(1):56–62. doi: 10.1016/j.cellimm.2008.06.012
139. Zhu B, Suzuki K, Goldberg HA, Rittling SR, Denhardt DT, McCulloch CA, et al. Osteopontin modulates CD44-dependent chemotaxis of peritoneal macrophages through G-protein-coupled receptors: evidence of a role for an intracellular form of osteopontin. J Cell Physiol (2004) 198(1):155–67. doi: 10.1002/jcp.10394
140. Lund SA, Wilson CL, Raines EW, Tang J, Giachelli CM, Scatena M. Osteopontin mediates macrophage chemotaxis via alpha4 and alpha9 integrins and survival via the alpha4 integrin. J Cell Biochem (2013) 114(5):1194–202. doi: 10.1002/jcb.24462
141. Chen Y, Colonna M. Microglia in alzheimer's disease at single-cell level. are there common patterns in humans and mice? J Exp Med (2021) 218(9):e20202717. doi: 10.1084/jem.20202717
142. Shen X, Qiu Y, Wight AE, Kim HJ, Cantor H. Definition of a mouse microglial subset that regulates neuronal development and proinflammatory responses in the brain. Proc Natl Acad Sci U S A. (2022) 119(8):e2116241119. doi: 10.1073/pnas.2116241119
143. Shinohara ML, Kim HJ, Kim JH, Garcia VA, Cantor H. Alternative translation of osteopontin generates intracellular and secreted isoforms that mediate distinct biological activities in dendritic cells. Proc Natl Acad Sci U S A. (2008) 105(20):7235–9. doi: 10.1073/pnas.0802301105
144. Inoue M, Shinohara ML. Intracellular osteopontin (iOPN) and immunity. Immunol Res (2011) 49(1-3):160–72. doi: 10.1007/s12026-010-8179-5
145. Yamamoto N, Sakai F, Kon S, Morimoto J, Kimura C, Yamazaki H, et al. Essential role of the cryptic epitope SLAYGLR within osteopontin in a murine model of rheumatoid arthritis. J Clin Invest. (2003) 112(2):181–8. doi: 10.1172/JCI17778
146. Weber GF, Zawaideh S, Hikita S, Kumar VA, Cantor H, Ashkar S. Phosphorylation-dependent interaction of osteopontin with its receptors regulates macrophage migration and activation. J Leukoc Biol (2002) 72(4):752–61. doi: 10.1189/jlb.72.4.752
147. Muzio L, Viotti A, Martino G. Microglia in neuroinflammation and neurodegeneration: From understanding to therapy. Front Neurosci (2021) 15:742065. doi: 10.3389/fnins.2021.742065
148. Kwon HS, Koh SH. Neuroinflammation in neurodegenerative disorders: the roles of microglia and astrocytes. Transl Neurodegener. (2020) 9(1):42. doi: 10.1186/s40035-020-00221-2
149. Maduna T, Audouard E, Dembele D, Mouzaoui N, Reiss D, Massotte D, et al. Microglia express mu opioid receptor: Insights from transcriptomics and fluorescent reporter mice. Front Psychiatry (2018) 9:726. doi: 10.3389/fpsyt.2018.00726
150. Campbell GR, To RK, Spector SA. TREM-1 protects HIV-1-Infected macrophages from apoptosis through maintenance of mitochondrial function. mBio. (2019) 10(6):e02638-19. doi: 10.1128/mBio.02638-19
151. Yuan Z, Fan X, Staitieh B, Bedi C, Spearman P, Guidot DM, et al. HIV-Related proteins prolong macrophage survival through induction of triggering receptor expressed on myeloid cells-1. Sci Rep (2017) 7:42028. doi: 10.1038/srep42028
152. Avey D, Sankararaman S, Yim AKY, Barve R, Milbrandt J, Mitra RD. Single-cell RNA-seq uncovers a robust transcriptional response to morphine by glia. Cell Rep (2018) 24(13):3619–29.e4. doi: 10.1016/j.celrep.2018.08.080
153. Karagiannis TT, Cleary JP Jr., Gok B, Henderson AJ, Martin NG, Yajima M, et al. Single cell transcriptomics reveals opioid usage evokes widespread suppression of antiviral gene program. Nat Commun (2020) 11(1):2611. doi: 10.1038/s41467-020-16159-y
154. Farhadian SF, Mehta SS, Zografou C, Robertson K, Price RW, Pappalardo J, et al. Single-cell RNA sequencing reveals microglia-like cells in cerebrospinal fluid during virologically suppressed HIV. JCI Insight (2018) 3(18):e121718. doi: 10.1172/jci.insight.121718
Keywords: neuroHIV, microglia, HIV, SIV, opioid, opiate, nonhuman primate, central nervous system
Citation: Fox HS, Niu M, Morsey BM, Lamberty BG, Emanuel K, Periyasamy P, Callen S, Acharya A, Kubik G, Eudy J, Guda C, Dyavar SR, Fletcher CV, Byrareddy SN and Buch S (2022) Morphine suppresses peripheral responses and transforms brain myeloid gene expression to favor neuropathogenesis in SIV infection. Front. Immunol. 13:1012884. doi: 10.3389/fimmu.2022.1012884
Received: 06 August 2022; Accepted: 31 October 2022;
Published: 16 November 2022.
Edited by:
Thomas Rogers, Temple University, United StatesReviewed by:
Sylvia Fitting, University of North Carolina at Chapel Hill, United StatesKurt F. Hauser, Virginia Commonwealth University, United States
Copyright © 2022 Fox, Niu, Morsey, Lamberty, Emanuel, Periyasamy, Callen, Acharya, Kubik, Eudy, Guda, Dyavar, Fletcher, Byrareddy and Buch. 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.
*Correspondence: Howard S. Fox, hfox@unmc.edu
†Present address: Shetty Ravi Dyavar, ADICET Bio, Redwood City, CA, United States
‡ Deceased