Immune-metabolic trajectories delineate subgroups in paediatric long COVID

Abstract Most children and adolescents recover rapidly from SARS-CoV-2 infection, yet a subset develops paediatric long COVID (LC). How immune ontogeny shapes LC biology and heterogeneity remains unclear. We deeply phenotype a two-visit cohort with severe LC ( n = 74) and controls ( n = 27) spanning up to 3.2 years post index infection. Symptom burden remains high and neurofilament light chain (NfL) percentiles inversely associate with functional status (Bell score; $$r$$ r = −0.3536, P = 0.0060). Cardiopulmonary assessment and serology are unremarkable. Conventional autoantibodies are not enriched, whereas anti-DFS70 supports subgrouping. Immune features are temporally structured; SARS-CoV-2–associated mediators decline within 1 year, while innate-weighted, Th2-skewed cytokines persist. Metabolomics (43 metabolites) recapitulate the identified subgroups and align with EBV serostatus, disease phase (