Abstract
COVID-19 has caused more than 7 million deaths worldwide, and according to the World Health Organization, it continues to result in more than 1000 reported deaths per month at the time of this writing. It is crucial to understand the immune response to COVID-19 since the virus continues to persist. Natural killer (NK) cells play a critical role in the immune defense against viral infections, including COVID-19. While it is well documented that infected patients have a reduction in lymphocytes and NK cells, gaps in knowledge exist regarding the function of NK cells. To study the function of NK cells in patients hospitalized with COVID-19, peripheral blood was obtained from patients admitted to the medical (non-ICU) wards at a large tertiary hospital. We demonstrated a decrease in the mature cytotoxic subset of NK cells within the peripheral blood of patients hospitalized with COVID-19. We also observed a notable reduction in the cytotoxic function of NK cells against tumor targets. We examined the mechanisms leading to NK cell killing. We found reductions in the intracellular levels of effector molecules, the degranulation of cytotoxic granules, and the extracellular concentrations of released effector molecules. We identified dysfunctional intracellular granule trafficking required to position the granules for degranulation, which would be consistent with the reduced release of effector molecules. We found clusters of inhibitory receptors were upregulated in subsets of NK cells, in keeping with inhibition of cytotoxicity. Additionally, males with COVID-19 showed NK cell defects compared to healthy males, while no significant differences were observed in females. Our findings highlight defects in cytolytic effector molecules, granule trafficking and release, and increased expression of inhibitory receptors on NK cells in patients hospitalized with COVID-19, in addition to a sex difference in cytolytic function, which contributes to defective NK cell function in COVID-19.
Introduction
Importance of COVID-19
The COVID-19 pandemic began in December 2019 with the emergence of a novel member of the coronavirus family, SARS-CoV-2. As of July 14, 2024, more than 778 million people have been infected globally, resulting in more than seven million deaths worldwide (World Health Organization, 2024). The infection spans the spectrum from asymptomatic to severe respiratory disease, leading to hospitalization, admission to the ICU and death (Ziehr et al., 2020; ; ). The devastating impact of COVID-19 has led to intense interest in the host’s defense against the infection, involving organ-specific defense and innate and adaptive immunity, including natural killer (NK) cells.
The role of NK cells in host defense and control of activation
Natural killer (NK) cells are innate lymphoid cells that serve as one of the first lines of defense against tumor cells and microbial pathogens, including bacteria, fungi, parasites, and viruses (; ; ; ; ). Unlike other cytotoxic lymphocytes, natural killer cells do not require previous exposure to recognize and deploy their cytotoxic products to attack and kill pathogens, pathogen-infected or malignantly transformed cells (). Viruses and tumors can cause a reduction in inhibitory MHC-I molecules on the target cells, making them susceptible to NK cell targeting through activating receptors ().
Effect of viruses on NK cell activity
Some viruses stimulate NK cell activity, while others suppress their function, leading to contrasting effects on host defense. For instance, the human immunodeficiency virus and respiratory syncytial virus have been shown to enhance NK cell activity (; van Erp et al., 2019). By contrast, other viruses, such as influenza A, human pegivirus, and vaccinia, reduce NK cell function in infected individuals (; ; ). Different viruses can manipulate various pathways of activation and inhibition in NK cells. For example, human cytomegalovirus and human T lymphotropic virus type 1 enhance NK cell function through augmented degranulation and increased IFN-γ production. However, they can also inhibit NK cell function by altering receptor expression, including checkpoint proteins (; Tschan-Plessl et al., 2016; ). Consequently, it is necessary to understand the sequence of events during NK cell function required to predict the effect of viruses under different conditions.
Pathways leading to cytotoxicity
When an NK cell encounters a tumor or virus-infected cell, it kills the target cell by releasing cytotoxic effector molecules. These molecules, such as perforin, granulysin, and granzymes, are released from the NK cells’ granules, also called vesicles or secretory lysosomes. This degranulation event requires the receptor activation signal to be stronger than the inhibitory signal. This triggers cytoskeletal rearrangement, which is necessary for the proper trafficking of vesicles within the NK cells toward the immunological synapse, where they can release their cytotoxic proteins.
Perforin damages or creates pores in the target cell membrane, allowing the entry of granzymes (; ). Additionally, granulysin associates with the plasma membrane of tumor cells, causing mitochondrial damage and calcium overload, leading to cell death (). Granzymes, including Granzyme A, B and K, are a family of enzymes that induce caspase-dependent cell death (; ; ; Zhao et al., 2007). While studies have described the intracellular levels of granzymes A, B and perforin in patients with COVID-19, it still needs to be clarified whether the extracellular concentration of these molecules is altered upon stimulation.
NK cell activity in COVID-19
NK cells are dysfunctional in severe COVID-19. In hospitalized patients with severe COVID-19 (grade 6 to 9 on the WHO progression scale), the number of NK cells is reduced, and the expression of IFN-stimulated genes and genes involved in IFN-α signaling is increased, leading to defective NK cell cytotoxicity (; Zheng et al., 2020; ). IFN-α stimulates NK cell function and can downregulate IFN-γ production (; ). When considering intracellular effector molecule content, there is less consensus about whether granzyme B and perforin expression increase or decrease (; Zheng et al., 2020; Witkowski et al., 2021). There is also an impairment in cytotoxic degranulation assessed by CD107a presentation (; Witkowski et al., 2021).
Independent studies have looked at limited panels of activating and inhibitory receptors, and the conclusions drawn so far show an exhausted phenotype based on the expression of NKG2D, NKG2A, PD-1, CTLA-4, and TIM-3 (Zheng et al., 2020; Varchetta et al., 2021; Ustiuzhanina et al., 2023). Transcriptional analysis has proposed that these defects are due to a transforming growth factor-β (TGF-β) signature leading to reduced cell-cell adhesion, granule exocytosis and cytotoxicity (Witkowski et al., 2021). However, a systematic assessment of the sequence of intracellular events leading to NK cell cytotoxicity and the expression of inhibitory receptors as possible mechanisms of defective cytotoxicity in patients hospitalized with COVID-19 (PHWC-19) has not been performed.
Demographic and clinical factors affecting NK cell responses in COVID-19
Demographic and clinical factors can affect NK cell intracellular mechanisms. Sex differences are present across many aspects of the immune system and are essential to understand and treat infections and immunologic disorders fully. Males are more likely to become hospitalized when infected with COVID-19 and comprise the majority of all age and ethnic groups in hospitalized patients (). There are differences in receptor and cytokine expression, which could be the basis for sex differences in the immune response to SARS-CoV-2. Despite the relationship between sex differences and the severity of COVID-19, there is a lack of information about NK cell function in males and females hospitalized with COVID-19.
The ancestral strain of SARS-CoV-2 predates the subsequent variants of concern (VOCs) that evolved later in the pandemic. These variants, including Alpha, Beta, Gamma, Delta, and Omicron, are progressively more transmissible due to genetic adaptation and external pressures such as vaccine-induced and natural immunity (; ). It has been shown that the CD56dim CD16bright NK cells from patients infected with Omicron BA.2 demonstrated greater IFN-γ secretion and cytotoxicity in previously vaccinated individuals compared to unvaccinated controls (). However, there was no comparison of NK cell functioning in patients infected with earlier VOCs.
Dexamethasone, a corticosteroid, was introduced as an anti-inflammatory medication to reduce the progression of organ injury and failure that could result in death (). Moderate administration of dexamethasone benefits PHWC-19 with respiratory support (supplemental oxygen or ventilation), and no advantages have been found in patients without respiratory aids (). There have been studies regarding the effects of corticosteroids, with some studies showing an enhancement of proliferation and cytokine production of NK cells while others showing reduced cytotoxicity (; ; ). The role of dexamethasone as an independent variable on NK cell function during COVID-19 infection has not been investigated. One study found no difference in the cytotoxicity of NK cells from patients before and after the establishment of dexamethasone as a standard treatment for COVID-19, but did not investigate further (Witkowski et al., 2021).
The significance and approach
We utilized antibody labelling and mass cytometry to determine the percentages of different blood NK cell subsets and other lymphocytes reduced in PHWC-19. NK cells from PHWC-19 were co-cultured with a tumor cell target to compare cytotoxicity with NK cells from healthy subjects. This led to our investigation of the mechanisms involved in NK cell cytotoxicity. Using flow and mass cytometry, we found that NK cells from PHWC-19 had reduced effector molecule expression, degranulation, and release using flow cytometry and ELISA. We used live cell imaging to establish that granule trafficking and polarization were defective in NK cells from PHWC-19, which is critical for granule release and effectiveness. We used a panel of antibodies of activating and inhibitory receptors to show that inhibitory receptor expression was enhanced, which could affect downstream granule trafficking, granule content, and degranulation. Lastly, we explored the role of demographic and clinical factors like age, sex, VOCs, prior vaccination and dexamethasone treatment on NK cell functioning to understand how NK cells are affected by COVID-19.
Materials and methods
Study participants
Peripheral blood samples were collected from forty-four patients hospitalized in a medical ward at Foothills Medical Centre, Calgary, Canada, after informed consent was obtained between November 2020 and September 2023. The patients were categorized within the moderate range of 4 or 5 on the WHO clinical progression scale for COVID-19 disease severity (). These patients presented with clinical features requiring hospitalization but were not admitted to the intensive care unit. Patients were excluded if they had major comorbidities that might affect study outcomes, were immunosuppressed, including having an active malignancy or receiving immunosuppressive drugs. Characteristics such as age, sex, vaccination status, dexamethasone treatment, and variants at the time of hospitalization are summarized in Table 1. When patient samples were obtained, peripheral blood samples were collected on the same day from forty-five healthy donors at the Cumming School of Medicine to act as the paired negative controls.
Table 1
| Characteristic1 | Non-ICU hospitalized patients (N = 44) | Healthy controls (N = 45) |
|---|---|---|
| Age group – No. | ||
| Mean ± SD | 63 ± 19 | 36 ± 19 |
| 20-40: | 6 | 32 |
| 41-60: | 13 | 5 |
| 61-80+: | 25 | 8 |
| Male Sex – No. (%) | 24 (54%) | 21 (47%) |
| Dexamethasone Treatment – No. | ||
| Yes: | 22 | N/A |
| No: | 22 | N/A |
| Vaccination Status – No. | ||
| Yes: | 19 | 34 |
| No: | 21 | 7 |
| Unknown: | 4 | 4 |
| Presumed and Confirmed Variants – No. | ||
| Ancestral: | 9 | N/A |
| Alpha: | 7 | N/A |
| Delta: | 15 | N/A |
| Omicron: | 13 | N/A |
Characteristics of study subjects.
Participant demographics in the healthy and non-ICU hospitalized COVID-19 cohort. Information includes age, sex, dexamethasone treatment, vaccination status, and confirmed and presumed variant status based on the timing of diseases.
Primary natural killer cell isolation from whole blood
Peripheral blood mononuclear cells (PBMCs) were isolated from whole blood using the Ficoll-Paque® separation technique (Millipore Sigma, Germany). PBMCs were used for flow cytometry studies, while primary NK cells (pNK) were isolated by negative selection using magnetic beads (Miltenyi Biotec, Germany). Three mL of whole blood was put aside for mass cytometry.
NK cellular cytotoxicity
NK cell cytotoxicity was assessed using 721.221 tumor cells (a gift from Dr. Nicole Bernard, McGill University, Montreal). The 721.221 cells were cultured in RPMI 1640 (Gibco, Thermo Fisher Scientific Inc, Waltham, MA) with 10% fetal bovine serum (FBS) and 1% Penicillin-Streptomycin-Glutamine (PSQ). The 721.221 cells were stained with the carboxyfluorescein succinimidyl ester (CFSE) (Millipore Sigma, Germany) as a cell tracking dye to distinguish them from the primary cells. The cells were co-cultured with NK cells in RPMI 1640 with 10% FBS, 1% non-essential amino acids, 1% Na pyruvate, 1% PSQ at 5:1 and 10:1 effector-to-target (E:T) ratios for 4 hours before adding 7-amino actinomycin D (7AAD) (Invitrogen, Thermo Fisher Scientific Inc, Waltham) and analyzed by flow cytometry (Guava EasyCyte, EMD Millipore, Danvers, MA). The positive control for cell death was obtained by fixing 721.221 cells with 4% paraformaldehyde (PFA) followed by Triton-X100 (0.1% in PBS). The percentage of cell death is analyzed using the FlowJo LLC program by examining the percentage of cells that are double positive for the CFSE and 7AAD stains.
Antibody labelling and flow cytometry for cell surface markers, intracellular cytotoxic effector molecules, and checkpoint inhibitor receptors
PBMCs were fixed with 4% PFA and permeabilized with saponin (0.1%). Intracellular effector proteins granulysin, perforin, and granzymes A, B, and K, and checkpoint proteins LAG-3, PD-1, TIGIT, and CD94 were labelled with fluorescent antibodies or their respective isotype control antibodies for 30 minutes along with anti-CD56 and anti-CD3 to define NK cells (Supplementary Table 1). Geometric mean fluorescence intensity was assessed by flow cytometry (LSRII, BD Biosciences, New Jersey). The healthy and COVID-19 samples were collected and processed simultaneously to limit interexperimental variability. Initially, PBMCs were gated on CD56+ CD3- cell population for NK cells. Histograms were then utilized to determine the percentage of NK cells positive for each marker. Flow cytometry data was analyzed using the FlowJo LLC program.
NK cell degranulation
NK cell degranulation was assessed as previously described (). The Anti-CD107a fluorescent antibody and the intracellular protein transport inhibitor cocktail Brefeldin A/Monensin were added to the NK cell plus 721.221 cell co-culture for four hours. Cells were fixed with 4% PFA, labelled with CD56 and CD3, and analyzed by flow cytometry (Guava EasyCyte, EMD Millipore, Danvers, MA). These experiments were conducted concurrently in the healthy and COVID-19 samples to pair the results.
Effector molecule release ELISA assay and analyses
Supernatants from the co-cultures of NK and tumor cells were collected at 4 hours and stored at -80°C prior to measuring levels of released perforin, granzyme A and granzyme B by ELISA. Effector molecule release was measured using the manufacturer’s reagents and instructions (Abcam, Cambridge).
Granule trafficking
Live-cell imaging of primary NK cells in the presence of 721.221 target cells was used to assess spatial and temporal changes in conjugate formation and granule trafficking. The NK cells were treated with Cresyl violet to label granules (Millipore Sigma, Germany) and tubulin tracker green to label microtubules (Invitrogen, Thermo Fisher Scientific Inc, Waltham, MA). A Ti Eclipse microscope (Nikon Instruments Inc., Melville, NY) was used with an objective lens with a magnification of 60X, a numerical aperture of 1.2 and oil immersion. Images were analyzed using the ImageJ Fiji program following the protocol described previously to measure granule distances to the microtubule-organizing center (MTOC) and the synapse (). The synapse is defined as the center of the area of contact between the NK cell and the 721.221 tumor cell. The width of the synapse was determined when the plasma membrane of the NK cell and the tumor cell was deformed by the contact and defined as the linear distance between the points of contact of both cells’ plasma membranes. The granule distance is defined as the mean distance from the granules to the closest point on the immunological synapse. Distances of the MTOC to synapse and granule to synapse were determined over time using the Trackmate plugin in ImageJ (). Imaging in the healthy donor was performed at a frame rate of 1 image every minute for 60 minutes, and imaging in the PHWC-19 was performed at a frame rate of 1 image every minute for 80 minutes. The time to initial polarization was defined as the time between the start of a sustained contact between the NK cell and the tumor cell and the initial decrease in distance of the granules to the synapse.
Surface receptors assessed by mass cytometry (CyTOF)
The whole blood was diluted at a ratio of 1:1.4 of blood to the PROT1 stabilizer for ten minutes and stored at -80°C until the samples were prepared for mass cytometry (HELIOS, Standard BioTools Inc, San Francisco). The cells were thawed, and the red blood cells were lysed using the PROT1 lysis buffer and washed with a cell staining media (CSM) consisting of phosphate-buffered saline (PBS) (Gibco, Thermo Fisher Scientific Inc, Waltham, MA) with 1% bovine serum albumin (BSA). The samples were fixed, permeabilized and barcoded with the Cell-ID 20-Plex Pd Barcoding Kit (Standard BioTools Inc, San Francisco) according to the manufacturer’s recommendations. After barcoding, the cells were adjusted to have 5 x 105 cells per sample for 3 x 106 cells per tube. Cells were washed, and samples were combined and stained with a surface staining cocktail of lanthanide-labelled antibodies (Supplementary Table 2) for 30 mins at room temperature. Cells were reconstituted in 0.3% Saponin/PBS with iridium and 1.6% paraformaldehyde in the fridge overnight. The next day, the cells were washed with CSM and PBS before being aspirated into a pellet and brought to the HELIOS machine for acquisition. All the data was uploaded to the Cytobank program to generate all gates and run all analysis methods, including viSNE for dimensionality reduction and FlowSom for clustering cells. GraphPad Prism (GraphPad Software Inc, CA) was used to conduct the PCA analysis of the data from mass cytometry. Gating strategies for cell populations have been provided in Supplementary Figure 1.
Statistical analysis
Statistical analysis was performed using GraphPad Prism (GraphPad Software Inc, CA). A paired t-test was used to assess cytotoxicity and degranulation after ensuring normal distribution based on a Shapiro-Wilk test. The remaining analysis was done with Mann-Whitney, Wilcoxon matched-pairs signed rank, or Kruskal-Wallis ANOVA based on paired or unpaired samples and the number of comparisons. P<0.05 was considered statistically significant.
Results
PHWC-19 have reduced numbers of peripheral blood CD56+ CD3- NK cells
NK cells play an essential role in the innate immune system during many respiratory diseases. These cells can kill tumors, virus-infected cells, and bacterial and fungal pathogens (; ; ). Since the number of NK cells could affect the severity and degree of immune dysfunction in PHWC-19, the peripheral whole blood of PHWC-19 was analyzed using mass cytometry. We found a significant decrease in the percentage of NK cells and, indeed, all mononuclear cells in PHWC-19 compared to healthy donors (Figure 1A, gating strategy in Supplementary Figure 1). Using flow cytometry, we found the percentage of CD3- CD56+ NK cells in PHWC-19 was lower compared to healthy controls processed on the same day (9.8 ± 3.8% in PHWC-19, 22.8 ± 8.0% in healthy subjects) (Figure 1B). We further assessed different NK cell subsets based on the expression of CD56 and CD57. We found that the mature, differentiated, and cytotoxic NK cell subset, CD56dim CD57+, was decreased in PHWC-19, while the immature cytotoxic NK cells, CD56dim CD57- increased (Figure 1C) (; Vivier et al., 2008). The less abundant cytokine-producing CD56bright CD57+/- population remained unaltered due to COVID-19 (Figure 1C) (; ). Although the reduced numbers and maturity of NK cells may contribute to the observed reduction in cytotoxicity, we were interested in determining whether there is a defect at the cellular level of NK cell function.
Figure 1
PHWC-19 has decreased cytotoxic abilities in peripheral CD56+ CD3- NK cells
To test the cytotoxic ability of NK cells from PHWC-19, we used an established method of tumor target killing (; ; ). NK cells were co-cultured with the 721.221 cell line, and killing was assessed using 7-AAD uptake and flow cytometry. There was a decrease in 721.221 cell death in response to NK cells from PHWC-19 at E:T ratios of 5:1 and 10:1 compared to healthy control subjects (Figure 2A).
Figure 2
Since the mean age of PHWC-19 was older than that of our healthy control subjects, and aging has been reported to affect some aspects of immune cell function, we tested the effect of age on NK cell cytotoxicity (; ). Linear regression analysis of the healthy donors showed no decline in the percentage of 721.221 cells killed when co-cultured with NK cells from older healthy subjects (Figure 2B). Surprisingly, there was a significant increase in the percentage of 721.221 cells killed with increasing age of the PHWC-19 (Figure 2B). Although there were fewer young PHWC-19, results were significant, revealing that COVID-19 caused a defect in cytotoxicity that was greater in younger patients than in older patients who required hospitalization.
NK cell effector molecules and degranulation were altered in PHWC-19
After establishing that cytotoxicity was reduced in NK cells from PHWC-19, we were interested in the sequence of intracellular events that might account for the reduced cytotoxicity. Effector molecules are preformed and stored in granules before their release. Granzymes A, B, and K are serine proteases, while perforin and granulysin are membrane-disrupting or pore-forming proteins involved in the defense against tumors and viruses (; van Domselaar and Bovenschen, 2011; ). Using intracellular antibody labelling and flow cytometry, we compared the relative abundance of intracellular granzymes A, B, K, perforin and granulysin in NK cells from PHWC-19 and healthy adults. We found that only the intracellular levels of granzyme B were significantly decreased in the NK cells of PHWC-19 compared to those from healthy subjects (Figure 3A), and the fluorescent intensity of granzyme B correlated with cytotoxicity (Supplementary Figure 2A), revealing that the level of granzyme B was important for cytotoxicity. Intracellular granzyme A showed a lower trend in NK cells from PHWC-19 but did not reach statistical significance using a Wilcoxon matched-pairs signed rank test (Figure 3A).
Figure 3
In contrast, granzyme K showed no significant differences between the healthy and COVID-19 groups (Figure 3A). Among the pore-forming effector molecules, intracellular perforin and granulysin were not significantly different between the healthy and COVID-19 samples (Figure 3A). The intracellular content of granzymes and pore-forming proteins reflects the stores available. However, since only a fraction of the stores are utilized during each round of killing, the effect would potentially manifest only after multiple rounds of killing, leading to the depletion of the stores (). Because of this, we were interested in whether degranulation was affected.
Degranulation occurs when the effector-containing granules fuse with the NK cell plasma membrane and release their contents. To assess degranulation in NK cells from patients hospitalized with COVID-19, flow cytometry was used to detect the cell surface membrane expression of CD107a (LAMP1) during the killing of tumor cells. There was a significantly lower expression of CD107a after co-culture of NK cells from PHWC-19 with 721.221 cells, and the fluorescent intensity of CD107a correlated with cytotoxicity (Supplementary Figure 2B), revealing the importance of reduced degranulation of NK cells from hospitalized SARS-CoV-2 infected patients compared to healthy NK cells (Figure 3B).
The release of effector proteins during degranulation increases the concentration of effectors in the extracellular compartment. Since degranulation was reduced in NK cells from patients hospitalized with COVID-19, we then determined if there was a difference in the amount of effector molecules released. The supernatants of the co-cultures were collected, and the effector molecules, granzyme A, B, and perforin, were measured by ELISA. There was a significant decrease in the concentration of granzymes A and B released from the NK cells obtained from PHWC-19 (Figure 3C). Meanwhile, extracellular concentrations of perforin were not significantly different between the healthy and COVID-19 cohorts (Figure 3C). This data reveals that intracellular stores of critical effector molecules were reduced, and degranulation and the release of granzymes were impaired in NK cells from PHWC-19, which was associated with reduced cytotoxicity.
Granule trafficking was impaired in NK cells from PHWC-19
To release the effector molecules at the target cell interface, the intracellular cytolytic granules must move along the microtubules and reach the area of contact with the target cell. We sought to determine if there was impaired trafficking of the granules and the microtubule-organizing center (MTOC). Granule trafficking requires the granules to converge on the MTOC via microtubule motor proteins and for the converged granules and MTOC to polarize to the NK cell synapse with the target cell (). Since enveloped viruses disrupt many forms of cellular transport (), we asked whether the defect in granule release was due to deficiencies in trafficking or failure to reach the plasma membrane. To examine granule movement, we used live-cell imaging of NK cells from healthy subjects and PHWC-19 in the presence of 721.221 target cells. NK cells from healthy subjects rapidly mobilized the granules and MTOC to be adjacent to the tumor cell, indicating successful granule trafficking (Figure 4A, Supplementary Video 1A). By contrast, the cytolytic granules and MTOC in NK cells from PHWC-19 did not polarize to the synapse with the target cell, indicating failed granule trafficking, which would lead to defective granule release and cytotoxicity (Figure 4A, Supplementary Video 1B). The imaging of NK cells from hospitalized patients with COVID-19 was continued for an additional 20 minutes to capture any events that may have taken longer to occur. Quantifying the mean granule-to-synapse distance indicated decreased localization of granules towards the synapse in NK cells from PHWC-19 (Figure 4B), although the time to initial polarization was similar in healthy and PHWC-19 (Supplementary Figure 3). Similarly, reduced localization of the MTOC was also seen in the NK cells of PHWC-19 (Figure 4C).
Figure 4
NK cells from healthy subjects formed a broad, concave base of attachment between the NK cell and the tumor cell. By contrast, the NK cells from PHWC-19 formed a synapse with the tumor cells, but the synapse was not as intimate as the synapse formed by NK cells from healthy control subjects, and the width of the synapse was greater in NK cells from healthy control subjects than in NK cells from PHWC-19 (Figures 4A, D, Supplementary Video 1B). Thus, NK cells from PHWC-19 demonstrated impaired granule trafficking and a less intimate synapse formation that would contribute to defective cytotoxicity.
Expression of inhibitory and activating NK cell surface receptors
Because of the suboptimal conjugate formation and granule polarization in NK cells from PHWC-19, we hypothesized that there might be changes in inhibitory receptor expression, which could affect granule trafficking and degranulation. The inhibitory receptors include a subset called checkpoint inhibitory proteins. These proteins regulate the NK cell function by reducing their cytolytic activity. We performed experiments using samples of NK cells obtained from healthy and PHWC-19 on the same day to assess PD-1, LAG-3, TIGIT, and CD94. Only PD-1 was significantly increased in NK cells from PHWC-19 (Figure 5A). Meanwhile, LAG-3 expression was reduced in PHWC-19, and TIGIT and CD94 showed no difference between the healthy and COVID-19 groups (Figure 5A). Since different cell subsets are present, we examined the expression patterns of inhibitory and activating receptors on different NK cell subsets defined by CD56 and CD57 expression. Using mass cytometry, we studied the expression of the inhibitory and checkpoint receptors PD-1, LAG-3, TIGIT, CTLA4, TIM-3, Siglec 9, NKG2A, NKB1, CD161, CD158b, CD47 and the activating receptor NKG2D on these subsets. The differential expression of each receptor within the subsets was visualized via a log2 heatmap (Figure 5B). The more mature and most cytotoxic CD56dim CD57+ NK cells in PHWC-19 showed a higher expression of PD-1, CTLA-4, and CD158b and a decrease in CD47 compared to healthy controls, suggesting that the net effect is to dampen cytotoxicity (Figure 5C). Similarly, the more abundant and slightly less cytotoxic CD56dim CD57- NK cells showed a significant upregulation of PD-1, LAG3, CTLA4, and CD158b, along with a downregulation of NKG2D, NKG2A, and CD161 in PHWC-19 compared to healthy controls. This also suggests a net negative effect (Figure 5C). Overall, these cytotoxic NK cell populations exhibit an increased expression of inhibitory receptors in PHWC-19, which is consistent with reduced killing. In the less prevalent CD56bright CD57- NK cells, which are associated with cytokine production and lower cytotoxicity, there was an increased expression of PD-1, CTLA4, CD158b, and decreased expression of NKG2D, CD161, and CD47 in PHWC-19 (Figure 5C), revealing a receptor profile consistent with a negative effect on signaling.
Figure 5
We identified major clusters within CD56dim CD57- and CD56dim CD57+ NK cells based on the co-expression of activating and inhibitory receptors (Figure 6A). This aided us in identifying alterations in NK cells with specific receptor co-expression between healthy donors and PHWC-19. The major clusters for CD56dim CD57- NK cells were 1, 4, 7, 9, and 10 (Figure 6B). Among these, cluster 4 was significantly increased in PHWC-19, while cluster 7 was significantly decreased in PHWC-19 (Figure 6B). Cluster 4 of CD56dim CD57- NK cells was characterized by higher co-expression of CTLA4, CD161, and CD47, while cluster 7 showed higher co-expression of NKG2D, CTLA4, and CD161 (Figure 6C). Dendrogram clustering revealed three major groups of receptors co-expressed on CD56dim CD57- NK cells, which showed similarities between healthy and PHWC-19 groups (Figure 6C).
Figure 6
The major clusters for CD56dim CD57+ NK cells were clusters 1, 9, and 10 (Figure 6B). Cluster 9 was significantly decreased in PHWC-19 (Figure 6B). This cluster showed higher co-expression of NKG2D, TIGIT, CTLA4, CD161, and CD47 (Figure 6C). Dendrogram clustering revealed four major groups of receptors co-expressed on CD56dim CD57+ NK cells, which were similar between healthy and PHWC-19 donors (Figure 6C). We also investigated the clusters within CD56bright CD57- and CD56bright CD57+ NK cells, but these populations did not show major differences between healthy donors and PHWC-19 (Supplementary Figures 4A–C). The dendrograms of CD56bright CD57- and CD56bright CD57+ NK cells showed similar groupings in the healthy and PHWC-19 groups (Supplementary Figure 4C). Additional PCA analysis was performed on CD56bright and dim and CD57+ and – subsets (Supplementary Figure 5). These results reveal that specific combinations of inhibitory receptors cluster into major subsets and that these clusters differ in their expression in PHWC-19 compared to control subjects.
Effect of variants, dexamethasone, sex, and vaccine status on NK cell killing in patients hospitalized with COVID-19
After establishing the intracellular mechanistic defects responsible for impaired cytotoxicity in PHWC-19, we wondered how other biological, therapeutic and virologic factors might influence the defect. These characteristics include the different variants of the SARS-CoV-2 virus, the treatment of patients with dexamethasone, and the sex-related differences. The SARS-CoV-2 virus underwent genetic mutations, producing variants of concern. The severity of the disease and immune responses differed depending on the variant causing the infection (; ). The major variants of concern sequentially evolved from Alpha to Beta, Gamma, Delta, and Omicron. These variants have shown increased transmissibility, decreased virulence, and antigenicity (defined by antibody binding) (). We sought to determine the cytotoxicity of NK cells obtained from PHWC-19 infected with different SARS-CoV-2 variants compared to healthy control subjects. Unfortunately, the sample size was small, and while the differences in the degree of impaired NK cell killing did not reach statistical significance, it is worth noting that there is a trend suggesting that the reduction in NK cell cytotoxicity was greater with the Alpha variant compared to NK cells from PHWC-19 infected with Delta and Omicron (Figure 7A), which correlates with the lack of effective host defense and worsened pathology in patients infected with earlier variants of concern.
Figure 7
During the transition from the time when the Alpha variant was dominant to the time when the Delta variant was dominant, dexamethasone was routinely added to therapy for the management of hypoxemic hospitalized patients (). It is possible that dexamethasone diminished the function and the number of NK cells (). To address this possibility, we performed experiments to determine whether dexamethasone affected cytotoxicity in PHWC-19. In these experiments, multiple comparisons were made. Control PHWC-19 who were enrolled before the use of dexamethasone or where the care team chose not to use it served as the comparator for PHWC-19 that received dexamethasone, and NK cells from healthy subjects were compared to their NK cells treated in vitro with dexamethasone. Previous studies suggested that the serum concentration after a dose of 6mg would be 0.1μM (Weijtens et al., 1998; ; ). To evaluate the effect of the equivalent of 6mg dexamethasone on NK cells in vivo, NK cells from healthy control subjects were pretreated with 0.1 μM in vitro. There was no significant reduction in NK cell cytotoxicity between the NK cells from PHWC-19 treated or not treated with dexamethasone, and no reduction in cytotoxicity among healthy subjects treated with dexamethasone compared to the control in vitro (Figure 7B).
Differences in biological sex influence many aspects of the immune system, and we are beginning to learn their effect on infections. Males were more likely to become hospitalized when infected with SARS-CoV-2, and males comprise the majority of all age and ethnic groups in hospitalized patients (). To determine whether biological sex affected NK cell cytotoxicity, the responses from NK cells from male PHWC-19 were compared to those of healthy male control subjects, and responses of NK cells from female PHWC-19 were compared to those of healthy female control subjects. There was a significant decrease in the killing of NK cells in male PHWC-19 compared to the healthy male subjects (Figure 7C). By contrast, we did not detect a difference in NK cell cytotoxicity in female PHWC-19 compared to healthy control females.
A multivariate analysis was performed to examine the influence of age, healthy vs PHWC-19, the influence of dexamethasone, and sex. The regression formulation for an E:T ratio of 10:1 was Cytotoxicity = 31.898 + 0.340*Age + 22.722*Status (healthy) + 9.143*Dexamethasone - 14.046*Sex (males). Where P = 0.010 (Age), 0.001 (Status-healthy), 0.015 (Sex-Male).
To determine whether vaccine status affected NK cell function in patients who had been hospitalized with COVID-19, the cytotoxicity was compared between NK cells from PHWC-19 who had been vaccinated and those who had not. This demonstrated no difference in the cytotoxicity (Supplementary Figure 6), suggesting that once hospitalized with COVID-19, vaccine status did not influence the reduction in NK cell function.
Assessment of sex differences in cytotoxic mechanisms in healthy and hospitalized patients positive for COVID-19
Having shown that the defects in cytotoxicity resulting from COVID-19 disease were correlated with the male but not female sex, we were interested in exploring whether the intracellular defects also exhibited a similar pattern of sex-related differences. Unfortunately, the sample size was small, and while trends emerged, the intracellular expression of effector molecules, degranulation or extracellular free effector molecules released into the supernatant failed to demonstrate statistically significant differences between males and females (Supplementary Figures 7A–C).
We also assessed the expression of activating and inhibitory receptors in male versus female participants to determine whether there were differences. NK cells from male PHWC-19 showed an increased expression of PD-1 compared to healthy subjects, while females did not (Figure 8A). By contrast, there were no differences in TIGIT and CD94 expression in males or females (Figure 8A).
Figure 8
Further examination of the NK cell subsets revealed similar changes in males and females, including upregulation of PD-1, CTLA4, and CD158b, compared to sex-matched controls for the CD56dim CD57- NK cells (Figure 8B). Meanwhile, CD56dim CD57+ NK cells from male PHWC-19 showed a higher expression of CTLA4 and CD158b and NK cells from female PHWC-19 had a higher expression of CTLA4, which was balanced by a lowered CD47 expression (Figure 8B). Lastly, in the CD56bright CD57- NK cells, male PHWC-19 exhibited more of an inhibitory phenotype with increased CD158b and decreased NKG2D expression (Figure 8B). In females, more inhibitory receptors are expressed in the healthy donors or evenly expressed between the healthy and PHWC-19 cohorts (Figure 8B). In summary, there were significant differences in the male patients, including reduced activating receptor expression and increased inhibitory receptors.
Two studies have shown the sexually dimorphic nature of dexamethasone treatment that is more beneficial in men than in women (; ). We were interested in whether this effect extended beyond neutrophils and involved NK cells. While there was no difference at a 5:1 ratio, at the 10:1 ratio, males who did not receive dexamethasone had a trend towards a more pronounced NK cell cytotoxicity defect than male PHWC-19 who received dexamethasone (Supplementary Figure 8). While dexamethasone treatment may not show overall effects on NK cell function, dexamethasone may have a role in recovering NK cell cytotoxicity specific to males.
Discussion
NK cells are crucial in defending against viral infections, including killing SARS-CoV-2-infected cells. While other viruses affect NK cell responses, gaps remain in understanding NK cell function in SARS-CoV-2-infected patients. Most studies enrolled severe hospitalized COVID-19 cases; our study focused on patients with moderate disease not requiring ICU care. Few studies have explored NK cell function across different COVID-19 variants, vaccination statuses, and dexamethasone use. Studies have not previously investigated the sequence of events leading to NK cell cytotoxicity to completely understand how the infection changes NK cell responses over time. Our study is novel in examining the interactions between receptor subset clustering, stores of effector molecules, granule trafficking, polarization of granules, and effector molecule release in the context of a viral infection, and the roles of sex and dexamethasone on NK cell cytotoxicity and mechanisms. Our research will drive NK cell-dependent therapies for PHWC-19 by targeting the possible dysfunctions we outlined ().
Our study builds upon previous work in the field, showing a reduction in the percentages of all mononuclear cells and, more specifically, a reduction in the proportion of NK cells from patients hospitalized with COVID-19 (Zheng et al., 2020; ; ). We also found that CD56dim NK cells comprise most NK cells, and the ratio of CD56dim to CD56bright was preserved in PHWC-19 (). However, our studies revealed that CD56dim CD57- NK cells were increased in PHWC-19, while CD56dim CD57+ cells were reduced. This change contrasts with HIV, in which CD57- cells are reduced, while the number of CD57+ cells remains the same (). In addition to reducing the number of NK cells, it indicates a shift to a more immature, less cytotoxic phenotype in PHWC-19 that could contribute to the observed defect in NK cell cytotoxicity ().
While the number of NK cells is important, their cytotoxic function is crucial to maintaining host defense. Previous studies showed reduced NK cell degranulation in severe COVID-19 patients against K562 and SW480 tumor cells but did not assess cytotoxicity or moderately ill patients (; Varchetta et al., 2021). Reduced CD56dim CD57+ NK cells may impair cytotoxicity, as they are more cytotoxic than increased immature CD56dim CD57- NK cells ().
We examined age-related changes, noting that although the average age of the COVID-19 group was older, younger PHWC-19 showed a relatively greater defect in NK cell cytotoxicity. Increasing age did not adversely affect NK cell cytotoxicity in either cohort. Surprisingly, advanced age was associated with improved NK cell function in PHWC-19, consistent with studies linking NK cell function more to health status and comorbidities than age (; ; ). The increase in NK cell function was even more apparent in healthy control subjects, making it unlikely that inflammation or NK cell memory contributed to the increase. We further investigated the mechanisms behind the impaired NK cell cytotoxicity.
We began investigating NK cell processes sequentially by examining the effector molecules and granule-mediated release. NK cells deploy stored lytic granules to release the effector molecules during cytotoxic responses. Granzymes and pore-forming effector molecules are stored within the NK cell and, when needed, are released upon activation. Perforin is crucial in helping granzymes enter target cells (). Without perforin, granzymes are less effective at entering target cells, directly affecting cytotoxicity. While some studies indicate elevated perforin in severe COVID-19 patients, others have reported reduced or unchanged levels (; ; ; ; Witkowski et al., 2021; Zenarruzabeitia et al., 2021). Our findings align with the latter, showing no difference in perforin levels between healthy and COVID-19 groups, possibly due to differences in disease severity (Zenarruzabeitia et al., 2021).
We observed an overall reduction in granzyme B in NK cells taken from moderate illness patients, which reflects findings in patients with severe illness from COVID-19 (; Zheng et al., 2020). We extended these findings by showing lowered degranulation and extracellular concentrations of granzyme B in the supernatant post-stimulation. By contrast, other reports observed an early upregulation of granzyme B, primarily evident in patients with severe COVID-19 (; Witkowski et al., 2021; Zenarruzabeitia et al., 2021). Though less studied, granzyme A is critical for NK cell cytotoxic function (; ; van Domselaar and Bovenschen, 2011). Previous studies noted increased granzyme A+ NK cells in PHWC-19 but did not measure extracellular levels (; ; ). Our study revealed a significant reduction in the extracellular concentrations of granzyme B and a trend toward a decrease in granzyme A from NK cells of PHWC-19. This discrepancy may depend on disease severity. We also demonstrated a positive correlation between intracellular granzyme B levels and cytotoxicity, which supports the idea that a reduction in granzyme B levels in PHWC-19 would contribute to decreased cytotoxic abilities of NK cells.
Intracellular effector molecule levels reflect cellular readiness, but assessing their release is crucial for understanding impaired cytotoxicity. We observed decreased degranulation of NK cells from PHWC-19 in response to 721.221 cells, and found a positive correlation between degranulation and the cytotoxic ability of NK cells, extending previous reports of reduced CD107a expression in response to K562 tumor cells in PHWC-19 (Zheng et al., 2020; Witkowski et al., 2021). Other viruses, such as influenza A virus, exhibited inhibition of granule release against tumor cells, while human T-lymphotropic virus-1 and human cytomegalovirus upregulate degranulation (; Tschan-Plessl et al., 2016; ). Extracellular levels of perforin and granzymes A and B revealed decreased granzyme release from PHWC-19 NK cells, contributing to reduced cytotoxicity, especially in line with reduced degranulation. However, no change was observed in extracellular levels of perforin, indicating that reduced degranulation does not always mean a reduction in effector molecule availability. Nevertheless, our findings led us to conclude that alterations in NK cell effector molecules occur at all three stages, including stores, degranulation and effector release.
Having demonstrated that the release of effector molecules was reduced, we were interested in determining whether a defect in the transport process failed to bring these molecules to the synapse with the target cell. Lytic granules containing effector molecules are mobilized along the NK cell microtubules to position them at the immunological synapse before release. We assessed conjugate formation with the target and polarization of granules and the MTOC to the synapse. The findings demonstrate that the initial phase of conjugation progresses normally while the later stages of granule polarization are impaired. This study extends previous observations of defective granule polarization in NK cells from HIV patients by demonstrating similar defects in NK cells from PHWC-19 (). Although certain tumors resist NK cell killing by preventing granule polarization despite proper conjugation, we found that conjugation was defective, as indicated by the width and morphology of the synapse (). This critical process explains why the cells may not degranulate or release effector proteins, even if their stores were sufficient.
Due to interrupted granule trafficking and impaired synapse formation in NK cells from PHWC-19, we examined whether receptor expression was affected, potentially influencing signaling and cytotoxic responses. NK cell subsets were predominantly skewed towards an inhibited phenotype based on NKG2D, PD-1, LAG3, CTLA4, and CD158b (KIR2DL2/DL3) expression in PHWC-19. Previous studies focused on decreased NKG2D and increased PD-1 expression across all NK cells rather than specific subsets (Varchetta et al., 2021; ; Ustiuzhanina et al., 2023). Our study also investigated LAG3 and CTLA-4, finding increased expression favoring an exhausted phenotype in PHWC-19 (Wilk et al., 2020; Zheng et al., 2020; ). We found significant increases in CD158b in immature cytotoxic CD56dim CD57- NK cells and the immunoregulatory CD56bright CD57- subset (). These changes are similar to human cytomegalovirus, which showed an increase in inhibitory receptors and a reduction in the activating receptors (; ; Tschan-Plessl et al., 2016), but dissimilar to hepatitis C virus, which upregulates the expression of activating receptors, including NKG2D, while lowering the expression of the inhibitory receptor KIR3DL1 (; ). Different viruses will target different receptors to modulate the functioning of NK cells.
By examining inhibitory receptors, we found that subsets of NK cells are more likely to express various families of these receptors in healthy donors and PHWC-19, indicating that individual NK cells expressing a particular receptor are more likely to express another specific receptor. Our study showed that the receptor clusters were similar between healthy donors and PHWC-19 within each NK cell subset, revealing that the clusters were maintained throughout the disease process. However, we also saw that clusters differed between CD57- and CD57+ cell subsets, suggesting that the maturation process alters these clusters. While previous studies have not examined these specific receptors, clusters can form due to differentiation or viral ligand exposure (; ; ). The latter concept states that NK cell diversity may result from interactions with more infected host cells as they mature ().
After examining the mechanisms related to NK cell function, we investigated correlations between demographic and clinical factors and NK cell function in PHWC-19 across different VOCs: Ancestral, Alpha, Delta, and Omicron. While studies have reported that NK cell dysfunction correlates with disease severity, none have examined the impact of different VOCs (Zheng et al., 2020). Despite small sample sizes, we found less severe defects in NK cell function in patients infected with Omicron than Delta and Alpha. The Alpha variant, dominant when vaccine administration was low, potentially resulted in more severe NK defects than in vaccinated Omicron patients. Vaccine-induced antibody-dependent cellular cytotoxicity (ADCC) might compensate for classical pathway defects in NK cell cytotoxicity in Omicron infections (). NK cells are capable of vaccine-induced ADCC, as seen in infections like HIV, malaria, and influenza (; ; ). These observations may contribute to NK cell cytotoxic functions being less impaired in Omicron infections compared to Alpha and Delta (; Wolter et al., 2022; ). Our study is limited in testing this hypothesis because we did not have adequate sample sizes across the various experiments.
Various interventions influencing immune responses were used throughout the pandemic, including vaccination and dexamethasone. We observed no difference in NK cell cytotoxicity between vaccinated and unvaccinated PHWC-19, suggesting that once hospitalized with COVID-19, vaccination status did not influence cytotoxicity. Alternatively, the responses of NK cells from patients who received dexamethasone were compared to those of patients who had not received dexamethasone and to healthy controls treated in vitro with serum concentrations equivalent to oral administration of 6 mg/day. Our findings confirm that defects in NK cell cytotoxicity in PHWC-19 are independent of dexamethasone administration (Witkowski et al., 2021). One study showed that dexamethasone impacts neutrophils in a sex-dependent manner, with males responding better than females during severe COVID-19 (; ). Our study is the first to explore the sex-specific effects of dexamethasone on NK cell function. Previous research indicates that glucocorticosteroids like dexamethasone generally restore immune responses in males but do not improve them in females (; ). NK cells from male PHWC-19 were more responsive to dexamethasone, exhibiting a trend of less defective cytotoxicity than untreated males.
Our study on sex-related differences in NK cell cytotoxicity was motivated by findings that males tend to have more severe COVID-19 outcomes than females (; ). We found a greater defect in NK cell cytotoxicity in males hospitalized with COVID-19 compared to females. Females undergo inactivation of a single X chromosome to balance gene expression from two X chromosomes (). Escaping X inactivation provides women with an advantage in tackling immune challenges because of higher expression levels of immunological factors encoded on the X chromosome, such as TLR7, TLR8, IL2R, CXCR3, and IRAK1 (Wang et al., 2016; ). Similar patterns could be present during COVID-19, where women have shown an upregulation of genes regulating innate antiviral responses, TLR8, TLR9, and IFIT1-3 (). This prompted further investigation into NK cell functional mechanisms. Despite small sample sizes, we observed differences in activating and inhibitory receptor expression. More inhibitory receptors are expressed in CD56dim CD57+ and CD56bright CD57- NK cells from male PHWC-19. These alterations might contribute to the cytotoxicity defect, but further studies are needed to understand the mechanisms behind these defects in male PHWC-19. Further studies would also be required to determine whether these results have implications for long COVID or post-viral remodeling.
COVID-19 disease, in addition to vaccine status and degree of exposure, as public health measures evolved, has continuously changed, making it hard to control extraneous factors. We aimed to collect and investigate various factors, but we were unable to control some elements, such as the timing in the hospital before the blood draw, which could influence NK cell responses, other medications like remdesivir and comorbidities. We tried to address the medications and comorbidities by excluding immunosuppressed individuals and those taking immunosuppressive drugs. Our study did not collect information to determine whether remdesivir influenced our results, although remdesivir has not been shown to affect NK cell functioning in isolation (Zhou et al., 2020; ).
Our study found that NK cells in PHWC-19 have impaired cytotoxic function. We investigated several mechanisms related to NK cell cytotoxicity and discovered that receptor expression skews towards an inhibited phenotype and are co-expressed in specific clusters, possibly affecting granule polarization, effector molecule production, and degranulation. The defect stems from various concurrent malfunctions in mechanisms. Variant and sex-dependent differences were observed, but further studies are needed to understand the underlying mechanisms.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by the Conjoint Health Research Ethics Board of the University of Calgary. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written or oral informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
AD: Project administration, Visualization, Formal analysis, Data curation, Methodology, Writing – review & editing, Conceptualization, Writing – original draft, Investigation. MJP: Conceptualization, Writing – review & editing, Formal analysis, Project administration, Data curation, Writing – original draft, Methodology, Visualization, Investigation. AM: Investigation, Writing – review & editing, Data curation, Formal analysis. DF: Conceptualization, Methodology, Investigation, Writing – review & editing. MP: Investigation, Writing – review & editing. AT: Resources, Project administration, Writing – review & editing. CM: Writing – review & editing, Methodology, Supervision, Funding acquisition, Writing – original draft, Conceptualization, Visualization, Resources, Project administration.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. Canadian Institutes of Health Research Grant numbers 365812 and 447325. Funded material and operating supplies, fees for core facilities, technician and student salaries. Publications costs were funded by The Jesse Bowden Lloyd Professorship the Department of Microbiology, Immunology and Infectious Diseases and C Mody.
Acknowledgments
We gratefully acknowledge Rommy Koetzel for technical support and Zahra Kiani for helpful discussions. Additionally, we would like to thank Jared Schlechte and Anowara Islam for their help with CyTOF and Karen Poon from the Nicole Perkins Microbial Communities (NPMC) Core Labs for her help with flow cytometry.
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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/fcimb.2025.1635043/full#supplementary-material
Supplementary Figure 1Gating strategy for mass cytometry. The HELIOS machine was used for acquisition, and the data was uploaded to the Cytobank program, where manual hierarchical gating was performed after barcoding and labelling with lanthanide–labelled antibodies. All cell populations were initially gated on CD45+ to identify white blood cells via a histogram. Cellular populations were further identified with the following definitions: B cells (CD45+ CD19+), NKT Cells (CD45+ CD3+ CD56+), CD8+ T cells (CD45+ CD3+ CD8+), CD4+ T cells (CD45+ CD3+ CD4+), NK cells (CD45+ NKp46+), and monocytes (CD45+ CD14+). NK cells were further stratified into subsets of CD56dim CD57-, CD56dim CD57+, CD56bright CD57-, and CD56bright CD57 +. No signal was detected below zero, which caused the abrupt cutoff.
Supplementary Figure 2The correlation of granzyme B expression and degranulation with cytotoxicity. (A) 721.221 tumor cell killing by NK cells increases with increasing intracellular granzyme B levels in NK cells from healthy individuals and PHWC-19 (NHealthy=8 and NCOVID=8) (B) 721.221 tumor cell killing by NK cells increases with increasing degranulation in NK cells from healthy individuals and PHWC-19 (NHealthy=10 and NCOVID=10). Circles indicate healthy individuals, while triangles indicate PHWC-19. Experiments were conducted at two effector-to-target ratios (E:T): 5:1 and 10:1. Each point represents the average of three triplicate results. Granzyme B and CD107a levels were measured using fluorescent antibodies and detected by flow cytometry.
Supplementary Figure 3Comparison of time to initial polarization of granules between NK cells from healthy donors and PHWC-19. Time to initial polarization was defined as the time between the start of a sustained contact between the NK cell and the tumor cell and the initial movement of the granules to the synapse. N = 10 conjugates for healthy and PHCW-19. Statistics were calculated using an unpaired T-test after testing for a normal distribution. P<0.05.
Supplementary Figure 4Mass cytometry profiling of NK receptor co-expression on specific CD56bright NK cell populations. (A) Representative t-SNE dimensionality reduction plots of the clusters within the CD56bright CD57- (left) and CD56bright CD57+ (right) NK cell populations from healthy donors (N = 15) and PHWC-19 (N = 17). (B) Pooled relative percentages of each cluster within the CD56bright CD57- (left) and CD56bright CD57+ (right) NK cell populations from healthy donors (N = 15) and PHWC-19 (N = 17) detected by single-cell mass cytometry. The legend applies to sections A and B. (C) Heatmap and dendrogram of log2 expression of activating and inhibitory receptors within the clusters of CD56bright NK cell subset from whole blood of healthy donors (N = 15) and PHWC-19 (N = 17). Data was collected by single-cell mass cytometry. Lighter colors (larger numbers) on the legend indicate greater receptor expression than the darker colors (smaller numbers). Hierarchical clustering heat map was drawn using ChiPlot (https://www.chiplot.online/). * = P<0.05, ** = P<0.01, *** = P<0.001. Statistics for B were calculated using the Mann-Whitney test.
Supplementary Figure 5Principal component analysis (PCA) of receptor expression in NK cell subsets. (A) Comparison of the receptor expression correlations in CD56dim CD57- NK cells in healthy donors and PHWC-19. (B) Comparison of the receptor expression correlations in CD56dim CD57+ NK cells in healthy donors and PHWC-19. (C) Comparison of the receptor expression correlations in CD56bright CD57- NK cells in healthy donors and PHWC-19. (D) Comparison of the receptor expression correlations in CD56bright CD57+ NK cells in healthy donors and PHWC-19. The NK cells were collected from the whole blood of healthy donors (N = 15) and PHWC-19 (N = 17). Data was collected by single-cell mass cytometry.
Supplementary Figure 6Vaccine status and NK cell cytotoxicity in PHWC-19. Comparison of 721.221 tumor cell killing by NK cells from PHWC-19 that were known to be vaccinated to those that were not (NVax=6–7 and NNo Vax=4-5). Experiments were conducted at two effector-to-target ratios (E:T): 5:1 and 10:1.
Supplementary Figure 7Comparisons of sex-related differences in granule-mediated mechanisms. (A) Paired comparison of the expression of intracellular effector molecules within CD56+CD3- PBMCs from healthy individuals and patients with COVID-19 for differences between sexes. The levels of intracellular effector molecules like granzymes A (N = 4), B (N = 10), K (N = 4), perforin (N = 9), and granulysin (N = 10) were measured using fluorescent antibodies and detected by flow cytometry. (B) Comparison of the degranulation levels of intracellular granules from NK cells from healthy individuals and patients with COVID-19 for differences between sexes. The levels of CD107a were tagged using fluorescent antibodies (anti-CD107a) after the NK cells were co-cultured with 721.221 tumor cells and detected by flow cytometry. N = 6. Experiments were conducted at two effector-to-target ratios (E:T): 5:1 and 10:1. (C) Comparison of the level of extracellular effector molecules from NK cells from healthy individuals and patients with COVID-19 for differences between sexes. The granzyme A (N = 6), granzyme B (N = 10), and perforin (N = 4) levels in the supernatants collected from a co-culture of the pNKs with 721.221 tumor cells were detected using ELISA. Experiments were conducted at one effector-to-target ratio (E:T) of 10:1. Males are on the left side, while females are on the right side of each panel. Circles indicate healthy donors, while triangles indicate patients hospitalized with COVID-19. * = P<0.05, ** = P<0.01, *** = P<0.001; calculated using the Kruskal-Wallis ANOVA with Dunn’s multiple comparisons test.
Supplementary Figure 8NK cell cytotoxicity defects were less pronounced in male PHWC-19 treated with dexamethasone than in male PHWC-19 not treated with dexamethasone. Circles indicate healthy donors, and triangles indicate patients hospitalized with COVID-19. To determine the effect of dexamethasone on NK cells from healthy subjects, dexamethasone was added to the cultured media (0.1 μM) (NHealthy=24 and NCOVID=14). Experiments were conducted at two effector-to-target ratios (E:T): 5:1 and 10:1. Each point represents a triplicate result’s average for each subject. Statistics were calculated using the Mann-Whitney test; * = P < 0.05.
Supplementary Table 1Flow cytometry antibody panel.
Supplementary Table 2CYTOF antibody panel.
References
1
AbelA. M.YangC.ThakarM. S.MalarkannanS. (2018). Natural killer cells: development, maturation, and clinical utilization. Front. Immunol.9. doi: 10.3389/fimmu.2018.01869
2
AhlenstielG.TiterenceR. H.KohC.EdlichB.FeldJ. J.RotmanY.et al. (2010). Natural killer cells are polarized toward cytotoxicity in chronic hepatitis C in an interferon-alfa-dependent manner. Gastroenterology138, 325-35.e1–2. doi: 10.1053/j.gastro.2009.08.066
3
AlefishatE.JelinekH. F.MousaM.TayG. K.AlsafarH. S. (2022). Immune response to SARS-CoV-2 variants: A focus on severity, susceptibility, and preexisting immunity. J. Infect. Public Health15, 277–288. doi: 10.1016/j.jiph.2022.01.007
4
AmorimC. F.CarvalhoN. B.NetoJ. A.SantosS. B.GrassiM. F. R.CarvalhoL. P.et al. (2019). The role of NK cells in the control of viral infection in HTLV-1 carriers. J. Immunol. Res.2019, 1–9. doi: 10.1155/2019/6574828
5
ArabiM.Al-NajjarY.MhaimeedN.SalamehM. A.PaulP.AlAnniJ.et al. (2023). Severity of the Omicron SARS-CoV-2 variant compared with the previous lineages: A systematic review. J. Cell Mol. Med.27, 1443–1464. doi: 10.1111/jcmm.17747
6
AroraG.HartG. T.Manzella-LapeiraJ.DoritchamouJ. Y.NarumD. L.ThomasL. M.et al. (2018). NK cells inhibit Plasmodium falciparum growth in red blood cells via antibody-dependent cellular cytotoxicity. Elife7, e36806. doi: 10.7554/eLife.36806
7
BeresfordP. J.ZhangD.OhD. Y.FanZ.GreerE. L.RussoM. L.et al. (2001). Granzyme A activates an endoplasmic reticulum-associated caspase-independent nuclease to induce single-stranded DNA nicks. J. Biol. Chem.276, 43285–43293. doi: 10.1074/jbc.M108137200
8
BéziatV.DuffyD.QuocS. N.Le Garff-TavernierM.DecocqJ.CombadièreB.et al. (2011). CD56brightCD16+ NK cells: A functional intermediate stage of NK cell differentiation. J. Immunol.186, 6753–6761. doi: 10.4049/jimmunol.1100330
9
BéziatV.LiuL. L.MalmbergJ.-A.IvarssonM. A.SohlbergE.BjörklundA. T.et al. (2013). NK cell responses to cytomegalovirus infection lead to stable imprints in the human KIR repertoire and involve activating KIRs. Blood121, 2678–2688. doi: 10.1182/blood-2012-10-459545
10
BozzanoF.DentoneC.PerroneC.Di BiagioA.FenoglioD.ParodiA.et al. (2021). Extensive activation, tissue trafficking, turnover and functional impairment of NK cells in COVID-19 patients at disease onset associates with subsequent disease severity. PloS Pathog.17, e1009448. doi: 10.1371/journal.ppat.1009448
11
CarabelliA. M.PeacockT. P.ThorneL. G.HarveyW. T.HughesJ.de SilvaT. I.et al. (2023). SARS-CoV-2 variant biology: immune escape, transmission and fitness. Nat. Rev. Microbiol. 21, 162–177. doi: 10.1038/s41579-022-00841-7
12
ChenL.JondalM.YakimchukK. (2018). Regulatory effects of dexamethasone on NK and T cell immunity. Inflammopharmacology26, 1331–1338. doi: 10.1007/s10787-017-0418-0
13
ChiveroE. T.BhattaraiN.McLindenJ. H.XiangJ.StapletonJ. T. (2015). Human Pegivirus (HPgV; formerly known as GBV-C) inhibits IL-12 dependent natural killer cell function. Virology485, 116–127. doi: 10.1016/j.virol.2015.07.008
14
ChowdhuryD.LiebermanJ. (2008). Death by a thousand cuts: Granzyme pathways of programmed cell death. Annu. Rev. Immunol.26, 389–420. doi: 10.1146/annurev.immunol.26.021607.090404
15
CianciR.MassaroM. G.De SantisE.TottiB.GasbarriniA.GambassiG.et al. (2023). Changes in lymphocyte subpopulations after remdesivir therapy for COVID-19: A brief report. Int. J. Mol. Sci.24, 14973. doi: 10.3390/ijms241914973
16
ClausM.PierisN.UrlaubD.BrödeP.SchaafB.DurakD.et al. (2023). Early expansion of activated adaptive but also exhausted NK cells during acute severe SARS-CoV-2 infection. Front. Cell Infect. Microbiol.13. doi: 10.3389/fcimb.2023.1266790
17
EitlerJ.WotschelN.MillerN.BoisselL.KlingemannH. G.WelsW.et al. (2021). Inability of granule polarization by NK cells defines tumor resistance and can be overcome by CAR or ADCC mediated targeting. J. Immunother. Cancer9, e001334. doi: 10.1136/jitc-2020-001334
18
Epling-BurnetteP. K.BaiF.PainterJ. S.RollisonD. E.SalihH. R.KruschM.et al. (2007). Reduced natural killer (NK) function associated with high-risk myelodysplastic syndrome (MDS) and reduced expression of activating NK receptors. Blood109, 4816–4824. doi: 10.1182/blood-2006-07-035519
19
ErshovD.PhanM.-S.PylvänäinenJ. W.RigaudS. U.Le BlancL.Charles-OrszagA.et al. (2022). TrackMate 7: integrating state-of-the-art segmentation algorithms into tracking pipelines. Nat. Methods19, 829–832. doi: 10.1038/s41592-022-01507-1
20
FanY.YangB.WuC. (2008). Phenotypically and functionally distinct subsets of natural killer cells in human PBMCs. Cell Biol. Int.32, 188–197. doi: 10.1016/j.cellbi.2007.08.025
21
FarkouhA.BaumgärtelC.GottardiR.HemetsbergerM.CzejkaM.Kautzky-WillerA. (2021). Sex-related differences in drugs with anti-inflammatory properties. J. Clin. Med.10, 1441. doi: 10.3390/jcm10071441
22
Garcia-PeñarrubiaP.KosterF. T.KelleyR. O.McDowellT. D.BankhurstA. D. (1989). Antibacterial activity of human natural killer cells. J. Exp. Med.169, 99–113. doi: 10.1084/jem.169.1.99
23
GoodridgeJ. P.ÖnfeltB.MalmbergK. (2015). Newtonian cell interactions shape natural killer cell education. Immunol. Rev.267, 197–213. doi: 10.1111/imr.12325
24
GrossmanW. J.VerbskyJ. W.TollefsenB. L.KemperC.AtkinsonJ. P.LeyT. J. (2004). Differential expression of granzymes A and B in human cytotoxic lymphocyte subsets and T regulatory cells. Blood104, 2840–2848. doi: 10.1182/blood-2004-03-0859
25
GwalaniL. A.OrangeJ. S. (2018). Single degranulations in NK cells can mediate target cell killing. J. Immunol.200, 3231–3243. doi: 10.4049/jimmunol.1701500
26
HagemannK.RieckenK.JungJ. M.HildebrandtH.MenzelS.BundersM. J.et al. (2022). Natural killer cell-mediated ADCC in SARS-CoV-2-infected individuals and vaccine recipients. Eur. J. Immunol.52, 1297–1307. doi: 10.1002/eji.202149470
27
HassanZ.KumarN. D.ReggioriF.KhanG. (2021). How viruses hijack and modify the secretory transport pathway. Cells10, 2535. doi: 10.3390/cells10102535
28
HaynesB. F.GilbertP. B.McElrathM. J.Zolla-PaznerS.TomarasG. D.AlamS. M.et al. (2012). Immune-correlates analysis of an HIV-1 vaccine efficacy trial. New Engl. J. Med.366, 1275–1286. doi: 10.1056/NEJMoa1113425
29
HeldW.CoudertJ. D.ZimmerJ. (2003). The NK cell receptor repertoire: formation, adaptation and exploitation. Curr. Opin. Immunol.15, 233–237. doi: 10.1016/S0952-7915(02)00031-6
30
HermensJ. M.KesmirC. (2023). Role of T cells in severe COVID-19 disease, protection, and long term immunity. Immunogenetics75, 295–307. doi: 10.1007/s00251-023-01294-9
31
Hernández-BlancoC.Al-Akioui-SanzK.HerreraL.Aguirre-PortolésC.Lozano-OjalvoD.Pérez-RodriguezL.et al. (2025). The phase I RELEASE clinical trial to evaluate the safety of NK cells in COVID-19. iScience28, 111698. doi: 10.1016/j.isci.2024.111698
32
HongH. S.EberhardJ. M.KeudelP.BollmannB. A.BallmaierM.BhatnagarN.et al. (2010). HIV infection is associated with a preferential decline in less-differentiated CD56dim CD16+ NK cells. J. Virol.84, 1183–1188. doi: 10.1128/JVI.01675-09
33
HsuH.-T.CariseyA. F.OrangeJ. S. (2017). “Measurement of lytic granule convergence after formation of an NK cell immunological synapse,” in Methods in Molecular Biology, Humana Press, New York, NY. 497–515. doi: 10.1007/978-1-4939-6881-7_31
34
IslamA.LiS. S.OykhmanP.Timm-McCannM.HustonS. M.StackD.et al. (2013). An Acidic Microenvironment Increases NK Cell Killing of Cryptococcus neoformans and Cryptococcus gattii by Enhancing Perforin Degranulation. PloS Pathog.9, e1003439. doi: 10.1371/journal.ppat.1003439
35
JacobsJ. L.HaidarG.MellorsJ. W. (2023). COVID-19: challenges of viral variants. Annu. Rev. Med.74, 31–53. doi: 10.1146/annurev-med-042921-020956
36
JiangY.WeiX.GuanJ.QinS.WangZ.LuH.et al. (2020). COVID-19 pneumonia: CD8+ T and NK cells are decreased in number but compensatory increased in cytotoxic potential. Clin. Immunol.218, 108516. doi: 10.1016/j.clim.2020.108516
37
KavianN.HachimA.PoonL. L. M.ValkenburgS. A. (2020). Vaccination with ADCC activating HA peptide epitopes provides partial protection from influenza infection. Vaccine38, 5885–5890. doi: 10.1016/j.vaccine.2020.07.008
38
KiesslingR.KleinE.WigzellH. (1975). Natural” killer cells in the mouse. I. Cytotoxic cells with specificity for mouse Moloney leukemia cells. Specificity and distribution according to genotype. Eur. J. Immunol.5, 112–117. doi: 10.1002/eji.1830050208
39
KirwanS.MerriamD.BarsbyN.McKinnonA.BurshtynD. N. (2006). Vaccinia virus modulation of natural killer cell function by direct infection. Virology347, 75–87. doi: 10.1016/j.virol.2005.11.037
40
KrämerB.KnollR.BonaguroL.ToVinhM.RaabeJ.Astaburuaga-GarcíaR.et al. (2021). Early IFN-α signatures and persistent dysfunction are distinguishing features of NK cells in severe COVID-19. Immunity54, 2650–2669.e14. doi: 10.1016/j.immuni.2021.09.002
41
KrenskyA. M.ClaybergerC. (2009). Biology and clinical relevance of granulysin. Tissue Antigens73, 193–198. doi: 10.1111/j.1399-0039.2008.01218.x
42
KyeiS. K.OgbomoH.LiS.Timm-McCannM.XiangR. F.HustonS. M.et al. (2016). Mechanisms by which interleukin-12 corrects defective NK cell anticryptococcal activity in HIV-infected patients. mBio7, 10.1128/mbio.00878-16. doi: 10.1128/mBio.00878-16
43
LeeA. J.MianF.PoznanskiS. M.StackarukM.ChanT.ChewM. V.et al. (2019). Type I interferon receptor on NK cells negatively regulates interferon-γ Production. Front. Immunol.10. doi: 10.3389/fimmu.2019.01261
44
LeeG.SchaunerR.BurkeJ.BoroczJ.VasanaS.SobierajL.et al. (2023). NK cells from COVID-19 positive patients exhibit enhanced cytotoxic activity upon NKG2A and KIR2DL1 blockade. Front. Immunol.14. doi: 10.3389/fimmu.2023.1022890
45
LiM.GuoW.DongY.WangX.DaiD.LiuX.et al. (2020). Elevated exhaustion levels of NK and CD8+ T cells as indicators for progression and prognosis of COVID-19 disease. Front. Immunol.11. doi: 10.3389/fimmu.2020.580237
46
LiebermanJ. (2003). The ABCs of granule-mediated cytotoxicity: New weapons in the arsenal. Nat. Rev. Immunol.3, 361–370. doi: 10.1038/nri1083
47
LinL.LiuY.TangX.HeD. (2021). The disease severity and clinical outcomes of the SARS-coV-2 variants of concern. Front. Public Health9. doi: 10.3389/fpubh.2021.775224
48
LisovskyI.IsitmanG.BruneauJ.BernardN. F. (2015). Functional analysis of NK cell subsets activated by 721.221 and K562 HLA-null cells. J. Leukoc. Biol.97, 761–767. doi: 10.1189/jlb.4AB1014-499R
49
LodoenM. B.LanierL. L. (2006). Natural killer cells as an initial defense against pathogens. Curr. Opin. Immunol.18, 391–398. doi: 10.1016/j.coi.2006.05.002
50
Lopez-VergèsS.MilushJ. M.PandeyS.YorkV. A.Arakawa-HoytJ.PircherH.et al. (2010). CD57 defines a functionally distinct population of mature NK cells in the human CD56dimCD16+ NK-cell subset. Blood116, 3865–3874. doi: 10.1182/blood-2010-04-282301
51
LyonM. F. (1962). Sex chromatin and gene action in the mammalian X-chromosome. Am. J. Hum. Genet.14, 135–148.
52
MacFarlaneA. W.JillabM.SmithM. R.AlpaughR. K.ColeM. E.LitwinS.et al. (2017). NK cell dysfunction in chronic lymphocytic leukemia is associated with loss of the mature cells expressing inhibitory killer cell Ig-like receptors. Oncoimmunology6(7), e1330235. doi: 10.1080/2162402X.2017.1330235
53
MansourabadiA. H.AghamajidiA.DorfakiM.KeshavarzF.ShafeghatZ.MoazzeniA.et al. (2023). B lymphocytes in COVID-19: a tale of harmony and discordance. Arch. Virol.168, 148. doi: 10.1007/s00705-023-05773-y
54
MaoH.TuW.LiuY.QinG.ZhengJ.ChanP.-L.et al. (2010). Inhibition of human natural killer cell activity by influenza virions and hemagglutinin. J. Virol.84, 4148–4157. doi: 10.1128/JVI.02340-09
55
MarshallJ. C.MurthyS.DiazJ.AdhikariN. K.AngusD. C.ArabiY. M.et al. (2020). A minimal common outcome measure set for COVID-19 clinical research. Lancet Infect. Dis.20, e192–e197. doi: 10.1016/S1473-3099(20)30483-7
56
MartinvaletD.ZhuP.LiebermanJ. (2005). Granzyme A induces caspase-independent mitochondrial damage, a required first step for apoptosis. Immunity22, 355–370. doi: 10.1016/j.immuni.2005.02.004
57
MaucourantC.FilipovicI.PonzettaA.AlemanS.CornilletM.HertwigL.et al. (2020). Natural killer cell immunotypes related to COVID-19 disease severity. Sci. Immunol.5, eabd6832. doi: 10.1126/sciimmunol.abd6832
58
MazzoniA.SalvatiL.MaggiL.CaponeM.VanniA.SpinicciM.et al. (2020). Impaired immune cell cytotoxicity in severe COVID-19 is IL-6 dependent. J. Clin. Invest.130, 4694–4703. doi: 10.1172/JCI138554
59
Montaño MendozaV. M.Mendez CortinaY. A.Rodríguez-PereaA. L.FernandezG. J.RugelesM. T.Velilla HernandezP. A.et al. (2023). Biological sex and age-related differences shape the antiviral response to SARS-CoV-2 infection. Heliyon9, e13045. doi: 10.1016/j.heliyon.2023.e13045
60
MorganD. J.DavisD. M. (2017). Distinct effects of dexamethasone on human natural killer cell responses dependent on cytokines. Front. Immunol.8. doi: 10.3389/fimmu.2017.00432
61
MuntasellA.VilchesC.AnguloA.López-BotetM. (2013). Adaptive reconfiguration of the human NK-cell compartment in response to cytomegalovirus: A different perspective of the host-pathogen interaction. Eur. J. Immunol.43, 1133–1141. doi: 10.1002/eji.201243117
62
Nanwani-NanwaniK.López-PérezL.Giménez-EsparzaC.Ruiz-BarrancoI.CarrilloE.ArellanoM. S.et al. (2022). Prevalence of post-intensive care syndrome in mechanically ventilated patients with COVID-19. Sci. Rep.12, 7977. doi: 10.1038/s41598-022-11929-8
63
NielsenC. M.WhiteM. J.GoodierM. R.RileyE. M. (2013). Functional significance of CD57 expression on human NK cells and relevance to disease. Front. Immunol.4. doi: 10.3389/fimmu.2013.00422
64
OgataK.AnE.ShioiY.NakamuraK.LuoS.YokoseN.et al. (2001). Association between natural killer cell activity and infection in immunologically normal elderly people. Clin. Exp. Immunol.124, 392–397. doi: 10.1046/j.1365-2249.2001.01571.x
65
OgbomoH.Timm-McCannM.BarnesT.XiangR. F.JamilK.GangulyA.et al. (2018). Granule-dependent NK cell killing of cryptococcus requires kinesin to reposition the cytolytic machinery for directed cytotoxicity. Cell Rep.24, 3017–3032. doi: 10.1016/j.celrep.2018.08.027
66
OkadaS.LiQ.WhitinJ. C.ClaybergerC.KrenskyA. M. (2003). Intracellular mediators of granulysin-induced cell death. J. Immunol.171, 2556–2562. doi: 10.4049/jimmunol.171.5.2556
67
OlivieroB.VarchettaS.PaudiceE.MicheloneG.ZaramellaM.MavilioD.et al. (2009). Natural killer cell functional dichotomy in chronic hepatitis B and chronic hepatitis C virus infections. Gastroenterology137, 1151–1160.e7. doi: 10.1053/j.gastro.2009.05.047
68
PatoneM.ThomasK.HatchR.TanP. S.CouplandC.LiaoW.et al. (2021). Mortality and critical care unit admission associated with the SARS-CoV-2 lineage B.1.1.7 in England: an observational cohort study. Lancet Infect. Dis.21, 1518–1528. doi: 10.1016/S1473-3099(21)00318-2
69
PengH.XiangT.XuF.JiangY.ZhongL.PengY.et al. (2023). Redistribution and Activation of CD16brightCD56dim NK Cell Subset to Fight against Omicron Subvariant BA.2 after COVID-19 Vaccination. Microorganisms11, 940. doi: 10.3390/microorganisms11040940
70
PiccolellaE.LombardiG.VismaraD.Del GalloF.ColizziV.DoleiA.et al. (1986). Effects of dexamethasone on human natural killer cell cytotoxicity, interferon production, and interleukin-2 receptor expression induced by microbial antigens. Infect. Immun.51, 712–714. doi: 10.1128/iai.51.2.712-714.1986
71
PoliA.MichelT.ThérésineM.AndrèsE.HentgesF.ZimmerJ. (2009). CD56bright natural killer (NK) cells: an important NK cell subset. Immunology126, 458–465. doi: 10.1111/j.1365-2567.2008.03027.x
72
QiS.NgwaC.Morales ScheihingD. A.Al MamunA.AhnstedtH. W.FingerC. E.et al. (2021). Sex differences in the immune response to acute COVID-19 respiratory tract infection. Biol. Sex Differ12, 66. doi: 10.1186/s13293-021-00410-2
73
RichardJ.SindhuS.PhamT. N. Q.BelzileJ.-P.CohenÉ.A. (2010). HIV-1 Vpr up-regulates expression of ligands for the activating NKG2D receptor and promotes NK cell–mediated killing. Blood115, 1354–1363. doi: 10.1182/blood-2009-08-237370
74
RoccaY.PouxvielhK.MarotelM.BenezechS.JaegerB.AllatifO.et al. (2022). Combinatorial expression of NK cell receptors governs cell subset reactivity and effector functions but not tumor specificity. J. Immunol.208, 1802–1812. doi: 10.4049/jimmunol.2100874
75
SchmidtS.ZimmermannS.-Y.TramsenL.KoehlU.LehrnbecherT. (2013). Natural killer cells and antifungal host response. Clin. Vaccine Immunol.20, 452–458. doi: 10.1128/CVI.00606-12
76
SchurzH.SalieM.TrompG.HoalE. G.KinnearC. J.MöllerM. (2019). The X chromosome and sex-specific effects in infectious disease susceptibility. Hum. Genomics13, 2. doi: 10.1186/s40246-018-0185-z
77
SeeD. M.KhemkaP.SahlL.BuiT.TillesJ. G. (1997). The role of natural killer cells in viral infections. Scand. J. Immunol.46, 217–224. doi: 10.1046/j.1365-3083.1997.d01-121.x
78
SinhaS.RosinN. L.AroraR.LabitE.JafferA.CaoL.et al. (2022). Dexamethasone modulates immature neutrophils and interferon programming in severe COVID-19. Nat. Med.28, 201–211. doi: 10.1038/s41591-021-01576-3
79
SolanaR.MarianiE. (2000). NK and NK/T cells in human senescence. Vaccine18, 1613–1620. doi: 10.1016/S0264-410X(99)00495-8
80
SolanaR.PawelecG.TarazonaR. (2006). Aging and innate immunity. Immunity24, 491–494. doi: 10.1016/j.immuni.2006.05.003
81
SongD.DuBoisD. C.AlmonR. R.JuskoW. J. (2018). Modeling sex differences in anti-inflammatory effects of dexamethasone in arthritic rats. Pharm. Res.35, 203. doi: 10.1007/s11095-018-2483-5
82
SpoorenbergS. M. C.DeneerV. H. M.GruttersJ. C.PullesA. E.VoornG. P. P.RijkersG. T.et al. (2014). Pharmacokinetics of oral vs intravenous dexamethasone in patients hospitalized with community-acquired pneumonia. Br. J. Clin. Pharmacol.78, 78–83. doi: 10.1111/bcp.12295
83
TakahashiT.EllingsonM. K.WongP.IsraelowB.LucasC.KleinJ.et al. (2020). Sex differences in immune responses that underlie COVID-19 disease outcomes. Nature588, 315–320. doi: 10.1038/s41586-020-2700-3
84
The Recovery Collaborative Group (2021). Dexamethasone in hospitalized patients with covid-19. New Engl. J. Med.384, 693–704. doi: 10.1056/NEJMoa2021436
85
ThieryJ.LiebermanJ. (2014). Perforin: A key pore-forming protein for immune control of viruses and cancer. Subcell Biochem.80, 197–220. doi: 10.1007/978-94-017-8881-6_10
86
TrapaniJ. A.BirdP. I. (2008). A renaissance in understanding the multiple and diverse functions of granzymes? Immunity29, 665–667. doi: 10.1016/j.immuni.2008.10.002
87
Tschan-PlesslA.SternM.SchmiedL.RetièreC.HirschH. H.GarzoniC.et al. (2016). Human cytomegalovirus infection enhances NK cell activity in vitro. Transplant. Direct2, e89. doi: 10.1097/TXD.0000000000000605
88
UstiuzhaninaM. O.VavilovaJ. D.BoykoA. A.StreltsovaM. A.KustS. A.KanevskiyL. M.et al. (2023). Coordinated loss and acquisition of NK cell surface markers accompanied by generalized cytokine dysregulation in COVID-19. Int. J. Mol. Sci.24, 1996. doi: 10.3390/ijms24031996
89
van DomselaarR.BovenschenN. (2011). Cell death-independent functions of granzymes: hit viruses where it hurts. Rev. Med. Virol.21, 301–314. doi: 10.1002/rmv.697
90
van ErpE. A.FeyaertsD.DuijstM.MulderH. L.WichtO.LuytjesW.et al. (2019). Respiratory syncytial virus infects primary neonatal and adult natural killer cells and affects their antiviral effector function. J. Infect. Dis.219, 723–733. doi: 10.1093/infdis/jiy566
91
VarchettaS.MeleD.OlivieroB.MantovaniS.LudovisiS.CerinoA.et al. (2021). Unique immunological profile in patients with COVID-19. Cell Mol. Immunol.18, 604–612. doi: 10.1038/s41423-020-00557-9
92
VivierE.TomaselloE.BaratinM.WalzerT.UgoliniS. (2008). Functions of natural killer cells. Nat. Immunol.9, 503–510. doi: 10.1038/ni1582
93
WangJ.SyrettC. M.KramerM. C.BasuA.AtchisonM. L.AngueraM. C. (2016). Unusual maintenance of X chromosome inactivation predisposes female lymphocytes for increased expression from the inactive X. Proc. Natl. Acad. Sci.113, E2029-38. doi: 10.1073/pnas.1520113113
94
WeijtensO.SchoemakerR. C.CohenA. F.RomijnF. P. H. T. M.LentjesE. G. W. M.van RooijJ.et al. (1998). Dexamethasone concentration in vitreous and serum after oral administration. Am. J. Ophthalmol.125, 673–679. doi: 10.1016/S0002-9394(98)00003-8
95
WilkA. J.RustagiA.ZhaoN. Q.RoqueJ.Martínez-ColónG. J.McKechnieJ. L.et al. (2020). A single-cell atlas of the peripheral immune response in patients with severe COVID-19. Nat. Med.26, 1070–1076. doi: 10.1038/s41591-020-0944-y
96
WitkowskiM.TizianC.Ferreira-GomesM.NiemeyerD.JonesT. C.HeinrichF.et al. (2021). Untimely TGFβ responses in COVID-19 limit antiviral functions of NK cells. Nature600, 295–301. doi: 10.1038/s41586-021-04142-6
97
WolterN.JassatW.WalazaS.WelchR.MoultrieH.GroomeM.et al. (2022). Early assessment of the clinical severity of the SARS-CoV-2 omicron variant in South Africa: a data linkage study. Lancet399, 437–446. doi: 10.1016/S0140-6736(22)00017-4
98
World Health Organization (2024). WHO COVID-19 dashboard: Number of COVID-19 deaths reported to WHO (Geneva, CH: WHO Health Emergencies Programme). Available online at: https://data.who.int/dashboards/covid19/deaths?n=c (Accessed March 29, 2024).
99
ZenarruzabeitiaO.Astarloa-PandoG.TerrénI.OrrantiaA.Pérez-GarayR.Seijas-BetolazaI.et al. (2021). T cell activation, highly armed cytotoxic cells and a shift in monocytes CD300 receptors expression is characteristic of patients with severe COVID-19. Front. Immunol.12. doi: 10.3389/fimmu.2021.655934
100
ZhaoT.ZhangH.GuoY.ZhangQ.HuaG.LuH.et al. (2007). Granzyme K cleaves the nucleosome assembly protein SET to induce single-stranded DNA nicks of target cells. Cell Death Differ14, 489–499. doi: 10.1038/sj.cdd.4402040
101
ZhengM.GaoY.WangG.SongG.LiuS.SunD.et al. (2020). Functional exhaustion of antiviral lymphocytes in COVID-19 patients. Cell Mol. Immunol.17, 533–535. doi: 10.1038/s41423-020-0402-2
102
ZhouL.HuntingtonK.ZhangS.CarlsenL.SoE.-Y.ParkerC.et al. (2020). Natural Killer cell activation, reduced ACE2, TMPRSS2, cytokines G-CSF, M-CSF and SARS-CoV-2-S pseudovirus infectivity by MEK inhibitor treatment of human cells. doi: 10.1101/2020.08.02.230839
103
ZiehrD. R.AlladinaJ.PetriC. R.MaleyJ. H.MoskowitzA.MedoffB. D.et al. (2020). Respiratory pathophysiology of mechanically ventilated patients with COVID-19: A cohort study. Am. J. Respir. Crit. Care Med.201, 1560–1564. doi: 10.1164/rccm.202004-1163LE
Summary
Keywords
COVID-19, cytotoxic effector molecules, granule trafficking, natural killer (Nk) cell, cellular cytotoxicity, immunology, sex differences, SARS-CoV-2
Citation
Dagar A, Polyak MJ, Mok ACH, Feehan D, Potemkin M, Tremblay A and Mody CH (2025) The effect of COVID-19 and sex differences on natural killer cell cytotoxicity. Front. Cell. Infect. Microbiol. 15:1635043. doi: 10.3389/fcimb.2025.1635043
Received
25 May 2025
Accepted
05 September 2025
Published
22 September 2025
Volume
15 - 2025
Edited by
John Hiscott, Istituto Pasteur Italia Cenci Bolognetti Foundation, Italy
Reviewed by
Kai Qin, Fred Hutchinson Cancer Center, United States
Bingyu Yan, Stanford University, United States
Updates
Copyright
© 2025 Dagar, Polyak, Mok, Feehan, Potemkin, Tremblay and Mody.
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: Christopher H. Mody, cmody@ucalgary.ca
Disclaimer
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.