Genome-wide mapping of native co-localized G4s and R-loops in living cells
eLife assessment
This study describes a useful antibody-free method to map both G-quadruplexes and R-loops in vertebrate cells independently of the BG4 and S9.6 antibodies. It also reveals that the helicase Dhx9 can affect the self-renewal and differentiation capacities of mESCs, perhaps by regulating co-localized G4s and R-loops. The datasets provided might constitute a good starting point for future functional studies, and although the strength of the evidence that DHX9 interferes with the ability of mESCs to differentiate by regulating directly the stability of either G4s or R-loops has been improved compared to a previous version, it is still incomplete.
https://doi.org/10.7554/eLife.99026.3.sa0Useful: Findings that have focused importance and scope
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Incomplete: Main claims are only partially supported
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
The interplay between G4s and R-loops are emerging in regulating DNA repair, replication, and transcription. A comprehensive picture of native co-localized G4s and R-loops in living cells is currently lacking. Here, we describe the development of HepG4-seq and an optimized HBD-seq methods, which robustly capture native G4s and R-loops, respectively, in living cells. We successfully employed these methods to establish comprehensive maps of native co-localized G4s and R-loops in human HEK293 cells and mouse embryonic stem cells (mESCs). We discovered that co-localized G4s and R-loops are dynamically altered in a cell type-dependent manner and are largely localized at active promoters and enhancers of transcriptional active genes. We further demonstrated the helicase Dhx9 as a direct and major regulator that modulates the formation and resolution of co-localized G4s and R-loops. Depletion of Dhx9 impaired the self-renewal and differentiation capacities of mESCs by altering the transcription of co-localized G4s and R-loops -associated genes. Taken together, our work established that the endogenous co-localized G4s and R-loops are prevalently persisted in the regulatory regions of active genes and are involved in the transcriptional regulation of their linked genes, opening the door for exploring broader roles of co-localized G4s and R-loops in development and disease.
Introduction
Genomic DNA can form various types of non-B secondary structures, including G-quadruplexes (G4s), R-loops, Z-DNA, i-motifs, Cruciform, and others (Matos-Rodrigues et al., 2023). Among them, G4s and R-loops are two of the best studied. G4s are built by stacked guanine tetrads connected via Hoogsteen hydrogen bonds and can be formed by intra- or inter-molecular folding of the tetramers (Panyutin et al., 1989; Sen and Gilbert, 1988; Sundquist and Klug, 1989; Williamson et al., 1989). R-loops are three-stranded structures containing a DNA-RNA hybrid and a displaced single-stranded DNA (Aguilera and García-Muse, 2012; Xu and Clayton, 1996). Both G4s and R-loops are involved in key biological processes, including transcription, replication, genomic instability, class switch recombination in B cells, DNA damage and repair, and telomere maintenance (García-Muse and Aguilera, 2019; Hänsel-Hertsch et al., 2017; Varshney et al., 2020; Yang et al., 2023).
R-loops appear to have a strong sequence preference with high G/C ratios (Ginno et al., 2013; Ginno et al., 2012). Reports about the interplay between R-loops and G4s are emerging. Specific G4 ligands stabilized G4s and simultaneously increase R-loop levels within minutes in human cancer cells, which finally induced DNA damage (De Magis et al., 2019). Reactive oxygen species (ROS) have been reported to induce G4 and R-loop formation at transcriptionally active sites, and their inter-regulation is essential for the DNA repair (Tan et al., 2020). In a reconstituted eukaryotic DNA replication system, the interplay of R-loops and G4s was shown to impact replication fork progression by inducing fork stalling (Kumar et al., 2021). Single-molecule fluorescence studies showed the existence of a positive feedback mechanism of G4 and R-loop formation during transcription, where the transcription-induced R-loop precedes and facilitates G4 formation in the non-template strand, and in turn, G4 promotes the R-loop formation in the following rounds of transcription (Lee et al., 2020; Lim and Hohng, 2020). Wulfridge et al. reported that the architectural protein CCCTC binding factor (CTCF)-bound sites are enriched for R-loops and G4s which facilitate CTCF binding to promote chromatin looping interactions (Wulfridge et al., 2023).
Detection of G4s and R-loops have been largely based on the use of a single-chain variable fragment (scFv) BG4 for G4s and a monoclonal antibody S9.6 for R-loops. In recent years, these two antibodies have been coupled with deep sequencing to genome-widely detect G4s and R-loops (Galli et al., 2022; Ginno et al., 2012; Hänsel-Hertsch et al., 2016; Jiang et al., 2023; Lyu et al., 2021). Using BG4 and S9.6 -based CUT & Tag, G4s, and R-loops showed a high degree of co-occurrence in mESCs (Lyu et al., 2021). However, given that a group of helicases, RNA-binding factors, endonucleases, and DNA topoisomerases cooperate to actively dissolve G4s and R-loops restoring B-formed DNA duplexes (Hänsel-Hertsch et al., 2017; Varshney et al., 2020; Yang et al., 2023), a steady state equilibrium is generally set at low levels in living cells under physiological conditions (Miglietta et al., 2020) and thus the addition of high-affinity antibodies may pull the equilibrium towards folded states. Additionally, the specificity of the S9.6 antibody on R-loops has been questioned recently for accurate quantification and mapping of R-loops (Hartono et al., 2018; König et al., 2017; Phillips et al., 2013).
To understand the co-localized G4s and R-loops in living cells under physiological conditions, we sought to develop an in vivo strategy for G4 profiling based on the G4-hemin complex-induced proximal labeling and R-loop profiling based on the N-terminal hybrid-binding domain (HBD) of RNase H1. Recent studies showed that G4s could tightly form a complex with the cellular cofactor hemin both in vitro and in living cells, where hemin binds by end-stacking on the terminal G-quartets of G4s without affecting the folding of G4s (Gray et al., 2019; Stadlbauer et al., 2021). The G4-hemin complex has been shown to act as a peroxidase to catalyze oxidation reactions in the presence of hydrogen peroxide (H2O2) (Cheng et al., 2009; Lat et al., 2020; Yang et al., 2011). The H2O2-activated G4-hemin complex oxidizes the biotin tyramide to phenoxyl redicals that covalently conjugate biotin to G4 itself and its proximal DNA within 10 nm (equivalent to approximately 31 bp) in vitro and in vivo (Einarson and Sen, 2017; Lat et al., 2020). Here, we have utilized the G4-hemin-mediated proximal biotinylation rection to develop a new method HepG4-seq (for high throughput sequencing of hemin-induced proximal labeled G4s) to map the genomic native G4s under physiological conditions. The HBD domain of RNase H1 mediates the specific recognition of DNA/RNA hybrid in a sequence-independent manner, which is a gold standard for R-loop recognition in the cell (Nowotny et al., 2008). The catalytically inactive RNase H1 or its HBD domain -based methods have been successfully used to identify genome-wide native R-loops (Chen et al., 2017; Wang et al., 2021). We have adapted the ‘GST-His6-2xHBD’-mediated CUT&Tag protocol (Wang et al., 2021) to develop the HBD-seq protocol in this study.
We have combined the HepG4-seq and HBD-seq to profile the genome-wide native co-localized G4s and R-loops with high signal-to-noise ratios in HEK293 cells and mESCs. We observed that the co-localized G4s and R-loops exhibit cell type-dependent distributions and are largely localized at active promoters and enhancers of transcriptionally active genes. We further showed that ~70% of the co-localized G4s and R-loops in mESCs were directly bound by the helicase Dhx9 and that depletion of Dhx9 significantly altered the levels of ~6200 co-localized G4s and R-loops bound by Dhx9. Furthermore, depletion of Dhx9 was shown to impair the self-renewal and differentiation capacities of mESCs by altering the transcription of co-localized G4s and R-loops -associated genes.
Results
Mapping of the native DNA G4 through the G4-hemin-mediated proximal biotinylation
The DNA G4-hemin complex could act as a mimic peroxidase to oxidize the biotin tyramide to phenoxyl radicals that can covalently conjugate biotin to G4 itself and its proximal DNA within ~30 bp in the presence of H2O2 (Cheng et al., 2009; Einarson and Sen, 2017; Lat et al., 2020). However, the efficiency of peroxidase-mediated biotinylation on DNA is limited using the substrate biotin tyramide (Zhou et al., 2019a). Recently, biotin aniline (Bio-An) has been shown to have superior labeling efficiency on DNA than biotin tyramide, when catalyzed by the engineered peroxidase APEX2 (Zhou et al., 2019a). The free heme concentration in normal human erythrocytes is 21±2 µM (Aich et al., 2015). To explore the G4-hemin-mediated biotinylation in the living cells, we treated HEK293 cells with 25 µM hemin and 500 µM Bio-An for 2 hr prior to activation with 1 mM H2O2 for 1 min, and then quenched the labeling reaction and performed the immunofluorescence staining using Alexa Fluor 647 conjugated recombinant streptavidin (Strep-647) that specifically recognizes biotin. As shown in Figure 1A, cells treated with hemin and Bio-An exhibited a robust fluorescence signal, while the absence of either hemin or Bio-An almost completely abolished the biotinylation signals, suggesting a specific and active biotinylation activity. To understand whether the addition of hemin disturbs the formation of G4s, we performed the BG4 CUT&Tag-seq using the recombinant BG4 and Tn5 on HEK293 cells treated with and without 25 µM hemin (Figure 1—figure supplement 1A, B). The heatmap and profile plot analysis showed similar BG4 CUT&Tag signals between the hemin-treated and control samples (Figure 1—figure supplement 1B). There were only 174 BG4 CUT&Tag peaks with significantly differential signals between the hemin-treated and control samples (Figure 1—figure supplement 1C). These data suggest that the hemin treatment condition we used does not significantly affect G4 folding. Therefore, hemin-induced proximal biotinylation of G4s could be utilized to mark the native G4s in living cells.
The recombinant streptavidin monomer (mSA) combines the streptavidin and rhizavidin sequences to achieve specific monovalent detection of biotin or biotinylated molecules with a high affinity (Kd = 2.8 nM) (Lim et al., 2013). The Moon-tag system consists of a 15-amino acid peptide GP41-tag and a 123-amino acid anti-GP41-tag nanobody with an affinity of ~30 nM in vitro (Boersma et al., 2019). We expressed and purified the recombinant mSA fused with the anti-GP41 nanobody (mSA-scFv), and the recombinant Tn5 fused with the GP41-tag and protein G (GP41-pG-Tn5) from E. coli (Figure 1—figure supplement 1A). To map the hemin-induced biotinylated G4s, we developed a new method HepG4-seq (Figure 1B), where mSA-scFv recognizes the biotinylated G4s and recruits the transposase Tn5 to achieve ‘Cleavage Under Targets and Tagmentation (CUT&Tag).’ Deep sequencing analysis of the biotinylated G4s fragments identified 6,799 consensus peaks from two independent biological repeats in HEK293 cells, where the signals were dramatically diminished in HEK293 cells without treatment with hemin and Bio-An, suggesting the specificity of HepG4-seq (Figure 1C, Supplementary file 1). Several representative HepG4s-seq-identified G4 peaks were shown in Figure 1D. Genomic distribution analysis showed that G4s are mainly localized in promoters (38.7%) and gene bodies (47.1%) (Figure 1E).
We also evaluated the HepG4-seq-identified peaks using a G4-forming sequences (PQS) predication tool pqsfinder which has been shown to have 96% accuracy on ~400 known and experimentally observed G4 structures (Hon et al., 2017). The peaks identified by HepG4-seq overlap quite well with the center of pqsfinder max Scores that report the PQS quality (Figure 1C and D). The motif enrichment analysis by HOMER (Heinz et al., 2010) revealed a high prevalence of G-rich sequences in HepG4-seq peaks (Figure 1F). All above further validated the specificity of HepG4-seq in capturing G4s.
Induction of DNA G4s by inhibiting G4 resolving helicase
The RecQ-like helicases Bloom syndrome protein (BLM) and Werner syndrome ATP-dependent helicase (WRN) are the first recognized and the best-characterized DNA G4-resolving mammalian helicases (Fry and Loeb, 1999; Mendoza et al., 2016; Mohaghegh et al., 2001). The small molecules ML216 and NSC617145 are selective and cell-permeable inhibitors of BLM and WRN, respectively, by inhibiting their ATPase activity (Aggarwal et al., 2013a; Aggarwal et al., 2013b; Nguyen et al., 2013). To investigate the effect of BLM or WRN inhibition on native G4s, we treated HEK293 cells with ML216 or NSC617145 for 16 hr and then labeled cellular G4s by hemin-G4-induced biotinylation in living cells. Immunofluorescence staining of the treated cells using Strep-647 showed that the treatment of ML216 or NSC617145 remarkably elevated signals of native G4s (Figure 1G). Furthermore, we performed HepG4-seq on HEK293 cells treated with ML216 or NSC617145. Notably, HepG4-seq identified 77,003 peaks from ML216- or NSC617145- treated HEK293 cells, and ~70,000 new G4 peaks were induced by inhibition of BLM or WRN (Supplementary file 1). The signals of G4s detected by HepG4-seq were significantly increased after inhibiting BLM or WRN (Figure 1H and I) (Mann-Whitney test, p<2.2E-16). Representative G4s peaks are shown in Figure 1J. In addition, the top enriched motifs in the 70000 extra HepG4-seq peaks are G-riched (Figure 1—figure supplement 1E). Taken together, these data suggest that HepG4-seq is able to efficiently detect dynamic native G4s.
Mapping of native co-localized G4s and R-loops in HEK293 cells
The HBD domain of RNase H1 has been demonstrated as a DNA/RNA hybrid recognition sensor and applied to identify genome-wide native DNA/RNA hybrids using the recombinant GST-His6-2xHBD coupling with Tn5-based CUT&Tag (Nowotny et al., 2008; Wang et al., 2021). Given that the GST-fusion proteins are prone to form variable high molecular-weight aggregates and these aggregates often undermine the reliability of the fusion proteins (Deceglie et al., 2014; Ki and Pack, 2020), we produced the recombinant two copies of HBDs fused with EGFP and V5-tag (HBD-V5) (Figure 1—figure supplement 1A) and used the anti-V5 tag antibody instead of the anti-His tag antibody for the CUT&Tag-seq. We call this modified protocol as HBD-seq for mapping the native R-loops in cells (Figure 2A). We performed the HBD-seq on HEK293 cells, and revealed 42,488 consensus native R-loops peaks with a high signal-to-noise ratio while the HBD-seq signals were dramatically diminished in HEK293 cells treated with the RNases prior to HBD-seq (Figure 2B), suggesting the specificity of HBD-seq in detecting native R-loops.
We then analyzed the regions co-occupied by both HepG4-seq-identified G4s and HBD-seq-identified R-loops, and revealed 5030 native co-localized peaks in HEK293 cells, ranging in size from 100 bp to ~1.5 kb (Figure 2C and D, Figure 1—figure supplement 1D). 73.8% of these co-localized peaks are localized at promoters, 5’UTR, exon1, and intron 1 (Figure 2E). When we performed a metagene analysis of these co-localized peaks, a distinct peak was detected around the transcription start site (TSS) (Figure 2F). Representative co-localized peaks are shown in the Figure 2G. The motifs enrichment analysis by HOMER (Heinz et al., 2010) showed that G-rich sequences are highly enriched in the co-localized peaks (Figure 2H).
The co-localized G4s and R-loops-mediated transcriptional regulation in HEK293 cells
The predominant distribution of co-localized peaks around TSS implies that they may participate in the transcriptional regulation of their associated genes. RNA-seq analysis revealed that the RNA levels of co-localized G4s and R-loops-associated genes are significantly higher than all genes, G4s or R-loops- associated genes with the Mann-Whitney test p<2.2E-16 (Figure 3A). Different from G4s and R-loops, the co-localized G4s and R-loops are mainly localized within 1 kb of the TSS of transcriptionally active genes (~60% peaks with FPKM ≥ 5) (Figure 3B). To investigate the transcriptional regulation of co-localized G4s and R-loops in living cells, we performed the RNA-seq on HEK293 cells treated with and without ML216 or NSC617145 and then analyzed the differential gene expression using DESeq2 (Love et al., 2014). As a result, hundreds of genes were linked to co-localized G4s and R-loops with increased G4 signals (at least 1.5 foldchange) and at the same time exhibited significant changes in expression levels upon inhibition of BLM or WRN in HEK293 cells (Figure 3C), suggesting that co-localized G4s and R-loops could regulate the transcription of their associated genes. Distribution analysis showed that these differential co-localized G4s and R-loops are mainly localized in the promoter-TSS (Figure 3D). Among the differential genes, 125 genes were co-regulated by both ML216 and NSC617145 (Figure 3E), suggesting that BLM and WRN could co-regulate the transcription of genes by resolving G4s. Gene ontology (GO) analysis showed that co-localized G4s & R-loops-regulated genes in HEK293 cells are mainly involved in cell cycle regulation, DNA/mRNA metabolic regulation, DNA damage response, chromatin binding, kinase binding, cell-substrate junction, et al (Figure 3F).
Mapping of native co-localized G4s and R-loops in mESCs
mESCs are pluripotent stem cells that c differentiate into various types of cells of three germ lineages (Murry and Keller, 2008; Young, 2011). To understand the regulatory roles of co-localized G4s and R-loops in mESCs, we performed the HepG4-seq and HBD-seq on mESCs and finally uncovered 68,482 native overlapping peaks in mESCs, ranging in size from 100 bp to ~2 kb (Figure 4A–D, Figure 4—figure supplement 1, Supplementary file 2). Notably, unlike HEK293 cells, a large number of native G4s (95,128) were identified by HepG4-seq in mESCs (Figure 4A, Supplementary file 2), which well overlap with the PQS predicted by pqsfinder (Figure 4B), suggesting that native G4s exhibit obvious cell type-specific distribution.
For the genomic distribution of co-localized G4s and R-loops, unlike HEK293 cells, only 34.2% peaks are localized in promoters, exon1, and intron1 while 25.9% peaks in intergenic regions (Figure 4E). The distinct number and localization feature of co-localized G4s and R-loops in HEK293 cells and mESCs shows the cell type-specific distribution. The metagene analysis of overlapping peaks exhibited a distinct peak around TSS (Figure 4F), suggesting the potential for transcriptional regulation. The motifs enrichment analysis found that G-rich sequences are highly enriched in these overlapping peaks, similar to those in HEK293 cells (Figure 4G). Representative peaks were found in several key regulatory genes of mESCs (Figure 4H).
Characterization of native co-localized G4s and R-loops in mESCs
Similar to HEK293 cells, the RNA levels of co-localized G4s and R-loops-associated genes were seen to be significantly higher in mESCs (Mann-Whitney test p<2.2E-16) (Figure 5A). However, unlike in HEK293 cells (Figure 3B), the overlapping peaks in mESCs are mainly localized in the proximal promoters (1 kb from TSS, 9063 peaks with FPKM ≥ 5) and the region 5–50 kb from the TSS of transcriptionally active genes (32,690 peaks with FPKM ≥ 5) (Figure 5B), suggesting that co-localized G4s and R-loops are possibly distributed in active promoters or enhancers. To test this idea, we analyzed the co-localization of G4 (HepG4-seq), R-loop (HBD-seq), ChIP-seq signals of multiple chromatin markers and RNA polymerase II with the phosphorylated serine 5 at its CTD domain (RNAP) that marks the transcriptionally initiated RNA Polymerase II (Hsin and Manley, 2012). As a result, co-localized G4s and R-loops were observed to well overlap with active chromatin markers (H3K4me3, H3K27ac, H3K36me3, H3K4me1) and RNAP but not the repressed chromatin marker H3K27me3 (Figure 5C).
Extensive studies define promoters into active, bivalent, and repressed states based on patterns of H3K4me3 and H3K27me3 (Figure 5D); enhancers are defined as active, poised, and unmarked states based on patterns of H3K27ac and H3K4me1 (Figure 5D; Atlasi and Stunnenberg, 2017; Bernstein et al., 2006; Bibikova et al., 2008; Calo and Wysocka, 2013; Heintzman et al., 2007). Notably, 20,741 and 19,726 co-localized G4s and R-loops are found in promoters and enhancers, respectively; 18,496 peaks are seen in active promoters; 15,787 peaks are seen in active enhancers (Figure 5D, Supplementary file 2). The co-localized G4s and R-loops in active promoters show high and almost equal signals of G4s and R-loops, and enrich H3K4me3, H3K27ac, and RNAP (Figure 5E). The co-localization of G4s, R-loops and RNAP at active promoters suggests that co-localized G4s and R-loops are likely linked to promoter-associated nascent RNAs (Core et al., 2008; Li and Fu, 2019; Preker et al., 2008; Seila et al., 2008). Interestingly, a medium level of H3K4me1 is present in a bimodal pattern beside co-localized G4s and R-loops at active promoters (Figure 5E). The co-localized G4s and R-loops in bivalent promoters exhibit low signals of both G4s and R-loops, and overlap with a high level of H3K27me3, a medium level of H3K4me1, and a low level of H3K4m3 (Figure 5E). The co-localized G4s and R-loops in active enhancers exhibit sharp peaks and well overlap with H3K27ac, H3K4me1, H3K4me3, and RNAP; the co-localized G4s and R-loops in poised enhancers enrich the active histone marks H3K4me1 and H3K4me3, and the repressive mark H3K27me3 while the signal of RNAP is very low (Figure 5F). The unmarked enhancers-associated co-localized G4s and R-loops show low levels of all marks tested. Given that enhancer RNAs (eRNAs) have been widely identified as non-coding RNAs in enhancers and are functionally important for enhancer activity (Andersson et al., 2014; Kim et al., 2010; Sigova et al., 2015), the co-occupacy of eRNAs are likely involved in the formation of co-localized G4s and R-loops in enhancers and shed light on the new regulatory mechanism of eRNA action.
Modulation of co-localized G4s and R-loops by the helicase Dhx9
Dhx9 (also known as RNA Helicase A) is a versatile helicase capable of directly resolving R-loops and G4s or promoting R-loop formation by unwinding secondary structures in the nascent RNA strand (Chakraborty and Grosse, 2011; Chakraborty et al., 2018; Cristini et al., 2018; Matsui et al., 2020; Tang et al., 2022; Yuan et al., 2021), and has been reported to play roles in DNA replication, transcription, translation, RNA processing and transport and maintenance of genomic stability (Aktaş et al., 2017; Aratani et al., 2001; Chellini et al., 2022; Jain et al., 2013; Tang et al., 2022). Thus, Dhx9 is a promising regulator of co-localized G4s and R-loops.
To investigate the role of Dhx9 in modulating co-localized G4s and R-loops, we generated the Dhx9 knockout mESCs (dhx9KO) by CRISPR/Cas9-medidated gene editing (Ran et al., 2013). The depletion of Dhx9 in the dhx9KO mESC clone was confirmed by western blot assay (Figure 6A) and immunofluorescence staining (Figure 6B). We next examined the G4 an R-loop levels in the dhx9KO mESCs by performing HepG4-seq and HBD-seq. Notably, compared to the wildtype mESCs, a large amount of G4s or R-loops within co-localized G4s and R-loops exhibited significantly up-regulated or down-regulated signals in the dhx9KO mESCs (Figure 6C, Figure 6—figure supplement 1A, Supplementary file 2), suggesting that Dhx9 could unwind or promote the formation of co-localized G4s and R-loops in mESCs. Interestingly, only a small proportion of co-localized G4s and R-loops displayed differential G4s and R-loops at the same time in the dhx9KO mESCs (Figure 6D, Figure 6—figure supplement 1B), suggesting that Dhx9 cannot simultaneously unwind or promote G4s and R-loops within co-localized G4 and R-loop regions and that multiple helicases or regulators are required for modulating these regions.
Given that co-localized G4s and R-loops have been shown to be enriched in active promoters and enhancers (Figure 5D–F), the loss-of-Dhx9-induced differential co-localized G4s and R-loops preferentially localize in the active and bivalent promoters and all three types of enhancers (Figure 6E). To explore the effect of Dhx9 on the transcription of co-localized G4s and R-loops -associated genes, we performed the RNA-seq on wild-type and dhx9KO mESCs. Differential gene expression analysis revealed 1647 significantly up-regulated genes and 1916 significantly down-regulated genes in the absence of Dhx9 (Figure 6—figure supplement 1C, D). Importantly, loss of Dhx9 resulted in hundreds of G4-, R-loop-, and G4&R-loop-associated genes with significantly differential expression (Figure 6F), suggesting that Dhx9 could regulate transcription by modulating G4s, R-loops, and co-localized G4s and R-loops. Representative Dhx9-regulated locus are shown in Figure 6G. GO analysis showed that co-localized G4s and R-loops-associated genes that showed differential expression after knocking out Dhx9 are mainly involved in negative regulation of cell differentiation, head development, positive regulation of cell motility, DNA-binding transcription activator activity, pattern specification process, embryonic organ morphogenesis, et al., suggesting that Dhx9 may regulate the cell fate of mESCs by modulating co-localized G4s and R-loops (Figure 6H). Coinciding with GO analysis, Dhx9 knockout in the mouse causes embryonic lethality, and Dhx9 knockdown leads to large structural changes in chromatin and eventually cell death (He et al., 2008; Zhang et al., 2004). Heterozygous loss-of-function variants of DHX9 are associated with neurodevelopmental disorders in humans (Yamada et al., 2023).
Characterization of co-localized G4s and R-loops directly bound by Dhx9
Tens of helicases or regulators have been reported to directly resolve or stabilize G4s or R-loops (Mendoza et al., 2016; Varshney et al., 2020; Yang et al., 2023). Interestingly, the loss of Dhx9 caused 30 of these helicases/regulators to be significantly differentially expressed (Figure 7A) and Dhx9 physically interacts with at least 10 of them based on the STRING protein-protein interaction network database (Szklarczyk et al., 2023; Figure 6—figure supplement 1E). These data suggest that Dhx9 could also indirectly modulate G4s and R-loops by affecting other helicases or regulators. Thus, to explore the direct target co-localized G4s and R-loops of Dhx9, we performed the CUT&Tag-seq using the Dhx9 antibody (Kaya-Okur et al., 2019) and revealed 54,982 Dhx9 binding peaks in wild-type mESCs (Figure 7B). Notably, 65.5% Dhx9 binding peaks well overlapped with 69.9% co-localized G4s and R-loops in mESCs (Figure 7C and D), suggesting that Dhx9 is a direct and major regulator of co-localized G4s and R-loops in mESCs. Motif analysis showed that G-rich sequences are highly enriched in the Dhx9 binding peaks overlapping with co-localized G4s and R-loops in mESCs (Figure 7E), further demonstrating that Dhx9 directly bind to co-localized G4s and R-loops that harbor G-rich sequences as shown in Figure 4G.
We next compared the Dhx9-bound co-localized G4s and R-loops in wild-type and dhx9KO mESCs, and identified 1382 significantly up-regulated peaks (823 with increased G4s and 559 with increased R-loops) and 4789 significantly down-regulated peaks (2278 with decreased G4s and 2511 with decreased R-loops) (Figure 7F, Figure 6—figure supplement 1F), accounting for ~50–75% of differential co-localized G4s and R-loops in absence of Dhx9 (Figure 6C). Analysis of the genomic distribution of Dhx9-bound differential co-localized G4s and R-loops found that these peaks are mainly localized in the active and bivalent promoters and all three types of enhancers (Figure 7G). The Dhx9-bound differential co-localized G4s and R-loops were linked to 852 genes with significantly differential expression (Figure 6—figure supplement 1G, Supplementary file 2), which are enriched in GO terms related to pattern specification, cell junction organization, brain development, negative regulation of cell differentiation, mesenchyme/mesoderm development, embryonic morphogenesis, et al. (Figure 7H). Several key regulators of mouse embryonic stem cell and embryonic development, such as Nanog, Lin28a, Bmp4, Wnt8a, Gata2, and Lef1, were shown to be transcriptionally regulated by Dhx9 through direct modulation of their associated G4s and R-loops (Figure 7I). These data suggest that Dhx9 significantly contributes to transcriptional regulation of co-localized G4s and R-loops-associated genes.
Dhx9 regulates the cell fate of mESCs
To understand the role of Dhx9 in regulating the cell fate of mESCs, we first examined the RNA levels of several key genes that maintain the pluripotency of mESCs by quantitative RT-PCR (qRT-PCR) and found that the RNA levels of Lin28a and Oct4 were significantly decreased and the RNA level of Nanog was significantly increased when Dhx9 was knocked out (Figure 8A). Western blot assay showed that dhx9KO mESCs produce obviously lower level of Lin28a protein than wild-type mESCs, consistent with its RNA level, whereas in contrast to the RNA level, the protein level of Nanog was significantly decreased in dhx9KO mESCs, suggesting that Dhx9 directly or indirectly modulates the translation of Nanog (Figure 8B). In line with the western blot assay, immunofluorescence staining of WT and dhx9KO mESCs showed that the loss of Dhx9 leads to reduced protein level of Nanog, but not Oct4, while dhx9KO mESCs exhibited normal morphology (Figure 8C). The mESCs can be maintained in a proliferative state for prolonged periods, which was known as ‘self-renewal’ (Liang and Zhang, 2013; Murry and Keller, 2008). Nanog and Lin28a have been reported to promote embryonic stem self-renewal (Chambers et al., 2003; Mitsui et al., 2003; Xu et al., 2009). Coinciding with reduced levels of Nanog and Lin28a proteins, dhx9KO mESCs was shown to modestly arrest at the S phase of the cell cycle (Figure 8D) and exhibited significantly attenuated proliferation capacity (Figure 8E), suggesting that Dhx9 regulates the self-renewal of mESCs.
mESCs are pluripotent stem cells which are able to differentiate into three germ lineages (Murry and Keller, 2008; Young, 2011). In the absence of differentiation inhibitor LIF, mESCs cultured in suspension spontaneously form three-dimensional aggregates called embryoid bodies (EBs), which could recapitulate many aspects of early embryogenesis, including the induction of three early germ lineages (Simunovic and Brivanlou, 2017). To understand the role of Dhx9 in regulating the pluripotency of mESCs, we performed the EB assay using the wild-type and dhx9KO mESCs. As the EB differentiation progressed, loss of Dhx9 resulted in apparently smaller and fewer EBs than wild-type cells (Figure 8F). At the same time, we collected EBs at different days of EB differentiation and examined the RNA levels of well-known markers of three germ lineages by qRT-PCR. As shown in Figure 8G, all maker genes tested displayed significantly differential expression, suggesting that Dhx9 regulates the pluripotency of mESCs, which is in line with the GO enrichment results of Dhx9-regulated co-localized G4s and R-loops-associated genes in Figure 7H. Taken together, Dhx9 regulates the self-renewal and differentiation capacities of mESCs.
Comparisons of HepG4-seq and HBD-seq with previous methods
To compare the performance of HepG4-seq and BG4 CUT&Tag, we analyzed the data using the same bioinformatic pipeline. As shown in Figure 9A, 80% (5,459) and 71% (67,935) HepG4-seq peaks overlap with BG4 CUT&Tag peaks in HEK293 cells and mESCs, respectively. SEACR is a highly selective peak caller for CUT&Tag and CUT&Run (Meers et al., 2019). With the help of SEACR, BG4, CUT&Tag obtained a much larger number of peaks in HEK293 cells and mESCs than previously reported MACS2-identified peaks. Surprisingly, 70,954 and 103,303 consensus peaks were identified by BG4 CUT&Tag alone in HEK293 cells and mESCs, respectively (Figure 9A), suggesting that BG4 antibody may promote the folding of G4 sequences during BG4-seq or that BG4 antibody is more sensitive than G4-hemin-induced proximal biotinylation in recognizing native G4s. As shown in Figure 9B, the patterns of BG4 CUT&Tag peaks overlap well with those of PQS while HepG4-seq peaks exhibit a much sharper shape and better resolution around the center of PQS patterns, indicating that BG4 may promote the folding of G4 motifs. Interestingly, the signal intensities of the HepG4-seq peaks are significantly lower than those of BG4 CUT&Tag in HEK293 cells, but obviously higher than those of BG4 CUT&Tag in mESCs. In the in vitro hemin-G4-induced self-biotinylation assay, parallel G4s exhibit higher peroxidase activities than anti-parallel G4s (Lat et al., 2020). The signal intensity changes of HepG4-seq in different cells possibly reflect the dynamics of G4 conformation. In the future, people may need to combine HepG4-seq and BG4s-eq to carefully explain the endogenous G4s.
During the past decade, R-loops were mapped using either S9.6 monoclonal antibody or catalytically inactive ribonuclease H1 for specific DNA-RNA hybrids. First, we compared the HBD-seq with S9.6 CUT&Tag in mESCs using the same bioinformatic pipeline and SEACR for peak calling. As a result, HBD-seq and S9.6 CUT&Tag identified similar numbers of peaks (125,311 for HBD-seq, 154,993 for S9.6 CUT&Tag), and 45,989 HBD-seq peaks overlap with 88,036 S9.6 peaks, indicating that S9.6 peaks are possibly narrower than HBD-seq peaks. As shown in Figure 9D, the overlapping peaks showed similar heatmap patterns and signal intensities between HBD-seq and S9.6 CUT&Tag while, notably, peaks of S9.6 CUT&Tag exhibit a much sharper signals at the center of peaks, which reflect the difference of binding modules of HBD and S9.6. GST-2xHBD and S9.6 mediated CUT&Tag were shown to generate highly similar native R-loop profiles (Wang et al., 2021). HBD-seq, as an optimized version of GST-2xHBD CUT&Tag, still generated ~80,000 extra peaks compared to S9.6 CUT&Tag, suggesting that the status of the recombinant 2xHBD protein could substantially affect R-loop mapping. To understand the difference among inactive ribonuclease H1-based methods, we compared the HBD-seq, MapR (inactive ribonuclease H1-mediated CUT&Run) (Yan et al., 2019) and GST-2xHBD CUT&Tag (Wang et al., 2021) and re-analyzed the data using the same pipeline and peak caller SEACR. In HEK293 cells, HBD-seq, MapR, and GST-2xHBD CUT&Tag identified 42,488, 12,520, and 17,570 peaks, respectively; 5752 MaR peaks and 11,701 GST-2xHBD CUT&Tag peaks overlap with HBD-seq peaks (Figure 9E). Notably, the heatmap of all HBD-seq peaks showed an obviously stronger signal intensities and sharper peak shapes than MapR and GST-2xHBD CUT&Tag, suggesting that HBD-seq provides superior R-loop signals. In consideration of the difference in DNA-RNA hybrid binding affinity and/or specificity between HBD and S9.6, the HBD-seq and S9.6 CUT&Tag may need to be combined together to carefully explain the native R-loops.
Discussion
In this study, we developed the new method ‘HepG4-seq’ and optimized the RNase H1 HBD domain-based HBD-seq to robustly map endogenous G4s and R-loops, respectively, in living cells with high specificity. Using the HepG4-seq and HBD-seq, we systematically characterized the native co-localized G4s and R-loops in HEK293 cells and mESCs, and revealed that co-localized G4s and R-loops are dynamically altered in a cell type-dependent way and largely localized at active promoters and enhancers of transcriptional active genes. Small molecules-induced inhibition of helicases BLM or WRN resulted in a significant accumulation of G4s within co-localized G4s and R-loops and at the same time led to genes with significantly differential expression in HEK293 cells that are enriched in the processes related to cell cycle, DNA metabolic, DNA damage response, chromatin binding, et al. Furthermore, we characterized the helicase Dhx9 as a key regulator of co-localized G4s and R-loops which efficiently unwinds or promotes co-localized G4s and R-loops, and illustrated that depletion of Dhx9 significantly altered the transcription of co-localized G4s and R-loops-associated genes that are enriched in embryonic development, cell differentiation, and germ lineage development, et al. Therefore, the loss of Dhx9 apparently impaired the self-renewal and pluripotency of mESCs.
In this study, we utilized a low dosage of hemin, similar to the physiological concentration in normal human erythrocytes, to spark the peroxidase activity of endogenous G4s without significantly altering the levels of native G4s (Figure 1—figure supplement 1B, C) and further robustly biotinylated G4s themselves by G4-hermin complex-mediated proximity labeling in just 1 min in living cells (Cheng et al., 2009; Einarson and Sen, 2017; Lat et al., 2020; Li et al., 2016; Stadlbauer et al., 2021; Yang et al., 2011). In consideration of the high affinity and specificity, the recombinant streptavidin monomer (Lim et al., 2013) is able to recognize the biotinylated G4s with high sensitivity and specificity and thereby yield robust CUT&Tag signals with the help of Moon-tag system (Boersma et al., 2019). Therefore, our HepG4-seq strategy is able to robustly and specifically capture native G4s. In HEK293 cells, HepG4-seq uncovered 6799 consensus G4 peaks under wild-type status and 77,003 G4 peaks in the presence of BLM/WRN inhibitors, suggesting that HepG4-seq is capable of detecting endogenous G4s with high sensitivity. Yet, the hemin-G4 complex-induced peroxidase activities were reported to be variable between parallelled and anti-parallelled G4s in the in vitro assay (Lat et al., 2020), although only 2–4 synthesized short DNA fragments were used for the test. Given that a method of genome-widely distinguishing parallelled and anti-parallelled G4s is lacking, HepG4-seq and BG4-seq may be combined to characterize the conformation of G4s.
Notably, we also discovered that the native co-localized G4s and R-loops landscape is altered in a cell-dependent manner and that approximately 10 folds more peaks were observed in mESCs than in HEK293 cells, which reflects that co-localized G4s and R-loops have the potential to regulate the complex pluripotency network in mESCs. CTCF is a key regulator of genome organization and gene expression (Ong and Corces, 2014). Recently, Wulfridge et al. reported that CTCF-bound regions are enriched for both R-loops, and G4s and G4s associated with R-loops promote CTCF binding (Wulfridge et al., 2023). Interestingly, the enriched motif with the most significant p-value in mESCs co-localized G4s and R-loops (Figure 4G) well matches the motif of CTCF ChIP-seq, suggesting that co-localized G4s and R-loops may be able to modulate the CTCF binding. To be noticed, the co-localized G4s and R-loops were identified by bioinformatically intersection analysis. The heterogenous distribution between cells will give false positive co-occurrence patterns. Recently, multi-CUT&Tag, a method for simultaneous mapping of multiple chromatin proteins, has been reported (Bartosovic and Castelo-Branco, 2023; Gopalan et al., 2021; Meers et al., 2023). As the multi-CUT&Tag, we could combine the HepG4 and HBD-seq to simultaneously map the co-localized G4s and R-loops even in the single cell level.
Furthermore, while 47,857 co-localized G4s and R-loops are directly bound by Dhx9 in the wild-type mESCs (Supplementary file 2), only 4060 of them display significantly differential signals in the absence of Dhx9. In addition, a limited number of co-localized G4s and R-loops showed differential G4s and R-loops at the same time in the dhx9KO mESCs. These data suggest that redundant regulators exist. It is worth noting that depletion of Dhx9 significantly altered the transcription of 30 known G4s and/or R-loops helicases/regulators (Figure 7A) and that half of these helicases/regulators are able to establish physical interaction network (Figure 6—figure supplement 1E). Due to the redundancy and compensatory roles of G4 and R-loop regulators, a comprehensively mechanistic study of G4 and R-loop interplay in living cells is difficult and lacking in the field. A degron system-mediated simultaneous and/or stepwise degradation system of multiple regulators will help us elucidate the interplaying effects between G4s and R-loops.
Taken together, our study provides new insights into exploring the regulatory roles of co-localized G4s and R-loops in development and disease.
Materials and methods
Expression and purification of recombinant proteins
Request a detailed protocolTo obtain the plasmid expressing the recombinant protein GP41-pG-Tn5, we synthesized a DNA fragment encoding a single copy of GP41 peptide (Boersma et al., 2019) and then followed the Gibson assembly protocol (Gibson et al., 2009) to insert this fragment between 10xHis-tag and pG (GB1) of the plasmid pET28a_His-GB1-Tn5 (Xu et al., 2021) to get the plasmid pET28a_His-GP41-GB1-Tn5. For the plasmid expressing the recombinant protein HBD-V5, we synthesized the DNA fragment encoding 10xHis-tag, EGFP, two copies of HBD domain of RNase H1 with optimized codon, and 2xV5-tag (10xHis-EGFP-2xHBD-2xV5) and cloned it into pET28a to get the plasmid pET28a_His-EGFP-2xHBD-2xV5. For the plasmid expressing the recombinant protein BG4-V5, we synthesized the DNA fragment encoding BG4, 2xV5-tag, and 6xHis-tag (BG4-2xV5-His) and cloned it into pET28a to get the plasmid pET28a_ BG4-2xV5-His. For the plasmid expressing the recombinant protein mSA-anti-GP41scFv, we synthesized the DNA fragment encoding mSA (Lim et al., 2013), anti-GP41scFv (Boersma et al., 2019) and 10xHis-tag (mSA-GP41scFv-His) and cloned it into pET28a to get the plasmid pET28a_ mSA-GP41scFv-His.
All plasmids were transformed into E. coli BL21(DE3) cells. A single colony was cultured in a 5 mL LB medium containing 50 µg/ml kanamycin at 37 °C with 200 rpm shaking for overnight. The cells were then transferred to 500 mL of the LB medium for continued incubation at 37 °C until the optical density at 600 nm (OD600) reached 0.8. IPTG was added to the final concentration of 0.2 mM for induction at 20 °C for 20 hr. Cells were collected and resuspended in Buffer A (50 mM Tris pH 7.4, 800 mM NaCl, and 20 mM imidazole). Cells were then lysed by sonication and centrifuged at 4 °C for 30 min at 12,000 rpm. Solubilized proteins in the supernatant were purified by the Ni-NTA resin (Cube biotech #31103) and eluted with Buffer B (50 mM Tris pH 7.4, 350 mM NaCl, and 500 mM imidazole). The eluted proteins were concentrated to remove imidazole and RNase through a 30 kDa cut-off Amicon Ultra centrifugal filter unit (Millipore #UFC903096). Glycerol was added to the concentrated protein to get the final concentration of 50%. The proteins with 50% glycerol were frozen at –80 °C. The purity of the recombinant proteins was analyzed by SDS-PAGE with Coomassie blue staining. Protein concentration was determined by OD280.
Cell lines and cell culture
The HEK293 and mESC cell lines V6.5 were verified and obtained from Laurie Boyer’s lab at the Massachusetts Institute of Technology (Boyer et al., 2006). Both the HEK293 and V6.5 cell lines are authenticated by the STR profiling and are negative for mycoplasma contamination testing by PCR method. HEK293 cells were maintained in the medium Dulbecco’s Modified Eagle’s Medium (DMEM, Sigma #D5671), 10% FBS (Gibco #10099–141), 1 x non-essential amino acid solution (Sigma #M7145), 1 x GlutaMAX (Sigma #G8541), 1 mM sodium pyruvate (Sigma #S8636), 0.1 mM β-mercaptoethanol (Sigma #M3148), 100 U/mL Penicillin-Streptomycin (Gibco #15140–122) under the condition of 37 °C and 5% CO2. V6.5 cells were maintained on mitomycin C-inactivated mouse embryonic fibroblasts (MEFs) in the mESC medium DMEM, 10% FBS, 2000 U/ml LIF, 1 non-essential amino acid solution, 1x GlutaMAX, 1 mM sodium pyruvate, 0.1 mM β-mercaptoethanol, 100 U/mL Penicillin-Streptomycin under the condition of 37 °C and 5% CO2. The cell lines were routinely checked free of mycoplasma contamination by PCR.
Small molecule treatment
Request a detailed protocolHEK293 cells were treated with DMSO or 25 μM ML216 (Selleck # S0469, inhibitor of BLM helicase) or 3 μM NSC 617145 (MCE # HY-110185, inhibitor of WRN helicase) for 16 hr, and then were collected for different assays.
Generation of Dhx9KO ESCs with CRISPR/Cas9 genome editing
Generation of Cas9 and sgRNAs expressing vectors
Request a detailed protocolWe synthesized the piggyBac 5’ and 3’ inverted repeat and cloned them into pUC57 to obtain pUC57.piggyBac. We then amplified the EF1a promoter (from Addgene# #26777), eSpCas9(1.1) (from Addgene #71814), and IRES-NeoR-WPRE cassette (from Addgene#50917) and ligated them in order to pUC57.piggyBac to get the PiggyBac_EF1a-eSpCas9-IRES-NeoR-WPRE by Gibson Assembly (Gibson et al., 2009). We synthesized the PGK promoter-PuroR-SV40 polyA cassette and two copies of the U6 promoter-guide RNA scaffold (one copy with two Bbs I sites and the other with two Bsa I sites for inserting CRISPR targeting sequence) and ligated them in order to pUC57 to obtain pUC57_sgRNAduo-Puro by Gibson Assembly (Gibson et al., 2009). Two sgRNAs were designed using http://crispor.tefor.net/. Two pairs of oligos for sgRNAs (Key resources table) were annealed and phosphorylated, and cloned into pUC57_sgRNAduo-Puro separately via Bbs I and Bsa I sites to obtain Dhx9_sgRNA1&2.
Generation of dhx9KO ESCs
Request a detailed protocolWe transfected the PiggyBac_EF1a-eSpCas9-IRES-NeoR-WPRE plasmid and Supper PiggyBac Transposase plasmid (Beijing Zoman Biotech) into V6.5 mESCs using the Lonza Nucleofector 2b with the mouse ES cell nucleofector Kit and then treated the transfected cells with G418 (400 µg/ml) for five days and pick up clones individually under microscope to obtain the single clone constitutively expressing eSpCas9 (Cas9 mESCs), which was considered as wildtype cells in the case of comparing with dhx9KO cells. Next, we transfected the plasmid Dhx9_sgRNA1&2 into Cas9 mESCs using the Lonza Nucleofector 2b with the mouse ES cell nucleofector Kit and treated cells with 1 µg/mL puromycin for two days 48 hr after transfection. The clones were picked up individually under a microscope and screened for dhx9KO cells by western blot analysis using an anti-Dhx9 antibody (ABclonal, A4563).
Western blot
Request a detailed protocolmESCs were dissociated using trypsin and removed MEF cells, and then lysed in the RIPA buffer. Proteins in cell lysate were collected by centrifuging at 15,000 g for 15 min at 4 °C, and were separated by SDS-PAGE and electro-transferred to the 0.45 µm PVDF membrane (Millipore #IPVH00010). The membrane was blotted under 5% milk prepared by nonfat-dried milk and PBST (1 x PBS +0.1% Tween-20) and washed by PBST. The primary antibodies were incubated with membranes for overnight at 4 °C, and the HRP-conjugated secondary antibodies were incubated for 30 min at room temperature. Antibodies were listed in the key resources able.
mESC differentiation into embryoid bodies
Request a detailed protocolmESCs were cultured in mESC medium for two passages before differentiation. To initiate differentiation, cells were removed MEFs and diluted to 10,000 cells/mL in EB-Diff medium (standard mESCs medium without LIF). The 2.5 mL diluted cells were plated on the ultra-low-attachment six-well cell culture plate and resuspended cells by shaking the plate to avoid cell adherence. The medium was replaced by EB-Diff medium with 50 µg/mL ascorbic acid (Sigma #A4544) every other day. EBs were allowed to grow for 5 days and analyzed at various times for gene expression.
RNA isolation and quantitative real-time PCR
Request a detailed protocolTotal RNAs were isolated by the acid guanidinium thiocyanate-phenol-chloroform extraction protocol from cells (Chomczynski and Sacchi, 2006). For RT-qPCR analysis, 1 μg total RNA was reverse-transcribed by MonScript RTIII All-in-One Mix (Monad, #MR05101M) following the instruction of the manufacturer. The RT product was subjected to real-time PCR using the SYBR Green PCR master mix and the Roche LC96 Lightcycler. Relative expression values were calculated by the ΔΔCT method using GAPDH for normalization. All reactions were performed in triplicate.
Cell cycle analysis
Request a detailed protocolCells were fixed overnight at 4 °C in 100% ethanol, washed with PBS buffer containing 1% BSA, and stained with 1 µg/mL of DAPI and 1 mg/mL RNase A in PBS for 1 hr at 37 °C. Samples were analyzed on an Attune NxT flow cytometer. The cell cycle was analyzed using ModFit LT.
Cell proliferation assay
Request a detailed protocolBoth WT and dhx9KO ESCs were seeded in 96-well cell culture plates at 5000, 7500, and 10,000 cells/well, respectively (six replicates per condition). Six blank wells were used as controls. Fresh medium was changed after 24 hr. After another 18 hr, cell proliferation was measured using the CCK-8 Cell Proliferation and Cytotoxicity Assay Kit (Solarbio, # CA1210).
Immunofluorescence
Request a detailed protocolCells were fixed with 4% paraformaldehyde for 15 min at room temperature (Figures 6 and 8) or 100% cold methanol for 15 min at 4 °C (Figure 1). The fixed cells were washed three times with PBS and permeabilized in PBS containing 0.2% Triton-X100 and 0.1% Tween-20 for 10 min at room temperature. The permeabilized cells were washed twice with PBS and incubated in PBS containing 10% NBS and 0.1% Tween-20 for 1 hr at room temperature. Cells were stained with primary antibody overnight at 4 °C, and washed three times with PBS containing 0.1% Tween-20. Next, cells were stained with secondary antibody for 30 min at room temperature and washed three times with PBS containing 0.1% Tween-20. Finally, cells were stained with 10 ug/mL Hoechst33342 for 10 min at room temperature, and imaged using the Nikon Eclipse Ti2 microscope or Olympus SpinSR. Antibodies were listed in the key resources table.
RNA sequencing (RNA-seq)
RNA-seq
Request a detailed protocolThe mRNA was first isolated from 1 μg total RNAs using the poly(A) mRNA Capture Module (ABclonal, #RK20340). RNA-seq libraries were prepared from the isolated mRNAs using the Fast RNA-seq Lib Prep Kit V2 (ABclonal, #RK20306). The libraries were finally sequenced on NovaSeq 6000 platform (Illumina) by HaploX Genomics Center.
Bioinformatic analysis
Request a detailed protocolThe pair-end (PE) sequencing reads were analyzed for quality control, filtered, and trimmed off adaptors to get the clean reads using fastp (Chen et al., 2018). Ribosomal RNAs (rRNAs) were removed from clean reads using SortMeRNA (Kopylova et al., 2012). Reads without rRNAs were mapped to the hg38 human genome or mm10 mouse genome using STAR (Dobin et al., 2013) with default settings to generate the bam files. RNA levels of each gene and biotypes of RNAs were quantified using RSEM (Li and Dewey, 2011) and normalized by the FPKM method. The bigwig file were converted from bam files by deepTools (Ramírez et al., 2014) (v2.0) and visualized using IGV (Robinson et al., 2011) genome browser. The differential gene expression analysis was performed using DESeq2 (Love et al., 2014).
HepG4-seq, BG4-seq, HBD-seq, and Dhx9 CUT&Tag
HepG4-seq
Request a detailed protocolHEK293 cells or mESCs were treated with or without 25 μM Hemin (Non-label control) (Sigma, #51280) and 0.5 mM Biotin-aniline (Iris Biotech, #LS-3970) for 2 hr. Next, the H2O2 was added to a final concentration of 1 mM to activate the biotinylation reaction. After 1 min of very gentle shaking, all medium was replaced with quench buffer (1 x PBS, 10 mM Trolox,20 mM sodium ascorbate). After 2 min of very gentle shaking, cells were washed three times with cold PBS and dissociated by TrypLE digestion. Next, cells were washed twice with Buffer 1 (10 mM HEPES pH 7.5, 150 mM NaCl, 0.5 mM spermidine, 1x Protease inhibitor cocktail). Trypan blue was used to detect cell viability above 90%. Live cells (1×105) were taken and conjugated to concanavalin A–coated magnetic beads (Epizyme Biotech, #YJ012) through incubation at room temperature for 15 min. The bead-bound cells were washed three times with Buffer 1 and resuspended in 50 μL AB buffer (Buffer 1 with 1% BSA, 2 mM EDTA, and 0.025% digitonin) containing 4.5 ug recombinant mSA-scFv at 4 °C by rotating overnight. Cells were then washed with DB-1 buffer (Buffer 1 with 0.025% digitonin) twice, and resuspended in 100 µL DB-300 buffer (Buffer 1 with 150 mM NaCl and 0.025% digitonin) containing 2 uL GP41-pG-Tn5 adapter complex and incubated at room temperature for 1 hr with rotation. To prepare the GP41-pG-Tn5 adapter complex, we incubated 12 µl 20 µM recombinant GP41-pG-Tn5 protein and 2.4 ul 100 µM annealed adaptor in 20 µl reaction as described (Kaya-Okur et al., 2019). Cells were washed three times with DB-300 buffer to remove unbound GP41-pG-Tn5 protein. Next, biotinylated genomic DNA was fragmented by activating Tn5 with Mg2+ and the sequencing library of biotinylated DNA fragments was built as described in the CUT&Tag study (Kaya-Okur et al., 2019). The library was finally sequenced on NovaSeq 6000 platform (Illumina) by HaploX Genomics Center.
BG4-seq, HBD-seq
Request a detailed protocolCells were harvested using TrpLE and washed with Buffer 1. Live cells (1×105) were taken and conjugated to concanavalin A–coated magnetic beads (Epizyme Biotech, #YJ012) through incubation at room temperature for 15 min. The bead-bound cells were washed three times with Buffer 1 and resuspended in 50 μL AB buffer containing 2 ug recombinant BG4-V5 (BG4-seq) or 3 μg recombinant HBD-V5 (HBD-seq) at 4 °C by rotating overnight. For the HBD-seq control samples, 1 uL of 20 mg/mL RNase A (NEB, #T3018-2) and 2 µL of 5 U/µL RNase H (NEB, #M0297S) were added to 50 µL of AB buffer and incubated at 37 °C for 4 hr with rotation prior to addition of the recombinant HBD-V5. After overnight incubation at 4 °C, cells were washed once with DB-1 buffer and resuspended in 50 µL DB-1 buffer containing 2 ug anti-V5 (Sino Biological, #100378-T36). After incubation at room temperature for 1 hr with rotation, cells were washed twice with DB-1 buffer and resuspended in 50 µL DB-1 buffer containing 0.25 ug anti-Rabbit IgG (Sino Biological, #SSA018). After incubation at room temperature for 30 min with rotation, cells were washed twice with DB-1 buffer and tagmented with GP41-pG-Tn5 adapter complex as HepG4-seq. The sequencing library was prepared as CUT&Tag (Kaya-Okur et al., 2019) using the tagmented fragments. The library was finally sequenced on NovaSeq 6000 platform (Illumina) by HaploX Genomics Center.
Dhx9 CUT&Tag
Request a detailed protocolCells were harvested using TrpLE and washed with Buffer 1. Live cells (1×105) were taken and conjugated to concanavalin A–coated magnetic beads (Epizyme Biotech, #YJ012) through incubation at room temperature for 15 min. Cells were fixed with 0.1% formaldehyde (Cell Signaling, #50-00-0) for 2 min as described (Kaya-Okur et al., 2019). The bead-bound cells were washed three times with Buffer 1 and resuspended in 50 μL AB buffer containing 0.04 ug Dhx9 antibody (ABclonal, #A4563) at 4 °C by rotating overnight. After overnight incubation at 4 °C, cells were washed once with DB-1 buffer and resuspended in 50 µL DB-1 buffer 0.25 ug anti-Rabbit IgG (Sino Biological, #SSA018). After incubation at room temperature for 30 mins with rotation, cells were washed twice with DB-1 buffer and tagmented with pG-Tn5 adapter complex (Xu et al., 2021). The sequencing library was prepared as CUT&Tag (Kaya-Okur et al., 2019) using the tagmented fragments. The library was finally sequenced on NovaSeq 6000 platform (Illumina) by HaploX Genomics Center.
Bioinformatic analysis
Request a detailed protocolThe PE sequencing reads were analyzed for quality control, filtered, and trimmed off adaptors to get the clean reads using fastp (Chen et al., 2018) and Cutadapt (Martin, 2011). Clean Reads were mapped to hg38 human genome or mm10 mouse genome using Bowie2 (Langmead and Salzberg, 2012) with the settings (--end-to-end --very-sensitive --no-mixed --no-discordant --phred33 -I 10 X 500) to generate the bam files. High-quality mapping reads were extracted by setting parameter samtools -q to 40 (Cheng et al., 2009). Duplicate reads were removed using the Picard toolkit (RRID:SCR_006525,https://broadinstitute.github.io/picard/). Only uniquely mapping and de-duplicated reads (quality score >40) were retained. To normalize reads coverage, scale factors were calculated by ChIPseqSpikeInFree software (Jin et al., 2020) with default parameters. Given the difference between DMSO and inhibitor treatment is too big, ChIPseqSpikeInFree failed to calculate the scale factor and thereby we calculated the scale factor using the total numbers of uniquely mapping reads. To create coverage plots or heatmaps, the reads coverage were normalized by setting the scaleFactor parameters and converted to the bigwig files from bam files using deepTools (Ramírez et al., 2014). The heatmaps and profile plots were generated from normalized bed files using deepTools (Ramírez et al., 2014). The bigwig files were visualized using the IGV genome browser (Robinson et al., 2011). The bam files were converted to normalized bedgraph files using the BEDTools genomecov command and the scale factor calculated above. Peaks were called from normalized bedgraph files using SEACR (Meers et al., 2019) with the settings (0.05 non-stringent). High confident peaks were consensus peaks between two replicates generated with BEDTools (Quinlan and Hall, 2010) intersect command. Co-localized G4s and R-loops peaks are overlapping G4 and DNA/RNA hybrid peaks generated with BEDTools (Quinlan and Hall, 2010) intersect command. To detect the differential peaks, peaks from different groups were first re-centered and trimmed based on calculating their reads coverage within a fix bin (400 bp) using R package DiffBind (v3.12); the read coverage of each peak was normalized using the scale factor calculated above; the significantly differential peaks were defined as peaks with the significance cut-off p-value <0.05 and fold change (FC) ≥1.5; BEDTools were then used to identify overlaps between the significantly differential peaks and original peaks from SEACR to retrieve original peak sites with differential levels in different groups that were used subsequently as input file of HOMER software (Heinz et al., 2010). Genomic localization of peaks was annotated by HOMER annotatePeaks.pl (Heinz et al., 2010). Motif enrichment of peaks was analyzed by HOMER findMotifsGenome.pl (Heinz et al., 2010). G4/DNA/RNA hybrid/Co-localized G4s and R-loops-associated genes are the nearest genes annotated by HOMER to G4/DNA/RNA hybrid/Co-localized G4s and R-loops peaks. MapR, GST-2xHBD CUT&Tag, S9.6 CUT&Tag were analyzed using the same pipeline above.
Analysis of histone marks, promoters, and enhancers
Request a detailed protocolThe bioinformatic analysis pipeline of these ChIP-seq data is the same as the HepG4-seq except that ChIP-seq peaks were called from bam files using MACS2 with default settings and the input DNAs were used as the control for MACS2. For histone marks, promoters are defined as a 2 kb region centered on the TSS of each Ensemble transcript as Wamstad et al., 2012. The enhancers was generated by combining all H3K4me1 or H3K27ac enriched regions and further filtered by removing any interval that overlapped with regions +/-1 kb of TSS annotated by mm10 genome build as Wamstad et al., 2012. If an extended region overlapped the +/-1 kb window around a TSS, the region was trimmed +/-5 kb of the TSS, leaving the distal portion in either direction as independent enhancer regions (Wamstad et al., 2012).
GO enrichment analysis and STRING protein-protein interaction analysis
Request a detailed protocolGO enrichment analysis was performed using the Metascape web interface (Zhou et al., 2019b) and G4/DNA/RNA hybrid/Co-localized G4s and R-loops-associated genes’ ID as input. The STRING protein-protein physical interaction analysis was performed using the STRING web interface (Szklarczyk et al., 2023) with the settings (physical subnetwork, network edges indicate confidence, medium confidence 0.4).
Data availability
The HepG4-seq, HBD-seq, BG4-seq, Dhx9 CUT&Tag and RNA-seq data have been deposited to the Gene Expression Omnibus (accession code GSE254764 and GSE254763). The ChIP-seq data of histone markers and RNAP are openly available in GEO (GSE47950; Wamstad et al., 2012). Data of BG4 CUT&Tag in mESCs, S9.6 CUT&Tag in mESCs, GST-2xHBD CUT&Tag in HEK293 cells and MapR in HEK293 cells are openly available in GEO (GSE173103 and GSE156400).
-
NCBI Gene Expression OmnibusID GSE254764. Genome-wide mapping of native G-loops in living cells.
-
NCBI Gene Expression OmnibusID GSE254763. RNA-seq analysis of wild-type HEK293, HEK293 treated with DMSO, ML216 or NSC617145, wild-type mESCs and dhx9KO mESCs.
-
NCBI Gene Expression OmnibusID GSE47950. Dynamic and Coordinated Epigenetic Regulation of Developmental Transitions in the Cardiac Lineage.
-
NCBI Gene Expression OmnibusID GSE173103. Genome-wide mapping of G-quadruplex structures with CUT&Tag.
-
NCBI Gene Expression OmnibusID GSE156400. Genomic Profiling of Native R-Loops with a DNA-RNA Hybrid Recognition Sensor Protein.
References
-
The free heme concentration in healthy human erythrocytesBlood Cells, Molecules & Diseases 55:402–409.https://doi.org/10.1016/j.bcmd.2015.09.003
-
Dual roles of RNA helicase A in CREB-dependent transcriptionMolecular and Cellular Biology 21:4460–4469.https://doi.org/10.1128/MCB.21.14.4460-4469.2001
-
The interplay of epigenetic marks during stem cell differentiation and developmentNature Reviews. Genetics 18:643–658.https://doi.org/10.1038/nrg.2017.57
-
Multimodal chromatin profiling using nanobody-based single-cell CUT&TagNature Biotechnology 41:794–805.https://doi.org/10.1038/s41587-022-01535-4
-
Unraveling epigenetic regulation in embryonic stem cellsCell Stem Cell 2:123–134.https://doi.org/10.1016/j.stem.2008.01.005
-
Modification of enhancer chromatin: what, how, and why?Molecular Cell 49:825–837.https://doi.org/10.1016/j.molcel.2013.01.038
-
The DNA/RNA helicase DHX9 contributes to the transcriptional program of the androgen receptor in prostate cancerJournal of Experimental & Clinical Cancer Research 41:178.https://doi.org/10.1186/s13046-022-02384-4
-
fastp: an ultra-fast all-in-one FASTQ preprocessorBioinformatics 34:i884–i890.https://doi.org/10.1093/bioinformatics/bty560
-
BookExpression and purification of large active GST fusion enzymesIn: Labrou NE, editors. Protein Downstream Processing: Design, Development and Application of High and Low-Resolution Methods. Humana Press. pp. 169–180.https://doi.org/10.1007/978-1-62703-977-2_15
-
STAR: ultrafast universal RNA-seq alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Self-biotinylation of DNA G-quadruplexes via intrinsic peroxidase activityNucleic Acids Research 45:9813–9822.https://doi.org/10.1093/nar/gkx765
-
Human werner syndrome DNA helicase unwinds tetrahelical structures of the fragile X syndrome repeat sequence d(CGG)nThe Journal of Biological Chemistry 274:12797–12802.https://doi.org/10.1074/jbc.274.18.12797
-
DNA G-Quadruplex recognition in vitro and in live cells by a structure-specific nanobodyJournal of the American Chemical Society 144:23096–23103.https://doi.org/10.1021/jacs.2c10656
-
G-quadruplexes sequester free heme in living cellsCell Chemical Biology 26:1681–1691.https://doi.org/10.1016/j.chembiol.2019.10.003
-
G-quadruplex structures mark human regulatory chromatinNature Genetics 48:1267–1272.https://doi.org/10.1038/ng.3662
-
DNA G-quadruplexes in the human genome: detection, functions and therapeutic potentialNature Reviews Molecular Cell Biology 18:279–284.https://doi.org/10.1038/nrm.2017.3
-
The affinity of the s9.6 antibody for double-stranded rnas impacts the accurate mapping of r-loops in fission yeastJournal of Molecular Biology 430:272–284.https://doi.org/10.1016/j.jmb.2017.12.016
-
The RNA polymerase II CTD coordinates transcription and RNA processingGenes & Development 26:2119–2137.https://doi.org/10.1101/gad.200303.112
-
DHX9 helicase is involved in preventing genomic instability induced by alternatively structured DNA in human cellsNucleic Acids Research 41:10345–10357.https://doi.org/10.1093/nar/gkt804
-
Fusion tags to enhance heterologous protein expressionApplied Microbiology and Biotechnology 104:2411–2425.https://doi.org/10.1007/s00253-020-10402-8
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Chromatin-associated RNAs as facilitators of functional genomic interactionsNature Reviews. Genetics 20:503–519.https://doi.org/10.1038/s41576-019-0135-1
-
Stable, high-affinity streptavidin monomer for protein labeling and monovalent biotin detectionBiotechnology and Bioengineering 110:57–67.https://doi.org/10.1002/bit.24605
-
Peak calling by sparse enrichment analysis for CUT&RUN chromatin profilingEpigenetics & Chromatin 12:42.https://doi.org/10.1186/s13072-019-0287-4
-
G-quadruplexes and helicasesNucleic Acids Research 44:1989–2006.https://doi.org/10.1093/nar/gkw079
-
G-quadruplex-R-loop interactions and the mechanism of anticancer G-quadruplex bindersNucleic Acids Research 48:11942–11957.https://doi.org/10.1093/nar/gkaa944
-
The Bloom’s and Werner’s syndrome proteins are DNA structure-specific helicasesNucleic Acids Research 29:2843–2849.https://doi.org/10.1093/nar/29.13.2843
-
CTCF: an architectural protein bridging genome topology and functionNature Reviews. Genetics 15:234–246.https://doi.org/10.1038/nrg3663
-
The sub-nanomolar binding of DNA-RNA hybrids by the single-chain Fv fragment of antibody S9.6Journal of Molecular Recognition 26:376–381.https://doi.org/10.1002/jmr.2284
-
deepTools: a flexible platform for exploring deep-sequencing dataNucleic Acids Research 42:W187–W191.https://doi.org/10.1093/nar/gku365
-
Genome engineering using the CRISPR-Cas9 systemNature Protocols 8:2281–2308.https://doi.org/10.1038/nprot.2013.143
-
Divergent transcription from active promotersScience 322:1849–1851.https://doi.org/10.1126/science.1162253
-
Insights into G-Quadruplex-Hemin dynamics using atomistic simulations: Implications for reactivity and foldingJournal of Chemical Theory and Computation 17:1883–1899.https://doi.org/10.1021/acs.jctc.0c01176
-
Polymerase η Recruits DHX9 helicase to promote replication across guanine quadruplex structuresJournal of the American Chemical Society 144:14016–14020.https://doi.org/10.1021/jacs.2c05312
-
The regulation and functions of DNA and RNA G-quadruplexesNature Reviews. Molecular Cell Biology 21:459–474.https://doi.org/10.1038/s41580-020-0236-x
-
G-quadruplexes associated with R-loops promote CTCF bindingMolecular Cell 83:3064–3079.https://doi.org/10.1016/j.molcel.2023.07.009
-
RNA-DNA hybrid formation at the human mitochondrial heavy-strand origin ceases at replication start sites: an implication for RNA-DNA hybrids serving as primersThe EMBO Journal 15:3135–3143.
-
Method for efficient soluble expression and purification of recombinant hyperactive Tn5 transposaseProtein Expression and Purification 183:105866.https://doi.org/10.1016/j.pep.2021.105866
-
Heterozygous loss-of-function DHX9 variants are associated with neurodevelopmental disorders: Human genetic and experimental evidencesEuropean Journal of Medical Genetics 66:104804.https://doi.org/10.1016/j.ejmg.2023.104804
-
Characterization of G‐Quadruplex/Hemin peroxidase: substrate specificity and inactivation kineticsChemistry – A European Journal 17:14475–14484.https://doi.org/10.1002/chem.201101941
-
Helicases in R-loop formation and resolutionThe Journal of Biological Chemistry 299:105307.https://doi.org/10.1016/j.jbc.2023.105307
-
TDRD3 promotes DHX9 chromatin recruitment and R-loop resolutionNucleic Acids Research 49:8573–8591.https://doi.org/10.1093/nar/gkab642
-
Expanding APEX2 substrates for proximity‐dependent labeling of nucleic acids and proteins in living cellsAngewandte Chemie International Edition 58:11763–11767.https://doi.org/10.1002/anie.201905949
Article and author information
Author details
Funding
National Key Research and Development Program of China Stem Cell and Translational Research (2018YFA0109200)
- Zhihong Xue
Project for Cultivating Leading Talent in Science and Technology of Sichuan University (0082604151357)
- Zhihong Xue
Postdoctoral Research Foundation of China (2018M643467)
- Lin Yu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Laurie Boyer at the Massachusetts Institute of Technology for providing the HEK293 and V6.5 cell lines. We thank Yi Liu (University of Texas Southwestern Medical Center, USA) and Fan Lai (Yunnan University, China) for insightful discussions and for critical evaluation of the manuscript. Z X is supported by the National Key Research and Development Program of China, Stem Cell and Translational Research (2018YFA0109200), and Project for Cultivating Leading Talent in Science and Technology of Sichuan University (0082604151357). LY is supported by the Postdoctoral Research Foundation of China (2018M643467).
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.99026. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Liu, Shen, Ren et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 893
- views
-
- 56
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Cell Biology
Myosin 10 (Myo10) is a motor protein known for its role in filopodia formation. Although Myo10-driven filopodial dynamics have been characterized, there is no information about the absolute number of Myo10 molecules during the filopodial lifecycle. To better understand molecular stoichiometries and packing restraints in filopodia, we measured Myo10 abundance in these structures. We combined SDS-PAGE densitometry with epifluorescence microscopy to quantitate HaloTag-labeled Myo10 in U2OS cells. About 6% of total intracellular Myo10 localizes to filopodia, where it enriches at opposite cellular ends. Hundreds of Myo10s are in a typical filopodium, and their distribution across filopodia is log-normal. Some filopodial tips even contain more Myo10 than accessible binding sites on the actin filament bundle. Live-cell movies reveal a dense cluster of over a hundred Myo10 molecules that initiates filopodial elongation. Hundreds of Myo10 molecules continue to accumulate during filopodial growth, but accumulation ceases when retraction begins. Rates of filopodial elongation, second-phase elongation, and retraction are inversely related to Myo10 quantities. Our estimates of Myo10 molecules in filopodia provide insight into the physics of packing Myo10, its cargo, and other filopodia-associated proteins in narrow membrane compartments. Our protocol provides a framework for future work analyzing Myo10 abundance and distribution upon perturbation.
-
- Biochemistry and Chemical Biology
Ion channels are biological transistors that control ionic flux across cell membranes to regulate electrical transmission and signal transduction. They are found in all biological membranes and their conductive state kinetics are frequently disrupted in human diseases. Organelle ion channels are among the most resistant to functional and pharmacological interrogation. Traditional channel protein reconstitution methods rely upon exogenous expression and/or purification from endogenous cellular sources which are frequently contaminated by resident ionophores. Here, we describe a fully synthetic method to assay functional properties of polycystin channels that natively traffic to primary cilia and endoplasmic reticulum organelles. Using this method, we characterize their oligomeric assembly, membrane integration, orientation, and conductance while comparing these results to their endogenous channel properties. Outcomes define a novel synthetic approach that can be applied broadly to investigate channels resistant to biophysical analysis and pharmacological characterization.