ASAR15, A cis-acting locus that controls chromosome-wide replication timing and stability of human chromosome 15

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

Monoallelically expressed noncoding RNAs form nucleolar territories on NOR-containing chromosomes and regulate rRNA expression

Out of the several hundred copies of rRNA genes arranged in the nucleolar organizing regions (NOR) of the five human acrocentric chromosomes, ~50% remain transcriptionally inactive. NOR-associated sequences and epigenetic modifications contribute to the differential expression of rRNAs. However, the mechanism(s) controlling the dosage of active versus inactive rRNA genes within each NOR in mammals is yet to be determined. We have discovered a family of ncRNAs, SNULs (Single NUcleolus Localized RNA), which form constrained sub-nucleolar territories on individual NORs and influence rRNA expression. Individual members of the SNULs monoallelically associate with specific NOR-containing chromosomes. SNULs share sequence similarity to pre-rRNA and localize in the sub-nucleolar compartment with pre-rRNA. Finally, SNULs control rRNA expression by influencing pre-rRNA sorting to the DFC compartment and pre-rRNA processing. Our study discovered a novel class of ncRNAs influencing rRNA expression by forming constrained nucleolar territories on individual NORs.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

In Vivo Clonal Analysis Reveals Random Monoallelic Expression in Lymphocytes That Traces Back to Hematopoietic Stem Cells

Evaluating the epigenetic landscape in the stem cell compartment at the single-cell level is essential to assess the cells’ heterogeneity and predict their fate. Here, using a genome-wide transcriptomics approach in vivo, we evaluated the allelic expression imbalance in the progeny of single hematopoietic cells (HSCs) as a read-out of epigenetic marking. After 4 months of extensive proliferation and differentiation, we found that X-chromosome inactivation (XCI) is tightly maintained in all single-HSC derived hematopoietic cells. In contrast, the vast majority of the autosomal genes did not show clonal patterns of random monoallelic expression (RME). However, a persistent allele-specific autosomal transcription in HSCs and their progeny was found in a rare number of cases, none of which has been previously reported. These data show that: 1) XCI and RME in the autosomal chromosomes are driven by different mechanisms; 2) the previously reported high frequency of genes under RME in clones expanded in vitro (up to 15%) is not found in clones undergoing multiple differentiation steps in vivo; 3) prior to differentiation, HSCs have stable patterns of autosomal RME. We propose that most RME patterns in autosomal chromosomes are erased and established de novo during cell lineage differentiation.

Three-dimensional genome organization via triplex-forming RNAs

An increasing number of long noncoding RNAs (lncRNAs) have been proposed to act as nuclear organization factors during interphase. Direct RNA-DNA interactions can be achieved by the formation of triplex helix structures where a single-stranded RNA molecule hybridizes by complementarity into the major groove of double-stranded DNA. However, whether and how these direct RNA-DNA associations influence genome structure in interphase chromosomes remain poorly understood. Here we theorize that RNA organizes the genome in space via a triplex-forming mechanism. To test this theory, we apply a computational modeling approach of chromosomes that combines restraint-based modeling with polymer physics. Our models suggest that colocalization of triplex hotspots targeted by lncRNAs could contribute to large-scale chromosome compartmentalization cooperating, rather than competing, with architectural transcription factors such as CTCF.

Differential Allelic Expression among Long Non-Coding RNAs

Long non-coding RNAs (lncRNA) comprise a diverse group of non-protein-coding RNAs >200 bp in length that are involved in various normal cellular processes and disease states, and can affect coding gene expression through mechanisms in cis or in trans. Since the discovery of the first functional lncRNAs transcribed by RNA Polymerase II, H19 and Xist, many others have been identified and noted for their unusual transcriptional pattern, whereby expression from one chromosome homolog is strongly favored over the other, also known as mono-allelic or differential allelic expression. lncRNAs with differential allelic expression have been observed to play critical roles in developmental gene regulation, chromosome structure, and disease. Here, we will focus on known examples of differential allelic expression of lncRNAs and highlight recent research describing functional lncRNAs expressed from both imprinted and random mono-allelic expression domains.

X-Chromosome Inactivation and Autosomal Random Monoallelic Expression as "Faux Amis"

X-chromosome inactivation (XCI) and random monoallelic expression of autosomal genes (RMAE) are two paradigms of gene expression regulation where, at the single cell level, genes can be expressed from either the maternal or paternal alleles. X-chromosome inactivation takes place in female marsupial and placental mammals, while RMAE has been described in mammals and also other species. Although the outcome of both processes results in random monoallelic expression and mosaicism at the cellular level, there are many important differences. We provide here a brief sketch of the history behind the discovery of XCI and RMAE. Moreover, we review some of the distinctive features of these two phenomena, with respect to when in development they are established, their roles in dosage compensation and cellular phenotypic diversity, and the molecular mechanisms underlying their initiation and stability.

Role of mammalian long non-coding RNAs in normal and neuro oncological disorders

Long non-coding RNAs (lncRNAs) are expressed at lower levels than protein-coding genes but have a crucial role in gene regulation. LncRNA is distinct, they are being transcribed using RNA polymerase II, and their functionality depends on subcellular localization. Depending on their niche, they specifically interact with DNA, RNA, and proteins and modify chromatin function, regulate transcription at various stages, forms nuclear condensation bodies and nucleolar organization. lncRNAs may also change the stability and translation of cytoplasmic mRNAs and hamper signaling pathways. Thus, lncRNAs affect the physio-pathological states and lead to the development of various disorders, immune responses, and cancer. To date, ~40% of lncRNAs have been reported in the nervous system (NS) and are involved in the early development/differentiation of the NS to synaptogenesis. LncRNA expression patterns in the most common adult and pediatric tumor suggest them as potential biomarkers and provide a rationale for targeting them pharmaceutically. Here, we discuss the mechanisms of lncRNA synthesis, localization, and functions in transcriptional, post-transcriptional, and other forms of gene regulation, methods of lncRNA identification, and their potential therapeutic applications in neuro oncological disorders as explained by molecular mechanisms in other malignant disorders.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Regulatory and Functional Involvement of Long Non-Coding RNAs in DNA Double-Strand Break Repair Mechanisms

Protection of genome integrity is vital for all living organisms, particularly when DNA double-strand breaks (DSBs) occur. Eukaryotes have developed two main pathways, namely Non-Homologous End Joining (NHEJ) and Homologous Recombination (HR), to repair DSBs. While most of the current research is focused on the role of key protein players in the functional regulation of DSB repair pathways, accumulating evidence has uncovered a novel class of regulating factors termed non-coding RNAs. Non-coding RNAs have been found to hold a pivotal role in the activation of DSB repair mechanisms, thereby safeguarding genomic stability. In particular, long non-coding RNAs (lncRNAs) have begun to emerge as new players with vast therapeutic potential. This review summarizes important advances in the field of lncRNAs, including characterization of recently identified lncRNAs, and their implication in DSB repair pathways in the context of tumorigenesis.

Gene regulation by long non-coding RNAs and its biological functions

Evidence accumulated over the past decade shows that long non-coding RNAs (lncRNAs) are widely expressed and have key roles in gene regulation. Recent studies have begun to unravel how the biogenesis of lncRNAs is distinct from that of mRNAs and is linked with their specific subcellular localizations and functions. Depending on their localization and their specific interactions with DNA, RNA and proteins, lncRNAs can modulate chromatin function, regulate the assembly and function of membraneless nuclear bodies, alter the stability and translation of cytoplasmic mRNAs and interfere with signalling pathways. Many of these functions ultimately affect gene expression in diverse biological and physiopathological contexts, such as in neuronal disorders, immune responses and cancer. Tissue-specific and condition-specific expression patterns suggest that lncRNAs are potential biomarkers and provide a rationale to target them clinically. In this Review, we discuss the mechanisms of lncRNA biogenesis, localization and functions in transcriptional, post-transcriptional and other modes of gene regulation, and their potential therapeutic applications.

Weak interactions in higher-order chromatin organization

The detailed principles of the hierarchical folding of eukaryotic chromosomes have been revealed during the last two decades. Along with structures composing three-dimensional (3D) genome organization (chromatin compartments, topologically associating domains, chromatin loops, etc.), the molecular mechanisms that are involved in their establishment and maintenance have been characterized. Generally, protein-protein and protein-DNA interactions underlie the spatial genome organization in eukaryotes. However, it is becoming increasingly evident that weak interactions, which exist in biological systems, also contribute to the 3D genome. Here, we provide a snapshot of our current understanding of the role of the weak interactions in the establishment and maintenance of the 3D genome organization. We discuss how weak biological forces, such as entropic forces operating in crowded solutions, electrostatic interactions of the biomolecules, liquid-liquid phase separation, DNA supercoiling, and RNA environment participate in chromosome segregation into structural and functional units and drive intranuclear functional compartmentalization.

Reciprocal monoallelic expression of ASAR lncRNA genes controls replication timing of human chromosome 6

DNA replication occurs on mammalian chromosomes in a cell-type distinctive temporal order known as the replication timing program. We previously found that disruption of the noncanonical lncRNA genes ASAR6 and ASAR15 results in delayed replication timing and delayed mitotic chromosome condensation of human chromosomes 6 and 15, respectively. ASAR6 and ASAR15 display random monoallelic expression and display asynchronous replication between alleles that is coordinated with other random monoallelic genes on their respective chromosomes. Disruption of the expressed allele, but not the silent allele, of ASAR6 leads to delayed replication, activation of the previously silent alleles of linked monoallelic genes, and structural instability of human chromosome 6. In this report, we describe a second lncRNA gene (ASAR6-141) on human chromosome 6 that when disrupted results in delayed replication timing in cisASAR6-141 is subject to random monoallelic expression and asynchronous replication and is expressed from the opposite chromosome 6 homolog as ASAR6 ASAR6-141 RNA, like ASAR6 and ASAR15 RNAs, contains a high L1 content and remains associated with the chromosome territory where it is transcribed. Three classes of cis-acting elements control proper chromosome function in mammals: origins of replication, centromeres, and telomeres, which are responsible for replication, segregation, and stability of all chromosomes. Our work supports a fourth type of essential chromosomal element, the "Inactivation/Stability Center," which expresses ASAR lncRNAs responsible for proper replication timing, monoallelic expression, and structural stability of each chromosome.

Replication timing and nuclear structure

DNA replication proceeds along spatially and temporally coordinated patterns within the nucleus, thus protecting the genome during the synthesis of new genetic material. While we have been able to visualize replication patterns on DNA fibers for 50 years, recent developments and discoveries have provided a greater insight into how DNA replication is controlled. In this review, we highlight many of these discoveries. Of great interest are the physiological role of the replication timing program, cis and trans-acting factors that modulate replication timing and the effects of chromatin structure on the replication timing program. We also discuss future directions in the study of replication timing.

Allele-specific control of replication timing and genome organization during development

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

Replication timing kept in LINE

O'Neill and O'Neill discuss Platt et al.’s findings that LINE1 elements are key to control of replication timing by ASAR long noncoding RNAs.

Accurate and synchronous replication timing between chromosome homologues is essential for maintaining chromosome stability, yet how this is achieved has remained a mystery. In this issue, Platt et al. (2018. J. Cell Biol.https://doi.org/10.1083/jcb.201707082) identify antisense LINE (L1) transcripts within long noncoding RNAs as the critical factor in maintaining synchronous chromosome-wide replication timing.

L1 retrotransposon antisense RNA within ASAR lncRNAs controls chromosome-wide replication timing

Proper chromosome duplication is critical for genome integrity and normal cellular function. Platt et al. show that the lncRNA genes ASAR6 and ASAR15 control chromosome-wide replication timing via the antisense strand of L1 retrotransposons located within ASAR6 and ASAR15 RNAs.

Mammalian cells replicate their chromosomes via a temporal replication program. The ASAR6 and ASAR15 genes were identified as loci that when disrupted result in delayed replication and condensation of entire human chromosomes. ASAR6 and ASAR15 are monoallelically expressed long noncoding RNAs that remain associated with the chromosome from which they are transcribed. The chromosome-wide effects of ASAR6 map to the antisense strand of an L1 retrotransposon within ASAR6 RNA, deletion or inversion of which delayed replication of human chromosome 6. Furthermore, ectopic integration of ASAR6 or ASAR15 transgenes into mouse chromosomes resulted in delayed replication and condensation, an increase in H3K27me3, coating of the mouse chromosome with ASAR RNA, and a loss of mouse Cot-1 RNA expression in cis. Targeting the antisense strand of the L1 within ectopically expressed ASAR6 RNA restored normal replication timing. Our results provide direct evidence that L1 antisense RNA plays a functional role in chromosome-wide replication timing of mammalian chromosomes.

Programming asynchronous replication in stem cells

Many regions of the genome replicate asynchronously and are expressed monoallelically. It is thought that asynchronous replication may be involved in choosing one allele over the other, but little is known about how these patterns are established during development. We show that, unlike somatic cells, which replicate in a clonal manner, embryonic and adult stem cells are programmed to undergo switching, such that daughter cells with an early-replicating paternal allele are derived from mother cells that have a late-replicating paternal allele. Furthermore, using ground-state embryonic stem (ES) cells, we demonstrate that in the initial transition to asynchronous replication, it is always the paternal allele that is chosen to replicate early, suggesting that primary allelic choice is directed by preset gametic DNA markers. Taken together, these studies help define a basic general strategy for establishing allelic discrimination and generating allelic diversity throughout the organism.

Monoallelic Gene Expression in Mammals

Monoallelic expression not due to cis-regulatory sequence polymorphism poses an intriguing problem in epigenetics because it requires the unequal treatment of two segments of DNA that are present in the same nucleus and that can indeed have absolutely identical sequences. Here, I focus on a few recent developments in the field of monoallelic expression that are of particular interest and raise interesting questions for future work. One development is regarding analyses of imprinted genes, in which recent work suggests the possibility that intriguing networks of imprinted genes exist and are important for genetic and physiological studies. Another issue that has been raised in recent years by a number of publications is the question of how skewed allelic expression should be for it to be designated as monoallelic expression and, further, what methods are appropriate or inappropriate for analyzing genomic data to examine allele-specific expression. Perhaps the most exciting recent development in mammalian monoallelic expression is a clever and carefully executed analysis of genetic diversity of autosomal genes subject to random monoallelic expression (RMAE), which provides compelling evidence for distinct evolutionary forces acting on random monoallelically expressed genes.

Order from clutter: selective interactions at mammalian replication origins

Mammalian chromosome duplication progresses in a precise order and is subject to constraints that are often relaxed in developmental disorders and malignancies. Molecular information about the regulation of DNA replication at the chromatin level is lacking because protein complexes that initiate replication seem to bind chromatin indiscriminately. High-throughput sequencing and mathematical modelling have yielded detailed genome-wide replication initiation maps. Combining these maps and models with functional genetic analyses suggests that distinct DNA-protein interactions at subgroups of replication initiation sites (replication origins) modulate the ubiquitous replication machinery and supports an emerging model that delineates how indiscriminate DNA-binding patterns translate into a consistent, organized replication programme.

Regulation of Replication Origins

In eukaryotes, genome duplication starts concomitantly at many replication initiation sites termed replication origins. The replication initiation program is spatially and temporally coordinated to ensure accurate, efficient DNA synthesis that duplicates the entire genome while maintaining other chromatin-dependent functions. Unlike in prokaryotes, not all potential replication origins in eukaryotes are needed for complete genome duplication during each cell cycle. Instead, eukaryotic cells vary the use of initiation sites so that only a fraction of potential replication origins initiate replication each cell cycle. Flexibility in origin choice allows each eukaryotic cell type to utilize different initiation sites, corresponding to unique nuclear DNA packaging patterns. These patterns coordinate replication with gene expression and chromatin condensation. Budding yeast replication origins share a consensus sequence that marks potential initiation sites. Metazoan origins, on the other hand, lack a consensus sequence. Rather, they are associated with a collection of structural features, chromatin packaging features, histone modifications, transcription, and DNA-DNA/DNA-protein interactions. These features confer cell type-specific replication and expression and play an essential role in maintaining genomic stability.

A replicator-specific binding protein essential for site-specific initiation of DNA replication in mammalian cells

Mammalian chromosome replication starts from distinct sites; however, the principles governing initiation site selection are unclear because proteins essential for DNA replication do not exhibit sequence-specific DNA binding. Here we identify a replication-initiation determinant (RepID) protein that binds a subset of replication-initiation sites. A large fraction of RepID-binding sites share a common G-rich motif and exhibit elevated replication initiation. RepID is required for initiation of DNA replication from RepID-bound replication origins, including the origin at the human beta-globin (HBB) locus. At HBB, RepID is involved in an interaction between the replication origin (Rep-P) and the locus control region. RepID-depleted murine embryonic fibroblasts exhibit abnormal replication fork progression and fewer replication-initiation events. These observations are consistent with a model, suggesting that RepID facilitates replication initiation at a distinct group of human replication origins.

Origins of mammalian DNA replication are poorly characterised because they lack an Identifiable consensus sequence. Here the authors identify RepID, a protein that binds to a subset of G-rich replication origins and facilitates initiation from those origins.

RNA as a fundamental component of interphase chromosomes: could repeats prove key?

Beginning with the precedent of XIST RNA as a 'chromosomal RNA' (cRNA), there is growing interest in the possibility that a diversity of non-coding RNAs may function in chromatin. We review findings which lead us to suggest that RNA is essentially a widespread component of interphase chromosomes. Further, RNA likely contributes to architecture and regulation, with repeat-rich 'junk' RNA in euchromatin (ecRNA) promoting a more open chromatin state. Thousands of low-abundance nuclear RNAs have been reported, however it remains a challenge to determine which of these may function in chromatin. Recent findings indicate that repetitive sequences are enriched in chromosome-associated non-coding RNAs, and repeat-rich RNA shows unusual properties, including localization and stability, with similarities to XIST RNA. We suggest two frontiers in genome biology are emerging and may intersect: the broad contribution of RNA to interphase chromosomes and the distinctive properties of repeat-rich intronic or intergenic junk sequences that may play a role in chromosome structure and regulation.

Replication timing and transcriptional control: beyond cause and effect-part III

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

The Properties of Long Noncoding RNAs That Regulate Chromatin

Beyond coding for proteins, RNA molecules have well-established functions in the posttranscriptional regulation of gene expression. Less clear are the upstream roles of RNA in regulating transcription and chromatin-based processes in the nucleus. RNA is transcribed in the nucleus, so it is logical that RNA could play diverse and broad roles that would impact human physiology. Indeed, this idea is supported by well-established examples of noncoding RNAs that affect chromatin structure and function. There has been dramatic growth in studies focused on the nuclear roles of long noncoding RNAs (lncRNAs). Although little is known about the biochemical mechanisms of these lncRNAs, there is a developing consensus regarding the challenges of defining lncRNA function and mechanism. In this review, we examine the definition, discovery, functions, and mechanisms of lncRNAs. We emphasize areas where challenges remain and where consensus among laboratories has underscored the exciting ways in which human lncRNAs may affect chromatin biology.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.