Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 30, Number 6—June 2024
Research

SARS-CoV-2 Disease Severity and Cycle Threshold Values in Children Infected during Pre-Delta, Delta, and Omicron Periods, Colorado, USA, 2021–2022

On This Page
Methods
Results
Discussion
Cite This Article
Figures
Figure 1
Figure 2
Figure 3
Figure 4
Tables
Table 1
Table 2
Table 3
Downloads
Article
Appendix
RIS [TXT - 2 KB]
Article Metrics
Metric Details
Related Articles
S. pyogenes Disease, Germany, 2005–2023
Facial Palsy after COVID-19 Vaccination
SARS-CoV-2 RBD and Neutralizing Antibody Titers
More articles on Coronavirus, COVID-19
Author affiliations: Colorado Department of Public Health and Environment, Denver, Colorado, USA (L. Bankers, S.C. O’Brien, D.M. Tapay, I. Armistead, A. Burakoff, S.R. Matzinger); University of Colorado, Aurora, Colorado, USA (E. Ho, S.R. Dominguez); Children’s Hospital Colorado, Aurora (E. Ho, S.R. Dominguez)

Cite This Article

Abstract

In adults, viral load and disease severity can differ by SARS-CoV-2 variant, patterns less understood in children. We evaluated symptomatology, cycle threshold (Ct) values, and SARS-CoV-2 variants among 2,299 pediatric SARS-CoV-2 patients (0–21 years of age) in Colorado, USA, to determine whether children infected with Delta or Omicron had different symptom severity or Ct values than during earlier variants. Children infected during the Delta and Omicron periods had lower Ct values than those infected during pre-Delta, and children <1 year of age had lower Ct values than older children. Hospitalized symptomatic children had lower Ct values than asymptomatic patients. Compared with pre-Delta, more children infected during Delta and Omicron were symptomatic (75.4% pre-Delta, 95.3% Delta, 99.5% Omicron), admitted to intensive care (18.8% pre-Delta, 39.5% Delta, 22.9% Omicron), or received oxygen support (42.0% pre-Delta, 66.3% Delta, 62.3% Omicron). Our data reinforce the need to include children, especially younger children, in pathogen surveillance efforts.

SARS-CoV-2 is characterized by diverse variants (1) with differing transmissibility and disease severity (2). The rapid evolution and spread of new variants has required a nimble public health response to understand dynamics and clinical implications (3), but most work thus far has focused on adults.

The highly transmissible SARS-CoV-2 Delta variant (B.1.617.2) was detected and became the predominant variant in the United States over an 8-week period during May–June 2021 (4,5). Adults infected by Delta exhibited higher viral loads (6) and potentially increased disease severity than those infected by previous variants (3,5). Taylor et al. (5) found a significantly higher proportion of hospitalized patients after Delta became predominant than in earlier pandemic phases. However, in-hospital outcomes did not differ between Delta and pre-Delta variants. Similarly, Twohig et al. (3) observed higher risk for hospitalization and emergency care for those infected by Delta than by Alpha (B.1.1.7).

During December 2021, Delta was swiftly overtaken by the more transmissible Omicron (B.1.1.529), which became predominant in the United States over a 2-week period (4). Adults infected by Omicron tended to exhibit similar viral loads but lower disease severity compared with Delta-infected adults (7). Ulloa et al. (7) found reduced risk for hospitalization, intensive care unit (ICU) admission, and death among persons infected by Omicron compared with Delta. Taylor et al. (8) observed smaller proportions of hospitalized patients admitted to the ICU or requiring invasive mechanical ventilation and lower rates of in-hospital death with Omicron than Delta. In addition, some studies found that vaccine effectiveness against hospitalization and visits to emergency department or urgent care was lower against Omicron than against Delta (9).

Less is known about symptom severity among pediatric patients infected with different variants. Some work suggests little or no difference in disease severity among children infected by Delta compared with previous variants (10). However, much of that work has been limited to hospitalized patients. In addition, although the incidence rate of detected infections with Omicron was higher for young children compared with Delta, some studies showed that clinical outcomes of infections with Omicron tended to be less severe than Delta (11). However, the large increase in the number of infections during Omicron’s predominance might increase the number of severe outcomes (12).

The relationship between cycle threshold (Ct) values and quantitative viral load is tightly inversely correlated (13), such that relatively lower Ct values are indicative of higher viral loads. Similar to adults, children who have symptomatic SARS-CoV-2 infections might have lower Ct values (suggesting higher viral loads) compared with those who have asymptomatic infections (14), and those infected with Delta or Omicron might have lower Ct values than those infected with other variants (15). However, the interplay between SARS-CoV-2 variant, Ct, and symptom severity has not been well studied among pediatric cases, especially among nonhospitalized children.

We evaluated hospitalized and nonhospitalized SARS-CoV-2 pediatric cases in Colorado, USA, during January 2021–January 2022. We used clinical surveillance data, Ct values, and whole-genome sequencing (WGS) to determine whether children infected with either Delta or early Omicron variants had different symptom severity or Ct values from children infected with earlier SARS-CoV-2 variants.

Methods

Ethics Statement

The data used in this study were generated for public health surveillance purposes. This activity was determined to be consistent with enhanced disease surveillance activities, not human subjects research, by the Colorado Department of Public Health and Environment’s (CDPHE) Communicable Disease branch. Institutional review board approval was provided by the Colorado Multiple Institutional Review Board. Informed consent was waived.

Study Population

The study population consisted of hospitalized (hereafter inpatient) and nonhospitalized (hereafter outpatient) cases from the Children’s Hospital Colorado (CHCO) hospitals and outpatient clinics who were Colorado residents <21 years of age and tested positive for SARS-CoV-2 by PCR during January 1, 2021–January 31, 2022. We collected demographic and clinical information including age, race/ethnicity, symptomatology, number of days between symptom onset and positive test (hereafter symptom onset date), and vaccination status. Data were extracted from a state communicable disease database (Colorado Electronic Disease Reporting System), in which communicable diseases are entered as part of public health surveillance and investigation activities.

We performed inpatient chart abstraction to collect data about comorbidities, admission and discharge dates, whether patients were admitted to the hospital because of COVID-19 as opposed to with COVID-19, whether patients were symptomatic because of COVID-19, whether they were admitted to the pediatric ICU (PICU) because of COVID-19, whether they received oxygen support (and type of support) because of COVID-19, and vaccination status. We limited comparisons among inpatients and between inpatients and outpatients to those hospitalized because of COVID-19. We classified patients as admitted because of COVID-19 if their primary complaint symptoms at time of admission to CHCO were consistent with COVID-19 and a COVID-19 test was positive at admission, or if symptoms consistent with COVID-19 developed during hospitalization (with positive test upon admission or later during hospitalization) that would have resulted in admission if they were not already hospitalized. We classified patients as admitted with COVID-19 if they sought care at CHCO with a non–COVID-19 primary diagnosis (e.g., trauma, psychiatry, social reasons, surgery, non–COVID-19 medical diagnoses) and tested positive upon admission on routine surveillance testing. Those persons might also have had >1 symptoms consistent with COVID-19 that did not require admission solely for those symptoms. We determined classification by reviewing clinical provider documentation of assessments and medical decision-making, because International Classification of Diseases, 10th Revision, codes were not consistently available. We considered all common symptoms of COVID-19, including fever, congestion/rhinorrhea, cough, shortness of breath, vomiting/diarrhea, fatigue, headache, and loss of sense of taste/smell.

RNA Extraction and Quantitative Reverse Transcription PCR

SARS-CoV-2–positive nasopharyngeal swab specimens collected at CHCO hospitals (inpatient) and outpatient clinics (outpatients) were sent to the CDPHE Laboratory. We extracted RNA using the Applied Biosystems MagMAX Viral/Pathogen II Nucleic Acid Isolation Kit on the KingFisher Flex System for automated extraction (ThermoFisher Scientific, https://www.thermofisher.com). We used the Applied Biosystems TaqPath COVID-19 Combo Kit multiplexed reverse transcription PCR (RT-PCR) (ThermoFisher Scientific) to obtain Ct values for open reading frame 1ab (ORF1ab), spike (S), and nucleocapsid (N) gene targets. We used N gene Ct values as a correlate of viral load (13).

Library Preparation and WGS

We performed WGS on either GridION (Oxford Nanopore Technologies, https://www.nanoporetech.com) or NextSeq550 (Illumina, https://www.illumina.com) on all samples that met sequencing criteria (N gene Ct <28). For samples sequenced on GridION, we performed library preparation following the ARTIC tiled PCR amplicon sequencing protocol (16,17), and for samples sequenced on the NextSeq550, we performed library preparation and single-end Illumina sequencing following the Illumina COVID-Seq assay (18) (Appendix).

Bioinformatic Analysis of WGS Data

We performed assembly and analysis of WGS data on the terra.bio platform (19), using our custom, publicly available workflows for SARS-CoV-2 (20,21) (Appendix). Sequenced samples with >50% coverage are publicly available in the GISAID (https://www.gisaid.org) database and the National Center for Biotechnology Information Sequence Read Archive (BioProject accession no. PRJNA686984) (Appendix Table 1).

Data Analysis and Statistics

We merged WGS results and deidentified patient data using internal identifiers in Tableau version 2021.4 (22). To include data for samples that were not successfully sequenced, we performed analyses in 2 ways.

Figure 1

Variant counts of sequencing-confirmed lineages by week of sample collection in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Gray boxes indicate time periods of potentially mixed lineage that were excluded from time period analyses.

Figure 1. Variant counts of sequencing-confirmed lineages by week of sample collection in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Gray...

Our first approach mitigates the risk that Ct value analyses could be biased toward lower Ct values because of ability to sequence. We generated a lineage distribution over the course of the study and assigned all samples with Ct values to variant periods. We defined 3 periods on the basis of when Delta or Omicron were at a prevalence of >80% among sequenced Colorado pediatric samples: pre-Delta during January 1–May 1, 2021; Delta during June 15–September 29, 2021; and Omicron during December 26, 2021–January 31, 2022. We excluded samples collected between the pre-Delta and Delta periods and between the Delta and Omicron periods (in which variants were mixed) to improve the accuracy of variant period assignments (Figure 1). Second, to evaluate sensitivity, we performed analyses only on samples that were successfully sequenced and assigned a lineage.

We visualized analyses of Ct values, variant period and lineage, and patient data using Tableau. We performed statistical analyses (analysis of covariance, 1- or 2-way analysis of variance with Tukey test, or χ2 test, as appropriate) in Rstudio base packages version 1.4.1106 (RStudio, http://www.rstudio.com). We calculated descriptive characteristics of inpatients including counts and proportions with 95% CI using R version 4.1.1 (The R Foundation for Statistical Computing, https://www.r-project.org). We considered a p value <0.05 statistically significant.

Results

Characteristics of Patient Population

The study population included 2,299 persons, 1,629 (70.9%) outpatients and 670 inpatients (29.1%). Of the 670 inpatients, 395 were hospitalized because of COVID-19 and 275 were hospitalized with COVID-19) (Table 1; Appendix Tables 2, 3). Among the entire population, 40.1% were female, 32.6% were non-Hispanic White, 24.4% were Hispanic, and the median age was 6.5 years. A total of 1,724 cases occurred during the pre-Delta (n = 429), Delta (n = 752), and Omicron (n = 543) periods. Proportions of male versus female patients and those with >1 known comorbidity did not differ by variant period (p = 0.44 by χ2 test for sex and p = 0.62 by χ2 test for comorbidities). However, proportions of racial/ethnic groups differed significantly across variant periods (p<0.0001 by χ2 test).

Among inpatients, a larger proportion infected during the Delta or Omicron periods were symptomatic (p<0.0001 by χ2 test), admitted to the PICU (p = 0.002 by χ2 test), or received oxygen support (p = 0.04 by χ2 test) because of COVID-19 compared with persons infected during the pre-Delta period (Table 2; Appendix Tables 4, 5). Mean hospital stay durations were 5.72 (SD +8.30) days for the pre-Delta period, 7.84 (SD +20.6) days for the Delta period, and 3.32 (SD +5.79) days for the Omicron period and were not significantly different (p = 0.1).

We obtained Ct values for 1,796 (78.1%) of 2,299 persons; of those, 362 were inpatients and 1,434 were outpatients. Of those 1,796 persons, 1,240 were collected during the 3 variant periods: pre-Delta period (n = 307), Delta period (n = 582), and Omicron period (n = 351). We successfully sequenced 1,276 (55.5%) of 2,299 samples (219 inpatients and 1,057 outpatients) to >50% coverage and obtained lineage calls for 1,177 (51.2%) samples (183 inpatients and 1,053 outpatients; 654 Delta, 334 Omicron, and 189 other).

Ct Value Patterns across All Variants

Figure 2

Ct value patterns across all variants in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Boxplots indicate overall Ct value patterns across categories of patient characteristics, regardless of variant period. A) Age group; B) vaccination status (unvaccinated vs. vaccinated with any number of doses); C) vaccination status (unvaccinated versus vaccinated by number of doses); D) patient type (outpatient/not hospitalized, hospitalized because of COVID-19, hospitalized but not because of COVID-19); E) symptomatic versus asymptomatic; F) number of days between symptom onset and positive test (symptom onset groups); G) hospitalized but not admitted to PICU versus admitted to PICU; H) any type of supplemental oxygen support versus no oxygen support received; I) highest level of supplemental oxygen support received (none, noninvasive oxygen support, invasive oxygen support). Significance was determined using 1-way analysis of variance with Tukey test. Brackets indicate which comparisons correspond to the significance codes, and connected brackets indicate comparisons that have the same significance code. Ct, cycle threshold; PICU, pediatric intensive care unit.

Figure 2. Ct value patterns across all variants in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Boxplots indicate overall Ct value...

We first evaluated overall patterns of Ct values among our dataset categories. Patients <1 year of age had significantly lower Ct values than the other age groups (adjusted p<0.0001 for 1–4 years, 5–11 years, and >12 years) (Table 3; Figure 2). We observed the same pattern among patients <1 year of age when limited to sequenced samples (adjusted p = 0.0003 for patients 1–4 years of age, adjusted p = 0.0002 for 5–11 years, and adjusted p = 0.0002 for >12 years) (Appendix Tables 6–8, Figure 1). Patients who received any number of vaccine doses had significantly higher Ct values than unvaccinated patients (adjusted p = 0.003). In addition, patients who had received a booster had significantly higher Ct values than unvaccinated patients (adjusted p = 0.0081).

Regarding potential disease severity indicators, we found that inpatients hospitalized with COVID-19 had significantly higher Ct values than both inpatients hospitalized because of COVID-19 (adjusted p<0.0001) and outpatients (adjusted p<0.0001). Ct values for inpatients hospitalized because of COVID-19 and outpatients did not differ significantly from each other (Table 3; Figure 2). Symptomatic patients had significantly lower Ct values than asymptomatic patients (adjusted p = 0.0081). Those tested soon after symptom onset had lower Ct values than patients tested 8–14 days after symptom onset (adjusted p<0.0001 for 0–3 days and 4–7 days after symptom onset). Similarly, Ct values were significantly lower in the groups tested 0–3 days and 4–7 days after symptom onset than in the 8–14 days group when limited to sequenced samples (adjusted p<0.0001 for 0–3 days, adjusted p = 0.0055 for 4–7 days) (Appendix Tables 6–8, Figure 1).

Among inpatients, we did not observe a significant difference in Ct values according to PICU admission status. However, those who received any type of supplemental oxygen support (invasive or noninvasive) and those who received noninvasive supplemental oxygen had lower Ct values than those who did not receive any form of supplemental oxygen (adjusted p = 0.0002 for any type of supplemental oxygen and adjusted p = 0.0012 for noninvasive supplemental oxygen).

Ct Value Patterns by Variant Period

Figure 3

Ct value patterns of patient characteristics across variant periods in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Boxplots indicate overall Ct value patterns among patient characteristics across variant periods. A) Variant period; B) age group; C) vaccination status (unvaccinated vs. vaccinated with any number of doses); D) vaccination status (unvaccinated vs. vaccinated by number of doses). Red = pre-Delta, Blue = Delta, Orange = Omicron. Significance was determined using 1-way analysis of variance (A) or 2-way (B–D) with Tukey test. Brackets indicate which comparisons correspond to the significance codes, and connected brackets indicate comparisons that have the same significance code. Ct, cycle threshold.

Figure 3. Ct value patterns of patient characteristics across variant periods in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Boxplots indicate...

Figure 4

Ct value patterns among markers of disease severity across variant periods in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022. Boxplots indicate overall Ct value patterns across disease severity markers by variant period. A) Patient type (outpatient/not hospitalized, hospitalized because of COVID-19); B) symptomatic versus asymptomatic; C) number of days between symptom onset and positive test (symptom onset groups); D) hospitalized but not admitted to PICU versus admitted to PICU; E) any type of supplemental oxygen support versus no oxygen support received; F) highest level of supplemental oxygen support received (none, noninvasive oxygen support, invasive oxygen support). Red = pre-Delta, Blue = Delta, Orange = Omicron. Significance was determined using 2-way analysis of variance with Tukey test. Brackets indicate which comparisons correspond to the significance codes, and connected brackets indicate comparisons that have the same significance code. Ct, cycle threshold; PICU, pediatric intensive care unit.

Figure 4. Ct value patterns among markers of disease severity across variant periods in study of SARS-CoV-2 disease severity in children during pre-Delta, Delta, and Omicron periods, Colorado, USA, January 2021–January 2022....

We analyzed Ct value patterns for each dataset category across variant periods (Table 3; Figures 3, 4) as well as among lineages for sequenced samples (Appendix Tables 6–8, Figures 2 and 3). Ct values were significantly lower during the Delta and Omicron periods than in the pre-Delta period (adjusted p value = 0.0004 for Delta period and adjusted p = 0.0003 for Omicron period) but did not differ significantly between Delta and Omicron (Table 3; Figure 3). We did not observe this pattern in the subset of sequenced samples (Appendix Tables 6–8, Figure 2). We observed a significant interaction between variant period and patient age (adjusted p = 0.001). Among patients infected during the Delta period, children <1 year of age had significantly lower Ct values than those >1 year of age (adjusted p = 0.0021 for children 1–4 years of age, adjusted p<0.0001 for 5–11 years, and adjusted p<0.0001 for >12 years) (Table 3; Figure 3). In addition, patients <1 year of age infected during Delta had significantly lower Ct values than children of the same age group infected during the Omicron period (adjusted p = 0.0156). Patients infected during the Delta period who were 1–4 years of age had lower Ct values than children in the same age group who were infected during the pre-Delta period (adjusted p = 0.031). We did not observe a significant interaction between lineage and patient age among sequenced samples (Appendix Tables 6–8, Figure 2). Furthermore, we did not observe a significant interaction between either variant period or sequenced lineage and vaccination status (Table 3; Figure 3; Appendix Tables 6–8, Figure 2).

When we evaluated potential indicators of disease severity across variant periods, we did not observe a statistically significant interaction between the variant period and patient type (inpatient because of COVID-19 vs. outpatient), the presence or absence of symptoms, days from symptom onset, or whether inpatients were admitted to the PICU (Table 3; Figure 4). Among patients who did not receive supplemental oxygen, those infected during the pre-Delta period had higher Ct values than those infected during the Delta (adjusted p = 0.0017) or Omicron (adjusted p = 0.0238) periods, in addition to having higher Ct values than children infected during pre-Delta who received supplemental oxygen (adjusted p = 0.0121). However, when compared by the highest level of supplemental oxygen support received (none, invasive, or noninvasive), results were not statistically significant. Among sequenced samples, we did not observe any statistically significant interactions between lineage and disease severity indicators, likely because of Ct value bias and smaller sample sizes (Appendix Tables 6–8, Figure 3).

Discussion

We found that children infected with SARS-CoV-2 during the Delta and Omicron periods had lower Ct values (suggesting higher viral load) than children infected during the pre-Delta period. Patients who were <1 year of age had lower Ct values than children in other age groups, a pattern potentially driven by Delta. In addition, inpatients with symptomatic infections had lower Ct values than asymptomatic inpatients, but we did not observe significant differences in Ct values among disease severity markers among variant periods. Therefore, Ct value might not be a strong correlate of disease severity outcomes across variant periods in children. During the Delta and Omicron periods, a greater proportion of inpatients were symptomatic, hospitalized in the PICU, or received oxygen support because of COVID-19 than during the pre-Delta period, but those proportions did not differ significantly between Delta and Omicron. This finding suggests that, for hospitalized children, Delta and Omicron infections might be more severe than earlier variants, regardless of Ct value.

Although many studies have evaluated relationships between SARS-CoV-2 variants and Ct values or viral load (6,7), transmission (23), and symptom severity (3,5,7) in adults, fewer studies have focused on children. Therefore, the effects of SARS-CoV-2 variant on Ct values or viral load and symptom severity in this population remain poorly understood.

Previous studies in adults observed that those infected by Delta exhibit higher viral loads (6) than those infected by previous variants and those infected by Omicron exhibit similar viral loads to those infected by Delta (7). Consistent with our results, previous studies in children observed lower Ct values in hospitalized children infected with Delta and Omicron than for previous variants (15). On the other hand, studies regarding the relationship between viral load and patient age yielded mixed results (24). Consistent with our results, others have found that infants (25) and young children (26,27) exhibit higher viral loads than older children. One study observed that older children exhibit similar viral loads to adults (25), whereas others found no significant effect of age on viral load (28–30). Therefore, more work with larger sample sizes or meta-analyses will be required to describe the relationship between age and Ct values or viral load in pediatric cases.

Although not directly addressed in this study, differences in viral loads among age groups or among SARS-CoV-2 variants have the potential to influence SARS-CoV-2 transmission dynamics. Early household transmission studies found that children transmit SARS-CoV-2 at lower rates than adults (23,31). However, more recent studies showed increased transmission by children in schools (32) and athletic facilities (33) during Delta compared with earlier variants. In addition, Omicron is more contagious than previous variants, including Delta (34). Finally, Costa et al. (29) found faster RNA clearance in children than in adults, suggesting a shorter viral shedding period, which might contribute to differing transmission dynamics between adults and children. This finding suggests that children might have a more significant effect on SARS-CoV-2 transmission dynamics than previously thought.

Consistent with our study, Chung et al. (28) found lower Ct values in symptomatic children than asymptomatic children. However, we did not observe a significant relationship between Ct values and markers of disease severity across variant periods among hospitalized children. Instead, higher proportions of hospitalized children exhibited disease severity markers during the Delta and Omicron periods than during pre-Delta, regardless of Ct values. Similarly, Quintero et al. (15) found that children infected by Delta were more likely to be symptomatic than those infected by earlier variants, and Mitchell et al. (35) found the highest proportion of PICU admissions after the emergence of Omicron. A review by Khemiri et al. (36) observed similar clinical manifestations in children and adults infected by Delta and Omicron. Furthermore, children, particularly those <5 years of age, were at higher risk for hospitalization (8,37–39) and had a higher incidence of croup (40) during Omicron than during earlier variants. In sum, those studies suggest that Delta and Omicron infections might be more severe in children than previous variants.

The first limitation of our study is that analyses of the relationship between WGS-based SARS-CoV-2 lineage calls and Ct values might be biased toward lower Ct values because the CDPHE Laboratory uses a cutoff of <28 Ct for WGS. As such, we performed our analyses using variant periods to include all samples for which we had Ct data, regardless of whether they were successfully sequenced, with the understanding that a small number of samples might be assigned to the incorrect variant period. We found that our analysis lost some power when only looking at samples that were successfully sequenced, but the results largely agree with the full dataset that used variant periods. Second, although the relationship between Ct values and quantitative viral load are not exact, they are tightly inversely correlated (13). Because we used a single quantitative RT-PCR protocol and platform for all samples, we expect the Ct patterns, as they relate to viral load, to be consistent. Third, whether samples with high Ct values represent the beginning or end of infection is unclear. However, it is assumed that most patients are tested early in the course of infection. Fourth, because we concluded data collection in January 2022, our analyses of Omicron only include the sublineages BA.1.1.529, BA.1, and BA.1.1. Because of the rapid evolution of SARS-CoV-2, numerous Omicron sublineages have emerged since and are not included here. Those other sublineages could exhibit different characteristics, such as different Ct values and altered transmissibility or virulence (34). Future work could address those differences among Omicron sublineages. Finally, our analyses only included Colorado residents and therefore might not be generalizable to other jurisdictions or geographic locations.

Our study provides insight into the relationship between Ct values, the presence of symptoms, and disease severity markers across variant periods. A strength of our study is the inclusion of in-depth data for both inpatient and outpatient populations. Similar studies are limited to performing in-depth analyses on Ct values to inpatient populations (37). In addition, our study is unique in its integration of clinical and laboratory data. Furthermore, rather than relying on RT-PCR of a small number of mutations to delineate variant periods, we were able to use entire genomes, enabling more accurate lineage determination.

In conclusion, our results suggest that children infected with either Delta or Omicron had lower Ct values (suggesting higher viral load) and potentially greater disease severity than children infected with earlier variants. Those patterns were particularly pronounced in the youngest children (<1 year of age). Our data highlight the importance of monitoring children, particularly in younger age groups, as new variants emerge and including pediatric cases in surveillance efforts for current and future SARS-CoV-2 variants to understand the potential differences in Ct values and viral loads, transmission, and disease severity across the age spectrum.

Dr. Bankers is the Bioinformatics and Genomic Analysis supervisor at the Colorado Department of Public Health and Environment Laboratory in Denver. Her research interests focus on the intersection of evolutionary genomics, pathogen genomic surveillance, and public health.

Top

Acknowledgments

We thank the Microbiology Laboratory at CHCO for submitting samples to the CDPHE Lab, the accessioning team at the CDPHE lab for logging in all samples, and the sequencing team for sample processing. We thank Brian Erly for helpful comments and edits on an earlier draft.

The data included in this work were generated for public health surveillance purposes and was funded by 6 NU50CK000552-02-01 Epidemiology and Laboratory Capacity for Infectious Diseases (ELC) and NU50CK000483 Emerging Infections Program (EIP) Cooperative Agreements from the Centers for Disease Control and Prevention.

Top

References

  1. Harvey  WT, Carabelli  AM, Jackson  B, Gupta  RK, Thomson  EC, Harrison  EM, et al.; COVID-19 Genomics UK (COG-UK) Consortium. SARS-CoV-2 variants, spike mutations and immune escape. Nat Rev Microbiol. 2021;19:409–24. DOIPubMedGoogle Scholar
  2. Telenti  A, Hodcroft  EB, Robertson  DL. The evolution and biology of SARS-CoV-2 variants. Cold Spring Harb Perspect Med. 2022;12:a041390. DOIPubMedGoogle Scholar
  3. Twohig  KA, Nyberg  T, Zaidi  A, Thelwall  S, Sinnathamby  MA, Aliabadi  S, et al.; COVID-19 Genomics UK (COG-UK) consortium. Hospital admission and emergency care attendance risk for SARS-CoV-2 delta (B.1.617.2) compared with alpha (B.1.1.7) variants of concern: a cohort study. Lancet Infect Dis. 2022;22:35–42. DOIPubMedGoogle Scholar
  4. Lambrou  AS, Shirk  P, Steele  MK, Paul  P, Paden  CR, Cadwell  B, et al.; Strain Surveillance and Emerging Variants Bioinformatic Working Group; Strain Surveillance and Emerging Variants NS3 Working Group. Genomic Surveillance for SARS-CoV-2 Variants: Predominance of the Delta (B.1.617.2) and Omicron (B.1.1.529) Variants - United States, June 2021-January 2022. MMWR Morb Mortal Wkly Rep. 2022;71:206–11. DOIPubMedGoogle Scholar
  5. Taylor  CA, Patel  K, Pham  H, Whitaker  M, Anglin  O, Kambhampati  AK, et al.; COVID-NET Surveillance Team. Severity of disease among adults hospitalized with laboratory-confirmed COVID-19 before and during the period of SARS-CoV-2 B.1.617.2 (Delta) predominance— COVID-NET, 14 states, January–August 2021. MMWR Morb Mortal Wkly Rep. 2021;70:1513–9. DOIPubMedGoogle Scholar
  6. Wang  Y, Chen  R, Hu  F, Lan  Y, Yang  Z, Zhan  C, et al. Transmission, viral kinetics and clinical characteristics of the emergent SARS-CoV-2 Delta VOC in Guangzhou, China. EClinicalMedicine. 2021;40:101129. DOIPubMedGoogle Scholar
  7. Ulloa  AC, Buchan  SA, Daneman  N, Brown  KA. Estimates of SARS-CoV-2 Omicron variant severity in Ontario, Canada. JAMA. 2022;327:1286–8. DOIPubMedGoogle Scholar
  8. Taylor  CA, Whitaker  M, Anglin  O, Milucky  J, Patel  K, Pham  H, et al.; COVID-NET Surveillance Team. COVID-19–associated hospitalizations among adults during SARS-CoV-2 Delta and Omicron variant predominance, by race/ethnicity and vaccination status—COVID-NET, 14 states, July 2021–January 2022. MMWR Morb Mortal Wkly Rep. 2022;71:466–73. DOIPubMedGoogle Scholar
  9. Ferdinands  JM, Rao  S, Dixon  BE, Mitchell  PK, DeSilva  MB, Irving  SA, et al. Waning 2-dose and 3-dose effectiveness of mRNA vaccines against COVID-19–associated emergency department and urgent care encounters and hospitalizations among adults during periods of Delta and Omicron variant predominance—VISION Network, 10 states, August 2021–January 2022. MMWR Morb Mortal Wkly Rep. 2022;71:255–63. DOIPubMedGoogle Scholar
  10. Delahoy  MJ, Ujamaa  D, Whitaker  M, O’Halloran  A, Anglin  O, Burns  E, et al.; COVID-NET Surveillance Team. Hospitalizations associated with COVID-19 among children and adolescents—COVID-NET, 14 states, March 1, 2020–August 14, 2021. MMWR Morb Mortal Wkly Rep. 2021;70:1255–60. DOIPubMedGoogle Scholar
  11. Wang  L, Berger  NA, Kaelber  DC, Davis  PB, Volkow  ND, Xu  R. Incidence rates and clinical outcomes of SARS-CoV-2 infection with the Omicron and Delta variants in children younger than 5 years in the US. JAMA Pediatr. 2022;176:811–3. DOIPubMedGoogle Scholar
  12. Belay  ED, Godfred-Cato  S. SARS-CoV-2 spread and hospitalisations in paediatric patients during the omicron surge. Lancet Child Adolesc Health. 2022;6:280–1. DOIPubMedGoogle Scholar
  13. Walker  AS, Pritchard  E, House  T, Robotham  JV, Birrell  PJ, Bell  I, et al.; COVID-19 Infection Survey team. Ct threshold values, a proxy for viral load in community SARS-CoV-2 cases, demonstrate wide variation across populations and over time. eLife. 2021;10:e64683. DOIPubMedGoogle Scholar
  14. Rostad  CA, Kamidani  S, Anderson  EJ. Implications of SARS-CoV-2 viral load in children: getting back to school and normal. JAMA Pediatr. 2021;175:e212022. DOIPubMedGoogle Scholar
  15. Quintero  AM, Eisner  M, Sayegh  R, Wright  T, Ramilo  O, Leber  AL, et al. Differences in SARS-CoV-2 clinical manifestations and disease severity in children and adolescents by infecting variant. Emerg Infect Dis. 2022;28:2270–80. DOIPubMedGoogle Scholar
  16. DNA Pipelines R&D. Farr B, Rajan D, Betteridge E, Shirley L, Quail M, et al. COVID-19 ARTIC v3 Illumina library construction and sequencing protocol. 2020 Nov 4 [cited 2022 Feb 15]. https://www.protocols.io/view/covid-19-artic-v3-illumina-library-construction-an-j8nlke665l5r/v5
  17. Quick  J. nCoV-2019 sequencing protocol v3 (LoCost). 2020 Aug 25 [cited 2022 Feb 15]. https://www.protocols.io/view/ncov-2019-sequencing-protocol-v3-locost-bp2l6n26rgqe/v3
  18. Illumina CovidSeq Test. Instructions for use [cited 2022 Feb 15]. https://support.illumina.com/content/dam/illumina-support/documents/documentation/chemistry_documentation/Illumina-COVIDSeq-Test/illumina-covidseq-test-instructions-for-use-1000000128490-03.PDF
  19. The Broad Institute. Terra [cited 2022 Feb 22]. https://app.terra.bio
  20. Colorado Department of Public Health and Environment. CDPHE-bioinformatics/CDPHE-SARS-CoV-2 [cited 2022 Feb 22]. https://github.com/CDPHE-bioinformatics/CDPHE-SARS-CoV-2
  21. Colorado Department of Public Health and Environment. CDPHEsarscov2 [cited 2022 Feb 22]. https://dockstore.org/organizations/CDPHE/collections/CDPHEsarscov2
  22. Tableau (version. 9.1). J Med Libr Assoc. 2016;104:182–3. DOIGoogle Scholar
  23. Silverberg  SL, Zhang  BY, Li  SNJ, Burgert  C, Shulha  HP, Kitchin  V, et al. Child transmission of SARS-CoV-2: a systematic review and meta-analysis. BMC Pediatr. 2022;22:172. DOIPubMedGoogle Scholar
  24. Bayhan  Gİ, Özkubat Korkmaz  I, Şahiner  ES, Tekeli  N, Akdağ  AG, Uyan Erten  AZ, et al. Are SARS-CoV-2 viral loads in children lower than in adults? J Infect. 2023;86:e13–4. DOIPubMedGoogle Scholar
  25. Ochoa  V, Díaz  FE, Ramirez  E, Fentini  MC, Carobene  M, Geffner  J, et al.; INBIRS COVID-19 Study Group. Infants younger than 6 months infected with SARS-CoV-2 show the highest respiratory viral loads. J Infect Dis. 2022;225:392–5. DOIPubMedGoogle Scholar
  26. Cendejas-Bueno  E, Romero-Gómez  MP, Escosa-García  L, Jiménez-Rodríguez  S, Mingorance  J, García-Rodríguez  J; SARS-CoV-2 Working Group. Lower nasopharyngeal viral loads in pediatric population. The missing piece to understand SARS-CoV-2 infection in children? J Infect. 2021;83:e18–9. DOIPubMedGoogle Scholar
  27. Heald-Sargent  T, Muller  WJ, Zheng  X, Rippe  J, Patel  AB, Kociolek  LK. Age-related differences in nasopharyngeal severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) levels in patients with mild to moderate coronavirus disease 2019 (COVID-19). JAMA Pediatr. 2020;174:902–3. DOIPubMedGoogle Scholar
  28. Chung  E, Chow  EJ, Wilcox  NC, Burstein  R, Brandstetter  E, Han  PD, et al. Comparison of symptoms and RNA levels in children and adults with SARS-CoV-2 infection in the community setting. JAMA Pediatr. 2021;175:e212025–212025. DOIPubMedGoogle Scholar
  29. Costa  R, Bueno  F, Albert  E, Torres  I, Carbonell-Sahuquillo  S, Barrés-Fernández  A, et al. Upper respiratory tract SARS-CoV-2 RNA loads in symptomatic and asymptomatic children and adults. Clin Microbiol Infect. 2021;27:1858.e1–7. DOIPubMedGoogle Scholar
  30. Yonker  LM, Boucau  J, Regan  J, Choudhary  MC, Burns  MD, Young  N, et al. Virologic features of severe acute respiratory syndrome coronavirus 2 infection in children. J Infect Dis. 2021;224:1821–9. DOIPubMedGoogle Scholar
  31. Maltezou  HC, Vorou  R, Papadima  K, Kossyvakis  A, Spanakis  N, Gioula  G, et al. Transmission dynamics of SARS-CoV-2 within families with children in Greece: A study of 23 clusters. J Med Virol. 2021;93:1414–20. DOIPubMedGoogle Scholar
  32. Torjesen  I. Covid-19: Delta variant is now UK’s most dominant strain and spreading through schools. BMJ. 2021;373:n1445. DOIPubMedGoogle Scholar
  33. Sheikh  A, McMenamin  J, Taylor  B, Robertson  C; Public Health Scotland and the EAVE II Collaborators. SARS-CoV-2 Delta VOC in Scotland: demographics, risk of hospital admission, and vaccine effectiveness. Lancet. 2021;397:2461–2. DOIPubMedGoogle Scholar
  34. Setiabudi  D, Sribudiani  Y, Hermawan  K, Andriyoko  B, Nataprawira  HM. The Omicron variant of concern: The genomics, diagnostics, and clinical characteristics in children. Front Pediatr. 2022;10:898463. DOIPubMedGoogle Scholar
  35. Mitchell  R, Cayen  J, Thampi  N, Frenette  C, Bartoszko  J, Choi  KB, et al. Trends in severe outcomes among adult and pediatric patients hospitalized with COVID-19 in the Canadian Nosocomial Infection Surveillance Program, March 2020 to May 2022. JAMA Netw Open. 2023;6:e239050. DOIPubMedGoogle Scholar
  36. Khemiri  H, Ayouni  K, Triki  H, Haddad-Boubaker  S. SARS-CoV-2 infection in pediatric population before and during the Delta (B.1.617.2) and Omicron (B.1.1.529) variants era. Virol J. 2022;19:144. DOIPubMedGoogle Scholar
  37. Acker  KP, Levine  DA, Varghese  M, Nash  KA, RoyChoudhury  A, Abramson  EL, et al. Indications for hospitalization in children with SARS-CoV-2 infection during the Omicron wave in New York City. Children (Basel). 2022;9:1043. DOIPubMedGoogle Scholar
  38. Iuliano  AD, Brunkard  JM, Boehmer  TK, Peterson  E, Adjei  S, Binder  AM, et al. Trends in disease severity and health care utilization during the early Omicron variant period compared with previous SARS-CoV-2 high transmission periods—United States, December 2020–January 2022. MMWR Morb Mortal Wkly Rep. 2022;71:146–52. DOIPubMedGoogle Scholar
  39. Marks  KJ, Whitaker  M, Agathis  NT, Anglin  O, Milucky  J, Patel  K, et al.; COVID-NET Surveillance Team. Hospitalization of infants and children aged 0–4 years with laboratory-confirmed COVID-19—COVID-NET, 14 states, March 2020–February 2022. MMWR Morb Mortal Wkly Rep. 2022;71:429–36. DOIPubMedGoogle Scholar
  40. Iijima  H, Kubota  M, Ogimi  C. Clinical characteristics of pediatric patients with COVID-19 between Omicron era vs. pre-Omicron era. J Infect Chemother. 2022;28:1501–5. DOIPubMedGoogle Scholar

Top

Figures
Tables

Top

Cite This Article

DOI: 10.3201/eid3006.231427

Original Publication Date: May 14, 2024

Table of Contents – Volume 30, Number 6—June 2024

EID Search Options
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.

Top

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Laura Bankers, Colorado Department of Public Health and Environment, 8100 E Lowry Blvd, Denver, CO 80230, USA

Send To

10000 character(s) remaining.

Top

Page created: April 08, 2024
Page updated: May 22, 2024
Page reviewed: May 22, 2024
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external