Skip to main content
Advertisement

Main menu

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research

User menu

  • My alerts

Search

  • Advanced search
Life Science Alliance
  • Other Publications
    • EMBO Press
    • The EMBO Journal
    • EMBO reports
    • EMBO Molecular Medicine
    • Molecular Systems Biology
    • Rockefeller University Press
    • Journal of Cell Biology
    • Journal of Experimental Medicine
    • Journal of General Physiology
    • Cold Spring Harbor Laboratory Press
    • Genes & Development
    • Genome Research
  • My alerts
Life Science Alliance

Advanced Search

  • Home
  • Articles
    • Newest Articles
    • Current Issue
    • Archive
    • Subjects
  • Submit
    • Submit a Manuscript
    • Author Guidelines
    • License, Copyright, Fee
    • FAQ
    • Why Submit
  • About
    • About Us
    • Editors & Staff
    • Board Members
    • Licensing and Reuse
    • Reviewer Guidelines
    • Privacy Policy
    • Advertise
    • Contact Us
    • LSA LLC
  • Alerts
  • Follow lsa Template on Twitter
Resource
Transparent Process
Open Access

Genome-wide R-loop analysis defines unique roles for DDX5, XRN2, and PRMT5 in DNA/RNA hybrid resolution

View ORCID ProfileOscar D Villarreal, Sofiane Y Mersaoui, Zhenbao Yu, Jean-Yves Masson, View ORCID ProfileStéphane Richard  Correspondence email
Oscar D Villarreal
1Segal Cancer Center, Lady Davis Institute for Medical Research and Gerald Bronfman Department of Oncology and Departments of Biochemistry, Human Genetics and Medicine, McGill University, Montréal, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Oscar D Villarreal
Sofiane Y Mersaoui
1Segal Cancer Center, Lady Davis Institute for Medical Research and Gerald Bronfman Department of Oncology and Departments of Biochemistry, Human Genetics and Medicine, McGill University, Montréal, Canada
2Genome Stability Laboratory, Centre Hospitalier Universitaire de Québec Research Center, Oncology Axis; Department of Molecular Biology, Medical Biochemistry and Pathology; Laval University Cancer Research Center, Québec City, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhenbao Yu
1Segal Cancer Center, Lady Davis Institute for Medical Research and Gerald Bronfman Department of Oncology and Departments of Biochemistry, Human Genetics and Medicine, McGill University, Montréal, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jean-Yves Masson
2Genome Stability Laboratory, Centre Hospitalier Universitaire de Québec Research Center, Oncology Axis; Department of Molecular Biology, Medical Biochemistry and Pathology; Laval University Cancer Research Center, Québec City, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stéphane Richard
1Segal Cancer Center, Lady Davis Institute for Medical Research and Gerald Bronfman Department of Oncology and Departments of Biochemistry, Human Genetics and Medicine, McGill University, Montréal, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stéphane Richard
  • For correspondence: stephane.richard@mcgill.ca
Published 3 August 2020. DOI: 10.26508/lsa.202000762
  • Article
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF
Loading

Article Figures & Data

Figures

  • Tables
  • Supplementary Materials
  • Figure 1.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 1. Genome-wide R-loops in DDX5, XRN2, and PRMT5 knockdown cells.

    (A) U2OS cells were transfected with siRNAs for control (CTL), DDX5, XRN2, and PRMT5 and cell lysates were separated by SDS–PAGE and immunoblotted with the indicated antibodies to confirm successful knockdowns. β-Actin was used as a loading control. (B) Total read counts within R-loop peaks in each knockdown condition, normalized to library size, and averaged for the two biological replicates. Peaks called by the MACS algorithm v2.2.6 in broad mode (q-value < 0.1) for each replicate were merged into a consensus list across all treatments through DiffBind v2.14.0. Error bars denote SD of the replicates. (C) Read coverage of a representative peak profile with a gain in R-loop signal relative to siCTL at the RFNG-GPS1-DUS1L gene loci for all cells, generated by Integrative Genomic Viewer v2.8.0. DNA extracted the control DNA treated with RNase-H (black); siCTL cells (cyan); siDDX5 (red); siPRMT5 (magenta); and siXRN2 (green).

  • Figure S1.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure S1.

    R-loop signal reproducibility and comparison within consensus and representative peaks. (A) Scatter plots of read counts normalized to library size within the R-loop consensus peaks showing the reproducibility and Pearson correlation coefficient between the biological replicates (termed “A” and “B”) of siCTL, siDDX5, siXRN2, and siPRMT5 U2OS cells. (B) Read coverage of representative peak profiles with R-loop signal gain, showing the FOS (left), SSTR5 (center), and SPIB-MYBPC2 (right) gene loci, generated by Integrative Genomic Viewer v2.8.0. Control cells are cyan; cells were treated with RNase H (black), siDDX5 (red), siPRMT5 (magenta), or siXRN2 (green). (C) Volcano plots of differential R-loop signal relative to siCTL in the consensus list of peaks, highlighting significant gains or losses upon knockdown treatment as reported by DESeq2 v1.26.0 (absolute log2 fold change > 1 and false discovery rate < 0.1, Wald test). (D) Scatter plot of log2 normalized mean read concentration within R-loop consensus peaks for each comparison, highlighting significant gains or losses in signal upon the treatment.

  • Figure S2.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure S2.

    DRIP-qPCR analysis of gain R-loops identified in the DRIP-seq screen. (A) Validation of the hybrid signal (R-loop) obtained in the DRIP-seq analysis using DRIP-qPCR at the indicated loci in siCTL (red) and siDDX5 (blue) U2OS. Upper panels show R-loop signals at each locus obtained by sequencing (Integrative Genomic Viewer) and lower panels indicated the R-loop signals relative to input obtained by DRIP-qPCR at the indicated locus. The bar graphs show the average and SEM from three independent experiments. Statistical significance was assessed using t test. *P < 0.05, **P < 0.01, and ****P < 0.0001. (A, B) Western blotting of protein extract of siCTL and siDDX5 transfected U2OS used for panel (A) for experiment #1 and #2. The migration of DDX5 and β-ACTIN (ACTB, loading control) is shown with an arrowhead.

  • Figure 2.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 2. R-loop gains and losses in siDDX5, siXRN2, and siPRMT5 cells.

    (A) Total amount of R-loop gain and loss consensus peaks called for each knockdown condition relative to siCTL by DESeq2 v1.26.0 (absolute log2 fold change > 1 and false discovery rate < 0.1, Wald test). (B) Venn diagrams showing the overlaps among consensus peaks with gain (top) or loss (bottom) in R-loop signal upon each knockdown condition. (C) Distribution of log2 normalized read concentration within R-loop gain, loss or unchanged consensus peaks at the control (red) and treated (siDDX5, green; siXRN2, cyan; and siPRMT5, purple) cells. Peaks are split into three panels for each knockdown treatment.

  • Figure S3.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure S3.

    RNA-seq gene expression analysis in the knockdown treatments. (A) Principal component analysis plot showing the first two principal components of the RNA-seq gene expression results for each knockdown treatment in biological triplicates, quantified through HOMER v4.11.1 and normalized through the regularized logarithm (rlog) transformation of DESeq2 v1.26.0. (B) Same as (A) but displayed as a heat map of the Pearson correlation coefficients between the normalized gene expression of the samples. (C) Volcano plot of the RNA-seq differential expression for the genes lying nearest to the DRIP-seq consensus peaks, highlighting genes that lie near peaks with a gain (red) or loss (blue) in R-loop signal upon the corresponding knockdown treatment. Dashed lines indicate an false discovery rate of 0.05 and absolute log fold change of one. (D) The RNA-seq gene expression for the peaks with a gain in R-loop signal upon knockdown shows no significant difference in distribution among treatments. (E) Venn diagrams for the overlaps between genes lying nearest the R-loop gain peaks and the genes whose RNA-seq expression is up-regulated or down-regulated upon the knockdown treatment relative to the control, as reported by DESeq2 v1.26.0 (absolute log2 fold change > 1 and false discovery rate < 0.05, Wald test).

  • Figure 3.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 3. R-loop gains are elevated close to neighboring genes on chromosomes with high gene density.

    (A) Distribution of distance to the neighboring gene from the DRIP-seq consensus peaks with a gain (red), loss (green), or unchanged (blue) R-loop signal upon each siRNA condition, measured through Bedtools v2.26.0 with Ensembl gene annotation. Definition of neighboring gene is the second nearest gene to the peak. (B) Ratio between the total amount of R-loop gain or loss peaks and the amount of unchanged peaks upon knockdown treatment, as a function of the mean distance from all consensus peaks to their corresponding neighboring (i.e., second nearest) genes, measured separately for each chromosome.

  • Figure 4.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 4. Distribution and coverage of DRIP-seq peaks across genomic locations.

    (A) Distribution of the peaks with no change in R-loop signal upon any of the knockdown treatments (left; Unchanged) and of the gain peaks upon each of the siDDX5, siXRN2, and siPRMT5 cells. The percentage of the R-loop gain peaks are distributed as 5′-UTR and 3′-UTR; promoter and transcription start site (promoter-TSS); transcription termination site (TTS), and noncoding, intronic, and intergenic peaks. (B) Distribution in percentage of the gain peaks in the intersection of two of three of the knockdown treatments (left) or in all of the knockdown treatments simultaneously (right). (C, D, E, F, G, H) Genome-wide R-loop signal profile for each treatment compared with control (C, E, G) and normalized to siCTL (D, F, H). (C, D, E, F, G, H) DRIP-seq coverage was measured using all reads (C, D), reads overlapping with intergenic regions (E, F), and reads overlapping with introns (G, H) through NGS.PLOT v2.63. (I) Histogram of peaks lying nearer to the TSS or to the TTS of the nearest gene using the full list of consensus peaks (left) or the gain peaks relative to the control (right) measured through HOMER v4.11.1. (J) Percentage of peaks lying nearer to the TSS or to the TTS of the nearest gene. From left to right: full list of consensus peaks; peaks with a gain in siDDX5 only; peaks with a gain in both siDDX5 and siXRN2; peaks with a gain in both siDDX5 and PRMT5; and peaks with a gain in siDDX5, siXRN2, and siPRMT5 cells.

  • Figure S4.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure S4.

    R-loop peak enrichment near TSS/TTS or GC-rich areas. (A, B) Distribution of the distances nearest transcription start site or transcription termination site for the R-loop peaks with a gain, loss, or no change upon the knockdown treatments. (C) Histogram of position-dependent nucleotide frequencies for the 4 kb relative to the center of DRIP-seq peaks using the full list of consensus peaks (left) or the gain peaks relative to the control (right) measured through HOMER v4.11.1. (D, E) Distribution of GC-rich (D) or CpG (E) percentage for the genomic regions covering R-loop peak gains, losses, or no change for each of the knockdown treatments.

  • Figure 5.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure 5. R-loops in siDDX5 induce antisense transcription.

    (A) Read coverage of DRIP-seq and RNA-seq signal centered at the R-loop gain peak associated to the EGR1, ACTG1, RHOB, RB1CC1, SOGA1, STIL, and UBALD1 genes and the downstream of IER2 gene loci in siCTL and siDDX5 (absolute log2 fold change > 1 and false discovery rate < 0.1, Wald test). The bottom two tracks show the intergenic RNA-seq coverage of the pooled replicates. Whereas the gene expression decreased upon the knockdown treatment, the read coverage at the intergenic region adjacent to the transcription start site increased, thus indicating possible antisense transcription. Control cells are blue, siDDX5 cells are red. (B) Antisense expression was quantified by reverse transcription quantitative PCR (RT-qPCR) using cDNAs transcribed with random primers. The amplified fragments are located at the peak regions upstream of EGR1, ACTG1, RHOB, RB1CC1, SOGA1, STIL, and UBALD1 genes and the downstream of IER2 gene. The relative expression was normalized with GAPDH. The graph shows the average and SEM from four independent experiments. Statistical significance was assessed using t test. *P < 0.05, ** P < 0.01, and ***P < 0.001. (C) Agarose gel image of the RT-PCR products amplified at the promoter region of the EGR1, ACTG1, RHOB, RB1CC1, SOGA1, STIL1, and UBALD1 genes and downstream of IER2 gene. The RT primers used are sense primers, which hybrid with the antisense strand corresponding to the gene; antisense primer, annealing to the sense strand; random primers; and no primer. DNA markers are shown on the left are in base pairs.

  • Figure S5.
    • Download figure
    • Open in new tab
    • Download powerpoint
    Figure S5. R-loops in siXRN2 and siPRMT5 induce antisense transcription. Read coverage of DRIP-seq and RNA-seq signal centered at the R-loop gain peak associated to the RPLP1, RPP25L, FBXW5, and PMEPA1 gene loci in siCTL and siXRN2 (absolute log2 fold change > 1 and false discovery rate < 0.1, Wald test) as well as to CUL3, NYAP1, CPT1A, and NACC2 gene loci in siCTL and siPRMT5 (absolute log2 fold change > 1 and false discovery rate < 0.1, Wald test).

    The bottom two tracks show the intergenic RNA-seq coverage of the pooled replicates. Whereas the gene expression decreased upon in the knockdowns, the read coverage at the intergenic region adjacent to the transcription start site increased, thus indicating possible antisense transcription. Control cells are blue, siXRN2 cells are green, and siPRMT5 cells are magenta.

Tables

  • Figures
  • Supplementary Materials
    • View popup
    Table 1.

    R-loop genomic coverage.

    siRNATypeSpace (Mb)
    DDX5Gain5.33
    DDX5Loss0.43
    DDX5Unchanged129.26
    XRN2Gain5.70
    XRN2Loss0.18
    XRN2Unchanged129.14
    PRMT5Gain15.29
    PRMT5Loss0.40
    PRMT5Unchanged119.33
    • View popup
    Table 2.

    Primers for reverse transcription (RT)-PCR analysis.

    EGR1RT sense: 5′-CCCTGTTCGCGTTCGGCCCC-3′
    RT antisense: 5′-GCTCGGTGCTGCCCCCTGGAG-3′
    PCR forward: 5′-CACCCCCTGCTTCCTTCTCC-3′
    PCR reverse: 5′-CGACGCAGTGAGCACGAACT-3′
    ACTG1RT sense: 5′-CGGAGCAGAACGTAG-3′
    RT antisense: 5′-GCCCAGAATCTCCGG-3′
    PCR forward: 5′-GTGTCCCTCGGTGTGTGACG-3′
    PCR reverse: 5′-CGGGCAAGGCTGTCAGGTAT-3′
    RHOBRT sense: 5′-CGGGACTTGGAAGAG-3′
    RT antisense: 5′-GCTCTGGCGGTACCC-3′
    PCR forward: 5′-GGGGCCCTAAACCACAGGAG-3′
    PCR reverse: 5′-GCCCCTCTTCCTGGCAAACT-3′
    RB1CC1RT sense: 5′-CGGGACTTGGAAGAG-3′
    RT antisense: 5′-GCTTGTTCCCCTCAG-3′
    PCR forward: 5′-TCCCAACCATTAGGGTGCTCA-3′
    PCR reverse: 5′-GCGGCACCATTTCTCAGACC-3′
    SOGA1RT sense: 5′-GAGATGGAGTCTAGC-3′
    RT antisense: 5′-CAGGAGTTCGAGACC-3′
    PCR forward: 5′-ACCTCGGCTCACTGCAACCT-3′
    PCR reverse: 5′-CCAACATGATGAAACCCCGTCT-3′
    STILRT sense: 5′-TTGAACTCGGGAGGC-3′
    RT antisense: 5′-CGCGCTCGACCAATC-3′
    PCR forward: 5′-GTTCTTCGGGTGTCCGCTTC-3′
    PCR reverse: 5′-CGGCGCTCCAGGATCAAG-3′
    UBALD1RT sense: 5′-TAGAGACGGTTTGAC-3′
    RT antisense: 5′-TTCCTGGCCCTGACC-3′
    PCR forward: 5′-GTCCTGGGCCTAGGCAATCC-3′
    PCR reverse: 5′-GGGAGCGAATTTCGGAAACC-3′
    IER2RT sense: 5′-CCGGTTACCACGTGG-3′
    RT antisense: 5′-TGATACTGTAGGGCC-3′
    PCR forward: 5′-CGGGCATTCCCTAACTGGTG-3′
    PCR reverse: 5′-GTGCAATCGATCCCCAGCTC-3′
    • View popup
    Table 3.

    Primers pairs for antisense RT-qPCR analysis.

    EGR15′-AGGCTCGGGGTGAGGAGTGT-3′
    5′-CGACGCAGTGAGCACGAACT-3′
    ACTG15′-GTGTCCCTCGGTGTGTGACG-3′
    5′-CAACAGACCCACCCGGACTC-3′
    RHOB5′-GCCAGGAAGAGGGGCAATTC-3′
    5′-GTCCGGGAGCTGGCTGTCT-3′
    RB1CC15′-TCCCAACCATTAGGGTGCTCA-3′
    5′-CGCCACAACCACGTTTTCAG-3′
    SOGA15′-ACCTCGGCTCACTGCAACCT-3′
    5′-CAAATTAGCCGGGCGTGGTA-3′
    STIL5′-GTTCTTCGGGTGTCCGCTTC-3′
    5′-CGCAATGGAAAGCCCAGCTA-3′
    UBALD15′-TCCTCGGACCCCGAGTAGGT-3′
    5′-GGGAGCGAATTTCGGAAACC-3′
    IER25′-CGGGCATTCCCTAACTGGTG-3′
    5′-AAAGCCCCGATCTCCCTGTC-3′
    • View popup
    Table 4.

    Primers for DRIP-qPCR validation.

    FOS5′-CCTGCAAGATCCCTGATGACCT-3′
    5′-AGGGTGAAGGCCTCCTCAGACT-3′
    KLF25′-GACAACAGTGGGGAGTGGACCTT-3′
    5′-CTGAGGGATCCTTGCCCTACATC-3′
    JUNB5′-CCGGATGTGCACTAAAATGGAAC-3′
    5′-AGTCGTGTAGAGAGAGGCCACCA-3′
    CTNNB15′-GCCATTTTAAGCCTCTCGGTCTG-3′
    5′-CTCCTCAGACCTTCCTCCGTCTC-3′
    LY6E5′-GAAGGCTGCTGAGTTTCCTCCTC-3′
    5′-GCTTCTCTCCTGACCCACTCCTC-3′
    SNHG125′-CTGGGACTATAAGCACGCACCAC-3′
    5′-TTGGGGTCAGGAGTTCAAGACTG-3′
    SOWAHC5′-GCTAGCCTTCTGGGAAAAGTGGA-3′
    5′-GAAGTGGAGGGCAGAGAAGAGGT-3′
    RPS23-15′-TTAGTCGGTTCAGGGCAACTTGA-3′
    5′-CTAAGACACTCGCCTCACCTGGA-3′
    RPS23-25′-GTTCATGCCTGTAATCCCAGCAC-3′
    5′-GTATGACTTTGCTGCCCAGGATG-3′

Supplementary Materials

  • Figures
  • Tables
  • Supplemental Data 1.

    Consensus R-loop peaks among the siLuciferase (siCTL), siDDX5, siXRN2, and siPRMT5 conditions.[LSA-2020-00762_Supplemental_Data_1.xlsx]

  • Supplemental Data 2.

    RNA-seq gene expression in siLuciferase (siCTL), siDDX5-, siXRN2-, and siPRMT5-transfected U2OS cells.[LSA-2020-00762_Supplemental_Data_2.xlsx]

  • Supplemental Data 3.

    Antisense intergenic transcription hits identified in the neighborhood of R-loop gain peaks.[LSA-2020-00762_Supplemental_Data_3.xlsx]

PreviousNext
Back to top
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Life Science Alliance.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Genome-wide R-loop analysis defines unique roles for DDX5, XRN2, and PRMT5 in DNA/RNA hybrid resolution
(Your Name) has sent you a message from Life Science Alliance
(Your Name) thought you would like to see the Life Science Alliance web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
DRIP-seq analysis of siDDX5, siXRN2 and siPRMT5 U2OS cells
Oscar D Villarreal, Sofiane Y Mersaoui, Zhenbao Yu, Jean-Yves Masson, Stéphane Richard
Life Science Alliance Aug 2020, 3 (10) e202000762; DOI: 10.26508/lsa.202000762

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
DRIP-seq analysis of siDDX5, siXRN2 and siPRMT5 U2OS cells
Oscar D Villarreal, Sofiane Y Mersaoui, Zhenbao Yu, Jean-Yves Masson, Stéphane Richard
Life Science Alliance Aug 2020, 3 (10) e202000762; DOI: 10.26508/lsa.202000762
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
Issue Cover

In this Issue

Volume 3, No. 10
October 2020
  • Table of Contents
  • Cover (PDF)
  • About the Cover
  • Masthead (PDF)
Advertisement

Jump to section

  • Article
    • Abstract
    • Introduction
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgements
    • References
  • Figures & Data
  • Info
  • Metrics
  • Reviewer Comments
  • PDF

Subjects

  • Systems & Computational Biology

EMBO Press LogoRockefeller University Press LogoCold Spring Harbor Logo

Content

  • Home
  • Newest Articles
  • Current Issue
  • Archive
  • Subject Collections

For Authors

  • Submit a Manuscript
  • Author Guidelines
  • License, copyright, Fee

Other Services

  • Alerts
  • Twitter
  • RSS Feeds

More Information

  • Editors & Staff
  • Reviewer Guidelines
  • Feedback
  • Licensing and Reuse
  • Privacy Policy

ISSN: 2575-1077
© 2021 Life Science Alliance LLC

Life Science Alliance is registered as a trademark in the U.S. Patent and Trade Mark Office and in the European Union Intellectual Property Office.