A methylation study implicates the rewiring of brain neural circuits during puberty in the emergence of sex differences in depression symptoms.

BACKGROUND: Women are 1.5-3 times more likely to suffer from depression than men. This sex bias first emerges during puberty and then persists across the reproductive years. As the cause remains largely elusive, we performed a methylation-wide association study (MWAS) to generate novel hypotheses. METHODS: We assayed nearly all 28 million possible methylation sites in blood in 595 blood samples from 487 participants aged 9-17. MWASs were performed to identify methylation sites associated with increasing sex differences in depression symptoms as a function of pubertal stage. Epigenetic deconvolution was applied to perform analyses on a cell-type specific level. RESULTS: In monocytes, a substantial number of significant associations were detected after controlling the FDR at 0.05. These results could not be explained by plasma testosterone/estradiol or current/lifetime trauma. Our top results in monocytes were significantly enriched (ratio of 2.48) for genes in the top of a large genome-wide association study (GWAS) meta-analysis of depression and neurodevelopment-related Gene Ontology (GO) terms that remained significant after correcting for multiple testing. Focusing on our most robust findings (70 genes overlapping with the GWAS meta-analysis and the significant GO terms), we find genes coding for members of each of the major classes of axon guidance molecules (netrins, slits, semaphorins, ephrins, and cell adhesion molecules). Many of these genes were previously implicated in rodent studies of brain development and depression-like phenotypes, as well as human methylation, gene expression and GWAS studies. CONCLUSIONS: Our study suggests that the emergence of sex differences in depression may be related to the differential rewiring of brain circuits between boys and girls during puberty.

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.

Test-statistic inflation in methylome-wide association studies.

Recent years have seen a surge of methylome-wide association studies (MWAS). We observed that many of these studies suffer from test statistic inflation that is most likely caused by commonly used quality control (QC) pipelines not going far enough to remove technical artefacts. To support this claim, we reanalysed GEO datasets with an improved QC pipeline that reduced test-statistic inflation parameter lambda from the original mean/median of 20.16/15.17 to 3.07/1.14. Furthermore, the mean/median number of methylome-wide significant findings was reduced by 65,688/57,805 loci after more thorough QC. To avoid such false positives we argue for more extensive QC and that reporting the test-statistic inflation parameter lambda become standard for all MWAS allowing readers to better assess the risk of false discoveries.

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.

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.

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.

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.

A Whole Methylome Study of Ethanol Exposure in Brain and Blood: An Exploration of the Utility of Peripheral Blood as Proxy Tissue for Brain in Alcohol Methylation Studies.

BACKGROUND: Recent reviews have highlighted the potential use of blood-based methylation biomarkers as diagnostic and prognostic tools of current and future alcohol use and addiction. Due to the substantial overlap that often exists between methylation patterns across different tissues, including blood and brain, blood-based methylation may track methylation changes in brain; however, little work has explored the overlap in alcohol-related methylation in these tissues. METHODS: To study the effects of alcohol on the brain methylome and identify possible biomarkers of these changes in blood, we performed a methylome-wide association study in brain and blood from 40 male DBA/2J mice that received either an acute ethanol (EtOH) or saline intraperitoneal injection. To investigate all 22 million CpGs in the mouse genome, we enriched for the methylated genomic fraction using methyl-CpG binding domain (MBD) protein capture followed by next-generation sequencing (MBD-seq). We performed association tests in blood and brain separately followed by enrichment testing to determine whether there was overlapping alcohol-related methylation in the 2 tissues. RESULTS: The top result for brain was a CpG located in an intron of Ttc39b (p = 5.65 × 10-08 ), and for blood, the top result was located in Espnl (p = 5.11 × 10-08 ). Analyses implicated pathways involved in inflammation and neuronal differentiation, such as CXCR4, IL-7, and Wnt signaling. Enrichment tests indicated significant overlap among the top results in brain and blood. Pathway analyses of the overlapping genes converge on MAPKinase signaling (p = 5.6 × 10-05 ) which plays a central role in acute and chronic responses to alcohol and glutamate receptor pathways, which can regulate neuroplastic changes underlying addictive behavior. CONCLUSIONS: Overall, we have shown some methylation changes in brain and blood after acute EtOH administration and that the changes in blood partly mirror the changes in brain suggesting the potential for DNA methylation in blood to be biomarkers of alcohol use.

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.

Epigenetic Aging in Major Depressive Disorder.

OBJECTIVE: Major depressive disorder is associated with an increased risk of mortality and aging-related diseases. The authors examined whether major depression is associated with higher epigenetic aging in blood as measured by DNA methylation (DNAm) patterns, whether clinical characteristics of major depression have a further impact on these patterns, and whether the findings replicate in brain tissue. METHOD: DNAm age was estimated using all methylation sites in blood of 811 depressed patients and 319 control subjects with no lifetime psychiatric disorders and low depressive symptoms from the Netherlands Study of Depression and Anxiety. The residuals of the DNAm age estimates regressed on chronological age were calculated to indicate epigenetic aging. Major depression diagnosis and clinical characteristics were assessed with questionnaires and psychiatric interviews. Analyses were adjusted for sociodemographic characteristics, lifestyle, and health status. Postmortem brain samples of 74 depressed patients and 64 control subjects were used for replication. Pathway enrichment analysis was conducted using ConsensusPathDB to gain insight into the biological processes underlying epigenetic aging in blood and brain. RESULTS: Significantly higher epigenetic aging was observed in patients with major depression compared with control subjects (Cohen's d=0.18), with a significant dose effect with increasing symptom severity in the overall sample. In the depression group, epigenetic aging was positively and significantly associated with childhood trauma score. The case-control difference was replicated in an independent data set of postmortem brain samples. The top significantly enriched Gene Ontology terms included neuronal processes. CONCLUSIONS: As compared with control subjects, patients with major depression exhibited higher epigenetic aging in blood and brain tissue, suggesting that they are biologically older than their corresponding chronological age. This effect was even more profound in the presence of childhood trauma.

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.

Methyl-CpG-Binding Domain Sequencing: MBD-seq.

Detailed biological knowledge about the potential importance of the methylome is typically lacking for common diseases. Therefore, methylome-wide association studies (MWAS) are critical to detect disease relevant methylation sites. Methyl-CpG-binding domain sequencing (MBD-seq) offers potential advantages compared to antibody-based enrichment, but performance depends critically on using an optimal protocol. Using an optimized protocol, MBD-seq can approximate the sensitivity/specificity obtained with whole-genome bisulfite sequencing, but at a fraction of the costs and time to complete the project. Thus, MBD-seq offers a comprehensive first pass at the CpG methylome and is economically feasible with the samples sizes required for MWAS.

Association of Childhood Trauma Exposure With Adult Psychiatric Disorders and Functional Outcomes.

Importance: Being exposed to trauma is a common childhood experience associated with symptoms and impairments in childhood. Objective: To assess the association between cumulative childhood trauma exposure and adult psychiatric and functional outcomes. Design, Setting, and Participants: Prospective, population-based cohort study of 1420 participants. A community representative sample of participants was assessed with structured Child and Adolescent Psychiatric Assessment interviews up to 8 times in childhood (ages 9-16 years; 6674 observations; 1993-2000) for lifetime trauma exposure as defined by the Diagnostic and Statistical Manual of Mental Disorders. Participants were followed up 4 times in adulthood (ages 19, 21, 25, and 30 years; 4556 observations of 1336 participants; 1999-2015) with the structured Young Adult Psychiatric Assessment Interview for psychiatric outcomes, functional outcomes, and evidence of a disrupted transition to adulthood. Analysis was completed in 2018. Exposure: Participants were assessed with the structured Child and Adolescent Psychiatric Assessment interview (parent and self-report) up to 8 times in childhood for lifetime trauma exposure (ages 9-16 years; 6674 observations; 1993-2000). Main Outcomes and Measures: Participants were assessed up to 4 times with the structured Young Adult Psychiatric Assessment interview (self-report) in adulthood (ages 19, 21, 25, and 30 years; 4556 observations of 1336 participants; 1999-2015) for psychiatric outcomes, functional outcomes, and evidence of a disrupted transition to adulthood. Results: Among the 1420 study participants, 630 (49.0%) were female and 983 (89.4%) were white. By age 16 years, 30.9% of children (n = 451) were exposed to 1 traumatic event, 22.5% (n = 289) were exposed to 2 such events, and 14.8% (n = 267) were exposed to 3 or more. Cumulative childhood trauma exposure to age 16 years was associated with higher rates of adult psychiatric disorders (odds ratio for any disorder, 1.2; 95% CI, 1.0-1.4) and poorer functional outcomes, including key outcomes that indicate a significantly disrupted transition to adulthood (eg, failure to hold a job and social isolation). Childhood trauma exposure continued to be associated with higher rates of adult psychiatric and functional outcomes after adjusting for a broad range of childhood risk factors, including psychiatric functioning and family adversities and hardships (adjusted odds ratio for any disorder, 1.3; 95% CI, 1.0-1.5). Conclusions and Relevance: Cumulative childhood trauma exposure was associated with poor adult outcomes even after accounting for many of the childhood and family factors associated with both trauma exposure and poor adult outcomes. Childhood trauma exposures are common, but often preventable, thus providing a clear target for child-focused public health efforts to ameliorate long-term morbidity.

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

We recently showed that, after optimization, our methyl-CpG binding domain sequencing (MBD-seq) application approximates the methylome-wide coverage obtained with whole-genome bisulfite sequencing (WGB-seq), but at a cost that enables adequately powered large-scale association studies. A prior drawback of MBD-seq is the relatively large amount of genomic DNA (ideally >1 µg) required to obtain high-quality data. Biomaterials are typically expensive to collect, provide a finite amount of DNA, and may simply not yield sufficient starting material. The ability to use low amounts of DNA will increase the breadth and number of studies that can be conducted. Therefore, we further optimized the enrichment step. With this low starting material protocol, MBD-seq performed equally well, or better, than the protocol requiring ample starting material (>1 µg). Using only 15 ng of DNA as input, there is minimal loss in data quality, achieving 93% of the coverage of WGB-seq (with standard amounts of input DNA) at similar false/positive rates. Furthermore, across a large number of genomic features, the MBD-seq methylation profiles closely tracked those observed for WGB-seq with even slightly larger effect sizes. This suggests that MBD-seq provides similar information about the methylome and classifies methylation status somewhat more accurately. Performance decreases with <15 ng DNA as starting material but, even with as little as 5 ng, MBD-seq still achieves 90% of the coverage of WGB-seq with comparable genome-wide methylation profiles. Thus, the proposed protocol is an attractive option for adequately powered and cost-effective methylome-wide investigations using (very) low amounts of DNA.

Enrichment methods provide a feasible approach to comprehensive and adequately powered investigations of the brain methylome.

Methylome-wide association studies are typically performed using microarray technologies that only assay a very small fraction of the CG methylome and entirely miss two forms of methylation that are common in brain and likely of particular relevance for neuroscience and psychiatric disorders. The alternative is to use whole genome bisulfite (WGB) sequencing but this approach is not yet practically feasible with sample sizes required for adequate statistical power. We argue for revisiting methylation enrichment methods that, provided optimal protocols are used, enable comprehensive, adequately powered and cost-effective genome-wide investigations of the brain methylome. To support our claim we use data showing that enrichment methods approximate the sensitivity obtained with WGB methods and with slightly better specificity. However, this performance is achieved at <5% of the reagent costs. Furthermore, because many more samples can be sequenced simultaneously, projects can be completed about 15 times faster. Currently the only viable option available for comprehensive brain methylome studies, enrichment methods may be critical for moving the field forward.

Genomewide Association Study of Alcohol Dependence Identifies Risk Loci Altering Ethanol-Response Behaviors in Model Organisms.

BACKGROUND: Alcohol dependence (AD) shows evidence for genetic liability, but genes influencing risk remain largely unidentified. METHODS: We conducted a genomewide association study in 706 related AD cases and 1,748 unscreened population controls from Ireland. We sought replication in 15,496 samples of European descent. We used model organisms (MOs) to assess the role of orthologous genes in ethanol (EtOH)-response behaviors. We tested 1 primate-specific gene for expression differences in case/control postmortem brain tissue. RESULTS: We detected significant association in COL6A3 and suggestive association in 2 previously implicated loci, KLF12 and RYR3. None of these signals are significant in replication. A suggestive signal in the long noncoding RNA LOC339975 is significant in case:control meta-analysis, but not in a population sample. Knockdown of a COL6A3 ortholog in Caenorhabditis elegans reduced EtOH sensitivity. Col6a3 expression correlated with handling-induced convulsions in mice. Loss of function of the KLF12 ortholog in C. elegans impaired development of acute functional tolerance (AFT). Klf12 expression correlated with locomotor activation following EtOH injection in mice. Loss of function of the RYR3 ortholog reduced EtOH sensitivity in C. elegans and rapid tolerance in Drosophila. The ryanodine receptor antagonist dantrolene reduced motivation to self-administer EtOH in rats. Expression of LOC339975 does not differ between cases and controls but is reduced in carriers of the associated rs11726136 allele in nucleus accumbens (NAc). CONCLUSIONS: We detect association between AD and COL6A3, KLF12, RYR3, and LOC339975. Despite nonreplication of COL6A3, KLF12, and RYR3 signals, orthologs of these genes influence behavioral response to EtOH in MOs, suggesting potential involvement in human EtOH response and AD liability. The associated LOC339975 allele may influence gene expression in human NAc. Although the functions of long noncoding RNAs are poorly understood, there is mounting evidence implicating these genes in multiple brain functions and disorders.

Correcting for cell-type effects in DNA methylation studies: reference-based method outperforms latent variable approaches in empirical studies.

Based on an extensive simulation study, McGregor and colleagues recently recommended the use of surrogate variable analysis (SVA) to control for the confounding effects of cell-type heterogeneity in DNA methylation association studies in scenarios where no cell-type proportions are available. As their recommendation was mainly based on simulated data, we sought to replicate findings in two large-scale empirical studies. In our empirical data, SVA did not fully correct for cell-type effects, its performance was somewhat unstable, and it carried a risk of missing true signals caused by removing variation that might be linked to actual disease processes. By contrast, a reference-based correction method performed well and did not show these limitations. A disadvantage of this approach is that if reference methylomes are not (publicly) available, they will need to be generated once for a small set of samples. However, given the notable risk we observed for cell-type confounding, we argue that, to avoid introducing false-positive findings into the literature, it could be well worth making this investment.Please see related Correspondence article: https://genomebiology.biomedcentral.com/articles/10/1186/s13059-017-1149-7 and related Research article: https://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-0935-y.

Evaluation of Methyl-Binding Domain Based Enrichment Approaches Revisited.

Methyl-binding domain (MBD) enrichment followed by deep sequencing (MBD-seq), is a robust and cost efficient approach for methylome-wide association studies (MWAS). MBD-seq has been demonstrated to be capable of identifying differentially methylated regions, detecting previously reported robust associations and producing findings that replicate with other technologies such as targeted pyrosequencing of bisulfite converted DNA. There are several kits commercially available that can be used for MBD enrichment. Our previous work has involved MethylMiner (Life Technologies, Foster City, CA, USA) that we chose after careful investigation of its properties. However, in a recent evaluation of five commercially available MBD-enrichment kits the performance of the MethylMiner was deemed poor. Given our positive experience with MethylMiner, we were surprised by this report. In an attempt to reproduce these findings we here have performed a direct comparison of MethylMiner with MethylCap (Diagenode Inc, Denville, NJ, USA), the best performing kit in that study. We find that both MethylMiner and MethylCap are two well performing MBD-enrichment kits. However, MethylMiner shows somewhat better enrichment efficiency and lower levels of background "noise". In addition, for the purpose of MWAS where we want to investigate the majority of CpGs, we find MethylMiner to be superior as it allows tailoring the enrichment to the regions where most CpGs are located. Using targeted bisulfite sequencing we confirmed that sites where methylation was detected by either MethylMiner or by MethylCap indeed were methylated.