Abnormal X chromosome inactivation and sex-specific gene dysregulation after ablation of FBXL10
© The Author(s) 2016
Received: 8 March 2016
Accepted: 10 May 2016
Published: 31 May 2016
Almost all CpG-rich promoters in the mammalian genome are bound by the multidomain FBXL10 protein (also known as KDM2B, JHDM1B, CXXC2, and NDY1). FBXL10 is expressed as two isoforms: FBXL10-1, a longer form that contains an N-terminal histone demethylase domain with C-terminal F-box, CXXC, PHD, RING, and leucine-rich repeat domains, and FBXL10-2, a shorter form that initiates at an alternative internal exon and which lacks the histone demethylase domain but retains all other annotated domains. Selective deletion of Fbxl10-1 had been reported to produce a low penetrance and variable phenotype; most of the mutant animals were essentially normal. We constructed mutant mouse strains that were either null for Fbxl10-2 but wild type for Fbxl10-1 or null for both Fbxl10-1 and Fbxl10-2.
Deletion of Fbxl10-2 (in a manner that does not perturb expression of Fbxl10-1) produced a phenotype very different from the Fbxl10-1 mutant, with craniofacial abnormalities, neural tube defects, and increased lethality, especially in females. Mutants that lacked both FBXL10-1 and FBXL10-2 showed embryonic lethality and even more extreme sexual dimorphism, with more severe gene dysregulation in mutant female embryos. X-linked genes were most severely dysregulated, and there was marked overexpression of Xist in mutant females although genes that encode factors that bind to Xist RNA were globally downregulated in mutant female as compared to male embryos.
FBXL10 is the first factor shown to be required both for the normal expression and function of the Xist gene and for normal expression of proteins that associate with Xist RNA; it is proposed that FBXL10 coordinates the expression of Xist RNA with proteins that associate with this RNA. The function of FBXL10 is largely independent of the histone demethylase activity of the long form of the protein.
FBXL10 is a multidomain chromosomal protein that is localized to essentially all CpG-rich promoters via a CXXC domain that binds to unmethylated CpG dinucleotides [1, 3]. FBXL10 is present in a variant form of Polycomb Repressive Complex 1 (PRC1), but most FBXL10 is deposited at promoters that are not associated with Polycomb factors. Many biological functions have been ascribed to FBXL10; among these are the regulation of nucleolar function , organization of Polycomb Repressive Complexes [3, 4], self-renewal of stem cells , and regulation of histone H2A ubiquitylation . However, these studies were conducted prior to the demonstration that FBXL10 is bound to all CpG-rich promoters and none utilized mouse strains null for Fbxl10. FBXL10 has also been reported to be an oncogene and a tumor suppressor in leukemia [7, 8]. The biological activities of FBXL10 are usually attributed to the JmjC histone demethylase domain, although FBXL10 also contains CXXC, PHD, RING, F-box, and leucine-rich repeat (LRR) domains. It was recently reported that FBXL10 is required for the protection of sequences bound by Polycomb Repressive Complexes from de novo methylation, while promoters bound by FBXL10 but not by Polycomb Repressive Complexes did not undergo de novo methylation after removal of FBXL10 .
Ablation of the long form of the protein (the JmjC-containing FBXL10-1) was reported to produce a fully recessive phenotype in which most homozygous mice were of normal phenotype but some displayed coloboma, exencephaly, and rare tail kinks . A single conservative amino acid substitution that affects both FBXL10-1 and FBXL10-2 causes paunch calf syndrome with craniofacial dysmorphia and perinatal lethality in Romagnolo cattle; this mutation is also fully recessive . A deletion that removed the CXXC domain (which is common to both FBXL10-1 and FBXL10-2) was reported to produce increased lethality and developmental abnormalities in heterozygous mice as a result of haploinsufficiency, but this mutant allele produced large amounts of an abnormal protein that could have exerted neomorphic effects . The semidominant character of this allele contrasts with the fully recessive nature of all other mutant alleles of Fbxl10 (including alleles described here).
We report here that a mutant allele that deleted Fbxl10-2 in a manner that did not affect expression of Fbxl10-1 was fully recessive and produced a phenotype that differed markedly from the Fbxl10-1 deletion. The Fbxl10-2 mutation caused craniofacial abnormalities, with high rates of perinatal lethality and a strain-dependent cleft palate and eyes-open-at-birth phenotype. We observed that female homozygotes had a more severe phenotype and that very few survived to adulthood. We further investigated the sexual dimorphism in embryos that lacked both FBXL10-1 and FBXL10-2 and found dysregulation of many genes in female embryos, with pronounced dysregulation of X-linked genes and strong overexpression of Xist specifically in female embryos, with a net decrease in X-linked gene expression that was not detectable in male embryos. Male embryos showed greatly decreased Tsix expression but did not detectably reactivate Xist. The data indicate that removal of FBXL10 results in widespread dysregulation of gene expression, with much larger effects in female embryos and even greater dysregulation of X-linked genes. FBXL10 is a novel factor that is required for normal expression of Xist and Tsix and for the coordinated expression of Xist RNA and the protein factors that associate with this RNA.
Two biological replicates for each genotype of both sexes were analyzed. Total RNA was extracted from single embryos dissected at embryonic day E9.5. The sex of the embryos was determined by PCR detection of the Sry gene using genomic DNA extracted from the yolk sac. RNA-seq reads were mapped to the mouse genome (mm9) using TopHat (v1.4.1; ). Mapped reads were assigned to RefSeq annotated genes using htseq-count (version 0.5.3p9; http://www-huber.embl.de/users/anders/HTSeq/doc/count.html). Differential expression was assessed by edgeR (http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2796818/) with TMM normalization. Log2 fold-changes were computed after applying a floor of 1 cpm to ensure that genes with very low expression but large fold-changes did not affect the analysis. Pearson’s correlation coefficient was used to compare the global variation in expression levels in +/+ and T/T mice. Distributions of log2 fold expression changes between autosomal and ChrX genes were compared using the Mann–Whitney U test.
Gene set enrichment analysis was performed as in  using the log2 fold-change to rank genes. A custom gene signature from 81 Xist interacting factors  was used for the analysis. The density of Xist bound factors was computed as the number of genes in a 2000 gene window normalized by the window size as the window slides across the ranked gene list. P values are those reported from GSEA for the enrichment of Xist bound factors.
FBXL10 ChIP-Seq data from mouse ES cells are from GSE41316 .
The RNA-seq data were deposited at the Gene Expression Omnibus under the Accession number GSE76552.
Selective deletion of FBXL10-2: sexually dimorphic phenotype distinct from FBXL10-1 mutation
As shown in Fig. 1a, b, FBXL10-2 differs from FBXL10-1 by the absence of an N-terminal domain that contains a JmjC histone demethylase domain. As mentioned previously, a deletion that selectively removed FBXL10-1 was reported to have a weakly penetrant recessive phenotype characterized by low-frequency exencephaly, coloboma, and tail kinks, although most homozygous mutants were without visible phenotype . We selectively removed FBXL10-2 by introduction of loxP sites followed by Cre-mediated excision of the alternative internal promoter and first exon that drive expression of FBXL10-2 (Fig. 1b, Additional file 1: Fig. S1). As shown in Fig. 1b, c, the deletion removed all detectable FBXL10-2 protein without measurable effect on FBXL10-1.
The FBXL10-1 and FBXL10-2 deletions had markedly different phenotypes. The FBXL10-2 deletion did not show exencephaly, coloboma, or kinked tail as did the FBXL10-1 mutants. In contrast, surviving animals lacking FBXL10-2 were runted and displayed craniofacial abnormalities (Additional file 1: Fig. S2a), while mutants that died prenatally showed craniofacial and neural tube abnormalities (Fig. 1d). When backcrossed onto the FVB genetic background, the Fbxl10 ∆−2/∆−2 mutation additionally caused cleft pallet and open eyelids at birth (Fig. 1d, Additional file 1: Fig. S2b). We performed all gene expression experiments on the C57Bl6 strain background.
Abnormal development and gene expression in Fbxl10-null embryos
Additional file 3: Table S2 shows differential gene expression in male versus female Fbxl10 null mutants, and Fig. 3c shows examples of genes that showed sex-specific dysregulation. The imprinted Meg3 and Gtl2 genes  were notably upregulated in mutant XX embryos relative to XY embryos, although other nearby genes such as Rian were expressed at normal levels (Fig. 3c). Expression of Meg3 and Gtl2 was increased by a factor of greater than 2, which suggests that activation of the normally silent imprinted alleles was unlikely to be responsible for the increased expression. Abnormalities of gene expression were widespread and more apparent in female than in male embryos.
Abnormal expression of X-linked genes in the absence of FBXL10
Normal reactivation of the paternally imprinted X chromosome in Fbxl10 null mutants
The observed global downregulation of the active X chromosome could arise from the failure of the paternally imprinted X chromosome to be reactivated. A defective counting in cells that have a paternally imprinted X chromosome would cause the inactivation of two X chromosomes in half the cells (the cells that randomly inactivate the X chromosome of maternal origin). We reasoned that a failure to reactivate the paternally imprinted X chromosome would lead to a skewed X inactivation pattern. We tested for potential allelic imbalance of two X-linked transcripts using single nucleotide polymorphisms (SNPs) to trace the parental origin of the transcripts (Additional file 1: Fig. S4). Haus7 and Sh3bgr1 are X-linked genes with reduced expression in mutant females (Additional file 1: Fig. S4a, c). Our SNP analysis showed that transcripts from each parental allele were equally represented in female mutant embryos (Additional file 1: Fig. S4b, d). Therefore, X inactivation is not skewed in female embryos that lack FBXL10, which indicates that choice in somatic X inactivation is normal in the absence of FBXL10.
Abnormal expression of Xist RNA and Xist RNA-binding factors in Fbxl10-null female embryos
Individual ablations of FBXL10-1 and FBXL-2 produced non-overlapping phenotypes. The long form deletion produced low-penetrance exencephaly, colobomo, and rare tail kinks with low mortality , while deletion of the short form caused high-penetrance anterior–posterior craniofacial contraction with strain-dependent cleft palate and an eyes-open-at-birth phenotype, with high perinatal lethality. This is an unusual case in which selective deletion of two isoforms of the same gene gives rise to highly dissimilar phenotypes. The data also indicate that histone demethylation via the JmjC domain of FBXL10-1 is largely dispensable for the biological activity of the protein. There are other cases in which the histone modifying domains of chromosomal proteins are largely dispensable; the histone demethylase activity of KDM6B is dispensable for the biological function of the protein , as is the E3 ubiquitin ligase activity of RING1B/RNF2 .
We observed sexual dimorphism in the form of more severe developmental defects in female mice that lacked FBXL10-2; we then tested for sexually dimorphic phenotypes in mouse embryos that lacked both the FBXL10-1 and FBXL10-2. The Fbxl0 null mutation was found to cause much more severe developmental abnormalities and earlier death in female embryos; at E10.5, the most advanced female embryos resembled the least advanced male embryos. Analysis of gene expression abnormalities by RNA-seq on single E9.5 embryos showed much greater dysregulation of gene expression in female embryos. This was unexpected given that FBXL10 is localized to nearly all CpG-rich promoters [1, 3], which account for ~76 % of all promoters.
The sexually dimorphic phenotypes were apparent well prior to the onset of sexual differentiation, which suggested that abnormalities of sex chromosome function might be involved. Inspection of RNA-seq data showed that Xist was markedly overexpressed in Fbxl10 null mutant female embryos, while Tsix was reduced in mutant male embryos. Deficiencies of the MYST1-(MOF)-containing MSL and NSL complexes have been reported to cause downregulation of Tsix and biallelic expression of Xist with Xist RNA coating of both X chromosomes , but examination of RNA-seq data showed that all of the ~20 components of the MSL and NSL complexes were expressed normally in male and female Fbxl10 mutants (data not shown). Deletion of Tsix has also been reported to cause biallelic expression of Xist in female cells , but the marked downregulation of Tsix in male Fbxl10 mutant embryos or ES cells did not result in detectable expression of Xist from the male X chromosome, and expression of Tsix was not detectable in mutant female embryos. Furthermore, FBXL10 is not detectably enriched at the promoters of the Xist or Tsix genes .
A recent study identified a set of ~81 proteins that interact with Xist RNA . We found that many of these factors were expressed at reduced levels in Fbxl10-null female embryos but were expressed at normal levels in mutant male embryos. Hnrnpu, which is required for the association of the Xist complex with the inactive X chromosome , is reduced by more than twofold, as is the histone ubiquitin transferase RING1B/RNF2. The pronounced imbalance created by excess Xist RNA combined to the reduced concentrations of Xist binding factors will create a state in which highly heterogeneous and ectopic Xist–protein complexes in abnormally large amounts can disrupt patterns of gene expression. Genes on the X chromosome encode numerous transcription regulators; these include Tbx22, Dax1, Sox3, and Pou3F4. We suggest that ectopic Xist silencing complexes that result from the overexpression of Xist coupled to the underexpression of Xist binding factors leads to a pervasive dysregulation of gene expression driven by repression of genes on the active X chromosome. This gene dysregulation is independent of the methylation abnormalities induced by the deprivation of FBXL10, which affects only genes that are bound by Polycomb Repressive Complexes 1 and 2 .
We show that removal of FBXL10 results in overexpression of Xist in females, reduced expression of proteins reported to complex with Xist RNA, and greater perturbations of global gene expression profiles in females as compared to males; the data indicate that FBXL10 is a factor required for the coordinated expression of Xist, Tsix, and Xist RNA-associated factors. In the absence of FBXL10 abnormal X inactivation leads to a global reduction of X-linked gene expression and dysregulated expression of autosomal genes specifically in female embryos.
MB, JRE, and THB designed the study, MB carried out the experimental procedures, JRE performed computational analysis, and MB, JRE, and THB wrote the paper. All authors read and approved the final manuscript.
This study was supported by Grants from the NIH to J. R. E. and T. H. B. and by a Grant from the US Department of Defense to J. R. E. We thank V. Bardwell (University of Minnesota) for her gift of antiserum to FBXL10.
The authors declare that they have no competing interests.
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.
- Farcas AM, Blackledge NP, Sudbery I, Long HK, McGouran JF, Rose NR, et al. KDM2B links the polycomb repressive complex 1 (PRC1) to recognition of CpG islands. eLife. 2012;1:e00205.View ArticlePubMedPubMed CentralGoogle Scholar
- He J, Shen L, Wan M, Taranova O, Wu H, Zhang Y. Kdm2b maintains murine embryonic stem cell status by recruiting PRC1 complex to CpG islands of developmental genes. Nat Cell Biol. 2013;15:373–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Frescas D, Guardavaccaro D, Bassermann F, Koyama-Nasu R, Pagano M. JHDM1B/FBXL10 is a nucleolar protein that represses transcription of ribosomal RNA genes. Nature. 2007;450(7167):309–13.View ArticlePubMedGoogle Scholar
- Blackledge NP, Farcas AM, Kondo T, King HW, McGouran JF, Hanssen LL, et al. Variant PRC1 complex-dependent H2A ubiquitylation drives PRC2 recruitment and polycomb domain formation. Cell. 2014;157(6):1445–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Kottakis F, Foltopoulou P, Sanidas I, Keller P, Wronski A, Dake BT, et al. NDY1/KDM2B functions as a master regulator of polycomb complexes and controls self-renewal of breast cancer stem cells. Cancer Res. 2014;74(14):3935–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu X, Johansen JV, Helin K. Fbxl10/Kdm2b recruits polycomb repressive complex 1 to CpG Islands and regulates H2A ubiquitylation. Mol Cell. 2014;49(6):1134–46.View ArticleGoogle Scholar
- Pasqualucci L, Trifonov V, Fabbri G, Ma J, Rossi D, Chiarenza A, et al. Analysis of the coding genome of diffuse large B-cell lymphoma. Nat Genet. 2011;43(9):830–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Suzuki T, Minehata K-I, Akagi K, Jenkins NA, Copeland NG. Tumor suppressor gene identification using retroviral insertional mutagenesis in Blm-deficient mice. EMBO J. 2006;25(14):3422–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Boulard M, Edwards JR, Bestor TH. FBXL10 protects Polycomb-bound genes from hypermethylation. Nat Genet. 2015;47(5):479–85.View ArticlePubMedGoogle Scholar
- Fukuda T, Tokunaga A, Sakamoto R, Yoshida N. Fbxl10/Kdm2b deficiency accelerates neural progenitor cell death and leads to exencephaly. Mol Cell Neurosci. 2011;46(3):614–24.View ArticlePubMedGoogle Scholar
- Testoni S, Bartolone E, Rossi M, Patrignani A, Bruggmann R, Lichtner P, et al. KDM2B is implicated in bovine lethal multi-organic developmental dysplasia. PLoS ONE. 2012;7(9):e45634.View ArticlePubMedPubMed CentralGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-seq. Bioinformatics. 2009;25(9):1105–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Chu C, Zhang QC, da Rocha ST, Flynn RA, Bharadwaj M, Calabrese JM, et al. Systematic discovery of Xist RNA binding proteins. Cell. 2015;161(2):404–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Takada S, Tevendale M, Baker J, Georgiades P, Campbell E, Freeman T, Johnson MH, Paulsen M, Ferguson-Smith AC. Delta-like and gtl2 are reciprocally expressed, differentially methylated linked imprinted genes on mouse chromosome 12. Curr Biol. 2000;10(18):1135–8.View ArticlePubMedGoogle Scholar
- Shpargel KB, Sengoku T, Yokoyama S, Magnuson T. UTX and UTY demonstrate histone demethylase-independent function in mouse embryonic development. PLoS Genet. 2012;8(9):e1002964.View ArticlePubMedPubMed CentralGoogle Scholar
- Illingworth RS, Moffat M, Mann AR, Read D, Hunter CJ, Pradeepa MM, Adams IR, Bickmore WA. The E3 ubiquitin ligase activity of RING1B is not essential for early mouse development. Genes Dev. 2015;29(18):1897–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Chelmicki T, Dündar F, Turley MJ, Khanam T, Aktas T, Ramírez F, et al. MOF-associated complexes ensure stem cell identity and Xist repression. Elife. 2014;3:e02024.View ArticlePubMedPubMed CentralGoogle Scholar
- Sado T, Wang Z, Sasaki H, Li E. Regulation of imprinted X-chromosome inactivation in mice by Tsix. Development. 2001;128(8):1275–86.PubMedGoogle Scholar
- Hacisuleyman E, Goff LA, Trapnell C, Williams A, Henao-Mejia J, Sun L, et al. Topological organization of multichromosomal regions by the long intergenic noncoding RNA Firre. Nat Struct Mol Biol. 2014;21(2):198–206.View ArticlePubMedPubMed CentralGoogle Scholar