Medicine, Faculty of
Other UBC
Medical Genetics, Department of
DSpace
Epigenetics & Chromatin. 2021 Feb 17;14(1):12
Children's Hospital (Vancouver, B.C.). Research Institute
University of British Columbia. Centre for Molecular Medicine and Therapeutics
The Author(s)
Balaton, Bradley P.
Fornes, Oriol
Wasserman, Wyeth W.
Brown, Carolyn Janet
Background: X-chromosome inactivation (XCI) in eutherian mammals is the epigenetic inactivation of one of the two X chromosomes in XX females in order to compensate for dosage differences with XY males. Not all genes are inactivated, and the proportion escaping from inactivation varies between human and mouse (the two species that have been extensively studied). Results: We used DNA methylation to predict the XCI status of X-linked genes with CpG islands across 12 different species: human, chimp, bonobo, gorilla, orangutan, mouse, cow, sheep, goat, pig, horse and dog. We determined the XCI status of 342 CpG islands on average per species, with most species having 80–90% of genes subject to XCI. Mouse was an outlier, with a higher proportion of genes subject to XCI than found in other species. Sixteen genes were found to have discordant X-chromosome inactivation statuses across multiple species, with five of these showing primate-specific escape from XCI. These discordant genes tended to cluster together within the X chromosome, along with genes with similar patterns of escape from XCI. CTCF-binding, ATAC-seq signal and LTR repeats were enriched at genes escaping XCI when compared to genes subject to XCI; however, enrichment was only observed in three or four of the species tested. LINE and DNA repeats showed enrichment around subject genes, but again not in a consistent subset of species. Conclusions: In this study, we determined XCI status across 12 species, showing mouse to be an outlier with few genes that escape inactivation. Inactivation status is largely conserved across species. The clustering of genes that change XCI status across species implicates a domain-level control. In contrast, the relatively consistent, but not universal correlation of inactivation status with enrichment of repetitive elements or CTCF binding at promoters demonstrates gene-based influences on inactivation state. This study broadens enrichment analysis of regulatory elements to species beyond human and mouse.
Balaton et al. Epigenetics & Chromatin (2021) 14:12 https://doi.org/10.1186/s13072-021-00386-8RESEARCHCross-species examination of X-chromosome inactivation highlights domains of escape from silencingBradley P. Balaton1, Oriol Fornes1,2,3, Wyeth W. Wasserman1,2,3 and Carolyn J. Brown1* Abstract Background: X-chromosome inactivation (XCI) in eutherian mammals is the epigenetic inactivation of one of the two X chromosomes in XX females in order to compensate for dosage differences with XY males. Not all genes are inactivated, and the proportion escaping from inactivation varies between human and mouse (the two species that have been extensively studied).Results: We used DNA methylation to predict the XCI status of X-linked genes with CpG islands across 12 different species: human, chimp, bonobo, gorilla, orangutan, mouse, cow, sheep, goat, pig, horse and dog. We determined the XCI status of 342 CpG islands on average per species, with most species having 80–90% of genes subject to XCI. Mouse was an outlier, with a higher proportion of genes subject to XCI than found in other species. Sixteen genes were found to have discordant X-chromosome inactivation statuses across multiple species, with five of these show-ing primate-specific escape from XCI. These discordant genes tended to cluster together within the X chromosome, along with genes with similar patterns of escape from XCI. CTCF-binding, ATAC-seq signal and LTR repeats were enriched at genes escaping XCI when compared to genes subject to XCI; however, enrichment was only observed in three or four of the species tested. LINE and DNA repeats showed enrichment around subject genes, but again not in a consistent subset of species.Conclusions: In this study, we determined XCI status across 12 species, showing mouse to be an outlier with few genes that escape inactivation. Inactivation status is largely conserved across species. The clustering of genes that change XCI status across species implicates a domain-level control. In contrast, the relatively consistent, but not uni-versal correlation of inactivation status with enrichment of repetitive elements or CTCF binding at promoters demon-strates gene-based influences on inactivation state. This study broadens enrichment analysis of regulatory elements to species beyond human and mouse.Keywords: X-chromosome inactivation, Cross-species, DNA methylation, Escape from X-chromosome inactivation, CpG islands, Dosage compensation, Mammals, ATAC-seq, CTCF, Repetitive elements© The Author(s) 2021. 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://creat iveco mmons .org/licen ses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creat iveco mmons .org/publi cdoma in/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.BackgroundHuman and mouse differ in both the initiation and com-pleteness of X-chromosome inactivation (XCI) [1, 2]. In contrast to human, mouse has imprinted XCI early in development, which is maintained in extraembryonic (placental) tissues [3–5]. In placenta, rat [6] and vole [7] also have imprinted XCI while horse/donkey hybrids [8] and pig [9] have random XCI. The story is unclear in cow, Open AccessEpigenetics & Chromatin*Correspondence: carolyn.brown@ubc.ca1 Department of Medical Genetics, The University of British Columbia, Vancouver, CanadaFull list of author information is available at the end of the articlePage 2 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 where both random [10] and imprinted [11] XCI have been reported. At the blastocyst stage, human as well as rabbit express XIST (the RNA that initiates the silencing cascade) from both alleles, while mouse has exclusively paternal Xist expression [1]. Cow has been observed to upregulate XIST at a similar stage to human and rab-bit [12]. Human and rabbit also showed later inactiva-tion timing than mouse [1]. See [13] for a review of XCI across species.Not all genes are subject to XCI, and here again, there is a substantial difference between human and mouse. Escape from XCI is generally defined as having an inac-tive X (Xi) expression of at least 10% of active X (Xa) expression [14]. Around 12% of X chromosome genes are escaping XCI in human [15], while in mouse the pro-portion of genes escaping from XCI is only 3–7% [16]. In human, an additional 15% of genes variably escape from XCI, differing in their XCI status between different tis-sues, populations, individuals or studies [15, 17]. Large-scale studies have not been reported in species outside of human and mouse, and the studies in mouse generally report only on the genes escaping from XCI. The varia-tion between species highlights the importance of study-ing XCI across a range of species; particularly as the most common model organism, mouse, appears quite different from human.There are various methods to examine the XCI sta-tus of genes, with the above numbers being determined using a combination of allelic expression and DNA methylation (DNAme). Additional methods to assess XCI status are reviewed in [18]. For allelic expression to be used to examine XCI escape status, the samples ana-lyzed must be skewed so that the majority of cells in the sample have the same Xi. Skewing of XCI > 90% occurs infrequently in human, but at elevated incidence in blood [19] and cancer due to its monoclonal origin [20]. Cell lines that have undergone clonal selection or which are skewed due to X-linked diseases have also been used [14]. Mouse lines with the gene that controls initiation of XCI, Xist, knocked out on one allele exclusively inactivate the X chromosome with functional Xist [16]; and selectable markers such as fluorescent proteins can also be inserted on one of the X chromosomes in order to select for cell populations with a consistent Xa [21]. Trophoblast cells in mouse have imprinted XCI, and have also been used to determine XCI status [22]. Overall, the requirement for skewing of XCI dramatically limits the datasets that can be used to analyze escape from XCI using allelic expression.DNAme-based analyses circumvent this challenge. DNAme of CpG islands at promoters is strongly predic-tive of XCI escape status [23]. CpG islands are regions of at least 200 bp with high GC content and limited depletion of CG dinucleotides, and are often associated with the promoters of genes, particularly house-keeping genes [24]. Males have low DNAme of promoter CpG islands on the X chromosome, while females, with one Xa and one Xi, will have one relatively unmethylated chromosome and one methylated chromosome, for an average methylation level around 50%. DNAme in gene bodies differs between genes escaping from and subject to XCI, but these differences are subtler and may be tis-sue-specific [23, 25–27]).Knowing the XCI status of genes is important, as genes that escape from XCI often have sex-biased expression, being higher in males if a gametolog is also present on the Y, and higher in females if not [17]. Furthermore, hav-ing two active copies of a gene has been argued to pro-tect females from cancers as both copies will need to be mutated in order to have loss of function [28]. In indi-vidual species, knowing which genes escape from XCI will be useful for mapping the effect of X-linked genes to various traits, and understanding XCI within a species is important for genomic selection strategies in breed-ing for agriculture [29]. Additionally, the knowledge of which genes escape from XCI across species can further our understanding of the underlying mechanism allowing some genes to escape XCI and give insight into the evolu-tionary development of XCI.Here, we compared the XCI status of human and mouse, first examining allelic expression and DNAme in human and mouse to establish robust thresholds of DNAme as an indicator of XCI. We then used DNAme data across two separate groups, one of nine different mammalian species, and one of five different primate species, to examine conservation of XCI escape status across species. Finally, we performed an analysis testing elements previously seen enriched at genes with various XCI statuses (repetitive elements, CTCF and ATAC-seq) for enrichment with our XCI status calls across species.ResultsXCI status calls from allelic expressionTo obtain DNAme thresholds separating genes escap-ing XCI from genes subject to XCI, we first needed to establish which genes were escaping versus subject to XCI using allelic expression data. Allelic expression data requires skewed Xi choice and thus was only available for two species: human and mouse (Fig. 1, Additional file 1: Figures S1, S2). Expression-based XCI status calls were determined using a binomial model as previously described [16], with genes having an Xi/Xa expression ratio significantly over 0.1 being called as escaping XCI and those with Xi/Xa significantly under 0.1 being called as subject to XCI. For humans, we obtained data for eight skewed samples from cancer-related samples and Page 3 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 we identified 44 genes escaping XCI, 262 genes subject to XCI and 21 genes variably escaping from XCI in them (Additional file 2: Table S1). We called genes as variably escaping if they had at least 33% of informative samples with each XCI status. The majority of these XCI status calls agreed with previous studies, with discordance for only 53 genes, (17% of genes with an XCI status call in both), 39 of which were reported to variably escape from XCI here or previously [15]. We attribute the low num-ber of genes variably escaping in our current study to the limited number of samples available and the frequency of informative, heterozygous SNPs per sample, result-ing in a mean of 3.5 informative samples per gene. With more samples, we would expect to observe more variably escaping genes.In mouse we classified 16 genes as escaping XCI, 662 genes subject to XCI and 10 genes variably escaping from XCI (Additional file 3: Table S2). We used three different mouse expression datasets (Keown et al., Berletch et al. and Wu et al.) and results were 97%, 90% and 87% con-cordant when datasets were compared with each other [16, 21, 26]. Most of the discordance in our results arises from identifying more genes variably escaping in the Wu dataset than the other two datasets. Additionally, our use of a threshold of 0.1 rather than 0 to call escape from XCI and the inclusion of a variable escape category resulted in more discordant calls relative to those assigned by Berletch [16]. Figure 1 shows a clear DNAme difference between genes with an Xi/Xa expression ratio under this 0.1 threshold and genes with an Xi/Xa expression ratio over the threshold.Establishing thresholds for calling XCI status from DNAmeDNAme data have also been used to call XCI status [23], and is now available from a number of species where expression in individuals with skewed Xi choice is not available. Our search of GEO [30] for DNAme data across eutherian species found datasets with females for 12 dif-ferent species: human, chimp, bonobo, gorilla, orangutan, mouse, cow, sheep, pig, horse, goat and dog (Additional file 4: Table S3). Most of the datasets used whole genome bisulfite sequencing (WGBS), while horse was limited to a reduced representation bisulfite sequencing (RRBS) dataset and many of the primates and dog were processed on the Illumina Infinium Human Methylation450 Bead-Chip array (450k array), with probes that did not map well to the species in question being filtered out by the source publications. Plotting male versus female DNAme at promoter CpG islands on the X chromosome showed similar trends across species (Additional file 1: Figure S3) with a cluster of sites with less than 10% methyla-tion in both, the bulk of sites showing higher female and low male methylation, and the cluster that is over 70% XCI status of past studies:Escapes XCI Subject to XCI XCI status variesbaFig. 1 Using Xi/Xa expression ratio to establish thresholds of DNAme for XCI status calls. Two species are featured: human (a) and mouse (b). Each point is a SNP with Xi/Xa expression data, matched to the closest CpG island within 2 kb of the closest TSS (accounting for which splice variants would include the SNP) in order to have matched DNAme values. Lines are drawn at 0.1 Xi/Xa expression and at 10, 15 and 60% DNAme as they were used as thresholds to call XCI escape status subsequently. Points are colored based on their XCI status calls. For human, previously published XCI status calls were used [15], while in mouse, which did not have studies calling genes as subject to XCI, they were colored based on their Xi/Xa expression-based XCI status calls featured here. Genes in the pseudoautosomal region, which matches to the Y chromosome, were filtered out. CEMT30, a leukemia cancer sample was used for A, while the Keown et al. data were used for BPage 4 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 methylated in both sexes being under-represented on the array data. There are some differences in the amount of male hemi-methylated islands and the female DNAme average across species, which could be due to differences across species or due to the different tissues and methods of assessing DNAme used.DNAme levels for human and mouse were compared to Xi/Xa expression in order to establish thresholds of DNAme for calling escape from XCI (Fig. 1, Addi-tional file 1: Figures S1, S2). There was good correlation between XCI status calls made using Xi/Xa expression and DNAme with a 10% DNAme threshold. An uncall-able zone between 10 and 15% DNAme was added to lower the chance of miscalling genes, as most discordan-cies between Xi/Xa expression-based calls and DNAme-based calls had DNAme levels in this range. DNAme at genes subject to XCI was lower than expected if the Xi was completely hypermethylated, with an average DNAme of 38% and 27% in human and mouse, respec-tively (Table 1). This shows that the DNAme on the Xi is not complete at these CpG islands. Looking at auto-somal imprinted genes, the expected 50% DNAme ratio was found, demonstrating that lower methylation is not a problem inherent with this analysis or datasets, rather it reflects the DNAme levels of the Xi (Additional file 1: Figure S4).XCI status calls from DNAmeApplying our DNAme thresholds across species to make XCI status calls generated between 26 and 567 XCI status calls per species, with a median of 342 calls per species (Additional file 2: Table S1, Additional file 3: Table S2). Most species had 80–90% of genes identified as subject to XCI by DNAme (Fig. 2), while mouse had 95% of genes subject to XCI and horse only had 76% of genes subject to XCI. The decreased number of genes subject to XCI in horse may be due to the data being generated using RRBS, which provides sparser data and, unlike 450k array data, the sparse CpGs assessed are not the same across samples. In other species the average DNAme at genes subject to XCI ranged from 31% in sheep to 41% in the chimp 450k array data. The 450k array data tended to have higher DNAme than WGBS data, with values between 38 and 41%. Comparison between human and chimp WGBS and 450k array data at the same genes showed that the WGBS and 450k array data differ in DNAme levels, with R2 values of 0.04 in chimp and 0.59 in human (Additional file 1: Figure S5). Differences may be due to having more CpG sites averaged in the WGBS data. Of the genes that had XCI status calls from both DNAme determining methods, 98% of human genes had the same XCI status calls when analyzed with WGBS or Table 1 Mean DNAme for genes subject to XCI per datasetThe mean DNAme of CpG islands at genes found subject to XCI was calculated per datasetSpecies Data type Average DNAme (%)Human WGBS 38450k array 41Chimp WGBS 35450k array 41Bonobo 450k array 38Gorilla 450k array 39Orangutan 450k array 39Mouse WGBS 27Cow WGBS 37Sheep WGBS 31Goat WGBS 33Pig WGBS 38Horse RRBS 37Dog 450k array 3920040060008 11 161 1 12 4 6 6 2 11 6 6 3 4 3 5 14 2 11 17 5 10 12 5 10 14 12 11 14 14 11 116 1 2 0 0 2 2 6 2 10 4 2 0 4 0 180 96 87 83 95 88 86 89 87 76 85 88 86 82 89 89XCI Status:Species% of callsHumanChimpMouseCowPigSheepGoatDogHorseHumanGorillaBonoboChimpOrangutanHumanMousesllac sutats ICX fo rebmuNabNEscapes XCI Subject to XCI Variably escapes XCIData Type: Ex’n DNAmeXi/Xa WGBS RRBS 450k arrayFig. 2 The number and type of XCI status calls per species. The number of XCI status calls per dataset (a) and the percentage of calls with each XCI status per dataset (b) are shown. Datasets (columns) were sorted by technique used to generate the data. Species names are colored by the type of data used to generate XCI status callsPage 5 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 the 450k array, as did 92% of chimp genes. The largest impact of using the 450k array instead of WGBS was at genes escaping from XCI, which occasionally crossed the threshold to being called subject to XCI, particularly in chimp, likely due to the low sample size in WGBS (only one sample). Many genes were not assigned a call in one of the datasets as they were hypermethylated. XCI status calls made using our DNAme thresholds were generally consistent so we did not discard the 450k array datasets.Horse had elevated numbers of variably escaping genes (10%), which were close to that seen previously in human, while other species (including human) only had 0–5% of genes found variably escaping from XCI. The variation in proportion of variable escape genes seen here could be due to low sample size (in everything except human WGBS), or from our methods of calling variable escape genes being more stringent than previous studies. We required at least 33% of informative samples to have each XCI status before calling a gene as variably escap-ing from XCI, similar to the initial survey of human XCI status by Carrel and Willard [14]. Reducing this require-ment to only 10% of samples increased the number of variably escaping genes found in human to 63—almost a quarter of informative genes. These include 37 new genes called which did not have enough informative samples to be called as escaping or subject to XCI with our initial thresholds, as well as 15 genes which changed from an initial call of escaping XCI (12 genes) or subject to XCI (three genes). Although this lower threshold called more genes, we used our 33% threshold of variable escape calls for subsequent studies as we wished to focus on genes that we were confident changed their XCI status between species, rather than differing levels of variable escape from XCI.Overall, we saw that calls of XCI status using DNAme agreed well with those made using allelic expression, and provided an opportunity to examine XCI across multi-ple species. While WGBS resulted in the most XCI sta-tus calls, 450k array DNAme-based calls were generally concordant. These studies showed an average of 11% of genes escaping from XCI across 12 different species, with mouse being an outlier with only 5% of genes escaping from XCI.Conservation of XCI status calls across speciesXCI status calls per gene were compared across species, focusing on genes that were informative in 4 + species. We observed 267 genes being completely conserved across all informative species, with only eight of these genes escaping from XCI and the rest being subject to XCI. Of the eight conserved XCI escapees, two (DDX3X and KDM6A) have Y homologues across eutherian mam-mals [31], five have Y pseudogenes in human (ARSD, STS, PNPLA4, EIF2S3 and MED14) [32], and one has no known Y homology (CTPS2) (Fig. 3a). To avoid biasing the analysis with the more conserved primates, the spe-cies were grouped into two groups: primates with 450k array data, and other datasets (including the human and chimp WGBS data). A clear difference in conservation of status was seen between these two groups, with 97% of genes having completely conserved XCI status across primates, while only 75% of genes had conserved XCI status across all mammals (Additional file 2: Table S1). Of the genes which were usually subject to XCI (> 75% of informative species subject to XCI), 79% of these had all informative species subject to XCI. Genes that usually escaped from XCI were less concordant, with only 61% of these genes having entirely conserved XCI status across all informative species. A similar trend was seen in the all primates group.There were 16 genes that varied frequently (2 + spe-cies escaping XCI and 2 + species subject to XCI) in the all mammals group and none that varied greatly across primates, again showing the higher similarity in XCI status across closely related species (Fig. 3). Of these 16 genes, four showed primate-specific escape from XCI (RPS4X, CDK16, EIF1AX and GEMIN8) and one showed artiodactyla-specific (cow, sheep, goat, pig) XCI (KDM5C). The pattern of conservation of the other genes variably escaping across species did not match any phylogenetic patterns. The primate-specific escape genes RPS4X and EIF1AX have been shown to have primate-specific retention of their Y homolog while KDM5C, the gene that is subject to XCI only in artio-dactyla has lost its Y homolog in bulls, while retaining it in mouse and primates [31]. We show the WGBS data surrounding the CpG island at the transcription start (See figure on next page.)Fig. 3 Concordant and discordant escape genes across species. Eight genes escape XCI in all informative species (a), while 259 genes were subject to XCI in all informative species (not shown). Discordant genes in two different groups of species were examined, only primates (b) and all mammals (c, limited to only 2 primate species). The intersection of a gene and species is colored based on that gene’s XCI status call in that species. Genes that did not have an XCI status call in a species are colored grey. Only escape genes informative in at least 4 + species were selected for a. Genes were selected for b if they had at least one discordant primate species while genes in c required two XCI statuses with two or more species. To match best across species within groups, 450k array data were prioritized in b and WGBS data were prioritized in c. Genes are organized based on their position on the human X chromosome with a horizontal black line denoting the centromere. Green boxes highlight domains of adjacent genes with similar changes to XCI statuses across speciesPage 6 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 TCEAL1DIAPH2RPS4XRIBC1SMC1AKDM5CCDK16UBA1BCOREIF1AXTXLNGSYAP1GEMIN8OFD1TRAPPC2TCEANChumanchimpmousecowsheepgoat pighorse doggeneXCI statusescapes XCIno callsubject to XCIvariably escapes XCIKDM6ADDX3XMED14EIF2S3CTPS2PNPLA4STSARSDhumanchimpbonobogorillaorangutanmousecowsheepgoat pighorse dogspeciesgeneRAP2CSTAG2DIAPH2SPIN4CXorf38CDKL5OFD1TRAPPC2TCEANChumanchimpbonobogorillaorangutanspeciesgenerouprdant escape genesc all mammals ga conco b primate groupPage 7 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 site (TSS) of the ubiquitous escape gene KDM6A, the artiodactyla-specific subject gene KDM5C and the pri-mate-specific escape gene RPS4X (Fig. 4).CDKL5 was the only gene seen to have more than one discordant species in primates (Fig. 3b), being subject to XCI in the human WGBS data, variable in orangu-tan and the human 450k array data and escaping in chimp and bonobo. In gorilla, CDKL5 appeared sub-ject to XCI, but half of the data were in the uncallable region between 10 and 15% DNAme so it was not called as subject to XCI. Other genes had only one species of primates discordant from the rest, usually gorilla or bonobo.Role for alternative promoter usage in escape from XCIUBA1 was particularly interesting as it has been shown previously in human to have two different TSSs with dif-fering XCI statuses [33]. This pattern of multiple TSSs with differing XCI status was seen also in chimp and horse (although data are sparse in horse) (Fig. 5). In cow, the upstream TSS and CpG island are not annotated, but the region homologous to the human upstream TSS showed a DNAme pattern consistent with a promoter subject to XCI, and in pig the CpG islands are annotated but the gene is not. Similarly, in mouse both TSSs (which are annotated but lack CpG island definition) had female-specific DNAme. Mouse has been shown to have fewer CpG islands than human, with CpG island loss from KDM5CKDM6A RPS4Xa b cKDM6A KDM5C IQSEC2HumanPIN4 RPS4XChimpKDM6A KDM5C IQSEC2 PIN4 RPS4XMouseKdm6a Kdm5cIqsec2 Rps4xCowKDM6A KDM5C RPS4XSex:FemaleMaleGoatKDM5C IQSEC2 RPS4XCpG IslandXCI status:Subject to XCIEscapes XCIPigKDM5CKDM6A (not annotated) RPS4XFig. 4 Featured genes compared across species. Male and female DNAme values are graphed by gene and dataset. KDM6A is featured as it is concordantly escaping across species (a). KDM5C is featured because it is known to escape XCI across species but is here shown to be subject to XCI in artiodactyla (cow, sheep, pig and goat) (b). RPS4X is featured because it is a well-known primate-specific escape gene (c). Male methylation is shown in blue and female in red. Annotated CpG islands are shown under the methylation data in purple. Genes are shown colored by their XCI status with arrows at the TSS pointing in the direction of transcription. All of the methylation data shown are from WGBS. Pig did not have KDM6A annotated, but predictions from other species show it located at this CpG island. Goat did not have a CpG island or hypomethylated region at the annotated KDM6A Page 8 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 the ancestral genome being four times as high in mouse as human [34]. The island is still large enough to see hypomethylation on the Xa so the cutoff for minimum island size may be too high in some species. Overall, the alternative TSSs are conserved across species; however, the XCI status of the downstream TSS changes from escaping from XCI in human, chimp and horse to being subject to XCI in mouse and cow. In humans, both TSSs HumanSex:FemaleMaleCpG IslandChimpMouseCowXCI status:Subject to XCIEscapes XCIHorsePigFig. 5 DNAme across the variably escaping gene UBA1. UBA1 is featured as it has multiple different TSSs with CpG islands that have different XCI statuses. Male methylation is shown in blue and female in red. Annotated CpG islands are shown under the methylation data in orange. Genes are shown colored by their XCI status with arrows at the TSS pointing in the direction of transcription. All of the methylation data shown, except for horse are from WGBS. Horse used RRBS data, which is why the data are so sparsePage 9 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 were always found within the same topologically associ-ated domain (TAD) and sub-TAD. Examining TSS usage in the other genes featured in Fig. 3C, we were able to map the TSS and CpG islands using either the University of California Santa Cruz Genome Browser (UCSC) [35] for that species or using the UCSC liftover tool across species, suggesting that the change in XCI status across species was not due to differences in TSS usage between species.Domains of escape from XCI across speciesLooking at the position of genes escaping XCI along the human X chromosome, we saw that most genes escap-ing XCI clustered into domains on the short arm of the X chromosome, similar to what has been described previ-ously [14]. Ten of the 23 transitions between clusters of genes escaping or variably escaping from XCI and genes subject to XCI fell near TAD boundaries in human [36], again similar to what has been seen previously [37]. These clusters of genes escaping from XCI often matched across species. Genes discordant in more than one species were also often clustered, while the genes discordant in only one species were generally scattered by themselves. Some of the genes within discordant clusters were not featured in Fig. 3 as they were missing data in some species. Only two of the strongly discordant genes featured in Fig. 3 are located on the long arm of the X chromosome and they did not form a cluster.We investigated these domains of changing XCI status further by examining whether the discordant species had altered the chromosomal arrangement of these genes. For the primate-specific region of genes escaping XCI spanning the genes TCEANC to GEMIN8, most species had the same gene order, orientation and flanking genes as observed for human (Additional file 1: Figure S6), although some small changes were observed in gorilla, mouse, cow and sheep. In human and mouse, the two species with Hi-C data, there is a TAD spanning from EGFL6 (which neighbors TCEANC) to GEMIN8, which may coordinate the regulation of this region, although if regulated as a domain, EGFL6 would be expected to also escape XCI in primates. There was no data here giv-ing an XCI status for EGFL6, but a previous study had seen it as subject to XCI in human [38]. Gorilla was the only primate that did not demonstrate escape from XCI across this domain, with only the gene GEMIN8 escap-ing XCI. A small insertion was present in gorilla, but it was outside of the TAD which cast doubt about whether it could be the cause of this discordance from the other primates. None of the structural differences in this region were conserved across species with concordant XCI status; thus, we found no detectable genomic correlate underpinning the change in XCI status. Similar results were found for the other discordant regions.Correlation of features with XCI status across speciesThese genes that transition their inactivation status across species provided a dataset to interrogate for fac-tors underlying establishment of silencing or escape from silencing. We considered various factors pertaining to CpG islands in addition to enrichment of various classes of DNA repeats. No differences were seen in CpG island size, nor CpG and GC content between species with discordant XCI status at specific genes. Differences in islands between all genes escaping from versus subject to XCI per species were seen in some species, but no char-acteristic was seen to be significant after multiple testing correction or in more than one species.Different classes of repeats were tested for correla-tion with genes escaping from versus subject to XCI in human, chimp, mouse, cow, sheep, pig and horse. There were significantly more LINE repeats within 15 kb upstream of genes subject to XCI than for genes escap-ing from XCI in chimp, mouse, sheep and horse (Fig. 6a, Additional file 5: Table S4, t-test, corrected p-val-ues < 0.01). Other repeat classes found enriched across multiple species include LTR, DNA and snRNA repeats, which were enriched at genes escaping XCI in 3 species (Additional file 1: Figure S7). SINE repeats, which have previously been seen enriched at genes escaping from XCI [39], were only found significant in horse, which unexpectedly had more SINE repeats near genes subject to XCI than at genes escaping from XCI. Human still had more SINE repeats near genes escaping XCI than subject to XCI on average, but this difference failed to reach sig-nificance in this study.We compared CTCF-binding signal between genes found escaping vs subject to XCI across species. For this, we predicted the probability of CTCF binding across species by using a DanQ model [40] trained on human CTCF ChIP data from ENCODE [41] and validated on mouse (Additional file 1: Figure S8). There were signifi-cant differences in the amount of CTCF-binding signal within 4 kb of TSSs escaping vs subject to XCI in chimp, bonobo, gorilla, and horse but not in human, gorilla, mouse, cow, sheep, goat or pig (Fig. 6b, Additional file 5: Table S4). All of the species with significant differences had more CTCF-binding signal near genes escaping XCI. We also examined whether there were significant regions in the TCEANC to GEMIN8 cluster of discordant genes which correlated with a change in XCI status across spe-cies, but did not find any differences consistent across species (Additional file 6: Table S5).ATAC-seq is an assay for accessible chromatin [42]. Comparing ATAC-seq signal 250 bp up and Page 10 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 downstream of TSSs across species revealed significant differences in the mean female/male ratio across genes that were escaping vs subject to XCI in human, mouse and pig but not in cow or goat (Fig. 6c, Additional file 5: Table S4). ATAC-seq signal had a higher female/male ratio in genes escaping XCI than genes subject to XCI, as seen previously in human [43], and the same trend existed in species where the differences failed to reach significance. In the species with significant differences in ATAC-seq signal with XCI status, we did not see all tissues showing significant differences (Additional file 1: Figure S9). The differences were significant in the only tissue examined in human, two of the three exam-ined in pig, and one out of ten examined in mouse.Across all species examined, mouse genes appeared uniquely well-silenced. We clustered all species based on their XCI status calls (Additional file 1: Figure S10). The bovids (cow, sheep and goat) as a group clustered together, although mouse clusters with them for an unknown reason. Dog has very sparse data which may explain it clustering as an outlier, but we are unsure of the reason why pig clustered with dog instead of with the more closely related bovids. We observed clear abc* * ** * * ** ** ** ** * **FCTChtiwsnibNSSTfobk4nihtiwCATA-M/FqesSSTfobk51nihtiwdnalsiGpCfobk51nihtiwstaeperforebmuNFig. 6 Enrichment of elements which may be related to XCI status. a The number of repetitive elements of each class within 15 kb of each CpG island, sorted by XCI status. See Figure S7 for the repeat classes not shown here. b CTCF binding in overlapping 200-bp bins was predicted using a DanQ model [40]. The Y axis shows the number of bins with > 50% predicted probability of having CTCF binding within 4 kb of each TSS. c Female/male ATAC-seq signal averaged across samples within 250 bp of each TSS. F/M is female over male. Species with a * have significant differences between genes escaping XCI and those subject to XCI (t-test, adjusted p-value < 0.01). P-values are listed in Additional file 5: Table S4, along with the number of CpG islands or TSSs per XCI status in each species used for each analysisPage 11 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 separation of the primates from most other species due to the large number of primate-specific escape genes.DiscussionEscape from XCI is an important contributor to sex dif-ferences in expression and has even been argued to underlie a male predisposition to cancer [17, 28]. In addi-tion, genes subject to XCI can also have unique effects on phenotype, with some mutations having phenotypic effects only when separate cell populations are expressing two different alleles [44, 45]. Mutations that are delete-rious at the cellular level or affect the region controlling choice of Xi can lead to skewed Xi choice, leaving the individual vulnerable to recessive mutations on the oppo-site X chromosome [46, 47]. Knowing the XCI status of genes is also important for estimating the effect of an X-linked allele in genome- or epigenome-wide associa-tion studies [48, 49] and is important for genetic selection of X-linked genes in agriculture [29].To validate our use of DNAme to call XCI status, we compared expression-based calls with DNAme in human and mouse. The human Xi/Xa expression-based calls had 83% agreement with previous calls, with the discrepan-cies largely in genes variably escaping from XCI [15]. As cancer samples were used to allow Xi/Xa analysis, some epigenetic dysregulation may have occurred [20]. We took human DNAme data from IHEC which included multiple consortia, one of which was mostly cancer sam-ples while the other two were not. DNAme-based XCI status calls were quite similar between the consortia with only one gene being called as escaping in one con-sortia and subject to XCI in another (Additional file 7: Table S6). Our study was further limited by the need for heterozygous polymorphisms, thus with only 8 sam-ples, any mis-regulation may not have been noticeable, or led to false or missed calls of variable escape from XCI. Our human DNAme calls were 94% (WGBS) and 91% (450k array) concordant with previous XCI calls, and the two datasets analyzed here gave calls that were 97% concordant with each other. Of the few XCI status calls that were inconsistent with previous studies, 80% were in genes called as variably escaping from XCI, and are likely due to differences in the population or tissues sampled. While our mouse Xi/Xa expression-based calls had a median 90% concordancy across datasets, we only identified 60–86% of previously identified mouse escape genes, likely due to differences in thresholds between studies. There were no discordancies between our mouse DNAme calls and previous mouse studies; however the genes discordant between our Xi/Xa expression calls and previous mouse studies were not informative in our DNAme calls due to lack of CpG islands. Comparing our mouse DNAme calls to a previous study by Keown et al., which examined DNAme on the X chromosome in mouse brain, revealed no discordancies in genes called as escaping XCI, but there were differences in which genes were informative [26].In this study, we have made an average of 342 XCI sta-tus calls per species, for 12 different species. The propor-tion of genes subject to XCI differs, with most species having 80–90% of genes subject to XCI. The only species with more genes subject to XCI is mouse at 95%, and the only species with fewer was horse at 76%. Additionally, horse had elevated numbers of genes variably escap-ing from XCI (10), while other species only had 0–5% of genes variably escaping from XCI. A meta-analysis in human found 8% of genes variably escaping from XCI and a further 7% as varying between studies [15], while our current study identified 6% variable escape in human by expression and only 2% by DNAme. Our study is con-sistent with a previous study using DNAme to make XCI status calls that did not see many genes consistently vari-ably escaping from XCI [23]. Of the genes previously pre-dicted to variably escape from XCI [15], 69% had no data in this study due to lack of a CpG island and another 10% were hypermethylated in males or females and therefore XCI status could not be determined.Our DNAme analysis found that human genes subject to XCI have promoter CpG DNAme between 38% (in WGBS) and 41% (in 450k array analysis) which agrees with a previous analysis using the 450k DNAme array which showed genes subject to XCI having an average DNAme around 40% [23] (Table 1). Mouse had a lower 27% DNAme average for genes subject to XCI; other mouse studies have not examined genes which are sub-ject to XCI. Other species had DNAme averages in a range between human and mouse, but most were closer to human than mouse. Our DNAme thresholds to call genes as escaping from or subject to XCI were consistent across human and mouse WGBS, but as our data were from different studies using different techniques on dif-ferent tissues in different species there may be variation unaccounted for with our thresholds. However, WGBS and 450k array-based XCI status calls were consistent in both human and chimp and, with a few notable excep-tions, genes had concordant XCI status calls across spe-cies. Past studies of XCI status calls using DNAme in human did not see many differences in DNAme-based XCI status across tissues [23], so different tissues ana-lyzed may not cause many discordancies. Having male DNAme as a control and an upper threshold for calling genes as subject to XCI should reduce the chance of call-ing a gene as subject to XCI if it is instead silenced on both copies of the X in a tissue-specific manner. For the primate and dog samples which used the human 450k DNAme array, only probes which mapped consistently Page 12 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 between the species were kept by the source publica-tions [50, 51], and so these species may be enriched for genes with a conserved XCI status. Utilizing datasets from different studies confounds the species differences with other experimental differences including sample size as well as inclusion of male samples. The lack of male samples in some species prohibited us from filtering out genes that are methylated on the Xa and therefore would never be seen to escape XCI by DNAme.Many of the genes escaping from XCI have previously been seen grouped in domains [37], and here we see these domains conserved across species. Furthermore, we see that many of the genes that change XCI status across species are clustered into domains and many of these domains coincide with TADs in human. These domains suggest escape from XCI may be regulated at a domain level; however, we also see some genes being reg-ulated individually and even separate TSSs for the same gene can have opposite XCI statuses. Individual escape genes are often discordant in a few species. Coincidence of changes in XCI status with loss of Y homology empha-sizes the importance of dosage for determining genes whose escape from XCI is vital to survival. Generally, the TSS is seen to be conserved, even when a gene changes XCI status. Previous studies have suggested that CTCF and YY1 may be enriched near genes escaping from XCI [16, 53, 54]. CTCF has also been seen enriched at bound-aries between domains of genes with opposite XCI sta-tuses [56]. Repeat elements (SINE for genes escaping XCI and LINEs for genes subject to XCI) have also been seen enriched in 100-kb windows around TSSs as well as win-dows 15 kb upstream [39, 52].Our XCI status calls across species also allow us to check conservation of elements that may control XCI. A region escaping XCI in human was still able to escape from XCI when inserted at a mouse region which is nor-mally subject to XCI, showing that the mechanisms con-trolling escape from XCI are conserved and functional across species [55]. We suspect that any elements found to be important in human or mouse research will be con-served across species with the same XCI status; having a variety of mammalian species with XCI status calls gives us a platform to test this hypothesis.We compared DNA repeats and CpG island charac-teristics with XCI status within and across species and found none varied significantly across species per dis-cordant gene, few varied between XCI statuses within a species and none varied between XCI statuses in all species. Previous studies have examined enrichment of repetitive elements across differently sized regions ranging from 15 to 100 kb. The enrichment closer to the promoter may reflect gene-specific control, whereas enrichment across a broader range suggests regulation at the level of domains. These studies have seen enrich-ment of LINE and LTR MLT1K repeats at genes subject to XCI and SINE and MER33 repeats at genes escap-ing from XCI [39, 52]. Here, with a window of 15 kb, we replicated the enrichment for LINE repeats, with SINE repeats failing to reach significance and LTR and DNA repeats (which MLT1K and MER33 belong to) showing the opposite trend of previous studies. However, no ele-ment was consistently found across all species. We also predicted CTCF binding and observed that some species have more CTCF-binding signal around genes escaping XCI than genes subject to XCI as has been seen previ-ously [16, 53, 54]. ATAC-seq signal, which has previously been seen enriched at genes escaping XCI, was also seen enriched here, but again, only in some species [43]. A deeper bioinformatic analysis comparing our XCI status calls to features which differ across species with differing XCI status but are conserved in species with conserved XCI status might identify important regulatory features which control the XCI status of nearby genes or control XCI in general.These XCI status calls may be improved in the future through new techniques such as single-cell RNA-seq (scRNA-seq) which can make expression-based XCI status calls without the need for samples with skewed Xi choice. Cells can be analyzed individually or their Xi choice can be identified and then all of the cells with the same Xi can be pooled. scRNA-seq has also identified variable escape at the cellular level within a tissue [17], with most genes varying based on their Xi choice and one gene (TIMP1) seen to vary randomly with no observed difference in Xi choice between cells with different XCI status. Current scRNA-seq datasets have a limitation of low read depth per cell, which limits the ability to exam-ine lowly expressed genes [57]. Methods to enrich for the 3′ end of genes, such as the Chromium Next GEM Single Cell pipeline, are useful for quantifying expression per gene, but further limits the number of polymorphisms available for study. As sequencing becomes cheaper and scRNA-seq technology continues to develop, scRNA-seq may become the new gold standard for making XCI sta-tus calls.Non-CpG DNAme may allow us to use DNAme to examine XCI status in genes without CpG islands, as this mark is seen enriched in the gene body of transcribed genes [25]. Brain and pluripotent cells have the most abundant non-CpG DNAme, with other tissues having less than 1% non-CpG DNAme [58]. A study across mul-tiple tissues in human found 18% of genes (109 of 612) had female-specific non-CpG DNAme in at least one tis-sue, but of these 66% (72 genes) were only significant in one tissue (usually brain) [27]. Another study, in brain only, found 20% of genes escaping from XCI [25]. These Page 13 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 numbers are higher than other reports of escape, likely due to many of these genes variably escaping from XCI and only escaping from XCI in brain.Improved gene and genome annotations in some of the less well-studied species would enhance our XCI status calls across species. Many of the species examined here had their gene annotations generated bioinformatically using CESAR [59] mapping of human genes instead of being annotated with mRNA from that species. This may not have captured the correct TSS, and if transcription was no longer close to the same CpG island these XCI status calls would be invalid. With better annotations in the future, these datasets could be reprocessed to pro-vide more up-to-date XCI status calls with improved confidence.As mouse has considerably fewer genes escaping from XCI than other species, there may be a better species to use as a model for research related to which genes escape from XCI. Unfortunately, none of the species other than mouse examined here are small or make affordable model systems. Rabbit, for which there was no DNAme data available, has been shown to be more similar to human than mouse in aspects of XCI and may be a good species for further examination [1].ConclusionsOur study has created reference XCI status calls for 12 species, so that labs working with diverse mammalian species will have improved understanding of how their genes of interest are expressed in their species of inter-est. We have again confirmed that mouse has substan-tially fewer genes escaping from XCI than human, and shown that other mammals are more similar to human in this regard. Additionally, we have shown conservation of XCI status across the majority of X-linked genes and highlighted some genes of interest which are discordant across species. Interestingly, many of these discordant genes occur in domains of similarly regulated genes. In the future, we hope to use these XCI status calls to iden-tify elements which are controlling escape from XCI and which are conserved across species, and these discordant genes are ideal candidate regions to investigate.MethodsXi/Xa expression‑based XCI status callsHuman whole genome seq and RNA-seq data were obtained for 11 samples, from the Center for Epigenome Mapping Technologies. This data is from cancer samples, and because cancer has a clonal origin, we anticipated they would show skewing of XCI. Eight of the samples had skewed Xi choice, as could be seen by the majority of genes having an Xi/Xa ratio below 0.1. These samples were from brain, blood, breast and thyroid, however neither of the brain samples had fully skewed Xi choice and could be used in this analysis. Mouse RNA-seq data were obtained from two studies using crosses between two distantly related mouse strains, one of which used an Xist knockout to skew Xi selection [16] and another which used fluorescent markers expressed on each X chromosome to separate cells by Xi choice [21]. These mouse datasets have previously been used to find genes escaping XCI, but most mouse studies do not call genes which are subject to XCI, so they were reanalyzed here.The different species were processed differently due to different starting file types. The human data were pre-aligned, starting as DNA VCF files and RNA bam files. The DNA VCF files were indexed and then filtered to only heterozygous SNPs in exons using the bcftools view tool [60]. A BCF file was made for the expression data using samtools mpileup with the -t DP,AD options, fol-lowed by bcftools filter to filter for depth 30 or higher [61]. The RNA BCF file was then indexed and then bcftools call used to find indels and bcftools view used to filter for quality 30 + calls. In mouse, the data were avail-able as fastq files and were aligned using the MEA pipe-line [62]. The resulting unnormalized big wig files were then quantified at known polymorphisms to determine the number of reads on the Xi and Xa.The levels of each allele in the RNA were then extracted using R and compared at all the heterozygous sites found in the DNA analysis [63]. The ratio between alleles was used for graphing and the error rate determined using a binomial model with an α of 0.05 [16]. Genes were assigned XCI status calls per SNP, with a ratio of 0.1 being used as a threshold between genes escaping and subject to XCI and not giving an XCI status for genes who cross this threshold with their error rates.SNPs were mapped to splice variants which include the SNP and the closest TSS of these was used to con-nect DNAme and Xi/Xa expression for Fig. 1, Additional file 1: Figures S1 and S2.DNAme‑based XCI status callsGEO was searched for all WGBS, RRBS or 450k array data that was in eutherian mammals other than mouse and human. Human data were downloaded from the International Human Epigenomics Consortium (IHEC) [64], while a single mouse dataset with a high number of samples was downloaded [65]. Data were downloaded for Homo sapiens (human), Pan troglodytes (chimp), Pan paniscus (bonobo), Gorilla and Gorilla beringei (gorilla), Pongo pygmaeus and Pongo abelii (orangutan), Mus musculus (mouse), Bos Taurus (cow), Ovis aries (sheep), Capra aegagrus hircus (goat), Sus scrofa (pig), Equus ferus caballus (horse) and Canis familiaris (dog). When processed bigwig files were available they were chosen Page 14 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 over processing from raw data. Relevant genomes were downloaded from UCSC (Additional file 4: Table S3) and raw reads were aligned to them using BISMARK [66]. BISMARK methylation extractor was used to get bedGraph files and then UCSC tools bedGraphToBig-Wig tool used to make bigwig files. Gene and CpG island maps were downloaded from UCSC, and the UCSC tools bigWigAverageOverBed tool was used to quantify the mean methylation level across CpG islands. R was then used to annotate CpG islands within 2 kb of a gene’s TSS as belonging to that gene and XCI status calls were made, with islands with a mean DNAme below 10% being called as escaping XCI and islands with between 15 and 60% DNAme being called as subject to XCI. Islands for which over half of males had 15% DNAme or higher were discarded as having male hypermethylation and being uninformative. The mean DNAme across each sex was also calculated and compared per CpG island. The lack of TSSs mapped within each species precluded robust examination of non-CpG island promoter regions, as we were unsure of the exact location of the TSS.For datasets generated on the human 450k DNAme array, data were downloaded and filtered for promoter-associated probes. The mean DNAme of probes sharing an annotated CpG island were matched to their anno-tated genes and this was used for making XCI status calls as above.ClusteringXCI calls per species were transformed into numeric val-ues, with escape as 0, variable escape as 0.5 and subject to XCI as 1. The daisy function from the cluster package in R was used to compute distance and then hclust with the gower metric and complete method were used to per-form the clustering. The phylogenetic tree was generated using the online interactive Tree of Life tool [67].Conservation analysisR was used to collect and match all the XCI status calls across species. Genes were matched based on their name, controlling only for capitalization changes across species. Genes with XCI status calls in four or more species were included in further analysis. Datasets analyzed were split into two different groups: all mammals (human, chimp, mouse, cow, pig, sheep, and goat WGBS data, with horse RRBS and dog 450k array data) and primates (human, chimp, bonobo, gorilla and orangutan 450k array data). The two separate groups allowed us to examine conserva-tion of genes without our analyses being biased toward primate-specific calls.Statistical testsStatistical tests comparing enrichment of CpG island sta-tistics and various repeat classes between genes subject to or escaping from XCI were done using R. We used a t-test with the Benjamini–Hochberg method for multiple testing correction [68].Domain analysisDomains were identified based on conservation calls above and examined using the UCSC browser to com-pare the arrangement of genes. TAD boundaries were taken from Dixon, 2012 [36] and were annotated to genes if they were between it and the next gene or were within the gene body. Additionally, to confirm that UBA1 TSSs were within the same TAD, we used a larger set of TADs in the 3D genome browser [69].ATAC‑seq analysisATAC-seq data were downloaded, see Additional file 4: Table S3 for data sources. If bigwig files were available they were used, but if not we downloaded raw data and aligned it using HISAT2 [70]. The bamcoverage tool from the deepTools package [71] was used to generate bigwig files (normalized using RPKM) and bigWigAverageOver-Bed from UCSC utilities was used to determine the mean coverage in 250 bp up and downstream of each TSS. Each TSS was matched to the closest CpG island within 2 kb and any XCI status call from that island used for the TSS.CTCF predictionsCTCF binding was predicted using a strand-specific DanQ model [40]. The model was trained on human CTCF ChIP-seq data (i.e., positive sequences) and DNase I hypersensitive sites (i.e., negative sequences) from ENCODE [72]. Presence of a CTCF-binding site on the forward strand was required for positive sequences. Negative sequences were required to match the distribution of %GC content of positive sequences. To evaluate the ability of the model to make CTCF-binding predictions in species other than human, it was validated on mouse CTCF ChIP-seq data (also from ENCODE). We used this model to predict the prob-ability of having a CTCF-bound region per overlap-ping 200 bp bins in all species but dog (which had very few XCI status calls to compare to). The CTCF model, the data used to train and validate it, and the cross-species CTCF-binding predictions on the X chromo-somes of the studied species have been deposited on GitHub (https ://githu b.com/wasse rmanl ab/CTCF/). For the purpose of quantifying CTCF-binding signal per TSS, we counted the number of bins with an over Page 15 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 50% predicted probability of being a CTCF-bound region within 4 kb of each TSS. For our analysis of the TCEANC to GEMIN8 region, we counted the number of bins with over 50% probability of CTCF within each region.Supplementary InformationThe online version contains supplementary material available at https ://doi.org/10.1186/s1307 2-021-00386 -8.Additional file 1: Figure S1. The Xi/Xa expression ratio vs promoter DNAme level in individual human samples. Figure S2. The Xi/Xa expres-sion ratio vs promoter DNAme level in individual mouse samples. Figure S3. Male vs female DNAme across species. Figure S4. A comparison of imprinted genes and genes subject to XCI. Figure S5. Comparison of DNAme data generated using WGBS and the 450 k array. Figure S6. Cross-species comparison of a primate-specific escape domain. Figure S7. Number of repeats within 15kb per TSS for genes subject or escaping XCI across species. Figure S8. Tests on mouse CTCF of our model trained on human CTCF. Figure S9. Mean female/male ATAC-seq signal across samples within 250 bp of TSSs, separated by tissue. Figure S10. Clustering of species by XCI status calls.Additional file 2: Table S1. All XCI status calls made in this study com-pared to human.Table S1. All XCI status calls made in this study compared to human.Additional file 3: Table S2. Individual XCI status calls per dataset. Each sheet is a separate dataset analyzed.Additional file 4: Table S3.The sources of data used in this study. See separate sheets for different types of data used.Additional file 5: Table S4. Enrichment of repeats, CTCF and ATAC-seq at genes escaping vs subject to XCI.Additional file 6: Table S5. The number of predicted CTCF binding sites between genes in a discordant region.Additional file 7: Table S6. DNAme based XCI status calls compared across IHEC consortia.Abbreviations450k array: Illumina Infinium Human Methylation450 BeadChip array; DNAme: DNA methylation; RRBS: Reduced representation bisulfite sequencing; scRNA-seq: Single-cell RNA sequencing; TAD: Topologically associated domain; TSS: Transcription start site; UCSC: University of California Santa Cruz Genome Browser; WGBS: Whole genome bisulfite sequencing; Xa: Active X chromo-some; XCI: X-chromosome inactivation; Xi: Inactive X chromosome.AcknowledgementsWe thank the other members of the Brown and Wasserman labs for helpful comments during the development of this project.The human expression data and some of the human WGBS data were generated by The Canadian Epigenetics, Epigenomics, Environment and Health Research Consortium (CEEHRC) initiative funded by the Canadian Institutes of Health Research (CIHR), Genome BC, and Genome Quebec. Information about CEEHRC and the participating investigators and institutions can be found at http://www.cihr-irsc.gc.ca/e/43734 .html. We would also like to thank the research groups which generated the other sources of data used in this analysis.Authors’ contributionsBPB conducted the analyses, OF provided the CTCF prediction modeling, and all authors contributed to the interpretation of data and have approved the manuscript for publication.FundingBPB was supported by a CGS-D award from NSERC. Research was supported by CIHR project grant (PJT-16120).Availability of data and materialsSee Additional file 4: Table S3 and references for data sources used. All sources are publicly available.Ethics approval and consent to participateThe human expression data and some of the human WGBS data were gener-ated by The Canadian Epigenetics, Epigenomics, Environment and Health Research Consortium (CEEHRC) initiative funded by the Canadian Institutes of Health Research (CIHR), Genome BC, and Genome Quebec. Ethics approval for data access was provided by the University of British Columbia Clinical Research Ethics Board (H17-01363).Consent for publicationNot applicable.Competing interestsThe authors declare that they have no competing interests.Author details1 Department of Medical Genetics, The University of British Columbia, Vancou-ver, Canada. 2 BC Children’s Hospital Research Institute, Vancouver, Canada. 3 Centre for Molecular Medicine and Therapeutics, The University of British Columbia, Vancouver, Canada. Received: 4 December 2020 Accepted: 1 February 2021References 1. Okamoto I, Patrat C, Thépot D, Peynot N, Fauque P, Daniel N, et al. Euthe-rian mammals use diverse strategies to initiate X-chromosome inactiva-tion during development. Nature. 2011;472:370–4. 2. Carrel L, Brown CJ. When the Lyon(ized chromosome) roars: ongoing expression from an inactive X chromosome. Philos Trans R Soc Lond B Biol Sci. 2017;372:12. https ://doi.org/10.1098/rstb.2016.0355. 3. Mak W, Nesterova TB, de Napoles M, Appanah R, Yamanaka S, Otte AP, et al. Reactivation of the paternal X chromosome in early mouse embryos. Science. 2004;303:666–9. 4. Okamoto I, Otte AP, Allis CD, Reinberg D, Heard E. Epigenetic dynamics of imprinted X inactivation during early mouse development. Science. 2004;303:644–9. 5. Moreira de Mello JC, de Araújo ESS, Stabellini R, Fraga AM, de Souza JES, Sumita DR, et al. Random X inactivation and extensive mosaicism in human placenta revealed by analysis of allele-specific gene expression along the X chromosome. PLoS ONE. 2010;5:e10947. 6. Wake N, Takagi N, Sasaki M. Non-random inactivation of X chromosome in the rat yolk sac. Nature. 1976;262:580–1. 7. Shevchenko AI, Malakhova AA, Elisaphenko EA, Mazurok NA, Nesterova TB, Brockdorff N, et al. Variability of sequence surrounding the Xist gene in rodents suggests taxon-specific regulation of X chromosome inactiva-tion. PLoS ONE. 2011;6:e22771. 8. Wang X, Miller DC, Clark AG, Antczak DF. Random X inactivation in the mule and horse placenta. Genome Res. 2012;22:1855–63. 9. Zou H, Yu D, Du X, Wang J, Chen L, Wang Y, et al. No imprinted XIST expression in pigs: biallelic XIST expression in early embryos and random X inactivation in placentas. Cell Mol Life Sci. 2019;76:4525–38. 10. Chen Z, Hagen DE, Wang J, Elsik CG, Ji T, Siqueira LG, et al. Global assess-ment of imprinted gene expression in the bovine conceptus by next generation sequencing. Epigenetics. 2016;11:501–16. 11. Xue F, Tian XC, Du F, Kubota C, Taneja M, Dinnyes A, et al. Aberrant patterns of X chromosome inactivation in bovine clones. Nat Genet. 2002;31:216–20. 12. Yu B, van Tol HTA, Stout TAE, Roelen BAJ. Initiation of X Chromosome Inactivation during Bovine embryo development. Cells. 2020;9:2. https ://doi.org/10.3390/cells 90410 16. 13. Shevchenko AI, Dementyeva EV, Zakharova IS, Zakian SM. Diverse developmental strategies of X chromosome dosage compensation in eutherian mammals. Int J Dev Biol. 2019;63:223–33. 14. Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature. 2005;434:400–4.Page 16 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 15. Balaton BP, Cotton AM, Brown CJ. Derivation of consensus inactivation status for X-linked genes from genome-wide studies. Biol Sex Differ. 2015;6:35. 16. Berletch JB, Ma W, Yang F, Shendure J, Noble WS, Disteche CM, et al. Escape from X inactivation varies in mouse tissues. PLoS Genet. 2015;11:e1005079. 17. Tukiainen T, Villani A-C, Yen A, Rivas MA, Marshall JL, Satija R, et al. Landscape of X chromosome inactivation across human tissues. Nature. 2017;550:244–8. 18. Balaton BP, Brown CJ. Escape Artists of the X Chromosome. Trends Genet. 2016;32:348–59. 19. Vacca M, Della Ragione F, Scalabrì F, D’Esposito M. X inactivation and reactivation in X-linked diseases. Semin Cell Dev Biol. 2016;56:78–87. 20. Larson NB, Fogarty ZC, Larson MC, Kalli KR, Lawrenson K, Gayther S, et al. An integrative approach to assess X-chromosome inactivation using allele-specific expression with applications to epithelial ovarian cancer. Genet Epidemiol. 2017;41:898–914. 21. Wu H, Luo J, Yu H, Rattner A, Mo A, Wang Y, et al. Cellular resolution maps of X chromosome inactivation: implications for neural development, function, and disease. Neuron. 2014;81:103–19. 22. Calabrese JM, Sun W, Song L, Mugford JW, Williams L, Yee D, et al. Site-specific silencing of regulatory elements as a mechanism of X inactiva-tion. Cell. 2012;151:951–63. 23. Cotton AM, Price EM, Jones MJ, Balaton BP, Kobor MS, Brown CJ. Land-scape of DNA methylation on the X chromosome reflects CpG density, functional chromatin state and X-chromosome inactivation. Hum Mol Genet. 2015;24:1528–39. 24. CpG Islands in vertebrate genomes. J Mol Biol Academic Press. 1987;196:261–82. 25. Lister R, Mukamel EA, Nery JR, Urich M, Puddifoot CA, Johnson ND, et al. Global epigenomic reconfiguration during mammalian brain develop-ment. Science. 2013;341:1237905. 26. Keown CL, Berletch JB, Castanon R, Nery JR, Disteche CM, Ecker JR, et al. Allele-specific non-CG DNA methylation marks domains of active chromatin in female mouse brain. Proc Natl Acad Sci U S A. 2017;114:E2882–90. 27. Schultz MD, He Y, Whitaker JW, Hariharan M, Mukamel EA, Leung D, et al. Human body epigenome maps reveal noncanonical DNA methylation variation. Nature. 2015;523:212–6. 28. Dunford A, Weinstock DM, Savova V, Schumacher SE, Cleary JP, Yoda A, et al. Tumor-suppressor genes that escape from X-inactivation contribute to cancer sex bias. Nat Genet. 2017;49:10–6. 29. Couldrey C, Johnson T, Lopdell T, Zhang IL, Littlejohn MD, Keehan M, et al. Bovine mammary gland X chromosome inactivation. J Dairy Sci. 2017;100:5491–500. 30. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 2013;41:D991–5. 31. Bellott DW, Hughes JF, Skaletsky H, Brown LG, Pyntikova T, Cho T-J, et al. Mammalian Y chromosomes retain widely expressed dosage-sensitive regulators. Nature. 2014;508:494–9. 32. Wilson Sayres MA, Makova KD. Gene survival and death on the human Y chromosome. Mol Biol Evol. 2013;30:781–7. 33. Goto Y, Kimura H. Inactive X chromosome-specific histone H3 modifica-tions and CpG hypomethylation flank a chromatin boundary between an X-inactivated and an escape gene. Nucleic Acids Res. 2009;37:7416–28. 34. Jiang C, Han L, Su B, Li W-H, Zhao Z. Features and trend of loss of promoter-associated CpG islands in the human and mouse genomes. Mol Biol Evol. 2007;24:1991–2000. 35. Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The Human Genome Browser at UCSC. Genome Res. 2002;12:996–1006. 36. Dixon JR, Selvaraj S, Yue F, Kim A, Li Y, Shen Y, et al. Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature. 2012;485:376–80. 37. Marks H, Kerstens HHD, Barakat TS, Splinter E, Dirks RAM, van Mierlo G, et al. Dynamics of gene silencing during X inactivation using allele-specific RNA-seq. Genome Biol. 2015;16:149. 38. Cotton AM, Ge B, Light N, Adoue V, Pastinen T, Brown CJ. Analysis of expressed SNPs identifies variable extents of expression from the human inactive X chromosome. Genome Biol. 2013;14:R122. 39. Cotton AM, Chen C-Y, Lam LL, Wasserman WW, Kobor MS, Brown CJ. Spread of X-chromosome inactivation into autosomal sequences: role for DNA elements, chromatin features and chromosomal domains. Hum Mol Genet. 2014;23:1211–23. 40. Quang D, Xie X. DanQ: a hybrid convolutional and recurrent deep neural network for quantifying the function of DNA sequences. Nucleic Acids Res. 2016;44:e107. 41. Davis CA, Hitz BC, Sloan CA, Chan ET, Davidson JM, Gabdank I, et al. The Encyclopedia of DNA elements (ENCODE): data portal update. Nucleic Acids Res. 2018;46:D794-801. 42. Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Meth-ods. 2013;10:1213–8. 43. Qu K, Zaba LC, Giresi PG, Li R, Longmire M, Kim YH, et al. Individuality and variation of personal regulomes in primary human T cells. Cell Syst. 2015;1:51–61. 44. Twigg SRF, Babbs C, van den Elzen MEP, Goriely A, Taylor S, McGowan SJ, et al. Cellular interference in craniofrontonasal syndrome: males mosaic for mutations in the X-linked EFNB1 gene are more severely affected than true hemizygotes. Hum Mol Genet. 2013;22:1654–62. 45. Centerwall WR, Benirschke K. An animal model for the XXY Klinefelter’s syndrome in man: tortoiseshell and calico male cats. Am J Vet Res. 1975;36:1275–80. 46. Mitterbauer G, Winkler K, Gisslinger H, Geissler K, Lechner K, Mannhalter C. Clonality analysis using X-chromosome inactivation at the human androgen receptor gene (Humara). Evaluation of large cohorts of patients with chronic myeloproliferative diseases, secondary neutrophilia, and reactive thrombocytosis. Am J Clin Pathol. 1999;112:93–100. 47. Naumova AK, Olien L, Bird LM, Smith M, Verner AE, Leppert M, et al. Genetic mapping of X-linked loci involved in skewing of X chromosome inactivation in the human. Eur J Hum Genet. 1998;6:552–62. 48. Xu W, Hao M. A unified partial likelihood approach for X-chromosome association on time-to-event outcomes. Genet Epidemiol. 2018;42:80–94. 49. Chen B, Craiu RV, Sun L. Bayesian model averaging for the X-chromo-some inactivation dilemma in genetic association study. Biostatistics. 2020;21:319–35. 50. Epiphanio TMF, Fernandes NC, de Oliveira TF, Lopes PA, Réssio RA, Gon-çalves S, et al. Global DNA methylation of peripheral blood leukocytes from dogs bearing multicentric non-Hodgkin lymphomas and healthy dogs: A comparative study. PLoS ONE. 2019;14:e0211898. 51. Hernando-Herraez I, Prado-Martinez J, Garg P, Fernandez-Callejo M, Heyn H, Hvilsom C, et al. Dynamics of DNA methylation in recent human and great ape evolution. PLoS Genet. 2013;9:e1003763. 52. Wang Z, Willard HF, Mukherjee S, Furey TS. Evidence of influence of genomic DNA sequence on human X chromosome inactivation. PLoS Comput Biol. 2006;2:e113. 53. Giorgetti L, Lajoie BR, Carter AC, Attia M, Zhan Y, Xu J, et al. Structural organization of the inactive X chromosome in the mouse. Nature. 2016;535:575–9. 54. Chen C-Y, Shi W, Balaton BP, Matthews AM, Li Y, Arenillas DJ, et al. YY1 binding association with sex-biased transcription revealed through X-linked transcript levels and allelic binding analyses. Sci Rep. 2016;6:37324. 55. Peeters SB, Korecki AJ, Simpson EM, Brown CJ. Human cis-acting ele-ments regulating escape from X-chromosome inactivation function in mouse. Hum Mol Genet. 2018;27:1252–62. 56. Filippova GN, Cheng MK, Moore JM, Truong J-P, Hu YJ, Nguyen DK, et al. Boundaries between chromosomal domains of X inactivation and escape bind CTCF and lack CpG methylation during early development. Dev Cell. 2005;8:31–42. 57. Moreira de Mello JC, Fernandes GR, Vibranovski MD, Pereira LV. Early X chromosome inactivation during human preimplantation development revealed by single-cell RNA-sequencing. Sci Rep. 2017;7:10794. 58. He Y, Ecker JR. Non-CG Methylation in the Human Genome. Annu Rev Genomics Hum Genet. 2015;16:55–77.Page 17 of 17Balaton et al. Epigenetics & Chromatin (2021) 14:12 • fast, convenient online submission • thorough peer review by experienced researchers in your field• rapid publication on acceptance• support for research data, including large and complex data types• gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year • At BMC, research is always in progress.Learn more biomedcentral.com/submissionsReady to submit your research ? Choose BMC and benefit from: 59. Sharma V, Elghafari A, Hiller M. Coding exon-structure aware realigner (CESAR) utilizes genome alignments for accurate comparative gene annotation. Nucleic Acids Res. 2016;44:e103. 60. bcftools. http://samto ols.githu b.io/bcfto ols/bcfto ols.html. Accessed 9 Apr 2020 61. samtools. samtools/samtools. GitHub. https ://githu b.com/samto ols/samto ols. Accessed 9 Apr 2020 62. Richard Albert J, Koike T, Younesy H, Thompson R, Bogutz AB, Karimi MM, et al. Development and application of an integrated allele-specific pipeline for methylomic and epigenomic analysis (MEA). BMC Genomics. 2018;19:463. 63. Website. R Core Team (2019). R: A language and environment for statisti-cal computing. R Foundation for Statistical Computing, Vienna, Austria. https ://www.R-proje ct.org/. Accessed 9 Apr 2020 64. Bujold D, Morais DA, Gauthier C, Côté C, Caron M, Kwan T, et al. The international human epigenome consortium data portal. Cell Syst. 2016;3:496–9. 65. Duncan CG, Grimm SA, Morgan DL, Bushel PR, Bennett BD, et al. Dosage compensation and DNA methylation landscape of the X chromosome in mouse liver. Sci Rep. 2018;8:10138. 66. Krueger F, Andrews SR. Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics. 2011;27:1571–2. 67. Letunic I, Bork P. Interactive Tree Of Life (iTOL): an online tool for phyloge-netic tree display and annotation. Bioinformatics. 2007;23:127–8. 68. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;12:289–300. https ://doi.org/10.1111/j.2517-6161.1995.tb020 31.x. 69. Wang Y, Song F, Zhang B, Zhang L, Xu J, et al. The 3D genome browser: a web-based browser for visualizing 3D genome organization and long-range chromatin interactions. Genome Biol. 2018;19:151. 70. Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Bio-technol. 2019;37:907–15. 71. Ramírez F, Ryan DP, Grüning B, Bhardwaj V, Kilpert F, Richter AS, et al. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res. 2016;44:W160–5. 72. ENCODE Project Consortium, Moore JE, Purcaro MJ, Pratt HE, Epstein CB, Shoresh N, et al. Expanded encyclopaedias of DNA elements in the human and mouse genomes. Nature. 2020;583:699–710.Publisher’s NoteSpringer Nature remains neutral with regard to jurisdictional claims in pub-lished maps and institutional affiliations.
Article
eng
Reviewed
Vancouver : University of British Columbia Library
BioMed Central
Attribution 4.0 International (CC BY 4.0)
Faculty
X-chromosome inactivation
Cross-species
DNA methylation
Escape from X-chromosome inactivation
CpG islands
Dosage compensation
Mammals
ATAC-seq
CTCF
Repetitive elements
Cross-species examination of X-chromosome inactivation highlights domains of escape from silencing
Text