Infant RSV immunoprophylaxis changes nasal epithelial DNA methylation at six years of age
=========================================================================================

* Cheng-Jian Xu
* Nienke M. Scheltema
* Cancan Qi
* Rolf Vedder
* Laura B.C. Klein
* Elisabeth E. Nibbelke
* Cornelis K van der Ent
* Louis J. Bont
* Gerard H. Koppelman

## Abstract

**Background** Respiratory syncytial virus (RSV) infection has been associated with childhood wheeze and asthma, and potential mechanisms include persistent epigenetic effects.

**Methods** In the randomized, placebo-controlled MAKI trial, 429 preterm infants randomly received RSV immunoprophylaxis with palivizumab or placebo during their first RSV season. Children were followed until age 6 for asthma evaluation. DNA methylation in cells obtained by nasal brushes at age 6 was measured by Illumina MethylationEPIC array.

**Results** RSV immunoprophylaxis in infancy had significant impact on global methylation patterns in nasal cells at age 6. The first principal component related to the immunoprophylaxis intervention was enriched for the pathway “positive regulation of defense response to virus by host” and “antigen processing and presentation” and driven by methylation changes in *NOD2, DGKG, MSH3*, and *ITPR2*. Three CpGsites, cg18040241, cg08243963, cg19555973 were differentially methylated at genome-wide significance, but were not associated with asthma. Differential methylation region analysis identified regions near genes that were previously implicated in the development of asthma and allergy such as *HLA-DPA1, HLA-DPB1, FASLG*, and *CHI3L1*.

**Conclusions** The study provides the first proof of concept that RSV immunoprophylaxis during infancy has long-term effects on nasal epigenetic signatures at age 6, relating to host antiviral defense pathways.

Keywords
*   RSV infection
*   epigenetics
*   RSV immunoprophylaxis
*   DNA methylation
*   asthma

## Background

Human Respiratory syncytial virus (RSV)-related acute lower respiratory tract infection is a leading cause of severe respiratory morbidity and mortality in children [1]. RSV infection in early life is associated with increased risk of early childhood wheezing [2][3][4][5][6]. Recurrent wheezing and asthma are thought to result from alterations in early life immune development following RSV infection, when the immune system is immature. Premature born children are at higher risk for the adverse effects of early life RSV lower respiratory tract infection. However, how RSV infection would have a long-term effect on the host immune system and how this would relate to the development of asthma remains unknown.

Epigenetics refers to (heritable) changes in gene expression that are not encoded by genetic variation. Recent epigenome-wide association studies (EWAS) of allergic diseases have reported encouraging results [7][8][9][10], and epigenetics might explain the high degree of plasticity of the immune response throughout life [11]. For example, in a recent EWAS, we reported 14 CpG methylation sites in whole blood were associated with childhood asthma, that related to distinct gene expression signatures reflecting early life activation of eosinophils and effector/memory CD8 T cells and NK cells at the expense of naïve T cell subsets, providing support for a role of early life antiviral immunity in asthma inception [7]. Indeed, RSV infection may have epigenetic consequences. RSV infection of Dendritic Cells leads to alterations in histone methylation by the H3K4 demethylase KDM5B, resulting in decreased pro-inflammatory cytokine production and a subsequent increase of Th2 cytokines from T cells[12]. Moreover, RSV has been related to a nasal airway microRNA profile which is predicted to have a downregulated NF-kB signaling pathway [13]. Among these epigenetic markers, DNA methylation provides a stable marker which is suitable for investigating mechanisms for long-term regulation of gene expression upon environmental stimulation [14]. Many studies relate early life exposures to DNA methylation signatures later in life, yet these observational studies may subject to confounding and causation remains to be verified by future studies. A possible solution to this would be to nest epigenetic analyses in a randomized controlled clinical trial, where the exposure can be controlled by the randomized intervention.

We hypothesized that RSV infection in early life causes changes in DNA methylation in the nasal airway epithelium in childhood. To test this hypothesis, we investigated the randomized clinical trial MAKI [2] (ISRCTN registry, number [ISRCTN73641710](http://medrxiv.org/external-ref?link_type=ISRCTN&access_num=ISRCTN73641710)) to investigate the direct link between RSV prophylaxis during infancy and DNA methylation changes in the nasal epithelium at age 6. In this trial, we observed significant protective effects of early life RSV prophylaxis on having RSV infection and asthma symptoms in the first years of life [2][3]. We performed an epigenome-wide association study (EWAS) to identify differentially methylated patterns in relation to palivizumab prevention at the global, individual and regional CpG site level.

## Method

### Study design and participants

In the MAKI trial, 429 otherwise healthy late preterm infants were enrolled between 2008 and 2010, and they were born at 33 to 35 weeks of gestation. Infants were randomly assigned to receive palivizumab (n=214) or placebo (n=215) during their first RSV season. Baseline characteristics at randomization are available at NEJM.org[2]. Details about the design, definitions, protocol of the primary study, follow-up study and clinical assessment have been previously described[2][3]. Clinical assessment at age 6 years included questionnaires for respiratory symptoms (parent-reported current asthma) as well as a parent report of a physician’ s diagnosis of asthma, lung function tests (FEV0.5, FVC), and assessment of specific IgE to a panel of aeroallergens as well as total serum IgE. Positive specific IgE was defined as an allergen-specific IgE concentration of at least 0·35 kU/L for house dust mite, birch, mugwort, timothy, Aspergillus, fumigatus, dog, and cat[3].

### Nasal DNA methylation measurements

Nasal epithelial cells were collected by nasal brushing around age 6 years of age at a home visit by a trained investigator. Briefly, the right nostril of the subjects was examined and the inferior turbinate was located using a speculum and penlight. The lateral area underneath the inferior turbinate was then brushed for 3 seconds with two brushes (Copan, 56380CS01 FLOQswabs) and these were placed in a 2 ml screw-cap Eppendorf tube and put into a freezer at – 80°C until further processing.

DNA was extracted from nasal brushes using the DNA investigator kit (Qiagen, Benelux BV, Venlo, the Netherlands). This was followed by precipitation-based purification and concentration using GlycoBlue (Ambion). 500 ng of DNA was bisulphite-converted using the EZ 96-DNA methylation kit (Zymo Research), following the manufacturer’s standard protocol. After verification of the bisulphite conversion step using Sanger Sequencing, DNA concentration was normalized and the samples were randomized to avoid batch effects. One standard DNA sample per chip was included in this step for quality control. In total, 296 nasal epithelium samples with sufficient DNA quality and quantity were hybridized to the Infinium HumanMethylationEPIC BeadChip array (Illumina, San Diego, CA, USA).

### Data pre-processing and statistical analysis

DNA methylation data were pre-processed in R with the Bioconductor package Minfi[15], using the original IDAT files extracted from the HiScanSQ scanner. We implemented sample filtering to remove 6 bad quality samples (call rate <99%) and 16 samples with gender mismatch. During processing, bad quality probes which failed more than 10% of the samples, the probes on sex chromosomes, the probes that mapped to multiple loci, and the probes containing SNPs at the target CpG sites with a MAF>5% in European populations were excluded14. We subsequently implemented stratified quantile normalization[17]. To remove the effect of extreme outliers in data, we trimmed the methylation set using: (25th percentile −3*IQR) and (75th percentile+3*IQR), where IQR=interquartile range. In our predefined analysis plan, we followed the trial design by relating the DNA methylation signatures at age 6 years to the palivizumab intervention. The association of DNA methylation with palivizumab intervention was conducted using robust linear regression corrected for sex, age at nasal brush sampling, maternal smoking, batch and surrogate variable analysis (SVA)[18] which represent unknown confounders such as cell type heterogeneity. Principal component analysis (PCA) and nonmetric multidimensional scaling (NMDS)[19] were utilized to check the global effects of palivizumab prophylaxis. R package vegan was used for nonmetric multidimensional scaling analysis (NMDS). Principal components associated with the intervention were analyzed, and the top 20 CpG sites contributing to the first two principal components were annotated using illumina Infinium MethylationEPIC Manifest File (version 1.0), with significant pathways detected by gene network ([www.genenetwork.nl](http://www.genenetwork.nl))[20]. For single CpG site analysis, adjustment for multiple testing was done using Benjamini & Hochberg method. Significant differentially methylated CpG site were considered when the P-value after False Discovery rate correction (FDR) was less than 0.05. For regional analysis, comb-p[21] was used to identify differentially methylated regions (DMRs) by using the summery statistics from models of single CpG site analysis. All statistical analyses were performed using the computing environment R (version 3.32).

## Results

### Sample characteristics

Figure 1 illustrated the flow chart of initial enrollment, follow-up and quality control of DNA methylation data in the MAKI study. After quality control, 274 (63.8%) of the collected samples and 790,437 probes remained for further analyses. The baseline characteristics of 274 participating children (132 palivizumab; 142 placebo), representing 63.9 % of the original study population, are presented in Table 1. Median age at the time of nasal brush collection was 6.6 years. Current asthma was defined as parent reported wheeze or asthma medication in the past 12 months[3] ; current allergic asthma was defined as having asthma and specific IgE to aeroallergens at age 6 years. Asthma prevalence was somewhat lower in the palivizumab group compared to the placebo group (14 vs 24%, *P*=0.065) (Table 1), and no differences were observed in lung function (FEV0.5) and total serum IgE levels, as reported previously in the full study population[3].

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T1)

Table 1. Sample characteristics from MAKI methylation study

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/04/09/2021.04.06.21254883/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/F1)

Figure 1. Flow chart of enrollment at age 0, follow up at age 6, and quality control of DNA methylation data in MAKI study.

### Differential methylation profiling analysis

After QC, 790,437 CpG sites were available for downstream analyses. The principal components analysis (PCA) plot (Figure 2a) and nonmetric multidimensional scaling (NMDS) plot (Figure 3) illustrated an overall methylation difference at age 6 years between the RSV immunoprophylaxis group and placebo group. The first two principal components (Pvaluepc1=8.1×10−3; Pvaluepc2=1.0×10−2) were significantly associated with RSV immunoprophylaxis (Table E1). NMDS confirmed this overall methylation difference between the two groups (adonis, R2=0.107, *P*=0.0055, with 1999 permutations) (Figure 3).

View this table:
[Table E1](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T2)

Table E1 The association of the first ten principal components with Infant RSV immunoprophylaxis.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/04/09/2021.04.06.21254883/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/F2)

Figure 2. 
Principal component analysis (PCA) of global methylation patterns. **A**, PCA plot of global methylation patterns in DNA obtained nasal epithelial brushing in relation to Palivizumab intervention. Two PCs were significantly associated with the intervention: PC1 (p=8.1×10−3) and PC2 (p=1.0×10−2). **B**, Gene set analysis of PC1 using Go-biological process databases. **C**, Gene network of genes that significantly contribute to PC1. These genes annotate to immune related genes that were differentially methylated after RSV immunoprophylaxis in the first year of life. There are two clusters (blue and green) in this network.

![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/04/09/2021.04.06.21254883/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/F3)

Figure 3: 
Non metric dimensional biplot. NMDS biplots depicting the individual overall methylation levels (data points, n=274), colored by treatment group: placebo (blue, n=132) and Palivizumab (red, n=142). The difference in overall methylation levels is significant (adonis, R2= 0.0195, P value= 0.0115, with 1999 permutation a) NMDS1-NMDS2 b) NMDS1-NMDS3.

We next performed an Epigenome Wide Association Analysis (EWAS), and three CpG sites (cg18040241, cg08243963 and cg19555973) which are localized in Galactosidase *Beta 1 Like 2* (*GLB1L2)* and *SC5DL* (encoding lathosterol oxidase involved in cholesterol biosynthesis) were genome wide significantly associated with RSV immunoprophylaxis versus placebo (Figure 4a, Table E2, *P*<0.05 after FDR correction). These genes have not been previously implicated in viral immune responses or asthma development. The lambda value was 1.123, which indicated no strong inflation in the model after SVA correction (Figure 4b). In differential methylation region (DMR) analysis, we identified 43 DMRs in relation to RSV immunoprophylaxis (Table E3).

View this table:
[Table E2](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T3)

Table E2 Associations between DNA methylation and RSV immunoprophylaxis with adjustment for surrogate variables.

View this table:
[Table E3](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T4)

Table E3 DMRs associated with RSV immunoprophylaxis at an SLK-Sidak P value of less than 0.05.

![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2021/04/09/2021.04.06.21254883/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/F4)

Figure 4. 
a) Manhattan plot of 790,437 CpG sites for differentially methylation with early life RSV prophylaxis by palivizumab. P values (y axis) correspond to the significance of the difference in methylation between prevention and no prevention subject. The red line correspond to the P value FDR=0.05 threshed. b) A quantile-quantile (QQ) plot displays the experimentally observed Pvalues (vertical axis) as well as the expected Pvalues of a null distribution (horizontal axis). The observed pvalue is obtained from the association between CpG sites methylation and RSV prophylaxis. The model to test the association is methylation ∼ RSV prophylaxis + age at nasal brush sample + sex + maternal smoking + batch, without correction of surrogate variables. The grey area represents the 95% concentration band.

### Pathway analysis

To identify pathways explaining these overall methylation differences after RSV immunoprophylaxis, we selected the top 20 CpG sites contributing most to the first two PCs (Table E4 and E5). We next performed pathway analysis by Gene Network[22] on genes located closely to these CpG sites. The pathway analysis of PC1 indicated enrichment of Go biological pathways involved in “positive regulation of defense response to virus by host” driven by *NOD2* (a known innate immunity gene), *DGKG, MSH3* and *ITPR2*, and “viral budding via host endosomal sorting complexes required for transport (ESCRT)” driven by *CANX, SEMA4D, ITPR2* and *ADK* (Figure 2b, Table E6). Based on gene network clustering, Figure 2c depicts two clusters in this gene network: “positive regulation of interleukin-1 beta secretion” and “antigen processing and presentation of peptide antigen via MHC class I” (Table E7 and E8). The pathway analysis on genes linked to PC2 showed relation to innate immune response genes (Table E9).

View this table:
[Table E4](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T5)

Table E4 The top 20 CpG sites that contributed most in loadings to the first principal component (PC1).

View this table:
[Table E5](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T6)

Table E5 The top 20 CpGs that contributed most in loadings to the first principal component (PC2).

View this table:
[Table E6:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T7)

Table E6: Gene set analysis of PC1 using Go-Biological process databases.

View this table:
[Table E7:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T8)

Table E7: Gene set analysis by using Go-Biological process databases of green cluster of Figure 2c.

View this table:
[Table E8:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T9)

Table E8: Gene set analysis by using Go-Biological process databases of blue cluster of Figure 2c.

View this table:
[Table E9:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T10)

Table E9: Gene set analysis of PC2 by using Go-Biological process databases.

### Link of RSV prophylaxis epigenetic markers to asthma and allergy

We next related PC1 of RSV immunoprophylaxis associated methylation to patient reported asthma at age 6 years, but did not observe an association (Point-Biserial correlation −0.077, *P*=0.20). The three CpG sites which were associated with RSV immunoprophylaxis were also not associated with patient reported asthma (Table E10). We linked RSV prevention related epigenetic markers to asthma and allergy associated CpG sites at regional level as many well-known asthma and allergy genes have been identified (Table E3). This differential methylation region analysis identified regions near genes that were previously implicated in the development of asthma and allergy such as *HLA-DPA1, HLA-DPB1, FASLG*, and *CHI3L1*. [24] We further tested 30 CpG sites which were reported to be associated with allergic asthma from a recent nasal epigenomics study[23], yet none of these were associated with RSV immunoprophylaxis (Table E11). However, we did replicate 16 of these 30 CpG sites associated with allergic asthma in nasal epithelium in our study (Table E12).

View this table:
[Table E10:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T11)

Table E10: Associations between 3 CpG sites associated with RSV immunoprophylaxis with patient reported asthma.

View this table:
[Table E11:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T12)

Table E11: 
Associations between 30 CpG sites, which are reported to associate with allergic asthma in a previous study (Forno et al) with RSV immunoprophylaxis.

View this table:
[Table E12:](http://medrxiv.org/content/early/2021/04/09/2021.04.06.21254883/T13)

Table E12: 
Associations between 30 CpG sites, which are reported to associate with allergic asthma with allergic asthma in MAKI data and 16 out of 29 could be replicated (P value<0.05).

## Discussion

This work shows that RSV prophylaxis by palivizumab at infancy has global and persistent effects on nasal DNA methylation patterns in childhood. We could observe such changes at global, single CpG site and regional level. Using pathway analysis, we linked these immunoprophylaxis associated global DNA methylation changes to viral response and viral budding genes. We propose that prevention of RSV respiratory tract infection in early childhood leads to different timing and severity of early life RSV respiratory tract infection, which in turn affects local DNA methylation of relevant innate immune response genes. This was further evidenced by differential methylation of cg16834953 in the *NOD2* gene, which functions as a viral pattern recognition receptor and participates in inducing antiviral signaling [24].

An important question is whether the epigenetic patterns we identified would lead to functional differences of these nasal cells, for example in response to other respiratory infections. An alternative explanation would be a direct effect of palivizumab on global gene methylation patterns, or microbiome changes in relation to early life RSV infections that in turn could affect nasal DNA methylation. To our knowledge, this is the first time that long-term methylation effects could be related to the timing of early life RSV infection. We performed our study within the framework of the randomized controlled MAKI trial, enabling us to infer causality rather than reporting an association, as can be done in observational cohorts. As we expect all children to have had RSV respiratory tract infections at the age of six years, we interpret our findings on changes in palivizumab related DNA methylation as an effect of timing and/or severity of early life RSV infection. However, this hypothesis needs further validation in future studies, focused on the timing and severity of early life RSV infection, preferentially in a longitudinal design. The low prevalence of documented RSV infection in the first year of life did not enable us to perform these analyses in this cohort (see table 1).

We did not identify a direct link between the epigenetic markers associated with RSV immunoprophylaxis and asthma at age 6 years in global and single site analyses, but identified some overlap in regional analysis. Although some asthma genes were identified, there was no strong overlap between the RSV immunoprophylaxis associated gene set and previously published asthma associated CpG sites. This is in line with the clinical findings in this cohort, reporting no strong association with early life palivizumab prevention and doctor’s diagnosed asthma at age 6 years[3]. Much larger studies are needed to identify a potential link between RSV immunoprophylaxis, epigenetic changes and asthma[25].

Despite the overall robustness and biological plausibility of our study findings, there are some limitations to consider. Firstly, to the best of our knowledge, there are no suitable data or cohorts available for replication of our findings. The study by O’Brien *et al*.[26] did not include epigenetic studies in the long-term follow up [personal communication]. Second, there is no gold standard to define cell composition in nasal brushes, so we used surrogates variables[18] to represent cell types. Third, our study was single blinded. Fourth, we could not directly confirm that these methylation changes in nasal epithelium regulated gene expression, which should be part of future studies. Fifth, although we studied over 800,00 CpG sites of the human methylome, this still represents about 3 % of all human CpG sites. Further work confirming these differentially methylated sites and detailed assessment of the associated regions using pyrosequencing should be performed. Sixth, other factors may also associate with DNA methylations, such as age[27], cigarette smoke[28], air pollution[29], or other respiratory infections. However, given that our analysis was embedded in a randomized controlled clinical trial, it is very reasonable to assume that potential confounders are distributed equally across the two treatment arms; and that the results were driven by the RSV Immunoprohylaxis.

In conclusion, infant RSV prophylaxis persistently changes nasal epithelial DNA methylation at least until age 6 years, and these methylation profiles are linked to regulation of anti-viral host defense responses and viral budding process.

## Data Availability

The epigenetic datasets generated and analyzed are available from the corresponding author on reasonable request.

## Authors’ contributions

The study was conceived and designed by LJB and GHK. CJX, CQ and RV performed statistical analyses, CJX drafted the manuscript and interpreted the data. NMS, LJB and GHK contributed to the manuscript writing. NMS, LBCK, EEN, CK and LJB obtained informed consent and clinical information. All authors reviewed and provided comments to the final manuscript (CJX, NMS, CQ, RV, LBCK, EEN, CK, LJB, GHK). All authors read and approved the final manuscript.

## Availability of data and materials

The epigenetic datasets generated and analyzed are available from the corresponding author on reasonable request.

## Ethics approval and consent to participate

The MAKI trial was approved by the ethical committee Medisch Ethische Toetsingscommissie (METC) of the Universitair Medisch Centrum Utrecht(21/07/2008), and all study participants provided written consent.

## Funding

This study was supported by MedImmune grant ESR-14-10006

## Potential conflict of interest

This study was supported by a grant from MedImmune to LJB and GHK. Funders had no role in design, interpretation and reporting of this study. L.J.B. reports grants from AbbVie during the conduct of the study and grants from MedImmune, Janssen, MeMed, and the Bill & Melinda Gates Foundation. None of the other authors report competing interests relating to the work presented in this manuscript.

## Acknowledgements

The authors are indebted to all the participating children and their families for their commitment and participation. We thank all members of the MAKI research team, including the staff of the pediatric lung function laboratory, and the laboratory staff. We also thank the staff of the Genomics Analysis Facility for providing technical assistance with the methylation analyses of this study.

## Abbreviations

RSV
:   Respiratory syncytial virus;
EWAS
:   Epigenome-wide association studies;
PCA
:   Principal components analysis;
NMDS
:   Nonmetric multidimensional scaling;
DMR
:   Differential methylation region;
ESCRT
:   Endosomal sorting complexes required for transport.

*   Received April 6, 2021.
*   Revision received April 6, 2021.
*   Accepted April 9, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Shi T, McAllister DA, O’Brien KL, Simoes EAF, Madhi SA, Gessner BD, Polack FP, Balsells E, Acacio S, Aguayo C, Alassani I, Ali A, Antonio M, Awasthi S, Awori JO, Azziz-Baumgartner E, Baggett HC, Baillie VL, Balmaseda A, Barahona A, Basnet S, Bassat Q, Basualdo W, Bigogo G, Bont L, Breiman RF, Brooks WA, Broor S, Bruce N, Bruden D, et al. Global, regional, and national disease burden estimates of acute lower respiratory infections due to respiratory syncytial virus in young children in 2015: a systematic review and modelling study. The Lancet 2017; 390: 946–958.
    
    

2.  2.Blanken MO, Rovers MM, Molenaar JM, Winkler-Seinstra PL, Meijer A, Kimpen JLL, Bont L. Respiratory Syncytial Virus and Recurrent Wheeze in Healthy Preterm Infants. N. Engl. J. Med. 2013; 368: 1791–1799.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1211917&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23656644&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000318540000007&link_type=ISI) 

3.  3.Scheltema NM, Nibbelke EE, Pouw J, Blanken MO, Rovers MM, Naaktgeboren CA, Mazur NI, Wildenbeest JG, van der Ent CK, Bont LJ. Respiratory syncytial virus prevention and asthma in healthy preterm infants: a randomised controlled trial. Lancet Respir. Med. 2018; 6: 257–264.
    
    

4.  4.Sigurs N, Bjarnason R, Sigurbergsson F, Kjellman B. Respiratory syncytial virus bronchiolitis in infancy is an important risk factor for asthma and allergy at age 7. Am. J. Respir. Crit. Care Med. 2000; 161: 1501–1507.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/ajrccm.161.5.9906076&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10806145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086945400019&link_type=ISI) 

5.  5.Stein RT, Sherrill D, Morgan WJ, Holberg CJ, Halonen M, Taussig LM, Wright AL, Martinez FD. Respiratory syncytial virus in early life and risk of wheeze and allergy by age 13 years. The Lancet 1999; 354: 541–545.
    
    

6.  6.Rosas-Salazar C, Shilts MH, Tovchigrechko A, Schobel S, Chappell JD, Larkin EK, Gebretsadik T, Halpin RA, Nelson KE, Moore ML, Anderson LJ, Peebles RS, Das SR, Hartert TV. Nasopharyngeal Lactobacillus is associated with a reduced risk of childhood wheezing illnesses following acute respiratory syncytial virus infection in infancy. J. Allergy Clin. Immunol. 2018; 142: 1447-1456.e9.
    
    

7.  7.Xu C-J, Söderhäll C, Bustamante M, Baïz N, Gruzieva O, Gehring U, Mason D, Chatzi L, Basterrechea M, Llop S, Torrent M, Forastiere F, Fantini MP, Carlsen KCL, Haahtela T, Morin A, Kerkhof M, Merid SK,  Rijkom B van, Jankipersadsing SA, Bonder MJ, Ballereau S, Vermeulen CJ, Aguirre-Gamboa R,  Jongste JC de, Smit HA, Kumar A, Pershagen G, Guerra S, Garcia-Aymerich J, et al. DNA methylation in childhood asthma: an epigenome-wide meta-analysis. Lancet Respir. Med. 2018; 6: 379–388.
    
    

8.  8.Liang L, Willis-Owen SAG, Laprise C, Wong KCC, Davies GA, Hudson TJ, Binia A, Hopkin JM, Yang IV, Grundberg E, Busche S, Hudson M, Rönnblom L, Pastinen TM, Schwartz DA, Lathrop GM, Moffatt MF, Cookson WOCM. An epigenome-wide association study of total serum immunoglobulin E concentration. Nature 2015; 520: 670–674.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14125&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25707804&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

9.  9.Reese SE, Xu C-J, den Dekker HT, Lee MK, Sikdar S, Ruiz-Arenas C, Merid SK, Rezwan FI, Page CM, Ullemar V, Melton PE, Oh SS, Yang IV, Burrows K, Söderhäll C, Jima DD, Gao L, Arathimos R, Küpers LK, Wielscher M, Rzehak P, Lahti J, Laprise C, Madore A-M, Ward J, Bennett BD, Wang T, Bell DA, Vonk JM, Håberg SE, et al. Epigenome-wide meta-analysis of DNA methylation and childhood asthma. J. Allergy Clin. Immunol. 2019; 143: 2062–2074.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2018.11.043&link_type=DOI) 

10. 10.Peng C, Cardenas A, Rifas-Shiman SL, Hivert M-F, Gold DR, Platts-Mills TA, Lin X, Oken E, Avila L, Celedón JC, Weiss ST, Baccarelli AA, Litonjua AA, DeMeo DL. Epigenetic age acceleration is associated with allergy and asthma in children in Project Viva. J. Allergy Clin. Immunol. 2019; 143: 2263-2270.e14.
    
    

11. 11.Tost J. A translational perspective on epigenetics in allergic diseases. J. Allergy Clin. Immunol. 2018; 142: 715–726.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2018.07.009&link_type=DOI) 

12. 12.Ptaschinski C, Mukherjee S, Moore ML, Albert M, Helin K, Kunkel SL, Lukacs NW. RSV-Induced H3K4 Demethylase KDM5B Leads to Regulation of Dendritic Cell-Derived Innate Cytokines and Exacerbates Pathogenesis In Vivo. PLOS Pathog. 2015; 11: e1004978.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1004978&link_type=DOI) 

13. 13.Hasegawa K, Pérez-Losada M, Hoptay CE, Epstein S, Mansbach JM, Teach SJ, Piedra PA, Camargo CA, Freishtat RJ. RSV vs. rhinovirus bronchiolitis: difference in nasal airway microRNA profiles and NFκB signaling. Pediatr. Res. 2018; 83: 606–614.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29244796&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

14. 14.Chen R, Xia L, Tu K, Duan M, Kukurba K, Li-Pook-Than J, Xie D, Snyder M. Longitudinal personal DNA methylome dynamics in a human with a chronic condition. Nat. Med. 2018; : 1.
    
    

15. 15.Fortin J-P,  Triche J Timothy J, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics 2017; 33: 558–560.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw691&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

16. 16.McCartney DL, Walker RM, Morris SW, McIntosh AM, Porteous DJ, Evans KL. Identification of polymorphic and off-target probe binding sites on the Illumina Infinium MethylationEPIC BeadChip. Genomics Data 2016; 9: 22–24.
    
    

17. 17.Touleimat N, Tost J. Complete pipeline for Infinium® Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics 2012; 4: 325–341.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2217/epi.12.21&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22690668&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000306530200011&link_type=ISI) 

18. 18.Leek JT, Storey JD. Capturing Heterogeneity in Gene Expression Studies by Surrogate Variable Analysis. PLOS Genet. 2007; 3: e161.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.0030161&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17907809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

19. 19.Kenkel NC, Orloci L. Applying Metric and Nonmetric Multidimensional Scaling to Ecological Studies: Some New Results. Ecology 1986; 67: 919–928.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/1939814&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986D369900010&link_type=ISI) 

20. 20.Deelen P, van Dam S, Herkert JC, Karjalainen JM, Brugge H, Abbott KM, van Diemen CC, van der Zwaag PA, Gerkes EH, Zonneveld-Huijssoon E, Boer-Bergsma JJ, Folkertsma P, Gillett T, van der Velde KJ, Kanninga R, van den Akker PC, Jan SZ, Hoorntje ET, te Rijdt WP, Vos YJ, Jongbloed JDH, van Ravenswaaij-Arts CMA, Sinke R, Sikkema-Raddatz B, Kerstjens-Frederikse WS, Swertz MA, Franke L. Improving the diagnostic yield of exome- sequencing by predicting gene–phenotype associations using large-scale gene expression analysis. Nat. Commun. 2019; 10: 2837.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-10649-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31253775&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

21. 21.Pedersen BS, Schwartz DA, Yang IV, Kechris KJ. Comb-p: software for combining, analyzing, grouping and correcting spatially correlated P-values. Bioinformatics 2012; 28: 2986–2988.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bts545&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22954632&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000311303500019&link_type=ISI) 

22. 22.Fehrmann RSN, Karjalainen JM, Krajewska M, Westra H-J, Maloney D, Simeonov A, Pers TH, Hirschhorn JN, Jansen RC, Schultes EA, van Haagen HHHBM, de Vries EGE, te Meerman GJ, Wijmenga C, van Vugt MATM, Franke L. Gene expression analysis identifies global gene dosage sensitivity in cancer. Nat. Genet. 2015; 47: 115.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3173&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25581432&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

23. 23.Forno E, Wang T, Qi C, Yan Q, Xu C-J, Boutaoui N, Han Y-Y, Weeks DE, Jiang Y, Rosser F, Vonk JM, Brouwer S, Acosta-Perez E, Colón-Semidey A, Alvarez M, Canino G, Koppelman GH, Chen W, Celedón JC. DNA methylation in nasal epithelium, atopy, and atopic asthma in children: a genome-wide study. Lancet Respir. Med. 2019; 7: 336–346.
    
    

24. 24.Sabbah A, Chang TH, Harnack R, Frohlich V, Tominaga K, Dube PH, Xiang Y, Bose S. Activation of innate immune antiviral responses by Nod2. Nat. Immunol. 2009; 10: 1073–1080.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ni.1782&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19701189&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000270006900010&link_type=ISI) 

25. 25.Riddell CA, Bhat N, Bont LJ, Dupont WD, Feikin DR, Fell DB, Gebretsadik T, Hartert TV, Hutcheon JA, Karron RA, Nair H, Reiner RC, Shi T, Sly PD, Stein RT, Wu P, Zar HJ, Ortiz JR. Informing randomized clinical trials of respiratory syncytial virus vaccination during pregnancy to prevent recurrent childhood wheezing: A sample size analysis. Vaccine 2018; 36: 8100–8109.
    
    

26. 26.O’Brien KL, Chandran A, Weatherholtz R, Jafri HS, Griffin MP, Bellamy T, Millar EV, Jensen KM, Harris BS, Reid R, Moulton LH, Losonsky GA, Karron RA, Santosham M. Efficacy of motavizumab for the prevention of respiratory syncytial virus disease in healthy Native American infants: a phase 3 randomised double-blind placebocontrolled trial. Lancet Infect. Dis. 2015; 15: 1398–1408.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(15)00247-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26511956&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

27. 27.Xu C-J, Bonder MJ, Söderhäll C, Bustamante M, Baïz N, Gehring U, Jankipersadsing SA, van der Vlies P, van Diemen CC, van Rijkom B, others. The emerging landscape of dynamic DNA methylation in early childhood. BMC Genomics 2017; 18: 25.
    
    

28. 28.Joubert BR, Felix JF, Yousefi P, Bakulski KM, Just AC, Breton C, Reese SE, Markunas CA, Richmond RC, Xu C-J, Küpers LK, Oh SS, Hoyo C, Gruzieva O, Söderhäll C, Salas LA, Baïz N, Zhang H, Lepeule J, Ruiz C, Ligthart S, Wang T, Taylor JA, Duijts L, Sharp GC, Jankipersadsing SA, Nilsen RM, Vaez A, Fallin MD, Hu D, et al. DNA Methylation in Newborns and Maternal Smoking in Pregnancy: Genome-wide Consortium Meta-analysis. Am. J. Hum. Genet. 2016; 98: 680–696.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2016.02.019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27040690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F09%2F2021.04.06.21254883.atom) 

29. 29.Gruzieva O, Xu C-J, Breton CV, Annesi-Maesano I, Antó JM, Auffray C, Ballereau S, Bellander T, Bousquet J, Bustamante M, Charles M-A, de Kluizenaar Y, den Dekker HT, Duijts L, Felix JF, Gehring U, Guxens M, Jaddoe VVW, Jankipersadsing SA, Merid SK, Kere J, Kumar A, Lemonnier N, Lepeule J, Nystad W, Page CM, Panasevich S, Postma D, Slama R, Sunyer J, et al. Epigenome-Wide Meta-Analysis of Methylation in Children Related to Prenatal NO2 Air Pollution Exposure. Environ. Health Perspect. 2017; 125: 104–110.