- 1Cardiovascular Research Center, Lewis Katz School of Medicine at Temple University, Philadelphia, PA, United States
- 2Center for Inflammation and Lung Research, Lewis Katz School of Medicine at Temple University, Philadelphia, PA, United States
- 3Center for Metabolic Disease Research, Lewis Katz School of Medicine at Temple University, Philadelphia, PA, United States
- 4Department of Thoracic Medicine and Surgery, Lewis Katz School of Medicine at Temple University, Philadelphia, PA, United States
There is a high incidence of tobacco use among intravenous opioid drug users. It is well established that opioids and tobacco smoke induce a degree of immune activation, and recent work suggests that the combination of these drugs promotes further activation of the immune system. Our approach involved the treatment of wild-type mice with cigarette smoke (SM) for a period of eight weeks, and the chronic continuous administration of morphine (M) via mini-pumps for the final four weeks. In an effort to examine the responses of CD4+CD25highCD127low regulatory T (Treg) cells, the major immune suppressive cell type, to the combined chronic administration of SM and M, we determined the frequency of these cells in the spleen, lymph nodes and lungs. Flow cytometric analyses showed that SM and M individually, and the combination (SM + M) have differential effects on the numbers of Treg in the spleen, lymph node, and lung. Either SM or M alone increased Treg cell numbers in the spleen, but SM+M did not. Furthermore, SM + M decreased Treg cell numbers in the lymph node and lung. We then performed RNA-Seq on Treg cells from mice treated with SM, M, or SM + M, and we found that the S + M induced a number of significant changes in the transcriptome, that were not as apparent following treatment with either SM or M alone. This included an activation of TWEAK, PI3K/AKT and OXPHOS pathways and a shift to Th17 immunity. Our results have provided novel insights on tissue Treg cell changes, which we suggest are the result of transcriptomic reprogramming induced by SM, M, and SM + M, respectively. We believe these results may lead to the identification of novel therapeutic targets for suppressing smoke and opioid induced Treg cell impairment.
Introduction
Tobacco use is a major cause of death from cancer, cardiovascular disease, and pulmonary disease (1) mediated in part via 11 oncogenic pathways (2), signals from 29 particulate elements (3), a variety of DNA adducts (4, 5), and more than 5,000 chemicals (6, 7). In chronic obstructive pulmonary disease (COPD) patients, higher numbers of Treg cells have been found in the pulmonary lymphocyte follicles but not in the lung parenchyma (8). Furthermore, Foxp3+ Treg are upregulated in large airways but downregulated in small airways in patients with COPD which correlates with limited airflow (9, 10). SM extracts promote DNA methylatransferase 3a expression in dendritic cells, resulting in increased Th17 cell numbers and fewer Treg cells (11). In the early development of COPD in animal models, exposure to SM decreases the numbers of signal transducer and activator of transcription 5 (STAT5) and phospho-STAT5 (pSTAT5) expressing cells and the expression levels of Treg cytokines transforming growth factor- β (TGF-β) and interleukin-10 (IL-10). However SM exposure later increases the numbers of STAT3+ and pSTAT3+ cells and Th17 cytokine IL-17 (12). We previously reported that the number of Treg cells in the blood of patients with severe COPD was significantly reduced (13), suggesting that chronic inflammation may be promoted by the lack of Treg cells. However, the molecular basis for the change in Treg numbers at the level of the transcriptome remains unknown.
In addition, people who use illicit opiate drugs are reported to have high rates of tobacco use. In a 24-year follow-up of narcotics addicts, the death rate of smokers were four times than of non‐smokers (14, 15). Morphine, a breakdown product from heroin, has been found to impair innate immune responses, T cell activation, and shift toward type 2 T helper cell (Th2) differentiation with increased Th1 cell death (16). We have reported that opioids upregulate the expression of several important pro-inflammatory chemokines as well as the chemokine receptor C-C Motif Chemokine Receptor 5 (CCR5) on T cells, monocytes and macrophages (17). In addition, morphine can increase IL-10 levels (18) and concomitantly reduce IL-17 secretion (19) from cultured CD4+ T cells (20). However, molecular mechanisms underlying morphine effects on Treg homeostasis and function remain poorly defined.
In this study, we report the Treg changes in lung, lymph node and spleen, and accompanying transcriptomic reprogramming of Treg cells upon the administration of SM, M, and the combination of SM with M. The combination of SM plus M decreases Treg cells in lymph node and lung potentially by modulating cytokine expression, and cell death pathways. Since Treg cells are the major immunosuppressive cell type and play essential roles in maintaining immune homeostasis, and inhibiting inflammation, autoimmune diseases and anti-tumor immune response, novel therapeutic strategies could be developed to maintain Treg survival (21, 22).
Materials and Methods
Animals
B6.129S4-Arg1tm1.1Lky/J mice were purchased from the Jackson Laboratory (Stock No: 015857, Bar Harbor, ME). All mice were weaned at 3 weeks of age and maintained on chow diet until 8 weeks old. All animal experiments were performed in accordance with the Institutional Animal Care and Use Committee (IACUC) Guidelines and Authorization for the use of Laboratory Animals and were approved by the IACUC of Temple University Lewis Katz School of Medicine.
Cigarette Smoke Administration
SM from filtered 3R4F Cigarettes (University of Kentucky, Lexington, KY) was delivered using a TED-10 smoking system (Teaque Enterprises, Woodland, CA) using standard parameters (23). Mice were exposed in whole body chambers to mainstream SM as previously described (15). Mice were exposed to SM for two hours, 5 days a week for a total of 8 weeks. Control mice, not exposed to cigarette smoke, were exposed to filtered room air in identical whole body chambers for the same time period. An average total particulate matter (TPM) for the first week was 150 mg/day, and this was elevated to 225 mg/day beginning with the second week and for the remainder of the experiment. This dosing regimen is sufficient to establish an inflammatory response in the lung tissue.
Morphine Administration via Mini−Pumps
Morphine administration was conducted using 28 day Alzet mini-pumps (Durect Corp., Cupertino, CA). These pumps deliver a consistent and continuous dose of morphine over a 28-day-period. The animals received a dose of 1 mg/kg/day.
Cell Preparation and Flow Cytometry
To make single cell suspension from spleen and lymph node (LN), freshly collected spleen and LN were washed with PBS and homogenized gently between the frosted ends of the slides in FACS buffer. The suspended cells were transferred into a 15 ml conical tube and immediately centrifugated at 600 g for 5 minutes at 4°C, then the supernatant was removed. The red blood cells (RBC) in splenic cell pellets were lyzed by adding 1 ml ACK lysis buffer and incubated for 5 minutes, then 8 ml of FACS buffer was added to stop the RBC lysis reaction. Cells were spinned at 600 g for 5 minutes and the resuspended cells were passed through the cell strainer (40 µm) into 50 ml tubes and the cell concentrations were adjusted to 1-5x106 cells/ml.
To make single cell suspension from lung, lung lobs were collected and kept in DMEM medium-low containing 10% fetal bovine serum (FBS) on ice. Next, lung was washed twice with PBS and then micro-dissected into 2 mm pieces and digested in 2 ml PBS containing 20 mM HEPES, 5% FBS, 60 U/ml DNAse1 (Sigma, # 11284932001), and 450 U/ml collagenase type I (Sigma, # SCR103) at 37°C for 45 minutes with shake every 15 minutes. Then the mixture was passed through a 70 µm strainer into 50 ml tube. Cells were spinned at 600 g for 7 minutes and the lung cells were resuspend in in ice-cold FACS buffer.
Cells were stained with anti-CD3, CD4, CD25, CD127, CD8, and CD45 (BD Bioscience, San Jose, CA) for 15 minutes at 4°C, and washed twice with FACS buffer before running on BD™ LSR II flow cytometer. For cytokine detection, isolated splenocytes and lung cells were incubated in RPMI-1640 complete medium (supplemented with 10% FBS, 50 μM 2-ME, 10 mM L-glutamine, 10 mM HEPES, 1 mM sodium pyruvate, 100 U/mL pen/strep, and β-mercaptoethanol) at 37°C with cell stimulation cocktail (plus protein transport inhibitors) (#00-4975-93, Invitrogen) for 5 hours. Cells were stained with previously mentioned surface markers followed by intracellular staining with anti-IL-10, and anti-IL-17A. Data was analyzed on Flowjo v10.
CD4+CD25highCD127low Treg Purification and RNA-Seq
Splenocytes were stained with anti-CD4, CD25, and CD127 (BD Bioscience, San Jose, CA) for 15 minutes at 4°C and washed twice with FACS buffer. Cell-sorting experiments were performed using an Aria Cell Sorter (BD Biosciences) at Temple University Lewis Katz School of Medicine Flow Cytometry Core. The CD4+CD25highCD127low cells were sorted directly into TRIzol, and total RNAs from four groups (n = 2) were extracted using miRNeasy Mini kit (Qiagen) following the manufacturer instruction.
RNA-Seq was performed by BGI (North America) (24). Using the fragmentation buffer, the mRNAs were fragmented into short fragments (about 200–500 nucleotides), then the first-strand cDNA was synthesized by random hexamer primer using the fragments as templates, and dTTP was substituted by dUTP during the synthesis of the second strand. Short fragments were purified and resolved with elution buffer for end repair and single nucleotide A (adenine) addition. After that, the short fragments were connected with adaptors, and then the second strand was degraded finally using uracil N-glycosylase (2). After agarose gel electrophoresis, the suitable fragments were selected for polymerase chain reaction (PCR) amplification as templates. During the quality control steps, an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System (Thermo Fisher) were used for quantification and qualification of the sample library. At last, the library was sequenced with an Illumina HiSeq4000 using the PE100 strategy. Primary sequencing data produced by the Illumina Hiseq4000, called raw reads, were filtered into clean reads by removing adaptor containing and low-quality reads by BGI in-house software. A reference annotation–based assembly method was used to reconstruct the transcripts by TopHat (v2.0.10) and Cufflinks (v2.1.1), and background noise was reduced using fragments per kilobase million and coverage threshold.
Statistical Analysis of RNA-Seq Data
Data analysis was carried out using the statistical computing environment R, the Bioconductor suite of packages for R, and RStudio (24). Raw data were background-subtracted, variance-stabilized, and normalized by robust spline normalization. Differentially expressed genes were identified by linear modeling and Bayesian statistics using the Limma package (25). For comparisons between two groups, moderated t test was used for evaluation of statistical significance. All original RNA-seq data were deposited in the NCBI’s Gene Expression Omnibus database (GSE198210).
Ingenuity Pathway Analysis
We utilized Ingenuity Pathway Analysis IPA, Qiagen, https://www.qiagenbioinformatics.com/products/ingenuity-pathway-analysis/) to characterize the clinical relevance and molecular and cellular functions related to the identified genes in our microarray analysis. Differentially expressed genes were identified and uploaded into IPA for analysis. The core and pathways analysis was used to identify molecular and cellular pathways, as we previously reported (26–28). P value < 0.05, and |Z-score|≥1 was set as cut-off in this study. Of note, pathways with |Z-score|≥2 were designated as significantly enfluenced (29).
Results
SM + M Treatment Decreased CD4+CD25highCD127low Treg Cells in the Lung and Lymph Nodes
In an effort to examine the effects of SM and M on the regulatory T cell component of the imimune system, we first profiled T cell subsets by flow cytometry in the lung in animals subjected to SM, or chronic treatment with M, or the combination as described in Figure 1A. The results show (Figure 1B and Suplementery S1) that SM, M and SM + M decreased numbers of CD45+ total leukocytes in the lung compared to the control group but did not change the proportion of T cells in this organ. Further characterization of the T cell compartment, showed that each of the treatments increased the proportion of the CD8+ T cell subset, with a coinciding decrease in the CD4+ subset in the SM and SM + M but not in M group. Next, we examined the CD25highCD127low Treg cell frequency in both subsets and found that SM, M and SM + M did not significantly change Treg in the CD4+ or CD8+ T cell subsets in the lung compared to the control group. However, SM + M decreased CD4+CD25highCD127low Treg cells compared to either M treatment alone or SM treatment.
Figure 1 SM + M treatment decreased CD4+CD25highCD127low Treg cells in the lung and lymph nodes. (A) Schematic representation of stimulations in four groups of male mice: 1) non-stimulation controls (NS); 2) smoking (SM); 3) morphine (M); and 4) smoking plus morphine (SM + M). B-D. The flow cytometry detections were performed for analyzing the populations of CD45 (lung only), CD3, CD4, CD8, CD25highCD127low Treg in (B) lung (n = 4), (C) spleen, and (D) lymph node (LN) (n = 8). Bar graphs indicate the percentage of each subset in its parent population. Percentage of CD25highCD127low Treg in all cells from lung, spleen and LN were calculated and showed in the last row. (T-test, *p < 0.05, **p < 0.01, ***p < 0.001, * in black represents the difference between SM, M or SM+M and NS control, * in red represents the difference between SM+M and SM or M).
We also examined the T cell compartment in spleen and the results show (Figure 1C and Suplementary S1) that SM treatment alone slightly but significantly increased CD8+ subset compared to the control. In addition, SM + M slightly but significantly increased the CD4+ subset compared to the SM treated group. However, no significant changes were observed in the M treated group. Furthermore, we found that SM and M treatments alone significantly increased the percentage of CD25highCD127low Treg cells in in the spleen, in contrast to the SM + M treatment. Moreover, we examined the LN, and here again the SM and SM + M modestly, but significantly, increased the percentage of CD3+ T cells, but decreased the proportionof CD4+ T cell subset with a modest increase in the percentage of the CD8+ subset. In contrast to the lung, there was a significant reduction in the percentage of CD4+CD25highCD127low Treg cells (Figure 1D and Suplementary S1), and both the SM alone and M alone failed to exhibit this effect. As with the lung, the exposure to SM + M caused an decrease in CD8+CD25highCD127low cells in LN.
Taken together, these results have demonstrated that SM, M and SM + M did not change CD4+CD25highCD127low Treg cells in lung but SM + M decreased CD4+CD25highCD127low Treg cells compared to treatments with either SM or M alone. In addition, SM and M did, but SM + M did not, increase CD4+CD25highCD127low Treg cells in spleen. Finally, SM increased, SM + M decreased, but M did not change, the proportion of Treg cells in LN.
Morphine Is Stronger Than Cigarette Smoke in Re-Shaping Treg Cell Transcriptomes, and SM + M Have Synergy in Modulating Treg Cell Transcriptomes
Previuos studies reported that Treg cells are weakened in patients with COPD (10). Cigarette smokers with COPD also have significantly fewer Treg cells and a lower level of Foxp3 mRNA in the lung than healthy smokers (30). We previously reported that chronic M administration increased numbers of circulating Treg cells and functional Th17 cells (31), suggesting that even with increased Treg cell numbers, M may weaken Treg function in suppressing Th17. We hypothesized that SM, M, or SM + M may re-shape Treg cell transcriptomes. To test this hypothesis, we isolated splenic CD4+CD25highCD127low Treg cells and performed RNA-Seq on the isolated cells. Of note, to preserve Treg cell RNA quality, we did not use Treg cell-specific transcription factor Foxp3 as an intracellular staining marker for the Treg. Based on the results (Table 1), we found that SM alone significantly modulated 648 genes with 211 genes upregulated and 437 genes downregulated. The analysis also showed that morphine alone significantly modulated 902 genes with 233 genes upregulated and 669 genes downregulated, and SM + M significantly modulated 1,630 genes with 630 genes upregulated and 1,000 genes downregulated. Moreover, SM + M significantly modulated 1507 genes with 743 genes upregulated and 764 genes downregulated compared to the M treated alone. Finally, SM + M significantly modulated 710 genes with 349 genes upregulated and 361 genes downregulated compared to SM alone.
Table 1 CD4+CD25highCD127low regulatory T cells (Treg) RNA-Seq experiments led to identification of numerous differentially expressed genes (DEGs) in Treg in four groups of mice (n = 2).
These results have demonstrated that first, the 4 week chronic M administration is stronger than the 8 week SM exposure in re-shaping Treg cell transcriptomes, potentially due to the more focused receptor-signaling pathways of M than SM; and second, SM plus M have synergistic effects in modulating Treg transcriptomes.
SM, M, and SM + M Downregulate IL-10 Expression, and SM + M Increase Th17 Transcription Factor (TF) Rorc Expression in Treg Cells
We hypothesize that pathological conditions changes Treg characteristics, result in generating distinct cytokines, interactors/signaling pathways and chemokines. To test this hypothesis, we analyzed changes in cytokine mRNA expression in Treg cells from SM, M and SM + M treated animals. The results show (Figure 2A) that SM, M and SM + M treatments significanly downregulated IL-10 expression in Treg cells compared to the control. In addition, Treg cells from M-treated animals increased expression of the Th1 cytokine Ifng, the Th2 inducing cytokine Il4 (32), chemokine Ccl5 and receptor Ccr5, but decreased Ccr1, Ccr3, and Ccr6 expression. Treg cells from the SM + sM treated animals exhibited lower Ifng, but increased proinflammatory cytokine tumor necrosis factor (Tnf), Ccl1, Ccl9, and notably, increased expression of Th17 transcription factor Rorc.
Figure 2 Morphine (M) is stronger than cigarette smoke (SM) in re-shaping Treg cell transcriptomes and SM + M have synergy in modulating Treg cell transcriptomes. (A) Volcano plots illustrate log 10 (p-value) in relation to the log 2 (fold change) for the treatments versus control (n = 2). The differentially expressed genes with -log10 (p. value) >1.3 and |log2FC| >1 are displayed in red. Cytokines, chemokines and receptors with transcriptomic changes are labeled. (B) Bar graph show the level of transcripts of Foxp3 and cofactors, Treg effector molecules, and immune checkpoint molecules in Treg from four groups. (C) Bar graph show the level of transcripts of 15 nicotinic receptors and 4 opioid receptors in Treg from four groups. TPM, Transcripts per Million mapped reads (T test, *p < 0.05, **p < 0.01, ***p < 0.001, * in black represents the difference between SM, M or SM+M and NS control, * in red represents the difference between SM+M and SM or M).
In adition, we found that Foxp3 and its cofactors Hdac9, Gata1, Ikzf2 (encodes Helios), Ikzf4, Irf4, and Plagl1 were specifically inhibited following treatmen with M (Figure 2B), which was consistent with the previous finding that chronic M administration weakened Treg function (31). Since SM, M and SM + M significantly suppressed IL-10 production, we next investigated the levels of other Treg effector molecules, which are critical for Treg cell function (33). As expected, EGF-like growth factor amphiregulin (Areg), which confers Treg cell suppressive function via the EGFR/GSK-3β/Foxp3 axis, was significantly downregulated in SM, M and SM + M Treg cells. Lrrc32 (encodes glycoprotein-A repetitions predominant, GARP), which is critical for tethering TGF-β to the cell surface, was specifically decreased with M treatment. No change was observed in the expression of Gzmb (encodes granzyme B) and Tgfb1 (encodes TGF-β1). Moreover, we examined the expression of immune checkpoint molecules, and the results show that SM + M increased two co-inhibitory checkpoints Cd96 (93) and Pdcd1 (PD1) and decreased both co-inhibitory Ctla4 and co-stimulatory CD28. Morphine treatment alone decreased three TNF receptor superfamily members, Tnfrsf4 (OX40), Tnfrsf18 (GITR), Tnfrsf9 (4-1BB), but SM + M did not. Finally, SM increased co-stimulatory checkpoint Cd40lg and the induction was absent with the combined administration of SM plus M.
SM, M and SM + M Downregulate the Expression of Nicotinic Acetylcholine Receptor Chrnb1
In order to assess the expression of the nicotinic receptor and opioid receptor on Treg cells, we analyzed the expressions of three muscarinic acetylcholine receptor subunits, 12 nicotinic acetylcholine receptor (nAChR) subunits (34) and four opioid receptors (35) in SM, M, and SM + M treated Treg cells. The results show (Figure 2C) that among the nicotinic receptor subunits examined, nAChR Chrna1, Chrna2, Chrnb1, Chrmb2 have comparably higher expression levels than other subunits, suggesting that these subunits might be the dominant structural conponents on Treg cells. Among these, Chrnb1 was significantly downregulated with SM, M and SM + M treatment. In addition, M treatment specifically decreases the expression of muscarinic acetylcholine receptor Chrm3. On the other hand, analysis of the opioid receptor transcripts shows that Oprm1 exhibits a higher basal level among these receptors on Treg cells, but no significant change were observed between the treatment groups.
SM + M Have Synergistic Effects to Activate TWEAK, PI3K/AKT and OXPHOS Pathways on Treg Cells
In order to have a more comprehensive assessment of the modulation that SM, M and SM + M induced in Treg cells, we used IPA to determine the signaling pathways related to all of the transcriptomic reprogramming of Treg cells treated by SM, M, and SM + M. As shown in Figure 3A, SM + M preferentially activated TWEAK (TNF-like weak inducer of apoptosis) and phosphatidylinositol 3-kinase/protein kinase B (PI3K/AKT) (Z-score =1.964) signaling. In contrast with Treg cells from SM + M treated animals, pathways for IL-2 and IL-6 signaling, Th2 activation, and TREM-1 signaling were all inhibited in M treated Treg cells. SM treatment selectively activated LXR/RXR signaling. Moreover, several significantly modulated pathways were shared by SM + M, SM, and/or M alone. For example, SM + M and M both activated mTOR signaling. SM + M and SM both activated the PPARa/RXRa pathway, and inhibited both PKCθ signaling, and the STAT3 pathway.
Figure 3 SM + M have synergistic effects to activate TWEAK, PI3K/AKT and OXPHOS pathways on Treg cells. (A) Ingenuity Pathway Analysis (IPA) identify the signaling pathways based on the differentially expressed genes (with P value < 0.05) from each comparisons. Bar graph (left) show the top 10 pathways (ordered by |Z-score|) only modulated upon SM, M and SM + M treatment versus control (P value < 0.05, and |Z-score|≥1 was set as cut-off in this study). Plot graph (middle) show overlapped pathways modulated in SM + M and SM or/and M. Bar graph (right) show top pathways modulated upon SM + M versus SM or M alone. Bar in dark orange (Z-score ≥2) represent activation of the pathway and in dark blue (Z-score ≤ -2) represent inhibition. (B) Ingenuity Pathway Analysis (left) identify the metabolic pathways modulated upon treatments. Gene set enrichment assay (GSEA, middle) also identify the expression dataset of SM + M are significantly enriched in the class of genes involved in KEGG OXPHOS pathway. Heatmaps on the right display the gene change associate with the activation of OXPHOS in SM + M when compared with NS or M alone. (NES, normalized enrichment score; FDR, false discovery rate).
Sterol metabolism has been linked to the anti-inflammatory response in CD4+ T cells, therefore, we examined the metabolic pathways significantly affected by SM, M, and SM + M (Figure 3B). Our results showed that SM + M treatment activated oxidative phosphoration (OXPHOS), and superpathway of cholesterol biosynthesis, while suppressed phosphatidylcholine biosynthesis. In addition, using gene set enrichment analysis (GSEA), we also found genes activating OXPHOS were significantly enriched in SM + M Treg when compared to control or M Treg.
Treatment With S+M Stimulated IL-17 Production in CD4+CD25highCD127low Treg Cells From Spleen and Lung
To test the hypothesis that SM, M and SM + M weakened Treg function and SM + M-Treg cells exhibited IL-17-biased plasticity, which was suggested by the RNA-seq results, we determined the numbers of cells producing IL-10 and IL-17A in splenocytes and lung cells using flow cytometry. The results show (Figure 4A) that splenic IL-10 production was significantly decreased in animals following SM and M treatments, but was not significantly altered by SM + M. However, in the splenic CD4+CD25highCD127low Treg subsets, IL-10 production was inhibited not only by SM, and M, but also by SM + M treatment, which was consistent with the RNA-Seq results (Figure 2A). Moreover, we also verified that IL-17A production in the splenic CD4+CD25highCD127low Treg subsets was induced in SM + M Treg cells but not in the Treg cells from SM or M treated mice. Of note, M treatment alone significantly suppressed both IL-10 and IL-17A production in the spleen and in CD3+CD4+ splenocytes.
Figure 4 SM + M stimulated IL-17A production in CD4+CD25highCD127low Treg cells from spleen and lung. (A) Representative flow cytometry gating and quantitative bar graph show the IL-17A and IL-10 production in splenocytes upon treatments. Dot plot indicate the frequency of IL-10 and IL-17A producing cells in CD3+CD4+ T cells and Treg cells. (B) Representative flow cytometry gating and quantitative bar graph show the IL-17A and IL-10 production in lung cells upon treatments. Dot plot indicate the frequency of IL-10 and IL-17A producing cells in CD3+CD4+ T cells and Treg cells. (n = 8) (T-test, *p < 0.05, **p < 0.01,***p < 0.001, * in black represents the difference between SM, M or SM+M and NS control, * in red represents the difference between SM+M and SM or M).
Next, we examined IL-10 and IL-17A cytokine production in lung cells, and found that SM + M treatments significantly increased IL-10 production in lung cells compared to SM alone (Figure 4B). However, elevated levels of IL-10 were not observed in the Treg subsets in the lung. Furthermore, M, but not SM + M, treatment induced IL-17A production in lung cells. However, Treg cells from SM + M treated mice showed a slight but not signicant increase in IL-17 production compared to in the CD4+CD25highCD127low subset from M treated animals. It is worth noting that we did not analyze cells that produce both IL-10 and IL-17A in these experiments because of the low level of double-positive events.
Discussion
Treg cells have a great therapeutic potential for many diseases including systemic lupus erythematosus, organ-specific autoimmune diseases (36) (type I diabetes, psoriasis, myasthenia gravis, inflammatory bowel disease, and multiple sclerosis), transplantation, cancers (37, 38), and cardiovascular disease (39–41). Altered Treg cell numbers have been reported in tobacco use related diseases including increased Treg cell numbers in pulmonary lymphocyte follicles of moderate COPD patients (8), increased Foxp3+ Treg cells in large airways but decreased in small airways (main pathological sites in COPD) (9, 10), and decreased Treg cell numbers in the blood of patients with severe COPD (13), suggesting that the chronic inflammation in these patients may be promoted by lack of this regulatory cell population. In addition, studies revealed that M use alters T cell compartmentation and function, including the Treg cell component (20, 42). Our study observed that both SM and M alone increased the proportion of Treg cells in spleen, and SM increased the percentage of Treg cells in LN. Furthermore, we identified that SM + M decreased Treg cells in lung, spleen, and LN compared to either M alone or SM alone. The greater strength of the combination of SM and M is consistent with our previous studies suggesting that these treatments are more potent inhibitors of Treg cells than either SM or M alone (15).
In an effort to understand the molecular basis for the effects of SM and M, and the combination, we conducted RNAseq analysis of the Treg cells. Our results revealed that these treatments have a number of important influences on the transcripts, and molecular pathways that control Treg cell function. First, while we found that M is potent in re-shaping Treg transcriptomes, the combination of SM + M exhibit synergistic activity in modulating the Treg cell transcriptome. For example, we observed and verified that SM + M, as well as M and SM alone, downregulated IL-10 gene expression and protein production. Given the important role of this cytokine in the function of Treg cells (43–46), these results suggest that SM and M down-regulate both the numbers and functional activity of these important regulatory cells. In addition, we found that SM + M increased Th17 transcription factor Rorc and the IL-17A generation, which suggested that SM + M Treg cells exhibited IL-17-biased plasticity (47).
In addition, data from RNAseq analysis showed that M Treg cells decreased the expression of Foxp3 and its cofactors, suggesting a defective regulatory network, which may accelerate the expression of pro-inflammatory chemokine expression by T cells, monocytes and macrophages (17). In addition, M decreased the expression of co-stimulatory checkpoints Tnfrsf4, Tnfrsf9, Tnfrsf18, but increased Cd40lg in Treg. It has been reported that M inhibited antigen-presenting abilities of dendritic cells and anti-CD40 antibodies could restore the function (48). Taking together, further work is needed to assess the CD40/CD40L interaction between dendritic cells and Treg following chronic M treatment. Meanwhile, while M decreases Ctla-4 expression, SM + M Treg cells exhibited downregulation of both Ctla-4 and CD28, the co-inhibitory and co-stimulatory checkpoint couple that compete to bind to the antigen presenting cell (APC) ligands CD80/CD86. These results suggest a dysfunction in the binding properties of APC ligands CD80/CD86 in SM + M treated animals (49).
The greater strength of the combination of SM and M was also apparent in signaling and metabolic pathways. SM + M activated TWEAK signaling in Treg, which has been reported to amplify inflammation and promote IL-17 secrection (50), leading to tissue damage and potentially impeding endogenous repair mechanisms (51). Treatment with SM + M also resulted in the activation of PI3K-AKT singaling, which is reported to deliver a cell-intrinsic negative signal to restrain Treg cell suppressive activity (52). Moreover, SM + M and M both activated mTOR signaling, which, on the one hand, triggers Treg proliferation, but also weakens the functions and stability of Treg (53, 54). Besides, mTOR also directs Treg cell polarization into other functionally plastic subpopulations (55). Metabolically, SM + M activated the OXPHOS pathway. Compared to the Th subsets, Treg cells are less reliant on glycolysis and use OXPHOS for energy production (56). In general, Treg show a decrease in glycolysis and an increase in OXPHOS, fatty acid oxidation/synthesis during their phase of suppressive function (57). However, a very recent study reported that under Th17 conditions, OXPHOS fine-tunes pathogenic Th17 and Treg cell fate decision and mitochondrial respiration determines the lineage specification of pathogenic Th17 over immunosuppressive Treg cell (58).
In conclusion, S + M treatment downregulated the Treg population in the lung and peripheral lymphoid organs more than cigarrette smoke or morphine, and re-programs Treg transcritptomes via the cell death related TWEAK pathway, PI3K-AKT/mTOR pathway and metabolic-driven OXPHOS pathway. In addition, S + M treatment induced Treg cells to exhibit Th17-biased plasticity (Figure 5). Our findings provide novel insights on the roles of cigarette smoke, morphine and the combination in reprogramming Treg transcriptomes as well as novel targets for the future therapeutic interventions involving immunosuppression, including human immunodeficiency virus (HIV) pathologies, drug abuse pathologies, COPD, cardiovascular diseases, autoimmune diseases, transplantation, cancers and tissue repair.
Figure 5 Working model: Chronic exposure to the combination of cigarette smoke and morphine decreases CD4+ regulatory T cell numbers by reprogramming the Treg cell transcriptome.
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: NCBI, accession ID: GSE198210.
Ethics Statement
The animal study was reviewed and approved by Animal Care and Use Committee (IACUC) of Temple University Lewis Katz School of Medicine.
Author Contributions
YS, WC and AK performed animal works and YS carried out flow cytometry and RNAseq data gathering, data analysis and prepared tables and figures. WC, KX, AK, FS, YL, XJ, GC, HW aided with analysis of the data. XY and TR supervised the experimental design and manuscript writing. All authors read and approved the final manuscript.
Funding
Our research activities are supported by grants from the National Institutes of Health (NIH) (HL131460, HL132399, HL138749, HL147565, P30 DA13429, RO1 DA040619, RO1 DA049745). The content in this article is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
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.887681/full#supplementary-material
References
1. Benowitz NL. Nicotine Addiction. N Engl J Med (2010) 362(24):2295–303. doi: 10.1056/NEJMra0809890
2. Ma Y, Li MD. Establishment of a Strong Link Between Smoking and Cancer Pathogenesis Through DNA Methylation Analysis. Sci Rep (2017) 7(1):1811. doi: 10.1038/s41598-017-01856-4
3. Behera SN, Xian H, Balasubramanian R. Human Health Risk Associated With Exposure to Toxic Elements in Mainstream and Sidestream Cigarette Smoke. Sci Total Environ (2014) 472:947–56. doi: 10.1016/j.scitotenv.2013.11.063
4. Ma B, Stepanov I, Hecht SS. Recent Studies on DNA Adducts Resulting From Human Exposure to Tobacco Smoke. Toxics (2019) 7(1):219–46. doi: 10.3390/toxics7010016
5. Paiano V, Maertens L, Guidolin V, Yang J, Balbo S, Hecht SS. Quantitative Liquid Chromatography-Nanoelectrospray Ionization-High-Resolution Tandem Mass Spectrometry Analysis of Acrolein-DNA Adducts and Etheno-DNA Adducts in Oral Cells From Cigarette Smokers and Nonsmokers. Chem Res Toxicol (2020) 33(8):2197–207. doi: 10.1021/acs.chemrestox.0c00223
6. Talhout R, Schulz T, Florek E, van Benthem J, Wester P, Opperhuizen A. Hazardous Compounds in Tobacco Smoke. Int J Environ Res Public Health (2011) 8(2):613–28. doi: 10.3390/ijerph8020613
7. Turesky RJ, Lu K. Biomarkers of Environmental Toxicants: Exposure and Biological Effects. Toxics (2020) 8(2):1–4. doi: 10.3390/toxics8020037
8. Plumb J, Smyth LJ, Adams HR, Vestbo J, Bentley A, Singh SD. Increased T-Regulatory Cells Within Lymphocyte Follicles in Moderate COPD. Eur Respir J (2009) 34(1):89–94. doi: 10.1183/09031936.00100708
9. Isajevs S, Taivans I, Strazda G, Kopeika U, Bukovskis M, Gordjusina V, et al. Decreased FOXP3 Expression in Small Airways of Smokers With COPD. Eur Respir J (2009) 33(1):61–7. doi: 10.1183/09031936.00145307
10. Hou J, Sun Y. Role of Regulatory T Cells in Disturbed Immune Homeostasis in Patients With Chronic Obstructive Pulmonary Disease. Front Immunol (2020) 11:723. doi: 10.3389/fimmu.2020.00723
11. Huang D, Zhou B, Luo ZZ, Yu SC, Tang B. Cigarette Smoke Extract Promotes DNA Methyltransferase 3a Expression in Dendritic Cells, Inducing Th-17/Treg Imbalance via the C-Jun/allograft Inflammatory Factor 1 Axis. Kaohsiung J Med Sci (2021) 37(7):594–603. doi: 10.1002/kjm2.12367
12. Silva LEF, Lourenco JD, Silva KR, Santana FPR, Kohler JB, Moreira AR, et al. Th17/Treg Imbalance in COPD Development: Suppressors of Cytokine Signaling and Signal Transducers and Activators of Transcription Proteins. Sci Rep (2020) 10(1):15287. doi: 10.1038/s41598-020-72305-y
13. Chatila WM, Criner GJ, Hancock WW, Akimova T, Moldover B, Chang JK, et al. Blunted Expression of miR-199a-5p in Regulatory T Cells of Patients With Chronic Obstructive Pulmonary Disease Compared to Unaffected Smokers. Clin Exp Immunol (2014) 177(1):341–52. doi: 10.1111/cei.12325
14. Hser YI, McCarthy WJ, Anglin MD. Tobacco Use as a Distal Predictor of Mortality Among Long-Term Narcotics Addicts. Prev Med (1994) 23(1):61–9. doi: 10.1006/pmed.1994.1009
15. Cornwell WD, Sriram U, Seliga A, Zuluaga-Ramirez V, Gajghate S, Rom S, et al. Tobacco Smoke and Morphine Alter Peripheral and CNS Inflammation Following HIV Infection in a Humanized Mouse Model. Sci Rep (2020) 10(1):13977. doi: 10.1038/s41598-020-70374-7
16. Plein LM, Rittner HL. Opioids and the Immune System - Friend or Foe. Br J Pharmacol (2018) 175(14):2717–25. doi: 10.1111/bph.13750
17. Steele AD, Henderson EE, Rogers TJ. Mu-Opioid Modulation of HIV-1 Coreceptor Expression and HIV-1 Replication. Virology (2003) 309(1):99–107. doi: 10.1016/S0042-6822(03)00015-1
18. Vianello E, Kalousova M, Dozio E, Tacchini L, Zima T, Corsi Romanelli MM. Osteopontin: The Molecular Bridge Between Fat and Cardiac-Renal Disorders. Int J Mol Sci (2020) 21(15):1–13. doi: 10.3390/ijms21155568
19. Mai J, Nanayakkara G, Lopez-Pastrana J, Li X, Li YF, Wang X, et al. Interleukin-17a Promotes Aortic Endothelial Cell Activation via Transcriptionally and Post-Translationally Activating P38 Mitogen-Activated Protein Kinase (MAPK) Pathway. J Biol Chem (2016) 291(10):4939–54. doi: 10.1074/jbc.M115.690081
20. Abo-Elnazar S, Moaaz M, Ghoneim H, Molokhia T, El-Korany W. Th17/Treg Imbalance in Opioids and Cannabinoids Addiction: Relationship to NF-kappaB Activation in CD4+ T Cells. Egypt J Immunol (2014) 21(2):33–47.
21. Yang Y, Xiong Z, Zhang S, Yan Y, Nguyen J, Ng B, et al. Bcl-xL Inhibits T-Cell Apoptosis Induced by Expression of SARS Coronavirus E Protein in the Absence of Growth Factors. Biochem J (2005) 392(Pt 1):135–43. doi: 10.1042/BJ20050698
22. Yang Y, Yang F, Xiong Z, Yan Y, Wang X, Nishino M, et al. An N-Terminal Region of Translationally Controlled Tumor Protein is Required for its Antiapoptotic Activity. Oncogene (2005) 24(30):4778–88. doi: 10.1038/sj.onc.1208666
23. Lin CR, Bahmed K, Tomar D, Marchetti N, Criner GJ, Bolla S, et al. The Relationship Between DJ-1 and S100A8 in Human Primary Alveolar Type II Cells in Emphysema. Am J Physiol Lung Cell Mol Physiol (2019) 317(6):L791-804.
24. Li X, Wang L, Fang P, Sun Y, Jiang X, Wang H, et al. Lysophospholipids Induce Innate Immune Transdifferentiation of Endothelial Cells, Resulting in Prolonged Endothelial Activation. J Biol Chem (2018) 293(28):11033–45. jbc.RA118.002752. doi: 10.1074/jbc.RA118.002752
25. Beiting DP, Hidano S, Baggs JE, Geskes JM, Fang Q, Wherry EJ, et al. The Orphan Nuclear Receptor TLX Is an Enhancer of STAT1-Mediated Transcription and Immunity to Toxoplasma Gondii. PloS Biol (2015) 13(7):e1002200. doi: 10.1371/journal.pbio.1002200
26. Wang L, Fu H, Nanayakkara G, Li Y, Shao Y, Johnson C, et al. Novel Extracellular and Nuclear Caspase-1 and Inflammasomes Propagate Inflammation and Regulate Gene Expression: A Comprehensive Database Mining Study. J Hematol Oncol (2016) 9(1):122. doi: 10.1186/s13045-016-0351-5
27. Li YF, Nanayakkara G, Sun Y, Li X, Wang L, Cueto R, et al. Analyses of Caspase-1-Regulated Transcriptomes in Various Tissues Lead to Identification of Novel IL-1beta-, IL-18- and Sirtuin-1-Independent Pathways. J Hematol Oncol (2017) 10(1):40. doi: 10.1186/s13045-017-0406-2
28. Zhang R, Saredy J, Shao Y, Yao T, Liu Lu, Saaoud F, et al. End-Stage Renal Disease is Different From Chronic Kidney Disease in Upregulating ROS-Modulated Proinflammatory Secretome in PBMCs - A Novel Multiple-Hit Model for Disease Progression. Redox Biol (2020) 34:178–203. doi: 10.1016/j.redox.2020.101460
29. Kramer A, Green J, Pollard J Jr., Tugendreich S. Causal Analysis Approaches in Ingenuity Pathway Analysis. Bioinformatics (2014) 30(4):523–30. doi: 10.1093/bioinformatics/btt703
30. Lee SH, Goswami S, Grudo A, Song LZ, Bandi V, Goodnight-White S, et al. Antielastin Autoimmunity in Tobacco Smoking-Induced Emphysema. Nat Med (2007) 13(5):567–9. doi: 10.1038/nm1583
31. Cornwell WD, Lewis MG, Fan X, Rappaport J, Rogers TJ. Effect of Chronic Morphine Administration on Circulating T Cell Population Dynamics in Rhesus Macaques. J Neuroimmunol (2013) 265(1-2):43–50. doi: 10.1016/j.jneuroim.2013.09.013
32. Xu K, Yang WY, Nanayakkara GK, Shao Y, Yang F, Hu W, et al. GATA3, HDAC6, and BCL6 Regulate FOXP3+ Treg Plasticity and Determine Treg Conversion Into Either Novel Antigen-Presenting Cell-Like Treg or Th1-Treg. Front Immunol (2018) 9:45. doi: 10.3389/fimmu.2018.00045
33. Zemmour D, Zilionis R, Kiner E, Klein AM, Mathis D, Benoist C. Single-Cell Gene Expression Reveals a Landscape of Regulatory T Cell Phenotypes Shaped by the TCR. Nat Immunol (2018) 19(3):291–301. doi: 10.1038/s41590-018-0051-0
34. Friedman JR, Richbart SD, Merritt JC, Brown KC, Nolan NA, Akers AT, et al. Acetylcholine Signaling System in Progression of Lung Cancers. Pharmacol Ther (2019) 194:222–54. doi: 10.1016/j.pharmthera.2018.10.002
35. Waldhoer M, Bartlett SE, Whistler JL. Opioid Receptors. Annu Rev Biochem (2004) 73:953–90. doi: 10.1146/annurev.biochem.73.011303.073940
36. Yang XF, Fang P, Meng S, Jan M, Xiong X, Yin Y, et al. The Forkhead Transcription Factors Play Important Roles in Vascular Pathology and Immunology. Adv Exp Med Biol (2009) 665:90–105. doi: 10.1007/978-1-4419-1599-3_7
37. Sharabi A, Tsokos MG, Ding Y, Malek TR, Klatzmann D, Tsokos GC. Regulatory T Cells in the Treatment of Disease. Nat Rev Drug Discov (2018) 17(11):823–44. doi: 10.1038/nrd.2018.148
38. Togashi Y, Shitara K, Nishikawa H. Regulatory T Cells in Cancer Immunosuppression - Implications for Anticancer Therapy. Nat Rev Clin Oncol (2019) 16(6):356–71. doi: 10.1038/s41571-019-0175-7
39. Meng X, Yang J, Dong M, Zhang K, Tu E, Gao Q, et al. Regulatory T Cells in Cardiovascular Diseases. Nat Rev Cardiol (2016) 13(3):167–79. doi: 10.1038/nrcardio.2015.169
40. Pastrana JL, Sha X, Virtue A, Mai J, Cueto R, Lee IA, et al. Regulatory T Cells and Atherosclerosis. J Clin Exp Cardiolog (2012) 2012(Suppl 12):2. doi: 10.4172/2155-9880.S12-002
41. Xiong Z, Yan Y, Song J, Fang P, Yin Y, Yang Y, et al. Expression of TCTP Antisense in CD25(high) Regulatory T Cells Aggravates Cuff-Injured Vascular Inflammation. Atherosclerosis (2009) 203(2):401–8. doi: 10.1016/j.atherosclerosis.2008.07.041
42. Peng Y, Yang J, Guo D, Zheng C, Sun H, Zhang Q, et al. Sufentanil Postoperative Analgesia Reduce the Increase of T Helper 17 (Th17) Cells and FoxP3(+) Regulatory T (Treg) Cells in Rat Hepatocellular Carcinoma Surgical Model: A Randomised Animal Study. BMC Anesthesiol (2020) 20(1):212. doi: 10.1186/s12871-020-01129-0
43. Li X, Shao Y, Sha X, Fang P, Kuo YM, Andrews AJ, et al. IL-35 (Interleukin-35) Suppresses Endothelial Cell Activation by Inhibiting Mitochondrial Reactive Oxygen Species-Mediated Site-Specific Acetylation of H3K14 (Histone 3 Lysine 14). Arterioscler Thromb Vasc Biol (2018) 38(3):599–609. doi: 10.1161/ATVBAHA.117.310626
44. Shao Y, Cheng Z, Li X, Chernaya V, Wang H, Yang XF. Immunosuppressive/anti-Inflammatory Cytokines Directly and Indirectly Inhibit Endothelial Dysfunction–a Novel Mechanism for Maintaining Vascular Function. J Hematol Oncol (2014) 7:80. doi: 10.1186/s13045-014-0080-6
45. Li X, Fang P, Sun Y, Shao Y, Yang WY, Jiang X, et al. Anti-Inflammatory Cytokines IL-35 and IL-10 Block Atherogenic Lysophosphatidylcholine-Induced, Mitochondrial ROS-Mediated Innate Immune Activation, But Spare Innate Immune Memory Signature in Endothelial Cells. Redox Biol (2020) 28:101373. doi: 10.1016/j.redox.2019.101373
46. Shao Y, Yang WY, Saaoud F, Drummer Ct, Sun Y, Xu K, et al. IL-35 Promotes CD4+Foxp3+ Tregs and Inhibits Atherosclerosis via Maintaining CCR5-Amplified Treg-Suppressive Mechanisms. JCI Insight (2021) 6(19):127-37. doi: 10.1172/jci.insight.152511
47. Sakaguchi S, Vignali DA, Rudensky AY, Niec RE, Waldmann H. The Plasticity and Stability of Regulatory T Cells. Nat Rev Immunol (2013) 13(6):461–7. doi: 10.1038/nri3464
48. Chang MC, Chen YL, Chiang YC, Cheng YJ, Jen YW, Chen CA, et al. Anti-CD40 Antibody and Toll-Like Receptor 3 Ligand Restore Dendritic Cell-Mediated Anti-Tumor Immunity Suppressed by Morphine. Am J Cancer Res (2016) 6(2):157–72.
49. Shen H, Wu N, Nanayakkara G, Fu H, Yang Q, Yang WY, et al. Co-Signaling Receptors Regulate T-Cell Plasticity and Immune Tolerance. Front Biosci (Landmark Ed) (2019) 24:96–132. doi: 10.2741/4710
50. Park JS, Park MK, Lee SY, Oh HJ, Lim MA, Cho WT, et al. TWEAK Promotes the Production of Interleukin-17 in Rheumatoid Arthritis. Cytokine (2012) 60(1):143–9. doi: 10.1016/j.cyto.2012.06.285
51. Kumar M, Makonchuk DY, Li H, Mittal A, Kumar A. TNF-Like Weak Inducer of Apoptosis (TWEAK) Activates Proinflammatory Signaling Pathways and Gene Expression Through the Activation of TGF-Beta-Activated Kinase 1. J Immunol (2009) 182(4):2439–48. doi: 10.4049/jimmunol.0803357
52. Sauer S, Bruno L, Hertweck A, Finlay D, Leleu M, Spivakov M, et al. T Cell Receptor Signaling Controls Foxp3 Expression via PI3K, Akt, and mTOR. Proc Natl Acad Sci USA (2008) 105(22):7797–802. doi: 10.1073/pnas.0800928105
53. Butcher MJ, Filipowicz AR, Waseem TC, McGary CM, Crow KJ, Magilnick N, et al. Atherosclerosis-Driven Treg Plasticity Results in Formation of a Dysfunctional Subset of Plastic IFNgamma+ Th1/Tregs. Circ Res (2016) 119(11):1190–203. doi: 10.1161/CIRCRESAHA.116.309764
54. Li J, McArdle S, Gholami A, Kimura T, Wolf D, Gerhardt T, et al. Ccr5+T-Bet+FoxP3+ Effector CD4 T Cells Drive Atherosclerosis. Circ Res (2016) 118(10):1540–52. doi: 10.1161/CIRCRESAHA.116.308648
55. Shi H, Chi H. Metabolic Control of Treg Cell Stability, Plasticity, and Tissue-Specific Heterogeneity. Front Immunol (2019) 10:2716. doi: 10.3389/fimmu.2019.02716
56. Michalek RD, Gerriets VA, Jacobs SR, Macintyre AN, MacIver NJ, Mason EF, et al. Cutting Edge: Distinct Glycolytic and Lipid Oxidative Metabolic Programs Are Essential for Effector and Regulatory CD4+ T Cell Subsets. J Immunol (2011) 186(6):3299–303. doi: 10.4049/jimmunol.1003613
57. Kempkes RWM, Joosten I, Koenen H, He X. Metabolic Pathways Involved in Regulatory T Cell Functionality. Front Immunol (2019) 10:2839. doi: 10.3389/fimmu.2019.02839
Keywords: CD4+CD25highCD127low regulatory T cells (Treg), cigarette smoke, morphine, smoke plus morphine, RNA-Seq transcriptomes
Citation: Shao Y, Cornwell W, Xu K, Kirchhoff A, Saasoud F, Lu Y, Jiang X, Criner GJ, Wang H, Rogers TJ and Yang X (2022) Chronic Exposure to the Combination of Cigarette Smoke and Morphine Decreases CD4+ Regulatory T Cell Numbers by Reprogramming the Treg Cell Transcriptome. Front. Immunol. 13:887681. doi: 10.3389/fimmu.2022.887681
Received: 01 March 2022; Accepted: 24 March 2022;
Published: 20 April 2022.
Edited by:
Junji Xing, Houston Methodist Research Institute, United StatesReviewed by:
Mei-Zhen Cui, University of Texas of the Permian Basin, United StatesDevendra K Agrawal, Western University of Health Sciences, United States
Copyright © 2022 Shao, Cornwell, Xu, Kirchhoff, Saasoud, Lu, Jiang, Criner, Wang, Rogers and Yang. 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: Xiaofeng Yang, eGZ5YW5nQHRlbXBsZS5lZHU=; Thomas J. Rogers, dGhvbWFzLnJvZ2Vyc0B0ZW1wbGUuZWR1