Skip to main content

Active EBV infection in children: associations between DNA load, infection status, immune status, and disease severity

Abstract

Background

This study investigated active Epstein-Barr virus (EBV) infection in children and examined the associations among EBV deoxyribonucleic acid (DNA) load, infection types, disease severity, and immune characteristics.

Methods

A total of 35,956 pediatric patients who underwent EBV DNA load testing were included. Patients were categorized based on their EBV DNA levels and infection status.

Results

Spearman’s rank correlation analysis revealed a positive association between EBV DNA levels and the mortality rate, as well as the incidence rates of acute kidney injury (AKI), respiratory failure, cardiovascular complications, coagulation abnormalities, and liver injury. Mortality risk significantly increased when EBV DNA exceeded 1 × 105 copies/mL (adjusted odds ratio: 10.53, 95% confidence interval: 2.38–46.59, P < 0.05). As EBV DNA levels increase, the rise in mortality rate during activation- immunoglobulin G (IgG+) was more pronounced than that observed during primary infections. Gaussian mixture model clustering identified two immune clusters. Cluster 0 exhibited elevated pro-inflammatory indicators (IFN-γ, IL-6) and anti-inflammatory indicator (IL-10) levels, along with reduced immune cell counts. This cluster showed higher activation-IgG+ and mortality rates compared with Cluster 1.

Conclusions

An elevated EBV DNA load (> 1 × 105 copies/mL) in children is associated with increased mortality risk. High pro-inflammatory and anti-inflammatory states, coupled with low immune cell numbers, indicate critical condition. Simultaneous examinations of EBV DNA, antibodies, and immune status are recommended, especially for children with EBV DNA > 1 × 105 copies/mL, emphasizing the need for caution in those with activation-IgG+ and immune dysregulation.

Background

Epstein-Barr virus (EBV), a member of the gammaherpesvirus family, is a widely prevalent herpesvirus [1] infecting over 90% of adults worldwide. EBV infection commonly occurs during childhood [2,3,4]. Once infected, EBV establishes persistent infection in B lymphocytes and epithelial cells, typically in a latent form. A significant characteristic of EBV is its ability to reactivate multiple times after the initial infection, leading to secondary or recurrent infections. Both primary infection and reactivation of EBV can result in severe diseases, such as infectious mononucleosis (IM), lymphohistiocytosis (HLH), sepsis, and lymphoma [46]. EBV infection can also cause severe damage to distant organs, such as the liver, kidney, respiratory, and nervous system. This organ damage primarily arises from immune-mediated damage caused by lymphocyte infiltration and cytokine release post-EBV infection [4, 6,7,8], indicating disease progression and increasing the risk of life-threatening complications. Therefore, identifying an effective method to recognize the critical status of EBV infection is crucial for timely intervention to prevent and reduce the risk of severe conditions associated with infection.

EBV is typically diagnosed by detecting EBV-specific antibodies and measuring deoxyribonucleic acid (DNA) loads in body fluids. An elevated EBV DNA load is often associated with severe EBV-related diseases, such as HLH [9], lymphoma [10], and nasopharyngeal carcinoma [11, 12]. Patients with high EBV DNA loads often exhibit more severe disease symptoms and poorer prognosis. A study on lymphoma [13] has shown that patients who are EBV DNA positive have lower progression-free survival and overall survival rates. Reportedly, high EBV DNA levels negatively impact the immune system, contributing to poor outcomes. For instance, a study on the coronavirus disease 2019 (COVID-19) has shown that higher EBV DNA loads correlate with lower immune cells, particularly CD8+ T cells and natural killer (NK) cells, and increased disease severity [14]. However, research on the clinical value of EBV DNA load in pediatric EBV-related diseases remains limited. Although EBV DNA levels have been linked to prognosis in diseases such as HLH, lymphoma, and nasopharyngeal carcinoma, key questions remain unanswered. Specifically, the direct association between EBV DNA load and the risk of organ damage or poor prognosis remains unclear. Additionally, it is unclear whether a specific EBV DNA threshold can reliably guide clinical assessments regarding disease severity. In populations with activated EBV DNA, disease severity differs between patients with primary infection and those with viral reactivation. EBV infection has a synergistic pathogenic effect, leading to immune system damage. Conversely, preexisting immune disorders can trigger EBV reactivation and worsen the clinical condition. This raises the question: Is specific immune status associated with adverse outcomes in individuals with EBV activation?

Therefore, we conducted a study on hospitalized children who underwent EBV testing results to investigate the clinical significance of EBV DNA load in pediatric EBV-related diseases. However, the definition of EBV reactivation remains unclear. In this study, positive EBV immunoglobulin G (IgG) antibodies were used as indicates of EBV reactivation. We analyzed the relationship between EBV DNA levels and organ injury or adverse outcomes during hospitalization in children and explored the differences between primary infection and activation-IgG+. Additionally, we explored whether a specific immune profile was associated with adverse outcomes in individuals with EBV activation. A comprehensive analysis of these factors provides novel scientific evidence to enhance diagnostic accuracy and treatment strategies.

Methods

Study population

A retrospective analysis was conducted on clinical data from pediatric patients admitted to Hunan Children’s Hospital between January 2017 and August 2023 who undergone EBV DNA testing. The inclusion criteria included (1) age < 18 years, and (2) underwent EBV DNA testing with complete clinical data. The study was approved by the Medical Ethics Committee of Hunan Children’s Hospital (HCHLL-2022-50) and adhered to the ethical standards outlined in the Helsinki Declaration of 1964 and its subsequent amendments or similar ethical standards. The requirement for written informed consent was waived by the Medical Ethics Committee of the Hunan Children’s Hospital.

Measurement of EBV DNA, EBV antibodies, cytokines, lymphocyte subgroups, and liver function markers

Plasma EBV DNA levels were assessed by real-time fluorescence-based quantitative polymerase chain reaction (PCR) using an EBV DNA kit (Sansure Biotech Inc., Ltd of Changsha, Hunan Province, China). Plasma EB viral capsid antigen (VCA) immunoglobulin M (IgM), VCA IgG, early antigen (EA) IgM, and nuclear antigen (NA) IgG levels were measured via chemiluminescence immunoassay using an EBV antibody kit (Shenzhen YHLO Biotech Co, Ltd, China). Serum interleukin (IL)-2, IL-4, IL-6, IL-10, tumor necrosis factor (TNF)-α, and interferon (IFN)-γ levels were measured using multiplex bead-based flow cytometry with the cytokine detection kit (Qindao Raisecare Biological Technology Co., Ltd, China). Lymphocyte subsets, including NK cells (CD3/CD16+CD56+), total B lymphocytes (CD3CD19+), total T lymphocytes (CD3+CD19), helper/inducible T lymphocytes (CD3+CD4+), inhibitory/cytotoxic T lymphocytes (CTLs) (CD3+CD8+) and T lymphocytes + B lymphocytes + NK cells were measured by cytometry using a lymphocyte subsets kit (Shenzhen Mindray Bio-Medical Electronics Co., Ltd, China). Liver function markers—serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels—were measured using an automated biochemical analyzer with a liver function test kit (Beijing Jiuqiang Biotech, China). The testing procedures were performed by the Laboratory Department of Hunan Children’s Hospital following a standardized protocol.

Definition of EBV infection status

EBV infection states were classified based on established criteria [6, 15]: primary infection was defined as (1) the presence of VCA IgM, (2) low-affinity VCA IgG, or (3) a negative result in the EBV antibody test for both VCA IgM and VCA IgG with detection of EBV-DNA level > 400 copies/mL in plasma. Past infection was defined as the presence of high-affinity VCA IgG and/or IgG antibodies against EBV NA (EBNA). Activation-IgG+ was identified as a positive EBV antibody test for VCA IgG, indicating a past infection, along with an EBV-DNA level exceeding 400 copies/mL or the presence of EBV-VCA IgM in plasma. “Unknown status” was defined as patients who had not undergone EBV antibody testing. No-infection was defined as a negative result in the EBV antibody test for both VCA IgM and VCA IgG, with an EBV-DNA level ≤ 400 copies/mL in plasma.

Group definitions

All patients were categorized into six groups based on EBV DNA levels: 1) < 400 copies/mL; 2) > 400 copies/mL and ≤ 1 × 103 copies/mL; 3) > 1 × 103 copies/mL and ≤ 1 × 104 copies/mL; 4) > 1 × 104 copies/mL and ≤ 1 × 105 copies/mL; 5) > 1 × 105 copies/mL < and ≤ 1 × 106 copies/mL; and 6) > 1 × 106 copies/mL. Additionally, patients were further classified into three groups based on their infection status: primary infection group, activated infection (activation-IgG+ group), and the unknown status. Patients were also stratified into survival and death groups based on post-discharge outcomes.

Data collection

Information on underlying and accompanying diseases during hospitalization was extracted from the medical records. Underlying diseases included IM, sepsis, HLH, pneumonia, chronic active EBV (CAEBV), lymphoma, and immunodeficiency. Accompanying diagnoses included toxic encephalopathy, acute kidney injury (AKI), respiratory failure, acute respiratory distress syndrome (ARDS), myocardial injury, heart failure, shock, coagulation dysfunction, disseminated intravascular coagulation (DIC), and liver dysfunction. The diagnosis of CAEBV followed [16]: (1) persistent or recurrent IM-like symptoms for over 3 months; (2) detection of EBV antibodies (EBV-CA and EBV-EA) in tissues or peripheral blood samples, positive cells for EBV-encoded small RNA (EBER) in tissues, or EBV-DNA levels > 102.5 copies/mL in plasma and whole blood; and (3) absence of identifiable underlying immunodeficiency disease. The patients were required to meet al.l three criteria for the diagnosis of CAEBV. Lymphoma diagnosis was based on histopathological examination and immunohistochemical analysis according to the 2001 World Health Organization classification of lymphomas [17]. Other diseases were defined as patients who were not diagnosed with any of the following: IM, sepsis, HLH, pneumonia, CAEBV, or lymphoma. The diagnosis of AKI followed the 2012 Kidney Disease: Improving Global Outcomes (KDIGO) clinical practice guidelines, defined as the occurrence of any of the following [18]: (1) an increase in serum creatinine (SCr) by ≥ 0.3 mg/dL (26.5 µmol/L) within 48 h; (2) a known or presumed increase in SCr to 1.5 times the baseline within 7 days; (3) urine output < 0.5 mL/kg/h for 6 h. Respiratory failure was defined as a partial pressure of oxygen (PaO2) < 60 mmHg in ambient air, diagnosed according to the Berlin definition for ARDS [19]. Myocardial injury was defined as any elevation in cardiac troponin T, I or creatine kinase (CK)-MB. Heart failure met the criteria outlined in the 2016 European Society of Cardiology guidelines for the diagnosis and treatment of acute and chronic heart failure [20]. Shock was defined by low blood pressure requiring either vasopressor, inotropic support, or inadequate tissue perfusion [21, 22]. Coagulopathy was defined as a platelet count < 80 × 109/L or a 50% decrease in platelet count from the highest value within 3 days (for patients with chronic hematologic/tumor) or an international normalized ratio (INR) > 2. Liver dysfunction was defined as total bilirubin ≥ 4 mg/dL (not applicable for newborns) or ALT ≥ 2 times the upper limit of normal for age, in accordance with established criteria [23]. Accompanying symptoms were categorized into six groups based on the affected system: (1) renal injury: AKI; (2) respiratory system injury: respiratory failure, ARDS; (3) circulatory system injury: myocardial injury, heart failure, and shock; (4) coagulation dysfunction: DIC; (5) liver injury: liver dysfunction; and (6) immunological dysfunction: cytokine abnormalities and lymphocyte subgroup imbalances at the time of patient admission.

Statistical analysis

Categorical variables were presented as frequencies and percentages. Between-group comparisons of categorical variables were performed using the chi-square test or Fisher’s exact test, depending on data distribution. Logistic regression was used to estimate odds ratios (ORs) and 95% confidence intervals (CIs) for the association among EBV DNA levels, infection status, and mortality rate. Spearman’s rank correlation test was used to analyze the correlation between changes in EBV DNA levels, mortality rates, and organ damage. All hypothesis tests were two-tailed, with a type 1 error rate set at 5%. Multiple comparisons were adjusted using the Bonferroni method. All statistical analyses were performed using IBM SPSS Statistics for Windows version 25 (IBM Corp., Armonk, N.Y., USA) and Python software version 3.10 (Python Software Foundation, Beaverton, OR, USA), with data visualizations created using Python (Python Software Foundation) and GraphPad version 9.0 (San Diego, CA, USA).

To explore potential associations between immune indicators, prognosis, and infection status, we conducted clustering analysis on patients who underwent comprehensive immune parameter testing, including IFN-γ, IL-10, IL-6, IL-4, IL-2, and TNF-α levels; NK counts (CD3/CD16+CD56+); total B lymphocytes (CD3CD19+); total T lymphocytes (CD3+CD19); helper/inducible T lymphocytes (CD3+CD4+); inhibitory/CTLs (CD3+CD8+); and T lymphocytes + B lymphocytes + NK cells. For patients with complete immune indicators, a random split into training and validation sets was conducted at a ratio of 0.5: 0.5. Gaussian Mixture Model (GMM) clustering analysis was applied, with results visualized using principal component analysis. Clinical characteristics, infection status, and prognosis were compared among different clusters.

Results

Patients characteristics

This study included 35,956 pediatric patients (Fig. 1). The distribution of EBV DNA copy numbers was as follows: 90.34% had < 400 copies/mL, 1.84% had 400–1 × 103 copies/mL, 3.87% had 1 × 103 − 1 × 104 copies/mL, 2.12% had 1 × 104 − 1 × 105 copies/mL, 1.21% had 1 × 105 − 1 × 106 copies/mL, and 0.61% had > 1 × 106 copies/mL. The male-to-female ratio was 1.55: 1, and the median age at admission was 4.41 years (ranging: 0.02–17.98 years). Table 1 shows the characteristics of the study population according to EBV DNA copy numbers. The incidence of HLH and IM increased with higher EBV DNA loads (P < 0.001). Regarding prognosis, the distribution of patients with a high viral load (> 105 copies/mL) was higher in the mortality group compared with the survival group (P < 0.001). Patients with EBV DNA > 1 × 104 copies/mL had a higher pediatric intensive care unit admission rate than those with lower loads; however, the difference was not statistically significant (Table 1). The distributions of infection status in different age and sex groups are presented in Additional Table 1.

Fig. 1
figure 1

Study flow chart

Table 1 General characteristics of the population, N (%)

EBV DNA levels in different age groups and disease groups

Further analysis was conducted on patients with EBV DNA levels > 400 copies/mL. Patients aged 5–12 years and 12–18 years exhibited lower EBV DNA levels than those aged 1–5 years (P < 0.001) (Fig. 2A). Among the different disease groups, the CAEBV group had the highest EBV DNA levels, with significant statistical differences compared with the pneumonia and other diagnostic groups (P < 0.001). The HLH group also exhibited statistically statistical differences compared with the IM, pneumonia, and other diagnostic groups (P < 0.001) (Fig. 2B).

Fig. 2
figure 2

Distributions of EBV DNA load, infection status, and mortality. (A) EBV DNA levels in different age groups. 1–5 years group vs. 5–12 years group, P < 0.001; 1–5 years group vs. 12–18 years group, P < 0.001. (B) EBV DNA levels in different disease types. IM vs. HLH, P < 0.001; IM vs. pneumonia, P < 0.001; IM vs. others diagnosis, P < 0.001; sepsis vs. pneumonia, P < 0.001; sepsis vs. others diagnosis, P < 0.001; HLH vs. pneumonia, P < 0.001; HLH vs. others diagnosis, P < 0.001; pneumonia vs. CAEBV, P < 0.001; pneumonia vs. others diagnosis, P < 0.001; and CAEBV vs. others diagnosis, P < 0.001. (C) The proportion of primary infections and Activation-IgG+ across different disease types. (D) Mortality rates across different EBV DNA levels and infection states. (E) Violin plot of ALT levels at different levels of EBV DNA. (F) Violin plot of AST levels at different levels of EBV DNA

Distribution of infection status in different diseases

Figure 2C shows the distribution of infection status among the different disease groups. The proportion of activation-IgG+ decreased sequentially in CAEBV (100%), IM (50.53%), HLH (40.05%), sepsis (15.06%), pneumonia (13.37%), immunodeficiency (7.41%), other diagnosis (4.62%), and lymphoma (4.55%) (Fig. 2C). The pathological types of lymphoma are presented in Additional Table 2.

Table 2 Risk analysis of death by EBV DNA load and infection status

ALT and AST levels in patients with different EBV DNA loads

Significant differences in AST and ALT levels were observed between patients with different EBV DNA load (AST: H = 1,188.2, P < 0.001; ALT: H = 1,291.75, P < 0.001). Both AST and ALT levels increased with increasing EBV DNA levels. Stratified by EBV DNA levels (< 400, 400–1 × 10³, 1 × 10³ − 1 × 10⁴, 1 × 10⁴ − 1 × 10⁵, 1 × 10⁵ − 1 × 10⁶, > 1 × 10⁶ copies/mL), the median AST levels [M (P₂₅, P₇₅), U/L] were 29.0 (21.6, 42.8), 36.35 (25.7, 69.05), 40.1 (25.9, 75.5), 45.5 (29.4, 93.0), 57.25 (35.7, 129.55), and 88.4 (45.5, 170.3), while the ALT levels [M (P₂₅, P₇₅), U/L] were 18.3 (12.7, 30.7), 27.65 (15.12, 85.52), 29.0 (16.0, 95.41), 36.0 (17.3, 115.25), 61.55 (26.32, 183.0), and 108.8 (34.5, 243.0) (Fig. 2D, E).

Associations between EBV DNA load levels, infection status, and outcomes

Patients with EBV activation-IgG+ (1.45%) exhibited a higher mortality rate than those with a primary infection (0.62%). Compared with the primary infection group, the adjusted OR (adOR) for the activation-IgG+ group was 2.48 (95%CI: 1.04–5.89, P = 0.04). Children with EBV DNA levels of > 1 × 105 copies/mL had a significantly higher risk of mortality. Compared with the 400–1 × 103 copies/mL group, the adOR for the 1 × 105 − 1 × 106 copies/mL group was 11.02 (95%CI: 2.49–48.72, P = 0.002), while for the > 1 × 106 copies/mL group, the adOR was 12.53 (95%CI: 2.64–59.47, P = 0.002) (Table 2).

Regarding the correlation between EBV DNA load levels and mortality among patients with EBV DNA activation, a sharp increase in mortality was observed when the EBV DNA was > 1 × 105 copies/mL (Fig. 2F). Compared with patients with EBV DNA 400 − 10³ copies/mL, those with EBV DNA 1 × 10⁵ − 1 × 10⁶ copies/mL had an adOR of 10.53 (95% CI: 2.38–46.59, P = 0.002) (Table 2). Moreover, the increased rate of mortality in the activation-IgG⁺ group was higher than that in the primary infection group, with an adOR of 2.48 (95% CI: 1.04–5.89, P = 0.04) (Table 2).

Mortality rates and incidence of organ injury across different EBV DNA load groups and infection status

Among children with EBV DNA copy > 400 copies/mL, the incidence rates of liver injury, cardiovascular system injury, coagulation system injury, respiratory system injury and AKI were 18.83%, 4.38%, 4.32%, 2.16% and 0.92%, respectively. Spearman’s rank correlation trend test analysis indicated that the mortality rate and incidences of respiratory, cardiovascular, coagulation, and liver injuries were positively correlated with EBV DNA load levels (P < 0.05) (Fig. 3). Among children with primary infection, the mortality rate and coagulation system injury positively correlated with EBV DNA load (P < 0.05). Respiratory system, cardiovascular system, and liver injuries also showed an increasing trend with higher EBV DNA load levels, although the P-values were insignificant (P > 0.05). Among children with EBV activation-IgG+, the mortality rate and incidences of AKI, respiratory system injury, and coagulation system injury positively correlated with EBV DNA load levels (P < 0.05). Cardiovascular and liver injuries also showed an increasing trend with higher EBV DNA load levels, although the P-values were insignificant (P > 0.05) (Fig. 3).

Fig. 3
figure 3

Heatmap of mortality rates and organ injury across levels of EBV DNA and infection states. # The correlation coefficient of Spearman’s rank correlation trend test between the incidence rate and EBV DNA levels among patients with EBV DNA copy number > 400 copies/mL. * The P-value of Spearman’s rank correlation trend test between the incidence rate and EBV DNA levels among patients with EBV DNA copy number > 400 copies/mL

GMM clustering analysis of immunological indicators

This study included 323 patients with complete data on immunological indicators. These patients were split into a training set (50%) and a validation set (50%) for GMM clustering analysis of immunological indicators to explore infection status and prognosis in populations with different immune characteristics. Two clusters, Cluster 0 and 1, were identified in the training set. Table 3 summarizes the general characteristics of the two clusters. In Cluster 0 of the training set, the immune status presented elevated pro-inflammatory indicators (IFN-γ, IL-6) and anti-inflammatory indicator (IL-10) levels compared with Cluster 1. The immune cell counts in Cluster 0, including NK cells, total B lymphocytes, helper/inducer T lymphocytes, inhibitory/CTLs, total T lymphocytes, and T + B + NK cells, were lower than those in Cluster 1 (Fig. 4). The immune profiles of both clusters in the validation set resembled those in the training set (Fig. 4).

Table 3 General characteristics of different clusters
Fig. 4
figure 4

GMM clustering analysis of immune indicators, PCA visualization, and immune indicators among different clusters. (A) Line graph of standardized mean values for different immune indicators across different clusters. (B) PCA visualization of the training and validation sets. (C) Bar and scatter plots of immune indicators for different clusters in the training set and validation set. * P < 0.05, ** P < 0.01, *** P < 0.005, **** P < 0.001

Compared with Cluster 1, Cluster 0 demonstrated higher mortality rates in both the training and validation sets (Table 3). In the training set, the activation-IgG+ rates for Cluster 0 and 1 were 82.61% and 44.35%, respectively (P < 0.05), and the mortality rates were 23.91% and 0.87%, respectively (P < 0.05). In the validation set, the activation-IgG+ rates for Cluster 0 and 1 were 59.26% and 57.04%, respectively (P = 0.775), and the mortality rates were 37.04% and 2.22%, respectively (P < 0.05) (Table 3). In the training set, compared with the primary infection group, the activation-IgG+ group had an adOR of 2.48 (95% CI: 1.04–5.89) for developing the immune status of Cluster 1 (P < 0.05). No statistically significant differences were observed in this correlation in the validation set (Table 4). In both the training and validation sets, Cluster 0 was associated with a higher risk of death than Cluster 1 (P < 0.05, Table 5).

Table 4 Associations between clusters and infection status
Table 5 Associations between clusters and outcomes

In all infection status groups, Cluster 0 exhibited higher mortality rates than Cluster 1. In the primary infection group, the mortality rates for Cluster 0 and 1 were 18.18% and 0%, respectively (P = 0.016). In the group with unknown infection status, the mortality rates for Cluster 0 and 1 were 37.5% and 0%, respectively (P < 0.05). In the activation-IgG+ group, the mortality rates for Cluster 0 and 1 were 29.63% and 3.12%, respectively (P < 0.05) (Table 6).

Table 6 Associations between immune clusters and outcomes according to infection status

Discussion

Our study focuses on the correlation between EBV DNA load, infection status, and disease severity. We investigated the differences between infection states to characterize their distinct EBV DNA loads and clinical manifestations. Additionally, we explored the immune characteristics of the EBV DNA-activated pediatric patients, enhancing our understanding of the clinical presentation of pediatric EBV infection and providing insights for improved diagnosis and treatment.

EBV, also known as human herpesvirus 4, is a DNA virus in the herpesvirus family. Its life cycle includes latent and lytic phases. Following infection, EBV enters a latent state within memory B cells. However, under specific pathological conditions, it may reactivate, leading to cell lysis. EBV activation poses significant health risks. Our study, which includes the largest number of cases of activated EBV infection, observed a higher incidence of EBV infection in males compared with females, consistent with the findings of other researches [24]. However, EBV DNA copy numbers did not exhibit significant differences between sexes, suggesting comparable levels of DNA replication across sexes. In our study, the highest incidence of EBV activation occurred in children aged 1–5 years, corresponding to the peak age for IM in China, whereas in Western countries, primary EBV infection tends to occur later [3].

EBV is associated with various pediatric diseases. Primary infection often manifests as IM, while reactivation is more diverse and frequently linked to immunosuppression-related conditions such as sepsis, HLH, or chemotherapy [7]. Our findings indicate that patients with HLH and sepsis exhibit higher EBV DNA loads and increased rates of activation-IgG+. The mechanism of EBV reactivation is a complex process involving cellular immune responses to viral infection, host immune status, and other potential factors, including T-cell and NK-cell functional defects, macrophage activation, and impaired immune surveillance [25,26,27]. Specifically, pathological stimuli lead to abnormal proliferation of T cells and NK cells, resulting in elevated cytokine levels that promote the EBV reactivation. Inflammatory responses are crucial features of pediatric infectious diseases, such as sepsis and HLH. The release of inflammatory mediators, such as IL-6, creates a conducive microenvironment, allowing EBV to re-enter an active state. In this process, IL-6 activates pathways related to EBV, promoting viral replication [28]. In critically ill patients, prolonged pathological stress can impair cellular immune function [29, 30], and reduce CTL activity, leading to suppressed EBV clearance [31]. Additionally, an increase in regulatory T cells (Treg) may result in an insufficient response to EBV infection [32], further creating favorable conditions for EBV reactivation. Therefore, close monitoring for EBV reactivation is crucial in pediatric patients with conditions that affect immune regulation.

Reactivation of EBV reflects disruptions in host immune function and is associated with an unfavorable prognosis. Moreover, EBV DNA level serves as a key marker of EBV activity. Our study revealed a correlation between high EBV DNA load and in-hospital mortality, with a significant increase in the risk coefficient when the EBV DNA copy number increased to > 1 × 105 copies/mL. The in-hospital mortality rate showed a noticeable inflection point as EBV DNA levels increased, particularly among reactivated patients. Additionally, elevated EBV DNA levels were associated with organ injury, including the liver, coagulation system, circulatory system, respiratory system, and kidneys. Prior research suggests that EBV pathogenesis involves B, T, and NK cell activation, along with cytokine production [1]. During the initial infection, EBV first infects B cells, triggering the activation of CD8+ T cells. These activated CD8+ T cells destroy virus-infected cells through cytotoxic effects, leading to necrosis or apoptosis of infected B cells, ultimately leading to reduced B cell count and an increase in CD8+ T cells. Activated T-cells also release cytokines that cause clinical symptoms and injury to various organs [33,34,35].

To explore the relationship between immune status and post-EBV infection outcomes, we conducted unsupervised clustering analysis on EBV-activated cases using immune indicators. The analysis categorized the patients into two clusters. Based on the analysis of the immune status, EBV infection status, and prognostic characteristics of the two clusters, patients in Cluster 0 exhibited elevated levels of both pro-inflammatory and anti-inflammatory cytokines, such as IFN-γ, IL-10, and IL-6, compared with those in Cluster 1. However, Cluster 0 showed lower expressions of immune cells, including NK cells, total B lymphocytes, helper/inducer T lymphocytes, inhibitory/CTLs, total T lymphocytes and T lymphocytes + B lymphocytes + NK cells compared with Cluster 1. Notably, Cluster 0 had a higher mortality and activation-IgG+ rates, suggesting that the coexistence of high pro-inflammatory and anti-inflammatory states, along with low immune cell numbers, might indicate a critical condition and potentially unfavorable prognosis. Previous studies have suggested that the simultaneous elevation of pro- and anti-inflammatory indicators in response to a sustained immune response is indicative of compensatory anti-inflammatory response syndrome (CARS). Prolonged immune dysregulation may lead to a mixed anti-inflammatory response syndrome (MARS), and sustained immune dysregulation can further result in immune paralysis and apoptosis of immune cells [30, 36]. While excessive inflammatory factors can cause tissue damage, a strong anti-inflammatory response may hinder pathogen clearance, leading to worsening conditions. In this study, the immune profile of Cluster 0 resembled of the MARS state, with a higher proportion of reactivated children in Cluster 0 (pooled proportion, 73.97%) compared to Cluster 1 (pooled proportion, 51.2%). This suggests that reactivated patients may be more prone to this immune imbalance, which is characterized by a high pro-inflammatory and anti-inflammatory coexistence, along with low immune cell numbers. However, regardless of the primary infection or activation-IgG+, Cluster 0 exhibited higher mortality rates compared with Cluster 1. These findings highlight the importance of monitoring immune indicators in patients with activated EBV. Additionally, close monitoring is necessary for an immune state characterized by high pro-inflammatory and anti-inflammatory states, along with low immune cell numbers.

This study has some limitations. The sample size for comprehensive immune indicators was limited. Consequently, we could not dynamically assess the relationship between changes in EBV DNA copy number and the immune system, a facet that will be addressed in future research. Additionally, this study was conducted at a single-center design, necessitating further validation to confirm the generalizability of our conclusions.

Conclusion

In children with EBV activation-IgG+, an EBV DNA copy number exceeding 105 copies/mL was associated with a significantly increased mortality risk. Activation-IgG+ posed a higher mortality risk than primary infection, and an elevated EBV DNA copy number correlated with complications involving the liver, coagulation system, circulatory system, respiratory system, and AKI. In children with EBV activation-IgG+, a combination of high pro-inflammatory and anti-inflammatory activity, along with low immune cell numbers, may be indicative of a critical condition and poor prognosis.

Data availability

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

EBV:

Epstein-Barr virus

AKI:

Acute kidney injury

IM:

Infectious mononucleosis

HLH:

Lymphohistiocytosis

COVID-19:

Coronavirus disease 2019

NK:

Natural killer

PCR:

Polymerase chain reaction

VCA:

Viral capsid antigen, TNF, tumor necrosis factor

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

CAEBV:

Chronic active EBV

ARDS:

Acute respiratory distress syndrome

KDIGO:

Kidney Disease: Improving Global Outcomes

ORs:

Odds ratios

Cis:

Confidence intervals

GMM:

Gaussian mixture model

adOR:

Adjusted OR

CARS:

Compensatory anti-inflammatory response syndrome

MARS:

Mixed anti-inflammatory response syndrome

References

  1. Damania B, Kenney SC, Raab-Traub N. Epstein-Barr virus: biology and clinical disease. Cell. 2022;185(20):3652–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Lunn RM, Jahnke GD, Rabkin CS. Tumour virus epidemiology. Philosophical Trans Royal Soc B: Biol Sci 2017;372(1732).

  3. Chan KH, Dowd JB, Palermo T, Brite J, McDade TW, Aiello A. Seroprevalence of Epstein-Barr virus infection in U.S. Children ages 6–19, 2003–2010. PLoS ONE 2013;8(5).

  4. Ye Z, Chen L, Zhong H, Cao L, Fu P, Xu J. Epidemiology and clinical characteristics of Epstein-Barr virus infection among children in Shanghai, China, 2017–2022. Front Cell Infect Microbiol 2023;13.

  5. Lin W, Xiong G, Zhang B et al. Epstein-Barr virus (EBV) infection in Chinese children: A retrospective study of Age-Specific prevalence. PLoS ONE 2014;9(6).

  6. Shi T, Huang L, Chen Z, Tian J. Characteristics of primary Epstein-Barr virus infection disease spectrum and its reactivation in children, in Suzhou, China. J Med Virol. 2021;93(8):5048–57.

    Article  CAS  PubMed  Google Scholar 

  7. Yang Y, Gao F. Clinical characteristics of primary and reactivated Epstein-Barr virus infection in children. J Med Virol 2020.

  8. Taylor GS, Long HM, Brooks JM, Rickinson AB, Hislop AD. The immunology of Epstein-Barr Virus - Induced disease. Annu Rev Immunol. 2015;33(1):787–821.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang J, Qin S, Jin Z, et al. The clinical significance and prognostic role of Whole-Blood Epstein-Barr virus DNA in Lymphoma-Associated hemophagocytic lymphohistiocytosis. J Clin Immunol. 2023;43(6):1302–10.

    Article  CAS  PubMed  Google Scholar 

  10. Okamoto A, Yanada M, Inaguma Y, et al. The prognostic significance of EBV DNA load and EBER status in diagnostic specimens from diffuse large B-cell lymphoma patients. Hematol Oncol. 2015;35(1):87–93.

    Article  PubMed  Google Scholar 

  11. Li H, Cao D, Li S et al. Synergistic association of hepatitis B surface antigen and plasma Epstein-Barr virus DNA load on distant metastasis in patients with nasopharyngeal carcinoma. JAMA Netw Open 2023;6(2).

  12. Chan KCA, Woo JKS, King A, et al. Analysis of plasma Epstein - Barr virus DNA to screen for nasopharyngeal Cancer. N Engl J Med. 2017;377(6):513–22.

    Article  CAS  PubMed  Google Scholar 

  13. Qiu L, Si J, Kang J, et al. A retrospective analysis of EBV-DNA status with the prognosis of lymphoma. J Cell Mol Med. 2022;26(20):5195–201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Paolucci S, Cassaniti I, Novazzi F, et al. EBV DNA increase in COVID-19 patients with impaired lymphocyte subpopulation count. Int J Infect Dis. 2021;104:315–9.

    Article  CAS  PubMed  Google Scholar 

  15. Naughton P, Healy M, Enright F, Lucey B. Infectious mononucleosis: diagnosis and clinical interpretation. Br J Biomed Sci. 2021;78(3):107–16.

    Article  CAS  PubMed  Google Scholar 

  16. Okano M, Kawa K, Kimura H, et al. Proposed guidelines for diagnosing chronic active Epstein-Barr virus infection. Am J Hematol. 2005;80(1):64–9.

    Article  PubMed  Google Scholar 

  17. Chan JKC. The new world health organization classification of lymphomas: the past, the present and the future. Hematol Oncol. 2001;19(4):129–50.

    Article  CAS  PubMed  Google Scholar 

  18. Palevsky PM, Liu KD, Brophy PD, et al. KDOQI US commentary on the 2012 KDIGO clinical practice guideline for acute kidney injury. Am J Kidney Dis. 2013;61(5):649–72.

    Article  PubMed  Google Scholar 

  19. Force ADT, Ranieri VM, Rubenfeld GD, et al. Acute respiratory distress syndrome: the Berlin definition. JAMA. 2012;307(23):2526–33.

    Google Scholar 

  20. Ponikowski P, Voors AA, Anker SD, et al. 2016 ESC guidelines for the diagnosis and treatment of acute and chronic heart failure. Eur Heart J. 2016;37(27):2129–200.

    Article  PubMed  Google Scholar 

  21. Rhodes A, Evans LE, Alhazzani W, et al. Surviving Sepsis campaign: international guidelines for management of Sepsis and septic shock: 2016. Intensive Care Med. 2017;43(3):304–77.

    Article  PubMed  Google Scholar 

  22. Aneja RK, Carcillo JA. Differences between adult and pediatric septic shock. Minerva Anestesiol. 2011;77(10):986–92.

    CAS  PubMed  Google Scholar 

  23. Goldstein B, Giroir B, Randolph A. International pediatric sepsis consensus conference: definitions for sepsis and organ dysfunction in pediatrics*. Pediatr Crit Care Med. 2005;6(1):2–8.

    Article  PubMed  Google Scholar 

  24. Shi J, Ma W, Li W. Epidemiologic features of children with Epstein-Barr virus associated diseases in Hangzhou, China. J Med Virol. 2019;92(8):1277–82.

    Article  PubMed  Google Scholar 

  25. Ford M, Orlando E, Amengual JE. EBV reactivation and lymphomagenesis: more questions than answers. Curr Hematol Malig Rep. 2023;18(6):226–33.

    Article  PubMed  Google Scholar 

  26. Mertowska P, Mertowski S, Smolak K, et al. Could immune checkpoint disorders and EBV reactivation be connected in the development of hematological malignancies in immunodeficient patients?? Cancers. 2023;15:19.

    Article  Google Scholar 

  27. Murata T, Sugimoto A, Inagaki T et al. Molecular basis of Epstein - Barr virus latency establishment and lytic reactivation. Viruses 2021;13(12).

  28. Tao Q, Zhang G, Tsang CM et al. Enhanced IL-6/IL-6R signaling promotes growth and malignant properties in EBV-Infected premalignant and cancerous nasopharyngeal epithelial cells. PLoS ONE 2013;8(5).

  29. Duggal NA, Snelson C, Shaheen U, Pearce V, Lord JM. Innate and adaptive immune dysregulation in critically ill ICU patients. Sci Rep 2018;8(1).

  30. Liu D, Huang S-Y, Sun J-H et al. Sepsis-induced immunosuppression: mechanisms, diagnosis and current treatment options. Military Med Res 2022;9(1).

  31. Antia R, Hawkins JB, Delgado-Eckert E, Thorley-Lawson DA, Shapiro M. The cycle of EBV infection explains persistence, the sizes of the infected cell populations and which come under CTL regulation. PLoS Pathog 2013;9(10).

  32. Li J, Qian C-N, Zeng Y-X. Regulatory T cells and EBV associated malignancies. Int Immunopharmacol. 2009;9(5):590–2.

    Article  CAS  PubMed  Google Scholar 

  33. Freeman ML, Burkum CE, Jensen MK, Woodland DL, Blackman MA. γ-Herpesvirus reactivation differentially stimulates Epitope-Specific CD8 T cell responses. J Immunol. 2012;188(8):3812–9.

    Article  CAS  PubMed  Google Scholar 

  34. Farina A, Rosato E, York M, Gewurz BE, Trojanowska M, Farina GA. Innate immune modulation induced by EBV lytic infection promotes endothelial cell inflammation and vascular injury in scleroderma. Front Immunol 2021;12.

  35. Deng Y, Chatterjee B, Zens K, et al. CD27 is required for protective lytic EBV antigen - specific CD8 + T-cell expansion. Blood. 2021;137(23):3225–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Bone RC. Immunologic dissonance: a continuing evolution in our Understanding of the systemic inflammatory response syndrome (SIRS) and the multiple organ dysfunction syndrome (MODS). Ann Intern Med. 1996;125(8):680–7.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This study was supported by the Hunan Provincial Natural Science Foundation of China (No. 2022JJ40203, grant to HY; 2023JJ40351, grant to LXie), the National Natural Science Foundation of China (Young Scientists Fund, No. 82102285, grant to XLi), the Hunan Provincial Science and Technology Department Project (2020SK1014-3 and 2020SK2114, grant to XLu), the Hunan Provincial Key Laboratory of Emergency Medicine for Children (No. 2018TP1028, grant to ZX), and the Scientific Research Project of Hunan Provincial Health Commission (No. 202217012704, grant to HY; B202217018341, grant to LXie). The study sponsors have no role in the study design, data collection, data analysis, data interpretation, or writing of the report.

Author information

Authors and Affiliations

Authors

Contributions

HY and XLi contributed to the study’s conception, analysed and interpreted the data, and wrote the manuscript. ZX and XLu contributed to the study design, interpreted the data, and revised the manuscript. XHZ, TL, LX, LXie,XW, YY and LL performed chart reviews, interpreted the data, and revised the manuscript. ZX and XLu designed the study, interpreted the data, and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xiulan Lu or Zhenghui Xiao.

Ethics declarations

Ethics approval and consent to participate

The study protocol was reviewed and approved by the Medical Ethics Committee of the Hunan Children’s Hospital (HCHLL-2022-50). Informed consent was waved because of the retrospective design. The authors had no access to information that could identify individual participants during and after data collection.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yan, H., Li, X., Zeng, X. et al. Active EBV infection in children: associations between DNA load, infection status, immune status, and disease severity. Virol J 22, 113 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-025-02741-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s12985-025-02741-7

Keywords