- Research
- Open access
- Published:
X chromosome dosage and the genetic impact across human tissues
Genome Medicine volume 15, Article number: 21 (2023)
Abstract
Background
Sex chromosome aneuploidies (SCAs) give rise to a broad range of phenotypic traits and diseases. Previous studies based on peripheral blood samples have suggested the presence of ripple effects, caused by altered X chromosome number, affecting the methylome and transcriptome. Whether these alterations can be connected to disease-specific tissues, and thereby having clinical implication for the phenotype, remains to be elucidated.
Methods
We performed a comprehensive analysis of X chromosome number on the transcriptome and methylome in blood, fat, and muscle tissue from individuals with 45,X, 46,XX, 46,XY, and 47,XXY.
Results
X chromosome number affected the transcriptome and methylome globally across all chromosomes in a tissue-specific manner. Furthermore, 45,X and 47,XXY demonstrated a divergent pattern of gene expression and methylation, with overall gene downregulation and hypomethylation in 45,X and gene upregulation and hypermethylation in 47,XXY. In fat and muscle, a pronounced effect of sex was observed. We identified X chromosomal genes with an expression pattern different from what would be expected based on the number of X and Y chromosomes. Our data also indicate a regulatory function of Y chromosomal genes on X chromosomal genes.
Fourteen X chromosomal genes were downregulated in 45,X and upregulated in 47,XXY, respectively, in all three tissues (AKAP17A, CD99, DHRSX, EIF2S3, GTPBP6, JPX, KDM6A, PP2R3B, PUDP, SLC25A6, TSIX, XIST, ZBED1, ZFX). These genes may be central in the epigenetic and genomic regulation of sex chromosome aneuploidies.
Conclusion
We highlight a tissue-specific and complex effect of X chromosome number on the transcriptome and methylome, elucidating both shared and non-shared gene-regulatory mechanism between SCAs.
Background
In females, monosomy X results in Turner syndrome (TS; 45,X) and in males an additional X chromosome gives rise to Klinefelter syndrome (KS, 47,XXY) [1, 2]. Both syndromes share overlapping phenotypic traits such as infertility, hypergonadotropic hypgonadism, an increased risk of the metabolic syndrome, type II diabetes, osteoporosis, autoimmune disorders, and neuropsychiatric and neurodevelopmental disorders [3,4,5,6], all adding to an increased overall morbidity and mortality [3, 5, 7]. However, syndrome-specific traits are also prevalent in both TS and KS, such as cardiac anomalies in TS and venous thrombosis in KS [8, 9]. The overlapping phenotypic traits, in addition to individual groupwise stigmata, indicates that X chromosome number impacts central biological pathways involved in disease development. Elucidating how X chromosome number impacts the methylome, transcriptome and proteome, and regulates genomic pathways and networks, is essential in deciphering the molecular underpinnings of the phenotype in TS and KS as well as other sex chromosome aneuplodies (SCAs).
During the past decade, several studies have widened the knowledge about sex chromosome biology [10,11,12,13,14,15,16]. Collectively, they demonstrate a significant impact of especially X chromosome number on the human transciptome and methylome. These studies have primarily used peripheral blood or cells derived from blood for genomic analysis [10,11,12,13,14,15,16]. However, gene expression as well as DNA methylation is tissue-specific. Thus, blood may not always be a representative and relevant tissue to study to elucidate the molecular underpinnings of altered sex chromosome number ultimately leading to the phenotype seen in SCAs. Presently, a significant gap exists as to how sex chromosome number affects the transcriptome and methylome in specific tissues.
To this end, we collected blood samples, fat and muscle tissue biopsies from adult individuals with diverse karyotypes including 45,X (n = 36), 46,XX (n = 34), 46,XY (n = 16) and 47,XXY (n = 22) and performed gene expression analysis (RNA-seq) and DNA methylation analysis (850 k- Illumina Infinium assay), allowing comprehensive integrative analysis of the X chromosome number effect on both the transcriptome and methylome across three different target tissues (Fig. 1).
Methods
Design
The samples in this study were derived from four cross-sectional studies on patients with TS and KS and controls (NCT00624949, NCT00999310, NCT02526628, NCT01678261). The primary outcomes of the clinical trials have been published [10,11,12,13,14], and the circular RNA profile has been characterized in a subset of these samples [15]. The first participant was enrolled in 2012 and the last was enrolled in 2021. 47,XXY males (n = 22, median age (95%CI) 40.6 (37.0–44.3) were recruited through endocrine and clinical genetic hospital outpatient clinics as well as fertility clinics in Denmark, whereas 45,X women (n = 36, median age (95%CI) 46.1 (42.8–49.4) were recruited through the Danish National Society of Tuner Syndrome Contact Group and an endocrine outpatient clinic in Denmark. Healthy, age-matched men (n = 16, median age (95%CI): 40.2 (33.7–46.7) and women (n = 34, median age (95%CI) 43.9 (40.2–47.5) were recruited through advertisement to serve as controls. We obtained blood samples and abdominal subcutaneous fat tissue biopsies from all 22 47,XXY, and vastus lateralis muscle biopsies from 21 47,XXY. Of the 36 patients with TS, all had a blood sample drawn, and 10 had an abdominal subcutaneous fat tissue biopsy and a vastus lateralis muscle biopsy taken. Sixteen male controls contributed with a blood sample, a subcutaneous fat tissue biopsy, and a vastus lateralis muscle biopsy. Of the 34 female controls, all had a blood sample drawn, 11 had an abdominal subcutaneous fat tissue biopsy, and 10 had a vastus lateralis muscle biopsy taken. DNA methylation were not performed on four blood samples (45,X, n = 3; female control, n = 1), two fat tissue samples (45,X, n = 1; female control, n = 2), and in one muscle tissue sample (47,XXY, n = 1) due to technical reasons. RNA sequencing was not performed on 6 blood samples (45,X, n = 4, 47,XXY, n = 2), two fat samples (45,X, n = 1; male controls, n = 1), and 6 muscle tissue samples (45,X. n = 1; female controls, n = 5, male controls, n = 1) due to technical reasons (Fig. 1).
Sample collection
Blood samples for RNA extraction were drawn from the antecubital vein into PAXgene Blood RNA Tubes, placed 2 h at room temperature, and then stored overnight at − 21°, before storage at − 80°. Blood samples for DNA extraction were drawn into EDTA-treated tubes, and the peripheral blood samples were stored immediately at − 80 °C. Abdominal subcutaneous fat tissue biopsies were obtained from the lower abdomen by liposuction technique. The procedure was performed under sterile conditions and local anesthesia (10 ml lidocaine, 10 mg/ml). The biopsies were immediately cleaned and snap-frozen in liquid nitrogen and stored at − 80°. Muscle biopsies were taken from the vastus lateralis. The procedure was performed under sterile conditions and local anesthesia (10 ml lidocaine, 10 mg/ml). A small incision was made in the skin parallel to the orientation of the muscle fibers below. The muscle tissue was extracted using a 5-mm Bergström needle. The biopsies were immediately snap-frozen in liquid nitrogen and stored at − 80°.
RNA and DNA extraction
Blood
Blood total RNA was purified using the PaxGene blood Kit 262,174 (Qiagen) followed by on-chip electrophoresis on a Tapestation 4200 RNA Screen Tape System (Agilent) and by UV measurements on a Lunatic (Unchained Labs). Synthesis of directional RNA-seq libraries was conducted using the KAPA RNA HyperPrep with RiboErase Globin (HMR) (Roche) following the recommended procedure. Library preparation was automated on a Sciclone NGS (Caliper, Perkin Elmer) liquid handling robot. RNA-seq library quality was estimated by on-chip electrophoresis on a Tapestation 4200 D100 Screen Tape System (Agilent). The concentrations of the libraries were estimated using the Qubit dsDNA HS Assay (Thermo Fisher). Five hundred nanograms of total RNA was used as input.
Genomic DNA from the blood samples was extracted using the QIAmp® Mini Kit (Qiagen, Germany).
Fat tissue
The frozen biopsies were dissociated using a CovarisSP02, and DNA and RNA was extracted by the AllPrep DNA/RNA/PROTEIN Mini kit 8004 (Qiagen). Synthesis of directional RNA-seq libraries were conducted as described above using the KAPA RNA HyperPrep with RiboErase kit (HMR) (Roche).
Muscle tissue
The biopsies were sliced into 10-µm sections using a cryostat. RNA was extracted from all sections by the RNeasy Fibrouse tissue Mini kit (Qiagen). Synthesis of directional RNA-seq libraries were conducted as described above using the KAPA RNA HyperPrep with RiboErase kit (HMR) (Roche). DNA from muscle tissue was extracted using Gentra Puregene tissue kit catalog number 158667.
The isolated DNA from all tissues was quantified using a Qubit fluorimeter (manufacturer and model).
RNA sequencing and bioinformatics
The RNA-seq libraries were multiplexed paired-end sequenced on an Illumina NovaSeq 6000 (100 bp). Paired de-multiplexed fastq files were quality controlled using FastQC (Babraham Bioinformatics). Adaptor removal, in addition to trimming of low-quality ends, was then carried out using Trim Galore (Babraham Bioinformatics). Gene expression was quantified by quasi-mapping using Salmon [16]. A decoy-aware transcriptome index was built based on the hg38 transcriptome and selective alignment was run using the fastq pairs as input. Transcript abundancies were summarized to a gene level using the R package Tximeta [17]. Differential expression was then assessed using the R Bioconductor package DESeq2 using adjusted p-value < 0.05 (FDR < 0.05) [18]. All plots were made using ggplot2 [19]. The division of X chromosome genes into classes was based on Tukiainen et al. [20].
Enrichment analysis
To investigate possible biological functions of the RNA expression changes in our cohort, we performed overrepresentation analysis (ORA) using Clusterprofiler [21]. Differentially expressed genes were used as input. Disease associations were identified by Disease Ontology (DO) enrichment analysis, while Gene Ontology Biological Processes (GOBP) were used for functional enrichment [22, 23]. The top-enriched terms for each enrichment analysis were sorted according to p-value and presented as barplots.
Weighted gene co-expression network analysis
To identify co-expressed genes from the RNA-seq data and relate these to the number of sex chromosomes, weighted correlation network analysis (WGCNA, v1.70.3) was applied [24]. For each tissue, outlier samples were detected using hierarchical sample clustering, removing one 47,XXY blood sample. A signed co-expression network was constructed for each tissue, using a one-step approach, calculating adjacency choosing an appropriate soft thresholding power with approximate scale-free topology. Gene clustering was performed on the signed Topology Overlap Matrix by hierarchical clustering, identifying modules via the blockwiseModules function with a minModuleSize of 30 and a mergeCutHeight of 0.25. The module eigengenes were calculated via the moduleEigengenes function, and eigengene significance and corresponding p-value were obtained for each module-trait association in each tissue. Intramodular connectivity and gene significance for the number of X or Y chromosomes were extracted for each module of interest. For modules of interest, hub genes were identified using the chooseTopHubInEachModule function.
Methylation analysis (Infinium MethylationEPIC)
Infinium MethylationEPIC
1 μg of genomic DNA from all three tissues was bisulfide converted using the Illumina iScan (platform). The methylation level was measured using the Infinium MethylationEPIC DNA Analysis BeadChip (Array type) at Eurofins Genomics AS. Shortly, the DNA was bisulfide converted and then analyzed using Illumina’s Infinium assay that relies on direct hybridization of genomic targets to array-bound sequences. Single base extension was followed by florescence staining, signal amplification, and scanning on an iScan instrument. The raw intensity values obtained from the Infinium MethylationEPIC chip was read into R (v. 4.1.1) and further processed using the R package Minfi [25]. We filtered out cross-reactive probes and poor performance probes with a detection p-value < 0.01. Next, we applied preprocessFunnorm normalization, a between-array normalization method that removes variation by regressing out variability inferred by the control probes [26]. Cross-reactive probes were then removed, and the methylation values were calculated as M-values (logit [beta]) (Equation (I)).
Multidimensional scaling plots, using the top 1000 CpG sites as input, were evaluated to identify clusters of samples behaving differently than expected. Finally, the probes were annotated to the human genome version38 (https://github.com/achilleasNP/IlluminaHumanMethylationEPICanno.ilm10b5.hg38.git). For differential methylation analysis, the M-values were analyzed using LIMMA (~ 0 + karyotype) [27]. Differentially methylated positions were defined as adj. p-value < 0.05 and delta M > 1.0. DMRcate [28] were used to identify DNA regions that were differentially methylated, using a false discovery rate = 0.05, adjusted p-value < 0.05 and a |mean M-value|> 0.1. We defined a DMR as a minimum of 2 consecutive CpGs. If there were more than 1000 nucleotides between significant CpG, the sites were divided into separate DMRs. DNA methylation has been shown to play an important role in modulating gene expression. Therefore, we correlated changes in methylation and gene expression. For the correlation analysis, all significant DMRs (adj. P-value < 0.05) were correlated with gene expression, also including non-significant expressed genes, in order to evaluate trends in methylation status and gene expression. These were plotted in a scatter plot based on log2FC and ∆M-value. In the second part of the analysis, we used Spearman’s rank correlation to analyze the correlation between the specific methylation level (M-value) and RNA expression level (normalized counts) for shared genes. We considered the correlations to be interesting (either negative or positive) if p < 0.1.
Results
X chromosome number affects X chromosomal transcription signature across tissue types
First, we investigated if X chromosome number influenced X chromosomal gene expression in a tissue-specific way. Using unsupervised clustering, a clear clustering of 45,X, 46,XX, 46,XY, and 47,XXY was evident based on X chromosomal gene expression in fat and muscle (Fig. 2B, C), whereas in blood 45,X clustered with 46,XY and 47,XXY clustered with 46,XX (Fig. 2A). The observed pattern in fat and muscle indicated a possible effect of sex on X chromosomal gene expression. However, we cannot clearly separate the effect of sex from the presence of a Y chromosome, and the numbers of X chromosomes, as our cohort did not include samples from sex reversal conditions.
The impact of sex chromosome number on the X chromosomal expression pattern in blood (Additional file 3: Tables S1-S6) confirmed previous findings of predominantly gene downregulation in 45,X and gene upregulation in 47,XXY compared to controls with the same sex (Table 1) [14, 29, 30]. In fat and muscle (Additional file 3: Tables S7-S18), we observed the same pattern, both when comparing with same sex controls and opposite sex controls (Table 1).
X chromosomal genes are classified into three categories according to their X inactivation status: (1) escape genes (genes that escape X inactivation), (2) inactive genes (genes inactivated by X inactivation), and (3) genes that show variable inactivation [20]. The escape genes also include the genes located in the pseudoautosomal region 1 (PAR1) of the X chromosome. Among differentially expressed X chromosomal genes (xDEGs), we found that escape genes and PAR1 genes were overrepresented in all three tissues, as expected, and that the majority of xDEGs were located on Xp (Additional file 1: Fig. S1).
We identified several overlapping xDEGs between 45,X vs. 46,XX and 47,XXY vs. 46,XY (Table 1). 14 xDEGs overlapped in all three tissues (Table 2, Additional file: Fig. S2). Of these, three (JPX, TSIX, and XIST) are known to be involved in X chromosome inactivation, seven (AKAP17A, CD99, DHRSX, GTPBP6, PPP2R3B, SLC25A6, ZBED1) are located in PAR1, and four (EIF2S3, KDM6A, PUDP, ZFX) are annotated as escape genes. Notably, 6 of these 11 genes have widespread regulatory genomic functions (Table 2).
Most genes annotated as escape genes were differentially expressed in 45,X vs. 46,XX (34 out of 46) and 47,XXY vs. 46,XY (26 out of 46) in blood (Fig. 2D). As expected, these genes were downregulated in 45,X vs. 46,XX and upregulated in 47,XXY vs. 46,XY, except AP1S2 (Fig. 2D). The observed changes in gene expression agreed with previous studies showing that escape genes tend to be expressed from the inactive X chromosome (Xi) at a level that is between 10 and 95% the expression level seen from the active X chromosome (Xa) [31,32,33]. The same pattern was observed in fat and muscle tissue (Fig. 2D, Additional file 2: Text S1).
Of the genes annotated as inactive (blood, n = 254; fat, n = 315; muscle, n = 279), we identified 40 xDEGs in blood, six xDEGs in fat and eight xDEGs in muscle (differentially expressed in either 45,X vs. 46,XX, or 47,XXY vs. 46,XY, or both contrasts) (Fig. 2D). In blood, 18 of the 40 genes were downregulated and 22 upregulated in 45,X vs. 46,XX (Fig. 2D) with eight of these genes (ACSL4, ARHGAP6, ARMCX6, CLCN5, MOSPD2, RP2, SEPTIN6, ZNF275) showing the same expression pattern as demonstrated by Zhang et al. [34]. Only one gene, SEPTIN6, was found to be differentially expressed between 47,XXY vs. 46,XY in blood (Fig. 2D) [34].
To further examine the expression profile of X chromosomal genes in relation to X and Y chromosome number, we analyzed the change in expression level of PAR1 genes, escape genes, and inactive genes in the four karyotype groups (Fig. 2E, Additional file 1: Fig. S3). The mean change in expression of PAR1 genes between the karyotype groups in all three tissue types were as we would expect based on sex chromosome number and previous studies demonstrating a minor spread of inactivation into the PAR1 region of Xi in individuals with two X chromosomes, leading to higher expression of PAR1 genes among males compared with females [31]. Also, the mean changes in expression of escape genes between karyotype groups in all three tissues were as we would expect based on X chromosome number and evidence that escape genes are not fully expressed form the Xi [31, 32]. A tissue-specific-expression pattern was seen for inactive genes (Fig. 2E). In blood, the mean expression level of inactive genes in groups with more than one X chromosome (e.g., 46,XX and 47,XXY) were less than the mean expression in 45,X, indicating that the expression of genes on the Xa is to some extent also repressed. Notably, we found that the gene expression changes in 45,X vs. 46,XX and 45,X vs. 46,XY had the same level (Fig. 2E), indicating that Y chromosomal genes may affect the expression of inactive genes on the active X chromosome (Xa). Furthermore, the log2FC of 45,X vs. 47,XXY was higher than both 45,X vs. 46,XX and 45,X vs. 46,XY (Fig. 2E), indicating that as the number of sex chromosome increases, the expression of genes annotated as inactive on the Xa decreases—perhaps as a compensatory mechanism. In muscle tissue, we observed the same pattern as in blood for 45,X vs. 46,XX. However, when making comparisons where one of the groups had a Y chromosome (e.g., 46,XY or 47,XXY), we observed that the overall expression of inactive genes was increased in groups with only X chromosome(s) (e.g., 45,X or 46,XY), indicating that in muscle, Y chromosomal genes may impact the expression of inactive genes. In fat tissue, we observed the expected pattern of inactive genes, with no significant differences in mean expression between groups with different X chromosome numbers.
The effect of X chromosome number is genome-wide
Previous studies have documented a genome-wide effect of X chromosome number in blood [14, 29, 34, 35]. Our differential expression analysis of autosomal genes revealed the same genome-wide effect in blood as well as in fat and muscle (Table 1, Additional file 1: Figs. S4 and S5, Additional file 3: Tables S19-S36). Unsupervised clustering revealed a clear genotype-specific clustering in all three tissues based on autosomal genes (Fig. 3A–C). Four autosomal differentially expressed genes (aDEGs) in fat, 6 aDEGs in muscle, and 27 aDEGs in blood were shared between 45X vs. 46,XX and 47,XXY vs. 46,XY (Fig. 3D). Notably, all 37 aDEGs, except one, displayed an inverse expression pattern in 45,X vs. 46,XX compared to 47,XXY vs. 46,XY, being either upregulated in 45,X vs. 46,XX and downregulated in 47,XXY vs. 46,XY or vice versa (Fig. 3D). Only few aDEGs in 45,X vs. 46,XX overlapped across all three tissues (CCDC28A, OVCH1-AS1) and in 47,XXY vs. 46,XY only one overlapping aDEG was found (DDX43) (Fig. 3D, Additional file 1: Fig. S6).
Differentially expressed genes are associated with comorbidities in 45,X and 47,XXY and reveal distinct co-expression patterns
To gain insight into the biological function of the changes seen in the transcriptome of 45,X and 47,XXY, we performed enrichment analysis using all DEGs as input (padj < 0.05) (aDEGs and xDEGs). Enrichment was found for several terms related to known comorbidities in 45,X and 47,XXY (immune system, inflammatory diseases, skin diseases, neurodegenerative disorders, language disorders, abnormalities of dentition, hearing loss, congenital malformations, coagulation) with few overlaps between tissue types (Fig. 4, Additional file 2: Text S2).
To further study the transcriptional impact of SCAs, weighted gene correlation network analysis (WGCNA) was applied to identify co-expressed groups of genes (Additional file 1: Figs. S7-S10). Gene expression data from each tissue was used as input, firstly from all chromosomes and secondly restricted to the autosomes. We looked for strong correlations between each module of co-expressed genes and (1) karyotype and (2) the number of sex chromosomes (Additional file 1: Figs. S7 and S9). Here, the number of Y chromosomes represented sex (46,XY and 47,XXY vs. 45,X and 46,XX), whereas the number of X chromosomes represented the X number, independent of sex (46,XX and 47,XXY vs. 45,X and 46,XY). Within strongly correlated modules (highlighted in Additional file 1: Fig. S7), we further investigated the correlation between each gene and trait (Additional file 1: Fig. S8).
When including genes from both autosomes and sex chromosomes, one module from each tissue type was strongly associated with the number of X chromosomes (“royalblue” for blood, “skyblue3” for fat and “darkolivegreen” for muscle) (Additional file 1: Fig. S7A, B, C), indicating that the expression of these sets of genes were dependent on X chromosome number. As expected, these modules contained many X-linked genes, primarily from PAR1 or with Y-linked homologs. Common genes for the three tissue modules were XIST, TSIX, AKAP17A, KDM6A, PUDP, KDM5C, and ZFX (Additional file 1: Fig. S8). One module in blood was correlated to the number of Y chromosomes (“darkorange” module) (Additional file 1: Fig. S7A), whereas several modules in fat and muscle were both positively and negatively associated with the number of Y chromosomes (Additional file 1: Fig. S7B, C). Again, this showed a stronger and more complex influence of sex on gene expression in fat and muscle compared to blood. As expected, the modules that were positively correlated to the number of X or Y chromosomes were positively correlated to 47,XXY and negatively correlated to 45,X (Additional file 1: Fig. S7), in line with the sex chromosome number.
We then identified expression patterns solely from autosomes to investigate potential ripple effects of SCAs on autosomal gene expression. A single module was associated with the number of X chromosomes in blood and muscle (“greenyellow” for blood, “lightgreen” for muscle) (Additional file 1: Figs. S9A, C and S10). However, these correlations were less significant and no convincing associations to the X count were detected in fat (Additional file 1: Fig. S9B). Again, several modules in fat and muscle were associated with the presence of a Y chromosome (Additional file 1: Fig. S9B, C), showing that the sex influence on fat and muscle is stronger than in blood when only considering autosomal genes.
The impact of X chromosome number and sex on DNA methylation across tissue types
To investigate the impact of X chromosome number and sex on the DNA methylation profile, we performed unsupervised clustering based on all X chromosomal and autosomal CpG methylation sites. In blood, based on X chromosomal methylation sites in blood, a distinct clustering based on X chromosome number was revealed. 45,X and 46,XY clustered together and 46,XX and 47,XXY clustered together (Fig. 5A). In fat and muscle, a possible effect of sex and X chromosome number was observed (Fig. 5B, C). Based on autosomal methylation sites, a 45,X cluster and a 47,XXY cluster were observed, whereas 46,XX and 46,XY clustered closely together (Fig. 5D), in agreement with previous studies [14, 29]. When extending our analysis to fat and muscle tissue, we also observed a possible effect of sex and X chromosome number (Fig. 5E, F).
To further characterize these tissue differences, we investigated the distribution of methylation values in the four karyotype groups across tissues. On the X chromosome, 45,X had a higher density of hypomethylated sites and a lower density of hypermethylated sites in blood compared to 46,XY (Fig. 5G). To our surprise, the opposite was true for the distribution of methylation values in fat and muscle tissue. Here, we observed a lower density of hypomethylated and higher density of hypermethylated sites in 45,X compared to 46,XY (Fig. 5H). When comparing the distribution of methylation values on the X chromosome between 47,XXY and 46,XX, the same density of hypomethylated and hypermethylated sites were seen in blood (Fig. 5H). In muscle and fat, 47,XXY had a lower density of hypermethylated sites and a higher density of sites with an intermediary methylation level compared to 46,XX (Fig. 5H). On autosomes, the distribution of methylation values in all three tissues were similar for 45,X and 46,XX, showing a normal biphasic distribution (Fig. 5I). The same pattern was observed for 47,XXY and 46,XY, but with a tendency for 47,XXY to have a higher density of hypermethylated sites in muscle (Fig. 5J). These data illustrate that the impact of X chromosome number on the methylome is highly tissue-specific and differentially affected by sex.
The number of differentially methylated positions (DMPs) revealed a genome-wide effect of altered X chromosome number on methylation, across all three tissues (Table 3, Additional file 1: Fig. S11). Comparing 45,X vs. 46,XX, the majority of DMPs were hypomethylated across all tissues, while comparing 47,XXY vs. 46,XY the majority of DMPs were hypermethylated (Table 3), consistent with previous findings of an inverse methylation profile from peripheral blood [34]. The few autosomal DMPs (aDMPs) that overlapped between tissues were hypomethylated in 45,X and hypermethylated in 47,XXY, while the overlapping X chromosomal DMPs (xDMPs) were hypermethylated in 45,X and hypomethylated in 47,XXY (Additional file 1: Fig. S12). This inverse methylation pattern reflected the number of sex chromosomes as the comparison used both 46,XX and 46,XY as controls, thereby minimizing the effect of sex.
Analysis of DMPs in relation to genic location adds further evidence of an inverse methylation pattern between 45,X and 47,XXY
It is well established that CpG islands are usually methylated on the inactive X genes, while gene bodies are unmethylated. Furthermore, hypermethylation of X chromosomal gene bodies are correlated with gene expression [36,37,38]. We analyzed the identified DMPs in relation to genic locations. The DMPs were annotated according to their relation to gene, using eight annotation categories (1st Exon, 3′UTR, 5′ UTR, Body, Exon boundaries, TSS1500, TSS200, Unknown). The genic location of hypermethylated xDMPs in 45,X vs. 46,XY and hypomethylated xDMPs in 47,XXY vs. 46,XX were primarily localized in gene bodies in all three tissues (Additional file 1: Fig. S13A, B, C). Hypomethylated xDMPs in 45,X vs. 46,XY were primarily located at 5′UTR, TSS200, and TSS1500 in fat and muscle, but primarily in gene body in blood (Additional file 1: Fig. S13A, B, C).
Hypermethylated xDMPs in 47,XXY vs. 46,XX were also localized at 5′UTR and TS200 in fat, muscle, and blood. The distribution of aDMPs in relation to genic location was similar across all three tissues (Additional file 1: Fig. S13D, E). In conclusion, we demonstrate a methylation landscape in relation to genic location that presents as a mirrored image between 45,X and 47,XXY suggesting a common mechanism affecting methylation at specific CpG sites in both 45,X and 47,XXY. Furthermore, the methylation changes observed in muscle and fat may have larger functional effect because of their location within 5′UTR, TSS200, and TSS1500.
45,X and 47,XXY impacts the landscape of differentially methylated regions
Since changes in regional methylation have been shown to play a functional role in gene regulation [39], we investigated if 45,X and 47,XXY had any contrasting differentially methylated regions (DMRs) (Table 4).
Groups with unequal number of X chromosomes (47,XXY vs. 46,XY and 45,X vs. 46,XX) exhibited the largest difference in X chromosomal DMRs, while 47,XXY vs. 46,XX and 45,X vs. 46,XY had similar numbers of DMRs (Table 4). The number of autosomal DMRs in the same sex comparisons 47,XXY vs. 46,XY and 45,X vs. 46,XX were similar across all tissues (Table 4). This indicates that both X chromosome number and sex may impact regional methylation—with a more pronounced effect on methylation when lacking an X chromosome (45,X) compared with an additional X chromosome (47,XXY) (Tables 3 and 4).
DMPs with methylation patterns in opposite directions in 45,X and 47,XXY have a stronger impact on autosomal gene expression in 45,X
Epigenetic mechanisms including DNA methylation can regulate gene expression, mediating gene silencing, or activation [40]. To investigate if the inverse DNA methylation pattern between 45,X and 47,XXY was reflected in the transcriptome, we correlated aDMPs and xDMPs in both 45,X and 47,XXY to gene expression (adjp < 0.05, absolute delta M-value > 0.3) (Additional file 1: Fig. S14, Additional file 2: Text S3). Autosomal DNA methylation affected autosomal gene expression more severely in 45,X compared to 47,XXY, as a larger proportion of the autosomal genes annotated to the aDMPs were differentially expressed in 45,X vs. 46,XX compared with the 47,XXY vs. 46,XY contrast (Additional file 1: Fig. S14).
In blood, three autosomal genes (OR2C3, TRIM2, and HCG11) and one X chromosomal gene (KDM5C) were hypermethylated in 47,XXY and hypomethylated in 45,X and also differentially expressed in both 45,X vs 46,XX and 47,XXY vs 46,XY. Six autosomal genes (NOVA1, ANKRD55, DIP2C, SORBS2, SERINC5, SRGAP1) and one X chromosomal gene (ARSD) were differentially expressed and hypomethylated in 47,XXY and hypermethylated in 45,X (Additional file 1: Fig. S14). In muscle, only one autosomal gene, HENMT1, hypermethylated in 47,XXY and hypomethylated in 45,X, and one X chromosomal gene, EIF2S3, hypomethylated in 47,XXY and hypermethylated in 45,X was observed (Additional file 1: Fig. S14).
Correlation between DMRs and DEGs across tissues identified more X chromosomal DMR/DEGs than autosomal
We defined a DMR/DEG pair as a gene that was differentially expressed (adjp < 0.05) and associated with a significant differentially methylated region annotation (p-value < 0.05, absolute mean M-value > 0.1 and CpG > 1). Overall, the correlation between DMRs and DEGs was sparse across all comparisons (Table 5, Additional file 1: Fig. S15, Additional file 2: Text S4). However, one gene was of particular interest, KDM6A, identified as a xDMR/xDEG in both 45,X vs. 46,XX and 47,XXY vs. 46,XY in both blood, fat, and muscle. This DMR/DEG pair was hypomethylated and upregulated in 47,XXY, while hypermethylated and downregulated in 45,X.
Discussion
In this study, we present a comprehensive and integrative analysis of the transcriptome and methylome in three different phenotypically relevant tissues, across four karyotype groups; 45,X, 46,XX, 46,XY, and 47,XXY. To our knowledge, this is the first study analyzing the impact of X chromosome number on the transcriptome and methylome across multiple tissues from the same individuals. Our results significantly add to the emerging literature of the genetic and epigenetic architecture of SCAs.
Collectively, we demonstrate that the inverse genome-wide transcriptome profile in 45,X and 47,XXY, previously described in blood, can be extended to other tissue types of the same embryonic origin (e.g., mesoderm). However, our data also revealed that the effect of X chromosome number is tissue-specific and further indicate a possible effect of sex and presence or absence of the Y chromosome on the transcriptome in fat and muscle.
The higher number of DEGs between 45,X and 46,XX in blood is partly reflected by the larger number of blood samples. Using random down-sizing of the number of blood samples, before differential expression analysis, to reflect the number of fat and muscle samples resulted in 469 DEGs and 412 DEGs, respectively.
The similar expression patterns in the different tissue types indicate the existence of a universal X chromosome number dependent gene-regulatory mechanism in SCAs. Such a gene-regulatory mechanism may be initiated by escape genes having transcriptional regulatory functions and being able to start a cascade of regulatory pathways which in the end lead to the expression pattern illustrated. Evidence of a common gene-regulatory mechanism in SCAs comes from the finding of overlapping xDEGs between 45,X and 47,XXY across tissues displaying an inverse expression pattern. These genes—AKAP17A, CD99, DHRSX, EIF2S3, GTPBP6, JPX, KDM6A, PP2R3B, PUDP, SLC25A6, TSIX, XIST, ZBED1, ZFX—may be central in the epigenetic regulation of sex chromosome aneuploidies. Zhang et al. demonstrated that ZFX binding motifs were significantly enriched in autosomal promoter and enhancer regions of aDEGs in 45,X, but not in 47,XXY [34]. ZFX, encoding a transcriptional factor, has previously been proven differentially expressed in blood from 45,X and 47,XXY [14, 29, 34]. Thus, ZFX could be a central gene initiating the ripple effect of SCAs on the genome. Moreover, we also identified several shared aDEGs between 45,X and 47,XXY in blood, fat, and muscle, respectively. All, except one of these, displayed an inverse expression profile, further indicating a shared genomic regulatory mechanism of SCAs. The interplay identified between X chromosomes and autosomes was supported by the WGCNA analysis, showing that gene networks from the autosomes could be associated to with the number of X chromosome and that altered X chromosome number drives more subtle, tissue-dependent, changes to a larger set of autosomal genes.
As anticipated, most PAR1 genes, across all three tissues, were downregulated in 45,X, while upregulated in 47,XXY. This agrees with previous studies performed on blood cells [34, 35]. In addition, we also demonstrated that PAR1 genes display an expression difference between 45,X vs. 46,XX and 47,XXY vs. 46,XY that followed the expected pattern based on sex chromosome number and previous studies demonstrating a minor spread of inactivation into the PAR1 region of Xi [31]. Notably, in muscle, SHOX was upregulated in 45,X vs. 47,XXY, and a trend towards upregulation in 45,X vs. 46,XX and downregulation in 47,XXY vs. 46,XY was observed. SHOX is involved in growth and expression changes in SHOX drive differences in height [41, 42]. SHOX is normally expressed in skeletal muscle and is typically downregulated in females compared with males [20]. Since counts in all 4 karyotypes were low (< 200), this finding should be interpreted with caution, but could suggest that some PAR1 genes can exhibit tissue-specific expression pattern. However, it is also possible that temporal variation in SHOX expression may exist and that high levels of SHOX is only expressed in temporal windows during embryonic development and during rapid growth in childhood [43, 44].
In addition, we showed that some escape genes are not differentially expressed between 45,X vs. 46,XX and 47,XXY vs. 46,XY. This could indicate that regulatory mechanisms are at play in 45,X and 47,XXY to diminish the gene imbalance of an altered X chromosome number. Moreover, we showed that some escape genes and inactive genes were differentially expressed in 45,X but not in 47,XXY, or vice versa, indicating that some regulatory mechanism may be specific for the specific SCA subtypes. Previous findings on peripheral blood showed a similar expression pattern in 24 out of 28 escape genes in 45,X vs. 46,XX and of 27 out of 32 escape genes in 47,XXY vs. 46,XY [34]. Our data is also in agreement with a study performed on lymphoblastoid cell lines, demonstrating that some escape genes show a sublinear increase in expression in relation to increased X chromosome number in blood [35].
Based on our analyses, some genes annotated as inactivated were differentially expressed between 45,X and 46,XX, 47,XXY and 46,X or both, adding them to the pool of possible X chromosomal candidate genes involved in the phenotypic traits of 45,X and 47,XXY. Interestingly, our results indicate that not only genes annotated as inactive are inactivated on the Xi, but that these genes are also repressed on the Xa in individuals with more than one X chromosome, as previously reported [35].
In addition, we also found evidence that Y chromosomal genes may have the ability to repress the expression of genes on Xa and increase expression of genes from the Xi. However, these findings need further validation, before any conclusions can be drawn.
Collectively, our findings elucidate a much more complex relationship between X chromosomal gene expression and X chromosome number. These alterations may be caused by altered expression level of transcription factors and could be a compensatory mechanism to dampen the genomic impact of altered X chromosome number. However, these changes could also play an important role in the phenotype of SCAs.
Our analysis also revealed a tissue-specific and genome-wide impact of X chromosome number on the methylome. 45,X and 47,XXY demonstrated a unique methylation profile in all three tissues. Furthermore, our analysis indicated a possible impact of sex on the methylome in fat and muscle. We demonstrated that the general hypomethylation in 45,X and hypermethylation in 47,XXY in blood, described previously [14, 29, 34], could not be extended to fat and muscle tissue. In muscle tissue, 45,X was associated with autosomal hypermethylation compared to both 46,XX and 46,XY, while 47,XXY was associated with autosomal hypomethylation compared to both 46,XX and 46,XY. This pattern of autosomal methylation in 45,X and 47,XXY in muscle, is in agreement with the finding that 46,XY males present with hypomethylation of autosomal genes involved in muscle contraction and substrate metabolism compared with 46,XX females [45]. However, 47,XXY were associated with autosomal hypomethylation and 45,X with autosomal hypermethylated when compared to controls of the same sex. These altered autosomal methylation patterns in muscle compared to controls of the same sex may be linked to the increased risk of the metabolic syndrome and type 2 diabetes frequently seen in 45,X and 47,XXY [3, 46, 47].
Our finding that loss of an X chromosome had a more severe impact on DNA methylation compared with an additional X chromosome is in line with a previous study [30]. This may reflect a comparatively effective X inactivation of the additional X in 47,XXY and/or the lack of two copies of PAR genes in 45,X necessary for obtaining epigenetic balance [48]. Nevertheless, it is evident that common epigenetic regulatory mechanisms are at play in 45,X and 47,XXY, as 45,X and 47,XXY share many DMPs that display an inverse methylation pattern. This pattern has been demonstrated for blood [30, 34], but here we also show the existence of such a pattern in fat and muscle.
Previous studies have documented changes in the methylome and transcriptome in 45,X and 47,XXY to be complementary rather than overlapping [14, 34]. Here, we demonstrate several X chromosomal and autosomal genes that are both differentially methylated at single CpG sites and differentially expressed. We also demonstrate several X chromosomal and autosomal DMR/DEGs across all three tissues, although mostly prevalent in blood and mostly on the X chromosome. Some of these genes could warrant further attention, as they may be related to the phenotype of SCAs. CACNA2D1 is of interest as it was found to be hypomethylated and upregulated in 45,X in fat tissue. CACNA2D1 encodes a protein in the voltage-dependent calcium channel complex, which is essential in neurotransmission, and is associated with epilepsy [49], intellectual disability [50], and cardiac arrhythmias [51] which are conditions with an increased prevalence in 45,X [3, 12, 52]. STS and APLN, both hypomethylated and upregulated in 47,XXY in fat tissue, also seem to be of special interest. STS belongs to the sulfatase family that catalyzes the conversion of estrogen, cholesterol, and androgen precursors into their active form [53]. It is associated with increased blood pressure [54], which is a common comorbidity in sex chromosome abnormalities [46, 55]. APLN encodes a peptide hormone secreted by adipocytes, involved in the regulation of various metabolic functions and upregulated by insulin and obesity [56]. APLN expression has been positively correlated with total cholesterol and LDL cholesterol [57]. In muscle tissue, PWWP2B and THPO could be of special interest. PWWP2B was hypermethylated and downregulated in 45,X. Lack of PWWP2B promotes adipocyte thermogenesis [58]. THPO was hypermethylated and downregulated in 45,X. THPO encodes a cytokine that regulates platelet production and proliferation and controls thrombopoiesis [59]. Dysregulation of THPO has been implicated in cardiovascular diseases, which are prevalent in 45,X and 47,XXY [11, 60].
Of the observed DMR/DEG pairs, PRRG1 was of particular interest. PRRG1 encodes a vitamin K-dependent transmembrane protein and diseases associated with PRRG1 are involved in coagulation defects. Furthermore, a decrease in PRRG1 (as seen in 47,XXY) is associated with intima thickening in atherosclerosis [61]. SLC16A2 was hypomethylated and upregulated in blood in 45,X while hypermethylated and downregulated in 47,XXY in both blood and muscle tissue. SLC16A2 has been proposed involved in development of the central nervous system, while loss of function mutations has been associated with psychomotor retardation in males while females exhibit no neurological defects [62]. Mutations in SLC16A2 are the cause of Allan-Herndon-Dudley syndrome, associated with impaired thyroid metabolism, presenting with neurodevelopmental delay, hypertonia, and muscle weakness [63]. Both 47,XXY and 45,X suffer from thyroid disorders, and especially in 47,XXY neurocognition is affected, correlating well with downregulation of this gene [6]. MID2 was hypomethylated and upregulated in 45,X in both blood and muscle tissue. MID2 is associated with Opitz syndrome, a primary midline development disease characterized by congenital facial malformation, abnormalities of the central nervous system (including motor skill defects and developmental delay), and congenital heart defects [64, 65]. This gene was only affected in 45,X and could perhaps be involved in cardiac septal defects and coarctation of the aorta in 45,X [66].
KDM6A, has previously been found differentially expressed and methylated in blood from 45,X [14, 34]. Here, we demonstrate that KDM6A is a DMR/DEG pair in both 45,X and 47,XXY across all three tissues. KDM6A is thought to affect cardiac looping, neural tube development [67], germ cell development [68], and neurocognition and global growth and development of seizures [69]. MCF2 was hypermethylated and upregulated in 47,XXY and hypomethylated in 45,X across all three tissues. Differential methylation related to MCF2 is associated with several psychiatric disorders [70].
In muscle, one aDMR/DEG was of particular interest; PRKAG2, which was both hypermethylated and downregulated in 45,X. PRKAG2 encodes a subunit of the AMP activated protein kinase (AMPK) regulating muscle glucose metabolism and glycogen storage [71]. In response to reduced intracellular ATP levels, AMPK activates energy-producing pathways and inhibits biosynthesis of protein, carbohydrate, and lipids [72]. It is expressed in both skeletal muscle and cardiac muscle and has been shown to be involved in cardiac conduction disorders [73]. Furthermore, dysfunctional mutations in PRKAG2 have been proposed in liver cirrhosis [74]. Both liver disease and cardiac arrhythmias are prevalent in 45,X [3, 60].
Conclusion
The present study on the underlying genomics of 45,X and 47,XXY syndromes extends previous studies done on blood to other phenotypically relevant tissues and shows that there are both global and tissue-specific changes in both the methylome and transcriptome. We describe an inverse pattern of methylation and gene expression between 45,X and 47,XXY, a pattern most pronounced for X chromosomal genes. Furthermore, we point towards several new candidate genes. Our results add significantly to the emerging literature stating evidence of a comprehensive and complex X chromosome number effect on the human genome [34, 35].
Availability of data and materials
All sequence and methylation data from this study are deposited at the EGA repository (EGAS00001006996, EGAS00001007020) and will be available upon request within the General Data Protection Regulation.
Abbreviations
- Adjp:
-
Adjusted p-value
- SCAs:
-
Sex chromosome aneuploidies
- TS:
-
Turner syndrome
- KS:
-
Klinefelter syndrome
- DMR:
-
Differentially methylated region
- aDEG:
-
Autosomal differentially expressed genes
- xDEG:
-
X chromosomal differentially expressed genes
- aDMP:
-
Autosomal differentially methylated positions
- xDMP:
-
X chromosomal differentially methylated positions
- Fig:
-
Figure
- RNA seq:
-
Ribonucleic acid sequencing
- Xa:
-
Active X chromosome
- Xi:
-
Inactive X chromosome
- log2FC:
-
Log2fold change
References
Gravholt CH, Andersen NH, Conway GS, Dekkers OM, Geffner ME, Klein KO, et al. Clinical practice guidelines for the care of girls and women with Turner syndrome: proceedings from the 2016 Cincinnati International Turner Syndrome Meeting. Eur J Endocrinol. 2017;177(3):G1–70.
Zitzmann M, Aksglaede L, Corona G, Isidori AM, Juul A, T’Sjoen G, et al. European academy of andrology guidelines on Klinefelter Syndrome Endorsing Organization: European Society of Endocrinology. Andrology. 2021;9(1):145–67.
Viuff MH, Berglund A, Juul S, Andersen NH, Stochholm K, Gravholt CH. Sex hormone replacement therapy in Turner syndrome: impact on morbidity and mortality. J Clin Endocrinol Metab. 2020;105(2):468–78.
Jørgensen KT, Rostgaard K, Bache I, Biggar RJ, Nielsen NM, Tommerup N, et al. Autoimmune diseases in women with Turner’s syndrome. Arthritis Rheum. 2010;62(3):658–66.
Bojesen A, Gravholt CH. Morbidity and mortality in Klinefelter syndrome (47, XXY). Acta Paediatr. 2011;100(6):807–13.
Skakkebæk A, Moore PJ, Pedersen AD, Bojesen A, Kristensen MK, Fedder J, et al. The role of genes, intelligence, personality, and social engagement in cognitive performance in Klinefelter syndrome. Brain Behav. 2017;7(3):e00645.
Schoemaker MJ, Swerdlow AJ, Higgins CD, Wright AF, Jacobs PA. United Kingdom Clinical Cytogenetics G. Mortality in women with turner syndrome in Great Britain: a national cohort study. J Clin Endocrinol Metab. 2008;93(12):4735–42.
Zöller B, Ji J, Sundquist J, Sundquist K. High Risk of Venous Thromboembolism in Klinefelter Syndrome. J Am Heart Assoc. 2016;5(5):e003567.
Viuff MH, Trolle C, Wen J, Jensen JM, Norgaard BL, Gutmark EJ, et al. Coronary artery anomalies in Turner syndrome. J Cardiovasc Comput Tomogr. 2016;10(6):480–4.
Skakkebæk A, Gravholt CH, Rasmussen PM, Bojesen A, Jensen JS, Fedder J, et al. Neuroanatomical correlates of Klinefelter syndrome studied in relation to the neuropsychological profile. Neuroimage Clin. 2014;4:1–9.
Chang S, Biltoft D, Skakkebæk A, Fedder J, Bojesen A, Bor MV, et al. Testosterone treatment and association with thrombin generation and coagulation inhibition in Klinefelter syndrome: A cross-sectional study. Thromb Res. 2019;182:175–81.
Trolle C, Mortensen KH, Pedersen LN, Berglund A, Jensen HK, Andersen NH, et al. Long QT interval in Turner syndrome–a high prevalence of LQTS gene mutations. PLoS ONE. 2013;8(7):e69614.
Trolle C, Mortensen KH, Bjerre M, Hougaard DM, Cohen A, Andersen NH, et al. Osteoprotegerin in Turner syndrome - relationship to aortic diameter. Clin Endocrinol (Oxf). 2015;82(3):397–403.
Trolle C, Nielsen MM, Skakkebaek A, Lamy P, Vang S, Hedegaard J, et al. Widespread DNA hypomethylation and differential gene expression in Turner syndrome. Sci Rep. 2016;6:34220.
Johannsen EB, Just J, Viuff MH, Okholm TLH, Pedersen SB, Meyer Lauritsen K, et al. Sex chromosome aneuploidies give rise to changes in the circular RNA profile: A circular transcriptome-wide study of Turner and Klinefelter syndrome across different tissues. Front Genet. 2022;13:928874.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.
Love MI, Soneson C, Hickey PF, Johnson LK, Pierce NT, Shepherd L, et al. Tximeta: Reference sequence checksums for provenance identification in RNA-seq. PLoS Comput Biol. 2020;16(2):e1007664.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer; 2009. p. 213 (VIII).
Tukiainen T, Villani AC, Yen A, Rivas MA, Marshall JL, Satija R, et al. Landscape of X chromosome inactivation across human tissues. Nature. 2017;550(7675):244–8.
Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation. 2021;2(3):100141.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
The Gene Ontology Consortium. The Gene Ontology resource: enriching a GOld mine. Nucleic Acids Res. 2021;49(D1):D325-34.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9(1):559.
Fortin J-P, Triche TJ Jr, Hansen KD. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics. 2016;33(4):558–60.
Fortin J-P, Labbe A, Lemire M, Zanke BW, Hudson TJ, Fertig EJ, et al. Functional normalization of 450k methylation array data improves replication in large cancer studies. Genome Biol. 2014;15(11):503.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Peters TJ, Buckley MJ, Statham AL, Pidsley R, Samaras K, Lord RV, et al. De novo identification of differentially methylated regions in the human genome. Epigenetics Chromatin. 2015;8(1):6.
Skakkebæk A, Nielsen MM, Trolle C, Vang S, Hornshøj H, Hedegaard J, et al. DNA hypermethylation and differential gene expression associated with Klinefelter syndrome. Sci Rep. 2018;8(1):13740.
Sharma A, Jamil MA, Nuesgen N, Schreiner F, Priebe L, Hoffmann P, et al. DNA methylation signature in peripheral blood reveals distinct characteristics of human X chromosome numerical aberrations. Clin Epigenetics. 2015;7(1):76.
Balaton BP, Cotton AM, Brown CJ. Derivation of consensus inactivation status for X-linked genes from genome-wide studies. Biol Sex Differ. 2015;6(1):35.
Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature. 2005;434(7031):400–4.
Balaton BP, Brown CJ. Contribution of genetic and epigenetic changes to escape from X-chromosome inactivation. Epigenetics Chromatin. 2021;14(1):30.
Zhang X, Hong D, Ma S, Ward T, Ho M, Pattni R, et al. Integrated functional genomic analyses of Klinefelter and Turner syndromes reveal global network effects of altered X chromosome dosage. Proc Natl Acad Sci. 2020;117(9):4864–73.
Raznahan A, Parikshak NN, Chandran V, Blumenthal JD, Clasen LS, Alexander-Bloch AF, et al. Sex-chromosome dosage effects on gene expression in humans. Proc Natl Acad Sci USA. 2018;115(28):7398–403.
Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK. On the presence and role of human gene-body DNA methylation. Oncotarget. 2012;3(4):462–74.
Hellman A, Chess A. Gene body-specific methylation on the active X chromosome. Science. 2007;315(5815):1141–3.
Cotton AM, Lam L, Affleck JG, Wilson IM, Peñaherrera MS, McFadden DE, et al. Chromosome-wide DNA methylation analysis predicts human tissue-specific X inactivation. Hum Genet. 2011;130(2):187–201.
Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Gräf S, et al. An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res. 2008;18(9):1518–29.
Jaenisch R, Bird A. Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003;33(3):245–54.
Rao E, Weiss B, Fukami M, Rump A, Niesler B, Mertz A, et al. Pseudoautosomal deletions encompassing a novel homeobox gene cause growth failure in idiopathic short stature and Turner syndrome. Nat Genet. 1997;16(1):54–63.
Ottesen AM, Aksglaede L, Garn I, Tartaglia N, Tassone F, Gravholt CH, et al. Increased number of sex chromosomes affects height in a nonlinear fashion: a study of 305 patients with sex chromosome aneuploidy. Am J Med Genet A. 2010;152A(5):1206–12.
Marchini A, Rappold G, Schneider KU. SHOX at a glance: from gene to protein. Arch Physiol Biochem. 2007;113(3):116–23.
Clement-Jones M, Schiller S, Rao E, Blaschke RJ, Zuniga A, Zeller R, et al. The short stature homeobox gene SHOX is involved in skeletal abnormalities in Turner syndrome. Hum Mol Genet. 2000;9(5):695–702.
Landen S, Jacques M, Hiam D, Alvarez-Romero J, Harvey NR, Haupt LM, et al. Skeletal muscle methylome and transcriptome integration reveals profound sex differences related to muscle function and substrate metabolism. Clin Epigenetics. 2021;13(1):202.
Bojesen A, Høst C, Gravholt CH. Klinefelter’s syndrome, type 2 diabetes and the metabolic syndrome: the impact of body composition. Mol Hum Reprod. 2010;16(6):396–401.
Gravholt CH, Hjerrild BE, Mosekilde L, Hansen TK, Rasmussen LM, Frystyk J, et al. Body composition is distinctly altered in Turner syndrome: relations to glucose metabolism, circulating adipokines, and endothelial adhesion molecules. Eur J Endocrinol. 2006;155(4):583–92.
De Bonis ML, Cerase A, Matarazzo MR, Ferraro M, Strazzullo M, Hansen RS, et al. Maintenance of X- and Y-inactivation of the pseudoautosomal (PAR2) gene SPRY3 is independent from DNA methylation and associated to multiple layers of epigenetic modifications. Hum Mol Genet. 2006;15(7):1123–32.
Vergult S, Dheedene A, Meurs A, Faes F, Isidor B, Janssens S, et al. Genomic aberrations of the CACNA2D1 gene in three patients with epilepsy and intellectual disability. Eur J Human Genet EJHG. 2015;23(5):628–32.
Kessi M, Chen B, Peng J, Yan F, Yang L, Yin F. Calcium channelopathies and intellectual disability: a systematic review. Orphanet J Rare Dis. 2021;16(1):219.
Burashnikov E, Pfeiffer R, Barajas-Martinez H, Delpón E, Hu D, Desai M, et al. Mutations in the cardiac L-type calcium channel associated with inherited J-wave syndromes and sudden cardiac death. Heart Rhythm. 2010;7(12):1872–82.
Viuff MH, Stochholm K, Juul S, Gravholt CH. Disorders of the eye, ear, skin, and nervous system in women with Turner syndrome -a nationwide cohort study. Eur J Hum Genet. 2022;30:228-36.
Hobkirk R. Steroid sulfotransferases and steroid sulfate sulfatases: characteristics and biological roles. Can J Biochem Cell Biol. 1985;63(11):1127–44.
Valigora SD, Lib PK, Dunphy G, Turner M, Ely DL. Steroid sulfatase inhibitor alters blood pressure and steroid profiles in hypertensive rats. J Steroid Biochem Mol Biol. 2000;73(3–4):113–22.
Wen J, Trolle C, Viuff MH, Ringgaard S, Laugesen E, Gutmark EJ, et al. Impaired aortic distensibility and elevated central blood pressure in Turner syndrome: a cardiovascular magnetic resonance study. J Cardiovasc Magn Reson. 2018;20(1):80.
Boucher J, Masri B, Daviaud D, Gesta S, Guigné C, Mazzucotelli A, et al. Apelin, a newly identified adipokine up-regulated by insulin and obesity. Endocrinology. 2005;146(4):1764–71.
Sinitsky MY, Dyleva YA, Uchasova EG, Belik EV, Yuzhalin AE, Gruzdeva OV, et al. Adipokine gene expression in adipocytes isolated from different fat depots of coronary artery disease patients. Arch Physiol Biochem. 2022;128(1):261–9.
Yan L, Jin W, Zhao Q, Cui X, Shi T, Xu Y, et al. PWWP2B fine-tunes adipose thermogenesis by stabilizing HDACs in a NuRD subcomplex. Adv Sci (Weinh). 2021;8(16):e2102060.
Kaushansky K, Broudy VC, Lin N, Jorgensen MJ, McCarty J, Fox N, et al. Thrombopoietin, the Mp1 ligand, is essential for full megakaryocyte development. Proc Natl Acad Sci USA. 1995;92(8):3234–8.
Viuff MH, Stochholm K, Grønbaek H, Berglund A, Juul S, Gravholt CH. Increased occurrence of liver and gastrointestinal diseases and anaemia in women with Turner syndrome - a nationwide cohort study. Aliment Pharmacol Ther. 2021;53(7):821–9.
Schweighofer N, Haudum CW, Schmidt A, Kolesnik E, Colantonio CR, Pieske B, Pieber TR, Obermayer-Pietsch B. Diabetes-, sex-, and BMI specific associations of genetic variants in PRRG1 with cardiovascular surrogates in a large cohort at CV risk. 21st European Congress of Endocrinology. Lyon: Endocrine Abstracts; 2019.
SLC16A2 solute carrier family 16 member 2 [Homo sapiens (human)]. 2022. Available from: https://www.ncbi.nlm.nih.gov/gene/6567.
Sarret C, Oliver Petit I, Tonduti D, et al. Allan-Herndon-Dudley Syndrome. In: Adam MP, Ardinger HH, Pagon RA, Wallace SE, Bean LJH, Gripp KW, et al., editors. GeneReviews(®). Seattle: University of Washington, Seattle; 1993. (Copyright © 1993-2022).
Cox TC, Allen LR, Cox LL, Hopwood B, Goodwin B, Haan E, et al. New mutations in MID1 provide support for loss of function as the cause of X-linked Opitz syndrome. Hum Mol Genet. 2000;9(17):2553–62.
Robin NH, Opitz JM, Muenke M. Opitz G/BBB syndrome: clinical comparisons of families linked to Xp22 and 22q, and a review of the literature. Am J Med Genet. 1996;62(3):305–17.
Mortensen KH, Andersen NH, Gravholt CH. Cardiovascular phenotype in Turner syndrome–integrating cardiology, genetics, and endocrinology. Endocr Rev. 2012;33(5):677–714.
Lan F, Bayliss PE, Rinn JL, Whetstine JR, Wang JK, Chen S, et al. A histone H3 lysine 27 demethylase regulates animal posterior development. Nature. 2007;449(7163):689–94.
Mansour AA, Gafni O, Weinberger L, Zviran A, Ayyash M, Rais Y, et al. The H3K27 demethylase Utx regulates somatic and germ cell epigenetic reprogramming. Nature. 2012;488(7411):409–13.
Lindgren AM, Hoyos T, Talkowski ME, Hanscom C, Blumenthal I, Chiang C, et al. Haploinsufficiency of KDM6A is associated with severe psychomotor retardation, global growth restriction, seizures and cleft palate. Hum Genet. 2013;132(5):537–52.
Piton A, Gauthier J, Hamdan FF, Lafrenière RG, Yang Y, Henrion E, et al. Systematic resequencing of X-chromosome synaptic genes in autism spectrum disorder and schizophrenia. Mol Psychiatry. 2011;16(8):867–80.
Pena JLB, Santos WC, Siqueira MHA, Sampaio IH, Moura ICG, Sternick EB. Glycogen storage cardiomyopathy (PRKAG2): diagnostic findings of standard and advanced echocardiography techniques. Eur Heart J Cardiovasc Imaging. 2021;22(7):800–7.
Winder WW. Energy-sensing and signaling by AMP-activated protein kinase in skeletal muscle. J Appl Physiol (1985). 2001;91(3):1017–28.
Hu J, Tang B, Wang J, Huang K, Wang Y, Lu S, et al. Familial atrial enlargement, conduction disorder and symmetric cardiac hypertrophy are early signs of PRKAG2 R302Q. Curr Med Sci. 2020;40(3):486–92.
Beyzaei Z, Ezgu F, Geramizadeh B, Alborzi A, Shojazadeh A. Novel PRKAG2 variant presenting as liver cirrhosis: report of a family with 2 cases and review of literature. BMC Med Genomics. 2021;14(1):33.
Acknowledgements
We thank Lone Kvist and Pamela Celis for their technical assistance.
Funding
A.S. was supported by the Independent Research Fund Denmark (0134-00130B). M.H.V. was supported by Aarhus University. C.H.G. was supported by Novo Nordisk Foundation (NNF15OC0016474, NNF20OC0060610), “Fonden til lægevidenskabens fremme”, the Familien Hede Nielsen foundation, and the Independent Research Fund Denmark (0134-00406A).
Author information
Authors and Affiliations
Contributions
M.V., A.S., J.J., and C.H.G. designed research. A.S., S.C., S.B.P., K.M.L., M.G.B.P., and C.T. performed research. M.V., A.S., E.B.J., and J.J. analyzed data. M.V., A.S., E.B.J., J.J., and C.H.G. wrote the paper. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The study was approved by The Danish Data Protection Agency and the local and central ethics.committee (Region Midtjylland, Denmark number M-20080238 and M-20010248, Central
Denmark Regional Committee on Health Research Ethics number 1–10-72–131-15) and registered at ClinicalTrials.gov (NCT00624949, NCT00999310, NCT02526628, NCT01678261). The study was conducted according to the principles expressed in the Declaration of Helsinki. All participants provided informed consent.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Figure S1.
Manhattan plots of X chromosomal differentially expressed genes. Figure S2. Venn diagrams of X chromosomal differentially expressed genes. Figure S3. Circular plots of PAR genes. Figure S4. Manhattan plots of autosomal differentially expressed genes. Figure S5. Densityplot of autosomal differentially expressed genes (aDEGs). Figure S6. Venn diagrams of autosomal differentially expressed genes. Figure S7. WGCNA analysis module 1–3. Figure S8. WGCNA analysis trait 1–3. Figure S9. WGCNA analysis module 4–6. Figure S10. WGCNA analysis trait 4–5. Figure S11. Venn diagrams of X chromosomal and autosomal differentially methylated sites. Figure S12. Venn diagram of X chromosomal and autosomal differentially methylated sites. Figure S13. Distribution of X chromosomal differentially methylated sites (xDMPs) and autosomal differentially methylated sites (aDMPs) between different contrasts in relation to genic region. Figure S14. Correlation plots of significant differentially methylated positions (DMPs) in correlation with gene expression. Figure S15. Correlation plots of significantly differentially methylated regions (DMRs) in correlation with gene expression.
Additional file 2: Text S1.
Expression profile of escape genes in fat and muscle. Text S2. Enrichment analysis of DEGs. Text S3. Correlation between DMPs and gene expression. Text S4. Correlation between DMRs and DEGs across tissues.
Additional file 3: Table S1.
X chromosomal DEGs between 45,X vs 46,XX in blood. Table S2. X chromosomal DEGs between 45,X vs 47,XXY in blood. Table S3. X chromosomal DEGs between 45,X vs 46,XY in blood. Table S4. X chromosomal DEGs between 46,XX vs 46,XY in blood. Table S5. X chromosomal DEGs between 47,XXY vs 46,XX in blood. Table S6. X chromosomal DEGs between 47,XXY vs 46,XY in blood. Table S7. X chromosomal DEGs between 45,X vs 46,XX in fat. Table S8. X chromosomal DEGs between 45,X vs 47,XXY in fat. Table S9. X chromosomal DEGs between 45,X vs 46,XY in fat. Table S10. X chromosomal DEGs between 46,XX vs 46,XY in fat. Table S11. X chromosomal DEGs between 47,XXY vs 46,XX in fat. Table S12. X chromosomal DEGs between 47,XXY vs 46,XY in fat. Table S13. X chromosomal DEGs between 45,X vs 46,XX in muscle. Table S14. X chromosomal DEGs between 45,X vs 47,XXY in muscle. Table S15. X chromosomal DEGs between 45,X vs 46,XY in muscle. Table S16. X chromosomal DEGs between 46,XX vs 46,XY in muscle. Table S17. X chromosomal DEGs between 47,XXY vs 46,XX in muscle. Table S18. X chromosomal DEGs between 47,XXY vs 46,XY in muscle. Table S19. Autosomal DEGs between 45,X vs 46,XX in blood. Table S20. Autosomal DEGs between 45,X vs 47,XXY in blood. Table S21. Autosomal DEGs between 45,X vs 46,XY in blood. Table S22. Autosomal DEGs between 46,XX vs 46,XY in blood. Table S23. Autosomal DEGs between 47,XXY vs 46,XX in blood. Table S24. Autosomal DEGs between 47,XXY vs 46,XY in blood. Table S25. Autosomal DEGs between 45,X vs 46,XX in fat. Table S26. Autosomal DEGs between 45,X vs 47,XXY in fat. Table S27. Autosomal DEGs between 45,X vs 46,XY in fat. Table S28. Autosomal DEGs between 46,XX vs 46,XY in fat. Table S29. Autosomal DEGs between 47,XXY vs 46,XX in fat. Table S30. Autosomal DEGs between 47,XXY vs 46,XY in fat. Table S31. Autosomal DEGs between 45,X vs 46,XX in muscle. Table S32. Autosomal DEGs between 45,X vs 47,XXY in muscle. Table S33. Autosomal DEGs between 45,X vs 46,XY in muscle. Table S34. Autosomal DEGs between 46,XX vs 46,XY in muscle. Table S35. Autosomal DEGs between 47,XXY vs 46,XX in muscle. Table S36. Autosomal DEGs between 47,XXY vs 46,XY in muscle.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Viuff, M., Skakkebæk, A., Johannsen, E.B. et al. X chromosome dosage and the genetic impact across human tissues. Genome Med 15, 21 (2023). https://doi.org/10.1186/s13073-023-01169-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13073-023-01169-4