Search

Abstract Background A large number of analysis strategies are available for DNA methylation (DNAm) array and RNA-seq datasets, but it is unclear which strategies are best to use. We compare commonly used strategies and report how they influence results in large cohort studies. Results We tested the associations of DNAm and RNA expression with age, BMI, and smoking in four different cohorts (n = ~ 2900). By comparing strategies against the base model on the number and percentage of replicated CpGs for DNAm analyses or genes for RNA-seq analyses in a leave-one-out cohort replication approach, we find the choice of the normalization method and statistical test does not strongly influence the results for DNAm array data. However, adjusting for cell counts or hidden confounders substantially decreases the number of replicated CpGs for age and increases the number of replicated CpGs for BMI and smoking. For RNA-seq data, the choice of the normalization method, gene expression inclusion threshold, and statistical test does not strongly influence the results. Including five principal components or excluding correction of technical covariates or cell counts decreases the number of replicated genes. Conclusions Results were not influenced by the normalization method or statistical test. However, the correction method for cell counts, technical covariates, principal components, and/or hidden confounders does influence the results.

Abstract Background Cutaneous squamous cell carcinoma (cSCC) occurs 65–200 times more in immunosuppressed organ transplant patients than in the general population. T cells, which are targeted by the given immunosuppressive drugs, are involved in anti-tumor immune surveillance and are functionally regulated by DNA methylation. Prior to kidney transplantation, we aim to discover differentially methylated regions (DMRs) in T cells involved in de novo post-transplant cSCC development. Methods We matched 27 kidney transplant patients with a future de novo cSCC after transplantation to 27 kidney transplant patients without cSCC and studied genome-wide DNA methylation of T cells prior to transplantation. From 11 out of the 27 cSCC patients, the DNA methylation of T cells after transplantation was also examined to assess stability of the observed differences in DNA methylation. Raw methylation values obtained with the 450k array were confirmed with pyrosequencing. Results We found 16 DMRs between patients with a future cSCC and those who do not develop this complication after transplantation. The majority of the DMRs were located in regulatory genomic regions such as flanking bivalent transcription start sites and bivalent enhancer regions, and most of the DMRs contained CpG islands. Examples of genes annotated to the DMRs are ZNF577, coding for a zinc-finger protein, and FLOT1, coding for a protein involved in T cell migration. The longitudinal analysis revealed that DNA methylation of 9 DMRs changed significantly after transplantation. DNA methylation of 5 out of 16 DMRs was relatively stable, with a variation in beta-value lower than 0.05 for at least 50% of the CpG sites within that region. Conclusions This is the first study demonstrating that DNA methylation of T cells from patients with a future de novo post-transplant cSCC is different from patients without cSCC. These results were obtained before transplantation, a clinically relevant time point for cSCC risk assessment. Several DNA methylation profiles remained relatively stable after transplantation, concluding that these are minimally affected by the transplantation and possibly have a lasting effect on post-transplant cSCC development.