A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility.

We applied a combinatorial indexing assay, sci-ATAC-seq, to profile genome-wide chromatin accessibility in ∼100,000 single cells from 13 adult mouse tissues. We identify 85 distinct patterns of chromatin accessibility, most of which can be assigned to cell types, and ∼400,000 differentially accessible elements. We use these data to link regulatory elements to their target genes, to define the transcription factor grammar specifying each cell type, and to discover in vivo correlates of heterogeneity in accessibility within cell types. We develop a technique for mapping single cell gene expression data to single-cell chromatin accessibility data, facilitating the comparison of atlases. By intersecting mouse chromatin accessibility with human genome-wide association summary statistics, we identify cell-type-specific enrichments of the heritability signal for hundreds of complex traits. These data define the in vivo landscape of the regulatory genome for common mammalian cell types at single-cell resolution.

High-Throughput Single-Cell Sequencing with Linear Amplification.

Conventional methods for single-cell genome sequencing are limited with respect to uniformity and throughput. Here, we describe sci-L3, a single-cell sequencing method that combines combinatorial indexing (sci-) and linear (L) amplification. The sci-L3 method adopts a 3-level (3) indexing scheme that minimizes amplification biases while enabling exponential gains in throughput. We demonstrate the generalizability of sci-L3 with proof-of-concept demonstrations of single-cell whole-genome sequencing (sci-L3-WGS), targeted sequencing (sci-L3-target-seq), and a co-assay of the genome and transcriptome (sci-L3-RNA/DNA). We apply sci-L3-WGS to profile the genomes of >10,000 sperm and sperm precursors from F1 hybrid mice, mapping 86,786 crossovers and characterizing rare chromosome mis-segregation events in meiosis, including instances of whole-genome equational chromosome segregation. We anticipate that sci-L3 assays can be applied to fully characterize recombination landscapes, to couple CRISPR perturbations and measurements of genome stability, and to other goals requiring high-throughput, high-coverage single-cell sequencing.

X chromosome regulation: diverse patterns in development, tissues and disease.

Genes on the mammalian X chromosome are present in one copy in males and two copies in females. The complex mechanisms that regulate the X chromosome lead to evolutionary and physiological variability in gene expression between species, the sexes, individuals, developmental stages, tissues and cell types. In early development, delayed and incomplete X chromosome inactivation (XCI) in some species causes variability in gene expression. Additional diversity stems from escape from XCI and from mosaicism or XCI skewing in females. This causes sex-specific differences that manifest as differential gene expression and associated phenotypes. Furthermore, the complexity and diversity of X dosage regulation affect the severity of diseases caused by X-linked mutations.

Allele-specific non-CG DNA methylation marks domains of active chromatin in female mouse brain.

DNA methylation at gene promoters in a CG context is associated with transcriptional repression, including at genes silenced on the inactive X chromosome in females. Non-CG methylation (mCH) is a distinct feature of the neuronal epigenome that is differentially distributed between males and females on the X chromosome. However, little is known about differences in mCH on the active (Xa) and inactive (Xi) X chromosomes because stochastic X-chromosome inactivation (XCI) confounds allele-specific epigenomic profiling. We used whole-genome bisulfite sequencing in a mouse model with nonrandom XCI to examine allele-specific DNA methylation in frontal cortex. Xi was largely devoid of mCH, whereas Xa contained abundant mCH similar to the male X chromosome and the autosomes. In contrast to the repressive association of DNA methylation at CG dinucleotides (mCG), mCH accumulates on Xi in domains with transcriptional activity, including the bodies of most genes that escape XCI and at the X-inactivation center, validating this epigenetic mark as a signature of transcriptional activity. Escape genes showing CH hypermethylation were the only genes with CG-hypomethylated promoters on Xi, a well-known mark of active transcription. Finally, we found extensive allele-specific mCH and mCG at autosomal imprinted regions, some with a negative correlation between methylation in the two contexts, further supporting their distinct functions. Our findings show that neuronal mCH functions independently of mCG and is a highly dynamic epigenomic correlate of allele-specific gene regulation.

Escape from X inactivation varies in mouse tissues.

X chromosome inactivation (XCI) silences most genes on one X chromosome in female mammals, but some genes escape XCI. To identify escape genes in vivo and to explore molecular mechanisms that regulate this process we analyzed the allele-specific expression and chromatin structure of X-linked genes in mouse tissues and cells with skewed XCI and distinguishable alleles based on single nucleotide polymorphisms. Using a binomial model to assess allelic expression, we demonstrate a continuum between complete silencing and expression from the inactive X (Xi). The validity of the RNA-seq approach was verified using RT-PCR with species-specific primers or Sanger sequencing. Both common escape genes and genes with significant differences in XCI status between tissues were identified. Such genes may be candidates for tissue-specific sex differences. Overall, few genes (3-7%) escape XCI in any of the mouse tissues examined, suggesting stringent silencing and escape controls. In contrast, an in vitro system represented by the embryonic-kidney-derived Patski cell line showed a higher density of escape genes (21%), representing both kidney-specific escape genes and cell-line specific escape genes. Allele-specific RNA polymerase II occupancy and DNase I hypersensitivity at the promoter of genes on the Xi correlated well with levels of escape, consistent with an open chromatin structure at escape genes. Allele-specific CTCF binding on the Xi clustered at escape genes and was denser in brain compared to the Patski cell line, possibly contributing to a more compartmentalized structure of the Xi and fewer escape genes in brain compared to the cell line where larger domains of escape were observed.

Female bias in Rhox6 and 9 regulation by the histone demethylase KDM6A.

The Rhox cluster on the mouse X chromosome contains reproduction-related homeobox genes expressed in a sexually dimorphic manner. We report that two members of the Rhox cluster, Rhox6 and 9, are regulated by de-methylation of histone H3 at lysine 27 by KDM6A, a histone demethylase with female-biased expression. Consistent with other homeobox genes, Rhox6 and 9 are in bivalent domains prior to embryonic stem cell differentiation and thus poised for activation. In female mouse ES cells, KDM6A is specifically recruited to Rhox6 and 9 for gene activation, a process inhibited by Kdm6a knockdown in a dose-dependent manner. In contrast, KDM6A occupancy at Rhox6 and 9 is low in male ES cells and knockdown has no effect on expression. In mouse ovary where Rhox6 and 9 remain highly expressed, KDM6A occupancy strongly correlates with expression. Our study implicates Kdm6a, a gene that escapes X inactivation, in the regulation of genes important in reproduction, suggesting that KDM6A may play a role in the etiology of developmental and reproduction-related effects of X chromosome anomalies.

CTCF-mediated insulation and chromatin environment modulate Car5b escape from X inactivation.

Background: The number and escape levels of genes that escape X chromosome inactivation (XCI) in female somatic cells vary among tissues and cell types, potentially contributing to specific sex differences. Here we investigate the role of CTCF, a master chromatin conformation regulator, in regulating escape from XCI. CTCF binding profiles and epigenetic features were systematically examined at constitutive and facultative escape genes using mouse allelic systems to distinguish the inactive X (Xi) and active X (Xa) chromosomes. Results: We found that escape genes are located inside domains flanked by convergent arrays of CTCF binding sites, consistent with the formation of loops. In addition, strong and divergent CTCF binding sites often located at the boundaries between escape genes and adjacent neighbors subject to XCI would help insulate domains. Facultative escapees show clear differences in CTCF binding dependent on their XCI status in specific cell types/tissues. Concordantly, deletion but not inversion of a CTCF binding site at the boundary between the facultative escape gene Car5b and its silent neighbor Siah1b resulted in loss of Car5b escape. Reduced CTCF binding and enrichment of a repressive mark over Car5b in cells with a boundary deletion indicated loss of looping and insulation. In mutant lines in which either the Xi-specific compact structure or its H3K27me3 enrichment was disrupted, escape genes showed an increase in gene expression and associated active marks, supporting the roles of the 3D Xi structure and heterochromatic marks in constraining levels of escape. Conclusion: Our findings indicate that escape from XCI is modulated both by looping and insulation of chromatin via convergent arrays of CTCF binding sites and by compaction and epigenetic features of the surrounding heterochromatin.

KDM6A facilitates Xist upregulation at the onset of X inactivation.

X 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 lncRNA Xist on 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 gene Kdm6a , which encodes a histone demethylase that removes methyl groups at H3K27 to facilitate gene expression. Here, we investigate the role of KDM6A in the regulation of Xist . We observed impaired upregulation of Xist during early stages of differentiation in hybrid mouse ES cells following CRISPR/Cas9 knockout of Kdm6a . This is associated with reduced Xist RNA coating of the Xi, suggesting diminished XCI potency. Indeed, Kdm6a knockout results in aberrant overexpression of genes from the Xi after differentiation. KDM6A binds to the Xist promoter and knockout cells show an increase in H3K27me3 at Xist . These results indicate that KDM6A plays a role in the initiation of XCI through histone demethylase-dependent activation of Xist during early differentiation.

Sex-biased and parental allele-specific gene regulation by KDM6A.

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.

Genes that escape from X inactivation.

To achieve a balanced gene expression dosage between males (XY) and females (XX), mammals have evolved a compensatory mechanism to randomly inactivate one of the female X chromosomes. Despite this chromosome-wide silencing, a number of genes escape X inactivation: in women about 15% of X-linked genes are bi-allelically expressed and in mice, about 3%. Expression from the inactive X allele varies from a few percent of that from the active allele to near equal expression. While most genes have a stable inactivation pattern, a subset of genes exhibit tissue-specific differences in escape from X inactivation. Escape genes appear to be protected from the repressive chromatin modifications associated with X inactivation. Differences in the identity and distribution of escape genes between species and tissues suggest a role for these genes in the evolution of sex differences in specific phenotypes. The higher expression of escape genes in females than in males implies that they may have female-specific roles and may be responsible for some of the phenotypes observed in X aneuploidy.

X Inactivation and Escape: Epigenetic and Structural Features.

X inactivation represents a complex multi-layer epigenetic mechanism that profoundly modifies chromatin composition and structure of one X chromosome in females. The heterochromatic inactive X chromosome adopts a unique 3D bipartite structure and a location close to the nuclear periphery or the nucleolus. X-linked lncRNA loci and their transcripts play important roles in the recruitment of proteins that catalyze chromatin and DNA modifications for silencing, as well as in the control of chromatin condensation and location of the inactive X chromosome. A subset of genes escapes X inactivation, raising questions about mechanisms that preserve their expression despite being embedded within heterochromatin. Escape gene expression differs between males and females, which can lead to physiological sex differences. We review recent studies that emphasize challenges in understanding the role of lncRNAs in the control of epigenetic modifications, structural features and nuclear positioning of the inactive X chromosome. Second, we highlight new findings about the distribution of genes that escape X inactivation based on single cell studies, and discuss the roles of escape genes in eliciting sex differences in health and disease.

Epigenetic and genetic mechanisms contribute to telomerase inhibition by EGCG.

The ends of human chromosomes are protected from the degradation associated with cell division by 15-20 kb long segments of hexameric repeats of 5'-TTAGGG-3' termed telomeres. In normal cells telomeres lose up to 300 bp of DNA per cell division that ultimately leads to senescence; however, most cancer cells bypass this lifespan restriction through the expression of telomerase. hTERT, the catalytic subunit essential for the proper function of telomerase, has been shown to be expressed in approximately 90% of all cancers. In this study we investigated the hTERT inhibiting effects of (-)-epigallocatechin-3-gallate (EGCG), the major polyphenol found in green tea catechins, in MCF-7 breast cancers cells and HL60 promyelocytic leukemia cells. Exposure to EGCG reduced cellular proliferation and induced apoptosis in both MCF-7 and HL60 cells in vitro, although hTERT mRNA expression was decreased only in MCF-7 cells when treated with EGCG. Furthermore, down-regulation of hTERT gene expression in MCF-7 cells appeared to be largely due to epigenetic alterations. Treatment of MCF-7 cells with EGCG resulted in a time-dependent decrease in hTERT promoter methylation and ablated histone H3 Lys9 acetylation. In conjunction with demethylation, further analysis showed an increase in hTERT repressor E2F-1 binding at the promoter. From these findings, we propose that EGCG is effective in causing cell death in both MCF-7 and HL60 cancer cell lines and may work through different pathways involving both anti-oxidant effects and epigenetic modulation.

Epigenetic regulation of telomerase in retinoid-induced differentiation of human leukemia cells.

Changes in the promoter methylation of hTERT, the gene that encodes telomerase, a ribonucleoprotein responsible for replacing telomeric repeats, have been demonstrated in differentiating cells where hTERT is inhibited, suggesting epigenetic regulation of hTERT. All-trans retinoic acid (ATRA) induces differentiation in human leukemia cells and has had significant clinical success treating promyelocytic leukemia in what is termed 'differentiation therapy'. It is thought that the inhibition of telomerase is a target of retinoids and is closely tied to the differentiated phenotype. This study demonstrates the epigenetic changes associated with ATRA-induced inhibition of telomerase activity, including the hypoacetylation and hypermethylation of the hTERT promoter. Further, we have found changes in the differential expression of the three DNA methyltransferases during ATRA-induced differentiation of HL60 human leukemia cells. These results suggest that alteration of DNA methylation may play a role in the activation of telomerase in cancer cells and that epigenetic mechanisms may represent a target for differentiation therapy mechanisms. We propose that epigenetic changes in the hTERT promoter represent a stable locking mechanism in the retionoid-induced suppression of telomerase activity.

X-Chromosome Inactivation and Escape from X Inactivation in Mouse.

X chromosome inactivation silences one X chromosome in female mammals. However, this silencing is incomplete, and some genes escape X inactivation. We describe methods to determine the chromosome-wide X inactivation status of genes in tissues or cell lines derived from mice using a combination of skewing of X inactivation and allele-specific analyses of gene expression based on RNA-seq.

A method to detect DNA methyltransferase I gene transcription in vitro in aging systems.

Epigenetic alterations of DNA play key roles in determining gene structure and expression. Methylation of the 5-position of cytosine is thought to be the most common modification of the genome in mammals. Studies have generally shown that hypermethylation in gene regulatory regions is associated with inactivation and reduced transcription and that alteration in established methylation patterns during development can affect embryonic viability. Changes in methylation have also been associated with aging and cellular senescence as well as tumorogenesis. DNA methyltransferase 1 (DNMT1) is thought to play an important role in maintaining already established methylation patterns during DNA replication and catalyzes the transfer of a methyl moiety from S-adenosyl-L-methionine (SAM) to the 5-position of cytosines in the CpG dinucleotide. Several studies illustrate changes in activity and transcription of DNMT1 during aging and here we show a comprehensive method of detection of DNMT1 mRNA transcription from senescing cells in culture.

Aging cell culture: methods and observations.

Culturing and subcultivation of normal human diploid fibroblasts have advanced our understanding of the molecular events involved in aging. This progress is leading to the development of therapies that slow or ablate the adverse physiological and pathological changes associated with aging. It has been established that normal human diploid fibroblasts can proliferate in culture for only finite periods of time. Hayflick and Moorhead and others have described numerous types of cells, ranging from fetal to adult, that were incapable of indefinite proliferation. There are many ways to study aging in vitro, and this chapter summarizes some laboratory procedures.

A method to study the expression of DNA methyltransferases in aging systems in vitro.

The methylation of CpG dinucleotides located in key protein binding sites within gene regulatory regions often leads to gene silencing. A mechanism of aging is proposed whereby an accumulation of methylation at gene regulatory sites contributes to cellular senescence. DNA methyltransferases (DNMTs) are enzymes that catalyze the transfer of a methyl moiety from S-adenosyl-L-methionine (SAM) to the cytosine of a CpG dinucleotide and are responsible for establishing and maintaining methylation patterns in the genome. It is important to study not only transcription of the DNMTs, but also their protein expression because studies illustrate that it is possible for the enzymes to undergo posttranslational physical changes in response to stimulation even though gene transcription remains unchanged. Here, we discuss an in vitro method to study protein expression of DNMTs in aging systems.

Telomerase inhibition by retinoids precedes cytodifferentiation of leukemia cells and may contribute to terminal differentiation.

Human promyelocytic leukemia HL60 cells display high telomerase activity, a phenotype related to their immortal status. All-trans retinoic acid (ATRA) is a clinically effective cytodifferentiating agent. To understand the mechanism underlying ATRA-induced cytodifferentiation, we did a kinetic analysis of the role of ATRA in inhibiting telomerase in HL60 cells. Our studies indicate that telomerase inhibition by ATRA occurred relatively early after treatment of HL60 cells due to a rapid decrease in hTERT gene expression. More importantly, however, we found through monitoring the expression of CD11b, a marker for granulocytic differentiation of HL60 cells, that down-regulation of telomerase preceded the differentiation of HL60 cells. These observations suggest that the hTERT gene may be a primary target of ATRA regulation of cellular differentiation and the antileukemia activity of ATRA may be mediated by its ability to induce the differentiation of the promyelocytic leukemia cells through down-regulation of the hTERT gene.

X-chromosome inactivation and escape.

X-chromosome inactivation, which was discovered by Mary Lyon in 1961 results in random silencing of one X chromosome in female mammals. This review is dedicated to Mary Lyon, who passed away last year. She predicted many of the features of X inactivation, for e.g., the existence of an X inactivation center, the role of L1 elements in spreading of silencing and the existence of genes that escape X inactivation. Starting from her published work here we summarize advances in the field.

Identification of genes escaping X inactivation by allelic expression analysis in a novel hybrid mouse model.

X chromosome inactivation (XCI) is a female-specific mechanism that serves to balance gene dosage between the sexes whereby one X chromosome in females is inactivated during early development. Despite this silencing, a small portion of genes escape inactivation and remain expressed from the inactive X (Xi). Little is known about the distribution of escape from XCI in different tissues in vivo and about the mechanisms that control tissue-specific differences. Using a new binomial model in conjunction with a mouse model with identifiable alleles and skewed X inactivation we are able to survey genes that escape XCI in vivo. We show that escape from X inactivation can be a common feature of some genes, whereas others escape in a tissue specific manner. Furthermore, we characterize the chromatin environment of escape genes and show that expression from the Xi correlates with factors associated with open chromatin and that CTCF co-localizes with escape genes. Here, we provide a detailed description of the experimental design and data analysis pipeline we used to assay allele-specific expression and epigenetic characteristics of genes escaping X inactivation. The data is publicly available through the GEO database under ascension numbers GSM1014171, GSE44255, and GSE59779. Interpretation and discussion of these data are included in a previously published study (Berletch et al., 2015) [1].