- Research article
- Open Access
- Open Peer Review
Within-pair differences of DNA methylation levels between monozygotic twins are different between male and female pairs
BMC Medical Genomicsvolume 9, Article number: 55 (2016)
DNA methylation levels will be important for detection of epigenetic effects. However, there are few reports showing sex-related differences in the sensitivity to DNA methylation. To evaluate their sex-related individual differences in the sensitivity to methylation rigorously, we performed a systematic analysis of DNA methylation in monozygotic twins, an optimal model to evaluate them because the genetic backgrounds are the same.
We examined 30 male and 43 female older monozygotic twin pairs recruited from the registry established by the Center for Twin Research, Osaka University. Their methylation levels were determined using the Infinium HumanMethylation450 BeadChip Kit (Illumina), which interrogated 485577 highly informative CpG sites at the single-nucleotide resolution, and the median methylation level was calculated for each of the 25657 CpG islands. Within-pair differences of methylation levels (WPDMs) were greater in male pairs than female pairs for 86.0 % of autosomal CpG islands, but were higher in female pairs than male pairs for 76.7 % of X chromosomal CpG islands. Mean WPDMs of CpG islands in each autosomal chromosome were significantly higher in male pairs than in female whereas that in X chromosome was significantly higher in female pairs than in male. Multiple comparison indicated that WPDMs in three autosomal and two X-chromosomal CpG islands were significantly greater in male pairs, whereas those in 22 X-chromosomal CpG islands were significantly greater in female pairs.
Sex-related differences were present in the WPDMs of CpG islands in individuals with the same genetic background. These differences may be associated with the sexual influences in susceptibility of some diseases.
Human phenotypes, such as physical characteristics, abilities, and disease susceptibility, are determined by both genetic and environmental factors [1–4]. Environmental factors affect human phenotypes by changing the epigenetic modification of the genome, such as by DNA methylation and histone modification . Epigenetic modification changes impact cellular behavior by regulating the chromatin status and gene expression  and so the evaluation of epigenetic changes will be used as new laboratory tests. One of the most important epigenomic modifications is the methylation of genomic DNA, which is the covalent addition of a methyl group to the cytosine at CpG dinucleotides. The CpG sites present in the regions containing high numbers of CpG dinucleotides (CpG islands) are generally unmethylated, although those in the majority of other genomic regions are methylated. CpG islands overlap the promoter regions of 60–70 % of genes and are generally protected from methylation, allowing for the expression of downstream genes, the transcription of which is further regulated by histone modification .
Many reports show the within-pair differences of methylation levels (WPDMs) in discordant monozygotic twins for several disorders and traits because the aberrant DNA methylation of CpG islands may be an important epigenetic change that affects the developmental process of diseases or traits [8–19]. To identify the association of DNA methylation with the development of disease, general WPDMs in monozygotic twin pairs should be assessed. However, they have not yet been elucidated.
In this study, we examined the methylation levels of CpG islands in 113 monozygotic twins, calculated the WPDMs of genomic DNA, and compared the WPDMs between men and women to identify the sex difference in the WPDMs. WPDM of monozygotic twins can reflect the difference of the sensitivity to DNA methylation under the condition of the same genetic background. This study will be able to clarify the sex-related differences in the sensitivity to DNA methylation.
Subjects and Methods
A total of 113 healthy Japanese monozygotic twin volunteers (35 male and 78 female pairs) were recruited from the registry established by the Center for Twin Research, Osaka University (Table 1) . Blood was sampled at 9 am after a 12 h fast. A clinical examination was performed, and the twins completed health-related questionnaires. The twins in each pair were examined on the same day. Genomic DNA was isolated from peripheral blood mononuclear cells using a commercial kit (QIAamp DNA Mini Kit, QIAGEN, Germany). The zygosity of subjects was confirmed by the perfect matching of 15 short tandem repeat (STR) loci using the PowerPlex® 16 System (Promega, Madison, WI, USA).
Methylation level of CpG islands
Analysis of the methylation level was performed using an Infinium HumanMethylation450 BeadChip Kit (Illumina), which interrogated 485577 highly informative CpG sites at the single-nucleotide resolution for each sample using the standard manufacturer's protocol. The experiment was performed with 0.5 μg of high-quality genomic DNA. There were 2 bead types for each CpG site per locus on the chip. The raw data were analyzed using the Genome Studio software (Illumina), and the fluorescence intensity ratios between the 2 bead types were calculated. A ratio value of 0 was equal to the nonmethylation of the locus, and a ratio of 1 was equal to total methylation. These raw data were corrected to normalize the differences in detection ranges between the two probes of the Infinium Assay using a peak-based correction method . Normalized data were filtered to exclude invalid probes, such as null probes and probes with low reliability. After filtering, the data were categorized to each of 25657 CpG islands according to the registration of UCSC [22, 23], and a median methylation level was calculated when there were two or more probes in a CpG island. We used the statistical software R (ver.2.15.1) to perform these data analyses.
Within-pair differences of the methylation level (WPDM)
We calculated the absolute values of differences in each CpG island methylation level between individuals in each pair as follows:
where ML1 is the methylation level of one of each twin pair and ML2 is that of the other twin.
We also calculated the gender difference index of WPDMs in each CpG island as follows
This index is positive when the mean WPDM of a CpG island is higher in a male pair than a female pair.
Student’s t test was used to compare WPDMs between males and females. Statistical analysis was performed using the JMP10 software (SAS Institute, Inc., Tokyo, Japan).
Within-pair differences in the methylation levels (WPDMs) of CpG islands
As shown in Additional file 1: Figure S1, we could find that the WPDMs were larger in many autosomal CpG islands for male pairs than female pairs, whereas the WPDM in many X chromosomal CpG islands were larger in female pairs than male pairs. When we performed the same analysis using only an older subset (>55 years old) (Table 1), we obtained similar results (Fig. 1). As shown in Table 2, means WPDM of CpG islands in each autosomal chromosome were significantly higher in male than in female pairs, whereas that in X chromosome was significantly higher in female than in male pairs. In addition, median of WPDM were also showed the same significances (Table 2).
The WPDMs of CpG islands in older male and female pairs are shown in Additional file 2: Table S1 in ranking order. Table 3 shows the top-rank 50 CpG islands, which have large WPDMs in older male and female pairs, and the common CpG islands, which are included in the top-rank 50 CpG islands of both genders. These are shown in Table 4.
Gender difference index of WPDMs
As shown in Additional file 3: Figure S2, the gender difference indices of WPDMs were positive for 86.0 % (21439/24932) of autosomal CpG islands, but negative for 76.7 % (556/725) of X-chromosomal CpG islands.
Comparison of each WPDM between older male and female pairs
Of the 25657 CpG islands analyzed, 11461 CpG islands showed low P values (<0.05) for WPDMs between male and female pairs using Student’s t test. Among these significant CpG islands, WPDMs in the male pairs were higher in 11027 CpG islands (10975 were autosomal and 52 were X chromosome), whereas those in female pairs were higher in the other 434 islands (51 were autosomal and 383 were X chromosome) (Additional file 4: Table S2). To perform multiple comparisons, we corrected the P values using the Bonferroni method and found 27 significant CpG islands. Of them, 3 were in autosomal chromosomes (2, 8, 12 chromosomes) and 24 were in the X chromosomes (Table 5). The WPDM in male pairs was significantly higher in all 3 autosomal CpG islands (Fig. 2a–c) and 2 of 24 X chromosomal island (Figs. 2d, 2e). Those in the female pairs were significantly higher in 22 of 24 X chromosomal CpG islands (Figs. 3a-v).
We clarified in this study that some CpG islands show large WPDMs both in men and women (Table 4), WPDMs of autosomal CpG islands are generally large in men and those of X-chromosomal CpG islands are generally large in women (Fig. 1, Additional file 1: Figure S1 and Table 2), and multiple comparison indicated the significant differences in WPDMs of some CpG islands between men and women (Table 5) (Figs. 2 and 3). We suppose that these may be caused by the sex-related differences in sensitivity to the DNA methylation or the sex-related difference in the exposure to environment. Therefore, it will be required extra attention to sex-related individual differences when we analyze DNA methylation.
According to the UCSC database [22, 23], the CpG islands with large WPDMs common to both male and female pairs (Table 4) are located near the genes encoding the MC3R (melanocortin 3 receptor), KDM4B (lysine (K)-specific demethylase 4B), ADRBK1 (adrenergic beta receptor kinase 1, also known as GRK2), CCDC28A (coiled-coil domain containing 28A), C9orf171 (chromosome 9 open reading frame 171), LHX8 (LIM homeobox 8), NCOR2 (nuclear receptor corepressor 2), and so on (Table 4). Two of the genes, MC3R and ADRBK1, are related to the regulation of energy homeostasis [24, 25]. Such genes may be susceptible epigenetic changes by environmental factors in both men and women. In addition, these results will serve the data as controls when interpreting the biological relevance of sex-related CpG islands.
In the present study, we found that the WPDMs of most X chromosomal CpG islands are larger in female pairs. This may be due to the random inactivation of the X chromosome, which is specific for females . Interestingly, the WPDMs of most autosomal CpG islands were larger in male pairs. We confirmed these data using older twins because the WPDMs increase with age [27–30]. These indicate that individual differences in most autosomal methylation levels are greater in men than women and suggest that epigenetic changes of DNA in autosomal chromosomes may be more dynamic in men, indicating that men may be more sensitive to environmental factors or may encounter more opportunities to interact with environmental factors compared to women.
It is possible that the large differences in WPDMs of particular gene between men and women may be related to the sex differences in the disease susceptibility of acquired diseases which affected by DNA methylation in that gene. In the present study, statistical analyses indicate that WPDMs were significantly greater in 3 autosomal (Figs. 2a-c) and 2 X chromosomal CpG islands in men (Figs. 2d and e), but were significantly greater in 22 X chromosomal CpG islands in women (Figs. 3a-v). Two of these autosomal CpG islands are located near known genes, ADGRB1 (adhesion G protein-coupled receptor B1) and SLC6A12 (solute carrier family 6 (neurotransmitter transporter) member 12) (Table 5). Interestingly, glioblastoma , gastric cancer , and colorectal cancer , which are dominant in males [34–36], are associated with ADGRB1, and schizophrenia  and autism , which are also dominant in males [39, 40], are associated with SLC6A12.
By contrast, although the WPDMs of the majority of CpG islands in the X chromosome are greater in women, the WPDMs of the two CpG islands in the X chromosome were significantly greater in male pairs. These CpG islands are located near known genes, including ARSD (arylsulfatase D), KCNE1L known as KCNE5 (potassium channel voltage gated subfamily E regulatory beta subunit 5), GYG2 (glycogenin 2), and IRS4 (insulin receptor substrate 4) (Table 5). KCNE1L and ARSD are associated with atrial fibrillation  and gastric dilatation , respectively, both of which are also male dominant [43, 44]. GYG2 is involved in blood glucose homeostasis  and IRS4 encodes the insulin receptor substrate. The CpG sites in such glucose-related genes may be easily influenced by glucose levels, which are higher in men than in women . On the other hand, HCFC1 (host cell factor C1), which has a higher WPDMs in women, is associated with herpes simplex infection , which is female dominant .
Because one of the limitations of this study may be the sample size, which is not enough for high statistical power, there may be some other minor significances we could not find. Another limitation may be a lack of replication study because it is difficult to collect healthy twin data for another cohort. It will be important to analyze the age as co-factor to explore whether the pattern of sex difference changes with age although we could not because of the small sample size. In future, when DNA methylation levels are used as new laboratory tests, our data will be important to know the physiological difference and may also supply significances for diagnosis or prognosis of some sex-related disorders.
In conclusion, sex-related differences were present in the WPDMs of autosomal and X-chromosomal CpG islands, which were greater in men and women, respectively for individuals with the same genetic background. These differences may be associated with the sexual influences in susceptibility of some diseases.
Short tandem repeat
Within-pair differences of the methylation level
Bouchard Jr TJ, McGue M. Genetic and environmental influences on human psychological differences. J Neurobiol. 2003;54:4–45.
Ellis JA, Kemp AS, Ponsonby AL. Gene-environment interaction in autoimmune disease. Expert Rev Mol Med. 2014;16:e4.
Martin GM. Constitutional, somatic genetic and environmental aspects of the phenotypic diversity of aging in human subjects. Basic Life Sci. 1988;43:183–90.
McGue M, Bouchard Jr TJ. Genetic and environmental influences on human behavioral differences. Annu Rev Neurosci. 1998;21:1–24.
Meyer P. Epigenetic variation and environmental change. J Exp Bot. 2015;66:3541–8.
Esteller M. Epigenetics in cancer. N Engl J Med. 2008;358:1148–59.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128:683–92.
Castellani CA, Melka MG, Diehl EJ, Laufer BI, O'Reilly RL, Singh SM. DNA methylation in psychosis: insights into etiology and treatment. Epigenomics. 2015;7:67–74.
Dempster EL, Wong CC, Lester KJ, Burrage J, Gregory AM, Mill J, Eley TC. Genome-wide methylomic analysis of monozygotic twins discordant for adolescent depression. Biol Psychiatry. 2014;76:977–83.
Kratz CP, Edelman DC, Wang Y, Meltzer PS, Greene MH. Genetic and epigenetic analysis of monozygotic twins discordant for testicular cancer. Int J Mol Epidemiol Genet. 2014;5:135–9.
Liu F, Sun Q, Wang L, Nie S, Li J. Bioinformatics analysis of abnormal DNA methylation in muscle samples from monozygotic twins discordant for type 2 diabetes. Mol Med Rep. 2015;12:351–6.
Nilsson E, Jansson PA, Perfilyev A, Volkov P, Pedersen M, Svensson MK, et al. Altered DNA methylation and differential expression of genes influencing metabolism and inflammation in adipose tissue from subjects with type 2 diabetes. Diabetes. 2014;63:2962–76.
Ollikainen M, Ismail K, Gervin K, Kyllonen A, Hakkarainen A, Lundbom J, et al. Genome-wide blood DNA methylation alterations at regulatory elements and heterochromatic regions in monozygotic twins discordant for obesity and liver fat. Clin Epigenetics. 2015;7:39.
Ottini L, Rizzolo P, Siniscalchi E, Zijno A, Silvestri V, Crebelli R, Marcon F. Gene promoter methylation and DNA repair capacity in monozygotic twins with discordant smoking habits. Mutat Res Genet Toxicol Environ Mutagen. 2015;779:57–64.
Roos L, Spector TD, Bell CG. Using epigenomic studies in monozygotic twins to improve our understanding of cancer. Epigenomics. 2014;6:299–309.
Schreiner F, Gohlke B, Stutte S, Bartmann P, Hecher K, Oldenburg J, et al. 11p15 DNA-methylation analysis in monozygotic twins with discordant intrauterine development due to severe twin-to-twin transfusion syndrome. Clin Epigenetics. 2014;6:6.
Selmi C, Cavaciocchi F, Lleo A, Cheroni C, De Francesco R, Lombardi SA, et al. Genome-wide analysis of DNA methylation, copy number variation, and gene expression in monozygotic twins discordant for primary biliary cirrhosis. Front Immunol. 2014;5:128.
Tan Q, Frost M, Heijmans BT, von Bornemann HJ, Tobi EW, Christensen K, Christiansen L. Epigenetic signature of birth weight discordance in adult twins. BMC Genomics. 2014;15:1062.
Yuan W, Xia Y, Bell CG, Yet I, Ferreira T, Ward KJ, et al. An integrated epigenomic analysis for type 2 diabetes susceptibility loci in monozygotic twins. Nat Commun. 2014;5:5719.
Hayakawa K, Iwatani Y and Osaka Twin Research G. An overview of multidisciplinary research resources at the Osaka University Center for Twin Research. Twin Res Hum Genet. 2013;16:217–20.
Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F. Evaluation of the Infinium Methylation 450 K technology. Epigenomics. 2011;3:771–84.
The University of California. The Genome Browser. https://genome.ucsc.edu/. Accessed 22 Dec 2015
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D. The human genome browser at UCSC. Genome Res. 2002;12:996–1006.
Begriche K, Girardet C, McDonald P, Butler AA. Melanocortin-3 receptors and metabolic homeostasis. Prog Mol Biol Transl Sci. 2013;114:109–46.
Vila-Bedmar R, Garcia-Guerra L, Nieto-Vazquez I, Mayor Jr F, Lorenzo M, Murga C, Fernandez-Veledo S. GRK2 contribution to the regulation of energy expenditure and brown fat function. FASEB J. 2012;26:3503–14.
Gendrel AV, Heard E. Noncoding RNAs and epigenetic mechanisms during X-chromosome inactivation. Annu Rev Cell Dev Biol. 2014;30:561–80.
Hannum G, Guinney J, Zhao L, Zhang L, Hughes G, Sadda S, et al. Genome-wide methylation profiles reveal quantitative views of human aging rates. Mol Cell. 2013;49:359–67.
Heyn H, Li N, Ferreira HJ, Moran S, Pisano DG, Gomez A, et al. Distinct DNA methylomes of newborns and centenarians. Proc Natl Acad Sci U S A. 2012;109:10522–7.
Horvath S. DNA methylation age of human tissues and cell types. Genome Biol. 2013;14:R115.
Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, Beyan H, et al. Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010;20:434–9.
Zhu D, Hunter SB, Vertino PM, Van Meir EG. Overexpression of MBD2 in glioblastoma maintains epigenetic silencing and inhibits the antiangiogenic function of the tumor suppressor gene BAI1. Cancer Res. 2011;71:5859–70.
Lee JH, Koh JT, Shin BA, Ahn KY, Roh JH, Kim YJ, Kim KK. Comparative study of angiostatic and anti-invasive gene expressions as prognostic factors in gastric cancer. Int J Oncol. 2001;18:355–61.
Fukushima Y, Oshika Y, Tsuchida T, Tokunaga T, Hatanaka H, Kijima H, et al. Brain-specific angiogenesis inhibitor 1 expression is inversely correlated with vascularity and distant metastasis of colorectal cancer. Int J Oncol. 1998;13:967–70.
Kleihues P, Ohgaki H. Primary and secondary glioblastomas: from concept to clinical diagnosis. Neuro Oncol. 1999;1:44–51.
Buas MF, Vaughan TL. Epidemiology and risk factors for gastroesophageal junction tumors: understanding the rising incidence of this disease. Semin Radiat Oncol. 2013;23:3–9.
Koo JH, Leong RW. Sex differences in epidemiological, clinical and pathological characteristics of colorectal cancer. J Gastroenterol Hepatol. 2010;25:33–42.
Park HJ, Kim JW, Lee SK, Kim SK, Park JK, Cho AR, et al. Association between the SLC6A12 gene and negative symptoms of schizophrenia in a Korean population. Psychiatry Res. 2011;189:478–9.
Silva IM, Rosenfeld J, Antoniuk SA, Raskin S, Sotomaior VS. A 1.5Mb terminal deletion of 12p associated with autism spectrum disorder. Gene. 2014;542:83–6.
Markham JA. Sex steroids and schizophrenia. Rev Endocr Metab Disord. 2012;13:187–207.
Werling DM, Geschwind DH. Sex differences in autism spectrum disorders. Curr Opin Neurol. 2013;26:146–53.
Ravn LS, Aizawa Y, Pollevick GD, Hofman-Bang J, Cordeiro JM, Dixen U, et al. Gain of function in IKs secondary to a mutation in KCNE5 associated with atrial fibrillation. Heart Rhythm. 2008;5:427–35.
Haskins ME, Desnick RJ, DiFerrante N, Jezyk PF, Patterson DF. Beta-glucuronidase deficiency in a dog: a model of human mucopolysaccharidosis VII. Pediatr Res. 1984;18:980–4.
Benjamin EJ, Levy D, Vaziri SM, D'Agostino RB, Belanger AJ, Wolf PA. Independent risk factors for atrial fibrillation in a population-based cohort. The Framingham Heart Study. JAMA. 1994;271:840–4.
Jennings Jr PB, Butzin CA. Epidemiology of gastric dilatation-volvulus in the military working dog program. Mil Med. 1992;157:369–71.
Mu J, Skurat AV, Roach PJ. Glycogenin-2, a novel self-glucosylating protein involved in liver glycogen biosynthesis. J Biol Chem. 1997;272:27589–97.
Soeters MR, Sauerwein HP, Groener JE, Aerts JM, Ackermans MT, Glatz JF, et al. Gender-related differences in the metabolic response to fasting. J Clin Endocrinol Metab. 2007;92:3646–52.
Vogel JL, Kristie TM. The dynamics of HCF-1 modulation of herpes simplex virus chromatin during initiation of infection. Viruses. 2013;5:1272–91.
Liesegang TJ. Herpes simplex virus epidemiology and ocular importance. Cornea. 2001;20:1–13.
The authors are thankful to all of the consultants from the Osaka Twin Research Group (Yoshinori Iwatani (lead author, firstname.lastname@example.org), Shiro Yorifuji, Hiroyasu Iso, Kei Kamide, Jun Hatazawa, Shinji Kihara, Norio Sakai, Hiroko Watanabe, Kiyoko Makimoto, Mikio Watanabe, and Chika Honda, Center for Twin Research, Osaka University Graduate School of Medicine) and all of the technical and secretarial staff of the Center for Twin Research, Osaka University Graduate School of Medicine. The authors are also thankful to Beckman Coulter, Inc. (Tokyo, JAPAN) for collaborative studies.
This project was supported by University Grants from the Japanese Ministry of Education, Culture, Sports, Science and Technology, and also by JSPS KAKENHI Grant Number 24590695, and also by the Charitable Trust Laboratory Medicine Research Foundation of Japan.
Availability of data and material
The data supporting our findings can be found in Additional file 5: Table S3.
MW, CH, and YI conceived and designed the experiments. MW analyzed the data. MW and YI interpreted the results and wrote the paper. All authors reviewed and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Written consent for publication was obtained from all of the twins.
Ethics approval and consent to participate
Written informed consent was obtained from all of the twins, and the Ethics Committee of Osaka University approved the study protocol (No. 506).
Within-pair differences for the methylation levels (WPDMs) of each CpG island. Red circles indicate male pairs, and blue circles indicate female pairs. Within-pair differences in male pairs are greater in most autosomal CpG islands. (TIF 2145 kb)
Rank order within-pair differences in methylation levels (WPDMs) of CpG islands in older male and female pairs (in descending order). (XLSX 2065 kb)
Comparison of within-pair differences in methylation levels (WPDMs) between older male and older female pairs. The gender difference index is positive when the mean within-pair differences in the methylation levels are higher in male pairs than female pairs. (TIF 298 kb)
Results of statistical test comparing within-pair difference of methylation levels (WPDMs) between older male and female twin pairs. WPDMs in the male pairs were higher in 11027 CpG islands (10975 were autosomal and 52 were X chromosome), whereas those in female pairs were higher in the other 434 islands (51 were autosomal and 383 were X chromosome). (XLSX 1758 kb)
Data supporting the findings in this study. (XLSX 6571 kb)