Skip to main content

Advertisement

You are viewing the new article page. Let us know what you think. Return to old version

Research | Open | Published:

Sex influences DNA methylation and gene expression in human skeletal muscle myoblasts and myotubes

Abstract

Background

Sex differences are known to impact muscle phenotypes, metabolism, and disease risk. Skeletal muscle stem cells (satellite cells) are important for muscle repair and to maintain functional skeletal muscle. Here we studied, for the first time, effects of sex on DNA methylation and gene expression in primary human myoblasts (activated satellite cells) before and after differentiation into myotubes.

Method

We used an array-based approach to analyse genome-wide DNA methylation and gene expression in myoblasts and myotubes from 13 women and 13 men. The results were followed up with a reporter gene assay.

Results

Genome-wide DNA methylation and gene expression differences between the sexes were detected in both myoblasts and myotubes, on the autosomes as well as the X-chromosome, despite lack of exposure to sex hormones and other factors that differ between sexes. Pathway analysis revealed higher expression of oxidative phosphorylation and other metabolic pathways in myoblasts from women compared to men. Oxidative phosphorylation was also enriched among genes with higher expression in myotubes from women. Forty genes in myoblasts and 9 in myotubes had differences in both DNA methylation and gene expression between the sexes, including LAMP2 and SIRT1 in myoblasts and KDM6A in myotubes. Furthermore, increased DNA methylation of LAMP2 promoter had negative effects on reporter gene expression. Five genes (CREB5, RPS4X, SYAP1, XIST, and ZRSR2) showed differential DNA methylation and gene expression between the sexes in both myoblasts and myotubes. Interestingly, differences in DNA methylation and expression between women and men were also found during differentiation (myoblasts versus myotubes), e.g., in genes involved in energy metabolism. Interestingly, more DNA methylation changes occur in women compared to men on autosomes.

Conclusion

All together, we show that epigenetic and transcriptional differences exist in human myoblasts and myotubes as well as during differentiation between women and men. We believe that these intrinsic differences might contribute to sex dependent differences in muscular phenotypes.

Background

Sex contributes to differences in many aspects of metabolism and diseases such as obesity and diabetes [1], and it contributes to differences in skeletal muscle morphology and metabolism [2, 3]. Additionally, age-related muscle decline and remodelling differ in women versus men [4].

Adult skeletal muscle stem cells, so-called satellite cells, are responsible for regeneration and maintenance of skeletal muscle [5], thereby contributing to a healthy muscle phenotype. The satellite cells are activated in response to stress, after, e.g. injury or exercise, which initiates proliferation. Asymmetric cell division gives rise to new stem cells as well as muscle progenitors, called myoblasts, that eventually withdraw from the cell cycle, differentiate, fuse into myotubes and later mature into myofibres (Fig. 1) [6].

Fig. 1
figure1

Schematic figure of myogenesis as well as the study design. The leftmost part of the figure shows the adult myogenic process, from activation of skeletal muscle stem cells to a mature myofibre, as it occurs naturally in vivo. The rightmost part shows our study design covering the different analyses of genome-wide DNA methylation and gene expression in primary human myoblasts and in vitro differentiated myotubes from women and men. BMI, body mass index; Wilcoxon, Wilcoxon signed-rank test

Studies in mice have identified sex differences in skeletal muscle regeneration. Female myoblasts transplanted with higher efficiency regardless sex of the host, and transcriptional sex differences in the myoblasts were detected [7]. However, there is limited knowledge about the impact of sex on human myogenesis. Nevertheless, sex-specific differences in human skeletal muscle tissue transcriptome have been shown at baseline as well as in response to exercise and with age [8, 9]. Sex-specific muscle phenotypes may arise due to, e.g. hormonal effects of androgens or oestrogen, metabolic programming, or genetic differences [1]. Furthermore, sex-biased differentially expressed genes in skeletal muscle are located both on the autosomal chromosomes and the X-chromosome [8, 9].

Epigenetic modifications, such as DNA methylation, are important for inactivation of one X-chromosome in women, regulation of gene expression and differentiation of stem cells [10]. DNA methylations can be influenced by the environment and change with age [11]. We have previously shown genome-wide DNA methylation changes during differentiation of primary human myoblasts into myotubes, and obesity-associated epigenetic programming of these cells [12]. Sex-specific DNA methylation has been observed in other tissues, such as pancreatic islets, blood and liver among others [13,14,15,16,17]. However, it has not yet, to our knowledge, been studied in human muscle stem cells.

We hypothesise that sex influences DNA methylation and gene expression in human skeletal muscle stem cells and thereby affects myogenesis and muscle phenotypes. Our aim is to compare the genome-wide DNA methylome in primary human myoblasts and in vitro differentiated myotubes from women versus men, and relate epigenetic differences to gene transcription. We also studied the direct effect of increased methylation on transcriptional activity in cultured myoblasts.

Methods

Study subjects’ characteristics

Women and men were selected from a study population which has been previously described in detail elsewhere [12, 18, 19]. Thirteen women and 13 men without any known disease (except 6 individuals of each sex with obesity) were selected based on matching body mass index (BMI) and age from a subset of the cohort where DNA methylation and gene expression array data already existed (part of these data have been published previously [12] (Table 1). Data were not available for one man for the following phenotypes: glucose 2 h, Insulin 0 h, Insulin 2 h and homeostasis model assessment of insulin resistance and β-cell function (HOMA-IR and HOMA-B). Furthermore, mRNA expression data were not available for two women.

Table 1 Phenotypes of individuals included in the study

Human muscle stem cell isolation and culture

A detailed description of the method can be found in previous publications [12, 18]. In short, satellite cells were isolated from vastus lateralis. Myoblasts were seeded on BD Matrigel™ (BD Biosciences, NJ, USA), cultured in HAM/F10 supplied with 20% fetal bovine serum (FBS) and 1% penicillin/streptomycin (PS) and harvested when less than 50% confluent. At 70–80% confluency, growth medium was switched to DMEM (1.0 g/L glucose) supplied with 10% FBS and 1% PS until cells were confluent and aligned. Differentiation was then enhanced by switching to DMEM (4.5 g/L glucose) supplied with 2% horse serum and 1% PS for 5 days for myotube formation. Full differentiation was determined via visual confirmation where approximately 70% of the myotubes should have > 3 nuclei. DNA and RNA was extracted from the cells using the DNeasy blood and tissue kit (Qiagen, Hilding, Germany) and Trizol in combination with RNeasy minElute Cleanup kit (Qiagen), respectively.

Purity of the cultures in this cohort has been assessed by flow cytometry and published elsewhere [12]. All cultures expressed the myogenic marker CD56, while the endothelial and haematopoietic markers CD31 and CD45 were not detected.

DNA methylation analysis

DNA methylation was analysed genome-wide using the Infinium HumanMethylation450K BeadChip (Illumina, Inc., CA, USA). The method has been described in detail elsewhere [12]. Myoblasts and myotubes from the same individual were analysed on the same chip. Women and men were equally distributed on the different chips. Y-chromosome annotated probes were removed from initial analysis, together with rs-probes, non-CpG probes and cross-reactive probes with at least 49-base pair (bp) match [20]. Mean detection p value cut off was set to 0.01. Beta-mixture quantile normalisation (BMIQ) [21] was used to normalise between Infinium I and Infinium II assays. COMBAT [22] was used to correct for batch in non-paired analysis, i.e. women versus men. Methylation for each analysed site is presented as β value ranging from 0 to 1 (0–100% methylation). In total 458,475 probes were analysed for all chromosomes, 447,867 for autosomes and 10,608 for the X-chromosome. Y-chromosome DNA methylation data were normalised (including all chromosomes) and analysed separately.

The Illumina 450k array annotates the probes to different genomic regions [23]; in relation to the transcription start sites (TSS) either at 0 to 200 bp (TSS200) or between 200 and 1500 bp upstream of the TSS (TSS1500), to 5′ untranslated region (5′UTR), 1st exon, gene body (introns and exons, except 1st exon), 3′UTR or intergenic (not annotated to any gene/gene region). The probes are also annotated to CpG island regions; inside the CpG island itself, defined as 500 bp with GC content higher than 50% and observed/expected CpG ratio > 0.6 [23]. The 2-kb regions directly flanking the CpG island, either upstream (northern (N)) or downstream (southern (S)) of the islets, are called shores and the 2-kb regions flanking the shores are called the northern and southern shelf.

mRNA expression analysis

mRNA expression was analysed genome-wide using HumanHT-12 Expression BeadChip (Illumina). The array contains 47,231 probes and those with mean detection p value > 0.01 for more than 60% of the samples were filtered out. Data were background corrected, log2 transformed and quantile normalised using lumi package [24] in R (Additional file 1). In total, 16,955 transcripts, corresponding to 10,925 unique genes, were found expressed after quality control (QC). Of these genes, 9766 overlapped with genes from the methylation array.

Gene set enrichment analysis (GSEA)

GSEA was analysed using GSEA software [25, 26] to find enriched KEGG pathways. This software enabled us to analyse the complete expression data set using a pre-ranked list, based on t-statistics, where duplicate gene symbols were removed based on p value (lowest p value was kept in the data set). Default setting was used, except for minimal size of gene sets, which was set to two genes.

Luciferase assay

The luciferase assay method to measure the direct DNA methylation effect on gene expression has been described in detail elsewhere [12, 13]. In short, 1500-bp fragments of human promoters from LAMP2 and RPS4X were cloned into a CpG-free luciferase vector (pCpGL-basic) [27] and amplified by GenScript (GenScript USA Inc., Piscataway, NJ, USA). The constructs where then either mock methylated or methylated in vitro using the methyltransferases SssI, HhaI or HpaII (2.5 U/μg DNA) (New England Biolabs), which methylate cytosines in the following context: CG, GCGC and CCGG, respectively. Mouse myoblasts (C2C12) cultured in DMEM (4.5 g/L glucose) supplied with 10% FBS and 1% PS in a 96-well plate were co-transfected with either 150 ng DNA/well (LAMP2) or 50 ng DNA/well (RPS4X) of the pCpGL vector constructs together with a Renilla luciferase control reporter vector (Promega, Madison, WI, USA) using FuGeneHD (Promega). Luciferase and Renilla activity were measured 48 h post transfection in cell lysate using Dual-Luciferase Reporter Assay System (Promega). Luciferase activity was calculated as the ratio between the reporter gene firefly luciferase and the control vector Renilla luciferase.

DNA methylation and mRNA expression data in skeletal muscle

DNA methylation and gene expression data from non-diabetic individuals in a previously published study were used to study differences between sexes in skeletal muscle [28]. DNA methylation data were analysed by Infinium HumanMethylation450K BeadChip (Illumina) and gene expression data using GeneChip Human Gene 1.0 ST array by Applied Biosystems (Foster City, CA, USA). DNA methylation data for myotubes and skeletal muscle were merged based on Illumina TargetID and the expression data based on gene symbol.

Statistical analysis

t-tests were used to analyse differences in phenotypes between men and women and the data are presented as mean ± standard deviation (SD).

For DNA methylation and gene expression analyses, p values were calculated with linear regression corrected for age and BMI in the comparisons between women and men, and with a paired non-parametric test (Wilcoxon) for comparisons between myoblasts and myotubes from the same individuals. Frequencies of significant CpG sites in different genomic regions and regions in relation to CpG island were analysed with chi-square tests against the expected frequency (5% for p < 0.05) of all analysed sites. A false discovery rate (FDR) (Benjamini–Hochberg procedure) was performed to correct for multiple testing and to reduce the amount of false positives in all the array data (autosomal chromosomes and X-chromosome together) and GSEA. Separate FDR analyses were run for the Y-chromosome in men. FDR analyses of muscle data were performed after filtering on CpG sites/genes that overlapped with significant results in myotubes. Principal component (PC) analyses were performed after batch correction and the top PCs correlated with sex. Correlations between expression data from all samples were calculated using Pearson correlation analyses.

Luciferase assay results were analysed using paired t-test against the control.

Results

Differences in clinical phenotypes between women and men

Human myoblasts and myotubes derived from human satellite cells from 13 healthy women and 13 healthy men were included in the study (Fig. 1). Their characteristics are shown in Table 1. There were no differences in age or BMI between the groups. Women had higher gynoid fat mass and heart rate (HR), while men were taller and had higher whole-body fat-free mass, systolic blood pressure, waist-hip ratio and VO2 max (Table 1).

Sex-specific differences in DNA methylation and gene expression in human myoblasts

We started to analyse DNA methylation in cells harvested as proliferating myoblasts using Illumina 450k array (Fig. 1). Methylation data for a total number of 458,475 CpG sites, 447,867 on autosomal chromosomes and 10,608 on the X-chromosome were obtained from all individuals.

To examine potential sources of variation, we performed PC analysis of the methylation data set in myoblasts and correlated the top PCs with sex. The samples clustered based on sex (Additional file 2), as expected when methylation of the X-chromosome is included [13, 17]. Furthermore, sex was significantly correlated with the third PC (p = 1.47 × 10− 5), suggesting that sex has a relatively large impact on DNA methylation in myoblasts.

We continued to analyse the impact of sex on methylation of individual CpG sites in the myoblasts using linear regression. Based on FDR of less than 5% (q < 0.05), 12,177 CpG sites had significant differences, in methylation between women and men (Fig. 2a and Additional file 3). Of these, 5762 CpG sites were located on the autosomal chromosomes (q < 0.05) with absolute differences in methylation up to 46.1%. The majority (62%) had higher methylation in women compared to men. On the X-chromosome, 6415 CpG sites had different methylation (q < 0.05) between women and men with absolute differences up to 60.2%. Here, 78% had a higher methylation level in women, and notably 22% had a higher methylation level in men. More differences on the X-chromosome compared to autosomal chromosomes are expected due to the X-chromosome inactivation in women [13, 17, 29].

Fig. 2
figure2

Differential DNA methylation and gene expression in human myoblasts from women and men. a Manhattan plot of all analysed CpG sites for DNA methylation in myoblasts from women versus men. The blue line indicates the threshold for significance (q < 0.05). Green dots indicated with arrows are selected CpG sites for some genes of interest. Manhattan plot with black points are CpG sites only on the X-chromosome. N = 13 for both women and men. b, c A selection of relevant significant KEGG pathways (q < 0.05) from GSEA of expression data in myoblasts from women versus men (N = 11–13). Total, total number of genes in the pathway in the analysed dataset; Observed, observed number of genes in the dataset contributing to the enrichment; Val, valine; Leu, leucine; Ile, isoleucine; Gnrh, gonadotropin-releasing hormone; TGF-beta, transforming growth factor-beta. d Gene expression of selected genes in metabolic pathways and e DNA methylation of significant CpG sites annotated to these genes in myoblasts from women and men. *, q < 0.05, N = 11–13. Data are presented as mean ± SD

The CpG sites on the array have been annotated to different functional genomic regions as well as regions based on CpG content [23], see Additional file 4 for details. DNA methylation in these different regions might have diverse effects on transcription of their corresponding genes [10]. Hence, we evaluated where CpG sites with significant differences (q < 0.05) in DNA methylation between women and men are located. In general, we found more CpG sites than expected by chance close to TSS and in CpG islands on both autosomal chromosomes and the X-chromosome (Additional file 4).

To investigate if sex also affects gene expression in myoblasts, we used a microarray approach with mRNA from the same samples as for DNA methylation. Correlations of data from all samples were high (Additional file 1) and a PC analysis revealed significant correlation between sex and the third PC (p = 0.01) as well as the fourth PC (p = 8.64 × 10− 4).

We then performed GSEA [26] to analyse enrichment of KEGG pathways in the whole expression dataset. Twenty pathways were enriched (q < 0.05) with higher expression levels in women compared to men, including several pathways related to the cell cycle as well as metabolism of energy, proteins and fatty acids (Fig. 2b and Additional file 5). In contrast, only four pathways were enriched with higher expression in men compared to women (q < 0.05), mostly related to cell-cell communication, e.g. transforming growth factor-beta (TGF-beta) signalling (Fig. 2c and Additional file 5).

Next, we analysed expression of individual genes in myoblasts and found 137 unique genes with significantly different expression (q < 0.05) between women and men (Additional file 6). Importantly, 124 of these genes were located on autosomal chromosomes and include genes known to be involved in myogenesis and energy metabolism, in line with GSEA. For example, CTF1 had higher expression in men and MAMSTR, TFB1M and LDHB had higher expression in women.

DNA methylation is known to affect gene expression [10], and we therefore combined our gene expression data with the DNA methylation data. Forty unique genes had significant changes (q < 0.05) in both mRNA expression and methylation (Additional file 7), whereas 10 were located on the X-chromosome. Also here, we found several genes involved in metabolic pathways. For example, SIRT1 had higher expression in men, while ATP6V1B1, CPS1 and COQ3 had higher expression in women (Fig. 2d, e). Genes on the X-chromosome with both differential gene expression and DNA methylation included, e.g. XIST, involved in X-chromosome inactivation [30], LAMP2 and RPS4X (Fig. 3a–d).

Fig. 3
figure3

Genes with both differential DNA methylation and gene expression in human myoblasts from women versus men. ad Bar graphs showing mRNA expression and tables showing DNA methylation of a, b XIST and c, d LAMP2 and RPS4X in myoblasts from women versus men. Only data for CpG sites with significant differences between women and men are included. Below the tables are UCSC genome browser tracks for b XIST and d LAMP2 showing the location of the promoter region, CpG sites in the table above, the Refseq gene and histone modifications from ENCODE data in skeletal muscle myoblasts (HSMM). Red line on the chromosome above indicates location on the X-chromosome. N = 11–13, *q < 0.05; H3K4me1, Histone 3 lysine 4 (H3K4) methylation; H3K4me3, H3K4 trimethylation; H3K27ac, H3K27 acetylation; TSS200/TSS1500, 0–200 bp and 200–1500 bp upstream of the transcription start site (TSS), respectively. e In vitro methylation of LAMP2 and RPS4X promoters reduced expression (activity) of a luciferase reporter gene. Numbers above the bars indicate the number of cytosines targeted by each enzyme in the promoter sequences (SssI: all CG, HhaI: GCGC and HpaII: CCGG). *p < 0.05, N = 5–6. Mean of the control for each gene was set to 1. Data are presented as mean ± SD

DNA methylation of LAMP2 and RPS4X promoter regions directly affects transcriptional activity in myoblasts

DNA methylation in promoter regions is associated with low gene expression [10]. In order to study direct effects of promoter DNA methylation on gene expression level in myoblasts, we designed a luciferase assay by cloning the LAMP2 and RPS4X promoter regions, respectively, into a CpG-free vector. LAMP2 encodes lysosomal associated membrane protein 2 (LAMP-2) involved in autophagy, a process important for proper myogenesis [31,32,33]. LAMP2 had a higher expression level, accompanied with lower methylation level in the proximal promoter (TSS200), in myoblasts from men (Fig. 3c, d). RPS4X encodes a protein of the 40S ribosomal subunit and is known to escape X-chromosome inactivation [34]. As seen in Fig. 3c, d, RPS4X exhibited higher expression and lower promoter methylation in women. The luciferase assay clearly showed that increased methylation in the promoter region of these genes reduced expression of the reporter gene (Fig. 3e). This result supports that methylation in these promoters directly regulates the expression of their genes.

Sex-specific differences in DNA methylation and gene expression in human myotubes

We continued to analyse DNA methylation in myotubes derived from myoblasts of the 13 women and 13 men (Fig. 1). The samples clustered based on sex, as was also seen for the myoblasts, (Additional file 2), and sex was again significantly correlated with the third PC (p = 1.47 × 10–5).

Between women and men, 11,097 CpG sites had significantly different DNA methylation levels (q < 0.05) (Fig. 4a and Additional file 3). Of these CpG sites, 4918 were located on the autosomal chromosomes with absolute differences in methylation up to 39.2%, and 6179 sites were located on the X-chromosome, with absolute methylation differences up to 58.9%. Similar to the myoblasts, many significant sites on the X-chromosome (19%) had higher methylation level in men. We also found a similar pattern between myoblasts and myotubes when we analysed the distribution of these significant sites on the X-chromosome in relation to gene regions and CpG island regions (Additional file 4). However, these consistent patterns were not seen for the autosomal chromosomes (Additional file 4), indicating that sex-specific methylation differences on the X-chromosome are more persistent than on autosomal chromosomes during differentiation.

Fig. 4
figure4

Differential DNA methylation and gene expression in human myotubes from women and men. a Manhattan plot of all analysed CpG sites for DNA methylation between women and men in myotubes. The blue line indicates the threshold for significance (q < 0.05). Green dots indicated with arrows are selected CpG sites for some genes of interest. Manhattan plot with black points are CpG sites only on the X-chromosome. N = 13 for both women and men. b A selection of relevant significant KEGG-pathways (q < 0.05) from GSEA of expression data for women versus men in myotubes (N = 11–13). *, pathways enriched in both myoblasts and myotubes; Total, total number of genes in the pathway in the analysed dataset; Observed, observed number of genes in the dataset contributing to the enrichment

We proceeded to study the impact of sex on gene expression in myotubes. Correlations of data from all samples were high (Additional file 1), and sex correlated significantly with PC 3 to 5 (p = 0.003, 0.005 and 0.03) in a PC analysis.

GSEA of gene expression data in myotubes revealed 15 significant KEGG pathways (q < 0.05) that were upregulated in women compared to men (Fig. 4b and Additional file 5). Of note, 12 of these were also upregulated in myoblasts, including oxidative phosphorylation and DNA replication. However, no gene set was downregulated in myotubes from women compared to men.

Twenty-two unique genes were differently expressed (q < 0.05) between women and men in myotubes (Additional file 6). Of these, 13 genes were located on autosomal chromosomes and 9 genes on the X-chromosome.

Further, 9 unique genes in the myotubes had significant changes (q < 0.05) in both expression and DNA methylation between women and men (Additional file 7). Seven of these genes were located on the X-chromosome, including KDM6A (Fig. 4a) which is a histone demethylase important for myogenesis [35] and previously shown to be differentially expressed and methylated between women and men in other tissues [13, 17].

Sex differences in DNA methylation and gene expression in human myoblasts remain and arise in myotubes

To better understand which of the sex differences seen in myoblasts that persist after differentiation into myotubes, we compared our significant (q < 0.05) DNA methylation and gene expression hits in myoblasts and myotubes (Fig. 5a). Interestingly, more significant methylation and expression differences were seen in myoblasts compared to myotubes, both on autosomal chromosomes and on the X-chromosome. While 25% (1453) of the CpG sites that were differentially methylated in myoblasts on autosomal chromosomes remained significant in myotubes, 92% (5912 CpG sites) overlapped on the X-chromosome. Again, this suggests that sex differences in DNA methylation on the X-chromosome, compared to the autosomal chromosomes, are more persistent during myogenesis.

Fig. 5
figure5

Differences in expression and methylation between women and men in both myoblasts and myotubes. a Overlap of number of unique genes with significant differences (q < 0.05) in both gene expression and DNA methylation in myoblasts and myotubes from women versus men. b Gene expression and c DNA methylation in myoblasts and myotubes of CREB5, SYAP1, XIST, ZRSR2 and RPS4X. Only CpG sites with significant differences in DNA methylation in both myoblasts and myotubes are shown. d Genome browser track for ZRSR2 showing the location of the promoter region, CpG sites in the table above, the Refseq gene and histone modifications from ENCODE data in skeletal muscle myoblasts (HSMM). Red line on the chromosome above indicates location on the X-chromosome. N = 11–13, *q < 0.05; H3K4me1, Histone 3 lysine 4 (H3K4) methylation; H3K4me3, H3K4 trimethylation; H3K27ac, H3K27 acetylation; TSS200/TSS1500, 0–200 bp and 200–1500 bp upstream of the transcription start site (TSS), respectively. Data are presented as mean ± SD

Nine genes were differentially expressed (q < 0.05) between women and men in both myoblast and myotubes, and 5 of them (CREB5, SYAP1, XIST, ZRSR2 and RPS4X) also showed differences in DNA methylation (q < 0.05) (Fig. 5a–d). All of these genes, except CREB5, are located on the X-chromosome. They are involved in, e.g. splicing, X-chromosome inactivation and transcriptional activation [30, 36,37,38], suggesting that they may further regulate expression of other genes. Together, the data show that many sex differences in DNA methylation and gene expression, both on the X-chromosome and autosomal chromosomes, in human myoblasts remain after differentiation into myotubes. The findings also indicate that sex differences arise during cell specification, since several differences between women and men in myotubes were not seen in the myoblasts.

Sex-specific differences in global DNA methylation in human myoblasts and myotubes

Previous studies have shown sex differences in DNA methylation of specific functional genomic elements and CpG island regions [13, 17]. We therefore analysed the average DNA methylation levels in the myoblasts and myotubes from women and men for each of these regions (Additional file 8). Interestingly, men had higher methylation than women in 3′ untranslated regions (UTR) and shelves on the X-chromosome. Women had higher methylation in the other regions, as expected due to the role of DNA methylation for X-chromosome inactivation [13, 17, 29]. On the autosomal chromosomes, methylation levels in myotubes, but not myoblasts, were close to significantly higher in women compared to men in gene bodies, 3′ UTR, intergenic regions and CpG islands (q = 0.051). These data support that global methylation differences may arise during muscle differentiation.

Changes in DNA methylation and gene expression during myogenesis in women and men

We have previously shown that 3.7 times more DNA methylation changes occur during differentiation of myoblasts from obese compared with non-obese subjects [12]. We therefore studied whether sex also influences changes in DNA methylation and gene expression during myogenesis. Here, we analysed methylation and expression in myoblasts versus myotubes from women and men separately and compared the results (Fig. 6a, b and Additional file 3).

Fig. 6
figure6

Changes in DNA methylation and gene expression during myogenesis in women and men. a, b Overlap of number of a CpG sites with significantly different methylation (q < 0.05) and b unique genes with significantly different gene expression (q < 0.05) in myoblasts versus myotubes from women (N = 11) and men (N = 13). c mRNA expression of selected genes with expression difference in myoblasts versus myotubes from either women or men. In women, p = 0.58 and q = 0.69 for NDUFC1, p = 0.07 and q = 0.14 for NDUFB9. In men, p = 0.08 and q = 0.14 for MYOD1, p = 0.89 and q = 0.92 for LDHB. *q < 0.05. Data are presented as mean ± SD

In women, methylation status was altered in 24,723 CpG sites in the transition from myoblast to myotube (q < 0.05) (Fig. 6a). In men, methylation was altered in 16,789 sites in the same transition. Only 4258 of these sites were the same between women and men. Subsequently, most methylation changes seen during differentiation seem to be sex-specific. This was not driven by changes on the X-chromosome where relatively few changes occurred (212 in women and 262 in men, with 17 overlapping CpG sites). Interestingly, CpG sites that changed DNA methylation in women were enriched in gene bodies and open sea, while significant CpG sites in men were enriched close to TSS and in CpG islands (Additional file 4).

No CpG sites with significant difference in methylation were found on the Y chromosome in a separate analysis of DNA methylation in myoblasts versus myotubes from men (Additional file 9).

4333 unique genes in women and 5717 in men were differentially expressed between myoblasts and myotubes (q < 0.05) (Fig. 6b and Additional file 6). Expected expression changes of myogenic regulatory factors, muscle-specific genes and cell cycle genes during differentiation were seen for both women and men (Additional file 10) [12]. These data confirm that the myoblasts exit the cell cycle and differentiate into functional myotubes.

3742 genes with differential expression in myoblasts versus myotubes (q < 0.05) overlapped between women and men (Fig. 6b and Additional file 6). In line with the analysis in myoblasts and myotubes where we found differences in metabolic pathways between women and men, some genes involved in energy metabolism changed expression in either men or women during differentiation. For example, LDHB expression decreased in myoblasts versus myotubes from only women, while NDUFB9 and NDUFC1 expression increased in myoblasts versus myotubes from only men (Fig. 6c). Additionally, the expression of one myogenic regulatory factor, MYOD1, was only downregulated significantly in women (Fig. 6c).

Of note, three times more unique genes on the X-chromosome (68 versus 22) changed expression in men compared to women during myogenesis (Additional file 6).

Moreover, two genes on the Y-chromosome, KDM5D and EIF1AY, showed significant difference (q < 0.05) in gene expression between myoblasts and myotubes. Both had higher expression in myotubes (Additional file 9).

Sex-specific DNA methylation and gene expression in human skeletal muscle

Finally, we analysed if sex-specific differences in DNA methylation in myotubes can also be found in skeletal muscle tissue. Here, we used skeletal muscle from the non-diabetic individuals in a previously described cohort [28]. Methylation data were available from muscle of 7 women and 10 men. More than 50% of the differentially methylated CpG sites observed in myotubes (6539 out of 11,056) could be confirmed in muscle (q < 0.05) (Additional file 11). These included for example CpG sites annotated to KDM6A, RPS4X, SYAP1, XIST and ZRSR2 (Fig. 7a). A great majority of the overlapping CpG sites between myotubes and skeletal muscle showed higher methylation in women, and most were found on the X-chromosome.

Fig. 7
figure7

Differences in expression and methylation between women and men in and myotubes and skeletal muscle. a Table showing DNA methylation of significant CpG sites annotated to selected genes with differential DNA methylation and b bar graphs showing mRNA expression of genes with differential expression, between women and men in both myotubes and skeletal muscle. (N = 7–13), *q < 0.05. Data are presented as mean ± SD

Expression data in muscle from the same cohort as the methylation data were available (7 women and 12 men) for 22 of the genes with differential gene expression between women and men in myotubes (Additional file 11). Four of these (KDM6A, TXLNG, UBA1 and ZFP1) were also differentially expressed (q < 0.05) between women and men in whole muscle (Fig. 7b). ZFP1 is the only autosomal gene among these (chromosome 16), and the only one with higher expression in men compared to women.

Discussion

This study casts light on epigenetic sex differences present in human myoblasts and myotubes. Previous research has shown that sex differences are present from whole body metabolism, to tissues and on a cellular and genetic level [1,2,3, 13]. It is therefore important to bear in mind that sex might significantly influence the results of biomedical research. However, this factor is still often overlooked in different studies and in study designs.

Here, we found evidence that the sex origin of human myoblasts and myotubes influences their epigenome and transcriptome, which in turn can influence their differentiation and metabolism [12, 35, 39]. DNA methylation and gene expression differences between women and men that persist after myoblast differentiation into myotubes were found especially on the X-chromosome, but also on autosomal chromosomes. In addition, sex differences in DNA methylation emerged during differentiation at thousands of CpG sites, and were replicated in skeletal muscle tissue from another cohort.

Our study design, where cells from women and men were cultured under the same conditions, made it possible to study sex differences without acute effects from external factors, such as sex hormones, which might influence the results considerably [40,41,42]. Thus, the differences observed here are due to strictly genetic factors (XX or XY) or demonstrate an epigenetic memory from the in vivo environment. Indeed, obesity reprograms the DNA methylome of muscle stem cells [12], and elements such as sex hormones, exercise, diet and inflammatory components induce a long-term memory in the skeletal muscle [43,44,45,46]. Individuals in this study had an average BMI above normal weight (BMI > 25), which is also seen at a population level in many European countries [47]. Nevertheless, we adjusted for BMI in the statistical analyses, and PC analysis plots of methylation data did not cluster samples based on BMI or obesity. Methylation of only two CpG sites (cg02359186 and cg05801817) and expression of just one gene (MAMSTR) in myoblasts were associated with both sex and BMI. None of the sites/genes in myotubes significant for sex were also significantly associated with BMI. In addition, there were no difference in average BMI between women and men, and the number of obese individuals was the same in both groups. However, blood pressure and fat distribution differ between women and men in our study as expected [48,49,50,51,52]. These traits may affect whole body homeostasis and impact on the muscle stem cells [6]. We did not see a significant difference in VO2 max, an estimate of fitness levels, between the sexes after adjusting for weight. On the other hand, this has been observed previously and may require larger cohorts for detection [48, 53]. It is likely that a combination of genetic and environmental factors contributes to intrinsic sex differences in muscle stem cells. Another study has shown that knockout of the androgen receptor affects skeletal muscle of female and male mice differently [54]. Furthermore, better engraftment of muscle stem cells from female mice seems not to be caused by hormones or the immune response in the host, but innate factors in the stem cells [7]. This emphasises the importance of taking sex into account in biomedical research, as it can be of further relevance for future medicine related to muscular diseases and stem cell transplantation.

DNA methylation is known to differ between women and men, in particular due to the X-chromosome inactivation in women [13, 17, 29]. As expected, expression of XIST, a master regulator of X-chromosome inactivation, was significantly higher in both myoblasts and myotubes from women compared to men [13, 30]. Several CpG sites annotated to XIST also exhibited significant sex differences in methylation. This gene could serve as a positive control in studies regarding sex differences in gene expression and DNA methylation. Indeed, this gene has previously been suggested to serve as an expression marker, independent of tissue, for sex together with four Y-chromosome genes [55].

Epigenetic remodelling is important for cell specification and differentiation of muscle stem cells [12, 35, 39]. DNA methylation and histone modifications are known to cooperate during this process [56, 57]. We found lower methylation in KDM6A together with higher expression in myotubes and skeletal muscle from women versus men, in agreement with studies in other tissues [13, 17]. KDM6A is of great importance for removing repressive histone marks at muscle-specific genes during myogenesis [35]. SIRT1, a histone deacetylase, was also differentially methylated in myoblasts from women versus men and the expression was higher in men. A metabolic switch from fatty acid oxidation towards glycolysis during muscle stem cell activation leads to reduced NAD+ levels, which lowers SIRT1 activity [39]. This in turn leads to higher levels of H4K16 acetylation and expression of myogenic genes. Together, these indicate further differences in epigenetic regulation between women and men in myoblasts and myotubes. In future studies, it would be of interest to study whether histone modifications also differ between sexes in these cells, especially since other epigenetic enzymes changed expression in either men or women during differentiation.

The activity of epigenetic enzymes is closely linked to the metabolic state of the cell as indicated for SIRT1, and their substrate/products are mainly metabolic products [58]. GSEA showed higher expression of genes in oxidative phosphorylation in both myoblasts and myotubes from women. Oxidative metabolism enhances myoblast differentiation and muscle regeneration compared to glycolytic metabolism. In myotubes, high expression of genes in this pathway can indicate a better metabolic efficiency [39, 59, 60]. In line with this, expression of genes encoding parts of the respiratory chain (ATP6V1B1 and COQ3) and genes important for pyruvate generation (LDHB) were higher in myoblasts from women. Women also showed higher expression of genes in pathways related to metabolism of fatty acids and amino acids in myoblasts. These pathways provide metabolites to the oxidative metabolic pathways and are expected to be upregulated during differentiation [12]. Moreover, men had higher expression of TGF-beta signalling pathway in myoblasts. TGF-beta signalling promotes myoblast proliferation and inhibits muscle differentiation [61]. Treatment of human myoblasts with TGF-β1 reduces electron transport chain capacity and complex IV abundance [62]. It is tempting to speculate that metabolic regulation during myogenesis may be sex-specific, influence epigenetic mechanisms and thereby also the differentiation itself. In support of this hypothesis, we found differential DNA methylation and expression of, e.g. MAMSTR and CTF1 between women and men. MAMSTR showed higher expression in myoblasts from women while CTF1 showed higher expression in men. Mamstr acts as transcriptional regulator of MyoD in mice [63] and MYOD1 is a key myogenic regulatory factor that was differentially regulated during differentiation in women and men. Downregulation of MYOD1 in only women can be an indication of shorter differentiation processes in women. However, MYOD1 is normally not completely repressed in myotubes and has also been suggested to influence, e.g. fibre-type switch [64]. CTF1 encodes a cytokine, CT-1, which binds to the leukaemia inhibitory factor receptor (LIFR), necessary for myoblast proliferation [65]. CT-1 also plays a role in lipid homeostasis and improves glucose uptake [66]. Overall, these data support the notion that myoblasts from male mice proliferate better, while differentiation is enhanced in females [7].

Another cellular process, autophagy, is closely related to metabolism as well as crucial for myogenesis and maintenance of muscle homeostasis [33, 67]. However, the autophagic system needs to be fine-tuned. In the present study, we found that LAMP2, encoding a key regulator of lysosomal function and autophagy [68], had higher expression in men compared to women in the myoblasts. This was accompanied by lower methylation levels in the proximal promoter region in men. With the luciferase assay, we demonstrated that higher DNA methylation in the LAMP2 promoter directly reduces expression of the reporter gene. This implies that the methylation differences we found between women and men drive the difference observed in expression. Higher expression of LAMP2 in myoblasts from men can be an indicator of altered autophagy compared to women, which in turn might affect differentiation and myotube function [32, 33]. There is a clear evidence of autophagic vacuoles accumulating with LAMP-2 deficiency in muscle cells [69]. Thus, DNA methylation of the LAMP2 promoter can potentially be a way to help regulate autophagy. Of note, LAMP2 expression in monocytes from hypogonadal men is affected by testosterone replacement therapy [70].

Increased DNA methylation in the RPS4X promoter also affected transcriptional activity of the reporter gene directly. RPS4X escapes X-chromosome inactivation and higher expression levels in women are therefore expected [34]. Indeed, we found higher RPS4X expression in both myoblasts and myotubes from women, together with lower methylation levels in the promoter. Thus, genes on the X-chromosome that escape inactivation in women can have lower DNA methylation than in men, and this can directly affect gene expression. This pattern with regions of lower DNA methylation on the X-chromosome in women was also seen in the global analysis of DNA methylation.

One limitation in our study is the lack of adjustment for different numbers of X-chromosomes in women and men. However, we still believe that these data are of interest since it shows differences in DNA methylation between myoblasts and myotubes from women and men, respectively. Furthermore, we find regions with higher methylation in men on the X-chromosome, despite differences in the number of this chromosome, and several of these have been described in other metabolically important tissues before [13, 17]. DNA methylation is generally tissue-specific [10]. Differences between women and men in multiple tissues therefore suggest a common regulatory mechanism and point to an importance. For example, interplays of sex chromosomes with both autosomes and mitochondria have been demonstrated and may explain sex-specific phenotypes and risk of disease [71, 72].

Conclusions

We have shown that sex differences exist at the DNA methylome and transcriptome level in both myoblasts and myotubes, and these differences are likely to contribute to the phenotypical differences in myocytes and muscle tissue that exist between women and men.

Abbreviations

BMI:

Body mass index

BMIQ:

Beta-mixture quantile normalisation

bp:

Base pair

FBS:

Fetal bovine serum

FDR:

False discovery rate

Gnrh:

Gonadotropin-releasing hormone

GSEA:

Gene set enrichment analysis

H3K27ac:

Histone 3 acetylation on lysine 27

H3K4me1:

Histone 3 methylation on lysine 4

H3K4me3:

Histone 3 trimethylation on lysine 4

HOMA-B:

Homeostasis model assessment of β-cell function

HOMA-IR:

Homeostasis model assessment of insulin resistance

HSMM:

Skeletal muscle myoblasts

Ile:

Isoleucine

Leu:

Leucine

N shelf:

Northern shelf

N shore:

Northern shore

PC:

Principal component

PS:

Penicillin/streptomycin

QC:

Quality control

S shelf:

Southern shelf

S shore:

Southern shore

SD:

Standard deviation

TGF-beta:

Transforming growth factor-beta

TSS:

Transcription start site

Val:

Valine

References

  1. 1.

    Mauvais-Jarvis F. Sex differences in metabolic homeostasis, diabetes, and obesity. Biol Sex Differ. 2015;6:14.

  2. 2.

    Haizlip KM, Harrison BC, Leinwand LA. Sex-based differences in skeletal muscle kinetics and fiber-type composition. Physiology (Bethesda). 2015;30(1):30–9.

  3. 3.

    Lundsgaard AM, Kiens B. Gender differences in skeletal muscle substrate metabolism - molecular mechanisms and insulin sensitivity. Front Endocrinol (Lausanne). 2014;5:195.

  4. 4.

    Gheller BJ, Riddle ES, Lem MR, Thalacker-Mercer AE. Understanding age-related changes in skeletal muscle metabolism: differences between females and males. Annu Rev Nutr. 2016;36:129–56.

  5. 5.

    Saini A, Mastana S, Myers F, Lewis MP. ‘From death, lead me to immortality’ - mantra of ageing skeletal muscle. Curr Genomics. 2013;14(4):256–67.

  6. 6.

    Yin H, Price F, Rudnicki MA. Satellite cells and the muscle stem cell niche. Physiol Rev. 2013;93(1):23–67.

  7. 7.

    Deasy BM, Lu A, Tebbets JC, Feduska JM, Schugar RC, Pollett JB, et al. A role for cell sex in stem cell-mediated skeletal muscle regeneration: female cells have higher muscle regeneration efficiency. J Cell Biol. 2007;177(1):73–86.

  8. 8.

    Welle S, Tawil R, Thornton CA. Sex-related differences in gene expression in human skeletal muscle. PLoS One. 2008;3(1):e1385.

  9. 9.

    Lindholm ME, Huss M, Solnestam BW, Kjellqvist S, Lundeberg J, Sundberg CJ. The human skeletal muscle transcriptome: sex differences, alternative splicing, and tissue homogeneity assessed with RNA sequencing. FASEB J. 2014;28(10):4571–81.

  10. 10.

    Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.

  11. 11.

    Davegardh C, Garcia-Calzon S, Bacos K, Ling C. DNA methylation in the pathogenesis of type 2 diabetes in humans. Mol Metab. 2018;14:12–25.

  12. 12.

    Davegardh C, Broholm C, Perfilyev A, Henriksen T, Garcia-Calzon S, Peijs L, et al. Abnormal epigenetic changes during differentiation of human skeletal muscle stem cells from obese subjects. BMC Med. 2017;15(1):39.

  13. 13.

    Hall E, Volkov P, Dayeh T, Esguerra JL, Salo S, Eliasson L, et al. 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.

  14. 14.

    El-Maarri O, Becker T, Junen J, Manzoor SS, Diaz-Lacava A, Schwaab R, et al. 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.

  15. 15.

    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.

  16. 16.

    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.

  17. 17.

    Garcia-Calzon S, Perfilyev A, de Mello VD, Pihlajamaki J, Ling C. Sex differences in the Methylome and transcriptome of the human liver and circulating HDL-cholesterol levels. J Clin Endocrinol Metab. 2018;103(12):4395–408.

  18. 18.

    Green CJ, Pedersen M, Pedersen BK, Scheele C. Elevated NF-kappaB activation is conserved in human myocytes cultured from obese type 2 diabetic patients and attenuated by AMP-activated protein kinase. Diabetes. 2011;60(11):2810–9.

  19. 19.

    Moller JB, Pedersen M, Tanaka H, Ohsugi M, Overgaard RV, Lynge J, et al. Body composition is the Main determinant for the difference in type 2 diabetes pathophysiology between Japanese and Caucasians. Diabetes Care. 2014;37(3):796–804.

  20. 20.

    Chen YA, Lemire M, Choufani S, Butcher DT, Grafodatskaya D, Zanke BW, et al. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics. 2013;8(2):203–9.

  21. 21.

    Teschendorff AE, Marabita F, Lechner M, Bartlett T, Tegner J, Gomez-Cabrero D, et al. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics. 2013;29(2):189–96.

  22. 22.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

  23. 23.

    Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, et al. High density DNA methylation array with single CpG site resolution. Genomics. 2011;98(4):288–95.

  24. 24.

    Du P, Kibbe WA, Lin SM. lumi: a pipeline for processing Illumina microarray. Bioinformatics. 2008;24(13):1547–8.

  25. 25.

    Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, et al. PGC-1 alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.

  26. 26.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

  27. 27.

    Klug M, Rehli M. Functional analysis of promoter CpG methylation using a CpG-free luciferase reporter vector. Epigenetics. 2006;1(3):127–30.

  28. 28.

    Dayeh T, Tuomi T, Almgren P, Perfilyev A, Jansson PA, de Mello VD, et al. DNA methylation of loci within ABCG1 and PHOSPHO1 in blood DNA is associated with future type 2 diabetes risk. Epigenetics. 2016;11(7):482–8.

  29. 29.

    Hellman A, Chess A. Gene body-specific methylation on the active X chromosome. Science. 2007;315(5815):1141–3.

  30. 30.

    Cerase A, Pintacuda G, Tattermusch A, Avner P. Xist localization and function: new insights from multiple levels. Genome Biol. 2015;16:166.

  31. 31.

    Garcia-Prat L, Martinez-Vicente M, Perdiguero E, Ortet L, Rodriguez-Ubreva J, Rebollo E, et al. Autophagy maintains stemness by preventing senescence. Nature. 2016;529(7584):37–42.

  32. 32.

    McMillan EM, Quadrilatero J. Autophagy is required and protects against apoptosis during myoblast differentiation. Biochem J. 2014;462:267–77.

  33. 33.

    Fortini P, Ferretti C, Iorio E, Cagnin M, Garribba L, Pietraforte D, et al. The fine tuning of metabolism, autophagy and differentiation during in vitro myogenesis. Cell Death Dis. 2016;7:e2168.

  34. 34.

    Fisher EM, Beer-Romero P, Brown LG, Ridley A, McNeil JA, Lawrence JB, et al. Homologous ribosomal protein genes on the human X and Y chromosomes: escape from X inactivation and possible implications for Turner syndrome. Cell. 1990;63(6):1205–18.

  35. 35.

    Faralli H, Wang CC, Nakka K, Benyoucef A, Sebastian S, Zhuang LN, et al. UTX demethylase activity is required for satellite cell-mediated muscle regeneration. J Clin Investig. 2016;126(4):1555–65.

  36. 36.

    Nomura N, Zu YL, Maekawa T, Tabata S, Akiyama T, Ishii S. Isolation and characterization of a novel member of the gene family encoding the cAMP response element-binding protein CRE-BP1. J Biol Chem. 1993;268(6):4259–66.

  37. 37.

    Tronchere H, Wang J, Fu XD. A protein related to splicing factor U2AF35 that interacts with U2AF65 and SR proteins in splicing of pre-mRNA. Nature. 1997;388(6640):397–400.

  38. 38.

    Zinn AR, Alagappan RK, Brown LG, Wool I, Page DC. Structure and function of ribosomal protein S4 genes on the human and mouse sex chromosomes. Mol Cell Biol. 1994;14(4):2485–92.

  39. 39.

    Ryall JG, Dell'Orso S, Derfoul A, Juan A, Zare H, Feng X, et al. The NAD(+)-dependent SIRT1 deacetylase translates a metabolic switch into regulatory epigenetics in skeletal muscle stem cells. Cell Stem Cell. 2015;16(2):171–83.

  40. 40.

    Herbst KL, Bhasin S. Testosterone action on skeletal muscle. Curr Opin Clin Nutr Metab Care. 2004;7(3):271–7.

  41. 41.

    Sculthorpe N, Solomon AM, Sinanan ACM, Bouloux PMG, Grace F, Lewis MP. Androgens affect Myogenesis in vitro and increase local IGF-1 expression. Med Sci Sport Exer. 2012;44(4):610–5.

  42. 42.

    Kitajima Y, Ono Y. Estrogens maintain skeletal muscle and satellite cell functions. J Endocrinol. 2016;229(3):267–75.

  43. 43.

    Egner IM, Bruusgaard JC, Eftestol E, Gundersen K. A cellular memory mechanism aids overload hypertrophy in muscle long after an episodic exposure to anabolic steroids. J Physiol. 2013;591(24):6221–30.

  44. 44.

    Jacobsen SC, Brons C, Bork-Jensen J, Ribel-Madsen R, Yang B, Lara E, et al. Effects of short-term high-fat overfeeding on genome-wide DNA methylation in the skeletal muscle of healthy young men. Diabetologia. 2012;55(12):3341–9.

  45. 45.

    Sharples AP, Stewart CE, Seaborne RA. Does skeletal muscle have an ‘epi’-memory? The role of epigenetics in nutritional programming, metabolic disease, aging and exercise. Aging Cell. 2016;15:603–16.

  46. 46.

    Seaborne RA, Strauss J, Cocks M, Shepherd S, O'Brien TD, van Someren KA, et al. Human skeletal muscle possesses an epigenetic memory of hypertrophy. Sci Rep. 2018;8:1898.

  47. 47.

    Global Health Observatory (GHO) data - Overweight and obesity. WHO; (2017). http://www.who.int/gho/ncd/risk_factors/overweight/en/. Accessed 4 Oct 2018.

  48. 48.

    Wheatley CM, Snyder EM, Johnson BD, Olson TP. Sex differences in cardiovascular function during submaximal exercise in humans. Springerplus. 2014;3:445.

  49. 49.

    Xing CY, Tarumi T, Liu J, Zhang Y, Turner M, Riley J, et al. Distribution of cardiac output to the brain across the adult lifespan. J Cereb Blood Flow Metab. 2017;37(8):2848–56.

  50. 50.

    Wiklund P, Toss F, Weinehall L, Hallmans G, Franks PW, Nordstrom A, et al. Abdominal and gynoid fat mass are associated with cardiovascular risk factors in men and women. J Clin Endocrinol Metab. 2008;93(11):4360–6.

  51. 51.

    Gallagher D, Visser M, Wang Z, Harris T, Pierson RN Jr, Heymsfield SB. Metabolically active component of fat-free body mass: influences of age, adiposity, and gender. Metabolism. 1996;45(8):992–7.

  52. 52.

    Merz B, Nothlings U, Wahl S, Haftenberger M, Schienkiewitz A, Adamski J, et al. Specific metabolic markers are associated with future waist-gaining phenotype in women. PLoS One. 2016;11(6):e0157733.

  53. 53.

    Sharma HB, Kailashiya J. Gender difference in aerobic capacity and the contribution by body composition and Haemoglobin concentration: a study in young Indian National Hockey Players. J Clin Diagn Res. 2016;10(11):Cc9–Cc13.

  54. 54.

    MacLean HE, Chiu WSM, Notini AJ, Axell AM, Davey RA, McManus JF, et al. Impaired skeletal muscle development and function in male, but not female, genomic androgen receptor knockout mice. FASEB J. 2008;22(8):2676–89.

  55. 55.

    Staedtler F, Hartmann N, Letzkus M, Bongiovanni S, Scherer A, Marc P, et al. Robust and tissue-independent gender-specific transcript biomarkers. Biomarkers. 2013;18(5):436–45.

  56. 56.

    Conerly ML, Teves SS, Diolaiti D, Ulrich M, Eisenman RN, Henikoff S. Changes in H2A.Z occupancy and DNA methylation during B-cell lymphomagenesis. Genome Res. 2010;20(10):1383–90.

  57. 57.

    Ooi SK, Qiu C, Bernstein E, Li K, Jia D, Yang Z, et al. DNMT3L connects unmethylated lysine 4 of histone H3 to de novo methylation of DNA. Nature. 2007;448(7154):714–7.

  58. 58.

    Feil R, Fraga MF. Epigenetics and the environment: emerging patterns and implications. Nat Rev Genet. 2012;13(2):97–109.

  59. 59.

    Ryall JG. Metabolic reprogramming as a novel regulator of skeletal muscle development and regeneration. FEBS J. 2013;280(17):4004–13.

  60. 60.

    Cerletti M, Jang YC, Finley LW, Haigis MC, Wagers AJ. Short-term calorie restriction enhances skeletal muscle stem cell function. Cell Stem Cell. 2012;10(5):515–9.

  61. 61.

    Delaney K, Kasprzycka P, Ciemerych MA, Zimowska M. The role of TGF-beta 1 during skeletal muscle regeneration. Cell Biol Int. 2017;41(7):706–15.

  62. 62.

    Hoffmann C, Hockele S, Kappler L, de Angelis MH, Haring HU, Weigert C. The effect of differentiation and TGF beta on mitochondrial respiration and mitochondrial enzyme abundance in cultured primary human skeletal muscle cells. Sci Rep. 2018;8:737.

  63. 63.

    Mokalled MH, Johnson AN, Creemers EE, Olson EN. MASTR directs MyoD-dependent satellite cell differentiation during skeletal muscle regeneration. Genes Dev. 2012;26(2):190–202.

  64. 64.

    Ehlers ML, Celona B, Black BL. NFATc1 controls skeletal muscle fiber type and is a negative regulator of MyoD activity. Cell Rep. 2014;8(6):1639–48.

  65. 65.

    Broholm C, Laye MJ, Brandt C, Vadalasetty R, Pilegaard H, Pedersen BK, et al. LIF is a contraction-induced myokine stimulating human myocyte proliferation. J Appl Physiol. 2011;111(1):251–9.

  66. 66.

    Lopez-Yoldi M, Moreno-Aliaga MJ, Bustos M. Cardiotrophin-1: a multifaceted cytokine. Cytokine Growth Factor Rev. 2015;26(5):523–32.

  67. 67.

    Sandri M. Autophagy in skeletal muscle. FEBS Lett. 2010;584(7):1411–6.

  68. 68.

    Kaushik S, Cuervo AM. Chaperone-mediated autophagy: a unique way to enter the lysosome world. Trends Cell Biol. 2012;22(8):407–17.

  69. 69.

    Endo Y, Furuta A, Nishino I. Danon disease: a phenotypic expression of LAMP-2 deficiency. Acta Neuropathol. 2015;129(3):391–8.

  70. 70.

    Corrales JJ, Almeida M, Martin-Martin L, Miralles JM, Orfao A. Testosterone replacement therapy in hypogonadal men is associated with increased expression of LAMP-2 (CD107b) by circulating monocytes and dendritic cells. Clin Endocrinol. 2014;80(4):577–84.

  71. 71.

    Silkaitis K, Lemos B. Sex-biased chromatin and regulatory cross-talk between sex chromosomes, autosomes, and mitochondria. Biol Sex Differ. 2014;5:2.

  72. 72.

    Arnold AP, Chen X, Itoh Y. What a difference an X or Y makes: sex chromosomes, gene dose, and epigenetics in sexual differentiation. Handb Exp Pharmacol. 2012;214:67–88.

Download references

Acknowledgements

Thanks to SCIBLU genomics facility at Lund University for performing the DNA methylation array and mRNA expression array analysis. Thanks to Petr Volkov and Alexander Perfilyev for bioinformatics support.

Funding

This study was supported with grants from Diabetes Wellness Sweden (781-2015PG), Swedish Research Council, Region Skåne (ALF), Knut and Alice Wallenberg Foundation, Novo Nordisk Foundation, Swedish Foundation for Strategic Research for IRC15-0067 and Albert Påhlsson foundation. The Centre of Inflammation and Metabolism is supported by the Danish National Research Foundation (DNRF55). The Centre for Physical Activity Research is supported by Trygfonden. The Novo Nordisk Foundation Center for Basic Metabolic Research is supported by the Novo Nordisk Foundation to University of Copenhagen. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Author information

BKP designed the clinical study, previously published elsewhere [19] and MP carried out the clinical study including obtaining biopsies and collected phenotypes from the cohort. CB, TIH and CS isolated and cultured the cells. CB, CD and CL designed the study. CD and EHW performed bioinformatics analyses, analysed and interpreted data. CD and EHW wrote the manuscript. CD and EHW made all figures. CD, CL and EHW edited the manuscript. CS and BKP are responsible for the cohort, CL, CS and BKP contributed with resources and acquired financial support. All authors read and approved the final manuscript.

Correspondence to Cajsa Davegårdh or Charlotte Ling.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the local ethics committees (National Committee on Health Research Ethics (DNVK) KF 01-141/04) and conducted in accordance with the principles of the Helsinki Declaration. All study participants had provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

QC metrics and Pearson correlation of gene expression data in human myoblasts and myotubes from 11 women and 13 men. A-B Cumulative density plotted against log 2 intensity of expression data (A) before and (B) after normalisation. C-D Quantile normalisation of signal intensities (C) before and (D) after normalisation. E-F Person correlation between expression data of all samples in E) myoblasts and F) myotubes. (PDF 301 kb)

Additional file 2:

Clustering of samples based on PC analyses of methylation data in human myoblasts and myotubes from 13 women and 13 men. x-axis shows PC2 and y-axis shows PC1. (PDF 48 kb)

Additional file 3:

DNA methylation data from the Infinium HumanMethylation450K BeadChip (Illumina) from 13 women and 13 men. Different Excel sheets correspond to different comparison, i.e. myoblasts (women vs. men), myotubes (women vs. men), women (myoblasts vs. myotubes) and men (myoblasts vs. myotubes). (XLSX 1020 kb)

Additional file 4:

Distribution and frequency of CpG sites with significant difference in DNA methylation between women and men. A-D Frequencies of CpG sites that exhibit significant (q < 0.05) DNA methylation differences in women versus men in myoblasts and myotubes on the X-chromosome and the autosomal chromosomes in functional gene regions and CpG island regions. E-F) Frequencies of CpG sites that exhibit significant (q < 0.05) DNA methylation differences in myoblasts versus myotubes in women and men in functional gene regions and CpG island regions. Frequencies are compared to all analysed CpG sites using chi-square tests. *q < 0.05 in comparison to all analysed, ¤q < 0.05 for women versus men; TSS200/TSS1500: proximal promoter at 0–200 bp and 200–1500 bp, respectively, upstream of the transcription start site (TSS); UTR: untranslated region; CpG island: 200 bp (or more) stretch of DNA with GC content > 50% and CpG observed/expected ratio higher than 0.6; Shore: region of 2000 bp directly flanking the CpG island upstream (northern (N)) or southern (downstream (S)); Shelf: regions of 2000 bp flanking the island shores. (PDF 48 kb)

Additional file 5:

Results from GSEA of enriched KEGG pathway. Results from GSEA of KEGG pathways based on expression data for women versus men in myoblasts and myotubes, respectively. FDR q values were used to determine pathway enrichment (q < 0.05). (XLSX 13 kb)

Additional file 6:

mRNA expression data from the HumanHT-12 Expression BeadChip (Illumina) in myoblasts and myotubes from 11 women and 13 men. Different Excel sheets correspond to different comparisons, i.e. myoblasts (women vs. men), myotubes (women vs. men), women (myoblasts vs. myotubes) and men (myoblasts vs. myotubes). (XLSX 2430 kb)

Additional file 7:

Overlap between significant DNA methylation data and mRNA expression data in myoblasts and myotubes. DNA methylation and mRNA expression merged on annotated gene symbol for genes with significant differences in gene expression. Different Excel sheets correspond to different comparison, i.e. myoblasts (women vs men) and myotubes (women vs men), women (myoblasts vs. myotubes) and men (myoblasts vs. myotubes). (XLSX 2880 kb)

Additional file 8:

Average DNA methylation in myoblasts and myotubes from 13 women and 13 men. Average DNA methylation of CpG sites annotated to functional gene regions and CpG island regions on autosomal chromosomes and the X-chromosome, respectively. Data are presented as mean ± SD. *q < 0.05. (PDF 107 kb)

Additional file 9:

DNA methylation and mRNA expression data on the Y-chromosome in myoblasts versus myotubes from 13 men. Different Excel sheets correspond to DNA methylation data and mRNA expression data respectively. (XLSX 61 kb)

Additional file 10:

mRNA expression in myoblasts versus myotubes from 11 women and 13 men of A myogenic regulatory factors, B myocyte enhancer factor-2 genes, C-D muscle-specific genes and E-F cell cycle genes. Data are presented as mean ± SEM. *q < 0.05. (PDF 171 kb)

Additional file 11:

DNA methylation and mRNA expression data in skeletal muscle from women (N = 7) and men (N = 10–12). Different Excel sheets correspond to DNA methylation data and mRNA expression data respectively. (XLSX 1230 kb)

Rights and permissions

Open Access This 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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Gender
  • Muscle stem cells
  • Myoblasts
  • Myotubes
  • Skeletal muscle
  • Epigenetics
  • mRNA expression
  • Transcriptome
  • Epigenome
  • Myogenesis