Progressive accumulation of hyperinflammatory NKG2Dlow NK cells in early life defines a novel endotype of severe atopic dermatitis
==================================================================================================================================

* David E. Ochayon
* Stanley B. DeVore
* Wan-Chi Chang
* Durga Krishnamurthy
* Harsha Seelamneni
* Brittany Grashel
* Daniel Spagna
* Sandra Andorf
* Lisa J. Martin
* Jocelyn M. Biagini
* Stephen N. Waggoner
* Gurjit K. Khurana Hershey

## Abstract

Atopic dermatitis (AD) is a chronic inflammatory skin disease that often precedes the development of food allergy, asthma, and allergic rhinitis. The prevailing paradigm holds that a reduced frequency and function of natural killer (NK) cell contributes to AD pathogenesis, yet the underlying mechanisms and contributions of NK cells to allergic co-morbidities remain ill-defined. Herein, analysis of circulating NK cells in a longitudinal early life cohort of children with AD revealed a progressive accumulation of NK cells with low expression of the activating receptor NKG2D, which was linked to more severe AD and sensitivity to allergens. This was most notable in children co-sensitized to food and aero allergens, a risk factor for development of asthma. Individual-level longitudinal analysis in a subset of children revealed co-incident reduction of NKG2D on NK cells with acquired or persistent sensitization, and this was associated with impaired barrier function. Low expression of NKG2D on NK cells was paradoxically associated with depressed cytolytic function but exaggerated release of the proinflammatory cytokine TNF-α. These observations provide important new insights into a potential pathophysiological mechanism of atopic march involving altered NK-cell functional responses and define a novel endotype of severe AD.

*   innate lymphoid cells
*   eczema
*   atopic march
*   skin
*   hypersensitivity
*   cytotoxicity
*   innate immunity
*   tumor necrosis factor
*   allergic rhinitis
*   epidermal
*   phenotype

## Introduction

Atopic dermatitis (AD) is an inflammatory skin disease affecting up to 25% of children that often precedes the development of other allergic disorders, including food allergy and allergic rhinitis (*1, 2*). Moreover, nearly half of pediatric patients with AD will develop asthma (*3*). The mechanisms by which AD favors progression to allergies and asthma remain unclear.

Type 2 cytokines (e.g., IL-13) and effector cells (e.g., eosinophils) are generally considered drivers of AD and allergic disease, whereas type 1 immune cells such as natural killer (NK) cells are implicated in a counter regulatory role (*4*). Consistent with this, several studies have contributed a paradigm holding that reduced frequencies and/or diminished cytolytic functionality of circulating NK cells contributes to the pathogenesis of AD (*5-8*). Indeed, IL-15-based therapies that can expand NK cells provided some measure of benefit in mouse models of AD (*5*). NK cells are thought to limit eczema and allergic inflammation via IFN-*γ* counteraction of type 2 inflammation (*9*) and/or cytolytic regulation of type 2 immune cells (*10-12*). Nevertheless, transcriptomic analyses reveal a progressive enrichment for NK cells in lesional skin of patients with atopic dermatitis (*13-15*). Potential changes in NK-cell function that may contribute to development of allergic disease in children with AD have not been explored.

The Mechanisms of Progression of Atopic Dermatitis to Asthma in Children (MPAACH) cohort (N=651), the first US based prospective early life cohort of children with AD, was designed to elucidate the mechanistic underpinnings of AD and the atopic march (*16-19*). MPAACH is a longitudinal study where children are seen for annual visits that include questionnaires, a physical examination, and extensive biospecimen collections. Herein, we investigated peripheral NK cells in a subset of MPAACH children longitudinally over sequential annual visits to define phenotypic and functional changes that are related to clinical outcomes. Our findings reveal that more severe AD and sensitivity to allergens is associated with selective loss of NKG2D on NK cells but not other immune cells. Correspondingly, there is a progressive accumulation of NKG2Dlow CD56dim NK cells that exhibit reduced cytolytic function but elevated secretion of proinflammatory cytokines such as TNF-α. Thus, accumulation of hyperinflammatory NKG2Dlow NK cells defines an endotype of severe AD characterized by increased sensitivity to allergens.

## Results

### Altered repertoire of NK cells in allergen-sensitized children with AD

Peripheral blood samples from any of the first 3 annual visits of a randomly selected subset (n=124, 39-44 subjects per visit) of MPAACH subjects were utilized to assess NK cells. The demographic and clinical characteristics of the included subjects in comparison to the full MPAACH cohort are shown in **Table 1**. Compared to the full cohort, the randomly selected subset was slightly older (2.2 vs. 2.0 y; *p* < 0.001), less likely to be black (57% vs. 68%; *p* = 0.037), more likely to be co-sensitized to both food and aeroallergens (24% vs. 14%; *p* = 0.003), and more likely to be diagnosed with food allergy (23% vs. 14%; *p* = 0.014). The incidence of co-morbid conditions such as wheezing and allergies/allergic rhinitis significantly increased at later clinical visits in this cohort of children (**Supplemental Table 1**). **Table 1** also shows a comparison of complete blood counts (CBC) in the overall MPAACH cohort and sub-cohort evaluated in this study. The frequencies of different peripheral immune cells were similar in the overall and sub-study cohorts (**Table 1**). Comparisons of cohort demographics and CBCs by visit are shown in **Supplemental Tables 1 and 2**, respectively.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/T1)

Table 1. 
Characteristics of the MPAACH cohort and the Study Subset at Visit I1

Reduced frequencies of NK cells are commonly observed in the blood of patients with AD (*5-8*), but relationships between NK-cell frequencies and co-morbidities of AD are undefined. As expected, most NK cells in the peripheral blood of children with AD display a cytolytically mature CD16+CD56dim phenotype, with less common populations of more immature CD16lowCD56bright and CD16+CD56neg NK cells (**Supplemental Figure 1A+B**). Of note, the frequency of total NK cells (**Figure 1A**) and the ratio of CD56dim to CD56bright NK cells (**Figure 1B**) in blood was similar (p>0.1) in children with AD across all time points analyzed regardless of allergic sensitivity.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F1)

Figure 1. Altered repertoire of NK cells in allergen sensitized children with AD.
PBMC from children with AD who lack or exhibit allergic sensitivity (V1=44, V2=39, V3, n=41) were analyzed by flow cytometry at different annual clinical visits to ascertain (**A**) the percentages of CD56+ NK cells among total lymphocytes and (**B**) ratio of CD16+CD56dim to CD16lowCD56bright NK cell subsets. Differences in median values between sensitized and non-sensitized children were not statistically significant in a Mann-Whitney U test (p<0.1). (**C-F**) NK cells from V2 or V3 for 4 non-sensitized (grey) and 4 sensitized (purple) patients were subjected to (**C**) UMAP dimensional analysis and (**D**) FlowSOM clustering. (**E**) marker expression of 8 distinct subsets of NK cells within the UMAP distribution. (**F**) The proportional contribution of NK cells from sensitized and non-sensitized children to each FlowSOM cluster, Standard error of means within each cluster were compared using Mann-Whitney U test, p <0.1 values are reported.

Despite these similar frequencies of NK cells, changes in the expression levels of germline encoded NK-cell receptors in different individuals could contribute to altered functional responses to inflammatory or environmental (e.g., allergen) cues. We used flow cytometry and dimensionality reduction analysis (UMAP) in conjunction with unsupervised clustering (FlowSOM) to project differences in NK-cell population phenotypes in a subset (described in **Supplemental Table 3**) of allergen sensitized (n=4, purple) and non-sensitized (n=4, gray) children with AD based on the expression of CD16, NKG2D, CD69, CD62L, CCR4, Tim-3, NKp44, and CD45RO on NK cells (**Figure 1C**). We chose these receptors based on evidence that their expression changes or contributes to NK-cell functional activity in allergic or asthmatic disease (*12, 24-28*). Eight distinct populations of NK cells were apparent in this FlowSOM analysis (**Figure 1D+E** and **Supplemental Figure 2**). The CD16dim CD56bright NK cells cluster (#7, **Figure 1D**) was similarly (*p*>0.1) populated by cells from the sensitized and non-sensitized individuals (**Figure 1F**), consistent with conventional flow cytometry analysis (**Figure 1B**). Among the eight different NK-cell clusters, only two populations (Clusters 1 and 2) showed significant enrichment (*p*=0.029) in non-sensitized compared to sensitized children (**Figure 1F**). These two clusters were comprised of CD16+CD56dim NK cells characterized by high expression levels of NKG2D, low expression of CCR4 and CD62L, and elevated expression of CD69 (**Figure 1E)**. These results suggest that allergic sensitization in patients with AD is associated with alteration of the circulating NK cell repertoire.

### Expression of NKG2D on NK cells is inversely associated with severity of AD and with allergic co-morbidities

The diminished presence of NKG2D+ NK cells in this subset of allergen-sensitized individuals (**Figure 1**) may represent a specific early-life endotype of children with AD that reflects altered functional responsiveness to activating signals via NKG2D (*29*). Indeed, the reduced expression levels of NKG2D on CD56dim NK cells in sensitized compared to non-sensitized MPAACH children (**Figure 2A+B**) facilitated stratification and quantification of NKG2Dlow and NKG2Dhigh NK cells (**Supplemental Figure 1C**). The phenotype of reduced NKG2D expression was unique to NK cells, as expression levels of NKG2D on non-NK cells (CD3+ or CD14+ or CD19+) did not differ between the sensitized and non-sensitized groups (**Figure 2B**). The percentages of NKG2Dlow CD56dim (**Figure 2C**) NK cells were increased in allergen-sensitized patients relative to non-sensitized children with AD at V2 (*p*=0.024) and V3 (*p*=0.007). CD56bright NK cells showed a similar trend that did not reach statistical significance (**Figure 2D**). The magnitude of difference in NKG2D expression levels was greatest at V3 (**Figure 2C+D**, far right plots), suggesting a progressive aggravation of the NKG2D pathway.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F2)

Figure 2. Progressive longitudinal decline of NKG2D expression on CD56dim NK in allergen sensitized children.
(**A**) Representative plot of NKG2Dlow and NKG2Dhigh proportions on CD56dim NK cells determined by flow cytometry. (**B**) Mean fluorescence intensity (MFI) of NKG2D expression on CD56dim NK cells and non-NK cells in non-sensitized (closed black circles, n = 18) and sensitized (open red circles, n = 23) children at V3. Proportion of (**C**) CD56dim and (**D**) CD56bright NK cells gated as NKG2Dlow (defined in A) at V1, V2 and V3. Statistical differences between median values of non-sensitized (closed black circles, V1 n=23; V2, n=16; V3, n = 18) and sensitized (open red circles, V1, n =21; V2, n=23; V3, n =23) determined by two-tailed Mann-Whitney U test, p<0.1 values are reported.

Expression levels of the activating Fc receptor CD16 showed a similar progressive decline on CD56dim NK cells in sensitized MPAACH children relative to those who are non-sensitized (**Supplemental Figure 3A**). There were no apparent trends or differences in the expression of CCR4, CD62L, and Tim-3 by CD56dim NK cells between the two groups of children over time, while CD69 exhibited a consistent, non-significant reduction in expression on NK cells in sensitized compared to non-sensitized children (**Supplemental Figure 3B-D**). Similar non-significant trends were observed across the range of receptors analyzed in both CD56bright (**Supplemental Figure 4**) and CD56neg (**Supplemental Figure 5**) NK cells between sensitized and non-sensitized children.

Given the increasing trajectory of NKG2Dlow NK cells over time in children with AD who had persistent or developed allergic sensitization, we examined the association of this NK-cell phenotype with the severity of AD. Sensitized children exhibited significantly higher SCORAD than non-sensitized children (25.7 vs 19.7, *p* < 0.001). Moreover, there was a positive association (r=0.316, *p*=0.010) between the proportions of NKG2Dlow CD56dim NK cells and SCORAD in sensitized children (**Figure 3**) that was not observed in non-sensitized children (r=-0.148, *p*=0.27). Thus, reduced NKG2D expression on CD56dim NK cells is linked to more severe eczema in allergen-sensitized children.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F3)

Figure 3. Association between NKG2Dlow percentages to increased AD severity in sensitized children.
Spearman correlation test showed inverse relationships between SCORAD and proportion of NKG2Dlow in sensitized (open red circles, n=67) but not in non-sensitized (closed black circles, n=57) children (Combined from V1-V3) in CD56dim NK cells.

### NK cell phenotypic changes are most evident in children with co-sensitization to food and aero allergens

The type of allergic sensitization is important in both clinical outcome and treatment decision making. Specifically, co-sensitization to multiple allergens is a stronger risk factor for development of rhinitis and asthma compared to mono-sensitization (*30-32*). Since reduced NKG2D (**Figure 2**) and CD16 (**Supplemental Figure 3A**) expression were characteristics of NK cells of AD patients with allergen sensitivity, we compared the expression of these receptors on NK cells in patients with sensitivity to one type of allergen (food or aero mono-sensitization) to patients with sensitization to at least one food and one aero allergen (co-sensitization). Similar to observations in **Figure 2**, sensitization type had no effect on proportions of NKG2Dlow CD56dim NK cells at V1 (**Figure 4A**, left column). A significant enrichment of NKG2Dlow CD56dim NK cells was observed at V2 (*p*=0.013) and V3 (*p*=0.007) in co-sensitized individuals but was not evident until V3 (*p*=0.025) in mono-sensitized children (**Figure 4A**, middle and right columns). A similar trend was evident for CD16 expression, with significantly (*p*<0.05) reduced expression levels of CD16 observed on CD56dim NK cells in co-sensitized individuals at every clinical timepoint (**Figure 4B**). This effect was not observed until V2 and was of a reduced magnitude in mono-sensitized children (**Figure 4B**). Thus, co-sensitization is more strongly linked to reduced activating receptor expression on NK cells.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F4)

Figure 4. Reduced NK cell activating receptor expression in CD56dim NK cells exacerbated by allergen co-sensitization.
Flow cytometry was used to determine the proportion of (**A**) NKG2Dlow CD56dim NK cells or the (**B**) Mean fluorescence intensity of CD16 expression on CD56dim NK cells at each clinical visit for children with AD who are non-sensitized (closed black circles. V1, n=23, V2, n=16, V3, n = 18), sensitized to either food or aero allergens (monosensitization, open red circles. V1, n=14; V2, n=12; V3, n = 11), or sensitized to both food and aero allergens (co-sensitization, closed red circles. V1, n=7; V2, n=11; V3, n = 12). Significant differences between group medians were determined by one-tailed Mann-Whitney U test based on expectation of singular direction of effect in Figure 2, p<0.1 values are reported.

### NK cells with reduced NKG2D expression display decreased cytolytic function

Numerous studies measuring cytolytic activity of peripheral NK cells against the K562 cells have shown that this activity is reduced in individuals with AD relative to healthy controls (*5-8*). Given the importance of NKG2D to NK cell killing of K562 cells (*33*), we evaluated whether cytolytic function would be more substantially impacted in sensitized individuals whose NK cells exhibit lower NKG2D expression. Prior to use in killing assays, PBMC were either unstimulated (**Figure 5A**) or activated with IL-12 and IL-18 (**Figure 5B**), which are established stimulants of NK cell activity (*34*) that are implicated in the pathogenesis of AD (*35-37*) and associated allergic diseases (*38-40*). Spontaneous killing activity observed in unstimulated NK cells from sensitized individuals was significantly reduced compared to NK cells from non-sensitized children (*p*=0.016), which were similar to children without AD (non-AD, **Figure 5A**). In addition, the capacity of IL-12 and IL-18 to amplify NK cell cytolytic activity was reduced more than two-fold (*p*<0.0001) in sensitized children relative to non-sensitized children with AD and non-AD controls (**Figure 5B**). The CD56dim NK cells from the sensitized individuals used in the functional assays depicted in Figure 5, were confirmed to exhibit lower NKG2D expression compared to non-AD and non-sensitized children (**Supplemental Figure 6**). Thus, reduced expression of NKG2D in sensitized children with AD is linked to depressed cytolytic functionality of NK cells.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F5)

Figure 5. Reduced cytolytic function but enhanced NKG2D-induced cytokine production associated with low NKG2D expression.
**(A-B)** Peripheral blood (V1 or V2) from children without AD (closed gray circles, n=4) or children with AD with allergen sensitivity (closed red circles,V1 or V2, n=5) or non-sensitized (closed black circles, V1 or V2, n=5) were stimulated with (**A**) media or (**B**) media containing IL-12 and IL-18 prior to 4-hour co-culture with 51Cr-labeled K562 cells. Chromium release in culture supernatant was measured using a gamma reader. Comparison of the area under the curve was performed using ANOVA. (**C-D**) V3 Blood cells from n=21 patients were assigned into two groups based on (**C**) the median of mean fluorescence intensity of NKG2D expression on CD56dim NK cells and (**D-E**) IFN-γ and TNF-α expression was measured by flow cytometry following stimulation with plate bound MICA (10 μg/ml). Mann-Whitney U test was used to determine significant difference between the groups with higher (closed black squares, n=10) and lower (open white squares, n=11) median NKG2D expression level, p<0.1 values are reported.

### Peripheral NK cells from sensitized children with AD express higher levels of inflammatory cytokines

To directly assess the impact of reduced NKG2D expression on cytokine production by NK cells, we stimulated a subset of PBMC with plate-bound MICA, a known ligand for NKG2D (*29*). We utilized V3 samples (n=21) stratified as NKG2Dhigh or NKG2Dlow (**Figure 5C**) based on the observed expression level of NKG2D on NK cells by MFI (NKG2Dlow: below median; NKG2Dhigh: above median). On average, patient-derived NK cells with lower than the median of NKG2D MFI responded to stimulation with MICA with greater IFN-γ (**Figure 5D**) and TNF-α (**Figure 5E**) than patient-derived NK cells with higher than the median of NKG2D MFI. In addition, blood CD56dim NK cells from allergen sensitized children with AD at V1 and V2 exhibited 45% greater spontaneous expression of TNF-α than non-sensitized children even in the absence of *ex vivo* restimulation (**Figure 6A+B**). A similar trend for increased TNF-α expression was observed in CD56bright NK cells (**Figure 6C+D**). The expression of TNF-α remained elevated in these NK cell populations following stimulation with IL-12 and IL-18 (**Figure 6B+D**). In contrast, IFN-γ expression by peripheral NK cells was only evident following stimulation with IL-12 and IL-18 (**Supplemental Figure 7**) and was not significantly different amongst the groups. We examined NKG2D expression on the NK cells examined in Figure 6 and there was no observed association between NKG2D expression and inflammatory cytokine expression (**Supplemental Figure 8**). Interestingly and, consistent with progressive aggravation of NKG2D expression at later clinical visits (**Figure 2**), the MFI of NKG2D expression was inversely correlated with TNF-α (r=-0.44, *p=*0.048) production in samples derived from V2 but not those from V1 (r=0.07 *p*=0.77).

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F6)

Figure 6. CD56dim and CD56bright NK cells in allergen-sensitized patients show increased expression of TNF-α.
Percent TNFα-producing NK cells following culture of V1 and V2 peripheral blood (**A-B**) CD56dim or (**C-D**) CD56bright NK cells in media in the absence or presence of IL-12 and IL-18. NK cells derived from PBMC from MPAACH children who were: non-AD (closed gray circles, n = 7), non-sensitized (closed black circles, n=15), or allergen sensitized (open red circles, n=26) were evaluated for cytokine expression by flow cytometry. Statistical differences between group medians were determined by Mann-Whitney U test, p<0.1 values are reported.

### Individual-level longitudinal analysis of NKG2DlowCD56dim and NKG2DlowCD56bright NK cells reveals their association with persistent and acquired sensitization

In a subset of 10 children for whom individual-level longitudinal NK cell and sensitization data was available across multiple time points, we examined NK cells and allergic sensitization (non/transient vs. acquired/persistent) within each child over time. There were no differences in age, sex, race or public insurance in those who were non/transiently sensitized compared to those that were acquired/persistently sensitized (**Supplemental Table 4)**. At V3, children with acquired or persistent sensitization had increased proportions of NKG2DlowCD56dim cells (*p=*0.033, **Figure 7A**) compared to children who were non-sensitized or had transient sensitization, and a similar trend was observed in NKG2DlowCD56bright NK cells (*p=*0.057, **Figure 7B**). Over time from V1/2 to V3, children with acquired/persistent sensitization had a mean increase of 28.6% in NKG2DlowCD56dim cells and 29.6% in NKG2DlowCD56bright NK cells compared to only 2.7% and 2.9% for the non/transiently sensitized children (*p=*0.13 and 0.06, respectively (**Figure 7 C+D**). We next evaluated whether the change in the percentage of NK cells over time from V1/V2 to V3 was associated with AD severity and/or skin barrier function. Lesional TEWL was significantly positively associated with an increase in both NKG2DlowCD56dim cells (r=0.59, *p=*0.09, **Supplemental Figure 9A**) and CD56bright NK cells (r=0.55, *p=*0.10, **Supplemental Figure 9B**). When we stratified by sensitization phenotype, we observed that these associations were driven by the children with longitudinal acquired/persistent sensitization (**Supplemental Figure 9 C+D**). There were no associations between NK cell proportion over time and SCORAD (**Supplemental Figure 9E-H**). Collectively, these data suggest that the progressive enrichment in NKG2Dlow NK cells over time is co-incident with worsening skin barrier function and persistent or late sensitization.

![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F7.medium.gif)

[Figure 7.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F7)

Figure 7. Individual-level longitudinal analysis of proportions of NKG2Dlow CD56dim and CD56bright NK cells in MPAACH participants with different longitudinal sensitization patterns as indicated.
NKG2Dlow measurements in participants at V1/V2 and V3 were plotted individually to provide summary of non-(NS-NS) or transient (S-NS) sensitized (closed black circles) vs persistent (S-S) or acquired sensitized (NS-S) (open red circles) at V3 for (**A**) CD56dim and (**B**) CD56bright NK cells. Significant differences between group medians were determined by one-tailed Mann-Whitney U test based on expectation of singular direction of effect in Figure 2. Trajectories of proportions of (**C**) NKG2DlowCD56dim and (**D**) NKG2Dlow CD56bright NK cells within individual children over time who have the following longitudinal sensitization patterns: non-sensitization (NS-NS, closed black circles), transient sensitization (S-NS, open black circles), acquired sensitization (NS-S, closed red circles) and persistent sensitization (S-S, open red circles) were analyzed using paired t-test.

## Discussion

Numerous publications over the past four decades have compared patients with AD to healthy controls and established a paradigm that loss of NK-cell activity contributes to pathogenesis of eczema. As a result, several translational studies are focused on therapeutic expansion or restoration of NK cells in patients with AD. Our analyses in a well-defined longitudinal cohort of children with AD surprisingly revealed a progressive accumulation of hyperinflammatory NKG2Dlow NK cells that is associated with more severe AD and allergen sensitivity. The loss of NKG2D expression by NK cells was more pronounced in patients co-sensitized to both aero and food allergens, a known risk factor for development of asthma. An individual-level analysis of longitudinal paired sensitization and NK cell data over time further supported these fundings. Strikingly, there was a marked increase in the proportion of NKG2Dlow NK cells over time in children who were persistently sensitized or acquired sensitization over time in contrast to those who were non-sensitized or only transiently sensitized early. Further, this magnitude of the increase in NKG2Dlow NK cells was correlated with skin barrier dysfunction in these individuals. Collectively, our data indicate that dysfunctional hyperinflammatory NKG2Dlow NK cells accumulate co-incidentally with worsening sensitization and increasing barrier dysfunction. These observations provide important new insights into a potential pathophysiological mechanism of atopic march involving altered NK cell functional responses.

Atopic dermatitis and associated allergic co-morbidities are correlated with heightened expression of type 2 cytokines, including IL-4 (*41*). Therefore, one potential mechanism of reduced NKG2D expression and diminished cytolytic function of NK cells relates to potential effects of the type 2 cytokines on NK cells. Indeed, IL-4 is capable of down-regulating NKG2D expression on NK cells (*42*). In contrast, type 1 cytokines like IL-12 can promote increased expression of NKG2D (*43*), although this may be counteracted by suppressive effects of IL-4 on IL-12 responses (*44*). A potential source of type 2 cytokines are T cells (*45*). A recent study by Zhou et al. dysregulation of NK-cell surface markers expression and function was shown to be downstream of Th2 allergic T cells responses (*46*). Whether other type 2 cytokines such as IL-13, IL-5, IL-33, and TSLP can influence expression of NKG2D remains to be determined.

A second mechanism potentially contributing to downregulation of NKG2D on NK cells in allergen A second mechanism potentially contributing to downregulation of NKG2D on NK cells in allergen sensitized patients with AD lies in the effects of allergens themselves. Children sensitized to house dust mite exhibit increased sera levels of soluble ligands (e.g. MICA and ULBP2) for NKG2D compared to non-sensitized children (*47*). One effect of elevated expression of soluble NKG2D ligands is the downregulation of NKG2D expression on NK cells (*48, 49*). Thus, reduced NKG2D expression in allergen sensitized patients might be caused by allergen-induced increases in the levels of circulating soluble NKG2D ligands. Allergens often also possess proteolytic activity. Since metalloproteases like ADAM10 play a role in cleavage of both NKG2D ligands and the receptor CD16 from the surface of NK cells (*50, 51*), similar functionality of allergen-associated proteases or proteolytic activation of endogenous metalloproteinases may contribute to loss of NKG2D and CD16 expression. These observations may also explain the exacerbation of NKG2D loss in co-sensitized patients, who respond to a greater number of potential allergens and increased frequency of allergen exposure to amplify these effects. Moreover, aero allergens are more difficult to effectively avoid compared to food allergens and might favor chronic stimulation of these putative pathways influencing NK cell phenotype and function.

NK cells sense stressed cells via NKG2D to detect a repertoire of surface expressed stress ligands **(*52, 53*)**. Oncologic transformation, microbial infections or environmental stressors promotes the expression of NKG2D ligands (*53-55*). Here we show that NKG2D interaction with its ligands promotes NK cells activation (Figure 5). Moreover, the above-described interactions allow NK cells to eliminate virally transfected cells and stimulate myeloid cells (*56, 57*). Thus, we can speculate that reduced expression of NKG2D would be reflected in reduced ability to combat pathogens. Several studies showed increased susceptibility to viral and bacterial infections in children with AD (*58, 59*). Recently, a study conducted by Medeleanu et al. suggested that poly-sensitized children are more susceptible to lower respiratory tract infection in early life in addition to an increased association to AD comorbidities (*60*). Thus, our data supports the idea that NKG2D and CD16 are downstream to allergen sensitization, which potentially affects antimicrobial response in AD.

One curious observation in our report is the positive relationship between the severity of AD and the proportion of NKG2Dlow NK cells in sensitized children. Our analysis as well as other studies showed an association between allergen sensitization and breach in the skin barrier (*61*). Thus, interpretation of these data suggests that reduced integrity of skin barriers combined with sensitization may contribute to the downregulation of NKG2D. This may relate to internalization of allergens or inflammatory stimulants through breaches in the epidermal barrier. Indeed, introduction of *Staphylococcus aureus* through the skin barrier at sites of eczema lesions would provide additional ligands for NKG2D on monocytes (*62*) that could subsequently stimulate NK cells and promote loss of NKG2D. A second interpretation of our data holds that accumulation of hyperinflammatory NKG2Dlow NK cells in certain patients could exacerbate epidermal integrity disruption resulting in increased transdermal exposure to allergens. This hypothesis is supported by the progressive decline in NKG2D expression by NK cells over subsequent clinical visits.

TNF-α impairs skin barrier function (*63, 64*) and could therefore contribute to worsening of AD as well as increased translocation of potential allergens. Of note, elevated TNF-α expression by NK cells is also a characteristic of severe asthmatic patients (*26*). The increased capacity for IFN-γ and TNF-α secretion by NKG2Dlow NK cells present in sensitized patients in response to MICA may also contribute to the pathophysiology of atopic march. Indeed, an IFN-γ signature in the context of elevated type 2 cytokines is a characteristic of severe asthma (*26, 27, 65, 66*). In addition, activated NK cells can produce cytokines including TNF-α *(16)* and GM-CSF *(17)* that may aggravate allergic disease mechanisms (*18-21*). Indeed, disease pathogenesis may depend on a feedback loop between skin barrier disruptions (*22*) affecting NK-cell functionality that in turn exacerbates skin barrier disruption to amplify allergen exposure (*67*).

Reduced NKG2D- and CD16-mediated cytolytic functionality of NK cells in allergen-sensitized patients with AD may also undermine key immunoregulatory mechanisms that limit development of allergy, allergic rhinitis, and asthma. NK cell killing of dendritic cells (*10*), CD4 T cells (*68*), and eosinophils (*69*) is likely to limit development of allergic disease. Therefore, the diminished capacity of NKG2Dlow NK cells in allergen sensitized patients with AD may prime the development of more severe allergic diseases and asthma. Of note, the development of asthma is typically established early in life and has been shown to be associated with increased microbial infections (*70, 71*).

Our data show that a progressive accumulation of hyperinflammatory NKG2Dlow NK cells defines a novel endotype of severe eczema associated with allergen sensitivity. These observations reveal a pathophysiological mechanism of NK cell dysregulation in AD that potentially contributes to atopic march. Therefore, the current paradigm of NK-cell deficiency in AD and associated translational focus on therapeutic restoration of NK cell activity may need to be re-evaluated in light of possible pathogenic contributions of NK-cell cytokine production to skin barrier disruption and development of allergic co-morbidities in this patient population.

## Materials and Methods

### Study Subjects

MPAACH is a prospective longitudinal early life cohort of children with AD that has been previously described (*18*). Recruitment began in December of 2016 and is ongoing. Briefly, children in the greater Cincinnati area (southwest Ohio and northern Kentucky) were identified through public advertising or the medical chart at Cincinnati Children’s Hospital Medical Center. Eligible children were aged ≤2 years upon enrollment, had a gestation of ≥36 weeks and either a diagnosis of AD (based on the Hanifin and Rajka Criteria for Atopic Dermatitis (*20*)) or the parent(s)/legal authorized representative indicated a positive response to each of the 3 questions from the Children’s Eczema Questionnaire (CEQ) (*21*). Each year following enrollment, MPAACH children undergo annual visits during which clinical data is obtained through questionnaires, extensive biospecimens are collected, and the children undergo skin prick testing (described below). At age 7, pulmonary function testing is performed. In a subset of experiments, children without AD were also recruited to serve as controls for baseline NK-cell functional responses. A set of adult peripheral blood mononuclear cells (PBMC) were used as a phenotyping control. This study was approved by the Institutional Review Board at Cincinnati Children’s Hospital Medical Center, and all subjects signed informed consent/parental permission prior to participation.

### Clinical Outcome Definitions

At each annual visit, MPAACH children undergo assessment of the severity of their AD using Scoring for Atopic Dermatitis (SCORAD) (*22*), and the quality of their skin barrier function is assessed by quantification of transepidermal water loss (TEWL) at both lesional and never-lesional sites at each visit using the DermaLab TEWL probe (Cortex Technolog, Hadsund, Denmark) as previously described (*18*). All children undergo annual skin prick testing (SPT) to 11 aeroallergens and 6 foods. The aeroallergen panel includes mold mix 1 (*Alternaria tenuis, Hormodendrum Cladosporioides, Helminthosporium interseminatum, Aspergillus fumigatus, Aspergillus Niger, Penicillium notatum*), mold mix 2 (*Rhizopus nigricans, Pullularia pullulans, Fusarium vasinfectum, Mucor racemosus*), ragweed (*Ambrosia trifida, Ambrosia artemisiifolia*), grass (Kentucky Bluegrass, Orchard, Redtop, Timothy, Sweet Vernon), tree mix 1 (Pecan, Maple BHR, Oak RVW, American Sycamore, Black Willow), tree mix 2 (White Ash, Birch Mix (Red&White), Black Walnut, Common Cottonwood, & American Elm), weeds (Kochia, English Plantain, Lamb’s Quarters, Marshelder/Poverty (BPT Mix) Common Cocklebur, & Careless/Pigweed (CR Mix)), mite mix (*Dermatophagoides pteronyssinus* and *D. farina*), cockroach (*Periplaneta americana* and *Blattella germanica*), cat and dog (HollisterStier, Spokane, WA). The food allergen panel includes cow’s milk, egg white, egg yolk, soy, wheat, tree nuts (almond, brazil nut, cashew, hazel nut, pecan, pistachio, and walnut), and peanut (Greer, Lenoir, NC). A positive SPT is defined as a wheal measurement 3mm or larger than the diluent control. Co-sensitization was defined as a positive SPT to at least one or more aeroallergens and one or more foods. Four longitudinal individual-level sensitization phenotypes were defined based on sensitization data from V1, V2 and V3 in the same child: *persistently non-sensitized* (non-sensitized at V1, V2 and V3); *transient* (sensitized at V1 or V2, but not at V3); *acquired* (non-sensitized at V1 and/or V2 but sensitized at V3); and *persistently sensitized* (sensitized at V1 and/or V2 and V3).

### Blood Sample Collection and Processing

Whole blood samples collected from MPAACH children were used to isolate human PBMC using Ficoll®Paque Plus (Sigma-Aldrich, Inc. St. Louis). The concentration of white blood cells in PBMC was measured using DxH 520 Hematology Analyzer (Beckman Coulter, CA). PBMC were resuspended in 1 mL of RPMI 1640 medium with 20 % FBS and stored at 4°C for immediate analysis or at -80°C for future studies.

### Antibodies

The following conjugated antibodies (Abs) were used in the described studies: Brilliant violet (BV) 421 anti-CD56 (5.1H11), BV711 anti-CD16 (3G8), Alexa-Fluor (AF) 647 anti-NKp44 (9E2), PE anti-Tim-3 (P44-8), PE/Cy7 anti-CCR4, APC anti-IFN-γ (4S.B3), PE anti-TNF-α (MAb11), PE/Cy7 anti-CD107a (H4A3) and AF488 anti-CD45RO were purchased from Biolegend (San Diego, CA); Brilliant Ultra violet (BUV) 395 anti-CD3 (SK7), V500 anti-CD3 (SP34-2), V500 anti-CD19, V500 anti-CD14, BUV737 anti-CD69 (FN50), and BUV395 anti-CD62L (SK11) were purchased from Becton and Dickinson (BD, San Diego, CA); eFluor710 anti-NKG2D (1D11) from Invitrogen (Carlsbad, CA); Abs were used at doses titrated in our lab.

### Sample analysis and internal controls

All phenotypic analyses were performed on fresh samples directly *ex vivo* to avoid complications of freeze/thaw. To control for batch effects across the cohort, we: (1) used the same antibody clones (as mentioned above), flow cytometer (LSR Fortessa instrument, BD Biosciences), staining protocol and reagents: Brilliant Violet staining buffer (BD Biosciences), running buffer (MACSQuant®, Miltenyi Biotec, Bergisch Gladbach, Germany) and fixation buffer (eBioscience™ IC Fixation Buffer, Invitrogen) through the entire study, (2) samples analysis was performed in weekly batches, and (3) included an internal control of healthy adult PBMC cultured with IL-2 (200U/mL) in each experimental batch.

### NK-cell activation assays

Thawed PBMC were cultured at 2×106 cells per well in the presence of 200 U/mL of human IL-2 (Peprotech®, Rocky Hill, NJ) for 18 hours. After 18 hours, incubated cells were stimulated in the presence or absence of plate bound major histocompatibility complex class I chain-related A (MICA) (10 μg/mL) for 6 hours. Golgistop™, GolgiPlug™ (BD), and PE/Cy7 anti-CD107a antibody (1:200) were added for the final five hours of the incubation (according to manufacturer guidelines). At harvest, cells were stained first for surface markers and then for the expression of intracellular cytokines by utilizing eBioscience™ Intracellular Fixation & Permeabilization buffer set (Invitrogen). Flow cytometric data for phenotypic and functional analysis was obtained using an LSR Fortessa instrument and analyzed via FlowJo_v10.8.1 software (Treestar). Samples were selected for NKG2D functional assays from MPAACH children at V3 based on NKG2D-MFI expression as outlined in Figure 5C.

### NK-cell cytokine production assessment

Thawed PBMC were cultured at 2×106 cells per well in the presence of 200 U/mL of human IL-2 (Peprotech®, Rocky Hill, NJ) for 18 hours. After 18 hours of incubation, cells were stimulated in the presence of IL-12 and IL-18 (2.5 ng/mL and 10 ng/mL) or media for 36 hours. At the last five hours of incubation, Golgistop™ and GolgiPlug™ (BD) were added (according to manufacturer’s guidelines). Cells were harvested and stained first for surface markers and then for the expression of intracellular cytokines by utilizing eBioscience™ Intracellular Fixation & Permeabilization buffer set (Invitrogen). Flow cytometric data for all phenotypic and functional analysis were obtained using an LSR Fortessa instrument (BD Biosciences) and analyzed via FlowJo_v10.8.1 software (Treestar).

### Flow cytometry dimensionality reduction analysis

Dimensionality reduction analysis was performed by utilizing Flowjo packages (Down sample, Uniform Manifold Approximation and Projection (UMAP) (arXiv:1802.03426v3) and unsupervised clustering was performed via FlowSOM, we instructed FlowSOM to generate 8 clusters of NK cells (*23*) [https://www.flowjo.com/exchange/#/](https://www.flowjo.com/exchange/#/). After utilizing Down-sample on the gated NK cells, only samples that contained exatctly 10,000 NK cells were utilized, this ensured sufficient numbers of cells for generation of the concatenated file and for downstream analysis by UMAP and FlowSOM.

### NK-cell cytotoxicity assay

K562 target cells were labeled for 2 hours with 2 μCi of sodium chromate (Na251CrO4) per 2 × 104 target cells at 37 °C, 5% CO2 and then washed to remove excess 51Cr. Labeled target cells were added to 96-well round-bottom plates (1 × 104 cells/well). Isolated NK cells were added to the plates with E:T ratios ranging between 16:1 and 1:1. The amount of 51Cr released, which corresponds to target cell death, was measured by a gamma scintillation counter. The percent cytotoxicity against target cells was calculated as: ((experimental lysis – spontaneous lysis)/ (maximal lysis – spontaneous lysis)) ×100. To determine maximal lysis, 51Cr-labeled target cells were treated with 3% Triton X for 4 hours. To determine spontaneous release, target cells incubated without effector cells were used.

### Statistical analyses

Statistical analyses were performed in GraphPad Prism 9.5 (GraphPad software, CA) and R version 4.1.0 (R core team 2021). Prior to the analysis, Shapiro-Wilk tests were performed to examine the distribution of the data for each continuous variable. Descriptive statistics are presented in median (25th-75th percentile) or frequency (%) due to the non-normal distributions. Comparisons between the full cohort and the subset were conducted using Mann-Whitney U and Fisher’s exact tests for continuous and categorical variables, respectively. The frequency of each NK-cell subset was compared using Kruskal-Wallis test across visits. Spearman correlation coefficients were used to measure the strength of association between NKG2D and SCORAD and TEWL. For a subset of children for whom individual-level longitudinal NK cell data was available (n=10), analysis of the proportions of NKG2DlowCD56dim and NKG2DlowCD56bright NK cells expression was analyzed across multiple visits within the same child to determine their trajectories. The longitudinal analysis of NKG2D measurements on CD56dim and CD56bright NK cells was performed for each individual child within each sensitization phenotype. P-values < 0.1 are reported as significant in the manuscript. Details of statistical testing are indicated in each figure legend.

## Data Availability

All data produced in the present work are contained in the manuscript and raw data will be made available upon publication or reasonable request to the authors

## SUPPLEMENTAL MATERIAL

![Supplemental Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F8.medium.gif)

[Supplemental Figure 1.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F8)

Supplemental Figure 1. Age-dependent changes in blood NK cells subset composition in children in the MPAACH cohort.
**(A)** Representative flow cytometry analysis with NK cells sequentially gated as singlet, live lymphocyte, and lineage negative (CD3neg CD19neg CD14neg) events prior to discrimination of CD16+CD56neg, CD16+CD56dim, or CD16neg CD56bright NK cells. (**B**) Median values of the frequency of each NK-cell subsets at each of three annual study visits, V1 (light grey circles, n=44), V2 (grey circles, n=39), V3 (dark grey circles, n=41). (**C**) Strategies to gate NKG2Dlow and NKG2Dhigh NK cells in non-sensitized (black) and sensitized (red) MPAACH relative to IL-2-treated positive control NK cells (dark gray) and FMO control (light gray).

![Supplemental Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F9.medium.gif)

[Supplemental Figure 2.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F9)

Supplemental Figure 2. Characterization of Uniform Manifold Approximation and Projection clusters.
Analysis of surface marker expression on the various clusters of NK cells presented in Figure 1E.

![Supplemental Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F10.medium.gif)

[Supplemental Figure 3.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F10)

Supplemental Figure 3. Allergen sensitization dependent changes in CD56dim NK cell receptor expression across three annual visits.
The expression of various surface markers on non-sensitized (closed black circles, V1, n=23; V2, n=16; V3, n=18) and sensitized (open red circles, V1, n=21; V2, n=23; V3, n=23) (**A**) Representative histograms of mean fluorescence intensity showing reduced CD16 expression on NK cells from sensitized patients at V3 and comparison between medians of CD16 per visit. (**B-D**) Proportions of CD56dim NK cells in each group of patients expressing CCR4, CD62L, CD69, and Tim-3 at V1-V3. Significant differences between group medians were determined by Mann-Whitney U test, p<0.1 values are reported.

![Supplemental Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F11.medium.gif)

[Supplemental Figure 4.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F11)

Supplemental Figure 4. Allergen sensitization dependent changes in CD56bright NK cell receptor expression across three annual visits.
The expression of various surface markers on non-sensitized (closed black circles, V1, n=23; V2, n=16; V3, n=18) and sensitized (open red circles, V1, n=21; V2, n=23; V3, n=23) (**A-C**) Proportions of CD56bright NK cells in each group of patients expressing CCR4, CD62L, CD69, and Tim-3 at V1-V3. Significant differences between group medians were determined by Mann-Whitney U test. No significant difference between medians was observed.

![Supplemental Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F12.medium.gif)

[Supplemental Figure 5.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F12)

Supplemental Figure 5. Allergen sensitization dependent changes in CD56neg NK cell receptor expression across three annual visits.
The expression of various surface markers on non-sensitized (closed black circles, V1, n=23; V2, n=16; V3, n=18) and sensitized (open red circles, V1, n=21; V2, n=23; V3, n=23) (**A-C**) Proportions of CD56neg NK cells in each group of patients expressing CD16, CCR4, CD62L, CD69, and Tim-3 at V1-V3. Significant differences between group medians were determined by Mann-Whitney U test. No significant difference between medians was observed.

![Supplemental Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F13.medium.gif)

[Supplemental Figure 6.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F13)

Supplemental Figure 6. NK cells from participants with reduced cytotoxicity show reduced NKG2D expression.
NK cells derived from V1 and V2 for participants evaluated in Figure 4A+B without AD (closed gray circles, n=4) or MPAACH children with allergen sensitivity (open red circles, n=5) or non-sensitized (closed black circles, n=5) or were evaluated for MFI of NKG2D expression on (**A**) CD56dim and (**B**) CD56bright NK cells. Significant differences between group medians were determined by one-tailed Mann-Whitney U test, p<0.1 values are reported.

![Supplemental Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F14.medium.gif)

[Supplemental Figure 7.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F14)

Supplemental Figure 7. CD56dim and CD56bright NK cells in allergen-sensitized patients show increased expression of IFN-γ.
Following culture of V1 and V2 PBMC in media (unstimulated) or IL-12 and IL-18, NK cells derived from non-AD (closed gray circles, n = 7), non-sensitized (closed black circles, n=15) and allergen sensitized (open red circles, n=26) children were evaluated for cytokine expression by flow cytometry. Expression of IFN-γ on (**A**) CD56dim (**B**) CD56bright NK cells. Statistical differences between groups were determined by Mann-Whitney U test. No significant difference between medians was observed.

![Supplemental Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F15.medium.gif)

[Supplemental Figure 8.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F15)

Supplemental Figure 8. NKG2D expression on NK cells examined in cytokine stimulation assays.
NK cells derived from Figure 6 participants (V1 + V2) including children without AD (closed gray circles, n=6) or MPAACH children without (closed black circles, n=15) or with (open red circles, n=26) allergen sensitivity, were evaluated for MFI of NKG2D expression on (**A**) CD56dim and (**B**) CD56bright NK cells. Significant differences between group medians were determined by one-tailed Mann-Whitney U test. No significant difference between medians was observed.

![Supplemental Figure 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/05/2023.06.02.23290884/F16.medium.gif)

[Supplemental Figure 9.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/F16)

Supplemental Figure 9. Association between changes in NKG2D expression over visits and skin barrier function and AD severity.
**(A)** Association between L-TEWL and delta NKG2Dlow CD56dim and **(B)** delta NKG2Dlow CD56bright NK cells in children with longitudinal analysis of proportions of NK cells (closed blue circles, n = 10), **(C-D)** stratified by sensitization phenotype (non-sensitized children: closed black circles, n = 4; sensitized children: open red circles, n = 6). **(E)** Association between SCORAD and delta NKG2Dlow CD56dim and **(F)** delta NKG2Dlow CD56bright NK cells, **(G-H)** stratified by sensitization phenotype. Association was assessed using Spearman correlation test. Changes (delta) in NKG2Dlow CD56dim and NKG2Dlow CD56dim NK cells were calculated as V3 – V1 or V2 NK cells expression.

View this table:
[Supplemental Table 1.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/T2)

Supplemental Table 1. 
Sample size and demographics in children by visit1

View this table:
[Supplemental Table 2.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/T3)

Supplemental Table 2. 
Complete blood count and cytokine level by sensitization1

View this table:
[Supplemental Table 3.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/T4)

Supplemental Table 3. 
Characteristics of the MPAACH cohort and the subset used in UMAP at visit 11

View this table:
[Supplemental Table 4.](http://medrxiv.org/content/early/2023/06/05/2023.06.02.23290884/T5)

Supplemental Table 4. 
Demographics by sensitization phenotype1

## Acknowledgements

We thank the participating families and staff on the MPAACH study. Funding support provided by 5U19AI70235 (GKKH, JBM, LJM); American Heart Association (DEO); AADCRC (SA); R01AI148080, R01AR073228 (SNW), 5T32GM063483 (SBD).

## Footnotes

*   E-mail: Stephen.Waggoner{at}cchmc.org or Gurjit.Hershey{at}cchmc.org

*   Received June 2, 2023.
*   Revision received June 2, 2023.
*   Accepted June 5, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1. S. Weidinger,  L. A. Beck,  T. Bieber,  K. Kabashima,  A. D. Irvine, Atopic dermatitis. Nat Rev Dis Primers 4, 1 (2018).
    
    

2.  2. S. C. Dharmage et al., Atopic dermatitis and the atopic march revisited. Allergy 69, 17–27 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/all.12268&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24117677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

3.  3. J. M. Spergel, From atopic dermatitis to asthma: the atopic march. Ann Allergy Asthma Immunol 105, 99–106; quiz 107-109, 117 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.anai.2009.10.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20674819&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285832800002&link_type=ISI) 

4.  4. S. Romagnani, The role of lymphocytes in allergic disease. J Allergy Clin Immunol 105, 399–408 (2000).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1067/mai.2000.104575&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10719286&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088946700001&link_type=ISI) 

5.  5. M. R. Mack et al., Blood natural killer cell deficiency reveals an immunotherapy strategy for atopic dermatitis. Sci Transl Med 12, (2020).
    
    

6.  6. R. S. Lever,  M. J. Lesko,  R. M. MacKie,  D. M. Parrott, Natural killer cell activity in atopic dermatitis: a sequential study. Clin Allergy 15, 479–486 (1985).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3863720&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

7.  7. J. R. Jensen,  T. T. Sand,  A. S. Jorgensen,  K. Thestrup-Pedersen, Modulation of natural killer cell activity in patients with atopic dermatitis. J Invest Dermatol 82, 30–34 (1984).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1523-1747.ep12259055&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6418828&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1984RX67900009&link_type=ISI) 

8.  8. N. T. Kusaimi,  J. J. Trentin, Natural cell-mediated cytotoxic activity in the peripheral blood of patients with atopic dermatitis. Arch Dermatol 118, 568–571 (1982).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archderm.1982.01650200036012&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7103526&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

9.  9. J. Bi et al., NK Cells Alleviate Lung Inflammation by Negatively Regulating Group 2 Innate Lymphoid Cells. J Immunol 198, 3336–3344 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE5OC84LzMzMzYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

10. 10. K. Walwyn-Brown et al., Human NK Cells Lyse Th2-Polarizing Dendritic Cells via NKp30 and DNAM-1. J Immunol 201, 2028–2041 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjIwMS83LzIwMjgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

11. 11. A. Awad et al., Natural killer cells induce eosinophil activation and apoptosis. PLoS One 9, e94492 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0094492&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24727794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

12. 12. O. Haworth,  M. Cernadas,  B. D. Levy, NK cells are effectors for resolvin E1 in the timely resolution of allergic airway inflammation. J Immunol 186, 6129–6135 (2011).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4Ni8xMS82MTI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

13. 13. L. Mobus et al., Elevated NK-cell transcriptional signature and dysbalance of resting and activated NK cells in atopic dermatitis. J Allergy Clin Immunol 147, 1959–1965 e1952 (2021).
    
    

14. 14. L. Mobus et al., Atopic dermatitis displays stable and dynamic skin transcriptome signatures. J Allergy Clin Immunol 147, 213–223 (2021).
    
    

15. 15. F. Scholz,  S. Naik,  F. S. Sutterwala,  D. H. Kaplan, Langerhans Cells Suppress CD49a+ NK Cell-Mediated Skin Inflammation. J Immunol 195, 2335–2342 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE5NS81LzIzMzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

16. 16. M. G. Sherenian et al., Sensitization to peanut, egg or pets is associated with skin barrier dysfunction in children with atopic dermatitis. Clin Exp Allergy 51, 666–673 (2021).
    
    

17. 17. T. Gonzalez et al., Biofilm propensity of Staphylococcus aureus skin isolates is associated with increased atopic dermatitis severity and barrier dysfunction in the MPAACH pediatric cohort. Allergy 76, 302–313 (2021).
    
    

18. 18. J. M. Biagini Myers et al., Events in Normal Skin Promote Early-Life Atopic Dermatitis-The MPAACH Cohort. J Allergy Clin Immunol Pract 8, 2285–2293 e2286 (2020).
    
    

19. 19. J. W. Kroner et al., Atopic dermatitis independently increases sensitization above parental atopy: The MPAACH study. J Allergy Clin Immunol 145, 1464–1466 (2020).
    
    

20. 20. J. M. Hanifin,  G. Rajka, Diagnostic features of atopic dermatitis. Acta Derm Venereol 92 (uppl), 44–47 (1980).
    
    

21. 21. L. B. von Kobyletzki et al., Validation of a parental questionnaire to identify atopic dermatitis in a population-based sample of children up to 2 years of age. Dermatology 226, 222–226 (2013).
    
    

22. 22.Severity scoring of atopic dermatitis: the SCORAD index. Consensus Report of the European Task Force on Atopic Dermatitis. Dermatology 186, 23–31 (1993).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000247298&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8435513&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993KK23100006&link_type=ISI) 

23. 23. S. Van Gassen et al., FlowSOM: Using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A 87, 636–645 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cyto.a.22625&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25573116&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

24. 24. J. Devulder et al., Aberrant anti-viral response of natural killer cells in severe asthma. Eur Respir J 55, (2020).
    
    

25. 25. C. A. Anderson et al., A degradatory fate for CCR4 suggests a primary role in Th2 inflammation. J Leukoc Biol 107, 455–466 (2020).
    
    

26. 26. M. G. Duvall et al., Natural killer cell-mediated inflammation resolution is disabled in severe asthma. Sci Immunol 2, (2017).
    
    

27. 27. N. Farhadi et al., Natural killer cell NKG2D and granzyme B are critical for allergic pulmonary inflammation. J Allergy Clin Immunol 133, 827–835 e823 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2013.09.048&link_type=DOI) 

28. 28. T. Carbone et al., CD56highCD16-CD62L-NK cells accumulate in allergic contact dermatitis and contribute to the expression of allergic responses. J Immunol 184, 1102–1110 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE4NC8yLzExMDIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

29. 29. D. H. Raulet, Roles of the NKG2D immunoreceptor and its ligands. Nat Rev Immunol 3, 781–790 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri1199&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14523385&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185970700011&link_type=ISI) 

30. 30. M. Migueres et al., Types of sensitization to aeroallergens: definitions, prevalences and impact on the diagnosis and treatment of allergic respiratory disease. Clin Transl Allergy 4, 16 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/2045-7022-4-16&link_type=DOI) 

31. 31. V. Siroux et al., The asthma-rhinitis multimorbidity is associated with IgE polysensitization in adolescents and adults. Allergy 73, 1447–1458 (2018).
    
    

32. 32. G. Moscato et al., Importance of Polysensitization in Allergic Rhinitis and Asthma. Journal of Allergy and Clinical Immunology 117, S297 (2006).
    
    

33. 33. D. S. Bae,  Y. K. Hwang,  J. K. Lee, Importance of NKG2D-NKG2D ligands interaction for cytolytic activity of natural killer cell. Cell Immunol 276, 122–127 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cellimm.2012.04.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22613008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

34. 34. S. Ushio et al., Cloning of the cDNA for human IFN-gamma-inducing factor, expression in Escherichia coli, and studies on the biologic activities of the protein. J Immunol 156, 4274–4279 (1996).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE1Ni8xMS80Mjc0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

35. 35. J. H. Lee,  D. H. Cho,  H. J. Park, IL-18 and Cutaneous Inflammatory Diseases. Int J Mol Sci 16, 29357–29369 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms161226172&link_type=DOI) 

36. 36. K. L. Hon et al., Serum concentration of IL-18 correlates with disease extent in young children with atopic dermatitis. Pediatr Dermatol 21, 619–622 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.0736-8046.2004.21600.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15575842&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

37. 37. Q. Hamid et al., In vivo expression of IL-12 and IL-13 in atopic dermatitis. J Allergy Clin Immunol 98, 225–231 (1996).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0091-6749(96)70246-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8765838&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VA69200027&link_type=ISI) 

38. 38. I. Meyts et al., IL-12 contributes to allergen-induced airway inflammation in experimental asthma. J Immunol 177, 6460–6470 (2006).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE3Ny85LzY0NjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

39. 39. B. Verhaeghe et al., Up-regulation of IL-18 in allergic rhinitis. Allergy 57, 825–830 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1398-9995.2002.23539.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12169180&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177359800013&link_type=ISI) 

40. 40. J. S. Wild et al., IFN-gamma-inducing factor (IL-18) increases allergic sensitization, serum IgE, Th2 cytokines, and airway eosinophilia in a mouse model of allergic asthma. J Immunol 164, 2701–2710 (2000).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTA6IjE2NC81LzI3MDEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

41. 41. P. M. Brunner,  E. Guttman-Yassky,  D. Y. Leung, The immunology of atopic dermatitis and its reversibility with broad-spectrum and targeted therapies. J Allergy Clin Immunol 139, S65–S76 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2017.01.011&link_type=DOI) 

42. 42. J. Brady et al., The interactions of multiple cytokines control NK cell maturation. J Immunol 185, 6679–6688 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4NS8xMS82Njc5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

43. 43. C. Zhang et al., Interleukin-12 improves cytotoxicity of natural killer cells via upregulated expression of NKG2D. Hum Immunol 69, 490–500 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.humimm.2008.06.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18619507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

44. 44. S. Romagnani, The Th1/Th2 paradigm. Immunol Today 18, 263–266 (1997).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0167-5699(97)01070-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9190109&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XD23800004&link_type=ISI) 

45. 45. T. A. Wynn, Type 2 cytokines: mechanisms and therapeutic strategies. Nat Rev Immunol 15, 271–282 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri3831&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25882242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

46. 46. X. Zhou et al., Cytometric analysis reveals an association between allergen-responsive natural killer cells and human peanut allergy. J Clin Invest 132, (2022).
    
    

47. 47. G. Matusali et al., Plasma levels of soluble MICA and ULBP2 are increased in children allergic to dust mites. J Allergy Clin Immunol 130, 1003–1005 (2012).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22921872&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

48. 48. G. Matusali et al., Soluble ligands for the NKG2D receptor are released during HIV-1 infection and impair NKG2D expression and cytotoxicity of NK cells. FASEB J 27, 2440–2450 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1096/fj.12-223057&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23395909&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319667600031&link_type=ISI) 

49. 49. M. Jinushi et al., Impairment of natural killer cell and dendritic cell functions by the soluble form of MHC class I-related chain A in advanced human hepatocellular carcinomas. J Hepatol 43, 1013–1020 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jhep.2005.05.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16168521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234059700015&link_type=ISI) 

50. 50. R. Romee et al., NK cell CD16 surface expression and function is regulated by a disintegrin and metalloprotease-17 (ADAM17). Blood 121, 3599–3608 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czoxMToiMTIxLzE4LzM1OTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

51. 51. I. Waldhauer et al., Tumor-associated MICA is shed by ADAM proteases. Cancer Res 68, 6368–6376 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjEwOiI2OC8xNS82MzY4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

52. 52. S. Bauer et al., Activation of NK cells and T cells by NKG2D, a receptor for stress-inducible MICA. Science 285, 727–729 (1999).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyODUvNTQyOC83MjciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. 53. D. Cosman et al., ULBPs, novel MHC class I-related molecules, bind to CMV glycoprotein UL16 and stimulate NK cytotoxicity through the NKG2D receptor. Immunity 14, 123–133 (2001).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1074-7613(01)00095-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11239445&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000167092500003&link_type=ISI) 

54. 54. A. Soriani et al., Reactive oxygen species- and DNA damage response-dependent NK cell activating ligand upregulation occurs at transcriptional levels and requires the transcriptional factor E2F1. J Immunol 193, 950–960 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6OToiMTkzLzIvOTUwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

55. 55. T. Ebihara et al., Induction of NKG2D ligands on human dendritic cells by TLR ligand stimulation and RNA virus infection. Int Immunol 19, 1145–1155 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/intimm/dxm073&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17878262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000250681700001&link_type=ISI) 

56. 56. T. T. Greene et al., A Herpesviral induction of RAE-1 NKG2D ligand expression occurs through release of HDAC mediated repression. Elife 5, (2016).
    
    

57. 57. P. Eissmann et al., Multiple mechanisms downstream of TLR-4 stimulation allow expression of NKG2D ligands to facilitate macrophage/NK cell crosstalk. J Immunol 184, 6901–6909 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6MTE6IjE4NC8xMi82OTAxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDYvMDUvMjAyMy4wNi4wMi4yMzI5MDg4NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

58. 58. P. Y. Ong,  D. Y. Leung, Bacterial and Viral Infections in Atopic Dermatitis: a Comprehensive Review. Clin Rev Allergy Immunol 51, 329–337 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12016-016-8548-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

59. 59. M. Bohme,  E. Lannero,  M. Wickman,  S. L. Nordvall,  C. F. Wahlgren, Atopic dermatitis and concomitant disease patterns in children up to two years of age. Acta Derm Venereol 82, 98–103 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00015550252948112&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12125961&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176014600006&link_type=ISI) 

60. 60. M. Medeleanu et al., Moderate-to-severe lower respiratory tract infection in early life is associated with increased risk of polysensitization and atopic dermatitis: Findings from the CHILD Study. Journal of Allergy and Clinical Immunology: Global 1, 73–79 (2022).
    
    

61. 61. F. De Marchi et al., Correlation of skin barrier impairment in atopic dermatitis with aeroallergen sensitization. Allergy Asthma Proc 36, e127–133 (2015).
    
    

62. 62. M. Mellergaard et al., Staphylococcus aureus induces cell-surface expression of immune stimulatory NKG2D ligands on human monocytes. J Biol Chem 295, 11803–11821 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyOTUvMzMvMTE4MDMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

63. 63. M. O. Danso et al., TNF-alpha and Th2 cytokines induce atopic dermatitis-like features on epidermal differentiation proteins and stratum corneum lipids in human skin equivalents. J Invest Dermatol 134, 1941–1950 (2014).
    
    

64. 64. R. L. Riding,  J. E. Harris, The Role of Memory CD8(+) T Cells in Vitiligo. J Immunol 203, 11–19 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiamltbXVub2wiO3M6NToicmVzaWQiO3M6ODoiMjAzLzEvMTEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wNS8yMDIzLjA2LjAyLjIzMjkwODg0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

65. 65. J. A. Wisniewski et al., TH1 signatures are present in the lower airways of children with severe asthma, regardless of allergic status. J Allergy Clin Immunol 141, 2048–2060 e2013 (2018).
    
    

66. 66. J. W. Steinke et al., Bronchoalveolar lavage cytokine patterns in children with severe neutrophilic and paucigranulocytic asthma. J Allergy Clin Immunol 147, 686–693 e683 (2021).
    
    

67. 67. K. S. Lee et al., The prevention of TNF-alpha/IFN-gamma mixture-induced inflammation in human keratinocyte and atopic dermatitis-like skin lesions in Nc/Nga mice by mineral-balanced deep sea water. Biomed Pharmacother 97, 1331–1340 (2018).
    
    

68. 68. S. N. Waggoner,  M. Cornberg,  L. K. Selin,  R. M. Welsh, Natural killer cells act as rheostats modulating antiviral T cells. Nature 481, 394–398 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10624&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22101430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom) 

69. 69. S. Pesce et al., The Innate Immune Cross Talk between NK Cells and Eosinophils Is Regulated by the Interaction of Natural Cytotoxicity Receptors with Eosinophil Surface Ligands. Front Immunol 8, 510 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fimmu.2017.00510&link_type=DOI) 

70. 70. S. C. Dharmage,  J. L. Perret,  A. Custovic, Epidemiology of Asthma in Children and Adults. Front Pediatr 7, 246 (2019).
    
    

71. 71. A. Sharma et al., Associations between fungal and bacterial microbiota of airways and asthma endotypes. J Allergy Clin Immunol 144, 1214–1227 e1217 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2019.06.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31279011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F05%2F2023.06.02.23290884.atom)