Impact of flanking chromosomal sequences on localization and silencing by the human non-coding RNA XIST
- Angela D. Kelsey†1,
- Christine Yang†1,
- Danny Leung2, 3,
- Jakub Minks1,
- Thomas Dixon-McDougall1,
- Sarah E.L. Baldry1,
- Aaron B. Bogutz1,
- Louis Lefebvre1 and
- Carolyn J. Brown1Email authorView ORCID ID profile
© Kelsey et al. 2015
Received: 17 April 2015
Accepted: 10 September 2015
Published: 2 October 2015
X-chromosome inactivation is a striking example of epigenetic silencing in which expression of the long non-coding RNA XIST initiates the heterochromatinization and silencing of one of the pair of X chromosomes in mammalian females. To understand how the RNA can establish silencing across millions of basepairs of DNA we have modelled the process by inducing expression of XIST from nine different locations in human HT1080 cells.
Localization of XIST, depletion of Cot-1 RNA, perinuclear localization, and ubiquitination of H2A occurs at all sites examined, while recruitment of H3K9me3 was not observed. Recruitment of the heterochromatic features SMCHD1, macroH2A, H3K27me3, and H4K20me1 occurs independently of each other in an integration site-dependent manner. Silencing of flanking reporter genes occurs at all sites, but the spread of silencing to flanking endogenous human genes is variable in extent of silencing as well as extent of spread, with silencing able to skip regions. The spread of H3K27me3 and loss of H3K27ac correlates with the pre-existing levels of the modifications, and overall the extent of silencing correlates with the ability to recruit additional heterochromatic features.
The non-coding RNA XIST functions as a cis-acting silencer when expressed from nine different locations throughout the genome. A hierarchy among the features of heterochromatin reveals the importance of interaction with the local chromatin neighborhood for optimal spread of silencing, as well as the independent yet cooperative nature of the establishment of heterochromatin by the non-coding XIST RNA.
To avoid a functional gene dosage imbalance between the sexes, one of the two X chromosomes in female placental mammals is transcriptionally silenced . This process of X-chromosome inactivation (XCI) occurs early in development and is generally random in all human tissues with either the paternal or maternal X chromosome becoming the inactive X (Xi). The X-inactivation centre (XIC), which is located at Xq13 in humans, is the region of the X that is necessary for the chromosome to be inactivated, and contains the XIST gene (X-inactive specific transcript) [2–4]. Remarkably, the approximately 17 kb spliced and polyadenylated long non-coding XIST RNA uniquely localizes to the chromosome from which it is transcribed . The coating of the Xi by the XIST RNA results in a substantial epigenetic transformation, losing epigenetic modifications associated with active chromatin (notably histone acetylation) and gaining modifications associated with inactive chromatin (including H3K27me3, H3K9me2/3, H4K20me1, and H2AK119u1). The Xi also becomes enriched in several other proteins, including the histone variant macrohistone H2A (macroH2A), the nuclear matrix protein hnRNPU, and the epigenetic regulators SMCHD1 and ASH2L ([6, 7], reviewed in ). In addition, the Xi is peripherally or perinucleolarly localized  with perinucleolar targeting during S phase suggested to be important for maintenance of silencing. A further feature of the Xi is silencing of repetitive elements, as visualized by loss of RNA hybridization with a Cot-1 probe for repetitive DNA  resulting in what has been termed a ‘Cot-1 hole’.
The timing of acquisition of these features has been best studied in mouse, where the differentiation of embryonic stem cells (ESCs) provides an in vitro model for the events of XCI (reviewed in ). Early studies in mouse suggested the presence of a developmental window beyond which Xist was unable to induce X-chromosome inactivation , although macroH2A could be recruited . After this stage, Xist expression was no longer required for maintenance of silencing , consistent with studies in human showing maintenance of silencing in the absence of XIST [15, 16]. SATB1 has been suggested to be involved in defining such a window for Xist function ; however, Satb1/Satb2-deficient mice are able to undergo X inactivation . More recent studies have shown that induction of Xist can recruit H3K27me3 in mouse somatic cells , and XIST induction recruits multiple features of the Xi in human somatic cells [10, 20]. In addition, an ongoing role for Xist in stable silencing of the Xi has been shown by loss of Xist resulting in gene reactivation, loss of perinucleolar association and loss of H3K27me3 . Loss of X-linked gene silencing is enhanced by disruption of DNA methylation and other pathways that cooperate with Xist, thus reactivation of X-linked markers has been used to identify additional players in the pathway [22–24]. Characterization of the ongoing role for XIST in somatic cells has important implications for disease, as highlighted by a recent study showing that deletion of Xist results in hematological malignancies in female mice due to reactivation of X-linked genes .
There is substantial evidence for cooperativity of multiple silencing pathways in the initial silencing of the chromosome, with XCI able to proceed in the absence of key components of the silencing machinery such as PRC2 [26–28], PRC1 , or macroH2A . In addition to multiple factors cooperating in the process of XCI, different X-linked genes may be silenced (or maintained silent) by different players. For example, mutation of the Smchd1 gene results in loss of DNA methylation and partial to full reactivation of approximately 20 % of the X-linked genes in mouse [31, 32]. Surprisingly, many of the marks of an Xi can be recruited by a mouse transgene containing a deletion that makes the Xist RNA defective in silencing, although often the recruitment is not as effective as seen with a full-length Xist (reviewed in ). This silencing defective Xist RNA is also able to form a Cot-1 hole , consistent with the Cot-1 hole not reflecting X-linked gene silencing , but rather a core of silenced non-coding DNA .
Spread of silencing to autosomal genes has been observed in unbalanced X/autosome translocations; however, the extent of autosomal silencing is highly variable in both humans (reviewed in ) and mice . Silencing of autosomal genes has also been observed upon integration of Xist/XIST transgenes into autosomes, and localization of the RNA to the autosome is able to induce many features of the Xi including nucleolar localization [10, 20, 21, 38]. Recently, an XIST transgene was integrated into chromosome 21 in induced pluripotent stem cells from an individual with Down’s syndrome and corrected gene expression from chromosome 21 to near normal disomic levels . Together, these studies demonstrate the XIST RNA is able not only to spread along autosomal material but also to recruit some of the heterochromatic features associated with XCI to autosomes. The silencing of the trisomic chromosome 21 in induced pluripotent stem cells was proposed as a first step towards ‘chromosome therapy’ , and for such uses of XIST a better understanding of the influence of the chromatin neighborhood is necessary.
XIST expression is required to induce the cascade of changes that cooperatively silence the X, but relatively little is known about the process by which the non-coding RNA recruits these changes. In human somatic cells we have reported the recruitment of several features of XCI, including gene silencing, following XIST expression from an inducible transgene [20, 40]. The separation of XIST expression from the myriad of changes that occur during differentiation provides an opportunity to dissect the role of XIST in XCI. By using this inducible XIST transgene to examine the influence of XIST expression induced from nine different integration sites we aimed to establish a hierarchy to the features that are recruited by XIST and determine if any were influenced by the genomic context of the XIST integration. We were able to identify features (Cot-1 hole formation, perinucleolar localization, proximal reporter silencing, and H2AK119u1) that are recruited to all integrations and thus are compatible with any of the sites/genomic contexts tested, a feature that is not recruited to any site examined (H3K9me3) while another set of features (macroH2A, H3K27me3, SMCHD1, and H4K20me1 recruitment) are dependent on genomic context, but appear independent of each other.
XIST RNA expression results in depletion of Cot-1 RNA and increased perinucleolar location
In the HT1080 cells an inducible XIST cDNA has been previously reported to localize in cis to the chromosome from which it is transcribed, resulting in the recruitment of some chromatin modifications and repression of both flanking reporter genes as well as flanking endogenous genes [20, 40]. To test whether there are differences in the functionality of the XIST RNA due to influences of local DNA sequences and chromatin environment we integrated a single copy of the full-length inducible XIST cDNA into nine different FRT sites in HT1080 cells containing a constitutive Tet-repressor transgene to allow induction of XIST by treatment of the cells with the tetracycline analog doxycycline (DOX). The sites of the FRT integrations were identified by inverse PCR and sequencing to be 1p, 3q, 4q, 7p, 7q, 8p, 12q, 15q, and an X chromosome FRT site at Xq23 was previously reported . The integration sites (described in Additional file 1) are in both G-dark and G-light regions (four and five, respectively), and three integrations occurred within a gene.
Increased nucleolar association of chromosomes expressing transgenic XIST
XIST RNA expression silences nearby reporter and endogenous genes
To distinguish differences between clones that are attributable to the integration site rather than clone-to-clone variation, we analyzed clones that were independent integrations of the XIST cDNA into the same genomic location. Intriguingly, for the 12q integration site, we observed a significant correlation between the amount of silencing and the level of XIST expression (Fig. 3c). This correlation was observed for both cultures of the same single-cell clone and additional independent clones at the same integration site. Such a correlation was not observed for genes silenced by XIST expression from the 8p integration site using multiple cultures and two independent clones (Fig. 3d). Thus, while variation in the level of XIST expression occurred within a clone, the impact of this variation depended on the integration site, and with similar XIST expression a similar extent of silencing was observed between independent integrations into the same FRT site. A second clone from the Xq integration site showed very similar silencing, while a second clone from the 7q integration site showed lower XIST expression levels and failed to silence (Additional files 1 and 3), suggesting that the silencing ability of 7q integrations, like 12q integrations, might be influenced by XIST expression levels. Since the clones at the 8p integration site silenced across a range of XIST levels, it seemed possible that, unlike silencing of the reporter genes, the maintenance of silencing of endogenous genes might not be XIST-dependent, so we analyzed silencing of endogenous flanking genes after removal of DOX for the 8p and Xq integration sites. All genes examined showed partial to complete reactivation (Fig. 3e) 5 days after induction of XIST had ceased, indicating that the silencing observed requires ongoing XIST expression.
RNA-seq confirms differences in silencing capacity of XIST at different integrations
The significant differences between the integration sites in their capacity to silence endogenous genes suggested an important impact of the genomic context of the integration; however, it was possible that by chance the genes chosen to be tested were non-random in their ability to silence. Therefore, we chose to examine three clones with XIST at different integration sites using RNA-seq to generate a detailed view of any variability in silencing capacity. We chose to examine the 8p and 12q integration sites, which had shown the most silencing, but different sensitivities to XIST levels, as well as the 1p integration site, which had shown limited silencing.
Chromatin features of regions silenced by XIST in somatic cells
Recruitment of heterochromatic modifications to the site of XIST localization
Co-localization of histone modifications with induced XIST RNA signal
The only mark examined that was not detected at any site was the H3K9me3 mark, which is often associated with constitutive heterochromatin. H2AK119u1 was seen at all of the four integrations examined, suggesting it is less dependent on local chromatin structure than the other marks which were heavily dependent upon the integration site and surprisingly independent of each other. The X chromosome showed the best ability to recruit all features examined, with the next best integration site also being the one that showed the greatest gene silencing (8p).
We have examined the impact of genomic location of XIST on its ability to alter nuclear ultrastructure, chromatin state, and gene expression in a somatic cell line. We observe considerable heterogeneity in the extent of silencing and the recruitment of chromatin marks, depending upon the integration site, providing us with an opportunity to dissect the interactions between features in a system in which complete recruitment and silencing does not occur. The ability of XIST to localize to all integration sites supports that there is limited sequence specificity for the RNA with the chromosome, in agreement with the spread of inactivation reported in X/autosome translocations, and various mouse Xist transgenes (reviewed in [36, 50]). The lack of sequence specificity is also in line with the recently described proximity-transfer model which suggests that Xist first associates with sequences in physical proximity to the Xist locus and then transfers to gene-rich regions which are topologically associated with the integration site .
Spread of silencing to autosomes in these HT1080 cells is not seen to the extent that is observed for X/autosome translocations [48, 49] or an autosomal XIST transgene in iPS cells , suggesting that bypassing differentiation reduces the ability to inactivate a chromosome. While early mouse studies suggested the presence of a limited developmental window during which silencing could be induced  we see XIST-dependent silencing of endogenous genes in these HT1080 fibrosarcoma cells up to almost 50 Mb from the integration site, although endogenous gene silencing is most prominent closer to the integration site. It is possible that the cancerous origin of these cells has reactivated critical developmental gatekeepers, as specific mouse cancers have been shown to allow Xist-induced silencing [17, 25]. The ability of XIST to induce gene silencing in somatic cells has important implications for cancer cells where rearrangements or reactivation of XIST may bring previously active genes under the influence of XIST.
The other feature that we observed at all integrations examined was H2AK119u1, which is established by the PRC1 complex. We observed H2AK119u1 in the absence of SMCHD1, macroH2A, and H3K27me3 in the 1p integration cells, suggesting that XIST may be able to directly recruit the PRC1 complex, in agreement with the PRC2-independent PRC1 recruitment previously suggested in mice . As the 1p integration site demonstrated limited silencing in the presence of H2AK119u1 recruitment, we conclude that H2AK119u1 is insufficient for the spread of gene silencing. No integration site examined showed recruitment of H3K9me3, although the antibody clearly hybridized to the native Xi in female cells. Therefore it appears that facultative heterochromatic marks can be recruited by XIST in these cells, but the establishment of a ‘locked-in’ silent state required additional layers of developmentally-regulated chromatin condensation. Without such locks on silencing, the continued silencing of both reporter genes and endogenous genes was XIST-dependent, undergoing reactivation upon removal of XIST induction.
Interestingly, no single factor was seen to associate with more robust silencing (Fig. 7), demonstrating considerable redundancy in XIST-inducible silencing pathways. Consistent with previous studies in mice, silencing of endogenous genes in these human cells was observed in the absence of recruitment of macroH2A (7p) , SMCHD1 (3q) , or H3K27me3 (12q) . MacroH2A has been reported to be recruited by induction of mouse Xist expression , and is also lost upon deletion of Xist ; yet intriguingly we observe macroH2A recruitment to be site-dependent, rather than solely XIST-dependent. We observed that macroH2A is not required for XIST-induced silencing, which is consistent with the ability of mice that are knocked-out for both macroH2A1 and macroH2A2 to continue to undergo XCI . The recruitment of SMCHD1 and macroH2A seems to be independent of each other as, for example, the 3q integration site recruits macroH2A but not SMCHD1, while the 7q integration shows some SMCHD1 recruitment, but no macroH2A recruitment. There was a positive correlation between the amount of SMCHD1 recruitment and the final perinucleolar localization (Spearman r = 0.8; P = 0.01), with the proportion of cells positive for SMCHD1 always equal or lower to the percent seen to be perinucleolar. Thus perinucleolar localization may be required but not sufficient to recruit SMCHD1. The strong ability of the X integration to recruit all features confounded the ability to detect correlations; however, it appears that the majority of the features are independently recruited to the chromosome that expresses XIST. For several of the integration sites two or more independent clones with similar XIST levels behaved similarly supporting that the variation we observe is due to the integration site. In addition, the 4q and 12q clones were derived from a different subclone of HT1080 from the other autosomal integration sites; however, they did not appear to be more similar to each other, arguing against the variation arising during subclone generation.
The differential recruitment of chromatin marks depending on integration site suggests that multiple silencing pathways not only work in parallel to promote silencing, but that their recruitment is favored by different underlying DNA sequences, consistent with previous studies using genes that escape from inactivation to identify 12 features of the DNA sequence of the X chromosome that may be involved in the spread or maintenance of XCI .
Overall, we saw variability in the recruitment of chromatin marks between the integration sites, highlighting the importance of the integration site in modulating XIST function. Localization of XIST with a concomitant depletion of Cot-1 RNA, recruitment of H2AK119u1, and a shift to perinucleolar location was seen at all integrations, and thus reflect features established by XIST independent of the local chromatin environment. In contrast, recruitment of SMCHD1, macroH2A, H3K27me3, and H4K20me1 appeared to be strongly influenced by the site of XIST expression. As some silencing was observed at all integrations, this study demonstrates that silencing in human somatic cells can occur in the absence of macroH2A, SMCHD1, and H3K27me3/H4K20me1 recruitment, underscoring the independent but cooperative nature of the X-chromosome inactivation process. The X chromosome demonstrated the most consistent ability to recruit the heterochromatic marks of XCI, consistent with an evolutionary accumulation of DNA features enabling the recruitment of heterochromatic marks to the X chromosome.
Generation and culture of cell lines and identification of the transgene integration site
HT1080 HH1 cells were transfected with pcDNA6/TR, and two subclones (HT1080HH1-2-3 or 2–12) expressing high RNA levels of the Tet-repressor were subsequently transfected with pFRT/lacZeo (Life Technologies) at low concentrations. The full-length inducible XIST cDNA construct  was co-transfected with the pOG44 plasmid expressing Flp recombinase for site-specific recombination into the FRT site followed by Hygromycin selection and confirmation of loss of Zeomycin resistance. Cells were grown at 37 °C with 5 % CO2 in DMEM supplemented with penicillin/streptomycin, non-essential amino acids, and 10 % V/V fetal bovine serum. XIST expression was induced with the addition of 1 μg/mL doxycycline to the culture medium.
Southern blotting identified eight FRT integrations that appeared single copy. Inverse PCR utilizing primers complementary to a sequence within the integrated pFRT plasmid was used to identify the precise integration site of the transgenes in the HT1080 cell lines. The ends of linearized plasmids are subject to exonuclease activity, thus the actual integrated transgene often lacks several hundred of base pairs on each end which was first identified by a series of PCR assays prior to restriction endonuclease digestion with a frequently-cutting restriction endonucleases identified to cut in the remaining plasmid, followed by ligation with T4 DNA ligase (Invitrogen) to create circular DNA molecules. The captured genomic DNA was amplified by nested PCR with primers facing outward from the plasmid fragment for sequencing and the genomic location was identified using the BLAT algorithm . No X-linked integration was identified, so the F55 HT1080 clone from Yan and Boyd was used .
While cancer derived, the HT1080 cells remain diploid with four structural rearrangements detected by spectral karyotyping (46,XY,del (1)(p21), i(3)(p10), i(3)(q10), der(4)t(1;4)(p21;p16), der(5)t(5;5)(p15;?), der(11)t(3;11)(q11;q25). Using our allelic pyrosequencing assays we observed instability of chromosome 3 in two of nine clones (see also ) and homozygosity for assays on 4q in one clone, while assays on chromosomes 1, 7, 8, and 15 remained diploid in the clones tested, suggesting that the unbalanced rearranged chromosomes were the least stable in these cells. We analyzed two independent clones at the Xq, 7q, and 8p integration site and five independent clones of the 8p integration site. These clones were individual single-cell clones following Flp-mediated recombination into the FRT site. Each clone was selected for Hyg-resistance and assessed for loss of Zeomycin sensitivity.
As fluorescent reporters allow for efficient screening, we created a plasmid that carries both the inducible repeat A and a DsRed-Express2, driven by the mouse Pgk1 promoter (Fig. 2b). To test whether the ability of repeat A to silence the reporter depends on the genomic integration site, we inserted the repeat A – DsRed-Express2 transgene into six of the HT1080 cell lines with a known chromosomal location of the FRT integration site for assessment of silencing by flow cytometry as previously performed .
RNA FISH and immunofluorescence
Cells were grown on glass coverslips. Upon removal from cell culture the coverslips were first rinsed in ice-cold CSK buffer (0.3 M sucrose, 100 μM NaCl, 10 μM PIPES, 3 μM MgCl2), then permeabilized with 0.5 % Triton-X 100 in CSK for 8 min on ice and then fixed in 4 % paraformaldehyde for 8 min at room temperature. Coverslips were stored at 4 °C in 70 % ethanol. Just prior to RNA FISH, the coverslips were immersed in 100 % ethanol for 5 min and left to air dry. FISH was performed with two probes: an XIST probe and a Cot-1 (Invitrogen) probe. Both probes had been directly fluorescently labeled using the Nick Translation Reagent Kit (Abbott Molecular, Inc.) with Spectrum red-UTP (Vysis) for Cot-1 DNA probes and Spectrum green-UTP (Vysis) for XIST probes. Approximately 150 ng of each probe was mixed together along with 20 μg salmon testes DNA then air dried in a speed vacuum, resuspended in 10 μL deionized formamide, denatured at 80 °C for 10 min, and then mixed with 10 μL hybridization buffer (20 mg/mL BSA, 4XSSC, 20 %). This was pipetted onto a small square of Parafilm and the coverslip was placed on top of the probe mixture. Another piece of Parafilm was then placed on top and the edges were sealed to prevent the drying out of the coverslip. Hybridization took place overnight in a humidified chamber at 37 °C. The next day the coverslips were rinsed as follows: 20 min in 50 % formamide/50 % 4XSSC at 37 °C, 20 min in 2X SSC at 37 °C, and 20 min in 1X SSC at room temperature. Coverslips were then stained with DAPI and mounted onto microscope slides with Vectashield (Vector Laboratories). Cells were observed on a Leica inverted microscope (DMI 6000B) at 100X magnification and images were obtained using a Retiga 4000R (Q-Imaging) camera with Openlab software (PerkinElmer). Images were processed using Adobe Photoshop CS4 to reduce background and correct for variation in FISH efficiency between different images. A one-way ANOVA test in GraphPad was used to determine significantly different signal sizes. Line scans were generated using Image J software (NIH) by drawing a line through the area of interest and plotting the RGB intensities across the line.
To determine the nuclear location, XIST signals were scored visually in Photoshop as being either ‘perinuclear only’ (for example, adjacent to and in contact with the nuclear periphery), ‘perinucleolar only’ (that is, adjacent to and in contact with a Cot-1 negative nucleolus), ‘both’ (that is, adjacent to and in contact with both the nuclear periphery and a Cot-1 negative nucleolus), or ‘neither’. XIST signals scored as ‘both’ are included in the ‘perinuclear’/‘perinucleolar’ percentages in Table 1. For the ‘5d DOX’ counts, results for each integration site are the average of three independent experiments performed on different coverslips, by at least two independent observers, with a minimum of 50 cells counted each time. The ‘No DOX’ counts were done once with a minimum of 60 cells counted per integration site.
For combined RNA FISH and immunofluorescence, coverslips, which had been stored at 4 °C in 70 % ethanol, were first rinsed in PBS then placed onto a small amount of PBT (PBS with 1 % BSA and 0.1 % Tween 20) containing 0.4 U/μL Ribolock RNase inhibitor. Coverslips were sealed between two layers of Parafilm and left in the blocking buffer at room temperature for 20 min, then transferred from blocking buffer to PBT containing 1:100 primary antibody and 0.4 U/μL Ribolock, sealed between two layers of Parafilm and left at room temperature for 4–6 h. Coverslips were then washed three times, for 5 min each time, at room temperature in PBS containing 0.1 % Tween 20, then put onto a small amount of PBT containing 1:250 fluorescently labeled secondary antibody and 0.4 U/μL Ribolock, sealed between two layers of Parafilm and left at room temperature in the dark for 45 min. Coverslips were then washed three times, for 5 min each time, at room temperature in the dark in PBS containing 0.1 % Tween 20. Coverslips were then fixed in 4 % PFA in PBS for 10 min at room temperature in the dark, and washed for 5 min in PBS before continuing on to RNA FISH, making sure that the coverslips remained in the dark throughout the RNA FISH procedure. Antibodies used in immunofluorescence include: anti-H3K27me3 (07–449 from Millipore); anti-macroH2A (07–219 from Millipore); anti-SMCHD1 (ab31865 from Abcam); anti-H2AK119u1 (05–678 from Millipore); anti-H3K9me3 (07–442 from Millipore); anti-H4K20me1 (07–440 from Millipore).
RNA isolation, reverse transcription Q-PCR and allelic discrimination by pyrosequencing
RNA was isolated from cell pellets stored at −70 °C using TRIZOL (Invitrogen) according to the manufacturer’s instructions and then treated with DNase1. cDNA was generated in the range of 0.5–2.5 μg RNA using M-MLV reverse transcriptase for qPCR on a StepOnePlusTM Real-Time PCR System (Applied Biosystems, Darmstadt, Germany), using Maxima Hot Start Taq (Thermo Scientific) and EvaGreen dye (Biotium). The following conditions were used: 95° for 5 min, followed by 40 cycles of (95° for 15 s, 60° for 30 s, 72° for 1 min), and a melt curve stage of (95° for 15 s, 60° for 1 min, increase of 0.3° until 95°). The expression levels of genes of interest were normalized to the expression of ACTB or PGK1. Primer sequences are found in Additional file 6.
Pyrosequencing of cDNA before and after DOX induction of DNA, and of clones containing integrations on alternate chromosomes were examined. Each 25 μL PCR was performed with 1x PCR Buffer (Invitrogen), 0.2 mM dNTPs, 0.625 U Taq DNA polymerase (Invitrogen), 0.5 μM forward and reverse primers, and 50–100 ng of cDNA for (94C for 30 s, 58.3C for 30 s, 72C for 1 min) × 35 cycles, and 72 °C for 10 min for final extension. One of the forward and reverse primers was biotinylated for template isolation during pyrosequencing preparation. Universal M13 primer was also used for some assays (see Additional file 6), where the primer to be biotinylated instead contained the M13 sequence at the 5′end (5′CGC CAG GGT TTT CCC AGT CAC GAC3′). Nested PCR was run for the assays that utilized the universal primers: the first round of PCR was performed with the same cycling conditions but only for 15x cycles, with the M13-tagged primer (without the biotin) and its paired forward or reverse primer; the second round of PCR was run with 1 μL of PCR product from the first round of PCR as the template, as well as the biotinylated M13 primer and the non-M13-tagged primer that was used in the first round, under the same cycling condition but for 20 cycles. Pyrosequencing was performed on the PyroMark MD machine (Qiagen). Template preparation for pyrosequencing was done according to manufacturer’s protocol. For each assay, 10–15 μL of PCR products was used as template and CDT tips were used to dispense the dNTPs.
For both ChIP-seq and RNA-seq, library preparation and sequencing was done according to Illumina protocols, and reads were aligned to hg18 reference genome using Tophat. RNA-seq was performed as previously described , first on two 8p clones following 5 days DOX and the 8p clone with No DOX producing 50 bp paired-end reads, with a Pearson correlation between replicate DOX FPKM values of r = 0.9934 (P <0.0001). A second RNA-seq of the 12q and 1p integration sites with No DOX and following 5 days DOX produced 76 bp paired-end reads. Total RNA expression was quantified in FPKM using Cufflinks. ChIP of H3K27me3 (Active Motif antibody 61017) and H3K27ac (Active Motif antibody 39133) was performed on the 8p clone with and without DOX treatment, as previously done  using 20 μg of chromatin and 3 μg of antibody per IP. ChIP-seq produced 36 bp single-end reads, and the ChIP-seq data were analyzed using SeqMonk to identify enriched H3K27ac peaks (MACS, P value = 1 × 10−3) and to quantify the level of H3K27me3 in 2 kb and 1 Mb non-overlapping windows.
Allelic information for both RNA-seq and ChIP-seq was obtained with a custom workflow in Galaxy [56–58]. Briefly, mapped reads were run through Samtools mpileup to obtain SNPs relative to the reference genome in the control and experimental (integration of interest) samples, and allelic ratio of reads was calculated at each variant. In addition, only the variants with a biallelic ratio of 0.3–0.7 in the No DOX (control) sample were further examined. In order to phase the variants, the allele with lower reads for each variant site in the DOX (experimental) sample was considered to be silenced in cis with the XIST transgene. Variant sites were next combined based on regions of interest, either gene location or genomic regions for RNA-seq and ChIP-seq, respectively. Variant sites that had an allelic ratio 0.15 greater or less than the overall gene allelic ratio for the sample were removed and a new overall gene ratio was calculated. For ChIP-seq, variants with reads >2-fold different from the average number of reads in the region of interest were also excluded. Percent allelic silencing for a given gene or region of interest was defined as (control allelic ratio – experimental allelic ratio)/control allelic ratio. The calculated ratio was the frequency of the allele with lower reads in the experimental sample for each gene or region of interest. We also required that there were at least 4 reads and 5 reads with allelic information per gene (RNA-seq) and region of interest (ChIP-seq), respectively, in both the control and experimental samples to be considered for analysis.
We obtained normalized Hi-C data for human female fibroblast line IMR90  and the Hi-C analysis was done as previous . We calculated the 1 Mb Hi-C counts by summing the counts in all 40 kb bins within each 1 Mb bins across the chromosome, with the anchor being the 1 Mb bin containing the XIST transgene (AGPAT5 for 8p clone).
The sequence data described are available in GSE68109. Additional microscopy images are available in Figshare (http://dx.doi.org/10.6084/m9.figshare.1529822).
This work was supported by CIHR grants to CJB (MOP-13680) and LL (MOP-119357). The authors thank Jen Chow for developing the inducible XIST clones, and Jordan Henriksen, Irene Qi and Allison M. Cotton for assistance and advice on FISH.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lyon MF. Gene action in the X-chromosome of the mouse (Mus musculus L.). Nature. 1961;190:372–3.View ArticlePubMedGoogle Scholar
- Brown CJ, Lafreniere RG, Powers VE, Sebastio G, Ballabio A, Pettigrew AL, et al. Localization of the X inactivation centre on the human X chromosome in Xq13. Nature. 1991;349:82–4.View ArticlePubMedGoogle Scholar
- Brown CJ, Ballabio A, Rupert JL, Lafreniere RG, Grompe M, Tonlorenzi R, et al. A gene from the region of the human X inactivation centre is expressed exclusively from the inactive X chromosome. Nature. 1991;349:38–44.View ArticlePubMedGoogle Scholar
- Brown CJ, Hendrich BD, Rupert JL, Lafreniere RG, Xing Y, Lawrence J, et al. The human XIST gene: analysis of a 17 kb inactive X-specific RNA that contains conserved repeats and is highly localized within the nucleus. Cell. 1992;71:527–42.View ArticlePubMedGoogle Scholar
- Clemson CM, McNeil JA, Willard HF, Lawrence JB. XIST RNA paints the inactive X chromosome at interphase: evidence for a novel RNA involved in nuclear/chromosome structure. J Cell Biol. 1996;132:259–75.View ArticlePubMedGoogle Scholar
- Blewitt ME, Vickaryous NK, Hemley SJ, Ashe A, Bruxner TJ, Preis JI, et al. An N-ethyl-N-nitrosourea screen for genes involved in variegation in the mouse. Proc Natl Acad Sci U S A. 2005;102:7629–34.PubMed CentralView ArticlePubMedGoogle Scholar
- Pullirsch D, Hartel R, Kishimoto H, Leeb M, Steiner G, Wutz A. The Trithorax group protein Ash2l and Saf-A are recruited to the inactive X chromosome at the onset of stable X inactivation. Development. 2010;137:935–43.PubMed CentralView ArticlePubMedGoogle Scholar
- Sado T, Brockdorff N. Advances in understanding chromosome silencing by the long non-coding RNA Xist. Phil Trans Royal Society London Series B. 2013;368:20110325.View ArticleGoogle Scholar
- Barr ML, Carr DH. Correlations between sex chromatin and sex chromosomes. Acta Cytol. 1962;6:34–45.PubMedGoogle Scholar
- Hall LL, Byron M, Sakai K, Carrel L, Willard HF, Lawrence JB. An ectopic human XIST gene can induce chromosome inactivation in postdifferentiation human HT-1080 cells. Proc Natl Acad Sci U S A. 2002;99:8677–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Chow J, Heard E. X inactivation and the complexities of silencing a sex chromosome. Curr Op Cell biol. 2009;21:359–66.View ArticlePubMedGoogle Scholar
- Wutz A, Jaenisch R. A shift from reversible to irreversible X inactivation is triggered during ES cell differentiation. Mol Cell. 2000;5:695–705.View ArticlePubMedGoogle Scholar
- Rasmussen TP, Wutz A, Pehrson JR, Jaenisch R. Expression of Xist RNA is sufficient to initiate macrochromatin body formation. Chromosoma. 2001;110:411–20.View ArticlePubMedGoogle Scholar
- Csankovszki G, Panning B, Bates B, Pehrson JR, Jaenisch R. Conditional deletion of Xist disrupts histone macroH2A localization but not maintenance of X inactivation. Nat Genet. 1999;22:323–4.View ArticlePubMedGoogle Scholar
- Rack KA, Chelly J, Gibbons RJ, Rider S, Benjamin D, Lafreniere RG, et al. Absence of the XIST gene from late-replicating isodicentric X chromosomes in leukemia. Hum Mol Genet. 1994;3:1053–9.View ArticlePubMedGoogle Scholar
- Brown CJ, Willard HF. The human X inactivation center is not required for maintenance of X inactivation. Nature. 1994;368:154–6.View ArticlePubMedGoogle Scholar
- Agrelo R, Souabni A, Novatchkova M, Haslinger C, Leeb M, Komnenovic V, et al. SATB1 defines the developmental context for gene silencing by Xist in lymphoma and embryonic cells. Dev Cell. 2009;16:507–16.PubMed CentralView ArticlePubMedGoogle Scholar
- Nechanitzky R, Davila A, Savarese F, Fietze S, Grosschedl R. Satb1 and Satb2 are dispensable for X chromosome inactivation in mice. Dev Cell. 2012;23:866–71.View ArticlePubMedGoogle Scholar
- Jeon Y, Lee JT. YY1 tethers Xist RNA to the inactive X nucleation center. Cell. 2011;146:119–33.PubMed CentralView ArticlePubMedGoogle Scholar
- Chow JC, Hall LL, Baldry SE, Thorogood NP, Lawrence JB, Brown CJ. Inducible XIST-dependent X-chromosome inactivation in human somatic cells is reversible. Proc Natl Acad Sci U S A. 2007;104:10104–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang LF, Huynh KD, Lee JT. Perinucleolar targeting of the inactive X during S phase: evidence for a role in the maintenance of silencing. Cell. 2007;129:693–706.View ArticlePubMedGoogle Scholar
- Csankovszki G, Nagy A, Jaenisch R. Synergism of Xist RNA, DNA methylation, and histone hypoacetylation in maintaining X chromosome inactivation. J Cell Biol. 2001;153:773–83.PubMed CentralView ArticlePubMedGoogle Scholar
- Bhatnagar S, Zhu X, Ou J, Lin L, Chamberlain L, Zhu LJ, et al. Genetic and pharmacological reactivation of the mammalian inactive X chromosome. Proc Natl Acad Sci U S A. 2014;111:12591–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Minkovsky A, Sahakyan A, Rankin-Gee E, Bonora G, Patel S, Plath K. The Mbd1-Atf7ip-Setdb1 pathway contributes to the maintenance of X chromosome inactivation. Epigenetics Chromatin. 2014;7:12.PubMed CentralView ArticlePubMedGoogle Scholar
- Yildirim E, Kirby JE, Brown DE, Mercier FE, Sadreyev RI, Scadden DT, et al. Xist RNA is a potent suppressor of hematologic cancer in mice. Cell. 2013;152:727–42.View ArticlePubMedGoogle Scholar
- Silva J, Mak W, Zvetkova I, Appanah R, Nesterova TB, Webster Z, et al. Establishment of histone h3 methylation on the inactive X chromosome requires transient recruitment of Eed-Enx1 polycomb group complexes. Dev Cell. 2003;4:481–95.View ArticlePubMedGoogle Scholar
- Schoeftner S, Sengupta AK, Kubicek S, Mechtler K, Spahn L, Koseki H, et al. Recruitment of PRC1 function at the initiation of X inactivation independent of PRC2 and silencing. EMBO J. 2006;25:3110–22.PubMed CentralView ArticlePubMedGoogle Scholar
- Kalantry S, Magnuson T. The Polycomb Group Protein EED Is Dispensable for the Initiation of Random X-Chromosome Inactivation. PLoS Genet. 2006;2, e66.PubMed CentralView ArticlePubMedGoogle Scholar
- Leeb M, Wutz A. Ring1B is crucial for the regulation of developmental control genes and PRC1 proteins but not X inactivation in embryonic cells. J Cell Biol. 2007;178:219–29.PubMed CentralView ArticlePubMedGoogle Scholar
- Tanasijevic B, Rasmussen TP. X chromosome inactivation and differentiation occur readily in ES cells doubly-deficient for macroH2A1 and macroH2A2. PLoS One. 2011;6, e21512.PubMed CentralView ArticlePubMedGoogle Scholar
- Blewitt ME, Gendrel AV, Pang Z, Sparrow DB, Whitelaw N, Craig JM, et al. SmcHD1, containing a structural-maintenance-of-chromosomes hinge domain, has a critical role in X inactivation. Nat Genet. 2008;40:663–9.View ArticlePubMedGoogle Scholar
- Mould AW, Pang Z, Pakusch M, Tonks ID, Stark M, Carrie D, et al. Smchd1 regulates a subset of autosomal genes subject to monoallelic expression in addition to being critical for X inactivation. Epigenetics Chromatin. 2013;6:19.PubMed CentralView ArticlePubMedGoogle Scholar
- Chaumeil J, Le Baccon P, Wutz A, Heard E. A novel role for Xist RNA in the formation of a repressive nuclear compartment into which genes are recruited when silenced. Genes Dev. 2006;20:2223–37.PubMed CentralView ArticlePubMedGoogle Scholar
- Thorogood NP, Brown CJ. Active chromatin marks are retained on X chromosomes lacking gene or repeat silencing despite XIST/Xist expression in somatic cell hybrids. PLoS One. 2010;5, e10787.PubMed CentralView ArticlePubMedGoogle Scholar
- Clemson CM, Hall LL, Byron M, McNeil J, Lawrence JB. The X chromosome is organized into a gene-rich outer rim and an internal core containing silenced nongenic sequences. Proc Natl Acad Sci U S A. 2006;103:7688–93.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang C, Chapman AG, Kelsey AD, Minks J, Cotton AM, Brown CJ. X-chromosome inactivation: molecular mechanisms from the human perspective. Hum Genet. 2011;130:175–85.View ArticlePubMedGoogle Scholar
- Duthie SM, Nesterova TB, Formstone EJ, Keohane AM, Turner BM, Zakian SM, et al. Xist RNA exhibits a banded localization on the inactive X chromosome and is excluded from autosomal material in cis. Hum Mol Genet. 1999;8:195–204.View ArticlePubMedGoogle Scholar
- Lee JT, Strauss WM, Dausman JA, Jaenisch R. A 450 kb transgene displays properties of the mammalian X-inactivation center. Cell. 1996;86:83–94.View ArticlePubMedGoogle Scholar
- Jiang J, Jing Y, Cost GJ, Chiang JC, Kolpa HJ, Cotton AM, et al. Translating dosage compensation to trisomy 21. Nature. 2013;500:296–300.View ArticlePubMedGoogle Scholar
- Minks J, Baldry SE, Yang C, Cotton AM, Brown CJ. XIST-induced silencing of flanking genes is achieved by additive action of repeat a monomers in human somatic cells. Epigenetics Chromatin. 2013;6:23.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan C, Boyd DD. Histone H3 acetylation and H3 K4 methylation define distinct chromatin regions permissive for transgene expression. Mol Cell Biol. 2006;26:6357–71.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall LL, Carone DM, Gomez AV, Kolpa HJ, Byron M, Mehta N, et al. Stable C0T-1 repeat RNA is abundant and is associated with euchromatic interphase chromosomes. Cell. 2014;156:907–19.PubMed CentralView ArticlePubMedGoogle Scholar
- Engreitz JM, Pandya-Jones A, McDonel P, Shishkin A, Sirokman K, Surka C, et al. The Xist lncRNA exploits three-dimensional genome architecture to spread across the X chromosome. Science. 2013;341:1237973.PubMed CentralView ArticlePubMedGoogle Scholar
- Simon MD, Pinter SF, Fang R, Sarma K, Rutenberg-Schoenberg M, Bowman SK, et al. High-resolution Xist binding maps reveal two-step spreading during X-chromosome inactivation. Nature. 2013;504:465–9.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Nora EP, Dekker J, Heard E. Segmental folding of chromosomes: a basis for structural and regulatory chromosomal neighborhoods? BioEssays. 2013;35:818–28.PubMed CentralView ArticlePubMedGoogle Scholar
- Lyon MF. X-chromosome inactivation: a repeat hypothesis. Cytogenet Cell Genet. 1998;80:133–7.View ArticlePubMedGoogle Scholar
- Cotton AM, Chen CY, 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Bala Tannan N, Brahmachary M, Garg P, Borel C, Alnefaie R, Watson CT, et al. DNA methylation profiling in X;autosome translocations supports a role for L1 repeats in the spread of X chromosome inactivation. Hum Mol Genet. 2014;23:1224–36.PubMed CentralView ArticlePubMedGoogle Scholar
- Gendrel AV, Heard E. Noncoding RNAs and epigenetic mechanisms during X-chromosome inactivation. Annu Rev Cell Dev Biol. 2014;30:561–80.View ArticlePubMedGoogle Scholar
- Gendrel AV, Apedaile A, Coker H, Termanis A, Zvetkova I, Godwin J, et al. Smchd1-dependent and -independent pathways determine developmental dynamics of CpG island methylation on the inactive X chromosome. Dev Cell. 2012;23:265–79.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- UCSC Genome Bioinformatics. Available at: http://genome.ucsc.edu/cgi-bin/hgBlat.
- Xie W, Schultz Matthew D, Lister R, Hou Z, Rajagopal N, Ray P, et al. Epigenomic analysis of multilineage differentiation of human embryonic stem cells. Cell. 2013;153:1134–48.PubMed CentralView ArticlePubMedGoogle Scholar
- Hawkins RD, Hon GC, Lee LK, Ngo Q, Lister R, Pelizzola M, et al. Distinct epigenomic landscapes of pluripotent and lineage-committed human cells. Cell Stem Cell. 2010;6:479–91.PubMed CentralView ArticlePubMedGoogle Scholar
- Goecks J, Nekrutenko A, Taylor J. The Galaxy Team. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11:R86.PubMed CentralView ArticlePubMedGoogle Scholar
- Blankenberg D, Von Kuster G, Coraor N, Ananda G, Lazarus R, Mangan M, et al. Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol. 2010; Chapter 19:Unit 19.10.1-21.Google Scholar
- Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, et al. Galaxy: a platform for interactive large-scale genome analysis. Genome Res. 2005;15:1451–5.PubMed CentralView ArticlePubMedGoogle Scholar