A MBD-seq protocol for large-scale methylome-wide studies with (very) low amounts of DNA

Insights from a methylome-wide association study of antidepressant exposure

This study tests the association of whole-blood DNA methylation and antidepressant exposure in 16,531 individuals from Generation Scotland (GS), using self-report and prescription-derived measures. We identify 8 associations and a high concordance of results between self-report and prescription-derived measures. Sex-stratified analyses observe nominally significant increased effect estimates in females for four CpGs. There is observed enrichment for genes expressed in the Amygdala and annotated to synaptic vesicle membrane ontology. Two CpGs (cg15071067; DGUOK-AS1 and cg26277237; KANK1) show correlation between DNA methylation with the time in treatment. There is a significant overlap in the top 1% of CpGs with another independent methylome-wide association study of antidepressant exposure. Finally, a methylation profile score trained on this sample shows a significant association with antidepressant exposure in a meta-analysis of eight independent external datasets. In this large investigation of antidepressant exposure and DNA methylation, we demonstrate robust associations which warrant further investigation to inform on the design of more effective and tolerated treatments for depression.

Cell-type specific methylation changes in the newborn child associated to obstetric pain relief

Although it is widely known that various pharmaceuticals affect the methylome, the knowledge of the effects from anesthesia is limited, and nearly nonexistent regarding the effects of obstetric anesthesia on the newborn child. Using sequencing based-methylation data and a reference-based statistical deconvolution approach we performed methylome-wide association studies (MWAS) of neonatal whole blood, and for each cell-type specifically, to detect methylation variations that are associated with the pain relief administered to the mother during delivery. Significant findings were replicated in a different dataset and followed-up with gene ontology analysis to pinpoint biological functions of potential relevance to these neonatal methylation alterations. The MWAS analyses detected methylome-wide significant (q<0.1) alterations in the newborn for laughing gas in granulocytes (two CpGs, p<5.50x10-9, q = 0.067), and for pudendal block in monocytes (five CpGs across three loci, p<1.51 x10-8, q = 0.073). Suggestively significant findings (p<1.00x10-6) were detected for both treatments for bulk and all cell-types, and replication analyses showed consistent significant enrichment (odds ratios ranging 3.47-39.02; p<4.00×10-4) for each treatment, suggesting our results are robust. In contrast, we did not observe any overlap across treatments, suggesting that the treatments are associated with different alterations of the neonatal blood methylome. Gene ontology analyses of the replicating suggestively significant results indicated functions related to, for example, cell differentiation, intracellular membrane-bound organelles and calcium transport. In conclusion, for the first time, we investigated and detected effect of obstetric pain-relief on the blood methylome in the newborn child. The observed differences suggest that anesthetic treatment, such as laughing gas or pudendal block, may alter the neonatal methylome in a cell-type specific manner. Some of the observed alterations are part of gene ontology terms that previously have been suggested in relation to anesthetic treatment, supporting its potential role also in obstetric anesthesia.

The Impact of Childhood Mental Health and Substance Use on Methylation Aging Into Adulthood

OBJECTIVE

To test whether childhood mental health symptoms, substance use, and early adversity accelerate the rate of DNA methylation (DNAm) aging from adolescence to adulthood.

METHOD

DNAm was assayed from blood samples in 381 participants in both adolescence (mean [SD] age = 13.9 [1.6] years) and adulthood (mean [SD] age = 25.9 [2.7] years). Structured diagnostic interviews were completed with participants and their parents at multiple childhood observations (1,950 total) to assess symptoms of common mental health disorders (attention-deficit/hyperactivity disorder, oppositional defiant disorder, conduct disorder, anxiety, and depression) and common types of substance use (alcohol, cannabis, nicotine) and early adversities.

RESULTS

Neither childhood mental health symptoms nor substance use variables were associated with DNAm aging cross-sectionally. In contrast, the following mental health symptoms and substance variables were associated with accelerated DNAm aging from adolescence to adulthood: depressive symptoms (b = 0.314, SE = 0.127, p = .014), internalizing symptoms (b = 0.108, SE = 0.049, p = .029), weekly cannabis use (b =1.665, SE = 0.591, p = .005), and years of weekly cannabis use (b = 0.718, SE = 0.283, p = .012). In models testing all individual variables simultaneously, the combined effect of the variables was equivalent to a potential difference of 3.17 to 3.76 years in DNAm aging. A final model tested a variable assessing cumulative exposure to mental health symptoms, substance use, and early adversities. This cumulative variable was strongly associated with accelerated aging (b = 0.126, SE = 0.044, p = .005).

CONCLUSION

Mental health symptoms and substance use accelerated DNAm aging into adulthood in a manner consistent with a shared risk mechanism.

PLAIN LANGUAGE SUMMARY

Using data from 381 participants in the Great Smoky Mountains Study, the authors examined whether childhood mental health symptoms, substance use, and early adversity accelerate biological aging, as measured by DNA methylation age, from adolescence to adulthood. Depressive symptoms and cannabis use were found to significantly accelerate biological aging. Models that tested the combined effect of mental health symptoms, substance use, and early adversity demonstrated that there was a shared effect across these types of childhood problems on accelerated aging.

DNA methyltransferase 1 (DNMT1) promotes cyst growth and epigenetic age acceleration in autosomal dominant polycystic kidney disease

Alteration of DNA methylation leads to diverse diseases, and the dynamic changes of DNA methylation (DNAm) on sets of CpG dinucleotides in mammalian genomes are termed "DNAm age" and "epigenetic clocks" that can predict chronological age. However, whether and how dysregulation of DNA methylation promotes cyst progression and epigenetic age acceleration in autosomal dominant polycystic kidney disease (ADPKD) remains elusive. Here, we show that DNA methyltransferase 1 (DNMT1) is upregulated in cystic kidney epithelial cells and tissues and that knockout of Dnmt1 and targeting DNMT1 with hydralazine, a safe demethylating agent, delays cyst growth in Pkd1 mutant kidneys and extends life span of Pkd1 conditional knockout mice. With methyl-CpG binding domain (MBD) protein-enriched genome sequencing (MBD-seq), DNMT1 chromatin immunoprecipitation (ChIP)-sequencing and RNA-sequencing analysis, we identified two novel DNMT1 targets, PTPRM and PTPN22 (members of the protein tyrosine phosphatase family). PTPRM and PTPN22 function as mediators of DNMT1 and the phosphorylation and activation of PKD-associated signaling pathways, including ERK, mTOR and STAT3. With whole-genome bisulfide sequencing in kidneys of patients with ADPKD versus normal individuals, we found that the methylation of epigenetic clock-associated genes was dysregulated, supporting that epigenetic age is accelerated in the kidneys of patients with ADPKD. Furthermore, five epigenetic clock-associated genes, including Hsd17b14, Itpkb, Mbnl1, Rassf5 and Plk2, were identified. Thus, the diverse biological roles of these five genes suggest that their methylation status may not only predict epigenetic age acceleration but also contribute to disease progression in ADPKD.

Investigating neonatal health risk variables through cell-type specific methylome-wide association studies

Adverse neonatal outcomes are a prevailing risk factor for both short- and long-term mortality and morbidity in infants. Given the importance of these outcomes, refining their assessment is paramount for improving prevention and care. Here we aim to enhance the assessment of these often correlated and multifaceted neonatal outcomes. To achieve this, we employ factor analysis to identify common and unique effects and further confirm these effects using criterion-related validity testing. This validation leverages methylome-wide profiles from neonatal blood. Specifically, we investigate nine neonatal health risk variables, including gestational age, Apgar score, three indicators of body size, jaundice, birth diagnosis, maternal preeclampsia, and maternal age. The methylomic profiles used for this research capture data from nearly all 28 million methylation sites in human blood, derived from the blood spot collected from 333 neonates, within 72 h post-birth. Our factor analysis revealed two common factors, size factor, that captured the shared effects of weight, head size, height, and gestational age and disease factor capturing the orthogonal shared effects of gestational age, combined with jaundice and birth diagnosis. To minimize false positives in the validation studies, validation was limited to variables with significant cumulative association as estimated through an in-sample replication procedure. This screening resulted in that the two common factors and the unique effects for gestational age, jaundice and Apgar were further investigated with full-scale cell-type specific methylome-wide association analyses. Highly significant, cell-type specific, associations were detected for both common effect factors and for Apgar. Gene Ontology analyses revealed multiple significant biologically relevant terms for the five fully investigated neonatal health risk variables. Given the established links between adverse neonatal outcomes and both immediate and long-term health, the distinct factor effects (representing the common and unique effects of the risk variables) and their biological profiles confirmed in our work, suggest their potential role as clinical biomarkers for assessing health risks and enhancing personalized care.

Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Methylome-wide association study of anxiety disorders

Anxiety Disorders (ANX) such as panic disorder, generalized anxiety disorder, and phobias, are highly prevalent conditions that are moderately heritable. Evidence suggests that DNA methylation may play a role, as it is involved in critical adaptations to changing environments. Applying an enrichment-based sequencing approach covering nearly 28 million autosomal CpG sites, we conducted a methylome-wide association study (MWAS) of lifetime ANX in 1132 participants (618 cases/514 controls) from the Netherlands Study of Depression and Anxiety. Using epigenomic deconvolution, we performed MWAS for the main cell types in blood: granulocytes, T-cells, B-cells and monocytes. Cell-type specific analyses identified 280 and 82 methylome-wide significant associations (q-value < 0.1) in monocytes and granulocytes, respectively. Our top finding in monocytes was located in ZNF823 on chromosome 19 (p = 1.38 × 10-10) previously associated with schizophrenia. We observed significant overlap (p < 1 × 10-06) with the same direction of effect in monocytes (210 sites), T-cells (135 sites), and B-cells (727 sites) between this Discovery MWAS signal and a comparable replication dataset from the Great Smoky Mountains Study (N = 433). Overlapping Discovery-Replication MWAS signal was enriched for findings from published GWAS of ANX, major depression, and post-traumatic stress disorder. In monocytes, two specific sites in the FZR1 gene showed significant replication after Bonferroni correction with an additional 15 nominally replicated sites in monocytes and 4 in T-cells. FZR1 regulates neurogenesis in the hippocampus, and its knockout leads to impairments in associative fear memory and long-term potentiation in mice. In the largest and most extensive methylome-wide study of ANX, we identified replicable methylation sites located in genes of potential relevance for brain mechanisms of psychiatric conditions.

Genes implicated by a methylome-wide schizophrenia study in neonatal blood show differential expression in adult brain samples

Schizophrenia is a disabling disorder involving genetic predisposition in combination with environmental influences that likely act via dynamic alterations of the epigenome and the transcriptome but its detailed pathophysiology is largely unknown. We performed cell-type specific methylome-wide association study of neonatal blood (N = 333) from individuals who later in life developed schizophrenia and controls. Suggestively significant associations (P < 1.0 × 10-6) were detected in all cell-types and in whole blood with methylome-wide significant associations in monocytes (P = 2.85 × 10-9-4.87 × 10-9), natural killer cells (P = 1.72 × 10-9-7.82 × 10-9) and B cells (P = 3.8 × 10-9). Validation of methylation findings in post-mortem brains (N = 596) from independent schizophrenia cases and controls showed significant enrichment of transcriptional differences (enrichment ratio = 1.98-3.23, P = 2.3 × 10-3-1.0 × 10-5), with specific highly significant differential expression for, for example, BDNF (t = -6.11, P = 1.90 × 10-9). In addition, expression difference in brain significantly predicted schizophrenia (multiple correlation = 0.15-0.22, P = 3.6 × 10-4-4.5 × 10-8). In summary, using a unique design combining pre-disease onset (neonatal) blood methylomic data and post-disease onset (post-mortem) brain transcriptional data, we have identified genes of likely functional relevance that are associated with schizophrenia susceptibility, rather than confounding disease associated artifacts. The identified loci may be of clinical value as a methylation-based biomarker for early detection of increased schizophrenia susceptibility.

Early adversities accelerate epigenetic aging into adulthood: a 10-year, within-subject analysis

BACKGROUND

Longitudinal studies are needed to clarify whether early adversities are associated with advanced methylation age or if they actually accelerate methylation aging. This study test whether different dimensions of childhood adversity accelerate biological aging from childhood to adulthood, and, if so, via which mechanisms.

METHODS

381 participants provided one blood sample in childhood (average age 15.0; SD = 2.3) and another in young adulthood (average age 23.1; SD = 2.8). Participants and their parents provided a median of 6 childhood assessments (total = 1,950 childhood observations), reporting exposures to different types of adversity dimensions (i.e. threat, material deprivation, loss, unpredictability). The blood samples were assayed to estimate DNA methylation age in both childhood and adulthood and also change in methylation age across this period.

RESULTS

Cross-sectional associations between the childhood adversity dimensions and childhood measures of methylation age were non-significant. In contrast, multiple adversity dimensions were associated with accelerated within-person change in methylation age from adolescence to young adulthood. These associations attenuated in model testing all dimensions at the same time. Accelerated aging increased with increasing number of childhood adversities: Individuals with highest number of adversities experienced 2+ additional years of methylation aging compared to those with no exposure to childhood adversities. The association between total childhood adversity exposure and accelerated aging was partially explained by childhood depressive symptoms, but not anxiety or behavioral symptoms.

CONCLUSIONS

Early adversities accelerate epigenetic aging long after they occur, in proportion to the total number of such experiences, and in a manner consistent with a shared effect that crosses multiple early dimensions of risk.

Editorial: Accelerated epigenetic ageing as a consequence of early environmental adversity

In this editorial, I discuss findings and implications of a hypothesis-generating study by Copeland and colleagues showing that early adverse events related to unpredictability as well as a cumulative score summarising environmental adversity in childhood and adolescence are linked to accelerated epigenetic ageing. The setting of this study is the longitudinal Great Smoky Mountains Study, a richly phenotyped longitudinal cohort with biological information. In my discussion, I focus on potential mechanisms behind the reported findings and opportunities for next steps, for example provided by the richness of data in the cohort used in this study.

Circulating Cell-Free DNA-Based Methylation Pattern in Saliva for Early Diagnosis of Head and Neck Cancer

Head and neck cancer (HNC) remains one of the leading causes of mortality worldwide due to tumor diagnosis at a late stage, loco-regional aggression, and distant metastases. A standardized diagnostic procedure for HNC is a tissue biopsy that cannot faithfully portray the in-depth tumor dynamics. Therefore, there is an urgent need to develop simple, accurate, and non-invasive methods for cancer detection and follow-up. A saliva-based liquid biopsy allows convenient, non-invasive, and painless collection of high volumes of this biofluid, with the possibility of repetitive sampling, all enabling real-time monitoring of the disease. No approved clinical test for HNC has yet been established. However, epigenetic changes in saliva circulating cell-free DNA (cfDNA) have the potential for a wide range of clinical applications. Therefore, the aim of this review is to present an overview of cfDNA-based methylation patterns in saliva for early detection of HNC, with particular attention to circulating tumor DNA (ctDNA). Due to advancements in isolation and detection technologies, as well as next- and third-generation sequencing, recent data suggest that salivary biomarkers may be successfully applied for early detection of HNC in the future, but large prospective clinical trials are still warranted.

DNA methylation signatures of childhood trauma predict psychiatric disorders and other adverse outcomes 17 years after exposure

Childhood trauma is robustly linked to a broad range of adverse outcomes with consequences persisting far into adulthood. We conducted a prospective longitudinal study to predict psychiatric disorders and other adverse outcomes from trauma-related methylation changes 16.9 years after trauma exposure in childhood. Methylation was assayed using a sequencing-based approach that provides near-complete coverage of all 28 million sites in the blood methylome. Methylation data involved 673 assays from 489 participants aged 13.6 years (SD = 1.9) with outcomes measures collected at age 30.4 (SD = 2.26). For a subset of 303 participants we also generated methylation data in adulthood. Trauma-related methylation risk scores (MRSs) significantly predicted adult depression, externalizing problems, nicotine dependence, alcohol use disorder, serious medical problems, social problems and poverty. The predictive power of the MRSs was higher than that of reported trauma and could not be explained by the reported trauma, correlations with demographic variables, or a continuity of the predicted health problems from childhood to adulthood. Rather than measuring the occurrence of traumatic events, the MRSs seemed to capture the subject-specific impact of trauma. The majority of predictive sites did not remain associated with the outcomes suggesting the signatures of trauma do not become biologically embedded in the blood methylome. Instead, the long-term effects of trauma therefore seemed more consistent with a developmental mechanism where the initial subject-specific impacts of trauma are magnified over time. The MRSs have the potential to be a novel clinical biomarker for the assessment of trauma-related health risks.

SEEMLIS: a flexible semi-automated method for enrichment of methylated DNA from low-input samples

Background

DNA methylation alterations have emerged as hallmarks of cancer and have been proposed as screening, prognostic, and predictive biomarkers. Traditional approaches for methylation analysis have relied on bisulfite conversion of DNA, which can damage DNA and is not suitable for targeted gene analysis in low-input samples. Here, we have adapted methyl-CpG-binding domain protein 2 (MBD2)-based DNA enrichment for use on a semi-automated exclusion-based sample preparation (ESP) platform for robust and scalable enrichment of methylated DNA from low-input samples, called SEEMLIS.

Results

We show that combining methylation-sensitive enzyme digestion with ESP-based MBD2 enrichment allows for single gene analysis with high sensitivity for GSTP1 in highly impure, heterogenous samples. We also show that ESP-based MBD2 enrichment coupled with targeted pre-amplification allows for analysis of multiple genes with sensitivities approaching the single cell level in pure samples for GSTP1 and RASSF1 and sensitivity down to 14 cells for these genes in highly impure samples. Finally, we demonstrate the potential clinical utility of SEEMLIS by successful detection of methylated gene signatures in circulating tumor cells (CTCs) from patients with prostate cancer with varying CTC number and sample purity.

Conclusions

SEEMLIS is a robust assay for targeted DNA methylation analysis in low-input samples, with flexibility at multiple steps. We demonstrate the feasibility of this assay to analyze DNA methylation in prostate cancer cells using CTCs from patients with prostate cancer as a real-world example of a low-input analyte of clinical importance. In summary, this novel assay provides a platform for determining methylation signatures in rare cell populations with broad implications for research as well as clinical applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-022-01252-4.

Dual methylation and hydroxymethylation study of alcohol use disorder

Using an integrative, multi-tissue design, we sought to characterize methylation and hydroxymethylation changes in blood and brain associated with alcohol use disorder (AUD). First, we used epigenomic deconvolution to perform cell-type-specific methylome-wide association studies within subpopulations of granulocytes/T-cells/B-cells/monocytes in 1132 blood samples. Blood findings were then examined for overlap with AUD-related associations with methylation and hydroxymethylation in 50 human post-mortem brain samples. Follow-up analyses investigated if overlapping findings mediated AUD-associated transcription changes in the same brain samples. Lastly, we replicated our blood findings in an independent sample of 412 individuals and aimed to replicate published alcohol methylation findings using our results. Cell-type-specific analyses in blood identified methylome-wide significant associations in monocytes and T-cells. The monocyte findings were significantly enriched for AUD-related methylation and hydroxymethylation in brain. Hydroxymethylation in specific sites mediated AUD-associated transcription in the same brain samples. As part of the most comprehensive methylation study of AUD to date, this work involved the first cell-type-specific methylation study of AUD conducted in blood, identifying and replicating a finding in DLGAP1 that may be a blood-based biomarker of AUD. In this first study to consider the role of hydroxymethylation in AUD, we found evidence for a novel mechanism for cognitive deficits associated with AUD. Our results suggest promising new avenues for AUD research.

DNA methylation of the KLK8 gene in depression symptomatology

Background

Depression is a common, complex, and debilitating mental disorder estimated to be under-diagnosed and insufficiently treated in society. Liability to depression is influenced by both genetic and environmental risk factors, which are both capable of impacting DNA methylation (DNAm). Accordingly, numerous studies have researched for DNAm signatures of this disorder. Recently, an epigenome-wide association study of monozygotic twins identified an association between DNAm status in the KLK8 (neuropsin) promoter region and severity of depression symptomatology.

Methods

In this study, we aimed to investigate: (i) if blood DNAm levels, quantified by pyrosequencing, at two CpG sites in the KLK8 promoter are associated with depression symptomatology and depression diagnosis in an independent clinical cohort and (ii) if KLK8 DNAm levels are associated with depression, postpartum depression, and depression symptomatology in four independent methylomic cohorts, with blood and brain DNAm quantified by either MBD-seq or 450 k methylation array.

Results

DNAm levels in KLK8 were not significantly different between depression cases and controls, and were not significantly associated with any of the depression symptomatology scores after correction for multiple testing (minimum p value for KLK8 CpG1 = 0.12 for ‘Depressed mood,’ and for CpG2 = 0.03 for ‘Loss of self-confidence with other people’). However, investigation of the link between KLK8 promoter DNAm levels and depression-related phenotypes collected from four methylomic cohorts identified significant association (p value < 0.05) between severity of depression symptomatology and blood DNAm levels at seven CpG sites.

Conclusions

Our findings suggest that variance in blood DNAm levels in KLK8 promoter region is associated with severity of depression symptoms, but not depression diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01184-5.

A targeted solution for estimating the cell-type composition of bulk samples

Background

To avoid false-positive findings and detect cell-type specific associations in methylation and transcription investigations with bulk samples, it is critical to know the proportions of the major cell-types.

Results

We present a novel approach that allows for precise estimation of cell-type proportions using only a few highly informative methylation markers. The most reliable estimates were obtained with 17 amplicons (34 CpGs) using the MuSiC estimator, for which the average correlations between the estimated and the true cell-type proportions were 0.889. Furthermore, the estimates were not significantly different from the true values (P = 0.95) indicating that the estimator is unbiased and the standard deviation of the estimates further indicate high precision. Moreover, the overall variability of the estimates as measured by the Root Mean Squared Error (RMSE), which is a function of both bias and precision, was low (mean RMSE = 0.038). Taken together, these results indicate that the approach produced reliable estimates that are both unbiased and highly precise.

Conclusion

This cost-effective approach for estimating cell-type proportions in bulk samples allows for enhanced targeted analysis, which in turn will minimize the risk of reporting false-positive findings and allowing for detection of cell-type specific associations. The approach is applicable across platforms and can be extended to assess cell-type proportions for various tissues.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04385-0.

The Detection of Cancer Epigenetic Traces in Cell-Free DNA

Nucleic acid fragments found in blood circulation originate mostly from dying cells and carry signs pointing to specific features of the parental cell types. Deciphering these clues may be transformative for numerous research and clinical applications but strongly depends on the development and implementation of robust analytical methods. Remarkable progress has been achieved in the reliable detection of sequence alterations in cell-free DNA while decoding epigenetic information from methylation and fragmentation patterns requires more sophisticated approaches. This review discusses the currently available strategies for detecting and analyzing the epigenetic marks in the liquid biopsies.

Cell-free DNA methylome profiling by MBD-seq with ultra-low input

Methylation signatures in cell-free DNA (cfDNA) have shown great sensitivity and specificity in the characterization of tumour status and classification of tumour types, as well as the response to therapy and recurrence. Currently, most cfDNA methylation studies are based on bisulphite conversion, especially targeted bisulphite sequencing, while enrichment-based methods such as cfMeDIP-seq are beginning to show potential. Here, we report an enrichment-based ultra-low input cfDNA methylation profiling method using methyl-CpG binding proteins capture, termed cfMBD-seq. We optimized the conditions for cfMBD capture by adjusting the amount of MethylCap protein along with using methylated filler DNA. Our data show high correlation between low input cfMBD-seq and standard MBD-seq (>1000 ng input). When compared to cfMEDIP-seq, cfMBD-seq demonstrates higher sequencing data quality with more sequenced reads passed filter and less duplicate rate. cfMBD-seq also outperforms cfMeDIP-seq in the enrichment of CpG islands. This new bisulphite-free ultra-low input methylation profiling technology has great potential in non-invasive and cost-effective cancer detection and classification.

Liquid Biopsy of Methylation Biomarkers in Cell-Free DNA

Liquid biopsies, in particular, analysis of cell-free DNA (cfDNA), have emerged as a promising noninvasive diagnostic approach in oncology. Abnormal distribution of DNA methylation is one of the hallmarks of many cancers and methylation changes occur early during carcinogenesis. Systemic analysis of cfDNA methylation profiles is being developed for cancer early detection, monitoring for minimal residual disease (MRD), predicting treatment response and prognosis, and tracing the tissue origin. This review highlights the advantages and disadvantages of ctDNA profiling for noninvasive diagnosis of early-stage cancers and explores recent advances in the clinical application of ctDNA methylation assays. We also summarize the technologies for ctDNA methylation analysis and provide a brief overview of the bioinformatic approaches for analyzing DNA methylation sequencing data.

Cell-Free DNA-Methylation-Based Methods and Applications in Oncology

Liquid biopsy based on cell-free DNA (cfDNA) enables non-invasive dynamic assessment of disease status in patients with cancer, both in the early and advanced settings. The analysis of DNA-methylation (DNAm) from cfDNA samples holds great promise due to the intrinsic characteristics of DNAm being more prevalent, pervasive, and cell- and tumor-type specific than genomics, for which established cfDNA assays already exist. Herein, we report on recent advances on experimental strategies for the analysis of DNAm in cfDNA samples. We describe the main steps of DNAm-based analysis workflows, including pre-analytics of cfDNA samples, DNA treatment, assays for DNAm evaluation, and methods for data analysis. We report on protocols, biomolecular techniques, and computational strategies enabling DNAm evaluation in the context of cfDNA analysis, along with practical considerations on input sample requirements and costs. We provide an overview on existing studies exploiting cell-free DNAm biomarkers for the detection and monitoring of cancer in early and advanced settings, for the evaluation of drug resistance, and for the identification of the cell-of-origin of tumors. Finally, we report on DNAm-based tests approved for clinical use and summarize their performance in the context of liquid biopsy.

The story of the lost twins: decoding the genetic identities of the Kumhar and Kurcha populations from the Indian subcontinent

Background

The population structure of the Indian subcontinent is a tapestry of extraordinary diversity characterized by the amalgamation of autochthonous and immigrant ancestries and rigid enforcement of sociocultural stratification. Here we investigated the genetic origin and population history of the Kumhars, a group of people who inhabit large parts of northern India. We compared 27 previously published Kumhar SNP genotype data sampled from Uttar Pradesh in north India to various modern day and ancient populations.

Results

Various approaches such as Principal Component Analysis (PCA), Admixture, TreeMix concurred that Kumhars have high ASI ancestry, minimal Steppe component and high genomic proximity to the Kurchas, a small and relatively little-known population found ~ 2500 km away in Kerala, south India. Given the same, biogeographical mapping using Geographic Population Structure (GPS) assigned most Kumhar samples in areas neighboring to those where Kurchas are found in south India.

Conclusions

We hypothesize that the significant genomic similarity between two apparently distinct modern-day Indian populations that inhabit well separated geographical areas with no known overlapping history or links, likely alludes to their common origin during or post the decline of the Indus Valley Civilization (estimated by ALDER). Thereafter, while they dispersed towards opposite ends of the Indian subcontinent, their genomic integrity and likeness remained preserved due to endogamous social practices. Our findings illuminate the genomic history of two Indian populations, allowing a glimpse into one or few of numerous of human migrations that likely occurred across the Indian subcontinent and contributed to shape its varied and vibrant evolutionary past.

A methylation study of long-term depression risk

Recurrent and chronic major depressive disorder (MDD) accounts for a substantial part of the disease burden because this course is most prevalent and typically requires long-term treatment. We associated blood DNA methylation profiles from 581 MDD patients at baseline with MDD status 6 years later. A resampling approach showed a highly significant association between methylation profiles in blood at baseline and future disease status (P = 2.0 × 10^-16). Top MWAS results were enriched specific pathways, overlapped with genes found in GWAS of MDD disease status, autoimmune disease and inflammation, and co-localized with eQTLS and (genic enhancers of) of transcription sites in brain and blood. Many of these findings remained significant after correction for multiple testing. The major themes emerging were cellular responses to stress and signaling mechanisms linked to immune cell migration and inflammation. This suggests that an immune signature of treatment-resistant depression is already present at baseline. We also created a methylation risk score (MRS) to predict MDD status 6 years later. The AUC of our MRS was 0.724 and higher than risk scores created using a set of five putative MDD biomarkers, genome-wide SNP data, and 27 clinical, demographic and lifestyle variables. Although further studies are needed to examine the generalizability to different patient populations, these results suggest that methylation profiles in blood may present a promising avenue to support clinical decision making by providing empirical information about the likelihood MDD is chronic or will recur in the future.

Methylome-wide association findings for major depressive disorder overlap in blood and brain and replicate in independent brain samples

We present the first large-scale methylome-wide association studies (MWAS) for major depressive disorder (MDD) to identify sites of potential importance for MDD etiology. Using a sequencing-based approach that provides near-complete coverage of all 28 million common CpGs in the human genome, we assay methylation in MDD cases and controls from both blood (N = 1132) and postmortem brain tissues (N = 61 samples from Brodmann Area 10, BA10). The MWAS for blood identified several loci with P ranging from 1.91 × 10^-8 to 4.39 × 10^-8 and a resampling approach showed that the cumulative association was significant (P = 4.03 × 10^-10) with the signal coming from the top 25,000 MWAS markers. Furthermore, a permutation-based analysis showed significant overlap (P = 5.4 × 10^-3) between the MWAS findings in blood and brain (BA10). This overlap was significantly enriched for a number of features including being in eQTLs in blood and the frontal cortex, CpG islands and shores, and exons. The overlapping sites were also enriched for active chromatin states in brain including genic enhancers and active transcription start sites. Furthermore, three loci located in GABBR2, RUFY3, and in an intergenic region on chromosome 2 replicated with the same direction of effect in the second brain tissue (BA25, N = 60) from the same individuals and in two independent brain collections (BA10, N = 81 and 64). GABBR2 inhibits neuronal activity through G protein-coupled second-messenger systems and RUFY3 is implicated in the establishment of neuronal polarity and axon elongation. In conclusion, we identified and replicated methylated loci associated with MDD that are involved in biological functions of likely importance to MDD etiology.

The DNA hypermethylation phenotype of colorectal cancer liver metastases resembles that of the primary colorectal cancers

BACKGROUND

Identifying molecular differences between primary and metastatic colorectal cancers-now possible with the aid of omics technologies-can improve our understanding of the biological mechanisms of cancer progression and facilitate the discovery of novel treatments for late-stage cancer. We compared the DNA methylomes of primary colorectal cancers (CRCs) and CRC metastases to the liver. Laser microdissection was used to obtain epithelial tissue (10 to 25 × 10⁶ μm²) from sections of fresh-frozen samples of primary CRCs (n = 6), CRC liver metastases (n = 12), and normal colon mucosa (n = 3). DNA extracted from tissues was enriched for methylated sequences with a methylCpG binding domain (MBD) polypeptide-based protocol and subjected to deep sequencing. The performance of this protocol was compared with that of targeted enrichment for bisulfite sequencing used in a previous study of ours.

RESULTS

MBD enrichment captured a total of 322,551 genomic regions (249.5 Mb or ~ 7.8% of the human genome), which included over seven million CpG sites. A few of these regions were differentially methylated at an expected false discovery rate (FDR) of 5% in neoplastic tissues (primaries: 0.67%, i.e., 2155 regions containing 279,441 CpG sites; liver metastases: 1%, i.e., 3223 regions containing 312,723 CpG sites) as compared with normal mucosa samples. Most of the differentially methylated regions (DMRs; 94% in primaries; 70% in metastases) were hypermethylated, and almost 80% of these (1882 of 2396) were present in both lesion types. At 5% FDR, no DMRs were detected in liver metastases vs. primary CRC. However, short regions of low-magnitude hypomethylation were frequent in metastases but rare in primaries. Hypermethylated DMRs were far more abundant in sequences classified as intragenic, gene-regulatory, or CpG shelves-shores-island segments, whereas hypomethylated DMRs were equally represented in extragenic (mainly, open-sea) and intragenic (mainly, gene bodies) sequences of the genome. Compared with targeted enrichment, MBD capture provided a better picture of the extension of CRC-associated DNA hypermethylation but was less powerful for identifying hypomethylation.

CONCLUSIONS

Our findings demonstrate that the hypermethylation phenotype in CRC liver metastases remains similar to that of the primary tumor, whereas CRC-associated DNA hypomethylation probably undergoes further progression after the cancer cells have migrated to the liver.

Cell Type-Specific Methylome-wide Association Studies Implicate Neurotrophin and Innate Immune Signaling in Major Depressive Disorder

BACKGROUND

We sought to characterize methylation changes in brain and blood associated with major depressive disorder (MDD). As analyses of bulk tissue may obscure association signals and hamper the biological interpretation of findings, these changes were studied on a cell type-specific level.

METHODS

In 3 collections of human postmortem brain (n = 206) and 1 collection of blood samples (N = 1132) of MDD cases and controls, we used epigenomic deconvolution to perform cell type-specific methylome-wide association studies within subpopulations of neurons/glia for the brain data and granulocytes/T cells/B cells/monocytes for the blood data. Sorted neurons/glia from a fourth postmortem brain collection (n = 58) were used for validation purposes.

RESULTS

Cell type-specific methylome-wide association studies identified multiple findings in neurons/glia that were detected across brain collections and were reproducible in physically sorted nuclei. Cell type-specific analyses in blood samples identified methylome-wide significant associations in T cells, monocytes, and whole blood that replicated findings from a past methylation study of MDD. Pathway analyses implicated p75 neurotrophin receptor/nerve growth factor signaling and innate immune toll-like receptor signaling in MDD. Top results in neurons, glia, bulk brain, T cells, monocytes, and whole blood were enriched for genes supported by genome-wide association studies for MDD and other psychiatric disorders.

CONCLUSIONS

We both replicated and identified novel MDD-methylation associations in human brain and blood samples at a cell type-specific level. Our results provide mechanistic insights into how the immune system may interact with the brain to affect MDD susceptibility. Importantly, our findings involved associations with MDD in human samples that implicated many closely related biological pathways. These disease-linked sites and pathways represent promising new therapeutic targets for MDD.

Independent Methylome-Wide Association Studies of Schizophrenia Detect Consistent Case-Control Differences

Methylome-wide association studies (MWASs) are promising complements to sequence variation studies. We used existing sequencing-based methylation data, which assayed the majority of all 28 million CpGs in the human genome, to perform an MWAS for schizophrenia in blood, while controlling for cell-type heterogeneity with a recently generated platform-specific reference panel. Next, we compared the MWAS results with findings from 3 existing large-scale array-based schizophrenia methylation studies in blood that assayed up to ~450 000 CpGs. Our MWAS identified 22 highly significant loci (P < 5 × 10-8) and 852 suggestively significant loci (P < 1 × 10-5). The top finding (P = 5.62 × 10-11, q = 0.001) was located in MFN2, which encodes mitofusin-2 that regulates Ca2+ transfer from the endoplasmic reticulum to mitochondria in cooperation with DISC1. The second-most significant site (P = 1.38 × 10-9, q = 0.013) was located in ALDH1A2, which encodes an enzyme for astrocyte-derived retinoic acid-a key neuronal morphogen with relevance for schizophrenia. Although the most significant MWAS findings were not assayed on the arrays, we observed significant enrichment of overlapping findings with 2 of the 3 array datasets (P = 0.0315, 0.0045, 0.1946). Overrepresentation analysis of Gene Ontology terms for the genes in the significant overlaps suggested high similarity in the biological functions detected by the different datasets. Top terms were related to immune and/or stress responses, cell adhesion and motility, and a broad range of processes essential for neurodevelopment.

Performance of Mapping Approaches for Whole-Genome Bisulfite Sequencing Data in Crop Plants

DNA methylation is involved in many different biological processes in the development and well-being of crop plants such as transposon activation, heterosis, environment-dependent transcriptome plasticity, aging, and many diseases. Whole-genome bisulfite sequencing is an excellent technology for detecting and quantifying DNA methylation patterns in a wide variety of species, but optimized data analysis pipelines exist only for a small number of species and are missing for many important crop plants. This is especially important as most existing benchmark studies have been performed on mammals with hardly any repetitive elements and without CHG and CHH methylation. Pipelines for the analysis of whole-genome bisulfite sequencing data usually consists of four steps: read trimming, read mapping, quantification of methylation levels, and prediction of differentially methylated regions (DMRs). Here we focus on read mapping, which is challenging because un-methylated cytosines are transformed to uracil during bisulfite treatment and to thymine during the subsequent polymerase chain reaction, and read mappers must be capable of dealing with this cytosine/thymine polymorphism. Several read mappers have been developed over the last years, with different strengths and weaknesses, but their performances have not been critically evaluated. Here, we compare eight read mappers: Bismark, BismarkBwt2, BSMAP, BS-Seeker2, Bwameth, GEM3, Segemehl, and GSNAP to assess the impact of the read-mapping results on the prediction of DMRs. We used simulated data generated from the genomes of Arabidopsis thaliana, Brassica napus, Glycine max, Solanum tuberosum, and Zea mays, monitored the effects of the bisulfite conversion rate, the sequencing error rate, the maximum number of allowed mismatches, as well as the genome structure and size, and calculated precision, number of uniquely mapped reads, distribution of the mapped reads, run time, and memory consumption as features for benchmarking the eight read mappers mentioned above. Furthermore, we validated our findings using real-world data of Glycine max and showed the influence of the mapping step on DMR calling in WGBS pipelines. We found that the conversion rate had only a minor impact on the mapping quality and the number of uniquely mapped reads, whereas the error rate and the maximum number of allowed mismatches had a strong impact and leads to differences of the performance of the eight read mappers. In conclusion, we recommend BSMAP which needs the shortest run time and yields the highest precision, and Bismark which requires the smallest amount of memory and yields precision and high numbers of uniquely mapped reads.

MBD-seq - realities of a misunderstood method for high-quality methylome-wide association studies

The majority of methylome-wide association studies (MWAS) have been performed using commercially available array-based technologies such as the Infinium Human Methylation 450K and the Infinium MethylationEPIC arrays (Illumina). While these arrays offer a convenient and relatively robust assessment of the probed sites they only allow interrogation of 2-4% of all CpG sites in the human genome. Methyl-binding domain sequencing (MBD-seq) is an alternative approach for MWAS that provides near-complete coverage of the methylome at similar costs as the array-based technologies. However, despite publication of multiple positive evaluations, the use of MBD-seq for MWAS is often fiercely criticized. Here we discuss key features of the method and debunk misconceptions using empirical data. We conclude that MBD-seq represents an excellent approach for large-scale MWAS and that increased utilization is likely to result in more discoveries, advance biological knowledge, and expedite the clinical translation of methylome-wide research findings.

Cell-Free DNA Methylation Profiling Analysis-Technologies and Bioinformatics

Analysis of circulating nucleic acids in bodily fluids, referred to as “liquid biopsies”, is rapidly gaining prominence. Studies have shown that cell-free DNA (cfDNA) has great potential in characterizing tumor status and heterogeneity, as well as the response to therapy and tumor recurrence. DNA methylation is an epigenetic modification that plays an important role in a broad range of biological processes and diseases. It is well known that aberrant DNA methylation is generalizable across various samples and occurs early during the pathogenesis of cancer. Methylation patterns of cfDNA are also consistent with their originated cells or tissues. Systemic analysis of cfDNA methylation profiles has emerged as a promising approach for cancer detection and origin determination. In this review, we will summarize the technologies for DNA methylation analysis and discuss their feasibility for liquid biopsy applications. We will also provide a brief overview of the bioinformatic approaches for analysis of DNA methylation sequencing data. Overall, this review provides informative guidance for the selection of experimental and computational methods in cfDNA methylation-based studies.

RaMWAS: fast methylome-wide association study pipeline for enrichment platforms

Motivation

Enrichment-based technologies can provide measurements of DNA methylation at tens of millions of CpGs for thousands of samples. Existing tools for methylome-wide association studies cannot analyze datasets of this size and lack important features like principal component analysis, combined analysis with SNP data and outcome predictions that are based on all informative methylation sites.

Results

We present a Bioconductor R package called RaMWAS with a full set of tools for large-scale methylome-wide association studies. It is free, cross-platform, open source, memory efficient and fast.

Availability and implementation

Release version and vignettes with small case study at bioconductor.org/packages/ramwas Development version at github.com/andreyshabalin/ramwas.

Supplementary information

Supplementary data are available at Bioinformatics online.

Application of the geographic population structure (GPS) algorithm for biogeographical analyses of wild and captive gorillas

Background

The utilization of high resolution genome data has important implications for the phylogeographical evaluation of non-human species. Biogeographical analyses can yield detailed understanding of their population biology and facilitate the geo-localization of individuals to promote their efficacious management, particularly when bred in captivity. The Geographic Population Structure (GPS) algorithm is an admixture based tool for inference of biogeographical affinities and has been employed for the geo-localization of various human populations worldwide. Here, we applied the GPS tool for biogeographical analyses and localization of the ancestral origins of wild and captive gorilla genomes, of unknown geographic source, available in the Great Ape Genome Project (GAGP), employing Gorillas with known ancestral origin as the reference data.

Results

Our findings suggest that GPS was successful in recapitulating the population history and estimating the geographic origins of all gorilla genomes queried and localized the wild gorillas with unknown geographical origin < 150 km of National Parks/Wildlife Reserves within the political boundaries of countries, considered as prominent modern-day abode for gorillas in the wild. Further, the GPS localization of most captive-born gorillas was congruent with their previously presumed ancestral homes.

Conclusions

Currently there is limited knowledge of the ancestral origins of most North American captive gorillas, and our study highlights the usefulness of GPS for inferring ancestry of captive gorillas. Determination of the native geographical source of captive gorillas can provide valuable information to guide breeding programs and ensure their appropriate management at the population level. Finally, our findings shine light on the broader applicability of GPS for protecting the genetic integrity of other endangered non-human species, where controlled breeding is a vital component of their conservation.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2568-5) contains supplementary material, which is available to authorized users.

DNA Methylation Analysis

DNA methylation is a process by which methyl groups are added to cytosine or adenine. DNA methylation can change the activity of the DNA molecule without changing the sequence. Methylation of 5-methylcytosine (5mC) is widespread in both eukaryotes and prokaryotes, and it is a very important epigenetic modification event, which can regulate gene activity and influence a number of key processes such as genomic imprinting, cell differentiation, transcriptional regulation, and chromatin remodeling. Profiling DNA methylation across the genome is critical to understanding the influence of methylation in normal biology and diseases including cancer. Recent discoveries of 5-methylcytosine (5mC) oxidation derivatives including 5-hydroxymethylcytosine (5hmC), 5-formylcytsine (5fC), and 5-carboxycytosine (5caC) in mammalian genome further expand our understanding of the methylation regulation. Genome-wide analyses such as microarrays and next-generation sequencing technologies have been used to assess large fractions of the methylome. A number of different quantitative approaches have also been established to map the DNA epigenomes with single-base resolution, as represented by the bisulfite-based methods, such as classical bisulfite sequencing, pyrosequencing etc. These methods have been used to generate base-resolution maps of 5mC and its oxidation derivatives in genomic samples. The focus of this chapter is to provide the methodologies that have been developed to detect the cytosine derivatives in the genomic DNA.

How to Design a Whole-Genome Bisulfite Sequencing Experiment

Aside from post-translational histone modifications and small RNA populations, the epigenome of an organism is defined by the level and spectrum of DNA methylation. Methyl groups can be covalently bound to the carbon-5 of cytosines or the carbon-6 of adenine bases. DNA methylation can be found in both prokaryotes and eukaryotes. In the latter, dynamic variation is shown across species, along development, and by cell type. DNA methylation usually leads to a lower binding affinity of DNA-interacting proteins and often results in a lower expression rate of the subsequent genome region, a process also referred to as transcriptional gene silencing. We give an overview of the current state of research facilitating the planning and implementation of whole-genome bisulfite-sequencing (WGBS) experiments. We refrain from discussing alternative methods for DNA methylation analysis, such as reduced representation bisulfite sequencing (rrBS) and methylated DNA immunoprecipitation sequencing (MeDIPSeq), which have value in specific experimental contexts but are generally disadvantageous compared to WGBS.

Convergence of evidence from a methylome-wide CpG-SNP association study and GWAS of major depressive disorder

DNA methylation is an epigenetic modification that provides stability and diversity to the cellular phenotype. It is influenced by both genetic sequence variation and environmental factors, and can therefore potentially account for variation of heritable phenotypes and disorders. Therefore, methylome-wide association studies (MWAS) are promising complements to genome-wide association studies (GWAS) of genetic variants. Of particular interest are methylation sites (CpGs) that are created or destroyed by the alleles of single-nucleotide polymorphisms (SNPs), as these so-called CpG-SNPs may show variation in methylation levels on top of what can be explained by the sequence variation. Using sequencing-based data from 1132 major depressive disorder (MDD) cases and controls, we performed a MWAS of 970,414 common CpG-SNPs. The analysis identified 27 suggestively significant (P < 1.00 × 10^-5) CpG-SNPs associations. Furthermore, the MWAS results were over-represented (odds ratios ranging 1.36-5.00; P ranging 4.9 × 10^-3-8.1 × 10^-2) among findings from three recent GWAS for MDD-related phenotypes. Overlapping loci included, e.g., ROBO2, ASIC2, and DCC. As the CpG-SNP analysis accounts for the number of alleles that creates CpGs, the methylation differences could not be explained by differences in allele frequencies. Thus, the results show that the MWAS and GWASs provide independent lines of evidence for the involvement of these loci in MDD. In conclusion, our methylation study of MDD contributes novel information about loci of relevance that complements previous findings and generates new hypothesis about MDD etiology, such as that the functional effects of genetic association may be partly mediated and/or enhanced by the methylation status in these loci.