Robin F. Chan

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.

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

Contrasting influences of Drosophila white/mini-white on ethanol sensitivity in two different behavioral assays

BACKGROUND The fruit fly Drosophila melanogaster has been used extensively to investigate genetic mechanisms of ethanol (EtOH)-related behaviors. Many past studies in flies, including studies from our laboratory, have manipulated gene expression using transposons carrying the genetic-phenotypic marker mini-white(mini-w), a derivative of the endogenous gene white(w). Whether the mini-w transgenic marker or the endogenous w gene influences behavioral responses to acute EtOH exposure in flies has not been systematically investigated. METHODS We manipulated mini-w and w expression via (i) transposons marked with mini-w, (ii) RNAi against mini-w and w, and (iii) a null allele of w. We assessed EtOH sensitivity and tolerance using a previously described eRING assay (based on climbing in the presence of EtOH) and an assay based on EtOH-induced sedation. RESULTS In eRING assays, EtOH-induced impairment of climbing correlated inversely with expression of the mini-w marker from a series of transposon insertions. Additionally, flies harboring a null allele of w or flies with RNAi-mediated knockdown of mini-w were significantly more sensitive to EtOH in eRING assays than controls expressing endogenous w or the mini-w marker. In contrast, EtOH sensitivity and rapid tolerance measured in the EtOH sedation assay were not affected by decreased expression of mini-w or endogenous w in flies. CONCLUSIONS EtOH sensitivity measured in the eRING assay is noticeably influenced by w and mini-w, making eRING problematic for studies on EtOH-related behavior in Drosophila using transgenes marked with mini-w. In contrast, the EtOH sensitivity assay described here is a suitable behavioral paradigm for studies on EtOH sensitivity and rapid tolerance in Drosophila including those that use widely available transgenes marked with mini-w.

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

Abstract 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.

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

ABSTRACT 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.

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.

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.

An Inexpensive, Scalable Behavioral Assay for Measuring Ethanol Sedation Sensitivity and Rapid Tolerance in Drosophila

Alcohol use disorder (AUD) is a serious health challenge. Despite a large hereditary component to AUD, few genes have been unambiguously implicated in their etiology. The fruit fly, Drosophila melanogaster, is a powerful model for exploring molecular-genetic mechanisms underlying alcohol-related behaviors and therefore holds great promise for identifying and understanding the function of genes that influence AUD. The use of the Drosophila model for these types of studies depends on the availability of assays that reliably measure behavioral responses to ethanol. This report describes an assay suitable for assessing ethanol sensitivity and rapid tolerance in flies. Ethanol sensitivity measured in this assay is influenced by the volume and concentration of ethanol used, a variety of previously reported genetic manipulations, and also the length of time the flies are housed without food immediately prior to testing. In contrast, ethanol sensitivity measured in this assay is not affected by the vigor of fly handling, sex of the flies, and supplementation of growth medium with antibiotics or live yeast. Three different methods for quantitating ethanol sensitivity are described, all leading to essentially indistinguishable ethanol sensitivity results. The scalable nature of this assay, combined with its overall simplicity to set-up and relatively low expense, make it suitable for small and large scale genetic analysis of ethanol sensitivity and rapid tolerance in Drosophila.

Cell-type-specific methylome-wide association studies implicate neurodegenerative processes and neuroimmune communication in major depressive disorder

We studied the methylome in three collections of human postmortem brain (N=206) and blood samples (N=1,132) of subjects with major depressive disorder (MDD) and controls. Using an epigenomic deconvolution approach we performed cell-type-specific methylome-wide association studies (MWAS) within sub-populations of neurons/glia and granulocytes/T-cells/B-cells/monocytes for bulk brain and blood data, respectively. Multiple MWAS findings in neurons/glia replicated across brain collections (ORs=509-538, P-values<1×10−5) and were reproducible in an array-based MWAS of sorted neurons/glia from a fourth brain collection (N=58). Pathway analyses implicated p75NTR/VEGF signaling, neurodegeneration, and blood-brain barrier perturbation. Cell-type-specific analysis in blood identified associations in CD14+ monocytes -- a cell type strongly linked to neuroimmune processes and stress. Top results in neurons/glia/bulk and monocytes were enriched for genes supported by GWAS for MDD (ORs=2.02-2.87, P-values=0.003 to <1×10−5), neurodegeneration and other psychiatric disorders. In summary, we identified novel MDD-methylation associations by using epigenomic deconvolution that provided important mechanistic insights for the disease.

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.