RNA polymerase II speed: a key player in controlling and adapting transcriptome composition.

RNA polymerase II (RNA Pol II) speed or elongation rate, i.e., the number of nucleotides synthesized per unit of time, is a major determinant of transcriptome composition. It controls co-transcriptional processes such as splicing, polyadenylation, and transcription termination, thus regulating the production of alternative splice variants, circular RNAs, alternatively polyadenylated transcripts, or read-through transcripts. RNA Pol II speed itself is regulated in response to intra- and extra-cellular stimuli and can in turn affect the transcriptome composition in response to these stimuli. Evidence points to a potentially important role of transcriptome composition modification through RNA Pol II speed regulation for adaptation of cells to a changing environment, thus pointing to a function of RNA Pol II speed regulation in cellular physiology. Analyzing RNA Pol II speed dynamics may therefore be central to fully understand the regulation of physiological processes, such as the development of multicellular organisms. Recent findings also raise the possibility that RNA Pol II speed deregulation can be detrimental and participate in disease progression. Here, we review initial and current approaches to measure RNA Pol II speed, as well as providing an overview of the factors controlling speed and the co-transcriptional processes which are affected. Finally, we discuss the role of RNA Pol II speed regulation in cell physiology.

Very long intergenic non-coding (vlinc) RNAs directly regulate multiple genes in cis and trans.

BACKGROUND: The majority of the human genome is transcribed in the form of long non-coding (lnc) RNAs. While these transcripts have attracted considerable interest, their molecular mechanisms of function and biological significance remain controversial. One of the main reasons behind this lies in the significant challenges posed by lncRNAs requiring the development of novel methods and concepts to unravel their functionality. Existing methods often lack cross-validation and independent confirmation by different methodologies and therefore leave significant ambiguity as to the authenticity of the outcomes. Nonetheless, despite all the caveats, it appears that lncRNAs may function, at least in part, by regulating other genes via chromatin interactions. Therefore, the function of a lncRNA could be inferred from the function of genes it regulates. In this work, we present a genome-wide functional annotation strategy for lncRNAs based on identification of their regulatory networks via the integration of three distinct types of approaches: co-expression analysis, mapping of lncRNA-chromatin interactions, and assaying molecular effects of lncRNA knockdowns obtained using an inducible and highly specific CRISPR/Cas13 system. RESULTS: We applied the strategy to annotate 407 very long intergenic non-coding (vlinc) RNAs belonging to a novel widespread subclass of lncRNAs. We show that vlincRNAs indeed appear to regulate multiple genes encoding proteins predominantly involved in RNA- and development-related functions, cell cycle, and cellular adhesion via a mechanism involving proximity between vlincRNAs and their targets in the nucleus. A typical vlincRNAs can be both a positive and negative regulator and regulate multiple genes both in trans and cis. Finally, we show vlincRNAs and their regulatory networks potentially represent novel components of DNA damage response and are functionally important for the ability of cancer cells to survive genotoxic stress. CONCLUSIONS: This study provides strong evidence for the regulatory role of the vlincRNA class of lncRNAs and a potentially important role played by these transcripts in the hidden layer of RNA-based regulation in complex biological systems.

Circular ANRIL isoforms switch from repressors to activators of p15/CDKN2B expression during RAF1 oncogene-induced senescence.

Long non-coding RNAs (ncRNAs) are major regulators of gene expression and cell fate. The INK4 locus encodes the tumour suppressor proteins p15INK4b, p16INK4a and p14ARF required for cell cycle arrest and whose expression increases during senescence. ANRIL is a ncRNA antisense to the p15 gene. In proliferative cells, ANRIL prevents senescence by repressing INK4 genes through the recruitment of Polycomb-group proteins. In models of replicative and RASval12 oncogene-induced senescence (OIS), the expression of ANRIL and Polycomb proteins decreases, thus allowing INK4 derepression. Here, we found in a model of RAF1 OIS that ANRIL expression rather increases, due in particular to an increased stability. This led us to search for circular ANRIL isoforms, as circular RNAs are rather stable species. We found that the expression of two circular ANRIL increases in several OIS models (RAF1, MEK1 and BRAF). In proliferative cells, they repress p15 expression, while in RAF1 OIS, they promote full induction of p15, p16 and p14ARF expression. Further analysis of one of these circular ANRIL shows that it interacts with Polycomb proteins and decreases EZH2 Polycomb protein localization and H3K27me3 at the p15 and p16 promoters, respectively. We propose that changes in the ratio between Polycomb proteins and circular ANRIL isoforms allow these isoforms to switch from repressors of p15 gene to activators of all INK4 genes in RAF1 OIS. Our data reveal that regulation of ANRIL expression depends on the senescence inducer and underline the importance of circular ANRIL in the regulation of INK4 gene expression and senescence.

Functional annotation of the vlinc class of non-coding RNAs using systems biology approach.

Functionality of the non-coding transcripts encoded by the human genome is the coveted goal of the modern genomics research. While commonly relied on the classical methods of forward genetics, integration of different genomics datasets in a global Systems Biology fashion presents a more productive avenue of achieving this very complex aim. Here we report application of a Systems Biology-based approach to dissect functionality of a newly identified vast class of very long intergenic non-coding (vlinc) RNAs. Using highly quantitative FANTOM5 CAGE dataset, we show that these RNAs could be grouped into 1542 novel human genes based on analysis of insulators that we show here indeed function as genomic barrier elements. We show that vlinc RNAs genes likely function in cisto activate nearby genes. This effect while most pronounced in closely spaced vlinc RNA-gene pairs can be detected over relatively large genomic distances. Furthermore, we identified 101 vlinc RNA genes likely involved in early embryogenesis based on patterns of their expression and regulation. We also found another 109 such genes potentially involved in cellular functions also happening at early stages of development such as proliferation, migration and apoptosis. Overall, we show that Systems Biology-based methods have great promise for functional annotation of non-coding RNAs.

High-resolution profiling of gammaH2AX around DNA double strand breaks in the mammalian genome.

Chromatin acts as a key regulator of DNA-related processes such as DNA damage repair. Although ChIP-chip is a powerful technique to provide high-resolution maps of protein-genome interactions, its use to study DNA double strand break (DSB) repair has been hindered by the limitations of the available damage induction methods. We have developed a human cell line that permits induction of multiple DSBs randomly distributed and unambiguously positioned within the genome. Using this system, we have generated the first genome-wide mapping of gammaH2AX around DSBs. We found that all DSBs trigger large gammaH2AX domains, which spread out from the DSB in a bidirectional, discontinuous and not necessarily symmetrical manner. The distribution of gammaH2AX within domains is influenced by gene transcription, as parallel mappings of RNA Polymerase II and strand-specific expression showed that gammaH2AX does not propagate on active genes. In addition, we showed that transcription is accurately maintained within gammaH2AX domains, indicating that mechanisms may exist to protect gene transcription from gammaH2AX spreading and from the chromatin rearrangements induced by DSBs.

Cell cycle control of centromeric repeat transcription and heterochromatin assembly.

Heterochromatin in eukaryotic genomes regulates diverse chromosomal processes including transcriptional silencing. However, in Schizosaccharomyces pombe RNA polymerase II (RNAPII) transcription of centromeric repeats is essential for RNA-interference-mediated heterochromatin assembly. Here we study heterochromatin dynamics during the cell cycle and its effect on RNAPII transcription. We describe a brief period during the S phase of the cell cycle in which RNAPII preferentially transcribes centromeric repeats. This period is enforced by heterochromatin, which restricts RNAPII accessibility at centromeric repeats for most of the cell cycle. RNAPII transcription during S phase is linked to loading of RNA interference and heterochromatin factors such as the Ago1 subunit of the RITS complex and the Clr4 methyltransferase complex subunit Rik1 (ref. 7). Moreover, Set2, an RNAPII-associated methyltransferase that methylates histone H3 lysine 36 at repeat loci during S phase, acts in a pathway parallel to Clr4 to promote heterochromatin assembly. We also show that phosphorylation of histone H3 serine 10 alters heterochromatin during mitosis, correlating with recruitment of condensin that affects silencing of centromeric repeats. Our analyses suggest at least two distinct modes of heterochromatin targeting to centromeric repeats, whereby RNAPII transcription of repeats and chromodomain proteins bound to methylated histone H3 lysine 9 mediate recruitment of silencing factors. Together, these processes probably facilitate heterochromatin maintenance through successive cell divisions.

TDP1 mutation causing SCAN1 neurodegenerative syndrome hampers the repair of transcriptional DNA double-strand breaks.

TDP1 removes transcription-blocking topoisomerase I cleavage complexes (TOP1ccs), and its inactivating H493R mutation causes the neurodegenerative syndrome SCAN1. However, the molecular mechanism underlying the SCAN1 phenotype is unclear. Here, we generate human SCAN1 cell models using CRISPR-Cas9 and show that they accumulate TOP1ccs along with changes in gene expression and genomic distribution of R-loops. SCAN1 cells also accumulate transcriptional DNA double-strand breaks (DSBs) specifically in the G1 cell population due to increased DSB formation and lack of repair, both resulting from abortive removal of transcription-blocking TOP1ccs. Deficient TDP1 activity causes increased DSB production, and the presence of mutated TDP1 protein hampers DSB repair by a TDP2-dependent backup pathway. This study provides powerful models to study TDP1 functions under physiological and pathological conditions and unravels that a gain of function of the mutated TDP1 protein, which prevents DSB repair, rather than a loss of TDP1 activity itself, could contribute to SCAN1 pathogenesis.

Distinct roles of HDAC complexes in promoter silencing, antisense suppression and DNA damage protection.

Histone acetylation is important in regulating DNA accessibility. Multifunctional Sin3 proteins bind histone deacetylases (HDACs) to assemble silencing complexes that selectively target chromatin. We show that, in fission yeast, an essential HDAC, Clr6, exists in two distinct Sin3 core complexes. Complex I contains an essential Sin3 homolog, Pst1, and other factors, and predominantly targets gene promoters. Complex II contains a nonessential Sin3 homolog, Pst2, and several conserved proteins. It preferentially targets transcribed chromosomal regions and centromere cores. Defects in complex II abrogate global protective functions of chromatin, causing increased accessibility of DNA to genotoxic agents and widespread antisense transcripts that are processed by the exosome. Notably, the two Clr6 complexes differentially repress forward and reverse centromeric repeat transcripts, suggesting that these complexes regulate transcription in heterochromatin and euchromatin in similar manners, including suppression of spurious transcripts from cryptic start sites.

Small Interfering RNAs Targeting a Chromatin-Associated RNA Induce Its Transcriptional Silencing in Human Cells.

Transcriptional gene silencing by small interfering RNAs (siRNAs) has been widely described in various species, including plants and yeast. In mammals, its extent remains somewhat debated. Previous studies showed that siRNAs targeting gene promoters could induce the silencing of the targeted promoter, although the involvement of off-target mechanisms was also suggested. Here, by using nascent RNA capture and RNA polymerase II chromatin immunoprecipitation, we show that siRNAs targeting a chromatin-associated noncoding RNA induced its transcriptional silencing. Deletion of the sequence targeted by one of these siRNAs on the two alleles by genome editing further showed that this silencing was due to base-pairing of the siRNA to the target. Moreover, by using cells with heterozygous deletion of the target sequence, we showed that only the wild-type allele, but not the deleted allele, was silenced by the siRNA, indicating that transcriptional silencing occurred only in cis. Finally, we demonstrated that both Ago1 and Ago2 are involved in this transcriptional silencing. Altogether, our data demonstrate that siRNAs targeting a chromatin-associated RNA at a distance from its promoter induce its transcriptional silencing. Our results thus extend the possible repertoire of endogenous or exogenous interfering RNAs.

Impact of COVID-19 on the association between pulse oximetry and arterial oxygenation in patients with acute respiratory distress syndrome.

Managing patients with acute respiratory distress syndrome (ARDS) requires frequent changes in mechanical ventilator respiratory settings to optimize arterial oxygenation assessed by arterial oxygen partial pressure (PaO2) and saturation (SaO2). Pulse oxymetry (SpO2) has been suggested as a non-invasive surrogate for arterial oxygenation however its accuracy in COVID-19 patients is unknown. In this study, we aimed to investigate the influence of COVID-19 status on the association between SpO2 and arterial oxygenation. We prospectively included patients with ARDS and compared COVID-19 to non-COVID-19 patients, regarding SpO2 and concomitant arterial oxygenation (SaO2 and PaO2) measurements, and their association. Bias was defined as mean difference between SpO2 and SaO2 measurements. Occult hypoxemia was defined as a SpO2 ≥ 92% while concomitant SaO2 < 88%. Multiple linear regression models were built to account for confounders. We also assessed concordance between positive end-expiratory pressure (PEEP) trial-induced changes in SpO2 and in arterial oxygenation. We included 55 patients, among them 26 (47%) with COVID-19. Overall, SpO2 and SaO2 measurements were correlated (r = 0.70; p < 0.0001), however less so in COVID-19 than in non-COVID-19 patients (r = 0.55, p < 0.0001 vs. r = 0.84, p < 0.0001, p = 0.002 for intergroup comparison). Bias was + 1.1%, greater in COVID-19 than in non-COVID-19 patients (2.0 vs. 0.3%; p = 0.02). In multivariate analysis, bias was associated with COVID-19 status (unstandardized ß = 1.77, 95%CI = 0.38-3.15, p = 0.01), ethnic group and ARDS severity. Occult hypoxemia occurred in 5.5% of measurements (7.7% in COVID-19 patients vs. 3.4% in non-COVID-19 patients, p = 0.42). Concordance rate between PEEP trial-induced changes in SpO2 and SaO2 was 84%, however less so in COVID-19 than in non-COVID-19 patients (69% vs. 97%, respectively). Similar results were observed for PaO2 regarding correlations, bias, and concordance with SpO2 changes. In patients with ARDS, SpO2 was associated with arterial oxygenation, but COVID-19 status significantly altered this association.

Ataxia telangiectasia mutated activation by transcription- and topoisomerase I-induced DNA double-strand breaks.

Ataxia telangiectasia mutated (ATM), the deficiency of which causes a severe neurodegenerative disease, is a crucial mediator for the DNA damage response (DDR). As neurons have high rates of transcription that require topoisomerase I (TOP1), we investigated whether TOP1 cleavage complexes (TOP1cc)-which are potent transcription-blocking lesions-also produce transcription-dependent DNA double-strand breaks (DSBs) with ATM activation. We show the induction of DSBs and DDR activation in post-mitotic primary neurons and lymphocytes treated with camptothecin, with the induction of nuclear DDR foci containing activated ATM, gamma-H2AX (phosphorylated histone H2AX), activated CHK2 (checkpoint kinase 2), MDC1 (mediator of DNA damage checkpoint 1) and 53BP1 (p53 binding protein 1). The DSB-ATM-DDR pathway was suppressed by inhibiting transcription and gamma-H2AX signals were reduced by RNase H1 transfection, which removes transcription-mediated R-loops. Thus, we propose that Top1cc produce transcription arrests with R-loop formation and generate DSBs that activate ATM in post-mitotic cells.

Integrated analysis of H2A.Z isoforms function reveals a complex interplay in gene regulation.

The H2A.Z histone variant plays major roles in the control of gene expression. In human, H2A.Z is encoded by two genes expressing two isoforms, H2A.Z.1 and H2A.Z.2 differing by three amino acids. Here, we undertook an integrated analysis of their functions in gene expression using endogenously-tagged proteins. RNA-Seq analysis in untransformed cells showed that they can regulate both distinct and overlapping sets of genes positively or negatively in a context-dependent manner. Furthermore, they have similar or antagonistic function depending on genes. H2A.Z.1 and H2A.Z.2 can replace each other at Transcription Start Sites, providing a molecular explanation for this interplay. Mass spectrometry analysis showed that H2A.Z.1 and H2A.Z.2 have specific interactors, which can mediate their functional antagonism. Our data indicate that the balance between H2A.Z.1 and H2A.Z.2 at promoters is critically important to regulate specific gene expression, providing an additional layer of complexity to the control of gene expression by histone variants.

Balance between acetylation and methylation of histone H3 lysine 9 on the E2F-responsive dihydrofolate reductase promoter.

Epigenetic marks that specify silent heterochromatic domains in eucaryotic genomes include methylation of histone H3 lysine 9. Strikingly, active loci in the vicinity of silent domains are sometimes characterized by acetylation of histone H3 lysine 9, suggesting that the balance between these two competitive modifications is important for the establishment of specific chromatin structures. Some euchromatic genes, targeted by the retinoblastoma protein Rb, are also believed to be regulated by histone H3 lysine 9 methylation. Here, we study the dihydrofolate reductase promoter, which is repressed in G0 and at the beginning of G1 by p107 or p130, two Rb-related proteins. We found that these two pocket proteins share with Rb the ability to associate with the histone methyl transferase SUV39H1. SUV39H1 can be recruited to the E2F transcription factor and functions as a transcriptional corepressor. With ChIP assays followed by real-time PCR, we showed that K9 of histone H3 evolves from a hypermethylated state in G0 to a hyperacetylated state at the G1/S transition. Taken together, these results indicate that the temporal regulation of euchromatic promoters may involve controlling the balance between methylation and acetylation of histone H3 lysine 9, a feature previously described for the spatial regulation of chromatin function.

Control of Gene Expression in Senescence through Transcriptional Read-Through of Convergent Protein-Coding Genes.

Antisense RNAs are non-coding RNAs that can regulate their corresponding sense RNAs and are generally produced from specific promoters. We uncover here a family of antisense RNAs, named START RNAs, produced during cellular senescence by transcriptional read-through at convergent protein-coding genes. Importantly, START RNAs repress the expression of their corresponding sense RNAs. In proliferative cells, we found that the Pol II elongation rate is limited downstream of TTS at START loci, allowing transcription termination to occur before Pol II reaches the convergent genes, thus preventing antisense RNA production and interference with the expression of the convergent genes. START RNAs are repressed by H2A.Z histone variant, whose local occupancy decreases in senescence. Our results thus uncover a mechanism of gene expression regulation relying on read-through antisense transcript production at convergent genes, underlining the functional importance of chromatin regulation in the control of RNA pol II elongation rate at intergenic regions.

Functional and physical interaction between the histone methyl transferase Suv39H1 and histone deacetylases.

The histone methyl transferase Suv39H1 is involved in silencing by pericentric heterochromatin. It specifically methylates K9 of histone H3, thereby creating a high affinity binding site for HP1 proteins. We and others have shown recently that it is also involved in transcriptional repression by the retinoblastoma protein Rb. Strikingly, both HP1 localisation and repression by Rb also require, at least in part, histone deacetylases. We found here that repression of a heterologous promoter by Suv39H1 is dependent on histone deacetylase activity. However, the enzymatic activity of Suv39H1 is not required, since the N-terminal part is by itself a transcriptional repression domain. Coimmunoprecipitation experiments indicated that Suv39H1 can physically interact with HDAC1, -2 and -3, therefore suggesting that transcriptional repression by Suv39H1 could be the consequence of histone deacetylases recruitment. Consistent with this interpretation, the N-terminal transcriptional repression domain of Suv39H1 bound the so-called 'core histone deacetylase complex', composed of HDAC1, HDAC2 and the Rb-associated proteins RbAp48 and RbAp46. Taken together, our results suggest that a complex containing both the Suv39H1 histone methyl transferase and histone deacetylases could be involved in heterochromatin silencing or transcriptional repression by Rb.

A vlincRNA participates in senescence maintenance by relieving H2AZ-mediated repression at the INK4 locus.

Non-coding RNAs (ncRNAs) play major roles in proper chromatin organization and function. Senescence, a strong anti-proliferative process and a major anticancer barrier, is associated with dramatic chromatin reorganization in heterochromatin foci. Here we analyze strand-specific transcriptome changes during oncogene-induced human senescence. Strikingly, while differentially expressed RNAs are mostly repressed during senescence, ncRNAs belonging to the recently described vlincRNA (very long intergenic ncRNA) class are mainly activated. We show that VAD, a novel antisense vlincRNA strongly induced during senescence, is required for the maintenance of senescence features. VAD modulates chromatin structure in cis and activates gene expression in trans at the INK4 locus, which encodes cell cycle inhibitors important for senescence-associated cell proliferation arrest. Importantly, VAD inhibits the incorporation of the repressive histone variant H2A.Z at INK4 gene promoters in senescent cells. Our data underline the importance of vlincRNAs as sensors of cellular environment changes and as mediators of the correct transcriptional response.

VlincRNAs controlled by retroviral elements are a hallmark of pluripotency and cancer.

BACKGROUND: The function of the non-coding portion of the human genome remains one of the most important questions of our time. Its vast complexity is exemplified by the recent identification of an unusual and notable component of the transcriptome - very long intergenic non-coding RNAs, termed vlincRNAs. RESULTS: Here we identify 2,147 vlincRNAs covering 10 percent of our genome. We show they are present not only in cancerous cells, but also in primary cells and normal human tissues, and are controlled by canonical promoters. Furthermore, vlincRNA promoters frequently originate from within endogenous retroviral sequences. Strikingly, the number of vlincRNAs expressed from endogenous retroviral promoters strongly correlates with pluripotency or the degree of malignant transformation. These results suggest a previously unknown connection between the pluripotent state and cancer via retroviral repeat-driven expression of vlincRNAs. Finally, we show that vlincRNAs can be syntenically conserved in humans and mouse and their depletion using RNAi can cause apoptosis in cancerous cells. CONCLUSIONS: These intriguing observations suggest that vlincRNAs could create a framework that combines many existing short ESTs and lincRNAs into a landscape of very long transcripts functioning in the regulation of gene expression in the nucleus. Certain types of vlincRNAs participate at specific stages of normal development and, based on analysis of a limited set of cancerous and primary cell lines, they appear to be co-opted by cancer-associated transcriptional programs. This provides additional understanding of transcriptome regulation during the malignant state, and could lead to additional targets and options for its reversal.

Hyperphosphorylation of RNA polymerase II in response to topoisomerase I cleavage complexes and its association with transcription- and BRCA1-dependent degradation of topoisomerase I.

The progression of RNA polymerase II can be blocked by lesions on the DNA template. In this study, we focused on the modifications of the largest subunit of RNA polymerase II, Rpb1, in response to stabilized topoisomerase I (Top1)-DNA cleavage complexes. In addition to DNA modifications (base damages and strand breaks), Top1 cleavage complexes can be trapped by camptothecin (CPT) and its derivatives used in cancer treatment. We found that, within a few minutes, CPT produces the complete hyperphosphorylation of Rpb1 in both primary and transformed cancer cells. Hyperphosphorylation is rapidly reversible following CPT removal. Hyperphosphorylation occurs selectively on the serine 5 residue of the conserved heptapeptide repeats in the Rpb1 carboxy-terminal domain and is mediated principally by the transcription factor IIH-associated cyclin-dependent kinase Cdk7. Hyperphosphorylated Rpb1 is not primarily targeted for proteosomal degradation and instead is subjected to cycles of phosphorylation and dephosphorylation as long as Top1 cleavage complexes are trapped by CPT. Finally, we show that transcription-induced degradation of Top1 is Brca1 dependent, suggesting a role for Brca1 in the repair or removal of transcription-blocking Top1-DNA cleavage complexes.

Fission yeast homologs of human histone H3 lysine 4 demethylase regulate a common set of genes with diverse functions.

Schizosaccharomyces pombe contains two proteins, SWIRM1 and SWIRM2, with close homology to human histone H3 lysine 4 demethylase. Both proteins contain the amino oxidase catalytic domain and a recently described DNA interaction SWIRM domain. Here we describe the biochemical isolation and the functional characterization of SWIRM1 and SWIRM2. Our results indicate that while SWIRM2 is an essential gene, cells lacking SWIRM1 are viable. We found that SWIRM1 and SWIRM2 are stably associated in a multiprotein complex, but intriguingly, unlike their human counterpart, S. pombe SWIRM complex contains neither a histone deacetylase nor any detectable demethylase activity. Genome-wide chromatin immunoprecipitation unexpectedly showed the absence of both SWIRM proteins from heterochromatic domains. Instead, consistent with biochemical analyses, SWIRM1 and SWIRM2 co-localize to a common set of target gene promoters whose functions are implicated in diverse processes including mitochondrial metabolism and transcriptional regulation. Importantly, we show that SWIRM1 is not only required for optimum transcription of its target genes but also display a global role in regulation of antisense transcription.