Permutation-based significance analysis reduces the type 1 error rate in bisulfite sequencing data analysis of human umbilical cord blood samples

Essi Laajala*, Viivi Halla-aho, Toni Grönroos, Ubaid Ullah Kalim, Mari Vähä-Mäkilä, Mirja Nurmio, Henna Kallionpää, Niina Lietzén, Juha Mykkänen, Omid Rasool, Jorma Toppari, Matej Orešič, Mikael Knip, Riikka Lund, Riitta Lahesmaa, Harri Lähdesmäki*

*Corresponding author for this work

Research output: Contribution to journalArticleScientificpeer-review

2 Citations (Scopus)
31 Downloads (Pure)


DNA methylation patterns are largely established in-utero and might mediate the impacts of in-utero conditions on later health outcomes. Associations between perinatal DNA methylation marks and pregnancy-related variables, such as maternal age and gestational weight gain, have been earlier studied with methylation microarrays, which typically cover less than 2% of human CpG sites. To detect such associations outside these regions, we chose the bisulphite sequencing approach. We collected and curated clinical data on 200 newborn infants; whose umbilical cord blood samples were analysed with the reduced representation bisulphite sequencing (RRBS) method. A generalized linear mixed-effects model was fit for each high coverage CpG site, followed by spatial and multiple testing adjustment of P values to identify differentially methylated cytosines (DMCs) and regions (DMRs) associated with clinical variables, such as maternal age, mode of delivery, and birth weight. Type 1 error rate was then evaluated with a permutation analysis. We discovered a strong inflation of spatially adjusted P values through the permutation analysis, which we then applied for empirical type 1 error control. The inflation of P values was caused by a common method for spatial adjustment and DMR detection, implemented in tools comb-p and RADMeth. Based on empirically estimated significance thresholds, very little differential methylation was associated with any of the studied clinical variables, other than sex. With this analysis workflow, the sex-associated differentially methylated regions were highly reproducible across studies, technologies, and statistical models.

Original languageEnglish
Issue number12
Publication statusPublished - 2022
MoE publication typeA1 Journal article-refereed


  • analysis workflow
  • bisulphite sequencing
  • differential methylation
  • DNA methylation
  • pregnancy
  • RRBS
  • sex
  • spatial correlation
  • type 1 error
  • umbilical cord blood


Dive into the research topics of 'Permutation-based significance analysis reduces the type 1 error rate in bisulfite sequencing data analysis of human umbilical cord blood samples'. Together they form a unique fingerprint.

Cite this