Respiratory viral infections increase risk of asthma in infants and children. Infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus can cause severe lung inflammation and prolonged respiratory symptoms. We sought to determine whether SARS-CoV-2 infection modified pediatric incident asthma risk.
This retrospective cohort study examined children ages 1 to 16 within the Children’s Hospital of Philadelphia Care Network who received polymerase chain reaction (PCR) testing for SARS-CoV-2 between March 1, 2020 and February 28, 2021. Multivariable Cox regression models assessed the hazard ratio of new asthma diagnosis between SARS-CoV-2 PCR positive and SARS-CoV-2 PCR negative groups within an 18-month observation window. Models were adjusted for demographic characteristics, socioeconomic variables, and atopic comorbidities.
There were 27 423 subjects included in the study. In adjusted analyses, SARS-CoV-2 PCR positivity had no significant effect on the hazard of new asthma diagnosis (hazard ratio [HR]: 0.96; P = .79). Black race (HR: 1.49; P = .004), food allergies (HR: 1.26; P = .025), and allergic rhinitis (HR: 2.30; P < .001) significantly increased the hazard of new asthma diagnosis. Preterm birth (HR: 1.48; P = .005) and BMI (HR: 1.13; P < .001) significantly increased the hazard of new asthma diagnosis for children <5 years old.
SARS-CoV-2 PCR positivity was not associated with new asthma diagnosis in children within the observation period, although known risk factors for pediatric asthma were confirmed. This study informs the prognosis and care of children with a history of SARS-CoV-2 infection.
Respiratory viral infections, particularly rhinovirus and respiratory syncytial virus, are risk factors for the onset of pediatric asthma. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can cause severe lower respiratory tract inflammation, but little is known about if SARS-CoV-2 infection influences development of asthma.
By establishing that a history of SARS-CoV-2 polymerase chain reaction positivity is not associated with new asthma diagnosis, this study broadens our understanding of the long-term respiratory consequences of SARS-CoV-2 infection and informs the prognosis and care of millions of US children.
Respiratory viral infections during childhood have been identified as a potential risk factor for the development of pediatric asthma. Numerous longitudinal studies have demonstrated associations between acute wheezing illnesses in infancy secondary to infection with rhinovirus or respiratory syncytial virus (RSV) and subsequent progression to asthma in childhood and adolescence.1 ,2 Rhinovirus and RSV have been shown to induce type II inflammatory pathways that may lead to persistent airway inflammation, allergic sensitization, and remodeling seen in pediatric asthma.1 ,3 Smaller studies have isolated coronavirus strains from nasopharyngeal aspirates of infants and children during acute wheezing episodes or asthma exacerbations.4 –6 However, each of these studies identified aspirates containing human coronaviruses 229E and OC43, coronavirus strains with low pathogenicity and mostly responsible for benign upper respiratory tract infections.7 Thus, although respiratory viral infections have been associated with incident asthma, less is known about the impact of coronaviruses on the development of asthma in children.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in 2019 in Wuhan, China, is a highly pathogenic coronavirus strain with the potential for more severe lower respiratory symptoms including pneumonia or acute respiratory distress syndrome.7 As a severe inflammatory disease of the lower respiratory tract, SARS-CoV-2 infection may lead to prolonged respiratory symptoms even after resolution of infection.8 A small study from Iran found a higher incidence of chronic cough and asthma-like symptoms in children after hospitalization for SARS-CoV-2 infection, particularly in children with a personal history of atopy and family history of asthma.9 However, there are no studies to date that have investigated the incidence of asthma diagnosis after SARS-CoV-2 infection in a pediatric population. To address this knowledge gap, we used a large longitudinal pediatric cohort to ascertain whether children who tested positive for SARS-CoV-2 had altered rates of subsequent asthma diagnosis compared with those who tested negative for SARS-CoV-2 during the same interval.
Methods
Setting and Participants
Data were abstracted from the electronic health records (EHR) of children ages 1 to 16 who received care within the Children’s Hospital of Philadelphia (CHOP) Care Network from 2019 to 2022. We limited our analysis to subjects who obtained polymerase chain reaction (PCR) testing for SARS-CoV-2 between March 1, 2020 and February 28, 2021. Of note, SARS-CoV-2 PCR testing became available at CHOP in March 2020, when the first incident case was reported.10 The interval between March 1, 2020 and February 28, 2021 preceded both the introduction of routine rapid antigen testing and the introduction of vaccination for pediatric patients. All subjects were considered unvaccinated during the exposure window.
All subjects were required to have at least 1 ambulatory well child visit in the CHOP Care Network (CCN) in the year preceding the PCR test and at least 1 CCN visit at any time during an 18-month follow-up period after their first positive or last negative PCR test (if all test results were negative). A CCN follow-up visit could include primary care, urgent care, emergency department visits, telemedicine, or inpatient admission. Additionally, this 18-month follow-up period began 30 days after the corresponding PCR test. A 30-day lag time between the exposure and the start of the observation period was chosen to provide a reasonable separation between measurement of the exposure and the outcome. This excluded study subjects who may have received an asthma diagnosis on the same day or very close to the day of their SARS-CoV-2 PCR test.
Measures
Our primary definition of asthma diagnosis was the presence of at least 1 corresponding International Classification of Diseases (ICD) code for asthma (codes 493.×, J45.×) at the CCN visit, along with a prescription of an asthma-related medication. ICD codes related to “reactive airway disease” and “postviral wheeze” were excluded. An asthma-related medication included any of the following: β agonists, inhaled corticosteroids, combined inhaled corticosteroids and long-acting β agonists, leukotriene modulators, or biologics. For a sensitivity analysis, a secondary definition of asthma diagnosis was used that required a corresponding ICD code for asthma at 2 or more follow-up visits at least 6 months apart, as well as prescription of an asthma-related medication.
Covariates
Sex, age, race, and insurance status of study subjects were included for analysis. Although we acknowledge that race and ethnicity are social constructs, race was included as a variable in our analyses given prior research showing that Black and Latinx children have higher asthma prevalence and greater asthma morbidity, most likely because of a combination of factors, including structural racism, healthcare inequities, differential environmental exposures, and the effects of chronic, long-term stress.11 The Childhood Opportunity Index is a validated approach for estimating social determinants of health across various neighborhood-level indices,12 and this was used in our study to approximate socioeconomic status, housing quality, and the built environment. Gestational age, BMI, and atopic comorbidities (specifically atopic dermatitis, food allergies, and allergic rhinitis), were also included as potential confounders.
Statistical Analyses
The covariates of interest were summarized and compared between SARS-CoV-2 positive and SARS-CoV-2 negative groups using Pearson’s χ-square test for categorical variables and 2-sample t test for continuous variables (Table 1). Multivariable Cox regression models assessed the hazard of new asthma diagnosis comparing SARS-CoV-2 positive and SARS-CoV-2 negative groups. Models were adjusted for the covariates outlined in Table 1. Subgroup analyses were performed based on the age groups defined by the National Heart, Lung, and Blood Institute guidelines (ages 0–4, 5–11, and ≥12 years) given the phenotypic similarities within each group.13 The subgroup analyses used the primary definition of asthma for the outcome variable.
Variable . | SARS-CoV-2 PCR+, (n = 3147) . | SARS-CoV-2 PCR-, (n = 24 276) . | Total, (n = 27 423) . | P . |
---|---|---|---|---|
Sex | .88 | |||
Male, n (%) | 1634 (52) | 12 680 (52) | 14 314 (52) | |
Agea, n (%) | <.001 | |||
1–4 y | 885 (28.1) | 9351 (38.5) | 10 236 (37.3) | |
5–11 y | 1272 (40.4) | 9416 (38.8) | 10 688 (39.0) | |
≥12 y | 990 (31.5) | 5509 (22.7) | 6499 (23.7) | |
Racea, n (%) | <.001 | |||
White | 1755 (55.8) | 14 889 (61.3) | 16 644 (60.7) | |
Black | 780 (24.9) | 5048 (20.8) | 5828 (21.3) | |
Asian | 96 (3.05) | 796 (3.28) | 892 (3.25) | |
Multiple races | 119 (3.78) | 1012 (4.17) | 1131 (4.12) | |
Other | 397 (12.6) | 2531 (10.4) | 2928 (10.7) | |
Insurancea, n (%) | <.001 | |||
Non-Medicaid | 2170 (69.0) | 18 125 (74.7) | 20 295 (74.0) | |
Childhood Opportunity Index classification,a n (%) | <.001 | |||
Very low | 733 (23.4) | 4371 (18.2) | 5104 (18.6) | |
Low | 214 (6.8) | 1401 (5.82) | 1615 (5.89) | |
Moderate | 350 (11.2) | 2372 (9.86) | 2722 (9.93) | |
High | 625 (20.0) | 4763 (19.8) | 5388 (19.6) | |
Very high | 1206 (38.6) | 11 160 (46.4) | 12 366 (45.1) | |
Gestational agea, n (%) | <.001 | |||
Full term | 1796 (57.1) | 14 592 (60.1) | 16 388 (59.8) | |
Preterm | 290 (9.22) | 2584 (10.6) | 2874 (10.5) | |
Post-term | 214 (6.80) | 1788 (7.37) | 2002 (7.30) | |
Unknown | 847 (26.9) | 5312 (21.9) | 6159 (22.5) | |
Mean BMI at first visit (SD)a | 19.3 (4.8) | 18.4 (11.3) | 18.5 (10.8) | <.001 |
Atopic dermatitis, n (%) | 372 (11.8) | 3118 (12.8) | 3490 (12.7) | .11 |
Positive diagnosis | ||||
Food allergya, n (%) | 432 (13.7) | 3871 (16.0) | 4303 (15.7) | .001 |
Positive diagnosis | ||||
Allergic rhinitisa, n (%) | 1474 (46.8) | 10 496 (43.2) | 11 970 (43.6) | <.001 |
Positive diagnosis |
Variable . | SARS-CoV-2 PCR+, (n = 3147) . | SARS-CoV-2 PCR-, (n = 24 276) . | Total, (n = 27 423) . | P . |
---|---|---|---|---|
Sex | .88 | |||
Male, n (%) | 1634 (52) | 12 680 (52) | 14 314 (52) | |
Agea, n (%) | <.001 | |||
1–4 y | 885 (28.1) | 9351 (38.5) | 10 236 (37.3) | |
5–11 y | 1272 (40.4) | 9416 (38.8) | 10 688 (39.0) | |
≥12 y | 990 (31.5) | 5509 (22.7) | 6499 (23.7) | |
Racea, n (%) | <.001 | |||
White | 1755 (55.8) | 14 889 (61.3) | 16 644 (60.7) | |
Black | 780 (24.9) | 5048 (20.8) | 5828 (21.3) | |
Asian | 96 (3.05) | 796 (3.28) | 892 (3.25) | |
Multiple races | 119 (3.78) | 1012 (4.17) | 1131 (4.12) | |
Other | 397 (12.6) | 2531 (10.4) | 2928 (10.7) | |
Insurancea, n (%) | <.001 | |||
Non-Medicaid | 2170 (69.0) | 18 125 (74.7) | 20 295 (74.0) | |
Childhood Opportunity Index classification,a n (%) | <.001 | |||
Very low | 733 (23.4) | 4371 (18.2) | 5104 (18.6) | |
Low | 214 (6.8) | 1401 (5.82) | 1615 (5.89) | |
Moderate | 350 (11.2) | 2372 (9.86) | 2722 (9.93) | |
High | 625 (20.0) | 4763 (19.8) | 5388 (19.6) | |
Very high | 1206 (38.6) | 11 160 (46.4) | 12 366 (45.1) | |
Gestational agea, n (%) | <.001 | |||
Full term | 1796 (57.1) | 14 592 (60.1) | 16 388 (59.8) | |
Preterm | 290 (9.22) | 2584 (10.6) | 2874 (10.5) | |
Post-term | 214 (6.80) | 1788 (7.37) | 2002 (7.30) | |
Unknown | 847 (26.9) | 5312 (21.9) | 6159 (22.5) | |
Mean BMI at first visit (SD)a | 19.3 (4.8) | 18.4 (11.3) | 18.5 (10.8) | <.001 |
Atopic dermatitis, n (%) | 372 (11.8) | 3118 (12.8) | 3490 (12.7) | .11 |
Positive diagnosis | ||||
Food allergya, n (%) | 432 (13.7) | 3871 (16.0) | 4303 (15.7) | .001 |
Positive diagnosis | ||||
Allergic rhinitisa, n (%) | 1474 (46.8) | 10 496 (43.2) | 11 970 (43.6) | <.001 |
Positive diagnosis |
Significant test results given α = .05.
All analyses were completed using Stata 15.1 (Stata Corp, College Station, TX). A significance level of P < .05 was chosen a priori for all analyses. The CHOP Institutional Review Board determined this study was not human subject’s research.
Results
A description of the study sample (n = 27 423) appears in Table 1. Of the subjects that met inclusion criteria, 3147 subjects (11.5%) were in the SARS-CoV-2 positive group and 24 276 (88.5%) were in the SARS-CoV-2 negative group. SARS-CoV-2 positive individuals tended to be older, Black, Medicaid insured, among lower childhood opportunity index quintiles, and have a higher BMI (Table 1). Regarding atopic comorbidities, SARS-CoV-2 positive individuals were less likely to have a food allergy diagnosis but more likely to have an allergic rhinitis diagnosis than SARS-CoV-2 negative individuals (Table 1).
There were 573 total subjects who received an asthma diagnosis in the 18-month follow-up period. In our primary analysis of these subjects, 1.81% (n = 57) of SARS-CoV-2 positive subjects were subsequently diagnosed with asthma compared with 2.13% (n = 516) of SARS-CoV-2 negative subjects. Overall, SARS-CoV-2 positive subjects trended toward lower incidence of new asthma diagnosis at most time points over the 18-month follow-up period (Fig 1). In the regression analyses (Table 2), SARS-CoV-2 PCR positivity was not associated with the hazard of new asthma diagnosis in the 18-month follow-up period (hazard ratio [HR]: 0.96; 95% confidence interval [CI]: 0.73–1.27). Black race (HR: 1.49; 95% CI: 1.13–1.95), comorbid food allergy (HR: 1.26; 95% CI: 1.03–1.55), and allergic rhinitis (HR: 2.30; 95% CI: 1.93–2.74) were associated with a significantly higher hazard of new asthma diagnosis, whereas older age (ages 5–11: HR: 0.27; 95% CI: 0.22–0.34; ages 12+: HR: 0.16; 95% CI: 0.12–0.22) was associated with a significantly lower hazard of new asthma diagnosis (Table 2).
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 24 276 | (Ref) | (Ref) | (Ref) |
Positive | 3147 | 0.96 | 0.73–1.27 | .79 |
Sex | ||||
Male | 14 314 | (Ref) | (Ref) | (Ref) |
Female | 13 109 | 1.03 | 0.87–1.21 | .73 |
Age | ||||
1–4 y | 10 236 | (Ref) | (Ref) | (Ref) |
5–11 ya | 10 688 | 0.27 | 0.22–0.34 | <.001a |
≥12 ya | 6499 | 0.16 | 0.12–0.22 | <.001a |
Race | ||||
White | 16 644 | (Ref) | (Ref) | (Ref) |
Blacka | 5828 | 1.49 | 1.13–1.95 | .004a |
Asian | 892 | 1.20 | 0.75–1.91 | .45 |
Multiple racesa | 1131 | 1.46 | 1.01–2.12 | .044a |
Other | 2928 | 1.30 | 0.98–1.71 | .067 |
Insurance | ||||
Non-Medicaid | 20 295 | (Ref) | (Ref) | (Ref) |
Medicaid | 7128 | 1.00 | 0.81–1.23 | .995 |
Childhood Opportunity Index classification | ||||
Very low | 5104 | (Ref) | (Ref) | (Ref) |
Low | 1615 | 0.77 | 0.52–1.14 | .19 |
Moderate | 2722 | 0.82 | 0.58–1.17 | .28 |
High | 5388 | 1.03 | 0.75–1.41 | .86 |
Very high | 12 366 | 1.02 | 0.75–1.37 | .91 |
Gestational age | ||||
Full term | 16 388 | (Ref) | (Ref) | (Ref) |
Preterm | 2874 | 1.25 | 0.98–1.56 | .074 |
Post-term | 2002 | 0.90 | 0.64–1.26 | .54 |
Unknown | 6159 | 1.24 | 0.99–1.57 | .064 |
BMI | 27 280 | 1.00 | 0.99–1.01 | .935 |
Atopic dermatitis | ||||
No atopic dermatitis | 23 933 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 3490 | 1.03 | 0.82–1.29 | .82 |
Food allergies | ||||
No food allergies | 23 120 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 4303 | 1.26 | 1.03–1.55 | .025a |
Allergic rhinitis | ||||
No allergic rhinitis | 15 453 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 11 970 | 2.30 | 1.93–2.74 | <.001a |
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 24 276 | (Ref) | (Ref) | (Ref) |
Positive | 3147 | 0.96 | 0.73–1.27 | .79 |
Sex | ||||
Male | 14 314 | (Ref) | (Ref) | (Ref) |
Female | 13 109 | 1.03 | 0.87–1.21 | .73 |
Age | ||||
1–4 y | 10 236 | (Ref) | (Ref) | (Ref) |
5–11 ya | 10 688 | 0.27 | 0.22–0.34 | <.001a |
≥12 ya | 6499 | 0.16 | 0.12–0.22 | <.001a |
Race | ||||
White | 16 644 | (Ref) | (Ref) | (Ref) |
Blacka | 5828 | 1.49 | 1.13–1.95 | .004a |
Asian | 892 | 1.20 | 0.75–1.91 | .45 |
Multiple racesa | 1131 | 1.46 | 1.01–2.12 | .044a |
Other | 2928 | 1.30 | 0.98–1.71 | .067 |
Insurance | ||||
Non-Medicaid | 20 295 | (Ref) | (Ref) | (Ref) |
Medicaid | 7128 | 1.00 | 0.81–1.23 | .995 |
Childhood Opportunity Index classification | ||||
Very low | 5104 | (Ref) | (Ref) | (Ref) |
Low | 1615 | 0.77 | 0.52–1.14 | .19 |
Moderate | 2722 | 0.82 | 0.58–1.17 | .28 |
High | 5388 | 1.03 | 0.75–1.41 | .86 |
Very high | 12 366 | 1.02 | 0.75–1.37 | .91 |
Gestational age | ||||
Full term | 16 388 | (Ref) | (Ref) | (Ref) |
Preterm | 2874 | 1.25 | 0.98–1.56 | .074 |
Post-term | 2002 | 0.90 | 0.64–1.26 | .54 |
Unknown | 6159 | 1.24 | 0.99–1.57 | .064 |
BMI | 27 280 | 1.00 | 0.99–1.01 | .935 |
Atopic dermatitis | ||||
No atopic dermatitis | 23 933 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 3490 | 1.03 | 0.82–1.29 | .82 |
Food allergies | ||||
No food allergies | 23 120 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 4303 | 1.26 | 1.03–1.55 | .025a |
Allergic rhinitis | ||||
No allergic rhinitis | 15 453 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 11 970 | 2.30 | 1.93–2.74 | <.001a |
There were 27 423 total individuals, and 27 054 were used in the Cox regression analysis.
Significant test results given α = .05.
In our sensitivity analysis utilizing our secondary asthma definition that required an asthma ICD code at 2 or more follow-up visits, there were 418 total subjects who received and maintained an asthma diagnosis in the 18-month follow-up period. Of those, 1.12% (n = 35) of SARS-CoV-2 positive subjects were diagnosed with asthma compared with 1.59% (n = 383) of SARS-CoV-2 negative subjects. In the regression analyses (Table 3), SARS-CoV-2 PCR positivity was not associated with a new asthma diagnosis using this stricter definition. However, Black race, atopic comorbidities, and age were again associated with new asthma diagnosis when using the secondary asthma definition (Table 3).
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 24 143 | (Ref) | (Ref) | (Ref) |
Positive | 3125 | 0.80 | 0.56–1.13 | .20 |
Sex | ||||
Male | 14 234 | (Ref) | (Ref) | (Ref) |
Female | 13 034 | 1.05 | 0.86–1.27 | .63 |
Age | ||||
1–4 y | 10 162 | (Ref) | (Ref) | (Ref) |
5–11 ya | 10 639 | 0.21 | 0.16–0.27 | <.001a |
≥12 ya | 6467 | 0.077 | 0.048–0.12 | <.001a |
Race | ||||
White | 16 553 | (Ref) | (Ref) | (Ref) |
Blacka | 5790 | 1.71 | 1.24–2.34 | .001a |
Asian | 886 | 1.19 | 0.67–2.09 | .55 |
Multiple racesa | 1124 | 1.56 | 1.01–2.39 | .044a |
Othera | 2915 | 1.54 | 1.13–2.12 | .007a |
Insurance | ||||
Non-Medicaid | 20 180 | (Ref) | (Ref) | (Ref) |
Medicaid | 7088 | 1.10 | 0.86–1.41 | .43 |
Childhood Opportunity Index classification | ||||
Very low | 5069 | (Ref) | (Ref) | (Ref) |
Low | 1607 | 0.83 | 0.53–1.29 | .42 |
Moderate | 2705 | 0.80 | 0.53–1.21 | .29 |
High | 5355 | 1.10 | 0.77–1.59 | .59 |
Very high | 12 305 | 1.18 | 0.84–1.67 | .34 |
Gestational age | ||||
Full term | 16 297 | (Ref) | (Ref) | (Ref) |
Preterm | 2856 | 1.28 | 0.97–1.69 | .077 |
Post-term | 1996 | 1.02 | 0.71–1.47 | .91 |
Unknown | 6119 | 1.15 | 0.86–1.53 | .34 |
BMI | 27 127 | 1.00 | 0.99–1.01 | .83 |
Atopic dermatitis | ||||
No atopic dermatitis | 23 804 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 3464 | 0.94 | 0.72–1.23 | .67 |
Food allergies | ||||
No food allergies | 22 990 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 4278 | 1.34 | 1.09–1.72 | .008a |
Allergic rhinitis | ||||
No allergic rhinitis | 15 373 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 11 895 | 2.83 | 2.30–3.49 | <.001a |
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 24 143 | (Ref) | (Ref) | (Ref) |
Positive | 3125 | 0.80 | 0.56–1.13 | .20 |
Sex | ||||
Male | 14 234 | (Ref) | (Ref) | (Ref) |
Female | 13 034 | 1.05 | 0.86–1.27 | .63 |
Age | ||||
1–4 y | 10 162 | (Ref) | (Ref) | (Ref) |
5–11 ya | 10 639 | 0.21 | 0.16–0.27 | <.001a |
≥12 ya | 6467 | 0.077 | 0.048–0.12 | <.001a |
Race | ||||
White | 16 553 | (Ref) | (Ref) | (Ref) |
Blacka | 5790 | 1.71 | 1.24–2.34 | .001a |
Asian | 886 | 1.19 | 0.67–2.09 | .55 |
Multiple racesa | 1124 | 1.56 | 1.01–2.39 | .044a |
Othera | 2915 | 1.54 | 1.13–2.12 | .007a |
Insurance | ||||
Non-Medicaid | 20 180 | (Ref) | (Ref) | (Ref) |
Medicaid | 7088 | 1.10 | 0.86–1.41 | .43 |
Childhood Opportunity Index classification | ||||
Very low | 5069 | (Ref) | (Ref) | (Ref) |
Low | 1607 | 0.83 | 0.53–1.29 | .42 |
Moderate | 2705 | 0.80 | 0.53–1.21 | .29 |
High | 5355 | 1.10 | 0.77–1.59 | .59 |
Very high | 12 305 | 1.18 | 0.84–1.67 | .34 |
Gestational age | ||||
Full term | 16 297 | (Ref) | (Ref) | (Ref) |
Preterm | 2856 | 1.28 | 0.97–1.69 | .077 |
Post-term | 1996 | 1.02 | 0.71–1.47 | .91 |
Unknown | 6119 | 1.15 | 0.86–1.53 | .34 |
BMI | 27 127 | 1.00 | 0.99–1.01 | .83 |
Atopic dermatitis | ||||
No atopic dermatitis | 23 804 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 3464 | 0.94 | 0.72–1.23 | .67 |
Food allergies | ||||
No food allergies | 22 990 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 4278 | 1.34 | 1.09–1.72 | .008a |
Allergic rhinitis | ||||
No allergic rhinitis | 15 373 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 11 895 | 2.83 | 2.30–3.49 | <.001a |
There were 27 268 total individuals, and 26 899 were used in the Cox regression analysis.
Significant test results given α = .05.
For each subgroup analysis by age, there was no significant relationship between SARS-CoV-2 positivity and new asthma diagnosis in each of the 3 age categories (Tables 4–6). Among toddlers, race and atopic comorbidities were again associated with a higher hazard of new asthma diagnosis in the subsequent 18 months. Additionally, toddlers who were born preterm or had a higher BMI also had a significantly higher hazard of new asthma diagnosis in the subsequent 18 months (Table 4). Female sex significantly increased the hazard of new asthma diagnosis among school-age children and teenagers (Tables 5 and 6). All other risk factors (except for allergic rhinitis in school-age children) had no significant relationship with new asthma diagnosis among school-age children and teenagers (Tables 5 and 6).
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 9351 | (Ref) | (Ref) | (Ref) |
Positive | 885 | 0.92 | 0.64–1.31 | .63 |
Sex | ||||
Male | 5531 | (Ref) | (Ref) | (Ref) |
Female | 4705 | 0.83 | 0.68–1.02 | .079 |
Race | ||||
White | 5854 | (Ref) | (Ref) | (Ref) |
Blacka | 2408 | 1.78 | 1.28–2.46 | .001a |
Asian | 341 | 1.29 | 0.72–2.34 | .39 |
Multiple racesa | 502 | 1.58 | 1.02–2.49 | .043a |
Othera | 1131 | 1.52 | 1.08–2.12 | .015a |
Insurance | ||||
Non-Medicaid | 7055 | (Ref) | (Ref) | (Ref) |
Medicaid | 3181 | 1.00 | 0.77–1.28 | .98 |
Childhood Opportunity Index classification | ||||
Very low | 2238 | (Ref) | (Ref) | (Ref) |
Low | 688 | 0.88 | 0.57–1.37 | .58 |
Moderate | 1179 | 0.78 | 0.51–1.19 | .25 |
High | 2031 | 1.05 | 0.72–1.53 | .81 |
Very high | 4010 | 1.18 | 0.83–1.68 | .36 |
Gestational age | ||||
Full term | 7005 | (Ref) | (Ref) | (Ref) |
Preterma | 1262 | 1.48 | 1.13–1.93 | .005a |
Post-term | 875 | 1.01 | 0.69–1.47 | .98 |
Unknown | 1094 | 1.06 | 0.76–1.50 | .72 |
BMIa | 10 208 | 1.13 | 1.07–1.19 | <.001a |
Atopic dermatitis | ||||
No atopic dermatitis | 8724 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 1512 | 1.01 | 0.77–1.31 | .97 |
Food allergies | ||||
No food allergies | 8679 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 1557 | 1.38 | 1.08–1.75 | .009a |
Allergic rhinitis | ||||
No allergic rhinitis | 7081 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 3155 | 2.65 | 2.15–3.27 | <.001a |
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 9351 | (Ref) | (Ref) | (Ref) |
Positive | 885 | 0.92 | 0.64–1.31 | .63 |
Sex | ||||
Male | 5531 | (Ref) | (Ref) | (Ref) |
Female | 4705 | 0.83 | 0.68–1.02 | .079 |
Race | ||||
White | 5854 | (Ref) | (Ref) | (Ref) |
Blacka | 2408 | 1.78 | 1.28–2.46 | .001a |
Asian | 341 | 1.29 | 0.72–2.34 | .39 |
Multiple racesa | 502 | 1.58 | 1.02–2.49 | .043a |
Othera | 1131 | 1.52 | 1.08–2.12 | .015a |
Insurance | ||||
Non-Medicaid | 7055 | (Ref) | (Ref) | (Ref) |
Medicaid | 3181 | 1.00 | 0.77–1.28 | .98 |
Childhood Opportunity Index classification | ||||
Very low | 2238 | (Ref) | (Ref) | (Ref) |
Low | 688 | 0.88 | 0.57–1.37 | .58 |
Moderate | 1179 | 0.78 | 0.51–1.19 | .25 |
High | 2031 | 1.05 | 0.72–1.53 | .81 |
Very high | 4010 | 1.18 | 0.83–1.68 | .36 |
Gestational age | ||||
Full term | 7005 | (Ref) | (Ref) | (Ref) |
Preterma | 1262 | 1.48 | 1.13–1.93 | .005a |
Post-term | 875 | 1.01 | 0.69–1.47 | .98 |
Unknown | 1094 | 1.06 | 0.76–1.50 | .72 |
BMIa | 10 208 | 1.13 | 1.07–1.19 | <.001a |
Atopic dermatitis | ||||
No atopic dermatitis | 8724 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 1512 | 1.01 | 0.77–1.31 | .97 |
Food allergies | ||||
No food allergies | 8679 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosisa | 1557 | 1.38 | 1.08–1.75 | .009a |
Allergic rhinitis | ||||
No allergic rhinitis | 7081 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 3155 | 2.65 | 2.15–3.27 | <.001a |
There were 10 236 total individuals, and 10 118 were used in the Cox regression analysis.
Significant test results given α = .05.
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 9416 | (Ref) | (Ref) | (Ref) |
Positive | 1272 | 1.29 | 0.78–2.08 | .30 |
Sex | ||||
Male | 5595 | (Ref) | (Ref) | (Ref) |
Femalea | 5092 | 1.71 | 1.20–2.42 | .003a |
Race | ||||
White | 6536 | (Ref) | (Ref) | (Ref) |
Black | 2123 | 1.16 | 0.64–2.07 | .63 |
Asian | 343 | 1.21 | 0.49–3.02 | .68 |
Multiple races | 461 | 1.44 | 0.69–3.03 | .33 |
Other | 1225 | 0.70 | 0.36–1.37 | .30 |
Insurance | ||||
Non-Medicaid | 8102 | (Ref) | (Ref) | (Ref) |
Medicaid | 2586 | 0.80 | 0.50–1.26 | .33 |
Childhood Opportunity Index classification | ||||
Very low | 1822 | (Ref) | (Ref) | (Ref) |
Low | 601 | 0.56 | 0.21–1.49 | .25 |
Moderate | 953 | 0.99 | 0.48–2.06 | .98 |
High | 2104 | 0.94 | 0.48–1.84 | .87 |
Very high | 5127 | 0.76 | 0.40–1.46 | .41 |
Gestational age | ||||
Full term | 6554 | (Ref) | (Ref) | (Ref) |
Preterm | 1112 | 0.58 | 0.28–1.20 | .14 |
Post-term | 779 | 0.85 | 0.41–1.75 | .65 |
Unknowna | 2243 | 1.55 | 1.04–2.30 | .029a |
BMI | 10 620 | 0.99 | 0.95–1.04 | .80 |
Atopic dermatitis | ||||
No atopic dermatitis | 9350 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 1338 | 1.01 | 0.61–1.69 | .96 |
Food allergies | ||||
No food allergies | 8949 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosis | 1739 | 1.19 | 0.78–1.83 | .42 |
Allergic rhinitis | ||||
No allergic rhinitis | 5527 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 5161 | 2.29 | 1.57–3.35 | <.001a |
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 9416 | (Ref) | (Ref) | (Ref) |
Positive | 1272 | 1.29 | 0.78–2.08 | .30 |
Sex | ||||
Male | 5595 | (Ref) | (Ref) | (Ref) |
Femalea | 5092 | 1.71 | 1.20–2.42 | .003a |
Race | ||||
White | 6536 | (Ref) | (Ref) | (Ref) |
Black | 2123 | 1.16 | 0.64–2.07 | .63 |
Asian | 343 | 1.21 | 0.49–3.02 | .68 |
Multiple races | 461 | 1.44 | 0.69–3.03 | .33 |
Other | 1225 | 0.70 | 0.36–1.37 | .30 |
Insurance | ||||
Non-Medicaid | 8102 | (Ref) | (Ref) | (Ref) |
Medicaid | 2586 | 0.80 | 0.50–1.26 | .33 |
Childhood Opportunity Index classification | ||||
Very low | 1822 | (Ref) | (Ref) | (Ref) |
Low | 601 | 0.56 | 0.21–1.49 | .25 |
Moderate | 953 | 0.99 | 0.48–2.06 | .98 |
High | 2104 | 0.94 | 0.48–1.84 | .87 |
Very high | 5127 | 0.76 | 0.40–1.46 | .41 |
Gestational age | ||||
Full term | 6554 | (Ref) | (Ref) | (Ref) |
Preterm | 1112 | 0.58 | 0.28–1.20 | .14 |
Post-term | 779 | 0.85 | 0.41–1.75 | .65 |
Unknowna | 2243 | 1.55 | 1.04–2.30 | .029a |
BMI | 10 620 | 0.99 | 0.95–1.04 | .80 |
Atopic dermatitis | ||||
No atopic dermatitis | 9350 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 1338 | 1.01 | 0.61–1.69 | .96 |
Food allergies | ||||
No food allergies | 8949 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosis | 1739 | 1.19 | 0.78–1.83 | .42 |
Allergic rhinitis | ||||
No allergic rhinitis | 5527 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosisa | 5161 | 2.29 | 1.57–3.35 | <.001a |
There were 10 688 total individuals, and 10 541 were used in the Cox regression analysis.
Significant test results given α = .05.
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 5509 | (Ref) | (Ref) | (Ref) |
Positive | 990 | 0.47 | 0.17–1.31 | .15 |
Sex | ||||
Male | 3188 | (Ref) | (Ref) | (Ref) |
Femalea | 3311 | 2.04 | 1.14–3.66 | .016a |
Race | ||||
White | 4254 | (Ref) | (Ref) | (Ref) |
Black | 1297 | 0.87 | 0.30–2.51 | .80 |
Asian | 208 | 1.13 | 0.27–4.75 | .87 |
Multiple races | 168 | 0.78 | 0.11–5.79 | .81 |
Other | 572 | 1.84 | 0.82–4.13 | .14 |
Insurance | ||||
Non-Medicaid | 5138 | (Ref) | (Ref) | (Ref) |
Medicaid | 1361 | 1.16 | 0.54–2.48 | .71 |
Childhood Opportunity Index classification | ||||
Very low | 1044 | (Ref) | (Ref) | (Ref) |
Low | 326 | 0.36 | 0.043–2.97 | .34 |
Moderate | 590 | 0.78 | 0.20–2.99 | .72 |
High | 1253 | 1.17 | 0.38–3.63 | .79 |
Very high | 3229 | 0.84 | 0.28–2.56 | .76 |
Gestational age | ||||
Full term | 2829 | (Ref) | (Ref) | (Ref) |
Preterm | 500 | 1.32 | 0.50–3.51 | .58 |
Post-term | 348 | — | — | — |
Unknown | 2822 | 1.10 | 0.61–1.97 | .75 |
BMI | 6452 | 0.95 | 0.88–1.01 | .10 |
Atopic dermatitis | ||||
No atopic dermatitis | 5859 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 640 | 1.13 | 0.44–2.89 | .80 |
Food allergies | ||||
No food allergies | 5492 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosis | 1007 | 0.80 | 0.34–1.91 | .62 |
Allergic rhinitis | ||||
No allergic rhinitis | 2845 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosis | 3654 | 0.82 | 0.47–1.43 | .48 |
Variable . | N . | HR . | 95% CI . | P . |
---|---|---|---|---|
SARS-CoV-2 PCR result | ||||
Negative | 5509 | (Ref) | (Ref) | (Ref) |
Positive | 990 | 0.47 | 0.17–1.31 | .15 |
Sex | ||||
Male | 3188 | (Ref) | (Ref) | (Ref) |
Femalea | 3311 | 2.04 | 1.14–3.66 | .016a |
Race | ||||
White | 4254 | (Ref) | (Ref) | (Ref) |
Black | 1297 | 0.87 | 0.30–2.51 | .80 |
Asian | 208 | 1.13 | 0.27–4.75 | .87 |
Multiple races | 168 | 0.78 | 0.11–5.79 | .81 |
Other | 572 | 1.84 | 0.82–4.13 | .14 |
Insurance | ||||
Non-Medicaid | 5138 | (Ref) | (Ref) | (Ref) |
Medicaid | 1361 | 1.16 | 0.54–2.48 | .71 |
Childhood Opportunity Index classification | ||||
Very low | 1044 | (Ref) | (Ref) | (Ref) |
Low | 326 | 0.36 | 0.043–2.97 | .34 |
Moderate | 590 | 0.78 | 0.20–2.99 | .72 |
High | 1253 | 1.17 | 0.38–3.63 | .79 |
Very high | 3229 | 0.84 | 0.28–2.56 | .76 |
Gestational age | ||||
Full term | 2829 | (Ref) | (Ref) | (Ref) |
Preterm | 500 | 1.32 | 0.50–3.51 | .58 |
Post-term | 348 | — | — | — |
Unknown | 2822 | 1.10 | 0.61–1.97 | .75 |
BMI | 6452 | 0.95 | 0.88–1.01 | .10 |
Atopic dermatitis | ||||
No atopic dermatitis | 5859 | (Ref) | (Ref) | (Ref) |
Atopic dermatitis diagnosis | 640 | 1.13 | 0.44–2.89 | .80 |
Food allergies | ||||
No food allergies | 5492 | (Ref) | (Ref) | (Ref) |
Food allergy diagnosis | 1007 | 0.80 | 0.34–1.91 | .62 |
Allergic rhinitis | ||||
No allergic rhinitis | 2845 | (Ref) | (Ref) | (Ref) |
Allergic rhinitis diagnosis | 3654 | 0.82 | 0.47–1.43 | .48 |
There were 6499 total individuals, and 6395 were used in the Cox regression analysis. —, regression output unavailable as sample size insufficient.
Significant test results given α = .05.
Discussion
In this longitudinal cohort study of children observed 1 year before and 18-months after receiving SARS-CoV-2 PCR testing in the first year of the pandemic, we found that SARS-CoV-2 PCR positivity was not associated with subsequent asthma diagnosis. Although contrary to the literature regarding other respiratory viral infections, particularly RSV and rhinovirus, our findings are consistent with contemporary SARS-CoV-2 literature. For example, recent descriptive epidemiology studies have demonstrated a similar incidence of asthma diagnoses between pre and postpandemic levels, despite widespread transmission of SARS-CoV-2.14 ,15 Additionally, the Long COVID Research Group demonstrated that long-term pulmonary function is not affected in both symptomatic and asymptomatic children infected with the SARS-CoV-2 virus.16
Defining asthma phenotypes in a pediatric cohort is challenging, especially among younger children who may present with a self-limited acute wheezing episode. Our investigation found a higher incidence of asthma among young children, with older age being associated with a significantly lower likelihood of new asthma diagnosis. It is worth noting that prior epidemiologic analyses of our cohort have demonstrated a high cumulative incidence of asthma (15.9%) in the 0 to 5 age group.17 This value approximates that which was found in a recent multistate analysis (17.5%) of the 0 to 5 age group,18 indicating that asthma diagnosis rates in our region are similar to those observed nationally. To minimize the likelihood of transient wheeze confounding our current study, we excluded ICD codes relating to reactive airway and postviral wheeze from our asthma definition. Further, we performed a sensitivity analysis that required more than 1 encounter with an asthma diagnosis. The results of this analysis did not differ significantly from our primary analysis, indicating that transient wheeze episodes likely did not substantially alter our findings.
The lack of observed association between SARS-CoV-2 and incident asthma is also supported by the consistency of the association between other asthma risk factors and new asthma diagnosis in the current study. In the primary analysis, atopic history (specifically food allergies and allergic rhinitis) was associated with a significantly higher risk of asthma diagnosis. Atopy is a well-studied risk factor for the development of asthma given its biologic predisposition to other allergic conditions.19 As demonstrated by our subgroup analyses, these relationships are especially relevant for infants and toddlers when lung development is most vulnerable to complex differences in both biologic and environmental factors. Notably, the subgroup analysis of toddlers yielded significant relationships between preterm gestation and high BMI, both of which increased the risk of developing asthma. These are also well-studied risk factors across the literature, particularly among young children.20 ,21
This study further benefits from a relatively unique period in history in which there was more uniformity in and control over environmental risk factors for asthma development. Principally, the prevalence of other respiratory viral illnesses (including RSV, influenza, and rhinovirus) plummeted during the first year of the pandemic.10 This fact provided a natural opportunity to isolate the effects of SARS-CoV-2 from other respiratory viral infections. Additionally, global improvement in pollution trends and air quality, coupled with greater isolation indoors, may have helped minimize exposure to outdoor toxicants and allergens while also minimizing variability between urban and rural settings.22 Further, the pandemic generated a widespread practice of PCR testing of both symptomatic and asymptomatic individuals (for example, in the cases of preprocedural and preadmission testing). Typically, asymptomatic children are not tested for viral infections, making it difficult to define an epidemiologic comparison group that lacks the exposure of interest. In other words, the fact that healthcare systems were performing the gold-standard test on asymptomatic children during the first year of the pandemic created a natural experimental design that facilitated the current study.
Recent investigation has shed light on possible pathophysiologic mechanisms that may better inform the understanding of our results. SARS-CoV-2 infection tends to elicit an immunophenotype notable for a robust type I inflammatory response characterized by type I interferon and inflammasome activation and subsequent CD8 T cell response.23 This is well demonstrated in the context of the multisystem inflammatory syndrome in children, an immunopathologic response characterized by overactivation of these pathways.23 Infection with the SARS-CoV-2 virus has also been shown to induce eosinopenia regardless of disease severity.24 ,25 This is contrary to infection with RSV or rhinovirus, which have been linked to increases in certain atopic biomarkers, such as peripheral eosinophilia.1 Taken altogether, infection with the SARS-CoV-2 virus may drive the cellular immunophenotype at the level of the airway epithelium to a more type I, as opposed to type II, inflammatory state.26 ,27 Chronic asthma, on the other hand, is primarily mediated by type II inflammation characterized by IL-4, IL-5, IL-13, and eosinophilia.27 ,28
There are inherent limitations to our study. First, the use of EHR data from 1 institution may impact the generalizability of our findings. However, the CCN covers a 3-state region in one of the most populated parts of the country, including urban and suburban settings. Further, the CCN patient population is largely reflective of allergic disease patterns across the country, at least in the younger age group.18 Nonetheless, EHR data are also limited in its reliance upon provider-entered information. For example, our EHR data does not consistently capture certain risk factors, such as family history of asthma and secondhand smoke exposure. In addition, our study relies upon SARS-CoV-2 PCR test results for a 1-year exposure window, which may have misclassified individuals who later obtained a positive PCR or rapid antigen test after the exposure window. However, this misclassification is unlikely to change the direction of the association given the consistent relationship between positive and negative SARS-CoV-2 exposure groups over time in our analyses. Studying the first year of the pandemic was an important feature of our experimental design as it isolated a time interval before the deployment of home rapid antigen SARS-CoV-2 testing. However, patients seeking care during this period may have been less able to delay care until pandemic conditions improved. As such, our study population may have had a higher medical complexity than the US pediatric population as a whole. Because of the limited number of children with severe primary SARS-CoV-2 infection, low subject numbers precluded considering SARS-CoV-2 severity in our analyses. Therefore, future studies may identify subcohorts of subjects where SARS-CoV-2 infection may modify incident asthma risk. Finally, as our exposure window preceded the evolution of several SARS-CoV-2 variants (eg, ο), it is possible that these later SARS-CoV-2 strains may influence asthma risk differently.
Conclusions
Taken altogether, we find a consistent finding across all pediatric age groups that SARS-CoV-2 PCR positivity does not confer an additional risk for asthma diagnosis, at least within the first 18 months after PCR test. This information may be useful for families and providers alike in the prognostication of the long-term respiratory effects after SARS-CoV-2 infection in children. As one of the first population-based studies to investigate SARS-CoV-2 infection in asthma-naïve children, this study has several implications. This study refines our knowledge regarding the long-term respiratory effects of SARS-CoV-2 infection in children and adds to our growing understanding of how SARS-CoV-2 may influence asthma development in pediatric patients. Future work should consider SARS-CoV-2 infection severity and stratify between asymptomatic and symptomatic children and between those with mild, moderate, or severe symptoms. Future work should also consider asthma by the persistence and severity of symptoms in children with both new and prior asthma diagnoses.
Drs Senter and Aisenberg designed the study, coordinated and supervised data collection, conducted the initial analyses, and drafted the initial manuscript; Mr Dudley designed the study and data collection instruments, collected data, conducted the initial analyses, and critically reviewed and revised the manuscript; Mr Luan and Dr Huang coordinated and supervised data collection, and critically reviewed and revised the manuscript; Drs Kenyon and Hill conceptualized and designed the study, coordinated and supervised data collection, and critically reviewed and revised the manuscript; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.
FUNDING: Dr Kenyon is supported by the National Institutes of Health (NIH, grant K23HL136842). Asthma research in the Hill Laboratory is supported by the National Institutes of Health (grant R01 HL162715) and the Children’s Hospital of Philadelphia Research Institute. The other authors received no additional funding. This information or content and conclusions are those of the authors and do not represent the official views of, nor an endorsement by, the funders. None of the funding sources had a role in the design or conduct of the study.
CONFLICT OF INTEREST DISCLOSURES: The authors have no conflicts of interest relevant to this article to disclose.
Comments