Identification of G1-regulated genes in normally cycling human cells

Single-cell nascent RNA sequencing unveils coordinated global transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1,2. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations3. However, fundamental questions about the temporal regulation of transcription and enhancer-gene coordination remain unanswered, primarily because of the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq-a new single-cell nascent RNA sequencing assay that uses click chemistry-and unveil coordinated transcription throughout the genome. We demonstrate the episodic nature of transcription and the co-transcription of functionally related genes. scGRO-seq can estimate burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells and can leverage replication-dependent non-polyadenylated histone gene transcription to elucidate cell cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq enables the identification of networks of enhancers and genes. Our results suggest that the bursting of transcription at super-enhancers precedes bursting from associated genes. By imparting insights into the dynamic nature of global transcription and the origin and propagation of transcription signals, we demonstrate the ability of scGRO-seq to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Single-cell nascent RNA sequencing using click-chemistry unveils coordinated transcription

Transcription is the primary regulatory step in gene expression. Divergent transcription initiation from promoters and enhancers produces stable RNAs from genes and unstable RNAs from enhancers1-5. Nascent RNA capture and sequencing assays simultaneously measure gene and enhancer activity in cell populations6-9. However, fundamental questions in the temporal regulation of transcription and enhancer-gene synchrony remain unanswered primarily due to the absence of a single-cell perspective on active transcription. In this study, we present scGRO-seq - a novel single-cell nascent RNA sequencing assay using click-chemistry - and unveil the coordinated transcription throughout the genome. scGRO-seq demonstrates the episodic nature of transcription, and estimates burst size and frequency by directly quantifying transcribing RNA polymerases in individual cells. It reveals the co-transcription of functionally related genes and leverages the replication-dependent non-polyadenylated histone genes transcription to elucidate cell-cycle dynamics. The single-nucleotide spatial and temporal resolution of scGRO-seq identifies networks of enhancers and genes and indicates that the bursting of transcription at super-enhancers precedes the burst from associated genes. By imparting insights into the dynamic nature of transcription and the origin and propagation of transcription signals, scGRO-seq demonstrates its unique ability to investigate the mechanisms of transcription regulation and the role of enhancers in gene expression.

Disrupting the Cdk9/Cyclin T1 heterodimer of 7SK snRNP for the Brd4 and AFF1/4 guided reconstitution of active P-TEFb

P-TEFb modulates RNA polymerase II elongation through alternative interaction with negative and positive regulation factors. While inactive P-TEFbs are mainly sequestered in the 7SK snRNP complex in a chromatin-free state, most of its active forms are in complex with its recruitment factors, Brd4 and SEC, in a chromatin-associated state. Thus, switching from inactive 7SK snRNP to active P-TEFb (Brd4/P-TEFb or SEC/P-TEFb) is essential for global gene expression. Although it has been shown that cellular signaling stimulates the disruption of 7SK snRNP, releasing dephosphorylated and catalytically inactive P-TEFb, little is known about how the inactive released P-TEFb is reactivated. Here, we show that the Cdk9/CycT1 heterodimer released from 7SK snRNP is completely dissociated into monomers in response to stress. Brd4 or SEC then recruits monomerized Cdk9 and CycT1 to reassemble the core P-TEFb. Meanwhile, the binding of monomeric dephosphorylated Cdk9 to either Brd4 or SEC induces the autophosphorylation of T186 of Cdk9. Finally, the same mechanism is employed during nocodazole released entry into early G1 phase of cell cycle. Therefore, our studies demonstrate a novel mechanism by which Cdk9 and CycT1 monomers are reassembled on chromatin to form active P-TEFb by its interaction with Brd4 or SEC to regulate transcription.

Joint changes in RNA, RNA polymerase II, and promoter activity through the cell cycle identify non-coding RNAs involved in proliferation

Proper regulation of the cell cycle is necessary for normal growth and development of all organisms. Conversely, altered cell cycle regulation often underlies proliferative diseases such as cancer. Long non-coding RNAs (lncRNAs) are recognized as important regulators of gene expression and are often found dysregulated in diseases, including cancers. However, identifying lncRNAs with cell cycle functions is challenging due to their often low and cell-type specific expression. We present a highly effective method that analyses changes in promoter activity, transcription, and RNA levels for identifying genes enriched for cell cycle functions. Specifically, by combining RNA sequencing with ChIP sequencing through the cell cycle of synchronized human keratinocytes, we identified 1009 genes with cell cycle-dependent expression and correlated changes in RNA polymerase II occupancy or promoter activity as measured by histone 3 lysine 4 trimethylation (H3K4me3). These genes were highly enriched for genes with known cell cycle functions and included 57 lncRNAs. We selected four of these lncRNAs-SNHG26, EMSLR, ZFAS1, and EPB41L4A-AS1-for further experimental validation and found that knockdown of each of the four lncRNAs affected cell cycle phase distributions and reduced proliferation in multiple cell lines. These results show that many genes with cell cycle functions have concomitant cell-cycle dependent changes in promoter activity, transcription, and RNA levels and support that our multi-omics method is well suited for identifying lncRNAs involved in the cell cycle.

Involvement of CTCF in transcription regulation of EGR1 at early G1 phase as an architecture factor

Early growth response 1 (EGR1) is a transcription factor and regulates cellular processes such as proliferation, differentiation, and apoptosis. The expression of EGR1 is rapidly induced in response to several stimuli, and it activates the expression of downstream target genes involved in signaling cascades. EGR1 gene is also known to be transcribed in early G1 phase. However, the regulation of EGR1 transcription in early G1 phase is not clarified well. Here we found that CCCTC-binding factor (CTCF), a chromatin binding protein, is required to transcribe EGR1 gene at the onset of early G1 phase. We found that CTCF mediated the formation of higher-order chromatin structures among CTCF binding sites located in the EGR1 locus. Disruption of the CTCF-dependent higher-order chromatin structure using nuclease-dead Cas9 (dCas9)-mediated interference reduced the EGR1 transcription in early G1 phase. Collectively, we propose that CTCF has functional roles for the temporal expression of EGR1 in early G1 phase through regulation of higher-order chromatin structure organization.

Transcriptional Output Transiently Spikes Upon Mitotic Exit

The pulsatile nature of gene activity has recently emerged as a general property of the transcriptional process. It has been shown that the frequency and amplitude of transcriptional bursts can be subjected to extrinsic regulation. Here we have investigated if these parameters were constant throughout the cell cycle using the single molecule RNA FISH technique. We found evidence of transcriptional spikes upon mitotic exit in three different human cell lines. Recording of cell growth prior to hybridization and immuno-RNA FISH analysis revealed that these spikes were short-lived and subsided before completion of cytokinesis. The transient post-mitotic increase in transcriptional output was found to be the result of cells displaying a higher number of active alleles and/or an increased number of nascent transcripts per active allele, indicating that both the burst fraction and the amplitude of individual bursts can be increased upon mitotic exit. Our results further suggest that distinct regulatory mechanisms are at work shortly after mitotic exit and during the rest of interphase. We speculate that transcriptional spikes are associated with chromatin decondensation, a hallmark of post-mitotic cells that might alter the dynamics of transcriptional regulators and effectors.

A hyperactive transcriptional state marks genome reactivation at the mitosis-G1 transition

Hsiung et al. tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. During the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression.

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis–G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer–promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis–G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis–G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis–G1 transition might predispose cells to diverge in gene expression states.

Meta-analysis reveals conserved cell cycle transcriptional network across multiple human cell types

BACKGROUND

Cell division is central to the physiology and pathology of all eukaryotic organisms. The molecular machinery underpinning the cell cycle has been studied extensively in a number of species and core aspects of it have been found to be highly conserved. Similarly, the transcriptional changes associated with this pathway have been studied in different organisms and different cell types. In each case hundreds of genes have been reported to be regulated, however there seems to be little consensus in the genes identified across different studies. In a recent comparison of transcriptomic studies of the cell cycle in different human cell types, only 96 cell cycle genes were reported to be the same across all studies examined.

RESULTS

Here we perform a systematic re-examination of published human cell cycle expression data by using a network-based approach to identify groups of genes with a similar expression profile and therefore function. Two clusters in particular, containing 298 transcripts, showed patterns of expression consistent with cell cycle occurrence across the four human cell types assessed.

CONCLUSIONS

Our analysis shows that there is a far greater conservation of cell cycle-associated gene expression across human cell types than reported previously, which can be separated into two distinct transcriptional networks associated with the G₁/S-S and G₂-M phases of the cell cycle. This work also highlights the benefits of performing a re-analysis on combined datasets.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.

Identification of preferentially reactivated genes during early G1 phase using nascent mRNA as an index of transcriptional activity

During mammalian mitosis, transcription is silenced due to dissociation of transcription factors from DNA and chromosome condensation. At the end of mitosis, transcription is reactivated through chromosome relaxation and reloading of these factors to the DNA. Early G1 genes, which are preferentially reactivated during M/G1 transition, are important for maintenance of cellular function and are known to be strictly regulated. As only few early G1 genes have been identified to date, screening for early G1 genes by genome-wide analysis using nascent mRNA could contribute to the elucidation of the regulatory mechanisms during early G1. Here, we performed a detailed expression analysis for the M/G1 transition of mammalian cells by microarray analysis of nascent mRNA, and identified 298 early G1 genes. Analysis of these genes provides two important insights. Firstly, certain motifs are enriched in the upstream sequences of early G1 genes; from this we could predict candidate cognate transcription factors, including Sp1, which is recruited to the DNA in the early G1 phase. Secondly, we discovered that neighboring genes of early G1 genes were also frequently up-regulated in the G1 phase. Information about these numerous newly identified early G1 genes will likely contribute to an understanding of the regulatory mechanisms of the early G1 genes.

Activation transcription factor-3 activation and the development of spinal cord degeneration in a rat model of amyotrophic lateral sclerosis

It has been reported that an early activation of glial fibrillary acid protein (GFAP) in astroglial cells occurs simultaneously in peripheral nerves and spinal cord from the G93A SOD1 mouse model of amyotrophic lateral sclerosis (ALS), an invariably fatal neurodegenerative disorder. In ALS, the contribute to the pathological process of different cell types varies according to the disease stage, with a florid immune response in spinal cord at end stage disease. In this study, we have mapped in different anatomical sites the process of disease-induced functional perturbation from a pre-symptomatic stage using a marker of cellular distress expressed in neurons and glial cells, the activating transcription factor 3 (ATF-3), and applied large-scale gene expression analysis to define the pattern or transcriptional changes occurring in spinal cord from the G93A SOD1 rat model of ALS in parallel with ATF-3 neuronal activation. From the disease onset onward, transgenic lumbar spinal cord displayed ATF-3 transcriptional regulation and motor cells immunostaining in association with the over-expression of genes promoting cell growth, the functional integrity of cell organelles and involved in the modulation of immune responses. While spinal cord from the pre-symptomatic rat showed no detectable ATF-3 transcriptional regulation, ATF-3 activation was appreciated in large size neurofilament-rich, small size non-peptidergic and parvalbumin-positive neurons within the dorsal root ganglia (DRG), and in ventral roots Schwann cells alongside macrophages infiltration. This pattern of peripheral ATF-3 activation remained detectable throughout the disease process. In the G93A SOD1 rat model of ALS, signs of roots and nerves subtle distress preceded overt clinical-pathological changes, involving both glial cells and neurons that function as receptors of peripheral sensory stimuli from the muscle. In addition, factors previously described to be linked to ATF-3 activation under various experimental conditions of stress, become switched on in spinal cord from the end-stage transgenic rat model of ALS.