BackgroundRespiratory syncytial virus (RSV) infection has been associated with the subsequent development of recurrent wheezing and asthma, although the mechanisms involved are still unknown. We investigate the role of epigenetics in the respiratory morbidity after infection by comparing methylation patterns from children who develop recurrent wheezing (RW-RSV), subsequent asthma (AS-RVS), and those experiencing complete recovery (CR-RSV).MethodsProspective, observational study of infants aged < 2 years with RSV respiratory infection admitted to hospital and followed-up after discharge for at least three years. According to their clinical course, patients were categorized into subgroups: RW-RSV (n = 36), AS-RSV (n = 9), and CR-RSV (n = 32). The DNA genome-wide methylation pattern was analyzed in whole blood samples, collected during the acute phase of the infection, using the Illumina Infinium Methylation EPIC BeadChip (850K CpG sites). Differences in methylation were determined through a linear regression model adjusted for age, gender and cell composition.ResultsPatients who developed respiratory sequelae showed a statistically significant higher proportion of NK and CD8T cells (inferred through a deconvolution approach) than those with complete recovery. We identified 5,097 significant differentially methylated positions (DMPs) when comparing RW-RSV and AS-RVS together against CR-RSV. Methylation profiles affect several genes involved in airway inflammation processes. The most significant DMPs were found to be hypomethylated in cases and therefore generally leading to overexpression of affected genes. The lead CpG position (cg24509398) falls at the gene body of EYA3 (P-value = 2.77×10^-10), a tyrosine phosphatase connected with pulmonary vascular remodeling, a key process in the asthma pathology. Logistic regression analysis resulted in a diagnostic epigenetic signature of 3-DMPs (involving genes ZNF2698, LOC102723354 and RPL15/NKIRAS1) that allows to efficiently differentiate sequelae cases from CR-RSV patients (AUC = 1.00). Enrichment pathway analysis reveals the role of the cell cycle checkpoint (FDR P-value = 4.71×10^-2), DNA damage (FDP-value = 2.53×10^-2), and DNA integrity checkpoint (FDR P-value = 2.56×10^-2) in differentiating sequelae from CR-RSV patients.ConclusionsEpigenetic mechanisms might play a fundamental role in the long-term sequelae after RSV infection, contributing to explain the different phenotypes observed.