- Open Access
X chromosome dosage and presence of SRY shape sex-specific differences in DNA methylation at an autosomal region in human cells
Biology of Sex Differencesvolume 9, Article number: 10 (2018)
Sexual dimorphism in DNA methylation levels is a recurrent epigenetic feature in different human cell types and has been implicated in predisposition to disease, such as psychiatric and autoimmune disorders. To elucidate the genetic origins of sex-specific DNA methylation, we examined DNA methylation levels in fibroblast cell lines and blood cells from individuals with different combinations of sex chromosome complements and sex phenotypes focusing on a single autosomal region––the differentially methylated region (DMR) in the promoter of the zona pellucida binding protein 2 (ZPBP2) as a reporter.
Our data show that the presence of the sex determining region Y (SRY) was associated with lower methylation levels, whereas higher X chromosome dosage in the absence of SRY led to an increase in DNA methylation levels at the ZPBP2 DMR. We mapped the X-linked modifier of DNA methylation to the long arm of chromosome X (Xq13-q21) and tested the impact of mutations in the ATRX and RLIM genes, located in this region, on methylation levels. Neither ATRX nor RLIM mutations influenced ZPBP2 methylation in female carriers.
We conclude that sex-specific methylation differences at the autosomal locus result from interaction between a Y-linked factor SRY and at least one X-linked factor that acts in a dose-dependent manner.
Recent genome-wide studies of DNA methylation in human cells have detected hundreds to thousands of autosomal loci where methylation levels differ between the sexes [1,2,3,4,5]. It has been proposed that sex-specific differences in methylation levels contribute to sexual dimorphism in complex phenotypes such as early cognitive processes , pancreatic islet function , predisposition to asthma , and psychiatric disorders . While differences in methylation of X-linked genes are explained by X-chromosome inactivation, the cause of sex-specific methylation of autosomal loci remains elusive with the two main suspects: the sex chromosome complement and gonadal hormones. The goal of this study was to determine if the sex chromosomes contribute to sexual dimorphism in autosomal DNA methylation in human cells.
Mammalian males and females have different sex chromosome complements (XX in females and XY in males). The sex chromosome complement determines gonadal sex, which in turn defines production of gonadal hormones. The ensemble of current evidence suggests that both major factors, the sex chromosome complement and gonadal hormones either acting separately or in concert lead to sexual dimorphism in a wide array of phenotypes (reviewed in [7, 8]).
Sex-specific differences in phenotypes and gene expression patterns have been documented in different species from Drosophila to humans ([9,10,11,12]. Sexually dimorphic expression is seen in genes residing on the sex chromosomes, as well as autosomal genes. Expression of certain autosomal genes depends on the sex chromosome dosage, whereas other autosomal genes are controlled by sex hormones, as demonstrated in mice and humans [13,14,15,16]. In mice, about 10% of autosomal genes have different expression levels between males and females in somatic cells . Most of these differences depend on the sex of the mouse; however, expression levels of hundreds of genes depend on the sex chromosome complement . Sexual dimorphism in autosomal gene expression may result from different doses of transcription factors that are encoded by the sex chromosome-linked genes, especially those X-linked genes that escape X-inactivation. It may also arise as a consequence of sex-specific differences in methylation levels of regulatory elements. In humans, most sexually dimorphic autosomal differentially methylated regions have higher methylation levels in females [1, 2, 18], whereas repetitive elements have higher methylation in males [2, 19,20,21,22,23]. In females with Turner syndrome and monosomy X, global DNA methylation levels are lower than in 46,XX females, but higher than in 46,XY males [24, 25]. This suggests that both sex chromosomes as well as gonadal hormones contribute to the establishment or maintenance of global DNA methylation patterns.
Variation in DNA methylation often correlates with variation in chromatin modifications. Enhancer-specific active chromatin marks, H3K27ac and H3K4me1, co-localize with loci that show sex-specific differences in methylation in the mouse liver . Moreover, the analysis of combined data from the ENCODE and Roadmap Epigenomics projects suggests that sex-specific differences in epigenetic marks are more likely to occur at enhancer regions and “bivalent” regions that carry both active and repressive chromatin marks .
To elucidate the genetic mechanisms contributing to sexual dimorphism in DNA methylation, we asked whether the sex chromosome complement influenced DNA methylation levels in human cells (Additional file 1: Figure S1). To address this question, we focused on a single autosomal differentially methylated region (DMR) located within the promoter/enhancer region of the human zona pellucida binding protein 2 gene (ZPBP2) and used it as a marker DMR. DNA methylation of the ZPBP2 DMR shows robust differences between males and females in different tissues/cell types [2, 4, 6, 28]. Here, we show that both the Y and X chromosomes control ZPBP2 methylation levels.
DNA samples from the peripheral blood of five carriers of Xq13-q21 duplications with known boundaries of duplicated regions [29, 30] were used for methylation analysis. Three carriers of ATRX, chromatin remodeler (ATRX) mutations c.7366_7367insA , eight heterozygous carriers of ring finger protein, LIM domain interacting (RLIM) gene mutations p.Arg387Cys, and six non-carrier family members were also included in the study .
The DNA samples from males and females from the Saguenay-Lac-Saint-Jean (SLSJ) asthma familial collection were used as additional controls to compare methylation levels between males, females, and female carriers of X chromosomal duplications and mutations in X-linked genes [6, 28]. We therefore used ZPBP2 methylation data from previously published works as a reference for methylation levels in peripheral blood cells (PBC).
All experiments were conducted in accordance with the Canadian Tri-Council Policy Statement of Ethical Conduct for Research Involving Humans and approved by the REB of the MUHC.
DNA samples from fibroblast cell lines
The DNA samples extracted from untransformed human fibroblast cell lines were purchased from the Coriell Cell Repository (Camden, New Jersey, USA). Information on karyotype, genomic positions of translocation/deletion breakpoints, sex, age, number of passages, and the presence or absence of sex-determining region Y (SRY) were obtained from the Coriell Institute for Medical Research website . If no SRY information was available, the samples were genotyped for SRY by PCR (primers are listed in Additional file 2: Table S1). Cell lines and their characteristics are listed in Additional files 3 and 4: Tables S2 and S3.
DNA methylation assays
The DNA samples were treated with sodium bisulfite using the EpiTect Bisulfite Kit (Qiagen) per manufacturer’s protocol with modifications. The incubation step was extended by one cycle of 5 min at 95 °C and 2 h at 60 °C before holding at room temperature overnight until purification. The samples were purified according to the manufacturer’s protocol, dried at 56 °C for 5 min, and eluted in 40 μl (2 × 20 μl) of the elution buffer supplied with the kit.
Pyrosequencing methylation assays for the ZPBP2 promoter/enhancer region were carried as previously described [6, 28]. The assay interrogates 11 CGs located near the transcriptional start site of ZPBP2 (Additional file 5: Figure S2). This region is referred to as ZPBP2 DMR from this point on. Methylation was assayed using the PyroMark Q24 Advanced platform (Qiagen) and PyroMark Q24 Advanced CpG Reagents. The assay was conducted using separate sequencing primers (Additional file 2: Table S1). The results were analyzed by the PyroMark Q24 Advanced 3.0.0 software (Qiagen). Sodium bisulfite sequencing assays to interrogate methylation of 51 CGs of the ZPBP2 promoter region were conducted as previously described  (Additional file 5: Figure S2). The ZPBP2 DMR targeted by pyrosequencing overlaps with CGs 26–35 of the sodium bisulfite sequencing ZPBP2 methylation assay. Sanger sequencing of clones was performed by the McGill University and Genome Quebec sequencing platform.
To study the influences of phenotypic sex and sex chromosome complement on the methylation patterns in these cell lines, ZPBP2 DMR methylation was analyzed in 36 fibroblast cell lines: 13 cell lines from 46,XX females, 7 from 46,XY males, 3 from 46,XX males, 4 from 46,XY females, 7 from 45, X females, one from a 49,XYYYY male, and one from a 48/49,XXXY/XXXXY male (Additional file 3: Table S2). The ZPBP2 pyrosequencing methylation assay interrogates 11 CGs . Inconsistent quality measures were found at CG1, and therefore, 10 CGs were used for further analyses. The data were missing for at least one of the 10 CGs in seven of the 36 assayed cell lines. Therefore, only 29 fibroblast cell lines were included in the principal component (PC) analysis (Additional file 3: Table S2 and Additional file 6: Table S4).
To improve the signal to noise ratio, the first PC of the 10 CGs was calculated to summarize the methylation data. This PC was plotted against methylation data and multiplied by −1 to ensure that the PC was positively correlated with general methylation trends; the resulting oriented first PC was used as the dependent variable in subsequent linear models. Since the methylation levels at the 10 CGs are correlated, to perform separate analyses would lead to unnecessary multiple tests whereas the principal component captures the primary trend across the region.
Considering the available distribution of X and Y chromosome numbers, observations were grouped to compare samples from cell lines with one X chromosome (coded as “0”) to samples with two or more X chromosomes (coded as “1”). The presence of the SRY locus was coded as “1,” and absence as “0” (see Table 1). Relationships between the first PC of methylation and sex phenotype, presence of the Y chromosome, presence of the SRY region, number of X chromosomes, donor’s age, and number of passages were first explored descriptively by fitting separate linear regressions. Then, stepwise regression was used to obtain our final model, with model selection criterion based on likelihood ratio tests, with p = 0.10 for variable removal and p = 0.05 for variable inclusion. Likelihood ratio tests were used instead of Wald tests or score tests since the asymptotic convergence is often better in small sample sizes . The starting model used for stepwise regression included covariates: sex, age, passage number, number of Y chromosomes, as well as main effects, and the interaction for the presence of SRY and number of X chromosomes.
The presence of the SRY region and the number of X chromosomes influence ZPBP2 methylation levels.
To test the impact of the sex-chromosome complement on autosomal DNA methylation levels in human cells, as a first step, we used human non-transformed fibroblast cell lines derived from donors with different combinations of sex-chromosome complement and sex that would permit separating the effects of sex chromosomes from those of sex phenotype (Additional file 3: Table S2 and Additional file 1: Figure S1). We selected the ZPBP2 DMR located within the ZPBP2 promoter as a marker of sex-specific differences in methylation levels. Previous studies showed robust sex-specific differences in DNA methylation in PBC, higher methylation levels in females, and differences in methylation greater than 5% in this region [2, 4, 6, 28]. We analyzed ZPBP2 DMR methylation levels in 36 fibroblast cell lines with different combinations of sex chromosome complement and sex phenotype (Fig. 1a). Of these, 29 cell lines had complete data for the 10 CGs of the ZPBP2 DMR and were used for PCA: these included 11 cell lines from 46,XX females; 4––from 46,XY males; 3––from 46,XX males; 3––from 46,XY females, 6––from 45,X individuals; one from a 49,XYYYY male, and one from a 48/49,XXXY/XXXXY male (Additional file 7: Figure S3 and Additional file 6: Table S4).
Methylation levels of the ZPBP2 DMR in fibroblast cell lines, summarized by their first PC, were analyzed to examine different combinations of sex phenotype and sex chromosome complement (Fig. 1b). The first PC captures 59.2% of the total variation (Additional file 7: Figure S3). After grouping the independent variables of interest, there were 14 samples with one X chromosome and 15 with two or more X chromosomes, and there were 15 samples with SRY and 14 samples without SRY.
A number of initial linear models were fit to provide a general description of the effects of potential covariates on methylation. When sex phenotype, passage number, number of Y chromosomes, and age were, one at a time, regressed on the first PC of methylation, none of them showed a significant association (p = 0.88 for sex, p = 0.25 for passage number, p = 0.21 for number of Y chromosomes, p = 0.65 for age; see also Additional file 8: Table S5). In contrast, a linear model that regressed, individually, the number of X chromosomes and SRY on the first PC of methylation showed significant associations (p = 0.019 for number of X chromosomes, p = 0.037 for SRY).
Stepwise regression, using likelihood ratio tests, to choose variables for inclusion and exclusion selected a final model with covariates sex, SRY, number of X chromosomes, and an interaction between the latter two variables (Table 2 and Fig. 1b). All multiple regression models excluded one sample with ambiguous sex, and therefore, sample size was 28 in these models. The significant interaction between the number of X chromosomes and SRY presence on the methylation patterns is such that the relationship between methylation and the number of X chromosomes is mainly evident when SRY is absent (Fig. 1b). In the absence of SRY, methylation was higher in the samples with two or more X chromosomes, when compared to a reference group of samples with one X chromosome (Fig. 1b). In the presence of SRY, however, the number of X chromosomes appears to have little influence on ZPBP2 methylation (Table 3), although only three samples with SRY and two or more X chromosomes were available for analysis. We realize that due to the numbers of non-transformed fibroblast samples available for this analysis that power is limited (see Table 1 and Table S5); nevertheless, given that only one multiple linear model was fit, no adjustments for multiple comparisons were performed.
Stepwise multiple linear regression results, chosen by comparing model fits with likelihood ratio tests, between the first PC of methylation data at ZPBP2 and covariates. Estimates and standard errors are in units of methylation percentages. The number of samples in all analyses is 28. This model has an adjusted R-squared of 0.30 and a naïve R-squared of 0.41. Note that the p value of 0.054 for “sex” is obtained from a Wald test; this variable was retained in the model based on a likelihood ratio test with p = 0.032.
The results are estimated from the multiple linear regression in Table 2 and averaged over the two sexes.
The X-linked modifier of DNA methylation maps to chromosomal region Xq13-Xq21
Since the X chromosome dosage influenced methylation, the next step was to map the region of the X chromosome harboring the gene(s) responsible for this influence. We reasoned that the X-linked modifier of DNA methylation should satisfy the following criteria: (i) two or more copies of the modifier would result in higher methylation levels compared to one copy and (2) loss of function mutations in the modifier in 46,XX females would result in lower methylation levels at the reporter locus.
We first asked if the X-linked modifier resided on the Xp or Xq arms and compared ZPBP2 methylation in fibroblast cell lines from four Turner syndrome patients with monosomy X (45, X), four individuals with monosomy for most of the p-arm and trisomy for the q-arm of X chromosome (referred to as 46,X,i(Xq)) and three females with trisomy X (47 or 48,XXX) (Fig. 2a). Since the SRY region may reduce methylation levels in samples with two or more X chromosomes, we selected cell lines in which SRY was absent. The methylation levels at individual CGs within the ZPBP2 DMR were correlated (Pearson correlation coefficient ranging from 0.36 to 0.85). Hence, we used the average methylation levels across the 10 CGs as a read-out. Fibroblast cell lines from Turner syndrome patients with monosomy X had significantly lower methylation levels than those with three copies of the Xq (p = 0.00013, one-tailed Student’s t test) (Fig. 2b). A similar trend was observed in the extended ZPBP2 promoter region (Additional file 5: Figure S2B). Therefore, at least one of the genes that influence ZPBP2 DMR methylation resides on the long arm of the X chromosome.
To test the possibility that the number of inactive X chromosomes rather than the dose of a particular gene influenced ZPBP2 DMR methylation levels, we compared females with two inactive X chromosomes (Xi) and three copies of Xq (47/48,XXX) to those with one Xi and three copies of Xq (46,X,i(Xq)). No significant differences in methylation levels were found between these two groups (p = 0.13, one-tailed Student’s t test). Methylation in 46,X,i(Xq) samples was not statistically different from that in 46,XX females (p = 0.49, two-tailed Student’s t test).
Next, to further reduce the candidate region for the Xq-linked methylation modifier(s), we examined four fibroblast cell lines from females with different Xq deletions (Fig. 2a). ZPBP2 methylation levels in these cell lines were higher than in 45,X Turner syndrome patients (Fig. 2b). We hypothesized that all four cell lines with Xq deletions contained two or more copies of the X-linked modifier. Two Xq regions, referred to as region 1 and region 2 from this point on, were present in all four cell lines in more than one copy (Fig. 2a). Therefore, if there was only one methylation modifier in Xq, it had to reside in one of these two regions. Region 1 spans more than 10 Mb between the centromere and genomic position 78.6 Mb (hg19) (Xq11-q21) and includes the X-inactivation center. Region 2 spans a 2.2 Mb genomic interval between positions 106.6 and 108.8 Mb (Xq22.3) of the X chromosome and includes 18 annotated genes.
We next hypothesized that if region 1 harbored the X-linked modifier of methylation, female carriers of chromosomal duplications of this region would have higher methylation levels than females without duplications. There were no fibroblast cell lines with duplications of the Xq11-q21 region available to us. We therefore assayed ZPBP2 methylation in peripheral blood of three females who carried duplications of chromosomal region Xq13.2-q21.1 [29, 30] and compared their methylation levels to those observed in the blood of control females and males from the SLSJ familial collection  (Fig. 2c). Methylation in the duplication carriers was at the higher end of the range observed in non-carrier females, and the difference reached statistical significance in our sample (p = 0.0045, two-tailed Student’s t test) (Fig. 2c). Based on the boundaries of these duplications, we speculate that the candidate gene may reside within a 5.2 Mb genomic interval between genomic positions 73.1 and 78.4 (hg19) (Fig. 2a); however, this conclusion must be validated using a larger sample size.
Two genes, the ATRX, chromatin remodeler (ATRX) and ring finger protein, LIM domain interacting (RLIM) reside in region 1 and are of particular interest because they encode transcription factors that are involved in epigenetic silencing and X chromosome inactivation [36,37,38,39,40]. These genes are subject to X-inactivation , but could act before its onset in early embryos. In such a case, the mutation would have an effect independent of whether it resided on the X chromosome that was preferentially inactivated. We therefore hypothesized that if ATRX or RLIM influenced ZPBP2 methylation, their mutations would reduce methylation levels in heterozygous carrier females. ZPBP2 methylation levels were assayed in two female carriers of the disease-causing ATRX mutation c.7366_7367insA  and eight heterozygous carrier females of the disease-causing mutation p.Arg387Cys in the zinc-finger domain of RLIM . The samples from six female relatives of the RLIM mutation carriers were also assayed. ZPBP2 methylation levels in female carriers of RLIM and ATRX mutations were not different from those found in controls (Fig. 2c). Thus, neither ATRX nor RLIM mutations reduced the methylation levels of ZPBP2 in heterozygous females in our study.
Here, we report that X-chromosome dosage and sex determining factor SRY region contribute to sexual dimorphism in DNA methylation at an autosomal region, the promoter/enhancer of the ZPBP2 gene located in chromosomal region 17q12, whereas phenotypic sex and age have no significant effect on methylation. The presence of the whole Y-chromosome or its part harboring the SRY gene were associated with lower methylation levels, independent of the number of X chromosomes. To the best of our knowledge, this is the first study of the effect of sex chromosomes and sex phenotype on autosomal methylation in human cells.
In the DNA samples from Turner syndrome patients that carried only one X, the ZPBP2 DMR methylation levels were lower than in individuals who carried two or more copies of intact X chromosomes, or the Xq arms. Our findings are consistent with hypomethylation of the ZPBP2 promoter region in the blood cells of Turner syndrome patients compared to 46,XX females [24, 25]. The caveat here is that our conclusions are based on the assumption that rearrangements and aneuploidy for chromosomes other than the sex chromosomes are unlikely to affect ZPBP2 methylation levels. This assumption has not been tested.
To identify the X-linked genes responsible for the levels of DNA methylation at the ZPBP2 locus, we compared methylation in fibroblast cell lines from Turner syndrome patients with karyotypes 45,X and 46,X,i(Xq), females with trisomy X, as well as females with deletions within the Xq arm. To confirm the effect of the Xq13-q21 dosage on ZPBP2 methylation, we used the DNA samples from peripheral blood cells of female carriers of Xq13-q21 duplications. The striking differences in methylation between Turner syndrome patients with one and three copies of Xq suggest that at least one X-linked modifier of ZPBP2 methylation levels resides on the q-arm of the X chromosome. ZPBP2 methylation in the samples from Xq13-q21 duplication carriers was at the higher end of the range of controls, but the difference was not as dramatic as the difference between males and females. Therefore, using samples from female carriers of X chromosomal deletions may prove to be a more informative mapping approach.
The results of testing two candidate genes from region 1, ATRX and RLIM, did not support the hypothesis that the dosage of either one of these genes was responsible for the higher methylation of the ZPBP2 DMR in females. However, it is possible that a missense mutation in ATRX or RLIM may impact DNA methylation differently than a deletion of the entire gene would.
In principle, the X-linked gene(s) that is responsible for the dose-dependent effect on DNA methylation should either escape X-chromosome inactivation or act early in embryonic development before the onset of X-inactivation (reviewed in ). In humans, sex differences in DNA methylation are present in newborns  whereas the earliest developmental stage with sex differences in human gene expression is represented by embryonic stem cells . Sex differences in gene expression are detected in embryonic stem cells derived from mouse blastocysts and in bovine blastocysts and as early as 8-cell stage in mouse embryos [44,45,46]. Hence, one cannot rule out the possibility that the causal X-linked gene contributes to the establishment of sexually dimorphic DNA methylation patterns before the blastocyst stage, and silencing of one of its copies later in development does not alter DNA methylation patterns. Alternatively, if the sexual dimorphism in DNA methylation is established after the onset of X-inactivation or the X-linked factor is required for the maintenance of higher methylation levels in females, the causal X-linked gene must escape X-inactivation to act in a dose-dependent fashion. Indeed, the X-linked gene for lysine demethylase 5C (KDM5C) that escapes X-inactivation and its Y-homolog KDM5D, both influence DNA methylation at certain autosomal regions in a dose-dependent fashion . KDM5C resides on the short arm of the X chromosome and in principle may be responsible for the tendency to lower methylation in 46,X,i(Xq) compared to trisomy X females.
None of region 2 genes are known to escape X-inactivation. Region 1 is larger, and several of its genes escape X-inactivation including several long non-coding RNAs (e.g., X-inactive specific transcript (XIST), JPX transcript, XIST activator (JPX)) [41, 48]. Region 1 also harbors several micro RNAs whose inactivation status is not known [41, 48]. This warrants further studies of the role of region 1 genes in the establishment of autosomal methylation patterns.
X-chromosome inactivation requires epigenetic silencing of hundreds of genes on the same X chromosome. It has been suggested that the presence of an inactive X per se may influence the genomic distribution of heterochromatic proteins, acting as an epigenetic “sink” or “tank” [13, 49,50,51]. The data from the studies of meiotic silencing in the mammalian germ line show that, indeed, sex chromosomes and autosomes may compete for silencing factors in certain cell types [52,53,54,55] landing additional support for this hypothesis. If this were the mechanism by which X chromosome dosage affected autosomal methylation, it would cause differences in DNA methylation levels between individuals who carry an inactive X and those who do not. Moreover, one would predict that the number of inactive X chromosomes rather than number of copies of a particular X-linked gene would be critical. That is, the presence of an inactive isochromosome Xq would have the same impact on methylation as would one Xi or an Xi with a duplication have. Although our data are consistent with the possibility that the inactive X promotes higher DNA methylation levels at the ZPBP2 DMR (acting as a “supply center” rather than a “sink” for heterochromatic factors), higher methylation in duplication carriers who have only one Xi, favors the scenario with genes on the Xp and Xq influencing methylation rather than an effect of an inactive X being present.
In contrast to the X chromosome, the presence of the Y-chromosome reduces the level of DNA methylation at the ZPBP2 DMR and makes it insensitive to the X chromosome dosage. The SRY region recapitulates the effect of the whole Y in our small sample of fibroblast cell lines. SRY is a transcription factor that plays a critical role in sex determination (reviewed in ). Moreover, the mouse Sry has been shown to contribute to sexual dimorphism in autosomal gene expression levels . It has been also implicated as the cause of sexual dimorphism in human brain development  and Hirschsprung disease . The mechanism of the SRY effect on methylation has yet to be discovered; however, it is worth noting that our reporter gene ZPBP2 is highly expressed in testis and implicated in male fertility . We therefore hypothesize that SRY binds the ZPBP2 promoter and shields it from DNA methylation, thereby preserving its activity in males, in the germ line, and in somatic cells. Such a mechanism would be beneficial for male fertility and explain the enrichment of spermatogenesis-associated genes among those with sex-specific promoter methylation [2, 4]. Thus, although SRY resides within a gene-rich region of the Y chromosome and, in our study, its effect cannot be separated from the effects of closely linked genes, it represents a good candidate whose role needs further investigation.
Interestingly, in Drosophila, sexual dimorphism in body size is regulated by the dosage of the X-linked gene Myc, which is not subject to dosage compensation, and the sex-determining gene transformer (tra) . Mouse studies dissecting the mechanism underlying sexual dimorphism in gene expression levels suggest contribution of X- and Y- encoded factors as well as sex hormones . Thus, the sum of current data suggests that the paradigm implicating interaction between X-chromosome-linked and sex-determining genes in sexual dimorphism in insects [10, 60, 61] and mice  may also explain some instances of sexually dimorphic autosomal DNA methylation in humans. We focused our study of sexual dimorphism in DNA methylation at a single autosomal locus. However, our findings may reflect a more general rule and larger genomic studies may find similar regulation of methylation at other autosomal loci.
Our data demonstrate that interaction between the Y chromosome and the X chromosome dosage defines DNA methylation levels at the autosomal locus ZPBP2. Furthermore, our findings suggest that at least one X-linked gene that influences ZPBP2 DMR methylation levels resides on the long arm of chromosome X. This is the first study that attempts dissecting the genetic mechanisms underlying sex-specific differences in methylation levels in a human autosomal region, and our findings may be applicable to other loci.
ATRX, chromatin remodeler
Differentially methylated region
JPX transcript, XIST activator
Peripheral blood cells
Principal component analysis
Ring finger protein, LIM domain interacting
Sex determining region Y
X-inactive specific transcript
zona pellucida binding protein 2
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.
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.
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.
Yousefi P, Huen K, Dave V, Barcellos L, Eskenazi B, Holland N. Sex differences in DNA methylation assessed by 450 K BeadChip in newborns. BMC Genomics. 2015;16:911.
Maschietto M, Bastos LC, Tahira AC, Bastos EP, Euclydes VL, Brentani A, Fink G, de Baumont A, Felipe-Silva A, Francisco RP, et al. Sex differences in DNA methylation of the cord blood are related to sex-bias psychiatric diseases. Sci Rep. 2017;7:44547.
Naumova AK, Al Tuwaijri A, Morin A, Vaillancout VT, Madore AM, Berlivet S, Kohan-Ghadr HR, Moussette S, Laprise C. Sex- and age-dependent DNA methylation at the 17q12-q21 locus associated with childhood asthma. Hum Genet. 2013;132(7):811–22.
Burgoyne PS, Arnold AP. A primer on the use of mouse models for identifying direct sex chromosome effects that cause sex differences in non-gonadal tissues. Biol Sex Differ. 2016;7:68.
Arnold AP, Cassis LA, Eghbali M, Reue K, Sandberg K. Sex hormones and sex chromosomes cause sex differences in the development of cardiovascular diseases. Arterioscler Thromb Vasc Biol. 2017;37(5):746–56.
Ghosh S, Klein RS. Sex drives dimorphic immune responses to viral infections. J Immunol. 2017;198(5):1782–90.
Mathews KW, Cavegn M, Zwicky M. Sexual dimorphism of body size is controlled by dosage of the X-chromosomal gene Myc and by the sex-determining gene tra in drosophila. Genetics. 2017;205(3):1215–28.
Ecker S, Chen L, Pancaldi V, Bagger FO, Fernandez JM, Carrillo de Santa Pau E, Juan D, Mann AL, Watt S, Casale FP, et al. Genome-wide analysis of differential transcriptional and epigenetic variability across human immune cell types. Genome Biol. 2017;18(1):18.
Warnefors M, Mossinger K, Halbert J, Studer T, VandeBerg JL, Lindgren I, Fallahshahroudi A, Jensen P, Kaessmann H. Sex-biased microRNA expression in mammals and birds reveals underlying regulatory mechanisms and a role in dosage compensation. Genome Res. 2017;27:1961.
Wijchers PJ, Yandim C, Panousopoulou E, Ahmad M, Harker N, Saveliev A, Burgoyne PS, Festenstein R. Sexual dimorphism in mammalian autosomal gene regulation is determined not only by Sry but by sex chromosome complement as well. Dev Cell. 2010;19(3):477–84.
Miller VM, Garovic VD, Kantarci K, Barnes JN, Jayachandran M, Mielke MM, Joyner MJ, Shuster LT, Rocca WA. Sex-specific risk of cardiovascular disease and cognitive decline: pregnancy and menopause. Biol Sex Differ. 2013;4(1):6.
Irizar H, Munoz-Culla M, Sepulveda L, Saenz-Cuesta M, Prada A, Castillo-Trivino T, Zamora-Lopez G, Lopez de Munain A, Olascoaga J, Otaegui D. Transcriptomic profile reveals gender-specific molecular mechanisms driving multiple sclerosis progression. PLoS One. 2014;9(2):e90482.
Jansen R, Batista S, Brooks AI, Tischfield JA, Willemsen G, van Grootheest G, Hottenga JJ, Milaneschi Y, Mbarek H, Madar V, et al. Sex differences in the human peripheral blood transcriptome. BMC Genomics. 2014;15:33.
van Nas A, Guhathakurta D, Wang SS, Yehya N, Horvath S, Zhang B, Ingram-Drake L, Chaudhuri G, Schadt EE, Drake TA, et al. Elucidating the role of gonadal hormones in sexually dimorphic gene coexpression networks. Endocrinology. 2009;150(3):1235–49.
Mamrut S, Avidan N, Staun-Ram E, Ginzburg E, Truffault F, Berrih-Aknin S, Miller A. Integrative analysis of methylome and transcriptome in human blood identifies extensive sex- and immune cell-specific differentially methylated regions. Epigenetics. 2015;10:943.
El-Maarri O, Becker T, Junen J, Manzoor SS, Diaz-Lacava A, Schwaab R, Wienker T, Oldenburg J. Gender specific differences in levels of DNA methylation at selected loci from human total blood: a tendency toward higher methylation levels in males. Hum Genet. 2007;122(5):505–14.
El-Maarri O, Walier M, Behne F, van Uum J, Singer H, Diaz-Lacava A, Nusgen N, Niemann B, Watzka M, Reinsberg J, et al. Methylation at global LINE-1 repeats in human blood are affected by gender but not by age or natural hormone cycles. PLoS One. 2011;6(1):e16252.
Vaissiere T, Hung RJ, Zaridze D, Moukeria A, Cuenin C, Fasolo V, Ferro G, Paliwal A, Hainaut P, Brennan P, et al. Quantitative analysis of DNA methylation profiles in lung cancer identifies aberrant DNA methylation of specific genes and its association with gender and cancer risk factors. Cancer Res. 2009;69(1):243–52.
Pilsner JR, Hall MN, Liu X, Ilievski V, Slavkovich V, Levy D, Factor-Litvak P, Yunus M, Rahman M, Graziano JH, et al. Influence of prenatal arsenic exposure and newborn sex on global methylation of cord blood DNA. PLoS One. 2012;7(5):e37147.
Kippler M, Engstrom K, Mlakar SJ, Bottai M, Ahmed S, Hossain MB, Raqib R, Vahter M, Broberg K. Sex-specific effects of early life cadmium exposure on DNA methylation and implications for birth weight. Epigenetics. 2013;8(5):494–503.
Sharma A, Jamil MA, Nuesgen N, Schreiner F, Priebe L, Hoffmann P, Herns S, Nothen MM, Frohlich H, Oldenburg J, et al. DNA methylation signature in peripheral blood reveals distinct characteristics of human X chromosome numerical aberrations. Clin Epigenetics. 2015;7:76.
Trolle C, Nielsen MM, Skakkebaek A, Lamy P, Vang S, Hedegaard J, Nordentoft I, Orntoft TF, Pedersen JS, Gravholt CH. Widespread DNA hypomethylation and differential gene expression in Turner syndrome. Sci Rep. 2016;6:34220.
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.
Yen A, Kellis M. Systematic chromatin state comparison of epigenomes associated with diverse properties including sex and tissue type. Nat Commun. 2015;6:7973.
Al Tuwaijri A, Gagne-Ouellet V, Madore AM, Laprise C, Naumova AK. Local genotype influences DNA methylation at two asthma-associated regions, 5q31 and 17q21, in a founder effect population. J Med Genet. 2016;53(4):232–41.
Martinez F, Rosello M, Mayo S, Monfort S, Oltra S, Orellana C. Duplication at Xq13.3-q21.1 with syndromic intellectual disability, a probable role for the ATRX gene. Am J Med Genet A. 2014;164A(4):918–23.
Linhares ND, Valadares ER, da Costa SS, Arantes RR, de Oliveira LR, Rosenberg C, Vianna-Morgante AM, Svartman M. Inherited Xq13.2-q21.31 duplication in a boy with recurrent seizures and pubertal gynecomastia: clinical, chromosomal and aCGH characterization. Meta Gene. 2016;9:185–90.
Giorgio E, Brussino A, Biamino E, Belligni EF, Bruselles A, Ciolfi A, Caputo V, Pizzi S, Calcia A, Di Gregorio E, et al. Exome sequencing in children of women with skewed X-inactivation identifies atypical cases and complex phenotypes. Eur J Paediatr Neurol. 2017;21(3):475–84.
Hu H, Haas SA, Chelly J, Van Esch H, Raynaud M, de Brouwer AP, Weinert S, Froyen G, Frints SG, Laumonnier F, et al. X-exome sequencing of 405 unresolved families identifies seven novel intellectual disability genes. Mol Psychiatry. 2016;21(1):133–48.
Coriell Institute for Medical Research. https://www.coriell.org/. Accessed 1 Nov 2017.
Berlivet S, Moussette S, Ouimet M, Verlaan DJ, Koka V, Al Tuwaijri A, Kwan T, Sinnett D, Pastinen T, Naumova AK. Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines. Hum Genet. 2012;131(7):1161–71.
Barndorff-Nilesen OE, Cox DR. Inference and asymptotics. London: Chapman & Hall; 1994.
Baumann C, Schmidtmann A, Muegge K, De La Fuente R. Association of ATRX with pericentric heterochromatin and the Y chromosome of neonatal mouse spermatogonia. BMC Mol Biol. 2008;9:29.
Sarma K, Cifuentes-Rojas C, Ergun A, Del Rosario A, Jeon Y, White F, Sadreyev R, Lee JT. ATRX directs binding of PRC2 to Xist RNA and Polycomb targets. Cell. 2014;159(4):869–83.
Watson LA, Goldberg H, Berube NG. Emerging roles of ATRX in cancer. Epigenomics. 2015;7(8):1365–78.
Fukuda A, Mitani A, Miyashita T, Sado T, Umezawa A, Akutsu H. Maintenance of Xist imprinting depends on chromatin condensation state and Rnf12 dosage in mice. PLoS Genet. 2016;12(10):e1006375.
Wang F, Shin J, Shea JM, Yu J, Boskovic A, Byron M, Zhu X, Shalek AK, Regev A, Lawrence JB, et al. Regulation of X-linked gene expression during early mouse development by Rlim. elife. 2016;5.
Balaton BP, Cotton AM, Brown CJ. Derivation of consensus inactivation status for X-linked genes from genome-wide studies. Biol Sex Differ. 2015;6:35.
Payer B. Epigenetic regulation of X-chromosome inactivation. In: Naumova AK, Taketo T, editors. Epigenetics in human reproduction and development. Singapore: World Scientific Publishing; 2016. p. 113–58.
Ronen D, Benvenisty N. Sex-dependent gene expression in human pluripotent stem cells. Cell Rep. 2014;8(4):923–32.
Bermejo-Alvarez P, Rizos D, Rath D, Lonergan P, Gutierrez-Adan A. Sex determines the expression level of one third of the actively expressed genes in bovine blastocysts. Proc Natl Acad Sci U S A. 2010;107(8):3394–9.
Gross N, Kropp J, Khatib H. Sexual dimorphism of miRNAs secreted by bovine in vitro-produced embryos. Front Genet. 2017;8:39.
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.
Grafodatskaya D, Chung BH, Butcher DT, Turinsky AL, Goodman SJ, Choufani S, Chen YA, Lou Y, Zhao C, Rajendram R, et al. Multilocus loss of DNA methylation in individuals with mutations in the histone H3 lysine 4 demethylase KDM5C. BMC Med Genet. 2013;6:1.
Cotton AM, Ge B, Light N, Adoue V, Pastinen T, Brown CJ. Analysis of expressed SNPs identifies variable extents of expression from the human inactive X chromosome. Genome Biol. 2013;14(11):R122.
Blewitt ME, Vickaryous NK, Hemley SJ, Ashe A, Bruxner TJ, Preis JI, Arkell R, Whitelaw E. An N-ethyl-N-nitrosourea screen for genes involved in variegation in the mouse. Proc Natl Acad Sci U S A. 2005;102(21):7629–34.
Juriloff DM, Harris MJ. Hypothesis: the female excess in cranial neural tube defects reflects an epigenetic drag of the inactivating X chromosome on the molecular mechanisms of neural fold elevation. Birth Defects Res A Clin Mol Teratol. 2012;94(10):849–55.
Robinson WH, Cotton AM, Penaherrera MS, Peeters SB, Brown CJ. X chromosome inactivation. In: Naumova AK, Greenwood CM, editors. Epigenetics and complex traits. New York: Springer Science + Business; 2013. p. 63–85.
Forejt J, Gregorova S, Goetz P. XY pair associates with the synaptonemal complex of autosomal male-sterile translocations in pachytene spermatocytes of the mouse (Mus musculus). Chromosoma. 1981;82(1):41–53.
Homolka D, Ivanek R, Capkova J, Jansa P, Forejt J. Chromosomal rearrangement interferes with meiotic X chromosome inactivation. Genome Res. 2007;17(10):1431–7.
Turner JM. Meiotic sex chromosome inactivation. Development. 2007;134(10):1823–31.
Fayer S, Yu Q, Kim J, Moussette S, Camerini-Otero RD, Naumova AK. Robertsonian translocations modify genomic distribution of gammaH2AFX and H3.3 in mouse germ cells. Mamm Genome. 2016;27(5–6):225–36.
Larney C, Bailey TL, Koopman P. Switching on sex: transcriptional regulation of the testis-determining gene Sry. Development. 2014;141(11):2195–205.
Sekido R. The potential role of SRY in epigenetic gene regulation during brain sexual differentiation in mammals. Adv Genet. 2014;86:135–65.
Li Y, Kido T, Garcia-Barcelo MM, Tam PK, Tabatabai ZL, Lau YF. SRY interference of normal regulation of the RET gene suggests a potential role of the Y-chromosome gene in sexual dimorphism in Hirschsprung disease. Hum Mol Genet. 2015;24(3):685–97.
Lin YN, Roy A, Yan W, Burns KH, Matzuk MM. Loss of zona pellucida binding proteins in the acrosomal matrix disrupts acrosome biogenesis and sperm morphogenesis. Mol Cell Biol. 2007;27(19):6794–805.
Verhulst EC, van de Zande L. Double nexus––Doublesex is the connecting element in sex determination. Brief Funct Genomics. 2015;14(6):396–406.
Ledon-Rettig CC, Zattara EE, Moczek AP. Asymmetric interactions between doublesex and tissue- and sex-specific target genes mediate sexual dimorphism in beetles. Nat Commun. 2017;8:14593.
We thank all participants. We are grateful to Teruko Taketo and Mary Fujiwara for critically reading the manuscript.
This work was supported by the operating grants from the Canadian Institutes for Health Research (CIHR) (to A.K.N. and C.L.). CL is a Canada Research Chair in Environment and Genetics of Respiratory Disorders and Allergies, director of the Asthma Strategic Group of the Respiratory Health Network (RHN) of FRSQ. CL builds, coordinates, and manages the SLSJ study. AA is a recipient of the National Guard Health Affairs Saudi Arabia scholarship.
Availability of data and materials
All data are presented in the paper and additional files.
Ethics approval and consent to participate
All experiments were conducted in accordance with the Canadian Tri-Council Policy Statement of Ethical Conduct for Research Involving Humans and approved by the Research Ethics Board of the McGill University Health Center. All participants gave their consents.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Experimental design. (DOCX 18 kb)
Table S1. Primers used in the study. (DOCX 12 kb)
Table S2. Fibroblast cell lines used for exploratory analysis. (DOCX 16 kb)
Table S3. Fibroblast cell lines used for mapping of the X-linked modifier of methylation. (DOCX 14 kb)
Figure S2. Sodium bisulfite sequencing methylation analysis of 51 CGs across the ZPBP2 promoter region shows higher methylation levels in fibroblast cell lines with two and three copies of the Xq arm. A: location of the 51 CGs interrogated using the sodium bisulfite sequencing assay shown in the context of the UCSC browser (hg19). B: heatmap representing mean methylation levels for each of the 51 CGs in cell lines with one X chromosome (two fibroblast cell lines with karyotype 45,X and no SRY region) and two or three copies of Xq (data from four fibroblast cell lines with karyotype 46,XX and two fibroblast cell lines from Turner syndrome patients with karyotype 46,i(Xq)). The black box beneath the heatmap shows the location of the 10 CGs analyzed by the pyrosequencing assay. The color scale for percent methylation is shown on the right. (DOCX 292 kb)
Table S4. Methylation percentages at the analyzed ZPBP2 DMR CGs from human non-transformed fibroblast cell lines, and the first PC of these values. (DOCX 14 kb)
Figure S3. Scatter plot of the first (PC1) and second (PC2) eigenvectors from a principal component analysis of methylation levels at 10 CpGs near the ZPBP2 transcriptional start site. The first principal component explains 59.2% of the variance in methylation levels. Points are colored by groups defined by presence of SRY and the number of X chromosomes. (DOCX 118 kb)
Table S5. Univariate analysis results, full model multivariate analysis results, as well as post hoc power calculations. (DOCX 14 kb)