- Research article
- Open Access
Insights into the origin of DNA methylation differences between monozygotic twins discordant for schizophrenia
Journal of Molecular Psychiatryvolume 3, Article number: 7 (2015)
DNA methylation differences between monozygotic twins discordant for schizophrenia have been previously reported. However, the origin of methylation differences between monozygotic twins discordant for schizophrenia is not clear. The findings here argue that all DNA methylation differences may not necessarily represent the cause of the disease; rather some may result from the effect of antipsychotics.
Methylation differences in rat brain regions and also in two pairs of unrelated monozygotic twins discordant for schizophrenia have been studied using genome-wide DNA methylation arrays at Arraystar Inc. (Rockville, Maryland, USA). The identified gene promoters showing significant alterations to DNA methylation were then further characterized using ingenuity pathway analysis (Ingenuity System Inc, CA, USA).
Pathway analysis of the most significant gene promoter hyper/hypomethylation revealed a significant enrichment of DNA methylation changes in biological networks and pathways directly relevant to neural development and psychiatric disorders. These included HIPPO signaling (p = 3.93E-03) and MAPK signaling (p = 4.27E-03) pathways involving hypermethylated genes in schizophrenia-affected patients as compared to their unaffected co-twins. Also, a number of significant pathways and networks involving genes with hypomethylated gene promoters have been identified. These included CREB signaling in neurons (p = 1.53E-02), Dopamine-DARPP32 feedback in cAMP signaling (p = 7.43E-03) and Ephrin receptors (p = 1.13E-02). Further, there was significant enrichment for pathways involved in nervous system development and function (p = 1.71E-03-4.28E-02).
The findings highlight the significance of antipsychotic drugs on DNA methylation in schizophrenia patients. The unique pathways affected by DNA methylation in the two pairs of monozygotic twins suggest that patient-specific pathways are responsible for the disease; suggesting that patient-specific treatment strategies may be necessary in treating the disorder. The study reflects the need for developing personalized medicine approaches that take into consideration epigenetic variations between patients.
Epigenetic modifications may affect gene expression without the existence of changes in DNA sequence . They are known to affect a number of neurological disorders [2, 3] including psychosis . Specifically, DNA methylation is known to affect gene expression via transcriptional regulation . Further, a number of environmental factors may affect gene-specific DNA methylation, and this alteration may be maintained during genetic transmission through mitosis and meiosis . These dynamic changes may explain the discordance of monozygotic twins  that are considered genetically identical. In fact, several studies have reported methylation differences between monozygotic twins discordant for a variety of disorders including schizophrenia [8, 9]. Interestingly, hypermethylation has been reported to affect genes that have been implicated in psychiatric disorders . These results have been used to argue that aberrations in gene specific DNA methylation must represent, at least, one of the contributing mechanisms in the manifestation of psychiatric disorders . The epigenetic float and drift during development can be stochastic, caused by environmental factors or (most likely) the result of a combination of the two . Although the involvement of stochastic processes have been difficult to establish, a number of environmental agents are known to affect epigenetic-dependent disease phenotypes in a number of disorders. An example of this includes the hypermethylation of tumor suppressor genes in a variety of cancers . Such results question the cause of differentially methylated psychosis-related genes in monozygotic twins discordant for psychosis. There is a possibility that the observed methylation difference is caused by (still unknown) environment factors, such as drugs. This suggests that some of the observed differential methylation between monozygotic twins discordant for this disorder may not necessarily be the cause of the disorder in the affected twin.
In this report, we discuss two complementary sets of results. The first deals with identification of hypermethylation/hypomethylation of a set of genes in the affected individuals of two sets of unrelated monozygotic twins discordant for schizophrenia. The results support the literature that monozygotic twins do differ in the DNA methylation status of their genes and the hypermethylated/hypomethylation genes may take part in pathways affecting psychiatric disorders. Our second set of results shows that an antipsychotic drug (Olanzapine) treatment may cause its effect via hypermethylation/hypomethylation of a number of genes that have been implicated in psychosis . Interestingly, several genes affected by hyper/hypomethylation in the second study are psychosis related genes. Taken together, the results offer novel insight into the epigenetic cause(s) and treatment of psychotic disorders. It forms the focus of this report.
This study involves two psychosis-affected monozygotic twins and their healthy co-twins (n = 4). All subjects gave written informed consent to participate in the study and contributed whole blood samples for DNA analysis. The demographic information and the medication status of the affected twins are described elsewhere . Genomic DNA was extracted from the blood samples using the PerfectPure DNA Blood Kit, following the manufacture’s protocol (http://www.5prime.com). Samples of genomic DNA were sonicated to random fragments in a size range of about 200–1000 bp. Immunoprecipitation of methylated DNA fragments was performed using BiomagTM magnetic beads coupled mouse monoclonal antibody against 5-methylcytidine. The precipitated DNA was eluted and purified by phenol chloroform extraction and ethanol precipitation. The total input and immunoprecipitated DNA were labeled with Cy3- and Cy5-labeled random 9-mers, respectively, and hybridized to the NimbleGen Human DNA Methylation 3x720K Promoter Plus CpG Island Array, which is a multiplex slide with 3 identical arrays per slide and each array contains 27,728 CpG Islands annotated by UCSC and 22,532 well-characterized RefSeq promoter regions (from about −2,440 to +610 bp of the TSSs) totally covered by ~720,000 probes. Scanning was performed with the Axon GenePix 4000B microarray scanner. This methylation protocol was carried out at ArrayStar Inc (Rockville, MD, USA).
The details of olanzapine treatment of experimental rats are described fully elsewhere . Briefly, adult male Sprague–Dawley rats of 12 weeks of age (250 – 300 g) were purchased from Charles River, PQ, Canada. Rats were weighed and divided into two treatment groups (olanzapine treated and control) with comparable means of weight. Rats received injections of olanzapine (Zyprexa, Lilly, IN, USA; 2.5 mg/kg, i.m.; n = 8) or vehicle (phosphate buffered saline (PBS); n = 8) between 1:30 and 3:00 pm daily for 19 days.
Methylation enrichment and peak-finding
Raw data generated from scans were normalized (log2-ratio). It followed Median centering, quantile normalization and linear smoothing using Bioconductor packages Ringo, limma, and MEDME. The results were used in a sliding-window (750 bp) peak-finding algorithm provided by NimbleScan v2.5 (Roche-NimbleGen). A one-sided Kolmogorov-Smirnov (KS) test is applied to determine whether the probes are drawn from a significantly more positive distribution of intensity log2-ratios than those in the rest of the array. Each probe receives a log10 p-value score from the windowed KS test around that probe. If several adjacent probes rise significantly above a set threshold, the region is assigned to an enrichment peak (EP). The peak data files are generated from the p-value data files. NimbleScan detects peaks by searching for at least two probes above a p-value minimum cutoff (−log10) of 2. Peaks within 500 bp of each other are merged. Further details of the quality control and normalization procedures in the twins and rat datasets are presented elsewhere [14, 15]. The genome-wide methylation analysis was applied to two unrelated pairs of monozygotic twins discordant for schizophrenia. For microarray data of replicate samples, a T-test was performed between different groups to calculate the p values for each probe. Probes with p < 0.05 were identified as differentially methylated probes, which were further analyzed to find Differentially Methylated Regions (DMR). Finally, we have applied a more stringent threshold of significant differences in DNA methylation (p < 0.001) to filter the genes affected by DNA methylation changes in their promoter regions. It allowed identification of biologically relevant genes showing hypermethylation in schizophrenia-affected twins as compared to their unaffected co-twin.
DEP analysis using M’ method
When comparing two groups’ differentially enriched regions, the average of the log2-ratio values for each group (Experiment and Control) was obtained and M’ value (defined by the following equation) was calculated for each probe. Then, the NimbleScan sliding-window peak-finding algorithm was rerun on the data to find the differential enrichment peaks (DEP).
The differential enrichment peaks (DEP) called by the NimbleScan algorithm were filtered according to the following criteria:
At least one of the two groups has a median (log2 MeDIP/Input) > =0.3 and M’ > 0.
At least half of probes in a peak may have coefficient of variability (CV) < = 0.8 in both two groups.
Genes affected by DMRs were analyzed with Ingenuity Pathway Analysis (IPA) for each patient. Also, the genes identified were matched with genes in the schizophrenia gene database (www.szgene.org). Top pathways and networks related to the two pairs are reported here. Finally, we have assessed whether the genes affected by promoter hyper/hypomethylation were among the psychosis related genes listed in Gemma (www.chibi.ubc.ca/Gemma/home.html). Furthermore, IPA was also used to identify the most significant pathways and networks affected by hyper/hypomethylation in the rat brain regions. Similarities and differences between pathways and networks affected by hyper/hypomethylation in the twins and in rat brain regions (hippocampus and cerebellum) exposed to a therapeutic dose of an antipsychotic drug (olanzapine) were assessed (n = 8). Moreover, psychosis related genes showing methylation differences between schizophrenia-affected twins and healthy co-twins were identified.
A. Monozygotic twins discordant for schizophrenia – Twin Pair 1
Pathways and networks affected by hypermethylation in the patient
Hypermethylated genes in the affected as compared to unaffected twin of twin pair 1 were analyzed with Ingenuity Pathway Analysis (IPA). This analysis identified a number of relevant pathways and networks. These included HIPPO signaling (p = 3.93E-03) and MAPK signaling (p = 4.27E-03) canonical pathways (Table 1). Also, the diseases and disorders affected by hypomethylation of genes in this patient included developmental and hereditary disorders (p = 3.56E-03–4.72E-02). The most significant network affected by hypermethylation of genes in the patient was a cellular growth, proliferation and cellular movement and development network (score = 131) (Fig. 1).
Pathways and networks affected by hypomethylation in the patient
Among the most significant canonical pathways identified using the list of genes that showed promoter hypomethylation in the patient in twin pair 1 were dopamine-DARPP32 feedback in cAMP signaling (p = 7.43E-03), calcium signaling (p = 1.08E-02) and ephrin receptor signaling (p = 1.13E-02) (Fig. 2). Also, the most significant molecular and cellular functions revealed using the same list of genes were gene expression (p = 3.75E-04–4.97E-02), cellular assembly and organization (p = 1.15E-03–4.61E-02), lipid metabolism (p = 3.39E-03–3.69E-02), and small molecule biochemistry (p = 3.39E-03–4.51E-02). There was also clear enrichment for a pathway involved in nervous system development and function (p = 1.71E-03–4.28E-02) (Table 1).
B. Monozygotic twins discordant for schizophrenia – Twin Pair 2
Pathways and networks affected by hypermethylation in the patient
Table 2 presents that the triacylglycerol biosynthesis (p = 2.94E-03) pathway was affected by genes that showed promoter hypermethylation in the patient as compared to their unaffected co-twin. These differences also affected gene expression (p = 5.79E-06–3.62E-02) and molecular transport and protein trafficking (p = 4.01E-03–4.39E-03) functions. The most significant disorders implicated, similarly to twin pair 1, were developmental and hereditary disorders (p = 9.84E-04–4.10E-02), organismal injury and abnormalities (p = 4.41E-04–4.40E-02), and involved aberrations in gene expression (p = 5.79E-06–3.62E-02). The top network affected by hypermethylation of genes included cell death and survival, cellular growth and proliferation, and cell morphology (peak score = 47) (Fig. 3).
Pathways and networks affected by hypomethylation in the patient
It is apparent from the results included in Table 2 that promoter hypomethylation of genes in the patient significantly affects canonical pathways involving glutamate receptor signaling and CREB signaling in neurons (p = 1.53E-02) (Fig. 4), as well as molecular and cellular functions, specifically cell signalling (p = 4.34E-06–4.51E-02), nucleic acid metabolism (p = 4.34E-06–4.51E-02), and small molecule biochemistry (p = 4.34E-06–4.51E-02).
Shared pathways between olanzapine treated rat brains and affected twins
A therapeutic dose of olanzapine caused hyper- and hypo-methylation of genes involved in several psychosis related pathways and networks. Further details on the data are available on Melka et al. . Intriguingly, some of the most significant pathways identified in monozygotic twins discordant for schizophrenia were also affected in rat brain regions. The common pathways and networks affected by olanzapine-induced hypermethylation and the twin patients (exposed to antipsychotic drugs) included the MAPK signaling pathway (p = 1.59E-03), a cell death, cellular growth and proliferation network (peak score = 24) (Table 1, Additional file 1: Tables S1 and S2), and a molecular transport and protein trafficking (peak score = 12) network (Table 2, Additional file 1: Table S2). Similarly, hypomethylation of genes in both sets of results (i.e. olanzapine treated rat brains and the twin study) affected calcium signaling (p = 5.92E-03) (Table 1, Additional file 1: Table S1). We have noted, despite the fact that there are some commonly affected pathways and networks between the two sets of results, that the individual genes involved in hyper/hypomethylation are not identical in rat brains and twin patients. However, both sets of results show methylation changes on psychosis related genes (Fig. 5).
The genome-wide methylation results on two pairs of monozygotic twins discordant for schizophrenia offer a number of insights. First, the twins do differ in methylation patterns. It follows similar conclusions in a number of reports [16, 17]. Interestingly, these differences appear to accumulate over time, and may account for the changing transcriptome  as well as the discordance of monozygotic twins for a variety of phenotypes and disorders, particularly psychosis.
Second, the genes affected by hyper/hypomethylation of promoters are not random. Several genes with promoter hyper/hypomethylation in schizophrenia twins were found to be psychosis related (Fig. 5, Additional file 2: Table S3). Predominantly, these genes appear to play a role in pathways and networks that are relevant in the context of the disease in question. Such conclusions follow a number of reports where discordance for psychosis has been suggested to involve differences in methylation . This study has shown hyper- and hypo-methylation of genes affecting psychotic relevant pathways such as HIPPO signaling and Ephrin signaling, respectively, in twin pair 1. Similarly, hypermethylation of genes affecting Triacylglycerol Biosynthesis which had been associated with psychosis [19, 20]; and hypomethylation of genes affecting glutamate and CREB signaling pathways in twin pair 2, which are also associated with psychosis . Interestingly, neural stem cell proliferation has been previously associated with altered methylation of CREB . In both cases, the identified pathways are involved in disorders that may affect the development and course of psychosis. Further, it has been shown that epigenetic aberrations including methylation, may lead to disturbance of inhibitory activity and cortical associative functions . Interestingly, hippocampal functions and cognitive performance could be influenced by epigenetic regulation of genes implicated in glutamate signaling . The results support the relevance of developing personalized medicine not only with respect to genetic variations but also to epigenetic variations.
Third, the two twin pairs representing discordance for schizophrenia in two unrelated families show overlaps in pathways and interactions identified in the two patients. Of special interest to our findings are networks affected by hypermethylation of genes affecting Cell Death and Survival, Cellular Growth and Cellular Movement (Tables 1 and 2) in both schizophrenia affected patients. Similarly, genome-wide DNA methylation variations between adolescent monozygotic twins were associated with enrichment in pathways and networks involved in cellular growth and proliferation  and these significant differences in methylation of monozygotic twins have been associated with psychosis .
Finally, it is logical to ask the question regarding the cause of these differences. Is DNA methylation causing the disorder or is this due to other factors? To the best of our understanding this question has not been addressed satisfactorily, which is the focus of this section of the discussion. As mentioned, we have generated data on rats to further address this question. It involves the use of a rat model where a therapeutic dose of olanzapine is administered . The results show that olanzapine affects methylation of a large number of genes . They affect pathways and interactions that are compatible with psychosis. For example, we report epigenetic alterations of the dopaminergic system [27, 28] and serotonin transporter gene promoters , which have been evident in major psychiatric disorders, including schizophrenia . In particular, hyper/hypomethylation in rat brain regions support the dopamine hypothesis of schizophrenia reflecting the effect of olanzapine . This argument has been supported by two observations. First, common networks such as cellular growth and cellular movement were affected by hyper/hypomethylation in both patients. Second, psychotic-relevant pathways such as the Dopamine-DARPP32 feedback in cAMP signaling were also affected by olanzapine-induced methylation . Interestingly, those pathways have also been observed in the current twin study. We acknowledge that these results are based on only two pairs of MZ twins and there may be a need to confirm the results using additional methods such as bisulphite sequencing. However, the uniformity of the results being reported prompted us to ask the question, could methylation changes in the affected twin be caused by antipsychotic drugs that have been administered during the course of the management of their psychotic symptoms? We conclude by suggesting that this possibility cannot be ruled out. The implication is that studies on DNA methylation on psychiatric disorders where patients are often on medication should be interpreted carefully. More importantly, it may be premature for published reports to conclude that epigenetic changes are causes of psychosis in patients.
First, DNA methylation differences between monozygotic twins discordant for schizophrenia revealed common networks affected in both unrelated twin pairs suggesting that those networks may underlie the cause of the disease. Although there exist common pathways, this study also showed that unique genes were affected in each pair of twins suggesting that the aetiology and/or the pathophysiology of the diseases may involve unique (patient-specific) genes entailing the need for personalized medicine in treating the disorder. Importantly, these patient-specific genes together with the observed effect of antipsychotics on genes in rat brains as compared to monozygotic twins corroborate the known heterogeneous nature of the disorder. Second, the results also showed patient-specific pathways, which support previous findings highlighting that the effects of antipsychotic drugs administered to patients may have caused at least some of the hyper/hypomethylation differences. These results reflect the fact that personalized medicine is not only necessary for patients of wider genetic variations but also for patients of epigenetic variations such as DNA methylation. The findings could be relevant for future therapeutic investigation of complex disorders in addition to schizophrenia.
Availability of supporting data
Additional files, supporting the observations and discussion, are presented in Table S1, Table S2 and Table S3.
Bosia M, Pigoni A, Cavallaro R. Genomics and epigenomics in novel schizophrenia drug discovery: translating animal models to clinical research and back. Expert Opin Drug Discov. 2015;10(2):125–39.
Abdolmaleky HM, Smith CL, Zhou JR, Thiagalingam S. Epigenetic alterations of the dopaminergic system in major psychiatric disorders. Methods Mol Biol. 2008;448:187–212.
Rosales-Reynoso MA, Ochoa-Hernandez AB, Juarez-Vazquez CI, Barros-Nunez P. Epigenetic mechanisms in the development of memory and their involvement in certain neurological diseases. Neurologia. 2014. doi:10.1016/j.nrl.2014.02.004.
Gavin DP, Akbarian S. Epigenetic and post-transcriptional dysregulation of gene expression in schizophrenia and related disease. Neurobiol Dis. 2012;46(2):255–62.
Razin A, Kantor B. DNA methylation in epigenetic control of gene expression. Prog Mol Subcell Biol. 2005;38:151–67.
Abdolmaleky HM, Smith CL, Faraone SV, et al. Methylomics in psychiatry: modulation of gene-environment interactions may be through DNA methylation. Am J Med Genet B Neuropsychiatr Genet. 2004;127B(1):51–9.
Singh SM, Murphy B, O’Reilly R. Epigenetic contributors to the discordance of monozygotic twins. Clin Genet. 2002;62(2):97–103.
Petronis A, Gottesman II, Kan P, et al. Monozygotic twins exhibit numerous epigenetic differences: clues to twin discordance? Schizophr Bull. 2003;29(1):169–78.
Petronis A. Epigenetics and twins: three variations on the theme. Trends Genet. 2006;22(7):347–50.
Dempster EL, Pidsley R, Schalkwyk LC, et al. Disease-associated epigenetic changes in monozygotic twins discordant for schizophrenia and bipolar disorder. Hum Mol Genet. 2011;20(24):4786–96.
Ollikainen M, Smith KR, Joo EJ, et al. DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Hum Mol Genet. 2010;19(21):4176–88.
Fraga MF, Esteller M. Towards the human cancer epigenome: a first draft of histone modifications. Cell Cycle. 2005;4(10):1377–81.
Melka Melkaye G, Castellani Christina A, Laufer Benjamin I, Rajakumar N, O’Reilly R, Singh SM. Olanzapine induced DNA methylation changes support the dopamine hypothesis of psychosis. J Mol Psych. 2013;1(19):1–7.
Castellani CA, Laufer BI, Melka MG, Diehl EJ, O’Reilly RL, Singh SM. DNA methylation differences in monozygotic twin pairs discordant for schizophrenia identifies psychosis related genes and networks. BMC Med Genomics. 2015;8, 17-015-0093-1. doi:10.1186/s12920-015-0093-1.
Melka MG, Laufer BI, McDonald P, et al. The effects of olanzapine on genome-wide DNA methylation in the hippocampus and cerebellum. Clin Epigenetics.2014; 6(1), 1-7083-6-1. doi:10.1186/1868-7083-6-1.
Dempster EL, Wong CC, Lester KJ, et al. Genome-wide methylomic analysis of monozygotic twins discordant for adolescent depression. Biol Psychiatry. 2014;76(12):977–83.
van Dongen J, Ehli EA, Slieker RC, et al. Epigenetic variation in monozygotic twins: a genome-wide analysis of DNA methylation in buccal cells. Genes (Basel). 2014;5(2):347–65.
Bonsch D, Wunschel M, Lenz B, Janssen G, Weisbrod M, Sauer H. Methylation matters? Decreased methylation status of genomic DNA in the blood of schizophrenic twins. Psychiatry Res. 2012;198(3):533–7.
Foley DL, Mackinnon A, Watts GF, et al. Cardiometabolic risk indicators that distinguish adults with psychosis from the general population, by age and gender. PLoS One. 2013;8(12):e82606.
Chen SF, Hu TM, Lan TH, Chiu HJ, Sheen LY, Loh EW. Severity of psychosis syndrome and change of metabolic abnormality in chronic schizophrenia patients: severe negative syndrome may be related to a distinct lipid pathophysiology. Eur Psychiatry. 2014;29(3):167–71.
Lipina TV, Fletcher PJ, Lee FH, Wong AH, Roder JC. Disrupted-in-schizophrenia-1 Gln31Leu polymorphism results in social anhedonia associated with monoaminergic imbalance and reduction of CREB and beta-arrestin-1,2 in the nucleus accumbens in a mouse model of depression. Neuropsychopharmacology. 2013;38(3):423–36.
Yu M, Li W, Luo S, et al. Folic acid stimulation of neural stem cell proliferation is associated with altered methylation profile of PI3K/Akt/CREB. J Nutr Biochem. 2014;25(4):496–502.
Guidotti A, Auta J, Chen Y, et al. Epigenetic GABAergic targets in schizophrenia and bipolar disorder. Neuropharmacology. 2011;60(7–8):1007–16.
Bagot RC, Zhang TY, Wen X, et al. Variations in postnatal maternal care and the epigenetic regulation of metabotropic glutamate receptor 1 expression and hippocampal function in the rat. Proc Natl Acad Sci U S A. 2012;109 Suppl 2:17200–7.
Levesque ML, Casey KF, Szyf M, et al. Genome-wide DNA methylation variability in adolescent monozygotic twins followed since birth. Epigenetics. 2014;9(10):1410–22.
Byrne EM, Carrillo-Roa T, Henders AK, et al. Monozygotic twins affected with major depressive disorder have greater variance in methylation than their unaffected co-twin. Transl Psychiatry. 2013;3:e269.
Grace AA. Dopamine system dysregulation by the hippocampus: implications for the pathophysiology and treatment of schizophrenia. Neuropharmacology. 2012;62(3):1342–8.
Madras BK. History of the discovery of the antipsychotic dopamine D2 receptor: a basis for the dopamine hypothesis of schizophrenia. J Hist Neurosci. 2013;22(1):62–78.
Abdolmaleky HM, Nohesara S, Ghadirivasfi M, et al. DNA hypermethylation of serotonin transporter gene promoter in drug naive patients with schizophrenia. Schizophr Res. 2014;152(2–3):373–80.
Rivollier F, Lotersztajn L, Chaumette B, Krebs MO, Kebir O. Epigenetics of schizophrenia: a review. Encephale. 2014, doi:10.1016/j.encep.2014.06.005.
This research was supported by grants from the Canadian Institute of Health Research and Dr. Shiva M. Singh held a Senior Research Fellowship from the Ontario Mental Health Foundation during the course of this study.
The authors have declared that they have no competing interests.
MGM and CAC analyzed and interpreted the data. MGM wrote the first draft of the manuscript. CAC, RLO and SMS critically reviewed the manuscript for its intellectual contents; SMS supervised the study and co-wrote the manuscript. All authors have read and approved the final version of the manuscript.
Most significant pathways and networks identified by Ingenuity pathway analysis using genes that showed hypermethylation (a) or hypomethylation (b) in hippocampus, following a therapeutic dose of olanzapine treatment in a rat model in vivo. Table S2. Most significant pathways and networks identified by Ingenuity pathway analysis using genes that showed hypermethylation (a) or hypomethylation (b) in cerebellum, following therapeutic dose of olanzapine treatment in a rat model in vivo.
List of genes hypermethylated or hypomethylated in patient 1 and patient 2 relative to their respective healthy co-twins.