Continuous transcription initiation guarantees robust repair of all transcribed genes and regulatory regions

Liakos, Anastasios; Konstantopoulos, Dimitris; Lavigne, Matthieu D.; Fousteri, Maria

doi:10.1038/s41467-020-14566-9

Download PDF

Article
Open access
Published: 14 February 2020

Continuous transcription initiation guarantees robust repair of all transcribed genes and regulatory regions

Nature Communications volume 11, Article number: 916 (2020) Cite this article

6614 Accesses
20 Citations
20 Altmetric
Metrics details

Subjects

Abstract

Inhibition of transcription caused by DNA damage-impaired RNA polymerase II (Pol II) elongation conceals a local increase in de novo transcription, slowly progressing from Transcription Start Sites (TSSs) to gene ends. Although associated with accelerated repair of Pol II-encountered lesions and limited mutagenesis, it is still unclear how this mechanism is maintained during genotoxic stress-recovery. Here we uncover a widespread gain in chromatin accessibility and preservation of the active H3K27ac mark after UV-irradiation. The concomitant increase in Pol II escape from promoter-proximal pause (PPP) sites of most active genes, PROMPTs and enhancer RNAs favors unrestrained initiation, as evidenced by the synthesis of nascent RNAs including start RNAs. Accordingly, drug-inhibition of PPP-release replenishes levels of pre-initiating Pol II at TSSs after UV. Our data show that such continuous engagement of Pol II molecules ensures maximal transcription-driven repair throughout expressed genes and regulatory loci. Importantly, revealing this unanticipated regulatory layer of UV-response provides physiological relevant traction to the emerging concept that Pol II initiation rate is determined by pause-release dynamics.

A CSB-PAF1C axis restores processive transcription elongation after DNA damage repair

Article Open access 26 February 2021

Diana van den Heuvel, Cornelia G. Spruijt, … Martijn S. Luijsterburg

Active mRNA degradation by EXD2 nuclease elicits recovery of transcription after genotoxic stress

Article Open access 20 January 2023

Jérémy Sandoz, Max Cigrang, … Frédéric Coin

Elongation factor ELOF1 drives transcription-coupled repair and prevents genome instability

Article 09 June 2021

Marit E. Geijer, Di Zhou, … Jurgen A. Marteijn

Introduction

Pol II initiation at TSSs and release into productive elongation from PPP are ubiquitous and crucial steps regulating transcription of protein-coding genes and long non-coding RNAs^1,2 (together called mRNAs in this paper). Similarly, Pol II is regulated for transcribing regulatory non-coding regions: enhancer RNAs (eRNAs) are expressed bidirectionally from eTSSs^3,4,5, while PROMoter uPstream transcripts (or upstream antisense RNAs, collectively called PROMPTs here) are produced in the opposite direction to mRNA when two stable transcripts are not initiated in very close proximity and in opposite directions (bidirectional promoters)⁶. However, contrary to mRNAs, eRNAs and PROMPTs are short and unstable due to high early termination rates and increased susceptibility to degradation by the RNA exosome^6,7, thus making their detection technically challenging.

Initiation of transcription in all the above regions depends on the efficient assembly of the pre-initiation complex (PIC) upstream of TSS and serine 5 phosphorylation (S5P) of Pol II C-terminal domain (CTD)⁸. After elongation of ~30–60 nucleotides of initiation-associated RNAs (start-RNAs)^7,9, Pol II is paused at PPP sites by negative elongation factors DSIF and NELF^2,10. Signal-regulated phosphorylation of these factors and of Pol II CTD serine 2 (S2P) by P-TEFb is required for productive elongation^11,12,13. It recently emerged that, if this step does not occur rapidly, start-RNAs are terminated^12,14, implying that Pol II turnover at PPP sites is high at steady state, and that replenishment of Pol II engaged in early transcription (initiation to PPP) is achieved through the continuous re-entry of pre-initiating Pol II into PICs^14,15.

The integrity of the genetic information encoded in DNA sequence is persistently challenged by a variety of genotoxic perturbations¹⁶. A plethora of DNA damage response (DDR) mechanisms have evolved to guarantee the detection and removal of different types of DNA lesions, limiting the probability of mutagenesis by adjusting to the cell’s status and need for efficient recovery from DNA damage^16,17,18. Nucleotide excision repair (NER) plays a vital role in sensing and removing a large panel of helix-distorting DNA adducts such as cyclobutane pyrimidine dimers (CPDs) induced by ultraviolet (UV) light, and benzo[a]pyrene guanine adducts induced by cigarette smoke¹⁷. Transcription-coupled NER (TC-NER) is promptly triggered by elongating Pol II molecules encountering DNA adducts and speeds up excision and repair in expressed loci^19,20,21. In comparison, the second NER subpathway, Global Genome-NER (GG-NER), operates through the entire genome, but recognizes more stochastically helix distortions^17,21,22,23. Importantly, given all the classes of transcripts defined above, it is estimated that the coverage of transcribed regions²⁴ potentially scanned by TC-NER expands to more than 50% of the genome, thus qualifying transcription as a major driving force in safeguarding genomic stability.

Although TC-NER depends on lesion-sensing potential by elongating Pol II molecules, transcription elongation has been shown to be transiently inhibited after UV irradiation^25,26,27 due to a proportion of Pol II molecules stalling at encountered DNA damages^25,28. Moreover, depletion of the pre-initiating hypo-phosphorylated Pol II(hypo) isoform from chromatin shortly after UV irradiation^25,29,30 has led to the assumption that new transcription initiation events are transiently and globally repressed^{21,29,30,31,32,33}. On the other hand, recent reports^25,26,34 have revealed a functionally essential stress-dependent increase in 5′ nascent RNA (nRNA) activity globally that depends on the UV-induced increase in active P-TEFb levels^35,36, and on the rapid dissociation of the NELF complex³⁷. The ensuing fast and global release of de novo Pol II elongation waves from PPP sites into gene bodies boosts lesion-sensing activity and accelerates removal of DNA adducts by TC-NER in virtually all active mRNA genes²⁵. Together, these findings substantiate the possibility that initiation of transcription might not be as severely affected by UV as previously believed.

Taking also into consideration recent evidence that supports the model of disengagement of a given Pol II molecule from DNA template after damage recognition^21,31,38, it is tempting to assume that ensuring continuity in transcription initiation may bring advantages in the repair process. We thus hypothesized that the apparent loss of pre-initiating RNAPII may not be due to the absence of RNAPII recruitment at TSSs, but rather due to a decrease in the dwell time of Pol II-hypo isoform at TSS, as justified by the concomitant increase in Pol II-ser5P and Pol II-ser2P downstream of TSS²⁵. In this way, cells would be able to uninterruptedly feed the global release of scanning enzymes into transcribable sequences, and guarantee the detection of more lesions along genes’ template strand (TS). Herein, we decipher chromatin dynamics genome-wide upon UV damage, and find a significant gain in accessibility (ATAC-seq) at the TSSs of virtually all active regulatory regions controlling mRNAs, PROMPTs, and eRNA expression. This phenomenon is underlined by the maintenance of active histone marks (H3K27ac), the lack of deposition of transcriptional silencing modifications (H3K27me3) at transcribed loci, and is correlated with the influx of Pol II into productive elongation. The paradoxical decrease in pre-initiating Pol II-hypo at these TSSs upon UV is elucidated by revealing that the presence of Pol II-hypo is rescued when PPP release is drug-inhibited. Accordingly, preserved production of start-RNAs after UV stress lies under the increased production of nRNA, and is prevented only after inhibition of transcription initiation. The identified genome-wide dependence of initiation rate on promoter-proximal pause-release dynamics explains the seamless recruitment/initiation of Pol II upon UV, in turn enabling efficient repair of the totality of the sequences encoding active regulatory regions and mRNAs.

Results

Chromatin accessibility increases at active regulatory regions upon UV

To characterize the impact that UV irradiation might have on the chromatin landscape of transcriptional regulatory regions during the early recovery times (from 0.5 to 4 h), and how this could be linked to the widespread PPP release of elongating Pol II and the local increase in nRNA production downstream of TSS^25,26,34, we first determined the genome-wide changes in chromatin accessibility. The omni-ATAC-seq protocol³⁹ was implemented in our system involving UV-C irradiation (15 J/m²) of human skin fibroblasts synchronized in G1 (see “Methods” and also ref. ²⁵). We reproducibly measured chromatin accessibility before (NO UV) and 2 h after (+UV) irradiation (Supplementary Fig. 1a, b), performed peak calling and mapped a total of 106,052 accessible regions (ARs) across conditions (Supplementary Fig. 1b), and combined the replicates (see “Methods” for details). ARs were enriched at promoters and intragenic or intergenic regions with transcriptional regulatory function (TSSs, TSS flanks, and enhancers according to ChromHMM annotation, Fig. 1a–c; Supplementary Fig. 1c–e, “Methods”). Interestingly, we reveal a widespread increase (up to 1.71 average fold change (FC)) in chromatin accessibility after stress at 97.9% of promoter-, 94.6% of intragenic-, and 94.4% of intergenic ARs (Fig. 1b–d; Supplementary Fig. 1e, f).

**Fig. 1: Increase in chromatin accessibility in response to mild doses of UV irradiation.**

We then selected differentially accessible regions (DARs) by applying stringent thresholds both in terms of FC (Log₂ FC > 1) and P value (P < 0.001), and found that 6410 loci showed particularly increased chromatin accessibility upon UV (DAR-gain) (Fig. 1e, top panel). DAR-gain found at promoter regions represented 13.3% of all promoter ARs (Fig. 1e, lower panel), thus pinpointing towards a potentially functionally relevant chromatin opening at TSS regions. DAR-gain located at intragenic and intergenic loci (Fig. 1e) were linked to genes if they overlapped functional enhancers defined in FANTOM5 (see “Methods”). We found that genes associated with DAR-gain loci (either identified on their promoter or enhancers) were representative (adjusted (adj.) P < 0.05) of a number of biological pathways previously associated with DDR processes, including cellular response to stress, DNA repair, transcription regulation by TP53, and cell cycle checkpoints (Supplementary Fig. 2). In addition, we identified a broad range of many other significant Gene Ontology (GO) categories (163 in total, Supplementary Fig. 2), a result in line with the previously reported global PPP release of elongating Pol II waves at all active gene bodies upon UV irradiation²⁵.

Chromatin marks linked to transcription remain stable after UV

A number of studies have demonstrated that the turnover, modification, and/or degradation of histones around damage sites represent essential steps in conserved pathways that help cells deal with genotoxic stress^40,41. However, especially in the case of UV-C-induced DNA damage, little is known about the post-translational modifications (PTMs) of histones around transcriptional regulatory regions. To better interpret the increase in chromatin accessibility and clarify its possible impact on genome-wide transcription dynamics, we studied the differential presence of two histone PTMs representative for the transcription status of associated chromatin: the silencing mark H3K27me3 and the activation mark H3K27ac^42,43.

We conducted ChIP-seq experiments with antibodies specific for both H3K27ac and H3K27me3 in NO UV and +UV (+2 h) conditions, and focused our analysis on TSSs of mRNAs and on a robust set of eTSSs, which are known to be functional and potentially transcribed in the investigated cell type according to the FANTOM5 database (see “Methods”). We used the ChIP-seq data generated in this study (H3K27ac and H3K27me3), as well as previously published ChIP-seq data (Pol II-ser2P²⁵) from the steady state (NO UV) condition, to determine subsets of active (presence of H3K27ac and Pol II-ser2P peaks over TSS), repressed (presence of H3K27me3 peaks over TSS), and inactive loci (no peak detected over TSS for H3K27ac, H3K27me3, and Pol II) (Fig. 2a, see “Methods”) in our cell system. We associated the changes in histone marks and Pol II observed in these regions upon UV with ATAC-seq results. The increase in chromatin accessibility was detected at all active TSSs, which correspond largely to the promoters identified above (compare Figs. 1a and 2a, and see “Methods”), as well as FANTOM5-annotated active eTSSs upon UV (Fig. 2b, ATAC, 95% confidence interval (CI) excludes 0). This opening was in sharp contrast to the UV-induced global loss of Pol II-hypo at TSSs and eTSSs (Fig. 2b, Pol II-hypo, 95% CI excludes 0) observed 1.5 h after UV irradiation (8 J/m²) (Fig. 2a, b; Supplementary Fig. 3).

**Fig. 2: Histone modifications remain virtually stable upon UV damage.**

Strikingly, we also found preservation (slight but not significant increase) of H3K27ac levels (Fig. 2a, b, 95% CI includes 0; Supplementary Fig. 3), and we observed no exchange of H3K27ac for H3K27me3 in response to UV at these active TSSs and eTSSs. Reciprocally, there was no loss of H3K27me3 for H3K27ac, and no gain of Pol II at repressed loci (Fig. 2a; Supplementary Fig. 3). Accordingly, the results of our genome-wide analysis were consistent with biochemical evidence obtained by histone acetic extraction followed by western blot analysis, showing that the global levels of H3K27me3 or H3K27ac remain fairly stable during the early period of recovery from UV stress (15 J/m²) (Supplementary Fig. 3c, d).

We therefore conclude that depletion of detectable Pol II-hypo at TSSs and eTSSs does not occur due to repression of these loci by trimethylation of H3K27⁴³, or because of the loss of the activating histone mark H3K27ac⁴².

Chromatin opening parallels Pol II transition into elongation

To elucidate the functional advantage associated with increased chromatin accessibility in response to UV, we performed a thorough integrative analysis of our data in relation to previously published data sets (Pol II-ser2P from ref. ²⁵ and Cap Analysis Gene Expression (CAGE)-seq from ref. ⁴, see “Methods”). First, we customized a genome annotation, which unambiguously pinpoints to the TSSs of mRNAs, PROMPTs, and eRNAs that do not overlap with regions possibly being transcribed through from neighboring/overlapping genes, promoters, or enhancers (see “Methods”). We then established three categories (Fig. 3a–c), as per previously suggested models⁴⁴: first, active bidirectional promoter regions, which include the TSSs of mRNA–mRNA pairs transcribed in opposite directions (Fig. 3a); second, active unidirectional promoters, which include the TSS of a given mRNA gene (+ or −) for which we could associate an expressed PROMPT in the antisense direction (Fig. 3b); third, active intergenic—as opposed to intragenic—enhancers to avoid potential contamination by interfering reads that derive from overlapping transcription of other active elements (Fig. 3c). Importantly, PROMPT and enhancer transcriptional activity were defined from available CAGE data for the skin and dermal fibroblasts (FANTOM5 consortium, see “Methods”) that accurately determine transcript starting position (5′ end), abundance, and directionality of Pol II transcription in our model (Fig. 3a–c, CAGE). TSS loci were sorted by inter-TSS distance, which we defined as the distance separating TSSs and/or the summits of CAGE signals detected on the (+) and (−) strands (Fig. 3a, b; “Methods”). This allowed us to identify regions with overlapping (convergent, CONV) or non-overlapping (divergent, DIV) transcription (Fig. 3a, b). By focusing on the latter category, we could study the dynamics of transcription at play only in each direction, without having to deal with potential interferences from overlapping regions.

**Fig. 3: Release of Pol II from pausing sites in coding and non-coding transcribed regions upon UV stress.**

Using this setup, we discovered that the UV-dependent increase in chromatin accessibility (Fig. 3a–c, ATAC) was paralleled to the transition of Pol II into active elongation (Fig. 3a–c, Pol II-ser2P), not only at flanking mRNAs (Fig. 3d, e) but also at adjacent PROMPTs and eRNA sequences (Fig. 3e, f), as shown by the loss in Pol II reads at TSSs and the gain of reads in downstream regions. These results were confirmed quantitatively by showing that escape index (EI) of elongating Pol II increased in the +UV condition in comparison with NO UV for 90.1% of bidirectional promoters (Fig. 3g, Chi-square test P = 5.1 × 10⁻²⁶⁶), as well as for 70.1% of PROMPTs (Fig. 3h, Chi-square test P = 4.5 × 10⁻⁸⁹) and 68.6% of eRNAs (Fig. 3i, Chi-square test P = 2.5 × 10⁻⁴⁴). We conclude that the PPP release of Pol II upon genotoxic stress is synchronously triggered at all active transcription units, and coincides with increased chromatin breathing. These data extend the previously characterized transcription-driven genome surveillance mechanism²⁵ to essentially all active gene-regulatory regions, and give mechanistic insights into the synergy between the increase in chromatin accessibility and the transcriptional response observed upon UV.

DRB rescues post-UV detection of Pol II in PIC

We noted that although 63.65% of the transcribed genome shows reduction in transcription activity (coverage of the transcriptome with Log₂ FC (+UV/NO UV)) < 0, see “Methods”) in line with earlier published data⁴⁵, a local increase in nRNA synthesis downstream of TSS of all active genes is detected during the UV-recovery phase^25,26,27,34. This observation combined with the above findings on the UV-induced chromatin opening around virtually all active TSSs, PROMPTs, and eTSSs are hardly compatible with the previously suggested model of UV-induced global inhibition of transcription initiation. We thus searched for alternative reasons that could explain reduction of Pol II-hypo levels at active TSSs/eTSSs despite increased accessibility after UV.

We performed a set of experiments aiming to determine whether Pol II was actually recruited to TSSs upon UV (Pol II-hypo isoform used as a proxy). First, as depicted in Fig. 4a, we irradiated cells with a UV (15 J/m²) and left them to recover for 2 h, when the levels of Pol II-hypo are known to be severely depleted^25,29,30 (Supplementary Fig. 4a, b (DMSO NO UV vs DMSO +UV +2 h)). We then applied, or not, an inhibitor of Pol II release into elongation from PPP sites (DRB, see “Methods”). Cells were crosslinked (X) 2 h after the addition of DRB (or DMSO for the control cells). In accordance with the above mentioned previous reports, in cells that were crosslinked 2 h after UV irradiation in the absence of DRB (+UV/X 2 h), or in cells that were crosslinked 4 h after UV irradiation, and had been incubated with DMSO for the last 2 h (+UV/−DRB/ X 4 h), we detected only minimal levels of pre-initiating Pol II in total chromatin extracts or at TSSs, PROMPTs, and eTSSs, as revealed by Western blot analysis (Fig. 4b) and ChIP-seq (Fig. 4c, d), respectively. In contrast, when cells had been incubated with DRB for the last 2 h before being crosslinked at 4 h after UV irradiation (+UV/+DRB/X 4 h), we observed a significant rescue of pre-initiating Pol II (hypo) levels in total chromatin (Fig. 4b, two-sided Student’s t test P = 0.0055 compared with “+UV/-DRB /X 4 h” and P = 0.0156 compared with “+UV/X 2 h”). The restoration of pre-initiating Pol II levels was even more pronounced when we focused on the occupancy on active TSSs, PROMPTs, and eTSSs, where average read densities detected by Pol II-hypo ChIP-seq after DRB treatment (+UV/+DRB/X 4 h) matched the control NO UV levels (NO UV/+DRB/X 4 h) (Fig. 4c, d). Therefore, even by blocking the stress-triggered transition of Pol II molecules from PPP sites into elongation at 2 h post UV, when the prior-to-UV Pol II-hypo levels were almost completely depleted, we were able to reveal the underlying continuous de novo recruitment of Pol II-hypo molecules in PICs.

**Fig. 4: DRB rescues pre-initiating Pol II levels following DNA damage.**

We also applied DRB just before and for 2 h after UV (15 J/m², Supplementary Fig. 4a), and found a limited loss of pre-initiating Pol II in chromatin extracts upon UV (Supplementary Fig. 4b, c, two-sided Student’s t test P = 0.0145). This result was corroborated by ChIP-qPCR experiments (performed on the same chromatin extracts used above), as DRB prevented the UV-induced reduction in occupancy of Pol II-hypo at promoter/TSS-proximal regions of six active genes (Supplementary Fig. 4d, two-sided Student’s t test P = 0.003 for DMSO, while P = 0.2876 (non-significant) for DRB).

We thus conclude that the genome-wide UV-induced PPP release of Pol II molecules into elongation accelerates the transition into initiation of the next-to-be-recruited Pol II-hypo molecules, limiting the dwell time of this isoform at essentially all active TSSs, PROMPTs, and eTSSs.

Increased RNA synthesis from active TSSs upon UV

As UV irradiation does not inhibit neither the recruitment of Pol II-hypo into PICs nor Pol II escape into elongation, we next verified that newly synthesized RNA molecules were detected at the beginning of transcribed regions. We took advantage of our and other nRNA-seq data^25,27, and we examined if the previously characterized global increase in EU- or Bru-labeled RNA reads at the beginning of genes (see Supplementary Fig. 4 in ref. ²⁵) could originate from increased Pol II initiation at active TSSs (Fig. 5a, b; Supplementary Fig. 5a, b), as suggested before²⁷. In particular, at unidirectional promoters, we confirmed that nRNA synthesis was increased in the mRNA direction, but we also found a concomitant increase of nRNA production in the antisense, PROMPT direction. Similarly, we found widespread gains in eRNAs synthesis, which emanate equally in both directions from active eTSSs (Fig. 5a, b; Supplementary Fig. 5a, b). Identifying labeled nRNA even at short transcripts such as PROMPTs and eRNAs confirms active labeling close to TSSs, and validates the fact that regions directly downstream of TSSs get de novo transcribed during the post-UV period. Taken together, these data demonstrate that the continuous recruitment of Pol II-hypo molecules (see Fig. 4) and their fast transition into initiation/productive elongation (see Fig. 3), during the recovery from genotoxic stress, is accompanied by de novo synthesis of RNA directly downstream of TSSs.

**Fig. 5: De novo and increased RNA synthesis from all TSSs upon UV exposure.**

To further verify initiation activity during UV recovery, we exploited the possibility to track start-RNAs, which directly inform on the amount of dynamically engaged Pol II located within the initially transcribed sequence (approximately the first 100 nucleotides⁷). We followed the experimental procedure depicted in Fig. 5c, and applied, or not, transcription elongation (DRB) or initiation (triptolide, TRP) inhibitors 2 h post UV (15 J/m²). For each condition, we isolated small RNAs by size selection (<200 nucleotides), and we ligated an RNA–DNA linker to their 3′ ends. Reverse transcription (RT) was performed using a universal primer annealing to the linker sequence as previously described⁷. Subsequently, locus-specific qPCR reactions were performed in order to compare, in a quantitative way, the levels of start-RNAs at representative active loci for which we had identified Pol II-ser2P ChIP-seq or nRNA-seq signal (see “Methods”). Our results revealed that start-RNAs could be detected after UV treatment, validating the fact that initiation still occurs during the UV-recovery phase (Fig. 5d, +UV/−DRB). Similar results were obtained in the presence of the transcription elongation inhibitor (Fig. 5d, +UV/+DRB). However, the opposite was found after inhibiting transcription initiation by TRP, which as expected led to a clear reduction of start-RNAs (Fig. 5d, +UV/+TRP, two-sided Student’s t test, P = 0.0037 compared with “NO UV/+DRB”, P = 0.0016 compared with “+UV/−DRB”, and P = 0.0009 compared with “+UV/+DRB”). Together, these results consolidate further the evidence of the non-stop recruitment and initiation of Pol II at TSSs after UV irradiation.

Equal levels of Pol II-hypo at PICs prime for uniform TC-NER

Next, we took advantage of XR-seq data (eXcision-Repair sequencing)²⁰, which precisely and exclusively pinpoint the location and levels of transcription-dependent repair (TC-NER pathway) when the assay is performed in GG-NER-deficient cells (xeroderma pigmentosum (XP)-C cells). Given the strand specificity of the assay, we considered only the excision of CPD damages from template (non-coding) strand (TS) for mRNAs, PROMPTs, and eRNAs, which corresponded to the + (blue) or the − (red) strand of the genome (Fig. 6a, b), depending on the transcript orientation. Upon correlation with CAGE, we found that onset of TC-NER coincided with CAGE reads location, confirming the fact that TC-NER (triggered by damage-arrested Pol II molecules²⁰) and CAGE⁴ accurately locate active TSSs (Fig. 6b (compare with Fig. 3a–c, left), Fig. 6c, d). As expected, repair efficiency was equal in each direction for bidirectional active promoters (Fig. 6b–e). This result was also in line with Pol II-hypo ChIP-seq data showing equivalent amounts of Pol II recruitment at PICs in both directions (Supplementary Fig. 6a, b), and CAGE data indicating strand-balanced production of capped mRNAs (Fig. 6e, CAGE, boxes centered around Log₂ FC = 0). Nevertheless, we note that the variability between both directions was strikingly less for TC-NER and Pol II-hypo than for CAGE (Fig. 6e (proportion of non-significant F-Tests: P = 0), Supplementary Fig. 6b (top panel, proportion of non-significant F-Tests: P = 0)).

**Fig. 6: TC-NER is homogeneous at all transcribed regions.**

Next, we further investigated repair of PROMPTs and enhancers, a phenomenon previously observed, but hardly explained^20,46. We quantified strand-specific repair upstream and downstream of unidirectional promoters, and found that repair activity at unambiguously resolved divergent PROMPTs was stronger than expected from CAGE levels (Fig. 6b–d, DIV). Indeed, XR-seq read density was not correlated to the steady state levels of CAGE at those loci (Pearson correlation coefficient (PCC) = 0.1343). Also, FCs of TC-NER reads between mRNAs and PROMPTs were much smaller than those for CAGE (Fig. 6f, 95% CI excludes 0), thus matching the UV-independent Pol II-hypo uniformity (Supplementary Fig. 6a, b). Similarly, TC-NER levels on TS of eRNAs were higher than anticipated. Indeed, the densities of eRNA XR-seq reads were similar to those of mRNAs (Fig. 6b–d), and contrasted with the very low CAGE signal detected at these loci (Fig. 3a–c, left). Therefore, balanced Pol II-hypo loading in PICs at all classes of transcripts, during steady state or upon stress, allows for equal initiation events and mirrors the homogeneous levels of XR-seq detected in these regions. Taken together, our results demonstrate that the widespread continual initiation and release into productive elongation of Pol II waves maximize repair activity, regardless of prior-to-UV transcript expression level at all kinds of active regulatory regions (mRNAs, PROMPTs, and enhancers).

Continuous initiation drives TC-NER to completion

We next assessed the biological relevance of continuous transcription initiation from active regulatory regions during the UV-recovery period. We have reported previously²⁵ that in the absence of a UV-triggered PPP release of elongating Pol II waves, the scanning activity of pri-elongating (e.g., already elongating prior to UV) Pol II molecules is not sufficient to enable recognition of damaged thymidine dimers (TpT sequence, abbreviated as TT) in the totality of the transcribed genome. In accordance, a recent study reports that inhibition of transcription elongation with DRB limits the extent of excision activity detected by CPD XR-seq³⁸. Thus, sending Pol II molecules from TSS into gene bodies after UV is of pivotal importance to maximize the detection of TTs on active TS.

We quantified the effect of various regimes of DRB (see scheme in Fig. 7a) on the dynamics and distribution of TC-NER (CPD XR-seq data from XP-C cells³⁸) around TTs and around TSSs of mRNA genes, PROMPTs, and enhancers (Fig. 7a–d, respectively, see “Methods” and ref. ²⁵). We observe that when DRB was applied immediately after UV, repair of TTs located from TSS to upstream of the +1 -h wave front (WF) of UV-released Pol II (as obtained in the non-DRB control) are seriously affected (Fig. 7a–c, DRB + UV + 1 h, clusters 0-II).

**Fig. 7: TC-NER activity monitored with various regimes of DRB reveals the need for continuous transcription initiation.**

Accordingly, allowing only a limited number of pri-elongating molecules to be launched just before the irradiation, while preventing de novo Pol II release after stress (Fig. 7a–c, DRB2 + UV + 1 h) highlights how an insufficient feed in Pol II impairs TC-NER activity at both distal and proximal transcribed regions (compare signal before and after asterisk positions in Fig. 7 a, c).

We reasoned that the high extent of ongoing repair activity detected without DRB (UV + DMSO) is the result of the continuity in Pol II initiation. To confirm this directly, we mapped XR-seq reads around TTs located on TS 1, 4, and 8 h after UV in WT cells (data obtained from ref. ⁴⁶, see “Methods”). We found that significant levels of transcription-dependent excision activity were maintained at lesions located directly downstream of active genes TSSs at 4 and 8 h into the recovery process (compare Fig. 8a–c and Supplementary Fig. 7a–c). Notably, it also appeared that the bulk of the excision activity on TS of active genes shifts overtime (+8 h) from the proximal to the distal part of long genes (Fig. 8a–c, clusters III–VI; Supplementary Fig. 7d).

**Fig. 8: Continuous initiation maintains TC-NER efficiency during the whole recovery period and across the whole transcribed region.**

Overall, these results reveal that a large extent of the transcription-driven repair activity detected after UV is due to the ongoing entry of Pol II molecules at TSSs, which can scan farther and farther lesions along the gene bodies. Our analysis highlights the advantage of a continuously supplied transcription-dependent repair process (Fig. 8 and Supplementary Fig. 7d, TS of active genes) over slower CPD detection capabilities of GG-NER, which was detected at significantly lower levels across all regions on NTS of active genes and on both strands of inactive genes (Supplementary Fig. 7d).

Discussion

In this study, we provide quantitative insights into the molecular processes underlying the major transcription-coordinated cellular response that is activated in human cells upon genotoxic stress^{25,26,27,34,35,37}. The establishment of precise maps of chromatin state helped us to query in detail the impact of transcription on DNA repair activities at important functional regions, including PROMPT and eRNA loci. Our results support a model of continuous transcription initiation that can feed the widespread UV-triggered escape of Pol II into the elongation, enabling long-lasting efficient DNA lesion scanning of the whole transcribed genome.

The finding that an increase in chromatin accessibility parallels the conservation of H3K27ac-modified nucleosomes at the flanks of already opened regions in response to UV irradiation is compatible with reports showing that there can be a significant gain in DNA accessibility without changes in nucleosome occupancy during rapid transcriptional induction⁴⁷. Notably, the maintenance of H3K27ac at these sites prevents the imposition of repressive trimethylation at active loci (see Fig. 2), in accordance with the rule that H3K27ac and H3K27me3 are mutually exclusive⁴⁸. Moreover, finding more active transcription at these loci complies with prior reports, suggesting that an increase in gene expression is associated with surges in chromatin accessibility⁴⁹, and that the presence of nRNA inhibits the recruitment of H3K27me3-catalyzing polycomb-repressive complex 2 (PRC2) at active genes⁵⁰. Our results obtained using a range of mild UV doses (8–20 J/m²), and focusing on the early phase of recovery (1-4 h), support the fact that chromatin relaxation and maintenance of H3K27ac has an active role in the repair of transcription-blocking lesions, and give substance to previously observed rather low increase in acetylation of histones in similar conditions⁵¹. In contrast, drastic chromatin remodeling observed in murine cells at a later time during recovery (6 h) when much higher doses (80 J/m²) were used⁵² fits with findings signifying that cellular response to UV depends on the exposure dose^53,54. In other words, when cells deal with unmanageable levels of damages, they need to implement more radical expression changes, which are required for the associated fate of programmed death^55,56, a protective mechanism limiting the risk of malignant transformation.

Our analysis takes advantage of a high-resolution strand-specific map of TSSs for coding and non-coding (enhancers and PROMPTs) loci, and supports the idea that bidirectional transcription of divergent RNAs arises from two distinct hubs of transcription initiation (PICs), located within a single nucleosome-depleted region (NDR)^8,57,58. Indeed, for bidirectional mRNAs and mRNA-PROMPTs, the binding of Pol II-hypo occurs at both edges of highly accessible regions (see ATAC-seq vs Pol II-hypo in Supplementary Fig. 6a, c), which correspond to single NDRs flanked by H3K27ac nucleosomes (see arrows in Supplementary Fig. 6a, c). These observations also extend the evidence supporting the claim that enhancers and PROMPT PICs are organized in a similar manner to gene PICs^8,13. Finally, the observed differences in transcript levels between PROMPTs and mRNAs (see Fig. 5) are probably not due to differences in Pol II-hypo recruitment (see Supplementary Fig. 6a, b, bottom), but rather due to differences in the frequency of premature termination at PPP sites and/or differences in degradation of PROMPT RNAs by the RNA exosome, which is known to be inhibited upon UV stress^59,60.

By uncoupling TSSs of mRNA genes from those of PROMPTs and enhancers, we reveal that P-TEFb-dependent release of elongating Pol II from PPP sites extends to all actively transcribed regions (see Fig. 3). Interestingly, a growing number of studies have reported data to suggest that (i) UV irradiation preferentially inhibits elongation, rather than transcription initiation^25,26,27,34, (ii) P-TEFb and NELF are important regulators of UV response^35,37,61, and (iii) although elongation gradually decelerates due to the encounter of Pol II with DNA lesions, significant initiation/early elongation activity is observed in the first thousand bases of actively transcribed regions^25,26,34, a characteristic that has also been used for the identification of active TSSs genome-wide after UV²⁷. These features are consistent with our finding that new Pol II molecules are constantly recruited to PICs post-UV (see Fig. 4), suggesting that a pool of non-engaged polymerases is still available, and that they promptly proceed into initiation of start-RNAs and subsequently into elongation of longer nRNAs (see Fig. 5). Considering that Pol II ChIP-seq density depends, among others, on the epitope residence time at a given genomic locus⁶², and that Pol II molecules recruited in the PIC are readily phosphorylated upon UV^25,29,30, we propose that the rapid exchange of Pol II isoforms after UV irradiation represents a perfectly plausible cause for the decreased ability to detect Pol II-hypo molecules at TSSs upon UV. It also accounts for the increase in EI, and could explain the gain in accessibility around TSS (see Figs. 2–4, Supplementary Fig. 6, and ref. ²⁵). Such a model explains previously published data concerning the persistence of sufficiently important levels to be worthy of attention, of (i) PIC/basal transcription factors in nuclear extracts³⁰ or upstream of genes’ TSSs (TFIIB)³³ and (ii) nRNAs at the beginning of genes^25,26,27,34 upon UV.

We note that excision fragments (from XR-seq) are distributed more homogeneously at sense (mRNA) and antisense (PROMPT) strands of unidirectional TSSs, and at enhancers, than it could be predicted from the CAGE levels (see Fig. 6). This finding reinforces the possibility that efficient repair at stable and unstable transcripts is primed by the uniform recruitment of Pol II-hypo at all classes of PICs in the steady state (see Supplementary Fig. 6). Remarkably, UV induces continuous and uniform transition into initiation (see Fig. 4), and constantly feeds a long-lasting PPP release of lesion-sensing transcribing Pol II into the transcriptome (see Figs. 3 and 5).

This concept was further validated by re-analyzing an experiment mapping excision repair of CPDs (XR-seq) after inhibition of PPP release³⁸, and thus a fortiori preventing the beneficial effects of Pol II initiation. Indeed, upon DRB treatment, we find a drastic impairment of CPD excisions at the beginning of all classes of active transcripts (see Fig. 7d), and in particular at TTs located within the distance normally covered by de novo released Pol II molecules, without inhibition, at 1 h post-UV (low XR-seq signal for DRB vs DMSO in cluster 0-II, see Fig. 7a, c). In addition, the decrease in the percentage of excision in cluster II (between 10 kb and WF) in DRB2 experiment (see Fig. 7 a–c, DRB2 vs DMSO) indicates that continuity in initiation post UV is crucial for the repair of sequences located toward the 3′ end of longer genes. Taking into consideration that only one Pol II molecule can be accommodated per PPP site at each actively transcribed loci at the time of irradiation, applying DRB after UV mainly restricts the benefits related to the continuous recruitment, initiation, and release into elongation of multiple polymerases throughout the recovery period.

Critically, such a scheme is compatible with the evidence that continuous nRNA-seq signal can be detected uninterruptedly between 2 and 12 h post UV from the start to the end of genes²⁶. We propose that keeping initiation active maximizes the probability to repair quickly all lesions located on the TS. Our model is compatible with the idea that even if a Pol II molecule is ubiquitylated and degraded after the recruitment of the NER complex^21,31,38, the amount of Pol II still normally recycling after termination from all the short genes (which contain less or no lesions) will provide sufficient trailing polymerases for scanning lession-containing genes. In this way, this mechanism can sense efficiently one after the other the more distal lesions even in longer genes, and gives an alternative biological perspective to presently favored mode of action²¹. In agreement, we detect TC-NER activity on TS directly downstream of TSSs 4 and 8 h after UV. In addition, repair of distal damages in long genes is detected only when DRB is not applied after UV (compare clusters III–VI with DRB2 and DMSO in Fig. 7a–c), and intensifies at 8 h in normal conditions (see Fig. 8 clusters III-VI; Supplementary Fig. 7d, distal).

Increase in TC-NER at regulatory regions has also been observed in E. coli⁶³, and is compatible with the idea that the act of antisense transcription over regulatory regions exerts a meaningful biological function⁶⁴ conserved through evolution. Indeed, these DNA sequences may serve as binding sites for transcription factors, or encode target sites for RNA-binding proteins, enabling accurate regulation of topologically associated mRNA genes^58,65. Given the effect of DNA repair on the landscape of somatic mutations in cancer tissues^17,18, surveillance of these vital sequences impacts on cell’s fitness. We propose that our model could account for the low levels of substitutions recently observed upstream of genes’ TSSs and around DNAse-hypersensitive (DHS) sites^66,67. As such, transcription of non-coding regulatory loci could serve to keep their transcription factor-binding sequences (TFBS) in check.

Recent advances in the field of transcription regulation point to the fact that activation of paused genes is mediated through switching from a premature termination state of Pol II at PPP sites to a processive elongation state^11,14,15, implying that continuous cycles of initiation followed by rapid premature termination are required for fast transcriptional induction¹⁴. Our results, showing that persistent initiation guarantees a prolonged transcription-coupled NER across the whole transcriptome, are functionally linked to the fact that DNA damage-triggered widespread PPP release of a given Pol II is sufficient to drive immediate initiation of the next Pol II (see Figs. 4 and 5; Supplementary Figs. 4 and 5). In other words, the clearance rate of Pol II from TSSs highly depends on PPP status. These findings are favorable to the emerging concept that Pol II promoter-proximal pausing has an inhibitory effect on initiation^68,69,70, and highlight how this mechanism can function across the whole transcriptome. At the same time, our results provide a compelling physiological relevance to why cells could gain from firing initiation continuously, as a balance between promoter-proximal termination and escape into elongation allows efficient dynamic responses to stimuli or genotoxic stress.

Methods

Cell culture and treatments

Cells used in this study were VH10 hTERT-immortalized human skin fibroblasts, and were cultured, synchronized by low-serum starvation, and released in full medium as described previously²⁵, unless stated differently. When applied, 5,6-dichloro-1-β-D-ribofuranosylbenzimidazole (DRB, Calbiochem) and triptolide (called TRP, Invivogen) were used in a final concentration of 100 μΜ and 125 nM, respectively, and they were added directly in growth media at the indicated times. Cells were irradiated with UV-C (254 nm, TUV Lamp, Philips) (8–20 J/m², as indicated) and left to recover for indicated times (see figure legends and below).

ChIP-seq

ChIP-seq was performed as previously described²⁵ with minor changes (see Supplementary Methods for details). Cells were mock-treated (NO UV) or treated with UV (+UV) with 15 J/m² except if otherwise stated. Treated cells were left to recover for 2 h before harvesting or as indicated on the timeline (Fig. 4).

The antibodies used for ChIP were the following: H3K27ac (ab4729, Abcam), H3K27me3 (07-449, Millipore), and 8WG16 (Pol II-hypo) (05-952, Millipore). The primers used for ChIP-qPCR experiments were the following (5′ to 3′, F: forward, R: reverse, ChIA neg was the negative primer): SSBP1_F: GTGAGGGAGGAAGGGATAGC, SSBP1_R: AGGGCCAGACACCTACACAG, OSBPL9_F: ATTGGCGGCTCCCAAGAT, OSBPL9_R: GCATTGTAGTCCAGCACGAA, TRPM7_F: CCCAGGGAAACCTTCTCAG, TRPM7_R: TCGCACAATTATGAAAGACTCG, MYC_F: ACTCAGTCTGGGTGGAAGGTATC, MYC_R: GGAGGAATGATAGAGGCATAAGGAG, AKNA_F: CCGTTCCAATCCCTTACC, AKNA_R: TGGAACAAAGAATTCACAGG, APRT_F: GCCTTGACTCGCACTTTTGT, APRT_R: TAGGCGCCATCGATTTTAAG, ChIA_neg_F: AGTCTGAGCTTTGTGGACAGC, and ChIA_neg_R: CCCTCCCAGTATACAGTCTTGC. qPCR, library preparation, and next-generation sequencing were performed as previously described²⁵. Values of all qPCR replicates are supplied in the Source Data file.

Western blot analysis

Western blot analysis of equal amounts of crosslinked chromatin extracts or of histone extracts (see Supplemetary Methods) was performed as described²⁵. Antibodies used for western blot analysis are the following: anti-H3K27ac (ab4729, Abcam), anti-H3K27me3 (07-449, Millipore), 8WG16 (05-952, Millipore), anti-elongating RNA pol II (ab5095, Abcam), anti-Lamin B1 (ab65986, Abcam), anti-histone 4 (ab10158, Abcam), and anti-histone3 (ab1791, Abcam). Dilutions of antibodies were performed according to the manufacturer's guidelines. Time for analysis is indicated in the figures (Fig. 4b; Supplementary Fig. 4a, b; Supplementary Fig. 3c, d). Uncropped scans of all western blot figures are supplied in the Source Data file.

Assay for transposase-accessible chromatin (ATAC)-seq

ATAC-seq method (nuclei preparation, transposition, and amplification of transposed fragments for library preparation) was performed using Nextera DNA Library Prep Kit (Illumina, Inc.) and primers as described by Corces et al.³⁹ with minor modifications: (i) 70,000 cells were used per experimental condition and (ii) the DNase treatment of cells in culture medium, before the transposition reaction, was skipped. The UV dose applied for ATAC-seq experiments was 15 J/m², and treated cells were left to recover for 2 h before harvesting.

Start RNA isolation and qPCRs

To isolate small RNAs (smaller than 200 nucleotides), we used Qiagen miRNeasy Mini Kit and RNeasy MinElute Cleanup Kit according to the manufacturer's instructions. In order to monitor the efficiency of the different enzymatic reactions, we included in our experiments a spike-in RNA oligonucleotide of known sequence (oGAB11: rArGrUrCrArCrUrUrArGrCrGrArUrGrUrArCrArCrUrGrArCrUrGrUrG, synthesized and purified by IDT). After purification, small RNAs and spike-in molecules were ligated to the IDT DNA linker 1 (/5rApp/CTGTAGGCACCATCAAT/3ddC/). Specifically, samples were denatured for 2 min at 80 °C and then placed immediately on ice. Ligation mix (4.8 μl of 50% PEG, 2 μl of 10× RNA ligase buffer, linker and RNase-free H₂O, and 0.5 μl of truncated RNA ligase) (NEB, Cat No. M0351S) was added in a final volume of 20 μl. The reaction was incubated for 3 h at 37 °C. After H₂O was added to a final volume of 100 μl, ethanol precipitation (three volumes of 100% EtOH, with 1/10th volume of 3M NaAc, pH 5.2, and 10 μg of glycogen (ThermoFisher Scientific, Cat No. AM9510) was performed overnight at −80 °C. RNA was purified at 10 μl, and reverse transcription (RT) was performed using primer oLSC003: /5Phos/TCGTATGCCGTCTTCTGCTTG/iSp18/CACTCA/iSp18/AATGATACGGCGACCACCGATCCGACGATCATTGATGGTGCCTACAG according to Invitrogen Superscript II (Cat No. 18064014) instructions. qPCR was performed using gene-specific forward primers (sequences 5′ to 3′ for OSBPL9: ATTGGCGGCTCCCAAGAT, SSBP1: GTGAGGGAGGAAGGGATAGC, IFIT1: TCTCAGAGGAGCCTGGCTAA, and KPNA6: ATTTGGCGAGAGCCTGTCT) and one common reverse primer (oNTI230: 5′-AATGATACGGCGACCACCGA-3′), which anneals to RT primer oLSC003 sequence. Quantitative PCR results were obtained from two independent biological experiments. Values of all qPCR replicates are supplied in the Source Data file.

Read alignment, normalization, peak calling, and differential accessibility analysis

For all next-generation sequencing (NGS) data analyses, in-house scripts and pipelines were developed to automate and analyze the data consistently (see below for details). Code is available upon request. Sequenced data and generated wig profiles are available on Gene Expression Omnibus (GEO) (Accession ID: GSE125181). ATAC-seq and ChIP-seq reads were subjected to quality control, data filtering, and alignment, and wig profile generation was performed essentially as described previously²⁵ with minor modifications. Chip-seq data for Pol II-ser2P and Pol II-hypo (Figs. 2 and 3; Supplementary Figs. 3, 6) were obtained from ref. ²⁵ for NO UV and 2 or 1.5 h post UV (8 J/m²), respectively. nRNA-seq data were generated in our lab previously²⁵ (GSE83763) and obtained from ref. ²⁷ (GSE75398). Downloaded data were processed as described in ref. ²⁵ (see Fig. 5; Supplementary Fig. 5).

For H3K27ac and H3K27me3 ChIP-seq alignment files, peak calling was performed using SICER version 1.1⁷¹ with window parameter = 400 bp and gap parameter = 1, while false discovery rate (fdr) and log₂ fold change cutoffs were set to 0.01 and 1.5, respectively.

For ATAC-seq alignment files, peak calling was performed using MACS2⁷². Because of the variability of ATAC-seq fragment lengths, several runs of the peak-calling algorithm were performed, using different parameters per run, in an attempt to maximize the sensitivity of the detection of open-chromatin regions. In particular, --nomodel --shift 100 --extsize 200, --broad --shift 100 --extsize 200, --nomodel --shift 37 --extsize 73, --broad --shift 37 --extsize 73, --broad --nomodel --shift 37 --extsize 73 --keep-dup all, --broad --nomodel --shift 100 --extsize 200 --keep-dup all and --nomodel --shift 75 --extsize 150 --keep-dup all runs were combined, and detected peaks were filtered using fdr < 0.05 and fold change > 1. Only peaks present in five out of seven methods were kept for further analysis. Although the majority of peaks detected were common between conditions, a number of peaks were also detected only in NO UV and +UV. As these peaks were less in number than the common ones, and showed a rather low density of reads in another attempt to maximize the sensitivity of the detection of open-chromatin regions, and to unbiasedly investigate the changes occurring upon irradiation, we considered the union of the peaks in each condition to perform the rest of the analysis (Supplementary Table 1). To conduct differential accessibility analysis, diffBind R package (https://www.bioconductor.org/packages//2.10/bioc/html/DiffBind.html) was used, with the merged ATAC-seq peak set as a reference. Differential accessibility regions were detected and filtered by applying fold change (Log₂ FC ≥ 1) and P-value (p-val ≤ 0.001) thresholds.

Read-density plots

ATAC-seq, ChIP-seq, nRNA-seq, CAGE-seq, and XR-seq data were subjected to read-density analysis after read depth normalization of all samples per experiment. Heatmaps and average-density profiling were computed as described previously²⁵ around genomic regions of interest, as indicated in the figures. Heatmaps were generated directly using the software, from matrices of binned read densities (bin size is indicated in the figures) for all considered individual (n) items (metagenes). Read- density matrices were also imported in R and python custom scripts for (i) plotting average-density profiles (smoothing achieved by a moving window of the bin size as indicated) and (ii) for determination of read densities per genomic category.

Construction of mRNA-TSSs, PROMPT-TSSs, and eTSS annotation

To annotate transcription start sites (TSSs), all known protein-coding and non-coding RNA hg19 RefSeq transcripts release 86 was downloaded from UCSC table browser (http://genome-euro.ucsc.edu/cgi-bin/hgTables). For each transcript, a biotype was assigned using BioMart (www.biomart.org), and all the small non-coding RNAs were excluded. For all the gene models containing multiple alternative transcripts, TSS neighborhoods of a 100-bp window were clustered together, and only the longest transcript was kept, resulting in 30,473 transcripts. Transcripts were then separated into three groups, based on their transcriptional activity. TSS coordinates were extended to 2 kb on each direction, and were tested for overlap with the Pol II-ser2P-UV, H3K27ac-UV, and H3K27me3-UV peak sets. Regions overlapping with at least one Pol II-ser2P-UV and H3K27ac -UV peak were characterized as active, those overlapping with a H3K27me3-UV peak, but not with a Pol II-ser2P-UV or with a H3K27ac-UV peak, were characterized as repressed, and those that did not overlap with any of the above peak sets were characterized as inactive. Any region overlapping with both H3K27ac-UV and H3K27me3-UV peaks was excluded from the rest of the analysis. This resulted in 15,819 active, 2943 repressed, and 7608 inactive transcripts (Supplementary Table 1). To further classify the active TSSs in terms of transcription directionality, the annotation was split up into unidirectional and bidirectional references. All active transcript pairs with opposite direction of transcription (forward(+) vs reverse(−)), where −2 kb ≤ TSS_distance ≤ +2 kb, TSS_distance = TSS _{coordinate_forward strand} − TSS_{coordinate reverse strand} (inter-TSS distance), were characterized as bidirectional, while the rest of the annotations were characterized as unidirectional. Bidirectional pairs were further categorized into two groups of annotations: convergent bidirectional transcript pairs with TSS_distance ≤ 100 bp, and divergent bidirectional transcript pairs with TSS_distance> 100 bp. To optimize the categorization of convergent and divergent transcript pairs, TSS coordinates were redefined by scanning at a radius of 250 bp, to detect the nucleotide occupied by the maximum-sense CAGE signal. Any bidirectional pair with a non-significant CAGE peak in the aforementioned region was excluded from the analysis. This finally resulted in 12,859 unidirectional transcripts and 2822 active bidirectional TSS pairs, 1806 of which were characterized as divergent and 1016 as convergent (Supplementary Table 1).

To gain a complete overview of the non-coding antisense transcription events occurring around mRNA-TSSs, we also annotated upstream antisense (uaRNA) and downstream antisense (daRNA) transcripts (referred as an ensemble to PROMPTs in this paper for convenience). Only the active unidirectional mRNA-TSSs were used. For all the genes annotated with more than one mRNA transcript, only the leftmost TSS (for + strand genes), and rightmost TSS (for − strand genes) were considered for the rest of the analysis. The antisense CAGE peak with the highest summit in the region that ranged from −2 kb upstream to +1 kb downstream of each unidirectional TSS was considered to be the main PROMPT TSS for further analyses (inter-TSS distance = mRNA TSS - CAGE PROMPT). The above procedure was also repeated for the inactive transcript set to estimate the highest CAGE summit background distribution. The putative active PROMPT CAGE summits, which were higher than the average of the summit background distribution, were considered as active. This resulted in 5366 pairs of active unidirectional—PROMPT-TSSs, which were categorized to 1444 divergent and 3922 convergent pairs, as described above (Supplementary Table 1). By focusing on the divergent loci, the dynamics of transcription could be studied at play in each direction, without having to deal with interference from either direction. Therefore, analysis was focused on upstream antisense RNA, which corresponds to the original definition of PROMPTs⁶.

To annotate enhancer transcription start sites (eTSSs), all 65,423 human enhancers from phases 1 and 2 of the FANTOM5 project from http://fantom.gsc.riken.jp/5/datafiles/phase2.2/extra/Enhancers/human_permissive_enhancers_phase_1_and_2.bed.gz, and the center of each annotation, were considered as the corresponding transcription start site. Enhancers were separated into 6766 active, 4730 repressed, and 39,227 inactive following the pipeline described above. Active intergenic enhancers were further analyzed, and all the eTSSs within a distance of 10 kb from nearby active transcripts, or neighbor eTSSs within a distance of 2 kb, were excluded. The rest of the intergenic eTSSs were extended to 1 kb in both directions, and sense and antisense maximum CAGE summit heights were detected for each reference. This procedure was also repeated for the inactive enhancer set, and inactive sense and antisense highest CAGE summit background distributions were estimated as described above. Finally, the putative active intergenic sense and antisense CAGE summits, which were higher averages of the summit background distributions, were considered as active. This resulted in 1228 active intergenic eTSSs (Supplementary Table 1).

Promoter escape indices analysis

Promoter escape analysis was performed for a subset of active unidirectional and bidirectional transcripts, PROMPTs and active enhancers. In particular, to avoid the inclusion of Pol II-ser2P reads mapped in overlapping promoters and gene bodies, only active divergent unidirectional transcript—PROMPT pairs were considered, where TSS_distance> 100 bp, TSS_distance = TSS _{coordinate forward reference} − TSS _{coordinate reverse reference}, active divergent bidirectional transcript pairs with TSS_distance> 100 bp, and active intergenic enhancers with no nearby transcripts within 10 kb and no nearby eTSSs within 2 kb. TSSs and PROMPT-TSS promoter escape indexes (EI, inverse of pausing index) were calculated as previously defined²⁵, by taking the average coverage in rpm in the gene body (density in gene body was abbreviated as Db and ranged from 101 bp to 2 kb downstream of TSS or 101 bp downstream of TSS to TTS for genes larger or smaller than 2 kb, respectively) divided by the average coverage on the promoter-proximal region (Dp) ranged from 250 bp upstream to 100 bp downstream of TSS.

For enhancer escape analysis, EI was calculated as above, where density of reads at enhancer flanks (Df) is calculated for the regions ranging from −2 kb to −100 bp upstream of eTSS and from +100 bp to +2 kb downstream of eTSS, while density of reads on enhancer TSS (De) is calculated for the regions ranging from 100 bp upstream to 100 bp downstream of eTSS.

Nucleotide excision repair data meta-analysis

The strand-specific genome-wide maps of nucleotide excision repair of the UV-induced DNA damage (CPDs), available for XP-C mutants lacking the global genome nucleotide excision repair mechanism (GG-NER-deficient, TC-NER-proficient), were obtained from Hu et al.²⁰ (UV: 20 J/m², data used in Fig. 6, Gene Expression Omnibus (GEO) accession number GSE67941) and Chiou et al.³⁸ (UV: 20 J/m² data used in Fig. 7, GEO accession number GSE106823). XR-seq data for wild-type (WT) cells (UV: 10 J/m² used for Fig. 8 and Supplementary Fig. 7, GEO accession number GSE76391) were obtained from Adar et al.⁴⁶. Sequence read archive (SRA) data sets were downloaded from Gene Expression Omnibus using the sra toolkit prefetch (https://www.ncbi.nlm.nih.gov/sra/docs/sradownload/) command, and converted to fastq files using fastq-dump. Fastq quality control, data filtering, and short read alignment were performed as above. Meta-analysis involved that read counts were normalized to equal read depth. Heatmap read-density matrices and average read-density plots of CPD XR-seq read around potential pyrimidine dimers. For XP-C data, we focus on TpT sequences (TT) as explained previously²⁵, and filtered the TTs overlapping with enhancer regions defined above to avoid signals generated from eTSSs. For WT cells, the list of TTs was further curated to only consider TTs in active genes (or inactive genes) located on plus strand and filtering out TTs in the region between TSS and +2 kb for CONV and DIV bidirectional and unidirectional genes with inter-TSS < 100 bp (Supplementary Table 1). XR-seq reads were computed as described above. Read-density matrices were calculated for both strands (TS and NTS) separately when indicated. For WT cells, only minus reads were used as these correspond to excised DNA from TS for plus genes. The ratio of XR-seq reads between directions, and calculation of variability between directions was performed as described in the legend of Fig. 6. “S–F” scores and quantification of reads around TT loci was performed as described in Lavigne et al.²⁵, with cluster borders defined previously.

FANTOM5 Cap-analysis of gene expression (CAGE) sequencing data meta-analysis

The FANTOM5 strand-specific CAGE-seq alignment files of normal dermal fibroblast primary cells (six donors with source codes: 11269-116G9, 11346-117G5, 11418-118F5, 11450-119A1, 11454-119A5, and 11458-119A9) and normal skin fibroblasts (two donors with source codes: 11553-120C5 and 11561-120D4) were downloaded from ftp://ftp.biosciencedbc.jp/archive/fantom5/datafiles/phase2.2/basic/human.primary_cell.hCAGE and were combined. Heatmap read-density matrices and average read-density plots were computed as described in the section “Read densities heatmaps and average plots”. Read-density matrices were calculated for both strands separately.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

The data reported in this paper have been deposited with the Gene Expression Omnibus under accession code GSE125181.

References

Haberle, V. & Stark, A. Eukaryotic core promoters and the functional basis of transcription initiation. Nat. Rev. Mol. Cell Biol. 19, 621–637 (2018).
Article CAS PubMed PubMed Central Google Scholar
Jonkers, I. & Lis, J. T. Getting up to speed with transcription elongation by RNA polymerase II. Nat. Rev. Mol. Cell Biol. 16, 167–177 (2015).
Article CAS PubMed PubMed Central Google Scholar
Lai, F. & Shiekhattar, R. Enhancer RNAs: the new molecules of transcription. Curr. Opin. Genet Dev. 25, 38–42 (2014).
Article CAS PubMed PubMed Central Google Scholar
The FANTOM Consortium. et al. An atlas of active enhancers across human cell types and tissues. Nature 507, 455–461 (2014).
Article ADS CAS Google Scholar
Kim, T.-K. et al. Widespread transcription at neuronal activity-regulated enhancers. Nature 465, 182–187 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Preker, P. et al. RNA exosome depletion reveals transcription upstream of active human promoters. Science 322, 1851–1854 (2008).
Article ADS CAS PubMed Google Scholar
Henriques, T. et al. Stable pausing by RNA polymerase II provides an opportunity to target and integrate regulatory signals. Mol. Cell 52, 517–528 (2013).
Article CAS PubMed Google Scholar
Lai, W. K. M. & Pugh, B. F. Genome-wide uniformity of human ‘open’ pre-initiation complexes. Genome Res. 27, 15–26 (2017).
Article CAS PubMed PubMed Central Google Scholar
Nechaev, S. et al. Global analysis of short RNAs reveals widespread promoter-proximal stalling and arrest of Pol II in Drosophila. Science 327, 335–338 (2010).
Article ADS CAS PubMed Google Scholar
Cheng, B. & Price, D. H. Properties of RNA polymerase II elongation complexes before and after the P-TEFb-mediated transition into productive elongation. J. Biol. Chem. 282, 21901–21912 (2007).
Article CAS PubMed Google Scholar
Core, L. & Adelman, K. Promoter-proximal pausing of RNA polymerase II: a nexus of gene regulation. Genes Dev. https://doi.org/10.1101/gad.325142.119 (2019).
Article CAS Google Scholar
Henriques, T. et al. Widespread transcriptional pausing and elongation control at enhancers. Genes Dev. 32, 26–41 (2018).
Article MathSciNet CAS PubMed PubMed Central Google Scholar
Core, L. J. et al. Analysis of nascent RNA identifies a unified architecture of initiation regions at mammalian promoters and enhancers. Nat. Genet. 46, 1311–1320 (2014).
Article CAS PubMed PubMed Central Google Scholar
Krebs, A. R. et al. Genome-wide single-molecule footprinting reveals high RNA polymerase II turnover at paused promoters. Mol. Cell 67, 411–422.e4 (2017).
Article CAS PubMed PubMed Central Google Scholar
Steurer, B. et al. Live-cell analysis of endogenous GFP-RPB1 uncovers rapid turnover of initiating and promoter-paused RNA polymerase II. Proc. Natl Acad. Sci. USA 115, E4368–E4376 (2018).
Article CAS Google Scholar
Jackson, S. P. & Bartek, J. The DNA-damage response in human biology and disease. Nature 461, 1071–1078 (2009).
Article ADS CAS PubMed PubMed Central Google Scholar
Liakos, A., Lavigne, M. D. & Fousteri, M. Nucleotide excision repair: from neurodegeneration to cancer. in Personalised Medicine: Lessons from Neurodegeneration to Cancer (ed El-Khamisy, S.) 17–39 (Springer International Publishing, 2017) https://doi.org/10.1007/978-3-319-60733-7_2.
Google Scholar
Helleday, T., Eshtad, S. & Nik-Zainal, S. Mechanisms underlying mutational signatures in human cancers. Nat. Rev. Genet. 15, 585–598 (2014).
Article CAS PubMed PubMed Central Google Scholar
Vermeulen, W. & Fousteri, M. Mammalian transcription-coupled excision repair. Cold Spring Harb. Perspect. Biol. 5, a012625–a012625 (2013).
Article PubMed PubMed Central CAS Google Scholar
Hu, J., Adar, S., Selby, C. P., Lieb, J. D. & Sancar, A. Genome-wide analysis of human global and transcription-coupled excision repair of UV damage at single-nucleotide resolution. Genes Dev. 29, 948–960 (2015).
Article CAS PubMed PubMed Central Google Scholar
Lans, H., Hoeijmakers, J. H. J., Vermeulen, W. & Marteijn, J. A. The DNA damage response to transcription stress. Nat. Rev. Mol. Cell Biol. https://doi.org/10.1038/s41580-019-0169-4 (2019).
Article CAS PubMed Google Scholar
Scharer, O. D. Nucleotide excision repair in eukaryotes. Cold Spring Harb. Perspect. Biol. 5, a012609–a012609 (2013).
Article PubMed PubMed Central CAS Google Scholar
Cleaver, J. E., Lam, E. T. & Revet, I. Disorders of nucleotide excision repair: the genetic and molecular basis of heterogeneity. Nat. Rev. Genet. 10, 756–768 (2009).
Article CAS PubMed Google Scholar
Djebali, S. et al. Landscape of transcription in human cells. Nature 489, 101–108 (2012).
Article ADS CAS PubMed PubMed Central Google Scholar
Lavigne, M. D., Konstantopoulos, D., Ntakou-Zamplara, K. Z., Liakos, A. & Fousteri, M. Global unleashing of transcription elongation waves in response to genotoxic stress restricts somatic mutation rate. Nat. Commun. 8, 2076 (2017).
Article ADS PubMed PubMed Central CAS Google Scholar
Williamson, L. et al. UV irradiation induces a non-coding RNA that functionally opposes the protein encoded by the same gene. Cell 168, 843–855.e13 (2017).
Article CAS PubMed PubMed Central Google Scholar
Magnuson, B. et al. Identifying transcription start sites and active enhancer elements using BruUV-seq. Sci. Rep. 5, 17978 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Donahue, B. A., Yin, S., Taylor, J. S., Reines, D. & Hanawalt, P. C. Transcript cleavage by RNA polymerase II arrested by a cyclobutane pyrimidine dimer in the DNA template. Proc. Natl Acad. Sci. USA 91, 8502–8506 (1994).
Article ADS CAS PubMed PubMed Central Google Scholar
Rockx, D. A. et al. UV-induced inhibition of transcription involves repression of transcription initiation and phosphorylation of RNA polymerase II. Proc. Natl Acad. Sci. USA 97, 10503–10508 (2000).
Article ADS CAS PubMed PubMed Central Google Scholar
Heine, G. F., Horwitz, A. A. & Parvin, J. D. Multiple mechanisms contribute to inhibit transcription in response to DNA damage. J. Biol. Chem. 283, 9555–9561 (2008).
Article CAS PubMed PubMed Central Google Scholar
Gregersen, L. H. & Svejstrup, J. Q. The cellular response to transcription-blocking DNA damage. Trends Biochem. Sci. 43, 327–341 (2018).
Article CAS PubMed PubMed Central Google Scholar
Epanchintsev, A. et al. Cockayne’s syndrome A and B proteins regulate transcription arrest after genotoxic stress by promoting ATF3 degradation. Mol. Cell 68, 1054–1066.e6 (2017).
Article CAS PubMed Google Scholar
Gyenis, Á. et al. UVB induces a genome-wide acting negative regulatory mechanism that operates at the level of transcription initiation in human cells. PLoS Genet. 10, e1004483 (2014).
Article PubMed PubMed Central CAS Google Scholar
Andrade-Lima, L. C., Veloso, A., Paulsen, M. T., Menck, C. F. M. & Ljungman, M. DNA repair and recovery of RNA synthesis following exposure to ultraviolet light are delayed in long genes. Nucleic Acids Res. 43, 2744–2756 (2015).
Article CAS PubMed PubMed Central Google Scholar
Bugai, A. et al. P-TEFb activation by RBM7 shapes a pro-survival transcriptional response to genotoxic stress. Mol. Cell 74, 254–267.e10 (2019).
Article CAS PubMed PubMed Central Google Scholar
Chen, R. et al. PP2B and PP1 cooperatively disrupt 7SK snRNP to release P-TEFb for transcription in response to Ca2+ signaling. Genes Dev. 22, 1356–1368 (2008).
Article CAS PubMed PubMed Central Google Scholar
Borisova, M. E. et al. p38-MK2 signaling axis regulates RNA metabolism after UV-light-induced DNA damage. Nat. Commun. 9, 1017 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Chiou, Y.-Y., Hu, J., Sancar, A. & Selby, C. P. RNA polymerase II is released from the DNA template during transcription-coupled repair in mammalian cells. J. Biol. Chem. https://doi.org/10.1074/jbc.RA117.000971 (2017).
Article PubMed PubMed Central Google Scholar
Corces, M. R. et al. An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues. Nat. Methods 14, 959 (2017).
Article CAS PubMed PubMed Central Google Scholar
Polo, S. E. & Almouzni, G. Chromatin dynamics after DNA damage: the legacy of the access-repair-restore model. DNA Repair 36, 114–121 (2015).
Article CAS PubMed PubMed Central Google Scholar
Misteli, T. & Soutoglou, E. The emerging role of nuclear architecture in DNA repair and genome maintenance. Nat. Rev. Mol. Cell Biol. 10, 243–254 (2009).
Article CAS PubMed PubMed Central Google Scholar
Creyghton, M. P. et al. Histone H3K27ac separates active from poised enhancers and predicts developmental state. PNAS 107, 21931–21936 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Di Croce, L. & Helin, K. Transcriptional regulation by Polycomb group proteins. Nat. Struct. Mol. Biol. 20, 1147–1155 (2013).
Article PubMed CAS Google Scholar
Chen, Y. et al. Principles for RNA metabolism and alternative transcription initiation within closely spaced promoters. Nat. Genet. 48, 984–994 (2016).
Article CAS PubMed PubMed Central Google Scholar
Mayne, L. V. & Lehmann, A. R. Failure of RNA synthesis to recover after UV irradiation: an early defect in cells from individuals with Cockayne’s syndrome and xeroderma pigmentosum. Cancer Res. 42, 1473–1478 (1982).
CAS PubMed Google Scholar
Adar, S., Hu, J., Lieb, J. D. & Sancar, A. Genome-wide kinetics of DNA excision repair in relation to chromatin state and mutagenesis. Proc. Natl Acad. Sci. USA 113, E2124–E2133 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Mueller, B. et al. Widespread changes in nucleosome accessibility without changes in nucleosome occupancy during a rapid transcriptional induction. Genes Dev. 31, 451–462 (2017).
Article CAS PubMed PubMed Central Google Scholar
Karlić, R., Chung, H.-R., Lasserre, J., Vlahoviček, K. & Vingron, M. Histone modification levels are predictive for gene expression. PNAS 107, 2926–2931 (2010).
Article ADS PubMed PubMed Central Google Scholar
Ucar, D. et al. The chromatin accessibility signature of human immune aging stems from CD8+ T cells. J. Exp. Med. https://doi.org/10.1084/jem.20170416 (2017).
Article CAS PubMed PubMed Central Google Scholar
Beltran, M. et al. The interaction of PRC2 with RNA or chromatin is mutually antagonistic. Genome Res. 26, 896–907 (2016).
Article CAS PubMed PubMed Central Google Scholar
Ramanathan, B. & Smerdon, M. J. Changes in nuclear protein acetylation in u.v.-damaged human cells. Carcinogenesis 7, 1087–1094 (1986).
Article CAS PubMed Google Scholar
Schick, S. et al. Dynamics of chromatin accessibility and epigenetic state in response to UV damage. J. Cell Sci. 128, 4380–4394 (2015).
Article CAS PubMed Google Scholar
Li, G. & Ho, V. C. p53-dependent DNA repair and apoptosis respond differently to high- and low-dose ultraviolet radiation. Br. J. Dermatol. 139, 3–10 (1998).
Article CAS PubMed Google Scholar
Farrell, A. W., Halliday, G. M. & Lyons, J. G. Chromatin structure following UV-induced DNA damage—repair or death? Int. J. Mol. Sci. 12, 8063–8085 (2011).
Article CAS PubMed PubMed Central Google Scholar
Rehemtulla, A., Hamilton, C. A., Chinnaiyan, A. M. & Dixit, V. M. Ultraviolet radiation-induced apoptosis is mediated by activation of CD-95 (Fas/APO-1). J. Biol. Chem. 272, 25783–25786 (1997).
Article CAS PubMed Google Scholar
Caricchio, R., McPhie, L. & Cohen, P. L. Ultraviolet B radiation-induced cell death: critical role of ultraviolet dose in inflammation and lupus autoantigen redistribution. J. Immunol. 171, 5778–5786 (2003).
Article CAS PubMed Google Scholar
Ibrahim, M. M. et al. Determinants of promoter and enhancer transcription directionality in metazoans. Nat. Commun. 9, 4472 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Scruggs, B. S. et al. Bidirectional transcription arises from two distinct hubs of transcription factor binding and active chromatin. Mol. Cell 58, 1101–1112 (2015).
Article CAS PubMed PubMed Central Google Scholar
Blasius, M., Wagner, S. A., Choudhary, C., Bartek, J. & Jackson, S. P. A quantitative 14-3-3 interaction screen connects the nuclear exosome targeting complex to the DNA damage response. Genes Dev. 28, 1977–1982 (2014).
Article CAS PubMed PubMed Central Google Scholar
Tiedje, C. et al. p38 ^MAPK /MK2-mediated phosphorylation of RBM7 regulates the human nuclear exosome targeting complex. RNA 21, 262–278 (2015).
Article PubMed PubMed Central Google Scholar
Boeing, S. et al. Multiomic analysis of the UV-induced DNA damage response. Cell Rep. 15, 1597–1610 (2016).
Article CAS PubMed Google Scholar
Ehrensberger, A. H., Kelly, G. P. & Svejstrup, J. Q. Mechanistic interpretation of promoter-proximal peaks and RNAPII density maps. Cell 154, 713–715 (2013).
Article CAS PubMed Google Scholar
Adebali, O., Chiou, Y.-Y., Hu, J., Sancar, A. & Selby, C. P. Genome-wide transcription-coupled repair in Escherichia coli is mediated by the Mfd translocase. PNAS 114, E2116–E2125 (2017).
Article CAS PubMed PubMed Central Google Scholar
Murray, S. C. & Mellor, J. Using both strands: the fundamental nature of antisense transcription. Bioarchitecture 6, 12–21 (2016).
Article PubMed PubMed Central Google Scholar
Seila, A. C., Core, L. J., Lis, J. T. & Sharp, P. A. Divergent transcription: a new feature of active promoters. Cell Cycle 8, 2557–2564 (2009).
Article CAS PubMed Google Scholar
Haradhvala, N. J. et al. Mutational strand asymmetries in cancer genomes reveal mechanisms of DNA damage and repair. Cell 164, 538–549 (2016).
Article CAS PubMed PubMed Central Google Scholar
Sabarinathan, R., Mularoni, L., Deu-Pons, J., Gonzalez-Perez, A. & López-Bigas, N. Nucleotide excision repair is impaired by binding of transcription factors to DNA. Nature 532, 264–267 (2016).
Article ADS CAS PubMed Google Scholar
Shao, W. & Zeitlinger, J. Paused RNA polymerase II inhibits new transcriptional initiation. Nat. Genet. 49, 1045–1051 (2017).
Article CAS PubMed Google Scholar
Fitz, J., Neumann, T. & Pavri, R. Regulation of RNA polymerase II processivity by Spt5 is restricted to a narrow window during elongation. EMBO J. 37, e97965 (2018).
Article PubMed PubMed Central CAS Google Scholar
Gressel, S. et al. CDK9-dependent RNA polymerase II pausing controls transcription initiation. eLife 6, e29736 (2017).
Article PubMed PubMed Central Google Scholar
Xu, S., Grullon, S., Ge, K. & Peng, W. Spatial clustering for identification of ChIP-enriched regions (SICER) to map regions of histone methylation patterns in embryonic stem cells. Methods Mol. Biol. 1150, 97–111 (2014).
Article CAS PubMed PubMed Central Google Scholar
Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).
Article PubMed PubMed Central CAS Google Scholar

Download references

Acknowledgements

We thank Pantelis Hatzis, Mihalis Verykokakis, and members of the Fousteri lab for critical discussions and reading of the paper. We thank Panagiotis Moulos for help and discussion on the underdevelopment computational pipelines. We thank Vladimir Benes and the Genecore facility (EMBL, Germany) for the special care they use in sequencing our NGS libraries. This work was funded by a European Research Council grant to M.F., Agreement-309612 (TransArrest) and <Matching Funds> to MF funded by National sources.

Author information

These authors contributed equally: Anastasios Liakos, Dimitris Konstantopoulos.

Authors and Affiliations

Institute for Fundamental Biomedical Research, BSRC ‘Alexander Fleming’, 34 Fleming st., 16672, Vari, Athens, Greece
Anastasios Liakos, Dimitris Konstantopoulos, Matthieu D. Lavigne & Maria Fousteri
Department of Biology, School of Science, National & Kapodistrian University of Athens, Athens, Greece
Anastasios Liakos
Department of Biology, University of Crete, 70013, Herakleion, Greece
Dimitris Konstantopoulos

Authors

Anastasios Liakos
View author publications
You can also search for this author in PubMed Google Scholar
Dimitris Konstantopoulos
View author publications
You can also search for this author in PubMed Google Scholar
Matthieu D. Lavigne
View author publications
You can also search for this author in PubMed Google Scholar
Maria Fousteri
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

M.F. and M.D.L. designed the study and were responsible for interpretation of the results. M.F. directed the study and obtained financial support. M.D.L., M.F., and A.L. wrote the paper and all authors edited the paper. A.L. performed the experimental part of the study. D.K. performed the statistical and bioinformatics analyses. M.D.L. contributed significantly in the bioinformatics analysis of the data. All authors discussed the results, reviewed, commented on, and approved the final version of the paper. A.L and D.K. contributed equally to the paper.

Corresponding authors

Correspondence to Matthieu D. Lavigne or Maria Fousteri.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Reporting Summary

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Liakos, A., Konstantopoulos, D., Lavigne, M.D. et al. Continuous transcription initiation guarantees robust repair of all transcribed genes and regulatory regions. Nat Commun 11, 916 (2020). https://doi.org/10.1038/s41467-020-14566-9

Download citation

Received: 11 July 2019
Accepted: 16 January 2020
Published: 14 February 2020
DOI: https://doi.org/10.1038/s41467-020-14566-9

This article is cited by

Temporal-iCLIP captures co-transcriptional RNA-protein interactions
- Ross A. Cordiner
- Yuhui Dou
- Torben Heick Jensen
Nature Communications (2023)
Cellular fractionation reveals transcriptome responses of human fibroblasts to UV-C irradiation
- Jiena Liu
- Zhenzhen Wu
- Yuming Wang
Cell Death & Disease (2022)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.