Isogenic mice exhibit sexually-dimorphic DNA methylation patterns across multiple tissues
© The Author(s). 2017
Received: 5 May 2017
Accepted: 27 November 2017
Published: 13 December 2017
Cytosine methylation is a stable epigenetic modification of DNA that plays an important role in both normal physiology and disease. Most diseases exhibit some degree of sexual dimorphism, but the extent to which epigenetic states are influenced by sex is understudied and poorly understood. To address this deficit we studied DNA methylation patterns across multiple reduced representation bisulphite sequencing datasets (from liver, heart, brain, muscle and spleen) derived from isogenic male and female mice.
DNA methylation patterns varied significantly from tissue to tissue, as expected, but they also varied between the sexes, with thousands of sexually dimorphic loci identified. The loci affected were largely autonomous to each tissue, even within tissues derived from the same germ layer. At most loci, differences between genders were driven by females exhibiting hypermethylation relative to males; a proportion of these differences were independent of the presence of testosterone in males. Loci harbouring gender differences were clustered in ontologies related to tissue function.
Our findings suggest that gender is underwritten in the epigenome in a tissue-specific and potentially sex hormone-independent manner. Gender-specific epigenetic states are likely to have important implications for understanding sexually dimorphic phenotypes in health and disease.
Sexual dimorphism, in which the two sexes exhibit different characteristics, affects a range of traits in animals. Almost all human diseases exhibit some component of sexual dimorphism, which can manifest as differences in disease incidence, age of onset, or severity of symptoms . The discordance cannot be completely explained by genetic (i.e. sex chromosome) differences or the actions of sex hormones . Despite increasing efforts in understanding the epigenetic contribution to complex disease , the contribution of epigenetic factors (such as DNA methylation) to sexual dimorphism is relatively underexplored.
There is evidence for gender differences in DNA methylation in various tissues in eutherian mammals. Unsurprisingly, the majority of differentially methylated loci occur on the X-chromosome, where they are hypermethylated in females relative to males . However, sex-specific methylation also occurs on autosomes. A recent meta-analysis of Infinium 450 K array data from 76 individual studies identified gender-specific methylation at about 200 autosomal CpG sites in peripheral blood . Gender-specific methylation differences have also been reported using 450 K array in human prefrontal cortex , saliva , and pancreatic islets , and by reduced representation bisulphite sequencing (RRBS) in mouse liver . A notable recent study used RRBS to find gender-specific methylation at 160 autosomal loci in mouse liver . In this study, a bias towards hypomethylation in males was dependent on testosterone exposure during puberty, which was coincident with loss of methylation. This study suggests that there is an interaction between sex hormones and epigenetic marks, but whether gender-specific marks might exist independent of hormonal or other postnatal factors is unknown.
Given that the bulk of de novo DNA methylation in mammals takes place during the very early stages of embryogenesis, it is possible that at least some gender methylation differences are specified early in development, or are even innate. Such gender differences would be consistent across different tissues, possibly even tissues derived from different germ layers. Here we have explored this possibility by performing RRBS on male and female mouse tissues derived from each of the three germ layers: liver (endoderm), heart (mesoderm), and brain (ectoderm). We find sexually dimorphic, differentially methylated loci in all tissues. By combining our analyses with publically available data, we find that in the liver, a proportion of loci exhibit gender-specific methylation in the absence of testosterone. Very few gender-specific differences are shared among tissues, even those from the same germ layer.
Gender is predicted by DNA methylation patterns in the mouse liver
We generated snapshots of genome-wide cytosine methylation patterns with reduced representation bisulphite sequencing (RRBS). RRBS captures only ~1% of the mammalian genome, but the captured fraction is highly enriched for functional regions such as CpG islands, shores, and other gene regulatory elements. Here we used an enhanced RRBS protocol to produce methylomes from the livers of six female and six male isogenic C57BL/6 J mice. Overall, data quality was very high, as were inter-sample correlations (Additional file 1: Figure S1). For all analyses we considered only those CpGs that were present at >10× coverage in all 12 samples. To combine statistical evidence of neighbouring CpGs, we calculated DNA methylation levels in tiles of 100 bp across the genome (resulting in 820,388 tiles covered across all 12 samples).
We chose a subset of ten DMTs to validate by COBRA . All but one of the ten loci showed a difference in methylation levels between males and females as predicted by the RRBS signal (Additional file 3: Figure S2). This experimental validation indicates that our RRBS and informatics strategy detects gender-specific differences with high confidence.
Gender DMTs in the liver can be testosterone independent
Heart and brain also harbour gender DMTs
Gender DMTs are tissue specific
While this analysis suggests that gender DMTs do not arise in the germline, we considered whether gender DMTs might be common to tissues of the same germ layer origin (i.e. set early in development). We identified DMTs from additional mesodermal tissues studied by Reizel et al. (skeletal muscle and spleen from adult mice) to compare with our mesodermal (heart) tissue. We identified 354 and 48 gender DMTs in muscle and spleen respectively (Additional file 2: Tables S6 and S7), but these DMTs were exclusive to both each other and to the heart (Fig. 4d); this is despite almost all tiles from the muscle and spleen datasets being represented in our heart dataset. Taken together, these comparisons indicate that differential cytosine methylation between genders is almost entirely tissue-specific.
Here we confirm the existence of sex-specific cytosine methylation in the mouse liver and find that tissues derived from all three embryonic germ layers also exhibit sexually dimorphic patterns of methylation that are essentially idiosyncratic to a tissue. In the majority of instances, sex differences were manifest as a strong female bias towards hypermethylation, and this was the case in the liver even when the demethylating action of testosterone in males was removed. Our findings suggest that the DNA methylome undergoes gender differentiation in multiple tissues at some point after lineage specification, in response to tissue-dependent mechanisms.
Tissues of different origins displayed varying extents of gender-specific methylation: the brain harboured more than a thousand differentially methylated regions whereas the heart had only a few hundred. Sex-specific methylation in the brain was remarkable also because it occurred on a very high background of inter-individual variation. Extensive inter-individual variation in the methylome of the human brain has been recognised only recently, and attributed largely to a secondary influence of genetic variation . However in our study all samples were derived from isogenic mice, indicating that at least some proportion of inter-individual variation in methylation patterns in the brain is purely epigenetic in nature, and sex-specific epigenetic signatures are overlaid on this. As our criteria for calling sex-specific methylation was strict (q < 0.01, methylation difference ≥ 25%), it is probable that thousands more sex-specific differences exist in the brain, albeit smaller in magnitude or restricted to specific cell subtypes.
Despite identifying many thousands of gender differences in cytosine methylation overall, we found very few that were in common among tissues. Even tissues originating from the same germ layer displayed a completely autonomous set of gender methylation differences. Taken together with the ontological analysis showing enrichment for tissue-specific pathways, this strongly suggests that gender-specific epigenetic differences are likely to have a tissue-specific function. This idea is supported by our finding that gender-specific methylation differences overlap significantly with histone markers of enhancer regions. It is also consistent with the recent finding that gender-specific methylation in the liver correlates with hepatic gene expression, particularly when the methylation differences occur in regions corresponding to tissue-specific enhancers . The expression of more than a thousand genes in the mammalian liver is sex-dependent, which is not surprising given that sexual dimorphism in both steroid and drug metabolism underpin normal liver physiology . Whether gender biases in methylation in heart and brain also influence their respective transcriptomes is not known, however the tissue-specific functional pathways we find suggest that this may well be the case.
The mechanisms underpinning gender-specific methylation are generally unexplored. However two recent studies have implicated testosterone in the process. In the mouse forebrain, male-specific methylation patterns can be induced in females by neonatal administration of exogenous testosterone . In the liver, testosterone can trigger male-specific demethylation of certain regions; like females, castrated males maintain these regions as methylated . However, not all gender-specific methylation can be attributed to the actions of testosterone: in an unbiased comparison we found more than 300 gender differences between the livers of females and castrated males. While these testosterone-independent gender differences account for a minority of all liver differences, they show the same strong tendency for hypermethylation in females, suggesting that other sex-specific trans-acting factors are also capable of effecting methylation dimorphism. Such factors may involve early and indirect actions of the sex chromosomes, as has been shown for early embryonic gene expression . Female hormones such as estradiol are also capable of skewing gene expression patterns , and so may also be a contributor.
While the weight of our evidence supports widespread tissue-specific methylation differences between genders, our study is not without limitations and there are many unknowns that warrant further investigation. Comparative analysis of many tissues across multiple developmental stages (including embryonic, fetal, neonatal, and adult) would be required to establish the precise timing of establishment of gender differences, and the relationship to the changing hormonal milieu. Furthermore, RRBS allows the interrogation of methylation levels at only a representative fraction of the genome, and while this proportion is enriched for functionally relevant regions, whole genome bisulphite sequencing would be required to capture the full extent of gender-specific DNA methylation. Such studies will only become more feasible with the increasing affordability of whole genome sequencing.
Tissue-specific gender-bias in DNA methylation represents yet another aspect of sexual dimorphism but its function is currently unclear. While likely to reflect normal tissue physiology, such gender bias in epigenetic states may also have broader implications. Intrinsic gender bias in epigenetic state holds the potential to influence disease susceptibility and disease course, and modulate the response to environmental stressors. Gender differences in environmental epigenetic programming are pervasive  and our findings here suggest that at least some of the gender bias in induced phenotypes derives from baseline gender differences in epigenetic state.
Our study of the DNA methylomes of multiple tissues of male and female mice indicates that sex significantly influences DNA methylation patterns in a tissue-specific manner. These findings provide a platform to better understand the role of DNA methylation in health and disease and have important implications for the study of complex and programmed phenotypes. Our findings underscore the need to consider both genders in epigenome-wide association studies, and reinforce the requirements related to the choice of tissue to study.
Mice and tissues
All animals in this study were generated contemporaneously. Mice were handled in accordance with good practice as defined by the National Health and Medical Research Council (Australia) Statement on Animal Experimentation, and requirements of state government legislation. The study was approved by the Garvan/St Vincent’s Animal Ethics Committee (#13/35). Six male and six female C57BL/6 J littermate pairs were selected at weaning and maintained on NIH-31 chow until 24 weeks of age. Tissues (brain, heart, liver) were collected and snap frozen prior to storage at −80°C until DNA extraction. DNA was extracted from a 5 mm coronal section of the rostral end of the brain, a 5 mm apical section of the heart, and the extreme caudal section of the left lobe of the liver.
DNA extraction and RRBS
Frozen tissue sections were homogenised in lysis buffer (50 mM Tris-HCl pH 8.0, 100 mM EDTA, 1% SDS, 100 mM NaCl) and incubated with 400 μg/ml Proteinase K (Roche) overnight at 55°C, followed by phenol:chloroform extraction and ethanol precipitation. Genomic DNA was used for RRBS library preparation and sequencing through the Methyl-MiniSeq service of Zymo Research (Irvine, CA, USA).
Differential methylation analysis of RRBS data
RRBS data from liver heart and brain generated for this study have been deposited in the Gene Expression Omnibus (GEO) under accession number GSE84573. eRRBS and mapping was performed by Zymo Research using proprietary methods. Bed data supplied was reformatted, and identification of differentially methylated cytosines and annotation was carried out in R with the methylKit package, v0.99.2 . To combine statistical evidence of neighbouring CpGs, we calculated DNA methylation levels in tiles of 100 bp across the genome. For all analyses we considered only tiles that were present at > 10× coverage in all replicate samples and < 99.8th percentile of coverage values across all samples in the comparison. Hierarchical clustering was performed using correlation and the Ward’s minimum variance method. Identification of differentially methylated tiles was carried out with methylKit; thresholds for calling a difference were set at a methylation difference of ≥ 25%, and q ≤ 0.01. Genomic annotations of differentially methylated tiles were performed using the mm10 refgene table from UCSC Table Browser. Gene ontology analysis was performed with GREAT  using all default parameters with the exception that distal gene regulatory domains were set at a maximum of 100 kb. Histone modification enrichment analysis was performed by permutation of overlaps with replicated H3K4me1 and H3K27Ac peak files from ENCODE (liver, ENCSR000CDH and ENCSR000CAO; brain, ENCSR000CDF and ENCSR000CAE; heart, ENCSR000CDD and ENCSR000CAI); the overlap of DMTs was compared with the overlap observed using the same number of tiles randomly chosen from the dataset, over 10,000 iterations.
Differential methylation analysis of RRBS data obtained from Reizel et al.
Datasets from liver, castrated liver, muscle and spleen from adult mice (20 weeks) were obtained from GEO Accession GSE60012. Sequencing reads were aligned to the mm10 genome with Bismark  and methylation percentage calls for each CpG site calculated using MethylKit. The sex chromosomes were removed. Our quality control based on pairwise correlations of all samples identified five liver outliers (three female and two male samples) that were excluded from further analysis. Similarly to the analysis of our RBBS data, to combine statistical evidence of neighbouring CpGs, we calculated DNA methylation levels in tiles of 100 bp across the genome, and considered only tiles that were present at > 10× coverage in replicate samples and < 99.8th percentile of coverage values. Given the different number of replicate samples used by Reizel et al. compared to our experimental design, we set this coverage restriction in at least across six samples in each group for liver female vs. male samples comparison, five samples in each group for females vs. castrated liver samples, six samples in each group for muscle female vs. male samples and three samples in each group for the spleen female vs. male samples comparison. Identification of differentially methylated tiles and annotations was carried out with methylKit; thresholds were set at a methylation difference of ≥ 25%, and q ≤ 0.01.
Combined bisulphite restriction analysis (COBRA) validations
Ten loci identified as differentially methylated by RRBS were randomly chosen for independent verification by COBRA, with the additional criterion that they were in the vicinity of coding genes, long non-coding RNA or retrotransposon repeats. Primers were designed using MethPrimer  and sequences are available on request. Following bisulphite PCR and digestion with the relevant restriction enzyme, amplicons representing the fraction of methylated and unmethylated template were quantified using densitometry (Multi Gauge V2.3).
This work was supported by an Australian Research Council Discovery Project grant (DP120100825). HM was supported by an Australian Postgraduate Award. JEC is supported by an Australian Research Council DECRA (DE120100723). CMS was supported by an Australian Research Council Future Fellowship (FT120100097).
Availability of data and materials
The datasets generated in this study are available in the Gene Expression Omnibus (GEO) under accession number GSE84573. Data from Reizel et al. analysed in this study were obtained from GEO under accession GSE60012.
JEC and CMS conceived and oversaw the study; JEC, EG and CMS designed the analytical strategy; KB, HC, JEC and HM performed experiments; EG, HM, PEY analysed data; JEC, EG, HM, SH and CMS interpreted data; JEC, EG, SH, HM and CS drafted the manuscript. All authors read and approved the final manuscript.
All animals in this study were handled in accordance with good practice as defined by the National Health and Medical Research Council (Australia) Statement on Animal Experimentation, and requirements of state government legislation. The study was approved by the Garvan/St Vincent’s Animal Ethics Committee (#13/35).
Consent for publication
Keith Booher and Hunter Chung are employees of Zymo Research; The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ober C, Loisel DA, Gilad Y. Sex-specific genetic architecture of human disease. Nat Rev Genet. 2008;9(12):911–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Morrow EH. The evolution of sex differences in disease. Biol Sex Differ. 2015;6:5.View ArticlePubMedPubMed CentralGoogle Scholar
- Rakyan VK, Down TA, Balding DJ, Beck S. Epigenome-wide association studies for common human diseases. Nat Rev Genet. 2011;12(8):529–41.View ArticlePubMedPubMed CentralGoogle Scholar
- McCarthy NS, Melton PE, Cadby G, Yazar S, Franchina M, Moses EK, Mackey DA, Hewitt AW. Meta-analysis of human methylation data for evidence of sex-specific autosomal patterns. BMC Genomics. 2014;15:981.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu H, Wang F, Liu Y, Yu Y, Gelernter J, Zhang H. Sex-biased methylome and transcriptome in human prefrontal cortex. Hum Mol Genet. 2014;23(5):1260–70.View ArticlePubMedGoogle Scholar
- Liu J, Morgan M, Hutchison K, Calhoun VD. A study of the influence of sex on genome wide methylation. PLoS One. 2010;5(4):e10028.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall E, Volkov P, Dayeh T, Esguerra JL, Salo S, Eliasson L, Ronn T, Bacos K, Ling C. Sex differences in the genome-wide DNA methylation pattern and impact on gene expression, microRNA levels and insulin secretion in human pancreatic islets. Genome Biol. 2014;15(12):522.View ArticlePubMedPubMed CentralGoogle Scholar
- Orozco LD, Rubbi L, Martin LJ, Fang F, Hormozdiari F, Che N, Smith AD, Lusis AJ, Pellegrini M. Intergenerational genomic DNA methylation patterns in mouse hybrid strains. Genome Biol. 2014;15(5):R68.View ArticlePubMedPubMed CentralGoogle Scholar
- Reizel Y, Spiro A, Sabag O, Skversky Y, Hecht M, Keshet I, Berman BP, Cedar H. Gender-specific postnatal demethylation and establishment of epigenetic memory. Genes Dev. 2015;29(9):923–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Akalin A, Kormaksson M, Li S, Garrett-Bakelman FE, Figueroa ME, Melnick A, Mason CE. methylKit: a comprehensive R package for the analysis of genome-wide DNA methylation profiles. Genome Biol. 2012;13(10):R87.View ArticlePubMedPubMed CentralGoogle Scholar
- Xiong Z, Laird PW. COBRA: a sensitive and quantitative DNA methylation assay. Nucleic Acids Res. 1997;25(12):2532–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Illingworth RS, Gruenewald-Schneider U, De Sousa D, Webb S, Merusi C, Kerr AR, James KD, Smith C, Walker R, Andrews R, et al. Inter-individual variability contrasts with regional homogeneity in the human brain DNA methylome. Nucleic Acids Res. 2015;43(2):732–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Waxman DJ, Holloway MG. Sex differences in the expression of hepatic drug metabolizing enzymes. Mol Pharmacol. 2009;76(2):215–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghahramani NM, Ngun TC, Chen PY, Tian Y, Krishnan S, Muir S, Rubbi L, Arnold AP, de Vries GJ, Forger NG, et al. The effects of perinatal testosterone exposure on the DNA methylome of the mouse brain are late-emerging. Biol Sex Differ. 2014;5:8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lowe R, Gemma C, Rakyan VK, Holland ML. Sexually dimorphic gene expression emerges with embryonic genome activation and is dynamic throughout development. BMC Genomics. 2015;16:295.View ArticlePubMedPubMed CentralGoogle Scholar
- Sakakibara M, Uenoyama Y, Minabe S, Watanabe Y, Deura C, Nakamura S, Suzuki G, Maeda K, Tsukamura H. Microarray analysis of perinatal-estrogen-induced changes in gene expression related to brain sexual differentiation in mice. PLoS One. 2013;8(11):e79437.View ArticlePubMedPubMed CentralGoogle Scholar
- Gabory A, Attig L, Junien C. Sexual dimorphism in environmental epigenetic programming. Mol Cell Endocrinol. 2009;304(1-2):8–18.View ArticlePubMedGoogle Scholar
- McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G. GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010;28(5):495–501.View ArticlePubMedPubMed CentralGoogle Scholar
- Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27(11):1571–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Li LC, Dahiya R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002;18(11):1427–31.View ArticlePubMedGoogle Scholar