Replication and transcription: shaping the landscape of the genome

Globin vector regulatory elements are active in early hematopoietic progenitor cells

The globin genes are archetypal tissue-specific genes that are silent in most tissues but for late-stage erythroblasts upon terminal erythroid differentiation. The transcriptional activation of the β-globin gene is under the control of proximal and distal regulatory elements located on chromosome 11p15.4, including the β-globin locus control region (LCR). The incorporation of selected LCR elements in lentiviral vectors encoding β and β-like globin genes has enabled successful genetic treatment of the β-thalassemias and sickle cell disease. However, recent occurrences of benign clonal expansions in thalassemic patients and myelodysplastic syndrome in patients with sickle cell disease call attention to the non-erythroid functions of these powerful vectors. Here we demonstrate that lentivirally encoded LCR elements, in particular HS1 and HS2, can be activated in early hematopoietic cells including hematopoietic stem cells and myeloid progenitors. This activity is position-dependent and results in the transcriptional activation of a nearby reporter gene in these progenitor cell populations. We further show that flanking a globin vector with an insulator can effectively restrain this non-erythroid activity without impairing therapeutic globin expression. Globin lentiviral vectors harboring powerful LCR HS elements may thus expose to the risk of trans-activating cancer-related genes, which can be mitigated by a suitable insulator.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Liquid-Liquid Phase Separation in Chromatin

In eukaryotic cells, protein and RNA factors involved in genome activities like transcription, RNA processing, DNA replication, and repair accumulate in self-organizing membraneless chromatin subcompartments. These structures contribute to efficiently conduct chromatin-mediated reactions and to establish specific cellular programs. However, the underlying mechanisms for their formation are only partly understood. Recent studies invoke liquid-liquid phase separation (LLPS) of proteins and RNAs in the establishment of chromatin activity patterns. At the same time, the folding of chromatin in the nucleus can drive genome partitioning into spatially distinct domains. Here, the interplay between chromatin organization, chromatin binding, and LLPS is discussed by comparing and contrasting three prototypical chromatin subcompartments: the nucleolus, clusters of active RNA polymerase II, and pericentric heterochromatin domains. It is discussed how the different ways of chromatin compartmentalization are linked to transcription regulation, the targeting of soluble factors to certain parts of the genome, and to disease-causing genetic aberrations.

Computational methods for predicting 3D genomic organization from high-resolution chromosome conformation capture data

The advent of high-resolution chromosome conformation capture assays (such as 5C, Hi-C and Pore-C) has allowed for unprecedented sequence-level investigations into the structure-function relationship of the genome. In order to comprehensively understand this relationship, computational tools are required that utilize data generated from these assays to predict 3D genome organization (the 3D genome reconstruction problem). Many computational tools have been developed that answer this need, but a comprehensive comparison of their underlying algorithmic approaches has not been conducted. This manuscript provides a comprehensive review of the existing computational tools (from November 2006 to September 2019, inclusive) that can be used to predict 3D genome organizations from high-resolution chromosome conformation capture data. Overall, existing tools were found to use a relatively small set of algorithms from one or more of the following categories: dimensionality reduction, graph/network theory, maximum likelihood estimation (MLE) and statistical modeling. Solutions in each category are far from maturity, and the breadth and depth of various algorithmic categories have not been fully explored. While the tools for predicting 3D structure for a genomic region or single chromosome are diverse, there is a general lack of algorithmic diversity among computational tools for predicting the complete 3D genome organization from high-resolution chromosome conformation capture data.

Nanoscopic subcellular imaging enabled by ion beam tomography

Multiplexed ion beam imaging (MIBI) has been previously used to profile multiple parameters in two dimensions in single cells within tissue slices. Here, a mathematical and technical framework for three-dimensional (3D) subcellular MIBI is presented. Ion-beam tomography (IBT) compiles ion beam images that are acquired iteratively across successive, multiple scans, and later assembled into a 3D format without loss of depth resolution. Algorithmic deconvolution, tailored for ion beams, is then applied to the transformed ion image series, yielding 4-fold enhanced ion beam data cubes. To further generate 3D sub-ion-beam-width precision visuals, isolated ion molecules are localized in the raw ion beam images, creating an approach coined as SILM, secondary ion beam localization microscopy, providing sub-25 nm accuracy in original ion images. Using deep learning, a parameter-free reconstruction method for ion beam tomograms with high accuracy is developed for low-density targets. In cultured cancer cells and tissues, IBT enables accessible visualization of 3D volumetric distributions of genomic regions, RNA transcripts, and protein factors with 5 nm axial resolution using isotope-enrichments and label-free elemental analyses. Multiparameter imaging of subcellular features at near macromolecular resolution is implemented by the IBT tools as a general biocomputation pipeline for imaging mass spectrometry.

Superresolution imaging reveals spatiotemporal propagation of human replication foci mediated by CTCF-organized chromatin structures

Mammalian DNA replication is initiated at numerous replication origins, which are clustered into thousands of replication domains (RDs) across the genome. However, it remains unclear whether the replication origins within each RD are activated stochastically or preferentially near certain chromatin features. To understand how DNA replication in single human cells is regulated at the sub-RD level, we directly visualized and quantitatively characterized the spatiotemporal organization, morphology, and in situ epigenetic signatures of individual replication foci (RFi) across S-phase at superresolution using stochastic optical reconstruction microscopy. Importantly, we revealed a hierarchical radial pattern of RFi propagation dynamics that reverses directionality from early to late S-phase and is diminished upon caffeine treatment or CTCF knockdown. Together with simulation and bioinformatic analyses, our findings point to a "CTCF-organized REplication Propagation" (CoREP) model, which suggests a nonrandom selection mechanism for replication activation at the sub-RD level during early S-phase, mediated by CTCF-organized chromatin structures. Collectively, these findings offer critical insights into the key involvement of local epigenetic environment in coordinating DNA replication across the genome and have broad implications for our conceptualization of the role of multiscale chromatin architecture in regulating diverse cell nuclear dynamics in space and time.

An intergenic non-coding RNA promoter required for histone modifications in the human β-globin chromatin domain

Transcriptome analyses show a surprisingly large proportion of the mammalian genome is transcribed; much more than can be accounted for by genes and introns alone. Most of this transcription is non-coding in nature and arises from intergenic regions, often overlapping known protein-coding genes in sense or antisense orientation. The functional relevance of this widespread transcription is unknown. Here we characterize a promoter responsible for initiation of an intergenic transcript located approximately 3.3 kb and 10.7 kb upstream of the adult-specific human β-globin genes. Mutational analyses in β-YAC transgenic mice show that alteration of intergenic promoter activity results in ablation of H3K4 di- and tri-methylation and H3 hyperacetylation extending over a 30 kb region immediately downstream of the initiation site, containing the adult δ- and β-globin genes. This results in dramatically decreased expression of the adult genes through position effect variegation in which the vast majority of definitive erythroid cells harbor inactive adult globin genes. In contrast, expression of the neighboring ε- and γ-globin genes is completely normal in embryonic erythroid cells, indicating a developmentally specific variegation of the adult domain. Our results demonstrate a role for intergenic non-coding RNA transcription in the propagation of histone modifications over chromatin domains and epigenetic control of β-like globin gene transcription during development.

ATR-Chk1 activation mitigates replication stress caused by mismatch repair-dependent processing of DNA damage

The mismatch repair pathway (MMR) is essential for removing DNA polymerase errors, thereby maintaining genomic stability. Loss of MMR function increases mutation frequency and is associated with tumorigenesis. However, how MMR is executed at active DNA replication forks is unclear. This has important implications for understanding how MMR repairs O6-methylguanine/thymidine (MeG/T) mismatches created upon exposure to DNA alkylating agents. If MeG/T lesion recognition by MMR initiates mismatch excision, the reinsertion of a mismatched thymidine during resynthesis could initiate futile repair cycles. One consequence of futile repair cycles might be a disruption of overall DNA replication in the affected cell. Herein, we show that in MMR-proficient HeLa cancer cells, treatment with a DNA alkylating agent slows S phase progression, yet cells still progress into the next cell cycle. In the first S phase following treatment, they activate ataxia telangiectasia and Rad3-related (ATR)-Checkpoint Kinase 1 (Chk1) signaling, which limits DNA damage, while inhibition of ATR kinase activity accelerates DNA damage accumulation and sensitivity to the DNA alkylating agent. We also observed that exposure of human embryonic stem cells to alkylation damage severely compromised DNA replication in a MMR-dependent manner. These cells fail to activate the ATR-Chk1 signaling axis, which may limit their ability to handle replication stress. Accordingly, they accumulate double-strand breaks and undergo immediate apoptosis. Our findings implicate the MMR-directed response to alkylation damage as a replication stress inducer, suggesting that repeated MMR processing of mismatches may occur that can disrupt S phase progression.

Laser Microirradiation to Study In Vivo Cellular Responses to Simple and Complex DNA Damage

DNA damage induces specific signaling and repair responses in the cell, which is critical for protection of genome integrity. Laser microirradiation became a valuable experimental tool to investigate the DNA damage response (DDR) in vivo. It allows real-time high-resolution single-cell analysis of macromolecular dynamics in response to laser-induced damage confined to a submicrometer region in the cell nucleus. However, various laser conditions have been used without appreciation of differences in the types of damage induced. As a result, the nature of the damage is often not well characterized or controlled, causing apparent inconsistencies in the recruitment or modification profiles. We demonstrated that different irradiation conditions (i.e., different wavelengths as well as different input powers (irradiances) of a femtosecond (fs) near-infrared (NIR) laser) induced distinct DDR and repair protein assemblies. This reflects the type of DNA damage produced. This protocol describes how titration of laser input power allows induction of different amounts and complexities of DNA damage, which can easily be monitored by detection of base and crosslinking damages, differential poly (ADP-ribose) (PAR) signaling, and pathway-specific repair factor assemblies at damage sites. Once the damage conditions are determined, it is possible to investigate the effects of different damage complexity and differential damage signaling as well as depletion of upstream factor(s) on any factor of interest.

Identification and characterization of long non-coding RNA in prenatal and postnatal skeletal muscle of sheep

lncRNAs are a class of transcriptional RNA molecules of >200 nucleotides in length. However, the overall expression pattern and function of lncRNAs in sheep muscle is not clear. Here, we identified 1566 lncRNAs and 404 differentially expressed lncRNAs in sheep muscle from prenatal (110 days of fetus) and postnatal (2 to 3 years old of adult sheep) developmental stages by using RNA-seq technology. Several lncRNAs were identified by using RT-PCR and DNA sequencing. The expression levels of several lncRNAs were confirmed by qRT-PCR. We analyzed the effect of lncRNAs that act cis to the target genes. lncRNA targeting genes were involved in signaling pathways associated with growth and development of muscle by GO and KEGG enrichment analysis. Through our study, we provide a comprehensive expression profile of muscle lncRNAs in sheep, which provides valuable resources for further understanding genetic regulation of muscle growth and development from the perspective of lncRNA.

Interaction of the epigenetic integrator UHRF1 with the MYST domain of TIP60 inside the cell

Background

The nuclear epigenetic integrator UHRF1 is known to play a key role with DNMT1 in maintaining the DNA methylation patterns during cell division. Among UHRF1 partners, TIP60 takes part in epigenetic regulations through its acetyltransferase activity. Both proteins are involved in multiple cellular functions such as chromatin remodeling, DNA damage repair and regulation of stability and activity of other proteins. The aim of this work was to investigate the interaction between UHRF1 and TIP60 in order to elucidate the dialogue between these two proteins.

Methods

Biochemical (immunoprecipitation and pull-down assays) and microscopic (confocal and fluorescence lifetime imaging microscopy; FLIM) techniques were used to analyze the interaction between TIP60 and UHRF1 in vitro and in vivo. Global methylation levels were assessed by using a specific kit. The results were statistically analyzed using Graphpad prism and Origin.

Results

Our study shows that UHRF1, TIP60 and DNMT1 were found in the same epigenetic macro-molecular complex. In vitro pull-down assay showed that deletion of either the zinc finger in MYST domain or deletion of whole MYST domain from TIP60 significantly reduced its interaction with UHRF1. Confocal and FLIM microscopy showed that UHRF1 co-localized with TIP60 in the nucleus and confirmed that both proteins interacted together through the MYST domain of TIP60. Moreover, overexpression of TIP60 reduced the DNA methylation levels in HeLa cells along with downregulation of UHRF1 and DNMT1.

Conclusion

Our data demonstrate for the first time that TIP60 through its MYST domain directly interacts with UHRF1 which might be of high interest for the development of novel oncogenic inhibitors targeting this interaction.

Orientation and repositioning of chromosomes correlate with cell geometry-dependent gene expression

Extracellular matrix signals from the microenvironment regulate gene expression patterns and cell behavior. Using a combination of experiments and geometric models, we demonstrate correlations between cell geometry, three-dimensional (3D) organization of chromosome territories, and gene expression. Fluorescence in situ hybridization experiments showed that micropatterned fibroblasts cultured on anisotropic versus isotropic substrates resulted in repositioning of specific chromosomes, which contained genes that were differentially regulated by cell geometries. Experiments combined with ellipsoid packing models revealed that the mechanosensitivity of chromosomes was correlated with their orientation in the nucleus. Transcription inhibition experiments suggested that the intermingling degree was more sensitive to global changes in transcription than to chromosome radial positioning and its orientations. These results suggested that cell geometry modulated 3D chromosome arrangement, and their neighborhoods correlated with gene expression patterns in a predictable manner. This is central to understanding geometric control of genetic programs involved in cellular homeostasis and the associated diseases.

Long noncoding RNAs in psychiatric disorders

Long noncoding RNAs (lncRNAs) are nonprotein coding transcripts longer than 200 nucleotides. Many of these lncRNAs have regulatory functions and have recently emerged as major players in governing fundamental biological processes. Here, we review the definition, distribution, identification, databases, analysis, classification, and functions of lncRNAs. We also discuss the potential roles of lncRNAs in the etiological processes of psychiatric disorders and the implications for clinical diagnosis and treatment.

Crystal structure of the DNA binding domain of the transcription factor T-bet suggests simultaneous recognition of distant genome sites

The transcription factor T-bet (Tbox protein expressed in T cells) is one of the master regulators of both the innate and adaptive immune responses. It plays a central role in T-cell lineage commitment, where it controls the T_H1 response, and in gene regulation in plasma B-cells and dendritic cells. T-bet is a member of the Tbox family of transcription factors; however, T-bet coordinately regulates the expression of many more genes than other Tbox proteins. A central unresolved question is how T-bet is able to simultaneously recognize distant Tbox binding sites, which may be located thousands of base pairs away. We have determined the crystal structure of the Tbox DNA binding domain (DBD) of T-bet in complex with a palindromic DNA. The structure shows a quaternary structure in which the T-bet dimer has its DNA binding regions splayed far apart, making it impossible for a single dimer to bind both sites of the DNA palindrome. In contrast to most other Tbox proteins, a single T-bet DBD dimer binds simultaneously to identical half-sites on two independent DNA. A fluorescence-based assay confirms that T-bet dimers are able to bring two independent DNA molecules into close juxtaposition. Furthermore, chromosome conformation capture assays confirm that T-bet functions in the direct formation of chromatin loops in vitro and in vivo. The data are consistent with a looping/synapsing model for transcriptional regulation by T-bet in which a single dimer of the transcription factor can recognize and coalesce distinct genetic elements, either a promoter plus a distant regulatory element, or promoters on two different genes.

Novel genomic rearrangements mediated by multiple genetic elements in Streptococcus pyogenes M23ND confer potential for evolutionary persistence

Symmetric genomic rearrangements around replication axes in genomes are commonly observed in prokaryotic genomes, including Group A Streptococcus (GAS). However, asymmetric rearrangements are rare. Our previous studies showed that the hypervirulent invasive GAS strain, M23ND, containing an inactivated transcriptional regulator system, covRS, exhibits unique extensive asymmetric rearrangements, which reconstructed a genomic structure distinct from other GAS genomes. In the current investigation, we identified the rearrangement events and examined the genetic consequences and evolutionary implications underlying the rearrangements. By comparison with a close phylogenetic relative, M18-MGAS8232, we propose a molecular model wherein a series of asymmetric rearrangements have occurred in M23ND, involving translocations, inversions and integrations mediated by multiple factors, viz., rRNA-comX (factor for late competence), transposons and phage-encoded gene segments. Assessments of the cumulative gene orientations and GC skews reveal that the asymmetric genomic rearrangements did not affect the general genomic integrity of the organism. However, functional distributions reveal re-clustering of a broad set of CovRS-regulated actively transcribed genes, including virulence factors and metabolic genes, to the same leading strand, with high confidence (p-value ~10-10). The re-clustering of the genes suggests a potential selection advantage for the spatial proximity to the transcription complexes, which may contain the global transcriptional regulator, CovRS, and other RNA polymerases. Their proximities allow for efficient transcription of the genes required for growth, virulence and persistence. A new paradigm of survival strategies of GAS strains is provided through multiple genomic rearrangements, while, at the same time, maintaining genomic integrity.

Integrating epigenomic data and 3D genomic structure with a new measure of chromatin assortativity

BACKGROUND

Network analysis is a powerful way of modeling chromatin interactions. Assortativity is a network property used in social sciences to identify factors affecting how people establish social ties. We propose a new approach, using chromatin assortativity, to integrate the epigenomic landscape of a specific cell type with its chromatin interaction network and thus investigate which proteins or chromatin marks mediate genomic contacts.

RESULTS

We use high-resolution promoter capture Hi-C and Hi-Cap data as well as ChIA-PET data from mouse embryonic stem cells to investigate promoter-centered chromatin interaction networks and calculate the presence of specific epigenomic features in the chromatin fragments constituting the nodes of the network. We estimate the association of these features with the topology of four chromatin interaction networks and identify features localized in connected areas of the network. Polycomb group proteins and associated histone marks are the features with the highest chromatin assortativity in promoter-centered networks. We then ask which features distinguish contacts amongst promoters from contacts between promoters and other genomic elements. We observe higher chromatin assortativity of the actively elongating form of RNA polymerase 2 (RNAPII) compared with inactive forms only in interactions between promoters and other elements.

CONCLUSIONS

Contacts among promoters and between promoters and other elements have different characteristic epigenomic features. We identify a possible role for the elongating form of RNAPII in mediating interactions among promoters, enhancers, and transcribed gene bodies. Our approach facilitates the study of multiple genome-wide epigenomic profiles, considering network topology and allowing the comparison of chromatin interaction networks.

Detection of RNA-DNA association by a proximity ligation-based method

We describe a proximity ligation assay (PLA)-based method of assessing association of DNA and RNA in single cells during the cell cycle. Pulse-labeling of DNA with EdU and RNA with BrU and testing their close proximity by PLA demonstrates that RNA synthesis in individual cells resumes about 30–45 min after DNA replication. Consistent with this conclusion, RNA Pol II phosphorylated at Ser2 of its CTD is detected at the same time as RNA transcripts on nascent DNA. Our results also show that RNA is associated with DNA foci during all stages of mitosis.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Joint annotation of chromatin state and chromatin conformation reveals relationships among domain types and identifies domains of cell-type-specific expression

The genomic neighborhood of a gene influences its activity, a behavior that is attributable in part to domain-scale regulation. Previous genomic studies have identified many types of regulatory domains. However, due to the difficulty of integrating genomics data sets, the relationships among these domain types are poorly understood. Semi-automated genome annotation (SAGA) algorithms facilitate human interpretation of heterogeneous collections of genomics data by simultaneously partitioning the human genome and assigning labels to the resulting genomic segments. However, existing SAGA methods cannot integrate inherently pairwise chromatin conformation data. We developed a new computational method, called graph-based regularization (GBR), for expressing a pairwise prior that encourages certain pairs of genomic loci to receive the same label in a genome annotation. We used GBR to exploit chromatin conformation information during genome annotation by encouraging positions that are close in 3D to occupy the same type of domain. Using this approach, we produced a model of chromatin domains in eight human cell types, thereby revealing the relationships among known domain types. Through this model, we identified clusters of tightly regulated genes expressed in only a small number of cell types, which we term “specific expression domains.” We found that domain boundaries marked by promoters and CTCF motifs are consistent between cell types even when domain activity changes. Finally, we showed that GBR can be used to transfer information from well-studied cell types to less well-characterized cell types during genome annotation, making it possible to produce high-quality annotations of the hundreds of cell types with limited available data.

Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Epigenetic regulation of the replication program during mammalian cell differentiation remains poorly understood. We performed an integrative analysis of eleven genome-wide epigenetic profiles at 100 kb resolution of Mean Replication Timing (MRT) data in six human cell lines. Compared to the organization in four chromatin states shared by the five somatic cell lines, embryonic stem cell (ESC) line H1 displays (i) a gene-poor but highly dynamic chromatin state (EC4) associated to histone variant H2AZ rather than a HP1-associated heterochromatin state (C4) and (ii) a mid-S accessible chromatin state with bivalent gene marks instead of a polycomb-repressed heterochromatin state. Plastic MRT regions (≲ 20% of the genome) are predominantly localized at the borders of U-shaped timing domains. Whereas somatic-specific U-domain borders are gene-dense GC-rich regions, 31.6% of H1-specific U-domain borders are early EC4 regions enriched in pluripotency transcription factors NANOG and OCT4 despite being GC poor and gene deserts. Silencing of these ESC-specific “master” replication initiation zones during differentiation corresponds to a loss of H2AZ and an enrichment in H3K9me3 mark characteristic of late replicating C4 heterochromatin. These results shed a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and lineage commitment.

During development, embryonic stem cell (ESC) enter a program of cell differentiation eventually leading to all the necessary differentiated cell types. Understanding the mechanisms responsible for the underlying modifications of the gene expression program is of fundamental importance, as it will likely have strong impact on the development of regenerative medicine. We show that besides some epigenetic regulation, ubiquitous master replication origins at replication timing U-domain borders shared by 6 human cell types are transcriptionally active open chromatin regions specified by a local enrichment in nucleosome free regions encoded in the DNA sequence suggesting that they have been selected during evolution. In contrast, ESC specific master replication origins bear a unique epigenetic signature (enrichment in CTCF, H2AZ, NANOG, OCT4, …) likely contributing to maintain ESC chromatin in a highly dynamic and accessible state that is refractory to polycomb and HP1 heterochromatin spreading. These ESC specific master origins thus appear as key genomic regions where epigenetic control of chromatin organization is at play to maintain pluripotency of stem cell lineages and to guide lineage commitment to somatic cell types.

Chromatin associations in Arabidopsis interphase nuclei

The arrangement of chromatin within interphase nuclei seems to be caused by topological constraints and related to gene expression depending on tissue and developmental stage. In yeast and animals it was found that homologous and heterologous chromatin association are required to realize faithful expression and DNA repair. To test whether such associations are present in plants we analyzed Arabidopsis thaliana interphase nuclei by FISH using probes from different chromosomes. We found that chromatin fiber movement and variable associations, although in general relatively seldom, may occur between euchromatin segments along chromosomes, sometimes even over large distances. The combination of euchromatin segments bearing high or low co-expressing genes did not reveal different association frequencies probably due to adjacent genes of deviating expression patterns. Based on previous data and on FISH analyses presented here, we conclude that the global interphase chromatin organization in A. thaliana is relatively stable, due to the location of its 10 centromeres at the nuclear periphery and of the telomeres mainly at the centrally localized nucleolus. Nevertheless, chromatin movement enables a flexible spatial genome arrangement in plant nuclei.

Large transcription units unify copy number variants and common fragile sites arising under replication stress

Copy number variants (CNVs) resulting from genomic deletions and duplications and common fragile sites (CFSs) seen as breaks on metaphase chromosomes are distinct forms of structural chromosome instability precipitated by replication inhibition. Although they share a common induction mechanism, it is not known how CNVs and CFSs are related or why some genomic loci are much more prone to their occurrence. Here we compare large sets of de novo CNVs and CFSs in several experimental cell systems to each other and to overlapping genomic features. We first show that CNV hotpots and CFSs occurred at the same human loci within a given cultured cell line. Bru-seq nascent RNA sequencing further demonstrated that although genomic regions with low CNV frequencies were enriched in transcribed genes, the CNV hotpots that matched CFSs specifically corresponded to the largest active transcription units in both human and mouse cells. Consistently, active transcription units >1 Mb were robust cell-type-specific predictors of induced CNV hotspots and CFS loci. Unlike most transcribed genes, these very large transcription units replicated late and organized deletion and duplication CNVs into their transcribed and flanking regions, respectively, supporting a role for transcription in replication-dependent lesion formation. These results indicate that active large transcription units drive extreme locus- and cell-type-specific genomic instability under replication stress, resulting in both CNVs and CFSs as different manifestations of perturbed replication dynamics.

Screening of lymph nodes metastasis associated lncRNAs in colorectal cancer patients

AIM: To screen lymph nodes metastasis associated long noncoding RNAs (lncRNAs) in colorectal cancer through microarray analysis.

METHODS: Metastatic lymph node (MLN), normal lymph node (NLN) and tumor tissues of 3 colorectal cancer (CRC) patients were collected during the operation and validated by pathological examinations. RNAs were extracted from MLN, NLN, and cancer tissues separately. RNA quantity and quality were measured with a NanoDrop ND-1000 spectrophotometer and RNA integrity was assessed by standard denaturing agarose electrophoresis. Agilent Feature Extraction Software (Version 11.0.1.1) was used to analyze acquired array images. Four differently expressed lncRNAs were confirmed by quantitative real-time polymerase chain reaction (qRT-PCR) in 26 subsets of MLN, NLN, and tumor tissues.

RESULTS: Of 33045 lncRNAs, 1133 were differentially expressed in MLN compared with NLN, of which 260 were up-regulated and 873 down-regulated (≥ 2 fold-change). Five hundred and forty-five lncRNAs were differentially expressed in MLN compared with tumor tissues, of which 460 were up-regulated and 85 down-regulated (≥ 2 fold-change). Compared with NLN and cancer tissues, 14 lncRNAs were specifically up-regulated and 5 specifically down-regulated in MLN. AK307796, ENST00000425785, and AK021444 were confirmed to be specifically up-regulated in MLN and ENST00000465846 specifically down-regulated in MLN by qRT-PCR in 26 CRC patients.

CONCLUSION: The specifically expressed lncRNAs in MLN may exert a partial or key role in the progress of lymph nodes metastasis of CRC.

Chromatin structure and replication origins: determinants of chromosome replication and nuclear organization

The DNA replication program is, in part, determined by the epigenetic landscape that governs local chromosome architecture and directs chromosome duplication. Replication must coordinate with other biochemical processes occurring concomitantly on chromatin, such as transcription and remodeling, to insure accurate duplication of both genetic and epigenetic features and to preserve genomic stability. The importance of genome architecture and chromatin looping in coordinating cellular processes on chromatin is illustrated by two recent sets of discoveries. First, chromatin-associated proteins that are not part of the core replication machinery were shown to affect the timing of DNA replication. These chromatin-associated proteins could be working in concert, or perhaps in competition, with the transcriptional machinery and with chromatin modifiers to determine the spatial and temporal organization of replication initiation events. Second, epigenetic interactions are mediated by DNA sequences that determine chromosomal replication. In this review, we summarize recent findings and current models linking spatial and temporal regulation of the replication program with epigenetic signaling. We discuss these issues in the context of the genome's three-dimensional structure with an emphasis on events occurring during the initiation of DNA replication.

Cell-cycle specific association of transcription factors and RNA polymerase ii with the human β-globin gene locus

The human β-globin genes are regulated by a locus control region (LCR) and are expressed at extremely high levels in erythroid cells. How transcriptional fidelity of highly expressed genes is regulated and maintained during the cell cycle is not completely understood. Here, we analyzed the association of transcription factor USF, the co-activator CBP, topoisomerase I (Topo I), basal transcription factor TFIIB, and RNA polymerase II (Pol II) with the β-globin gene locus at specific cell-cycle stages. The data demonstrate that while association of Pol II with globin locus associated chromatin decreased in mitotically arrested cells, it remained bound at lower levels at the γ-globin gene promoter. During early S-phase, association of CBP, USF, and Pol II with the globin gene locus decreased. The re-association of CBP and USF2 with the LCR preceded re-association of Pol II, suggesting that these proteins together mediate recruitment of Pol II to the β-globin gene locus during S-phase. Finally, we analyzed the association of Topo I with the globin gene locus during late S-phase. In general, Topo I association correlated with the binding of Pol II. Inhibition of Topo I activity reduced Pol II binding at the LCR and intergenic regions but not at the γ-globin gene promoter. The data demonstrate dynamic associations of transcription factors with the globin gene locus during the cell cycle and support previous results showing that specific components of transcription complexes remain associated with highly transcribed genes during mitosis.

Human genome replication proceeds through four chromatin states

Advances in genomic studies have led to significant progress in understanding the epigenetically controlled interplay between chromatin structure and nuclear functions. Epigenetic modifications were shown to play a key role in transcription regulation and genome activity during development and differentiation or in response to the environment. Paradoxically, the molecular mechanisms that regulate the initiation and the maintenance of the spatio-temporal replication program in higher eukaryotes, and in particular their links to epigenetic modifications, still remain elusive. By integrative analysis of the genome-wide distributions of thirteen epigenetic marks in the human cell line K562, at the 100 kb resolution of corresponding mean replication timing (MRT) data, we identify four major groups of chromatin marks with shared features. These states have different MRT, namely from early to late replicating, replication proceeds though a transcriptionally active euchromatin state (C1), a repressive type of chromatin (C2) associated with polycomb complexes, a silent state (C3) not enriched in any available marks, and a gene poor HP1-associated heterochromatin state (C4). When mapping these chromatin states inside the megabase-sized U-domains (U-shaped MRT profile) covering about 50% of the human genome, we reveal that the associated replication fork polarity gradient corresponds to a directional path across the four chromatin states, from C1 at U-domains borders followed by C2, C3 and C4 at centers. Analysis of the other genome half is consistent with early and late replication loci occurring in separate compartments, the former correspond to gene-rich, high-GC domains of intermingled chromatin states C1 and C2, whereas the latter correspond to gene-poor, low-GC domains of alternating chromatin states C3 and C4 or long C4 domains. This new segmentation sheds a new light on the epigenetic regulation of the spatio-temporal replication program in human and provides a framework for further studies in different cell types, in both health and disease.

Previous studies revealed spatially coherent and biological-meaningful chromatin mark combinations in human cells. Here, we analyze thirteen epigenetic mark maps in the human cell line K562 at 100 kb resolution of MRT data. The complexity of epigenetic data is reduced to four chromatin states that display remarkable similarities with those reported in fly, worm and plants. These states have different MRT: (C1) is transcriptionally active, early replicating, enriched in CTCF; (C2) is Polycomb repressed, mid-S replicating; (C3) lacks of marks and replicates late and (C4) is a late-replicating gene-poor HP1 repressed heterochromatin state. When mapping these states inside the 876 replication U-domains of K562, the replication fork polarity gradient observed in these U-domains comes along with a remarkable epigenetic organization from C1 at U-domain borders to C2, C3 and ultimately C4 at centers. The remaining genome half displays early replicating, gene rich and high GC domains of intermingled C1 and C2 states segregating from late replicating, gene poor and low GC domains of concatenated C3 and/or C4 states. This constitutes the first evidence of epigenetic compartmentalization of the human genome into replication domains likely corresponding to autonomous units in the 3D chromatin architecture.

Robust 3D DNA FISH using directly labeled probes

3D DNA FISH has become a major tool for analyzing three-dimensional organization of the nucleus, and several variations of the technique have been published. In this article we describe a protocol which has been optimized for robustness, reproducibility, and ease of use. Brightly fluorescent directly labeled probes are generated by nick-translation with amino-allyldUTP followed by chemical coupling of the dye. 3D DNA FISH is performed using a freeze-thaw step for cell permeabilization and a heating step for simultaneous denaturation of probe and nuclear DNA. The protocol is applicable to a range of cell types and a variety of probes (BACs, plasmids, fosmids, or Whole Chromosome Paints) and allows for high-throughput automated imaging. With this method we routinely investigate nuclear localization of up to three chromosomal regions.

PfSec13 is an unusual chromatin-associated nucleoporin of Plasmodium falciparum that is essential for parasite proliferation in human erythrocytes

In Plasmodium falciparum, the deadliest form of human malaria, the nuclear periphery has drawn much attention due to its role as a sub-nuclear compartment involved in virulence gene expression. Recent data have implicated components of the nuclear envelope in regulating gene expression in several eukaryotes. Special attention has been given to nucleoporins that compose the nuclear pore complex (NPC). However, very little is known about components of the nuclear envelope in Plasmodium parasites. Here we characterize PfSec13, an unusual nucleoporin of P. falciparum, which shows unique structural similarities suggesting that it is a fusion between Sec13 and Nup145C of yeast. Using super resolution fluorescence microscopy (3D-SIM) and in vivo imaging, we show that the dynamic localization of PfSec13 during parasites' intra-erythrocytic development corresponds with that of the NPCs and that these dynamics are associated with microtubules rather than with F-actin. In addition, PfSec13 does not co-localize with the heterochormatin markers HP1 and H3K9me3, suggesting euchromatic location of the NPCs. The proteins associated with PfSec13 indicate that this unusual Nup is involved in several cellular processes. Indeed, ultrastructural and chromatin immunoprecipitation analyses revealed that, in addition to the NPCs, PfSec13 is found in the nucleoplasm where it is associated with chromatin. Finally, we used peptide nucleic acids (PNA) to downregulate PfSec13 and show that it is essential for parasite proliferation in human erythrocytes.

Insulators target active genes to transcription factories and polycomb-repressed genes to polycomb bodies

Polycomb bodies are foci of Polycomb proteins in which different Polycomb target genes are thought to co-localize in the nucleus, looping out from their chromosomal context. We have shown previously that insulators, not Polycomb response elements (PREs), mediate associations among Polycomb Group (PcG) targets to form Polycomb bodies. Here we use live imaging and 3C interactions to show that transgenes containing PREs and endogenous PcG-regulated genes are targeted by insulator proteins to different nuclear structures depending on their state of activity. When two genes are repressed, they co-localize in Polycomb bodies. When both are active, they are targeted to transcription factories in a fashion dependent on Trithorax and enhancer specificity as well as the insulator protein CTCF. In the absence of CTCF, assembly of Polycomb bodies is essentially reduced to those representing genomic clusters of Polycomb target genes. The critical role of Trithorax suggests that stable association with a specialized transcription factory underlies the cellular memory of the active state.

We have studied the nuclear localization of genes that are regulated by Polycomb mechanisms. The genomes of higher eukaryotes contain hundreds of genes that are regulated by Polycomb mechanisms. Once repressed by Polycomb complexes, they tend to stay repressed; but, when activated, they bind Trithorax protein and tend to maintain the active state epigenetically. Polycomb repression has been reported to make these genes associate in the nucleus to form “Polycomb bodies.” We find that this association is not caused by Polycomb complexes but by insulator elements accompanying the genes. We show that, when these genes are in the active state, the binding of Trithorax targets them to other nuclear regions where transcription occurs, so-called transcription factories. In these nuclear re-positionings the insulator provides the associative power while the state of activity determines whether a gene goes to a Polycomb body or to a transcription factory. The strong effect of Trithorax suggests the possibility that the stable association with a transcription factory it produces may account for the epigenetic memory of the active state.

Transcriptional regulation by histone modifications: towards a theory of chromatin re-organization during stem cell differentiation

Chromatin-related mechanisms, as e.g. histone modifications, are known to be involved in regulatory switches within the transcriptome. Only recently, mathematical models of these mechanisms have been established. So far they have not been applied to genome-wide data. We here introduce a mathematical model of transcriptional regulation by histone modifications and apply it to data of trimethylation of histone 3 at lysine 4 (H3K4me3) and 27 (H3K27me3) in mouse pluripotent and lineage-committed cells. The model describes binding of protein complexes to chromatin which are capable of reading and writing histone marks. Molecular interactions of the complexes with DNA and modified histones create a regulatory switch of transcriptional activity. The regulatory states of the switch depend on the activity of histone (de-) methylases, the strength of complex-DNA-binding and the number of nucleosomes capable of cooperatively contributing to complex-binding. Our model explains experimentally measured length distributions of modified chromatin regions. It suggests (i) that high CpG-density facilitates recruitment of the modifying complexes in embryonic stem cells and (ii) that re-organization of extended chromatin regions during lineage specification into neuronal progenitor cells requires targeted de-modification. Our approach represents a basic step towards multi-scale models of transcriptional control during development and lineage specification.

Regulatory elements associated with paternally-expressed genes in the imprinted murine Angelman/Prader-Willi syndrome domain

The Angelman/Prader-Willi syndrome (AS/PWS) domain contains at least 8 imprinted genes regulated by a bipartite imprinting center (IC) associated with the SNRPN gene. One component of the IC, the PWS-IC, governs the paternal epigenotype and expression of paternal genes. The mechanisms by which imprinting and expression of paternal genes within the AS/PWS domain – such as MKRN3 and NDN – are regulated by the PWS-IC are unclear. The syntenic region in the mouse is organized and imprinted similarly to the human domain with the murine PWS-IC defined by a 6 kb interval within the Snrpn locus that includes the promoter. To identify regulatory elements that may mediate PWS-IC function, we mapped the location and allele-specificity of DNase I hypersensitive (DH) sites within the PWS-IC in brain cells, then identified transcription factor binding sites within a subset of these DH sites. Six major paternal-specific DH sites were detected in the Snrpn gene, five of which map within the 6 kb PWS-IC. We postulate these five DH sites represent functional components of the murine PWS-IC. Analysis of transcription factor binding within multiple DH sites detected nuclear respiratory factors (NRF's) and YY1 specifically on the paternal allele. NRF's and YY1 were also detected in the paternal promoter region of the murine Mrkn3 and Ndn genes. These results suggest that NRF's and YY1 may facilitate PWS-IC function and coordinately regulate expression of paternal genes. The presence of NRF's also suggests a link between transcriptional regulation within the AS/PWS domain and regulation of respiration. 3C analyses indicated Mkrn3 lies in close proximity to the PWS-IC on the paternal chromosome, evidence that the PWS-IC functions by allele-specific interaction with its distal target genes. This could occur by allele-specific co-localization of the PWS-IC and its target genes to transcription factories containing NRF's and YY1.

Chromatin structure and organization: the relation with gene expression during development and disease

The elementary level of chromatin fiber, namely the nucleofilament, is known to undergo a hierarchical compaction leading to local chromatin loops, then chromatin domains and ultimately chromosome territories. These successive folding levels rely on the formation of chromatin loops ranging from few kb to some Mb. In addition to a packaging and structural role, the high-order organization of genomes functionally impacts on gene expression program. This review summarises to which extent each level of chromatin compaction does affect gene regulation. In addition, we point out the structural and functional changes observed in diseases. Emphasis will be mainly placed on the large-scale organization of the chromatin.

Chromosome conformation capture assays in bacteria

Bacterial chromosomes must be compacted by three-orders of magnitude to fit within the cell. While such compaction could in theory yield disordered structures, it is becoming increasingly clear that bacterial chromosomes are in fact arranged in regular and reproducible fashions and that their configurations are tightly connected to fundamental processes such as chromosome segregation. Nonetheless, due to throughput and resolution limitations associated with traditional assays, many question regarding bacterial chromosome structure and its relation to genome function remain. Here, I review the related technologies, chromosome conformation capture (3C) and chromosome conformation capture carbon copy (5C), which my collaborators and I recently introduced as tools to probe the high-resolution folding of entire bacterial genomes. These technologies utilize covalent cross-linking and proximity ligation to facilitate the measurement of the spatial positioning of hundreds of genomic loci, thereby opening the door to high-throughput studies of bacterial chromosome structure. Hence, 3C and 5C represent powerful new tools for assaying the three-dimensional architecture of bacterial genomes.

LINC complexes in health and disease

The cell nucleus communicates with the rest of the cell via nucleo/cytoplasmic transport of proteins and RNA through the nuclear pores. Direct mechanical links between the nucleus and the cytoplasm have recently emerged in the form of LINC (Linkers of the nucleoskeleton to the cytoskeleton) protein complexes. A LINC complex consists of four components. At its core are an inner nuclear membrane (INM) transmembrane protein and an outer nuclear membrane (ONM) transmembrane protein which physically interact with each other in the lumen of the NE. The INM LINC component interacts on the nucleoplasmic side with either the lamina or with an INM-associated protein. The ONM LINC component on the other hand contacts on the cytoplasmatic side a component of the cytoskeleton. This review highlights the components of LINC complexes and their emerging roles in mechanotransduction, nuclear migration, chromosome positioning, signaling, meiosis, cytoskeletal organization and human disease.

Nuclear positional control of HIV transcription in 4D

Retroviruses integrate their genome into the chromatin of the host cell and are subject to the same control mechanisms governing transcription in the nucleus. There is increasing evidence that the spatial position of a gene within the nucleus in time affects its activity. Therefore it becomes important to study the chromatin environment in space and time of the HIV-1 provirus, particularly in cells where a tight transcriptional control allows the virus to hide away from antiviral treatment and immune response. We recently showed that the HIV-1 provirus is found at the nuclear periphery of latently infected lymphocytes associated in trans with centromeric heterochromatin. After induction of transcription, this association was lost, although the location of the transcribing provirus remained peripheral. Our results reveal a novel mechanism of transcriptional silencing involved in HIV-1 post-transcriptional latency and open wider perspectives for the general organization of chromatin in the nucleus.

Super-resolution fluorescence imaging of chromosomal DNA

Super-resolution microscopy is a powerful tool for understanding cellular function. However one of the most important biomolecules - DNA - remains somewhat inaccessible because it cannot be effectively and appropriately labeled. Here, we demonstrate that robust and detailed super-resolution images of DNA can be produced by combining 5-ethynyl-2'-deoxyuridine (EdU) labeling using the 'click chemistry' approach and direct stochastic optical reconstruction microscopy (dSTORM). This method can resolve fine chromatin structure, and - when used in conjunction with pulse labeling - can reveal the paths taken by individual fibers through the nucleus. This technique should provide a useful tool for the study of nuclear structure and function.

ON SIMULATING THE GENERATION OF MOSAICISM DURING MAMMALIAN CEREBRAL CORTICAL DEVELOPMENT

Renewed calls for a systems biology reflect the hope hat enduring biological questions at single-cell and cell-population scales will be resolved as modern molecular biology, with its reductionist program, approaches a nearly-complete characterization of the molecular mechanisms of specific cellular processes. Due to the confounding complexity of biological organization across these scales, computational science is sought to complement the intuition of experimentalists. However, with respect to the molecular basis of cellular processes during development and disease, a gulf between feasible simulations and realistic biology persists. Formidable are the mathematical and computational challenges to conducting and validating cell population-scale simulations, drawn from single-cell level and molecular level details. Nonetheless, in some biological contexts, a focus on core processes crafted by evolution can yield coarse-grained mathematical models that retain explanatory potential despite drastic simplification of known biochemical kinetics.

In this article, we bring this modeling philosophy to bear on the nature of neural progenitor cell decision making during mammalian cerebral cortical development. Specifically, we present the computational component to a research program addressing developmental links between (i) the cellular response to endogenous DNA damage, (ii) primary mechanisms of neuronal genetic heterogeneity, or mosaicism, and (iii) the cell fate decision making that defines the population kinetics of neurogenesis.

Unravelling global genome organization by 3C-seq

Eukaryotic genomes exist in the cell nucleus as an elaborate three-dimensional structure which reflects various nuclear processes such as transcription, DNA replication and repair. Next-generation sequencing (NGS) combined with chromosome conformation capture (3C), referred to as 3C-seq in this article, has recently been applied to the yeast and human genomes, revealing genome-wide views of functional associations among genes and their regulatory elements. Here, we compare the latest genomic approaches such as 3C-seq and ChIA-PET, and provide a condensed overview of how eukaryotic genomes are functionally organized in the nucleus.

Ectopic gene expression and organogenesis in Arabidopsis mutants missing BRU1 required for genome maintenance

Chromatin reconstitution after DNA replication and repair is essential for the inheritance of epigenetic information, but mechanisms underlying such a process are still poorly understood. Previously, we proposed that Arabidopsis BRU1 functions to ensure the chromatin reconstitution. Loss-of-function mutants of BRU1 are hypersensitive to genotoxic stresses and cause release of transcriptional gene silencing of heterochromatic genes. In this study, we show that BRU1 also plays roles in gene regulation in euchromatic regions. bru1 mutations caused sporadic ectopic expression of genes, including those that encode master regulators of developmental programs such as stem cell maintenance and embryogenesis. bru1 mutants exhibited adventitious organogenesis, probably due to the misexpression of such developmental regulators. The key regulatory genes misregulated in bru1 alleles were often targets of PcG SET-domain proteins, although the overlap between the bru1-misregulated and PcG SET-domain-regulated genes was limited at a genome-wide level. Surprisingly, a considerable fraction of the genes activated in bru1 were located in several subchromosomal regions ranging from 174 to 944 kb in size. Our results suggest that BRU1 has a function related to the stability of subchromosomal gene regulation in the euchromatic regions, in addition to the maintenance of chromatin states coupled with heritable epigenetic marks.

Chromatin structure of pluripotent stem cells and induced pluripotent stem cells

Pluripotent embryonic stem (ES) cells are specialized cells with a dynamic chromatin structure, which is intimately connected with their pluripotency and physiology. In recent years somatic cells have been reprogrammed to a pluripotent state through over-expression of a defined set of transcription factors. These cells, known as induced pluripotent stem (iPS) cells, recapitulate ES cell properties and can be differentiated to apparently all cell lineages, making iPS cells a suitable replacement for ES cells in future regenerative medicine. Chromatin modifiers play a key function in establishing and maintaining pluripotency, therefore, elucidating the mechanisms controlling chromatin structure in both ES and iPS cells is of utmost importance to understanding their properties and harnessing their therapeutic potential. In this review, we discuss recent studies that provide a genome-wide view of the chromatin structure signature in ES cells and iPS cells and that highlight the central role of histone modifiers and chromatin remodelers in pluripotency maintenance and induction.

Characterization of the HIV-1 RNA associated proteome identifies Matrin 3 as a nuclear cofactor of Rev function

Background

Central to the fully competent replication cycle of the human immunodeficiency virus type 1 (HIV-1) is the nuclear export of unspliced and partially spliced RNAs mediated by the Rev posttranscriptional activator and the Rev response element (RRE).

Results

Here, we introduce a novel method to explore the proteome associated with the nuclear HIV-1 RNAs. At the core of the method is the generation of cell lines harboring an integrated provirus carrying RNA binding sites for the MS2 bacteriophage protein. Flag-tagged MS2 is then used for affinity purification of the viral RNA. By this approach we found that the viral RNA is associated with the host nuclear matrix component MATR3 (Matrin 3) and that its modulation affected Rev activity. Knockdown of MATR3 suppressed Rev/RRE function in the export of unspliced HIV-1 RNAs. However, MATR3 was able to associate with Rev only through the presence of RRE-containing viral RNA.

Conclusions

In this work, we exploited a novel proteomic method to identify MATR3 as a cellular cofactor of Rev activity. MATR3 binds viral RNA and is required for the Rev/RRE mediated nuclear export of unspliced HIV-1 RNAs.

Early in vitro differentiation of mouse definitive endoderm is not correlated with progressive maturation of nuclear DNA methylation patterns

The genome organization in pluripotent cells undergoing the first steps of differentiation is highly relevant to the reprogramming process in differentiation. Considering this fact, chromatin texture patterns that identify cells at the very early stage of lineage commitment could serve as valuable tools in the selection of optimal cell phenotypes for regenerative medicine applications. Here we report on the first-time use of high-resolution three-dimensional fluorescence imaging and comprehensive topological cell-by-cell analyses with a novel image-cytometrical approach towards the identification of in situ global nuclear DNA methylation patterns in early endodermal differentiation of mouse ES cells (up to day 6), and the correlations of these patterns with a set of putative markers for pluripotency and endodermal commitment, and the epithelial and mesenchymal character of cells. Utilizing this in vitro cell system as a model for assessing the relationship between differentiation and nuclear DNA methylation patterns, we found that differentiating cell populations display an increasing number of cells with a gain in DNA methylation load: first within their euchromatin, then extending into heterochromatic areas of the nucleus, which also results in significant changes of methylcytosine/global DNA codistribution patterns. We were also able to co-visualize and quantify the concomitant stochastic marker expression on a per-cell basis, for which we did not measure any correlation to methylcytosine loads or distribution patterns. We observe that the progression of global DNA methylation is not correlated with the standard transcription factors associated with endodermal development. Further studies are needed to determine whether the progression of global methylation could represent a useful signature of cellular differentiation. This concept of tracking epigenetic progression may prove useful in the selection of cell phenotypes for future regenerative medicine applications.