Monocytic-Myeloid Derived Suppressor Cells Suppress T-Cell Responses in Recovered SARS CoV2-Infected Individuals.

Coronavirus disease 2019 (COVID-19) caused by SARS Coronavirus 2 (CoV2) is associated with massive immune activation and hyperinflammatory response. Acute and severe CoV2 infection is characterized by the expansion of myeloid derived suppressor cells (MDSC) because of cytokine storm, these MDSC suppress T cell functions. However, the presence of MDSC and its effect on CoV2 antigen specific T cell responses in individuals long after first detection of CoV2 and recovery from infection has not been studied. We and others have previously shown that CD11b+CD33+CD14+HLA-DR-/lo monocytic MDSC (M-MDSC) are present in individuals with clinical recovery from viral infection. In this study, we compared the frequency, functional and transcriptional signatures of M-MDSC isolated from CoV2 infected individuals after 5-months of the first detection of the virus (CoV2+) and who were not infected with CoV2 (CoV2-). Compared to CoV2- individuals, M-MDSC were present in CoV2+ individuals at a higher frequency, the level of M-MDSC correlated with the quantity of IL-6 in the plasma. Compared to CoV2-, increased frequency of PD1+, CD57+ and CX3CR1+ T effector memory (TEM) cell subsets was also present in CoV2+ individuals, but these did not correlate with M-MDSC levels. Furthermore, depleting M-MDSC from peripheral blood mononuclear cells (PBMC) increased T cell cytokine production when cultured with the peptide pools of immune dominant spike glycoprotein (S), membrane (M), and nucleocapsid (N) antigens of CoV2. M-MDSC suppressed CoV2 S- antigen-specific T cell in ROS, Arginase, and TGFβ dependent manner. Our gene expression, RNA-seq and pathway analysis studies further confirm that M-MDSC isolated from CoV2+ individuals are enriched in pathways that regulate both innate and adaptive immune responses, but the genes regulating these functions (HLA-DQA1, HLA-DQB1, HLA-B, NLRP3, IL1β, CXCL2, CXCL1) remained downregulated in M-MDSC isolated from CoV2+ individuals. These results demonstrate that M-MDSC suppresses recall responses to CoV2 antigens long after recovery from infection. Our findings suggest M-MDSC as novel regulators of CoV2 specific T cell responses, and should be considered as target to augment responses to vaccine.

Monocytic MDSC inhibit CoV2 antigen specific T cell responses and maintain a divergent transcriptomic profile long after recovery from infection.

Introduction

The coronavirus disease 2019 (COVID-19) pandemic caused by SARS CoV2 (CoV2) virus remains a global threat and cause of massive worldwide morbidity and mortality (https://www.who.int/emergencies/diseases/novel-coronavirus-2019). Even though the pathogenesis of COVID-19 remains elusive, ranging from asymptomatic to multiorgan failure and death, CoV2 infection is unarguably associated with massive immune activation and hyperinflammatory response because of the cytokine storm (1–6). A decline in T cell numbers (lymphopenia) in severe COVID-19 disease suggests the important role of these cells during CoV2 infection (3, 7–11). Furthermore, an association of early detection of CoV2 specific T cells with mild disease and late induction with severe disease support that T cells are indispensable for viral clearance (12–15). Experimental evidence shows that the induction of T cell responses to multiple immune dominant CoV2 antigens is a critical parameter of antiviral efficacy (12, 16–20). Although the kinetics and duration of cellular immunity in CoV2 infection are highly heterogeneous, CoV2 specific T cells persist for at least 6-8 months after natural infection (17, 21–27). The regulation of T cell responses/effector function during the contraction phase of T cells remains less explored.

Myeloid derived suppressor cells (MDSC) are a heterogeneous population of cells derived from immature myeloid progenitors and have immunosuppression potential (28–34). Although much extensively studied in cancer, increasing evidence demonstrate their importance in microbial infections including viral respiratory infections (29, 32, 35–37). In humans, MDSC expresses ommon myeloid markers (CD11b+CD33+HLA DR-/lo) and, depending on the presence of CD15 or CD14, are divided into granulocytic or monocytic subsets, respectively (28, 38–40). These cells utilize multiple mechanisms such as 1) depleting the microenvironment of important nutrients for T cell activation by expressing high levels of arginase-1 (Arg-1) and indoleamine 2,3 dioxygenase (IDO), 2) producing oxygen species such as nitric oxide (NO), reactive oxygen species (ROS) and peroxynitrite (PNT), 3) producing immunosuppressive cytokines such as TGF-β and IL-10, and 4) inducing expansion of regulatory T-cell (Treg) to suppress adaptive response (41–45). Expansion and function of MDSC is regulated by STAT family of transcription factors primarily STAT3 and STAT5, which induces the expression of anti-apoptotic genes and prevents differentiation of myeloid progenitor cells into mature myeloid cells (28, 42, 46–48). Considerable research, predominantly performed in animal models, has demonstrated inhibition of antitumor and antimicrobial activity by MDSC. Multiple recent studies have also shown activation and expansion of MDSC in severe CoV2 infection result in decreased numbers and impaired function of natural killer and T cells (35, 49–53). However, information regarding the presence of MDSC, and their effect on T cell function in individuals with past CoV2 infection is limited.

In this study, we investigated CD11b+CD33+CD14+HLA-DR-/lo monocytic myeloid derived suppressor cells (M-MDSC) in individuals with a history of CoV2 infection (CoV2+) confirmed by RT PCR and individuals without CoV2 (CoV2-). We found an increased frequency of M-MDSC in the peripheral blood of CoV2+ as compared to CoV2-. Importantly, depleting M-MDSC augmented T cell cytokine production in response to the peptide pool of immune dominant CoV2 antigens. Our RNA sequencing (RNAseq) data reveals differential expression of genes in CoV2+ M-MDSC. These results provide insight that targeting M-MDSC may improve the cellular immune response to CoV2 vaccines.

Materials and Methods

Patient Population

Blood was obtained after receipt of written informed consent from BCG un-vaccinated, HIV-, CoV2– and CoV2+ persons. CoV2- and CoV2+ participants with clinical signs of CoV2 infection (temperature of at least 100.4°F, cough, shortness of breath, chills, sore throat, muscle pain, or new loss of taste or smell) or in close contact with anyone showing these symptoms at the time of sample collection were not enrolled for this study. All studies were conducted in accordance with the Declaration of Helsinki guidelines and approved by the institutional review board of the University of Georgia, Athens.

Cell Isolation and Culture

Peripheral blood mononuclear cells (PBMC) were isolated from freshly obtained blood by Ficoll density centrifugation (GE Healthcare) and cultured in RPMI1640 supplemented with AB human serum. For some experiments, PBMC and MDSC depleted PBMC were cultured with the PepTivator® peptide pools for surface glycoprotein (S), nucleocapsid phosphoprotein (N), or membrane glycoprotein (M) of SARS-CoV2 as per the manufacturer’s instructions (Miltenyi Biotec). Briefly, 1.5 x106 cells were cultured in a flat-bottom 96-well plate in a volume of 150 µl culture media, and peptide pools were added at a final concentration of 0.6 nmol of each peptide. To determine the effect of PD1, PD-L1, galectin, ROS, arginase, and iNOS cells were cultured with the peptide pool of S-protein in the presence or absence of neutralizing (PD1, PD-L1, and galectin) or isotype control antibodies, or chemical inhibitors (ROS, arginase, and iNOS).

Antibodies and Other Reagents

Antibodies used for flow cytometry were Alexa Fluor488-anti-CD11b, PE-Dazzel594-anti-HLA DR, PE/Cy7-anti-CD14, APC-anti-CD33, PerCP-Cy5.5-anti-CD66b, BV421-anti-CD19, BV510/PE/Cy7-anti-CD3, APC-eF780/APC-anti-CD4, PE-Fire700/APC-eF780-anti-CD8, AF488-anti-CD45RA, BV421/BV605-anti-CD62L, AF700-anti-CD27, BV421-anti-CD28, PerCP-Cy5.5-anti-PD1, PE-Dazzel594-anti-CD57, PE-anti-CX3CR1, PE-anti-IL-2; LIVE/DEAD fixable aqua stain dye. For neutralization studies, monoclonal antibodies to PD-1, PD-L1, galectin or isotype (10 µg/ml) control were used (all from Biolegend). Chemical inhibitors used were ROS inhibitor catalase (100 U/mL), arginase inhibitor NG-monomethyl-L-arginineacetate (0.5 mM) (both from Sigma Aldrich), and iNOS inhibitor N( ϖ)-hydroxy-nor-L-arginine (0.5 mM; EMD Millipore Corp.).

Depletion of MDSC From Fresh PBMC

Freshly isolated PBMC were stained for CD14, HLA DR, and LIVE/Dead Aqua stain; Aqua+ dead cells were excluded and CD14+HLA DR-/lo MDSC depleted using flow cytometry were collected separately. MDSC depleted PBMC were cultured for flow cytometry and measurement of cytokines. Total RNA was isolated from MDSC for transcriptomic and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) analysis.

Immunolabeling and Flow Cytometry

Cells were stained for surface markers using respective antibodies and cell staining buffer. For intracellular IL-2, surface-stained cells were fixed and permeablized using a Fixation/Permeabilization kit and anti-IL2 antibody. Flow cytometry was done on Quanteon and data was analyzed using Flowjo. A minimum of 300,000 Aqua- LIVE cells were collected for each sample. Controls for each experiment included unstained cells and fluorescence minus one (FMO).

Quantification of Cytokines

Supernatants collected from PBMC and MDSC depleted PBMC cultures were stored at -80°C. Quantity of IFNγ, TNFα, IL-10, IL-17, IL-12p70 were determined at 24-hours post-stimulation with Peptivator Peptide pools using Bioplex analyte detection kit (Biorad). Plasma levels of IL-6 and IL-8 were determined using IL-6 Quantikine HS (sensitivity 0.09 pg/ml) and IL-8 Quantikine HS (sensitivity 0.4 pg/ml), respectively (both from R&D Systems)

Quantitative Reverse-Transcription Polymerase Chain Reaction (qRT-PCR)

qRT-PCR was done as previously described (54). Total RNA was isolated from sorted CD14+HLA DR-/lo MDSC from CoV- and CoV2+ individuals using TRIzol™ reagent (Thermo Fisher Scientific) according to the manufacturer’s protocol; 250 ng RNA was used for cDNA synthesis using iScript™ cDNA Synthesis Kit (BioRad, USA) according to the manufacturer’s instructions. qRT-PCR was performed using the PowerUp™ SYBR™ Green Master Mix (Thermo Fisher Scientific, USA) and data were acquired with the Applied Biosystems™ StepOne™ Real-Time PCR System (Applied Biosystems/Life Technologies, Germany). The Human acidic ribosomal protein (HuPO) was used as a housekeeping gene. The list of genes and primers used in this study is in . Data were analyzed to calculate the relative quantification of the genes in comparison to the HuPO gene by comparative Ct method (2-ΔCt) (55).

Statistical Analysis

Data are expressed as mean values ± standard error mean (SEM). Paired Student t-tests were used to determine the statistical significance of in vitro experiments. Comparisons between CoV2- and CoV2+ subjects was made by non-parametric Mann Whitney U test; comparisons between different parameters were analyzed using Spearman correlation. Statistical analysis was performed using Graphpad Prism 9 (La Jolla, CA); p-values of <0.05 were considered statistically significant.

RNA-Seq Data Generation

CD14+HLA DR-/lo cells from freshly isolated PBMC were isolated after excluding dead cells by flow cytometry. RNA was isolated using RNeasy Mini Kit (QIAGEN, Germany) and quality was assessed using Bioanalyzer (Agilent Technologies, USA). All samples were checked for RNA integrity number (RIN). The mean RIN for 10 samples was 6.9 with a standard deviation 2.9. RNA-Seq libraries were prepared in two batches using 10 pg -1 ng total RNA and SMART-Seq HT PLUS kit kit (Takara Bio USA, Inc.). All libraries were sequenced to the depth of average of 34294989 (range from 16694996 to 46416998) using Illumina 500 mid-output PE75 sequencer at the Georgia Genomics andBioinformatics Core (GGBC), University of Georgia, Athens.

Differential Gene Expression and Pathway Analysis

Data were available from the GGBC, University of Georgia, Athens Core in the.fastq format. Data pre-processing included concatenating reads for one sample obtained on different lanes of the sequencer into one file. Filtering low quality reads and removal of the 3’ adapter sequences were further performed using the Trim Galore tool, which utilizes the Catadapt program (56). Reads were mapped to the latest version of the human genome hg38 (GRCh38.p13) using HISAT2 (57). Mapped reads were counted against the human GENCODE annotation (v37) (58) using HT-Seq (59). The EdgeR library (60) in the R computing environment was used for quality control of the RNA-Seq data, and ComBat-seq method (61) for correction of batch effects. Differential gene expression analysis was conducted using EdgeR. EdgeR uses empirical Bayes estimation and exact tests based on the negative binomial distribution of the RNA-Seq data, followed by false discovery rate (FDR) correction using the Benjamini-Hochberg method (62). Genes were considered differentially expressed when FDR-corrected p-values were <0.05, or nominal p-value <0.05 in an exploratory approach. Pathways over-represented for DEGs were identified using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 (63, 64) and KEGG (65, 66) and Reactome (67) databases. For pathway analysis, all genes with a nominal p-value <0.05 were used as input. DAVID uses a one-tailed Fisher’s Exact test for gene set enrichment analysis. Pathways with nominal p-values <0.1 were considered significantly enriched for DEGs.

Results

Patient Characteristics

All the participants were unvaccinated for COVID-19. CoV2- participants self-reported their COVID-19 status; CoV2+ participants who self-reported positive RT PCR in past were included in the study. Participants were enrolled at a median of 163 days (5.4 months) from the first detection of CoV2 infection ( ). Our participant enrollment form included self-report for Type 2 diabetes (T2D), hypertension, tuberculosis, smoking and chewing tobacco, any respiratory complications, or any other ongoing medical condition. Of all the CoV2- and CoV2+ participants, 2 CoV2- individuals and 1 CoV2+ reported T2D and/or hypertension, and 1 CoV2+ participant self-reported suffering from anxiety. None of the participants reported any other co-morbidity.

Table 1

Demographic and clinical variables in study particpants.

	CoV2- (N = 9)	CoV2+ (N = 13)	p-value
Age in years Median (Range)	49 (24-61)	33 (22-51)	0.59
Gender distribution (Number)  Female  Male	45	68
Time to CoV2+ (Days) Median (Range)	-	163 (18-313)	-

Increased Frequencies of M-MDSC in CoV2+ Individuals

Severe COVID-19 is associated with increased circulating MDSC, which declines upon clinical recovery from the disease (35, 49, 50, 52). Here we sought to investigate the presence of M-MDSC at a median of 163 days (approx. 5 months) from the first detection of CoV2 infection by RT PCR. We found an increased frequency of M-MDSC in peripheral blood of CoV2+ individuals as compared to CoV2- controls (0.2±0.09 vs 0.9±0.2; p=0.03) ( ). Of note, the M-MDSC frequency in peripheral blood of one CoV2- participant who self-reported as hypertensive was relatively higher (0.4%) as compared to other CoV2- participants. We have previously shown that M- MDSC expansion during viral infection is dependent on IL-6 (32, 68). Additionally, COVID-19 is associated with elevated levels of IL-6. Therefore, we quantified IL-6 in the plasma of CoV2- and CoV2+ individuals. Despite a long time after recovery from CoV2 infection, we found higher levels of IL-6 in plasma of CoV2+ individuals, compared to CoV2- controls (0.76±0.2 vs 1.3±0.3; p=0.04) ( ). Furthermore, a positive correlation was found between the quantity of plasma IL-6 and the circulating frequency of MDSC (Pearson r= 0.62; p=0.02) ( ). Of note, IL-8 was also elevated in CoV2+ individuals, but its level did not correlate with M-MDSC ( ). The plasma IL-6 quantity of the hypertensive CoV2- participant was not higher than the remaining participants.

Figure 1

M-MDSC expansion in CoV2+ individuals is dependent on IL-6: (A–C) Heparinized blood obtained from CoV2- and CoV2+ individuals was stained with anti-CD3, -CD19, -CD11b, -CD33, -CD14, CD66b, -HLA DR antibodies, cells were analyzed as CD3-CD19-CD66b-CD11b+CD33+CD14+HLA DR-/lo by flow cytometry. (A) Gating strategy for MDSC is shown (B) A representative dot plot with HLA DR-/lo region from CoV2(-)and CoV2(+)is shown. (C) Percentages of MDSC are shown. (D) The quantity of cytokine IL-6 in the plasma of CoV2- and CoV2+ individuals was measured by ELISA, as in Methods. (E) Plasma IL-6 quantity was correlated with the circulating frequency of M-MDSC in CoV2+ individuals. (C, D) Each dot in the plots depicts data of each individual donor, the plots include observations from 25th to 75th percentile. The horizontal line represents the median value. (D) Each dot in the plot depicts data of each individual donor; black solid and dotted lines, model-estimated values, and their 95% confidence intervals. *p < 0.05.

Phenotypic Evidence of T Cell Exhaustion and Senescence in CoV2+ Individuals

Effector memory T (TEM) cells remain long-term after an infection is eliminated. These T cell subsets are critical to eradicating virus by their ability to produce anti-viral cytokines to control viral replication. Apart from the expression of lymph node homing receptors CCR7 and CD62L, CD27 and CD28 along with CD45RA can be used to discriminate naïve and TEM cells with CD45RA+CD27+CD28+ (naïve T cells), CD45RA-CD62L-CD27+CD28+ (TEM) (69–72). Memory T cells with senescent and exhaustion phenotype have been reported during severe CoV2 infection, displaying a reduced capacity for antiviral cytokine production (73–75). In this research, we sought to determine if phenotypic T cell abnormalities persist 5 months from the first detection of CoV2 infection. For this, initially, we compared the frequency of various CD45RA-CD62L- TEM subsets present in the CD4+ and CD8+ compartments of the whole blood of CoV2- and CoV2+ individuals. We found a comparable frequency of CD27-CD28+, CD27+CD28+, CD27+CD28-, and CD27-CD28- TEM cell subset in the CD4+ and CD8+ compartments of the two groups ( , and , ). The hypertensive CoV2- participant with a higher frequency of MDSC also exhibited elevated frequencies of CD4+ and CD8+ CD27-CD28+ cells.

Figure 2

CoV2+ individuals exhibit an aberrant T cell phenotype 5 months after the first detection of infection: Heparinized blood obtained from CoV2- and CoV2+ individuals was stained with anti-CD3, -CD4, -CD8, -CD45RA, -CD62L, CD27, -CD28, -CD57, -PD-1 and -CX3CR1 antibodies. (A–C) Cells were analyzed for various T memory cell subsets as indicated. (A) Gating strategy for T cells is shown (B, C) Percentages of T memory cells are shown. (D–F) Percentages of PD-1+, CD57+ and CX3CR1+ cells in CD3+CD4+CD45RA-CD62L-CD27+CD28+ T memory cell subset (G–I) Percentages of PD-1+, CD57+ and CX3CR1+ cells in CD3+CD8+CD45RA-CD62L-CD27+CD28+ T memory cell subset is shown. For all, each dot in the plots depicts data of each individual donor, the plots include observations from 25th to 75th percentile. The horizontal line represents the median value. *p < 0.05; **p < 0.005 NS, Non-significant.

Table 2

Percentage of CD4 T effector memory (TEM) cell subsets in the peripheral blood of study participants.

	CoV2- (N = 9)	CoV2+ (N = 13)	p-value
CD27-CD28+	33.6 ± 7.7	26.8 ± 5.9	0.2
CD27+CD28-	12.6 ± 3.4	8.3 ± 1.8	0.5
CD27-CD28-	8.0 ± 2.6	7.2 ± 1.8	0.9
CD27+CD28+ PD1+ CD57+ CX3CR1+	52.3 ± 4.9	57.7 ± 4.9	0.4
	15.2 ± 6.4	26.5 ± 3.5	0.05
	13.8 ± 1.6	30.6 ± 6	0.036
	10.7 ± 1.7	19.7 ± 3.3	0.05

Statistical significant values are italicized.

Table 3

Percentage of CD8 T effector memory (TEM) cell subsets in the peripheral blood of study participants.

	CoV2- (N = 9)	CoV2+ (N = 13)	p-value
CD27-CD28+	24.4 ± 6.5	22.2 ± 5.3	0.7
CD27+CD28-	14.2 ± 3	13.7 ± 2.7	0.9
CD27-CD28-	8.9 ± 3	11.5 ± 2.5	0.7
CD27+CD28+ PD1+ CD57+ CX3CR1+	53.2 ± 4.8	52.6 ± 5.5	0.9
	10.9 ± 4.6	25.9 ± 3	0.005
	18 ± 3.5	30.9 ± 5.7	0.1
	12.4 ± 1.8	21.6 ± 2.9	0.02

Statistical significant values are italicized.

Of the four TEM subsets identified, the frequency of CD27+CD28+ was highest in CD4+ and CD8+ T cells. Therefore subsequently, we compared PD1+, CD57+, and CX3CR1+ cells in these TEM subsets of CD4+ and CD8+ compartments. As compared to CoV2- individuals, CoV2+ individuals exhibited increased frequency of CD4+CD57+ TEM cells (13.8±1.6 vs 30.6±6; p=0.04). Even though the frequency of CD4+PD1+ and CD4+CX3CR1+ TEM cells was also higher in CoV2+ individuals, it was less stringent than CD57+ cell frequency ( and ). In contrast to CD4+ TEM cells, CoV2+ individuals exhibited increased frequency of CD8+PD1+ and CD8+CX3CR1+ cells (10.9±4.6 vs 25.9±3; p=0.005 and 12.4±1.8 vs 21.6±2.9; p=0.02, respectively) ( and ). Of note, the frequency of CD4+CD57+ TEM cells, and CD8+PD1+ or CD8+CX3CR1+ TEM cells did not correlate with the circulating MDSC (Data not shown). Collectively these studies suggest persistence of T cells linked to the suppression of their response for long-duration post CoV2 infection.

M-MDSC Regulates CoV2 Specific T Cell Cytokine Production

We found increased MDSC five months post-CoV2 infection ( ). In the next set of experiments, we sought to investigate if these MDSC regulate T cell function in response to CoV2 antigens. For this, we depleted CD14+HLA DRlo/- MDSC from freshly isolated PBMC of CoV2- and CoV2+ individuals, and stimulated whole PBMC and MDSC depleted PBMC with peptide pools of S, N, and M antigens for 16-18 hours. The frequency of CD3+CD45RA-CD62L- CD4+ and CD8+ TEM cells was determined by flow cytometry. We have previously established that depleting CD14+HLA DRlo/- MDSC does not change CD3+ T cell percentage or expression of CD69, CD38, and HLA DR T cell activation markers in MDSC depleted PBMC fraction (76). In this research, we found that the net CoV2 specific frequency of CD3+CD45RA-CD62L-CD4+IL-2+ was greater in MDSC depleted PBMC cultures when compared to whole PBMC cultures of individuals with past CoV2 infection (p-values 0.007 for S peptide pool, 0.03 for N and M peptide pools) ( ). Similarly, the net CoV2 specific frequency of CD3+CD45RA-CD62L-CD8+IL-2+ was greater in MDSC depleted PBMC cultures when compared to whole PBMC cultures of individuals with past CoV2 infection (p-values 0.008 for S peptide pool, 0.03 for N and M peptide pools) ( ). The cultures of CoV2- individuals did not produce significant IL-2 in response to CoV2 antigens (data not shown).

Figure 3

M-MDSC regulates CoV2 antigen-specific IL-2 production: Freshly isolated PBMC from CoV2+ individuals were stained with anti-CD14 and -HLA DR antibodies; CD14+HLA DR-/lo M-MDSC were depleted from PBMC by flow cytometry. Whole PBMC (PBMC) and MDSC depleted PBMC (PBMC-MDSC) were cultured in the absence or presence of peptide pools of S-, N-, and M- antigens of CoV2 for 20-24 hours. Cells were stained with anti-CD3, -CD4, -CD8, -CD45RA, -CD62L, IL-2 antibodies, and LIVE/DEAD fixable stain. (A, B) Percentages of CD3+CD4+CD45RA-CD62L-IL-2+ cells was determined. (C, D) Percentages of CD3+CD8+CD45RA-CD62L-IL-2+ cells was determined. (A, C) Representative dot plots showing CD4+IL-2+ gated on Live CD3+CD4+CD45RA-CD62L- (A), and CD8+IL-2+ gated on Live CD3+CD8+CD45RA-CD62L- (C) are shown. (B, D) Each dot represents an individual donor. *p < 0.05, **p < 0.005, FMO, Fluorescence minus one.

Subsequently, we quantified T helper cell (TH) cytokines, cytokines that activate T cells, and immune suppressive cytokine IL-10 in the culture supernatants of whole PBMC and MDSC depleted PBMC stimulated for 48-72 hours with peptide pool of CoV2 antigen. The net IFNγ (p-values 0.004 for S and N peptide pools, 0.008 for M peptide pool), TNFα (p-values 0.004 for S, N, and M peptide pools), IL-17 (p-values 0.03 for S, N, and M peptide pools), and IL-12p70 (p-values 0.04 for S peptide pool, 0.03 N and M peptide pools) produced by MDSC depleted PBMC cultures was greater when compared to whole PBMC cultures ( ). In contrast, MDSC depleted PBMC cultures when stimulated with CoV2 antigen peptides produced less IL-10 when compared to whole PBMC cultures stimulated identically (p-values 0.03 for S peptide pool, 0.008 for N and M peptide pools) ( ). These data establish that CD14+HLA DRlo/- MDSC regulates the immunologic response to CoV2 and suppresses TH1 responses even after 5-6 months from initial exposure to the virus. The cultures of CoV2- individuals did not produce significant IFNγ in response to CoV2 antigens ( ).

Figure 4

M-MDSC regulates CoV2 antigen-specific T cell cytokine production: (A–E) Freshly isolated PBMC from CoV2+ individuals were stained with anti-CD14 and -HLA DR antibodies; CD14+HLA DR-/lo M-MDSC were depleted from PBMC by flow cytometry. Whole PBMC (PBMC) and MDSC depleted PBMC (PBMC-MDSC) were cultured in the absence or presence of peptide pools of S-, N-, and M- antigens of CoV2 for 48-72 hours. Culture supernatant was stored at -80°C until further use. The quantity of cytokines in the culture supernatant was measured using Bio-Plex multiplex immunoassay system. For all, each dot represents an individual donor. *p < 0.05, **p < 0.005.

M-MDSC Mediated T Cell Inhibition Is Dependent on Arg and ROS

MDSC utilize multiple mechanisms to confer their suppressor activity. These include biochemical and metabolic mediators ROS, iNOS, and arginase or cell surface inhibitory receptors PD-1, PD-L1, and α-galectin (28, 31, 35, 42, 43, 47). Here we sought to investigate if these mediators regulate T cell response, and delineate which of these mechanisms is utilized by MDSC to control T cell function in CoV2+ individuals. For this, initially, we depleted CD14+HLA DRlo/- MDSC from freshly isolated PBMC of CoV2+ individuals, and cultured whole PBMC and MDSC depleted PBMC with peptide pool of S-antigen in the presence or absence of blocking antibodies for PD-1, PD-L1, or galectin-9 or iNOS, arginase, or ROS inhibitor. The quantity of IFNγ was quantified in the culture supernatants at 48-72 hours by ELISA. Compared to unstimulated controls, PBMC cultured with S-antigen produced more IFNγ (undetectable vs 3.8 ± 1.4 ng/ml; p=0.0001). This was further increased upon blocking PD-1 or PD-L1 (3.8 ± 1.4 vs 4.8 ± 1.5 ng/ml; p=0.01 and 3.8 ± 1.4 vs 4.9 ± 1.6 ng/ml; p=0.02, respectively), or inhibiting iNOS, arginase, or ROS (3.8 ± 1.4 vs 5.6 ± 2.1 ng/ml; p=0.05, 3.8 ± 1.4 vs 5.5 ± 1.9 ng/ml; p=0.05, and 3.8 ± 1.4 vs 6.5 ± 1.4 ng/ml; p=0.04). Of note, blocking galectin-9 did not affect IFNγ production (3.8 ± 1.4 vs 3.5 ± 1.5 ng/ml; p=0.59) ( ).

Figure 5

M-MDSC mediated T cell inhibition is dependent on Arg and ROS: Freshly isolated PBMC from CoV2- and CoV2+ individuals were stained with anti-CD14 and -HLA DR antibodies; CD14+HLA DR-/lo M-MDSC were depleted from PBMC by flow cytometry. (A) Whole PBMC (PBMC) and MDSC depleted PBMC (PBMC-MDSC) of CoV2+ individuals were cultured in the absence or presence of the peptide pool of S-antigen in presence of blocking anti-PD1, -PD-L1, -galectin antibodies or isotype-matched control antibody, or in presence of iNOS inhibitor (L-NMMA), arginase inhibitor (nor-NOHA) or ROS inhibitor (Catalase). Supernatants were collected at 48-72 hours and the quantity of IFNγ was measured by ELISA. (B, C) RNA was isolated from sorted CD14+HLA DR-/lo MDSC using Trizol. The expression of (B) p47 (p47) and (C) Arg2 was determined by qRT PCR as detailed in Methods. For (A) Mean values +/- SEM are shown from 4 separate donors; (B, C) each dot in the plots depicts data of each individual donor, the plots include observations from 25th to 75th percentile. The horizontal line represents the median value. *p < 0.05, **p < 0.005; ***p < 0.0005; NS, Non-significant.

As previously observed ( ) compared to whole PBMC, IFNγ production was more in the culture supernatants of MDSC depleted PBMC cultures (3.8 ± 1.4 vs 11.8 ± 3 ng/ml; p=0.04). To explore the relative contribution of cell surface inhibitory receptors and biochemical mediators on the function of CoV2+ MDSC, blocking antibodies or biochemical inhibitors were added to MDSC depleted PBMC cultures, and the quantity of IFNγ was compared to MDSC depleted PBMC cultured without these inhibitors. We hypothesized that functionally active pathways in MDSC of CoV2+ individuals that may inhibit IFNγ production will be absent in MDSC depleted PBMC cultures. Therefore, the addition of antibodies or chemical inhibitors will have no effect on the IFNγ production in response to CoV2 antigen. However, similar to the whole PBMC culture, blocking PD-1, PD-L1, or inhibiting iNOS also increased IFNγ production in MDSC depleted PBMC cultures (11.8 ± 3 vs 13.4 ± 3.4 ng/ml; p=0.03, 11.8 ± 3 vs 14.5 ± 2.9 ng/ml; p=0.05, and 11.8 ± 3 vs 13.6 ± 3.4 ng/ml; p=0.03), suggesting that none of these pathways contribute to MDSC mediated T cell inhibition. Of note, inhibiting arginase or ROS did not affect IFNγ production (11.8 ± 3 vs 9.8 ± 4.7; p=0.41, and 11.8 ± 3 vs 10.6 ± 3.8; p=0.42, respectively) ( ), suggesting that M-MDSC are primarily responsible for arginase and ROS mediated IFNγ inhibition in PBMC cultures.

Concomitantly, we determined the gene expression of PD-1, PD-L1, NOS2, Arg 1, Arg 2, and p47 ox subunit of ROS producing enzyme NADPH oxidase in MDSC isolated from CoV2+ and CoV2- individuals. Compared to CoV2- controls, MDSC isolated from CoV2+ individuals exhibited an increased expression of p47 and Arg 2 (11.4 ± 3.7 vs 168.7 ± 37.8; p=0.009, and 0.2 ± 0.04 vs 2.5 ± 0.9; p=0.03, respectively) ( ). Of note, the expression of PD-1, PD-L1, NOS2, and Arg 1 remained undetectable. These results corroborated with our ELISA results. Specifically, neutralizing PD-1, PD-L1, and iNOS increased, whereas the addition of arginase or ROS inhibitors to MDSC depleted PBMC cultures do not affect IFNγ production. Collectively, these findings suggest that: (1) T cells utilize multiple mechanisms to control TH1 responses, and (2) MDSC utilizes arginase and ROS pathways to inhibit T cell function in CoV2+ individuals.

CoV2 Infection Induces Long-Term Changes in Gene Expression Profiles

To determine whether infection with CoV2 has a long-term effect on gene expression in M-MDSC, we conducted an RNA-Seq study using samples from ten study participants (5 CoV2- and 5 CoV2+), that were profiled in two separate batches. Our quality control procedures identified a strong batch effect in the data ( ). Following batch effect correction and normalization, two outlier samples, one in each batch were identified, both of which from the CoV2- group. None of these outliers reported any co-morbidities. These outliers were removed based on their relative log expression distributions, resulting in a dataset of eight samples with a similar distribution pattern ( ), which was used for the subsequent differential gene expression analysis. Results from the differential gene expression analysis are shown in .

First, we were interested in assessing the expression of genes quantified in a group of independent study participants by qRT PCR ( , ). Gene expression results of qRT PCR demonstrate significantly increased expression of p47, Arg2 ( ), and TGFβ ( ) in M-MDSC isolated from CoV2+ as compared to CoV2- participants. Additionally, we considered expression levels of transcription factors STAT1, STAT3, STAT5b, HIF-1α, and Nrf2 reported to regulate the expansion and function of MDSC. The median expression of all the genes was greater in M-MDSC isolated from CoV2+ participants, but the differences did not achieve statistical significance ( ). Of note, the expression of STAT3 in M-MDSC of CoV2+ when compared to CoV2- participants approached significance (0.4 ± 0.3 vs 15.9 ± 10.6; p=0.07). Collectively, our findings of positive correlation between IL-6 and frequency of M-MDSC ( ) and STAT3 expression ( ) suggest that the IL-6-STAT3 axis drives M-MDSC expansion in CoV2 infection. Contrary to the results obtained using qRT PCR, differences in the expression of genes between CoV2- and CoV2+ samples did not achieve statistical significance in the RNA-seq study. Nevertheless, some trends in expression remained the same, for example, p47, Arg2 and TGFβ were upregulated when quantified by both qRT PCR and RNA-seq ( ), and PD-1, PD-L1, NOS2 and Arg1 were undetectable by either method.

Figure 6

Expression of genes associated with expansion and function of M-MDSC: Freshly isolated PBMC from CoV2- and CoV2+ individuals were stained with anti-CD14 and -HLA DR antibodies; CD14+HLA DR-/lo M-MDSC were sorted by flow cytometry and stored in Trizol. (A–F) RNA was isolated and the expression of indicated genes was determined by qRT PCR as detailed in Methods. (G–I) Expression of indicated genes was plotted from RNA-Seq data, as counts per million (cpm). For all, each dot in the plots depicts data of each individual donor, the plots include observations from 25th to 75th percentile. The horizontal line represents the median value. p-values are shown, **p < 0.005; nominal p-values are indicated for the RNA-Seq data.

We next took the advantage of the RNA-seq data set to explore additional long-term effects of CoV2 infection on gene expression profiles. This analysis identified 12 new DEGs (FDR <0.05), of which 2 were upregulated and 10 were downregulated ( ); and a total of 188 DEGs with nominal p-value cut off of 0.05, of which 63 were up- and 125 downregulated ( ). Remarkably, despite prolonged time after detected CoV2 (median 163 days/5 months), the absolute fold changes in the expression of most of these genes were quite high, 2- to 378-fold ( ). Many of these genes were previously implicated in severe CoV2 infection; including some genes that did not reach statistical significance in our study. We show that: 1) some of these genes remained upregulated 5 months after infection (LYZ, PLAC8 [fold increase 1.6; p=0.07)], 2) were upregulated during severe disease and downregulated 5 months after recovery (PFKB3, CXCL8 (fold decrease 1.93; p=0.06), TREML4, THBS1, EREG), or 3) remained downregulated long-term (IL1B, NLRP3) (77–81). Our RNA-seq findings combined with cytokine studies ( , ) suggest that after recovery from COVID-19 disease, M-MDSC may contribute to control cytokine induced hyper inflammation, but retain their capacity to suppress TH1 cell function, critical for the anti-viral immune mechanism.

Figure 7

CoV2 infection induces long-term changes in global gene expression profiles: Differentially expressed genes in M-MDSC isolated from CoV2- and CoV2+ individuals were identified using EdgeR. (A) Heatmap was constructed using function heatmap.2 from the library gplots in R. For the heatmap construction, samples were arranged by group (Cov2- followed by CoV2+), and genes were arranged by the average fold change in expression CoV2+/CoV2- (greatest to least); construction of dendrograms was omitted. Gene expression was scaled by row; scale bar represents scaled gene expression values. (B) Volcano plot was constructed using EnhancedVolcano function from the library EnhancedVolcano in R. The input for the volcano plot consisted of gene symbols along with logFC and FDR from the EdgeR differential gene expression analysis output. Vertical dash lines represent the fold change cut off (greater than 1.5 or less than -1.5); horizontal dash line represents the FDR cut off (less than 0.05). (C) Heatmaps were constructed for genes (nominal p < 0.05) from selected pathways (indicated on the left of the heatmaps). Samples were arranged by group (Cov2- followed by CoV2+). Genes were arranged by pathway; genes mapped to two pathways were added to the heatmap once at their overlap.

Immune Response Pathways Are Enriched for Differentially Expressed Genes (DEGs) Modulated Long-Term After CoV2 Infection

We next explored the long-term effects of CoV2 infection on biological processes by conducting pathway enrichment analysis using KEGG and the Reactome databases. In total, this analysis identified 38 pathways significantly enriched for DEGs (nominal p-value < 0.1, ). Of these, ten pathways are directly related to the regulation of immune function, in particular regulating T cell responses ( ). Of note, the genes responsible for upregulated T cell function such as HLA-B, HLA DQB1, and innate immune signaling EREG, IL1B, NLRP3 remained downregulated in M-MDSC of CoV2+ participants ( ). These findings suggest M-MDSC retains the capacity to regulate anti-CoV2 immune function long after CoV2 infection.

Table 4

Pathways enriched for DEGs identified between CoV2+ and CoV2- groups (nominal p-value < 0.05).

Category	Term	Genes	p-value
KEGG_PATHWAY	hsa04612:Antigen processing and presentation	HSPA8, CD74, HLA-B, PSME2, HLA-DQA1, CTSB, HLA-DQB1	0.0002
REACTOME_PATHWAY	R-HSA-114608:R-HSA-114608:Platelet degranulation	TGFB1, HSPA5, SERPING1, TAGLN2, PPBP, THBS1, F5	0.0019
KEGG_PATHWAY	hsa04621:NOD-like receptor signaling pathway	RIPK2, IL1B, NLRP3, CXCL1, CXCL2	0.0033
KEGG_PATHWAY	hsa04668:TNF signaling pathway	RPS6KA4, IL1B, CXCL1, PTGS2, CXCL2, CREB5	0.0064
KEGG_PATHWAY	hsa04145:Phagosome	FCGR3A, TFRC, HLA-B, THBS1, HLA-DQA1, HLA-DQB1	0.0247
KEGG_PATHWAY	hsa05164:Influenza A	IVNS1ABP, HSPA8, IL1B, NLRP3, HLA-DQA1, HLA-DQB1	0.0428
REACTOME_PATHWAY	R-HSA-1257604:R-HSA-1257604:PIP3 activates AKT signaling	IRS2, NRG1, EREG, HBEGF	0.0475
KEGG_PATHWAY	hsa05168:Herpes simplex infection	CD74, EEF1D, IL1B, HLA-B, HLA-DQA1, HLA-DQB1	0.0512
REACTOME_PATHWAY	R-HSA-877300:R-HSA-877300:Interferon gamma signaling	HLA-B, TRIM21, HLA-DQA1, HLA-DQB1	0.0582
REACTOME_PATHWAY	R-HSA-202424:R-HSA-202424:Downstream TCR signaling	RIPK2, PSME2, HLA-DQA1, HLA-DQB1	0.0807

Pathway enrichment analysis was conducted using the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 and KEGG and Reactome databases. Pathways with nominal p-value < 0.1 were considered significant. Pathways identified from both databases were combined and sorted by the nominal p-value.

Discussion

Over the past couple of years, CoV2 infection has been a cause of mortality worldwide. With the development of various vaccines, hospitalization is reduced significantly in vaccinated individuals, yet the factors that may influence prolonged vaccine efficacy remain ill-defined (82–86). In this research, we demonstrate that elevated frequency of M-MDSC persists in peripheral blood of CoV2+ individuals even after 5-months post infection, as a result of persistently elevated levels of IL-6. We have also shown that these M-MDSC are functionally active, as depleting them increased the T cell cytokine production in response to CoV2 peptide pools. We show that M-MDSC suppresses T cell cytokine production through mechanisms involving ROS and arginase. Using the transcriptomic and pathway analysis approach, we further show M-MDSC from CoV2+ individuals are enriched in pathways that regulate both innate and adaptive immune control mechanisms. Importantly, the genes in these pathways remain downregulated long after recovery from COVID-19. These support the previous findings that MDSC regulates CD4 and CD8 T cell responses in CoV2 infection (35, 49). To our knowledge, this is the first demonstration of suppression of T cell recall responses by M-MDSC even after recovery from infection.

MDSC expand during various pathological conditions including microbial infection, due to acute or chronic inflammation. Independent studies on laboratory-confirmed CoV2 infected individuals with varying disease severity establish a higher frequency of both CD14+ M- and CD15+ G-MDSC in individuals with severe disease (35, 49–52). Of note, circulating MDSC inversely correlate with T cell count suggesting the contribution of these cells to T cell dysfunction (49, 50). In the present study, we provide evidence of M-MDSC mediated suppression of CoV2 immunogenic antigen-specific T cell responses long after the first exposure to the virus. COVID-19 is associated with massive immune activation, with elevated levels of proinflammatory cytokine IL-6. Therapeutic modulation of IL-6 levels by anti-IL-6 receptor monoclonal antibodies (mAbs) (tocilizumab, sarilumab) and anti-IL6 mAbs reduces the duration and/or severity of COVID-19 (87–90). Longitudinal studies on proinflammatory cytokine trajectory in COVID-19 remain unexplored, in this study, we found elevated levels of IL-6 and IL-8 in CoV2+ individuals even after 5-months of the first detection of infection. However, only IL-6 correlated with the frequency of M-MDSC. To our knowledge, this is the first demonstration of a direct correlation of plasma IL-6 level with circulating frequency of M-MDSC in CoV2+ individuals.

Although the mechanism of IL-8 (CXCL8) mediated M-MDSC expansion is not established, IL-8 facilitates differentiation and subsequent mobilization of MDSC in patients diagnosed with several types of cancers (91, 92). IL-6, through gp130/JAK/STAT pathway, phosphorylates STAT3, which translocates to the nucleus and regulates differentiation, mobilization, and survival of M-MDSC in human and mice disease conditions including COVID-19 (39, 41, 46, 93). An increased gene expression of STAT3 in M-MDSC of CoV2+ individuals suggests that IL-6-STAT3 is the critical mediator of M-MDSC expansion following CoV2 infection. In our study, STAT3 expression did not reach significance; this could be due to 1) a small sample size, and 2) lower STAT3 expression in CoV2+ individuals exhibiting a low level of M-MDSC ( ). Severe COVID-19 causes a dysregulated myeloid cell compartment resulting in an increase in dysfunctional HLA DRloS100Ahi monocytes cells resembling MDSC (80). In contrast, our transcriptomic data showed a trend of lower S100A9/A8 expression in CD14+HLA DRlo M-MDSC isolated 5 months after the first detection of CoV2 ( ). Even though the markers differentiating G- and M-MDSC in COVID-19 infection are not well established, S100A9/A8 proteins are majorly expressed by G-MDSC, and our transcriptomic analysis did not show the expression of CECAM8 (CD66b), a marker exclusively expressed on G-MDSC (28, 31, 44). It is also possible that the downregulation of S100A proteins in our study is due to the control in hyper inflammation over a period of time post-infection. Our findings resemble the transcriptional signature of M-MDSC reported by Kvedaraite et al, demonstrating an increased expression of IFITM2, IFITM3, and lower expression of CD83 ( ) in cells isolated from CoV2+ individuals as compared to CoV2- controls (94). We show here that the transcriptional signatures of M-MDSC and T cell suppression due to M-MDSC is maintained as long as after 5 months of the first detection of the infection. It appears that aberrant myelopoiesis continues after recovery from COVID-19, and IL-8 and IL-6 remain important mediators for the mobilization and expansion of STAT3 expressing M-MDSC.

STAT3 directly regulates the p47phox subunit of ROS producing NADPH oxidase (NOX2) complex (42, 95). Our finding that M-MDSC isolated from CoV2+ individuals express more p47, and ROS neutralization in M-MDSC depleted cultures did not augment IFNγ production in response to the peptide pool of immune dominant S-antigen is similar to other viral infections (96, 97). Efficient T cell signaling and subsequent immune effector function require proper interaction between T cell receptors (TCR) on APC and CD3 on T cells. ROS contributes to the generation of peroxynitrite that causes nitration of a tyrosine residue at TCR-CD3 complex in CD8+ cells, thereby weakening CD8-TCR interaction (98, 99). Whether ROS in M-MDSC of CoV2+ individuals also modifies tyrosine or any other amino acid at the TCR-T cell junction remains to be explored. Additionally, amino acid L-arginine is essential for T cell activation and function. Its depletion in the cell culture leads to a rapid decrease in CD3ζ levels, and this is reversed by L-arginine supplementation in vitro or arginase (ARG) inhibition when co-culture with ARG-producing cells is used (100–103). The two isoforms ARG1 (encoded by ARG1) and ARG2 (encoded by ARG2) of the enzyme ARG hydrolyze L-arginine (104–106). We found increased expression of ARG2 gene in M-MDSC of CoV2+ individuals, and ARG inhibition by nor-NOHA in M-MDSC depleted cultures did not affect IFNγ production support the findings of Falck-Jones et al. demonstrating lower expression of CD3ζ-chain in CD4+ and CD8+ T cells of COVID-19 patients, and an increased T cell function in PBMC-M-MDSC co-cultures supplemented with L-arginine (35). However, our transcriptomic analysis ( ) did not reveal the expression of ARG1 gene, and ARG1 expression in M-MDSC measured by qRT PCR remained undetectable suggesting a dominant role of ARG2 dependent depletion of L-arginine by M-MDSC in our study cohort. A similar prominent role of ARG2 in macrophages has been observed in Helicobacter pylori infection, and CD8+ T cells from ARG2-deficient mice show superior antitumor activity (107, 108). It is highly like that similar to human dendritic cells and other myeloid cell types, both ARG1 and ARG2 contribute to the immune suppressive function of M-MDSC in CoV2 infection; nonetheless, the relative contribution of the two isoforms to deplete L-arginine might be dependent on the disease stage (109, 110). Here we propose that ROS and Arginase produced by M-MDSC are the mediators of CoV2 antigen-specific T cell suppression.Increasing evidence suggest that both CD4 and CD8 T cell responses are important for COVID-19 outcome and maintenance of CoV2 immunity, even in the absence of humoral response. T cell responses to CoV2 has been shown to recognize epitopes across multiple viral proteins. Similar to the findings of Saini et al. and Tarke et al, we also found that CoV2+ individuals produced significant quantities of TH1 cytokines in response to the peptide pools of S, N and M proteins (20, 111). Other studies have shown that CoV2 S, N and M, together with non-structural protein (nsp)3, nsp4 and ORF3 are recognized by T cells (12, 19, 20). Although epitope mapping was beyond the scope of this study, and our sample size was small to establish immune dominant antigen for T cell responses in our cohort, we observed cultures stimulated with peptide pool of N antigen produced the maximum quantities of TH1 cytokine, and the T cell response was durable up to 5-months. Independent studies have shown that the durability and stability of CoV2 specific T cell responses are maintained from 3 to 12 months post infection (17). Our findings followed the trend previously observed to corresponding peptide pools, and that T cell recognition of multiple epitopes commonly occurs in CoV2 infection. Both G- and M- MDSC have been shown to suppress T cell proliferation and cytokine production during severe and mild CoV2 infection (35, 49, 52). Previous studies of Agrati et al. and Falck-Jones et al, show that depleting MDSC from PBMC of CoV2+ individuals increased the T cell cytokine production in response to Staphylococcus enterotoxin B (SEB) polyclonal stimulation, thus establishing that M-MDSC isolated from these individuals had a potent suppressive effect on T cells (35, 49). Here we show M-MDSC inhibits T cell cytokine production in response to CoV2 antigen, thus in an antigen specific manner. Our findings and those of others also show that MDSC isolated from CoV2+ individuals suppress T cell cytokine production via arginase and TGFβ dependent mechanisms. Importantly, we provide the evidence that M-MDSC identified in CoV2+ individuals long after recovery from infection are suppressive and functionally active. To our knowledge, this is the first report showing CoV2 antigen-specific immune suppressive activity of M-MDSC, and at a much later time point.

During chronic viral infection, persistent antigen disrupts memory cell development and leading to impairment in T cell function characterized by upregulation of exhaustion (PD-1, TIM-3) and senescence markers (CD57) (69, 112). Similar to other chronic viral infections we observed an increase in CD8+ TEM cells expressing PD-1, and CD4+ TEM cells expressing CD57 in the peripheral blood of CoV2+ individuals. However, neither these correlate to the frequency of M-MDSC in peripheral blood, nor does depleting M-MDSC from PBMC affect TEM cell subsets (data not shown). Based on our blocking studies ( ) and RNA-seq data ( ), in this study, we provide evidence PD-1-related inhibitory pathways do not contribute to M-MDSC mediated immune suppressive function of T cells in CoV2 infection. Even though the expression of PD-1 during the acute and early convalescent phase of CoV2 infection correlates more with activation state rather than functional exhaustion, future studies on understanding the role of these T cell subsets after clinical recovery from the infection will be important (18, 113). Additionally, COVID-19 vaccines are also administered to COVID-19 convalescent individuals. Given that, exhausted T cells lose their potential to differentiate into memory T cells, T cell exhaustion in individuals with past CoV2 infection can impede vaccine-induced T cell memory development.

Antigen presentation is a vital immune process essential for triggering T cell response and is a fundamental element of host defense (114). Classical antigen-presenting cells are of myeloid (macrophages, monocytes, dendritic cells) origin, which in association with MHC class II (for exogenous antigens) and MHC class I (for endogenous antigens) molecules activate CD4+ and CD8+ T cells, respectively (114). We show here that M-MDSC isolated from CoV2+ individuals suppress T cell cytokine production in an antigen-specific manner. Although we studied IL-2 production, but not other cytokines exclusively produced by CD4+ and CD8+ T cells, in line with this data, we found the genes HLA DQB1, HLA DQA1, and HLA A critical for antigen presentation are down-regulated in M-MDSC isolated from all CoV2+ individuals. Among the various genes categorized in the KEGG enriched antigen processing and presentation pathway, CD74, HSPA8, and PSME2 genes involved in the processing of exogenous antigens exhibited a variable expression pattern among CoV2+ participants (115). Similarly, PSME2 implicated in immunoproteasome assembly and required for efficient antigen processing of endogenous antigens also exhibited variable expression among CoV2+ participants (116). It appears that other participant-dependent factors may regulate antigen processing in M-MDSC, but inefficient antigen presentation is the dominant factor driving suppressed T cell responses in CoV2+ participants. Another pathway in the KEGG database related to pathogen/antigen uptake and its presentation to T cells is Phagosome. Although the expression of FCGR3A responsible for FcR mediated phagocytosis was variable among the CoV2+ participants, genes coding for proteins associated with Phagolysosome Biogenesis (THBS1) and microtubule activity of phagosome pathway (TFRC) along with HLA DQB1, HLA DQA1, and HLA A were consistently down-regulated (117, 118). Together, these data suggest that CoV2 infection induces long-term changes that negatively influence anti-CoV2 T cell response.

Additionally, compared to CoV2-, multiple pathways related to innate immune signaling were also enriched in M-MDSC of Cov2+ individuals. However, the genes driving these pathways were majorly downregulated. IL-1β-NLRP3 axis has a double-sword function in anti-viral immunity which may facilitate viral eradication but also results in the proinflammatory cell death known as pyroptosis (119). Chemokines CXCL1 and CXCL2 direct the recruitment of various immune effector cells and also activate the NLRP3 inflammasome (120). PI3K/AKT is one of the most critical pathways in innate immune signaling and cell survival. Here we found that genes IRS2, NRG1, HBEGF, and EREG regulating AKT signaling are downregulated in M-MDSC isolated from CoV2+ participants long after recovery from infection (121). These data suggest that during recovery from infection M-MDSC may have developed a mechanism to countermeasure hyper inflammation by inhibiting innate immune signaling, but retain their capacity to suppress T cell responses.

Limitations of Study

The pathological and immunological consequences of M-MDSC after recovery from CoV2 remain unclear. Given the immune regulatory role of M-MDSC in various infections and immunological abnormalities, it is highly likely that M-MDSC in individuals with critical COVID-19 illness contributes to immunosuppression (35, 49–52). While our study design provides robust and reproducible results suggesting an expansion of M-MDSC even after 5-months of the first report of infection, which exhibits the dual function of suppressing T cells responses and downregulated inflammatory responses, the small sample size in our study is a major limitation and reason for some findings unable to reach statistical significance. Thus, other mechanisms driving the response of M-MDSC, such as genetics, comorbidities, age, gender, or initial viral load remain unexplored. Due to the increased rate of vaccination in this community, the identification of non-vaccinated participants with a history of CoV2 infection, and non-vaccinated uninfected controls became a challenge. Thus this study was performed on a small cohort of non-vaccinated CoV2- and CoV2+ individuals. We performed RNA-seq studies in two batches to compare the global gene profile of M-MDSC in CoV2- and CoV2+ individuals; however, we observed a huge batch effect ( , ) with less coverage of reads across human transcriptome in Batch 2. This resulted in the identification of limited numbers of DEGs (12 DEGs with FDR <0.05) in M-MDSC from CoV2+ individuals, and the fold change expression of many genes (CSTA, SELL (CD62L), S100A8/A9, FCN1, DUSP6) important for characterization and function of M-MDSC did not reach statistical significance. Additionally, two CoV2- samples were outliers and excluded from the analysis. Despite these limitations, our study clearly demonstrates that M-MDSC retains their capacity to suppress T cell function even after recovery from CoV2 infection. Indeed, future studies performed on large cohorts of individuals receiving the COVID-19 vaccine will be important to establish the effect of M-MDSC on long-term vaccine efficacy.

Conclusion

Collectively, our ex vivo and in vitro data establish that CoV2 infected individuals have an increased quantity of M-MDSC even after 5-months of the first report of infection. Importantly, we show that M-MDSC controls T cell recall responses to peptides of major antigenic proteins of the virus. In addition, we show transcriptional changes in M-MDSC of CoV2 + individuals as compared to those of CoV2- individuals. Our findings suggest that M-MDSC are important regulators of immune responses, and provide attractive targets to augment vaccine responses. Our ongoing studies are centered to study differences in M-MDSC population in response to various CoV2 vaccines. An improved understanding of MDSC biology in CoV2 infection can contribute to novel approaches to prolong durability of vaccine response and/or restore normal immune function.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The name of the repository and accession number can be found below: (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE199286).

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional review board of the University of Georgia, Athens. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

AG developed the study, performed experiments and analyzed the data; KM performed the experiments; RU and LS performed the preliminary analysis of RNA-seq data; NBB analyzed RNA-seq data; AG and NBB wrote the manuscript.

Funding

AG was supported in part by AI127132 from the National Institute of Allergy and Infectious Diseases (NIAID), IUSSTF/VN-COVID/011/2020 from INDO-US Science and Technology Forum, and The University of Georgia Research Foundation. A part of this work was supported by the National Center for Advancing Translational Sciences of the National Institutes of Health under Award Number UL1TR002378. NB-B was supported in part by 1I01 BX005285 from the Office of Research and Development, Veterans Health Administration and R56 AI157755 (NIH/NIAID). LS and RX were supported by LS’s startup funds from the University of Georgia Office of Research.

Author Disclaimer

The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

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.