Methylome-wide Association Studies Research Articles

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.

Implication of DNA methylation during lifestyle mediated weight loss.

Over the past 50 years, the number of overweight/obese people increased significantly, making obesity a global public health challenge. Apart from rare monogenic forms, obesity is a multifactorial disease, most likely resulting from a concerted interaction of genetic, epigenetic and environmental factors. Although recent studies opened new avenues in elucidating the complex genetics behind obesity, the biological mechanisms contributing to individual's risk to become obese are not yet fully understood. Non-genetic factors such as eating behaviour or physical activity are strong contributing factors for the onset of obesity. These factors may interact with genetic predispositions most likely via epigenetic mechanisms. Epigenome-wide association studies or methylome-wide association studies are measuring DNA methylation at single CpGs across thousands of genes and capture associations to obesity phenotypes such as BMI. However, they only represent a snapshot in the complex biological network and cannot distinguish between causes and consequences. Intervention studies are therefore a suitable method to control for confounding factors and to avoid possible sources of bias. In particular, intervention studies documenting changes in obesity-associated epigenetic markers during lifestyle driven weight loss, make an important contribution to a better understanding of epigenetic reprogramming in obesity. To investigate the impact of lifestyle in obesity state specific DNA methylation, especially concerning the development of new strategies for prevention and individual therapy, we reviewed 19 most recent human intervention studies. In summary, this review highlights the huge potential of targeted interventions to alter disease-associated epigenetic patterns. However, there is an urgent need for further robust and larger studies to identify the specific DNA methylation biomarkers which influence obesity.

Relationship between 233 colorectal cancer risk loci and survival in 1926 patients with advanced disease

BackgroundGenome, transcriptome and methylome-wide association studies have identified single-nucleotide polymorphisms (SNPs) or genes at 258 loci associated with colorectal cancer (CRC) risk. We studied the relationship between these and patient outcome.MethodsWe studied 1926 unrelated patients with advanced CRC from COIN and COIN-B. Of 205 CRC-risk SNPs, 19 were directly genotyped and 162 were imputed, and of 53 risk genes, 52 were tested. An additive Cox model for overall survival was adjusted for known prognostic factors. For nominally significant SNPs or genes, we considered a recessive model with a Bonferroni corrected threshold of P = 2.1 × 10−4. We examined SNPs as expression quantitative trait loci (eQTL) and the relationship between gene expression in colorectal tumours and survival in 597 unrelated patients.ResultsEleven SNPs or genes were nominally associated with survival under an additive model. Only rs117079142 mapping to UTP23 and EIF3H (Hazard Ratio [HR] = 2.79, 95% Confidence Intervals [CI] = 1.70–4.58, P = 4.7 × 10−5) and rs9924886 mapping to CDH1 and CDH3 (HR = 1.24, 95% CI = 1.12–1.38, P = 5.2 × 10−5) passed the multiple testing threshold under a recessive model. rs117079142 was an eQTL for UTP23 and rs9924886 for CDH1, CDH3 and ZFP90. Decreased CDH1 expression in CRCs was associated with worse survival (HR = 2.18, 95% CI = 1.3–3.5, P = 1.8 × 10−3).Conclusionrs117079142 and rs9924886 may represent potential prognostic biomarkers for CRC.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

The Cannabis-Induced Epigenetic Regulation of Genes Associated with Major Depressive Disorder.

The prevalence of depression is increasing worldwide, as is the number of people suffering from treatment-resistant depression; these patients constitute 30% of those treated. Unfortunately, there have not been significant advances in the treatment of this disorder in the past few decades. Exposure to cannabis and cannabis-derived compounds impacts depression symptomatology in different ways, with evidence indicating that cannabidiol has antidepressant effects; there have been mixed results with medical cannabis. Even though the exact molecular mechanisms of the action underlying changes in depression symptomatology upon exposure to cannabis and cannabis-derived compounds are still unknown, there is strong evidence that these agents have a widespread impact on epigenetic regulation. We hypothesized that exposure to cannabis or cannabis-derived compounds changes the DNA methylation levels of genes associated with depression. To test this hypothesis, we first performed a literature search to identify genes that are differentially methylated upon exposure to cannabis and cannabis-derived compounds, as reported in methylome-wide association studies. We next checked whether genes residing in loci associated with depression, as identified in the largest currently available genome-wide association study of depression, were reported to be epigenetically regulated by cannabis or cannabis-related compounds. Multiple genes residing in loci associated with depression were found to be epigenetically regulated by exposure to cannabis or cannabis-derived compounds. This epigenomic regulation of depression-associated genes by cannabis or cannabis-derived compounds was reported across diverse organisms, tissues, and developmental stages and occurred in genes crucial for neuronal development, functioning, survival, and synapse functioning, as well as in genes previously implicated in other mental disorders.

Integrated methylome and phenome study of the circulating proteome reveals markers pertinent to brain health

Characterising associations between the methylome, proteome and phenome may provide insight into biological pathways governing brain health. Here, we report an integrated DNA methylation and phenotypic study of the circulating proteome in relation to brain health. Methylome-wide association studies of 4058 plasma proteins are performed (N = 774), identifying 2928 CpG-protein associations after adjustment for multiple testing. These are independent of known genetic protein quantitative trait loci (pQTLs) and common lifestyle effects. Phenome-wide association studies of each protein are then performed in relation to 15 neurological traits (N = 1,065), identifying 405 associations between the levels of 191 proteins and cognitive scores, brain imaging measures or APOE e4 status. We uncover 35 previously unreported DNA methylation signatures for 17 protein markers of brain health. The epigenetic and proteomic markers we identify are pertinent to understanding and stratifying brain health.

Complex trait methylation scores in the prediction of major depressive disorder.

SummaryBackgroundDNA methylation (DNAm) is associated with time-varying environmental factors that contribute to major depressive disorder (MDD) risk. We sought to test whether DNAm signatures of lifestyle and biochemical factors were associated with MDD to reveal dynamic biomarkers of MDD risk that may be amenable to lifestyle interventions.MethodsHere, we calculated methylation scores (MS) at multiple p-value thresholds for lifestyle (BMI, smoking, alcohol consumption, and educational attainment) and biochemical (high-density lipoprotein (HDL) and total cholesterol) factors in Generation Scotland (GS) (N=9,502) and in a replication cohort (ALSPACadults, N=565), using CpG sites reported in previous well-powered methylome-wide association studies. We also compared their predictive accuracy for MDD to a MDD MS in an independent GS sub-sample (N=4,432).FindingsEach trait MS was significantly associated with its corresponding phenotype in GS (βrange=0.089–1.457) and in ALSPAC (βrange=0.078–2.533). Each MS was also significantly associated with MDD before and after adjustment for its corresponding phenotype in GS (βrange=0.053–0.145). After accounting for relevant lifestyle factors, MS for educational attainment (β=0.094) and alcohol consumption (MSp-value<0.01–0.5; βrange=-0.069–0.083) remained significantly associated with MDD in GS. Smoking (AUC=0.569) and educational attainment (AUC=0.585) MSs could discriminate MDD from controls better than the MDD MS (AUC=0.553) in the independent GS sub-sample. Analyses implicating MDD did not replicate across ALSPAC, although the direction of effect was consistent for all traits when adjusting for the MS corresponding phenotypes.InterpretationWe showed that lifestyle and biochemical MS were associated with MDD before and after adjustment for their corresponding phenotypes (pnominal<0.05), but not when smoking, alcohol consumption, and BMI were also included as covariates. MDD results did not replicate in the smaller, female-only independent ALSPAC cohort (NALSPAC=565; NGS=9,502), potentially due to demographic differences or low statistical power, but effect sizes were consistent with the direction reported in GS. DNAm scores for modifiable MDD risk factors may contribute to disease vulnerability and, in some cases, explain additional variance to their observed phenotypes.FundingWellcome Trust.

P60. Meta-Analysis of Cell Type-Specific DNA Methylation of Childhood Attention-Deficit/Hyperactivity Disorder Symptoms

The interplay of genetic factors and environmental adversities in ADHD can be reflected by cell type-specific DNA methylation (DNAm). Previous methylome-wide association studies (MWAS) for ADHD have been performed in bulk tissue only, where meaningful DNAm signal could be left undetected. We aim to identify cell type-specific DNAm profiles associated with childhood ADHD symptoms.

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 changes induced by prenatal toxic metal exposure: An overview of epidemiological evidence.

Accumulating evidence suggests that exposure to unfavorable conditions early in life can substantially contribute to the risk of chronic disorders later in life (‘developmental programming’ phenomenon). The mechanistic basis for this phenomenon remains poorly understood so far, although epigenetic mechanisms such as DNA methylation, histone modifications and microRNA-mediated gene regulation apparently play a crucial role. The key role of epigenetic modifications triggered by unfavorable environmental cues during sensitive developmental periods in linking adverse early-life events to later-life health outcomes is evident from a large body of studies, including methylome-wide association studies and research of candidate genes. Toxic metals (TMs), such as heavy metals, including lead, chromium, cadmium, arsenic, mercury, etc., are among environmental contaminants currently most significantly impacting human health status. Since TMs can cross the placental barrier and accumulate in fetal tissues, exposure to high doses of these xenobiotics early in development is considered to be among important factors contributing to the developmental programming of adult-life diseases in modern societies. In this mini-review, we summarize epidemiological findings indicating that prenatal TM exposure can induce epigenetic dysregulation, thereby potentially affecting adult health outcomes.

Methylome-wide association study of early life stressors and adult mental health.

The environment and events that we are exposed to in utero, during birth and in early childhood influence our future physical and mental health. The underlying mechanisms that lead to these outcomes are unclear, but long-term changes in epigenetic marks, such as DNA methylation, could act as a mediating factor or biomarker. DNA methylation data were assayed at 713 522 CpG sites from 9537 participants of the Generation Scotland: Scottish Family Health Study, a family-based cohort with extensive genetic, medical, family history and lifestyle information. Methylome-wide association studies of eight early life environment phenotypes and two adult mental health phenotypes (major depressive disorder and brief resilience scale) were conducted using DNA methylation data collected from adult whole blood samples. Two genes involved with different developmental pathways (PRICKLE2, Prickle Planar Cell Polarity Protein 2 and ABI1, Abl-Interactor-1) were annotated to CpG sites associated with preterm birth (P < 1.27 × 10−9). A further two genes important to the development of sensory pathways (SOBP, Sine Oculis Binding Protein Homolog and RPGRIP1, Retinitis Pigmentosa GTPase Regulator Interacting Protein) were annotated to sites associated with low birth weight (P < 4.35 × 10−8). The examination of methylation profile scores and genes and gene-sets annotated from associated CpGs sites found no evidence of overlap between the early life environment and mental health conditions. Birth date was associated with a significant difference in estimated lymphocyte and neutrophil counts. Previous studies have shown that early life environments influence the risk of developing mental health disorders later in life; however, this study found no evidence that this is mediated by stable changes to the methylome detectable in peripheral blood.

Methylomic Investigation of Problematic Adolescent Cannabis Use and Its Negative Mental Health Consequences

The impact of adolescent cannabis use is a pressing public health question owing to the high rates of use and links to negative outcomes. This study considered the association between problematic adolescent cannabis use and methylation. Using an enrichment-based sequencing approach, a methylome-wide association study (MWAS) was performed of problematic adolescent cannabis use in 703 adolescent samples from the Great Smoky Mountain Study. Using epigenomic deconvolution, MWASs were performed for the main cell types in blood: granulocytes, T cells, B cells, and monocytes. Enrichment testing was conducted to establish overlap between cannabis-associated methylation differences and variants associated with negative mental health effects of adolescent cannabis use. Whole-blood analyses identified 45 significant CpGs, and cell type-specific analyses yielded 32 additional CpGs not identified in the whole-blood MWAS. Significant overlap was observed between the B-cell MWAS and genetic studies of education attainment and intelligence. Furthermore, the results from both T cells and monocytes overlapped with findings from an MWAS of psychosis conducted in brain tissue. In one of the first methylome-wide association studies of adolescent cannabis use, several methylation sites located in genes of importance for potentially relevant brain functions were identified. These findings resulted in several testable hypotheses by which cannabis-associated methylation can impact neurological development and inflammation response as well as potential mechanisms linking cannabis-associated methylation to potential downstream mental health effects.

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

ABSTRACTThe 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 Type–Specific Methylome-wide Association Studies Implicate Neurotrophin and Innate Immune Signaling in Major Depressive Disorder

BackgroundWe 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. MethodsIn 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. ResultsCell 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. ConclusionsWe 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.

Differential DNA methylation in blood as a mediator of the association between cigarette smoking and bladder cancer risk among postmenopausal women

ABSTRACTSmoking accounts for approximately 52% of bladder cancer incidence among postmenopausal women, but the underlying mechanism is poorly understood. Our study investigates whether changes in DNA methylation, as measured in blood, mediate the impact of smoking on bladder cancer risk among postmenopausal women. We conducted analyses among 206 cases and 251 controls that were current or never smokers at baseline from a previous case-control study of bladder cancer and genome-wide DNA methylation nested within the Women’s Health Initiative. Separate mediation analyses were conducted for three CpG sites demonstrating robust associations with smoking in prior methylome-wide association studies: cg05575921 (AhRR), cg03636183 (F2RL3), and cg19859270 (GPR15). We estimated causal effects using the regression-based, four-way decomposition approach, which addresses the interaction between smoking and each CpG site. The overall proportion of the excess relative risk mediated by cg05575921 was 92% (p-value = 0.004) and by cg19859270 was 79% (p-value = 0.02). The largest component of the excess relative risk of bladder cancer due to 30 pack-years of smoking history in current smokers was the mediated interaction for both cg05575921 (72%, p = 0.02) and cg19859270 (72%, p-value = 0.04), where the mediated interaction is the effect of smoking on bladder cancer that both acts through differential methylation and depends on smoking history. There was little evidence that smoking was mediated through cg03636183. Our results suggest that differential methylation of cg05575921 and cg19859270 mediate the effects of smoking on bladder cancer, potentially revealing downstream effects of smoking relevant for carcinogenesis.

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.

47 - METHYLOME-WIDE ASSOCIATION STUDIES FOR MAJOR DEPRESSIVE DISORDER IN BLOOD OVERLAP WITH METHYLATION RESULTS FROM BRAIN AND LARGE-SCALE GWAS

Background Epigenetic modifications such as DNA methylation provide stability and diversity to the cellular phenotype and aberrant methylation has been implicated in processes underlying psychiatric disorders. Therefore, studies combining DNA methylation and genotype information provide a promising approach to study disorders where genotype information alone has failed to reveal the full etiology. Methods We applied an optimized MBD-seq protocol to assay the complete CpG methylome in cases with Major Depressive Disorder (MDD) and controls using blood samples (N=1,132) from Netherlands Study of Depression and Anxiety and brain samples (N=64) from the Victorian Brain Bank Network. Data were analyzed with RaMWAS, a novel Bioconductor package specifically designed for Methylome-Wide Association Studies (MWAS). To study the overlap between top MWAS findings in blood and brain, we used a permutation based enrichment test (shiftR) that accounted for the dependency between adjacent CpG sites. Furthermore, we utilized the methylation data in combination with existing genotype information from the same individuals in a MWAS of CpGs created or destroyed by SNPs. Next, we tested whether top results from this CpG-SNP MWAS overlapped with recent large-scale GWAS to identify robust associations with genomic loci of importance for MDD etiology. Results The MWAS in blood identified five methylome-wide significant sites (P Discussion This work involves the largest MWAS for MDD performed to date. Our integrated analysis with brain tissue, genotype information and GWAS results highlighted biological functions of potential value for MDD etiology. Part of the associated methylation marks in blood overlapped with MWAS finding in brain. As blood can easily be collected in a clinical setting, these loci may be of direct value as potential diagnostic biomarkers for MDD.

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=1,132) and postmortem brain tissues (N=61 samples from Brodmann Area 10, BA10). The MWAS for blood identified several loci with P 1.91×10−8-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 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.

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.

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

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. 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. Release version and vignettes with small case study at bioconductor.org/packages/ramwas Development version at github.com/andreyshabalin/ramwas. Supplementary data are available at Bioinformatics online.