Vaccine Effectiveness Against Long COVID in Children: A Report from the RECOVER EHR Cohort
==========================================================================================

* Hanieh Razzaghi
* Christopher B. Forrest
* Kathryn Hirabayashi
* Qiong Wu
* Andrea Allen
* Suchitra Rao
* Yong Chen
* H. Timothy Bunnell
* Elizabeth A. Chrischilles
* Lindsay G. Cowell
* Mollie R. Cummins
* David A Hanauer
* Miranda Higginbotham
* Benjamin D. Horne
* Carol R. Horowitz
* Ravi Jhaveri
* Susan Kim
* Aaron Mishkin
* Jennifer A Muszynski
* Susanna Naggie
* Nathan M. Pajor
* Anuradha Paranjape
* Hayden T Schwenk
* Marion R. Sills
* Yacob G. Tedla
* David A. Williams
* Charles Bailey
* the RECOVER Consortium* A complete list of group members appears in the supplementary materials

## Abstract

**Objective** Vaccination reduces the risk of acute COVID-19 in children, but it is less clear whether it protects against long COVID. We estimated vaccine effectiveness (VE) against long COVID in children aged 5-17 years.

**Methods** This retrospective cohort study used data from 17 health systems in the RECOVER PCORnet electronic health record (EHR) Program for visits between vaccine availability, and October 29, 2022. Conditional logistic regression was used to estimate VE against long COVID with matching on age group (5-11, 12-17) and time period and adjustment for sex, ethnicity, health system, comorbidity burden, and pre-exposure health care utilization. We examined both probable (symptom-based) and diagnosed long COVID in the year following vaccination.

**Results** The vaccination rate was 56% in the cohort of 1,037,936 children. The incidence of probable long COVID was 4.5% among patients with COVID-19, while diagnosed long COVID was 0.7%. Adjusted vaccine effectiveness within 12 months was 35.4% (95 CI 24.5 – 44.5) against probable long COVID and 41.7% (15.0 – 60.0) against diagnosed long COVID. VE was higher for adolescents 50.3% [36.3 – 61.0]) than children aged 5-11 (23.8% [4.9 – 39.0]). VE was higher at 6 months (61.4% [51.0 – 69.6]) but decreased to 10.6% (–26.8 – 37.0%) at 18-months.

**Discussion** This large retrospective study shows a moderate protective effect of SARS-CoV-2 vaccination against long COVID. The effect is stronger in adolescents, who have higher risk of long COVID, and wanes over time. Understanding VE mechanism against long COVID requires more study, including EHR sources and prospective data.

**Article Summary** Vaccination against COVID-19 has a protective effect against long COVID in children and adolescents. The effect wanes over time but remains significant at 12 months.

**What’s Known on This Subject** Vaccines reduce the risk and severity of COVID-19 in children. There is evidence for reduced long COVID risk in adults who are vaccinated, but little information about similar effects for children and adolescents, who have distinct forms of long COVID.

**What This Study Adds** Using electronic health records from US health systems, we examined large cohorts of vaccinated and unvaccinated patients <18 years old and show that vaccination against COVID-19 is associated with reduced risk of long COVID for at least 12 months.

**Contributors’ Statement** Drs. Hanieh Razzaghi and Charles Bailey conceptualized and designed the study, supervised analyses, drafted the initial manuscript, and critically reviewed and revised the manuscript.

Drs. Christopher Forrest and Yong Chen designed the study and critically reviewed and revised the manuscript.

Ms. Kathryn Hirabayashi, Ms. Andrea Allen, and Dr. Qiong Wu conducted analyses, and critically reviewed and revised the manuscript.

Drs. Suchitra Rao, H Timothy Bunnell, Elizabeth A. Chrischilles, Lindsay G. Cowell, Mollie R. Cummins, David A. Hanauer, Benjamin D. Horne, Carol R. Horowitz, Ravi Jhaveri, Susan Kim, Aaron Mishkin, Jennifer A. Muszynski, Susanna Nagie, Nathan M. Pajor, Anuradha Paranjape, Hayden T. Schwenk, Marion R. Sills, Yacob G. Tedla, David A. Williams, and Ms. Miranda Higginbotham critically reviewed and revised the manuscript.

All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

**Authorship statement** Authorship has been determined according to ICMJE recommendations.

## Introduction

SARS-CoV-2 has infected over 676 million people, and has been associated with 6.9 million deaths1. While severity has been lower in children than adults, statistics for the United States (US) indicate >2100 deaths in children aged <21 years with COVID-192 and up to 2.1% of cases classified as severe disease3. Further, children may experience long-term health burdens after infection, collectively referred to as post-acute sequelae of SARS-CoV-2 (PASC) or “long COVID” 4 5,6. These may represent effects of host response to the original infection in addition to or instead of direct pathogen effects7–9. Long COVID is heterogeneous and likely underdiagnosed, but symptoms such as brain fog, dyspnea, gastrointestinal dysfunction, generalized pain, and fatigue can cause significant burden for children, even after mild COVID-194,5,10,11. Recent research has identified specific post-acute symptoms consistent with long COVID 12. Prevalence estimates vary widely10,12–14. Symptoms of long COVID have been reported less often and for shorter duration in children compared to adults 11, and higher rates have been reported for adolescents than younger children13. It is difficult to establish how much this results from differential reporting of symptoms at different ages, greater difficulty distinguishing long COVID from other childhood illnesses or effects of the pandemic (*e.g.* disruption of seasonal viral patterns, or of school progress)13. Despite the paucity of data about long COVID in children, developing mitigation strategies is an important concern not only for affected children but in societal discussions15 and for shaping public health policy16–18.

Vaccines provide the best opportunity for both reducing the risk of severe disease19–24 as well as short-term complications of acute infection. In adults, some reports25–29 identified a modest benefit associated with vaccination in preventing long COVID in adults, while others did not30,31. Little data are available on the effect of vaccination on long COVID in children; in one UK study including 6500 children with COVID-19, no association was seen4. However, to our knowledge, there are no studies assessing clinical data to address this question. Efforts to learn from clinical care at a large enough scale to identify underrepresented patterns provides an important complement to multicenter prospective cohorts. The NIH Researching COVID to Enhance Recovery (RECOVER) Initiative, which seeks to understand, treat, and prevent the post-acute sequelae of SARS-CoV-2 infection (PASC, or long COVID), is employing both strategies to better elucidate the condition, identify effective ways to treat patients, and support families impacted by long COVID. We report results from a large-scale, national collaboration of health systems from the PCORnet® network contributing electronic health record (EHR) data as part of RECOVER’s real-world program46.

## Methods

### Data sources

We conducted a retrospective cohort study based on EHR data from 41 health systems for children with COVID-19, a diagnosis of acute respiratory illness, fever, dyspnea, or cough, or a SARS-CoV-2 vaccine. Data, extracted quarterly, were transformed to either the PCORnet (v6.0) or Observational Medical Outcomes Partnership (OMOP) (v5.3) data models, and forwarded to the RECOVER-PCORnet Coordinating Center. For this study, we used the s8 version of the data, collected in April 2023, which comprises 25,389,893 patients.

Institutional Review Board oversight was provided by the Biomedical Research Alliance of New York, protocol # 21-08-508-380.

### Study Sample

Patients were eligible if they had an in-person visit between vaccine eligibility and October 29, 2022, and were between 5 and <18 years at the date of vaccine eligibility. The date of SARS-CoV-2 vaccine eligibility was defined as January 01, 2021, for patients aged 12 – <18 years, and October 29, 2021, for patients aged 5 – <12 years, corresponding to availability of vaccines. A small number of patients were excluded for documented COVID-19 before 28 days after their first vaccine dose.

For vaccinated patients, cohort entry was the date of their first vaccine dose, while for unvaccinated patients, the date of an in-person visit after vaccine eligibility was selected randomly, resulting in matched distribution of entry dates between cohorts. In addition, we required 1 contact between 8 days and 2 years prior to cohort entry to assess baseline health characteristics. To determine long COVID outcomes, we also required at least 1 in-person visit between 28 and 179 days after cohort entry.

### Outcomes

The main outcomes were *diagnosed* or *probable* long COVID. A patient had *diagnosed* long COVID if they had two or more visits with diagnosis codes specific for long COVID; a single diagnosis code classified children as having *probable* long COVID. In addition, to account for incomplete availability or use of long-COVID-specific diagnosis codes, especially for patients with early signs of long COVID, we also classified a patient as having *probable* long COVID if they had COVID-19 (SARS-CoV-2 PCR or antigen positive or COVID-19-specific diagnosis codes) plus at least two long-COVID-compatible diagnoses 28-179 days after infection. These diagnoses were identified in prior work as post-acute associations with COVID-1912,32 (Supplemental Table S3).

Some results are reported for pre-omicron and omicron eras. This applies only to adolescents, as vaccines for younger children were not available substantially before the omicron era. To maintain consistency across analyses, November 29, 2021, the first day a younger child could have had a qualifying post-vaccine infection, was used as the dividing line between eras; at that point omicron was circulating in the US, and became the dominant strain a few weeks later. For each patient, the date of a SARS-CoV-2 infection following cohort entry, if one occurred, was used to place their results in either of the two eras; if the patient never developed COVID-19, the date of cohort entry was used.

### Exposure

The primary exposure was receipt of at least one dose of SARS-CoV-2 vaccine. Vaccine doses were ascertained using a vaccine code (CVX or RxNorm) designating an administered or patient-reported dose, a procedure code indicating SARS-CoV-2 vaccination, or a vaccine description containing the string “COVID” or “SARS”. Where available, data was retrieved from state/regional vaccine information systems via hospital operational procedures. Adequate vaccine data capture was defined as a site-reported vaccination rate at least 60% of the CDC-estimated rate2 for its coverage area, based on addresses of patients seen during the pandemic. Where dose counts were used, all recorded doses within 14 days of a base dose were considered a single dose.

### Demographic and clinical data

Age, sex, and ethnicity were taken from EHRs. Comorbidity data used the body system taxonomy of the Pediatric Medical Complexity Algorithm (PMCA)14. Diagnoses received up to 3 years before cohort entry were considered. To classify patients as complex-chronic, chronic, or no chronic condition, the “more conservative” version of the algorithm was applied.

### Main analysis

We matched vaccinated and unvaccinated patients on age group at cohort entry (5-11, 12-17) and time period (six-month cohort entry period). We limited long COVID observations to 12 months following cohort entry to account for waning of vaccine effectiveness (VE) over time.

### Sensitivity analyses

To evaluate potential impact of attrition due to matching, we performed a sensitivity analysis using inverse probability of treatment weighting that incorporated age group, time of cohort entry, sex, ethnicity, health system, number of baseline visits to the health system, PMCA-defined body system(s) affected by chronic illness, progressive disease, and malignancy. Further, we estimated VE after changing the window for outcomes to 6 months and 18 months post cohort entry and restricting the vaccinated cohort to those with at least 2 recorded doses.

### Cohort Comparison

For unadjusted analyses, follow-up time was computed from cohort entry until the end of observation. We calculated an incidence rate ratio (IRR) comparing vaccinated to unvaccinated groups and derived 95% CIs using the estimation of Wald33. Corresponding VE was computed as 100 × (1 IRR).

For adjusted analyses, we used conditional logistic regression to account for stratification. Covariates included sex, ethnicity, health system, baseline utilization (0, 1-5, 6-10, 11-24, 25-49, 50-99, or 100+ encounters in the 3 years prior to cohort entry), the number of PMCA body systems affected, and the presence of progressive or malignant condition. Since the prevalence of long COVID was low, IRR ≈ OR, so the corresponding VE was computed as 100 × (1-OR).

### Mediation Analyses

To understand the direct and indirect impact of vaccination on long COVID, we conducted causal mediation analyses, where the treatment was pre-infection vaccination, the mediator was SARS-CoV-2 infection, and the outcome was long COVID. We fit two logistic regression models, the first of treatment on the mediator, and the second of treatment and mediator on the outcome, adjusting for the same set of covariates. From these models, we estimated the effect and CIs of vaccination on long COVID independent of infection, and of vaccination on long COVID mediated by infection.

### Statistical Tools

Data analysis was performed using R17 version 4.0.2. Exact matching used the matchIt package 34 version 4.4.0, while IPTW were computed via twang 35 version 2.5. Conditional logistic regression used the clogit() algorithm implemented in the survival package 36 version 3.5. Significant effects were identified based on a threshold of *a*<0.05.

## Results

### Study Cohort

Immunization data from 17 health systems was sufficiently robust for inclusion. Cohort formation is described in Figure 1: the 5–11-year-old group comprised 480,498 children, while the 12-17 year-old group comprised 719,519 children/adolescents. Overall, 55% received at least one SARS-CoV-2 vaccine, and 84% of vaccinated children received two or more doses. Compared to unvaccinated children, vaccinated patients were older, had lower baseline healthcare utilization, and fewer chronic conditions. Girls and children who identified as Asian or Hispanic were more common in the vaccinated group, while those identifying as Black and White were less common (Table 1). Descriptive statistics for the matched cohorts are provided in Supplemental Table S1. Vaccinated patients had slightly fewer comorbidities and healthcare utilization than unvaccinated patients. The prevalence of probable long COVID was 0.3% in the cohort overall; for children with COVID-19 after cohort entry, prevalence of diagnosed long COVID was 0.8%, and adding probable cases raised the prevalence of long COVID to 4.5%.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/28/2023.09.27.23296100/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/09/28/2023.09.27.23296100/F1)

Figure 1: 
Study Sample Selection. Steps in the construction of the overall study sample are shown, including demographic and vaccine-related eligibility criteria. At each step, the number of children satisfying all criteria to that point is shown, as is the percentage of children in the previous step who were retained. a Vaccination rates ≥70% of CDC regional estimate

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/09/28/2023.09.27.23296100/T1)

Table 1. Cohort Description for Full Cohort

### Effect of Vaccination on long COVID in Children and Adolescents

While the likelihood that a recorded vaccine is incorrectly attributed to a patient is very low, the risk that dose count is incomplete in reported data is well-described37, and discontinuation before a second dose is uncommon. It is therefore more accurate in this study design to include patients with evidence of vaccination, and, because any compensatory bias introduced thereby will tend to reduce estimates, more responsible to report results that do not overestimate vaccine effectiveness. We analyzed all children with evidence of vaccination prior to SARS-CoV-2 infection. Results were very similar across approaches: VE was 35% (95% CI 24– 45%) for the combined age groups (Table 2), with greater effect in adolescents than younger children. For children known to be completely vaccinated (*i.e.* 2+ doses) prior to SARS-CoV-2 infection, VE was 35% (95% CI 35– 53%) against probable long COVID within 12 months (Figure 2), with consistent findings within age strata, supporting the validity of our approach. For comparison, VE estimates against SARS-CoV-2 infection are shown in Supplemental Figure S1. We observed decreased likelihood of long COVID among boys and children with non-Hispanic Black or Asian ethnicity (Supplementary Table S2). The presence of chronic conditions had minimal effect, though likelihood of long COVID rose with increasing baseline utilization.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/28/2023.09.27.23296100/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/09/28/2023.09.27.23296100/F2)

Figure 2. 
Vaccine effectiveness estimates for children receiving at least 2 doses of mRNA vaccine in their primary series. Analyses were otherwise structured identically to the main analysis.

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/09/28/2023.09.27.23296100/T2)

Table 2. SARS-CoV-2 Vaccine Effectiveness Against Long COVID

The protective effect of vaccination appeared to wane over time. We observed higher VE against long COVID within 6 months of vaccination than 12 months, while extending the observation period to 18 months revealed a further diminution (Figure 3). Protection was also seen in children who started vaccination after recuperating from COVID-19 (Supplemental Figure S2).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/28/2023.09.27.23296100/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/09/28/2023.09.27.23296100/F3)

Figure 3. 
Durability of vaccine effectiveness of SARS-CoV-2 vaccine against long COVID. Results are shown for PASC occurring within 6, 12, or 18 months of vaccination.

One-to-one matching resulted in loss of 27% of unvaccinated and 64% of vaccinated patients due to the high overall vaccination rate. To assess the effect of this on our VE estimates, we used IPTW to adjust for multiple potential confounders without decreasing cohort size. Once again, a protective effect was seen, though with a smaller magnitude of 24% (95% CI 20 – 28%) for vaccination on long COVID up to a year later (Supplemental Figure S3).

To assess where in the path to long COVID the vaccine might have impact, we performed mediation analysis, estimating effect of vaccine on long COVID via change in COVID-19 risk (indirect effect) separately from other effects on long COVID (direct effect). For 5–11-year-olds, who became vaccine eligible just before the omicron variant became prevalent, the long COVID-associated OR for indirect effect on probable long COVID was 0.67 (95% CI 0.64 – 0.70), while for direct effect the OR was 1.34 (95% CI 1.13 – 1.59). For adolescents, results were similar: in the pre-omicron era the indirect effect OR was 0.46 (95% CI 0.45 – 0.48) and direct effect OR was 1.32 (95% CI 1.12 – 1.55); in the omicron era, the analogous ORs were 0.45 (95% CI 0.41 – 0.49) and 0.99 (95% CI 0.73 – 1.34). These results suggest that throughout the period of pediatric use, the effectiveness of vaccine against long COVID is closely tied to its effectiveness against the antecedent COVID-19 episode.

## Discussion

To our knowledge, this national-scale study, with a diverse cohort of over 1 million children, is the first to estimate the VE against long COVID in children. Our findings suggest that SARS-CoV-2 vaccine is associated with reduced risk of Long COVID in children and adolescents. The protective effect was greater among adolescents than younger children, reflecting activity against delta and omicron strains, while use for younger children is limited to the omicron era. Since diagnostic codes for long COVID (*e.g.* ICD-10-CM U09.9) are rare in children, we have examined both specific diagnosis as an outcome, and a broader definition based on the repeated presence from 1-6 months after infection of symptoms that occur more frequently following COVID-19 than in uninfected children12,32. The effect of vaccine is similar for both outcomes and is also seen in children vaccinated after an initial episode of COVID-19, suggesting that vaccination confers some benefit against future late effects. Mediation analysis suggests that the protective effect of vaccination may be primarily via reduction in the occurrence or severity of infection.

Our finding of a protective effect against long COVID in children are consistent with those observed in adults, while the VE estimates in our study are lower than studies assessing the short-term effect of vaccines in the one week to one month after vaccination 29. Other reasons are varied definitions of long COVID between these studies, and different methodological approaches. Our observed higher VE for adolescents is similar to studies examining acute infection 38. Further analyses are required to explore whether other factors contribute to this difference, such as long COVID symptoms being more frequently reported or documented in adolescents.

Waning effects of vaccine against long COVID are to be expected 39, similar to studies against acute infection 40 41 42. The reasons include reduced neutralization efficiency against newer strains and waning of pre-existing antibodies over time 43. It is also possible that long COVID starting several months after acute infection, or delayed diagnosis or recognition of long COVID, may contribute to this phenomenon. Further, waning effect may reflect the impact of successive infections rather than late-onset long COVID. This is difficult to assess analytically, given fewer COVID-19 episodes are being documented in the EHR during the omicron era. Evaluations of VE over time are needed to address whether these waning effects can be overcome through booster or annual vaccine doses.

VE studies using EHR data rely on accurate capture of clinical and vaccine data, with adjustment for important confounders to minimize potential bias 44 45. Our focus on data quality, including comparison to the CDC’s county-level vaccination data, reduced our cohort size but yielded more reliable vaccination data. Our assessment of VE against acute SARS-CoV-2 was similar to VE estimates from other pediatric studies, further validating our findings and approach. Finally, we performed matching on age and time period to limit secular effects and conducted sensitivity analyses including inverse probability of treatment weighting, with similar findings.

Nonetheless, several challenges exist with the identification of long COVID in children using EHR data, which may confound assessment of our outcome. There has been low uptake of the long COVID (ICD10-CM U09.9) diagnosis code in pediatrics; it is more often seen in older children and associated with respiratory, cardiovascular, or fatigue subtypes of long COVID 32. For this reason, we supplemented our outcome definition with a symptom-based computable phenotype definition of long COVID we have developed, which includes features related to long COVID based on statistical association and chart review validation12,32. However, long COVID presentations are heterogeneous, and symptoms and conditions may overlap with other childhood illnesses, with the potential for misclassification of long COVID in our study. These limitations on classification accuracy are expected to affect both vaccinated and unvaccinated cohorts, similarly, so would not confound the comparison of the two.

There are other limitations that warrant discussion. First, our study relies on secondary use of EHR data, and is subject to bias due to differential access to and utilization of healthcare during the pandemic, as well as clinical practice variation across sites. We required both pre– and post-vaccine visits to establish a baseline of utilization, which may create bias toward sicker patients, but importantly we did not observe major imbalances in utilization rates between cohorts. Next, testing practices changed over time, with an increased reliance on home antigen testing in later phases of the pandemic, potentially leading to misclassification of symptom based long COVID. In addition, a large proportion of patients in our data with a long COVID diagnosis code do not have documentation in the EHR of the antecedent COVID-19 episode, which prevented use of severity as an adjustment to vaccine effect.

## Conclusion

Our results provide substantial evidence in a large and diverse cohort of children receiving health care for protective effect in children 5 years and older. This study adds to the growing body of knowledge about mitigating effect of vaccines on COVID-19, while demonstrating the need for further research utilizing a range of designs to examine the protective effect of these vaccines against subsequent strains and to help guide vaccine policy.

## Supporting information

Supplement [[supplements/296100_file05.docx]](pending:yes)

RECOVER Consortium [[supplements/296100_file06.docx]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

[https://redcap.chop.edu/surveys/?s=88WJ8WAFPJ](https://redcap.chop.edu/surveys/?s=88WJ8WAFPJ) 

## Acknowledgments

We are grateful to the health systems, clinicians, and patients at member sites without whose trust research of this type would not be possible. We are also grateful to clinical data teams at each RECOVER site involved in data extraction and integration, as well as colleagues at the RECOVER Coordinating Center, particularly Susan Hague, for maintaining the RECOVER data infrastructure. Additionally, we would like to thank the National Community Engagement Group (NCEG), all patient, caregiver and community Representatives, and all the participants enrolled in the RECOVER Initiative

## Footnotes

*   Disclosures: Dr. Mollie Cummins is employed by Doxy.me Inc., a commercial telemedicine platform provider. Dr. Benjamin Horne is a member of the advisory boards of Opsis Health and Lab Me Analytics, and a consultant to Pfizer (regarding clinical risk scores; funds paid to Intermountain). Dr. Susanna Naggie reports research grants from Gilead Sciences and AbbVie, scientific advisor/stock options from Vir Biotechnologies, consulting with no financial payment from Pardes Biosciences and Silverback Therapeutics, DSMB fees from Personal Health Insights, Inc, event adjudication committee fees from BMS/PRA outside the submitted work. Dr. Mishkin receives Grant support from Pfizer paid directly to Institution Advisory Board for Takeda. Dr. Jhaveri is a consultant for AstraZeneca, Seqirus, Dynavax, receives an editorial stipend from Elsevier and Pediatric Infectious Diseases Society. Dr. Rao reports prior grant support from GSK and Biofire and is a consultant for Sequiris. The remaining authors have no disclosures

*   Funding/Support: This research was funded by the National Institutes of Health (NIH) Agreement OT2HL161847-01 as part of the Researching COVID to Enhance Recovery (RECOVER) program of research

*   Disclaimer: The views and conclusions contained in this document are those of the authors and should not be interpreted as representing the official policies, either expressed or implied, of the NIH

*   Abbreviations: PCR – polymerase chain reaction; EHR – electronic health record; CI – 95% confidence intervals; ICD-10-CM – International Classification of Diseases, version 10, clinical modification; PMCA – Pediatric Medical Complexity Algorithm; US – United States of America; VE – vaccine effectiveness; VIS – vaccine information system

*   This revised version of the manuscript includes minor formatting changes and an accurate author list.

*   Received September 27, 2023.
*   Revision received September 28, 2023.
*   Accepted September 28, 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.Johns Hopkins University. COVID-19 dashboard by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU). Accessed March 04, 2023. [https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html)
    
    

2.  2.Centers for Disease Control and Prevention (CDC). CDC COVID Data Tracker. Accessed March 04, 2023. [https://covid.cdc.gov/covid-data-tracker/#demographics](https://covid.cdc.gov/covid-data-tracker/#demographics)
    
    

3.  3.Forrest CB, Burrows EK, Mejias A, et al. Severity of Acute COVID-19 in Children <18 Years Old March 2020 to December 2021. Pediatrics. Apr 1 2022;149(4)doi:10.1542/peds.2021-055765
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1542/peds.2021-055765&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35322270&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

4.  4.Pinto Pereira SM, Nugawela MD, Rojas NK, et al. Post-COVID-19 condition at 6 months and COVID-19 vaccination in non-hospitalised children and young people. Arch Dis Child. Apr 2023;108(4):289–295. doi:10.1136/archdischild-2022-324656
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImFyY2hkaXNjaGlsZCI7czo1OiJyZXNpZCI7czo5OiIxMDgvNC8yODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8yOC8yMDIzLjA5LjI3LjIzMjk2MTAwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

5.  5.Trapani G, Verlato G, Bertino E, et al. Long COVID-19 in children: an Italian cohort study. Ital J Pediatr. Jun 3 2022;48(1):83. doi:10.1186/s13052-022-01282-x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13052-022-01282-x&link_type=DOI) 

6.  6.Davis HE, McCorkell L, Vogel JM, Topol EJ. Long COVID: major findings, mechanisms and recommendations. Nat Rev Microbiol. Mar 2023;21(3):133–146. doi:10.1038/s41579-022-00846-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41579-022-00846-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36639608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

7.  7.Kuczborska K, Buda P, Ksiazyk J. Long-COVID in immunocompromised children. Eur J Pediatr. Sep 2022;181(9):3501–3509. doi:10.1007/s00431-022-04561-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00431-022-04561-1&link_type=DOI) 

8.  8.Pierce CA, Herold KC, Herold BC, et al. COVID-19 and children. Science. Sep 9 2022;377(6611):1144–1149. doi:10.1126/science.ade1675
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.ade1675&link_type=DOI) 

9.  9.Molloy EJ, Nakra N, Gale C, Dimitriades VR, Lakshminrusimha S. Multisystem inflammatory syndrome in children (MIS-C) and neonates (MIS-N) associated with COVID-19: optimizing definition and management. Pediatr Res. May 2023;93(6):1499–1508. doi:10.1038/s41390-022-02263-w
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41390-022-02263-w&link_type=DOI) 

10. 10.Say D, Crawford N, McNab S, Wurzel D, Steer A, Tosif S. Post-acute COVID-19 outcomes in children with mild and asymptomatic disease. Lancet Child Adolesc Health. Jun 2021;5(6):e22–e23. doi:10.1016/S2352-4642(21)00124-3
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2352-4642(21)00124-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33891880&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

11. 11.Mizrahi B, Sudry T, Flaks-Manov N, et al. Long covid outcomes at one year after mild SARS-CoV-2 infection: nationwide cohort study. BMJ. Jan 11 2023;380:e072529. doi:10.1136/bmj-2022-072529
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjIwOiIzODAvamFuMTFfMTAvZTA3MjUyOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzI4LzIwMjMuMDkuMjcuMjMyOTYxMDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

12. 12.Rao S, Lee GM, Razzaghi H, et al. Clinical Features and Burden of Postacute Sequelae of SARS-CoV-2 Infection in Children and Adolescents. JAMA Pediatr. Aug 22 2022;doi:10.1001/jamapediatrics.2022.2800
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapediatrics.2022.2800&link_type=DOI) 

13. 13.Adler L, Israel M, Yehoshua I, et al. Long COVID symptoms in Israeli children with and without a history of SARS-CoV-2 infection: a cross-sectional study. BMJ Open. Feb 21 2023;13(2):e064155. doi:10.1136/bmjopen-2022-064155
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiMTMvMi9lMDY0MTU1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMjgvMjAyMy4wOS4yNy4yMzI5NjEwMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Funk AL, Kuppermann N, Florin TA, et al. Post-COVID-19 Conditions Among Children 90 Days After SARS-CoV-2 Infection. JAMA Netw Open. Jul 1 2022;5(7):e2223253. doi:10.1001/jamanetworkopen.2022.23253
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2022.23253&link_type=DOI) 

15. 15.Katika A, Zoulias E, Koufi V, Malamateniou F. Mining Greek Tweets on Long COVID Using Sentiment Analysis and Topic Modeling. Stud Health Technol Inform. Jun 29 2023;305:545–548. doi:10.3233/SHTI230554
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3233/SHTI230554&link_type=DOI) 

16. 16.Kasi SG, Dhir SK, Shah A, et al. Coronavirus Disease 2019 (COVID-19) Vaccination for Children: Position Statement of Indian Academy of Pediatrics Advisory Committee on Vaccination and Immunization Practices. Indian Pediatr. Jan 15 2022;59(1):51–57. doi:10.1007/s13312-022-2421-9
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s13312-022-2421-9&link_type=DOI) 

17. 17.Barie PS, Brindle ME, Khadaroo RG, Klassen TL,  Huston JM. Omicron, Long-COVID, and the Safety of Elective Surgery for Adults and Children: Joint Guidance from the Therapeutics and Guidelines Committee of the Surgical Infection Society and the Surgery Strategic Clinical Network, Alberta Health Services. Surg Infect (Larchmt*)*. Feb 2023;24(1):6–18 doi:10.1089/sur.2022.274
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/sur.2022.274&link_type=DOI) 

18. 18.Paul S, Mishra CM. Do we need to vaccinate every child against COVID-19: What evidence suggests-A systematic review of opinions. Front Public Health. 2022;10:1002992. doi:10.3389/fpubh.2022.1002992
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpubh.2022.1002992&link_type=DOI) 

19. 19.Baden LR, El Sahly HM, Essink B, et al. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. N Engl J Med. Feb 4 2021;384(5):403–416. doi:10.1056/NEJMoa2035389
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2035389&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33378609&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

20. 20. El Sahly HM, Baden LR, Essink B, et al. Efficacy of the mRNA-1273 SARS-CoV-2 Vaccine at Completion of Blinded Phase. N Engl J Med. Nov 4 2021;385(19):1774–1785. doi:10.1056/NEJMoa2113017
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2113017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34551225&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

21. 21.Thomas SJ, Moreira ED, Jr., Kitchin N, et al. Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine through 6 Months. N Engl J Med. Nov 4 2021;385(19):1761–1773. doi:10.1056/NEJMoa2110345
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2110345&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34525277&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

22. 22.Elliott P, Haw D, Wang H, et al. Exponential growth, high prevalence of SARS-CoV-2, and vaccine effectiveness associated with the Delta variant. Science. Dec 17 2021;374(6574):eabl9551. doi:10.1126/science.abl9551
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1126/science.abl9551&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34726481&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

23. 23.Gray G, Collie S, Goga A, et al. Effectiveness of Ad26.COV2.S and BNT162b2 Vaccines against Omicron Variant in South Africa. N Engl J Med. Jun 9 2022;386(23):2243–2245. doi:10.1056/NEJMc2202061
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2202061&link_type=DOI) 

24. 24.Tseng HF, Ackerson BK, Luo Y, et al. Effectiveness of mRNA-1273 against SARS-CoV-2 Omicron and Delta variants. Nat Med. May 2022;28(5):1063–1071. doi:10.1038/s41591-022-01753-y
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-022-01753-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35189624&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

25. 25.Al-Aly Z, Bowe B, Xie Y. Long COVID after breakthrough SARS-CoV-2 infection. Nat Med. Jul 2022;28(7):1461–1467. doi:10.1038/s41591-022-01840-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-022-01840-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35614233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

26. 26.Brannock MD, Chew RF, Preiss AJ, et al. Long COVID risk and pre-COVID vaccination in an EHR-based cohort study from the RECOVER program. Nat Commun. May 22 2023;14(1):2914. doi:10.1038/s41467-023-38388-7
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-023-38388-7&link_type=DOI) 

27. 27.Taquet M, Dercon Q, Harrison PJ. Six-month sequelae of post-vaccination SARS-CoV-2 infection: A retrospective cohort study of 10,024 breakthrough infections. Brain Behav Immun. Jul 2022;103:154–162. doi:10.1016/j.bbi.2022.04.013
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2022.04.013&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35447302&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

28. 28.Ayoubkhani D, Bermingham C, Pouwels KB, et al. Trajectory of long covid symptoms after covid-19 vaccination: community based cohort study. BMJ. May 18 2022;377:e069676. doi:10.1136/bmj-2021-069676
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjIwOiIzNzcvbWF5MThfMTEvZTA2OTY3NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzI4LzIwMjMuMDkuMjcuMjMyOTYxMDAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

29. 29.Notarte KI, Catahay JA, Velasco JV, et al. Impact of COVID-19 vaccination on the risk of developing long-COVID and on existing long-COVID symptoms: A systematic review. EClinicalMedicine. Nov 2022;53:101624. doi:10.1016/j.eclinm.2022.101624
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.eclinm.2022.101624&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=36051247&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

30. 30.Arjun MC, Singh AK, Pal D, et al. Characteristics and predictors of Long COVID among diagnosed cases of COVID-19. PLoS One. 2022;17(12):e0278825. doi:10.1371/journal.pone.0278825
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0278825&link_type=DOI) 

31. 31.Wynberg E, Han AX, Boyd A, et al. The effect of SARS-CoV-2 vaccination on post-acute sequelae of COVID-19 (PASC): A prospective cohort study. Vaccine. Jul 30 2022;40(32):4424–4431. doi:10.1016/j.vaccine.2022.05.090
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2022.05.090&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

32. 32.Lorman V, Rao S, Jhaveri R, et al. Understanding pediatric long COVID using a tree-based scan statistic approach: an EHR-based cohort study from the RECOVER Program. JAMIA Open. Apr 2023;6(1):ooad016. doi:10.1093/jamiaopen/ooad016
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jamiaopen/ooad016&link_type=DOI) 

33. 33.Sahai H, Khurshid A. Statistics in epidemiology: methods, techniques, and applications. CRC Press; 1996:321 p.
    
    

34. 34.Ho D, Imai K, King G, Stuart EA. MatchIt: Nonparametric Preprocessing for Parametric Causal Inference. Journal of Statistical Software. 06/14 2011;42(8):1–28. doi:10.18637/jss.v042.i08
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v042.i08&link_type=DOI) 

35. 35.Ridgeway GM, D. F.; Morral, A. R.; Cefalu, M.; Burgette, L. F.; Pane, J. D.; Griffin, B. A. Toolkit for Weighting and Analysis of Nonequivalent Groups: A Tutorial for the R TWANG Package. 2022.
    
    

36. 36.Therneau TM. A Package for Survival Analysis in R. Accessed 2023-08-07, 2023.
    
    

37. 37.Centers for Diseease Control and Prevention Data (CDC). Definitions for COVID-19 Vaccinations in the United States. Accessed August 25, 2022. [https://www.cdc.gov/coronavirus/2019-ncov/vaccines/reporting-vaccinations.html#anchor_1649867074684](https://www.cdc.gov/coronavirus/2019-ncov/vaccines/reporting-vaccinations.html#anchor_1649867074684)
    
    

38. 38.Fowlkes AL, Yoon SK, Lutrick K, et al. Effectiveness of 2-Dose BNT162b2 (Pfizer BioNTech) mRNA Vaccine in Preventing SARS-CoV-2 Infection Among Children Aged 5-11 Years and Adolescents Aged 12-15 Years – PROTECT Cohort, July 2021-February 2022. MMWR Morb Mortal Wkly Rep. Mar 18 2022;71(11):422–428. doi:10.15585/mmwr.mm7111e1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm7111e1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

39. 39.Yewdell JW. Individuals cannot rely on COVID-19 herd immunity: Durable immunity to viral disease is limited to viruses with obligate viremic spread. PLoS Pathog. Apr 2021;17(4):e1009509. doi:10.1371/journal.ppat.1009509
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1009509&link_type=DOI) 

40. 40.Klein NP, Demarco M, Fleming-Dutra KE, et al. Effectiveness of BNT162b2 COVID-19 Vaccination in Children and Adolescents. Pediatrics. May 1 2023;151(5)doi:10.1542/peds.2022-060894
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1542/peds.2022-060894&link_type=DOI) 

41. 41.Price AM, Olson SM, Newhams MM, et al. BNT162b2 Protection against the Omicron Variant in Children and Adolescents. N Engl J Med. May 19 2022;386(20):1899–1909. doi:10.1056/NEJMoa2202826
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2202826&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35353976&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

42. 42.Tan SHX, Cook AR, Heng D, Ong B, Lye DC, Tan KB. Effectiveness of BNT162b2 Vaccine against Omicron in Children 5 to 11 Years of Age. N Engl J Med. Aug 11 2022;387(6):525–532. doi:10.1056/NEJMoa2203209
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2203209&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35857701&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

43. 43.Lipsitch M, Krammer F, Regev-Yochay G, Lustig Y, Balicer RD. SARS-CoV-2 breakthrough infections in vaccinated individuals: measurement, causes and impact. Nat Rev Immunol. Jan 2022;22(1):57–65. doi:10.1038/s41577-021-00662-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41577-021-00662-4&link_type=DOI) 

44. 44.Sullivan SG, Tchetgen Tchetgen EJ, Cowling BJ. Theoretical Basis of the Test-Negative Study Design for Assessment of Influenza Vaccine Effectiveness. Am J Epidemiol. Sep 1 2016;184(5):345–53. doi:10.1093/aje/kww064
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kww064&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27587721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

45. 45.Lipsitch M, Jha A, Simonsen L. Observational studies and the difficult quest for causality: lessons from vaccine effectiveness and impact studies. Int J Epidemiol. Dec 1 2016;45(6):2060–2074. doi:10.1093/ije/dyw124
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw124&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27453361&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F28%2F2023.09.27.23296100.atom) 

46. 46.RECOVER: Researching COVID to Enhance Recovery. (n.d.). RECOVER: Researching COVID to Enhance Recovery. [https://recovercovid.org/](https://recovercovid.org/)