Genome-wide sequence and functional analysis of early replicating DNA in normal human fibroblasts

Replication Timing Aberration of KIF14 and MDM4 / 2 β Alleles and Aneuploidy as Markers of Chromosomal Instability and Poor Treatment Response in Ewing Family Tumor Patients

Replication timing of allelic gene pairs is strictly regulated according to expression, genome stability, and epigenetic changes, and tumorigenesis may be associated with changes in the allelic replication in various tumors. Our aim was to determine whether such alterations had a prognostic value in Ewing's family tumor (EFT) patients. The KIF14 and MDM4 / PI3KC 2β and the centromeric satellite sequence of chromosomes 8 and 12 were used for replication timing assessments. Aneuploidy was assessed by enumerating the copy numbers of chromosomes 8 and 12. Replication timing and aneuploidy were detected cytogenetically using multicolors fluorescence in situ hybridization assay applied in 135 EFT. Patients with trisomy 8 presented an association with an asynchronous replication pattern (SD) of MDM4 / PI3KC 2β genes ( p  = 0.013). Trisomy 12 was associated with a synchronous pattern (DD) of KIF14 probe signals ( p  = 0.04). The DD synchronous replication pattern of KIF14 showed a correlation with age ( p  < 0.0001), and the SS synchronous replication pattern of the same locus showed a correlation with lung metastatic ( p  = 0.012). The subgroup of patients presenting with multiplet signals of MDM4 / PI3KC 2β showed an association with treatment response ( p  = 0.045) and age ( p  = 0.033). Replication pattern of KIF14 may, significantly, be associated with chromosomal instability as MDM4 / PI3KC 2β may be a considerably new marker of poor treatment response in EFT patients.

Submegabase copy number variations arise during cerebral cortical neurogenesis as revealed by single-cell whole-genome sequencing

Reports of copy number variations (CNVs) within single human brain cells have been limited to megabase-scale alterations in relatively few cells, leaving unclear when CNVs first arise and whether their generation is regulated. Answering these questions has been limited by an absence of experimentation with model organisms that allow developmental assessments infeasible with human samples. Here, we identify the existence and developmental dynamics of cerebral cortical CNVs in mouse, showing that their prevalence increases through midneurogenesis. Our improved sequencing approach also allowed characterization of previously undocumented neural CNVs below 1 Mb in size, comprising half of all alterations. These data demonstrate the existence of myriad CNVs, which genomically diversify neural cells before incorporation into the mature organization of the brain.

Somatic copy number variations (CNVs) exist in the brain, but their genesis, prevalence, forms, and biological impact remain unclear, even within experimentally tractable animal models. We combined a transposase-based amplification (TbA) methodology for single-cell whole-genome sequencing with a bioinformatic approach for filtering unreliable CNVs (FUnC), developed from machine learning trained on lymphocyte V(D)J recombination. TbA–FUnC offered superior genomic coverage and removed >90% of false-positive CNV calls, allowing extensive examination of submegabase CNVs from over 500 cells throughout the neurogenic period of cerebral cortical development in Mus musculus. Thousands of previously undocumented CNVs were identified. Half were less than 1 Mb in size, with deletions 4× more common than amplification events, and were randomly distributed throughout the genome. However, CNV prevalence during embryonic cortical development was nonrandom, peaking at midneurogenesis with levels triple those found at younger ages before falling to intermediate quantities. These data identify pervasive small and large CNVs as early contributors to neural genomic mosaicism, producing genomically diverse cellular building blocks that form the highly organized, mature brain.

Distinct mutational behaviors differentiate short tandem repeats from microsatellites in the human genome

A tandem repeat's (TR) propensity to mutate increases with repeat number, and can become very pronounced beyond a critical boundary, transforming it into a microsatellite (MS). However, a clear understanding of the mutational behavior of different TR classes and motifs and related mechanisms is lacking, as is a consensus on the existence of a boundary separating short TRs (STRs) from MSs. This hinders our understanding of MSs' mutational properties and their effective use as genetic markers. Using indel calls for 179 individuals from 1000 Genomes Pilot-1 Project, we determined polymorphism incidence for four major TR classes, and formalized its varying relationship with repeat number using segmented regression. We observed a biphasic regime with a transition from a faster to a slower exponential growth at 9, 5, 4, and 4 repeats for mono-, di-, tri-, and tetranucleotide TRs, respectively. We used an in vitro mutagenesis assay to evaluate the contribution of strand slippage errors to mutability. STRs and MSs differ in their absolute polymorphism levels, but more importantly in their rates of mutability growth. Although strand slippage is a major factor driving mononucleotide polymorphism incidence, dinucleotide polymorphism incidence is greater than that expected due to strand slippage alone, indicating that additional cellular factors might be driving dinucleotide mutability in the human genome. Leveraging on hundreds of human genomes, we present the first comprehensive, genome-wide analysis of TR mutational behavior, encompassing several motif sizes and compositions.

Temporal and functional analysis of DNA replicated in early S phase

In summary, recently developed technologies have begun to draw back the curtain of mystery that obscures some of the basic mechanisms of DNA replication at multiple levels. Studies using extended DNA and chromatin fiber techniques have proven valuable for identifying the location of origins of replication at specific genomic sites and determining their temporal order of replication, for identifying and quantifying sites of DNA damage and localizing chromatin proteins in relation to sites of DNA replication. The future potential of these methods include further discoveries in functional genomics and contributions to the elucidation of the histone code. Such studies could prove very valuable in studies of the mechanisms of cancer development, aging, and other processes of disordered genomic functioning.

Replication factory activation can be decoupled from the replication timing program by modulating Cdk levels

Cdk activity can differentially regulate the number of active replication factories, replication rates, and the rate of progression through the timing program during S phase.

In the metazoan replication timing program, clusters of replication origins located in different subchromosomal domains fire at different times during S phase. We have used Xenopus laevis egg extracts to drive an accelerated replication timing program in mammalian nuclei. Although replicative stress caused checkpoint-induced slowing of the timing program, inhibition of checkpoint kinases in an unperturbed S phase did not accelerate it. Lowering cyclin-dependent kinase (Cdk) activity slowed both replication rate and progression through the timing program, whereas raising Cdk activity increased them. Surprisingly, modest alteration of Cdk activity changed the amount of DNA synthesized during different stages of the timing program. This was associated with a change in the number of active replication factories, whereas the distribution of origins within active factories remained relatively normal. The ability of Cdks to differentially effect replication initiation, factory activation, and progression through the timing program provides new insights into the way that chromosomal DNA replication is organized during S phase.

Temporal differences in DNA replication during the S phase using single fiber analysis of normal human fibroblasts and glioblastoma T98G cells

We have recently shown that replication forks pause near origins in normal human fibroblasts (NHF1-hTERT) but not glioblastoma T98G cells. This observation led us to question whether other differences in the replication program may exist between these cell types that may relate to their genetic integrity. To identify differences, we detected immunoflourescently the sequential incorporation of the nucleotide analogs IdU and CldU into replicating DNA at the start of every hour of a synchronized S phase. We then characterized the patterns of labeled replicating DNA tracks and quantified the percentages and lengths of the tracks found at these hourly intervals. From the directionality of labeling in single extended replicating DNA fibers, tracks were categorized as single bidirectional origins, unidirectional elongations, clusters of origins firing in tandem, or merging forks (terminations). Our analysis showed that the start of S phase is enriched in single bidirectional origins in NHF1-hTERT cells, followed by an increase in clustering during mid S phase and an increase in merging forks during late S phase. Early S phase in T98G cells also largely consisted of single bidirectional origin initiations; however, an increase in clustering was delayed until an hour later, and clusters were shorter in mid/late S phase than in NHF1-hTERT cells. The spike in merging forks also did not occur until an hour later in T98G cells. Our observations suggest models to explain the temporal replication of single and clustered origins, and suggest differences in the replication program in a normal and cancer cell line.

EpiGRAPH: user-friendly software for statistical analysis and prediction of (epi)genomic data

The EpiGRAPH web service http://epigraph.mpi-inf.mpg.de/ enables biologists to uncover hidden associations in vertebrate genome and epigenome datasets. Users can upload sets of genomic regions and EpiGRAPH will test multiple attributes (including DNA sequence, chromatin structure, epigenetic modifications and evolutionary conservation) for enrichment or depletion among these regions. Furthermore, EpiGRAPH learns to predictively identify similar genomic regions. This paper demonstrates EpiGRAPH's practical utility in a case study on monoallelic gene expression and describes its novel approach to reproducible bioinformatic analysis.

Replication-timing-correlated spatial chromatin arrangements in cancer and in primate interphase nuclei

Using published high-resolution data on S-phase replication timing, we determined the three-dimensional (3D) nuclear arrangement of 33 very-early-replicating and 31 very-late-replicating loci. We analyzed diploid human, non-human primate and rearranged tumor cells by 3D fluorescence in situ hybridization with the aim of investigating the impact of chromosomal structural changes on the nuclear organization of these loci. Overall, their topology was found to be largely conserved between cell types, species and in tumor cells. Early-replicating loci were localized in the nuclear interior, whereas late-replicating loci showed a broader distribution with a higher preference for the periphery than for late-BrdU-incorporation foci. However, differences in the spatial arrangement of early and late loci of chromosome 2, as compared with those from chromosome 5, 7 and 17, argue against replication timing as a major driving force for the 3D radial genome organization in human lymphoblastoid cell nuclei. Instead, genomic properties, and local gene density in particular, were identified as the decisive parameters. Further detailed comparisons of chromosome 7 loci in primate and tumor cells suggest that the inversions analyzed influence nuclear topology to a greater extent than the translocations, thus pointing to geometrical constraints in the 3D conformation of a chromosome territory.

DNA replication in early S phase pauses near newly activated origins

During the S phase of the cell cycle, the entire genome is replicated. There is a high level of orderliness to this process through the temporally and topologically coordinated activation of many replication origins situated along chromosomes. We investigated the program of replication from origins initiating in early S phase by labeling synchronized normal human fibroblasts (NHF1) with nucleotide analogs for various pulse times and measuring labeled tracks in combed DNA fibers. Our analysis showed that replication forks progress 9-35 kilobases from newly initiated origins, followed by a pause in synthesis before replication resumes. Pausing was not observed near origins that initiated in the middle of S phase. No evidence for pausing near origins was found at the beginning of the S phase in glioblastoma T98G cells. Treatment with the S phase checkpoint inhibitor caffeine abrogated pausing in NHF1 cells in early S phase. This suggests that pausing may comprise a novel aspect of the intra-S phase checkpoint pathway or a related new early S checkpoint. Further, it is possible that the loss of this regulatory process in cancer cells such as T98G could be a contributing factor in the genetic instability that typifies cancers.

Early replication and the apoptotic pathway

In higher eukaryotes there is a link between time of replication and transcription. It is generally accepted that genes that are actively transcribed are replicated in the first half of S phase while inactive genes replicate in the second half of S phase. We have recently reported that in normal human fibroblasts there are some functionally related genes that replicate at the same time in S phase. This had been previously reported for functionally related genes that are located in clusters, for example the alpha- and beta-globin complexes. We have shown, however, that this also occurs with some functionally related genes that are not organized in a cluster, but rather are distributed throughout the genome. For example, using GOstat analysis of data from our and other groups, we found an overrepresentation of genes involved in the apoptotic process among sequences that are replicated very early (approximately in the first hour of S phase) in both fibroblasts and lymphoblastoid cells. This finding leads us to question how and why the replication of genes in the apoptotic pathway is temporally organized in this manner. Here we discuss the possible explanations and implications of this observation.

Vascular injury during elevated glucose can be mitigated by erythropoietin and Wnt signaling

Impacting a significant portion of the world's population with increasing incidence in minorities, the young, and the physically active, diabetes mellitus (DM) and its complications affect approximately 20 million individuals in the United States and over 100 million individuals worldwide. In particular, vascular disease from DM may lead to some of the most serious complications that can extend into both the cardiac and nervous systems. Unique strategies that can prevent endothelial cell (EC) demise and elucidate novel cellular mechanisms for vascular cytoprotection become vital for the prevention and treatment of vascular DM complications. Here, we demonstrate that erythropoietin (EPO), an agent that has recently been shown to extend cell viability in a number of systems extending beyond hematopoietic cells, prevents EC injury and apoptotic nuclear DNA degradation during elevated glucose exposure. More importantly, EPO employs Wnt1, a cysteine-rich glycosylated protein involved in gene expression, cell differentiation, and cell apoptosis, to confer EC cytoprotection and maintains the integrity of Wnt1 expression during elevated glucose exposure. In addition, application of anti-Wnt1 neutralizing antibody abrogates the protective capacity of both EPO and Wnt1, illustrating that Wnt1 is an important component in the cytoprotection of ECs during elevated glucose exposure. Intimately linked to this cytoprotection is the downstream Wnt1 pathway of glycogen synthase kinase (GSK-3beta) that requires phosphorylation of GSK-3beta and inhibition of its activity by EPO. Interestingly, inhibition of GSK-3beta activity during elevated glucose leads to enhanced EC survival, but does not synergistically improve protection by EPO or Wnt1, suggesting that EPO and Wnt1 are closely tied to the blockade of GSK-3beta activity. Our work exemplifies an exciting potential application for EPO in regards to the treatment of DM vascular disease complications and highlights a previously unrecognized role for Wnt1 and the modulation of the downstream pathway of GSK-3beta to promote vascular cell viability during DM.

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.