COVID-19 Vaccination and Incidence of Pediatric SARS-CoV-2 Infection and Hospitalization – JAMA Network
							April 26, 2024
							Key Points  
    Question Was implementation of the    pediatric COVID-19 immunization program of California    associated with reductions in the reported pediatric COVID-19    incidence and hospitalizations?  
    Finding In this case series    including 3.9 million children, pediatric vaccination was    estimated to avert 146210 cases of COVID-19 among adolescents    aged 12 to 15 years during a 141-day postvaccine evaluation    period and 230134 cases among children aged 5 to 11 years    during a 199-day postvaccine evaluation period. In addition, an    estimated 168 hospitalizations were averted among children aged    6 to 59 months during a 225-day evaluation period.  
    Meaning The findings of this study    suggest that vaccination against SARS-CoV-2 was associated with    significant reductions in COVID-19 incidence and    hospitalizations among children in California.  
    Importance A SARS-CoV-2 vaccine    was approved for adolescents aged 12 to 15 years on May 10,    2021, with approval for younger age groups following    thereafter. The population level impact of the pediatric    COVID-19 vaccination program has not yet been    established.  
    Objective To identify whether    California's pediatric COVID-19 immunization program was    associated with changes in pediatric COVID-19 incidence and    hospitalizations.  
    Design, Setting, and Participants    A case series on COVID-19 vaccination including children    aged 6 months to 15 years was conducted in California. Data    were obtained on COVID-19 cases in California between April 1,    2020, and February 27, 2023.  
    Exposure Postvaccination    evaluation periods spanned 141 days (June 10 to October 29,    2021) for adolescents aged 12 to 15 years, 199 days (November    29, 2021, to June 17, 2022) for children aged 5 to 11 years,    and 225 days (July 17, 2022, to February 27, 2023) for those    aged 6 to 59 months. During these periods, statewide vaccine    coverage reached 53.5% among adolescents aged 12 to 15 years,    34.8% among children aged 5 to 11 years, and 7.9% among those    aged 6 to 59 months.  
    Main Outcomes and Measures    Age-stepped implementation of COVID-19 vaccination was    used to compare observed county-level incidence and    hospitalization rates during periods when each age group became    vaccine eligible to counterfactual rates predicted from    observations among other age groups. COVID-19 case and    hospitalization data were obtained from the California    reportable disease surveillance system.  
    Results Between April 1, 2020, and    February 27, 2023, a total of 3913063 pediatric COVID-19    cases and 12740 hospitalizations were reported in California.    Reductions of 146210 cases (95% prediction interval [PI],    136056-158948) were estimated among adolescents aged 12 to 15    years, corresponding to a 37.1% (35.5%-39.1%) reduction from    counterfactual predictions. Reductions of 230134    (200170-265149) cases were estimated among children aged 5 to    11 years, corresponding to a 23.7% (20.6%-27.3%) reduction from    counterfactual predictions. No evidence of reductions in    COVID-19 cases statewide were found among children aged 6 to 59    months (estimated averted cases, 259; 95% PI, 1938 to 1019),    although low transmission during the evaluation period may have    limited the ability to do so. An estimated 168 hospitalizations    (95% PI, 42-324) were averted among children aged 6 to 59    months, corresponding to a 24.4% (95% PI, 6.1%-47.1%)    reduction. In meta-analyses, county-level vaccination coverage    was associated with averted cases for all age groups. Despite    low vaccination coverage, pediatric COVID-19 immunization in    California averted 376 085 (95% PI, 348355-417328) reported    cases and 273 (95% PI, 77-605) hospitalizations among children    aged 6 months to 15 years over approximately 4 to 7 months    following vaccination availability.  
    Conclusions and Relevance The    findings of this case series analysis of 3913063 cases    suggest reduced pediatric SARS-CoV-2 transmission following    immunization. These results support the use of COVID-19    vaccines to reduce COVID-19 incidence and hospitalization in    pediatric populations.  
    Vaccination is among the most important interventions to reduce    the public health impact of infectious diseases.1 SARS-CoV-2    mRNA vaccines, including mRNA1273 (Moderna) and BNT162b2    (Pfizer BioNTech), were approved for adult use in December    2020.2 On May 10, 2021, the first    mRNA COVID-19 vaccine was approved for use in adolescents aged    12 to 15 years. Vaccines were subsequently approved for    children aged 5 to 11 years on October 29, 2021, and for    children aged 6 to 59 months on June 17, 2022 (Figure    1).2  
    COVID-19 vaccines are safe for children.3 However,    concerns over vaccine-related adverse events, lower vaccine    effectiveness against illness in children, and perceptions of a    milder disease course in children have resulted in high rates    of parental vaccine hesitancy4-6 and resistance to pediatric    vaccine mandates.7 While California has among    the highest rates of vaccination in the US,8 pediatric    vaccination coverage lags that of adults substantially, with    only 8.2% of children younger than 5 years and 37.8% of    children aged 5 to 11 years fully vaccinated as of May    2023.8 Severe manifestations of    COVID-19 are rare among children, but can occur.9 Understanding    the population-level impact of COVID-19 vaccinations in    SARS-CoV-2 infections and hospitalizations in pediatric    populations would aid in public health decision-making on    pediatric vaccine and booster policy and provide    pediatric-specific information on vaccine outcomes that could    be applied to future SARS-CoV-2 variants.  
    Herein, we analyze data on 3913063 pediatric cases of    COVID-19 and 12740 hospitalizations from California. Using the    phased introduction of the vaccine to individuals aged 12 to 15    years, 5 to 11 years, and 6 to 59 months, we estimated    statewide and county-specific outcomes associated with    vaccination on pediatric incidence and hospitalizations in    California.  
    We obtained deidentified information on all pediatric COVID-19    cases reported in California between April 1, 2020, and    February 27, 2023, from the California COVID-19 Reporting    System, along with the patients county of residence, age, and    hospitalization status. Each case was confirmed using a nucleic    acid amplification test. Because the research constitutes a    public health surveillance activity, the study did not    constitute human research and does not require institutional    review board review or exemption according to the Common Rule    (45 CFR 46). We followed the reporting    guideline for case series studies.  
    Daily cases were aggregated by county and age groups based on    dates of vaccination eligibility: 0 to 5 months (vaccine    ineligible), 6 to 59 months, 5 to 11 years, 12 to 15 years, and    older than 16 years (nonpediatric). To remove variation due to    differential health care seeking by day of week, we calculated    7-day moving averages of case counts per county and age group.    Due to small counts for pediatric hospitalizations, we    aggregated hospitalizations by week and age group within 5    California-designated regions (eFigure 1 in Supplement    1). Descriptions of other covariate data are covered in the    eMethods in Supplement 1).  
      Training and Prediction Periods    
    For each age group of interest (6-59 months, 5-11 years, and    12-15 years), we split data into age-eligible and    age-ineligible periods. The prevaccine eligibility period    encompassed data up to the date of vaccine eligibility. The    evaluation period for the outcome associated with the vaccine    lasted from 1 month following the date of vaccine eligibility    (to allow time to complete 2 doses) until the date that the    next age group became eligible or until the time of analysis    (February 27, 2023) (Figure 1). Ending the evaluation    period at the vaccine eligibility date of the next eligible age    group permitted that age group to be selected as a control time    series in our predictive models.  
      Candidate Model Generation and Selection    
    We developed a set of candidate predictive generalized linear    models, which are described in greater detail in the eMethods    in Supplement 1. Predictors eligible for    selection within candidate models included (1) log-incidence    series for other age groups (ie, <6 months, 6-59 months,    5-11 years, 12-15 years, and 16 years, omitting the group    being modeled) included as either covariates or as an offset    term for any 1 age group; (2) an indicator of vaccine age    eligibility for other age groups; (3) an indicator for    in-person school being in session; (4) interactions between    school and vaccine introduction indicators and time series for    other age groups, aiming to account for differences in constant    proportionality during school periods or when 1 age group    became vaccinated; and (5) seasonal controls. Eligible    predictors are summarized in eTable 1 in Supplement    1. Quasi-Poisson distributions were fit for the outcome to    account for overdispersion. Models were developed separately    for each California county.  
    Candidate models for hospitalizations included similar eligible    predictors, with 2 main differences: log weekly case incidence    series for other age groups was lagged by 2 weeks in accordance    with the expected lag between infection and    hospitalization10 and unlagged weekly    hospitalizations across other age groups were included as    possible predictors. Models were developed separately for each    of the 5 regions.  
    We used a time series with a 1-year gap cross-validation    approach (eMethods and eFigure 2 in Supplement    1) to select the best predictive model for each age group    and geographic area (county or region) within the prevaccine    period.11-13 For each area-age group    combination, we selected the model with the lowest    out-of-sample mean square error across holdout folds. For this    model, we also calculated the coefficient of determination,    r2, a goodness-of-fit metric. The selected    predictors varied by area and age group. Selected models for    each area-age group combination are included in eTable 2 (for    cases) and eTable 3 (for hospitalizations) in Supplement    1.  
      Calculation of Vaccine Outcomes and Association Between      Averted Cases and Vaccination Coverage    
    Selected models were fit to prevaccine eligibility data for    their age group and geographic area and then used to predict    counterfactual incidence or hospitalization in the postvaccine    period or the expected case or hospitalization counts had    vaccination not occurred. For inference, we computed 95%    prediction intervals (PIs) around the counterfactual    predictions, using a sandwich estimator to account for    overdispersion when computing SEs (eMethods in Supplement    1).14 Prediction intervals, which    are wider than CIs, capture the uncertainty around each future    predicted value. Statewide estimates were obtained by summing    predictions across geographies (eMethods in Supplement    1).  
    We estimated the absolute and relative differences between    predicted counterfactual values and observed values for each    county or region during the postvaccine evaluation period. To    understand the association between vaccination coverage and    averted cases, we fit regression models relating the reduction    in cases within each age and county to county-level vaccination    coverage within the same age group, using a fixed-effects    meta-analysis with weights equal to the inverse estimated SE of    the estimates per county. We used segmented regression models    (eMethods in Supplement 1) to examine whether there    were coverages below which reductions in cases could not be    identified or above which diminishing returns on vaccination    were observed.15  
    To examine whether postvaccine predictions from a different,    but well predictive model, yielded similar estimates of    vaccination outcomes, we repeated model selection using the    mean absolute error instead of the mean square error in our    cross-validation algorithm. We conducted jackknife analyses to    examine whether postvaccine predictions from any one county    were driving observed effects, dropping each county in turn    from the overall pool of counties and recalculating the primary    analytic end point of cases averted.  
    All analyses were conducted in R, version 3.6.0 (R Foundation    for Statistical Computing).16  
    Between April 1, 2020, and February 27, 2023, a total of    3913063 COVID-19 cases were reported in California among    individuals aged 18 years or younger. Of these, 47174 cases    (1.2%) were among children younger than 6 months, 517447    (13.2%) in children aged 6 to 59 months, 1590806 (40.7%) in    children aged 5 to 11 years, and 1511690 (38.6%) in    adolescents aged 12 to 15 years. A total of 12740    hospitalizations were reported: 1443 (11.3%) were among    children younger than 6 months, 3428 (26.9%) in children aged 6    to 59 months, 2536 (19.9%) in children aged 5 to 11 years, and    3921 (30.8%) in adolescents aged 12 to 15 years.  
      Vaccine-Attributable Averted Cases and Hospitalizations by      Pediatric Age Group    
    As shown in eFigure 3 in Supplement    1, r2 values for models fit to daily case    data were 0.92 (IQR, 0.79-0.96) for children aged 6 to 59    months, 0.89 (IQR, 0.78-0.95) for children aged 5 to 11 years,    and 0.79 (IQR, 0.62-0.90) for adolescents aged 12 to 15 years.    eFigure 4 in Supplement 1 shows the model fit for    hospitalizations. More details on model fit are included in the    eResults in Supplement 1.  
      Adolescents Aged 12 to 15 Years    
    Individuals aged 12 to 15 years were eligible to be vaccinated    against SARS-CoV-2 as of May 10, 2021. By October 29, 2021,    when the next age group became eligible, 53.5% of this    population had completed the 2-dose primary series of the    vaccine, corresponding to 1712686 individuals. County-level    vaccination rates ranged from 11.5% to 85.7%.8 During the    141 days spanning June 10 to October 29, 2021, 247700 COVID-19    cases were observed among individuals aged 12 to 15 years. We    estimated that 394506 (95% PI, 392545-396467) cases of    COVID-19 would have occurred absent vaccination, corresponding    to 146210 (95% PI, 136056-158948) cases averted statewide or    37.1% (95% PI, 34.5%-40.3%) of expected cases (Table, Figure 2C).    Incidence plots from all counties are included as eFigures 9-14    in Supplement 1.  
    During this same 141-day period, 688 hospitalizations were    observed among adolescents. We estimated that 59 (95% PI, 65    to 244) hospitalizations were averted or a reduction of 7.9%    (95% PI, 8.7% to 32.7%) from expectation (Table, Figure 3C).    Hospitalization plots from all regions are included as eFigures    12-14 in Supplement 1.  
      Children Aged 5 to 11 Years    
    Children aged 5 to 11 years were eligible for vaccination on    October 29, 2021. By June 17, 2022, 1219432 individuals    (34.8% of this population) had completed a primary series of    the vaccine, with a range of 10.0% to 74.7% by    county.8 During the 199-day period    following November 29, 2021, we estimated that 230134 (95% PI,    200170-265149) cases were averted due to the vaccine    corresponding to a reduction of 23.7% (95% PI, 20.6%-27.3%)    from counterfactual expectations (Table, Figure 2B). During this same    period, we estimated that 46 (95% PI, 79 to 221)    hospitalizations were averted, corresponding to 5.8% (95% PI,    10.2% to 28.6%) of expected hospitalizations (Table, Figure 3B).  
      Children Aged 6 to 59 Months    
    Children aged 6 to 59 months were eligible for vaccination on    June 17, 2022. By February 27, 2023, 177087 (7.9%) individuals    had received both doses of the primary series, with a range of    0.7% to 38.5% across counties.8 In the 225 days following    July 17, 2022, we did not detect any significant changes in    cases from counterfactual expectations in the postvaccine    period (estimated averted cases: 259; 95% PI, 1938  1019)    (Table). The    postvaccine evaluation period for this age group did not    include a surge in COVID-19 cases as it did for the other age    groups (Figure 2A). However, we estimated    that 168 (95% PI, 42-324) hospitalizations were averted    following vaccination, or a reduction of 24.4% (95% PI,    6.1%-47.1%) from counterfactual expectations (Table and Figure 3A). Summing    across all age groups, we estimated that pediatric vaccination    was associated with reductions of 376 085 (95% PI, 348 355-417    328) reported cases and 273 (95% PI, 77-605) hospitalizations    among children aged 6 months to 15 years during the 4 to 7    months following vaccine availability. This represents a    reduction of 26.3% of the number of cases and 12.4% of the    hospitalizations that would have been seen in this population    absent the vaccine.  
    As indicated in the eResults and eFigures 5 and 6 in Supplement    1, results for individuals aged 5 to 15 years were not    sensitive to the inclusion of any single county, although    results for children aged 6 to 59 months were sensitive to the    inclusion of Los Angeles (eFigure 7 in Supplement    1). Estimated cases (eTable 4 in Supplement    1) and hospitalizations (eTable 5 in Supplement    1) were consistent when model selection was done using mean    absolute error as the loss function for children aged 5 to 15    years. Estimated averted cases in children aged 6 to 59 months    were slightly lower using mean absolute error, but    hospitalization results were consistent (eResults in Supplement    1). Estimates of cases averted (eTable 6 in Supplement    1) and hospitalizations averted (eTable 7 in Supplement    1) made using the mean absolute error as the loss function    for each county or region are available, along with plots of    observed and counterfactual case and hospitalization series for    all geographic areas (eFigures 9-14 in Supplement    1).  
      Association Between Averted Cases and Vaccination    
    County-level vaccination coverage explained 26% of variation of    cases averted for children aged 6 to 59 months, 28% for    children aged 5 and 11 years, and 12% for adolescents aged 12    to 15 years (Figure 4). On average, every    increase of 10 vaccinations per 1000 children corresponded to a    reduction of 0.9 (95% CI, 0.3-1.4) cases per 1000 children for    individuals aged 6 to 59 months, 3.5 (95% CI, 1.9-5.1) cases    per 1000 children for those aged 5 and 11 years, and 2.0 (95%    CI, 0.6-3.4) cases per 1000 children for adolescents aged 12 to    15 years. Linear model fits had lower Akaike information    criterion and bayesian information criterion values than    segmented regression model fits for all age groups. Across all    age groups, pediatric vaccination rates in California were    generally highest among Bay Area counties (eFigure 8 in    Supplement    1), which also ranked highest for averted cases due to    vaccination (eResults in Supplement    1).  
    We provide evidence that Californias pediatric COVID-19    immunization program averted 376 085 (95% PI, 348355-417328)    reported cases and 273 (95% PI, 77-605) hospitalizations among    children aged 6 months to 15 years during the 4 to 7 months    following vaccine availability. This represents a reduction of    26.3% of the number of cases that would have been seen in this    population absent the vaccine. Prior work has similarly    reported a high impact of widespread administration of mRNA    vaccines in adult populations. In California, COVID-19 vaccines    were estimated to avert more than 1.5 million cases, 72 000    hospitalizations, and 19 000 deaths statewide during the first    10 months of vaccination (through October 16,    2021).17 In the US, each 10% increase    in vaccination coverage among individuals aged 18 years or    older at the county level was associated with an 8% reduction    in mortality and a 7% reduction in incidence.18 Similarly, a    study in Israel estimated that nearly 650 000 cases of COVID-19    were averted in the first 2 months following vaccination    introduction,19  
    Earlier studies have estimated vaccine effectiveness in    pediatric populations by comparing incidence rates among    vaccinated children with those in unvaccinated children using    test-negative designs,20-22 or    retrospective23,24 or prospective cohort    studies.25 Our counterfactual case    series approach, which has been used in other studies to    estimate the population-level impact of interventions with a    clearly specified rollout time,26,27 enables calculation of    vaccine program impact at the population level, without    information on individual vaccine status.  
    The cumulative effect of vaccination at the population level    may be meaningful even if individual vaccine effectiveness is    low. While influenza vaccine effectiveness was estimated at 29%    in 2017-2018,28 it was estimated that    widespread vaccination averted more than 3.1 million cases of    influenza in the US.29 Nevertheless, overall impact    depends on vaccine coverage. We identified positive    associations between county-level vaccination coverage and    averted cases in each age group, whereby each 10 additional    vaccinations per 1000 children corresponded to an average    reduction of 0.9 to 3.5 cases per 1000 children. Segmented    regression models associating vaccine coverage with averted    cases did not identify break points, suggesting that over the    range of vaccination coverages examined (0%-85%), we saw    neither diminishing returns on increased coverage owing to the    acquisition of sufficient population-level immunity nor a    threshold below which vaccination has limited public health    impact. This is consistent with the persistence of SARS-CoV-2    circulation in populations with high vaccination coverage and    resulting value of direct protection.  
    Results for individuals aged 6 to 59 months differed from those    of older age groups in that we found a significant reduction in    hospitalizations, but not cases, following vaccination. One    explanation for this discrepancy could be that postvaccine    evaluation period for children aged 6 to 59 months did not    include a surge in COVID-19 cases as it did for the other age    groups (Figure 2), potentially making it    difficult to detect statistically significant reductions from    the counterfactual. However, vaccine effectiveness of early    mRNA vaccines was lower against Omicron variants compared with    Alpha and Delta variants,30,31 and the Omicron variant    dominated during the postvaccine period for children aged 6    months to 11 years (Figure 1). The detection of    significant reductions in hospitalization in this age group,    but not others, may be due, in part, to the fact that COVID-19    mortality disproportionately affects very young children    compared with older children.32 For older age groups, we    also estimated reductions in hospitalizations, although the 95%    PI spans 0. However, we note that 95% CIs are narrower than PIs    and may not have encompassed the null.  
    This study has limitations. Case data represented individuals    who sought testing, which may be differential across    unvaccinated and vaccinated groups, geographies, and time.    Access to at-home testing likely resulted in further case    underascertainment. If individuals were, on average, less    likely to seek care for mild illness following vaccination, our    analysis could have overestimated the absolute effect of the    vaccine on cases averted. Overestimation of the relative effect    of the vaccine may have resulted if vaccine recipients were    disproportionately represented in the surveillance record both    before and after vaccine eligibility compared with    never-vaccinated individuals being more connected to care. Data    on hospitalizations are less likely to be subjected to biases    from differential case ascertainment. We estimated significant    reductions in hospitalizations following vaccine introduction    compared with counterfactual predictions.  
    Several considerations could lead to underestimates of the    association between vaccination and child long-term health.    First, asymptomatic cases are less likely to be reported, yet    remain an important outcome, as postCOVID-19 condition    symptoms may present after asymptomatic    infections.33-35 Second, we were unable to    estimate indirect outcomes associated with the vaccine in other    age groups or control for social contacts. If children    increased social contacts following receipt of the vaccine, as    has been shown elsewhere,36 they may be challenged more    frequently with SARS-CoV-2. Third, we assessed the outcomes of    the vaccine over a short postvaccination period, limiting our    ability to examine vaccine responses under waning immunity.  
    Two important limitations relate to model functional form.    First, attributing differences between the observed and the    predicted counterfactual cases to the vaccine assumes that the    associations between incidence in the age group being modeled    and incidence in the age groups selected as model predictors    would, absent the vaccine, remain constant over the pre-to-post    vaccine periods. This would not occur if one age group    developed increased immunity or if different variants had    differential age-disease associations. This is especially    salient for the 5- to 11-year age group, as the models were    primarily trained on data from the period when the Delta    variant predominated, yet the Omicron variant, which is less    reliant on angiotensin-converting enzyme 2 binding for    entry37 and disproportionately    influenced children younger than 5 years, prevailed in the    evaluation period. Accordingly, the effect of vaccination may    have been overestimated for this age group in counties where    the incidence in children younger than 5 years was selected as    a predictor (eTable 2 in the Supplement).  
    Second, there is potential for unstable predictions in the    evaluation period if the predictive model was faced with values    of selected predictors that fell outside the range of data used    to fit the model. Our time series with a gap cross-validation    approach guards against both of these limitations by    prioritizing selection of generalized linear models that do    well predicting values in periods that follow the training    period, and in periods where the predictors may fall outside    the range of what they were during the training    period.11-13,38 Moreover, generalized linear    models selected using different loss functions resulted in    similar model predictions during the postevaluation period,    suggesting that results are robust to differences in the nature    of the association between incidence in the modeled age group    and incidence in the predictor age groups.  
    In this case series analysis of 3913063 pediatric cases, we    provide evidence suggesting that programmatic vaccination    against SARS-CoV-2 was associated with significant reductions    in COVID-19 incidence among children in California in the 4 to    7 months following vaccine eligibility. At the county level, we    found associations of higher vaccine coverage with greater    reductions in pediatric cases. Our results support the use of    COVID-19 vaccines to reduce COVID-19 incidence and    hospitalization in pediatric populations.  
    Accepted for Publication: February 23, 2024.  
    Published: April 23, 2024.    doi:10.1001/jamanetworkopen.2024.7822  
    Open Access: This is an open access article    distributed under the terms of the CC-BY    License.  2024 Head JR et al. JAMA Network Open.  
    Corresponding Author: Justin V. Remais, PhD,    2121 Berkeley Way, #5302, Berkeley, CA 94720 (jvr@berkeley.edu).  
    Author Contributions: Dr Head had full access    to all of the data in the study and takes responsibility for    the integrity of the data and the accuracy of the data    analysis.  
    Concept and design: Head, Len, Lewnard, Remais.  
    Acquisition, analysis, or interpretation of data: Head,    Collender, Len, White, Sud, Camponuri, Lee, Remais.  
    Drafting of the manuscript: Head, Remais.  
    Critical review of the manuscript for important intellectual    content: All authors.  
    Statistical analysis: Head, Collender, Len, Camponuri,    Lee.  
    Obtained funding: Head, Remais.  
    Administrative, technical, or material support:    Collender, Len, White, Sud, Camponuri, Remais.  
    Supervision: Remais.  
    Conflict of Interest Disclosures: Dr Len    reported having been an employee of the California Department    of Public Health (CDPH). No other disclosures were reported.  
    Funding/Support: This project was supported by    a grant from the CDPH through the University of California    Health & CDPH COVID Modeling Consortium. Dr Head was supported    by the National Institute of Allergy and Infectious Diseases    (NIAID) National Institutes of Health (NIH) award K01AI173529.    Dr Remais was supported by NIAID NIH award R01AI148336.  
    Role of the Funder/Sponsor: The funding    organizations did not play a role in the design and conduct of    the study; collection, management, analysis, and interpretation    of the data; preparation, review, or approval of the    manuscript; and decision to submit the manuscript for    publication.  
    Data Sharing Statement: See Supplement    2.  
Visit link: 
COVID-19 Vaccination and Incidence of Pediatric SARS-CoV-2 Infection and Hospitalization - JAMA Network