ABSTRACT
Background Epstein-Barr virus-associated Infectious Mononucleosis (EBV-IM) is a common disease following primary EBV infection in children and adolescents. While EBV-IM is mostly self-limiting, symptoms like fatigue may persist over several months or even result in myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). This large clinical observational study aimed at identifying risk factors for protracted courses of EBV-IM in young people.
Methods A cohort of N=200 children, adolescents, and young adults with acute primary EBV infection was recruited from hospitals and private practices. Data on the patients’ medical history as well as clinical and laboratory parameters were collected at a baseline visit (V1) within four weeks after symptom onset (Tonset) and at two follow-up visits (V2 and V3) one and six months after Tonset. Risk factors for protracted symptoms at V3 were modeled using multivariable logistic regressions.
Results Protracted symptoms were observed in 55/183 (30.1%) and protracted fatigue in 34/181 (18.8%) patients at V3. A medical history indicating an increased susceptibility to infectious diseases as well as distinct severe IM symptoms, e.g. severe gastrointestinal symptoms, were significantly associated with protracted disease [OR: 2.31; P=0.011 and OR: 3.42; P=0.027] and with chronic fatigue [OR: 2.98; P=0.006 and OR: 3.54; P=0.034], respectively. Occurrence of twelve or more clinical and laboratory parameters until and including V1 discriminated between fatigue and no fatigue at V3 [OR 2.43, P=0.033].
Conclusion A clinical history of immune dysregulation as well as distinct severe IM symptoms might predict protracted post-viral disease and thus help in the identification of young patients at risk.
Highlights
Severe gastrointestinal symptoms are associated with protracted course of Epstein-Barr virus-associated Infectious Mononucleosis (EBV-IM).
Signs of immune dysregulation prior to EBV-IM can indicate an increased risk of protracted symptoms.
Greater number of initial symptoms helps to identify patients developing postviral fatigue.
INTRODUCTION
Infectious mononucleosis (IM) develops in up to 75% of adolescents and young adults after primary infection with Epstein-Barr virus (EBV) [1]. Symptoms typically include fever, lymphadenopathy, tonsillopharyngitis, splenomegaly, and fatigue, and usually are self-limiting within two weeks to three months [1]. However, acute complications such as kissing tonsils, splenic rupture, or life-threatening hemophagocytic lymphohistiocytosis (HLH) [1, 2] can occur, and lymphoproliferative disease can manifest in immunocompromised patients [3]. Moreover, post-viral symptoms may persist for weeks or years and may fulfill the clinical criteria of myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) [4, 5], significantly compromising daily activities [6].
Risk factors for long-term IM sequelae are still under investigation, and various scores have been published to assess IM severity, complexity, and protraction [7-13]. In a large cohort of college students, physical symptoms and immune dysregulation have been identified as risk factors for ME/CFS following IM [14], and autonomic symptoms along with distinct immune markers were suggested as predictors of severe ME/CFS [15].
Here, we report on a novel IM scoring system and on candidate risk factors for protracted symptoms identified in the Munich Infectious Mononucleosis (IMMUC) study. This observational, prospective, longitudinal study investigates a comprehensive set of routine and experimental data in 200 children, adolescents, and young adults with recent IM to identify biomarkers and causative factors of complicated disease trajectories, including chronic fatigue.
MATERIAL & METHODS
Ethics, Study Setting, and Inclusion Criteria
The IMMUC study was approved by the TUM University Hospital Ethics Committee (No. 112/14) and adheres to the 1964 Declaration of Helsinki, its later amendments, and the ICH-GCP regulations. All participants and/or legal guardians gave written informed consent. The study was registered at ClinicalTrials.gov (NCT06002802). 200 participants with recent IM onset were prospectively recruited from hospitals or private practices in or around Munich between February 2016 and June 2019 and observed up to 12 months after IM onset. General, clinical, and virological inclusion criteria are summarized in Tab. 1.
Study Visits, Study Periods and Study Data
Patients were seen at visit 1 (V1) within 28 days after onset (DAO) of IM (DAO 0=Tonset). V2 was scheduled at four to six weeks after Tonset and V3 at four to 12 months after Tonset. A defined set of virological (IMVir), clinical chemistry (IMLab), and clinical (IMClin) features of IM was investigated at each visit. Additionally, medical history data and external laboratory data were collected for the time periods before Tonset and between Tonset, V1, and V2, to generate combined datasets for the history period (HP) before Tonset, the very early period (VEP) from Tonset up to and including V1, and the early period (EP) from Tonset up to and including V2 (Fig. 1). All visits included a structured interview, a physical examination, and peripheral blood collection for routine (reported here) and experimental (reported elsewhere) laboratory analyses. Splenic sonography was performed at V1 and V2 and in case of persistent splenomegaly at V3. For a detailed account of the collected study data see Suppl. Methods, Suppl. Material M1, and Suppl. Material M2.
IM Scoring
A score was developed to grade severity (Smax), complexity (Cmax), and protraction of individual IM features and overall IM at each visit and for each period (Tab. 2). Up to five grades of severity (S0-5) were evaluated for each IM feature, and the highest severity grade of all distinct features was used to grade the overall IM (Suppl. Mat. M2). By calculating the number of IM features, complexity grades were determined for the overall IM (C0-C23), IMClin (C0-15), and IMLab (C0-8) IM. Protraction grades were defined as P0 (feature/IM not present anymore at V2 and V3), P1 (feature/IM present at V2 but not at V3), or P2 (feature/IM present at V3). Since the clinical significance of mild lymphadenopathy and splenomegaly at V3 was unclear, only S2-4 of the latter at V3 were counted as P2. The grading of fatigue at V3 differed from that at V1 and V2 due to an increasing clinical impact of fatigue over time.
Outcome Measures, Predictor Variables in Multivariable Logistic Regression Models
Long-term IM protraction (P2-IM vs. P0/1-IM) or chronic fatigue (P2-fatigue vs. P0/1-fatigue) were the binary outcome measures to be modelled. Predictor variables were either measured repetitively at visits V1 and V2 (quantitative and categorized IMVir, IMLab and IMClin features) and thus, included into the visit-dependent analyses (see below) or only once during the periods HP, VEP, and/or EP (quantitative and categorized baseline patient characteristics as well as Cmax-IM, -IMLab, -IMClin, Smax-IM, -IMLab, and -IMClin features) and thus included into the period-dependent analyses (see below).
Sample Size Calculation
We expected approximately 10% of patients with P2-IM, and therefore, set the sample size to N=200 patients (Suppl. Methods).
Visit-Dependent Analyses
Each multivariable logistic regression model included one of the predictor variables measured at V1 and V2, and sex and age as independent variables. Since individual patients were measured repetitively, an interaction effect between the predictor variable and DAO was included in each model. The effect of a predictor variable can thus be described in dependence on DAO. These models also determined a time interval at which the OR of the predictor variable was significantly different from 1. The robust Huber-White estimator accounted for minor misspecifications of the model assumption due to repeated measurements.
Period-Dependent Analyses
Each multivariable logistic regression model included one of the predictor variables measured once for either of the three periods: HP, VEP, or EP and sex and age as independent variables. A grid search was used to find the optimal cut-off value for Cmax-IM and Cmax-IMclin in terms of the area under the curve (AUC).
For the visit- and period-dependent analyses, models with less than five observations per independent variable per category variable of the outcome variable were discarded to prevent overfitting. Only significant results are presented here due to the abundance of variables investigated.
Other Statistical Analyses
Hypothesis tests were performed at exploratory 5% significance levels. P-values are exploratory and thus not corrected for multiple testing. Analyses were performed using R, version 4.1.3 (The R Foundation for Statistical Computing, Vienna, Austria) [16].
RESULTS
Of 200 study participants assessed at V1, 194 and 183 attended V2 and V3, respectively (Fig. 1A). Study visits (V1, V2, V3) and study periods (HP, VEP, EP) are depicted in Fig. 1B.
Baseline Patient Characteristics
Detailed data are shown in Suppl. Table S1A to S1F. The cohort included 92 children (1-11 years), 84 adolescents (12-17 years), and 24 adults (18-33 years). 54.5% patients were female. The mean age was 11.99±6.08 years (range: 1–33) (male: 11.76±6.89, female: 12.18±5.33; t-test: P=0.632). A history of increased susceptibility to infections and/or immune dysregulation (ELVIS and/or GARFIELD criteria [17]) prior to IM was reported by 47.5% of patients (Suppl. Table 1E). When asked for any fever or antibiotic therapy during the six months prior to IM, 37.0% and 13.6% of patients provided positive answers, respectively. 35.7% and 6.0% of patients remembered ≥1 or ≥2 family members with a history of IM.
Severity, Complexity, and Protraction of IM
Comprehensive data on IM severity, complexity, and protraction are provided in Fig. 3A-C, Suppl. Fig. S1, Suppl. Fig. S2 and Suppl. Tab. S2. IM severity and complexity decreased from V1 to V3 as expected. At least one IM feature persisted in 52.3% at V2 (P1) and in 27.9 % at V3 (P2) (Fig. 3C). 34 (18.8%) patients suffered from fatigue at V3, with mild fatigue (S1) reported by 82.4% and moderate fatigue (S2) (“reduced education/work/club sports”) by 17.6% of these patients (Fig. 3A,C).
Predictors for Long-Term Protraction (P2-IM)
Significant predictors of any IM feature at V3 (P2-IM) of the visit-dependent analyses are shown in Fig. 4A. Early positivity for VCAp18 IgG was the only risk factor, with decreasing effect over time (negative slope in Fig. 4A). Protective effects were also observed for early atypical lymphocytosis (until 6 DAO), tonsillopharyngitis (until 9 DAO) and leukocytosis (until 4 DAO). Complexity of IM and IMClin were significant risk factors for P2-IM from DAO=6 and DAO=5 onwards, respectively. Additional IMClin features predicted P2-IM, again with a delay of few days after Tonset and depending on severity (positive slope in Fig. 4A). Considering the AUC, the best predictor variables were complexity of IMClin (AUC 0.68), complexity of IM (AUC 0.66) and presence of other airway symptoms [S1-4 vs. S0] (AUC 0.62).
Significant predictors of any IM feature at V3 (P2-IM) of the period-dependent analyses are shown in Fig. 5A. IMClin complexity in the VEP and EP were the best predictors in terms of the AUC. The ROC analysis indicated that a cut-off of ≥10 and ≥9 IMClin features during the VEP (OR 2.61, 95%-CI: [1.27, 5.37]) and EP (OR 2.63, 95%-CI: [1.37, 5.03]) provided the greatest AUCs for predicting P2-IM, respectively. Analogously, a cut-off of ≥12 and ≥13 IM features provided the greatest AUCs for the VEP (OR 2.47, 95%-CI: [1.26, 4.87]) and EP (OR 2.99, 95%-CI: [1.49, 5.98]), respectively. Interestingly, ≥1 positive ELVIS criterion was also identified as a significant risk factor for P2-IM.
Predictors for Chronic Fatigue (P2-fatigue)
Significant predictors of fatigue at V3 (P2-fatigue) of the visit-dependent analyses are shown in Fig. 4B. Fatigue was the most frequent IM feature at V3 (Suppl. Fig. S3E). It was, therefore, not unexpected to find overlaps between predictors for P2-fatigue and P2-IM. The best predictor for P2-fatigue regarding AUC was other airway symptoms which were significant risk factors after DAO 17, with subsequent OR increase over time. Moreover, moderate to severe GI symptoms and VCAp18 IgG detection in immunoblots were significant predictors with an increasing and decreasing OR over time, respectively (Fig. 4B). We did not find a significant visit-dependent effect of (very) early fatigue for P2-fatigue.
Significant predictors of fatigue at V3 (P2-fatigue) of the period-dependent analyses are shown in Fig. 5B. Early exanthema, early musculoskeletal symptoms, and ≥8 minor infections per year (ELVIS criterion) prior to IM significantly predicted P2-fatigue. Early IMClin complexity was the best predictor (AUC 0.74) of the period-dependent analyses. ROC analysis indicated that a cut-off of ≥12 IM features in the VEP (OR 2.43, 95%-CI: [1.07, 5.49]), and of ≥11 features in the EP (OR 3.39, 95%-CI: [1.41, 8.17]), as well as ≥9 IMClin features in both periods (VEP: OR 2.47, 95%-CI: [1.14, 5.33]; EP: OR 4.41, 95%-CI: [1.90, 10.21]) resulted in the greatest AUC for predicting fatigue at V3 (P2-fatigue), respectively. Early fatigue was not a significant predictor of P2-fatigue.
DISCUSSION
Long-term sequelae of IM, such as chronic fatigue, manifest in a significant number of cases [19, 20] and hit young people in a very vulnerable stage of their lives, indicating a high medical need for research on risk factors, prevention, and therapeutic strategies [18, 19]. To identify biomarkers and predictors of acute and long-term IM complications we implemented the longitudinal IMMUC study and report here on the first results from a cohort of 200 pediatric and young adult patients.
In line with published data, 72.1% of the patients had fully recovered five to eight months after IM onset [14]. Remarkably, 27.9% of the study participants showed protracted disease for more than four months after IM onset, including 18.8% cases with chronic fatigue. Several clinical and laboratory parameters were identified that may predict protracted IM and thus identify patients at risk.
Of note, IM complexity was a significant predictor of long-term protraction. Within the VEP and EP, 50.5% and 94% study participants presented with maximum IM complexity. Ten or more very early IMClin features predicted protracted disease at V3 and nine or more predicted chronic fatigue. (Very) early IMLab features did not add any essential value to prediction accuracy.
Our findings are consistent with those reported by Cervia and colleagues, who identified the complexity of coronavirus disease 2019 (COVID-19) in combination with other clinical and laboratory characteristics as risk factors for post-acute COVID-19 syndromes [20]. It was hypothesized that a greater number of symptoms during the acute infection, including autonomic dysregulation, might reflect an initial immune (over)activation that results in prolonged impaired health [15, 21]. Our results are also in line with reports indicating an increased risk of ME/CFS in patients with more or more severe physical symptoms during the acute infection [13, 14].
Patients with distinct inborn errors of immunity (IEI) carry an increased susceptibility to EBV-associated HLH or lymphoproliferative disorders [3]. However, little is known about clinical IEI warning signs as possible predictors of protracted IM. In a previous study, 60% of ME/CFS patients reported a history of recurrent infections [22]. Our IMMUC data demonstrate that ≥1 major infection and ≥ 8 minor infections per year predicted any IM feature or fatigue at V3, respectively. These findings suggest that patients with an increased susceptibility to infection should be closely monitored and offered timely advice after IM diagnosis.
In addition, chronic fatigue was predicted by (very) early exanthema and facial edema, possibly indicative of IM-associated vascular inflammation [23-25], and was associated with early musculoskeletal symptoms and (very) early moderate to severe GI symptoms. Interestingly, GI symptoms prior to IM have been associated with the development of ME/CFS [15]. Taken together, these findings suggest that preexisting or early autonomic symptoms in the context of IM might predict post-viral disease.
Remarkably, the detection of EBV VCAp18 IgG considered a late marker of primary EBV infection [26], was identified as a significant predictor of any long-term protraction when detected within the first three weeks after IM onset. This might indicate that individuals with a long EBV incubation period prior to IM may face an increased risk of protracted disease, and early immunoblot analysis might help in prediction. Early detection of EBV VCAp18 IgG has been associated with a lower IM severity score by others but was not investigated as a predictor of protracted disease [27]. Moreover, EBV VCAp18 IgM-mediated complement deposition and EBV VCAp18 IgG-dependent cellular phagocytosis were positively associated with IM severity [28], and complement activation was discussed in ME/CFS-associated immune dysregulation [22].
Limitations
First, the results of the IMMUC scores might have been prone to interviewer bias. Second, most patients were recruited when initial IM symptoms had declined. Therefore, laboratory parameters at V1 may reflect a subacute stage of IM and remembered IM symptoms of the VEP might be influenced by recall bias. Third, the time period between V1 and follow-up visits varied to match best with school schedules and family plans.
Conclusion
In this IMMUC study, we implemented a novel scoring system to evaluate IM severity, complexity, and protraction and identified several candidate predictors of chronic postviral disease that are easily assessed during routine diagnostics. This score may help identifying patients at risk and facilitate prompt medical provisioning. In addition, this score forms the basis for the analysis of diverse immunological parameters investigated in this cohort that might further improve prediction for long-term sequelae in young people with IM.
Data Availability
All data produced in the present study is not openly available.
Author contributions
U.B. acquired funding, supervised and coordinated the study together with K.G.. U.B., M.B., J.G., A.H., S.S, N.K., T.B., J.M, D.H., P.L., S.E-S., H-J.D., S.D., F.H., C.F., T.S., M-M.S., A.B., and K.G. designed the study. M.B., J.G., L.S-H., R.P., F.H., C.F., M-M.S., A.B., K.G., and U.B. recruited and visited the patients and acquired the data. M.B., J.G., J.L.L., J.W., A.H., E.N., F.H., M-M.S., A.B., A.N., L.M., K.G., and U.B. analysed and interpreted the data. M.B., J.G., L.M., J.L.L., A.H., A.N., K.G., and U.B. wrote the first version of the manuscript. All authors read, edited, and commented on the draft. All authors approved the final version of the manuscript to be published and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Availability of data and materials
M.B., J.G., L.M., J.L.L., A.H., J.W., K.G., and U.B. have full access to all data and take responsibility for their integrity and accuracy.
Ethical considerations
The IMMUC study was approved by the Ethics Committee of the School of Medicine and Health of the Technical University of Munich (reference number: 112/14).
Funding
This work was funded by the German Center of Infection Research (DZIF) (grant numbers 07_905 and 07_909). Funding was acquired by U.B. and received by P.L., S.E-S., A.H., D.H., F.H., J.M., T.B., H-J.D., C.F., and T.S..
Potential conflicts of interest
A.B. received lecture fees from Siemens. U.B. received research grants for ME/CFS or EBV studies from the Federal Ministry of Education and Research, the Federal Ministry of Health, the Bavarian State Ministry of Health and Care, the Bavarian State Ministry of Science and the Arts, the German Center for Infection Research, the People for Children (Menschen für Kinder) foundation, the Weidenhammer-Zoebele Foundation, the Lost Voices Foundation, and the ME/CFS research foundation. All other authors have no conflict of interest to declare.
SUPPLEMENTARY METHODS
Study Data
Baseline patient characteristics included age, sex, body mass index (BMI), body height, serum titers (anti-diphtheria and anti-tetanus toxoid IgG, total IgA, IgM, IgG, IgE, vitamin D), and data from the general medical history of patients and families (Suppl. Tab. 1). The patient’s medical history addressed clinical signs of pathological susceptibility to infections (ELVIS) [17] and immune dysregulation (GARFIELD) [17], failure to thrive, developmental delay, vaccination status, infectious diseases (influenza, herpes labialis, varicella, measles), allergies, asthma, cancer, transplantation, transfusions, exposure to farm animals, as well as specific events within 6 months before Tonset (fever, use of antibiotics/analgesics/ antiphlogistics, traveling to certain regions). The family’s medical history addressed parent consanguinity, the number of siblings in the same household, as well as a history of IM, known inborn errors of immunity (IEI), cancer, neurodermatitis, autoimmune disease, allergies, and/or asthma (Suppl. Mat. M1). IMClin features included fever, fatigue, lymphadenopathy, tonsillopharyngitis, splenomegaly, facial edema, exanthema, other airway symptoms, other gastrointestinal (GI) symptoms, urogenital symptoms, bleeding and/or treatment for bleeding or anemia, neurological symptoms, musculoskeletal symptoms, hospitalization, and treatment for HLH (Suppl. Mat. M2); they were evaluated retrospectively and prospectively from Tonset until V3. IMLab features included lymphocytosis, neutropenia, anemia, thrombocytopenia, hyperferritinemia, hepatitis, nephropathy, and CRP elevation defined by age-adapted cut-off levels of blood count, transaminases (ASAT, ALAT), CRP, ferritin, creatinine, and (cystatin C if creatinine levels were elevated) in serum. IMVir features included EBV DNA in plasma (EBVPlasma), cell fraction (EBVcells), and oral washes (EBVoral) evaluated by PCR (TaqManTM 7500, Thermo Fisher Scientific, Waltham, USA), as well as anti-EBV VCA IgM, anti-EBV VCA IgG, and anti-EBNA-1 IgG determined in serum by CMIA [ARCHITECT i1000SRTM, Abbot Laboratories, Chicago, USA]), and serum IgG against EAp54, EAp138, VCAp18, VCAp23, BZLF1 or EBNA-1 detected by immunoblot (ProfiBlotTM, Tecan Group Ltd., Männedorf, Switzerland).
Sample Size Calculation
With 20 expected P2 cases and a required minimum number of 5-10 observations in the less frequent outcome class, we could consistently estimate multivariable logistic regression models with 2-4 parameters [31, 32]. An unpaired two-sided t-test with a 5% significance level has 80% power to detect an effect size of 0.7 when the sample sizes in two groups are N=18 (P2) and N=180 (P0-1) in a total sample size of N=198, as expected (nQuery Advisor 7.0).
SUPPLEMENTARY TABLES
SUPPLEMENTARY FIGURE LEGENDS
Acknowledgments
The authors would like to thank the IMMUC Study Group team members F. Fischer, M. Burggraf, P. Wallraven, R. Weggel, K. Bartl, L. Werny, K. Rautter, T. Hofberger, Y. Müller, H. Zietemann, E. Planatscher, J. Nückel, K. Michel, S. Strunz., J. Mücke, K. Mittelstraß, M. Bach, F. Martin, L. Kramer, C. Richter, A. Yardim, and L. Beckert for their great study assistance. Furthermore we would like to thank patients and their families for participating in the study, as well as all participating physicians for announcing the study and transferring interested patients. Without their tremendous support this work would have not been possible.
Footnotes
↵*** IMMUC Study Group team members are listed in the Acknowledgments.