Fast and precise single-cell data analysis using a hierarchical autoencoder.

A primary challenge in single-cell RNA sequencing (scRNA-seq) studies comes from the massive amount of data and the excess noise level. To address this challenge, we introduce an analysis framework, named single-cell Decomposition using Hierarchical Autoencoder (scDHA), that reliably extracts representative information of each cell. The scDHA pipeline consists of two core modules. The first module is a non-negative kernel autoencoder able to remove genes or components that have insignificant contributions to the part-based representation of the data. The second module is a stacked Bayesian autoencoder that projects the data onto a low-dimensional space (compressed). To diminish the tendency to overfit of neural networks, we repeatedly perturb the compressed space to learn a more generalized representation of the data. In an extensive analysis, we demonstrate that scDHA outperforms state-of-the-art techniques in many research sub-fields of scRNA-seq analysis, including cell segregation through unsupervised learning, visualization of transcriptome landscape, cell classification, and pseudo-time inference.

Interhomolog Homologous Recombination in Mouse Embryonic Stem Cells.

Homologous recombination is a critical mechanism for the repair of DNA double-strand breaks (DSBs). It occurs predominantly between identical sister chromatids and at lower frequency can also occur between homologs. Interhomolog homologous recombination (IH-HR) has the potential lead to substantial loss of genetic information, i.e., loss of heterozygosity (LOH), when it is accompanied by crossing over. In this chapter, we describe a system to study IH-HR induced by a defined DSB in mouse embryonic stem cells derived from F1 hybrid mice. This system is based on the placement of mutant selectable marker genes, one of which contains an I-SceI endonuclease cleavage site, on the two homologs such that repair of the I-SceI-generated DSB from the homolog leads to drug resistance. Loss of heterozygosity arising during IH-HR is analyzed using a PCR-based approach. Finally, we present a strategy to analyze the role of BLM helicase in this system.

Measurement of Homologous Recombination at Stalled Mammalian Replication Forks.

Site-specific replication fork barriers (RFBs) have proven valuable tools for studying mechanisms of repair at sites of replication fork stalling in prokaryotes and yeasts. We adapted the Escherichia coli Tus-Ter RFB for use in mammalian cells and used it to trigger site-specific replication fork stalling and homologous recombination (HR) at a defined chromosomal locus in mammalian cells. By comparing HR responses induced at the Tus-Ter RFB with those induced by a site-specific double-strand break (DSB), we have begun to uncover how the mechanisms of mammalian stalled fork repair differ from those underlying the repair of a replication-independent DSB. Here, we outline how to transiently express the Tus protein in mES cells, how to use flow cytometry to score conservative and aberrant repair outcomes, and how to quantify distinct repair outcomes in response to replication fork stalling at the inducible Tus-Ter chromosomal RFB.

Pluripotency of embryonic stem cells lacking clathrin-mediated endocytosis cannot be rescued by restoring cellular stiffness.

Mouse embryonic stem cells (mESCs) display unique mechanical properties, including low cellular stiffness in contrast to differentiated cells, which are stiffer. We have previously shown that mESCs lacking the clathrin heavy chain (Cltc), an essential component for clathrin-mediated endocytosis (CME), display a loss of pluripotency and an enhanced expression of differentiation markers. However, it is not known whether physical properties such as cellular stiffness also change upon loss of Cltc, similar to what is seen in differentiated cells, and if so, how these altered properties specifically impact pluripotency. Using atomic force microscopy (AFM), we demonstrate that mESCs lacking Cltc display higher Young's modulus, indicative of greater cellular stiffness, compared with WT mESCs. The increase in stiffness was accompanied by the presence of actin stress fibers and accumulation of the inactive, phosphorylated, actin-binding protein cofilin. Treatment of Cltc knockdown mESCs with actin polymerization inhibitors resulted in a decrease in the Young's modulus to values similar to those obtained with WT mESCs. However, a rescue in the expression profile of pluripotency factors was not obtained. Additionally, whereas WT mouse embryonic fibroblasts could be reprogrammed to a state of pluripotency, this was inhibited in the absence of Cltc. This indicates that the presence of active CME is essential for the pluripotency of embryonic stem cells. Additionally, whereas physical properties may serve as a simple readout of the cellular state, they may not always faithfully recapitulate the underlying molecular fate.

Millefy: visualizing cell-to-cell heterogeneity in read coverage of single-cell RNA sequencing datasets.

BACKGROUND: Read coverage of RNA sequencing data reflects gene expression and RNA processing events. Single-cell RNA sequencing (scRNA-seq) methods, particularly "full-length" ones, provide read coverage of many individual cells and have the potential to reveal cellular heterogeneity in RNA transcription and processing. However, visualization tools suited to highlighting cell-to-cell heterogeneity in read coverage are still lacking. RESULTS: Here, we have developed Millefy, a tool for visualizing read coverage of scRNA-seq data in genomic contexts. Millefy is designed to show read coverage of all individual cells at once in genomic contexts and to highlight cell-to-cell heterogeneity in read coverage. By visualizing read coverage of all cells as a heat map and dynamically reordering cells based on diffusion maps, Millefy facilitates discovery of "local" region-specific, cell-to-cell heterogeneity in read coverage. We applied Millefy to scRNA-seq data sets of mouse embryonic stem cells and triple-negative breast cancers and showed variability of transcribed regions including antisense RNAs, 3 ' UTR lengths, and enhancer RNA transcription. CONCLUSIONS: Millefy simplifies the examination of cellular heterogeneity in RNA transcription and processing events using scRNA-seq data. Millefy is available as an R package (https://github.com/yuifu/millefy) and as a Docker image for use with Jupyter Notebook (https://hub.docker.com/r/yuifu/datascience-notebook-millefy).

Evaluation of 3D Chromatin Interactions Using Hi-C.

The invention of Hi-C has greatly facilitated 3D genome research through an unbiased probing of 3D chromatin interactions. It produces enormous amount of sequencing data that capture multiscale chromatin conformation structures. In the last decade, numerous computational methods have been developed to analyze Hi-C data and predict A/B compartments, topologically associating domains (TADs), and significant chromatin contacts. This chapter introduced the iHiC package that provides several utilities to facilitate Hi-C data analysis with public software and demonstrated its application to a Hi-C dataset generated for mouse embryonic stem (ES) cells.

A gene regulatory network controls the balance between mesendoderm and ectoderm at pluripotency exit.

During embryogenesis, differentiation of pluripotent cells into somatic cell types depends both on signaling cues and intrinsic gene expression programs. While the molecular underpinnings of pluripotency are well mapped, much less is known on how mouse embryonic stem cells (mESCs) differentiate. Using RNA-Seq profiling during specification to the three germ layers, we showed that mESCs switched on condition-specific gene expression programs from the onset of the differentiation procedure and that primed pluripotency did not constitute an obligatory intermediate state. After inferring the gene network controlling mESC differentiation, we tested the role of the highly connected nodes by deleting them in a triple knock-in Sox1-Brachyury-Eomes mESC line reporting on ectoderm, mesoderm, and endoderm fates. This led to the identification of regulators of mESC differentiation that acted at several levels: Sp1 as a global break on differentiation, Nr5a2 controlling ectoderm specification, and notably Fos:Jun and Zfp354c as opposite switches between ectoderm and mesendoderm fate.

Aptamer-Based Fluorescent Biosensing of Adenosine Triphosphate and Cytochrome c via Aggregation-Induced Emission Enhancement on Novel Label-Free DNA-Capped Silver Nanoclusters/Graphene Oxide Nanohybrids.

Four fluorescent DNA-stabilized fluorescent silver nanoclusters (DNA-AgNCs) were designed and synthesized with differences in lengths of cytosine-rich DNA strand (as the stabilizing agent) and target-specific strand DNA aptamers for adenosine triphosphate (ATP) and cytochrome c (Cyt c). After their nanohybrid formation with graphene oxide (GO), it was unexpectedly found that, depending on the composition of the base and length of the strand DNA aptamer, the fluorescence intensity of three of the nanohybrids significantly enhanced. Our experimental observations and quantum mechanical calculations provided an insight into the mechanisms underlying the behavior of DNA-AgNCs/GO nanohybrids. The enhanced fluorescence was found to be attributed to the aggregation-induced emission enhancement (AIE) characteristic of the DNA-AgNCs adsorbed on the GO surface, as confirmed evidently by both fluorescence and transmission electron microscopies. The AIE is a result of hardness and oxidation properties of GO, which lead to enhanced argenophilic interaction and thus to increased Ag(I)-DNA complex shell aggregation. Consequently, two of the DNA-AgNCs/GO nanohybrids were successfully extended to construct highly selective, sensitive, label-free, and simple aptasensors for biosensing of ATP (LOD = 0.42 nM) and Cyt c (LOD = 2.3 nM) in lysed Escherichia coli DH5 α cells and mouse embryonic stem cells, respectively. These fundamental findings are expected to significantly influence the designing and engineering of new AgNCs/GO-based AIE biosensors.

Glycine cleavage system determines the fate of pluripotent stem cells via the regulation of senescence and epigenetic modifications.

Metabolic remodelling has emerged as critical for stem cell pluripotency; however, the underlying mechanisms have yet to be fully elucidated. Here, we found that the glycine cleavage system (GCS) is highly activated to promote stem cell pluripotency and during somatic cell reprogramming. Mechanistically, we revealed that the expression of Gldc, a rate-limiting GCS enzyme regulated by Sox2 and Lin28A, facilitates this activation. We further found that the activated GCS catabolizes glycine to fuel H3K4me3 modification, thus promoting the expression of pluripotency genes. Moreover, the activated GCS helps to cleave excess glycine and prevents methylglyoxal accumulation, which stimulates senescence in stem cells and during reprogramming. Collectively, our results demonstrate a novel mechanism whereby GCS activation controls stem cell pluripotency by promoting H3K4me3 modification and preventing cellular senescence.

Transcription Restart Establishes Chromatin Accessibility after DNA Replication.

DNA replication is highly disruptive to chromatin, leading to eviction of nucleosomes, RNA polymerase, and regulatory factors. When and how transcription resumes on DNA following DNA replication is unknown. Here we develop a replication-coupled assay for transposase-accessible chromatin (repli-ATAC-seq) to investigate active chromatin restoration post-replication in mouse embryonic stem cells. We find that nascent chromatin is inaccessible and transcriptionally silenced, with accessibility and RNA polymerase occupancy re-appearing within 30 minutes. Chromatin accessibility restores differentially genome wide, with super enhancers regaining transcription factor occupancy faster than other genomic features. We also identify opportunistic and transiently accessible chromatin within gene bodies after replication. Systematic inhibition of transcription shows that transcription restart is required to re-establish active chromatin states genome wide and resolve opportunistic binding events resulting from DNA replication. Collectively, this establishes a central role for transcription in overcoming the genome-wide chromatin inaccessibility imposed by DNA replication every cell division.

DamC reveals principles of chromatin folding in vivo without crosslinking and ligation.

Current understanding of chromosome folding is largely reliant on chromosome conformation capture (3C)-based experiments, where chromosomal interactions are detected as ligation products after chromatin crosslinking. To measure chromosome structure in vivo, quantitatively and without crosslinking and ligation, we implemented a modified version of DNA adenine methyltransferase identification (DamID) named DamC, which combines DNA methylation-based detection of chromosomal interactions with next-generation sequencing and biophysical modeling of methylation kinetics. DamC performed in mouse embryonic stem cells provides the first in vivo validation of the existence of topologically associating domains (TADs), CTCF loops and confirms 3C-based measurements of the scaling of contact probabilities. Combining DamC with transposon-mediated genomic engineering shows that new loops can be formed between ectopic and endogenous CTCF sites, which redistributes physical interactions within TADs. DamC provides the first crosslinking- and ligation-free demonstration of the existence of key structural features of chromosomes and provides novel insights into how chromosome structure within TADs can be manipulated.

5-Formylcytosine organizes nucleosomes and forms Schiff base interactions with histones in mouse embryonic stem cells.

Nucleosomes are the basic unit of chromatin that help the packaging of genetic material while controlling access to the genetic information. The underlying DNA sequence, together with transcription-associated proteins and chromatin remodelling complexes, are important factors that influence the organization of nucleosomes. Here, we show that the naturally occurring DNA modification, 5-formylcytosine (5fC) is linked to tissue-specific nucleosome organization. Our study reveals that 5fC is associated with increased nucleosome occupancy in vitro and in vivo. We demonstrate that 5fC-associated nucleosomes at enhancers in the mammalian hindbrain and heart are linked to elevated gene expression. Our study also reveals the formation of a reversible-covalent Schiff base linkage between lysines of histone proteins and 5fC within nucleosomes in a cellular environment. We define their specific genomic loci in mouse embryonic stem cells and look into the biological consequences of these DNA-histone Schiff base sites. Collectively, our findings show that 5fC is a determinant of nucleosome organization and plays a role in establishing distinct regulatory regions that control transcription.

Characterization of a core region in the A2UCOE that confers effective anti-silencing activity.

We have previously shown that reliability of the A2UCOE in driving transgene expression can be attributed to its resistance to DNA methylation, and its ability to confer this property to linked regulatory sequences. In order to gain a better understanding of how resistance to DNA methylation from the A2UCOE is conferred, and whether the anti-silencing effect from the A2UCOE is confined within a core region, we evaluated the anti-silencing effect of different sub-domains. We found that maximal epigenetic regulatory activity was contained within a 455 bp element derived from the CBX3 region when tested in the context of a lentiviral vector in murine embryonic stem (ES) cells and human inducible pluripotent stem (iPS) cells. This region possessed an active chromatin signature, and operated effectively in cis to protect linked heterologous regulatory elements from methylation, thereby conferring stable transgene expression. Defined UCOE elements may be particularly useful for use in vectors where gene expression is desired in methylation-prone chromatin environments such as those encountered in pluripotent stem cells.

miR-290 contributes to the low abundance of cyclin D1 protein in mouse embryonic stem cells.

Mouse miR-290 cluster miRNAs are expressed specifically in early embryos and embryonic germ cells. These miRNAs play critical roles in the maintenance of pluripotency and self-renewal. Here, we showed that Cyclin D1 is a direct target gene of miR-290 cluster miRNAs. Negative relationships between the expression of Cyclin D1 protein and miR-290 cluster miRNAs in pluripotent and non-pluripotent cells, as well as in differentiating CGR8 cells were observed. Inhibition of miR-290 cluster miRNAs could arrest cells at the G1 phase and slow down the cell proliferation in CGR8 mouse stem cells. Since miR-290 cluster miRNAs are the most dominant stem-cell-specific miRNAs, our results revealed an important cause for the absence of Cyclin D1 in mouse embryonic stem cells.

Mouse intestinal niche cells express a distinct α1,2-fucosylated glycan recognized by a lectin from Burkholderia cenocepacia.

In this study, we examined the distribution of fucosylated glycans in mouse intestines using a lectin, BC2LCN (N-terminal domain of the lectin BC2L-C from Burkholderia cenocepacia), as a probe. BC2LCN is specific for glycans with a terminal Fucα1,2Galß1,3-motif and it is a useful marker for discriminating the undifferentiated status of human induced/embryonic stem cells. Apparent BC2LCN reactivity was detected in the secretory granules of goblet cells in the ileum but not those in the colon. We also found distinctive reactivity in the crypt bottom, which is known as the stem cell zone, of the colon and the ileum. Other lectins for fucosylated glycans, including Ulex europaeus agglutinin-I, Pholiota squarrosa lectin and Aleuria aurantia lectin, did not exhibit similar reactivity in the crypt bottom. Remarkably, BC2LCN-positive epithelial cells could be labeled with a niche cell marker, c-Kit/CD117. Overall, our results indicate that intestinal niche cells express distinct fucosylated glycans recognized by BC2LCN. Increasing evidence suggests that the self-renewal and proliferation of stem cells depend on specific signals derived from niche cells. Our results highlight novel molecular properties of intestinal niche cells in terms of their glycosylation, which may help to understand the regulation of intestinal stem cells. The distinct expression of glycans may reflect the functional roles of niche cells. BC2LCN is a valuable tool for investigating the functional significance of protein glycosylation in stem cell regulation.

Thermoresponsive hydrogel maintains the mouse embryonic stem cell "naïve" pluripotency phenotype.

A chemically defined thermoresponsive hydrogel, poly(AEtMA-Cl-co-DEAEA) cross-linked with N,N'-methylenebisacrylamide, which allows enzyme-free passaging, was used as a substrate to culture murine embryonic stem cells (mESCs) under defined and undefined conditions. Analysis of 14 stem cell markers showed that the mESCs remained in a "naïve" state of pluripotency with differentiation potential to form endoderm, mesoderm, and ectoderm derived lineages. These results validate the use of a chemically defined hydrogel for standardised and inexpensive mESC culture.