Search for: All records

Abstract BackgroundX chromosome inactivation (XCI) is a female-specific process in which one X chromosome is silenced to balance X-linked gene expression between the sexes. XCI is initiated in early development by upregulation of the lncRNAXiston the future inactive X (Xi). A subset of X-linked genes escape silencing and thus have higher expression in females, suggesting female-specific functions. One of these genes is the highly conserved geneKdm6a, which encodes a histone demethylase that removes methyl groups at H3K27 to facilitate gene expression.KDM6Amutations have been implicated in congenital disorders such as Kabuki Syndrome, as well as in sex differences in development and cancer. MethodsKdm6awas knocked out (KO) using CRISPR/Cas9 gene editing in hybrid female mouse embryonic stem (ES) cells derived either from a 129 × Mus castaneus(cast) cross or a BL6 xcastcross. In one of the lines a transcriptional stop signal inserted inTsixresults in completely skewed X silencing upon differentiation. The effects of both homozygous and heterozygousKdm6aKO onXistexpression during the onset of XCI were measured by RT-PCR and RNA-FISH. Changes in gene expression and in H3K27me3 enrichment were investigated using allele-specific RNA-seq and Cut&Run, respectively. KDM6A binding to theXistgene was characterized by Cut&Run. ResultsWe observed impaired upregulation ofXistand reduced coating of the Xi during early stages of differentiation inKdm6aKO cells, both homozygous and heterozygous, suggesting a threshold effect of KDM6A. This was associated with aberrant overexpression of genes from the Xi after differentiation, indicating loss of X inactivation potency. Consistent with KDM6A having a direct role inXistregulation, we found that the histone demethylase binds to theXistpromoter and KO cells show an increase in H3K27me3 atXist, consistent with reduced expression. ConclusionsThese results reveal a novel female-specific role for the X-linked histone demethylase, KDM6A in the initiation of XCI through histone demethylase-dependent activation ofXistduring early differentiation. Plain language summaryX chromosome inactivation is a female-specific mechanism that evolved to balance sex-linked gene dosage between females (XX) and males (XY) by silencing one X chromosome in females. X inactivation begins with the upregulation of the long noncoding RNAXiston the future inactive X chromosome. While most genes become silenced on the inactive X chromosome some genes escape inactivation and thus have higher expression in females compared to males, suggesting that escape genes may have female-specific functions. One such gene encodes the histone demethylase KDM6A which function is to turn on gene expression by removing repressive histone modifications. In this study, we investigated the role of KDM6A in the regulation ofXistexpression during the onset of X inactivation. We found that KDM6A binds to theXistgene to remove repressive histone marks and facilitate its expression in early development. Indeed, depletion of KDM6A prevents upregulation ofXistdue to abnormal persistence of repressive histone modifications. In turn, this results in aberrant overexpression of genes from the inactive X chromosome. Our findings point to a novel mechanism ofXistregulation during the initiation of X inactivation, which may lead to new avenues of treatment to alleviate congenital disorders such as Kabuki syndrome and sex-biased immune disorders where X-linked gene dosage is perturbed. 

Abstract SummaryHere, we presented the scHiCDiff software tool that provides both nonparametric tests and parametirc models to detect differential chromatin interactions (DCIs) from single-cell Hi-C data. We thoroughly evaluated the scHiCDiff methods on both simulated and real data. Our results demonstrated that scHiCDiff, especially the zero-inflated negative binomial model option, can effectively detect reliable and consistent single-cell DCIs between two conditions, thereby facilitating the study of cell type-specific variations of chromatin structures at the single-cell level. Availability and implementationscHiCDiff is implemented in R and freely available at GitHub (https://github.com/wmalab/scHiCDiff). 

Abstract SummaryHiCube is a lightweight web application for interactive visualization and exploration of diverse types of genomics data at multiscale resolutions. Especially, HiCube displays synchronized views of Hi-C contact maps and 3D genome structures with user-friendly annotation and configuration tools, thereby facilitating the study of 3D genome organization and function. Availability and implementationHiCube is implemented in Javascript and can be installed via NPM. The source code is freely available at GitHub (https://github.com/wmalab/HiCube). 

Recent advances in high-throughput chromosome conformation capture (Hi-C) techniques have allowed us to map genome-wide chromatin interactions and uncover higher-order chromatin structures, thereby shedding light on the principles of genome architecture and functions. However, statistical methods for detecting changes in large-scale chromatin organization such as topologically associating domains (TADs) are still lacking. Here, we proposed a new statistical method, DiffGR, for detecting differentially interacting genomic regions at the TAD level between Hi-C contact maps. We utilized the stratum-adjusted correlation coefficient to measure similarity of local TAD regions. We then developed a nonparametric approach to identify statistically significant changes of genomic interacting regions. Through simulation studies, we demonstrated that DiffGR can robustly and effectively discover differential genomic regions under various conditions. Furthermore, we successfully revealed cell type-specific changes in genomic interacting regions in both human and mouse Hi-C datasets, and illustrated that DiffGR yielded consistent and advantageous results compared with state-of-the-art differential TAD detection methods. The DiffGR R package is published under the GNU General Public License (GPL) ≥ 2 license and is publicly available at https://github.com/wmalab/DiffGR. 

Abstract Background KDM6A is a demethylase encoded by a gene with female-biased expression due to escape from X inactivation. Its main role is to facilitate gene expression through removal of the repressive H3K27me3 mark, with evidence of some additional histone demethylase-independent functions. KDM6A mutations have been implicated in congenital disorders such as Kabuki Syndrome, as well as in sex differences in cancer. Methods Kdm6a was knocked out using CRISPR/Cas9 gene editing in F1 male and female mouse embryonic stem cells (ES) derived from reciprocal crosses between C57BL6 x Mus castaneus . Diploid and allelic RNA-seq analyses were done to compare gene expression between wild-type and Kdm6a knockout (KO) clones. The effects of Kdm6a KO on sex-biased gene expression were investigated by comparing gene expression between male and female ES cells. Changes in H3K27me3 enrichment and chromatin accessibility at promoter regions of genes with expression changes were characterized by ChIP-seq and ATAC-seq followed by diploid and allelic analyses. Results We report that Kdm6a KO in male and female embryonic stem (ES) cells derived from F1 hybrid mice cause extensive gene dysregulation, disruption of sex biases, and specific parental allele effects. Among the dysregulated genes are candidate genes that may explain abnormal developmental features of Kabuki syndrome caused by KDM6A mutations in human. Strikingly, Kdm6a knockouts result in a decrease in sex-biased expression and in preferential downregulation of the maternal alleles of a number of genes. Most promoters of dysregulated genes show concordant epigenetic changes including gain of H3K27me3 and loss of chromatin accessibility, but there was less concordance when considering allelic changes. Conclusions Our study reveals new sex-related roles of KDM6A in the regulation of developmental genes, the maintenance of sex-biased gene expression, and the differential expression of parental alleles. 

Abstract The recently developed Hi-C technique has been widely applied to map genome-wide chromatin interactions. However, current methods for analyzing diploid Hi-C data cannot fully distinguish between homologous chromosomes. Consequently, the existing diploid Hi-C analyses are based on sparse and inaccurate allele-specific contact matrices, which might lead to incorrect modeling of diploid genome architecture. Here we present ASHIC, a hierarchical Bayesian framework to model allele-specific chromatin organizations in diploid genomes. We developed two models under the Bayesian framework: the Poisson-multinomial (ASHIC-PM) model and the zero-inflated Poisson-multinomial (ASHIC-ZIPM) model. The proposed ASHIC methods impute allele-specific contact maps from diploid Hi-C data and simultaneously infer allelic 3D structures. Through simulation studies, we demonstrated that ASHIC methods outperformed existing approaches, especially under low coverage and low SNP density conditions. Additionally, in the analyses of diploid Hi-C datasets in mouse and human, our ASHIC-ZIPM method produced fine-resolution diploid chromatin maps and 3D structures and provided insights into the allelic chromatin organizations and functions. To summarize, our work provides a statistically rigorous framework for investigating fine-scale allele-specific chromatin conformations. The ASHIC software is publicly available at https://github.com/wmalab/ASHIC.